Implementing an empirical scalar tertiary anisotropic rheology (ESTAR) into large-scale ice sheet models

Abstract

Abstract. The microstructural evolution that occurs in polycrystalline ice during deformation leads to the development of anisotropic rheological properties that are not adequately described by the most common, isotropic, ice flow relation used in large-scale ice sheet models -the Glen flow relation. We present a preliminary assessment of the implementation in the Ice Sheet System Model (ISSM) of a computationally-efficient, empirical, scalar, tertiary, anisotropic rheology (ESTAR). The effect of this anisotropic rheology on ice flow dynamics is investigated by comparing idealised simulations using ESTAR with 5 those using the isotropic Glen flow relation, where the latter includes a flow enhancement factor. For an idealised embayed ice shelf, the Glen flow relation overestimates velocities by up to 17% when using an enhancement factor equivalent to the maximum value prescribed by ESTAR. Importantly, no single Glen enhancement factor can accurately capture the spatial variations in flow over the ice shelf. For flow-line studies of idealised grounded flow over a bumpy topography or a sticky base -both scenarios dominated at depth by bed-parallel shear -the differences between simulated velocities using ESTAR and the 10 Glen flow relation vary according to the value of the enhancement factor used to calibrate the Glen flow relation. These results demonstrate the importance of describing the anisotropic rheology of ice in a physically realistic manner, and have implications for simulations of ice sheet evolution used to reconstruct paleo-ice sheet extent and predict future ice sheet contributions to sea level.