Projecting Antarctic ice discharge using response functions from SeaRISE ice-sheet models

Abstract

Abstract. The largest uncertainty in projections of future sea-level change results from the potentially changing dynamical ice discharge from Antarctica. Basal ice-shelf melting induced by a warming ocean has been identified as a major cause for additional ice flow across the grounding line. Here we derive dynamic ice-sheet response functions for basal ice-shelf melting for four different Antarctic drainage regions using experiments from the Sea-level Response to Ice Sheet Evolution (SeaRISE) intercomparison project with five different Antarctic ice-sheet models. Under the assumptions of linear-response theory we project future ice-discharge for each model, each region and each of the four Representative Concentration Pathways (RCP) using oceanic temperatures from 19 comprehensive climate models of the Coupled Model Intercomparison Project, CMIP-5, and two ocean models from the EU-project Ice2Sea. The uncertainty in the climatic forcing, the oceanic response and the ice-model response is combined into an uncertainty range of future Antarctic ice-discharge induced from basal ice-shelf melt. The uncertainty range we derived for the Antarctic contribution to global sea-level rise from 1992 to 2011 is in full agreement with the observed contribution for this period if we use the three ice-sheet models with an explicit representation of ice-shelf dynamics and account for the time delayed warming of the oceanic subsurface compared with the surface air temperature. The median of the additional ice-loss for the 21st century (Table 6) is 0.07 m (66%-range: 0.02–0.14 m; 90%-range: 0.0–0.23 m) of global sea-level equivalent for the low-emission RCP-2.6 scenario and 0.09 m (66%-range: 0.04–0.21 m; 90%-range: 0.01–0.37 m) for the strongest RCP-8.5 if models with explicit ice-shelf representation are applied. These results were obtained using a time delay between the surface warming signal and the subsurface oceanic warming as observed in the CMIP-5 models. Without this time delay the values increase to 0.09 m (66%-range: 0.04–0.17 m; 90%-range: 0.02–0.25 m) for RCP-2.6 and 0.15 m (66%-range: 0.07–0.28 m; 90%-range: 0.04–0.43 m) for RCP-8.5. Our results are scenario dependent which is most visible in the upper percentiles of the distribution, i.e. highest contributions to sea level rise. All probability distributions, as provided in Fig. 12, are highly skewed towards high values. The applied ice-sheet models are coarse-resolution with limitations in the representation of grounding-line motion. However, we find the main uncertainty to be introduced by the external forcing to the ice-sheets, i.e. the climatic and oceanic uncertainty dominate. The scaling coefficients for the four different drainage basins provide valuable information for further assessments of future Antarctic ice discharge.