Antarctic ice dynamics amplified by Northern Hemisphere sea-level forcing.

Sea-level rise due to ice loss in the Northern Hemisphere in response to insolation and greenhouse gas forcing is thought to have caused grounding-line retreat of marine-based sectors of the Antarctic Ice Sheet (AIS)1-3. Such interhemispheric sea-level forcing may explain the synchronous evolution of global ice sheets over ice-age cycles. Recent studies that indicate that the AIS experienced substantial millennial-scale variability during and after the last deglaciation4-7 (roughly 20,000 to 9,000 years ago) provide further evidence of this sea-level forcing. However, global sea-level change as a result of mass loss from ice sheets is strongly nonuniform, owing to gravitational, deformational and Earth rotational effects8, suggesting that the response of AIS grounding lines to Northern Hemisphere sea-level forcing is more complicated than previously modelled1,2,6. Here, using an ice-sheet model coupled to a global sea-level model, we show that AIS dynamics are amplified by Northern Hemisphere sea-level forcing. As a result of this interhemispheric interaction, a large or rapid Northern Hemisphere sea-level forcing enhances grounding-line advance and associated mass gain of the AIS during glaciation, and grounding-line retreat and mass loss during deglaciation. Relative to models without these interactions, the inclusion of Northern Hemisphere sea-level forcing in our model increases the volume of the AIS during the Last Glacial Maximum (about 26,000 to 20,000 years ago), triggers an earlier retreat of the grounding line and leads to millennial-scale variability throughout the last deglaciation. These findings are consistent with geologic reconstructions of the extent of the AIS during the Last Glacial Maximum and subsequent ice-sheet retreat, and with relative sea-level change in Antarctica3-7,9,10.

Oceanic forcing of penultimate deglacial and last interglacial sea-level rise.

Sea-level histories during the two most recent deglacial-interglacial intervals show substantial differences1-3 despite both periods undergoing similar changes in global mean temperature4,5 and forcing from greenhouse gases6. Although the last interglaciation (LIG) experienced stronger boreal summer insolation forcing than the present interglaciation7, understanding why LIG global mean sea level may have been six to nine metres higher than today has proven particularly challenging2. Extensive areas of polar ice sheets were grounded below sea level during both glacial and interglacial periods, with grounding lines and fringing ice shelves extending onto continental shelves8. This suggests that oceanic forcing by subsurface warming may also have contributed to ice-sheet loss9-12 analogous to ongoing changes in the Antarctic13,14 and Greenland15 ice sheets. Such forcing would have been especially effective during glacial periods, when the Atlantic Meridional Overturning Circulation (AMOC) experienced large variations on millennial timescales16, with a reduction of the AMOC causing subsurface warming throughout much of the Atlantic basin9,12,17. Here we show that greater subsurface warming induced by the longer period of reduced AMOC during the penultimate deglaciation can explain the more-rapid sea-level rise compared with the last deglaciation. This greater forcing also contributed to excess loss from the Greenland and Antarctic ice sheets during the LIG, causing global mean sea level to rise at least four metres above modern levels. When accounting for the combined influences of penultimate and LIG deglaciation on glacial isostatic adjustment, this excess loss of polar ice during the LIG can explain much of the relative sea level recorded by fossil coral reefs and speleothems at intermediate- and far-field sites.

Centennial-scale Holocene climate variations amplified by Antarctic Ice Sheet discharge.

Proxy-based indicators of past climate change show that current global climate models systematically underestimate Holocene-epoch climate variability on centennial to multi-millennial timescales, with the mismatch increasing for longer periods. Proposed explanations for the discrepancy include ocean-atmosphere coupling that is too weak in models, insufficient energy cascades from smaller to larger spatial and temporal scales, or that global climate models do not consider slow climate feedbacks related to the carbon cycle or interactions between ice sheets and climate. Such interactions, however, are known to have strongly affected centennial- to orbital-scale climate variability during past glaciations, and are likely to be important in future climate change. Here we show that fluctuations in Antarctic Ice Sheet discharge caused by relatively small changes in subsurface ocean temperature can amplify multi-centennial climate variability regionally and globally, suggesting that a dynamic Antarctic Ice Sheet may have driven climate fluctuations during the Holocene. We analysed high-temporal-resolution records of iceberg-rafted debris derived from the Antarctic Ice Sheet, and performed both high-spatial-resolution ice-sheet modelling of the Antarctic Ice Sheet and multi-millennial global climate model simulations. Ice-sheet responses to decadal-scale ocean forcing appear to be less important, possibly indicating that the future response of the Antarctic Ice Sheet will be governed more by long-term anthropogenic warming combined with multi-centennial natural variability than by annual or decadal climate oscillations.

