Magnetic resonance imaging of nonaqueous phase liquid during soil vapor extraction in heterogeneous porous media

Soil vapor extraction (SVE) is commonly used to remediate nonaqueous phase liquids (NAPLs) from the vadose zone. This paper aims to determine the effect of grain size heterogeneity on the removal of NAPL in porous media during SVE. Magnetic resonance imaging (MRI) was used to observe and quantify the amount and location of NAPL in flow-through columns filled with silica gel grains. MRI is unique because it is nondestructive, allowing three-dimensional images to be taken of the phases as a function of space and time. Columns were packed with silica gel in three ways: coarse grains (250-550 microm) only, fine grains (32-63 microm) only, and a core of fine grains surrounded by a shell of coarse grains. Columns saturated with water were drained under a constant suction head, contaminated with decane, and then drained to different decane saturations. Each column was then continuously purged with water-saturated nitrogen gas and images were taken intermittently. Results showed that at residual saturation, a sharp volatilization front moved through the columns filled with either coarse-grain or fine-grain silica gel. In the heterogeneous columns, the volatilization front in the core lagged just behind the shell because gas flow was greater through the shell and decane in the core diffused outward to the shell. When decane saturation in the core was above residual saturation, decane volatilization occurred near the inlet, the relative decane saturation throughout the core dropped uniformly, and decane in the core flowed in the liquid phase to the shell to replenish volatilized decane. These results indicate that NAPL trapped in low-permeability zones can flow to replenish areas where NAPL is lost due to SVE. However, when residual NAPL saturation is reached, NAPL flow no longer occurs and diffusion limits removal from low-permeability zones.     

Effect of soil moisture dynamics on dense nonaqueous phase liquid (DNAPL) spill zone architecture in heterogeneous porous media

The amount, location, and form of NAPL in contaminated vadose zones are controlled by the spatial distribution of water saturation and soil permeability, the NAPL spill scenario, water infiltration events, and vapor transport. To evaluate the effects of these processes, we used the three-phase flow simulator STOMP, which includes a new permeability-liquid saturation-capillary pressure (k-S-P) constitutive model. This new constitutive model considers three NAPL forms: free, residual, and trapped. A 2-D vertical cross-section with five stratigraphic layers was assumed, and simulations were performed for seven cases. The conceptual model of the soil heterogeneity was based upon the stratigraphy at the Hanford carbon tetrachloride (CT) spill site. Some cases considered co-disposal of NAPL with large volumes of wastewater, as also occurred at the Hanford CT site. In these cases, the form and location of NAPL were most strongly influenced by high water discharge rates and NAPL evaporation to the atmosphere. In order to investigate the impact of heterogeneity, the hydraulic conductivity within the lower permeability layer was modeled as a realization of a random field having three different classes. For six extreme cases of 100 realizations, the CT mass that reached the water table varied by a factor of two, and was primarily controlled by the degree of lateral connectivity of the low conductivity class within the lowest permeability layer. The grid size at the top boundary had a dramatic impact on NAPL diffusive flux just after the spill event when the NAPL was present near the ground surface. NAPL evaporation with a fine grid spacing at the top boundary decreased CT mass that reached the water table by 74%, compared to the case with a coarse grid spacing, while barometric pumping had a marginal effect for the case of a continuous NAPL spill scenario considered in this work. For low water infiltration rate scenarios, the distribution of water content prior to a NAPL spill event decreased CT mass that reached the water table by 98% and had a significant impact on the formation of trapped NAPL. For all cases simulated, use of the new constitutive model that allows the formation of residual NAPL increased the amount of NAPL retained in the vadose zone. Density-driven advective gas flow from the ground surface controlled vapor migration in strongly anisotropic layers, causing NAPL mass flux to the lower layer to be reduced. These simulations indicate that consideration of the formation of residual and trapped NAPLs and dynamic boundary conditions (e.g., areas, rates, and periods of different NAPL and water discharge and fluctuations of atmospheric pressure) in the context of full three-phase flow are needed, especially for NAPL spill events at the ground surface. In addition, NAPL evaporation, density-driven gas advection, and NAPL vertical movement enhanced by water flow must be considered in order to predict NAPL distribution and migration in the vadose zone.     

