Simulating the fate and transport of TCE from groundwater to indoor air

This work provides an exploratory analysis on the relative importance of various factors controlling the fate and transport of volatile organic contaminants (in this case, TCE) from a DNAPL source zone located below the water table and into the indoor air. The analysis is conducted using the multi-phase compositional model CompFlow Bio, with the base scenario problem geometry reminiscent of a field experiment conducted by Rivett [Rivett, M.O., (1995), Soil–gas signatures from volatile chlorinated solvents: Borden field experiments. Groundwater, 33(1), 84–98.] at the Borden aquifer where groundwater was observed to transport a contaminant plume a substantial distance without vertical mass transport of the contaminant across the capillary fringe and into the vadose zone. Results for the base scenario model indicate that the structure of the permeability field was largely responsible for deflecting the groundwater plume upward towards the capillary fringe, permitting aqueous phase diffusion to transport the TCE into the vadose zone. Alternative permeability realizations, generated as part of a Monte Carlo simulation process, at times deflected the groundwater plume downwards causing the extended thickness of the saturated zone to insulate the vadose zone from exposure to the TCE by upward diffusive transport. Comparison of attenuation coefficients calculated using the CompFlow Bio and Johnson and Ettinger [Johnson, P.C. and Ettinger, R.A., (1991), Heuristic model for predicting the intrusion rate of contaminant vapors into buildings. Environmental Science and Technology, 25, 1445–1452.] heuristic model exhibited fortuitous agreement for the base scenario problem geometry, with this agreement diverging for the alternative permeability realizations as well as when parameters such as the foundation slab fracture aperture, the indoor air pressure drop, the capillary fringe thickness, and the infiltration rate were varied over typical ranges.     

Mathematical Modeling of Column and Field Dense Nonaqueous Phase Liquid Tracer Tests

Mathematical models for the simulation of dense nonaqueous phase liquid tracer tests (DTTs) in laboratory columns and in the field are developed and examined. The DTT technique is a means of estimating the quantity of dense nonaqueous phase liquid (DNAPL) in a domain of interest in an aquifer. The two-dimensional field DTT model uses the Method of Principal Directions and an asymmetrical upwind algorithm for describing advective transport. Both models include diffusion transport of tracer into and from low-permeability porous structures such as clay lenses, as well as the mass transport kinetics of partitioning tracer to and from the DNAPL droplets. The dependence of the effluent tracer concentration curves on the parameters of the models is explored, and conclusions are drawn regarding the applicability of, and several possible problems with, the DTT technique. Model results indicate that the DTT performs well at locating distributed droplets of DNAPL, but is unlikely to be useful in the assessment of pooled DNAPL.     

Self-inhibition can limit biologically enhanced TCE dissolution from a TCE DNAPL

Biodegradation of trichloroethene (TCE) near a Dense Non Aqueous Phase Liquid (DNAPL) can enhance the dissolution rate of the DNAPL by increasing the concentration gradient at the DNAPL-water interface. Two-dimensional flow-through sand boxes containing a TCE DNAPL and inoculated with a TCE dechlorinating consortium were set up to measure this bio-enhanced dissolution under anaerobic conditions. The total mass of TCE and daughter products in the effluent of the biotic boxes was 3-6 fold larger than in the effluent of the abiotic box. However, the mass of daughter products only accounted for 19-55% of the total mass of chlorinated compounds in the effluent, suggesting that bio-enhanced dissolution factors were maximally 1.3-2.2. The enhanced dissolution most likely primarily resulted from variable DNAPL distribution rather than biodegradation. Specific dechlorination rates previously determined in a stirred liquid medium were used in a reactive transport model to identify the rate limiting factors. The model adequately simulated the overall TCE degradation when predicted resident microbial numbers approached observed values and indicated an enhancement factor for TCE dissolution of 1.01. The model shows that dechlorination of TCE in the 2D box was limited due to the short residence time and the self-inhibition of the TCE degradation. A parameter sensitivity analysis predicts that the bio-enhanced dissolution factor for this TCE source zone can only exceed a value of 2 if the TCE self-inhibition is drastically reduced (when a TCE tolerant dehalogenating community is present) or if the DNAPL is located in a low-permeable layer with a small Darcy velocity.     

A three-layer diffusion-cell to examine bio-enhanced dissolution of chloroethene dense non-aqueous phase liquid

