Porous Aquifer Research Articles

Mechanisms controlling spatial variability of geogenic ammonium in coastal aquifers: Insights from Holocene sedimentary evolution.

The contamination of groundwater with geogenic ammonium (NH4+) across various geological backgrounds has garnered significant attention, particularly in coastal aquifer systems. However, there remains a gap in our understanding of the mechanisms governing the spatial variability of NH4+ in coastal groundwater at a macroscopic scale. In this study, we collected the sediment samples from two boreholes corresponding to high-NH4+-N and low-NH4+-N groundwater. We analyzed the age, physicochemical properties, and soluble organic matter (SOM) characteristics of these sediments. The aim was to reconstruct the sedimentary evolutionary history of the Pearl River Delta and establish a link between sedimentary evolutionary processes and organic matter (OM) to further elucidate the mechanism underlying the formation of spatial heterogeneity of NH4+ in groundwater. The results suggested that the studied Quaternary shallow confined porous aquifer system was predominantly formed during the Holocene and comprised three depositional stages, including fluvial facies, estuarine-tidal flat facies, and deltaic plain facies. The depositional environment significantly controlled the physicochemical and SOM characteristics of sediments. In the paleo-channel area, the aquifer was covered by estuarine-tidal flat facies sediments abundant in OM and exhibited considerable SOM degradation. Consequently, a substantial amount of ion-exchange form NH4+-N (IEF) was liberated through compaction and diffused into the aquifer. In the paleo-interfluve area, the aquifer was covered by fluvial sediments characterized by extensive weathering and low OM content, resulting in the limited production of significant amounts of IEF that could infiltrate into the aquifer. This study provided an inaugural elucidation of the control mechanism of sedimentary evolution on the spatial variability of NH4+ in coastal groundwater at a macroscale, thereby enhancing the scientific substantiation for both the exploitation and protection of groundwater resources.

Modelling approach integrating climate projections for coastal groundwater management

Climate change and excessive groundwater extraction are major contributors to rising groundwater salinization due to seawater intrusion in coastal aquifers. This study aims to define a wide-applicable approach in which hydrological balance, boundary conditions, and irrigation water demand, defined over time considering climate change predictions, can integrated into a numerical model of the groundwater system. The approach was tested in a selected coastal aquifer. The approach spans from the past, used to define steady or almost natural conditions for calibration purposes (1950–2000 in the test), to the future (2100), divided in decade steps. The water balance analysis is based on an inverse hydrogeological water balance approach. The future climate change predictions are used to assess variations in boundary conditions of the groundwater model concerning salinity and sea level, recharge, and inflow from upstream aquifers. The approach considers changes in agricultural activities, groundwater demand, and river stage. The regional model is generated using the MODFLOW code for the groundwater flow model and the SEAWAT code for the salt transport model. The test concerns the Metaponto coastal plain, in which a porous aquifer is at salinization risk due to seawater intrusion. In this way, different influences of climate change and human activities are combined to define a 3d view of groundwater depletion and salinization effects. Quantifying these potential effects or risks, adaptation scenarios with numerical assessments are outlined in this study.

Land surface response to groundwater drawdown and recovery in Taiyuan city, Northern China, analyzed with a long-term elevation change measurements from leveling and multi-sensor InSAR

