Description: The kinetics of particulate organic carbon (POC) mineralization in marine surface sediments is not well constrained. This creates considerable uncertainties when benthic processes are considered in global biogeochemical or Earth system circulation models to simulate climate-ocean interactions and biogeochemical tracers in the ocean. In an attempt to improve our understanding of the rate and depth distribution of organic carbon mineralization in bioturbated (0–10 cm) sediments, we parameterized a 1-D diagenetic model that simulates the reactivity of three discrete POC pools at global scale (a "multi-G" model). The rate constants of the three reactive classes (highly reactive, reactive, refractory) are fixed and determined to be 70 yr−1, 0.5 yr−1, and ~0.001 yr−1, respectively, based on the Martin curve model for pelagic POC degradation. In contrast to previous approaches, the reactivity of the organic material degraded in the seafloor is continuous with, and set by, the apparent reactivity of material sinking through the water column. The model is able to simulate a global database (185 stations) of benthic oxygen and nitrate fluxes across the sediment-water interface in addition to porewater oxygen and nitrate distributions and organic carbon burial efficiencies. It is further consistent with degradation experiments of fresh phytoplankton. We propose that an important yet mostly overlooked consideration in previous upscaling approaches is the proportion of the relative reactive POC classes reaching the seafloor in addition to their reactivity. The approach presented is applicable to both steady-state and non-steady state scenarios, and links POC degradation kinetics in sedimentary environments to water depth and the POC rain rate to the seafloor.

Description: The benthic environment is a crucial component of marine systems in the provision of ecosystem services, sustaining biodiversity and in climate regulation, and therefore important to human society. With the contemporary increase in computational power, model resolution and technological improvements in quality and quantity of benthic data, it is necessary to ensure that benthic systems are appropriately represented in coupled benthic-pelagic biogeochemical and ecological modelling studies. In this paper we focus on five topical challenges related to various aspects of modelling benthic environments: organic matter reactivity, dynamics of benthic-pelagic boundary layer, microphytobenthos, biological transport and small-scale heterogeneity, and impacts of episodic events. We discuss current gaps in their understanding and indicate plausible ways ahead. Further, we propose a three-pronged approach for the advancement of benthic and benthic-pelagic modelling, essential for improved understanding, management and prediction of the marine environment. This includes: (A) development of a traceable and hierarchical framework for benthic-pelagic models, which will facilitate integration among models, reduce risk of bias, and clarify model limitations; (B) extended cross-disciplinary approach to promote effective collaboration between modelling and empirical scientists of various backgrounds and better involvement of stakeholders and end-users; (C) a common vocabulary for terminology used in benthic modelling, to promote model development and integration, and also to enhance mutual understanding.

Description: Oxygen minimum zones (OMZs) that impinge on continental margins favor the release of phosphorus (P) from the sediments to the water column, enhancing primary productivity and the maintenance or expansion of low-oxygen waters. A comprehensive field program in the Peruvian OMZ was undertaken to identify the sources of benthic P at six stations, including the analysis of particles from the water column, surface sediments, and pore fluids, as well as in situ benthic flux measurements. A major fraction of solid-phase P was bound as particulate inorganic P (PIP) both in the water column and in sediments. Sedimentary PIP increased with depth in the sediment at the expense of particulate organic P (POP). The ratio of particulate organic carbon (POC) to POP exceeded the Redfield ratio both in the water column (202 ± 29) and in surface sediments (303 ± 77). However, the POC to total particulate P (TPP = POP + PIP) ratio was close to Redfield in the water column (103 ± 9) and in sediment samples (102 ± 15). This suggests that the relative burial efficiencies of POC and TPP are similar under low-oxygen conditions and that the sediments underlying the anoxic waters on the Peru margin are not depleted in P compared to Redfield. Benthic fluxes of dissolved P were extremely high (up to 1.04 ± 0.31 mmol m−2 d−1), however, showing that a lack of oxygen promotes the intensified release of dissolved P from sediments, whilst preserving the POC / TPP burial ratio. Benthic dissolved P fluxes were always higher than the TPP rain rate to the seabed, which is proposed to be caused by transient P release by bacterial mats that had stored P during previous periods when bottom waters were less reducing. At one station located at the lower rim of the OMZ, dissolved P was taken up by the sediments, indicating ongoing phosphorite formation. This is further supported by decreasing porewater phosphate concentrations with sediment depth, whereas solid-phase P concentrations were comparatively high.