Microbial reductive dechlorination of trichloroethene (TCE) and perchloroethene (PCE) in the vicinity of their dense non-aqueous phase liquid (DNAPL) has been shown to accelerate DNAPL dissolution. A three-layer diffusion-cell was developed to quantify this bio-enhanced dissolution and to measure the conditions near the DNAPL interface. The 12 cm long diffusion-cell setup consists of a 5.5 cm central porous layer (sand), a lower 3.5 cm DNAPL layer and a top 3 cm water layer. The water layer is frequently refreshed to remove chloroethenes at the upper boundary of the porous layer, while the DNAPL layer maintains the saturated chloroethene concentration at the lower boundary. Two abiotic and two biotic diffusion-cells with TCE DNAPL were tested. In the abiotic diffusion-cells, a linear steady state TCE concentration profile between the DNAPL and the water layer developed beyond 21 d. In the biotic diffusion-cells, TCE was completely converted into cis-dichloroethene (cis-DCE) at 2.5 cm distance of the DNAPL. Dechlorination was likely inhibited up to a distance of 1.5 cm from the DNAPL, as in this part the TCE concentration exceeded the culture’s maximum tolerable concentration (2.5 mM). The DNAPL dissolution fluxes were calculated from the TCE concentration gradient, measured at the interface of the DNAPL layer and the porous layer. Biotic fluxes were a factor 2.4 (standard deviation 0.2) larger than abiotic dissolution fluxes. This diffusion-cell setup can be used to study the factors affecting the bio-enhanced dissolution of DNAPL and to assess bioaugmentation, pH buffer addition and donor delivery strategies for source zones.     

A multi-objective optimization framework for surfactant-enhanced remediation of DNAPL contaminations

The occurrence of Dense Non-Aqueous Phase Liquid (DNAPL) contaminations in the subsurface is a threat for drinkwater resources in the western world. Surfactant-Enhanced Aquifer Remediation (SEAR) is widely considered as one of the most promising techniques to remediate DNAPL contaminations in-situ, be it with considerable additional costs compared to classical pump-and-treat remediations. A cost-effective design of the remediation set-up is therefore essential. In this work, a pilot SEAR test is executed at a DNAPL contaminated site in Belgium in order to collect data for the calibration of a multi-phase multi-component model. The calibrated model is used to assess a series of scenario-analyses for the full-scale remediation of the site. The remediation variables that were varied were the injection and extraction rate, the injection and extraction duration, and the surfactant injection concentrations. A constrained multi-objective optimization of the model was applied to obtain a Pareto set of optimal remediation strategies with different weights for the two objectives of the remediation: (i) the maximal removal of DNAPL and (ii) a total minimal cost. These Pareto curves can help decision makers to select an optimal remediation strategy in terms of cost and remediation efficiency. The Pareto front shows a considerable trade-off between the total remediation cost and the removed DNAPL mass.     

Three-fluid retention in porous media involving water, PCE and air

A classical way to obtain three-fluid retention curves in porous media from measured two-fluid retention curves is based on the Leverett concept, which states that the total volumetric liquid content in a water-wet porous medium, containing water, a nonaqueous-phase liquid (NAPL) and air, is a function of the capillary pressure across the interface between the continuous NAPL and air. This functional relationship results from the assumed condition that in a three-fluid porous medium, the intermediate wetting fluid spreads over the water-air interface. Application of Leverett's concept may not be valid, however, for nonspreading NAPLs like perchloroethylene (PCE). This paper discusses measurements of both PCE-air and water-PCE-air retention curves using a long vertical column in conjunction with a dual-energy gamma radiation system. The data indicate that the Leverett concept was applicable only until a critical PCE saturation had been reached.     

Controlled release, blind test of DNAPL remediation by ethanol flushing

A dense nonaqueous phase liquid (DNAPL) source zone was established within a sheet-pile isolated cell through a controlled release of perchloroethylene (PCE) to evaluate DNAPL remediation by in-situ cosolvent flushing. Ethanol was used as the cosolvent, and the main remedial mechanism was enhanced dissolution based on the phase behavior of the water-ethanol-PCE system. Based on the knowledge of the actual PCE volume introduced into the cell, it was estimated that 83 L of PCE were present at the start of the test. Over a 40-day period, 64% of the PCE was removed by flushing the cell with an alcohol solution of approximately 70% ethanol and 30% water. High removal efficiencies at the end of the test indicated that more PCE could have been removed had it been possible to continue the demonstration. The ethanol solution extracted from the cell was recycled during the test using activated carbon and air stripping treatment. Both of these treatment processes were successful in removing PCE for recycling purposes, with minimal impact on the ethanol content in the treated fluids. Results from pre- and post-flushing partitioning tracer tests overestimated the treatment performance. However, both of these tracer tests missed significant amounts of the PCE present, likely due to inaccessibility of the PCE. The tracer results suggest that some PCE was inaccessible to the ethanol solution which led to the inefficient PCE removal rates observed. The flux-averaged aqueous PCE concentrations measured in the post-flushing tracer test were reduced by a factor of 3 to 4 in the extraction wells that showed the highest PCE removal compared to those concentrations in the pre-flushing tracer test.     

