Modeling and assessing the impacts of climate change on groundwater recharge in endorheic basins of Northwest China

Climate change imposing additional stressors on groundwater resources globally, thereby predicting groundwater recharge (GR) changes is crucial to sustainably managing water resources, especially in the arid endorheic basins. Groundwater in the Endorheic Basins of Northwest China (NWEB) is potentially impacting regional socio-economic output and ecosystem stability due to the imbalance between supply and extraction exacerbated by climate change. Hence, recognizing the impacts of climate change on past and future GR is imperative for groundwater supply and sustainable groundwater management in the NWEB. Here, the impact of historical (1971-2020) and projected (2021-2100) climate changes on GR across the entire NWEB and three distinctive landscape regions (i.e., mountainous, oasis, and desert) were assessed. A coupled distributed hydrologic model (CWatM-HBV model), which integrates the Community Water Model (CWatM) and the HBV model, was run with three shared socioeconomic pathways (SSP1-2.6, SSP2-4.5, and SSP3-7.0) forcing from 10 general circulation models (GCMs) to simulate and analyze the interannual and seasonal variations of GR, along with their driving factors. Over the past 50 years, both precipitation and runoff have undergone significant increases, and leading to a dramatic rise in GR (0.09 mm yr-1). The future annual growth rate of GR is projected to range from 0.01 to 0.09 mm yr-1 from SSP1-2.6 to SSP3-7.0 across the entire NWEB, with the majority of the increase expected during the spring and summer seasons, driven by enhanced precipitation. GR from the mountainous region is the primary source (accounting for approximately 56-59 %) throughout the NWEB with the greatest increase anticipated. Precipitation and runoff have significant influences on GR in mountainous areas, and the impact of precipitation on GR is expected to increase over time. Changes in GR in oasis and desert areas are mainly limited by precipitation variation and increase in the SSP2-4.5 and SSP3-7.0 scenario. Additionally, the processes of glacial retreat and permafrost degradation will complicate the GR dynamics although the process is largely interfered with by anthropogenic environmental changes, especially in oasis-desert systems. The average annual recharge in the NWEB was 8.9 mm in the historical period and 13.6 ± 4.1 mm in the future. Despite an increase in GR due to climate change, groundwater storage is likely to continue to decline due to complex water demands in the NWEB. This study highlights the significance of future precipitation changes for GR and contributes to the understanding of the influence of climate change on groundwater systems and advances the sustainable management of water resources.

Evolution of groundwater hydrochemical characteristics and formation mechanism during groundwater recharge: A case study in the Hutuo River alluvial-pluvial fan, North China Plain

A pilot project for groundwater recharge from rivers is currently being carried out in North China Plain. To investigate the influence of river recharge on groundwater hydrochemical characteristics, dynamic monitoring and analysis of groundwater samples were conducted at a typical recharge site in the Hutuo River alluvial-pluvial fan in the North China Plain from 2019 to 2021. Hydrochemical, isotopic, and multivariate statistical analyses were used to systematically reveal the spatiotemporal variation of groundwater chemistry and its driving factors during groundwater recharge process. The results showed that the groundwater hydrochemical types and characteristics in different recharge areas and recharge periods exhibited obvious spatiotemporal differences. The groundwater type varied from HCO3·SO4-Na·Mg to HCO3·SO4-Ca·Mg in an upstream ecological area, while the groundwater type changed from SO4·HCO3-Mg·Ca to HCO3·SO4-Ca·Mg in the downstream impacted by reclaimed water. Changes in the contents of Ca2+, Mg2+ and HCO3- were mostly controlled by the water-rock interactions and mixing-dilution of recharge water, while the increases in Na+, NO3-, Cl-, SO42- and NO3- contents were mainly due to the infiltration of reclaimed water. Nitrogen and oxygen isotope (δ15N and δ18O) tests and the Bayesian isotope mixing model results further demonstrated that nitrate pollution mainly originated from anthropogenic sources, and the major contribution came from manure and sewage, with an average proportion of 64.6 %. Principal component analysis indicated that water-rock interactions, river-groundwater mixing and redox environment alternation were dominant factors controlling groundwater chemical evolution in groundwater recharge process.

Water, sanitation, and hygiene implications of large-scale recycling of treated municipal wastewater in semi-arid regions

Access to water, sanitation, and hygiene (WaSH) is crucial for national development, as it improves human health and fulfills a fundamental need. This study examines the impact of a large-scale groundwater (GW) recharge scheme using secondary treated wastewater (STW) on WaSH characteristics and identifies the major determinants of improved WaSH charecteristics in drought-hit regions of Kolar district, southern India. The study quantifies improved WaSH practices by comparing WaSH characteristics between impacted areas (influenced by STW) and non-impacted areas (not influenced by STW) of Kolar, using household survey data. Pearson's chi-square and student's t-test are used to verify differences between WaSH characteristics. Furthermore, a composite WaSH score is formulated, and a hierarchical stepwise multiple linear regression model is constructed to identify major determinants of improved WaSH scores. The results show that impacted areas have better WaSH characteristics, including daily water supply by gram panchayat, enhanced toilet uses among all family members, bathing patterns, cloth washing practices, toilet cleaning patterns, and water consumption per capita per day. The maximum and minimum WaSH scores of impacted areas were 17.50 and 6.50, respectively, while those of non-impacted areas were 14 and 4.5. This study finds that improved water availability, quality, and security due to daily water supply at the household level are the major determinants of improved WaSH practices. These results can inform policymakers in designing sanitation and hygiene improvement policies that integrate water recycling projects in drought-hit areas.

Active-distributed temperature sensing dataset beneath a braided river

Braided rivers play a significant role in replenishing groundwater, but our understanding of how these recharge rates fluctuate over time remains limited. Traditional techniques for gauging groundwater recharge are ineffective for studying complex braided river systems due to their insufficient spatiotemporal resolution. To address this gap, active-distributed temperature sensing (A-DTS) was used. This method combines fiber optic temperature measurements with an active heat source, enabling quantification of groundwater fluxes. In this study, twelve consecutive A-DTS surveys were conducted on a 100 m long hybrid fiber optic cable to a depth of 5 m beneath the Waikirikiri Selwyn River. This experiment was conducted during a period of relatively stable river stage and flow, highlighting the effectiveness of using A-DTS to measure temporal changes in groundwater recharge.

Numerical analysis of surface hydrogeological water budget to estimate unconfined aquifers recharge

Under changing climate, groundwater resources are the main drivers of socioeconomic development and ecosystem sustainability. This study assessed the contribution of two adjacent watersheds, Muse Street (MS) and West Wood (WW), with low and high urban development, to the Memphis aquifer recharge process in central Jackson, Tennessee, USA. The numerical MODFLOW model was created using data from 2017 to 2019 and calibrated using reported water budget components derived from in-situ data. The calibrated MODFLOW model was then used to investigate the impact of high and low urban developments on the recharge rate. The hydraulic parameters and recharge rates were optimized by adjusting the groundwater level based on the observed water level using PEST. The stochastic modeling was also carried out using the Latin Hypercube approach to reduce the uncertainty. The calibration results were satisfactory, with RMSE of 0.124 and 0.63 obtained in the WW and MS watersheds, respectively, indicating accurate estimation of the input parameters, precisely the hydrodynamic coefficients. The study results indicate that, per unit area, the MS watershed contributes 119% more to recharge and 186% more to riverbed leakage compared to the WW watershed. However, regarding total recharge and riverbed leakage, the WW watershed contributed more than the MS watershed. The results of this study have enhanced the knowledge of the impact of urbanization on hydrology and the recharge process in watersheds with diverse land uses.

Mechanism of groundwater recharge in the thick loess deposits by multiple environmental tracers

The loess-covered region accounts for ∼10 % of global land surface. Because of dry climate and thick vadose zones, the subsurface water flux is low but the water storage is relatively large . As a result, the groundwater recharge mechanism is complicated and currently controversial (e.g., piston flow or dual mode with piston and preferential flow). Taking typical tablelands in China's Loess Plateau as example study area, this study aims to qualitatively and quantitively evaluate the forms/rates and controls of groundwater recharge considering space and time. We collected 498 precipitation, soil water and groundwater samples in 2014-2021 for hydrochemical and isotopic analysis (Cl-, NO3-, δ18O, δ2H, 3H and 14C). A graphical method was employed to determine appropriate model to correct 14C age. Dual model exhibited in the recharge: regional-scale piston flow and local-scale preferential flow. Piston flow dominated groundwater recharge with a proportion of 77 %-89 %. Preferential flow gradually declined with increasing water table depths, and the upper depth limit may be <40 m. The dynamics of tracers proved that mixing and dispersion effects of aquifers limited the ability of tracers to capture preferential flow at short-time scales. Long-term average potential recharge (79 ± 49 mm/year) was close to actual recharge (85 ± 41 mm/year) at the regional scale, indicating the hydraulic equilibrium between unsaturated and saturated zones. The thickness of vadose zone controlled recharge forms, and precipitation dominated the potential and actual recharge rates. Land-use change can also affect the potential recharge rates at point and field scales but maintain the dominance of piston flow. The revealed spatially-varied recharge mechanism is useful for groundwater modeling and the method can be referred for studying recharge mechanism in thick aquifers.

