Geostatistical analysis of precipitation in the island of Crete (Greece) based on a sparse monitoring network

Based on the predictions of General Circulation Models, significant reduction of precipitation in Mediterranean areas is a possible scenario. Hence, better understanding of the spatial and temporal precipitation patterns is necessary in order to quantify desertification risks and design suitable mitigation measures. This study uses monthly precipitation measurements from a sparse network of 54 monitoring stations on the Mediterranean island of Crete (Greece). The study period extends from 1948 to 2012. The data reveal strong correlations between the western and eastern parts of the island. However, the average annual precipitation in the West is about 450 mm higher than that in the East. We construct a spatial model of average annual precipitation in Crete. The model involves a topographic trend and residuals with anisotropic spatial correlations which are fitted with a recently developed variogram function. We use regression kriging to generate annual precipitation maps and to identify locations of high estimation uncertainty. To our knowledge, this is the most detailed spatial analysis of precipitation in Crete to date. We present the analysis in detail for the year 1971. The trend accounts for ≈ 74% of the total variance. The highest precipitation estimate is 2331 mm in the West and 1781 mm in the East. The highest relative estimation uncertainty (≈ 20%) is observed along the southeastern coastline of the island, where the lowest values of annual precipitation are observed. This region includes one of the major agricultural areas of the island. The same overall patterns are found for other years in the study. Finally, we find no statistical evidence for a decrease in the global (over the entire island) annual precipitation during the study period.

[1]  Konstantine P. Georgakakos,et al.  A hydrologically useful station precipitation model: 1. Formulation , 1984 .

[2]  L. D. Monache,et al.  Bayesian Exploration of Multivariate Orographic Precipitation Sensitivity for Moist Stable and Neutral Flows , 2015 .

[3]  Patrick L Brezonik,et al.  Analysis and predictive models of stormwater runoff volumes, loads, and pollutant concentrations from watersheds in the Twin Cities metropolitan area, Minnesota, USA. , 2002, Water research.

[4]  David I. F. Grimes,et al.  Geostatistical analysis of rainfall , 2010 .

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

[6]  Phaedon C. Kyriakidis,et al.  Geostatistical Integration of Coarse Resolution Satellite Precipitation Products and Rain Gauge Data to Map Precipitation at Fine Spatial Resolutions , 2017, Remote. Sens..

[7]  J. Friesen,et al.  Adjustment of global precipitation data for enhanced hydrologic modeling of tropical Andean watersheds , 2017, Climatic Change.

[8]  Andreas N. Angelakis,et al.  Water resources management in the Island of Crete, Greece, with emphasis on the agricultural use , 2001 .

[9]  Pierre Ailliot,et al.  Stochastic weather generators: an overview of weather type models , 2015 .

[10]  U. Haberlandt,et al.  Spatial interpolation of hourly rainfall - effect of additional information, variogram inference and storm properties , 2010 .

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

[12]  Dionissios T. Hristopulos,et al.  Spartan Gibbs Random Field Models for Geostatistical Applications , 2003, SIAM J. Sci. Comput..

[13]  Robert M. Hirsch,et al.  Has the magnitude of floods across the USA changed with global CO2 levels? , 2012 .

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

[15]  M. Borga,et al.  Precipitation and temperature space–time variability and extremes in the Mediterranean region: evaluation of dynamical and statistical downscaling methods , 2013, Climate Dynamics.

[16]  J. Chilès,et al.  Geostatistics: Modeling Spatial Uncertainty , 1999 .

[17]  A. Dégre,et al.  Geostatistical interpolation of daily rainfall at catchment scale: the use of several variogram models in the Ourthe and Ambleve catchments, Belgium , 2011 .

[18]  Halit Apaydin,et al.  Evaluation of topographical and geographical effects on some climatic parameters in the Central Anatolia Region of Turkey , 2011 .

[19]  D. Allard,et al.  A stochastic daily weather generator for skewed data , 2010 .

[20]  P. Burman.  A comparative study of ordinary cross-validation, v-fold cross-validation and the repeated learning-testing methods , 1989 .

