ALBERT

Description: Until relatively recently, it was assumed that Arctic ice masses would respond to climatic/oceanic forcing over millennia, but observations made during the past two decades have radically altered this viewpoint and have demonstrated that marine-terminating outlet glaciers can undergo dramatic dynamic change at annual timescales. This paper reviews the substantial progress made in our understanding of the links between marine-terminating Arctic outlet glacier behaviour and the ocean-climate system during the past 20 years, when many ice masses have rapidly lost mass. Specifically, we assess three primary climatic/oceanic controls on outlet glacier dynamics, namely air temperature, ocean temperature and sea ice concentrations, and discuss key linkages between them. Despite recent progress, significant uncertainty remains over the response of marine-terminating outlet glaciers to these forcings, most notably: (1) the spatial variation in the relative importance of each factor; (2), the contribution of glacier-specific factors to glacier dynamics; and (3) the limitations in our ability to accurately model marine-terminating outlet glacier behaviour. Our present understanding precludes us from identifying patterns of outlet glacier response to forcing that are applicable across the Arctic and we underscore the potential danger of extrapolating rates of mass loss from a small sample of study glaciers.

Description: Variability in tidewater outlet glacier behavior under similar external forcing has been attributed to differences in outlet shape (i.e., bed elevation and width), but this dependence has not been investigated in detail. Here we use a numerical ice flow model to show that the dynamics of tidewater outlet glaciers under external forcing are highly sensitive to width and bed topography. Our sensitivity tests indicate that for glaciers with similar discharge, the trunks of wider glaciers and those grounded over deeper basal depressions tend to be closer to flotation, so that less dynamically induced thinning results in rapid, unstable retreat following a perturbation. The lag time between the onset of the perturbation and unstable retreat varies with outlet shape, which may help explain intra-regional variability in tidewater outlet glacier behavior. Further, because the perturbation response is dependent on the thickness relative to flotation, varying the bed topography within the range of observational uncertainty can result in either stable or unstable retreat due to the same perturbation. Thus, extreme care must be taken when interpreting the future behavior of actual glacier systems using numerical ice flow models that are not accompanied by comprehensive sensitivity analyses.

Description: Variability in tidewater outlet glacier behavior under similar external forcing has been attributed to differences in outlet shape (i.e. bed elevation and width), but this dependence has not been investigated in detail. Here we use a numerical ice flow model to show that the dynamics of tidewater outlet glaciers under external forcing are highly sensitive to width and bed topography. Our sensitivity tests indicate that for glaciers with similar discharge, the trunks of wider glaciers and those grounded over deeper basal depressions tend to be closer to flotation, so that less dynamically induced thinning results in rapid, unstable retreat following a perturbation. The lag time between the onset of the perturbation and unstable retreat varies with outlet shape, which may help explain intra-regional variability in tidewater outlet glacier behavior. Further, because the perturbation response is dependent on the thickness relative to flotation, varying the bed topography within the range of observational uncertainty can result in either stable or unstable retreat due to the same perturbation. Thus, extreme care must be taken when interpreting the future behavior of actual glacier systems using numerical ice flow models that are not accompanied by comprehensive sensitivity analyses.

Description: Predictions of marine ice-sheet behaviour require models that are able to robustly simulate grounding line migration. We present results of an intercomparison exercise for marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no effects of lateral buttressing). Unique steady-state grounding line positions exist for ice sheets on a downward sloping bed, while hysteresis occurs across an overdeepened bed, and stable steady state grounding line positions only occur on the downward-sloping sections. Models based on the shallow ice approximation, which does not resolve extensional stresses, do not reproduce the approximate analytical results unless appropriate parameterizations for ice flux are imposed at the grounding line. For extensional-stress resolving "shelfy stream" models, differences between model results were mainly due to the choice of spatial discretization. Moving grid methods were found to be the most accurate at capturing grounding line evolution, since they track the grounding line explicitly. Adaptive mesh refinement can further improve accuracy, including in fixed-grid models that generally perform poorly at coarse resolution. Fixed grid models with nested grid representations of the grounding line are able to generate accurate steady-state positions, but can be inaccurate over transients. Only one full Stokes model was included in the intercomparison, and consequently the accuracy of shelfy stream models as approximations of full Stokes models remains to be determined in detail, especially during transients.

