Showing papers on "Groundwater flow published in 2021"

Poor groundwater quality and high potential health risks in the Datong Basin, northern China: research from published data.

Abstract: Datong Basin in China is a typical arid–semiarid inland basin, with high levels and wide distributions of arsenic (As), fluoride (F−), and iodine (I). To better understand the presence of low-quality groundwater in Datong Basin and assess the health risks for local residents, groundwater samples were collected from the shallow aquifer and in medium-deep groundwater and analyzed for As, F−, I, and nitrate (NO3−). Maxima of 1932 μg/L for As, 80.89 mg/L for F−, 2300 μg/L for I, and 3854.74 mg/L for NO3− were detected in shallow groundwater, which greatly exceeded the WHO limits for drinking purpose. High-As groundwater was present in both shallow and medium-deep aquifers. High-F− and high-NO3− groundwater was widely distributed in the shallow aquifer, and high-I groundwater was mainly present in the medium-deep aquifers. Poor-quality groundwater in the Datong Basin is mainly caused by local geological and climatic conditions, which are characterized by strong evaporation, active water–rock interactions, thick lacustrine sediment, low groundwater flow rate, and reducing and weak alkaline environments. However, groundwater quality was further impacted by agricultural activities in some areas, as shallow groundwater was also polluted by nitrate. Datong Basin inhabitants face high health risk caused by high concentrations of As, F−, I, and NO3−. The mean noncarcinogenic risk values (HQtotal) were 18.40 for children, 10.94 for adult females, and 9.47 for adult males due to exposure to contaminants in shallow groundwater; and 13.76 for children, 8.18 for adult females, and 7.08 for adult males because of exposure to medium-deep groundwater. Further, the carcinogenic risks (CR) caused by exposure to As were very high for local inhabitants, with the mean and median CR values of 4.20×10−3 and 4.13×10−4 in shallow groundwater and 3.44×10−3 and 1.71×10−4 in medium-deep groundwater, respectively.

Finite Difference Modelings of Groundwater Flow for Constructing Artificial Recharge Structures

Abstract: Rainfall and deep percolation of irrigation are the most important sources of recharge in Khatoon-Abad plain, Kerman province, Iran. Artificial recharge is a practical strategy to increase the storage capacity of groundwater in drought condition. Modeling of groundwater flow for a steady-state period (April 2017 to October 2017) and an unsteady state period from October 2012 to August 2019 was simulated in monthly stress steps. In this study, root mean square error and coefficient of residual mass have been used to evaluate the results of calibration and verification process obtained by PMWIN software. Moreover, two strategies of artificial recharge were incorporated for considering the water table fluctuations. The results show that artificial recharge by direct injection of water into the aquifer prevents evaporation, runoff and outflow as runoff from the plain and increases the aquifer water availability. In different scenarios, the different amounts of recharge were applied to the aquifer, and the simulation results showed that the effect of recharge in the study period varies from 3 km to more than 7 km.

Application of physics-informed neural networks to inverse problems in unsaturated groundwater flow

Abstract: This paper investigates the application of Physics-Informed Neural Networks (PINNs) to inverse problems in unsaturated groundwater flow. PINNs are applied to the types of unsaturated groundwater fl...

Characterization of Groundwater Recharge and Flow in California's San Joaquin Valley From InSAR-Observed Surface Deformation.

Abstract: Surface deformation in California's Central Valley (CV) has long been linked to changes in groundwater storage. Recent advances in remote sensing have enabled the mapping of CV deformation and associated changes in groundwater resources at increasingly higher spatiotemporal resolution. Here, we use interferometric synthetic aperture radar (InSAR) from the Sentinel-1 missions, augmented by continuous Global Positioning System (cGPS) positioning, to characterize the surface deformation of the San Joaquin Valley (SJV, southern two-thirds of the CV) for consecutive dry (2016) and wet (2017) water years. We separate trends and seasonal oscillations in deformation time series and interpret them in the context of surface and groundwater hydrology. We find that subsidence rates in 2016 (mean -42.0 mm/yr; peak -345 mm/yr) are twice that in 2017 (mean -20.4 mm/yr; peak -177 mm/yr), consistent with increased groundwater pumping in 2016 to offset the loss of surface-water deliveries. Locations of greatest subsidence migrated outwards from the valley axis in the wetter 2017 water year, possibly reflecting a surplus of surface-water supplies in the lowest portions of the SJV. Patterns in the amplitude of seasonal deformation and the timing of peak seasonal uplift reveal entry points and potential pathways for groundwater recharge into the SJV and subsequent groundwater flow within the aquifer. This study provides novel insight into the SJV aquifer system that can be used to constrain groundwater flow and subsidence models, which has relevance to groundwater management in the context of California's 2014 Sustainable Groundwater Management Act (SGMA).