Effect of water content on transient nonequilibrium NAPL-gas mass transfer during soil vapor extraction

The effect of water content on the volatilization of nonaqueous phase liquid (NAPL) in unsaturated soils was characterized by one-dimensional venting experiments conducted to evaluate the lumped mass transfer coefficient. An empirical correlation based upon the modified Sherwood number, Peclet number, and normalized mean grain size was used to estimate initial lumped mass transfer coefficients over a range of water content. The effects of water content on the soil vapor extraction SVE process have been investigated through experimentation and mathematical modeling. The experimental results indicated that a rate-limited NAPL-gas mass transfer occurred in water-wet soils. A severe mass transfer limitation was observed at 61.0% water saturation where the normalized effluent gas concentrations fell below 1.0 almost immediately, declined exponentially from the initiation of venting, and showed long tailing. This result was attributed to the reduction of interfacial area between the NAPL and mobile gas phases due to the increased water content. A transient mathematical model describing the change of the lumped mass transfer coefficient was used. Simulations showed that the nonequilibrium mass transfer process could be characterized by the exponent beta, a parameter which described the reduction of the specific area available for NAPL volatilization. The nonequilibrium mass transfer limitations were controlled by the soil mean grain size and pore gas velocity, were well described by beta values below 1.0 at low water saturation, and were well predicted with beta values greater than 1.0 at high water saturation.     

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.     

Investigation of surfactant-enhanced mass removal and flux reduction in 3D correlated permeability fields using magnetic resonance imaging

Magnetic resonance imaging (MRI) was used to visualize the NAPL source zone architecture before and after surfactant-enhanced NAPL dissolution in three-dimensional (3D) heterogeneously packed flowcells characterized by different longitudinal correlation lengths: 2.1 cm (aquifer 1) and 1.1 cm (aquifer 2). Surfactant flowpaths were determined by imaging the breakthrough of a paramagnetic tracer (MnCl(2)) analyzed by the method of moments. In both experimental aquifers, preferential flow occurred in high permeability materials with low NAPL saturations, and NAPL was preferentially removed from the top of the aquifers with low saturation. Alternate flushing with water and two surfactant pulses (5-6 pore volumes each) resulted in approximately 63% of NAPL mass removal from both aquifers. However, overall reduction in mass flux (Mass Flux 1) exiting the flowcell was lower in aquifer 2 (68%) than in aquifer 1 (81%), and local effluent concentrations were found to increase by as high as 120 times at local sampling ports from aquifer 2 after surfactant flushing. 3D MRI images of NAPL revealed that NAPL migrated downward and created additional NAPL source zones in previously uncontaminated areas at the bottom of the aquifers. The additional NAPL source zones were created in the direction transverse to flow in aquifer 2, which explains the higher mass flux relative to aquifer 1. Analysis using a total trapping number indicates that mobilization of NAPL trapped in the two coarsest sand fractions is possible when saturation is below 0.5 and 0.4, respectively. Results from this study highlight the potential impacts of porous media heterogeneity and NAPL source zone architecture on advanced in-situ flushing technologies.     

Reaction-enhanced mass transfer and transport from non-aqueous phase liquid source zones

