Abstract
Many hydrological variables usually show the presence of spatial drifts. Most often they are accounted for by either universal or residual kriging and usually assuming low order polynomials. The Universal Kriging matrix (M), which includes in it the values of the polynomials at data locations (matrix F), in some cases may have a too large condition number and can even be nearly singular due to the fact that some columns are close to be linearly dependent. These problems are usually caused by a combination of pathological data locations and an inadequate choice of the coordinate system. As suggested by others, we have found that an appropriate scaling can alleviate the problem by significantly reducing the condition number of M (cond(M)). This scaling, however, does not affect the linear independence of this matrix. We show that a QR factorization of matrix F leads to a significant improvement on cond(M). An alternative to drift polynomials consists on using a set of functions derived from the eigenvectors of the variogram matrix (Γ). Although their potential usefulness as interpolating functions remains to be ascertained, they are optimal from the point of view of optimizing cond(M). In fact we are able to provide a rigorous proof for an upper bound of cond(M). Applications of the theoretical developments to hydraulic head data from an alluvial aquifer are also presented.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Preview
Unable to display preview. Download preview PDF.
References
An introduction to numerical analysis, John Wiley and Sons 692 pp. New Jersey, 573 p. Dahlquist, G. and A. Bjorck 1974 Numerical Methods. Prentice-Hall, Inc. Englewood Cliffs,, n. 5, p. 721–739
Diamond P. and Armstrong, M. 1984 Robustness of Variograms and Conditioning of Kriging Matrices. Math. GeoL, (16): 809–822.
ENRESA, 1992 Estudio Geoestadistico de la Piezometria del Acuffero de la Fâbrnca de Uranio de Andújar (Internal Report), 93 pp.
Juanes, R. and J. Samper 1996 Estudio geoestadístico de los niveles piezométricos en el acuffero de la F.U.A. (Internal Report )
ODowd, R.J., 1991 Conditioning of coefficient matrices of ordinary kriging, Math. Geol., 23 (5): 721–739.
Posa, D., 1989 Conditioning of the stationary kriging matrices for some well-known covariance models. Math. Geol., 21 (7): 755–765
Samper, J. and S. P. Neuman, 1989 Estimation of Spatial Covariance Structures by Adjoint State Maximum Likelihood Cross Validation. 1. Theory. Water Resour. Res. 25 (3): 351–362
Samper, J. and J. Carrera, 1990 GeostadLara: aplicaciones a la Hidrologfa Subterrbnea Ed. CIMNE. Barcelona Spain, 484 pp.
Samper, J. and J. Carrera, 1995 Numerical Modeling of Remedial Actions for an Uranium Mill tailings in Spain: from model construction to prediction uncertainty. In: Groundwater Quality: Remediation and Protection. Praga, IAHS Pub. N. 225, pp. 299–310.
Author information
Authors and Affiliations
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 1997 Springer Science+Business Media Dordrecht
About this paper
Cite this paper
López, C., Samper, J. (1997). Numerical Aspects of the Universal Kriging Method for Hydrological Applications. In: Soares, A., Gómez-Hernandez, J., Froidevaux, R. (eds) geoENV I — Geostatistics for Environmental Applications. Quantitative Geology and Geostatistics, vol 9. Springer, Dordrecht. https://doi.org/10.1007/978-94-017-1675-8_6
Download citation
DOI: https://doi.org/10.1007/978-94-017-1675-8_6
Publisher Name: Springer, Dordrecht
Print ISBN: 978-90-481-4861-5
Online ISBN: 978-94-017-1675-8
eBook Packages: Springer Book Archive