Review of Modeling Approaches to Groundwater Flow in Deformed Carbonate Aquifers

Abstract: We discuss techniques to represent groundwater flow in carbonate aquifers using the three existing modeling approaches: equivalent porous medium, conduit network, and discrete fracture network. Fractures in faulted stratigraphic successions are characterized by dominant sets of sub-vertical joints. Grid rotation is recommended using the equivalent porous medium to match higher hydraulic conductivity with the dominant orientation of the joints. Modeling carbonate faults with throws greater than approximately 100 m is more challenging. Such faults are characterized by combined conduit-barrier behavior. The barrier behavior can be modeled using the Horizontal Flow Barrier Package with a low-permeability vertical barrier inserted to represent the impediment of horizontal flow in faults characterized by sharp drops of the piezometric surface. Cavities can occur parallel to the strike of normal faults generating channels for the groundwater. In this case, flow models need to account for turbulence using a conduit network approach. Channels need to be embedded in an equivalent porous medium due to cavities a few centimeters large, which are present in carbonate aquifers even in areas characterized by low hydraulic gradients. Discrete fracture network modeling enables representation of individual rock discontinuities in three dimensions. This approach is used in non-heavily karstified aquifers at industrial sites and was recently combined with the equivalent porous medium to simulate diffusivity in the matrix. Following this review, we recommend that the future research combines three practiced modeling approaches: equivalent porous medium, discrete fracture network, and conduit network, in order to capture structural and flow aspects in the modeling of groundwater in carbonate rocks.

Groundwater geochemical signatures and implication for sustainable development in a typical endorheic watershed on Tibetan plateau

Abstract: Groundwater resource is significantly important for sustainable development of the world, especially for arid endorheic watersheds. A total of 28 groundwaters were collected for hydrogeochemical analysis from the arid Chaka watershed on Tibetan plateau to illustrate the hydrochemical evolution, formation mechanisms and feasibility of groundwater in small arid endorheic watersheds where groundwater is much scarcer. The results showed groundwater has a slightly alkaline nature, and varies from soft fresh HCO3-Ca type to hard brackish/saline Cl-Na type along the groundwater flow path in the watershed with the total hardness in the range of 270–2,127 mg/L and the total dissolved solids in the range of 282–41,770 mg/L. Nitrogen and fluoride in phreatic water are found sporadically exceeding the permissible limits with the maximum value of 118 mg/L for nitrate, 1.2 mg/L for ammonia and 1.2 mg/L for fluoride. Hydrochemistry of phreatic and confined groundwater is naturally governed by water-rock interactions including minerals (halite, gypsum and anhydrite) dissolution, silicate weathering and cation-exchange reaction. The salinity of phreatic water is also dominantly controlled by the strong evaporation. Human activity is one of the important mechanisms influencing the hydrochemical signature of groundwater regardless of the depth. Groundwater has a great hydrogeochemical discrepancy spatially across the watershed and varies from excellent to extremely poor quality in phreatic aquifers. A better water quality that under the good to medium categories was observed in the confined aquifers with 80% of samples having the EWQI value less than 100 and others in the range of 100–150. Phreatic groundwater away from the river and in the downstream area has a relatively poor quality for domestic and agricultural purposes, and should be avoided to direct utilization. This research can improve the understanding of groundwater hydrogeochemical feature, genesis, and its constraints on the availability and feasibility of groundwater resources in small arid watersheds worldwide.

Investigation of groundwater occurrence using gravity and electrical resistivity methods: a case study from Wadi Sar, Hijaz Mountains, Saudi Arabia

