Groundwater Model Research Articles

A novel groundwater monitoring network design framework for long-term and economical data monitoring

The groundwater monitoring network (GMN) design balancing prediction accuracy with cost minimization has paramount importance in depth to groundwater (DTW) monitoring, modeling, and management. This study demonstrates a robust and transferable GMN design framework for the underexplored statistically homogenous DTW to facilitate appropriate, long-term, cost-effective, and regional-scale groundwater monitoring. The inherent homogeneity in regional scale DTW data was investigated to cluster the observation wells into 18 strongly structured clusters/strata using optimal cluster number (k = 18), corresponding to the highest mean silhouette score of 0.72 and the K-means clustering algorithm. The clustered observation wells were used to evaluate the cluster random sampling (CRS) and stratified random sampling (SRS) for GMN design. The sampling techniques were evaluated to capture the spatial variability of DTW over the whole study region through spatial modeling at an adequate accuracy by the inverse distance weighting interpolation tools. The performance metrics for SRS recommended it as a more robust and appropriate technique for regional-scale GMN design than CRS. Then, the GMN was designed with 540 wells corresponding to 70% sampling under SRS for accurate, cost-effective, and long-term groundwater monitoring. By reducing 30% of observation wells, the cost-saving can contribute to cutting-edge monitoring infrastructure development. Additionally, the need-based balancing of the accuracy in spatial modeling and monitoring cost may further moderate the number of observation wells using the performance metrics of SRS corresponding to sampling percentages. Conclusively, this novel and transferable framework can be adopted globally with internally homogenous DTW data for effective GMN design.

Drivers of barrier island water-table fluctuations and groundwater salinization

Barrier islands are threatened by climate change as sea-level rise and higher frequency storm surge lead to more flooding and saltwater intrusion. Vegetation plays a vital role in preventing erosion of barrier islands due to aeolian and hydrological forces. However, vegetation on barrier islands is threatened by rising water tables causing hypoxic conditions and storm-surge overwash introducing saline water to the root zone. To better protect barrier island ecosystems, it is critical to identify the relative influence of different hydrological drivers on water table elevation and salinity, and understand how this influence varies spatially and temporally. In this study, three barrier island sites were instrumented with groundwater wells monitoring water level and specific conductance. Using these data, a set of transfer function noise models were calibrated and used to determine the relative influence of hydrologic drivers including precipitation, evapotranspiration, bay and ocean water levels, and wave height on groundwater levels and specific conductance. We found that drivers of water-level change and specific conductance vary strongly among sites, depending primarily on the surface water connectivity and the geology of the island. Sites with close connection to inlets showed more salinization and responded to a larger number of drivers, while sites that were poorly connected to the ocean responded to fewer drivers.

Quantifying Seepage‐Face Evaporation and Its Effects on Groundwater Flow and Solute Transport in Small‐Slope Tidal Flat

AbstractLarge‐scale seepage faces occur on tidal flats with gentle slope, which are widely distributed worldwide. Evaporation on these seepage faces, leading to salt retention and accumulation, may significantly impact the density‐dependent groundwater flow beneath the tidal flats. However, due to nonlinear complexities of the groundwater flow and solute transport on seepage faces, explicit boundary conditions and numerical models to quantify these processes are lacking. In this study, we present both mathematical and numerical models to quantify these processes. Compared to the results of our previous study, this paper shows that seepage‐face evaporation can (a) significantly increase the groundwater salinity in the upper intertidal zone, and form multiple groundwater circulation cells in the intertidal zone, (b) cause the disappearance of multiple seepage‐faces and reduce the spatial extent of seepage faces notably, (c) and intensify the groundwater and salt exchange as well as the seawater‐groundwater circulation through the intertidal zone.

Scrutinizing different predictive modeling validation methodologies and data-partitioning strategies: new insights using groundwater modeling case study

