iBet uBet web content aggregator. Adding the entire web to your favor.
iBet uBet web content aggregator. Adding the entire web to your favor.



Link to original content: https://doi.org/10.1038/nature23646
Very large release of mostly volcanic carbon during the Palaeocene–Eocene Thermal Maximum | Nature
Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Letter
  • Published:

Very large release of mostly volcanic carbon during the Palaeocene–Eocene Thermal Maximum

Abstract

The Palaeocene–Eocene Thermal Maximum1,2 (PETM) was a global warming event that occurred about 56 million years ago, and is commonly thought to have been driven primarily by the destabilization of carbon from surface sedimentary reservoirs such as methane hydrates3. However, it remains controversial whether such reservoirs were indeed the source of the carbon that drove the warming1,3,4,5. Resolving this issue is key to understanding the proximal cause of the warming, and to quantifying the roles of triggers versus feedbacks. Here we present boron isotope data—a proxy for seawater pH—that show that the ocean surface pH was persistently low during the PETM. We combine our pH data with a paired carbon isotope record in an Earth system model in order to reconstruct the unfolding carbon-cycle dynamics during the event6,7. We find strong evidence for a much larger (more than 10,000 petagrams)—and, on average, isotopically heavier—carbon source than considered previously8,9. This leads us to identify volcanism associated with the North Atlantic Igneous Province10,11, rather than carbon from a surface reservoir, as the main driver of the PETM. This finding implies that climate-driven amplification of organic carbon feedbacks probably played only a minor part in driving the event. However, we find that enhanced burial of organic matter seems to have been important in eventually sequestering the released carbon and accelerating the recovery of the Earth system12.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Figure 1: Stable-isotope data from DSDP Site 401.
Figure 2: δ13C evolution and pH reconstruction based on analysis of boron isotopes in M. subbotinae from Site 401.
Figure 3: Assimilation of data from our Earth system model.

Similar content being viewed by others

