Analysis and finite element approximation of a nonlinear stationary stokes problem arising in glaciology

The aim of this paper is to study a nonlinear stationary Stokes problem with mixed boundary conditions that describes the ice velocity and pressure fields of grounded glaciers under Glen's flow law. Using convex analysis arguments, we prove the existence and the uniqueness of a weak solution. A finite element method is applied with approximation spaces that satisfy the inf-sup condition, and a priori error estimates are established by using a quasinorm technique. Several algorithms (including Newton's method) are proposed to solve the nonlinearity of the Stokes problem and are proved to be convergent. Our results are supported by numerical convergence studies