Simulation and parameterization of superimposed ice formation

In cold Arctic snowpacks, meltwater retention is a significant factor controlling the timing and magnitude of runoff. Meltwater percolates vertically through the snowpack until it reaches an impermeable horizon, whereupon a saturated zone is established. If the underlying media is below the freezing point, accretive ice formation takes place. This process has previously been crudely parameterized or modelled numerically. Such ice is called either superimposed ice on glaciers or basal ice on bare land. Using theory derived from sea‐ice formation, an analytical solution to basal ice growth is proposed. Results are compared against growth rates derived from numerical modelling. In addition, model results are compared to field observations of ice temperatures. The analytical solution is further extended to account for the temperature gradient inside the underlying media and the variable thermal properties of the underlying media. In the analysis, observations and references have predominantly relied on knowledge from glaciers. However, the process of accretive ice growth is equally important in seasonal snow packs with a cold snow‐ground interface and on Arctic sea ice where the ice‐snow interface is well below freezing point. The simplification of this accretive ice growth problem makes the solution attractive for incorporation in large‐scale cryospheric models. Copyright © 2007 John Wiley & Sons, Ltd.     

Modeling preferential flow in reactive barriers: implications for performance and design

Reactive barriers are passive and in situ ground water treatment systems. Heterogeneities in hydraulic conductivity (K) within the aquifer-reactive barrier system will result in higher flux rates, and reduced residence times, through portions of the barrier. These spatial variations in residence time will affect the treatment capacity of the barrier. A numerical flow model was used to evaluate the effects of spatial variations in K on preferential flow through barriers. The simulations indicate that the impact of heterogeneities in K will be a function of their location and distribution; the more localized the high K zone, the greater the preferential flow. The geometry of the reactive barrier will also strongly influence flow distribution. Aquifer heterogeneities will produce greater preferential flow in thinner barriers compared to thicker barriers. If the barrier K is heterogeneous, greater preferential flow will occur in thicker barriers. The K of the barrier will affect the flow distribution; decreasing the K of the barrier can result in more even distribution of flow. Results indicate that less variable flow will be attained utilizing thicker, homogeneous barriers. The addition of homogeneous zones to thinner barriers will be effective at redistributing flow only if installed immediately adjacent to both the up- and downgradient faces of the barrier.     

Laboratory Column Experiments Using Polymer Mats to Remove Selected VOCs, PAHs, and Pesticides from Ground Water

Large-scale column experiments were undertaken to evaluate the potential of polymer mats to remove selected volatile organic compounds, polycyclic aromatic hydrocarbons, and pesticides (atrazine and fenamiphos) from ground water and potentially to act as permeable reactive barriers in contaminated ground water environments. The polymer mats, composed of interwoven silicone (dimethylsiloxane) tubes and purged with air, were installed in 2 m long flow-through columns. The polymer mats proved efficient in physically removing (stripping) benzene and naphthalene from contaminated water. Removal efficiencies for both these compounds from an aqueous phase flowing past a polymer mat were 75% or greater. However, for atrazine and fenamiphos, removal efficiencies were 5% or less, probably as a result of their lower Henry's law constants and possibly lower polymer diffusion coefficients.
These experiments indicate that, at least for relatively volatile compounds, polymer mats can provide a remediation technique for the removal of organic compounds from contaminated water. Application of this technique may be well suited as a longer-term, semipassive strategy to remediate contaminated ground water, using natural ground water flow to deliver contaminated ground water to polymer mats engineered as sorption-stripping barriers.
Additional benefits of this technique may include targeted delivery of gaseous chemical amendments, such as oxygen, to enhance aerobic biodegradation and to further reduce any residual concentrations of contaminants.     

Pumping-induced drawdown and stream depletion in a leaky aquifer system

