Variational Inference at Glacier Scale

Abstract

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.

References in corpus (9)

• Variational Inference: A Review for Statisticians
• MCMC using Hamiltonian dynamics
• Weight Uncertainty in Neural Networks
• B-PINNs: Bayesian Physics-Informed Neural Networks for Forward and Inverse PDE Problems with Noisy Data
• Automatic Differentiation Variational Inference
• Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet
• Hilbert Space Methods for Reduced-Rank Gaussian Process Regression
• Approximate Bayesian Model Inversion for PDEs with Heterogeneous and State-Dependent Coefficients
• Bayesian and Variational Bayesian approaches for flows in heterogenous random media