scholarly journals Full-Stokes modeling of grounding line dynamics, ice melt and iceberg calving for Thwaites Glacier, West Antarctica

2016 ◽  
Author(s):  
Hongju Yu ◽  
Eric Rignot ◽  
Mathieu Morlighem ◽  
Helene Seroussi

Abstract. Thwaites Glacier (TG), West Antarctica, has been losing mass and retreating rapidly in the past three decades. Here we present a two-dimensional, Full-Stokes (FS) modeling study of the grounding line dynamics and iceberg calving of TG. First, we compare FS with two simplified models, the higher-order (HO) model and the shallow-shelf approximation (SSA) model, to determine the impact of changes in ice shelf basal melt rate on grounding line dynamics. Second, we combine FS with the Linear Elastic Fracture Mechanics (LEFM) theory to simulate crevasse propagation and iceberg calving. In the first experiment, we find that FS requires basal melt rate consistent with remote sensing observations to reach steady state at TG’s current geometry while HO and SSA require unrealistically high basal melt rate. The grounding line of FS is also more sensitive to changes in basal melt rate than HO and SSA. In the second experiment, we find that only FS can produce surface and bottom crevasses that match radar sounding observations of crevasse width and height. We attribute the difference to the non- hydrostatic conditions of ice near the grounding line, which facilitate crevasse formation and are not accounted for in HO and SSA. Additional experiments using FS indicate that iceberg calving is significantly enhanced when surface crevasses exist near the grounding line, when ice shelf is shortened, or when the ice shelf front is undercut. We conclude that FS yields substantial improvements in the description of ice flow dynamics at the grounding line under high basal melt rate and in constraining crevasse formation and iceberg calving.

2016 ◽  
Author(s):  
Hongju Yu ◽  
Eric Rignot ◽  
Mathieu Morlighem ◽  
Helene Seroussi

Abstract. Thwaites Glacier (TG), West Antarctica, has been losing mass and retreating rapidly in the past few decades. Here, we present a study of its calving dynamics combining a two-dimensional flowband Full Stokes (FS) model of its viscous flow with linear elastic fracture mechanics (LEFM) theory to model crevasse propagation and ice fracturing. We compare the results with those obtained with the higher-order (HO) and the shallow-shelf approximation (SSA) models coupled with LEFM. We find that FS/LEFM produces surface and bottom crevasses that match the distribution of crevasse depth and width observed from NASA's Operation IceBridge radar depth sounders, whereas HO/LEFM and SSA/LEFM do not generate crevasses that match observations. We attribute the difference to the non-hydrostatic condition of ice near the grounding line, which facilitates crevasse formation, and is accounted for by the FS model but not by the HO or SSA model. We also find that calving is enhanced when pre-existing surface crevasses are present, when the ice shelf is shortened or when the ice shelf front is undercut. The role of undercutting depends on the time scale of calving events. It is more prominent for glaciers with rapid calving rates than glaciers with slow calving rates. Glaciers extending into a shorter ice shelf are more vulnerable to calving than glaciers developing a long ice shelf, especially as the ice front retreats close to the grounding line region, which leads to a positive feedback. We conclude that the FS/LEFM combination yields substantial improvements in capturing the stress field near the grounding line for constraining crevasse formation and iceberg calving.


2017 ◽  
Vol 11 (3) ◽  
pp. 1283-1296 ◽  
Author(s):  
Hongju Yu ◽  
Eric Rignot ◽  
Mathieu Morlighem ◽  
Helene Seroussi

