Contrasting the modelled sensitivity of the Amundsen Sea Embayment ice streams

Present-day mass loss from the West Antarctic ice sheet is centred on the Amundsen Sea Embayment (ASE), primarily through ice streams, including Pine Island, Thwaites and Smith glaciers. To understand the differences in response of these ice streams, we ran a perturbed parameter ensemble, using a vertically-integrated ice flow model with adaptive mesh refinement. We generated 71 sets of three physical parameters (basal traction coefficient, ice viscosity stiffening factor and sub-shelf melt rate), which we used to simulate the ASE for 50 years. We also explored the effects of different bed geometries and basal sliding laws. The mean rate of sea-level rise across the ensemble of simulations is comparable with current observed rates for the ASE. We found evidence that grounding line dynamics are sensitive to features in the bed geometry: simulations using BedMap2 geometry resulted in a higher rate of sea-level rise than simulations using a rougher geometry, created using mass conservation. Modelled grounding-line retreat of all the three ice streams was sensitive to viscosity and basal traction, while the melt rate was more important in Pine Island and Smith glaciers, which flow through more confined ice shelves than Thwaites, which has a relatively unconfined shelf

Proper orthogonal decomposition of ice velocity identifies drivers of flow variability at Sermeq Kujalleq (Jakobshavn Isbr AE)

Abstract. The increasing volume and spatio-temporal resolution of satellite-derived ice velocity data has created new exploratory opportunities for the quantitative analysis of glacier dynamics. One potential technique, Proper Orthogonal Decomposition (POD), also known as Empirical Orthogonal Functions, has proven to be a powerful and flexible technique for revealing coherent structures in a wide variety of environmental flows. In this study we investigate the applicability of POD to an openly available TanDEM-X/TerraSAR-X derived ice velocity dataset from Sermeq Kujalleq (Jakobshavn Isbræ), Greenland. We find three dominant modes with annual periodicity that we argue are explained by glaciological processes. Mode 1 is interpreted as relating to the stress-reconfiguration at the glacier terminus, known to be an important control on the glacier’s dynamics. Modes 2 and 3 together relate to the development of the spatially heterogenous glacier hydrological system and are primarily driven by the pressurisation and efficiency of the subglacial hydrological system. During the melt season, variations in the velocity shown in Modes 2 and 3 are explained by the drainage of nearby supraglacial melt ponds, as identified with a Google Earth Engine MODIS dynamic thresholding technique. By isolating statistical structures within velocity datasets, and through their comparison to glaciological theory and complementary datasets POD indicates which glaciological processes are responsible for the changing bulk velocity signal, as observed from space. With the proliferation of optical and radar derived velocity products (e.g. MEaSUREs/ESA CCI/PROMICE) we suggest POD, and potentially other modal decomposition techniques, will become increasingly useful in future studies of ice dynamics. </jats:p

Experimental protocol for sea level projections from ISMIP6 stand-alone ice sheet models

Projection of the contribution of ice sheets to sea level change as part of the Coupled Model Intercomparison Project Phase 6 (CMIP6) takes the form of simulations from coupled ice sheet–climate models and stand-alone ice sheet models, overseen by the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). This paper describes the experimental setup for process-based sea level change projections to be performed with stand-alone Greenland and Antarctic ice sheet models in the context of ISMIP6. The ISMIP6 protocol relies on a suite of polar atmospheric and oceanic CMIP-based forcing for ice sheet models, in order to explore the uncertainty in projected sea level change due to future emissions scenarios, CMIP models, ice sheet models, and parameterizations for ice–ocean interactions. We describe here the approach taken for defining the suite of ISMIP6 stand-alone ice sheet simulations, document the experimental framework and implementation, and present an overview of the ISMIP6 forcing to be used by participating ice sheet modeling groups

Future Sea Level Change Under Coupled Model Intercomparison Project Phase 5 and Phase 6 Scenarios From the Greenland and Antarctic Ice Sheets

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

Modeling the Greenland Ice Sheet's Committed Contribution to Sea Level During the 21st Century

Mass loss from the Greenland Ice Sheet can be partitioned between surface mass balance and discharge due to ice dynamics through its marine-terminating outlet glaciers. A perturbation to a glacier terminus (e.g., a calving event) results in both an instantaneous response in velocity and mass loss and a diffusive response due to the evolution of ice thickness over time. This diffusive response means the total impact of a retreat event can take decades to be fully realized. Here we model the committed response of the Greenland Ice Sheet by applying perturbations to the marine-terminating glacier termini that represent recent observed changes, and simulating the response over the 21st century, while holding the climate forcing constant. The sensitivity of the ice sheet response to model parameter uncertainty is explored within an ensemble framework, and Gravity Recovery and Climate Experiment data is used to constrain the results using a Bayesian calibration approach. We find that the Greenland Ice Sheet's committed contribution to 21st century sea level rise is at least 33.5 [17.5 52.4] mm (25th and 75th percentiles in brackets), with at least 6 mm being attributable directly to terminus retreat that occurred between 2007 and 2015. The spread in our projections is driven by uncertainty in the basal friction coefficient. Our results complement the ISMIP6 Greenland projections, which report the ice sheet response to future forcing, excluding the background response. In this way, we can obtain estimates of Greenland's total contribution to sea level rise in 2100

