Heuristic space–time design of monitoring wells for contaminant plume characterization in stochastic flow fields

An optimization methodology for designing groundwater quality monitoring networks applicable to stochastic flow fields is presented and evaluated. The approach sets itself apart from previous techniques by incorporating the time dimension directly into the objective function. This function is extremized using a directed partial enumeration strategy guided by physical considerations related to transport processes. The result is a set of monitoring well locations and a sampling schedule that minimizes plume characterization error while satisfying constraints on the maximum number of wells and allowable number of active wells. The method is evaluated using hypothetical plumes with varying degrees of heterogeneity. Results indicate that the proposed approach is successful in generating near-optimal sampling networks that satisfy all imposed constraints. Monitoring networks with as little as three active wells and a total of 12 wells are found to provide adequate plume characterization for low toxicity contaminants.     

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.     

Optimal design of pump-and-treat systems under uncertain hydraulic conductivity and plume distribution

In this work, we present a stochastic optimal control framework for assisting the management of the cleanup by pump-and-treat of polluted shallow aquifers. In the problem being investigated, hydraulic conductivity distribution and dissolved contaminant plume location are considered as the uncertain variables. The framework considers the subdivision of the cleanup horizon in a number of stress periods over which the pumping policy implemented until that stage is dynamically adjusted based upon new information that has become available in the previous stages. In particular, by following a geostatistical approach, we study the idea of monitoring the cumulative contaminant mass extracted from the installed recovery wells, and using these measurements to generate conditional realizations of the hydraulic conductivity field. These realizations are thus used to obtain a more accurate evaluation of the initial plume distribution, and modify accordingly the design of the pump-and-treat system for the remainder of the remedial process. The study indicates that measurements of contaminant mass extracted from pumping wells retain valuable information about the plume location and the spatial heterogeneity characterizing the hydraulic conductivity field. However, such an information may prove quite soft, particularly in the instances where recovery wells are installed in regions where contaminant concentration is low or zero. On the other hand, integrated solute mass measurements may effectively allow for reducing parameter uncertainty and identifying the plume distribution if more recovery wells are available, in particular in the early stages of the cleanup process.     

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.     

Assessing the natural attenuation of organic contaminants in aquifers using plume-scale electron and carbon balances: model development with analysis of uncertainty and parameter sensitivity

A quantitative methodology is described for the field-scale performance assessment of natural attenuation using plume-scale electron and carbon balances. This provides a practical framework for the calculation of global mass balances for contaminant plumes, using mass inputs from the plume source, background groundwater and plume residuals in a simplified box model. Biodegradation processes and reactions included in the analysis are identified from electron acceptors, electron donors and degradation products present in these inputs. Parameter values used in the model are obtained from data acquired during typical site investigation and groundwater monitoring studies for natural attenuation schemes. The approach is evaluated for a UK Permo-Triassic Sandstone aquifer contaminated with a plume of phenolic compounds. Uncertainty in the model predictions and sensitivity to parameter values was assessed by probabilistic modelling using Monte Carlo methods. Sensitivity analyses were compared for different input parameter probability distributions and a base case using fixed parameter values, using an identical conceptual model and data set. Results show that consumption of oxidants by biodegradation is approximately balanced by the production of CH4 and total dissolved inorganic carbon (TDIC) which is conserved in the plume. Under this condition, either the plume electron or carbon balance can be used to determine contaminant mass loss, which is equivalent to only 4% of the estimated source term. This corresponds to a first order, plume-averaged, half-life of > 800 years. The electron balance is particularly sensitive to uncertainty in the source term and dispersive inputs. Reliable historical information on contaminant spillages and detailed site investigation are necessary to accurately characterise the source term. The dispersive influx is sensitive to variability in the plume mixing zone width. Consumption of aqueous oxidants greatly exceeds that of mineral oxidants in the plume, but electron acceptor supply is insufficient to meet the electron donor demand and the plume will grow. The aquifer potential for degradation of these contaminants is limited by high contaminant concentrations and the supply of bioavailable electron acceptors. Natural attenuation will increase only after increased transport and dilution.     

Assessing impacts of partial mass depletion in DNAPL source zones: II. Coupling source strength functions to plume evolution