Groundwater salinity is a critical factor affecting water quality and ecosystem health, with implications for various sectors including agriculture, industry, and public health. Hence, the reliability and accuracy of groundwater salinity predictive models are paramount for effective decision-making in managing groundwater resources. This pioneering study presents the validation of a predictive model aimed at forecasting groundwater salinity levels using three different validation methods and various data partitioning strategies. This study tests three different data validation methodologies with different data-partitioning strategies while developing a group method of data handling (GMDH)-based model for predicting groundwater salinity concentrations in a coastal aquifer system. The three different methods are the hold-out strategy (last and random selection), k-fold cross-validation, and the leave-one-out method. In addition, various combinations of data-partitioning strategies are also used while using these three validation methodologies. The prediction model’s validation results are assessed using various statistical indices such as root mean square error (RMSE), means squared error (MSE), and coefficient of determination (R2). The results indicate that for monitoring wells 1, 2, and 3, the hold-out (random) with 40% data partitioning strategy gave the most accurate predictive model in terms of RMSE statistical indices. Also, the results suggested that the GMDH-based models behave differently with different validation methodologies and data-partitioning strategies giving better salinity predictive capabilities. In general, the results justify that various model validation methodologies and data-partitioning strategies yield different results due to their inherent differences in how they partition the data, assess model performance, and handle sources of bias and variance. Therefore, it is important to use them in conjunction to obtain a comprehensive understanding of the groundwater salinity prediction model's behavior and performance.

Cation exchange and leakage as dominant processes in controlling salinity and strontium in sandy and argillaceous coastal aquifer: Insights from hydrochemistry and multi-isotopes

Coastal groundwater is crucial for supporting ecosystems and meeting the water needs of coastal communities. However, over-exploitation of groundwater resources, including fresh water and deep brine, leads to groundwater depletion and exacerbates saline intrusion, posing significant challenges to sustainable development. This study focuses on the south bank of Laizhou Bay, a typical area threatened by saline intrusion due to the salt industry’s extraction of deep brine and over-exploitation of fresh groundwater. To assess the impact on salinization processes, we employed hydrochemistry and multi-isotope techniques to study strontium (Sr) distribution and salinization in coastal groundwater, unraveling their origins and explaining hydrochemical processes. The hydrochemical analysis revealed that evaporite dissolution primarily contributes to salinity, especially in deep brine and saline water (TDS > 30 g/L). The relationship between 87Sr/86Sr and chlor-alkali index (CAI) and saturation index (SI) indicated distinct sources of Sr. In fresh and brackish water (<10 g/L), Sr originates from mineral weathering and dissolution (SIs of albite and calcite increase with salinity), while in deep brine and saline water, Sr is influenced by cation exchange (with CAI-Ⅰ and CAI-Ⅱ greater than 0). Additionally, the interplay of anthropogenic input (nitrate), Sr, and EC highlights the differential anthropogenic impact on shallow groundwater, emphasizing the relatively closed environment of deep brine. Multi-isotope analysis showed that the average 87Sr/86Sr of shallow groundwater samples in the cone depression is 0.7112, between that of upstream surface water and shallow groundwater (0.7119) and deep brine and saline water (0.7101), indicating groundwater extraction induces deep brine leakage into shallow groundwater, elevating salinity and Sr concentrations. This poses a direct threat to local water security, including irrigation suitability and health risks of drinking water.In conclusion, we categorized the relationships between 87Sr/86Sr and Sr under different hydrochemical processes into different end-members based on their origins. This offers valuable insights in comparable regions for identifying groundwater salinization and saline intrusion, and the use of Sr and its isotopes can also help trace cation exchange and deep brine leakage in coastal aquifers.

Effects of hydraulic conductivity on simulating groundwater–land surface interactions over a typical endorheic river basin