In recent years, a comprehensive program of groundwater management was adopted in Taiyuan city, Northern China, with the aim to prevent the lowering of groundwater level and manage land subsidence due to excessive groundwater pumping in the past decades. Groundwater recovery has been observed in this city following the termination of groundwater pumping. While the link between land subsidence and groundwater pumping is well documented, how the land surface responds to the recovery of groundwater has not been studied in detail. In this study, we analyzed the leveling data from 1956 to 2000 to investigate changes in ground elevation in response to groundwater drawdown prior to the termination of groundwater pumping and combined them with InSAR (Interferometric Synthetic Aperture Radar) displacement data during 2003–2020 to evaluate the effect of groundwater recharge after the termination of groundwater pumping. A method was proposed to merge InSAR time series of displacement derived from multiple satellites (ENVISAT, COSMO-SkyMed, TerraSAR-X, and Sentinel-1) to achieve a consistent, long-term continuous deformation. We found that previous groundwater depression cones in Taiyuan (Xizhang, Wanbailin, Wujiabao, and Xiayuan) have turned into groundwater recovery and the corresponding subsidence centers have switched to land uplift. The subsidence center of Xizhang, north of Taiyuan city, experienced a total subsidence of 749 mm by Aug 2003; it then turned into uplift with an amount of rebound of 9.2 % of previous subsidence by the end of 2020. In central area of Taiyuan, the cumulated subsidence is 1127 mm in Wanbailin, 1817 mm in Xiayuan, and 3343 mm in Wujiabao. Beginning in Aug 2010, they all turned into land uplift, with a rebound of 7.9 %, 4.8 % and 3.6 % of the previous subsidence in each center, respectively. An analysis of the correlation between groundwater level and ground displacement suggests that the ground uplift is related to the poro-elastic rebound mechanism in the porous aquifer system, which is caused by the rise in pore pressure with groundwater recovery, and the amount of subsidence/uplift is controlled by the stratigraphic properties. The adverse impacts on the built environment due to the rising groundwater level were discussed. The findings improve our understanding of anthropogenic ground deformation in complex urban aquifers influenced by groundwater pumping and recharge and provides essential information that can contribute to future groundwater management and groundwater-related hazard mitigation over the Taiyuan city.

Phenomenal study of microbial impact on hydrogen storage in aquifers: A coupled multiphysics modelling

Underground hydrogen storage (UHS) has been considered as an integral part of energy transition from fossil fuels to renewable sources. Porous aquifers can serve as typical sites for this purpose due to their worldwide distribution and huge storage capacity. However, the diverse microbial species that inhabit the aquifers are known to be able to catalyze in-situ biochemical reactions within the hosting rock. These reactions can normally lead to microbial consumption of hydrogen, microbial clogging of pore space and thus affect hydrogen injection and withdrawal rates as reported in the literature. So far, these phenomena have been widely reported but rarely quantified. In this study, we build a coupled hydrological-mechanical-chemical-biological (HMCB) multiphysics model to simulate these microbe-related processes during UHS in aquifers. The model consists of a complete set of partial differential equations to describe: (1) rock deformation; (2) water-hydrogen two-phase flow; (3) microbes and dissolved hydrogen transport; (4) mineral dissolution/precipitation; and (5) microbial activities involving adsorption/desorption and growth/decay. All these processes are linked together through the porosity/permeability models which consider the joint impacts of microbial clogging, mineral dissolution/precipitation and effective stress. This multiphysics model is verified against laboratory biochemical reaction data and microbial transport data. Then, the verified model is used to investigate the impacts of iron-reduction bacteria (IRB) activities on UHS in aquifers. Based on the simulation results, it can be concluded that (1) hydrogen saturation at the top surface of the aquifer is the greatest while microbial activities surrounding the injection well is stimulated the most; (2) microbial activities influence the initial few cycles of hydrogen injection and withdrawal but the impacts gradually diminish with the dissolution of Fe2O3; (3) hydrogen recovery efficiency is degraded due to the combined effects of hydrogen consumption, water production and microbial clogging with the microbial clogging impact being the most significant; and (4) effective stress impacts aquifer permeability throughout UHS operations while microbial clogging influences it in the initial few cycles. For mineral dissolution/precipitation, the impact can be neglected.

New Criteria to Estimate Local Thermal Nonequilibrium Conditions for Heat Transport in Porous Aquifers

