Simulation of the effect of remediation on EDB and 1,2-DCA plumes at sites contaminated by leaded gasoline

An analytical model is used to simulate the effects of partial source removal and plume remediation on ethylene dibromide (EDB) and 1,2-dichloroethane (1,2-DCA) plumes at contaminated underground storage tank (UST) sites. The risk posed by EDB, 1,2-DCA, and commingled gasoline hydrocarbons varies throughout the plume over time. Dissolution from the light nonaqueous phase liquid (LNAPL) determines the concentration of each contaminant near the source, but biological decay in the plume has a greater influence as distance downgradient from the source increases. For this reason, compounds that exceed regulatory standards near the source may not in downgradient plume zones. At UST sites, partial removal of a residual LNAPL source mass may serve as a stand alone remedial technique if dissolved concentrations in the source zone are within several orders of magnitude of the applicable government or remedial standards. This may be the case with 1,2-DCA; however, EDB is likely to be found at concentrations that are orders of magnitude higher than its low Maximum Contaminant Level (MCL) of 0.05 μg/L (micrograms per liter). For sites with significant EDB contamination, even when plume remediation is combined with source depletion, significant timeframes may be required to mitigate the impact of this compound. Benzene and MTBE are commonly the focus of remedial efforts at UST sites, but simulations presented here suggest that EDB, and to a lesser extent 1,2-DCA, could be the critical contaminants to consider in the remediation design process at many sites.     

Plume persistence caused by back diffusion from thin clay layers in a sand aquifer following TCE source-zone hydraulic isolation

This paper concludes that back diffusion from one or a few thin clayey beds in a sand aquifer can cause contaminant persistence above MCLs in a sand aquifer long after the source zone initially causing the plume is isolated or removed. This conclusion is based on an intensive case study of a TCE contaminated site in Florida, with the processes evaluated using numerical modeling. At this site, the TCE DNAPL zone formed decades ago, and was hydraulically isolated by means of an innovative system performing groundwater extraction, treatment and re-injection. Treated water is re-injected in a row of injection wells situated a short distance downgradient of the extraction wells, creating a clean-water displacement front to efficiently flush the downgradient plume. This scheme avoids the creation of stagnation zones typical of most groundwater pump-and-treat systems, thereby minimizing the time for aquifer flushing and therefore downgradient cleanup. The system began operation in August 2002 and although the performance monitoring shows substantial declines in concentrations, detectable levels of TCE and degradation products persist downgradient of the re-injection wells, long after the TCE should have disappeared based on calculations assuming a nearly homogenous sand aquifer. Three hypotheses were assessed for this plume persistence: 1) incomplete source-zone capture, 2) DNAPL occurrence downgradient of the re-injection wells, and 3) back diffusion from one or more thin clay beds in the aquifer. After careful consideration, the first two hypotheses were eliminated, leaving back diffusion as the only plausible hypothesis, supported by detailed measurements of VOC concentrations within and near the clay beds and also by numerical model simulations that closely represent the field site hydrogeologic conditions. The model was also used to simulate a more generalized, hypothetical situation where more thin clayey beds occur in a sand aquifer with an underlying aquitard. While there is no doubt that DNAPL source mass reduction can eventually improve downgradient groundwater quality, the magnitude and time scale over which the improvement occurs is the major uncertainty given current characterization approaches. This study shows that even one thin clay bed, less than 0.2 m thick, can cause plume persistence due to back diffusion for several years or even decades after the flux from the source is completely isolated. Thin clay beds, which have a large storage capacity for dissolved and sorbed contaminant mass, are common in many types of sandy aquifers. However, without careful inspection of continuous cores and sampling, such thin clay beds, and their potential for causing long-term back-diffusion effects, can easily go unnoticed during site characterization.     

Impact of nonaqueous phase liquid (NAPL) source zone architecture on mass removal mechanisms in strongly layered heterogeneous porous media during soil vapor extraction

