A new rate law describing microbial respiration

Closed ecosystems extract energy through self-organized nutrient cycles

Our planet is a self-sustaining ecosystem powered by light energy from the sun, but roughly closed to matter. Many ecosystems on Earth are also approximately closed to matter and recycle nutrients by self-organizing stable nutrient cycles, e.g., microbial mats, lakes, open ocean gyres. However, existing ecological models do not exhibit the self-organization and dynamical stability widely observed in such planetary-scale ecosystems. Here, we advance a conceptual model that explains the self-organization, stability, and emergent features of closed microbial ecosystems. Our model incorporates the bioenergetics of metabolism into an ecological framework. By studying this model, we uncover a crucial thermodynamic feedback loop that enables metabolically diverse communities to almost always stabilize nutrient cycles. Surprisingly, highly diverse communities self-organize to extract [Formula: see text]10[Formula: see text] of the maximum extractable energy, or [Formula: see text]100 fold more than randomized communities. Further, with increasing diversity, distinct ecosystems show strongly correlated fluxes through nutrient cycles. However, as the driving force from light increases, the fluxes of nutrient cycles become more variable and species-dependent. Our results highlight that self-organization promotes the efficiency and stability of complex ecosystems at extracting energy from the environment, even in the absence of any centralized coordination.

Consider the Anoxic Microsite: Acknowledging and Appreciating Spatiotemporal Redox Heterogeneity in Soils and Sediments

Reduction-oxidation (redox) reactions underlie essentially all biogeochemical cycles. Like most soil properties and processes, redox is spatiotemporally heterogeneous. However, unlike other soil features, redox heterogeneity has yet to be incorporated into mainstream conceptualizations of soil biogeochemistry. Anoxic microsites, the defining feature of redox heterogeneity in bulk oxic soils and sediments, are zones of oxygen depletion in otherwise oxic environments. In this review, we suggest that anoxic microsites represent a critical component of soil function and that appreciating anoxic microsites promises to advance our understanding of soil and sediment biogeochemistry. In sections 1 and 2, we define anoxic microsites and highlight their dynamic properties, specifically anoxic microsite distribution, redox gradient magnitude, and temporality. In section 3, we describe the influence of anoxic microsites on several key elemental cycles, organic carbon, nitrogen, iron, manganese, and sulfur. In section 4, we evaluate methods for identifying and characterizing anoxic microsites, and in section 5, we highlight past and current approaches to modeling anoxic microsites. Finally, in section 6, we suggest steps for incorporating anoxic microsites and redox heterogeneities more broadly into our understanding of soils and sediments.

Functional convergence in slow-growing microbial communities arises from thermodynamic constraints

The dynamics of microbial communities is complex, determined by competition for metabolic substrates and cross-feeding of byproducts. Species in the community grow by harvesting energy from chemical reactions that transform substrates to products. In many anoxic environments, these reactions are close to thermodynamic equilibrium and growth is slow. To understand the community structure in these energy-limited environments, we developed a microbial community consumer-resource model incorporating energetic and thermodynamic constraints on an interconnected metabolic network. The central element of the model is product inhibition, meaning that microbial growth may be limited not only by depletion of metabolic substrates but also by accumulation of products. We demonstrate that these additional constraints on microbial growth cause a convergence in the structure and function of the community metabolic network-independent of species composition and biochemical details-providing a possible explanation for convergence of community function despite taxonomic variation observed in many natural and industrial environments. Furthermore, we discovered that the structure of community metabolic network is governed by the thermodynamic principle of maximum free energy dissipation. Our results predict the decrease of functional convergence in faster growing communities, which we validate by analyzing experimental data from anaerobic digesters. Overall, the work demonstrates how universal thermodynamic principles may constrain community metabolism and explain observed functional convergence in microbial communities.

Sulfate reduction accelerates groundwater arsenic contamination even in aquifers with abundant iron oxides

Groundwater contamination by geogenic arsenic is a global problem affecting nearly 200 million people. In South and Southeast Asia, a cost-effective mitigation strategy is to use oxidized low-arsenic aquifers rather than reduced high-arsenic aquifers. Aquifers with abundant oxidized iron minerals are presumably safeguarded against immediate arsenic contamination, due to strong sorption of arsenic onto iron minerals. However, preferential pumping of low-arsenic aquifers can destabilize the boundaries between these aquifers, pulling high-arsenic water into low-arsenic aquifers. We investigate this scenario in a hybrid field-column experiment in Bangladesh where naturally high-arsenic groundwater is pumped through sediment cores from a low-arsenic aquifer, and detailed aqueous and solid-phase measurements are used to constrain reactive transport modelling. Here we show that elevated groundwater arsenic concentrations are induced by sulfate reduction and the predicted formation of highly mobile, poorly sorbing thioarsenic species. This process suggests that contamination of currently pristine aquifers with arsenic can occur up to over 1.5 times faster than previously thought, leading to a deterioration of urgently needed water resources.

Redox Properties of Solid Phase Electron Acceptors Affect Anaerobic Microbial Respiration under Oxygen-Limited Conditions in Floodplain Soils

Mountain floodplain soils often show spatiotemporal variations in redox conditions that arise due to changing hydrology and resulting biogeochemistry. Under oxygen-depleted conditions, solid phase terminal electron acceptors (TEAs) can be used in anaerobic respiration. However, it remains unclear to what degree the redox properties of solid phases limit respiration rates and hence organic matter degradation. Here, we assess such limitations in soils collected across a gradient in native redox states from the Slate River floodplain (Colorado, U.S.A.). We incubated soils under anoxic conditions and quantified CO₂ production and microbial Fe(III) reduction, the main microbial metabolic pathway, as well as the reactivity of whole-soil solid phase TEAs toward mediated electrochemical reduction. Fe(III) reduction occurred together with CO₂ production in native oxic soils, while neither Fe(II) nor CO₂ production was observed in native anoxic soils. Initial CO₂ production rates increased with increasing TEA redox reactivity toward mediated electrochemical reduction across all soil depths. Low TEA redox reactivity appears to be caused by elevated Fe(II) concentrations rather than crystallinity of Fe(III) phases. Our findings illustrate that the buildup of Fe(II) in systems with long residence times limits the thermodynamic viability of dissimilatory Fe(III) reduction and thereby limits the mineralization of organic carbon.

