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.     

Factors Governing Sustainable Groundwater Pumping near a River

The objective of this paper was to provide new insights into processes affecting riverbank filtration (RBF). We consider a system with an inflatable dam installed for enhancing water production from downstream collector wells. Using a numerical model, we investigate the impact of groundwater pumping and dam operation on the hydrodynamics in the aquifer and water production. We focus our study on two processes that potentially limit water production of an RBF system: the development of an unsaturated zone and riverbed clogging. We quantify river clogging by calibrating a time‐dependent riverbed permeability function based on knowledge of pumping rate, river stage, and temperature. The dynamics of the estimated riverbed permeability reflects clogging and scouring mechanisms. Our results indicate that (1) riverbed permeability is the dominant factor affecting infiltration needed for sustainable RBF production; (2) dam operation can influence pumping efficiency and prevent the development of an unsaturated zone beneath the riverbed only under conditions of sufficient riverbed permeability; (3) slow river velocity, caused by dam raising during summer months, may lead to sedimentation and deposition of fine‐grained material within the riverbed, which may clog the riverbed, limiting recharge to the collector wells and contributing to the development of an unsaturated zone beneath the riverbed; and (4) higher river flow velocities, caused by dam lowering during winter storms, scour the riverbed and thus increase its permeability. These insights can be used as the basis for developing sustainable water management of a RBF system.     

Analytical Solutions for Determination of Non-Steady- State and Steady-State Capture Zones

This paper presents analytical solutions for determining non-steady-state capture zones produced by a single recovery well and steady-state capture zones produced by multiple recovery wells. Analysis of non-steady-slate capture zones is based on the lime-dependent location of caplure zone stagnation points and the geometric similarity between steady-slate and non-steady-state capture zones. The analytical solution of steady-state capture zones is obtained from spatial variations of discharge potential across the capture zone boundary. Both capture zone analyses are based on the assumptions of uniform flow field with a constant hydraulic conductivity, the Dupuit assumption of insignificant vertical flow, a negligible delayed yield, and a fully penetrating well with a constant pumping rate. For a ground water pump-and-trcat remediation program, the pumping rate and well location design variables can be adjusted to ensure containment of the ground water contaminant plume.     

Analysis of steady ground water flow toward wells in a confined-unconfined aquifer

A confined aquifer may become unconfined near the pumping wells when the water level falls below the confining unit in the case where the pumping rate is great and the excess hydraulic head over the top of the aquifer is small. Girinskii's potential function is applied to analyze the steady ground water flow induced by pumping wells with a constant-head boundary in a mixed confined-unconfined aquifer. The solution of the single-well problem is derived, and the critical radial distance at which the flow changes from confined to unconfined condition is obtained. Using image wells and the superposition method, an analytic solution is presented to study steady ground water flow induced by a group of pumping wells in an aquifer bounded by a river with constant head. A dimensionless function is introduced to determine whether a water table condition exists or not near the pumping wells. An example with three pumping wells is used to demonstrate the patterns of potentiometric surface and development of water table around the wells.     

Interaction of Aquifer and River‐Canal Network near Well Field

The article presents semi‐analytical mathematical models to asses (1) enhancements of seepage from a canal and (2) induced flow from a partially penetrating river in an unconfined aquifer consequent to groundwater withdrawal in a well field in the vicinity of the river and canal. The nonlinear exponential relation between seepage from a canal reach and hydraulic head in the aquifer beneath the canal reach is used for quantifying seepage from the canal reach. Hantush's (1967) basic solution for water table rise due to recharge from a rectangular spreading basin in absence of pumping well is used for generating unit pulse response function coefficients for water table rise in the aquifer. Duhamel's convolution theory and method of superposition are applied to obtain water table position due to pumping and recharge from different canal reaches. Hunt's (1999) basic solution for river depletion due to constant pumping from a well in the vicinity of a partially penetrating river is used to generate unit pulse response function coefficients. Applying convolution technique and superposition, treating the recharge from canal reaches as recharge through conceptual injection wells, river depletion consequent to variable pumping and recharge is quantified. The integrated model is applied to a case study in Haridwar (India). The well field consists of 22 pumping wells located in the vicinity of a perennial river and a canal network. The river bank filtrate portion consequent to pumping is quantified.     

Modeling of Contaminant Movement Near Pumping Wells: Saturated-Unsaturated Flow with Particle Tracking

