A parallel high-order accurate finite element nonlinear Stokes ice sheet model and benchmark experiments

The numerical modeling of glacier and ice sheet evolution is a subject of growing interest, in part because of the potential for models to inform estimates of global sea level change. This paper focuses on the development of a numerical model that determines the velocity and pressure fields within a...

Full description

Bibliographic Details
Published in:Journal of Geophysical Research: Earth Surface
Main Authors: Leng, Wei, Ju, Lili, Gunzburger, Max, Price, Stephen, Ringler, Todd
Other Authors: United States. Department of Energy., United States. Department of Energy. Office of Biological and Environmental Research.
Format: Article in Journal/Newspaper
Language:English
Published: University of South Carolina 2012
Subjects:
Online Access:https://doi.org/10.1029/2011JF001962
https://digital.library.unt.edu/ark:/67531/metadc838667/
Description
Summary:The numerical modeling of glacier and ice sheet evolution is a subject of growing interest, in part because of the potential for models to inform estimates of global sea level change. This paper focuses on the development of a numerical model that determines the velocity and pressure fields within an ice sheet. Our numerical model features a high-fidelity mathematical model involving the nonlinear Stokes system and combinations of no-sliding and sliding basal boundary conditions, high-order accurate finite element discretizations based on variable resolution grids, and highly scalable parallel solution strategies, all of which contribute to a numerical model that can achieve accurate velocity and pressure approximations in a highly efficient manner. We demonstrate the accuracy and efficiency of our model by analytical solution tests, established ice sheet benchmark experiments, and comparisons with other well-established ice sheet models.