Ice flow localisation enhanced by composite ice rheology

Ice’s predominantly viscous rheology exhibits a significant temperature and strain-rate dependence, commonly captured as a single deformation mechanism by Glen's flow law. However, Glen’s power-law relationship may fail to capture accurate stress levels at low and elevated strain-rates ultimately leading to velocity over- and under-estimates, respectively. Alternative more complex flow laws such as Goldsby rheology combine various creep mechanisms better accounting for micro-scale observations resulting in enhanced localisation of ice flow at glacier scales and internal sliding.

The challenge in implementing Goldsby rheology arises with the need of computing an accurate partitioning of the total strain-rate among the active creep mechanisms. Some of these mechanisms exhibit grain-size evolution sensitivity potentially impacting the larger scale ice dynamics.

We here present a consistent way to compute the effective viscosity of the ice using Goldsby rheology for temperature and strain-rate ranges relevant to ice flow. We implement a local iteration procedure to ensure accurate implicit partitioning of the total strain-rate among the active creep mechanisms including grain-size evolution. We discuss the composite deformation maps and compare the results against Glen's flow law. We incorporate our implicit rheology solver into an implicit 2D thermo-mechanical ice flow solver to investigate localisation of ice flow over variable topography and in shear margin configurations. We quantify discrepancies  in surface velocity patterns when using Goldsby rheology instead of Glen's flow law.