Analytical solutions, describing the time-dependent DNAPL source-zone mass and contaminant discharge rate, derived previously in Part I [Falta, R.W., Rao, P.S., Basu, N., this issue. Assessing the impacts of partial mass depletion in DNAPL source zones: I. Analytical modeling of source strength functions and plume response. J. Contam. Hydrol.] are used as a flux-boundary condition in a semi-analytical contaminant transport model. These analytical solutions assume a power relationship between the flow-averaged source concentration, and the source DNAPL mass; the empirical exponent (gamma) is a function of the flow field heterogeneity, DNAPL architecture, and the correlation between them. The DNAPL source strength terms can account for partial source remediation, either at time zero, or at some later time after the DNAPL release. The transport model considers advection, retardation, three-dimensional dispersion, and sequential first-order decay/production of several species. A separate solution is used to compute the time-dependent mass of each contaminant in the plume. A series of examples using different values of gamma shows how the benefits of partial DNAPL source remediation can vary with site conditions. In general, when gamma>1, relatively large short-term reductions in the plume concentrations and mass occur, but the source longevity is not strongly affected. Conversely, when gamma<1, the short-term reductions in the plume concentrations and mass are smaller, but the source longevity can be greatly reduced. In either case, the source remediation effort is much more effective if it is undertaken at an early time, before much contaminant mass has entered the plume. If the remediation effort is significantly delayed, the leading parts of the plume are not affected by the source remediation, and additional control or remediation of the plume itself is required.     

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.     

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.     

Optimization-Based Atmospheric Plume Source Identification

This work applies optimization and an Eulerian inversion approach presented by Bagtzoglou and Baun in 2005 in order to reconstruct contaminant plume time histories and to identify the likely source of atmospheric contamination using data from a real test site for the first time. Present-day distribution of an atmospheric contaminant plume as well as data points reflecting the plume history allow the reconstruction and provide the plume velocity, distribution, and probable source. The method was tested to a hypothetical case and with data from the Forest Atmosphere Transfer and Storage (FACTS) experiment in the Duke experimental forest site. In the scenarios presented herein, as well as in numerous cases tested for verification purposes, the model conserved mass, successfully located the peak of the plume, and managed to capture the motion of the plume well but underestimated the contaminant peak.     

Theoretical analysis of deviations from local equilibrium during sorbing solute transport through idealized stratified aquifers

This paper studies the spreading characteristics of reactive solute plumes in idealized stratified aquifers. The aquifer consists of two layers having different permeabilities with flow parallel to the stratification. The solute is assumed to adsorb onto the aquifer solids according to a first-order reversible kinetic rate law; the adsorption parameters are spatially uniform. We use the Aris moment method to examine analytically the time evolution of the lower-order spatial moments of the depth-averaged contaminant plume for an instantaneous input of mass. The results demonstrate that sorption kinetics cause the total dissolved mass and average velocity of the contaminant plume to decrease with increasing travel time. The plume variance is shown to depend upon three factors: intra-layer longitudinal dispersion, intra-layer kinetics, and vertical averaging. The results indicate that the relative importance of sorption kinetics diminishes as the permeability contrast between the layers increases. We present a simple criterion that can be used to assess the applicability of the local equilibrium assumption in idealized stratified systems.     

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.     

Reliability of different sampling densities for estimating and mapping lichen diversity in biomonitoring studies

Sampling requirements related to lichen biomonitoring include optimal sampling density for obtaining precise and unbiased estimates of population parameters and maps of known reliability. Two available datasets on a sub-national scale in Italy were used to determine a cost-effective sampling density to be adopted in medium-to-large-scale biomonitoring studies. As expected, the relative error in the mean Lichen Biodiversity (Italian acronym: BL) values and the error associated with the interpolation of BL values for (unmeasured) grid cells increased as the sampling density decreased. However, the increase in size of the error was not linear and even a considerable reduction (up to 50%) in the original sampling effort led to a far smaller increase in errors in the mean estimates (<6%) and in mapping (<18%) as compared with the original sampling densities. A reduction in the sampling effort can result in considerable savings of resources, which can then be used for a more detailed investigation of potentially problematic areas. It is, however, necessary to decide the acceptable level of precision at the design stage of the investigation, so as to select the proper sampling density.     

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.     

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.     

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.     

Comparing estimates of fugitive landfill methane emissions using inverse plume modeling obtained with Surface Emission Monitoring (SEM), Drone Emission Monitoring (DEM), and Downwind Plume Emission Monitoring (DWPEM)

ABSTRACT

As part of the global effort to quantify and manage anthropogenic greenhouse gas emissions, there is considerable interest in quantifying methane emissions in municipal solid waste landfills. A variety of analytical and experimental methods are currently in use for this task. In this paper, an optimization-based estimation method is employed to assess fugitive landfill methane emissions. The method combines inverse plume modeling with ambient air methane concentration measurements. Three different measurement approaches are tested and compared. The method is combined with surface emission monitoring (SEM), above ground drone emission monitoring (DEM), and downwind plume emission monitoring (DWPEM). The methodology is first trialed and validated using synthetic datasets in a hand-generated case study. A field study is also presented where SEM, DEM and DWPEM are tested and compared. Methane flux during two-days measurement campaign was estimated to be between 228 and 350 g/s depending on the type of measurements used. Compared to SEM, using unmanned aerial systems (UAS) allows for a rapid and comprehensive coverage of the site. However, as showed through this work, advancement of DEM-based methane sampling is governed by the advances that could be made in UAS-compatible measurement instrumentations. Downwind plume emission monitoring led to a smaller estimated flux compared with SEM and DEM without information about positions of major leak points in the landfill. Even though, the method is simple and rapid for landfill methane screening. Finally, the optimization-based methodology originally developed for SEM, shows promising results when it is combined with the drone-based collected data and downwind concentration measurements. The studied cases also discovered the limitations of the studied sampling strategies which is exploited to identify improvement strategies and recommendations for a more efficient assessment of fugitive landfill methane emissions.

