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.R1 Manuscript Type: Letter Date Submitted by the Author: 08-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 multiplicative 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 15 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-law8 relations. These formulations include multiplicative coefficients – the sliding9 coefficient and rate factor – whose values and units depend on the respec-10 tive exponents. This dependency complicates the systematic exploration of11 parameter space, especially in ensemble simulations. To address this, we pro-12 pose dimensionless formulations of both sliding and flow laws, in which the13 coefficients are of order unity and decoupled from the exponents. This separa-14 tion simplifies sensitivity studies and parameter variations. The dimensionless15 laws are straightforward to implement in existing models; we demonstrate16 this with the SICOPOLIS ice-sheet model using three test simulations in an17 idealized set-up. These simulations illustrate that independent variation of18 exponents and coefficients is feasible and practical, supporting the use of di-19 mensionless laws in efforts to better constrain ice dynamics in past and future20 climate scenarios.21 1 INTRODUCTION22 Ice sheets and glaciers flow due to two different processes, namely basal sliding and internal deformation.23 Basal sliding describes the sliding of glacier ice on the underlying substrate, which can be either hard24 bedrock or a deformable sediment layer between ice and bedrock. Internal deformation results from the25 Page 2 of 15 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 2 viscous properties of “hot” polycrystalline ice (that is, with a homologous temperature T{Tmnear unity,26 where Tis the absolute temperature and Tmthe pressure melting point).27 In a dynamic/thermodynamic ice sheet or glacier model, both processes must be included. Basal sliding,28 which in reality is a complex process that depends on a multitude of factors such as the basal temperature,29 roughness of the bedrock, softness of the sediment layer (if existing) and hydrological conditions, is usually30 parameterized by a sliding law that relates the sliding velocity to the basal stresses. Internal deformation31 can be modelled by a non-linear viscous flow law that describes the relation between the macroscopic32 deformation (strain rate) and internal stresses (e.g., Hooke, 2005; Greve and Blatter, 2009; Cuffey and33 Paterson, 2010).34 Popular forms for such relations are the Weertman–Budd sliding law and the Nye–Glen flow law (see35 below for references). They have in common that they are expressed as power laws with some exponents,36 of which the optimal values are debated, and contain a factor to close the respective equation. This37 factor, the “sliding coefficient” in case of the sliding law and the “rate factor” in case of the flow law,38 may contain remaining dependencies, such as on the temperature. In a dimensional formulation, the units39 and numerical values of these factors depend strongly on the choice of the exponents, which makes it40 cumbersome to vary the exponents, for instance, within an ensemble of simulations for a given scenario.41 To overcome this obstacle, we propose fully or partly dimensionless versions of the sliding and flow laws,42 which have in common that the respective factor is dimensionless and generally of order unity. These43 formulations decouple the value of the exponents from the value of the factor, so that an independent44 variation of the exponents and factors becomes possible. We demonstrate this useful feature by some45 simple simulations with the ice-sheet model SICOPOLIS (SImulation COde for POLythermal Ice Sheets;46 SICOPOLIS Authors, 2025).47 2 BASAL SLIDING LAWS48 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 15 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.49 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 the50 normal stress Nb, the above forms are often referred to as the Weertman sliding law.51 In principle, vband τbare vector quantities. For simplicity, we formulate the sliding laws only with the52 respective magnitudes. However, to interpret the results correctly, it must be kept in mind that vband τb 53 are anti-parallel to each other due to the nature of friction.54 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 basal sliding is most relevant: Typical thickness rHs “ 1 km , typical surface slope ε“10´2, typical normal stress rNbs “ ρgrHs “ 107Pa , typical shear stress rτbs “ εrNbs “ 105Pa , typical sliding velocity rvbs “ 100 m a´1, (5) Page 4 of 15 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 4 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, τb“ rτbs˜τb, Nb“ rNbs˜ Nb, Cb“ rCbs˜ Cb,(8) where the quantities marked by the tilde symbol are the non-dimensional basal sliding velocity, shear stress, normal stress and sliding coefficient, respectively. Inserting Eq. (8) in the sliding law (3) yields its fully non-dimensional form, ˜vb“˜ Cb ˜τp b ˜ Nq b ,(9) in which all quantities are supposed to be of order unity.55 A dimensional form of Eq. (3) can be kept by only making use of the scaling (8)4of 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 quantities56 requires only minimal adaptations.57 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) Page 5 of 15 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 5 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. 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 glacier58 model? In Table 1 we have compiled some parameter combinations that were used along with Weertman59 or Weertman–Budd sliding laws in the literature. The impossibility of comparing the various dimensional60 sliding coefficients Cbfor different exponents pp, qqbecomes immediately evident. They do not even have61 a common unit, and the respective numerical value tells nothing about the actual strength of basal sliding.62 In the second case, pp, qq“p1,2q, the numerical value of Cbis greater than 109; however, the small63 dimensionless value means that it produces only very little basal sliding ( ˜ Cb«0.04). By contrast, in the64 third case, pp, qq“p3,0q, the numerical value of Cbis merely 10´12; however, it corresponds to pronounced65 basal sliding ( ˜ Cb“10). The dimensionless sliding coefficients ˜ Cbgive a much better idea about what the66 respective value means physically, and allow comparing values across different sliding laws.67 To further strengthen our point, suppose that we wish to test a sliding law with a new set of exponents,68 for instance pp, qq “ p3,1.5q. Working with the dimensional sliding coefficient Cb, we would not have any69 idea which order of magnitude may be suited for its numerical value, and which range of values mean70 Page 6 of 15 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 6 strong or weak sliding. However, if the dimensionless sliding coefficient ˜ Cbis used, we can immediately71 start with an initial guess ˜ Cb“1and, from there on, refine the sliding law by, e.g., tuning to observed72 flow speeds. According to the scaling (7), (8)4, the dimensional equivalent of ˜ Cb“1would be Cb“73 10´2.5m a´1Pa´1.5“3.162 ˆ10´3m a´1Pa´1.5.74 We have only discussed cases with a constant sliding parameter; however, the non-dimensionalization75 method is of course not limited to this. It can also be applied to a spatially variable sliding coefficient, which76 may arise from an inversion procedure (e.g., Morlighem and others, 2013). Alternative sliding laws, such as77 the Coulomb-limited rules discussed by Cornford and others (2020), allow similar non-dimensionalization,78 although we refrain from working out the details here.79 3 FLOW LAWS80 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) 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 (“Einstein summation convention”: summation over the twice-appearing81 index iimplied, thus tii is the trace of the stress tensor), and δij is the Kronecker symbol, in other words,82 the unit tensor in index notation.83 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), Page 7 of 15 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 7 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, the84 deformability of ice differs markedly from that behaviour, and the value is frequently chosen as n“3,85 or within the range from 1.5 to 4.2 (Cuffey and Paterson, 2010) (while recent evidence from laboratory86 experiments actually supports n“1for temperate ice; Schohn and others, 2025).87 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 and88 Blatter, 2009).89 Similar to the procedure in Sect. 2, we now introduce scales (typical values) for the relevant quantities, considered suitable for areas where rather large deformations take place: Typical deviatoric stress rτs “ 105Pa , typical strain rate rds “ 2.5ˆ10´2a´1, (20) where the scale rτsis deemed appropriate for both tD ij and τe. Using Eq. (18) entails the choice for the scale of the rate factor: rAs “ rds rτsn.(21) We introduce the dimensionless quantities tD ij “ rτs˜ tD ij , τe“ rτs˜τe, dij “ rds˜ dij , A “ rAs˜ A , (22) Page 8 of 15 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 14 Getraer B and Morlighem M (2025) Increasing the Glen-Nye power-law exponent accelerates ice-loss projections198 for the Amundsen Sea Embayment, West Antarctica. Geophysical Research Letters,52(7), e2024GL112516 (doi:199 10.1029/2024GL112516)200 Glen JW (1955) The creep of polycrystalline ice. Proceedings of the Royal Society A,228(1175), 519–538 (doi:201 10.1098/rspa.1955.0066)202 Goldberg DN (2011) A variationally derived, depth-integrated approximation to a higher-order glaciological flow203 model. Journal of Glaciology,57(201), 157–170 (doi: 10.3189/002214311795306763)204 Grandadam F (2024) Implementation of the Depth Integrated Viscosity Approximation in SICOPOLIS. Internship205 Report, Claude Bernard University Lyon 1, France, and Hokkaido University, Sapporo, Japan (doi: 10.5281/zen-206 odo.14732938)207 Greve R (2005) Relation of measured basal temperatures and the spatial distribution of the geothermal heat flux for208 the Greenland ice sheet. Annals of Glaciology,42, 424–432 (doi: 10.3189/172756405781812510)209 Greve R and Blatter H (2009) Dynamics of Ice Sheets and Glaciers. Springer, Berlin, Germany etc., ISBN 978-3-210 642-03414-5 (doi: 10.1007/978-3-642-03415-2)211 Greve R and Blatter H (2016) Comparison of thermodynamics solvers in the polythermal ice sheet model SICOPOLIS.212 Polar Science,10(1), 11–23 (doi: 10.1016/j.polar.2015.12.004)213 Hindmarsh RCA (2009) Consistent generation of ice-streams via thermo-viscous instabilities modulated by membrane214 stresses. Geophysical Research Letters,36(6), L06502 (doi: 10.1029/2008GL036877)215 Hindmarsh RCA and Le Meur E (2001) Dynamical processes involved in the retreat of marine ice sheets. Journal of216 Glaciology,47(157), 271–282 (doi: 10.3189/172756501781832269)217 Hooke RL (2005) Principles of Glacier Mechanics. Cambridge University Press, Cambridge, UK and New York, NY,218 USA, 2nd edition219 Hutter K and Jöhnk K (2004) Continuum Methods of Physical Modeling. Springer, Berlin, Germany etc. (doi:220 10.1007/978-3-662-06402-3)221 Lipscomb WH, Price SF, Hoffman MJ, Leguy GR, Bennett AR, Bradley SL, Evans KJ, Fyke JG, Kennedy JH,222 Perego M, Ranken DM, Sacks WJ, Salinger AG, Vargo LJ and Worley PH (2019) Description and evaluation of the223 Community Ice Sheet Model (CISM) v2.1. Geoscientific Model Development,12(1), 387–424 (doi: 10.5194/gmd-224 12-387-2019)225 Page 15 of 15 Cambridge University Press Journal of Glaciology For Peer Review Greve: Non-dimensional sliding and flow laws 15 Millstein JD, Minchew BM and Pegler SS (2022) Ice viscosity is more sensitive to stress than commonly assumed.226 Communications Earth & Environment,3, 57 (doi: 10.1038/s43247-022-00385-x)227 Morlighem M, Seroussi H, Larour E and Rignot E (2013) Inversion of basal friction in Antarctica using exact and228 incomplete adjoints of a higher-order model. Journal of Geophysical Research: Earth Surface,118(3), 1746–1753229 (doi: 10.1002/jgrf.20125)230 Nye JF (1957) The distribution of stress and velocity in glaciers and ice sheets. Proceedings of the Royal Society A,231 239(1216), 113–133 (doi: 10.1098/rspa.1957.0026)232 Payne AJ and Dongelmans PW (1997) Self-organization in the thermomechanical flow of ice sheets. Journal of233 Geophysical Research: Solid Earth,102(B6), 12219–12233 (doi: 10.1029/97JB00513)234 Payne AJ, Huybrechts P, Abe-Ouchi A, Calov R, Fastook JL, Greve R, Marshall SJ, Marsiat I, Ritz C, Tarasov L235 and Thomassen MPA (2000) Results from the EISMINT model intercomparison: the effects of thermomechanical236 coupling. Journal of Glaciology,46(153), 227–238 (doi: 10.3189/172756500781832891)237 Rückamp M, Greve R and Humbert A (2019) Comparative simulations of the evolution of the Greenland ice sheet238 under simplified Paris Agreement scenarios with the models SICOPOLIS and ISSM. Polar Science,21, 14–25 (doi:239 10.1016/j.polar.2018.12.003)240 Saito F, Abe-Ouchi A, Takahashi K and Blatter H (2016) SeaRISE experiments revisited: potential sources of spread241 in multi-model projections of the Greenland ice sheet. The Cryosphere,10(1), 43–63 (doi: 10.5194/tc-10-43-2016)242 Sayag R and Tziperman E (2008) Spontaneous generation of pure ice streams via flow instability: Role of longi-243 tudinal shear stresses and subglacial till. Journal of Geophysical Research: Solid Earth,113(B5), B05411 (doi:244 10.1029/2007JB005228)245 Schohn CM, Iverson NR, Zoet LK, Fowler JR and Morgan-Witts N (2025) Linear-viscous flow of temperate ice.246 Science,387(6730), 182–185 (doi: 10.1126/science.adp7708)247 SICOPOLIS Authors (2025) SICOPOLIS v25. GitHub, URL https://github.com/sicopolis/sicopolis248 Weertman J (1957) On the sliding of glaciers. Journal of Glaciology,3(21), 33–38 (doi: 10.3189/S0022143000024709)249 Page 16 of 15 Cambridge University Press Journal of Glaciology