Variational Inference at Glacier Scale

• Main Author: Brinkerhoff, Douglas J.
• Format: Article in Journal/Newspaper
• Language: unknown
• Published: arXiv 2021
• Subjects: Computational Physics physics.comp-ph; Machine Learning cs.LG; FOS Physical sciences; FOS Computer and information sciences; Greenland; glacier; Ice Sheet
• Online Access: https://dx.doi.org/10.48550/arxiv.2108.07263; https://arxiv.org/abs/2108.07263

We characterize the complete joint posterior distribution over spatially-varying basal traction and and ice softness parameters of an ice sheet model from observations of surface speed by using stochastic variational inference combined with natural gradient descent to find an approximating variational distribution. By placing a Gaussian process prior over the parameters and casting the problem in terms of eigenfunctions of a kernel, we gain substantial control over prior assumptions on parameter smoothness and length scale, while also rendering the inference tractable. In a synthetic example, we find that this method recovers known parameters and accounts for mutual indeterminacy, both of which can influence observed surface speed. In an application to Helheim Glacier in Southeast Greenland, we show that our method scales to glacier-sized problems. We find that posterior uncertainty in regions of slow flow is high regardless of the choice of observational noise model.