Factors Influencing the Stream‐Aquifer Flow Exchange Coefficient

Knowledge of river gain from or loss to a hydraulically connected water table aquifer is crucial in issues of water rights and also when attempting to optimize conjunctive use of surface and ground waters. Typically in groundwater models this exchange flow is related to a difference in head between the river and some point in the aquifer, through a “coefficient.” This coefficient has been defined differently as well as the location for the head in the aquifer. This paper proposes a new coefficient, analytically derived, and a specific location for the point where the aquifer head is used in the difference. The dimensionless part of the coefficient is referred to as the SAFE (stream‐aquifer flow exchange) dimensionless conductance. The paper investigates the factors that influence the value of this new conductance. Among these factors are (1) the wetted perimeter of the cross‐section, (2) the degree of penetration of the cross‐section, and (3) the shape of the cross‐section. The study shows that these factors just listed are indeed ordered in their respective level of importance. In addition the study verifies that the analytical correct value of the coefficient is matched by finite difference simulation only if the grid system is sufficiently fine. Thus the use of the analytical value of the coefficient is an accurate and efficient alternative to ad hoc estimates for the coefficient typically used in finite difference and finite element methods.     

Impacts of Seasonal Pumping on Stream‐Aquifer Interactions in Miryang,Korea

Large agricultural fields in South Korea are located mostly on alluvial plains, where a significant amount of groundwater is used for heating of water‐curtain insulated greenhouses. Such greenhouses are commonly used for crop cultivation during the winter dry season from November to March. After use the groundwater is discharged directly into streams, causing groundwater depletion. A hydrogeological study was carried out in a typical agricultural area of this type, located on an alluvial aquifer near the Nakdong River. Groundwater levels, chemical characteristics, and temperatures from 68 observation wells were analyzed to determine the impacts of seasonal groundwater pumping on the groundwater system and stream‐aquifer interactions. Our results show that the groundwater system has not yet reached a state of dynamic equilibrium. Decades of excessive seasonal pumping have caused a gradual decline of groundwater levels, leading to groundwater depletion, especially in areas further from the river. Seasonal pumping has also significantly affected groundwater quality in the aquifer near the river. Groundwater temperature is decreasing (in this case a disadvantage), and saline groundwater is being diluted by induced recharge. The results of this study provide a basic outline for effective integrated water management that is widely applicable in South Korea.     

Spatio‐Temporal Variations in Stream–Aquifer Interactions Following Construction of Weirs in Korea

The “Four Major Rivers Restoration Project” was conducted to secure sufficient water resources, introduce comprehensive flood control measures, improve water quality, and restore river ecosystems in Korea. As a part of the project, 16 sites were dredged and weirs were installed in the Han, Geum, Yeongsan, and Nakdong Rivers from late 2010 to early 2012. Groundwater data were obtained from 213 groundwater monitoring wells near the four major rivers to analyze the impacts of weir construction on the nearby groundwater flow system. The groundwater level and chemical characteristics were analyzed to investigate how the groundwater flow system and water quality changed following weir construction. Our results show that the groundwater level immediately increased with increased river levels following weir construction. In addition, the hydrologic condition of some rivers upstream of the weirs was changed from gaining to losing streams. Consequently, the direction of groundwater flow changed from perpendicular to parallel to the river, and groundwater downstream of the weir became recharged from the area upstream of the weir. This should affect groundwater quality, which should become similar to the river water; however, this change has not yet been observed. Therefore, both further monitoring of the groundwater quality and further hydrogeochemical analysis are required for quantitative evaluation of the effects of weir construction in the study area.     

Regional Flow and Deformation Analysis of Basin‐Fill Aquifer Systems Using Stress‐Dependent Parameters

Changes in effective stress due to water pressure variations modify the intrinsic hydrodynamic properties of aquifers and aquitards. Overexploited groundwater systems, such as basins with heavy pumping, are subject to nonrecoverable modifications. This results in loss of permeability, porosity, and specific storage due to system consolidation. This paper presents (1) the analytical development of model functions relating effective stress to hydrodynamic parameters for aquifers and aquitards constituted of unconsolidated granular sediments, and (2) a modeling approach for the analysis of aquifer systems affected by effective stress variations, taking into account the aforementioned dependency. The stress‐dependent functions were fit to laboratory data, and used in the suggested modeling approach. Based on only few unknowns, this approach is computationally simple, efficiently captures the hydromechanical processes that are active in regional aquifer systems under stress, and readily provides an estimate of their consolidation.     

A Stream‐Based Methane Monitoring Approach for Evaluating Groundwater Impacts Associated with Unconventional Gas Development

Gaining streams can provide an integrated signal of relatively large groundwater capture areas. In contrast to the point‐specific nature of monitoring wells, gaining streams coalesce multiple flow paths. Impacts on groundwater quality from unconventional gas development may be evaluated at the watershed scale by the sampling of dissolved methane (CH₄) along such streams. This paper describes a method for using stream CH₄ concentrations, along with measurements of groundwater inflow and gas transfer velocity interpreted by 1‐D stream transport modeling, to determine groundwater methane fluxes. While dissolved ionic tracers remain in the stream for long distances, the persistence of methane is not well documented. To test this method and evaluate CH₄ persistence in a stream, a combined bromide (Br) and CH₄ tracer injection was conducted on Nine‐Mile Creek, a gaining stream in a gas development area in central Utah. A 35% gain in streamflow was determined from dilution of the Br tracer. The injected CH₄ resulted in a fivefold increase in stream CH₄ immediately below the injection site. CH₄ and δ¹³C_CH4 sampling showed it was not immediately lost to the atmosphere, but remained in the stream for more than 2000 m. A 1‐D stream transport model simulating the decline in CH₄ yielded an apparent gas transfer velocity of 4.5 m/d, describing the rate of loss to the atmosphere (possibly including some microbial consumption). The transport model was then calibrated to background stream CH₄ in Nine‐Mile Creek (prior to CH₄ injection) in order to evaluate groundwater CH₄ contributions. The total estimated CH₄ load discharging to the stream along the study reach was 190 g/d, although using geochemical fingerprinting to determine its source was beyond the scope of the current study. This demonstrates the utility of stream‐gas sampling as a reconnaissance tool for evaluating both natural and anthropogenic CH₄ leakage from gas reservoirs into groundwater and surface water.     

Hydrogeophysical Methods for Analyzing Aquifer Storage and Recovery Systems

Hydrogeophysical methods are presented that support the siting and monitoring of aquifer storage and recovery (ASR) systems. These methods are presented as numerical simulations in the context of a proposed ASR experiment in Kuwait, although the techniques are applicable to numerous ASR projects. Bulk geophysical properties are calculated directly from ASR flow and solute transport simulations using standard petrophysical relationships and are used to simulate the dynamic geophysical response to ASR. This strategy provides a quantitative framework for determining site‐specific geophysical methods and data acquisition geometries that can provide the most useful information about the ASR implementation. An axisymmetric, coupled fluid flow and solute transport model simulates injection, storage, and withdrawal of fresh water (salinity ～500 ppm) into the Dammam aquifer, a tertiary carbonate formation with native salinity approximately 6000 ppm. Sensitivity of the flow simulations to the correlation length of aquifer heterogeneity, aquifer dispersivity, and hydraulic permeability of the confining layer are investigated. The geophysical response using electrical resistivity, time‐domain electromagnetic (TEM), and seismic methods is computed at regular intervals during the ASR simulation to investigate the sensitivity of these different techniques to changes in subsurface properties. For the electrical and electromagnetic methods, fluid electric conductivity is derived from the modeled salinity and is combined with an assumed porosity model to compute a bulk electrical resistivity structure. The seismic response is computed from the porosity model and changes in effective stress due to fluid pressure variations during injection/recovery, while changes in fluid properties are introduced through Gassmann fluid substitution.     

In‐Situ Sonication for Enhanced Recovery of Aquifer Microbial Communities

