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: 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: The accumulation of gas hydrates in marine sediments is essentially controlled by the accumulation of particulate organic carbon (POC) which is microbially converted into methane, the thickness of the gas hydrate stability zone (GHSZ) where methane can be trapped, the sedimentation rate (SR) that controls the time that POC and the generated methane stays within the GHSZ, and the delivery of methane from deep-seated sediments by ascending pore fluids and gas into the GHSZ. Recently, Wallmann et al. (2012) presented transfer functions to predict the gas hydrate inventory in diffusion-controlled geological systems based on SR, POC and GHSZ thickness for two different scenarios: normal and full compacting sediments. We apply these functions to global data sets of bathymetry, heat flow, seafloor temperature, POC input and SR, estimating a global mass of carbon stored in marine methane hydrates from 3 to 455 Gt of carbon (GtC) depending on the sedimentation and compaction conditions. The global sediment volume of the GHSZ in continental margins is estimated to be 60–67 × 1015 m3, with a total of 7 × 1015 m3 of pore volume (available for GH accumulation). However, seepage of methane-rich fluids is known to have a pronounced effect on gas hydrate accumulation. Therefore, we carried out a set of systematic model runs with the transport-reaction code in order to derive an extended transfer function explicitly considering upward fluid advection. Using averaged fluid velocities for active margins, which were derived from mass balance considerations, this extended transfer function predicts the enhanced gas hydrate accumulation along the continental margins worldwide. Different scenarios were investigated resulting in a global mass of sub-seafloor gas hydrates of ~ 550 GtC. Overall, our systematic approach allows to clearly and quantitatively distinguish between the effect of biogenic methane generation from POC and fluid advection on the accumulation of gas hydrate, and hence, provides a simple prognostic tool for the estimation of large-scale and global gas hydrate inventories in marine sediments.

Description: This study presents benthic data from 12 samplings from February to December 2010 in a 28 m deep channel in the southwest Baltic Sea. In winter, the distribution of solutes in the porewater was strongly modulated by bioirrigation which efficiently flushed the upper 10 cm of sediment, leading to concentrations which varied little from bottom water values. Solute pumping by bioirrigation fell sharply in the summer as the bottom waters became severely hypoxic (〈 2 μM O2). At this point the giant sulfide-oxidizing bacteria Beggiatoa was visible on surface sediments. Despite an increase in O2 following mixing of the water column in November, macrofauna remained absent until the end of the sampling. Contrary to expectations, metabolites such as dissolved inorganic carbon, ammonium and hydrogen sulfide did not accumulate in the upper 10 cm during the hypoxic period when bioirrigation was absent, but instead tended toward bottom water values. This was taken as evidence for episodic bubbling of methane gas out of the sediment acting as an abiogenic irrigation process. Porewater–seawater mixing by escaping bubbles provides a pathway for enhanced nutrient release to the bottom water and may exacerbate the feedback with hypoxia. Subsurface dissolved phosphate (TPO4) peaks in excess of 400 μM developed in autumn, resulting in a very large diffusive TPO4 flux to the water column of 0.7 ± 0.2 mmol m−2 d−1. The model was not able to simulate this TPO4 source as release of iron-bound P (Fe–P) or organic P. As an alternative hypothesis, the TPO4 peak was reproduced using new kinetic expressions that allow Beggiatoa to take up porewater TPO4 and accumulate an intracellular P pool during periods with oxic bottom waters. TPO4 is then released during hypoxia, as previous published results with sulfide-oxidizing bacteria indicate. The TPO4 added to the porewater over the year by organic P and Fe–P is recycled through Beggiatoa, meaning that no additional source of TPO4 is needed to explain the TPO4 peak. Further experimental studies are needed to strengthen this conclusion and rule out Fe–P and organic P as candidate sources of ephemeral TPO4 release. A measured C/P ratio of 〈 20 for the diffusive flux to the water column during hypoxia directly demonstrates preferential release of P relative to C under oxygen-deficient bottom waters. This coincides with a strong decrease in dissolved inorganic N/P ratios in the water column to ~ 1. Our results suggest that sulfide oxidizing bacteria could act as phosphorus capacitors in systems with oscillating redox conditions, releasing massive amounts of TPO4 in a short space of time and dramatically increasing the internal loading of TPO4 to the overlying water.

