We use Recurrence Quantification Analysis (RQA) to study features of Generalized Synchronization (GS) between El Niño-Southern Oscillation (ENSO) and monthly hydrological anomalies (HyAns) of rainfall and streamflows in Colombia. To that end, we check the sensitivity of the RQA concerning diverse HyAns estimation methods, which constitutes a fundamental procedure for any climatological analysis at inter-annual timescales. In general, the GS and its sensitivity to HyAns methods are quantified by means of time-lagged joint recurrence analysis. Then, we link the GS results with the dynamics of major physical mechanisms that modulate Colombia's hydroclimatology, including the Caribbean, the CHOCO and the Orinoco Low-Level Jets (LLJs), and the Cross-Equatorial Flow (CEF) over northwestern Amazonia (southern Colombia). Our findings show that RQA exhibits significant differences depending on the HyAns methods. GS results are similar for the HyAns methods with variable annual cycle but the time-lags seem to be sensitive. On the other hand, our results make evident that HyAns in the Pacific, Caribbean, and Andean regions of Colombia exhibit strong (weak) GS with the ENSO signal during La Niña (El Niño), when hydrological anomalies are positive (negative). Results from the GS analysis allow us to identify spatial patterns of non-linear dependence between ENSO and the Colombian's climatology. The mentioned moisture transport sources constitute the interdependence mechanism and contribute to explain hydrological anomalies in Colombia during the phases of ENSO. During La Niña (El Niño), GS is strong (weak) for the Caribbean and the CHOCO LLJs whereas GS is moderate (strong) for the Orinoco LLJ. Moreover, moisture advection by the Caribbean and CHOCO LLJs exhibit synchrony with HyAns at 0–2 (2–4) months-lags over north-western Colombia and the Orinoco LLJ moisture advection synchronizes with HyAns at similar month-lags over the Amazon region of Colombia. Furthermore, our results suggest a strong (weak) GS between negative (positive) Sea Surface Temperatures (SST) anomalies in the Eastern Pacific and rainfall anomalies in Colombia. In contrast, GS is strong (weak) for positive (negative) SST anomalies in the Central Pacific. Our GS results contribute to advance our understanding on the regional effects of both phases of ENSO in Colombia, whose socio-economical, environmental and ecological impacts cannot be overstated. This work provides a novel approach that reveals new insights into the impact of ENSO on northern South America.

The synoptic mode of variability (SMV) refers to changes in atmospheric conditions over periods ranging from 2 to 10 days. In tropical regions, this variability is driven by tropical waves that have a clear signal on the wavenumber–frequency power spectra of precipitation. This study uses the ensemble empirical mode decomposition (EEMD) method to identify the SMV in daily precipitation and streamflows in 47 river basins over northern South America. We found the presence of the frequency bands with periods of 3–12 days and 6–18 days, which agrees with the SMV associated with tropical waves that modulate precipitation over the region. Furthermore, our results reveal that variance explained by the SMV in rainfall over each catchment is greater than the variance explained by those SMV in streamflows, which suggests that catchments efficiently filter out this variability. We found that SMV explains from 5% to 20% of streamflow variability for catchments ranging from 1000 km${}^2$ to 5000 km${}^2$. Additionally, the variance explained by SMV decreases as a power fit with the catchment area. Thus, this study characterizes the SMV for potential applications on regional hydrology, diagnosis, modeling, short-time forecasting, prediction, and management of water resources.

