Articles | Volume 18, issue 15
Research article
13 Aug 2021
Research article |  | 13 Aug 2021

New insights into large-scale trends of apparent organic matter reactivity in marine sediments and patterns of benthic carbon transformation

Felipe S. Freitas, Philip A. Pika, Sabine Kasten, Bo B. Jørgensen, Jens Rassmann, Christophe Rabouille, Shaun Thomas, Henrik Sass, Richard D. Pancost, and Sandra Arndt

Constraining the mechanisms controlling organic matter (OM) reactivity and, thus, degradation, preservation, and burial in marine sediments across spatial and temporal scales is key to understanding carbon cycling in the past, present, and future. However, we still lack a detailed quantitative understanding of what controls OM reactivity in marine sediments and, consequently, a general framework that would allow model parametrization in data-poor areas. To fill this gap, we quantify apparent OM reactivity (i.e. OM degradation rate constants) by extracting reactive continuum model (RCM) parameters (a and v, which define the shape and scale of OM reactivity profiles, respectively) from observed benthic organic carbon and sulfate dynamics across 14 contrasting depositional settings distributed over five distinct benthic provinces. We further complement the newly derived parameter set with a compilation of 37 previously published RCM a and v estimates to explore large-scale trends in OM reactivity. Our analysis shows that the large-scale variability in apparent OM reactivity is largely driven by differences in parameter a (10−3–107) with a high frequency of values in the range 100–104 years. In contrast, and in broad agreement with previous findings, inversely determined v values fall within a narrow range (0.1–0.2). Results also show that the variability in parameter a and, thus, in apparent OM reactivity is a function of the whole depositional environment, rather than traditionally proposed, single environmental controls (e.g. water depth, sedimentation rate, OM fluxes). Thus, we caution against the simplifying use of a single environmental control for predicting apparent OM reactivity beyond a specific local environmental context (i.e. well-defined geographic scale). Additionally, model results indicate that, while OM fluxes exert a dominant control on depth-integrated OM degradation rates across most depositional environments, apparent OM reactivity becomes a dominant control in depositional environments that receive exceptionally reactive OM. Furthermore, model results show that apparent OM reactivity exerts a key control on the relative significance of OM degradation pathways, the redox zonation of the sediment, and rates of anaerobic oxidation of methane. In summary, our large-scale assessment (i) further supports the notion of apparent OM reactivity as a dynamic ecosystem property, (ii) consolidates the distributions of RCM parameters, and (iii) provides quantitative constraints on how OM reactivity governs benthic biogeochemical cycling and exchange. Therefore, it provides important global constraints on the most plausible range of RCM parameters a and v and largely alleviates the difficulty of determining OM reactivity in RCM by constraining it to only one variable, i.e. the parameter a. It thus represents an important advance for model parameterization in data-poor areas.

1 Introduction

Organic matter (OM) buried in marine sediments represents the largest reactive reservoir of reduced carbon on Earth (e.g. Hedges, 1992). The majority of OM that reaches modern global ocean seafloors originates from contemporary primary productivity (PP) in terrestrial (global net primary productivity, NPP = 59 Pg C yr−1) or marine environments (global NPP = 49 Pg C yr−1) (Hedges et al., 1997; Hedges and Oades, 1997; Kharbush et al., 2020). Additionally, it comprises OM derived from chemoautotrophy (Lengger et al., 2019; Vasquez-Cardenas et al., 2020) and ancient PP, the latter delivered through rock weathering, soil erosion, and thawing permafrost (Hedges et al., 1997; Kusch et al., 2021; Grotheer et al., 2020). Since OM production slightly exceeds degradation, a small fraction of the photosynthetically produced organic carbon escapes heterotrophic degradation and is buried in sediments (Berner, 2003; Burdige, 2006). This relatively small imbalance connects the short- and long-term carbon cycles, controls the long-term evolution of atmospheric CO2, and has enabled the accumulation of oxygen in the atmosphere (Berner, 2003; Hedges, 1992). Nevertheless, several data compilations reveal that rates of OM degradation and preservation vary significantly in space (Arndt et al., 2013; Burdige, 2006; Emerson and Hedges, 1988; Middelburg, 2019) and time (Berner, 2003). Significant progress has been made over the past decades in identifying the potential environmental controls of this spatio-temporal variability (Arndt et al., 2013; Bogus et al., 2012; LaRowe et al., 2020b; Zonneveld et al., 2010). It is becoming increasingly clear that the susceptibility of OM to heterotrophic degradation is not solely an inherent characteristic of the OM itself, but also depends on its environmental context (Kayler et al., 2019; LaRowe et al., 2020b; Mayer, 1995; Schmidt et al., 2011; Zonneveld et al., 2010). In fact, OM reactivity is determined by the dynamic interplay of a plethora of different factors (Arndt et al., 2013; Burdige, 2006; Middelburg, 2019; LaRowe et al., 2020b) (Table 1 and references therein). However, the relative significance of each of these factors remains poorly quantified. Consequently, we currently lack a general framework that allows quantification of the reactivity of OM in a given environmental context. This knowledge gap seriously compromises our ability to constrain OM reactivity in both space and time (Arndt et al., 2013; Lessin et al., 2018). The current lack of such a general framework ultimately limits our ability to quantify pathways and rates of biogeochemical processes over different temporal and spatial scales and, thus, predict the impact and feedbacks of past, present, and future climate change on global biogeochemical cycles and climate (Hülse et al., 2018).

Table 1Environmental factors influencing organic matter reactivity in marine sediments.

Download Print Version | Download XLSX

Traditionally, apparent (i.e. estimated representation of) OM reactivity has been determined through laboratory-based incubation experiments in shallow (few centimetres deep) and mixed sediments (e.g. Dai et al., 2009; Grossi et al., 2003; Westrich and Berner, 1984) or integrated model–data approaches in deeply buried sediments (hundreds of centimetres to metres deep and hundreds to thousands and millions of years old) (e.g. Boudreau and Ruddick, 1991; Middelburg, 1989; Wehrmann et al., 2013), using reaction–transport models (RTM). The latter are powerful tools since they allow the extraction of quantitative measures of OM reactivity from comprehensive multi-species porewater datasets integrated over greater sediment depths and ages (e.g. Middelburg, 2019). However, due to our limited quantitative understanding of the environmental controls on OM reactivity, the application of RTM approaches to data-poor settings is still complicated by the difficulty in constraining model parameters.

The need for predictive frameworks that constrain OM degradation model parameters has been recognized since the early days of RTM. Based on the rationale that the reactivity of OM settling onto the sediment is controlled by the degree of degradation, and thus its residence time in the oxygenated water column, early efforts have explored correlations between OM reactivity and water depth, sedimentation rate, or OM fluxes (Boudreau, 1997; Boudreau and Ruddick, 1991; Müller and Suess, 1979; Tromp et al., 1995). However, the derived relationships were based on limited datasets, and a more recent compilation of global data did not reveal statistically significant correlations between OM reactivity and single depositional environmental descriptors (Arndt et al., 2013). However, the latter combined OM degradation model parameters derived from structurally different models and used observational data of different complexity, thus somewhat compromising the comparability of their findings. Nevertheless, it emphasized the role of a complex interplay of environmental controls in reactivity (Kayler et al., 2019; LaRowe et al., 2020b; Mayer, 1995; Schmidt et al., 2011; Zonneveld et al., 2010).

The efforts by Arndt et al. (2013) to develop a general predictive framework highlighted two critical issues: OM reactivity is variable within and across settings with no significative global trends, and there is a crucial need to resolve the lack of comparative studies. Therefore, we here use a common, integrated model–data approach to inversely determine the model parameters that control apparent OM reactivity from contemporary sediment depth profiles (Sect. 3) at 14 different sites from five different depositional environments (Sect. 2). OM degradation is formulated according to the reactive continuum model (RCM) (Boudreau and Ruddick, 1991), which represents the most suitable degradation model type to explore links between reactivity parameters and environmental controls (see Sect. 3.2 for details). Derived RCM parameters are discussed in the environmental context of depositional settings. We further complement the newly derived parameter set with a compilation of previously published RCM parameter values to explore the distribution of model parameters and provide first-order constraints on model parametrization in data-poor areas. Finally, we evaluate the control of OM reactivity on benthic cycling and fluxes by linking RCM parameters to estimated OM degradation rates and relative contributions of metabolic pathways, anaerobic oxidation of methane (AOM) coupled to sulfate reduction, and fluxes of dissolved species across the sediment–water interface (SWI) at each of the 14 studied sites.

Figure 1Depositional environment geographic locations. (a) Global scale; (b) European margin; (c) Arabian Sea region. Map produced with Ocean Data View (Schlitzer, 2019).

2 Study sites

Our data–model study covers a wide range of depositional environments, from coastal to shelf and slope settings characterized by widely different environmental conditions (Fig. 1, Table 2). It comprises the following benthic provinces (Seiter et al., 2004): northern European margin (EUR1), NW Mediterranean margin off Rhone delta (EUR2), western Arabian Sea (WARAB), Bering Sea (NWPAC), and SW Atlantic off La Plata River (RIOPLATA). Most of those regions are characterized by predominantly temperate climate regimes, except for the Arabian Sea region (arid) and the Bering Sea area (snow/polar) (Chen and Chen, 2013). OM characteristics also cover a broad spectrum, ranging from organic-rich (TOC > 2 wt  %) to lean (TOC < 1 wt  %) sediments (Premuzic et al., 1982; Seiter et al., 2004) and from predominantly terrestrial (coast), to mixed terrestrial and marine (coast/shelf), to predominantly marine (shelf/slope) OM (Hedges et al., 1997; Meyers, 1997). The pelagic ecosystem structure reflects local and regional characteristics. Overall, as a common characteristic they display the occurrence of spring–summer phytoplankton blooms (e.g. Cowie, 2005; Rixen et al., 2019; Fleming-Lehtinen and Laamanen, 2012; Jensen et al., 1990; Coyle et al., 2008). This is also reflected by the composition of the phytoplankton community and thus OM production at the photic zone and vertical transport efficiency to sediments (e.g. Bach et al., 2019). This wide spectrum of depositional environments (Tables 2 and S1) offers a unique opportunity to quantify OM reactivity and explore its multidisciplinary links with environmental drivers.

Table 2Oceanographic context of depositional environments studied in the large-scale organic matter reactivity quantification.

NA: not available.

Download Print Version | Download XLSX

3 Methods

Here, we adopt an integrated model–data approach to determine inversely OM reactivity parameters from contemporary total organic carbon (TOC) and sulfate (SO42-) profiles measured at 14 sites across the five different depositional environments described in Sect. 2. The following sections provide an overview of the observational data (Sect. 3.1), a brief description of the applied RTM (Sect. 3.2.1; for further details see Supplements, Sect. S2), a detailed description of the inverse model approach used to extract RCM parameters from the observed depth profiles (Sect. 3.2.4), and a summary of the analysed model output (Sect. 3.2.5).

3.1 Observational data

We benefit from previously published datasets to develop our large-scale OM reactivity assessments. TOC and porewater SO42- depth profiles as well as site-specific physical characteristics and supporting data are available for all sites and were compiled from the PANGAEA database (, last access: 20 November 2020) (except from Rhone delta and Severn estuary datasets), as well as corresponding published literature (see below). Specific details about sampling strategies and analytical methods can be found in the respective publications. The dataset for the Severn estuary (Table 2; row A) is from Thomas (2014). Rhone delta data (Table 2; rows B and C) are from the DICASE project (Rassmann et al., 2016). Data for the northern European margin sites Aarhus Bay, Arkona Basin, and the Skagerrak transect (Table 2; rows D, E, G, H, and I) are part of the Europe-funded project METROL (Aquilina et al., 2010; Dale et al., 2008a, b; Knab et al., 2008; Mogollón et al., 2012). Data from the Helgoland mud area, North Sea (Table 2; row F) (Henkel and Kulkarni, 2015), the Arabian Sea oxygen minimum zone (OMZ) transect (Table 2; rows J, K, and L) (Bohrmann and cruise participants, 2010a, b, c), Bering Sea (Table 2; row M) (Gersonde, 2009), and Argentine Basin (Table 2; row N) (Henkel et al., 2011) are presented in the PANGAEA database (see above). We use those site-specific datasets to inform our data–model analysis and to constrain OM reactivity and related benthic processes.

3.2 Model description

3.2.1 General RTM framework

The Biogeochemical Reaction Network Simulator (BRNS) (Aguilera et al., 2005; Regnier et al., 2002) is an adaptive simulation environment that has been successfully employed to reproduce and quantify diagenetic processes in marine sediments across a wide range of depositional environments and timescales (Dale et al., 2008a; Thullner et al., 2009; Wehrmann et al., 2013). The BRNS is suitable for large, mixed kinetic–equilibrium reaction networks (Dale et al., 2009; Jourabchi, 2005; Thullner et al., 2009). The concentration depth profiles of solid and dissolved species in marine sediments are calculated according to the vertically resolved mass conservation equation of solid and dissolved species in porous media (Berner, 1980; Boudreau, 1997):

(1) σ C i t = z D bio σ C i z + D i σ C i z - σ ω C i z + α i σ C i 0 - C i + n s i n r n .

The first three terms on the right-hand side represent the transport process (bioturbation and molecular diffusion, advection, and bioirrigation), whereas the last term denotes the sum of all reactions (production and consumption) affecting species i. Table 2 summarizes all symbols employed here.

Briefly, the implemented reaction network encompasses the most pertinent primary and secondary redox reactions found in the upper layers of marine sediments (e.g. Aguilera et al., 2005; Thullner et al., 2009; Van Cappellen and Wang, 1996; Wang and Van Cappellen, 1996). It explicitly accounts for the heterotrophic degradation of OM coupled to the consumption of oxygen (aerobic OM degradation), nitrate (denitrification), sulfate (organoclastic sulfate reduction), and methanogenesis. Additionally, the reaction network accounts for secondary redox reactions, i.e. re-oxidation of reduced species produced during primary redox reactions. It explicitly resolves nitrification, sulfide re-oxidation by O2, anaerobic oxidation of methane (AOM) coupled to sulfate reduction, and CH4 reoxidation by O2.