Abstract: An integrated approach of Gravity Recovery and Climate Experiment (GRACE) and vertical electrical resistivity sounding (VES) technique has been carried out to investigate regionally and locally the groundwater potentialities of Wadi Sar in the Hijaz Mountains. Our findings are (1) the terrestrial water storage variations (ΔTWS) are estimated at −2.06±0.34 mm/year; (2) the Global Land Data Assimilation System–derived soil moisture storage variations (ΔSMS) are estimated at −0.067±0.005 mm/year; (3) the groundwater storage variations (ΔGWS) show a negative trend estimated at −2.00±0.34 mm/year during the period April 2002–July 2017; and (4) the average annual precipitation (AAP) rate is estimated at 115 mm during the period 2002–2018. Three geoelectrical layers are identified from the inversion of the electrical resistivity data: (5) the surface layer of high resistivity values is consisted of dry unconsolidated Quaternary deposits; (6) the second layer represents the fractured groundwater aquifer of low resistivity values and variable thickness; (7) the third layer is composed of fractured basement rocks of higher resistivity values. The streams are draining the surface water toward the Najd Pedi plain aquifer; (8) the faults are acting as conduits for groundwater flow away from the Wadi. The current study indicates the occurrence of groundwater at the downstream zones of the Wadi Sar, but with a general decrease trend. The integrated approach provides a better understanding of the groundwater potentialities in the arid regions.

Hydrochemical characteristics and functions of groundwater in southern Laizhou Bay based on the multivariate statistical analysis approach

Abstract: Influenced by paleoclimate, paleogeography and hydrogeology, groundwater with different total dissolved solids (TDS) levels, such as freshwater, brackish water, saline water and brine, widely occurs across southern Laizhou Bay (SLB) . Due to groundwater overexploitation, several groundwater depression cones have gradually formed since the 1970s, resulting in saline water intrusion. In this study, the multivariate statistical analysis (MSA) approach was adopted to identify the hydrochemical characteristics of shallow groundwater (SG) and deep groundwater (DG) in the SLB. Moreover, we applied Gibbs diagrams, hydrochemical facies evolution diagrams (HFE-Diagram) and ion proportion relations to analyze the factors influencing the hydrochemical characteristics. The results showed that hydrochemical zonation of the various water types (ranging from Ca-HCO3 to Na-Cl) from south to north was mainly controlled by migration and mixing of saline water with regional groundwater. Based on the MSA results and hydrochemical characteristics, SG functional areas were divided into the groundwater flow zone, intense groundwater mixing zone and groundwater-seawater interaction zone, and DG functional areas were divided into the groundwater flow zone, groundwater mixing zone and underground brine storage zone. Finally, SG and DG groundwater flow evolution models were developed . The comprehensive method proposed in this study combines the MSA approach with regional geological characteristics and provides a reference for the evaluation of complex multiaquifer systems and analysis of the groundwater evolution characteristics in other regions.

Numerical and Experimental Validation of the Applicability of Active-DTS Experiments to Estimate Thermal Conductivity and Groundwater Flux in Porous Media

Abstract: Groundwater flow depends on the heterogeneity of hydraulic properties whose field characterization is challenging. Recently developed active‐Distributed Temperature Sensing (DTS) experiments offer the possibility to directly measure groundwater fluxes resulting from heterogeneous flow fields. Here, based on fundamental principles and numerical simulations, two interpretation methods of active‐DTS experiments are proposed to estimate both the porous media thermal conductivities and the groundwater fluxes in sediments. These methods rely on the interpretation of the temperature increase measured along a single heated fiber optic (FO) cable and consider heat transfer processes occurring both through the FO cable itself and through the porous media. The first method relies on the Moving Instantaneous Line Source (MILS) model that reproduces the temperature increase and provides estimates of thermal conductivity and groundwater flux as well as an evaluation of the temperature rise due to the FO cable. The second method, based on the graphical identification of three characteristic times, provides complementary estimates of flux, fully independent of the effect of the FO cable. Sandbox experiments provide an experimental validation of the interpretation methods, demonstrate the excellent accuracy of groundwater flux estimates (< 5%) and highlight the complementarity of both methods. Active‐DTS experiments allow investigating groundwater fluxes over a large range spanning 1x10‐6 to 5x10‐2 m/s, depending on the duration of the experiment. Considering the applicability of active‐DTS experiments in different contexts, we propose a general experimental framework for the application of both interpretation methods in the field, making active‐DTS field experiments especially promising for many subsurface applications.

PFAS soil and groundwater contamination via industrial airborne emission and land deposition in SW Vermont and Eastern New York State, USA.

