Evaluating alternate biokinetic models for trace pollutant cometabolism.

Mathematical models of cometabolic biodegradation kinetics can improve our understanding of the relevant microbial reactions and allow us to design in situ or in-reactor applications of cometabolic bioremediation. A variety of models are available, but their ability to describe experimental data has not been systematically evaluated for a variety of operational/experimental conditions. Here five different models were considered: first-order; Michaelis-Menten; reductant; competition; and combined models. The models were assessed on their ability to fit data from simulated batch experiments covering a realistic range of experimental conditions. The simulated observations were generated by using the most complex model structure and parameters based on the literature, with added experimental error. Three criteria were used to evaluate model fit: ability to fit the simulated experimental data, identifiability of parameters using a colinearity analysis, and suitability of the model size and complexity using the Bayesian and Akaike Information criteria. Results show that no single model fits data well for a range of experimental conditions. The reductant model achieved best results, but required very different parameter sets to simulate each experiment. Parameter nonuniqueness was likely to be due to the parameter correlation. These results suggest that the cometabolic models must be further developed if they are to reliably simulate experimental and operational data.

Review of reactive kinetic models describing reductive dechlorination of chlorinated ethenes in soil and groundwater.

Reductive dechlorination is a major degradation pathway of chlorinated ethenes in anaerobic subsurface environments, and reactive kinetic models describing the degradation process are needed in fate and transport models of these contaminants. However, reductive dechlorination is a complex biological process, where many microbial populations including dechlorinating, fermentative, methanogenic, iron and sulfate reducing, interact. In this article the modeling approaches and the experimental data needed to calibrate them are reviewed, classified, and discussed. Model approaches considered include first order kinetics, Monod kinetics to describe sequential reductive dechlorination and bacterial growth, and metabolic models which simulate fermentation and redox processes interacting with reductive dechlorination processes. The review shows that the estimated kinetic parameters reported vary over a wide range, and that experimental microbial data are scarce. Very few studies have been performed evaluating the influence of sulfate and iron reduction, and contradictory conclusions on the interaction of redox processes with reductive dechlorination have been reported. The modeling approaches for metabolic reductive dechlorination employing different descriptions of the interaction between redox and dechlorination processes and competition for hydrogen are classified. The current concepts lead to different results, suggesting a need for further investigations on the interactions between the microbial communities performing dechlorination and redox processes, including the establishment of biomarkers quantifying dechlorination, and on geochemical characterization. Finally, the relevance of laboratory data and the development of practical modeling tools for field applications are discussed.

The gompertz function can coherently describe microbial mineralization of growth-sustaining pesticides.

Mineralization of (14)C-labeled tracers is a common way of studying the environmental fate of xenobiotics, but it can be difficult to extract relevant kinetic parameters from such experiments since complex kinetic functions or several kinetic functions may be needed to adequately describe large data sets. In this study, we suggest using a two-parameter, sigmoid Gompertz function for parametrizing mineralization curves. The function was applied to a data set of 252 normalized mineralization curves that represented the potential for degradation of the herbicide MCPA in three horizons of an agricultural soil. The Gompertz function fitted most of the normalized curves, and trends in the data set could be visualized by a scatter plot of the two Gompertz parameters (rate constant and time delay). For agricultural topsoil, we also tested the effect of the MCPA concentration on the mineralization kinetics. Reduced initial concentrations lead to shortened lag-phases, probably due to reduced need for bacterial growth. The effect of substrate concentration could be predicted by simply changing the time delay of the Gompertz curves. This delay could to some extent also simulate concentration effects for 2,4-D mineralization in agricultural soil and aquifer sediment and 2,6-dichlorobenzamide mineralization in single-species, mineral medium.

Is there an environmental benefit from remediation of a contaminated site? Combined assessments of the risk reduction and life cycle impact of remediation.