Description: Porewater profiles in sediment cores from mangrove-dominated coastal lagoons (Celestún and Chelem) on the Yucatán Peninsula, Mexico, reveal the widespread coexistence of dissolved methane and sulfate. This observation is interesting since dissolved methane in porewaters is typically oxidized anaerobically by sulfate. To explain the observations we used a numerical transport-reaction model that was constrained by the field observations. The model suggests that methane in the upper sediments is produced in the sulfate reduction zone at rates ranging between 0.012 and 31 mmolm-2 d-1, concurrent with sulfate reduction rates between 1.1 and 24 mmol SO2- 4 m-2 d-1. These processes are supported by high organic matter content in the sediment and the use of non-competitive substrates by methanogenic microorganisms. Indeed sediment slurry incubation experiments show that non-competitive substrates such as trimethylamine (TMA) and methanol can be utilized for microbial methanogenesis at the study sites. The model also indicates that a significant fraction of methane is transported to the sulfate reduction zone from deeper zones within the sedimentary column by rising bubbles and gas dissolution. The shallow depths of methane production and the fast rising methane gas bubbles reduce the likelihood for oxidation, thereby allowing a large fraction of the methane formed in the sediments to escape to the overlying water column.

Description: Highlights • An empirical diagenetic model includes redox-dependent P storage by microorganisms. • Sediment mixing and burrowing by animals strongly enhances P burial. • Supporting evidence for a decrease in oceanic P inventory in the early Paleozoic. Abstract A diagenetic model is used to simulate the diagenesis and burial of particulate organic carbon (Corg) and phosphorus (P) in marine sediments underlying anoxic versus oxic bottom waters. The latter are physically mixed by animals moving through the surface sediment (bioturbation) and ventilated by burrowing, tube-dwelling organisms (bioirrigation). The model is constrained using an empirical database including burial ratios of Corg with respect to organic P (Corg:Porg) and total reactive P (Corg:Preac), burial efficiencies of Corg and Porg, and inorganic carbon-to-phosphorus regeneration ratios. If Porg is preferentially mineralized relative to Corg during aerobic respiration, as many previous studies suggest, then the simulated Porg pool is found to be completely depleted. A modified model that incorporates the redox-dependent microbial synthesis of polyphosphates and Porg (termed the microbial P pump) allows preferential mineralization of the bulk Porg pool relative to Corg during both aerobic and anaerobic respiration and is consistent with the database. Results with this model show that P burial is strongly enhanced in sediments hosting fauna. Animals mix highly labile Porg away from the aerobic sediment layers where mineralization rates are highest, thereby mitigating diffusive PO43− fluxes to the bottom water. They also expand the redox niche where microbial P uptake occurs. The model was applied to a hypothetical shelf setting in the early Paleozoic; a time of the first radiation of benthic fauna. Results show that even shallow bioturbation at that time may have had a significant impact on P burial. Our model provides support for a recent study that proposed that faunal radiation in ocean sediments led to enhanced P burial and, possibly, a stabilization of atmospheric O2 levels. The results also help to explain Corg:Porg ratios in the geological record and the persistence of Porg in ancient marine sediments.

