Assessing the transformation of chlorinated ethenes in aquifers with limited potential for natural attenuation: added values of compound-specific carbon isotope analysis and groundwater dating

The evaluation of biotransformation of chlorinated ethenes (CEs) in contaminated aquifers is challenging when variable redox conditions and groundwater flow regime are limiting factors. By using compound-specific stable carbon isotope analysis (C-CSIA) and ³H-³He based groundwater dating, we assessed three CE-contaminated field sites that differed in groundwater flow velocities, redox conditions, and level of contamination. CE isotopic signatures and carbon isotopic mass balances were applied to quantify CE transformation, whereas groundwater dating allowed determining degradation timescales and assessing hydrodynamic regimes. The combination of these techniques enabled at all field sites to indicate zones within the aquifers where CE dechlorination preferably occurred, sometimes even to metabolites of no toxic concern. However, the natural transformation processes were insufficient to mitigate the entire CE contamination at the studied sites. Such situations of limited transformation are worldwide far more common than sites where optimal natural (mainly redox) conditions are enabling complete CEs degradation. Despite such constraints for natural transformation, this study showed that even under non-favorable biogeochemical CEs degradation, the combination of CSIA and groundwater dating provide valuable information to the understanding of the fate of the CEs, thus, being an important contribution in the definition of efficient remediation measures at any given biogeochemical conditions.     

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.     

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

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

Occurrence and attenuation of specific organic compounds in the groundwater plume at a former gasworks site

The changing contaminant pattern with travelled distance was investigated in the anaerobic groundwater plume downstream from an extended zone containing residual NAPL at a former gas manufacturing plant. With increasing distance, O- and N-heterocyclic aromatic compounds are enriched in the plume relative to the usually assessed coal tar constituents (poly- and monocyclic aromatic compounds). In a first approximation, the overall concentration decrease of the investigated compounds follows a first order overall decay. The half life distance in the plume downgradient from the source varied between 20 m for benzene and up to 167-303 m for alkyl-naphthalenes. Acenaphthene is degraded only within about 50 m downstream from the source area, then its concentration remains constant (ca. 180 microg/l) and far above the legal limit. Dimethyl-benzofurans were the most recalcitrant among all compounds which could be quantified with the analytical method available. The overall groundwater contamination in the plume is seriously underestimated if only BTEX and 16-EPA-PAHs are monitored.     

Attenuation of landfill leachate by UK Triassic sandstone aquifer materials: 2. Sorption and degradation of organic pollutants in laboratory columns

The sorption and degradation of dissolved organic matter (DOM) and 13 organic micropollutants (BTEX, aromatic hydrocarbons, chloro-aromatic and -aliphatic compounds, and pesticides) in acetogenic and methanogenic landfill leachate was studied in laboratory columns containing Triassic sandstone aquifer materials from the English Midlands. Solute sorption and degradation relationships were evaluated using a simple transport model. Relative to predictions, micropollutant sorption was decreased up to eightfold in acetogenic leachate, but increased up to sixfold in methanogenic leachate. This behaviour reflects a combination of interactions between the micropollutants, leachate DOM and aquifer mineral fraction. Sorption of DOM was not significant. Degradation of organic fractions occurred under Mn-reducing and SO₄-reducing conditions. Degradation of some micropollutants occurred exclusively under Mn-reducing conditions. DOM and benzene were not significantly degraded under the conditions and time span (up to 280 days) of the experiments. Most micropollutants were degraded immediately or after a lag phase (32–115 days). Micropollutant degradation rates varied considerably (half-lives of 8 to >2000 days) for the same compounds (e.g., TeCE) in different experiments, and for compounds (e.g., naphthalene, DCB and TeCA) within the same experiment. Degradation of many micropollutants was both simultaneous and sequential, and inhibited by the utilisation of different substrates. This mechanism, in combination with lag phases, controls micropollutant degradation potential in these systems more than the degradation rate. These aquifer materials have a potentially large capacity for in situ bioremediation of organic pollutants in landfill leachate and significant degradation may occur in the Mn-reducing zones of leachate plumes. However, degradation of organic pollutants in acetogenic leachate may be limited in aquifers with low pH buffering capacity and reducible Mn oxides. Contaminants in this leachate present a greater risk to groundwater resources in these aquifers than methanogenic leachate.     

