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.     

PCE dissolution and simultaneous dechlorination by nanoscale zero-valent iron particles in a DNAPL source zone

While the capability of nanoscale zero-valent iron (NZVI) to dechlorinate organic compounds in aqueous solutions has been demonstrated, the ability of NZVI to remove dense non-aqueous phase liquid (DNAPL) from source zones under flow-through conditions similar to a field scale application has not yet been thoroughly investigated. To gain insight on simultaneous DNAPL dissolution and NZVI-mediated dechlorination reactions after direct placement of NZVI into a DNAPL source zone, a combined experimental and modeling study was performed. First, a DNAPL tetrachloroethene (PCE) source zone with emplaced NZVI was built inside a small custom-made flow cell and the effluent PCE and dechlorination byproducts were monitored over time. Second, a model for rate-limited DNAPL dissolution and NZVI-mediated dechlorination of PCE to its three main reaction byproducts with a possibility for partitioning of these byproducts back into the DNAPL was formulated. The coupled processes occurring in the flow cell were simulated and analyzed using a detailed three-dimensional numerical model. It was found that subsurface emplacement of NZVI did not markedly accelerate DNAPL dissolution or the DNAPL mass-depletion rate, when NZVI at a particle concentration of 10g/L was directly emplaced in the DNAPL source zone. To react with NZVI the DNAPL PCE must first dissolve into the groundwater and the rate of dissolution controls the longevity of the DNAPL source. The modeling study further indicated that faster reacting particles would decrease aqueous contaminant concentrations but there is a limit to how much the mass removal rate can be increased by increasing the dechlorination reaction rate. To ensure reduction of aqueous contaminant concentrations, remediation of DNAPL contaminants with NZVI should include emplacement in a capture zone down-gradient of the DNAPL source.     

Surfactant enhanced recovery of tetrachloroethylene from a porous medium containing low permeability lenses. 1. Experimental studies

A matrix of batch, column and two-dimensional (2-D) box experiments was conducted to investigate the coupled effects of rate-limited solubilization and layering on the entrapment and subsequent recovery of a representative dense NAPL, tetrachloroethylene (PCE), during surfactant flushing. Batch experiments were performed to determine the equilibrium solubilization capacity of the surfactant, polyoxyethylene (20) sorbitan monooleate (Tween 80), and to measure fluid viscosity, density and interfacial tension. Results of one-dimensional column studies indicated that micellar solubilization of residual PCE was rate-limited at Darcy velocities ranging from 0.8 to 8.2 cm/h and during periods of flow interruption. Effluent concentration data were used to develop effective mass transfer coefficient (Ke) expressions that were dependent upon the Darcy velocity and duration of flow interruption. To simulate subsurface heterogeneity, 2-D boxes were packed with layers of F-70 Ottawa sand and Wurtsmith aquifer material within 20-30 mesh Ottawa sand. A 4% Tween 80 solution was then flushed through PCE-contaminated boxes at several flow velocities, with periods of flow interruption. Effluent concentration data and visual observations indicated that both rate-limited solubilization and pooling of PCE above the fine layers reduced PCE recovery to levels below those anticipated from batch and column measurements. These experimental results demonstrate the potential impact of both mass transfer limitations and subsurface layering on the recovery of PCE during surfactant enhanced aquifer remediation.     

Modeling field-scale cosolvent flooding for DNAPL source zone remediation

A three-dimensional, compositional, multiphase flow simulator was used to model a field-scale test of DNAPL removal by cosolvent flooding. The DNAPL at this site was tetrachloroethylene (PCE), and the flooding solution was an ethanol/water mixture, with up to 95% ethanol. The numerical model, UTCHEM accounts for the equilibrium phase behavior and multiphase flow of a ternary ethanol-PCE-water system. Simulations of enhanced cosolvent flooding using a kinetic interphase mass transfer approach show that when a very high concentration of alcohol is injected, the DNAPL/water/alcohol mixture forms a single phase and local mass transfer limitations become irrelevant. The field simulations were carried out in three steps. At the first level, a simple uncalibrated layered model is developed. This model is capable of roughly reproducing the production well concentrations of alcohol, but not of PCE. A more refined (but uncalibrated) permeability model is able to accurately simulate the breakthrough concentrations of injected alcohol from the production wells, but is unable to accurately predict the PCE removal. The final model uses a calibration of the initial PCE distribution to get good matches with the PCE effluent curves from the extraction wells. It is evident that the effectiveness of DNAPL source zone remediation is mainly affected by characteristics of the spatial heterogeneity of porous media and the variable (and unknown) DNAPL distribution. The inherent uncertainty in the DNAPL distribution at real field sites means that some form of calibration of the initial contaminant distribution will almost always be required to match contaminant effluent breakthrough curves.     