References

  1. McInerney, F. A. & Wing, S. L. The Paleocene–Eocene Thermal Maximum: a perturbation of carbon cycle, climate, and biosphere with implications for the future. Annu. Rev. Earth Planet. Sci. 39, 489–516 (2011)

    Article  ADS  CAS  Google Scholar 

  2. Dunkley Jones, T. et al. Climate model and proxy data constraints on ocean warming across the Paleocene–Eocene Thermal Maximum. Earth Sci. Rev. 125, 123–145 (2013)

    Article  ADS  CAS  Google Scholar 

  3. Dickens, G. R., O’Neil, J. R., Rea, D. K. & Owen, R. M. Dissociation of oceanic methane hydrate as a cause of the carbon isotope excursion at the end of the Paleocene. Paleoceanography 10, 965–971 (1995)

    Article  ADS  Google Scholar 

  4. DeConto, R. M. et al. Past extreme warming events linked to massive carbon release from thawing permafrost. Nature 484, 87–91 (2012)

    Article  ADS  CAS  PubMed  Google Scholar 

  5. Higgins, J. A. & Schrag, D. P. Beyond methane: towards a theory for the Paleocene–Eocene Thermal Maximum. Earth Planet. Sci. Lett. 245, 523–537 (2006)

    Article  ADS  CAS  Google Scholar 

  6. Ridgwell, A. & Schmidt, D. N. Past constraints on the vulnerability of marine calcifiers to massive carbon dioxide release. Nat. Geosci. 3, 196–200 (2010)

    Article  ADS  CAS  Google Scholar 

  7. Cui, Y. et al. Slow release of fossil carbon during the Palaeocene–Eocene Thermal Maximum. Nat. Geosci. 4, 481–485 (2011)

    Article  ADS  CAS  Google Scholar 

  8. Zeebe, R. E., Zachos, J. C. & Dickens, G. R. Carbon dioxide forcing alone insufficient to explain Palaeocene–Eocene Thermal Maximum warming. Nat. Geosci. 2, 576–580 (2009)

    Article  ADS  CAS  Google Scholar 

  9. Frieling, J. et al. Thermogenic methane release as a cause for the long duration of the PETM. Proc. Natl Acad. Sci. USA 113, 12059–12064 (2016)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  10. Svensen, H. et al. Release of methane from a volcanic basin as a mechanism for initial Eocene global warming. Nature 429, 542–545 (2004)

    Article  ADS  CAS  PubMed  Google Scholar 

  11. Storey, M., Duncan, R. A. & Swisher, C. C. Paleocene–Eocene Thermal Maximum and the opening of the Northeast Atlantic. Science 316, 587–589 (2007)

    Article  ADS  CAS  PubMed  Google Scholar 

  12. Bowen, G. J. & Zachos, J. C. Rapid carbon sequestration at the termination of the Palaeocene–Eocene Thermal Maximum. Nat. Geosci. 3, 866–869 (2010)

    Article  ADS  CAS  Google Scholar 

  13. Rohling, E. J. et al. Making sense of palaeoclimate sensitivity. Nature 491, 683–691 (2012)

    Article  ADS  CAS  Google Scholar 

  14. Gibbs, S. J. et al. Ocean warming, not acidification, controlled coccolithophore response during past greenhouse climate change. Geology 44, 59–62 (2016)

    Article  ADS  Google Scholar 

  15. Hönisch, B. et al. The geological record of ocean acidification. Science 335, 1058–1063 (2012)

    Article  ADS  PubMed  CAS  Google Scholar 

  16. Charles, A. J. et al. Constraints on the numerical age of the Paleocene-Eocene boundary. Geochem. Geophys. Geosyst. 12, Q0AA17 (2011)

    Article  CAS  Google Scholar 

  17. Le Quéré, C. et al. Global carbon budget 2016. Earth Syst. Sci. Data 8, 605–649 (2016)

    Article  ADS  Google Scholar 

  18. Schaller, M. F., Fung, M. K., Wright, J. D., Katz, M. E. & Kent, D. V. Impact ejecta at the Paleocene–Eocene boundary. Science 354, 225–229 (2016)

    Article  ADS  CAS  PubMed  Google Scholar 

  19. Wieczorek, R., Fantle, M. S., Kump, L. R. & Ravizza, G. Geochemical evidence for volcanic activity prior to and enhanced terrestrial weathering during the Paleocene Eocene Thermal Maximum. Geochim. Cosmochim. Acta 119, 391–410 (2013)

    Article  ADS  CAS  Google Scholar 

  20. Penman, D. E., Hönisch, B., Zeebe, R. E., Thomas, E. & Zachos, J. C. Rapid and sustained surface ocean acidification during the Paleocene–Eocene Thermal Maximum. Paleoceanography 29, 357–369 (2014)

    Article  ADS  Google Scholar 

  21. Anagnostou, E. et al. Changing atmospheric CO2 concentration was the primary driver of early Cenozoic climate. Nature 533, 380–384 (2016)

    Article  ADS  CAS  PubMed  Google Scholar 

  22. Schubert, B. A. & Jahren, A. H. Reconciliation of marine and terrestrial carbon isotope excursions based on changing atmospheric CO2 levels. Nat. Commun. 4, 1653 (2013)

    Article  ADS  PubMed  CAS  Google Scholar 

  23. Penman, D. E. et al. An abyssal carbonate compensation depth overshoot in the aftermath of the Palaeocene–Eocene Thermal Maximum. Nat. Geosci. 9, 575–580 (2016)

    Article  ADS  CAS  Google Scholar 

  24. Turner, S. K. & Ridgwell, A. Development of a novel empirical framework for interpreting geological carbon isotope excursions, with implications for the rate of carbon injection across the PETM. Earth Planet. Sci. Lett. 435, 1–13 (2016)

    Article  ADS  CAS  Google Scholar 

  25. Röhl, U., Westerhold, T., Bralower, T. J. & Zachos, J. C. On the duration of the Paleocene–Eocene thermal maximum (PETM). Geochem. Geophys. Geosyst. 8, Q12002 (2007)

    Article  ADS  CAS  Google Scholar 

  26. Svensen, H., Planke, S. & Corfu, F. Zircon dating ties NE Atlantic sill emplacement to initial Eocene global warming. J. Geol. Soc. Lond. 167, 433–436 (2010)

    Article  CAS  Google Scholar 

  27. Saunders, A. D. Two LIPs and two Earth-system crises: the impact of the North Atlantic Igneous Province and the Siberian Traps on the Earth-surface carbon cycle. Geol. Mag. 153, 201–222 (2016)

    Article  ADS  CAS  Google Scholar 

  28. Rampino, M. R. Peraluminous igneous rocks as an indicator of thermogenic methane release from the North Atlantic Volcanic Province at the time of the Paleocene–Eocene Thermal Maximum (PETM). Bull. Volcanol. 75, 1–5 (2013)

    Article  Google Scholar 

  29. Ma, Z. et al. Carbon sequestration during the Palaeocene–Eocene Thermal Maximum by an efficient biological pump. Nat. Geosci. 7, 382–388 (2014)

    Article  ADS  CAS  Google Scholar 

  30. Dickson, A. J., Cohen, A. S. & Coe, A. L. Seawater oxygenation during the Paleocene–Eocene Thermal Maximum. Geology 40, 639–642 (2012)

    Article  ADS  CAS  Google Scholar 

  31. Jenkyns, H. C. Cretaceous anoxic events—from continents to oceans. J. Geol. Soc. Lond. 137, 171–188 (1980)

    Article  Google Scholar 

  32. Turner, S. K., Sexton, P. F., Charles, C. D. & Norris, R. D. Persistence of carbon release events through the peak of early Eocene global warmth. Nat. Geosci. 7, 748–751 (2014)

    Article  ADS  CAS  Google Scholar 

  33. Payne, J. L. et al. Calcium isotope constraints on the end-Permian mass extinction. Proc. Natl Acad. Sci. USA 107, 8543–8548 (2010)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  34. Pardo, A., Keller, G., Molina, E. & Canudo, J. I. Planktic foraminiferal turnover across the Paleocene–Eocene transition at DSDP site 401, Bay of Biscay, North Atlantic. Mar. Micropaleontol. 29, 129–158 (1997)

    Article  ADS  Google Scholar 

  35. Bornemann, A. et al. Persistent environmental change after the Paleocene–Eocene Thermal Maximum in the eastern North Atlantic. Earth Planet. Sci. Lett. 394, 70–81 (2014)

    Article  ADS  CAS  Google Scholar 

  36. Sexton, P. F., Wilson, P. A. & Pearson, P. N. Microstructural and geochemical perspectives on planktic foraminiferal preservation: “glassy” versus “frosty”. Geochem. Geophys. Geosyst. 7, Q12P19 (2006)

    Article  CAS  Google Scholar 

  37. Foster, G. L. Seawater pH, pCO2 and [CO32−] variations in the Caribbean Sea over the last 130 kyr: a boron isotope and B/Ca study of planktic forminifera. Earth Planet. Sci. Lett. 271, 254–266 (2008)

    Article  ADS  CAS  Google Scholar 

  38. Foster, G. L. et al. Interlaboratory comparison of boron isotope analyses of boric acid, seawater and marine CaCO3 by MC-ICPMS and NTIMS. Chem. Geol. 358, 1–14 (2013)

    Article  ADS  CAS  Google Scholar 

  39. Barker, S., Greaves, M. & Elderfield, H. A study of cleaning procedures used for foraminiferal Mg/Ca paleothermometry. Geochem. Geophys. Geosyst. 4, 8407 (2003)

    Article  ADS  CAS  Google Scholar 

  40. Henehan, M. J. et al. Calibration of the boron isotope proxy in the planktonic foraminifera Globigerinoides ruber for use in palaeo-CO2 reconstruction. Earth Planet. Sci. Lett. 364, 111–122 (2013)

    Article  ADS  CAS  Google Scholar 

  41. Nunes, F. & Norris, R. D. Abrupt reversal in ocean overturning during the Palaeocene/Eocene warm period. Nature 439, 60–63 (2006)

    Article  ADS  CAS  PubMed  Google Scholar 

  42. Sanyal, A., Bijma, J., Spero, H. & Lea, D. W. Empirical relationship between pH and the boron isotopic composition of Globigerinoides sacculifer: implications for the boron isotope paleo-pH proxy. Paleoceanography 16, 515–519 (2001)

    Article  ADS  Google Scholar 

  43. Martínez-Boti, M. A. et al. Boron isotope evidence for oceanic carbon dioxide leakage during the last deglaciation. Nature 518, 219–222 (2015)

    Article  ADS  PubMed  CAS  Google Scholar 

  44. Zeebe, R. E., Wolf-Gladrow, D. A., Bijma, J. & Hönisch, B. Vital effects in foraminifera do not compromise the use of delta B-11 as a paleo-pH indicator: evidence from modeling. Paleoceanography 18, 1043 (2003)

    Article  ADS  Google Scholar 

  45. Hönisch, B. et al. The influence of symbiont photosynthesis on the boron isotopic composition of foraminifera shells. Mar. Micropaleontol. 49, 87–96 (2003)

    Article  ADS  Google Scholar 

  46. Foster, G. L. & Rae, J. W. B. Reconstructing ocean pH with boron isotopes in foraminifera. Annu. Rev. Earth Planet. Sci. 44, 207–237 (2016)

    Article  ADS  CAS  Google Scholar 

  47. Klochko, K., Kaufman, A. J., Yao, W. S., Byrne, R. H. & Tossell, J. A. Experimental measurement of boron isotope fractionation in seawater. Earth Planet. Sci. Lett. 248, 276–285 (2006)

    Article  ADS  CAS  Google Scholar 

  48. Henehan, M. J. et al. A new boron isotope-pH calibration for Orbulina universa, with implications for understanding and accounting for ‘vital effects’. Earth Planet. Sci. Lett. 454, 282–292 (2016)

    Article  ADS  CAS  Google Scholar 

  49. Kim, S.-T. & O’Neil, J. R. Equilibrium and nonequilibrium oxygen isotope effects in synthetic carbonates. Geochim. Cosmochim. Acta 61, 3461–3475 (1997)

    Article  ADS  CAS  Google Scholar 

  50. Tindall, J. et al. Modelling the oxygen isotope distribution of ancient seawater using a coupled ocean-atmosphere GCM: implications for reconstructing early Eocene climate. Earth Planet. Sci. Lett. 292, 265–273 (2010)

    Article  ADS  CAS  Google Scholar 

  51. Evans, D. & Müller, W. Deep time foraminifera Mg/Ca paleothermometry: nonlinear correction for secular change in seawater Mg/Ca. Paleoceanography 27, PA4205 (2012)

    Article  ADS  Google Scholar 

  52. Spivack, A. J. & Edmond, J. M. Boron isotope exchange between seawater and the oceanic crust. Geochim. Cosmochim. Acta 51, 1033–1043 (1987)

    Article  ADS  CAS  Google Scholar 

  53. Lemarchand, D., Gaillardet, J., Lewin, E. & Allegre, C. J. Boron isotope systematics in large rivers: implications for the marine boron budget and paleo-pH reconstruction over the Cenozoic. Chem. Geol. 190, 123–140 (2002)

    Article  ADS  CAS  Google Scholar 

  54. Paillard, D., Labeyrie, L. & Yiou, P. Macintosh program performs time-series analysis. EOS 77, 379 (1996)

    Article  ADS  Google Scholar 

  55. Giusberti, L. et al. Mode and tempo of the Paleocene–Eocene thermal maximum in an expanded section from the Venetian pre-Alps. Geol. Soc. Am. Bull. 119, 391–412 (2007)

    Article  ADS  CAS  Google Scholar 

  56. Röhl, U., Bralower, T. J., Norris, R. D. & Wefer, G. New chronology for the late Paleocene thermal maximum and its environmental implications. Geology 28, 927–930 (2000)

    Article  ADS  Google Scholar 

  57. Farley, K. A. & Eltgroth, S. F. An alternative age model for the Paleocene–Eocene Thermal Maximum using extraterrestrial He-3. Earth Planet. Sci. Lett. 208, 135–148 (2003)

    Article  ADS  CAS  Google Scholar 

  58. Wright, J. D. & Schaller, M. F. Evidence for a rapid release of carbon at the Paleocene–Eocene thermal maximum. Proc. Natl Acad. Sci. USA 110, 15908–15913 (2013)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  59. Zeebe, R. E., Dickens, G. R., Ridgwell, A., Sluijs, A. & Thomas, E. Onset of carbon isotope excursion at the Paleocene–Eocene thermal maximum took millennia, not 13 years. Proc. Natl Acad. Sci. USA 111, E1062–E1063 (2014)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  60. Pearson, P. N. & Nicholas, C. J. Layering in the Paleocene/Eocene boundary of the Millville core is drilling disturbance. Proc. Natl Acad. Sci. USA 111, E1064–E1065 (2014)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  61. Stassen, P., Speijer, R. P. & Thomas, E. Unsettled puzzle of the Marlboro clays. Proc. Natl Acad. Sci. USA 111, E1066–E1067 (2014)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  62. Wright, J. D. & Schaller, M. F. Reply to Pearson and Nicholas, Stassen et al., and Zeebe et al.: Teasing out the missing piece of the PETM puzzle. Proc. Natl Acad. Sci. USA 111, E1068–E1071 (2014)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  63. Pearson, P. N. & Thomas, E. Drilling disturbance and constraints on the onset of the Paleocene–Eocene boundary carbon isotope excursion in New Jersey. Clim. Past 11, 95–104 (2015)

    Article  Google Scholar 

  64. Zeebe, R. E., Ridgwell, A. & Zachos, J. C. Anthropogenic carbon release rate unprecedented during the past 66 million years. Nat. Geosci. 9, 325–329 (2016)

    Article  ADS  CAS  Google Scholar 

  65. Brady, P. V. The effect of silicate weathering on global temperature and atmospheric CO2 . J. Geophys. Res. 96, 18101–18106 (1991)

    Article  ADS  CAS  Google Scholar 

  66. Edwards, N. R. & Marsh, R. Uncertainties due to transport-parameter sensitivity in an efficient 3-D ocean-climate model. Clim. Dyn. 24, 415–433 (2005)

    Article  Google Scholar 

  67. Ridgwell, A. et al. Marine geochemical data assimilation in an efficient Earth system model of global biogeochemical cycling. Biogeosciences 4, 87–104 (2007)

    Article  ADS  CAS  Google Scholar 

  68. Ridgwell, A. & Hargreaves, J. C. Regulation of atmospheric CO2 by deep-sea sediments in an Earth system model. Global Biogeochem. Cycles 21, GB2008 (2007)

    Article  ADS  CAS  Google Scholar 

  69. Colbourn, G., Ridgwell, A. & Lenton, T. M. The time scale of the silicate weathering negative feedback on atmospheric CO2 . Glob. Biogeochem. Cycles 29, 583–596 (2015)

    Article  ADS  CAS  Google Scholar 

  70. Lord, N. S., Ridgwell, A., Thorne, M. C. & Lunt, D. J. An impulse response function for the ‘long tail’ of excess atmospheric CO2 in an Earth system model. Global Biogeochem. Cycles 30, 2–17 (2016)

    Article  ADS  CAS  Google Scholar 

  71. Cui, Y. & Kump, L. R. Global warming and the end-Permian extinction event: proxy and modeling perspectives. Earth Sci. Rev. 149, 5–22 (2015)

    Article  ADS  CAS  Google Scholar 

  72. IPCC. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (eds Stocker, T. F. et al.) (Cambridge Univ. Press, 2013)

  73. Ridgwell, A. Glacial–Interglacial Perturbations in the Global Carbon Cycle. PhD thesis, Univ. East Anglia (2001)

  74. Cao, L. et al. The role of ocean transport in the uptake of anthropogenic CO2 . Biogeosciences 6, 375–390 (2009)

    Article  ADS  CAS  Google Scholar 

