Fact-checked by Grok 2 weeks ago

Kriging

Kriging is a geostatistical method of spatial interpolation that estimates values of a spatial random field at unsampled locations using a weighted average of known observations, incorporating the spatial autocorrelation structure to provide optimal unbiased predictions along with associated uncertainty measures.^[1] Originating in the field of mining engineering, it was pioneered in the 1950s by South African engineer Danie G. Krige through empirical techniques for estimating gold ore grades from borehole samples.^[2] The approach was formalized in 1963 by French mathematician Georges Matheron, who developed the theoretical framework of geostatistics and named the method "kriging" in Krige's honor. At its core, kriging models the spatial dependence of data through the variogram, a function that quantifies how dissimilarity between observations increases with distance, enabling the calculation of weights that minimize prediction error.^[2] It assumes stationarity in the statistical properties of the field and serves as an exact interpolator, meaning predictions at observed locations match the data exactly.^[1] Common variants include ordinary kriging, which estimates a constant unknown mean; simple kriging, assuming a known mean; and universal kriging, which incorporates deterministic trends like spatial coordinates.^[1] These extensions allow kriging to handle non-stationarities and provide the best linear unbiased estimator under Gaussian assumptions, often viewed as a form of Gaussian process regression.^[3] Kriging has become a foundational tool in geostatistics, applied across diverse domains such as environmental monitoring for mapping soil contaminants, hydrology for groundwater modeling, and petroleum engineering for reservoir characterization.^[1] Its ability to generate probabilistic estimates distinguishes it from deterministic interpolators like inverse distance weighting, supporting decision-making in resource estimation and risk assessment.^[2] Since its inception, advancements in computational tools have expanded its use to large-scale datasets, including remote sensing and climate modeling, while maintaining its emphasis on rigorous uncertainty quantification.^[1]

Introduction

Definition and Purpose

Kriging is a geostatistical interpolation technique that estimates values of a spatial random field at unsampled locations using observed data points, under the assumption of spatial autocorrelation among nearby observations.^[3] This method originated in the field of mining geology and is now widely applied in environmental science, hydrology, and resource estimation to predict continuous spatial phenomena such as pollutant concentrations or soil properties.^[2] By modeling the covariance structure of the data, kriging produces predictions that account for the inherent variability and dependence in spatial datasets.^[3] The primary purpose of kriging is to deliver unbiased estimates with the minimum possible variance, making it the best linear unbiased estimator (BLUE) for spatial prediction.^[3] Unlike simpler deterministic methods like inverse distance weighting (IDW), which assign weights based solely on Euclidean distances and often oversmooth data, kriging incorporates a probabilistic model of spatial continuity to yield more precise and reliable interpolations, particularly in heterogeneous environments.^[4] This optimality ensures that predictions are not systematically biased and have the lowest prediction error among linear combinations of the observed data.^[3] In essence, kriging functions as a tailored form of Gaussian process regression for geospatial applications, where the spatial field is treated as a realization of a multivariate Gaussian distribution conditioned on the observations.^[5] A practical example is in mineral resource evaluation, where kriging predicts ore grades at untested sites within a mine based on limited drill core samples, enabling informed decisions on extraction feasibility and reserve volumes.^[6] The method models spatial dependence through tools like the variogram to quantify how similarity decreases with distance.^[3]

Historical Development

The origins of kriging trace back to the work of South African mining engineer Danie G. Krige, who in 1951 developed an empirical method for estimating gold ore grades in Witwatersrand mines using weighted averages of nearby drill hole samples to account for spatial variability.^[7] This approach, detailed in Krige's MSc thesis and subsequent publication, addressed practical challenges in mine valuation by improving the accuracy of reserve estimates over traditional methods like polygonal estimation.^[8] Krige's technique was initially applied routinely in South African gold mines during the early 1950s, marking the first systematic use of spatial interpolation in mineral resource evaluation.^[9] The theoretical formalization of kriging occurred in 1963 through the efforts of French mathematician Georges Matheron at the École des Mines de Paris, who built upon Krige's empirical results to establish a rigorous geostatistical framework.^[10] Matheron coined the term "kriging" in 1962 as a tribute to Krige, integrating concepts like the variogram—introduced to quantify spatial continuity and dependence in data—to derive unbiased linear predictors for unsampled locations.^[11] This development transformed Krige's practical tool into a statistically grounded method, emphasizing best linear unbiased estimation under assumptions of spatial stationarity. Kriging evolved rapidly from its mining roots in the 1950s to a foundational geostatistical technique by the 1970s, with Matheron's seminal two-volume "Traité de Géostatistique Appliquée" (1962–1963) providing the comprehensive theoretical basis and applications that popularized it beyond ore evaluation.^[12] By the 1980s, kriging had gained traction in environmental sciences for tasks such as mapping pollutant distributions and groundwater contamination, further facilitated in the 1990s by accessible software like the U.S. Environmental Protection Agency's Geo-EAS package, which implemented kriging routines for spatial analysis of environmental data.^[13] This period solidified kriging's role as a versatile interpolation method across disciplines, bridging empirical mining practices with broader statistical applications.^[14]

Mathematical Foundations

Spatial Dependence and Stationarity

Spatial dependence is a fundamental concept in geostatistics, positing that values of a spatial process at nearby locations exhibit greater similarity than those at distant locations. This principle, known as Tobler's First Law of Geography, states that "everything is related to everything else, but near things are more related than distant things," providing the theoretical foundation for interpolation methods like kriging, where predictive weights are derived from spatial proximity and correlation structures. Kriging relies on this dependence to model spatial autocorrelation, assuming that the influence of observed data points diminishes with increasing distance, thereby enabling unbiased predictions at unsampled locations. In practice, this dependence is quantified through tools like the variogram, which measures dissimilarity as a function of separation distance, though detailed modeling is addressed elsewhere. Stationarity assumptions underpin the validity of these models by ensuring that statistical properties of the spatial process remain consistent across the domain. Strict stationarity requires that the joint probability distribution of the process is invariant under spatial translations, implying uniformity in all moments and higher-order dependencies. Second-order stationarity, a weaker and more commonly invoked condition in geostatistics, assumes a constant mean throughout the domain and a covariance function that depends solely on the lag vector between points, allowing the process variance to be well-defined and separable from location. Intrinsic stationarity relaxes these further by focusing on the stationarity of increments rather than the process itself, where the variance of differences between values at points separated by a fixed lag is constant, facilitating the use of variograms even when means or variances vary spatially. This form is particularly useful for processes exhibiting local homogeneity in fluctuations but global trends.^[15] For more complex non-stationary scenarios, intrinsic random functions of order k (IRF-k) generalize the framework by assuming that the k-th order differences of the process are intrinsically stationary, accommodating polynomial drifts or trends of order up to k-1. Introduced by Matheron, these functions extend geostatistical inference to datasets with underlying smooth variations, such as geological formations, by stabilizing higher-order increments.^[15] Examples illustrate these concepts in real-world applications. Climate variables like daily temperature in a flat agricultural region often satisfy second-order stationarity over moderate scales, with constant means and lag-dependent covariances reflecting atmospheric mixing. In contrast, elevation data across varied topography typically exhibit non-stationarity due to systematic trends from underlying geology, necessitating IRF approaches or trend removal to model residual fluctuations effectively.^[16]

Variogram and Covariance Functions