Intermediate-scale 2D experimental investigation of in situ chemical oxidation using potassium permanganate for remediation of complex DNAPL source zones

In situ chemical oxidation is a technology that has been applied to speed up remediation of a contaminant source zone by inducing increased mass transfer from DNAPL sources into the aqueous phase for subsequent destruction. The DNAPL source zone can consist of one or more individual sources that may be present as an interconnected pool of high saturation, as a region of disconnected ganglia at residual saturation, or as combinations of these two morphologies. Potassium permanganate (KMnO(4)) is a commonly employed oxidant that has been shown to rapidly destroy DNAPL compounds like PCE and TCE following second-order kinetics in an aqueous system. During the oxidation of a target DNAPL compound, or naturally occurring reduced species in the subsurface, manganese oxide (MnO(2)) solids are produced. Research has shown that these manganese oxide solids may result in permeability reductions in the porous media thus reducing the ability for oxidant to be transported to individual DNAPL sources. It can also occur at the DNAPL-water interface, decreasing contact of the oxidant with the DNAPL. Additionally, MnO(2) formation at the DNAPL-water interface, and/or flow-bypassing as a result of permeability reductions around the source, may alter the mass transfer from the DNAPL into the aqueous phase, potentially diminishing the magnitude of any DNAPL mass depletion rate increase induced by oxidation. An experiment was performed in a two-dimensional (2D) sand-filled tank that included several discrete DNAPL source zones. Spatial and temporal monitoring of aqueous PCE, chloride, and permanganate concentrations was used to relate changes in mass depletion of, and mass flux, from DNAPL residual and pool source zones to chemical oxidation performance and MnO(2) formation. During the experiment, permeability changes were monitored throughout the 2D tank and these were related to MnO(2) deposition as measured through post-oxidation soil coring. Under the conditions of this experiment, MnO(2) formation was found to reduce permeability in and around DNAPL source zones resulting in changes to the overall flow pattern, with the effects depending on source zone configuration. A pool with little or no residual around it, in a relatively homogeneous flow field, appeared to benefit from resulting MnO(2) pore-blocking that substantially reduced mass transfer from the pool even though there was relatively little PCE mass removed from the pool. In contrast, a pool with residual around it (in a more typical heterogeneous flow field) appeared to undergo increased mass transfer as MnO(2) reduced permeability, altering the water flow and increasing the mixing at the DNAPL-water interface. Further, the magnitude of increased PCE mass depletion during oxidation appeared to depend on the PCE source configuration (pool versus ganglia) and decreased as MnO(2) was formed and deposited at the DNAPL-water interface. Overall, the oxidation of PCE mass appeared to be rate-limited by the mass transfer from the DNAPL to aqueous phase.     

Variations in expression of carbon isotope fractionation of chlorinated ethenes during biologically enhanced PCE dissolution close to a source zone

The stable carbon isotope values of tetrachloroethene (PCE) and its degradation products were monitored during studies of biologically enhanced dissolution of PCE dense nonaqueous phase liquid (DNAPL) to determine the effect of PCE dissolution on observed isotope values. The degradation of PCE was monitored in a 2-dimensional model aquifer and in a pilot test cell (PTC) at Dover Air Force Base, both with emplaced PCE DNAPL sources. Within the plume down gradient from the source, the isotopic fractionation of dissolved PCE and its degradation products were consistent with those observed in biodegradation laboratory studies. However, close to the source zone significant shifts in the isotope values of dissolved PCE were not observed in either the model aquifer or PTC due to the constant input of newly dissolved, non fractionated PCE, and the small isotopic fractionation associated with PCE reductive dechlorination by the mixed microbial culture used. Therefore the identification of reductive dechlorination in the presence of PCE DNAPL was based upon the appearance of daughter products and the isotope values of those daughter products. An isotope model was developed to simulate isotope values of PCE during the dissolution and degradation of PCE adjacent to a DNAPL source zone. With the exception of very high degradation rate constants (>1/day) stable carbon isotope values of PCE estimated by the model remained within error of the isotope value of the PCE DNAPL, consistent with measured isotope values in the model aquifer and in the PTC.     