A transient axisymmetric saturated-unsaturated numerical flow model was coupled with a particle tracking model to investigate the movement of contaminants when a shallow unconfined aquifer is pumped at a constant rate. The particle tracking model keeps track of locations and masses of solutes in the aquifer, and the time of capture by the well. At the end of each time-step the flow model solves the Richard's equation for the hydraulic head distribution from which elemental velocities are calculated. Solutes are then displaced for a period equivalent to the time-step using both the magnitude and direction of the elemental velocities. Numerical experiments were performed to investigate effluent concentrations in wells with screens of different length and in different positions relative to zones of stratified contamination. At early times of pumping the effluent concentrations were similar to the concentrations adjacent to the well screen, but at late times, the concentrations approached the vertically averaged concentration in the aquifer. Time to attain the vertically averaged concentration was determined by the well geometry, initial location of the contaminant plume in relation to the well screen, and hydraulic properties of the aquifer. The results are consistent with the hydraulics of flow to a pumping well and of particular importance, they demonstrate that short-term pump tests could give erroneous design concentrations for pump-and-treat systems. The model provides a means of quantifying arrival times and mixing ratios. It could therefore provide a useful means of designing production wells in aquifers with stratified contamination and more efficient recovery systems for aquifer remediation.     

Groundwater flow path determination during riverbank filtration affected by groundwater exploitation: a case study of Liao River,Northeast China

The groundwater flow path plays an important role in maintaining hydrological and ecological quality and security, which are important in the comprehensive management and use of both groundwater and surface water. In this study, an integrated multi-tracer-constrained framework was used to determine the groundwater flow path. The results show that there are shallow and deep flow paths in riverbank filtration, controlled by the different permeabilities of riverbed sediments and aquifers at different depths. The contribution of river water to shallow groundwater is less than that to deep groundwater because of the low permeability of the riverbed sediment in the dense muddy layer in the shallow slope of the river valley. This contribution decreases with increasing distance from the Liao River. The shallow groundwater quality is better than the deep groundwater quality because of its longer residence time.     

Modeling Stream-Aquifer Interactions Under Seasonal Groundwater Pumping and Managed Aquifer Recharge

In South Korea, a significant amount of groundwater is used for the heating of water-curtain insulated greenhouses during the winter dry season, which had led to problems of groundwater depletion. A managed aquifer recharge (MAR) project is currently underway with the goal of preventing such groundwater depletion in a typical cultivation area, located on an alluvial aquifer near the Nam River. In the present study, FEFLOW, a three-dimensional finite element model, was used to evaluate different strategies for MAR of the cultivation areas. A conceptual model was developed to simulate the stream-aquifer dynamics under the influence of seasonal groundwater pumping and MAR. The optimal rates and duration of MAR were assessed by analyzing the recovery of the groundwater levels and the change in the groundwater temperature. The simulation results indicate that a MAR rate of 8000 m³/d effectively restores the groundwater level when the injection wells are located inside the groundwater depletion area. It is also demonstrated that starting the MAR before the beginning of the seasonal pumping is more effective. Riverbank filtration is preferable for securing the injection water owing to plentiful source of induced recharge from the river. Locating the pumping wells adjacent to the river where there are thick permeable layers could be a good strategy for minimizing decreases in the groundwater level and temperature.     

Vertical Contaminant Profiling of Volatile Organic* in a Deep Fractured Basalt Aquifer

Volatile organic compounds delected in ground water from wells at Test Area North (TAN) at the Idaho National Engineering Laboratory (INEL) prompted RCRA facility investigations in 1989 and 1990 and a CERCLA-driven RI/FS in 1992. In order to address ground water treatment feasibility, one of the main objectives, of the 1992 remedial investigation was to determine the vertical extent of ground water contamination, where the principle contaminant, of concern is trichloroethylene (TCE). It was hypothesized that a sedimentary interbed at depth in the fractured basalt aquifer could be inhibiting vertical migration of contaminants to lower aquifers. Due to the high cost of drilling and installation of ground water monitoring wells at this facility (greater than $100,000 per well), a real time method was proposed for obtaining and analyzing ground water samples during drilling to allow accurate placement of well screens in zones of predicted VOC contamination. This method utilized an inflatable pump packer pressure transducer system interfaced with a datalogger and PC at land surface. This arrangement allowed for real lime monitoring of hydraulic head above and below the packer to detect leakage around the packer during pumping and enabled collection of head data during pumping for estimating hydrologic properties. Analytical results were obtained in about an hour from an on-site mobile laboratory equipped with a gas chromalograplvmass spectrometer (GC/MS). With the hydrologic and analytical results in hand, a decision was made to either complete the well or continue drilling to the next test zone. In almost every case, analytical results of ground water samples taken from the newly installed wells closely replicated the water quality of ground water samples obtained through the pump packer system.     

