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.     

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.     

Biogeochemical characterisation of a coal tar distillate plume

The distillation of acidified coal tars for up to 50 years has given rise to a phenol plume approximately 500 m long, 50 m deep and containing up to 15 g l(-1) dissolved organic carbon (DOC) in the Triassic Sandstones aquifer. A conceptual biogeochemical model based on chemical and microbiological analysis of groundwater samples has been developed as a preliminary to more detailed studies of the controls on natural attenuation. While the development of redox zones and the production of methane and carbon dioxide provide evidence of natural attenuation, it appears that degradation is slow. The existence of sulphate in the plume indicates that this electron acceptor has not been depleted and that consequently methanogenesis is probably limited. Based on a simple estimate of sulphate input concentration, a half-life of about 15 years has been estimated for sulphate reduction. Geochemical modelling predicts that increased alkalinity within the plume has not led to carbonate precipitation, and thus within the limits of accuracy of the measurement, alkalinity may reflect the degree of biodegradation. This implies a loss of around 18% of the DOC over a 30-year period. Despite limited degradation, microbial studies show that there are diverse microbial communities in the aquifer with the potential for both anaerobic and aerobic biodegradation. Microbial activity was found to be greatest at the leading edge of the plume where DOC concentrations are 60 mg l(-1) or less, but activity could still be observed in more contaminated samples even though cells could not be cultured. The study suggests that degradation may be limited by the high phenol concentrations within the core of the plume, but that once diluted by dispersion, natural attenuation may proceed. More detailed studies to confirm these initial findings are identified and form the basis of associated papers.     

Variations in expression of carbon isotope fractionation of chlorinated ethenes during biologically enhanced PCE dissolution close to a source zone

The stable carbon isotope values of tetrachloroethene (PCE) and its degradation products were monitored during studies of biologically enhanced dissolution of PCE dense nonaqueous phase liquid (DNAPL) to determine the effect of PCE dissolution on observed isotope values. The degradation of PCE was monitored in a 2-dimensional model aquifer and in a pilot test cell (PTC) at Dover Air Force Base, both with emplaced PCE DNAPL sources. Within the plume down gradient from the source, the isotopic fractionation of dissolved PCE and its degradation products were consistent with those observed in biodegradation laboratory studies. However, close to the source zone significant shifts in the isotope values of dissolved PCE were not observed in either the model aquifer or PTC due to the constant input of newly dissolved, non fractionated PCE, and the small isotopic fractionation associated with PCE reductive dechlorination by the mixed microbial culture used. Therefore the identification of reductive dechlorination in the presence of PCE DNAPL was based upon the appearance of daughter products and the isotope values of those daughter products. An isotope model was developed to simulate isotope values of PCE during the dissolution and degradation of PCE adjacent to a DNAPL source zone. With the exception of very high degradation rate constants (>1/day) stable carbon isotope values of PCE estimated by the model remained within error of the isotope value of the PCE DNAPL, consistent with measured isotope values in the model aquifer and in the PTC.     

In situ modification of an aquifer material by a cationic surfactant to enhance retardation of organic contaminants

Sorption of hexadecyltrimethylammonium chloride (HDTMA), a cationic surfactant, on aquifer material from Columbus AFB, Mississippi, U.S.A., was examined. Transport studies using flow-through columns and a box model aquifer showed that an almost stationary zone of HDTMA-modified aquifer material could be produced in situ without a significant decrease in hydraulic conductivity.Perchloroethylene (PCE) and naphthalene sorption isotherms on the HDTMA-modified aquifer material were linear, and sorption coefficients were increased by over two orders of magnitude relative to the unmodified material. The retardation of PCE by insitu emplaced HDTMA zones within a column was examined. Agreement between batch- and column-derived sorption coefficients and breakthrough curve symmetry indicates that local equilibrium was attained. Significant retardation of a naphthalene plume by an in situ emplaced surfactant zone was demonstrated in the box model aquifer system.The experimental results indicate that it is feasible to create in situ a sorbent zone within an aquifer using cationic surfactants. In most situations, the sorbent zone concept needs to be coupled with contaminant degradation processes for sorbent emplacement to be a practical tool in the remediation of groundwater contamination sites. Sorbent zones may be of benefit in the engineering of suitable environments for microbial or abiotic degradation reactions and by providing time slow reactions to occur.     

