Abstract:The distribution of centered Gaussian random fields (GRFs) indexed by compacta such as smooth, bounded Euclidean domains or smooth, compact and orientable manifolds is determined by their covariance operators. We consider centered GRFs given as variational solutions to coloring operator equations driven by spatial white noise, with an elliptic self-adjoint pseudodifferential coloring operator from the Hörmander class. This includes the Matérn class of GRFs as a special case. Using biorthogonal multiresolution analyses on the manifold, we prove that the precision and covariance operators, respectively, may be identified with bi-infinite matrices and finite sections may be diagonally preconditioned rendering the condition number independent of the dimension p of this section. We prove that a tapering strategy by thresholding applied on finite sections of the bi-infinite precision and covariance matrices results in optimally numerically sparse approximations. That is, asymptotically only linearly many nonzero matrix entries are sufficient to approximate the original section of the bi-infinite covariance or precision matrix using this tapering strategy to arbitrary precision. The locations of these nonzero matrix entries can be determined a priori. The tapered covariance or precision matrices may also be optimally diagonally preconditioned. Analysis of the relative size of the entries of the tapered covariance matrices motivates novel, multilevel Monte Carlo (MLMC) oracles for covariance estimation, in sample complexity that scales log-linearly with respect to the number p of parameters. In addition, we propose and analyze novel compressive algorithms for simulating and kriging of GRFs. The complexity (work and memory vs. accuracy) of these three algorithms scales near-optimally in terms of the number of parameters p of the sample-wise approximation of the GRF in Sobolev scales.

Constrained realizations and minimum variance reconstruction of non-Gaussian random fields

Convergence of the Centered Maximum of Log-Correlated Gaussian Fields

Wiener Reconstruction of The Large Scale Structure

Sampling Gaussian stationary random fields: A stochastic realization approach

Multilevel approximation of Gaussian random fields: Covariance compression, estimation, and spatial prediction

A Unifying Framework for Global Gaussianization: Asymptotic Equivalence of Locally Stationary Processes and Bivariate White Noise

Simultaneous Inference for Non-Stationary Random Fields, with Application to Gridded Data Analysis

Fixed-domain asymptotic properties of maximum composite likelihood estimators for max-stable Brown-Resnick random fields

The exact non-Gaussian weak lensing likelihood: A framework to calculate analytic likelihoods for correlation functions on masked Gaussian random fields

Dynamical Gaussian, lognormal, and reverse lognormal maximum‐likelihood ensemble filter

$G$-Gaussian random fields and stochastic quantization under nonlinear expectation

Theoretical and Computational Guarantees of Mean Field Variational Inference for Community Detection

Second Maximum of a Gaussian Random Field and Exact (t-)Spacing test

Adaptive Bayesian estimation using a Gaussian random field with inverse Gamma bandwidth

Nonlinear Realization of the General Covariance Group Revisited

An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach

Construction of non-Gaussian random fields with any given correlation structure

Finite range decompositions of Gaussian fields with applications to level-set percolation

Fast estimation of the look-elsewhere effect using Gaussian random fields

Asymptotic properties of the maximum likelihood and cross validation estimators for transformed Gaussian processes

Simulation of Gaussian random field in a ball