Attributing long-term sea-level rise to Paris Agreement emission pledges.

The main contributors to sea-level rise (oceans, glaciers, and ice sheets) respond to climate change on timescales ranging from decades to millennia. A focus on the 21st century thus fails to provide a complete picture of the consequences of anthropogenic greenhouse gas emissions on future sea-level rise and its long-term impacts. Here we identify the committed global mean sea-level rise until 2300 from historical emissions since 1750 and the currently pledged National Determined Contributions (NDC) under the Paris Agreement until 2030. Our results indicate that greenhouse gas emissions over this 280-y period result in about 1 m of committed global mean sea-level rise by 2300, with the NDC emissions from 2016 to 2030 corresponding to around 20 cm or 1/5 of that commitment. We also find that 26 cm (12 cm) of the projected sea-level-rise commitment in 2300 can be attributed to emissions from the top 5 emitting countries (China, United States of America, European Union, India, and Russia) over the 1991-2030 (2016-2030) period. Our findings demonstrate that global and individual country emissions over the first decades of the 21st century alone will cause substantial long-term sea-level rise.

Asynchronous warming and δ18O evolution of deep Atlantic water masses during the last deglaciation.

The large-scale reorganization of deep ocean circulation in the Atlantic involving changes in North Atlantic Deep Water (NADW) and Antarctic Bottom Water (AABW) played a critical role in regulating hemispheric and global climate during the last deglaciation. However, changes in the relative contributions of NADW and AABW and their properties are poorly constrained by marine records, including δ18O of benthic foraminiferal calcite (δ18Oc). Here, we use an isotope-enabled ocean general circulation model with realistic geometry and forcing conditions to simulate the deglacial water mass and δ18O evolution. Model results suggest that, in response to North Atlantic freshwater forcing during the early phase of the last deglaciation, NADW nearly collapses, while AABW mildly weakens. Rather than reflecting changes in NADW or AABW properties caused by freshwater input as suggested previously, the observed phasing difference of deep δ18Oc likely reflects early warming of the deep northern North Atlantic by ∼1.4 °C, while deep Southern Ocean temperature remains largely unchanged. We propose a thermodynamic mechanism to explain the early warming in the North Atlantic, featuring a strong middepth warming and enhanced downward heat flux via vertical mixing. Our results emphasize that the way that ocean circulation affects heat, a dynamic tracer, is considerably different from how it affects passive tracers, like δ18O, and call for caution when inferring water mass changes from δ18Oc records while assuming uniform changes in deep temperatures.

Northern Hemisphere forcing of Southern Hemisphere climate during the last deglaciation.

According to the Milankovitch theory, changes in summer insolation in the high-latitude Northern Hemisphere caused glacial cycles through their impact on ice-sheet mass balance. Statistical analyses of long climate records supported this theory, but they also posed a substantial challenge by showing that changes in Southern Hemisphere climate were in phase with or led those in the north. Although an orbitally forced Northern Hemisphere signal may have been transmitted to the Southern Hemisphere, insolation forcing can also directly influence local Southern Hemisphere climate, potentially intensified by sea-ice feedback, suggesting that the hemispheres may have responded independently to different aspects of orbital forcing. Signal processing of climate records cannot distinguish between these conditions, however, because the proposed insolation forcings share essentially identical variability. Here we use transient simulations with a coupled atmosphere-ocean general circulation model to identify the impacts of forcing from changes in orbits, atmospheric CO(2) concentration, ice sheets and the Atlantic meridional overturning circulation (AMOC) on hemispheric temperatures during the first half of the last deglaciation (22-14.3 kyr BP). Although based on a single model, our transient simulation with only orbital changes supports the Milankovitch theory in showing that the last deglaciation was initiated by rising insolation during spring and summer in the mid-latitude to high-latitude Northern Hemisphere and by terrestrial snow-albedo feedback. The simulation with all forcings best reproduces the timing and magnitude of surface temperature evolution in the Southern Hemisphere in deglacial proxy records. AMOC changes associated with an orbitally induced retreat of Northern Hemisphere ice sheets is the most plausible explanation for the early Southern Hemisphere deglacial warming and its lead over Northern Hemisphere temperature; the ensuing rise in atmospheric CO(2) concentration provided the critical feedback on global deglaciation.

Global warming preceded by increasing carbon dioxide concentrations during the last deglaciation.

The covariation of carbon dioxide (CO(2)) concentration and temperature in Antarctic ice-core records suggests a close link between CO(2) and climate during the Pleistocene ice ages. The role and relative importance of CO(2) in producing these climate changes remains unclear, however, in part because the ice-core deuterium record reflects local rather than global temperature. Here we construct a record of global surface temperature from 80 proxy records and show that temperature is correlated with and generally lags CO(2) during the last (that is, the most recent) deglaciation. Differences between the respective temperature changes of the Northern Hemisphere and Southern Hemisphere parallel variations in the strength of the Atlantic meridional overturning circulation recorded in marine sediments. These observations, together with transient global climate model simulations, support the conclusion that an antiphased hemispheric temperature response to ocean circulation changes superimposed on globally in-phase warming driven by increasing CO(2) concentrations is an explanation for much of the temperature change at the end of the most recent ice age.

The multimillennial sea-level commitment of global warming.

Global mean sea level has been steadily rising over the last century, is projected to increase by the end of this century, and will continue to rise beyond the year 2100 unless the current global mean temperature trend is reversed. Inertia in the climate and global carbon system, however, causes the global mean temperature to decline slowly even after greenhouse gas emissions have ceased, raising the question of how much sea-level commitment is expected for different levels of global mean temperature increase above preindustrial levels. Although sea-level rise over the last century has been dominated by ocean warming and loss of glaciers, the sensitivity suggested from records of past sea levels indicates important contributions should also be expected from the Greenland and Antarctic Ice Sheets. Uncertainties in the paleo-reconstructions, however, necessitate additional strategies to better constrain the sea-level commitment. Here we combine paleo-evidence with simulations from physical models to estimate the future sea-level commitment on a multimillennial time scale and compute associated regional sea-level patterns. Oceanic thermal expansion and the Antarctic Ice Sheet contribute quasi-linearly, with 0.4 m °C(-1) and 1.2 m °C(-1) of warming, respectively. The saturation of the contribution from glaciers is overcompensated by the nonlinear response of the Greenland Ice Sheet. As a consequence we are committed to a sea-level rise of approximately 2.3 m °C(-1) within the next 2,000 y. Considering the lifetime of anthropogenic greenhouse gases, this imposes the need for fundamental adaptation strategies on multicentennial time scales.

Younger Dryas cooling and the Greenland climate response to CO2.

Greenland ice-core δ(18)O-temperature reconstructions suggest a dramatic cooling during the Younger Dryas (YD; 12.9-11.7 ka), with temperatures being as cold as the earlier Oldest Dryas (OD; 18.0-14.6 ka) despite an approximately 50 ppm rise in atmospheric CO(2). Such YD cooling implies a muted Greenland climate response to atmospheric CO(2), contrary to physical predictions of an enhanced high-latitude response to future increases in CO(2). Here we show that North Atlantic sea surface temperature reconstructions as well as transient climate model simulations suggest that the YD over Greenland should be substantially warmer than the OD by approximately 5 °C in response to increased atmospheric CO(2). Additional experiments with an isotope-enabled model suggest that the apparent YD temperature reconstruction derived from the ice-core δ(18)O record is likely an artifact of an altered temperature-δ(18)O relationship due to changing deglacial atmospheric circulation. Our results thus suggest that Greenland climate was warmer during the YD relative to the OD in response to rising atmospheric CO(2), consistent with sea surface temperature reconstructions and physical predictions, and has a sensitivity approximately twice that found in climate models for current climate due to an enhanced albedo feedback during the last deglaciation.

Global climate evolution during the last deglaciation.

Deciphering the evolution of global climate from the end of the Last Glacial Maximum approximately 19 ka to the early Holocene 11 ka presents an outstanding opportunity for understanding the transient response of Earth's climate system to external and internal forcings. During this interval of global warming, the decay of ice sheets caused global mean sea level to rise by approximately 80 m; terrestrial and marine ecosystems experienced large disturbances and range shifts; perturbations to the carbon cycle resulted in a net release of the greenhouse gases CO(2) and CH(4) to the atmosphere; and changes in atmosphere and ocean circulation affected the global distribution and fluxes of water and heat. Here we summarize a major effort by the paleoclimate research community to characterize these changes through the development of well-dated, high-resolution records of the deep and intermediate ocean as well as surface climate. Our synthesis indicates that the superposition of two modes explains much of the variability in regional and global climate during the last deglaciation, with a strong association between the first mode and variations in greenhouse gases, and between the second mode and variations in the Atlantic meridional overturning circulation.

Ice-shelf collapse from subsurface warming as a trigger for Heinrich events.

Episodic iceberg-discharge events from the Hudson Strait Ice Stream (HSIS) of the Laurentide Ice Sheet, referred to as Heinrich events, are commonly attributed to internal ice-sheet instabilities, but their systematic occurrence at the culmination of a large reduction in the Atlantic meridional overturning circulation (AMOC) indicates a climate control. We report Mg/Ca data on benthic foraminifera from an intermediate-depth site in the northwest Atlantic and results from a climate-model simulation that reveal basin-wide subsurface warming at the same time as large reductions in the AMOC, with temperature increasing by approximately 2 °C over a 1-2 kyr interval prior to a Heinrich event. In simulations with an ocean model coupled to a thermodynamically active ice shelf, the increase in subsurface temperature increases basal melt rate under an ice shelf fronting the HSIS by a factor of approximately 6. By analogy with recent observations in Antarctica, the resulting ice-shelf loss and attendant HSIS acceleration would produce a Heinrich event.

Global and regional temperature change over the past 4.5 million years.

Much of our understanding of Cenozoic climate is based on the record of δ18O measured in benthic foraminifera. However, this measurement reflects a combined signal of global temperature and sea level, thus preventing a clear understanding of the interactions and feedbacks of the climate system in causing global temperature change. Our new reconstruction of temperature change over the past 4.5 million years includes two phases of long-term cooling, with the second phase of accelerated cooling during the Middle Pleistocene Transition (1.5 to 0.9 million years ago) being accompanied by a transition from dominant 41,000-year low-amplitude periodicity to dominant 100,000-year high-amplitude periodicity. Changes in the rates of long-term cooling and variability are consistent with changes in the carbon cycle driven initially by geologic processes, followed by additional changes in the Southern Ocean carbon cycle.

Assessing population exposure to coastal flooding due to sea level rise.

The exposure of populations to sea-level rise (SLR) is a leading indicator assessing the impact of future climate change on coastal regions. SLR exposes coastal populations to a spectrum of impacts with broad spatial and temporal heterogeneity, but exposure assessments often narrowly define the spatial zone of flooding. Here we show how choice of zone results in differential exposure estimates across space and time. Further, we apply a spatio-temporal flood-modeling approach that integrates across these spatial zones to assess the annual probability of population exposure. We apply our model to the coastal United States to demonstrate a more robust assessment of population exposure to flooding from SLR in any given year. Our results suggest that more explicit decisions regarding spatial zone (and associated temporal implication) will improve adaptation planning and policies by indicating the relative chance and magnitude of coastal populations to be affected by future SLR.

Rapid postglacial rebound amplifies global sea level rise following West Antarctic Ice Sheet collapse.

