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.     

Relative velocities of DNAPL and aqueous phase plume migration

Numerical simulation is used to examine the relative velocities of DNAPL and aqueous phase plumes in sandy aquifers where lateral spreading of DNAPL has occurred at the base of the aquifer. The scenario being modeled is one where a permeable aquifer is underlain by a sloping aquitard, which results in lateral migration of the DNAPL down the slope, in addition to lateral migration of an aqueous phase plume subject to a specified hydraulic gradient. A sensitivity analysis is presented to the impacts of both DNAPL properties and geologic properties. The most important chemical properties governing the relative velocities of the DNAPL and the shallow aqueous phase plume are the DNAPL viscosity and the aqueous component soil-water partition coefficient (Kd). The dip of the underlying aquitard was found to be relatively unimportant, at least for the range of values studied. The scenario under consideration can be important in conceptual model development and remedial design, as in certain cases DNAPL could be migrating in areas without the evidence of a well-developed aqueous phase plume. The implication of this work is that the absence of a shallow aqueous phase plume directly downgradient of a DNAPL source zone does not rule out the possibility of deep occurrences of DNAPL beyond the shallow monitoring well network. A further finding of this study is that the occurrence of a highly sorbing compound in groundwater at virtually any concentration may indicate the immediate upgradient presence of residual or pooled DNAPL.     

Multiphase flow and transport in fractured clay/sand sequences.

