Robert M. DeConto

Rapid Cenozoic glaciation of Antarctica induced by declining atmospheric CO2.

The sudden, widespread glaciation of Antarctica and the associated shift towards colder temperatures at the Eocene/Oligocene boundary (∼34 million years ago) (refs 1–4) is one of the most fundamental reorganizations of global climate known in the geologic record. The glaciation of Antarctica has hitherto been thought to result from the tectonic opening of Southern Ocean gateways, which enabled the formation of the Antarctic Circumpolar Current and the subsequent thermal isolation of the Antarctic continent. Here we simulate the glacial inception and early growth of the East Antarctic Ice Sheet using a general circulation model with coupled components for atmosphere, ocean, ice sheet and sediment, and which incorporates palaeogeography, greenhouse gas, changing orbital parameters, and varying ocean heat transport. In our model, declining Cenozoic CO2 first leads to the formation of small, highly dynamic ice caps on high Antarctic plateaux. At a later time, a CO2 threshold is crossed, initiating ice-sheet height/mass-balance feedbacks that cause the ice caps to expand rapidly with large orbital variations, eventually coalescing into a continental-scale East Antarctic Ice Sheet. According to our simulation the opening of Southern Ocean gateways plays a secondary role in this transition, relative to CO2 concentration.

Modelling West Antarctic ice sheet growth and collapse through the past five million years

The West Antarctic ice sheet (WAIS), with ice volume equivalent to ∼5 m of sea level, has long been considered capable of past and future catastrophic collapse. Today, the ice sheet is fringed by vulnerable floating ice shelves that buttress the fast flow of inland ice streams. Grounding lines are several hundred metres below sea level and the bed deepens upstream, raising the prospect of runaway retreat. Projections of future WAIS behaviour have been hampered by limited understanding of past variations and their underlying forcing mechanisms. Its variation since the Last Glacial Maximum is best known, with grounding lines advancing to the continental-shelf edges around ∼15 kyr ago before retreating to near-modern locations by ∼3 kyr ago. Prior collapses during the warmth of the early Pliocene epoch and some Pleistocene interglacials have been suggested indirectly from records of sea level and deep-sea-core isotopes, and by the discovery of open-ocean diatoms in subglacial sediments. Until now, however, little direct evidence of such behaviour has been available. Here we use a combined ice sheet/ice shelf model capable of high-resolution nesting with a new treatment of grounding-line dynamics and ice-shelf buttressing to simulate Antarctic ice sheet variations over the past five million years. Modelled WAIS variations range from full glacial extents with grounding lines near the continental shelf break, intermediate states similar to modern, and brief but dramatic retreats, leaving only small, isolated ice caps on West Antarctic islands. Transitions between glacial, intermediate and collapsed states are relatively rapid, taking one to several thousand years. Our simulation is in good agreement with a new sediment record (ANDRILL AND-1B) recovered from the western Ross Sea, indicating a long-term trend from more frequently collapsed to more glaciated states, dominant 40-kyr cyclicity in the Pliocene, and major retreats at marine isotope stage 31 (∼1.07 Myr ago) and other super-interglacials.

Contribution of Antarctica to past and future sea-level rise

Polar temperatures over the last several million years have, at times, been slightly warmer than today, yet global mean sea level has been 6–9 metres higher as recently as the Last Interglacial (130,000 to 115,000 years ago) and possibly higher during the Pliocene epoch (about three million years ago). In both cases the Antarctic ice sheet has been implicated as the primary contributor, hinting at its future vulnerability. Here we use a model coupling ice sheet and climate dynamics—including previously underappreciated processes linking atmospheric warming with hydrofracturing of buttressing ice shelves and structural collapse of marine-terminating ice cliffs—that is calibrated against Pliocene and Last Interglacial sea-level estimates and applied to future greenhouse gas emission scenarios. Antarctica has the potential to contribute more than a metre of sea-level rise by 2100 and more than 15 metres by 2500, if emissions continue unabated. In this case atmospheric warming will soon become the dominant driver of ice loss, but prolonged ocean warming will delay its recovery for thousands of years.

Global cooling during the eocene-oligocene climate transition.

