# Widespread energy limitation to life in global subseafloor sediments

J. Bradley, S. Arndt, J. Amend, E. Burwicz, A. Dale, M. Egger, D. Larowe

► To cite this version:

J. Bradley, S. Arndt, J. Amend, E. Burwicz, A. Dale, et al. Widespread energy limitation to life in global subseafloor sediments. Science Advances, 2020, 6 (32), 10.1126/sciadv.aba0697. hal-04450513

HAL Id: hal-04450513

https://hal.science/hal-04450513

Submitted on 21 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# Widespread energy limitation to life in global subseafloor sediments

J. Bradley, S. Arndt, J. Amend, E. Burwicz, A. Dale, M. Egger, D. Larowe

► To cite this version:

J. Bradley, S. Arndt, J. Amend, E. Burwicz, A. Dale, et al. Widespread energy limitation to life in global subseafloor sediments. Science Advances, 2020, 6 (32), 10.1126/sciadv.aba0697. hal-04450513

HAL Id: hal-04450513

https://hal.science/hal-04450513v1

Submitted on 21 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
ENVIRONMENTAL STUDIES

# Widespread energy limitation to life in global subseafloor sediments

J. A. Bradley<fnref n="1" />,²*, S. Arndt³, J. P. Amend⁴,⁵, E. Burwicz⁶, A. W. Dale⁶, M. Egger⁷, D. E. LaRowe⁴

Microbial cells buried in subseafloor sediments comprise a substantial portion of Earth's biosphere and control global biogeochemical cycles; however, the rate at which they use energy (i.e., power) is virtually unknown. Here, we quantify organic matter degradation and calculate the power utilization of microbial cells throughout Earth's Quaternary-age subseafloor sediments. Aerobic respiration, sulfate reduction, and methanogenesis mediate 6.9, 64.5, and 28.6% of global subseafloor organic matter degradation, respectively. The total power utilization of the subseafloor sediment biosphere is 37.3 gigawatts, less than 0.1% of the power produced in the marine photic zone. Aerobic heterotrophs use the largest share of global power (54.5%) with a median power utilization of 2.23 × 10⁻¹⁸ watts per cell, while sulfate reducers and methanogens use 1.08 × 10⁻¹⁹ and 1.50 × 10⁻²⁰ watts per cell, respectively. Most subseafloor cells subsist at energy fluxes lower than have previously been shown to support life, calling into question the power limit to life.

Copyright © 2020
The Authors, some
rights reserved;
exclusive licensee
American Association
for the Advancement
of Science. No claim to
original U.S. Government
Works. Distributed
under a Creative
Commons Attribution
License 4.0 (CC BY).

INTRODUCTION

Marine sediments harbor a vast number of microorganisms across a wide range of depositional settings (1). These microorganisms have been shown to survive for extraordinarily long periods of time (2). Contrary to life at Earth's surface, subsurface organisms are severely energy limited (3). Many of these sediment-dwelling microbes are novel and uncultured (4, 5) and survive under extreme energy limitation for millennia, thus calling into question the limit for life (2, 6). They are also responsible for degrading Earth's largest pool of organic C (7), exerting a major control on global climate and bio-geochemical cycles (8). The availability of energy to microorganisms in subseafloor sediments varies considerably from one environment to another (9) and imposes limits on their rates of activity, as well as their survival (3, 10). However, outside of measurements from specific sites (11–13) and laboratory experiments (14), the rate of energy utilization, i.e., power, of subseafloor life is not well known. Therefore, questions pertaining to microbial activity levels, community assembly, adaptive evolution, and even the lower power limit to life remain unanswered (3, 15, 16). Here, we quantify the power utilization for all microorganisms contained in subseafloor sediments that were deposited during the past 2.59 million years (i.e., the Quaternary period). A bioenergetic ecosystem model was developed and implemented on a 0.25° × 0.25° resolution global grid. Our model calculates the volumetric distribution, rate, and thermodynamic properties of particulate organic carbon (POC) degradation, as well as the global distribution of cells and electron acceptors, in Quaternary-age subseafloor sediments. On this basis, we calculate, in three dimensions, the cell-specific power of microbial life in the global marine sediment habitat.

RESULTS AND DISCUSSION

Metabolism and organic carbon degradation in global marine sediments

We designate POC degraded in oxic sediments to be mediated by O₂, POC degraded in the sulfate-reducing zone to be mediated by SO₄²⁻, and attribute POC degraded beneath the sulfate-methane transition to methanogenesis (reactions 1 to 3; table S1). The thickness and depth of Holocene- and Pleistocene-aged sediments, as well as the depth of oxygen penetration and the depth of the sulfate-methane transition, are highly variable across the seafloor (fig. S1). Oxic sediments are mostly restricted to deep-water abyssal zones (fig. S2). Together, we estimate that 2.7% of the Quaternary subseafloor sediment volume is oxic and that it contains 6.80 × 10³ Pg of organic C (4.6% of Quaternary subseafloor sediment POC). Our mechanistically derived estimate for the global mass of organic carbon in oxic marine sediment (6.80 × 10³ Pg of C) is very closely aligned with an independent empirically derived estimate (6.7 × 10³ to 16.0 × 10³ Pg of C) (17). Sulfate-reducing and methanogenic sediments are more widespread (33.0 and 64.3% of sediment volume, respectively), particularly underlying shallow-water shelf and margin settings, and in deeper and thicker sedimentary layers (fig. S2). We calculate that sulfate-reducing sediments contain 8.31 × 10⁴ Pg of organic C (56.9% of Quaternary subseafloor sediment POC) and methanogenic sediments contain 5.63 × 10⁴ Pg of organic C (38.5% of Quaternary subseafloor sediment POC). We find that rates of POC degradation are mostly reflective of the volumetric POC distribution across subseafloor sediments. The majority of POC degradation in Quaternary sediments globally is mediated by sulfate reduction (64.5%; 0.171 Pg of C year⁻<fnref n="1" />). We calculate that oxic sediments mediate 6.9% of global subseafloor POC degradation (0.018 Pg of C year⁻<fnref n="1" />), while methanogenic sediments mediate 28.6% (0.076 Pg of C year⁻<fnref n="1" />) (Fig. 1). Our results underpin the understanding that sulfate plays a major role in the global degradation of organic C in marine sediments (18–20), while simultaneously verifying the importance of methanogenesis (21).

Power utilization in global marine sediments

The flux of energy, i.e., power, that is supplied to and used by
microorganisms is calculated using the rates of POC degradation

<footnote n="1">¹School of Geography, Queen Mary University of London, London, UK. ²Department of Geochemistry, GFZ, German Research Centre for Geosciences, Potsdam, Germany. ³Department of Geosciences, Environment and Society, Université Libre de Bruxelles, Brussels, Belgium. ⁴Department of Earth Sciences, University of Southern California, Los Angeles, CA, USA. ⁵Department of Biological Sciences, University of Southern California, Los Angeles, CA, USA. ⁶GEOMAR, Helmholtz Centre for Ocean Research, Kiel, Germany. ⁷The Ocean Cleanup Foundation, Rotterdam, Netherlands.</footnote>

<footnote>*Corresponding author. Email: jbradley.earth@gmail.com</footnote>
# Widespread energy limitation to life in global subseafloor sediments

J. A. Bradley<fnref n="1" />,²*, S. Arndt³, J. P. Amend⁴,⁵, E. Burwicz⁶, A. W. Dale⁶, M. Egger⁷, D. E. LaRowe⁴

Microbial cells buried in subseafloor sediments comprise a substantial portion of Earth's biosphere and control global biogeochemical cycles; however, the rate at which they use energy (i.e., power) is virtually unknown. Here, we quantify organic matter degradation and calculate the power utilization of microbial cells throughout Earth's Quaternary-age subseafloor sediments. Aerobic respiration, sulfate reduction, and methanogenesis mediate 6.9, 64.5, and 28.6% of global subseafloor organic matter degradation, respectively. The total power utilization of the subseafloor sediment biosphere is 37.3 gigawatts, less than 0.1% of the power produced in the marine photic zone. Aerobic heterotrophs use the largest share of global power (54.5%) with a median power utilization of 2.23 × 10⁻¹⁸ watts per cell, while sulfate reducers and methanogens use 1.08 × 10⁻¹⁹ and 1.50 × 10⁻²⁰ watts per cell, respectively. Most subseafloor cells subsist at energy fluxes lower than have previously been shown to support life, calling into question the power limit to life.

Copyright © 2020
The Authors, some
rights reserved;
exclusive licensee
American Association
for the Advancement
of Science. No claim to
original U.S. Government
Works. Distributed
under a Creative
Commons Attribution
License 4.0 (CC BY).

## INTRODUCTION

Marine sediments harbor a vast number of microorganisms across a wide range of depositional settings (1). These microorganisms have been shown to survive for extraordinarily long periods of time (2). Contrary to life at Earth's surface, subsurface organisms are severely energy limited (3). Many of these sediment-dwelling microbes are novel and uncultured (4, 5) and survive under extreme energy limitation for millennia, thus calling into question the limit for life (2, 6). They are also responsible for degrading Earth's largest pool of organic C (7), exerting a major control on global climate and bio-geochemical cycles (8). The availability of energy to microorganisms in subseafloor sediments varies considerably from one environment to another (9) and imposes limits on their rates of activity, as well as their survival (3, 10). However, outside of measurements from specific sites (11–13) and laboratory experiments (14), the rate of energy utilization, i.e., power, of subseafloor life is not well known. Therefore, questions pertaining to microbial activity levels, community assembly, adaptive evolution, and even the lower power limit to life remain unanswered (3, 15, 16). Here, we quantify the power utilization for all microorganisms contained in subseafloor sediments that were deposited during the past 2.59 million years (i.e., the Quaternary period). A bioenergetic ecosystem model was developed and implemented on a 0.25° × 0.25° resolution global grid. Our model calculates the volumetric distribution, rate, and thermodynamic properties of particulate organic carbon (POC) degradation, as well as the global distribution of cells and electron acceptors, in Quaternary-age subseafloor sediments. On this basis, we calculate, in three dimensions, the cell-specific power of microbial life in the global marine sediment habitat.

## RESULTS AND DISCUSSION

### Metabolism and organic carbon degradation in global marine sediments

We designate POC degraded in oxic sediments to be mediated by O₂, POC degraded in the sulfate-reducing zone to be mediated by SO₄²⁻, and attribute POC degraded beneath the sulfate-methane transition to methanogenesis (reactions 1 to 3; table S1). The thickness and depth of Holocene- and Pleistocene-aged sediments, as well as the depth of oxygen penetration and the depth of the sulfate-methane transition, are highly variable across the seafloor (fig. S1). Oxic sediments are mostly restricted to deep-water abyssal zones (fig. S2). Together, we estimate that 2.7% of the Quaternary subseafloor sediment volume is oxic and that it contains 6.80 × 10³ Pg of organic C (4.6% of Quaternary subseafloor sediment POC). Our mechanistically derived estimate for the global mass of organic carbon in oxic marine sediment (6.80 × 10³ Pg of C) is very closely aligned with an independent empirically derived estimate (6.7 × 10³ to 16.0 × 10³ Pg of C) (17). Sulfate-reducing and methanogenic sediments are more widespread (33.0 and 64.3% of sediment volume, respectively), particularly underlying shallow-water shelf and margin settings, and in deeper and thicker sedimentary layers (fig. S2). We calculate that sulfate-reducing sediments contain 8.31 × 10⁴ Pg of organic C (56.9% of Quaternary subseafloor sediment POC) and methanogenic sediments contain 5.63 × 10⁴ Pg of organic C (38.5% of Quaternary subseafloor sediment POC). We find that rates of POC degradation are mostly reflective of the volumetric POC distribution across subseafloor sediments. The majority of POC degradation in Quaternary sediments globally is mediated by sulfate reduction (64.5%; 0.171 Pg of C year⁻<fnref n="1" />). We calculate that oxic sediments mediate 6.9% of global subseafloor POC degradation (0.018 Pg of C year⁻<fnref n="1" />), while methanogenic sediments mediate 28.6% (0.076 Pg of C year⁻<fnref n="1" />) (Fig. 1). Our results underpin the understanding that sulfate plays a major role in the global degradation of organic C in marine sediments (18–20), while simultaneously verifying the importance of methanogenesis (21).

### Power utilization in global marine sediments

The flux of energy, i.e., power, that is supplied to and used by microorganisms is calculated using the rates of POC degradation

<footnote n="1">¹School of Geography, Queen Mary University of London, London, UK. ²Department of Geochemistry, GFZ, German Research Centre for Geosciences, Potsdam, Germany. ³Department of Geosciences, Environment and Society, Université Libre de Bruxelles, Brussels, Belgium. ⁴Department of Earth Sciences, University of Southern California, Los Angeles, CA, USA. ⁵Department of Biological Sciences, University of Southern California, Los Angeles, CA, USA. ⁶GEOMAR, Helmholtz Centre for Ocean Research, Kiel, Germany. ⁷The Ocean Cleanup Foundation, Rotterdam, Netherlands.</footnote>

<footnote>*Corresponding author. Email: jbradley.earth@gmail.com</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Volumetric distribution and flux of sediment, cells, organic carbon, and energy across major marine sediment catabolic zones in global Quaternary subseafloor sediment. (A) The proportion of total sediment volume (cubic centimeters), number of microbial cells, POC stored (petagrams of C) and degraded (petagrams of C per year), and power used (gigawatts) in each catabolic zone in Quaternary sediment. (B) Frequency distribution of cell-specific power utilization in Quaternary-age sediment. Green denotes oxic, lilac signifies sulfate-reducing, and orange represents methanogenic sediment.</figcaption></figure>

combined with the Gibbs energies of O₂, SO₄²⁻, and CH₄-mediated catabolism under representative in situ conditions (tables S1 and S2). The combined total power used by all microorganisms in subseafloor sediment is 37.3 GW. This is equivalent to only 0.08% of the power produced photosynthetically in the marine photic zone (22, 23). Microbial power utilization is the highest, per cubic centimeter of sediment, in shallow coastal settings (up to 10⁻⁹ W cm⁻³), and it is uniformly low (~10⁻¹⁴ W cm⁻³) in the oldest Pleistocene sediments (2.59 million years) (fig. S2). High power (per cubic centimeter) in sediments underlying northern high latitudes and Southeast Asian shelves can be attributed to high concentrations of reactive POC (fig. S3) that degrade relatively rapidly.

Notably, we find that O₂-mediated POC degradation constitutes the major power source in subseafloor sediments (20.3 GW; 54.5% of global subseafloor power), despite oxic sediment comprising a minor fraction of total sediment volume and POC. Microorganisms in sulfate-reducing sediment use 14.6 GW of power globally (39.1% of the global total), while microorganisms in methanogenic sediments use only 2.4 GW, equivalent to only 6.4% of the global subseafloor power utilization. The dichotomy between both the volume and amount of POC stored within each sedimentary layer and the power realized by POC degradation for each catabolic strategy underscores the need to consider not only POC concentrations and other geochemical constituents but also the rate and energetics of major catabolic pathways occurring in the subseafloor biosphere.

## Cell-specific power utilization

The power per individual microbial cell is calculated by dividing the total power used by the number of microorganisms present in a sedimentary horizon. Our results confirm that subseafloor sediment microorganisms are characterized by extreme and widespread energy limitation. We calculate that 84.0% of microorganisms subsist at an energy flux less than the lowest previously calculated estimate of the basal power limit to life for natural settings [1.9 × 10⁻¹⁹ W per cell (13)]. The median rate of microbial energy utilization in Quaternary-age marine sediments, i.e., cell-specific power, is 3.32 × 10⁻²⁰ W per cell, which is more than two orders of magnitude lower than the lowest experimentally derived maintenance power for microorganisms

(24) and five orders of magnitude lower than the lowest culture-based estimates of maintenance power (25).

