Reactive transport modelling of biogeochemical processes and carbon isotope geochemistry inside a landfill leachate plume

Process-Based and Probabilistic Quantification of Co and Ni Mobilization Risks Induced by Managed Aquifer Recharge

Managed aquifer recharge (MAR) is an increasingly used water management technique that enhances water availability while commonly generating water quality benefits. However, MAR activities may also trigger adverse geochemical reactions, especially during the injection of oxidant-enriched waters into reducing aquifers. Where this occurs, the environmental risks and the viability of mitigating them must be well understood. Here, we develop a rigorous approach for assessing and managing the risks from MAR-induced metal mobilization. First, we develop a process-based reactive transport model to identify and quantify the main hydrogeochemical drivers that control the release of metals and their mobility. We then apply a probabilistic framework to interrogate the inherent uncertainty associated with adjustable model parameters and consider this uncertainty (i) in long-term predictions of groundwater quality changes and (ii) in scenarios that investigate the effectiveness of modifications in the water treatment process to mitigate metal release and mobility. The results suggested that Co, Ni, Zn, and Mn were comobilized during pyrite oxidation and that metal mobility was controlled (i) by the sediment pH buffering capacity and (ii) by the sorption capacity of the native aquifer sediments. Both tested mitigation strategies were shown to be effective at reducing the risk of elevated metal concentrations.

Reactive Transport: A Review of Basic Concepts with Emphasis on Biochemical Processes

Reactive transport (RT) couples bio-geo-chemical reactions and transport. RT is important to understand numerous scientific questions and solve some engineering problems. RT is highly multidisciplinary, which hinders the development of a body of knowledge shared by RT modelers and developers. The goal of this paper is to review the basic conceptual issues shared by all RT problems, so as to facilitate advancement along the current frontier: biochemical reactions. To this end, we review the basic equations to indicate that chemical systems are controlled by the set of equilibrium reactions, which are easy to model, but whose rate is controlled by mixing. Since mixing is not properly represented by the standard advection-dispersion equation (ADE), we conclude that this equation is poor for RT. This leads us to review alternative transport formulations, and the methods to solve RT problems using both the ADE and alternative equations. Since equilibrium is easy, difficulties arise for kinetic reactions, which is especially true for biochemistry, where numerous challenges are open (how to represent microbial communities, impact of genomics, effect of biofilms on flow and transport, etc.). We conclude with the basic eleven conceptual issues that we consider fundamental for any conceptually sound RT effort.

Evaluation of MNA in A Chlorinated Solvents-Contaminated Aquifer Using Reactive Transport Modeling Coupled with Isotopic Fractionation Analysis

Groundwater contamination by chlorinated hydrocarbons is a worldwide problem that poses important challenges in remediation processes. In Italy, the Legislative Decree 152/06 defines the water quality limits to be obtained during the cleanup process. In situ bioremediation techniques are becoming increasingly important due to their affordability and, under the right conditions, because they can be more effective than conventional methodologies. In the initial feasibility study phase, the numerical modeling supports the reliability of each technique. Two different codes, BIOCHLOR and PHREEQC were discussed and compared assuming different field conditions. Isotopic Fractionation-Reactive Transport Models were then developed in one synthetic and one simple field case. From the results, the two codes were in agreement and also able to demonstrate the Monitored Natural Attenuation processes occurring at the dismissed site located in Italy. Finally, the PHREEQC model was used to forecast the remediation time frame by MNA, hypothesizing a complete source cleanup: a remediation time frame of about 10–11 years was achieved by means of natural attenuation processes.

Coupling Flow, Heat, and Reactive Transport Modeling to Reproduce In Situ Redox Potential Evolution: Application to an Infiltration Pond

Redox potential (Eh) measurements are widely used as indicators of the dominant reduction-oxidation reactions occurring underground. Yet, Eh data are mostly used in qualitative terms, as actual values cannot be used to distinguish uniquely the dominant redox processes at a sampling point and should therefore be combined with a detailed geochemical characterization of water samples. In this work, we have intensively characterized the redox potential of the first meter of soil in an infiltration pond recharged with river water using a set of in situ sensors measuring every 12 min during a 1 year period. This large amount of data combined with hydrogeochemical campaigns allowed developing a reactive transport model capable of reproducing the redox potential in space and time together with the site hydrochemistry. Our results showed that redox processes were mainly driven by the amount of sedimentary organic matter in the system as well as by seasonal variation of temperature. As a subsidiary result, our work emphasizes the need to use a fully coupled model of flow, heat transport, solute transport, and the geochemical reaction network to fully reproduce the Eh observations in the topsoil.

Evaluation of geochemical processes and nitrate pollution sources at the Ljubljansko polje aquifer (Slovenia): A stable isotope perspective

The Ljubljansko polje aquifer, which is the main supply of drinking water for the local population in Ljubljana, Slovenia is highly vulnerable to anthropogenic pollution. In this study, the geochemistry of major constituents including nitrate concentrations and the dual isotopes of nitrate were used to ascertain the spatial distribution of processes and nitrate sources in the groundwater from seven wells at three different water supplies: Kleče, Hrastje and Jarški prod. The groundwater is of the Ca²⁺-Mg²⁺-HCO₃^- type approaching equilibrium with respect to dolomite and are moderately supersaturated with calcite. The groundwater nitrate concentrations ranged from 5.32 to 50.1 mg L^-1 and are well above the threshold value for anthropogenic activity (3 mg L^-1). The δ¹⁵N_NO3 values ranged from 1.4 to 9.7‰, while δ¹⁸O_NO3 values were from 6.3 to 34.6‰. Based on isotope mixing model three sources of nitrate were identified: atmospheric deposition, fertilizers and soil nitrogen. At Kleče 8, Kleče 12 and Jarški prod 3 the low δ¹⁵N_NO3 and high δ¹⁸O_NO3 values result from pristine nitrate sources, while in Hrastje 3 and Kleče 11 equal amounts of nitrate derived from soils with mixed fertilization and sewage. The data also indicate that the main sources of high nitrate concentrations in groundwater are from fertilizers and sewage-manure (comprising up to 64%). Such levels occurred in the Hrastje and Kleče 11 wells where precipitation is the main source of groundwater. Nitrate derived from atmospheric deposition accounted for approximately 10% of the nitrate in the groundwater. The message from this study is that to reduce the nitrogen load and improve water quality will involve containment and the careful management of sources from urban and agriculture inputs such as sewage-manure and fertilizers.