About 34 million years ago, Earths climate shifted from a relatively ice-free world to one with glacial conditions on Antarctica characterized by substantial ice sheets. How Earths temperature changed during this climate transition remains poorly understood, and evidence for Northern Hemisphere polar ice is controversial. Here, we report proxy records of sea surface temperatures from multiple ocean localities and show that the high-latitude temperature decrease was substantial and heterogeneous. High-latitude (45 degrees to 70 degrees in both hemispheres) temperatures before the climate transition were ∼20°C and cooled an average of ∼5°C. Our results, combined with ocean and ice-sheet model simulations and benthic oxygen isotope records, indicate that Northern Hemisphere glaciation was not required to accommodate the magnitude of continental ice growth during this time.

Obliquity-paced Pliocene West Antarctic ice sheet oscillations

Thirty years after oxygen isotope records from microfossils deposited in ocean sediments confirmed the hypothesis that variations in the Earth’s orbital geometry control the ice ages, fundamental questions remain over the response of the Antarctic ice sheets to orbital cycles. Furthermore, an understanding of the behaviour of the marine-based West Antarctic ice sheet (WAIS) during the ‘warmer-than-present’ early-Pliocene epoch (∼5–3 Myr ago) is needed to better constrain the possible range of ice-sheet behaviour in the context of future global warming. Here we present a marine glacial record from the upper 600 m of the AND-1B sediment core recovered from beneath the northwest part of the Ross ice shelf by the ANDRILL programme and demonstrate well-dated, ∼40-kyr cyclic variations in ice-sheet extent linked to cycles in insolation influenced by changes in the Earth’s axial tilt (obliquity) during the Pliocene. Our data provide direct evidence for orbitally induced oscillations in the WAIS, which periodically collapsed, resulting in a switch from grounded ice, or ice shelves, to open waters in the Ross embayment when planetary temperatures were up to ∼3 °C warmer than today and atmospheric CO2 concentration was as high as ∼400 p.p.m.v. (refs 5, 6). The evidence is consistent with a new ice-sheet/ice-shelf model that simulates fluctuations in Antarctic ice volume of up to +7 m in equivalent sea level associated with the loss of the WAIS and up to +3 m in equivalent sea level from the East Antarctic ice sheet, in response to ocean-induced melting paced by obliquity. During interglacial times, diatomaceous sediments indicate high surface-water productivity, minimal summer sea ice and air temperatures above freezing, suggesting an additional influence of surface melt under conditions of elevated CO2.

Thresholds for Cenozoic bipolar glaciation.

The long-standing view of Earth’s Cenozoic glacial history calls for the first continental-scale glaciation of Antarctica in the earliest Oligocene epoch (∼33.6 million years ago), followed by the onset of northern-hemispheric glacial cycles in the late Pliocene epoch, about 31 million years later. The pivotal early Oligocene event is characterized by a rapid shift of 1.5 parts per thousand in deep-sea benthic oxygen-isotope values (Oi-1) within a few hundred thousand years, reflecting a combination of terrestrial ice growth and deep-sea cooling. The apparent absence of contemporaneous cooling in deep-sea Mg/Ca records, however, has been argued to reflect the growth of more ice than can be accommodated on Antarctica; this, combined with new evidence of continental cooling and ice-rafted debris in the Northern Hemisphere during this period, raises the possibility that Oi-1 represents a precursory bipolar glaciation. Here we test this hypothesis using an isotope-capable global climate/ice-sheet model that accommodates both the long-term decline of Cenozoic atmospheric CO2 levels and the effects of orbital forcing. We show that the CO2 threshold below which glaciation occurs in the Northern Hemisphere (∼280 p.p.m.v.) is much lower than that for Antarctica (∼750 p.p.m.v.). Therefore, the growth of ice sheets in the Northern Hemisphere immediately following Antarctic glaciation would have required rapid CO2 drawdown within the Oi-1 timeframe, to levels lower than those estimated by geochemical proxies and carbon-cycle models. Instead of bipolar glaciation, we find that Oi-1 is best explained by Antarctic glaciation alone, combined with deep-sea cooling of up to 4 °C and Antarctic ice that is less isotopically depleted (-30 to -35‰) than previously suggested. Proxy CO2 estimates remain above our model’s northern-hemispheric glaciation threshold of ∼280 p.p.m.v. until ∼25 Myr ago, but have been near or below that level ever since. This implies that episodic northern-hemispheric ice sheets have been possible some 20 million years earlier than currently assumed (although still much later than Oi-1) and could explain some of the variability in Miocene sea-level records.

2.8 Million Years of Arctic Climate Change from Lake El’gygytgyn, NE Russia

