Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

Early Last Interglacial ocean warming drove substantial ice mass loss from Antarctica

View ORCID ProfileChris S. M. Turney, View ORCID ProfileChristopher J. Fogwill, View ORCID ProfileNicholas R. Golledge, Nicholas P. McKay, Erik van Sebille, Richard T. Jones, David Etheridge, View ORCID ProfileMauro Rubino, David P. Thornton, View ORCID ProfileSiwan M. Davies, Christopher Bronk Ramsey, View ORCID ProfileZoë A. Thomas, View ORCID ProfileMichael I. Bird, View ORCID ProfileNiels C. Munksgaard, Mika Kohno, John Woodward, Kate Winter, Laura S. Weyrich, Camilla M. Rootes, Helen Millman, Paul G. Albert, Andres Rivera, View ORCID ProfileTas van Ommen, Mark Curran, Andrew Moy, Stefan Rahmstorf, Kenji Kawamura, View ORCID ProfileClaus-Dieter Hillenbrand, Michael E. Weber, Christina J. Manning, Jennifer Young, and Alan Cooper
  1. aPalaeontology, Geobiology and Earth Archives Research Centre, School of Biological, Earth and Environmental Sciences, University of New South Wales, Kensington NSW 2033, Australia;
  2. bAustralian Research Council Centre of Excellence in Australian Biodiversity and Heritage, School of Biological, Earth and Environmental Sciences, University of New South Wales, Kensington NSW 2033, Australia;
  3. cChronos 14Carbon-Cycle Facility, University of New South Wales, Sydney NSW 2052, Australia;
  4. dSchool of Geography, Geology and the Environment, Keele University, Staffordshire ST5 5BG, United Kingdom;
  5. eAntarctic Research Centre, Victoria University of Wellington, Wellington 6140, New Zealand;
  6. fEnvironment and Climate, GNS Science, Avalon, Lower Hutt 5011, New Zealand;
  7. gSchool of Earth and Sustainability, Northern Arizona University, Flagstaff, AZ 86011;
  8. hGrantham Institute, Imperial College London, London SW7 2AZ, United Kingdom;
  9. iDepartment of Physics, Imperial College London, London SW7 2AZ, United Kingdom;
  10. jInstitute for Marine and Atmospheric Research Utrecht, Utrecht University, 3584 CS Utrecht, The Netherlands;
  11. kDepartment of Geography, Exeter University, Devon EX4 4RJ, United Kingdom;
  12. lClimate Science Centre, Commonwealth Scientific and Industrial Research Organisation Ocean and Atmosphere, Aspendale, VIC 3195 Australia;
  13. mDipartimento di Matematica e Fisica, Università della Campania “Luigi Vanvitelli,” 81100 Caserta, Italy;
  14. nDepartment of Geography, Swansea University, Swansea SA2 8PP, United Kingdom;
  15. oResearch Laboratory for Archaeology and the History of Art, University of Oxford, Oxford OX1 3TG, United Kingdom;
  16. pCentre for Tropical Environmental and Sustainability Science, College of Science and Engineering, James Cook University, Cairns, QLD 4870, Australia;
  17. qAustralian Research Council Centre of Excellence in Australian Biodiversity and Heritage, James Cook University, Cairns, QLD 4870, Australia;
  18. rResearch Institute for the Environment and Livelihoods, Charles Darwin University, Darwin NT 0909, Australia;
  19. sDepartment of Geochemistry, Geoscience Center, University of Göttingen, 37077 Göttingen, Germany;
  20. tDepartment of Geography and Environmental Sciences, Faculty of Engineering and Environment, Northumbria University, Newcastle upon Tyne NE1 8ST, United Kingdom;
  21. uAustralian Centre for Ancient DNA, University of Adelaide, Adelaide SA 5005, Australia;
  22. vAustralian Research Council Centre of Excellence in Australian Biodiversity and Heritage, University of Adelaide, Adelaide SA 5005, Australia;
  23. wDepartment of Geography, University of Sheffield, Sheffield S3 7ND, United Kingdom;
  24. xDepartamento de Geografia, Universidad de Chile, 8331051 Santiago, Chile;
  25. yDepartment of the Environment and Energy, Australian Antarctic Division, Kingston, TAS 7050, Australia;
  26. zAntarctic Climate and Ecosystems Cooperative Research Centre, University of Tasmania, Hobart, TAS 7001, Australia;
  27. aaEarth System Analysis, Potsdam Institute for Climate Impact Research, D-14412 Potsdam, Germany;
  28. bbInstitute of Physics and Astronomy, University of Potsdam, 14476 Potsdam, Germany;
  29. ccResearch Organizations of Information and Systems, National Institute of Polar Research, Tachikawa, Tokyo 190-8518, Japan;
  30. ddDepartment of Polar Science, Graduate University for Advanced Studies, Tachikawa, Tokyo 190-8518, Japan;
  31. eeInstitute of Biogeosciences, Japan Agency for Marine–Earth Science and Technology, Yokosuka 237-0061, Japan;
  32. ffPalaeo Environments, Ice Sheets and Climate Change, British Antarctic Survey, Cambridge CB3 0ET, United Kingdom;
  33. ggSteinmann Institute, University of Bonn, 53115 Bonn, Germany;
  34. hhDepartment of Earth Sciences, Royal Holloway University of London, Surrey TW20 OEX, United Kingdom;
  35. iiSouth Australian Museum, Adelaide, South Australia 5005, Australia

See allHide authors and affiliations

PNAS February 25, 2020 117 (8) 3996-4006; first published February 11, 2020; https://doi.org/10.1073/pnas.1902469117
Chris S. M. Turney
aPalaeontology, Geobiology and Earth Archives Research Centre, School of Biological, Earth and Environmental Sciences, University of New South Wales, Kensington NSW 2033, Australia;
bAustralian Research Council Centre of Excellence in Australian Biodiversity and Heritage, School of Biological, Earth and Environmental Sciences, University of New South Wales, Kensington NSW 2033, Australia;
cChronos 14Carbon-Cycle Facility, University of New South Wales, Sydney NSW 2052, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Chris S. M. Turney
  • For correspondence: c.turney@unsw.edu.au
Christopher J. Fogwill
aPalaeontology, Geobiology and Earth Archives Research Centre, School of Biological, Earth and Environmental Sciences, University of New South Wales, Kensington NSW 2033, Australia;
dSchool of Geography, Geology and the Environment, Keele University, Staffordshire ST5 5BG, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Christopher J. Fogwill
Nicholas R. Golledge
eAntarctic Research Centre, Victoria University of Wellington, Wellington 6140, New Zealand;
fEnvironment and Climate, GNS Science, Avalon, Lower Hutt 5011, New Zealand;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Nicholas R. Golledge
Nicholas P. McKay
gSchool of Earth and Sustainability, Northern Arizona University, Flagstaff, AZ 86011;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Erik van Sebille
hGrantham Institute, Imperial College London, London SW7 2AZ, United Kingdom;
iDepartment of Physics, Imperial College London, London SW7 2AZ, United Kingdom;
jInstitute for Marine and Atmospheric Research Utrecht, Utrecht University, 3584 CS Utrecht, The Netherlands;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Richard T. Jones
kDepartment of Geography, Exeter University, Devon EX4 4RJ, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David Etheridge
lClimate Science Centre, Commonwealth Scientific and Industrial Research Organisation Ocean and Atmosphere, Aspendale, VIC 3195 Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mauro Rubino
dSchool of Geography, Geology and the Environment, Keele University, Staffordshire ST5 5BG, United Kingdom;
lClimate Science Centre, Commonwealth Scientific and Industrial Research Organisation Ocean and Atmosphere, Aspendale, VIC 3195 Australia;
mDipartimento di Matematica e Fisica, Università della Campania “Luigi Vanvitelli,” 81100 Caserta, Italy;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Mauro Rubino
David P. Thornton
lClimate Science Centre, Commonwealth Scientific and Industrial Research Organisation Ocean and Atmosphere, Aspendale, VIC 3195 Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Siwan M. Davies
nDepartment of Geography, Swansea University, Swansea SA2 8PP, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Siwan M. Davies
Christopher Bronk Ramsey
oResearch Laboratory for Archaeology and the History of Art, University of Oxford, Oxford OX1 3TG, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zoë A. Thomas
aPalaeontology, Geobiology and Earth Archives Research Centre, School of Biological, Earth and Environmental Sciences, University of New South Wales, Kensington NSW 2033, Australia;
bAustralian Research Council Centre of Excellence in Australian Biodiversity and Heritage, School of Biological, Earth and Environmental Sciences, University of New South Wales, Kensington NSW 2033, Australia;
cChronos 14Carbon-Cycle Facility, University of New South Wales, Sydney NSW 2052, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Zoë A. Thomas
Michael I. Bird
pCentre for Tropical Environmental and Sustainability Science, College of Science and Engineering, James Cook University, Cairns, QLD 4870, Australia;
qAustralian Research Council Centre of Excellence in Australian Biodiversity and Heritage, James Cook University, Cairns, QLD 4870, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Michael I. Bird
Niels C. Munksgaard
qAustralian Research Council Centre of Excellence in Australian Biodiversity and Heritage, James Cook University, Cairns, QLD 4870, Australia;
rResearch Institute for the Environment and Livelihoods, Charles Darwin University, Darwin NT 0909, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Niels C. Munksgaard
Mika Kohno
sDepartment of Geochemistry, Geoscience Center, University of Göttingen, 37077 Göttingen, Germany;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John Woodward
tDepartment of Geography and Environmental Sciences, Faculty of Engineering and Environment, Northumbria University, Newcastle upon Tyne NE1 8ST, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kate Winter
tDepartment of Geography and Environmental Sciences, Faculty of Engineering and Environment, Northumbria University, Newcastle upon Tyne NE1 8ST, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Laura S. Weyrich
uAustralian Centre for Ancient DNA, University of Adelaide, Adelaide SA 5005, Australia;
vAustralian Research Council Centre of Excellence in Australian Biodiversity and Heritage, University of Adelaide, Adelaide SA 5005, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Camilla M. Rootes
wDepartment of Geography, University of Sheffield, Sheffield S3 7ND, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Helen Millman
aPalaeontology, Geobiology and Earth Archives Research Centre, School of Biological, Earth and Environmental Sciences, University of New South Wales, Kensington NSW 2033, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paul G. Albert
oResearch Laboratory for Archaeology and the History of Art, University of Oxford, Oxford OX1 3TG, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andres Rivera
xDepartamento de Geografia, Universidad de Chile, 8331051 Santiago, Chile;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tas van Ommen
yDepartment of the Environment and Energy, Australian Antarctic Division, Kingston, TAS 7050, Australia;
zAntarctic Climate and Ecosystems Cooperative Research Centre, University of Tasmania, Hobart, TAS 7001, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Tas van Ommen
Mark Curran
yDepartment of the Environment and Energy, Australian Antarctic Division, Kingston, TAS 7050, Australia;
zAntarctic Climate and Ecosystems Cooperative Research Centre, University of Tasmania, Hobart, TAS 7001, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrew Moy
yDepartment of the Environment and Energy, Australian Antarctic Division, Kingston, TAS 7050, Australia;
zAntarctic Climate and Ecosystems Cooperative Research Centre, University of Tasmania, Hobart, TAS 7001, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stefan Rahmstorf
aaEarth System Analysis, Potsdam Institute for Climate Impact Research, D-14412 Potsdam, Germany;
bbInstitute of Physics and Astronomy, University of Potsdam, 14476 Potsdam, Germany;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kenji Kawamura
ccResearch Organizations of Information and Systems, National Institute of Polar Research, Tachikawa, Tokyo 190-8518, Japan;
ddDepartment of Polar Science, Graduate University for Advanced Studies, Tachikawa, Tokyo 190-8518, Japan;
eeInstitute of Biogeosciences, Japan Agency for Marine–Earth Science and Technology, Yokosuka 237-0061, Japan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Claus-Dieter Hillenbrand
ffPalaeo Environments, Ice Sheets and Climate Change, British Antarctic Survey, Cambridge CB3 0ET, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Claus-Dieter Hillenbrand
Michael E. Weber
ggSteinmann Institute, University of Bonn, 53115 Bonn, Germany;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Christina J. Manning
hhDepartment of Earth Sciences, Royal Holloway University of London, Surrey TW20 OEX, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jennifer Young
uAustralian Centre for Ancient DNA, University of Adelaide, Adelaide SA 5005, Australia;
vAustralian Research Council Centre of Excellence in Australian Biodiversity and Heritage, University of Adelaide, Adelaide SA 5005, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alan Cooper
iiSouth Australian Museum, Adelaide, South Australia 5005, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  1. Edited by Johannes Sutter, University of Bern, Bern, Switzerland, and accepted by Editorial Board Member Jean Jouzel December 27, 2019 (received for review February 11, 2019)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Fifty years ago, it was speculated that the marine-based West Antarctic Ice Sheet is vulnerable to warming and may have melted in the past. Testing this hypothesis has proved challenging due to the difficulty of developing in situ records of ice sheet and environmental change spanning warm periods. We present a multiproxy record that implies loss of the West Antarctic Ice Sheet during the Last Interglacial (129,000 to 116,000 y ago), associated with ocean warming and the release of greenhouse gas methane from marine sediments. Our ice sheet modeling predicts that Antarctica may have contributed several meters to global sea level at this time, suggesting that this ice sheet lies close to a “tipping point” under projected warming.

