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
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 51 print issues and online access
$199.00 per year
only $3.90 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
References
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)
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)
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)
DeConto, R. M. et al. Past extreme warming events linked to massive carbon release from thawing permafrost. Nature 484, 87–91 (2012)
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)
Ridgwell, A. & Schmidt, D. N. Past constraints on the vulnerability of marine calcifiers to massive carbon dioxide release. Nat. Geosci. 3, 196–200 (2010)
Cui, Y. et al. Slow release of fossil carbon during the Palaeocene–Eocene Thermal Maximum. Nat. Geosci. 4, 481–485 (2011)
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)
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)
Svensen, H. et al. Release of methane from a volcanic basin as a mechanism for initial Eocene global warming. Nature 429, 542–545 (2004)
Storey, M., Duncan, R. A. & Swisher, C. C. Paleocene–Eocene Thermal Maximum and the opening of the Northeast Atlantic. Science 316, 587–589 (2007)
Bowen, G. J. & Zachos, J. C. Rapid carbon sequestration at the termination of the Palaeocene–Eocene Thermal Maximum. Nat. Geosci. 3, 866–869 (2010)
Rohling, E. J. et al. Making sense of palaeoclimate sensitivity. Nature 491, 683–691 (2012)
Gibbs, S. J. et al. Ocean warming, not acidification, controlled coccolithophore response during past greenhouse climate change. Geology 44, 59–62 (2016)
Hönisch, B. et al. The geological record of ocean acidification. Science 335, 1058–1063 (2012)
Charles, A. J. et al. Constraints on the numerical age of the Paleocene-Eocene boundary. Geochem. Geophys. Geosyst. 12, Q0AA17 (2011)
Le Quéré, C. et al. Global carbon budget 2016. Earth Syst. Sci. Data 8, 605–649 (2016)
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)
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)
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)
Anagnostou, E. et al. Changing atmospheric CO2 concentration was the primary driver of early Cenozoic climate. Nature 533, 380–384 (2016)
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)
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)
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)
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)
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)
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)
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)
Ma, Z. et al. Carbon sequestration during the Palaeocene–Eocene Thermal Maximum by an efficient biological pump. Nat. Geosci. 7, 382–388 (2014)
Dickson, A. J., Cohen, A. S. & Coe, A. L. Seawater oxygenation during the Paleocene–Eocene Thermal Maximum. Geology 40, 639–642 (2012)
Jenkyns, H. C. Cretaceous anoxic events—from continents to oceans. J. Geol. Soc. Lond. 137, 171–188 (1980)
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)
Payne, J. L. et al. Calcium isotope constraints on the end-Permian mass extinction. Proc. Natl Acad. Sci. USA 107, 8543–8548 (2010)
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)
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)
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)
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)
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)
Barker, S., Greaves, M. & Elderfield, H. A study of cleaning procedures used for foraminiferal Mg/Ca paleothermometry. Geochem. Geophys. Geosyst. 4, 8407 (2003)
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)
Nunes, F. & Norris, R. D. Abrupt reversal in ocean overturning during the Palaeocene/Eocene warm period. Nature 439, 60–63 (2006)
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)
Martínez-Boti, M. A. et al. Boron isotope evidence for oceanic carbon dioxide leakage during the last deglaciation. Nature 518, 219–222 (2015)
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)
Hönisch, B. et al. The influence of symbiont photosynthesis on the boron isotopic composition of foraminifera shells. Mar. Micropaleontol. 49, 87–96 (2003)
Foster, G. L. & Rae, J. W. B. Reconstructing ocean pH with boron isotopes in foraminifera. Annu. Rev. Earth Planet. Sci. 44, 207–237 (2016)
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)
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)
Kim, S.-T. & O’Neil, J. R. Equilibrium and nonequilibrium oxygen isotope effects in synthetic carbonates. Geochim. Cosmochim. Acta 61, 3461–3475 (1997)
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)
Evans, D. & Müller, W. Deep time foraminifera Mg/Ca paleothermometry: nonlinear correction for secular change in seawater Mg/Ca. Paleoceanography 27, PA4205 (2012)
Spivack, A. J. & Edmond, J. M. Boron isotope exchange between seawater and the oceanic crust. Geochim. Cosmochim. Acta 51, 1033–1043 (1987)
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)
Paillard, D., Labeyrie, L. & Yiou, P. Macintosh program performs time-series analysis. EOS 77, 379 (1996)
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)
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)
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)
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)
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)
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)
Stassen, P., Speijer, R. P. & Thomas, E. Unsettled puzzle of the Marlboro clays. Proc. Natl Acad. Sci. USA 111, E1066–E1067 (2014)
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)
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)
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)
Brady, P. V. The effect of silicate weathering on global temperature and atmospheric CO2 . J. Geophys. Res. 96, 18101–18106 (1991)
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)
Ridgwell, A. et al. Marine geochemical data assimilation in an efficient Earth system model of global biogeochemical cycling. Biogeosciences 4, 87–104 (2007)
Ridgwell, A. & Hargreaves, J. C. Regulation of atmospheric CO2 by deep-sea sediments in an Earth system model. Global Biogeochem. Cycles 21, GB2008 (2007)
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)
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)
Cui, Y. & Kump, L. R. Global warming and the end-Permian extinction event: proxy and modeling perspectives. Earth Sci. Rev. 149, 5–22 (2015)
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)
Ridgwell, A. Glacial–Interglacial Perturbations in the Global Carbon Cycle. PhD thesis, Univ. East Anglia (2001)
Cao, L. et al. The role of ocean transport in the uptake of anthropogenic CO2 . Biogeosciences 6, 375–390 (2009)
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
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
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). a–h, 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. i–p, Maximum carbon emission rate per time interval. q–x, Cumulative carbon emission per time interval in Eg of carbon. y, z, aa–af, 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. c–f, 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 c–f.
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). c–f, 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. c–f, 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.
a–p, 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. a–p, 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.
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)
Rights and permissions
About this article
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
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nature23646