Natural attenuation of chlorinated ethene compounds: model development and field-scale application at the Dover site

A multi-dimensional and multi-species reactive transport model was developed to aid in the analysis of natural attenuation design at chlorinated solvent sites. The model can simulate several simultaneously occurring attenuation processes including aerobic and anaerobic biological degradation processes. The developed model was applied to analyze field-scale transport and biodegradation processes occurring at the Area-6 site in Dover Air Force Base, Delaware. The model was calibrated to field data collected at this site. The calibrated model reproduced the general groundwater flow patterns, and also, it successfully recreated the observed distribution of tetrachloroethene (PCE), trichloroethene (TCE), dichloroethylene (DCE), vinyl chloride (VC) and chloride plumes. Field-scale decay rates of these contaminant plumes were also estimated. The decay rates are within the range of values that were previously estimated based on lab-scale microcosm and field-scale transect analyses. Model simulation results indicated that the anaerobic degradation rate of TCE, source loading rate, and groundwater transport rate are the important model parameters. Sensitivity analysis of the model indicated that the shape and extent of the predicted TCE plume is most sensitive to transmissivity values. The total mass of the predicted TCE plume is most sensitive to TCE anaerobic degradation rates. The numerical model developed in this study is a useful engineering tool for integrating field-scale natural attenuation data within a rational modeling framework. The model results can be used for quantifying the relative importance of various simultaneously occurring natural attenuation processes.     

Reactive transport modeling of processes controlling the distribution and natural attenuation of phenolic compounds in a deep sandstone aquifer

Reactive solute transport modeling was utilized to evaluate the potential for natural attenuation of a contaminant plume containing phenolic compounds at a chemical producer in the West Midlands, UK. The reactive transport simulations consider microbially mediated biodegradation of the phenolic compounds (phenols, cresols, and xylenols) by multiple electron acceptors. Inorganic reactions including hydrolysis, aqueous complexation, dissolution of primary minerals, formation of secondary mineral phases, and ion exchange are considered. One-dimensional (1D) and three-dimensional (3D) simulations were conducted. Mass balance calculations indicate that biodegradation in the saturated zone has degraded approximately 1-5% of the organic contaminant plume over a time period of 47 years. Simulations indicate that denitrification is the most significant degradation process, accounting for approximately 50% of the organic contaminant removal, followed by sulfate reduction and fermentation reactions, each contributing 15-20%. Aerobic respiration accounts for less than 10% of the observed contaminant removal in the saturated zone. Although concentrations of Fe(III) and Mn(IV) mineral phases are high in the aquifer sediment, reductive dissolution is limited, producing only 5% of the observed mass loss. Mass balance calculations suggest that no more than 20-25% of the observed total inorganic carbon (TIC) was generated from biodegradation reactions in the saturated zone. Simulations indicate that aerobic biodegradation in the unsaturated zone, before the contaminant entered the aquifer, may have produced the majority of the TIC observed in the plume. Because long-term degradation is limited to processes within the saturated zone, use of observed TIC concentrations to predict the future natural attenuation may overestimate contaminant degradation by a factor of 4-5.     

The variability and intrinsic remediation of a BTEX plume in anaerobic sulphate-rich groundwater

Data from long-term groundwater sampling, limited coring, and associated studies are synthesised to assess the variability and intrinsic remediation/natural attenuation of a dissolved hydrocarbon plume in sulphate-rich anaerobic groundwater. Fine vertical scale (0.25- and 0.5-m depth intervals) and horizontal plume-scale (>400 m) characteristics of the plume were mapped over a 5-year period from 1991 to 1996. The plume of dissolved BTEX (benzene, toluene, ethylbenzene, xylene) and other organic compounds originated from leakage of gasoline from a subsurface fuel storage tank. The plume was up to 420 m long, less than 50 m wide and 3 m thick. In the first few years of monitoring, BTEX concentrations near the point of leakage were in approximate equilibrium with non-aqueous phase liquid (NAPL) gasoline. NAPL composition of core material and long-term trends in ratios of BTEX concentrations in groundwater indicated significant depletion (water washing, volatilisation and possibly biodegradation) of benzene from residual NAPL after 1992. Large fluctuations in BTEX concentrations in individual boreholes were shown to be largely attributable to seasonal groundwater flow variations. A combination of temporal and spatial groundwater quality data was required to adequately assess the stationarity of plumes, so as to allow inference of intrinsic remediation. Contoured concentration data for the period 1991 to 1996 indicated that plumes of toluene and o-xylene were, at best, only partially steady state (pseudo-steady state) due to seasonal groundwater flow changes. From this analysis, it was inferred that significant remediation by natural biodegradation was occurring for BTEX component plumes such as toluene and o-xylene, but provided no conclusive evidence of benzene biodegradation. Issues associated with field quantification of intrinsic remediation from groundwater sampling are highlighted. Preferential intrinsic biodegradation of selected organic compounds within the BTEX plume was shown to be occurring, in parallel with sulphate reduction and bicarbonate production. Ratios of average hydrocarbon concentrations to benzene for the period 1991 to 1992 were used to estimate degradation rates (half-lives) at various distances along the plume. The estimates varied with distance, the narrowest range being, for toluene, 110 to 260 days. These estimates were comparable to rates determined previously from an in situ tracer test and from plume-scale modelling.     

