Influence of temporally variable groundwater flow conditions on point measurements and contaminant mass flux estimations

Monitoring of contaminant concentrations, e.g., for the estimation of mass discharge or contaminant degradation rates, often is based on point measurements at observation wells. In addition to the problem, that point measurements may not be spatially representative, a further complication may arise due to the temporal dynamics of groundwater flow, which may cause a concentration measurement to be not temporally representative. This paper presents results from a numerical modeling study focusing on temporal variations of the groundwater flow direction. “Measurements” are obtained from point information representing observation wells installed along control planes using different well frequencies and configurations. Results of the scenario simulations show that temporally variable flow conditions can lead to significant temporal fluctuations of the concentration and thus are a substantial source of uncertainty for point measurements. Temporal variation of point concentration measurements may be as high as the average concentration determined, especially near the plume fringe, even when assuming a homogeneous distribution of the hydraulic conductivity. If a heterogeneous hydraulic conductivity field is present, the concentration variability due to a fluctuating groundwater flow direction varies significantly within the control plane and between the different realizations. Determination of contaminant mass fluxes is also influenced by the temporal variability of the concentration measurement, especially for large spacings of the observation wells. Passive dosimeter sampling is found to be appropriate for evaluating the stationarity of contaminant plumes as well as for estimating average concentrations over time when the plume has fully developed. Representative sampling has to be performed over several periods of groundwater flow fluctuation. For the determination of mass fluxes at heterogeneous sites, however, local fluxes, which may vary considerably along a control plane, have to be accounted for. Here, dosimeter sampling in combination with time integrated local water flux measurements can improve mass flux estimates under dynamic flow conditions.     

Stochastic evaluation of mass discharge from pointlike concentration measurements

The contaminant mass discharge crossing a control plane is an important metric in the assessment of natural attenuation at contaminated sites. For risk-assessment purposes, the mass discharge must be estimated together with a level of uncertainty. We present a conditional Monte Carlo approach that allows estimating the statistical distribution of mass discharge. The approach is based on conditioning multiple realizations of the hydraulic conductivity field on all data available. We jointly determine a first-order decay coefficient in each realization, leading to conditional statistical distributions of all estimated parameters and the total mass discharge. The resulting statistical distribution of contaminant mass discharges can be used in the assessment of risks at the contaminated site. The method is applied to data of hypothetical test cases, which gives the opportunity to compare estimation results to the true field. As concentration data, we account for pointlike measurements obtained in multi-level sampling wells. The obtained empirical distribution of mass discharge crossing the multi-level sampling fence could be well fitted by a log-normal distribution.     

Contaminant degradation in physically and chemically heterogeneous aquifers

This paper examines the importance of the correlation between hydraulic conductivity (K) and degradation rate constant (k) during the transport of reactive contaminants in heterogeneous aquifers. We simulated reactive transport in an ensemble of two-dimensional heterogeneous aquifers. Two sets of transport simulations were conducted: one in which a perfect positive correlation was assumed between ln(K) and ln(k), and one in which a perfect negative correlation was assumed. We found that the sign of the correlation has important consequences for the contaminant transport. Qualitatively, a negative correlation leads to significantly more pronounced "fingering" of the contaminant plume than does a positive correlation, with potentially important consequences for downgradient receptors. Quantitatively, the expected behavior (as quantified by the contaminant mass remaining in the aquifer) is statistically different between the positive and negative cases: on average, more contaminant mass persists when K and k are negatively correlated. Also, the negative correlation leads to more variability between realizations of the ensemble, whereas a positive correlation induces relatively little variability between realizations. We discuss the implications of these findings for the management of contaminated aquifers.     

Plume persistence caused by back diffusion from thin clay layers in a sand aquifer following TCE source-zone hydraulic isolation