Temporal evolution of DNAPL source and contaminant flux distribution: impacts of source mass depletion

We investigated, using model simulations, the changes occurring in the distribution of dense non-aqueous phase liquid (DNAPL) mass (Sn) within the source zone during depletion through dissolution, and the resulting changes in the contaminant flux distribution (J) at the source control plane (CP). Two numerical codes (ISCO3D and T2VOC) were used to simulate selected scenarios of DNAPL dissolution and transport in three-dimensional, heterogeneous, spatially correlated, random permeability fields with emplaced sources. Data from the model simulations were interpreted based on population statistics (mean, standard deviation, coefficient of variation) and spatial statistics (centroid, second moments, variograms). The mean and standard deviation of the Sn and J distributions decreased with source mass depletion by dissolution. The decrease in mean and standard deviation was proportional for the J distribution resulting in a constant coefficient of variation (CV), while for the Sn distribution, the mean decreased faster than the standard deviation. The spatial distributions exhibited similar behavior as the population distribution, i.e., the CP flux distribution was more stable (defined by temporally constant second moments and range of variograms) than the Sn distribution. These observations appeared to be independent of the heterogeneity of the permeability (k) field (variance of the log permeability field=1 and 2.45), correlation structure (positive vs. negative correlation between the k and Sn domains) and the DNAPL dissolution model (equilibrium vs. rate-limited), for the cases studied. Analysis of data from a flux monitoring field study (Hill Air Force Base, Utah) at a DNAPL source CP before and after source remediation also revealed temporal invariance of the contaminant flux distribution. These modeling and field observations suggest that the temporal evolution of the contaminant flux distribution can be estimated if the initial distribution is known. However, the findings are preliminary and broader implications to sampling strategies for remediation performance assessment need to be evaluated in additional modeling and experimental studies.     

Effects of source zone heterogeneity on surfactant-enhanced NAPL dissolution and resulting remediation end-points

The effectiveness of removal of nonaqueous phase liquids (NAPLs) from the entrapment source zone of the subsurface has been limited by soil heterogeneity and the inability to locate all entrapped sources. The goal of this study was to demonstrate the uncertainty of degree of source removal associated with aquifer heterogeneity. In this demonstration, source zone NAPL removal using surfactant-enhanced dissolution was considered. Model components that simulate the processes of natural dissolution in aqueous phase and surfactant-enhanced dissolution were incorporated into an existing code of contaminant transport. The dissolution modules of the simulator used previously developed Gilland-Sherwood type phenomenological models of NAPL dissolution to estimate mass transfer coefficients that are upscaleable to multidimensional flow conditions found at field sites. The model was used to simulate the mass removal from 10 NAPL entrapment zone configurations based on previously conducted two-dimensional tank experiments. These entrapment zones represent the NAPL distribution in spatially correlated random fields of aquifer hydraulic conductivity. The numerical simulations representing two-dimensional conditions show that effectiveness of mass removal depends on the aquifer heterogeneity that controls the NAPL entrapment and delivery of the surfactant to the locations of entrapped NAPLs. Flow bypassing resulting from heterogeneity and the reduction of relative permeability due to NAPL entrapment reduces the delivery efficiency of the surfactant, thus prolonging the remediation time to achieve desired end-point NAPL saturations and downstream dissolved concentrations. In some extreme cases, the injected surfactant completely bypassed the NAPL source zones. It was also found that mass depletion rates for different NAPL source configurations vary significantly. The study shows that heterogeneity result in uncertainties in the mass removal and achievable end-points that are directly related to dissolved contaminant plume development downstream of the NAPL entrapment zone.     

Effectiveness of nanoscale zero-valent iron for treatment of a PCE-DNAPL source zone