AbstractA fundamental assumption in numerous studies of heat transfer in porous media is local thermal equilibrium (LTE), which assumes that the temperature of the porous media at the fluid and solid interface is in instantaneous equilibrium. Although significant efforts have been made to quantify the occurrence and consequences of local thermal nonequilibrium (LTNE), where the temperatures of the fluid and adjacent solid phases differ, there is no simple expression for quantifying the occurrence and effects of local thermal disequilibrium. Using a numerical model combining LTE and LTNE models, we develop here two simple general criteria based on Darcian velocities (q) and particle sizes (dp) of porous media for determining when LTNE effects occur (denoted as g(dp, q)) and when they become significant (denoted as f(dp, q)). Results show that using an LTE model can result in an underestimation of effective thermal diffusivity and the unaffected Darcian velocities when g(dp, q) > 0. It is possible that using the LTE model can result in an underestimation of the effective thermal diffusivity by more than 200 times within Darcian velocities ranging from 0 to 60 m/d. In the case of g(dp, q) < 0, the use of the LTE model can result in an overestimation of effective thermal diffusivity and Darcian velocities. The performances of the newly developed general criteria are demonstrated using three typical data sets and corresponding numerical models. These data sets include new heat tracer tests conducted in the laboratory and the field, as well as temperature‐time series collected in streambed sediments from a previous study by Shanafield et al. (2012, https://doi.org/10.5194/hessd‐9‐4305‐2012). The potential LTNE effects should be considered when using heat as a tracer to characterize flow and heat transport in porous media in the presence of Darcian velocities less than 2 m/d and particle sizes larger than 10 mm.

Sedimentary and hydraulic characterization to identify hydrofacies in Barreiras Formation, Rio de Janeiro state - Brazil

The use of analogous outcrops in characterizing hydrogeological heterogeneities is widespread worldwide. It is a powerful tool for sediment characterization and determination of hydraulic parameters in rural and urban areas lacking hydrogeological information. This study investigates and characterizes the depositional and post-depositional factors that control the quality of the porous aquifer of the Barreiras Formation, Campos Basin, Rio de Janeiro, Brazil. The integrated analysis between petrological and hydraulic parameters in two outcrops analogous to the Barreiras aquifer demonstrated the role of post-depositional processes, which altered the original deposition conditions. The two outcrops are predominantly composed of sandy and muddy facies, with hydraulic conductivity ranging between 10−4 to 10−5 cm/s and 10−5 to 10−8 cm/s, respectively. This allowed the definition of two hydrofacies, formed by sandy layers as reservoirs and by muddy layers acting as hydraulic barriers. In this context, the Barreiras Formation can be characterized as a poor aquifer, differing from the typical braided fluvial deposits. This evidence draws attention to the indiscriminate use of typical sedimentary models in environments with a high rate of chemical-biological weathering, not only in underground flow models but also in studies of underground resources protection and exploitation, such as the Barreiras aquifer, which consists of one of the main water sources for public, industrial and agricultural supply.

A Novel Method to Calculate Water Influx Parameters and Geologic Reserves for Fractured-Vuggy Reservoirs with Bottom/Edge Water

In practical oilfield production, the phenomenon of water influx typically shortens the water-free recovery period of wells, leading to water flooding and causing a sharp decline in the production well yields, bringing great harm to production. Water invasion usually occurs as a result of the elastic expansion of the water as well as the compaction of the aquifer pore space. However, it can be due to the special characteristics of fractured-vuggy reservoirs such as non-homogeneity and the discrete distribution of the pore spaces. It is challenging to use traditional seepage flow theories to analyze the characteristics of water influx. Also, reservoir numerical simulation methods require numerous parameters which are difficult to obtain, which significantly reduces the accuracy of the results. In this study, considering the driving energy for water influx, a water influx characteristic model was obtained by fitting a graph plate. Subsequently, an iterative calculation method was used to simultaneously obtain water influx volume and OOIP. The aquifer to hydrocarbon ratio was determined by fitting the water influx curve with the graphic plate. Results show that the calculation method is sensitive to the values of reservoir pressure and the crude oil formation volume factor. After applying the method to one field case, it was discovered that water influx performance can be characterized into two types, i.e., linear water influx and logarithmic water influx. In the early stages, the water influx rate of logarithmic water influx is greater compared to linear water influx. However, the volume and energy of waterbody are limited, and the water invasion phenomenon occurs almost exclusively within a short period after the invasion. On the other hand, the volume of waterbody invaded by linear water influx is larger, and it can maintain a stable rate of water influx. The results of the study can provide theoretical support for the waterbody energy evaluation and dynamic analysis of water influx, as well as the control and management of water in these types of reservoirs.

Groundwater Depletion. Are Environmentally Friendly Energy Recharge Dams a Solution?

Groundwater is a primary source of drinking water; however, groundwater depletion constitutes a common phenomenon worldwide. The present research aims to quantify groundwater depletion in three aquifers in Greece, including the porous aquifers in the Eastern Thermaikos Gulf, Mouriki, and the Marathonas basin. The hypothesis is to reverse the phenomenon by adopting an environmentally acceptable methodology. The core of the suggested methodology was the simulation of groundwater using MODFLOW-NWT and the application of managed aquifer recharge (MAR) by using water from small dams after the generation of hydropower. Surface run-off and groundwater recharge values were obtained from the ArcSWAT simulation. The predicted future climatic data were obtained from the Coordinated Regional Climate Downscaling Experiment (CORDEX), considering the Representative Concentration Pathway (RCP) 4.5 and the climate model REMO2009. Groundwater flow simulations from 2010 to 2020 determined the existing status of the aquifers. The simulation was extended to the year 2030 to forecast the groundwater regime. In all three sites, groundwater depletion occurred in 2020, while the phenomenon will be exacerbated in 2030, as depicted in the GIS maps. During 2020, the depletion zones extended 11%, 28%, and 23% of the aquifers in Mouriki, the Eastern Thermaikos Gulf, and the Marathonas basin, respectively. During 2030, the depletion zones will increase to 50%, 42%, and 44% of the aquifers in Mouriki, the Eastern Thermaikos Gulf, and the Marathonas basin, respectively. The simulation was extended to 2040 by applying MAR with the water from the existing dams as well as from additional dams. In all sites, the application of MAR contributed to the reversal of groundwater depletion, with a significant amount of hydropower generated. Until 2040, the application of MAR will reduce the depletion zones to 0.5%, 9%, and 12% of the aquifers in Mouriki, the Eastern Thermaikos Gulf, and the Marathonas basin, respectively. Apart from over-pumping, climatic factors such as long periods of drought have exacerbated groundwater depletion. The transformation of dams to mini-scale hydropower facilities combined with MAR will benefit clean energy production, save CO2 emissions, and lead to an economically feasible strategy against groundwater depletion.

Unveiling subsurface heterogeneity in porous aquifers: Insights from hydrogeophysics and derivative analysis

Groundwater is a crucial resource, particularly in urban areas where the need for alternative water sources is rising. Securing groundwater resources is paramount, requiring a concentrated effort to assess these resources and mitigate the increasing water shortages in urban regions. Despite its significance, the application of hydrogeophysics and derivative analysis in understanding aquifer dynamics is often overlooked and underused. Consequently, this study argues that neglecting the integration of hydrogeophysics dataset and derivative analysis in aquifer characterization leads to developing models that lack solution-oriented approaches, hindering effective groundwater management. This study aimed to enhance understanding of aquifer dynamics for solution-based modeling while emphasizing the importance of integrating hydrogeophysics dataset and derivative analysis to highlight aquifer system heterogeneities. The electrical resistivity tomography and derivative analysis of pumping test were used in this study to image the subsurface and amplify aquifer flow regimes respectively. The results of the electrical resistivity survey revealed a distinct layer of fine to medium-grain sand at depths of approximately 60 m in certain areas intercalated with medium-to-coarse-grain sand and thin layers of peat. Derivative analysis plots indicated that the predominant flow regime is linear and bilinear, with evidence of fracture dewatering during pumping cycles suggesting an unconfined aquifer. Additionally, aquifer heterogeneity and a no-flow boundary were evident during the pumping cycle. This study underscores the efficacy of combining hydrogeophysics and derivative analysis of pumping tests as a robust approach for imaging and validating subsurface conditions. The implications of these findings extend beyond the specific case study, offering valuable insights for groundwater utilization, monitoring, and management on a broader scale. By elucidating effective methodologies, this research enhances groundwater security and meets the escalating demand for sustainable water sources in urban areas. In conclusion, electrical resistivity tomography is an effective tool for characterizing subsurface aquifer systems, while derivative analysis of pumping tests is crucial for highlighting aquifer system heterogeneities. These methods offer a more practical interpretation compared to traditional drawdown versus time curve approaches, making the results invaluable for groundwater usage monitoring and management.

Mixing and transport of CO2 across a monolayer-covered surface in an open cylinder driven by a rotating knife edge

The dissolution of CO2 into water has attracted much attention for the modeling of air-sea gas exchange, CO2 sequestration into porous aquifers, and gas transport in bioreactors. In nearly all circumstances, the surface of water is covered by surfactants that make the gas-liquid interface behave elastically and exhibit surface (excess) shear viscosity. Here, transport and mixing of CO2 across the interface is studied in an axisymmetric system with inertia. A rotating knife edge at the interface with negligible contact surface area takes advantage of surface shear viscosity to drive a strong overturning flow in the bulk via viscous coupling. When CO2 is dissolved into water, this system has a forced convection nature that enhances the transport of CO2 into water. This transport leads to intense radial gradients in the dissolved CO2 resulting in baroclinic production of azimuthal vorticity which alters the underlying meridional flow and further enhances mixing by breaking down transport barriers. Flow visualization experiments reveal baroclinic instability and numerical simulations detail the hydrodynamics, including the effects of varying the CO2 partial pressure above the water.

Assessment of contamination dispersion in a porous aquifer environment using explicit and implicit methods in the MODFLOW–MODPATH model

This study presents a meticulously edited article that assesses the uncertainty associated with contamination dispersion in a porous aquifer environment. The article employs a conceptual model incorporating finite difference flow modeling and particle tracking methods grounded in particle movement theory within groundwater flow. The groundwater flow model utilizes an automatic calibration technique, PES, and independent verification, to minimize statistical discrepancies in the optimized parameters. The statistical summary of the model output reveals the average pollution concentration in the study area, specifically the landfill aquifer in northern Iran. The results indicate that the initial contamination value at the boundaries is zero, whereas 91,802 units are discharged into the sea. Moreover, approximately 3,317,290,003 pollution units have been extracted from pumping wells along the contamination plume over 10 years. Furthermore, the drainage network releases 1,778,680,001 pollution units, resulting in a net outflow of 5,082,740,007 units from the same drainage network. Additionally, 9,912,180,003 units were discharged from the leaky boundaries. It should be noted that these statistics can be reversed depending on the specified period. The difference between the input and output pollution values in the study area indicates a net accumulation of 8192 units. To analyze the sensitivity of contamination dispersion in the porous aquifer environment and the variations in coastal drainage patterns obtained from the finite difference simulation, the explicit and implicit methods, as well as the random walk particle tracking method, were investigated in a portion of the basin area using the South Side Implicit method. One of the main objectives of this sensitivity analysis is to evaluate the results of the uncalibrated qualitative model and examine the particle dispersion in the adjacent area of the contaminated landfill region to assess the error of the proposed model by developing a limited test. The calculated error or computational difference between the explicit and implicit methods and the applet model is considered negligible. For example, the calculated error between the explicit and backward methods and the applet model ranges from 0.086 to 0.372 units for a concentration interval of 0.512 units and 0.680 units for a time interval 0.537. Therefore, the computational differences are not significant. Furthermore, based on theoretical sensitivity analysis of contamination dispersion, the laboratory model demonstrates that the flow direction in this modeling is from west to east, as expected.

Hydrodynamic Characterization of Carbonate Aquifers Using Atypical Pumping Tests without the Interruption of the Drinking Water Supply

The Gran Sasso carbonate aquifer is the largest and most productive in the Apennines. Its hydrogeological structure has been studied since the middle of the last century for the springs’ characterization for drinking purposes and for a motorway tunnel. Meanwhile, its hydrodynamic parametrization is less developed and has been limited to monitoring the discharge and chemical and isotopic parameters. Secondary porosity characterizes the aquifer, and an underlying impermeable marly complex represents the basal aquiclude. It might appear inappropriate to characterize the hydraulic properties via pumping tests, as their reliability has been proven in homogeneous and isotropic media. However, the high extent of the aquifer, the wells’ location, the scarcity of information available and the lack of alternatives has forced the estimation of hydrodynamic parameters as in porous aquifers and the experimental testing of the aquifer, especially in maximum pumping conditions, for a possible exploitation increase. Since aquifer testing was performed during the normal well field’s activities, it was not possible to perform typical tests. Therefore, the step-drawdown test was conducted by turning on an increasing number of wells over time and keeping the observation points fixed. As results, a mean hydraulic conductivity of 5 × 10−3 m/s and a mean transmissivity of 0.3 m2/s were established without interrupting the water supply; meanwhile, the influence radius and flow directions were also estimated.

Study on Migration Characteristics of Pollutants in Groundwater at a Proposed Hazardous Waste Landfill

Objective: The paper aims to analyze the hydrogeological conditions of a proposed hazardous waste landfill and the migration characteristics of lead, zinc, and nickel in fractured aquifers and porous aquifers under accident conditions and provide a reference for the influence of the proposed landfill on groundwater. Methods: In this study, based on a 1:50000 regional hydrogeological survey and 1:2000 site hydrogeological mapping, the hydrogeological conceptual model was established. Finite difference software GMS was used to analyze the migration characteristics. Results: The study demonstrated that when the pollutants in the hazardous waste landfill leaked, they migrated from northeast to southwest along the gully. The pollutants in the porous aquifer migrated quickly, and the polluted area expanded rapidly from point to surface. The pollutants migration in fractured aquifers was slow, and the groundwater quality was deteriorating continuously. During the simulation period, the pollutants of lead, zinc and nickel all polluted the aquifer. Among them, the lead pollution range w reported to be the largest, with an exceeding distance of 216.7 m; the zinc pollution range was the smallest, with an exceeding distance of 33.3 m, and the exceeding distance of nickel was 165.1 m. Conclusion: In order to ensure the safety of the groundwater environment in the simulated area, the impervious treatment must be carried out according to the requirements of the proposed hazardous waste landfill. Meantime, an emergency plan should be formulated.

Assessing the salinization mechanisms of coastal brackish springs

Brackish springs have been extensively reported in Mediterranean coastal carbonate formations. The phenomenon is puzzling because these springs may discharge high flow rates with significant salinities at elevations of several meters above sea level. Although these springs have been studied for millennia (since the ancient Greeks), controversy persists. In essence, they are typically assumed to consist of a karstic conduit bringing fresh groundwater that meets a saltwater conduit at a branching point, where the brackish mixture flows up to the spring mouth. Here, we analyze the hydraulics of the system for two simple cases, depending on whether seawater flow occurs through an open conduit or a porous aquifer. To this end, we derive the equations governing variable-density turbulent flow through the conduits. These turn out to be similar to the traditional ones, except that (1) Bernouilli’s head (energy per unit weight) is substituted by the energy per unit volume, and (2) this energy is not a potential, as a flow path dependent term needs to be added. We solve these equations to obtain the characteristic functions relating fresh to seawater discharge depending on spring concentration. These functions are specific to every brackish spring and representative of the karst system and salinization mechanism. They allow us to show that the spring salinity is mostly controlled by the weight of the water column flowing towards the spring mouth (for low flow rates) and by energy dissipation (for high flow rates). When concentration and flow rate data are available, it is possible to characterize: (1) the resistance of the upwards conduit from the response at high flow rate; (2) the depth of the branching point from the concentration at low discharges; and (3) the hydraulic resistance of the seawater conduit from the slope of concentration vs. flow rate at low spring discharges. These patterns are compared to field data from three different springs, which yields insights on the conceptual model governing every particular spring.

Anaerobic dihydrogen consumption of nutrient-limited aquifer sediment microbial communities examined by stable isotope analysis

ABSTRACT The biogeochemical consequences of dihydrogen (H2) underground storage in porous aquifers are poorly understood. Here, the effects of nutrient limitations on anaerobic H2 oxidation of an aquifer microbial community in sediment microcosms were determined in order to evaluate possible responses to high H2 partial pressures. Hydrogen isotope analyses of H2 yielded isotope depletion in all biotic setups indicating microbial H2 consumption. Carbon isotope analyses of carbon dioxide (CO2) showed isotope enrichment in all H2-supplemented biotic setups indicating H2-dependent consumption of CO2 by methanogens or homoacetogens. Homoacetogenesis was indicated by the detection of acetate and formate. Consumption of CO2 and H2 varied along the differently nutrient-amended setups, as did the onset of methane production. Plotting carbon against hydrogen isotope signatures of CH4 indicated that CH4 was produced hydrogenotrophically and fermentatively. The putative hydrogenotrophic Methanobacterium sp. was the dominant methanogen. Most abundant phylotypes belonged to typical ferric iron reducers, indicating that besides CO2, Fe(III) was an important electron acceptor. In summary, our study provides evidence for the adaptability of subsurface microbial communities under different nutrient-deficient conditions to elevated H2 partial pressures.

Uranium anomaly in groundwater of the hard rock aquifer system in southeast Brazil

This study presents a hydrogeological characterization that fills knowledge gaps in an area where a complex geological scenario connects two fractured (granite/metamorphic rocks and diabase) and one porous (sedimentary rocks) aquifer systems. Previous groundwater studies from the four local deep supply wells suggested the sedimentary aquifer as the source of uranium (∼60 μg/L) in one of the wells. Despite the water extraction interruption from this well, this finding needed to be supported by evidence provided by the geological and structural characteristics of the aquifers. Data collected through video inspection, caliper, optical (OPTV), and acoustic borehole televiewers (ATV) revealed the predominance of the crystalline aquifer in the area. As recorded by video inspection, the well's borehole with high U intersects the crystalline rocks at several open fractures and presents the highest gamma radiation values. Therefore, we infer that the origin of dissolved U is the weathering of primary minerals locally occurring in those rocks. The U leaching and complexation as calcium-uranyl-carbonates favors its solubility in that specific aqueous medium, as shown by hydrogeochemical modeling using PHREEQC. Additionally, the major ions and trace elements concentrations allow the division of the wells into two groups: a) two with higher U concentrations, influenced by the granite/metamorphic aquifer and possibly longer residence times, and b) the other two wells with lower U concentrations, higher nitrate concentrations, influenced by the diabase aquifer, possibly have younger residence times.

Estimating the vulnerability of groundwater resources to diffuse pollution in highlands areas: review of the literature and critical analysis (highlands of Cameroon)

Abstract Groundwater is a major resource for drinking water, especially in developing countries, where it is less expensive to treat than surface water. Today, the resource is highly susceptible to pollution, particularly as a result of human activity. This review was based on a literature review and critical analysis of models for estimating the groundwater vulnerability. The results show that the deepest porous aquifers are the least susceptible to pollution, whereas those in karstic and fissured environments are susceptible, whatever their depth. Pollution usually arises from human activity. Critical analysis of the literature shows that existing methods are developed in specific environmental contexts. Given the variability of factors in space and time, these methods do not take the intrinsic realities of all natural settings into account adequately and are not perfectly applicable in all environments. This highlights the need to develop appropriate models for each environment, such as that of the highlands in countries such as Cameroon.

Evaluating the impact of muon-induced cosmogenic 39Ar and 37Ar underground production on groundwater dating with field observations and numerical modeling

Groundwater dating by radioactive cosmogenic tracers such as 39Ar relies on the decay rate from a known initial atmospheric activity (100%modern). Thereby, it is assumed that cosmogenic 39Ar production in the subsurface is negligible at depths below the water table and that contributions from natural rock radioactivity are minor or missing. Here we present 39Ar data from aquifers located in quaternary glacial sediments and tertiary limestones in Denmark, which unequivocally demonstrate that cosmogenic production can induce considerable age biases. 39Ar values larger than 100%modern are observed at relatively shallow groundwater depths in non-radiogenic rocks. These activities are compared to calculations based on previously assessed depth-dependent production rates in rocks and realistic estimates of the emanated fractions to the water phase. The water residence time distribution with depth, which was determined by numerical flow modeling and particle tracking, underpinned the significance of muon-induced 39Ar production. The short-lived isotope 37Ar is produced by similar processes as 39Ar and demonstrated its usefulness as an indicator of local underground production in an aquifer. The significance of cosmogenic underground production in other possible recharge scenarios was then assessed by explicitly simulating the radioargon accumulation and decay in a 2D synthetical numerical model. These simulations demonstrated that underground production is negligible when the water infiltrates freely in a porous aquifer. However, in the presence of a confining layer impeding the infiltration at shallow depths (<30 m), as is the case in our study site in Denmark for instance, over-modern 39Ar activities (>100%modern) may occur. The age concluded from the dissolved activities is then possibly biased towards young values. Special attention should thus be paid to the recharge rates when using 39Ar for dating groundwater. 37Ar activities provide complementary information about the strength and mechanisms of underground production.

Understanding groundwater storage and drainage dynamics of a high mountain catchment with complex geology using a semi-distributed process-based modelling approach

Mountainous groundwater systems are usually characterized by great geological complexity and highly heterogeneous hydraulic properties. For that reason, proper system characterization, monitoring and modelling remain challenging. In this study, we investigated a geologically complex alpine catchment (from 1400 to 2900 m asl) in the Dolomites (Italian Alps) by combining hydrogeological field investigation, hydrological monitoring and modelling. A semi-distributed process-based hydrological model was applied to simulate the continuous measured catchment discharge in a period of three years, which covers a large variation of hydrodynamic conditions. The model structure couples the sequential hydrogeological units within the studied catchment: (1) the fractured dolomitic rocks as bedrock aquifer and 2) the unconsolidated deposits accumulating on the slopes and at the valley floor as porous aquifer. In order to evaluate the model structure and parameterization in depth, we applied a multi-step evaluation approach considering both parameter sensitivity and uncertainty. The modelling results demonstrate that the newly developed model can reproduce most discharge behavior of aquifers. The model indicates a seasonal dynamic linkage between surface and subsurface storage units during different flow conditions. Besides the matrix and conduit flow in fractured dolomitic aquifer, it highlights the important role of unconsolidated sediments (porous aquifer) to the storage and discharge behavior of the entire groundwater system. Furthermore, with the comprehensive model evaluation we learned the model performance deficit during extreme high flow conditions and proposed a more detailed hydrogeological conceptual model. We believe that the proposed modelling approach can be transferred to other high mountain catchments with similar hydro(geo)logical characteristics to characterize the collective behavior of aquifer systems, explore conceptual model weakness and optimize catchment monitoring work.