This paper concludes that back diffusion from one or a few thin clayey beds in a sand aquifer can cause contaminant persistence above MCLs in a sand aquifer long after the source zone initially causing the plume is isolated or removed. This conclusion is based on an intensive case study of a TCE contaminated site in Florida, with the processes evaluated using numerical modeling. At this site, the TCE DNAPL zone formed decades ago, and was hydraulically isolated by means of an innovative system performing groundwater extraction, treatment and re-injection. Treated water is re-injected in a row of injection wells situated a short distance downgradient of the extraction wells, creating a clean-water displacement front to efficiently flush the downgradient plume. This scheme avoids the creation of stagnation zones typical of most groundwater pump-and-treat systems, thereby minimizing the time for aquifer flushing and therefore downgradient cleanup. The system began operation in August 2002 and although the performance monitoring shows substantial declines in concentrations, detectable levels of TCE and degradation products persist downgradient of the re-injection wells, long after the TCE should have disappeared based on calculations assuming a nearly homogenous sand aquifer. Three hypotheses were assessed for this plume persistence: 1) incomplete source-zone capture, 2) DNAPL occurrence downgradient of the re-injection wells, and 3) back diffusion from one or more thin clay beds in the aquifer. After careful consideration, the first two hypotheses were eliminated, leaving back diffusion as the only plausible hypothesis, supported by detailed measurements of VOC concentrations within and near the clay beds and also by numerical model simulations that closely represent the field site hydrogeologic conditions. The model was also used to simulate a more generalized, hypothetical situation where more thin clayey beds occur in a sand aquifer with an underlying aquitard. While there is no doubt that DNAPL source mass reduction can eventually improve downgradient groundwater quality, the magnitude and time scale over which the improvement occurs is the major uncertainty given current characterization approaches. This study shows that even one thin clay bed, less than 0.2 m thick, can cause plume persistence due to back diffusion for several years or even decades after the flux from the source is completely isolated. Thin clay beds, which have a large storage capacity for dissolved and sorbed contaminant mass, are common in many types of sandy aquifers. However, without careful inspection of continuous cores and sampling, such thin clay beds, and their potential for causing long-term back-diffusion effects, can easily go unnoticed during site characterization.     

Pulsed pumping process optimization using a potential flow model

A computational model is applied to the optimization of pulsed pumping systems for efficient in situ remediation of groundwater contaminants. In the pulsed pumping mode of operation, periodic rather than continuous pumping is used. During the pump-off or trapping phase, natural gradient flow transports contaminated groundwater into a treatment zone surrounding a line of injection and extraction wells that transect the contaminant plume. Prior to breakthrough of the contaminated water from the treatment zone, the wells are activated and the pump-on or treatment phase ensues, wherein extracted water is augmented to stimulate pollutant degradation and recirculated for a sufficient period of time to achieve mandated levels of contaminant removal. An important design consideration in pulsed pumping groundwater remediation systems is the pumping schedule adopted to best minimize operational costs for the well grid while still satisfying treatment requirements. Using an analytic two-dimensional potential flow model, optimal pumping frequencies and pumping event durations have been investigated for a set of model aquifer-well systems with different well spacings and well-line lengths, and varying aquifer physical properties. The results for homogeneous systems with greater than five wells and moderate to high pumping rates are reduced to a single, dimensionless correlation. Results for heterogeneous systems are presented graphically in terms of dimensionless parameters to serve as an efficient tool for initial design and selection of the pumping regimen best suited for pulsed pumping operation for a particular well configuration and extraction rate. In the absence of significant retardation or degradation during the pump-off phase, average pumping rates for pulsed operation were found to be greater than the continuous pumping rate required to prevent contaminant breakthrough.     

Analysis of groundwater contamination using concentration-time series recorded during an integral pumping test: bias introduced by strong concentration gradients within the plume

When only few monitoring wells are available to assess the extent and level of groundwater contamination, inversion of concentration breakthrough curves acquired during an integral pumping test can be used as an alternative quantification method. The idea is to use concentration-time series recorded during integral pumping tests through an inversion technique to estimate contaminant mass fluxes crossing a control plane. In this paper, we examine how a longitudinal concentration gradient along a contaminant plume length scale affects the estimated inversed-concentration distribution and its associated mass flux. The analytically inversed-concentration distribution at the imaginary control plane (ICP) is compared to a numerically generated concentration distribution, treating the latter one as a "real contaminant plume" characterized by the presence of a longitudinal concentration gradient. It is found that the analytically inversed-concentration can lead to overestimation or underestimation of concentration distribution values depending on the transport time period and dispersivity values. At lower dispersivity values, with shorter transport time periods, the analytically inversed-concentration distribution overestimates the "real" concentration distribution. A better fit of the estimated concentration distribution to the "real" one is observed when the transport time period increases, i.e. when the advective front has already crossed the ICP. However, for higher dispersivity values, underestimation of the real concentration distribution is observed. Deviation of the inversed-concentration distribution from the "real" one is assessed for a site-specific concentration gradient term. A concentration gradient adjusted contaminant mass flux is thus formulated to evaluate groundwater contamination levels at a given time period through an ICP. This concentration gradient ratio can indicate whether the ICP is well positioned to evaluate accurately contaminant mass fluxes which are representative of groundwater contamination levels.     

