Description: Establishing whether past millennial-scale climate change affected the stability of marine methane hydrate is important for our understanding of climatic change and determining the fate of marine hydrates in a future warmer world. We show, using three-dimensional seismic data offshore of Mauritania, that episodic, millennial-scale shifts of the base of the hydrate stability zone can be imaged below the ocean floor. Process modeling suggests that the base of the hydrate stability zone should have shallowed and deepened in response to climate change over the past ∼150,000 yr. Specifically, there is seismic evidence for millennial-scale shifts during the Holocene (∼11,700 yr ago to present) at a temporal resolution that has previously been unrealized. This is the first evidence that millennial-scale climatic cycles caused hydrate formation and dissociation and that hydrate instability should be expected in a warming world.

Description: In spite of significant progress in paleoclimate reconstructions and modelling of different aspects of the past glacial cycles, the mechanisms which transform regional and seasonal variations in solar insolation into long-term and global-scale glacial–interglacial cycles are still not fully understood – in particular, in relation to CO2 variability. Here using the Earth system model of intermediate complexity CLIMBER-2 we performed simulations of the co-evolution of climate, ice sheets, and carbon cycle over the last 400000 years using the orbital forcing as the only external forcing. The model simulates temporal dynamics of CO2, global ice volume, and other climate system characteristics in good agreement with paleoclimate reconstructions. These results provide strong support for the idea that long and strongly asymmetric glacial cycles of the late Quaternary represent a direct but strongly nonlinear response of the Northern Hemisphere ice sheets to orbital forcing. This response is strongly amplified and globalised by the carbon cycle feedbacks. Using simulations performed with the model in different configurations, we also analyse the role of individual processes and sensitivity to the choice of model parameters. While many features of simulated glacial cycles are rather robust, some details of CO2 evolution, especially during glacial terminations, are sensitive to the choice of model parameters. Specifically, we found two major regimes of CO2 changes during terminations: in the first one, when the recovery of the Atlantic meridional overturning circulation (AMOC) occurs only at the end of the termination, a pronounced overshoot in CO2 concentration occurs at the beginning of the interglacial and CO2 remains almost constant during the interglacial or even declines towards the end, resembling Eemian CO2 dynamics. However, if the recovery of the AMOC occurs in the middle of the glacial termination, CO2 concentration continues to rise during the interglacial, similar to the Holocene. We also discuss the potential contribution of the brine rejection mechanism for the CO2 and carbon isotopes in the atmosphere and the ocean during the past glacial termination.

Description: Glacial cycles of the late Quaternary are controlled by the asymmetrically varying mass balance of continental ice sheets in the Northern Hemisphere. Surface mass balance is governed by processes of ablation and accumulation. Here two ablation schemes, the positive-degree-day (PDD) method and the surface energy balance (SEB) approach, are compared in transient simulations of the last glacial cycle with the Earth system model of intermediate complexity CLIMBER-2. The standard version of the CLIMBER-2 model incorporates the SEB approach and simulates ice volume variations in reasonable agreement with paleoclimate reconstructions during the entire last glacial cycle. Using results from the standard CLIMBER-2 model version, we simulated ablation with the PDD method in offline mode by applying different combinations of three empirical parameters of the PDD scheme. We found that none of the parameter combinations allow us to simulate a surface mass balance of the American and European ice sheets that is similar to that obtained with the standard SEB method. The use of constant values for the empirical PDD parameters led either to too much ablation during the first phase of the last glacial cycle or too little ablation during the final phase. We then substituted the standard SEB scheme in CLIMBER-2 with the PDD scheme and performed a suite of fully interactive (online) simulations of the last glacial cycle with different combinations of PDD parameters. The results of these simulations confirmed the results of the offline simulations: no combination of PDD parameters realistically simulates the evolution of the ice sheets during the entire glacial cycle. The use of constant parameter values in the online simulations leads either to a buildup of too much ice volume at the end of glacial cycle or too little ice volume at the beginning. Even when the model correctly simulates global ice volume at the last glacial maximum (21 ka), it is unable to simulate complete deglaciation during the Holocene. According to our simulations, the SEB approach proves superior for simulations of glacial cycles.

