Abstract:Highlights • An efficient stochastic realization approach is proposed for sampling Gaussian stationary random fields. • A rational spectrum is constructed using covariance extension techniques for a rational ARMA model. • Samples of random fields are generated via the ARMA recursion with a low computational cost. • The proposed approach is capable of performing multiscale simulations. Generating large-scale samples of stationary random fields is of great importance in the fields such as geomaterial modeling and uncertainty quantification. Traditional methodologies based on covariance matrix decomposition have the difficulty of being computationally expensive, which is even more serious when the dimension of the random field is large. This paper proposes an efficient stochastic realization approach for sampling Gaussian stationary random fields from a systems and control point of view. Specifically, we take the exponential and squared exponential covariance functions as examples and make a decoupling assumption when there are multiple dimensions. Then a rational spectral density is constructed in each dimension using techniques from covariance extension, and the corresponding autoregressive moving-average (ARMA) model is obtained via spectral factorization. As a result, samples of the random field with a specific covariance function can be generated very efficiently in the space domain by implementing the ARMA recursion using a Gaussian white noise input. Such a procedure is computationally cheap due to the fact that the constructed ARMA model has a low order. Furthermore, the same method is integrated to multiscale simulations where interpolations of the generated samples are achieved when one zooms into finer scales. Both theoretical analysis and simulation results show that our approach performs favorably compared with covariance matrix decomposition methods.

Rational Covariance Functions for Nonstationary Random Fields.

The moduli of non-differentiability for Gaussian random fields with stationary increments

Stationary Random Fields in Space and Time with Rational Spectral Densities.

Generation of Two Correlated Stationary Gaussian Processes

Covariance Matrix Functions of Isotropic Vector Random Fields

Vector Random Fields with Compactly Supported Covariance Matrix Functions

Covariance Matrix Functions of Vector Χ^2 Random Fields in Space and Time

Vector Random Fields with Long-Range Dependence

Time-Varying Isotropic Vector Random Fields on Compact Two-Point Homogeneous Spaces

Hybrid Parametric Classes of Isotropic Covariance Functions for Spatial Random Fields

Stochastic Representations of Isotropic Vector Random Fields on Spheres

Isotropic Covariance Matrix Functions on Compact Two-Point Homogeneous Spaces

Nonstationary covariance functions that model space–time interactions

The Circular Matern Covariance Function and its Link to Markov Random Fields on the Circle

Covariance Matrices for Second-Order Vector Random Fields in Space and Time.

Asymptotic covariances for functionals of weakly stationary random fields

Sampling Gaussian stationary random fields: A stochastic realization approach

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

Isotropic Covariance Matrix Polynomials on Spheres

ℓ1-symmetric vector random fields

Isotropic Covariance Matrix Functions on All Spheres