The impact of ground water pumping on nearby streams is often estimated using analytic models of the interconnected stream-aquifer system. A common assumption of these models is that the pumped aquifer is underlain by an impermeable formation. A new semianalytic solution for drawdown and stream depletion has been developed that does not require this assumption. This solution shows that pumping-induced flow (leakage) through an underlying aquitard can be an important recharge mechanism in many stream-aquifer systems. The relative importance of this source of recharge increases with the distance between the pumping well and the stream. The distance at which leakage becomes the primary component of the pumping-induced recharge depends on the specific properties of the aquifer, aquitard, and streambed. Even when the aquitard is orders of magnitude less transmissive than the aquifer, leakage can be an important recharge mechanism because of the large surface area over which it occurs. Failure to consider aquitard leakage can lead to large overestimations of both the drawdown produced by pumping and the contribution of stream depletion to the pumping-induced recharge. The ramifications for water resources management and water rights adjudication can be significant. A hypothetical example helps illustrate these points and demonstrates that more attention should be given to estimating the properties of aquitards underlying stream-aquifer systems. The solution presented here should serve as a relatively simple but versatile tool for practical assessments of pumping-induced stream-aquifer interactions. However, this solution should not be used for such assessments without site-specific data that indicate pumping has induced leakage through the aquitard.     

Sub-Saharan African Ground Water Protection—Building on International Experience

Sub-Saharan Africa faces significant challenges in dealing with ground water pollution. These countries can look to successes and missteps on other continents to help choose their own individual paths to ensuring reliable and clean supplies of ground water. In the large view, sub-Saharan Africa can define specific levels of acceptable risk in water quality that drive cleanup efforts and are amenable to acceptance across national and geographic boundaries. Ground water quality databases must be expanded, and data must be available in an electronic form that is flexible, expandable, and uniform, and that can be used over wide geographic areas. Guidance from other continents is available on well construction, sampling and monitoring, interim remediation, technical impracticability, monitored natural attenuation, and many specific issues such as how to deal with small waste generators and septic contamination of water supply wells. It is important to establish a common African view on the appropriateness of other nations’ ground water quality guidance for African issues, economic conditions, and community circumstances. Establishing numerical, concentration-based, water quality action levels for pollutants in ground water, which many neighboring African nations could hold comparable, would set the stage for risk-based remediation of contaminated sites. Efforts to gain public, grass-roots understanding and support for stable and balanced enforcement of standards are also key. Finally, effective capacity building in the region could be an eventual solution to ground water quality problems; with increased numbers of trained environmental professionals, ground water throughout the region can be protected and contaminated sites cleaned up.     

Modeling of Plume Capture by Continuous, Low-Permeability Barriers

Modeling was performed to simulate ground water flow through reactive barriers of lower hydraulic conductivity than the surrounding aquifer to determine the plume capture widths. As a plume approaches such a barrier, it spreads laterally. Therefore, to intercept an entire plume, the barrier must be wider than the upgradient width of the undisturbed plume. The results indicate that, for practical values of barrier thickness and plume width, hydraulic conductivities ten-fold less than that of the aquifer can be accommodated by making the width of the barrier approximately 20% greater than the upgradient width of the plume. Barrier hydraulic conductivities one-hundred-fold less than that of the aquifer may require barrier widths up to twice the width of the upgradient plume for plumes 100 feet wide (33 m) and as little as 1.1 times for plumes 1000 feet wide (325 m). The results presented here lend support to the view that novel emplacement methods that create zones of slightly lower hydraulic conductivity than the native aquifer may be viable alternatives to the excavation-and-backfill approaches which have thus far been used for installing permeable reactive barriers.     

Stream depletion in alluvial valleys using the SDF semianalytical model

