Fast sampling of parameterised Gaussian random fields

Abstract

Gaussian random fields are popular models for spatially varying uncertainties, arising for instance in geotechnical engineering, hydrology or image processing. A Gaussian random field is fully characterised by its mean function and covariance operator. In more complex models these can also be partially unknown. In this case we need to handle a family of Gaussian random fields indexed with hyperparameters. Sampling for a fixed configuration of hyperparameters is already very expensive due to the nonlocal nature of many classical covariance operators. Sampling from multiple configurations increases the total computational cost severely. In this report we employ parameterised Karhunen-Loève expansions for sampling. To reduce the cost we construct a reduced basis surrogate built from snapshots of Karhunen-Loève eigenvectors. In particular, we consider Matérn-type covariance operators with unknown correlation length and standard deviation. We suggest a linearisation of the covariance function and describe the associated online-offline decomposition. In numerical experiments we investigate the approximation error of the reduced eigenpairs. As an application we consider forward uncertainty propagation and Bayesian inversion with an elliptic partial differential equation where the logarithm of the diffusion coefficient is a parameterised Gaussian random field. In the Bayesian inverse problem we employ Markov chain Monte Carlo on the reduced space to generate samples from the posterior measure. All numerical experiments are conducted in 2D physical space, with non-separable covariance operators, and finite element grids with $\sim 10^4$ degrees of freedom.

35 pages

References in corpus (12)

• Numerical solution of fractional elliptic stochastic PDEs with spatial white noise
• The ROMES method for statistical modeling of reduced-order-model error
• Coordinate Transformation and Polynomial Chaos for the Bayesian Inference of a Gaussian Process with Parametrized Prior Covariance Function
• Fast spatial Gaussian process maximum likelihood estimation via skeletonization factorizations
• Survey of multifidelity methods in uncertainty propagation, inference, and optimization
• Weak convergence of Galerkin approximations for fractional elliptic stochastic PDEs with spatial white noise
• Parallel Markov Chain Monte Carlo
• Goal-oriented adaptive surrogate construction for stochastic inversion
• MCMC Methods for Functions: Modifying Old Algorithms to Make Them Faster
• Multilevel ensemble Kalman filtering for spatially extended models
• Low-rank computation of posterior covariance matrices in Bayesian inverse problems
• Linear-Cost Covariance Functions for Gaussian Random Fields

Cited by in corpus (6)

• On the well-posedness of Bayesian inverse problems
• Regularity and convergence analysis in Sobolev and Hölder spaces for generalized Whittle-Matérn fields
• Multilevel adaptive sparse Leja approximations for Bayesian inverse problems
• Nesterov Acceleration for Ensemble Kalman Inversion and Variants
• Gaussian Random Fields as an Abstract Representation of Patient Metadata for Multimodal Medical Image Segmentation
• Hierarchical surrogate-based Approximate Bayesian Computation for an electric motor test bench