Trace Element Behavior in Methane-Rich and Methane-Free Groundwater in North and East Texas

There is concern about adverse impacts of natural gas (primarily methane) production on groundwater quality; however, data on trace element concentrations are limited. The objective of this study was to compare the distribution of trace elements in groundwater samples with and without dissolved methane in aquifers overlying the Barnett Shale (Hood and Parker counties, 207 samples) and the Haynesville Shale (Panola County, 42 samples). Both shales have been subjected to intensive hydraulic fracturing for gas production. Well clusters with high dissolved methane were previously found in these counties and are thought to be of natural origin. Overall, groundwater in these counties is of excellent quality with typically low elemental concentrations. Several statistical analyses strongly suggest that most trace element concentrations, generally at low background levels, are no higher and even reduced when dissolved methane is present. In addition, trace element concentrations are not correlated with distance to gas wells. The reduction in trace element concentrations is attributed to anaerobic microbial degradation of methane, is associated with a higher pH (>8.5), and, likely, with precipitation of carbonates and pyrite and formation of clays. Trace and other elements are likely incorporated within the precipitating mineral crystalline network or sorbed. High pH values are found throughout these high-methane clusters (e.g., Parker-Hood cluster), even in subregions where methane is not present, which is consistent with a pervasive natural origin of dissolved methane rather than a limited gas well source.

Applying the Rayleigh Approach for Stable Isotope-Based Analysis of VOC Biodegradation in Diffusion-Dominated Systems

Compound-specific stable isotope analysis (CSIA) has become an established tool for assessing biodegradation in the subsurface. Diffusion-dominated vapor phase transport thereby is often excluded from quantitative assessments due to the problem of diffusive mixing of concentrations with different isotopic signatures for CSIA interpretation. In soils and other unsaturated porous media volatile organic compounds (VOCs) however, are mainly transported via gas-phase diffusion and may thus prohibit a CSIA-based quantitative assessment of the fate of VOCs. The present study presents and verifies a concept for the assessment of biodegradation-induced stable isotope fractionation along a diffusive transport path of VOCs in unsaturated porous media. For this purpose data from batch and column toluene biodegradation experiments in unsaturated porous media were combined with numerical reactive transport simulations; both addressing changes of concentration and stable isotope fractionation of toluene. The numerical simulations are in good agreement with the experiment data, and our results show that the presented analytically derived assessment concept allows using the slope of the Rayleigh plot to obtain reasonable estimates of effective in situ fractionation factors in spite of diffusion-dominated transport. This enlarges the application range of CSIA and provides a mean for a better understanding of VOC fate in the unsaturated subsurface.

Environmental monitoring of a landfill area through the application of carbon stable isotopes, chemical parameters and multivariate analysis

Leachate produced during an organic matter decomposition process has a complex composition and can cause contamination of surface and groundwaters adjacent to a landfill area. The monitoring of these areas is extremely important for the characterization of the leachate produced and to avoid or mitigate environmental damages. Thus, the present study has the objective of monitoring the area of a Brazilian landfill using conventional parameters (dissolved metals and anions in water) and alternative, stable carbon isotopes parameters (δ¹³C of dissolved organic and inorganic carbons in water) in addition to multivariate analysis techniques. The use of conventional and alternative parameters together with multivariate analysis showed that cells of the residues are at different phases of stabilization of the organic matter and probably already at C3 of the methanogenic phase of decomposition. In addition, the data showed that organic matter stabilization ponds present in the landfill are efficient and improve the quality of the leachate. Enrichment of the heavy ¹³C isotope in both surface and groundwater suggested contamination in two sampling sites.

Impact of an historic underground gas well blowout on the current methane chemistry in a shallow groundwater system

The rapid increase in shale gas production in recent years has led to increased attention to its potential negative environmental effects, including the risks of contaminating groundwater with methane and other substances. In this context, the uncontrolled gas migration that is triggered during well blowouts is an understudied environmental hazard. We show that the methane chemistry in shallow groundwater overlying the site of a catastrophic underground blowout continues to be impacted 50 y later. The occurrence of anaerobic methane oxidation limits the spatial extent to which the dissolved thermogenic methane plume could be observed and discerned from local biogenic methane sources. However, it also highlights the requirement to carry out monitoring in close proximity to potential gas leakage sources.

Blowouts present a small but genuine risk when drilling into the deep subsurface and can have an immediate and significant impact on the surrounding environment. Nevertheless, studies that document their long-term impact are scarce. In 1965, a catastrophic underground blowout occurred during the drilling of a gas well in The Netherlands, which led to the uncontrolled release of large amounts of natural gas from the reservoir to the surface. In this study, the remaining impact on methane chemistry in the overlying aquifers was investigated. Methane concentrations higher than 10 mg/L (n = 12) were all found to have δ¹³C-CH₄ values larger than −30‰, typical of a thermogenic origin. Both δ¹³C-CH₄ and δD-CH₄ correspond to the isotopic composition of the gas reservoir. Based on analysis of local groundwater flow conditions, this methane is not a remnant but most likely the result of ongoing leakage from the reservoir as a result of the blowout. Progressive enrichment of both δ¹³C-CH₄ and δD-CH₄ is observed with increasing distance and decreasing methane concentrations. The calculated isotopic fractionation factors of ε_C = 3 and ε_D = 54 suggest anaerobic methane oxidation is partly responsible for the observed decrease in concentrations. Elevated dissolved iron and manganese concentrations at the fringe of the methane plume show that oxidation is primarily mediated by the reduction of iron and manganese oxides. Combined, the data reveal the long-term impact that underground gas well blowouts may have on groundwater chemistry, as well as the important role of anaerobic oxidation in controlling the fate of dissolved methane.

Identification of groundwater redox process induced by landfill leachate based on sensitive factor method

Landfill site is a significant source of groundwater pollution. To ensure that the groundwater contamination of landfills can be controlled and repaired scientifically, the identification of groundwater pollution process is needed. On the basis of biogeochemical process of leachate pollutants in the groundwater environment, a sensitive factor method for the identification of groundwater redox process from landfills was established in this research. The method encompasses four phases, including sensitive factors selection, redox zone characterization, weight calculation, and redox zone identification. In the sensitive factor index system employed here, five indicators involving dissolved oxygen (DO), nitrite, Fe²⁺, sulfide, and CO₂ were selected. The boundary of each redox zones was determined by the quantitative method, and the weight of each indicator was calculated by combined weight method. This method was applied to a landfill site in the northeast of China. The result showed that there were five redox zones that appeared in pollution plume, including methanogenic zone (MGZ), sulfate reduction zone (SRZ), iron reduction zone (IRZ), nitrate reduction zone (NRZ), and oxygen reduction zone (ORZ). The results were consistent with the actual situation of the site. The sensitive factor method was scientific and effective to identify the groundwater redox process in landfill and can provide reference data related to investigation and remediation of groundwater pollution in landfill sites.

Combination of a crop model and a geochemical model as a new approach to evaluate the sustainability of an intensive agriculture system

By combining a crop model (STICS) and a geochemical model (PHREEQC), a new approach to assess the sustainability of agrosystems is proposed. It is based upon aqueous geochemistry and the stepwise modifications of soil solution during its transfer from the surface till aquifer. Meadows of Crau (SE France), irrigated since the 16th century, were field monitored (2012-2015) and modelled. Except for N, the mineral requirements of hay are largely covered by dissolved elements brought by irrigation water with only slight deficits in K and P, which are compensated by P-K fertilizers and the winter pasture by sheep. N cycle results in a very small nitrate leakage. The main determinants of the chemical composition changes of water are: concentration by evaporation, equilibration with soil pCO₂, mineral nutrition of plants, input of fertilizers, sheep grazing, mineral-solution interactions in superficial formations till the aquifer, including ion exchange. Inverse modelling with PHREEQC allows for quantifying these processes. For groundwater, measured composition fit statistically very well with those computed, validating thus this approach. This long-term established agrosystem protects both soil and water resources: soil nutritional status remains constant with even some P and (minor) K fixation in soils; long-term decarbonatation occurs but it is greatly slowed by saturation of irrigation water by carbonate; P fixation in soil protects groundwater from eutrophication.

A preliminary report of indigenous fungal isolates from contaminated municipal solid waste site in India

Municipal solid waste (MSW) containing harmful substances is a major concern in waste management and can cause adverse effects on diversity of fungi in soil. The main objective was to evaluate the fungal diversity inhabiting in the soil nearby MSW disposal site. The fungal strains were isolated in potato dextrose agar (PDA), media at temperatures 28 ± 1 °C by using standard serial dilution pour plate method, and appeared fungal colonies identified based on morphological characteristics. The overall most fungal diversity was found in soil sample collected from S5, followed by S4, S3, S1, and least in S2 site. A total of 24 fungal isolates recovered from the different MSW sites and Aspergillus sp., Fusarium sp., and Curvularia sp. genus has isolated from all the samples. In addition, the metal tolerance index performed because it needs to classify the fungus for their best use as potential agent for environmental protection. The metal tolerance outcomes revealed that both metals (cadmium and chromium) has appeared as the highest growth inhibitor for most strains and even fungal colonies did not propagate very well on the surface of media. Therefore, these findings suggest that the pre-adapted indigenous fungal isolates have proven remarkable tolerance ability to both metals. Furthermore, these highly metal-tolerant fungal strains are recommended for detail research or can use in pilot-scale bioremediation application to treat contaminated site.

Dual Mechanism Conceptual Model for Cr Isotope Fractionation during Reduction by Zerovalent Iron under Saturated Flow Conditions

Chromium isotope analysis is rapidly becoming a valuable complementary tool for tracking Cr(VI) treatment in groundwater. Evaluation of various treatment materials has demonstrated that the degree of isotope fractionation is a function of the reaction mechanism, where reduction of Cr(VI) to Cr(III) induces the largest fractionation. However, it has also been observed that uniform flow conditions can contribute complexity to isotope measurements. Here, laboratory batch and column experiments were conducted to assess Cr isotope fractionation during Cr(VI) reduction by zerovalent iron under both static and saturated flow conditions. Isotope measurements were accompanied by traditional aqueous geochemical measurements (pH, Eh, concentrations) and solid-phase analysis by scanning electron microscopy and X-ray absorption spectroscopy. Increasing δ(53)Cr values were associated with decreasing Cr(VI) concentrations, which indicates reduction; solid-phase analysis showed an accumulation of Cr(III) on the iron. Reactive transport modeling implemented a dual mechanism approach to simulate the fractionation observed in the experiments. The faster heterogeneous reaction pathway was associated with minimal fractionation (ε=-0.2‰), while the slower homogeneous pathway exhibited a greater degree of fractionation (ε=-0.9‰ for the batch experiment, and ε=-1.5‰ for the column experiment).

Determination of the origin of dissolved inorganic carbon in groundwater around a reclaimed landfill in Otwock using stable carbon isotopes

Chemical and isotopic analyses of groundwater from piezometers located around a reclaimed landfill in Otwock (Poland) were performed in order to trace the origin of dissolved inorganic carbon (DIC) in the groundwater. Due to differences in the isotopic composition of carbon from different sources, an analysis of stable carbon isotopes in the groundwater, together with the Keeling plot approach and a two-component mixing model allow us to evaluate the relative contributions of carbon from these sources in the groundwater. In the natural (background) groundwater, DIC concentrations and the isotopic composition of DIC (δ(13)CDIC) comes from two sources: decomposition of organic matter and carbonate dissolution within the aquifer sediments, whereas in the leachate-contaminated groundwater, DIC concentrations and δ(13)CDIC values depend on the degradation of organic matter within the aquifer sediments and biodegradation of organic matter stored in the landfill. From the mixing model, about 4-54% of the DIC pool is derived from organic matter degradation and 96-46% from carbonate dissolution in natural conditions. In the leachate-contaminated groundwater, about 20-53% of the DIC is derived from organic matter degradation of natural origin and 80-47% from biodegradation of organic matter stored in the landfill. Partial pressure of CO2 (P CO2) was generally above the atmospheric, hence atmospheric CO2 as a source of carbon in DIC pool was negligible in the aquifer. P CO2 values in the aquifer in Otwock were always one to two orders of magnitude above the atmospheric P CO2, and thus CO2 escaped directly into the vadose zone.

