GSTAT: a program for geostatistical modelling, prediction and simulation

Abstract Gstat is a computer program for variogram modelling, and geostatistical prediction and simulation. It provides a generic implementation of the multivariable linear model with trends modelled as a linear function of coordinate polynomials or of user–defined base functions, and independent or dependent, geostatistically modelled, residuals. Simulation in gstat comprises conditional or unconditional (multi-) Gaussian sequential simulation of point values or block averages, or (multi-) indicator sequential simulation. Besides many of the popular options found in other geostatistical software packages, gstat offers the unique combination of (i) an interactive user interface for modelling variograms and generalized covariances (residual variograms), that uses the device-independent plotting program gnuplot for graphical display, (ii) support for several ascii and binary data and map file formats for input and output, (iii) a concise, intuitive and flexible command language, (iv) user customization of program defaults, (v) no built-in limits, and (vi) free, portable ansi-c source code. This paper describes the class of problems gstat can solve, and addresses aspects of efficiency and implementation, managing geostatistical projects, and relevant technical details.

[1]  Tony Mason,et al.  Lex & Yacc , 1992 .

[2]  Peter K. Kitanidis,et al.  Orthonormal residuals in geostatistics: Model criticism and parameter estimation , 1991 .

[3]  Donald E. Myers Vector Conditional Simulation , 1989 .

[4]  Derek Karssenberg,et al.  Integrating dynamic environmental models in GIS: The development of a Dynamic Modelling language , 1996, Trans. GIS.

[5]  M. Stein Large sample properties of simulations using latin hypercube sampling , 1987 .

[6]  J. Carr,et al.  Revisiting the accurate calculation of block-sample covariances using Gauss quadrature , 1993 .

[7]  Norman R. Draper,et al.  Applied regression analysis (2. ed.) , 1981, Wiley series in probability and mathematical statistics.

[8]  X. F. Zhang,et al.  On the weighted least-squares method for fitting a semivariogram model , 1995 .

[9]  D. Myers Pseudo-cross variograms, positive-definiteness, and cokriging , 1991 .

[10]  Jeffrey H. Kingston The design and implementation of the lout document formatting language , 1993, Softw. Pract. Exp..

[11]  David E. Stewart,et al.  Meschach : matrix computations in C , 1994 .

[12]  K. Mardia,et al.  Maximum likelihood estimation of models for residual covariance in spatial regression , 1984 .

[13]  A. Journel Fundamentals of geostatistics in five lessons , 1991 .

[14]  Mark E. Johnson Multivariate Statistical Simulation: Johnson/Multivariate , 1987 .

[15]  Gene H. Golub,et al.  Matrix computations , 1983 .

[16]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

[17]  Ronald Christensen,et al.  Quadratic covariance estimation and equivalence of predictions , 1993 .

[18]  Edzer Pebesma,et al.  Mapping groundwater quality in the Netherlands , 1997 .

[19]  Margaret Armstrong,et al.  Problems with universal kriging , 1984 .

[20]  Michael Edward Hohn,et al.  An Introduction to Applied Geostatistics: by Edward H. Isaaks and R. Mohan Srivastava, 1989, Oxford University Press, New York, 561 p., ISBN 0-19-505012-6, ISBN 0-19-505013-4 (paperback), $55.00 cloth, $35.00 paper (US) , 1991 .

[21]  P. Burrough Principles of Geographical Information Systems for Land Resources Assessment , 1986 .

[22]  P. Kitanidis Generalized covariance functions in estimation , 1993 .

[23]  Peter K. Kitanidis,et al.  Statistical estimation of polynomial generalized covariance functions and hydrologic applications , 1983 .

[24]  Ronald Christensen Plane Answers to Complex Questions , 1987 .

[25]  M. E. Johnson,et al.  Multivariate Statistical Simulation , 1988 .

[26]  N. Draper,et al.  Applied Regression Analysis , 1966 .

[27]  J. Gómez-Hernández,et al.  Joint Sequential Simulation of MultiGaussian Fields , 1993 .

[28]  E. J. Englund,et al.  Geo-EAS (Geostatistical Environmental Assessment Software) user's guide , 1988 .

[29]  Carol A. Gotway,et al.  The use of conditional simulation in nuclear-waste-site performance assessment , 1994 .

[30]  John R. Levine,et al.  Lex & yacc, 2nd edition , 1992 .

[31]  Noel A Cressie,et al.  Statistics for Spatial Data. , 1992 .

[32]  Alfred V. Aho,et al.  Compilers: Principles, Techniques, and Tools , 1986, Addison-Wesley series in computer science / World student series edition.

[33]  Noel A Cressie,et al.  Multivariable spatial prediction , 1993 .

[34]  Yvan Pannatier,et al.  Variowin: Software for Spatial Data Analysis in 2D , 1996 .

[35]  D. Myers Matrix formulation of co-kriging , 1982 .

[36]  N. Cressie Fitting variogram models by weighted least squares , 1985 .

[37]  R. Christensen Linear Models for Multivariate, Time Series, and Spatial Data , 1997 .