A numerical model (Queen's University Multi-Phase Flow Simulator, QUMPFS) was used to assess the rate of trichloroethylene (TCE) dense, non-aqueous phase liquid (DNAPL) migration through fractured clay, with special attention focused on the influence of interbedded sand lenses. The presence of these sand lenses was found to increase the time required for the non-wetting phase to migrate through the full 30 m vertical extent of the clay sequence from a few days to several years. Applied vertical hydraulic gradients were found to be moderately influential in systems consisting solely of fractured clays, yet one of the dominant factors controlling speed of vertical migration when sand lenses were present. Larger displacement pressure of the sands relative to that of the fractures leads to slower DNAPL migration rates, due to the delays that occur during build-up of capillary pressures. Dissolution of DNAPL and subsequent matrix diffusion of the aqueous phase has little effect on the rate of DNAPL migration through systems consisting of fractured clay only, yet slows the rate of migration in systems containing sand lenses. In all cases examined, the rate of DNAPL loading to the lower aquifer far exceeded the rate of aqueous phase mass loading. It was also found that DNAPL reaches the lower aquifer at approximately the same time as the aqueous phase plumes even for systems experiencing downward groundwater flow due to the attenuation of the aqueous phase through matrix diffusion.     

Effect of source variability and transport processes on carbon isotope ratios of TCE and PCE in two sandy aquifers

Chlorinated ethenes often migrate over extended distances in aquifers and may originate from different sources. The aim of this study was to determine whether stable carbon isotope ratios remain constant during dissolution and transport of chlorinated ethenes and whether the ratios can be used to link plumes to their sources. Detailed depth-discrete delineation of the carbon isotope ratio in a tetrachloroethene (PCE) plume and in a trichloroethene (TCE) plume was done along cross-sections orthogonal to groundwater flow in two sandy aquifers in the Province of Ontario, Canada. At the TCE site, TCE concentrations up to solubility were measured in one high concentration zone close to the bottom of the aquifer from where dense non-aqueous phase liquid (DNAPL) was collected. A laboratory experiment using the DNAPL indicated that only very small carbon isotope fractionation occurs during dissolution of TCE (0.26 per thousand), which is consistent with field observations. At most sampling points, the delta(13)C of dissolved TCE was similar to that of the DNAPL except for a few sampling points at the bottom of the aquifer close to the underlying aquitard. At these points, a (13)C enrichment of up to 2.4 per thousand was observed, which was likely due to biodegradation and possibly preferential diffusion of TCE with (12)C into the aquitard. In contrast to the TCE site, several distinct zones of high concentration were observed at the PCE site and from zones to zone, the delta(13)C values varied substantially from -24.3 per thousand to -33.6 per thousand. Comparison of the delta(13)C values in the high concentration zones made it possible to divide the plume in the three different domains, each probably representing a different episode and location of DNAPL release. The three different zones could still be distinguished 220 m from the DNAPL sources. This demonstrates that carbon isotope ratios can be used to differentiate between different zones in chlorinated ethene plumes and to link plume zones to their sources. In addition, subtle variations in delta(13)C at plume fringes provided insight into mechanisms of plume spreading in transverse vertical direction. These variations were identified because of the high-resolution provided by the monitoring network.     

Co-injection of air and steam for the prevention of the downward migration of DNAPLs during steam enhanced extraction: an experimental evaluation of optimum injection ratio predictions

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 a previous work, a theoretical model was developed to predict the optimum injection ratio of air to steam that would eliminate accumulation of DNAPL ahead of the temperature front and thus minimize the potential for downward migration. In this work, the theoretical model is summarized, and an experiment is presented in order to evaluate the optimum injection ratio prediction. In the experiment, a two-dimensional water saturated sand pack is contaminated with a known mass of TCE (DNAPL). The system is then remediated by co-injecting air and steam at the predicted optimum injection ratio, calculated based on the average contaminant soil concentration in the sand pack. Results for the co-injection of air and steam are compared to results for the injection of pure steam or pure air. Injection at the predicted optimum injection ratio for a volumetric average NAPL saturation, reduced accumulation of the contaminant ahead of the condensation front by over 90%, as compared to steam injection alone. This indicates that the optimum injection ratio prediction is a valuable tool for limiting the spreading of DNAPL during steam-enhanced extraction. Injection at the optimum injection ratio resulted in earlier recovery of contaminant than for steam injection alone. Co-injection of steam and air is also shown to result in much higher recovery rates than air injection alone.     

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.     

A new concept: the use of neutrally-buoyant microemulsions for DNAPL remediation

Even in the absence of mobilization of dense nonaqueous phase liquid (DNAPL), the microemulsion that forms when the surfactant solubilizes a dense contaminant such as trichloroethylene will be more dense than water and tends to migrate downward. This paper addresses the issue of migration with a new concept: surfactant enhanced aquifer remediation at neutral buoyancy. Laboratory results of surfactant remediation in two-dimensional model aquifers show that downward migration of microemulsion containing solubilized dense contaminants can be reduced to an acceptable level, even in the absence of capillary barriers in the aquifer. One model experiment was designed to exhibit a small degree of vertical migration and full capture of the microemulsion at the extraction well. The second experiment was designed to demonstrate the effect of large buoyancy forces that lead to excessive downward migration of the microemulsion. Density measurements of aqueous solutions containing sodium dihexyl sulfosuccinate surfactant, isopropanol, trichloroethylene, and sodium chloride are presented. A companion paper presents the results of the flow and transport calculations needed for this approach to surfactant flooding.     

Removal of DNAPL contamination from the saturated zone by the combined effect of vertical upward flushing and density reduction

A common aspect of innovative remediation techniques is that they tend to reduce the interfacial tension between the aqueous and non-aqueous phase liquids, resulting in mobilization of the organic contaminant. This complicates the remediation of aquifers, contaminated with Dense Non-Aqueous Phase Liquids (DNAPLs), as they are likely to migrate downwards, deeper into the aquifer and into finer layers. A possible solution is the use of swelling alcohols, which tend to reduce the density difference between the aqueous phase and the DNAPL. To avoid premature mobilization upon the initial contact between the DNAPL and the alcohol, several researchers have proposed the use of vertical upward flow of the alcohol. In this paper, we present an equation, which describes the upward mobilization of both continuous and discontinuous DNAPLs and so the important parameters governing the upward controlled mobilization of the DNAPL. The need and required magnitude of this specific discharge was investigated by conducting four column experiments in which the initial density of the DNAPL and the permeability was varied. It was shown that the required flow velocities increase with the permeability of the porous medium and the initial density difference between the aqueous phase and the DNAPL. Whenever the specific discharge falls below the critical value, the DNAPL moves downward. A second set of column experiments looked at the impact of permeability of porous medium on the solubilization and mobilization of DNAPL during alcohol flooding. Columns, packed with coarse or fine sand, containing a residual trichloroethylene (TCE) or perchloroethylene (PCE) saturation were flushed with the alcohol mixture at a fixed specific discharge rate. The induced pressure gradients in the aqueous phase, which were higher in the fine sand, resulted for this porous medium in extensive mobilization of the DNAPL against the direction of the buoyancy force. The density of the first NAPL coming out of the top of the fine sand was close to that of the pure DNAPL. In the coarser sand, the pressure gradients were sufficient to prevent downward migration of the DNAPL, but upward mobilization was minimal. The predominant removal mechanism in this case was the much slower solubilization.     

Analytical solution of two-dimensional solute transport in an aquifer–aquitard system

This study deals with two-dimensional solute transport in an aquifer–aquitard system by maintaining rigorous mass conservation at the aquifer–aquitard interface. Advection, longitudinal dispersion, and transverse vertical dispersion are considered in the aquifer. Vertical advection and diffusion are considered in the aquitards. The first-type and the third-type boundary conditions are considered in the aquifer. This study differs from the commonly used averaged approximation (AA) method that treats the mass flux between the aquifer and aquitard as an averaged volumetric source/sink term in the governing equation of transport in the aquifer. Analytical solutions of concentrations in the aquitards and aquifer and mass transported between the aquifer and upper or lower aquitard are obtained in the Laplace domain, and are subsequently inverted numerically to yield results in the real time domain (the Zhan method). The breakthrough curves (BTCs) and distribution profiles in the aquifer obtained in this study are drastically different from those obtained using the AA method. Comparison of the numerical simulation using the model MT3DMS and the Zhan method indicates that the numerical result differs from that of the Zhan method for an asymmetric case when aquitard advections are at the same direction. The AA method overestimates the mass transported into the upper aquitard when an upward advection exists in the upper aquitard. The mass transported between the aquifer and the aquitard is sensitive to the aquitard Peclet number, but less sensitive to the aquitard diffusion coefficient.     

Evaluation of trichloroethene recovery processes in heterogeneous aquifer cells flushed with biodegradable surfactants

The ability of two biodegradable surfactants, polyoxyethylene (20) sorbitan monooleate (Tween 80) and sodium dihexyl sulfosuccinate (Aerosol MA), to recover a representative dense non-aqueous-phase liquid (DNAPL), trichloroethene (TCE), from heterogeneous porous media was evaluated through a combination of batch and aquifer cell experiments. An aqueous solution containing 3.3% Aerosol MA, 8% 2-propanol and 6 g/l CaCl(2) yielded a weight solubilization ratio (WSR) of 1.21 g TCE/g surfactant, with a corresponding liquid-liquid interfacial tension (IFT) of 0.19 dyn/cm. Flushing of aquifer cells containing a TCE-DNAPL source zone with approximately two pore volumes of the AMA formulation resulted in substantial (>30%) mobilization of TCE-DNAPL. However, a TCE mass recovery of 81% was achieved when the aqueous-phase flow rate was sufficient to displace the mobile TCE-DNAPL toward the effluent well. Aqueous solutions of Tween 80 exhibited a greater capacity to solubilize TCE (WSR=1.74 g TCE/g surfactant) and exerted markedly less reduction in IFT (10.4 dyn/cm). These data contradict an accepted empirical correlation used to estimate IFT values from solubilization capacity, and indicate a unique capacity of T80 to form concentrated TCE emulsions. Flushing of aquifer cells with less than 2.5 pore volumes of a 4% T80 solution achieved TCE mass recoveries ranging from 66 to 85%, with only slight TCE-DNAPL mobilization (<5%) occurring when the total trapping number exceeded 2 x 10(-5). These findings demonstrate the ability of Tween 80 and Aerosol MA solutions to efficiently recover TCE from a heterogeneous DNAPL source zone, and the utility of the total trapping number as a design parameter for a priori prediction of DNAPL mobilization and bank angle formation when flushing with low-IFT solutions. Given their potential to stimulate microbial reductive dechlorination at low concentrations, these surfactants are well-suited for remedial action plans that couple aggressive mass removal followed by enhanced bioremediation to treat chlorinated solvent source zones.     

Changes in contaminant mass discharge from DNAPL source mass depletion: evaluation at two field sites

Changes in contaminant fluxes resulting from aggressive remediation of dense nonaqueous phase liquid (DNAPL) source zone were investigated at two sites, one at Hill Air Force Base (AFB), Utah, and the other at Ft. Lewis Military Reservation, Washington. Passive Flux Meters (PFM) and a variation of the Integral Pumping Test (IPT) were used to measure fluxes in ten wells installed along a transect down-gradient of the trichloroethylene (TCE) source zone, and perpendicular to the mean groundwater flow direction. At both sites, groundwater and contaminant fluxes were measured before and after the source-zone treatment. The measured contaminant fluxes (J; ML(-2)T(-1)) were integrated across the well transect to estimate contaminant mass discharge (M(D); MT(-1)) from the source zone. Estimated M(D) before source treatment, based on both PFM and IPT methods, were approximately 76 g/day for TCE at the Hill AFB site; and approximately 640 g/day for TCE, and approximately 206 g/day for cis-dichloroethylene (DCE) at the Ft. Lewis site. TCE flux measurements made 1 year after source treatment at the Hill AFB site decreased to approximately 5 g/day. On the other hand, increased fluxes of DCE, a degradation byproduct of TCE, in tests subsequent to remediation at the Hill AFB site suggest enhanced microbial degradation after surfactant flooding. At the Ft. Lewis site, TCE mass discharge rates subsequent to remediation decreased to approximately 3 g/day for TCE and approximately 3 g/day for DCE approximately 1.8 years after remediation. At both field sites, PFM and IPT approaches provided comparable results for contaminant mass discharge rates, and show significant reductions (>90%) in TCE mass discharge as a result of DNAPL mass depletion from the source zone.     

Refinement of the density-modified displacement method for efficient treatment of tetrachloroethene source zones

A novel method to remediate dense nonaqueous phase liquid (DNAPL) source zones that incorporates in situ density conversion of DNAPL via alcohol partitioning followed by displacement with a low interfacial tension (IFT) surfactant flood has been developed. Previous studies demonstrated the ability of the density-modified displacement (DMD) method to recover chlorobenzene (CB) and trichloroethene (TCE) from heterogeneous porous media without downward migration of the dissolved plume or free product. However, the extent of alcohol (n-butanol) partitioning required for in situ density conversion of high-density NAPLs, such as tetrachloroethene (PCE), could limit the utility of the DMD method. Hence, the objective of this study was to compare the efficacy of two n-butanol delivery approaches: an aqueous solution of 6% (wt) n-butanol and a surfactant-stabilized macroemulsion containing 15% (vol) n-butanol in water, to achieve density reduction of PCE-NAPL in two-dimensional (2-D) aquifer cells. Results of liquid-liquid equilibrium studies indicated that density conversion of PCE relative to water occurred at an n-butanol mole fraction of 0.56, equivalent to approximately 5 ml n-butanol per 1 ml of PCE when in equilibrium with an aqueous solution. In 2-D aquifer cell studies, density conversion of PCE was realized using both n-butanol preflood solutions, with effluent NAPL samples exhibiting density reductions ranging from 0.51 to 0.70 g/ml. Although the overall PCE mass recoveries were similar (91% and 93%) regardless of the n-butanol delivery method, the surfactant-stabilized macroemulsion preflood removed approximately 50% of the PCE mass. In addition, only 1.2 pore volumes of the macroemulsion solution were required to achieve in situ density conversion of PCE, compared to 6.4 pore volumes of the 6% (wt) n-butanol solution. These findings demonstrate that use of the DMD method with a surfactant-stabilized macroemulsion containing n-butanol holds promise as an effective source zone remediation technology, allowing for efficient recovery of PCE-DNAPL while mitigating downward migration of the dissolved plume and free product.     

A controlled field experiment on groundwater contamination by a multicomponent DNAPL: creation of the emplaced-source and overview of dissolved plume development

A unique field experiment has been undertaken at the CFB Borden research site to investigate the development of dissolved chlorinated solvent plumes from a residual dense non-aqueous phase liquid (DNAPL) source. The "emplaced-source" tracer test methodology involved a controlled emplacement of a block-shaped source of sand containing chlorinated solvents below the water table. The gradual dissolution of this residual DNAPL solvent source under natural aquifer conditions caused dissolved solvent plumes of trichloromethane (TCM), trichloroethene (TCE) and perchloroethene (PCE) to continuously develop down gradient. Source dissolution and 3-D plume development were successfully monitored via 173 multilevel samplers over a 475-day tracer test period prior to site remediation research being initiated. Detailed groundwater level and hydraulic conductivity data were collected. Development of plumes with concentrations spanning 1-700,000 micrograms/1 is described and key processes controlling their migration identified. Plumes were observed to be narrow due to the weakness of transverse dispersion processes and long due to advection and significant longitudinal dispersion, very limited sorptive retardation and negligible, if any, attenuation due to biodegradation or abiotic reaction. TCM was shown to be essentially conservative, TCE very nearly conservative and PCE, consistent with its greater hydrophobicity, more retarded yet having a greater mobility than observed in previous Borden field tests. The absence of biodegradation was ascribed to the prevailing aerobic conditions and lack of any additional biodegradable carbon substrates. The transient groundwater flow regime caused significant transverse lateral plume movement, plume asymmetry and was likely responsible for most of the, albeit limited, transverse horizontal plume spreading. In agreement with the widespread incidence of extensive TCE and PCE plumes throughout the industrialized world, the experiment indicates such solvent plumes are likely to be highly mobile and persistent, at least in aquifers that are aerobic and have low sorption potential (low foc content).     

Dip-angle influence on areal DNAPL recovery by co-solvent flooding with and without pre-flooding

A two-dimensional (2D) laboratory model was used to study effects of gravity on areal recovery of a representative dense non-aqueous phase liquid (DNAPL) contaminant by an alcohol pre-flood and co-solvent flood in dipping aquifers. Recent studies have demonstrated that injection of alcohol and co-solvent solutions can be used to reduce in-situ the density of DNAPL globules and displace the contaminant from the source zone. However, contact with aqueous alcohol reduces interfacial tension and causes DNAPL swelling, thus facilitating risk of uncontrolled downward DNAPL migration. The 2D laboratory model was operated with constant background gradient flow and a DNAPL spill was simulated using tetrachloroethene (PCE). The spill was dispersed to a trapped, immobile PCE saturation by a water flood. Areal PCE recovery was studied using a double-triangle well pattern to simulate a remediation scheme consisting of an alcohol pre-flood using aqueous isobutanol ( approximately 10% vol.) followed by a co-solvent flood using a solution of ethylene glycol (65%) and 1-propanol (35%). Experiments were conducted with the 2D model oriented in the horizontal plane and compared to experiments at the 15 degrees and 30 degrees dip-angle orientations. Injection was applied either in the downward or upward direction of flow. Experimental results were compared to theoretical predictions for flood front stability and used to evaluate effects of gravity on areal PCE recovery. Sensitivity experiments were performed to evaluate effects of the alcohol pre-flood on PCE areal recovery. For experiments conducted with the alcohol pre-flood and the 2D model oriented in the horizontal plane, results indicate that 89-93% of source zone PCE was recovered. With injection oriented downward, results indicate that areal PCE recovery was 70-77% for a 15 degrees dip angle and 57-59% for a 30 degrees dip angle. With injection oriented upward, results indicate that areal PCE recovery was 57-60% at the 30 degrees dip angle, which was similar to PCE recovery for injection in the downward flow direction. Lower areal PCE recovery at greater dip angles in either direction of flow was attributed to DNAPL swelling and migration, flood front instabilities and bypassing of the displaced fluid past the extraction wells during the alcohol pre-flood. Additional results demonstrate that the use of an alcohol pre-flood can be beneficial in improving DNAPL recovery in the horizontal orientation, but pre-flooding may reduce areal recovery efficiency in dip-angle orientations. This study also demonstrates the use of theoretical perturbation (fingering) analysis in predicting NAPL recovery efficiency for flooding processes in remediating aquifers with dip angles.     

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.     

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.     

Influence of mass transfer characteristics for DNAPL source depletion and contaminant flux in a highly characterized glaciofluvial aquifer

The transfer of contaminant mass between the nonaqueous- and aqueous-phases is a process of central importance for the remediation of sites contaminated by dense nonaqueous-phase liquids (DNAPLs). This paper describes a comparison of the results obtained with various alternative DNAPL-aqueous-phase mass transfer models contained in the literature for predicting DNAPL source-zone depletion times in groundwater systems. These dissolution models were largely developed through laboratory column experiments. To gain insight into the implications of various representations of the local-scale kinetic as well as equilibrium DNAPL dissolution processes, aquifer heterogeneity and the complex architecture of a DNAPL source-zone, the aqueous-phase contaminant concentrations and mass fluxes arriving at a down-gradient compliance boundary are analyzed in a conditional stochastic framework. The hydrogeologic setting is a heterogeneous fluvial aquifer in Southwest Germany, referred to as the aquifer analog dataset, that was intensively characterized in three dimensions for hydrogeological parameters that include permeability, effective porosity, grain size, mineralogy and sorption coefficients. By embedding the various dissolution models into the compositional, multiphase flow model, CompFlow, the relative times predicted for complete depletion of a released DNAPL source due to natural dissolution are explored. Issues related to achieving environmental benefits through, for example, partial DNAPL-zone source removal via enhanced remedial technologies are also discussed. In this context, performance metrics in the form of peak aqueous-phase contaminant concentrations and mass fluxes arriving at a down-gradient compliance boundary are compared to each other. This is done for each of the alternative mass transfer models. A significant reduction in the fractional flux at a downstream location from the DNAPL source can be achieved by partial source-zone mass reduction; however, peak concentration levels at the same location remain much higher than the United States Environment Protection Agency (US-EPA) drinking water limits. Although groundwater quality was found to improve more rapidly for the equilibrium dissolution model, it is also shown that dissolution models that promote rapid DNAPL disappearance produce greater prediction uncertainty in the aqueous-phase flux reduction.     

The significance of heterogeneity on mass flux from DNAPL source zones: an experimental investigation

Understanding the process of mass transfer from source zones of aquifers contaminated with organic chemicals in the form of dense non-aqueous phase liquids (DNAPL) is of importance in site management and remediation. A series of intermediate-scale tank experiments was conducted to examine the influence of aquifer heterogeneity on DNAPL mass transfer contributing to dissolved mass emission from source zone into groundwater under natural flow before and after remediation. A Tetrachloroethylene (PCE) spill was performed into six source zone models of increasing heterogeneity, and both the spatial distribution of the dissolution behavior and the net effluent mass flux were examined. Experimentally created initial PCE entrapment architecture resulting from the PCE migration was largely influenced by the coarser sand lenses and the PCE occupied between 30 and 60% of the model aquifer depth. The presence of DNAPL had no apparent effect on the bulk hydraulic conductivity of the porous media. Up to 71% of PCE mass in each of the tested source zone was removed during a series of surfactant flushes, with associated induced PCE mobilization responsible for increasing vertical DNAPL distributions. Effluent mass flux due to water dissolution was also found to increase progressively due to the increase in NAPL-water contact area even though the PCE mass was reduced. Doubling of local groundwater flow velocities showed negligible rate-limited effects at the scale of these experiments. Thus, mass transfer behavior was directly controlled by the morphology of DNAPL within each source zone. Effluent mass flux values were normalized by the up-gradient DNAPL distributions. For the suite of aquifer heterogeneities and all remedial stages, normalized flux values fell within a narrow band with mean of 0.39 and showed insensitivity to average source zone saturations.     

Three-dimensional density-dependent flow and multicomponent reactive transport modeling of chlorinated solvent oxidation by potassium permanganate

A popular method for the treatment of aquifers contaminated with chlorinated solvents is chemical oxidation based on the injection of potassium permanganate (KMnO₄). Both the high density (1025 gL^− 1) and reactivity of the treatment solution influence the fate of permanganate (MnO₄) in the subsurface and affect the degree of contaminant treatment. The MIN3P multicomponent reactive transport code was enhanced to simulate permanganate-based remediation, to evaluate the pathways of MnO₄ utilization, and to assess the role of density contrasts for the delivery of the treatment solution. The modified code (MIN3P-D) provides a direct coupling between density-dependent fluid flow, solute transport, contaminant treatment, and geochemical reactions. The model is used to simulate a field trial of TCE oxidation in a sandy aquifer that is underlain by an aquitard. Three-dimensional simulations are conducted for a coupled reactive system comprised of ten aqueous components, two mineral phases, TCE (dissolved, adsorbed, and NAPL), reactive organic matter, and including ion exchange reactions. Model parameters are constrained by literature data and a detailed data set from the field site under investigation. The general spatial and transient evolution in observed concentrations of the oxidant, dissolved TCE, and reaction products are adequately reproduced by the simulations. The model elucidates the important role of density-induced flow and transport on the distribution of the treatment solution into NAPL containing regions located at the aquifer–aquitard interface. Model results further suggest that reactions that do not directly affect the stability of MnO₄ have a negligible effect on solution density and MnO₄ delivery.     

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.