The objective of the following research is to theoretically quantify the enhancement of interphase mass transfer of dissolved non-aqueous phase liquid (NAPL) compounds from the non-aqueous phase to the aqueous phase and the enhancement of dispersive mass transport from a NAPL zone due to destruction of dissolved NAPL compounds. For relatively slow reaction rates, such as for permanganate and perchloroethene (PCE), local-scale mass transfer enhancement is expected to be small. Dispersive mass transport with reaction from a horizontal NAPL zone can be quantified using equations derived for a mathematically equivalent falling film reactor system. In contrast to local-scale interphase mass transfer, dispersive mass transport from NAPL zones may be significantly increased by reaction. Enhancement factors due to destruction of the NAPL compound(s) are mainly dependent on NAPL solubility and oxidant concentration and to a lesser extent on reaction rate, stoichiometry, and transverse dispersion coefficients. Higher NAPL solubility and/or lower oxidant concentration reduces the maximum expected enhancement factor. Reaction enhancement factors for mass transport from NAPL zones are expected to be in the range of 5-50 for permanganate and chlorinated solvents. Theoretical results suggest that assuming instantaneous reaction rates may be appropriate for dispersive mass transport from NAPL zones.     

Petroleum mass removal from low permeability sediment using air sparging/soil vapor extraction: impact of continuous or pulsed operation

Air sparging and soil vapor extraction (AS/SVE) are innovative remediation techniques that utilize volatilization and microbial degradation to remediate petroleum spills from soils and groundwater. This in situ study investigated the use of AS/SVE to remediate a gasoline spill from a leaking underground storage tank (UST) in the low permeability, clayey soil of the Appalachian Piedmont. The objectives of this study were to evaluate AS/SVE in low permeability soils by quantifying petroleum mass removal rates, monitoring vadose zone contaminant levels, and comparing the mass extraction rates of continuous AS/SVE to 8 and 24 h pulsed operation. The objectives were met by collecting AS/SVE exhaust gas samples and vadose zone air from multi-depth soil vapor probes. Samples were analyzed for O₂, CO₂, BTEX (benzene, toluene, ethylbenzene, xylene), and total combustible hydrocarbon (TCH) concentrations using portable hand meters and gas chromatography. Continuous AS/SVE was effective in removing 608 kg of petroleum hydrocarbons from low permeability soil in 44 days (14.3 kg day⁻¹). Mass removal rates ranged from 2.6 times higher to 5.1 times lower than other AS/SVE studies performed in sandy sediments. BTEX levels in the vadose zone were reduced from about 5 ppm to 1 ppm. Ten pulsed AS/SVE tests removed 78 kg in 23 days and the mean mass removal rate (17.6 kg day⁻¹) was significantly higher than the last 15 days of continuous extraction. Pulsed operation may be preferable to continuous operation because of increased mass removal and decreased energy consumption.     

Non-aqueous phase liquid spreading during soil vapor extraction

Many non-aqueous phase liquids (NAPLs) are expected to spread at the air-water interface, particularly under non-equilibrium conditions. In the vadose zone, this spreading should increase the surface area for mass transfer and the efficiency of volatile NAPL recovery by soil vapor extraction (SVE). Observations of spreading on water wet surfaces led to a conceptual model of oil spreading vertically above a NAPL pool in the vadose zone. Analysis of this model predicts that spreading can enhance the SVE contaminant recovery compared to conditions where the liquid does not spread. Experiments were conducted with spreading volatile oils hexane and heptane in wet porous media and capillary tubes, where spreading was observed at the scale of centimeters. Within porous medium columns up to a meter in height containing stagnant gas, spreading was less than ten centimeters and did not contribute significantly to hexane volatilization. Water film thinning and oil film pinning may have prevented significant oil film spreading, and thus did not enhance SVE at the scale of a meter. The experiments performed indicate that volatile oil spreading at the field scale is unlikely to contribute significantly to the efficiency of SVE.     

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.     

Mass flux from a non-aqueous phase liquid pool considering spontaneous expansion of a discontinuous gas phase