Multi-component reactive transport modeling of natural attenuation of an acid groundwater plume at a uranium mill tailings site

Natural attenuation of an acidic plume in the aquifer underneath a uranium mill tailings pond in Wyoming, USA was simulated using the multi-component reactive transport code PHREEQC. A one-dimensional model was constructed for the site and the model included advective-dispersive transport, aqueous speciation of 11 components, and precipitation-dissolution of six minerals. Transport simulation was performed for a reclamation scenario in which the source of acidic seepage will be terminated after 5 years and the plume will then be flushed by uncontaminated upgradient groundwater. Simulations show that successive pH buffer reactions with calcite, Al(OH)3(a), and Fe(OH)3(a) create distinct geochemical zones and most reactions occur at the boundaries of geochemical zones. The complex interplay of physical transport processes and chemical reactions produce multiple concentration waves. For SO4(2-) transport, the concentration waves are related to advection-dispersion, and gypsum precipitation and dissolution. Wave speeds from numerical simulations compare well to an analytical solution for wave propagation.     

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.     

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.     

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.     

Mixing-controlled biodegradation in a toluene plume--results from two-dimensional laboratory experiments

Various abiotic and biotic processes such as sorption, dilution, and degradation are known to affect the fate of organic contaminants, such as petroleum hydrocarbons in saturated porous media. Reactive transport modeling of such plumes indicates that the biodegradation of organic pollutants is, in many cases, controlled by mixing and therefore occurs locally at the plume's fringes, where electron donors and electron-acceptors mix. Herein, we aim to test whether this hypothesis can be verified by experimental results obtained from aerobic and anaerobic degradation experiments in two-dimensional sediment microcosms. Toluene was selected as a model compound for oxidizable contaminants. The two-dimensional microcosm was filled with quartz sand and operated under controlled flow conditions simulating a contaminant plume in otherwise uncontaminated groundwater. Aerobic degradation of toluene by Pseudomonas putida mt-2 reduced a continuous 8.7 mg L(-1) toluene concentration by 35% over a transport distance of 78 cm in 15.5 h. In comparison, under similar conditions Aromatoleum aromaticum strain EbN1 degraded 98% of the toluene infiltrated using nitrate (68.5+/-6.2 mg L(-1)) as electron acceptor. A major part of the biodegradation activity was located at the plume fringes and the slope of the electron-acceptor gradient was steeper during periods of active biodegradation. The distribution of toluene and the significant overlap of nitrate at the plume's fringe indicate that biokinetic and/or microscale transport processes may constitute additional limiting factors. Experimental data is corroborated with results from a reactive transport model using double Monod kinetics. The outcome of the study shows that in order to simulate degradation in contaminant plumes, detailed data sets are required to test the applicability of models. These will have to deal with the incorporation of existing parameters coding for substrate conversion kinetics and microbial growth.     

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.     

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.     

Microcosm studies of microbial degradation in a coal tar distillate plume

Investigation of a groundwater plume containing up to 24 g l(-1) phenolic compounds suggested that over a period of nearly 50 years, little degradation had occurred despite the presence of a microbial community and electron acceptors within the core of the plume. In order to study the effect of contaminant concentration on degradation behaviour, laboratory microcosm experiments were performed under aerobic and anaerobic conditions at four different concentrations obtained by diluting contaminated with uncontaminated groundwater. The microcosms contained groundwater with total phenols at ca. 200, 250, 660 and 5000 mg l(-1), and aquifer sediment that had been acclimatised within the plume for several months. The microcosms were operated for a period of 390-400 days along with sterile controls to ascertain whether degradation was microbially mediated or abiotic. Under aerobic conditions, degradation only occurred at concentrations up to 660 mg l(-1) total phenols. At phenol concentrations below 250 mg l(-1) a benzoquinone intermediate, thought to originate from the degradation of 2,5-dimethylphenol, was isolated and identified. This suggested an unusual degradative pathway for this compound; its aerobic degradation more commonly proceeding via catecholic intermediates. Under anaerobic conditions, degradation only occurred in the most dilute microcosm (total phenols 195 mg l(-1)) with a loss of p-cresol accompanied by a nonstoichiometric decrease in nitrate and sulphate. By inference, iron(III) from the sediment may also have been used as a terminal electron acceptor, in which case the amount of biologically available iron released was calculated as 1.07 mg Fe(III)/g of sediment. The study shows that natural attenuation is likely to be stimulated by dilution of the plume.     