The semivariogram, denoted as \gamma(\mathbf{h}), quantifies the average dissimilarity between values of a spatial random field Z(\mathbf{x}) separated by a lag vector \mathbf{h}, and is defined as \gamma(\mathbf{h}) = \frac{1}{2} \mathbb{E} \left[ (Z(\mathbf{x}) - Z(\mathbf{x} + \mathbf{h}))^2 \right].^[17] This measure arises under the assumption of intrinsic stationarity, where the first two moments of the increments are translation-invariant.^[17] The empirical semivariogram \hat{\gamma}(\mathbf{h}) is estimated from observed data pairs \{Z(\mathbf{x}_i), Z(\mathbf{x}_i + \mathbf{h})\} using the method-of-moments estimator \hat{\gamma}(\mathbf{h}) = \frac{1}{2 |N(\mathbf{h})|} \sum_{N(\mathbf{h})} [Z(\mathbf{x}_i) - Z(\mathbf{x}_i + \mathbf{h})]^2, where N(\mathbf{h}) is the set of pairs separated by \mathbf{h}. Estimation considers potential anisotropy by directional binning of lags, allowing separate models for different orientations if spatial dependence varies with direction. The nugget effect, appearing as \hat{\gamma}(0^+), captures microscale variability or measurement error unresolved at the sampling scale. Theoretical semivariogram models are fitted to the empirical values to ensure validity for kriging. Common isotropic models include the spherical \gamma(\mathbf{h}) = \begin{cases} c_0 + c \left( \frac{3}{2} \frac{|\mathbf{h}|}{a} - \frac{1}{2} \left( \frac{|\mathbf{h}|}{a} \right)^3 \right) & |\mathbf{h}| \leq a \\ c_0 + c & |\mathbf{h}| > a \end{cases}, exponential \gamma(\mathbf{h}) = c_0 + c \left(1 - e^{-|\mathbf{h}|/a}\right), and Gaussian \gamma(\mathbf{h}) = c_0 + c \left(1 - e^{-(|\mathbf{h}|/a)^2}\right), where c_0 is the nugget, c the sill, and a the range. Fitting proceeds via weighted least squares, minimizing \sum_w [\hat{\gamma}(\mathbf{h}_k) - \gamma(\mathbf{h}_k; \boldsymbol{\theta})]^2 with weights accounting for estimation variance, or maximum likelihood, maximizing the Gaussian log-likelihood under the model parameters \boldsymbol{\theta}.^[18] Under second-order stationarity, the covariance function C(\mathbf{h}) = \mathrm{Cov}(Z(\mathbf{x}), Z(\mathbf{x} + \mathbf{h})) relates to the semivariogram by C(\mathbf{h}) = C(0) - \gamma(\mathbf{h}), where C(0) = \sigma^2 is the variance (sill). Valid covariance functions must be positive definite, ensuring non-negative variances in any linear combination of the field, a property verified via Bochner's theorem for continuous models or spectral analysis. Cross-validation assesses model adequacy by omitting each data point, predicting it via kriging, and evaluating errors such as mean error (near zero for unbiasedness) and mean squared error (minimized for accuracy). This technique guides selection among candidate models by comparing prediction performance across the dataset.

The Kriging Predictor

General Formulation

Kriging provides a method for predicting the value of a spatial random process Z(\mathbf{x}) at an unsampled location \mathbf{x}_0 based on observed values at known locations \mathbf{x}_1, \dots, \mathbf{x}_n. The general kriging predictor is formulated as a linear combination of the observations:

\hat{Z}(\mathbf{x}_0) = \sum_{i=1}^n \lambda_i Z(\mathbf{x}_i),

where the weights \lambda_i are chosen to minimize the prediction error variance while satisfying the unbiasedness condition \sum_{i=1}^n \lambda_i = 1.^[19] This condition ensures that the expected value of the predictor equals the true value at \mathbf{x}_0, i.e., E[\hat{Z}(\mathbf{x}_0) - Z(\mathbf{x}_0)] = 0, assuming the process has a constant mean.^[19] To determine the optimal weights, the kriging system solves a constrained optimization problem that minimizes the prediction variance subject to the unbiasedness constraint. For ordinary kriging, this leads to the augmented linear system

\begin{pmatrix} \boldsymbol{\Gamma} & \mathbf{1} \\ \mathbf{1}^T & 0 \end{pmatrix} \begin{pmatrix} \boldsymbol{\lambda} \\ \mu \end{pmatrix} = \begin{pmatrix} \boldsymbol{\gamma} \\ 1 \end{pmatrix},

where \boldsymbol{\Gamma} is the n \times n variogram matrix with entries \Gamma_{ij} = \gamma(\mathbf{x}_i - \mathbf{x}_j), \boldsymbol{\gamma} is the vector of variograms between \mathbf{x}_0 and the observation points with \gamma_k = \gamma(\mathbf{x}_0 - \mathbf{x}_k), \gamma(\cdot) denotes the semivariogram function, \mathbf{1} is a vector of ones, and \mu is the Lagrange multiplier.^[19] The variogram matrix \boldsymbol{\Gamma} captures the spatial dependence structure derived from the second-order properties of the process.^[19] The associated prediction variance, or kriging variance, quantifies the uncertainty in the estimate and is given by

\sigma^2(\mathbf{x}_0) = \mu - \boldsymbol{\gamma}^T \boldsymbol{\lambda},

where \mu is the Lagrange multiplier from the kriging system. Note that \gamma(0) = 0 by definition of the semivariogram, while the sill of the variogram, \lim_{\|\mathbf{h}\| \to \infty} \gamma(\mathbf{h}), is equivalent to the process variance (or nugget plus sill if microscale variation is present).^[20] This variance is always non-negative and achieves zero at observed locations, reflecting exact interpolation.^[21] In the context of Gaussian processes, simple kriging (with known zero mean) corresponds exactly to the conditional mean of a zero-mean Gaussian process given the observations, providing exact interpolation where the posterior mean passes through the data points and the posterior variance vanishes there. Ordinary kriging approximates this under an unknown constant mean, providing similar exact interpolation properties and serving as a form of empirical Bayesian inference under Gaussian assumptions with a constant mean.^[21]^[22]

Best Linear Unbiased Estimation

Kriging provides the best linear unbiased estimator (BLUE) for predicting the value of a spatial random field Z(\mathbf{x}) at an unobserved location \mathbf{x}_0, based on observations Z(\mathbf{x}_i) for i = 1, \dots, n. The estimator takes the linear form \hat{Z}(\mathbf{x}_0) = \sum_{i=1}^n \lambda_i Z(\mathbf{x}_i), where the weights \lambda = (\lambda_1, \dots, \lambda_n)^T are chosen to satisfy two key criteria: unbiasedness, meaning E[\hat{Z}(\mathbf{x}_0) - Z(\mathbf{x}_0)] = 0 (which implies \sum_{i=1}^n \lambda_i = 1 under stationarity assumptions), and minimum variance of the prediction error among all linear unbiased estimators. This optimality follows from the Gauss-Markov theorem applied to spatial processes, ensuring the kriging predictor has the smallest mean squared error in the class of linear estimators.^[3] The derivation of the BLUE weights proceeds by minimizing the prediction error variance \text{Var}(\hat{Z}(\mathbf{x}_0) - Z(\mathbf{x}_0)) subject to the unbiasedness constraint \sum_{i=1}^n \lambda_i = 1. This variance is expressed as

\text{Var}(\hat{Z}(\mathbf{x}_0) - Z(\mathbf{x}_0)) = \sum_{i=1}^n \sum_{j=1}^n \lambda_i \lambda_j C(\mathbf{x}_i, \mathbf{x}_j) - 2 \sum_{i=1}^n \lambda_i C(\mathbf{x}_i, \mathbf{x}_0) + C(\mathbf{x}_0, \mathbf{x}_0),

where C(\cdot, \cdot) denotes the covariance function. To solve this constrained optimization, the method of Lagrange multipliers is employed, forming the Lagrangian

L(\boldsymbol{\lambda}, \mu) = \sum_{i=1}^n \sum_{j=1}^n \lambda_i \lambda_j C(\mathbf{x}_i, \mathbf{x}_j) - 2 \sum_{i=1}^n \lambda_i C(\mathbf{x}_i, \mathbf{x}_0) + C(\mathbf{x}_0, \mathbf{x}_0) + 2\mu \left(1 - \sum_{i=1}^n \lambda_i\right).

Taking partial derivatives with respect to each \lambda_k and setting them to zero yields

\sum_{j=1}^n \lambda_j C(\mathbf{x}_j, \mathbf{x}_k) - C(\mathbf{x}_k, \mathbf{x}_0) + \mu = 0, \quad k = 1, \dots, n,

along with the constraint equation. In matrix form, this results in the kriging system

\begin{pmatrix} \mathbf{C} & \mathbf{1} \\ \mathbf{1}^T & 0 \end{pmatrix} \begin{pmatrix} \boldsymbol{\lambda} \\ \mu \end{pmatrix} = \begin{pmatrix} \mathbf{c}_0 \\ 1 \end{pmatrix},

