Regional Groundwater System Research Articles

Model simplification to simulate groundwater recharge from a perched gravel-bed river

Gravel-bed rivers are an important source of groundwater recharge in some regions of the world. Their interactions with groundwater are complex and highly variable in space and time, with considerable water storage in the riverbed sediments. In losing river sections, where most of the groundwater recharge occurs, the river can be separated from the regional groundwater system by an unsaturated zone (i.e., perched). The complexity of groundwater–surface water interactions in these environments calls for the use of 3D fully integrated hydrological models to represent them, but their computational intensity limits their practicality for parameter inference, uncertainty quantification and regional scale problems. On the other hand, the simple groundwater–surface water exchange functions currently implemented in regional scale groundwater models are not suited to represent complex gravel-bed river systems such as braided rivers. There is therefore a need for developing groundwater–surface water exchange functions tailored to gravel-bed rivers that can be used in regional scale models.To address this issue, we developed a model simplification framework that combines a 3D integrated surface and subsurface hydrological model, a 2D cross-sectional river-aquifer model and a 1D conductance-based analytical model. We aim at broadly simplifying the 3D model while ensuring the appropriate simulation of groundwater recharge. We demonstrate our modelling approach on the Selwyn River (New Zealand) using piezometric data and groundwater recharge estimates derived from field observations and satellite imagery.Our results indicates that groundwater recharge from this river can be simulated using a simple 1D analytical model, which can easily be implemented in regional groundwater models (e.g., MODFLOW models). However, to represent properly the time variability of groundwater recharge, it is essential to use the groundwater level in the shallow aquifer associated with the river as input to the regional groundwater model. Our approach is generally transferable to other gravel-bed rivers but requires some observations of river losses for proper calibration.

Mobility and persistence of pesticides and emerging contaminants in age-dated and redox-classified groundwater under a range of land use types

Understanding groundwater contamination patterns is hampered by the heterogeneous groundwater age and redox status over the depth range typically sampled for identifying pesticides and emerging contaminants threats. This study explores depth patterns of groundwater age and redox status across various land use types, unraveling spatial and temporal trends of pesticides and emerging contaminants using data from groundwater quality monitoring in the south of the Netherlands. The Netherlands is an ideal testing ground due to its high population density and widespread groundwater contamination from multiple sources. 146 multi-level observation wells were age-dated using 3H/3He, and contaminant concentrations were analyzed based on recharge year, land use type, and redox conditions, mitigating uncertainties from spatial and depth-dependent variations in both groundwater age and redox status. Redox-recharge year diagrams were developed to visually evaluate contaminant patterns in relation to these factors and to assess concentration patterns in relation to contamination history. Most detections of pesticides, metabolites, and emerging contaminants occurred in the youngest recharge periods (2000–2010 and 2010–2020) and in agricultural areas. However, certain contaminants, including BAM, desphenyl-chloridazon, short-chain PFCAs, PFOA, and EDTA, were consistently found in older water and Fe- or SO4-reduced conditions, indicating their mobility and persistence in the regional groundwater system. Comparing the presence of contaminants in specific redox classes and recharge periods with known application or leaching history provides insights into retardation (e.g., PFOS) and degradation (e.g., 2-hydroxy-atrazine, benzotriazole), explaining lower detection frequencies in earlier recharge periods. Identifying recharge years from age-dated groundwater helps relate contaminants to farmland application or river water recharge periods, revealing leaching history and contamination origins. The presented framework has the potential to enhance the interpretation of large groundwater datasets from dedicated, short-screened observation wells, such as those from the Danish GRUMO network, the Dutch monitoring networks, and parts of the US National Water Quality Program.

Deep desert aquifers as an archive for Mid- to Late Pleistocene hydroclimate: An example from the southeastern Mediterranean