Denitrification and phenol degradation in a contaminated aquifer

A natural groundwater system modified by pollutant phenols and agricultural nitrate has been modelled in the laboratory by a series of sacrificial microcosm experiments. Samples of aquifer sediment and groundwater from the margin of the phenol plume were used to inoculate anaerobic microcosms enriched in nitrate and pollutant phenols. Rapid degradation of phenol and p-cresol was observed over a 35-day period leading to the generation of inorganic carbon and a number of transient intermediates. O-cresol proved to be recalcitrant on the experimental time-scale. A mass balance calculation shows that, during degradation, carbon was conserved in the aqueous phase. Groundwater-sediment interactions were monitored using carbon stable isotope data. A mass balance for solution TIC indicates thatp-cresol degradation stimulated the dissolution of sedimentary carbonate phases due to the formation of carbonic acid. Compound-specific carbon isotope analysis (GC-IRMS) was used to search for 13C enrichment in residual p-cresol. A slight enrichment trend (epsilon = -2.5/1000) was tentatively identified. The potential of this fractionation effect for obtaining in situ degradation rates is discussed. Results from the microcosm experiments help to explain the observed distribution of nitrate and phenols within the polluted aquifer.     

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.     

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.     

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.     

Identification of a reactive degradation zone at a landfill leachate plume fringe using high resolution sampling and incubation techniques

Vertical small-scale variation in phenoxy acid herbicide degradation across a landfill leachate plume fringe was studied using laboratory degradation experiments. Sediment cores (subdivided into 5 cm segments) were collected in the aquifer and the sediment and porewater were used for microcosm experiments (50 experiments) and for determination of solid organic carbon, solid-water partitioning coefficients, specific phenoxy acid degraders and porewater chemistry. Results from a multi-level sampler installed next to the cores provided information on the plume position and oxygen concentration in the groundwater. Oxygen concentration was controlled individually in each microcosm to mimic the conditions at their corresponding depths. A highly increased degradation potential existed at the narrow plume fringe (37.7 to 38.6 masl), governed by the presence of phenoxy acids and oxygen. This resulted in the proliferation of a microbial population of specific phenoxy acid degraders, which further enhanced the degradation potential for phenoxy acids at the fringe. The results illustrate the importance of fringe degradation processes in contaminant plumes. Furthermore, they highlight the relevance of using high-resolution sampling techniques as well as controlled microcosm experiments in the assessment of the natural attenuation capacity of contaminant plumes in groundwater.     

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.     

An analytical quantification of mass fluxes and natural attenuation rate constants at a former gasworks site

A new integral groundwater investigation approach was used for the first time to quantify natural attenuation rates at field scale. In this approach, pumping wells positioned along two control planes were operated at distances of 140 and 280 m downstream of a contaminant source zone at a former gasworks site polluted with BTEX- (benzene, toluene, ethyl-benzene, o-, p-xylene) and PAH- (polycyclic aromatic hydrocarbons) compounds. Based on the quantified changes in total contaminant mass fluxes between the control planes, first-order natural attenuation rate constants could be estimated. For BTEX-compounds, these ranged from 1.4e-02 to 1.3e-01 day(-1) whereas for PAH-compounds natural attenuation rate constants of 3.7e-04 to 3.1e-02 day(-1) were observed. Microbial degradation activity at the site was indicated by an increase in dissolved iron mass flux and a reduction in sulphate mass flux between the two investigated control planes. In addition to information about total contaminant mass fluxes and average concentrations, an analysis of the concentration-time series measured at the control planes also allowed to semi-quantitatively delineate the aquifer regions most likely contaminated by the BTEX- and PAH-compounds.