Natural attenuation of sulfonamides and metabolites in contaminated groundwater - Review, advantages and challenges of current documentation techniques

Sulfonamides are applied worldwide as antibiotics. They are emerging contaminants of concern, as their presence in the environment may lead to the spread of antibiotic resistance genes. Sulfonamides are present in groundwater systems, which suggest their persistence under certain conditions, highlighting the importance of understanding natural attenuation processes in groundwater. Biodegradation is an essential process, as degradation of sulfonamides reduces the risk of antibiotic resistance spreading. In this review, natural attenuation, and in particular assessment of biodegradation, is evaluated for sulfonamides in groundwater systems. The current knowledge level on biodegradation is reviewed, and a scientific foundation is built based on sulfonamide degradation processes, pathways, metabolites and toxicity. An overview of bacterial species and related metabolites is provided. The main research effort has focused on aerobic conditions while investigations under anaerobic conditions are lacking. The level of implementation in research is laboratory scale; here we strived to bridge towards field application and assessment, by assessing approaches commonly used in monitored natural attenuation. Methods to document contaminant mass loss are assessed to be applicable for sulfonamides, while the approach is limited by a lack of reference standards for metabolites. Furthermore, additional information is required on relevant metabolites in order to improve risk assessments. Based on the current knowledge on biodegradation, it is suggested to use the presence of substituent-containing metabolites from breakage of the sulfonamide bridge as specific indicators of degradation. Microbial approaches are currently available for assessment of microbial community's capacities, however, more knowledge is required on indigenous bacteria capable of degrading sulfonamides and on the impact of environmental conditions on biodegradation. Compound specific stable isotope analysis shows great potential as an additional in situ method, but further developments are required to analyse for sulfonamides at environmentally relevant levels. Finally, in a monitored natural attenuation scheme it is assessed that approaches are available that can uncover some processes related to the fate of sulfonamides in groundwater systems. Nevertheless, there are still unknowns related to relevant bacteria and metabolites for risk assessment as well as the effect of environmental settings such as redox conditions. Alongside, uncovering the fate of sulfonamides in future research, the applicability of the natural attenuation documentation approaches will advance, and provide a step towards in situ remedial concepts for the frequently detected sulfonamides.

Quantification of contaminant mass discharge from point sources in aquitard/aquifer systems based on vertical concentration profiles and 3D modeling

Point sources with contaminants, such as chlorinated solvents, per- and polyfluoroalkyl substances (PFAS), or pesticides, are often located in low-permeability aquitards, where they can act as long-term sources and threaten underlying groundwater resources. We demonstrate the use of a 3D numerical model integrating comprehensive hydrogeological and contamination data to determine the contaminant mass discharge (CMD) from an aquitard into the underlying aquifer. A mature point source with a dissolved chlorinated solvent in a clayey till is used as an example. The quantitative determination is facilitated by model calibration to high-resolution vertical concentration profiles obtained by direct-push sampling techniques in the aquifer downgradient of the contaminant source zone. The concentration profiles showed a plume sinking with distance from the source characteristic for such aquitard/aquifer settings. The sinking is caused by the interplay between infiltrating water and horizontal groundwater flow. The application of 3D solute transport modeling on high-resolution profiles allowed for determining the infiltration rate, the hydraulic conductivity in the aquitard, and, ultimately, the CMD. Different source zone conceptualizations demonstrate the potential effects of fractures and sorption in source zones in aquitards on CMD development. Fractures in the aquitard had a minor influence on the current CMD determined with the presented approach. Still, fractures with hydraulic apertures larger than 10 μm were crucial for the temporal development of the CMD and plume. A thorough characterization of the source zone conditions combined with high-resolution concentration profiles and detailed modeling is valuable for shedding light on the probable future development of groundwater contamination arising from sources in aquitard/aquifer settings and evaluating remedial actions.

A novel concept for estimating the contaminant mass discharge of chlorinated ethenes emanating from clay till sites

Interest in using contaminant mass discharge (CMD) for risk assessment of contaminated sites has increased over the years, as it accounts for the contaminant mass that is moving and posing a risk to water resources and receptors. The most common investigation of CMD involves a transect of multilevel wells; however, this is an expensive undertaking, and it is difficult to place it in the right position in a plume. Additionally, infrastructure at the site needs to be considered. To derive an initial CMD estimate at a contaminated site and to allow for the prioritization of further investigations and remedial actions, the ProfileFlux method has been developed. It is targeted at former industrial sites with a source zone in a low conductivity layer with primarily vertical flow overlying an aquifer with primarily horizontal groundwater flow. The ProfileFlux method was developed for mature chlorinated solvent plumes, typically originating from more than 30 to 50-year-old spills, as the usage of chlorinated solvents is primarily historical. Thus, it is assumed that the contaminant had time to distribute in the low conductivity layer by mainly diffusive processes. Today the contamination is continuously released to the underlying aquifer, where advection and dispersive (other than diffusive) processes are of higher importance. The approach combines high-resolution, depth-discrete vertical concentration profiles and a simple 2D flow and transport model to estimate CMD by comparing measured and simulated concentration profiles. The study presented herein includes a global sensitivity analysis, in order to identify crucial field parameters, and of particular importance in this regard are source length, groundwater flux and infiltration. The ProfileFlux method was tested at a well-examined industrial site primarily contaminated with trichloroethylene, thereby allowing a comparison between CMD from the ProfileFlux method and the traditional transect method. CMD was estimated at 117-170 g/year, when using the ProfileFlux method, against 143 g/year with the transect method, thus validating ProfileFlux method's ability to estimate CMD. In addition, applying the method identified weak points in the conceptual site model. The method will be incorporated into a user-friendly online tool directed at environmental consultants and decision-makers working on the risk assessment and prioritization of contaminated sites with the specific hydrogeological conditions of an aquifer with an overlying low permeability layer.