Download references

Acknowledgements

This study was funded by a UK Ocean Acidification Research Program NERC/DEFRA/DECC grant (NE/H017518/1) to P.N.P., G.L.F. and P.F.S. (also supporting M.G.). A.R. was supported by a Heising–Simons Foundation award, and by EU grant ERC 2013-CoG-617313. E.T. was in part supported by the National Science Foundation Division of Ocean Sciences (grant no. NSF OCE 1536611). H.P. was in part supported by ERC grant 2013-CoG-617462. This study used samples provided by the International Ocean Discovery Program. We thank A. Milton at the University of Southampton for maintaining the mass spectrometers used in this study, and M. Davies at The Open University for assistance with sample preparation. We thank L. Haxhiaj and D. Nürnberg at GEOMAR Kiel and H. Kuhnert at MARUM Bremen for their help with carbon and oxygen isotope analyses.

Author information

Authors and Affiliations

Authors

Contributions

G.L.F., P.F.S. and P.N.P. developed the concept and designed the study. M.G. and E.A. carried out the preparation of chemical samples, as well as elemental and isotopic analyses. P.F.S. performed foraminifer taxonomy and prepared foraminifer samples for the analyses. R.D.N. and E.T. supplied washed coarse-fraction samples. P.F.S. developed the age model. A.R. devised and conducted the Earth system modelling and analysis. H.P. carried out the carbon and oxygen isotopic analyses. M.G., A.R., G.L.F. and P.F.S. led the writing of the manuscript. All authors contributed to the interpretation of results and writing of the final text.