Stepwise decreases in microbial power utilization largely follow the redox state of the sediment (Fig. 1). The delineation of microbial power utilization according to catabolic zones is consistent with a continuous decrease in reactivity and degradation rate of organic matter with increasing sediment depth and sediment age (26). Microorganisms situated within oxic sediments exhibit the highest power per cell, with a median power utilization of 2.23 × 10⁻¹⁸ W per cell. This is broadly consistent with field data—derived estimates of the maintenance power of microorganisms from oxygen-bearing sediments [5.30 × 10⁻¹⁸ W per cell (6)]. The power utilization of microorganisms situated in anoxic sediments is generally lower than that in oxic sediments (27, 28). This is clearly reflected by our data: We find that cell-specific power in anoxic sediments is one to two orders of magnitude lower than that in oxic sediments, with a median power utilization of 1.08 × 10⁻¹⁹ and 1.50 × 10⁻²⁰ W per cell in sulfate-reducing and methanogenic sediments, respectively. Globally, we find that extremely few cells (<0.02%) subsist at less than 1 zW (10⁻²¹ W), substantiating a theoretically posed “minimum decay prevention power” for microbial life (13).

Distinct high- and low-power zones arise according to catabolic strategy and geography (Fig. 2). It can be shown that the somewhat multimodal distribution of cellular power utilization (Fig. 1) arises from log-normal distributions within each domain (i.e., shelf, margin, and abyssal settings) (fig. S4). Microorganisms contained in continental shelf sediments, particularly in northern high latitudes and near Southeast Asia, use comparatively higher power than elsewhere (Fig. 2). This highlights that sediment POC with relatively short depositional pathways (i.e., situated in shallow water, <200 m) is tilted toward higher reactivity (fig. S3) (29) and supports microorganisms at a higher activity state (Fig. 2). High cell-specific power is also prevalent in oxic sediments underlying the ocean rises of the East and Central Pacific, the Southwest and Southeast Indian Ridges, and the Mid-Atlantic Ridge (2000 to 3000 m). These regions are characterized by low cell abundance (1, 2), where microorganisms are sustained for millennia with extremely low mortality (2, 30) on ancient POC (17), but via a highly exergonic terminal electron
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Volumetric distribution and flux of sediment, cells, organic carbon, and energy across major marine sediment catabolic zones in global Quaternary subseafloor sediment. (A) The proportion of total sediment volume (cubic centimeters), number of microbial cells, POC stored (petagrams of C) and degraded (petagrams of C per year), and power used (gigawatts) in each catabolic zone in Quaternary sediment. (B) Frequency distribution of cell-specific power utilization in Quaternary-age sediment. Green denotes oxic, lilac signifies sulfate-reducing, and orange represents methanogenic sediment.</figcaption></figure>

combined with the Gibbs energies of O₂, SO₄²⁻, and CH₄-mediated catabolism under representative in situ conditions (tables S1 and S2). The combined total power used by all microorganisms in subseafloor sediment is 37.3 GW. This is equivalent to only 0.08% of the power produced photosynthetically in the marine photic zone (22, 23). Microbial power utilization is the highest, per cubic centimeter of sediment, in shallow coastal settings (up to 10⁻⁹ W cm⁻³), and it is uniformly low (~10⁻¹⁴ W cm⁻³) in the oldest Pleistocene sediments (2.59 million years) (fig. S2). High power (per cubic centimeter) in sediments underlying northern high latitudes and Southeast Asian shelves can be attributed to high concentrations of reactive POC (fig. S3) that degrade relatively rapidly.

Notably, we find that O₂-mediated POC degradation constitutes the major power source in subseafloor sediments (20.3 GW; 54.5% of global subseafloor power), despite oxic sediment comprising a minor fraction of total sediment volume and POC. Microorganisms in sulfate-reducing sediment use 14.6 GW of power globally (39.1% of the global total), while microorganisms in methanogenic sediments use only 2.4 GW, equivalent to only 6.4% of the global subseafloor power utilization. The dichotomy between both the volume and amount of POC stored within each sedimentary layer and the power realized by POC degradation for each catabolic strategy underscores the need to consider not only POC concentrations and other geochemical constituents but also the rate and energetics of major catabolic pathways occurring in the subseafloor biosphere.

## Cell-specific power utilization

The power per individual microbial cell is calculated by dividing the total power used by the number of microorganisms present in a sedimentary horizon. Our results confirm that subseafloor sediment microorganisms are characterized by extreme and widespread energy limitation. We calculate that 84.0% of microorganisms subsist at an energy flux less than the lowest previously calculated estimate of the basal power limit to life for natural settings [1.9 × 10⁻¹⁹ W per cell (13)]. The median rate of microbial energy utilization in Quaternary-age marine sediments, i.e., cell-specific power, is 3.32 × 10⁻²⁰ W per cell, which is more than two orders of magnitude lower than the lowest experimentally derived maintenance power for microorganisms

(24) and five orders of magnitude lower than the lowest culture-based estimates of maintenance power (25).

Stepwise decreases in microbial power utilization largely follow the redox state of the sediment (Fig. 1). The delineation of microbial power utilization according to catabolic zones is consistent with a continuous decrease in reactivity and degradation rate of organic matter with increasing sediment depth and sediment age (26). Microorganisms situated within oxic sediments exhibit the highest power per cell, with a median power utilization of 2.23 × 10⁻¹⁸ W per cell. This is broadly consistent with field data—derived estimates of the maintenance power of microorganisms from oxygen-bearing sediments [5.30 × 10⁻¹⁸ W per cell (6)]. The power utilization of microorganisms situated in anoxic sediments is generally lower than that in oxic sediments (27, 28). This is clearly reflected by our data: We find that cell-specific power in anoxic sediments is one to two orders of magnitude lower than that in oxic sediments, with a median power utilization of 1.08 × 10⁻¹⁹ and 1.50 × 10⁻²⁰ W per cell in sulfate-reducing and methanogenic sediments, respectively. Globally, we find that extremely few cells (<0.02%) subsist at less than 1 zW (10⁻²¹ W), substantiating a theoretically posed “minimum decay prevention power” for microbial life (13).

Distinct high- and low-power zones arise according to catabolic strategy and geography (Fig. 2). It can be shown that the somewhat multimodal distribution of cellular power utilization (Fig. 1) arises from log-normal distributions within each domain (i.e., shelf, margin, and abyssal settings) (fig. S4). Microorganisms contained in continental shelf sediments, particularly in northern high latitudes and near Southeast Asia, use comparatively higher power than elsewhere (Fig. 2). This highlights that sediment POC with relatively short depositional pathways (i.e., situated in shallow water, <200 m) is tilted toward higher reactivity (fig. S3) (29) and supports microorganisms at a higher activity state (Fig. 2). High cell-specific power is also prevalent in oxic sediments underlying the ocean rises of the East and Central Pacific, the Southwest and Southeast Indian Ridges, and the Mid-Atlantic Ridge (2000 to 3000 m). These regions are characterized by low cell abundance (1, 2), where microorganisms are sustained for millennia with extremely low mortality (2, 30) on ancient POC (17), but via a highly exergonic terminal electron
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Cell-specific power utilization across major marine sediment catabolic zones. Power per cell (watts) calculated on a global scale and depth-integrated for the (A) oxic, (B) sulfate-reducing, and (C) methanogenic sedimentary layers. White areas denote absence of the corresponding catabolic zone.</figcaption></figure>

transfer step (table S2). Sulfate-reducing sediments underlie many of the major ocean basins (>4000-m water depth), including the North Pacific, North Atlantic, Argentine, Canary, Angola, Cape, Southwest Pacific, and Mid-Indian Ocean basins. These sediments contain cells using very little power (<10⁻²⁰ W per cell). Methanogenesis, which is prevalent in shallow-water sediments near continental land masses, is shown to provide even less power to individual organisms (10⁻²¹ to 10⁻²⁰ W per cell). These extremely low-power settings, largely in deeper and thicker sedimentary layers underlying shallow-water shelf and margin settings, are typically characterized by high rates of mortality among subseafloor microbial communities (31). Nevertheless, the sustained presence of microbes, even in

deeper, low-power sedimentary layers, indicates that at least a fraction of the microbial population are capable of surviving at powers far lower than has been calculated for any natural environment previously (13).

Uncertainty in our calculations arises from (i) the rate of POC degradation, (ii) the Gibbs energy of the reaction, and (iii) the number of cells carrying out the reaction. Uncertainty is evaluated for each of these factors (i to iii) individually and concurrently (i.e., together) (see Supplementary Text). The reactivity of POC deposited onto the seafloor and its evolution during burial are accounted for by a set of adjustable parameters (see Materials and Methods) using a theoretically and empirically well-established reactive continuum modeling approach, validated using a large body of carefully dated organic matter and pore water depth profiles, rate measurements, isotopic signatures, and authigenic mineral distributions from a wide range of marine settings across contrasting depositional environments and temporal scales. We prescribe parameter values that are consistent with decades of observations, laboratory experiments, and modeling, in which the decrease in organic matter reactivity with increasing sediment depth and age, as well as from continental to abyssal environments, is widely observed (29). Nevertheless, to address the uncertainty in our results associated with (i) POC degradation rates, we evaluate a range of simulations encapsulating all possible variability in the reactivity of POC deposited onto the seafloor and its evolution during burial, constrained by a global parameter compilation (see Supplementary Text). To address the uncertainty associated with variability in (ii) the Gibbs energy arising from spatial heterogeneities in the concentrations of reactants and products, we run additional simulations, where the concentrations of the reactants and products are varied over realistic ranges to represent high-energy and low-energy scenarios. The calculation of power per cell also relies on an estimate of the number of cells carrying out the reaction, for which we use the best available spatially resolved global model for subseafloor biomass distribution (1) and simplify the reaction network to account for the most pertinent biogeochemical processes and the major catabolic pathways of POC oxidation. A large body of evidence suggests that the vast majority of microorganisms in subseafloor sediments are involved in the degradation of organic matter and overwhelmingly hold the genetic potential to degrade organic molecules, especially in organic-rich, coastal, and/or anoxic settings (see Supplementary Text). Nevertheless, to account for uncertainty in (iii) the number of microbial cells present in the sediment and carrying out the metabolism prescribed to a certain redox zone, we carry out additional simulations where cell abundance is varied by up to 50%. These factors (i to iii), tested individually and simultaneously, produce relatively minor deviations in the median-per-cell power (1.42 × 10⁻²⁰ to 10.20 × 10⁻²⁰ W per cell; see Supplementary Text). Under all tested scenarios, the percentage of microorganisms subsisting at less than 1.9 × 10⁻¹⁹ W per cell is at least 65.7%— underscoring the highly energy-limited nature of microbial life in subseafloor sediments.

Quantifying the energy turnover of the subseafloor biosphere is critical to determining the cell-specific minimum power requirement for life, as well as to predict the habitable boundaries of life on Earth. This study advances previous knowledge of the abundance, distribution, and activity of microorganisms in the subsurface (1, 32, 33) by providing the catabolic rates and specific activities of these microorganisms on a global scale. Considering that a small number of cells use a disproportionately large fraction of the total energy in
<figure><img src="image_3.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 2.** Cell-specific power utilization across major marine sediment catabolic zones. Power per cell (watts) calculated on a global scale and depth-integrated for the (A) oxic, (B) sulfate-reducing, and (C) methanogenic sedimentary layers. White areas denote absence of the corresponding catabolic zone.</figcaption></figure>

transfer step (table S2). Sulfate-reducing sediments underlie many of the major ocean basins (>4000-m water depth), including the North Pacific, North Atlantic, Argentine, Canary, Angola, Cape, Southwest Pacific, and Mid-Indian Ocean basins. These sediments contain cells using very little power (<10<sup>-20</sup> W per cell). Methanogenesis, which is prevalent in shallow-water sediments near continental land masses, is shown to provide even less power to individual organisms (10<sup>-21</sup> to 10<sup>-20</sup> W per cell). These extremely low-power settings, largely in deeper and thicker sedimentary layers underlying shallow-water shelf and margin settings, are typically characterized by high rates of mortality among subseafloor microbial communities (31). Nevertheless, the sustained presence of microbes, even in

deeper, low-power sedimentary layers, indicates that at least a fraction of the microbial population are capable of surviving at powers far lower than has been calculated for any natural environment previously (13).

Uncertainty in our calculations arises from (i) the rate of POC degradation, (ii) the Gibbs energy of the reaction, and (iii) the number of cells carrying out the reaction. Uncertainty is evaluated for each of these factors (i to iii) individually and concurrently (i.e., together) (see Supplementary Text). The reactivity of POC deposited onto the seafloor and its evolution during burial are accounted for by a set of adjustable parameters (see Materials and Methods) using a theoretically and empirically well-established reactive continuum modeling approach, validated using a large body of carefully dated organic matter and pore water depth profiles, rate measurements, isotopic signatures, and authigenic mineral distributions from a wide range of marine settings across contrasting depositional environments and temporal scales. We prescribe parameter values that are consistent with decades of observations, laboratory experiments, and modeling, in which the decrease in organic matter reactivity with increasing sediment depth and age, as well as from continental to abyssal environments, is widely observed (29). Nevertheless, to address the uncertainty in our results associated with (i) POC degradation rates, we evaluate a range of simulations encapsulating all possible variability in the reactivity of POC deposited onto the seafloor and its evolution during burial, constrained by a global parameter compilation (see Supplementary Text). To address the uncertainty associated with variability in (ii) the Gibbs energy arising from spatial heterogeneities in the concentrations of reactants and products, we run additional simulations, where the concentrations of the reactants and products are varied over realistic ranges to represent high-energy and low-energy scenarios. The calculation of power per cell also relies on an estimate of the number of cells carrying out the reaction, for which we use the best available spatially resolved global model for subseafloor biomass distribution (1) and simplify the reaction network to account for the most pertinent biogeochemical processes and the major catabolic pathways of POC oxidation. A large body of evidence suggests that the vast majority of microorganisms in subseafloor sediments are involved in the degradation of organic matter and overwhelmingly hold the genetic potential to degrade organic molecules, especially in organic-rich, coastal, and/or anoxic settings (see Supplementary Text). Nevertheless, to account for uncertainty in (iii) the number of microbial cells present in the sediment and carrying out the metabolism prescribed to a certain redox zone, we carry out additional simulations where cell abundance is varied by up to 50%. These factors (i to iii), tested individually and simultaneously, produce relatively minor deviations in the median-per-cell power (1.42 x 10<sup>-20</sup> to 10.20 x 10<sup>-20</sup> W per cell; see Supplementary Text). Under all tested scenarios, the percentage of microorganisms subsisting at less than 1.9 x 10<sup>-19</sup> W per cell is at least 65.7%— underscoring the highly energy-limited nature of microbial life in subseafloor sediments.

Quantifying the energy turnover of the subseafloor biosphere is critical to determining the cell-specific minimum power requirement for life, as well as to predict the habitable boundaries of life on Earth. This study advances previous knowledge of the abundance, distribution, and activity of microorganisms in the subsurface (1, 32, 33) by providing the catabolic rates and specific activities of these microorganisms on a global scale. Considering that a small number of cells use a disproportionately large fraction of the total energy in
marine sediments, we find that extraordinarily low power is a characteristic of almost all subseafloor sediment-hosted life. We report energy fluxes that fall short of previously calculated cell maintenance requirements for natural settings, verifying the overwhelmingly low-power nature of the deep biosphere (10). Considering that cells are subsisting at or near basal maintenance power requirements (13), it is unlikely that growth (i.e., cellular division), and by extension evolution, can be sustained widely throughout subseafloor sediments, validating what has shown to be the case for specific sites (15). The prevalence of dormant and nongrowing microorganisms, however, affords a survival advantage permitting ecological coexistence (34) and maintenance of biodiversity (35) in subseafloor sediments. Overall, our data represent evidence for widespread energy limitation in the marine subsurface, and that the overwhelming majority of subseafloor sediment microorganisms may subsist at close to the power limit for life.

## MATERIALS AND METHODS

### Experimental design

Here, we present and apply a new bioenergetic modeling approach that is used to quantify the distribution and rate of POC degradation, as well as the global distribution of cells, electron acceptors, and the thermodynamic properties of POC degradation. On this basis, we calculate the cell-specific power (i.e., energy per unit time per individual microbial cell), P (watts) according to

$$P = \frac{r \cdot \Delta G_r}{Y} \qquad (1)$$

