Spatial interpolation of monthly climate data for Finland: comparing the performance of kriging and generalized additive models

The Finnish Meteorological Institute has calculated statistics for the new reference period of 1981–2010. During this project, the grid size has been reduced from 10 to 1 km, the evaluation of the interpolation has been improved, and comparisons between different methods has been performed. The climate variables of interest were monthly mean temperature and mean precipitation, for which the spatial variability was explained using auxiliary information: mean elevation, sea percentage, and lake percentage. We compared three methods for spatial prediction: kriging with external drift (KED), generalized additive models (GAM), and GAM combined with residual kriging (GK). Every interpolation file now has attached statistical key figures describing the bias and the normality of the prediction error. According to the cross-validation results, GAM was the best method for predicting mean temperatures, with only very small differences relative to the other methods. For mean precipitation, KED produced the most accurate predictions, followed by GK. In both cases, there was notable seasonal variation in the statistical skill scores. For the new reference period and future interpolations, KED was chosen as the primary method due to its robustness and accuracy.

[1]  Eric F. Wood,et al.  Correction of Global Precipitation Products for Orographic Effects , 2006 .

[2]  Gudmund Høst,et al.  Kriging by local polynomials , 1999 .

[3]  M. R. Holdaway,et al.  Spatial modeling and interpolation of monthly temperature using kriging , 2006 .

[4]  Tomislav Hengl,et al.  Comparison of kriging with external drift and regression-kriging. Technical report , 2003 .

[5]  Trevor Hastie,et al.  Generalized linear and generalized additive models in studies of species distributions: setting the scene , 2002 .

[6]  A. Vajda,et al.  The influence of natural conditions on the spatial variation of climate in Lapland, northern Finland , 2003 .

[7]  L. Kullman.  A Richer, Greener and Smaller Alpine World: Review and Projection of Warming-Induced Plant Cover Change in the Swedish Scandes , 2010, AMBIO.

[8]  M. Perry,et al.  The generation of monthly gridded datasets for a range of climatic variables over the UK , 2005 .

[9]  H. Tuomenvirta,et al.  Climate change projections for Finland during the 21st century , 2004 .

[10]  M. Hutchinson,et al.  A comparison of two statistical methods for spatial interpolation of Canadian monthly mean climate data , 2000 .

[11]  Guangqing Chi,et al.  Applied Spatial Data Analysis with R , 2015 .

[12]  N. Cressie.  The origins of kriging , 1990 .

[13]  C. Willmott.  Some Comments on the Evaluation of Model Performance , 1982 .

[14]  Christopher Daly,et al.  Guidelines for assessing the suitability of spatial climate data sets , 2006 .

[15]  S. Vicente‐Serrano,et al.  Comparative analysis of interpolation methods in the middle Ebro Valley (Spain): application to annual precipitation and temperature , 2003 .

[16]  Miska Luoto,et al.  Potential effect of climate change on the distribution of palsa mires in subarctic Fennoscandia , 2006 .

[17]  Holdaway.  Spatial modeling and interpolation of monthly temperature using kriging , 1996 .

[18]  David Ginsbourger,et al.  A Note on the Choice and the Estimation of Kriging Models for the Analysis of Computer Experiments , 2007 .

[19]  G. Heuvelink,et al.  Spatio-temporal prediction of daily temperatures using time-series of MODIS LST images , 2013, Theoretical and Applied Climatology.

[20]  Geoffrey M. Laslett,et al.  Kriging and Splines: An Empirical Comparison of their Predictive Performance in Some Applications , 1994 .

[21]  S. Wood,et al.  GAMs with integrated model selection using penalized regression splines and applications to environmental modelling , 2002 .

[22]  S. Wood.  Generalized Additive Models: An Introduction with R , 2006 .

[23]  S. Wood.  Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models , 2011 .

[24]  T. Yee,et al.  Generalized additive models in plant ecology , 1991 .

[25]  Gwo-Fong Lin,et al.  A spatial interpolation method based on radial basis function networks incorporating a semivariogram model , 2004 .