Fluctuations of electrical conductivity as a natural tracer for bank filtration in a losing stream

A key parameter used in the assessment of bank filtration is the travel time of the infiltrated river water during the passage through groundwater. We analyze time series of electrical conductivity (EC) in the river and adjacent groundwater observation wells to investigate travel times of young hyporheic groundwater in adjoining channelized and restored sections of River Thur in North-East Switzerland. To quantify mixing ratios and mean residence times we perform cross-correlation analysis and non-parametric deconvolution of the EC time series. Measurements of radon-222 in the groundwater samples validate the calculated residence times. A simple relationship between travel time and distance to the river has not been observed. Therefore, we speculate that the lateral position and depth of the thalweg as well as the type of bank stabilization might control the infiltration processes in losing rivers. Diurnal oscillations of EC observed in the river and in nearby observation wells facilitate analyzing the temporal variation of infiltration. The diurnal oscillations are particularly pronounced in low flow situations, while the overall EC signal is dominated by individual high-flow events. Differences in travel times derived from diurnal and overall EC signals thus reflect different infiltration regimes.     

Empirical relationships for estimating stream depletion by a well pumping near a gaining stream

Siting wells near streams requires an accurate estimate of the quantity of water derived from the river due to pumping. A number of hydrogeological and hydraulic parameters influence this value. This study estimates stream depletion under steady-state conditions for a variety of hydrogeological systems. A finite differences model was used to analyze several hydrogeological situations, and for each of these the stream depletion was estimated using an advective transport method. An empirical equation for stream depletion was obtained for the case of a stream that partially penetrates the aquifer and a pumping well that is screened over a portion of the aquifer. The derived equation, which is valid for both isotropic and anisotropic conditions, expresses stream depletion as a function of the unit inflow to the river, the discharge of the pumping well, the well screen length, the distance between the river and pumping well, the wetted perimeter, and a new parameter called "overlap," which is defined to be the distance between the riverbed and the top of well screen. The overlap parameter makes it possible to consider indirectly the vertical component of flow, which is accentuated when the well is screened below the streambed. The formula proposed here should be useful in deciding where to locate a pumping well and to decide the appropriate length of its screen.     

Pestizide in Oberflächenwässern — Befunde in Talsperren und in Uferfiltraten Pesticides in Surface Water — Analytical Results for Drinking Water Reservoirs and Bank Filtrate Waters

The results of a series of investigations are presented, performed by the governmental chemical laboratory Karlsruhe (CLUA) in order to determine pesticides in drinking water for purposes of food monitoring. 231 samples of 8 municipal water supplies (groundwater wells/no bank filtrate) and 34 private water supplies close to the bank of the river Rhine (groundwater wells/partly bank filtrate) were analysed. The sampling sites were located between Mannheim and Greffern, Germany. In 1986 high amounts of atrazine and simazine were determined in the drinking water wells of several private water supplies, especially at the sampling sites in the south. The likewise high chloride content of the samples was taken as a proof of bank filtrate contributing to the water. The decrease of atrazine load in the river water of the Rhine since 1987 has its parallel in the analytical results for the concerned drinking water wells. Obviously the pesticide contamination of water from bank filtrate pumping wells is mainly depending on the pesticide load of the river water. The application of fertilizers and pesticides in the fields seems not to be a significant source of the pesticide contamination of dringing waters derived from bank filtrate water. For reasons of public health a monitoring of pesticide content in bank filtrate pumping wells of private water supplies is recommended.     

Optimal pumping locations of skimming wells

Abstract

A real-life problem involving pumping of groundwater from a series of existing wells along a river flood plain underlain with geologically saline water is examined within a conceptual framework. Unplanned pumping results in upconing of saline water. Therefore, it is necessary to determine optimal locations of fixed capacity pumping wells in space and time from a set of pre-selected candidate wells that minimize total salinity concentration in space and time. The nonlinear, non-convex, combinatorial problem involving zero—one decision variables is solved in a simulation—optimization (S/O) framework. Optimization is accomplished by using simulated annealing (SA)—a search algorithm. The computational burden is primarily managed by replacing the numerical model with a surrogate simulator—artificial neural network (ANN). The computational burden is further reduced through intuitive algorithmic guidance. The model results suggest that the skimming wells must be operated from optimal locations such that they are staggered in space and time to obtain least saline water.     