Hydraulic conductivity is the key variable in subsurface configuration that is important in modeling groundwater and its interactions with land surface processes. However, an acquisition of quantitatively distributed hydraulic conductivity is challenging due to limitations on data availability, accuracy and resolution, which inhibits an advanced understanding of hydrologic and energy cycles at the basin scale. This study investigated the role of hydraulic conductivity in simulating groundwater-land surface interactions over a typical endorheic river basin in northwestern China. A 3D variably saturated groundwater model coupled to a land surface model was configured to account for interactions between groundwater dynamics and land energy fluxes. Four scenarios of hydraulic conductivity were tested and inter-compared, adopting both homogeneous and heterogeneous parameterization strategies based on Global Soil Dataset for use in Earth System Models (GSDE) and HYdrogeology MaPS (GLHYMPS). The overall objectives are to diagnose effects of hydraulic conductivity on simulating groundwater–land surface interactions, with a particular focus on differences in spatial distributions and statistical characteristics between scenarios, as well as the dynamics and feedback mechanisms between groundwater and land surface energy fluxes. Comprehensive simulations and analyses were conducted to discern spatial distributions and statistical characteristics of key hydrological variables (streamflow, water table depth, latent heat flux, sensible heat flux) under different scenarios. Results showed that the hydraulic conductivity of the saturated zone had a significant control over water table depth, while the combine hydraulic conductivity in the vadose zone and the saturated zone governed the evolution of streamflow. Moreover, the water table critical range—where the land energy budget is highly sensitive to groundwater table—was 0.1 m–5 m and the area might expand approximately 4.5 % after considering the heterogeneity of hydraulic conductivity. Compared to the homogeneous scenario, the spatial heterogeneity of hydraulic conductivity amplified water table depth variability while diminishing land energy fluxes variability.

Feasibility of sustainability separation of drinking water network by groundwater desalination Bajestan urban catchment, Iran

Fresh water supply in arid land is facing many challenges that require proper management of water resources, especially in drinking water. Management of water projects depends on the conditions of each country or region. Bajestan is located in Razavi Khorasan province, Iran. The certain resource of drinking water supply is desalinated of high salinity groundwater from Playa. Government has implemented the project “Separation of Drinking and sanitary water networks by groundwater desalination” since of 2014, by creating “Government Drinking Water Withdrawal Smart Stations (GDWWSS)" in this city. This study evaluates the feasibility of different methods of separation drinking water from sanitary water by considering sustainability dimensions these methods. These methods including of Home Water Desalination Machines, Mobile Packaged Water, GDWWSS (with new design), Private Drinking Water Withdrawal Stations (PDWWS) and Dual Drinking Water Distribution Network (DDWDN). Findings showed, that if just the economic dimension is considered in the drinking water separation project, PDWWS would be the best method but since the social dimension is an integral part of a sustainable project, the most compatible method is DDWDN. The environmental evaluation showed that production of salt water and salt (NaCl) is the most common adverse effect of all the proposed methods, but the DDWDN will produce more salt wastewater. On the other hand, Mobile packaged water method will produce 17,000 tons of PET annually, which will cause biological problems if it is not recycled. Therefore, the best way is to use the water resources of the plain and not use the desalination of salt groundwater resources. Also, to increase people's satisfaction and reduce water consumption with good quality water, DDWDN is suggested.

Managed aquifer recharge in island aquifer under thermal influences on the fresh-saline water interface

Storing and utilizing freshwater in coastal aquifers with managed aquifer recharge (MAR) presents significant challenges, particularly in small island settings where seawater intrusion is a major source of groundwater salinization. Existing studies usually evaluate density-dependent flow under isothermal condition, which neglects the potential impact of temperature on fluid density. To fill the knowledge gap, we performed both laboratory experiments and numerical modeling to investigate the impact of warmer water injection on freshwater storage and extracted water quality for island aquifer. It is found that warmer water (70 °C) injection compromise the “barrier effect” achieved by normal water injection, resulting in a reduced ability to prevent coastal seawater intrusion by, creating an upward hydraulic gradient beneath the injection well. In addition, warm water injection result in a 5 % decrease in freshwater storage capacity for island aquifer. However, by employing a combination of receding tide cycles and intermittent injection-extraction operations, the duration of lowered saline concentration within the extracted water can be extended by approximately 10 % with an 11 % reduction in the overall extracted saline concentration. This study provides valuable insights into the influence of warmer water injection on coastal aquifer quality. The findings from both the laboratory experiments and numerical simulations enhance the understanding of thermal density-dependent dynamics of coastal aquifers, offering valuable guidance for effective freshwater storage and usage in island settings.

