Humic acid enhanced remediation of an emplaced diesel source in groundwater. 2. Numerical model development and application

A pilot scale experiment for humic acid-enhanced remediation of diesel fuel, described in Part 1 of this series, is numerically simulated in three dimensions. Groundwater flow, enhanced solubilization of the diesel source, and reactive transport of the dissolved contaminants and humic acid carrier are solved with a finite element Galerkin approach. The model (BIONAPL) is calibrated by comparing observed and simulated concentrations of seven diesel fuel components (BTEX and methyl-, dimethyl- and trimethylnaphthalene) over a 1500-day monitoring period. Data from supporting bench scale tests were used to estimate contaminant-carrier binding coefficients and to simulate two-site sorption of the carrier to the aquifer sand. The model accurately reproduced the humic acid-induced 10-fold increase in apparent solubility of trimethylnaphthalene. Solubility increases on the order of 2-5 were simulated for methylnaphthalene and dimethylnaphthalene, respectively. Under the experimental and simulated conditions, the residual 500-ml diesel source was almost completely dissolved and degraded within 5 years. Without humic acid flushing, the simulations show complete source dissolution would take about six times longer.     

Mechanisms of electron acceptor utilization: implications for simulating anaerobic biodegradation

Simulation of biodegradation reactions within a reactive transport framework requires information on mechanisms of terminal electron acceptor processes (TEAPs). In initial modeling efforts, TEAPs were approximated as occurring sequentially, with the highest energy-yielding electron acceptors (e.g. oxygen) consumed before those that yield less energy (e.g., sulfate). Within this framework in a steady state plume, sequential electron acceptor utilization would theoretically produce methane at an organic-rich source and Fe(II) further downgradient, resulting in a limited zone of Fe(II) and methane overlap. However, contaminant plumes often display much more extensive zones of overlapping Fe(II) and methane. The extensive overlap could be caused by several abiotic and biotic processes including vertical mixing of byproducts in long-screened monitoring wells, adsorption of Fe(II) onto aquifer solids, or microscale heterogeneity in Fe(III) concentrations. Alternatively, the overlap could be due to simultaneous utilization of terminal electron acceptors. Because biodegradation rates are controlled by TEAPs, evaluating the mechanisms of electron acceptor utilization is critical for improving prediction of contaminant mass losses due to biodegradation. Using BioRedox-MT3DMS, a three-dimensional, multi-species reactive transport code, we simulated the current configurations of a BTEX plume and TEAP zones at a petroleum-contaminated field site in Wisconsin. Simulation results suggest that BTEX mass loss due to biodegradation is greatest under oxygen-reducing conditions, with smaller but similar contributions to mass loss from biodegradation under Fe(III)-reducing, sulfate-reducing, and methanogenic conditions. Results of sensitivity calculations document that BTEX losses due to biodegradation are most sensitive to the age of the plume, while the shape of the BTEX plume is most sensitive to effective porosity and rate constants for biodegradation under Fe(III)-reducing and methanogenic conditions. Using this transport model, we had limited success in simulating overlap of redox products using reasonable ranges of parameters within a strictly sequential electron acceptor utilization framework. Simulation results indicate that overlap of redox products cannot be accurately simulated using the constructed model, suggesting either that Fe(III) reduction and methanogenesis are occurring simultaneously in the source area, or that heterogeneities in Fe(III) concentration and/or mineral type cause the observed overlap. Additional field, experimental, and modeling studies will be needed to address these questions.     

Progression of natural attenuation processes at a crude-oil spill site: I. Geochemical evolution of the plume

A 16-year study of a hydrocarbon plume shows that the extent of contaminant migration and compound-specific behavior have changed as redox reactions, most notably iron reduction, have progressed over time. Concentration changes at a small scale, determined from analysis of pore-water samples drained from aquifer cores, are compared with concentration changes at the plume scale, determined from analysis of water samples from an observation well network. The small-scale data show clearly that the hydrocarbon plume is growing slowly as sediment iron oxides are depleted. Contaminants, such as ortho-xylene that appeared not to be moving downgradient from the oil on the basis of observation well data, are migrating in thin layers as the aquifer evolves to methanogenic conditions. However, the plume-scale observation well data show that the downgradient extent of the Fe2+ and BTEX plume did not change between 1992 and 1995. Instead, depletion of the unstable Fe (III) oxides near the subsurface crude-oil source has caused the maximum dissolved iron concentration zone within the plume to spread at a rate of approximately 3 m/year. The zone of maximum concentrations of benzene, toluene, ethylbenzene and xylene (BTEX) has also spread within the anoxic plume. In monitoring the remediation of hydrocarbon-contaminated ground water by natural attenuation, subtle concentration changes in observation well data from the anoxic zone may be diagnostic of depletion of the intrinsic electron-accepting capacity of the aquifer. Recognition of these subtle patterns may allow early prediction of growth of the hydrocarbon plume.     