[26]  S. Wood.  Stable and Efficient Multiple Smoothing Parameter Estimation for Generalized Additive Models , 2004 .

[27]  David Ginsbourger,et al.  A Note on the Choice and the Estimation of Kriging Models for the Analysis of Computer Experiments , 2007 .

[28]  G. Matheron.  Principles of geostatistics , 1963 .

[29]  P. McCullagh,et al.  Generalized Linear Models , 1992 .

[30]  Edzer Pebesma,et al.  Applied Spatial Data Analysis with R. Springer , 2008 .

[31]  R. Olea.  Geostatistics for Natural Resources Evaluation By Pierre Goovaerts, Oxford University Press, Applied Geostatistics Series, 1997, 483 p., hardcover, $65 (U.S.), ISBN 0-19-511538-4 , 1999 .

[32]  T. Hastie,et al.  Comparative performance of generalized additive models and multivariate adaptive regression splines for statistical modelling of species distributions , 2006 .

[33]  M. Wand,et al.  Geoadditive models , 2003 .

[34]  Marc Voltz,et al.  A comparison of kriging, cubic splines and classification for predicting soil properties from sample information , 1990 .

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

[36]  Ten-Lin Liu,et al.  Interpolating Soil Properties Using Kriging Combined with Categorical Information of Soil Maps , 2006 .

[37]  R. Tibshirani,et al.  Generalized Additive Models , 1991 .

[38]  M. Luoto,et al.  Novel theoretical insights into geomorphic process–environment relationships using simulated response curves , 2011 .

[39]  Francis W. Zwiers,et al.  Detecting and attributing external influences on the climate system: a review of recent advances , 2005 .

[40]  P. McCullagh,et al.  Monograph on Statistics and Applied Probability , 1989 .

[41]  G. Heuvelink,et al.  Optimization of sample patterns for universal kriging of environmental variables , 2007 .

[42]  P. Goovaerts.  Ordinary Cokriging Revisited , 1998 .

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

[44]  P. Legendre.  Spatial Autocorrelation: Trouble or New Paradigm? , 1993 .

[45]  S. P. Neuman,et al.  Effects of kriging and inverse modeling on conditional simulation of the Avra Valley Aquifer in southern Arizona , 1982 .

[46]  A. Journel,et al.  When do we need a trend model in kriging? , 1989 .

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

[48]  M. Hutchinson,et al.  Spatial interpolation of monthly mean climate data for China , 2005 .

[49]  A. Burnicki,et al.  Stochastic Simulation of Land-Cover Change Using Geostatistics and Generalized Additive Models , 2002 .

[50]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[51]  Jennifer A. Miller,et al.  Modeling the distribution of four vegetation alliances using generalized linear models and classification trees with spatial dependence , 2002 .

[52]  G. Heuvelink,et al.  A generic framework for spatial prediction of soil variables based on regression-kriging , 2004 .

[53]  C. Rolland.  Spatial and Seasonal Variations of Air Temperature Lapse Rates in Alpine Regions , 2003 .

[54]  P. Goovaerts.  Geostatistics in soil science: state-of-the-art and perspectives , 1999 .

[55]  C. Frei,et al.  Comparison of six methods for the interpolation of daily, European climate data , 2008 .

[56]  Gerhard Tutz,et al.  Variable Selection and Model Choice in Geoadditive Regression Models , 2009, Biometrics.

[57]  Budiman Minasny,et al.  Spatial prediction of topsoil salinity in the Chelif Valley, Algeria, using local ordinary kriging with local variograms versus whole-area variogram , 2001 .

[58]  D. Rizzo,et al.  Characterization of aquifer properties using artificial neural networks: Neural kriging , 1994 .

[59]  M. Araújo,et al.  The importance of biotic interactions for modelling species distributions under climate change , 2007 .

[60]  G. Robertson.  Geostatistics in Ecology: Interpolating With Known Variance , 1987 .

[61]  Peter Haggett,et al.  Trend-Surface Mapping in Geographical Research , 1965 .

[62]  David W. S. Wong,et al.  An adaptive inverse-distance weighting spatial interpolation technique , 2008, Comput. Geosci..