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 The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) is the primary effort of CMIP6 (Coupled Model Intercomparison Project–Phase 6) focusing on ice sheets, designed to provide an ensemble of process‐based projections of the ice‐sheet contribution to sea‐level rise over the twenty‐first century. However, the behavior of the Antarctic Ice Sheet beyond 2100 remains largely unknown: several instability mechanisms can develop on longer time scales, potentially destabilizing large parts of Antarctica. Projections of Antarctic Ice Sheet evolution until 2300 are presented here, using an ensemble of 16 ice‐flow models and forcing from global climate models. Under high‐emission scenarios, the Antarctic sea‐level contribution is limited to less than 30 cm sea‐level equivalent (SLE) by 2100, but increases rapidly thereafter to reach up to 4.4 m SLE by 2300. Simulations including ice‐shelf collapse lead to an additional 1.1 m SLE on average by 2300, and can reach 6.9 m SLE. Widespread retreat is observed on that timescale in most West Antarctic basins, leading to a collapse of large sectors of West Antarctica by 2300 in 30%–40% of the ensemble. While the onset date of retreat varies among ice models, the rate of upstream propagation is highly consistent once retreat begins. Calculations of sea‐level contribution including water density corrections lead to an additional ∼10% sea level and up to 50% for contributions accounting for bedrock uplift in response to ice loading. Overall, these results highlight large sea‐level contributions from Antarctica and suggest that the choice of ice sheet model remains the leading source of uncertainty in multi‐century projections.

The last deglaciation (21 to 8 ka) of the Eurasian ice sheet (EIS) is thought to have been responsible for a sea level rise of about 20 m. While many studies have examined the timing and rate of the EIS retreat during this period, many questions remain about the key processes that triggered the EIS deglaciation 21 kyr ago. Due to its large marine-based parts in the Barents–Kara (BKIS) and British Isles sectors, the BKIS is often considered to be a potential analogue of the current West Antarctic ice sheet (WAIS). Identifying the mechanisms that drove the EIS evolution might provide a better understanding of the processes at play in the West Antarctic destabilization. To investigate the relative impact of key drivers on the EIS destabilization, we used the three-dimensional ice sheet model GRISLI (GRenoble Ice Shelf and Land Ice) (version 2.0) forced by climatic fields from five Paleoclimate Modelling Intercomparison Project phases 3 and 4 (PMIP3, PMIP4) Last Glacial Maximum (LGM) simulations. In this study, we performed sensitivity experiments to test the response of the simulated Eurasian ice sheets to surface climate, oceanic temperatures (and thus basal melting under floating ice tongues), and sea level perturbations. Our results highlight that the EIS retreat simulated with the GRISLI model is primarily triggered by atmospheric warming. Increased atmospheric temperatures further amplify the sensitivity of the ice sheets to sub-shelf melting. These results contradict those of previous modelling studies mentioning the central role of basal melting on the deglaciation of the marine-based Barents–Kara ice sheet. However, we argue that the differences with previous works are mainly related to differences in the methodology followed to generate the initial LGM ice sheet. Due to the strong sensitivity of EIS to the atmospheric forcing highlighted with the GRISLI model and the limited extent of the confined ice shelves during the LGM, we conclude by questioning the analogy between EIS and the current WAIS. However, because of the expected rise in atmospheric temperatures, the risk of hydrofracturing is increasing and could ultimately put the WAIS in a configuration similar to the past Eurasian ice sheet.

The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from −5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty.