An existing multiphase flow simulator was modified in order to determine the effects of four mechanisms on NAPL mass removal in a strongly layered heterogeneous vadose zone during soil vapor extraction (SVE): a) NAPL flow, b) diffusion and dispersion from low permeability zones, c) slow desorption from sediment grains, and d) rate-limited dissolution of trapped NAPL. The impacts of water and NAPL saturation distribution, NAPL-type (i.e., free, residual, or trapped) distribution, and spatial heterogeneity of the permeability field on these mechanisms were evaluated. Two different initial source zone architectures (one with and one without trapped NAPL) were considered and these architectures were used to evaluate seven different SVE scenarios. For all runs, slow diffusion from low permeability zones that gas flow bypassed was a dominant factor for diminished SVE effectiveness at later times. This effect was more significant at high water saturation due to the decrease of gas-phase relative permeability. Transverse dispersion contributed to fast NAPL mass removal from the low permeability layer in both source zone architectures, but longitudinal dispersion did not affect overall mass removal time. Both slow desorption from sediment grains and rate-limited mass transfer from trapped NAPL only marginally affected removal times. However, mass transfer from trapped NAPL did affect mass removal at later time, as well as the NAPL distribution. NAPL flow from low to high permeability zones contributed to faster mass removal from the low permeability layer, and this effect increased when water infiltration was eliminated. These simulations indicate that if trapped NAPL exists in heterogeneous porous media, mass transfer can be improved by delivering gas directly to zones with trapped NAPL and by lowering the water content, which increases the gas relative permeability and changes trapped NAPL to free NAPL.     

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.     

In situ fluidization for peat bed rupture,and preliminary economic analysis

This study concerns in situ fluidization (ISF), a new remediation method with potential application to the remediation of NAPL and heavy metal contaminants, by their release from the fluidized zone generated by a water jet. The present study examines the effect of ISF on layers of peat, of significance owing to its role as an important NAPL and metal contaminant trap. Once trapped, such contaminants are not readily accessible by most remedial methods, due to the low permeability and diffusivity of the peat. A simple tank experiment is used to demonstrate rupture of a peat layer by ISF, with removal of the peat as elutriated fines and segregated peat chunks. The application of ISF in the field is then examined by three field trials in uncontaminated sands, in both saturated and unsaturated conditions. Fluidized depths of up to 1.9 m in the saturated zone (with refusal on a peat layer) and 2.5 m in the unsaturated zone (no refusal) were attained, using a 1.9-m-long, 50 mm diameter jet operated at 5-13 1 s(-1). Pulses of dark turbidity and shell fragments in the effluent indicated the rupture of peat and shelly layers. The experiments demonstrate the hydraulic viability of ISF in the field, and its ability to remove peat-based contaminants. The issues of appropriate jet design and water generation during ISF are discussed, followed by a preliminary economic analysis of ISF relative to existing remediation methods.     

A laboratory study of sediment and contaminant release during gas ebullition

Significant quantities of gas are generated from labile organic matter in contaminated sediments. The implications for the gas generation and subsequent release of contaminants from sediments are unknown but may include enhanced direct transport such as pore water advection and diffusion. The behavior of gas in sediments and the resulting migration of a polyaromatic hydrocarbon, viz phenanthrene, were investigated in an experimental system with methane injection at the base of a sediment column. Hexane above the overlying water layer was used to trap any phenanthrene migrating out of the sediment layer. The rate of suspension of solid particulate matter from the sediment bed into the overlying water layer was also monitored. The experiments indicated that significant amounts of both solid particulate matter and contaminant can be released from a sediment bed by gas movement with the amount of release related to the volume of gas released. The effective mass transfer coefficient of gas bubble-facilitated contaminant release was estimated under field conditions, being around three orders of magnitude smaller than that of bioturbation. A thin sand-capping layer (2 cm) was found to dramatically reduce the amount of contaminant or particles released with the gas because it could prevent or at least reduce sediment suspension. Based on the experimental observations, gas bubble-facilitated contaminant transport pathways for both uncapped and capped systems were proposed. Sediment cores were sliced to obtain phenanthrene concentration. X-ray computed tomography (CT) was used to investigate the void space distribution in the sediment penetrated by gas bubbles. The results showed that gas bubble migration could redistribute the sediment void spaces and may facilitate pore water circulation in the sediment.     

