Groundwater-surface water interactions are important in controlling lake water residence time, biogeochemistry, and water availability for downstream communities in tropical volcanic catchments. To better understand these complex seasonal interactions, a multi-tracer approach including water and inorganic carbon stable isotopes (δ2H, δ18O, δ13CDIC), hydrochemistry, and 222Rn was applied in Lake Hule, northern Costa Rica. Seasonal isotope mass balance calculations using lake, stream, precipitation, and groundwater isotope compositions were supplemented with local hydrometeorological information. Evaporation to inflow ratios (E/I) revealed a small variability between the dry (December-April) and wet seasons (May-November), with relatively low evaporation losses, 2.9±1.0 % and 3.2±1.8 %, respectively. Bayesian end-member analysis indicated that annual inputs from groundwater, precipitation, and runoff represented 61.3±8.1%, 24.4±8.4, and 14.3±5.9% of total inflow, respectively. Temporal variations of δ13CDIC also confirmed the key role carbonate buffering plays in this lake and indicated greater CO2 degassing from groundwater sources in the wet season. This first tracer-aided assessment in a volcanic lake maar of northern Costa Rica provides evidence of previously unknown groundwater-surface water interactions and poses a promising tool for estimating seasonal variability of groundwater discharge into natural lakes across the volcanic front of Central America.

A document by Fabian Quichimbo. Click on the document to view its contents.