During the last deglaciation, the climate evolves from a cold state at the Last Glacial Maximum (LGM) at 21 ka (thousand years ago) with large ice sheets to the warm Holocene at ∼9 ka with reduced ice sheets. The deglacial ice sheet melt can impact the climate through multiple ways: changes of topography and albedo, bathymetry and coastlines, and freshwater fluxes (FWFs). In the PMIP4 (Paleoclimate Modelling Intercomparison Project – Phase 4) protocol for deglacial simulations, these changes can be accounted for or not depending on the modelling group choices. In addition, two ice sheet reconstructions are available (ICE-6G_C and GLAC-1D). In this study, we evaluate all these effects related to ice sheet changes on the climate using the iLOVECLIM model of intermediate complexity. We show that the two reconstructions yield the same warming to a first order but with a different amplitude (global mean temperature of 3.9 ∘C with ICE-6G_C and 3.8 ∘C with GLAC-1D) and evolution. We obtain a stalling of temperature rise during the Antarctic Cold Reversal (ACR, from ∼14 to ∼12 ka) similar to proxy data only with the GLAC-1D ice sheet reconstruction. Accounting for changes in bathymetry in the simulations results in a cooling due to a larger sea ice extent and higher surface albedo. Finally, freshwater fluxes result in Atlantic meridional overturning circulation (AMOC) drawdown, but the timing in the simulations disagrees with proxy data of ocean circulation changes. This questions the causal link between reconstructed freshwater fluxes from ice sheet melt and recorded AMOC weakening.

This study presents the application of an interactive downscaling in Europe using iLOVECLIM (a model of intermediate complexity), increasing its atmospheric resolution from 5.56 to 0.25∘ kilometric. A transient simulation using the appropriate climate forcings for the entire Holocene (11.5–0 ka BP) was done for both the standard version of the model and with an interactive downscaling applied. Our results show that simulations from downscaling present spatial variability that agrees better with proxy-based reconstructions and other climate models as compared to the standard model. The downscaling scheme simulates much higher (by at least a factor of 2) precipitation maxima and provides detailed information in mountainous regions. We focus on examples from the Scandes mountains, the Alps, the Scottish Highlands, and the Mediterranean. The higher spatial resolution of the downscaling provides a more realistic overview of the topography and gives local climate information, such as precipitation and temperature gradient, that is important for paleoclimate studies. With downscaling, we simulate similar trends and spatial patterns of the precipitation changes reconstructed by other proxy studies (for example in the Alps) as compared to the standard version. Our downscaling tool is numerically cheap, implying that our model can perform kilometric, multi-millennial simulations and is suitable for future studies.

Crotti et al. reconstructed the dynamics of the Wilkes Subglacial Basin (Antarctica) during the past 350,000 years. Their study reveals that a portion of the East Antarctic ice sheet experienced an extensive retreat 330,000 years ago.

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.

The Last Glacial Maximum (LGM, ∼ 21 000 years ago) has been a major focus for evaluating how well state-of-the-art climate models simulate climate changes as large as those expected in the future using paleoclimate reconstructions. A new generation of climate models has been used to generate LGM simulations as part of the Paleoclimate Modelling Intercomparison Project (PMIP) contribution to the Coupled Model Intercomparison Project (CMIP). Here, we provide a preliminary analysis and evaluation of the results of these LGM experiments (PMIP4, most of which are PMIP4-CMIP6) and compare them with the previous generation of simulations (PMIP3, most of which are PMIP3-CMIP5). We show that the global averages of the PMIP4 simulations span a larger range in terms of mean annual surface air temperature and mean annual precipitation compared to the PMIP3-CMIP5 simulations, with some PMIP4 simulations reaching a globally colder and drier state. However, the multi-model global cooling average is similar for the PMIP4 and PMIP3 ensembles, while the multi-model PMIP4 mean annual precipitation average is drier than the PMIP3 one. There are important differences in both atmospheric and oceanic circulations between the two sets of experiments, with the northern and southern jet streams being more poleward and the changes in the Atlantic Meridional Overturning Circulation being less pronounced in the PMIP4-CMIP6 simulations than in the PMIP3-CMIP5 simulations. Changes in simulated precipitation patterns are influenced by both temperature and circulation changes. Differences in simulated climate between individual models remain large. Therefore, although there are differences in the average behaviour across the two ensembles, the new simulation results are not fundamentally different from the PMIP3-CMIP5 results. Evaluation of large-scale climate features, such as land–sea contrast and polar amplification, confirms that the models capture these well and within the uncertainty of the paleoclimate reconstructions. Nevertheless, regional climate changes are less well simulated: the models underestimate extratropical cooling, particularly in winter, and precipitation changes. These results point to the utility of using paleoclimate simulations to understand the mechanisms of climate change and evaluate model performance.

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.

