The numerical stability of linear systems arising in kriging, estimation, and simulation of random fields, is studied analytically and numerically. In the state-space formulation of kriging, as developed here, the stability of the kriging system depends on the condition number of the prior, stationary covariance matrix. The same is true for conditional random field generation by the superposition method, which is based on kriging, and the multivariate Gaussian method, which requires factoring a covariance matrix. A large condition number corresponds to an ill-conditioned, numerically unstable system. In the case of stationary covariance matrices and uniform grids, as occurs in kriging of uniformly sampled data, the degree of ill-conditioning generally increases indefinitely with sampling density and, to a limit, with domain size. The precise behavior is, however, highly sensitive to the underlying covariance model. Detailed analytical and numerical results are given for five one-dimensional covariance models: (1) hole-exponential, (2) exponential, (3) linear-exponential, (4) hole-Gaussian, and (5) Gaussian. This list reflects an approximate ranking of the models, from “best” to “worst” conditioned. The methods developed in this work can be used to analyze other covariance models. Examples of such representative analyses, conducted in this work, include the spherical and periodic hole-effect (hole-sinusoidal) covariance models. The effect of small-scale variability (nugget) is addressed and extensions to irregular sampling schemes and higher dimensional spaces are discussed.
[1]
M. Ekstrom.
A spectral characterization of the ill-conditioning in numerical deconvolution
,
1973,
IEEE Transactions on Audio and Electroacoustics.
[2]
J. P. Delhomme,et al.
Spatial variability and uncertainty in groundwater flow parameters: A geostatistical approach
,
1979
.
[3]
A. Mantoglou,et al.
The Turning Bands Method for simulation of random fields using line generation by a spectral method
,
1982
.
[4]
Eric F. Wood,et al.
On the relationship between kriging and state estimation
,
1982
.
[5]
G. Dagan.
Stochastic modeling of groundwater flow by unconditional and conditional probabilities: 1. Conditional simulation and the direct problem
,
1982
.
[6]
Erik H. Vanmarcke,et al.
Random Fields: Analysis and Synthesis.
,
1985
.
[7]
Gene H. Golub,et al.
Matrix computations
,
1983
.
[8]
G. Christakos.
On the Problem of Permissible Covariance and Variogram Models
,
1984
.
[9]
John P. Lewis,et al.
Generalized stochastic subdivision
,
1987,
TOGS.
[10]
W. Press,et al.
Numerical Recipes: The Art of Scientific Computing
,
1987
.
[11]
R. Ababou,et al.
Implementation of the three‐dimensional turning bands random field generator
,
1989
.
[12]
D. Posa.
Conditioning of the stationary kriging matrices for some well-known covariance models
,
1989
.
[13]
R. J. O'dowd.
Conditioning of coefficient matrices of Ordinary Kriging
,
1991
.