Abstract: In order to understand the extent to which airborne PFAS emission can impact soil and groundwater, we conducted a sampling campaign in areas of conserved forest lands near Bennington, VT/Hoosick Falls, NY. This has been home to sources of PFAS air-emissions from Teflon-coating operations for over 50 years. Since 2015, the Vermont and New York Departments of Environmental Conservation have documented ∼1200 residential wells and two municipal water systems across a 200 km2 area contaminated with perfluorooctanoic acid (PFOA). Given the large areal extent of the plume, and the fact that much of the contaminated area lies up-gradient and across rivers from manufactures, we seek to determine if groundwater contamination could have resulted primarily from air-emission, land deposition, and subsequent leaching to infiltrating groundwater. Sampling of soils and groundwater in the Green Mountain National Forest (GMNF) downwind of factories shows that both soil and groundwater PFOA contamination extend uninterrupted from inhabited areas into conserved forest lands. Groundwater springs and seeps in the GMNF located 8 km downwind, but >300 meters vertically above factories, contain up to 100 ppt PFOA. Our results indicate that air-emitted PFAS can contaminate groundwater and soil in areas outside of those normally considered down-gradient of a source with respect to regional groundwater flow.

Prediction of shallow landslides in pyroclastic-covered slopes by coupled modeling of unsaturated and saturated groundwater flow

Abstract: Slopes covered with unsaturated shallow pyroclastic deposits, lying upon fractured limestone bedrock, are widespread in the mountains around Naples (southern Italy). Rainfall infiltration, reducing soil suction, eventually triggers shallow landslides. While drastic reduction of suction is unanimously recognized as the triggering mechanism, there is still debate about the hydrological processes controlling slope drainage and causing the establishment of landslide predisposing conditions. Field observations at the slope of Cervinara suggested that temporary storage of water in a perched aquifer, in the upper part of the fractured bedrock, may affect the leakage through the soil-bedrock interface. Hence, a physically based model, coupling flows in the unsaturated soil cover and in the perched aquifer, has been applied to three large rainfall events which occurred in December 1999 (when a landslide was triggered), January 2009, and November 2012. The results highlight that the different responses of soil and aquifer to precipitations, related not only to rainfall event characteristics (i.e., duration and mean intensity) but also to the initial conditions of the slope, determined by antecedent precipitations, can play a prominent role in the triggering of landslides. In fact, further simulations with synthetic rainfall events and different initial conditions provide a possible interpretation of the triggering of the landslide of December 1999, indicating that a soil profile with dry conditions at the base and a low level in the perched aquifer, typical of late autumn, can impede the drainage of infiltrating water through the soil-bedrock interface, thus favoring the build-up of pore pressure within the soil cover.

Origin and geochemical evolution of groundwater in the Abaya Chamo basin of the Main Ethiopian Rift: application of multi-tracer approaches

Abstract: The fractured volcanic aquifer of the Abaya Chamo basin in the southern Ethiopian Rift represents an important source for water supply. This study investigates the geochemical evolution of groundwater and the groundwater flow system in this volcanic aquifer system using hydrochemistry and environmental tracers. Water types of groundwater were found to transform from Ca-Mg-HCO3 (western part of Lake Abaya area) to Na-HCO3 (northwestern part), from the highland down to the Rift Valley. Silicate hydrolysis and Ca/Na ion exchange are the major geochemical processes that control groundwater chemistry along the flow path. Groundwaters are of meteoric origin. The δ18O and δD content of groundwater ranges from −4.9 to −1.1‰ and –27 to 5‰, respectively. The δ18O and δD values that lie on the summer local meteoric water line indicate that the groundwater was recharged mainly by summer rainfall. δ13CDIC values of cold groundwater range from −12 to −2.7‰, whereas δ13CDIC of thermal groundwater ranges from −8.3 to +1.6‰. The calculated δ13CCO2(g) using δ13CDIC and DIC species indicates the uptake of soil CO2 for cold groundwater and the influx of magmatic CO2 through deep-seated faults for thermal groundwater. In the western part of Lake Abaya area, the shallow and deep groundwater are hydraulically connected, and the uniform water type is consistent with a fast flow of large gradient. In contrast, in the northern part of Lake Abaya area, water underwent deep circulation and slow flow, so the water types—e.g. high F− (up to 5.6 mg/L) and Na+—varied laterally and vertically.

A Basin-Scale Groundwater Flow Model of the Columbia Plateau Regional Aquifer System in the Palouse (USA): Insights for Aquifer Vulnerability Assessment