The Antarctic ice sheet's contribution to global sea level rise over the 21st century is of primary societal importance and remains largely uncertain as of yet. In particular, in the recent literature, the contribution of the Antarctic ice sheet by 2100 can be negative (sea level fall) by a few centimetres or positive (sea level rise), with some estimates above 1 m. The Ice Sheet Model Intercomparison Project for the Coupled Model Intercomparison Project – phase 6 (ISMIP6) aimed at reducing the uncertainties in the fate of the ice sheets in the future by gathering various ice sheet models in a common framework. Here, 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 show that our model is strongly sensitive to the climate forcing used, with a contribution of the Antarctic ice sheet to global sea level rise by 2100 that ranges from −50 to +150 mm sea level equivalent (SLE). Future oceanic warming leads to a decrease in thickness of the ice shelves, resulting in grounding-line retreat, while increased surface mass balance partially mitigates or even overcompensates the dynamic ice sheet contribution to global sea level rise. Most of the ice sheet changes over the next century are dampened under low-greenhouse-gas-emission scenarios. Uncertainties related to sub-ice-shelf melt rates induce large differences in simulated grounding-line retreat, confirming the importance of this process and its representation in ice sheet models for projections of the Antarctic ice sheet's evolution.

Antarctica's ice shelves modulate the grounded ice flow, and weakening of ice shelves due to climate forcing will decrease their ‘buttressing’ effect, causing a response in the grounded ice. While the processes governing ice-shelf weakening are complex, uncertainties in the response of the grounded ice sheet are also difficult to assess. The Antarctic BUttressing Model Intercomparison Project (ABUMIP) compares ice-sheet model responses to decrease in buttressing by investigating the ‘end-member’ scenario of total and sustained loss of ice shelves. Although unrealistic, this scenario enables gauging the sensitivity of an ensemble of 15 ice-sheet models to a total loss of buttressing, hence exhibiting the full potential of marine ice-sheet instability. All models predict that this scenario leads to multi-metre (1–12 m) sea-level rise over 500 years from present day. West Antarctic ice sheet collapse alone leads to a 1.91–5.08 m sea-level rise due to the marine ice-sheet instability. Mass loss rates are a strong function of the sliding/friction law, with plastic laws cause a further destabilization of the Aurora and Wilkes Subglacial Basins, East Antarctica. Improvements to marine ice-sheet models have greatly reduced variability between modelled ice-sheet responses to extreme ice-shelf loss, e.g. compared to the SeaRISE assessments.

Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between −7.8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between −6.1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.

The Greenland ice sheet is one of the largest contributors to global mean sea-level rise today and is expected to continue to lose mass as the Arctic continues to warm. The two predominant mass loss mechanisms are increased surface meltwater run-off and mass loss associated with the retreat of marine-terminating outlet glaciers. In this paper we use a large ensemble of Greenland ice sheet models forced by output from a representative subset of the Coupled Model Intercomparison Project (CMIP5) global climate models to project ice sheet changes and sea-level rise contributions over the 21st century. The simulations are part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). We estimate the sea-level contribution together with uncertainties due to future climate forcing, ice sheet model formulations and ocean forcing for the two greenhouse gas concentration scenarios RCP8.5 and RCP2.6. The results indicate that the Greenland ice sheet will continue to lose mass in both scenarios until 2100, with contributions of 90±50 and 32±17 mm to sea-level rise for RCP8.5 and RCP2.6, respectively. The largest mass loss is expected from the south-west of Greenland, which is governed by surface mass balance changes, continuing what is already observed today. Because the contributions are calculated against an unforced control experiment, these numbers do not include any committed mass loss, i.e. mass loss that would occur over the coming century if the climate forcing remained constant. Under RCP8.5 forcing, ice sheet model uncertainty explains an ensemble spread of 40 mm, while climate model uncertainty and ocean forcing uncertainty account for a spread of 36 and 19 mm, respectively. Apart from those formally derived uncertainty ranges, the largest gap in our knowledge is about the physical understanding and implementation of the calving process, i.e. the interaction of the ice sheet with the ocean.

