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

Abstract

Abstract. 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.