Quantifying chemical reactions by using mixing analysis

This work is motivated by a sound understanding of the chemical processes that affect the organic pollutants in an urban aquifer. We propose an approach to quantify such processes using mixing calculations. The methodology consists of the following steps: (1) identification of the recharge sources (end-members) and selection of the species (conservative and non-conservative) to be used, (2) identification of the chemical processes and (3) evaluation of mixing ratios including the chemical processes. This methodology has been applied in the Besòs River Delta (NE Barcelona, Spain), where the River Besòs is the main aquifer recharge source. A total number of 51 groundwater samples were collected from July 2007 to May 2010 during four field campaigns. Three river end-members were necessary to explain the temporal variability of the River Besòs: one river end-member is from the wet periods (W1) and two are from dry periods (D1 and D2). This methodology has proved to be useful not only to compute the mixing ratios but also to quantify processes such as calcite and magnesite dissolution, aerobic respiration and denitrification undergone at each observation point.

Reactive modelling of 1,2-DCA and DOC near the shoreline

1,2-Dichloroethane (1,2-DCA) was found to be the most abundant compound among chlorinated hydrocarbons detected in a petrochemical plant in southern Italy. This site is located near the coastline, and it is set above an unconfined coastal aquifer, where seawater intrusion is present. The presence of organic and inorganic contaminants at this site has required the implementation of remediation strategies, consisting of pumping wells (hydraulic barrier) and a horizontal flow barrier. The purpose of this work was to assess the influence of salt water intrusion on the degradation rate of 1,2-DCA. This was done on a three-dimensional domain relative to a limited portion of a well characterized field site, accounting for density-dependent flow and reactive transport modelling of 1,2-DCA and Dissolved Organic Carbon (DOC). The modelling procedure was performed employing SEAWAT-4.0 and PHT3D, to reproduce the complex three-dimensional flow and transport domain. In order to determine the fate of 1,2-DCA, detailed field investigations provided intensive depth profile information. Different, kinetically controlled degradation rates were simulated to explain the observed, selective degradation of pollutants in groundwater. Calibration of the model was accomplished by comparison with the two different sets of measurements obtained from the MLS devices and from pumping wells. With the calibrated model, it was possible to distinguish between dispersive non-reactive processes and bacterially mediated reactions. In the non-reactive model, 1,2-DCA sorption was simulated using linear sorption coefficient determined with field data and 1,2-DCA degradation was simulated using a first order decay coefficient using literature data as initial guess. Finally, on the reactive transport model, where a two-step approach with partial equilibrium approach was implemented, the effects of neglecting the cation exchange capacity, omitting density-dependent flow, and refining the vertical discretization of the model were investigated. Comparison of results from various scenarios shows that geochemical changes in inorganic constituents can be used to improve the site's conceptual model, and establishes that natural degradation processes can be suitable for 1,2-DCA as a remediation option.

Stochastic analytical modeling of the biodegradation of steady plumes

We present a stochastic analytical framework to assess the contaminant concentration of a steady plume undergoing biodegradation. The method is focused on heterogeneous formations, and it embeds both fringe and core degradation. The Lagrangian concentration approach of Fiori (2001) was employed, which is suited for describing the interplay between the large scale advection caused by heterogeneity and the local dispersion processes. The principal scope of the model is to provide a relatively simple tool for a quick assessment of the contamination level in aquifers, as function of a few relevant, physically based dimensionless parameters. The solution of the analytical model is relatively simple and generalizes previous approaches developed for homogeneous formations. It is found that heterogeneity generally enhances mixing and degradation; in fact, the plume shear and distortion operated by the complex, heterogeneous velocity field facilitates local dispersion in diluting the contaminant and mixing it with the electron acceptor. The decay of the electron donor concentration, and so the plume length, is proportional to the transverse pore-scale dispersivity, which is indeed the parameter ruling mixing and hence degradation. While the theoretical plume length is controlled by the fringe processes, the core degradation may determine a significant decay of concentration along the mean flow direction, thus affecting the length of the plume. The method is applied to the crude oil contamination event at the Bemijdi site, Minnesota (USA).

Reactive transport modeling of chromium isotope fractionation during Cr(VI) reduction

Chromium isotope fractionation is indicative of mass-transfer processes, such as reduction of Cr(VI) to Cr(III) during groundwater remediation. Laboratory experiments comparing batch and column treatment of Cr(VI) using organic carbon suggest that the associated isotope fractionation may be influenced by solute-transport mechanisms. These batch and column experiments were simulated using the reactive transport model MIN3P to further evaluate the effects of Cr reduction and transport on isotope fractionation under saturated flow conditions. Simulation of the batch experiment provided a good fit to the experimental data, where a fractionation factor (α₅₃) of 0.9965 was attributed to a single, dominant Cr(VI) removal mechanism. Calibration of the column simulations to the experimental results suggested the presence of a second, more rapid Cr(VI) removal mechanism with α₅₃ = 0.9992. Results from this study demonstrate that the interpretation of Cr isotope fractionation during reduction can be complex, particularly where multiple removal mechanisms are evident. Reactive transport modeling of Cr isotope fractionation can provide a quantitative assessment of the contaminant removal mechanisms, thus improving the application of Cr isotope measurements as a tool to track Cr(VI) migration and attenuation in groundwater.

Transverse hydrodynamic dispersion effects on isotope signals in groundwater chlorinated solvents' plumes

