Stable isotopic and geochemical identification of groundwater evolution and recharge sources in the arid Shule River Basin of Northwestern China

Stable isotopic (δD_VSMOW and δ¹⁸O_VSMOW) and geochemical signatures were employed to constrain the geochemical evolution and sources of groundwater recharge in the arid Shule River Basin, Northwestern China, where extensive groundwater extraction occurs for agricultural and domestic supply. Springs in the mountain front of the Qilian Mountains, the Yumen‐Tashi groundwater (YTG), and the Guazhou groundwater (GZG) were Ca‐HCO₃, Ca‐Mg‐HCO₃‐SO₄ and Na‐Mg‐SO₄‐Cl type waters, respectively. Total dissolved solids (TDS) and major ion (Mg²⁺, Na⁺, Ca²⁺, K⁺, SO₄^2?, Cl^? and NO₃^?) concentrations of groundwater gradually increase from the mountain front to the lower reaches of the Guazhou Basin. Geochemical evolution in groundwater was possibly due to a combination of mineral dissolution, mixing processes and evapotranspiration along groundwater flow paths. The isotopic and geochemical variations in melt water, springs, river water, YTG and GZG, together with the end‐member mixing analysis (EMMA) indicate that the springs in the mountain front mainly originate from precipitation, the infiltration of melt water and river in the upper reaches; the lateral groundwater from the mountain front and river water in the middle reaches are probably effective recharge sources for the YTG, while contribution of precipitation to YTG is extremely limited; the GZG is mainly recharged by lateral groundwater flow from the Yumen‐Tashi Basin and irrigation return flow. The general characteristics of groundwater in the Shule River Basin have been initially identified, and the results should facilitate integrated management of groundwater and surface water resources in the study area. Copyright © 2015 John Wiley & Sons, Ltd.     

A two‐dimensional analytical solution for groundwater flow in a leaky confined aquifer system near open tidal water

Groundwater in coastal areas is commonly disturbed by tidal fluctuations. A two‐dimensional analytical solution is derived to describe the groundwater fluctuation in a leaky confined aquifer system near open tidal water under the assumption that the groundwater head in the confined aquifer fluctuates in response to sea tide whereas that of the overlying unconfined aquifer remains constant. The analytical solution presented here is an extension of the solution by Sun for two‐dimensional groundwater flow in a confined aquifer and the solution by Jiao and Tang for one‐dimensional groundwater flow in a leaky confined aquifer. The analytical solution is compared with a two‐dimensional finite difference solution. On the basis of the analytical solution, the groundwater head distribution in a leaky confined aquifer in response to tidal boundaries is examined and the influence of leakage on groundwater fluctuation is discussed. Copyright © 2001 John Wiley & Sons, Ltd.     

Using hydrochemical,stable isotope,and river water recharge data to identify groundwater flow paths in a deeply buried karst system

The deeply buried river‐connected Xishan karst aquifer (XKA) in western Beijing, China, has been suffering from diminishing recharge for several decades, which in turn leads to the disappearing of spring water outflows and continuously lowering of groundwater level in the area. Thus, it is important to correctly recognize the groundwater recharge and flow paths for the sustainable development of the XKA. To investigate these issues, the hydrochemical and isotopic compositions are analysed for both surface water and groundwater samples collected over an area of about 280 km². Results show that (a) the river water is characterized by high Na contents; (b) the δ²H and δ¹⁸O values in the river water are distinctively higher than those of groundwater samples, after experiencing the long‐time evaporative enrichment in the upstream reservoir; (c) the Sr concentrations and ⁸⁷Sr/⁸⁶Sr ratios of groundwater clearly indicated the interaction between water and carbonate minerals but excluded the water–silicate interaction; and (d) the groundwater samples in the direct recharge area of the XKA have the lowest Na concentrations and the δ²H and δ¹⁸O values. Based on the large differences in the Na contents and ¹⁸O values of groundwater and surface water, a simple two‐component mixing model is developed for the study area and the fractions of the river water are estimated for groundwater samples. We find that the distribution pattern of the river water fractions in the XKA clearly shows a change of directions in the preferential flow path of the groundwater from its source zone to the discharge area. Overall, our results suggest that the recharged surface water can be a useful evidence for delineating the groundwater flow path in river‐connected karst aquifer. This study improves our understanding of the heterogeneity in karst groundwater systems.     

Groundwater mixing between different aquifer types in a complex structural setting discerned by elemental and stable isotope geochemistry

Karst aquifers are well known for their intricate stratigraphy and geologic structures, which make groundwater characterization challenging because flowpaths and recharge sources are complex and difficult to evaluate. Geochemical data, collected from ten closely spaced production wells constructed in two karst aquifers (Bangor Limestone (Mb) and Tuscumbia Limestone/Fort Payne Chert (Mftp)) in Trussville, north‐central Alabama, illustrate two distinctive groundwater end‐members: (1) higher major ion, dissolved inorganic carbon, conductivity, alkalinity concentrations, heavier δ¹³C ratios (max: −10.2 ± 0.2‰ Vienna Pee Dee Belemnite (PDB)) and lower residence times (mean: 19.5 ± 2 years, n = 2) of groundwater in the Mb aquifer and (2) lower constituent concentrations, lighter δ¹³C ratios (min: −13.4 ± 0.2‰ PDB) and longer residence times of groundwater (mean: 23.6 ± 2 years, n = 4) in the Mftp aquifer. Summer and fall data and the binary mixing model show aquifer inter‐flow mixing along solution fractures and confirms the distinctive groundwater geochemistry of the two aquifers. Lowering of static water levels over the summer (drawdown from 2 to 5.2 m) leads to more reducing groundwater conditions (lower Eh values) and slightly enriched δ¹⁸O and δD ratios during the fall [δ¹⁸O: −4.8 ± 0.1 to −5.4 ± 0.1‰ Vienna Standard Mean Oceanic Water (VSMOW), n = 9; δD: −25.4 ± 1 to −27.4 ± 1‰ VSMOW, n = 9] when compared with summer season samples (δ¹⁸O: −5.1 ± 0.1 to −5.7 ± 0.1‰ VSMOW, n = 11; δD: −25.0 ± 1 to −30.6 ± 1‰ VSMOW, n = 11). GIS analyses confirm the localized origin of recharge to the investigated aquifers. The combination of GIS, field parameters and geochemistry analyses can be successfully used to identify recharge sources, evaluate groundwater flow and transport pathways and to improve understanding of how groundwater withdrawals impact the sustainability and susceptibility to contamination of karst aquifers. Copyright © 2015 John Wiley & Sons, Ltd.     

Isotopes and geochemistry in a managed aquifer recharge scheme: a case study of fresh water injection at the Damascus University Campus,Syria

A combination of stable isotopes (¹⁸O and ²H) and hydrochemistry has been applied to investigate storage processes in relation to aquifer storage and recovery (ASR) of the shallow alluvial Quaternary aquifer in Damascus basin. The stored water, entirely taken from the Figeh springs during flood periods, was injected in a single well having a brackish groundwater. Water samples were collected from four observation wells drilled in the Damascus University Campus (DUC) site during a 3‐year period (2006–2008). The injectant water, which deviates in its chemical and isotopic signatures from that of the ambient groundwater, shows that the stored water plume remains within close proximity to the injection well (IW) (<≈ 100 m). Thus, only two wells (W13 and W14) located at a distance less than 80 m from the injection point were affected by this injection. The observation wells located at longer distances from the IW (≈145 m and ≈ 600 m for wells W15 and WHz, respectively) were completely unaffected by the injection. Although most of the chemical and isotopic parameters usefully reflected the mixing process that occurs between the injectant water and ambient groundwater, the stable isotope (¹⁸O) and chloride (Cl⁻) were the most sensitive parameters that quickly reflect this signature. Using a simple mass balance, the calculated proportion of injectant water reaching the well W13 was in the range of 50–90%. This proportion was even lower (30–55%) in the case of well W14. Although the drought event prevailing during this study did not much help to inject further amounts of water, higher than the injected volume (0·2416 M m³) and also not favourable to better evaluate the fate and subsurface hydrological processes, these findings offer encouragement to continue the ASR activities, as an alternative way for better management of water resources in this basin facing intensive problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Identification of groundwater recharge sources and processes in a heterogeneous alluvial aquifer: results from multi‐level monitoring of hydrochemistry and environmental isotopes in a riverside agricultural area in Korea

