Description: The fluids emanating from active submarine hydrothermal vent chimneys provide a window into subseafloor processes and, through mixing with seawater, are responsible for steep thermal and compositional gradients that provide the energetic basis for diverse biological communities. Although several models have been developed to better understand the dynamic interplay of seawater, hydrothermal fluid, minerals and microorganisms inside chimney walls, none provide a fully integrated approach to quantifying the biogeochemistry of these hydrothermal systems. In an effort to remedy this, a fully coupled biogeochemical reaction-transport model of a hydrothermal vent chimney has been developed that explicitly quantifies the rates of microbial catalysis while taking into account geochemical processes such as fluid flow, solute transport and oxidation–reduction reactions associated with fluid mixing as a function of temperature. The metabolisms included in the reaction network are methanogenesis, aerobic oxidation of hydrogen, sulfide and methane and sulfate reduction by hydrogen and methane. Model results indicate that microbial catalysis is generally fastest in the hottest habitable portion of the vent chimney (77–102 °C), and methane and sulfide oxidation peak near the seawater-side of the chimney. The fastest metabolisms are aerobic oxidation of H2 and sulfide and reduction of sulfate by H2 with maximum rates of 140, 900 and 800 pmol cm−3 d−1, respectively. The maximum rate of hydrogenotrophic methanogenesis is just under 0.03 pmol cm−3 d−1, the slowest of the metabolisms considered. Due to thermodynamic inhibition, there is no anaerobic oxidation of methane by sulfate (AOM). These simulations are consistent with vent chimney metabolic activity inferred from phylogenetic data reported in the literature. The model developed here provides a quantitative approach to describing the rates of biogeochemical transformations in hydrothermal systems and can be used to constrain the role of microbial activity in the deep subsurface.

Description: This study investigates the controls on organic carbon and molybdenum (Mo) accumulation in sediments deposited within the Western Interior Seaway across the Cenomanian–Turonian boundary interval (94.34–93.04 Ma) including Oceanic Anoxic Event 2 (OAE2). Carbon fluxes to the sediment–water interface (reflecting changes in primary productivity) and bottom-water oxygen concentrations (reflecting preservation effects) are reconstructed from field data and used to constrain a benthic model that simulates the geochemistry of unconsolidated sediments as they were deposited. The results show that increased availability of reactive iron prevents Mo sequestration as thiomolybdate (MoS42 −) during OAE2 (O2 ~ 105 μM) by (i) inhibiting sulfate reduction, and (ii) buffering any free sulfide that becomes available. In the post-OAE2 period (O2 ~ 50 μM), Mo accumulation is favored by a large reduction in iron flux. Importantly, this occurs in parallel with oxygenated bottom waters and high rates of aerobic carbon degradation in the surface sediments, implying that elevated Mo burial fluxes in ancient marine facies do not necessarily reflect euxinic or even anoxic conditions within the water column. Our findings suggest that both an increase in production and preservation lead to enrichment in organic carbon in the Western Interior Seaway. More generally, the results demonstrate that a careful consideration of the coupling between iron, carbon and oxygen cycles during the early stages of diagenesis is critical for interpreting geochemical proxies in modern and ancient settings.

