the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Timescales of outletglacier flow with negligible basal friction: theory, observations and modeling
Johannes Feldmann
Anders Levermann
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 sealevel 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 outletglacier 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 widthtolength ratio of its confinement, (2) the third power of the confinement depth and (3) the temperaturedependent 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 icesheet instabilities (MISIs) as found in Antarctica. Assuming that changes in the iceflow velocity due to icedynamic imbalance are proportional to the equilibrium velocity, we combine the scaling relation with a statistical analysis of the topography of 13 MISIprone 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 outletglacier flow, augmenting the commonly exploited and comparatively computationally expensive approach of numerical modeling.
 Article
(5753 KB) 
Supplement
(68044 KB)  BibTeX
 EndNote
Global sealevel rise is increasingly driven by the contributions of the ice sheets of Greenland and Antarctica. Both ice sheets are observed to be out of balance (The IMBIE team, 2018; Rignot et al., 2019; Mouginot et al., 2019a; The IMBIE Team, 2020; Diener et al., 2021). This imbalance is mainly due to the increased discharge of ice into the ocean, which accounts for more than half of the cumulative net mass loss from the Greenland Ice Sheet and virtually all of the loss from the Antarctic Ice Sheet. Their discharge is projected to raise global sea level by several decimeters within the 21st century (Bamber et al., 2019; Levermann et al., 2020; Pattyn and Morlighem, 2020; Edwards et al., 2021; DeConto et al., 2021). Fed by snowfall, the ice sheets spread under their own weight towards the ocean. The flow of the ice accelerates from the interior to the coast, where the ice sheets are drained by numerous outlet glaciers that often stream through topographic confinements, e.g., valleyshaped bed troughs, before they discharge into the ocean (Morlighem et al., 2017, 2020; Maier et al., 2019). Such confinements tend to induce stabilizing iceinternal stresses that are transmitted upstream, exerting a buttressing effect on the hinterland (MacAyeal, 1989; Dupont and Alley, 2005; Schoof, 2007; Gudmundsson et al., 2012; Haseloff and Sergienko, 2018; Pegler, 2018), regulating the flow speed and thus the rate of ice discharge into the ocean. Note that we refer to the term “buttressing” as a backforce induced by laterally confined ice flow that is either floating (often referred to as “iceshelf buttressing”) or grounded.
In both ice sheets, the ice discharge mainly takes place through fast flowing outlet glaciers like Jakobshavn, Helheim or Kangerlussuaq Glaciers in Greenland (Catania et al., 2020; The IMBIE Team, 2020) and Pine Island or Thwaites Glaciers in West Antarctica (Mouginot et al., 2014; The IMBIE team, 2018; Diener et al., 2021). In Antarctica, large parts of the ice are grounded below sea level with a bed topography that is sloping down in the inland direction (retrograde slope; Pattyn and Morlighem, 2020; Morlighem et al., 2020). These two preconditions make the ice sheet prone to the socalled marine ice sheet instability (MISI), a positive feedback between ice discharge and icesheet retreat down the retrograde bed slope (Weertman, 1974; Schoof, 2007; Pattyn and Morlighem, 2020). The buttressing effect of laterally confined ice flow can slow down or even prevent MISItype destabilization (Gudmundsson et al., 2012; Haseloff and Sergienko, 2018; Pegler, 2018). The development of a MISI was hypothesized for the West Antarctic Ice Sheet more than 40 years ago (Hughes, 1973; Mercer, 1978) and meanwhile could indeed have been initiated, as suggested by observations and icedynamic modeling (Favier et al., 2014; Joughin et al., 2014; Seroussi et al., 2017; The IMBIE Team, 2020; Diener et al., 2021).
Determining the quantitative influence of buttressing on the timescales of confined ice flow and the related discharge, in general, as well as on MISI timescales, in particular, requires detailed numerical computations (Favier et al., 2014; Joughin et al., 2014; Pattyn and Morlighem, 2020). By comparing similar ice dynamics, it is however possible to derive scaling relations that comprise the complex stress and velocity field without the need of computationally expensive simulations (Feldmann and Levermann, 2016; Levermann and Feldmann, 2019). Such relations provide a link between the characteristic scales of a glacier, e.g., between its iceflow timescale and its geometric dimensions. While scaling relations are valid for each glacier separately, they also allow the comparison of similar glaciers, e.g., the comparison between the timescales of two glaciers based on their individual geometric scales.
To derive the scaling relation presented here, we conduct a scaling analysis of the governing dynamic equations of fast, shallow and laterally confined ice flow under the assumption of small basal friction. To this end, we employ the concept of similitude (Rayleigh, 1915), which is a common and simple procedure in hydrodynamics and engineering to investigate the scaling behavior of physical systems (Szücs, 1980). A famous example for its application is the Reynolds number for the Navier–Stokes equation (Reynolds, 1883; Kundu et al., 2012). A fixed Reynolds number provides a scaling relation that ensures similar flow patterns of a fluid under variation of its characteristic properties. The scaling relation derived in the present study links the timescale of the flow of an outlet glacier to its geometry, i.e., to the spatial dimensions of its lateral confinement, and to the temperature of the ice. Earlier work that dealt with the topic of confined outlet glacier flow includes analyses of the dependency of ice flux and groundingline dynamics on the confinement width and indepth investigations of the role of lateral shear margins (e.g., Nye, 1965; Raymond, 1996; Joughin et al., 2004; Gudmundsson et al., 2012; Pegler et al., 2013; Haseloff et al., 2015; Schoof et al., 2017; Pegler, 2018). Laterally unconfined ice flow and its response time to small perturbations have been examined for a regime of low basal shear and driving stress in Sergienko and Wingham (2019).
In a previous study, we applied the similitude principle to the onedimensional flow problem in combination with observational data to compare the timescales of geometrically similar, MISIprone Antarctic outlet glaciers (Levermann and Feldmann, 2019). This onedimensional approach, however, neglected the important buttressing effect of the outlet's lateral confinement. Our present analysis respects the buttressing effect of the grounded ice portion through the width of the confinement as one of the dominant parameters determining the timescale of the ice dynamics. The derivation of the associated scaling relation is outlined in Sect. 2. It is then tested in Sect. 3, examining the scaling of (1) observed Greenland outletglacier flow speeds and (2) simulated Antarctictype MISI retreat rates. In Sect. 4, we apply the scaling relation to a set of MISIprone Antarctic outlet glaciers to compare their potential retreat timescales, making the assumption that the retreat speed of an outlet glacier is proportional to its flow speed. The method we use to obtain the outletspecific timescales is principally similar to but more sophisticated than the one used in Levermann and Feldmann (2019), as we here extend it to both horizontal dimensions and carry out a more profound statistical analysis of the obtained outletspecific scales. The limitations of our approach are discussed in Sect. 5 where we also draw our conclusions.
Our derivation of the scaling relation (detailed in Appendix A1) is based on the shallowshelf approximation (SSA) of the Stokes stress balance of the ice (MacAyeal, 1989; Bueler and Brown, 2009). The SSA captures the dynamics of a fast and shallow ice flow and neglects vertical shearing inside the ice. It describes a flow which is controlled by the balance between horizontal iceinternal stresses, basal shear stresses between ice and bed, and the driving stress due to the surface gradient of the ice. Many of the outlet glaciers around Greenland and Antarctica are ice streams (Joughin et al., 2018; Mouginot et al., 2019a) that slide over a slippery bed (Morlighem et al., 2013; Maier et al., 2019). While the fast ice streaming can generally produce high basal shear stresses, it has been shown that in the well lubricated, marine groundingline regions, the outlets are typically characterized by comparatively low basal shear (Joughin et al., 2004; Sergienko and Hindmarsh, 2013; Sergienko et al., 2014). The lateral confinements of the ice streams can induce high lateral iceinternal shear stresses with a strong influence on the ice dynamics, especially if the ice flows through deep and narrow bed troughs (Dupont and Alley, 2005; Gudmundsson et al., 2012; Leguy, 2015). Here we consider the basal stresses in the groundingline region to be small compared to such iceinternal stresses and thus neglect them in the stress balance.
Typically, the ice streams in Greenland and Antarctica are much longer than wide, thus exhibiting a small widthtolength ratio, $R=W/L$, as the characteristic alongflow length scale L is typically observed to be a multiple of the width scale W (Morlighem et al., 2017; Mouginot et al., 2019a; Morlighem et al., 2020). Without loss of generality we place the two horizontal axes of the coordinate system of our analysis into the main flow direction (x) and perpendicular (y) to it (Fig. 1). Using this information in a simple dimensional analysis of the SSA, in which all quantities have been nondimensionalized (Appendix A1), the stress balance reduces to
with
where v_{x} is the ice speed in x direction and thereby along the bed trough, H is the ice thickness, h represents the surface elevation of the ice, n is Glen's flow law exponent determining the ice rheology (Glen, 1955), and ∂_{x} and ∂_{y} denote the spatial derivatives in the x and y directions. Notably, Eq. (1) includes a single dimensionless coefficient, λ, in which A is the ice softness, D and T represent the characteristic vertical length scale and the timescale of the system, respectively, ρ denotes the ice density and g is the gravitational acceleration. The coefficient λ represent the balance between (1) the leadingorder iceinternal stress, i.e., the horizontal shearing of the ice flow in acrosstrough direction and (2) the driving stress of the ice, which is determined by the surface gradient along the trough. As long as λ remains constant, the ratio of these two stresses will remain the same, while the characteristic scales of the system are allowed to vary. In other words, the ice dynamics of the considered system will remain similar under a scaling of the system that leaves the system's governing dimensionless number λ unchanged (concept of similitude). Considering a reference system and applying a scaling to it (denoted by a dash), the requirement of ${\mathit{\lambda}}^{\prime}=\mathit{\lambda}$ then yields a time–scaling relation that assures the similitude of the two systems as follows:
Here, $\mathit{\tau}={T}^{\prime}/T$ is the ratio of the timescales, $\mathit{\alpha}={A}^{\prime}/A$ denotes the scaling of the ice softness, and $\mathit{\beta}={D}^{\prime}/D$, $\mathit{\omega}={W}^{\prime}/W$ and $\mathit{\delta}={L}^{\prime}/L$ represent the scaling of the confinement depth, width and length, respectively (Fig. 1, Table 1). Laboratory studies yield a range of possible values for n ranging from 2 to 4 (e.g., Duval, 2013). Using the most commonly applied value of n=3, we obtain
where $\mathit{\gamma}=\mathit{\omega}/\mathit{\delta}={R}^{\prime}/R$ expresses the scaling of the horizontal aspect ratio (widthtolength ratio) of the confinement. The timescale is thus a strongly nonlinear function of the spatial dimensions of the confinement, i.e., decreasing with the fourth power of its widthtolength ratio and with the third power of the bed depth. Furthermore, the timescale is inversely related to the ice softness, which means that higher ice temperatures lead to a faster timescale. This is due to the fact that a higher softness enhances the deformation rate of the ice (Glen, 1955), reducing the restraint due to lateral shear and promoting alongflow extension of the ice, both inducing faster ice flow and discharge. The different characteristics of the geometry will generally have influences along the following lines: a deeper bed trough might imply a stronger ice gradient at the grounding line for similar glaciers, and thereby a stronger driving stress. In case of a MISItype retreat, a shorter horizontal length scale, i.e., a steeper retrograde bed slope, leads to a higher flux increase for each increment of groundingline retreat down the slope, hence amplifying the positive feedback loop between ice discharge and retreat (Schoof, 2007; Pattyn and Morlighem, 2020). A wider confinement diminishes the buttressing strength of the system (Dupont and Alley, 2005; Gudmundsson et al., 2012; Leguy, 2015), which in turn provides less moderation to the unstable retreat.
Using the confinement geometry, the timescaling τ (Eq. 3) can be translated into a scaling of the glacier outflow μ, i.e., the discharge through the confinement. For this purpose, the velocity scaling $\mathit{\nu}={V}^{\prime}/V=\mathit{\delta}/\mathit{\tau}$ (based on the velocity scale $V=L/T$) is multiplied with the scaling of the confinement depth β and width ω, respectively, yielding
This serves as a measure for the relative outflow between individual outlets, i.e., the average velocity times the vertical outflow area near the grounding line. Introducing the ratio of the depth gradients, $\mathrm{\Gamma}\equiv \mathit{\beta}/\mathit{\delta}$, this becomes
i.e., the outgoing volume flux scales linearly with the ice softness and the bed depth, respectively, and with the fourth power of the gradient of the bed topography and with the fifth power of the confinement width. Please note that the above relations only hold under the listed assumptions, especially the assumption of geometric similarity, i.e., for glaciers that have qualitatively the same topography only in different sizes.
The analytically derived timescaling relation (Eq. 4) is tested in two very different settings. First, we analyze observational geometric and velocity data from strongly confined outlet glaciers in Greenland. Second, we investigate the time scaling of simulated Antarctictype MISI.
3.1 Observed velocity scaling of Greenland outlet glaciers
The Greenland Ice Sheet drains through many very narrow, fjordtype outlets, suggesting a particularly strong control of the bed topography on the iceflow timescale (Morlighem et al., 2017; Catania et al., 2020). We use a dataset from Beckmann et al. (2019) which covers 12 outlet glaciers around Greenland (Figs. 2 and S1 in the Supplement) to compare the observed relationship between the iceflow timescale and the confinement width of the outlets to the results predicted by theory (Eq. 3). The comparison with the theory is only valid for observed glaciers that are approximately dynamically and geometrically similar, which excludes two glaciers from the analysis (see Appendix B for details). Observations allow us to directly use the velocity scale $V=L/T$ as a timescale measure, with velocity scaling $\mathit{\nu}={V}^{\prime}/V=\mathit{\delta}/\mathit{\tau}$. The geometry dependence of the velocity scaling according to Eq. (3) reduces to
for n=3 if differences in the width of the glaciers dominate the velocity field ($\mathit{\alpha}=\mathit{\beta}=\mathit{\delta}=\mathrm{1}$). This quartic relationship turns out to approximate the observed velocitywidth scaling reasonably well (Fig. 3). Deviations from the predicted curve likely occur due to the fact that the comparison of the complex realworld glaciers breaks the underlying theoretical assumption of perfect similarity. Deviations also result from differences between the individual glaciers in their ice softness fields, in the degree of nonlinearity of their ice rheologies (flowlaw exponent n) or in the vertical and alongflow geometry of their topographic confinements, which here are all assumed to be identical.
3.2 Modelled MISI time scaling with idealized geometry
The scaling relation (Eq. 4) applies to the timescale of glacier flow, as supported by the comparison with observations of Greenland outlet glaciers of similar geometry in the previous section. This flow timescale is linked to the ice discharge across the grounding line, since a faster coastal flow, i.e., a smaller τ, implies a larger groundingline discharge μ (Eq. 5). Groundingline retreat depends on the divergence of groundingline discharge, which involves the flow speed at the grounding line. If we were to seek a relation for the groundingline retreat rate, we could make the assumption that the changes in the iceflow speed in the course of retreat are proportional the equilibrium iceflow speed. This is a strong assumption and it might not be justified in all cases. Under this assumption, the timescale of glacier retreat would be proportional to the timescale of glacier flow (see Appendix A2 for the mathematical examination). This would mean that Eq. (4) would also be applicable to the time scaling of glacier retreat.
To test this, we simulate MISItype outletglacier retreat timescales under a systematic variation of each of the variables involved in Eq. (4) in an ensemble of idealized, threedimensional numerical experiments. For this purpose we employ the Parallel Ice Sheet Model (PISM; Bueler and Brown, 2009; Winkelmann et al., 2011; PISM authors, 2021). The model solves the full SSA equations (Eq. A1) of the stress balance from which the above scaling relation is derived, thus accounting for all the horizontal compressive, tensile and shear stresses that act within the ice as well as basal shearing.
The bed topography prescribed in the simulations is obtained and modified from the MISMIP+ intercomparison project, which is the current benchmark for idealized simulations of a buttressed marine icesheet–shelf system (AsayDavis et al., 2016; Cornford et al., 2020; Lhermitte et al., 2020; Zhang et al., 2020; Leguy et al., 2021; Feldmann et al., 2022). It shapes a trough with a retrograde bed section (Fig. 4), resembling the geometric characteristics of a deeply incised, MISIprone marine Antarctic outlet such as Pine Island Glacier. Ice flow in this setup mainly takes place from the interior through the bed trough, where an ice stream feeds a bayshaped ice shelf that calves into the ocean (Fig. S3). There is also a perpendicular flow component, i.e., from the channel's lateral ridges down into the trough. Resulting convergent flow and associated horizontal shearing (Fig. S4) enable the emergence of buttressing. Details on the model setup and the experimental design are given in Appendix C.
Cornford et al. (2020)Cornford et al. (2020)Cornford et al. (2020)In our simulations, we vary separately the four relevant parameters which are found to determine the timescale of the ice dynamics according to Eq. (4). That is, the depth D, the width W and the length L of the bed trough, and the ice temperature T_{ice} (ice softness A(T_{ice}); see Eq. C1) are altered one at a time (see Table 2 for parameter ranges). The timescale T is measured as the time per km of groundingline retreat (inverse retreat speed) down the central part of the retrograde bed section (marked grey in Fig. 5). This way, we obtain four sets of simulations, i.e., for the icesoftness scaling α, the confinementdepth scaling β, the length scaling δ and the width scaling ω, each of them yielding a curve for the timescaling τ (large panels in Fig. 6). To examine the scaling behavior of the measured retreat timescales, we normalize each set of simulations with respect to the same reference softness/confinement depth/confinement length/confinement width and the associated reference timescale (blue axes in Fig. 6), and compare the data to the theoretically expected scaling behavior according to Eq. (4) (blue curves in Fig. 6). For each scaling parameter, the numerical results are in good agreement with our analytically derived timescaling relation of Eq. (4). Linear regression of the data in doublelogarithmic plots (insets in Fig. 6) yields that the scaling exponents of the geometric scales D, W and L inferred from our simulations slightly underestimate the theoretical values (deviation by ≈13 %). Note that the theoretical scaling relation assumes negligible basal stresses, while in our simulations the basal stresses are not negligible (due to the prescribed nonnegligible value of the basal friction coefficient C in the sliding law), which is likely one reason for the deviations between simulations and theory. Additional simulations that vary the friction coefficient over 1 order of magnitude show that the retreat speed increases with decreasing friction with the time scaling, following Eq. (4) to a good approximation (Fig. S5). Deviations from theory increase with decreasing friction, since reducing C flattens the originally parabolic glacier profile from the reference simulation, leading to a violation of the similarity principle.
An additional set of experiments of altered accumulation rate a is conducted to analyze how a scaling of the surface mass balance affects the time scaling of grounding line movement. The simulations yield a slightly larger or smaller ice sheet depending on whether the accumulation rate is increased or decreased. This in turn slightly changes the depth, length and width measurements of our analysis, and thereby the timescale. Figure S5 shows that the changes induced by the relative accumulation changes have a relatively weak influence on the timescale.
Our analysis of observed Greenland outletglacier velocities revealed that their scaling follows the theoretical prediction from the scaling relation (Sect. 3.1; Fig. 3). The validity of the scaling relation is further supported by the conducted ensemble of idealized, numerical simulations of Antarctictype outletglacier retreat (Sect. 3.2), as simulated MISI retreat timescales are in good agreement with the ones predicted from theory (Fig. 6). This leads us to apply the theory to observations of Antarctica's MISIprone outlet glaciers to compare their retreat timescales after a potential destabilization. We base our approach on the fact that the currently observed geometry and physical conditions of the considered outlet glaciers represent the balance of the relevant forces that act on the ice. This static situation carries the information for their initial retreat timescale after a potential destabilization.
Analyzing the major outlets around Antarctica, we find 13 glaciers that show a clear retrograde bed slope upstream of their grounding line and have a generally similar bed profile. Please refer to the Figs. S9–S23 for maps of the topographic features of each glacier, and Appendix D1 for a definition of the applied similarity criteria. These glaciers are situated in the Amundsen Sea Sector, East Antarctica and the Filchner–Ronne region (Fig. 7). The three geometric scales required by Eq. (4), i.e., the depth, width and length scales, are obtained from the recent BedMachine Antarctica compilation (Morlighem et al., 2020). The ice softness (Fig. S9), entering the scaling relation as the fourth parameter, is obtained from a presentday simulation of the Antarctic Ice Sheet, which was conducted in the course of the recent ISMIP 6 Antarctica intercomparison project (Seroussi et al., 2020). Pine Island Glacier (PIG) is selected as the scaling reference for our analysis as one of the most studied glaciers in Antarctica. The scales are extracted along the streamlines of each outlet (Figs. S10–S24, panels a, c and d), based on the MEaSUREs Antarctic velocity map from Mouginot et al. (2019a). From this, we obtain probability distributions of the respective scaling ratios (Figs. S10–S24e), representing the spatial variability of the extracted data. The distributions of the retreat time scaling (Fig. 8) are calculated according to Eq. (4) using a Monte Carlo method. The medians of these distributions, $<\mathit{\tau}>$, represent the best estimate of our calculations. Details on our method and the underlying datasets are given in Appendix D.
The obtained timescales range over 3 orders of magnitude (Fig. 7 and Table 3), mainly due to differences in the bedslope magnitude (length scale), the width of the confinement and its depth, as these scales enter Eq. (4) by the power of 3 and 4, respectively. Below, we will refer to the inverse values of the timescaling estimates, i.e., $<\mathit{\tau}{>}^{\mathrm{1}}$, providing a measure of the glaciers' potential retreat speed relative to the reference PIG.
By far, the fastest retreat speed after a potential destabilization is detected for Mellor Glacier (MG; $<{\mathit{\tau}}_{\mathrm{MG}}{>}^{\mathrm{1}}=\mathrm{67}$). Here, the slowing influence of its relatively narrow confinement is strongly overcompensated by the combination of a deep trough (about 2000 m; together with Foundation Ice Stream (FIS), the deepest of the ensemble, see Table 3 and Fig. S8) and a steep retrograde bed slope (with about −55 m km^{−1}, the steepest of the entire ensemble). The fast retreat speed of Ninnis Glacier (NG; $<{\mathit{\tau}}_{\mathrm{NG}}{>}^{\mathrm{1}}=\mathrm{20}$) is primarily caused by its relatively steep retrograde slope. Similar retreat rates are obtained for the outlets of Thwaites Glacier (TG; $<{\mathit{\tau}}_{\mathrm{TG}}{>}^{\mathrm{1}}=\mathrm{19}$) and Cook Glacier (CG; $<{\mathit{\tau}}_{\mathrm{CG}}{>}^{\mathrm{1}}=\mathrm{16}$), which is for different reasons. While CG has the widest confinement among all investigated outlets (100 km), being almost twice as wide as TG’s confinement, its retrograde slope is only half as steep compared to the slope of TG. Both effects roughly balance out, and since the scaling ratios for the ice softness and the trough depth are very similar, respectively, this leads to the similar timescale. The reference glacier of our analysis, PIG, is found to have a much slower retreat speed. The effect of having the second softest ice of all examined outlets is outdone by the glacier’s comparatively gentle slope and narrow trough, respectively. Slower retreat rates are only found for Recovery Glacier (RG; $<{\mathit{\tau}}_{\mathrm{RG}}{>}^{\mathrm{1}}=\mathrm{0.8}$; slightly narrower, flatter bed trough) and FIS ($<{\mathit{\tau}}_{\mathrm{FIS}}{>}^{\mathrm{1}}=\mathrm{0.1}$; much flatter slope). Smith Glacier and Denman Glacier show a pronounced bimodal timescale distribution (highlighted in Fig. 8), which is an imprint of their bimodal distributions: for both glaciers, the inferred depth scales D span a relatively large range and cluster around the two endings of it, revealing two distinct depth regimes at the grounding line (Figs. S11 and S14). For some glaciers, the likely ranges of τ are very large, which is a combined effect of (1) our method of extracting the characteristic scales, capturing the full spatial variability inside the confinement and (2) the high degree of nonlinearity of the scaling relation, amplifying the uncertainty in the obtained scaling ratios. Note that we do not account for the uncertainty of the underlying topographic data itself. As the density of measurements is relatively high over many of the Antarctic outlet glaciers, the vertical accuracy over the regions considered in this study is typically on the order of ∼10 m, lying at the very low end of the overall Antarctic uncertainty range (Morlighem et al., 2020).
Despite having the fastest timescale, the sealevel contribution of a potential retreat of MG could be comparatively small (Fig. 7). For instance, the volume of marine ice in the drainage basin of Slessor Glacier (SLG)/RG and CG/NG is more than 2.5 times as large compared to the basin of MG. The consequences of a destabilization could thus be more serious for CG, NG and Support Force Glacier (SFG), as these tributaries show both a high retreat speed and a large sealevel potential.
Exploiting the similitude principle, we derived a timescaling relation for the case of fast, shallow and laterally confined outletglacier flow with negligible basal friction (Eq. 4). According to this relation, the timescale of an outlet is proportional to the inverse of its ice softness and highly nonlinearly related to the width, length and depth scales of its confinement. That is, deeper and wider outlets with a steeper retrograde slope have a faster timescale than shallower, more confined outlets with a flatter slope (Fig. 1). We showed that the derived scaling relation is able to predict the scaling of (1) observed flow velocities of Greenland outlet glaciers (Sect. 3.1; Fig. 2) and (2) Antarctictype MISI retreat speed as simulated in idealized experiments (Sect. 3.2; Fig. 6). Motivated by these simulations, we applied the scaling relation to the observed topography of 13 MISIprone Antarctic outlets to compare their time scaling of potential unstable retreat (Sect. 4; Fig. 7). Within this set of outlets, Mellor, Ninnis, Thwaites and Cook glaciers have the fastest timescales – about 16 to 67 times faster than the reference Pine Island Glacier (Fig. 8; Table 3) – due to their particularly wide and/or deep confinements and/or steep slopes (Figs. S8 and S10–S24). The underlying approach and its application have several limitations which restrict the results with respect to the following:

