Storm-damaged saline-contaminated boreholes as a means of aquifer contamination

Saline water from a storm surge can flow down storm-damaged submerged water supply wells and contaminate boreholes and surrounding aquifers. Using data from conventional purging techniques, aquifer test response analysis, chemical analysis, and regression analysis of chloride/silica (Cl/Si) ratio, equations were derived to estimate the volume of saline water intrusion into a well and a porous media aquifer, the volume of water needed to purge a well shortly following an intrusion event, and the volume of water needed after delay of several or more months, when the saline plume has expanded. Purging time required is a function of volume of water and pumping rate. The study site well is located within a shoreline community of Lake Pontchartrain, St. Tammany Parish, in southeastern Louisiana, United States, which was impacted by two hurricane storm surges and had neither been rehabilitated nor chlorinated prior to our study. Chemical analysis of water samples in fall 2005 and purging of well and aquifer in June 6, 2006, indicated saline water had intruded the well in 2005 and the well and aquifer in 2006. The volume of water needed to purge the study well was approximately 200 casing volumes, which is significantly greater than conventionally used during collection of water samples for water quality analyses.     

Evaluation of Bias Associated with Capture Maps Derived from Nonlinear Groundwater Flow Models

The impact of groundwater withdrawal on surface water is a concern of water users and water managers, particularly in the arid western United States. Capture maps are useful tools to spatially assess the impact of groundwater pumping on water sources (e.g., streamflow depletion) and are being used more frequently for conjunctive management of surface water and groundwater. Capture maps have been derived using linear groundwater flow models and rely on the principle of superposition to demonstrate the effects of pumping in various locations on resources of interest. However, nonlinear models are often necessary to simulate head‐dependent boundary conditions and unconfined aquifers. Capture maps developed using nonlinear models with the principle of superposition may over‐ or underestimate capture magnitude and spatial extent. This paper presents new methods for generating capture difference maps, which assess spatial effects of model nonlinearity on capture fraction sensitivity to pumping rate, and for calculating the bias associated with capture maps. The sensitivity of capture map bias to selected parameters related to model design and conceptualization for the arid western United States is explored. This study finds that the simulation of stream continuity, pumping rates, stream incision, well proximity to capture sources, aquifer hydraulic conductivity, and groundwater evapotranspiration extinction depth substantially affect capture map bias. Capture difference maps demonstrate that regions with large capture fraction differences are indicative of greater potential capture map bias. Understanding both spatial and temporal bias in capture maps derived from nonlinear groundwater flow models improves their utility and defensibility as conjunctive‐use management tools.     

Boundary condition effects on maximum groundwater withdrawal in coastal aquifers

Prevention of sea water intrusion in coastal aquifers subject to groundwater withdrawal requires optimization of well pumping rates to maximize the water supply while avoiding sea water intrusion. Boundary conditions and the aquifer domain size have significant influences on simulating flow and concentration fields and estimating maximum pumping rates. In this study, an analytical solution is derived based on the potential-flow theory for evaluating maximum groundwater pumping rates in a domain with a constant hydraulic head landward boundary. An empirical correction factor, which was introduced by Pool and Carrera (2011) to account for mixing in the case with a constant recharge rate boundary condition, is found also applicable for the case with a constant hydraulic head boundary condition, and therefore greatly improves the usefulness of the sharp-interface analytical solution. Comparing with the solution for a constant recharge rate boundary, we find that a constant hydraulic head boundary often yields larger estimations of the maximum pumping rate and when the domain size is five times greater than the distance between the well and the coastline, the effect of setting different landward boundary conditions becomes insignificant with a relative difference between two solutions less than 2.5%. These findings can serve as a preliminary guidance for conducting numerical simulations and designing tank-scale laboratory experiments for studying groundwater withdrawal problems in coastal aquifers with minimized boundary condition effects.     

Concentration Profiles in Screened Wells under Static and Pumped Conditions

