Modeling kinetic processes controlling hydrogen and acetate concentrations in an aquifer-derived microcosm

Microbiological hydrogen (H2 ) thresholds in anaerobic continuous-flow systems: Effects of system characteristics

Hydrogen (H₂ ) concentrations that were associated with microbiological respiratory processes (RPs) such as sulfate reduction and methanogenesis were quantified in continuous-flow systems (CFSs) (e.g., bioreactors, sediments). Gibbs free energy yield (ΔǴ ~ 0) of the relevant RP has been proposed to control the observed H₂ concentrations, but most of the reported values do not align with the proposed energetic trends. Alternatively, we postulate that system characteristics of each experimental design influence all system components including H₂ concentrations. To analyze this proposal, a Monod-based mathematical model was developed and used to design a gas-liquid bioreactor for hydrogenotrophic methanogenesis with Methanobacterium bryantii M.o.H. Gas-to-liquid H₂ mass transfer, microbiological H₂ consumption, biomass growth, methane formation, and Gibbs free energy yields were evaluated systematically. Combining model predictions and experimental results revealed that an initially large biomass concentration created transients during which biomass consumed [H₂ ]_L rapidly to the thermodynamic H₂ -threshold (≤1 nM) that triggerred the microorganisms to stop H₂ oxidation. With no H₂ oxidation, continuous gas-to-liquid H₂ transfer increased [H₂ ]_L to a level that signaled the methanogens to resume H₂ oxidation. Thus, an oscillatory H₂ -concentration profile developed between the thermodynamic H₂ -threshold (≤1 nM) and a low [H₂ ]_L (~10 nM) that relied on the rate of gas-to-liquid H₂ -transfer. The transient [H₂ ]_L values were too low to support biomass synthesis that could balance biomass losses through endogenous oxidation and advection; thus, biomass declined continuously and disappeared. A stable [H₂ ]_L (1807 nM) emerged as a result of abiotic H₂ -balance between gas-to-liquid H₂ transfer and H₂ removal via advection of liquid-phase.

Vertical Hydrologic Exchange Flows Control Methane Emissions from Riverbed Sediments

CH₄ emissions from inland waters are highly uncertain in the current global CH₄ budget, especially for streams, rivers, and other lotic systems. Previous studies have attributed the strong spatiotemporal heterogeneity of riverine CH₄ to environmental factors such as sediment type, water level, temperature, or particulate organic carbon abundance through correlation analysis. However, a mechanistic understanding of the basis for such heterogeneity is lacking. Here, we combine sediment CH₄ data from the Hanford reach of the Columbia River with a biogeochemical-transport model to show that vertical hydrologic exchange flows (VHEFs), driven by the difference between river stage and groundwater level, determine CH₄ flux at the sediment-water interface. CH₄ fluxes show a nonlinear relationship with the magnitude of VHEFs, where high VHEFs introduce O₂ into riverbed sediments, which inhibit CH₄ production and induce CH₄ oxidation, and low VHEFs cause transient reduction in CH₄ flux (relative to production) due to reduced advective CH₄ transport. In addition, VHEFs lead to the hysteresis of temperature rise and CH₄ emissions because high river discharge caused by snowmelt in spring leads to strong downwelling flow that offsets increasing CH₄ production with temperature rise. Our findings reveal how the interplay between in-stream hydrologic flux besides fluvial-wetland connectivity and microbial metabolic pathways that compete with methanogenic pathways can produce complex patterns in CH₄ production and emission in riverbed alluvial sediments.

Release of trace elements during bioreductive dissolution of magnetite from metal mine tailings: Potential impact on marine environments