Effect of water-table fluctuation on dissolution and biodegradation of a multi-component, light nonaqueous-phase liquid

Light nonaqueous-phase liquids (LNAPLs) such as gasoline and diesel fuel are among the most common causes of soil and groundwater contamination. Dissolution and subsequent advective transport of LNAPL components can negatively impact water supplies, while biodegradation is thought to be an important sink for this class of contaminants. We present a laboratory investigation of the effect of a water-table fluctuation on dissolution and biodegradation of a multi-component LNAPL (85% hexadecane, 5% toluene, 5% ethylbenzene, and 5% 2-methylnapthalene on a molar basis) in a pair of similar model aquifers (80 cm x 50 cm x 3 cm), one of which was subjected to a water-table fluctuation. Water-table fluctuation resulted in LNAPL and air entrapment below the water table, an increase in the vertical extent of the LNAPL source zone (by factor 6.7), and an increase in the volume of water passing through the source zone (by factor ~18). Effluent concentrations of dissolved LNAPL components were substantially higher and those of dissolved nitrate lower in the model aquifer where a fluctuation had been induced. Thus, water-table fluctuation led to enhanced biodegradation activity (28.3 mmol of nitrate consumed compared to 16.3 mmol in the model without fluctuation) as well as enhanced dissolution of LNAPL components. Despite the increased biodegradation, fluctuation led to increased elution of dissolved LNAPL components from the system (by factors 10-20). Hence, water-table fluctuations in LNAPL-contaminated aquifers might be expected to result in increased exposure of downgradient receptors to LNAPL components. Accordingly, water-table fluctuations in contaminated aquifers are probably undesirable unless the LNAPL is of minimal solubility or the dissolved-phase plume is not expected to reach a receptor due to distance or the presence of some form of containment.     

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.     

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.     

Fingerprinting and Groundwater Model Application to Evaluate Hydraulic Barrier Efficiency (Italy)

Remediation actions at contaminated sites are based on multiple numerical model scenarios considering different parameter distributions, source positions and contaminant transport paths. In some cases the excess of scenarios is due to uncertainties in the conceptual model as a result of the spread of contamination through heterogeneities in the physical system. Reduction of project hypotheses and conceptual model uncertainty is therefore needed. This result can be achieved by coupling hydrogeological investigations with environmental forensic techniques, better localization of the source and understanding of contamination history. In this respect, in the present study, compositional fingerprinting and groundwater flow modeling were applied to a former oil storage facility where, even though a hydraulic barrier had been built to stop the hydrocarbon plume, the presence of some hydrocarbons was still found in downgradient monitoring wells. The final aim was to evaluate the efficacy of the hydraulic barrier and identify of the source of pollution. Fingerprinting results indicated pollution with a gasoline-diesel mixture much altered by water washing and/or biodegradation. Comparison of seven groundwater samples collected in wells and monitoring wells was performed by analyzing the volatile fraction (BTEX) and the total ion chromatogram (TIC), focusing attention on: n-alkanes (m/z 85), alkylcyclohexanes (m/z 83), isoprenoids (m/z 113), C4-alkylbenzenes (m/z 134), C3-C6 alkylbenzenes and polycyclic aromatic hydrocarbons (PAHs). The most probable scenario was then identified by combining the results of fingerprinting with different contaminant paths obtained using the numerical model.     

Processes controlling the distribution and natural attenuation of dissolved phenolic compounds in a deep sandstone aquifer