Potential impacts of saline groundwater pumping on seawater intrusion in a coastal aquifer system

ABSTRACT Desalination plants employing feed water sourced from saline groundwater (SGW) are a feasible approach to tackle freshwater scarcity and seawater intrusion (SWI). Available studies lacked comprehensive analysis quantifying the impact of pumping from a coastal aquifer, on changes in solute transport mechanism. This study investigates the impact of different rates of SGW pumping on the solute transport in a coastal aquifer. A regional-scale model resembling the Upper Floridan aquifer was modelled using SEAWAT to quantify SWI variables. The study found that the length of intrusion of toe (Ltoe), width of mixing zone (WMZ) and residual salt mass (RSM) varied with dispersivity values. The ratio of saltwater to freshwater influx was unaffected by dispersivity values. Significant decrease in the toe-length occurred when the pumping rate was highest and the dispersivity was lowest. The WMZ was found to be the most sensitive SWI variable towards varying dispersivity, with 65% difference for higher pumping rates. The mid-bottom region was experiencing substantial modifications, leading to an inclined response at increased pumping rates, which suggested that the diluting process was primarily taking place in the bottom portion. This study gives insight into saline wedge movement in a coastal aquifer system which will benefit in the design of SGW well system for desalination plants.

Delineating the mechanisms controlling groundwater salinization using chemo-isotopic data and meta-heuristic clustering algorithms (case study: Saguenay-Lac-Saint-Jean region in the Canadian Shield, Quebec, Canada).

This study employed meta-heuristic clustering algorithms to determine the source and mechanism of groundwater salinization in Quebec's Saguenay-Lac-Saint-Jean (SLSJ) region, utilizing hydrogeochemical (38 inorganic constituents, including minor, major, and trace elements) and isotopic data (δ18O and δ2H). A total of 382 groundwater and precipitation samples were examined. Among the meta-heuristic algorithms, Artificial Bee Colony K-Means (ABCKM), Differential Evolution K-Means (DEKM), Harmony Search K-Means (HSKM), Particle Swarm Optimization K-Means (PSOKM), and Genetic K-Means (GKM) were used and investigated, and finally, PSOKM displayed superior performance and was chosen for further investigation. Analysis of diverse plots and hydrogeochemical modeling unveiled the impact of the Laflamme Sea invasion on groundwater chemistry. PSOKM1, PSOKM4, and PSOKM5 exhibited notable carbonate and silicate dissolution, with PSOKM4 demonstrating predominant carbonate dissolution. Cation exchange was identified through binary plots and Chloro Alkaline Index (CAI), with reverse cation exchange predominantly observed in most PSOKM4 samples, while positive values suggested direct cation exchange in other clusters. Spatial dynamics analysis using HFE-D indicated that salinization occurs as groundwater flows through crystalline bedrock aquifers, resulting in a transition from HCO3- dominance in PSOKM4 to Cl- dominance in the remaining clusters. Interaction between groundwater and rock along this path facilitated a transformation towards a Na-Cl end-member. The closely aligned stable isotopes with the Global Meteoric Water Line (GMWL) indicated a blend of meteoric water and seawater as the groundwater source.

An integrated modeling approach to simulate human-crop-groundwater interactions in intensively irrigated regions

Global food security is threatened by aquifer depletion in arid and semi-arid regions. Ensuring aquifer sustainability while minimizing farmers' costs requires effective groundwater management policies, which are informed by hydro-economic models. However, the reliability of these models is typically compromised by oversimplifications. This study introduces a hydro-economic modeling framework (ABM-MODFLOW), applied in the U.S. High Plains, overcoming previous models’ shortcomings. It employs an interdisciplinary approach, simulating real-world farming decisions and the resultant impacts on groundwater systems. Model validation indicates satisfactory performance in reproducing historical data and trends, outperforming the standalone groundwater model. Results show that agents with high irrigation capacities and fine soils earn the highest profits, whereas those with low irrigation capacities and coarse soils achieve the lowest profits. This modeling framework provides an advanced approach for simulating interactions between human decisions, crops, and groundwater in heavily stressed aquifers globally, supporting the development of effective groundwater management policies.