Understanding the ocean circulation changes associated with abrupt climate events is key to better assessing climate variability and understanding its different natural modes. Sedimentary Pa∕Th, benthic δ13C and Δ14C are common proxies used to reconstruct past circulation flow rate and ventilation. To overcome the limitations of each proxy taken separately, a better approach is to produce multiproxy measurements on a single sediment core. Yet, different proxies can provide conflicting information about past ocean circulation. Thus, modelling them in a consistent physical framework has become necessary to assess the geographical pattern and the timing and sequence of the multiproxy response to abrupt circulation changes. We have implemented a representation of the 231Pa and 230Th tracers into the model of intermediate complexity iLOVECLIM, which already included δ13C and Δ14C. We have further evaluated the response of these three ocean circulation proxies to a classical abrupt circulation reduction obtained by freshwater addition in the Nordic Seas under preindustrial boundary conditions. The proxy response is shown to cluster in modes that resemble the modern Atlantic water masses. The clearest and most coherent response is obtained in the deep (> 2000 m) northwest Atlantic, where δ13C and Δ14C significantly decrease, while Pa∕Th increases. This is consistent with observational data across millennial-scale events of the last glacial. Interestingly, while in marine records, except in rare instances, the phase relationship between these proxies remains unclear due to large dating uncertainties, in the model the bottom water carbon isotope (δ13C and Δ14C) response lags behind the sedimentary Pa∕Th response by a few hundred years.

The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.

Climate model simulations are inherently biased. It is a notably difficult problem when dealing with climate impact assessments and model-data integration. This is especially true when looking at derived quantities such as biomes, where not only climate but also vegetation dynamics biases come into play. To overcome such difficulties, we evaluate the performance of an existing methodology to correct climate model outputs, applied here for the first time to long past climate conditions. The proposed methodology relies on the ‘Cumulative Distribution Function-transform’ (CDF-t) technique, which allows to account for climate change within the bias-correction procedure. The results are evaluated in two independent ways: (i) using forward modelling, so that model results are directly comparable to reconstructed vegetation distribution; (ii) using climatic reconstructions based on an inverse modelling approach. The modelling is performed using the intermediate complexity model iLOVECLIM in the standard global and interactively downscaled over the Europe version. The combined effects of dynamical downscaling and bias correction resulted in significantly stronger agreement between the simulated results and pollen-based biome reconstructions (BIOME6000) for the pre-industrial (0.18 versus 0.44) and mid-Holocene (MH) (0.31 versus 0.40). Higher correlation is also observed between statistically modelled global gridded potential natural distribution and modelled biomes (0.36 versus 0.41). Similarly, we find higher correlation between the reconstructed and the modelled temperatures for the MH (0.02 versus 0.21). No significant difference is found for the Last Glacial Maximum when using temperature reconstructions, due to the low number of data points available. Our findings show that the application of the CDF-t method on simulated climate variables enables us to simulate palaeoclimate and vegetation distribution in better agreement with independent reconstructions.

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.

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