The effects of transverse hydrodynamic dispersion on altering transformation-induced compound-specific isotope analysis (CSIA) signals within groundwater pollution plumes have been assessed with reactive transport modeling accommodating diffusion-induced isotope fractionation (DIF) and implementing different parameterizations of local transverse dispersion. The model reproduced previously published field data showing a negative carbon isotope pattern (-2 ‰) at the fringes of a nondegrading PCE plume. We extended the study to reactive transport scenarios considering vinyl chloride as a model compound and assessing, through a detailed sensitivity analysis, the coupled effects of transverse hydrodynamic dispersion (with and without DIF) and aerobic fringe degradation on the evolution of carbon and chloride isotope ratios. Transformation-induced positive isotope signals were increasingly attenuated with distance from the source and higher degradation rate. The effect of DIF on the overall isotope signal attenuation was greatest near the source and for low values of groundwater flow velocity, transverse dispersion coefficient, molecular weight, rate constant, and isotope fractionation factor, α, of the degradation reaction. Models disregarding DIF underestimate the actual α. The approximately twice larger DIF effect for chlorine than for carbon together with the low α for oxidation resulted in strong chlorine CSIA depletions for VC at the plume fringe.

Pollution profiles and physicochemical parameters in old uncontrolled landfills

The long-term effectiveness of the geological barrier beneath municipal-waste landfills is a critical issue for soil and groundwater protection. This study examines natural clayey soils directly in contact with the waste deposited in three landfills over 12 years old in Spain. Several physicochemical and geological parameters were measured as a function of depth. Electrical conductivity (EC), water-soluble organic carbon (WSOC), Cl(-), NH(4)(+), Na(+) and exchangeable NH(4)(+) and Na(+) were used as parameters to measure the penetration of landfill leachate pollution. Mineralogy, specific surface area and cationic-exchange capacities were analyzed to characterize the materials under the landfills. A principal component analysis, combined with a Varimax rotation, was applied to the data to determine patterns of association between samples and variables not evident upon initial inspection. The main factors explaining the variation in the data are related to waste composition and local geology. Although leachates have been in contact with clays for long time periods (13-24 years), WSOC and EC fronts are attenuated at depths of 0.2-1.5m within the clay layer. Taking into account this depth of the clayey materials, these natural substrata (>45% illite-smectite-type sheet silicates) are suitable for confining leachate pollution and for complying with European legislation. This paper outlines the relevant differences in the clayey materials of the three landfills in which a diffusive flux attenuation capacity (A(c)) is defined as a function (1) of the rate of decrease of the parameters per meter of material, (2) of the age and area of the landfill and (3) of the quantity and quality of the wastes.

Biogeochemical evolution of a landfill leachate plume, Norman, Oklahoma

Leachate from municipal landfills can create groundwater contaminant plumes that may last for decades to centuries. The fate of reactive contaminants in leachate-affected aquifers depends on the sustainability of biogeochemical processes affecting contaminant transport. Temporal variations in the configuration of redox zones downgradient from the Norman Landfill were studied for more than a decade. The leachate plume contained elevated concentrations of nonvolatile dissolved organic carbon (NVDOC) (up to 300 mg/L), methane (16 mg/L), ammonium (650 mg/L as N), iron (23 mg/L), chloride (1030 mg/L), and bicarbonate (4270 mg/L). Chemical and isotopic investigations along a 2D plume transect revealed consumption of solid and aqueous electron acceptors in the aquifer, depleting the natural attenuation capacity. Despite the relative recalcitrance of NVDOC to biodegradation, the center of the plume was depleted in sulfate, which reduces the long-term oxidation capacity of the leachate-affected aquifer. Ammonium and methane were attenuated in the aquifer relative to chloride by different processes: ammonium transport was retarded mainly by physical interaction with aquifer solids, whereas the methane plume was truncated largely by oxidation. Studies near plume boundaries revealed temporal variability in constituent concentrations related in part to hydrologic changes at various time scales. The upper boundary of the plume was a particularly active location where redox reactions responded to recharge events and seasonal water-table fluctuations. Accurately describing the biogeochemical processes that affect the transport of contaminants in this landfill-leachate-affected aquifer required understanding the aquifer's geologic and hydrodynamic framework.

Metal speciation in landfill leachates with a focus on the influence of organic matter

This study characterises the heavy-metal content in leachates collected from eight landfills in France. In order to identify heavy metal occurrence in the different size fractions of leachates, a cascade filtration protocol was applied directly in the field, under a nitrogen gas atmosphere to avoid metal oxidation. The results of analyses performed on the leachates suggest that most of the metals are concentrated in the <30 kDa fraction, while lead, copper and cadmium show an association with larger particles. Initial speciation calculations, without considering metal association with organic matter, suggest that leachate concentrations in lead, copper, nickel and zinc are super-saturated with respect to sulphur phases. Speciation calculations that account for metal complexation with organic matter, considered as fulvic acids based on C1(s) NEXAFS spectroscopy, show that this mechanism is not sufficient to explain such deviation from equilibrium conditions. It is therefore hypothesized that the deviation results also from the influence of biological activity on the kinetics of mineral phase precipitation and dissolution, thus providing a dynamic system. The results of chemical analyses of sampled fluids are compared with speciation calculations and some implications for the assessment of metal mobility and natural attenuation in a context of landfill risk assessment are discussed.

34S/32S fractionation during sulfate reduction in groundwater treatment systems: reactive transport modeling

Isotope ratio measurements provide a tool for indicating the relative significance of biogeochemical reactions and for constraining estimates of the extent and rate of reactions in passive treatment systems. In this paper, the reactive transport model MIN3P is used to evaluate sulfur isotope fractionation in column experiments designed to simulate treatment of contaminated water by microbially mediated sulfate reduction occurring within organic carbon-based and iron and carbon-based permeable reactive barriers. A mass dependent fractionation model was used to determine reaction rates for 32S and 34S compounds during reduction, precipitation, and dissolution reactions and to track isotope-dependent mass transfer during SO4 removal. The δ34S values obtained from the MIN3P model were similar to those obtained from the Rayleigh equation, indicating that there was not a significant difference between the conceptual models. Differences between the MIN3P derived α value and the Rayleigh equation derived value were attributed to minor changes in the dissolution and precipitation rate of gypsum and mathematical differences in the fitting models. The results indicated that the prediction of δ34S was fairly insensitive to differences in the fractionation factor at the concentration ranges measured in the current study. However, more significant differences would be expected at low sulfate conditions.