A mechanistic derivation of the Monod bioreaction equation for a limiting nutrient

We present a quasi-steady state mechanistic derivation of the Monod bioreaction equation based upon a conceptual model involving aqueous phase diffusive transport of substrate towards a spherical microbe; transport of the substrate across its surface membrane; and reaction depleting the substrate within the microbe. The resulting Monod coefficients [Formula: see text] and [Formula: see text] are dependent upon substrate-species pairs and the mass transfer properties of the system. Two substrate transport scenarios are investigated: (1) a constant rate model that is a function of a constant flux across the surface of the microbe; and (2) a linear rate model that is the product of a constant transport velocity and the concentration of substrate in contact with the surface of the microbe. The model is verified and parameterized using benzene, toluene, and phenol depletion and biomass growth data obtained from Reardon et al. (Biotechnol Bioeng: 385-400, 2000). Calibration results indicate a normalized surface to bulk concentration ratio of nearly unity in all simulations for benzene, toluene, and phenol when paired with P. putida F1, implying that the process is not aqueous phase diffusion limited.

Physiological Acclimation Extrapolates the Kinetics and Thermodynamics of Methanogenesis From Laboratory Experiments to Natural Environments

Chemotrophic microorganisms face the steep challenge of limited energy resources in natural environments. This observation has important implications for interpreting and modeling the kinetics and thermodynamics of microbial reactions. Current modeling frameworks treat microbes as autocatalysts, and simulate microbial energy conservation and growth with fixed kinetic and thermodynamic parameters. However, microbes are capable of acclimating to the environment and modulating their parameters in order to gain competitive fitness. Here we constructed an optimization model and described microbes as self-adapting catalysts by linking microbial parameters to intracellular metabolic resources. From the optimization results, we related microbial parameters to the substrate concentration and the energy available in the environment, and simplified the relationship between the kinetics and the thermodynamics of microbial reactions. We took as examples Methanosarcina and Methanosaeta – the methanogens that produce methane from acetate – and showed how the acclimation model extrapolated laboratory observations to natural environments and improved the simulation of methanogenesis and the dominance of Methanosaeta over Methanosarcina in lake sediments. These results highlight the importance of physiological acclimation in shaping the kinetics and thermodynamics of microbial reactions and in determining the outcome of microbial interactions.

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

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

Temperature dependence of nitrate-reducing Fe(II) oxidation by Acidovorax strain BoFeN1 - evaluating the role of enzymatic vs. abiotic Fe(II) oxidation by nitrite

Fe(II) oxidation coupled to nitrate reduction is a widely observed metabolism. However, to what extent the observed Fe(II) oxidation is driven enzymatically or abiotically by metabolically produced nitrite remains puzzling. To distinguish between biotic and abiotic reactions, we cultivated the mixotrophic nitrate-reducing Fe(II)-oxidizing Acidovorax strain BoFeN1 over a wide range of temperatures and compared it to abiotic Fe(II) oxidation by nitrite at temperatures up to 60°C. The collected experimental data were subsequently analyzed through biogeochemical modeling. At 5°C, BoFeN1 cultures consumed acetate and reduced nitrate but did not significantly oxidize Fe(II). Abiotic Fe(II) oxidation by nitrite at different temperatures showed an Arrhenius-type behavior with an activation energy of 80±7 kJ/mol. Above 40°C, the kinetics of Fe(II) oxidation were abiotically driven, whereas at 30°C, where BoFeN1 can actively metabolize, the model-based interpretation strongly suggested that an enzymatic pathway was responsible for a large fraction (ca. 62%) of the oxidation. This result was reproduced even when no additional carbon source was present. Our results show that at below 30°C, i.e., at temperatures representing most natural environments, biological Fe(II) oxidation was largely responsible for overall Fe(II) oxidation, while abiotic Fe(II) oxidation by nitrite played a less important role.

Importance of thermodynamics dependent kinetic parameters in nitrate-based souring mitigation studies

Souring is the unwanted formation of hydrogen sulfide (H₂S) by sulfate-reducing microorganisms (SRM) in sewer systems and seawater flooded oil reservoirs. Nitrate treatment (NT) is one of the major methods to alleviate souring: The mechanism of souring remediation by NT is stimulation of nitrate reducing microorganisms (NRM) that depending on the nitrate reduction pathway can outcompete SRM for common electron donors, or oxidize sulfide to sulfate. However, some nitrate reduction pathways may challenge the efficacy of NT. Therefore, a precise understanding of souring rate, nitrate reduction rate and pathways is crucial for efficient souring management. Here, we investigate the necessity of incorporating two thermodynamic dependent kinetic parameters, namely, the growth yield (Y), and F_T, a parameter related to the minimum catabolic energy production required by cells to utilize a given catabolic reaction. We first show that depending on physiochemical conditions, Y and F_T for SRM change significantly in the range of [0-0.4] mole biomass per mole electron donor and [0.0006-0.5], respectively, suggesting that these parameters should not be considered constant and that it is important to couple souring models with thermodynamic models. Then, we highlight this further by showing an experimental dataset that can be modeled very well by considering variable F_T. Next, we show that nitrate based lithotrophic sulfide oxidation to sulfate (lNRM₃) is the dominant nitrate reduction pathway. Then, arguing that thermodynamics would suggest that S° consumption should proceed faster than S⁰ production, we infer that the reason for frequently observed S⁰ accumulation is its low solubility. Last, we suggest that nitrate based souring treatment will suffer less from S⁰ accumulation if we (i) act early, (ii) increase temperature and (iii) supplement stoichiometrically sufficient nitrate.

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.

The view of microbes as energy converters illustrates the trade-off between growth rate and yield

