Search

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.

Earlier large-scale Greenland ice sheet sea-level projections (e.g., those run during the ice2sea and SeaRISE initiatives) have shown that ice sheet initial conditions have a large effect on the projections and give rise to important uncertainties. The goal of the initMIP-Greenland intercomparison exercise is to compare, evaluate and improve the initialisation techniques used in the ice sheet modelling community and to estimate the associated uncertainties in modelled mass changes. initMIP-Greenland is the first in a series of ice sheet model intercomparison activities within ISMIP6 (the Ice Sheet Model Intercomparison Project for CMIP6), which is the primary activity within the Coupled Model Intercomparison Project - phase 6 (CMIP6) focusing on the ice sheets. Two experiments for the large-scale Greenland ice sheet have been designed to allow intercomparison between participating models of 1) the initial present-day state of the ice sheet and 2) the response in two idealised forward experiments. The forward experiments serve to evaluate the initialisation in terms of model drift (forward run without additional forcing) and in response to a large perturbation (prescribed surface mass balance anomaly), and should not be interpreted as sea-level projections. We present and discuss results that highlight the diversity of data sets, boundary conditions and initialisation techniques used in the community to generate initial states of the Greenland ice sheet. We find good agreement across the ensemble for the dynamic response to surface mass balance changes in areas where the simulated ice sheets overlap, but differences arising from the initial size of the ice sheet. The model drift in the control experiment is reduced for models that participated in earlier intercomparison exercises.

Estimating the future evolution of the Antarctic Ice Sheet (AIS) is critical for improving future sea level rise (SLR) projections. Numerical ice sheet models are invaluable tools for bounding Antarctic vulnerability; yet, few continental-scale projections of century-scale AIS SLR contribution exist, and those that do vary by up to an order of magnitude. This is partly because model projections of future sea level are inherently uncertain and depend largely on the model's boundary conditions and climate forcing, which themselves are unknown due to the uncertainty in the projections of future anthropogenic emissions and subsequent climate response. Here, we aim to improve the understanding of how uncertainties in model forcing and boundary conditions affect ice sheet model simulations. With use of sampling techniques embedded within the Ice Sheet System Model (ISSM) framework, we assess how uncertainties in snow accumulation, ocean-induced melting, ice viscosity, basal friction, bedrock elevation, and the presence of ice shelves impact continental-scale 100-year model simulations of AIS future sea level contribution. Overall, we find that AIS sea level contribution is strongly affected by grounding line retreat, which is driven by the magnitude of ice shelf basal melt rates and by variations in bedrock topography. In addition, we find that over 1.2 m of AIS global mean sea level contribution over the next century is achievable, but not likely, as it is tenable only in response to unrealistically large melt rates and continental ice shelf collapse. Regionally, we find that under our most extreme 100-year warming experiment generalized for the entire ice sheet, the Amundsen Sea sector is the most significant source of model uncertainty (1032 mm 6σ spread) and the region with the largest potential for future sea level contribution (297 mm). In contrast, under a more plausible forcing informed regionally by literature and model sensitivity studies, the Ronne basin has a greater potential for local increases in ice shelf basal melt rates. As a result, under this more likely realization, where warm waters reach the continental shelf under the Ronne ice shelf, it is the Ronne basin, particularly the Evans and Rutford ice streams, that are the greatest contributors to potential SLR (161 mm) and to simulation uncertainty (420 mm 6spread).

The increasing contribution of the Antarctic Ice Sheet to sea level rise is linked to reductions in ice shelf buttressing, driven in large part by basal melting of ice shelves. These ocean-driven buttressing losses are being compounded as ice shelves weaken and fracture. To date, model projections of ice sheet evolution have not accounted for weakening ice shelves. Here we present the first constitutive framework for ice deformation that explicitly includes mechanical weakening, based on observations of the progressive degradation of the remnant Larsen B Ice Shelf from 2000 to 2015. We implement this framework in an ice sheet model and are able to reproduce most of the observed weakening of the ice shelf. In addition to predicting ice shelf weakening and reduced buttressing, this new framework opens the door for improved understanding and predictions of iceberg calving, meltwater routing and hydrofracture, and ice shelf collapse.