Adverse impacts of mine tailings on water and sediments quality are major worldwide environmental problems. Due to the environmental issues associated with the deposition of mine tailings on land, a controversial discussed alternative is submarine tailings disposal (STD). However, Fe(III) bioreduction of iron oxides (e.g., magnetite) in the tailings disposed might cause toxic effects on coastal environments due to the release of different trace elements (TEs) contained in the oxides. To study the extent and kinetics of magnetite bioreduction under marine conditions and the potential release of TEs, a number of batch experiments with artificial seawater (pH 8.2) and a marine microbial strain (Shewanella loihica) were performed using several magnetite ore samples from different mines and a mine tailings sample. The elemental composition of the magnetite determined in the tailings showed relatively high amounts of TEs (e.g., Mn, Zn, Co) compared with those of the magnetite ore samples (LA-ICP-MS and EMPA analyses). The experiments were conducted at 10 °C in the dark for up to 113 days. Based on the consumption of lactate and production of acetate and aqueous Fe(II) over time, the magnitude of Fe(III) bioreduction was calculated using a geochemical model including Monod kinetics. Model simulations reproduced the release of iron and TEs observed throughout the experiments, e.g., Mn (up to 203 μg L^-1), V (up to 79 μg L^-1), As (up to 17 μg L^-1) and Cu (up to 328 μg L^-1), suggesting a potential contamination of pore water by STD. Therefore, the results of this study can help to better evaluate the potential impacts of STD.

Multiscale Modeling of Uranium Bioreduction in Porous Media by One-Dimensional Biofilms

We formulate a multiscale mathematical model that describes the bioreduction of uranium in porous media. On the mesoscale we describe the bioreduction of uranium [VI] to uranium [IV] using a multispecies one-dimensional biofilm model with suspended bacteria and thermodynamic growth inhibition. We upscale the mesoscopic (colony scale) model to the macroscale (reactor scale) and investigate the behavior of substrate utilization and production, attachment and detachment processes, and thermodynamic effects not usually considered in biofilm growth models. Simulation results of the reactor model indicate that thermodynamic inhibition quantitatively alters the dynamics of the model and neglecting thermodynamic effects may over- or underestimate chemical concentrations in the system. Furthermore, we numerically investigate uncertainties related to the specific choice of attachment and detachment rate coefficients and find that while increasing the attachment rate coefficient or decreasing the detachment rate coefficient leads to thicker biofilms, performance of the reactor remains largely unaffected.

Thermodynamic Inhibition in Chemostat Models : With an Application to Bioreduction of Uranium

We formulate a mathematical model of bacterial populations in a chemostat setting that also accounts for thermodynamic growth inhibition as a consequence of chemical reactions. Using only elementary mathematical and chemical arguments, we carry this out for two systems: a simple toy model with a single species, a single substrate, and a single reaction product, and a more involved model that describes bioreduction of uranium[VI] into uranium[IV]. We find that in contrast to most traditional chemostat models, as a consequence of thermodynamic inhibition the equilibria concentrations of nutrient substrates might depend on their inflow concentration and not only on reaction parameters and the reactor's dilution rate. Simulation results of the uranium degradation model indicate that thermodynamic growth inhibition quantitatively alters the solution of the model. This suggests that neglecting thermodynamic inhibition effects in systems where they play a role might lead to wrong model predictions and under- or over-estimate the efficacy of the process under investigation.

Advanced redox zonation of the San Pedro Sula alluvial aquifer (Honduras) using data fusion and multivariate geostatistics

The incorrect wastewater management and the land use distribution lead to severe environmental problems, creating heavy eutrophication condition in surface-water; when surface-water/groundwater relationships exist, the organic matter transferred to the aquifer oxidizes and triggers redox processes (i.e. Terminal Electron Accepting Processes, TEAPs), that provoke severe groundwater quality modifications and complicate its exploitation and management. For this reason, the definition of the redox zonation within an aquifer is an effective tool for the identification of the contamination sources and for the conceptual model refinement, when remediation strategies need to be planned. Although the redox processes are dynamic reactions, the aquifer redox zonation is generally aimed to identify homogenous zones, characterized by a predominant TEAP. To overcome this methodological approach, the Multi-Collocated Factorial Kriging (MCFK) was applied to redox-related physico-chemical parameters, which allowed identifying their spatial relationships at different scales, transferring this method from precision agriculture and soil science to hydrogeochemistry. The selected study area is the San Pedro Sula aquifer (Honduras), a multi-layer alluvial aquifer characterized by well-known surface-water/groundwater interactions and heavy eutrophicated streams. Here, high concentrations of Mn and Fe were found in the aquifer. The MCFK results identified a short-range (2300 m) factor, highlighting a strong relation between Mn concentrations and anoxic conditions, due to the organic matter transfer from eutrophicated surface-water into the aquifer. Simultaneously, the relationship between Fe and turbidity is related to a fine Fe(III) oxi-hydroxide colloidal phase, developed when different redox conditions of groundwater mix up in the wells. The long-range (6000 m) factor points out that Fe is related to redox processes at a wider scale, especially in the northern San Pedro Sula alluvial plain. These results are supported by both the Principal Component Analysis and the hydrogeochemical numerical modeling. As a result, different TEAPs occur simultaneously in contaminated areas, acting at multiple scales.