A Monitoring Network Design Procedure for Three-Dimensional (3D) Groundwater Contaminant Source Identification

Finding the location and concentration of contaminant sources is an important step in groundwater remediation and management. This discovery typically requires the solution of an inverse problem. This inverse problem can be formulated as an optimization problem where the objective function is the sum of the square of the errors between the observed and predicted values of contaminant concentration at the observation wells. Studies show that the source identification accuracy is dependent on the observation locations (i.e., network geometry) and frequency of sampling; thus, finding a set of optimal monitoring well locations is very important for characterizing the source. The objective of this study is to propose a sensitivity-based method for optimal placement of monitoring wells by incorporating two uncertainties: the source location and hydraulic conductivity. An optimality metric called D-optimality in combination with a distance metric, which tends to make monitoring locations as far apart from each other as possible, is developed for finding optimal monitoring well locations for source identification. To address uncertainty in hydraulic conductivity, an integration method of multiple well designs is proposed based on multiple hydraulic conductivity realizations. Genetic algorithm is used as a search technique for this discrete combinatorial optimization problem. This procedure was applied to a hypothetical problem based on the well-known Borden Site data in Canada. The results show that the criterion-based selection proposed in this paper provides improved source identification performance when compared to uniformly distributed placement of wells.     

Analysis of groundwater contamination using concentration–time series recorded during an integral pumping test: Bias introduced by strong concentration gradients within the plume

When only few monitoring wells are available to assess the extent and level of groundwater contamination, inversion of concentration breakthrough curves acquired during an integral pumping test can be used as an alternative quantification method. The idea is to use concentration–time series recorded during integral pumping tests through an inversion technique to estimate contaminant mass fluxes crossing a control plane. In this paper, we examine how a longitudinal concentration gradient along a contaminant plume length scale affects the estimated inversed-concentration distribution and its associated mass flux. The analytically inversed-concentration distribution at the imaginary control plane (ICP) is compared to a numerically generated concentration distribution, treating the latter one as a “real contaminant plume” characterized by the presence of a longitudinal concentration gradient. It is found that the analytically inversed-concentration can lead to overestimation or underestimation of concentration distribution values depending on the transport time period and dispersivity values. At lower dispersivity values, with shorter transport time periods, the analytically inversed-concentration distribution overestimates the “real” concentration distribution.A better fit of the estimated concentration distribution to the “real” one is observed when the transport time period increases, i.e. when the advective front has already crossed the ICP. However, for higher dispersivity values, underestimation of the real concentration distribution is observed. Deviation of the inversed-concentration distribution from the “real” one is assessed for a site-specific concentration gradient term. A concentration gradient adjusted contaminant mass flux is thus formulated to evaluate groundwater contamination levels at a given time period through an ICP. This concentration gradient ratio can indicate whether the ICP is well positioned to evaluate accurately contaminant mass fluxes which are representative of groundwater contamination levels.     

Monitoring groundwater contamination and delineating source zones at industrial sites: uncertainty analyses using integral pumping tests