Description: Palaeo data suggest that Greenland must have been largely ice free during Marine Isotope Stage 11 (MIS-11). However, regional summer insolation anomalies were modest during this time compared to MIS-5e, when the Greenland ice sheet likely lost less volume. Thus it remains unclear how such conditions led to an almost complete disappearance of the ice sheet. Here we use transient climate–ice sheet simulations to simultaneously constrain estimates of regional temperature anomalies and Greenland’s contribution to the MIS-11 sea-level highstand. We find that Greenland contributed 6.1 m (3.9–7.0 m, 95% credible interval) to sea level, ∼7 kyr after the peak in regional summer temperature anomalies of 2.8 °C (2.1–3.4 °C). The moderate warming produced a mean rate of mass loss in sea-level equivalent of only around 0.4 m per kyr, which means the long duration of MIS-11 interglacial conditions around Greenland was a necessary condition for the ice sheet to disappear almost completely.

Description: Multiple proxy data reveal that the early to middle Holocene (ca. 8–6 kyr B.P.) was warmer than the preindustrial period in most regions of the Northern Hemisphere. This warming is presumably explained by the higher summer insolation in the Northern Hemisphere, owing to changes in the orbital parameters. Subsequent cooling in the late Holocene was accompanied by significant changes in vegetation cover and an increase in atmospheric CO2 concentration. The essential question is whether it is possible to explain these changes in a consistent way, accounting for the orbital parameters as the main external forcing for the climate system. We investigate this problem using the computationally efficient model of climate system, CLIMBER‐2, which includes models for oceanic and terrestrial biogeochemistry. We found that changes in climate and vegetation cover in the northern subtropical and circumpolar regions can be attributed to the changes in the orbital forcing. Explanation of the atmospheric CO2 record requires an additional assumption of excessive CaCO3 sedimentation in the ocean. The modeled decrease in the carbonate ion concentration in the deep ocean is similar to that inferred from CaCO3 sediment data [Broecker et al., 1999]. For 8 kyr B.P., the model estimates the terrestrial carbon pool ca. 90 Pg higher than its preindustrial value. Simulated atmospheric δ13C declines during the course of the Holocene, similar to δ13C data from the Taylor Dome ice core [Indermühle et al., 1999]. Amplitude of simulated changes in δ13C is smaller than in the data, while a difference between the model and the data is comparable with the range of data uncertainty.

Description: The surface energy and mass balance of ice sheets strongly depends on the amount of solar radiation absorbed at the surface, which is mainly controlled by the albedo of snow and ice. Here, using an Earth system model of intermediate complexity, we explore the role played by surface albedo for the simulation of glacial cycles. We show that the evolution of the Northern Hemisphere ice sheets over the last glacial cycle is very sensitive to the representation of snow albedo in the model. It is well known that the albedo of snow depends strongly on snow grain size and the content of light-absorbing impurities. Excluding either the snow aging effect or the dust darkening effect on snow albedo leads to an excessive ice build-up during glacial times and consequently to a failure in simulating deglaciation. While the effect of snow grain growth on snow albedo is well constrained, the albedo reduction due to the presence of dust in snow is much more uncertain because the light-absorbing properties of dust vary widely as a function of dust mineral composition. We also show that assuming slightly different optical properties of dust leads to very different ice sheet and climate evolutions in the model. Conversely, ice sheet evolution is less sensitive to the choice of ice albedo in the model. We conclude that a proper representation of snow albedo is a fundamental prerequisite for a successful simulation of glacial cycles.

