Hydroclimatic Variability and Floodplain Wetland Dynamics in the Magdalena River: A Case Study of Zambrano, Colombia

Understanding the interactions between vegetation dynamics and surface water availability is essential for assessing the resilience of tropical floodplain ecosystems under increasing hydroclimatic variability. This study analyzes the spatio-temporal dynamics of vegetation cover, surface water, and land use in Zambrano, a floodplain-dominated sector of the lower Magdalena River basin (Colombian Caribbean), over the period 1990–2025. Multi-temporal Landsat imagery was used to derive the Normalized Difference Vegetation Index (NDVI) and the Normalized Difference Water Index (NDWI), enabling the evaluation of seasonal and interannual ecohydrological variability under contrasting dry and rainy conditions. In addition, land-use classification was performed using a CORINE Land Cover methodology adapted for Colombia (CLC-C) to characterize the spatial organization of the landscape and its influence on vegetation–water interactions. Results show that vegetation dynamics are strongly controlled by hydroclimatic seasonality. Dense vegetation consistently expands during rainy periods, while dry seasons promote the expansion of open and sparse vegetation, reflecting seasonal vegetation stress rather than long-term degradation. NDWI patterns indicate that surface water and soil moisture are highly seasonal and spatially constrained, with open water largely confined to the Magdalena River channel and localized floodplain depressions. Extreme hydroclimatic events associated with the El Niño–Southern Oscillation (ENSO) produce abrupt but temporary changes in vegetation structure and surface moisture distribution. A strong inverse correlation between NDVI and NDWI reflects the contrasting spectral responses of vegetation and water surfaces resulting from the shared near-infrared (NIR) band in both indices. This spectral relationship is consistent with the observed seasonal variations in vegetation greenness and surface moisture across the floodplain. Land-use analysis reveals the progressive consolidation of the landscape, where the agropastoral matrix expanded from ~18,000 ha in 1990 to over 22,000 ha by 2025, driving a systematic reduction in natural and semi-natural forest structures. Forest conservation areas serve as critical ecological buffers, exhibiting lower seasonal variability in vegetation greenness. Overall, the results indicate that the Zambrano floodplain functions as a structurally stable yet highly responsive ecohydrological system, where vegetation dynamics and surface water availability are predominantly governed by interannual hydroclimatic pulses rather than long-term directional degradation. These findings demonstrate that while the structural matrix of the floodplain exhibits strong baseline resilience, its ecological functioning remains critically coupled with, and vulnerable to, the extreme phase shifts in ENSO cycles.