Purge and pump samples from screened wells reflect concentration averaging and contaminant redistribution by wellbore flow. These issues were assessed in a screened well at the Hanford Site by investigating the vertical profile of a technetium-99 plume in a conventional well under static and pumped conditions. Specific conductance and technetium-99 concentrations were well correlated, and this enabled measurement of specific conductance to be used as a surrogate for technetium-99 concentration. Time-series measurements were collected during purging from three specific conductance probes installed in the well at 1.2, 3.1, and 4.9 m below the static water level in a 7.7-m-deep screened well. The vertical contaminant profile adjacent to the well in the aquifer was calculated using the concentration profile in the well during pumping, the pumping flow rate, and a wellbore flow and mixing model. The plume was found to be stratified in the aquifer—the highest concentrations occurred adjacent to the upper part of the screened interval. The purge and pump sample concentrations were 41% to 58% of the calculated peak concentration in the aquifer. Plume stratification in the aquifer adjacent to the well screen became more pronounced as pumping continued. Extended pumping may have partially reversed the effect of contaminant redistribution in the aquifer by wellbore flow and allowed the stratification of the plume to be more observable. It was also found that the vertical profile of contamination in the well under static (i.e., nonpumping conditions) was not representative of the profile in the aquifer. Thus, passive or micropurge sampling techniques, which sample the wellbore water at different depths, would not yield results representative of the aquifer in this well.     

Estimating Flow Using Tracers and Hydraulics in Synthetic Heterogeneous Aquifers

Regional ground water flow is most usually estimated using Darcy's law, with hydraulic conductivities estimated from pumping tests, but can also be estimated using ground water residence times derived from radioactive tracers. The two methods agree reasonably well in relatively homogeneous aquifers but it is not clear which is likely to produce more reliable estimates of ground water flow rates in heterogeneous systems. The aim of this paper is to compare bias and uncertainty of tracer and hydraulic approaches to assess ground water flow in heterogeneous aquifers. Synthetic two-dimensional aquifers with different levels of heterogeneity (correlation lengths, variances) are used to simulate ground water flow, pumping tests, and transport of radioactive tracers. Results show that bias and uncertainty of flow rates increase with the variance of the hydraulic conductivity for both methods. The bias resulting from the nonlinearity of the concentration–time relationship can be reduced by choosing a tracer with a decay rate similar to the mean ground water residence time. The bias on flow rates estimated from pumping tests is reduced when performing long duration tests. The uncertainty on ground water flow is minimized when the sampling volume is large compared to the correlation length. For tracers, the uncertainty is related to the ratio of correlation length to the distance between sampling wells. For pumping tests, it is related to the ratio of correlation length to the pumping test's radius of influence. In regional systems, it may be easier to minimize this ratio for tracers than for pumping tests.     

The impact of well drawdowns on the mixing process of river water and groundwater and water quality in a riverside well field,Northeast China

Riverbank filtration (RBF) has been widely used throughout the world as an effective means to regulate surface water and groundwater resources and pretreat raw water for municipal water supply. The quality of the water from a riverside well field and the mixing ratios of surface water and groundwater is primarily impacted by the hydrodynamic processes in the RBF system. The RBF system is largely controlled by the water exploiting system in addition to the natural hydrologic condition of the river–aquifer system. As one of the most important design parameters of the riverside well field, the drawdown of groundwater level greatly determines the water head differences between the river water and groundwater as well as the field flow net, which subsequently impacts the mixing of river water and groundwater and water quality significantly. This study aimed to improve the understanding of the mixing process between the surface water and groundwater and estimate the impact of the RBF on the mixing ratio of surface water–groundwater and water quality quantitatively. A set of field pumping tests with various groundwater level drawdowns were carried out independently and successively at a riverside field with a single pumping well near the Songhua River in Northeast China in August 2017. During these tests, the water levels and hydrochemical parameters of the Songhua River, the adjacent aquifer, and the pumping well were monitored. The dynamic mixing process of the river water and groundwater and water quality under various drawdown conditions were analysed systematically using analytical methods. The results obtained from Dupuit method and the Mirror Image method in conjunction with the Hydrochemical Tracing method showed that the pumping water directly from the river water reached 60% ± 10% after a steady flow net was established. The larger the proportion of the pumping water captured from the river, the better quality of the pumping water was, because the quality of the river water (only limited to some water quality parameters monitored which were minority) was better than that of the groundwater. The results also showed that total Fe, TDS, total hardness, COD_Mn, and K⁺ were relatively sensitive to the changes of groundwater drawdown, and their concentrations decreased with an increase in the groundwater drawdown. It can be concluded that both the mixing ratio of the surface water and the groundwater and the water quality of the riverside well field can be regulated through adjusting the designed drawdown of the groundwater level, which is helpful for the design and the optimization of the riverside well water intake engineering.     