Environmental pollution indicated by the presence of $P{M}_{2.5}$ particulate matter varies based on prevailing atmospheric conditions described by certain meteorological variables. Consequently, it is important to understand atmospheric behavior in areas such as the Aburrá Valley, which experiences recurrent pollution events twice a year. This study examines the behavior of specific meteorological variables and $P{M}_{2.5}$ particulate matter in the Aburrá Valley. By using statistical analysis tools such as correlation coefficients, principal component analysis (PCA), and multiple linear regression models, the research identifies relationships between $P{M}_{2.5}$ and daily cycles of temperature, rainfall, radiation, and wind speed and direction. Datasets were analyzed considering periods before and after the COVID-19 lockdown (pre-pandemic and pandemic, respectively), and specific pollution events were also analyzed. Furthermore, this work considers the relationships between $P{M}_{2.5}$ and meteorological variables, contrasting the pre-pandemic and pandemic periods. This study characterizes diurnal cycles of meteorological variables and their relationship with $P{M}_{2.5}$. There are consistent patterns among temperature, atmospheric boundary layer (ABL) height, and solar radiation, whereas precipitation and relative humidity show the opposite behavior. $P{M}_{2.5}$ exhibits similar relative frequency functions during both daytime and nighttime, regardless of rainfall. An inverse relationship is noted between $P{M}_{2.5}$ levels and ABL height at different times of the day. Moreover, the PCA results show that the first principal component explains around 60% of the total variance in the hydrometeorological data. The second PC explains 10%, and the rest of the variance is distributed among the other three to eight PCs. In this sense, there is no significant difference between the two PCAs with hydrometeorological data from a pre-pandemic period and a COVID-19 pandemic period. Multiple regression analysis indicates a significant and consistent dependence of $P{M}_{2.5}$ on temperature and solar radiation across both analyzed periods. The application of Generalized Additive Models (GAMs) to our dataset yielded promising results, reflecting the complex relationship between meteorological variables and $P{M}_{2.5}$ concentrations. The metrics obtained from the GAM were as follows: Mean Squared Error (MSE) of 98.04, Root Mean Squared Error (RMSE) of 9.90, R-squared (R²) of 0.24, Akaike Information Criterion (AIC) of 110,051.34, and Bayesian Information Criterion (BIC) of 110,140.63. In comparison, the linear regression model exhibited slightly higher MSE (100.49), RMSE (10.02), and lower R-squared (0.22), with AIC and BIC values of 110,407.45 and 110,460.67, respectively. Although the improvement in performance metrics from GAM over the linear model is not conclusive, they indicate a better fit for the complexity of atmospheric dynamics influencing $P{M}_{2.5}$ levels. These findings underscore the intricate interplay of meteorological factors and particulate matter concentration, reinforcing the necessity for advanced modeling techniques in environmental studies. This work presents new insights that enhance the diagnosis, understanding, and modeling of environmental pollution, thereby supporting informed decision-making and strengthening management efforts.

We use a Generalized Watershed Loading Function (GWLF) model to simulate streamflow in the Gualí River Basin. The model's performance is assessed using metrics such as percentage of bias (PBIAS), Nash–Sutcliffe efficiency (NSE), RMSE–observations standard deviation ratio (RSR), and Kling–Gupta efficiency (KGE) to indicate good performance. Furthermore, we analyze projections of precipitation and streamflow using several global climate models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6) and three shared socioeconomic pathways (SSP1-2.6, SSP2-4.5, SSP5-8.5). Despite the uncertainties and coarse resolution, our results show that increases in the mean streamflow and significantly decreasing trends in projected precipitation and streamflow are observed from 2015 to 2099 under the SSP5-8.5 scenario. Furthermore, our findings suggest an increase in long-term mid-flow and low-flow. Moreover, this work provides a methodological framework for hydrological modeling in small tropical river basins, by incorporating data from GCMs while raising concerns and caveats. This study offers valuable insights into the potential effects of climate change on streamflow in an Andean river basin characterized by volcanic activity and significant human impacts. The findings reported here provide useful information for future decisions related to water supply for the social, environmental, and productive sectors in the seven towns within the catchment.

We study spatial scaling and complexity properties of Amazonian radar rainfall fields using the Beta-Lognormal Model (BL-Model) with the aim to characterize and model the process at a broad range of spatial scales. The Generalized Space q-Entropy Function (GSEF), an entropic measure defined as a continuous set of power laws covering a broad range of spatial scales, $S_q\left(\lambda \right)\sim {\lambda }^{\Omega \left(q\right)}$, is used as a tool to check the ability of the BL-Model to represent observed 2-D radar rainfall fields. In addition, we evaluate the effect of the amount of zeros, the variability of rainfall intensity, the number of bins used to estimate the probability mass function, and the record length on the GSFE estimation. Our results show that: (i) the BL-Model adequately represents the scaling properties of the q-entropy, $S_q$, for Amazonian rainfall fields across a range of spatial scales $\lambda$ from 2 km to 64 km; (ii) the q-entropy in rainfall fields can be characterized by a non-additivity value, ${\mathit{q}}_{sat}$, at which rainfall reaches a maximum scaling exponent, ${\Omega }_{sat}$; (iii) the maximum scaling exponent ${\Omega }_{sat}$ is directly related to the amount of zeros in rainfall fields and is not sensitive to either the number of bins to estimate the probability mass function or the variability of rainfall intensity; and (iv) for small-samples, the GSEF of rainfall fields may incur in considerable bias. Finally, for synthetic 2-D rainfall fields from the BL-Model, we look for a connection between intermittency using a metric based on generalized Hurst exponents, $M(q_1,q_2)$, and the non-extensive order (q-order) of a system, ${\Theta }_q$, which relates to the GSEF. Our results do not exhibit evidence of such relationship.