Dataset for laboratory treatability experiment with activated carbon and bioamendments to enhance biodegradation of chlorinated ethenes

This dataset describes the outcome of a laboratory trichloroethene (TCE) treatability experiment with liquid activated carbon and bioamendments. The treatability experiment included unamended microcosms, bioamended microcosms with a Dehalococcoides containing culture and electron donor, and bioamended microcosms including liquid activated carbon (PlumeStop®). Data were collected frequently over an 85-day experimental period. Data were collected for the following parameters: redox sensitive species, chlorinated ethenes, non-chlorinated end-products, electron donors, compound specific isotopes, specific bacteria and functional genes. The reductive dechlorination of TCE could be described by a carbon isotope enrichment factor (ε_C) of -7.1 ‰. In the amended systems, the degradation rates for the TCE degradation were 0.08–0.13 d⁻¹ and 0.05–0.09 d⁻¹ determined by concentrations and isotope fractionation, respectively. Dechlorination of cis-DCE was limited. This dataset assisted in identifying the impact of different bioamendments and activated carbon on biodegradation of chlorinated ethenes. The dataset is useful in optimising design and setup for future laboratory and field investigations. This study provides novel information on the effect of low dose liquid activated carbon on chlorinated ethenes degradation by applying isotopic and microbial techniques, and by linking the outcome to a field case study. The data presented in this article are related to the research article “Assessment of chlorinated ethenes degradation after field scale injection of activated carbon and bioamendments: Application of isotopic and microbial analyses” (Ottosen et al., 2021).

Assessment of chlorinated ethenes degradation after field scale injection of activated carbon and bioamendments: Application of isotopic and microbial analyses

Over the last decade, activated carbon amendments have successfully been applied to retain chlorinated ethene subsurface contamination. The concept of this remediation technology is that activated carbon and bioamendments are injected into aquifer systems to enhance biodegradation. While the scientific basis of the technology is established, there is a need for methods to characterise and quantify the biodegradation at field scale. In this study, an integrated approach was applied to assess in situ biodegradation after the establishment of a cross sectional treatment zone in a TCE plume. The amendments were liquid activated carbon, hydrogen release donors and a Dehalococcoides containing culture. The integrated approach included spatial and temporal evaluations on flow and transport, redox conditions, contaminant concentrations, biomarker abundance and compound-specific stable isotopes. This is the first study applying isotopic and microbial techniques to assess field scale biodegradation enhanced by liquid activated carbon and bioamendments. The injection enhanced biodegradation from TCE to primarily cis-DCE. The Dehalococcoides abundances facilitated characterisation of critical zones with insufficient degradation and possible explanations. A conceptual model of isotopic data together with distribution and transport information improved process understanding; the degradation of TCE was insufficient to counteract the contaminant input by inflow into the treatment zone and desorption from the sediment. The integrated approach could be used to document and characterise the in situ degradation, and the isotopic and microbial data provided process understanding that could not have been gathered from conventional monitoring tools. However, quantification of degradation through isotope data was restricted for TCE due to isotope masking effects. The combination of various monitoring tools, applied frequently at high-resolution, with system understanding, was essential for the assessment of biodegradation in the complex, non-stationary system. Furthermore, the investigations revealed prospects for future research, which should focus on monitoring contaminant fate and microbial distribution on the sediment and the activated carbon.

Natural attenuation of a chlorinated ethene plume discharging to a stream: Integrated assessment of hydrogeological, chemical and microbial interactions

Attenuation processes of chlorinated ethenes in complex near-stream systems result in site-specific outcomes of great importance for risk assessment of contaminated sites. Additional interdisciplinary and comprehensive field research is required to enhance process understanding in these systems. In this study, several methods were combined in a multi-scale interdisciplinary in-situ approach to assess and quantify the near-stream attenuation of a chlorinated ethene plume, mainly consisting of cis-dichloroethene (cis-DCE) and vinyl chloride (VC), discharging to a lowland stream (Grindsted stream, Denmark) over a monitoring period of seven years. The approach included: hydrogeological characterisation, reach scale contaminant mass balance analysis, quantification of contaminant mass discharge, streambed fluxes of chlorinated ethenes quantified using Sediment Bed Passive Flux Meters (SBPFMs), assessment of redox conditions, temporal assessment of contaminant concentrations, microbial analysis, and compound-specific isotope analysis (CSIA). This study site exhibits a special attenuation behaviour not commonly encountered in field studies: the conversion from an initially limited degradation case (2012-16), despite seemingly optimal conditions, to one presenting notable levels of degradation (2019). Hence, this study site provides a new piece to the puzzle, as sites with different attenuation behaviours are required in order to acquire the full picture of the role groundwater-surface water interfaces have in risk mitigation. In spite of the increased degradation in the near-stream plume core, the contaminant attenuation was still incomplete in the discharging plume. A conceptualization of flow, transport and processes clarified that hydrogeology was the main control on the natural attenuation, as short residence times of 0.5-37 days restricted the time in which dechlorination could occur. This study reveals the importance of: taking an integrated approach to understand the influence of all attenuation processes in groundwater - surface water interactions; considering the scale and domain of interest when determining the main processes; and monitoring sufficiently both spatially and temporally to cover the transient conditions.

Transport of citrate and polymer coated gold nanoparticles (AuNPs) in porous media: Effect of surface property and Darcy velocity

