Numerical simulation of a managed aquifer recharge system designed to supply drinking water to the city of Amsterdam,The Netherlands

Managed aquifer recharge (MAR) is increasingly used to secure drinking water supply worldwide. The city of Amsterdam (The Netherlands) depends largely on the MAR in coastal dunes for water supply. A new MAR scheme is proposed for the production of 10 × 10⁶ m³/year, as required in the next decade. The designed MAR system consists of 10 infiltration ponds in an artificially created sandbank, and 25 recovery wells placed beneath the ponds in a productive aquifer. Several criteria were met for the design, such as a minimum residence time of 60 days and maximum drawdown of 5 cm. Steady-state and transient flow models were calibrated. The flow model computed the infiltration capacity of the ponds and drawdowns caused by the MAR. A hypothetical tracer transport model was used to compute the travel times from the ponds to the wells and recovery efficiency of the wells. The results demonstrated that 98% of the infiltrated water was captured by the recovery wells which accounted for 65.3% of the total abstraction. Other sources include recharge from precipitation (6.7%), leakages from surface water (13.1%), and natural groundwater reserve (14.9%). Sensitivity analysis indicated that the pond conductance and hydraulic conductivity of the sand aquifer in between the ponds and wells are important for the infiltration capacity. The temperature simulation showed that the recovered water in the wells has a stable temperature of 9.8–12.5 °C which is beneficial for post-treatment processes. The numerical modelling approach is useful and helps to gain insights for implementation of the MAR.

     

Aquifer residence times for recycled water estimated using chemical tracers and the propagation of temperature signals at a managed aquifer recharge site in Australia

A prerequisite for minimizing contamination risk whilst conducting managed aquifer recharge (MAR) with recycled water is estimating the residence time in the zone where pathogen inactivation and biodegradation processes occur. MAR in Western Australia’s coastal aquifers is a potential major water source. As MAR with recycled water becomes increasingly considered in this region, better knowledge of applied and incidental tracer-based options from case studies is needed. Tracer data were collected at a MAR site in Floreat, Western Australia, under a controlled pumping regime over a distance of 50 m. Travel times for bromide-spiked groundwater were compared with two incidental tracers in recycled water: chloride and water temperature. The average travel time using bromide was 87?±?6 days, whereas the estimates were longer based on water temperature (102?±?17 days) and chloride (98?±?60 days). The estimate of average flow velocity based on water temperature data was identical to the estimate based on bromide within a 25-m section of the aquifer (0.57?±?0.04 m day^?1). This case study offers insights into the advantages, challenges and limitations of using incidental tracers in recycled water as a supplement to a controlled tracer test for estimating aquifer residence times.     

Tracing the decomposition of dissolved organic carbon in artificial groundwater recharge using carbon isotope ratios

Reducing the concentration of dissolved organic C (DOC) in water is one of the main challenges in the process of artificial groundwater recharge. At the Tuusula waterworks in southern Finland, surface water is artificially recharged into an esker by pond infiltration and an equal amount of groundwater is daily pumped from the aquifer. This groundwater study was conducted to consider the role of redox processes in the decomposition of DOC. The isotopic composition of dissolved inorganic C (δ¹³C_DIC) in the recharged water was used as a tracer for redox reactions. The isotopic composition of O and H in water was determined in order to calculate mixing ratios between the local groundwater and the infiltrated surface water. Three distinct processes in the reduction of the DOC content were traced using isotopic methods and concentration analyses of DIC and DOC: (1) the decomposition of DOC, (2) adsorption of DOC on mineral matter, and (3) the dilution of artificially recharged water by mixing with local groundwater. The largest decrease (44%) in the DOC content occurred during the early stage of subsurface flow, within 350 m of the infiltration ponds. The reduction of DOC was accompanied by an equal increase in DIC and a significant drop in δ¹³C_DIC. This change is attributed to the oxidative decomposition of DOC. A further 23% decrease in DOC is attributed to adsorption and a final drop of 14% to dilution with local groundwater.     

Hydrogeochemical and isotopic evidence of groundwater evolution and recharge in aquifers in Beijing Plain, China