Similitude requirement. The application of Eq. (4) to observed glaciers is limited by the requirement of geometric similarity. That is to say that the approach can only be used to compare existing outlet glaciers if they can be scaled onto each other as shown in the Supplement figures (Figs. S10–S24a). Effects of differences in iceshelf buttressing or smaller topographic differences such as ice rises and details in the bed topography can alter the timescale of different glaciers in a way that cannot be captured by the approach because it weakens the similitude assumption. Perfect geometric similarity between two systems can only be achieved in a synthetic environment, i.e., in idealized experiments such as the ones carried out in this study to test the scaling relation (Sect. 3.2). The results from our numerical simulations suggest that the scaling relation is valid for strictly similar geometries (Fig. 6). The realworld glaciers considered here show geometric similarity in a broader sense.

Region of applicability. Our analysis is restricted to regions in which basal shear stresses are small compared to the iceinternal stresses (including the buttressinginducing lateral shear stresses). This assumption generally holds for coastal regions of most Antarctic outlet glaciers where basal shear stresses are found to be close to zero (Joughin et al., 2004; Morlighem et al., 2013; Sergienko and Hindmarsh, 2013; Sergienko et al., 2014) – with the possibility of the sporadic occurrence of confined and comparatively small riblike patterns of high basal resistance in some areas (Sergienko and Hindmarsh, 2013; Sergienko et al., 2014). The region of low basal friction is likely to follow the grounding line in case of retreat, as the lubrication of the bed is largely caused by the abundance of liquid water at the base of the ice sheet (Schoof, 2005; Martin et al., 2011; Leguy et al., 2014).