Abstract

The future response of the Antarctic ice sheet to rising temperatures remains highly uncertain. A useful period for assessing the sensitivity of Antarctica to warming is the Last Interglacial (LIG) (129 to 116 ky), which experienced warmer polar temperatures and higher global mean sea level (GMSL) (+6 to 9 m) relative to present day. LIG sea level cannot be fully explained by Greenland Ice Sheet melt (∼2 m), ocean thermal expansion, and melting mountain glaciers (∼1 m), suggesting substantial Antarctic mass loss was initiated by warming of Southern Ocean waters, resulting from a weakening Atlantic meridional overturning circulation in response to North Atlantic surface freshening. Here, we report a blue-ice record of ice sheet and environmental change from the Weddell Sea Embayment at the periphery of the marine-based West Antarctic Ice Sheet (WAIS), which is underlain by major methane hydrate reserves. Constrained by a widespread volcanic horizon and supported by ancient microbial DNA analyses, we provide evidence for substantial mass loss across the Weddell Sea Embayment during the LIG, most likely driven by ocean warming and associated with destabilization of subglacial hydrates. Ice sheet modeling supports this interpretation and suggests that millennial-scale warming of the Southern Ocean could have triggered a multimeter rise in global sea levels. Our data indicate that Antarctica is highly vulnerable to projected increases in ocean temperatures and may drive ice–climate feedbacks that further amplify warming.

  • Antarctic ice sheets
  • marine ice sheet instability (MISI)
  • paleoclimatology
  • polar amplification
  • tipping element

The projected contribution of the Antarctic ice sheet to 21st-century global mean sea level (GMSL) ranges from negligible (1) to several meters (2, 3). Valuable insights into the response of ice sheets to warming may be gained from the Last Interglacial (LIG) (or Marine Isotope Stage [MIS] 5e in marine sediment records; 129,000 to 116,000 y before present or 129 to 116 ky) (4⇓⇓⇓⇓–9). This period experienced warmer polar temperatures and higher GMSL (+6 to 9 m, possibly up to 11 m) (4, 10⇓⇓–13) relative to present day, and was the most geographically widespread expression of high sea level during a previous warm period (4, 10). LIG sea level cannot be fully explained by Greenland Ice Sheet melt (∼2 m) (8), ocean thermal expansion, and melting mountain glaciers (∼1 m) (4), implying substantial Antarctic mass loss (3, 4, 14, 15). Half a century ago, John Mercer was the first to propose that the marine-based West Antarctic Ice Sheet (WAIS) is vulnerable to a warming atmosphere through loss of buttressing ice shelves and may have made a significant contribution to global sea level during the LIG (5⇓–7). Recent work has further demonstrated that extensive deep, marine-based sectors of the East Antarctic Ice Sheet (EAIS) may have accelerated melting, thus contributing to higher LIG sea levels (14). While an isotopic signature of a relatively cool LIG climate preserved in the Mount Moulton blue ice field (16) may be explained by substantial WAIS mass loss (17), no direct physical evidence has yet been identified (4, 18). Temperature estimates derived from climate model simulations provide an indirect measure of change but typically suggest ∼1 °C less warming than proxy-based reconstructions (4, 8, 19). When used to drive ice sheet models, these climate anomalies are not sufficient to remove the floating ice shelves that buttress ice flow from central Antarctica (20). In an attempt to bypass these problems, ice sheet models have been driven by a wide range of prescribed climate scenarios; however, these suggest widely different sensitivities dependent on model physics and parameterization (21, 22), with >2 °C (and in some instances >4 °C) ocean warming required for the loss of the WAIS, exceeding paleoclimate estimates (3, 9, 20, 23) and different sensitivities of Antarctic ice sheet sectors (18, 24, 25).