A semianalytical method commonly used for quantifying stream depletion caused by ground water pumping was reviewed for applicability in narrow alluvial aquifers. This stream depletion factor (SDF) method is based on the analytic Glover model, but uses a numerical model-derived input parameter, called the SDF, to partly account for mathematically nonideal conditions such as variable transmissivity and nearby aquifer boundaries. Using the SDF can improve and simplify depletion estimates. However, the method's approximations introduce error that increases with proximity to the impermeable aquifer boundary. This article reviews the history of the method and its assumptions. New stream depletion response curves are presented as functions of well position within bounded aquifers. A simple modification to modeled SDF values is proposed that allows the impermeable boundary to be accounted for with image wells, but without overaccounting for boundary effects that are already reflected in modeled SDFs. It is shown that SDFs for locations closer to the river than to the aquifer boundary do not reflect impermeable-boundary effects, and thus need no modification, and boundary effects in the other portion of the aquifer follow a predictable removable pattern. This method is verified by comparing response curves using modified SDFs with response curves from an extensively calibrated numerical model of a managed ground water recharge site. The modification improves SDF-based stream depletion estimates in bounded aquifers while still benefiting from the additional information contained in SDF maps and retaining their value as standardized references for water rights administration.     

Injection-extraction treatment well pairs: an alternative to permeable reactive barriers

Two of the biggest drawbacks of using permeable reactive barriers (PRBs) to treat contaminated ground water are the high capital cost of installation, particularly when the contaminated ground water is deep below ground surface, and the uncertainty of whether or not PRBs remain effective for the long time scales (e.g., decades) needed for many contaminant plumes. The use of an injection-extraction treatment well pair (IETWP) for capture and treatment of contaminated ground water can circumvent these difficulties, while still providing many of the same advantages offered by PRBs. In this paper, the hydraulics of IETWPs and PRBs are compared, focusing primarily on the width of the captured plume. It is demonstrated that IETWPs act as hydraulic barriers in a manner similar to PRBs, and that IETWPs provide excellent plume capture. A mathematical expression is presented for the plume capture width of an IETWP oriented perpendicular to the ground water flow direction in a homogeneous aquifer. Also discussed are other practical considerations that might determine whether an IETWP is better suited than a PRB for a particular contaminated site; these considerations include operating and maintenance costs, and the conditions under which an IETWP system can be used for in situ remediation.     

The Effect of Undetected Barriers on Groundwater Drawdown and Recovery

In large-scale pumping projects, such as mine dewatering, predictions are often made about the rate of groundwater level recovery after pumping has ceased. However, these predictions may be impacted by geological uncertainty—including the presence of undetected impermeable barriers. During pumping, an impermeable barrier may be undetected if it is located beyond the maximum extent of the cone of depression; yet it may still control drawdown during the recovery phase. This has implications for regional-scale modeling and monitoring of groundwater level recovery. In this article, non-dimensional solutions are developed to show the conditions under which a barrier may be undetected during pumping but still significantly impact groundwater level recovery. The magnitude of the impact from an undetected barrier will increase as the ratio of pumping rate to aquifer transmissivity increases. The results are exemplified for a hypothetical aquifer with an unknown barrier 3 km from a pumping well. The difference in drawdown between a model with and without a barrier may be <1 m in the 10 years while pumping is occurring, but up to 50 m after pumping has ceased.     

B.—UNESCO

Summary

Methods are outlined for determining the hydraulic properties of an aquifer, namely, the transmissivity and storativity from a well located near one or more impermeable barriers. The procedures are based on the theory of the nonsteady radial flow toward a well steadily discharging from an aquifer of semiinfinite areal extent. Two different flow systems have been considered, namely a well near an impermeable barrier and a well located between two impermeable barriers.     

Multiobjective analysis of a public wellfield using artificial neural networks

As competition for increasingly scarce ground water resources grows, many decision makers may come to rely upon rigorous multiobjective techniques to help identify appropriate and defensible policies, particularly when disparate stakeholder groups are involved. In this study, decision analysis was conducted on a public water supply wellfield to balance water supply needs with well vulnerability to contamination from a nearby ground water contaminant plume. With few alternative water sources, decision makers must balance the conflicting objectives of maximizing water supply volume from noncontaminated wells while minimizing their vulnerability to contamination from the plume. Artificial neural networks (ANNs) were developed with simulation data from a numerical ground water flow model developed for the study area. The ANN-derived state transition equations were embedded into a multiobjective optimization model, from which the Pareto frontier or trade-off curve between water supply and wellfield vulnerability was identified. Relative preference values and power factors were assigned to the three stakeholders, namely the company whose waste contaminated the aquifer, the community supplied by the wells, and the water utility company that owns and operates the wells. A compromise pumping policy that effectively balances the two conflicting objectives in accordance with the preferences of the three stakeholder groups was then identified using various distance-based methods.     