Corresponding author

Correspondence to Marcus Gutjahr.

Ethics declarations

Competing interests

The authors declare no competing financial interests.

Additional information

Reviewer Information Nature thanks T. Bralower, K. Meissner and the other anonymous reviewer(s) for their contribution to the peer review of this work.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

Extended Data Figure 1 Elemental and stable-isotope cross-plots for M. subbotinae.

a–f, We screened samples for chemical consistency by checking various elemental ratios (Al/Ca, B/Ca and Mg/Ca) as well as measured δ11B and δ13C.

Extended Data Figure 2 Stable-isotope data for foraminifera and bulk carbonate.

Foraminifera-based stable-isotope compositions were generated from identical samples after splitting the δ13C/δ18O fraction from the δ11B fraction. For both foraminifera and bulk carbonate, the data are plotted against the depth of the core sample. The same data are plotted against age relative to CIE in Fig. 1.

Extended Data Figure 3 Illustration of the results of δ11B-to-pH conversion, and differences in age models.

a, We used our δ11B measurements as a proxy to calculate the evolution of pH at Site 401 during the PETM CIE, using either the borate ion47 (red) or the T. sacculifer43 (green) calibration. The age scale used follows ref. 25. b, Direct comparison of our two age models, plotting the reconstructed pH evolution of Site 401 using either the age model of ref. 57 or our preferred age model25. c, Expanded view of b.