Time period of applicability. In principle, the applicability of our analysis is also limited by possible changes in surface mass balance. These are, however, likely slower than the icedynamic discharge. Estimates of the past decades show a background snowfall in coastal regions of at most 1 m yr^{−1} (Wessem et al., 2014; Agosta et al., 2019). Changes in this quantity can be estimated by a 5–7 % increase for every degree of warming in the atmosphere (Frieler et al., 2015), and are thus on the order of 0.05–0.07 m yr^{−1} K^{−1}. Differences in the surface mass balance between individual outlets are modeled to be on the order of 0.1 m yr^{−1} (Wessem et al., 2014; Agosta et al., 2019). In contrast, the highest observed dynamic thinning rates of the destabilized outlets in West Antarctica's Amundsen Sea Sector typically reach up to 5–10 m yr^{−1} (Pritchard et al., 2009; Konrad et al., 2018; Shepherd et al., 2019), making changes in the surface mass balance a secondorder effect with respect to the retreat timescale. Observed average groundingline retreat rates of these outlets reach up to 1 km yr^{−1} (Rignot et al., 2014; Konrad et al., 2018), from which we estimate that our results may be valid for at least 2 decades after a potential destabilization. After that, the analysis might have to be updated with the altered topography. A comparison of our numerical scaling results (insets of Fig. 6) to an additional set of simulations of altered surface accumulation yields that a scaling of the surface mass balance has a relatively weak influence on the time scaling of groundingline retreat (Fig. S6).