The effects of vertical barrier walls on the hydraulic control of contaminated groundwater

We present explicit analytic solutions describing the hydraulic head and discharge vector for two-dimensional, steady groundwater flow past an impermeable barrier embedded in a regional flow field. We use the solution to investigate the effects of open vertical barriers on the flow field; in particular, we examine the hydraulic containment of contaminant plumes or source zones by combination of a vertical barrier wall and extraction wells. We quantify the local reduction in discharge rates due to the barrier wall and the local increase in the size of the capture zone of an extraction well near an open, up-gradient barrier. We find that the combination of an open vertical barrier with down-gradient extraction wells can be very effective in decreasing the well discharge rate necessary to control a contaminant plume or source area. Design charts are presented for quantifying the effects of the barrier wall on the hydraulic control of the groundwater flow field and for estimating the jump in head across a barrier. The charts are appropriate for use in the preliminary design and cost estimating of remedial systems, and for the design of dewatering systems.     

A Study on Confined Flow of Ground Water Through a Tunnel

Confined flow of ground water through a tunnel, which might be encountered in tunneling under the bottom of a sea or river, is numerically analyzed by a reductive finite element method formulated in our research. That is, the rate and potential distribution of the confined flow of ground water through an opening are obtained in connection with the permeability of rock masses, the thickness of covered ground, the location of impermeable bedrock, and other variables. In addition, flow through an opening in the ground with highly permeable masses and discharge of ground water through a tunnel in grouted masses are illustrated, and some useful results for the practical application of tunneling work are obtained.     

Screening of one-well hydraulic barrier design alternatives

This study develops a robust method for screening one-well hydraulic barrier design alternatives that can be easily computed without a numerical simulation model. The paper outlines the general method and shows its implementation with hydraulic barriers using a single pumping well. For such barriers, the method is easily computable with spreadsheets and/or charts depicted within the paper and posted online. The method applies the potential flow theory, which leads to using a curvilinear coordinate system for all types of calculations. For contaminant transport calculations, the method applies the boundary layer theory. For calculations of aquifer remediation, the method refers to bulk characteristics of the domain. As an example, the method has been applied to calculate the possible containment of a wide part of the coastal plain aquifer in Israel, which is contaminated by entrapped kerosene (a light nonaqueous phase liquid).     

Design Screening Tools for Passive Funnel and Gate Systems

The funnel and gate remediation concept (Star and Cherry 1993) represents a promising, yet relatively under-developed, technology for the passive control and in situ remediation of contaminated ground water. Effective design and implementation of such a system may, however, prove difficult under conditions of large or unpredictable variations in contaminant migration or ground water flow.
Numerical modeling of two-dimensional ground water flow has been used to predict the hydraulic performance of passive, straight, or winged funnel and gate configurations over a range of hydrogeologic and ambient ground water flow conditions. The results of these analyses were used to construct generic correlation diagrams relating upstream capture zone or gale through put to the barrier, gale, and aquifer characteristics. These diagrams serve as useful screening tools to (1) quantitatively estimate the capture zone of pre-determined funnel and gale configurations, or (2) develop preliminary funnel and gale designs that will yield a desired capture zone, independent of aquifer characteristics.     

MNA as a Remedy for Arsenic Mobilized by Anthropogenic Inputs of Organic Carbon