The partitioning of non-aqueous phase liquid (NAPL) compounds to a discontinuous gas phase results in the repeated spontaneous expansion, snap-off, and vertical mobilization of the gas phase. This mechanism has the potential to significantly affect the mass transfer processes that control the dissolution of NAPL pools by increasing the vertical transport of NAPL mass and increasing the total mass transfer rate from the surface of the pool. The extent to which this mechanism affects mass transfer from a NAPL pool depends on the rate of expansion and the mass of NAPL compound in the gas phase. This study used well-controlled bench-scale experiments under no-flow conditions to quantify for the first time the expansion of a discontinuous gas phase in the presence of NAPL. Air bubbles placed in glass vials containing NAPL increased significantly in volume, from a radius of 1.0 mm to 2.0 mm over 215 days in the presence of tetrachloroethene (PCE), and from a radius of 1.2 mm to 2.3 mm over 22 days in the presence of trans-1,2-dichloroethene (tDCE). A one-dimensional mass transfer model, fit to the experimental data, showed that this expansion could result in a mass flux from the NAPL pool that was similar in magnitude to the mass flux expected for the dissolution of a NAPL pool in a two-fluid (NAPL and water) system. Conditions favouring the significant effect of a discontinuous gas phase on mass transfer were identified as groundwater velocities less than ~0.01 m/day, and a gas phase that covers greater than ~10% of the pool surface area and is located within ~0.01 m of the pool surface. Under these conditions the mass transfer via a discontinuous gas phase is expected to affect, for example, efforts to locate NAPL source zones using aqueous concentration data, and predict the lifetime and risk associated with NAPL source zones in a way that is not currently included in the common conceptual models used to assess NAPL-contaminated sites.     

Assessing performance and closure for soil vapor extraction: integrating vapor discharge and impact to groundwater quality

Soil vapor extraction (SVE) is typically effective for removal of volatile contaminants from higher-permeability portions of the vadose zone. However, contamination in lower-permeability zones can persist due to mass transfer processes that limit the removal effectiveness. After SVE has been operated for a period of time and the remaining contamination is primarily located in lower-permeability zones, the remedy performance needs to be evaluated to determine whether the SVE system should be optimized, terminated, or transitioned to another technology to replace or augment SVE. Numerical modeling of vapor-phase contaminant transport was used to investigate the correlation between measured vapor-phase mass discharge, MF(r), from a persistent, vadose-zone contaminant source and the resulting groundwater contaminant concentrations. This relationship was shown to be linear, and was used to directly assess SVE remediation progress over time and to determine the level of remediation in the vadose zone necessary to protect groundwater. Although site properties and source characteristics must be specified to establish a unique relation between MF(r) and the groundwater contaminant concentration, this correlation provides insight into SVE performance and support for decisions to optimize or terminate the SVE operation or to transition to another type of treatment.     

Hot water flushing for immiscible displacement of a viscous NAPL

Thermal remediation techniques, such as hot water flooding, are emerging technologies that have been proposed for the removal of nonaqueous phase liquids (NAPLs) from the subsurface. In this study a combined laboratory and modeling investigation was conducted to determine if hot water flooding techniques would improve NAPL mass removal compared to ambient temperature water flushing. Two experiments were conducted in a bench scale two-dimensional sandbox (55 cmx45 cmx1.3 cm) and NAPL saturations were quantified using a light transmission apparatus. In these immiscible displacement experiments the aqueous phase, at 22 degrees C and 50 degrees C, displaced a zone with initial NAPL saturations on the order of 85%. The interfacial tension and viscosity of the selected light NAPL, Voltesso 35, are strongly temperature-dependent. Experimental results suggest that hot water flooding reduced the size of the high NAPL saturation zone, in comparison to the cold water flood, and yielded greater NAPL mass recovery (75% NAPL removal vs. 64%). Hot water flooding did not, however, result in lower residual NAPL saturations. A numerical simulator was modified to include simultaneous flow of water and organic phases, energy transport, temperature and pressure. Model predictions of mass removal and NAPL saturation profiles compared well with observed behavior. A sensitivity analysis indicates that the utility of hot water flooding improves with the increasing temperature dependence of NAPL hydraulic properties.     

Axi-symmetric simulation of soil vapor extraction influenced by soil fracturing