A comparative life cycle assessment is presented for four different management options for a trichloroethene-contaminated site with a contaminant source zone located in a fractured clay till. The compared options are (i) long-term monitoring (ii) in-situ enhanced reductive dechlorination (ERD), (iii) in-situ chemical oxidation (ISCO) with permanganate and (iv) long-term monitoring combined with treatment by activated carbon at the nearby waterworks. The life cycle assessment included evaluation of both primary and secondary environmental impacts. The primary impacts are the local human toxic impacts due to contaminant leaching into groundwater that is used for drinking water, whereas the secondary environmental impacts are related to remediation activities such as monitoring, drilling and construction of wells and use of remedial amendments. The primary impacts for the compared scenarios were determined by a numerical risk assessment and remedial performance model, which predicted the contaminant mass discharge over time at a point of compliance in the aquifer and at the waterworks. The combined assessment of risk reduction and life cycle impacts showed that all management options result in higher environmental impacts than they remediate, in terms of person equivalents and assuming equal weighting of all impacts. The ERD and long-term monitoring were the scenarios with the lowest secondary life cycle impacts and are therefore the preferred alternatives. However, if activated carbon treatment at the waterworks is required in the long-term monitoring scenario, then it becomes unfavorable because of large secondary impacts. ERD is favorable due to its low secondary impacts, but only if leaching of vinyl chloride to the groundwater aquifer can be avoided. Remediation with ISCO caused the highest secondary impacts and cannot be recommended for the site.

Development and sensitivity analysis of a fully kinetic model of sequential reductive dechlorination in groundwater.

A fully kinetic biogeochemical model of sequential reductive dechlorination (SERD) occurring in conjunction with lactate and propionate fermentation, iron reduction, sulfate reduction, and methanogenesis was developed. Production and consumption of molecular hydrogen (H(2)) by microorganisms have been modeled using modified Michaelis-Menten kinetics and has been implemented in the geochemical code PHREEQC. The model have been calibrated using a Shuffled Complex Evolution Metropolis algorithm to observations of chlorinated solvents, organic acids, and H(2) concentrations in laboratory batch experiments of complete trichloroethene (TCE) degradation in natural sediments. Global sensitivity analysis was performed using the Morris method and Sobol sensitivity indices to identify the most influential model parameters. Results show that the sulfate concentration and fermentation kinetics are the most important factors influencing SERD. The sensitivity analysis also suggests that it is not possible to simplify the model description if all system behaviors are to be well described.

Increasing urban water self-sufficiency: new era, new challenges.

Urban water supplies are traditionally based on limited freshwater resources located outside the cities. However, a range of concepts and techniques to exploit alternative water resources has gained ground as water demands begin to exceed the freshwater available to cities. Based on 113 cases and 15 in-depth case studies, solutions used to increase water self-sufficiency in urban areas are analyzed. The main drivers for increased self-sufficiency were identified to be direct and indirect lack of water, constrained infrastructure, high quality water demands and commercial and institutional pressures. Case studies demonstrate increases in self-sufficiency ratios to as much as 80% with contributions from recycled water, seawater desalination and rainwater collection. The introduction of alternative water resources raises several challenges: energy requirements vary by more than a factor of ten amongst the alternative techniques, wastewater reclamation can lead to the appearance of trace contaminants in drinking water, and changes to the drinking water system can meet tough resistance from the public. Public water-supply managers aim to achieve a high level of reliability and stability. We conclude that despite the challenges, self-sufficiency concepts in combination with conventional water resources are already helping to reach this goal.

Environmental impacts of remediation of a trichloroethene-contaminated site: life cycle assessment of remediation alternatives.

The environmental impacts of remediation of a chloroethene-contaminated site were evaluated using life cycle assessment (LCA). The compared remediation options are (i) in situ bioremediation by enhanced reductive dechlorination (ERD), (ii) in situ thermal desorption (ISTD), and (iii) excavation of the contaminated soil followed by off-site treatment and disposal. The results showed that choosing the ERD option will reduce the life-cycle impacts of remediation remarkably compared to choosing either ISTD or excavation, which are more energy-demanding. In addition to the secondary impacts of remediation, this study includes assessment of local toxic impacts (the primary impact) related to the on-site contaminant leaching to groundwater and subsequent human exposure via drinking water. The primary human toxic impacts were high for ERD due to the formation and leaching of chlorinated degradation products, especially vinyl chloride during remediation. However, the secondary human toxic impacts of ISTD and excavation are likely to be even higher, particularly due to upstream impacts from steel production. The newly launched model, USEtox, was applied for characterization of primary and secondary toxic impacts and combined with a site-dependent fate model of the leaching of chlorinated ethenes from the fractured clay till site.

