The KRIGE2D Procedure

The Nugget Effect

For all the semivariogram models considered previously, the following property holds:

gamma Subscript z Baseline left-parenthesis 0 right-parenthesis equals limit Underscript h down-arrow 0 Endscripts gamma Subscript z Baseline left-parenthesis h right-parenthesis equals 0

However, a plot of the experimental semivariogram might indicate a discontinuity at h = 0; that is, gamma Subscript z Baseline left-parenthesis h right-parenthesis right-arrow c Subscript n Baseline greater-than 0 as h right-arrow 0, while gamma Subscript z Baseline left-parenthesis 0 right-parenthesis equals 0. The quantity c Subscript n is called the nugget effect; this term is from mining geostatistics where nuggets literally exist, and it represents variations at a much smaller scale than any of the measured pairwise distances—that is, at distances h much-less-than h Subscript m i n, where

h Subscript m i n Baseline equals min Underscript i comma j Endscripts h Subscript i j Baseline equals min Underscript i comma j Endscripts bar bold-italic s Subscript i Baseline minus bold-italic s Subscript j bar

Nonzero nugget effects have been associated with conceptual and theoretical difficulties; see Cressie (1993, section 2.3.1) and Christakos (1992, section 7.4.3) for details. There is no practical difficulty, however; you simply visually extrapolate the experimental semivariogram as h right-arrow 0. The importance of availability of data at small lag distances is again illustrated.

As an example, an exponential semivariogram with a nugget effect c Subscript n has the form

gamma Subscript z Baseline left-parenthesis h right-parenthesis equals c Subscript n Baseline plus sigma 0 squared left-bracket 1 minus exp left-parenthesis minus StartFraction h Over a 0 EndFraction right-parenthesis right-bracket comma h greater-than 0

and

gamma Subscript z Baseline left-parenthesis 0 right-parenthesis equals 0

where the factor sigma 0 squared is called the partial sill and the sill c 0 equals c Subscript n Baseline plus sigma 0 squared.

This is illustrated in Figure 11 for the parameters a 0 equals 1, sigma 0 squared equals 4, and nugget effect c Subscript n Baseline equals 1.5.

You can specify the nugget effect in PROC KRIGE2D with the NUGGET= option in the MODEL statement. It is a separate, additive term independent of direction; that is, it is isotropic. The way to approximate an anisotropic nugget effect is described in the following section.

Figure 11: Exponential Semivariogram Model with a Nugget Effect c Subscript n Baseline equals 1.5

Exponential Semivariogram Model with a Nugget Effect cn=1.5


Last updated: December 09, 2022