Laboratory studies to characterize the efficacy of sand capping a coal tar-contaminated sediment

Placement of a microbial active sand cap on a coal tar-contaminated river sediment has been suggested as a cost effective remediation strategy. This approach assumes that the flux of contaminants from the sediment is sufficiently balanced by oxygen and nutrient fluxes into the sand layer such that microbial activity will reduce contaminant concentrations within the new benthic zone and reduce the contaminant flux to the water column. The dynamics of such a system were evaluated using batch and column studies with microbial communities from tar-contaminated sediment under different aeration and nutrient inputs. In a 30-d batch degradation study on aqueous extracts of coal tar sediment, oxygen and nutrient concentrations were found to be key parameters controlling the degradation rates of polycyclic aromatic hydrocarbons (PAHs). For the five PAHs monitored (naphthalene, fluorene, phenanthrene, anthracene, and pyrene), degradation rates were inversely proportional to molecular size. For the column studies, where three columns were packed with a 20-cm sand layer on the top of a 5 cm of sediment layer, flow was established to sand layers with (1) aerated water, (2) N(2) sparged water, or (3) HgCl(2)-sterilized N(2) sparged water. After steady-state conditions, PAH concentrations in effluents were the lowest in the aerated column, except for pyrene, whose concentration was invariant with all effluents. These laboratory scale studies support that if sufficient aeration can be achieved in the field through either active and passive means, the resulting microbially active sand layer can improve the water quality of the benthic zone and reduce the flux of many, but not all, PAHs to the water column.     

An analytical solution to contaminant diffusion in semi-infinite clayey soils with piecewise linear adsorption

The adsorption of contaminants onto soil particles typically is nonlinear if the contaminant concentration is sufficiently high. A simplified piecewise linear adsorption isotherm consistent with experimental results is proposed as an approximation for nonlinear adsorption behavior. This approximation allows for the use of analytical solution to model solute diffusion of contaminants that exhibit nonlinear adsorption. A moving boundary is introduced to represent significant changes in the retardation factor of clay with an increase in solute concentration. The proposed analytical solutions were validated using experimental data presented in the literature. There is negligible difference between the results obtained by the proposed analytical solution and those obtained by the linear model when C_m/C₀ reached 0.5. The results also show that the model based on linear adsorption using the initial secant of the Freundlich isotherm leads to significantly lower estimated breakthrough time for the contaminant of interest than that obtained using the proposed model. The earlier breakthrough is due to an under-estimation of the amount of adsorption. The proposed method is relatively simple to apply and can be used for evaluating experimental results and verifying more complex numerical models.     

An analytical solution for vertical transport of volatile chemicals in the vadose zone

An analytical solution is presented for one-dimensional vertical transport of volatile chemicals through the vadose zone to groundwater. The solution accounts for the important transport mechanisms of the steady advection of water and gas, diffusion and dispersion in water and gas, as well as adsorption, and first-order degradation. By assuming a linear, equilibrium partitioning between water, gas and the adsorbed chemical phases, the dependent variable in the mathematical model becomes the total resident concentration. The general solution was derived for cases having a constant initial total concentration over a discrete depth interval and a zero initial total concentration elsewhere. A zero concentration gradient is assumed at the groundwater table. Examples are given to demonstrate the application of the new solution for calculating the case of a non-uniform initial source concentration, and estimating the transport of chemicals to the groundwater and the atmosphere. The solution was also used to verify a numerical code called VLEACH. We discovered an error in VLEACH, and found that the new solution agreed very well with the numerical results by corrected VLEACH. A simplified solution to predict the migration of volatile organic chemical due to the gas density effect has shown that a high source concentration may lead to significant downward advective gas-phase transport in a soil with a high air-permeability.     

