Introduction to Choosing a Kriging Plan

Clayton Deutsch1 and Jared Deutsch2

1 University of Alberta
2 University of Alberta

Learning Objectives
• Classify estimates according to their purpose: interim estimation, final decisions, visual‐
ization, and probabilistic prediction.
• Know best practice in choosing a fit‐for‐purpose kriging plan.
• Understand stationarity, kriging performance measures, and other concerns when choos‐
ing a kriging search plan.

1 Introduction
This lesson relates primarily to mining where estimation is still widely practiced for resource and
reserve estimation. Kriging is the primary technique for the estimation of grades. Kriging is a linear
unbiased estimator that minimizes the estimation variance using a site‐specific variogram model of
spatial variability accounting for anisotropy and other spatial features (Journel & Huijbregts, 1978).
The number of data and locations of data used to inform the estimate compose the kriging search
plan. The kriging search plan is typically restricted due to concerns about computational time, sta‐
tionarity, conditional bias, and histogram reproduction. These restrictions and concerns are the
subject of this lesson.
There is no doubt that the prerequisites to kriging are more important than the details of the
kriging plan. These prerequisites include ensuring data quality, compositing, managing outliers,
subdividing the data into reasonable subsets (domains), establishing an appropriate coordinate sys‐
tem for estimation, choosing the correct block size for estimation, and managing large scale trends
and contacts between estimation domains. Nevertheless, details of the kriging plan are important
provided these prerequisites have been met in a reasonable fashion.
The choice of a kriging search plan is dictated by the purpose of the estimate; there is no uni‐
versal best kriging plan. Best practice in choosing a kriging search plan depends on the estimate
purpose. This lesson addresses the selection of a kriging search plan for four types of estimates,
and discusses concerns of stationarity, kriging measures, and other common practices.

2 Purpose of the Estimate

We classify the purpose of estimates into four categories: 1) interim estimates, 2) final estimates, 3)
visualization and trend models, and 4) probabilistic predictions. Each of these estimation purposes
has a different set of criteria for choosing the best kriging search plan.

Interim Estimates
Interim estimates are estimates awaiting more data. Early in the life cycle of a mineral deposit there
are relatively few data available for estimation. At the time of mining, many more data are available
in the form of blast holes and drill holes. This disparity in available data is termed the informa‐
tion effect (Rossi & Deutsch, 2013). One goal of interim estimates is to provide close estimates of
ore tonnage, ore grade and waste tonnage within reasonably large production volumes (sometimes
called mining panels). A too‐large search will result in over smoothing of the grades and inaccurate

GeostatisticsLessons.com ©2015 C. Deutsch and J. Deutsch 1

estimation of tonnage and grade. The search could be restricted to improve the estimate of ton‐
nage and grade; however, this will come at a cost of conditional bias. This cost is not severe since
interim estimates are not final estimates; these may be used in the mine plan and for technical and
economic evaluation; however, additional information will be available prior to mining and final
grade control decisions.

Final Estimates
Final estimates are an interesting time in a mining context; the classification of material as ore or
waste (grade control) will have a direct and final economic impact on the mine. The emphasis is on
the minimization of conditional bias, minimization of mean squared error and minimization of Type
I and Type II errors (Isaaks, 2005). Final estimates are constructed with the goal of making the best
possible estimate according to our mean squared error criterion. Final estimates are made with all
data that will be available; no additional data is expected prior to use of the estimate.
Interim estimates in an underground mining context may be considered as near‐final estimates
depending on the flexibility of the mining method. If there is little future flexibility, then concerns
of conditional bias and Type I/II errors are important and estimates could be considered as final.

Visualization and Trend Model Estimates

Visualization and trend models include all models where the goal is a smooth, interpretable map
of the properties of interest. The purpose is to understand the entire spatial distribution of the
property being estimated, not discrete estimates one at a time as in resource and reserve estima‐
tion. Visualization estimates are commonly used with geologic mapping to evaluate our conceptual
geologic model, and draw conclusions on the presence and nature of trends in the data. The require‐
ment for these models is that they be free of numerical artifacts to permit reliable interpretation.

Probabilistic Prediction Estimates