Water Level in Observation Wells Simulated From Fracture and Matrix Water Heads Outputted by Dual‐Continuum Hydrogeological Models: POWeR‐FADS

AbstractDo observation wells in fractured porous aquifers measure water head in the fracture network, water head in the matrix, or some combination of both? This question necessarily arises when calibrating dual‐continuum hydrogeological models against on‐field data. One can assume that observation wells measure fracture water head, because matrix permeability is negligible compared to fracture permeability. Nevertheless, this reasoning is invalid for wells poorly‐connected to the fractures. Yet, the possibility of such a poor connection at given depths has never been implemented in a physics‐based manner when comparing matrix and fracture water heads simulated by dual‐continuum models to on‐field data. To fill this knowledge gap, a physically based, easy to calibrate, open‐source postprocessing tool, POWeR‐FADS (Program for Observation Well Representation in Fractured Aquifer Dual‐continuum Simulations), available at https://github.com/BJeannot1/POWeR-FADS, has been developed. It introduces as parameters well geometry and the altitude of lowest interception of the fractures by the well. From these, POWeR‐FADS nonintrusively postprocesses time series of matrix and fracture water heads at the well, as simulated by any planar, bidimensional dual‐continuum hydrogeological model, to calculate water exchanges involving the observation well and thus the evolution of water level in the well. Synthetic test cases show that POWeR‐FADS makes it possible to simulate peculiar behaviors that are similar to patterns actually observed by the authors in on‐site observation wells of a fractured porous aquifer, like “floors” in observed water levels, delayed but sharp rises at the beginning of recharge events, or inflexion points accelerating the drawdown velocity during the recession phase.