where \mathbf{C} is the n \times n covariance matrix with entries C(\mathbf{x}_i, \mathbf{x}_j), \mathbf{c}_0 is the vector with entries C(\mathbf{x}_i, \mathbf{x}_0), and \mathbf{1} is a vector of ones. Solving this system provides the optimal weights that minimize the error variance while ensuring unbiasedness.^[23]^[3] The BLUE property of kriging exhibits invariance under affine transformations of the data. Specifically, if the observed values are transformed as Z'(\mathbf{x}_i) = a Z(\mathbf{x}_i) + b for constants a \neq 0 and b, the kriging estimator transforms accordingly to \hat{Z}'(\mathbf{x}_0) = a \hat{Z}(\mathbf{x}_0) + b, preserving unbiasedness and minimum variance in the transformed space. This robustness stems from the linear structure and the use of second-moment properties.^[3] Compared to other linear estimators, such as those based on geometric distances (e.g., inverse distance weighting), kriging achieves superior variance reduction by explicitly accounting for the spatial covariance structure, which captures dependencies beyond mere proximity. This incorporation of the covariance function ensures that the prediction error variance is lower, providing more reliable uncertainty quantification in spatially correlated data.^[23]

Kriging Variants

Simple Kriging

Simple kriging is a foundational geostatistical method for interpolating values at unsampled locations in a spatial field, assuming the underlying process has a known constant mean and exhibits second-order stationarity. Developed as part of the early geostatistical framework, it minimizes the mean squared prediction error among linear unbiased estimators under these conditions.^[24] This approach is particularly suited to scenarios where prior knowledge of the global mean allows for straightforward application without additional estimation steps.^[25] The method relies on two primary assumptions: the spatial random field possesses a constant known mean \mu, and it is second-order stationary, such that the expected value is constant and the covariance between observations depends solely on their spatial separation h.^[24] These assumptions ensure that the covariance structure fully captures the spatial dependence, enabling reliable predictions without needing to model local variations in the mean.^[26] Given n observations Z(x_1), \dots, Z(x_n) at sampled locations \mathbf{x} = (x_1, \dots, x_n), the simple kriging predictor \hat{Z}(x_0) at an unsampled location x_0 is formulated as the known mean plus a weighted sum of deviations from that mean:

\hat{Z}(x_0) = \mu + \mathbf{c}^T \mathbf{C}^{-1} (\mathbf{Z} - \mu \mathbf{1}),

where \mathbf{Z} is the vector of observations, \mathbf{1} is the vector of ones, \mathbf{c} is the n \times 1 vector of covariances between x_0 and the sampled locations, \mathbf{C} is the n \times n covariance matrix among the sampled locations, and the weights \boldsymbol{\lambda} = \mathbf{C}^{-1} \mathbf{c} for the centered process satisfy \sum \lambda_i = 1.^[25] This is equivalent to \hat{Z}(x_0) = \sum_{i=1}^n \lambda_i Z(x_i). The associated prediction variance, which quantifies uncertainty, is

\sigma^2_{SK}(x_0) = C(0) - \mathbf{c}^T \mathbf{C}^{-1} \mathbf{c},

with C(0) denoting the process variance at lag zero.^[24] Simple kriging offers advantages in simplicity and efficiency compared to more general variants when the mean is reliably known from prior data, reducing computational demands by avoiding mean estimation.^[26] Additionally, when the spatial process is Gaussian, simple kriging provides the exact conditional expectation, making it optimal beyond mere linearity. A representative application involves interpolating rainfall in a uniform climatic region, where the global mean is derived from historical averages across the field, allowing simple kriging to leverage the covariance structure for precise spatial predictions.^[27]

Ordinary Kriging

Ordinary kriging is a geostatistical technique for predicting values at unsampled locations in a spatial process where the mean is assumed to be constant but unknown across the domain. It achieves best linear unbiased prediction by solving for optimal weights that minimize the prediction variance while enforcing an unbiasedness constraint through a Lagrange multiplier. This method was formalized as part of the foundational geostatistical framework developed by Georges Matheron in the 1960s and 1970s.^[28]^[29] The primary assumptions of ordinary kriging are that the spatial process has an unknown constant mean and satisfies second-order stationarity, implying that the mean is constant and the covariance between any two points depends solely on their separation distance. These assumptions ensure that the spatial dependence structure can be modeled reliably using a covariance or variogram function, allowing for the estimation of the process without prior information on the mean value.^[30]^[31] To obtain the kriging weights \boldsymbol{\lambda}, the method solves an augmented linear system that incorporates the unbiasedness constraint \sum \lambda_i = 1:

\begin{bmatrix} \boldsymbol{\Gamma} & \mathbf{1} \\ \mathbf{1}^T & 0 \end{bmatrix} \begin{bmatrix} \boldsymbol{\lambda} \\ \mu \end{bmatrix} = \begin{bmatrix} \boldsymbol{\gamma} \\ 1 \end{bmatrix},

where \boldsymbol{\Gamma} is the variogram matrix among the sampled locations (with \Gamma_{ij} = \gamma(\mathbf{x}_i - \mathbf{x}_j)), \mathbf{1} is the vector of ones, \boldsymbol{\gamma} is the vector of variogram values between the sampled locations and the prediction point \mathbf{x}_0 (with \gamma_i = \gamma(\mathbf{x}_0 - \mathbf{x}_i)), \boldsymbol{\lambda} contains the weights, and \mu is the Lagrange multiplier enforcing the constraint. The prediction at \mathbf{x}_0 is then \hat{Z}(\mathbf{x}_0) = \boldsymbol{\lambda}^T \mathbf{Z}, where \mathbf{Z} is the vector of observed values.^[30]^[31]^[29] The associated prediction variance, or kriging variance, quantifies the uncertainty and is calculated as \sigma^2_{OK}(\mathbf{x}_0) = \boldsymbol{\gamma}^T \boldsymbol{\lambda} + \mu, where \gamma(0) = 0. This variance expression accounts for the spatial structure and the enforced constraint, providing a measure of reliability for the prediction.^[30]^[31] Ordinary kriging approximates simple kriging when the sample mean is close to the true mean, as the Lagrange multiplier \mu then aligns closely with the known mean assumption in simple kriging, reducing the additional complexity of the constraint.^[32]^[30] A classic application of ordinary kriging is in ore grade estimation within mining, where sample data from drill holes are used to predict metal concentrations at unsampled sites without assuming prior knowledge of the average grade, thereby supporting resource evaluation and mine planning. For instance, in gold deposit modeling, ordinary kriging weights nearby samples to estimate grades while ensuring the predictions are unbiased and variance-minimized, leading to more accurate block models for extraction decisions.^[29]^[30]

Universal Kriging

Universal kriging extends the kriging framework to handle spatial processes with non-stationary means, where the expected value varies systematically according to a known functional form rather than being constant. Developed by Georges Matheron, this method models the process as Z(\mathbf{x}) = \mu(\mathbf{x}) + Y(\mathbf{x}), with \mu(\mathbf{x}) = \mathbf{X}(\mathbf{x}) \boldsymbol{\beta} representing the deterministic trend—typically a linear combination of basis functions such as polynomials in coordinates (e.g., constant, linear, or quadratic terms)—and Y(\mathbf{x}) denoting the zero-mean stochastic residuals. The coefficients \boldsymbol{\beta} are unknown and estimated from the data, while the residuals are assumed to be second-order stationary with a known covariance structure or variogram, ensuring the covariance of residuals depends only on spatial separation. The core of universal kriging involves solving a generalized linear system that simultaneously determines the interpolation weights and the trend parameters to achieve the best linear unbiased predictor (BLUP). In matrix notation using the covariance form, the system is given by

\begin{pmatrix} \boldsymbol{\Sigma} & \mathbf{F} \\ \mathbf{F}^T & \mathbf{0} \end{pmatrix} \begin{pmatrix} \boldsymbol{\omega} \\ \boldsymbol{\mu} \end{pmatrix} = \begin{pmatrix} \mathbf{c}_0 \\ \mathbf{f}(\mathbf{x}_0) \end{pmatrix},

where \boldsymbol{\Sigma} is the n \times n covariance matrix of the observations Z(\mathbf{s}_1), \dots, Z(\mathbf{s}_n), \mathbf{F} is the n \times p design matrix with rows \mathbf{X}(\mathbf{s}_i)^T ( p basis functions), \mathbf{c}_0 is the n \times 1 vector of covariances between Z(\mathbf{x}_0) and the observations, \mathbf{f}(\mathbf{x}_0) = \mathbf{X}(\mathbf{x}_0), \boldsymbol{\omega} are the weights, and \boldsymbol{\mu} are the Lagrange multipliers associated with the trend constraints. The predictor is then \hat{Z}(\mathbf{x}_0) = \boldsymbol{\omega}^T \mathbf{Z} + \mathbf{f}(\mathbf{x}_0)^T \hat{\boldsymbol{\beta}}, where \hat{\boldsymbol{\beta}} = -\boldsymbol{\mu}, ensuring unbiasedness by satisfying \mathbf{F}^T \boldsymbol{\omega} = \mathbf{f}(\mathbf{x}_0). An equivalent variogram-based system replaces covariances with the variogram matrix \boldsymbol{\Gamma} and vector \boldsymbol{\gamma}_0. The kriging variance for universal kriging, which quantifies prediction uncertainty, is