Revisiting Antarctic ice loss due to marine ice cliff instability

Predictions for sea-level rise from Antarctica this century range from zero to over one metre. The highest are driven by the controversial ‘marine ice cliff instability’ (MICI) hypothesis, where coastal ice cliffs rapidly collapse after ice shelves disintegrate from surface and sub-shelf melting caused by global warming. Here we quantify ice sheet modelling uncertainties for the MICI study and show the probability distributions are skewed towards lower values (most likely value: 45 cm under very high greenhouse gas concentrations). However, MICI is not required to reproduce sea-level changes in the mid-Pliocene, Last Interglacial or 1992-2017, and without it the results agree with previous studies (all 95th percentiles are less than 43 cm). We therefore find previous interpretations of the MICI projections over-estimate sea-level rise this century. The hypothesis is still poorly-constrained: confidence in projections with MICI would require much greater diversity in models of ice shelf vulnerability and ice cliff collapse

Choice of observation type affects Bayesian calibration of Greenland Ice Sheet model simulations

Abstract. Determining reliable probability distributions for ice sheet mass change over the coming century is critical to refining uncertainties in sea-level rise projections. Bayesian calibration, a method for constraining projection uncertainty using observations, has been previously applied to ice sheet projections but the impact of the chosen observation type on the calibrated posterior probability distributions has not been quantified. Here, we perform three separate Bayesian calibrations to constrain uncertainty in Greenland Ice Sheet (GrIS) simulations of the committed mass loss in 2100 under the current climate, using observations of velocity change, dynamic ice thickness change, and mass change. Comparing the posterior probability distributions shows that the median ice sheet mass change can differ by 119 % for the particular model ensemble that we used, depending on the observation type used in the calibration. More importantly for risk-averse sea-level planning, posterior probabilities of high-end mass change scenarios are highly sensitive to the observation selected for calibration. Furthermore, we show that using mass change observations alone may result in model simulations that overestimate flow acceleration and underestimate dynamic thinning around the margin of the ice sheet. Finally, we look ahead and present ideas for ways to improve Bayesian calibration of ice sheet projections. </jats:p

Experimental protocol for sea level projections from ISMIP6 stand-alone ice sheet models

Projection of the contribution of ice sheets to sea level change as part of the Coupled Model Intercomparison Project Phase 6 (CMIP6) takes the form of simulations from coupled ice sheet-climate models and stand-alone ice sheet models, overseen by the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). This paper describes the experimental setup for process-based sea level change projections to be performed with stand-alone Greenland and Antarctic ice sheet models in the context of ISMIP6. The ISMIP6 protocol relies on a suite of polar atmospheric and oceanic CMIP-based forcing for ice sheet models, in order to explore the uncertainty in projected sea level change due to future emissions scenarios, CMIP models, ice sheet models, and parameterizations for ice-ocean interactions. We describe here the approach taken for defining the suite of ISMIP6 stand-alone ice sheet simulations, document the experimental framework and implementation, and present an overview of the ISMIP6 forcing to be used by participating ice sheet modeling groups

Understanding of Contemporary Regional Sea-Level Change and the Implications for the Future

Global sea level provides an important indicator of the state of the warming climate, but changes in regional sea level are most relevant for coastal communities around the world. With improvements to the sea-level observing system, the knowledge of regional sea-level change has advanced dramatically in recent years. Satellite measurements coupled with in situ observations have allowed for comprehensive study and improved understanding of the diverse set of drivers that lead to variations in sea level in space and time. Despite the advances, gaps in the understanding of contemporary sea-level change remain and inhibit the ability to predict how the relevant processes may lead to future change. These gaps arise in part due to the complexity of the linkages between the drivers of sea-level change. Here we review the individual processes which lead to sea-level change and then describe how they combine and vary regionally. The intent of the paper is to provide an overview of the current state of understanding of the processes that cause regional sea-level change and to identify and discuss limitations and uncertainty in our understanding of these processes. Areas where the lack of understanding or gaps in knowledge inhibit the ability to provide the needed information for comprehensive planning efforts are of particular focus. Finally, a goal of this paper is to highlight the role of the expanded sea-level observation network—particularly as related to satellite observations—in the improved scientific understanding of the contributors to regional sea-level change