where r denotes the rate of reaction (grams of C degraded per cubic centimeter per second), ΔG<sub>r</sub> (joules per gram of C) represents the Gibbs energy of the reaction, and Y (cells per cubic centimeter) is the number of microbial cells carrying out the reaction. The modeling framework is implemented on a 0.25° × 0.25° resolution global grid using the following global datasets and models.

### Reaction-transport model

We use a one-dimensional reaction-transport model (RTM) to calculate the distribution and degradation rate of POC in subseafloor sediments throughout the Quaternary period, following the approach described in (36). Similar to Wallmann et al. (37), constant boundary conditions and parameters are used to characterize the Holocene (0 to 11,700 years old) and Pleistocene (11,700 years old to 2.59 million years old) depositional environments.

The one-dimensional conservation equation for POC in porous media is given by (38–39)

$$\frac{\partial(1 - \Phi)\text{POC}}{\partial t} = \frac{\partial}{\partial z}\left(D_b(1 - \Phi)\frac{\partial\text{POC}}{\partial z}\right) - \frac{\partial(1 - \Phi)\omega\text{POC}}{\partial z} + (1 - \Phi)R_{\text{POC}} \qquad (2)$$

where POC (grams of C per cubic centimeter dry sediment) is the concentration of POC, t refers to time, D<sub>b</sub> (square centimeters per year) stands for the bioturbation coefficient, ω (centimeters per year) represents the sedimentation rate, and R<sub>POC</sub> (grams of C per cubic centimeter per year) denotes the rate of organic matter degradation. The porosity, Φ, of marine sediments in the shelf, margin, and abyss

domains was calculated as a function of depth, z (meters), assuming steady-state compaction, according to (40)

$$\Phi_z = \Phi_0 \cdot e^{-c_0 z} \qquad (3)$$

where Φ<sub>0</sub> denotes the porosity at the sediment-water interface (SWI) and c<sub>0</sub> (per meter) stands for the compaction length scale, which characterizes how a given sediment type will compact under its own weight.

Quaternary sediments are divided into three layers: Holocene bioturbated [top 10 cm (39)], nonbioturbated Holocene (<11,700 years), and Pleistocene (<2.59 million years). Sediment mixing was assumed to be constant over the bioturbated layer and nonexistent immediately below it. The rate of organic matter degradation, R<sub>POC</sub>, was described using a reactive continuum model (RCM). The RCM assumes a continuous yet dynamic distribution of organic compounds comprising a range of reactivities and reproducing the often-observed decrease in apparent POC reactivity with increasing sediment depth (and thus age) (41). Within the RCM, R<sub>POC</sub> is given by

$$R_{\text{POC}} = - \int_{0}^{\infty} k \cdot om(k, t) dk \qquad (4)$$

where om(k,t) represents a probability density function that determines the concentration of organic matter having a degradability between k and k + dk at time t, with k being analogous to a reaction rate constant. The initial distribution or organic compounds [om(k,0)] cannot be inferred by observations and may take different mathematical forms. We use a gamma function (36, 41–43), assuming first-order degradation kinetics, whereby the initial (t = 0) distribution of om over k is given by

$$om(k, 0) = \frac{\text{POC}_0 \cdot a^v \cdot k^{v-1} \cdot e^{-a \cdot k}}{\Gamma(v)} \qquad (5)$$

where POC<sub>0</sub> is the initial organic matter content (at the SWI), Γ is the gamma function, a (years) is the average lifetime of the reactive components of the POC, and v is a dimensionless parameter determining the shape of the distribution near k = 0. Assuming steady-state conditions (∂POC/∂t = 0) and a known organic C content at the sediment water interface, POC<sub>0</sub>, the change in the POC concentration as a function of depth, POC(z), is given by (41)

$$\text{POC}(z) = \text{POC}_0 \cdot \left( \frac{a}{a + \text{age}(z)} \right)^v \qquad (6)$$

where age(z) refers to the age of the sediment layer at depth z. Following the approaches described in (36, 44), we use a multi-G approximation of the RCM for the bioturbated sediment. Within the bioturbated zone, POC is represented by 500 distinct fractions that are degraded according to a first-order organic matter degradation rate law with a degradation rate constant, k<sub>i</sub>

$$R_{\text{POC}} = \sum_{i=1}^{500} k_i \cdot \text{POC}_i(z) \qquad (7)$$

where

$$\text{POC}_i(0) = F_i \cdot \text{POC} \qquad (8)$$

The initial proportion of total organic matter in fraction i, F<sub>i</sub>, as well as its respective reactivity, k<sub>i</sub>, can be determined through the initial probability density function that determines the concentration
marine sediments, we find that extraordinarily low power is a characteristic of almost all subseafloor sediment-hosted life. We report energy fluxes that fall short of previously calculated cell maintenance requirements for natural settings, verifying the overwhelmingly low-power nature of the deep biosphere (10). Considering that cells are subsisting at or near basal maintenance power requirements (13), it is unlikely that growth (i.e., cellular division), and by extension evolution, can be sustained widely throughout subseafloor sediments, validating what has shown to be the case for specific sites (15). The prevalence of dormant and nongrowing microorganisms, however, affords a survival advantage permitting ecological coexistence (34) and maintenance of biodiversity (35) in subseafloor sediments. Overall, our data represent evidence for widespread energy limitation in the marine subsurface, and that the overwhelming majority of subseafloor sediment microorganisms may subsist at close to the power limit for life.

## MATERIALS AND METHODS

### Experimental design

Here, we present and apply a new bioenergetic modeling approach that is used to quantify the distribution and rate of POC degradation, as well as the global distribution of cells, electron acceptors, and the thermodynamic properties of POC degradation. On this basis, we calculate the cell-specific power (i.e., energy per unit time per individual microbial cell), P (watts) according to

$$P = \frac{r \cdot \Delta G_r}{Y} \qquad (1)$$

where r denotes the rate of reaction (grams of C degraded per cubic centimeter per second), ΔG<sub>r</sub> (joules per gram of C) represents the Gibbs energy of the reaction, and Y (cells per cubic centimeter) is the number of microbial cells carrying out the reaction. The modeling framework is implemented on a 0.25° × 0.25° resolution global grid using the following global datasets and models.

### Reaction-transport model

We use a one-dimensional reaction-transport model (RTM) to calculate the distribution and degradation rate of POC in subseafloor sediments throughout the Quaternary period, following the approach described in (36). Similar to Wallmann et al. (37), constant boundary conditions and parameters are used to characterize the Holocene (0 to 11,700 years old) and Pleistocene (11,700 years old to 2.59 million years old) depositional environments.

The one-dimensional conservation equation for POC in porous media is given by (38–39)

$$\frac{\partial(1 - \Phi)\text{POC}}{\partial t} = \frac{\partial}{\partial z}\left(D_b(1 - \Phi)\frac{\partial\text{POC}}{\partial z}\right) - \frac{\partial(1 - \Phi)\omega\text{POC}}{\partial z} + (1 - \Phi)R_{\text{POC}} \qquad (2)$$

where POC (grams of C per cubic centimeter dry sediment) is the concentration of POC, t refers to time, D<sub>b</sub> (square centimeters per year) stands for the bioturbation coefficient, ω (centimeters per year) represents the sedimentation rate, and R<sub>POC</sub> (grams of C per cubic centimeter per year) denotes the rate of organic matter degradation. The porosity, Φ, of marine sediments in the shelf, margin, and abyss

domains was calculated as a function of depth, z (meters), assuming steady-state compaction, according to (40)

$$\Phi_z = \Phi_0 \cdot e^{-c_0 z} \qquad (3)$$

where Φ<sub>0</sub> denotes the porosity at the sediment-water interface (SWI) and c<sub>0</sub> (per meter) stands for the compaction length scale, which characterizes how a given sediment type will compact under its own weight.

Quaternary sediments are divided into three layers: Holocene bioturbated [top 10 cm (39)], nonbioturbated Holocene (<11,700 years), and Pleistocene (<2.59 million years). Sediment mixing was assumed to be constant over the bioturbated layer and nonexistent immediately below it. The rate of organic matter degradation, R<sub>POC</sub>, was described using a reactive continuum model (RCM). The RCM assumes a continuous yet dynamic distribution of organic compounds comprising a range of reactivities and reproducing the often-observed decrease in apparent POC reactivity with increasing sediment depth (and thus age) (41). Within the RCM, R<sub>POC</sub> is given by

$$R_{\text{POC}} = - \int_{0}^{\infty} k \cdot om(k, t) dk \qquad (4)$$

where om(k,t) represents a probability density function that determines the concentration of organic matter having a degradability between k and k + dk at time t, with k being analogous to a reaction rate constant. The initial distribution or organic compounds [om(k,0)] cannot be inferred by observations and may take different mathematical forms. We use a gamma function (36, 41–43), assuming first-order degradation kinetics, whereby the initial (t = 0) distribution of om over k is given by

$$om(k, 0) = \frac{\text{POC}_0 \cdot a^v \cdot k^{v-1} \cdot e^{-a \cdot k}}{\Gamma(v)} \qquad (5)$$

where POC<sub>0</sub> is the initial organic matter content (at the SWI), Γ is the gamma function, a (years) is the average lifetime of the reactive components of the POC, and v is a dimensionless parameter determining the shape of the distribution near k = 0. Assuming steady-state conditions (∂POC/∂t = 0) and a known organic C content at the sediment water interface, POC<sub>0</sub>, the change in the POC concentration as a function of depth, POC(z), is given by (41)

$$\text{POC}(z) = \text{POC}_0 \cdot \left( \frac{a}{a + \text{age}(z)} \right)^v \qquad (6)$$

where age(z) refers to the age of the sediment layer at depth z. Following the approaches described in (36, 44), we use a multi-G approximation of the RCM for the bioturbated sediment. Within the bioturbated zone, POC is represented by 500 distinct fractions that are degraded according to a first-order organic matter degradation rate law with a degradation rate constant, k<sub>i</sub>

$$R_{\text{POC}} = \sum_{i=1}^{500} k_i \cdot \text{POC}_i(z) \qquad (7)$$

where

$$\text{POC}_i(0) = F_i \cdot \text{POC} \qquad (8)$$

The initial proportion of total organic matter in fraction i, F<sub>i</sub>, as well as its respective reactivity, k<sub>i</sub>, can be determined through the initial probability density function that determines the concentration
of organic matter having a degradability between k and k + dk at time 0 (Eq. 4). The initial fraction of total POC characterized by a distinct reactivity k is given by

$$f(k, 0) = \frac{om(k, 0)}{POC_0} = \frac{ia^v \cdot k^{v-1} \cdot e^{-a \cdot k}}{\Gamma(v)} \quad (9)$$

The initial fraction of POC within the reactivity range between 0 and k, i.e., having a reactivity ≤k, is then given by integrating Eq. 9, assuming a, v, and k > 0

$$\begin{aligned} F(k, 0) &= \int_{0}^{k} f(0, k) dk = \int_{0}^{k} \frac{a^v \cdot k^{v-1} \cdot e^{-a \cdot k}}{\Gamma(v)} dk \\ &= \frac{a^v \cdot k^v \cdot (a \cdot k)^{-v} (\Gamma(v) - \Gamma(v, a \cdot k))}{\Gamma(v)} \quad (10) \\ &= \left( \frac{1 - \Gamma(v, a \cdot k)}{\Gamma(v)} \right) \end{aligned}$$

where Γ(v,a·k) denotes the inverse gamma function.

In the bioturbated sediment, the RCM is approximated by dividing the reactivity range k = [10⁻¹⁵, 10⁻(log(a)+2)] into 500 equal reactivity bins, kᵢ, thus ensuring a comprehensive approximation of the gamma function defined by the respective a and v values. The initial fraction, Fᵢ, of total POC within the reactivity bin kᵢ-1 and kᵢ (and thus with reactivity kᵢ = kᵢ-1 + (kᵢ - kᵢ-1)/2) in the 500G model can then be calculated as

$$F_i = F(k_j, 0) - F(k_{j-1}, 0) \quad (11)$$

The most reactive fraction, F₅₀₀, with reactivity k₅₀₀ = 10⁻log(a)+2 year⁻¹ is calculated on the basis of the upper incomplete gamma function

$$F_{500} = \int_{k_{500}}^{\infty} f(k_{500}, 0) dk = \frac{\Gamma(v, a \cdot k_{500})}{\Gamma(v)} \quad (12)$$

The derived rate constants are used in Eq. 2 by expressing RₚOC according to Eq. 7 to determine POC concentrations and degradation rates in the Holocene bioturbated layer (<10 cm). Assuming steady-state conditions, the general solution of Eq. 2 for each organic matter fraction i is given by

$$\mathrm{POC}_i(z) = A_i e^{(\alpha_i z)} + B_i e^{(b_i z)} \quad (13)$$

where

$$\alpha_i = \frac{\omega - \sqrt{\omega^2 + 4D_b \cdot k_i}}{2D_b} \quad (14)$$

$$b_i = \frac{\omega + \sqrt{\omega^2 + 4D_b \cdot k_i}}{2D_b} \quad (15)$$

and

$$\mathrm{POC}(z) = \sum_{i=1}^{500} \mathrm{POC}_i(z) \quad (16)$$

The integration constants Aᵢ and Bᵢ are defined by the chosen boundary conditions. Here, we apply a known POC concentration at the sediment-water interface [POC(0) = POC₀] and assume

continuity (equal flux and concentration) across the bottom of the bioturbated layer, i.e.

$$\mathrm{POC}(\mathrm{zbio}) = \mathrm{POC}_{\mathrm{zbio}}; \quad \left. \frac{-D_b}{dz} \mathrm{POC} \right|_{\mathrm{zbio}} = 0 \quad (17)$$

Below the Holocene bioturbated zone, the values of age(z) that are required to evaluate Eq. 6 are calculated using burial rates, ω(z), porosity depth profiles, Φ(z), and the apparent age of organic matter at the lower limit of the Holocene bioturbated zone, age<sub>zbio</sub>. By inserting POC<sub>bio</sub> and age<sub>zbio</sub> into Eq. 6 for POC(z) and age(z), respectively, age<sub>zbio</sub> can be solved

$$\mathrm{age}_{\mathrm{zbio}} = \frac{-a \cdot (\exp(\ln(\mathrm{POC}_{\mathrm{bio}}/\mathrm{POC}_0)/v) - 1)}{\exp(\ln(\mathrm{POC}_{\mathrm{bio}}/\mathrm{POC}_0)/v)} \quad (18)$$

Assuming an exponentially decreasing porosity, Eq. 2, and steady-state compaction, the burial velocity, ω, at depth z is then [e.g., (38)]

$$\omega(z) = \left( \frac{1 - \Phi_0}{1 - \Phi(z)} \right) \omega_0 \quad (19)$$

where ω₀ corresponds to the burial velocity at the SWI. The age of a given sediment layer at depth z below the Holocene bioturbated zone, age(z), is given by [e.g., (38)]

$$\mathrm{age}(z) = \int_{0}^{z} \omega^{-1} dz \quad (20)$$

Substituting Eq. 19 into Eq. 20 results in

$$\mathrm{age}(z) = \frac{1}{(1 - \Phi_0) \omega_0} \int_{0}^{z} (1 - \Phi) dz \quad (21)$$

which, upon integration, leads to

$$\mathrm{age}(z) = \frac{z + \frac{\Phi_0}{c_0} \cdot (\exp(-c_0 \cdot z) - 1)}{\omega_0 \cdot (1 - \Phi_0)} \quad (22)$$

The age of POC below the Holocene bioturbated zone is thus given by

$$\mathrm{age}(z) = \mathrm{age}_{\mathrm{zbio}} + \frac{z + \frac{\Phi_0}{c_0} \cdot (\exp(-c_0 \cdot z) - 1)}{\omega_0 \cdot (1 - \Phi_0)} \quad (23)$$

The depth distribution of organic matter in marine sediments deposited since the beginning of the Quaternary period can thus be calculated with knowledge of the sedimentation rate, level of bioturbation, porosity structure, bulk organic matter concentration at the SWI, and the distribution of organic compounds across the reactivity range at the sediment-water interface.

