Wettability hysteresis and its implications for DNAPL source zone distribution

Subsurface heterogeneity at sites contaminated with nonaqueous phase liquids (NAPLs) reduces the effectiveness of traditional remediation measures. One cause may be the increased proportion of NAPL that is hydraulically isolated due to capillary trapping in heterogeneously-wetted materials. This study examines the wettability of ten materials, ranging from minerals, such as calcite and dolomite, to carbonaceous materials, such shale and coal, in air and water, NAPL and air, and NAPL and water systems. The wettability differed depending on which phase the solid material was initially immersed in: the less crystalline solids, if initially contacted by water were water-wet, but if initially contacted by NAPL were NAPL-wet. This difference, termed here wettability hysteresis, was observed for a suite of halogenated NAPLs and was independent of equilibration time. The degree of wettability hysteresis was greatest in the NAPL and water systems, with the magnitude of the difference increasing with the carbonaceous materials. Since the degree of capillary trapping in subsurface materials is related to wettability, the phenomenon of wettability hysteresis suggests that system history is a factor that may increase the heterogeneity of NAPL source zones.     

Laboratory investigation of flux reduction from dense non-aqueous phase liquid (DNAPL) partial source zone remediation by enhanced dissolution

This study investigated the benefits of partial removal of dense nonaqueous phase liquid (DNAPL) source zones using enhanced dissolution in eight laboratory scale experiments. The benefits were assessed by characterizing the relationship between reductions in DNAPL mass and the corresponding reduction in contaminant mass flux. Four flushing agents were evaluated in eight controlled laboratory experiments to examine the effects of displacement fluid property contrasts and associated override and underride on contaminant flux reduction (R(j)) vs. mass reduction (R(m)) relationships (R(j)(R(m))): 1) 50% ethanol/50% water (less dense than water), 2) 40% ethyl-lactate/60% water (more dense than water), 3) 18% ethanol/26% ethyl-lactate/56% water (neutrally buoyant), and 4) 2% Tween-80 surfactant (also neutrally buoyant). For each DNAPL architecture evaluated, replicate experiments were conducted where source zone dissolution was conducted with a single flushing event to remove most of the DNAPL from the system, and with multiple shorter-duration floods to determine the path of the R(j)(R(m)) relationship. All of the single-flushing experiments exhibited similar R(j)(R(m)) relationships indicating that override and underride effects associated with cosolvents did not significantly affect the remediation performance of the agents. The R(j)(R(m)) relationship of the multiple injection experiments for the cosolvents with a density contrast with water tended to be less desirable in the sense that there was less R(j) for a given R(m). UTCHEM simulations supported the observations from the laboratory experiments and demonstrated the capability of this model to predict R(j)(R(m)) relationships for non-uniformly distributed NAPL sources.     

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

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

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

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

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.     

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.     

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.     

Dissolved plume attenuation with DNAPL source remediation, aqueous decay and volatilization--analytical solution, model calibration and prediction uncertainty

A vertically-integrated analytical model for dissolved phase transport is described that considers a time-dependent DNAPL source based on the upscaled dissolution kinetics model of Parker and Park with extensions to consider time-dependent source zone biodecay, partial source mass reduction, and remediation-enhanced source dissolution kinetics. The model also considers spatial variability in aqueous plume decay, which is treated as the sum of aqueous biodecay and volatilization due to diffusive transport and barometric pumping through the unsaturated zone. The model is implemented in Excel/VBA coupled with (1) an inverse solution that utilizes prior information on model parameters and their uncertainty to condition the solution, and (2) an error analysis module that computes parameter covariances and total prediction uncertainty due to regression error and parameter uncertainty. A hypothetical case study is presented to evaluate the feasibility of calibrating the model from limited noisy field data. The results indicate that prediction uncertainty increases significantly over time following calibration, primarily due to propagation of parameter uncertainty. However, differences between the predicted performance of source zone partial mass reduction and the known true performance were reasonably small. Furthermore, a clear difference is observed between the predicted performance for the remedial action scenario versus that for a no-action scenario, which is consistent with the true system behavior. The results suggest that the model formulation can be effectively utilized to assess monitored natural attenuation and source remediation options if careful attention is given to model calibration and prediction uncertainty issues.     

表面活性剂强化抽出处理含水层中DNAPL污染物的去除特征

