The Paris Climate Agreement and future sea-level rise from Antarctica.

The Paris Agreement aims to limit global mean warming in the twenty-first century to less than 2 degrees Celsius above preindustrial levels, and to promote further efforts to limit warming to 1.5 degrees Celsius1. The amount of greenhouse gas emissions in coming decades will be consequential for global mean sea level (GMSL) on century and longer timescales through a combination of ocean thermal expansion and loss of land ice2. The Antarctic Ice Sheet (AIS) is Earth's largest land ice reservoir (equivalent to 57.9 metres of GMSL)3, and its ice loss is accelerating4. Extensive regions of the AIS are grounded below sea level and susceptible to dynamical instabilities5-8 that are capable of producing very rapid retreat8. Yet the potential for the implementation of the Paris Agreement temperature targets to slow or stop the onset of these instabilities has not been directly tested with physics-based models. Here we use an observationally calibrated ice sheet-shelf model to show that with global warming limited to 2 degrees Celsius or less, Antarctic ice loss will continue at a pace similar to today's throughout the twenty-first century. However, scenarios more consistent with current policies (allowing 3 degrees Celsius of warming) give an abrupt jump in the pace of Antarctic ice loss after around 2060, contributing about 0.5 centimetres GMSL rise per year by 2100-an order of magnitude faster than today4. More fossil-fuel-intensive scenarios9 result in even greater acceleration. Ice-sheet retreat initiated by the thinning and loss of buttressing ice shelves continues for centuries, regardless of bedrock and sea-level feedback mechanisms10-12 or geoengineered carbon dioxide reduction. These results demonstrate the possibility that rapid and unstoppable sea-level rise from Antarctica will be triggered if Paris Agreement targets are exceeded.

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.

Understanding of Contemporary Regional Sea-Level Change and the Implications for the Future.

Global sea level provides an important indicator of the state of the warming climate, but changes in regional sea level are most relevant for coastal communities around the world. With improvements to the sea-level observing system, the knowledge of regional sea-level change has advanced dramatically in recent years. Satellite measurements coupled with in situ observations have allowed for comprehensive study and improved understanding of the diverse set of drivers that lead to variations in sea level in space and time. Despite the advances, gaps in the understanding of contemporary sea-level change remain and inhibit the ability to predict how the relevant processes may lead to future change. These gaps arise in part due to the complexity of the linkages between the drivers of sea-level change. Here we review the individual processes which lead to sea-level change and then describe how they combine and vary regionally. The intent of the paper is to provide an overview of the current state of understanding of the processes that cause regional sea-level change and to identify and discuss limitations and uncertainty in our understanding of these processes. Areas where the lack of understanding or gaps in knowledge inhibit the ability to provide the needed information for comprehensive planning efforts are of particular focus. Finally, a goal of this paper is to highlight the role of the expanded sea-level observation network-particularly as related to satellite observations-in the improved scientific understanding of the contributors to regional sea-level change.

Impact of climate change on New York City's coastal flood hazard: Increasing flood heights from the preindustrial to 2300 CE.

The flood hazard in New York City depends on both storm surges and rising sea levels. We combine modeled storm surges with probabilistic sea-level rise projections to assess future coastal inundation in New York City from the preindustrial era through 2300 CE. The storm surges are derived from large sets of synthetic tropical cyclones, downscaled from RCP8.5 simulations from three CMIP5 models. The sea-level rise projections account for potential partial collapse of the Antarctic ice sheet in assessing future coastal inundation. CMIP5 models indicate that there will be minimal change in storm-surge heights from 2010 to 2100 or 2300, because the predicted strengthening of the strongest storms will be compensated by storm tracks moving offshore at the latitude of New York City. However, projected sea-level rise causes overall flood heights associated with tropical cyclones in New York City in coming centuries to increase greatly compared with preindustrial or modern flood heights. For the various sea-level rise scenarios we consider, the 1-in-500-y flood event increases from 3.4 m above mean tidal level during 1970-2005 to 4.0-5.1 m above mean tidal level by 2080-2100 and ranges from 5.0-15.4 m above mean tidal level by 2280-2300. Further, we find that the return period of a 2.25-m flood has decreased from ∼500 y before 1800 to ∼25 y during 1970-2005 and further decreases to ∼5 y by 2030-2045 in 95% of our simulations. The 2.25-m flood height is permanently exceeded by 2280-2300 for scenarios that include Antarctica's potential partial collapse.

