Groundwater recharge in highly-fractured volcanic aquifers remains poorly understood in the humid tropics, whereby rapid demographic growth and unregulated land use change are resulting in extensive surface water pollution and a large dependency on groundwater extraction. Here we present a multi-tracer approach including δ18O-δ2H, 3H/3He, and noble gases within the most prominent multi-aquifer system of central Costa Rica, with the objective to assess dominant groundwater recharge characteristics and age distributions. We sampled wells and large springs across an elevation gradient from 868 to 2,421 m asl. Our results suggest relatively young apparent ages ranging from 0.0±3.2 up to 76.6±9.9 years. Helium isotopes R/RA (0.99 to 5.4) indicate a dominant signal from the upper mantle across the aquifer. Potential recharge elevations ranged from ~1,400 to 2,650 m asl, with recharge temperatures varying from ~11°C to 19°C with a mean value of 14.5±1.9°C. Recharge estimates ranged from 129±78 to 1,605±196 mm/yr with a mean value of 642±117 mm/yr, representing 20.1±4.0% of the total mean annual rainfall as effective recharge. The shallow unconfined aquifer is characterised by young and rapidly infiltrating waters, whereas the deeper aquifer units have relatively older waters. These results are intended to guide the delineation and mapping of critical recharge areas in mountain headwaters to enhance water security and sustainability in the most important headwater dependent systems of Costa Rica.

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.