Field-scale characterisations of contaminant plumes in groundwater, as well as source zone delineations, are associated with uncertainties that can be considerable. A major source of uncertainty in environmental datasets is due to variability of sampling results, as a direct consequence of the heterogeneity of environmental matrices. We develop a methodology for quantifying uncertainties in field-scale mass flow and average concentration estimations, using integral pumping tests (IPTs), where the contaminant concentration is measured as a function of time in a pumping well. This procedure increases the sampling volume and reduces the effect of small-scale variability that may bias point-scale measurements. In particular, using IPTs, the interpolation uncertainty of conventional point-scale measurements is transformed to a quantifiable uncertainty related to the (unknown) plume position relative to the pumping well. We show that this plume position uncertainty generally influenced the predicted mass flows and average concentrations (of acenapthene, benzene and CHCs) to a greater extent than a boundary condition uncertainty related to the local water balance, considering 19 control planes at a highly heterogeneous industrial site in southwest Germany. Furthermore, large (order of magnitude) uncertainties only occurred if the conditions were strongly heterogeneous in the nearest vicinity of the well. We also develop a consistent methodology for an assessment of the combined effect of uncertainty in hydraulic conditions and uncertainty in reactive transport parameters for delimiting of both contaminant source zones and zones absent of source, based on (downgradient) IPTs.     

Quantification of groundwater contamination in an urban area using integral pumping tests

In this paper, the integral groundwater investigation method is used for the quantification of PCE and TCE mass flow rates at an industrialized urban area in Linz, Austria. In this approach, pumping wells positioned along control planes perpendicular to the groundwater flow direction are operated for a time period on the order of days and sampled for contaminants. The concentration time series of the contaminants measured during operation of the pumping wells are then used to determine contaminant mass flow rates, mean concentrations and the plume shapes and positions at the control planes. The three control planes used in Linz were positioned downstream of a number of potential source zones, which are distributed over the field site. By use of the integral investigation method, it was possible to identify active contaminant sources, quantify the individual source strength in terms of mass flow rates at the control planes and estimate the contaminant plume position relative to the control planes. The source zones emitting the highest PCE and TCE mass flow rates could be determined, representing the areas where additional investigation and remediation activities will be needed. Additionally, large parts of the area investigated could be excluded from further investigation and remediation activities.     

Predicting natural attenuation of xylene in groundwater using a numerical model.

The aquifer beneath an abandoned refinery in the Lower Rhine area, Germany, was contaminated with a number of different mineral oil products. Groundwater sampling in the area around the former xylene plant revealed that a xylene plume had developed in the underlying groundwater, and moreover, that there is strong evidence for in situ microbial xylene degradation with oxygen, nitrate, sulfate and ferric iron as electron acceptors. In order to prevent further xylene spreading, three pumping wells extracting contaminated water were installed downgradient of the spill zone. The numerical reactive transport code Transport Biochemisty Chemistry (TBC) was applied to this situation to quantify the relation of microbial degradation to xylene removal by the pumping wells. It could be shown that the unamended in situ degradation was an appreciable xylene removal process that contributed to about one-third to the total xylene removal (degradation plus extraction). A further objective of the model application was to predict xylene spreading under regional flow conditions, i.e. without operation of the three pumping wells, to consider the possible effects of natural xylene attenuation. To accomplish this, the model calibrated for the situation with operating wells was transferred to the hydraulic situation of regional flow while retaining the parameters of the biochemical model. It turned out that the xylene plume that is expected to develop downgradient of the source area will be limited to an extension of not more than 1000 m. An interesting feature of the simulations results was that xylene degradation under iron-reducing conditions, which was of minor importance for the situation with operating pumping wells, becomes the dominant degradation mechanism under regional flow conditions. Moreover, iron reduction will be the key process in controlling plume evolution. The model application illustrates that multi-species reactive transport models are needed to adequately transfer reactive processes from one hydraulic situation to another, while single species models are not suited for this predictive task.     

Probabilistic study of well capture zones distribution at the Lauswiesen field site