Crater Core The high-northern latitudes of the Arctic have an important influence on climate and constitute a region with a unique array of complex feedbacks that make it difficult to understand the workings of its climate. Melles et al. (p. 315, published online 21 June) developed a 2.8-million-year record of Arctic climate, using a sediment core from a lake in northeastern Russia that was formed more than 3.5 million years ago by a meteorite impact. Pronounced glacial episodes began 2.6 million years ago but did not achieve orbital pacing for another 700,000 years. A sediment core from a Russian lake provides a high-latitude climate record where prior terrestrial records have been sparse. The reliability of Arctic climate predictions is currently hampered by insufficient knowledge of natural climate variability in the past. A sediment core from Lake El’gygytgyn in northeastern (NE) Russia provides a continuous, high-resolution record from the Arctic, spanning the past 2.8 million years. This core reveals numerous “super interglacials” during the Quaternary; for marine benthic isotope stages (MIS) 11c and 31, maximum summer temperatures and annual precipitation values are ~4° to 5°C and ~300 millimeters higher than those of MIS 1 and 5e. Climate simulations show that these extreme warm conditions are difficult to explain with greenhouse gas and astronomical forcing alone, implying the importance of amplifying feedbacks and far field influences. The timing of Arctic warming relative to West Antarctic Ice Sheet retreats implies strong interhemispheric climate connectivity.

Past extreme warming events linked to massive carbon release from thawing permafrost

Between about 55.5 and 52 million years ago, Earth experienced a series of sudden and extreme global warming events (hyperthermals) superimposed on a long-term warming trend. The first and largest of these events, the Palaeocene–Eocene Thermal Maximum (PETM), is characterized by a massive input of carbon, ocean acidification and an increase in global temperature of about 5 °C within a few thousand years. Although various explanations for the PETM have been proposed, a satisfactory model that accounts for the source, magnitude and timing of carbon release at the PETM and successive hyperthermals remains elusive. Here we use a new astronomically calibrated cyclostratigraphic record from central Italy to show that the Early Eocene hyperthermals occurred during orbits with a combination of high eccentricity and high obliquity. Corresponding climate–ecosystem–soil simulations accounting for rising concentrations of background greenhouse gases and orbital forcing show that the magnitude and timing of the PETM and subsequent hyperthermals can be explained by the orbitally triggered decomposition of soil organic carbon in circum-Arctic and Antarctic terrestrial permafrost. This massive carbon reservoir had the potential to repeatedly release thousands of petagrams (1015 grams) of carbon to the atmosphere–ocean system, once a long-term warming threshold had been reached just before the PETM. Replenishment of permafrost soil carbon stocks following peak warming probably contributed to the rapid recovery from each event, while providing a sensitive carbon reservoir for the next hyperthermal. As background temperatures continued to rise following the PETM, the areal extent of permafrost steadily declined, resulting in an incrementally smaller available carbon pool and smaller hyperthermals at each successive orbital forcing maximum. A mechanism linking Earth’s orbital properties with release of soil carbon from permafrost provides a unifying model accounting for the salient features of the hyperthermals.

Sea-level rise due to polar ice-sheet mass loss during past warm periods

