Search

Glacial terminations are marked by a re-organisation of the different components of the climate system. In particular, rapid ice sheet disintegration leads to multiple complex feedback loops that are still poorly understood. To further investigate this aspect, we use here a fully coupled Northern Hemisphere ice sheet–climate model to perform numerical experiments of the last two glacial terminations. We show that even if the first-order climate trajectory is similar for the two terminations, the difference in terms of solar insolation leads to important changes for the ice sheet–climate system. Warmer temperatures during the penultimate termination are compatible with higher sea level during the last interglacial period with respect to the Holocene. We simulate a last interglacial Greenland contribution to sea level rise of about 2 m of sea level equivalent. We also simulate warmer subsurface Southern Ocean, compatible with an additional contribution from the Antarctic ice sheet. In addition, even without considering freshwater flux to the ocean resulting from ice sheet melting, the two terminations display different Atlantic overturning circulation sensitivity, this circulation being more prone to collapses during the penultimate termination. Finally, with additional sensitivity experiments we show that, for the two terminations, the Northern Hemisphere insolation is the main driver for the ice sheet retreat even if vegetation changes have also to be taken into account to simulate the full deglaciation. Conversely, even though it impacts the temperature, greenhouse gas concentration change alone does not explain the amplitude of ice sheet retreat and only modulates its timing. [ABSTRACT FROM AUTHOR]

During the Quaternary, ice sheets experienced several retreat-advanced cycles, strongly influencing climate patterns. In order to properly simulate these phenomena, it is preferable to use physics-based models instead of parameterizations to estimate surface mass balance (SMB) which has a strong influence on the ice sheet evolution. To further investigate the potential of these SMB models, this work evaluates BESSI (BErgen Snow Simulator), a multi-layer snow model with high computational efficiency, as an alternative to providing SMB for paleo studies. First, we validate the snow model using the regional climate model MAR (Modεave;le Atmosphérique Régional) as forcing and reference for the present-day climate over Greenland and Antarctic Ice Sheets. The evolution of SMB over the Last Interglacial period (LIG) (130–116 kaBP) is computed by forcing BESSI with transient climate forcing obtained from an Earth system model i LOVECLIM for both ice sheets. For present-day climate conditions, BESSI exhibits good performance compared to MAR despite a much simpler model set-up. The model also captures well the variation of SMB and its components during the LIG. Compared to the current simple melt estimation scheme of i LOVECLIM (ITM), BESSI is able to capture different SMB patterns for two particular ice sheet climate conditions thanks to its higher physical constraints while ITM displays a strong sensitivity to its parameters and input fields (temperature). The findings suggest that BESSI can provide more reliable SMB estimations for the i LOVECLIM framework to improve the model simulations of the ice sheet evolution and interactions with climate. [ABSTRACT FROM AUTHOR]

Glacial terminations are marked by a re-organisation of the different components of the climate system. In particular, rapid ice sheet disintegration leads to multiple complex feedback loops that are still poorly understood. To further investigate this aspect, we use here a fully coupled Northern Hemisphere ice sheet – climate model to perform numerical experiments of the last two glacial terminations. We show that even if the first-order climate trajectory is similar for the two terminations, the difference in terms of solar insolation leads to important changes for the ice sheet – climate system. Warmer temperatures during the penultimate termination are compatible with higher sea level during the last interglacial period with respect to the Holocene. We simulate a last interglacial Greenland contribution to sea level rise of about 2 m of sea level equivalent. We also simulate warmer subsurface Southern Ocean, compatible with an additional contribution from the Antarctic ice sheet. In addition, even without considering freshwater flux to the ocean resulting from ice sheet melting, the two terminations display different Atlantic overturning circulation sensitivity, this circulation being more prone to collapses during the penultimate termination. Finally, with additional sensitivity experiments we show that, for the two terminations, the Northern Hemisphere insolation is the main driver for the ice sheet retreat even if vegetation changes have also to be taken into account to simulate the full deglaciation. Conversely, even though it impacts the temperature, greenhouse gas concentration change alone does not explain the amplitude of ice sheet retreat and only modulate its timing. [ABSTRACT FROM AUTHOR]

The last deglaciation offers an unique opportunity to understand the climate–ice-sheet interactions in a global warming context. In this paper, to tackle this question, we use an Earth system model of intermediate complexity coupled to an ice sheet model covering the Northern Hemisphere to simulate the last deglaciation and the Holocene (26–0 ka). We use a synchronous coupling every year between the ice sheet and the rest of the climate system and we ensure a closed water cycle considering the release of freshwater flux to the ocean due to ice sheet melting. Our reference experiment displays a gradual warming in response to the forcings, with no abrupt changes. In this case, while the amplitude of the freshwater flux to the ocean induced by ice sheet retreat is realistic, it is sufficient to shut down the Atlantic meridional overturning circulation from which the model does not recover within the time period simulated. However, with reduced freshwater flux we are nonetheless able to obtain different oceanic circulation evolutions, including some abrupt transitions between shut-down and active circulation states in the course of the deglaciation. The inclusion of a parameterisation for the sinking of brines around Antarctica also produces an abrupt recovery of the Atlantic meridional overturning circulation, absent in the reference experiment. The fast oceanic circulation recoveries lead to abrupt warming phases in Greenland. Our simulated ice sheet geometry evolution is in overall good agreement with available global reconstructions, even though the abrupt sea level rise at 14.6 ka is underestimated, possibly because the climate model underestimates the millennial-scale temperature variability. In the course of the deglaciation, large-scale grounding line instabilities are simulated both for the Eurasian and North American ice sheets. The first instability occurs in the Barents–Kara seas for the Eurasian ice sheet at 14.5 ka. A second grounding line instability occurs ca. 12 ka in the proglacial lake that formed at the southern margin of the North American ice sheet. With additional asynchronously coupled experiments, we assess the sensitivity of our results to different ice sheet model choices related to surface and sub-shelf mass balance, ice deformation and grounding line representation. While the ice sheet evolutions differ within this ensemble, the global climate trajectory is only weakly affected by these choices. In our experiments, only the abrupt shifts in the oceanic circulation due to freshwater fluxes are able to produce some millennial-scale variability since no self-generating abrupt transitions are simulated without these fluxes. [ABSTRACT FROM AUTHOR]

