Description: Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and inform on the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimated the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes and the forcings employed. This study presents results from 18 simulations from 15 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100, forced with different scenarios from the Coupled Model Intercomparison Project Phase 5 (CMIP5) representative of the spread in climate model results. The contribution of the Antarctic ice sheet in response to increased warming during this period varies between −7.8 and 30.0 cm of Sea Level Equivalent (SLE). The evolution of the West Antarctic Ice Sheet varies widely among models, with an overall mass loss up to 21.0 cm SLE in response to changes in oceanic conditions. East Antarctica mass change varies between −6.5 and 16.5 cm SLE, 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 mass loss of 8 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from 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 AOGCMs show an overall mass loss of 10 mm SLE compared to simulations done under present-day conditions, with limited mass gain in East Antarctica.

Description: Ice-penetrating radar1,2,3 and ice core drilling4 have shown that large parts of the north-central Greenland ice sheet are melting from below. It has been argued that basal ice melt is due to the anomalously high geothermal flux1,4 that has also influenced the development of the longest ice stream in Greenland1. Here we estimate the geothermal flux beneath the Greenland ice sheet and identify a 1,200-km-long and 400-km-wide geothermal anomaly beneath the thick ice cover. We suggest that this anomaly explains the observed melting of the ice sheet’s base, which drives the vigorous subglacial hydrology3 and controls the position of the head of the enigmatic 750-km-long northeastern Greenland ice stream5. Our combined analysis of independent seismic, gravity and tectonic data6,7,8,9 implies that the geothermal anomaly, which crosses Greenland from west to east, was formed by Greenland’s passage over the Iceland mantle plume between roughly 80 and 35 million years ago. We conclude that the complexity of the present-day subglacial hydrology and dynamic features of the north-central Greenland ice sheet originated in tectonic events that pre-date the onset of glaciation in Greenland by many tens of millions of years.

Description: We introduce the coupled model of the Green- land glacial system IGLOO 1.0, including the polythermal ice sheet model SICOPOLIS (version 3.3) with hybrid dy- namics, the model of basal hydrology HYDRO and a param- eterization of submarine melt for marine-terminated outlet glaciers. The 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, impos- ing 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 atmo- sphere model MAR with ERA reanalysis boundary condi- tions. For the palaeo-part of the spin-up, we add the temper- ature anomaly derived from the GRIP ice core to the years 1961–1990 average surface temperature field. For our pro- jections, we apply surface temperature and surface mass bal- ance anomalies derived from RCP 4.5 and RCP 8.5 scenar- ios created by MAR with boundary conditions from simula- tions 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. Fur- ther on, the impact of elevation–surface mass balance feed- back, introduced via the MAR data, on future sea level rise is inspected. In our projections, we found the Greenland ice sheet to contribute between 1.9 and 13.0 cm to global sea level rise until the year 2100 and between 3.5 and 76.4 cm 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.7 cm, and in the year 2300 it ranges from 1.7 to 21.8 cm. Additionally, taking the Helheim and Store glaciers as examples, we inves- tigate the role of ocean warming and surface runoff change for the melting of outlet glaciers. It shows that ocean temper- ature and subglacial discharge are about equally important for the melting of the examined outlet glaciers.

Description: Here, we define a test bed for fast flow regions and its vicinity embedded in an ice sheet. This test bed is designed for outlet glaciers and ice streams of the Greenland ice sheet. It consists of a fine resolution part with a manufactured basal trough over which the professional software COMSOL (Multiphysics Modeling Software) operates as a full- Stokes model. Results by COMSOL are compared with coarse resolution simulations with the ice-sheet model SICOPOLIS operating in shallow-ice-approximation mode and using parameterizations of the fast flow effects. For simplification, in this preliminary approach, both models run in isothermal mode. Definition of surface mass balance follows the EISMINT intercomparison project with parameters adapted to the Greenland ice sheet. In particular, we inspect with this test bed upstream and lateral flow effects of ice streams and outlet glaciers. We present first simulations with this approach, although presentation of the test bed itself is the main emphasis of this presentation.