Anthropogenic influences on groundwater in the vicinity of a long‐lived radioactive waste repository

The groundwater flow system in the Culebra Dolomite Member (Culebra) of the Permian Rustler Formation is a potential radionuclide release pathway from the Waste Isolation Pilot Plant (WIPP), the only deep geological repository for transuranic waste in the United States. In early conceptual models of the Culebra, groundwater levels were not expected to fluctuate markedly, except in response to long‐term climatic changes, with response times on the order of hundreds to thousands of years. Recent groundwater pressures measured in monitoring wells record more than 25 m of drawdown. The fluctuations are attributed to pumping activities at a privately owned well that may be associated with the demand of the Permian Basin hydrocarbon industry for water. The unprecedented magnitude of drawdown provides an opportunity to quantitatively assess the influence of unplanned anthropogenic forcings near the WIPP. Spatially variable realizations of Culebra saturated hydraulic conductivity and specific storage were used to develop groundwater flow models to estimate a pumping rate for the private well and investigate its effect on advective transport. Simulated drawdown shows reasonable agreement with observations (average Model Efficiency coefficient = 0.7). Steepened hydraulic gradients associated with the pumping reduce estimates of conservative particle travel times across the domain by one half and shift the intersection of the average particle track with the compliance boundary by more than 2 km. The value of the transient simulations conducted for this study lies in their ability to (a) improve understanding of the Culebra groundwater flow system and (b) challenge the notion of time‐invariant land use in the vicinity of the WIPP.     

Monitoring of drawdown pattern during pumping in an unconfined physical aquifer model

In this study, we attempted to analyse a drawdown pattern around a pumping well in an unconfined sandy gravelly aquifer constructed in a laboratory tank by means of both experimental and numerical modelling of groundwater flow. The physical model consisted of recharge, aquifer and discharge zones. Permeability and specific yield of the aquifer material were determined by Dupuit approximation under steady‐state flow and stepwise gravitational drainage of groundwater, respectively. The drawdown of water table in pumping and neighbouring observation wells was monitored to investigate the effect of no‐flow boundary on the drawdown pattern during pumping for three different boundary conditions: (i) no recharge and no discharge with four no‐flow boundaries (Case 1); (ii) no recharge and reservoir with three no‐flow boundaries (Case 2); (iii) recharge and discharge with two no‐flow boundaries (Case 3). Based on the aquifer parameters, numerical modelling was also performed to compare the simulated drawdown with that observed. Results showed that a large difference existed between the simulated drawdown and that observed in wells for all cases. The reason for the difference could be explained by the formation of a curvilinear type water table between wells rather than a linear one due to a delayed response of water table in the capillary fringe. This phenomenon was also investigated from a mass balance study on the pumping volume. The curvilinear type of water table was further evidenced by measurement of water contents at several positions in the aquifer between wells using time domain reflectometry (TDR). This indicates that the existing groundwater flow model applicable to an unconfined aquifer lacks the capacity to describe a slow response of water table in the aquifer and care should be taken in the interpretation of water table formation in the aquifer during pumping. Copyright © 2001 John Wiley & Sons, Ltd.     

Two-stage aquifer pumping subject to slow desorption and persistent sources