## Total POC budget

The fraction of total POC preserved in a layer of sediment that accumulated over a given time interval, PEᵢ, is given as the ratio of the total amount of POC stored in the ith sediment layer [i = Holocene bioturbated (0 to z<sub>bio</sub>), Holocene nonbioturbated (z<sub>bio</sub> to z<sub>holo</sub>), and Pleistocene (z<sub>holo</sub> to z<sub>pleis</sub>)], POCᵢ (grams of C per square centimeter),
of organic matter having a degradability between k and k + dk at time 0 (Eq. 4). The initial fraction of total POC characterized by a distinct reactivity k is given by

$$f(k, 0) = \frac{om(k, 0)}{POC_0} = \frac{ia^v \cdot k^{v-1} \cdot e^{-a \cdot k}}{\Gamma(v)} \quad (9)$$

The initial fraction of POC within the reactivity range between 0 and k, i.e., having a reactivity ≤k, is then given by integrating Eq. 9, assuming a, v, and k > 0

$$\begin{aligned} F(k, 0) &= \int_{0}^{k} f(0, k) dk = \int_{0}^{k} \frac{a^v \cdot k^{v-1} \cdot e^{-a \cdot k}}{\Gamma(v)} dk \\ &= \frac{a^v \cdot k^v \cdot (a \cdot k)^{-v} (\Gamma(v) - \Gamma(v, a \cdot k))}{\Gamma(v)} \quad (10) \\ &= \left( \frac{1 - \Gamma(v, a \cdot k)}{\Gamma(v)} \right) \end{aligned}$$

where $\Gamma(v, a \cdot k)$ denotes the inverse gamma function.

In the bioturbated sediment, the RCM is approximated by dividing the reactivity range $k = [10^{-15}, 10^{-(\log(a)+2)}]$ into 500 equal reactivity bins, $k_j$, thus ensuring a comprehensive approximation of the gamma function defined by the respective a and v values. The initial fraction, $F_i$, of total POC within the reactivity bin $k_{j-1}$ and $k_j$ (and thus with reactivity $k_i = k_{j-1} + (k_j - k_{j-1})/2$) in the 500G model can then be calculated as

$$F_i = F(k_j, 0) - F(k_{j-1}, 0) \quad (11)$$

The most reactive fraction, $F_{500}$, with reactivity $k_{500} = 10^{-\log(a)+2}$ year$^{-1}$ is calculated on the basis of the upper incomplete gamma function

$$F_{500} = \int_{k_{500}}^{\infty} f(k_{500}, 0) dk = \frac{\Gamma(v, a \cdot k_{500})}{\Gamma(v)} \quad (12)$$

The derived rate constants are used in Eq. 2 by expressing $R_{POC}$ according to Eq. 7 to determine POC concentrations and degradation rates in the Holocene bioturbated layer (<10 cm). Assuming steady-state conditions, the general solution of Eq. 2 for each organic matter fraction i is given by

$$POC_i(z) = A_i e^{(\alpha_i z)} + B_i e^{(b_i z)} \quad (13)$$

where

$$\alpha_i = \frac{\omega - \sqrt{\omega^2 + 4D_b \cdot k_i}}{2D_b} \quad (14)$$

$$b_i = \frac{\omega + \sqrt{\omega^2 + 4D_b \cdot k_i}}{2D_b} \quad (15)$$

and

$$POC(z) = \sum_{i=1}^{500} POC_i(z) \quad (16)$$

The integration constants $A_i$ and $B_i$ are defined by the chosen boundary conditions. Here, we apply a known POC concentration at the sediment-water interface [POC(0) = POC$_0$] and assume

continuity (equal flux and concentration) across the bottom of the bioturbated layer, i.e.

$$\mathrm{POC}(\mathrm{zbio}) = \mathrm{POC}_{\mathrm{zbio}}; \left. \frac{-D_b d\mathrm{POC}}{dz} \right|_{\mathrm{zbio}} = 0 \quad (17)$$

Below the Holocene bioturbated zone, the values of age(z) that are required to evaluate Eq. 6 are calculated using burial rates, $\omega(z)$, porosity depth profiles, $\Phi(z)$, and the apparent age of organic matter at the lower limit of the Holocene bioturbated zone, age$_{zbio}$. By inserting POC$_{bio}$ and age$_{zbio}$ into Eq. 6 for POC(z) and age(z), respectively, age$_{zbio}$ can be solved

$$\mathrm{age}_{\mathrm{zbio}} = \frac{-a \cdot (\exp(\ln(\mathrm{POC}_{\mathrm{bio}}/\mathrm{POC}_0)/v) - 1)}{\exp(\ln(\mathrm{POC}_{\mathrm{bio}}/\mathrm{POC}_0)/v)} \quad (18)$$

Assuming an exponentially decreasing porosity, Eq. 2, and steady-state compaction, the burial velocity, $\omega$, at depth z is then [e.g., (38)]

$$\omega(z) = \left( \frac{1 - \Phi_0}{1 - \Phi(z)} \right) \omega_0 \quad (19)$$

where $\omega_0$ corresponds to the burial velocity at the SWI. The age of a given sediment layer at depth z below the Holocene bioturbated zone, age(z), is given by [e.g., (38)]

$$\mathrm{age}(z) = \int_{0}^{z} \omega^{-1} dz \quad (20)$$

Substituting Eq. 19 into Eq. 20 results in

$$\mathrm{age}(z) = \frac{1}{(1 - \Phi_0) \omega_0} \int_{0}^{z} (1 - \Phi) dz \quad (21)$$

which, upon integration, leads to

$$\mathrm{age}(z) = \frac{z + \frac{\Phi_0}{c_0} \cdot (\exp(-c_0 \cdot z) - 1)}{\omega_0 \cdot (1 - \Phi_0)} \quad (22)$$

The age of POC below the Holocene bioturbated zone is thus given by

$$\mathrm{age}(z) = \mathrm{age}_{\mathrm{zbio}} + \frac{z + \frac{\Phi_0}{c_0} \cdot (\exp(-c_0 \cdot z) - 1)}{\omega_0 \cdot (1 - \Phi_0)} \quad (23)$$

The depth distribution of organic matter in marine sediments deposited since the beginning of the Quaternary period can thus be calculated with knowledge of the sedimentation rate, level of bioturbation, porosity structure, bulk organic matter concentration at the SWI, and the distribution of organic compounds across the reactivity range at the sediment-water interface.

## Total POC budget

The fraction of total POC preserved in a layer of sediment that accumulated over a given time interval, $\overline{PE}_i$, is given as the ratio of the total amount of POC stored in the $i$th sediment layer [$i = Holocene bioturbated (0 to z_{bio}), Holocene nonbioturbated (z_{bio} to z_{holo}), and Pleistocene (z_{holo} to z_{pleis})$], $\overline{POC}_i$ (grams of C per square centimeter),
and the total steady-state input of POC to that respective layer, $\overline{T}_i$, (grams of C per square centimeter)

$$ \overline{PE}_i = \overline{POC}_i / \overline{T}_i \qquad (24) $$

where

$$ \overline{POC}_i = \int_{z_i}^{z_{i-1}} POC(z) dz \qquad (25) $$

and

$$ \overline{T}_i = POC(z_i - 1) \cdot \Delta z_i \qquad (26) $$

The amount of POC degraded in the ith layer, $\overline{R}_i$, is given by

$$ \overline{R}_i = \int_{z_i}^{z_{i-1}} k(z) \cdot POC(z) dz \qquad (27) $$

Two different global datasets (0.25° × 0.25°) of values for $\omega$, the sedimentation rate, are used: one for Holocene bioturbated and Holocene sediments and the other for Pleistocene sediments.

## Parameters and forcings

The concentration of Holocene and Pleistocene POC at the sediment-water interface, $POC_0$, and sedimentation rates, $\omega$, are constrained according to (37), which used data from (45, 46) and an algorithm that correlates water depth and sedimentation rate (47). Following (36), the bioturbation coefficient is calculated as a function of water depth based on a compilation of empirical data collected in (48). Its values range spatially from 0.59 to 27 cm² year⁻¹, decreasing in magnitude as water depth increases. It is constant throughout the depth of Holocene bioturbated zone and immediately drops to zero beneath it. Values of $\Phi_0$ and $c_0$ are chosen to describe the shelf, margin, and abyss based on sediments that are representative of these domains (see table S3) (49). We partitioned the ocean floor into shelf, margin, and abyss domains to specify values for some of the model parameters (fig. S5 and table S3). We choose a constant v parameter of 0.125 for all three sediment domains, which is characteristic of fresh organic matter and based on observations that the v parameter values do not vary much between sites (29, 41). We link values of the a parameter with sedimentation rates, based on global observations (29, 41). This approach accounts for order-of-magnitude changes in a due to factors that control organic matter transit times from its source to deposition. These parameters thus reflect typically observed RCM parameter variability across various depositional environments. Uncertainties arising from parameter values are addressed in Supplementary Text.

## Global reaction network and catabolic strategy

We designate all sediments shallower than the maximum $O_2$ penetration depth as oxic. To map $O_2$ penetration depth in marine sediments on a global scale (fig. S1), we combine (i) global datasets of oxygen concentrations in marine sediments and (ii) previous modeling results (2). First, we use a global dataset of oxygen concentrations and corresponding water depths (table S4) to formulate a regression model [$R^2$ = 0.64; based on (50)]. This model is extrapolated over the entire ocean floor. For data in columns 81 to 91 (table S4), oxygen penetration depth was determined following the methods described in (51). Second, we superimpose (onto our global map) the results from a modeling study (2), which designated regions of the seafloor where $O_2$ is known to penetrate through all sedimentary

layers to basement rock. We define the sulfate reduction zone as the horizon between the maximum $O_2$ penetration depth and the depth of the sulfate-methane transition. We use data from (52) to map the global sulfate-methane transition depth onto our 0.25° × 0.25° global grid. The reaction zones designated within our model are based on the best available published datasets and include not only global compilations of marine sediment core data (50, 52) but also global models (2, 52). Data from the numerous core samples (n = 1704) documenting the distribution of sulfate-reducing and methanogenic sediments included in our study are available via (52) and the International Ocean Discovery Program. We designate POC degraded in oxic sediments to be mediated by $O_2$, POC degraded in the sulfate-reducing zone to be mediated by $SO_4^{2-}$, and prescribe POC degraded beneath the sulfate-methane transition to methanogenesis (reactions 1 to 3; table S1), thus simplifying the reaction network to account for the most pertinent biogeochemical process, based on widely accepted modeling and experimental studies (18, 21, 53). Although other oxidants are used by microorganisms for degrading organic matter in marine sediments, such as nitrate, and metal oxides, these have been shown to be quantitatively unimportant metabolic pathways on a global scale (29, 44, 54–56)—mainly because of their low concentration compared to $SO_4^{2-}$, the additional consumption of Fe and Mn oxides in a range of rapid secondary redox reactions (57), and the rapid exchange of $O_2$ in shallow sediment. $O_2$, $SO_4^{2-}$, and $CH_4$-mediated POC degradation are considered here because the vast majority of marine sedimentary organic carbon is degraded via these pathways (18, 21, 54). Energy may also be provided by fermentation, although its importance is not well known. We use acetate ($CH_3COO^-$) as a proxy for POC. To undertake a study of this scale, such an assumption is essential, since the identity of the thousands of organic molecules consumed by microorganisms in natural settings is spatially disparate and is rarely known. Acetate is a fair choice of molecule, being a regular constituent of marine sediment pore water (26). Fractional changes in Gibbs energies associated with the choice of organic compound do not meaningfully affect our results (see “Uncertainty” section in Supplementary Text). In addition, the Gibbs energies of organic matter degradation are, on a per-electron basis, much more sensitive to the identity of the electron acceptor than that of the organic compound (58, 59). By focusing on the oxidant, we are capturing the first-order energetic differences of organic matter degradation in different environmental settings.

## Gibbs energy calculations

The amount of energy available from the oxidation of organic matter by aerobic, sulfate-reducing, and methanogenic pathways is calculated on the basis of the Gibbs energy function

$$ \Delta G_r = \Delta G_r^0 + RT \ln Q_r \qquad (28) $$

where $\Delta G_r^0$ and $Q_r$ refer to the standard molal Gibbs energy and the reaction quotient of the indicated reaction, respectively, R represents the gas constant, and T denotes temperature in kelvin. Values of $\Delta G_r^0$ were calculated using the revised Helgeson-Kirkham-Flowers equations of state (60–62), the SUPCRT92 software package (63), and thermodynamic data taken from (64–68). Individual values of $Q_r$ are calculated for each reaction using

$$ Q_r = \Pi a_i^{v_i} \qquad (29) $$
and the total steady-state input of POC to that respective layer, $\overline{T}_i$, (grams of C per square centimeter)

$$ \overline{PE}_i = \overline{POC}_i / \overline{T}_i \qquad (24) $$

where

$$ \overline{POC}_i = \int_{z_i}^{z_{i-1}} POC(z) dz \qquad (25) $$

and

$$ \overline{T}_i = POC(z_i - 1) \cdot \Delta z_i \qquad (26) $$

The amount of POC degraded in the ith layer, $\overline{R}_i$, is given by

$$ \overline{R}_i = \int_{z_i}^{z_{i-1}} k(z) \cdot POC(z) dz \qquad (27) $$

Two different global datasets (0.25° × 0.25°) of values for $\omega$, the sedimentation rate, are used: one for Holocene bioturbated and Holocene sediments and the other for Pleistocene sediments.

## Parameters and forcings

The concentration of Holocene and Pleistocene POC at the sediment-water interface, $POC_0$, and sedimentation rates, $\omega$, are constrained according to (37), which used data from (45, 46) and an algorithm that correlates water depth and sedimentation rate (47). Following (36), the bioturbation coefficient is calculated as a function of water depth based on a compilation of empirical data collected in (48). Its values range spatially from 0.59 to 27 cm² year⁻¹, decreasing in magnitude as water depth increases. It is constant throughout the depth of Holocene bioturbated zone and immediately drops to zero beneath it. Values of $\Phi_0$ and $c_0$ are chosen to describe the shelf, margin, and abyss based on sediments that are representative of these domains (see table S3) (49). We partitioned the ocean floor into shelf, margin, and abyss domains to specify values for some of the model parameters (fig. S5 and table S3). We choose a constant v parameter of 0.125 for all three sediment domains, which is characteristic of fresh organic matter and based on observations that the v parameter values do not vary much between sites (29, 41). We link values of the a parameter with sedimentation rates, based on global observations (29, 41). This approach accounts for order-of-magnitude changes in a due to factors that control organic matter transit times from its source to deposition. These parameters thus reflect typically observed RCM parameter variability across various depositional environments. Uncertainties arising from parameter values are addressed in Supplementary Text.

## Global reaction network and catabolic strategy

We designate all sediments shallower than the maximum $O_2$ penetration depth as oxic. To map $O_2$ penetration depth in marine sediments on a global scale (fig. S1), we combine (i) global datasets of oxygen concentrations in marine sediments and (ii) previous modeling results (2). First, we use a global dataset of oxygen concentrations and corresponding water depths (table S4) to formulate a regression model [$R^2$ = 0.64; based on (50)]. This model is extrapolated over the entire ocean floor. For data in columns 81 to 91 (table S4), oxygen penetration depth was determined following the methods described in (51). Second, we superimpose (onto our global map) the results from a modeling study (2), which designated regions of the seafloor where $O_2$ is known to penetrate through all sedimentary