We evaluated sources and pathways of groundwater recharge for a heterogeneous alluvial aquifer beneath an agricultural field, based on multi‐level monitoring of hydrochemistry and environmental isotopes of a riverside groundwater system at Buyeo, Korea. Two distinct groundwater zones were identified with depth: (1) a shallow oxic groundwater zone, characterized by elevated concentrations of NO₃^? and (2) a deeper (>10–14 m from the ground surface) sub‐oxic groundwater zone with high concentrations of dissolved Fe, silica, and HCO₃^?, but little nitrate. The change of redox zones occurred at a depth where the aquifer sediments change from an upper sandy stratum to a silty stratum with mud caps. The δ¹⁸O and δ²H values of groundwater were also different between the two zones. Hydrochemical and δ¹⁸O? δ²H data of oxic groundwater are similar to those of soil water. This illustrates that recharge of oxic groundwater mainly occurs through direct infiltration of rain and irrigation water in the sandy soil area where vegetable cropping with abundant fertilizer use is predominant. Oxic groundwater is therefore severely contaminated by agrochemical pollutants such as nitrate. In contrast, deeper sub‐oxic groundwater contains only small amounts of dissolved oxygen (DO) and NO₃^?. The ³H contents and elevated silica concentrations in sub‐oxic groundwater indicate a somewhat longer mean residence time of groundwater within this part of the aquifer. Sub‐oxic groundwater was also characterized by higher δ¹⁸O and δ²H values and lower d‐excess values, indicating significant evaporation during recharge. We suggest that recharge of sub‐oxic groundwater occurs in the areas of paddy rice fields where standing irrigation and rain water are affected by strong evaporation, and that reducing conditions develop during subsequent sub‐surface infiltration. This study illustrates the existence of two groundwater bodies with different recharge processes within an alluvial aquifer. Copyright © 2009 John Wiley & Sons, Ltd.     

Characterizing the climatic water balance dynamics and different runoff components in a poorly gauged tropical forested catchment,Nicaragua

ABSTRACT

The water balance dynamics and runoff components of a tropical forested catchment (46?km²) on the southwestern Pacific coast of Nicaragua were studied combining hydrometry, geological characterization and hydrochemical and isotopic tracers (three-component hydrograph separation). The climatic water balance was estimated for 2010/11, 2011/12 and 2012/13 with net values of 811?mm year^-1, 782?mm year^-1 and –447?mm year^-1, respectively. Runoff components were studied at different spatial and temporal scales, demonstrating that different sources and temporal contributions are controlled by dominant landscape elements and antecedent rainfall. In forested sub-catchments, permeable soils, stratigraphy and steep slopes favour subsurface stormflow generation contributing 50% and 53% to total discharge. At catchment scale, landscape elements such as smooth slopes, wide valleys, deeper soils and water table allow groundwater recharge during rainfall events. Groundwater dominates the hydrograph (50% of total discharge) under dry prior conditions. However, low soil infiltration capacity generates a larger surface runoff component (42%) under wet prior conditions which dominates total discharge. Our results show that forested areas are important to reduce surface runoff and thus soil degradation, which is relevant for the design of water management plans.

Editor D. Koutsoyiannis Associate editor D. Gerten     

Significance of isotopic and geochemical methods to determine the evolution of inland brackish and bitter water: An example from the Zuli river in the upper reaches of the Yellow River,China