Oxygenated gasoline release in the unsaturated zone, Part 2: Downgradient transport of ethanol and hydrocarbons

In the event of a gasoline spill containing oxygenated compounds such as ethanol and MTBE, it is important to consider the impacts these compounds might have on subsurface contamination. One of the main concerns commonly associated with ethanol is that it might decrease the biodegradation of aromatic hydrocarbon compounds, leading to an increase in the hydrocarbon dissolved plume lengths. The first part of this study (Part 1) showed that when gasoline containing ethanol infiltrates the unsaturated zone, ethanol is likely to partition to and be retained in the unsaturated zone pore water. In this study (Part 2), a controlled field test is combined with a two-dimensional laboratory test and three-dimensional numerical modelling to investigate how ethanol retention in the unsaturated zone affects the downgradient behaviour of ethanol and aromatic hydrocarbon compounds. Ethanol transport downgradient was extremely limited. The appearance of ethanol in downgradient wells was delayed and the concentrations were lower than would be expected based on equilibrium dissolution. Oscillations in the water table resulted in minor flushing of ethanol, but its effect could still be perceived as an increase in the groundwater concentrations downgradient from the source zone. Ethanol partitioning to the unsaturated zone pore water reduced its mass fraction within the NAPL thus reducing its anticipated impact on the fate of the hydrocarbon compounds. A conceptual numerical simulation indicated that the potential ethanol-induced increase in benzene plume length after 20 years could decrease from 136% to 40% when ethanol retention in the unsaturated zone is considered.     

Combined application of conservative transport modelling and compound-specific carbon isotope analyses to assess in situ attenuation of benzene, toluene, and o-xylene

In recent years, compound specific isotope analyses (CSIA) have developed into one of the most powerful tools for the quantification of in situ biodegradation of organic contaminants. In this approach, the calculation of the extent of biodegradation of organic contaminants in aquifers is usually based on the Rayleigh equation, and thus neglects physical transport processes such as dispersion that contribute to contaminant dilution in aquifers. Here we combine compound specific isotope analyses with a conservative transport model to study the attenuation of aromatic hydrocarbons at a former gasworks site. The conservative transport model was first used to simulate concentration reductions caused by dilution at wells downgradient of a BTEX source. In a second step, the diluted concentrations, together with the available stable carbon isotope ratios and carbon fractionation factors for benzene, toluene and o-xylene were applied in the Rayleigh equation to quantify the degree of biodegradation at each of those wells. At the investigated site, where other attenuation processes such as sorption and volatilisation were proven to be negligible, the combined approach is recommended for benzene, which represents a compound for which the effect of biodegradation is comparable to or less than the effect of dilution. As demonstrated for toluene and o-xylene, the application of the Rayleigh equation alone is sufficient if dilution can be proved to be insignificant in comparison to biodegradation. The analysis also suggests that the source width and the position of the observation wells relative to the plume center line are significantly related to the degree of dilution.     

Persistence of LNAPL sources: relationship between risk reduction and LNAPL recovery