An investigation was conducted in Beijing to identify the groundwater evolution and recharge in the quaternary aquifers. Water samples were collected from precipitation, rivers, wells, and springs for hydrochemical and isotopic measurements. The recharge and the origin of groundwater and its residence time were further studied. The groundwater in the upper aquifer is characterized by Ca-Mg-HCO₃ type in the upstream area and Na-HCO₃ type in the downstream area of the groundwater flow field. The groundwater in the lower aquifer is mainly characterized by Ca-Mg-HCO₃ type in the upstream area and Ca-Na-Mg-HCO₃ and Na-Ca-Mg-HCO₃ type in the downstream area. The δD and δ¹⁸O in precipitation are linearly correlated, which is similar to WMWL. The δD and δ¹⁸O values of river, well and spring water are within the same ranges as those found in the alluvial fan zone, and lay slightly above or below LMWL. The δD and δ¹⁸O values have a decreasing trend generally following the precipitation → surface water → shallow groundwater → spring water → deep groundwater direction. There is evidence of enrichment of heavy isotopes in groundwater due to evaporation. Tritium values of unconfined groundwater give evidence for ongoing recharge in modern times with mean residence times <50 a. It shows a clear renewal evolution along the groundwater flow paths and represents modern recharge locally from precipitation and surface water to the shallow aquifers (<150 m). In contrast, according to ¹⁴C ages in the confined aquifers and residence time of groundwater flow lines, the deep groundwater is approximately or older than 10 ka, and was recharged during a period when the climate was wetter and colder mainly from the piedmont surrounding the plain. The groundwater exploitation is considered to be “mined unsustainably” because more water is withdrawn than it is replenished.     

屋顶雨水回灌裂隙岩溶含水层连通示踪试验

雨水回灌后对岩溶水的影响分析和效果需要对回灌井和下游的观测井监测数据说明。由于岩溶含水层中岩溶发育极不均匀,在测定地下水位并判定出地下水流向的基础上,连通示踪试验是选择与回灌井相对应的观测井的有效办法。用氯化钠作为示踪剂,监测不同测井不同埋深地下水电导率随时间变化的方法确定了岩溶水观测井选址、岩溶含水层的结构及水文地质条件。示踪试验结果表明,3#观测井(西院井)可用来监测雨水回灌效果,并计算出该地区的地下水势流速在9 m/h～20 m/h之间。该成果对北方地区岩溶含水层回灌和环境评价具有重要的参考价值。     

Hydrogeochemical transport modeling of the infiltration of tertiary treated wastewater in a dune area, Belgium

Managed artificial recharge (MAR) is a well-established practice for augmentation of depleted groundwater resources or for environmental benefit. At the St-André MAR site in the Belgian dune area, groundwater resources are optimised through re-use of highly treated wastewater by means of infiltration ponds. The very high quality of the infiltration water sets this system apart from other MAR systems. The low total dissolved solid (TDS) content in the infiltration water (less than 50 mg/L) compared to the dune aquifer (500 mg/L) triggers a number of reactions, increasing the TDS through soil-aquifer passage. Multi-component reactive transport modelling was applied to analyse the geochemical processes that occur. Carbonate dissolution is the main process increasing the TDS of the infiltration water. Oxic aquifer conditions prevail between the infiltration ponds and the extraction wells. This is driven by the high flow velocities, leaving no time to consume O₂ between the ponds and extraction wells. Cation exchange is important when infiltration water is replaced by native dune water or when significant changes in infiltration-water quality occur. The seasonal variation of O₂ and temperature in the infiltration water are the main drivers for seasonal changes in the concentration of all major ions.     

Isotopic characteristics of precipitation, surface and ground waters in the Yinchuan plain, Northwest China