Soil Conservation Service-Curve Number method-based historical analysis of long-term (1936-2016) temporal evolution of city-scale potential natural groundwater recharge from precipitation: case study Algiers (Algeria)

Managing groundwater resources in urban areas requires an adequate understanding and assessment of urban hydrogeological systems (structure, components, connections, and imposed conditions) as a part of a larger, dynamically evolving environment. Urbanization and climate change are amongst the widely recognized signs of such a continuous evolution. Within this context, the present study gives a quantitative assessment of the impact of these two factors threatening water resources in urban environments. The Soil Conservation Service-Curve Number (SCS-CN) method is used to conduct a long-term quantitative analysis of the temporal evolution of the potential natural groundwater recharge from precipitation at the scale of Algiers city for an 80-year-long period (1936-2016). The length of the study period allowed us to account for and analyze important changes in urban settings and climatic conditions within the study zone. Overall, two trend shifts over three distinct periods were found to characterize the temporal evolution of precipitation, several climate change indicators defined for the study, and the potential natural aquifer recharge. A strong, approximately 1:4, linear correlation between the estimated city-scale potential natural aquifer recharge and precipitation was observed for the studied period (R2 = 0.748). Moreover, even though the urban area has known a rapid (2nd order polynomial) increase from 1936 to 2016, climate change (accounted for via the changes in precipitation regime) impacted the city-scale potential natural groundwater recharge with higher magnitudes than urbanization. Finally, the computed climate change indicators show that starting in the mid-1980s, Algiers has started receiving less precipitations, with fewer heavy rain events and longer dry condition periods.

Identifying the suitable managed aquifer recharge (MAR) strategy in an overexploited and contaminated river basin

Managed aquifer recharge (MAR) is a promising adaptation measure to reduce vulnerability to climate change and hydrological variability. However, in areas where the basin is highly polluted, densely populated, and intensely cultivated, implementing suitable MAR strategies is a significant challenge. This study used a geographic information system-based multicriteria decision analysis (GIS-MCDA) approach to delineate the MAR potential sites using seven thematic layers describing surface and subsurface features. Further, basin-specific MAR approach was developed using information such as polluted water areas, canal network distribution for water supply, and cropping patterns. The results of this study indicate that only 17% of the area is highly suitable, while 54% and 29% were found moderately suitable and unsuitable for the MAR approach. Since most highly and moderately suitable sites were falling in the agricultural areas, agricultural-based MAR (AgMAR) was considered a preferred option. AquaCrop model for sugarcane was developed considering excess canal water supply during the grand growth stage to understand the AgMAR potential in the study area. It was observed that the potential recharge under normal irrigation scenarios varies from 135.5 to 272 mm/year, which can be increased through AgMAR up to 545 mm/year depending on the water availability for excess irrigations. This study provides an improved understanding of the parameters that should be considered for MAR site selection and post-GIS-MCDA analysis to assess the basin-specific MAR strategy.

Coupling hydrogeochemistry and stable isotopes (δ2H, δ18O and δ13C) to identify factors affecting arsenic enrichment of surface water and groundwater in Precambrian sedimentary rocks, eastern salt range, Punjab, Pakistan

The study area is a part of the Salt Range, where water quality is being deteriorated by natural and anthropogenic sources. This research integrates water quality assessment, arsenic enrichment, hydrogeochemical processes, groundwater recharge and carbon sources in aquifer. Total dissolved solid (TDS) contents in springs water, lake water and groundwater are in range of 681-847 mg/L, 2460-5051 mg/L and 513-7491 mg/L, respectively. The higher concentrations of magnesium and calcium in water bodies next to sodium are because of carbonates, sulfates, halite and silicates dissolution. The average concentrations of ions in groundwater are in order of HCO3- > SO42- > Cl- > Na+  > Mg2+ > Ca2+ > K+ > NO3-, virtually analogous to springs water, but different from lake water, categorized as poor quality and unfit for drinking purposes. Based on major ions hydrochemistry, NaCl and mixed Ca-Mg-Cl type hydrochemical facies are associated with concentration of arsenic (4.2-39.5 µg/L) in groundwater. Groundwater samples (70%) having arsenic concentration (11 ≤ As ≤ 39.5 µg/L) exceeded from World Health Organization (WHO) guideline (As ≤ 10 µg/L) in near neutral to slightly alkaline (6.7 ≤ pH ≤ 8.3), positive Eh(6 ≤ Eh ≤ 204 mV), signifying its oxic condition. Eh-pH diagrams for arsenic and iron indicate that 80% of groundwater for arsenic and iron were in compartments of HAsO42- and Fe(OH)3, unveil oxic environment. Arsenic is moderately positive correlated with TDS, sodium, chloride, bicarbonate, nitrate, sulfate and weak negative with δ13CDIC in surface and groundwater, forecasting multiple sources of arsenic to aquifer. Stable isotopes of waters show recharge of groundwater from local rain and lake water. The lower δ13CDIC values of groundwater are modified by influx of CO2 produced during biological oxidation of soil natural organic matter.

Impact of climate change on groundwater recharge in shallow young glacial aquifers in northern Poland

We investigated the influence of climate change in the period 1951-2020 on shallow aquifers in the Brda and Wda outwash plains (Pomeranian Region, Northern Poland). There was a significant temperature rise (0.3 °C/10 years), which accelerated after 1980 (0.66 °C/10 years). Precipitation became increasingly irregular - extremely rainy years occurred right after or before extremely dry years, and intensive rainfall events became more frequent after 2000. The groundwater level decreased over the last 20 years, even though the average annual precipitation was higher than in the previous 50 years. We carried out numerical simulations of water flow in representative soil profiles for the years 1970-2020 using the HYDRUS-1D model, developed and calibrated during our earlier work at an experimental site in the Brda outwash plain (Gumuła-Kawęcka et al., 2022). We used a relationship between the water head and flux at the bottom of the soil profiles (the third-type boundary condition) to reproduce groundwater table fluctuations caused by recharge variability in time. The calculated daily recharge showed a decreasing linear trend for the last 20 years (0.05-0.06 mm d^-1/10 years), and dropping trends in water table level and soil water content in the entire profile of vadose zone. Field tracer experiments were performed to estimate impact of extremely rain events on water flux in vadose zone. The results suggest that tracer travel times are strongly determined by water content in the unsaturated zone which is determined by precipitation amount in span of weeks, rather than extremely high precipitation events.

Field variation of groundwater recharge and its uncertainty via multiple tracers' method in deep loess vadose zone

Accurate estimation of groundwater recharge is a precondition for assessing its spatial variation at different scales, especially field scale. In the field, the limitations and uncertainties of different methods are first evaluated based on site-specific conditions. In this study, we evaluated field variation in groundwater recharge via multiple tracers in the deep vadose zone on the Chinese Loess Plateau. Five deep soil profiles (approximately 20 m deep) were collected in the field. Soil water content and particle compositions were measured to analyse soil variation, and soil water isotope (³H, ¹⁸O, and ²H) and anion (NO₃^- and Cl^-) profiles were used to estimate recharge rates. Distinct peaks in soil water isotope and nitrate profiles indicated a vertical one-dimensional water flow in the vadose zone. Although the soil water content and particle composition were moderately variable, no significant differences were observed in recharge rates among the five sites (p > 0.05) owing to the identical climate and land use. The recharge rates did not show a significant difference (p > 0.05) between different tracers' methods. However, recharge estimates by the chloride mass balance method indicated higher variations (23.5 %) than those by the peak depth method (11.2 % to 18.7 %) among five sites. Moreover, if considering the contribution of immobile water in vadose zone, groundwater recharge would be overestimated (25.4 % to 37.8 %) using the peak depth method. This study provides a favourable reference for accurate groundwater recharge and its variation evaluated using different tracers' methods in deep vadose zone.

Determination of groundwater origins and vulnerability based on multi-tracer investigations: New contributions from passive sampling and suspect screening approach

Knowledge about groundwater origins and their interactions with surface water is fundamental to assess their vulnerability. In this context, hydrochemical and isotopic tracers are useful tools to investigate water origins and mixing. More recent studies examined the relevance of contaminants of emerging concern (CECs) as co-tracers to distinguish sources contributing to groundwater bodies. However, these studies focused on known and targeted CECs a priori selected regarding their origin and/or concentrations. This study aimed to improve these multi-tracer approaches using passive sampling and qualitative suspect screening by exploring a larger variety of historical and emerging concern contaminants in combination with hydrochemistry and water molecule isotopes. With this objective, an in-situ study was conducted in a drinking water catchment area located in an alluvial aquifer recharged by several water sources (both surface and groundwater sources). CECs determined by passive sampling and suspect screening allowed to provide in-depth chemical fingerprints of groundwater bodies by enabling the investigation of >2500 compounds with an increased analytical sensitivity. Obtained cocktails of CECs were discriminating enough to be used as chemical tracer in combination with hydrochemical and isotopic tracers. In addition, the occurrence and type of CECs contributed to a better understanding of groundwater-surface water interactions and highlighted short-time hydrological processes. Furthermore, the use of passive sampling with suspect screening analysis of CECs lead to a more realistic assessment and mapping of groundwater vulnerability.

Geochemical and isotopic tracers to define the aquifer's vulnerability: the case study of the alluvial multi-aquifer system of the Friulian plain