The model transport and reaction equations were solved sequentially according to Regnier et al. (1998). The size of the model domain was fixed at 1000 cm for all sites, except for the Bering Sea, in which the model domain was extended to 1500 cm, due to the low sedimentation rate assumed for this site (Table 2). This choice is based on initial tests and ensures that the model domain covers the diagenetically most active zone, thus reducing the influence of biogeochemical dynamics in underlying sediments on biogeochemical dynamics within the model domain. The model equations were solved on an uneven grid and with a time-adapting time step (Δt< 0.01) run until steady state was reached. The Supplement includes a detailed description of the general RTM framework used here, its parametrization, and its numerical solution (Sect. S2).

3.2.2 OM degradation model

OM is composed of a complex and dynamic mixture of compounds that are distributed over a wide, continuous spectrum of reactivities. Thus, OM degradation is described by the reactive continuum model (RCM) (Boudreau and Ruddick, 1991), which assumes a continuous distribution of OM compounds over the entire reactivity spectrum. The RCM assumes a gamma distribution to describe the probability density function of OM distribution, om(k,t) (Aris, 1968; Boudreau and Ruddick, 1991; Ho and Aris, 1987). The overall rate of OM degradation RTOC is given by the integral

(2) r n = R TOC = - 0 k o m k , t d k .

Here, om(k,t) determines the concentration of OM having a degradability between k and k+dk at time t, with k being equivalent to the first-order degradation rate constant. The initial distribution (t=0) of OM is given by

(3) o m ( k , 0 ) = TOC ( 0 ) a v k v - 1 e - a k Γ ( v ) ,

where TOC(0) denotes the initial OM content, Γ(v) is the gamma distribution, a is the average lifetime of the more reactive components of bulk OM, and v represents the dimensionless scaling parameter of the distribution near k=0. The free, positive parameters a and v delineate the shape of the initial distribution of OM compounds along the range of k, and thus the overall reactivity of bulk OM. As such, the RCM approach requires the definition of two parameters that will define the shape of the OM distribution over reactivity k (Sect. 3.2.5, Eq. 5).

Due to the rapid depletion of the most reactive compounds, the reactivity of the bulk material decreases during degradation, reflecting the widely observed reactivity decrease with burial time, depth, and age (Boudreau and Ruddick, 1991; Middelburg, 1989). This indicates that degradation of OM proceeds at different rates in parallel. Interactions between different compounds or transformations of compounds can change the reactivity of a given compound. While the RCM does not explicitly account for such interactions and transformations, the overall OM profiles take these interactions implicitly into account. Thus, OM compounds are continuously and dynamically distributed over a range of reactivities that capture the decrease in apparent reactivity with burial age and depth as the most reactive compounds are successively degraded. The interplay of a and v drives the bulk OM reactivity depth profiles and consequently the yielded profiles from the reaction network implemented in the BRNS. The evolution of OM concentration as a function of depth z is given by

(4) TOC z = TOC ( 0 ) a a + age ( z ) v ,

where a and v are the free RCM parameters (Arndt et al., 2013; Boudreau and Ruddick, 1991) and age(z) denotes the age of the sediment layer at depth z. For non-bioturbated sediments (z> zbio) the burial age(z) can be calculated as a function of the burial velocity (Sect. S2.3, Eq. S8). However, within the bioturbated upper sediment layers, the age distribution of reactive species is controlled by both sedimentation, bioturbation, and the reactivity k of reactive species (Meile and Van Cappellen, 2005). In such cases, we apply a multi-G approximation for the RCM in the bioturbated sediments (Sect. S2.3, Eqs. S9–S18).

3.2.3 Boundary conditions

Boundary conditions place the BRNS in the environmental context of each of the study sites (Table 4). The boundary conditions are constrained based on either site-specific measurements or alternatively published data if direct observations were not available. Here, we assume a fixed boundary concentration of OM at the SWI. Concentrations of O2 are based on measurements in Skagerrak (Canfield et al., 1993) and the Rhone delta (Rassmann et al., 2016) sediments. Similarly, concentrations of NO3- are based on measurements in Skagerrak sediments (Canfield et al., 1993). Because of the lack of site-specific information, these O2 and NO3- measurements are adopted at all oxic sites (where SO42->28 mM) since they are representative for coastal and shallow slope sites, and their exact value exerts no influence on deep SO42- and CH4 concentrations. For anoxic sites (where SO42-< 28 mM), O2 and NO3- bottom water concentrations are set to zero. SO42- concentrations were constrained based on measurements for each simulated site. Concentrations of CH4 and HS were set to zero at SWI, as those species are rapidly oxidized in the overlying bottom water. A no-flux Cz=0 condition was implemented for all species at the lower boundary, assuming that biogeochemical dynamics in underlying sediments exert no influence on diagenetic processes in the model domain (e.g. Thullner et al., 2009).

3.2.4 Inverse modelling

We use an inverse model approach to extract the optimal OM reactivity parameter set a and v (i.e. the apparent OM reactivity) by assuming that the rank of a parameter set depends on the similarity between simulated and measured data. In general, all inverse modelling approaches are ultimately affected by a few limitations. Among others, the quality of the inverse modelling results is directly affected by the quality of the observational data. In particular, the lack of observations and data at the SWI due to core top loss during sediment sampling by gravity and piston corers can impact results. Additionally, processes that are not explicitly described in the applied model can impact inverse model results because reaction rate parameters implicitly account for those processes. Furthermore, downcore changes in concentration–depth profiles may also be influenced by transient features in bottom water conditions and/or fluxes, and the validity of the steady-state assumption thus needs to be carefully evaluated. Finally, multiple parameter sets might fit the observations equally well. Because both a and v exert an influence on the apparent OM reactivity and its evolution with depth (see below, Eq. 5), they are not completely independent parameters. For instance, a decrease in v (decreasing reactivity) can be compensated for by a simultaneous decrease in parameter a (increasing reactivity). Consequently, different a and v parameter couples might result in equally statistically satisfying fits between the simulated and observed profiles.

These limitations can be alleviated by using comprehensive, multi-component observational datasets, because the inclusion of additional information adds further constraints. However, such datasets are often not available. For instance, the first efforts to quantify a and v parameter values, and thus reactivity, on the basis of observations solely focused on downcore TOC profiles (Boudreau and Ruddick, 1991). Such an approach was instrumental for the development of the RCM and for producing the first large-scale dataset of a and v. However, the increasing availability and accessibility of multi-species datasets allow for the use of additional information that is available across different sites. Because OM heterotrophic degradation is the driver behind biogeochemical dynamics in marine sediments, downcore profiles of TEAs (e.g. SO42-) (e.g. Bowles et al., 2014; Jørgensen et al., 2019b) or reduced products of OM degradation (such as CH4 and NH4+) could be used as further constraints. However, data availability is often still limited; therefore, it is important to identify a minimum set of observational data that are widely available and comparably easy to measure. We tested the suitability of different artificial porewater datasets for the inverse determination of apparent OM reactivity. Using an artificially generated dataset (see Supplement Sect. S3), we ran the BRNS in the same fashion as we did for the site-specific data–model approach (see below). Our sensitivity analysis confirms that when TOC is considered as a single constraint, multiple pairs of a and v produce equivalent TOC depth profiles (Supplement Fig. S1). Model results show that including SO42- depth profiles as an additional constraint (Supplement Fig. S2) facilitates the identification of a best-fit a and v parameter set (Supplement Fig. S3). CH4 profiles were used as an additional qualitative constraint when the data were available, although measured ex situ CH4 profiles are often not reliable due to sampling and measurement uncertainties in deeper sediment layers where concentrations exceed millimolar levels (Dale et al., 2008a; Hensen et al., 2003; Hilligsøe et al., 2018). The dynamics of SO42- and CH4 are solely controlled by OM degradation and AOM (e.g. Bowles et al., 2014; Egger et al., 2018; Regnier et al., 2011), although their distributions can also be affected by bioirrigation in particularly shallow sulfate–methane transition zones (SMTZs) (e.g. Dale et al., 2019). In anoxic settings and deeply buried sediments, SO42- is the dominant TEA and CH4 is the most common reduced species (Jørgensen et al., 2019a, b). Thus, a combination of TOC, SO42-, and CH4 (if available to verify the depths of the SMTZ) depth profiles incorporates the information contained in the observed benthic sulfur and carbon dynamics and is sufficient to extract robust estimates of apparent OM reactivity and its evolution from the sediment–water interface down to the SMTZ. In addition, because changes in av exert different effects on TOC and SO42- depth profiles (see Figs. S1 and S2), including these two species reduces the impact of the av correlation on the uniqueness of fit. Thus, here we perform a site-specific data–model fit based on TOC and SO42- (and CH4) depth profiles.

The optimal parameter set was then determined for each site by assuming that the rank of a parameter set depends on the similarity between simulated and measured data. Best-fit parameter couples were found by minimizing the misfit between model results and simulations. We quantified the similarity between the simulated and observed TOC and SO42- depth profiles in terms of correlation coefficient (r), centred root-mean-square difference (RMSD), and standard deviations (SD). All these measures can be graphically represented in a single Taylor diagram (Taylor, 2001) that statistically quantifies the misfit between the observed and simulated data and visualizes how closely the simulated depth profile resembles the observed one.

The best-fit RCM parameters were inversely determined by first running the model for each site with a set of av couples over the entire range of previously published a (a=10-3–107 year) and v (v=10-2–100) values (Arndt et al., 2013, and references therein). The approach adopted for the sampling of the parameter space is a critical step in inverse modelling since the chosen statistical measures may reveal, in addition to a global optimum, multiple local optima. Because of the comparably low computational cost of the forward modelling, the low number of unknown parameters (here, a and v), and the weak nonlinearity of the problem, we used a simple two-step nested regular sampling of the two-dimensional parameter space. First, the entire plausible av space (see above) was sampled on a coarse regular grid (Δ log(a) = 1, Δv=0.02), and, for each of these runs, r, RMSD, and SD were calculated for both TOC and SO42- profiles. Based on these measures, a new, more finely resolved av grid (10×10) was defined around the av couples that showed the smallest combined misfit between simulated and observed TOC and SO42- profiles. On this grid, av couples were again sampled regularly, and the optimized av couple was determined by the lowest combined misfit (Table 5).

3.2.5 Quantification of organic matter reactivity and associated benthic processes

The apparent reactivity of bulk OM k for each site is calculated based on the RCM reactivity parameters a and v (Boudreau and Ruddick, 1991) that yield the best data–model fit and is given by

(5) k ( z ) = v a + age ( z ) .

The apparent reactivity of bulk OM at the sediment–water interface (z = 0) is thus given by

(6) k ( 0 ) = v a .

The BRNS simulates downcore concentration profiles for each species i that is explicitly resolved in the reaction network, as well as reaction rates rn (Table S4). From those reaction rates, the total rate of organic carbon oxidation, i.e. the depth integral of OM degradation for each primary redox reaction n, can be calculated by integrating the rate over the entire model domain L (e.g. Thullner et al., 2009):

(7) R TOC , total = n = 1 - 4 0 L r n d x = - 0 k o m k , t d k ,

where n=1, 2, 3, and 4 denote aerobic OM degradation, denitrification, organoclastic sulfate reduction, and methanogenesis, respectively (Table S4). The OM fluxes (Jin) at the SWI are calculated as the sum of OM out of the sediment (JOM,out) and the burial fluxes at the base of the model domain (JOM,bur) (e.g. Burdige, 2006). Here, we assume that JOM,out is equivalent to the depth-integrated rates of OM degradation (Eq. 7) and JOM,bur is solely governed by advective processes at depth z=1000 cm. As such


Additionally, the fluxes of species i at the SWI, i.e. the benthic–pelagic fluxes, are calculated from the model-derived depth profiles of species i according to e.g. Burdige (2006):

(11) J i = J i , Advection + J i , Diffusion + J i , Bioturbation + J i , Bioirrigation .

Positive Ji values indicate upward fluxes from the sediment, whereas negative values represent downward fluxes into the sediment.

4 Results and discussion

The integrated data–model analysis yielded a comprehensive picture of OM reactivity parameters (a, v, and k(z)), as well as OM degradation dynamics and benthic–pelagic coupling for our large-scale compilation of depositional environments (Table 6). Because model parameters implicitly account for all the processes that are not explicitly described in the model, this variability provides important insights into the environmental controls on OM reactivity.

Figure 2Total organic carbon data–model best fit: (a) Severn estuary; (b) Rhone pro-delta; (c) Rhone shelf; (d) Aarhus Bay; (e) Arkona Basin; (f) Helgoland mud area, North Sea; (g) Skagerrak – S10; (h) Skagerrak – S11; (i) Skagerrak – S13; (j) Arabian Sea – OMZ core; (k) Arabian Sea – OMZ transition; (l) Arabian Sea – below OMZ; (m) Bering Sea; (n) Argentine Basin.


4.1 Constraining apparent OM reactivity

4.1.1 Inverse modelling