Description: We present sedimentary geochemical data and in situ benthic flux measurements of dissolved inorganic nitrogen (DIN: NO3−, NO2−, NH4+) and oxygen (O2) from 7 sites with variable sand content along 18°N offshore Mauritania (NW Africa). Bottom water O2 concentrations at the shallowest station were hypoxic (42 μM) and increased to 125 μM at the deepest site (1113 m). Total oxygen uptake rates were highest on the shelf (−10.3 mmol O2 m−2 d−1) and decreased quasi-exponentially with water depth to −3.2 mmol O2 m−2 d−1. Average denitrification rates estimated from a flux balance decreased with water depth from 2.2 to 0.2 mmol N m−2 d−1. Overall, the sediments acted as net sink for DIN. Observed increases in δ15NNO3 and δ18ONO3 in the benthic chamber deployed on the shelf, characterized by muddy sand, were used to calculate apparent benthic nitrate fractionation factors of 8.0‰ (15εapp) and 14.1‰ (18εapp). Measurements of δ15NNO2 further demonstrated that the sediments acted as a source of 15N depleted NO2−. These observations were analyzed using an isotope box model that considered denitrification and nitrification of NH4+ and NO2−. The principal findings were that (i) net benthic 14N/15N fractionation (εDEN) was 12.9 ± 1.7‰, (ii) inverse fractionation during nitrite oxidation leads to an efflux of isotopically light NO2− (−22 ± 1.9‰), and (iii) direct coupling between nitrification and denitrification in the sediment is negligible. Previously reported εDEN for fine-grained sediments are much lower (4–8‰). We speculate that high benthic nitrate fractionation is driven by a combination of enhanced porewater–seawater exchange in permeable sediments and the hypoxic, high productivity environment. Although not without uncertainties, the results presented could have important implications for understanding the current state of the marine N cycle.

Description: Two ∼6 m long sediment cores were collected along the ∼300 m isobath on the Alaskan Beaufort Sea continental margin. Both cores showed distinct sulfate-methane transition zones (SMTZ) at 105 and 120 cm below seafloor (cmbsf). Sulfate was not completely depleted below the SMTZ but remained between 30 and 500 μM. Sulfate reduction and anaerobic oxidation of methane (AOM) determined by radiotracer incubations were active throughout the methanogenic zone. Although a mass balance could not explain the source of sulfate below the SMTZ, geochemical profiles and correlation network analyses of biotic and abiotic data suggest a cryptic sulfur cycle involving iron, manganese and barite. Inhibition experiments with molybdate and 2-bromoethanesulfonate (BES) indicated decoupling of sulfate reduction and AOM and competition between sulfate reducers and methanogens for substrates. While correlation network analyses predicted coupling of AOM to iron reduction, the addition of manganese or iron did not stimulate AOM. Since none of the classical archaeal anaerobic methanotrophs (ANME) were abundant, the involvement of unknown or unconventional phylotypes in AOM is conceivable. The resistance of AOM activity to inhibitors implies deviation from conventional enzymatic pathways. This work suggests that the classical redox cascade of electron acceptor utilization based on Gibbs energy yields does not always hold in diffusion-dominated systems, and instead biotic processes may be more strongly coupled to mineralogy.

Description: Quantification of global biogeochemical cycles requires knowledge of the rates at which microorganisms catalyze chemical reactions. In order for models that describe these processes to capture global patterns of change, the underlying formulations in them must account for biogeochemical transformations over seasonal and millennial time scales in environments characterized by different energy levels. Building on existing models, a new thermodynamic limiting function is introduced. With only one adjustable parameter, this function that can be used to model microbial metabolism throughout the range of conditions in which organisms are known to be active. The formulation is based on a comparison of the amount of energy available from any redox reaction to the energy required to maintain a membrane potential, a proxy for the minimum amount of energy required by an active microorganism. This function does not require species- or metabolism-specific parameters, and can be used to model metabolisms that capture any amount of energy. The utility of this new thermodynamic rate limiting term is illustrated by applying it to three low-energy processes: fermentation, methanogenesis and sulfate reduction. The model predicts that the rate of fermentation will be reduced by half once the Gibbs energy of the catalyzed reaction reaches −12 kJ (mol e−)−1, and then slowing exponentially until the energy yield approaches zero. Similarly, the new model predicts that the low energy yield of methanogenesis, −4 to −0.5 kJ (mol e−)−1, for a partial pressure of H2 between 11 and 0.6 Pa decreases the reaction rate by 95–99%. Finally, the new function’s utility is illustrated through its ability to accurately model sulfate concentration data in an anoxic marine sediment.