\sigma^2_{UK}(\mathbf{x}_0) = C(0) - \mathbf{c}_0^T \boldsymbol{\omega} - \boldsymbol{\mu}^T \mathbf{f}(\mathbf{x}_0),

where C(0) is the process variance; in variogram terms, it becomes \sigma^2_{UK}(\mathbf{x}_0) = \boldsymbol{\gamma}_0^T \boldsymbol{\omega} + \boldsymbol{\mu}^T \mathbf{f}(\mathbf{x}_0) under intrinsic stationarity assumptions for residuals. This variance accounts for both the spatial correlation of residuals and the uncertainty in the trend estimate, generally exceeding that of ordinary kriging due to the additional variability from \boldsymbol{\beta}. An alternative implementation involves trend removal: first fit the trend via generalized least squares using the residual covariance structure, subtract it from observations to yield residuals, apply ordinary kriging to the residuals, and add the fitted trend back at \mathbf{x}_0; this is computationally equivalent under the model assumptions. A practical example arises in topographic modeling, where elevation data exhibit a linear trend due to underlying geology. The mean is modeled as \mu(x, y) = \beta_0 + \beta_1 x + \beta_2 y, with basis functions \mathbf{X}(x, y) = (1, x, y)^T, and residuals assumed to follow a spherical variogram with nugget effect, sill, and range fitted from data. For a dataset of elevations at irregular points, universal kriging yields predictions that capture both the regional slope and local fluctuations, with variances reflecting sparser trend information in peripheral areas.

Other Variants

Indicator kriging is a geostatistical technique designed for handling categorical, binary, or non-normal continuous data by transforming the variable into a set of binary indicators at specified thresholds, allowing estimation of conditional cumulative distribution functions via ordinary kriging applied to each indicator variogram.^[33] This approach enables probabilistic modeling of spatial distributions without assuming normality, making it suitable for applications like ore grade estimation where high and low values exhibit asymmetric continuity.^[33] Multiple indicators can be combined to recover the full conditional distribution at unsampled locations, providing both point estimates and uncertainty measures for decision-making in resource evaluation.^[33] Co-kriging extends the kriging framework to multivariate spatial data where variables are correlated, incorporating cross-covariances between primary and secondary variables to improve prediction accuracy beyond univariate kriging.^[29] By solving a system of equations that includes auto- and cross-variograms, co-kriging leverages auxiliary data—such as densely sampled secondary variables—to enhance estimates of the target variable, particularly when direct observations are sparse.^[29] This method is widely applied in scenarios like environmental monitoring, where correlated pollutants or geochemical indicators are jointly analyzed for more precise spatial interpolation.^[29] Bayesian kriging integrates prior distributions on model parameters, such as the mean, trend, and covariance structure, to derive posterior predictions, offering a probabilistic framework that contrasts with the frequentist approach of classical kriging by explicitly accounting for uncertainty in hyperparameters.^[34] Markov chain Monte Carlo (MCMC) methods are commonly employed to sample from the posterior, enabling flexible incorporation of expert knowledge or historical data into the spatial prediction process.^[34] This variant is particularly useful in complex settings with limited data, where it provides full posterior distributions for predictions rather than point estimates and variances.^[34] Log-normal kriging addresses positively skewed data, such as mineral grades or rainfall amounts, by applying a logarithmic transformation to achieve approximate normality before performing kriging, followed by a bias-corrected back-transformation to the original scale.^[35] Disjunctive kriging, a related nonlinear extension, further handles skewed distributions using expansions like Hermite polynomials to model the variable as a sum of uncorrelated factors, allowing estimation of nonlinear functions of the spatial variable without transformation.^[36] These methods preserve the log-normal or multimodal characteristics of the data, improving accuracy in applications involving multiplicative processes or bounded positive values.^[35]^[36] Block kriging modifies the standard point kriging system to estimate average values over finite areas or volumes, such as mine blocks, by adjusting the variogram to account for the support size and shape of the block relative to point samples. This involves computing block-to-point covariances or using regularization in the variogram model, which reduces smoothing effects and provides variance estimates tailored to the larger support, essential for resource recovery planning. It is routinely used in mining to assess tonnage and grade over practical exploitation units, ensuring unbiased predictions that reflect the averaging inherent in block sampling.

Computation and Implementation

Estimation Steps

The estimation of spatial predictions using kriging follows a structured workflow that ensures the method's assumptions are met and its results are reliable. This process begins with preparing the input data and culminates in generating predictions along with measures of uncertainty. The steps are iterative, often requiring validation and adjustment to account for the spatial structure of the data.^[24] Step 1: Data Exploration and Declustering. Initial data preparation involves exploratory analysis to identify outliers, trends, and sampling irregularities, such as uneven density that can bias estimates. For uneven sampling, declustering techniques are applied to weight observations inversely proportional to local density, providing a more representative global mean and reducing overrepresentation from clustered points. This step ensures stationarity assumptions are reasonably satisfied before proceeding.^[24]^[25] Step 2: Variogram Modeling. Next, the spatial dependence structure is modeled using an empirical variogram, calculated from pairwise differences in observed values as a function of separation distance. A theoretical variogram model (e.g., spherical or exponential) is fitted to the empirical points, often through least-squares or maximum likelihood methods. Validation occurs via cross-validation, where each point is temporarily removed, predicted from neighbors, and errors assessed for bias and variance; this confirms the model's adequacy for capturing spatial continuity. As detailed in the variogram and covariance functions section, this modeling is crucial for determining covariance weights.^[24]^[25] Step 3: Choosing the Kriging Variant and Setting Up the System Matrix. The appropriate kriging variant is selected based on knowledge of the mean: ordinary kriging is typically chosen when the mean is unknown but assumed constant, as it incorporates a Lagrange multiplier to enforce unbiasedness without requiring a prior mean estimate. The system of equations is then formulated as a matrix, where the covariance matrix between observed points and the target location defines the weights for linear combination of observations.^[24]^[25] Step 4: Solving the Kriging Equations. The kriging weights are obtained by solving the linear system, typically via direct matrix inversion for small datasets (e.g., fewer than 100 points), which yields exact solutions. For larger datasets, iterative methods like conjugate gradient solvers are employed to avoid computational instability and high memory demands.^[24]^[25] Step 5: Generating Predictions and Variance Maps; Diagnostic Checks. Predictions at unsampled locations are computed as the weighted sum of neighboring observations, while the kriging variance—derived from the model—provides a map of prediction uncertainty, highlighting areas of high reliability or extrapolation risk. Diagnostic checks include verifying that the mean prediction error approximates zero and that the standardized errors follow a standard normal distribution, often through leave-one-out cross-validation to assess overall performance.^[24]^[25] For large datasets, approximations such as kriging with moving neighborhoods are used, where predictions at each location incorporate only a local subset of nearby points (e.g., 10-50) within a defined search radius and anisotropy, balancing accuracy with computational efficiency. This approach mitigates the quadratic scaling of full kriging while preserving local spatial structure.^[24]^[25]

Software Tools