The Friuli-Venezia Giulia Region (north of Italy) is characterized by the presence of high-quality freshwater resources which benefit local citizens, animals, environmental habitats, and also agriculture and production activities. Waters from wells, canal, and wastewater selected in the Fiume Veneto area, through a detailed lithological modeling, were sampled and analyzed to characterize them from a geochemical point of view. The chemical and isotopic characterization made it possible to establish provenance, and the average age of water used, making available the estimation of the relationships between recharge capacity and water use in the Fiume Vento area. The focus of this study is to define the average age of the resources based on the time required for the recharge contributions to compensate the losses induced by exploitation. The results made it possible to support the plans for a water balance using the provenance and average age of water sources for the protection of water reserves formed by the multi-aquifer system of the high and medium Friuli plain. The methodology applied has followed the legislation of the water directive considering the overexploitation due to unauthorized withdrawals of the sampling area.

Nitrogen fate during agricultural managed aquifer recharge: Linking plant response, hydrologic, and geochemical processes

Agricultural managed aquifer recharge (Ag-MAR, on-farm recharge), where farmland is flooded with excess surface water to intentionally recharge groundwater, has received increasing attention by policy makers and researchers in recent years. However, there remain concerns about the potential for Ag-MAR to exacerbate nitrate (NO₃^-) contamination of groundwater, and additional risks, such as greenhouse gas emissions and crop tolerance to prolonged flooding. Here, we conducted a large-scale, replicated winter groundwater recharge experiment to quantify the effect of Ag-MAR on soil N biogeochemical transformations, potential NO₃^- leaching to groundwater, soil physico-chemical conditions, and crop yield. The field experiment was conducted in two grapevine vineyards in the Central Valley of California, which were each flooded for 2 weeks and 4 weeks, respectively, with 1.31 and 1.32 m³ m^-2 of water. Hydrologic, geochemical, and microbial results indicate that NO₃^- leaching from the first 1 m of the vadose zone was the dominant N loss pathway during flooding. Based on pore water sample and N₂O emission data, denitrification played a lesser role in decreasing NO₃^- in the root zone but prolonged anoxic conditions resulted in a significant 29 % yield decrease in the 4-week flooded vineyard. The results from this research, combined with data from previous studies, are summarized in a new conceptual model for integrated water-N dynamics under Ag-MAR. The proposed model can be used to determine best Ag-MAR management practices.

Mapping potential groundwater accumulation zones for Karachi city using GIS and AHP techniques

Karachi is the largest industrial metropolitan of Pakistan facing an acute water shortage which is leading to an overdraft of groundwater resources in the city. Groundwater is an important freshwater resource for the city as millions of people depend for sustenance. However, over-exploitation of groundwater has led to decreased groundwater levels within the city leading to environmental issues of depleting aquifers, deteriorating groundwater quality, land subsidence, and harm to groundwater-dependent ecosystems. The objective of the study was to assess the potential groundwater accumulation zones by integrating hydrogeological aspects of the city through nine thematic layers using the Geographic Information System (GIS) based multi-criteria decision analysis (MCDA) technique. The potential groundwater accumulation map reveals that 20% of the area has a low potential, 70% has moderate potential, and around 10% of the area in the city is composed of a high potential accumulation zone. The upstream regions of the city have the highest recharge potential because of sandy soil and barren land use, which promote high infiltration rates. The urbanized downstream areas have the lowest recharge potential due to impervious fabric. The findings reveal that the MCDA technique can be used with confidence in data-scarce regions for groundwater resource assessment and management. The recharge potential map can help better manage groundwater resources in the city by helping explore groundwater extraction opportunities and could hint at areas suitable for artificial recharge wells/ponds.

Evaluating potential groundwater recharge in the unsteady state for deep-rooted afforestation in deep loess deposits

Groundwater recharge reduces due to high transpiration from shallow-rooted to deep-rooted afforestation. However, reaching a steady state in recharge process is challenging and no methods are available for assessing potential groundwater recharge under unsteady state. Hence, this study developed a new method to quantify groundwater recharge in the unsteady state by (1) calculating the water age (A₂) at maximum root depth (D₂) for deep-rooted afforestation using the chloride accumulative age method; (2) determining the soil depth (D₁) corresponding to A₂ under shallow-rooted vegetation using the multi-year average pore water velocity multiplied by A₂; (3) calculating the reduction in groundwater recharge (∆R) from shallow- to deep-rooted afforestation as the depth difference multiplied by the average water content between D₁ and D₂, divided by stand age. The average groundwater recharge for deep-rooted afforestation is equal to the average annual groundwater recharge under shallow-rooted vegetation minus ∆R. Soil cores with >25 m soil profiles below four land-use types of Hippophae rhamnoides Linn. (H. rhamnoides), Platycladus orientalis (L.) Franco (P. orientalis), Robinia pseudoacacia L. (R. pseudoacacia), and grassland were collected to measure soil water content, root distribution, and chloride and tritium contents. The results revealed that: (1) maximum root depths were 11.0 ± 0.5, 20.2 ± 1.2, and 22.6 ± 0.8 m, with soil water deficits of 373.48, 823.65, and 1847.92 mm under H. rhamnoides, P. orientalis, and R. pseudoacacia, respectively; (2) groundwater recharge following land-use change has not reached a steady state; (3) an average annual groundwater recharge was 89.12 mm yr^-1 under grassland, amounting to 16 % of the average annual precipitation; deep-rooted afforestation did not significantly differ, with 83.55, 84.91, and 85.65 mm yr^-1 under H. rhamnoides, P. orientalis, and R. pseudoacacia, respectively. This study contributes to a rational assessment of groundwater resources under unsteady state during land-use change.

Ground truthing global-scale model estimates of groundwater recharge across Africa

Groundwater is an essential resource for natural and human systems throughout the world and the rates at which aquifers are recharged constrain sustainable levels of consumption. However, recharge estimates from global-scale models regularly disagree with each other and are rarely compared to ground-based estimates. We compare long-term mean annual recharge and recharge ratio (annual recharge/annual precipitation) estimates from eight global models with over 100 ground-based estimates in Africa. We find model estimates of annual recharge and recharge ratio disagree significantly across most of Africa. Furthermore, similarity to ground-based estimates between models also varies considerably and inconsistently throughout the different landscapes of Africa. Models typically showed both positive and negative biases in most landscapes, which made it challenging to pinpoint how recharge prediction by global-scale models can be improved. However, global-scale models which reflected stronger climatic controls on their recharge estimates compared more favourably to ground-based estimates. Given this significant uncertainty in recharge estimates from current global-scale models, we stress that groundwater recharge prediction across Africa, for both research investigations and operational management, should not rely upon estimates from a single model but instead consider the distribution of estimates from different models. Our work will be of particular interest to decision makers and researchers who consider using such recharge outputs to make groundwater governance decisions or investigate groundwater security especially under the potential impact of climate change.

Integrated application of a Bayesian mixing model, numerical model, and environmental tracers to characterize groundwater recharge sources in a mountainous area

In this study, the combined use of a Bayesian mixing model (BMM), numerical model (random walk particle tracking-RWPT), and environmental tracers (δ¹⁸O-δD, ³H, and CFC) was applied to elucidate the probabilistic contribution of the recharge sources, flow path, and residence time of groundwater across the mountainous area of Jeju Island, South Korea. Especially, the BMM ability to estimate the variable recharge contributions to the aquifer by different elevations and seasons was investigated. The δ¹⁸O-δD isotopes showed that groundwater in the study area was primarily fed by precipitation during the wet season, and the BMM estimated that wet season recharge contributed to approximately 64% of the total. The BMM-based probabilistic estimation of recharge sources revealed a mixed contribution of source waters from different elevations. A notable difference in recharge flow path was observed between highland (>450 masl) and lowland (<400 masl) wells, where the inflow of source water from the regional flow was dominant in the former and both regional and local recharges served as significant groundwater sources in the latter. Evidence from age tracers (³H and CFC-12) also supported different recharge mechanisms between highland and lowland wells. A reasonable match between the BMM- and RWPT-derived recharge contributions (RMSE 0.02-0.06) was achieved within the uncertainty ranges, with RWPT being particularly useful for capturing different flow paths between highland and lowland wells. The dynamics revealed here provide important information for establishing an improved and informed groundwater management plan for the mountainous area of Jeju Island. Ultimately, this study highlights the advantageous integrated analysis of BMM, RWPT, and environmental tracer analyses to enhance the reliability of recharge area estimation and increase the collective understanding of complex hydrogeological systems in mountainous areas.

Influence of inter-aquifer leakage on well-injection capacity: Theory and aquifer-scale mapping for artificial recharge

Aquifer storage and recovery (ASR) is an important water resources management technique that involves the injection of a large volume of water underground. For the successful implementation of an ASR project, a target aquifer should have a sufficient injection capacity, which is the maximum volume of water that can be safely injected. In nature, no aquitard is perfectly impermeable, and inter-aquifer leakage may have a major impact on injection capacity. Despite the importance of determining the injection capacity for ASR planning, there is no quantitative methodology that estimates the injection capacity of leaky aquifers. In this study, we first develop a solution for injection capacity with inter-aquifer leakage based on the Hantush - Jacob solution, and conduct a comprehensive sensitivity analysis to elucidate the influence of inter-aquifer leakage on injection capacity. From the sensitivity analysis, we show that inter-aquifer leakage can impact injection capacity by more than one order of magnitude, depending on the hydrogeological and operational parameters. We then develop a practical mapping methodology that estimates the injection capacity of leaky aquifers. We demonstrate the proposed methodology by applying it to a potential ASR site in Minnesota, USA, where ASR is considered as a solution to alleviate groundwater contamination by PFAS chemicals. The case study results reveal significant spatial variability in injection capacity over the study area and show an average increase in the injection capacity of about 26% compared to that in the nonleaky scenario. We also analyze the uncertainty in the estimated injection capacity due to the variability of aquitard properties and show that the variability of aquitard vertical hydraulic conductivity leads to a larger uncertainty in the estimated injection capacity than does the variability of aquitard thickness. This study elucidates the effects of inter-aquifer leakage on injection capacity and provides a practical methodology for injection capacity mapping.

