Search

We present transient simulations of the last glacial inception using the Earth system model CLIMBER-X with dynamic vegetation, interactive ice sheets, and visco-elastic solid Earth responses. The simulations are initialized at the middle of the Eemian interglacial (125 kiloyears before present, ka) and run until 100 ka, driven by prescribed changes in Earth's orbital parameters and greenhouse gas concentrations from ice core data. CLIMBER-X simulates a rapid increase in Northern Hemisphere ice sheet area through MIS5d, with ice sheets expanding over northern North America and Scandinavia, in broad agreement with proxy reconstructions. While most of the increase in ice sheet area occurs over a relatively short period between 119 and 117 ka, the larger part of the increase in ice volume occurs afterwards with an almost constant ice sheet extent. We show that the vegetation feedback plays a fundamental role in controlling the ice sheet expansion during the last glacial inception. In particular, with prescribed present-day vegetation the model simulates a global sea level drop of only ∼ 20 m, compared with the ∼ 35 m decrease in sea level with dynamic vegetation response. The ice sheet and carbon cycle feedbacks play only a minor role during the ice sheet expansion phase prior to ∼ 115 ka but are important in limiting the deglaciation during the following phase characterized by increasing summer insolation. The model results are sensitive to climate model biases and to the parameterization of snow albedo, while they show only a weak dependence on changes in the ice sheet model resolution and the acceleration factor used to speed up the climate component. Overall, our simulations confirm and refine previous results showing that climate–vegetation–cryosphere feedbacks play a fundamental role in the transition from interglacial to glacial states characterizing Quaternary glacial cycles.

The origin of the 100 kyr cyclicity, which dominates ice volume variations and other climate records over the past million years, remains debatable. Here, using a comprehensive Earth system model of intermediate complexity, we demonstrate that both strong 100 kyr periodicity in the ice volume variations and the timing of glacial terminations during past 800 kyr can be successfully simulated as direct, strongly nonlinear responses of the climate-cryosphere system to orbital forcing alone, if the atmospheric CO2 concentration stays below its typical interglacial value. The existence of long glacial cycles is primarily attributed to the North American ice sheet and requires the presence of a large continental area with exposed rocks. We show that the sharp, 100 kyr peak in the power spectrum of ice volume results from the long glacial cycles being synchronized with the Earth's orbital eccentricity. Although 100 kyr cyclicity can be simulated with a constant CO2 concentration, temporal variability in the CO2 concentration plays an important role in the amplification of the 100 kyr cycles.

Using a new approach to force an ice sheet model, we performed an ensemble of simulations of the Greenland Ice Sheet evolution during the last two glacial cycles, with emphasis on the Eemian Interglacial. This ensemble was generated by perturbing four key parameters in the coupled regional climate-ice sheet model and by introducing additional uncertainty in the prescribed "background" climate change. The sensitivity of the surface melt model to climate change was determined to be the dominant driver of ice sheet instability, as reflected by simulated ice sheet loss during the Eemian Interglacial period. To eliminate unrealistic parameter combinations, constraints from present-day and paleo information were applied. The constraints include (i) the diagnosed present-day surface mass balance partition between surface melting and ice discharge at the margin, (ii) the modeled present-day elevation at GRIP; and (iii) the modeled elevation reduction at GRIP during the Eemian. Using these three constraints, a total of 360 simulations with 90 different model realizations were filtered down to 46 simulations and 20 model realizations considered valid. The paleo constraint eliminated more sensitive melt parameter values, in agreement with the surface mass balance partition assumption. The constrained simulations resulted in a range of Eemian ice loss of 0.4–4.4 m sea level equivalent, with a more likely range of about 3.7–4.4 m sea level if the GRIP δ18O isotope record can be considered an accurate proxy for the precipitation-weighted annual mean temperatures.

A new version of the Earth system model of intermediate complexity, CLIMBER-2, which includes the three-dimensional polythermal ice-sheet model SICOPOLIS, is used to simulate the last glacial cycle forced by variations of the Earth's orbital parameters and atmospheric concentration of major greenhouse gases. The climate and ice-sheet components of the model are coupled bi-directionally through a physically-based surface energy and mass balance interface. The model accounts for the time-dependent effect of aeolian dust on planetary and snow albedo. The model successfully simulates the temporal and spatial dynamics of the major Northern Hemisphere (NH) ice sheets, including rapid glacial inception and strong asymmetry between the ice-sheet growth phase and glacial termination. Spatial extent and elevation of the ice sheets during the last glacial maximum agree reasonably well with palaeoclimate reconstructions. A suite of sensitivity experiments demonstrates that simulated ice-sheet evolution during the last glacial cycle is very sensitive to some parameters of the surface energy and mass-balance interface and dust module. The possibility of a considerable acceleration of the climate ice-sheet model is discussed.

We investigate glacial inception and glacial thresholds in the climate-cryosphere system utilising the Earth system model of intermediate complexity CLIMBER-2, which includes modules for atmosphere, terrestrial vegetation, ocean and interactive ice sheets. The latter are described by the three-dimensional polythermal ice-sheet model SICOPOLIS. A bifurcation which represents glacial inception is analysed with two different model setups: one setup with dynamical ice-sheet model and another setup without it. The respective glacial thresholds differ in terms of maximum boreal summer insolation at 65° N (hereafter referred as Milankovitch forcing (MF)). The glacial threshold of the configuration without ice-sheet dynamics corresponds to a much lower value of MF compared to the full model. If MF attains values only slightly below the aforementioned threshold there is fast transient response. Depending on the value of MF relative to the glacial threshold, the transient response time of inland-ice volume in the model configuration with ice-sheet dynamics ranges from 10 000 to 100 000 years. Due to these long response times, a glacial threshold obtained in an equilibrium simulation is not directly applicable to the transient response of the climate-cryosphere system to time-dependent orbital forcing. It is demonstrated that in transient simulations just crossing of the glacial threshold does not imply large-scale glaciation of the Northern Hemisphere. We found that in transient simulations MF has to drop well below the glacial threshold determined in an equilibrium simulation to initiate glacial inception. Finally, we show that the asynchronous coupling between climate and inland-ice components allows one sufficient realistic simulation of glacial inception and, at the same time, a considerable reduction of computational costs.