Processes controlling the distribution and natural attenuation (NA) of phenol, cresols and xylenols released from a former coal-tar distillation plant in a deep Triassic sandstone aquifer are evaluated from vertical profiles along the plume centerline at 130 and 350 m from the site. Up to four groups of contaminants (phenols, mineral acids, NaOH, NaCl) form discrete and overlapping plumes in the aquifer. Their distribution reflects changing source history with releases of contaminants from different locations. Organic contaminant distribution in the aquifer is determined more by site source history than degradation. Contaminant degradation at total organic carbon (TOC) concentrations up to 6500 mg l(-1) (7500 mg l(-1) total phenolics) is occurring by aerobic respiration NO3-reduction, Mn(IV)-/Fe(III)-reduction, SO4-reduction, methanogenesis and fermentation, with the accumulation of inorganic carbon, organic metabolites (4-hydroxybenzaldehyde, 4-hydroxybenzoic acid), acetate, Mn(II), Fe(II), S(-II), CH4 and H2 in the plume. Aerobic and NO3-reducing processes are restricted to a 2-m-thick plume fringe but Mn(IV)-/Fe(II)-reduction, SO4-reduction, methanogenesis and fermentation occur concomitantly in the plume. Dissolved H2 concentrations in the plume vary from 0.7 to 110 nM and acetate concentrations reach 200 mg l(-1). The occurrence of a mixed redox system and concomitant terminal electron accepting processes (TEAPs) could be explained with a partial equilibrium model based on the potential in situ free energy (deltaGr) yield for oxidation of H2 by specific TEAPs. Respiratory processes rather than fermentation are rate limiting in determining the distribution of H2 and TEAPs and H2 dynamics in this system. Most (min. 90%) contaminant degradation has occurred by aerobic and NO3-reducing processes at the plume fringe. This potential is determined by the supply of aqueous O2 and NO3 from uncontaminated groundwater, as controlled by transverse mixing, which is limited in this aquifer by low dispersion. Consumption to date of mineral oxides and SO4 is, respectively, <0.15% and 0.4% of the available aquifer capacity, and degradation using these oxidants is <10%. Fermentation is a significant process in contaminant turnover, accounting for 21% of degradation products present in the plume, and indicating that microbial respiration rates are slow in comparison with fermentation. Under present conditions, the potential for degradation in the plume is very low due to inhibitory effects of the contaminant matrix. Degradation products correspond to <22% mass loss over the life of the plume, providing a first-order plume scale half-life >140 years. The phenolic compounds are biodegradable under the range of redox conditions in the aquifer and the aquifer is not oxidant limited, but the plume is likely to be long-lived and to expand. Degradation is likely to increase only after contaminant concentrations are reduced and aqueous oxidant inputs are increased by dispersion of the plume. The results imply that transport processes may exert a greater control on the natural attenuation of this plume than aquifer oxidant availability.     

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.     

Numerical experiments and field results on the size of steady state plumes

Contaminated groundwater poses a serious risk for drinking water supplies. Under certain conditions, however, groundwater contamination remains restricted to a tolerable extent because of natural attenuation processes. We present an innovative approach to evaluate the size of these so-called steady-state plumes by 2-D and 1-D modelling in homogeneous aquifers. If longitudinal mixing is negligible, scenarios can be modelled in a simplified way using a 1-D domain vertical to the direction of flow. We analysed the sensitivity of the plume length with respect to biodegradation kinetics, flow velocity, transverse vertical dispersivity alphat, the source and aquifer geometry and reaction stoichiometry. Our findings indicate that for many readily biodegradable compounds transverse-dispersive mixing rather than reaction kinetics is the limiting factor for natural attenuation. Therefore, if alphat, aquifer and source geometry and concentrations of electron acceptors and donors are known, the length of the steady state contaminant plume can be predicted. The approach is validated under field conditions for an ammonium plume at a former landfill site in SW Germany.     

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.     

Mineralogical compositions of aquifer matrix as necessary initial conditions in reactive contaminant transport models

Mineralogical compositions and their spatial distributions are important initial conditions for reactive transport modeling. However, popular Kd-based "reactive" transport models only require contaminant concentrations in the pore fluids as initial conditions, and minerals implicitly represent infinite sources and sinks in these models. That situation results in a general neglect of mineralogical characterization in site investigations. This study uses a coupled multi-component reactive mass transport model to predict the natural attenuation of a ground water plume at a uranium mill tailings site in western USA. Numerous ground water geochemistry data are available at this site, but mineralogical data are sketchy. Even given the well-defined pore fluid chemistry, variations of secondary mineral species and mineral abundances in the aquifer resulted in significantly different modeling outcomes. Results show that the amount of calcite in the aquifer determines the distances of plume migration. The possible presence of jurbanite, an aluminum sulfate phase, can store acidity temporarily but cause more severe contamination on a later date. The surfaces of iron oxyhydroxides can store significant amounts of sulfate and protons and serve as a second source for prolonged contamination. These simulations under field conditions illustrate that mineralogical compositions are an essential requirement for accurate prediction of contaminant fate and transport.     

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.     

