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
arXiv:1410.1221 · doi:10.1016/j.jcp.2015.04.047
Abstract
The majority of research on efficient and scalable algorithms in computational science and engineering has focused on the forward problem: given parameter inputs, solve the governing equations to determine output quantities of interest. In contrast, here we consider the broader question: given a (large-scale) model containing uncertain parameters, (possibly) noisy observational data, and a prediction quantity of interest, how do we construct efficient and scalable algorithms to (1) infer the model parameters from the data (the deterministic inverse problem), (2) quantify the uncertainty in the inferred parameters (the Bayesian inference problem), and (3) propagate the resulting uncertain parameters through the model to issue predictions with quantified uncertainties (the forward uncertainty propagation problem)? We present efficient and scalable algorithms for this end-to-end, data-to-prediction process under the Gaussian approximation and in the context of modeling the flow of the Antarctic ice sheet and its effect on sea level. The ice is modeled as a viscous, incompressible, creeping, shear-thinning fluid. The observational data come from InSAR satellite measurements of surface ice flow velocity, and the uncertain parameter field to be inferred is the basal sliding parameter. The prediction quantity of interest is the present-day ice mass flux from the Antarctic continent to the ocean. We show that the work required for executing this data-to-prediction process is independent of the state dimension, parameter dimension, data dimension, and number of processor cores. The key to achieving this dimension independence is to exploit the fact that the observational data typically provide only sparse information on model parameters. This property can be exploited to construct a low rank approximation of the linearized parameter-to-observable map.
References in corpus (1)
Cited by in corpus (28)
- Optimal design of acoustic metamaterial cloaks under uncertainty
- Goal-Oriented Optimal Design of Experiments for Large-Scale Bayesian Linear Inverse Problems
- Solution of nonlinear Stokes equations discretized by high-order finite elements on nonconforming and anisotropic meshes, with application to ice sheet dynamics
- Taylor approximation and variance reduction for PDE-constrained optimal control under uncertainty
- Hessian-based adaptive sparse quadrature for infinite-dimensional Bayesian inverse problems
- Weighted BFBT Preconditioner for Stokes Flow Problems with Highly Heterogeneous Viscosity
- Bayesian Poroelastic Aquifer Characterization from InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate
- Residual-based error correction for neural operator accelerated infinite-dimensional Bayesian inverse problems
- Bayesian inference of heterogeneous epidemic models: Application to COVID-19 spread accounting for long-term care facilities
- Variational Inference at Glacier Scale
- Bayesian Poroelastic Aquifer Characterization from InSAR Surface Deformation Data Part II: Quantifying the Uncertainty
- Randomized Truncated SVD Levenberg-Marquardt Approach to Geothermal Natural State and History Matching
- Stochastic Learning Approach to Binary Optimization for Optimal Design of Experiments
- A fast and scalable computational framework for large-scale and high-dimensional Bayesian optimal experimental design
- Bayesian Inference with Generative Adversarial Network Priors
- Iterative solution of the Lippmann-Schwinger equation in strongly scattering acoustic media by randomized construction of preconditioners
- Variational Bayes' method for functions with applications to some inverse problems
- A Robust Hierarchical Solver for Ill-conditioned Systems with Applications to Ice Sheet Modeling
- A Data-Scalable Randomized Misfit Approach for Solving Large-Scale PDE-Constrained Inverse Problems
- Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initializaiton
- Real-time Bayesian inference at extreme scale: A digital twin for tsunami early warning applied to the Cascadia subduction zone
- Variance-based sensitivity of Bayesian inverse problems to the prior distribution
- Real-time aerodynamic load estimation for hypersonics via strain-based inverse maps
- Taylor approximation for chance constrained optimization problems governed by partial differential equations with high-dimensional random parameters
- Goal-Oriented Real-Time Bayesian Inference for Linear Autonomous Dynamical Systems With Application to Digital Twins for Tsunami Early Warning
- A stochastic subspace approach to gradient-free optimization in high dimensions
- Randomized Strong Recursive Skeletonization: Simultaneous Compression and LU Factorization of Hierarchical Matrices using Matrix-Vector Products
- Stability of Gibbs Posteriors from the Wasserstein Loss for Bayesian Full Waveform Inversion