Boris M van Breukelen

Abstract: Groundwater systems are increasingly used for seasonal aquifer thermal energy storage (SATES) for periodic heating and cooling of buildings. Its use is hampered in contaminated aquifers because of the potential environmental risks associated with the spreading of contaminated groundwater, but positive side effects, such as enhanced contaminant remediation, might also occur. A first reactive transport study is presented to assess the effect of SATES on the fate of chlorinated solvents by means of scenario modeling, with emphasis on the effects of transient SATES pumping and applicable kinetic degradation regime. Temperature effects on physical, chemical, and biological reactions were excluded as calculations and initial simulations showed that the small temperature range commonly involved (Delta T<15 degrees C) only caused minor effects. The results show that a significant decrease of the contaminant mass and (eventually) plume volume occurs when degradation is described as sediment-limited with a constant rate in space and time, provided that dense non-aqueous phase liquid (DNAPL) is absent. However, in the presence of DNAPL dissolution, particularly when the dissolved contaminant reaches SATES wells, a considerably larger contaminant plume is created, depending on the balance between DNAPL dissolution and mass removal by degradation. Under conditions where degradation is contaminant-limited and degradation rates depend on contaminant concentrations in the aquifer, a SATES system does not result in enhanced remediation of a contaminant plume. Although field data are lacking and existing regulatory constraints do not yet permit the application of SATES in contaminated aquifers, reactive transport modeling provides a means of assessing the risks of SATES application in contaminated aquifers. The results from this study are considered to be a first step in identifying the subsurface conditions under which SATES can be applied in a safe or even beneficial manner. (C) 2013 Elsevier B.V. All rights reserved.

Abstract: This reactive transport modeling study presents a follow up to the mass balance-based identification and quantification of the main hydrogeochemical processes that occurred during an aquifer storage and recovery (ASR) trial in an anoxic sandy aquifer (Herten, the Netherlands). Kinetic rate expressions were used to simulate oxidation of pyrite, soil organic matter (SOM), and ferrous iron, and dissolution of calcite and Mn-siderite. Cation exchange, precipitation of Fe- and Mn-(hydr) oxides, and surface complexation were treated as equilibrium processes. The PHREEQC model was automatically calibrated with PEST to observations from the first ASR cycle, and was then allowed to run for all 14 cycles to evaluate its long term performance. A sensitivity analysis was conducted to find the most controlling model parameters. Pyrite was ranked as the most important reductant, followed by SOM, whereas Fe(II) was least important. Moreover, the pH and oxygen gradients were found to enhance the rate of pyrite over SOM oxidation with distance away from the ASR well. The increasing sorption capacity of precipitating Fe-hydroxides was reflected by the decreasing Fe(II) concentrations with subsequent cycles whereas Mn(II) showed a tendency to mobilize during recovery and remain above standards. Oxidation and dissolution rates were found to depend on travel time and injection rate as well as on the presence or absence of flow. Oxygen enrichment of the injection water increased oxidation rates and therefore accelerated the aquifer's leaching from its reactive species. We specifically focused on impeding the release of Mn(II) to the groundwater, a process that acted as a restraining factor for the feasibility of ASR application at this site. The undesirable side-effects of oxygen enrichment as well as the Mn(II) issues were found to be partly suppressed by enriching the source water with pH buffers according to scenario simulations. (C) 2013 Elsevier Ltd. All rights reserved.