为明确表面活性剂强化抽出处理含水层中DNAPL污染物过程中表面活性剂的增强修复效果,及DNAPL自身理化性质和介质孔径的影响,利用数码图像分析技术对1,2-二氯乙烷和四氯乙烯2种DNAPL在石英砂填充的二维砂箱中的抽取迁移过程进行了实验模拟研究,并对抽出水样中DNAPL的浓度进行了测试分析。结果表明,实验条件下加入低浓度(0.18%)的十二烷基苯磺酸钠(SDBS)大幅提高了对弱透水层截留的2种DNAPL聚集体的抽出处理效率。1,2-二氯乙烷在该表面活性剂溶液中的表观溶解度远高于四氯乙烯,因此其短时间内的绝对去除率更高。SDBS强化抽出处理DNAPL的作用机理以增溶作用为主,而其增流作用使DNAPL迁移流动后分布面积增大,增加了与表面活性剂溶液接触的面积,对增溶作用起到促进效果。细粒介质中DNAPL迁移后的最大分布面积较大,因此体系中DNAPL的溶解速率较高。在DNAPL聚集体质量与水力梯度固定的条件下,油水界面张力越低,DNAPL的密度越大,DNAPL垂向迁移的风险就越大。本研究为修复工程中如何依据DNAPL种类与场地多孔介质的情况选择表面活性剂提供了参考。     

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.     

Controlled release, blind test of DNAPL remediation by ethanol flushing

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

Manipulation of density and viscosity for the optimization of DNAPL recovery by alcohol flooding

Laboratory experiments demonstrate that in situ recovery of pooled tetrachloroethene (PCE) from porous media may be accomplished more efficiently using multiple-step alcohol floods than with single alcohol floods. To optimize flooding efficiency while maintaining a low risk of downward DNAPL mobilization, a three-step flooding process is developed employing an isobutanol preflood, a composite alcohol mainflood, and a polymer solution postflood. The density and viscosity of these solutions are manipulated to prevent the onset and propagation of viscous and gravitational fingers, while maintaining phase behavior critical for efficient miscible NAPL displacement. An aqueous partitioning preflood solution of 10% by volume (10% v) isobutanol reduces the NAPL density in situ to approximately 1.00 g/ml by swelling the NAPL prior to miscible displacement induced by the mainflood. The composite alcohol mainflood, containing 65% v ethylene glycol and 35% v 1-propanol maintains miscibility while achieving neutral buoyancy and near stable displacement of the NAPL. Aqueous solutions of xanthan gum polymer efficiently displace the mainflood, reducing viscous fingering associated with waterfloods. Two-dimensional experiments using the multiple-step technique achieve 99.8% DNAPL mass recovery using a total of 0.45 pore volumes of alcohol, illustrating greater recovery efficiency than previous alcohol flooding formulations under comparable conditions.     

DNAPL distribution in the source zone: effect of soil structure and uncertainty reduction with increased sampling density

This paper focuses on parameters describing the distribution of dense nonaqueous phase liquid (DNAPL) contaminants and investigates the variability of these parameters that results from soil heterogeneity. In addition, it quantifies the uncertainty reduction that can be achieved with increased density of soil sampling. Numerical simulations of DNAPL releases were performed using stochastic realizations of hydraulic conductivity fields generated with the same geostatistical parameters and conditioning data at two sampling densities, thus generating two simulation ensembles of low and high density (three-fold increase) of soil sampling. The results showed that DNAPL plumes in aquifers identical in a statistical sense exhibit qualitatively different patterns, ranging from compact to finger-like. The corresponding quantitative differences were expressed by defining several alternative measures that describe the DNAPL plume and computing these measures for each simulation of the two ensembles. The uncertainty in the plume features under study was affected to different degrees by the variability of the soil, with coefficients of variation ranging from about 20% to 90%, for the low-density sampling. Meanwhile, the increased soil sampling frequency resulted in reductions of uncertainty varying from 7% to 69%, for low- and high-uncertainty variables, respectively. In view of the varying uncertainty in the characteristics of a DNAPL plume, remedial designs that require estimates of the less uncertain features of the plume may be preferred over others that need a more detailed characterization of the source zone architecture.     

The influence of cosolvent and heat on the solubility and reactivity of organophosphorous pesticide DNAPL alkaline hydrolysis