An Idealized Ground-Water Flow and Chemical Transport Model (S-PATHS)

The number of studies on the actual and potential environmental consequences of contaminated ground water is growing. One means of studying these consequences is through an idealized flow and transport model, S-PATHS, which allows the hydrologist to determine the salient features of contaminant migration with a minimum of data. The transport of contaminants by ground water from many waste disposal sites can be geometrically idealized as flow between a line and a circle. The flow system adjacent to the disposal site can be represented as a contaminant line source, and a downgradient pumping well as a circular sink. To study waste disposal sites on a larger scale the model geometry is reversed and the disposal site is represented as a circular source, and a river or other convenient line of evaluation is represented as a line sink. This idealization allows S-PATHS to describe the flow and transport process directly by a single partial differential expression. S-PATHS considers transmissivity, effective porosity, sorption, source strength, source concentration, decay, potentiometric gradient, circle size, and distance to the line. Coding for the model is not lengthy and can be run on a large-capacity, hand-held calculator.     

Monitoring strategies at phreatic wellfields: a 3D travel time approach

Ground water quality networks for monitoring phreatic drinking water wellfields are generally established for two main purposes: (1) the short-term safeguarding of public water supply and (2) signaling and predicting future quality changes in the extracted ground water. Six monitoring configurations with different well locations and different screen depths and lengths were evaluated using a numerical model of the 3D ground water flow toward a partially penetrating pumping well in a phreatic aquifer. Travel times and breakthrough curves for observation and pumping wells were used to judge the effectiveness of different design configurations for three monitoring objectives: (1) early warning; (2) prediction of future quality changes; and (3) evaluation of protection measures inside a protection zone. Effectiveness was tested for scenarios with advective transport, first-order degradation, and linear sorption. It is shown that the location and especially the depth of the observation wells should be carefully chosen, taking into account the residence time from the surface to the observation well, the residual transit times to the extraction well, and the transformation and retardation rates. Shallow monitoring was most functional for a variety of objectives and conditions. The larger the degradation rates or retardation, the shallower should the monitoring be for effective early warning and prediction of future ground water quality. The general approach followed in the current study is applicable for many geohydrological situations, tuning specific monitoring objectives with residence times and residual transit times obtained from a site-specific ground water flow model.     

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.     

A modeling framework for the design of collector wells

We present results of a design study performed for the Saylorville Wellfield in Iowa, which is owned and operated by the Des Moines Water Works. The purpose of this study was to estimate wellfield capacity and provide a preliminary design for two radial collector wells to be constructed in the outwash aquifer along the Des Moines River near Saylorville, Iowa. After a field investigation to characterize the aquifer, regional two-dimensional and local three-dimensional, steady-state groundwater flow modeling was performed to locate and design the wells. This modeling was the foundation for design recommendations based on the relative performance of 12 collector well designs with varying lateral numbers, elevations, screen lengths, and orientations. For each site, alternate designs were evaluated based on model estimates of the capacity, the percent of surface water captured, and the production per unit length of screen. Many of our results are consistent with current design practices based on experience and intuition, but our methods allow for a quantitative approach for comparing alternate designs. Although the results are site-specific, the framework for evaluating the hydraulic design of the Saylorville radial collector wells is broadly applicable and could be used at other riverbank filtration sites. In addition, many of the conclusions from this design study may apply at other sites where construction of radial collector wells is being considered.     

Does Pumping Volume Affect the Concentration of Micropollutants in Groundwater Samples?

Information on concentrations of micropollutants (such as pharmaceuticals, pesticides, and industrial chemicals) in most highly dynamic riverbank filtration (RBF) systems is lacking, in contrast to data on standard chemical parameters. Sampling protocols have thus far been based on the stabilization of standard chemical parameters in relatively pristine environments. To determine whether groundwater samples for micropollutant analysis can be taken at a similar pumping volume as samples for testing standard chemical parameters in both environments, three groundwater monitoring wells in an RBF system were sampled at two points in time (after pumping of 3 well volumes and after pumping of 15 well volumes). Micropollutant concentrations were not significantly different between the two sampling points; therefore, appropriate samples can be drawn after pumping 3 well volumes. For a specific microbiological parameter (leucin incorporation), a statistically significant difference was found.