With the release of nanoparticles (NPs) into the subsurface, it is imperative to better understand the fate and transport of NPs in porous media. Three types of stable AuNPs were used as model NPs to investigate the impact of surface coatings (type and coverage) and water velocity on the NP transport in a porous media (column studies). The NPs were electrostatic stabilized citrate AuNPs and sterically stabilized AuNPs with amphiphilic block co-polymer (PVA-COOH) in two particle/polymer ratios (weak vs. strong stabilization). The citrate AuNPs transport was sensitive to ionic changes in the mixing front of the plume, where destabilization occurred, and will therefore depend on the size/type of release. Blocking of deposition sites by aggregates was seen to facilitate transport, whereby a higher flow velocity (larger shadow zone) also resulted in better transport. The polymeric surface coating had great impact with steric repulsion as a main force contributing to the transport of NPs in the porous media. Sufficient polymer coating was crucial to obtain highly unfavorable attachment conditions (very low α) where the enhanced NP mobility was independent of the water velocity (comparable to solute tracer). Without sufficient steric stabilization, the transport and recovery was significantly reduced compared to the solute tracer, but increased with increasing water velocity. This highlights the importance of sufficient surface coating to achieve enhanced mobility, but also the increased risk of spreading to down-gradient receptors. For the (weakly) sterically stabilized NPs, the loss of polymer through ligand exchange with the porous media negates transport.

Chlorinated ethene plume evolution after source thermal remediation: Determination of degradation rates and mechanisms

The extent, mechanism(s), and rate of chlorinated ethene degradation in a large tetrachloroethene (PCE) plume were investigated in an extensive sampling campaign. Multiple lines of evidence for this degradation were explored, including compound-specific isotope analysis (CSIA), dual C-Cl isotope analysis, and quantitative real-time polymerase chain reaction (qPCR) analysis targeting the genera Dehalococcoides and Dehalogenimonas and the genes vcrA, bvcA, and cerA. A decade prior to this sampling campaign, the plume source was thermally remediated by steam injection. This released dissolved organic carbon (DOC) that stimulated microbial activity and created reduced conditions within the plume. Based on an inclusive analysis of minor and major sampling campaigns since the initial site characterization, it was estimated that reduced conditions peaked 4 years after the remediation event. At the time of this study, 11 years after the remediation event, the redox conditions in the aquifer are returning to their original state. However, the DOC released from the remediated source zone matches levels measured 3 years prior and plume conditions are still suitable for biotic reductive dechlorination. Dehalococcoides spp., Dehalogenimonas spp., and vcrA, bvcA, and cerA reductive dehalogenase genes were detected close to the source, and suggest that complete, biotic PCE degradation occurs here. Further downgradient, qPCR analysis and enriched δ¹³C values for cis-dichloroethene (cDCE) suggest that cDCE is biodegraded in a sulfate-reducing zone in the plume. In the most downgradient portion of the plume, lower levels of specific degraders supported by dual C-Cl analysis indicate that the biodegradation occurs in combination with abiotic degradation. Additionally, 16S rRNA gene amplicon sequencing shows that organizational taxonomic units known to contain organohalide-respiring bacteria are relatively abundant throughout the plume. Hydraulic conductivity testing was also conducted, and local degradation rates for PCE and cDCE were determined at various locations throughout the plume. PCE degradation rates from sampling campaigns after the thermal remediation event range from 0.11 to 0.35 yr^-1. PCE and cDCE degradation rates from the second to the third sampling campaigns ranged from 0.08 to 0.10 yr^-1 and 0.01 to 0.07 yr^-1, respectively. This is consistent with cDCE as the dominant daughter product in the majority of the plume and cDCE degradation as the time-limiting step. The extensive temporal and spatial analysis allowed for tracking the evolution of the plume and the lasting impact of the source remediation and illustrates that the multiple lines of evidence approach is essential to elucidate the primary degradation mechanisms in a plume of such size and complexity.

Mobility of electrostatically and sterically stabilized gold nanoparticles (AuNPs) in saturated porous media

The stability of gold nanoparticles (AuNPs) stabilized electrostatically with citrate or (electro)sterically by commercially available amphiphilic block copolymers (PVP-VA or PVA-COOH) was studied under various physicochemical conditions. Subsequently, the mobility of the AuNPs in porous media (sand) was investigated in column studies under environmental relevant physicochemical conditions. Electrostatically stabilized AuNPs were unstable under most physicochemical conditions due to the compression of the electrical double layer. Consequently, aggregation and deposition rapidly immobilized the AuNPs. Sterically stabilized AuNPs showed significantly less sensitivity towards changes in the physicochemical conditions with high stability, high mobility with negligible retardation, and particle deposition rate coefficients ranging an order of magnitude (1.5 × 10^-3 to 1.5 × 10^-2 min^-1) depending on the type and amount of stabilizer, and thereby the surface coverage and attachment affinity. The transport of sterically stabilized AuNPs is facilitated by reversible deposition in shallow energy minima with continuous reentrainment and blocking of available attachment sites by deposited AuNPs. The stability and mobility of NPs in the environment will thereby be highly dependent on the specific stabilizing agent and variations in the coverage on the NP. Under the given experimental conditions, transport distances of the most mobile AuNPs of up to 20 m is expected. Due to their size-specific plasmonic properties, the easily detectable AuNPs are proposed as potential model or tracer particles for studying transport of various stabilized NPs under environmental conditions.

A modeling approach integrating microbial activity, mass transfer, and geochemical processes to interpret biological assays: An example for PCE degradation in a multi-phase batch setup

The rate at which organic contaminants can be degraded in aquatic environments is not only dependent upon specific degrading bacteria, but also upon the composition of the microbial community, mass transfer of the contaminant, and abiotic processes that occur in the environment. In this study, we present three-phase batch experiments of tetrachloroethene (PCE) degradation by a consortium of organohalide-respiring bacteria, cultivated alone or in communities with iron- and/or sulfate-reducers. We developed a modeling approach to quantitatively evaluate the experimental results, comprised of chemical and biomolecular time series data. The model utilizes the IPhreeqc module to couple multi-phase mass transfer between gaseous, organic and aqueous phases with microbial and aquatic geochemical processes described using the geochemical code PHREEQC. The proposed approach is able to capture the contaminant degradation, the microbial population dynamics, the effects of multi-phase kinetic mass transfer and sample removal, and the geochemical reactions occurring in the aqueous phase. The model demonstrates the importance of aqueous speciation and abiotic reactions on the bioavailability of the substrates. The model-based interpretation allowed us to quantify the reaction kinetics of the different bacterial guilds. The model further revealed that the inclusion of sulfate-reducing bacteria lowers the rate of PCE degradation and that this effect is moderated in the presence of iron-reducing bacteria.