Concurrent bioremediation of perchlorate and 1,1,1-trichloroethane in an emulsified oil barrier

A detailed field pilot test was conducted to evaluate the use of edible oil emulsions for enhanced in situ biodegradation of perchlorate and chlorinated solvents in groundwater. Edible oil substrate (EOS) was injected into a line of ten direct push injection wells over a 2-day period to form a 15-m-long biologically active permeable reactive barrier (bio-barrier). Field monitoring results over a 2.5-year period indicate the oil injection generated strongly reducing conditions in the oil-treated zone with depletion of dissolved oxygen, nitrate, and sulfate, and increases in dissolved iron, manganese and methane. Perchlorate was degraded from 3100 to 20,000 microg/L to below detection (<4 microg/L) in the injection and nearby monitor wells within 5 days following the injection. Two years after the single emulsion injection, perchlorate was less than 6 microg/L in every downgradient well compared to an average upgradient concentration of 13,100 microg/L. Immediately after emulsion injection, there were large shifts in concentrations of chlorinated solvents and degradation products due to injection of clean water, sorption to the oil and adaptation of the in situ microbial community. Approximately 4 months after emulsion injection, concentrations of 1,1,1-trichloroethane (TCA), perchloroethene (PCE), trichloroethene (TCE) and their degradation products appeared to reach a quasi steady-state condition. During the period from 4 to 18 months, TCA was reduced from 30-70 microM to 0.2-4 microM during passage through the bio-barrier. However, 1-9 microM 1,1-dichloroethane (DCA) and 8-14 microM of chloroethane (CA) remained indicating significant amounts of incompletely degraded TCA were discharging from the oil-treated zone. During this same period, PCE and TCE were reduced with concurrent production of 1,2-cis-dichloroethene (cis-DCE). However, very little VC or ethene was produced indicating reductive dechlorination slowed or stopped at cis-DCE. The incomplete removal of TCA, PCE and TCE is likely associated with the short (5-20 days) hydraulic retention time of contaminants in the oil-treated zone. The permeability of the injection wells declined by 39-91% (average=68%) presumably due to biomass growth and/or gas production. However, non-reactive tracer tests and detailed monitoring of the perchlorate plume demonstrated that the permeability loss did not result in excessive flow bypassing around the bio-barrier. Contaminant transport and degradation within the bio-barrier was simulated using an advection-dispersion-reaction model where biodegradation rate was assumed to be linearly proportional to the residual oil concentration (Soil) and the contaminant concentration. Using this approach, the calibrated model was able to closely match the observed contaminant distribution. The calibrated model was then used to design a full-scale barrier to treat both ClO4 and chlorinated solvents.     

The effect of surface-active solutes on water flow and contaminant transport in variably saturated porous media with capillary fringe effects

Organic contaminants that decrease the surface tension of water (surfactants) can have an effect on unsaturated flow through porous media due to the dependence of capillary pressure on surface tension. We used an intermediate-scale 2D flow cell (2.44 x 1.53 x 0.108 m) packed with a fine silica sand to investigate surfactant-induced flow perturbations. Surfactant solution (7% 1-butanol and dye tracer) was applied at a constant rate at a point source located on the soil surface above an unconfined synthetic aquifer with ambient groundwater flow and a capillary fringe of approximately 55 cm. A glass plate allowed for visual flow and transport observations. Thirty instrumentation stations consist of time domain reflectometry probes and tensiometers measured in-situ moisture content and pressure head, respectively. As surfactant solution was applied at the point source, a transient flow perturbation associated with the advance of the surfactant solution was observed. Above the top of the capillary fringe the advance of the surfactant solution caused a visible drainage front that radiated from the point source. Upon reaching the capillary fringe, the drainage front caused a localized depression of the capillary fringe below the point source because the air-entry pressure decreased in proportion to the decrease in surface tension caused by the surfactant. Eventually, a new capillary fringe height was established. The height of the depressed capillary fringe was proportional to height of the initial capillary fringe multiplied by the relative surface tension of the surfactant solution. The horizontal transport of surfactant in the depressed capillary fringe, driven primarily by the ambient groundwater flow, caused the propagation of a wedge-shaped drying front in the downgradient direction. Comparison of dye transport during the surfactant experiment to dye transport in an experiment without surfactant indicated that because surfactant-induced drainage decreased the storage capacity of the vadose zone, the dye breakthrough time to the water table was more than twice as fast when the contaminant solution contained surfactant. The extensive propagation of the drying front and the effect of vadose zone drainage on contaminant breakthrough time suggest the importance of considering surface tension effects on unsaturated flow and transport in systems containing surface-active organic contaminants or systems, where surfactants are used for remediation of the vadose zone or unconfined aquifers.     