We inversely determined the set of RCM parameters a and v that produces the best model fit to the observed TOC (Fig. 2) and SO42- (Fig. 3) depth profiles. Table 5 shows the quantitative assessments of these model fits and associated uncertainties. Overall, TOC profiles display a satisfying fit between model and observations. However, the quality of fit is often compromised by a lack of observations and data at the SWI due to core loss during sediment sampling by gravity and piston corers, poor downcore sampling resolution, and possible non-steady-state conditions. Generally, the latter may be problematic when trying to reproduce sedimentary conditions of dynamic shallow systems and long-term burial in shelf and deep-sea sites. At estuarine sites, steady-state conditions may be impacted by the highly dynamic nature of such settings. The Severn estuary exhibits strong tidal dynamics that result in continuous cycles of sediment resuspension (Manning et al., 2010). The Rhone delta experiences strong pulses of fresh water and sediments associated with flood events, as well as variable sedimentation rates (Antonelli et al., 2008; Cathalot et al., 2010; Zebracki et al., 2015). In shelf and slope settings, steady-state assumptions can be impacted by fluctuations on sedimentation rates and lateral transport of sediments across shelves (e.g. Henkel et al., 2011; Riedinger et al., 2014, 2017; Hensen et al., 2003; Cowie, 2005). Further, due to the possible loss of the upper few millimetres to centimetres during sampling with gravity and piston cores, it can be challenging to quantify OM contents at the SWI. Some of our sites (e.g. northern European sites; Aquilina et al., 2010; Knab et al., 2008) are likely to be impacted by such factors. Overall, such limitations in inverse model results are highlighted by the somewhat low correlation coefficients for TOC profiles in most cases (rTOC< 0.75; Table 3). We find that adding additional constraints to the model–data fit relieves such limitation. The SO42- depth profiles (Fig. 3) display a good agreement between present-day observations and model simulations, with the majority of sites exhibiting high correlation coefficients (rSO4>0.95; Table 3). Nevertheless, Fig. 3 shows cases where the model-derived sulfate depth profiles display some mismatches relative to observations (Arabian Sea below OMZ, Fig. 3l; Bering Sea, Fig. 3m). At the Bering Sea, the data–model discrepancy is likely associated with kinks in both TOC (Fig. 2m) and SO42- (Fig. 3m) profiles. They likely represent changes in the depositional regime over the studied timescale (e.g. Hensen et al., 2003) that are not captured by the model. The mismatch in the Arabian Sea at the site below the OMZ (Fig. 3l) could originate from SO42- recycling through H2S oxidation by MnO2 and Fe(OH)3 (e.g. Rassmann et al., 2020; Schulz et al., 1994); however we are unable to further explore this hypothesis due to a lack of metal oxide data to constrain those biogeochemical reactions.

Figure 3Sulfate data–model best fit: (a) Severn estuary; (b) Rhone pro-delta; (c) Rhone shelf; (d) Aarhus Bay; (e) Arkona Basin; (f) Helgoland mud area, North Sea; (g) Skagerrak – S10; (h) Skagerrak – S11; (i) Skagerrak – S13; (j) Arabian Sea – OMZ core; (k) Arabian Sea – OMZ transition; (l) Arabian Sea – below OMZ; (m) Bering Sea; (n) Argentine Basin.


Table 3Summary of model elements incorporated in the BRNS.

Download Print Version | Download XLSX

4.1.2 Apparent OM reactivity parameters v and a

The inversely determined OM reactivity parameters a and v (Sect. 3.2.4) derived from the best-fit model solutions (Table 6) reveal important differences in apparent OM reactivity and its distribution on a global scale. Alongside a and v parameters compiled from the literature (Arndt et al., 2009; Boudreau and Ruddick, 1991; Contreras et al., 2013; Henkel et al., 2012; Mewes et al., 2016; Mogollón et al., 2012, 2016; Wehrmann et al., 2013; Ye et al., 2016), they provide useful constraints on the plausible range of a and v, as well as their control on the global variability in apparent OM reactivity.

Parameter v exerts a scaling influence on OM reactivity distribution. High v values result in a high initial OM reactivity, whereas low v values produce low initial reactivity (Arndt et al., 2013). For instance, an extremely low v value of  0.01 for OM deposited during the late Cretaceous Oceanic Anoxic Event 2 at Demerara Rise (Arndt et al., 2009) has been linked to intense and rapid sulfurization of OM in the strongly euxinic subtropical Atlantic Ocean (Hülse et al., 2019). While previously published values for contemporary ocean sediments fall within the range of 0.1 to 1.08 (Arndt et al., 2013; Boudreau and Ruddick, 1991), inversely determined v values predominantly fall within the range of 0.1–0.2 (apparent 6th–11th order of reaction) – although lower and higher v values were also determined (Fig. 4a). Newly determined v parameters thus confirm the previously observed dominance of v values between 0.1 and 0.2 (Boudreau and Ruddick, 1991) and suggest that v is relatively constant across a wide range of different depositional environments encountered in the contemporary ocean. These findings also agree with the parametrization of another empirically derived continuum model description of OM degradation (Middelburg, 1989), which applies a globally constant parameter of 0.16:

(12) k z = 0.16 ( a + age ( z ) ) - 0.95 .

Although the parameters of the so-called power model cannot be directly compared to those of the RCM due to different exponents, the RCM is nevertheless equivalent to a general power law with an exponent of 1:

(13) k z = v ( a + age ( z ) ) - 1 .

Thus, the v parameter range of the RCM is in good agreement with the 0.16 value imposed in the power-law model.

Figure 4Frequency histogram of organic matter reactivity parameter distribution on a global scale: (a) scaling parameter – v; (b) shaping parameter – a.


Table 4Site-specific transport parameters and upper boundary conditions implemented in the reaction–transport model.

Download Print Version | Download XLSX

Nevertheless, our results also reveal exceptions. Particularly high v values (v> 0.25) were determined for the Arabian Sea and Bering Sea. Similarly, Boudreau and Ruddick (1991) also report exceptionally high values (v≈1) for the Peru Margin and north Philippine Sea but invoked the possibility of non-steady-state depositional conditions as a possible explanation. In our study, v>0.25 also coincides in part with the poorest data–model fits, such that the robustness of such high v values is uncertain (Table 5). In particular, the high v value in the Bering Sea (Table 6) might result from uncertainties (see above, Sect. 4.1.1) associated with boundary conditions and transient features (e.g. Henkel et al., 2011; Hensen et al., 2003). However, the model-derived v> 0.25 for the Arabian Sea agrees with previous diagenetic model results for that region. Luff et al. (2000) developed a 3G model to investigate OM dynamics and found that bulk OM contains exceptionally reactive compounds with first-order k of 0.2–30 yr−1. In the Arabian Sea, the OM flux to the sediment is supported by high primary production rates during the monsoon season (e.g. Cowie, 2005), complemented by deep carbon fixation by microorganisms (Lengger et al., 2019). Observations indicate an efficient vertical transport of diatom and haptophyte detritus to the sediment (Barlow et al., 1999; Cowie, 2005; Latasa and Bidigare, 1998; Shalapyonok et al., 2001). High transport rates and pronounced oxygen minimum zones might limit heterotrophic degradation of the sinking OM in the water column, thus delivering comparably unaltered OM to the sediment. Biomarker analysis of benthic OM indicates the deposition of highly reactive OM (Koho et al., 2013) and a reactivity decrease with extended exposure to oxic conditions in the water column (Vandewiele et al., 2009). Therefore, the somewhat atypical high v values determined (together with comparably low a values; see below) for the Arabian Sea OMZ transect are supported by interdisciplinary observations. Thus, v has a relatively narrow range of 0.1 to 0.2, with higher > 0.2 values generally restricted to depositional environments that are characterized by high primary production rates, anoxic water conditions, and enhanced vertical OM fluxes within the OMZ.

Parameter a determines both the overall apparent reactivity and the shape of the reactivity decrease with depth. High a values result in an overall low OM reactivity, but a slow decrease in reactivity with depth. In contrast, low a values result in a higher reactivity but also a more rapid decrease in reactivity with sediment depth (Arndt et al., 2013). Parameter a is often related to the degree of “pre-ageing” or degradation (e.g. Middelburg, 1989; Middelburg et al., 1993). Very fresh marine OM has been shown to be characterized by substantially low a values (a=3.1×10-4 years) (Boudreau and Ruddick, 1991), while the heavily sulfurized OM deposited during Cretaceous Oceanic Anoxic Event 2 in the strongly euxinic subtropical Atlantic Ocean is characterized by very high a values (a=104 years) (Arndt et al., 2009). Our inversely derived a values span several orders of magnitude (Table 6; Fig. 4b). They thus corroborate the high variability observed in previously published a-value compilations (e.g. Arndt et al., 2013; Boudreau and Ruddick, 1991) and suggest that the variability of parameter a exerts the major control on the spatial variability in OM reactivity. Additionally, inversely determined a values also broadly agree with previously observed global patterns. Low values (a< 102 years; i.e. high apparent OM reactivity) are generally determined for depositional environments that are characterized by the rapid deposition of predominantly fresh marine OM (e.g. Arabian Sea and shallow sites at the northern European margin), while depositional settings that receive pre-degraded OM and/or complex mixtures of OM (e.g. Skagerrak) exhibit high a (a> 103 years). However, despite the broad range of a values, the ensemble of inversely determined a values shows that most values fall into the range of 100104 years, thus narrowing the range of a values.

Table 5Summary of statistical tests for the best-fit av pair at each site derived from Taylor diagrams.

Download Print Version | Download XLSX

4.1.3 OM reactivity distributions

Inversely determined parameters a and v also allow the assessment of the apparent OM reactivity evolution during burial. Here, we calculate the probability distribution of OM over k at the time of deposition on the SWI, as well as after burial at 10, 100, and 1000 cm b.s.f. (centimetres below seafloor). Since v values are broadly similar across depositional settings, the differences in k distributions are mainly driven by the range of a values (Fig. 5). However, it is important to note that sedimentation rates ω and OM contents also influence the evolution of these distributions as they co-control degradation rates (Eqs. 5, S8). Therefore, the changes in OMreactivity distributions at the considered burial depths are not entirely driven by differences in reactivity, but more broadly reflect different exposure to heterotrophic degradation across sites.

Figure 5Apparent organic matter reactivity k initial probability distributions and downcore evolution as a function of burial depths: (a) Severn estuary; (b) Rhone pro-delta; (c) Rhone shelf; (d) Aarhus Bay; (e) Arkona Basin; (f) Helgoland mud area, North Sea; (g) Skagerrak – S10; (h) Skagerrak – S11; (i) Skagerrak – S13; (j) Arabian Sea – OMZ core; (k) Arabian Sea – OMZ transition; (l) Arabian Sea – below OMZ; (m) Bering Sea; (n) Argentine Basin.


The Arabian Sea region (Fig. 5j–l) is characterized by the deposition of highly reactive OM, the concentration of which rapidly decreases with sediment depth (high v and low a; Table 6). The OM reactivity distributions show a clear depletion of the most reactive OM compounds in the uppermost 10 cm. Below 100 cm b.s.f. the remaining bulk OM is dominated by less reactive compounds (k(z) < 10−5 yr−1). The rapid decrease in apparent OM reactivity reflects the often rapid deposition of predominantly marine OM (Cowie, 2005; Seiter et al., 2005), whose degradation during settling is limited by short residence times in combination with locally anoxic conditions (Luff et al., 2000). The shallow, sub-oxic and anoxic sediments of the northern European margin (Aarhus Bay – Fig. 5d; Arkona Basin – Fig. 5e) reveal a similar OM–k distribution and depth profile. However, the loss of the most reactive fractions within the uppermost 10 cm is less pronounced. Such a pattern may be attributed to the burial of terrestrial yet reactive OM (Aquilina et al., 2010). In contrast, Skagerrak sediments (S10 – Fig. 5g; S13 – Fig. 5i) are characterized by the most homogeneous, albeit comparably unreactive, OM mixture (high a values; Table 4). The OM–k distributions reveal small changes from the SWI down to 1000 cm b.s.f., reflecting a low overall reactivity, and thus low degradation rates. This somewhat surprising finding can be explained with the atypical character of the Skagerrak – a shallow environment that connects the brackish Baltic Sea with the saline North Sea and is characterized by an estuarine-like circulation. Due to this circulation pattern, the Skagerrak acts as the final sink for marine OM produced in the North Sea (Lohse et al., 1995; Van Weering et al., 1987). Consequently, more than 90 % of the OM deposited in the Skagerrak is derived from re-suspended and pre-degraded marine OM or terrestrial OM (ca. 20 % of the OM buried) (de Haas and van Weering, 1997) of generally low quality (Dauwe et al., 1999a, b).

Table 6Benthic organic matter dynamics derived from the coupled model–data assessments across the depositional settings.

n.d.: not determined (oxygen and nitrate benthic fluxes not determined for locations where SWI concentrations of those TEA were assumed to be zero).

Download Print Version | Download XLSX

Between those two endmember cases (high-reactivity Arabian Sea and low-reactivity Skagerrak), sediments from other sites generally reveal small changes in the continuous OM–k distributions within the uppermost 10 cm. At all sites, k values have converged to < 10−4 yr−1, below 100 cm b.s.f., where age(z) ≫ a, and changes in OM degradation rate constants with depth are thus largely controlled by age(z) (Eqs. 5, S8). As such, only the least reactive fractions are preserved, leading to an overall low apparent OM reactivity in deeper sediments that is largely independent of the initial apparent OM reactivity k(0). Sediments from the Bering Sea (Fig. 5m) are a notable exception. OM reactivity decreases gradually throughout the sediment column, reflecting both unusually high a and v values (Table 6). The Bering Sea is one of the most productive ocean regions, and OM deposition is mostly fuelled by spring and summer blooms (Coyle et al., 2008; Stabeno and Hunt, 2002; Stockwell et al., 2001). However, due to great water depths and low sedimentation rates, OM can experience extensive ageing and pre-processing prior to burial and result in the OM distribution pattern observed here (Fig. 5m).

4.1.4 Emerging environmental patterns in OM reactivity

Parameter a distributions exert the main control on the global variability in apparent OM reactivity. It also displays plausible values spanning several orders of magnitude (Fig. 4b, Table 6). Thus, our ability to quantify and predict benthic OM dynamics strongly depends on a general framework that allows constraint of parameter a when observational data are scarce or unavailable (e.g. global, future, and past model scenarios). Therefore, we further explore these links combining our newly generated dataset with those compiled from previously published RCM parameters (Fig. 6).

Figure 6Depositional environment controls on organic matter reactivity parameters a (a–d) and k (e–g): (a) v vs. a; (b, e) water depth h(c, f) sedimentation rate ω; (d, g) organic matter fluxes JOM,in.