Polar amplification will result in amplified temperature changes in the Arctic with respect to the rest of the globe, making the Greenland ice sheet particularly vulnerable to global warming. While the ice sheet has been showing an increased mass loss in the past decades, its contribution to global sea level rise in the future is of primary importance since it is at present the largest single-source contribution after the thermosteric contribution. The question of the fate of the Greenland and Antarctic ice sheets for the next century has recently gathered various ice sheet models in a common framework within the Ice Sheet Model Intercomparison Project for the Coupled Model Intercomparison Project – phase 6 (ISMIP6). While in a companion paper we present the GRISLI-LSCE (Grenoble Ice Sheet and Land Ice model of the Laboratoire des Sciences du Climat et de l'Environnement) contribution to ISMIP6-Antarctica, we present here the GRISLI-LSCE contribution to ISMIP6-Greenland. We show an important spread in the simulated Greenland ice loss in the future depending on the climate forcing used. The contribution of the ice sheet to global sea level rise in 2100 can thus be from as low as 20 mm sea level equivalent (SLE) to as high as 160 mm SLE. Amongst the models tested in ISMIP6, the Coupled Model Intercomparison Project – phase 6 (CMIP6) models produce larger ice sheet retreat than their CMIP5 counterparts. Low-emission scenarios in the future drastically reduce the ice mass loss. The oceanic forcing contributes to about 10 mm SLE in 2100 in our simulations. In addition, the dynamical contribution to ice thickness change is small compared to the impact of surface mass balance. This suggests that mass loss is mostly driven by atmospheric warming and associated ablation at the ice sheet margin. With additional sensitivity experiments we also show that the spread in mass loss is only weakly affected by the choice of the ice sheet model mechanical parameters. [ABSTRACT FROM AUTHOR]

Providing reliable projections of the ice sheet contribution to future sea-level rise has become one of the main challenges of the ice sheet modelling community. To increase confidence in future projections, a good knowledge of the present-day state of ice flow dynamics, which is critically dependent on basal conditions, is strongly needed. The main difficulty is tied to the scarcity of observations at the ice–bed interface at the scale of the whole ice sheet, resulting in poorly constrained parameterisations in ice sheet models. To circumvent this drawback, inverse modelling approaches can be developed to infer initial conditions for ice sheet models that best reproduce available data. Most often such approaches allow for a good representation of the mean present-day state of the ice sheet but are accompanied with unphysical trends. Here, we present an initialisation method for the Greenland ice sheet using the thermo-mechanical hybrid GRISLI (GRenoble Ice Shelf and Land Ice) ice sheet model. Our approach is based on the adjustment of the basal drag coefficient that relates the sliding velocities at the ice–bed interface to basal shear stress in unfrozen bed areas. This method relies on an iterative process in which the basal drag is periodically adjusted in such a way that the simulated ice thickness matches the observed one. The quality of the method is assessed by computing the root mean square errors in ice thickness changes. Because the method is based on an adjustment of the sliding velocities only, the results are discussed in terms of varying ice flow enhancement factors that control the deformation rates. We show that this factor has a strong impact on the minimisation of ice thickness errors and has to be chosen as a function of the internal thermal state of the ice sheet (e.g. a low enhancement factor for a warm ice sheet). While the method performance slightly increases with the duration of the minimisation procedure, an ice thickness root mean square error (RMSE) of 50.3 m is obtained in only 1320 model years. This highlights a rapid convergence and demonstrates that the method can be used for computationally expensive ice sheet models. [ABSTRACT FROM AUTHOR]

As pointed out by the forth assessment report of the Intergovernmental Panel on Climate Change, IPCC-AR4 (Meehl et al., 2007), the contribution of the two major ice sheets, Antarctica and Greenland, to global sea level rise, is a subject of key importance for the scientific community. By the end of the next century, a 3-5 °C warming is expected in Greenland. Similar temperatures in this region were reached during the last interglacial (LIG) period, 130-115 ka BP, due to a change in orbital configuration rather than to an anthropogenic forcing. Ice core evidence suggests that the Greenland ice sheet (GIS) survived this warm period, but great uncertainties remain about the total Greenland ice reduction during the LIG. Here we perform long-term simulations of the GIS using an improved ice sheet model. Both the methodologies chosen to reconstruct palaeoclimate and to calibrate the model are strongly based on proxy data. We suggest a relatively low contribution to LIG sea level rise from Greenland melting, ranging from 0.7 to 1.5m of sea level equivalent, contrasting with previous studies. Our results suggest an important contribution of the Antarctic ice sheet to the LIG highstand. [ABSTRACT FROM AUTHOR]

The article presents study which examines the sensitivity of Greenland ice sheet (GSI) model to atmospheric forcing fields. It states that modelled precipitation and near surface air temperature is provided as upper boundary conditions to the GRenoble Ice Shelf and Land Ice model (GRISLI) hybrid ice sheet model (ISM) in Greenland. Results show that the model is much more sensitive to surface mass balance than to the dynamic parameters.