Evaluation of bioaugmentation with entrapped degrading cells as a soil remediation technology.

Soil augmentation with microbial degraders immobilized on carriers is evaluated as a potential remediation technology using a mathematical model that includes degradation within spatially distributed carriers and diffusion or advection-dispersion as contaminant mass transfer mechanisms. The total volume of carriers is a critical parameter affecting biodegradation performance. In the absence of advection, 320 and 20 000 days are required to mineralize 90% of the herbicide linuron by Variovorax sp. SRS16 encapsulated in 2 mm beads with 5 and 20 mm spacings, respectively. Given that many pesticide degraders have low intrinsic degradation rates and that only limited carrier to soil volume ratios are practically feasible, bioaugmented soils are characterized by low effective degradation rates and can be considered fully mixed. A simple exponential model is then sufficient to predict biodegradation as verified by comparisons with published experimental data. By contrast, the full spatially distributed model is needed to adequately model the degradation of faster degrading contaminants such as naphthalene and benzene which can be mass-transfer limited. Dimensionless Damköhler numbers are proposed to determine whether the spatially distributed model is required. Results show that field scale applications of immobilized degraders will be limited by the amount of carriers required to reach acceptable degradation rates.

Assessing the Transport of Pharmaceutical Compounds in a Layered Aquifer Discharging to a Stream.

A groundwater plume containing high concentrations of pharmaceutical compounds, mainly sulfonamides, barbiturates, and ethyl urethane, in addition to chlorinated ethenes and benzene was investigated. The contamination originating from a former pharmaceutical industry discharges into a multilayered aquifer system and a downgradient stream. In this study, geological and hydrogeological data were integrated into a numerical flow model to examine identified trends using statistical approaches, including principal component analysis and hierarchal cluster analysis. A joint interpretation of the groundwater flow paths and contaminant concentrations in the different compartments (i.e., groundwater and hyporheic zone) provided insight on the transport processes of the different contaminant plumes to the stream. The analysis of historical groundwater concentrations of pharmaceutical compounds at the site suggested these compounds are slowly degrading. The pharmaceutical compounds migrate in both a deep semiconfined aquifer, as well as in the shallow unconfined aquifer, and enter the stream along a 2-km stretch. This contrasted with the chlorinated ethenes, which mainly discharge to the stream as a focused plume from the unconfined aquifer. The integrated approach developed here, combining groundwater flow modeling and statistical analyses of the contaminant concentration data collected in groundwater and the hyporheic zone, lead to an improved understanding of the observed distribution of contaminants in the unconfined and semiconfined aquifers, and thus to their discharge to the stream. This approach is particularly relevant for large and long-lasting contaminant sources and plumes, such as abandoned landfills and industrial production sites, where field investigations may be very expensive.

The valuation of water quality: effects of mixing different drinking water qualities.

As water supplies increasingly turn to use desalination technologies it becomes relevant to consider the options for remineralization and blending with mineral rich water resources. We present a method for analyzing economic consequences due to changes in drinking water mineral content. Included impacts are cardiovascular diseases, dental caries, atopic eczema, lifetime of dish and clothes washing machines, heat exchangers, distribution systems, bottled water consumption and soap usage. The method includes an uncertainty assessment that ranks the impacts having the highest influence on the result and associated uncertainty. Effects are calculated for a scenario where 50% of Copenhagen's water supply is substituted by desalinated water. Without remineralization the total impact is expected to be negative (euro -0.44+/-0.2/m(3)) and individual impacts expected in the range of euro 0.01-0.51/m(3) delivered water. Health impacts have the highest contribution to impact size and uncertainty. With remineralization it is possible to reduce several negative impacts and the total impact is expected to be positive (euro 0.14+/-0.08/m(3)).

A simple contaminant fate and transport modelling tool for management and risk assessment of groundwater pollution from contaminated sites.