The isotopic characteristics of hydrogen and oxygen of various water bodies in the Yinchuan plain were investigated. A total of 131 water samples were collected and another 99 water samples collected by other scholars and the International Atomic Energy Agency were referred to in the study. The stable isotopic compositions of precipitation and lake waters are influenced by dry climate and strong evaporation of the area, and the slope of local meteoric water line and lake water line are smaller than that of global meteoric water line. The isotopic compositions of the Yellow River water are significantly lower than the weighted averages of local atmospheric precipitation and are controlled by the runoff of upstream areas. The isotopes of phreatic water suggest that the single phreatic water is predominantly recharged by the bedrock fissure water of Helan Mountain, while the phreatic water in multilayer structure area is recharged by multiple sources. The confined waters in the area may be recharged under cooler climate conditions than the present, which makes the heavy isotopes depleted. The ³H contents of various water bodies are in the following order: the Yellow River water > multilayer phreatic water > single phreatic water > lake water > drainage ditch water > upper confined water > lower confined water. This is in accordance with the recharge, discharge and circulation conditions of the specific water bodies. The isotope compositions of waters in the Yinchuan plain are mainly affected by external water recharge rather than local precipitation. This study is meaningful and helpful in understanding the groundwater flow systems and water cycle in the area.     

Groundwater residence times and flow paths in fractured rock determined using environmental tracers in the Mission Tunnel; Santa Barbara County,California, USA

Multiple geochemical tracers [ion chemistry, stable isotopes of water, chlorofluorocarbons (CFC), tritium] and a 25-year-long record of discharge were used to understand residence times and flow paths of groundwater seeps in the fractured rock aquifer surrounding the Mission Tunnel, Santa Barbara, California. Tritium data from individual seeps indicate that seep waters are a mixture of >45-year-old (recharged prior to the nuclear bomb tests) and young groundwater. CFC data support this interpretation, however, a two-end member mixing model cannot completely explain the age tracer data. Microbial degradation and partial re-equilibration complicate the CFC signal. Spectral analysis of precipitation and groundwater seepage records shows that seepage lags precipitation by 3 months. This delay is related to the advancement of the wetting front and increasing the number of active flow paths. Additionally, the amount of seepage produced by precipitation is less during extended periods of drought than during normal or wet periods, suggesting antecedent conditions strongly affect flow through this fractured rock aquifer.     

Effects of rainwater-harvesting-induced artificial recharge on the groundwater of wells in Rajasthan, India

In light of the increasing deterioration of groundwater supplies in Rajasthan, India, rainwater harvesting practices in southern Rajasthan were studied to determine the effects of artificially recharged groundwater on the supply and quality of local groundwater. A physical and geochemical investigation utilizing environmental tracers (δ¹⁸O and Cl^–), groundwater level and groundwater quality measurements, and geological surveys was conducted with two objectives: (1) to quantify the proportion of artificially recharged groundwater in wells located near rainwater harvesting structures and (2) to examine potential effects of artificial recharge on the quality of groundwater in these wells. A geochemical mixing model revealed that the proportion of artificial recharge in these wells ranged from 0 to 75%. Groundwater tracer, water table, and geological data provided evidence of complex groundwater flow and were used to explain the spatial distribution of artificial recharge. Furthermore, wells receiving artificial recharge had improved groundwater quality. Statistical analysis revealed a significant difference between the water quality in these wells and wells determined not to receive artificial recharge, for electrical conductivity and SO ₄ ^– . The findings from this study provide quantitative evidence that rainwater harvesting structures in southern Rajasthan influence the groundwater supply and quality of nearby wells by artificially recharging local groundwater.     

Patterns in groundwater chemistry resulting from groundwater flow

 Groundwater flow influences hydrochemical patterns because flow reduces mixing by diffusion, carries the chemical imprints of biological and anthropogenic changes in the recharge area, and leaches the aquifer system. Global patterns are mainly dictated by differences in the flux of meteoric water passing through the subsoil. Within individual hydrosomes (water bodies with a specific origin), the following prograde evolution lines (facies sequence) normally develop in the direction of groundwater flow: from strong to no fluctuations in water quality, from polluted to unpolluted, from acidic to basic, from oxic to anoxic–methanogenic, from no to significant base exchange, and from fresh to brackish. This is demonstrated for fresh coastal-dune groundwater in the Netherlands. In this hydrosome, the leaching of calcium carbonate as much as 15 m and of adsorbed marine cations (Na⁺, K⁺, and Mg²⁺) as much as 2500 m in the flow direction is shown to correspond with about 5000 yr of flushing since the beach barrier with dunes developed. Recharge focus areas in the dunes are evidenced by groundwater displaying a lower prograde quality evolution than the surrounding dune groundwater. Artificially recharged Rhine River water in the dunes provides distinct hydrochemical patterns, which display groundwater flow, mixing, and groundwater ages. Received, May 1998 · Revised, August 1998 · Accepted, October 1998     