Determination of potential recharge zones and its validation against groundwater quality parameters through the application of GIS and remote sensing techniques in uMhlathuze catchment, KwaZulu-Natal, South Africa

Urbanization has accelerated changes in the quantity and stability of the water resources in the uMhlathuze watershed of KwaZulu-Natal. This study applied the use of GIS and remote sensing to demarcate groundwater potential recharge zones in the uMhlathuze catchment using AHP approach and Catastrophe theory by assigning weights to 10 parameters with their sub-criteria and the results were validated against groundwater quality data. It was discovered that 22.92% and 26.38% of the catchment is encompassed by 'Low' groundwater potential recharge zones, 0.37% and 0.08% by 'Very low' groundwater potential recharge zones, 9.42% and 10.26% by 'Good' groundwater potential recharge zones, 66.87% and 63.19% by 'Moderate', and 0.42% and 0.09% by 'Very good' for AHP and Catastrophe theory respectively. The resultant map demonstrated that recharge potential of groundwater is lowest in mountainous regions coupled with hard rock geology of low transmissivity, whereas the highest potential prevails in lower slopes and plains with more permeable soil. The findings of the validation revealed that the lowest and highest total dissolved solid, nitrate and groundwater level overlaps with the 'Good and Very good' and 'Low and Very low' groundwater potential recharge zones respectively. From the collective findings of this study, it is inferred that the convergence and use of GIS and remote sensing for delineating the groundwater potential recharge zones are effective. The study further recommends that this method can be applied in research/projects involving the implementation of artificial groundwater recharge structures for better groundwater planning and governance.

Formation mechanism of hydrogeochemical characterization of mineral water in Antu County, Changbai Mountain area

Antu County in the Changbai Mountains is an important source of mineral water, but there is a lack of research on the source of groundwater characteristic components, affecting the protection of water resources. This study obtained hydrochemical and isotopic data (28 groups in total, April and September in 2019) by summarizing research and sampling data in order to identify the formation process of characteristics. The formation mechanism of the characteristic components was revealed using geostatistical, isotopic, and hydrogeochemical inversion simulations. The results show that the metasilicic acid is a common component of groundwater water chemistry in the study area. The water body primarily receives stable recharge from low-mineralized precipitation with ages ranging from 27.7 to 38.4 years and recharge elevations ranging from 1160 to 2393 m, providing ample time for water-rock interaction. The dissolution of olivine, pyroxene, albite, and other siliceous minerals is the source of characteristic components, and deep faults and deep basalt heat flow are the key conditions for the formation of metasilicic acid. When low-mineralized precipitation recharges the underground aquifer, it dissolves the silica-aluminate and silicon-containing minerals in the surrounding rocks through the water-rock action under the effect of CO₂, causing a large amount of metasilic acid to dissolve into the groundwater and forming metasilic acid-type mineral water.

The impact of land use and land cover change on groundwater recharge in northwestern Bangladesh

Groundwater recharge is affected by various anthropogenic activities, land use and land cover (LULC) change among these. The long-term temporal and seasonal changes in LULC have a substantial influence on groundwater flow dynamics. Therefore, assessment of the impacts of LULC changes on recharge is necessary for the sustainable management of groundwater resources. The objective of this study is to examine the effects of LULC changes on groundwater recharge in the northwestern part of Bangladesh. Spatially distributed monthly groundwater recharge was simulated using a semi-physically based water balance model. Long-term temporal LULC change analysis was conducted using LULC maps from 2006 to 2016, while wet and dry LULC maps were used to examine seasonal changes. The results show that the impervious built-up area has increased by 80.3%, whereas vegetated land cover has decreased by 16.4% over the study period. As a result, groundwater recharge in 2016 has decreased compared to the level seen in 2006. However, the decrease in recharge due to long-term temporal LULC changes is very small at the basin scale (2.6 mm/year), although the impact on regional level is larger (17.1 mm/year) due to urbanization. Seasonal LULC variations also affect recharge due to the higher potential for dry seasonal LULC compared to the wet seasonal LULC, a substantial difference (20.6 mm/year). The results reveal important information about the groundwater system and its response to land cover changes in northwestern Bangladesh.

Impacts of climate and land use change on groundwater recharge under shared socioeconomic pathways: A case of Siem Reap, Cambodia

The rapid pace of urbanization blended with climate change has significantly altered surface and groundwater flows. In the context of tourism-driven economic potential areas, these drivers have greater effects, including threatening groundwater availability. This study assessed the combined impacts of climate and land use changes on the groundwater recharge (GWR) in Siem Reap, Cambodia utilizing Phase Six of the Coupled Model Intercomparison Project (CMIP6) global climate models (GCMs), DynaCLUE land-use model, and Soil Water Assessment Tool (SWAT). Three climate models CanESM5, EC_Earth3, and MIROC6, out of seven, best captured the observed data after performance evaluation through the entropy method, were bias-corrected linearly for two shared socioeconomic pathways (SSPs) - SSP2-4.5 and SSP5-8.5. The results indicate a general increase in precipitation under both SSPs, while the average annual maximum temperature is likely to increase by 0.024 °C/year and 0.049 °C/year under SSP2-4.5 and SSP5-8.5, respectively. A similar trend but relatively higher increase is expected for the minimum temperature. Furthermore, the historical land use change showed the expansion of urban settlement by 373% between 2004 and 2019 at the expense of forest and shrubland. Future land use projections from the DynaCLUE model show that the urban settlements in the study area are likely to expand, from their 2019 condition, by 55% in 2030, 209% in 2060, and 369% in 2090 under SSP2 and at double of these rates under SSP5 scenario. The GWR is expected to rise by 39-53% during the wet season and decrease by 13-29% during the dry season under both scenarios. Meanwhile, under constant land use, the GWR is likely to increase more compared to other scenarios, highlighting the importance of land use planning to policymakers and planners. Additionally, the study shall also be important to practitioners and researchers in understanding, planning, and evaluating the performance of multiple climate models in groundwater assessment.

Surface water and groundwater interaction in the Kosi River alluvial fan of the Himalayan Foreland

We report the isotopic composition of the surface water and groundwater of the Kosi River fan on the Himalayan Foreland, India. We have collected 65 water samples from surface water (Kosi River (n = 2), streams (n = 9), waterlogging (n = 29), and canal (n = 4)), and groundwater (n = 21) for δ¹⁸O and δ²H analysis during December 2019. We obtained groundwater level data measured at the observation wells from the Central Groundwater Board, India, for 1996 and 2017. The groundwater level varies from 1.0 to 8.1 m below ground level (bgl) and from 0.5 to 9.0 m bgl during 1996 and 2017, respectively. We have used water table fluctuation approach to estimate the recharge rate. The recharge rate in the Kosi Fan varies from 0.7 to 21.4 mm/year from 1996 to 2017. Further, we have used δ¹⁸O and δ²H values of water samples to identify the source and the interaction between surface water and groundwater. The δ¹⁸O value of groundwater shows a wide variation (from -9.3‰ to -5.6‰) compared to the surface water, i.e., streams (-7.8‰ to -6.4‰) and canals (-6.9‰ to -6.0‰), suggesting mixing in groundwater during recharge processes. Furthermore, we have used a two-component mixing model to assess the fraction contribution from streams and precipitation to groundwater. The estimated fraction contribution from stream water to groundwater ranges from 45 to 83%. We also suggest higher recharge is limited up to the depth of 6 m bgl. We suggest precipitation and surface water actively recharge groundwater. We conclude that marked spatial variation in the isotopic composition of groundwater is mainly due to the local recharge sources and interaction between surface water and groundwater.

Developing a pan-European high-resolution groundwater recharge map - Combining satellite data and national survey data using machine learning

Groundwater recharge quantification is essential for sustainable groundwater resources management, but typically limited to local and regional scale estimates. A high-resolution (1 km × 1 km) dataset consisting of long-term average actual evapotranspiration, effective precipitation, a groundwater recharge coefficient, and the resulting groundwater recharge map has been created for all of Europe using a variety of pan-European and seven national gridded datasets. As an initial step, the approach developed for continental scale mapping consists of a merged estimate of actual evapotranspiration originating from satellite data and the vegetation controlled Budyko approach to subsequently estimate effective precipitation. Secondly, a machine learning model based on the Random Forest regressor was developed for mapping groundwater recharge coefficients, using a range of covariates related to geology, soil, topography and climate. A common feature of the approach is the validation and training against effective precipitation, recharge coefficients and groundwater recharge from seven national gridded datasets covering the UK, Ireland, Finland, Denmark, the Netherlands, France and Spain, representing a wide range of climatic and hydrogeological conditions across Europe. The groundwater recharge map provides harmonised high-resolution estimates across Europe and locally relevant estimates for areas where this information is otherwise not available, while being consistent with the existing national gridded datasets. The Pan-European groundwater recharge pattern compares well with results from the global hydrological model PCR-GLOBWB 2. At country scale, the results were compared to a German recharge map showing great similarity. The full dataset of long-term average actual evapotranspiration, effective precipitation, recharge coefficients and groundwater recharge is available through the EuroGeoSurveys' open access European Geological Data Infrastructure (EGDI).