Many efforts have been made to illuminate the nature of past hydroclimates in semi-arid and arid regions, where current and future shifts in water availability have enormous consequences on human subsistence. Deep desert aquifers, where groundwater is stored for prolonged periods, might serve as a direct record of major paleo-recharge events. To date, groundwater-based paleoclimate reconstructions have mainly focused on a relatively narrow timescale (up to ∼40 kyr), limited by the relatively short half-life of the widely used radiocarbon (5.73 kyr). Here we demonstrate the usage of deep regional aquifers in the arid southeastern Mediterranean as a hydroclimate archive for earlier Mid-to-Late Pleistocene epochs. State-of-the-art dating tools, primarily the 81Kr radioisotope (t1/2 = 229 kyr), were combined with other atmosphere-derived tracers to illuminate the impact of four distinguishable wetter episodes over the past 400 kyr, with differences in climatic conditions and paleo-recharge locations. Variations in stable water isotope composition suggest moisture transport from more proximal (Mediterranean) and distal (Atlantic) sources to different parts of the region at distinct times. Large variability in the computed noble gas-based recharge temperature (NGT), ranging ~15–30 °C, cannot be explained by climate variations solely, and points to different recharge pathways, including geothermal heating in the deep unsaturated zone and recharge from high-elevation (colder) regions. The obtained groundwater record complements and enhances the interpretation of other terrestrial archives in the arid region, including a contribution of valuable information regarding the moisture source origin as reflected in the deuterium-excess values, which is unattainable from the common practice analysis of calcitic cave deposits. We conclude that similar applications in other deep (hundred-m-order) regional groundwater systems (e.g., the Sahara desert aquifers) can significantly advance our understanding of long-term (up to 1 Myr) paleo-hydroclimate in arid regions, including places where no terrestrial remnants, such as cave, lake, and spring sediments, are available.

Fluid upwelling and alluvial controls on spring localization: An example from Sri Lanka

In water-scarce regions, perennial springs can be a valuable source of drinking water. However, to identify unreported springs and shallow water upwelling zones, it is essential to understand the factors that control spring localization. In a crystalline basement, as in Sri Lanka, without a sedimentary cover, faults and fractures provide the only far-reaching fluid pathways and springs commonly emerge at fault/fracture intersections. While surveying cold and hot water springs in Sri Lanka, it was observed that all springs probed were located at the edge of alluvium. In order to gain insight into this relationship, we performed a topographic and geomorphological analysis was conducted utilizing remote sensing, geological and soil maps, and geological mapping in the field. The results of our analysis of 27 springs indicate that their localization is controlled by fault intersections, non-permeable clay in the alluvium and laterite, and the chemically weathered surface of the bedrock. Furthermore, the constant discharge rates observed over the years and isotope analysis suggest that the springs are part of a tens-of-kilometer-wide regional groundwater system. Based on these results, we propose a conceptual model in which water rises at fault intersections from depth until it reaches the base of the alluvium where up to several meters thick clay with low to zero permeability further inhibits vertical flow forcing the water to spread laterally. Along the alluvium clay boundary with the more permeable weathered bedrock, the water continues its path to the surface. The localization of springs differs from that of fault intersection by tens of meters, with the potential for mixing between shallow and deep groundwater. This observed effect of alluvium and their contact boundaries on spring localization has not been reported for Sri Lanka. Consequently, discharge rates may be significantly increased if the fault intersections are specifically targeted by shallow drilling.

Geophysical and physicochemical assessment of groundwater and the implication on the public health in rural and semi-urban areas of Northern Nigeria