Abstract. Thwaites Glacier (TG), West Antarctica, has been losing mass and retreating rapidly in the past few decades. Here, we present a study of its calving dynamics combining a two-dimensional flow-band full-Stokes (FS) model of its viscous flow with linear elastic fracture mechanics (LEFM) theory to model crevasse propagation and ice fracturing. We compare the results with those obtained with the higher-order (HO) and the shallow-shelf approximation (SSA) models coupled with LEFM. We find that FS/LEFM produces surface and bottom crevasses that are consistent with the distribution of depth and width of surface and bottom crevasses observed by NASA's Operation IceBridge radar depth sounder and laser altimeter, whereas HO/LEFM and SSA/LEFM do not generate crevasses that are consistent with observations. We attribute the difference to the nonhydrostatic condition of ice near the grounding line, which facilitates crevasse formation and is accounted for by the FS model but not by the HO or SSA models. We find that calving is enhanced when pre-existing surface crevasses are present, when the ice shelf is shortened or when the ice shelf front is undercut. The role of undercutting depends on the timescale of calving events. It is more prominent for glaciers with rapid calving rates than for glaciers with slow calving rates. Glaciers extending into a shorter ice shelf are more vulnerable to calving than glaciers developing a long ice shelf, especially as the ice front retreats close to the grounding line region, which leads to a positive feedback to calving events. We conclude that the FS/LEFM combination yields substantial improvements in capturing the stress field near the grounding line of a glacier for constraining crevasse formation and iceberg calving.


2011 ◽  
Vol 5 (4) ◽  
pp. 2035-2060 ◽  
Author(s):  
A. Luckman ◽  
D. Jansen ◽  
B. Kulessa ◽  
E. C. King ◽  
P. Sammonds ◽  
...  

Abstract. Basal crevasses extend upwards from the base of ice bodies and can penetrate more than halfway through the ice column under conditions found commonly on ice shelves. As a result, they may locally modify the exchange of mass and energy between ice shelf and ocean, and by altering the shelf's mechanical properties could play a fundamental role in ice shelf stability. Although early studies revealed that such features may be abundant on Antarctic ice shelves, their geometrical properties and spatial distribution has gained little attention. We investigate basal crevasses in Larsen C Ice Shelf using field radar survey, remote sensing and numerical modelling. We demonstrate that a group of features visible in MODIS imagery are the surface expressions of basal crevasses in the form of surface troughs, and find that basal crevasses can be generated as a result of stresses well downstream of the grounding line. We show that linear elastic fracture mechanics modelling is a good predictor of basal crevasse penetration height where stresses are predominantly tensile, and that measured surface trough depth does not always reflect this height, probably because of snow accumulation in the trough, marine ice accretion in the crevasse, or stress bridging from the surrounding ice. We conclude that all features visible in MODIS imagery of ice shelves and previously labelled simply as "crevasses", where they are not full thickness rifts, must be basal crevasse troughs, highlighting a fundamental structural property of many ice shelves that may have been previously overlooked.


2020 ◽  
Vol 14 (5) ◽  
pp. 1673-1683
Author(s):  
Bradley Paul Lipovsky

Abstract. Understanding the processes that govern ice shelf extent is important to improving estimates of future sea-level rise. In present-day Antarctica, ice shelf extent is most commonly determined by the propagation of through-cutting fractures called ice shelf rifts. Here, I present the first three-dimensional analysis of ice shelf rift propagation. I model rifts using the assumptions of linear elastic fracture mechanics (LEFM). The model predicts that rifts may be stabilized (i.e., stop propagating) when buoyant flexure results in the partial contact of rift walls. This stabilizing tendency may be overcome, however, by processes that act in the ice shelf margins. In particular, loss of marginal strength, modeled as a transition from zero tangential displacement to zero tangential shear stress, is shown to favor rift propagation. Rift propagation may also be triggered if a rift is carried with the ice flow (i.e., advected) out of an embayment and into a floating ice tongue. I show that rift stability is closely related to the transition from uniaxial to biaxial extension known as the compressive arch. Although the partial contact of rift walls is fundamentally a three-dimensional process, I demonstrate that it may be parameterized within more numerically efficient two-dimensional calculations. This study constitutes a step towards a first-principle description of iceberg calving due to ice shelf rift propagation.


2021 ◽  
Vol 15 (1) ◽  
pp. 113-132
Author(s):  
Jan De Rydt ◽  
Ronja Reese ◽  
Fernando S. Paolo ◽  
G. Hilmar Gudmundsson