Here, we report a high-resolution record of environmental change and ice flow dynamics from the Patriot Hills Blue Ice Area (BIA), exposed in Horseshoe Valley (Ellsworth Mountains; Methods) (Fig. 1A). Horseshoe Valley is a locally sourced compound glacier system (i.e., with negligible inflow) that is buttressed by, but ultimately coalesces with, the Institute Ice Stream via the Horseshoe Valley Trough, making the area sensitive to dynamic ice sheet changes across the broader Weddell Sea Embayment (WSE) (26). Due to strong prevailing katabatic airflow, an extensive BIA (more than 1,150 m across) has formed to the leeward side of the Patriot Hills, where ancient ice is drawn up from depth within Horseshoe Valley (Fig. 1E). Regional airborne and detailed local ground-penetrating radar (GPR) surveys show a remarkably coherent series of dipping (24 to 45°) layers, broken by two discontinuities, which represent isochrons across the Patriot Hills BIA, extending thousands of meters into Horseshoe Valley. A “horizontal ice core” across the BIA spans the time intervals 0 to 80 ky and 130 to 134 ky (Methods and SI Appendix, Fig. S5) constrained by analysis of trace gases and geochemically identified volcanic layers exposed across the transect, which have been Bayesian age modeled against the recently compiled continuous 156-ky global greenhouse gas time series (CO2, CH4, and N2O) (27) on the AICC2012 age scale (28) (Fig. 1B and Methods). The record is located 50 km inland from the modern grounding line of the Filchner–Ronne Ice Shelf in the WSE (29) and close to the Rutford Ice Stream, one of the largest methane hydrate reserves identified in Antarctica [total organic carbon estimated to be 21,000 Gt (30), equivalent to ∼2,000 y of the current carbon emission rate of 10 GtC/year (https://www.co2.earth/global-co2-emissions)]. Today, precipitation at the site is delivered via storms originating from the South Atlantic or Weddell Sea (31). Crucially, the Ellsworth Mountains also lie in a sector of the continent that is highly responsive to isostatic rebound under a scenario of substantial WAIS mass loss, potentially preserving ice from around the time of the LIG in small valley glaciers and higher ground areas (32).

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Location and age profile of the Patriot Hills BIA. (A) Location of Antarctic ice and marine records discussed in this study and austral spring–summer (October to March) SST trends (over the period 1981 to 2010; HadISST data). (B) Trace gas (circles), tephra (triangles), and boundary (square) age solutions for surface ice along transect B–B′ relative to an arbitrary datum along the transect (displayed in D). The dashed lines denote unconformities D0–D2 at their surface expression. (C) Basal topography of the Ellsworth Subglacial Highlands (West Antarctica) with the locations of airborne radio-echo sounding transect A–A′ (displayed in E) and Rutford Ice Stream (IS) (29). The Horseshoe Valley, Independence, and Ellsworth troughs are given by the initials HV, IT, and ET, respectively. (D) The location of Patriot Hills in Horseshoe Valley (LIMA background image) with the BIA climate line (marked by transect B–B′), dominant ice flow direction, and distance to grounding line. (E) Airborne radio-echo sounding cross-section of ice within Horseshoe Valley, Independence, and Ellsworth troughs (modified from ref. 29). Digitization highlights basal topography (brown), lower basal ice unit (gray), and upper basal ice unit (red) as well as internal stratigraphic features (black for observed, dashed for inferred, and purple for best estimate).

The Patriot Hills Record

The isotopic series of δD across the Patriot Hills BIA exhibits a coherent record of relatively low values between 18 and 80 ky, consistent with a glacial-age sequence (Fig. 2E). Below these layers and at the periphery of zones of higher ice flow (29), we find an older unit of ice exposed at the surface expressed by a step change to enriched (interglacial) isotopic values (Fig. 2E and SI Appendix, Fig. S7), implying proximal warmer conditions and reduced sea ice extent (33). Importantly, we identify a distinct tephra horizon near the boundary of this older unit of ice, which, based on major and trace element geochemical fingerprinting (Fig. 3 and SI Appendix, Fig. S11), is correlated to a volcanic ash from the penultimate deglaciation (Termination II) referred to as Tephra B in marine sediments on the West Antarctic continental margin (34) and identified at 1,785.14-m depth in the Dome Fuji ice core, where it is dated to 130.7 ± 1.8 ky (AICC2012 timescale) (28, 33, 34). The start of the oldest section of the sequence is dated here to 134.1 ± 2.2 ky, consistent with modeling studies, airborne radio-echo sounding lines, and GPR profiles, which imply older ice exists at depth in the Ellsworth Mountains (29, 32) (Fig. 1 B–E).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Climate, ocean circulation, and sea level changes over the past 140 ky. (A) δ18O record from the North Greenland (NGRIP) ice core (106, 107). (B) Bermuda Rise 231Pa/230Th data (reversed axis; 1σ uncertainty) with dashed horizontal line denoting production ratio of 0.093 marking sluggish/absent AMOC (42). Selected North Atlantic Heinrich (H) events and reduced AMOC shown. (C) Biogenic opal flux from ODP Site 1094 (53.2°S) as a measure of wind-driven upwelling in the Southern Ocean (45). (D) Comparison between the recently compiled global atmospheric methane time series (red line; 2σ envelope) (27) with the methane record from the West Antarctic Patriot Hills (black circles with 1σ uncertainty; open circles mark anomalously high-concentration data excluded from age model; Methods). (E) The Patriot Hills record. Pie chart representation (circle and segments) of percentage methane-utilizing bacteria in 16S rRNA samples from Patriot Hills; crosses denote absence of these bacteria (Methods). Triangles denote the presence of geochemically identified tephra layers in the Patriot Hills transect, with δD (and mean Holocene; blue envelope 1σ) values. The gray shading denotes the timing of the surface elevation change across the WSE as indicated by the hiatus in the Patriot Hills sequence and inferred substantial Antarctic ice mass loss, consistent with the reported divergence of the isotopic signal observed between the horizontal Mount Moulton ice core record from the WAIS and East Antarctic ice cores (16, 17, 33), and peak global sea level (10). (F) Temporal changes in ocean productivity with peak productivity (PP) (green shading) during interglacials and subsequent enhanced content of calcareous microfossils in Antarctic continental margin sediments (red shading) (34). The dashed black line shows position of tephra identified in the Patriot Hills (−340 m), Dome Fuji (1,785.14 m), and Tephra B in marine sediments from the West Antarctic continental margin. (G) East Antarctic Dome Fuji δ18O record (28, 33). (H) Reconstructed relative sea level curve with 2σ envelope (10). The yellow shading highlights the timing of iceberg-rafted Heinrich debris event 11 (H11), when large amounts of iceberg-rafted debris were deposited in the North Atlantic (43) and the 231Pa/230Th ratio on Bermuda Rise shifted toward the production ratio of 0.093, representative of sluggish or absent AMOC (42); the circled numbers 1 and 2 denote enhanced upwelling-induced warming in the Southern Ocean and Antarctic ice mass loss, respectively. (I) Close-up of Termination II and the onset of the LIG highlighting the high-precision correlation enabled by the Patriot Hills tephra (∼130 ky) and the carbon isotopic composition of benthic foraminifera from ODP Site 1089 and ODP Site 1094 (46) (Fig. 1A). The cream shading highlights the inferred collapse of the AABW reported from ODP 1094 (46). Dashed vertical line denotes LIG tephra in Patriot Hills, Dome Fuji, and West Antarctic continental margin.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Average trace element concentrations of Patriot Hills tephra at −340 m and Tephra B from marine sediment cores PC108 (4.65-m depth) and PC111 (6.86-m depth) (34) normalized to Primitive Mantle (108). (A) Biplots show comparison between selected trace element concentrations of the tephra in the different sequences. Error bars on plots show 2σ of replicate analyses of MPI-DING StHs6/80-G (87), but errors are typically smaller than the data symbols (B–E).

The combined tephra and trace gas analyses suggest a ∼50-ky hiatus after Termination II (130.1 ± 1.8 ky). Radio-echo sounding surveys across the WSE have identified a large subglacial basin comprising landforms reflecting restricted, dynamic, marine-proximal alpine glaciation, with hanging tributary valleys feeding an overdeepened Ellsworth Trough (35). The extensive nature of the subglacial features implies substantial and repeated mass loss of the marine sections of the WAIS (presumably through the Pleistocene), with the ice margin some 200 km inland of present day (35). However, the timing of most recent retreat is currently unknown. While previous surface exposure dating in the region has suggested that the WAIS contribution to global sea level rise during warmer periods was limited to 3.3 m above present (36), relatively short-duration interglacial periods may have resulted in near-complete deglaciation (35). Previous work has interpreted erosional features D1 and D2 in the Patriot Hills BIA to be a consequence of extensive ice surface lowering in Horseshoe Valley (up to ∼500 m since the Last Glacial Maximum, 21 ky) and more exposure of katabatic-enhancing nunataks, resulting in increased wind scour (26, 37). While this scenario may explain unconformity D0, previous work has demonstrated Horseshoe Valley and the wider WSE to be highly sensitive to periods of ice stream advance or retreat in the last glacial cycle and Holocene, with dramatic reductions in surface elevation (26, 37⇓–39), changes that may result in more than just increased wind scour. Importantly, the head of Horseshoe Valley is an overdeepened trough (down to ∼2,000 m below sea level), while toward the mouth of the valley, a subglacial ridge is found at ∼200 m below current sea level with an ice thickness of some 750 m (Methods and SI Appendix, Fig. S3), allowing the isolation and stagnation of ice in Horseshoe Valley over multiple millennia. Furthermore, glaciological investigations assessing the impact of ice shelf loss on glaciers along the Antarctic Peninsula provide important insights into the preservation of ice, albeit on a smaller scale. The 2002 Larsen B ice shelf collapse led to many of the tributary glaciers abruptly changing from a convex to a concave profile (cross-section) (40), with relict ice left isolated on the upper flanks of the valleys (41). These scenarios are consistent with extensive grounding line retreat across the inner shelf of the Weddell Sea and associated substantial ice loss across the wider WSE (29).

The ice at Patriot Hills therefore appears to preserve a record of glacier flow in Horseshoe Valley up to the moment when the Filcher–Ronne Ice Shelf collapsed, after which the sequence remained isolated due to regional ice flow reconfiguration for multiple millennia; a situation that persisted until the ice surface had risen sufficiently to enable the regional ice flow to recover sometime during late MIS 5. We cannot, however, discount the possibility that there were one or more cycles of ice mass gain and loss through MIS 5. The presence of a discrete older ice unit along the flanks of the Ellsworth Mountains (29) (Fig. 1 and SI Appendix, Fig. S2) and the subsequent inferred highly variable climate and/or sea ice extent across the wider WSE (SI Appendix, Figs. S7 and S13) imply the preservation of ice from MIS 6/5 (Termination II) and 5/4 transitions in Horseshoe Valley. Our data provide evidence for substantial mass loss across the WSE during the LIG.

Ocean Warming

What could be the cause of this ice loss in the South Atlantic sector of the Southern Ocean? Recent work has proposed that the iceberg-rafted Heinrich 11 event between 135 and 130 ky (during Termination II) may have significantly reduced North Atlantic Deep Water (NADW) formation and shut down the Atlantic meridional overturning circulation (AMOC) (42), resulting in net heat accumulation in the Southern Hemisphere (the bipolar seesaw pattern of northern cooling and southern warming) (43, 44) (Fig. 4A). Under this scenario, surface cooling during Heinrich 11 increased the northern latitudinal temperature gradient and caused a southward migration of the Intertropical Convergence Zone and midlatitude Southern Hemisphere westerly airflow (14, 45). Importantly, Heinrich 11 was probably one of the largest of the iceberg-rafting events over the last 140 ky (including H-1 and H-2) and during a time of likely weakened AMOC (42). In the Southern Ocean, the associated northward Ekman transport of cool surface waters (something akin to today; Fig. 1A) was likely compensated by increased delivery of relatively warm and nutrient-rich Circumpolar Deep Water (CDW) toward the Antarctic margin (14, 34, 43, 45, 46), potentially leading to enhanced thermal erosion of ice at exposed grounding lines (43, 47). This interpretation is supported by the enriched benthic foraminifera 13C values into the LIG (46), a proxy for the influence of NADW on CDW in the south, implying northern (warmer) waters were reaching far south for much of this period (and a cause of persistent loss of ice volume) (Fig. 2I). The unambiguous precise correlation between the Patriot Hills ice and West Antarctic marine records (34) afforded by the Termination II tephra demonstrates that the warming recorded in the BIA is coincident with a major, well-documented peak in marine temperatures and productivity around the Antarctic continent and in the Southern Ocean (34, 45, 46) (Fig. 2). The subsequent delivery of large volumes of associated freshwater into the Southern Ocean during the LIG would have reduced Antarctic Bottom Water (AABW) production (46), resulting in increased deepwater formation in the North Atlantic (43, 48, 49) (Fig. 4C). Recent modeling results suggest that increased heat transport beneath the ice shelves can drive extensive grounding-line retreat, triggering substantial drawdown of the Antarctic ice sheet (2, 14, 20) (Fig. 4B). Of concern, warming of the ocean cavity in the WSE is projected to increase during the 21st century (50).

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Ocean–atmospheric interactions during Termination II and the LIG. Panels show changing Atlantic meridional overturning circulation (AMOC) in response to iceberg discharge (A and B) in the North Atlantic (Heinrich event 11) during Termination II and (C) from the Antarctic Ice Sheet (AIS) during the LIG, with inferred shifts in atmospheric circulation including midlatitude Southern Hemisphere westerly (crossed circle) airflow and Intertropical Convergence Zone (ITCZ) (14, 43, 45, 46, 48). The vertical arrows denote CH4 and heat flux associated with Antarctic coastal easterly (dot in circle) and westerly (crossed circle) airflow (30, 47). AABW, AAIW, CDW, NAIW, and NADW define Antarctic Bottom Water, Antarctic Intermediate Water, Circumpolar Deep Water, North Atlantic Intermediate Water, and North Atlantic Deep Water, respectively.

With Southern Ocean warming and concurrent ice sheet retreat, the large methane reservoirs in Antarctic sedimentary basins (e.g., Rutford Ice Stream) could have become vulnerable to release (30) and may have contributed to elevated atmospheric levels through the LIG (8, 27) (Fig. 2D). High-latitude open water and sea ice are rich in microbial communities, components of which may be collected by passing storms and delivered onto the ice sheet (e.g., prokaryotes, DNA), offering insights into offshore environmental processes (51, 52). To investigate environmental changes prior to and after the ice sheet reconfiguration recorded in the Patriot Hills BIA, we applied an established ancient DNA methodology and sequencing to provide a description of ancient microbial species preserved within the ice (Methods). Methane-utilizing microorganisms were found in three samples along the Patriot Hills transect and were absent from other samples on the transect and laboratory controls. While such microbes are not obligate methylotrophs and can be present in nonmethane-dominated environments (53), they would be expected to be at very different abundances to what we find. The most striking feature of the Patriot Hills BIA genetic record was detected immediately prior to inferred ice loss, where Methyloversatilis microbes dominated the detectable microbial diversity (∼130 ky) (Fig. 2E and SI Appendix, Fig. S15). Methyloversatilis was only found in high abundance in this sample (with trace amounts identified at ∼22 ky). Crucially, Methyloversatilis are facultative methylotrophs and live on single and multicarbon sources (54), consistent with elevated levels of CH4 and active methane oxidation by Methyloversatilis or other methanotrophic taxa in marine sediments or in the water column during the end of Termination II (SI Appendix). More work is needed to explore the potential for microbial methane utilization in this unique environment.

Antarctic Ice Sheet Modeling

The inferred substantial mass loss across the WSE implies a major role for ocean warming during Termination II and the LIG. To provide a framework for interpreting ice sheet dynamics around the Patriot Hills and across Antarctica, we present a series of temperature sensitivity experiments using the Parallel Ice Sheet Model, version 0.6.3 (Fig. 5) (2). We report here nine different simulations that capture a range of ocean and atmospheric warming scenarios (0° to 3 °C). Importantly, the most comprehensive published high-latitude (≥40° S) network of quantified sea surface temperature (SST) estimates suggests an early LIG (∼130 ky) warming of 1.6 ± 0.9 °C relative to present day (9, 23), providing an upper limit on the sensitivity of the Antarctic ice sheet to ocean temperatures. The pattern of circum-Antarctic ocean warming during this time period is not well established so we assume a spatially uniform warming pattern relative to present day temperatures. Our model time series illustrates that the majority of ice loss takes place within the first two millennia, depending on the magnitude of the forcing (Fig. 5 and Table 1). This corresponds to the time period of inferred loss of marine-based sectors of the ice sheet (Fig. 2), primarily in West Antarctica. In contrast to some whole-continent models, our simulations do not include mechanisms by which a grounded ice cliff may collapse (3), a process that produces considerably faster and greater ice margin retreat than reported here.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Modeled Antarctic ice sheet evolution under idealized forcing scenarios consistent with range of inferred LIG temperatures. (A) Sea level equivalent mass loss for ice sheet simulations forced by a range of air and ocean temperature anomalies relative to present day. “dT” and “dOT” describe atmospheric and ocean temperature anomalies, respectively. B–D show Antarctic Ice Sheet extent and elevation with 2 °C warmer ocean temperatures over time intervals of 1, 2, and 5 ky, respectively (with no atmospheric warming); equivalent sea level contribution is given in the Bottom Left corner of each panel. Locations of Patriot Hills (Ellsworth Mountains, WAIS) and ice core records discussed in this study are shown in B. Inset box in B outlines region shown in Fig. 6.

View this table:
  • View inline
  • View popup
Table 1.

Sea level equivalent mass loss (meters) for Antarctic ice sheet simulations forced over 10,000 y by range of annual air and ocean temperature anomalies relative to present day

For the 2 °C warmer than present day ocean temperature scenario (comparable to reconstructed estimates) (9, 23), with no additional atmospheric warming, our model predicts a contribution to GMSL rise of 3.8 m in the first millennium of forcing (Fig. 5B). The loss of the Filchner–Ronne Ice Shelf within 200 y of warming triggers a nonlinear response by removing the buttressing force that stabilizes grounded ice across large parts of the WSE and the EAIS (most notably the Recovery Basin) (Fig. 6 and SI Appendix, Fig. S17). Ongoing slower ice loss subsequently occurs around the margins of East Antarctica, producing a sustained contribution to sea level rise. Even for relatively cool ocean-forced runs, we find the shelves collapse quickly between the 200-y intervals (SI Appendix, Fig. S18). Indeed, during the warmer ocean model runs, the shelves disappear too quickly to observe the relevant processes on the timescale covered by the snapshots. For instance, under the scenario of 2 °C linear warming, the ice shelves disappear within 600 y of forcing (when temperatures reached between +0.4 and +0.8 °C). Other modeling studies using a range of different setups have reported similar rapid losses of the ice shelves during the onset of the LIG (24, 25). Our results are therefore consistent with an increasing body of evidence that the stability of Antarctic ice shelves is vulnerable to a relatively low temperature threshold (2, 24, 25).

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

Bed (black line) and surface (blue) elevation changes at Patriot Hills (Ellsworth Mountains, WAIS) in response to 2 °C warmer ocean temperatures over a time interval of 5 ky (with no atmospheric warming) (A). Bed (black line) and surface (blue) elevation changes vs. time, with phases of the prevalence of particular processes, such as ice shelf collapse (mint shaded), regional uplift (gray shaded), and dynamic thinning (light-brown shaded), highlighted. (B–G) Selected time slices corresponding to dashed lines in A showing ice shelf extent and ice sheet elevation in the Weddell Sea Embayment (WSE) over the first 3 ky. Location of Patriot Hills is marked by the red square; the gray shaded areas are ice shelf covered, while the white areas are free of both grounded and floating glacial ice.

Recent work has suggested that the Ellsworth Mountains would have experienced a relatively large positive isostatic adjustment (>200 m) accompanying the loss of the WAIS (24, 25, 32), although the model outputs may be underestimated (25). To investigate how an evolving ice sheet geometry would manifest across the wider region, we extracted local ice surface and bed elevations for the WSE from the model simulation that uses a preindustrial ice sheet configuration with 2 °C ocean warming and no atmospheric warming. Fig. 6 A–G illustrates the sequence of events that take place as the ice sheet evolves. First, loss of the Filchner–Ronne Ice Shelf in the Weddell Sea triggers a nonlinear response, removing the buttressing force that stabilizes grounded ice across large parts of the WSE and the EAIS (most notably the Recovery Basin) (55). The loss of back-stress allows for an acceleration of grounded ice and a rapid but short-lived thinning episode (32). At the Patriot Hills, bedrock uplift of ∼30 m over this 0.2-ky period is outpaced by a surface lowering of ∼75 m, implying a net ice sheet thinning of around 105 m. Subsequently, regional-scale isostatic uplift elevates both the bed topography (∼250 m) and ice sheet surface (∼350 m) relative to the initial configuration. The difference between these two values reflects positive net mass balance of the ice sheet here (∼0.055 m/y). After around 2.5 ky, renewed dynamic thinning of the ice sheet in the Patriot Hills leads to a rapid thinning and lowering of the ice sheet surface, at a rate exceeding regional-scale bedrock subsidence (120 m over 0.4 ky, or 0.3 m/y, compared to ∼70 m over 3.2 ky, or 0.022 m/y, respectively) (Fig. 6). For the 1 and 3 °C warming scenarios, similar spatial losses are modeled, with GMSL rises of 2.2 and 4.7 m for the first millennium, respectively (Table 1). Atmospheric warming of the magnitude suggested by Antarctic cores (>4 °C) (16, 17, 56⇓–58) adds an additional meter of equivalent global sea level within the first millennium (SI Appendix, Fig. S19).

Previous work has highlighted the sensitivities of the Ronne–Filchner and Ross ice shelves to warming under a range of model setups (3, 18). Recently published transient ice sheet model simulations covering the last glacial–interglacial cycle have investigated a range of scenarios encompassing different geothermal heat fluxes, ice shelf calving heights, mantle viscosity values, temperature and sea level forcing scenarios, etc. (24, 25). Importantly, these studies recognize the loss of the WAIS forced by warming across what is relatively narrow LIG temperature peak, with a maximum bedrock elevation of ∼400 to 650 m, and surface elevation changes of >1,500 m, larger than that reported here. However, it is important to note that these relatively large estimates are likely influenced by the glacial loading that was experienced during MIS 6.

The Patriot Hills record is consistent with basin-scale mass loss early in the LIG (15, 32) as a consequence of regional ice dynamic changes and isostatically driven isolation of Horseshoe Valley from sustained ocean forcing. While some modeling studies have argued the loss of the Filchner–Ronne Ice Shelf does not display a strong marine ice sheet instability feedback (59) and that isostatically driven rebound may halt ice retreat (18), our results suggest otherwise. Indeed, recent work has proposed that if mass loss comparable to recent decades is maintained for as little as 60 y, the WAIS could be irrevocably destabilized over subsequent millennia through the collapse in the Amundsen Sea sector, overcoming any isostatically driven rebound (60). Future work will be required to undertake large ensembles of high-resolution ice sheet model simulations that capture the full range of ice dynamics, ice–ocean–atmosphere coupling, MIS 6 ice sheet configuration, and spatial and temporal temperature evolution across this period to fully capture the uncertainty associated with LIG mass loss. However, we consider our ice sheet modeling simulations to be comparable to previous studies (24, 25) in the magnitude of rate of change and mass loss, and support the interpretation of the Patriot Hill BIA record. Our results suggest substantial ice sheet mass loss and flow reconfiguration in response to ocean warming, outpacing any bedrock rebound that might have stabilized the ice sheet (Fig. 6). Furthermore, marine-based ice sheets are particularly vulnerable to hysteresis effects (61), which could explain the 50-ky hiatus in the Patriot Hills blue ice record, particularly given the relatively low modeled temperature threshold (0.5 to 0.7 °C ocean warming) for ice shelf loss (SI Appendix, Fig. S17).

The evidence for substantial mass loss from Antarctica in the early LIG has important implications for the future (4, 62). Our field-based reconstruction and modeling results support a growing body of evidence that the Antarctic ice sheet is highly sensitive to ocean temperatures. Driven by enhanced basal melt through increased heat transport into cavities beneath the ice shelves (2, 47), this process is projected to increase with a weakening AMOC during the 21st century (50, 63⇓–65), which may lead to other positive feedbacks such as destabilization of methane hydrate reserves (30).

Methods

Patriot Hills.

Site description and geomorphological context.

The Patriot Hills BIA (Horseshoe Valley, Ellsworth Mountains; 80°18′S, 81°21′W) is a slow flowing (<12 m⋅y−1) compound glacier system situated within an overdeepened catchment that coalesces with the Institute Ice Stream at the periphery of the WSE (29, 37, 66⇓–68) (Fig. 1 and SI Appendix, Figs. S1–S4). Airborne radio-echo sounding surveys across the Ellsworth Mountains have revealed several wide (up to 34 km across) and long (260 km) subglacial troughs containing ice up to 2,620 m thick (Fig. 1) (29), along the side of which, two radar zones have been interpreted to indicate layers of ice with contrasting physical properties, consistent with snow deposited during previous glacial/interglacial transitions. In contrast to the other troughs across the Ellsworth Mountains, contemporary ice within the Horseshoe Valley Trough maintains the slowest average flow speeds of all, at 12 m⋅a−1 (cf. the main trunk of the Institute Ice Stream reaches speeds up to 415 m⋅a−1). This is in large part due to the configuration of the Horseshoe Valley Trough where the ice thickness measures in excess of 2,000 m at the head of the valley and reduces to ∼1,400 m downstream; toward the mouth of the valley, a subglacial ridge is found at ∼200 m below sea level with the ice thickness some 750 m thick (SI Appendix, Fig. S3) (69). The new Digital Elevation Model data for the WSE is available at https://data.bas.ac.uk/full-record.php?id=GB/NERC/BAS/PDC/00937. The configuration of the bed and resulting slow flow in Horseshoe Valley has two major benefits for our study. It allows 1) a long record of ice to accumulate, and 2) the isolation and preservation of ice during periods of regional and Antarctic-wide mass loss.

In the lee of a small mountain chain at the end of Horseshoe Valley called Patriot Hills, strong local katabatic winds descend into the valley from the polar plateau, ablating the ice sheet surface by up to 170 kg⋅m−2⋅y−1 (68). As a result, ancient ice is drawn up from depth in the Horseshoe Valley Trough to form an extensive BIA (more than 1,150 m across; SI Appendix, Fig. S4) (31, 37, 38). High-resolution analysis using GPR (37) and isotopes identifies three distinct unconformities [surface distances relative to an arbitrary transect datum (31) set at zero]: 247 m (D1), 360 m (D2), and −339 m (D0). Based on the trace gas, tephra, and isotopic values of the surface ice beyond D0 (closest to Patriot Hills), we interpret this section of the record to be Termination II in age (see below). No glaciomarine sediments have been identified at any of the boundaries.

Previous work has interpreted erosional features D1 and D2 in the Patriot Hills BIA to be a consequence of extensive ice surface lowering in Horseshoe Valley (up to ∼500 m since the Last Glacial Maximum, 21 ky) and more exposure of katabatic-enhancing nunataks, resulting in increased wind scour (26, 37). While this scenario may explain unconformity D0, other studies have demonstrated Horseshoe Valley and the wider WSE to be highly sensitive to periods of rapid ice stream advance or retreat in the last glacial cycle and Holocene with dramatic reductions in surface elevation (26, 37⇓–39). Recent work investigating the impact of ice shelf loss on glaciers along the Antarctic Peninsula provides important insights, albeit on a smaller scale. The 2002 Larsen B ice shelf collapse led to many of the tributary glaciers abruptly changing from a convex to a concave profile (40), with relict ice left isolated on the upper flanks of the valleys (41). Under a scenario of extreme ice surface lowering arising from ocean warming during the early LIG, the ice at Patriot Hills preserves a record of glacier flow in the overdeepened Horseshoe Valley up to the moment when the Filcher–Ronne Ice Shelf collapsed, after which the sequence likely remained isolated for multiple millennia until the ice surface had risen sufficiently to reincorporate the isolated ice into the glacier sometime during late MIS 5. The relatively enriched deuterium and 18O stable isotope values, ancient DNA (notably the detection of Methyloversatilis microbes in the sample form −340 m in the Patriot Hills record), and ice sheet modeling are consistent with early offshore warming in the south Atlantic and substantial ice mass loss in the early LIG (34, 46, 62), preserving most (if not all) of the Termination II ice record during the period represented by the D0 unconformity (see below). We therefore consider D0 reflects a significant fall in surface elevation and change in flow direction due to isostatically driven isolation of the valley during a period of rapid drawdown of the ice streams across the WSE.

Chronology.

Chronological control across the transect is provided by a comprehensive suite of trace gas samples—carbon dioxide (CO2), methane (CH4) and nitrous oxides (N2O)—and volcanic tephra horizons. The trace gas measurements provide a range of possible age solutions against the recently published 156-ky smoothed global time series for these gas species (27), which together with the absolute constraints provided by the tephra horizons, allows the development of a robust chronological framework that can be tied directly to the isotopic series through high-resolution GPR (31, 37) (SI Appendix, Figs. S6 and S7). A Kovacs 9-cm-diameter ice corer was used to collect ice for gas and taken from >3-m depth to minimize modern air contamination and/or alteration (31). The samples were double bagged and sealed in the field, and transported frozen to the Commonwealth Scientific and Industrial Research Organisation (CSIRO) ICELAB facility in Melbourne for the extraction and measurement of trace gases using a modified dry extraction “cheese grater” and cryogenic trapping technique (70, 71). The trapped air samples were analyzed by gas chromatography, and the trace gas concentrations are reported against the calibration scales maintained by CSIRO GASLAB (72). Where sufficient material was available, duplicates were analyzed.

The presence of visible tephra layers (volcanic ash horizons) provides additional chronological control for the Patriot Hills BIA. Here, we report two new tephras from Patriot Hills at 10 and −340 m, both observed as ∼4-cm units of dispersed shards (SI Appendix, Fig. S9). Shards were extracted by centrifugation of the melted ice samples and put onto a glass slide for electron microprobe analysis. The slides were ground and polished using silica carbide paper and decreasing grades of diamond suspension to expose fresh sections of glass. Single-grain analyses of 10 oxides were performed on a Cameca SX-100 electron microprobe at the Tephrochronology Analytical Unit, University of Edinburgh. See SI Appendix for operating conditions (73); geochemical results are provided in SI Appendix, Table S1. The shards from 10 m are bimodal, with a basanitic and trachytic composition (SI Appendix, Fig. S10). The shards from −340 m are trachytic in composition and exhibit a tightly clustered population (SI Appendix, Fig. S11). Both were compared to published tephras from across Antarctica (34, 74⇓⇓⇓⇓⇓⇓⇓⇓⇓–84). The 10-m tephra has the closest match to be the basanite Tephra C from the WAIS Divide at 3,149.12 m (Similarity Coefficient or SC = 0.98), equivalent to 44.9 ± 0.3 ky (84). The −340-m tephra revealed the closest match to a tephra layer in the Dome Fuji ice core at 1,785.14-m depth [SC = 0.966; equivalent to 130.7 ± 1.8 ky on the AICC2012 timescale (28, 77, 85); data previously unpublished].

A widespread tephra found in marine sedimentary records on the West Antarctic continental margin (Tephra B) has been proposed to correlate to the tephra at Dome Fuji 1,785.14 m, but the correlation has until now remained only tentative in the absence of any reported geochemistry from the latter (34). Here, we find the major oxides from Tephra B have a close match to Patriot Hills −340 m (SC = 0.948), consistent with this interpretation. To test this correlation, we undertook trace element analysis of the glass shards from Patriot Hills at −340 m. Unfortunately, the Dome Fuji shards were too thin for analysis. However, we were able to undertake trace element analyses on Tephra B samples from two marine sediment cores from the West Antarctic continental margin: PC108 (4.65-m depth) and PC111 (6.86-m depth) (34). Trace element analysis of volcanic glass shards were performed using an Agilent 8900 triple-quadrupole inductively coupled plasma mass spectrometry (ICP-MS) (ICP-QQQ) coupled to a Resonetics 193-nm ArF excimer laser ablation in the Department of Earth Sciences, Royal Holloway, University of London. See SI Appendix for operating conditions (86). Accuracies of laser ablation ICP-MS analyses of ATHO-G and reference StHs6/80-G MPI-DING (87) glass were typically ≤5%. Identical trace element glass chemistries (Fig. 2 and SI Appendix, Table S2) strongly support the correlation of Patriot Hills −340-m tephra horizon and the marine West Antarctic Tephra B (34), which is in turn correlated to Dome Fuji 1,785.14 m (33, 34, 77, 85), and probably originates from the Marie Byrd Land volcanic province (West Antarctica) (34). The recognition of a widespread tephra horizon across a large sector of the Antarctic at the very onset of the LIG provides a time-parallel marker horizon crucial for future studies investigating Antarctic ice sheet mass loss.

To develop an age model, we undertook Bayesian age modeling using a Poisson process deposition model (P_sequence) in the software package OxCal, version 4.2.4 (https://c14.arch.ox.ac.uk/) (SI Appendix, Tables S3 and S4) (88, 89). Using Bayes theorem, the algorithms employed sample possible solutions with a probability that is the product of the prior and likelihood probabilities (90, 91). “Calibration curves” with 20-y resolution were developed for the three trace gas species using the 156-ky time series (27). Taking into account the deposition model, the reported ages of the tephra layers, and the common age solutions offered by the trace gas measurements, the posterior probability densities quantify the most probable age distributions. The available constraints suggest the 1,156-m-long Patriot Hills BIA transect spans time intervals from ∼134.2 to ∼1.3 ky comprising four key zones: 4 (−362 to −339 m, equivalent to 134.2 ± 2.2 to 130.1 ± 1.8 ky), 3 (−326 to 240 m, equivalent to 80 ± 6.1 to 22.7 ± 2.8 ky), 2 (240 to 360 m, equivalent to 22.7 ± 2.8 to 10.3 ± 0.4 ky), and 1 (360 to 800 m, 10.3 ± 0.4 to 1.3 ± 0.6 ky). The Agreement Index (a measure of the agreement between the model—prior—and the observational data—likelihood) for the Patriot Hills age model was 101.6% (Aoverall = 71.2%), exceeding the recommended rejection Agreement Index threshold of 60% (89) (Methods). Regardless of the relatively large uncertainty associated with the oldest section of ice (zone 4), the identification of the 130.7 ± 1.8 ky (AICC2012 timescale) Tephra B/Dome Fuji 1,785.14 m (28, 33, 34) within Patriot Hills at −340 m unambiguously demonstrates the presence of Termination II-age ice. Future age constraints will inevitably help improve the accuracy and precision of the age model.

Isotopes.

δD and δ18O isotopic measurements were performed between 1- and 3-m resolution at James Cook University using diffusion sampling–cavity ring-down spectrometry (International Atomic Energy WICO Laboratory ID 16139) (92). This system continuously converts liquid water into water vapor for real-time stable isotope analysis by laser spectroscopy (Picarro L2120-i). See SI Appendix for operating conditions. To ensure reproducibility, a subset of samples was rerun at University of New South Wales ICELAB for δD and δ18O using a Los Gatos Research Liquid Water Isotope Analyzer 24 d (International Atomic Energy WICO Laboratory ID 16117). Reported overall analytical precision on long-term ice core standards is <0.32‰ for δD and <0.13 for δ18O values. All isotopic values are expressed relative to the Vienna Standard Mean Ocean Water 2 (VSMOW2). The isotopic datasets generated in this study are available at the publicly accessible National Oceanic and Atmospheric Administration (NOAA) Paleoclimatology Database (93) and are available upon request.

Ancient DNA analysis.

BIAs offer the opportunity to process large-volume samples of continental Antarctic ice in the field (∼7 kg per temporal sample), creating the prospect of generating sufficient microbial concentrations to permit detailed genetic biodiversity surveys (51, 52) (Fig. 2). To obtain the samples, a Kovac corer was thoroughly cleaned with 1 to 3% bleach and wiped with 95% ethanol between core extractions to minimize cross-contamination. After coring, the top 1 m of ice was removed and discarded, before 1- to 2-m-long cores were collected in 50-cm sections and immediately placed into clean PFTE flexible plastic tubing. A heat sealer was used to close the tubing at the top and bottom of the core. The sealed core was then cut from the remaining tubing with a sterile blade, and the process was repeated to encase the core in a second layer of the plastic tubing for protection during transport. Within 1 to 6 h of extraction, the tubing-encased BIA cores were hung inside a large dome tent to melt via solar radiation over 12 to 24 h, using black plastic bin liners around the plastic tubing to speed up the process where necessary. The melted BIA sample was transferred from the inside layer of tubing directly into a hand-powered vacuum filtration system cleaned with 1 to 3% bleach and ethanol wipes between samples. For each sample, disposable, sterile, 0.45-µm nitrocellulose filters were used to filter and collect whole bacterial organisms trapped in the ice during its formation, and reduce noise caused by environmental DNA. Filters were stored in sterile plastic bags, frozen at −20 °C, and returned to the Australian Centre for Ancient DNA in Adelaide for ultraclean genetic analysis.

Strict ancient DNA methodologies designed to assess low-biomass microbial samples were applied (94) (see SI Appendix for detailed methodology and analysis). DNA from all ice samples as well as extensive sampling and laboratory controls were extracted using two methods to maximize species recovery, and 16S ribosomal RNA libraries were amplified in triplicate using published, universal bacterial and archaeal 16S ribosomal RNA (rRNA) primers. After DNA sequencing, all individually indexed 16S rRNA libraries were de-multiplexed, quality filtered, and imported into QIIME, version 1.8.0. Microbial taxa were identified by comparing sequences to the Geengenes, version 13, reference database and binning sequences with 97% similar to known species into operational taxonomic units using closed reference clustering in UCLUST. Sampling and laboratory contaminants were then filtered from ice samples, and an average of 30.8% of the reads for each sample were retained (SI Appendix, Table S5). Retained sequences were then pooled, and the resulting taxa present in each sample were explored as a proportion of the total filtered DNA sequencing reads. Alpha and beta diversity was explored in QIIME, and importantly, no statistically significant differences in diversity were detected across the samples. Ancient DNA sample data are available upon request. While the current sample numbers limit resolution, our study highlights the untapped potential of BIA genetic data to exploit cryosphere microbial communities to investigate glaciological and environmental change (52).

Ice Sheet Modeling.

To investigate former ice sheet dynamics around the Patriot Hills and across Antarctica, we take a range of values for polar ocean warming (1 to 3 °C) (9, 11, 23) and employ the Parallel Ice Sheet Model (PISM), version 0.6.3 (2), an open source 3D, thermomechanical coupled ice sheet/ice shelf model. PISM employs a stress balance that superposes solutions of the shallow-ice and shallow-shelf equations, and incorporates a pseudoplastic basal substrate rheology to allow for realistic sliding over meltwater saturated sediments, a bed deformation model that simulates mantle dissipation and rebound arising from spatial changes in ice loading through time (95), and a subgrid basal traction and driving stress interpolation scheme to allow realistic grounding-line motion (96, 97). We prescribe a mantle viscosity of 1 × 1020 Pa⋅s, which is lower than the PISM default (1 × 1021) and intended to capture more accurately the weaker mantle of West Antarctica, where the majority of mass loss takes place. In the experiments presented here, we chose not to implement the subgrid scale interpolated ice shelf basal melt component of this scheme (2, 98). Calving is parameterized using horizontal strain rates and a minimum thickness criterion (220 m) (99, 100). Our experimental methodology is identical to that described in detail elsewhere (101, 102). Climate and ocean temperature perturbations are applied as spatially uniform linear increments added to boundary distributions representing present-day conditions. Linear increases take place between 2,000 and 3,000 model years. The first 2,000 y (no forcing) allow any transient behavior associated with model initialization to take place in the absence of environmental perturbations, whereas the subsequent 1,000 y force the ice sheet to evolve slowly to changes in air and ocean temperature and precipitation. All experiments are run at a spatial resolution of 20 km.

Reconstructed summer SST anomalies relative to present day (the 1998 World Ocean Atlas) (9) were used to inform on a range of warmer air and ocean LIG conditions and applied to a stable modern configuration of the Antarctic Ice Sheet to help interpret the Patriot Hills record (Table 1). A limitation of this approach is that the transient history from the preceding glacial state is not simulated. However, for the response of the ice shelves, this colder prehistory should not be critical, and the experiments as performed are directly relevant for the future of the ice sheets. From these simulations, we extract data from the first 10 ky. The ice sheet modeling outputs support the view that ocean (rather than atmosphere) warming was the primary driver of ice shelf collapse and substantial early LIG mass loss in Patriot Hills and across large parts of Antarctica (SI Appendix, Fig. S18). With a surface ocean warming of 2 °C, our simulations suggest isolation and stagnation of ice in Horseshoe Valley and the loss of the Bungenstock Ice Rise within 400 y of warming (equivalent to 0.8 °C of warming as a result of the linear temperature increase over 2,000 to 3,000 model years) (Figs. 5 and 6) and ultimately restricted ice across the wider WSE (SI Appendix, Fig. S20).

We caution that, for the LIG, subsurface ocean warming is poorly constrained. While recent work has suggested that sea surface warming may propagate to depths important for ice shelves (including embayments) within a few decades (103, 104), proxy SSTs could instead record “bottom-up” warming (i.e., as a consequence of circulatory change) and may underestimate the magnitude of the warming. We however, consider, that warming of +2 °C is likely to be at the upper end of potential LIG warming scenarios (14) and the forcing used here in our simulations to be conservative. Recent work using PISM showed that substantial collapse of WAIS is possible within only a few centuries even under modest warming (105). Those simulations used a much stiffer bed parameterization and were run at 5-km resolution. Other studies have suggested with 5-km resolved PISM simulations that, if mass loss comparable to recent decades is maintained for as little as 60 y, the WAIS could be irrevocably destabilized over subsequent millennia through the collapse in the Amundsen Sea sector (60), overcoming any isostatically driven rebound. On the basis of these comparisons, we can be confident that our inference of substantial mass loss from WAIS under modest ocean/atmosphere warming is not especially dependent on the model used, the way that the bed is parameterized, or the resolution of the simulations. Modeled Antarctic ice sheet contributions to global sea level are provided in Table 1. The ice sheet model data are available upon request.

Acknowledgments

C.S.M.T., C.J.F., M.I.B., A.C., and N.R.G. are supported by their respective Australian Research Council (ARC) and Royal Society of New Zealand fellowships. Fieldwork was undertaken under ARC Linkage Project (LP120200724), supported by Linkage Partner Antarctic Logistics and Expeditions. J.W. and K.W. undertook GPR survey of the Patriot Hills record through the Natural Environment Research Council Project (NE/I027576/1) with logistic field support from the British Antarctic Survey. S.M.D. acknowledges financial support from Coleg Cymraeg Cenedlaethol, the European Research Council, and the Fulbright Commission (259253 and FP7/2007-2013). K.K. was supported by Japan Society for the Promotion of Science and the Ministry of Education, Culture, Sports, Science and Technology’s Grants-in-Aid for Scientific Research (15KK0027 and 17H06320). We thank Dr. Chris Hayward and Dr. Gwydion Jones for electron microprobe assistance; Kathryn Lacey and Gareth James for help with preparing the tephra samples; Drs. Nelia Dunbar, Nels Iverson, and Andrei Kurbatov for discussions on the tephra correlations; CSIRO GASLAB personnel for support of gas analysis; Prof. Bill Sturges and Dr. Sam Allin of the Centre for Ocean and Atmospheric Sciences (University of East Anglia, Norwich, UK) for performing the sulfur hexafluoride analyses; Levke Caesar (Potsdam Institute for Climate Impact Research) for preparing the recent trend in SSTs in Fig. 1; Vicki Taylor (British Ocean Sediment Core Research Facility, Southampton, UK) for assistance with marine core sampling; and Dr. Emilie Capron (British Antarctic Survey) for advice on reconstructing early southern LIG temperatures. We thank CSIRO’s contribution, which was supported in part by the Australian Climate Change Science Program, an Australian Government Initiative. We also acknowledge Johannes Sutter, Torsten Albrecht, and Jonathan Kingslake for advice and data on their model simulations. We also thank the editor and two anonymous reviewers for their insightful comments that helped improve this manuscript.

Footnotes

  • ↵1To whom correspondence may be addressed. Email: c.turney{at}unsw.edu.au.
  • ↵2Deceased February 6, 2018.

  • Author contributions: C.S.M.T., C.J.F., and N.R.G. designed research; C.S.M.T., C.J.F., N.R.G., N.P.M., E.v.S., R.T.J., D.E., M.R., D.P.T., S.M.D., C.B.R., Z.A.T., M.I.B., N.C.M., M.K., J.W., K.W., L.S.W., C.M.R., H.M., P.G.A., A.R., A.M., C.-D.H., M.E.W., C.J.M., J.Y., and A.C. performed research; L.S.W. and A.C. contributed new reagents/analytic tools; C.S.M.T., C.J.F., N.R.G., N.P.M., E.v.S., R.T.J., D.E., M.R., S.M.D., C.B.R., Z.A.T., M.I.B., N.C.M., M.K., J.W., K.W., L.S.W., P.G.A., T.v.O., M.C., A.M., S.R., K.K., C.-D.H., M.E.W., and A.C. analyzed data; and C.S.M.T. wrote the paper with input from all authors.

  • The authors declare no competing interest.

  • This article is a PNAS Direct Submission. J.S. is a guest editor invited by the Editorial Board.

  • Data deposition: The isotopic datasets generated in this study are available at the publicly accessible National Oceanic and Atmospheric Administration (NOAA) Paleoclimatology Database (https://www.ncdc.noaa.gov/paleo/study/28610).

  • This article contains supporting information online at https://www.pnas.org/lookup/suppl/doi:10.1073/pnas.1902469117/-/DCSupplemental.

  • Copyright © 2020 the Author(s). Published by PNAS.

This open access article is distributed under Creative Commons Attribution License 4.0 (CC BY).

References

  1. ↵
    1. C. Ritz et al
    ., Potential sea-level rise from Antarctic ice-sheet instability constrained by observations. Nature 528, 115–118 (2015).
    OpenUrlPubMed
  2. ↵
    1. N. R. Golledge et al
    ., The multi-millennial Antarctic commitment to future sea-level rise. Nature 526, 421–425 (2015).
    OpenUrl
  3. ↵
    1. R. M. DeConto,
    2. D. Pollard
    , Contribution of Antarctica to past and future sea-level rise. Nature 531, 591–597 (2016).
    OpenUrlCrossRefPubMed
  4. ↵
    1. A. Dutton et al
    ., Sea-level rise due to polar ice-sheet mass loss during past warm periods. Science 349, aaa4019 (2015).
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. J. H. Mercer
    , West Antarctic Ice Sheet and CO2 greenhouse effect: A threat of disaster. Nature 271, 321–325 (1978).
    OpenUrlCrossRef
  6. ↵
    1. J. H. Mercer,
    2. C. Emiliani
    , Antarctic ice and interglacial high sea levels. Science 168, 1605–1606 (1970).
    OpenUrlFREE Full Text
  7. ↵
    1. J. H. Mercer
    , Antarctic Ice and Sangamon Sea Level (International Association of Scientific Hydrology Publication, Commission of Snow and Ice, General Assembly of Bern, 1968), vol. 79, pp. 217–225.
    OpenUrl
  8. ↵
    1. NEEM Community Members
    , Eemian interglacial reconstructed from a Greenland folded ice core. Nature 493, 489–494 (2013).
    OpenUrlCrossRefPubMed
  9. ↵
    1. E. Capron,
    2. A. Govin,
    3. R. Feng,
    4. B. L. Otto-Bliesner,
    5. E. W. Wolff
    , Critical evaluation of climate syntheses to benchmark CMIP6/PMIP4 127 ka last interglacial simulations in the high-latitude regions. Quat. Sci. Rev. 168, 137–150 (2017).
    OpenUrl
  10. ↵
    1. K. M. Grant et al
    ., Sea-level variability over five glacial cycles. Nat. Commun. 5, 5076 (2014).
    OpenUrlCrossRefPubMed
  11. ↵
    1. C. S. M. Turney,
    2. R. T. Jones
    , Does the Agulhas Current amplify global temperatures during super-interglacials? J. Quat. Sci. 25, 839–843 (2010).
    OpenUrlCrossRef
  12. ↵
    1. C. Buizert et al
    ., Radiometric 81Kr dating identifies 120,000-year-old ice at Taylor Glacier, Antarctica. Proc. Natl. Acad. Sci. U.S.A. 111, 6876–6881 (2014).
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. E. J. Rohling et al
    ., Differences between the last two glacial maxima and implications for ice-sheet, δ18O, and sea-level reconstructions. Quat. Sci. Rev. 176, 1–28 (2017).
    OpenUrl
  14. ↵
    1. C. J. Fogwill et al
    ., Testing the sensitivity of the East Antarctic Ice Sheet to Southern Ocean dynamics: Past changes and future implications. J. Quaternary Sci. 29, 91–98 (2014).
    OpenUrl
  15. ↵
    1. E. J. Rohling et al
    ., Asynchronous Antarctic and Greenland ice-volume contributions to the last interglacial sea-level highstand. Nat. Commun. 10, 5040 (2019).
    OpenUrl
  16. ↵
    1. E. V. Korotkikh et al
    ., The last interglacial as represented in the glaciochemical record from Mount Moulton Blue Ice Area, West Antarctica. Quat. Sci. Rev. 30, 1940–1947 (2011).
    OpenUrl
  17. ↵
    1. E. J. Steig et al
    ., Influence of West Antarctic Ice Sheet collapse on Antarctic surface climate. Geophys. Res. Lett. 42, 4862–4868 (2015).
    OpenUrl
  18. ↵
    1. J. Kingslake et al
    ., Extensive retreat and re-advance of the West Antarctic Ice Sheet during the Holocene. Nature 558, 430–434 (2018).
    OpenUrlCrossRefPubMed
  19. ↵
    1. B. L. Otto-Bliesner et al
    ., How warm was the last interglacial? New model—data comparisons. Philos. Trans. R. Soc. A 371, 20130097 (2013).
    OpenUrlCrossRefPubMed
  20. ↵
    1. J. Sutter,
    2. P. Gierz,
    3. K. Grosfeld,
    4. M. Thoma,
    5. G. Lohmann
    , Ocean temperature thresholds for last interglacial West Antarctic Ice Sheet collapse. Geophys. Res. Lett. 43, 2675–2682 (2016).
    OpenUrl
  21. ↵
    1. A. A. Robel,
    2. A. F. Banwell
    , A speed limit on ice shelf collapse through hydrofracture. Geophys. Res. Lett. 46, 12092–12100 (2019).
    OpenUrl
  22. ↵
    1. F. Clerc,
    2. B. M. Minchew,
    3. M. D. Behn
    , Marine ice cliff instability mitigated by slow removal of ice shelves. Geophys. Res. Lett. 46, 12108–12116 (2019).
    OpenUrl
  23. ↵
    1. J. S. Hoffman,
    2. P. U. Clark,
    3. A. C. Parnell,
    4. F. He
    , Regional and global sea-surface temperatures during the last interglaciation. Science 355, 276–279 (2017).
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. T. Albrecht,
    2. R. Winkelmann,
    3. A. Levermann
    , Glacial cycles simulation of the Antarctic Ice Sheet with PISM—Part 1: Boundary conditions and climatic forcing. Cryosphere Discuss. 2019, 1–56 (2019).
    OpenUrl
  25. ↵
    1. J. Sutter et al
    ., Modelling the Antarctic Ice Sheet across the mid-Pleistocene transition—implications for oldest ice. Cryosphere 13, 2023–2041 (2019).
    OpenUrl
  26. ↵
    1. C. J. Fogwill et al
    ., Drivers of abrupt Holocene shifts in West Antarctic ice stream direction determined from combined ice sheet modelling and geologic signatures. Antarct. Sci. 26, 674–686 (2014).
    OpenUrl
  27. ↵
    1. P. Köhler,
    2. C. Nehrbass-Ahles,
    3. J. Schmitt,
    4. T. F. Stocker,
    5. H. Fischer
    , A 156 kyr smoothed history of the atmospheric greenhouse gases CO2, CH4, and N2O and their radiative forcing. Earth Syst. Sci. Data 9, 363–387 (2017).
    OpenUrl
  28. ↵
    1. D. Veres et al
    ., The Antarctic ice core chronology (AICC2012): An optimized multi-parameter and multi-site dating approach for the last 120 thousand years. Clim. Past 9, 1733–1748 (2013).
    OpenUrlCrossRef
  29. ↵
    1. K. Winter et al
    ., Airborne radar evidence for tributary flow switching in Institute Ice Stream, West Antarctica: Implications for ice sheet configuration and dynamics. J. Geophys. Res. Earth Surf. 120, 1611–1625 (2015).
    OpenUrl
  30. ↵
    1. J. L. Wadham et al
    ., Potential methane reservoirs beneath Antarctica. Nature 488, 633–637 (2012).
    OpenUrlCrossRefPubMed
  31. ↵
    1. C. Turney et al
    ., Late Pleistocene and early Holocene change in the Weddell Sea: A new climate record from the Patriot Hills, Ellsworth Mountains, West Antarctica. J. Quat. Sci. 28, 697–704 (2013).
    OpenUrl
  32. ↵
    1. S. L. Bradley,
    2. M. Siddall,
    3. G. A. Milne,
    4. V. Masson-Delmotte,
    5. E. Wolff
    , Where might we find evidence of a Last Interglacial West Antarctic Ice sheet collapse in Antarctic ice core records? Global Planet. Change 88–89, 64–75 (2012).
    OpenUrl
  33. ↵
    1. R. Uemura et al
    ., Asynchrony between Antarctic temperature and CO2 associated with obliquity over the past 720,000 years. Nat. Commun. 9, 961 (2018).
    OpenUrl
  34. ↵
    1. C.-D. Hillenbrand et al
    ., Volcanic time-markers for marine isotopic stages 6 and 5 in Southern Ocean sediments and Antarctic ice cores: Implications for tephra correlations between palaeoclimatic records. Quat. Sci. Rev. 27, 518–540 (2008).
    OpenUrl
  35. ↵
    1. N. Ross et al
    ., Steep reverse bed slope at the grounding line of the Weddell Sea sector in West Antarctica. Nat. Geosci. 5, 393–396 (2012).
    OpenUrlCrossRef
  36. ↵
    1. A. S. Hein et al
    ., Evidence for the stability of the West Antarctic Ice Sheet divide for 1.4 million years. Nat. Commun. 7, 10325 (2016).
    OpenUrl
  37. ↵
    1. K. Winter et al
    ., Assessing the continuity of the blue ice climate record at Patriot Hills, Horseshoe Valley, West Antarctica. Geophys. Res. Lett. 43, 2019–2026 (2016).
    OpenUrl
  38. ↵
    1. C. J. Fogwill et al
    ., Antarctic ice sheet discharge driven by atmosphere-ocean feedbacks at the last glacial termination. Sci. Rep. 7, 39979 (2017).
    OpenUrl
  39. ↵
    1. M. Siegert,
    2. N. Ross,
    3. H. Corr,
    4. J. Kingslake,
    5. R. Hindmarsh
    , Late Holocene ice-flow reconfiguration in the Weddell Sea sector of West Antarctica. Quat. Sci. Rev. 78, 98–107 (2013).
    OpenUrl
  40. ↵
    1. C. A. Shuman,
    2. E. Berthier,
    3. T. A. Scambos
    , 2001–2009 elevation and mass losses in the Larsen A and B embayments, Antarctic Peninsula. J. Glaciol. 57, 737–754 (2011).
    OpenUrlCrossRef
  41. ↵
    1. H. Rott et al
    ., Changing pattern of ice flow and mass balance for glaciers discharging into the Larsen A and B embayments, Antarctic Peninsula, 2011 to 2016. Cryosphere 12, 1273–1291 (2018).
    OpenUrl
  42. ↵
    1. E. Böhm et al
    ., Strong and deep Atlantic meridional overturning circulation during the last glacial cycle. Nature 517, 73–76 (2015).
    OpenUrlCrossRefPubMed
  43. ↵
    1. G. Marino et al
    ., Bipolar seesaw control on last interglacial sea level. Nature 522, 197–201 (2015).
    OpenUrlCrossRefPubMed
  44. ↵
    1. J. B. Pedro et al
    ., Beyond the bipolar seesaw: Toward a process understanding of interhemispheric coupling. Quat. Sci. Rev. 192, 27–46 (2018).
    OpenUrl
  45. ↵
    1. S. L. Jaccard et al
    ., Two modes of change in Southern Ocean productivity over the past million years. Science 339, 1419–1423 (2013).
    OpenUrlAbstract/FREE Full Text
  46. ↵
    1. C. T. Hayes et al
    ., A stagnation event in the deep South Atlantic during the Last Interglacial period. Science 346, 1514–1517 (2014).
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. S. R. Rintoul
    , The global influence of localized dynamics in the Southern Ocean. Nature 558, 209–218 (2018).
    OpenUrl
  48. ↵
    1. R. T. Jones et al
    ., Delayed maximum northern European summer temperatures during the Last Interglacial as a result of Greenland Ice Sheet melt. Geology 45, 23–26 (2017).
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. E. L. Deaney,
    2. S. Barker,
    3. T. van de Flierdt
    , Timing and nature of AMOC recovery across Termination 2 and magnitude of deglacial CO2 change. Nat. Commun. 8, 14595 (2017).
    OpenUrl
  50. ↵
    1. H. H. Hellmer,
    2. F. Kauker,
    3. R. Timmermann,
    4. J. Determann,
    5. J. Rae
    , Twenty-first-century warming of a large Antarctic ice-shelf cavity by a redirected coastal current. Nature 485, 225–228 (2012).
    OpenUrlCrossRefPubMed
  51. ↵
    1. P. A. Santibáñez et al
    ., Prokaryotes in the WAIS Divide ice core reflect source and transport changes between Last Glacial Maximum and the early Holocene. Glob. Change Biol. 24, 2182–2197 (2018).
    OpenUrl
  52. ↵
    1. A. Boetius,
    2. A. M. Anesio,
    3. J. W. Deming,
    4. J. A. Mikucki,
    5. J. Z. Rapp
    , Microbial ecology of the cryosphere: Sea ice and glacial habitats. Nat. Rev. Microbiol. 13, 677–690 (2015).
    OpenUrlCrossRefPubMed
  53. ↵
    1. L. Chistoserdova,
    2. M. G. Kalyuzhnaya,
    3. M. E. Lidstrom
    , The expanding world of methylotrophic metabolism. Annu. Rev. Microbiol. 63, 477–499 (2009).
    OpenUrlCrossRefPubMed
  54. ↵
    1. M. G. Kalyuzhnaya et al
    ., Methyloversatilis universalis gen. nov., sp. nov., a novel taxon within the Betaproteobacteria represented by three methylotrophic isolates. Int. J. Syst. Evol. Microbiol. 56, 2517–2522 (2006).
    OpenUrlCrossRefPubMed
  55. ↵
    1. G. H. Gudmundsson,
    2. F. S. Paolo,
    3. S. Adusumilli,
    4. H. A. Fricker
    , Instantaneous Antarctic ice-sheet mass loss driven by thinning ice shelves. Geophys. Res. Lett. 46, 13903–13909 (2019).
    OpenUrl
  56. ↵
    1. K. Kawamura et al
    ., Northern Hemisphere forcing of climatic cycles in Antarctica over the past 360,000 years. Nature 448, 912–916 (2007).
    OpenUrlCrossRefPubMed
  57. ↵
    1. R. Schneider,
    2. J. Schmitt,
    3. P. Köhler,
    4. F. Joos,
    5. H. Fischer
    , A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the last glacial inception. Clim. Past 9, 2507–2523 (2013).
    OpenUrlCrossRef
  58. ↵
    1. S. Eggleston,
    2. J. Schmitt,
    3. B. Bereiter,
    4. R. Schneider,
    5. H. Fischer
    , Evolution of the stable carbon isotope composition of atmospheric CO2 over the last glacial cycle. Paleoceanography 31, 434–452 (2016).
    OpenUrlCrossRef
  59. ↵
    1. M. Mengel,
    2. J. Feldmann,
    3. A. Levermann
    , Linear sea-level response to abrupt ocean warming of major West Antarctic ice basin. Nat. Clim. Chang. 6, 71–74 (2015).
    OpenUrl
  60. ↵
    1. J. Feldmann,
    2. A. Levermann
    , Collapse of the West Antarctic Ice Sheet after local destabilization of the Amundsen Basin. Proc. Natl. Acad. Sci. U.S.A. 112, 14191–14196 (2015).
    OpenUrlAbstract/FREE Full Text
  61. ↵
    1. C. Schoof
    , Ice sheet grounding line dynamics: Steady states, stability, and hysteresis. J. Geophys. Res. Earth Surf. 112, F03S28 (2007).
    OpenUrl
  62. ↵
    1. A. Dutton,
    2. J. M. Webster,
    3. D. Zwartz,
    4. K. Lambeck,
    5. B. Wohlfarth
    , Tropical tales of polar ice: Evidence of Last Interglacial polar ice sheet retreat recorded by fossil reefs of the granitic Seychelles Islands. Quat. Sci. Rev. 107, 182–196 (2015).
    OpenUrlCrossRef
  63. ↵
    1. A. Shepherd,
    2. H. A. Fricker,
    3. S. L. Farrell
    , Trends and connections across the Antarctic cryosphere. Nature 558, 223–232 (2018).
    OpenUrl
  64. ↵
    1. L. Caesar,
    2. S. Rahmstorf,
    3. A. Robinson,
    4. G. Feulner,
    5. V. Saba
    , Observed fingerprint of a weakening Atlantic Ocean overturning circulation. Nature 556, 191–196 (2018).
    OpenUrl
  65. ↵
    1. E. Darelius,
    2. I. Fer,
    3. K. W. Nicholls
    , Observed vulnerability of Filchner-Ronne Ice Shelf to wind-driven inflow of warm deep water. Nat. Commun. 7, 12300 (2016).
    OpenUrl
  66. ↵
    1. R. G. Bingham et al
    ., Ice-flow structure and ice dynamic changes in the Weddell Sea sector of West Antarctica from radar-imaged internal layering. J. Geophys. Res. Earth Surf. 120, 655–670 (2015).
    OpenUrl
  67. ↵
    1. E. Rignot,
    2. J. Mouginot,
    3. B. Scheuchl
    , MEaSUREs InSAR-Based Antarctica Ice Velocity Map (Version 2, NASA National Snow and Ice Data Center Distributed Active Archive Center, Boulder, CO, 2017).
  68. ↵
    1. G. Casassa et al
    ., Elevation change and ice flow at Horseshoe Valley, Patriot Hills, West Antarctica. Ann. Glaciol. 39, 20–28 (2004).
    OpenUrl
  69. ↵
    1. H. Jeofry et al
    ., A new bed elevation model for the Weddell Sea sector of the West Antarctic Ice Sheet. Earth Syst. Sci. Data 10, 711–725 (2018).
    OpenUrl
  70. ↵
    1. D. M. Etheridge et al
    ., Natural and anthropogenic changes in atmospheric CO2 over the last 1000 years from air in Antarctic ice and firn. J. Geophys. Res. D Atmos. 101, 4115–4128 (1996).
    OpenUrl
  71. ↵
    1. M. Rubino et al
    ., A revised 1000 year atmospheric δ13C-CO2 record from Law Dome and South Pole, Antarctica. J. Geophys. Res. 118, 8482–8499 (2013).
    OpenUrl
  72. ↵
    1. N. Tindale,
    2. N. Derek,
    3. P. Fraser
    1. R. Francey et al
    ., “The CSIRO (Australia) measurement of greenhouse gases in the global atmosphere” in Baseline Atmospheric Program Australia 1999–2000, N. Tindale, N. Derek, P. Fraser, Eds. (Bureau of Meteorology and CSIRO Atmospheric Research, Melbourne, VIC, Australia, 2003), pp. 42–53.
  73. ↵
    1. C. Hayward
    , High spatial resolution electron probe microanalysis of tephras and melt inclusions without beam-induced chemical modification. Holocene 22, 119–125 (2012).
    OpenUrlCrossRef
  74. ↵
    1. B. Narcisi,
    2. J. R. Petit,
    3. B. Delmonte,
    4. I. Basile-Doelsch,
    5. V. Maggi
    , Characteristics and sources of tephra layers in the EPICA-Dome C ice record (East Antarctica): Implications for past atmospheric circulation and ice core stratigraphic correlations. Earth Planet. Sci. Lett. 239, 253–265 (2005).
    OpenUrl
  75. ↵
    1. B. Narcisi,
    2. J. R. Petit,
    3. B. Delmonte
    , Extended East Antarctic ice-core tephrostratigraphy. Quat. Sci. Rev. 29, 21–27 (2010).
    OpenUrl
  76. ↵
    1. B. Narcisi,
    2. J. R. Petit,
    3. M. Tiepolo
    , A volcanic marker (92 ka) for dating deep east Antarctic ice cores. Quat. Sci. Rev. 25, 2682–2687 (2006).
    OpenUrl
  77. ↵
    1. M. Kohno,
    2. Y. Fujii,
    3. T. Hirata
    , Chemical composition of volcanic glasses in visible tephra layers found in a 2503 m deep ice core from Dome Fuji, Antarctica. Ann. Glaciol. 39, 576–584 (2004).
    OpenUrl
  78. ↵
    1. N. W. Dunbar,
    2. W. C. McIntosh,
    3. R. P. Esser
    , Physical setting and tephrochronology of the summit caldera ice record at Mount Moulton, West Antarctica. Geol. Soc. Am. Bull. 120, 796–812 (2008).
    OpenUrlAbstract/FREE Full Text
  79. ↵
    1. N. W. Dunbar,
    2. A. V. Kurbatov
    , Tephrochronology of the Siple Dome ice core, West Antarctica: Correlations and sources. Quat. Sci. Rev. 30, 1602–1614 (2011).
    OpenUrl
  80. ↵
    1. B. Narcisi,
    2. J. R. Petit,
    3. A. Langone,
    4. B. Stenni
    , A new Eemian record of Antarctic tephra layers retrieved from the Talos Dome ice core (Northern Victoria Land). Global Planet. Change 137, 69–78 (2016).
    OpenUrl
  81. ↵
    1. B. Narcisi,
    2. J. R. Petit,
    3. B. Delmonte,
    4. C. Scarchilli,
    5. B. Stenni
    , A 16,000-yr tephra framework for the Antarctic ice sheet: A contribution from the new Talos Dome core. Quat. Sci. Rev. 49, 52–63 (2012).
    OpenUrl
  82. ↵
    1. B. Narcisi,
    2. J. R. Petit,
    3. J. Chappellaz
    , A 70 ka record of explosive eruptions from the TALDICE ice core (Talos Dome, East Antarctic plateau). J. Quat. Sci. 25, 844–849 (2010).
    OpenUrl
  83. ↵
    1. I. Basile,
    2. J. R. Petit,
    3. S. Touron,
    4. F. E. Grousset,
    5. N. Barkov
    , Volcanic layers in Antarctic (Vostok) ice cores: Source identification and atmospheric implications. J. Geophys. Res. D Atmos. 106, 31915–31931 (2001).
    OpenUrl
  84. ↵
    1. N. A. Iverson et al
    ., The first physical evidence of subglacial volcanism under the West Antarctic Ice Sheet. Sci. Rep. 7, 11457 (2017).
    OpenUrl
  85. ↵
    1. S. Fujita,
    2. F. Parrenin,
    3. M. Severi,
    4. H. Motoyama,
    5. E. Wolff
    , Volcanic synchronization of Dome Fuji and Dome C Antarctic deep ice cores over the past 216 kyr. Clim. Past 11, 1395–1416 (2015).
    OpenUrl
  86. ↵
    1. E. Tomlinson,
    2. T. Thordarson,
    3. W. Müller,
    4. M. Thirlwall,
    5. M. Menzies
    , Microanalysis of tephra by LA-ICP-MS—strategies, advantages and limitations assessed using the Thorsmörk ignimbrite (Southern Iceland). Chem. Geol. 279, 73–89 (2010).
    OpenUrlCrossRef
  87. ↵
    1. K. P. Jochum et al
    ., MPI-DING reference glasses for in situ microanalysis: New reference values for element concentrations and isotope ratios. Geochem. Geophys. Geosyst. 7, Q02008 (2006).
    OpenUrl
  88. ↵
    1. C. Bronk Ramsey,
    2. S. Lee
    , Recent and planned developments of the program OxCal. Radiocarbon 55, 720–730 (2013).
    OpenUrlCrossRef
  89. ↵
    1. C. Bronk Ramsey
    , Dealing with outliers and offsets in radiocarbon dating. Radiocarbon 51, 1023–1045 (2009).
    OpenUrlCrossRef
  90. ↵
    1. T. Higham et al
    ., The timing and spatiotemporal patterning of Neanderthal disappearance. Nature 512, 306–309 (2014).
    OpenUrlCrossRefPubMed
  91. ↵
    1. C. Bronk Ramsey
    , Probability and dating. Radiocarbon 40, 461–474 (1998).
    OpenUrl
  92. ↵
    1. N. C. Munksgaard,
    2. C. M. Wurster,
    3. M. I. Bird
    , Continuous analysis of δ18O and δD values of water by diffusion sampling cavity ring-down spectrometry: A novel sampling device for unattended field monitoring of precipitation, ground and surface waters. Rapid Commun. Mass Spectrom. 25, 3706–3712 (2011).
    OpenUrlPubMed
  93. ↵
    1. C. S. M. Turney,
    2. C. J. Fogwill
    , Patriot Hills, West Antarctica 135,000 year blue ice isotope data. National Oceanic and Atmospheric Administration National Centers for Environmental Information Paleoclimate Database. https://www.ncdc.noaa.gov/paleo-search/study/28610. Deposited 24 December 2018.
  94. ↵
    1. C. J. Adler et al
    ., Sequencing ancient calcified dental plaque shows changes in oral microbiota with dietary shifts of the Neolithic and Industrial revolutions. Nat. Genet. 45, 450–455 (2013).
    OpenUrlCrossRefPubMed
  95. ↵
    1. E. Bueler,
    2. C. S. Lingle,
    3. J. Brown
    , Fast computation of a viscoelastic deformable Earth model for ice-sheet simulations. Ann. Glaciol. 46, 97–105 (2007).
    OpenUrl
  96. ↵
    1. J. Feldmann,
    2. T. Albrecht,
    3. C. Khroulev,
    4. F. Pattyn,
    5. A. Levermann
    , Resolution-dependent performance of grounding line motion in a shallow model compared with a full-Stokes model according to the MISMIP3d intercomparison. J. Glaciol. 60, 353–360 (2014).
    OpenUrlCrossRef
  97. ↵
    1. E. Bueler,
    2. J. Brown
    , Shallow shelf approximation as a “sliding law” in a thermomechanically coupled ice sheet model. J. Geophys. Res. 114, F03008 (2009).
    OpenUrl
  98. ↵
    1. M. A. Martin,
    2. A. Levermann,
    3. R. Winkelmann
    , Comparing ice discharge through West Antarctic gateways: Weddell vs. Amundsen Sea warming. Cryosphere Discuss. 9, 1705–1733 (2015).
    OpenUrl
  99. ↵
    1. T. Albrecht,
    2. A. Levermann
    , Fracture field for large-scale ice dynamics. J. Glaciol. 58, 165–176 (2012).
    OpenUrlCrossRef
  100. ↵
    1. A. Levermann et al
    ., Kinematic first-order calving law implies potential for abrupt ice-shelf retreat. Cryosphere 6, 273–286 (2012).
    OpenUrl
  101. ↵
    1. A. R. A. Aitken et al
    ., Repeated large-scale retreat and advance of Totten Glacier indicated by inland bed erosion. Nature 533, 385–389 (2016).
    OpenUrl
  102. ↵
    1. N. R. Golledge,
    2. R. H. Levy,
    3. R. M. McKay,
    4. T. R. Naish
    , East Antarctic Ice Sheet most vulnerable to Weddell Sea warming. Geophys. Res. Lett. 44, 2343–2351 (2017).
    OpenUrl
  103. ↵
    1. D. Roemmich et al
    ., Unabated planetary warming and its ocean structure since 2006. Nat. Clim. Chang. 5, 240–245 (2015).
    OpenUrl
  104. ↵
    1. W. Llovel,
    2. I. Fukumori,
    3. B. Meyssignac
    , Depth-dependent temperature change contributions to global mean thermosteric sea level rise from 1960 to 2010. Global Planet. Change 101, 113–118 (2013).
    OpenUrl
  105. ↵
    1. N. R. Golledge et al
    ., Global environmental consequences of twenty-first-century ice-sheet melt. Nature 566, 65–72 (2019).
    OpenUrl
  106. ↵
    1. S. J. Johnsen et al
    ., Oxygen isotope and palaeotemperature records from six Greenland ice-core stations: Camp Century, Dye-3, GRIP, GISP2, Renland and NorthGRIP. J. Quat. Sci. 16, 299–307 (2001).
    OpenUrlCrossRef
  107. ↵
    1. K. K. Andersen et al
    ., The Greenland Ice Core Chronology 2005, 15–42 ka. Part 1: Constructing the time scale. Quat. Sci. Rev. 25, 3246–3257 (2006).
    OpenUrlCrossRef
  108. ↵
    1. S.-s. Sun,
    2. W. F. McDonough
    , Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. Geol. Soc. Lond. Spec. Publ. 42, 313–345 (1989).
    OpenUrlCrossRef
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Early Last Interglacial ocean warming drove substantial ice mass loss from Antarctica
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Early Last Interglacial ocean warming drove substantial ice mass loss from Antarctica
Chris S. M. Turney, Christopher J. Fogwill, Nicholas R. Golledge, Nicholas P. McKay, Erik van Sebille, Richard T. Jones, David Etheridge, Mauro Rubino, David P. Thornton, Siwan M. Davies, Christopher Bronk Ramsey, Zoë A. Thomas, Michael I. Bird, Niels C. Munksgaard, Mika Kohno, John Woodward, Kate Winter, Laura S. Weyrich, Camilla M. Rootes, Helen Millman, Paul G. Albert, Andres Rivera, Tas van Ommen, Mark Curran, Andrew Moy, Stefan Rahmstorf, Kenji Kawamura, Claus-Dieter Hillenbrand, Michael E. Weber, Christina J. Manning, Jennifer Young, Alan Cooper
Proceedings of the National Academy of Sciences Feb 2020, 117 (8) 3996-4006; DOI: 10.1073/pnas.1902469117

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Early Last Interglacial ocean warming drove substantial ice mass loss from Antarctica
Chris S. M. Turney, Christopher J. Fogwill, Nicholas R. Golledge, Nicholas P. McKay, Erik van Sebille, Richard T. Jones, David Etheridge, Mauro Rubino, David P. Thornton, Siwan M. Davies, Christopher Bronk Ramsey, Zoë A. Thomas, Michael I. Bird, Niels C. Munksgaard, Mika Kohno, John Woodward, Kate Winter, Laura S. Weyrich, Camilla M. Rootes, Helen Millman, Paul G. Albert, Andres Rivera, Tas van Ommen, Mark Curran, Andrew Moy, Stefan Rahmstorf, Kenji Kawamura, Claus-Dieter Hillenbrand, Michael E. Weber, Christina J. Manning, Jennifer Young, Alan Cooper
Proceedings of the National Academy of Sciences Feb 2020, 117 (8) 3996-4006; DOI: 10.1073/pnas.1902469117
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

Article Classifications

  • Physical Sciences
  • Earth, Atmospheric, and Planetary Sciences
Proceedings of the National Academy of Sciences: 117 (8)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • The Patriot Hills Record
    • Ocean Warming
    • Antarctic Ice Sheet Modeling
    • Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Model of the Amazon forest
News Feature: A sea in the Amazon
Did the Caribbean sweep into the western Amazon millions of years ago, shaping the region’s rich biodiversity?
Image credit: Tacio Cordeiro Bicudo (University of São Paulo, São Paulo, Brazil), Victor Sacek (University of São Paulo, São Paulo, Brazil), and Lucy Reading-Ikkanda (artist).
Syrian archaeological site
Journal Club: In Mesopotamia, early cities may have faltered before climate-driven collapse
Settlements 4,200 years ago may have suffered from overpopulation before drought and lower temperatures ultimately made them unsustainable.
Image credit: Andrea Ricci.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Birds nestling on tree branches
Parent–offspring conflict in songbird fledging
Some songbird parents might improve their own fitness by manipulating their offspring into leaving the nest early, at the cost of fledgling survival, a study finds.
Image credit: Gil Eckrich (photographer).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490