Australian gridded chloride deposition-rate dataset

Chloride deposition-rate measurements at points within Australia are upscaled to the entire continent on a regular 0.05° grid. The upscaling uses a double-exponential correlation between deposition rate and distance to the coast, where the parameters in the double-exponential are spatially varying. These parameters are estimated using least-squares with Tikhonov regularisation to ensure minimal spatial variability. A calibration-constrained, null-space Monte-Carlo analysis is used to quantify uncertainty in the prediction. The resulting dataset consists of the best-fit chloride deposition rates across Australia, as well as estimates of uncertainty. The dataset can be used for various purposes including: estimating groundwater recharge through the use of the chloride mass-balance method; catchment salt balance estimates; regional investigations of groundwater hydrochemistry; and, corrosion prediction.

Assessment of the spatial-temporal distribution of groundwater recharge in data-scarce large-scale African river basin

The systematic assessment of spatial and temporal distribution of groundwater recharge (GWR) is crucial for the sustainable management of the water resources systems, especially in large-scale river basins. This helps in identifying critical zones in which GWR largely varies and thus leads to negative consequences. However, such analyses might not be possible when the models require detailed hydro-climate and hydrogeological data in data-scarce regions. Hence, this calls for alternate suitable modeling approaches that are applicable with the limited data and, however, includes the detailed assessment of the spatial-temporal distribution of different water balance components especially the GWR component. This paper aimed at investigating the spatial and temporal distribution of the GWR at monthly, seasonal and annual scales using the WetSpass-M physically distributed hydrological model, which is not requiring the detailed catchment information. In addition, the study conducted the sensitivity analysis of model parameters to assess the significant variation of GWR. The large-scale river basins such as the Omo river basin, Ethiopia, were chosen to demonstrate the potential of the WetSpass-M model under limited data conditions. From the modeling results, it was found that the maximum average monthly GWR of 13.4 mm occurs in July. The estimated average seasonal GWR is 32.5 mm/yr and 47.6 mm/yr in the summer and winter seasons, respectively. Further, it was found that GWR is highly sensitive to the parameter such as average rainfall intensity factor.

Processes explaining the origin and evolution of groundwater composition in the Andean Precordillera and Altiplano of the Tarapacá Region of northern Chile

In the arid area of northern Chile, groundwater resources in the Andean formations are essential for native populations, ecological services, mining, and other human activities. Validated conceptual hydrogeological models are required for current and future water and land management. This work aims to explain the processes controlling the origin and distribution of recharge and groundwater composition in the Andean Precordillera and Altiplano of the Tarapacá Region of northern Chile, using major solutes in spring, river, and well water, and the stable and radioactive isotopes of water oxygen, hydrogen, and dissolved inorganic carbon. The waters are mainly of the Na-Ca-SO₄ type. Processes controlling the chemical evolution of waters are atmospheric dust contribution, evapo-concentration, and enhanced volcanic rock weathering, as well as halite dissolution in some locations. The isotopic composition of Precordillera eastern flank water samples follows an evaporation line, while those in the western flank, in the Altiplano, follow a line that is parallel to the local meteoric line, suggesting unsaturated zone evaporation processes of infiltrated rainfall. δ¹³C_DIC contents (-2 to -27‰) indicate mixing processes, volcanic CO₂ in the Altiplano, and calcite dissolution in some sectors. In the western depression, the only recharge is due to water infiltration in creek channels. In the highland areas, 5-25% of precipitation produces recharge. The estimated groundwater renewal time in the Precordillera was 3-14 kyr. The piezometric elevation in the Precordillera due to low-permeability intrusive rocks and local recharge prevents the east-west groundwater transfer from the Altiplano to the western depression and explains why the volcanic CO₂ in the Altiplano basins is not observed on the western flank. These results provide new insights for the evolution of water quality in volcanic aquifers in arid environments and provide considerations for estimating groundwater residence times using radiocarbon in areas influenced by volcanic CO₂.

Appraisal of groundwater recharge in Neelum watershed (Upper Indus Basin) using geospatial water balance technique

Water availability is important for survival of millions of people living in the Himalayan region of Upper Indus Basin and adequate monitoring system is for better water resources management. In the present study, groundwater recharge appraisals in the Neelum watershed (Upper Indus Basin) were investigated by using water balance and geospatial modeling techniques on monthly time-scale climate data from 1989 to 2015. Results demonstrated that on an average out of total annual rainfall (i.e., 2028 mm), about 46% of the rainfall convert to surface runoff and 35% loss to atmosphere via evapo-transpiration (ET), while the remaining 18% contribute to infiltrate the groundwater recharge. Groundwater recharge enhanced during snow-melt from December to March and the rainfall infiltration increased during July and August months. Similarly, the infiltration ranges 106-177 mm from January to March and 45-51 mm from December to July. The groundwater discharge in the form of oozing from the spring occurred during the remaining six months, which ultimately contributed to the baseflow of the stream. Findings from the study revealed variations in groundwater recharge during the years and hence recommended more hydrological studies to predict future changes in climate and land use for sustainable development of freshwater resources in the Upper Indus Basin.

Salix psammophila afforestations can cause a decline of the water table, prevent groundwater recharge and reduce effective infiltration

Afforestation can reduce desertification and soil erosion. However, the hydrologic implications of afforestation are not well investigated, especially in arid and semi-arid regions. China has the largest area of afforestation in the world, with one-third of the world's total plantation forests. How the shrubs affect evapotranspiration, soil moisture dynamics, and groundwater recharge remains unclear. We designed two pairs of lysimeters, one being 1.2 m deep and the other one 4.2 m deep. Each pair consists of one lysimeter with bare soil, while on the other one a shrub is planted. The different water table depths were implemented to understand how depth to groundwater affects soil moisture and water table dynamics under different hydrological conditions. Soil moisture, water table depth, sap flow, and rainfall were measured concurrently. Our study confirms that for the current meteorological conditions in the Ordos plateau recharge is reduced or even prohibited through the large-scale plantation Salix psammophila. Shrubs also raise the threshold of precipitation required to increase soil moisture of the surface ground. For the conditions we analyzed, a minimum of 6 mm of precipitation was required for infiltration processes to commence. In addition to the hydrological analysis, the density of root distribution is assessed outside of the lysimeters for different water table depths. The results suggest that the root-density distribution is strongly affected by water table depth. Our results have important implications for the determination of the optimal shrub-density in future plantations, as well as for the conceptualization of plant roots in upcoming numerical models.

Managed aquifer recharge implementation criteria to achieve water sustainability

Depletion of groundwater is accelerated due to an increase in water demand for applications in urbanized areas, agriculture sectors, and energy extraction, and dwindling surface water during changing climate. Managed aquifer recharge (MAR) is one of the several methods that can help achieve long-term water sustainability by increasing the natural recharge of groundwater reservoirs with water from non-traditional supplies such as excess surface water, stormwater, and treated wastewater. Despite the multiple benefits of MAR, the wide-scale implementation of MAR is lacking, partly because of challenges to select the location for MAR implementation and identify the MAR type based on site conditions and needs. In this review, we provide an overview of MAR types with a basic framework to select and implement specific MAR at a site based on water availability and quality, land use, source type, soil, and aquifer properties. Our analysis of 1127 MAR projects shows that MAR has been predominantly implemented in sites with sandy clay loam soil (soil group C) and with access to river water for recharge. Spatial analysis reveals that many regions with depleting water storage have opportunities to implement MAR projects. Analyzing data from 34 studies where stormwater was used for recharge, we show that MAR can remove dissolved organic carbon, most metals, E. coli but not efficient at removing most trace organics, and enterococci. Removal efficiency depends on the type of MAR. In the end, we highlight potential challenges for implementing MAR at a site and additional benefits such as minimizing land subsidence, flood risk, augmenting low dry-season flow, and minimizing salt-water intrusion. These results could help identify locations in the water-stressed regions to implement specific MAR for water sustainability.

Shallow urban aquifers under hyper-recharge equatorial conditions and strong anthropogenic constrains. Implications in terms of groundwater resources potential and integrated water resources management strategies

Humid equatorial regions are recognized as the least documented in term of hydrogeological functioning of aquifers despite the fact that they house a lot of developing countries and that groundwater is often the main water resource. Regarding this aspect, a study was conducted in sub-Saharan Africa, focusing on the Mio-Pliocene aquifer in Douala megacity (Cameroon) which is the rainiest city in West-Africa (about 4000 mm/year) with one of the greatest demographic growth rate of the African continent. Firstly, groundwater recharge rate has been calculated through water balance and Water Table Fluctuation methods. Results show that the aquifer is characterized by a high recharge of 600-760 mm/year. Then infiltration process and groundwater flow conditions have been examined by combining hydrogeological and isotopic methods. Rainwater infiltrated is recycled in the vadose zone through plants roots transpiration and groundwater flows with a Darcy velocity of 5 m/day. From the recharge area to the estuary, the mineralization increases controlled by anthropogenic activities and water-rocks interactions which are amplified by the residence time and accelerated by the hot and humid climate of Douala. The paper ends with the determination of natural background levels (NBLs) and threshold values (TV) of chemical components in groundwater to assess the contamination for different flow paths. This multi-proxy study and the establishment of NBLs and TV can be beneficial to improve groundwater resources management. Moreover, the conceptual model provided in this study could be used as a reference for porous aquifers submitted to high rainfall amount.