Thermal tracer testing in a sedimentary aquifer: field experiment (Lauswiesen,Germany) and numerical simulation

An active and short-duration thermal tracer test (TTT) was conducted in a shallow sedimentary aquifer at the Lauswiesen test site, near Tübingen, Germany. By injecting 16  m³ of warm water at 22°C, a thermal anomaly was created, which propagated along the local groundwater flow direction. This was comprehensively monitored in five observation wells at a few meters distance. The purpose of this well-controlled experiment was to determine the practicability of such a TTT and its suitability to examine hydraulic characteristics of heterogeneous aquifers. The results showed that the thermal peak arrival times in the observation wells were consistent with previous observations from alternative field testing such as direct-push injection logging (DPIL). Combined analysis of depth-dependent temperatures and peak arrival times, and comparison with a numerical heat transport model, offers valuable insights into the natural flow field and spatial distribution of hydraulic conductivities. The study was able to identify vertical flow focusing and bypassing, which are attributed to preferential flow paths common in such sedimentary sand and gravel aquifers. These findings are fundamental for further development of experimental designs of active and short-duration TTTs and provide a basis for a more quantitative analysis of advective and conductive transport processes.     

Characterization of karst conduits by tracer tests for an artificial recharge scheme

Challenged by rapidly changing climate in combination with an increase in anthropogenic pressures, karst groundwater resources in the Old Town of Lijiang (OTLJ), SW China, are diminishing. Higher frequency and longer duration of dried-up periods have been observed at the Heilongtan Park (HP) Springs in recent years. Thus, there is an urgent need for an artificial recharge scheme, aimed at replenishing groundwater in the aquifer and increasing the outflow of the springs to ensure effective water resources management. Evaluation of the scheme feasibility, prior to its implementation, is important. In this study, tracer tests were conducted between the recharge area and receiving springs in order to gain insight into the transport mechanisms of karst groundwater and the structural characteristics of the aquifer. Multiple underground flow paths, exhibiting high conductivity between the recharge area and HP Springs, were revealed by the interpretation of tracer breakthrough curves. Three springs considered as the leakages of the scheme were identified. Moreover, the outflow of springs at HP and OTLJ were predicted to be increased by the artificially recharged water after 9.2 and 12.5 days, respectively. Quantitative analysis of tracer recoveries demonstrates that the springs to be recharged and the springs considered as leakages, respectively, share 45 and 55% of the increased outflow. The feasibility of the scheme has been confirmed by the tracer tests. This report provides references for the evaluation of artificial groundwater recharge and protection strategies, particularly in large and poorly investigated karst spring fields.

     

Use of artificial and natural tracers to assess groundwater transit-time distribution and flow systems in a high-alpine karst system (Wetterstein Mountains,Germany)

Groundwater in mountainous karst regions is vital for regional water budgets and freshwater supply. Owing to increasing water demand and climate change, detailed knowledge of the highly heterogeneous alpine aquifer systems is required. Multi-tracer analyses have been conducted in the steep karstic Wetterstein Mountains, which includes Germany’s highest summit, Zugspitze (2,962 m asl). Results of artificial tracer tests demonstrate well-developed flow paths through the unsaturated zone (up to 1,000 m thickness). Flow paths cross topographic divides and contribute to deep drainage systems underneath alpine valleys. Cross-formational flow has been identified. Quantitative analysis of tailing-dominated breakthrough curves and stable isotopes (¹⁸O) has enabled determination of the mean transit-time distribution. A fast-flow component with transit times between 3 and 13 days was found in karst conduits and open fissures, dependent on flow conditions. An intermediate-flow component, showing mean transit times of about 2.9–4.9 months, was found in well-drained fissures and fractures. A slow-flow component with mean transit times greater than 1 year is attributable to slow flow and low storage in the poorly drained fissures and rock matrix. The conceptual model enables a better understanding of drainage, water resources and vulnerability of the high-alpine karst system.     