layers to basement rock. We define the sulfate reduction zone as the horizon between the maximum $O_2$ penetration depth and the depth of the sulfate-methane transition. We use data from (52) to map the global sulfate-methane transition depth onto our 0.25° × 0.25° global grid. The reaction zones designated within our model are based on the best available published datasets and include not only global compilations of marine sediment core data (50, 52) but also global models (2, 52). Data from the numerous core samples (n = 1704) documenting the distribution of sulfate-reducing and methanogenic sediments included in our study are available via (52) and the International Ocean Discovery Program. We designate POC degraded in oxic sediments to be mediated by $O_2$, POC degraded in the sulfate-reducing zone to be mediated by $SO_4^{2-}$, and prescribe POC degraded beneath the sulfate-methane transition to methanogenesis (reactions 1 to 3; table S1), thus simplifying the reaction network to account for the most pertinent biogeochemical process, based on widely accepted modeling and experimental studies (18, 21, 53). Although other oxidants are used by microorganisms for degrading organic matter in marine sediments, such as nitrate, and metal oxides, these have been shown to be quantitatively unimportant metabolic pathways on a global scale (29, 44, 54–56)—mainly because of their low concentration compared to $SO_4^{2-}$, the additional consumption of Fe and Mn oxides in a range of rapid secondary redox reactions (57), and the rapid exchange of $O_2$ in shallow sediment. $O_2$, $SO_4^{2-}$, and $CH_4$-mediated POC degradation are considered here because the vast majority of marine sedimentary organic carbon is degraded via these pathways (18, 21, 54). Energy may also be provided by fermentation, although its importance is not well known. We use acetate ($CH_3COO^-$) as a proxy for POC. To undertake a study of this scale, such an assumption is essential, since the identity of the thousands of organic molecules consumed by microorganisms in natural settings is spatially disparate and is rarely known. Acetate is a fair choice of molecule, being a regular constituent of marine sediment pore water (26). Fractional changes in Gibbs energies associated with the choice of organic compound do not meaningfully affect our results (see “Uncertainty” section in Supplementary Text). In addition, the Gibbs energies of organic matter degradation are, on a per-electron basis, much more sensitive to the identity of the electron acceptor than that of the organic compound (58, 59). By focusing on the oxidant, we are capturing the first-order energetic differences of organic matter degradation in different environmental settings.

## Gibbs energy calculations

The amount of energy available from the oxidation of organic matter by aerobic, sulfate-reducing, and methanogenic pathways is calculated on the basis of the Gibbs energy function

$$ \Delta G_r = \Delta G_r^0 + RT \ln Q_r \qquad (28) $$

where $\Delta G_r^0$ and $Q_r$ refer to the standard molal Gibbs energy and the reaction quotient of the indicated reaction, respectively, R represents the gas constant, and T denotes temperature in kelvin. Values of $\Delta G_r^0$ were calculated using the revised Helgeson-Kirkham-Flowers equations of state (60–62), the SUPCRT92 software package (63), and thermodynamic data taken from (64–68). Individual values of $Q_r$ are calculated for each reaction using

$$ Q_r = \Pi a_i^{v_i} \qquad (29) $$
where $a_i$ stands for the activity of the ith species and $v_i$ corresponds to the stoichiometric coefficient of the ith species in the reaction of interest. Activity of the ith species is calculated according to

$$a_i = m_i \gamma_i \quad (30)$$

where $m_i$ and $\gamma_i$ denote the molality and individual activity coefficient of the ith species, respectively. Values of $\gamma_i$ are computed as a function of temperature and ionic strength using an extended version of the Debye-Hückel equation (69).

We have selected specific sets of representative conditions that are used to calculate Gibbs energies of each reaction (table S2). We define reactions at 5°C and 100 bars of pressure. For aerobic heterotrophy, we base our calculations on chemical data from oxic South Pacific Gyre sediments (2). For sulfate-reducing and methanogenic sediments, we use chemical data from anoxic Limfjorden and Peru Margin sediments (9, 70). To address the uncertainty in our results arising from spatial variability of these conditions, we ran additional simulations, where the concentrations of the reactants and products were varied to represent high-energy and low-energy scenarios (see Supplementary Text).

## Cell abundance

Global subseafloor cell abundance was estimated on the basis of a power law formulation from (1)

$$Y = b \cdot z^m \quad (31)$$

where Y (cells per cubic centimeter) is cell abundance, z (meters) denotes the depth below the seafloor, and b and m are parameters based on mean sedimentation rate and distance from landmasses greater than $10^5$ km$^2$. b and m were interpolated onto our 0.25° × 0.25° global grid based on data provided in (1). Cell concentration in the uppermost bioturbated layer ($Y_{bio}$, 0 to 10 cm) was uniformly described, following (1), as

$$Y_{bio} = b \cdot 0.1^m \quad (32)$$

## Global totals

We quantified total subseafloor sedimentary volume, POC degradation, and cell abundance by integrating across specific sediment depth and age horizons and by summing across the global grid representing the entire ocean floor.

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/32/eaba0697/DC1

View/request a protocol for this paper from Bio-protocol.

## REFERENCES AND NOTES

1. J. Kallmeyer, R. Pockalny, R. R. Adhikari, D. C. Smith, S. D'Hondt, Global distribution of microbial abundance and biomass in subseafloor sediment. *Proc. Natl. Acad. Sci. U.S.A.* **109**, 16213–16216 (2012).

2. S. D'Hondt, F. Inagaki, C. A. Zarikian, L. J. Abrams, N. Dubois, T. Engelhardt, H. Evans, T. Ferdelman, B. Gribsholt, R. N. Harris, B. W. Hoppie, J.-H. Hyun, J. Kallmeyer, J. Kim, J. E. Lynch, C. C. McKinley, S. Mitsunobu, Y. Morono, R. W. Murray, R. Pockalny, J. Sauvage, T. Shimono, F. Shiraishi, D. C. Smith, C. E. Smith-Duque, A. J. Spivack, B. O. Steinsbu, Y. Suzuki, M. Szpak, L. Toffin, G. Uramoto, Y. T. Yamaguchi, G.-I. Zhang, X.-H. Zhang, W. Ziebis, Presence of oxygen and aerobic communities from sea floor to basement in deep-sea sediments. *Nat. Geosci.* **8**, 299–304 (2015).

3. B. B. Jørgensen, I. P. G. Marshall, Slow microbial life in the seabed. *Ann. Rev. Mar. Sci.* **8**, 311–332 (2016).

4. K. G. Lloyd, L. Schreiber, D. G. Petersen, K. U. Kjeldsen, M. A. Lever, A. D. Steen, R. Stepanauskas, M. Richter, S. Kleindienst, S. Lenk, A. Schramm, B. B. Jørgensen, Predominant archaea in marine sediments degrade detrital proteins. *Nature* **496**, 215–218 (2013).

5. J. S. Lipp, Y. Morono, F. Inagaki, K.-U. Hinrichs, Significant contribution of Archaea to extant biomass in marine subsurface sediments. *Nature* **454**, 991–994 (2008).

6. H. Roy, J. Kallmeyer, R. R. Adhikari, R. Pockalny, B. B. Jørgensen, S. D'Hondt, Aerobic microbial respiration in 86-million-year-old deep-sea red clay. *Science* **336**, 922–925 (2012).

7. R. J. Parkes, B. Cragg, E. Roussel, G. Webster, A. Weightman, H. Sass, A review of prokaryotic populations and processes in sub-seafloor sediments, including biosphere:geosphere interactions. *Mar. Geol.* **352**, 409–425 (2014).

8. D. H. Rothman, Atmospheric carbon dioxide levels for the last 500 million years. *Proc. Natl. Acad. Sci. U.S.A.* **99**, 4167–4171 (2002).

9. S. D'Hondt, B. B. Jørgensen, D. J. Miller, A. Batzke, R. Blake, B. A. Cragg, H. Cypionka, G. R. Dickens, T. Ferdelman, K.-U. Hinrichs, N. G. Holm, R. Mitterer, A. Spivack, G. Wang, B. Bekins, B. Engelen, K. Ford, G. Gettemy, S. D. Rutherford, H. Sass, C. G. Skilbeck, I. W. Aiello, G. Guérin, C. H. House, F. Inagaki, P. Meister, T. Naehr, S. Niitsuma, R. J. Parkes, A. Schippers, D. C. Smith, A. Teske, J. Wiegel, C. N. Padilla, J. L. S. Acosta, Distributions of microbial activities in deep subseafloor sediments. *Science* **306**, 2216–2221 (2004).

10. B. B. Jørgensen, S. D'Hondt, A starving majority deep beneath the seafloor. *Science* **314**, 932–934 (2006).

11. E. Trembath-Reichert, Y. Morono, A. Ijiri, T. Hoshino, K. S. Dawson, F. Inagaki, V. J. Orphan, Methyl-compound use and slow growth characteristic microbial life in 2-km-deep subseafloor coal and shale beds. *Proc. Natl. Acad. Sci. U.S.A.* **114**, E9206–E9215 (2017).

12. F. Inagaki, K.-U. Hinrichs, Y. Kubo, M. W. Bowles, V. B. Heuer, W.-L. Hong, T. Hoshino, A. Ijiri, H. Imachi, M. Ito, M. Kaneko, M. A. Lever, Y.-S. Lin, B. A. Methé, S. Morita, Y. Morono, W. Tanikawa, M. Bihan, S. A. Bowden, M. Elvert, C. Glombitza, D. Gross, G. J. Harrington, T. Hori, K. Li, D. Limmer, C.-H. Liu, M. Murayama, N. Ohkouchi, S. Ono, Y.-S. Park, S. C. Phillips, X. Prieto-Mollar, M. Purkey, N. Riedinger, Y. Sanada, J. Sauvage, G. Snyder, R. Susilawati, Y. Takano, E. Tasumi, T. Terada, H. Tomaru, E. Trembath-Reichert, D. T. Wang, Y. Yamada, Exploring deep microbial life in coal-bearing sediment down to ~2.5 km below the ocean floor. *Science* **349**, 420–424 (2015).

13. D. E. LaRowe, J. P. Amend, Power limits for microbial life. *Front. Microbiol.* **6**, 718 (2015).

14. Y. Morono, T. Terada, M. Nishizawa, M. Ito, F. Hillion, N. Takahata, Y. Sano, F. Inagaki, Carbon and nitrogen assimilation in deep subseafloor microbial cells. *Proc. Natl. Acad. Sci. U.S.A.* **108**, 18295–18300 (2011).

15. P. Starnawski, T. Bataillon, T. J. G. Ettema, L. M. Jochum, L. Schreiber, X. Chen, M. A. Lever, M. F. Polz, B. B. Jørgensen, A. Schramm, K. U. Kjeldsen, Microbial community assembly and evolution in subseafloor sediment. *Proc. Natl. Acad. Sci. U.S.A.* **114**, 2940–2945 (2017).

16. W. D. Orsi, V. P. Edgcomb, G. D. Christman, J. F. Biddle, Gene expression in the deep biosphere. *Nature* **499**, 205–208 (2013).

17. E. R. Estes, R. Pockalny, S. D'Hondt, F. Inagaki, Y. Morono, R. W. Murray, D. Nordlund, A. J. Spivack, S. D. Wankel, N. Xiao, C. M. Hansel, Persistent organic matter in oxic subseafloor sediment. *Nat. Geosci.* **12**, 126–131 (2019).

18. D. Canfield, E. Kristensen, B. Thamdrup, *Advances in Marine Biology: Aquatic Geomicrobiology* (Academic Press, 2005).

19. B. B. Jørgensen, Mineralization of organic matter in the sea bed—The role of sulphate reduction. *Nature* **296**, 643–645 (1982).

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

21. B. B. Jørgensen, S. Kasten, Sulfur cycling and methane oxidation, in *Marine Geochemistry*, H. D. Schulz, M. Zabel, Eds. (Springer Berlin Heidelberg, 2006), pp. 271–309.

22. K. H. Nealson, P. G. Conrad, Life: Past, present and future. *Philos. Trans. R. Soc. Lond. B Biol. Sci.* **354**, 1923–1939 (1999).

23. C. B. Field, M. J. Behrenfeld, J. T. Randerson, P. Falkowski, Primary production of the biosphere: Integrating terrestrial and oceanic components. *Science* **281**, 237–240 (1998).

24. E. Marschall, M. Jogler, U. Henßge, J. Overmann, Large-scale distribution and activity patterns of an extremely low-light-adapted population of green sulfur bacteria in the Black Sea. *Environ. Microbiol.* **12**, 1348–1362 (2010).

25. L. Tijhuis, M. C. Van Loosdrecht, J. J. Heijnen, A thermodynamically based correlation for maintenance Gibbs energy requirements in aerobic and anaerobic chemotrophic growth. *Biotechnol. Bioeng.* **42**, 509–519 (1993).

26. F. Beulig, H. Roy, C. Glombitza, B. B. Jørgensen, Control on rate and pathway of anaerobic organic carbon degradation in the seabed. *Proc. Natl. Acad. Sci. U.S.A.* **115**, 367–372 (2018).
where $a_i$ stands for the activity of the ith species and $v_i$ corresponds to the stoichiometric coefficient of the ith species in the reaction of interest. Activity of the ith species is calculated according to

$$a_i = m_i \gamma_i \quad (30)$$

where $m_i$ and $\gamma_i$ denote the molality and individual activity coefficient of the ith species, respectively. Values of $\gamma_i$ are computed as a function of temperature and ionic strength using an extended version of the Debye-Hückel equation (69).

We have selected specific sets of representative conditions that are used to calculate Gibbs energies of each reaction (table S2). We define reactions at 5°C and 100 bars of pressure. For aerobic heterotrophy, we base our calculations on chemical data from oxic South Pacific Gyre sediments (2). For sulfate-reducing and methanogenic sediments, we use chemical data from anoxic Limfjorden and Peru Margin sediments (9, 70). To address the uncertainty in our results arising from spatial variability of these conditions, we ran additional simulations, where the concentrations of the reactants and products were varied to represent high-energy and low-energy scenarios (see Supplementary Text).

## Cell abundance

Global subseafloor cell abundance was estimated on the basis of a power law formulation from (1)

$$Y = b \cdot z^m \quad (31)$$

where Y (cells per cubic centimeter) is cell abundance, z (meters) denotes the depth below the seafloor, and b and m are parameters based on mean sedimentation rate and distance from landmasses greater than $10^5$ km$^2$. b and m were interpolated onto our 0.25° × 0.25° global grid based on data provided in (1). Cell concentration in the uppermost bioturbated layer ($Y_{bio}$, 0 to 10 cm) was uniformly described, following (1), as

$$Y_{bio} = b \cdot 0.1^m \quad (32)$$

## Global totals

We quantified total subseafloor sedimentary volume, POC degradation, and cell abundance by integrating across specific sediment depth and age horizons and by summing across the global grid representing the entire ocean floor.

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/32/eaba0697/DC1

View/request a protocol for this paper from Bio-protocol.

## REFERENCES AND NOTES

1. J. Kallmeyer, R. Pockalny, R. R. Adhikari, D. C. Smith, S. D'Hondt, Global distribution of microbial abundance and biomass in subseafloor sediment. *Proc. Natl. Acad. Sci. U.S.A.* **109**, 16213–16216 (2012).

2. S. D'Hondt, F. Inagaki, C. A. Zarikian, L. J. Abrams, N. Dubois, T. Engelhardt, H. Evans, T. Ferdelman, B. Gribsholt, R. N. Harris, B. W. Hoppie, J.-H. Hyun, J. Kallmeyer, J. Kim, J. E. Lynch, C. C. McKinley, S. Mitsunobu, Y. Morono, R. W. Murray, R. Pockalny, J. Sauvage, T. Shimono, F. Shiraishi, D. C. Smith, C. E. Smith-Duque, A. J. Spivack, B. O. Steinsbu, Y. Suzuki, M. Szpak, L. Toffin, G. Uramoto, Y. T. Yamaguchi, G.-I. Zhang, X.-H. Zhang, W. Ziebis, Presence of oxygen and aerobic communities from sea floor to basement in deep-sea sediments. *Nat. Geosci.* **8**, 299–304 (2015).

3. B. B. Jørgensen, I. P. G. Marshall, Slow microbial life in the seabed. *Ann. Rev. Mar. Sci.* **8**, 311–332 (2016).

4. K. G. Lloyd, L. Schreiber, D. G. Petersen, K. U. Kjeldsen, M. A. Lever, A. D. Steen, R. Stepanauskas, M. Richter, S. Kleindienst, S. Lenk, A. Schramm, B. B. Jørgensen, Predominant archaea in marine sediments degrade detrital proteins. *Nature* **496**, 215–218 (2013).