Using the Gasoline Additive MTBE in Forensic Environmental Investigations

A variety of additives are used in gasoline, and they can sometimes be used to help identify the source, timing, or number of gasoline spills at a site. The physicochemical characteristics of the additive MTBE, and its historical use pattern in the United States since 1979, make it a key compound to study when conducting forensic investigations of gasoline spills. MTBE's low octanol: water distribution coefficient and high solubility cause it to dissolve into groundwater more readily than other gasoline components. Thus, the initial appearance of MTBE in the groundwater is often a good indicator of a recent gasoline spill. MTBE's very low retardation and minimal biodegradation in groundwater can be used with transport rate calculations to establish relatively accurate estimates of spill timing. Because MTBE moves faster in groundwater than BTEX compounds, if a gasoline spill site has a BTEX plume that is longer than the MTBE plume, it is certain that at least two distinctly different gasoline releases have occurred. This allows for the identification of new gasoline spills, even when substantial subsurface petroleum contamination already exists. An example application is reviewed to demonstrate the use of MTBE data in forensic investigations.     

Modelling of physical and reactive processes during biodegradation of a hydrocarbon plume under transient groundwater flow conditions

Numerical experiments of non-reactive and reactive transport were carried out to quantify the influence of a seasonally varying, transient flow field on transport and natural attenuation at a hydrocarbon-contaminated field site. Different numerical schemes for solving advective transport were compared to assess their capability to model low transversal dispersivities in transient flow fields. For the field site, it is shown that vertical plume spreading is largely inhibited, particularly if sorption is taken into account. For the reactive simulations, a biodegradation reaction module for the geochemical transport model PHT3D was developed. Results of the reactive transport simulations show that under the site-specific conditions the temporal variations in groundwater flow do, to a modest extent, affect average biodegradation rates and average total (dissolved) contaminant mass in the aquifer. The model simulations demonstrate that the seasonal variability in groundwater flow only results in significantly enhanced biodegradation rates when a differential sorption of electron donor (toluene) and electron acceptor (sulfate) is assumed.     

Using the Gasoline Additive MTBE in Forensic Environmental Investigations

A variety of additives are used in gasoline, and they can sometimes be used to help identify the source, timing, or number of gasoline spills at a site. The physicochemical characteristics of the additive MTBE, and its historical use pattern in the United States since 1979, make it a key compound to study when conducting forensic investigations of gasoline spills. MTBE's low octanol : water distribution coefficient and high solubility cause it to dissolve into groundwater more readily than other gasoline components. Thus, the initial appearance of MTBE in the groundwater is often a good indicator of a recent gasoline spill. MTBE's very low retardation and minimal biodegradation in groundwater can be used with transport rate calculations to establish relatively accurate estimates of spill timing. Because MTBE moves faster in groundwater than BTEX compounds, if a gasoline spill site has a BTEX plume that is longer than the MTBE plume, it is certain that at least two distinctly different gasoline releases have occurred. This allows for the identification of new gasoline spills, even when substantial subsurface petroleum contamination already exists. An example application is reviewed to demonstrate the use of MTBE data in forensic investigations.     

Biodegradation potential of MTBE in a fractured chalk aquifer under aerobic conditions in long-term uncontaminated and contaminated aquifer microcosms

