Modelling the three-dimensional, diagnostic anisotropy field of an ice rise

TL;DR

This paper presents three-dimensional simulations of ice crystal orientation fabrics in Derwael Ice Rise, highlighting the importance of anisotropy modeling for interpreting radar data and improving ice flow understanding.

Contribution

It introduces a novel 3D simulation framework for ice fabric anisotropy, incorporating complex flow effects and providing a basis for better interpretation of observational data.

Findings

High horizontal divergence areas show significant fabric variation.

Model parameters greatly influence predicted fabric types.

Framework aids in planning observational campaigns.

Abstract

Polar ice develops anisotropic crystal orientation fabrics under deformation, yet ice is most often modelled as an isotropic fluid. We present three-dimensional simulations of the crystal orientation fabric of Derwael Ice Rise including the surrounding ice shelf using a crystal orientation tensor evolution equation corresponding to a fixed velocity field. We use a semi-Lagrangian numerical method that constrains the degree of crystal orientation evolution to solve the equations in complex flow areas. We perform four simulations based on previous studies, altering the rate of evolution of the crystal anisotropy and its dependence on a combination of the strain rate and deviatoric stress tensors. We provide a framework for comparison with radar observations of the anisotropy field, outlining areas where the assumption of one vertical eigenvector may not hold and provide resulting errors in measured eigenvalues. We recognise the areas of high horizontal divergence at the ends of the flow divide as important areas to make comparisons with observations. Here, poorly constrained model parameters result in the largest difference in fabric type. These results are important in the planning of future campaigns for gathering data to constrain model parameters and as a link between observations and computationally-efficient, simplified models of anisotropy.