Extended Data Figure 4 Selection of age-model tie points.

Comparison of bulk carbonate δ13C (top) and δ18O (bottom) values for Site 401 and Site 690 (ref. 25). Vertical lines show the age tie points used to derive the age model relative to the PETM CIE (see Methods).

Extended Data Figure 5 Key results of sensitivity experiments.

The figure shows the influence of uncertainties in the CIE onset duration on diagnosed total carbon release. For these idealized experiments, the CIE onset phase was assumed to occur linearly, with the decline in ocean surface dissolved inorganic carbon (DIC) δ13C (by 3.5‰) and pH (by 0.3 pH units) having a duration (Δτonset) that varied from 100 to 20,000 years; thereafter, the target pH and δ13C values were held constant until the end of the experiment (at 50,000 years). ah, The evolution with time of these target ocean surface variables, with pH on the left-hand y axis, and δ13C(DIC) on the right-hand axis. ip, Maximum carbon emission rate per time interval. qx, Cumulative carbon emission per time interval in Eg of carbon. y, z, aaaf, Average emitted δ13C per time interval.

Extended Data Figure 6 Spatial and temporal evolution of mean annual surface ocean pH in the Earth system model cGENIE, shown both for the PETM and for preindustrial and future times.

a, Black line, global mean surface ocean pH values across the PETM, from experiment R07sm_Corg (these are our main pH estimates, obtained using the inorganic borate ion calibration and the RH07 age model, and including an assumption of organic carbon burial after the peak PETM). Red circles represent the annual mean pH values at Site 401 (location shown in b) in the model, taken at times in the model simulation that have corresponding δ11B-derived pH data points (see Fig. 3b; note that we do not use all of the observed data points). b, Model-projected spatial pattern of annual mean surface ocean pH at time zero (that is, PETM onset). The star shows the location of Site 401. cf, Model-projected spatial patterns of the annual mean surface ocean pH anomaly compared with time 0, for the highlighted time points from a (5.0, 31.6, 58.2 and 71.5 kyr after onset). g, Model-projected spatial pattern of annual mean surface ocean pH in the modern ocean under pre-industrial (year 1765) atmospheric CO2 levels (278 p.p.m.). The model is configured as described in ref. 74 and driven with a CO2 emissions scenario that is consistent with RCP 6.0. h, i, Model-projected spatial pattern of the annual mean surface ocean pH anomaly compared with that for 1765, at years 2010 and 2050. The scale is as for cf.