The potential for aerobic biodegradation of MTBE in a fractured chalk aquifer is assessed in microcosm experiments over 450 days, under in situ conditions for a groundwater temperature of 10 °C, MTBE concentration between 0.1 and 1.0 mg/L and dissolved O₂ concentration between 2 and 10 mg/L. Following a lag period of up to 120 days, MTBE was biodegraded in uncontaminated aquifer microcosms at concentrations up to 1.2 mg/L, demonstrating that the aquifer has an intrinsic potential to biodegrade MTBE aerobically. The MTBE biodegradation rate increased three-fold from a mean of 6.6 ± 1.6 μg/L/day in uncontaminated aquifer microcosms for subsequent additions of MTBE, suggesting an increasing biodegradation capability, due to microbial cell growth and increased biomass after repeated exposure to MTBE. In contaminated aquifer microcosms which also contained TAME, MTBE biodegradation occurred after a shorter lag of 15 or 33 days and MTBE biodegradation rates were higher (max. 27.5 μg/L/day), probably resulting from an acclimated microbial population due to previous exposure to MTBE in situ. The initial MTBE concentration did not affect the lag period but the biodegradation rate increased with the initial MTBE concentration, indicating that there was no inhibition of MTBE biodegradation related to MTBE concentration up to 1.2 mg/L. No minimum substrate concentration for MTBE biodegradation was observed, indicating that in the presence of dissolved O₂ (and absence of inhibitory factors) MTBE biodegradation would occur in the aquifer at MTBE concentrations (ca. 0.1 mg/L) found at the front of the ether oxygenate plume. MTBE biodegradation occurred with concomitant O₂ consumption but no other electron acceptor utilisation, indicating biodegradation by aerobic processes only. However, O₂ consumption was less than the stoichiometric requirement for complete MTBE mineralization, suggesting that only partial biodegradation of MTBE to intermediate organic metabolites occurred. The availability of dissolved O₂ did not affect MTBE biodegradation significantly, with similar MTBE biodegradation behaviour and rates down to ca. 0.7 mg/L dissolved O₂ concentration. The results indicate that aerobic MTBE biodegradation could be significant in the plume fringe, during mixing of the contaminant plume and uncontaminated groundwater and that, relative to the plume migration, aerobic biodegradation is important for MTBE attenuation. Moreover, should the groundwater dissolved O₂ concentration fall to zero such that MTBE biodegradation was inhibited, an engineered approach to enhance in situ bioremediation could supply O₂ at relatively low levels (e.g. 2–3 mg/L) to effectively stimulate MTBE biodegradation, which has significant practical advantages. The study shows that aerobic MTBE biodegradation can occur at environmentally significant rates in this aquifer, and that long-term microcosm experiments (100s days) may be necessary to correctly interpret contaminant biodegradation potential in aquifers to support site management decisions.     

Hydrochemical and isotopic effects associated with petroleum fuel biodegradation pathways in a chalk aquifer

Hydrochemical data, compound specific carbon isotope analysis and isotopic enrichment trends in dissolved hydrocarbons and residual electron acceptors have been used to deduce BTEX and MTBE degradation pathways in a fractured chalk aquifer. BTEX compounds are mineralised sequentially within specific redox environments, with changes in electron acceptor utilisation being defined by the exhaustion of specific BTEX components. A zone of oxygen and nitrate exhaustion extends approximately 100 m downstream from the plume source, with residual sulphate, toluene, ethylbenzene and xylene. Within this zone complete removal of the TEX components occurs by bacterial sulphate reduction, with sulphur and oxygen isotopic enrichment of residual sulphate (epsilon(s) = -14.4 per thousand to -16.0 per thousand). Towards the plume margins and at greater distance along the plume flow path nitrate concentrations increase with delta15N values of up to +40 per thousand indicating extensive denitrification. Benzene and MTBE persist into the denitrification zone, with carbon isotope enrichment of benzene indicating biodegradation along the flow path. A Rayleigh kinetic isotope enrichment model for 13C-enrichment of residual benzene gives an apparent epsilon value of -0.66 per thousand. MTBE shows no significant isotopic enrichment (delta13C = -29.3 per thousand to -30.7 per thousand) and is isotopically similar to a refinery sample (delta13C = -30.1 per thousand). No significant isotopic variation in dissolved MTBE implies that either the magnitude of any biodegradation-induced isotopic fractionation is small, or that relatively little degradation has taken place in the presence of BTEX hydrocarbons. It is possible, however, that MTBE degradation occurs under aerobic conditions in the absence of BTEX since no groundwater samples were taken with co-existing MTBE and oxygen. Low benzene delta13C values are correlated with high sulphate delta34S, indicating that little benzene degradation has occurred in the sulphate reduction zone. Benzene degradation may be associated with denitrification since increased benzene delta13C is associated with increased delta15N in residual nitrate. Re-supply of electron acceptors by diffusion from the matrix into fractures and dispersive mixing is an important constraint on degradation rates and natural attenuation capacity in this dual-porosity aquifer.     