Warming climate, melting ice, rising seas We know that the sea level will rise as climate warms. Nevertheless, accurate projections of how much sea-level rise will occur are difficult to make based solely on modern observations. Determining how ice sheets and sea level have varied in past warm periods can help us better understand how sensitive ice sheets are to higher temperatures. Dutton et al. review recent interdisciplinary progress in understanding this issue, based on data from four different warm intervals over the past 3 million years. Their synthesis provides a clear picture of the progress we have made and the hurdles that still exist. Science, this issue 10.1126/science.aaa4019 Reconstructing past magnitudes, rates, and sources of sea-level rise can help project what our warmer future may hold. BACKGROUND Although thermal expansion of seawater and melting of mountain glaciers have dominated global mean sea level (GMSL) rise over the last century, mass loss from the Greenland and Antarctic ice sheets is expected to exceed other contributions to GMSL rise under future warming. To better constrain polar ice-sheet response to warmer temperatures, we draw on evidence from interglacial periods in the geologic record that experienced warmer polar temperatures and higher GMSLs than present. Coastal records of sea level from these previous warm periods demonstrate geographic variability because of the influence of several geophysical processes that operate across a range of magnitudes and time scales. Inferring GMSL and ice-volume changes from these reconstructions is nontrivial and generally requires the use of geophysical models. ADVANCES Interdisciplinary studies of geologic archives have ushered in a new era of deciphering magnitudes, rates, and sources of sea-level rise. Advances in our understanding of polar ice-sheet response to warmer climates have been made through an increase in the number and geographic distribution of sea-level reconstructions, better ice-sheet constraints, and the recognition that several geophysical processes cause spatially complex patterns in sea level. In particular, accounting for glacial isostatic processes helps to decipher spatial variability in coastal sea-level records and has reconciled a number of site-specific sea-level reconstructions for warm periods that have occurred within the past several hundred thousand years. This enables us to infer that during recent interglacial periods, small increases in global mean temperature and just a few degrees of polar warming relative to the preindustrial period resulted in ≥6 m of GMSL rise. Mantle-driven dynamic topography introduces large uncertainties on longer time scales, affecting reconstructions for time periods such as the Pliocene (~3 million years ago), when atmospheric CO2 was ~400 parts per million (ppm), similar to that of the present. Both modeling and field evidence suggest that polar ice sheets were smaller during this time period, but because dynamic topography can cause tens of meters of vertical displacement at Earth’s surface on million-year time scales and uncertainty in model predictions of this signal are large, it is currently not possible to make a precise estimate of peak GMSL during the Pliocene. OUTLOOK Our present climate is warming to a level associated with significant polar ice-sheet loss in the past, but a number of challenges remain to further constrain ice-sheet sensitivity to climate change using paleo–sea level records. Improving our understanding of rates of GMSL rise due to polar ice-mass loss is perhaps the most societally relevant information the paleorecord can provide, yet robust estimates of rates of GMSL rise associated with polar ice-sheet retreat and/or collapse remain a weakness in existing sea-level reconstructions. Improving existing magnitudes, rates, and sources of GMSL rise will require a better (global) distribution of sea-level reconstructions with high temporal resolution and precise elevations and should include sites close to present and former ice sheets. Translating such sea-level data into a robust GMSL signal demands integration with geophysical models, which in turn can be tested through improved spatial and temporal sampling of coastal records. Further development is needed to refine estimates of past sea level from geochemical proxies. In particular, paired oxygen isotope and Mg/Ca data are currently unable to provide confident, quantitative estimates of peak sea level during these past warm periods. In some GMSL reconstructions, polar ice-sheet retreat is inferred from the total GMSL budget, but identifying the specific ice-sheet sources is currently hindered by limited field evidence at high latitudes. Given the paucity of such data, emerging geochemical and geophysical techniques show promise for identifying the sectors of the ice sheets that were most vulnerable to collapse in the past and perhaps will be again in the future. Peak global mean temperature, atmospheric CO2, maximum global mean sea level (GMSL), and source(s) of meltwater. Light blue shading indicates uncertainty of GMSL maximum. Red pie charts over Greenland and Antarctica denote fraction (not location) of ice retreat. Interdisciplinary studies of geologic archives have ushered in a new era of deciphering magnitudes, rates, and sources of sea-level rise from polar ice-sheet loss during past warm periods. Accounting for glacial isostatic processes helps to reconcile spatial variability in peak sea level during marine isotope stages 5e and 11, when the global mean reached 6 to 9 meters and 6 to 13 meters higher than present, respectively. Dynamic topography introduces large uncertainties on longer time scales, precluding robust sea-level estimates for intervals such as the Pliocene. Present climate is warming to a level associated with significant polar ice-sheet loss in the past. Here, we outline advances and challenges involved in constraining ice-sheet sensitivity to climate change with use of paleo–sea level records.

The role of carbon dioxide during the onset of Antarctic glaciation

Antarctica glaciation began soon after a large decrease in the concentration of atmospheric carbon dioxide around 35 million years ago. Earth’s modern climate, characterized by polar ice sheets and large equator-to-pole temperature gradients, is rooted in environmental changes that promoted Antarctic glaciation ~33.7 million years ago. Onset of Antarctic glaciation reflects a critical tipping point for Earth’s climate and provides a framework for investigating the role of atmospheric carbon dioxide (CO2) during major climatic change. Previously published records of alkenone-based CO2 from high- and low-latitude ocean localities suggested that CO2 increased during glaciation, in contradiction to theory. Here, we further investigate alkenone records and demonstrate that Antarctic and subantarctic data overestimate atmospheric CO2 levels, biasing long-term trends. Our results show that CO2 declined before and during Antarctic glaciation and support a substantial CO2 decrease as the primary agent forcing Antarctic glaciation, consistent with model-derived CO2 thresholds.