Kriging implementations are available in various open-source software libraries, particularly in statistical programming environments. In R, the gstat package provides tools for variogram fitting and spatial predictions using methods such as ordinary and universal kriging, along with support for spatio-temporal extensions.^[37] The geoR package extends this capability with Bayesian geostatistical analysis, enabling uncertainty quantification in model parameters for Gaussian processes akin to kriging.^[38] In Python, scikit-learn offers Gaussian process regression, which mathematically aligns with kriging for spatial interpolation and prediction under stationarity assumptions.^[39] Additionally, PyKrige serves as a dedicated toolkit for 2D and 3D ordinary and universal kriging, incorporating standard variogram models like spherical and exponential. Commercial software caters to professional applications, especially in resource estimation and visualization. Surfer, developed by Golden Software, facilitates 2D and 3D mapping through kriging gridding, allowing users to create contoured surfaces from irregularly spaced data.^[40] Isatis.neo by Geovariances supports comprehensive mining workflows, including kriging for mineral resource estimation with features for domaining and simulation-based uncertainty assessment.^[41] Similarly, the Groundwater Modeling System (GMS) from Aquaveo integrates kriging routines based on the Geostatistical Software Library (GSLIB) for subsurface modeling in hydrogeological contexts.^[42] These tools commonly include automated variogram modeling to fit experimental semivariograms, streamlining the estimation of spatial covariance structures.^[43] Support for kriging variants, such as co-kriging, enables multivariate predictions by incorporating auxiliary variables with cross-variograms.^[44] Visualization of uncertainty is a key feature, often through prediction variance maps that highlight prediction reliability across the spatial domain.^[27] Open-source options like gstat, geoR, scikit-learn, and PyKrige are freely accessible to academics and researchers, promoting widespread adoption in educational and non-commercial settings. Commercial tools such as Surfer and Isatis.neo require licensing but offer robust support for industry-scale computations. Many integrate with geographic information systems (GIS); for instance, ArcGIS provides the Geostatistical Analyst extension for kriging within a broader spatial analysis framework.^[45] Recent trends emphasize cloud-based platforms for handling large datasets. Google Earth Engine supports kriging interpolation on feature collections, enabling scalable environmental applications like temperature mapping over vast regions without local computational limits.^[46]

Applications

Mining and Earth Sciences

In mining and earth sciences, kriging is extensively applied for ore reserve estimation, where block kriging is particularly valued for interpolating grades across discretized blocks that represent smelting units, thereby providing averaged estimates that reduce uncertainty in tonnage and grade calculations.^[47] This method accounts for spatial variability in mineral deposits by minimizing estimation variance through variogram-based weighting, enabling more reliable predictions of recoverable resources compared to simpler interpolation techniques.^[48] For instance, in iron ore deposits like Choghart in Iran, block kriging has been used to construct three-dimensional models that integrate drillhole data, yielding tonnage estimates with quantified errors that inform mine planning and economic viability assessments.^[47] Grade control during excavation relies heavily on ordinary kriging to deliver real-time estimates of mineral grades, allowing operators to optimize blasting patterns and selective mining that minimize dilution and ore loss.^[49] By applying ordinary kriging to closely spaced drillhole samples as mining progresses, estimators can delineate high-grade zones with sufficient precision to guide excavation decisions, such as adjusting blast designs to target ore selectively while avoiding waste.^[50] This approach is common in open-pit operations, where rapid kriging computations support on-site decisions that enhance recovery rates and reduce operational costs, as demonstrated in various gold and base metal mines.^[51] In petroleum exploration and reservoir characterization, co-kriging integrates sparse well data with densely sampled seismic attributes to estimate properties like porosity, improving the spatial continuity of models for hydrocarbon reservoirs.^[52] This multivariate extension of kriging leverages cross-correlations between primary (e.g., well-logged porosity) and secondary (e.g., seismic impedance) variables, resulting in more accurate interpolations across unsampled areas and better delineation of productive zones.^[53] Applications in fields like those in southwest Iran have shown co-kriging to outperform univariate methods by incorporating seismic trends, leading to refined reservoir simulations that support enhanced oil recovery strategies.^[52] A seminal case study traces kriging's origins to the Witwatersrand gold mines in South Africa during the 1950s, where D. G. Krige developed early distance-weighted averaging techniques to estimate gold grades from drillhole data, laying the foundation for modern geostatistics.^[54] This work evolved into routine applications by the early 1960s on large Anglovaal gold mines, marking one of the first industrial uses of kriging for reserve evaluation in the region. Today, these methods persist in constructing three-dimensional geological models for Witwatersrand operations, integrating kriging with seismic data to map complex reef structures and update reserves amid ongoing extraction.^[55] One key benefit of kriging in these contexts is its ability to quantify estimation risk through the kriging variance, which serves as a measure of local uncertainty and informs economic decision-making, such as cutoff grade optimization and investment risk assessment.^[56] By providing not only point or block estimates but also variance metrics, kriging enables mining engineers to evaluate the reliability of tonnage and grade projections, supporting probabilistic analyses that balance resource recovery against operational hazards.^[57] This risk quantification has become integral to standards like those from the Canadian Institute of Mining, ensuring that reserve statements reflect both mean values and their associated uncertainties for sustainable mine development.^[58]

Environmental and Hydrological Modeling

In environmental and hydrological modeling, kriging techniques are widely applied to interpolate sparse spatial data for simulating processes such as pollutant dispersion and water flow dynamics. Universal kriging, which accounts for underlying trends in the data, is particularly effective for estimating hydraulic conductivity fields in groundwater systems, where heterogeneity and regional drifts influence flow patterns. For instance, in aquifer studies, universal kriging has been used to generate continuous maps of hydraulic conductivity from borehole measurements, enabling more accurate simulations of groundwater movement and contaminant transport.^[59]^[60] For air quality assessment, indicator kriging is employed to map exceedance probabilities of pollutant thresholds, transforming concentration data into binary indicators (exceedance or non-exceedance) to better capture non-linear spatial relationships. This approach facilitates probabilistic risk maps for pollutants like PM_2.5 and NO₂, integrating monitoring station data to predict areas at risk of violating air quality standards. By providing variance estimates, indicator kriging quantifies uncertainty in these predictions, supporting regulatory decisions on emission controls.^[61]^[62]^[63] In climate modeling, simple kriging serves as a foundational method for interpolating temperature data across grids derived from weather station networks, assuming stationarity in the mean and variance. This technique produces smooth surfaces of monthly or annual temperatures, essential for downscaling climate models and assessing regional warming trends. Variogram models in simple kriging may incorporate anisotropy to reflect directional influences, such as prevailing winds affecting temperature gradients in hydrological catchments.^[64]^[65]^[66] A notable case study involves co-kriging to map nitrate contamination risks in vulnerable aquifers, where nitrate measurements from wells are jointly interpolated with auxiliary land-use data (e.g., agricultural intensity) to enhance prediction accuracy. In the Vega de Granada aquifer, Spain, compositional co-kriging revealed high-risk zones linked to intensive farming, with probability maps showing exceedance risks above 50 mg/L in 30-40% of the area, aiding targeted remediation efforts. This multivariate extension leverages correlations between nitrate levels and land-use patterns to reduce estimation errors in sparse datasets.^[67] Overall, kriging's advantages in these applications include its robustness to sparse monitoring networks typical of remote environmental sites, allowing reliable interpolations from limited observations, and its provision of kriging variance for uncertainty quantification, which is critical for risk-based hydrological assessments like flood or contamination forecasting.^[68]^[69]^[70]

Engineering and Computer Experiments

In engineering and computer experiments, kriging serves as a surrogate modeling technique within the design and analysis of computer experiments (DACE) framework to emulate computationally expensive simulations, such as computational fluid dynamics (CFD) analyses, thereby facilitating efficient optimization processes.^[71] Introduced in seminal work on DACE, kriging models the deterministic output of a computer code as a realization of a stochastic process, enabling interpolation and prediction across the input space with associated uncertainty estimates.^[71] This approach is particularly valuable in engineering design, where direct evaluations of high-fidelity simulations are prohibitive due to time and resource constraints, allowing for rapid exploration of design spaces.^[72] Space-filling designs, such as Latin hypercube sampling (LHS), are commonly employed to generate initial datasets for fitting kriging models efficiently in DACE. LHS stratifies the input space into equally probable intervals, ensuring that sample points are maximally dispersed and representative of the design domain, which enhances the kriging model's predictive accuracy with fewer evaluations. For instance, in mechanical engineering applications, LHS combined with kriging reduces the number of required simulation runs while maintaining robust coverage of multidimensional parameter spaces, as demonstrated in surrogate-based optimization workflows.^[73] In robust optimization, kriging's prediction variance is leveraged to propagate uncertainties in mechanical designs, enabling the identification of solutions that minimize sensitivity to input variations. By incorporating the kriging variance as a measure of model uncertainty and noise propagation, engineers can formulate objectives that balance mean performance with robustness, such as minimizing the variance of structural responses under parametric fluctuations.^[74] This statistics-based approach, using kriging metamodels to approximate both mean and variance, has been applied to optimize components like cantilever beams, where it efficiently handles noisy evaluations and achieves designs with improved reliability.^[74] A representative case study involves aerodynamic shape optimization, where kriging predicts lift coefficients from sparse CFD simulations of airfoil geometries. In optimizing transonic airfoils, such as the RAE 2822, kriging surrogates trained on a limited set of high-fidelity evaluations (e.g., 50-100 points) enable global search algorithms to maximize lift-to-drag ratios, reducing the overall computational cost by over 90% compared to direct evaluations while converging to near-optimal shapes.^[75] Extensions to multi-fidelity modeling via co-kriging further enhance efficiency by combining low- and high-resolution simulations. Co-kriging constructs a hierarchical surrogate that uses inexpensive low-fidelity data to inform the high-fidelity model, scaling the correlation structure across fidelity levels to improve predictions in engineering optimization. For example, in structural design problems, this approach integrates coarse mesh finite element analyses with refined ones, achieving accuracy comparable to full high-fidelity kriging with significantly fewer expensive evaluations.^[76]