In the context of global warming, growing attention is paid to the evolution of the Greenland ice sheet (GrIS) and its contribution to sea-level rise at the centennial timescale. Atmosphere–GrIS interactions, such as the temperature–elevation and the albedo feedbacks, have the potential to modify the surface energy balance and thus to impact the GrIS surface mass balance (SMB). In turn, changes in the geometrical features of the ice sheet may alter both the climate and the ice dynamics governing the ice sheet evolution. However, changes in ice sheet geometry are generally not explicitly accounted for when simulating atmospheric changes over the Greenland ice sheet in the future. To account for ice sheet–climate interactions, we developed the first two-way synchronously coupled model between a regional atmospheric model (MAR) and a 3-D ice sheet model (GRISLI). Using this novel model, we simulate the ice sheet evolution from 2000 to 2150 under a prolonged representative concentration pathway scenario, RCP8.5. Changes in surface elevation and ice sheet extent simulated by GRISLI have a direct impact on the climate simulated by MAR. They are fed to MAR from 2020 onwards, i.e. when changes in SMB produce significant topography changes in GRISLI. We further assess the importance of the atmosphere–ice sheet feedbacks through the comparison of the two-way coupled experiment with two other simulations based on simpler coupling strategies: (i) a one-way coupling with no consideration of any change in ice sheet geometry; (ii) an alternative one-way coupling in which the elevation change feedbacks are parameterized in the ice sheet model (from 2020 onwards) without taking into account the changes in ice sheet topography in the atmospheric model. The two-way coupled experiment simulates an important increase in surface melt below 2000 m of elevation, resulting in an important SMB reduction in 2150 and a shift of the equilibrium line towards elevations as high as 2500 m, despite a slight increase in SMB over the central plateau due to enhanced snowfall. In relation with these SMB changes, modifications of ice sheet geometry favour ice flux convergence towards the margins, with an increase in ice velocities in the GrIS interior due to increased surface slopes and a decrease in ice velocities at the margins due to decreasing ice thickness. This convergence counteracts the SMB signal in these areas. In the two-way coupling, the SMB is also influenced by changes in fine-scale atmospheric dynamical processes, such as the increase in katabatic winds from central to marginal regions induced by increased surface slopes. Altogether, the GrIS contribution to sea-level rise, inferred from variations in ice volume above floatation, is equal to 20.4 cm in 2150. The comparison between the coupled and the two uncoupled experiments suggests that the effect of the different feedbacks is amplified over time with the most important feedbacks being the SMB–elevation feedbacks. As a result, the experiment with parameterized SMB–elevation feedback provides a sea-level contribution from GrIS in 2150 only 2.5 % lower than the two-way coupled experiment, while the experiment with no feedback is 9.3 % lower. The change in the ablation area in the two-way coupled experiment is much larger than those provided by the two simplest methods, with an underestimation of 11.7 % (14 %) with parameterized feedbacks (no feedback). In addition, we quantify that computing the GrIS contribution to sea-level rise from SMB changes only over a fixed ice sheet mask leads to an overestimation of ice loss of at least 6 % compared to the use of a time variable ice sheet mask. Finally, our results suggest that ice-loss estimations diverge when using the different coupling strategies, with differences from the two-way method becoming significant at the end of the 21st century. In particular, even if averaged over the whole GrIS the climatic and ice sheet fields are relatively similar; at the local and regional scale there are important differences, highlighting the importance of correctly representing the interactions when interested in basin scale changes.

In this paper, we present the GRISLI (Grenoble ice sheet and land ice) model in its newest revision (version 2.0). Whilst GRISLI is applicable to any given ice sheet, we focus here on the Antarctic ice sheet because it highlights the importance of grounding line dynamics. Important improvements have been implemented in the model since its original version (Ritz et al., 2001). Notably, GRISLI now includes a basal hydrology model and an explicit flux computation at the grounding line based on the analytical formulations of Schoof (2007) or Tsai et al. (2015). We perform a full calibration of the model based on an ensemble of 300 simulations sampling mechanical parameter space using a Latin hypercube method. Performance of individual members is assessed relative to the deviation from present-day observed Antarctic ice thickness. To assess the ability of the model to simulate grounding line migration, we also present glacial–interglacial ice sheet changes throughout the last 400 kyr using the best ensemble members taking advantage of the capacity of the model to perform multi-millennial long-term integrations. To achieve this goal, we construct a simple climatic perturbation of present-day climate forcing fields based on two climate proxies: atmospheric and oceanic. The model is able to reproduce expected grounding line advances during glacial periods and subsequent retreats during terminations with reasonable glacial–interglacial ice volume changes.