Linking Spectral Induced Polarization (SIP) and Subsurface Microbial Processes: Results from Sand Column Incubation Experiments

Geophysical techniques, such as spectral induced polarization (SIP), offer potentially powerful approaches for in situ monitoring of subsurface biogeochemistry. The successful implementation of these techniques as monitoring tools for reactive transport phenomena, however, requires the deconvolution of multiple contributions to measured signals. Here, we present SIP spectra and complementary biogeochemical data obtained in saturated columns packed with alternating layers of ferrihydrite-coated and pure quartz sand, and inoculated with Shewanella oneidensis supplemented with lactate and nitrate. A biomass-explicit diffusion-reaction model is fitted to the experimental biogeochemical data. Overall, the results highlight that (1) the temporal response of the measured imaginary conductivity peaks parallels the microbial growth and decay dynamics in the columns, and (2) SIP is sensitive to changes in microbial abundance and cell surface charging properties, even at relatively low cell densities (<10⁸ cells mL^-1). Relaxation times (τ) derived using the Cole-Cole model vary with the dominant electron accepting process, nitrate or ferric iron reduction. The observed range of τ values, 0.012-0.107 s, yields effective polarization diameters in the range 1-3 μm, that is, 2 orders of magnitude smaller than the smallest quartz grains in the columns, suggesting that polarization of the bacterial cells controls the observed chargeability and relaxation dynamics in the experiments.

Biodegradation of phenolic compounds and their metabolites in contaminated groundwater using microbial fuel cells

This is the first study demonstrating the biodegradation of phenolic compounds and their organic metabolites in contaminated groundwater using bioelectrochemical systems (BESs). The phenols were biodegraded anaerobically via 4-hydroxybenzoic acid and 4-hydroxy-3-methylbenzoic acid, which were retained by electromigration in the anode chamber. Oxygen, nitrate, iron(III), sulfate and the electrode were electron acceptors for biodegradation. Electro-active bacteria attached to the anode, producing electricity (~1.8mW/m(2)), while utilizing acetate as an electron donor. Electricity generation started concurrently with iron reduction; the anode was an electron acceptor as thermodynamically favorable as iron(III). Acetate removal was enhanced by 40% in the presence of the anode. However, enhanced removal of phenols occurred only for a short time. Field-scale application of BESs for in situ bioremediation requires an understanding of the regulation and kinetics of biodegradation pathways of the parent compounds to relevant metabolites, and the syntrophic interactions and carbon flow in the microbial community.

An improved cellular automaton method to model multispecies biofilms

Biomass-spreading rules used in previous cellular automaton methods to simulate multispecies biofilm introduced extensive mixing between different biomass species or resulted in spatially discontinuous biomass concentration and distribution; this caused results based on the cellular automaton methods to deviate from experimental results and those from the more computationally intensive continuous method. To overcome the problems, we propose new biomass-spreading rules in this work: Excess biomass spreads by pushing a line of grid cells that are on the shortest path from the source grid cell to the destination grid cell, and the fractions of different biomass species in the grid cells on the path change due to the spreading. To evaluate the new rules, three two-dimensional simulation examples are used to compare the biomass distribution computed using the continuous method and three cellular automaton methods, one based on the new rules and the other two based on rules presented in two previous studies. The relationship between the biomass species is syntrophic in one example and competitive in the other two examples. Simulation results generated using the cellular automaton method based on the new rules agree much better with the continuous method than do results using the other two cellular automaton methods. The new biomass-spreading rules are no more complex to implement than the existing rules.

Diversity of planktonic and attached bacterial communities in a phenol-contaminated sandstone aquifer