Identification of abiotic and biotic reductive dechlorination in a chlorinated ethene plume after thermal source remediation by means of isotopic and molecular biology tools

Thermal tetrachloroethene (PCE) remediation by steam injection in a sandy aquifer led to the release of dissolved organic carbon (DOC) from aquifer sediments resulting in more reduced redox conditions, accelerated PCE biodegradation, and changes in microbial populations. These changes were documented by comparing data collected prior to the remediation event and eight years later. Based on the premise that dual C-Cl isotope slopes reflect ongoing degradation pathways, the slopes associated with PCE and TCE suggest the predominance of biotic reductive dechlorination near the source area. PCE was the predominant chlorinated ethene near the source area prior to thermal treatment. After thermal treatment, cDCE became predominant. The biotic contribution to these changes was supported by the presence of Dehalococcoides sp. DNA (Dhc) and Dhc targeted rRNA close to the source area. In contrast, dual C-Cl isotope analysis together with the almost absent VC (13)C depletion in comparison to cDCE (13)C depletion suggested that cDCE was subject to abiotic degradation due to the presence of pyrite, possible surface-bound iron (II) or reduced iron sulphides in the downgradient part of the plume. This interpretation is supported by the relative lack of Dhc in the downgradient part of the plume. The results of this study show that thermal remediation can enhance the biodegradation of chlorinated ethenes, and that this effect can be traced to the mobilisation of DOC due to steam injection. This, in turn, results in more reduced redox conditions which favor active reductive dechlorination and/or may lead to a series of redox reactions which may consecutively trigger biotically induced abiotic degradation. Finally, this study illustrates the valuable complementary application of compound-specific isotopic analysis combined with molecular biology tools to evaluate which biogeochemical processes are taking place in an aquifer contaminated with chlorinated ethenes.

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.

Stable carbon isotope analysis to distinguish biotic and abiotic degradation of 1,1,1-trichloroethane in groundwater sediments

The fate and treatability of 1,1,1-TCA by natural and enhanced reductive dechlorination was studied in laboratory microcosms. The study shows that compound-specific isotope analysis (CSIA) identified an alternative 1,1,1-TCA degradation pathway that cannot be explained by assuming biotic reductive dechlorination. In all biotic microcosms 1,1,1-TCA was degraded with no apparent increase in the biotic degradation product 1,1-DCA. 1,1,1-TCA degradation was documented by a clear enrichment in (13)C in all biotic microcosms, but not in the abiotic control, which suggests biotic or biotically mediated degradation. Biotic degradation by reductive dechlorination of 1,1-DCA to CA only occurred in bioaugmented microcosms and in donor stimulated microcosms with low initial 1,1,1-TCA or after significant decrease in 1,1,1-TCA concentration (after∼day 200). Hence, the primary degradation pathway for 1,1,1-TCA does not appear to be reductive dechlorination via 1,1-DCA. In the biotic microcosms, the degradation of 1,1,1-TCA occurred under iron and sulfate reducing conditions. Biotic reduction of iron and sulfate likely resulted in formation of FeS, which can abiotically degrade 1,1,1-TCA. Hence, abiotic degradation of 1,1,1-TCA mediated by biotic FeS formation constitute an explanation for the observed 1,1,1-TCA degradation. This is supported by a high 1,1,1-TCA (13)C enrichment factor consistent with abiotic degradation in biotic microcosms. 1,1-DCA carbon isotope field data suggest that this abiotic degradation of 1,1,1-TCA is a relevant process also at the field site.

The impact of bioaugmentation on dechlorination kinetics and on microbial dechlorinating communities in subsurface clay till

A molecular study on how the abundance of the dechlorinating culture KB-1 affects dechlorination rates in clay till is presented. DNA extracts showed changes in abundance of specific dechlorinators as well as their functional genes. Independently of the KB-1 added, the microbial dechlorinator abundance increased to the same level in all treatments. In the non-bioaugmented microcosms the reductive dehalogenase gene bvcA increased in abundance, but when KB-1 was added the related vcrA gene increased while bvcA genes did not increase. Modeling showed higher vinyl-chloride dechlorination rates and shorter time for complete dechlorination to ethene with higher initial concentration of KB-1 culture, while cis-dichloroethene dechlorination rates were not affected by KB-1 concentrations. This study provides high resolution abundance profiles of Dehalococcoides spp. (DHC) and functional genes, highlights the ecological behavior of KB-1 in clay till, and reinforces the importance of using multiple functional genes as biomarkers for reductive dechlorination.

Identification of chlorinated solvents degradation zones in clay till by high resolution chemical, microbial and compound specific isotope analysis

