文章基本信息

标题：Stable finite volume element schemes for the shallow-ice approximation
本地全文：下载
作者：ED BUELER
期刊名称：Journal of Glaciology
印刷版ISSN：0022-1430
电子版ISSN：1727-5652
出版年度：2016
卷号：62
期号：232
页码：230-242
DOI：10.1017/jog.2015.3
出版社：Cambridge University Press
摘要：ABSTRACT The isothermal, non-sliding shallow-ice approximation, combined with mass conservation, is a fundamental model for ice-sheet and glacier flow. It determines the ice extent, geometry and velocity by the solution of a free-boundary problem. In this paper, the steady-state form of this problem is solved directly, without time-stepping, thereby demonstrating a fully implicit scheme with no stability restrictions. The classical Mahaffy (1976) finite difference method is first reinterpreted as a ‘finite volume element’ scheme that has both an everywhere-defined approximate thickness function and an approximation of the conservation equation in flux integral form. From this reinterpretation an improved scheme is built by using better quadrature in the integral and upwinding on that part of the flux which is proportional to the bed gradient. The discrete equations are then solved by a parallel Newton scheme which respects the constraint that ice thickness is non-negative. The results show good accuracy on both flat-bed and bedrock-step exact solutions. The method is then applied at high resolution to model the steady-state geometry of the Greenland ice sheet, using only bedrock elevation and present-day surface mass balance as input data.
关键词：glacier modelling; ice; sheet modelling; free boundary; mass conservation; bedrock topography