This paper presents the inclusion of an online dynamical downscaling of temperature and precipitation within the model of intermediate complexity iLOVECLIM v1.1. We describe the following methodology to generate temperature and precipitation fields on a 40 km × 40 km Cartesian grid of the Northern Hemisphere from the T21 native atmospheric model grid. Our scheme is not grid specific and conserves energy and moisture in the same way as the original climate model. We show that we are able to generate a high-resolution field which presents a spatial variability in better agreement with the observations compared to the standard model. Although the large-scale model biases are not corrected, for selected model parameters, the downscaling can induce a better overall performance compared to the standard version on both the high-resolution grid and on the native grid. Foreseen applications of this new model feature include the improvement of ice sheet model coupling and high-resolution land surface models.

This study investigates the evolution of Greenland ice sheet flux focusing on five of the main fast flowing regions (Petermann glacier, North East Greenland Ice Stream, Kangerdlugssuaq glacier, Helheim glacier and Jakobshavn glacier) in response to 20th and 21st century climate change. A hybrid (shallow ice and shallow shelf) ice-sheet model (ISM) is forced with the combined outputs of a set of seven CMIP5 models and the regional climate model MAR. The ISM simulates the present-day ice velocity pattern, topography and surface mass balance (SMB) in good agreement with observations. Except for the Kangerdlugssuaq glacier, over the 21st century all the fast-flowing areas have exhibited a decrease in ice flux as a result of a negative SMB rather than dynamical changes. Only the fronts of Kangerdlugssuaq and Helheim glaciers have shown an interannual variability driven by dynamical rather than climate changes. Finally, the results predict a substantial inland ice margin retreat by the end of the 21st century, especially along the northern coasts.

The understanding of the nature and behavior of ice sheets in past warm periods is important for constraining the potential impacts of future climate change. The Pliocene warm period (between 3.264 and 3.025 Ma) saw global temperatures similar to those projected for future climates; nevertheless, Pliocene ice locations and extents are still poorly constrained. We present results from the efforts to simulate mid-Pliocene Greenland Ice Sheets by means of the international Pliocene Ice Sheet Modeling Intercomparison Project (PLISMIP). We compare the performance of existing numerical ice sheet models in simulating modern control and mid-Pliocene ice sheets with a suite of sensitivity experiments guided by available proxy records. We quantify equilibrated ice sheet volume on Greenland, identifying a potential range in sea level contributions from warm Pliocene scenarios. A series of statistical measures are performed to quantify the confidence of simulations with focus on inter-model and inter-scenario differences. We find that Pliocene Greenland Ice Sheets are less sensitive to differences in ice sheet model configurations and internal physical quantities than to changes in imposed climate forcing. We conclude that Pliocene ice was most likely to be limited to the highest elevations in eastern and southern Greenland as simulated with the highest confidence and by synthesizing available regional proxies; however, the extent of those ice caps needs to be further constrained by using a range of general circulation model (GCM) climate forcings.

Airborne and ground-based measurements of methane (CH4), carbon dioxide (CO2) and boundary layer thermodynamics were recorded over the Fennoscandian landscape (67–69.5° N, 20–28° E) in July 2012 as part of the MAMM (Methane and other greenhouse gases in the Arctic: Measurements, process studies and Modelling) field campaign. Employing these airborne measurements and a simple boundary layer box model, net regional-scale (~ 100 km) fluxes were calculated to be 1.2 ± 0.5 mg CH4 h−1 m−2 and −350 ± 143 mg CO2 h−1 m−2. These airborne fluxes were found to be relatively consistent with seasonally averaged surface chamber (1.3 ± 1.0 mg CH4 h−1 m−2) and eddy covariance (1.3 ± 0.3 mg CH4 h−1 m−2 and −309 ± 306 mg CO2 h−1 m−2) flux measurements in the local area. The internal consistency of the aircraft-derived fluxes across a wide swath of Fennoscandia coupled with an excellent statistical comparison with local seasonally averaged ground-based measurements demonstrates the potential scalability of such localised measurements to regional-scale representativeness. Comparisons were also made to longer-term regional CH4 climatologies from the JULES (Joint UK Land Environment Simulator) and HYBRID8 land surface models within the area of the MAMM campaign. The average hourly emission flux output for the summer period (July–August) for the year 2012 was 0.084 mg CH4 h−1 m−2 (minimum 0.0 and maximum 0.21 mg CH4 h−1 m−2) for the JULES model and 0.088 mg CH4 h−1 m−2 (minimum 0.0008 and maximum 1.53 mg CH4 h−1 m−2) for HYBRID8. Based on these observations both models were found to significantly underestimate the CH4 emission flux in this region, which was linked to the under-prediction of the wetland extents generated by the models.