The degradation of chlorinated ethenes and ethanes in clay till was investigated at a contaminated site (Vadsby, Denmark) by high resolution sampling of intact cores combined with groundwater sampling. Over decades of contamination, bioactive zones with degradation of trichloroethene (TCE) and 1,1,1-trichloroethane (1,1,1-TCA) to 1,2-cis-dichloroethene (cis-DCE) and 1,1-dichloroethane, respectively, had developed in most of the clay till matrix. Dehalobacter dominated over Dehalococcoides (Dhc) in the clay till matrix corresponding with stagnation of sequential dechlorination at cis-DCE. Sporadically distributed bioactive zones with partial degradation to ethene were identified in the clay till matrix (thickness from 0.10 to 0.22 m). In one sub-section profile the presence of Dhc with the vcrA gene supported the occurrence of degradation of cis-DCE and VC, and in another enriched δ(13)C for TCE, cis-DCE and VC documented degradation. Highly enriched δ(13)C for 1,1,1-TCA (25‰) and cis-DCE (-4‰) suggested the occurrence of abiotic degradation in a third sub-section profile. Due to fine scale heterogeneity the identification of active degradation zones in the clay till matrix depended on high resolution subsampling of the clay till cores. The study demonstrates that an integrated approach combining chemical analysis, molecular microbial tools and compound specific isotope analysis (CSIA) was required in order to document biotic and abiotic degradations in the clay till system.

A remediation performance model for enhanced metabolic reductive dechlorination of chloroethenes in fractured clay till

A numerical model of metabolic reductive dechlorination is used to describe the performance of enhanced bioremediation in fractured clay till. The model is developed to simulate field observations of a full scale bioremediation scheme in a fractured clay till and thereby to assess remediation efficiency and timeframe. A relatively simple approach is used to link the fermentation of the electron donor soybean oil to the sequential dechlorination of trichloroethene (TCE) while considering redox conditions and the heterogeneous clay till system (clay till matrix, fractures and sand stringers). The model is tested on lab batch experiments and applied to describe sediment core samples from a TCE-contaminated site. Model simulations compare favorably to field observations and demonstrate that dechlorination may be limited to narrow bioactive zones in the clay matrix around fractures and sand stringers. Field scale simulations show that the injected donor is expected to be depleted after 5 years, and that without donor re-injection contaminant rebound will occur in the high permeability zones and the mass removal will stall at 18%. Long remediation timeframes, if dechlorination is limited to narrow bioactive zones, and the need for additional donor injections to maintain dechlorination activity may limit the efficiency of ERD in low-permeability media. Future work should address the dynamics of the bioactive zones, which is essential to understand for predictions of long term mass removal.

Dissolution and transport of coal tar compounds in fractured clay-rich residuum

We investigated the dissolution and transport of organic contaminants from a crude coal tar mixture in a monolith of fractured clay-rich residuum. An electrolyte solution was eluted through the residuum monolith containing a small emplaced source of coal tar under biologically inhibited and mildly acidic conditions. Concentrations of 10 coal tar compounds, representing mono-, poly-, and heterocyclic aromatic hydrocarbons that constitute crude coal tar were monitored in the effluent over a period of 377 days. Most compounds appeared in the effluent within the first 0.1 pore volume eluted indicating the importance of rapid dissolution and transport through the fracture networks. The concentrations continued to rise but did not reach the corresponding effective solubility limit in most cases. Compounds that were less soluble and those that were more susceptible to sorption or matrix diffusion eluted at a much slower rate. Analysis of contaminant concentrations in microcore residuum samples indicated that all 10 compounds had spread throughout the entire monolith and had diffused into the fine-grained matrix between fractures. These data suggest that the predominantly fine pore structure did not appear to inhibit coal tar dissolution and subsequent transport, even though only a small portion of tar was in direct contact with fractures and macropores that control most flow.

Sorption of chlorinated solvents and degradation products on natural clayey tills

The sorption of chlorinated solvents and degradation products on seven natural clayey till samples from three contaminated sites was investigated by laboratory batch experiments in order to obtain reliable sorption coefficients (K(d) values). The sorption isotherms for all compounds were nearly linear, but fitted by Freundlich isotherms slightly better over the entire concentration range. For chloroethylenes, tetrachloroethylene (PCE) was most strongly sorbed to the clayey till samples (K(d)=0.84-2.45Lkg(-1)), followed by trichloroethylene (TCE, K(d)=0.62-0.96Lkg(-1)), cis-dichloroethylene (cis-DCE, K(d)=0.17-0.82Lkg(-1)) and vinyl chloride (VC, K(d)=0.12-0.36Lkg(-1)). For chloroethanes, 1,1,1-trichloroethane (1,1,1-TCA) was most strongly sorbed (K(d)=0.2-0.45Lkg(-1)), followed by 1,1-dichloroethane (1,1-DCA, K(d)=0.16-0.24Lkg(-1)) and chloroethane (CA, K(d)=0.12-0.18Lkg(-1)). This is consistent with the order of hydrophobicity of the compounds. The octanol-water coefficient (logK(ow)) correlated slightly better with logK(d) values than logK(oc) values indicating that the K(d) values may be independent of the actual organic carbon content (f(oc)). The estimated logK(oc) or logK(d) for chlorinated solvents and degradation products determined by regression of data in this study were significantly higher than values determined by previously published empirical relationships. The site specific K(d) values as well as the new empirical relationship compared well with calculations on water and soil core concentration for cis-DCE and VC from the Rugårdsvej site. In conclusion, this study with a wide range of chlorinated ethenes and ethanes--in line with previous studies on PCE and TCE--suggest that sorption in clayey tills could be higher than typically expected.

Assessing chlorinated ethene degradation in a large scale contaminant plume by dual carbon-chlorine isotope analysis and quantitative PCR

