Orthonormal residuals in geostatistics: Model criticism and parameter estimation

In linear geostatistics, models for the mean function (drift) and the variogram or generalized covariance function are selected on the basis of the modeler's understanding of the phenomenon studied as well as data. One can seldom be assured that the most appropriate model has been selected; however, analysis of residuals is helpful in diagnosing whether some important characteristic of the data has been neglected and, ultimately, in providing a reasonable degree of assurance that the selected model is consistent with the available information. The orthonormal residuals presented in this work are kriging errors constructed so that, when the correct model is used, they are uncorrelated and have zero mean and unit variance. It is suggested that testing of orthonormal residuals is a practical way for evaluating the agreement of the model with the data and for diagnosing model deficiencies. Their advantages over the usually employed standardized residuals are discussed. A set of tests are presented. Orthonormal residuals can also be useful in the estimation of the covariance (or variogram) parameters for a model that is considered correct.