Abstract: Aquifer thermal energy storage (ATES) uses groundwater to store energy for heating or cooling purposes in the built environment. This paper presents field and laboratory results aiming to elucidate the effects that ATES operation may have on chemical groundwater quality. Field data from an ATES site in the south of the Netherlands show that ATES results in chemical quality perturbations due to homogenisation of the initially present vertical water quality gradient. We tested this hypothesis by numerical modelling of groundwater flow and coupled SO4 transport during extraction and injection of groundwater by the ATES system. The modelling results confirm that extracting groundwater from an aquifer with a natural quality stratification, mixing this water in the ATES system, and subsequent injection in the second ATES well can adequately describe the observation data. This mixing effect masks any potential temperature effects in typical low temperature ATES systems (<25 degrees C) which was the reason to complement the field investigations with laboratory experiments focusing on temperature effects. The laboratory experiments indicated that temperature effects until 25 degrees C are limited; most interestingly was an increase in arsenic concentration. At 60 degrees C, carbonate precipitation, mobilisation of dissolved oxygen concentration, K and Li, and desorption of trace metals like As can occur.

Abstract: Carbon (C), chlorine (Cl), and hydrogen (H) isotope effects were determined during dechlorination of TCE to ethene by a mixed Dehalococcoides (Dhc) culture. The C isotope effects for the dechlorination steps were consistent with data published in the past for reductive dechlorination (RD) by Dhc. The Cl effects (combined with an inverse H effect in TCE) suggested that dechlorination proceeded through nucleophilic reactions with cobalamin rather than by an electron transfer mechanism. Depletions of 37Cl in daughter compounds, resulting from fractionation at positions away from the dechlorination center (secondary isotope effects), further support the nucleophilic dechlorination mechanism. Determination of C and Cl isotope ratios of the reactants and products in the reductive dechlorination chain offers a potential tool for differentiation of Dhc activity from alternative transformation mechanisms (e.g., aerobic degradation and reductive dechlorination proceeding via outer sphere mechanisms), in studies of in situ attenuation of chlorinated ethenes. Hydrogenation of the reaction products (DCE, VC, and ethene) showed a major preference for the 1H isotope. Detection of depleted dechlorination products could provide a line of evidence in discrimination between alternative sources of TCE (e.g., evolution from DNAPL sources or from conversion of PCE).

Abstract: Abstract Aquifers used for the production of drinking water are increasingly being used for the generation of shallow geothermal energy. This causes temperature perturbations far beyond the natural variations in aquifers and the effects of these temperature variations on groundwater quality, in particular trace elements, have not been investigated. Here, we report the results of column experiments to assess the impacts of temperature variations (5°C, 11°C, 25°C and 60°C) on groundwater quality in anoxic reactive unconsolidated sandy sediments derived from an aquifer system widely used for drinking water production in the Netherlands. Our results showed that at 5 °C no effects on water quality were observed compared to the reference of 11°C (in situ temperature). At 25°C, As concentrations were significantly increased and at 60 °C, significant increases were observed pH and DOC, P, K, Si, As, Mo, V, B, and F concentrations. These elements should therefore be considered for water quality monitoring programs of shallow geothermal energy projects. No consistent temperature effects were observed on Na, Ca, Mg, Sr, Fe, Mn, Al, Ba, Co, Cu, Ni, Pb, Zn, Eu, Ho, Sb, Sc, Yb, Ga, La, and Th concentrations, all of which were present in the sediment. The temperature-induced chemical effects were probably caused by (incongruent) dissolution of silicate minerals (K and Si), desorption from, and potentially reductive dissolution of, iron oxides (As, B, Mo, V, and possibly P and DOC), and mineralisation of sedimentary organic matter (DOC and P).

