Empirical semivariogram and covariance functions

The semivariogram and covariance functions are theoretical quantities that you cannot observe, so you estimate them from your data, using what is called the empirical semivariogram and empirical covariance functions. Often, you can gain insight into the quantities by looking at the way they are estimated. Suppose you take all pairs of data that are a similar distance and direction from each other, such as those connected by the blue lines in the following figure.

Empirical semivariogram

For all the pairs of locations s i and s j that are joined by lines, compute

average[(z(s i ) - z(s j ))2] ,

where z(s i ) is the measured value at location s i . If all the pairs of locations s i and s j are close to each other, it's expected that z(s i ) and z(s j ) will be similar in value, so when you take the differences and square them, the average should be small. As s i and s j get farther apart, it's expected that their values will become more dissimilar, so when you take their differences and square them, the average will get larger.

In the covariance function, for all the pairs of locations s i and s j that are joined by lines, the software computes

average [(Z(s i )-Element of the covariance function)(Z(s j )-Element of the covariance function)],

where z(s i ) is the measured value at location s i and Element of the covariance function is the mean of all of the data. Now, if all the pairs s i and s j are close to each other, it's expected either that both z(s i ) and z(s j ) will be above the mean Element of the covariance function or both will be below the mean. Either way, their product is positive, so when you average all of the products, you expect a positive value. If s i and s j are far apart, it's expected that about half the time the products will be negative and half the time they will be positive, so you expect their average to be near zero.

In Geostatistical Analyst, for all pairs that have a similar distance and angle, the average values calculated above are plotted on a semivariogram or covariance surface. For example, here is an empirical semivariogram surface:

Empirical semivariogram surface

Learn more about empirical semivariograms


7/11/2012