Nanoscale zero-valent iron (nZVI) has received considerable attention as a potential in situ remediation technology for treating chlorinated solvent source zones. Experimental and mathematical modeling studies were conducted to investigate the performance of nZVI in the transformation of tetrachloroethene (PCE) entrapped as a dense nonaqueous phase liquid (DNAPL). Injection of a 60 g/L suspension of nZVI into a column containing 20-30 mesh Ottawa sand and PCE-DNAPL at a residual saturation of 5.5% resulted in a uniform distribution of nZVI and minimal displacement of PCE. Subsequent flushing with 267 pore volumes of water containing 3mM CaCl(2) at a Darcy velocity of 0.75 m/day resulted in steady-state effluent concentrations of PCE near the solubility limit (ca. 200mg/L) and production of dissolved-phase ethene (10-30 mg/L). Over the duration of the experiment, approximately 30% of the initial PCE-DNAPL mass reacted to form ethene, 50% was eluted as dissolved-phase PCE, and 20% remained in the column as PCE-DNAPL. To further explore the implications of the nZVI column results, a multiphase transport model was developed that incorporated rate-limited PCE-DNAPL dissolution and reactions with nZVI. Using a fitted pseudo first-order transformation rate coefficient of 1.421/h, the model accurately captured observed trends in effluent concentrations of PCE and ethene and overall mass balance. A model sensitivity study reveals a strong dependence of treatment effectiveness on system characteristics. The sensitivity analysis suggests that an increase in the extent of PCE transformation is facilitated by decreasing flow rate, emplacement of nZVI down-gradient of the DNAPL source zone, and decreasing length of the DNAPL source zone. These findings indicate that, although emplacement of high concentrations of nZVI within a PCE-DNAPL source zone can result in substantial transformation of the parent compound, careful attention to design parameters (e.g. flow rate, location and amount nZVI delivered) will be required to achieve complete conversion to benign reaction products.     

Comparison of bioaugmentation and biostimulation for the enhancement of dense nonaqueous phase liquid source zone bioremediation.

Two 11.7-m(3) experimental controlled release systems (ECRS), packed with sandy model aquifer material and amended with tetrachloroethene (PCE) dense nonaqueous phase liquid (DNAPL) source zone, were operated in parallel with identical flow regimes and electron donor amendments. Hydrogen Releasing Compound (Regenesis Bioremediation Products, Inc., San Clemente, California), and later dissolved lactate, served as electron donors to promote dechlorination. One ECRS was bioaugmented with an anaerobic dechlorinating consortium directly into the source zone, and the other served as a control (biostimulated only) to determine the benefits of bioaugmentation. The presence of halorespiring bacteria in the aquifer matrix before bioaugmentation, shown by nested polymerase chain reaction with phylogenetic primers, suggests that dechlorinating catabolic potential may be somewhat widespread. Results obtained corroborate that source zone reductive dechlorination of PCE is possible at near field scale and that a system bioaugmented with a competent halorespiring consortium can enhance DNAPL dissolution and dechlorination processes at significantly greater rates than in a system that is biostimulated only.     

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

DNAPL source depletion: linking architecture and flux response

The relationship between dense non-aqueous phase liquid (DNAPL) mass reduction and contaminant mass flux was investigated experimentally in four model source zones. The flow cell design for the experiments featured a segmented extraction well that allowed for analysis of spatially resolved flux information. This flux information was coupled with image analysis of the NAPL spatial distribution to investigate the relationship between flux and the up-gradient NAPL architecture. Results indicate that in the systems studied, the relationship between DNAPL mass reduction and contaminant mass flux was primarily controlled by the NAPL architecture. A specific definition of NAPL architecture was employed where the source zone is resolved into a collection of streamtubes with spatial variability in NAPL saturation along each streamtube integrated and transformed into an effective NAPL content for each streamtube. The distribution of NAPL contents among the streamtubes (NAPL architecture) controlled dissolution dynamics. Two simplified models, a streamtube model and an effective Damkohler number model, were investigated for their ability to simulate dissolution dynamics.     

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.     

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.     

Multiphase flow and transport caused by spontaneous gas phase growth in the presence of dense non-aqueous phase liquid