Contaminated sites pose a significant threat to groundwater resources. The resources that can be allocated by water regulators for site investigation and cleanup are limited compared to the large number of contaminated sites. Numerical transport models of individual sites require large amounts of data and are labor intensive to set up, and thus they are likely to be too expensive to be useful in the management of thousands of contaminated sites. Therefore, simple tools based on analytical solutions of contaminant transport models are widely used to assess (at an early stage) whether a site might pose a threat to groundwater. We present a tool consisting of five different models, representing common geological settings, contaminant pathways, and transport processes. The tool employs a simplified approach for preliminary, conservative, fast and inexpensive estimation of the contamination levels of aquifers. This is useful for risk assessment applications or to select and prioritize the sites, which should be targeted for further investigation. The tool is based on steady-state semi-analytical models simulating different contaminant transport scenarios from the source to downstream groundwater, and includes both unsaturated and saturated transport processes. The models combine existing analytical solutions from the literature for vertical (from the source to the top of the aquifer) and horizontal (within the aquifer) transport. The effect of net recharge causing a downward migration and an increase of vertical dispersion and dilution of the plume is also considered. Finally, we illustrate the application of the tool for a preliminary assessment of two contaminated sites in Denmark and compare the model results with field data. The comparison shows that a first preliminary assessment with conservative, and often non-site specific parameter selection, is qualitatively consistent with broad trends in observations and provides a conservative estimate of contamination.

A calcite permeable reactive barrier for the remediation of Fluoride from spent potliner (SPL) contaminated groundwater.

The use of calcite (CaCO3) as a substrate for a permeable reactive barrier (PRB) for removing fluoride from contaminated groundwater is proposed and is illustrated by application to groundwater contaminated by spent potliner leachate (SPL), a waste derived from the aluminium smelting process. The paper focuses on two issues in the implementation of calcite permeable reactive barriers for remediating fluoride contaminated water: the impact of the groundwater chemical matrix and CO2 addition on fluoride removal. Column tests comparing pure NaF solutions, synthetic SPL solutions, and actual SPL leachate indicate that the complex chemical matrix of the SPL leachate can impact fluoride removal significantly. For SPL contaminant mixtures, fluoride removal is initially less than expected from idealized, pure, solutions. However, with time, the effect of other contaminants on fluoride removal diminishes. Column tests also show that pH control is important for optimizing fluoride removal with the mass removed increasing with decreasing pH. Barrier pH can be regulated by CO2 addition with the point of injection being critical for optimising the remediation performance. Experimental and model results show that approximately 99% of 2300 mg/L fluoride can be removed when CO2 is injected directly into the barrier. This can be compared to approximately 30-50% removal when the influent solution is equilibrated with atmospheric CO2 before contact with calcite.

Risk assessment and prioritisation of contaminated sites on the catchment scale.

Contaminated sites pose a significant threat to groundwater resources worldwide. Due to limited available resources a risk-based prioritisation of the remediation efforts is essential. Existing risk assessment tools are unsuitable for this purpose, because they consider each contaminated site separately and on a local scale, which makes it difficult to compare the impact from different sites. Hence a modelling tool for risk assessment of contaminated sites on the catchment scale has been developed. The CatchRisk screening tool evaluates the risk associated with each site in terms of its ability to contaminate abstracted groundwater in the catchment. The tool considers both the local scale and the catchment scale. At the local scale, a flexible, site specific leaching model that can be adjusted to the actual data availability is used to estimate the mass flux over time from identified sites. At the catchment scale, a transport model that utilises the source flux and a groundwater model covering the catchment is used to estimate the transient impact on the supply well. The CatchRisk model was tested on a groundwater catchment for a waterworks north of Copenhagen, Denmark. Even though data scarcity limited the application of the model, the sites that most likely caused the observed contamination at the waterworks were identified. The method was found to be valuable as a basis for prioritising point sources according to their impact on groundwater quality. The tool can also be used as a framework for testing hypotheses on the origin of contamination in the catchment and for identification of unknown contaminant sources.

Delay in catchment nitrogen load to streams following restrictions on fertilizer application.