Groundwater resource plays a central role in every sector of life globally. Unfortunately, the roles of groundwater are often unappreciated in many parts of the world, especially in developing nations. The paper applied a combined geo-electrical method and hydrochemical analysis alongside physical surveys to investigate aquifer promising zones, the quality and the source of water supply in the selected regions. The geo-electric resistivity data were used to select the most aquifer-promising zone for groundwater development with a key interest in groundwater potential and aquifer overburden protective capacity overlying the water-bearing unit. Seventeen (17) vertical electrical sounding (VES) acquired in four regions reveals that the regional weathered/fractured unit which constitutes the aquifer unity extends from 13.4 m to 41.3 m in depth with the subsurface resistivity value varying between 53 Ωm–1013 Ωm in terms of resistivity values. The overburdened thickness (protective soil layers) stretches from 2.5 m to 13 m with a resistive layer of resistivity values ranging from 111 Ωm to 2877 Ωm, which could provide high aquifer protective cover for the regional groundwater systems. The computed aquifer transmissivity and hydraulic conductivity range from 0.54 m2/h\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${m}^{2}/h$$\\end{document}–15.85 m2/h\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${m}^{2}/h$$\\end{document}, and 0.95 m/day–9.19 m/day, respectively. The estimated average longitudinal conductance (0.2 moh) reveals that the regional groundwater systems are moderately protective against surface contaminants. In addition, the physicochemical water analysis was examined, which exposes the level of heavy metals (HMs) compositions in the regional groundwater systems and the possible public health risks compared to WHO permissible drinking water limits. The sample water from the Rivers had the highest concentration of HMs (0.089 mg/I–2.012 mg/l), followed by the hand-dug wells (HWDs) (0.033 mg/l–1.555 mg/l) and the lowest in the borehole (0.011 mg/l–0.511 mg/l). In contrast, the borehole water sample shows a moderation where the HMs presence falls within the WHO standard for drinking water. The outcomes of the study show that 63% of the investigated region consumed water with a high concentration of HMs above the WHO permissible standard allowance for drinking water, while only 37% took water that falls within the WHO permissible standard. In conclusion, the study showed significant impacts of anthropogenic activities on the regional water quality, which provides vital information about the regional groundwater systems and the quality of water supply for adequate groundwater development and sustainability.

Complex Effects of Tides on Coastal Groundwater Revealed by High‐Resolution Integrated Flow Modeling

AbstractRiver deltas typically have high population density and support a wide range of intensive and prosperous socioeconomic activities. The hydrological processes in these regions are complex, primarily due to the interactions among the river, aquifer, and sea. However, a systematic and quantitative elaboration of the river‐aquifer‐sea interactions is still lacking. Here, we developed an integrated hydrological flow model for the Pearl River Delta (PRD), which contains the world's largest urban area in both size and population, to gain a deeper understanding of the complexities in the river‐aquifer‐sea interactions. The model performance was validated and cross‐checked via observations at gauging stations and independent remote‐sensing products (e.g., soil moisture, evapotranspiration, and total water storage anomalies). Based on the 10‐year simulation results (2004–2013), the major findings of this study are as follows: (a) accurate representation of the tidal effect is important not only for simulating short‐term flow dynamics but also for capturing the characteristics of long‐term hydrological fluxes and states; (b) the flow‐model‐computed average groundwater discharge rate per unit length of the coastline for the PRD is 3.01 m3/d/m, which is comparable with those derived from water budget approaches but 1–2 orders of magnitude lower than the total submarine groundwater discharge (SGD) estimated by using isotope tracer‐based methods; (c) the temporal variation of SGD is controlled by tidal forcing on an hourly time scale, but by terrestrial hydrological processes on monthly and annual time scales; and (d) an integrated hydrological flow model can be used to identify distinct and large subsurface zones sensitive to tidal fluctuations, quantifying the pivotal role of ocean tides in shaping the coastal groundwater system. This study represents a first step in using an integrated hydrological model to explore river‐aquifer‐sea interactions and their effects on the regional groundwater system simultaneously driven by meteorological and tidal forcings.

Control of Structural Landform Evolution on Karst Groundwater Cycle in a Large-Scale Anticlinorium.