Abstract. Pine Island Glacier in West Antarctica is among the fastest changing glaciers worldwide. Over the last 2 decades, the glacier has lost in excess of a trillion tons of ice, or the equivalent of 3 mm of sea level rise. The ongoing changes are thought to have been triggered by ocean-induced thinning of its floating ice shelf, grounding line retreat, and the associated reduction in buttressing forces. However, other drivers of change, such as large-scale calving and changes in ice rheology and basal slipperiness, could play a vital, yet unquantified, role in controlling the ongoing and future evolution of the glacier. In addition, recent studies have shown that mechanical properties of the bed are key to explaining the observed speed-up. Here we used a combination of the latest remote sensing datasets between 1996 and 2016, data assimilation tools, and numerical perturbation experiments to quantify the relative importance of all processes in driving the recent changes in Pine Island Glacier dynamics. We show that (1) calving and ice shelf thinning have caused a comparable reduction in ice shelf buttressing over the past 2 decades; that (2) simulated changes in ice flow over a viscously deforming bed are only compatible with observations if large and widespread changes in ice viscosity and/or basal slipperiness are taken into account; and that (3) a spatially varying, predominantly plastic bed rheology can closely reproduce observed changes in flow without marked variations in ice-internal and basal properties. Our results demonstrate that, in addition to its evolving ice thickness, calving processes and a heterogeneous bed rheology play a key role in the contemporary evolution of Pine Island Glacier.


2004 ◽  
Vol 39 ◽  
pp. 133-138 ◽  
Author(s):  
Robert Thomas ◽  
Eric Rignot ◽  
Pannirselvam Kanagaratnam ◽  
William Krabill ◽  
Gino Casassa

AbstractPine Island Glacier, flowing into the Amundsen Sea from West Antarctica, thinned substantially during the 1990s, its grounding line receded by several km, and its velocity increased by >10% to values approaching 3 km a–1. Here, we use these observations, together with estimates of ice thickness and surface strain rates, to estimate the perturbation in forces resisting ice flow compatible with the observations. The analysis assumes that such perturbations are transmitted far upstream from where they originate, and that creep response to the perturbations can be described by equations similar to those that govern ice-shelf creep. It indicates that observed acceleration between 1996 and 2000 could have been caused by progressive ungrounding within the most seaward 25 km ‘ice plain’ of the grounded glacier. Earlier retreat and thinning of the glacier’s floating ice shelf may have provided the conditions that initiated ungrounding of the ice plain. Our analysis indicates that continued ice-plain thinning at the current rate of about 2 ma–1 will result in a velocity increase by 1 km a–1 within the next 11 years as the ice plain becomes totally ungrounded.


1998 ◽  
Vol 27 ◽  
pp. 182-186 ◽  
Author(s):  
Christina L. Hulbe ◽  
Eric Rignot ◽  
Douglas R. Macayeal

Comparison between numerical model ice-shelf flow simulations and synthetic aperture radar (SAR) interferograms is used to study ice-flow dynamics at the Hemmen Ice Rise (HIR) and Lassiter Coast (LC) corners of the iceberg-calving front of the Filchnei—Ronne Ice Shelf, Antarctica. The interferograms are constructed from SAR images provided by the European Space Agency's remote-sensing satellites (ERS-1/2). Narrow bands of large shear strain rate are observed along the boundaries between fast-flowing ice-shelf ice and no-flow boundaries. Large rifts, opened where the ice shelf separates from the coast, appear to be filled with a melange of sea ice, ice-shelf fragments, and snow. Trial and error is used to find the best match between artificial interferograms, constructed from modelled ice flow, and the observed interferograms. We find that at both HIR and LC, ice with in the coastal boundary layers must be significantly softer than adjacent ice. At HIR the rift-filling ice melange transmits stress from one ice-shelf fragment to another; thus it must have mechanical competence and must moderate both separation of the ice shelf from the coast and the release of icebergs. However, the ice melange along the LC does not. The difference may be related to melange thickness, which could vary in the two locations due to differences in sub-ice-shelf oceanography or perhaps to regional atmospheric warming, currently under way along the Antarctic Peninsula. Future warming could weaken the melange ice around HIR as well, causing the ice shelf to lose contact with that shelf-front anchor.


Sign in / Sign up

Export Citation Format

Share Document