Climate change, inter-annual precipitation variability, recurrent droughts, and flash flooding, coupled with increasing water needs, are shaping the co-evolution of socioeconomic and cultural assemblages, water laws and regulations, and equitable drinking water access and allocation worldwide. Recognizing the need for mitigation strategies for drinking water availability in urban areas, the Isotope Hydrology Section of the International Atomic Energy Agency (IAEA) coordinated a state-of-the-art global assessment to evaluate water sources and distribution of drinking water supply in urban centers, an initiative entitled “Use of Isotope Techniques for the Evaluation of Water Sources for Domestic Supply in Urban Areas (2018-2023)”. Here, we report on a) current research trends for studying urban drinking water systems during the last two decades and b) the development, testing, and integration of new methodologies, aiming for a better assessment, mapping, and management of water resources used for drinking water supply in urban settings. Selected examples of water isotope applications (Canada, USA, Costa Rica, Ecuador, Morocco, Botswana, Romania, Slovenia, India, and Nepal) provide context to the insights and recommendations reported and highlight the versatility of water isotopes to underpin seasonal and temporal variations across various environmental and climate scenarios. The study revealed that urban areas depend on a large spectrum of water recharge across mountain ranges, extensive local groundwater extraction, and water transfer from nearby or distant river basins. The latter is reflected in the spatial isotope snapshot variability. High-resolution monitoring (hourly and sub-hourly) isotope sampling revealed large diurnal variations in the wet tropics (Costa Rica) (up to 1.5‰ in δ 18O) and more uniform diurnal variations in urban centers fed by groundwater sources (0.08 ‰ in δ 18O) ([Ljubljana](https://www.google.com/search?client=firefox-b-1-d&sca_esv=f5a20a2e9138d638&sca_upv=1&sxsrf=ADLYWIKR6-DvBtjaWqFYRhn6VgnegOa8kg:1717189104058&q=Ljubljana&stick=H4sIAAAAAAAAAONgVuLQz9U3SMrNNXnEaMwt8PLHPWEprUlrTl5jVOHiCs7IL3fNK8ksqRQS42KDsnikuLjgmngWsXL6ZJUm5WQl5iUCAAFa64FOAAAA&sa=X&ved=2ahUKEwjMrrz047iGAxWyG9AFHSVwCBgQzIcDKAB6BAgTEAE), Slovenia). Similarly, while d-excess was fairly close to the global mean value (+10 ‰) across all urban centers (10-15‰), reservoir-based drinking water systems show significantly lower values (up to ~ -20 ‰) (Arlington, TX, USA and Gaborone, Botswana), as a result of strong evapoconcentration processes. δ 18O time series and depth-integrated sampling highlighted the influence of the catchment damping ratio in the ultimate intake water composition. By introducing new, traceable spatial and temporal tools that span from the water source to the end-user and are linked to the engineered and socio-economic structure of the water distribution system, governmental, regional, or community-based water operators and practitioners could enhance drinking water treatment strategies (including more accurate surface water blending estimations) and improve urban water management and conservation plans in the light of global warming.

Extracting water from discrete xylem and soil samples and continuous (long-term) monitoring of water vapor across the soil-plant-atmosphere continuum remains challenging and under a large debate. Here, we present a detailed one-year study in the Fort Worth Botanic Garden (north-central Texas) to evaluate the analytical robustness of the centrifugation extraction method and understand water sourcing from three common urban tree species (Elderberry, Sambucus canadensis; Cherry Laurel, Prunus caroliniana; and Boxelder Maple, Acer negundo). Xylem (N=110) isotope ratios (δ 18O and δ 2H) are compared to local precipitation (N=498), throughfall (N=33), and soil water (N=105) at different depths (0-38 cm). Complementary soil water samples were obtained from cup suction lysimeters (N=42) (0-38 cm). Soil and xylem water extraction volumes ranged from 100 µL to 7.5 mL in plant samples and from 100 µL to 10.5 mL in soil samples. Extraction success rates were 68.8% and 75.2% for xylem and soil samples, respectively. The minimum sample total water content for effective extractions was determined as 10.6% (soil) and 17.8% (xylem). Xylem mean narrowband and broadband (proxy for organic contamination) were 0.23±0.40 (-) and 1.00±0.01 (-), respectively. These values agree with mean narrowband and broadband metrics from throughfall and soils, which highlight the non-invasive nature of centrifugated extractions. Annual mean soil δ 18O compositions (-3.6±1.7‰) corresponded with the throughfall input (-3.6±2.4‰). Xylem δ 18O compositions exhibited a strong temporal enrichment trend at the end of the winter, summer, and fall seasons. Mean spring xylem δ 18O (-2.85‰) was less variable and close to soil mean compositions (-2.82‰). For this season, Bayesian mixing analysis showed source water contributions from distinct soil depths: 0 cm to 12.7 cm for Boxelder Maple, 12.7 cm to 25.4 cm for Cherry Laurel, and 12.7 cm to 38.1 cm for Elderberry. Our results offer a standardized and effective protocol for centrifugation extractions and reveal plant water uptake preferences in a highly altered urban green space during an unprecedented warm year.

The spatial variation of soil water isotopes (SWI) - representing the baseline for investigating root water uptake (RWU) depths with water stable isotope techniques - has rarely been investigated. Here, we use spatial SWI depth profile sampling in combination with unmanned aerial vehicle (UAV) based land surface temperature estimates and vegetation indices (VI) in order to improving process understanding of the relationships between soil water content and isotope patterns with canopy status. We carried out a spatial sampling of ten SWI depth profiles in a tropical dry forest. UAV data were collected and analyzed to obtain detailed characterization of soil temperature and canopy status. We then performed a statistical analysis between the VI and land surface temperatures with soil water content and SWI values at different spatial resolutions (3 cm to 5 m). Best relationships were used for generating soil water isoscapes for the entire study area. Results suggest that soil water content and SWI values are strongly mediated by canopy parameters (VI). Various VI correlate strongly with soil water content and SWI values across all depths. SWI at the surface depend on land surface temperature (R² of 0.65 for δ18O and 0.57 for δ2H). Strongest overall correlations were found at a spatial resolution of 0.5 m. We speculate that this might be the ideal resolution for spatially characterizing SWI patterns and investigate RWU. Supporting spatial analyses of SWI with UAV-based approaches might be a future avenue for improving the spatial representation and credibility of such studies.

The El Niño-Southern Oscillation (ENSO) phenomena, originating in the tropical Pacific region, is an interannual climate variability driven by sea surface temperature and atmospheric pressure changes that affect weather patterns globally. In Mesoamerica, ENSO can cause significant changes in rainfall patterns with major impacts on water resources. This commentary presents results from a nearly 10-yr hydrometric and tracer monitoring network across north-central Costa Rica, a region known as a headwater-dependent system. This monitoring system has recorded different El Niño and La Niña events, as well as the direct/indirect effects of several hurricane and tropical storm passages. Our results show that ENSO exerts a significant but predictable impact on rainfall anomalies, groundwater recharge, and spring discharge, as evidenced by second-order water isotope parameters (e.g., line conditioned-excess or LC-excess). The Oceanic Niño Index (ONI) is correlated with a reduction in mean annual and cold front rainfall across the headwaters of north-central Costa Rica. During El Niño conditions, rainfall is substantially reduced (by up to 69.2%) during the critical cold fronts period, subsequently limiting groundwater recharge and promoting an early onset of baseflow conditions. In contrast, La Niña is associated with increased rainfall and groundwater recharge (by up to 94.7% during active cold front periods). During La Niña, the long-term mean spring discharge (39 Ls -1) is exceeded 63-80% of the time, whereas, during El Niño, the exceedance time ranges between 26% and 44%. These stark shifts in regional hydroclimatic variability are imprinted on the hydrogen and oxygen isotopic compositions of meteoric waters. Drier conditions favored lower LC-excess in rainfall (-17.3‰) and spring water (-6.5‰), whereas wetter conditions resulted in greater values (rainfall=+17.5‰; spring water=+10.7‰). The lower and higher LC-excess values in rainfall corresponded to the very strong 2014-16 El Niño and 2018 La Niña, respectively. During the recent triple-dip 2021-23 La Niña, LC-excess exhibited a significant and consistently increasing trend. These findings highlight the importance of combining hydrometric, synoptic, and isotopic monitoring as ENSO sentinels to advance our current understanding of ENSO impacts on hydrological systems across the humid Tropics. Such information is critical to constraining 21 st century projections of future water stress across this fragile region.

Nitrate contamination is affecting groundwater across the tropics. This study describes isotopic and ionic spatial trends across a tropical and volcanic multi-aquifer system in central Costa Rica in relation to land use change over four decades. Springs and wells (from 800 to 2,400 m asl) were sampled for NO3- and Cl- concentrations, δ18Owater, δ15NNO3, and δ18ONO3. A Bayesian isotope mixing model was used to estimate source contributions to the nitrate legacy in groundwater. Land use change was evaluated using satellite imagery from 1979 and 2019. The lower nitrate concentrations (< 1 mg/L) were reported in headwater springs near protected forested areas, while greater concentrations (up to ~63 mg/L) were reported in wells (mid- and low-elevation sites in the unconfined unit) and low-elevation springs. High-elevation springs were characterized by low Cl- concentrations and moderate NO3-/Cl- ratios, indicating the potential influence of soil nitrogen inputs. Wells and low-elevation springs exhibited greater NO3-/Cl- ratios and Cl- concentrations above 100 mg/L. A decreasing trend in NO3-/Cl- ratios coupled with greater Cl- values was also detected. Bayesian calculations suggest a mixture of sewage (domestic septic tanks), soil nitrogen (forested recharge areas), and chemical fertilizers (coffee plantations), as a direct result of abrupt land use change in the last 40 years. Our results confirm the incipient trend in increasing groundwater nitrogen and highlight the urgent need for a multi-municipal plan to transition from domestic septic tanks to regional sewage treatment and sustainable agricultural practices to prevent future groundwater quality degradation effectively.

Tracer-aided studies to understand source water partitioning in tropical ecosystems are limited. Here we report dry season source water partitioning in five unique ecosystems distributed across Costa Rica in altitudinal (<150-3,400 m asl) and latitudinal (Caribbean and Pacific slopes) gradients: evergreen and seasonal rainforests, cloud forest, Páramo, and dry forest. Soil and plant samples were collected during the dry season (2021). Plant and soil water extractions (triplicates) were conducted using controlled centrifugation. Stem water extraction efficiency and stem water content were calculated via gravimetric measurements. Water source contributions were estimated using a Bayesian mixing model. Isotope ratios in soil and stems exhibited a strong meteoric origin. Enrichment trends were detected mainly in stems and cactus samples within the dry forest ecosystem. Soil profiles revealed nearly uniform isotopic profiles; however, a depletion trend was observed in the Páramo ecosystem below 25 cm depth. More enriched compositions were reported in cactus samples for extracted water volumes above ~20% ( Adj. r2=0.34, p<0.01). The most prominent dry season water source in the evergreen rainforest (74.0%), seasonal rainforest (86.4%), and cloud forest (66.0%) corresponded with soil water. In the Páramo ecosystem, recent rainfall produced by trade wind incursions resulted in the most significant water source (61.9%), whereas in the dry forest, mean annual precipitation (38.6%) and baseflow (33.1%) were the dominant sources. The latter highlights the prevalence of distinct water uptake sources between recent cold front’s rainfall to more well-mixed soil moisture during the dry season.

There is still limited understanding of how waters mix, where waters come from and for how long they reside in tropical catchments. In this study, we used a tracer-aided model (TAM) and a gamma convolution integral model (GM) to assess runoff generation, mixing processes, water ages and transit times (TT) in the pristine humid tropical rainforest Quebrada Grande catchment in central Costa Rica. Models are based on a four-year data record (2016 to 2019) of continuous hydrometric and stable isotope observations. Both models agreed on a young water component of fewer than 95 days in age for 75% of the study period. The streamflow water ages ranged from around two months for wetter years (2017) and up to 9.5 months for drier (2019) years with a better agreement between the GM estimated TTs and TAM water ages for younger waters. Such short TTs and water ages result from high annual rainfall volumes even during drier years with 4,300 mm of annual precipitation (2019) indicating consistent quick near-surface runoff generation with limited mixing of waters and a supra-regional groundwater flow of likely unmeasured older waters. The TAM in addition to the GM allowed simulating streamflow (KGE > 0.78), suggesting an average groundwater contribution of less than 40% to streamflow. The model parameter uncertainty was constrained in calibration using stable water isotopes (δ2H), justifying the higher TAM model parameterization. We conclude that the multi-model analysis provided consistent water age estimates of a young water dominated catchment. This study represents an outlier compared to the globally predominant old water paradox, exhibiting a tropical rainforest catchment with higher new water fractions than older water.

Understanding how droughts propagate through the hydrological cycle from precipitation to streamflow and groundwater is important for improving water and risk management policies. At the catchment scale, the analysis of drought propagation and classification into drought types is usually done manually, which can be time consuming and difficult to replicate. Here, we developed an automated, objective procedure for classification of different drought types with the aim to study drought propagation in the tropics. The method was applied to the Savegre catchment in Costa Rica as a proof-of-concept. We first confirmed that drought events in the catchment could be classified into the process-based typology from the literature: classical rainfall deficit drought, wet-to-dry season drought, and composite drought. The automation algorithm was able to replicate the classification obtained with the manual typology with the exception of two events, and thus it is a development towards objective and time efficient hydrological drought analysis in tropical catchments. Most of the detected hydrological droughts (80% and 76% of all river discharge and baseflow droughts, respectively) were classical rainfall deficit droughts, which suggests that climate plays a more important role in drought development than catchment characteristics in this catchment. However, the importance of catchment characteristics was revealed by the presence of severe composite drought events and by the attenuation of significant precipitation droughts.

The páramos, a neotropical alpine grassland-peatland dominated biome of the northern Andes and Central America, play an essential role in regional and global cycles of water, carbon, and nutrients. They act as water towers, delivering water and ecosystem services mainly from the continental water divide at the Andean highland down to the Pacific and Amazon regions. The anthropogenic influence in form of the climate crisis exerts enormous pressures on these identified “hotspot” ecosystems and increases the vulnerability of nearby populations undermining the socio-economic and human development. Further, increasing pressures reduce the resilience to face climate shocks, and dramatically alters the hydro-climatic regime and shifts the páramos from long-term carbon sinks towards carbon sources. Despite their importance and vulnerability, only three decades ago, páramos, were globally among the least studied ecosystems. However, researchers have since identified them as ideal targets for solving water scarcity issues and to offset carbon emissions. Increasing awareness of the need for hydrological evidence to guide sustainable management of the páramos prompted action for generating data and to fill long-standing knowledge gaps. This has led to a remarkably successful community-research-policy effort to generate this knowledge. The combination of well-established and innovative approaches was used to data collection, processing and knowledge extraction. In this review, we provide a short overview of the state of knowledge of the hydrometeorology, flux dynamics, anthropogenic and the influence of extreme events in the regional páramos. Then, we present emerging technologies for hydrology and water resources research and management applied to páramos. Lastly, we discuss how converging science and policy efforts have leveraged traditional and new observational techniques to generate an evidence base that can support the sustainable management of the páramos. We conclude that this co-evolution of science policy was able to cover different spatial and temporal scales. Finally, we outline future research directions to showcase how sustainable long-term data collection can be sustained for the responsible development of páramo water towers.