Reversibility of seawater intrusion in a coastal aquifer: Insights from long-term field observation and numerical modeling

It is essential to understand basis processes affecting the reversibility of seawater intrusion (SWI) and the related timescale after groundwater pumping stopped for coastal groundwater management. Groundwater salinity variations presented by previous studies showed a fast SWI process and the following seawater retreat (SWR) process indicated by overall salinity decrease in the coast karst aquifer in the Dalian Peninsula, northeast China. Cross-sectional variable-density flow and transport simulations using equivalent porous medium (EPM) model, dual-domain (DDM) model and the EPM with a single discrete feature representative of conduit flow (EPM-DF) are conducted to understand the SWI reversibility related to the characteristics of the flow system, including aquifer parameters and seaside boundary conditions. Large dispersity combined with a low salinity concentration at the sea boundary are necessary in the EPM and DDM models to produce a fast SWI and low salinity concentration as observed. The enhanced dispersion in the DDM model enhances the mixing in the transition zone (TZ) and produces more obvious seasonal variation but longer SWR process compared to the EPM model. The EPM-DF model produces significant salinity seasonal variation during the SWI process and a fast decrease in the SWR process. The overall fast system response but low peak concentration during the SWI might be attributed to the highly heterogeneous characteristics of the karst aquifer system, which produces a wide front edge of TZ as observed. The dispersion-dominated EPM and DDM models produced higher degree of SWI reversibility compared to the EPM-DF model, as indicated by the weaker lag effects between groundwater levels and salinity concentration. Moreover, the SWI reversibility may arise not only from the heterogeneous characteristics of the karstic aquifer system but also from the rate for rise in the inland groundwater level during the seasonal variation cycle. An overall SWR process can be divided into two stages: the first stage with rapid decrease in salinity concentration but small changes in toe associated with TZ widening and the second stage with rapid toe retreat during aquifer flushing. The flushing of the brine leakage from the solar salt field may extend the SWR process. Considering this relatively longer water quality response time than the time span in general groundwater management plans, a reasonable objective and how maintain the hydraulic head during the SWI controlling have been of a great concern to the coastal groundwater management.

Assessment and mapping of soil salinity and groundwater quality in Nijhum Island, Hatiya Upazila, Noakhali

Crop yield is heavily influenced by soil salinity. The objective of the research is to assess the depth wise soil salinity distribution pattern and groundwater quality by measuring physico-chemical parameters in the agricultural land of Nijhum Island, Noakhali, Bangladesh. This study also explicates the spatial variation of soil salinity in the study area. Soil samples were collected from three different depths to estimate the Electrical Conductivity (EC) using a paste of 1:5 soil and distilled water suspension (1:5 weight-to-volume method). Spatial analyst tool of ArcGIS was used and observed that the land of Nijhum Island was affected with severe salinity in almost 85% of sampling plots and 15% of sampling plots were extremely saline. The soil salinity process has been noticed due to seawater intrusion from the Bay of Bengal through the Meghna estuary which threatens crop production and soil fertility, threatening the population's livelihood. This study provides a baseline understanding of soil salinity in Nijhum Island to help decision-makers and smallholder farmers improve their livelihoods. Thus, groundwater depth did not affect soil salinity significantly. The result contradicts most studies that found a strong influence on soil salinity. Because the study area is coastal, most of the salt comes from seawater inundation and capillary rise, which reduces groundwater availability.

Spatial modelling of present and future groundwater potentials in Nigeria; towards a sustainable water demand and supply