Structural landform evolution and hydrogeochemical analyses are crucial for understanding the characteristics of karst groundwater systems and the development of deep karst formed by complex aquifers in a tectonic collision zone. Detailed structural landform evolution analysis was carried out along the large-scale anticlinorium to investigate the temporal evolution of karst aquifer systems and karstification. Results showed that the tectonic activity included weak horizontal compression and slow vertical uplift during the Triassic to Middle Jurassic, forming a denuded clastic platform. This period was mainly preserved in the geological record as burial karst. From the Late Jurassic to the Early Cretaceous, the study area was strongly compressed by S-N-trending stress, and developed E-W-trending high-angle imbricate thrust structures, which controlled the formation of folded and fault-blocked mountains. Vertical multilayered strata underwent a strong horizontal extrusion, forming a large-scale anticlinorium with secondary folds and faults. With the exposure of carbonate rocks due to rapid crustal uplift, karst began to develop, forming a vertical multilayer karst aquifer system and controlling the distribution of karst groundwater. The Fangxian faulted basin was formed from the Late Cretaceous to the Paleogene, whereby landforms were dominated by intermountain basins. Slow crustal uplift caused the retreat of the denudation line to the east, leading to an increase in hydrodynamic conditions and karstification, and the initiation of early karst groundwater systems. Since the Neogene, intermittent and rapid crustal uplift has led to the deepening of rivers, resulting in the formation of peak clusters and canyons, the development of deep karst, and the complete formation of karst groundwater systems. Combined with hydrogeochemical and borehole data, local, intermediate, and regional karst groundwater systems were identified. It has vital significance to the geological route selection or construction of deep-buried tunnels and the utilization of karst groundwater.

Response of hydrogeological processes in a regional groundwater system to environmental changes: A modeling study of Yinchuan Basin, China

The sustainable development of groundwater resources in arid and semi-arid regions is a challenging task hindered by climate change and human activities. The rational utilization and management of groundwater resources is, therefore, dependent on an understanding of the influences of human and climatic factors on the spatial distribution of groundwater resources and their change over time. The thick Quaternary aquifers in the Yinchuan Basin, China were used herein as an example of how to quantitatively assess spatial and temporal trends in groundwater resources in response to human activities and climate change. A 3D transient groundwater flow model was constructed and used to simulate the evolution and spatial variability of hydrogeological processes from 1990 to 2020. By subsequently applying regime shift detection and correlation analysis to the simulation results, we found that: 1) groundwater storage was continuously depleted over the 30-year period, reaching a cumulative depletion of 1.89×109 m3; 2) human activities were mainly responsible for variations in regional hydrogeological processes for a period of up to 30 years. Climate only affected short-term interannual fluctuations in groundwater storage; 3) human activities (e.g., river water diversion and groundwater abstractions) were the decisive factors causing a continuous reduction of groundwater resources. A policy-driven reduction in water diversion from the Yellow River directly led to a significant drop in groundwater storage, which had a consequent effect on surface water and groundwater interactions and altered agricultural irrigation patterns (crop patterns and irrigation methods); 4) the amount of groundwater recharge from the Yellow River and local lakes increased from 1990 to 2020, whereas the discharge of groundwater to the Yellow River and lakes decreased.

Development of subsurface drainage systems: Discharge – retention – recharge

Sufficient freshwater is needed for water dependent sectors such as agriculture, nature, drinking water, and industry. However, even in low-lying, flood prone countries like the Netherlands, climate change, weather extremes, economic growth, urbanization, land subsidence and increased food production will make it more complex to guarantee sufficient freshwater for all sectors. Furthermore, the frequency and amplitude of extremely dry and wet weather conditions is expected to increase. The current Dutch water management system is not designed to anticipate these extremes. Over the last decades, drained Dutch agricultural fields, land consolidation and urbanization resulted in declining groundwater tables. Additionally, the fresh water demand of different sectors (agriculture, industry, drinking water) increased, causing an increased pressure on the regional groundwater system. As a consequence, the annual groundwater table in sandy soil areas dropped over time with the effect that, nowadays, fresh water is becoming scarce in dry periods. In this paper we provide insight in the shifting water management strategy in the Netherlands (1950–2020), with the corresponding drainage systems, developing from conventional drainage (approx. 1950–1990), to controlled drainage (1990’s onwards), climate adaptive drainage (2010 onwards) and subirrigation systems (2018 onwards). Furthermore, we provide insight in the effect of subirrigation on groundwater levels and crop yields, based on both international literature and measurements of Dutch field pilots. Although subirrigation can contribute to improved soil moisture conditions for crop growth on field scale, we show that the water volume needed for subirrigation can be large and could put a significant pressure on the available regional water sources. Therefore, efficient and responsible use of the available external water sources for subirrigation (e.g. surface water, treated waste water, or groundwater) is required. Finally, the implementation of controlled drainage with subirrigation asks for correct implementation in the regional balance: it requires an integral, catchment-wide approach.