Relation between timescales of ice flow and retreat. The theoretical scaling relation is derived for the time scaling of the ice flow, e.g., the scaling of the flow speed of similar outlet glacier as examined in Sect. 3.1. We use the same equation to compare the time scaling of grounding line retreat, which requires that the timescales of ice flow and retreat, T_{flow} and T_{retr}, are proportional to each other. This proportionality is given under the assumption that changes in the iceflow speed due to retreat are proportional to the iceflow speed in equilibrium (as shown in Appendix A2). Although confirmed in the idealized simulations (Sect. 3.2), this assumption generally puts further constraints on the applicability of Eq. (4) to the potential retreat timescales of realworld glaciers.
Within these very strong constraints, the scaling relation can capture dependencies of highly complex stress and flow fields. It resolves the buttressing effect of the grounded portion of an outlet explicitly, by extracting its dependence on the geometry of the outlets' grounded lateral confinement. In our analysis, we infer the associated characteristic scales upstream of the grounding line focusing on the retrograde bed section as the crucial region that determines the timescale of a potential instability, and thus we do not take into account the iceshelf geometry. In our simulations, the ice shelves and the upstream confinement share the same width scale (Figs. 4 and S3). The outlets feeding the large Filchner–Ronne Ice Shelf deviate from this assumption. They might respond on a slower timescale than inferred here, due to the stabilizing influence of the large ice rises and other smaller pinning points in the central region of the ice shelf and the convergent flow of at least 10 ice streams that enter the ice shelf, effects that cannot be accounted for by our approach.
For the majority of outlets analyzed here, a comparison of the bed topography and the iceflow pattern suggests a strong control of the bed topography on the width of the outlet (e.g., Pine Island Glacier, Denman Glacier, Mellor Glacier or Slessor Glacier; Figs. S10, S14, S16, S17; visualized schematically in Fig. 1). That is, the lateral margins of these glaciers are clearly defined by the lateral ridges of the underlying bed trough. However, a few glaciers (e.g., Thwaites Glacier and Ninnis Glacier; Figs. S10 and S13) reveal a much weaker topographic control. We base the delineation of the outlets' lateral margins on the efolding scale of their velocity field. While other methods of deducing the glacierwidth scale may lead to different results in cases of weak topographic influence, our approach provides a welldefined width measure that can be applied consistently to all outlets, independent of their individual degree of topographic control.
Outlets discarded in this study include Pope Glacier and Kohler Glacier (eastern and western branch). Though both troughs of Kohler Glacier have a retrograde bed section upstream of the grounding line, their strongly undulating bed topography further upstream (Figs. S22 and S23) deviates substantially from the rest of the ensemble. Pope Glacier (Fig. S24) is excluded from our analysis due to its very short retrograde slope section (3 km). The main outlets feeding Ross Ice Shelf, e.g., Bindschadler Ice Stream and MacAyeal Ice Stream, and also Totten Glacier in East Antarctica, are found to have a generally, i.e., fluxweighted, flat or prograde bed topography upstream of the grounding line.
The timescales obtained here differ from the ones estimated in a previous study in which we also applied the similitude principle to Antarctic outlets (Levermann and Feldmann, 2019). This is due to the fact that the two studies make different assumptions in the dimensional analysis of the SSA leading to different scaling relations which yield different results. For instance, the onedimensional analysis in Levermann and Feldmann (2019) does not account for the acrossflow direction (y direction here), and thus neglects the effect of horizontal shearing. In the present study, horizontal shear stresses turn out to represent the most important stresses compared to all other iceinternal stresses, originating from the assumption of a generally small widthtolength ratio of Antarctica's outlet glaciers. In fact, the timescale of TG estimated here is about 10 times faster than in the previous study, since the present study takes into account the much larger width of TG (and thus the smaller buttressing effect of its grounded portion) compared to the reference PIG. In contrast, in the present study, basal stresses are neglected, whereas in Levermann and Feldmann (2019), there is a separate scaling law that involves basal friction parameter. Furthermore, the present study uses a more sophisticated method of extracting the outlets' characteristic scales, i.e., they are extracted for each glacier along a set of streamlines spanning its confinement compared to a single linear central profile used in Levermann and Feldmann (2019). The method of extraction also has an influence on which glaciers are labeled as MISI prone, such that the ensembles of outlets for which we obtain timescaling values presented here and in our previous work are not the same. Such qualitative differences between the two studies naturally lead to different timescale estimates, which should be interpreted in the light of the underlying assumptions.
The scaling approach presented here provides a useful tool that in particular allows to compare (MISI) timescales between similar glaciers within the limits of the underlying assumptions and constraints – without the need for computationally expensive numerical simulations. Further studies are required in particular to explore the role of iceshelf buttressing on the time scaling, allowing for a more comprehensive scaling analysis.
A1 Timescales of ice flow
We apply the concept of similitude (Buckingham, 1914; Rayleigh, 1915; Szücs, 1980) to the dynamics of ice sheets based on the shallowshelf approximation (SSA; Morland, 1987; MacAyeal, 1989) of the Stokes stress balance (e.g., Le Meur et al., 2004). Neglecting the terms of vertical shearing in the stress balance and accounting for the small thicknesstolength ratio of ice sheets, the SSA represents the relevant dynamics of floating ice shelves and grounded ice streams, characterized by fast, pluglike flow. In particular, we assume isothermal ice, a condition that has been used to analyze icesheet dynamics in a number of previous studies (e.g., Dupont and Alley, 2005; Goldberg et al., 2009; Gudmundsson et al., 2012; Pattyn et al., 2013; Feldmann and Levermann, 2016; Seroussi and Morlighem, 2018; Cornford et al., 2020).
For an ice stream with horizontal velocity components v_{x} and v_{y}, ice thickness H, ice surface elevation h, and spatially uniform (isothermal) ice softness A, the SSA has the following two horizontal components:
where ${\dot{\mathit{\u03f5}}}_{\mathrm{e}}$ is the effective strain rate, τ_{b} is the basal shear stress, ρ is the ice density, g the gravitational acceleration and n denotes Glen's flowlaw exponent (Glen, 1955). The SSA implies that the horizontal velocities do not vary in vertical (z) direction. Equation (A1) describes the balance between (1) iceinternal stresses due to shearing, extension and compression (in large brackets), (2) the basal shear stress ${\mathit{\tau}}_{\mathrm{b}}=({\mathit{\tau}}_{\mathrm{b},x},{\mathit{\tau}}_{\mathrm{b},y})$ between ice and bed and (3) the driving stress due to the gradient in the icesurface elevation. Here we consider the case in which the driving stress is mainly balanced by the iceinternal stresses, and we will neglect the basal shear stress term in the following.
The effective strain rate ${\dot{\mathit{\u03f5}}}_{\mathrm{e}}$ in Eq. (A1) can be written as
Introducing the geometric scales S_{x} and S_{y} as the two characteristic horizontal dimensions of the system and a timescale T, we bring the components of the effective strain rate (Eq. A2) into dimensionless form (Buckingham, 1914) via ${x}^{*}=\frac{x}{{S}_{x}}$, ${y}^{*}=\frac{y}{{S}_{y}}$, ${v}_{x}^{*}=\frac{{v}_{x}T}{{S}_{x}}$ and ${v}_{y}^{*}=\frac{{v}_{y}T}{{S}_{y}}$, which yields
with
Defining the dimensionless variables ${H}^{*}=\frac{H}{D}$ and ${h}^{*}=\frac{h}{D}$ with the vertical geometric scale D, we now nondimensionalize the x component of the SSA (Eq. A1), writing
with
We specify that the main flow direction, i.e., along the topographic confinement, is in x direction and that the transversal, acrosstrough, flow component is in y direction. In this case, the two geometric horizontal scales can be associated with a length scale L=S_{x} and a width scale W=S_{y} of the confinement. Observations indicate that many of the outlet glaciers of Greenland and Antarctica stream through bed confinements with a characteristic length that is several times (in places 1 order of magnitude) larger than the confinement width. Consequently, we assume L≫W and thus R≪1.
Accounting for the leading order term in the effective strain rate (Eq. A3), which is on the order of $\mathcal{O}(\mathrm{1}/R)$ (using Landau notation, e.g., Hardy and Wright, 1979), the effective strain rate can be approximated by
Analogous, in the SSA (Eq. A5), only the $\mathrm{1}/{R}^{\mathrm{2}}$ term remains on the lefthand side of the equation. Using Eq. (A7), we can eventually express the approximated SSA as
with
From Eq. (A8), we see that the stress balance remains exactly the same as long as the dimensionless coefficient λ remains the same. In other words, the firstorder icesheet dynamics of our problem are expected to be similar under a transformation that leaves this coefficient unchanged. Applying a scaling to the system (denoted by a dash) and requiring similarity (${\mathit{\lambda}}^{\prime}=\mathit{\lambda}$) yields
Introducing the scaling ratios $\mathit{\alpha}=\frac{{A}^{\prime}}{A}$, $\mathit{\beta}=\frac{{D}^{\prime}}{D}$, $\mathit{\omega}=\frac{{W}^{\prime}}{W}$, $\mathit{\delta}=\frac{{L}^{\prime}}{L}$ and $\mathit{\tau}=\frac{{T}^{\prime}}{T}$ for the geometric dimensions, the ice softness and the timescale of the system, respectively, we arrive at
Since n is always positive, the timescale of the system will increase (τ>1) for a reduction of the icesoftness scale (α<1), a shortening of the vertical geometry (β<1), an elongation of the geometry in mainflow direction (δ>1) or a narrowing of the acrossmainflow geometry (ω<1). Defining $\mathit{\gamma}={R}^{\prime}/R=\mathit{\omega}/\mathit{\delta}$ and making the common choice of n=3 (Glen, 1955; Duval, 2013), we eventually obtain
A2 Relation between timescales of ice flow and groundingline retreat
Note that Eq. (A12) refers to the timescale of glacier flow, T_{flow}, which is associated with the velocity of ice flow in equilibrium (see Appendix B). In order to find a relation of this timescale to the timescale of groundingline retreat, T_{retr}, in the following, we analyze the time evolution of a change in the ice thickness close to the grounding line of a perturbed outlet. We define this change in ice thickness as
where H_{eq} is the ice thickness in equilibrium and H is the ice thickness in the perturbed state. The time evolution of δH then reads
since $\frac{\partial {H}_{\mathrm{eq}}}{\partial t}=\mathrm{0}$ per definition. We now introduce the ice thickness equation (e.g., Whillans, 1977), which relates changes in the ice thickness to the divergence of the ice flux Q and the surface accumulation a
Accounting for the fact that the fast retreat of Antarctic outlets such Pine Island Glacier and Thwaites Glacier is due to their icedynamic imbalance (Pritchard et al., 2009; Konrad et al., 2018; Shepherd et al., 2019), i.e., due to changes in the ice flux, we can neglect changes in a, and write
where we used Q=Hu with the iceflow velocity u in equilibrium. Note that we used $\mathit{\delta}\mathrm{\nabla}\cdot \left(Hu\right)=\mathrm{\nabla}\cdot \mathit{\delta}\left(Hu\right)$ due to the linearity of the nabla operator ∇. We now write Eq. (A16) in full, which yields
If we make the plausible assumption that the equilibrium velocity profile increases in downstream direction, i.e., $\mathrm{\nabla}\cdot u>\mathrm{0}$, then the first term on the righthand side of Eq. (A17) has a stabilizing character. That is, the change δH and its temporal evolution, $\frac{\partial \mathit{\delta}H}{\partial t}$, have opposite signs such that the change is dampened over time. Although this is an important term in the general case, we neglect it here, since it cannot be the decisive term in the case of selfsustained groundingline retreat, which we analyze here. If we consider the spatial derivatives of δH and δu to be of second order, then Eq. (A17), in combination with Eq. (A14), reduces to
Writing down the dimensions of Eq. (A18), its lefthand side reads
with icethickness scale D and timescale of retreat T_{retr}, which is associated with the thinning rate. The righthand side of Eq. (A18) yields
with length scale L. To find an expression for [δu], we make the assumption that the characteristic scale of a change in the iceflow velocity is proportional to the iceflow velocity itself, i.e.,
where T_{flow} is the timescale of the equilibrium iceflow velocity. The dimensional analysis of Eq. (A18) thus yields
Considering a scaled (dashed) and an unscaled system with ${\mathit{\tau}}_{\mathrm{retr}}=\frac{{T}_{\mathrm{retr}}^{\prime}}{{T}_{\mathrm{retr}}}$ and ${\mathit{\tau}}_{\mathrm{flow}}=\frac{{T}_{\mathrm{flow}}^{\prime}}{{T}_{\mathrm{flow}}}$, it follows
i.e., the time scaling of groundingline retreat is equal to the time scaling of ice flow.
The observational data used in our analysis are obtained from Beckmann et al. (2019) who inferred alongflow profiles of the bed elevation, ice thickness, ice velocity and glacier width for a selection of Greenland outlet glaciers, based on fluxweighted averages using geometric data from Morlighem et al. (2014) and velocity data from Rignot and Mouginot (2012). The 12 glaciers cover a variety of sizes, locations, discharge rates and climatic conditions across the Greenland Ice Sheet (Beckmann et al., 2019). At the same time, most of these outlets are qualitatively similar to each other, as their flow is strongly confined by a narrow lateral bed geometry and they are grounded on bedrock that is generally sloping up in landward direction (Figs. 2 and S1). Two glaciers deviating from these conditions, i.e., Daugaard–Jensen Glacier and Rink Glacier, are excluded from our analysis. Daugaard–Jensen Glacier is the only outlet in the ensemble that is grounded on a landward downsloping (instead of upsloping) bed topography. It has been shown that these two types of bed slope can have fundamentally different effects on the iceflow dynamics (Schoof, 2007). Rink Glacier exhibits a pronounced inversion of the iceflow acceleration near the terminus, which might be due to a local specific topographic feature, making the outlet's icespeed profile qualitatively very different compared to the rest of the outlets (except for Upernavik Isstrom N). Since the inversion for the case of Upernavik Isstrom N is much less pronounced, shorter in extent and closer to the calving front, we decided to include this glacier in our analysis. In fact, these choices involve some degree of subjectivity and a different approach of ensuring qualitative similarity between the glaciers might lead to a different set of valid glaciers to be analyzed.
To extract a characteristic scale for the confinement width W of each of the individual glaciers, we average the width profile over the first 60 km upstream of the glacier terminus. The same is done for the velocity scale V, though here leaving out the first 20 km to exclude the possible influence of ice–ocean interactions, i.e., melting and calving events, and their effects near the glacier terminus. The analysis of Gade Glacier, with its relatively short catchment, is limited to the first 40 km upstream of its terminus.
With the flowspeed scale $V=L/T$, the speed scaling ratio $\mathit{\nu}={V}^{\prime}/V$ can be written as
Assuming that both the time scaling and the velocity scaling of the considered glaciers are dominated by the width of their narrow topographic confinements ($\mathit{\alpha}=\mathit{\beta}=\mathit{\delta}=\mathrm{1}$), the analytically derived scaling relation (Eq. 4) reduces to
and consequently
which provides a good approximation to the scaling behavior deduced from the observational data (Fig. 2). Note that the results are affected by the concrete choice of averaging ranges along the glaciers' velocity and width profiles from which the scales of V and W are estimated (Fig. S2). For instance, taking the full range of the velocity profile to calculate the velocity average instead of leaving out the oceanward third of the profile leads to larger deviations from the theoretically predicted scaling curve (Fig. S2, compare panels A/B to C/D).
We apply PISM to a modified version of the MISMIP+ setup to simulate the unstable retreat of an inherently buttressed Antarctictype outlet glacier (Figs. 4 and S3). In the default MISMIP+ setting, the model spinup yields a steadystate outlet glacier with a stable groundingline position on the retrograde slope section of the bed (AsayDavis et al., 2016; Cornford et al., 2020). A perturbation of this equilibrium leads to a reversible retreat of the outlet. That is, the triggered transient retreat stops and the grounding line migrates back toward its original equilibrium location once the perturbation ceases, even for strong perturbation magnitudes of basal iceshelf melting or a reduction of the basal resistance (Gudmundsson et al., 2012; Cornford et al., 2020). In order to be able to extract timescales of irreversible, i.e., MISItype retreat, we widened and deepened the bed trough of the MISMIP+ setup. For this purpose, we doubled the extent of the computational domain in y direction (320 km instead of 160 km) and also extended it in x direction (800 km instead of 700 km). In MISMIP+, the bed topography results from the superposition of the two components B_{x}(x) and B_{y}(y), representing the two horizontal dimensions. In our experiments, the y component, B_{y}(y), is the same as in the MISMIP+ setup (Cornford et al., 2020, Eq. 1), but we generally apply larger values for the channel width w_{c} throughout our experiments (see Table 2 for the range of used values). The x component, B_{x}(x), is qualitatively very similar to the one used in MISMIP+, but it is generated using a piecewise cubic spline interpolation instead of prescribing a polynomial. This allows for a convenient prescription of the position and elevation of the endpoints of the retrograde bed section: at the nodes x_{0}=0 (location of the inland summit), x_{1}=200 and x_{2}=500 km (landward and oceanward ends of the retrograde bed section) and x_{3}=700 km, we prescribe the bed elevation (${B}_{x}\left({x}_{\mathrm{0}}\right)=\mathrm{400}\phantom{\rule{0.125em}{0ex}}\mathrm{m}$, ${B}_{x}\left({x}_{\mathrm{1}}\right)=\mathrm{750}\phantom{\rule{0.125em}{0ex}}\mathrm{m}$, ${B}_{x}\left({x}_{\mathrm{2}}\right)=\mathrm{600}\phantom{\rule{0.125em}{0ex}}\mathrm{m}$ and ${B}_{x}\left({x}_{\mathrm{3}}\right)=\mathrm{800}\phantom{\rule{0.125em}{0ex}}\mathrm{m}$) and set its first derivative to zero (since the nodes are locations of local extrema). Beyond x_{3}, the bed elevation is set to ${B}_{x}(x>{x}_{\mathrm{3}})=\mathrm{800}\phantom{\rule{0.125em}{0ex}}\mathrm{m}$. The shape of B_{x}(x) resulting from the interpolation is shown in Fig. 4b.
For a sufficiently wide trough, the outlet glacier that evolves after model initialization retreats down the entire retrograde slope, finding a stable equilibrium only on the landward upsloping bed section close to the ice divide (reddish lines in Fig. 5). Using this experimental design, all of our conducted simulations share the same initial conditions (except for the varied parameter), and thus have a high degree of similarity, which facilitates the application of the scaling to our numerical results.
The experiments are initiated from a block of ice of 2000 m thickness and are run into equilibrium for several 10 000 model years. Ice surface accumulation due to snowfall takes place at a spatially uniform rate (a=0.15 m yr^{−1}) and is constant in time. Basal and surface melting of the ice are neglected. The temperature within the ice body, T_{ice}, is chosen to be uniform (isothermal) and is translated into an icesoftness value A(T_{ice}) by an Arrhenius law (Glen, 1955) of the form
with constants A_{0}, Q and R_{r} (see Table 2).
For the grounded portion of the ice sheet, the basal shear stress in Eq. (A1), ${\mathit{\tau}}_{\mathrm{b}}=({\mathit{\tau}}_{\mathrm{b},x},{\mathit{\tau}}_{\mathrm{b},y})$, is given by a Weertmantype sliding law (e.g., Pattyn et al., 2013),
with constant friction coefficient C and basal sliding exponent m. The parameter values used in the simulations are given in Table 2. Like in several earlier studies, which involved simulations of icesheet dynamics in idealized settings (e.g., Feldmann and Levermann, 2015, 2016; Feldmann et al., 2022), we interpolate the basal friction at the grounding line according to a subgrid interpolation of the grounding line position (Gladstone et al., 2010; Feldmann et al., 2014). This approach attenuates the discontinuity in the basal stresses when going from the last grounded (nonzero basal friction) into the first floating cell (zero basal friction). Note that more complex sliding laws, such as a Coulomblimited law, intrinsically ensure that the basal stresses approach zero within a transition zone at the grounding line. The application of the Coulomblimited law from Tsai et al. (2015) in PISM in the course of the MISMIP+ intercomparison exercise (Cornford et al., 2020) revealed a narrow transition zone of 2 km (equal to the horizontal resolution used here). Thus, in both model realizations (Coulomblimited law vs. Weertman law with friction interpolation), the drop in basal stress takes place over the same short distance of 2 km upstream of the grounding line. Both model realizations lead to qualitatively the same response to the applied perturbations in the MISMIP+ experiments, with moderate deviations in the initial groundingline position and the response magnitude (Cornford et al., 2020).
The simulations are carried out on a regular horizontal grid of 2 km resolution, which assures adequate accuracy in modeling the rapid ice dynamics, as demonstrated by a convergence study (Fig. S7). In PISM, the grounding lines are diagnosed via the flotation criterion and thus evolve freely. Grounding line movement has been evaluated in the model intercomparison projects MISMIP3d (Pattyn et al., 2013; Feldmann et al., 2014) and MISMIP+ (Cornford et al., 2020).
D1 Method
To apply the above derived theory to Antarctica's outlet glaciers, their characteristic scales involved in Eq. (4) have to be extracted from observations. To this end, for each outlet first we define a groundingline section that spans the glacier's main trunk, i.e., the part of the grounding line across which the vast majority of the glacier's total ice discharge takes place. Then, based on the velocity field of the ice, we infer the glacier's streamlines that meet this groundingline section. The sampling resolution between and along the streamlines is 1 km. The fluxweighted average of the outlet's bed profile over all streamlines provides a general start point (p_{1}) and end point (p_{2}) of its retrograde bed section directly upstream of the grounding line (exemplary indicated in Fig. S10). For each streamline, the confinement depth D is then measured as the depth of the bed below zero at location p_{1}. The length scale L of the retrograde bed section is obtained from the bedslope magnitude S via $L=\mathrm{1}/S$ for each of the streamlines that have a negative mean slope between p_{1} and p_{2}. The confinement width W is measured in perpendicular direction to a central streamline of the outlet (indicated in black in Figs. S10–S24) as the distance at which the iceflow speed has declined to the eth fraction of the respective centerline value. This method of detecting the glacier's lateral shear margins facilitates the localization of its lateral boundary, especially if the outlet is rather weakly confined by the bed topography (e.g., Thwaites Glacier, Fig. S10). For outlets that are instead sharply confined by the bed (e.g., Pine Island Glacier), the velocitydeduced shear margins agree well with the lateral topographic confinement (Fig. S10). The values for W are obtained along the first 20 km upstream of p_{1}. This captures the nonlocal influence of the confinement width W (due to the nonlocal nature of the stress balance, Eq. 1) over the same distance upstream of p_{1} for each outlet, independently of the length of its retrograde slope section. Over the same range, the ice softness A is averaged for each streamline. The above approach yields a set of measurements for each characteristic scale, spanning the observational range of this scale. The number of measurements for the scales D,L and A is given by the number of the glacierspecific (retrograde) streamlines, whereas for W it is given by the 1 km sampling along the central streamline.
To ensure similarity between the measured outlets and comparability of their timescales via Eq. (4), we apply the following three criteria:

