Higher-order Numerical Simulation of Elasto-visco-plastic Sea Ice Models on Idealised Benchmarks

In this study, we numerically compare two elasto-visco-plastic sea ice models: the Maxwell Elasto Brittle (MEB) and the Brittle Bingham Maxwell (BBM). We examine the linear kinematic features and overall deformation of sea ice in two idealised scenarios: when the ice is compressed in one direction and when it is affected by a moving cyclone. We also provide a detailed analysis of the different parameters used in the models and their effect on the simulation results.

The models are solved using a higher-order explicit discontinuous Galerkin discretization, implemented in a software package neXtSIM_DG: next generation sea-ice model with DG. Numerical implementation is available at https://github.com/nextsimdg/nextsimdg.