Polluted aquifers contain indigenous microbial communities with the potential for in situ bioremediation. However, the effect of hydrogeochemical gradients on in situ microbial communities (especially at the plume fringe, where natural attenuation is higher) is still not clear. In this study, we used culture-independent techniques to investigate the diversity of in situ planktonic and attached bacterial communities in a phenol-contaminated sandstone aquifer. Within the upper and lower plume fringes, denaturing gradient gel electrophoresis profiles indicated that planktonic community structure was influenced by the steep hydrogeochemical gradient of the plume rather than the spatial location in the aquifer. Under the same hydrogeochemical conditions (in the lower plume fringe, 30 m below ground level), 16S rRNA gene cloning and sequencing showed that planktonic and attached bacterial communities differed markedly and that the attached community was more diverse. The 16S rRNA gene phylogeny also suggested that a phylogenetically diverse bacterial community operated at this depth (30 mbgl), with biodegradation of phenolic compounds by nitrate-reducing Azoarcus and Acidovorax strains potentially being an important process. The presence of acetogenic and sulphate-reducing bacteria only in the planktonic clone library indicates that some natural attenuation processes may occur preferentially in one of the two growth phases (attached or planktonic). Therefore, this study has provided a better understanding of the microbial ecology of this phenol-contaminated aquifer, and it highlights the need for investigating both planktonic and attached microbial communities when assessing the potential for natural attenuation in contaminated aquifers.

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

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

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

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

Microbial physiology-based model of ethanol metabolism in subsurface sediments

A biogeochemical reaction model was developed based on microbial physiology to simulate ethanol metabolism and its influence on the chemistry of anoxic subsurface environments. The model accounts for potential microbial metabolisms that degrade ethanol, including those that oxidize ethanol directly or syntrophically by reducing different electron acceptors. Out of the potential metabolisms, those that are active in the environment can be inferred by fitting the model to experimental observations. This approach was applied to a batch sediment slurry experiment that examined ethanol metabolism in uranium-contaminated aquifer sediments from Area 2 at the U.S. Department of Energy Field Research Center in Oak Ridge, TN. According to the simulation results, complete ethanol oxidation by denitrification, incomplete ethanol oxidation by ferric iron reduction, ethanol fermentation to acetate and H(2), hydrogenotrophic sulfate reduction, and acetoclastic methanogenesis: all contributed significantly to the degradation of ethanol in the aquifer sediments. The assemblage of the active metabolisms provides a frame work to explore how ethanol amendment impacts the chemistry of the environment, including the occurrence and levels of uranium. The results can also be applied to explore how diverse microbial metabolisms impact the progress and efficacy of bioremediation strategies.

The application of a simplified method to map the aerobic acetate mineralization rates at the groundwater table of the Netherlands

A simplified method is used to assess the microbial activity of subsoils and soils across a broad geographic scale. Acetate was selected because it is a major intermediate in catabolic biochemical pathways. In order to get minimal disturbance, only a small amount of tritium labelled acetate and water is added to the subsoil material. After an incubation time, the subsoil material is separated from the water by centrifugation and the formed tritium labelled water is separated from the remaining acetate by evaporation. The data of 128 locations in the Netherlands were plotted in a soil map and were also compared with the depth, dry weight, electric conductivity, pH and nitrate concentration. The peat areas consisted of limed meadows with a high groundwater level whereas the sand areas often showed deeper groundwater levels and a lower pH. The subsoils at the groundwater table of the peat areas, which are in contact with soil air, showed a higher mineralization rate compared with the surface soils in our study. In contrast, the mineralization rate of the subsoil at the groundwater table of sandy soils showed on average a factor 30 lower rate. Nevertheless, the self purification capacity of the subsoil can be vital under weather conditions where the surface soil becomes less active.

Modelling of an enhanced PAH attenuation experiment and associated biogeochemical changes at a former gasworks site in southern Germany

Former manufactured gas plant sites often form a widespread contaminant source in the subsurface, leading to large plumes that contain a wide variety of tar-oil related compounds. Although most of these compounds eventually degrade naturally, the relevant processes tend to be slow and inefficient, often leaving active remediation as the only viable option to eliminate the risks of toxic substances to reach potential receptors such as surface waters or drinking water wells. In this study we use a reactive transport model to analyse the fate of a contaminant plume containing acenaphthene, methylbenzofurans and dimethylbenzofurans (i) prior to the installation of an active remediation scheme and (ii) for an enhanced remediation experiment during which O(2) and H(2)O(2) were added to the contaminated groundwater through a recirculation well. The numerical model developed for this study considers the primary contaminant degradation reactions (i.e., microbially mediated redox reactions) as well as secondary and competing mineral precipitation/dissolution reactions that affect the site's hydrochemistry and/or contaminant fate. The model was calibrated using a variety of constraints to test the uncertainty on model predictions resulting from the undocumented presence of reductants such as pyrite. The results highlight the important role of reactive transport modelling for the development of a comprehensive process understanding.