A MIKE SHE hydrological-solute transport model including nitrate reduction is employed to evaluate the delayed response in nitrogen loads in catchment streams following the implementation of nitrogen mitigation measures since the 1980s. The nitrate transport lag times between the root zone and the streams for the period 1950-2011 were simulated for two catchments in Denmark and compared with observational data. Results include nitrogen concentration and mass discharge to streams. By automated baseflow separation, stream discharge was separated into baseflow and drain flow components, and the nitrogen concentration and mass discharge in baseflow and drain flow were determined. This provided insight on the development of stream nitrogen loads, with a short average lag time in drain flow and a long average lag time in baseflow. The long term effect of nitrogen mitigation measures was determined, with results showing that there is a 15 years long delay in the appearance of peak nitrogen loads in streams. This means that real time stream monitoring data cannot be used alone to assess the effect of nitrogen mitigation measures.

A Bayesian belief network approach for assessing uncertainty in conceptual site models at contaminated sites.

A key component in risk assessment of contaminated sites is in the formulation of a conceptual site model (CSM). A CSM is a simplified representation of reality and forms the basis for the mathematical modeling of contaminant fate and transport at the site. The CSM should therefore identify the most important site-specific features and processes that may affect the contaminant transport behavior at the site. However, the development of a CSM will always be associated with uncertainties due to limited data and lack of understanding of the site conditions. CSM uncertainty is often found to be a major source of model error and it should therefore be accounted for when evaluating uncertainties in risk assessments. We present a Bayesian belief network (BBN) approach for constructing CSMs and assessing their uncertainty at contaminated sites. BBNs are graphical probabilistic models that are effective for integrating quantitative and qualitative information, and thus can strengthen decisions when empirical data are lacking. The proposed BBN approach facilitates a systematic construction of multiple CSMs, and then determines the belief in each CSM using a variety of data types and/or expert opinion at different knowledge levels. The developed BBNs combine data from desktop studies and initial site investigations with expert opinion to assess which of the CSMs are more likely to reflect the actual site conditions. The method is demonstrated on a Danish field site, contaminated with chlorinated ethenes. Four different CSMs are developed by combining two contaminant source zone interpretations (presence or absence of a separate phase contamination) and two geological interpretations (fractured or unfractured clay till). The beliefs in each of the CSMs are assessed sequentially based on data from three investigation stages (a screening investigation, a more detailed investigation, and an expert consultation) to demonstrate that the belief can be updated as more information becomes available.

Characterization of chlorinated solvent contamination in limestone using innovative FLUTe® technologies in combination with other methods in a line of evidence approach.

Characterization of dense non-aqueous phase liquid (DNAPL) source zones in limestone aquifers/bedrock is essential to develop accurate site-specific conceptual models and perform risk assessment. Here innovative field methods were combined to improve determination of source zone architecture, hydrogeology and contaminant distribution. The FACT™ is a new technology and it was applied and tested at a contaminated site with a limestone aquifer, together with a number of existing methods including wire-line coring with core subsampling, FLUTe® transmissivity profiling and multilevel water sampling. Laboratory sorption studies were combined with a model of contaminant uptake on the FACT™ for data interpretation. Limestone aquifers were found particularly difficult to sample with existing methods because of core loss, particularly from soft zones in contact with chert beds. Water FLUTe™ multilevel groundwater sampling (under two flow conditions) and FACT™ sampling and analysis combined with FLUTe® transmissivity profiling and modeling were used to provide a line of evidence for the presence of DNAPL, dissolved and sorbed phase contamination in the limestone fractures and matrix. The combined methods were able to provide detailed vertical profiles of DNAPL and contaminant distributions, water flows and fracture zones in the aquifer and are therefore a powerful tool for site investigation. For the limestone aquifer the results indicate horizontal spreading in the upper crushed zone, vertical migration through fractures in the bryozoan limestone down to about 16-18m depth with some horizontal migrations along horizontal fractures within the limestone. Documentation of the DNAPL source in the limestone aquifer was significantly improved by the use of FACT™ and Water FLUTe™ data.

Modeling the Factors Impacting Pesticide Concentrations in Groundwater Wells.

