On the dynamics of shallow ice sheets. Modelling and analysis

Paolo Piersanti,Roger Temam
2024-07-06
Abstract:We formulate a model describing the evolution of thickness in a shallow ice sheet lying over a lithosphere. The model is thus governed by a set of variational inequalities that involve nonlinearities in the time derivative and in the elliptic term. Solutions are constructed via a finite difference scheme in time.
Analysis of PDEs,Mathematical Physics
What problem does this paper attempt to address?
The paper aims to address the mathematical modeling and analysis of the evolution of shallow ice sheet thickness. Specifically, the authors propose a model that describes the evolution of shallow ice sheet thickness and treat it as a time-dependent obstacle problem. In this model, the ice sheet thickness is constrained by non-negativity, meaning that the ice sheet thickness is always greater than or equal to the bedrock height. The paper describes this problem by introducing a variational inequality and proves that the model has a solution under certain assumptions. Additionally, the paper discusses how to prove the existence of the solution through penalty methods and semi-discrete schemes. This study fills a gap in the existing literature regarding the research on time-dependent obstacle problems in two-dimensional spatial domains.