Groundwater最新文献

A sophisticated understanding of the three-dimensional distribution of silt- and clay-rich bodies of strata (elements) in aquifers is critical given that they not only have the potential to act as aquitards or semi-confining units and vertically partition groundwater flow into separate aquifer zones, but also provide lateral barriers to groundwater flow, impacting contaminant distribution and groundwater flow dynamics. Additionally, when in prolonged contact with dense nonaqueous phase liquid (DNAPL) or contaminated groundwater, fine-grained elements may become storage zones for contaminant mass via matrix diffusion and thus serve as long-term secondary sources of contamination to groundwater that can confound remediation strategies and render remedy performance projections unreliable. The stratigraphic architecture of aquifer systems, including fine-grained facies architecture, is complex but is not random and can be effectively predicted through application of facies models. This paper reviews depositional models (“facies models”) for common depositional environments with a focus on alluvial end-members of braided fluvial, meandering fluvial, and alluvial fan facies models. We examine the facies models from the perspective of aquitards and present case studies to provide an overview of the expected aquitard dimensions and characteristics. The critical yet underappreciated role of the paleosol as a potential aquitard is also examined, and basic criteria for differentiating ancient floodplain clay units with high lateral continuity from other laterally discontinuous clay units are provided.

Leaching of marine salinity in the porewater of glaciomarine muds is one precursor to landslide hazard. In this study, groundwater modeling is used to quantify vertical groundwater flow, constrain paleosalinity, and characterize past and future progression of leaching with depth in Champlain Sea sediments. The Breckenridge Creek site, ~15 km northwest of Ottawa, Canada, was cored within a thick sequence (up to 98 m) of Champlain Sea muds that form a regional aquitard in the St. Lawrence Lowlands and Ottawa Valley. Porewater chloride concentrations ([Cl]), up to 12,250 mg/kg, and δ¹⁸O as high as −7.18‰, indicate remnant seawater. One-dimensional groundwater transport modeling simulates porewater [Cl] and δ¹⁸O with depth simultaneously and constrains specific discharge, q, from 2.40 to 2.51 mm/a. Groundwater transport modeling and three-component mixing of seawater, glacial meltwater and meteoric water constrain the range of initial [Cl] between 14,000 and 15,700 mg/kg (72–80% seawater) and initial δ¹⁸O between −5.99 and −5.61‰. The glacial meltwater component of Champlain Sea bottom waters at the Breckenridge site has a maximum δ¹⁸O value of −22.4‰. Downward leaching to the salinity threshold of <2 g/L for geotechnical sensitivity development reached a depth of 20.6 m. Modeling indicates the leaching front currently progresses at a rate of 2.5 m/1000 years, slower than advection of freshwater infiltration because of upward diffusion and dispersion of marine solutes. Notably for landslide hazard, the highest measurements of geotechnical sensitivity coincide with the leached zone.

The compaction simulation of compressible delay interbed is an important part of land subsidence simulation. Currently, the most widely used MODFLOW software has two modules, SUB and CSUB, both of which can simulate compressible delay interbed. The difference lies in that the head diffusion equation of the SUB module is based on the principle of head change, while CSUB can use either head change or geological stress variation principles. When based on the principle of geostress variation, the CSUB method is more physically reasonable. However, its limitation lies in the fact that, when solving the diffusion equation for compressible delay interbeds, it does not account for the effects of variations in the discrete nodal cell thickness and hydraulic conductivity of the interbed. This study improves the solution method for the head diffusion equation of compressible delay interbeds based on the principle of geostress variation. The Kozeny–Carman equation was introduced to establish a relationship between the hydraulic conductivity and porosity of the interbeds, while variations in the thickness of discrete nodal cells were also incorporated into the solution process. Collectively, these improvements lead to a more rigorous approach. To verify the effectiveness of the proposed simulation method, three representative test cases were developed and comprehensively compared with the CSUB results. The results indicate that notable discrepancies emerge between the two approaches when the interbed undergoes substantial compression, whereas the method proposed in this study effectively prevents the occurrence of “overcompaction” within the interbed.

Aquitards play critical roles in a variety of hydrogeologic processes. Despite their importance, coverage of aquitards in introductory hydrogeology textbooks is generally limited. This paper provides examples of classroom and field activities designed with an aquitard focus for instructors wishing to supplement textbook content. These activities emphasize high-resolution head profiles. Examining head profiles prompts students to think about how aquitards influence head with depth and conversely how these plots can be used to delineate and characterize aquitards. During a classroom activity, students explore the connection between changes in vertical gradient and changes in hydraulic conductivity by sketching conceptual head profiles based on given boundary conditions and several aquifer/aquitard scenarios. In a companion field exercise, students measure high-resolution head profiles using CMT multilevel systems at an outdoor learning laboratory. Students compare the high-resolution head profiles to lower resolution profiles they obtain from clusters of conventional wells. The field exercise provides students with a tactile experience that can help build intuition for vertical head changes, practice interpreting aquitards from head profiles, and an example of how lower resolution head profiles may create uncertainty in aquitard delineation and vertical gradient estimates. A paleosol at the site forms a prominent aquitard providing a unique basis for discussions about the geology of aquitards and characteristics influencing aquitard integrity. Regardless of the approach used, incorporating more aquitard content into hydrogeology courses at all levels will be beneficial for future hydrogeologists tackling a range of issues from sustainable water supplies to waste disposal.