Description: Depth-integrated (1-D) flowline models have been widely used to simulate fast-flowing tidewater glaciers and predict future change because their computational efficiency allows for continuous grounding line tracking, high horizontal resolution, and a physically-based calving criterion, which are all essential to realistic modeling of tidewater glaciers. As with all models, the values for parameters describing ice rheology and basal friction must be assumed and/or tuned based on observations. For prognostic studies, these parameters are typically tuned so that the glacier matches observed thickness and speeds at an initial state, to which a perturbation is applied. While it is well know that ice flow models are sensitive to these parameters, the sensitivity of tidewater glacier models has not been systematically investigated. Here we investigate the sensitivity of such flowline models of outlet glacier dynamics to uncertainty in three key parameters that influence a glacier's resistive stress components. We find that, within typical observational uncertainty, similar initial (i.e. steady-state) glacier configurations can be produced with substantially different combinations of parameter values, leading to differing transient responses after a perturbation is applied. In cases where the glacier is initially grounded near flotation across a basal overdeepening, as typically observed for rapidly changing glaciers, these differences can be dramatic owing to the threshold of stability imposed by the flotation criterion. The simulated transient response is particularly sensitive to the parameterization of ice rheology: differences in ice temperature of ∼ 2 °C can determine whether the glaciers thin to flotation and retreat unstably or remain grounded on a marine shoal. Due the highly non-linear dependence of tidewater glaciers on model parameters, we recommend that their predictions are accompanied by sensitivity tests that take parameter uncertainty into account.

Description: Depth-integrated (1-D) flowline models have been widely used to simulate fast-flowing tidewater glaciers and predict change because the continuous grounding line tracking, high horizontal resolution, and physically based calving criterion that are essential to realistic modeling of tidewater glaciers can easily be incorporated into the models while maintaining high computational efficiency. As with all models, the values for parameters describing ice rheology and basal friction must be assumed and/or tuned based on observations. For prognostic studies, these parameters are typically tuned so that the glacier matches observed thickness and speeds at an initial state, to which a perturbation is applied. While it is well know that ice flow models are sensitive to these parameters, the sensitivity of tidewater glacier models has not been systematically investigated. Here we investigate the sensitivity of such flowline models of outlet glacier dynamics to uncertainty in three key parameters that influence a glacier's resistive stress components. We find that, within typical observational uncertainty, similar initial (i.e., steady-state) glacier configurations can be produced with substantially different combinations of parameter values, leading to differing transient responses after a perturbation is applied. In cases where the glacier is initially grounded near flotation across a basal over-deepening, as typically observed for rapidly changing glaciers, these differences can be dramatic owing to the threshold of stability imposed by the flotation criterion. The simulated transient response is particularly sensitive to the parameterization of ice rheology: differences in ice temperature of ~ 2 °C can determine whether the glaciers thin to flotation and retreat unstably or remain grounded on a marine shoal. Due to the highly non-linear dependence of tidewater glaciers on model parameters, we recommend that their predictions are accompanied by sensitivity tests that take parameter uncertainty into account.

Description: Glacier calving can cause violent impulse waves which, upon landfall, can lead to destructive tsunami-like waves. Here we present data acquired during a calving event from Eqip Sermia, an ocean-terminating glacier in West Greenland. During an exceptionally well documented event, the collapse of 9 × 105 m3 ice from a 200 m high ice cliff caused an impulse wave of 50 m height, traveling at a speed of 25–30 m s-1. This wave was filmed from a tour boat in 800 m distance from the calving face, and simultaneously measured with a terrestrial radar interferometer and a tide gauge. Tsunami wave run-up height on the steep opposite shore in 4 km distance was 10–15 m, destroying infrastructure and eroding old vegetation. These observations indicate that such high tsunami waves are a recent phenomenon in the history of this glacier. Analysis of the data shows that only moderately bigger tsunami waves are to be expected in the future, even under rather extreme scenarios.

Description: Predictions of marine ice-sheet behaviour require models that are able to robustly simulate grounding line migration. We present results of an intercomparison exercise for marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no effects of lateral buttressing). Unique steady state grounding line positions exist for ice sheets on a downward sloping bed, while hysteresis occurs across an overdeepened bed, and stable steady state grounding line positions only occur on the downward-sloping sections. Models based on the shallow ice approximation, which does not resolve extensional stresses, do not reproduce the approximate analytical results unless appropriate parameterizations for ice flux are imposed at the grounding line. For extensional-stress resolving "shelfy stream" models, differences between model results were mainly due to the choice of spatial discretization. Moving grid methods were found to be the most accurate at capturing grounding line evolution, since they track the grounding line explicitly. Adaptive mesh refinement can further improve accuracy, including fixed grid models that generally perform poorly at coarse resolution. Fixed grid models, with nested grid representations of the grounding line, are able to generate accurate steady state positions, but can be inaccurate over transients. Only one full-Stokes model was included in the intercomparison, and consequently the accuracy of shelfy stream models as approximations of full-Stokes models remains to be determined in detail, especially during transients.