This work focuses on improving pump-and-treat remediation by optimizing a two-stage operational scheme to reduce volumes extracted when confronted with nonequilibrium desorption, low-permeability units, and continuous contaminant sources such as non-aqueous phase liquids (NAPL). Q1 and Q2 are the initial short-term high pumping rate and later long-term low pumping rate, respectively. A two-dimensional ground water flow and transport management model was used to test the proposed strategy for plumes developed from finite (NAPL-free) and continuous (NAPL-driven) contaminant sources in homogeneous and nonhomogeneous (zoned) aquifers. Remediation scenarios were simulated over durations of 2000, 6000, and 15,000 d to determine (1) the optimal time to switch from a preset Q1 to Q2 and (2) the value of Q2. The problem was constrained by mass removal requirements, maximum allowable downgradient concentrations, and practical bounds on Q2. Q1 was fixed at preset values 50% to 200% higher than the single-stage pumping rates (i.e., steady pumping rates during entire remediation period) necessary to achieve a desired cleanup level and capture the plume. Results for the NAPL-free homogeneous case under nonequilibrium desorption conditions achieved the same level of cleanup as single-stage pumping, while reducing extracted volumes by up to 36%. Comparable savings were obtained with NAPL-driven sources only when the source concentration was reduced by at least 2 orders of magnitude. For the zoned aquifer, the proposed strategy provided volume savings of up to 24% under NAPL-free and reduced source conditions.     

Estimation of Salt Water Upconing Using a Steady‐State Solution for Partial Completion of a Pumped Well

A new steady‐state analytical solution to the two‐dimensional radial‐flow equation was developed for drawdown (head) conditions in an aquifer with constant transmissivity, no‐flow conditions at the top and bottom, constant head conditions at a known radial distance, and a partially completed pumping well. The solution was evaluated for accuracy by comparison to numerical simulations using MODFLOW. The solution was then used to estimate the rise of the salt water‐fresh water interface (upconing) that occurs under a pumping well, and to calculate the critical pumping rate at which the interface becomes unstable, allowing salt water to enter the pumping well. The analysis of salt water‐fresh water interface rise assumed no significant effect on upconing by recharge; this assumption was tested and supported using results from a new steady‐state analytical solution developed for recharge under two‐dimensional radial‐flow conditions. The upconing analysis results were evaluated for accuracy by comparison to those from numerical simulations using SEAWAT for salt water‐fresh water interface positions under mild pumping conditions. The results from the equation were also compared with those of a published numerical sharp‐interface model applied to a case on Cape Cod, Massachusetts. This comparison indicates that estimating the interface rise and maximum allowable pumping rate using the analytical method will likely be less conservative than the maximum allowable pumping rate and maximum stable interface rise from a numerical sharp‐interface model.     

Fracture Flow Characterization with Low-Noise Spontaneous Potential Logging

Geophysical well logging has been applied for fracture characterization in crystalline terrains by physical properties measurements and borehole wall imaging. Some of these methods can be applied to monitor pumping tests to identify fractures contributing to groundwater flow and, with this, determine hydraulic conductivity and transmissivity along the well. We present a procedure to identify fractures contributing to groundwater flow using spontaneous potential measurements generated by electrokinetic processes when the borehole water head is lowered and then monitored while recovering. The electrokinetic model for flow through a tabular gap is used to interpret the measured data and determine the water head difference that drives the flow through the fracture. We present preliminary results at a test site in crystalline rocks on the campus of the University of São Paulo.     

Quantifying changes of effective springshed area and net recharge through recession analysis of spring flow

The average flow of Silver Springs, one of the largest magnitude springs in Central Florida, declined 32% from 2000 to 2012. The average groundwater head in the springshed declined 0.14 m, and the spring pool altitude increased 0.24 m during the same period. This paper presents a novel explanation of the spring flow recession curve for Silver Springs using the Torricelli model, which uses the groundwater head at a sentinel well, the spring pool altitude and the net recharge to groundwater. The effective springshed area and net recharge (defined as recharge minus groundwater pumping and evapotranspiration) were estimated based on the observed recession slopes for spring flow, groundwater head and spring pool altitude. The results indicate that the effective springshed area continuously declined since 1989 and the net recharge declined since the 1970s with a significant drop in 2002. Subsequent to 2002, the net recharge increased modestly but not to the levels prior to the 1990s. The reduction in net recharge was caused by changes in hydroclimatic conditions including precipitation and air temperature, along with groundwater withdrawals, which contributed to the declined spring flow. Copyright © 2016 John Wiley & Sons, Ltd.     

A local model for analysis of pump and treat systems with vertical barrier walls

We present a mathematical model of local, steady groundwater flow near a vertical barrier wall. Flow features represented in the model include an impermeable arc-shaped barrier wall and multiple wells; distant boundary conditions are not included explicitly, but their effects on the local flow field are modelled by specifying a uniform flow at infinity and a constant areal recharge within a local domain. We develop an explicit closed-form solution to the boundary-value problem using the analytic element method. The solution is an extension of a harmonic solution presented by Anderson and Mesa [Anderson EI, Mesa E. The effects of vertical barrier walls on the hydraulic control of contaminated groundwater. Adv Water Resourc 2006;29(1):89–98] which does not include the effects of recharge. We demonstrate that the general solution with recharge consists of the harmonic solution superposed on a special case of the harmonic solution along with two elementary one-dimensional flow solutions. The results are used to investigate the effects of areal recharge on the capture zone envelopes of the pumping wells and on the reduction in discharge that can be achieved by including a barrier wall in a pump and treat design. We find that the benefits of including an open barrier wall in a design, measured as a reduction in the pumping rate required to contain a plume, increase for higher recharge rates. Dimensionless plots of capture zone envelopes are presented for a practical well and barrier wall configuration.     

Influence of a thin veneer of low‐hydraulic‐conductivity sediment on modelled exchange between river water and groundwater in response to induced infiltration

A thin layer of fine‐grained sediment commonly is deposited at the sediment–water interface of streams and rivers during low‐flow conditions, and may hinder exchange at the sediment–water interface similar to that observed at many riverbank‐filtration (RBF) sites. Results from a numerical groundwater‐flow model indicate that a low‐permeability veneer reduces the contribution of river water to a pumping well in a riparian aquifer to various degrees, depending on simulated hydraulic gradients, hydrogeological properties, and pumping conditions. Seepage of river water is reduced by 5–10% when a 2‐cm thick, low‐permeability veneer is present on the bed surface. Increasing thickness of the low‐permeability layer to 0·1 m has little effect on distribution of seepage or percentage contribution from the river to the pumping well. A three‐orders‐of‐magnitude reduction in hydraulic conductivity of the veneer is required to reduce seepage from the river to the extent typically associated with clogging at RBF sites. This degree of reduction is much larger than field‐measured values that were on the order of a factor of 20–25. Over 90% of seepage occurs within 12 m of the shoreline closest to the pumping well for most simulations. Virtually no seepage occurs through the thalweg near the shoreline opposite the pumping well, although no low‐permeability sediment was simulated for the thalweg. These results are relevant to natural settings that favour formation of a substantial, low‐permeability sediment veneer, as well as central‐pivot irrigation systems, and municipal water supplies where river seepage is induced via pumping wells. Published in 2011 by John Wiley & Sons, Ltd.     

Spatial uncertainty in bias corrected climate change projections and hydrogeological impacts

The question of which climate model bias correction methods and spatial scales for correction are optimal for both projecting future hydrological changes as well as removing initial model bias has so far received little attention. For 11 climate models (CMs), or GCM/RCM – Global/Regional Climate Model pairing, this paper analyses the relationship between complexity and robustness of three distribution‐based scaling (DBS) bias correction methods applied to daily precipitation at various spatial scales. Hydrological simulations are forced by CM inputs to assess the spatial uncertainty of groundwater head and stream discharge given the various DBS methods. A unique metric is devised, which allows for comparison of spatial variability in climate model bias and projected change in precipitation. It is found that the spatial variability in climate model bias is larger than in the climate change signals. The magnitude of spatial bias seen in precipitation inputs does not necessarily correspond to the magnitude of biases seen in hydrological outputs. Variables that integrate basin responses over time and space are more sensitive to mean spatial biases and less so on extremes. Hydrological simulations forced by the least parameterized DBS approach show the highest error in mean and maximum groundwater heads; however, the most highly parameterised DBS approach shows less robustness in future periods compared with the reference period it was trained in. For hydrological impacts studies, choice of bias correction method should depend on the spatial scale at which hydrological impacts variables are required and whether CM initial bias is spatially uniform or spatially varying. Copyright © 2015 John Wiley & Sons, Ltd.     

测试水井导水性能的简便试验

这种测试水井含水层导水系数的新方法是利用井口水面空气压力阶变试验,计算出聊城井含水层的导水系数为14.4m2/d,而通过抽水试验得出的导水系数为43.6m2/d,二者有相同的数量级,但新方法比抽水试验简单易行和经济快速,用它来解释地下水微动态更合理,容易推广应用     