We apply a new parameterisation of the Greenland ice sheet (GrIS) feedback between surface mass balance (SMB: the sum of surface accumulation and surface ablation) and surface elevation in the MAR regional climate model (Edwards et al., 2014) to projections of future climate change using five ice sheet models (ISMs). The MAR (Modèle Atmosphérique Régional: Fettweis, 2007) climate projections are for 2000–2199, forced by the ECHAM5 and HadCM3 global climate models (GCMs) under the SRES A1B emissions scenario. The additional sea level contribution due to the SMB–elevation feedback averaged over five ISM projections for ECHAM5 and three for HadCM3 is 4.3% (best estimate; 95% credibility interval 1.8–6.9%) at 2100, and 9.6% (best estimate; 95% credibility interval 3.6–16.0%) at 2200. In all results the elevation feedback is significantly positive, amplifying the GrIS sea level contribution relative to the MAR projections in which the ice sheet topography is fixed: the lower bounds of our 95% credibility intervals (CIs) for sea level contributions are larger than the "no feedback" case for all ISMs and GCMs. Our method is novel in sea level projections because we propagate three types of modelling uncertainty – GCM and ISM structural uncertainties, and elevation feedback parameterisation uncertainty – along the causal chain, from SRES scenario to sea level, within a coherent experimental design and statistical framework. The relative contributions to uncertainty depend on the timescale of interest. At 2100, the GCM uncertainty is largest, but by 2200 both the ISM and parameterisation uncertainties are larger. We also perform a perturbed parameter ensemble with one ISM to estimate the shape of the projected sea level probability distribution; our results indicate that the probability density is slightly skewed towards higher sea level contributions.

We present a new parameterisation that relates surface mass balance (SMB: the sum of surface accumulation and surface ablation) to changes in surface elevation of the Greenland ice sheet (GrIS) for the MAR (Modèle Atmosphérique Régional: Fettweis, 2007) regional climate model. The motivation is to dynamically adjust SMB as the GrIS evolves, allowing us to force ice sheet models with SMB simulated by MAR while incorporating the SMB–elevation feedback, without the substantial technical challenges of coupling ice sheet and climate models. This also allows us to assess the effect of elevation feedback uncertainty on the GrIS contribution to sea level, using multiple global climate and ice sheet models, without the need for additional, expensive MAR simulations. We estimate this relationship separately below and above the equilibrium line altitude (ELA, separating negative and positive SMB) and for regions north and south of 77° N, from a set of MAR simulations in which we alter the ice sheet surface elevation. These give four "SMB lapse rates", gradients that relate SMB changes to elevation changes. We assess uncertainties within a Bayesian framework, estimating probability distributions for each gradient from which we present best estimates and credibility intervals (CI) that bound 95% of the probability. Below the ELA our gradient estimates are mostly positive, because SMB usually increases with elevation: 0.56 (95% CI: −0.22 to 1.33) kg m−3 a−1 for the north, and 1.91 (1.03 to 2.61) kg m−3 a−1 for the south. Above the ELA, the gradients are much smaller in magnitude: 0.09 (−0.03 to 0.23) kg m−3 a−1 in the north, and 0.07 (−0.07 to 0.59) kg m−3 a−1 in the south, because SMB can either increase or decrease in response to increased elevation. Our statistically founded approach allows us to make probabilistic assessments for the effect of elevation feedback uncertainty on sea level projections (Edwards et al., 2014).