Impact of major nearshore land reclamation project on offshore groundwater system

Offshore groundwater (OG) systems are widely distributed in the coastal areas, where nearshore land reclamation projects are constructed for seaports, airports, and urban other infrastructures. However, there is almost no study on the interaction between these projects and the OG system. This paper presents the first such an effort using the Hong Kong International Airport near Lantau Island, as a case study by developing a three-dimensional density-dependent flow and transport model to investigate the interaction between the reclaimed airport and the OG system. Results show that the offshore relatively fresh groundwater may flow from the nearby Lantau Island through the submarine confined aquifer. The higher pore water pressure in the basal aquifer beneath the airport is generated by the regional groundwater system of the nearby Lantau Island. Besides, the reclaimed airport also serves as extra aquifer space for the storage of fresh groundwater resources. The expansion of the 3rd runway will further increase the water table and move the center of the water mound toward the 3rd runway area. The groundwater system in the airport will reach a new steady state 50 years after the expansion, and the volume of the fresh groundwater lens is calculated to be ~5.2 × 107 m3 with a maximum water table of 1.4 m above mean sea level. The findings of this study have implications on the large-scale nearshore land reclamation which may be affected by the OG system in other coastal areas.

A novel analytical model for the transit time distributions in urban groundwater systems

With the ongoing rapid urbanization across the globe, its interference with groundwater resources is critical to freshwater sustainability. The groundwater transit time distribution (TTD) lumps the flow and transport processes of a regional groundwater system and therefore characterizes the aquifer’s resilience to nonpoint-source contamination. However, the influence of large-scale urban areas on the regional groundwater TTD is not clear. This study proposed a novel analytical model for groundwater TTDs accounting for the effects of impervious urban structures. After the verification against results from particle tracking, we apply this analytical expression to investigate how the position and spatial extent of the urban area change the TTDs from the pre-urban ones. The sensitivity analysis suggests that urban areas tend to increase both the mean and the variance of groundwater transit times. Moreover, modeling results highlight the importance of the spatial relationship between the urban area and the aquifer in determining the urban groundwater TTDs. For aquifers intersected by a local urban area, mean transit time (MTT) is dominated by the horizontal extent of the urban area (i.e., the impervious area), whereas for aquifers intersected by a regional urban area, MTT is strongly controlled by the vertical extent of the urban area in addition to its horizontal size. Being computationally efficient, the proposed analytical model can aid decision-making in urban freshwater resources management and urban planning.

Groundwater, biodiversity, and the role of flow system scale

AbstractGroundwater‐dependent ecosystems and species (GDEs) are found throughout watersheds at locations of groundwater discharge, yet not all GDEs are the same, nor are the groundwater systems supporting them. Groundwater moves along a variety of flow paths of different lengths and with different contributing areas, ranging from shorter local flow paths with low discharge and large seasonal variability to streams, springs and wetlands to longer regional flow paths with potentially larger discharge and low seasonal variability, commonly at low basin elevations. How does this variation in physical hydrology affect the type and distribution of GDEs? Using data on hypsographic position, groundwater‐dependent species distributions, groundwater pumping and streamflow from Oregon, USA, we provide a conceptual model and initial supporting evidence demonstrating that spatial variation in groundwater flow path scales, illustrated using basin hypsography, is a driver of non‐random distribution of GDEs across watersheds. Further, we posit that the spatial variation in primary stressors to groundwater (e.g. pumping and climate change) will differentially affect GDEs depending on their hypsographic position. Furthermore, because of their use for irrigation and municipal water supply, regional groundwater systems and associated species are more likely to be studied and receive regulatory protection. Our initial data point to a disproportionate focus on larger discharge, lower elevation GDEs, which leads to a bias in our understanding of the full suite of biodiversity associated with groundwater discharge as well as their stressors and potential mechanisms for protection.

