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: 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: Projections of the contribution of the Greenland ice sheet to sea level rise comprise uncertainties that arise from the imposed climate forcing and from the underlying mathematical and numerical description used by ice flow models. Here, we present a comparative modelling study with the models SICOPOLIS, using the shallow ice approximation (SIA) on a structured grid, and ISSM, using a higher-order (HO) approximation of the Stokes equation on an unstructured grid. Starting from a paleoclimatic spin-up produced by SICOPOLIS, the models are forced with two different, simplified warming scenarios based on RCP2.6 projections from climate models, which are in line with the limit of global warming negotiated for the Paris Agreement. ISSM/HO produces lower flow speeds at the glacier termini, but more acceleration in narrow outlet glaciers compared to SICOPOLIS/SIA. This leads to a larger elevation reduction for ISSM/HO, and thus a positive feedback on the surface mass balance (with that of ISSM/HO becoming ∼50 Gt a−1 more negative). Across the two models and scenarios, the projected mass loss by 2300 is ∼62–88 mm sea level equivalent.

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

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.