Despite the large scatter, broad trends emerge that generally agree with previously proposed relationships kω (Boudreau, 1997; Tromp et al., 1995), kJOM,in (Boudreau, 1997; Müller and Suess, 1979), and aω (Boudreau and Ruddick, 1991) (Fig. 6). Depositional environments characterized by low sedimentation rates, low OM fluxes, or great water depths generally display high a values (i.e. low apparent reactivity), whereas settings characterized by high sedimentation rates, high OM fluxes, or shallow water depth generally display lower a values (high apparent reactivity). While these trends are extremely weak on a global scale, they are relatively more pronounced within a given geographical setting (i.e. on a local/regional scale).

However, in agreement with the previous global assessment (Arndt et al., 2013), we do not find any statistically significant global or regional correlation between the reactivity parameters (a and k(0)) and master characteristics of depositional environments. Earlier proposed empirical relationships were based on a limited number of observations available at that time, and our results, as well as previously compiled parameter values, thus confirm cautious statements regarding the global validity of previously published relationships (Boudreau, 1997; Müller and Suess, 1979; Tromp et al., 1995). These results further emphasize that the spatial variability in a and thus apparent OM reactivity are not determined by one single characteristic of the depositional environment. Instead, it results from a complex and dynamic interplay of multiple environmental controls (e.g. Kayler et al., 2019; LaRowe et al., 2020b; Mayer, 1995; Zonneveld et al., 2010; Schmidt et al., 2011), and the ensemble of results indicates that sedimentation rates and OM deposition fluxes are poor proxies for the complex interplay of environmental controls on OM reactivity (Fig. 6c, f).

However, Fig. 6b indicates that water depth seems to be a useful and easily accessible first-order proxy for this complex interplay of environmental controls on OM reactivity. Thus, it may serve as an acceptable predictor for first-order global reactivity patterns in the absence of more suitable information. At the global scale, inversely determined apparent OM reactivity for the deep, well-oxygenated sites (e.g. Bering Sea and Argentine Basin) fall on the lower end of the reactivity spectrum, whereas shallow coastal environments, such as the northern European region, reveal generally higher apparent OM reactivities. In addition, within the northern European region and along the Rhone delta transect, parameter a generally increases (i.e. decrease in apparent reactivity) with increasing water depth (Fig. 6b, e). These trends can be explained with the broad large-scale patterns of environmental conditions along the global hypsometry. In deeper depositional environments, prolonged pelagic degradation (Cram et al., 2018; Maerz et al., 2020; Müller and Suess, 1979; Weber et al., 2016), exposure to oxygen (Hartnett et al., 1998), and the continuous ageing during settling (Griffith et al., 2010) result in a low apparent OM reactivity in underlying sediments. In contrast, several processes occurring in shallow environments, such as rapid settling, intense lateral transport, limited degradation in the water column, dynamic reworking through intense macrofaunal activity, and fluctuating redox conditions (e.g. Aller, 1994; Canfield, 1994; Hedges et al., 1997; Burdige, 2007) contribute to higher apparent OM reactivities. Local reactivity–water depth trends in shallow environments are driven by some of the same controls.

The weak water depth–reactivity trends and numerous exceptions at both global and regional level caution against the uncritical use of simplistic kh relationships. Further important controls often complicate and sometimes inhibit the use of water depth as a predictor for apparent OM reactivity. For instance, the Severn estuary is an exception to this observed local, shallow environmental trend. This region is characterized by one of the largest tidal ranges on Earth, resulting in pronounced environmental gradients along the dynamic land–ocean transition. OM of terrestrial, anthropogenic, and marine sources mixes along the gradient and interacts with high loads of fine, suspended particular material. Additionally, intense tidally driven deposition–resuspension cycles (Dyer, 1984; Jonas and Millward, 2010; Manning et al., 2010; Uncles, 2010) result in a continuous re-exposure of benthic OM, which might cause the lower apparent reactivity of the buried OM. The Arabian Sea is another example where the kh relationship deviates from the general global trends. Here, reactivity is 2–3 orders of magnitude higher than that derived for other deep-sea sites. Instead, reactivity is rather comparable to values determined for shallow sites in the northern European margin and the Rhone delta. In addition, across the Arabian Sea OMZ transect, only a modest regional depth trend is observed, with a decrease in reactivity (i.e. increase in a) from the OMZ to increasing water depths and bottom water oxygen concentrations (Cowie, 2005; Luff et al., 2000). While distinct from the overall kh discussed above, the Arabian Sea OMZ region exhibits relatively similar patterns to other highly productive, anoxic water column systems, e.g. Peru margin (Boudreau and Ruddick, 1991); Blake ridge (Marquardt et al., 2010; Wallmann et al., 2006); and California, Chile, Costa Rica, and Namibia margins (Marquardt et al., 2010). This suggests a general separation from anoxic and oxic water columns shaping benthic OM reactivity.

4.2 Influence of OM reactivity on benthic–pelagic coupling

Both quantity and quality of the deposited OM exerts an important control on benthic–pelagic coupling processes with potential implications for global biogeochemical cycles and climate. While the influence of OM quantity on benthic–pelagic exchange along the global hypsometry has been previously investigated (Krumins et al., 2013; Soetaert et al., 1996, 1998; Thullner et al., 2009), the role of reactivity on benthic–pelagic fluxes remains less constrained.

4.2.1 Control of OM reactivity on OM degradation rates and respiration pathways

In agreement with previous findings (Henrichs, 1992), depth-integrated OM degradation rates (RTOC, total; Eq. 7) show a strong correlation with OM fluxes (JOM,in) (Fig. 7a, r2=0.97; p<0.01; n=14). In contrast, they do not reveal a strong link with apparent OM reactivity (i.e. parameter a), thus indicating that the magnitude of the OM deposition flux rather than its quality often exerts a first-order control on depth-integrated rates of OM degradation. In our dataset, apparent OM reactivity only becomes a dominant control on RTOC, total in depositional environments that receive exceptionally reactive OM (i.e. a< 101 years), such as the Arabian Sea region (Table 6).

Figure 7Benthic organic matter degradation dynamics. (a) Organic matter fluxes vs. depth-integrated organic matter degradation rates; (b) relative contribution of each metabolic pathway to heterotrophic organic matter degradation.


OM reactivity does exert an important impact on the relative contributions of metabolic pathways (Fig. 7). Overall, rates of OM degradation coupled to consumption of oxygen (aerobic OM degradation pathway) and nitrate (denitrification) make a small relative contribution (<10 %) to total degradation rates in the investigated oxic depositional environments (Rhone, Skagerrak, Bering Sea, and Argentine Basin; Fig. 7b). The relative contributions quantified here are generally lower than previous estimates for coastal and shelf sediments (>20 %) (e.g. Long Island Sound, Mackin and Swider, 1989; Skagerrak, Canfield et al., 1993; eastern Canadian continental margin, Boudreau et al., 1998; Barents Sea, Freitas et al., 2020; and northern Gulf of Mexico, Owings et al., 2021). Additionally, estimates are consistently lower than values in the deep sea (>70 %) (e.g. Cape Hatteras continental rise; Heggie et al., 1987), as well as over the global-scale hypsometry (Thullner et al., 2009). However, the overall dynamic nature of the considered environments, which favours the dominance of organoclastic sulfate reduction (Bowles et al., 2014; Jørgensen and Kasten, 2006; Thullner et al., 2009), as well as the length of sediment column (1000 cm b.s.f.) considered here renders a direct comparison difficult. Here, the contribution of the aerobic OM degradation pathway further decreases with a decrease in OM apparent reactivity (i.e. increase in a) as a consequence of the resulting higher burial of OM to sub-oxic or anoxic sediments (e.g. Meister et al., 2013). Additionally, the enhanced production of reduced species at depth further decreases the relative contribution of the aerobic OM degradation pathway to overall OM degradation rates by channelling more oxygen consumption through the re-oxidation of upward-diffusing reduced species (e.g. ammonium and hydrogen sulfide) (Glud, 2008; Hensen et al., 2006).

It has been suggested that anoxic metabolic pathways represent significant pathways of OM degradation (e.g. Bowles et al., 2014; Jørgensen and Kasten, 2006; Thullner et al., 2009). Our results show that sulfate reduction is the main oxidative pathway in regions characterized by high apparent OM reactivity (a< 102 years) and intermediate OM fluxes (JOM,in 10−4–10−3 mol C cm−2 yr−1), e.g. Arabian Sea region (Cowie, 2005; Luff et al., 2000) and Arkona Basin and Aarhus Bay (Dale et al., 2008b; Jørgensen et al., 2019a; Mogollón et al., 2012). Additionally, high relative contributions of sulfate reduction are also simulated for environments that are characterized by a reduced supply of less reactive (a> 103 years) OM deposition fluxes, such as the Argentine Basin and Bering Sea sites (e.g. Hensen et al., 2003). While high apparent OM reactivity enhances sulfate consumption in the upper sediment, where sulfate is efficiently replenished by exchange with bottom waters, a low OM reactivity allows for a deeper sulfate penetration (Fig. 3m–n) that sustains sulfate reduction over a greater sediment volume (e.g. Bowles et al., 2014; Meister et al., 2013; Riedinger et al., 2005, 2014, 2017). In contrast, methanogenesis is an important degradation pathway in environments that are characterized by a high supply (JOM,in>10-3–10−2 mol C cm−2 yr−1) of both highly reactive (a=100–101 years) OM (Rhone proximal zone and Helgoland mud area, North Sea) and the least reactive (a>103 years) OM (Skagerrak). The fast burial of reactive OM drives the development of a shallow methanogenic zone in the sediment (shallow sulfate penetration; Fig. 3b, f) (e.g. Aromokeye et al., 2020; Oni et al., 2015a, b; Rassmann et al., 2020), whereas very low apparent reactivity OM escapes sulfate reduction and is degraded deeper in the sediment via methanogenesis (e.g. Dale et al., 2008a; Knab et al., 2008; Meister et al., 2013). In both cases, the produced methane diffuses up and supports an enhanced consumption of sulfate through AOM (Sect. 4.2.2). Consequently, the relative contribution of organoclastic sulfate reduction decreases due to the presence of methane (e.g. Bowles et al., 2014; Niewöhner et al., 1998; Regnier et al., 2011; Jørgensen and Kasten, 2006; Jørgensen et al., 2019a; Riedinger et al., 2005, 2014, 2017).

We do not consider the degradation pathways mediated by metal oxides due to a lack of SWI data that would allow constraint of boundary conditions of those processes. This might result in a slight overrepresentation of sulfate reduction (e.g. at the Skagerrak; Canfield et al., 1993; Rysgaard et al., 2001). It has been demonstrated that metal oxide pathways represent a minor contribution to total OM heterotrophic degradation on a global scale (Thullner et al., 2009), although in continental margin sediments iron hydroxide reduction may represent an important OM degradation pathway (Beckler et al., 2016).

4.2.2 Anaerobic oxidation of methane dynamics

AOM occurs at the sulfate–methane transition zone (SMTZ), which can vary in sediment depth from centimetres to hundreds of metres as a function of different environmental controls, such as OM quantity and quality, sedimentation rate, active fluid flow, or microbial growth dynamics (Puglini et al., 2020; Regnier et al., 2011). Reflecting the diversity of the depositional environments studied here, we observe a broad range of SMTZ depths (Fig. 8a–c). Generally and in agreement with previous findings (Egger et al., 2018; Regnier et al., 2011), we observe a deepening of the SMTZ with decreasing sedimentation rates in addition to an important control of OM reactivity on the depth of the SMTZ (Fig. 8b) and the depth-integrated AOM rates (ΣAOM, Fig. 8d).

Figure 8Anaerobic oxidation of methane (AOM) patterns as a function of apparent organic matter reactivity and depositional rates. (a) AOM rates vs. SMTZ depth; SMTZ depths vs. (b) reactivity parameter a and (c) sedimentation rates ω; AOM rates vs. (d) reactivity parameter a and (e) sedimentation rates ω.


Figure 9Benthic fluxes of dissolved species (a–e) and transport mechanism relative contributions (f–j): (a, f) oxygen, (b, g) nitrate, (c, h) sulfate, (d, i) ammonium, and (e, j) phosphate.


The lowest ΣAOM and deepest SMTZ depths are determined for sites that are characterized by low apparent OM reactivity (i.e. deep-sea, well oxygenated bottom waters regions; Table 6). In those regions, low OM reactivity results in low rates of heterotrophic degradation (Fig. 7a) and, thus, deep SO42- penetration (Fig. 3m–n) or no SO42- depletion at all. In contrast, high OM reactivity in combination with high depositional fluxes (Fig. 8b–e) result in the complete depletion of TEAs in the upper sediment (z< 50 cm), allowing the establishment of methanogenesis at relatively shallow sediment depth. For instance, sediments from the Helgoland mud area, North Sea, are characterized by a shallow SMTZ and intense ΣAOM (Oni et al., 2015a, b). Yet, high OM reactivity does not always result in shallow SMTZ and intense AOM rates. Associated with moderate depositional fluxes, high apparent OM reactivity reduces the amount of reactive OM that reaches the methanogenic zone, inducing a deepening of the SMTZ and a decrease in ΣAOM. For instance, sediments from the Aarhus Bay (Flury et al., 2016) and Arkona Basin (Mogollón et al., 2012) display high OM reactivity (a< 101 years), which supports high rates of degradation close to the SWI (Fig. 7a) where SO42- is constantly replenished from bottom waters and limits OM fluxes to the deeper methanogenic sediment (e.g. Flury et al., 2016; Meister et al., 2013; Regnier et al., 2011). For moderate depositional fluxes, associated with lower OM reactivity (a> 102 years), e.g. the Skagerrak sites, pre-aged OM buried in these sediments (Aquilina et al., 2010) is predominantly degraded via methanogenesis (Fig. 7b). Consequently, a shoaling of the SMTZ is observed (Dale et al., 2008a; Knab et al., 2008), and higher ΣAOM values are determined (Table 6).