Analytical Solutions of Travel Time to a Pumping Well with Variable Evapotranspiration

Analytical solutions of groundwater travel time to a pumping well in an unconfined aquifer have been developed in previous studies, however, the change in evapotranspiration was not considered. Here, we develop a mathematical model of unconfined flow toward a discharge well with redistribution of groundwater evapotranspiration for travel time analysis. Dependency of groundwater evapotranspiration on the depth to water table is described using a linear formula with an extinction depth. Analytical solutions of groundwater level and travel time are obtained. For a typical hypothetical example, these solutions perfectly agree with the numerical simulation results based on MODFLOW and MODPATH. As indicated in a dimensionless framework, a lumped parameter which is proportional to the pumping rate controls the distributions of groundwater evapotranspiration rate and the travel time along the radial direction.     

Identifying groundwater pumping source information using simulated annealing

This study develops an approach referred to as SA‐MF using simulated annealing (SA) and a three‐dimensional groundwater flow model (MODFLOW‐2000) to determine the pumping source information, namely the pumping source location, the pumping rate and the pumping period. Eight scenarios, each with various cases for aquifers with homogeneous or heterogeneous conductivities, are used to test the applicability and accuracy of SA‐MF in determining the pumping source information. The results show that at least five measured heads should be used to analyse the pumping source identification problem and at least three observation wells are required to effectively determine the pumping source information. The SA‐MF gives good estimated results in both synthetic and real problems even if the measured heads contain measurement errors. Copyright © 2008 John Wiley & Sons, Ltd.     

Tritium in ground water of the georgia piedmont: Implications for recharge and flow paths

Environmental tritium was measured in 33 natural water samples representative of precipitation, stream runoff, and groundwater (derived principally from production wells) within the Georgia Piedmont Province. Major ion analyses were used to assist in the interpretation of the tritium results. Tritium concentrations were significantly greater within shallow groundwater derived from the regolith (28–34 TU) and stream runoff (25–30 TU) than within recent rainfall (4–17 TU). Based upon the decay-corrected tritium input function, this probably indicates that at least some of the shallow water is stored within the regolith for a period of approximately 25 years. A ‘post-bomb’ component of recharge was present in all groundwater derived from production wells in the study area. Groundwater sampled from the bedrock aquifers was commonly less tritiated than either stream runoff or shallow water stored in the regolith. the lower tritium concentrations May, have resulted from the mixing of ‘pre-bomb’ water stored within the fractures or the transitional zone directly above the bedrock and modern water stored in the shallow regolith. the preponderance of modern water provides evidence that groundwater flow paths are areally restricted within this setting, probably confined to local surface water drainage basins. the residence time of groundwater in the Piedmont is limited by the lack of deep, gravity-driven regional flow and the localized vertical flow induced by pumping. the results of this study indicate that relatively small tritium concentration variations (10-20 TU) May, have regional hydrological significance in the southeastern Piedmont Province and similar settings.     

The Influence of Clay Zones on Land Subsidence from Groundwater Pumping

The objective of this article is to analyze the influence of clay zones on subsidence from groundwater pumping. Finite element analyses were conducted on a sand‐only aquifer and a sand aquifer with two clay zones located at different distances from the well face. A model that accounts for recoverable and nonrecoverable strains was used to simulate the sand and clay. This model couples the groundwater flow with the stress‐deformation response of the aquifer materials. Each aquifer was pumped from a single well for a period of 6 months, and then the groundwater level was lowered gradually to an elevation below the elevation of the clay zones and kept there for 10 years. The groundwater level was then raised gradually back to the original elevation over a period of 10 years. The results of the analyses show that the ground surface subsidence profile is strongly influenced by the presence of the clays zones. The ground surface sags where these clay zones are present resulting in a wavy ground surface profile. Subsidence continued when pumping is stopped, albeit at a much slower rate than during pumping, and when the groundwater level is below the elevation of the clay zones. Clay zones further away from the well face lag the subsidence of clay zones nearer the well face because of lower changes in hydrostatic head. Sags in ground surface subsidence profile from groundwater pumping are indicators of the presence of low hydraulic conductive geological materials.