Analytical modelling of fringe and core biodegradation in groundwater plumes

Biodegradation can be divided into two categories depending on the location at which it occurs within the plume: degradation at the plume fringes, and degradation in the interior (core). Available analytical solutions are limited to the consideration of either fringe or core degradation, which in turn limits the applicability of these solutions. Here, a new analytical approach to modelling plumes with both fringe and core degradation is presented. The approach relies on the use of readily available analytical solutions for solute transport. Using a well-known solution for three-dimensional solute transport from a planar source, an approximate solution is derived for the maximum plume length at steady-state conditions. This is verified through the use of a numerical solution. The solution suggests that the parameters controlling the plume length are: (i) the size of the contaminant source, (ii) electron acceptor to electron donor ratio, (iii) transverse dispersivities and (iv) the ratio between degradation rate constant and velocity (lambda/v). The latter term provides a simple check on the relative weights of transport to core degradation and can be used to estimate the importance of core degradation in the overall plume attenuation. The well-documented Bemidji field site has both fringe and core degradation. The new combined degradation model estimates the length of the plume with 10 m of the observed length; core only and fringe only solutions overestimate the length by more than a factor of 2.

Biogeochemistry at a wetland sediment-alluvial aquifer interface in a landfill leachate plume

The biogeochemistry at the interface between sediments in a seasonally ponded wetland (slough) and an alluvial aquifer contaminated with landfill leachate was investigated to evaluate factors that can effect natural attenuation of landfill leachate contaminants in areas of groundwater/surface-water interaction. The biogeochemistry at the wetland-alluvial aquifer interface differed greatly between dry and wet conditions. During dry conditions (low water table), vertically upward discharge was focused at the center of the slough from the fringe of a landfill-derived ammonium plume in the underlying aquifer, resulting in transport of relatively low concentrations of ammonium to the slough sediments with dilution and dispersion as the primary attenuation mechanism. In contrast, during wet conditions (high water table), leachate-contaminated groundwater discharged upward near the upgradient slough bank, where ammonium concentrations in the aquifer where high. Relatively high concentrations of ammonium and other leachate constituents also were transported laterally through the slough porewater to the downgradient bank in wet conditions. Concentrations of the leachate-associated constituents chloride, ammonium, non-volatile dissolved organic carbon, alkalinity, and ferrous iron more than doubled in the slough porewater on the upgradient bank during wet conditions. Chloride, non-volatile dissolved organic carbon (DOC), and bicarbonate acted conservatively during lateral transport in the aquifer and slough porewater, whereas ammonium and potassium were strongly attenuated. Nitrogen isotope variations in ammonium and the distribution of ammonium compared to other cations indicated that sorption was the primary attenuation mechanism for ammonium during lateral transport in the aquifer and the slough porewater. Ammonium attenuation was less efficient, however, in the slough porewater than in the aquifer and possibly occurred by a different sorption mechanism. A stoichiometrically balanced increase in magnesium concentration with decreasing ammonium and potassium concentrations indicated that cation exchange was the sorption mechanism in the slough porewater. Only a partial mass balance could be determined for cations exchanged for ammonium and potassium in the aquifer, indicating that some irreversible sorption may be occurring. Although wetlands commonly are expected to decrease fluxes of contaminants in riparian environments, enhanced attenuation of the leachate contaminants in the slough sediment porewater compared to the aquifer was not observed in this study. The lack of enhanced attenuation can be attributed to the fact that the anoxic plume, comprised largely of recalcitrant DOC and reduced inorganic constituents, interacted with anoxic slough sediments and porewaters, rather than encountering a change in redox conditions that could cause transformation reactions. Nevertheless, the attenuation processes in the narrow zone of groundwater/surface-water interaction were effective in reducing ammonium concentrations by a factor of about 3 during lateral transport across the slough and by a factor of 2 to 10 before release to the surface water. Slough porewater geochemistry also indicated that the slough could be a source of sulfate in dry conditions, potentially providing a terminal electron acceptor for natural attenuation of organic compounds in the leachate plume.

Eukaryotic diversity in an anaerobic aquifer polluted with landfill leachate

Eukaryotes may influence pollutant degradation processes in groundwater ecosystems by activities such as predation on bacteria and recycling of nutrients. Culture-independent community profiling and phylogenetic analysis of 18S rRNA gene fragments, as well as culturing, were employed to obtain insight into the sediment-associated eukaryotic community composition in an anaerobic sandy aquifer polluted with landfill leachate (Banisveld, The Netherlands). The microeukaryotic community at a depth of 1 to 5 m below the surface along a transect downgradient (21 to 68 m) from the landfill and at a clean reference location was diverse. Fungal sequences dominated most clone libraries. The fungal diversity was high, and most sequences were sequences of yeasts of the Basidiomycota. Sequences of green algae (Chlorophyta) were detected in parts of the aquifer close (<30 m) to the landfill. The bacterium-predating nanoflagellate Heteromita globosa (Cercozoa) was retrieved in enrichments, and its sequences dominated the clone library derived from the polluted aquifer at a depth of 5 m at a location 21 m downgradient from the landfill. The number of culturable eukaryotes ranged from 10(2) to 10(3) cells/g sediment. Culture-independent quantification revealed slightly higher numbers. Groundwater mesofauna was not detected. We concluded that the food chain in this polluted aquifer is short and consists of prokaryotes and fungi as decomposers of organic matter and protists as primary consumers of the prokaryotes.

Beyond the Rayleigh equation: reactive transport modeling of isotope fractionation effects to improve quantification of biodegradation