The outlets have to show icestream characteristics. The maximum flow speeds of the grounded part of the glacier have to be at least 100 m yr^{−1} (most of them are streaming at >500 m yr^{−1}).

MISI susceptibility. Within the first 20 km upstream of the grounding line, the calculated fluxaveraged bed profile has to be below sea level and it has to continuously slope down in landward direction over a distance of at least 5 km (the vast majority of the analyzed glaciers has an initial retrograde bed section of 10 km or more).

General geometric similarity. While each glacier's profile has of course its own characteristics, we require that fluctuations in the bedslope direction have to be modest, i.e., we discard glaciers with more than one pronounced peak in the bed profile within the first 20 km (which excludes the strongly undulating Kohler Glacier tributaries from our analysis).
Applying these criteria, we obtain an ensemble of 13 outlet glaciers around Antarctica that are suitable to be involved in our analysis. Note that since the choice of the similarity criteria is not unique a different set of criteria than the one used here may lead to a different ensemble of outlets.
We choose Pine Island Glacier as the reference system for calculating the scaling ratios α, β, δ and ω that enter the timescaling relation. That is, the scales of all outlets analyzed in this study are expressed relatively to the scales of this reference. Pine Island Glacier is one of the most prominent Antarctic outlets and it is relatively well observed. Due to the nature of the conducted scaling analysis, the scales calculated here could also be expressed relative to any other of the examined outlets without changing the results. Each scaling ratio of a specific glacier is considered as a distribution of results which accounts for the uncertainty of the two involved scales. That is, for each scaling ratio, all measurements of the scale of the specific glacier and the reference are combined with each other. The distribution for the timescaling ratio τ is then calculated according to Eq. (4), based on the likely range (range between the 17th and 83rd percentiles) of the distributions of α, β, δ and ω, by conducting Monte Carlo simulations with a sample size of 1000 (Fig. 8). This likely range helps to confine the resulting timescaling distributions: outliers in the measurements of the characteristic scales can generate strong uncertainties in the calculated scaling ratios which, in turn, are amplified through the nonlinearity of Eq. (4) in combination with the Monte Carlo method. Due to the large uncertainty in the observed alongflow bed shape (Morlighem et al., 2020) and the associated relatively large spread in the length scale, the δ scaling ratio is exempt from the Monte Carlo method, and the median of the distribution is used in Eq. (4) instead.
D2 Underlying datasets
The characteristic scales of the outlet glaciers are all obtained from datasets that represent presentday conditions of the Antarctic Ice Sheet. The bed topography stems from the Bedmachine Antarctica dataset (Morlighem et al., 2020), from which we infer the outletspecific geometric scales D, L and W. The regularly updated, continentwide compilation involves data from various sources, including satellite, airborne radar, oversnow radar and seismicsounding measurements. The novel dataset applies a mass conservation method to overcome limitations of the earlier BEDMAP2 compilation (Fretwell et al., 2013), revealing previously unknown bed features of Antarctica's outlets. with substantial implications regarding their susceptibility to MISI.
Due to the absence of observed Antarctic icesoftness data, we use results from a presentday, continentalscale simulation of the Antarctic Ice Sheet, conducted with PISM in the course of the recent ISMIP 6 Antarctica intercomparison project (Seroussi et al., 2020). The thermomechanically coupled model is run into equilibrium under presentday conditions at a horizontal resolution of 8 km and a vertical resolution that ranges from 13 m at the ice base to 100 m at the top of the numerical domain. As part of solving the stress balance of the ice, PISM computes a temperaturedependent ice softness field, from which the vertically averaged form serves as the underlying dataset for our inference of the icesoftness scale A (Fig. S9).
The ice flow field used to track the streamlines of the individual outlet glaciers is taken from Mouginot et al. (2019a). It is the most recent and precise dataset of the presentday Antarctic ice surface velocity, derived by applying a combination of interferometric phasemapping and speckletracking techniques to satellite measurements from the last 25 years.
Antarctic bed topography and ice thickness are taken from the Bedmachine Antarctica compilation (MEaSUREs BedMachine Antarctica, Version 1) and can be downloaded from https://doi.org/10.5067/C2GFER6PTOS4 (Morlighem, 2019). Antarctic surface velocities (MEaSUREs PhaseBased Antarctica Ice Velocity Map, Version 1) are available at https://doi.org/10.5067/PZ3NJ5RXRH10 (Mouginot et al., 2019b). The simulation data generated and analyzed in this study are available from https://doi.org/10.5281/zenodo.7395830 (Feldmann, 2022). PISM is freely available as opensource code from https://www.pism.io (PISM authors, 2021). The code version used in this study can be obtained from https://doi.org/10.5281/zenodo.7395020 (Feldmann and PISM authors, 2022).
The supplement related to this article is available online at: https://doi.org/10.5194/tc173272023supplement.
Both authors designed the study, developed the theory and wrote the paper. JF conducted the model simulations and prepared the figures.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We thank Camilla Schelpe and one anonymous referee for their valuable comments and suggestions which helped to improve the paper.
This work was supported by the Deutsche Forschungsgemeinschaft (DFG) in the framework of the priority programme “Antarctic Research with comparative investigations in Arctic ice areas” through grants LE 1448/81 and WI 4556/61. We acknowledge the European Regional Development Fund (ERDF), the German Federal Ministry of Education and Research, and the Land Brandenburg for supporting this project by providing resources on the highperformance computer system at the Potsdam Institute for Climate Impact Research. Development of PISM is supported by NSF grants PLR1644277 and PLR1914668 and NASA grants NNX17AG65G and 20CRYO20200052.
This paper was edited by Alexander Robinson and reviewed by Camilla Schelpe and one anonymous referee.
Agosta, C., Amory, C., Kittel, C., Orsi, A., Favier, V., Gallée, H., van den Broeke, M. R., Lenaerts, J. T. M., van Wessem, J. M., van de Berg, W. J., and Fettweis, X.: Estimation of the Antarctic surface mass balance using the regional climate model MAR (1979–2015) and identification of dominant processes, The Cryosphere, 13, 281–296, https://doi.org/10.5194/tc132812019, 2019. a, b
AsayDavis, X. S., Cornford, S. L., Durand, G., GaltonFenzi, B. K., Gladstone, R. M., Gudmundsson, G. H., Hattermann, T., Holland, D. M., Holland, D., Holland, P. R., Martin, D. F., Mathiot, P., Pattyn, F., and Seroussi, H.: Experimental design for three interrelated marine ice sheet and ocean model intercomparison projects: MISMIP v. 3 (MISMIP +), ISOMIP v. 2 (ISOMIP +) and MISOMIP v. 1 (MISOMIP1), Geosci. Model Dev., 9, 2471–2497, https://doi.org/10.5194/gmd924712016, 2016. a, b
Bamber, J. L., Oppenheimer, M., Kopp, R. E., Aspinall, W. P., and Cooke, R. M.: Ice Sheet Contributions to Future SeaLevel Rise from Structured Expert Judgment, P. Natl. Acad. Sci. USA, 116, 11195–11200, https://doi.org/10.1073/pnas.1817205116, 2019. a
Beckmann, J., Perrette, M., Beyer, S., Calov, R., Willeit, M., and Ganopolski, A.: Modeling the response of Greenland outlet glaciers to global warming using a coupled flow line–plume model, The Cryosphere, 13, 2281–2301, https://doi.org/10.5194/tc1322812019, 2019. a, b, c, d
Buckingham, E.: On Physically Similar Systems; Illustrations of the Use of Dimensional Equations, Phys. Rev., 4, 345–376, https://doi.org/10.1103/PhysRev.4.345, 1914. a, b
Bueler, E. and Brown, J.: Shallow Shelf Approximation as a “Sliding Law” in a Thermomechanically Coupled Ice Sheet Model, J. Geophys. Res., 114, F03008, https://doi.org/10.1029/2008JF001179, 2009. a, b
Catania, G. A., Stearns, L. A., Moon, T. A., Enderlin, E. M., and Jackson, R. H.: Future Evolution of Greenland's MarineTerminating Outlet Glaciers, J. Geophys. Res.Earth, 125, e2018JF004873, https://doi.org/10.1029/2018JF004873, 2020. a, b
Cornford, S. L., Seroussi, H., AsayDavis, X. S., Gudmundsson, G. H., Arthern, R., Borstad, C., Christmann, J., Dias dos Santos, T., Feldmann, J., Goldberg, D., Hoffman, M. J., Humbert, A., Kleiner, T., Leguy, G., Lipscomb, W. H., Merino, N., Durand, G., Morlighem, M., Pollard, D., Rückamp, M., Williams, C. R., and Yu, H.: Results of the third Marine Ice Sheet Model Intercomparison Project (MISMIP+), The Cryosphere, 14, 2283–2301, https://doi.org/10.5194/tc1422832020, 2020. a, b, c, d, e, f, g, h, i, j, k
DeConto, R. M., Pollard, D., Alley, R. B., Velicogna, I., Gasson, E., Gomez, N., Sadai, S., Condron, A., Gilford, D. M., Ashe, E. L., Kopp, R. E., Li, D., and Dutton, A.: The Paris Climate Agreement and Future SeaLevel Rise from Antarctica, Nature, 593, 83–89, https://doi.org/10.1038/s41586021034270, 2021. a
Diener, T., Sasgen, I., Agosta, C., Fürst, J. J., Braun, M. H., Konrad, H., and Fettweis, X.: Acceleration of Dynamic Ice Loss in Antarctica From Satellite Gravimetry, Front. Earth Sci., 9, 741789, https://doi.org/10.3389/feart.2021.741789, 2021. a, b, c
Dupont, T. K. and Alley, R. B.: Assessment of the Importance of IceShelf Buttressing to IceSheet Flow, Geophys. Res. Lett., 32, L04503, https://doi.org/10.1029/2004GL022024, 2005. a, b, c, d
Duval, P.: Creep Behavior of Ice in Polar Ice Sheets, in: The Science of Solar System Ices, edited by: Gudipati, M. S. and CastilloRogez, J., Astrophysics and Space Science Library, Springer, New York, NY, 227–251, https://doi.org/10.1007/9781461430766_8, 2013. a, b
Edwards, T. L., Nowicki, S., Marzeion, B., Hock, R., Goelzer, H., Seroussi, H., Jourdain, N. C., Slater, D. A., Turner, F. E., Smith, C. J., McKenna, C. M., Simon, E., AbeOuchi, A., Gregory, J. M., Larour, E., Lipscomb, W. H., Payne, A. J., Shepherd, A., Agosta, C., Alexander, P., Albrecht, T., Anderson, B., AsayDavis, X., Aschwanden, A., Barthel, A., Bliss, A., Calov, R., Chambers, C., Champollion, N., Choi, Y., Cullather, R., Cuzzone, J., Dumas, C., Felikson, D., Fettweis, X., Fujita, K., GaltonFenzi, B. K., Gladstone, R., Golledge, N. R., Greve, R., Hattermann, T., Hoffman, M. J., Humbert, A., Huss, M., Huybrechts, P., Immerzeel, W., Kleiner, T., Kraaijenbrink, P., Le Clec'h, S., Lee, V., Leguy, G. R., Little, C. M., Lowry, D. P., Malles, J.H., Martin, D. F., Maussion, F., Morlighem, M., O'Neill, J. F., Nias, I., Pattyn, F., Pelle, T., Price, S. F., Quiquet, A., Radić, V., Reese, R., Rounce, D. R., Rückamp, M., Sakai, A., Shafer, C., Schlegel, N.J., Shannon, S., Smith, R. S., Straneo, F., Sun, S., Tarasov, L., Trusel, L. D., Van Breedam, J., van de Wal, R., van den Broeke, M., Winkelmann, R., Zekollari, H., Zhao, C., Zhang, T., and Zwinger, T.: Projected Land Ice Contributions to TwentyFirstCentury Sea Level Rise, Nature, 593, 74–82, https://doi.org/10.1038/s4158602103302y, 2021. a
Favier, L., Durand, G., Cornford, S. L., Gudmundsson, G. H., Gagliardini, O., GilletChaulet, F., Zwinger, T., Payne, A. J., and Le Brocq, A. M.: Retreat of Pine Island Glacier Controlled by Marine IceSheet Instability, Nat. Clim. Change, 4, 117–121, https://doi.org/10.1038/nclimate2094, 2014. a, b
Feldmann, J.: Idealized simulations of marine ice sheet instability, Zenodo [data set], https://doi.org/10.5281/zenodo.7395830, 2022. a
Feldmann, J. and Levermann, A.: Interaction of marine icesheet instabilities in two drainage basins: simple scaling of geometry and transition time, The Cryosphere, 9, 631–645, https://doi.org/10.5194/tc96312015, 2015. a
Feldmann, J. and Levermann, A.: Similitude of ice dynamics against scaling of geometry and physical parameters, The Cryosphere, 10, 1753–1769, https://doi.org/10.5194/tc1017532016, 2016. a, b, c
Feldmann, J. and PISM authors: PISM code for modeling instability timescales (v1.0663ga7de0cb71), Zenodo [code], https://doi.org/10.5281/zenodo.7395020, 2022. a
Feldmann, J., Albrecht, T., Khroulev, C., Pattyn, F., and Levermann, A.: ResolutionDependent Performance of Grounding Line Motion in a Shallow Model Compared with a FullStokes Model According to the MISMIP3d Intercomparison, J. Glaciol., 60, 353–360, https://doi.org/10.3189/2014JoG13J093, 2014. a, b
Feldmann, J., Reese, R., Winkelmann, R., and Levermann, A.: Shearmargin melting causes stronger transient ice discharge than icestream melting in idealized simulations, The Cryosphere, 16, 1927–1940, https://doi.org/10.5194/tc1619272022, 2022. a, b
Fretwell, P., Pritchard, H. D., Vaughan, D. G., Bamber, J. L., Barrand, N. E., Bell, R., Bianchi, C., Bingham, R. G., Blankenship, D. D., Casassa, G., Catania, G., Callens, D., Conway, H., Cook, A. J., Corr, H. F. J., Damaske, D., Damm, V., Ferraccioli, F., Forsberg, R., Fujita, S., Gim, Y., Gogineni, P., Griggs, J. A., Hindmarsh, R. C. A., Holmlund, P., Holt, J. W., Jacobel, R. W., Jenkins, A., Jokat, W., Jordan, T., King, E. C., Kohler, J., Krabill, W., RigerKusk, M., Langley, K. A., Leitchenkov, G., Leuschen, C., Luyendyk, B. P., Matsuoka, K., Mouginot, J., Nitsche, F. O., Nogi, Y., Nost, O. A., Popov, S. V., Rignot, E., Rippin, D. M., Rivera, A., Roberts, J., Ross, N., Siegert, M. J., Smith, A. M., Steinhage, D., Studinger, M., Sun, B., Tinto, B. K., Welch, B. C., Wilson, D., Young, D. A., Xiangbin, C., and Zirizzotti, A.: Bedmap2: improved ice bed, surface and thickness datasets for Antarctica, The Cryosphere, 7, 375–393, https://doi.org/10.5194/tc73752013, 2013. a
Frieler, K., Clark, P. U., He, F., Buizert, C., Reese, R., Ligtenberg, S. R. M., van den Broeke, M. R., Winkelmann, R., and Levermann, A.: Consistent Evidence of Increasing Antarctic Accumulation with Warming, Nat. Clim. Change, 5, 348–352, https://doi.org/10.1038/nclimate2574, 2015. a
Gladstone, R. M., Payne, A. J., and Cornford, S. L.: Parameterising the grounding line in flowline ice sheet models, The Cryosphere, 4, 605–619, https://doi.org/10.5194/tc46052010, 2010. a
Glen, J. W.: The Creep of Polycrystalline Ice, P. Roy. Soc. Lond. A Mat., 228, 519–538, https://doi.org/10.1098/rspa.1955.0066, 1955. a, b, c, d, e
Goldberg, D., Holland, D. M., and Schoof, C.: Grounding Line Movement and Ice Shelf Buttressing in Marine Ice Sheets, J. Geophys. Res., 114, F04026, https://doi.org/10.1029/2008JF001227, 2009. a
Gudmundsson, G. H., Krug, J., Durand, G., Favier, L., and Gagliardini, O.: The stability of grounding lines on retrograde slopes, The Cryosphere, 6, 1497–1505, https://doi.org/10.5194/tc614972012, 2012. a, b, c, d, e, f, g
Hardy, G. H. and Wright, E. M.: Some Notations, in: An Introduction to the Theory of Numbers, Oxford University Press, 7–8, ISBN 9780198531715, 1979. a
Haseloff, M. and Sergienko, O. V.: The Effect of Buttressing on Grounding Line Dynamics, J. Glaciol., 64, 417–431, https://doi.org/10.1017/jog.2018.30, 2018. a, b
Haseloff, M., Schoof, C., and Gagliardini, O.: A Boundary Layer Model for Ice Stream Margins, J. Fluid Mech., 781, 353–387, https://doi.org/10.1017/jfm.2015.503, 2015. a
Hughes, T.: Is the West Antarctic Ice Sheet Disintegrating?, J. Geophys. Res., 78, 7884–7910, https://doi.org/10.1029/JC078i033p07884, 1973. a
Joughin, I., MacAyeal, D. R., and Tulaczyk, S.: Basal Shear Stress of the Ross Ice Streams from Control Method Inversions, J. Geophys. Res.Sol. Ea., 109, B09405, https://doi.org/10.1029/2003JB002960, 2004. a, b, c
Joughin, I., Smith, B. E., and Medley, B.: Marine Ice Sheet Collapse Potentially Under Way for the Thwaites Glacier Basin, West Antarctica, Science, 344, 735–738, https://doi.org/10.1126/science.1249055, 2014. a, b
Joughin, I., Smith, B. E., and Howat, I. M.: A Complete Map of Greenland Ice Velocity Derived from Satellite Data Collected over 20 Years, J. Glaciol., 64, 1–11, https://doi.org/10.1017/jog.2017.73, 2018. a
Konrad, H., Shepherd, A., Gilbert, L., Hogg, A. E., McMillan, M., Muir, A., and Slater, T.: Net Retreat of Antarctic Glacier Grounding Lines, Nat. Geosci., 11, 258–262, https://doi.org/10.1038/s415610180082z, 2018. a, b, c
Kundu, P., Cohen, I., and Hu, H.: Fluid Mechanics, Edition, Academic Press, New Delhi, https://doi.org/10.1016/C20090634103, 2012. a
Le Meur, E., Gagliardini, O., Zwinger, T., and Ruokolainen, J.: Glacier Flow Modelling: A Comparison of the Shallow Ice Approximation and the FullStokes Solution, C. R. Phys., 5, 709–722, https://doi.org/10.1016/j.crhy.2004.10.001, 2004. a
Leguy, G.: Chapter 7.1.5, in: The Effect of a BasalFriction Parameterization on GroundingLine Dynamics in IceSheet Models (PhD Thesis), New Mexico Institute of Mining and Technology, 3723078, https://www.proquest.com/openview/e50f74521532de35c3c3e5d8141ee50f/1.pdf?pqorigsite=gscholar&cbl=18750 (last access: 22 January 2023), 2015. a, b
Leguy, G. R., AsayDavis, X. S., and Lipscomb, W. H.: Parameterization of basal friction near grounding lines in a onedimensional ice sheet model, The Cryosphere, 8, 1239–1259, https://doi.org/10.5194/tc812392014, 2014. a
Leguy, G. R., Lipscomb, W. H., and AsayDavis, X. S.: Marine ice sheet experiments with the Community Ice Sheet Model, The Cryosphere, 15, 3229–3253, https://doi.org/10.5194/tc1532292021, 2021. a
Levermann, A. and Feldmann, J.: Scaling of instability timescales of Antarctic outlet glaciers based on onedimensional similitude analysis, The Cryosphere, 13, 1621–1633, https://doi.org/10.5194/tc1316212019, 2019. a, b, c, d, e, f, g
Levermann, A., Winkelmann, R., Albrecht, T., Goelzer, H., Golledge, N. R., Greve, R., Huybrechts, P., Jordan, J., Leguy, G., Martin, D., Morlighem, M., Pattyn, F., Pollard, D., Quiquet, A., Rodehacke, C., Seroussi, H., Sutter, J., Zhang, T., Van Breedam, J., Calov, R., DeConto, R., Dumas, C., Garbe, J., Gudmundsson, G. H., Hoffman, M. J., Humbert, A., Kleiner, T., Lipscomb, W. H., Meinshausen, M., Ng, E., Nowicki, S. M. J., Perego, M., Price, S. F., Saito, F., Schlegel, N.J., Sun, S., and van de Wal, R. S. W.: Projecting Antarctica's contribution to future sea level rise from basal ice shelf melt using linear response functions of 16 ice sheet models (LARMIP2), Earth Syst. Dynam., 11, 35–76, https://doi.org/10.5194/esd11352020, 2020. a
Lhermitte, S., Sun, S., Shuman, C., Wouters, B., Pattyn, F., Wuite, J., Berthier, E., and Nagler, T.: Damage Accelerates Ice Shelf Instability and Mass Loss in Amundsen Sea Embayment, P. Natl. Acad. Sci. USA, 117, 24735–24741, https://doi.org/10.1073/pnas.1912890117, 2020. a
MacAyeal, D. R.: LargeScale Ice Flow over a Viscous Basal Sediment: Theory and Application to Ice Stream B, Antarctica, J. Geophys. Res.Sol. Ea., 94, 4071–4087, https://doi.org/10.1029/JB094iB04p04071, 1989. a, b, c
Maier, N., Humphrey, N., Harper, J., and Meierbachtol, T.: Sliding Dominates SlowFlowing Margin Regions, Greenland Ice Sheet, Sci. Adv., 5, eaaw5406, https://doi.org/10.1126/sciadv.aaw5406, 2019. a, b
Martin, M. A., Winkelmann, R., Haseloff, M., Albrecht, T., Bueler, E., Khroulev, C., and Levermann, A.: The Potsdam Parallel Ice Sheet Model (PISMPIK) – Part 2: Dynamic equilibrium simulation of the Antarctic ice sheet, The Cryosphere, 5, 727–740, https://doi.org/10.5194/tc57272011, 2011. a
Mercer, J. H.: West Antarctic Ice Sheet and CO_{2} Greenhouse Effect: A Threat of Disaster, Nature, 271, 321–325, https://doi.org/10.1038/271321a0, 1978. a
Morland, L. W.: Unconfined IceShelf Flow, in: Dynamics of the West Antarctic Ice Sheet, edited by: Van der Veen, C. J. and Oerlemans, J., Glaciology and Quaternary Geology, Springer Netherlands, 99–116, https://doi.org/10.1007/9789400937451_6, 1987. a
Morlighem, M.: MEaSUREs BedMachine Antarctica, Version 1, NASA National Snow and Ice Data Center Distributed Active Archive Center [data set], https://doi.org/10.5067/C2GFER6PTOS4, 2019. a
Morlighem, M., Seroussi, H., Larour, E., and Rignot, E.: Inversion of Basal Friction in Antarctica Using Exact and Incomplete Adjoints of a HigherOrder Model: Antarctic basal friction inversion, J. Geophys. Res.Earth, 118, 1746–1753, https://doi.org/10.1002/jgrf.20125, 2013. a, b
Morlighem, M., Rignot, E., Mouginot, J., Seroussi, H., and Larour, E.: Deeply Incised Submarine Glacial Valleys beneath the Greenland Ice Sheet, Nat. Geosci., 7, 418–422, https://doi.org/10.1038/ngeo2167, 2014. a
Morlighem, M., Williams, C. N., Rignot, E., An, L., Arndt, J. E., Bamber, J. L., Catania, G., Chauché, N., Dowdeswell, J. A., Dorschel, B., Fenty, I., Hogan, K., Howat, I., Hubbard, A., Jakobsson, M., Jordan, T. M., Kjeldsen, K. K., Millan, R., Mayer, L., Mouginot, J., Noël, B. P. Y., O'Cofaigh, C., Palmer, S., Rysgaard, S., Seroussi, H., Siegert, M. J., Slabon, P., Straneo, F., van den Broeke, M. R., Weinrebe, W., Wood, M., and Zinglersen, K. B.: BedMachine v3: Complete Bed Topography and Ocean Bathymetry Mapping of Greenland From Multibeam Echo Sounding Combined With Mass Conservation, Geophys. Res. Lett., 44, 11051–11061, https://doi.org/10.1002/2017GL074954, 2017. a, b, c
Morlighem, M., Rignot, E., Binder, T., Blankenship, D., Drews, R., Eagles, G., Eisen, O., Ferraccioli, F., Forsberg, R., Fretwell, P., Goel, V., Greenbaum, J. S., Gudmundsson, H., Guo, J., Helm, V., Hofstede, C., Howat, I., Humbert, A., Jokat, W., Karlsson, N. B., Lee, W. S., Matsuoka, K., Millan, R., Mouginot, J., Paden, J., Pattyn, F., Roberts, J., Rosier, S., Ruppel, A., Seroussi, H., Smith, E. C., Steinhage, D., Sun, B., van den Broeke, M. R., van Ommen, T. D., van Wessem, M., and Young, D. A.: Deep Glacial Troughs and Stabilizing Ridges Unveiled beneath the Margins of the Antarctic Ice Sheet, Nat. Geosci., 13, 132–137, https://doi.org/10.1038/s4156101905108, 2020. a, b, c, d, e, f, g
Mouginot, J., Rignot, E., and Scheuchl, B.: Sustained Increase in Ice Discharge from the Amundsen Sea Embayment, West Antarctica, from 1973 to 2013, Geophys. Res. Lett., 41, 1576–1584, https://doi.org/10.1002/2013GL059069, 2014. a
Mouginot, J., Rignot, E., and Scheuchl, B.: ContinentWide, Interferometric SAR Phase, Mapping of Antarctic Ice Velocity, Geophys. Res. Lett., 46, 1–9, https://doi.org/10.1029/2019GL083826, 2019a. a, b, c, d, e
Mouginot, J., Rignot, E., and Scheuchl, B.: MEaSUREs PhaseBased Antarctica Ice Velocity Map, Version 1, NASA National Snow and Ice Data Center Distributed Active Archive Center [data set], https://doi.org/10.5067/PZ3NJ5RXRH10, 2019b. a
Nye, J. F.: The Flow of a Glacier in a Channel of Rectangular, Elliptic or Parabolic CrossSection, J. Glaciol., 5, 661–690, https://doi.org/10.3189/S0022143000018670, 1965. a
Pattyn, F. and Morlighem, M.: The Uncertain Future of the Antarctic Ice Sheet, Science, 367, 1331–1335, https://doi.org/10.1126/science.aaz5487, 2020. a, b, c, d, e
Pattyn, F., Perichon, L., Durand, G., Favier, L., Gagliardini, O., Hindmarsh, R. C., Zwinger, T., Albrecht, T., Cornford, S., Docquier, D., Fürst, J. J., Goldberg, D., Gudmundsson, G. H., Humbert, A., Hütten, M., Huybrechts, P., Jouvet, G., Kleiner, T., Larour, E., Martin, D., Morlighem, M., Payne, A. J., Pollard, D., Rückamp, M., Rybak, O., Seroussi, H., Thoma, M., and Wilkens, N.: GroundingLine Migration in PlanView Marine IceSheet Models: Results of the Ice2sea MISMIP3d Intercomparison, J. Glaciol., 59, 410–422, https://doi.org/10.3189/2013JoG12J129, 2013. a, b, c
Pegler, S. S.: Suppression of Marine Ice Sheet Instability, J. Fluid Mech., 857, 648–680, https://doi.org/10.1017/jfm.2018.742, 2018. a, b, c
Pegler, S. S., Kowal, K. N., Hasenclever, L. Q., and Worster, M. G.: Lateral Controls on GroundingLine Dynamics, J. Fluid Mech., 722, R1, https://doi.org/10.1017/jfm.2013.140, 2013. a
PISM authors: PISM, a Parallel Ice Sheet Model, https://www.pism.io (last access: 22 January 2023), 2021. a, b
Pritchard, H. D., Arthern, R. J., Vaughan, D. G., and Edwards, L. A.: Extensive Dynamic Thinning on the Margins of the Greenland and Antarctic Ice Sheets, Nature, 461, 971–975, https://doi.org/10.1038/nature08471, 2009. a, b
Rayleigh: The Principle of Similitude, Nature, 95, 66–68, https://doi.org/10.1038/095066c0, 1915. a, b
Raymond, C.: Shear Margins in Glaciers and Ice Sheets, J. Glaciol., 42, 90–102, https://doi.org/10.3189/S0022143000030550, 1996. a
Reynolds, O.: An Experimental Investigation of the Circumstances Which Determine Whether the Motion of Water Shall Be Direct or Sinuous, and of the Law of Resistance in Parallel Channels, Philos. T. R. Soc. Lond., 174, 935–982, https://doi.org/10.1098/rstl.1883.0029, 1883. a
Rignot, E. and Mouginot, J.: Ice Flow in Greenland for the International Polar Year 2008–2009, Geophys. Res. Lett., 39, L11501, https://doi.org/10.1029/2012GL051634, 2012. a
Rignot, E., Mouginot, J., Morlighem, M., Seroussi, H., and Scheuchl, B.: Widespread, Rapid Grounding Line Retreat of Pine Island, Thwaites, Smith, and Kohler Glaciers, West Antarctica, from 1992 to 2011, Geophys. Res. Lett., 41, 3502–3509, https://doi.org/10.1002/2014GL060140, 2014. a
Rignot, E., Mouginot, J., Scheuchl, B., van den Broeke, M., van Wessem, M. J., and Morlighem, M.: Four Decades of Antarctic Ice Sheet Mass Balance from 1979–2017, P. Natl. Acad. Sci. USA, 116, 1095–1103, https://doi.org/10.1073/pnas.1812883116, 2019. a
Schoof, C.: The Effect of Cavitation on Glacier Sliding, P. Roy. Soc. A, 461, 609–627, https://doi.org/10.1098/rspa.2004.1350, 2005. a
Schoof, C.: Ice Sheet Grounding Line Dynamics: Steady States, Stability, and Hysteresis, J. Geophys. Res., 112, F03S28, https://doi.org/10.1029/2006JF000664, 2007. a, b, c, d
Schoof, C., Davis, A. D., and Popa, T. V.: Boundary layer models for calving marine outlet glaciers, The Cryosphere, 11, 2283–2303, https://doi.org/10.5194/tc1122832017, 2017. a
Sergienko, O. V. and Hindmarsh, R. C. A.: Regular Patterns in Frictional Resistance of IceStream Beds Seen by Surface Data Inversion, Science, 342, 1086–1089, https://doi.org/10.1126/science.1243903, 2013. a, b, c
Sergienko, O. V. and Wingham, D. J.: Grounding Line Stability in a Regime of Low Driving and Basal Stresses, J. Glaciol., 65, 833–849, https://doi.org/10.1017/jog.2019.53, 2019. a
Sergienko, O. V., Creyts, T. T., and Hindmarsh, R. C. A.: Similarity of Organized Patterns in Driving and Basal Stresses of Antarctic and Greenland Ice Sheets beneath Extensive Areas of Basal Sliding, Geophys. Res. Lett., 41, 3925–3932, https://doi.org/10.1002/2014GL059976, 2014. a, b, c
Seroussi, H. and Morlighem, M.: Representation of basal melting at the grounding line in ice flow models, The Cryosphere, 12, 3085–3096, https://doi.org/10.5194/tc1230852018, 2018. a
Seroussi, H., Nakayama, Y., Larour, E., Menemenlis, D., Morlighem, M., Rignot, E., and Khazendar, A.: Continued Retreat of Thwaites Glacier, West Antarctica, Controlled by Bed Topography and Ocean Circulation: Iceocean modeling of thwaites glacier, Geophys. Res. Lett., 44, 6191–6199, https://doi.org/10.1002/2017GL072910, 2017. a
Seroussi, H., Nowicki, S., Payne, A. J., Goelzer, H., Lipscomb, W. H., AbeOuchi, A., Agosta, C., Albrecht, T., AsayDavis, X., Barthel, A., Calov, R., Cullather, R., Dumas, C., GaltonFenzi, B. K., Gladstone, R., Golledge, N. R., Gregory, J. M., Greve, R., Hattermann, T., Hoffman, M. J., Humbert, A., Huybrechts, P., Jourdain, N. C., Kleiner, T., Larour, E., Leguy, G. R., Lowry, D. P., Little, C. M., Morlighem, M., Pattyn, F., Pelle, T., Price, S. F., Quiquet, A., Reese, R., Schlegel, N.J., Shepherd, A., Simon, E., Smith, R. S., Straneo, F., Sun, S., Trusel, L. D., Van Breedam, J., van de Wal, R. S. W., Winkelmann, R., Zhao, C., Zhang, T., and Zwinger, T.: ISMIP6 Antarctica: a multimodel ensemble of the Antarctic ice sheet evolution over the 21st century, The Cryosphere, 14, 3033–3070, https://doi.org/10.5194/tc1430332020, 2020. a, b
Shepherd, A., Gilbert, L., Muir, A. S., Konrad, H., McMillan, M., Slater, T., Briggs, K. H., Sundal, A. V., Hogg, A. E., and Engdahl, M. E.: Trends in Antarctic Ice Sheet Elevation and Mass, Geophys. Res. Lett., 46, 8174–8183, https://doi.org/10.1029/2019GL082182, 2019. a, b
Szücs, E.: Fundamental Studies in Engineering II, Similitude and Modeling, Elsevier Scientific Publishing Co., Amsterdam, Netherlands, ISBN 9780080983783, 1980. a, b
The IMBIE team: Mass Balance of the Antarctic Ice Sheet from 1992 to 2017, Nature, 558, 219–222, https://doi.org/10.1038/s415860180179y, 2018. a, b
The IMBIE Team: Mass Balance of the Greenland Ice Sheet from 1992 to 2018, Nature, 579, 233–239, https://doi.org/10.1038/s4158601918552, 2020. a, b, c
Tsai, V. C., Stewart, A. L., and Thompson, A. F.: Marine IceSheet Profiles and Stability under Coulomb Basal Conditions, J. Glaciol., 61, 205–215, https://doi.org/10.3189/2015JoG14J221, 2015. a
Weertman, J.: Stability of the Junction of an Ice Sheet and an Ice Shelf, J. Glaciol., 13, 3–11, https://doi.org/10.3189/S0022143000023327, 1974. a
Wessem, J. M. V., Reijmer, C. H., Morlighem, M., Mouginot, J., Rignot, E., Medley, B., Joughin, I., Wouters, B., Depoorter, M. A., Bamber, J. L., Lenaerts, J. T. M., Berg, W. J. V. D., Broeke, M. R. V. D., and Meijgaard, E. V.: Improved Representation of East Antarctic Surface Mass Balance in a Regional Atmospheric Climate Model, J. Glaciol., 60, 761–770, https://doi.org/10.3189/2014JoG14J051, 2014. a, b
Whillans, I. M.: The Equation of Continuity and Its Application to the Ice Sheet Near “Byrd” Station, Antarctica, J. Glaciol., 18, 359–371, https://doi.org/10.3189/S0022143000021055, 1977. a
Winkelmann, R., Martin, M. A., Haseloff, M., Albrecht, T., Bueler, E., Khroulev, C., and Levermann, A.: The Potsdam Parallel Ice Sheet Model (PISMPIK) – Part 1: Model description, The Cryosphere, 5, 715–726, https://doi.org/10.5194/tc57152011, 2011. a
Zhang, T., Price, S. F., Hoffman, M. J., Perego, M., and AsayDavis, X.: Diagnosing the sensitivity of groundingline flux to changes in subiceshelf melting, The Cryosphere, 14, 3407–3424, https://doi.org/10.5194/tc1434072020, 2020. a
 Abstract
 Introduction
 Derivation of the scaling relation
 Testing of the scaling relation
 Application of scaling relation to Antarctic outlet glaciers
 Discussion and conclusions
 Appendix A: Similitude analysis of shallow ice dynamics
 Appendix B: Analysis of Greenland outlet glaciers
 Appendix C: Idealized simulations of MISItype retreat
 Appendix D: Measuring the Antarctic outlets
 Code and data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References
 Supplement
 Abstract
 Introduction
 Derivation of the scaling relation
 Testing of the scaling relation
 Application of scaling relation to Antarctic outlet glaciers
 Discussion and conclusions
 Appendix A: Similitude analysis of shallow ice dynamics
 Appendix B: Analysis of Greenland outlet glaciers
 Appendix C: Idealized simulations of MISItype retreat
 Appendix D: Measuring the Antarctic outlets
 Code and data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References
 Supplement