Dynamic Antarctic ice sheet during the early to mid-Miocene.

Geological data indicate that there were major variations in Antarctic ice sheet volume and extent during the early to mid-Miocene. Simulating such large-scale changes is problematic because of a strong hysteresis effect, which results in stability once the ice sheets have reached continental size. A relatively narrow range of atmospheric CO2 concentrations indicated by proxy records exacerbates this problem. Here, we are able to simulate large-scale variability of the early to mid-Miocene Antarctic ice sheet because of three developments in our modeling approach. (i) We use a climate-ice sheet coupling method utilizing a high-resolution atmospheric component to account for ice sheet-climate feedbacks. (ii) The ice sheet model includes recently proposed mechanisms for retreat into deep subglacial basins caused by ice-cliff failure and ice-shelf hydrofracture. (iii) We account for changes in the oxygen isotopic composition of the ice sheet by using isotope-enabled climate and ice sheet models. We compare our modeling results with ice-proximal records emerging from a sedimentological drill core from the Ross Sea (Andrill-2A) that is presented in a companion article. The variability in Antarctic ice volume that we simulate is equivalent to a seawater oxygen isotope signal of 0.52-0.66‰, or a sea level equivalent change of 30-36 m, for a range of atmospheric CO2 between 280 and 500 ppm and a changing astronomical configuration. This result represents a substantial advance in resolving the long-standing model data conflict of Miocene Antarctic ice sheet and sea level variability.

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 (10(15) 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.

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 approximately 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 approximately 15 kyr ago before retreating to near-modern locations by approximately 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 ( approximately 1.07 Myr ago) and other super-interglacials.

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 ( approximately 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 CO(2) levels and the effects of orbital forcing. We show that the CO(2) threshold below which glaciation occurs in the Northern Hemisphere ( approximately 280 p.p.m.v.) is much lower than that for Antarctica ( approximately 750 p.p.m.v.). Therefore, the growth of ice sheets in the Northern Hemisphere immediately following Antarctic glaciation would have required rapid CO(2) 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 degrees C and Antarctic ice that is less isotopically depleted (-30 to -35 per thousand) than previously suggested. Proxy CO(2) estimates remain above our model's northern-hemispheric glaciation threshold of approximately 280 p.p.m.v. until approximately 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.

Competing climate feedbacks of ice sheet freshwater discharge in a warming world.

Freshwater discharge from ice sheets induces surface atmospheric cooling and subsurface ocean warming, which are associated with negative and positive feedbacks respectively. However, uncertainties persist regarding these feedbacks' relative strength and combined effect. Here we assess associated feedbacks in a coupled ice sheet-climate model, and show that for the Antarctic Ice Sheet the positive feedback dominates in moderate future warming scenarios and in the early stage of ice sheet retreat, but is overwhelmed by the negative feedback in intensive warming scenarios when the West Antarctic Ice Sheet undergoes catastrophic collapse. The Atlantic Meridional Overturning Circulation is affected by freshwater discharge from both the Greenland and the Antarctic ice sheets and, as an interhemispheric teleconnection bridge, exacerbates the opposing ice sheet's retreat via the Bipolar Seesaw. These results highlight the crucial role of ice sheet-climate interactions via freshwater flux in future ice sheet retreat and associated sea-level rise.

Climate model differences contribute deep uncertainty in future Antarctic ice loss.

Future projections of ice sheets in response to different climate scenarios and their associated contributions to sea level changes are subject to deep uncertainty due to ice sheet instability processes, hampering a proper risk assessment of sea level rise and enaction of mitigation/adaptation strategies. For a systematic evaluation of the uncertainty due to climate model fields used as input to the ice sheet models, we drive a three-dimensional model of the Antarctic Ice Sheet (AIS) with the output from 36 climate models to simulate past and future changes in the AIS. Simulations show that a few climate models result in partial collapse of the West AIS under modeled preindustrial climates, and the spread in future changes in the AIS's volume is comparable to the structural uncertainty originating from differing ice sheet models. These results highlight the need for improved representations of physical processes important for polar climate in climate models.

Toward a Cenozoic history of atmospheric CO2.

The geological record encodes the relationship between climate and atmospheric carbon dioxide (CO2) over long and short timescales, as well as potential drivers of evolutionary transitions. However, reconstructing CO2 beyond direct measurements requires the use of paleoproxies and herein lies the challenge, as proxies differ in their assumptions, degree of understanding, and even reconstructed values. In this study, we critically evaluated, categorized, and integrated available proxies to create a high-fidelity and transparently constructed atmospheric CO2 record spanning the past 66 million years. This newly constructed record provides clearer evidence for higher Earth system sensitivity in the past and for the role of CO2 thresholds in biological and cryosphere evolution.

Numerical simulations of a kilometre-thick Arctic ice shelf consistent with ice grounding observations.

Recently obtained geophysical data show sets of parallel erosional features on the Lomonosov Ridge in the central Arctic Basin, indicative of ice grounding in water depths up to 1280 m. These features have been interpreted as being formed by an ice shelf-either restricted to the Amerasian Basin (the "minimum model") or extending across the entire Arctic Basin. Here, we use a numerical ice sheet-shelf model to explore how such an ice shelf could form. We rule out the "minimum model" and suggest that grounding on the Lomonosov Ridge requires complete Arctic ice shelf cover; this places a minimum estimate on its volume, which would have exceeded that of the modern Greenland Ice Sheet. Buttressing provided by an Arctic ice shelf would have increased volumes of the peripheral terrestrial ice sheets. An Arctic ice shelf could have formed even in the absence of a hypothesised East Siberian Ice Sheet.

Windblown Pliocene diatoms and East Antarctic Ice Sheet retreat.

Marine diatoms in tillites along the Transantarctic Mountains (TAMs) have been used to suggest a diminished East Antarctic Ice Sheet (EAIS) during Pliocene warm periods. Updated ice-sheet modelling shows significant Pliocene EAIS retreat, creating marine embayments into the Wilkes and Aurora basins that were conducive to high diatom productivity and rapid accumulation of diatomaceous sediments. Here we show that subsequent isostatic uplift exposed accumulated unconsolidated marine deposits to wind erosion. We report new atmospheric modelling utilizing Pliocene climate and derived Antarctic landscapes indicating that prevailing mid-altitude winds transported diatoms towards the TAMs, dominantly from extensive emerged coastal deposits of the Aurora Basin. This result unifies leading ideas from competing sides of a contentious debate about the origin of the diatoms in the TAMs and their link to EAIS history, supporting the view that parts of the EAIS are vulnerable to relatively modest warming, with possible implications for future sea-level rise.

Antarctic glacio-eustatic contributions to late Miocene Mediterranean desiccation and reflooding.

The Messinian Salinity Crisis (MSC) was a marked late Neogene oceanographic event during which the Mediterranean Sea evaporated. Its causes remain unresolved, with tectonic restrictions to the Atlantic Ocean or glacio-eustatic restriction of flow during sea-level lowstands, or a mixture of the two mechanisms, being proposed. Here we present the first direct geological evidence of Antarctic ice-sheet (AIS) expansion at the MSC onset and use a δ(18)O record to model relative sea-level changes. Antarctic sedimentary successions indicate AIS expansion at 6 Ma coincident with major MSC desiccation; relative sea-level modelling indicates a prolonged ∼50 m lowstand at the Strait of Gibraltar, which resulted from AIS expansion and local evaporation of sea water in concert with evaporite precipitation that caused lithospheric deformation. Our results reconcile MSC events and demonstrate that desiccation and refilling were timed by the interplay between glacio-eustatic sea-level variations, glacial isostatic adjustment and mantle deformation in response to changing water and evaporite loads.

2.8 million years of Arctic climate change from Lake El'gygytgyn, NE Russia.

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.

The role of carbon dioxide during the onset of Antarctic glaciation.

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 (CO(2)) during major climatic change. Previously published records of alkenone-based CO(2) from high- and low-latitude ocean localities suggested that CO(2) increased during glaciation, in contradiction to theory. Here, we further investigate alkenone records and demonstrate that Antarctic and subantarctic data overestimate atmospheric CO(2) levels, biasing long-term trends. Our results show that CO(2) declined before and during Antarctic glaciation and support a substantial CO(2) decrease as the primary agent forcing Antarctic glaciation, consistent with model-derived CO(2) thresholds.

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 (approximately 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.