This study examines the effect of pumping, hydrogeology, and pesticide characteristics on pesticide concentrations in production wells using a reactive transport model in two conceptual hydrogeologic systems; a layered aquifer with and without a stream present. The pumping rate can significantly affect the pesticide breakthrough time and maximum concentration at the well. The effect of the pumping rate on the pesticide concentration depends on the hydrogeology of the aquifer; in a layered aquifer, a high pumping rate resulted in a considerably different breakthrough than a low pumping rate, while in an aquifer with a stream the effect of the pumping rate was insignificant. Pesticide application history and properties have also a great impact on the effect of the pumping rate on the concentration at the well. The findings of the study show that variable pumping rates can generate temporal variability in the concentration at the well, which helps understanding the results of groundwater monitoring programs. The results are used to provide guidance on the design of pumping and regulatory changes for the long-term supply of safe groundwater. The fate of selected pesticides is examined, for example, if the application of bentazone in a region with a layered aquifer stops today, the concentration at the well can continue to increase for 20 years if a low pumping rate is applied. This study concludes that because of the rapid response of the pesticide concentration at the drinking water well due to changes in pumping, wellhead management is important for managing pesticide concentrations.

Dilution and volatilization of groundwater contaminant discharges in streams.

An analytical solution to describe dilution and volatilization of a continuous groundwater contaminant plume into streams is developed for risk assessment. The location of groundwater plume discharge into the stream (discharge through the side versus bottom of the stream) and different distributions of the contaminant plume concentration (Gaussian, homogeneous or heterogeneous distribution) are considered. The model considering the plume discharged through the bank of the river, with a uniform concentration distribution was the most appropriate for risk assessment due to its simplicity and limited data requirements. The dilution and volatilization model is able to predict the entire concentration field, and thus the mixing zone, maximum concentration and fully mixed concentration in the stream. It can also be used to identify groundwater discharge zones from in-stream concentration measurement. The solution was successfully applied to published field data obtained in a large and a small Danish stream and provided valuable information on the risk posed by the groundwater contaminant plumes. The results provided by the dilution and volatilization model are very different to those obtained with existing point source models, with a distributed source leading to a larger mixing length and different concentration field. The dilution model can also provide recommendations for sampling locations and the size of impact zones in streams. This is of interest for regulators, for example when developing guidelines for the implementation of the European Water Framework Directive.

Sources, occurrence and predicted aquatic impact of legacy and contemporary pesticides in streams.

We couple current findings of pesticides in surface and groundwater to the history of pesticide usage, focusing on the potential contribution of legacy pesticides to the predicted ecotoxicological impact on benthic macroinvertebrates in headwater streams. Results suggest that groundwater, in addition to precipitation and surface runoff, is an important source of pesticides (particularly legacy herbicides) entering surface water. In addition to current-use active ingredients, legacy pesticides, metabolites and impurities are important for explaining the estimated total toxicity attributable to pesticides. Sediment-bound insecticides were identified as the primary source for predicted ecotoxicity. Our results support recent studies indicating that highly sorbing chemicals contribute and even drive impacts on aquatic ecosystems. They further indicate that groundwater contaminated by legacy and contemporary pesticides may impact adjoining streams. Stream observations of soluble and sediment-bound pesticides are valuable for understanding the long-term fate of pesticides in aquifers, and should be included in stream monitoring programs.

Does microbial centimeter-scale heterogeneity impact MCPA degradation in and leaching from a loamy agricultural soil?

The potential for pesticide degradation varies greatly at the centimeter-scale in agricultural soil. Three dimensional numerical simulations were conducted to evaluate how such small-scale spatial heterogeneity may affect the leaching of the biodegradable pesticide 2-methyl-4-chlorophenoxyacetic acid (MCPA) in the upper meter of a variably-saturated, loamy soil profile. To incorporate realistic spatial variation in degradation potential, we used data from a site where 420 mineralization curves over 5 depths have been measured. Monod kinetics was fitted to the individual curves to derive initial degrader biomass values, which were incorporated in a reactive transport model to simulate heterogeneous biodegradation. Six scenarios were set up using COMSOL Multiphysics to evaluate the difference between models having different degrader biomass distributions (homogeneous, heterogeneous, or no biomass) and either matrix flow or preferential flow through a soil matrix with a wormhole. MCPA leached, within 250 days, below 1m only when degrader biomass was absent and preferential flow occurred. Both biodegradation in the plow layer and the microbially active lining of the wormhole contributed to reducing MCPA-leaching below 1m. The spatial distribution of initial degrader biomass within each soil matrix layer, however, had little effect on the overall MCPA-leaching.