Probabilistic predictions, including multigaussian kriging (Verly, 1984) and indicator kriging (Journel,
1983) within an estimation or simulation context constitute the fourth type of estimate. In this case,
the kriging equations are used to construct a conditional probability distribution at an unsampled
location in the presence of correlated data. For multigaussian kriging, the kriging equations are
typically referred to as the normal equations. The normal equation formalism is used throughout
Gaussian and probabilistic geostatistics. The goal of probabilistic predictions is to infer the most
accurate set of probabilities, free from any conditional bias or artifacts. Smoothing is not a concern
since the predicted distribution represents the uncertainty and variability.

3 Best Practice
As discussed earlier, best practice for all estimate types is to only apply estimation within a decided
stationary domain. All prerequisites including composite length selection, block size selection, and
stationary domain selection must be considered prior to estimation. The decision of stationarity
has typically been made long before the application of kriging; this decision is made as soon as data
are grouped together for analysis.

Best Practice: Interim Estimates

Best practice for interim estimates is to use a restricted search with ordinary kriging to compute
estimates that closely matches the anticipated grade tonnage curve (Parker, 1979). Ordinary krig‐
ing provides direct control on the smoothing; a small number of data in the search will lead to a
histogram of estimates with little smoothing. The spatial distribution of estimates will always be

GeostatisticsLessons.com ©2015 C. Deutsch and J. Deutsch 2

Figure 1: Effect of restricting the maximum number of search data used in kriging estimates com‐
pared to the target histogram. The maximum number of search data used for each histogram is an‐
notated (5, 10, 40). The desired histogram corresponds to the dispersion variance of blocks within
the domain.

smooth, but the histogram of estimates can be controlled. The search should be restricted by lim‐
iting the maximum number of data used to compute the kriging estimate. A known histogram of
values at the scale of interest (such as the selective mining unit size) is targeted. This histogram
may be calibrated from production data, if any are available, and by volume variance relations. The
impact of restricting the number of search data is sketched in the following figure.
There are additional considerations when restricting the search for ordinary kriging. The number
of data used from each drill hole may be restricted. This is common when making estimates in a 3D
domain with vertical drilling and a relatively short composite length; much more data is collected
in the vertical direction compared to the horizontal plane. A maximum per drill hole avoids using
data from only a single drill hole when making an estimate with very few data. The search range
may also be restricted; the range of the variogram is a common choice for the search range, but in
practice restricting the number of data used in an estimate is a more reliable method of restricting
the search.
The consequence of restricting the number of data used in the kriging estimate is conditional
bias. High grade areas are overestimated because local high grade values are given more weight.
Low grade areas are correspondingly underestimated because of additional weight given to low
grade values. There will be no global bias and there should be no bias in ore tonnes and ore grade
at the targeted cutoff. Conditional bias is not a problem for interim estimates as they will be updated
before final decisions are made; however, care should be taken not to use this type of estimate for
final decisions.

Best Practice: Final Estimates

The kriging plan for final estimates is chosen to minimize mean squared error and conditional bias.
Best practice for making a final estimate is to use ordinary kriging with a large number of search data
to minimize the estimation error. Simple kriging is typically not as effective due to an overly strong
dependence on the global mean. This is shown in the following figure plotting the mean squared
error from cross validation as a function of the maximum number of search data used in kriging

GeostatisticsLessons.com ©2015 C. Deutsch and J. Deutsch 3

Figure 2: Mean squared error in cross validation as a function of the number of search data included
for a copper porphyry deposit, using data from (Deutsch, Szymanski, & Deutsch, 2014).

a copper porphyry deposit. Using additional search data results in only a marginal improvement
past a certain amount of data, but does not result in a worse estimate. Specific concerns about
stationarity and performance measures are addressed later in this lesson.
The specific number of search data chosen when making a final estimate is a function of the
computational time required for the estimate. Choosing 40 search data in a 3D estimation context,
or 24 in a 2D context, is often a reasonable number to balance the computational time and desire
for the best estimate. The number of data used weakly depends on the support of the data. More
data of small support could be considered.
The search range should be set at the variogram range or even larger to include the required
number of data. Often, production data are reasonably closely spaced and the effective search
range that will get the required number of data is reasonably small. Limiting the kriging search
range to an arbitrary value that excludes data from the estimate is not best practice. The search
radius should consider the variogram anisotropy.