4.2.3 OM reactivity controls on sediment–water interface fluxes patterns

Model results suggest that apparent OM reactivity also influences nutrient recycling at the seafloor. Thus, we estimate benthic–pelagic fluxes of dissolved species (O2, NO3-, SO42-, NH4+, and PO43-). We emphasize that these flux calculations are based on simulated depth profiles at steady state and are thus only weakly affected by sampling limitations (e.g. possible loss of uppermost sediment layers associated with the choice of sampling devices; Hensen et al., 2006). Our estimates allow us to investigate the links between nutrient recycling and OM reactivity and to assess qualitative differences across sites.

Apparent OM reactivity broadly controls the spatial patterns of benthic–pelagic fluxes, as well as the relative importance of different transport processes (Glud, 2008; Bourgeois et al., 2017). Overall, molecular diffusion is the main transport pathway (Fig. 9f–j). However, the relative importance of bioirrigation increases at environments characterized by low deposition of OM and less reactive OM. This occurs because of the relative decrease in diffusive fluxes as a consequence of weak concentration gradients across the SWI. Bioturbation fluxes are generally low and reflect our assumption of depth-dependent Dbio (Middelburg et al., 1997). Advective fluxes become an important transport mechanism driving downward sulfate fluxes, especially in rapidly accumulating sediments and at sites characterized by low OM reactivity and moderate depositional rate.

Overall, high apparent OM reactivity drives high sedimentary uptake fluxes of dissolved TEA (Table 6; Fig. 9a–c). Sulfate largely dominates the benthic TEA uptake, particularly at sites where oxygen and nitrate are unavailable (Fig. 9c). At the Rhone proximal zone, unusually high sedimentation rates (Pastor et al., 2011) deliver fresh and highly reactive OM (Cathalot et al., 2010; Pruski et al., 2015), supporting high OM fluxes and degradation rates (Fig. 7a; Table 6). Consequently, intense oxygen consumption occurs at the SWI (Rassmann et al., 2016) driven by both aerobic degradation of OM and re-oxidation of reduced species. Similarly, in the Arabian Sea sediments below the OMZ the high input of fresh phytoplankton debris (Cowie, 2005; Koho et al., 2013; Rixen et al., 2019; Vandewiele et al., 2009) and potentially chemoautotroph biomass (Lengger et al., 2019) associated with intense degradation rates in the uppermost sediment layers (Luff et al., 2000) result in a high relative contribution of aerobic degradation of OM and denitrification to total OM oxidation. Consequently, both benthic oxygen and nitrate uptake fluxes are comparably high. In contrast, at regions characterized by deposition of less reactive OM (k< 10−4 yr−1; Skagerrak, Argentine Basin, and Bering Sea), benthic–pelagic fluxes of TEAs are generally low.

Benthic fluxes play a crucial role in nutrient recycling. In the Arabian Sea region, large phytoplankton blooms are associated with monsoon conditions, which result in upwelling of nutrient-rich bottom water (SW monsoon) and deepening of the mixed layer (NE monsoon) (Cowie, 2005; Luff et al., 2000; Rixen et al., 2019). Similarly, across the northern European margin, spring and summer diatom blooms are common (Fleming-Lehtinen and Laamanen, 2012; Jensen et al., 1990; Karlson et al., 1996; Lomstein et al., 1990; Wiltshire and Manly, 2004) and maintained by benthic nutrient fluxes. Additionally, bottom water upwelling and nutrient recycling are important mechanisms sustaining spring and summer PP at the Bering Sea (Coyle et al., 2008). Benthic ammonium and phosphate recycling fluxes mirror the reactivity and OM degradation patterns (Fig. 9d–e). The largest fluxes are determined for the Arabian Sea region, as well as the shallow Aarhus Bay and Arkona Basin (Table 6). In contrast, the lowest nutrient recycling fluxes are observed in regions characterized by the deposition of less reactive OM (phosphate fluxes are negligible at Skagerrak sites S10 and S13), where heterotrophic degradation rates are slow.

5 Conclusions

We developed and applied an inverse modelling approach to quantify apparent OM reactivity (i.e. parameters a and v) based on a common minimum set of benthic observational data, as well as a common RTM approach for 14 different sites across five different depositional environments. We provide important first-order constraints on the reactive continuum model (RCM) parameterization that can inform on the choice of parameter values in data-poor areas, including global-scale and timescale scenarios.

Our findings corroborate previous results (Arndt et al., 2013; Boudreau and Ruddick, 1991; Forney and Rothman, 2012; Middelburg, 1989) that the RCM parameter v is globally relatively constant (v=0.1–0.2). Exceptionally high v> 0.2 is often associated with deposition and burial of highly reactive phytoplankton debris in high-productivity regions associated with well-established OMZs. In contrast, in agreement with previous findings, RCM parameter a can span several orders of magnitude at a global scale, suggesting that the parameter a is the main parameter describing the variability of apparent OM reactivity. Consequently, future modelling efforts to quantify OM reactivity on a global scale could be reduced to one main reactivity variable. Based on our results and previous findings, we consolidate the range of predominant a distribution to 100–104 years. This is a valuable constraint when dealing with data-poor regions and timescales (e.g. Hülse et al., 2018), since it excludes extreme values at both ends of the a range.

Additionally, results indicate that the large variability in apparent OM reactivity is linked to a combination of multiple environmental drivers that control the quality of OM delivered to sediments and the timescale of settling and burial. Therefore, our findings contribute to the notion that apparent OM reactivity is controlled by a complex and dynamic interplay of environmental drivers, which are measurable and allow the quantification of a, v, and k. Our model results caution against the use of a single environmental master variable such as water depth, sedimentation rate, or OM deposition fluxes beyond well-defined geographical scales. Yet, results show that, if no other information is available, water depth can serve as a useful proxy for the complex and dynamic interplay of environmental drivers and can be used to predict first-order, large-scale OM reactivity patterns. However, more holistic, interdisciplinary exploration of OM reactivity in its entire environmental context is needed to mechanistically understand benthic carbon and nutrient dynamics.

Finally, results also show that OM reactivity exerts a dominant control on the redox zonation of the sediment, the depth of the SMTZ, depth-integrated AOM rates, and benthic–pelagic exchange fluxes, indicating that these processes could serve as predictor variables for OM reactivity. In contrast, depth-integrated OM degradation rates are largely controlled by the magnitude rather than the quality of OM deposition to the sediment.

Code and data availability

The model code is available at the GitHub repository (, Freitas, 2021). The data are available in the cited literature in the paper.


The supplement related to this article is available online at:

Author contributions

FSF, SA, and RDP designed this study. PAP developed and implemented changes to model description (multi-G approximation of the RCM). SK, BBJ, JR, CR, ST, and HS provided porewater and sediment data to inform model experiments and contributed to discussing and interpreting the data and model results. FSF compiled data, ran all model simulations, and performed data–model integration. FSF wrote the manuscript with significant contributions from SA. All co-authors edited and approved the manuscript.

Competing interests

The authors declare that they have no conflict of interest.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


Felipe S. Freitas was supported by a PhD scholarship from the Science without Borders programme (Ciência sem Fronteiras), sponsored by the CAPES Foundation within the Ministry of Education, Brazil, grant number 9999.009541/2013-06. Felipe S. Freitas also thanks the UKRI Natural Environment Research Council (NERC) for funding, grant number NE/P006493/1 (Changing Arctic Ocean Seafloor – ChAOS). Richard D. Pancost acknowledges the advanced ERC grant “The greenhouse earth system” (T-GRES, project reference 340923) and the Royal Society Wolfson Research Merit Award. Sandra Arndt and Philip A. Pika acknowledge funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement number 643052 745 (C-CASCADES). Shaun Thomas was funded by a President's Scholarship, Cardiff University. Christophe Rabouille and Jens Rassmann were supported by the INSU/EC2CO programme MissRhoDia and by the French state programme “Investissement d'avenir” run by the National Research Agency (AMORAD project ANR-11-RSNR-0002). The authors also thank Bernard Boudreau and the anonymous reviewer as well as associate editor Marilaure Grégoire for their constructive comments on early versions of the paper. Finally, we thank co-editor in chief Steve Bouillon for granting us a discount on the APC.

Financial support

This research has been supported by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (grant no. 9999.009541/2013-06).

Review statement

This paper was edited by Marilaure Grégoire and reviewed by Bernard Boudreau and one anonymous referee.


Aguilera, D. R., Jourabchi, P., Spiteri, C., and Regnier, P.: A knowledge-based reactive transport approach for the simulation of biogeochemical dynamics in Earth systems, Geochem. Geophy. Geosy., 6, 1–18,, 2005. 

Aller, R. C.: Diagenetic Processes Near the Sediment-Water Interface of Long Island Sound, I., in: Advances in Geophysics, Vol. 22, Elsevier, 237–350,, 1980. 

Aller, R. C.: Bioturbation and remineralization of sedimentary organic matter: effects of redox oscillation, Chem. Geol., 114, 331–345,, 1994. 

Aller, R. C. and Aller, J. Y.: The effect of biogenic irrigation intensity and solute exchange on diagenetic reaction rates in marine sediments, J. Mar. Res., 56, 905–936,, 1998. 

Aller, R. C. and Cochran, J. K.: The Critical Role of Bioturbation for Particle Dynamics, Priming Potential, and Organic C Remineralization in Marine Sediments: Local and Basin Scales, Front. Earth Sci., 7, 157,, 2019. 

Antonelli, C., Eyrolle, F., Rolland, B., Provansal, M., and Sabatier, F.: Suspended sediment and 137Cs fluxes during the exceptional December 2003 flood in the Rhone River, southeast France, Geomorphology, 95, 350–360,, 2008. 

Aquilina, A., Knab, N. J., Knittel, K., Kaur, G., Geissler, A., Kelly, S. P., Fossing, H., Boot, C. S., Parkes, R. J., Mills, R. A., Boetius, A., Lloyd, J. R., and Pancost, R. D.: Biomarker indicators for anaerobic oxidizers of methane in brackish-marine sediments with diffusive methane fluxes, Org. Geochem., 41, 414–426,, 2010. 

Aris, R.: Prolegomena to the rational analysis of systems of chemical reactions II. Some addenda, Arch. Ration. Mech. Anal., 27, 356–364,, 1968. 

Arndt, S., Hetzel, A., and Brumsack, H.-J.: Evolution of organic matter degradation in Cretaceous black shales inferred from authigenic barite: A reaction-transport model, Geochim. Cosmochim. Ac., 73, 2000–2022,, 2009. 

Arndt, S., Jørgensen, B. B., LaRowe, D. E., Middelburg, J. J., Pancost, R. D., and Regnier, P.: Quantifying the degradation of organic matter in marine sediments: A review and synthesis, Earth-Sci. Rev., 123, 53–86,, 2013. 

Aromokeye, D. A., Oni, O. E., Tebben, J., Yin, X., Richter-Heitmann, T., Wendt, J., Nimzyk, R., Littmann, S., Tienken, D., Kulkarni, A. C., Henkel, S., Hinrichs, K.-U., Elvert, M., Harder, T., Kasten, S., and Friedrich, M. W.: Crystalline iron oxides stimulate methanogenic benzoate degradation in marine sediment-derived enrichment cultures, ISME J., 15, 965–980,, 2020. 

Ausín, B., Bruni, E., Haghipour, N., Welte, C., Bernasconi, S. M., and Eglinton, T. I.: Controls on the abundance, provenance and age of organic carbon buried in continental margin sediments, Earth Planet. Sc. Lett., 558, 116759,, 2021. 

Bach, L. T., Stange, P., Taucher, J., Achterberg, E. P., Algueró-Muñiz, M., Horn, H., Esposito, M., and Riebesell, U.: The Influence of Plankton Community Structure on Sinking Velocity and Remineralization Rate of Marine Aggregates, Global Biogeochem. Cy., 33, 971–994,, 2019. 

Barlow, R. G., Mantoura, R. F. C., and Cummings, D. G.: Monsoonal influence on the distribution of phytoplankton pigments in the Arabian Sea, Deep-Sea Res. Pt. II, 46, 677–699,, 1999. 

Beckler, J. S., Kiriazis, N., Rabouille, C., Stewart, F. J., and Taillefert, M.: Importance of microbial iron reduction in deep sediments of river-dominated continental-margins, Mar. Chem., 178, 22–34,, 2016. 

Berner, R. A.: Early diagenesis: a theoretical approach, Princeton University Press, Princeton, NJ, 241 pp., 1980. 

Berner, R. A.: The long-term carbon cycle, fossil fuels and atmospheric composition, Nature, 426, 323–326,, 2003. 

Bianchi, T. S.: The role of terrestrially derived organic carbon in the coastal ocean: A changing paradigm and the priming effect, P. Natl. Acad. Sci. USA, 108, 19473,, 2011. 

Bianchi, T. S., Thornton, D. C. O., Yvon-Lewis, S. A., King, G. M., Eglinton, T. I., Shields, M. R., Ward, N. D., and Curtis, J.: Positive priming of terrestrially derived dissolved organic matter in a freshwater microcosm system, Geophys. Res. Lett., 42, 5460–5467,, 2015. 

Bianchi, T. S., Cui, X., Blair, N. E., Burdige, D. J., Eglinton, T. I., and Galy, V.: Centers of organic carbon burial and oxidation at the land-ocean interface, Org. Geochem., 115, 138–155,, 2018. 

Bogus, K. A., Zonneveld, K. A. F., Fischer, D., Kasten, S., Bohrmann, G., and Versteegh, G. J. M.: The effect of meter-scale lateral oxygen gradients at the sediment-water interface on selected organic matter based alteration, productivity and temperature proxies, Biogeosciences, 9, 1553–1570,, 2012. 

Bohrmann, G. and cruise participants: Documentation of sediment core GeoB12308-5, PANGAEA,, 2010a. 

Bohrmann, G. and cruise participants: Documentation of sediment core GeoB12309-3, PANGAEA,, 2010b. 

