2022, Feldmann, Johannes, Reese, Ronja, Winkelmann, Ricarda, Levermann, Anders

Basal ice-shelf melting is the key driver of Antarctica's increasing sea-level contribution. In diminishing the buttressing force of the ice shelves that fringe the ice sheet, the melting increases the ice discharge into the ocean. Here we contrast the influence of basal melting in two different ice-shelf regions on the time-dependent response of an isothermal, inherently buttressed ice-sheet-shelf system. In the idealized numerical simulations, the basal-melt perturbations are applied close to the grounding line in the ice-shelf's (1) ice-stream region, where the ice shelf is fed by the fastest ice masses that stream through the upstream bed trough and (2) shear margins, where the ice flow is slower. The results show that melting below one or both of the shear margins can cause a decadal to centennial increase in ice discharge that is more than twice as large compared to a similar perturbation in the ice-stream region. We attribute this to the fact that melt-induced ice-shelf thinning in the central grounding-line region is attenuated very effectively by the fast flow of the central ice stream. In contrast, the much slower ice dynamics in the lateral shear margins of the ice shelf facilitate sustained ice-shelf thinning and thereby foster buttressing reduction. Regardless of the melt location, a higher melt concentration toward the grounding line generally goes along with a stronger response. Our results highlight the vulnerability of outlet glaciers to basal melting in stagnant, buttressing-relevant ice-shelf regions, a mechanism that may gain importance under future global warming.

2023, Feldmann, Johannes, Levermann, Anders

The timescales of the flow and retreat of Greenland's and Antarctica's outlet glaciers and their potential instabilities are arguably the largest uncertainty in future sea-level projections. Here we derive a scaling relation that allows the comparison of the timescales of observed complex ice flow fields with geometric similarity. The scaling relation is derived under the assumption of fast, laterally confined, geometrically similar outlet-glacier flow over a slippery bed, i.e., with negligible basal friction. According to the relation, the time scaling of the outlet flow is determined by the product of the inverse of (1) the fourth power of the width-To-length ratio of its confinement, (2) the third power of the confinement depth and (3) the temperature-dependent ice softness. For the outflow at the grounding line of streams with negligible basal friction, this means that the volume flux is proportional to the ice softness and the bed depth, but goes with the fourth power of the gradient of the bed and with the fifth power of the width of the stream. We show that the theoretically derived scaling relation is supported by the observed velocity scaling of outlet glaciers across Greenland as well as by idealized numerical simulations of marine ice-sheet instabilities (MISIs) as found in Antarctica. Assuming that changes in the ice-flow velocity due to ice-dynamic imbalance are proportional to the equilibrium velocity, we combine the scaling relation with a statistical analysis of the topography of 13 MISI-prone Antarctic outlets. Under these assumptions, the timescales in response to a potential destabilization are fastest for Thwaites Glacier in West Antarctica and Mellor, Ninnis and Cook Glaciers in East Antarctica; between 16 and 67 times faster than for Pine Island Glacier. While the applicability of our results is limited by several strong assumptions, the utilization and potential further development of the presented scaling approach may help to constrain timescale estimates of outlet-glacier flow, augmenting the commonly exploited and comparatively computationally expensive approach of numerical modeling.

2022, Schlemm, Tanja, Feldmann, Johannes, Winkelmann, Ricarda, Levermann, Anders

Owing to global warming and particularly high regional ocean warming, both Thwaites and Pine Island Glaciers in the Amundsen region of the Antarctic Ice Sheet could lose their buttressing ice shelves over time. We analyse the possible consequences using the parallel ice sheet model (PISM), applying a simple cliff-calving parameterization and an ice mélange-buttressing model. We find that the instantaneous loss of ice-shelf buttressing, due to enforced ice-shelf melting, initiates grounding-line retreat and triggers marine ice sheet instability (MISI). As a consequence, the grounding line progresses into the interior of the West Antarctic Ice Sheet and leads to a sea level contribution of 0.6m within 100a. By subjecting the exposed ice cliffs to cliff calving using our simplified parameterization, we also analyse marine ice cliff instability (MICI). In our simulations it can double or even triple the sea level contribution depending on the only loosely constrained parameter that determines the maximum cliff-calving rate. The speed of MICI depends on this upper bound of the calving rate, which is given by the ice mélange buttressing the glacier. However, stabilization of MICI may occur for geometric reasons. Because the embayment geometry changes as MICI advances into the interior of the ice sheet, the upper bound on calving rates is reduced and the progress of MICI is slowed down. Although we cannot claim that our simulations bear relevant quantitative estimates of the effect of ice-mélange buttressing on MICI, the mechanism has the potential to stop the instability. Further research is needed to evaluate its role for the past and future evolution of the Antarctic Ice Sheet.