The development of water resources in Nigeria has not risen alongside population increase. Although groundwater is the available source of water at all locations in Nigeria, groundwater aquifer is becoming deeper due to increased groundwater harvest, especially in urban areas. There is a need for continuous groundwater modelling using geological, climatic, environmental data, and spatial tools especially under climate change and intense landcover conversions. The Spatial Multi-Criteria Analysis (SMCA) was adopted in this study to model the current groundwater potential in Nigeria. The multi-criteria analysis tool in ArcGIS was explored to overlay the nine factor maps to model historic and futuristic groundwater potentials. Our groundwater factor maps show an interesting pattern across Nigeria, the southern parts show more potential considering suitability factors like; surface water density, rainfall, temperature, soil, land cover and elevation in the region. Whereas the upland regions even with higher rainfall and lower temperatures suitable for groundwater recharge, are disadvantaged by geomorphological factors. Northern part of Nigeria shows high potentials, considering the geology, soil, lineament density and slope, but disadvantaged by other factors like lower rainfall and higher temperatures. This informed the final groundwater potential maps; results for the historic potential revealed that, no location in Nigeria has optimal (9–10), very poor or no (1–2) groundwater potential. The results further revealed that areas with higher groundwater potentials are largely within the corridors of major rivers in Nigeria (Niger and Benue), covering about 17.6% of the Nigerian landmass, while 2.6%, 33.7, 44% and 2.2% are occupy areas with very low, low, moderate and very high groundwater potentials respectively. Further analysis highlighs locations of concern due to climatic and environmental changes. Interestingly, our groundwater projection results show a persistent increase in groundwater potential from 2021 to 2100 if current landuses and environmental factors are maintained, and if the projected increase in rainfall is true. Despite these groundwater potentials, the recent issues of shallow aquifers have been cautioned by scientists, especially due to groundwater uncertainties in Nigeria especially with intense landcover conversions, combined with accelerated water demands due to increase in population, and incessant groundwater extractions. Therefore, there is a need to seek caution in the pattern of unregulated and incessant groundwater harvest in Nigeria and we recommend frequent updates of the groundwater potential using geospatial tools to inform governing policies on a centralised consolidated sustainable water supply in Nigeria.

Evaluating parameter inversion efficiency in Heterogeneous Groundwater models using Karhunen-Loève expansion: a comparative study of genetic algorithm, ensemble smoother, and MCMC

Evaluating parameter inversion efficiency in Heterogeneous Groundwater models using Karhunen-Loève expansion: a comparative study of genetic algorithm, ensemble smoother, and MCMC

Effects of Biochar and Straw Return on Soil Microbial Community Characteristics and Functional Differences in Saline Water Drip Irrigation Cotton Fields

In arid areas, fresh water resources are insufficient, and agricultural water mainly depends on shallow saline groundwater. However, long-term saline irrigation will cause soil salt accumulation and soil environment deterioration, which is not conducive to crop growth. In this study, based on the long-term irrigation of fresh water (0.35 dS·m-1, FW) and saline water (8.04 dS·m-1, SW), biochar (3.7 t·hm-2, BC) and straw (6 t·hm-2, ST) were added to the soil by an equal-carbon design. The aim was to clarify the effects of biochar and straw returning on the physical and chemical properties and microbial community structure of salinized soil. The results showed that saline irrigation significantly increased soil water content, electrical conductivity, available phosphorus, and total carbon content but significantly decreased pH value and available potassium content. The contents of available phosphorus, available potassium, and total carbon in soil were significantly increased by biochar and straw returning, but the conductivity value of soil irrigated with saline water was significantly decreased. The dominant bacteria in each treatment were Proteobacteria, Actinomycetes, Acidobacteria, Chloromycetes, and Blastomonas. Saline water irrigation significantly increased the relative abundance of Blastomonas and Proteobacteria but significantly decreased the relative abundance of Acidobacteria and Actinobacteria. Under the condition of fresh water irrigation, the relative abundance of Chlorocurvula was significantly reduced by the return of biochar. Straw returning significantly increased the relative abundance of Proteobacteria but significantly decreased the relative abundance of Acidobacteria, Actinomyces, Chloromyces, and Blastomonas. Under saline irrigation, the relative abundance of Chlorocurvula and Blastomonas were significantly reduced by biochar return to field. Straw returning significantly increased the relative abundance of Proteobacteria but significantly decreased the relative abundance of Acidobacteria, Actinomyces, Chloromyces, and Blastomonas. LEfSe analysis showed that saline irrigation decreased the potential markers and functional numbers of soil microorganisms.Under saline irrigation, biochar returning increased the number of potential markers and functions of soil microorganisms. Straw returning to field increases the number of potential markers of soil microorganisms. RDA results showed that soil microbial community and functional structure were significantly correlated with EC1:5, SWC, and pH. Saline water irrigation will deteriorate the soil environment, which is not conducive to agricultural production, among which EC1:5, SWC, and pH are important factors driving changes in soil microbial community and functional structure. Using biochar and straw to return to the field can reduce the harm of salt to soil and crops, laying a foundation for improving agricultural productivity.