Abstract: Bedrock aquifers are vulnerable to contamination due to the preferential movement of pollutants via rock discontinuities and porous layers. In this research, we propose an approach to assess vulnerability in three dimensions by combining stable isotope values and particle tracking in a vertically anisotropic aquifer of basaltic-fluvial origin at the basin-scale. A steady-state flow and particle tracking model is presented for the Columbia River Basalt aquifer in the South Fork Palouse Basin. Backward particle analysis combined with the distribution of δ2H, δ18O and 4He values vs. depth shows how the aquifer is characterized by two separated zones. A shallow ( 150 mBGL) aquifer zone is characterized by much higher particle travel times as well as a distinctive isotopic fingerprint. At such depths, penetration of particles is partially impeded by the low hydraulic conductivity of the sedimentary layers and recharge preferentially occurs in correspondence of the basin margin. Along this margin, the vulnerability is higher for the contaminants to enter the aquifer system and reach the pumping wells. Thus, following this research, we envisage efforts to combine stable isotope techniques with particle tracking analysis in three dimensions to define areas exposed to contamination risk in fluvio-volcanic bedrock aquifers. These research efforts can represent an approach to integrate with two-dimensional GIS tools that are commonly used to assess aquifer vulnerability.

A model ensemble generator to explore structural uncertainty in karst systems with unmapped conduits

Abstract: Karst aquifers are characterized by high-conductivity conduits embedded in a low-conductivity fractured matrix, resulting in extreme heterogeneity and variable groundwater flow behavior. The conduit network controls groundwater flow, but is often unmapped, making it difficult to apply numerical models to predict system behavior. This paper presents a multi-model ensemble method to represent structural and conceptual uncertainty inherent in simulation of systems with limited spatial information, and to guide data collection. The study tests the new method by applying it to a well-mapped, geologically complex long-term study site: the Gottesacker alpine karst system (Austria/Germany). The ensemble generation process, linking existing tools, consists of three steps: creating 3D geologic models using GemPy (a Python package), generating multiple conduit networks constrained by the geology using the Stochastic Karst Simulator (a MATLAB script), and, finally, running multiple flow simulations through each network using the Storm Water Management Model (C-based software) to reject nonbehavioral models based on the fit of the simulated spring discharge to the observed discharge. This approach captures a diversity of plausible system configurations and behaviors using minimal initial data. The ensemble can then be used to explore the importance of hydraulic flow parameters, and to guide additional data collection. For the ensemble generated in this study, the network structure was more determinant of flow behavior than the hydraulic parameters, but multiple different structures yielded similar fits to the observed flow behavior. This suggests that while modeling multiple network structures is important, additional types of data are needed to discriminate between networks.

Groundwater flow velocities in karst aquifers; importance of spatial observation scale and hydraulic testing for contaminant transport prediction.

Abstract: We review scale dependence of hydraulic conductivities and effective porosities for prediction of contaminant transport in four UK karst aquifers. Approaches for obtaining hydraulic parameters include core plug, slug, pumping and pulse tests, calibration of groundwater flow models and spring recession curves. Core plug and slug tests are unsuitable because they do not characterize a large enough volume to include a representative fracture network. Pumping test values match regional-scale hydraulic conductivities from flow modelling for the less intensively karstified aquifers: Magnesian Limestone, Jurassic Limestone and Cretaceous Chalks. Reliable bulk hydraulic conductivities were not available for the intensively karstified Carboniferous Limestone due to dominance of flow through pipe conduits in Mendips. Here, the only hydraulic conductivity value found from spring recession is one order of magnitude higher than that indicated by pumping tests. For all four carbonate aquifers, effective porosities assumed for transport modelling are two orders of magnitude higher than those found from tracer and hydrogeophysical tests. Thus, a combination of low hydraulic conductivities and assumed flowing porosities resulted in underestimated flow velocities. The UK karst aquifers are characterized by a range of hydraulic behaviours that fit those of karst aquifers worldwide. Indeed, underestimation of flow velocity due to inappropriate parameter selection is common to intensively karstified aquifers of southern France, north-western Germany and Italy. Similar issues arise for the Canadian Silurian carbonates where the use of high effective porosities (e.g. 5%) in transport models leads to underestimation of groundwater velocities. We recommend values in the range of 0.01–1% for such aquifers.