On non-dimensional forms of basal sliding laws and flow laws for ice-sheet and glacier modelling

For Peer Review THIS MANUSCRIPT HAS BEEN SUBMITTED TO THE JOURNAL OF GLACIOLOGY AND HAS NOT BEEN PEER-REVIEWED. On non-dimensional forms of basal sliding laws and flow laws for ice-sheet and glacier modelling Journal: Journal of Glaciology Manuscript ID JOG-2025-0089.R2 Manuscript Type: Letter Date Submitted by the Author: 22-Oct-2025 Complete List of Authors: Greve, Ralf; Hokkaido University, Institute of Low Temperature Science Keywords: Glacier flow, Glacier modelling, Ice dynamics, Ice-sheet modelling Abstract: Ice sheets and glaciers flow through basal sliding and internal deformation, each governed by physical laws commonly expressed as power-law relationships. These formulations include coefficients - the sliding coefficient and rate factor - whose values and units depend on the respective exponents. This dependency complicates the systematic exploration of parameter space, especially in ensemble simulations. To address this, we propose dimensionless formulations of both sliding and flow laws, in which the coefficients are of order unity and decoupled from the exponents. This separation simplifies sensitivity studies and parameter variations. The dimensionless laws are straightforward to implement in existing models; we demonstrate this with the SICOPOLIS ice-sheet model using three test simulations in an idealized set-up. These simulations illustrate that independent variation of exponents and coefficients is feasible and practical, supporting the use of dimensionless laws in efforts to better constrain ice dynamics in past and future climate scenarios. Cambridge University Press Journal of Glaciology For Peer Review Page 1 of 16 Cambridge University Press Journal of Glaciology For Peer Review Journal of Glaciology, Vol. 71, No. 1, 2025 1 On non-dimensional forms of basal sliding laws and flow1 laws for ice-sheet and glacier modelling2 Ralf GREVE1,2 3 1Institute of Low Temperature Science, Hokkaido University, Sapporo, Japan4 2Arctic Research Center, Hokkaido University, Sapporo, Japan5 Correspondence: Ralf Greve <gr[email protected]>6 ABSTRACT. Ice sheets and glaciers flow through basal sliding and internal7 deformation, each governed by physical laws commonly expressed as power-8 law relations. These formulations include coefficients – the sliding coefficient9 and rate factor – whose values and units depend on the respective exponents.10 This dependency complicates the systematic exploration of parameter space,11 especially in ensemble simulations. To address this, we propose dimensionless12 formulations of both sliding and flow laws, in which the coefficients are of order13 unity and decoupled from the exponents. This separation simplifies sensitivity14 studies and parameter variations. The dimensionless laws are straightforward15 to implement in existing models; we demonstrate this with the SICOPOLIS16 ice-sheet model using three test simulations in an idealized set-up. These17 simulations illustrate that independent variation of exponents and coefficients18 is feasible and practical, supporting the use of dimensionless laws in efforts to19 better constrain ice dynamics in past and future climate scenarios.20 1 INTRODUCTION21 Ice sheets and glaciers flow due to two different processes, namely basal sliding and internal deformation.22 Basal sliding describes the sliding of glacier ice on the underlying substrate, which can be either hard23 bedrock or a deformable sediment layer between ice and bedrock. Internal deformation is governed by the24 Page 2 of 16 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 2 non-linear viscous properties of “hot” polycrystalline ice (that is, with a homologous temperature T{Tm 25 near unity, where Tis the absolute temperature and Tmthe pressure melting point).26 In a dynamic/thermodynamic ice sheet or glacier model, both processes must be included. Basal sliding,27 which in reality is a complex process that depends on a multitude of factors such as the basal temperature,28 roughness of the bedrock, softness of the subglacial sediment layer (if existing) and hydrological conditions,29 is usually parameterized by a sliding law that relates the sliding velocity to the basal stresses. Internal30 deformation can be modelled by a non-linear viscous flow law that describes the relation between the31 macroscopic deformation (strain rate) and internal stresses (e.g., Hooke, 2005; Greve and Blatter, 2009;32 Cuffey and Paterson, 2010).33 Popular forms for such relations are the Weertman–Budd sliding law and the Nye–Glen flow law (see34 below for references). They have in common that they are expressed as power laws with some exponents,35 of which the optimal values are debated, and contain a factor to close the respective equation. This36 factor, the “sliding coefficient” in case of the sliding law and the “rate factor” in case of the flow law,37 may contain remaining dependencies, such as on the temperature. In a dimensional formulation, the units38 and numerical values of these factors depend strongly on the choice of the exponents, which makes it39 cumbersome to vary the exponents over their potential range of values, for instance, within an ensemble40 of simulations for a given scenario. To overcome this obstacle, we propose fully or partly dimensionless41 versions of the sliding and flow laws, which have in common that the respective factor is dimensionless42 and generally of order unity. These formulations decouple the value of the exponents from the value of the43 factor, so that the factors and exponents can be varied independently. We demonstrate this useful feature44 by some simple simulations with the ice-sheet model SICOPOLIS (SImulation COde for POLythermal Ice45 Sheets; SICOPOLIS Authors, 2025).46 2 BASAL SLIDING LAWS47 Basal sliding laws (aka basal friction laws) relate the shear stress (drag) at the base of an ice sheet or glacier, τb, to the basal normal stress, Nband the basal sliding velocity, vb. In a general, implicit form, a basal sliding law can be expressed as fpvb, τb, Nbq “ 0,(1) Page 3 of 16 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 3 where fis a function unspecified at this stage. The list of variables is not necessarily exhaustive as further dependencies, for instance on basal temperature or the presence of basal water, may be included. Note that, in the presence of subglacial water, the basal normal stress is often understood as the difference between the ice overburden stress, Nb,i, and the basal water pressure, pb,w, Nb“Nb,i´pb,w,(2) and then called the reduced normal stress or, alternatively, the effective pressure.48 A popular form of a basal sliding law is the Weertman–Budd sliding law, which results when assuming an explicit form of Eq. (1) solved for vb, with power-law dependencies on τband Nb: vb“Cb τp b Nq b ,(3) where Cbis the sliding coefficient and pp, qqare the non-negative sliding exponents (Weertman, 1957; Budd and others, 1979, 1984; Budd and Jenssen, 1987). Alternatively, Eq. (3) can be solved for the shear stress, τb“C‹ bv1{p bNq{p b,with C‹ b“C´1{p b,(4) where C‹ bis the basal friction coefficient. For the case q“0, that is, ignoring the dependence on the49 normal stress Nb, the above forms are often referred to as the Weertman sliding law.50 In principle, vband τbare vector quantities. For simplicity, we formulate the sliding laws only with the51 respective magnitudes. However, to interpret the results correctly, it must be kept in mind that vband τb 52 are anti-parallel to each other due to the nature of friction.53 Let us now non-dimensionalize the sliding law (3) by introducing scales (typical values) for the relevant quantities (e.g., Hutter and Jöhnk, 2004). We consider a situation near the edge of an ice sheet where Page 4 of 16 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 4 basal sliding is most relevant: rHs “ 1 km ptypical thicknessq,(5a) ε“10´2ptypical surface slopeq,(5b) rNbs “ ρgrHs “ 107Pa ptypical normal stressq,(5c) rτbs “ εrNbs “ 105Pa ptypical shear stressq,(5d) rvbs “ 100 m a´1ptypical sliding velocityq,(5e) where we have used approximate values ρ«103kg m´3for the ice density and g«10 m s´2for the acceleration due to gravity, which is sufficiently accurate for the sake of a scaling analysis. This scaling is consistent with the linear sliding law [Cb“10´3m a´1Pa´1,pp, qq “ p1,0q] used for the EISMINT Phase 2 Simplified Geometry Experiments (Payne and others, 2000): 100 m a´1 loooomoooon vb “10´3m a´1Pa´1 looooooooomooooooooon Cb ˆ105Pa loomoon τb .(6) An appropriate choice for the scale of the sliding coefficient results from Eq. (3) as rCbs “ rvbsrNbsq rτbsp.(7) We now use the above scales to introduce dimensionless quantities as follows: vb“ rvbs˜vb,(8a) τb“ rτbs˜τb,(8b) Nb“ rNbs˜ Nb,(8c) Cb“ rCbs˜ Cb,(8d) where the quantities marked by the tilde symbol are the non-dimensional basal sliding velocity, shear stress, normal stress and sliding coefficient, respectively. Inserting Eqs. (8) in the sliding law (3) yields its fully non-dimensional form, ˜vb“˜ Cb ˜τp b ˜ Nq b ,(9) Page 5 of 16 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 5 in which all quantities are supposed to be of order unity.54 A dimensional form of Eq. (3) can be kept by only making use of the scaling (8d) of the sliding coefficient, vb“ rCbs˜ Cb τp b Nq b ,(10) which has the advantage that its implementation in an existing model based on dimensional quantities55 requires only minimal adaptations.56 In order to obtain the fully or partly dimensionless counterparts of Eq. (4), we note the scaling and non-dimensionalization of the friction coefficient C‹ b: C‹ b“ rC‹ bs˜ C‹ b,with rC‹ bs “ rCbs´1{p“rτbs rvbs1{prNbsq{p.(11) The fully non-dimensional form of Eq. (4) results then as ˜τb“˜ C‹ b˜v1{p b˜ Nq{p b,(12) and the dimensional form in which only the scaling (11) of the friction coefficient is used reads τb“ rC‹ bs˜ C‹ bv1{p bNq{p b.(13) Why do we promote using Eqs. (10) or (13) instead of Eqs. (3) or (4) in an ice sheet or glacier57 model? In Table 1 we have compiled some parameter combinations that were used along with Weertman58 or Weertman–Budd sliding laws in the literature. The impossibility of comparing the various dimensional59 sliding coefficients Cbfor different exponents pp, qqbecomes immediately evident. They do not even have60 a common unit, and the respective numerical value tells nothing about the actual strength of basal sliding.61 In the second case, pp, qq“p1,2q, the numerical value of Cbis greater than 109; however, the small62 dimensionless value means that it produces only very little basal sliding ( ˜ Cb«0.04). By contrast, in the63 third case, pp, qq “ p3,0q, the numerical value of Cbis merely 10´12; however, it corresponds to pronounced64 basal sliding ( ˜ Cb“10). The dimensionless sliding coefficients ˜ Cbgive a much better idea about what the65 respective value means physically, and allow comparing values across different sliding laws.66 To further strengthen our point, suppose that we wish to test a sliding law with a new set of exponents,67 for instance pp, qq “ p3,1.5q. Working with the dimensional sliding coefficient Cb, we would not have any68 Page 6 of 16 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 6 pp, qqCbrCbs˜ CbReference p1,0q10´3m a´1Pa´110´3m a´1Pa´11Payne and others (2000) p1,2q3.985 ˆ109m a´1Pa 1011 m a´1Pa 0.03985 Budd and others (1984): p3,0q10´12 m a´1Pa´310´13 m a´1Pa´310 Cornford and others (2020) p3,1q1.607 ˆ10´6m a´1Pa´210´6m a´1Pa´21.607 Saito and others (2016): p3,2q6.72 m a´1Pa´110 m a´1Pa´10.672 Rückamp and others (2019) Table 1. Sliding exponents pp, qq, dimensional sliding coefficients Cb, scales rCbsand dimensionless sliding coefficients ˜ Cbfor several Weertman (q“0) or Weertman–Budd (qą0) sliding laws used in the literature. :: Rather than using the normal stress Nb, these sliding laws were formulated with the pressure head Z“Nb{pρgq. We converted the sliding coefficients given in these studies accordingly, using ρ“910 kg m´3and g“9.81 m s´2. idea which order of magnitude may be suited for its numerical value, and which range of values mean69 strong or weak sliding. However, if the dimensionless sliding coefficient ˜ Cbis used, we can immediately70 start with an initial guess ˜ Cb“1and, from there on, refine the sliding law by, e.g., tuning to observed flow71 speeds. According to the scaling (7), (8d), the dimensional equivalent of ˜ Cb“1would be Cb“ rCbs “72 10´2.5m a´1Pa´1.5“3.162 ˆ10´3m a´1Pa´1.5.73 We have only discussed cases with a constant sliding parameter; however, the non-dimensionalization74 method is of course not limited to this. It can also be applied to a spatially variable sliding coefficient, which75 may arise from an inversion procedure (e.g., Morlighem and others, 2013). Alternative sliding laws, such as76 the Coulomb-limited rules discussed by Cornford and others (2020), allow similar non-dimensionalization,77 although we refrain from working out the details here.78 3 FLOW LAWS79 A similar problem of units and hugely varying numerical values arises for the flow law of polycrystalline ice. It is a viscous flow law that relates the strain-rate (stretching) tensor dij to the stress deviator tD ij. The strain-rate tensor is defined as dij “1 2ˆBvi Bxj `Bvj Bxi˙pi, j “1,2,3q,(14) Page 7 of 16 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 7 where xidenotes the Cartesian coordinates (x1“x,x2“y,x3“z), and viis the velocity vector. The stress deviator is the traceless part of the full stress tensor tij, tij “ ´p δij `tD ij ,(15) where p“ ´tii{3is the pressure (we assume the Einstein summation convention: summation over the80 twice-appearing index iimplied, thus tii is the trace of the stress tensor), and δij is the Kronecker delta81 symbol, in other words, the unit tensor in index notation.82 For the flow law, usually collinearity between the symmetric tensors dij and tD ij is assumed. We note the form given by Greve and Blatter (2009), dij “AfpτeqtD ij ,(16) where Ais the rate factor, τe“ rptD ijtD ijq{2s1{2the effective stress (summation over both iand jimplied), and fpτeqis the creep function. The rate factor depends on the temperature relative to the pressure melting point, T1, via an Arrhenius law (e.g., Cuffey and Paterson, 2010), but it is sometimes chosen as a constant parameter for simplicity. In the Nye–Glen flow law (Glen, 1955; Nye, 1957), the creep function is expressed as a power law, fpτeq “ τn´1 e,(17) so that dij “Aτn´1 etD ij ,(18) where nis the stress exponent. A value of n“1would correspond to a Newtonian fluid; however, the83 deformability of ice differs markedly from that behaviour, and the value is frequently chosen as n“3,84 or within the range from 1.5 to 4.2 (Cuffey and Paterson, 2010) (while recent evidence from laboratory85 experiments actually supports n“1for temperate ice; Schohn and others, 2025).86 The Nye–Glen flow law (18) can also be inverted for the stress deviator, tD ij “A‹d´p1´1{nq edij ,with A‹“A´1{n,(19) where A‹is the associated rate factor and de“ rpdijdijq{2s1{2the effective strain rate (e.g., Greve and87 Page 8 of 16 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 14 REFERENCES175 Bons PD, Kleiner T, Llorens MG, Prior DJ, Sachau T, Weikusat I and Jansen D (2018) Greenland ice sheet: Higher176 nonlinearity of ice flow significantly reduces estimated basal motion. Geophysical Research Letters,45(13), 6542–177 6548 (doi: 10.1029/2018GL078356)178 Budd WF and Jenssen D (1987) Numerical modelling of the large-scale basal water flux under the West Antarctic ice179 sheet. In CJ van der Veen and J Oerlemans (eds.), Dynamics of the West Antarctic Ice Sheet, 293–320, D. Reidel180 Publishing Company, Dordrecht, The Netherlands (doi: 10.1007/978-94-009-3745-1)181 Budd WF, Keage PL and Blundy NA (1979) Empirical studies of ice sliding. Journal of Glaciology,23(89), 157–170182 (doi: 10.3189/S0022143000029804)183 Budd WF, Jenssen D and Smith IN (1984) A three-dimensional time-dependent model of the West Antarctic ice184 sheet. Annals of Glaciology,5, 29–36 (doi: 10.3189/1984AoG5-1-29-36)185 Cornford SL, Seroussi H, Asay-Davis XS, Gudmundsson GH, Arthern R, Borstad C, Christmann J, Dias dos Santos186 T, Feldmann J, Goldberg D, Hoffman MJ, Humbert A, Kleiner T, Leguy G, Lipscomb WH, Merino N, Durand187 G, Morlighem M, Pollard D, Rückamp M, Williams CR and Yu H (2020) Results of the third Marine Ice Sheet188 Model Intercomparison Project (MISMIP+). The Cryosphere,14(7), 2283–2301 (doi: 10.5194/tc-14-2283-2020)189 Cuffey KM and Paterson WSB (2010) The Physics of Glaciers. Elsevier, Amsterdam, The Netherlands etc., 4th190 edition, ISBN 978-0-12-369461-4191 Dunse T, Greve R, Schuler TV and Hagen JO (2011) Permanent fast flow versus cyclic surge behaviour:192 numerical simulations of the Austfonna ice cap, Svalbard. Journal of Glaciology,57(202), 247–259 (doi:193 10.3189/002214311796405979)194 Fowler AC and Johnson C (1996) Ice-sheet surging and ice-stream formation. Annals of Glaciology,23, 68–73 (doi:195 10.3189/S0260305500013276)196 Getraer B and Morlighem M (2025) Increasing the Glen-Nye power-law exponent accelerates ice-loss projections197 for the Amundsen Sea Embayment, West Antarctica. Geophysical Research Letters,52(7), e2024GL112516 (doi:198 10.1029/2024GL112516)199 Glen JW (1955) The creep of polycrystalline ice. Proceedings of the Royal Society A,228(1175), 519–538 (doi:200 10.1098/rspa.1955.0066)201 Goldberg DN (2011) A variationally derived, depth-integrated approximation to a higher-order glaciological flow202 model. Journal of Glaciology,57(201), 157–170 (doi: 10.3189/002214311795306763)203 Page 15 of 16 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 15 Grandadam F (2024) Implementation of the Depth Integrated Viscosity Approximation in SICOPOLIS. Internship204 Report, Claude Bernard University Lyon 1, France, and Hokkaido University, Sapporo, Japan (doi: 10.5281/zen-205 odo.14732938)206 Greve R (2005) Relation of measured basal temperatures and the spatial distribution of the geothermal heat flux for207 the Greenland ice sheet. Annals of Glaciology,42, 424–432 (doi: 10.3189/172756405781812510)208 Greve R and Blatter H (2009) Dynamics of Ice Sheets and Glaciers. Springer, Berlin, Germany etc., ISBN 978-3-209 642-03414-5 (doi: 10.1007/978-3-642-03415-2)210 Greve R and Blatter H (2016) Comparison of thermodynamics solvers in the polythermal ice sheet model SICOPOLIS.211 Polar Science,10(1), 11–23 (doi: 10.1016/j.polar.2015.12.004)212 Hindmarsh RCA (2009) Consistent generation of ice-streams via thermo-viscous instabilities modulated by membrane213 stresses. Geophysical Research Letters,36(6), L06502 (doi: 10.1029/2008GL036877)214 Hindmarsh RCA and Le Meur E (2001) Dynamical processes involved in the retreat of marine ice sheets. Journal of215 Glaciology,47(157), 271–282 (doi: 10.3189/172756501781832269)216 Hooke RL (2005) Principles of Glacier Mechanics. Cambridge University Press, Cambridge, UK and New York, NY,217 USA, 2nd edition, ISBN 9780511614231 (doi: 10.1017/CBO9780511614231)218 Hutter K and Jöhnk K (2004) Continuum Methods of Physical Modeling. Springer, Berlin, Germany etc. (doi:219 10.1007/978-3-662-06402-3)220 Lipscomb WH, Price SF, Hoffman MJ, Leguy GR, Bennett AR, Bradley SL, Evans KJ, Fyke JG, Kennedy JH,221 Perego M, Ranken DM, Sacks WJ, Salinger AG, Vargo LJ and Worley PH (2019) Description and evaluation of the222 Community Ice Sheet Model (CISM) v2.1. Geoscientific Model Development,12(1), 387–424 (doi: 10.5194/gmd-223 12-387-2019)224 Millstein JD, Minchew BM and Pegler SS (2022) Ice viscosity is more sensitive to stress than commonly assumed.225 Communications Earth & Environment,3, 57 (doi: 10.1038/s43247-022-00385-x)226 Morlighem M, Seroussi H, Larour E and Rignot E (2013) Inversion of basal friction in Antarctica using exact and227 incomplete adjoints of a higher-order model. Journal of Geophysical Research: Earth Surface,118(3), 1746–1753228 (doi: 10.1002/jgrf.20125)229 Nye JF (1957) The distribution of stress and velocity in glaciers and ice sheets. Proceedings of the Royal Society A,230 239(1216), 113–133 (doi: 10.1098/rspa.1957.0026)231 Page 16 of 16 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 16 Payne AJ and Dongelmans PW (1997) Self-organization in the thermomechanical flow of ice sheets. Journal of232 Geophysical Research: Solid Earth,102(B6), 12219–12233 (doi: 10.1029/97JB00513)233 Payne AJ, Huybrechts P, Abe-Ouchi A, Calov R, Fastook JL, Greve R, Marshall SJ, Marsiat I, Ritz C, Tarasov L234 and Thomassen MPA (2000) Results from the EISMINT model intercomparison: the effects of thermomechanical235 coupling. Journal of Glaciology,46(153), 227–238 (doi: 10.3189/172756500781832891)236 Rückamp M, Greve R and Humbert A (2019) Comparative simulations of the evolution of the Greenland ice sheet237 under simplified Paris Agreement scenarios with the models SICOPOLIS and ISSM. Polar Science,21, 14–25 (doi:238 10.1016/j.polar.2018.12.003)239 Saito F, Abe-Ouchi A, Takahashi K and Blatter H (2016) SeaRISE experiments revisited: potential sources of spread240 in multi-model projections of the Greenland ice sheet. The Cryosphere,10(1), 43–63 (doi: 10.5194/tc-10-43-2016)241 Sayag R and Tziperman E (2008) Spontaneous generation of pure ice streams via flow instability: Role of longi-242 tudinal shear stresses and subglacial till. Journal of Geophysical Research: Solid Earth,113(B5), B05411 (doi:243 10.1029/2007JB005228)244 Schohn CM, Iverson NR, Zoet LK, Fowler JR and Morgan-Witts N (2025) Linear-viscous flow of temperate ice.245 Science,387(6730), 182–185 (doi: 10.1126/science.adp7708)246 SICOPOLIS Authors (2025) SICOPOLIS v25. GitHub, URL https://github.com/sicopolis/sicopolis247 Weertman J (1957) On the sliding of glaciers. Journal of Glaciology,3(21), 33–38 (doi: 10.3189/S0022143000024709)248 Page 17 of 16 Cambridge University Press Journal of Glaciology