The application of thermodynamics to microbial growth has a long tradition that originated in the middle of the 20th century. This approach reflects the view that self-replication is a thermodynamic process that is not fundamentally different from mechanical thermodynamics. The key distinction is that a free energy gradient is not converted into mechanical (or any other form of) energy but rather into new biomass. As such, microbes can be viewed as energy converters that convert a part of the energy contained in environmental nutrients into chemical energy that drives self-replication. Before the advent of high-throughput sequencing technologies, only the most central metabolic pathways were known. However, precise measurement techniques allowed for the quantification of exchanged extracellular nutrients and heat of growing microbes with their environment. These data, together with the absence of knowledge of metabolic details, drove the development of so-called black-box models, which only consider the observable interactions of a cell with its environment and neglect all details of how exactly inputs are converted into outputs. Now, genome sequencing and genome-scale metabolic models (GEMs) provide us with unprecedented detail about metabolic processes inside the cell. However, mostly due to computational complexity issues, the derived modelling approaches make surprisingly little use of thermodynamic concepts. Here, we review classical black-box models and modern approaches that integrate thermodynamics into GEMs. We also illustrate how the description of microbial growth as an energy converter can help to understand and quantify the trade-off between microbial growth rate and yield.

The Energetic Potential for Undiscovered Manganese Metabolisms in Nature

Microorganisms are found in nearly every surface and near-surface environment, where they gain energy by catalyzing reactions among a wide variety of chemical compounds. The discovery of new catabolic strategies and microbial habitats can therefore be guided by determining which redox reactions can supply energy under environmentally-relevant conditions. In this study, we have explored the thermodynamic potential of redox reactions involving manganese, one of the most abundant transition metals in the Earth's crust. In particular, we have assessed the Gibbs energies of comproportionation and disproportionation reactions involving Mn²⁺ and several Mn-bearing oxide and oxyhydroxide minerals containing Mn in the +II, +III, and +IV oxidation states as a function of temperature (0-100°C) and pH (1-13). In addition, we also calculated the energetic potential of Mn²⁺ oxidation coupled to O₂, NO₂ ^-, NO₃ ^-, and FeOOH. Results show that these reactions-none of which, except O₂ + Mn²⁺, are known catabolisms-can provide energy to microorganisms, particularly at higher pH values and temperatures. Comproportionation between Mn²⁺ and pyrolusite, for example, can yield 10 s of kJ (mol Mn)^-1. Disproportionation of Mn³⁺ can yield more than 100 kJ (mol Mn)^-1 at conditions relevant to natural settings such as sediments, ferromanganese nodules and crusts, bioreactors and suboxic portions of the water column. Of the Mn²⁺ oxidation reactions, the one with nitrite as the electron acceptor is most energy yielding under most combinations of pH and temperature. We posit that several Mn redox reactions represent heretofore unknown microbial metabolisms.

Controls on Interspecies Electron Transport and Size Limitation of Anaerobically Methane-Oxidizing Microbial Consortia