Extended Data Figure 7 Spatial and temporal evolution of surface sedimentary calcium carbonate content in cGENIE during the PETM.

a, Black line, global mean surface levels (in wt%) of sedimentary calcium carbonate (CaCO3) across the PETM, from experiment R07sm_Corg. White circles, times from PETM onset onwards that correspond to the δ11B-derived pH data points in Fig. 3b and Extended Data Fig. 6. The white circles do not represent ‘values’, and simply mark specific time points. b, Model-projected spatial pattern of surface sedimentary wt% CaCO3 at time zero (PETM onset). Shown are the locations of sites for which surface ocean pH has been reconstructed (see Fig. 2) and at which detailed down-core model–data comparison is carried out (Extended Data Fig. 9). cf, Model-projected spatial patterns of the surface sedimentary wt% CaCO3 anomaly compared with time 0, for the time points highlighted in a. g, For reference, the assumed seafloor bathymetry in the model (and the locations of the four data-rich sites that are discussed in Supplementary Information).

Extended Data Figure 8 Spatial and temporal evolution of sea surface temperature in cGENIE during the PETM.

a, Black line, global and annual mean SSTs during the PETM, from experiment R07sm_Corg. Yellow circles, annual mean SST values at Site 401 in the model, at the times from PETM onset onwards that correspond to the δ11B-derived pH data points (see Fig. 3b). Blue and orange circles, δ18O- and Mg/Ca-derived SST estimates, respectively. b, Model-projected spatial pattern of annual mean SSTs at time 0. The star shows the location of Site 401. cf, Model-projected spatial patterns of the annual mean SST anomaly compared with time 0, for the time points shown in a.

