Analysis and numerical simulations of a landfast ice model

TL;DR

This paper provides a rigorous analytical and numerical investigation of a landfast ice model, revealing its well-posedness, equilibrium states, and differences from classical models, with implications for understanding Arctic ice behavior.

Contribution

The study offers the first comprehensive analytical and numerical analysis of a landfast ice model, including well-posedness and long-term behavior insights.

Findings

Existence of local and global strong solutions under certain conditions.

Numerical simulations show stationary equilibrium states with vanishing velocity.

Differences between landfast and classical sea-ice models are demonstrated.

Abstract

In this manuscript, we consider a common modeling framework for Arctic landfast ice based on the work of Lemieux et al. [27], which is designed for use in large-scale climate models. This approach extends the classical viscous-plastic sea-ice model introduced by Hibler [18], which remains the most used model for simulating large-scale sea-ice dynamics in climate science. In particular, landfast ice refers to sea-ice that is attached to the coastline or grounded and therefore exhibits nearly vanishing motion. We present a rigorous analytical and numerical study of this landfast ice model. The main analytical contributions are the local strong well-posedness, the global strong well-posedness in the absence of external forces and for initial data close to constant equilibrium solutions, and the existence of time-periodic solutions. Complementing the analysis, we perform numerical simulations that illustrate key qualitative differences between landfast ice and classical viscous-plastic sea-ice models. In particular, the simulations reveal the formation of stationary equilibrium states characterized by vanishing ice velocity. These observations are consistent with the global-in-time existence result close to equilibria established in Theorem 4.1 as well as the time-periodic result in Theorem 5.2. The combined analytical and numerical results provide new insight into the structure, stability, and long-term behavior of landfast ice dynamics.