Joint propagation of variability and imprecision in assessing the risk of groundwater contamination

Estimating risks of groundwater contamination often require schemes for representing and propagating uncertainties relative to model input parameters. The most popular method is the Monte Carlo method whereby cumulative probability distributions are randomly sampled in an iterative fashion. The shortcoming of the approach, however, arises when probability distributions are arbitrarily selected in situations where available information is incomplete or imprecise. In such situations, alternative modes of information representation can be used, for example the nested intervals known as “possibility distributions”. In practical situations of groundwater risk assessment, it is common that certain model parameters may be represented by single probability distributions (representing variability) because there are data to justify these distributions, while others are more faithfully represented by possibility distributions (representing imprecision) due to the partial nature of available information. This paper applies two recent methods, designed for the joint-propagation of variability and imprecision, to a groundwater contamination risk assessment. Results of the joint-propagation methods are compared to those obtained using both interval analysis and the Monte Carlo method with a hypothesis of stochastic independence between model parameters. The two joint-propagation methods provide results in the form of families of cumulative distributions of the probability of exceeding a certain value of groundwater concentration. These families are delimited by an upper cumulative distribution and a lower distribution respectively called Plausibility and Belief after evidence theory. Slight differences between the results of the two joint-propagation methods are explained by the different assumptions regarding parameter dependencies. Results highlight the point that non-conservative results may be obtained if single cumulative probability distributions are arbitrarily selected for model parameters in the face of imprecise information and the Monte Carlo method is used under the assumption of stochastic independence. The proposed joint-propagation methods provide upper and lower bounds for the probability of exceeding a tolerance threshold. As this may seem impractical in a risk-management context, it is proposed to introduce “a-posteriori subjectivity” (as opposed to the “a-priori subjectivity” introduced by the arbitrary selection of single probability distributions) by defining a single indicator of evidence as a weighted average of Plausibility and Belief, with weights to be defined according to the specific context.     

Comment on “Analysis of groundwater contamination using concentration-time series recorded during an integral pumping test: Bias introduced by strong concentration gradients within the plume” by Allelign Zeru and Gerhard Schäfer

We consider the results of a recent paper in this journal [Zeru, A. and Schäfer, G., 2005. Analysis of groundwater contamination using concentration–time series recorded during an integral pumping test: Bias introduced by strong concentration gradients within the plume. Journal of Contaminant Hydrology 81 (2005) 106–124], which addresses the field-scale characterisation of contaminant plumes in groundwater. There, it is concluded that contaminant concentration gradients can bias Integral Pumping Test (IPT) interpretations considerably, in particular if IPTs are conducted in advective fronts of contaminant plumes. We discuss implications of this setting and also argue that the longitudinal and transverse dispersivities used in the examples of Zeru and Schäfer (2005) of up to 30 m and 3 m, respectively, are generally very high for the here relevant capture zone scale (< 20 m). However, regardless of both longitudinal and transverse concentration gradients, we further show through a counter-example that IPT results are unbiased as long as the concentration attenuation along the flow direction is linear over the capture zone extent.     

Biogeochemical processes at the fringe of a landfill leachate pollution plume: potential for dissolved organic carbon, Fe(II), Mn(II), NH4, and CH4 oxidation