The mechanical properties of glacier beds play a fundamental role in regulating the sensitivity of glaciers to environmental forcing across a wide range of timescales. Glaciers are commonly underlain by deformable till whose mechanical properties and influence on ice flow are not well understood but are critical for reliable projections of future glacier states. Using synoptic-scale observations of glacier motion in different seasons to constrain numerical ice flow models, we study the mechanics of the bed beneath Hofsjökull, a land-terminating ice cap in central Iceland. Our results indicate that the bed deforms plastically and weakens following incipient summertime surface melt. Combining the inferred basal shear traction fields with a Coulomb-plastic bed model, we estimate the spatially distributed effective basal water pressure and show that changes in basal water pressure and glacier accelerations are non-local and non-linear. These results motivate an idealized physical model relating mean basal water pressure and basal slip rate wherein the sensitivity of glacier flow to changes in basal water pressure is inversely related to the ice surface slope.

Calving is a major mechanism of ice discharge of the Antarctic and Greenland ice sheets, and a change in calving front position affects the entire stress regime of marine terminating glaciers. The representation of calving front dynamics in a 2-D or 3-D ice sheet model remains non-trivial. Here, we present the theoretical and technical framework for a level-set method, an implicit boundary tracking scheme, which we implement into the Ice Sheet System Model (ISSM). This scheme allows us to study the dynamic response of a drainage basin to user-defined calving rates. We apply the method to Jakobshavn Isbræ, a major marine terminating outlet glacier of the West Greenland Ice Sheet. The model robustly reproduces the high sensitivity of the glacier to calving, and we find that enhanced calving triggers significant acceleration of the ice stream. Upstream acceleration is sustained through a combination of mechanisms. However, both lateral stress and ice influx stabilize the ice stream. This study provides new insights into the ongoing changes occurring at Jakobshavn Isbræ and emphasizes that the incorporation of moving boundaries and dynamic lateral effects, not captured in flow-line models, is key for realistic model projections of sea level rise on centennial timescales.

The Ice Sheet System Model (ISSM) team, which is funded primarily by NASA (Cryosphere and Modeling, Analysis and Prediction programs), as well as the Jet Propulsion Laboratory, University of California, Irvine, and National Science Foundation's Office of Polar Programs, organized a workshop in collaboration with the Bjerknes Centre for Climate Research (BCCR) at the University of Bergen in Norway, in June 2014.

Ice flow numerical models are essential for predicting the evolution of ice sheets in a warming climate. Recent research emphasizes the need for higher-order and even full-Stokes flow models, instead of the traditional shallow-ice approximation, whose assumptions are not valid in certain critical areas. These higher-order models are, however, computationally intensive and difficult to use at the continental scale. Here we present a new technique, the Tiling method, to couple ice flow models of varying orders of complexity. The goal of the method is to limit the spatial extent of where higherorder models are applied to reduce the computational cost, while maintaining the model precision. We apply this method on synthetic geometries to demonstrate its practical use. We first use a geometry for which all models yield the same results to check the consistency of the method. Then we apply our method to a geometry for which a full-Stokes model is required in the vicinity of the ice front. Our results show that the hybrid models present significant improvements over mono-model approaches and reduce computational times.

Interferometric radar images collected by ERS-1, ERS-2 and RADARSAT-1 are used to observe the rupture tip of rifts that propagate along Hemmen Ice Rise on the Ronne Ice Shelf, Antarctica. Interferograms generated in 1992 and 1997 allow for the observation of ice deformation accumulated over 9 and 24 days respectively. These interferograms are combined, in order to separate the continuous process of creep deformation from the more cyclic motion caused by variations in ocean tide. An examination of local gradients in creep deformation reveals the pattern of ice deformation around and near the rupture tips and rifts with great precision (up to 10 cm a-1). We compare the observations with a deformation model for ice and obtain the following results: (1) The tidal oscillation of the Ronne Ice Shelf only yields small deformations along the rifts and near the rupture tips. (2) Along the ice front, the rifts and at the rupture tips, vertical bending is observed which is well explained by a model of viscous deformation of ice. Furthermore, the model indicates that the deformation pattern observed at the rupture tips is a sensitive indicator of the propagation state of the rifts (i.e. active vs inactive). (3) The viscous adjustment of ice is the dominant mode of deformation, masking the deformation pattern predicted by linear elastic fracture mechanics (LEFM). (4) Yet, at a spatial scale equivalent to the length of a rift, the propagation rate is well predicted by LEFM.