The delineation of well capture zones is of utmost environmental and engineering relevance as pumping wells are commonly used both for drinking water supply needs, where protection zones have to be defined, and for investigation and remediation of contaminated aquifers. We analyze the probabilistic nature of well capture zones within the well field located at the "Lauswiesen" experimental site. The test site is part of an alluvial heterogeneous aquifer located in the Neckar river valley, close to the city of Tübingen in South-West Germany. We explore the effect of different conceptual models of the structure of aquifer heterogeneities on the delineation of three-dimensional probabilistic well catchment and time-related capture zones, in the presence of migration of conservative solutes. The aquifer is modeled as a three-dimensional, doubly stochastic composite medium, where distributions of geo-materials and hydraulic properties are uncertain. We study the relative importance of uncertain facies geometry and uncertain hydraulic conductivity and porosity on predictions of catchment and solute time of travel to the pumping well by focusing on cases in which (1) the facies distribution is random, but the hydraulic properties of each material are fixed, and (2) both facies geometry and material properties vary stochastically. The problem is tackled within a conditional numerical Monte Carlo framework. Results are provided in terms of probabilistic demarcations of the three-dimensional well catchment and time-related capture zones. Our findings suggest that the uncertainty associated with the prediction of the location of the outer boundary of well catchment at the "Lauswiesen" site is significantly affected by the conceptual model adopted to incorporate the heterogeneous nature of the aquifer domain in a predictive framework. Taking into account randomness of both lithofacies distribution and materials hydraulic conductivity allows recognizing the existence of preferential flow paths that influence the extent of the well catchment and the solute travel time distribution at the site.     

Estimating contaminant mass discharge: a field comparison of the multilevel point measurement and the integral pumping investigation approaches and their uncertainties

In this field study, two approaches to assess contaminant mass discharge were compared: the sampling of multilevel wells (MLS) and the integral groundwater investigation (or integral pumping test, IPT) that makes use of the concentration-time series obtained from pumping wells. The MLS approached used concentrations, hydraulic conductivity and gradient rather than direct chemical flux measurements, while the IPT made use of a simplified analytical inversion. The two approaches were applied at a control plane located approximately 40m downgradient of a gasoline source at Canadian Forces Base Borden, Ontario, Canada. The methods yielded similar estimates of the mass discharging across the control plane. The sources of uncertainties in the mass discharge in each approach were evaluated, including the uncertainties inherent in the underlying assumptions and procedures. The maximum uncertainty of the MLS method was about 67%, and about 28% for the IPT method in this specific field situation. For the MLS method, the largest relative uncertainty (62%) was attributed to the limited sampling density (0.63 points/m(2)), through a novel comparison with a denser sampling grid nearby. A five-fold increase of the sampling grid density would have been required to reduce the overall relative uncertainty for the MLS method to about the same level as that for the IPT method. Uncertainty in the complete coverage of the control plane provided the largest relative uncertainty (37%) in the IPT method. While MLS or IPT methods to assess contaminant mass discharge are attractive assessment tools, the large relative uncertainty in either method found for this reasonable well monitored and simple aquifer suggests that results in more complex plumes in more heterogeneous aquifers should be viewed with caution.     

An analytical solution for two-dimensional contaminant transport during groundwater extraction

We obtain an analytical solution for two-dimensional steady-state transport of conservative contaminant between injecting and pumping wells. Flow and transport are considered in the vertical cross-section. The Dupuit approximation and conformal mapping onto the complex potential domain are employed to determine the velocity and concentration distributions, respectively. We use this solution to derive a priori conditions under which widely used 1-D analytical solutions with constant velocity and dispersion coefficients provide accurate approximations. These conditions are formulated in terms of aquifer parameters, such as hydraulic conductivity, porosity and dispersivities, and remediation strategy, e.g., well spacing and pumping regimes.     

Flux-based assessment at a manufacturing site contaminated with trichloroethylene

Groundwater and contaminant fluxes were measured, using the passive flux meter (PFM) technique, in wells along a longitudinal transect passing approximately through the centerline of a trichloroethylene (TCE) plume at a former manufacturing plant located in the Midwestern US. Two distinct zones of hydraulic conductivity were identified from the measured groundwater fluxes; a 6-m-thick upper zone ( approximately 7 m to 13 m below the ground surface or bgs) with a geometric mean Darcy flux (q(0)) of 2 cm/day, and a lower zone ( approximately 13 m to 16.5m bgs) with a q(0) approximately 15 cm/day; this important hydrogeologic feature significantly impacts any remediation technology used at the site. The flux-averaged TCE concentrations estimated from the PFM results compared well with existing groundwater monitoring data. It was estimated that at least 800 kg of TCE was present in the source zone. The TCE mass discharge across the source control plane (85 m x 38 m) was used to estimate the "source strength" ( approximately 365 g/day), while mass discharges across multiple down-gradient control planes were used to estimate the plume-averaged, TCE degradation rate constant (0.52 year(-1)). This is close to the rate estimated using the conventional centerline approach (0.78 year(-1)). The mass discharge approach provides a more robust and representative estimate than the centerline approach since the latter uses only data from wells along the plume centerline while the former uses all wells in the plume.     