Effects of land use on groundwater recharge of a loess terrace under long-term irrigation

Understanding the relation between land use and groundwater recharge is of direct interest in areas that are prone to geohazards. This study was performed to characterize the effects of land use on groundwater table (GWT) distribution in the Heifangtai (HFT) loess terrace under long-term flood irrigation, as irrigation intensity was governed by the crop types. In light of the regional temperature-vegetable dryness index (TVDI) and optical images obtained in different seasons over four years, the dataset incorporated the growth cycles of the local crops, which in turn improved the workload and accuracy of land use detection. Irrigation intensities for different crops were recorded during field investigations. A total of 26 electrical resistivity tomography (ERT) profiles were conducted to estimate the GWT distribution of HFT terrace, which was further utilized to reveal the relation between land use (i.e. irrigation) and GWT fluctuations. The results indicated that high local GWT was associated with vegetable fields which had the highest irrigation intensity, and by contrast, low local GWT was attributed to the perennial fruit tree fields with the lowest irrigation intensity. The discharge was monitored intermittently for over three years to analyze the effect on the local GWT close to margin of HFT. The result showed a strong correlation between spring discharge and GWT fluctuations along the margin of HFT, indicating that groundwater recharge was affected by both irrigation and spring discharges. These results suggested that incorporating growth cycles of crops can effectively facilitate the interpretation of remote sensing data for land use detection, and the results can provide useful guidance for improving irrigation programs and for alleviating geohazards in regions that are under long-term irrigation.

Groundwater recharge and water table response to changing conditions for aquifers at different physiography: The case of a semi-humid river catchment, northwestern highlands of Ethiopia

Groundwater recharge estimation, aquifer response to meteorological variables, and evapotranspiration calculations have been performed on a semi-humid catchment, in northwestern Ethiopian plateau. The Soil Moisture Balance (SMB), WetSpass water balance model, Water Table Fluctuation (WTF), and Chloride Mass Balance (CMB) methods are applied to estimate the groundwater recharge. Accordingly, 431 mm, 462 mm, and 477 mm recharge amounts are estimated as mean annual value, respectively, using SMB, WetSpass, and CMB methods. Based on the WTF method, the annual recharge rates of the volcanic aquifers range from 157 mm to 760 mm. The SMB and WetSpass methods are less effective for the flat physiographic area, where the recharge rate is storage controlled rather than precipitation amount. The calculated high recharge for maintain-front aquifers using WTF is attributed to extra rising due to lateral groundwater flow, which restricts the reliability of the method for such aquifer geometries. High groundwater level rising rate (121 mm/day) has been observed for the steeply sloping, low rates (11 mm/day) for the flat floodplain, and intermediate rate (52 mm/day) for the gently sloping volcanic aquifers. Similarly, receding rates of 3.18 mm/day were found for the steeply sloping, 0.40 mm/day for the floodplain, and 1.14 mm/day for the gentle sloping aquifers. The recession, in all of the topographies, is happening with second-order polynomial decay function. A strong connection between the shallow and deep groundwater aquifers is noted. Storage change in the relatively deeper volcanic aquifers is due to vertical groundwater flow from the overlying alluvial aquifer. This indicates that the recharge mechanism is local, and may be the reason for the low aquifer productivity of the Dangila wellfield. Diurnal water table fluctuation is detrended from the receding trend of the dry period, and evapotranspiration from the groundwater is estimated at 28% of total ET.

Hydrological response of karst stream to precipitation variation recognized through the quantitative separation of runoff components

Various water transmitting media are related to highly variable water source compositions, which limit the understanding of the aquifer structure and hydrological processes in a karst catchment. This study aims to understand the variation in water contribution by matrix, fissure, and conduit flows during storm and seasonal scales based on discharge, electrical conductivity (EC), and nitrate measurements of stream water in a typical dolomite catchment during 2017-2018 and discusses the hydrological response mechanism of a karst aquifer to rainfall characteristics. Time-series analyses of discharge and EC indicated that the rapid response time (mean lag time < 1 h) was mainly controlled by rainfall intensity, and the lag time decreased significantly when the rainfall intensity was lower than 15 mm/h. However, the mean discharge was dominated by the rainfall amount and antecedent moisture state. Hydrograph separation based on nitrate indicates that the contribution of soil water was irrelevant for recharging the stream during a non-rain period, whereas epikarst water contributed more than 83.2% of the total flow during a rainfall event. As indicated by the EC frequency distribution analyses, the contribution ratios of the surface, conduit, fissure, and matrix flows were 1:1.8:2.1:7.1, 1:1.6:5.3:6.3, and 0:0:0:1 during stormy, heavy, and light rainfall events, respectively. These parameters indicate that the degree of karstification was low in the karst aquifer. Seasonal frequency distribution analyses of EC indicate that higher rainfall amounts and rainfall intensities during the wet season promoted the contribution of conduit flow to approximately 11.4% of the total flow; however, matrix flow dominated the recharge of the streamflow and its contribution was more than 55.6% during each season. Our results suggest that the permeability of the epikarst matrix dominates the storage and transfer functions in dolomite karst aquifers with low karstification.

Managed Aquifer Recharge in the Gulf Countries: A Review and Selection Criteria

The Gulf Cooperation Council (GCC) countries are arid with very limited availability of water resources. In recent years, these countries have started an intensive program to increase the storage of groundwater through various techniques of managed aquifer recharge (MAR). Water consisting of varying quantity and quality (derived from various sources) are used via MAR techniques to increase the groundwater storage and, if possible to enhance its quality, respectively. This paper presents a review of the MAR techniques practiced in GCC countries including the implementation strategies of the different structures. Generally, seven MAR techniques are utilized in GCC countries including dams, aquifer storage and recovery (ASR) technique, aquifer storage transfer and recovery (ASTR) technique, ponds, soil aquifer treatment (SAT) technique, rooftop rainwater harvesting, and Karez/Ain system. Results indicated that ASR using excess desalinated water or treated sewage effluent (TSE) is the most used MAR technique in GCC countries, followed by the use of ASTR, dams, and ponds. Based on this review, twelve different selection criteria have been developed for GCC countries for better MAR practice in the future.

Occurrence and risk assessment of fluoroquinolone antibiotics in reclaimed water and receiving groundwater with different replenishment pathways

Artificial recharge to groundwater with reclaimed water is considered a promising method to alleviate groundwater depletion and over-exploitation. However, the occurrence of fluoroquinolone antibiotics (FQs) was ubiquitous in wastewater, surface water, groundwater and even drinking water threating human health and ecology. In this study, the occurrence of six selected FQs in reclaimed water effluent and their removal by tertiary treatment units were investigated. The overall removal efficiencies in average of the tertiary treatment processes in Beijing and Changzhou were ranging from 21.2% to 55.2%. Activated carbon exhibited better performance for FQs removal than ozone and biological treatment such as membrane bioreactor, anaerobic-anoxic-oxic and biofilter. The results of two pilot study showed that the impact of reclaimed water to groundwater quality in terms of FQs concentration by direct injection in GBD was stronger than surface spreading in Changzhou, which might be due to the recharge strategy and the physical and chemical characteristics of sediment and aquifer soil. The hazard quotient (HQ) values of ofloxacin (OFL) in reclaimed water was up to 12.54, indicating the extreme eco-toxicological risk, while enrofloxacin (ENR) exhibited medium risk. After recharge with reclaimed water, the HQ values of OFL and ENR in groundwater ranged from low to medium ecological risk to the environment. Thus, the FQs in reclaimed water need to be paid more attention during their reuse for groundwater recharge, especially by direct injection. It is suggested that FQs should be considered in the priority substances lists in standards and guidelines of reclaimed water reuse for groundwater recharge to ensure the safety of groundwater.

Identifying groundwater end members by spatio-temporal isotopic and hydrogeochemical records

Intensive groundwater use has altered the local hydrological cycle within the Bajío Guanajuatense, Mexico. To improve the knowledge of this hydrogeological system and support water management in the area, groundwater end members were identified using multivariate statistical analysis. Pumped groundwater is composed of two well-mixed end members: (a) recent recharge, affected by a reuse cycle through irrigation where nitrate and chloride evolve and reach levels of 368 mg/L and greater than 100 mg/L, respectively, and (b) deep old groundwater. Mixing estimations show that most wells extract at least 70% of deep groundwater, and some of them extract more than 94%, posing a development and groundwater sustainability conundrum in the area.

Modelling the hydrologic effects of vegetation growth on the long-term trajectory of a reclamation watershed

