Optimal groundwater contamination monitoring using pumping wells.

This study presents a new method for selecting monitoring wells for optimal evaluation of groundwater quality. The basic approach of this work is motivated by difficulties in interpolating groundwater quality from information collected for only few sampled wells. The well selection relies on other existing data relevant to contaminant distribution in the sampling domain, e.g. predictions of models which rely on past measurements. The objective of this study is to develop a method of selecting the optimal wells, from which measurements could best serve some external model, e.g. a kriging system for characterizing the entire plume distribution, a flow-and-transport model for predicting a future distribution, or an inverse model for locating contaminant sources or estimating aquifer parameters. The decision variable at each sampling round determines the specific wells to be sampled. The study objective is accomplished through a spatially-continuous utility density function (UDF) which describes the utility of sampling at every point. The entire methodology which utilizes the UDF in conjunction with a sampling algorithm is entitled the UDF method. By applying calculations in steady and unsteady state sampling domains the effectiveness of the UDF method is demonstrated.

Implementing a method of screening one-well hydraulic barrier design alternatives.

This article provides details of applying the method developed by the authors (Rubin et al. 2008b) for screening one-well hydraulic barrier design alternatives. The present article with its supporting information (manual and electronic spreadsheets with a case history example) provides the reader complete details and examples of solving the set of nonlinear equations developed by Rubin et al. (2008b). It allows proper use of the analytical solutions and also depicting the various charts given by Rubin et al. (2008b). The final outputs of the calculations are the required position and the discharge of the pumping well. If the contaminant source is nonaqueous phase liquid (NAPL) entrapped within the aquifer, then the method provides an estimate of the aquifer remediation progress (which is a by-product) due to operating the hydraulic barrier.

Screening of one-well hydraulic barrier design alternatives.

Abstract 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).

Parameters that control the cleanup of fractured permeable aquifers.

This study develops a modeling approach for simulating and evaluating entrapped light nonaqueous-phase liquid (light NAPL-LNAPL) dissolution and transport of the solute in a fractured permeable aquifer (FPA). The term FPA refers to an aquifer made of porous blocks of high permeability that embed fractures. The fracture network is part of the domain characterized by high permeability and negligible storage. Previous studies show that sandstone aquifers often represent FPAs. The basic model developed in this study is a two-dimensional (2-D) model of permeable blocks that embed oblique equidistant fractures with constant aperture and orientation. According to this model, two major parameters govern NAPL dissolution and transport of the solute. These parameters are: 1) the dimensionless interphase mass transfer coefficient, K(f0), and 2) the mobility number, N(M0). These parameters represent measures of heterogeneity affecting flow, NAPL dissolution, and transport of the solute in the domain. The parameter K(f0) refers to the rate at which organic mass is transferred from the NAPL into the water phase. The parameter N(M0) represents the ratio of flow through the porous blocks to flow through the fracture network in regions free of entrapped NAPL. It also provides a measure of groundwater flow bypassing regions contaminated by entrapped NAPL. In regions contaminated by entrapped NAPL our simulations have often indicated very low permeability of the porous blocks, enabling a significant increase of the fracture flow at the expense of the permeable block flow. Two types of constitutive relationships also affect the rate of FPA cleanup: 1) the relationship between the saturation of the entrapped NAPL and the permeability of the porous blocks, and 2) the relationships representing effects of the entrapped NAPL saturation and the permeable block flow velocity on rates of interphase mass transfer. This study provides basic tools for evaluating the characteristics of pump-and-treat cleanup of FPAs by referring to sets of parameters and constitutive relationships typical of FPAs. The numerical simulations carried out in this study show that at high initial saturation of the entrapped NAPL, during initial stages of the FPA cleanup the contaminant concentration increases, but later it decreases. This phenomenon originates from significant groundwater bypassing the NAPL entrapped in the permeable blocks via the fracture network.

Groundwater contamination downstream of a contaminant penetration site. I. Extension-expansion of the contaminant plume.

This study concerns the possible use of boundary layer (BL) approach for the analysis and evaluation of contaminant transport in groundwater due to contaminant penetration into the groundwater aquifer through a site of limited size. The contaminant penetration may occur through either the upper (surface) or lower (bedrock) boundary of the aquifer. Two general cases of contaminant penetration mechanisms are considered: (1) the contaminant is transferred through an interface between a contaminating and freshwater fluid phases, and (2) the contaminant arrives at groundwater by leakage and percolation. For the purpose of BL evaluation the contaminant plume is divided into three different sections: (1) the penetration section, (2) the extension-expansion section, and (3) the spearhead section. In each section a different BL method approach yields simple analytical expressions for the description of the contaminant plume migration and contaminant transport. Previous studies of the BL method can be directly applied to the evaluation of contaminant transport at the contaminant penetration section. The present study extends those studies and concerns the contaminant transport in the two other sections, which are located downstream of the penetration section. This study shows that the contaminant concentration profiles in sections 2 and 3 incorporate two BLs: (1) an inner BL adjacent to the aquifer bottom or surface boundary, and (2) an outer BL, which develops above or below the inner one. The method developed in the present study has been applied to practical issues concerning salinity penetration into groundwater in south central Kansas.

Groundwater contamination downstream of a contaminant penetration site. II. Horizontal penetration of the contaminant plume.

Part I of this study (Rubin, H.; Buddemeier, R.W. Groundwater Contamination Downstream of a Contaminant Penetration Site Part 1: Extension-Expansion of the Contaminant Plume. J. of Environmental Science and Health Part A (in press).) addressed cases, in which a comparatively thin contaminated region represented by boundary layers (BLs) developed within the freshwater aquifer close to contaminant penetration site. However, at some distance downstream from the penetration site, the top of the contaminant plume reaches the top or bottom of the aquifer. This is the location of the "attachment point," which comprises the entrance cross section of the domain evaluated by the present part of the study. It is shown that downstream from the entrance cross section, a set of two BLs develop in the aquifer, termed inner and outer BLs. It is assumed that the evaluated domain, in which the contaminant distribution gradually becomes uniform, can be divided into two sections, designated: (a) the restructuring section, and (b) the establishment section. In the restructuring section, the vertical concentration gradient leads to expansion of the inner BL at the expense of the outer BL, and there is almost no transfer of contaminant mass between the two layers. In the establishment section, each of the BLs occupies half of the aquifer thickness, and the vertical concentration gradient leads to transfer of contaminant mass from the inner to the outer BL. By use of BL approximations, changes of salinity distribution in the aquifer are calculated and evaluated. The establishment section ends at the uniformity point, downstream from which the contaminant concentration profile is practically uniform. The length of the restructuring section, as well as that of the establishment section, is approximately proportional to the aquifer thickness squared, and is inversely proportional to the transverse dispersivity. The study provides a convenient set of definitions and terminology that are helpful in visualizing the gradual development of uniform contaminant concentration distribution in an aquifer subject to contaminant plume penetration. The method developed in this study can be applied to a variety of problems associated with groundwater quality, such as initial evaluation of field data, design of field data collection, the identification of appropriate boundary conditions for numerical models, selection of appropriate numerical modeling approaches, interpretation and evaluation of field monitoring results, etc.