The fate of chlorinated ethenes in a large contaminant plume originating from a tetrachloroethene (PCE) source in a sandy aquifer in Denmark was investigated using novel methods including compound-specific carbon and chlorine isotope analysis and quantitative real-time polymerase chain reaction (qPCR) methods targeting Dehaloccocoides sp. and vcrA genes. Redox conditions were characterized as well based on concentrations of dissolved redox sensitive compounds and sulfur isotopes in SO(4)(2-). In the first 400 m downgradient of the source, the plume was confined to the upper 20 m of the aquifer. Further downgradient it widened in vertical direction due to diverging groundwater flow reaching a depth of up to 50 m. As the plume dipped downward and moved away from the source, O(2) and NO(3)(-) decreased to below detection levels, while dissolved Fe(2+) and SO(4)(2-) increased above detectable concentrations, likely due to pyrite oxidation as confirmed by the depleted sulfur isotope signature of SO(4)(2-). In the same zone, PCE and trichloroethene (TCE) disappeared and cis-1,2-dichloroethene (cDCE) became the dominant chlorinated ethene. PCE and TCE were likely transformed by reductive dechlorination rather than abiotic reduction by pyrite as indicated by the formation of cDCE and stable carbon isotope data. TCE and cDCE showed carbon isotope trends typical for reductive dechlorination with an initial depletion of (13)C in the daughter products followed by an enrichment of (13)C as degradation proceeded. At 1000 m downgradient of the source, cDCE was the dominant chlorinated ethene and had reached the source δ(13)C value confirming that cDCE was not affected by abiotic or biotic degradation. Further downgradient (up to 1900 m), cDCE became enriched in (13)C by up to 8 ‰ demonstrating its further transformation while vinylchloride (VC) concentrations remained low (<1 μg/L) and ethene was not observed. The correlated shift of carbon and chlorine isotope ratios of cDCE by 8 and 3.9 ‰, respectively, the detection of Dehaloccocides sp genes, and strongly reducing conditions in this zone provide strong evidence for reductive dechlorination of cDCE. The significant enrichment of (13)C in VC indicates that VC was transformed further, although the mechanism could not be determined. The transformation of cDCE was the rate limiting step as no accumulation of VC occurred. In summary, the study demonstrates that carbon-chlorine isotope analysis and qPCR combined with traditional approaches can be used to gain detailed insight into the processes that control the fate of chlorinated ethenes in large scale plumes.

Field evaluation of biological enhanced reductive dechlorination of chloroethenes in clayey till

The performance of enhanced reductive dechlorination (ERD) for in situ remediation of cis-1,2-dichloroethene (cDCE) and vinyl chloride in clayey till was investigated in a pilot test. A dilute groundwater solution containing emulsified soybean oil and Dehalococcoides bacteria was injected into a sand-filled hydraulic fracture. Fermentation of the ERD solution caused the establishment of a dechlorinating bioactive zone in the fracture within 1 month of injection. By 148 days, all the cDCE in the fracture was dechlorinated to ethene. Analysis of a clay core from Day 150 indicated that electron donor and fermentation products diffused from the fracture at least 10 cm into clay and that stimulated dechlorination occurred in the clay in the presence of Dehalococcoides (7.9.10(4) cells g(-1)). Comparison of chloroethene profiles in the Day 150 core to modeled diffusion profiles indicated degradation occurred in a bioactive zone extending approximately 5 to 6 cm into the clay matrix. These data suggest that a bioactive zone established in a sand-filled fracture can expand into the adjacent clayey till matrix and facilitate mass transfer from the matrix to the bioactive zone. These findings offer promise for ERD and support further development of methods for deploying ERD in clayey till and other low-permeability deposits.

Modeling multi-component transport and enhanced anaerobic dechlorination processes in a single fracture-clay matrix system

Clayey tills contaminated with chlorinated solvents are a threat to groundwater and are difficult to remediate. A numerical model is developed for assessing leaching processes and for simulating the remediation via enhanced anaerobic dechlorination. The model simulates the transport of a contaminant in a single fracture-clay matrix system coupled with a reactive model for anaerobic dechlorination. The model takes into account microbially driven anaerobic dechlorination, where sequential Monod kinetics with competitive inhibition is used to model the reaction rates, and degradation is localized to account for potential pore size limitations on microbial entry to the clay matrix. The model is used to assess the distribution of TCE and its daughter products in the clay matrix and the concentration of the different compounds at the outlet of the fracture. The time frame for complete cleanup and the contaminant flux out of the clay system are assessed for different distributions of microbial degradation. Results from a set of scenarios show that time to remove 90% of the initial mass is halved when dechlorination occurs in a 5cm reaction zone in the clay at the fracture-matrix interface (from 419 to 195years) and decreases by an order of magnitude when dechlorination occurs in the entire matrix (to 32years). The fracture spacing and the microbial parameters are shown to be the critical parameter for estimation of time frames depending on the system in question. Generally, the system is more sensitive to the physical processes, mainly diffusion in the matrix, than to the biogeochemical processes, when dechlorination is assumed to take place in a limited reaction zone only. The inclusion of sequential dechlorination in clay fracture transport models is crucial, as the contaminant flux to the aquifer will increase as a result of degradation due to the higher mobility of the formed daughter products DCE and VC. The model is used to examine the relationship between flux reduction and mass removal for fractured clay systems.

Carbon isotope fractionation during diffusion and biodegradation of petroleum hydrocarbons in the unsaturated zone: field experiment at Vaerløse Airbase, Denmark, and modeling

A field experiment was conducted in Denmark in order to evaluate the fate of 13 volatile organic compounds (VOCs) that were buried as an artificial fuel source in the unsaturated zone. Compound-specific isotope analysis showed distinct phases in the 13C/12C ratio evolution in VOC vapors within 3 m from the source over 114 days. At day 3 and to a lesser extent at day 6, the compounds were depleted in 13C by up to -5.7% per hundred with increasing distance from the source compared to the initial source values. This trend can be explained by faster outward diffusion of the molecules with 12C only compared to molecules with a 13C. Then, the isotope profile leveled out, and several compounds started to become enriched in 13C by up to 9.5% per hundred with increasing distance from the source, due to preferential removal of the molecules with 12C only, through biodegradation. Finally, as the amount of a compound diminished in the source, a 13C enrichment was also observed close to the source. The magnitude of isotope fractionation tended to be larger the smaller the mass of the molecule was. This study demonstrates that, in the unsaturated zone, carbon isotope ratios of hydrocarbons are affected by gas-phase diffusion in addition to biodegradation, which was confirmed using a numerical model. Gas-phase diffusion led to shifts in delta(13)C >1% per hundred during the initial days after the spill, and again during the final stages of source volatilization after >75% of a compound had been removed. In between, diffusion has less of an effect, and thus isotope data can be used as an indicator for hydrocarbon biodegradation.