Reclamation watersheds that integrate fen peatlands into the design require the inclusion of uplands that are capable of supporting forest development while concurrently supplying sufficient groundwater recharge to downgradient wetland ecosystems. This necessitates selecting materials with suitable soil hydraulic properties and identifying the appropriate thickness and layering to fulfill the dual function of uplands as water storage, and water conveyance features. Currently, these systems incorporate tailings sand - a mine waste material - overlain by a cover soil of fine forest-floor material. The developmental pathway of these uplands is currently unknown, and it is unclear whether these landforms will provide enough groundwater recharge once a climax vegetation community establishes. Therefore, this research attempts to estimate the maximum density of vegetation, and associated water balance fluxes of a constructed upland integrated into a peatland watershed. The numerical modelling software HYDRUS-1D simulated soil moisture dynamics using a 65-year meteorological record, and a plant water stress algorithm was used to estimate the maximum sustainable leaf area index that the upland could support. Based on the thickness of the cover soil, the upland could support an average leaf area index of 1.2. Under this vegetation density, average annual groundwater recharge was 83 mm, and predominantly supplied by snowmelt (64%). Given this quantity of recharge, the model indicates that the upland will continue to provide enough groundwater to offset the anticipated water deficit in the downgradient fen ecosystem. However, by altering the design of the upland, specifically the spatial arrangement and thickness of cover soil, the same recharge could be supplied while also allowing for a higher average vegetation density. Such a design could allow for the creation of watersheds with a higher proportion of peatland.

Feasibility of rainwater harvesting for sustainable water management in urban areas of Egypt

Egypt's limited water resources, rapid population growth, and climate change are increasing the gap between water demand and supply. Meanwhile, significant amounts of rain fall in some regions in Egypt during specific storm events, which in some cases, lead to disasters like flash floods and inundations. Rainwater harvesting (RWH) can be considered as a sustainable promising solution to water shortage and inundation problems. In this work, the feasibility of RWH for urban areas was assessed over 22 cities throughout Egypt. Results show that the annual volume of rainwater harvested can reach 142.5 MCM in the considered cities, provided that all rain falling on the urban areas is collected. High potential of rainfall harvesting was found for cities that located on the North Coast, e.g., the potential water saving from the share of RWH in Alexandria can satisfy around 12% of its future supplementary domestic water needs. In contrast, rainfall over the cities located on the middle and the south of the country is insignificant to be harvested. A case study for the 5th settlement region in Cairo was discussed in terms of groundwater recharge and surface runoff estimation for two conditions: No-RWH and RWH systems by implementing recharge wells to store rainwater into the aquifer. Land cover classification maps of urban areas were created by using the ARCGIS software to estimate equivalent infiltration coefficients. The results demonstrate that the implementation of such RWH system has a significant impact on the regional water cycle, where the effective infiltration coefficient increased from 10% (No-RWH) to 75% (RWH) in the case study. Accordingly, the runoff coefficient decreased in the case study from 0.8 (No-RWH) to 0.15 (RWH), and the volume of runoff decreased in the case of RWH by around 82% lower than that of the No-RWH condition. Thus, direct infiltration of RWH into an aquifer can play an important role in sound water management for urban environments, as this may lead to a significant reduction in risks of flooding and expenses of municipal drainage systems installation and operation.

Calculating man-made depletion of a stressed multiple aquifer resource on a national scale

An inexorable depletion of groundwater occurs where groundwater abstraction exceeds the natural recharge, a typical state of (semi-)arid regions, which calls for sustainable management of groundwater resources. This study aims to assess the available storage and recharge rates on a national scale in time and space by modelling the natural recharge in combination with a method to evaluate changing groundwater volumes, which revealed measures to quantify the overdraft of the observed national groundwater resources in Jordan. Applying the combination of hydrological model and method to evaluate changing groundwater volumes, a climate-driven systematic decline of groundwater recharge was eliminated as responsible process, while overdraft leads to dropping groundwater tables. The major findings are, the intensity of groundwater abstraction from a basin becomes visible through the fact, that simulated baseflow exceeds significantly the observed baseflow. About 75% of Jordan's groundwater basins are subject to intense groundwater depletion, reaching annual rates of up to 1 m in some basins. The most affected areas are the basins Zarka, Azraq and the predominantly fossil groundwater reservoirs in Southern Jordan. Contrasting the past, when variable annual precipitation patterns did not negatively influence groundwater recharge, simulations show significantly reduced annual groundwater recharge all over Jordan. Particularly affected is the agricultural backbone in the Jordan Mountains, where recharge rates are predicted to vary between -30 mm/yr and +10 mm/yr in the coming decades, being reflected in the disappearance of freshwater springs and ascending saltwater. The applied methodology is relevant and transferable to other data- and water scarce areas worldwide, allowing (i) a fast estimation of groundwater reservoir development on a national scale and (ii) an investigation of long-term effects of overdraft.

Using machine learning algorithms to map the groundwater recharge potential zones

Groundwater recharge is indispensable for the sustainable management of freshwater resources, especially in the arid regions. Here we address some of the important aspects of groundwater recharge through machine learning algorithms (MLAs). Three MLAs including, SVM, MARS, and RF were validated for higher prediction accuracies in generating groundwater recharge potential maps (GRPMs). Accordingly, soil permeability samples were prepared and are arbitrarily grouped into training (70%) and validation (30%) samples. The GRPMs are generated using sixteen effective factors, such as elevation (denoted using a digital elevation model; DEM), aspect, slope angle, TWI (topographic wetness index), fault density, MRVBF (multiresolution index of valley bottom flatness), rainfall, lithology, land use, drainage density, distance from rivers, distance from faults, annual ETP (evapo-transpiration), minimum temperature, maximum temperature, and rainfall 24-hr. Subsequently, the VI (variables importance) is assessed based on the LASSO algorithm. The GRPMs of three MLAs were validated using the ROC-AUC (receiver operating characteristic-area under curve) and various techniques including true positive rate (TPR), false positive rate (FPR), F-measures, fallout, sensitivity, specificity, true skill statistics (TSS), and corrected classified instances (CCI). Based on the validation, the RF algorithm performed better (AUC = 0.987) than the SVM (AUC = 0.963) and the MARS algorithm (AUC = 0.962). Furthermore, the accuracy of these MLAs are included in excellent class, based on the ROC curve threshold. Our case study shows that the GRPMs are potential guidelines for decision-makers in drafting policies related to the sustainable management of the groundwater resources.

Can nature-based solutions contribute to water security in Bhopal?

Bhojtal, a large man-made lake bordering the city of Bhopal (Madhya Pradesh state, central India), is important for the city's water supply, connoted the lifeline of the city. Despite the dry though not arid and markedly seasonal climate, soil impermeability hampers infiltration into the complex geology underlying the Bhojtal catchment. Rural communities in the catchment are nonetheless high dependent on underlying aquifers. This paper develops baseline understanding of trends in the ecology, water quality and uses of Bhojtal, discussing their implications for the long-term wellbeing of the Bhopal city region. It highlights increasing dependency on water diverted from out-of-catchment sources, and also abstraction across the Bhojtal catchment in excess of replenishment that is depressing groundwater and contributing to reported declining lake level and water quality. Despite some nature-based management initiatives, evidence suggests little progress in haltering on-going groundwater depression and declines in lake water level and quality. Significant declines in ecosystem services produced by Bhojtal are likely without intervention, a major concern given the high dependency of people in the Bhopal region on Bhojtal for their water supply and socio-economic and cultural wellbeing. Over-reliance on appropriation of water from increasingly remote sources is currently compensating for lack of attention to measures protecting or regenerating local resources that may provide greater resilience and regional self-sufficiency. Improved knowledge of catchment hydrogeology on a highly localised scale could improve the targeting and efficiency of water harvesting and other management interventions in the Bhojtal catchment, and their appropriate hybridisation with engineered solutions, protecting the catchment from unintended impacts of water extraction or increasing its carrying capacity, and also providing resilience to rising population and climate change. Ecosystem service assessment provides useful insights into the breadth of benefits of improved management of Bhojtal and its catchment.

Impacts of land-use changes on the groundwater recharge in the Ho Chi Minh city, Vietnam

Ho Chi Minh City (HCMC), Vietnam has undergone tremendous transformation in land-use practices in the past few decades. The groundwater-related issues have also been a major concern in the fast-growing southern city of Vietnam. Quantitative prediction of the impact on groundwater recharge due to changes in the land-use pattern of a watershed is crucial in developing sound groundwater management schemes. This study aims to evaluate the impacts of change in land-use patterns on the quantity of groundwater recharge in HCMC. An empirical land-use projection model (Conversion of Land-use and its Effects, Dyna-CLUE) and a hydrological model (Soil and Water Assessment Tool, SWAT) was used for the study. Three future land-use scenarios of Low Urbanization Scenario (LU), Medium Urbanization Scenario (MU) and High Urbanization Scenario (HU) were developed in Dyna-CLUE focusing on the increase of built-up area to generate land-use maps of HCMC until the year 2100. The land-use maps for all three scenarios were then used in the calibrated hydrological model SWAT to get the future recharge in the near future (2016-2045), mid future (2046-2075) and far future (2076-2100). The recharge was observed to increase in the far future of LU by 10% while reduction of 30% and 52% in annual average recharge was observed in far future of MU and HU respectively. It was, thus, observed that change in built-up area has a significant effect on the groundwater recharge in HCMC.

Soil moisture dynamics modelling of a reclaimed upland in the early post-construction period