Sampling methods for characterization of microbial communities in aquifers should target both suspended and attached microorganisms (biofilms). We investigated the effectiveness and reproducibility of low‐frequency (200 Hz) sonication pulses on improving extraction efficiency and quality of microorganisms from a petroleum‐contaminated aquifer in Studen (Switzerland). Sonication pulses at different power levels (0.65, 0.9, and 1.1 kW) were applied to three different groundwater monitoring wells. Groundwater samples extracted after each pulse were compared with background groundwater samples for cell and adenosine tri‐phosphate concentration. Turbidity values were obtained to assess the release of sediment fines and associated microorganisms. The bacterial community in extracted groundwater samples was analyzed by terminal‐restriction‐fragment‐length polymorphism and compared with communities obtained from background groundwater samples and from sediment cores. Sonication enhanced the extraction efficiency up to 13‐fold, with most of the biomass being associated with the sediment fines extracted with groundwater. Consecutive pulses at constant power were decreasingly effective, while pulses with higher power yielded the best results both in terms of extraction efficiency and quality. Our results indicate that low‐frequency sonication may be a viable and cost‐effective tool to improve the extraction of microorganisms from aquifers, taking advantage of existing groundwater monitoring wells.     

Geostatistical Borehole Image‐Based Mapping of Karst‐Carbonate Aquifer Pores

Quantification of the character and spatial distribution of porosity in carbonate aquifers is important as input into computer models used in the calculation of intrinsic permeability and for next‐generation, high‐resolution groundwater flow simulations. Digital, optical, borehole‐wall image data from three closely spaced boreholes in the karst‐carbonate Biscayne aquifer in southeastern Florida are used in geostatistical experiments to assess the capabilities of various methods to create realistic two‐dimensional models of vuggy megaporosity and matrix‐porosity distribution in the limestone that composes the aquifer. When the borehole image data alone were used as the model training image, multiple‐point geostatistics failed to detect the known spatial autocorrelation of vuggy megaporosity and matrix porosity among the three boreholes, which were only 10 m apart. Variogram analysis and subsequent Gaussian simulation produced results that showed a realistic conceptualization of horizontal continuity of strata dominated by vuggy megaporosity and matrix porosity among the three boreholes.     

Sampling the Chemistry of Shallow Aquifer Systems — A Case Study

Pumped waters from 14 Pennsylvania wells, located in shallow sandstone, siltstone and shale aquifers, were continuously monitored for dissolved oxygen (D. O.), nitrate (NO₃), pH, electrical conductivity (EC) and water temperature in a discharge manifold at the well head. The amount of pumping or purging required to stabilize these parameter readings varied by well site and parameter being analyzed. However, the purging required was generally greatest for D. O. and least for water temperature where: D. O. < NO₃ pH < EC < water temperature. Wells located near the siltstone-shale interface generally required far more purging than did wells located elsewhere. Although parameter stability was often achieved within purging one bore volume, the complexity, diversity, and variability in the data and these well-ground water systems, suggest that no single purging rule is appropriate. Instead, the extent of purging required before sampling these shallow aquifers should be determined by incorporating on-site monitoring of target or related parameters into the purging process.
From a sampling perspective, the relationship between NO₃ and D. O. concentrations during purging were analyzed relative to aquifer type. For most wells located in sandstone or siltstone, NO₃ concentrations remained relatively constant during purging irrespective of changes in D. O. For most wells located in shale, these two were positively and similarly correlated, suggesting that a general relationship exists.     

A Leaky Aquifer below Champlain Sea Clay: Closed‐Form Solutions for Natural Seepage

Closed‐form solutions are proposed for natural seepage in semiconfined (leaky) aquifers such as those existing below the massive Champlain Sea clay layers in the Saint‐Lawrence River Valley. The solutions are for an ideal horizontal leaky aquifer below an ideal aquitard that may have either a constant thickness and a constant hydraulic head at its surface, or a variable thickness and a variable hydraulic head at its surface. A few simplifying assumptions were needed to obtain the closed‐form solutions. These have been verified using a finite element method, which did not make any of the assumptions but gave an excellent agreement for hydraulic heads and groundwater velocities. For example, the difference between the two solutions was smaller than 1 mm for variations in the 5 to 8 m range for the hydraulic head in the semiconfined aquifer. Note that fitting the hydraulic head data of monitoring wells to the theoretical solutions gives only the ratio of the aquifer and aquitard hydraulic conductivities, a clear case of multiple solutions for an inverse problem. Consequently, field permeability tests in the aquitard and the aquifer, and pumping tests in the aquifer, are still needed to determine the hydraulic conductivity values.     