With the increasing demand for water resources, the utilization of marginal water resources of poor-quality has become a focus of attention. The brackish water developed in the Loess Plateau is not only salty but also famous for its ‘bitterness’. In the present work, multi-isotope analysis (Sr, B) was combined with geochemical analysis to gain insight into the hydrogeochemical evolution and formation mechanisms of brackish water. These results demonstrate that groundwater in the headwater is influenced by carbonate weathering. After the confluence of several tributaries in the headwater, the total dissolved solids (TDS) of water is significantly increased. The dissolution of evaporates is shown to be the main source of salinity in brackish water, which also greatly affects the strontium isotopic composition of water. This includes the dissolution of Mg-rich minerals, which is the main cause of the bitterness. Furthermore, the release of calcium from the dissolution of gypsum may induce calcite precipitation and incongruent dissolution of dolomite, which also contributes to the enrichment of magnesium. The highly fractionated boron isotopic values observed in the upstream groundwater were explained by the absorption with clay minerals. The inflow of brackish groundwater is the source of river water. Then evaporation further aggravates the salinization of river water, with water quality evolving to saline conditions in the lower reach. When the river reaches the valley plain, the ⁸⁷Sr/⁸⁶Sr ratios decreases significantly, which is primarily related to erosion of the riverbanks during runoff. These results indicate that water resource sustainability could be enhanced by directing focus to mitigating salinization in the source area of the catchment.     

Linking tracers,water age and conceptual models to identify dominant runoff processes in a sparsely monitored humid tropical catchment

Assessing catchment runoff response remains a key research frontier because of limitations in current observational techniques to fully characterize water source areas and transit times in diverse geographical environments. Here, we report a study that combines empirical data with modelling to identify dominant runoff processes in a sparsely monitored humid tropical catchment. The analysis integrated isotope tracers into conceptual rainfall–runoff models of varying complexity (from 5 to 11 calibrated parameters) that are able to simulate discharge and tracer concentrations and track the evolving age of stream water exiting the catchment. The model structures can be seen as competing hypotheses of catchment functioning and were simultaneously calibrated against uncertain streamflow gaugings and a 2‐year daily isotope rainfall–runoff record. Comparison of the models was facilitated using global parameter sensitivity analysis and the resulting effect on calibration. We show that a variety of tested model structures reproduced water and tracer dynamics in stream, but the simpler models failed to adequately reproduce both. The resulting water age distributions of the tested models varied significantly with little similarity between the stream water age and stored water age distributions. The sensitivity analysis revealed that only some of the more complex models (from eight parameters) could be better constrained to infer more plausible water age distributions and catchment storage estimates. These models indicated that the age of water stored in the catchment is generally older compared with the age of water fluxes, with evapotranspiration age being younger compared with streamflow. However, the water age distributions followed a similar temporal behaviour dominated by climatic seasonality. Stream water ages increased during the dry season (greater than 1 year) and decreased with increased streamflow (a few weeks old) during the wet season. We further show that the ratios of the streamwater age to stored water age distribution and the water age distribution of actual evapotranspiration to the stored water age distribution from constrained models could potentially serve as useful hydrological indicators of catchment functioning. Copyright © 2016 John Wiley & Sons, Ltd.     

Groundwater recharge in the Akaki catchment,central Ethiopia: evidence from environmental isotopes (δ18O, δ2H and 3H) and chloride mass balance

Recharge patterns, possible flow paths and the relative age of groundwater in the Akaki catchment in central Ethiopia have been investigated using stable environmental isotopes δ¹⁸O and δ²H and radioactive tritium (³H) coupled with conservative chloride measurements. Stable isotopic signatures are encoded in the groundwater solely from summer rainfall. Thus, groundwater recharge occurs predominantly in the summer months from late June to early September during the major Ethiopian rainy season. Winter recharge is lost through high evaporation–evapotranspiration within the unsaturated zone after relatively long dry periods of high accumulated soil moisture deficits. Chloride mass balance coupled with the isotope results demonstrates the presence of both preferential and piston flow groundwater recharge mechanisms. The stable and radioactive isotope measurements further revealed that groundwater in the Akaki catchment is found to be compartmentalized into zones. Groundwater mixing following the flow paths and topography is complicated by the lithologic complexity. An uncommon, highly depleted stable isotope and zero‐³H groundwater, observed in a nearly east–west stretch through the central sector of the catchment, is coincident with the Filwoha Fault zone. Here, deep circulating meteoric water has lost its isotopic content through exchange reactions with CO₂ originating at deeper sources or it has been recharged with precipitation from a different rainfall regime with a depleted isotopic content. Copyright © 2006 John Wiley & Sons, Ltd.