Soil engineering in vivo: harnessing natural biogeochemical systems for sustainable, multi-functional engineering solutions

Carbon sequestration, infrastructure rehabilitation, brownfields clean-up, hazardous waste disposal, water resources protection and global warming-these twenty-first century challenges can neither be solved by the high-energy consumptive practices that hallmark industry today, nor by minor tweaking or optimization of these processes. A more radical, holistic approach is required to develop the sustainable solutions society needs. Most of the above challenges occur within, are supported on, are enabled by or grown from soil. Soil, contrary to conventional civil engineering thought, is a living system host to multiple simultaneous processes. It is proposed herein that 'soil engineering in vivo', wherein the natural capacity of soil as a living ecosystem is used to provide multiple solutions simultaneously, may provide new, innovative, sustainable solutions to some of these great challenges of the twenty-first century. This requires a multi-disciplinary perspective that embraces the science of biology, chemistry and physics and applies this knowledge to provide multi-functional civil and environmental engineering designs for the soil environment. For example, can native soil bacterial species moderate the carbonate cycle in soils to simultaneously solidify liquefiable soil, immobilize reactive heavy metals and sequester carbon-effectively providing civil engineering functionality while clarifying the ground water and removing carbon from the atmosphere? Exploration of these ideas has begun in earnest in recent years. This paper explores the potential, challenges and opportunities of this new field, and highlights one biogeochemical function of soil that has shown promise and is developing rapidly as a new technology. The example is used to propose a generalized approach in which the potential of this new field can be fully realized.

Energetic constraints on H2-dependent terminal electron accepting processes in anoxic environments: a review of observations and model approaches

Microbially mediated terminal electron accepting processes (TEAPs) to a large extent control the fate of redox reactive elements and associated reactions in anoxic soils, sediments, and aquifers. This review focuses on thermodynamic controls and regulation of H2-dependent TEAPs, case studies illustrating this concept, and the quantitative description of thermodynamic controls in modeling. Other electron transfer processes are considered where appropriate. The work reviewed shows that thermodynamics and microbial kinetics are connected near thermodynamic equilibrium. Free energy thresholds for terminal respiration are physiologically based and often near -20 kJ mol(-1), depending on the mechanism of ATP generation; more positive free energy values have been reported under "starvation conditions" for methanogenesis and lower values for TEAPs that provide more energy. H2-dependent methanogenesis and sulfate reduction are under direct thermodynamic control in soils and sediments and generally approach theoretical minimum energy thresholds. If H2 concentrations are lowered by thermodynamically more potent TEAPs, these processes are inhibited. This principle is also valid for TEAPS providing more free energy, such as denitrification and arsenate reduction, but electron donor concentration cannot be lowered so that the processes reach theoretical energy thresholds. Thermodynamics and kinetics have been integrated by combining traditional descriptions of microbial kinetics with the equilibrium constant K and reaction quotient Q of a process, taking into account process-specific threshold energies. This approach is dynamically evolving toward a general concept of microbially driven electron transfer in anoxic environments and has been used successfully in applications ranging from bioreactor regulation to groundwater and sediment biogeochemistry.

In silico Geobacter sulfurreducens metabolism and its representation in reactive transport models

Microbial activity governs elemental cycling and the transformation of many anthropogenic substances in aqueous environments. Through the development of a dynamic cell model of the well-characterized, versatile, and abundant Geobacter sulfurreducens, we showed that a kinetic representation of key components of cell metabolism matched microbial growth dynamics observed in chemostat experiments under various environmental conditions and led to results similar to those from a comprehensive flux balance model. Coupling the kinetic cell model to its environment by expressing substrate uptake rates depending on intra- and extracellular substrate concentrations, two-dimensional reactive transport simulations of an aquifer were performed. They illustrated that a proper representation of growth efficiency as a function of substrate availability is a determining factor for the spatial distribution of microbial populations in a porous medium. It was shown that simplified model representations of microbial dynamics in the subsurface that only depended on extracellular conditions could be derived by properly parameterizing emerging properties of the kinetic cell model.