Characterization of Regional Groundwater System Based on Aquifer Response to Recharge–Discharge Phenomenon and Hierarchical Clustering Analysis

The investigations of groundwater hydrograph reasonably reflect the aquifer response to recharge–discharge phenomenon and its characteristics. A better understanding of aquifer characteristics such as regional aquifer classification, recharge and discharge patterns, aquifer geology and flow patterns are the surface indicators that may be more effective and less costly for interpreting basic regional hydrogeological conditions and assessments. This study deals with the application of Hierarchical Clustering Analysis to understand the groundwater spatio-temporal patterns and to visualize/classify the nature of the aquifer in the regional area of Kaohsiung City, Taiwan. Groundwater level fluctuation patterns and slopes of rising and recession limbs are used to identify the pumping effects and classify aquifers. The results of clustering analysis show that the groundwater observation wells in the study area can be divided into five major characteristics along with the upstream to downstream of Kaoping River. The clusters are consistent with basic lithology distribution and age of sedimentary, which represents the characteristics of groundwater level fluctuation. The identified groundwater hydrographs patterns provide newer insights related to aquifer response to recharge–discharge phenomenon, types of aquifers and their behaviors. The knowledge of water level fluctuations in the observation wells provides a piece of prior information about the abstraction of groundwater. The proposed aquifer classification and pumping effect have great potential for applied use in groundwater management e.g., save drilling cost.

Analysis of coastal groundwater hydrochemistry evolution based on groundwater flow system division

Groundwater hydrochemistry evolution in carbonated coastal aquifers is impacted by seawater intrusion. A well-documented groundwater flow system division is the premise and the key of hydrochemistry evolution or seawater intrusion research under complex geological and topographical conditions. In this research, the regional groundwater system in the study area is divided into three groundwater flow systems, according to topographic, geological, and hydrodynamic conditions. The hydrogeochemistry evolution in the Daweijia coastal aquifer is characterized on the basis of the groundwater flow system division. The results show that the groundwater hydrochemistry evolution varies in different groundwater flow systems in the Daweijia area: (1) Hydrochemistry evolution in groundwater flow system I (SI) is controlled by lateral recharge, dissolution of evaporates, and seawater intrusion. The hydrochemistry presents characteristic HCO3-Ca·Mg, Cl·HCO3-Ca, and Cl-Na·Ca along the runoff path, and seawater intrusion is found 1.3 km inland. (2) Hydrochemistry composition in groundwater flow system II (SII) is influenced by water–rock interaction, freshwater-seawater mixing, evaporation, and reverse cation exchange. Additionally, human activities are potentially a primary factor in the hydrochemical evolution of SII and cause NO3 type groundwater. Seawater intrusion is found 2.8 km inland, with a mixing proportion of < 5%. (3) Hydrochemistry evolution in groundwater flow system III (SIII) is influenced by seawater intrusion, local precipitation, and a short runoff path. A classic seawater intrusion model is present in SIII. The hydrochemical distribution presents a regular changing pattern following the order: HCO3-Ca and HCO3·Cl-Ca (freshwater) and Cl-Mg·Ca/Na·Ca (brackish water). The degree of seawater intrusion in SIII is more serious than that in the SI and SII. Isotopic results generally reveal that the groundwater is primarily derived from the combination of local precipitation and lateral recharge from the system boundary. The study proposes an effective method to analyze groundwater hydrochemical evolution under complex geological conditions in coastal areas.

GRACE-based high-resolution propagation threshold from meteorological to groundwater drought