Bohrmann, G. and cruise participants: Documentation of sediment core GeoB12312-3, PANGAEA,, 2010c. 

Boudreau, B. P.: Is burial velocity a master parameter for bioturbation?, Geochim. Cosmochim. Ac., 58, 1243–1249,, 1994. 

Boudreau, B. P.: Diagenetic Models and Their Implementation: Modelling Transport and Reactions in Aquatic Sediments, Springer Berlin Heidelberg, Berlin, Heidelberg, 1997. 

Boudreau, B. P. and Ruddick, B. R.: On a reactive continuum representation of organic matter diagenesis, Am. J. Sci., 291, 507–538,, 1991. 

Boudreau, B. P., Mucci, A., Sundby, B., Luther, G. W., and Silverberg, N.: Comparative diagenesis at three sites on the Canadian continental margin, J. Mar. Res., 56, 1259–1284,, 1998. 

Bourgeois, S., Archambault, P., and Witte, U.: Organic matter remineralization in marine sediments: A Pan-Arctic synthesis: Pan-Arctic Benthic Remineralization, Global Biogeochem. Cy., 31, 190–213,, 2017. 

Bowles, M. W., Mogollon, J. M., Kasten, S., Zabel, M., and Hinrichs, K.-U.: Global rates of marine sulfate reduction and implications for sub-sea-floor metabolic activities, Science, 344, 889–891,, 2014. 

Bradley, J. A., Arndt, S., Amend, J. P., Burwicz, E., Dale, A. W., Egger, M., and LaRowe, D. E.: Widespread energy limitation to life in global subseafloor sediments, Sci. Adv., 6, eaba0697,, 2020. 

Burdige, D. J.: Geochemistry of marine sediments, Princeton University Press, Princeton, NJ, 609 pp., 2006. 

Burdige, D. J.: Preservation of Organic Matter in Marine Sediments: Controls, Mechanisms, and an Imbalance in Sediment Organic Carbon Budgets?, Chem. Rev., 107, 467–485,, 2007. 

Canfield, D. E.: Factors influencing organic carbon preservation in marine sediments, Chem. Geol., 114, 315–329,, 1994. 

Canfield, D. E., Jørgensen, B. B., Fossing, H., Glud, R., Gundersen, J., Ramsing, N. B., Thamdrup, B., Hansen, J. W., Nielsen, L. P., and Hall, P. O. J.: Pathways of organic carbon oxidation in three continental margin sediments, Mar. Geol., 113, 27–40,, 1993. 

Cathalot, C., Rabouille, C., Pastor, L., Deflandre, B., Viollier, E., Buscail, R., Grémare, A., Treignier, C., and Pruski, A.: Temporal variability of carbon recycling in coastal sediments influenced by rivers: assessing the impact of flood inputs in the Rhône River prodelta, Biogeosciences, 7, 1187–1205,, 2010. 

Chen, D. and Chen, H. W.: Using the Köppen classification to quantify climate variation and change: An example for 1901–2010, Environ. Dev., 6, 69–79,, 2013. 

Contreras, S., Meister, P., Liu, B., Prieto-Mollar, X., Hinrichs, K.-U., Khalili, A., Ferdelman, T. G., Kuypers, M. M. M., and Jorgensen, B. B.: Cyclic 100-ka (glacial-interglacial) migration of subseafloor redox zonation on the Peruvian shelf, P. Natl. Acad. Sci. USA, 110, 18098–18103,, 2013. 

Cowie, G.: The biogeochemistry of Arabian Sea surficial sediments: A review of recent studies, Prog. Oceanogr., 65, 260–289,, 2005. 

Coyle, K. O., Pinchuk, A. I., Eisner, L. B., and Napp, J. M.: Zooplankton species composition, abundance and biomass on the eastern Bering Sea shelf during summer: The potential role of water-column stability and nutrients in structuring the zooplankton community, Deep-Sea Res. Pt. II, 55, 1775–1791,, 2008. 

Cram, J. A., Weber, T., Leung, S. W., McDonnell, A. M. P., Liang, J.-H., and Deutsch, C.: The Role of Particle Size, Ballast, Temperature, and Oxygen in the Sinking Flux to the Deep Sea, Global Biogeochem. Cy., 32, 858–876,, 2018. 

Dai, J., Sun, M.-Y., Culp, R. A., and Noakes, J. E.: A laboratory study on biochemical degradation and microbial utilization of organic matter comprising a marine diatom, land grass, and salt marsh plant in estuarine ecosystems, Aquat. Ecol., 43, 825–841,, 2009. 

Dale, A. W., Regnier, P., Knab, N. J., Jørgensen, B. B., and Van Cappellen, P.: Anaerobic oxidation of methane (AOM) in marine sediments from the Skagerrak (Denmark): II. Reaction-transport modeling, Geochim. Cosmochim. Ac., 72, 2880–2894,, 2008a. 

Dale, A. W., Aguilera, D. R., Regnier, P., Fossing, H., Knab, N. J., and Jørgensen, B. B.: Seasonal dynamics of the depth and rate of anaerobic oxidation of methane in Aarhus Bay (Denmark) sediments, J. Mar. Res., 66, 127–155,, 2008b. 

Dale, A. W., Brüchert, V., Alperin, M., and Regnier, P.: An integrated sulfur isotope model for Namibian shelf sediments, Geochim. Cosmochim. Ac., 73, 1924–1944,, 2009. 

Dale, A. W., Flury, S., Fossing, H., Regnier, P., Røy, H., Scholze, C., and Jørgensen, B. B.: Kinetics of organic carbon mineralization and methane formation in marine sediments (Aarhus Bay, Denmark), Geochim. Cosmochim. Ac., 252, 159–178,, 2019. 

Dauwe, B., Middelburg, J. J., Van Rijswijk, P., Sinke, J., Herman, P. M. J., and Heip, C. H. R.: Enzymatically hydrolyzable amino acids in North Sea sediments and their possible implication for sediment nutritional values, J. Mar. Res., 57, 109–134,, 1999a. 

Dauwe, B., Middelburg, J. J., Herman, P. M. J., and Heip, C. H. R.: Linking diagenetic alteration of amino acids and bulk organic matter reactivity, Limnol. Oceanogr., 44, 1809–1814,, 1999b. 

de Haas, H. and van Weering, T. C. E.: Recent sediment accumulation, organic carbon burial and transport in the northeastern North Sea, Mar. Geol., 136, 173–187,, 1997. 

Dyer, K. R.: Sedimentation processes in the Bristol Channel/Severn Estuary, Mar. Pollut. Bull., 15, 53–57,, 1984. 

Egger, M., Riedinger, N., Mogollón, J. M., and Jørgensen, B. B.: Global diffusive fluxes of methane in marine sediments, Nat. Geosci., 11, 421–425,, 2018. 

Emerson, S. and Hedges, J. I.: Processes controlling the organic carbon content of open ocean sediments, Paleoceanography, 3, 621–634,, 1988. 

Faust, J. C., Tessin, A., Fisher, B. J., Zindorf, M., Papadaki, S., Hendry, K. R., Doyle, K. A., and März, C.: Millennial scale persistence of organic carbon bound to iron in Arctic marine sediments, Nat. Commun., 12, 275,, 2021. 

Fleming-Lehtinen, V. and Laamanen, M.: Long-term changes in Secchi depth and the role of phytoplankton in explaining light attenuation in the Baltic Sea, Estuar. Coast. Shelf Sci., 102/103, 1–10,, 2012. 

Flury, S., Røy, H., Dale, A. W., Fossing, H., Tóth, Z., Spiess, V., Jensen, J. B., and Jørgensen, B. B.: Controls on subsurface methane fluxes and shallow gas formation in Baltic Sea sediment (Aarhus Bay, Denmark), Geochim. Cosmochim. Ac., 188, 297–309,, 2016. 

Forney, D. C. and Rothman, D. H.: Inverse method for estimating respiration rates from decay time series, Biogeosciences, 9, 3601–3612,, 2012. 

Freitas, F. S., Hendry, K. R., Henley, S. F., Faust, J. C., Tessin, A. C., Stevenson, M. A., Abbott, G. D., März, C., and Arndt, S.: Benthic-pelagic coupling in the Barents Sea: an integrated data-model framework, Phil. Trans. R. Soc. A., 378, 1–15,, 2020. 

Freitas, F. S.: Model code for New insights into large-scale trends of apparent organic matter reactivity in marine sediments and patterns of benthic carbon transformation, Zenodo [data set], (last access: 9 August 2021), 2021. 

Gersonde, R.: The expedition of the research vessel “Sonne” to the subpolar North Pacific and the Bering Sea in 2009 (SO202-INOPEX) (Reports on polar and marine research No. 643), Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, 2009. 

Glud, R. N.: Oxygen dynamics of marine sediments, Mar. Biol. Res., 4, 243–289,, 2008. 

Griffith, D. R., Martin, W. R., and Eglinton, T. I.: The radiocarbon age of organic carbon in marine surface sediments, Geochim. Cosmochim. Ac., 74, 6788–6800,, 2010. 

Grossi, V., Caradec, S., and Gilbert, F.: Burial and reactivity of sedimentary microalgal lipids in bioturbated Mediterranean coastal sediments, Mar. Chem., 81, 57–69,, 2003. 

Grotheer, H., Meyer, V., Riedel, T., Pfalz, G., Mathieu, L., Hefter, J., Gentz, T., Lantuit, H., Mollenhauer, G., and Fritz, M.: Burial and Origin of Permafrost-Derived Carbon in the Nearshore Zone of the Southern Canadian Beaufort Sea, Geophys. Res. Lett., 47, e2019GL085897,, 2020. 

Guenet, B., Danger, M., Abbadie, L., and Lacroix, G.: Priming effect: bridging the gap between terrestrial and aquatic ecology, Ecology, 91, 2850–2861,, 2010. 

Hartnett, H. E., Keil, R. G., Hedges, J. I., and Devol, A. H.: Influence of oxygen exposure time on organic carbon preservation in continental margin sediments, Nature, 391, 572–575,, 1998. 

Hedges, J. I.: Global biogeochemical cycles: progress and problems, Mar. Chem., 39, 67–93,, 1992. 

Hedges, J. I. and Oades, J. M.: Comparative organic geochemistries of soils and marine sediments, Org. Geochem., 27, 319–361,, 1997. 

Hedges, J. I., Keil, R. G., and Benner, R.: What happens to terrestrial organic matter in the ocean?, Org. Geochem., 27, 195–212,, 1997. 

Heggie, D., Maris, C., Hudson, A., Dymond, J., Beach, R., and Cullen, J.: Organic carbon oxidation and preservation in NW Atlantic continental margin sediments, Geol. Soc. Lond. Spec. Publ., 31, 215,, 1987. 

Hemingway, J. D., Rothman, D. H., Grant, K. E., Rosengard, S. Z., Eglinton, T. I., Derry, L. A., and Galy, V. V.: Mineral protection regulates long-term global preservation of natural organic carbon, Nature, 570, 228–231,, 2019. 

Henkel, S. and Kulkarni, A.: Documentation of sediment core HE443/10-3,, 2015. 

Henkel, S., Strasser, M., Schwenk, T., Hanebuth, T. J. J., Hüsener, J., Arnold, G. L., Winkelmann, D., Formolo, M., Tomasini, J., Krastel, S., and Kasten, S.: An interdisciplinary investigation of a recent submarine mass transport deposit at the continental margin off Uruguay, Geochem. Geophy. Geosy., 12, 1–19,, 2011. 

Henkel, S., Mogollón, J. M., Nöthen, K., Franke, C., Bogus, K., Robin, E., Bahr, A., Blumenberg, M., Pape, T., Seifert, R., März, C., de Lange, G. J., and Kasten, S.: Diagenetic barium cycling in Black Sea sediments – A case study for anoxic marine environments, Geochim. Cosmochim. Ac., 88, 88–105,, 2012. 

Henrichs, S. M.: Early diagenesis of organic matter in marine sediments: progress and perplexity, Mar. Chem., 39, 119–149,, 1992. 

Hensen, C., Zabel, M., Pfeifer, K., Schwenk, T., Kasten, S., Riedinger, N., Schulz, H. D., and Boetius, A.: Control of sulfate pore-water profiles by sedimentary events and the significance of anaerobic oxidation of methane for the burial of sulfur in marine sediments, Geochim. Cosmochim. Ac., 67, 2631–2647,, 2003. 

Hensen, C., Zabel, M., and Schulz, H. N.: Benthic Cycling of Oxygen, Nitrogen and Phosphorus, in: Marine Geochemistry, edited by: Schulz, H. D. and Zabel, M., Springer-Verlag, Berlin/Heidelberg, 207–240, 6, 2006. 

Hilligsøe, K. M., Jensen, J. B., Ferdelman, T. G., Fossing, H., Lapham, L., Røy, H., and Jørgensen, B. B.: Methane fluxes in marine sediments quantified through core analyses and seismo-acoustic mapping (Bornholm Basin, Baltic Sea), Geochim. Cosmochim. Ac., 239, 255–274,, 2018. 

Ho, T. C. and Aris, R.: On apparent second-order kinetics, AIChE J., 33, 1050–1051,, 1987. 

Hoefs, M. J. L., Rijpstra, W. I. C., and Sinninghe Damsté, J. S.: The influence of oxic degradation on the sedimentary biomarker record I: evidence from Madeira Abyssal Plain turbidites, Geochim. Cosmochim. Ac., 66, 2719–2735,, 2002. 

Huguet, C., de Lange, G. J., Gustafsson, Ö., Middelburg, J. J., Sinninghe Damsté, J. S., and Schouten, S.: Selective preservation of soil organic matter in oxidized marine sediments (Madeira Abyssal Plain), Geochim. Cosmochim. Ac., 72, 6061–6068,, 2008. 

Hülse, D., Arndt, S., Daines, S., Regnier, P., and Ridgwell, A.: OMEN-SED 1.0: a novel, numerically efficient organic matter sediment diagenesis module for coupling to Earth system models, Geosci. Model Dev., 11, 2649–2689,, 2018. 