Various redox reactions may occur at the fringe of a landfill leachate plume, involving oxidation of dissolved organic carbon (DOC), CH4, Fe(II), Mn(II), and NH4 from leachate and reduction of O2, NO3 and SO4 from pristine groundwater. Knowledge on the relevance of these processes is essential for the simulation and evaluation of natural attenuation (NA) of pollution plumes. The occurrence of such biogeochemical processes was investigated at the top fringe of a landfill leachate plume (Banisveld, the Netherlands). Hydrochemical depth profiles of the top fringe were captured via installation of a series of multi-level samplers at 18, 39 and 58 m downstream from the landfill. Ten-centimeter vertical resolution was necessary to study NA within a fringe as thin as 0.5 m. Bromide appeared an equally well-conservative tracer as chloride to calculate dilution of landfill leachate, and its ratio to chloride was high compared to other possible sources of salt in groundwater. The plume fringe rose steadily from a depth of around 5 m towards the surface with a few meters in the period 1998-2003. The plume uplift may be caused by enhanced exfiltration to a brook downstream from the landfill, due to increased precipitation over this period and an artificial lowering of the water level of the brook. This rise invoked cation exchange including proton buffering, and triggered degassing of methane. The hydrochemical depth profile was simulated in a 1D vertical reactive transport model using PHREEQC-2. Optimization using the nonlinear optimization program PEST brought forward that solid organic carbon and not clay minerals controlled retardation of cations. Cation exchange resulted in spatial separation of Fe(II), Mn(II) and NH4 fronts from the fringe, and thereby prevented possible oxidation of these secondary redox species. Degradation of DOC may happen in the fringe zone. Re-dissolution of methane escaped from the plume and subsequent oxidation is an explanation for absence of previously present nitrate and anaerobic conditions in pristine groundwater above the plume. Stable carbon isotope (delta13C) values of methane confirm anaerobic methane oxidation immediately below the fringe zone, presumably coupled to reduction of sulfate, desorbed from iron oxide. Methane must be the principle reductant consuming soluble electron-acceptors in pristine groundwater, thereby limiting NA for other solutes including organic micro-pollutants at the fringe of this landfill leachate plume.     

Powdered activated carbon coupled with enhanced coagulation for natural organic matter removal and disinfection by-product control: application in a Western Australian water treatment plant

The removal of organic precursors of disinfection by-products (DBPs), i.e. natural organic matter (NOM), prior to disinfection and distribution is considered as the most effective approach to minimise the formation of DBPs. This study investigated the impact of the addition of powdered activated carbon (PAC) to an enhanced coagulation treatment process at an existing water treatment plant on the efficiency of NOM removal, the disinfection behaviour of the treated water, and the water quality in the distribution system. This is the first comprehensive assessment of the efficacy of plant-scale application of PAC combined with enhanced coagulation on an Australian source water. As a result of the PAC addition, the removal of NOM improved by 70%, which led to a significant reduction (80-95%) in the formation of DBPs. The water quality in the distribution system also improved, indicated by lower concentrations of DBPs in the distribution system and better maintenance of disinfectant residual at the extremities of the distribution system. The efficacy of the PAC treatment for NOM removal was shown to be a function of the characteristics of the NOM and the quality of the source water, as well as the PAC dose. PAC treatment did not have the capacity to remove bromide ion, resulting in the formation of more brominated DBPs. Since brominated DBPs have been found to be more toxic than their chlorinated analogues, their preferential formation upon PAC addition must be considered, especially in source waters containing high concentrations of bromide.     

Transport and degradation of propylene glycol in the vadose zone: model development and sensitivity analysis

Transport and degradation of de-icing chemical (containing propylene glycol, PG) in the vadose zone were studied with a lysimeter experiment and a model, in which transient water flow, kinetic degradation of PG and soil chemistry were combined. The lysimeter experiment indicated that aerobic as well as anaerobic degradation occurs in the vadose zone. Therefore, the model included both types of degradation, which was made possible by assuming advection-controlled (mobile) and diffusion-controlled (immobile) zones. In the mobile zone, oxygen can be transported by diffusion in the gas phase. The immobile zone is always water-saturated, and oxygen only diffuses slowly in the water phase. Therefore, the model is designed in a way that the redox potential can decrease when PG is degraded, and thus, anaerobic degradation can occur. In our model, manganese oxide (MnO₂, which is present in the soil) and NO \(_{3}^{-}\) (applied to enhance biodegradation) can be used as electron acceptors for anaerobic degradation. The application of NO \(_{3}^{-}\) does not result in a lower leaching of PG nor in a slower depletion of MnO₂. The thickness of the snowcover influences the leached fraction of PG, as with a high infiltration rate, transport is fast, there is less time for degradation and thus more PG will leach. The model showed that, in this soil, the effect of the water flow dominates over the effect of the degradation parameters on the leaching at a 1-m depth.     