The transboundary Hueco Bolson aquifer is the primary water supply for El Paso, Texas, and Ciudad Juárez, Mexico as well as U.S. Army Fort Bliss and smaller cities in Texas and Mexico. Binational groundwater data exchanges between the United States and Mexico that have focused on the aquifer have evolved over more than a century, shaped by scientific, political, and social dynamics. The history can be viewed through distinct periods. Efforts to understand the groundwater resource began with early reconnaissance surveys and have evolved through successive efforts to refine the characterization of groundwater flowpaths, residence times, surface water–groundwater interactions, and aquifer salinization. Lessons from these efforts highlight the importance of perseverance, mutual respect, and formal agreements, such as the 1999 memorandum between El Paso Water Utilities and the Junta Municipal de Agua y Saneamiento de Juárez, in sustaining long-term cooperation. The accumulated datasets chronicle the evolution of hydrogeologic conceptual and mathematical models while providing a foundation for ongoing research, sustainable water-use strategies, and the long-term stewardship of the Hueco Bolson aquifer shared by El Paso and Ciudad Juárez. The Hueco Bolson case demonstrates how long-term, cooperative data collection can improve scientific understanding and management of complex transboundary aquifer systems.

The concentration of fluoride (F⁻) in many geothermal waters worldwide exceeds the World Health Organization's drinking water guideline of 1.5 mg/L, highlighting a widespread water quality issue in geothermal systems. It is essential to predict the mobility of F⁻ in geothermal reservoirs compared to that of chloride (Cl⁻). Column experiments were conducted at 45°C under hydrostatic pressures from atmospheric pressure to 12 MPa to simulate geothermal reservoir conditions, with the computer code CXTFIT 2.1 used to fit the data and determine transport and adsorption parameters of Cl⁻ and F⁻. The effect of hydrostatic pressure on F⁻ transport was evaluated. Results showed that the Convection–Dispersion Equation (CDE) accurately fits the breakthrough curves of Cl⁻ and higher pressures enhance its dispersion transport. In contrast, F⁻ transport under atmospheric pressure shows rate-limited non-equilibrium, with the Two-Site Model (TSM) fitting F⁻ curves better. However, when pressure exceeds 6 MPa, equilibrium adsorption becomes more pronounced, and the equilibrium CDE better characterizes F⁻ transport. Higher pressure also increases the retardation factor (R) of F⁻, suggesting more instantaneous adsorption sites under these conditions. The release of hydroxide (OH⁻) during F⁻ adsorption increases pH values in the effluent notably. Moreover, greater pH fluctuations occur with increasing pressure due to the enhanced adsorption of F⁻ by the packed matrix. The findings would be helpful for understanding the solute transport and adsorption mechanism of fluoride under varying hydrostatic pressures in groundwater and geothermal systems.

The tracer flowmeter and depth-dependent sampler (TFDDS) has been used to characterize flow and chemistry in hundreds of public supply wells. The TFDDS surveys address production and water quality compliance issues for utilities facing mounting treatment costs. Surveyed wells often show a stratified distribution for constituents of concern that were underestimated and poorly discretized using conventional zone testing in the pilot hole. Elevated concentrations of metals, semimetals, and radionuclides typically occurring at or near lithologic boundaries between fine- and coarse-grained sediments are missed since the primary focus of the zone test is to estimate the yield of the permeable sediments that are centered within a coarse unit. To minimize the risk of constructing new wells that fail to produce compliant drinking water relative to regulatory standards, we propose a new approach that vertically stacks flow and mass balance chemistry profiling of long screened test wells (LSTWs) under pumping conditions prior to installing the more costly public supply well. This approach integrates the TFDDS with a spectrum of drilling methods to rapidly provide hydraulic and chemistry data. A detailed distribution of anthropogenic and geogenic constituents within the saturated zones is produced, characterizing both permeable sediments and contaminant-laden boundaries between finer- and coarser-grained deposits. Applying this high-resolution data to well design leads to informed predictions on whether the new well will be compliant. Two case studies are presented, one for public supply well design and the other for managed aquifer recharge.