Description: A simple prognostic tool for gas hydrate (GH) quantification in marine sediments is presented based on a diagenetic transport-reaction model approach. One of the most crucial factors for the application of diagenetic models is the accurate formulation of microbial degradation rates of particulate organic carbon (POC) and the coupled formation of biogenic methane. Wallmann et al. (2006) suggested a kinetic formulation considering the ageing effects of POC and accumulation of reaction products (CH4, CO2) in the pore water. This model is applied to data sets of several ODP sites in order to test its general validity. Based on a thorough parameter analysis considering a wide range of environmental conditions, the POC accumulation rate (POCar in g/m2/yr) and the thickness of the gas hydrate stability zone (GHSZ in m) were identified as the most important and independent controls for biogenic GH formation. Hence, depth-integrated GH inventories in marine sediments (GHI in g of CH4 per cm2 seafloor area) can be estimated as: GHI=a ·POCar·GHSZb ·exp(−GHSZc/POCar/d)+e with a = 0.00214, b = 1.234, c = −3.339, d = 0.3148, e = −10.265. The transfer function gives a realistic first order approximation of the minimum GH inventory in low gas flux (LGF) systems. The overall advantage of the presented function is its simplicity compared to the application of complex numerical models, because only two easily accessible parameters need to be determined.

Description: Carbon cycling in Peruvian margin sediments (11° S and 12° S) was examined at 16 stations from 74 m on the inner shelf down to 1024 m water depth by means of in situ flux measurements, sedimentary geochemistry and modeling. Bottom water oxygen was below detection limit down to ca. 400 m and increased to 53 μM at the deepest station. Sediment accumulation rates and benthic dissolved inorganic carbon fluxes decreased rapidly with water depth. Particulate organic carbon (POC) content was lowest on the inner shelf and at the deep oxygenated stations (〈 5%) and highest between 200 and 400 m in the oxygen minimum zone (OMZ, 15–20%). The organic carbon burial efficiency (CBE) was unexpectedly low on the inner shelf (〈 20%) when compared to a global database, for reasons which may be linked to the frequent ventilation of the shelf by oceanographic anomalies. CBE at the deeper oxygenated sites was much higher than expected (max. 81%). Elsewhere, CBEs were mostly above the range expected for sediments underlying normal oxic bottom waters, with an average of 51 and 58% for the 11° S and 12° S transects, respectively. Organic carbon rain rates calculated from the benthic fluxes alluded to a very efficient mineralization of organic matter in the water column, with a Martin curve exponent typical of normal oxic waters (0.88 ± 0.09). Yet, mean POC burial rates were 2–5 times higher than the global average for continental margins. The observations at the Peruvian margin suggest that a lack of oxygen does not affect the degradation of organic matter in the water column but promotes the preservation of organic matter in marine sediments.

