In simulations of groundwater flow through dipping aquifers, layers of model cells are often "deformed" to follow the top and bottom elevations of the aquifers. When this approach is used in MODFLOW, adjacent cells withi...In simulations of groundwater flow through dipping aquifers, layers of model cells are often "deformed" to follow the top and bottom elevations of the aquifers. When this approach is used in MODFLOW, adjacent cells within the same model layer are vertically offset from one another, and the standard conductance-based (two-point) formulation for flow between cells does not rigorously account for these offsets. The XT3D multi-point flow formulation in MODFLOW 6 is designed to account for geometric irregularities in the grid, including vertical offsets, and to provide accurate results for both isotropic and anisotropic groundwater flow. A recent study evaluated the performance of the standard formulation and XT3D using a simple, synthetic benchmark model of a steeply dipping aquifer. Although XT3D generally improved the accuracy of flow simulations relative to the standard formulation as expected, neither formulation produced accurate flows in cases that involved large vertical offsets. In this paper, we explain that the inability of XT3D to produce accurate flows in the steeply dipping aquifer benchmark was not due to an inherent limitation of the flow formulation, but rather to the limited cell connectivity inherent in the most commonly used discretization packages in MODFLOW 6. Furthermore, we demonstrate that XT3D is able to produce the expected accuracy when adequate cell connectivity is introduced using MODFLOW's unstructured grid type and the aquifer is discretized vertically using at least two model layers.
Wetlands, as crucial terrestrial carbon reservoirs, have recently suffered severe degradation due to intense human activities. Lacustrine sediments serve as vital indicators for understanding wetland environmental change...Wetlands, as crucial terrestrial carbon reservoirs, have recently suffered severe degradation due to intense human activities. Lacustrine sediments serve as vital indicators for understanding wetland environmental changes. In the current paper, porewater samples were extracted from lacustrine sediment in three boreholes with a depth of ~75 cm in the Huixian karst wetland, southwest China, to study the chemical and dissolved inorganic carbon (DIC) evolution under anthropogenic influence. Two boreholes are situated beneath the Mudong Lake, while the other one is in the degraded wetland area. The results show that porewater in the central region of Mudong Lake is natural HCO-Ca type water and recharged by karst groundwater as evidenced by depleted H -O isotopes. Methanogenesis prevails in this area, suggested by positive δC values ranging from 4.29‰ to 7.05‰. However, shallow porewater at the western edge of Mudong Lake and porewater in the degraded wetland exhibit significantly higher concentrations of NO and SO , resulting from the agricultural input and recharged groundwater influenced by oxidation of pyrite. These processes lead to a decrease in methane production and generate DIC through degradation of organic fertilizer and carbonate weathering by sulfuric acid, thereby significantly altering porewater δC values. Two types of DIC mixing processes were observed based on the increasing δC values with depth, which can be attributed to the unique karst groundwater subsystems. This work highlights the potential impact of human-induced porewater chemical variations on the fate of DIC, particularly in karst wetland environments.
Continental glaciations during the Pleistocene Epoch created complex systems of aquifers and aquitards across many northern regions of the Earth. The low hydraulic conductivities of glacial till aquitards suggest that li...Continental glaciations during the Pleistocene Epoch created complex systems of aquifers and aquitards across many northern regions of the Earth. The low hydraulic conductivities of glacial till aquitards suggest that limited recharge will reach the underlying aquifers, potentially preserving old groundwaters. Here, we characterize the recharge history in intertill and buried valley aquifers in Saskatchewan, Canada using C, H, He δH, δO, and major ions. Intertill aquifers with depths of <30 m had corrected C ages ranging from 0 to 15.5 ka. These aquifers also contained H and/or elevated NO in some locations, indicating that a component of modern recharge had mixed with older water. A single sample from the Judith River bedrock aquifer in the region had a corrected C age of 10.2 ka and elevated NO. Samples from buried valley aquifers with depths of 89 to 123 m contained older waters with ages >38 ka in some locations, indicating that recharge occurred before the last glacial advance over the region. While measuring tracers that cover a wide range of ages is necessary to understand these flow systems, δH and δO were less diagnostic because values of modern winter precipitation overlapped with groundwaters with a wide range of ages. The range of ages present in the intertill aquifers of the region indicates that these systems are currently being recharged, which indicates some development of groundwater resources is possible but also points to a need for groundwater protection measures.
Quantifying lacustrine groundwater discharge (LGD) is important for understanding the dynamics of lake ecosystems and their expansion. This study focuses on Lake Qinghai, employing radium isotope models to evaluate the c...Quantifying lacustrine groundwater discharge (LGD) is important for understanding the dynamics of lake ecosystems and their expansion. This study focuses on Lake Qinghai, employing radium isotope models to evaluate the contributions of both shallow and deep groundwater. The data indicate that the activity of Ra and Ra demonstrates a pronounced gradient, decreasing from the shoreline to the center of Lake Qinghai. Additionally, vertical stratification characteristics were observed. The spatial distribution of radium isotope activity suggests that there is discharge of both shallow and deep groundwater into the lake. Deep groundwater migrates slowly and its apparent age reflects the time elapsed since the water became enriched in Ra and was isolated from the source, in the study system this age is estimated to be 10.1 d. In contrast, shallow groundwater displayed varied apparent ages in different regions: 7.9 d in the north, 13.1 d in the south, and 7.4 d in the southeastern area of the lake. The LGDs of shallow groundwater discharge in the north, south, and southeast areas of Lake Qinghai were estimated by Ra as 1.89 × 10 to 2.69 × 10 m/d, 3.25 × 10 to 3.99 × 10 m/d, and 4.51 × 10 to 6.33 × 10 m/d, respectively. For deep groundwater, the LGD was 0.16 × 10 to 0.29 × 10 m/d. Annually, the total LGD fluxes of shallow and deep groundwater are 27.86 × 10 to 37.59 × 10 m/year and 0.58 × 10 to 1.06 × 10 m/year, respectively. This study is the first to evaluate shallow and deep groundwater discharge around the lake. Understanding these discharge dynamics is essential for developing effective management strategies to preserve lake environments.
Sea water intrusion (SWI) simulators are essential tools to assist the sustainable management of coastal aquifers. These simulators require the solution of coupled variable-density partial differential equations (PDEs),...Sea water intrusion (SWI) simulators are essential tools to assist the sustainable management of coastal aquifers. These simulators require the solution of coupled variable-density partial differential equations (PDEs), which reproduce the processes of groundwater flow and dissolved salt transport. The solution of these PDEs is typically addressed numerically with the use of density-dependent flow simulators, which are computationally intensive in most practical applications. To this end, model surrogates are generally developed as substitutes for full-scale aquifer models to trade off accuracy in exchange for computational efficiency. Surrogates represent an attractive option to support groundwater management situations in which fast simulators are required to evaluate large sets of alternative pumping strategies. Reduced-order models, a sub-category of surrogate models, are based on the original model equations and may provide quite accurate results at a small fraction of computational cost. In this study, a variable-density flow reduced-order model based on proper orthogonal decomposition (POD) and utilizing a fully coupled flow and solute-transport model is implemented with a finite-difference (FD) approach for simulating SWI in coastal aquifers. The accuracy and computational efficiency of the FD-POD approach for both homogeneous and-more realistic-heterogeneous systems are investigated using test cases based on the classic Henry's problem (Henry 1964). The findings demonstrate that the combined FD-POD approach is effective in terms of both accuracy and computational gain and can accommodate the output of the most popular variable-density flow models, such as those from USGS's MODFLOW family.
Due to increasing global demand for fresh water, it is increasingly necessary to understand how aquifer pumping affects groundwater chemistry. However, comprehensive predictive relationships between pumping and groundwat...Due to increasing global demand for fresh water, it is increasingly necessary to understand how aquifer pumping affects groundwater chemistry. However, comprehensive predictive relationships between pumping and groundwater quality have yet to be developed, as the available data, which are often collected over inconsistent time intervals, are poorly suited for long-term historical correlation studies. For example, we needed an adequate statistical method to better understand relationships between pumping rate and water quality in the City of Norman (OK, USA). Here we used the interval-scaled change in mean pumping rate combined with the Quadrant method to examine correlations between pumping rates and changes in trace metal concentrations. We found that correlations vary across the study area and are likely dependent on a variety of factors specific to each well. Comparing the Quadrant method to the commonly used Kendall's tau correlation, which requires different assumptions about aquifer behavior, the methods produced similar correlations when sample sizes were large and the time interval between samples was relatively short. Sample sizes were then artificially restricted to determine correlation reproducibility. Despite being less reproducible overall, the Quadrant method was more reproducible when there were large time intervals between samples and very small sample sizes (n ~ 4), but not as reproducible as significant (p ≤ 0.1) Kendall's tau correlations. Therefore, the Quadrant method may be useful for further investigating the effects of pumping in cases where Kendall's tau does not produce significant correlations.
Drilling wells in unconsolidated formations is commonly undertaken to extract drinking water and other applications, such as aquifer thermal energy storage (ATES). To increase the efficiency of an ATES system, the drilli...Drilling wells in unconsolidated formations is commonly undertaken to extract drinking water and other applications, such as aquifer thermal energy storage (ATES). To increase the efficiency of an ATES system, the drilling campaigns are targeting greater depths and enlarging the wellbore diameter in the production section to enhance the flow rates. In these cases, wells are more susceptible to collapse. Drilling fluids for shallow formations often have little strengthening properties and, due to single-string well design, come into contact with both the aquifer and the overburden. Drilling fluids and additives are experimentally investigated to be used to improve wellbore stability in conditions simulating field conditions in unconsolidated aquifers with a hydraulic conductivity of around 10 m/d. The impact on wellbore stability is evaluated using a new experimental setup in which the filtration rate is measured, followed by the use of a fall cone penetrometer augmented with an accelerometer to directly test the wellbore strengthening, and imaging with a scanning electron microscope (SEM) to investigate the (micro)structure of the filter cakes produced. Twelve drilling fluids are investigated with different concentrations of bentonite, polyanionic cellulose (PAC), Xanthan Gum, calcium carbonate (CaCO), and aluminum chloride hexahydrate ([Al(HO)]Cl). The filtration results indicate that calcium carbonate, average d <20 μm, provides pore throat bridging and filter cake formation after approximately 2 min, compared to almost instantaneous discharge when using conventional drilling fluids. The drilling fluid containing 2% [Al(HO)]Cl forms a thick (4 mm) yet permeable filter cake, resulting in high filtration losses. The fall cone results show a decrease of cone penetration depth up to 20.78%, and a 40.27% increase in deceleration time while penetrating the sample with CaCO compared with conventional drilling fluid containing bentonite and PAC, indicating a significant strengthening effect. The drilling fluids that contain CaCO, therefore, show high promise for field implementation.
Our recent steady-state mass-balance modeling suggests that most global carbonic-acid weathering of silicate rocks occurs in the vadose zone of aquifer systems not on the surface by atmospheric CO. That is, the weatherin...Our recent steady-state mass-balance modeling suggests that most global carbonic-acid weathering of silicate rocks occurs in the vadose zone of aquifer systems not on the surface by atmospheric CO. That is, the weathering solute flux is nearly equal to the total global continental riverine carbon flux, signifying little atmospheric weathering by carbonic acid. This finding challenges previous carbon models that utilize silicate weathering as a control of atmospheric CO levels. A robust analysis utilizing global estimates of groundwater carbon concentration generated by a geospatial machine learning algorithm was coupled with recharge flux in a geographic information system environment to yield a total global groundwater carbon flux of between 0.87 and 0.96 Pg C/year to the surface environment. On discharging to the surface, the carbon is speciated between 0.01 and 0.11 Pg C/year as CaCO; 0.35 and 0.38 Pg C/year as CO gas; and 0.49 and 0.51 Pg C/year as dissolved HCO . This total weathering carbon flux was calculated for direct ocean discharge (0.030 Pg C/year); endorheic basins (0.046 Pg C/year); cold-wet exorheic basins (0.058 Pg C/year); warm-dry exorheic basins (0.072 Pg C/year); cold-dry exorheic basins (0.115 Pg C/year); and warm-wet exorheic basins (0.448 Pg C/year).
The potential performance of a hypothetical colloidal-activated carbon (CAC) in situ remedy for perfluorooctanoic acid (PFOA) and perfluorooctane sulfonic acid (PFOS) in groundwater in coastal zones was evaluated using e...The potential performance of a hypothetical colloidal-activated carbon (CAC) in situ remedy for perfluorooctanoic acid (PFOA) and perfluorooctane sulfonic acid (PFOS) in groundwater in coastal zones was evaluated using estimated hydrogeologic and geochemical parameters for a coastal site in the United States. With these parameters, a reactive transport model (ISR-MT3DMS) was used to assess the effects of tidal fluctuations and near-shore geochemistry on CAC performance. The average near-shore ionic strength of 84 mM at the site was conservatively estimated to result in an increase in the adsorption of PFOA to CAC by about 50% relative to non-coastal sites with ionic strength <10 mM. The modeling also confirmed the hypothesis that tidally induced groundwater flow reversals near the shore would result in the accumulation of PFOA at the downgradient edge of the CAC zone. Slow desorption of PFOA from this downgradient CAC boundary may sustain downgradient plume concentrations above a strict cleanup criterion (e.g., USEPA MCL of 0.004 μg/L), for decades; however, there was still a large PFOA mass flux reduction (>99.9%) achieved after several decades at the shore. CAC longevity was substantially greater for PFOS with a similar source concentration; however, the higher PFOS distribution coefficient (K) in soil downgradient from the CAC zone resulted in substantially longer flushing times. It is recommended that short-term remedial action objectives for CAC remedies at coastal sites be based on mass flux reduction targets over a period of several decades, given the demonstrated challenges in trying to achieve very low cleanup criteria downgradient of a CAC zone in the short term.
Small island communities often rely on groundwater as their primary source of fresh water. However, the limited land area and high proportion of coastal zones pose unique challenges to groundwater management. A detailed...Small island communities often rely on groundwater as their primary source of fresh water. However, the limited land area and high proportion of coastal zones pose unique challenges to groundwater management. A detailed understanding of the subsurface structure can provide valuable insights into aquifer structure, groundwater vulnerability, saltwater intrusion, and the location of water resources. These insights can guide groundwater management strategies, for example, pollution regulation, promotion of sustainable agriculture, establishment of coastal buffer zones, and re-naturalization of land cover. Ordinarily, structural characterization relies on geological mapping and boreholes, however, such approaches can have insufficient spatial resolution to aid groundwater management. In this study, transient electromagnetic (TEM) methods are used to map the subsurface of a small, 13.2 km, Danish Island. The approach successfully identified two previously unknown paleochannels, where the interface between Quaternary aquifer units and an underlying Paleogene Clay aquiclude had maximum depths of 100 and 160 m below sea level. Before this, the interface was assumed to be 15 to 25 m below sea level: therefore, these paleochannels present substantial potential groundwater resources. Resolving geological heterogeneity within the Quaternary deposits was less successful and future work will focus on addressing these limitations. Nonetheless, in several locations, evidence of saltwater intrusion was observed within the Quaternary units. This work demonstrates how TEM mapping can identify water resources, define aquifer boundaries, and aid water management decisions. Such approaches could be applied in other areas, particularly small islands, where similar groundwater challenges exist.
The step-drawdown pumping test often experiences a transition from confined to unconfined conditions due to the continuously increasing pumping rate. However, the current well hydraulics model has not accurately interpre...The step-drawdown pumping test often experiences a transition from confined to unconfined conditions due to the continuously increasing pumping rate. However, the current well hydraulics model has not accurately interpreted this phenomenon. In this study, we developed an analytical solution to address the confined-unconfined conversion in step-drawdown pumping tests based on Girinskii's potential and superposition theory. Additionally, a field step-drawdown pumping test featuring confined-unconfined conversion was conducted to apply the proposed analytical solution. The particle swarm optimization algorithm was employed to simultaneously estimate multiple parameters. The results demonstrate that the newly proposed solution provides a better fit to the observed drawdown in the pumping well compared to previous models. The hydrogeological parameters (K, S), well loss coefficient (B), and critical time for confined-unconfined conversion (t) were estimated to be K = 7.15 m/d, S = 6.65 × 10, B = 7.48 × 10, and t = 1152 min, respectively. Neglecting the confined-unconfined conversion in step-drawdown pumping tests leads to underestimation of drawdown inside the pumping well due to an overestimation of the aquifer thickness. After the conversion from confined to unconfined conditions, the estimated well loss coefficient decreased by 88% compared to its pre-conversion value. This highlights the necessity of adjusting the well loss coefficient in the step-drawdown pumping test model to account for confined-unconfined conversion. In summary, this study introduces a new method for interpreting parameters in step-drawdown pumping tests and provides field validation for its effectiveness.
This study enhances the understanding of riverbank filtration and improves management of the Mississippi River valley alluvial (MRVA) aquifer during a managed aquifer recharge (MAR) pilot project at Shellmound, MS. Using...This study enhances the understanding of riverbank filtration and improves management of the Mississippi River valley alluvial (MRVA) aquifer during a managed aquifer recharge (MAR) pilot project at Shellmound, MS. Using high-resolution electrical resistivity tomography (ERT) and self-potential (SP) geophysical methods, we characterized the heterogeneous MRVA aquifer and monitored groundwater flow near a pumping well. ERT was used to provide detailed spatial characterization, filling gaps left by airborne electromagnetic (AEM) data and soil boring logs, while SP techniques were used to monitor groundwater flow, predict drawdown trends, and investigate surface-groundwater interactions. Results showed that SP signals were influenced by groundwater flow, river infiltration, and water mixing due to pumping disturbance of natural geochemical stratification, with significant river interaction observed after 1 h of pumping. The integration of ERT and SP methods revealed lithologic heterogeneity, explaining greater drawdowns on the northern side of the well and increased flow from the riverside. This comprehensive approach offers valuable insights into aquifer management and sustainability.
Research into land subsidence caused by groundwater withdrawal is hindered by the availability of measured heads, subsidence, and forcings. In this paper, a parsimonious, linked data-driven and physics-based approach is...Research into land subsidence caused by groundwater withdrawal is hindered by the availability of measured heads, subsidence, and forcings. In this paper, a parsimonious, linked data-driven and physics-based approach is introduced to simulate pumping-induced subsidence; the approach is intended to be applied at observation well nests. Time series analysis using response functions is applied to simulate heads in aquifers. The heads in the clay layers are simulated with a one-dimensional diffusion model, using the heads in the aquifers as boundary conditions. Finally, simulated heads in the layers are used to model land subsidence. The developed approach is applied to the city of Bangkok, Thailand, where relatively short time series of head and subsidence measurements are available at or near 23 well nests; an estimate of basin-wide pumping is available for a longer period. Despite the data scarcity, data-driven time series models at observation wells successfully simulate groundwater dynamics in aquifers with an average root mean square error (RMSE) of 2.8 m, relative to an average total range of 21 m. Simulated subsidence matches sparse (and sometimes very noisy) land subsidence measurements reasonably well with an average RMSE of 1.6 cm/year, relative to an average total range of 5.4 cm/year. Performance is not good at eight out of 23 locations, most likely because basin-wide pumping is not representative of localized pumping. Overall, this study demonstrates the potential of a parsimonious, linked data-driven, and physics-based approach to model pumping-induced subsidence in areas with limited data.
Numerical modeling offers a valuable alternative to analytical solutions for pumping test analysis. However, little is known about how discretization impacts results accuracy and runtime. This study presents a systematic...Numerical modeling offers a valuable alternative to analytical solutions for pumping test analysis. However, little is known about how discretization impacts results accuracy and runtime. This study presents a systematic method for defining the spatiotemporal discretization of pumping test numerical models based on dimensionless parameters. Two types of analysis are considered: one where observations are made in the pumping well, and another one where observations are made in different wells. The influence of the discretization parameters on results accuracy and runtime is investigated and an optimal set of parameters is determined that minimizes runtime while maintaining the maximum error under 1% for an "average" aquifer. Lower runtimes are achieved when the analysis focuses on the pumping well, which is attributed to the steady-state analytical solution approximating drawdown in the well in the numerical scheme employed. Additional tests demonstrate the robustness of the derived set of parameters in different configurations.
Automated water level measurements collected using vented pressure transducers in deep wells screened across the water table may exhibit a greater response to barometric pressure changes than the true water level. The ca...Automated water level measurements collected using vented pressure transducers in deep wells screened across the water table may exhibit a greater response to barometric pressure changes than the true water level. The cause was hypothesized to be disequilibrium in barometric pressure between the wellbores and land surface due to air exchange with the deep vadose zone. In this study, vented and nonvented pressure transducers were installed and operated simultaneously in two deep wells screened across the water table. A vent tube open to the atmosphere at land surface allowed for barometric compensation of the vented transducers. Two nonvented transducers were installed in each well, one submerged in the water and one above the water surface. The difference in readings allowed for barometric compensation. Manual measurements were also collected. It was confirmed that measurements from the vented transducers exhibited greater variability in response to barometric pressure changes than the nonvented transducers and manual measurements. Comparison of the downhole barometric pressure measurements to values from a nearby meteorology station showed the response in the wells to changes in barometric pressure was time-lagged and attenuated. Thus, the reference pressure from land surface supplied to the vented transducers was not representative of the air pressure within the wells. This caused fluctuations of the transducer readings in response to barometric pressure changes to be greater than the true water level change. This issue can be resolved by the use of nonvented pressure transducers.