Effect of sulfate on the simultaneous bioreduction of iron and uranium

The biogeochemistry related to iron- and sulfate-reducing conditions influences the fate of contaminants such as petroleum hydrocarbons, trace metals, and radionuclides (i.e., uranium) released into the subsurface. An understanding of these processes is imperative to successfully predict the fate of contaminants during bioremediation scenarios. A series of flow-through sediment column experiments were performed to determine if the commencement of sulfate-reducing conditions would occur while bioavailable Fe(III) was present and to determine how the bioreduction of a contaminant (uranium) was affected by the switch from iron-dominated to sulfate-dominated reducing conditions. The results presented herein demonstrated that, under biostimulation, sulfate reduction can commence even though a significant pool of bioavailable Fe(III) is present. In addition, the rate of U(VI) reduction was not negatively affected by the commencement of sulfate-reducing conditions.

Mixing-controlled biodegradation in a toluene plume--results from two-dimensional laboratory experiments

Various abiotic and biotic processes such as sorption, dilution, and degradation are known to affect the fate of organic contaminants, such as petroleum hydrocarbons in saturated porous media. Reactive transport modeling of such plumes indicates that the biodegradation of organic pollutants is, in many cases, controlled by mixing and therefore occurs locally at the plume's fringes, where electron donors and electron-acceptors mix. Herein, we aim to test whether this hypothesis can be verified by experimental results obtained from aerobic and anaerobic degradation experiments in two-dimensional sediment microcosms. Toluene was selected as a model compound for oxidizable contaminants. The two-dimensional microcosm was filled with quartz sand and operated under controlled flow conditions simulating a contaminant plume in otherwise uncontaminated groundwater. Aerobic degradation of toluene by Pseudomonas putida mt-2 reduced a continuous 8.7 mg L(-1) toluene concentration by 35% over a transport distance of 78 cm in 15.5 h. In comparison, under similar conditions Aromatoleum aromaticum strain EbN1 degraded 98% of the toluene infiltrated using nitrate (68.5+/-6.2 mg L(-1)) as electron acceptor. A major part of the biodegradation activity was located at the plume fringes and the slope of the electron-acceptor gradient was steeper during periods of active biodegradation. The distribution of toluene and the significant overlap of nitrate at the plume's fringe indicate that biokinetic and/or microscale transport processes may constitute additional limiting factors. Experimental data is corroborated with results from a reactive transport model using double Monod kinetics. The outcome of the study shows that in order to simulate degradation in contaminant plumes, detailed data sets are required to test the applicability of models. These will have to deal with the incorporation of existing parameters coding for substrate conversion kinetics and microbial growth.

Microbial processes as key drivers for metal (im)mobilization along a redox gradient in the saturated zone

Two sites representing different aquifer types, i.e., Dommel (sandy) and Flémalle (gravelly loam) along the Meuse River, have been selected to conduct microcosm experiments. Various conditions ranging from aerobic over nitrate- to sulphate reducing were imposed. For the sandy aquifer, nitrate reducing conditions predominated, which specifically in the presence of a carbon source led to pH increases and enhanced Zn removal. For the calcareous gravelly loam, sulphate reduction was dominant resulting in immobilization of both Zn and Cd. For both aquifer types and almost all redox conditions, higher arsenic concentrations were measured in the groundwater. Analyses of different specific microbial populations by polymerase chain reaction (PCR) revealed the dominance of denitrifiers for the Dommel site, while sulfate reducing bacteria (SRB) were the prevailing population for all redox conditions in the Flémalle samples.

In situ hydrogen consumption kinetics as an indicator of subsurface microbial activity