Microwave emissions at the L-band (1⁻2 GHz) in Antarctica are characterized by a significant contribution of ice layers at great depth, from hundreds to a thousand meters. Brightness temperatures, thus, could provide the internal temperature of the ice sheet. However, there are two difficulties to overcome in developing an accurate retrieval algorithm. First, it is difficult to know precisely from which depths waves are emanating because the ice-absorption coefficient is uncertain at the L-band, despite several formulations proposed in the literature over the past few decades. Second, emissivity potentially varies in Antarctica due to remnant scattering in firn (or ice), even at the Brewster angle, and despite the low frequency, limiting the accuracy of the estimate of the physical temperature. Here, we present a retrieval method able to disentangle the absorption and emissivity effects from brightness temperature over the whole continent. We exploit the fact that scattering and absorption are controlled by different physical parameters and phenomena that can be considered as statistically independent. This independence provides a constraint to the retrieval method, that is then well-conditioned and solvable. Our results show that (1) the retrieved absorption agrees with the permittivity model proposed by Mätzler et al. (2006), and (2) emissivity shows significant variations, up to 6% over the continent, which are correlated with wind speed and accumulation patterns. A possible cause of this latter point is density heterogeneity and sastrugi buried in the firn. These two results are an important step forward for the accurate retrieval of internal temperature using low-frequency microwave radiometers.

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.5 m of sea level equivalent, contrasting with previous studies. Our results suggest an important contribution of the Antarctic ice sheet to the LIG highstand.

Predicting the climate for the future and how it will impact ice sheet evolution requires coupling ice sheet models with climate models. However, before we attempt to develop a realistic coupled setup, we propose, in this study, to first analyse the impact of a model simulated climate on an ice sheet. We undertake this exercise for a set of regional and global climate models. Modelled near surface air temperature and precipitation are provided as upper boundary conditions to the GRISLI (GRenoble Ice Shelf and Land Ice model) hybrid ice sheet model (ISM) in its Greenland configuration. After 20 kyrs of simulation, the resulting ice sheets highlight the differences between the climate models. While modelled ice sheet sizes are generally comparable to the observed one, there are considerable deviations among the ice sheets on regional scales. These deviations can be explained by biases in temperature and precipitation near the coast. This is especially true in the case of global models. But the deviations between the climate models are also due to the differences in the atmospheric general circulation. To account for these differences in the context of coupling ice sheet models with climate models, we conclude that appropriate downscaling methods will be needed. In some cases, systematic corrections of the climatic variables at the interface may be required to obtain realistic results for the Greenland ice sheet (GIS).

A 2.5-dimensional climate model of intermediate complexity, CLIMBER-2, fully coupled with the GREMLINS 3-D thermo-mechanical ice sheet model is used to simulate the evolution of major Northern Hemisphere ice sheets during the last glacial-interglacial cycle and to investigate the ice sheets responses to both insolation and atmospheric CO2 concentration. This model reproduces the main phases of advance and retreat of Northern Hemisphere ice sheets during the last glacial cycle, although the amplitude of these variations is less pronounced than those based on sea level reconstructions. At the last glacial maximum, the simulated ice volume is 52.5×1015 m3 and the spatial distribution of both the American and Eurasian ice complexes is in reasonable agreement with observations, with the exception of the marine parts of these former ice sheets. A set of sensitivity studies has also been performed to assess the sensitivity of the Northern Hemisphere ice sheets to both insolation and atmospheric CO2. Our results suggest that the decrease of summer insolation is the main factor responsible for the early build up of the North American ice sheet around 120 kyr BP, in agreement with benthic foraminifera δ18O signals. In contrast, low insolation and low atmospheric CO2 concentration are both necessary to trigger a long-lasting glaciation over Eurasia.