Abstract: The hydrogeochemical processes that took place during an aquifer storage and recovery (ASR) trial in a confined anoxic sandy aquifer (Herten, the Netherlands) were identified and quantified, using observation wells at 0.1, 8 and 25 m distance from the ASR well. Oxic drinking water was injected in 14 ASR cycles in the period 2000-2009. The main reactions consisted of the oxidation of pyrite, sedimentary organic matter, and (adsorbed) Fe(II) and Mn(II) in all aquifer layers (A-D), whereas the dissolution of carbonates (Mg-calcite and Mn-siderite) occurred mainly in aquifer layer D. Extinction of the mobilization of SO4, Fe(II), Mn(II), As, Co, Ni, Ca and total inorganic C pointed at pyrite and calcite leaching in layer A, whereas reactions with Mn-siderite in layer D did not show a significant extinction over time. Iron(II) and Mn(II) removal during recovery was demonstrated by particle tracking and pointed at sorption to neoformed ferrihydrite. Part of the oxidants was removed by neoformed organic material in the ASR proximal zone (0 -ca. 5 m) where micro-organisms grow during injection and die away when storage exceeds about 1 month. Anoxic conditions during storage led to increased concentrations for a.o. Fe(II), Mn(II) and NH4 as noted for the first 50-200 m(3) of abstracted water during the recovery phase. With a mass balance approach the water-sediment reactions and leaching rate of the reactive solid phases were quantified. Leaching of pyrite and calcite reached completion at up to 8 m distance in layer A, but not in layer D. The mass balance approach moreover showed that Mn-siderite in layer D was probably responsible for the Mn(II) exceedances of the drinking water standard (0.9 mu mol/L) in the recovered water. Leaching of the Mn-siderite up to 8 m from the ASR well would take 1600 more pore volumes of drinking water injection (on top of the realized 460). (C) 2012 Elsevier Ltd. All rights reserved.

Abstract: This article describes the post audit and inverse modeling of a 1-D forward reactive transport model. The model simulates the changes in water quality following artificial recharge of pre-treated water from the river Rhine in the Amsterdam Water Supply Dunes using the PHREEQC-2 numerical code. One observation dataset is used for model calibration, and another dataset for validation of model predictions. The total simulation time of the model is 50 years, from 1957 10 2007, with recharge composition varying on a monthly basis and the post audit is performed 26 years after the former model simulation period. The post audit revealed that the original model could reasonably predict conservative transport and kinetic redox reactions (oxygen and nitrate reduction coupled to the oxidation of soil organic carbon), but showed discrepancies in the simulation of cation exchange. Conceptualizations of the former model were inadequate to accurately simulate water quality changes controlled by cation exchange, especially concerning the breakthrough of potassium and magnesium fronts. Changes in conceptualization and model design, including the addition of five flow paths, to a total of six, and the use of parameter estimation software (PEST), resulted in a better model to measurement fit and system representation. No unique parameter set could be found for the model, primarily due to high parameter correlations, and an assessment of the predictive error was made using a calibration constrained Monte-Carlo method, and evaluated against field observations. The predictive error was found to be low for Na+ and Ca2+, except for greater travel times, while the K+ and Mg2+ error was restricted to the exchange fronts at some of the flow paths. Optimized cation exchange coefficients were relatively high, especially for potassium, but still within the observed range in literature. The exchange coefficient for potassium agrees with strong fixation on illite, a main clay mineral in the area. Optimized CEC values were systematically lower than clay and organic matter contents indicated, possibly reflecting preferential flow of groundwater through the more permeable but less reactive aquifer parts. Whereas the artificial recharge initially acted as an intrusion of relatively saline water triggering Na+ for Ca2+ exchange, further increasing total hardness of the recharged water, the gradual long-term reduction in salinity of the river Rhine since the mid 1970s has shifted to an intrusion of fresher water causing Ca2+ for Na+ exchange. As a result, seasonal and longer term reversal of the initial cation exchange processes was observed adding to the general long-term reduction in total hardness of the recharged Rhine water. (C) 2012 Elsevier B.V. All rights reserved.

Abstract: 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 parts per thousand) 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, alpha, of the degradation reaction. Models disregarding DIF underestimate the actual alpha. The approximately twice larger DIF effect for chlorine than for carbon together with the low alpha for oxidation resulted in strong chlorine CSIA depletions for VC at the plume fringe.