Mapping Potential Groundwater‐Dependent Ecosystems for Sustainable Management

Ecosystems which rely on either the surface expression or subsurface presence of groundwater are known as groundwater‐dependent ecosystems (GDEs). A comprehensive inventory of GDE locations at an appropriate management scale is a necessary first‐step for sustainable management of supporting aquifers; however, this information is unavailable for most areas of concern. To address this gap, this study created a two‐step algorithm which analyzed existing geospatial and remote sensing data to identify potential GDEs at both state/province and aquifer/basin scales. At the state/province scale, a geospatial information system (GIS) database was constructed for Texas, including climate, topography, hydrology, and ecology data. From these data, a GDE index was calculated, which combined vegetative and hydrological indicators. The results indicated that central Texas, particularly the Edwards Aquifer region, had highest potential to host GDEs. Next, an aquifer/basin scale remote sensing‐based algorithm was created to provide more detailed maps of GDEs in the Edwards Aquifer region. This algorithm used Landsat ETM+ and MODIS images to track the changes of NDVI for each vegetation pixel. The NDVI dynamics were used to identify the vegetation with high potential to use groundwater—such plants remain high NDVI during extended dry periods and also exhibit low seasonal and inter‐annual NDVI changes between dry and wet seasons/years. The results indicated that 8% of natural vegetation was very likely using groundwater. Of the potential GDEs identified, 75% were located on shallow soil averaging 45 cm in depth. The dominant GDE species were live oak, ashe juniper, and mesquite.     

Incorporating Super‐Diffusion due to Sub‐Grid Heterogeneity to Capture Non‐Fickian Transport

Numerical transport models based on the advection‐dispersion equation (ADE) are built on the assumption that sub‐grid cell transport is Fickian such that dispersive spreading around the average velocity is symmetric and without significant tailing on the front edge of a solute plume. However, anomalous diffusion in the form of super‐diffusion due to preferential pathways in an aquifer has been observed in field data, challenging the assumption of Fickian dispersion at the local scale. This study develops a fully Lagrangian method to simulate sub‐grid super‐diffusion in a multidimensional regional‐scale transport model by using a recent mathematical model allowing super‐diffusion along the flow direction given by the regional model. Here, the time randomizing procedure known as subordination is applied to flow field output from MODFLOW simulations. Numerical tests check the applicability of the novel method in mapping regional‐scale super‐diffusive transport conditioned on local properties of multidimensional heterogeneous media.     

Patterns of Entrapped Air Dissolution in a Two‐Dimensional Pilot‐Scale Synthetic Aquifer

Past studies of entrapped air dissolution have focused on one‐dimensional laboratory columns. Here the multidimensional nature of entrapped air dissolution was investigated using an indoor tank (180 × 240 × 600 cm³) simulating an unconfined sand aquifer with horizontal flow. Time domain reflectometry (TDR) probes directly measured entrapped air contents, while dissolved gas conditions were monitored with total dissolved gas pressure (P_TDG) probes. Dissolution occurred as a diffuse wedge‐shaped front from the inlet downgradient, with preferential dissolution at depth. This pattern was mainly attributed to increased gas solubility, as shown by P_TDG measurements. However, compression of entrapped air at greater depths, captured by TDR and leading to lower quasi‐saturated hydraulic conductivities and thus greater velocities, also played a small role. Linear propagation of the dissolution front downgradient was observed at each depth, with both TDR and P_TDG, with increasing rates with depth (e.g, 4.1 to 5.7× slower at 15 cm vs. 165 cm depth). P_TDG values revealed equilibrium with the entrapped gas initially, being higher at greater depth and fluctuating with the barometric pressure, before declining concurrently with entrapped air contents to the lower P_TDG of the source water. The observed dissolution pattern has long‐term implications for a wide variety of groundwater management issues, from recharge to contaminant transport and remediation strategies, due to the persistence of entrapped air near the water table (potential timescale of years). This study also demonstrated the utility of P_TDG probes for simple in situ measurements to detect entrapped air and monitor its dissolution.