Description: Over the past two decades net mass loss from the Greenland ice sheet quadrupled, caused by enhanced surface melting and speedup of the marine-terminating outlet glaciers. This speedup has been related, among other factors, to enhanced submarine melting, which in turn is caused by warming of the surrounding ocean and by increased subglacial discharge. For the future and recent mass balance changes of the Greenland Ice Sheet, ice-ocean processes potentially play an important role, yet they are not properly represented in contemporary Greenland Ice Sheet models. In this work we performed numerical experiments with a one-dimensional plume model coupled to a one-dimensional model of outlet glacier. We investigate the response of a coupled ice-flow plume model to possible outcomes of climate change. In particularly, we examine the transient and equilibrium response of the outlet glaciers to changes in ocean temperature and subglacial discharge which affects both: glacier geometry and submarine melt rates.

Description: We present simulations with the dynamic/thermodynamic ice sheet model SICOPOLIS (version 3) coupled to HYDRO, a model of basal hydrology. SICOPOLIS describes the evolution of ice thickness, temperature and water content of ice sheets. Recently, the treatment of longitudinal and lateral stresses (“shelfy stream approximation”) for the dynamics and the enthalpy method as an alternative method for solving the energy equation were included into the model. HYDRO describes the basal water transport using the hydrological potential. In a bi-directional coupling, HYDRO receives the basal water fluxes from SICOPOLIS, while the basal water from HYDRO affects the basal sliding in SICOPOLIS. Here, we present offline simulations with HYDRO as well as simulations with SICOPOLIS-only and the coupled model SICOPOLIS-HYDRO. Several sensitivity studies highlight the importance of basal processes. In particular, we inspect the role of horizontal resolution. It shows that not only horizontal resolution plays an important role for resolving outlet glaciers, but also the coupled model better reproduces outlet glaciers compared to the uncoupled one; even the North-East-Greenland Ice Stream is modelled quite well without the need for special regional tuning.

Description: 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.

Description: The land ice contribution to global mean sea level rise has not yet been predicted1 using ice sheet and glacier models for the latest set of socio-economic scenarios, nor using coordinated exploration of uncertainties arising from the various computer models involved. Two recent international projects generated a large suite of projections using multiple models2,3,4,5,6,7,8, but primarily used previous-generation scenarios9 and climate models10, and could not fully explore known uncertainties. Here we estimate probability distributions for these projections under the new scenarios11,12 using statistical emulation of the ice sheet and glacier models. We find that limiting global warming to 1.5 degrees Celsius would halve the land ice contribution to twenty-first-century sea level rise, relative to current emissions pledges. The median decreases from 25 to 13 centimetres sea level equivalent (SLE) by 2100, with glaciers responsible for half the sea level contribution. The projected Antarctic contribution does not show a clear response to the emissions scenario, owing to uncertainties in the competing processes of increasing ice loss and snowfall accumulation in a warming climate. However, under risk-averse (pessimistic) assumptions, Antarctic ice loss could be five times higher, increasing the median land ice contribution to 42 centimetres SLE under current policies and pledges, with the 95th percentile projection exceeding half a metre even under 1.5 degrees Celsius warming. This would severely limit the possibility of mitigating future coastal flooding. Given this large range (between 13 centimetres SLE using the main projections under 1.5 degrees Celsius warming and 42 centimetres SLE using risk-averse projections under current pledges), adaptation planning for twenty-first-century sea level rise must account for a factor-of-three uncertainty in the land ice contribution until climate policies and the Antarctic response are further constrained.

Description: Projections of the sea level contribution from the Greenland and Antarctic ice sheets (GrIS and AIS) rely on atmospheric and oceanic drivers obtained from climate models. The Earth System Models participating in the Coupled Model Intercomparison Project phase 6 (CMIP6) generally project greater future warming compared with the previous Coupled Model Intercomparison Project phase 5 (CMIP5) effort. Here we use four CMIP6 models and a selection of CMIP5 models to force multiple ice sheet models as part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). We find that the projected sea level contribution at 2100 from the ice sheet model ensemble under the CMIP6 scenarios falls within the CMIP5 range for the Antarctic ice sheet but is significantly increased for Greenland. Warmer atmosphere in CMIP6 models results in higher Greenland mass loss due to surface melt. For Antarctica, CMIP6 forcing is similar to CMIP5 and mass gain from increased snowfall counteracts increased loss due to ocean warming.

Description: 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.