Surface conductivity of clays

SUMMARY Clay minerals are extensively used in a wide range of applications. In particular, clay-bearing formations are considered as suitable radioactive waste repository. Electrical resistivity tomography is an appropriate tool to monitor the properties of clay-bearing locations. However, an inherent drawback of a conventional resistivity survey is its ambiguity in distinguishing between the effects of groundwater salinity, clay content and porosity. A discrimination can be achieved on the basis of the induced polarization method that provides a complex conductivity. The main purpose of this study is the investigation of the complex conductivity of clay samples with a special focus on the contribution of surface conductivity produced by an excess of ions in the electrical double-layer coating the solid particles. Six clay mixtures were selected that include an almost pure kaolinite sample, a sample consisting of a mixture of kaolinite, illite and smectite, a crushed saponite breccia, a Ca-bentonite sample and two illite clay samples. Besides the enriched kaolinite, the other samples are natural geomaterials that contain more than 40 weight per cent clay minerals. The mineralogical compositions of the samples were determined by quantitative X-ray diffraction analysis. The clay powder was mixed with a varying volume of sodium chloride solution to get plastic state clay samples with varying water content. The samples were investigated by the spectral induced polarization method in a frequency range between 1 mHz and 1 kHz. The resulting complex conductivity spectra indicate a decrease of the real part of the electrical conductivity with rising water content for the illite, bentonite and saponite breccia samples. The overall conductivity of these clay samples is dominated by their surface conductivity. In contrast, the electrical conductivity of kaolinite and kaolinite–illite mixture does not show any significant changes with the water content. For all samples, the imaginary part of electrical conductivity increases at low water content. The real part of the surface conductivity indicates a linear dependence on the volumetric clay content. The slope of this linear relationship can be used to distinguish the types of clay. The ratio between imaginary conductivity and surface conductivity, which decreases with increasing clay content, proves to be a suitable parameter that characterizes the connectivity of clay aggregates in the sample. The surface conductivity of the pure kaolinite sample has been determined in an additional multisalinity experiment. The resulting surface conductivity is in good agreement with the experiment of varying water content. The multisalinity experiment has shown that the resulting petrophysical parameters depend on the procedure of sample packing, which may lead to anisotropy. The effect of anisotropy is attributed to the alignment of the plate-like kaolinite particles in the course of the packing and consolidation procedure.

Joint inversion of induced polarization and hydraulic tomography data for hydraulic conductivity imaging

SUMMARY For accurate modelling of groundwater flow and transport processes within an aquifer, precise knowledge about hydraulic conductivity K and its small-scale heterogeneities is fundamental. Methods based on pumping tests, such as hydraulic tomography (HT), allow for retrieving reliable K-estimates, but are limited in their ability to image structural features with high resolution, since the data from time-consuming hydraulic tests are commonly sparse. In contrast, geophysical methods like induced polarization (IP) can potentially yield structural images of much higher resolution, but depend on empirical petrophysical laws that may introduce significant uncertainties to the K-estimation. Therefore, this paper presents a joint inversion procedure for both HT and IP data, which allows for combining the complementary abilities of both methods. Within this approach, a traveltime inversion is applied to the HT data, while the IP inversion is based on a full-decay time-domain forward response, as well as a reparametrization of the Cole–Cole model to invert for K directly. The joint inversion is tested on a synthetic model mimicking horizontally layered sediments, and the results are compared with the individual HT and IP inversions. It is shown that jointly inverting both data sets consistently improves the results by combining the complementary sensitivities of the two methods, and that the inversion is more robust against changes in the experimental setups. Furthermore, we illustrate how a joint inversion approach can correct biases within the petrophysical laws by including reliable K-information from hydraulic tests and still preserving the high-resolution structural information from IP. The different inversion results are compared based on the structural similarity index (SSIM), which underlines the robustness of the joint inversion compared to using the data individually. Hence, the combined application of HT and IP within field surveys and a subsequent joint inversion of both data sets may improve our understanding of hydraulically relevant subsurface structures, and thus the reliability of groundwater modelling results.