[21]  P. Nastos,et al.  Spatial and temporal variability of Aridity Index in Greece , 2010 .

[22]  E. Xoplaki,et al.  A fast nonparametric spatio‐temporal regression scheme for generalized Pareto distributed heavy precipitation , 2014 .

[23]  B. Minasny,et al.  The Matérn function as a general model for soil variograms , 2005 .

[24]  G. Tsakiris,et al.  Regional Drought Assessment Based on the Reconnaissance Drought Index (RDI) , 2007 .

[25]  D. Easterling,et al.  Observations: Atmosphere and surface , 2013 .

[26]  M. Stone.  Cross‐Validatory Choice and Assessment of Statistical Predictions , 1976 .

[27]  T. W. Anderson,et al.  An Introduction to Multivariate Statistical Analysis , 1959 .

[28]  Statistical Tests of Taylor's Hypothesis: An Application to Precipitation Fields , 2009 .

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

[30]  D. Hristopulos.  Covariance functions motivated by spatial random field models with local interactions , 2014, Stochastic Environmental Research and Risk Assessment.

[31]  A. Longobardi,et al.  Spatial and temporal distribution of precipitation in a Mediterranean area (southern Italy) , 2016, Environmental Earth Sciences.

[32]  Peter M. Atkinson,et al.  Mapping precipitation in Switzerland with ordinary and indicator kriging. Special issue: Spatial Interpolation Comparison 97 , 1998 .

[33]  Cristina Portalés,et al.  Seasonal precipitation interpolation at the Valencia region with multivariate methods using geographic and topographic information , 2009 .

[34]  T. McKee,et al.  THE RELATIONSHIP OF DROUGHT FREQUENCY AND DURATION TO TIME SCALES , 1993 .

[35]  G. Corzo,et al.  Spatio-temporal analysis of annual rainfall in Crete, Greece , 2018, Acta Geophysica.

[36]  Ker-Chau Li,et al.  Asymptotic Optimality for $C_p, C_L$, Cross-Validation and Generalized Cross-Validation: Discrete Index Set , 1987 .

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

[38]  David Higdon,et al.  Non-Stationary Spatial Modeling , 2022, 2212.08043.

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

[40]  P. Jones,et al.  REPRESENTING TWENTIETH CENTURY SPACE-TIME CLIMATE VARIABILITY. , 1998 .

[41]  P. Guttorp,et al.  Studies in the history of probability and statistics XLIX On the Matérn correlation family , 2006 .

[42]  Dionissios T. Hristopulos,et al.  Improvement of groundwater level prediction in sparsely gauged basins using physical laws and local geographic features as auxiliary variables , 2013 .

[43]  Ricardo A. Olea,et al.  Geostatistics for Engineers and Earth Scientists , 1999, Technometrics.

[44]  Shaun Lovejoy,et al.  Universal Multifractals: Theory and Observations for Rain and Clouds , 1993 .

[45]  O. Tzoraki,et al.  Spatial Water Use efficiency Index towards resource sustainability: application in the island of Crete, Greece , 2015 .

[46]  Ioannis K. Tsanis,et al.  Temporal and spatial variation of annual rainfall on the island of Crete, Greece , 2003 .

[47]  C. W. Richardson,et al.  Comparison of the WGEN and LARS-WG stochastic weather generators for diverse climates , 1998 .

[48]  D. Rossiter,et al.  Multi-scale characterization of soil variability within an agricultural landscape mosaic system in southern Cameroon , 2005 .

[49]  L. V. Watrous.  Lasithi, a history of settlement on a highland plain in Crete , 1983 .

[50]  Dionissios T. Hristopulos,et al.  Improving kriging of groundwater level data using nonlinear normalizing transformations—a field application , 2012 .

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

[52]  Anja Vogler,et al.  An Introduction to Multivariate Statistical Analysis , 2004 .

[53]  B. Katsoulis,et al.  Spatial and Temporal Variation of Precipitation in Greece and Surrounding Regions Based on Global Precipitation Climatology Project Data , 2008 .

[54]  Luca Vogt.  Statistics For Spatial Data , 2016 .

[55]  Sylvain Arlot,et al.  A survey of cross-validation procedures for model selection , 2009, 0907.4728.

[56]  Uwe Haberlandt,et al.  Geostatistical interpolation of hourly precipitation from rain gauges and radar for a large-scale extreme rainfall event , 2007 .

[57]  B. Efron.  Estimating the Error Rate of a Prediction Rule: Improvement on Cross-Validation , 1983 .

[58]  I. Tsanis,et al.  Drought Assessment Based on Multi-Model Precipitation Projections for the Island of Crete , 2013 .

[59]  Jacques Rivoirard.  On the Structural Link Between Variables in Kriging with External Drift , 2002 .

[60]  Dionissios T. Hristopulos,et al.  Analytic Properties and Covariance Functions for a New Class of Generalized Gibbs Random Fields , 2006, IEEE Transactions on Information Theory.

[61]  Phaedon C. Kyriakidis,et al.  Geostatistical Mapping of Precipitation from Rain Gauge Data Using Atmospheric and Terrain Characteristics , 2001 .

[62]  Dionissios T. Hristopulos,et al.  Erratum: Spartan Gibbs Random Field Models for Geostatistical Applications , 2005, SIAM J. Sci. Comput..

[63]  A. Yaglom.  Correlation Theory of Stationary and Related Random Functions I: Basic Results , 1987 .

[64]  Dionissios T. Hristopulos,et al.  An application of Spartan spatial random fields in environmental mapping: focus on automatic mapping capabilities , 2008 .

[65]  Richard W. Katz,et al.  Precipitation as a Chain-Dependent Process , 1977 .

[66]  Jennifer L. Dungan,et al.  Comparison of regression and geostatistical methods for mapping Leaf Area Index (LAI) with Landsat ETM+ data over a boreal forest. , 2005 .

[67]  A. Bárdossy,et al.  Interpolation of precipitation under topographic influence at different time scales , 2013 .

[68]  Jin Li,et al.  Spatial interpolation methods applied in the environmental sciences: A review , 2014, Environ. Model. Softw..

[69]  F. Moral.  Comparison of different geostatistical approaches to map climate variables: application to precipitation , 2010 .

[70]  Hans Kaper,et al.  Mathematics and Climate , 2013 .

[71]  Gerard B. M. Heuvelink,et al.  Comparison of regression and kriging techniques for mapping the average annual precipitation of Turkey , 2012, Int. J. Appl. Earth Obs. Geoinformation.

[72]  P. Jones,et al.  Representing Twentieth-Century Space-Time Climate Variability. Part II: Development of 1901-96 Monthly Grids of Terrestrial Surface Climate , 2000 .

[73]  M. Zukovic,et al.  SPARTAN RANDOM PROCESSES IN TIME SERIES MODELING , 2007, 0709.3418.

[74]  D. T. Hristopulos,et al.  Comparison of stochastic and deterministic methods for mapping groundwater level spatial variability in sparsely monitored basins , 2012, Environmental Monitoring and Assessment.

[75]  Dionissios T. Hristopulos,et al.  Comparison of spatiotemporal variogram functions based on a sparse dataset of groundwater level variations , 2017 .

[76]  M. Parlange,et al.  Improved interpolation of meteorological forcings for hydrologic applications in a Swiss Alpine region , 2011 .

[77]  R. Adler.  The Geometry of Random Fields , 2009 .

[78]  A. Davison,et al.  Geostatistics of Dependent and Asymptotically Independent Extremes , 2013, Mathematical Geosciences.

[79]  M. Widmann,et al.  Statistical downscaling of GCM-simulated precipitation using Model Output Statistics , 2014 .

[80]  Giovanni Seni,et al.  Ensemble Methods in Data Mining: Improving Accuracy Through Combining Predictions , 2010, Ensemble Methods in Data Mining.

[81]  Anastassia Baxevani,et al.  A spatiotemporal precipitation generator based on a censored latent Gaussian field , 2015 .

[82]  F M Danson,et al.  Modelling the spatial distribution of Echinococcus multilocularis infection in foxes. , 2004, Acta tropica.