Fracturing, either pneumatic or hydraulic, is a method to improve the performance of soil vapor extraction (SVE) in relatively low permeability soils (< 10(-5) cm/s). A two-dimensional model is presented to simulate trichloroethylene (TCE) soil vapor extraction modified by fracturing. Flow and transport is modeled using mobile macropore and micropore networks, which also have been identified in the literature as dual porosity, dual permeability, or heterogeneous flow models. In this model, fluids can flow in both the macropore and micropore networks. This represents a more general model compared to immobile micropore, mobile macropore models presented thus far in the literature for vapor flow and transport in two dimensions. The model considers pressure- and concentration-driven exchange between the macropore and micropore networks, concentration-driven exchange between the gas and sorbed phases within each network, and equilibrium exchange between the gas and water and a sorbed phase within each network. The parameters employed in an example simulation are based on field measurements made at a fractured site. Considered in the simulations were the influence of the volume percentage of fractures, the length of fractures, the relative location of the water table, and the influence of pulsed pumping. For these simulations, internetwork concentration-driven exchange most significantly affected mass removal. The volume percentage of fractures more significantly influence flow and mass removal than the length of fractures. The depth of the water table below the contamination plume only significantly influenced flow and mass removal when the water table was within 60 cm of the bottom of the contaminated soil in the vadose zone for the parameters considered in this study. Pulsed pumping was not found to increase the amount of mass removed in this study.     

Parameters that control the cleanup of fractured permeable aquifers

This study develops a modeling approach for simulating and evaluating entrapped light nonaqueous-phase liquid (light NAPL-LNAPL) dissolution and transport of the solute in a fractured permeable aquifer (FPA). The term FPA refers to an aquifer made of porous blocks of high permeability that embed fractures. The fracture network is part of the domain characterized by high permeability and negligible storage. Previous studies show that sandstone aquifers often represent FPAs. The basic model developed in this study is a two-dimensional (2-D) model of permeable blocks that embed oblique equidistant fractures with constant aperture and orientation. According to this model, two major parameters govern NAPL dissolution and transport of the solute. These parameters are: 1) the dimensionless interphase mass transfer coefficient, K(f0), and 2) the mobility number, N(M0). These parameters represent measures of heterogeneity affecting flow, NAPL dissolution, and transport of the solute in the domain. The parameter K(f0) refers to the rate at which organic mass is transferred from the NAPL into the water phase. The parameter N(M0) represents the ratio of flow through the porous blocks to flow through the fracture network in regions free of entrapped NAPL. It also provides a measure of groundwater flow bypassing regions contaminated by entrapped NAPL. In regions contaminated by entrapped NAPL our simulations have often indicated very low permeability of the porous blocks, enabling a significant increase of the fracture flow at the expense of the permeable block flow. Two types of constitutive relationships also affect the rate of FPA cleanup: 1) the relationship between the saturation of the entrapped NAPL and the permeability of the porous blocks, and 2) the relationships representing effects of the entrapped NAPL saturation and the permeable block flow velocity on rates of interphase mass transfer. This study provides basic tools for evaluating the characteristics of pump-and-treat cleanup of FPAs by referring to sets of parameters and constitutive relationships typical of FPAs. The numerical simulations carried out in this study show that at high initial saturation of the entrapped NAPL, during initial stages of the FPA cleanup the contaminant concentration increases, but later it decreases. This phenomenon originates from significant groundwater bypassing the NAPL entrapped in the permeable blocks via the fracture network.     

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.     

Prediction of mass-transfer coefficient for solute transport in porous media