About 382 Tg yr⁻¹ of methane rising through the seafloor is oxidized anaerobically (W. S. Reeburgh, Chem Rev 107:486–513, 2007, https://doi.org/10.1021/cr050362v), preventing it from reaching the atmosphere, where it acts as a strong greenhouse gas. Microbial consortia composed of anaerobic methanotrophic archaea and sulfate-reducing bacteria couple the oxidation of methane to the reduction of sulfate under anaerobic conditions via a syntrophic process. Recent experimental studies and modeling efforts indicate that direct interspecies electron transfer (DIET) is involved in this syntrophy. Here, we explore a fluorescent in situ hybridization-nanoscale secondary ion mass spectrometry data set of large, segregated anaerobic oxidation of methane (AOM) consortia that reveal a decline in metabolic activity away from the archaeal-bacterial interface and use a process-based model to identify the physiological controls on rates of AOM. Simulations reproducing the observational data reveal that ohmic resistance and activation loss are the two main factors causing the declining metabolic activity, where activation loss dominated at a distance of <8 μm. These voltage losses limit the maximum spatial distance between syntrophic partners with model simulations, indicating that sulfate-reducing bacterial cells can remain metabolically active up to ∼30 μm away from the archaeal-bacterial interface. Model simulations further predict that a hybrid metabolism that combines DIET with a small contribution of diffusive exchange of electron donors can offer energetic advantages for syntrophic consortia.

Sulfate-dependent reversibility of intracellular reactions explains the opposing isotope effects in the anaerobic oxidation of methane

The anaerobic oxidation of methane (AOM) is performed by methanotrophic archaea (ANME) in distinct sulfate-methane interfaces of marine sediments. In these interfaces, AOM often appears to deplete methane in the heavy isotopes toward isotopic compositions similar to methanogenesis. Here, we shed light on this effect and its physiological underpinnings using a thermophilic ANME-1-dominated culture. At high sulfate concentrations, residual methane is enriched in both ¹³C and ²H (¹³α = 1.016 and ²α = 1.155), as observed previously. In contrast, at low sulfate concentrations, the residual methane is substantially depleted in ¹³C (¹³α = 0.977) and, to a lesser extent, in ²H. Using a biochemical-isotopic model, we explain the sulfate dependence of the net isotopic fractionation through the thermodynamic drive of the involved intracellular reactions. Our findings relate these isotopic patterns to the physiology and environment of the ANME, thereby explaining a commonly observed isotopic enigma.

Overlapping redox zones control arsenic pollution in Pleistocene multi-layer aquifers, the Po Plain (Italy)

Understanding the factors that control As concentrations in groundwater is vital for supplying safe groundwater in regions with As-polluted aquifers. Despite much research, mainly addressing Holocene aquifers hosting young (<100 yrs) groundwater, the source, transport, and fate of As in Pleistocene aquifers with fossil (>12,000 yrs) groundwaters are not yet fully understood and so are assessed here through an evaluation of the redox properties of the system in a type locality, the Po Plain (Italy). Analyses of redox-sensitive species and major ions on 22 groundwater samples from the Pleistocene arsenic-affected aquifer in the Po Plain shows that groundwater concentrations of As are controlled by the simultaneous operation of several terminal electron accepters. Organic matter, present as peat, is abundant in the aquifer, allowing groundwater to reach a quasi-steady-state of highly reducing conditions close to thermodynamic equilibrium. In this system, simultaneous reduction of Fe-oxide and sulfate results in low concentrations of As (median 7 μg/L) whereas As reaches higher concentrations (median of 82 μg/L) during simultaneous methanogenesis and Fe-reduction. The position of well-screens is an additional controlling factor on groundwater As: short screens that overlap confining aquitards generate higher As concentrations than long screens placed away from them. A conceptual model for groundwater As, applicable worldwide in other Pleistocene aquifers with reducible Fe-oxides and abundant organic matter is proposed: As may have two concentration peaks, the first after prolonged Fe-oxide reduction and until sulfate reduction takes place, the second during simultaneous Fe-reduction and methanogenesis.

Thermodynamic Inhibition in a Biofilm Reactor with Suspended Bacteria

We formulate a biofilm reactor model with suspended bacteria that accounts for thermodynamic growth inhibition. The reactor model is a chemostat style model consisting of a single replenished growth promoting substrate, a single reaction product, suspended bacteria, and wall attached bacteria in the form of a bacterial biofilm. We present stability conditions for the washout equilibrium using standard techniques, demonstrating that analytical results are attainable even with the added complexity from thermodynamic inhibition. Furthermore, we numerically investigate the longterm behaviour. In the computational study, we investigate model behaviour for select parameters and two commonly used detachment functions. We investigate the effects of thermodynamic inhibition on the model and find that thermodynamic inhibition limits substrate utilization/production both inside the biofilm and inside the aqueous phase, resulting in less suspended bacteria and a thinner biofilm.

Phytoplankton Temporal Strategies Increase Entropy Production in a Marine Food Web Model

We develop a trait-based model founded on the hypothesis that biological systems evolve and organize to maximize entropy production by dissipating chemical and electromagnetic free energy over longer time scales than abiotic processes by implementing temporal strategies. A marine food web consisting of phytoplankton, bacteria, and consumer functional groups is used to explore how temporal strategies, or the lack thereof, change entropy production in a shallow pond that receives a continuous flow of reduced organic carbon plus inorganic nitrogen and illumination from solar radiation with diel and seasonal dynamics. Results show that a temporal strategy that employs an explicit circadian clock produces more entropy than a passive strategy that uses internal carbon storage or a balanced growth strategy that requires phytoplankton to grow with fixed stoichiometry. When the community is forced to operate at high specific growth rates near 2 d⁻¹, the optimization-guided model selects for phytoplankton ecotypes that exhibit complementary for winter versus summer environmental conditions to increase entropy production. We also present a new type of trait-based modeling where trait values are determined by maximizing entropy production rather than by random selection.

Not Just Numbers: Mathematical Modelling and Its Contribution to Anaerobic Digestion Processes

Mathematical modelling of bioprocesses has a long and notable history, with eminent contributions from fields including microbiology, ecology, biophysics, chemistry, statistics, control theory and mathematical theory. This richness of ideas and breadth of concepts provide great motivation for inquisitive engineers and intrepid scientists to try their hand at modelling, and this collaboration of disciplines has also delivered significant milestones in the quality and application of models for both theoretical and practical interrogation of engineered biological systems. The focus of this review is the anaerobic digestion process, which, as a technology that has come in and out of fashion, remains a fundamental process for addressing the global climate emergency. Whether with conventional anaerobic digestion systems, biorefineries, or other anaerobic technologies, mathematical models are important tools that are used to design, monitor, control and optimise the process. Both highly structured, mechanistic models and data-driven approaches have been used extensively over half a decade, but recent advances in computational capacity, scientific understanding and diversity and quality of process data, presents an opportunity for the development of new modelling paradigms, augmentation of existing methods, or even incorporation of tools from other disciplines, to ensure that anaerobic digestion research can remain resilient and relevant in the face of emerging and future challenges.

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.

Thermodynamic modelling of synthetic communities predicts minimum free energy requirements for sulfate reduction and methanogenesis

Microbial communities are complex dynamical systems harbouring many species interacting together to implement higher-level functions. Among these higher-level functions, conversion of organic matter into simpler building blocks by microbial communities underpins biogeochemical cycles and animal and plant nutrition, and is exploited in biotechnology. A prerequisite to predicting the dynamics and stability of community-mediated metabolic conversions is the development and calibration of appropriate mathematical models. Here, we present a generic, extendable thermodynamic model for community dynamics and calibrate a key parameter of this thermodynamic model, the minimum energy requirement associated with growth-supporting metabolic pathways, using experimental population dynamics data from synthetic communities composed of a sulfate reducer and two methanogens. Our findings show that accounting for thermodynamics is necessary in capturing the experimental population dynamics of these synthetic communities that feature relevant species using low energy growth pathways. Furthermore, they provide the first estimates for minimum energy requirements of methanogenesis (in the range of −30 kJ mol⁻¹) and elaborate on previous estimates of lactate fermentation by sulfate reducers (in the range of −30 to −17 kJ mol⁻¹ depending on the culture conditions). The open-source nature of the developed model and demonstration of its use for estimating a key thermodynamic parameter should facilitate further thermodynamic modelling of microbial communities.

Consistent microbial dynamics and functional community patterns derived from first principles

Microbial communities are key engines that drive earth's biogeochemical cycles. However, existing ecosystem models have only limited ability to predict microbial dynamics and require the calibration of multiple population-specific empirical equations. In contrast, we build on a new kinetic "Microbial Transition State" (MTS) theory of growth derived from first principles. We show how the theory coupled to simple mass and energy balance calculations provides a framework with intrinsically important qualitative properties to model microbial community dynamics. We first show how the theory can simultaneously account for the influence of all the resources needed for growth (electron donor, acceptor, and nutrients) while still producing consistent dynamics that fulfill the Liebig rule of a single limiting substrate. We also show consistent patterns of energy-dependent microbial successions in mixed culture without the need for calibration of population-specific parameters. We then show how this approach can be used to model a simplified activated sludge community. To this end, we compare MTS-derived dynamics with those of a widely used activated sludge model and show that similar growth yields and overall dynamics can be obtained using two parameters instead of twelve. This new kinetic theory of growth grounded by a set of generic physical principles parsimoniously gives rise to consistent microbial population and community dynamics, thereby paving the way for the development of a new class of more predictive microbial ecosystem models.

Identifying and Quantifying the Intermediate Processes during Nitrate-Dependent Iron(II) Oxidation

Microbially driven nitrate-dependent iron (Fe) oxidation (NDFO) in subsurface environments has been intensively studied. However, the extent to which Fe(II) oxidation is biologically catalyzed remains unclear because no neutrophilic iron-oxidizing and nitrate reducing autotroph has been isolated to confirm the existence of an enzymatic pathway. While mixotrophic NDFO bacteria have been isolated, understanding the process is complicated by simultaneous abiotic oxidation due to nitrite produced during denitrification. In this study, the relative contributions of biotic and abiotic processes during NDFO were quantified through the compilation and model-based interpretation of previously published experimental data. The kinetics of chemical denitrification by Fe(II) (chemodenitrification) were assessed, and compelling evidence was found for the importance of organic ligands, specifically exopolymeric substances secreted by bacteria, in enhancing abiotic oxidation of Fe(II). However, nitrite alone could not explain the observed magnitude of Fe(II) oxidation, with 60-75% of overall Fe(II) oxidation attributed to an enzymatic pathway for investigated strains: Acidovorax ( A.) strain BoFeN1, 2AN, A. ebreus strain TPSY, Paracoccus denitrificans Pd 1222, and Pseudogulbenkiania sp. strain 2002. By rigorously quantifying the intermediate processes, this study eliminated the potential for abiotic Fe(II) oxidation to be exclusively responsible for NDFO and verified the key contribution from an additional, biological Fe(II) oxidation process catalyzed by NDFO bacteria.

Dynamics of soil biogeochemical gas emissions shaped by remolded aggregate sizes and carbon configurations under hydration cycles

Changes in soil hydration status affect microbial community dynamics and shape key biogeochemical processes. Evidence suggests that local anoxic conditions may persist and support anaerobic microbial activity in soil aggregates (or in similar hot spots) long after the bulk soil becomes aerated. To facilitate systematic studies of interactions among environmental factors with biogeochemical emissions of CO₂ , N₂ O and CH₄ from soil aggregates, we remolded silt soil aggregates to different sizes and incorporated carbon at different configurations (core, mixed, no addition). Assemblies of remolded soil aggregates of three sizes (18, 12, and 6 mm) and equal volumetric proportions were embedded in sand columns at four distinct layers. The water table level in each column varied periodically while obtaining measurements of soil GHG emissions for the different aggregate carbon configurations. Experimental results illustrate that methane production required prolonged inundation and highly anoxic conditions for inducing measurable fluxes. The onset of unsaturated conditions (lowering water table) resulted in a decrease in CH₄ emissions while temporarily increasing N₂ O fluxes. Interestingly, N₂ O fluxes were about 80% higher form aggregates with carbon placement in center (anoxic) core compared to mixed carbon within aggregates. The fluxes of CO₂ were comparable for both scenarios of carbon sources. These experimental results highlight the importance of hydration dynamics in activating different GHG production and affecting various transport mechanisms about 80% of total methane emissions during lowering water table level are attributed to physical storage (rather than production), whereas CO₂ emissions (~80%) are attributed to biological activity. A biophysical model for microbial activity within soil aggregates and profiles provides a means for results interpretation and prediction of trends within natural soils under a wide range of conditions.

Anaerobic microsites have an unaccounted role in soil carbon stabilization

Soils represent the largest carbon reservoir within terrestrial ecosystems. The mechanisms controlling the amount of carbon stored and its feedback to the climate system, however, remain poorly resolved. Global carbon models assume that carbon cycling in upland soils is entirely driven by aerobic respiration; the impact of anaerobic microsites prevalent even within well-drained soils is missed within this conception. Here, we show that anaerobic microsites are important regulators of soil carbon persistence, shifting microbial metabolism to less efficient anaerobic respiration, and selectively protecting otherwise bioavailable, reduced organic compounds such as lipids and waxes from decomposition. Further, shifting from anaerobic to aerobic conditions leads to a 10-fold increase in volume-specific mineralization rate, illustrating the sensitivity of anaerobically protected carbon to disturbance. The vulnerability of anaerobically protected carbon to future climate or land use change thus constitutes a yet unrecognized soil carbon–climate feedback that should be incorporated into terrestrial ecosystem models.

Mechanisms controlling soil carbon storage and feedbacks to the climate system remain poorly constrained. Here, the authors show that anaerobic microsites stabilize soil carbon by shifting microbial metabolism to less efficient anaerobic respiration and protecting reduced organic compounds from decomposition.

Atmospheric CH4 oxidation by Arctic permafrost and mineral cryosols as a function of water saturation and temperature

The response of methanotrophic bacteria capable of oxidizing atmospheric CH₄ to climate warming is poorly understood, especially for those present in Arctic mineral cryosols. The atmospheric CH₄ oxidation rates were measured in microcosms incubated at 4 °C and 10 °C along a 1-m depth profile and over a range of water saturation conditions for mineral cryosols containing type I and type II methanotrophs from Axel Heiberg Island (AHI), Nunavut, Canada. The cryosols exhibited net consumption of ~2 ppmv CH₄ under all conditions, including during anaerobic incubations. Methane oxidation rates increased with temperature and decreased with increasing water saturation and depth, exhibiting the highest rates at 10 °C and 33% saturation at 5 cm depth (260 ± 60 pmol CH₄ gdw^-1 d^-1 ). Extrapolation of the CH₄ oxidation rates to the field yields net CH₄ uptake fluxes ranging from 11 to 73 μmol CH₄  m^-2 d^-1 , which are comparable to field measurements. Stable isotope mass balance indicates ~50% of the oxidized CH₄ is incorporated into the biomass regardless of temperature or saturation. Future atmospheric CH₄ uptake rates at AHI with increasing temperatures will be determined by the interplay of increasing CH₄ oxidation rates vs. water saturation and the depth to the water table during summer thaw.

An oligotrophic deep-subsurface community dependent on syntrophy is dominated by sulfur-driven autotrophic denitrifiers

Subsurface lithoautotrophic microbial ecosystems (SLiMEs) under oligotrophic conditions are typically supported by H₂ Methanogens and sulfate reducers, and the respective energy processes, are thought to be the dominant players and have been the research foci. Recent investigations showed that, in some deep, fluid-filled fractures in the Witwatersrand Basin, South Africa, methanogens contribute <5% of the total DNA and appear to produce sufficient CH₄ to support the rest of the diverse community. This paradoxical situation reflects our lack of knowledge about the in situ metabolic diversity and the overall ecological trophic structure of SLiMEs. Here, we show the active metabolic processes and interactions in one of these communities by combining metatranscriptomic assemblies, metaproteomic and stable isotopic data, and thermodynamic modeling. Dominating the active community are four autotrophic β-proteobacterial genera that are capable of oxidizing sulfur by denitrification, a process that was previously unnoticed in the deep subsurface. They co-occur with sulfate reducers, anaerobic methane oxidizers, and methanogens, which each comprise <5% of the total community. Syntrophic interactions between these microbial groups remove thermodynamic bottlenecks and enable diverse metabolic reactions to occur under the oligotrophic conditions that dominate in the subsurface. The dominance of sulfur oxidizers is explained by the availability of electron donors and acceptors to these microorganisms and the ability of sulfur-oxidizing denitrifiers to gain energy through concomitant S and H₂ oxidation. We demonstrate that SLiMEs support taxonomically and metabolically diverse microorganisms, which, through developing syntrophic partnerships, overcome thermodynamic barriers imposed by the environmental conditions in the deep subsurface.

Thermodynamic and Kinetic Response of Microbial Reactions to High CO2

Geological carbon sequestration captures CO₂ from industrial sources and stores the CO₂ in subsurface reservoirs, a viable strategy for mitigating global climate change. In assessing the environmental impact of the strategy, a key question is how microbial reactions respond to the elevated CO₂ concentration. This study uses biogeochemical modeling to explore the influence of CO₂ on the thermodynamics and kinetics of common microbial reactions in subsurface environments, including syntrophic oxidation, iron reduction, sulfate reduction, and methanogenesis. The results show that increasing CO₂ levels decreases groundwater pH and modulates chemical speciation of weak acids in groundwater, which in turn affect microbial reactions in different ways and to different extents. Specifically, a thermodynamic analysis shows that increasing CO₂ partial pressure lowers the energy available from syntrophic oxidation and acetoclastic methanogenesis, but raises the available energy of microbial iron reduction, hydrogenotrophic sulfate reduction and methanogenesis. Kinetic modeling suggests that high CO₂ has the potential of inhibiting microbial sulfate reduction while promoting iron reduction. These results are consistent with the observations of previous laboratory and field studies, and highlight the complexity in microbiological responses to elevated CO₂ abundance, and the potential power of biogeochemical modeling in evaluating and quantifying these responses.

Diurnal dynamics of oxygen and carbon dioxide concentrations in shoots and rhizomes of a perennial in a constructed wetland indicate down-regulation of below ground oxygen consumption

Wetland plants actively provide oxygen for aerobic processes in submerged tissues and the rhizosphere. The novel concomitant assessment of diurnal dynamics of oxygen and carbon dioxide concentrations under field conditions tests the whole-system interactions in plant-internal gas exchange and regulation. Oxygen concentrations ([O2]) were monitored in-situ in central culm and rhizome pith cavities of common reed (Phragmites australis) using optical oxygen sensors. The corresponding carbon dioxide concentrations ([CO2]) were assessed via gas samples from the culms. Highly dynamic diurnal courses of [O2] were recorded, which started at 6.5-13 % in the morning, increased rapidly up to 22 % during midday and declined exponentially during the night. Internal [CO2] were high in the morning (1.55-17.5 %) and decreased (0.04-0.94 %) during the rapid increase of [O2] in the culms. The observed negative correlations between [O2] and [CO2] particularly describe the below ground relationship between plant-mediated oxygen supply and oxygen use by respiration and biogeochemical processes in the rhizosphere. Furthermore, the nocturnal declining slopes of [O2] in culms and rhizomes indicated a down-regulation of the demand for oxygen in the complete below ground plant-associated system. These findings emphasize the need for measurements of plant-internal gas exchange processes under field conditions because it considers the complex interactions in the oxic-anoxic interface.

Microbial diversity arising from thermodynamic constraints

The microbial world displays an immense taxonomic diversity. This diversity is manifested also in a multitude of metabolic pathways that can utilise different substrates and produce different products. Here, we propose that these observations directly link to thermodynamic constraints that inherently arise from the metabolic basis of microbial growth. We show that thermodynamic constraints can enable coexistence of microbes that utilise the same substrate but produce different end products. We find that this thermodynamics-driven emergence of diversity is most relevant for metabolic conversions with low free energy as seen for example under anaerobic conditions, where population dynamics is governed by thermodynamic effects rather than kinetic factors such as substrate uptake rates. These findings provide a general understanding of the microbial diversity based on the first principles of thermodynamics. As such they provide a thermodynamics-based framework for explaining the observed microbial diversity in different natural and synthetic environments.

The Thermodynamics of Marine Biogeochemical Cycles: Lotka Revisited

Nearly 100 years ago, Alfred Lotka published two short but insightful papers describing how ecosystems may organize. Principally, Lotka argued that ecosystems will grow in size and that their cycles will spin faster via predation and nutrient recycling so as to capture all available energy, and that evolution and natural selection are the mechanisms by which this occurs and progresses. Lotka's ideas have often been associated with the maximum power principle, but they are more consistent with recent developments in nonequilibrium thermodynamics, which assert that complex systems will organize toward maximum entropy production (MEP). In this review, we explore Lotka's hypothesis within the context of the MEP principle, as well as how this principle can be used to improve marine biogeochemistry models. We need to develop the equivalent of a climate model, as opposed to a weather model, to understand marine biogeochemistry on longer timescales, and adoption of the MEP principle can help create such models.

Power limits for microbial life

To better understand the origin, evolution, and extent of life, we seek to determine the minimum flux of energy needed for organisms to remain viable. Despite the difficulties associated with direct measurement of the power limits for life, it is possible to use existing data and models to constrain the minimum flux of energy required to sustain microorganisms. Here, a we apply a bioenergetic model to a well characterized marine sedimentary environment in order to quantify the amount of power organisms use in an ultralow-energy setting. In particular, we show a direct link between power consumption in this environment and the amount of biomass (cells cm^-3) found in it. The power supply resulting from the aerobic degradation of particular organic carbon (POC) at IODP Site U1370 in the South Pacific Gyre is between ∼10^-12 and 10^-16 W cm^-3. The rates of POC degradation are calculated using a continuum model while Gibbs energies have been computed using geochemical data describing the sediment as a function of depth. Although laboratory-determined values of maintenance power do a poor job of representing the amount of biomass in U1370 sediments, the number of cells per cm^-3 can be well-captured using a maintenance power, 190 zW cell^-1, two orders of magnitude lower than the lowest value reported in the literature. In addition, we have combined cell counts and calculated power supplies to determine that, on average, the microorganisms at Site U1370 require 50–3500 zW cell^-1, with most values under ∼300 zW cell^-1. Furthermore, we carried out an analysis of the absolute minimum power requirement for a single cell to remain viable to be on the order of 1 zW cell^-1.

Distribution of iron- and sulfate-reducing bacteria across a coastal acid sulfate soil (CASS) environment: implications for passive bioremediation by tidal inundation

Coastal acid sulfate soils (CASS) constitute a serious and global environmental problem. Oxidation of iron sulfide minerals exposed to air generates sulfuric acid with consequently negative impacts on coastal and estuarine ecosystems. Tidal inundation represents one current treatment strategy for CASS, with the aim of neutralizing acidity by triggering microbial iron- and sulfate-reduction and inducing the precipitation of iron-sulfides. Although well-known functional guilds of bacteria drive these processes, their distributions within CASS environments, as well as their relationships to tidal cycling and the availability of nutrients and electron acceptors, are poorly understood. These factors will determine the long-term efficacy of "passive" CASS remediation strategies. Here we studied microbial community structure and functional guild distribution in sediment cores obtained from 10 depths ranging from 0 to 20 cm in three sites located in the supra-, inter- and sub-tidal segments, respectively, of a CASS-affected salt marsh (East Trinity, Cairns, Australia). Whole community 16S rRNA gene diversity within each site was assessed by 454 pyrotag sequencing and bioinformatic analyses in the context of local hydrological, geochemical, and lithological factors. The results illustrate spatial overlap, or close association, of iron-, and sulfate-reducing bacteria (SRB) in an environment rich in organic matter and controlled by parameters such as acidity, redox potential, degree of water saturation, and mineralization. The observed spatial distribution implies the need for empirical understanding of the timing, relative to tidal cycling, of various terminal electron-accepting processes that control acid generation and biogeochemical iron and sulfur cycling.

Thermodynamic controls on the kinetics of microbial low-pH Fe(II) oxidation

Acid mine drainage (AMD) is a major worldwide environmental threat to surface and groundwater quality. Microbial low-pH Fe(II) oxidation could be exploited for cost-effective AMD treatment; however, its use is limited because of uncertainties associated with its rate and ability to remove Fe from solution. We developed a thermodynamic-based framework to evaluate the kinetics of low-pH Fe(II) oxidation. We measured the kinetics of low-pH Fe(II) oxidation at five sites in the Appalachian Coal Basin in the US and three sites in the Iberian Pyrite Belt in Spain and found that the fastest rates of Fe(II) oxidation occurred at the sites with the lowest pH values. Thermodynamic calculations showed that the Gibbs free energy of Fe(II) oxidation (ΔG(oxidation)) was also most negative at the sites with the lowest pH values. We then conducted two series of microbial Fe(II) oxidation experiments in laboratory-scale chemostatic bioreactors operated through a series of pH values (2.1-4.2) and found the same relationships between Fe(II) oxidation kinetics, ΔG(oxidation), and pH. Conditions that favored the fastest rates of Fe(II) oxidation coincided with higher Fe(III) solubility. The solubility of Fe(III) minerals, thus plays an important role on Fe(II) oxidation kinetics. Methods to incorporate microbial low-pH Fe(II) oxidation into active and passive AMD treatment systems are discussed in the context of these findings. This study presents a simplified model that describes the relationship between free energy and microbial kinetics and should be broadly applicable to many biogeochemical systems.

The potential for biologically catalyzed anaerobic methane oxidation on ancient Mars

This study examines the potential for the biologically mediated anaerobic oxidation of methane (AOM) coupled to sulfate reduction on ancient Mars. Seven distinct fluids representative of putative martian groundwater were used to calculate Gibbs energy values in the presence of dissolved methane under a range of atmospheric CO2 partial pressures. In all scenarios, AOM is exergonic, ranging from -31 to -135 kJ/mol CH4. A reaction transport model was constructed to examine how environmentally relevant parameters such as advection velocity, reactant concentrations, and biomass production rate affect the spatial and temporal dependences of AOM reaction rates. Two geologically supported models for ancient martian AOM are presented: a sulfate-rich groundwater with methane produced from serpentinization by-products, and acid-sulfate fluids with methane from basalt alteration. The simulations presented in this study indicate that AOM could have been a feasible metabolism on ancient Mars, and fossil or isotopic evidence of this metabolic pathway may persist beneath the surface and in surface exposures of eroded ancient terrains.

Microbial activity in the marine deep biosphere: progress and prospects

The vast marine deep biosphere consists of microbial habitats within sediment, pore waters, upper basaltic crust and the fluids that circulate throughout it. A wide range of temperature, pressure, pH, and electron donor and acceptor conditions exists—all of which can combine to affect carbon and nutrient cycling and result in gradients on spatial scales ranging from millimeters to kilometers. Diverse and mostly uncharacterized microorganisms live in these habitats, and potentially play a role in mediating global scale biogeochemical processes. Quantifying the rates at which microbial activity in the subsurface occurs is a challenging endeavor, yet developing an understanding of these rates is essential to determine the impact of subsurface life on Earth's global biogeochemical cycles, and for understanding how microorganisms in these “extreme” environments survive (or even thrive). Here, we synthesize recent advances and discoveries pertaining to microbial activity in the marine deep subsurface, and we highlight topics about which there is still little understanding and suggest potential paths forward to address them. This publication is the result of a workshop held in August 2012 by the NSF-funded Center for Dark Energy Biosphere Investigations (C-DEBI) “theme team” on microbial activity (www.darkenergybiosphere.org).

Development of microbial-enzyme-mediated decomposition model parameters through steady-state and dynamic analyses

We developed a microbial-enzyme-mediated decomposition (MEND) model, based on the Michaelis-Menten kinetics, that describes the dynamics of physically defined pools of soil organic matter (SOC). These include particulate, mineral-associated, dissolved organic matter (POC, MOC, and DOC, respectively), microbial biomass, and associated exoenzymes. The ranges and/or distributions of parameters were determined by both analytical steady-state and dynamic analyses with SOC data from the literature. We used an improved multi-objective parameter sensitivity analysis (MOPSA) to identify the most important parameters for the full model: maintenance of microbial biomass, turnover and synthesis of enzymes, and carbon use efficiency (CUE). The model predicted that an increase of 2 degrees C (baseline temperature 12 degrees C) caused the pools of POC-cellulose, MOC, and total SOC to increase with dynamic CUE and decrease with constant CUE, as indicated by the 50% confidence intervals. Regardless of dynamic or constant CUE, the changes in pool size of POC, MOC, and total SOC varied from -8% to 8% under +2 degrees C. The scenario analysis using a single parameter set indicates that higher temperature with dynamic CUE might result in greater net increases in both POC-cellulose and MOC pools. Different dynamics of various SOC pools reflected the catalytic functions of specific enzymes targeting specific substrates and the interactions between microbes, enzymes, and SOC. With the feasible parameter values estimated in this study, models incorporating fundamental principles of microbial-enzyme dynamics can lead to simulation results qualitatively different from traditional models with fast/slow/passive pools.

A theoretical reassessment of microbial maintenance and implications for microbial ecology modeling

We attempted to reconcile three microbial maintenance models (Herbert, Pirt, and Compromise) through a theoretical reassessment. We provided a rigorous proof that the true growth yield coefficient (Y(G)) is the ratio of the specific maintenance rate (a in Herbert) to the maintenance coefficient (m in Pirt). Other findings from this study include: (1) the Compromise model is identical to the Herbert for computing microbial growth and substrate consumption, but it expresses the dependence of maintenance on both microbial biomass and substrate; (2) the maximum specific growth rate in the Herbert (μ(max,H)) is higher than those in the other two models (μ(max,P) and μ(max,C)), and the difference is the physiological maintenance factor (m(q) = a); and (3) the overall maintenance coefficient (m(T)) is more sensitive to m(q) than to the specific growth rate (μ(G)) and Y(G). Our critical reassessment of microbial maintenance provides a new approach for quantifying some important components in soil microbial ecology models.

Modeling of mitochondria bioenergetics using a composable chemiosmotic energy transduction rate law: theory and experimental validation

Mitochondrial bioenergetic processes are central to the production of cellular energy, and a decrease in the expression or activity of enzyme complexes responsible for these processes can result in energetic deficit that correlates with many metabolic diseases and aging. Unfortunately, existing computational models of mitochondrial bioenergetics either lack relevant kinetic descriptions of the enzyme complexes, or incorporate mechanisms too specific to a particular mitochondrial system and are thus incapable of capturing the heterogeneity associated with these complexes across different systems and system states. Here we introduce a new composable rate equation, the chemiosmotic rate law, that expresses the flux of a prototypical energy transduction complex as a function of: the saturation kinetics of the electron donor and acceptor substrates; the redox transfer potential between the complex and the substrates; and the steady-state thermodynamic force-to-flux relationship of the overall electro-chemical reaction. Modeling of bioenergetics with this rate law has several advantages: (1) it minimizes the use of arbitrary free parameters while featuring biochemically relevant parameters that can be obtained through progress curves of common enzyme kinetics protocols; (2) it is modular and can adapt to various enzyme complex arrangements for both in vivo and in vitro systems via transformation of its rate and equilibrium constants; (3) it provides a clear association between the sensitivity of the parameters of the individual complexes and the sensitivity of the system's steady-state. To validate our approach, we conduct in vitro measurements of ETC complex I, III, and IV activities using rat heart homogenates, and construct an estimation procedure for the parameter values directly from these measurements. In addition, we show the theoretical connections of our approach to the existing models, and compare the predictive accuracy of the rate law with our experimentally fitted parameters to those of existing models. Finally, we present a complete perturbation study of these parameters to reveal how they can significantly and differentially influence global flux and operational thresholds, suggesting that this modeling approach could help enable the comparative analysis of mitochondria from different systems and pathological states. The procedures and results are available in Mathematica notebooks at http://www.igb.uci.edu/tools/sb/mitochondria-modeling.html.

Variation in energy available to populations of subsurface anaerobes in response to geological carbon storage

Microorganisms can strongly influence the chemical and physical properties of the subsurface. Changes in microbial activity caused by geological CO(2) storage, therefore, have the potential to influence the capacity, injectivity, and integrity of CO(2) storage reservoirs and ultimately the environmental impact of CO(2) injection. This analysis uses free energy calculations to examine variation in energy available to Fe(III) and SO(4)(2-) reducers and methanogens because of changes in the bulk composition of brine and shallow groundwater following subsurface CO(2) injection. Calculations were performed using data from two field experiments, the Frio Formation experiment and an experiment at the Zero Emission Research and Technology test site. Energy available for Fe(III) reduction increased significantly during CO(2) injection in both experiments, largely because of a decrease in pH from near-neutral levels to just below 6. Energy available to SO(4)(2-) reducers and methanogens varied little. These changes can lead to a greater rate of microbial Fe(III) reduction following subsurface CO(2) injection in reservoirs where Fe(III) oxides or oxyhydroxides are available and the rate of Fe(III) reduction is limited by energy available prior to injection.