Extended Data Figure 9 Down-core model–data evaluation at four data-rich sites.

ap, Comparisons for four Ocean Drilling Program Sites (401, 865, 1,209 and 1,263) for which surface ocean pH has been reconstructed across the PETM (Fig. 2; this study and ref. 20). q, The palaeo-locations of these sites in the cGENIE Earth system model. ap, Model–data comparisons are made for: wt% CaCO3 (a, e, i, m); δ13C values from bulk carbonate (b, f, j, n); and surface ocean pH (c, g, k, o). Panels d, h, l and p provide an orientation in time, showing the projected evolution of atmospheric δ13C from CO2 in the model. For wt% CaCO3 and δ13C of bulk carbonate, model points (resolved at 1-cm intervals) are plotted as filled yellow circles. Model-projected pH values (global and annual means, as in Fig. 3h and Extended Data Fig. 6a) and atmospheric δ13C values for CO2 are shown as continuous red lines. In all cases, observed data values are shown as asterisks. For Sites 865, 1,209 and 1,263, the age models—employing the original relative age-model constraints20 used to convert from model-simulated sediment depths (resolved at 1-cm intervals) at each location in cGENIE—were calculated using a constant detrital flux accumulation rate. The observed data are plotted on the respective Site-690-derived age models25. Both model-based and data-based age scales are synchronized to time 0 (PETM onset; horizontal line). See Supplementary Information for details.

Extended Data Table 1 Key results from individual model runs

Supplementary information

Supplementary Information

This file contains a detailed account about the Earth System Modelling approaches that were used and additional references. (PDF 331 kb)

Supplementary Table 1

This table contains foraminifera-based stable isotope results, relative sample ages, selected elemental ratios as well as the calculated mixed layer pH. (XLSX 59 kb)

Supplementary Table 2

This table contains bulk carbonate stable carbon and oxygen isotope results, presented alongside relative ages following our two alternative age models (see Methods). (XLSX 80 kb)

PowerPoint slides

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gutjahr, M., Ridgwell, A., Sexton, P. et al. Very large release of mostly volcanic carbon during the Palaeocene–Eocene Thermal Maximum. Nature 548, 573–577 (2017). https://doi.org/10.1038/nature23646

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/nature23646

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing