Highly restricted near-surface permafrost extent during the mid-Pliocene warm period

Accurate understanding of permafrost dynamics is critical for evaluating and mitigating impacts that may arise as permafrost degrades in the future; however, existing projections have large uncertainties. Studies of how permafrost responded historically during Earth's past warm periods are helpful in exploring potential future permafrost behavior and to evaluate the uncertainty of future permafrost change projections. Here, we combine a surface frost index model with outputs from the second phase of the Pliocene Model Intercomparison Project to simulate the near-surface (~3 to 4 m depth) permafrost state in the Northern Hemisphere during the mid-Pliocene warm period (mPWP, ~3.264 to 3.025 Ma). This period shares similarities with the projected future climate. Constrained by proxy-based surface air temperature records, our simulations demonstrate that near-surface permafrost was highly spatially restricted during the mPWP and was 93 ± 3% smaller than the preindustrial extent. Near-surface permafrost was present only in the eastern Siberian uplands, Canadian high Arctic Archipelago, and northernmost Greenland. The simulations are similar to near-surface permafrost changes projected for the end of this century under the SSP5-8.5 scenario and provide a perspective on the potential permafrost behavior that may be expected in a warmer world.

Abrupt climate changes in the last two deglaciations simulated with different Northern ice sheet discharge and insolation

There were significant differences between the last two deglaciations, particularly in Atlantic Meridional Overturning Circulation (AMOC) and Antarctic warming in the deglaciations and the following interglacials. Here, we present transient simulations of deglaciation using a coupled atmosphere–ocean general circulation model for the last two deglaciations focusing on the impact of ice sheet discharge on climate changes associated with the AMOC in the first part, and the sensitivity studies using a Northern Hemisphere ice sheet model in the second part. We show that a set of abrupt climate changes of the last deglaciation, including Bolling–Allerod warming, the Younger Dryas, and onset of the Holocene were simulated with gradual changes of both ice sheet discharge and radiative forcing. On the other hand, penultimate deglaciation, with the abrupt climate change only at the beginning of the last interglacial was simulated when the ice sheet discharge was greater than in the last deglaciation by a factor of 1.5. The results, together with Northern Hemisphere ice sheet model experiments suggest the importance of the transient climate and AMOC responses to the different orbital forcing conditions of the last two deglaciations, through the mechanisms of mass loss of the Northern Hemisphere ice sheet and meltwater influx to the ocean.

Glacial carbon cycle changes by Southern Ocean processes with sedimentary amplification

Recent paleo reconstructions suggest that increased carbon storage in the Southern Ocean during glacial periods contributed to low glacial atmospheric carbon dioxide concentration (pCO₂). However, quantifying its contribution in three-dimensional ocean general circulation models (OGCMs) has proven challenging. Here, we show that OGCM simulation with sedimentary process considering enhanced Southern Ocean salinity stratification and iron fertilization from glaciogenic dust during glacial periods improves model-data agreement of glacial deep water with isotopically light carbon, low oxygen, and old radiocarbon ages. The glacial simulation shows a 77-ppm reduction of atmospheric pCO₂, which closely matches the paleo record. The Southern Ocean salinity stratification and the iron fertilization from glaciogenic dust amplified the carbonate sedimentary feedback, which caused most of the increased carbon storage in the deep ocean and played an important role in pCO₂ reduction. The model-data agreement of Southern Ocean properties is crucial for simulating glacial changes in the ocean carbon cycle.

Antarctic surface temperature and elevation during the Last Glacial Maximum

Water-stable isotopes in polar ice cores are a widely used temperature proxy in paleoclimate reconstruction, yet calibration remains challenging in East Antarctica. Here, we reconstruct the magnitude and spatial pattern of Last Glacial Maximum surface cooling in Antarctica using borehole thermometry and firn properties in seven ice cores. West Antarctic sites cooled ~10°C relative to the preindustrial period. East Antarctic sites show a range from ~4° to ~7°C cooling, which is consistent with the results of global climate models when the effects of topographic changes indicated with ice core air-content data are included, but less than those indicated with the use of water-stable isotopes calibrated against modern spatial gradients. An altered Antarctic temperature inversion during the glacial reconciles our estimates with water-isotope observations.

Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison

Earlier large-scale Greenland ice sheet sea-level projections (e.g., those run during the ice2sea and SeaRISE initiatives) have shown that ice sheet initial conditions have a large effect on the projections and give rise to important uncertainties. The goal of the initMIP-Greenland intercomparison exercise is to compare, evaluate and improve the initialisation techniques used in the ice sheet modelling community and to estimate the associated uncertainties in modelled mass changes. initMIP-Greenland is the first in a series of ice sheet model intercomparison activities within ISMIP6 (the Ice Sheet Model Intercomparison Project for CMIP6), which is the primary activity within the Coupled Model Intercomparison Project - phase 6 (CMIP6) focusing on the ice sheets. Two experiments for the large-scale Greenland ice sheet have been designed to allow intercomparison between participating models of 1) the initial present-day state of the ice sheet and 2) the response in two idealised forward experiments. The forward experiments serve to evaluate the initialisation in terms of model drift (forward run without additional forcing) and in response to a large perturbation (prescribed surface mass balance anomaly), and should not be interpreted as sea-level projections. We present and discuss results that highlight the diversity of data sets, boundary conditions and initialisation techniques used in the community to generate initial states of the Greenland ice sheet. We find good agreement across the ensemble for the dynamic response to surface mass balance changes in areas where the simulated ice sheets overlap, but differences arising from the initial size of the ice sheet. The model drift in the control experiment is reduced for models that participated in earlier intercomparison exercises.