Implications: Fugitive landfill methane emission estimation is tackled in the present study. An optimization-based method combined with inverse plume modeling is employed to treat data from surface emission monitoring, drone-based emission monitoring and downwind plume emission monitoring. The study helped revealing the advantages and the limitations of the studied sampling strategies. Recommendations for an efficient assessment of landfill methane emissions are formulated. The method trialed in this study for fugitive landfill methane emission could also be appropriate for rapid screening of analogous greenhouse gas emission hotspots.     

A new mixed segment-puff approach for dispersion modeling

This paper presents a new mixed methodology for realistic and cost-effective simulation of shortterm air quality dispersion phenomena using the Gaussian formula. The method can be applied to shortrange, intermediate and, especially, long-range transport simulations. Pollutant dynamics are described by the temporal evolution of plume elements, treated as segments or puffs according to their size. While the segments provide a numerically fast simulation during transport conditions, the puffs allow a proper simulation of calm or low-wind situations.The methodology is incorporated into a computer package (AVACTA II, Release 3) that gives the user large flexibility in defining the computational domain, the three-dimensional meteorological and emission input, the receptor locations, and in selecting plume rise and sigma formulas. AVACTA II provides both pollutant concentration fields and dry/wet deposition patterns. The model uses linear chemistry and is applicable to any two-species reaction chain (e.g., SO₂ and SO²⁻₄) where this approximation is reasonable and an appropriate reaction rate is available.     

Reactive transport modeling of processes controlling the distribution and natural attenuation of phenolic compounds in a deep sandstone aquifer

Reactive solute transport modeling was utilized to evaluate the potential for natural attenuation of a contaminant plume containing phenolic compounds at a chemical producer in the West Midlands, UK. The reactive transport simulations consider microbially mediated biodegradation of the phenolic compounds (phenols, cresols, and xylenols) by multiple electron acceptors. Inorganic reactions including hydrolysis, aqueous complexation, dissolution of primary minerals, formation of secondary mineral phases, and ion exchange are considered. One-dimensional (1D) and three-dimensional (3D) simulations were conducted. Mass balance calculations indicate that biodegradation in the saturated zone has degraded approximately 1-5% of the organic contaminant plume over a time period of 47 years. Simulations indicate that denitrification is the most significant degradation process, accounting for approximately 50% of the organic contaminant removal, followed by sulfate reduction and fermentation reactions, each contributing 15-20%. Aerobic respiration accounts for less than 10% of the observed contaminant removal in the saturated zone. Although concentrations of Fe(III) and Mn(IV) mineral phases are high in the aquifer sediment, reductive dissolution is limited, producing only 5% of the observed mass loss. Mass balance calculations suggest that no more than 20-25% of the observed total inorganic carbon (TIC) was generated from biodegradation reactions in the saturated zone. Simulations indicate that aerobic biodegradation in the unsaturated zone, before the contaminant entered the aquifer, may have produced the majority of the TIC observed in the plume. Because long-term degradation is limited to processes within the saturated zone, use of observed TIC concentrations to predict the future natural attenuation may overestimate contaminant degradation by a factor of 4-5.     

Simultaneous optimization of dense non-aqueous phase liquid (DNAPL) source and contaminant plume remediation

A framework is developed for simultaneous, optimal design of groundwater contaminant source removal and plume remediation strategies. The framework allows for varying degrees of effort and cost to be dedicated to source removal versus plume remediation. We have accounted for the presence of physical heterogeneity in the DNAPL source, since source heterogeneity controls mass release into the plume and the efficiency of source removal efforts. We considered high and low estimates of capital and operating costs for chemical flushing removal of the source, since these are expected to vary form site to site. Using the lower chemical flushing cost estimates, it is found that the optimal allocation of funds to source removal or plume remediation is sensitive to the degree of heterogeneity in the source. When the time elapsed between the source release and the implementation of remediation was varied, it was found that, except for the longest elapsed time (50,000 days), a combination of partial source removal and plume remediation was most efficient. When first-order, dissolved contaminant degradation was allowed, source removal was found to be unnecessary for the cases where the degradation rate exceeded intermediate values of the first-order rate constant. Finally, it was found that source removal became more necessary as the degree of aquifer heterogeneity increased.     

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.