Bottom-up approach for the reaction of xenobiotics and their metabolites with model substances for natural organic matter by electrochemistry-mass spectrometry (EC-MS)

Risk assessment of xenobiotics requires a comprehensive understanding of their transformation in the environment. As most of the transformation processes usually involve a redox reaction or a hydrolysis as the first steps of the transformation, we applied an approach that uses an electrochemical cell to investigate model “redox” reactions in aqueous solutions for environmental processes. We investigated the degradation of a variety of xenobiotics from polar to nonpolar and analyzed their degradation products by on-line coupling of electrochemistry with mass spectrometry (EC-MS). Furthermore, we evaluated possible binding reactions with regard to the generation of non-extractable residues with some model substances (catechol, phthalic acid, γ-l-Glutamyl-l-cysteinyl-glycine (GSH) and l-histidine) deduced from a natural organic matter (NOM) structure model and identified possible binding-sites.Whereas typically investigations in soil/water-systems have been applied, we used to our knowledge for the first time a bottom-up approach, starting from the chemicals of interest and different model substances for natural organic matter to evaluate chemical binding mechanisms (or processes) in the EC-MS under redox conditions. Under oxidative conditions, bindings of the xenobiotics with catechol, GSH and histidine were found, but no reactions with the model compound phthalic acid were observed. In general, no chemical binding has yet been found under reductive conditions. In some cases (i.e. benzo[a]anthracene) the oxidation product only underwent a binding reaction, whereas the xenobiotic itself did not undergo any reactions.EC-MS is a promising fast and simple screening method to investigate the environmental behavior of xenobiotics and to evaluate the potential risks of newly synthesized substances.     

An inventory of nitrous oxide emissions from agriculture in the UK using the IPCC methodology: emission estimate,uncertainty and sensitivity analysis

Nitrous oxide emission from UK agriculture was estimated, using the IPCC default values of all emission factors and parameters, to be 87 Gg N₂O–N in both 1990 and 1995. This estimate was shown, however, to have an overall uncertainty of 62%. The largest component of the emission (54%) was from the direct (soil) sector. Two of the three emission factors applied within the soil sector, EF1 (direct emission from soil) and EF3_PRP (emission from pasture range and paddock) were amongst the most influential on the total estimate, producing a ±31 and +11% to −17% change in emissions, respectively, when varied through the IPCC range from the default value. The indirect sector (from leached N and deposited ammonia) contributed 29% of the total emission, and had the largest uncertainty (126%). The factors determining the fraction of N leached (Frac_LEACH) and emissions from it (EF5), were the two most influential. These parameters are poorly specified and there is great potential to improve the emission estimate for this component. Use of mathematical models (NCYCLE and SUNDIAL) to predict Frac_LEACH suggested that the IPCC default value for this parameter may be too high for most situations in the UK.Comparison with other UK-derived inventories suggests that the IPCC methodology may overestimate emission. Although the IPCC approach includes additional components to the other inventories (most notably emission from indirect sources), estimates for the common components (i.e. fertiliser and animals), and emission factors used, are higher than those of other inventories. Whilst it is recognised that the IPCC approach is generalised in order to allow widespread applicability, sufficient data are available to specify at least two of the most influential parameters, i.e. EF1 and Frac_LEACH, more accurately, and so provide an improved estimate of nitrous oxide emissions from UK agriculture.     

Impact analysis of carbon taxes on the steamcoal and natural gas market in the eastern United States: a linear complementarity model

The impact of a federal excise tax on Appalachian steamcoal on the consumption and production of steamcoal and natural gas is investigated. To examine the comprehensive effects of such taxes, the linear complementarity model is used to ensure consistent results. Predictions on coal and natural gas prices are also made for useful policy implementations.     

