the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Brief communication: Is vertical shear in an ice shelf (still) negligible?
Chris Miele
Timothy C. Bartholomaus
Ellyn M. Enderlin
Vertical shear is recognized today as a key component of the stress balance of ice shelves. However, the first ice shelf models were built on the neglect of vertical shear. Partly due to its historical treatment, it remains common to discuss vertical shear as though it were still considered negligible in ice shelf models. Here, we offer a historical perspective on the changing treatment of vertical shear over time, and we emphasize the term's nonnegligibility in current ice shelf modeling. We illustrate our discussion in the simplest context of an analytic, isothermal, shallowiceshelf model.
Please read the corrigendum first before continuing.

Notice on corrigendum
The requested paper has a corresponding corrigendum published. Please read the corrigendum first before downloading the article.

Article
(1336 KB)

The requested paper has a corresponding corrigendum published. Please read the corrigendum first before downloading the article.
 Article
(1336 KB)  Fulltext XML
 Corrigendum
 BibTeX
 EndNote
Analytic models of floating ice shelves date back to at least 1957, when Johannes Weertman derived expressions for the tension and velocity gradients within a uniformthickness ice shelf (Weertman, 1957). Weertman found that, for a shelf with uniform surface elevation h=h_{T}; uniform density ρ; and no lateral flow, the depthaveraged longitudinal deviatoric tension, $\stackrel{\mathrm{\u203e}}{{\mathit{\tau}}_{xx}}$, could be calculated via
Using a depthaveraged constitutive relation, expressions of this form permit the calculation of strain rates and velocities. Nearly 2 decades later, Thomas (1973) set out to generalize Weertman's expression for shelves of nonuniform thickness. Using the same underlying assumptions as Weertman but imposing no restrictions on the surface elevation h=h(x), Thomas obtained an expression nearly identical to Weertman's, wherein the depthaveraged deviatoric tension is
By Thomas' analysis, Weertman's solution is valid regardless of how h varies along a shelf. Thomas’ expression remains the generally accepted description of a nonuniformthickness shelf in longitudinal extension, and it is routinely cited or independently derived in the literature (Sanderson, 1979; Cuffey and Paterson, 2010; Gudmundsson, 2013; Oerlemans, 2021; Millstein et al., 2022). However, though Eq. (2) has persisted, the formulation of this model has quietly undergone a conceptual shift over the decades. This conceptual shift relates to the role of vertical shear in ice shelves – a topic which is sometimes incompletely communicated today, and on which we seek to provide clarification.
In originally deriving the nonuniformthickness model of Eq. (2), Thomas' “sole restriction [was] that of zero shear stresses in vertical planes.” The neglect of vertical shear stress, τ_{xz}, was universal in the formulation of ice shelf models at the time (Weertman, 1957; Thomas, 1973; Robin, 1975; Sanderson, 1979). However, it was understood by some authors to be theoretically suspect. In a seminal paper titled “Is Vertical Shear in an Ice Shelf Negligible?”, Sanderson and Doake (1979) argued that vertical shear is fundamentally linked with the thickness gradient of an ice shelf and that, strictly speaking, vertical shear “cannot be precisely zero” except in the case of uniform thickness. This observation did not challenge the practical utility of Eq. (2) (Sanderson and Doake, 1979, found vertical shear to be small enough that its neglect was, in fact, justified, answering their own titular question in the affirmative), but it highlighted a relationship that had been missed in Thomas' analysis.
The formulation of Thomas' model evolved with the development of the shallowshelf approximation (SSA) (Morland, 1987; MacAyeal, 1989; Schoof and Hindmarsh, 2010). The SSA, besides empowering a leap forward in computational glaciology, was accompanied by two key theoretical advances in ice shelf modeling: how the neglect of terms is justified and which terms are neglected.
The SSA is built on the fundamental assumption that the thicknesstolength aspect ratio, ϵ, of an ice shelf is small (this is the “shallowness” of the SSA). With ϵ≪1, larger powers of the aspect ratio attain smaller values. In contemporary terminology, an “nthorder approximation” is obtained by neglecting any term appearing as a coefficient of ϵ^{n+1} after nondimensionalization (i.e., setting those terms to zero in the approximate model, with the understanding that they are not exactly zero in real life). This dimensionalanalysis approach to excluding terms adds quantitative rigour to the approach of Sanderson and Doake (1979), for whom negligibility was more qualitatively assessed. As we sketch out in the next section and contrary to the postulates of the first ice shelf modelers, dimensional analysis does not lead to the wholesale neglect of vertical shear from shallowshelf models. Using the SSA as a starting point to derive an analytic model for a longitudinally extending ice shelf, Thomas' Eq. (2) results, albeit without the assumption of vanishing vertical shear. This is the modern approach to deriving Eq. (2).
However, the inclusion of vertical shear stress in the presentday interpretation of Thomas' model may somewhat clash with intuition, not least because the vertical shear term does not actually appear anywhere in Eq. (2). Even in literature postdating the development of the SSA, it is common to encounter language which, to a novice glaciologist, might seem to imply that vertical shear is still discarded entirely from shallowiceshelf models.^{1} This potential misstep is the primary motivation for the present paper. In the sections below, we briefly illustrate (a) that dimensional analysis does not entail the neglect of vertical shear from the stress balance of even a zeroorder shallowiceshelf model and (b) that, despite the term's absence in Eq. (2), the vertical shear stress of an ice shelf can be directly calculated. For compactness, we present this discussion in the simplest context of an isothermal, unconfined ice shelf in one horizontal dimension (as depicted in Fig. 1).
The typical balance of momentum for a 2D ice shelf crosssection, in x and z, can be expressed as
where each τ_{ij} term is a deviatoric stress and P is pressure. Both historically and today, simplification of these equations has typically been carried out by neglecting the “bridging term” ($\frac{\partial}{\partial x}{\mathit{\tau}}_{xz}$ of Eq. 3b). Possibly in an effort to maintain internal consistency, the pioneering authors discussed above tended to additionally neglect all other appearances of τ_{xz}, including the second term of Eq. (3a).
In contrast to their approach and with asterisks denoting appropriately scaled parameters, Eq. (3a)–(3b) can be nondimensionalized as in Weis et al. (1999) (see also MacAyeal, 1989, or Schoof and Hindmarsh, 2010) to become
By this approach, in the zeroth and firstorder approximations, the bridging term in Eq. (4b) will be neglected as a coefficient of ϵ^{2} (i.e., omitted for the purpose of further simplification), in agreement with Weertman, Thomas, and others. However, the vertical shear term in Eq. (4a) is retained, as it represents a coefficient of ϵ^{0}. Thus, dimensional analysis provides an internally consistent means of neglecting $\frac{\partial}{\partial x}{\mathit{\tau}}_{xz}$ while retaining $\frac{\partial}{\partial z}{\mathit{\tau}}_{xz}$.
Following the workflow presented in Sect. 5.2 of Greve and Blatter (2009), the neglect of $\frac{\partial}{\partial x}{\mathit{\tau}}_{xz}$ in Eq. (4b) yields the modified xmomentum equation shown below (where we have omitted asterisks for readability).
Although the most general shallowshelf models are constructed by depthintegrating equations of the above form (resulting in, for example, Eq. 6.55 of Greve and Blatter, 2009), we can provide an even simpler solution for an isothermal shelf, for which viscosity is depth invariant, and, consequently, τ_{xx} is depth invariant and equal to its depthaveraged value, $\stackrel{\mathrm{\u203e}}{{\mathit{\tau}}_{xx}}$. It can be verified, by direct substitution (see also Schoof and Hindmarsh, 2010), that Eq. (5) is solved by the system
where Eq. (6a) is simply the isothermal case of Thomas' Eq. (2), and Eq. (6b) is the (nonzero!) vertical shear stress which must accompany that solution. In fact, if we had set τ_{xz} to 0, Eq. (6a) would not actually solve Eq. (5), unless we additionally assumed that $\frac{\partial}{\partial x}h=\mathrm{0}$. This observation succinctly illustrates the work of Sanderson and Doake (1979): τ_{xz}=0 only to the extent that a shelf has uniform thickness.
In discussions of shallowiceshelf models, it is fairly common to hear vertical shear spoken of as “zero”, “neglected”, or otherwise unimportant. While this certainly was an approximation made by early ice shelf modelers, this language is at odds with current modeling practice. Indeed, as first shown by Sanderson and Doake (1979) in “Is Vertical Shear in an Ice Shelf Negligible?”, nonzero vertical shear stress is a fundamental requirement of a nonuniformthickness ice shelf. With dimensional analysis now enabling modelers to more rigorously define negligibility, the presentday answer to their question is a resounding “no”.
No data sets were used in this article.
CM devised the project and wrote the manuscript, with guidance from TCB and EME. TCB and EME both reviewed and edited the manuscript and acquired the financial support necessary for the project. TCB directly supervised CM.
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.
This work was made possible by grant nos. 1716865 and 1933105 from the US National Science Foundation and grant no. 80NSSC18K1477 from NASA.
This research has been supported by the National Science Foundation (grant nos. 1716865 and 1933105) and the National Aeronautics and Space Administration (grant no. 80NSSC18K1477).
This paper was edited by Reinhard Drews and reviewed by four anonymous referees.
Bondzio, J. H., Seroussi, H., Morlighem, M., Kleiner, T., Rückamp, M., Humbert, A., and Larour, E. Y.: Modelling calving front dynamics using a levelset method: application to Jakobshavn Isbræ, West Greenland, The Cryosphere, 10, 497–510, https://doi.org/10.5194/tc104972016, 2016. a
Bueler, E. and Brown, J.: Shallow shelf approximation as a “sliding law” in a thermomechanically coupled ice sheet model, J. Geophys. Res.Sol. Ea., 114, 1–21, https://doi.org/10.1029/2008JF001179, 2009. a
Cuffey, K. M. and Paterson, W. S. B.: The Physics of Glaciers, ButterworthHeinemann, Burlington, MA, 4th edn., https://doi.org/10.3189/002214311796405906, 2010. a, b
Greve, R. and Blatter, H.: Dynamics of Ice Sheets and Glaciers, Springer, Berlin, https://doi.org/10.3189/002214311798043717, 2009. a, b
Gudmundsson, G. H.: Iceshelf buttressing and the stability of marine ice sheets, The Cryosphere, 7, 647–655, https://doi.org/10.5194/tc76472013, 2013. a
Larour, E., Seroussi, H., Morlighem, M., and Rignot, E.: Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM), J. Geophys. Res.Earth, 117, F01022, https://doi.org/10.1029/2011JF002140, 2012. a
MacAyeal, D. R.: Largescale ice flow over a viscous basal sediment: theory and application to Ice Stream B, Antarctica, J. Geophys. Res., 94, 4071–4087, 1989. a, b
Millstein, J. D., Minchew, B. M., and Pegler, S. S.: Ice viscosity is more sensitive to stress than commonly assumed, Communications Earth & Environment, 3, 1–7, https://doi.org/10.1038/s4324702200385x, 2022. a
Morland, L. W.: Unconfined iceshelf flow, in: Dynamics of the West Antarctic Ice Sheet, 28, 5241–5244, https://doi.org/10.1007/9789400937451_6, 1987. a
Oerlemans, J.: Analytical Models of Ice Sheets and Ice Shelves, in: Glaciers and Ice Sheets in the Climate System, edited by: Fowler, A. and Ng, F., Springer International Publishing, 241–254, https://doi.org/10.1007/9783030425845_10, 2021. a
Pattyn, F. and Decleir, H.: Numerical simulation of Shirase Glacier, East Queen Maud Land, Antarctica, NIPR Symp. Polar Meteorol. Glaciol., 9, 87–109, https://doi.org/10.15094/00003882, 1995. a
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
Robin, G. D. Q.: Ice shelves and ice flow, Nature, 253, 168–172, https://doi.org/10.1038/253168a0, 1975. a
Rückamp, M., Neckel, N., Berger, S., Humbert, A., and Helm, V.: Calving Induced Speedup of Petermann Glacier, J. Geophys. Res.Earth, 124, 216–228, https://doi.org/10.1029/2018JF004775, 2019. a
Sanderson, T. J.: Equilibrium profile of ice shelves, J. Glaciol., 22, 435–460, https://doi.org/10.3189/S0022143000014453, 1979. a, b
Sanderson, T. J. and Doake, C. S.: Is Vertical Shear in an Ice Shelf Negligible?, J. Glaciol., 22, 285–292, https://doi.org/10.1017/S0022143000014271, 1979. a, b, c, d, e
Schoof, C. and Hindmarsh, R. C.: Thinfilm flows with wall slip: An asymptotic analysis of higher order glacier flow models, Q. J. Mech. Appl. Math., 63, 73–114, https://doi.org/10.1093/qjmam/hbp025, 2010. a, b, c
Thomas, R. H.: The creep of ice shelves: theory, J. Glaciol., 12, 45–53, https://doi.org/10.3189/S0022143000022693, 1973. a, b
Weertman, J.: Deformation of floating ice shelves, J. Glaciol., 3, 38–42, https://doi.org/10.3189/S0022143000024710, 1957. a, b
Weis, M., Greve, R., and Hutter, K.: Theory of shallow ice shelves, Continuum Mech. Therm., 11, 15–50, 1999. a
For example, in constructing the ice shelf model of Pattyn and Decleir (1995), “the [vertical] shear stress term in [the xmomentum equation] is omitted.” Bueler and Brown (2009) state that “$D(\mathit{v}{)}_{\mathrm{13}},D(\mathit{v}{)}_{\mathrm{31}},D(\mathit{v}{)}_{\mathrm{23}},D(\mathit{v}{)}_{\mathrm{32}}$ [i.e., the vertical shear strain rates] are all negligible in the SSA.” Cuffey and Paterson (2010) specify that, to construct a nonuniformthickness ice shelf model, the “assumption must be made that the slope at the bottom surface of the shelf is small so that the stress τ_{xz} will be negligible.” Larour et al. (2012) introduce the SSA as obtained by “assuming that vertical shear is negligible” and then specify that ${\dot{\mathit{\epsilon}}}_{xz}={\dot{\mathit{\epsilon}}}_{yz}$ = 0. In an ice shelf model intercomparison, Pattyn et al. (2013) write, “A further approximation, known as the shallowshelf approximation (SSA), is obtained by neglecting vertical shear.” Bondzio et al. (2016) describe the SSA as an approximation which “neglects all vertical shearing but includes membrane stresses”, and Rückamp et al. (2019) affirm that “the SSA neglects vertical shearing.”
 Abstract
 Introduction
 The conceptual evolution of Thomas' model
 Preconstruction: which terms are neglected from shallowiceshelf models by dimensional analysis?
 Constructing the simplest shallowshelf model
 Concluding remarks
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References
The requested paper has a corresponding corrigendum published. Please read the corrigendum first before downloading the article.
 Article
(1336 KB)  Fulltext XML
negligiblein ice shelf models. We address this miscommunication, providing conceptual guidance regarding this often misrepresented stress. Fundamentally, vertical shear is required to balance thickness gradients in ice shelves.
 Abstract
 Introduction
 The conceptual evolution of Thomas' model
 Preconstruction: which terms are neglected from shallowiceshelf models by dimensional analysis?
 Constructing the simplest shallowshelf model
 Concluding remarks
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References