Description: We have developed and employed an Earth system model to explore the forcings of atmospheric pCO2 change and the chemical and isotopic evolution of seawater over the last glacial cycle. Concentrations of dissolved phosphorus (DP), reactive nitrogen, molecular oxygen, dissolved inorganic carbon (DIC), total alkalinity (TA), 13C-DIC, and 14CDIC were calculated for 24 ocean boxes. The bi-directional water fluxes between these model boxes were derived from a 3-D circulation field of the modern ocean (Opa 8.2, NEMO) and tuned such that tracer distributions calculated by the box model were consistent with observational data from the modern ocean. To model the last 130 kyr, we employed records of past changes in sea-level, ocean circulation, and dust deposition. According to the model, about half of the glacial pCO2 drawdown may be attributed to marine regressions. The glacial sea-level low-stands implied steepened ocean margins, a reduced burial of particulate organic carbon, phosphorus, and neritic carbonate at the margin seafloor, a decline in benthic denitrification, and enhanced weathering of emerged shelf sediments. In turn, low-stands led to a distinct rise in the standing stocks of DIC, TA, and nutrients in the global ocean, promoted the glacial sequestration of atmospheric CO2 in the ocean, and added 13C- and 14C-depleted DIC to the ocean as recorded in benthic foraminifera signals. The other half of the glacial drop in pCO2 was linked to inferred shoaling of Atlantic meridional overturning circulation and more efficient utilization of nutrients in the Southern Ocean. The diminished ventilation of deep water in the glacial Atlantic and Southern Ocean led to significant 14C depletions with respect to the atmosphere. According to our model, the deglacial rapid and stepwise rise in atmospheric pCO2 was induced by upwelling both in the Southern Ocean and subarctic North Pacific and promoted by a drop in nutrient utilization in the Southern Ocean. The deglacial sea-level rise led to a gradual decline in nutrient, DIC, and TA stocks, a slow change due to the large size and extended residence times of dissolved chemical species in the ocean. Thus, the rapid deglacial rise in pCO2 can be explained by fast changes in ocean dynamics and nutrient utilization whereas the gradual pCO2 rise over the Holocene may be linked to the slow drop in nutrient and TA stocks that continued to promote an ongoing CO2 transfer from the ocean into the atmosphere.

Description: The eastern tropical South Pacific (ETSP) represents one of the most productive areas in the ocean that is characterized by a pronounced oxygen minimum zone (OMZ). Particulate organic matter (POM) that sinks out of the euphotic zone is supplied to the anoxic sediments and utilized by microbial communities. The degradation of POM is associated with dissolved organic matter (DOM) production and reworking. The release of recalcitrant DOM to the overlying waters may represent an important organic matter escape mechanism from remineralization within sediments but received little attention in OMZ regions so far. Here, we combine measurements of dissolved organic carbon (DOC) and dissolved organic nitrogen (DON) with DOM optical properties in the form of chromophoric (CDOM) and fluorescent (FDOM) DOM from pore waters and near-bottom waters of the ETSP off Peru. We evaluate diffusion–driven fluxes and net in situ fluxes of DOC and DON in order to investigate processes affecting DOM cycling at the sediment–water interface along a transect 12° S. To our knowledge, these are the first data for sediment release of DON and pore water CDOM and FDOM for the ETSP off Peru. Pore-water DOC and DON accumulated with increasing sediment depth, suggesting an imbalance between DOM production and remineralization within sediments. High DON accumulation resulted in very low pore water DOC / DON ratios (〉 1) which could be caused by either an "imbalance" in DOC and DON remineralization, or to the presence of an additional nitrogen source. Diffusion driven fluxes of DOC and DON exhibited high spatial variability. They varied from 0.2–0.1 mmol m−2 d−1 to 2.52–1.3 mmol m−2 d−1 and from −0.042–0.02 mmol m−2 d−1 to 3.32–1.7 mmol m−2 d−1, respectively. Generally low net in situ DOC and DON fluxes as well as steepening of spectral slope (S) of CDOM and accumulation of humic-like FDOM at the near-bottom waters over time indicated active microbial DOM utilization at the sediment–water interface, potentially stimulated by nitrate (NO3−) and nitrite (NO2−). The microbial DOC utilization rates, estimated in our study, may be sufficient to support denitrification rates of 0.2–1.4 mmol m−2 d−1, suggesting that sediment release of DOM contributes substantially to nitrogen loss processes in the ETSP off Peru.

