Ice Sheet Modeling: Accuracy of First-Order Stokes Model with Basal Sliding

• Main Author: Jonsson, Eskil
• Format: Bachelor Thesis
• Language: English
• Published: Uppsala universitet, Institutionen för geovetenskaper 2018
• Subjects: Ice sheet modelling; First-Order Stokes model; basal sliding; Stokes equations; Physical Geography; Naturgeografi; Ice Sheet; Ice Shelf
• Online Access: http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-360245

Some climate models are still lacking features such as dynamical modelling of ice sheets due to their computational cost which results in poor accuracy and estimates of e.g. sea level rise. The need for low-cost high-order models initiated the development of the First-Order Stokes (or Blatter-Pattyn) model which retains much of the accuracy of the full-Stokes model but is also cost-effective. This model has proven accurate for ice sheets and glaciers with frozen bedrocks, or no-slip basal boundary conditions. However, experimental evidence seems to be lacking regarding its accuracy under sliding, or stress-free, bedrock conditions (ice-shelf conditions). Hence, it became of interest to investigate this. Numerical experiments were set up by formulating the first-order Stokes equations as a variational finite element problem, followed by implementing them using the open-source FEniCS framework. Two types of geometries were used with both no-slip and slip basal boundary conditions. Specifically, experiments B and D from the Ice Sheet Model Intercomparison Project for Higher-Order ice sheet Models (ISMIP-HOM) were used to benchmark the model. Local model errors were investigated and a convergence analysis was performed for both experiments. The results yielded an inherent model error of about 0.06% for ISMIP-HOM B and 0.006% for ISMIPHOM D, mostly relating to the different types of geometries used. Errors in stress-free regions were greater and varied on the order of 1%. This was deemed fairly accurate, and probably enough justification to replace models such as the Shallow Shelf Approximation with the First-Order Stokes model in some regions. However, more rigorous tests with real-world geometries may be warranted. Also noteworthy were inconsistent results in the vertical velocity under slippery conditions (ISMIPHOM D) which could either be due to coding errors or an inherent problem with the decoupling of the horizontal and vertical velocities of the First-Order Stokes model. This should be further investigated. ...