Assessing the impact of source-zone remediation efforts at the contaminant-plume scale through analysis of contaminant mass discharge

The long-term impact of source-zone remediation efforts was assessed for a large site contaminated by trichloroethene. The impact of the remediation efforts (soil vapor extraction and in-situ chemical oxidation) was assessed through analysis of plume-scale contaminant mass discharge, which was measured using a high-resolution data set obtained from 23 years of operation of a large pump-and-treat system. The initial contaminant mass discharge peaked at approximately 7kg/d, and then declined to approximately 2kg/d. This latter value was sustained for several years prior to the initiation of source-zone remediation efforts. The contaminant mass discharge in 2010, measured several years after completion of the two source-zone remediation actions, was approximately 0.2kg/d, which is ten times lower than the value prior to source-zone remediation. The time-continuous contaminant mass discharge data can be used to evaluate the impact of the source-zone remediation efforts on reducing the time required to operate the pump-and-treat system, and to estimate the cost savings associated with the decreased operational period. While significant reductions have been achieved, it is evident that the remediation efforts have not completely eliminated contaminant mass discharge and associated risk. Remaining contaminant mass contributing to the current mass discharge is hypothesized to comprise poorly accessible mass in the source zones, as well as aqueous (and sorbed) mass present in the extensive lower-permeability units located within and adjacent to the contaminant plume. The fate of these sources is an issue of critical import to the remediation of chlorinated-solvent contaminated sites, and development of methods to address these sources will be required to achieve successful long-term management of such sites and to ultimately transition them to closure.     

Field demonstration and evaluation of the Passive Flux Meter on a CAH groundwater plume

This study comprises the first application of the Passive Flux Meter (PFM) for the measurement of chlorinated aliphatic hydrocarbon (CAH) mass fluxes and Darcy water fluxes in groundwater at a European field site. The PFM was originally developed and applied to measurements near source zones. The focus of the PFM is extended from near source to plume zones. For this purpose, 48 PFMs of 1.4 m length were constructed and installed in eight different monitoring wells in the source and plume zone of a CAH-contaminated field site located in France. The PFMs were retrieved, sampled, and analyzed after 3 to 11 weeks of exposure time, depending on the expected contaminant flux. PFM evaluation criteria include analytical, technical, and practical aspects as well as conditions and applicability. PFM flux data were compared with so-called traditional soil and groundwater concentration data obtained using active sampling methods. The PFMs deliver reasonable results for source as well as plume zones. The limiting factor in the PFM applicability is the exposure time together with the groundwater flux. Measured groundwater velocities at the field site range from 2 to 41 cm/day. Measured contaminant flux data raise up to 13 g/m²/day for perchloroethylene in the plume zone. Calculated PFM flux averaged concentration data and traditional concentration data were of similar magnitude for most wells. However, both datasets need to be compared with reservation because of the different sampling nature and time. Two important issues are the PFM tracer loss during installation/extraction and the deviation of the groundwater flow field when passing the monitoring well and PFM. The demonstration of the PFM at a CAH-contaminated field site in Europe confirmed the efficiency of the flux measurement technique for source as well as plume zones. The PFM can be applied without concerns in monitoring wells with European standards. The acquired flux data are of great value for the purpose of site characterization and mass discharge modeling, and can be used in combination with traditional soil and groundwater sampling methods.     

Toward enhanced subsurface intervention methods using chaotic advection

