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:

Constraints on global mean sea level during Pliocene warmth

Nature volume 574, pages 233–236 (2019)Cite this article

Subjects

Abstract

Reconstructing the evolution of sea level during past warmer epochs such as the Pliocene provides insight into the response of sea level and ice sheets to prolonged warming1. Although estimates of the global mean sea level (GMSL) during this time do exist, they vary by several tens of metres2,3,4, hindering the assessment of past and future ice-sheet stability. Here we show that during the mid-Piacenzian Warm Period, which was on average two to three degrees Celsius warmer than the pre-industrial period5, the GMSL was about 16.2 metres higher than today owing to global ice-volume changes, and around 17.4 metres when thermal expansion of the oceans is included. During the even warmer Pliocene Climatic Optimum (about four degrees Celsius warmer than pre-industrial levels)6, our results show that the GMSL was 23.5 metres above the present level, with an additional 1.6 metres from thermal expansion. We provide six GMSL data points, ranging from 4.39 to 3.27 million years ago, that are based on phreatic overgrowths on speleothems from the western Mediterranean (Mallorca, Spain). This record is unique owing to its clear relationship to sea level, its reliable U–Pb ages and its long timespan, which allows us to quantify uncertainties on potential uplift. Our data indicate that ice sheets are very sensitive to warming and provide important calibration targets for future ice-sheet models7.

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

Access options

Buy this article

  • Purchase on Springer Link
  • Instant access to full article PDF
Buy now

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

Fig. 1: POS in Coves d’Artà.
Fig. 2: GIA and other corrections.
Fig. 3: Pliocene sea-level and CO2 concentration changes.

Similar content being viewed by others

Data availability