5. J. S. Lipp, Y. Morono, F. Inagaki, K.-U. Hinrichs, Significant contribution of Archaea to extant biomass in marine subsurface sediments. *Nature* **454**, 991–994 (2008).

6. H. Roy, J. Kallmeyer, R. R. Adhikari, R. Pockalny, B. B. Jørgensen, S. D'Hondt, Aerobic microbial respiration in 86-million-year-old deep-sea red clay. *Science* **336**, 922–925 (2012).

7. R. J. Parkes, B. Cragg, E. Roussel, G. Webster, A. Weightman, H. Sass, A review of prokaryotic populations and processes in sub-seafloor sediments, including biosphere:geosphere interactions. *Mar. Geol.* **352**, 409–425 (2014).

8. D. H. Rothman, Atmospheric carbon dioxide levels for the last 500 million years. *Proc. Natl. Acad. Sci. U.S.A.* **99**, 4167–4171 (2002).

9. S. D'Hondt, B. B. Jørgensen, D. J. Miller, A. Batzke, R. Blake, B. A. Cragg, H. Cypionka, G. R. Dickens, T. Ferdelman, K.-U. Hinrichs, N. G. Holm, R. Mitterer, A. Spivack, G. Wang, B. Bekins, B. Engelen, K. Ford, G. Gettemy, S. D. Rutherford, H. Sass, C. G. Skilbeck, I. W. Aiello, G. Guérin, C. H. House, F. Inagaki, P. Meister, T. Naehr, S. Niitsuma, R. J. Parkes, A. Schippers, D. C. Smith, A. Teske, J. Wiegel, C. N. Padilla, J. L. S. Acosta, Distributions of microbial activities in deep subseafloor sediments. *Science* **306**, 2216–2221 (2004).

10. B. B. Jørgensen, S. D'Hondt, A starving majority deep beneath the seafloor. *Science* **314**, 932–934 (2006).

11. E. Trembath-Reichert, Y. Morono, A. Ijiri, T. Hoshino, K. S. Dawson, F. Inagaki, V. J. Orphan, Methyl-compound use and slow growth characteristic microbial life in 2-km-deep subseafloor coal and shale beds. *Proc. Natl. Acad. Sci. U.S.A.* **114**, E9206–E9215 (2017).

12. F. Inagaki, K.-U. Hinrichs, Y. Kubo, M. W. Bowles, V. B. Heuer, W.-L. Hong, T. Hoshino, A. Ijiri, H. Imachi, M. Ito, M. Kaneko, M. A. Lever, Y.-S. Lin, B. A. Methé, S. Morita, Y. Morono, W. Tanikawa, M. Bihan, S. A. Bowden, M. Elvert, C. Glombitza, D. Gross, G. J. Harrington, T. Hori, K. Li, D. Limmer, C.-H. Liu, M. Murayama, N. Ohkouchi, S. Ono, Y.-S. Park, S. C. Phillips, X. Prieto-Mollar, M. Purkey, N. Riedinger, Y. Sanada, J. Sauvage, G. Snyder, R. Susilawati, Y. Takano, E. Tasumi, T. Terada, H. Tomaru, E. Trembath-Reichert, D. T. Wang, Y. Yamada, Exploring deep microbial life in coal-bearing sediment down to ~2.5 km below the ocean floor. *Science* **349**, 420–424 (2015).

13. D. E. LaRowe, J. P. Amend, Power limits for microbial life. *Front. Microbiol.* **6**, 718 (2015).

14. Y. Morono, T. Terada, M. Nishizawa, M. Ito, F. Hillion, N. Takahata, Y. Sano, F. Inagaki, Carbon and nitrogen assimilation in deep subseafloor microbial cells. *Proc. Natl. Acad. Sci. U.S.A.* **108**, 18295–18300 (2011).

15. P. Starnawski, T. Bataillon, T. J. G. Ettema, L. M. Jochum, L. Schreiber, X. Chen, M. A. Lever, M. F. Polz, B. B. Jørgensen, A. Schramm, K. U. Kjeldsen, Microbial community assembly and evolution in subseafloor sediment. *Proc. Natl. Acad. Sci. U.S.A.* **114**, 2940–2945 (2017).

16. W. D. Orsi, V. P. Edgcomb, G. D. Christman, J. F. Biddle, Gene expression in the deep biosphere. *Nature* **499**, 205–208 (2013).

17. E. R. Estes, R. Pockalny, S. D'Hondt, F. Inagaki, Y. Morono, R. W. Murray, D. Nordlund, A. J. Spivack, S. D. Wankel, N. Xiao, C. M. Hansel, Persistent organic matter in oxic subseafloor sediment. *Nat. Geosci.* **12**, 126–131 (2019).

18. D. Canfield, E. Kristensen, B. Thamdrup, *Advances in Marine Biology: Aquatic Geomicrobiology* (Academic Press, 2005).

19. B. B. Jørgensen, Mineralization of organic matter in the sea bed—The role of sulphate reduction. *Nature* **296**, 643–645 (1982).

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

21. B. B. Jørgensen, S. Kasten, Sulfur cycling and methane oxidation, in *Marine Geochemistry*, H. D. Schulz, M. Zabel, Eds. (Springer Berlin Heidelberg, 2006), pp. 271–309.

22. K. H. Nealson, P. G. Conrad, Life: Past, present and future. *Philos. Trans. R. Soc. Lond. B Biol. Sci.* **354**, 1923–1939 (1999).

23. C. B. Field, M. J. Behrenfeld, J. T. Randerson, P. Falkowski, Primary production of the biosphere: Integrating terrestrial and oceanic components. *Science* **281**, 237–240 (1998).

24. E. Marschall, M. Jogler, U. Henßge, J. Overmann, Large-scale distribution and activity patterns of an extremely low-light-adapted population of green sulfur bacteria in the Black Sea. *Environ. Microbiol.* **12**, 1348–1362 (2010).

25. L. Tijhuis, M. C. Van Loosdrecht, J. J. Heijnen, A thermodynamically based correlation for maintenance of Gibbs energy requirements in aerobic and anaerobic chemotrophic growth. *Biotechnol. Bioeng.* **42**, 509–519 (1993).

26. F. Beulig, H. Roy, C. Glombitza, B. B. Jørgensen, Control on rate and pathway of anaerobic organic carbon degradation in the seabed. *Proc. Natl. Acad. Sci. U.S.A.* **115**, 367–372 (2018).
27. T. M. McCollom, J. P. Amend, A thermodynamic assessment of energy requirements for biomass synthesis by chemolithoautotrophic micro-organisms in oxic and anoxic environments. *Geobiology* **3**, 135–144 (2005).

28. M. A. Lever, K. L. Rogers, K. G. Lloyd, J. Overmann, B. Schink, R. K. Thauer, T. M. Hoehler, B. B. Jørgensen, Life under extreme energy limitation: A synthesis of laboratory- and field-based investigations. *FEMS Microbiol. Rev.* **39**, 688–728 (2015).

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

30. J. A. Bradley, J. P. Amend, D. E. La Rowe, Necromass as a limited source of energy for microorganisms in marine sediments. *J. Geophys. Res. Biogeosci.* **123**, 577–590 (2018).

31. S. Braun, S. S. Mhatre, M. Jaussi, H. Roy, K. U. Kjeldsen, C. Pearce, M.-S. Seidenkrantz, B. B. Jørgensen, B. A. Lomstein, Microbial turnover times in the deep seabed studied by amino acid racemization modelling. *Sci. Rep.* **7**, 5680 (2017).

32. Y. M. Bar-On, R. Phillips, R. Milo, The biomass distribution on Earth. *Proc. Natl. Acad. Sci. U.S.A.* **115**, 6506–6511 (2018).

33. C. Magnabosco, L.-H. Lin, H. Dong, M. Bomberg, W. Ghiorse, H. Stan-Lotter, K. Pedersen, T. L. Kieft, E. van Heerden, T. C. Onstott, The biomass and biodiversity of the continental subsurface. *Nat. Geosci.* **11**, 707–717 (2018).

34. S. E. Jones, J. T. Lennon, Dormancy contributes to the maintenance of microbial diversity. *Proc. Natl. Acad. Sci. U.S.A.* **107**, 5881–5886 (2010).

35. J. T. Lennon, S. E. Jones, Microbial seed banks: The ecological and evolutionary implications of dormancy. *Nat. Rev. Microbiol.* **9**, 119–130 (2011).

36. D. E. La Rowe, S. Arndt, J. A. Bradley, E. Burwicz, A. W. Dale, J. P. Amend, Organic carbon and microbial activity in marine sediments on a global scale throughout the Quaternary. *Geochimica et Cosmochimica Acta* **10**.1016/j.gca.2020.07.017 (2020).

37. K. Wallmann, E. Pinero, E. Burwicz, M. Haeckel, C. Hensen, A. Dale, L. Ruepke, The global inventory of methane hydrate in marine sediments: A theoretical approach. *Energies* **5**, 2449–2498 (2012).

38. R. A. Berner, Early diagenesis: A theoretical approach, in *Princeton Series in Geochemistry* (Princeton Univ. Press, 1980).

39. B. P. Boudreau, *Diagenetic Models and Their Implementation: Modelling Transport and Reactions in Aquatic Sediments* (Springer, 1997), vol. 171.

40. L. F. Athy, Density, porosity, and compaction of sedimentary rocks. *Am. Assoc. Pet. Geol. Bull.* **14**, 1–24 (1930).

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

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

43. T. C. Ho, R. Aris, On apparent second-order kinetics. *AIChE J.* **33**, 1050–1051 (1987).

44. A. W. Dale, L. Nickelsen, F. Scholz, C. Hensen, A. Oschlies, K. Wallmann, A revised global estimate of dissolved iron fluxes from marine sediments. *Global Biogeochem. Cycles* **29**, 691–707 (2015).

45. K. Seiter, C. Hensen, J. Schröter, M. Zabel, Organic carbon content in surface sediments—Defining regional provinces. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **51**, 2001–2026 (2004).

46. E. A. Romankevich, A. A. Vetrov, V. I. Peresypkin, Organic matter of the World Ocean. *Russ. Geol. Geophys.* **50**, 299–307 (2009).

47. E. B. Burwicz, L. H. Rüpke, K. Wallmann, Estimation of the global amount of submarine gas hydrates formed via microbial methane formation based on numerical reaction-transport modeling and a novel parameterization of Holocene sedimentation. *Geochim. Cosmochim. Acta* **75**, 4562–4576 (2011).

48. J. J. Middelburg, K. Soetaert, P. M. J. Herman, Empirical relationships for use in global diagenetic models. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **44**, 327–344 (1997).

49. T. Hantschel, A. I. Kauerauf, *Fundamentals of Basin and Petroleum Systems Modeling* (Springer-Verlag, 2009).

50. R. N. Glud, Oxygen dynamics of marine sediments. *Mar. Biol. Res.* **4**, 243–289 (2008).

51. R. N. Glud, J. K. Gundersen, O. Holby, Benthic in situ respiration in the upwelling area off central Chile. *Mar. Ecol. Prog. Ser.* **186**, 9–18 (1999).

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

53. D. Hülse, S. Arndt, S. Daines, P. Regnier, A. Ridgwell, 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).

54. M. Thullner, A. W. Dale, P. Regnier, Global-scale quantification of mineralization pathways in marine sediments: A reaction-transport modeling approach. *Geochem. Geophys. Geosyst.* **10**, Q10012 (2009).

55. R. J. Parkes, G. Webster, B. A. Cragg, A. J. Weightman, C. J. Newberry, T. G. Ferdelman, J. Kallmeyer, B. B. Jørgensen, I. W. Aiello, J. C. Fry, Deep sub-seafloor prokaryotes stimulated at interfaces over geological time. *Nature* **436**, 390–394 (2005).

56. J. J. Middelburg, L. A. Levin, Coastal hypoxia and sediment biogeochemistry. *Biogeosciences* **6**, 1273–1293 (2009).

57. P. Van Cappellen, Y. Wang, 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).

58. D. E. La Rowe, P. Van Cappellen, Degradation of natural organic matter: A thermodynamic analysis. *Geochim. Cosmochim. Acta* **75**, 2030–2042 (2011).

59. D. E. La Rowe, J. P. Amend, Catabolic rates, population sizes and doubling/replacement times of microorganisms in natural settings. *Am. J. Sci.* **315**, 167–203 (2015).

60. H. C. Helgeson, D. H. Kirkham, G. C. Flowers, Theoretical prediction of the thermodynamic behavior of aqueous electrolytes by high pressures and temperatures: IV. Calculation of activity coefficients, osmotic coefficients, and apparent molal and standard and relative partial molal properties to 600°C and 5 kb. *Am. J. Sci.* **281**, 1249–1516 (1981).

61. J. C. Tanger IV, H. C. Helgeson, Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Revised equations of state for the standard partial molal properties of ions and electrolytes. *Am. J. Sci.* **288**, 19–98 (1988).

62. E. L. Shock, E. H. Oelkers, J. W. Johnson, D. A. Sverjensky, H. C. Helgeson, Calculation of the thermodynamic properties of aqueous species at high pressures and temperatures. Effective electrostatic radii, dissociation constants and standard partial molal properties to 1000°C and 5 kbar. *J. Chem. Soc. Faraday Trans.* **88**, 803–826 (1992).

63. J. W. Johnson, E. H. Oelkers, H. C. Helgeson, SUPCRT92: A software package for calculating the standard molal thermodynamic properties of minerals, gases, aqueous species, and reactions from 1 to 5000 bar and 0 to 1000°C. *Comput. Geosci.* **18**, 899–947 (1992).

64. E. L. Shock, H. C. Helgeson, Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Correlation algorithms for ionic species and equation of state predictions to 5 kb and 1000°C. *Geochim. Cosmochim. Acta* **52**, 2009–2036 (1988).

65. E. L. Shock, H. C. Helgeson, D. A. Sverjensky, Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Standard partial molal properties of inorganic neutral species. *Geochim. Cosmochim. Acta* **53**, 2157–2183 (1989).

66. E. L. Shock, H. C. Helgeson, Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Standard partial molal properties of organic species. *Geochim. Cosmochim. Acta* **54**, 915–945 (1990).

67. D. A. Sverjensky, E. L. Shock, H. C. Helgeson, Prediction of the thermodynamic properties of aqueous metal complexes to 1000°C and 5 kb. *Geochim. Cosmochim. Acta* **61**, 1359–1412 (1997).

68. M. D. Schulte, E. L. Shock, R. H. Wood, The temperature dependence of the standard-state thermodynamic properties of aqueous nonelectrolytes. *Geochim. Cosmochim. Acta* **65**, 3919–3930 (2001).

69. H. C. Helgeson, Thermodynamics of hydrothermal systems at elevated temperatures and pressures. *Am. J. Sci.* **267**, 729–804 (1969).

70. B. B. Jørgensen, R. J. Parkes, Role of sulfate reduction and methane production by organic carbon degradation in eutrophic fjord sediments (Limfjorden, Denmark). *Limnol. Oceanogr.* **55**, 1338–1352 (2010).

71. A. Vion, L. Menot, Continental margins between 140m and 3500m depth (2009).

72. A. D. Steen, A. Crits-Christoph, P. Carini, K. M. De Angelis, N. Fierer, K. G. Lloyd, J. C. Thrash, High proportions of bacteria and archaea across most biomes remain uncultured. *ISME J.* **13**, 3126–3130 (2019).

73. L. A. Zinke, C. Glombitza, J. T. Bird, H. Roy, B. B. Jørgensen, K. G. Lloyd, J. P. Amend, B. K. Reese, Microbial organic matter degradation potential in Baltic Sea sediments influenced by depositional conditions and in situ geochemistry. *Appl. Environ. Microbiol.* **85**, e02164-18 (2019).

74. W. D. Orsi, T. A. Richards, W. R. Francis, Predicted microbial secretomes and their target substrates in marine sediment. *Nat. Microbiol.* **3**, 32–37 (2018).

75. B. J. Baker, C. S. Lazar, A. P. Teske, G. J. Dick, Genomic resolution of linkages in carbon, nitrogen, and sulfur cycling among widespread estuary sediment bacteria. *Microbiome* **3**, 14 (2015).