Disconnected bubbles or ganglia of trapped gas may occur below the top of the capillary fringe through a number of mechanisms. In the presence of dense non-aqueous phase liquid (DNAPL), the disconnected gas phase experiences mass transfer of dissolved gases, including volatile components from the DNAPL. The properties of the gas phase interface can also change. This work shows for the first time that when seed gas bubbles exist spontaneous gas phase growth can be expected to occur and can significantly affect water-gas-DNAPL distributions, fluid flow, and mass transfer. Source zone behaviour was observed in three different experiments performed in a 2-dimensional flow cell. In each case, a DNAPL pool was created in a zone of larger glass beads over smaller glass beads, which served as a capillary barrier. In one experiment effluent water samples were analyzed to determine the vertical concentration profile of the plume above the pool. The experiments effectively demonstrated a) a cycle of spontaneous gas phase expansion and vertical advective mobilization of gas bubbles and ganglia above the DNAPL source zone, b) DNAPL redistribution caused by gas phase growth and mobilization, and c) that these processes can significantly affect mass transport from a NAPL source zone.     

Long-term ground penetrating radar monitoring of a small volume DNAPL release in a natural groundwater flow field

An earlier field experiment at Canadian Forces Base Borden by Brewster and Annan [Geophysics 59 (1994) 1211] clearly demonstrated the capability of ground penetrating radar (GPR) reflection profiling to detect and monitor the formation of DNAPL layers in the subsurface. Their experiment involved a large volume release (770 L) of tetrachloroethylene into a portion of the sand aquifer that was hydraulically isolated from groundwater flow by sheet pile walls. In this study, we evaluated the ability of GPR profiling to detect and monitor much smaller volume releases (50 L). No subsurface confining structure was used in this experiment; hence, the DNAPL impacted zone was subjected to the natural groundwater flow regime. This condition allowed us to geophysically monitor the DNAPL mass loss over a 66 month period. Reflectivity variations on the GPR profiles were used to infer the presence and evolution of the solvent layers. GPR imaging found significant reflectivity increases due to solvent layer formation during the two week period immediately after the release. These results demonstrated the capacity of GPR profiling for the detection and monitoring of lesser volume DNAPL releases that are more representative of small-scale industrial spills. The GPR imaged solvent layers subsequently reduced in both areal extent and reflectivity after 29 months and almost completely disappeared by the end of the 66 month monitoring period. Total DNAPL mass estimates based on GPR profiling data indicated that the solvent mass was reduced to 34%-36% of its maximum value after 29 months; only 4%-9% of the solvent mass remained in the study area after 66 months. These results are consistent with independent hydrogeological estimates of remaining DNAPL mass based on the downgradient monitoring of the dissolved solvent phase. Hence, we have concluded that the long-term GPR reflectivity changes of the DNAPL layers are likely the result from the dissolution of chlorinated solvents residing in those layers. The long-term monitoring results demonstrated that GPR profiling is a promising non-invasive method for use at DNAPL contaminated sites in sandy aquifers where temporal information about immiscible contaminant mass depletion due to either natural flow or remediation is needed. However, our results also indicated that the GPR signature of older DNAPL impacted zones may not differ greatly from the uncontaminated background if significant mass reduction due to dissolution has occurred.     

Assessing the impacts of partial mass depletion in DNAPL source zones I. Analytical modeling of source strength functions and plume response

Analytical solutions are developed for approximating the time-dependent contaminant discharge from DNAPL source zones undergoing dissolution and other decay processes. The source functions assume a power relationship between source mass and chemical discharge and can consider partial DNAPL source remediation (depletion) at any time after the initial DNAPL release. The source functions are used as a time-dependent boundary condition in an idealized chemical transport model to develop leading order approximations of the plume response to DNAPL source removal. The results suggest that partial DNAPL remediation does not tend to have a dramatic impact on the maximum extent of the plume if very low concentration values are used to define the plume boundaries. However, the solutions show that partial DNAPL removal from the source zone is likely to lead to large reductions in plume concentrations and mass, and it reduces the longevity of the plume. When the mass discharge from the source zone is linearly related to the DNAPL mass, it is shown that partial DNAPL depletion leads to linearly proportional reductions in the plume mass and concentrations.     

Mass removal of chlorinated ethenes from rough-walled fractures using permanganate

In situ chemical oxidation (ISCO) employing permanganate is an emerging technology that has been successful at enhancing mass removal from DNAPL source zones in unconsolidated media at the pilot-scale. The focus of this study was to evaluate the applicability of flushing a permanganate solution across two single vertical fractures in a laboratory environment to remove free phase DNAPL. The fracture experiments were designed to represent a portion of a larger fractured aquifer system impacted by a near-surface DNAPL spill over a shallow fractured rock aquifer. Each fracture was characterized by hydraulic and tracer tests, and the aperture field for one of the fractures was mapped using a co-ordinate measurement machine. Following DNAPL emplacement, a series of water and permanganate flushes were performed. To support observations from the fracture experiments, a set of batch experiments was conducted. The data from both fracture experiments showed that the post-oxidation effluent concentration was not impacted by the oxidant flush; however, changes in the aperture distribution, flow field, and flow rate were observed. These changes resulted in a significant decrease to the mass loading from the fractures, and were attributed to the build-up of oxidation by-products (manganese oxides and carbon dioxide) within the fracture which was corroborated by the batch experiment data and visual examination of the walls of one fracture. These results provide insight into the potential impact that a permanganate solution and oxidation by-products can have on the aperture distribution within a fracture and on DNAPL mass transfer rates. A permanganate flush or injection completed within a fractured rock aquifer may lead to the development of an insoluble product adjacent to the DNAPL which results in the reduction or complete elimination of advective regions near the DNAPL and reduces mass transfer rates. This outcome would have significant implications on the plume generating potential of the remaining DNAPL.     

Investigation of surfactant-enhanced dissolution of entrapped nonaqueous phase liquid chemicals in a two-dimensional groundwater flow field.

Because of their low solubility, waste chemicals in the form of nonaqueous phase liquids (NAPLs) that are entrapped in subsurface formations act as long-term sources of groundwater contamination. In the design of remediation schemes that use surfactants, it is necessary to estimate the mass transfer rate coefficients under multi-dimensional flow fields that exit at field sites. In this study, we investigate mass transfer under a two-dimensional flow field to obtain an understanding of the basic mechanisms of surfactant-enhanced dissolution and to quantify the mass transfer rates. Enhanced dissolution experiments in a two-dimensional test cell were conducted to measure rates of mass depletion from entrapped NAPLs to a flowing aqueous phase containing a surfactant. In situ measurement of transient saturation changes using a gamma attenuation system revealed dissolution patterns that are affected by the dimensionality of the groundwater flow field. Numerical modeling of local flow fields that changed with time, due to depletion of NAPL sources, enabled the examination of the basic mechanisms of NAPL dissolution in complex groundwater systems. Through nonlinear regression analysis, mass transfer rates were correlated to porous media properties, NAPL saturation and aqueous phase velocity. Results from the experiments and numerical analyses were used to identify deficiencies in existing methods of analysis that uses assumptions of one-dimensional flow, homogeneity of aquifer properties, local equilibrium and idealized transient mass transfer.     

A pragmatic approach for estimation of source-zone emissions at LNAPL contaminated sites

When considering natural attenuation as a remediation strategy at a site contaminated by a light non-aqueous phase liquid (LNAPL), it is important to consider the emission of contaminants from the source zone. A quantification of source-zone emissions is essential both for comparison with down-gradient mass fluxes to provide an estimate of fractional mass flux reduction, as well as for estimating the source lifetime. Because the spatial distribution of LNAPL at a field site is strongly dependent on both the spill circumstances and the heterogeneity of the geologic materials, which can be problematic for in-situ determination, alternative methods for estimating source-zone emissions are needed. In this work, a three-dimensional multiphase flow and transport modelling approach is used to investigate the relationship between the lateral extent of an LNAPL body and the emission of contaminants to groundwater at a contaminated site. For simulations involving an LNAPL release in an aquifer comprised of heterogeneous porosity and permeability distributions that were generated geostatistically, it is shown that a simple linear relationship exists between the lateral extent of the LNAPL body in the capillary fringe and the emission to the aqueous phase. The parameters describing the relationship are found to be linear functions of the groundwater flow velocity and the vertical infiltration rate. This site-specific relationship provides a simple method to estimate contaminant emissions to groundwater at LNAPL contaminated sites.