Geodetic, seismic, and geological evidence indicates that West Antarctica is underlain by low-viscosity shallow mantle. Thus, as marine-based sectors of the West Antarctic Ice Sheet (WAIS) retreated during past interglacials, or will retreat in the future, exposed bedrock will rebound rapidly and flux meltwater out into the open ocean. Previous studies have suggested that this contribution to global mean sea level (GMSL) rise is small and occurs slowly. We challenge this notion using sea level predictions that incorporate both the outflux mechanism and complex three-dimensional viscoelastic mantle structure. In the case of the last interglacial, where the GMSL contribution from WAIS collapse is often cited as ~3 to 4 meters, the outflux mechanism contributes ~1 meter of additional GMSL change within ~1 thousand years of the collapse. Using a projection of future WAIS collapse, we also demonstrate that the outflux can substantially amplify GMSL rise estimates over the next century.

Regional and global sea-surface temperatures during the last interglaciation.

The last interglaciation (LIG, 129 to 116 thousand years ago) was the most recent time in Earth's history when global mean sea level was substantially higher than it is at present. However, reconstructions of LIG global temperature remain uncertain, with estimates ranging from no significant difference to nearly 2°C warmer than present-day temperatures. Here we use a network of sea-surface temperature (SST) records to reconstruct spatiotemporal variability in regional and global SSTs during the LIG. Our results indicate that peak LIG global mean annual SSTs were 0.5 ± 0.3°C warmer than the climatological mean from 1870 to 1889 and indistinguishable from the 1995 to 2014 mean. LIG warming in the extratropical latitudes occurred in response to boreal insolation and the bipolar seesaw, whereas tropical SSTs were slightly cooler than the 1870 to 1889 mean in response to reduced mean annual insolation.

Regional and global forcing of glacier retreat during the last deglaciation.

The ongoing retreat of glaciers globally is one of the clearest manifestations of recent global warming associated with rising greenhouse gas concentrations. By comparison, the importance of greenhouse gases in driving glacier retreat during the most recent deglaciation, the last major interval of global warming, is unclear due to uncertainties in the timing of retreat around the world. Here we use recently improved cosmogenic-nuclide production-rate calibrations to recalculate the ages of 1,116 glacial boulders from 195 moraines that provide broad coverage of retreat in mid-to-low-latitude regions. This revised history, in conjunction with transient climate model simulations, suggests that while several regional-scale forcings, including insolation, ice sheets and ocean circulation, modulated glacier responses regionally, they are unable to account for global-scale retreat, which is most likely related to increasing greenhouse gas concentrations.

Coherent changes of southeastern equatorial and northern African rainfall during the last deglaciation.

During the last deglaciation, wetter conditions developed abruptly ~14,700 years ago in southeastern equatorial and northern Africa and continued into the Holocene. Explaining the abrupt onset and hemispheric coherence of this early African Humid Period is challenging due to opposing seasonal insolation patterns. In this work, we use a transient simulation with a climate model that provides a mechanistic understanding of deglacial tropical African precipitation changes. Our results show that meltwater-induced reduction in the Atlantic meridional overturning circulation (AMOC) during the early deglaciation suppressed precipitation in both regions. Once the AMOC reestablished, wetter conditions developed north of the equator in response to high summer insolation and increasing greenhouse gas (GHG) concentrations, whereas wetter conditions south of the equator were a response primarily to the GHG increase.

A reconstruction of regional and global temperature for the past 11,300 years.

Surface temperature reconstructions of the past 1500 years suggest that recent warming is unprecedented in that time. Here we provide a broader perspective by reconstructing regional and global temperature anomalies for the past 11,300 years from 73 globally distributed records. Early Holocene (10,000 to 5000 years ago) warmth is followed by ~0.7°C cooling through the middle to late Holocene (<5000 years ago), culminating in the coolest temperatures of the Holocene during the Little Ice Age, about 200 years ago. This cooling is largely associated with ~2°C change in the North Atlantic. Current global temperatures of the past decade have not yet exceeded peak interglacial values but are warmer than during ~75% of the Holocene temperature history. Intergovernmental Panel on Climate Change model projections for 2100 exceed the full distribution of Holocene temperature under all plausible greenhouse gas emission scenarios.