The Rayleigh equation is commonly applied to evaluate the extent of degradation at contaminated sites for which compound-specific isotope analysis (CSIA) data are available. However, it was shown recently that (i) the Rayleigh equation systematically underestimates the extent of biodegradation in physically heterogeneous systems, while (ii) it overestimates biodegradation if sorption-based carbon isotope fractionation is relevant. This paper further explores these two isotope effects not captured by the Rayleigh equation by means of a numerical modeling approach. The reactive multicomponent transport simulations show that the systematic underestimation is considerably larger for fringe-controlled and Monod-type degradation reactions than for previously assumed redox-insensitive first-order degradation kinetics, while forthe nonsteady state front portion of plumes, the Rayleigh equation may falsely indicate the occurrence of and/or overestimate biodegradation. The latter anomaly results from carbon isotope fractionation during sorption. It occurs for both supply-controlled degradation at the plume fringe and slow, reaction-controlled degradation inside the plume core. The numerical model approach enables a more accurate interpretation of CSIA data and thereby improves the quantification of biodegradation processes.

One-dimensional model for biogeochemical interactions and permeability reduction in soils during leachate permeation

This paper uses the findings from a column study to develop a reactive model for exploring the interactions occurring in leachate-contaminated soils. The changes occurring in the concentrations of acetic acid, sulphate, suspended and attached biomass, Fe(II), Mn(II), calcium, carbonate ions, and pH in the column are assessed. The mathematical model considers geochemical equilibrium, kinetic biodegradation, precipitation-dissolution reactions, bacterial and substrate transport, and permeability reduction arising from bacterial growth and gas production. A two-step sequential operator splitting method is used to solve the coupled transport and biogeochemical reaction equations. The model gives satisfactory fits to experimental data and the simulations show that the transport of metals in soil is controlled by multiple competing biotic and abiotic reactions. These findings suggest that bioaccumulation and gas formation, compared to chemical precipitation, have a larger influence on hydraulic conductivity reduction.

Modelling assessment of regional groundwater contamination due to historic smelter emissions of heavy metals

Historic emissions from ore smelters typically cause regional soil contamination. We developed a modelling approach to assess the impact of such contamination on groundwater and surface water load, coupling unsaturated zone leaching modelling with 3D groundwater transport modelling. Both historic and predictive modelling were performed, using a mass balance approach for three different catchments in the vicinity of three smelters. The catchments differ in their hydrology and geochemistry. The historic modelling results indicate that leaching to groundwater is spatially very heterogeneous due to variation in soil characteristics, in particular soil pH. In the saturated zone, cadmium is becoming strongly retarded due to strong sorption at neutral pH, even though the reactivity of the sandy sediments is low. A comparison between two datasets (from 1990 to 2002) on shallow groundwater and modelled concentrations provided a useful verification on the level of statistics of "homogeneous areas" (areas with comparable land use, soil type and geohydrological situation) instead of comparison at individual locations. While at individual locations observations and the model varies up to two orders of magnitude, for homogeneous areas, medians and ranges of measured concentrations and the model results are similar. A sensitivity analysis on metal input loads, groundwater composition and sediment geochemistry reveals that the best available information scenario based on the median value of input parameters for the model predicts the range in observed concentrations very well. However, the model results are sensitive to the sediment contents of the reactive components (organic matter, clay minerals and iron oxides). Uncertainty in metal input loads and groundwater chemistry are of lesser importance. Predictive modelling reveals a remarkable difference in geochemical and hydrological controls on subsurface metal transport at catchment-scale. Whether the surface water load will peak within a few decades or continue to increase until after 2050 depends on the dominant land use functions in the areas, their hydrology and geochemical build-up.

The opposing effects of bacterial activity and gas production on anaerobic TCE degradation in soil columns

This laboratory study explores the effect of growth substrate concentration on the anaerobic degradation of trichloroethylene (TCE) in sand packed columns. In all columns the growth substrate rapidly degraded to gas, that formed a separate phase. Biomass accumulated in the 0-4.8 cm section of the columns in proportion to the influent growth substrate concentration and biomass concentrations in the remaining sections of all columns were similar to the column receiving the lowest substrate concentration. Increases in growth substrate concentration up to 3030 mg-CODl(-1) promoted TCE degradation, but a further increase to 14300 mg-CODl(-1) reduced the amount of TCE completely dechlorinated but did not affect the production of chlorinated TCE intermediates. The mathematical model developed here satisfactorily described the enhancement in TCE dehalogenation for substrate concentration up to 3030 mg-CODl(-1); reproducing TCE dehalogenation for 14300 mg-CODl(-1) required that the moisture content used in simulation be lowered to 0.1. The study shows that volatilization of TCE can be significant and volatilization losses should be taken into account when anaerobic activity in in-situ bioremediation applications is stimulated via addition of growth substrates. An implication of the modeling simulations is that maintaining a lower, but uniform, substrate concentration over the contaminated region may lead to faster contaminant degradation.

Modelling of iron cycling and its impact on the electron balance at a petroleum hydrocarbon contaminated site in Hnevice, Czech Republic

Over a period of several decades multiple leaks of large volumes from storage facilities located near Hnevice (Czech Republic) have caused the underlying Quaternary aquifer to be severely contaminated with nonaqueous phase liquid (NAPL) petroleum hydrocarbons. Beginning in the late 1980's the NAPL plume started to shrink as a consequence of NAPL dissolution exceeding replenishment and due to active remediation. The subsurface was classified geochemically into four different zones, (i) a contaminant-free zone never occupied by NAPL or dissolved contaminants, (ii) a re-oxidation zone formerly occupied by NAPL, (iii) a zone currently occupied by NAPL, and (iv) a lower fringe zone between the overlying NAPL and the deeper underlying contaminant-free zone. The study investigated the spatial and temporal variability of the redox zonation at the Hnevice site and quantified the influence of iron-cycling on the overall electron balance. As a first step inverse geochemical modelling was carried out to identify possible reaction models and mass transfer processes. In a subsequent step, two-dimensional (forward) multi-component reactive transport modelling was performed to evaluate and quantify the major processes that control the geochemical evolution at the site. The study explains the observed enrichment of the lower fringe zone with ferrihydrite as a result of the re-oxidation of ferrous iron. It suggests that once the NAPL zone started to shrink the dissolution of previously formed siderite and FeS by oxygen and nitrate consumed a significant part of the oxidation capacity for a considerable time period and therefore limited the penetration of electron acceptors into the NAPL contaminated zone.