76. I. P. G Marshall, S. M. Karst, P. H. Nielsen, B. B. Jørgensen, Metagenomes from deep Baltic Sea sediments reveal how past and present environmental conditions determine microbial community composition. *Mar. Genomics* **3**, 58–68 (2018).

77. L. A. Zinke, M. M. Mullis, J. T. Bird, I. P. G. Marshall, B. B. Jørgensen, K. G. Lloyd, J. P. Amend, B. K. Reese, Thriving or surviving? Evaluating active microbial guilds in Baltic Sea sediment. *Environ. Microbiol. Rep.* **9**, 528–536 (2017).

78. W. D. Orsi, Ecology and evolution of seafloor and subseafloor microbial communities. *Nat. Rev. Microbiol.* **16**, 671–683 (2018).

79. J. F. Diddle, J. S. Lipp, M. A. Lever, K. G. Lloyd, K. B. Sørensen, R. Anderson, H. F. Fredricks, M. Elvert, T. J. Kelly, D. P. Schrag, M. L. Sogin, J. E. Brenchley, A. Teske, C. H. House, K.-U. Hinrichs, Heterotrophic Archaea dominate sedimentary subsurface ecosystems off Peru. *Proc. Natl. Acad. Sci. U.S.A.* **103**, 3846–3851 (2006).

80. S. E. Ruff, J. F. Diddle, A. P. Teske, K. Knittel, A. Boetius, A. Ramette, Global dispersion and local diversification of the methane seep microbiome. *Proc. Natl. Acad. Sci. U.S.A.* **112**, 4015–4020 (2015).
27. T. M. McCollom, J. P. Amend, A thermodynamic assessment of energy requirements for biomass synthesis by chemolithoautotrophic micro-organisms in oxic and anoxic environments. *Geobiology* **3**, 135–144 (2005).

28. M. A. Lever, K. L. Rogers, K. G. Lloyd, J. Overmann, B. Schink, R. K. Thauer, T. M. Hoehler, B. B. Jørgensen, Life under extreme energy limitation: A synthesis of laboratory- and field-based investigations. *FEMS Microbiol. Rev.* **39**, 688–728 (2015).

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

30. J. A. Bradley, J. P. Amend, D. E. La Rowe, Necromass as a limited source of energy for microorganisms in marine sediments. *J. Geophys. Res. Biogeosci.* **123**, 577–590 (2018).

31. S. Braun, S. S. Mhatre, M. Jaussi, H. Røy, K. U. Kjeldsen, C. Pearce, M.-S. Seidenkrantz, B. B. Jørgensen, B. A. Lomstein, Microbial turnover times in the deep seabed studied by amino acid racemization modelling. *Sci. Rep.* **7**, 5680 (2017).

32. Y. M. Bar-On, R. Phillips, R. Milo, The biomass distribution on Earth. *Proc. Natl. Acad. Sci. U.S.A.* **115**, 6506–6511 (2018).

33. C. Magnabosco, L.-H. Lin, H. Dong, M. Bomberg, W. Ghiorse, H. Stan-Lotter, K. Pedersen, T. L. Kieft, E. van Heerden, T. C. Onstott, The biomass and biodiversity of the continental subsurface. *Nat. Geosci.* **11**, 707–717 (2018).

34. S. E. Jones, J. T. Lennon, Dormancy contributes to the maintenance of microbial diversity. *Proc. Natl. Acad. Sci. U.S.A.* **107**, 5881–5886 (2010).

35. J. T. Lennon, S. E. Jones, Microbial seed banks: The ecological and evolutionary implications of dormancy. *Nat. Rev. Microbiol.* **9**, 119–130 (2011).

36. D. E. La Rowe, S. Arndt, J. A. Bradley, E. Burwicz, A. W. Dale, J. P. Amend, Organic carbon and microbial activity in marine sediments on a global scale throughout the Quaternary. *Geochimica et Cosmochimica Acta* **10**.1016/j.gca.2020.07.017 (2020).

37. K. Wallmann, E. Pinero, E. Burwicz, M. Haeckel, C. Hensen, A. Dale, L. Ruepke, The global inventory of methane hydrate in marine sediments: A theoretical approach. *Energies* **5**, 2449–2498 (2012).

38. R. A. Berner, Early diagenesis: A theoretical approach, in *Princeton Series in Geochemistry* (Princeton Univ. Press, 1980).

39. B. P. Boudreau, *Diagenetic Models and Their Implementation: Modelling Transport and Reactions in Aquatic Sediments* (Springer, 1997), vol. 171.

40. L. F. Athy, Density, porosity, and compaction of sedimentary rocks. *Am. Assoc. Pet. Geol. Bull.* **14**, 1–24 (1930).

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

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

43. T. C. Ho, R. Aris, On apparent second-order kinetics. *AIChE J.* **33**, 1050–1051 (1987).

44. A. W. Dale, L. Nickelsen, F. Scholz, C. Hensen, A. Oschlies, K. Wallmann, A revised global estimate of dissolved iron fluxes from marine sediments. *Global Biogeochem. Cycles* **29**, 691–707 (2015).

45. K. Seiter, C. Hensen, J. Schröter, M. Zabel, Organic carbon content in surface sediments—Defining regional provinces. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **51**, 2001–2026 (2004).

46. E. A. Romankevich, A. A. Vetrov, V. I. Peresypkin, Organic matter of the World Ocean. *Russ. Geol. Geophys.* **50**, 299–307 (2009).

47. E. B. Burwicz, L. H. Rüpke, K. Wallmann, Estimation of the global amount of submarine gas hydrates formed via microbial methane formation based on numerical reaction-transport modeling and a novel parameterization of Holocene sedimentation. *Geochim. Cosmochim. Acta* **75**, 4562–4576 (2011).

48. J. J. Middelburg, K. Soetaert, P. M. J. Herman, Empirical relationships for use in global diagenetic models. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **44**, 327–344 (1997).

49. T. Hantschel, A. I. Kauerauf, *Fundamentals of Basin and Petroleum Systems Modeling* (Springer-Verlag, 2009).

50. R. N. Glud, Oxygen dynamics of marine sediments. *Mar. Biol. Res.* **4**, 243–289 (2008).

51. R. N. Glud, J. K. Gundersen, O. Holby, Benthic in situ respiration in the upwelling area off central Chile. *Mar. Ecol. Prog. Ser.* **186**, 9–18 (1999).

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

53. D. Hülse, S. Arndt, S. Daines, P. Regnier, A. Ridgwell, 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).

54. M. Thullner, A. W. Dale, P. Regnier, Global-scale quantification of mineralization pathways in marine sediments: A reaction-transport modeling approach. *Geochem. Geophys. Geosyst.* **10**, Q10012 (2009).

55. R. J. Parkes, G. Webster, B. A. Cragg, A. J. Weightman, C. J. Newberry, T. G. Ferdelman, J. Kallmeyer, B. B. Jørgensen, I. W. Aiello, J. C. Fry, Deep sub-seafloor prokaryotes stimulated at interfaces over geological time. *Nature* **436**, 390–394 (2005).

56. J. J. Middelburg, L. A. Levin, Coastal hypoxia and sediment biogeochemistry. *Biogeosciences* **6**, 1273–1293 (2009).

57. P. Van Cappellen, Y. Wang, 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).

58. D. E. La Rowe, P. Van Cappellen, Degradation of natural organic matter: A thermodynamic analysis. *Geochim. Cosmochim. Acta* **75**, 2030–2042 (2011).

59. D. E. La Rowe, J. P. Amend, Catabolic rates, population sizes and doubling/replacement times of microorganisms in natural settings. *Am. J. Sci.* **315**, 167–203 (2015).

60. H. C. Helgeson, D. H. Kirkham, G. C. Flowers, Theoretical prediction of the thermodynamic behavior of aqueous electrolytes by high pressures and temperatures: IV. Calculation of activity coefficients, osmotic coefficients, and apparent molal and standard and relative partial molal properties to 600°C and 5 kb. *Am. J. Sci.* **281**, 1249–1516 (1981).

61. J. C. Tanger IV, H. C. Helgeson, Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Revised equations of state for the standard partial molal properties of ions and electrolytes. *Am. J. Sci.* **288**, 19–98 (1988).

62. E. L. Shock, E. H. Oelkers, J. W. Johnson, D. A. Sverjensky, H. C. Helgeson, Calculation of the thermodynamic properties of aqueous species at high pressures and temperatures. Effective electrostatic radii, dissociation constants and standard partial molal properties to 1000°C and 5 kbar. *J. Chem. Soc. Faraday Trans.* **88**, 803–826 (1992).

63. J. W. Johnson, E. H. Oelkers, H. C. Helgeson, SUPCRT92: A software package for calculating the standard molal thermodynamic properties of minerals, gases, aqueous species, and reactions from 1 to 5000 bar and 0 to 1000°C. *Comput. Geosci.* **18**, 899–947 (1992).

64. E. L. Shock, H. C. Helgeson, Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Correlation algorithms for ionic species and equation of state predictions to 5 kb and 1000°C. *Geochim. Cosmochim. Acta* **52**, 2009–2036 (1988).

65. E. L. Shock, H. C. Helgeson, D. A. Sverjensky, Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Standard partial molal properties of inorganic neutral species. *Geochim. Cosmochim. Acta* **53**, 2157–2183 (1989).

66. E. L. Shock, H. C. Helgeson, Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Standard partial molal properties of organic species. *Geochim. Cosmochim. Acta* **54**, 915–945 (1990).

67. D. A. Sverjensky, E. L. Shock, H. C. Helgeson, Prediction of the thermodynamic properties of aqueous metal complexes to 1000°C and 5 kb. *Geochim. Cosmochim. Acta* **61**, 1359–1412 (1997).

68. M. D. Schulte, E. L. Shock, R. H. Wood, The temperature dependence of the standard-state thermodynamic properties of aqueous nonelectrolytes. *Geochim. Cosmochim. Acta* **65**, 3919–3930 (2001).

69. H. C. Helgeson, Thermodynamics of hydrothermal systems at elevated temperatures and pressures. *Am. J. Sci.* **267**, 729–804 (1969).

70. B. B. Jørgensen, R. J. Parkes, Role of sulfate reduction and methane production by organic carbon degradation in eutrophic fjord sediments (Limfjorden, Denmark). *Limnol. Oceanogr.* **55**, 1338–1352 (2010).

71. A. Vion, L. Menot, Continental margins between 140m and 3500m depth (2009).

72. A. D. Steen, A. Crits-Christoph, P. Carini, K. M. De Angelis, N. Fierer, K. G. Lloyd, J. C. Thrash, High proportions of bacteria and archaea across most biomes remain uncultured. *ISME J.* **13**, 3126–3130 (2019).

73. L. A. Zinke, C. Glombitza, J. T. Bird, H. Røy, B. B. Jørgensen, K. G. Lloyd, J. P. Amend, B. K. Reese, Microbial organic matter degradation potential in Baltic Sea sediments influenced by depositional conditions and in situ geochemistry. *Appl. Environ. Microbiol.* **85**, e02164-18 (2019).

74. W. D. Orsi, T. A. Richards, W. R. Francis, Predicted microbial secretomes and their target substrates in marine sediment. *Nat. Microbiol.* **3**, 32–37 (2018).

75. B. J. Baker, C. S. Lazar, A. P. Teske, G. J. Dick, Genomic resolution of linkages in carbon, nitrogen, and sulfur cycling among widespread estuary sediment bacteria. *Microbiome* **3**, 14 (2015).

76. I. P. G Marshall, S. M. Karst, P. H. Nielsen, B. B. Jørgensen, Metagenomes from deep Baltic Sea sediments reveal how past and present environmental conditions determine microbial community composition. *Mar. Genomics* **3**, 58–68 (2018).

77. L. A. Zinke, M. M. Mullis, J. T. Bird, I. P. G. Marshall, B. B. Jørgensen, K. G. Lloyd, J. P. Amend, B. K. Reese, Thriving or surviving? Evaluating active microbial guilds in Baltic Sea sediment. *Environ. Microbiol. Rep.* **9**, 528–536 (2017).

78. W. D. Orsi, Ecology and evolution of seafloor and subseafloor microbial communities. *Nat. Rev. Microbiol.* **16**, 671–683 (2018).

79. J. F. Diddle, J. S. Lipp, M. A. Lever, K. G. Lloyd, K. B. Sørensen, R. Anderson, H. F. Fredricks, M. Elvert, T. J. Kelly, D. P. Schrag, M. L. Sogin, J. E. Brenchley, A. Teske, C. H. House, K.-U. Hinrichs, Heterotrophic Archaea dominate sedimentary subsurface ecosystems off Peru. *Proc. Natl. Acad. Sci. U.S.A.* **103**, 3846–3851 (2006).

80. S. E. Ruff, J. F. Diddle, A. P. Teske, K. Knittel, A. Boetius, A. Ramette, Global dispersion and local diversification of the methane seep microbiome. *Proc. Natl. Acad. Sci. U.S.A.* **112**, 4015–4020 (2015).
81. B. K. Reese, L. A. Zinke, M. S. Sobol, D. E. La Rowe, B. N. Orcutt, X. Zhang, U. Jaekel, F. Wang, T. Dittmar, D. Defforey, B. Tully, A. Paytan, J. B. Sylvan, J. P. Amend, K. J. Edwards, P. Girguis, Nitrogen cycling of active bacteria within oligotrophic sediment of the mid- atlantic ridge flank. *Geomicrobiol. J.* **35**, 468–483 (2018).

82. B. J. Tully, J. F. Heidelberg, Potential mechanisms for microbial energy acquisition in oxic deep-sea sediments. *Appl. Environ. Microbiol.* **82**, 4232–4243 (2016).

83. A. Durbin, A. Teske, Microbial diversity and stratification of South Pacific abyssal marine sediments. *Environ. Microbiol.* **13**, 3219–3234 (2011).

84. J. Chun, A. Oren, A. Ventosa, H. Christensen, D. R. Arahal, M. S. da Costa, A. P. Rooney, H. Yi, X.-W. Xu, S. De Meyer, M. E. Trujillo, Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes. *Int. J. Syst. Evol. Microbiol.* **68**, 461–466 (2018).

85. A. Tkacz, M. Hortala, P. S. Poole, Absolute quantitation of microbiota abundance in environmental samples. *Microbiome* **6**, 110 (2018).

86. C. L. Myers, D. R. Barrett, M. A. Hibbs, C. Huttenhower, O. G. Troyanskaya, Finding function: Evaluation methods for functional genomic data. *BMC Genomics* **7**, 187 (2006).

87. J. M. Janda, S. L. Abbott, 16S rRNA gene sequencing for bacterial identification in the diagnostic laboratory: Pluses, perils, and pitfalls. *J. Clin. Microbiol.* **45**, 2761–2764 (2007).

88. K. H. Nealson, Sediment bacteria: Who's there, what are they doing, and what's new? *Annu. Rev. Earth Planet. Sci.* **25**, 403–434 (1997).

89. C. Arnosti, C. Bell, D. L. Moorhead, R. L. Sinsabaugh, A. D. Steen, M. Stromberger, M. Wallenstein, M. N. Weintraub, Extracellular enzymes in terrestrial, freshwater, and marine environments: Perspectives on system variability and common research needs. *Biogeochemistry* **117**, 5–21 (2014).

90. K. G. Lloyd, A. D. Steen, J. Ladau, J. Yin, L. Crosby, Phylogenetically novel uncultured microbial cells dominate Earth microbiomes. *mSystems* **3**, e00055-18 (2018).

91. D. J. Burdige, K. G. Gardner, Molecular weight distribution of dissolved organic carbon in marine sediment pore waters. *Mar. Chem.* **62**, 45–64 (1998).

92. C. A. Fox, H. A. Abdulla, D. J. Burdige, J. P. Lewicki, T. Komada, Composition of dissolved organic matter in pore waters of anoxic marine sediments analyzed by ¹H nuclear magnetic resonance spectroscopy. *Front. Mar. Sci.* **5**, 172 (2018).

93. D. E. LaRowe, E. Burwicz, S. Arndt, A. W. Dale, J. P. Amend, Temperature and volume of global marine sediments. *Geology* **45**, 275–278 (2017).

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

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