Light nonaqueous phase liquids (LNAPLs), such as fuels, are the source of much soil and groundwater contamination. Though the mobility of LNAPLs is limited in many environments, dissolved-phase components, such as benzene, can produce groundwater plumes that are more mobile than the LNAPL source. In such a setting, it is commonly assumed that recovery of the LNAPL will result in a reduction in risk associated with the dissolved phase. This paper synthesizes several existing multiphase and chemical transport solutions into a single linked methodology that predicts concentrations of soluble constituents within and downgradient of LNAPL source zones from dissolution of those constituents into groundwater flowing through and below LNAPL sources. This approach has been applied to a variety of LNAPL spill conditions. For biodegradable compounds, these analyses show that the period of time where the dissolved-phase plume is expanding is very small compared to the duration of most LNAPL sources, and that the downgradient extent is generally less than about 100 m for BTEX compounds. Therefore, the risk to receptors, as measured by the maximum downgradient extent of dissolved-phase plume or the maximum concentration of these compounds at a downgradient receptor, is generally unrelated to the longevity of the LNAPL sources. The maximum downgradient extent of the dissolved-phase plume is determined almost entirely by the groundwater velocity and the biodegradation rate. These analyses further demonstrate that recovery of LNAPL by hydraulic methods is often ineffective at reducing risk. Except in coarse-grained soils or intermediate soils with significant LNAPL saturations, free-product recovery approaches do not result in significant reductions in the longevity of downgradient dissolved-phase contamination. Further, for biodegradable constituents, remediation does not result in a near-term decrease in the downgradient extent of contamination. Cleanup methods that act to change the composition of the LNAPL source are more effective at reducing the downgradient concentrations, particularly for fine-grained soils or when LNAPL saturations are low.     

Assessing the natural attenuation of organic contaminants in aquifers using plume-scale electron and carbon balances: model development with analysis of uncertainty and parameter sensitivity

A quantitative methodology is described for the field-scale performance assessment of natural attenuation using plume-scale electron and carbon balances. This provides a practical framework for the calculation of global mass balances for contaminant plumes, using mass inputs from the plume source, background groundwater and plume residuals in a simplified box model. Biodegradation processes and reactions included in the analysis are identified from electron acceptors, electron donors and degradation products present in these inputs. Parameter values used in the model are obtained from data acquired during typical site investigation and groundwater monitoring studies for natural attenuation schemes. The approach is evaluated for a UK Permo-Triassic Sandstone aquifer contaminated with a plume of phenolic compounds. Uncertainty in the model predictions and sensitivity to parameter values was assessed by probabilistic modelling using Monte Carlo methods. Sensitivity analyses were compared for different input parameter probability distributions and a base case using fixed parameter values, using an identical conceptual model and data set. Results show that consumption of oxidants by biodegradation is approximately balanced by the production of CH4 and total dissolved inorganic carbon (TDIC) which is conserved in the plume. Under this condition, either the plume electron or carbon balance can be used to determine contaminant mass loss, which is equivalent to only 4% of the estimated source term. This corresponds to a first order, plume-averaged, half-life of > 800 years. The electron balance is particularly sensitive to uncertainty in the source term and dispersive inputs. Reliable historical information on contaminant spillages and detailed site investigation are necessary to accurately characterise the source term. The dispersive influx is sensitive to variability in the plume mixing zone width. Consumption of aqueous oxidants greatly exceeds that of mineral oxidants in the plume, but electron acceptor supply is insufficient to meet the electron donor demand and the plume will grow. The aquifer potential for degradation of these contaminants is limited by high contaminant concentrations and the supply of bioavailable electron acceptors. Natural attenuation will increase only after increased transport and dilution.     

Modelling of iron cycling and its impact on the electron balance at a petroleum hydrocarbon contaminated site in Hnevice, Czech Republic

Over a period of several decades multiple leaks of large volumes from storage facilities located near Hnevice (Czech Republic) have caused the underlying Quaternary aquifer to be severely contaminated with nonaqueous phase liquid (NAPL) petroleum hydrocarbons. Beginning in the late 1980's the NAPL plume started to shrink as a consequence of NAPL dissolution exceeding replenishment and due to active remediation. The subsurface was classified geochemically into four different zones, (i) a contaminant-free zone never occupied by NAPL or dissolved contaminants, (ii) a re-oxidation zone formerly occupied by NAPL, (iii) a zone currently occupied by NAPL, and (iv) a lower fringe zone between the overlying NAPL and the deeper underlying contaminant-free zone. The study investigated the spatial and temporal variability of the redox zonation at the Hnevice site and quantified the influence of iron-cycling on the overall electron balance. As a first step inverse geochemical modelling was carried out to identify possible reaction models and mass transfer processes. In a subsequent step, two-dimensional (forward) multi-component reactive transport modelling was performed to evaluate and quantify the major processes that control the geochemical evolution at the site. The study explains the observed enrichment of the lower fringe zone with ferrihydrite as a result of the re-oxidation of ferrous iron. It suggests that once the NAPL zone started to shrink the dissolution of previously formed siderite and FeS by oxygen and nitrate consumed a significant part of the oxidation capacity for a considerable time period and therefore limited the penetration of electron acceptors into the NAPL contaminated zone.     