Overestimate of committed warming

Palaeoclimate variations are an essential component in constraining future projections of climate change as a function of increasing anthropogenic greenhouse gases1 . The Earth System Sensitivity (ESS) describes the multi-millennial response of Earth (in terms of global mean temperature) to a doubling of CO2 concentrations. A recent study2 used a correlation of inferred temperatures and radiative forcing from greenhouse gases over the past 800,000 years3 to estimate the ESS from present day CO2 is about 9°C, and to imply a long-term commitment of 3–7°C even if greenhouse gas levels remain at present-day concentrations. However, we demonstrate that the methodology of ref. 2 does not reliably estimate the ESS in the presence of orbital forcing of ice age cycles and therefore conclude that the inferred2 present-day committed warming is considerably overestimated. There is a Reply to this Comment by Snyder, C. W. Nature 547, http://dx.doi.org/10.1038/nature22804 (2017).

State dependence of climatic instability over the past 720,000 years from Antarctic ice cores and climate modeling

Climatic variabilities on millennial and longer time scales with a bipolar seesaw pattern have been documented in paleoclimatic records, but their frequencies, relationships with mean climatic state, and mechanisms remain unclear. Understanding the processes and sensitivities that underlie these changes will underpin better understanding of the climate system and projections of its future change. We investigate the long-term characteristics of climatic variability using a new ice-core record from Dome Fuji, East Antarctica, combined with an existing long record from the Dome C ice core. Antarctic warming events over the past 720,000 years are most frequent when the Antarctic temperature is slightly below average on orbital time scales, equivalent to an intermediate climate during glacial periods, whereas interglacial and fully glaciated climates are unfavourable for a millennial-scale bipolar seesaw. Numerical experiments using a fully coupled atmosphere-ocean general circulation model with freshwater hosing in the northern North Atlantic showed that climate becomes most unstable in intermediate glacial conditions associated with large changes in sea ice and the Atlantic Meridional Overturning Circulation. Model sensitivity experiments suggest that the prerequisite for the most frequent climate instability with bipolar seesaw pattern during the late Pleistocene era is associated with reduced atmospheric CO2 concentration via global cooling and sea ice formation in the North Atlantic, in addition to extended Northern Hemisphere ice sheets.

Ice Sheet Model Intercomparison Project (ISMIP6) contribution to CMIP6

Reducing the uncertainty in the past, present and future contribution of ice sheets to sea-level change requires a coordinated effort between the climate and glaciology communities. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) is the primary activity within the Coupled Model Intercomparison Project - phase 6 (CMIP6) focusing on the Greenland and Antarctic Ice Sheets. In this paper, we describe the framework for ISMIP6 and its relationship to other activities within CMIP6. The ISMIP6 experimental design relies on CMIP6 climate models and includes, for the first time within CMIP, coupled ice sheet - climate models as well as standalone ice sheet models. To facilitate analysis of the multi-model ensemble and to generate a set of standard climate inputs for standalone ice sheet models, ISMIP6 defines a protocol for all variables related to ice sheets. ISMIP6 will provide a basis for investigating the feedbacks, impacts, and sea-level changes associated with dynamic ice sheets and for quantifying the uncertainty in ice-sheet-sourced global sea-level change.

Ice Sheet Model Intercomparison Project (ISMIP6) contribution to CMIP6

Reducing the uncertainty in the past, present and future contribution of ice sheets to sea-level change requires a coordinated effort between the climate and glaciology communities. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) is the primary activity within the Coupled Model Intercomparison Project - phase 6 (CMIP6) focusing on the Greenland and Antarctic Ice Sheets. In this paper, we describe the framework for ISMIP6 and its relationship to other activities within CMIP6. The ISMIP6 experimental design relies on CMIP6 climate models and includes, for the first time within CMIP, coupled ice sheet - climate models as well as standalone ice sheet models. To facilitate analysis of the multi-model ensemble and to generate a set of standard climate inputs for standalone ice sheet models, ISMIP6 defines a protocol for all variables related to ice sheets. ISMIP6 will provide a basis for investigating the feedbacks, impacts, and sea-level changes associated with dynamic ice sheets and for quantifying the uncertainty in ice-sheet-sourced global sea-level change.

Exposure age and ice-sheet model constraints on Pliocene East Antarctic ice sheet dynamics