Dynamics and anthropogenic impacts of multiple karst flow systems in a mountainous area of South China

The Xiangxi River basin, South China, is a steep terrane with well-developed karst features and an important Cambrian-Ordovician aquifer. Meteoric water in this mountainous area features a mean δ¹⁸O elevation gradient of –2.4?‰/km. This gradient was used to estimate mean recharge elevations of 760 m for Shuimoxi (SMX) spring, 1,060 m for Xiangshuidong (XSD) spring, and 1,430 m for drill hole ZK03, indicating multiple flow paths in the Cambrian-Ordovician karst aquifer. Mean residence times of 230 and 320 days and ～2 years were estimated for these features, respectively, using the damped running average model that predicts the isotopic variations in groundwater from those in precipitation. Groundwater in the regional karst flow system has the longest residence time, the highest recharge elevation, the longest flow paths, the lowest addition of anthropogenic components, and the greatest amount of water–rock interaction as indicated by its higher dissolved solids, Mg²⁺ concentrations and Mg/Ca ratios than the springs. In contrast, the local and shallow karst flow systems respond rapidly to recharge events. Artificial tracer tests prove that these shallow karst systems can also quickly transmit anthropogenic contaminants, indicating that they are highly vulnerable to human impacts, which include the enrichment of NO₃ ^–. The intensity of water–rock interaction and groundwater vulnerability are mainly determined by the structure and dynamics of the multiple karst flow systems.     

Residence times of shallow groundwater in West Africa: implications for hydrogeology and resilience to future changes in climate

Although shallow groundwater (<50 mbgl) sustains the vast majority of improved drinking-water supplies in rural Africa, there is little information on how resilient this resource may be to future changes in climate. This study presents results of a groundwater survey using stable isotopes, CFCs, SF₆, and ³H across different climatic zones (annual rainfall 400–2,000 mm/year) in West Africa. The purpose was to quantify the residence times of shallow groundwaters in sedimentary and basement aquifers, and investigate the relationship between groundwater resources and climate. Stable-isotope results indicate that most shallow groundwaters are recharged rapidly following rainfall, showing little evidence of evaporation prior to recharge. Chloride mass-balance results indicate that within the arid areas (<400 mm annual rainfall) there is recharge of up to 20 mm/year. Age tracers show that most groundwaters have mean residence times (MRTs) of 32–65 years, with comparable MRTs in the different climate zones. Similar MRTs measured in both the sedimentary and basement aquifers suggest similar hydraulic diffusivity and significant groundwater storage within the shallow basement. This suggests there is considerable resilience to short-term inter-annual variation in rainfall and recharge, and rural groundwater resources are likely to sustain diffuse, low volume abstraction.     

Aquifer system response to intensive pumping in urban areas of the Gangetic plains,India: the case study of Patna

A significant component of domestic demand for water of urban areas located in the Gangetic plains is met by heavy pumping of groundwater. The present study is focused on the Patna municipal area, inhabited by 17 million people and spanning over 134 km², where entire urban water demand is catered from pumping by wells of various capacities and designs. The present study examines the nature of the aquifer system within the urban area, the temporal changes in the water/piezometric level and the recharge mechanism of the deeper aquifers. The aquifer system is made up of medium-to-coarse unconsolidated sand, lying under a ~40-m-thick predominantly argillaceous unit holding 8- to 13-m-thick localised sand layers and continues up to 220 m below ground. Groundwater occurs under semi-confined condition, with transmissivity of aquifers in 5,500–9,200 m² day^?1 range. Hydraulic head of the deeper aquifer remains in 9–19 m range below ground, in contrast to 1–9 m range of that of the upper aquitard zone. The estimated annual groundwater extraction from the deeper aquifer is ~212.0 million m³, which has created a decline of 3.9 m in the piezometric level of the deeper aquifer during the past 30 years. Unregulated construction of deep tube wells with mushrooming of apartment culture may further exacerbate the problem. The sand layers within the aquitard zone are experiencing an annual extraction of 14.5 million m³ and have exhibited stable water level trend for past one and half decades. This unit is recharged from monsoon rainfall, besides contribution from water supply pipe line leakage and seepage from unlined storm water drains.     