A knowledge-based approach to environmental biomonitoring

This paper presents the design, development and implementation of an integrated GIS-controlled knowledge-based system for environmental monitoring applications, utilizing indigenous flora for assessing quality. The system gathers and combines geographical, ecological, and physicochemical data of organisms' response to pollution within an intelligent computer program that (a) recognises groups of indigenous species suitable for long-term monitoring of a specific pollutant or a combination of pollutants, (b) estimates the ambient concentration of pollutant(s) from the population of the species comprising the bioindicator group and (c) provides biomonitoring capacity indices at national and international/transboundary levels. Significantly, a novel system in the form of a rational framework at the conceptual design level has been developed, that actually contributes towards achieving a cost-effective long-term biomonitoring program, with the flexibility to counter on-course any (anticipated or not) variations/modifications of the surveillance environment: the scheme assumes a robust dynamic cooperation between instrumental and biomonitoring systems, with a view to minimise uncertainty and monitoring costs and increase reliability of pollution control and abatement, aiming eventually at the shifting, partially or totally, from instrumental to natural monitoring. The proposed approach is presently implemented at pilot-scale for establishing a biomonitoring network at a large industrial area in Greece. The results obtained indicate that a cost-effective program can be only attained and maintained under a suitable financial/organizational scheme at the macro level, whereas the micro level viability strongly depends upon careful management of human resources and fixed assets.

A steady state redox zone approach for modeling the transport and degradation of xenobiotic organic compounds from a landfill site

A redox zonation approach is used as a framework for obtaining biodegradation rate constants of xenobiotic compounds in a landfill plume (Grindsted, Denmark). The aquifer is physically heterogeneous in terms of a complex zonation of different geological units close to the landfill and biogeochemically heterogeneous in terms of a specified redox zonation. First-order degradation rates of six organic compounds (benzene, toluene, ethylbenzene, o-xylene, m/p-xylene, and naphthalene) were calculated in the methanogenic/sulfate- and Fe-reducing zones. The numerical simulations show that all compounds are anaerobically biodegraded, but at very different rates. High rates of biodegradation of most of the compounds (except benzene) were found in the Fe-reducing zone. These rates generally agree with previously published rates. Only o-xylene and toluene were significantly biodegraded in the methanogenic/sulfate-reducing environment. All rates in this redox zone are generally much lower than previously published rates.

Fringe-controlled natural attenuation of phenoxy acids in a landfill plume: integration of field-scale processes by reactive transport modeling

Data obtained from a field study of an aquifer contaminated by landfill leachate and related laboratory experiments suggest that natural attenuation of phenoxy acid herbicides such as mecoprop (MCPP) occurs in the transition zone between the anaerobic plume core and the overlying aerobic water body. The location of this transition zone is assumed to be primarily controlled by vertical transverse dispersion processes occurring downstream of the pollution source. A reactive transport modeling study was carried out to evaluate this conceptual model. The transport was simulated for a two-dimensional vertical cross section to quantify the combined physical, geochemical, and microbial processes that affect the fate of the phenoxy acid herbicides. The simulations, showing removal of phenoxy acids, an increase of phenoxy acid degraders in the fringe zone, and a dependency of the results on vertical transverse dispersivity, are compatible with the hypothesis of fringe-controlled aerobic biodegradation of the phenoxy acids.

Geobacteraceae community composition is related to hydrochemistry and biodegradation in an iron-reducing aquifer polluted by a neighboring landfill

Relationships between community composition of the iron-reducing Geobacteraceae, pollution levels, and the occurrence of biodegradation were established for an iron-reducing aquifer polluted with landfill leachate by using cultivation-independent Geobacteraceae 16S rRNA gene-targeting techniques. Numerical analysis of denaturing gradient gel electrophoresis (DGGE) profiles and sequencing revealed a high Geobacteraceae diversity and showed that community composition within the leachate plume differed considerably from that of the unpolluted aquifer. This suggests that pollution has selected for specific species out of a large pool of Geobacteraceae. DGGE profiles of polluted groundwater taken near the landfill (6- to 39-m distance) clustered together. DGGE profiles from less-polluted groundwater taken further downstream did not fall in the same cluster. Several individual DGGE bands were indicative of either the redox process or the level of pollution. This included a pollution-indicative band that dominated the DGGE profiles from groundwater samples taken close to the landfill (6 to 39 m distance). The clustering of these profiles and the dominance by a single DGGE band corresponded to the part of the aquifer where organic micropollutants and reactive dissolved organic matter were attenuated at relatively high rates.

Quantification of sequential chlorinated ethene degradation by use of a reactive transport model incorporating isotope fractionation

Compound-specific isotope analysis (CSIA) enables quantification of biodegradation by use of the Rayleigh equation. The Rayleigh equation fails, however, to describe the sequential degradation of chlorinated aliphatic hydrocarbons (CAHs) involving various intermediates that are controlled by simultaneous degradation and production. This paper shows how isotope fractionation during sequential degradation can be simulated in a 1D reactive transport code (PHREEQC-2). 12C and 13C isotopes of each CAH were simulated as separate species, and the ratio of the rate constants of the heavy to light isotope equaled the kinetic isotope fractionation factor for each degradation step. The developed multistep isotope fractionation reactive transport model (IF-RTM) adequately simulated reductive dechlorination of tetrachloroethene (PCE) to ethene in a microcosm experiment. Transport scenarios were performed to evaluate the effect of sorption and of different degradation rate constant ratios among CAH species on the downgradient isotope evolution. The power of the model to quantify degradation is illustrated for situations where mixed sources degrade and for situations where daughter products are removed by oxidative processes. Finally, the model was used to interpret the occurrence of reductive dechlorination at a field site. The developed methodology can easily be incorporated in 3D solute transport models to enable quantification of sequential CAH degradation in the field by CSIA.

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.