Best Practice: Visualization and Trend Model Estimates

Best practice for visualization and trend model estimates is similar to final estimates. The goal of
a smoothly varying interpretable map is best met by kriging with a large search. Typically simple
kriging is used to reduce the impact of local means in sparsely sampled and peripheral areas which
could influence our interpretation. Ideally global kriging is used, but a very large number of search
data may be used if there are too many data for global kriging (about 10,000). Limiting the number of
search data results in “shadows” and numerical artifacts, as shown in the following figure. Artifacts
are visible particularly in the upper portion of the map which is undesirable for interpretation.

Best Practice: Probabilistic Prediction Estimates

Probabilistic prediction estimates, including multigaussian and indicator kriging, are made with the
goal of inferring the best probability distribution. Best practice is to use a large search to include all
data that should influence the predicted probability distribution. Simple kriging is the correct kriging
type for prediction in a multivariate Gaussian context. This is true for both multigaussian kriging for

GeostatisticsLessons.com ©2015 C. Deutsch and J. Deutsch 4

Figure 3: Smoothly varying visualization map (left) and a restricted search map (right) with visible
shadows and artifacts.

local uncertainty characterization or within simulation. Ordinary kriging is more often applied for in‐
dicator kriging. Simple kriging is required to ensure covariance reproduction; however, stationarity
concerns are often considered more important and there are no strong theoretical requirements as
with the multivariate Gaussian distribution. Concerns such as over‐smoothing for interim estimates
are not applicable in a probabilistic estimation context where the parameters of probability distri‐
butions are being estimated. Limiting the number of search data results in suboptimal estimates of
the probability distributions.

4 Considerations
As kriging has been widely applied since the 1960s, a large number of specific considerations have
been discussed by geostatistical practitioners. Some of these considerations, including stationar‐
ity, kriging performance measures, and multipass kriging are discussed here as they relate to best
practices for estimation.

What About Stationarity?

One might think that the number of search data should be reduced to account for local fluctuations
in stationarity. The logic is that a limited search will produce a better estimate of the local mean
for estimation. This has not been observed in practice; using more search data does not result in
a worse (higher mean squared error) estimate. The apparent improvement due to locally adapting
to variations in the mean is not backed up in jackknife, cross validation, or when checking with
production data. Final estimates benefit from a reasonably large number of data despite concerns
about smoothing.

GeostatisticsLessons.com ©2015 C. Deutsch and J. Deutsch 5

What About Kriging Performance Measures?
Measures of kriging performance are variably popular. The slope of regression, kriging efficiencies,
the prevalence of negative weights, and weight to the simple kriging mean are often cited as impor‐
tant measures and used in the selection of kriging plans (Vann, Jackson, & Bertoli, 2003).
The slope of regression of estimates against true values is used as an approximation of the con‐
ditional bias of the kriging estimate. Simple kriging, which is an unbiased estimator, has a slope of
regression of exactly one. In ordinary kriging, where a Lagrange multiplier is used, the slope of re‐
gression is typically less than one indicating a conditional bias. Increasing deviation from a slope of
one indicates an increased conditional bias. Although the slope of regression is often documented;
it is of little utility when choosing a kriging search plan. In an interim estimation context, the most
important parameter is matching the desired histogram. In a final, visualization, or probabilistic
estimation context, a very large number of search data is used to minimize the mean squared error
and conditional bias is essentially zero. Ideally the conditional bias is eliminated for a final estimate,
but this may not always be possible (Isaaks, 2005).
Kriging efficiency (Krige, 1997) is a measure of the proportion of variance explained. The related
measure of the statistical efficiency of kriging (Deutsch, Szymanski, & Deutsch, 2014) provides a
measure of the effectiveness of an estimator relative to the optimal global simple kriging estimator.
These measures may be mapped over the domain of interest. Application of these efficiency mea‐
sures for decision making is challenging due to large differences between mineral deposits (Rossi
& Deutsch, 2013). When mapped, the kriging efficiency is relatively low in sparsely sampled areas,
and high in densely sampled areas. It is evident that estimates are worse when data are widely
spaced. The kriging variance is a robust measure of data spacing and may be used in place of the
kriging efficiency.