Hülse, D., Arndt, S., and Ridgwell, A.: Mitigation of Extreme Ocean Anoxic Event Conditions by Organic Matter Sulfurization, Paleoceanogr. Paleocl., 34, 476–489,, 2019. 

Jensen, M., Lomstein, E., and Sørensen, J.: Benthic NH4+ and NO3- flux following sedimentation of a spring phytoplankton bloom in Aarhus Bight, Denmark, Mar. Ecol. Prog. Ser., 61, 87–96,, 1990. 

Jonas, P. J. C. and Millward, G. E.: Metals and nutrients in the Severn Estuary and Bristol Channel: Contemporary inputs and distributions, Mar. Pollut. Bull., 61, 52–67,, 2010. 

Jørgensen, B. B. and Kasten, S.: Sulfur Cycling and Methane Oxidation, in: Marine Geochemistry, edited by: Schulz, H. D. and Zabel, M., Springer-Verlag, Berlin/Heidelberg, 271–309,, 2006. 

Jørgensen, B. B., Beulig, F., Egger, M., Petro, C., Scholze, C., and Røy, H.: Organoclastic sulfate reduction in the sulfate-methane transition of marine sediments, Geochim. Cosmochim. Ac., 254, 231–245,, 2019a. 

Jørgensen, B. B., Findlay, A. J., and Pellerin, A.: The Biogeochemical Sulfur Cycle of Marine Sediments, Front. Microbiol., 10, 849,, 2019b. 

Jourabchi, P.: Quantitative interpretation of pH distributions in aquatic sediments: A reaction-transport modeling approach, Am. J. Sci., 305, 919–956,, 2005. 

Karlson, B., Edler, L., Granéli, W., Sahlsten, E., and Kuylenstierna, M.: Subsurface chlorophyll maxima in the Skagerrak-processes and plankton community structure, J. Sea Res., 35, 139–158,, 1996. 

Kayler, Z. E., Premke, K., Gessler, A., Gessner, M. O., Griebler, C., Hilt, S., Klemedtsson, L., Kuzyakov, Y., Reichstein, M., Siemens, J., Totsche, K.-U., Tranvik, L., Wagner, A., Weitere, M., and Grossart, H.-P.: Integrating Aquatic and Terrestrial Perspectives to Improve Insights Into Organic Matter Cycling at the Landscape Scale, Front. Earth Sci., 7, 127,, 2019. 

Keil, R. G. and Cowie, G. L.: Organic matter preservation through the oxygen-deficient zone of the NE Arabian Sea as discerned by organic carbon:mineral surface area ratios, Mar. Geol., 161, 13–22,, 1999. 

Kharbush, J. J., Close, H. G., Van Mooy, B. A. S., Arnosti, C., Smittenberg, R. H., Le Moigne, F. A. C., Mollenhauer, G., Scholz-Böttcher, B., Obreht, I., Koch, B. P., Becker, K. W., Iversen, M. H., and Mohr, W.: Particulate Organic Carbon Deconstructed: Molecular and Chemical Composition of Particulate Organic Carbon in the Ocean, Front. Mar. Sci., 7, 518,, 2020. 

Knab, N. J., Cragg, B. A., Borowski, C., Parkes, R. J., Pancost, R., and Jørgensen, B. B.: Anaerobic oxidation of methane (AOM) in marine sediments from the Skagerrak (Denmark): I. Geochemical and microbiological analyses, Geochim. Cosmochim. Ac., 72, 2868–2879,, 2008. 

Koho, K. A., Nierop, K. G. J., Moodley, L., Middelburg, J. J., Pozzato, L., Soetaert, K., van der Plicht, J., and Reichart, G.-J.: Microbial bioavailability regulates organic matter preservation in marine sediments, Biogeosciences, 10, 1131–1141,, 2013. 

Krumins, V., Gehlen, M., Arndt, S., Van Cappellen, P., and Regnier, P.: Dissolved inorganic carbon and alkalinity fluxes from coastal marine sediments: model estimates for different shelf environments and sensitivity to global change, Biogeosciences, 10, 371–398,, 2013. 

Kusch, S., Rethemeyer, J., Ransby, D., and Mollenhauer, G.: Permafrost Organic Carbon Turnover and Export Into a High-Arctic Fjord: A Case Study From Svalbard Using Compound-specific 14C Analysis, J. Geophys. Res.-Biogeo., 126, e2020JG006008,, 2021. 

LaRowe, D. E., Arndt, S., Bradley, J. A., Burwicz, E., Dale, A. W., and Amend, J. P.: Organic carbon and microbial activity in marine sediments on a global scale throughout the Quaternary, Geochim. Cosmochim. Ac., 286, 227–247,, 2020a. 

LaRowe, D. E., Arndt, S., Bradley, J. A., Estes, E. R., Hoarfrost, A., Lang, S. Q., Lloyd, K. G., Mahmoudi, N., Orsi, W. D., Shah Walter, S. R., Steen, A. D., and Zhao, R.: The fate of organic carbon in marine sediments – New insights from recent data and analysis, Earth-Sci. Rev., 204, 103146,, 2020b. 

Latasa, M. and Bidigare, R. R.: A comparison of phytoplankton populations of the Arabian Sea during the Spring Intermonsoon and Southwest Monsoon of 1995 as described by HPLC-analyzed pigments, Deep-Sea Res. Pt. II, 45, 2133–2170,, 1998. 

Lengger, S. K., Rush, D., Mayser, J. P., Blewett, J., Schwartz-Narbonne, R., Talbot, H. M., Middelburg, J. J., Jetten, M. S. M., Schouten, S., Sinninghe Damsté, J. S., and Pancost, R. D.: Dark carbon fixation in the Arabian Sea oxygen minimum zone contributes to sedimentary organic carbon (SOM), Global Biogeochem. Cy., 2019GB006282,, 2019. 

Lessin, G., Artioli, Y., Almroth-Rosell, E., Blackford, J. C., Dale, A. W., Glud, R. N., Middelburg, J. J., Pastres, R., Queirós, A. M., Rabouille, C., Regnier, P., Soetaert, K., Solidoro, C., Stephens, N., and Yakushev, E.: Modelling Marine Sediment Biogeochemistry: Current Knowledge Gaps, Challenges, and Some Methodological Advice for Advancement, Front. Mar. Sci., 5, 1–8,, 2018. 

Lohse, L., Malschaert, J. F. P., Slomp, C. P., Helder, W., and van Raaphorst, W.: Sediment-water fluxes of inorganic nitrogen compounds along the transport route of organic matter in the North Sea, Ophelia, 41, 173–197,, 1995. 

Lomstein, E., Jensen, M. H., and Sørensen, J.: Intracellular NH4+ and NO3- pools associated with deposited phytoplankton in a marine sediment (Aarhus Bight, Denmark), Mar. Ecol. Prog. Ser., 61, 97–105, 1990. 

Luff, R., Wallmann, K., Grandel, S., and Schlüter, M.: Numerical modeling of benthic processes in the deep Arabian Sea, Deep-Sea Res. Pt. II, 47, 3039–3072,, 2000. 

Mackin, J. E. and Swider, K. T.: Organic matter decomposition pathways and oxygen consumption in coastal marine sediments, J. Mar. Res., 47, 681–716,, 1989. 

Maerz, J., Six, K. D., Stemmler, I., Ahmerkamp, S., and Ilyina, T.: Microstructure and composition of marine aggregates as co-determinants for vertical particulate organic carbon transfer in the global ocean, Biogeosciences, 17, 1765–1803,, 2020. 

Manning, A. J., Langston, W. J., and Jonas, P. J. C.: A review of sediment dynamics in the Severn Estuary: Influence of flocculation, Mar. Pollut. Bull., 61, 37–51,, 2010. 

Marquardt, M., Hensen, C., Piñero, E., Wallmann, K., and Haeckel, M.: A transfer function for the prediction of gas hydrate inventories in marine sediments, Biogeosciences, 7, 2925–2941,, 2010. 

Mayer, L. M.: Relationships between mineral surfaces and organic carbon concentrations in soils and sediments, Chem. Geol., 114, 347–363,, 1994a. 

Mayer, L. M.: Surface area control of organic carbon accumulation in continental shelf sediments, Geochim. Cosmochim. Ac., 58, 1271–1284,, 1994b. 

Mayer, L. M.: Sedimentary organic matter preservation: an assessment and speculative synthesis – a comment, Mar. Chem., 49, 123–126,, 1995. 

Meile, C. and Van Cappellen, P.: Particle age distributions and O 2 exposure times: Timescales in bioturbated sediments, Global Biogeochem. Cy., 19, 1–12,, 2005. 

Meister, P., Liu, B., Ferdelman, T. G., Jørgensen, B. B., and Khalili, A.: Control of sulphate and methane distributions in marine sediments by organic matter reactivity, Geochim. Cosmochim. Ac., 104, 183–193,, 2013. 

Mewes, K., Mogollón, J. M., Picard, A., Rühlemann, C., Eisenhauer, A., Kuhn, T., Ziebis, W., and Kasten, S.: Diffusive transfer of oxygen from seamount basaltic crust into overlying sediments: An example from the Clarion–Clipperton Fracture Zone, Earth Planet. Sc. Lett., 433, 215–225,, 2016. 

Meyers, P. A.: Organic geochemical proxies of paleoceanographic, paleolimnologic, and paleoclimatic processes, Org. Geochem., 27, 213–250,, 1997. 

Middelburg, J. J.: A simple rate model for organic matter decomposition in marine sediments, Geochim. Cosmochim. Ac., 53, 1577–1581,, 1989. 

Middelburg, J. J.: Reviews and syntheses: to the bottom of carbon processing at the seafloor, Biogeosciences, 15, 413–427,, 2018. 

Middelburg, J. J.: Marine Carbon Biogeochemistry: A Primer for Earth System Scientists, Springer International Publishing, Cham,, 2019. 

Middelburg, J. J., Vlug, T., Jaco, F., and van der Nat, W. A.: Organic matter mineralization in marine systems, Glob. Planet. Change, 8, 47–58,, 1993. 

Middelburg, J. J., Soetaert, K., and Herman, P. M. J.: Empirical relationships for use in global diagenetic models, Deep-Sea Res. Pt. I, 44, 327–344,, 1997. 

Mogollón, J. M., Dale, A. W., Fossing, H., and Regnier, P.: Timescales for the development of methanogenesis and free gas layers in recently-deposited sediments of Arkona Basin (Baltic Sea), Biogeosciences, 9, 1915–1933,, 2012. 

Mogollón, J. M., Mewes, K., and Kasten, S.: Quantifying manganese and nitrogen cycle coupling in manganese-rich, organic carbon-starved marine sediments: Examples from the Clarion-Clipperton fracture zone, Geophys. Res. Lett., 43, 7114–7123,, 2016. 

Mollenhauer, G., Eglinton, T. I., Ohkouchi, N., Schneider, R. R., Müller, P. J., Grootes, P. M., and Rullkötter, J.: Asynchronous alkenone and foraminifera records from the Benguela Upwelling System, Geochim. Cosmochim. Ac., 67, 2157–2171,, 2003. 

Mollenhauer, G., Inthorn, M., Vogt, T., Zabel, M., Sinninghe Damsté, J. S., and Eglinton, T. I.: Aging of marine organic matter during cross-shelf lateral transport in the Benguela upwelling system revealed by compound-specific radiocarbon dating, Geochem. Geophy. Geosy., 8, 1–16,, 2007. 

Müller, P. J. and Suess, E.: Productivity, sedimentation rate, and sedimentary organic matter in the oceans – I. Organic carbon preservation, Deep-Sea Res. Pt. A, 26, 1347–1362,, 1979. 

Niewöhner, C., Hensen, C., Kasten, S., Zabel, M., and Schulz, H. D.: Deep Sulfate Reduction Completely Mediated by Anaerobic Methane Oxidation in Sediments of the Upwelling Area off Namibia, Geochim. Cosmochim. Ac., 62, 455–464,, 1998. 

Ohkouchi, N., Eglinton, T. I., Keigwin, L. D., and Hayes, J. M.: Spatial and Temporal Offsets Between Proxy Records in a Sediment Drift, Science, 298, 1224,, 2002. 

Oni, O., Miyatake, T., Kasten, S., Richter-Heitmann, T., Fischer, D., Wagenknecht, L., Kulkarni, A., Blumers, M., Shylin, S. I., Ksenofontov, V., Costa, B. F. O., Klingelhöfer, G., and Friedrich, M. W.: Distinct microbial populations are tightly linked to the profile of dissolved iron in the methanic sediments of the Helgoland mud area, North Sea, Front. Microbiol., 6, 365,, 2015a. 

Oni, O. E., Schmidt, F., Miyatake, T., Kasten, S., Witt, M., Hinrichs, K.-U., and Friedrich, M. W.: Microbial Communities and Organic Matter Composition in Surface and Subsurface Sediments of the Helgoland Mud Area, North Sea, Front. Microbiol., 6, 1290,, 2015b. 

Owings, S. M., Bréthous, L., Eitel, E. M., Fields, B. P., Boever, A., Beckler, J. S., Bombled, B., Lansard, B., Metzger, E., Rabouille, C., and Taillefert, M.: Differential manganese and iron recycling and transport in continental margin sediments of the Northern Gulf of Mexico, Mar. Chem., 229, 103908,, 2021. 

Pastor, L., Deflandre, B., Viollier, E., Cathalot, C., Metzger, E., Rabouille, C., Escoubeyrou, K., Lloret, E., Pruski, A. M., Vétion, G., Desmalades, M., Buscail, R., and Grémare, A.: Influence of the organic matter composition on benthic oxygen demand in the Rhône River prodelta (NW Mediterranean Sea), Cont. Shelf Res., 31, 1008–1019,, 2011. 

Premuzic, E. T., Benkovitz, C. M., Gaffney, J. S., and Walsh, J. J.: The nature and distribution of organic matter in the surface sediments of world oceans and seas, Org. Geochem., 4, 63–77,, 1982. 