Limitations and Extensions

Assumptions and Challenges

Kriging methods rely on several foundational assumptions to ensure the validity of their predictions as best linear unbiased estimators (BLUE). A primary assumption is second-order stationarity of the underlying spatial process, which requires a constant mean across the domain and a covariance structure that depends solely on the lag distance between points, allowing the variogram or covariance function to model spatial dependence consistently.^[77] For simple kriging, an additional assumption of Gaussianity—or normality of the spatial random field—is often invoked to achieve optimality in the mean squared prediction error sense, though ordinary kriging maintains unbiasedness without full normality by estimating the mean from the data.^[78] Violations of these assumptions, such as non-stationarity or non-Gaussian distributions, can propagate errors through the interpolation process.^[2] Equally critical is the assumption of a correctly specified variogram model, which captures the spatial autocorrelation structure; this model directly determines the weights assigned to observed data points in the kriging predictor. Misspecification of the variogram—such as choosing an inappropriate functional form (e.g., exponential versus spherical) or inaccurate parameter estimation—leads to suboptimal weights and biased predictions, potentially underestimating variability or introducing systematic errors in the kriging surface.^[79] Kriging also implicitly assumes the absence of significant outliers in the dataset, as anomalous points can inflate the nugget effect in the variogram or skew weights, resulting in distorted predictions that fail to represent true spatial patterns.^[26] When these assumptions hold, kriging provides reliable interpolations, but real-world data often deviate, necessitating careful model diagnostics. Computational challenges further complicate kriging applications, particularly the O(n³) time complexity required for inverting the n × n covariance matrix during weight estimation, which becomes infeasible for datasets with thousands of points and demands substantial memory and processing resources.^[80] This issue is exacerbated by ill-conditioned matrices, often arising from collinear or densely clustered data points, which cause near-singularity and numerical instability, leading to unreliable or divergent solutions unless regularization techniques are applied.^[81] Non-ergodicity presents a theoretical hurdle, as the finite sample variogram may not ergodically represent the population variogram, resulting in over-smoothed kriging estimates that regress extremes toward the mean and underrepresent local heterogeneity or sharp gradients in the spatial field.^[82] Boundary effects in finite spatial domains introduce additional practical challenges, manifesting as bias near edges where fewer surrounding observations are available to inform predictions, often causing underprediction of variability or directional artifacts in the interpolated surface.^[83] To evaluate adherence to assumptions and overall model performance, cross-validation diagnostics are essential; leave-one-out cross-validation, for instance, assesses unbiasedness by verifying that the mean error (ME) approximates zero and optimizes accuracy by minimizing the root mean square error (RMSE), providing quantitative measures of prediction reliability without requiring independent validation data.^[84]

Modern Developments

Since the 2000s, kriging has seen significant advancements in scalability to handle big data, particularly through approximations like fixed rank kriging, which employs low-rank covariance structures to reduce computational complexity for datasets exceeding millions of points, enabling efficient predictions on massive spatial grids.^[85] Cluster-based kriging further enhances this by partitioning large datasets into subsets for local modeling, followed by aggregation, which has proven effective for high-dimensional problems while maintaining prediction accuracy.^[86] Sparse Gaussian process approximations, akin to kriging, leverage inducing points to approximate full covariance matrices, facilitating applications to satellite imagery analysis where data volumes from remote sensing exceed traditional limits.^[87] Hybrid approaches integrating kriging with machine learning have addressed non-stationarity in spatial data, such as random forest residuals kriging, where a random forest model captures nonlinear trends and kriging interpolates the residuals, improving predictions for complex environmental variables like soil organic matter.^[88] Deep Gaussian processes, an extension of kriging's Gaussian process foundation, stack multiple layers to model hierarchical non-stationarities, offering superior flexibility for geospatial regression tasks compared to single-layer kriging.^[89] Universal kriging has been applied in climate change modeling on global datasets to incorporate trends from general circulation models. Bayesian variants of kriging have been used to map COVID-19 spatial spread, employing empirical Bayesian kriging to interpolate incidence rates across regions while accounting for uncertainty in underreported data.^[90] For uncertainty quantification, ensemble kriging methods combine multiple kriging models with other surrogates like artificial neural networks, yielding robust variance estimates in AI-driven predictions and reducing overfitting in high-stakes scenarios.^[91] In the 2020s, trends emphasize kriging's integration with deep learning for multi-scale spatial prediction in environmental AI, such as deep kriging neural networks that fuse convolutional layers with Gaussian processes to enhance resolution in remote sensing tasks like land cover mapping.^[92]

References