Characterization and quantification of pneumatic fracturing effects at a clay till site

Environmental fracturing offers assistance to remediation efforts at contaminated, low-permeability sites via creation of active fracture networks, and hence, reduction of mass transport limitations set by diffusion in low-permeability matrices. A pilot study of pneumatic fracturing, focusing on direct documentation of fracture propagation patterns and spacing, was performed at a typical basal clay till site. The study applied a novel package of documentation methods, including injection of five tracers with different characteristics (bromide, uvitex, fluorescein, rhodamine WT, and brilliant blue), subsequent tracer-filled fracture documentation via direct and indirect methods, and geological characterization of the fractured site. The direct documentation methods consisted of Geoprobe coring, augering, and excavation. A mass balance and conceptual model have been established for the distribution of the injected tracers in the subsurface. They reveal that tracer was distributed within 2 m of the fracturing well, mainly in existing fractures above the redox boundary (2 to 4 m.b.s.; 5 to 10 cm spacing). Spacing of observed tracer-filled fractures was large (>1 m) at greater depths. The number of fractures induced/activated could possibly be increased via adjustments to the fracturing equipment design.

Role of diffusion in chemical oxidation of PCE in a dual permeability system

The effect of in situ chemical oxidation (ISCO) with permanganate (MnO4-) on tetrachloroethene (PCE) in a dual permeability system consisting of low permeability clay with high permeability sand lenses was investigated by two-dimensional laboratory experiments. The experiments imitate a field remediation at a former dry cleaning facility in Denmark. Results from laboratory experiments and field observations both show that after an application of MnO4- in the sand layer, the diffusion rate into the matrix is decreased due to reaction with PCE and the natural oxidant demand (NOD) related to the clayey till. A narrow but very efficient reaction zone is created in the clayey till. Initiallythe zone developed rapidlyfollowed by a slower expansion with time. PCE will counter diffuse into the reaction zone, where it will be degraded as long as MnO4- is present. A mass balance for the laboratory experiment revealed that the reaction between MnO4- and the clayey till was responsible for up to 90% of the total MnO4- consumption. Based on laboratory experiments, the high MnO4- consumption from reaction with clayey till appears to have been the limiting factor for the oxidation of PCE at the field site and is thus impairing the efficiency of ISCO as a remedy.

Quantification of potassium permanganate consumption and PCE oxidation in subsurface materials

A series of laboratory scale batch slurry experiments were conducted in order to establish a data set for oxidant demand by sandy and clayey subsurface materials as well as to identify the reaction kinetic rates of permanganate (MnO(4)(-)) consumption and PCE oxidation as a function of the MnO(4)(-) concentration. The laboratory experiments were carried out with 31 sandy and clayey subsurface sediments from 12 Danish sites. The results show that the consumption of MnO(4)(-) by reaction with the sediment, termed the natural oxidant demand (NOD), is the primary reaction with regards to quantification of MnO(4)(-) consumption. Dissolved PCE in concentrations up to 100 mg/l in the sediments investigated is not a significant factor in the total MnO(4)(-) consumption. Consumption of MnO(4)(-) increases with an increasing initial MnO(4)(-) concentration. The sediment type is also important as NOD is (generally) higher in clayey than in sandy sediments for a given MnO(4)(-) concentration. For the different sediment types the typical NOD values are 0.5-2 g MnO(4)(-)/kg dry weight (dw) for glacial meltwater sand, 1-8 g MnO(4)(-)/kg dw for sandy till and 5-20 g MnO(4)(-)/kg dw for clayey till. The long term consumption of MnO(4)(-) and oxidation of PCE can not be described with a single rate constant, as the total MnO(4)(-) reduction is comprised of several different reactions with individual rates. During the initial hours of reaction, first order kinetics can be applied, where the short term first order rate constants for consumption of MnO(4)(-) and oxidation of PCE are 0.05-0.5 h(-1) and 0.5-4.5 h(-1), respectively. The sediment does not act as an instantaneous sink for MnO(4)(-). The consumption of MnO(4)(-) by reaction with the reactive species in the sediment is the result of several parallel reactions, during which the reaction between the contaminant and MnO(4)(-) also takes place. Hence, application of low MnO(4)(-) concentrations can cause partly oxidation of PCE, as the oxidant demand of the sediment does not need to be met fully before PCE is oxidised.

Transport of hydrocarbons from an emplaced fuel source experiment in the vadose zone at Airbase Vaerløse, Denmark

An emplaced hydrocarbon source field experiment was conducted in the relatively homogeneous sandy geology of the vadose zone at Airbase Vaerløse, Denmark. The source (10.2 l of NAPL) consisted of 13 hydrocarbons (n-, iso- and cyclo-alkanes and aromates) and CFC-113 as a tracer. Monitoring in the 107 soil gas probes placed out to 20 m from the centre of the source showed spreading of all the compounds in the pore air and all compounds were measured in the pore air within a few hours after source emplacement. Seven of the fourteen compounds were depleted from the source within the 1 year of monitoring. The organic vapours in the pore air migrated radially from the source. The CFC-113 concentrations seemed to be higher in the deeper soil gas probes compared with the hydrocarbons, indicating a high loss of CFC-113 to the atmosphere and the lack of degradation of CFC-113. For the first days after source emplacement, the transport of CFC-113, hexane and toluene was successfully simulated using a radial gas-phase diffusion model for the unsaturated zone. Groundwater pollution caused by the vadose zone hydrocarbon vapours was only detected in the upper 30 cm of the underlying groundwater and only during the first 3 months of the experiment. Only the most water-soluble compounds were detected in the groundwater and concentrations decreased sharply with depth (approximately one order of magnitude within 10 cm depth) to non-detect at 30 cm depth. The groundwater table varied more than 1 m within the measurement period. However that did not influence the direction of the groundwater flow. Approximately 7 months after source emplacement the groundwater table rose more than 1 m within 1 month. That did not cause additional pollution of the groundwater.

