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.

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.

Mountain System Recharge (MRS) processes are the natural recharge pathways in arid and semi-arid mountainous regions. However, MSR processes are often poorly understood and characterized in hydrologic models. Mountains are the primary source of water supply to valley aquifers via multiple pathways including lateral groundwater flow from the mountain block (Mountain-block Recharge, MBR) and focused recharge from mountain streams contributing to mountain front recharge (MFR) at the piedmont zone. Here, we present a multi-tool isogeochemical approach to characterize mountain flow paths and MSR processes in the northern Tulare basin, California. We used groundwater chemistry data to delineate hydrochemical facies and explain the chemical evolution of groundwater from the Sierra Nevada to the Central Valley aquifer. Isotope tracers helped to validate MSR processes. Novel application of End-Member Mixing Analysis (EMMA) using conservative chemical components revealed three MSR end-members: (1) evaporated Ca-HCO3 water type associated with MFR, (2) non-evaporated Ca-HCO3 and Na-HCO3 water types with short residence times associated with shallow MBR, and (3) Na-HCO3 groundwater type with long residence time associated with deep MBR. We quantified the contribution of each MSR process to the valley aquifer using mixing ratio calculation (MIX). Our results show that deep MBR is a significant component of recharge representing more than 50% of the valley groundwater. Greater hydraulic connectivity between the Sierra Nevada and Central Valley has significant implications for parameterizing Central Valley groundwater flow models and improving groundwater management. Our framework is useful for understanding MSR processes in other snow-dominated mountain watersheds.

Numerous socio-economic activities depend on the seasonal rainfall and groundwater recharge cycle across the Central American Isthmus. Population growth and unregulated land use changes resulted in extensive surface water pollution and a large dependency on groundwater resources. This work combines stable isotope variations in rainfall, surface water, and groundwater of Costa Rica, Nicaragua, El Salvador, and Honduras to develop a regionalized rainfall isoscape, isotopic lapse rates, spatial-temporal isotopic variations, and air mass back trajectories determining potential mean recharge elevations, moisture circulation patterns, and surface water-groundwater interactions. Intra-seasonal rainfall modes resulted in two isotopically depleted incursions (W-shaped isotopic pattern) during the wet season and two enriched pulses during the Mid-Summer Drought and the months of the strongest trade winds. Notable isotopic sub-cloud fractionation and near-surface secondary evaporation were identified as common denominators within the Central American Dry Corridor. Groundwater and surface water isotope ratios depicted the strong orographic separation into the Caribbean and Pacific domains, mainly induced by the governing moisture transport from the Caribbean Sea, complex rainfall producing systems across the N-S mountain range, and the subsequent mixing with local evapotranspiration, and, to a lesser degree, the eastern Pacific Ocean fluxes. Groundwater recharge was characterized by a) depleted recharge in highland areas (72.3%), b) rapid recharge via preferential flow paths (13.1%), and enriched recharge due to near-surface secondary fractionation (14.6%). Median recharge elevation ranged from 1,104 to 1,979 m asl. These results are intended to enhance forest conservation practices, inform water protection regulations, and facilitate water security and sustainability planning in the Central American Isthmus.