[1]
A practical primer on geostatistics - USGS Publications Warehouse
Jul 6, 2009 · Geostatistics characterizes incompletely known spatial systems using numerical techniques and probabilistic models, using every measurement's ...
[2]
Kriging Interpolation Explanation | Columbia Public Health
Kriging is a method of spatial interpolation that originated in the field of mining geology as is named after South African mining engineer Danie Krige.
[3]
[PDF] "Kriging" in - UC Davis Statistics
Kriging, at its most fundamental level, is an interpola- tion method used to convert partial observations of a spatial field to predictions of that field at ...
[4]
An Experimental Comparison of Ordinary and Universal Kriging and ...
Among numerous findings, the most striking was that the two kriging methods were substantially superior to the inverse distance weighting methods over all ...
[5]
https://arxiv.org/abs/2408.02331
[6]
[PDF] Introduction to Choosing a Kriging Plan - Geostatistics Lessons
Kriging is the primary technique for the estimation of grades. Kriging is a linear unbiased estimator that minimizes the estimation variance using a site- ...
[7]
A statistical approach to some basic mine valuation problems on the ...
A statistical approach to some basic mine valuation problems on the Witwatersrand, by D.G. Krige, published in the Journal, December 1951 : introduction by the ...Missing: thesis | Show results with:thesis
[8]
Krige, D.G. (1951) A Statistical Approaches to Some Basic Mine ...
Krige, D.G. (1951) A Statistical Approaches to Some Basic Mine Valuation Problems on the Witwatersrand. Journal of the Chemical, Metallurgical and Mining ...
[9]
[PDF] IN MEMORY OF PROFESSOR DANIE KRIGE - Gold Fields
the first application of Kriging. It can be labeled 'simple elementary. Kriging', being based on the spatial correlation between the peripheral values and ...
[10]
GEORGES MATHERON: FOUNDER OF SPATIAL STATISTICS - jstor
Originally a geologist, Professor Georges Matheron (1930-2000) became increasingly mathematical in his approach to problem-solving.
[11]
Danie Krige - SciELO South Africa
The term 'géostatistique' was introduced by Matheron in 1962, as was the term 'kriging', although the latter appeared as a French word in a footnote.
[12]
Matheron, G. (1962) Trait de gostatistique applique, vol 14. Editions ...
Dec 18, 2024 · Matheron, G. (1962) Traité de géostatistique appliquée, vol 14. Editions Technip. ... ABSTRACT: Geostatistical Kriging is performed on hydrologic ...Missing: Georges | Show results with:Georges
[13]
[PDF] Geostatistics Without Tears
Jul 5, 2006 · environmental agencies. In the mid 1980's the Environmental Protection Agency (EPA) commissioned a geostatistical software package, GEO-EAS ...
[14]
Fifty Years of Kriging | SpringerLink
Jun 26, 2018 · Random function models and kriging constitute the core of the geostatistical methods created by Georges Matheron in the 1960s and further ...
[15]
The intrinsic random functions and their applications
Jul 1, 2016 · The intrinsic random functions and their applications. Published online by Cambridge University Press: 01 July 2016. G. Matheron.
[16]
A space and time scale‐dependent nonlinear geostatistical ...
Jul 21, 2015 · A geostatistical approach to downscaling climate model data is presented; Downscaled precipitation and temperature reproduce properties of ...<|separator|>
[17]
Non-stationary variogram models for geostatistical sampling ...
There are several possible approaches for implementing a non-stationary variogram. Three common approaches are (i) segmentation, (ii) locally adaptive Kriging ...
[18]
https://link.springer.com/article/10.1007/BF01032109
[19]
Fitting variogram models by weighted least squares
Carroll, R. J. and Ruppert, D., 1982, A comparison between maximum likelihood and generalized least squares in a heteroscedastic linear model:Jour. Amer. Stat.
[20]
[PDF] Interpolation, Kriging, Gaussian Processes - Duke People
The surfaces of equation (51), plotted in Figure 4, were interpolated with IDW inter- polation with α = 0.1 and q = 3. The Kriging and Gaussian Process ...
[21]
[PDF] A Note on Kriging and Gaussian Processes - DigitalCommons@USU
Gaussian Processes (GP) define a distribution over functions, while Kriging is a spatial interpolation method based on GP modeling.
[22]
[PDF] Best Linear Unbiased Estimation and Kriging - Alert Geomaterials
Best Linear Unbiased Estimation (BLUE). To make the estimator error as small as possible, its mean should be zero and its variance minimal. The mean is ...
[23]
[PDF] A Practical Primer on Geostatistics - USGS Publications Warehouse
These assumptions make simple kriging the most restricted form of kriging ... Ordinary kriging normal equations for optimal weights : Ordinary kriging ...
[24]
[PDF] Introduction to Geostatistics — Course Notes - University of Wyoming
fitting a permissible mathematical function to the experimental variogram; (3) conducting kriging interpolation based on this function. In the above example ...
[25]
Chapter 14 Kriging | Spatial Statistics for Data Science - Paula Moraga
Kriging (Matheron 1963) is a spatial interpolation method used to obtain predictions at unsampled locations based on observed geostatistical data.Missing: history | Show results with:history
[26]
How Kriging works—ArcGIS Pro
Kriging is an advanced geostatistical procedure that generates an estimated surface from a scattered set of points with z-values. Unlike other interpolation ...
[27]
[PDF] MATHERON - Paris
Economic Geology. Vol. 58, 1963, pp. 1246–1266. PRINCIPLES OF GEOSTATISTICS. G. MATHERON. ABSTRACT. Knowledge of ore grades and ore reserves as well as error ...
[28]
[PDF] Mining Geostatistics
The distribution of ore grades within a deposit is of mixed character, being partly structured and partly random. On one hand, the mineralizing process.
[29]
[PDF] Statistics for Spatial Data
Statistics for spatial and temporal data would provide dynamic models for phenomena distributed through space and evolving in time. Onward into the next decade!
[30]
[PDF] Kriging methods in spatial statistics - mediaTUM
Cressie, N. A. C. (1990). The Origins of Kriging. Mathematical Geology 22(3), 239–252. Cressie, N. A. C. (1993). Statistics for Spatial Data. Wiley Series ...
[31]
[PDF] 1.2 Kriging - University of Washington Department of Statistics
Gaussian process. µ(s)=EZ(s) Var Z(s) < ∞. Z is strictly stationary if. Z is ... when µ and C are known (simple kriging). The prediction variance is p(X) ...
[32]
Nonparametric estimation of spatial distributions
May 25, 1982 · Such rich structural information allows a nonparametric risk-qualified, estimation of local and global spatial distributions. Article PDF ...
[33]
A Bayesian Analysis of Kriging: Technometrics: Vol 35, No 4
Mar 12, 2012 · Technometrics Volume 35, 1993 - Issue 4 ... A Bayesian Analysis of Kriging. Mark S. Handcock Department of Statistics and Operations Research, ...
[34]
The lognormal approach to predicting local distributions of selective ...
Journel, A. G., 1977, Kriging in terms of projections: Math. Geol., v. 69 ... Journel, A.G. The lognormal approach to predicting local distributions of selective ...
[35]
A Simple Substitute for Conditional Expectation : The Disjunctive ...
In this paper, a new procedure for non linear estimation is proposed: it is better than the usual best linear estimation, and necessitates less ...
[36]
[PDF] gstat: Spatial and Spatio-Temporal Geostatistical Modelling ...
Function krigeST is a R implementation of the kriging function from gstat using spatio-temporal covariance models following the implementation of krige0.
[37]
[PDF] geoR.pdf
The geoR package is for geostatistical analysis, including variogram-based, likelihood-based and Bayesian methods.
[38]
1.7. Gaussian Processes - Scikit-learn
Gaussian Processes (GP) are a nonparametric supervised learning method used to solve regression and probabilistic classification problems.Gaussian Processes · GaussianProcessRegressor · 1.8. Cross decomposition · RBF
[39]
Surfer | Products - Golden Software
Surfer transforms complex geospatial data into 2D/3D models, with analysis tools, and creates maps quickly. It is used by scientists and engineers.Free Trial · Downloads & Installs · Grapher · PricingMissing: commercial Isatis GMS
[40]
Isatis.neo | Geostatistics made accessible - Geovariances
Isatis.neo is advanced geostatistical software for exploring, analyzing, and visualizing spatial data, simplifying complex processes and building robust models.
[41]
[PDF] GMS TUTORIALS - Geostatistics – 2D - Aquaveo
Kriging is a set of linear regression routines that minimize estimation variance from a predefined covariance model. 1. Select the Interpolation | Interpolation ...Missing: software | Show results with:software
[42]
Variogram modelling for kriging in Surfer - a tutorial
Aug 1, 2025 · Variogram modeling characterizes spatial continuity of data by fitting a model to an experimental variogram, which is calculated from the data.Missing: commercial Isatis GMS
[43]
[PDF] Co-kriging with the gstat package of the R environment for statistical ...
Jul 14, 2018 · This technical note shows how to perform co-kriging using the gstat geostatistical package [12] of the R environment for statistical comput-.Missing: paper | Show results with:paper
[44]
ArcGIS Geostatistical Analyst - Spatial Interpolation Methods - Esri
ArcGIS Geostatistical Analyst provides the essential tools needed to predict and model spatial phenomena, making it a key resource for professionals in fields ...Missing: free academics
[45]
ee.FeatureCollection.kriging - Earth Engine - Google for Developers
Oct 6, 2023 · The kriging method on a FeatureCollection returns an Image representing the results of sampling a Kriging estimator at each pixel. · The method ...
[46]
Reserve estimation of central part of Choghart north anomaly iron ...
In the kriging process the sum of weights assigned to the input data is adjusted to one, and the error of estimation is kept to a minimum. Kriging utilizes the ...
[47]
3.2.3: Overview of Reserve Estimation Methods | MNG 230
This estimation is often accomplished using a technique known as kriging. Kriging provides an optimal interpolation using the variogram; and the technique is ...
[48]
[PDF] Quantitative Kriging Neighbourhood Analysis for the Mining Geologist
Kriging (Matheron, 1962, 1963a,. 1963b; Journel and Huijbregts, 1978) is also a necessary step in the main methods of conditional simulation used in the mining.<|separator|>
[49]
Introduction to Choosing a Kriging Plan - Geostatistics Lessons
Oct 30, 2015 · Typically simple kriging is used to reduce the impact of local means in sparsely sampled and peripheral areas which could influence our ...
[50]
Real-Time Mining
May 10, 2019 · ordinary Kriging based on real-time mining monitoring data. ... • Block 1: the next blast block that has been characterized by grade control.
[51]
Estimation of reservoir porosity using analysis of seismic attributes in ...
Jan 28, 2020 · Consolidated dependent kriging is one of a variety of co-kriging methods. It is used when secondary data (seismic) exist in all parts of the ...
[52]
Porosity from seismic data: A geostatistical approach | Geophysics
Mar 2, 2017 · Using a geostatistical technique called cokriging, the areal distribution of porosity is estimated first in a numerically simulated ...
[53]
Kriging: Understanding allays intimidation
Matheron anglicized the term to kriging when he published a paper for English-speaking readers. France dominated the development and application of ...
[54]
[PDF] Geostatistical applications in petroleum reservoir modelling - SAIMM
This paper briefly discusses porosity modelling by using kriging and sequential Gaussian simulation, and permeability modelling by using collocated co-.
[55]
[PDF] Uncertainty of Mineral Resource Estimates - Geovariances
Linear interpolation techniques like kriging provide kriging variances, a first measure of uncertainty. The kriging variance is the variance of the error “true ...
[56]
Uncertainty Quantification in Mineral Resource Estimation
Aug 11, 2024 · Linear kriging methods are famous for quantifying optimal weights based on minimizing expected error variance. Linear kriging expresses ...Geostatistical Techniques · Uncertainty At Sampled... · Uncertainty At Unsampled...
[57]
[PDF] CIM Estimation of Mineral Resources and Mineral Reserves Best ...
• the kriging variance or standard deviation of the block ... At an early stage of the Mineral Reserve estimation process, various mining methods should.
[58]
Applications of Universal Kriging to an Aquifer Study in New Jersey
Kriging was used to (1) estimate the altitude of an aquifer surface, (2) estimate hydraulic conductivities from point data, and (3) estimate the associated ...Missing: modeling | Show results with:modeling
[59]
[PDF] Geostatistical Analysis of Hydraulic Conductivity in Heterogeneous ...
Figure 1.1 Generalized chain of events for a groundwater modeling effort 1-3 ... such as universal kriging automatically calculate a trend, it ...
[60]
A pragmatic approach to estimate the number of days in ...
We use a kriging model to combine surface observations and the CHIMERE model. · Daily probabilities of exceedance are computed with a Gaussian hypothesis. · A ...
[61]
Uncertainty assessment of PM2.5 contamination mapping using ...
Apr 12, 2016 · The uncertainty assessment methods currently in use include the sequential indicator simulation (SIS) and indicator kriging techniques. However, ...
[62]
[PDF] Modeling threshold exceedance probabilities of spatially correlated ...
Jan 29, 2009 · First method is indicator kriging, that is spatial interpolation of the ... Guidance report on preliminary assessment under. EC air quality ...
[63]
Spatial interpolation of temperature in the United States using ...
This paper contributes to the literature by developing a new kriging model for interpolating the air temperature in the mainland of the US in 2010.
[64]
Spatial modeling and interpolation of monthly temperature using ...
Jun 13, 2025 · Interpolation by kriging was applied to: (1) the basic monthly temperature values for all 35 stations from. 1900 through May 1993, (2) the ...
[65]
How do I take groundwater flow direction into consideration when ...
How do I take groundwater flow direction into consideration when 3D kriging? Kriging can't specifically take flow direction and magnitude as terms, but you can ...Missing: directional | Show results with:directional
[66]
Compositional cokriging for mapping the probability risk of ...
This aquifer is a highly nitrate vulnerable zone because agricultural land use is very important, which implies diffuse contamination (Chica-Olmo et al., 2014).
[67]
Developing Spatially Interpolated Surfaces and Estimating Uncertainty
In particular, kriging is a statistical model that produces both a spatial surface of predictions for the process of interest as well as the uncertainty ...
[68]
[PDF] Bayesian Kriging for Enhancing Copernicus Reanalysis Data and ...
Jul 23, 2025 · Overall, the study shows that Bayesian Kriging can generate high resolution, uncertainty maps from heterogeneous and sparse data records, making ...
[69]
Spatial measurement error and correction by spatial SIMEX in linear ...
In practice, spatial air pollution models are fit with sparse monitoring data. Hence, we examine the effects of estimation error in the Kriging model parameters ...3.1. Bias Analysis For... · 5. Simulation Study · 7. Discussion And...
[70]
Design and Analysis of Computer Experiments - Project Euclid
Jerome Sacks. William J. Welch. Toby J. Mitchell. Henry P. Wynn. "Design and Analysis of Computer Experiments." Statist. Sci. 4 (4) 409 - 423, November, 1989.
[71]
[PDF] dace.pdf - A MATLAB Kriging Toolbox - Omicron ApS
Aug 1, 2002 · This report describes the background for and use of the software package DACE. (Design and Analysis of Computer Experiments), ...
[72]
Application of Latin Hypercube Sampling Based Kriging Surrogate ...
Aug 10, 2025 · 2. Latin Hypercube Sampling Method · (1) · Correlated random variables are not considered because of · Principal Component Analysis and Nataf or ...Missing: DACE | Show results with:DACE
[73]
A robust optimization using the statistics based on kriging metamodel
The statistics such as mean and variance are obtained based on the reliable kriging model and the second-order statistical approximation method. Then, the ...
[74]
Surrogate-Based Aerodynamic Shape Optimization by Variable ...
A surrogate-based optimization algorithm for transonic airfoil design is presented. The approach replaces the direct optimization of an accurate, ...
[75]
Applications of multi-fidelity multi-output Kriging to engineering ...
May 15, 2023 · Multi-fidelity multi-output Kriging is used in engineering design optimization, including gas turbine combustor, vibrating truss, and airfoil ...
[76]
[PDF] Classical Geostatistical Methods - University of Iowa
One common stationarity assumption is that of second-order stationarity, which specifies that. Cov[e(s), e(t)] = C(st), for all s, t€ D. (3.2). In other words ...
[77]
None
### Summary of Assumptions and Challenges of Kriging/GPR from https://arxiv.org/pdf/2408.02331
[78]
Data interpolation with Kriging - Coastal Wiki
Feb 12, 2024 · The concept was developed by D. Krige (1951) and the theoretical foundation was given by G. Matheron (1969). Kriging refers ...
[79]
Advances in Kriging-Based Autonomous X-Ray Scattering ... - Nature
Jan 28, 2020 · We demonstrated the successful application of this method to exploring materials science problems using x-ray scattering measurements at a synchrotron beamline.Theory · Synthetic Test · Kriging Vs...
[80]
[PDF] Efficient kriging for real-time spatio-temporal interpolation
In this work, we formulate the kriging problem, to first reduce the computational cost to O(N3). We use an iterative solver (Saad, 2003), and further accelerate ...
[81]
Details of Ordinary Kriging
Additional problems with kriging, and with spatial estimation methods in general, are related to the necessary assumption of ergodicity of the spatial process.Missing: formulation | Show results with:formulation
[82]
[PDF] The Problem of Kriging when Estimating in a Finite Domain | CCG
While such weighting of the boundary samples is theoretically valid, we believe that it could lead to biased estimation of finite domain, especially if the data ...
[83]
Using cross validation to assess interpolation results—ArcGIS Pro
Cross validation is a leave-one-out resampling method that first uses all input points to estimate the parameters of an interpolation model.Missing: ME | Show results with:ME
[84]
Fixed Rank Kriging for Very Large Spatial Data Sets
Cressie (1993), section 3.4.5, gave a formula for the kriging predictor of Y ... N. (. 1993. ) Statistics for Spatial Data. , revised edn. New York: Wiley.
[85]
Cluster-based Kriging approximation algorithms for complexity ...
Sep 9, 2019 · In this paper, we propose a general methodology for the complexity reduction, called cluster Kriging, where the whole data set is partitioned into smaller ...
[86]
Efficient multi-scale Gaussian process regression for massive ... - GMD
Jul 31, 2020 · We design and implement a computationally efficient multi-scale Gaussian process (GP) software package, satGP, geared towards remote sensing applications.
[87]
An application of random forest plus residuals kriging approach
In this study, a hybrid approach, random forest plus residuals kriging (RFRK), was proposed to predict and map the spatial pattern of SOM for the rubber ...
[88]
[PDF] How Deep Are Deep Gaussian Processes?
Recent research has shown the potential utility of deep Gaussian processes. These deep structures are probability distributions, designed through ...
[89]
Towards annual updating of forced warming to date and constrained ...
Oct 17, 2025 · This study relies on the Kriging for Climate Change (KCC) statistical method, first introduced to constrain Global mean Surface Temperature (GST) ...Missing: universal | Show results with:universal<|separator|>
[90]
Ensemble of surrogates combining Kriging and Artificial Neural ...
In this work, two active learning approaches are proposed to combine Kriging and ANN models for reliability analysis.
[91]
DKNN: deep kriging neural network for interpretable geospatial ...
This study proposes a novel geospatial artificial intelligence (GeoAI) framework called deep kriging neural network (DKNN).Missing: 2020s | Show results with:2020s