Equilibrium nonaqueous phase liquid pool geometry in coarse soils with discrete textural interfaces

This paper presents a model for the geometry of nonaqueous phase liquid (NAPL) pools and mounds in homogeneous soils and soils with discrete textural interfaces. It is shown that the concepts of capillary pressure-saturation curve hysteresis and entry pressures are integral to the complete conceptualization of pool and mound geometry. Unless hysteresis is included in the analysis, light NAPL (LNAPL) in homogeneous soils cannot exist in pools at all, and dense NAPL (DNAPL) will not mound on horizontal textural interfaces unless lateral confining boundaries are present. The proposed model also implies that remobilization of DNAPL pools will occur at lower hydraulic gradients than those predicted with previous models. Comparing predicted and experimental DNAPL and LNAPL pool thicknesses and the location of an LNAPL lens with respect to the top of the capillary fringe validate the model.     

Laboratory evidence of natural remobilization of multicomponent DNAPL pools due to dissolution

Mixtures of dense non-aqueous phase liquids (DNAPLs) trapped in the subsurface can act as long-term sources of contamination by dissolving into flowing groundwater. In general, the components of higher solubility are removed more quickly, thus altering the composition of the remaining DNAPL, and possibly leading to changes in its physical properties. Through the development of a simple compositional model, Roy et al. [J. Contam. Hydrol. 2002 (59) 163] showed that preferential dissolution of a mixed DNAPL could potentially result in changes in density and interfacial tension that could subsequently lead to remobilization of an initially static DNAPL pool. The laboratory experiments presented in this next paper provide a proof-of-concept for the previously presented theory, demonstrating and quantifying this process of remobilization. In addition, the experiments provide a data set for evaluation of the model presented by Roy et al. [J. Contam. Hydrol. 2002 (59) 163]. In the four experiments, a DNAPL pool comprised of tetrachloroethene and benzene was created as an open pool overlying glass beads within a water-saturated 2-D flow box. Experiments included rectangular and triangular pools. In each of the experiments, remobilization (as breakthrough) was observed more than 2 weeks after formation of the initial pool. During each experiment, the pool height declined as mass was lost by dissolution, while sampling indicated a decrease in the mole fraction of benzene, the more soluble component. Small protuberances formed along the bottom of the pool as its composition changed with time and the displacement pressure was achieved for various pore throats. Eventually one of the protuberances extended further, forming a finger (breakthrough). In general, the pool emptied as the finger proceeded further into the beads. It was also shown theoretically and experimentally that remobilization will occur sooner for pools with a triangular (pointing down), rather than rectangular, shape. The experimental results were simulated using the model developed by Roy et al. [J. Contam. Hydrol. 2002 (59) 163]. The model matched the observations well, suggesting that it accurately represents the primary mechanisms involved with natural remobilization under the conditions of the study.     

A theoretical model of air and steam co-injection to prevent the downward migration of DNAPLs during steam-enhanced extraction

When steam is injected into soil containing a dense volatile non-aqueous phase liquid contaminant the DNAPL vaporized within the heated soil region condenses and accumulates ahead of the steam condensation front. If enough DNAPL accumulates, gravitational forces can overcome trapping forces allowing the liquid contaminant to flow downward. By injecting air with steam, a portion of the DNAPL vapor remains suspended in equilibrium with the air, decreasing liquid contaminant accumulation ahead of the steam condensation front, and thus reducing the possibility of downward migration. In this work, a one-dimensional theoretical model is developed to predict the injection ratio of air to steam that will prevent the accumulation of volatile DNAPLs. The contaminated region is modeled as a one-dimensional homogeneous porous medium with an initially uniform distribution of a single component contaminant. Mass and energy balances are combined to determine the injection ratio of air to steam that eliminates accumulation of the contaminant ahead of the steam condensation front, and hence reduces the possibility of downward migration. The minimum injection ratio that eliminates accumulation is defined as the optimum injection ratio. Example calculations are presented for three DNAPLs, carbon tetrachloride (CCl4), trichloroethylene (TCE), and perchloroethylene (PCE). The optimum injection ratio of air to steam is shown to depend on the initial saturation and the volatility of the liquid contaminant. Numerical simulation results are presented to validate the model, and to illustrate downward migration for ratios less than optimum. Optimum injection ratios determined from numerical simulations are shown to be in good agreement with the theoretical model.