Description: Oceanic Anoxic Events (OAEs) document major perturbations of the global carbon cycle with repercussions on the Earth’s climate and ocean circulation that are relevant to understand future climate trends. Here, we compare sedimentation patterns, nutrient cycling, organic carbon accumulation and carbon isotope variability across Cretaceous Oceanic Anoxic Events OAE1a and OAE2 in two drill cores with unusually high sedimentation rates from the Vocontian Basin (southern France) and Tarfaya Basin (southern Morocco). OAE1a and OAE2 exhibit remarkable similarities in the evolution of their δ13C excursion with long-lasting negative carbon isotope excursions preceding the onset of both anoxic events, supporting the view that OAEs were triggered by massive emissions of volcanic CO2 into the atmosphere. Based on analysis of cyclic sediment variations, we estimated the duration of the individual phases within the carbon isotope excursions. For both events, we identify: (1) a precursor phase lasting ~ 430 kyr and ~ 130 kyr, (2) an onset phase of ~ 390 and ~ 70 kyr, (3) a peak phase of ~ 600 and ~ 90 kyr, (4) a plateau phase of ~ 1400 and ~ 200 kyr and (5) a recovery phase of ~ 630 and ~ 440 kyr, respectively. The total duration of the positive carbon isotope excursion is estimated as 3400 kyr and 790 kyr and that of the main carbon accumulation phase as 980 kyr and 180 kyr, for OAE1a and OAE 2 respectively. The extended duration of the peak, plateau and recovery phases requires fundamental changes in global nutrient cycles either (1) through excess nutrient inputs to the oceans by increasing continental weathering and river discharge or (2) through nutrient-recycling from the marine sediment reservoir. We investigated the role of phosphorus on the development of carbon accumulation by analysing phosphorus speciation across OAE2 and the mid-Cenomanian Event (MCE) in the Tarfaya Basin. The ratios of organic carbon and total nitrogen to reactive phosphorus (Corg/Preact and Ntotal/Preact) prior to OAE2 and the MCE hover close to or below the Redfield ratio characteristic of marine organic matter. Decreases in reactive phosphorus resulting in Corg/Preact and Ntotal/Preact above the Redfield ratio during the later phase of OAE2 and the MCE indicate leakage from the sedimentary column into the water column under the influence of intensified and expanded oxygen minimum zones. These results suggest that a positive feedback loop, rooted in the benthic phosphorus cycle, contributed to increased marine productivity and carbon burial over an extended period of time during OAEs.

Description: Benthic fluxes of dissolved silica (Si) from sediments into the water column are driven by the dissolution of biogenic silica (bSiO2) and terrigenous Si minerals and modulated by the precipitation of authigenic Si phases. Each of these processes has a specific effect on the isotopic composition of silica dissolved in sediment pore waters such that the determination of pore water δ30Si values can help to decipher the complex Si cycle in surface sediments. In this study, the δ30Si signatures of pore fluids and bSiO2 in the Guaymas Basin (Gulf of California) were analyzed, which is characterized by high bSiO2 accumulation and hydrothermal activity. The δ30Si signatures were investigated in the deep basin, in the vicinity of a hydrothermal vent field, and at an anoxic site located within the pronounced oxygen minimum zone (OMZ). The pore fluid δ30Sipf signatures differ significantly depending on the ambient conditions. Within the basin, δ30Sipf is essentially uniform averaging +1.2 ± 0.1 ‰ (1SD). Pore fluid δ30Sipf values from within the OMZ are significantly lower (0.0 ± 0.5 ‰, 1SD), while pore fluids close to the hydrothermal vent field are higher (+2.0 ± 0.2 ‰, 1SD). Reactive transport modelling results show that the δ30Sipf is mainly controlled by silica dissolution (bSiO2 and terrigenous phases) and Si precipitation (authigenic aluminosilicates). Precipitation processes cause a shift to high pore fluid δ30Sipf signatures, most pronounced at the hydrothermal site. Within the OMZ however, additional dissolution of isotopically depleted Si minerals (e.g. clays) facilitated by high mass accumulation rates of terrigenous material (MARterr) is required to promote the low δ30Sipf signatures while precipitation of authigenic aluminosilicates seems to be hampered by high water / rock ratios. Guaymas OMZ δ30Sipf values are markedly different from those of the Peruvian OMZ, the only other marine setting where Si isotopes have been investigated to constrain early diagenetic processes. These differences highlight the fact that δ30Sipf signals in OMZs worldwide are not alike and each setting can result in a range of δ30Sipf values as a function of the environmental conditions. We conclude that the benthic silica cycle is more complex than previously thought and that additional Si isotope studies are needed to decipher the controls on Si turnover in marine sediment and the role of sediments in the marine silica cycle.