Comparison of ice dynamics using full-Stokes and Blatter–Pattyn approximation: application to the Northeast Greenland Ice Stream

Full-Stokes (FS) ice sheet models provide the most sophisticated formulation of ice sheet flow. However, their applicability is often limited due to the high computational demand and numerical challenges. To balance computational demand and accuracy, the so-called Blatter–Pattyn (BP) stress regime i...

Full description

Bibliographic Details
Published in:The Cryosphere
Main Authors: M. Rückamp, T. Kleiner, A. Humbert
Format: Article in Journal/Newspaper
Language:English
Published: Copernicus Publications 2022
Subjects:
geo
Online Access:https://doi.org/10.5194/tc-16-1675-2022
https://tc.copernicus.org/articles/16/1675/2022/tc-16-1675-2022.pdf
https://doaj.org/article/d9efccf55dff4833a40f104f25995cd1
Description
Summary:Full-Stokes (FS) ice sheet models provide the most sophisticated formulation of ice sheet flow. However, their applicability is often limited due to the high computational demand and numerical challenges. To balance computational demand and accuracy, the so-called Blatter–Pattyn (BP) stress regime is frequently used. Here, we explore the dynamic consequences of using simplified approaches by solving FS and the BP stress regime applied to the Northeast Greenland Ice Stream. To ensure a consistent comparison, we use one single ice sheet model to run the simulations under identical numerical conditions. A sensitivity study to the horizontal grid resolution (from 12.8 to a resolution of 0.1 km) reveals that velocity differences between the FS and BP solution emerge below ∼ 1 km horizontal resolution and continuously increase with resolution. Over the majority of the modelling domain both models reveal similar surface velocity patterns. At the grounding line of the 79∘ North Glacier the simulations show considerable differences whereby the BP model overestimates ice discharge of up to 50 % compared to FS. A sensitivity study to the friction type reveals that differences are stronger for a power-law friction than a linear friction law. Model differences are attributed to topographic variability and the basal drag, in which neglected stress terms in BP become important.