Groundwater drought could cause tremendous damage to the social-economy via land subsidence, seawater intrusion and permanent loss of aquifer storage capacity, and often show strong association with meteorological drought. To date, the threshold for meteorological drought triggering groundwater drought and its dominant factors have been not clarified, which inhibits the effective groundwater drought risk management based on preceding meteorological drought information. In this study, we used the Standardized precipitation index (SPI) and the drought severity index of groundwater storage anomalies (GWSA-DSI) to characterize meteorological and groundwater droughts in the Xijiang River Basin (XRB) of China, respectively. A probabilistic framework is proposed to identify the high-resolution propagation thresholds from meteorological to groundwater drought on 0.25° grid. Results show that GWSA-DSI can reliably identify groundwater drought events, and the propagation time from meteorological to groundwater drought ranges from 8 to 42 months. Although the XRB is located in a humid region with abundant precipitation, the probability of groundwater drought occurrence reached 43.8%, 54.8%, 61.2%, and 64.2% under a light, moderate, severe and extreme meteorological drought event, respectively. The propagation threshold triggering light groundwater drought is mainly dominated by moderate and severe meteorological droughts, which showed an increasing trend from central to southeast of XRB. Soil evaporation and watershed elevation are the main influencing factors on the propagation threshold. It is worth noting that anthropogenic overexploitation of groundwater not only destroy the dynamic balance of regional groundwater system, but also interfere with the propagation processes of meteorological to groundwater drought. The results have great implications for more reliably monitoring and predicting the dynamics of groundwater systems under drought stress, and our proposed framework can also be extended to other regions.

Induced polarization effects in airborne transient electromagnetic data collected in the McMurdo Dry Valleys, Antarctica

SUMMARYAirborne electromagnetics (EM) is a geophysical tool well suited to mapping glacial and hydrogeological structures in polar environments. This non-invasive method offers significant spatial coverage without requiring access to the ground surface, enabling the mapping of geological units to hundreds of metres depth over highly varied terrain. This method shows great potential for large-scale surveys in polar environments, as common targets such as permafrost, ice and brine-rich groundwater systems in these settings can be easily differentiated because of their significant contrasts in electrical properties. This potential was highlighted in a 2011 airborne EM survey in the McMurdo Dry Valleys that mapped the existence of a large-scale regional groundwater system in Taylor Valley. A more comprehensive airborne EM survey was flown in November 2018 to broadly map potential groundwater systems throughout the region. Data collected in this survey displayed significant perturbations from a process called induced polarization (IP), an effect that can greatly limit or prevent traditional EM workflows from producing reliable geological interpretations. Here, we present several examples of observed IP signatures over a range of conditions and detail how workflows explicitly designed to handle IP effects can produce reliable geological interpretations and data fits in these situations. Future polar EM surveys can be expected to encounter strong IP effects given the likely presence of geological materials (e.g. ice and permafrost) that can accentuate the influence of IP.

Chemical evolution of an inland sabkha: a case study from Sabkha Matti, Saudi Arabia

Sabkha Matti is a potential discharge point for regional groundwater systems in the Rub’ al Khali topographic basin of Saudi Arabia. The hydrogeochemical evolution of this sabkha (salt flat) and the underlying aquifers were assessed by using a combination of geological, hydraulic, hydrochemical, and isotopic approaches. A compilation of the geologic structure and lithology in combination with boron isotope data attributed the origin of salinity in the regional aquifers underlying Sabkha Matti to the entrapment of ancient seawater in a coastal lagoon environment during the depositional time of the late Oligocene-Miocene formations. Major ionic constituents and strontium isotope ratios from the trapped paleo-seawater were observed to resemble those obtained from the near-surface sabkha brine. These data, along with hydrostatic head measurements, suggest that the origin of the solutes in Sabkha Matti is associated with brines ascending from the underlying formations. The data also show that solute concentrations in Sabkha Matti are increasing over time through a combination of evaporation, mineral dissolution by recharge, and a density-driven convection mechanism that circulates the solutes. On the other hand, stable isotopes of water (δ18O, δ2H) and radio isotope data (14C and 3H) suggest that the existing waters in the sabkha and the underlying aquifers are relatively recent and were recharged during the wet phase in the late Pleistocene-Holocene between 3,200 and 14,000 BP. This study shows that the water and solutes in Sabkha Matti area are a combination from different sources.

Strategies to Overcome Intermediate Accumulation During in situ Nitrate Remediation in Groundwater by Hydrogenotrophic Denitrification.