A PCE groundwater plume discharging to a river: influence of the streambed and near-river zone on contaminant distributions

An investigation of a tetrachloroethene (PCE) groundwater plume originating at a dry cleaning facility on a sand aquifer and discharging to a river showed that the near-river zone strongly modified the distribution, concentration, and composition of the plume prior to discharging into the surface water. The plume, streambed concentration, and hydrogeology were extensively characterized using the Waterloo profiler, mini-profiler, conventional and driveable multilevel samplers (MLS), Ground Penetrating Radar (GPR) surveys, streambed temperature mapping (to identify discharge zones), drivepoint piezometers, and soil coring and testing. The plume observed in the shallow streambed deposits was significantly different from what would have been predicted based on the characteristics of the upgradient plume. Spatial and temporal variations in the plume entering the near-river zone contributed to the complex contaminant distribution observed in the streambed where concentrations varied by factors of 100 to 5000 over lateral distances of less than 1 to 3.5 m. Low hydraulic conductivity semi-confining deposits and geological heterogeneities at depth below the streambed controlled the pattern of groundwater discharge through the streambed and influenced where the plume discharged into the river (even causing the plume to spread out over the full width of the streambed at some locations). The most important effect of the near-river zone on the plume was the extensive anaerobic biodegradation that occurred in the top 2.5 m of the streambed, even though essentially no biodegradation of the PCE plume was observed in the upgradient aquifer. Approximately 54% of the area of the plume in the streambed consisted solely of PCE transformation products, primarily cis-1,2-dichloroethene (cDCE) and vinyl chloride (VC). High concentrations in the interstitial water of the streambed did not correspond to high groundwater-discharge zones, but instead occurred in low discharge zones and are likely sorbed or retarded remnants of past high-concentration plume discharges. The high-concentration areas (up to 5529 microg/l of total volatile organics) in the streambed are of ecological concern and represent potential adverse exposure locations for benthic and hyporheic zone aquatic life, but the effect of these exposures on the overall health of the river has yet to be determined. Even if the upgradient source of PCE is remediated and additional PCE is prevented from reaching the streambed, the high-concentration deposits in the streambed will likely take decades to hundreds of years to flush completely clean under natural conditions because these areas have low vertical groundwater flow velocities and high retardation factors. Despite high concentrations of contaminants in the streambed, PCE was detected in the surface water only rarely due to rapid dilution in the river and no cDCE or VC was detected. Neither the sampling of surface water nor the sampling of the groundwater from the aquifer immediately adjacent to the river gave an accurate indication of the high concentrations of PCE biodegradation products present in the streambed. Sampling of the interstitial water of the shallow streambed deposits is necessary to accurately characterize the nature of plumes discharging to rivers.     

Assessing impacts of partial mass depletion in DNAPL source zones: II. Coupling source strength functions to plume evolution

Analytical solutions, describing the time-dependent DNAPL source-zone mass and contaminant discharge rate, derived previously in Part I [Falta, R.W., Rao, P.S., Basu, N., this issue. Assessing the impacts of partial mass depletion in DNAPL source zones: I. Analytical modeling of source strength functions and plume response. J. Contam. Hydrol.] are used as a flux-boundary condition in a semi-analytical contaminant transport model. These analytical solutions assume a power relationship between the flow-averaged source concentration, and the source DNAPL mass; the empirical exponent (gamma) is a function of the flow field heterogeneity, DNAPL architecture, and the correlation between them. The DNAPL source strength terms can account for partial source remediation, either at time zero, or at some later time after the DNAPL release. The transport model considers advection, retardation, three-dimensional dispersion, and sequential first-order decay/production of several species. A separate solution is used to compute the time-dependent mass of each contaminant in the plume. A series of examples using different values of gamma shows how the benefits of partial DNAPL source remediation can vary with site conditions. In general, when gamma>1, relatively large short-term reductions in the plume concentrations and mass occur, but the source longevity is not strongly affected. Conversely, when gamma<1, the short-term reductions in the plume concentrations and mass are smaller, but the source longevity can be greatly reduced. In either case, the source remediation effort is much more effective if it is undertaken at an early time, before much contaminant mass has entered the plume. If the remediation effort is significantly delayed, the leading parts of the plume are not affected by the source remediation, and additional control or remediation of the plume itself is required.