The presented research concerned the compatibility of cosolvents with in situ alkaline hydrolysis (ISAH) for treatment of organophosphorous (OPP) pesticide contaminated sites. In addition, the influence of moderate temperature heat increments was studied as a possible enhancement method. A complex dense non-aqueous phase liquid (DNAPL) of primarily parathion (~50 %) and methyl parathion (~15 %) obtained from the Danish Groyne 42 site was used as a contaminant source, and ethanol and propan-2-ol (0, 25, and 50 v/v%) was used as cosolvents in tap water and 0.34 M NaOH. Both cosolvents showed OPP solubility enhancement at 50 v/v% cosolvent content, with slightly higher OPP concentrations reached with propan-2-ol. Data on hydrolysis products did not show a clear trend with respect to alkaline hydrolysis reactivity in the presence of cosolvents. Results indicated that the hydrolysis rate of methyl-parathion (MP3) decreased with addition of cosolvent, whereas the hydrolysis rate of ethyl-parathion (EP3) remained constant, and overall indications were that the hydrolysis reactions were limited by the rate of hydrolysis rather than NAPL dissolution. In addition to cosolvents, the influence of low-temperature heating on ISAH was studied. Increasing reaction temperature from 10 to 30 °C provided an average rate of hydrolysis enhancement by a factor of 1.4–4.8 dependent on the base of calculation. When combining 50 v/v% cosolvent addition and heating to 30 °C, EP3 solubility was significantly enhanced and results for O,O-diethyl-thiophosphoric acid (EP2 acid) showed a significant enhancement of hydrolysis as well. However, this could not be supported by para-nitrophenol (PNP) data indicating the instability of this product in the presence of cosolvent.

     

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.     

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.     

Performance of a field-scale permeable reactive barrier based on organic substrate and zero-valent iron for in situ remediation of acid mine drainage

A permeable reactive barrier (PRB) was installed in Aznalcóllar (Spain) in order to rehabilitate the Agrio aquifer groundwater severely contaminated with acid mine drainage after a serious mining accident. The filling material of the PRB consisted of a mixture of calcite, vegetal compost and, locally, Fe⁰ and sewage sludge. Among the successes of the PRB are the continuous neutralisation of pH and the removal of metals from groundwater within the PRB (removals of >95 %). Among the shortcomings are the improper PRB design due to the complexity of the internal structure of the Agrio alluvial deposits (which resulted in an inefficient capture of the contaminated plume), the poor degradability of the compost used and the short residence time within the PRB (which hindered a complete sulphate reduction), the clogging of a section of the PRB and the heterogeneities of the filling material (which resulted in preferential flows within the PRB). Undoubtedly, it is only through accumulated experience at field-scale systems that the potentials and limits of the PRB technology can be determined.     

Modeling in situ ozonation for the remediation of nonvolatile PAH-contaminated unsaturated soils

Mathematical models were developed to investigate the characteristics of gaseous ozone transport under various soil conditions and the feasibility of in situ ozone venting for the remediation of unsaturated soils contaminated with phenanthrene. On the basis of assumptions for the mass transfer and reactions of ozone, three approaches were considered: equilibrium, kinetic, and lump models. Water-saturation-dependent reactions of gaseous ozone with soil organic matter (SOM) and phenanthrene were employed. The models were solved numerically by using the finite-difference method, and the model parameters were determined by using the experimental data of Hsu [The use of gaseous ozone to remediate the organic contaminants in the unsaturated soils, PhD Thesis, Michigan State Univ., East Lansing, MI, 1995]. The transport of gas-phase ozone is significantly retarded by ozone consumption due to reactions with SOM and phenanthrene, in addition to dissolution. An operation time of 156 h was required to completely remove phenanthrene in a 5-m natural soil column. In actual situations, however, the operation time is likely to be longer than the ideal time because of unknown factors including heterogeneity of the porous medium and the distribution of SOM and contaminant. The ozone transport front length was found to be very limited (< 1 m). The sensitivity analysis indicated that SOM is the single most important factor affecting in situ ozonation for the remediation of unsaturated soil contaminated with phenanthrene. Models were found to be insensitive to the reaction mechanisms of phenathrene with either gas-phase ozone or dissolved ozone. More study is required to quantify the effect of OH* formation on the removal of contaminant and on ozone transport in the subsurface.     

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.