96. J. J. Middelburg, A simple rate model for organic matter decomposition in marine sediments. *Geochim. Cosmochim. Acta* **53**, 1577–1581 (1989).

97. J. T. Westrich, R. A. Berner, The role of sedimentary organic-matter in bacterial sulfate reduction: The G Model Tested. *Limnol. Oceanogr.* **29**, 236–249 (1984).

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

99. R. N. Glud, J. K. Gundersen, B. B. Jørgensen, N. P. Revsbech, H. D. Schulz, Diffusive and total oxygen uptake of deep-sea sediments in the eastern South Atlantic Ocean: In situ and laboratory measurements. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **41**, 1767–1788 (1994).

100. R. N. Glud, J. K. Gundersen, H. Roy, B. B. Jørgensen, Seasonal dynamics of benthic O₂ uptake in a semienclosed bay: Importance of diffusion and faunal activity. *Limnol. Oceanogr.* **48**, 1265–1276 (2003).

101. R. N. Glud, O. Holby, F. Hoffmann, D. Canfield, Benthic mineralization and exchange in Arctic sediments (Svalbard, Norway). *Mar. Ecol. Prog. Ser.* **173**, 237–251 (1998).

102. F. Wenzhöfer, O. Holby, O. Kohls, Deep penetrating benthic oxygen profiles measured in situ by oxygen optodes. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **48**, 1741–1755 (2001).

103. F. Wenzhöfer, M. Adler, O. Kohls, C. Hensen, B. Strottmann, S. Boehme, H. D. Schulz, Calcite dissolution driven by benthic mineralization in the deep-sea: In situ measurements of Ca²⁺, pH, pCO₂ and O₂. *Geochim. Cosmochim. Acta* **65**, 2677–2690 (2001).

104. R. A. Jahnke, S. R. Emerson, C. E. Reimers, J. Schuffert, K. Ruttenberg, D. Archer, Benthic recycling of biogenic debris in the eastern tropical Atlantic Ocean. *Geochim. Cosmochim. Acta* **53**, 2947–2960 (1989).

105. F. Wenzhöfer, R. N. Glud, Benthic carbon mineralization in the Atlantic: A synthesis based on in situ data from the last decade. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **49**, 1255–1279 (2002).

106. U. Witte, N. Aberle, M. Sand, F. Wenzhöfer, Rapid response of a deep-sea benthic community to POM enrichment: An in situ experimental study. *Mar. Ecol. Prog. Ser.* **251**, 27–36 (2003).

107. U. Witte, F. Wenzhöfer, S. Sommer, A. Boetius, P. Heinz, N. Aberle, M. Sand, A. Cremer, W.-R. Abraham, B. B. Jørgensen, O. Pfannkuche, In situ experimental evidence of the fate of a phytodetritus pulse at the abyssal sea floor. *Nature* **424**, 763–766 (2003).

108. K. S. Black, G. R. Fones, O. C. Peppe, H. A. Kennedy, I. Bentaleb, An autonomous benthic lander: Preliminary observations from the UK BENBO thematic programme. *Cont. Shelf Res.* **21**, 859–877 (2001).

109. R. N. Glud, H. Stahl, P. Berg, F. Wenzhöfer, K. Oguri, H. Kitazato, In situ microscale variation in distribution and consumption of 2: A case study from a deep ocean margin sediment (Sagami Bay, Japan). *Limnol. Oceanogr.* **54**, 1–12 (2009).

110. P. Giordani, W. Helder, E. Koning, S. Misercocchi, R. Danovaro, A. Malaguti, Gradients of benthic–pelagic coupling and carbon budgets in the Adriatic and Northern Ionian Sea. *J. Mar. Syst.* **33–34**, 365–387 (2002).

111. B. Lansard, C. Rabouille, L. Denis, C. Grenz, In situ oxygen uptake rates by coastal sediments under the influence of the Rhône River (NW Mediterranean Sea). *Cont. Shelf Res.* **28**, 1501–1510 (2008).

112. B. Lansard, C. Rabouille, L. Denis, C. Grenz, Benthic remineralization at the land-ocean interface: A case study of the Rhône River (NW Mediterranean Sea). *Estuar. Coast. Shell Sci.* **81**, 544–554 (2009).

113. O. Sachs, E. J. Sauter, M. Schlüter, M. M. R. van der Loeff, K. Jerosch, O. Holby, Benthic organic carbon flux and oxygen penetration reflect different plankton provinces in the Southern Ocean. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **56**, 1319–1335 (2009).

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

**Acknowledgments:** We thank R. Glud (University of Southern Denmark) and F. Wenzhöfer (Alfred Wegener Institute, Bremerhaven) for providing oxygen penetration data and for discussion related to the approach and D. Hüelse for discussion. We also thank J. Kallmeyer (German Research Centre for Geosciences, GFZ) and R. Pockalny (University of Rhode Island) for sharing published datasets on cell abundance and oxygen in global sediments. We also thank the anonymous reviewers of the manuscript, who provided many helpful comments and suggestions. **Funding:** This research was supported by the NSF-sponsored Center for Dark Energy Biosphere Investigations (C-DEBI, OCE-0939564), the Deep Carbon Observatory (the Deep Life Modelling and Visualization (DLMV) and Deep Energy programs), the NASA Astrobiology Institute–Life Underground (NAI-LU, NNA13AA92A), the NASA Astrobiology Program under the Joint NASA-NSF Ideas Lab on the “Origins of Life” (NSF Solicitation 16-570), the Alexander von Humboldt Foundation, a NERC grant (NE/T010967/1), the European Union Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement no. 643052 (C-CASCADES), and the USC Zumberge Fund Individual Grant Program. This is C-DEBI contribution 534. **Author contributions:** J.A.B. designed the study with contributions from D.E.L. and J.P.A. D.E.L., J.A.B., S.A., E.B., and A.W.D. designed and implemented the RTM. M.E. assisted with defining the sulfate-methane transition depth. J.A.B. performed additional calculations, analyzed data, and wrote the manuscript with significant contributions from all co-authors. **Competing interests:** The authors declare that they have no competing interests. **Data and materials availability:** All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper are available on https://www.bco-dmo.org/project/776336.

Submitted 1 November 2019

Accepted 23 June 2020

Published 5 August 2020

10.1126/sciadv.aba0697

**Citation:** J. A. Bradley, S. Arndt, J. P. Amend, E. Burwicz, A. W. Dale, M. Egger, D. E. LaRowe, Wide-spread energy limitation to life in global subseafloor sediments. *Sci. Adv.* **6**, eaba0697 (2020).
81. B. K. Reese, L. A. Zinke, M. S. Sobol, D. E. La Rowe, B. N. Orcutt, X. Zhang, U. Jaekel, F. Wang, T. Dittmar, D. Defforey, B. Tully, A. Paytan, J. B. Sylvan, J. P. Amend, K. J. Edwards, P. Girguis, Nitrogen cycling of active bacteria within oligotrophic sediment of the mid- atlantic ridge flank. *Geomicrobiol. J.* **35**, 468–483 (2018).

82. B. J. Tully, J. F. Heidelberg, Potential mechanisms for microbial energy acquisition in oxic deep-sea sediments. *Appl. Environ. Microbiol.* **82**, 4232–4243 (2016).

83. A. Durbin, A. Teske, Microbial diversity and stratification of South Pacific abyssal marine sediments. *Environ. Microbiol.* **13**, 3219–3234 (2011).

84. J. Chun, A. Oren, A. Ventosa, H. Christensen, D. R. Arahal, M. S. da Costa, A. P. Rooney, H. Yi, X.-W. Xu, S. De Meyer, M. E. Trujillo, Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes. *Int. J. Syst. Evol. Microbiol.* **68**, 461–466 (2018).

85. A. Tkacz, M. Hortala, P. S. Poole, Absolute quantitation of microbiota abundance in environmental samples. *Microbiome* **6**, 110 (2018).

86. C. L. Myers, D. R. Barrett, M. A. Hibbs, C. Huttenhower, O. G. Troyanskaya, Finding function: Evaluation methods for functional genomic data. *BMC Genomics* **7**, 187 (2006).

87. J. M. Janda, S. L. Abbott, 16S rRNA gene sequencing for bacterial identification in the diagnostic laboratory: Pluses, perils, and pitfalls. *J. Clin. Microbiol.* **45**, 2761–2764 (2007).

88. K. H. Nealson, Sediment bacteria: Who's there, what are they doing, and what's new? *Annu. Rev. Earth Planet. Sci.* **25**, 403–434 (1997).

89. C. Arnosti, C. Bell, D. L. Moorhead, R. L. Sinsabaugh, A. D. Steen, M. Stromberger, M. Wallenstein, M. N. Weintraub, Extracellular enzymes in terrestrial, freshwater, and marine environments: Perspectives on system variability and common research needs. *Biogeochemistry* **117**, 5–21 (2014).

90. K. G. Lloyd, A. D. Steen, J. Ladau, J. Yin, L. Crosby, Phylogenetically novel uncultured microbial cells dominate Earth microbiomes. *mSystems* **3**, e00055-18 (2018).

91. D. J. Burdige, K. G. Gardner, Molecular weight distribution of dissolved organic carbon in marine sediment pore waters. *Mar. Chem.* **62**, 45–64 (1998).

92. C. A. Fox, H. A. Abdulla, D. J. Burdige, J. P. Lewicki, T. Komada, Composition of dissolved organic matter in pore waters of anoxic marine sediments analyzed by ¹H nuclear magnetic resonance spectroscopy. *Front. Mar. Sci.* **5**, 172 (2018).

93. D. E. LaRowe, E. Burwicz, S. Arndt, A. W. Dale, J. P. Amend, Temperature and volume of global marine sediments. *Geology* **45**, 275–278 (2017).

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

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

96. J. J. Middelburg, A simple rate model for organic matter decomposition in marine sediments. *Geochim. Cosmochim. Acta* **53**, 1577–1581 (1989).

97. J. T. Westrich, R. A. Berner, The role of sedimentary organic-matter in bacterial sulfate reduction: The G Model Tested. *Limnol. Oceanogr.* **29**, 236–249 (1984).

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

99. R. N. Glud, J. K. Gundersen, B. B. Jørgensen, N. P. Revsbech, H. D. Schulz, Diffusive and total oxygen uptake of deep-sea sediments in the eastern South Atlantic Ocean: In situ and laboratory measurements. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **41**, 1767–1788 (1994).

100. R. N. Glud, J. K. Gundersen, H. Roy, B. B. Jørgensen, Seasonal dynamics of benthic O₂ uptake in a semienclosed bay: Importance of diffusion and faunal activity. *Limnol. Oceanogr.* **48**, 1265–1276 (2003).

101. R. N. Glud, O. Holby, F. Hoffmann, D. Canfield, Benthic mineralization and exchange in Arctic sediments (Svalbard, Norway). *Mar. Ecol. Prog. Ser.* **173**, 237–251 (1998).

102. F. Wenzhöfer, O. Holby, O. Kohls, Deep penetrating benthic oxygen profiles measured in situ by oxygen optodes. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **48**, 1741–1755 (2001).

103. F. Wenzhöfer, M. Adler, O. Kohls, C. Hensen, B. Strottmann, S. Boehme, H. D. Schulz, Calcite dissolution driven by benthic mineralization in the deep-sea: In situ measurements of Ca²⁺, pH, pCO₂ and O₂. *Geochim. Cosmochim. Acta* **65**, 2677–2690 (2001).

104. R. A. Jahnke, S. R. Emerson, C. E. Reimers, J. Schuffert, K. Ruttenberg, D. Archer, Benthic recycling of biogenic debris in the eastern tropical Atlantic Ocean. *Geochim. Cosmochim. Acta* **53**, 2947–2960 (1989).

105. F. Wenzhöfer, R. N. Glud, Benthic carbon mineralization in the Atlantic: A synthesis based on in situ data from the last decade. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **49**, 1255–1279 (2002).

106. U. Witte, N. Aberle, M. Sand, F. Wenzhöfer, Rapid response of a deep-sea benthic community to POM enrichment: An in situ experimental study. *Mar. Ecol. Prog. Ser.* **251**, 27–36 (2003).

107. U. Witte, F. Wenzhöfer, S. Sommer, A. Boetius, P. Heinz, N. Aberle, M. Sand, A. Cremer, W.-R. Abraham, B. B. Jørgensen, O. Pfannkuche, In situ experimental evidence of the fate of a phytodetritus pulse at the abyssal sea floor. *Nature* **424**, 763–766 (2003).

108. K. S. Black, G. R. Fones, O. C. Peppe, H. A. Kennedy, I. Bentaleb, An autonomous benthic lander: Preliminary observations from the UK BENBO thematic programme. *Cont. Shelf Res.* **21**, 859–877 (2001).

109. R. N. Glud, H. Stahl, P. Berg, F. Wenzhöfer, K. Oguri, H. Kitazato, In situ microscale variation in distribution and consumption of 2: A case study from a deep ocean margin sediment (Sagami Bay, Japan). *Limnol. Oceanogr.* **54**, 1–12 (2009).

110. P. Giordani, W. Helder, E. Koning, S. Misercocchi, R. Danovaro, A. Malaguti, Gradients of benthic-pelagic coupling and carbon budgets in the Adriatic and Northern Ionian Sea. *J. Mar. Syst.* **33–34**, 365–387 (2002).

111. B. Lansard, C. Rabouille, L. Denis, C. Grenz, In situ oxygen uptake rates by coastal sediments under the influence of the Rhône River (NW Mediterranean Sea). *Cont. Shelf Res.* **28**, 1501–1510 (2008).

112. B. Lansard, C. Rabouille, L. Denis, C. Grenz, Benthic remineralization at the land-ocean interface: A case study of the Rhône River (NW Mediterranean Sea). *Estuar. Coast. Shell Sci.* **81**, 544–554 (2009).

113. O. Sachs, E. J. Sauter, M. Schlüter, M. M. R. van der Loeff, K. Jerosch, O. Holby, Benthic organic carbon flux and oxygen penetration reflect different plankton provinces in the Southern Ocean. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **56**, 1319–1335 (2009).

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

**Acknowledgments:** We thank R. Glud (University of Southern Denmark) and F. Wenzhöfer (Alfred Wegener Institute, Bremerhaven) for providing oxygen penetration data and for discussion related to the approach and D. Hüelse for discussion. We also thank J. Kallmeyer (German Research Centre for Geosciences, GFZ) and R. Pockalny (University of Rhode Island) for sharing published datasets on cell abundance and oxygen in global sediments. We also thank the anonymous reviewers of the manuscript, who provided many helpful comments and suggestions. **Funding:** This research was supported by the NSF-sponsored Center for Dark Energy Biosphere Investigations (C-DEBI, OCE-0939564), the Deep Carbon Observatory (the Deep Life Modelling and Visualization (DLMV) and Deep Energy programs), the NASA Astrobiology Institute–Life Underground (NAI-LU, NNA13AA92A), the NASA Astrobiology Program under the Joint NASA-NSF Ideas Lab on the “Origins of Life” (NSF Solicitation 16-570), the Alexander von Humboldt Foundation, a NERC grant (NE/T010967/1), the European Union Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement no. 643052 (C-CASCADES), and the USC Zumberge Fund Individual Grant Program. This is C-DEBI contribution 534. **Author contributions:** J.A.B. designed the study with contributions from D.E.L. and J.P.A. D.E.L., J.A.B., S.A., E.B., and A.W.D. designed and implemented the RTM. M.E. assisted with defining the sulfate-methane transition depth. J.A.B. performed additional calculations, analyzed data, and wrote the manuscript with significant contributions from all co-authors. **Competing interests:** The authors declare that they have no competing interests. **Data and materials availability:** All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper are available on https://www.bco-dmo.org/project/776336.

Submitted 1 November 2019

Accepted 23 June 2020

Published 5 August 2020

10.1126/sciadv.aba0697

**Citation:** J. A. Bradley, S. Arndt, J. P. Amend, E. Burwicz, A. W. Dale, M. Egger, D. E. LaRowe, Wide-spread energy limitation to life in global subseafloor sediments. *Sci. Adv.* **6**, eaba0697 (2020).