Assessment of key transport parameters in a karst system under different dynamic conditions based on tracer experiments: the Jeita karst system,Lebanon

Artificial tracer experiments were conducted in the mature karst system of Jeita (Lebanon) under various flow conditions using surface and subsurface tracer injection points, to determine the variation of transport parameters (attenuation of peak concentration, velocity, transit times, dispersivity, and proportion of immobile and mobile regions) along fast and slow flow pathways. Tracer breakthrough curves (TBCs) observed at the karst spring were interpreted using a two-region nonequilibrium approach (2RNEM) to account for the skewness in the TBCs’ long tailings. The conduit test results revealed a discharge threshold in the system dynamics, beyond which the transport parameters vary significantly. The polynomial relationship between transport velocity and discharge can be related to the variation of the conduit’s cross-sectional area. Longitudinal dispersivity in the conduit system is not a constant value (α?=?7–10 m) and decreases linearly with increasing flow rate because of dilution effects. Additionally, the proportion of immobile regions (arising from conduit irregularities) increases with decreasing water level in the conduit system. From tracer tests with injection at the surface, longitudinal dispersivity values are found to be large (8–27 m). The tailing observed in some TBCs is generated in the unsaturated zone before the tracer actually arrives at the major subsurface conduit draining the system. This work allows the estimation and prediction of the key transport parameters in karst aquifers. It shows that these parameters vary with time and flow dynamics, and they reflect the geometry of the flow pathway and the origin of infiltrating (potentially contaminated) recharge.     

Geochemical characterization and origins of the thermal springs in southern Gaoligong Mountains,China

Hydrogeochemistry and environmental isotope data were utilized to understand origins and characteristics of the thermal springs in southern Gaoligong Mountains, China. The groundwater at the thermal springs has low values of total dissolved solids, and its main water types are Na-HCO₃. The thermal springs are mainly recharged from meteoric precipitations. The recharge areas are located near the springs at an approximate elevation of 1,800 m. The groundwater of the thermal springs is immature and partially equilibrated with a strong mixture of the shallow cold waters during the flow process. The shallow cold water accounts for more than 90 %. The temperatures of thermal reservoir that feed the springs are between 146 and 260 °C, and the calculated groundwater circulation depths range from 2,000 to 5,700 m below ground surface.     

Distribution and migration mechanism of fluoride in groundwater in the Manas River Basin,Northwest China

Elevated fluoride (F) concentration in groundwater is posing a public health risk in the Manas River Basin (MRB), Northwest China. Based on the characterization of regional groundwater flow, 90 groundwater samples from aquifers were analyzed, along with top-soil leachate and pore-water samples from aquitards. Stable oxygen (δ¹⁸O) and hydrogen isotopes, radiocarbon and hydrochemical analyses of the groundwater and pore-water samples were conducted to trace groundwater hydrological and hydrochemical processes and thereby understand the distribution and migration mechanism of F. The groundwater is recharged by meteoric precipitation through vapor condensation processes in the Tianshan Mountains. The F concentration in groundwater samples from this basin ranged from 0.11 to 48.15 mg/L (mean 2.56 mg/L). In 37 of the 90 groundwater samples, the F concentrations were above the safe level for drinking water. The F concentrations progressively increased with the residence time and well depths in the northwest of the alluvial-fluvial plain, where groundwater is overexploited for agricultural and domestic use. Positive correlations between F and sodium (Na)/calcium (Ca) indicate that the enrichment and migration of F are influenced by cation exchange processes under high-Na and alkaline pH conditions. The relationships between δ¹⁸O and F and chloride (Cl) concentrations were nonlinear due to leaching and mixing processes. This shows that vertical leaching by irrigation return flow and mixing with pore water are the dominant processes driving the migration of F in the groundwater flow system of MRB, in addition to geochemical processes.