Biodegradation of petroleum hydrocarbon vapors: laboratory studies on rates and kinetics in unsaturated alluvial sand

Predictions of natural attenuation of volatile organic compounds (VOCs) in the unsaturated zone rely critically on information about microbial biodegradation kinetics. This study aims at determining kinetic rate laws for the aerobic biodegradation of a mixture of 12 volatile petroleum hydrocarbons and methyl tert-butyl ether (MTBE) in unsaturated alluvial sand. Laboratory column and batch experiments were performed at room temperature under aerobic conditions, and a reactive transport model for VOC vapors in soil gas coupled to Monod-type degradation kinetics was used for data interpretation. In the column experiment, an acclimatization of 23 days took place before steady-state diffusive vapor transport through the horizontal column was achieved. Monod kinetic parameters Ks and vmax could be derived from the concentration profiles of toluene, m-xylene, n-octane, and n-hexane, because substrate saturation was approached with these compounds under the experimental conditions. The removal of cyclic alkanes, isooctane, and 1,2,4-trimethylbenzene followed first-order kinetics over the whole concentration range applied. MTBE, n-pentane, and chlorofluorocarbons (CFCs) were not visibly degraded. Batch experiments suggested first-order disappearance rate laws for all VOCs except n-octane, which decreased following zero-order kinetics in live batch experiments. For many compounds including MTBE, disappearance rates in abiotic batch experiments were as high as in live batches indicating sorption. It was concluded that the column approach is preferable for determining biodegradation rate parameters to be used in risk assessment models.     

Biodegradation of hydrocarbons vapors: Comparison of laboratory studies and field investigations in the vadose zone at the emplaced fuel source experiment, Airbase Vaerløse, Denmark

The natural attenuation of volatile organic compounds (VOCs) in the unsaturated zone can only be predicted when information about microbial biodegradation rates and kinetics are known. This study aimed at determining first-order rate coefficients for the aerobic biodegradation of 13 volatile petroleum hydrocarbons which were artificially emplaced as a liquid mixture during a field experiment in an unsaturated sandy soil. Apparent first-order biodegradation rate coefficients were estimated by comparing the spatial evolution of the resulting vapor plumes to an analytical reactive transport model. Two independent reactive numerical model approaches have been used to simulate the diffusive migration of VOC vapors and to estimate degradation rate coefficients. Supplementary laboratory column and microcosm experiments were performed with the sandy soil at room temperature under aerobic conditions. First-order kinetics adequately matched the lab column profiles for most of the compounds. Consistent compound-specific apparent first-order rate coefficients were obtained by the three models and the lab column experiment, except for benzene. Laboratory microcosm experiments lacked of sensitivity for slowly degrading compounds and underestimated degradation rates by up to a factor of 5. Addition of NH3 vapor was shown to increase the degradation rates for some VOCs in the laboratory microcosms. All field models suggested a significantly higher degradation rate for benzene than the rates measured in the lab, suggesting that the field microbial community was superior in developing benzene degrading activity.     

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

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

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.     

Analytical modelling of stable isotope fractionation of volatile organic compounds in the unsaturated zone

Analytical models were developed that simulate stable isotope ratios of volatile organic compounds (VOCs) near a point source contamination in the unsaturated zone. The models describe diffusive transport of VOCs, biodegradation and source ageing. The mass transport is governed by Fick's law for diffusion. The equation for reactive transport of VOCs in the soil gas phase was solved for different source geometries and for different boundary conditions. Model results were compared to experimental data from a one-dimensional laboratory column and a radial-symmetric field experiment. The comparison yielded a satisfying agreement. The model results clearly illustrate the significant isotope fractionation by gas phase diffusion under transient state conditions. This leads to an initial depletion of heavy isotopes with increasing distance from the source. The isotope evolution of the source is governed by the combined effects of isotope fractionation due to vaporisation, diffusion and biodegradation. The net effect can lead to an enrichment or depletion of the heavy isotope in the remaining organic phase, depending on the compound and element considered. Finally, the isotope evolution of molecules migrating away from the source and undergoing degradation is governed by a combined degradation and diffusion isotope effect. This suggests that, in the unsaturated zone, the interpretation of biodegradation of VOC based on isotopic data must always be based on a model combining gas phase diffusion and degradation.     

Tracer test for the measurement of gas diffusion and non-aqueous phase liquid (NAPL) saturation in soil

During soil bioremediation, the diffusion of oxygen into the soil is an important prerequisite for aerobic biodegradation, and the decrease of petroleum products is the ultimate goal. Both processes need to be monitored. The aim of this work was to develop a gas tracer test that yields information on both, gas diffusion and residual saturation with non-aqueous phase liquids (NAPLs) in unsaturated soil heaps. One conservative tracer (methane) and 4 partitioning gas tracers (diethylether, methyl tert-butyl ether, chloroform and n-heptane) were injected as vapors into laboratory columns filled with unsaturated sand with increasing NAPL saturation. Breakthrough curves of gaseous compounds were measured at two points and compared to analytical solutions of an analytical diffusive-reactive transport equation. By fitting of methane data, robust results for effective diffusivity (tortuosity) were obtained. NAPL saturation was most accurately measured by the moderately water soluble tracers (ethers and chloroform). The hydrophobic tracer n-heptane did not partition into water-immersed NAPL. An easy and accurate way to assess air-NAPL partitioning constants from gas chromatography retention times is furthermore reported. It is concluded that gas tracer tests have the potential for measuring two important properties in soil bioremediation systems easily and quickly.     

Inverse modeling of BTEX dissolution and biodegradation at the Bemidji, MN crude-oil spill site

The U.S. Geological Survey (USGS) solute transport and biodegradation code BIOMOC was used in conjunction with the USGS universal inverse modeling code UCODE to quantify field-scale hydrocarbon dissolution and biodegradation at the USGS Toxic Substances Hydrology Program crude-oil spill research site located near Bemidji, MN. This inverse modeling effort used the extensive historical data compiled at the Bemidji site from 1986 to 1997 and incorporated a multicomponent transport and biodegradation model. Inverse modeling was successful when coupled transport and degradation processes were incorporated into the model and a single dissolution rate coefficient was used for all BTEX components. Assuming a stationary oil body, we simulated benzene, toluene, ethylbenzene, m,p-xylene, and o-xylene (BTEX) concentrations in the oil and ground water, respectively, as well as dissolved oxygen. Dissolution from the oil phase and aerobic and anaerobic degradation processes were represented. The parameters estimated were the recharge rate, hydraulic conductivity, dissolution rate coefficient, individual first-order BTEX anaerobic degradation rates, and transverse dispersivity. Results were similar for simulations obtained using several alternative conceptual models of the hydrologic system and biodegradation processes. The dissolved BTEX concentration data were not sufficient to discriminate between these conceptual models. The calibrated simulations reproduced the general large-scale evolution of the plume, but did not reproduce the observed small-scale spatial and temporal variability in concentrations. The estimated anaerobic biodegradation rates for toluene and o-xylene were greater than the dissolution rate coefficient. However, the estimated anaerobic biodegradation rates for benzene, ethylbenzene, and m,p-xylene were less than the dissolution rate coefficient. The calibrated model was used to determine the BTEX mass balance in the oil body and groundwater plume. Dissolution from the oil body was greatest for compounds with large effective solubilities (benzene) and with large degradation rates (toluene and o-xylene). Anaerobic degradation removed 77% of the BTEX that dissolved into the water phase and aerobic degradation removed 17%. Although goodness-of-fit measures for the alternative conceptual models were not significantly different, predictions made with the models were quite variable.     

A set of constitutive relationships accounting for residual NAPL in the unsaturated zone

Although laboratory experiments show that non-aqueous phase liquid (NAPL) is retained in the unsaturated zone, no existing multiphase flow model has been developed to account for residual NAPL after NAPL drainage in the unsaturated zone. We developed a static constitutive set of saturation-capillary pressure relationships for water, NAPL and air that accounts for both this residual NAPL and entrapped NAPL. The set of constitutive relationships is formulated similarly to the set of scaled relationships that is frequently applied in continuum models. The new set consists of three fluid-phase systems: a three-phase system and a two-phase system, that both comply with the original constitutive model, and a newly introduced residual NAPL system. The new system can be added relatively easily to the original two- and three-phase systems. Entrapment is included in the model. The constitutive relationships of the non-drainable residual NAPL system are based on qualitative fluid behavior derived from a pore scale model. The pore scale model reveals that the amount of residual NAPL depends on the spreading coefficient and the water saturation. Furthermore, residual NAPL is history-dependent. At the continuum scale, a critical NAPL pressure head defines the transition from free, mobile NAPL to residual NAPL. Although the Pc-S relationships for water and total liquid are not independent in case of residual NAPL, two two-phase Pc-S relations can represent a three-phase residual system of Pc-S relations. A newly introduced parameter, referred to as the residual oil pressure head, reflects the mutual dependency of water and oil. Example calculations show consistent behavior of the constitutive model. Entrapment and retention in the unsaturated zone cooperate to retain NAPL. Moreover, the results of our constitutive model are in agreement with experimental observations.     

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.     

Groundwater-surface water interaction and its role on TCE groundwater plume attenuation

A field investigation of a TCE plume in a surficial sand aquifer shows that groundwater-surface water interactions strongly influence apparent plume attenuation. At the site, a former industrial facility in Connecticut, depth-discrete monitoring along three cross-sections (transects) perpendicular to groundwater flow shows a persistent VOC plume extending 700 m from the DNAPL source zone to a mid-size river. Maximum TCE concentrations along a transect 280 m from the source were in the 1000s of microg/L with minimal degradation products. Beyond this, the land surface drops abruptly to a lower terrace where a shallow pond and small streams occur. Two transects along the lower terrace, one midway between the facility and river just downgradient of the pond and one along the edge of the river, give the appearance that the plume has strongly attenuated. At the river, maximum TCE concentrations in the 10s of microg/L and similar levels of its degradation product cis-DCE show direct plume discharge from groundwater to the river is negligible. Although degradation plays a role in the strong plume attenuation, the major attenuation factor is partial groundwater plume discharge to surface water (i.e. the pond and small streams), where some mass loss occurs via water-air exchange. Groundwater and stream mass discharge estimates show that more than half of the plume mass discharge crossing the first transect, before surface water interactions occur, reaches the river directly via streamflow, although river concentrations were below detection due to dilution. This study shows that groundwater and surface water concentration measurements together provide greater confidence in identifying and quantifying natural attenuation processes at this site, rather than groundwater measurements alone.     

A functional relation for field-scale nonaqueous phase liquid dissolution developed using a pore network model

A pore network model with cubic chambers and rectangular tubes was used to estimate the nonaqueous phase liquid (NAPL) dissolution rate coefficient, Kdissai, and NAPL/water total specific interfacial area, ai. Kdissai was computed as a function of modified Peclet number (Pe') for various NAPL saturations (SN) and ai during drainage and imbibition and during dissolution without displacement. The largest contributor to ai was the interfacial area in the water-filled corners of chambers and tubes containing NAPL. When Kdissai was divided by ai, the resulting curves of dissolution coefficient, Kdiss versus Pe' suggested that an approximate value of Kdiss could be obtained as a weak function of hysteresis or SN. Spatially and temporally variable maps of Kdissai calculated using the network model were used in field-scale simulations of NAPL dissolution. These simulations were compared to simulations using a constant value of Kdissai and the empirical correlation of Powers et al. [Water Resour. Res. 30(2) (1994b) 321]. Overall, a methodology was developed for incorporating pore-scale processes into field-scale prediction of NAPL dissolution.     

Design and laboratory testing of a chamber device to measure total flux of volatile organic compounds from the unsaturated zone under natural conditions

To determine if an aquifer contaminated with volatile organic compounds (VOCs) has potential for natural remediation, all natural processes affecting the fate and transport of VOCs in the subsurface must be identified and quantified. This research addresses the quantification of air-phase volatile organic compounds (VOCs) leaving the unsaturated zone soil gas and entering the atmosphere-including the additional flux provided by advective soil-gas movement induced by barometric pumping. A simple and easy-to-use device for measuring VOC flux under natural conditions is presented. The vertical flux chamber (VFC) was designed using numerical simulations and evaluated in the laboratory. Mass-balance numerical simulations based on continuously stirred tank reactor equations (CSTR) provided information on flux measurement performance of several sampling configurations with the final chamber configuration measuring greater than 96% of model-simulated fluxes. A laboratory device was constructed to evaluate the flux chamber under both diffusion-only and advection-plus-diffusion transport conditions. The flux chamber measured an average of 82% of 15 diffusion-only fluxes and an average of 95% of 15 additional advection-plus-diffusion flux experiments. The vertical flux chamber has the capability of providing reliable measurement of VOC flux from the unsaturated zone under both diffusion and advection transport conditions.