Mine reclamation landscapes typically comprise layers of mine waste materials such as tailings sands, capped with a cover soil. In addition to the arrangement and placement of these materials, their hydraulic properties govern the performance of the built system. Soil evolution due to freeze-thaw cycling can result in dramatically altered soil hydraulic properties compared to the as-built material. Therefore, prediction of present and future hydrologic behaviour relies on understanding the nature and magnitude of this change and the elapsed time associated with stabilization. This research quantifies the transient hydraulic properties of mine reclamation materials at a constructed upland within a reclaimed watershed, and models the effect of this evolution on the partitioning of soil moisture between evaporation and groundwater recharge. Soil moisture dynamics were simulated using HYDRUS-1D for the ice-free period two, three, and five years after construction. A capillary barrier between the fine-grained cover soil and coarse-grained tailings sand regulated percolation past the interface. Soil evolution of the cover soil was responsible for an increase in saturated hydraulic conductivity by an order of magnitude, decrease in air-entry pressure by a factor of 4, and decrease in the van Genuchten n parameter by a factor of 2. The altered soil hydraulic properties associated with the weathered cover soil ultimately resulted in a 64% increase in groundwater recharge as a consequence of the capillary barrier weakening. The cover soil exhibited minor spatial heterogeneity in soil hydraulic properties, and did not contribute substantial uncertainty to the estimates of groundwater recharge and evaporation. Cover soil thickness exerted a strong influence on the partitioning of soil moisture. Reclaimed uplands will provide the most recharge to downgradient ecosystems in the period following the completion of soil evolution (~4 years) but preceding substantial vegetation development.

Rainfall, groundwater, and surface water isotope data from extreme tropical cyclones (2016-2019) within the Caribbean Sea and Atlantic Ocean basins

Under a changing climate, projections estimate that over the next thirty years, extreme Tropical Cyclones (TCs) will increase in frequency, with two to three times more Category 4 and 5 hurricanes in the Atlantic basin between 20°N and 40°N. In recent years, the Caribbean Sea and Atlantic Ocean basins have experienced several extreme TCs, resulting in extensive human, ecological, and economic damage [1], [2], [3]. To improve understanding of TCs and their potential impacts in the face of climate change, physically based understanding of past climate and modern TC dynamics is necessary. Despite the well-known Atlantic hurricane season, surface observations of the isotopic evolution of TC's moisture and the propagation of isotopically distinct pulses across surface and subsurface water reservoirs are lacking. In this data article, we provide novel high frequency sampling of surface rainfall isotope compositions (δ¹⁸O, δ²H, and d-excess in ‰) for Hurricanes Otto (Costa Rica, 2016), Nate (Costa Rica, 2017), Irma (Cuba and The Bahamas, 2017), Maria (Cuba and The Bahamas, 2017), and Dorian (The Bahamas, 2019). These five TCs were characterized by unprecedented impacts during continental and maritime landfalls and passages. In total, 161 surface rainfall samples were collected in passive devices [4] with event-based and daily frequencies, resulting in the first surface isotopic tempestology anatomy across the Caribbean Sea and Atlantic Ocean basins to date. Derived rainfall from TCs often results in large input amounts of isotopically distinct water over an area from few hours to several days, and therefore this unique isotope composition is propagated through surface and shallow subsurface reservoirs. Our data also include spring (N=338) and surface water (N=334) isotope compositions following the impact of Hurricane Otto and Tropical Storm Nate in central Costa Rica. As this region is well-known for its diverse rainfall dynamics and as a climate change 'hot spot' [5], [6], [7], our data provide an opportunity to improve and complement modern and past climate interpretations often derived from satellite products and calcite-δ¹⁸O paleoclimatic archives in light of climatic forcing, TC rainfall amounts and recharge rates, and the hypothesized climatic-induced decline of past Mesoamerican civilizations.

A global-scale dataset of direct natural groundwater recharge rates: A review of variables, processes and relationships

Groundwater recharge indicates the existence of renewable groundwater resources and is therefore an important component in sustainability studies. However, recharge is also one of the least understood, largely because it varies in space and time and is difficult to measure directly. For most studies, only a relatively small number of measurements is available, which hampers a comprehensive understanding of processes driving recharge and the validation of hydrogeological model formulations for small- and large-scale applications. We present a new global recharge dataset encompassing >5000 locations. In order to gain insights into recharge processes, we provide a systematic analysis between the dataset and other global-scale datasets, such as climatic or soil-related parameters. Precipitation rates and seasonality in temperature and precipitation were identified as the most important variables in predicting recharge. The high dependency of recharge on climate indicates its sensitivity to climate change. We also show that vegetation and soil structure have an explanatory power for recharge. Since these conditions can be highly variable, recharge estimates based only on climatic parameters may be misleading. The freely available dataset offers diverse possibilities to study recharge processes from a variety of perspectives. By noting the existing gaps in understanding, we hope to encourage the community to initiate new research into recharge processes and subsequently make recharge data available to improve recharge predictions.

Potential groundwater recharge from deep drainage of irrigation water

Knowledge of soil water dynamics in the deep vadose zone provides valuable information on the temporal and spatial variability of groundwater recharge. However, semi-arid climate can complicate how the input of water, such as irrigation, can contribute to potential groundwater recharge. This study assessed the recharge rates and their timing under irrigated cropland from a semi-arid region of northern Iran. A deep drainage (10 m) experiment was performed and in situ soil water content was measured to analyze the soil water dynamics and model hydraulic parameters using HYSDRUS-1D. The best parameters selected from inverse parameter optimization were used to calibrate model and estimate the long-term (20-year) average groundwater recharge and the influence of the root zone, unsaturated zone and the time scale on the recharge processes. The simulated annual flux ranged from 24 mm to 268 mm (mean of 110 mm) at 2-m depth and ranged between 26 mm to 207 mm (mean of 95 mm) at the 10-m depth. High fluxes, observed between December and April, may be the result of greater precipitation combined with the irrigation return flow. The May-October period showed a gradual decrease in flux at the depth of 2 m. At the depth of 10 m, the flux showed some continuity (base flux) during the long-term recharge simulation. In total, 12.7% of the input water contributed to the recharge of the groundwater. The annual soil water fluxes were almost similar irrespective of depth below the root zone and the flux rates did not show any clear relation between the different components of the water budget at any depth. This approach improved our understanding of the recharge process in the deep vadose zone in a semiarid region and can help for the development of effective management of groundwater resources.

Groundwater recharge sites and pollution sources in the wine-producing Guadalupe Valley (Mexico): Restrictions and mixing prior to transfer of reclaimed water from the US-México border

Rapid depletion of aquifers in semiarid and arid regions threatens water security. This holds true especially in emerging countries where insufficient knowledge about aquifer systems precludes the implementation of advanced management measures, such as managed aquifer recharge. This study deals with the generation of baseline knowledge for the assessment of aquifers in arid and semiarid regions where artificial recharge with reclaimed water gains increasing impetus. The Guadalupe aquifer in Baja California provides water to 57% of the Mexican wine industry. Recent plans foresee a partial replenishment of its depleted groundwater reserves by transferring treated waste water from the Mexico-USA border for irrigation. The aquifer demonstrated to have a rapid response by rising the water table of about +20 m in relation to natural recharge under an intense rainfall period of 236 mm. Two predominant recharge sources were identified based on a geochemical multi-tracer approach: (a) water of modern age (<5 yr, >1.8 TU) and mixed water of recent-submodern age (³H 0.8-1.8 TU), and (b) sub-modern waters that were recharged before 1952 (³H < 0.5 TU). Water of the first type originate in the main Guadalupe stream, which has a more depleted average δ¹⁸O isotope value (-7.8‰) than average local rainwater (-2.0‰). The stream water initially has a Na-HCO₃ composition and recharges the entire Calafia zone and most groundwater along the riverbed across the valley. Water of the second type is mostly derived from hill-slope groundwater that has a stable isotope composition of mixed local rainwater and a NaCl composition. High total dissolved solids >2 g l^-1 together with enriched NO₃^- and Se concentrations characterize groundwater in the downstream the Porvenir zone. The geochemical age of this older, hill-slope groundwater suggests that its replenishment takes at least several decades when it becomes exhausted.

Groundwater Management in Coastal Areas through Landscape Scale Planning: A Systematic Literature Review

Groundwater is one of the main resources for social-ecological systems. As part of the total water cycle and deeply connected with land use, groundwater management faces many challenges, especially in coastal areas. Landscape Scale Planning is an emerging approach for land use planning providing a framework for management based on evidence, given that landscapes have physical and information flows. Landscape Scale Planning embraces the following three dimensions: (i) the spatial dimension centres on the recognition of distinct landscape units; (ii) the temporal dimension entails past, current and future uses of a landscape; and (iii) the modification dimension involves the anthropogenic alterations that affected and will affect the landscape and its features along the spatial and temporal dimensions. Through a systematic literature review of 28 selected publications, this paper explores how groundwater management can be improved through a Landscape Scale Planning approach. The results show that Landscape Scale Planning can be applied as an integrative framework for groundwater management. Landscape units based on, but not limited to, geology, topography, cultural and socio-economic aspects can aid groundwater management to consider the differing spatial and temporal characteristics of the aquifer. Landscape Scale Planning can also favour the inclusion of land use change dynamics in groundwater management processes. To this end, the paper proposes guidelines for applying Landscape Scale Planning to inform groundwater management and consider land use changes.