There are few methods available for broadly assessing microbial community metabolism directly within a groundwater environment. In this study, hydrogen consumption rates were estimated from in situ injection/withdrawal tests conducted in two geochemically varying, contaminated aquifers as an approach towards developing such a method. The hydrogen consumption first-order rates varied from 0.002 nM h(-1) for an uncontaminated, aerobic site to 2.5 nM h(-1) for a contaminated site where sulfate reduction was a predominant process. The method could accommodate the over three orders of magnitude range in rates that existed between subsurface sites. In a denitrifying zone, the hydrogen consumption rate (0.02 nM h(-1)) was immediately abolished in the presence of air or an antibiotic mixture, suggesting that such measurements may also be sensitive to the effects of environmental perturbations on field microbial activities. Comparable laboratory determinations with sediment slurries exhibited hydrogen consumption kinetics that differed substantially from the field estimates. Because anaerobic degradation of organic matter relies on the rapid consumption of hydrogen and subsequent maintenance at low levels, such in situ measures of hydrogen turnover can serve as a key indicator of the functioning of microbial food webs and may be more reliable than laboratory determinations.

Control of hydrogen partial pressures on the rates of syntrophic microbial metabolisms: a kinetic model for butyrate fermentation

A new model describing the rate of syntrophic butyrate fermentation is constructed based on a thermodynamically consistent rate law and the metabolic pathway. This model takes into account the mechanism of reverse electron transfer and proposes that the net amount of energy saved by microorganisms as ATP depends on hydrogen partial pressures in the environment. Hydrogen partial pressures thus control not only the energy available in the environment but also the energy conserved by microorganisms. This new model predicts the rates of butyrate fermentation as a product of a kinetic factor and a thermodynamic potential factor: the kinetic factor describes how butyrate concentration controls the rates; the thermodynamic factor accounts for how the thermodynamic driving force controls the rates. Increases in hydrogen partial pressures decrease the energy available, lowering the driving force and fermentation rates. To maintain butyrate fermentation at significant rates, microorganisms decrease the amount of energy conserved, maximizing the driving force. Application of the new model demonstrates that the thermodynamic driving force is a dominant factor in controlling the rates of butyrate fermentation.

Ecological control analysis: being(s) in control of mass flux and metabolite concentrations in anaerobic degradation processes

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.

Anaerobic dechlorination and redox activities after full-scale Electrical Resistance Heating (ERH) of a TCE-contaminated aquifer

The effects of Electrical Resistance Heating (ERH) on dechlorination of TCE and redox conditions were investigated in this study. Aquifer and groundwater samples were collected prior to and after ERH treatment, where sediments were heated to approximately 100 degrees C. Sediment samples were collected from three locations and examined in microcosms for 250 to 400 days of incubation. Redox activities, in terms of consumed electron acceptors, were low in unamended microcosms with field-heated sediments, although they increased upon lactate-amendment. TCE was not dechlorinated or stalled at cDCE with field-heated sediments, which was similar or lower compared to the degree of dechlorination in unheated microcosms. However, in microcosms which were bioaugmented with a mixed anaerobic dechlorinating culture (KB-1) and lactate, dechlorination past cDCE to ethene was observed in field-heated sediments. Dechlorination and redox activities in microcosms with field-heated sediments were furthermore compared with controlled laboratory-heated microcosms, which were heated to 100 degrees C for 10 days and then slowly cooled to 10 degrees C. In laboratory-heated microcosms, TCE was not dechlorinated and redox activities remained low in unamended and lactate-amended sediments, although organic carbon was released to the aqueous phase. In contrast, in field-heated sediments, high aqueous concentrations of organic carbon were not observed in unamended microcosms, and TCE was dechlorinated to cDCE upon lactate amendment. This suggests that dechlorinating microorganisms survived the ERH or that groundwater flow through field-heated sediments carried microorganisms into the treated area and transported dissolved organic carbon downstream.

Modeling the dynamics of fermentation and respiratory processes in a groundwater plume of phenolic contaminants interpreted from laboratory- to field-scale

A biodegradation model with consecutive fermentation and respiration processes, developed from microcosm experiments and simulated mathematically with microbial growth kinetics, has been implemented into a field-scale reactive transport model of a groundwater plume of phenolic contaminants. Simulation of the anaerobic plume core with H2 and acetate as intermediate products of biodegradation allows the rates and parameter values forfermentation processes and individual respiratory terminal electron accepting processes (TEAPS) to be estimated using detailed, spatially discrete, hydrochemical field data. The modeling of field-scale plume development includes consideration of microbial acclimatization, substrate toxicity toward degradation, bioavailability of mineral oxides, and adsorption of biogenic Fe(ll) species in the aquifer, identified from complementary laboratory process studies. The results suggest that plume core processes, particularly fermentation and Fe(lll)-reduction, are more important for degradation than previously thought, possibly with a greater impact than plume fringe processes (aerobic respiration, denitrification, and SO4-reduction). The accumulation of acetate as a fermentation product within the plume contributes significantly to the mass balance for carbon. These results demonstrate the value of quantifying fermentation products within organic contaminant plumes and strongly suggest that the conceptual model selected for reactive processes plays a dominant role in the quantitative assessment of risk reduction by naturally occurring biodegradation processes.