Many intervention activities in the terrestrial subsurface involve the need to recover/emplace distributions of scalar quantities (e.g. dissolved phase concentrations or heat) from/in volumes of saturated porous media. These scalars can be targeted by pump-and-treat methods or by amendment technologies. Application examples include in-situ leaching for metals, recovery of dissolved contaminant plumes, or utilizing heat energy in geothermal reservoirs. While conventional pumping methods work reasonably well, costs associated with maintaining pumping schedules are high and improvements in efficiency would be welcome. In this paper we discuss how transient switching of the pressure at different wells can intimately control subsurface flow, generating a range of "programmed" flows with various beneficial characteristics. Some programs produce chaotic flows which accelerate mixing, while others create encapsulating flows which can isolate fluid zones for lengthy periods. In a simplified model of an aquifer subject to balanced pumping, chaotic flow topologies have been predicted theoretically and verified experimentally using Hele-Shaw cells. Here, a survey of the key characteristics of chaotic advection is presented. Mathematical methods are used to show how these characteristics may translate into practical situations involving regional flows and heterogeneity. The results are robust to perturbations, and withstand significant aquifer heterogeneity. It is proposed that chaotic advection may form the basis of new efficient technologies for groundwater interventions.     

Evaluation of groundwater flow patterns around a dual-screened groundwater circulation well

Dual-screened groundwater circulation wells (GCWs) can be used to remove contaminant mass and to mix reagents in situ. GCWs are so named because they force water in a circular pattern between injection and extraction screens. The radial extent, flux and direction of the effective flow of this circulation cell are difficult to measure or predict. The objective of this study is to develop a robust protocol for assessing GCW performance. To accomplish this, groundwater flow patterns surrounding a GCW are assessed using a suite of tools and data, including: hydraulic head, in situ flow velocity, measured hydraulic conductivity data from core samples, chemical tracer tests, contaminant distribution data, and numerical flow and transport models. The hydraulic head data show patterns that are consistent with pumping on a dual-screened well, however, many of the observed changes are smaller than expected. In situ thermal perturbation flow sensors successfully measured horizontal flow, but vertical flow could not be determined with sufficient accuracy to be useful in mapping flow patterns. Two types of chemical tracer tests were utilized at the site and showed that much of the flow occurs within a few meters of the GCW. Flow patterns were also assessed based on changes in contaminant (trichloroethylene, TCE) concentrations over time. The TCE data clearly showed treated water moving away from the GCW at shallow and intermediate depths, but the circulation of that water back to the well, except very close to the well, was less clear. Detailed vertical and horizontal hydraulic conductivities were measured on 0.3 m-long sections from a continuous core from the GCW installation borehole. The measured vertical and horizontal hydraulic conductivity data were used to construct numerical flow and transport models, the results of which were compared to the head, velocity and concentration data. Taken together, the field data and modeling present a fairly consistent picture of flow and transport around the GCW. However, the time and expense associated with conducting all of those tests would be prohibitive for most sites. As a consequence, a sequential protocol for GCW characterization is presented here in which the number of tools used can be adjusted to meet the needs of individual sites. While not perfect, we believe that this approach represents the most efficient means for evaluating GCW performance.     

Quantification of mass fluxes and natural attenuation rates at an industrial site with a limited monitoring network: a case study

At many "real world" field sites, the number of available monitoring wells is limited due to economic or geological reasons. Under such restricted conditions, it is difficult to perform a reliable field investigation and to quantify primary lines of evidence for natural attenuation (NA), like the documentation of a decrease of contaminant mass flux in flow direction. This study reports the results of a groundwater investigation at a former manufactured gas plant situated in a Quaternary river valley in southwest Germany. The location, infrastructure and aquifer setting are typical of many industrial sites in Germany. Due to difficult drilling conditions (coarse glaciofluvial gravel deposits and an anthropogenic fill above the aquifer), only 12 monitoring wells were available for the investigation and localisation of the contaminant plume. These wells were situated along three control planes (CP) downgradient from the contaminant source, with four wells along each plane. Based on the sparse set of monitoring wells, field scale mass fluxes and first-order natural attenuation rate constants of benzene, toluene, ethylbenzene, and o-xylene and p-xylene (BTEX) and low molecular weight polycyclic aromatic hydrocarbons (PAH) were estimated utilizing different point scale and also a new integral investigation method. The results show that even at a heterogeneous site with a sparse monitoring network point scale investigation methods can provide reliable information on field scale natural attenuation rates, if a dependable flow model or tracer test data is available. If this information is not available, only the new integral investigation method presented can yield adequate results for the quantification of contaminant mass fluxes under sparse monitoring conditions.