Hydrochemical evolution of groundwater in a river corridor: the compounded impacts of various environmental factors

Various environmental factors could induce groundwater hydrochemical evolution, which should be considered when addressing groundwater environmental issues. However, the complex interactions among the environmental controls remain unknown for a groundwater flow system spanning multiple geological units. To fill this gap, we conducted a study on the groundwater hydrochemical evolution in the fluvial corridor of the Fen River, Northern China, utilizing a combination of hydrogeochemical and multiple isotope methods. Results reveal that the groundwater in the corridor has significantly degraded due to high concentrations of SO42-\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ ext{SO}}_{4}^{2 - }$$\\end{document}, NO3-\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ ext{NO}}_{3}^{ - }$$\\end{document}, Cl−, or F−. We find an unordered evolution of the hydrochemical composition of groundwater along this corridor. These evolutions are driven by mineral dissolution/precipitation, dedolomitization, and cation exchange processes. Human activities play a significant role, with notable contributions including NO3-\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ ext{NO}}_{3}^{ - }$$\\end{document} fluxes from agricultural fertilizers, manure, and sewage, as well as SO42-\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ ext{SO}}_{4}^{2 - }$$\\end{document} fluxes arising from coal mining activities. The combination between karst spring areas and faults/uplifts, between coal-bearing strata and mining activities fosters the mixing of karst water/mine water with shallow groundwater, promoting dramatic hydrochemical evolution of groundwater. The flat terrain and the natural blockage formed by mountains significantly enhance water–rock interactions and groundwater evaporation by slowing groundwater flow. The contribution of evaporation on groundwater salinity ranges from ~ 0.2% to 4.8%, highlighting its importance in the groundwater hydrochemical evolution. This study unravels the multifaceted nature of groundwater hydrochemical evolution. It emphasizes the four types of environmental controls, including hydrogeochemical processes, climate factors, human activities, and variations in geological settings, can be equally important, which is usually ignored. The findings enrich our understanding of groundwater evolution and highlight the challenges encountered in regions composed of diverse geological units.

The Hydraulic Evolution of Groundwater-Fed Pit Lakes After Mine Closure.

Open pit mining frequently requires regional water tables to be lowered to access ore deposits. When mines close, dewatering ceases allowing the water table to recover. In arid and semi-arid mining regions, the developing pit lakes are predominantly fed by groundwater during this recovery phase and pit lakes develop first into "terminal sinks" for the surrounding groundwater system. With time, the re-establishment of regional hydraulic gradients can cause pit lakes to develop into throughflow systems, in which pit lake water outflows into adjacent aquifers. In this study, we use numerical groundwater modeling to aid process understanding of how regional hydraulic gradients, aquifer properties, net evaporation rates, and pit geometry determine the hydraulic evolution of groundwater-fed pit lakes. We find that before the recovery of the regional water table to its new equilibrium, pit lakes frequently transition to throughflow systems. Throughflow from pit lakes to downstream aquifers can develop within two decades following cessation of dewatering even under low hydraulic gradients (e.g., 5 × 10-4) or high net evaporation rates (e.g., 2.5 m/year). Pit lakes remain terminal sinks only under suitable combinations of high evaporation rates, low hydraulic gradients, and low hydraulic conductivities. In addition, we develop an approximate analytical solution for a rapid assessment of the hydraulic status of pit lakes under steady-state conditions. Understanding whether pit lakes remain terminal sinks or transition into throughflow systems largely determines the long-term water quality of pit lakes and downstream aquifers. This knowledge is fundamental for mine closure and planning post-mining land use.