Journal cover Journal topic
The Cryosphere An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

Journal metrics

  • IF value: 4.524 IF 4.524
  • IF 5-year value: 5.558 IF 5-year 5.558
  • CiteScore value: 4.84 CiteScore 4.84
  • SNIP value: 1.425 SNIP 1.425
  • SJR value: 3.034 SJR 3.034
  • IPP value: 4.65 IPP 4.65
  • h5-index value: 52 h5-index 52
  • Scimago H index value: 55 Scimago H index 55
Volume 12, issue 7 | Copyright
The Cryosphere, 12, 2229-2248, 2018
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.

Research article 11 Jul 2018

Research article | 11 Jul 2018

A Bayesian hierarchical model for glacial dynamics based on the shallow ice approximation and its evaluation using analytical solutions

Giri Gopalan1, Birgir Hrafnkelsson1, Guðfinna Aðalgeirsdóttir2, Alexander H. Jarosch2, and Finnur Pálsson2 Giri Gopalan et al.
  • 1Faculty of Physical Sciences, School of Engineering and Natural Sciences, University of Iceland, Reykjavik, Iceland
  • 2Institute of Earth Sciences, University of Iceland, Reykjavik, Iceland

Abstract. Bayesian hierarchical modeling can assist the study of glacial dynamics and ice flow properties. This approach will allow glaciologists to make fully probabilistic predictions for the thickness of a glacier at unobserved spatiotemporal coordinates, and it will also allow for the derivation of posterior probability distributions for key physical parameters such as ice viscosity and basal sliding. The goal of this paper is to develop a proof of concept for a Bayesian hierarchical model constructed, which uses exact analytical solutions for the shallow ice approximation (SIA) introduced by Bueler et al. (2005). A suite of test simulations utilizing these exact solutions suggests that this approach is able to adequately model numerical errors and produce useful physical parameter posterior distributions and predictions. A byproduct of the development of the Bayesian hierarchical model is the derivation of a novel finite difference method for solving the SIA partial differential equation (PDE). An additional novelty of this work is the correction of numerical errors induced through a numerical solution using a statistical model. This error-correcting process models numerical errors that accumulate forward in time and spatial variation of numerical errors between the dome, interior, and margin of a glacier.

Publications Copernicus
Short summary
Geophysical systems can often contain scientific parameters whose values are uncertain, complex underlying dynamics, and field measurements with errors. These components are naturally modeled together within what is known as a Bayesian hierarchical model (BHM). This paper constructs such a model for shallow glaciers based on an approximation of the underlying dynamics. The evaluation of this model is aided by the use of exact analytical solutions from the literature.
Geophysical systems can often contain scientific parameters whose values are uncertain, complex...