Abstract: There is increasing interest in combined carbon-chlorine compound-specific isotope analysis (CSIA) to differentiate between contaminant sources and to assess transformation processes, However, the significant abundance of polychlorinated molecules with several heavy chlorine isotopes complicates the evaluation of chlorine isotope trends. Therefore, the goal of this study was to develop a conceptual and mathematical framework that describes the expected chlorine isotope fractionation patterns during multistep transformation of chlorinated compounds. Reductive dechlorination of chlorinated ethenes served as an example. The study demonstrates that chlorine isotope trends can be simulated by reproducing the average behavior of light and heavy isotopes or by explicitly simulating molecules with different numbers of heavy isotopes (isotopologues). The calculations reveal that initial chlorine isotope ratios of products equal the isotope ratios of their parent compounds in the absence of secondary isotope effects, while steadily increasing during transformation. The slopes in dual isotope plots are linear for reactant and product during a one-step reaction. They become nonlinear for products that are degraded further but converge to characteristic slopes. Consideration of different scenarios reveals that combined carbon-chlorine isotope analysis bears high potential to differentiate between contaminant sources, to elucidate reaction mechanisms in laboratory studies, and to identify transformation processes in the field.

Abstract: Open waste dump systems are still widely used in Indonesia. The Jatibarang landfill receives 650-700 tons of municipal waste per day from the city of Semarang, Central Java. Some of the leachate from the landfill (lows via several natural and collection ponds to a nearby river. The objectives of the study were to identify seasonal landfill leachate characteristics in this surface water and to determine the occurrence of natural attenuation, in particular the potential for biodegradation, along the flow path. Monthly measurements of general landfill leachate parameters. organic matter-related factors and redox-related components revealed that leachate composition was influenced by seasonal precipitation. In the dry season, electrical conductivity and concentrations of BOD, COD, N-organic matter, ammonia, sulphate and calcium were significantly higher (1.1-2.3 fold) than during the wet season. Dilution was the major natural attenuation process acting on leachate. Heavy metals had the highest impact on river water quality. Between the landfill and the river, a fivefold dilution occurred during the dry season due to active springwater infiltration, while rainwater led to a twofold dilution in the wet season. Residence time of leachate in the surface leachate collection system was less than 70 days. Field measurements and laboratory experiments showed that during this period hardly any biodegradation of organic matter and ammonia occurred (less than 25%). However, the potential for biodegradation of organic matter and ammonia was clearly revealed during 700 days of incubation of leachate in the laboratory (over 65%). If the residence time of leachate discharge can be increased to allow for biodegradation processes and precipitation reactions, the polluting effects of leachate on the river can be diminished. (C) 2008 Elsevier Ltd. All rights reserved.

Abstract: 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 for the 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.

Abstract: Heterogeneity in eukaryotic and bacteria community structure in surface and subsurface sediment samples downgradient of the Banisveld landfill (The Netherlands) was studied using a culturing-independent molecular approach. Along a transect covering the part of the aquifer most polluted by landfill leachate, sediment was sampled at 1-m depth intervals, until a depth of 5.5 m, at four distances from the landfill. Two drillings were placed in a nearby clean area as a reference. Denaturing gradient gel electrophoresis banding patterns revealed high bacterial and eukaryotic diversity and complex community structures. Bacteria and eukaryotic community profiles in polluted samples grouped different from those in clean samples. Bacteria community profiles in surface samples clustered together and separately from subsurface community profiles. Subsurface bacteria profiles clustered in a location-specific manner. Eukaryotic community structure did not significantly relate to distance from the landfill or depth. No significant spatial autocorrelation of bacteria or eukaryotic communities was observed over 1-m depth intervals per sampling location. Spatial heterogeneity in sediment-associated bacterial communities appears to be much larger than in groundwater. We discuss how on the one hand, spatial heterogeneity may complicate the assessment of microbial community structure and functioning, while on the other it may provide better opportunities for natural attenuation.

Abstract: 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 I 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.

Abstract: Identification of the functional groups of microorganisms that are predominantly in control of fluxes through, and concentrations in, microbial networks would benefit microbial ecology and environmental biotechnology: the properties of those controlling microorganisms could be studied or monitored specifically or their activity could be modulated in attempts to manipulate the behaviour of such networks. Herein we present ecological control analysis (ECA) as a versatile mathematical framework that allows for the quantification of the control of each functional group in a microbial network on its process rates and concentrations of intermediates. In contrast to current views, we show that rates of flow of matter are not always limited by a single functional group; rather flux control can be distributed over several groups. Also, control over intermediate concentrations is always shared. Because of indirect interactions, through other functional groups, the concentration of an intermediate can also be controlled by functional groups not producing or consuming it. Ecological control analysis is illustrated by a case study on the anaerobic degradation of organic matter, using experimental data obtained from the literature. During anaerobic degradation, fermenting microorganisms interact with terminal electron-accepting microorganisms (e.g. halorespirers, methanogens). The analysis indicates that flux control mainly resides with fermenting microorganisms, but can shift to the terminal electron-accepting microorganisms under less favourable redox conditions. Paradoxically, halorespiring microorganisms do not control the rate of perchloroethylene and trichloroethylene degradation even though they catalyse those processes themselves.

Abstract: The Rayleigh equation relates the change in isotope ratio of an element in a substrate to the extent of substrate consumption via a single kinetic isotopic fractionation factor (alpha). Substrate consumption is, however, commonly distributed over several metabolic pathways each potentially having a different alpha. Therefore, extended Rayleigh-type equations were derived to account for multiple competing degradation pathways. The value of alpha as expressed in the environment appears a function of the alpha values and rate constants of the various involved degradation pathways. Remarkably, the environmental or apparent alpha value changes and shows non-Rayleigh behavior over a large and relevant concentration interval if Monod kinetics applies and the half-saturation constants of the competing pathways differ. Derived equations were applied to previously published data and enabled (i) quantification of the share that two competing degradation pathways had on aerobic 1,2-dichloroethane (1,2-DCA) biodegradation in laboratory batch experiments and (ii) calculation of the extent of methyl tert-butyl ether (MTBE) biodegradation shared over aerobic and anaerobic degradation at a field site by means of an improved solution to two-dimensional (carbon and hydrogen) compound-specific isotope analysis (CSIA).

Abstract: Using culture-independent 16S rRNA gene-based methods, we previously observed that Geobacteraceae were a major component of the microbial communities in the iron-reducing aquifer polluted by the Banisveld landfill, The Netherlands. However, phylogenetic information does not tell about the functional potential of the detected Geobacteraceae, nor can phylogenetic information easily be used to establish the presence of other iron-reducers. Therefore, we enriched for iron-reducing consortia using a range of culturing media, with various electron donors and acceptors and varying incubation conditions ( pH, temperature), and by applying dilution-to-extinction culturing. Enrichments and strains isolated from these enrichments were characterized by 16S rRNA gene-based methods. The number of culturable iron-reducers was less than 110 iron-reducing bacteria per gram of sediment. The Geobacter phylotype that was previously found to constitute a major part of the microbial communities in a part of the aquifer where organic matter was attenuated at a relatively high rate, was not isolated. The isolation of another Geobacter strain and Serratia, Clostridium, Rhodoferax and Desulfitobacterium strains suggest the presence of a diverse iron-reducing community. Physiological capabilities of the isolates are described and discussed in relation to the hydrogeochemistry and the high abundance of Geobacteraceae in the aquifer polluted by the Banisveld landfill.

Abstract: Quantifying the share of destructive and nondestructive processes to natural attenuation (NA) of groundwater pollution plumes is of high importance to the evaluation and acceptance of NA as remediation strategy. Dilution as consequence of hydrodynamic dispersion may contribute considerably to NA, however, without reducing the mass of pollution. Unfortunately, tracers to quantify dilution are usually lacking. Degradation though of low-molecular-weight organic chemicals such as BTEX, chlorinated ethenes, and MTBE is uniquely associated with increases in isotope ratios for steady-state plumes. Compound-specific isotope analysis (CSIA) data are commonly interpreted by means of the Rayleigh equation, originally developed for closed systems, to calculate the extent of degradation under open system field conditions. For that reason, the validity of this approach has been questioned. The Rayleigh equation was accordingly modified to account for dilution, and showed that dilution contributed several to many times more to NA than biodegradation at a groundwater benzene plume. Derived equations also (i) underlined that field-derived isotopic enrichment factors underestimate actual values operative as a consequence of dilution, and (ii) provided a check on the lower limit of isotopic fractionation, thereby resulting in more reliable predictions on the extent of degradation.

Abstract: 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 (PHREEGC-2). C-12 and C-13 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 tetra chloroethene (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.

Abstract: Managers of landfill sites are faced with enormous challenges when attempting to detect and delineate leachate plumes with a limited number of monitoring wells, assess spatial and temporal trends for hundreds of contaminants, and design long-term monitoring (LTM) strategies. Subsurface microbial ecology is a unique source of data that has been historically underutilized in LTM groundwater designs. This paper provides a methodology for utilizing qualitative and quantitative information (specifically, multiple water quality measurements and genome-based data) from a landfill leachate contaminated aquifer in Banisveld, The Netherlands, to improve the estimation of parameters of concern. We used a principal component analysis (PCA) to reduce nonindependent hydrochemistry data, Bacteria and Archaea community profiles from 16S rDNA denaturing gradient gel electrophoresis (DGGE), into six statistically independent variables, representing the majority of the original dataset variances. The PCA scores grouped samples based on the degree or class of contamination and were similar over considerable horizontal distances. Incorporation of the principal component scores with traditional subsurface information using cokriging improved the understanding of the contaminated area by reducing error variances and increasing detection efficiency. Combining these multiple types of data (e.g., genome-based information, hydrochemistry, borings) may be extremely useful at landfill or other LTM sites for designing cost-effective strategies to detect and monitor contaminants.

Abstract: 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.

Abstract: The biogeochemical processes goveming leachate attenuation inside a landfill leachate plume (Banisveld, the Netherlands) were revealed and quantified using the ID reactive transport model PHREEQC-2. Biodegradation of dissolved organic carbon (DOC) was simulated assuming first-order oxidation of two DOC fractions with different reactivity, and was coupled to reductive dissolution of iron oxide. The following secondary geochemical processes were required in the model to match observations: kinetic precipitation of calcite and siderite, cation exchange, proton buffering and degassing. Rate constants for DOC oxidation and carbonate mineral precipitation were determined, and other model parameters were optimized using the nonlinear optimization program PEST by means of matching hydrochemical observations closely (pH, DIC, DOC, Na, K, Ca, Mg, NH4, Fe(II), SO4, Cl, CH4, saturation index of calcite and siderite). The modelling demonstrated the relevance and impact of various secondary geochemical processes on leachate plume evolution. Concomitant precipitation of siderite masked the act of iron reduction. Cation exchange resulted in release of Fe(II) from the pristine anaerobic aquifer to the leachate. Degassing, triggered by elevated CO2 pressures caused by carbonate precipitation and proton buffering at the front of the plume, explained the observed downstream decrease in methane concentration.

Abstract: 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 O-2, 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 (delta(13)C) 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. (C) 2004 Elsevier B.V. All rights reserved.

Abstract: The biogeochemical processes were identified which improved the leachate composition in the flow direction of a landfill leachate plume (Banisveld, The Netherlands). Groundwater observation wells were placed at specific locations after delineating the leachate plume using geophysical tests to map subsurface conductivity. Redox processes were determined using the distribution of solid and soluble redox species, hydrogen concentrations, concentration of dissolved gases (N-2, Ar, and CH4), and stable isotopes (delta(15)N-NO3, delta(34)S-SO4, delta(13)C-CH4, delta(2)H-CH4, and delta(13)C of dissolved organic and inorganic carbon (DOC and DIC, respectively)). The combined application of these techniques improved the redox interpretation considerably. Dissolved organic carbon (DOC) decreased downstream in association with increasing delta(13)C-DOC values confirming the occurrence of degradation. Degradation of DOC was coupled to iron reduction inside the plume, while denitrification could be an important redox process at the top fringe of the plume. stable carbon and hydrogen isotope signatures of methane indicated that methane was formed inside the landfill and not in the plume. Total gas pressure exceeded hydrostatic pressure in the plume, and methane seems subject to degassing. Quantitative proof for DOC degradation under iron-reducing conditions could only be obtained if the geochemical processes cation exchange and precipitation of carbonate minerals (siderite and calcite) were considered and incorporated in an inverse geochemical model of the plume. Simulation of delta(13)C-DIC confirmed that precipitation of carbonate minerals happened. (C) 2003 Elsevier Science B.V. All rights reserved.

Abstract: Compound-specific carbon and hydrogen isotope analysis was used to investigate biodegradation of benzene and ethylbenzene in contaminated groundwater at Dow Benelux BV industrial Site. delta(13)C values for dissolved benzene and ethylbenzene in downgradient samples were enriched by up to 2 +/- 0.5parts per thousand in C-13, compared to the delta(13)C value of the source area samples. delta(2)H values for dissolved benzene and ethylbenzene in downgradient samples exhibited larger isotopic enrichments of up to 27 +/- 5parts per thousand for benzene and up to 50 +/- 5parts per thousand for ethylbenzene relative to the source area. The observed carbon and hydrogen isotopic fractionation in downgradient samples provides evidence of biodegradation of both benzene and ethylbenzene within the study area at Dow Benelux BV. The estimated extents of biodegradation of benzene derived from carbon and hydrogen isotopic compositions for each sample are in agreement, supporting the conclusion that biodegradation is the primary control on the observed differences in carbon and hydrogen isotope values. Combined carbon and hydrogen isotope analyses provides the ability to compare biodegradation in the field based on two different parameters, and hence provides a stronger basis for assessment of biodegradation of petroleum hydrocarbon contaminants.

Abstract: Knowledge about the relationship between microbial community structure and hydrogeochemistry (e.g., pollution, redox and degradation processes) in landfill leachate-polluted aquifers is required to develop tools for predicting and monitoring natural attenuation. In this study analyses of pollutant and redox chemistry were conducted in parallel with culture-independent profiling of microbial communities present in a well-defined aquifer (Banisveld, The Netherlands). Degradation of organic contaminants occurred under iron-reducing conditions in the plume of pollution, while upstream of the landfill and above the plume denitrification was the dominant redox process. Beneath the plume iron reduction occurred. Numerical comparison of 16S ribosomal DNA (rDNA)based denaturing gradient gel electrophoresis (DGGE) profiles of Bacteria and Archaea in 29 groundwater samples revealed a clear difference between the microbial community structures inside and outside the contaminant plume. A similar relationship was not evident in sediment samples. DGGE data were supported by sequencing cloned 16S rDNA. Upstream of the landfill members of the beta subclass of the class Proteobacteria (beta -proteobacteria) dominated. This group was not encountered beneath the landfill, where gram-positive bacteria dominated. Further downstream the contribution of gram-positive bacteria to the clone library decreased, while the contribution of delta -proteobacteria strongly increased and beta -proteobacteria reappeared. The beta -proteobacteria (Acidovorax, Rhodoferax) differed considerably from those found upstream (Gallionella, Azoarcus). Direct comparisons of cloned 16S rDNA with bands in DGGE profiles revealed that the data from each analysis were comparable. A relationship was observed between the dominant redox processes and the bacteria identified. In the iron-reducing plume members of the family Geobacteraceae made a strong contribution to the microbial communities. Because the only known aromatic hydrocarbon-degrading, iron-reducing bacteria are Geobacter spp., their occurrence in landfill leachate-contaminated aquifers deserves more detailed consideration.

Abstract: Databases containing information regarding presence and activity of microbial communities will be very useful for determination of the potential for intrinsic bioremediation in landfill leachate polluted aquifers. Simple analyses such as community-level physiological profiling (CLPP) and denaturing gradient gel electrophoresis (DGGE) of 16S rDNA fragments yield large sets of data for inclusion into such databases. In this study we describe the development of a method for anaerobic CLPP, using commercially available Biolog plates. Incubation at the in situ temperature of the aquifer (10 degreesC) for 28 days was optimal for obtaining a specific, reproducible physiological profile. Anaerobic incubation was essential for profiling anaerobic communities. The anaerobic cultivation-dependent CLPP method and cultivation-independent DGGE were applied to groundwater and sediment samples from the aquifer near the Coupepolder landfill in The Netherlands. A combination of computer-assisted CLPP and DGCE analysis of both groundwater and sediment samples yielded the best separating power for characterizing microbial communities in the aquifer. Communities in groundwater were significantly different from those in the corresponding sediment. Microbial communities present in subsamples from sediment cores usually were similar for the various sampling locations. Variation was observed for the heterogeneous sediment beneath the landfill. Both anaerobic CLPP and DGGE analysis clearly separated microbial communities from the polluted aquifer underneath the landfill from those in the less or not polluted aquifer downstream and upstream of the landfill.

Abstract: Previously, we observed that microbial community structure and functional diversity in aquifers might be enhanced by landfill leachate infiltration. To study this hypothesis, groundwater samples were taken near the Banisveld landfill, The Netherlands. Based on hydrochemical parameters, the samples clustered into two groups. One group corresponded to polluted samples from the plume of landfill leachate and the second group to clean samples from outside the plume. Most Probable Number-Biolog was used to select Eco Biolog plates with similar inoculum densities. Analysis of substrate utilization profiles of these plates revealed that anaerobic microbial communities in polluted samples clustered separately from those in clean samples. Especially substrates containing an aromatic nucleus were more utilized by microbial communities in the leachate plume. Both substrate richness and functional diversity were significantly enhanced in the plume of pollution. This study shows that community-level physiological profiling is a useful and simple tool to distinguish between anaerobic microbial communities in and near a landfill leachate plume.

Abstract: Water quality changes were modelled along a flowpath in a plume of artificially recharged, pretreated Rhine water in the dunes of the Amsterdam Water Supply, after 24 years of infiltration. The hydrogeochemical transport model PHREEQC was extended with dispersion/diffusion and kinetics for selected chemical reactions. In the model the following reactions were included: cation-exchange, calcite dissolution and precipitation, and kinetic oxygen consumption and denitrification by oxidation of organic matter. Monthly-averaged values were used for the infiltration water quality. Traveltimes from infiltration area to sampling points were determined with chloride and tritium, and used to place the 3D field-observations in the 1D column-model. Values for CEC were variable for seven layers in the model. Infiltration of pretreated Rhine water in the dune aquifer can be considered an intrusion of more saline water. It caused desorption of Ca2+, in exchange for Na+, K+ and Mg2+ from Rhine water. Because of variations in total solute concentrations in infiltration water, local small scale freshening fronts (Ca2+ sorption, Na+ desorption) were created by seasonally decreasing salt concentrations. The undersaturation with respect to calcite in the infiltration water, and the CO2 produced during consumption of oxygen, resulted in dissolution of calcite. Precipitation of calcite occurred in response to desorption of calcium from the exchanger in the downstream parts, overall, a net dissolution of calcite was simulated. Good results were generally achieved for all components: sulfate, nitrate, chloride, alkalinity, calcium, magnesium, potassium, sodium, H-3 and O-2. The contributions of the different geochemical reactions to the water quality are illustrated with computer simulations for the individual processes. (C) 1998 Elsevier Science B.V. All rights reserved.