Description: We introduce the coupled model of the Greenland glacial system IGLOO 1.0, including the polythermal ice sheet model SICOPOLIS (version 3.3) with hybrid dynamics, the model of basal hydrology HYDRO and a parameterization of submarine melt for marine-terminated outlet glaciers. Aim of this glacial system model is to gain a better understanding of the processes important for the future contribution of the Greenland ice sheet to sea level rise under future climate change scenarios. The ice sheet is initialized via a relaxation towards observed surface elevation, imposing the palaeo-surface temperature over the last glacial cycle. As a present-day reference, we use the 1961-1990 standard climatology derived from simulations of the regional atmosphere model MAR with ERA reanalysis boundary conditions. For the palaeo-part of the spin-up, we add the temperature anomaly derived from the GRIP ice core to the years 1961–1990 average surface temperature field. For our projections, we apply surface temperature and surface mass balance anomalies derived from RCP 4.5 and RCP 8.5 scenarios created by MAR with boundary conditions from simulations with three CMIP5 models. The hybrid ice sheet model is fully coupled with the model of basal hydrology. With this model and the MAR scenarios, we perform simulations to estimate the contribution of the Greenland ice sheet to future sea level rise until the end of the 21st and 23rd centuries. Further on, the impact of elevation-surface mass balance feedback, introduced via the MAR data, on future sea level rise is inspected. In our projections, we found the Greenland ice sheet to contribute to global sea level rise between 1.9 and 13.0cm until the year 2100 and between 3.5 and 76.4cm until the year 2300, including our simulated additional sea level rise due to elevation-surface mass balance feedback. Translated into additional sea level rise, the strength of this feedback in the year 2100 varies from 0.4 to 1.7cm, and in the year 2300 it ranges from 1.7 to 21.8cm. Additionally, taking Helheim and Store Glaciers as examples, we investigate the role of ocean warming and surface runoff change for the melting of outlet glaciers. It shows that ocean temperature and subglacial discharge are about equally important for the melting of the examined outlet glaciers.

Description: We discuss potential transitions of six climatic subsystems with large-scale impact on Europe, sometimes denoted as tipping elements. These are the ice sheets on Greenland and West Antarctica, the Atlantic thermohaline circulation, Arctic sea ice, Alpine glaciers and northern hemisphere stratospheric ozone. Each system is represented by co-authors actively publishing in the corresponding field. For each subsystem we summarize themechanism of a potential transition in a warmer climate along with its impact on Europe and assess the likelihood for such a transition based on published scientific literature. As a summary, the ‘tipping’ potential for each system is provided as a function of global mean temperature increase which required some subjective interpretation of scientific facts by the authors and should be considered as a snapshot of our current understanding.

Description: The geological record shows that abrupt changes in the Earth system can occur on timescales short enough to challenge the capacity of human societies to adapt to environmental pressures. In many cases, abrupt changes arise from slow changes in one component of the Earth system that eventually pass a critical threshold, or tipping point, after which impacts cascade through coupled climate–ecological–social systems. The chance of detecting abrupt changes and tipping points increases with the length of observations. The geological record provides the only long-term information we have on the conditions and processes that can drive physical, ecological and social systems into new states or organizational structures that may be irreversible within human time frames. Here, we use well-documented abrupt changes of the past 30 kyr to illustrate how their impacts cascade through the Earth system. We review useful indicators of upcoming abrupt changes, or early warning signals, and provide a perspective on the contributions of palaeoclimate science to the understanding of abrupt changes in the Earth system.

Description: Understanding the future fate of the Greenland Ice Sheet (GIS) in the context of anthropogenic CO2 emissions is crucial to predict sea level rise. With the fully coupled Earth system model of intermediate complexity CLIMBER-X, we study the stability of the GIS and its transient response to CO2 emissions over the next 10 Kyr. Bifurcation points exist at global temperature anomalies of 0.6 and 1.6 K relative to pre-industrial. For system states in the vicinity of the equilibrium ice volumes corresponding to these temperature anomalies, mass loss rate and sensitivity of mass loss to cumulative CO2 emission peak. These critical ice volumes are crossed for cumulative emissions of 1,000 and 2,500 GtC, which would cause long-term sea level rise by 1.8 and 6.9 m respectively. In summary, we find tipping of the GIS within the range of the temperature limits of the Paris agreement. Key Points Bifurcation points exist at global mean temperature anomalies of 0.6 and 1.6 K relative to pre-industrial Mass loss rate and sensitivity to cumulative CO2 emission peak near the equilibrium ice volumes belonging to these temperature anomalies Substantial long-term mass loss of the Greenland ice sheet for cumulative emissions larger than 1,000 Gt carbon