The data produced in this study are available at the NOAA (https://www.ncdc.noaa.gov/paleo/study/27530) and Pangaea (https://doi.pangaea.de/10.1594/PANGAEA.905851) data repositories.

Code availability

The computer code used to do the sea-level (GIA) calculation, written in MATLAB, is available on github: https://github.com/jaustermann/SLcode.

References

  1. Dutton, A. et al. Sea-level rise due to polar ice-sheet mass loss during past warm periods. Science 349, aaa4019 (2015).

    CAS  Google Scholar 

  2. Rohling, E. J. et al. Sea-level and deep-sea-temperature variability over the past 5.3 million years. Nature 508, 477–482 (2014); corrigendum 510, 432 (2014).

    ADS  CAS  PubMed  Google Scholar 

  3. Raymo, M. E., Kozdon, R., Evans, D., Lisiecki, L. & Ford, H. L. The accuracy of mid-Pliocene δ18O-based ice volume and sea level reconstructions. Earth Sci. Rev. 177, 291–302 (2018).

    ADS  CAS  Google Scholar 

  4. De Boer, B., van de Wal, R. S. W., Bintanja, R., Lourens, L. J. & Tuenter, E. Cenozoic global ice-volume and temperature simulations with 1-D ice-sheet models forced by benthic δ18O records. Ann. Glaciol. 51, 23–33 (2010).

    ADS  Google Scholar 

  5. Haywood, A. M. et al. Large-scale features of Pliocene climate: results from the Pliocene Model Intercomparison Project. Clim. Past 9, 191–209 (2013).

    Google Scholar 

  6. Fedorov, A. V. et al. Patterns and mechanisms of early Pliocene warmth. Nature 496, 43–49 (2013).

    ADS  CAS  Google Scholar 

  7. DeConto, R. M. & Pollard, D. Contribution of Antarctica to past and future sea-level rise. Nature 531, 591–597 (2016).

    ADS  CAS  Google Scholar 

  8. Zachos, J. C., Shackleton, N. J., Revenaugh, J. S., Palike, H. & Flower, B. P. Climate response to orbital forcing across the Oligocene-Miocene boundary. Science 292, 274–278 (2001).

    ADS  CAS  Google Scholar 

  9. Pagani, M., Liu, Z. H., LaRiviere, J. & Ravelo, A. C. High Earth-system climate sensitivity determined from Pliocene carbon dioxide concentrations. Nat. Geosci. 3, 27–30 (2010).

    ADS  CAS  Google Scholar 

  10. Lisiecki, L. E. & Raymo, M. E. A Pliocene-Pleistocene stack of 57 globally distributed benthic δ18O records. Paleoceanography 20, PA1003 (2005).

    ADS  Google Scholar 

  11. Dwyer, G. S. & Chandler, M. A. Mid-Pliocene sea level and continental ice volume based on coupled benthic Mg/Ca palaeotemperatures and oxygen isotopes. Phil. Trans. R. Soc. A 367, 157–168 (2009).

    ADS  CAS  Google Scholar 

  12. Naish, T. R. & Wilson, G. S. Constraints on the amplitude of Mid-Pliocene (3.6–2.4 Ma) eustatic sea-level fluctuations from the New Zealand shallow-marine sediment record. Phil. Trans. R. Soc. A 367, 169–187 (2009).

    ADS  CAS  Google Scholar 

  13. Rovere, A. et al. The Mid-Pliocene sea-level conundrum: glacial isostasy, eustasy and dynamic topography. Earth Planet. Sci. Lett. 387, 27–33 (2014).

    ADS  CAS  Google Scholar 

  14. Raymo, M. E., Mitrovica, J. X., O’Leary, M. J., DeConto, R. M. & Hearty, P. J. Departures from eustasy in Pliocene sea-level records. Nat. Geosci. 4, 328–332 (2011).

    ADS  CAS  Google Scholar 

  15. Rowley, D. B. et al. Dynamic topography change of the Eastern United States since 3 million years ago. Science 340, 1560–1563 (2013).

    ADS  CAS  Google Scholar 

  16. Edwards, T. L. et al. Revisiting Antarctic ice loss due to marine ice-cliff instability. Nature 566, 58–64 (2019).

    ADS  CAS  PubMed  Google Scholar 

  17. Miller, K. G. et al. High tide of the warm Pliocene: implications of global sea level for Antarctic deglaciation. Geology 40, 407–410 (2012).

    ADS  CAS  Google Scholar 

  18. Winnick, M. J. & Caves, J. K. Oxygen isotope mass-balance constraints on Pliocene sea level and East Antarctic Ice Sheet stability. Geology 43, 879–882 (2015).

    ADS  CAS  Google Scholar 

  19. Polyak, V. J. et al. A highly resolved record of relative sea level in the western Mediterranean Sea during the last interglacial period. Nat. Geosci. 11, 860–864 (2018).

    ADS  CAS  Google Scholar 

  20. Kendall, R. A., Mitrovica, J. X. & Milne, G. A. On post-glacial sea level—II. Numerical formulation and comparative results on spherically symmetric models. Geophys. J. Int. 161, 679–706 (2005).

    ADS  Google Scholar 

  21. Moucha, R. & Ruetenik, G. A. Interplay between dynamic topography and flexure along the U.S. Atlantic passive margin: insights from landscape evolution modeling. Global Planet. Change 149, 72–78 (2017).

    ADS  Google Scholar 

  22. Stap, L. B. et al. CO2 over the past 5 million years: continuous simulation and new δ11B-based proxy data. Earth Planet. Sci. Lett. 439, 1–10 (2016).

    ADS  CAS  Google Scholar 

  23. De Schepper, S., Gibbard, P. L., Salzmann, U. & Ehlers, J. A global synthesis of the marine and terrestrial evidence for glaciation during the Pliocene Epoch. Earth Sci. Rev. 135, 83–102 (2014).

    ADS  Google Scholar 

  24. Passchier, S. Linkages between East Antarctic Ice Sheet extent and Southern Ocean temperatures based on a Pliocene high-resolution record of ice-rafted debris off Prydz Bay, East Antarctica. Palaeoecanogr. Paleoclimatol. 26, PA4204 (2011).

    ADS  Google Scholar 

  25. Mudelsee, M. & Raymo, M. E. Slow dynamics of the Northern Hemisphere glaciation. Paleoceanography 20, PA4022 (2005).

    ADS  Google Scholar 

  26. Williams, T. et al. Evidence for iceberg armadas from East Antarctica in the Southern Ocean during the late Miocene and early Pliocene. Earth Planet. Sci. Lett. 290, 351–361 (2010).

    ADS  CAS  Google Scholar 

  27. Foster, G. L. & Rohling, E. J. Relationship between sea level and climate forcing by CO2 on geological timescales. Proc. Natl Acad. Sci. USA 110, 1209–1214 (2013).

    ADS  CAS  Google Scholar 

  28. Church, J. A. et al. Sea level change. Climate Change 2013: The Physical Science Basis (Cambridge Univ. Press, 2013).

  29. Morlighem, M. et al. BedMachine v3: complete bed topography and ocean bathymetry mapping of Greenland from multibeam echo sounding combined with mass conservation. Geophys. Res. Lett. 44, 11051–11061 (2017).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  30. Bamber, J. L., Riva, R. E., Vermeersen, B. L. & LeBrocq, A. M. Reassessment of the potential sea-level rise from a collapse of the West Antarctic Ice Sheet. Science 324, 901–903 (2009).

    ADS  CAS  Google Scholar 

  31. Fretwell, P. et al. Bedmap2: improved ice bed, surface and thickness datasets for Antarctica. Cryosphere 7, 375–393 (2013).

    ADS  Google Scholar 

  32. Shakun, J. D. et al. Minimal East Antarctic Ice Sheet retreat onto land during the past eight million years. Nature 558, 284–287 (2018).

    ADS  CAS  Google Scholar 

  33. Rovere, A. et al. The analysis of Last Interglacial (MIS 5e) relative sea-level indicators: reconstructing sea-level in a warmer world. Earth Sci. Rev. 159, 404–427 (2016).

    ADS  Google Scholar 

  34. Decker, D. D. et al. U–Pb dating of cave spar: a new shallow crust landscape evolution tool. Tectonics 37, 208–223 (2018).

    ADS  Google Scholar 

  35. Ludwig, K. R. & Titterington, D. M. Calculation of 230Th/U isochrons, ages, and errors. Geochim. Cosmochim. Acta 58, 5031–5042 (1994).

    ADS  CAS  Google Scholar 

  36. Ludwig, K. R. User's Manual for Isoplot 3.75: A Geochronological Toolkit for Microsoft Excel http://www.bgc.org/isoplot_etc/isoplot/isoplot3_75-4_15manual.pdf (Berkeley Geochronology Center Special Publication 5, 2012).

  37. Schoene, B. et al. Reassessing the uranium decay constants for geochronology using ID-TIMS U–Pb data. Geochim. Cosmochim. Acta 70, 426–445 (2006).

    ADS  CAS  Google Scholar 

  38. Cheng, H. et al. Improvements in 230Th dating, 230Th and 234U half-life values, and U–Th isotopic measurements by multi-collector inductively coupled plasma mass spectrometry. Earth Planet. Sci. Lett. 371–372, 82–91 (2013).

    ADS  CAS  Google Scholar 

  39. Yuan, H. et al. Evaluation of lead isotope compositions of NIST NBS 981 measured by thermal ionization mass spectrometer and multiple-collector inductively coupled plasma mass spectrometer. Solid Earth Sci. 1, 74–78 (2016).

    ADS  Google Scholar 

  40. Dziewonski, A. M. & Anderson, D. L. Preliminary reference Earth model. Phys. Earth Planet. Inter. 25, 297–356 (1981).

    ADS  Google Scholar 

  41. Elderfield, H. et al. Evolution of ocean temperature and ice volume through the Mid- Pleistocene Climate Transition. Science 337, 704–709 (2012).

    ADS  CAS  Google Scholar 

  42. Bierman, P. R., Shakun, J. D., Corbett, L. B., Zimmerman, S. R. & Rood, D. H. A persistent and dynamic East Greenland Ice Sheet over the past 7.5 million years. Nature 540, 256–260 (2016).

    ADS  CAS  Google Scholar 

  43. Shakun, J. D., Lea, D. W., Lisiecki, L. E. & Raymo, M. E. An 800-kyr record of global surface ocean δ18O and implications for ice volume-temperature coupling. Earth Planet. Sci. Lett. 426, 58–68 (2015).

    ADS  CAS  Google Scholar 

  44. Pomar, L. & Ward, W. C. in Sequence Stratigraphy and Depositional Response to Eustatic, Tectonic and Climatic Forcing (ed. Haq, B. U.) Ch. 4, 87–112 (Kluwer, 1995).

  45. Just, J., Hübscher, C., Betzler, C., Lüdmann, T. & Reicherter, K. Erosion of continental margins in the Western Mediterranean due to sea-level stagnancy during the Messinian Salinity Crisis. Geo-Mar. Lett. 31, 51–64 (2011).

    ADS  Google Scholar 

  46. Dorale, J. A. et al. Sea-level highstand 81,000 years ago in Mallorca. Science 327, 860–863 (2010).

    ADS  CAS  Google Scholar 

  47. Creveling, J. R., Mitrovica, J. X., Clark, P. U., Waelbroeck, C. & Pico, T. Predicted bounds on peak global mean sea level during marine isotope stages 5a and 5c. Quat. Sci. Rev. 163, 193–208 (2017).

    ADS  Google Scholar 

  48. Pomar, L. Reef geometries, erosion surfaces and high-frequency sea-level changes, upper Miocene Reef Complex, Mallorca, Spain. Sedimentology 38, 243–269 (1991).

    ADS  Google Scholar 

  49. Pomar, L. La evolución tectonosedimentaria de las Baleares: análisis crítico. Acta Geol. Hisp. 14, 293–310 (1979).

    Google Scholar 

  50. Miller, K. G. et al. The Phanerozoic record of global sea-level change. Science 310, 1293–1298 (2005).

    ADS  CAS  PubMed  Google Scholar 

  51. Flament, N. et al. A review of observations and models of dynamic topography. Lithosphere 5, 189–210 (2013).

    ADS  Google Scholar 

  52. Stocchi, P. et al. MIS 5e relative sea-level changes in the Mediterranean Sea: contribution of isostatic disequilibrium. Quat. Sci. Rev. 185, 122–134 (2018).

    ADS  Google Scholar 

  53. Meehl, G. A. et al. Climate Change 2007: The Physical Science Basis (Cambridge Univ. Press, 2007).

Download references

Acknowledgements

We thank the owner and personnel of Coves d’Artà for granting permission and offering logistic support during the field research conducted for this study. We thank F. L. Forray and G. Lucia for helping with the coring process. B.P.O. and V.J.P. are funded by a collaborative NSF grant (AGS 1602670). Additional research costs were covered by a NSF grant (EAR 0326902 to Y.A. and V.J.P.) and MINECO grants (CGL2013-48441-P and CGL2016-79246-P to J.J.F.). O.A.D. received student research grants from the Cave Research Foundation, the Geological Society of America, and the Fred L. and Helen M. Tharp Endowed Scholarship (School of Geosciences, University of South Florida). J.A. thanks the PALSEA working group and the NSF (grant OCE-0825293 “PLIOMAX”) for facilitating discussions at regular meetings, and the Vetlesen Foundation for support.

Author information

Author notes

  1. These authors contributed equally: Oana A. Dumitru, Jacqueline Austermann

Authors and Affiliations

  1. Karst Research Group, School of Geosciences, University of South Florida, Tampa, FL, USA

    Oana A. Dumitru & Bogdan P. Onac

  2. Department of Earth and Environmental Sciences, Columbia University, Lamont-Doherty Earth Observatory, Palisade, NY, USA

    Jacqueline Austermann

  3. Department of Earth and Planetary Sciences, University of New Mexico, Albuquerque, NM, USA

    Victor J. Polyak & Yemane Asmerom

  4. Earth Sciences Research Group, Universitat de les Illes Balears, Palma, Spain

    Joan J. Fornós, Joaquín Ginés & Angel Ginés

  5. Department of Geology, Babeş-Bolyai University, Cluj-Napoca, Romania

    Bogdan P. Onac

Authors
  1. Oana A. Dumitru
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Jacqueline Austermann
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Victor J. Polyak
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Joan J. Fornós
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yemane Asmerom
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Joaquín Ginés
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Angel Ginés
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Bogdan P. Onac
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

B.P.O., J.J.F. and O.A.D. designed the research. A.G. and J.G. first documented the presence of POS in the cave studied. B.P.O., J.J.F., J.G., A.G. and O.A.D. measured POS elevations and collected samples. V.J.P., O.A.D. and Y.A. measured U–Pb ages of POS samples. J.A. produced the GIA and long-term deformation corrections. O.A.D., J.A. and B.P.O. drafted and wrote the majority of the manuscript with input from all authors.

Corresponding author

Correspondence to Bogdan P. Onac.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

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

Peer review information Nature thanks Eelco Rohling, Jon Woodhead and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Extended data figures and tables

Extended Data Fig. 1 Schematic profile of a coastal cave in Mallorca hosting POS at different levels.

a, b, Standalone (a) and cave-wall (b) POS structures. The asterisked POS is an example of an asymmetric knob-like carbonate encrustation that forms when only the tip of the stalactite is submerged.

Extended Data Fig. 2 U–Pb isochrons and calculated age estimates for the six POS samples.

a, 235U–207Pb two-dimensional isochron for sample AR-02; b–f, Concordia-constrained linear three-dimensional isochron for samples AR-05 (b), AR-11 (c), AR-15 (d), AR-09 (e) and AR-03 (f). Error ellipses on individual ages are 2σ.

Extended Data Fig. 3 GIA contribution due to ongoing adjustment.

This GIA contribution is caused by the incomplete present-day adjustment to the late Pleistocene ice and ocean loading cycles. a, Model simulation using a viscosity structure of 5 × 1020 Pa s viscosity in the upper mantle, 5 × 1021 Pa s viscosity in the lower mantle, and an elastic lithospheric thickness of 96 km. b, Standard deviation of model predictions obtained using 36 different radial viscosity profiles, including varying the lithospheric thickness. The square marks the position of Coves d’Artà. The figures were produced using Matlab 2015b and the m_map plotting package (https://www.eoas.ubc.ca/~rich/map.html).

Extended Data Fig. 4 GIA contribution due to Pliocene ice age cycles.

The model simulation uses a viscosity structure of 5 × 1020 Pa s viscosity in the upper mantle, 5 × 1021 Pa s viscosity in the lower mantle, and an elastic lithospheric thickness of 96 km. a, Snapshot of sea level at 3.244 Ma (grey vertical line in b) assuming a GMSL curve based on the LR04 benthic record10. The colour scale is chosen to diverge around the GMSL value of 13 m. The red square marks the position of Coves d’Artà. The figure was produced using Matlab 2015b and the m_map plotting package (https://www.eoas.ubc.ca/~rich/map.html). b, Local sea level at Coves d’Artà based on Rohling et al.2 (blue), de Boer et al.4 (yellow), and the LR04 benthic record10 (red). Respective GMSL curves are shown in black and mostly coincide with local sea level at Coves d’Artà (note that for the estimates based on Rohling et al.2 the black GMSL curve is mostly behind the local sea-level curve in blue). Sea level is relative to the beginning of this run (4.9 Ma).

Extended Data Fig. 5 Total GIA correction.

a–c, Local sea-level change at Coves d’Artà, as calculated from the GIA model based on the GMSL curves by Rohling et al.2 (a), the LR04 benthic record10 (b) and de Boer et al.4 (c). Uncertainties due to Earth’s viscoelastic structure are denoted by grey bands. d–f, GIA correction colour coded by the GMSL value; standard deviations are shown as grey bands. Black markers indicate the GIA correction and its uncertainty for each POS. Results are for the GMSL curves by Rohling et al.2 (d), the LR04 benthic record10 (e) and de Boer et al.4 (f).

Extended Data Fig. 6 Uplift rate estimation.

Determining the amount of uplift based on the best fit of observed relative sea-level changes across the POS to other GMSL reconstructions over the same time interval. a–c, GMSL curves2,4,10; grey bars are 1σ uncertainties. Boxes indicate the age uncertainty for each POS and the 50th and 99th percentiles of the GMSL values that fall within this age range. We calculate synthetic sea-level changes relative to the youngest POS and compare them to the observed sea-level changes, assuming a range of uplift rates. d–f, Histograms of uplift rates in which we find a good fit between the observed and the synthetic data. Percentiles (16th, 50th and 84th) are shown by vertical lines (solid line is the median, dashed lines are uncertainty bounds). We conducted ten million iterations for this Monte Carlo search. g, Histogram combining all uplift rates that resulted in a good fit.

Extended Data Fig. 7 Joint uplift rate estimation from varying the upper and lower cutoff.

a–i, Joint histograms for a variety of lower and upper percentile cutoffs. The lower cutoff was varied between the 40th, 50th and 60th percentiles (different rows), whereas the upper cutoff was varied between the 90th, 95th and 99th percentiles (different columns). Percentiles in the histograms (16th, 50th and 84th) are shown by vertical lines (solid line is the median, dashed lines are uncertainty bounds). Panel f is identical to Extended Data Fig. 6g. j, Combination of all joint histograms to obtain our best-estimate uplift rate used for Table 1, Fig. 2c and Fig. 3.

Extended Data Fig. 8 Reconstructed elevation of the GMSL for each POS after all corrections have been applied.

a–f, PDFs for the GMSL estimate for AR-02, AR-05, AR-11, AR-15, AR-09 and AR-03, respectively. PDFs are constructed assuming Gaussian uncertainties for the measured elevation of POS, the respective indicative range, the GIA correction, thermal expansion and POS age (Table 1); a non-Gaussian distribution is obtained for the uplift correction (Table 1, Extended Data Figs. 6, 7). The mode (thick black line) and lower and upper uncertainty bounds (16th and 84th percentiles, thick dashed lines) are shown by vertical lines and correspond to the GMSL estimates reported in Table 1. We used a kernel with 1-m bandwidth to calculate the mode and PDF (thin black line).

Extended Data Table 1 GMSL estimates for different percentiles (10th, 33th, 50th, 66th and 90th) following IPCC’s likelihood scale
Full size table

Source data

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Dumitru, O.A., Austermann, J., Polyak, V.J. et al. Constraints on global mean sea level during Pliocene warmth. Nature 574, 233–236 (2019). https://doi.org/10.1038/s41586-019-1543-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41586-019-1543-2

This article is cited by

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Search

Advanced 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