Effect of pH control at the anode for the electrokinetic removal of phenanthrene from kaolin soil

Polycyclic aromatic hydrocarbon (PAH)-contaminated soils exist at numerous sites, and these sites may threaten public health and the environment because many PAH compounds are toxic, mutagenic, and/or carcinogenic. PAHs are also hydrophobic and persistent, so conventional remediation methods are often costly or inefficient, especially when the contaminants are present in low permeability and/or organic soils. An innovative technique, electrokinetically enhanced in situ flushing, has the potential to increase soil-solution-contaminant interaction and PAH removal efficiency for low permeability soils; however, the electrolysis reaction at the anode may adversely affect the remediation of low acid buffering capacity soils, such as kaolin. Therefore, the objective of this study was to improve the remediation of low acid buffering soils by controlling the pH at the anode to counteract the electrolysis reaction. Six bench-scale electrokinetic experiments were conducted, where each test employed one of three different flushing solutions, deionized water, a surfactant, or a cosolvent. For each of these solutions, tests were performed with and without a 0.01 M NaOH solution at the anode to control the pH. The test using deionized water with pH control generated a higher electroosmotic flow than the equivalent test performed without pH control, but the electroosmotic flow difference between the surfactant and cosolvent tests with and without pH control was minor compared to that observed with the deionized water tests. Controlling the pH was beneficial for increasing contaminant solubilization and migration from the soil region adjacent to the anode, but the high contaminant concentrations that resulted in the middle or cathode soil regions indicates that subsequent changes in the soil and/or solution chemistry caused contaminant deposition and low overall contaminant removal efficiency.     

Methane oxidation in a crude oil contaminated aquifer: Delineation of aerobic reactions at the plume fringes

High resolution direct-push profiling over short vertical distances was used to investigate CH(4) attenuation in a petroleum contaminated aquifer near Bemidji, Minnesota. The contaminant plume was delineated using dissolved gases, redox sensitive components, major ions, carbon isotope ratios in CH(4) and CO(2), and the presence of methanotrophic bacteria. Sharp redox gradients were observed near the water table. Shifts in δ(13)C(CH4) from an average of -57.6‰ (±1.7‰) in the methanogenic zone to -39.6‰ (±8.7‰) at 105m downgradient, strongly suggest CH(4) attenuation through microbially mediated degradation. In the downgradient zone the aerobic/anaerobic transition is up to 0.5m below the water table suggesting that transport of O(2) across the water table is leading to aerobic degradation of CH(4) at this interface. Dissolved N(2) concentrations that exceeded those expected for water in equilibrium with the atmosphere indicated bubble entrapment followed by preferential stripping of O(2) through aerobic degradation of CH(4) or other hydrocarbons. Multivariate and cluster analysis were used to distinguish between areas of significant bubble entrapment and areas where other processes such as the infiltration of O(2) rich recharge water were important O(2) transport mechanisms.     

An assessment of natural biotransformation of petroleum hydrocarbons and chlorinated solvents at an aquifer plume transect

Field biogeochemical characterization and laboratory microcosm studies were performed to assess the potential for future biotransformation of trichloroethylene (TCE) and toluene in a plume containing petroleum hydrocarbons and chlorinated solvents at the former Wurtsmith Air Force Base in Oscoda, MI. In situ terminal electron accepting processes (TEAPs), contaminant composition and microbial phylogeny were studied at a plume transect 100 m downgradient of the source. The presence of reduced electron acceptors, relevant microbial communities, and elevated dissolved methane and carbon dioxide concentrations at the transect, as well as downgradient accumulation of BTEX metabolites and dechlorination products, indicated that past or current reductive dechlorination at the transect was likely driven by BTEX biodegradation in the methanogenic zone. However, TCE and toluene mineralization in sediment-groundwater microcosms without added electron acceptors did not exceed 5% during 300 days of incubation and was nearly invariant with original sediment TEAP, even following amendments of nitrogen and phosphorus. Mineralization rates were on the order of 0.0015-0.03 mumol/g day. After 8 months, microcosms showed evidence of methanogenesis, but CH4 and CO2 production arose from the degradation of contaminants other than toluene. Cis-dichloroethylene was observed in only one methanogenic microcosm after more than 500 days. It appears likely that spatially and temporally dynamic redox zonation at the plume transect will prevent future sustained reductive dehalogenation of highly chlorinated solvents, for during the course of a year, the predominant TEAP at the highly contaminated water table shifted from methanogenesis to iron- and sulfate-reduction. It is recommended that biotransformation studies combine considerations of long-term, spatially relevant changes in redox zonation with laboratory-scale studies of electron donor utilization and cometabolic substrate transformation to yield a more accurate assessment of natural bioattenuation of specific pollutants in aquifers contaminated by undefined organic waste mixtures.     

Progression of natural attenuation processes at a crude oil spill site: II. Controls on spatial distribution of microbial populations

A multidisciplinary study of a crude-oil contaminated aquifer shows that the distribution of microbial physiologic types is strongly controlled by the aquifer properties and crude oil location. The microbial populations of four physiologic types were analyzed together with permeability, pore-water chemistry, nonaqueous oil content, and extractable sediment iron. Microbial data from three vertical profiles through the anaerobic portion of the contaminated aquifer clearly show areas that have progressed from iron-reduction to methanogenesis. These locations contain lower numbers of iron reducers, and increased numbers of fermenters with detectable methanogens. Methanogenic conditions exist both in the area contaminated by nonaqueous oil and also below the oil where high hydrocarbon concentrations correspond to local increases in aquifer permeability. The results indicate that high contaminant flux either from local dissolution or by advective transport plays a key role in determining which areas first become methanogenic. Other factors besides flux that are important include the sediment Fe(II) content and proximity to the water table. In locations near a seasonally oscillating water table, methanogenic conditions exist only below the lowest typical water table elevation. During 20 years since the oil spill occurred, a laterally continuous methanogenic zone has developed along a narrow horizon extending from the source area to 50-60 m downgradient. A companion paper [J. Contam. Hydrol. 53, 369-386] documents how the growth of the methanogenic zone results in expansion of the aquifer volume contaminated with the highest concentrations of benzene, toluene, ethylbenzene, and xylenes.     

Competing TCE and cis-DCE degradation kinetics by zero-valent iron-experimental results and numerical simulation

The successful dechlorination of mixtures of chlorinated hydrocarbons with zero-valent metals requires information concerning the kinetics of simultaneous degradation of different contaminants. This includes intraspecies competitive effects (loading of the reactive iron surface by a single contaminant) as well as interspecies competition of several contaminants for the reactive sites available. In columns packed with zero-valent iron, the degradation behaviour of trichloroethylene (TCE), cis-dichloroethylene (DCE) and mixtures of both was measured in order to investigate interspecies competition. Although a decreasing rate of dechlorination is to be expected, when several degradable substances compete for the reactive sites on the iron surface, TCE degradation is nearly unaffected by the presence of cis-DCE. In contrast, cis-DCE degradation rates decrease significantly when TCE is added. A new modelling approach is developed in order to identify and quantify the observed competitive effects. The numerical model TBC (Transport, Biochemistry and Chemistry, Sch?fer et al., 1998a) is used to describe adsorption, desorption and dechlorination in a mechanistic way. Adsorption and degradation of a contaminant based on a limited number of reactive sites leads to a combined zero- and first-order degradation kinetics for high and low concentrations, respectively. The adsorption of several contaminants with different sorption parameters to a limited reactive surface causes interspecies competition. The reaction scheme and the parameters required are successfully transferred from Arnold and Roberts (2000b) to the model TBC. The degradation behaviour of the mixed contamination observed in the column experiments can be related to the adsorption properties of TCE and cis-DCE. By predicting the degradation of the single substances TCE and cis-DCE as well as mixtures of both, the calibrated model is used to investigate the effects of interspecies competition on the design of permeable reactive iron barriers. Even if TCE is present in only small concentrations (>3% of molar cis-DCE concentration) it is the contaminant limiting the residence time and the required thickness of the iron barrier.     

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.     

Monitoring groundwater contamination and delineating source zones at industrial sites: uncertainty analyses using integral pumping tests

Field-scale characterisations of contaminant plumes in groundwater, as well as source zone delineations, are associated with uncertainties that can be considerable. A major source of uncertainty in environmental datasets is due to variability of sampling results, as a direct consequence of the heterogeneity of environmental matrices. We develop a methodology for quantifying uncertainties in field-scale mass flow and average concentration estimations, using integral pumping tests (IPTs), where the contaminant concentration is measured as a function of time in a pumping well. This procedure increases the sampling volume and reduces the effect of small-scale variability that may bias point-scale measurements. In particular, using IPTs, the interpolation uncertainty of conventional point-scale measurements is transformed to a quantifiable uncertainty related to the (unknown) plume position relative to the pumping well. We show that this plume position uncertainty generally influenced the predicted mass flows and average concentrations (of acenapthene, benzene and CHCs) to a greater extent than a boundary condition uncertainty related to the local water balance, considering 19 control planes at a highly heterogeneous industrial site in southwest Germany. Furthermore, large (order of magnitude) uncertainties only occurred if the conditions were strongly heterogeneous in the nearest vicinity of the well. We also develop a consistent methodology for an assessment of the combined effect of uncertainty in hydraulic conditions and uncertainty in reactive transport parameters for delimiting of both contaminant source zones and zones absent of source, based on (downgradient) IPTs.     

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.     

A hybrid method for inverse characterization of subsurface contaminant flux

The methods presented in this work provide a potential tool for characterizing contaminant source zones in terms of mass flux. The problem was conceptualized by considering contaminant transport through a vertical "flux plane" located between a source zone and a downgradient region where contaminant concentrations were measured. The goal was to develop a robust method capable of providing a statement of the magnitude and uncertainty associated with estimated contaminant mass flux values. In order to estimate the magnitude and transverse spatial distribution of mass flux through a plane, the problem was considered in an optimization framework. Two numerical optimization techniques were applied, simulated annealing (SA) and minimum relative entropy (MRE). The capabilities of the flux plane model and the numerical solution techniques were evaluated using data from a numerically generated test problem and a nonreactive tracer experiment performed in a three-dimensional aquifer model. Results demonstrate that SA is more robust and converges more quickly than MRE. However, SA is not capable of providing an estimate of the uncertainty associated with the simulated flux values. In contrast, MRE is not as robust as SA, but once in the neighborhood of the optimal solution, it is quite effective as a tool for inferring mass flux probability density functions, expected flux values, and confidence limits. A hybrid (SA-MRE) solution technique was developed in order to take advantage of the robust solution capabilities of SA and the uncertainty estimation capabilities of MRE. The coupled technique provided probability density functions and confidence intervals that would not have been available from an independent SA algorithm and they were obtained more efficiently than if provided by an independent MRE algorithm.