Redox processes and release of organic matter after thermal treatment of a TCE-contaminated aquifer

Redox conditions in heated and unheated microcosm experiments were studied to evaluate the effect of thermal remediation treatment on biogeochemical processes in subsurface environments. The results were compared to field-scale observations from thermal treatments of contaminated sites. Trichloroethene-contaminated aquifer material and groundwater from Ft. Lewis, WA were incubated for 200 days at ambient temperature (i.e., 10 degrees C) or heated to 100 degrees C for 10 days and cooled slowly over a period of 150 days to mimic a thermal treatment. Increases of up to 14 mM dissolved organic carbon were observed in the aqueous phase after heating. Redox conditions did generally not change during heating in the laboratory experiment, and only minor changes occurred as an effect of heat treatment in the field. The conditions were slightly manganese/iron-reducing in two sediments and possibly sulfate-reducing in the third sediment based on production of up to 0.20 mM dissolved iron and 0.15 mM dissolved manganese and consumption of 0.08 mM sulfate. The calculated energy gain of less than -20 kJ/mol H2 for iron and sulfate reduction as well as methane production indicated that these processes were thermodynamically favorable. Sulfate reduction and methane production occurred in the unheated microcosms upon lactate amendment. Little or no reduction of the redox level was identified in heated lactate-amended microcosms, possibly because of limited microbial activity. Because the redox conditions, pH, and alkalinity remained within normal aquifer levels upon heating, bioaugmentation may be feasible for stimulating anaerobic dechlorination in heated samples or in future field applications.

Mechanisms of electron acceptor utilization: implications for simulating anaerobic biodegradation

Simulation of biodegradation reactions within a reactive transport framework requires information on mechanisms of terminal electron acceptor processes (TEAPs). In initial modeling efforts, TEAPs were approximated as occurring sequentially, with the highest energy-yielding electron acceptors (e.g. oxygen) consumed before those that yield less energy (e.g., sulfate). Within this framework in a steady state plume, sequential electron acceptor utilization would theoretically produce methane at an organic-rich source and Fe(II) further downgradient, resulting in a limited zone of Fe(II) and methane overlap. However, contaminant plumes often display much more extensive zones of overlapping Fe(II) and methane. The extensive overlap could be caused by several abiotic and biotic processes including vertical mixing of byproducts in long-screened monitoring wells, adsorption of Fe(II) onto aquifer solids, or microscale heterogeneity in Fe(III) concentrations. Alternatively, the overlap could be due to simultaneous utilization of terminal electron acceptors. Because biodegradation rates are controlled by TEAPs, evaluating the mechanisms of electron acceptor utilization is critical for improving prediction of contaminant mass losses due to biodegradation. Using BioRedox-MT3DMS, a three-dimensional, multi-species reactive transport code, we simulated the current configurations of a BTEX plume and TEAP zones at a petroleum-contaminated field site in Wisconsin. Simulation results suggest that BTEX mass loss due to biodegradation is greatest under oxygen-reducing conditions, with smaller but similar contributions to mass loss from biodegradation under Fe(III)-reducing, sulfate-reducing, and methanogenic conditions. Results of sensitivity calculations document that BTEX losses due to biodegradation are most sensitive to the age of the plume, while the shape of the BTEX plume is most sensitive to effective porosity and rate constants for biodegradation under Fe(III)-reducing and methanogenic conditions. Using this transport model, we had limited success in simulating overlap of redox products using reasonable ranges of parameters within a strictly sequential electron acceptor utilization framework. Simulation results indicate that overlap of redox products cannot be accurately simulated using the constructed model, suggesting either that Fe(III) reduction and methanogenesis are occurring simultaneously in the source area, or that heterogeneities in Fe(III) concentration and/or mineral type cause the observed overlap. Additional field, experimental, and modeling studies will be needed to address these questions.