Estimating Variogram Uncertainty

The variogram is central to any geostatistical survey, but the precision of a variogram estimated from sample data by the method of moments is unknown. It is important to be able to quantify variogram uncertainty to ensure that the variogram estimate is sufficiently accurate for kriging. In previous studies theoretical expressions have been derived to approximate uncertainty in both estimates of the experimental variogram and fitted variogram models. These expressions rely upon various statistical assumptions about the data and are largely untested. They express variogram uncertainty as functions of the sampling positions and the underlying variogram. Thus the expressions can be used to design efficient sampling schemes for estimating a particular variogram. Extensive simulation tests show that for a Gaussian variable with a known variogram, the expression for the uncertainty of the experimental variogram estimate is accurate. In practice however, the variogram of the variable is unknown and the fitted variogram model must be used instead. For sampling schemes of 100 points or more this has only a small effect on the accuracy of the uncertainty estimate. The theoretical expressions for the uncertainty of fitted variogram models generally overestimate the precision of fitted parameters. The uncertainty of the fitted parameters can be determined more accurately by simulating multiple experimental variograms and fitting variogram models to these. The tests emphasize the importance of distinguishing between the variogram of the field being surveyed and the variogram of the random process which generated the field. These variograms are not necessarily identical. Most studies of variogram uncertainty describe the uncertainty associated with the variogram of the random process. Generally however, it is the variogram of the field being surveyed which is of interest. For intensive sampling schemes, estimates of the field variogram are significantly more precise than estimates of the random process variogram. It is important, when designing efficient sampling schemes or fitting variogram models, that the appropriate expression for variogram uncertainty is applied.