Bioremediation of polluted groundwater is one of the most difficult actions in environmental science. Nonetheless, the clean-up of nitrate polluted groundwater may become increasingly important as nitrate concentrations frequently exceed the EU drinking water limit of 50 mg L–1, largely due to intensification of agriculture and food production. Denitrifiers are natural catalysts that can reduce increasing nitrogen loading of aquatic ecosystems. Porous aquifers with high nitrate loading are largely electron donor limited and additionally, high dissolved oxygen concentrations are known to reduce the efficiency of denitrification. Therefore, denitrification lag times (time prior to commencement of microbial nitrate reduction) up to decades were determined for such groundwater systems. The stimulation of autotrophic denitrifiers by the injection of hydrogen into nitrate polluted regional groundwater systems may represent a promising remediation strategy for such environments. However, besides high costs other drawbacks, such as the transient or lasting accumulation of the cytotoxic intermediate nitrite or the formation of the potent greenhouse gas nitrous oxide, have been described. In this article, we detect causes of incomplete denitrification, which include environmental factors and physiological characteristics of the underlying bacteria and provide possible mitigation approaches.

Multi-method geochemical characterization of groundwater from a hypogene karst system

An approach, combining several geochemical methods, was used to determine the groundwater properties and components of a hypogene karst system, where sampling is restricted only to the spring sites, and with a limited number of available sampling locations. Radiogenic isotopes (3H, 14C) were used to constrain the groundwater mean residence time and separate different groundwater components. Noble gases, stable isotopes of water (δ18O, δ2H), dissolved inorganic carbon (δ13C) and dissolved sulfate (δ34S, δ18O), and major ion and trace element composition were used to identify the source of water, its chemical evolution and water–rock interactions, as well as to identify the contribution and composition of endogenic gases. This approach was applied to three low-temperature thermal springs located in Mariovo (North Macedonia) associated with fossil hypogene caves, previously identified by morphological and geochemical studies of caves and cave deposits. Based on the obtained results, the main studied springs represent an output part of a regional hypogene karst groundwater system with a deep-circulating (~1 km), old (~15 ka), thermal (≥60 °C) water, which mixes with young (<50 years), cold (<14 °C) and shallow epigene karst groundwater. The output parts are structurally controlled, at the interception of low topography and deep faults, along which the groundwater interacts with deep-seated gases, dominantly CO2 of metamorphic origin (δ13C of +4.5‰ VPDB), with some contribution of mantle helium. The thermal karst groundwater interacts at depth with volcanic rocks from the nearby Neogene-Quaternary volcanic complex, as well as with metamorphic basement rocks and granitoids.

Groundwater and Solute Budget (A Case Study from Sabkha Matti, Saudi Arabia)

Sabkha Matti is the largest inland sabkha (2950 km2) in the Arabian Peninsula. The drainage area supporting this sabkha is &gt;250,000 km2 and is the discharge point for part of the ten thousand meter thick regional groundwater systems ranging in age from Precambrian through Miocene in the Rub’ al Khali structural basin. A hydrologic budget was constructed for this sabkha, where water fluxes were calculated on the basis of hydraulic gradient and conductivities measured in both shallow and deep wells. The evaporation rates from the surface of the sabkha were estimated from the published data and indicate that almost all the annual rainfall is lost by surface evaporation. The water flux multiplied by its solute concentration showed that nearly all the solutes in the sabkha were derived by upward leakage from the underlying regional aquifers rather than the weathering of the aquifer framework, from precipitation, or from other sources. Steady-state estimates within a rectilinear control volume of the sabkha indicate that about 1 m3/year of water enters by lateral groundwater flow, 2 m3/year of water exits by lateral groundwater flow, 20 m3/year enters by upward leakage, 780 m3/year enters by recharge from rainfall, and 780 m3/year is lost by evaporation. The proposed conceptual model of the hydrology for sabkha Matti is assumed to apply to the rest of the inland sabkhas of the Arabia Peninsula and to many ancient environments of deposition observed in the geologic record.