Compositional evolution of the emplaced fuel source in the vadose zone field experiment at Airbase Verløse, Denmark

A field experiment was performed in a sandy vadose zone, studying the fate of an emplaced fuel-NAPL source, composed of 13 hydrocarbons and a tracer. The UNIFAC model was used to testthe nonideal behavior of the source, and the numerical model MIN3P was used for assessing the effect of biodegradation on source evolution. The diffusive loss to the surrounding vadose zone and the atmosphere created temporary gradients in mole fractions of the individual compounds within the source NAPL. The evolution of the source composition corresponded in general with expectations based on Raoult's Law, with the exception thatthe mole fractions of aromatic compounds in the source NAPL decreased faster than fractions of aliphatic compounds of similar volatility. Calculation of activity coefficients (y) using the UNIFAC model implied nonideal conditions, with composition-dependent gammas different from 1. Positive deviations were calculated for the aromatic compounds. The effect of biodegradation on source depletion, evaluated by numerical modeling, was greater for the aromatic as compared to the aliphatic compounds. Hence, the faster depletion of the aromatic relative to aliphatic compounds of similar volatility is both a result of the nonideality of the mixture and a result of partitioning and biodegradation in the pore-water. Vapor concentrations of the compounds in the source were in reasonable agreement with predictions based on the modified Raoult's Law with the UNIFAC predicted gammas and the NAPL composition for the most volatile compounds. For the less volatile compounds, the measured vapor concentrations were lower than predicted with the largest deviations for the least volatile compounds. This field experiment illustrated that nonideal behavior and bioenhanced source depletion need to be considered at multicomponent NAPL spill sites.

Degradation of the (R)- and (S)-enantiomers of the herbicides MCPP and dichlorprop in a continuous field-injection experiment

An aerobic field-injection experiment was performed to study the degradation and migration of different herbicides at trace levels in an aerobic aquifer at Vejen, Denmark. Mecoprop (MCPP) and dichlorprop monitored in a dense network of multilevel samplers were both degraded within a distance of 1 m after a period of 120 days. The study showed that no preferential degradation of the (R)- and (S)-enantiomers of MCPP and of dichlorprop took place as the enantiomeric forms of the phenoxy acids were degraded simultaneously in the aquifer.

Sorption and degradation of the herbicide 2-methyl-4,6-dinitrophenol under aerobic conditions in a sandy aquifer in Vejen, Denmark

A pulse (7 days) and a continuous (216 days), natural gradient field injection experiment with herbicides, including 2-methyl-4,6-dinitrophenol (4,6-dinitro-o-cresol, abbreviated DNOC), and a bromide tracer were conducted in a shallow, aerobic aquifer near Vejen, Denmark. The pulse and continuous plume were monitored in a dense, three-dimensional monitoring network installed in the aquifer downgradient of the injection. The sorption and degradation of DNOC were evaluated based on moment analysis of breakthrough curves, cross sections, and snapshots of the DNOC plume and supported by results from laboratory experiments. Significant and spatially variable sorption of DNOC (Kd range, 0.10-0.98 L/kg) was observed due to a specific binding of DNOC to clay minerals. The spatial variation was mainly a result of variation in pH, with stronger sorption at lower pH, whereas other factors such as cation composition on the solid matrix appeared to be negligible. Significant degradation of DNOC in the aquifer was revealed by moment analysis of data from the continuous field injection experiment. Degradation of DNOC in the field was slow and/or subject to long lag phases, and the data suggested spatially varying degradation potentials. This was supported by the laboratory experiments. The potential for natural attenuation of DNOC in aerobic aquifers appears promising.

Sorption of three heterocyclic compounds and their hydrocarbon analogue from coal-tar contaminated water in a natural clayey till

The natural attenuation of contaminants in fractured media is controlled by biodegradation and by diffusion into the matrix. The effective rate of diffusion and the mass of contaminants thereby removed from the aqueous phase in the fracture depends on the sorption of the contaminants. The sorption of three heterocyclic compounds and their hydrocarbon analogue from an aqueous solution of a complex mixture of 25 organic coal-tar compounds including all major groups of organic coal-tar compounds in clayey till was studied in batch experiments for a large range of concentrations including concentrations approaching the effective solubility of the mixture and for two different solid:liquid (s:l) ratios. For the high s:l ratio sorption experiment the isotherms were close to linear and best fitted by Freundlich isotherms. Sorption of the four compounds was correlated with Kow, but an order of magnitude greater than predicted from a Kow-Kow relationship for aquifer materials. Sorption was primarily dominated by hydrophobic sorption. However, at high surface density carbazole sorption appeared to be influenced by dipole-dipole attraction. For the low s: 1 ratio experiment at higher concentrations a dramatic increase in sorption was observed for all four compounds. The dramatic increase is likely a result of multi-layer formation/condensation of organic compounds on mineral surfaces at high surface density of organic compounds possibly induced by the complexity of solutes. The steepness of the increase and sensitivity to total compound surface density appears to be related to compound hydrophobicity. The impact of these four compounds on groundwater quality from sources of coal-tar in fractured clayey till overburden may be significantly lower than predicted from general sorption relationships.