Description: A kinetic-bioenergetic reaction model for the anaerobic oxidation of methane (AOM) in coastal marine sediments is presented. The model considers a fixed depth interval of sediments below the zone of bioturbation (the window-of-observation), subject to seasonal variations of temperature and inputs of organic substrates and sulfate. It includes (1) nine microbially-mediated reaction pathways involved in CH4 production/consumption; (2) an explicit representation of five functional microbial groups; and (3) bioenergetic limitations of the microbial metabolic pathways. Fermentation of organic substrates is assumed to produce hydrogen (H2) and acetate (Ac) as key reactive intermediates. Competition among the metabolic pathways is controlled by the relative kinetic efficiencies of the various microbial processes and by bioenergetic constraints. Model results imply that the functional microbial biomasses within the window-of-observation undergo little variation over the year, as a result of kinetic and thermodynamic buffering of the seasonal forcings. Furthermore, the microbial processes proceed at only small fractions of their maximum potential rates. These findings provide a theoretical justification for the approximation of steady-state microbial biomasses, which is frequently used in diagenetic modeling. In contrast, AOM rates show a strong seasonal evolution: AOM only becomes spontaneous in winter, when hydrogenotrophic sulfate reduction (hySR) sufficiently reduces the local H2 concentration. The bioenergetic limitation of AOM is thus a critical factor modulating this process in seasonally-forced nearshore marine sediments. A global sensitivity analysis based on a 2-level factorial design reveals that AOM rates are most sensitive to the kinetic parameters describing hySR and acetotrophic methanogenesis (acME). The growth and substrate uptake kinetics of AOM are unimportant, whereas the threshold value of ATP energy conservation for AOM is the most sensitive thermodynamic parameter. These results confirm that anaerobic methane oxidizing microorganisms are metabolizing close to their thermodynamic limit, with the energetic balance being controlled by the relative rates of hySR and acME. The removal of Ac by acME primarily allows more sulfate (SO42−) to be utilized for H2 oxidation, thereby promoting AOM.

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: Sediment incubation microcosms, multitrap apparatus and water column variables have been employed to describe the dynamic changes in benthic-pelagic coupling between nutrient pools in the Pontevedra ria, NW Spain, during spring and summer 1998. A comparison of the chemical characteristics of suspended and bed sediment together with sediment nutrient effluxes revealed that particulate organic nitrogen and carbon were progressively depleted upon transit through the nutrient pools. The main fate of particulate organic nitrogen reaching the bed sediment is denitrification, although resuspension is also important. An estimate of net denitrification (Dnet) was made at 2 muddy sites in the ria with a mass balance at the benthic boundary layer. First-order approximations calculate Dnet to be 178 and 182 µmolN m-2 h-1 at the 2 stations, and agree well with previous modelling estimates. Denitrification is highest when upwelling relaxes and the flux of organic matter to the sediment increases. Regular inputs of offshore water ensure water renewal and re-oxygenation of bottom waters, thus preventing anoxia, particularly in the summer. With upwelling, large quantities of ammonium are effluxed to the water column (250 µmolNH4+ m-2 h-1), probably as a result of the bed sediment resuspension engendered by upwelling and stirring of phytodetrital fluff held in suspension as neutrally buoyant material above the sediment surface. We hypothesise that hydrodynamical processes play an important role in determining the quantity of nutrients remineralised in the Pontevedra ria and, in the case of nitrogen, the rate of denitrification at the benthic boundary layer.

Description: Geochemical data (CH4, SO42−, I−, Cl−, particulate organic carbon (POC), δ13C-CH4, and δ13C-CO2) are presented from the upper 30 m of marine sediment on a tectonic submarine accretionary wedge offshore southwest Taiwan. The sampling stations covered three ridges (Tai-Nan, Yung-An, and Good Weather), each characterized by bottom simulating reflectors, acoustic turbidity, and different types of faulting and anticlines. Sulfate and iodide concentrations varied little from seawater-like values in the upper 1–3 m of sediment at all stations; a feature that is consistent with irrigation of seawater by gas bubbles rising through the soft surface sediments. Below this depth, sulfate was rapidly consumed within 5–10 m by anaerobic oxidation of methane (AOM) at the sulfate-methane transition. Carbon isotopic data imply a mainly biogenic methane source. A numerical transport-reaction model was used to identify the supply pathways of methane and estimate depth-integrated turnover rates at the three ridges. Methane gas ascending from deep layers, facilitated by thrusts and faults, was by far the dominant term in the methane budget at all sites. Differences in the proximity of the sampling sites to the faults and anticlines mainly accounted for the variability in gas fluxes and depth-integrated AOM rates. By comparison, methane produced in situ by POC degradation within the modeled sediment column was unimportant. This study demonstrates that the geochemical trends in the continental margins offshore SW Taiwan are closely related to the different geological settings.

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.