Modeling an air pollution episode in northwestern United States: Identifying the effect of nitrogen oxide and volatile organic compound emission changes on air pollutants formation using direct sensitivity analysis

Air quality impacts of volatile organic compound (VOC) and nitrogen oxide (NO_x) emissions from major sources over the northwestern United States are simulated. The comprehensive nested modeling system comprises three models: Community Multiscale Air Quality (CMAQ), Weather Research and Forecasting (WRF), and Sparse Matrix Operator Kernel Emissions (SMOKE). In addition, the decoupled direct method in three dimensions (DDM-3D) is used to determine the sensitivities of pollutant concentrations to changes in precursor emissions during a severe smog episode in July of 2006. The average simulated 8-hr daily maximum O₃ concentration is 48.9 ppb, with 1-hr O₃ maxima up to 106 ppb (40 km southeast of Seattle). The average simulated PM_2.5 (particulate matter with an aerodynamic diameter <2.5 μm) concentration at the measurement sites is 9.06 μg m^?3, which is in good agreement with the observed concentration (8.06 μg m^?3). In urban areas (i.e., Seattle, Vancouver, etc.), the model predicts that, on average, a reduction of NO_x emissions is simulated to lead to an increase in average 8-hr daily maximum O₃ concentrations, and will be most prominent in Seattle (where the greatest sensitivity is??0.2 ppb per % change of mobile sources). On the other hand, decreasing NO_x emissions is simulated to decrease the 8-hr maximum O₃ concentrations in remote and forested areas. Decreased NO_x emissions are simulated to slightly increase PM_2.5 in major urban areas. In urban areas, a decrease in VOC emissions will result in a decrease of 8-hr maximum O₃ concentrations. The impact of decreased VOC emissions from biogenic, mobile, nonroad, and area sources on average 8-hr daily maximum O₃ concentrations is up to 0.05 ppb decrease per % of emission change, each. Decreased emissions of VOCs decrease average PM_2.5 concentrations in the entire modeling domain. In major cities, PM_2.5 concentrations are more sensitive to emissions of VOCs from biogenic sources than other sources of VOCs. These results can be used to interpret the effectiveness of VOC or NO_x controls over pollutant concentrations, especially for localities that may exceed National Ambient Air Quality Standards (NAAQS).

Implications: The effect of NO_x and VOC controls on ozone and PM_2.5 concentrations in the northwestern United States is examined using the decoupled direct method in three dimensions (DDM-3D) in a state-of-the-art three-dimensional chemical transport model (CMAQ). NO_x controls are predicted to increase PM_2.5 and ozone in major urban areas and decrease ozone in more remote and forested areas. VOC reductions are helpful in reducing ozone and PM_2.5 concentrations in urban areas. Biogenic VOC sources have the largest impact on O₃ and PM_2.5 concentrations.     

Determination of chlorinated paraffins in sediments from the Firth of Clyde by gas chromatography with electron capture negative ionisation mass spectrometry and carbon skeleton analysis by gas chromatography with flame ionisation detection

Short chain chlorinated paraffins (SCCPs) are a group of persistent organic pollutants (POPs) of increasing concern, but are to date not widely investigated in the environment, largely due to the challenges involved in their quantification. Here, SCCPs were quantified in marine sediments from the Firth of Clyde, Scotland, by gas chromatography with electron capture negative ionisation mass spectrometry (GC-ECNIMS) and through carbon skeleton analysis by gas chromatography with flame ionisation detection (GC-FID), and the analytical challenges encountered are discussed. Concentrations in the sediments ranged from 0.4 to 69 μg kg⁻¹ when determined by GC-ECNIMS, and from 5.6 to 379 μg kg⁻¹ when determined by GC-FID. For 8 out of 11 samples, analysis by GC-FID gave higher results than analysis by GC-ECNIMS. Unexpected aspects of the analysis, such as the presence of high concentrations of longer chain chlorinated paraffins in the samples, are also presented.