Nonstationary models for exploring and mapping monthly precipitation in the United Kingdom

Spatial statistical algorithms are used widely for both the exploration and mapping of environmental variables such as precipitation amount. One limitation of standard approaches to characterization and spatial interpolation is the usual assumption of stationarity. In short, spatial variation is assumed to be constant across the region of interest. Much research effort has been extended in developing approaches that allow for local variation in spatial structure. Simple moving window approaches are examples of such developments. The purpose of this paper is to apply selected methods for exploring and mapping monthly precipitation amount in the UK. Global regression, moving window regression (MWR) and geographically weighted regression (GWR) are used to explore the relationship between altitude and precipitation amount. Inverse distance weighting (IDW) is used as a basic approach to spatial prediction. Global and local variogram models are estimated and modelled to assess variation in spatial structure of precipitation amount and to inform spatial prediction using ordinary kriging (OK), kriging with an external drift (KED) and simple kriging with local means (SKlm). In the latter case, local means are derived using global regression, MWR (using ordinary least squares (OLS) and generalized least squares (GLS)) and GWR. The importance of choice of algorithm for the estimation of the variogram and for spatial prediction is assessed. The benefits of local as against global approaches and multivariate as against univariate prediction procedures are considered. It is demonstrated that use of elevation data to inform the prediction process reduces prediction errors and that there is also a small reduction in prediction errors when local variogram models are used. In this study, KED based on local variogram models provides more accurate predictions (as judged by cross-validation statistics) than any other approach. The paper concludes by considering some key issues and possible avenues for future work. Copyright © 2009 Royal Meteorological Society

[1]  Eulogio Padro-Igúzquiza VARFIT: a fortran-77 program for fitting variogram models by weighted least squares , 1999 .

[2]  Donald E. Myers,et al.  Basic Linear Geostatistics , 1998, Technometrics.

[3]  Michael F. Hutchinson,et al.  Interpolation of Rainfall Data with Thin Plate Smoothing Splines - Part II: Analysis of Topographic Dependence , 2002 .

[4]  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 .

[5]  C. Lloyd Assessing the effect of integrating elevation data into the estimation of monthly precipitation in Great Britain , 2005 .

[6]  Richardo A. Olea,et al.  Semivariogram modeling by weighted least squares , 1996 .

[7]  Christel Prudhomme,et al.  MAPPING EXTREME RAINFALL IN A MOUNTAINOUS REGION USING GEOSTATISTICAL TECHNIQUES: A CASE STUDY IN SCOTLAND , 1999 .

[8]  Penelope A. Hancock,et al.  Spatial interpolation of large climate data sets using bivariate thin plate smoothing splines , 2006, Environ. Model. Softw..

[9]  Noel A. C. Cressie,et al.  Statistics for Spatial Data: Cressie/Statistics , 1993 .

[10]  Carol A. Gotway,et al.  Statistical Methods for Spatial Data Analysis , 2004 .

[11]  Edzer Pebesma,et al.  GSTAT: a program for geostatistical modelling, prediction and simulation , 1998 .

[12]  C. Lloyd Multivariate Interpolation of Monthly Precipitation Amount in the United Kingdom , 2010 .

[13]  A. Flint,et al.  Precipitation Estimation in Mountainous Terrain Using Multivariate Geostatistics. Part II: Isohyetal Maps , 1990 .

[14]  P. Atkinson,et al.  Deriving DSMs from LiDAR data with kriging , 2002 .

[15]  Eulogio Pardo-Igúzquiza,et al.  An automatic moving window approach for mapping meteorological data , 2005 .

[16]  R. Reese Geostatistics for Environmental Scientists , 2001 .

[17]  P. Kitanidis Introduction to Geostatistics: Applications in Hydrogeology , 1997 .

[18]  Geostatistical Interpolation of Rainfall in Mountainous Areas , 2001 .

[19]  Guy Delrieu,et al.  Rain Measurement by Raingage-Radar Combination: A Geostatistical Approach , 1988 .

[20]  G. McCabe,et al.  Precipitation interpolation in mountainous regions using multiple linear regression , 1998 .

[21]  J. Jaime Gómez-Hernández,et al.  Incorporating Information from a Digital Elevation Model for Improving the Areal Estimation of Rainfall , 2001 .

[22]  Michel Gevers,et al.  OPTIMAL ESTIMATION OF THE AVERAGE AREAL RAINFALL AND OPTIMAL SELECTION OF RAIN-GAUGE LOCATIONS , 1984 .

[23]  J. Salas,et al.  A COMPARATIVE ANALYSIS OF TECHNIQUES FOR SPATIAL INTERPOLATION OF PRECIPITATION , 1985 .

[24]  C. Lloyd Local Models for Spatial Analysis , 2006 .

[25]  Noel A Cressie,et al.  Statistics for Spatial Data, Revised Edition. , 1994 .

[26]  A. Flint,et al.  Precipitation Estimation in Mountainous Terrain Using Multivariate Geostatistics. Part I: Structural Analysis , 1992 .

[27]  T. C. Haas,et al.  Lognormal and Moving Window Methods of Estimating Acid Deposition , 1990 .

[28]  Lubos Mitas,et al.  Multivariate Interpolation of Precipitation Using Regularized Spline with Tension , 2002, Trans. GIS.

[29]  Michael F. Hutchinson,et al.  Interpolation of Rainfall Data with Thin Plate Smoothing Splines - Part I: Two Dimensional Smoothing of Data with Short Range Correlation , 1998 .

[30]  P. Burrough,et al.  Principles of geographical information systems , 1998 .

[31]  Peter Guttorp,et al.  Advances in Modeling and Inference for Environmental Processes with Nonstationary Spatial Covariance , 2001 .

[32]  C. Daly,et al.  A Statistical-Topographic Model for Mapping Climatological Precipitation over Mountainous Terrain , 1994 .

[33]  Hans Wackernagel,et al.  Multivariate Geostatistics: An Introduction with Applications , 1996 .

[34]  E. Pardo‐Igúzquiza COMPARISON OF GEOSTATISTICAL METHODS FOR ESTIMATING THE AREAL AVERAGE CLIMATOLOGICAL RAINFALL MEAN USING DATA ON PRECIPITATION AND TOPOGRAPHY , 1998 .

[35]  Edzer J. Pebesma,et al.  Multivariable geostatistics in S: the gstat package , 2004, Comput. Geosci..

[36]  Chris Brunsdon,et al.  Spatial variations in the average rainfall–altitude relationship in Great Britain: an approach using geographically weighted regression , 2001 .

[37]  Timothy C. Coburn,et al.  Geostatistics for Natural Resources Evaluation , 2000, Technometrics.

[38]  P. Goovaerts Geostatistical approaches for incorporating elevation into the spatial interpolation of rainfall , 2000 .

[39]  Eulogio Pardo-Igúzquiza,et al.  MLREML: a computer program for the inference of spatial covariance parameters by maximum likelihood and restricted maximum liklihood , 1997 .

[40]  Gerard B. M. Heuvelink,et al.  About regression-kriging: From equations to case studies , 2007, Comput. Geosci..

[41]  A. McBratney,et al.  Choosing functions for semi‐variograms of soil properties and fitting them to sampling estimates , 1986 .

[42]  M. F. Hutchinson,et al.  Interpolating Mean Rainfall Using Thin Plate Smoothing Splines , 1995, Int. J. Geogr. Inf. Sci..