The potential application of monitored natural attenuation (MNA) as a remedy for ground water contaminated with arsenic (As) is examined for a subset of contaminated sites, specifically those where naturally occurring As has been mobilized due to localized anthropogenic organic carbon (OC) releases. This includes sites subject to petroleum releases, exposure to landfill leachates, and OC additions for biostimulation of reductive dechlorination of chlorinated solvents. The key characteristic of these sites is that, under conditions prevailing before the anthropogenic OC introduction, the naturally occurring As in the subsurface was not mobile and did not adversely affect ground water quality. This suggests that, in the far-field (where background conditions are (re) established), As may be sequestered upon contact of the contaminated ground water with either or both the (uncontaminated) ambient ground water and the background aquifer minerals. The observed extents of elevated concentrations (or "footprints") of As and other chemical species, such as dissolved OC and iron (Fe), and related parameters, such as redox potential ( E _h) and dissolved oxygen, and their evolution over time can be used to assess the mobilization and sequestration of As and the potential feasibility of MNA as a remedial option. Ultimately, the capacity for As sequestration must be assessed in the context of the OC loading to the site, which may require "active" measures for source control. Monitoring is needed to confirm the continuing effectiveness of the MNA remedy or to indicate if contingency measures must be implemented.     

Influence of soil moisture hysteresis on the functioning of capillary barriers

Previous experimental studies of capillary barriers have identified highly hysteretic soil moisture retention characteristics in the materials used. In this study, numerical modelling is used to analyse the role of soil moisture hysteresis in capillary barrier functioning. Comparisons between published experimental results and model simulations indicate that soil moisture hysteresis was a necessary inclusion in the modelling approach to adequately reproduce pore water pressure distributions and the timing of breakthrough occurrences. Under hypothetical intermittent infiltration and evaporation conditions, the predicted volumetric water content in the moisture retention layer was significantly different for hysteretic and non‐hysteretic models. The hysteresis effect was found to be dependent on the nature of infiltration–evaporation cycling, although the predicted volume of flow through the hysteretic barrier was lower than that of the non‐hysteretic case, regardless of the nature of the cyclic upper boundary conditions. For practical engineering designs, where the water leakage through the barrier is the primary concern, the inclusion of soil moisture hysteresis in numerical modelling is needed. Copyright © 2009 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.     

The effects of rock fragments on hydrologic and hydraulic responses along a slope

A model describing the three‐dimensional matrix flow along a slope with rock fragments or impermeable blocks was developed. The model was combined with modified Picard's iteration to ensure mass conservation in the unsaturated flow. We found that rock fragments obstruct water flow along the slope. The groundwater table must be raised to provide a sufficient pore water pressure gradient to facilitate water flow, but higher pore water pressure may induce slope failure. We also conducted a bench‐scale laboratory flume experiment to examine the effects of impermeable blocks on downstream seepage flow. In addition, a numerical experiment was conducted to examine how different arrangements of impermeable blocks affect downstream seepage flow and pore water pressure. This research demonstrated that the hydraulic phenomena were affected when impermeable blocks were present, and pore water pressure increased as the position of impermeable blocks was lowered. Copyright © 2006 John Wiley & Sons, Ltd.     

Monitoring the Removal of Phosphate from Ground Water Discharging through a Pond-Bottom Permeable Reactive Barrier

Installation of a permeable reactive barrier to intercept a phosphate (PO₄) plume where it discharges to a pond provided an opportunity to develop and test methods for monitoring the barrier's performance in the shallow pond-bottom sediments. The barrier is composed of zero-valent-iron mixed with the native sediments to a 0.6-m depth over a 1100-m² area. Permanent suction, diffusion, and seepage samplers were installed to monitor PO₄ and other chemical species along vertical transects through the barrier and horizontal transects below and near the top of the barrier. Analysis of pore water sampled at about 3-cm vertical intervals by using multilevel diffusion and suction samplers indicated steep decreases in PO₄ concentrations in ground water flowing upward through the barrier. Samples from vertically aligned pairs of horizontal multiport suction samplers also indicated substantial decreases in PO₄ concentrations and lateral shifts in the plume's discharge area as a result of varying pond stage. Measurements from Lee-style seepage meters indicated substantially decreased PO₄ concentrations in discharging ground water in the treated area; temporal trends in water flux were related to pond stage. The advantages and limitations of each sampling device are described. Preliminary analysis of the first 2 years of data indicates that the barrier reduced PO₄ flux by as much as 95%.