What About Negative Kriging Weights?

Negative kriging weights are applied to data that are screened behind other data more highly corre‐
lated to the location being estimated. Most negative kriging weights improve estimation because
local trends in the data are more accurately interpolated. There are concerns, however, when the
negative weights become extreme and when they cause negative grade estimates. Clearly, nega‐
tive grade estimates should be reset to zero. A significant proportion of negative estimates in an
estimated model implies that the variogram may have been modeled too continuously.
Closely related to negative kriging weights is the so‐called string effect where samples at the
ends of a string of data receive more weight than data within the string. These strings of data are
caused by the search for data along a drill hole. Cross validation should identify this problem with
some large errors. Limiting the number of data used in the kriging plan per drill hole will mitigate
this problem.

What About…?
There are some other considerations. In general, except for interim estimates, the practitioner
should err on the side of using more data in the kriging. The kriging equations will sort out the
optimal weight for the data; data far from the unsampled location or screened by closer data will
not receive any significant weight.
The use of multiple search passes is aimed at the important problem of classification and at
restricting estimation to use nearby data in areas of more data. This practice introduces artifacts
where one search strategy ceases to be satisfied and another is considered. Choosing a fairly large
search and restricting the maximum number of data accounts for varying data density. The esti‐
mates may be better classified by data spacing calibrated to simulation‐based uncertainty measures.
There may be a need to carefully clip a kriged model to avoid estimating too far from the data,
below the deepest drill holes or at the margins. There may be a need to account for soft bound‐

GeostatisticsLessons.com ©2015 C. Deutsch and J. Deutsch 6

aries and many other practical details. These concerns warrant additional discussion, and will be
addressed in future lessons; this lesson forms an introduction to choosing a kriging plan.

5 References
Deutsch, J., Szymanski, J., & Deutsch, C. (2014). Checks and measures of performance for kriging
estimates. Journal of the Southern African Institute of Mining and Metallurgy, 114(3), 223–
223.
Isaaks, E. (2005). The kriging oxymoron: A conditionally unbiased and accurate predictor (2nd
edition). In O. Leuangthong & C. V. Deutsch (Eds.), Geostatistics banff 2004 (Vol. 14, pp.
363–374). Springer Netherlands.
Journel, A. G. (1983). Nonparametric estimation of spatial distributions. Journal of the Interna‐
tional Association for Mathematical Geology, 15(3), 445–468. Journal Article.
Journel, A. G., & Huijbregts, Ch. J. (1978). Mining geostatistics (p. 600). book, Blackburn Press.
Krige, D. (1997). A practical analysis of the effects of spatial structure and of data available and
accessed, on conditional biases in ordinary kriging. Geostatistics Wollongong, 96, 799–810.
Parker, H. (1979). The volume variance relationship: A useful tool for mine planning. Engineering
and Mining Journal, 180, 106–126.
Rossi, M. E., & Deutsch, C. V. (2013). Mineral resource estimation. Springer Science & Business
Media.
Vann, J., Jackson, S., & Bertoli, O. (2003). Quantitative kriging neighbourhood analysis for the
mining geologist–a description of the method with worked case examples. In 5th interna‐
tional mining geology conference (Vol. 8, pp. 215–223). Bedigo: AUSIMM.
Verly, G. (1984). The block distribution given a point multivariate normal distribution. In G. Verly,
M. David, A. G. Journel, & A. Marechal (Eds.), Geostatistics for natural resources characteri‐
zation (pp. 495–515). Springer Netherlands.

Citation
Deutsch, C. V., & Deutsch, J. L. (2015). Introduction to Choosing a Kriging Plan. In J. L. Deutsch (Ed.),
Geostatistics Lessons. Retrieved from http://geostatisticslessons.com/lessons/introkrigingplan

GeostatisticsLessons.com ©2015 C. Deutsch and J. Deutsch 7