Pruski, A. M., Buscail, R., Bourgeois, S., Vétion, G., Coston-Guarini, J., and Rabouille, C.: Biogeochemistry of fatty acids in a river-dominated Mediterranean ecosystem (Rhône River prodelta, Gulf of Lions, France): Origins and diagenesis, Org. Geochem., 83/84, 227–240,, 2015. 

Puglini, M., Brovkin, V., Regnier, P., and Arndt, S.: Assessing the potential for non-turbulent methane escape from the East Siberian Arctic Shelf, Biogeosciences, 17, 3247–3275,, 2020. 

Ragavan, P. and Kumar, S.: Potential role of priming effect in the open ocean oxygen minimum zones: an outlook, Hydrobiologia, 848, 2437–2448,, 2021. 

Rassmann, J., Lansard, B., Pozzato, L., and Rabouille, C.: Carbonate chemistry in sediment porewaters of the Rhône River deltadriven by early diagenesis (northwestern Mediterranean), Biogeosciences, 13, 5379–5394,, 2016. 

Rassmann, J., Eitel, E. M., Lansard, B., Cathalot, C., Brandily, C., Taillefert, M., and Rabouille, C.: Benthic alkalinity and dissolved inorganic carbon fluxes in the Rhône River prodelta generated by decoupled aerobic and anaerobic processes, Biogeosciences, 17, 13–33,, 2020. 

Regnier, P., Mouchet, A., Wollast, R., and Ronday, F.: A discussion of methods for estimating residual fluxes in strong tidal estuaries, Cont. Shelf Res., 18, 1543–1571,, 1998. 

Regnier, P., O'Kane, J. P., Steefel, C. I., and Vanderborght, J. P.: Modeling complex multi-component reactive-transport systems: towards a simulation environment based on the concept of a Knowledge Base, Appl. Math. Model., 26, 913–927,, 2002. 

Regnier, P., Dale, A. W., Arndt, S., LaRowe, D. E., Mogollón, J., and Van Cappellen, P.: Quantitative analysis of anaerobic oxidation of methane (AOM) in marine sediments: A modeling perspective, Earth-Sci. Rev., 106, 105–130,, 2011. 

Riedinger, N., Pfeifer, K., Kasten, S., Garming, J. F. L., Vogt, C., and Hensen, C.: Diagenetic Alteration of Magnetic Signals by Anaerobic Oxidation of Methane Related to a Change in Sedimentation Rate, Geochim. Cosmochim. Ac., 69, 4117–4126,, 2005. 

Riedinger, N., Formolo, M. J., Lyons, T. W., Henkel, S., Beck, A., and Kasten, S.: An inorganic geochemical argument for coupled anaerobic oxidation of methane and iron reduction in marine sediments, Geobiology, 12, 172–181,, 2014. 

Riedinger, N., Brunner, B., Krastel, S., Arnold, G. L., Wehrmann, L. M., Formolo, M. J., Beck, A., Bates, S. M., Henkel, S., Kasten, S., and Lyons, T. W.: Sulfur Cycling in an Iron Oxide-Dominated, Dynamic Marine Depositional System: The Argentine Continental Margin, Front. Earth Sci., 5, 33, 1–19,, 2017. 

Rixen, T., Gaye, B., and Emeis, K.-C.: The monsoon, carbon fluxes, and the organic carbon pump in the northern Indian Ocean, Prog. Oceanogr., 175, 24–39,, 2019. 

Rysgaard, S., Fossing, H., and Jensen, M. M.: Organic matter degradation through oxygen respiration, denitrification, and manganese, iron, and sulfate reduction in marine sediments (the Kattegat and the Skagerrak), Ophelia, 55, 77–91,, 2001. 

Salvadó, J. A., Tesi, T., Andersson, A., Ingri, J., Dudarev, O. V., Semiletov, I. P., and Gustafsson, Ö.: Organic carbon remobilized from thawing permafrost is resequestered by reactive iron on the Eurasian Arctic Shelf, Geophys. Res. Lett., 42, 8122–8130,, 2015. 

Schlitzer, R.: Ocean data view, available at: (last access: 20 November 2020), 2019. 

Schmidt, M. W. I., Torn, M. S., Abiven, S., Dittmar, T., Guggenberger, G., Janssens, I. A., Kleber, M., Kögel-Knabner, I., Lehmann, J., Manning, D. A. C., Nannipieri, P., Rasse, D. P., Weiner, S., and Trumbore, S. E.: Persistence of soil organic matter as an ecosystem property, Nature, 478, 49–56,, 2011. 

Schulz, H. D., Dahmke, A., Schinzel, U., Wallmann, K., and Zabel, M.: Early diagenetic processes, fluxes, and reaction rates in sediments of the South Atlantic, Geochim. Cosmochim. Ac., 58, 2041–2060,, 1994. 

Seiter, K., Hensen, C., Schröter, J., and Zabel, M.: Organic carbon content in surface sediments – defining regional provinces, Deep Sea Res. Part Oceanogr. Res. Pap., 51, 2001–2026,, 2004. 

Seiter, K., Hensen, C., and Zabel, M.: Benthic carbon mineralization on a global scale, Global Biogeochem. Cy., 19, 1–26,, 2005. 

Shalapyonok, A., Olson, R. J., and Shalapyonok, L. S.: Arabian Sea phytoplankton during Southwest and Northeast Monsoons 1995: composition, size structure and biomass from individual cell properties measured by flow cytometry, Deep-Sea Res. Pt. II, 48, 1231–1261,, 2001. 

Shields, M. R., Bianchi, T. S., Gélinas, Y., Allison, M. A., and Twilley, R. R.: Enhanced terrestrial carbon preservation promoted by reactive iron in deltaic sediments, Geophys. Res. Lett., 43, 1149–1157,, 2016. 

Sinninghe Damsté, J. S., Irene, W., Rijpstra, C., de Leeuw, J. W., and Schenck, P. A.: Origin of organic sulphur compounds and sulphur-containing high molecular weight substances in sediments and immature crude oils, Org. Geochem., 13, 593–606,, 1988. 

Sinninghe Damsté, J. S., Kok, M. D., Köster, J., and Schouten, S.: Sulfurized carbohydrates: an important sedimentary sink for organic carbon?, Earth Planet. Sc. Lett., 164, 7–13,, 1998. 

Sinninghe Damsté, J. S., Rijpstra, W. I. C., and Reichart, G.: The influence of oxic degradation on the sedimentary biomarker record II. Evidence from Arabian Sea sediments, Geochim. Cosmochim. Ac., 66, 2737–2754,, 2002. 

Soetaert, K., Herman, P. M. J., and Middelburg, J. J.: A model of early diagenetic processes from the shelf to abyssal depths, Geochim. Cosmochim. Ac., 60, 1019–1040,, 1996. 

Soetaert, K., Herman, P. M. J., Middelburg, J. J., and Heip, C.: Assessing organic matter mineralization, degradability and mixing rate in an ocean margin sediment (Northeast Atlantic) by diagenetic modeling, J. Mar. Res., 56, 519–534,, 1998. 

Stabeno, P. J. and Hunt, G. L.: Overview of the Inner Front and Southeast Bering Sea Carrying Capacity Programs, Deep-Sea Res. Pt. II, 49, 6157–6168,, 2002. 

Stevenson, M. A., Faust, J. C., Andrade, L. L., Freitas, F. S., Gray, N. D., Tait, K., Hendry, K. R., Hilton, R. G., Henley, S. F., Tessin, A., Leary, P., Papadaki, S., Ford, A., März, C., and Abbott, G. D.: Transformation of organic matter in a Barents Sea sediment profile: coupled geochemical and microbiological processes, Philos. T. R. Soc. Math. Phys. Eng. Sci., 378, 20200223,, 2020. 

Stockwell, D. A., Whitledge, T. E., Zeeman, S. I., Coyle, K. O., Napp, J. M., Brodeur, R. D., Pinchuk, A. I., and Hunt, G. L.: Anomalous conditions in the south-eastern Bering Sea, 1997: nutrients, phytoplankton and zooplankton, Fish. Oceanogr., 10, 99–116,, 2001. 

Sun, M.-Y. and Wakeham, S. G.: Molecular evidence for degradation and preservation of organic matter in the anoxic Black Sea Basin, Geochim. Cosmochim. Ac., 58, 3395–3406,, 1994. 

Taylor, K. E.: Summarizing multiple aspects of model performance in a single diagram, J. Geophys. Res.-Atmos., 106, 7183–7192,, 2001. 

Thomas, S.: The response of anaerobic prokaryotic processes and communities in the Severn Estuary sediments to environmental changes, Cardiff University, Cardiff, 2014. 

Thullner, M., Dale, A. W., and Regnier, P.: Global-scale quantification of mineralization pathways in marine sediments: A reaction-transport modeling approach, Geochem. Geophy. Geosy., 10, 1–24,, 2009. 

Tromp, T. K., Van Cappellen, P., and Key, R. M.: A global model for the early diagenesis of organic carbon and organic phosphorus in marine sediments, Geochim. Cosmochim. Ac., 59, 1259–1284,, 1995. 

Uncles, R. J.: Physical properties and processes in the Bristol Channel and Severn Estuary, Mar. Pollut. Bull., 61, 5–20,, 2010. 

Van Cappellen, P. and Wang, Y.: Cycling of iron and manganese in surface sediments; a general theory for the coupled transport and reaction of carbon, oxygen, nitrogen, sulfur, iron, and manganese, Am. J. Sci., 296, 197–243,, 1996. 

van Nugteren, P., Moodley, L., Brummer, G.-J., Heip, C. H. R., Herman, P. M. J., and Middelburg, J. J.: Seafloor ecosystem functioning: the importance of organic matter priming, Mar. Biol., 156, 2277–2287,, 2009. 

Van Weering, T. C. E., Berger, G. W., and Kalf, J.: Recent sediment accumulation in the Skagerrak, Northeastern North Sea, Neth. J. Sea Res., 21, 177–189,, 1987. 

Vandewiele, S., Cowie, G., Soetaert, K., and Middelburg, J. J.: Amino acid biogeochemistry and organic matter degradation state across the Pakistan margin oxygen minimum zone, Deep-Sea Res. Pt. II, 56, 376–392,, 2009. 

Vasquez-Cardenas, D., Meysman, F. J. R., and Boschker, H. T. S.: A Cross-System Comparison of Dark Carbon Fixation in Coastal Sediments, Global Biogeochem. Cy., 34, e2019GB006298,, 2020. 

Wakeham, S. G., Lee, C., Hedges, J. I., Hernes, P. J., and Peterson, M. J.: Molecular indicators of diagenetic status in marine organic matter, Geochim. Cosmochim. Ac., 61, 5363–5369,, 1997. 

Wallmann, K., Aloisi, G., Haeckel, M., Obzhirov, A., Pavlova, G., and Tishchenko, P.: Kinetics of organic matter degradation, microbial methane generation, and gas hydrate formation in anoxic marine sediments, Geochim. Cosmochim. Ac., 70, 3905–3927,, 2006. 

Wang, D., Zhu, M.-X., Yang, G.-P., and Ma, W.-W.: Reactive Iron and Iron-Bound Organic Carbon in Surface Sediments of the River-Dominated Bohai Sea (China) Versus the Southern Yellow Sea, J. Geophys. Res.-Biogeo., 124, 79–98,, 2019. 

Wang, Y. and Van Cappellen, P.: A multicomponent reactive transport model of early diagenesis: Application to redox cycling in coastal marine sediments, Geochim. Cosmochim. Ac., 60, 2993–3014,, 1996. 

Weber, T., Cram, J. A., Leung, S. W., DeVries, T., and Deutsch, C.: Deep ocean nutrients imply large latitudinal variation in particle transfer efficiency, P. Natl. Acad. Sci. USA, 113, 8606–8611,, 2016. 

Wehrmann, L. M., Arndt, S., März, C., Ferdelman, T. G., and Brunner, B.: The evolution of early diagenetic signals in Bering Sea subseafloor sediments in response to varying organic carbon deposition over the last 4.3 Ma, Geochim. Cosmochim. Ac., 109, 175–196,, 2013. 

Westrich, J. T. and Berner, R. A.: The role of sedimentary organic matter in bacterial sulfate reduction: The G model tested, Limnol. Oceanogr., 29, 236–249,, 1984. 

Wiltshire, K. H. and Manly, B. F. J.: The warming trend at Helgoland Roads, North Sea: phytoplankton response, Helgol. Mar. Res., 58, 269–273,, 2004. 

Ye, H., Yang, T., Zhu, G., Jiang, S., and Wu, L.: Pore water geochemistry in shallow sediments from the northeastern continental slope of the South China sea, Mar. Pet. Geol., 75, 68–82,, 2016.  

Zebracki, M., Eyrolle-Boyer, F., Evrard, O., Claval, D., Mourier, B., Gairoard, S., Cagnat, X., and Antonelli, C.: Tracing the origin of suspended sediment in a large Mediterranean river by combining continuous river monitoring and measurement of artificial and natural radionuclides, Sci. Total Environ., 502, 122–132,, 2015. 

Zonneveld, K. A. F., Versteegh, G. J. M., Kasten, S., Eglinton, T. I., Emeis, K.-C., Huguet, C., Koch, B. P., de Lange, G. J., de Leeuw, J. W., Middelburg, J. J., Mollenhauer, G., Prahl, F. G., Rethemeyer, J., and Wakeham, S. G.: Selective preservation of organic matter in marine environments; processes and impact on the sedimentary record, Biogeosciences, 7, 483–511,, 2010. 

Short summary
It remains challenging to fully understand what controls carbon burial in marine sediments globally. Thus, we use a model–data approach to identify patterns of organic matter reactivity at the seafloor across distinct environmental conditions. Our findings support the notion that organic matter reactivity is a dynamic ecosystem property and strongly influences biogeochemical cycling and exchange. Our results are essential to improve predictions of future changes in carbon cycling and climate.
Final-revised paper