The Late Pliocene epoch is a potential analogue for future climate in a warming world. Here we reconstruct Plio-Pleistocene East Antarctic Ice Sheet (EAIS) variability using cosmogenic nuclide exposure ages and model simulations to better understand ice sheet behaviour under such warm conditions. New and previously published exposure ages indicate interior-thickening during the Pliocene. An ice sheet model with mid-Pliocene boundary conditions also results in interior thickening and suggests that both the Wilkes Subglacial and Aurora Basins largely melted, offsetting increased ice volume. Considering contributions from West Antarctica and Greenland, this is consistent with the most recent IPCC AR5 estimate, which indicates that the Pliocene sea level likely did not exceed +20 m on Milankovitch timescales. The inception of colder climate since ∼3 Myr has increased the sea ice cover and inhibited active moisture transport to Antarctica, resulting in reduced ice sheet thickness, at least in coastal areas.

Insolation-driven 100,000-year glacial cycles and hysteresis of ice-sheet volume

The growth and reduction of Northern Hemisphere ice sheets over the past million years is dominated by an approximately 100,000-year periodicity and a sawtooth pattern (gradual growth and fast termination). Milankovitch theory proposes that summer insolation at high northern latitudes drives the glacial cycles, and statistical tests have demonstrated that the glacial cycles are indeed linked to eccentricity, obliquity and precession cycles. Yet insolation alone cannot explain the strong 100,000-year cycle, suggesting that internal climatic feedbacks may also be at work. Earlier conceptual models, for example, showed that glacial terminations are associated with the build-up of Northern Hemisphere 'excess ice', but the physical mechanisms underpinning the 100,000-year cycle remain unclear. Here we show, using comprehensive climate and ice-sheet models, that insolation and internal feedbacks between the climate, the ice sheets and the lithosphere-asthenosphere system explain the 100,000-year periodicity. The responses of equilibrium states of ice sheets to summer insolation show hysteresis, with the shape and position of the hysteresis loop playing a key part in determining the periodicities of glacial cycles. The hysteresis loop of the North American ice sheet is such that after inception of the ice sheet, its mass balance remains mostly positive through several precession cycles, whose amplitudes decrease towards an eccentricity minimum. The larger the ice sheet grows and extends towards lower latitudes, the smaller is the insolation required to make the mass balance negative. Therefore, once a large ice sheet is established, a moderate increase in insolation is sufficient to trigger a negative mass balance, leading to an almost complete retreat of the ice sheet within several thousand years. This fast retreat is governed mainly by rapid ablation due to the lowered surface elevation resulting from delayed isostatic rebound, which is the lithosphere-asthenosphere response. Carbon dioxide is involved, but is not determinative, in the evolution of the 100,000-year glacial cycles.

Habitable zone limits for dry planets

Most discussion of habitable planets has focused on Earth-like planets with globally abundant liquid water. For an "aqua planet" like Earth, the surface freezes if far from its sun, and the water vapor greenhouse effect runs away if too close. Here we show that "land planets" (desert worlds with limited surface water) have wider habitable zones than aqua planets. For planets at the inner edge of the habitable zone, a land planet has two advantages over an aqua planet: (i) the tropics can emit longwave radiation at rates above the traditional runaway limit because the air is unsaturated and (ii) the dry air creates a dry stratosphere that limits hydrogen escape. At the outer limits of the habitable zone, the land planet better resists global freezing because there is less water for clouds, snow, and ice. Here we describe a series of numerical experiments using a simple three-dimensional global climate model for Earth-sized planets. Other things (CO(2), rotation rate, surface pressure) unchanged, we found that liquid water remains stable at the poles of a low-obliquity land planet until net insolation exceeds 415 W/m(2) (170% that of modern Earth), compared to 330 W/m(2) (135%) for the aqua planet. At the outer limits, we found that a low-obliquity land planet freezes at 77%, while the aqua planet freezes at 90%. High-obliquity land and aqua planets freeze at 58% and 72%, respectively, with the poles offering the last refuge. We show that it is possible that, as the Sun brightens, an aqua planet like Earth can lose most of its hydrogen and become a land planet without first passing through a sterilizing runaway greenhouse. It is possible that Venus was a habitable land planet as recently as 1 billion years ago.

Deepwater formation in the North Pacific during the Last Glacial Termination

Between approximately 17,500 and 15,000 years ago, the Atlantic meridional overturning circulation weakened substantially in response to meltwater discharges from disintegrating Northern Hemispheric glacial ice sheets. The global effects of this reorganization of poleward heat flow in the North Atlantic extended to Antarctica and the North Pacific. Here we present evidence from North Pacific paleo surface proxy data, a compilation of marine radiocarbon age ventilation records, and global climate model simulations to suggest that during the early stages of the Last Glacial Termination, deep water extending to a depth of approximately 2500 to 3000 meters was formed in the North Pacific. A switch of deepwater formation between the North Atlantic and the North Pacific played a key role in regulating poleward oceanic heat transport during the Last Glacial Termination.