Several previously reported laboratory studies related to transport of solutes through packed columns were utilized to develop predictive relationships for mass-transfer rate coefficient. The data were classified into two groups: those obtained under rate-limited mass transfer between mobile and immobile water regions (physical nonequilibrium conditions), and those derived from rate-limited mass transfer between instantaneous and slow sorption sites (sorption nonequilibrium conditions). The mass-transfer coefficient in all these studies was obtained by fitting breakthrough data to a transport model employing a first-order rate limitations with a "constant" mass-transfer coefficient, independent of flow conditions. This study demonstrated that the mass-transfer coefficient in these models is dependent on system parameters including pore-water velocity, length-scale, retardation coefficient, and particle or aggregate size. Predictive relationships were developed, through regression analysis, relating mass-transfer coefficient to residence time. The developed relationships adequately estimated previously reported field mass-transfer values. Successful simulations of field desorption data reported by Bahr [J. Contam. Hydrol. 4 (1989) 205] further demonstrate the potential applicability of the developed relationships.     

A high-resolution non-invasive approach to quantify oxygen transport across the capillary fringe and within the underlying groundwater

Oxygen transport across the capillary fringe is relevant for many biogeochemical processes. We present a non-invasive technique, based on optode technology, to measure high-resolution concentration profiles of oxygen across the unsaturated/saturated interface. By conducting a series of quasi two-dimensional flow-through laboratory experiments, we show that vertical hydrodynamic dispersion in the water-saturated part of the capillary fringe is the process limiting the mass transfer of oxygen. A number of experimental conditions were tested in order to investigate the influence of grain size and horizontal flow velocity on transverse vertical dispersion in the capillary fringe. In the same setup, analogous experiments were simultaneously carried out in the fully water-saturated zone, therefore allowing a direct comparison with oxygen transfer across the capillary fringe. The outcomes of the experiments under various conditions show that oxygen transport in the two zones of interest (i.e., the unsaturated/saturated interface and the saturated zone) is characterized by very similar transverse dispersion coefficients. An influence of the capillary fringe morphology on oxygen transport has not been observed. These results may be explained by the narrow grain size distribution used in the experiments, leading to a steep decline in water saturation at the unsaturated/saturated interface and to the absence of trapped gas in this transition zone. We also modeled flow (applying the van Genuchten and the Brooks-Corey relationships) and two-dimensional transport across the capillary fringe, obtaining simulated profiles of equivalent aqueous oxygen concentration that were in good agreement with the observations.     

A large-scale experiment on mass transfer of trichloroethylene from the unsaturated zone of a sandy aquifer to its interfaces

A large-scale experiment was conducted to investigate the transport of trichloroethylene (TCE) vapors in the unsaturated zone and to determine the mass transfer to the groundwater and the atmosphere. The experiment involved injection of 5 1 of TCE in the unsaturated zone under controlled conditions, with multidepth sampling of gas and water through the unsaturated zone and across the capillary zone into underlying groundwater. The mass transfer of TCE vapors from the vadose zone to the atmosphere was quantified using a vertical flux chamber. A special soil water sampler was used to monitor transport across the capillary fringe. Experimental data indicated that TCE in the unsaturated zone was mainly transported to the atmosphere and this exchange reduced significantly the potential for groundwater pollution. The maximum measured TCE flux to the atmosphere was about 3 g/m(2)/day. Observed and calculated fluxes based on vertical TCE vapor concentration gradients and Fick's law were in good agreement. This confirms that TCE vapor transport under the experimental conditions was governed essentially by molecular diffusion. TCE vapors also caused a lower, but significant contamination of the underlying groundwater by dispersion across the capillary fringe with a corresponding maximum flux of about 0.1 g/m(2)/day. This mass transfer to groundwater is partly uncertain due to an inadvertent entry of some nonaqueous phase liquid (NAPL) from the source area into the saturated zone. Application of an analytical solution to estimate the TCE flux from the unsaturated zone to the groundwater indicated that this phenomenon is not only influenced by molecular diffusion but also by vertical dispersion. The mass balance indicates that, under the given experimental conditions (e.g. proximity of the source emplacement relative to the soil surface, relatively high permeable porous medium), nearly 95% of the initial TCE mass was transferred to the atmosphere.     

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.     

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.