Tidewater glaciers have been observed to experience instantaneous, stepwise increases in velocity during iceberg-calving events due to a loss of resistive stresses. These changes in stress can ...potentially impact tidewater glacier stability by promoting additional calving and affecting the viscous delivery of ice to the terminus. Using flow models and perturbation theory, we demonstrate that calving events and subsequent terminus readvance produce quasi-periodic, sawtooth oscillations in stress that originate at the terminus and propagate upstream. The stress perturbations travel at speeds much greater than the glacier velocities and, for laterally resisted glaciers, rapidly decay within a few ice thickness of the terminus. Consequently, because terminus fluctuations due to individual calving events tend to be much higher frequency than climate variations, individual calving events have little direct impact on the viscous delivery of ice to the terminus. This suggests that the primary mechanism by which calving events can trigger instability is by causing fluctuations in stress that weaken the ice and lead to additional calving and sustained terminus retreat. Our results further demonstrate a stronger response to calving events in simulations that include the full stress tensor, highlighting the importance of accounting for higher order stresses when developing calving parameterizations.
Marine ice-cliff instability could accelerate ice loss from Antarctica, and according to some model predictions could potentially contribute >1 m of global mean sea level rise by 2100 at current ...emission rates. Regions with over-deepening basins >1 km in depth (e.g., the West Antarctic Ice Sheet) are particularly susceptible to this instability, as retreat could expose increasingly tall cliffs that could exceed ice stability thresholds. Here, we use a suite of high-fidelity glacier models to improve understanding of the modes through which ice cliffs can structurally fail and derive a conservative ice-cliff failure retreat rate parameterization for ice-sheet models. Our results highlight the respective roles of viscous deformation, shear-band formation, and brittle-tensile failure within marine ice-cliff instability. Calving rates increase non-linearly with cliff height, but runaway ice-cliff retreat can be inhibited by viscous flow and back force from iceberg mélange.
The simple calving laws currently used in ice-sheet models do not adequately reflect the complexity and diversity of calving processes. To be effective, calving laws must be grounded in a sound ...understanding of how calving actually works. Here, we develop a new strategy for formulating calving laws, using (a) the Helsinki Discrete Element Model (HiDEM) to explicitly model fracture and calving processes, and (b) the continuum model Elmer/Ice to identify critical stress states associated with HiDEM calving events. A range of observed calving processes emerges spontaneously from HiDEM in response to variations in ice-front buoyancy and the size of subaqueous undercuts. Calving driven by buoyancy and melt under-cutting is under-predicted by existing calving laws, but we show that the location and magnitude of HiDEM calving events can be predicted in Elmer/Ice from characteristic stress patterns. Our results open the way to developing calving laws that properly reflect the diversity of calving processes, and provide a framework for a unified theory of the calving process continuum.
We review the two main approaches to estimating sea level rise over the coming century: physically plausible models of reduced complexity that exploit statistical relationships between sea level and ...climate forcing, and more complex physics‐based models of the separate elements of the sea level budget. Previously, estimates of future sea level rise from semiempirical models were considerably larger than those from process‐based models. However, we show that the most recent estimates of sea level rise by 2100 using both methods have converged, but largely through increased contributions and uncertainties in process‐based model estimates of ice sheets mass loss. Hence, we focus in this paper on ice sheet flow as this has the largest potential to contribute to sea level rise. Progress has been made in ice dynamics, ice stream flow, grounding line migration, and integration of ice sheet models with high‐resolution climate models. Calving physics remains an important and difficult modeling issue. Mountain glaciers, numbering hundreds of thousands, must be modeled by extensive statistical extrapolation from a much smaller calibration data set. Rugged topography creates problems in process‐based mass balance simulations forced by regional climate models with resolutions 10–100 times larger than the glaciers. Semiempirical models balance increasing numbers of parameters with the choice of noise model for the observations to avoid overfitting the highly autocorrelated sea level data. All models face difficulty in separating out non‐climate‐driven sea level rise (e.g., groundwater extraction) and long‐term disequilibria in the present‐day cryosphere‐sea level system.
Key Points
We review the two approaches to estimating sea level rise by 2100
We give an overview of glacier, ice sheet and semiempirical modeling
We discuss limitations of state of the art ice sheet and semiempirical models
Full text
Available for:
FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
Abstract
We investigate the physical basis of the crevasse-depth (CD) calving law by analysing relationships between glaciological stresses and calving behaviour at Sermeq Kujalleq (Store Glacier), ...Greenland. Our observations and model simulations show that the glacier has a stable position defined by a compressive arch between lateral pinning points. Ice advance beyond the arch results in calving back to the stable position; conversely, if melt-undercutting forces the ice front behind the stable position, it readvances because ice velocities exceed subaqueous melt rates. This behaviour is typical of self-organising criticality, in which the stable ice-front position acts as an attractor between unstable super-critical and sub-critical regimes. This perspective provides strong support for a ‘position-law’ approach to modelling calving at Sermeq Kujalleq, because any calving ‘rate’ is simply a by-product of how quickly ice is delivered to the critical point. The CD calving law predicts ice-front position from the penetration of surface and basal crevasse fields, and accurately simulates super-critical calving back to the compressive arch and melt-driven calving into the sub-critical zone. The CD calving law reflects the glaciological controls on calving at Sermeq Kujalleq and exhibits considerable skill in simulating its mean position and seasonal fluctuations.
Predictions of marine ice-sheet behaviour require models able to simulate grounding-line migration. We present results of an intercomparison experiment for plan-view 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 buttressing effects from lateral drag). Perturbation experiments specifying spatial variation in basal sliding parameters permitted the evolution of curved grounding lines, generating buttressing effects. The experiments showed regions of compression and extensional flow across the grounding line, thereby invalidating the boundary layer theory. Steady-state grounding-line positions were found to be dependent on the level of physical model approximation. Resolving grounding lines requires inclusion of membrane stresses, a sufficiently small grid size (<500 m), or subgrid interpolation of the grounding line. The latter still requires nominal grid sizes of <5 km. For larger grid spacings, appropriate parameterizations for ice flux may be imposed at the grounding line, but the short-time transient behaviour is then incorrect and different from models that do not incorporate grounding-line parameterizations. The numerical error associated with predicting grounding-line motion can be reduced significantly below the errors associated with parameter ignorance and uncertainties in future scenarios.
We investigate the temporal evolution and spatial distribution of mass balance on the glacier Midtre Lovénbreen, Svalbard. Running a diagnostic high-resolution full-stress ice flow model with ...geometries obtained from five digital elevation models (DEMs) in the period 1962–2005, we compute velocity fields and linearly interpolated volume change of the glacier. We evaluate the kinematic free surface equation using these model outputs to solve the surface mass balance (SMB). Monitoring data on Midtre Lovénbreen allows model results to be compared with point measurements from the glacier over several decades. This method allows us to estimate the mass balance over the entire glacier surface, beyond the spatially limited field measurements, and to derive past SMB over an extended time period.
Abstract
Retreat of calving glaciers worldwide has contributed substantially to sea-level rise in recent decades. Mass loss by calving contributes significantly to the uncertainty of sea-level rise ...projections. At Bowdoin Glacier, Northwest Greenland, most calving occurs by a few large events resulting from kilometre-scale fractures forming parallel to the calving front. High-resolution terrestrial radar interferometry data of such an event reveal that crevasse opening is fastest at low tide and accelerates during the final 36 h before calving. Using the ice flow model Elmer/Ice, we identify the crevasse water level as a key driver of modelled opening rates. Sea water-level variations in the range of local tidal amplitude (1 m) can reproduce observed opening rate fluctuations, provided crevasse water level is at least 4 m above the low-tide sea level. The accelerated opening rates within the final 36 h before calving can be modelled by additional meltwater input into the crevasse, enhanced ice cliff undercutting by submarine melt, ice damage increase due to tidal cyclic fatigue, crevasse deepening or a combination of these processes. Our results highlight the influence of surface meltwater and tides on crevasse opening leading to major calving events at grounded tidewater glaciers such as Bowdoin.
Iceberg calving accounts for between 30 % and 60 % of net mass loss from
the Greenland Ice Sheet, which has intensified and is now the single largest
contributor to global sea level rise in the ...cryosphere. Changes to calving
rates and the dynamics of calving glaciers represent a significant
uncertainty in projections of future sea level rise. A growing body of
observational evidence suggests that calving glaciers respond rapidly to
regional environmental change, but predictive capacity is limited by the lack
of suitable models capable of simulating calving mechanisms realistically.
Here, we use a 3-D full-Stokes calving model to investigate the environmental
sensitivity of Store Glacier, a large outlet glacier in West Greenland. We
focus on two environmental processes: undercutting by submarine melting and
buttressing by ice mélange, and our results indicate that Store Glacier
is likely to be able to withstand moderate warming perturbations in which the
former is increased by 50 % and the latter reduced by 50 %. However,
severe perturbation with a doubling of submarine melt rates or a complete
loss of ice mélange destabilises the calving front in our model runs.
Furthermore, our analysis reveals that stress and fracture patterns at
Store's terminus are complex and varied, primarily due to the influence of
basal topography. Calving style and environmental sensitivity vary greatly,
with propagation of surface crevasses significantly influencing iceberg
production in the northern side, whereas basal crevasses dominate in the
south. Any future retreat is likely to be initiated in the southern side by a
combination of increased submarine melt rates in summer and reduced
mélange strength in winter. The lateral variability, as well as the
importance of rotational and bending forces at the terminus, underlines the
importance of using the 3-D full-Stokes stress solution when modelling
Greenland's calving glaciers.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
Computer models are necessary for understanding and predicting marine ice sheet behaviour. However, there is uncertainty over implementation of physical processes at the ice base, both for grounded ...and floating glacial ice. Here we implement several sliding relations in a marine ice sheet flow-line model accounting for all stress components and demonstrate that model resolution requirements are strongly dependent on both the choice of basal sliding relation and the spatial distribution of ice shelf basal melting.Sliding relations that reduce the magnitude of the step change in basal drag from grounded ice to floating ice (where basal drag is set to zero) show reduced dependence on resolution compared to a commonly used relation, in which basal drag is purely a power law function of basal ice velocity. Sliding relations in which basal drag goes smoothly to zero as the grounding line is approached from inland (due to a physically motivated incorporation of effective pressure at the bed) provide further reduction in resolution dependence.A similar issue is found with the imposition of basal melt under the floating part of the ice shelf: melt parameterisations that reduce the abruptness of change in basal melting from grounded ice (where basal melt is set to zero) to floating ice provide improved convergence with resolution compared to parameterisations in which high melt occurs adjacent to the grounding line.Thus physical processes, such as sub-glacial outflow (which could cause high melt near the grounding line), impact on capability to simulate marine ice sheets. If there exists an abrupt change across the grounding line in either basal drag or basal melting, then high resolution will be required to solve the problem. However, the plausible combination of a physical dependency of basal drag on effective pressure, and the possibility of low ice shelf basal melt rates next to the grounding line, may mean that some marine ice sheet systems can be reliably simulated at a coarser resolution than currently thought necessary.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK