Surface elevation errors in finite element Stokes models for glacier evolution

内容摘要

The primary data which determine the evolution of glaciation are the bedrock elevation and the surface mass balance. From this data, which we assume is defined over a fixed land region, the glacier's geometry solves a free boundary problem which balances the time derivative of the surface elevation, the surface velocity from the Stokes flow of the ice, and the surface balance rate. This problem can be posed in weak form as a variational inequality over a cone of admissible surface elevation functions, those which are above the bedrock topography. After some preparatory theory for the Stokes problem, we conjecture that the corresponding continuous-space, implicit time-step variational inequality problem is well-posed if the surface kinematical equation is appropriately regularized. This conjecture is supported by physical arguments and numerical evidence. We then prove a general theorem which bounds the numerical error made by finite element approximations of nonlinear-operator variational inequalities in Banach spaces. This bound is a sum of error terms of different types, special to variational inequalities. When it is applied to the implicit time-step glacier problem there are th