# Terrain Segmentation of Greece Using the Spatial and Seasonal Variation of Reference Crop Evapotranspiration

Authors: Vassilis Aschonitis; George Miliaresis; Kleoniki Demertzi; Dimitris Papamichail
Journal: Advances in Meteorology
 (2016)

---

## Abstract

The study presents a combination of techniques for integrated analysis of reference crop evapotranspiration (ETo) in GIS environment. The analysis is performed for Greece and includes the use of (a) ASCE-standardized Penman-Monteith method for the estimation of 50-year mean monthly ETo, (b) cross-correlation and principal components analysis for the analysis of the spatiotemporal variability of ETo, (c) K-means clustering for terrain segmentation to regions with similar temporal variability of ETo, and (d) general linear models for the description of ETo based on clusters attributes. Cross-correlation revealed a negative correlation of ETo with both elevation and latitude and a week positive correlation with longitude. The correlation between ETo and elevation was maximized during the warm season, while the correlation with latitude was maximized during winter. The first two principal components accounted for the 97.9% of total variance of mean monthly ETo. K-means segmented Greece to 11 regions/clusters. The categorical factor of cluster number together with the parameters of elevation, latitude, and longitude described satisfactorily the ETo through general linear models verifying the robustness of the cluster analysis. This research effort can contribute to hydroclimatic studies and to environmental decision support in relation to water resources management in agriculture.


---

## Body

## 1. Introduction

The evapotranspiration rate, which can be achieved under no water restrictions from a well-watered ideal grass surface, is called reference crop evapotranspirationETo and it is one of the most important hydroclimatic parameters for the implementation of various hydrological and agricultural applications [1]. Various methods have been developed for the assessment of ETo [2, 3] and the ASCE-standardized Penman-Monteith method, which is an update of the FAO-56 Penman-Monteith, has been proposed by the ASCE-EWRI Task Committee as the most precise method for ETo estimations [2].The intra-annual and spatial variability of theETo for regions of Greece has been investigated using different methods by many authors [4–18]. The latest and the most detailed mean monthly ETo estimations at country scale based on FAO Penman-Monteith method [1] have been performed by Dalezios et al. [5] using 66 meteorological stations covering a period of 15 years. The authors used the kriging technique for the development of 50 km resolution grids. Mardikis et al. [11] tested different interpolation methods for ETo based on 93 meteorological stations and provided improved methods which include the effect of elevation.The last years, the use of climatic models succeeded to generate high resolution grids of various climatic parameters which enabled the development of globalETo maps. Significant works of long-term mean ETo estimations at global scales have been performed by many authors [19–22]. ETo maps using the FAO-56 Penman-Monteith and ASCE-standardized Penman-Monteith method, which are the most demanding methods in terms of climatic data and elaboration, have been developed at resolutions of 10 arc-min [20] and 0.5 degrees [21] (1 degree = 60 arc-min = 3600 arc-sec ≈ 111 km at equator). The highest spatial resolution ETo maps (30 arc-sec) have been developed [22] using the simple method of Hargreaves [23], which requires only temperature data.Grids of mean monthlyETo can be used in an attempt to capture the spatial and seasonal patterns and to segment the terrain to regions with distinct seasonal ETo variability. Common approaches for such analysis are the typical multivariate methods such as correlation analysis, principal components, factorial kriging analysis, multi-Gaussian co-kriging, regression analysis, and so forth, which have successfully been used for other environmental quality parameters [24–30]. In addition, terrain segmentation techniques have been applied to multitemporal land surface temperature (LST) data in an attempt to define subregions with different seasonal LST variability, to assess its sensitivity to climatic change by revealing thermal anomalies and to support environmental analysis [31–36]. Similar regionalization and segmentation approaches have also been applied successfully in order to define subregions with different seasonal precipitation patterns [37–39].Objectives of the study are (a) to develop high resolution grids (30 arc-sec) of 50-years mean monthlyETo for Greece using the ASCE-standardized Penman-Monteith, (b) to analyze the spatiotemporal variability of ETo, and (c) to provide a terrain segmentation scheme for Greece based on the spatiotemporal variation of ETo. The techniques of cross-correlation, principal components (PCA), cluster analysis, and general linear models (GLMs) were used. Cross-correlation and PCA are used to capture the dependence of spatiotemporal variation of mean monthly ETo on geographical attributes (latitude and longitude), topography (elevation), and seasonality. K-means clustering is used to segment the Greek territory to regions/clusters in which the mean monthly ETo observations approximate to representative mean values with similar temporal variability. GLMs are used to (a) to parameterize the effect of clusters and together with elevation, latitude, and longitude to build models which can provide estimations of mean monthly ETo (b) to verify the robustness of the cluster analysis and (c) to assess the contribution/effect of each cluster in the mean monthly ETo estimations. This effort will be valuable in assisting environmental decision support and agricultural planning in relation to water resources management.


## 2. Data and Methods

### 2.1. Study Area and Data

The study area is Greece (South-East Europe) which is confined between the 34° and 42° parallel N., with a meridional extent from 19° to 28° E. Greece has a typical Mediterranean climate: relatively cold and rainy winters, relatively warm and dry summers, and, generally, extended periods of sunshine. The spatial heterogeneity of climate is mainly attributed to the mountain range of Pindos located in the central part of Greece, while significant differences in the winter’s severity are observed between the central-north continental territory and the islands where in the second case the winter is milder (Hellenic National Meteorological Service:http://www.hnms.gr/hnms/english/climatology/climatology_html).The analysis was based on climatic data, which were obtained from the following databases:(i)
[40]: this database provides mean monthly values for the parameters of maximum, minimum, and mean temperature at 30 arc-sec spatial resolution (~1 × 1 km). The data are provided as grids of mean monthly values of the period 1950–2000 (http://www.worldclim.org/). The database also includes a revised version of the GTOPO30 DEM based on SRTM DEM at 30 arc-sec spatial resolution. The data were produced using observation-based datasets after elaboration with the thin-plate smoothing spline algorithm implemented in the ANUSPLIN package for interpolation, using latitude, longitude, and elevation as independent variables.(ii)
[41]: this database provides mean monthly values of parameters such as solar radiation, specific humidity, wind speed at 10 m above surface (which was converted to 2 m above surface), precipitation, and temperature for the period 1948–2006 at 0.5 degrees spatial resolution. The dataset was constructed by combining a suite of global observation-based datasets and it is available in the form of NetCDF files (Network Common Data Form) of monthly values of each year for the period 1948–2006 (http://hydrology.princeton.edu/data.pgf.php). Resampling to 30 arc-sec spatial resolution was performed using the bilinear resampling scheme. The data were converted to grids of mean monthly values of the period 1950–2000.


### 2.2. Reference Crop Evapotranspiration

The estimation ofETo using the ASCE-standardized Penman-Monteith method is performed by the following equation [2]:(1)ETo=0.408ΔRn-G+γu2es-eaCn/Tmean+273.16Δ+γ1+Cdu2,where ETo is the daily reference crop evapotranspiration (mm d−1), Rn is the daily net radiation at the crop surface (MJ m−2 d−1), u2 is the mean daily wind speed at 2 m height above the soil surface (m s−1), Tmean is the mean daily air temperature (°C), G is the daily soil heat flux density at the soil surface (MJ m−2 d−1), es is the mean daily saturation vapour pressure (kPa), ea is the mean daily actual vapour pressure (kPa), Δ is the slope of the saturation vapour pressure-temperature curve (kPa °C−1), γ is the psychometric constant (kPa °C−1), Cn and Cd are constants, which vary according to the time step and the reference crop type and describe the bulk surface resistance and aerodynamic roughness. The short reference crop corresponds to clipped grass of 12 cm height and surface resistance of 70 s m−1 where the constants Cn and Cd have the values 900 and 0.34, respectively. The tall reference crop corresponds to full cover alfalfa of 50 cm height and surface resistance of 45 s m−1, where the constants Cn and Cd have the values 1600 and 0.38, respectively [2]. In this study, the ETo is estimated using (1) for the commonly used short reference crop as follows: (a) mean monthly values of maximum, minimum, and mean temperature were obtained from the database of [40], while (b) wind speed, specific humidity, and solar radiation were obtained from the database of [41]. The specific humidity was converted to actual vapour pressure [42] before its use. The equations used for intercalculations in ASCE method are given in [2, 43, 44]. All the calculations were performed in ArcGIS 9.3 ESRI environment. Twelve rasters (one for each month) of mean monthly ETo plus one of mean annual ETo of the period 1950–2000 were developed for the Greek territory.


### 2.3. Methodology for Terrain Segmentation Based on the Annual and Seasonal Variation of a Hydroclimatic Parameter

A set of four statistical methods such as cross-correlation analysis, principal component analysis,K-means cluster analysis, and general linear models [45–48] are used in this study. Correlation analysis and principal component analysis (PCA) reveal the temporal and spatial pattern evident within a multitemporal dataset. PCA is a linear transformation technique that produces a set of images known as principal components PCs that are uncorrelated with one another while they are ordered in terms of the amount of variance they explain from the original image set [49, 50]. PCs are computed from the linear combination of eigenvectors and the corresponding pixel values of the initial images [51]. PCA has traditionally been used in remote sensing as a means of data compaction since it is usually observed that the first 2 or 3 components are able to explain the majority of the variability in data values. Later components thus tend to be dominated by noise effects. By rejecting these later components, the volume of data is reduced, with no appreciable loss of information [52]. Standardized principal components analysis [53] is applied (data per month is centered with mean 0.0 and standard deviation 1.0) and so each image is not weighed according to its variance.K-means cluster analysis was used to partition the multitemporal (12-dimensional) imagery of ETo into K exclusive clusters. It begins by initializing cluster centroids, then assigns each pixel to the cluster whose centroid is nearest, updates the cluster centroids, and then repeats the process until the stopping criteria are satisfied [51]. The analysis uses Euclidian distance for calculating the distances between pixels and cluster centroids. The underlying idea of cluster analysis is that the cluster centroids represent the mean expression of the derived clusters. So clustering of the multitemporal data sets is expected to define groups of pixels with a rather common centroid curve that expresses their average monthly variability [32]. Elevation, latitude, and longitude statistics per cluster are computed in order to assist interpretation. Finally the clusters were interpreted according to their centroid and their spatial arrangement [36].General multiple linear regression analysis was performed using Statgraphics Centurion software (StatPoint Technologies). The parameterization of segmentation effects was made using an integer cluster number CN ranging from 1 toK. The CN categorical factor together with the quantitative factors of latitude (Lat), longitude (Lon), and elevation (H) were used as independent variables to describe the mean monthly values of ETo. The values of the aforementioned parameters were extracted from 29765 randomly selected positions homogeneously distributed in the entire Greek territory. The selection was performed using the Sampling Design Tool developed by NOAA (National Oceanic and Atmospheric Administration of USA) based on the procedure of “stratified random sampling.” For each month, the mean values of the dependent variable (mean monthly ETo) of sampling positions belonging to each cluster were compared using ANOVA-LSD. This analysis was performed in order to verify that the sampling procedure retained the differences between clusters, which were derived by K-means using the whole number of pixels of each cluster. Square root transformation was used for the Lat, Lon, and H values while BoxCox transformation [54] was used for the mean monthly ETo values to avoid/reduce normality deviations of the dependent and independent variables. The selection of the parameters for the regression model was performed using their variance inflation factor (VIF) in order to avoid multi-co-linearity effects. VIF values above 10 are usually considered to indicate serious multicollinearity [45]. Thus, the following linear model was built to describe the mean monthly ETo (mm) for each month:(2)BoxCoxETo,i=a0+∑j=1Naj·ICN,j+aN+1·Lon+aN+2·Lat+aN+3·H,where a0 is intercept of the regression model and aj is the regression coefficients of the categorical factor ICN which is regulated by the cluster number CN, aN+1,N+2,N+3: the regression coefficients of the quantitative factors Lon, Lat (decimal degrees), and H (m), respectively, and i is the month (starting with i=1 for January and ending with i=12 for December). The value of N is equal to the number of clusters minus 1 (N=K-1) whereas the categorical factor ICN is regulated by the cluster number CN as follows:(3)ICN,j=1,ifCN=j,-1,ifCN=N+1=K,0,ifCN≠j,K.The general form of BoxCox transformation for the dependent variableETo for each month i is given according to the following formula [54]:(4)BoxCoxETo,i=1+K1ETo,i+λ2λ1-1forλ1≠0,(5a)K1=1λ1K2λ1-1,(5b)K2=∏k=1nETo,i,k+λ21/n,where λ1,2 is the power and shift parameters of BoxCox transformation, respectively, and n is the number of ETo observations (samples) for a month i. Note that K2 is the geometric mean of ETo+λ2. In order to improve the efficiency of the model, the power parameter λ1 was optimized while the shift parameter λ2 was set to 0 based on the default optimization procedure incorporated in the statistical software. Τhe optimal transformation is the one that minimizes the mean squared error of the transformed dependent variable [54].Outliers were not removed during the regression analysis because of the large number of data (there is no automatic procedure to remove outliers). Autocorrelation of the residuals was tested using the Durbin-Watson test [55, 56], which provides values between 0 and 4 where the optimum value for no autocorrelation is equal to 2. ANOVA was used to estimate the statistical significance of the regression model while square correlation coefficient R2 adjusted for degrees of freedom was used to evaluate the explanatory power of the model. Type III sums of squares analysis was used to estimate the statistical significance of the independent variables CN, Lat, Lon, and H.


## 2.1. Study Area and Data

The study area is Greece (South-East Europe) which is confined between the 34° and 42° parallel N., with a meridional extent from 19° to 28° E. Greece has a typical Mediterranean climate: relatively cold and rainy winters, relatively warm and dry summers, and, generally, extended periods of sunshine. The spatial heterogeneity of climate is mainly attributed to the mountain range of Pindos located in the central part of Greece, while significant differences in the winter’s severity are observed between the central-north continental territory and the islands where in the second case the winter is milder (Hellenic National Meteorological Service:http://www.hnms.gr/hnms/english/climatology/climatology_html).The analysis was based on climatic data, which were obtained from the following databases:(i)
[40]: this database provides mean monthly values for the parameters of maximum, minimum, and mean temperature at 30 arc-sec spatial resolution (~1 × 1 km). The data are provided as grids of mean monthly values of the period 1950–2000 (http://www.worldclim.org/). The database also includes a revised version of the GTOPO30 DEM based on SRTM DEM at 30 arc-sec spatial resolution. The data were produced using observation-based datasets after elaboration with the thin-plate smoothing spline algorithm implemented in the ANUSPLIN package for interpolation, using latitude, longitude, and elevation as independent variables.(ii)
[41]: this database provides mean monthly values of parameters such as solar radiation, specific humidity, wind speed at 10 m above surface (which was converted to 2 m above surface), precipitation, and temperature for the period 1948–2006 at 0.5 degrees spatial resolution. The dataset was constructed by combining a suite of global observation-based datasets and it is available in the form of NetCDF files (Network Common Data Form) of monthly values of each year for the period 1948–2006 (http://hydrology.princeton.edu/data.pgf.php). Resampling to 30 arc-sec spatial resolution was performed using the bilinear resampling scheme. The data were converted to grids of mean monthly values of the period 1950–2000.


## 2.2. Reference Crop Evapotranspiration

The estimation ofETo using the ASCE-standardized Penman-Monteith method is performed by the following equation [2]:(1)ETo=0.408ΔRn-G+γu2es-eaCn/Tmean+273.16Δ+γ1+Cdu2,where ETo is the daily reference crop evapotranspiration (mm d−1), Rn is the daily net radiation at the crop surface (MJ m−2 d−1), u2 is the mean daily wind speed at 2 m height above the soil surface (m s−1), Tmean is the mean daily air temperature (°C), G is the daily soil heat flux density at the soil surface (MJ m−2 d−1), es is the mean daily saturation vapour pressure (kPa), ea is the mean daily actual vapour pressure (kPa), Δ is the slope of the saturation vapour pressure-temperature curve (kPa °C−1), γ is the psychometric constant (kPa °C−1), Cn and Cd are constants, which vary according to the time step and the reference crop type and describe the bulk surface resistance and aerodynamic roughness. The short reference crop corresponds to clipped grass of 12 cm height and surface resistance of 70 s m−1 where the constants Cn and Cd have the values 900 and 0.34, respectively. The tall reference crop corresponds to full cover alfalfa of 50 cm height and surface resistance of 45 s m−1, where the constants Cn and Cd have the values 1600 and 0.38, respectively [2]. In this study, the ETo is estimated using (1) for the commonly used short reference crop as follows: (a) mean monthly values of maximum, minimum, and mean temperature were obtained from the database of [40], while (b) wind speed, specific humidity, and solar radiation were obtained from the database of [41]. The specific humidity was converted to actual vapour pressure [42] before its use. The equations used for intercalculations in ASCE method are given in [2, 43, 44]. All the calculations were performed in ArcGIS 9.3 ESRI environment. Twelve rasters (one for each month) of mean monthly ETo plus one of mean annual ETo of the period 1950–2000 were developed for the Greek territory.


## 2.3. Methodology for Terrain Segmentation Based on the Annual and Seasonal Variation of a Hydroclimatic Parameter

A set of four statistical methods such as cross-correlation analysis, principal component analysis,K-means cluster analysis, and general linear models [45–48] are used in this study. Correlation analysis and principal component analysis (PCA) reveal the temporal and spatial pattern evident within a multitemporal dataset. PCA is a linear transformation technique that produces a set of images known as principal components PCs that are uncorrelated with one another while they are ordered in terms of the amount of variance they explain from the original image set [49, 50]. PCs are computed from the linear combination of eigenvectors and the corresponding pixel values of the initial images [51]. PCA has traditionally been used in remote sensing as a means of data compaction since it is usually observed that the first 2 or 3 components are able to explain the majority of the variability in data values. Later components thus tend to be dominated by noise effects. By rejecting these later components, the volume of data is reduced, with no appreciable loss of information [52]. Standardized principal components analysis [53] is applied (data per month is centered with mean 0.0 and standard deviation 1.0) and so each image is not weighed according to its variance.K-means cluster analysis was used to partition the multitemporal (12-dimensional) imagery of ETo into K exclusive clusters. It begins by initializing cluster centroids, then assigns each pixel to the cluster whose centroid is nearest, updates the cluster centroids, and then repeats the process until the stopping criteria are satisfied [51]. The analysis uses Euclidian distance for calculating the distances between pixels and cluster centroids. The underlying idea of cluster analysis is that the cluster centroids represent the mean expression of the derived clusters. So clustering of the multitemporal data sets is expected to define groups of pixels with a rather common centroid curve that expresses their average monthly variability [32]. Elevation, latitude, and longitude statistics per cluster are computed in order to assist interpretation. Finally the clusters were interpreted according to their centroid and their spatial arrangement [36].General multiple linear regression analysis was performed using Statgraphics Centurion software (StatPoint Technologies). The parameterization of segmentation effects was made using an integer cluster number CN ranging from 1 toK. The CN categorical factor together with the quantitative factors of latitude (Lat), longitude (Lon), and elevation (H) were used as independent variables to describe the mean monthly values of ETo. The values of the aforementioned parameters were extracted from 29765 randomly selected positions homogeneously distributed in the entire Greek territory. The selection was performed using the Sampling Design Tool developed by NOAA (National Oceanic and Atmospheric Administration of USA) based on the procedure of “stratified random sampling.” For each month, the mean values of the dependent variable (mean monthly ETo) of sampling positions belonging to each cluster were compared using ANOVA-LSD. This analysis was performed in order to verify that the sampling procedure retained the differences between clusters, which were derived by K-means using the whole number of pixels of each cluster. Square root transformation was used for the Lat, Lon, and H values while BoxCox transformation [54] was used for the mean monthly ETo values to avoid/reduce normality deviations of the dependent and independent variables. The selection of the parameters for the regression model was performed using their variance inflation factor (VIF) in order to avoid multi-co-linearity effects. VIF values above 10 are usually considered to indicate serious multicollinearity [45]. Thus, the following linear model was built to describe the mean monthly ETo (mm) for each month:(2)BoxCoxETo,i=a0+∑j=1Naj·ICN,j+aN+1·Lon+aN+2·Lat+aN+3·H,where a0 is intercept of the regression model and aj is the regression coefficients of the categorical factor ICN which is regulated by the cluster number CN, aN+1,N+2,N+3: the regression coefficients of the quantitative factors Lon, Lat (decimal degrees), and H (m), respectively, and i is the month (starting with i=1 for January and ending with i=12 for December). The value of N is equal to the number of clusters minus 1 (N=K-1) whereas the categorical factor ICN is regulated by the cluster number CN as follows:(3)ICN,j=1,ifCN=j,-1,ifCN=N+1=K,0,ifCN≠j,K.The general form of BoxCox transformation for the dependent variableETo for each month i is given according to the following formula [54]:(4)BoxCoxETo,i=1+K1ETo,i+λ2λ1-1forλ1≠0,(5a)K1=1λ1K2λ1-1,(5b)K2=∏k=1nETo,i,k+λ21/n,where λ1,2 is the power and shift parameters of BoxCox transformation, respectively, and n is the number of ETo observations (samples) for a month i. Note that K2 is the geometric mean of ETo+λ2. In order to improve the efficiency of the model, the power parameter λ1 was optimized while the shift parameter λ2 was set to 0 based on the default optimization procedure incorporated in the statistical software. Τhe optimal transformation is the one that minimizes the mean squared error of the transformed dependent variable [54].Outliers were not removed during the regression analysis because of the large number of data (there is no automatic procedure to remove outliers). Autocorrelation of the residuals was tested using the Durbin-Watson test [55, 56], which provides values between 0 and 4 where the optimum value for no autocorrelation is equal to 2. ANOVA was used to estimate the statistical significance of the regression model while square correlation coefficient R2 adjusted for degrees of freedom was used to evaluate the explanatory power of the model. Type III sums of squares analysis was used to estimate the statistical significance of the independent variables CN, Lat, Lon, and H.


## 3. Results and Discussion

### 3.1. Spatial Variation of Mean AnnualETo

The 30 arc-sec resolution (~1 × 1 km) maps of elevationH and mean reference crop evapotranspiration ETo for the period 1950–2000 estimated by the ASCE standardized Penman-Monteith method are given in Figures 1(a) and 1(b), respectively, while the frequency distributions of pixel values are given, respectively, in Figures 2(a) and 2(b).Figure 1
Maps of (a) elevationH and (b) mean annual ETo for the period 1950–2000 estimated by the ASCE standardized Penman-Monteith method.
(a)
(b)Figure 2
Frequency analysis of the spatial variation of (a) elevationH and (b) mean annual ETo estimated by the ASCE standardized Penman-Monteith method for the period 1950–2000.
(a)
(b)TheETo map (Figure 1(b)) indicates that the lowland areas (below 250 m) with a distance less than 20–25 km from the shoreline and with latitude less than 39 degrees present annual values of ETo>1200 mm. Regions with extremely high mean annual ETo values >1400 mm occupy approximately 13.5% of the Greek territory and they are distributed mainly to lowlands along the coastlines of (a) Crete and Dodecanese islands (southern and eastern Aegean Sea), (b) eastern Attika, and (c) eastern Peloponnesus. These findings are in accordance to the study of [5] who used a 50 km resolution grid and the studies of [57, 58] who analyzed data from meteorological stations.


### 3.2. Spatial and Seasonal Variation of Mean MonthlyETo

The correlation matrix of the mean monthly values of reference evapotranspirationETo versus H, Lat, and Lon is quantified in Table 1. According to Table 1, the mean monthly ETo presents negative correlations with H and Lat while a rather week positive correlation between ETo and Lon exists. The absolute values of ETo-H correlations are maximized during the warm season while the ETo-Lat correlations are maximized during the winter. The absolute values of ETo-Lon correlations are generally low and they are maximized during the period of August–January.Table 1
Cross-correlation matrix of the mean monthlyETo (January to December) as well as the latitude, longitude, and elevation images. Correlations are performed between pairs of variables taking values from pixels of the same geographical position.
H
Lat
Lon
E
T
o
Jan
Feb
Mar
Apr
May
Jun
Jul
Aug
Sep
Oct
Nov
Dec
H
1
Lat
0.03
1
Lon
−0.30
−0.14
1
Jan
−0.55
−0.70
0.32
1
Feb
−0.63
−0.65
0.21
0.97
1
Mar
−0.70
−0.60
0.15
0.93
0.98
1
Apr
−0.77
−0.56
0.22
0.88
0.94
0.96
1
May
−0.76
−0.54
0.24
0.83
0.88
0.92
0.98
1
Jun
−0.72
−0.58
0.23
0.82
0.87
0.90
0.95
0.98
1
Jul
−0.70
−0.53
0.25
0.77
0.81
0.84
0.91
0.97
0.98
1
Aug
−0.72
−0.56
0.36
0.82
0.84
0.85
0.92
0.96
0.98
0.99
1
Sep
−0.72
−0.61
0.43
0.88
0.89
0.90
0.94
0.96
0.96
0.95
0.98
1
Oct
−0.72
−0.61
0.29
0.93
0.96
0.97
0.98
0.95
0.95
0.91
0.93
0.96
1
Nov
−0.64
−0.67
0.34
0.98
0.98
0.96
0.94
0.91
0.89
0.85
0.88
0.94
0.98
1
Dec
−0.59
−0.67
0.35
0.99
0.98
0.94
0.90
0.86
0.85
0.79
0.84
0.90
0.95
0.99
1PCA was implemented in order to interpret the monthly intercorrelations observed in Table1 and the corresponding eigenvalues and eigenvectors are presented in Table 2. According to the eigenvalues, the first two principal components (PC-1 and PC-2) account for the 97.9% of the total variance observed in the 12 monthly ETo images.Table 2
Eigenvalues and eigenvectors for the PCs of the mean monthlyETo.
Principal components (PCs)
Eigenvectors
1
2
3
4
5
6
7
8
9
10
11
12
Jan
0.28
−0.41
0.33
−0.19
−0.40
−0.03
0.00
0.04
−0.44
0.33
−0.34
−0.17
Feb
0.29
−0.32
−0.16
−0.36
0.06
0.06
0.16
−0.52
0.10
0.21
0.50
0.20
Mar
0.29
−0.20
−0.49
−0.19
0.11
0.61
−0.02
0.22
0.03
−0.28
−0.28
−0.14
Apr
0.29
0.05
−0.46
0.39
−0.14
−0.39
0.09
−0.43
−0.09
−0.12
−0.39
0.11
May
0.29
0.25
−0.25
0.19
−0.58
0.04
−0.26
0.26
−0.06
0.06
0.51
−0.16
Jun
0.29
0.31
−0.08
−0.30
0.07
−0.31
0.66
0.39
0.11
0.13
−0.03
−0.09
Jul
0.28
0.44
0.11
−0.46
0.00
−0.05
−0.47
−0.02
0.01
−0.04
−0.21
0.48
Aug
0.28
0.35
0.34
−0.04
0.11
0.12
−0.01
−0.44
0.02
−0.25
0.01
−0.63
Sep
0.29
0.17
0.34
0.46
0.02
0.49
0.33
−0.01
−0.06
0.10
0.00
0.44
Oct
0.30
−0.04
−0.08
0.21
0.67
−0.16
−0.28
0.18
−0.39
0.28
0.16
−0.10
Nov
0.29
−0.24
0.15
0.20
0.02
−0.10
−0.23
0.11
0.78
0.27
−0.16
−0.10
Dec
0.28
−0.36
0.28
0.04
0.01
−0.29
0.01
0.21
−0.03
−0.71
0.21
0.16
Eigenvalue
11.13
0.62
0.15
0.04
0.02
0.01
0.01
0.01
0.01
0.00
0.00
0.00
Variance %
92.71
5.17
1.25
0.32
0.18
0.12
0.08
0.06
0.04
0.04
0.02
0.02The PC-1 map (Figure3(a)) is composed of linear combinations of monthly ETo images which are not influenced by seasonality since the weights of eigenvectors slightly vary between 0.28 and 0.30 (Table 2). Taking into account the correlation coefficients of monthly ETo versus H and Lat from Table 1, it is easy to detect that PC-1 amplifies the effects of elevation and latitude. The effects of elevation and latitude indirectly express the effects of temperature, which is higher in the lowlands but also higher in lower latitudes. The PC-1 is spatially maximised in the lowlands of the southern continental Greece and in the islands.Figure 3
The first two principal components PC-1 and PC-2 that account for the 97.9% of theETo variance.
(a)
(b)For PC-2 map (Figure3(b)), the linear combinations of monthly ETo images amplify the difference between the periods of October–March and April–September which present negative and positive weights of eigenvectors, respectively (Table 2). The difference between these two periods is spatially maximised in central Greece, in central Macedonia (Thessaloniki plain and Chalkidiki regions), in eastern Thrace, in the islands of north-eastern Aegean sea and south Ionian sea. These patterns are in accordance with the spatiotemporal patterns given by [5, 8, 59].


### 3.3. Terrain Segmentation Based on the Spatial and Seasonal Variation ofETo

For the implementation of theK-means method, the following criteria were used: (a) small clusters with area extent (occurrence) less than 0.001% were eliminated by merging them with larger clusters that are closest to their centroids and the stopping criterion was defined as the percentage of the migrating pixels during a specific iteration (if it was less than 0.001% of the entire image pixels the clustering was terminated), (b) the maximum number of clusters K was defined by a trial and error procedure, and (c) the maximum number of iterations was set equal to 300. At the end of the analysis, the mean monthly values of ETo are given for each cluster (described as multipart polygons).Cluster analysis revealed 11 clusters (K=11) in which the mean monthly ETo observations approximate to representative mean values with similar temporal variability. The characteristics of their centroids (latitude, longitude), their percent area extent, mean elevation, and the variation of monthly ETo are given in Table 3. The spatial distribution of the 11 clusters is given in Figure 4.Table 3
Cluster centroids, percent area (%), meanH, mean monthly, and annual ETo.
CN
Centroid
Area %
H (m)
E
T
o (mm)
Lat
Lon
Jan
Feb
Mar
Apr
May
Jun
Jul
Aug
Sep
Oct
Nov
Dec
Annual
1
39.98
22.91
32.12
488
18
28
52
81
118
153
172
156
104
64
31
21
997
2
40.00
22.51
13.58
987
12
20
41
67
103
136
154
139
90
52
22
14
850
3
40.11
22.24
3.68
1427
9
15
32
53
88
119
136
122
77
41
16
9
718
4
40.05
22.57
0.85
1791
9
14
28
45
75
103
119
107
66
33
13
8
620
5
38.90
22.40
20.38
341
26
36
61
93
130
168
189
171
116
75
40
30
1134
6
37.88
22.67
13.54
234
38
46
71
103
141
179
202
184
128
86
51
41
1269
7
37.86
23.53
5.99
131
42
49
74
109
149
188
215
198
138
93
57
45
1356
8
38.68
24.70
2.83
413
25
31
54
90
134
175
203
186
126
74
38
28
1164
9
37.25
26.52
2.55
133
47
50
74
107
153
199
232
218
155
99
62
51
1447
10
35.90
24.70
3.33
218
57
60
83
113
145
185
202
191
140
101
68
61
1404
11
35.56
25.29
1.15
94
67
68
91
120
150
192
208
200
148
110
78
72
1504
Total Greece
25
34
57
87
125
161
182
165
112
71
38
28
1085Figure 4
The spatial distribution of the 11 clusters ofETo (on the bottom left corner of cluster figures are also given the CN values).In general,ETo is maximized during summer and minimized during winter where the 76.7% of annual ETo in Greece is observed during April to September (warm-dry period), while the 23.3% is observed during October to March (cold-wet period) (Table 3). Clusters of significant interest are the following:(i)
Clusters 9, 10, and 11 present mean annual values ofETo>1400 mm, cluster 9 presents the highest mean monthly ETo among all clusters during the period of May–September, and cluster 11 presents the highest mean monthly ETo among all clusters during the period of October–April.(ii)
Clusters 3 and 4 present the lower mean annual values ofETo and correspond to the higher elevation areas of Greece.(iii)
Clusters 1 and 5 cover the major part of the agricultural land in Greece.


### 3.4. Regression Models of Mean MonthlyETo

The eleven clusters (K=11) obtained by the K-means analysis were used to calibrate (2) for each month using the observations of monthly ETo, Lon, Lat, H, and CN values of the 29765 selected sampling positions. The sampling positions were distributed to each cluster based on the % cluster’s area extent (Table 3). Their number within each cluster is given in Table 4. The comparison of the mean monthly values of the dependent variable (the BoxCox transformed mean monthly ETo) of sampling positions belonging to each cluster based on ANOVA-LSD is also given in Table 4. This analysis verified that the sampling procedure retained the differences between clusters which were derived by K-means taking into account the fact that only 4 out of 660 combinations of pairs comparisons were found nonstatistically different (e.g., pairs of clusters 3 and 4 for January, clusters 9 and 10 in June, clusters 7 and 9 in August, and clusters 9 and 10 in October) (Table 4).Table 4
Number of sampling positions within each cluster and ANOVA-LSD comparisons of the mean values of the dependent variable between the 11 clusters for each month.
CN
Sampling positions
Jan
Feb
Mar
Apr
May
Jun
Jul
Aug
Sep
Oct
Nov
Dec
1
9628
c†
d
d
d
d
d
d
d
d
d
d
d
2
4071
b
c
c
c
c
c
c
c
c
c
c
c
3
1103
a
b
b
b
b
b
b
b
b
b
b
b
4
256
a
a
a
a
a
a
a
a
a
a
a
a
5
6096
e
f
f
f
f
f
f
f
f
f
f
f
6
4008
f
g
g
g
g
g
g
g
g
g
g
g
7
1737
g
h
h
h
h
h
h
h
h
h
h
h
8
847
d
e
e
e
e
e
e
e
e
e
e
e
9
741
h
i
i
i
i
i
i
h
i
i
i
i
10
961
i
j
j
j
j
i
j
i
j
i
j
j
11
317
j
k
k
k
k
j
k
j
k
j
k
k
†Different letters indicate statistically significant differences between clusters of each month at 95% confidence level (ANOVA-LSD).The variance inflation factors VIFs of the selected independent variables of (2) ranged between 1.26 and 3.3 (the maximum VIF value corresponds to the transformed value of elevation). The regression coefficients and the statistical tests of general multiple regression analysis are given in Table 5 (twelve models, one for each month) whereas the respective graphs of observed versus predicted values of the BoxCox transformed mean monthly ETo are given in Figure 5. For all models, type III sums of squares analysis showed that the independent variables CN, Lat, Lon, and H were statistically significant at 95% confidence level with p<0.001, while the Durbin-Watson test indicated no serial autocorrelation of the residuals at 95% confidence level. Square correlation coefficients R2 ranged between 0.88 and 0.96 indicating satisfactory accuracy of the twelve models.Table 5
Regression coefficients and statistical tests of the general linear model described by  (2) for each month.
Month
Jan
Feb
Mar
Apr
May
Jun
Regression coefficient
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
a
0 (intercept)
185.01 ± 2.24
193.20 ± 2.10
256.76 ± 2.37
238.60 ± 2.25
248.98 ± 2.32
315.45 ± 2.80
a
1 (ICN,1)
−7.82 ± 0.07
−7.34 ± 0.07
−7.24 ± 0.08
−7.02 ± 0.07
−6.87 ± 0.07
−9.63 ± 0.09
a
2 (ICN,2)
−11.67 ± 0.11
−12.27 ± 0.10
−13.52 ± 0.11
−15.10 ± 0.11
−16.37 ± 0.11
−20.87 ± 0.13
a
3 (ICN,3)
−13.85 ± 0.16
−16.10 ± 0.15
−19.25 ± 0.17
−23.78 ± 0.16
−27.12 ± 0.17
−32.94 ± 0.20
a
4 (ICN,4)
−14.18 ± 0.27
−16.68 ± 0.25
−20.39 ± 0.28
−27.84 ± 0.27
−35.58 ± 0.28
−43.84 ± 0.34
a
5 (ICN,5)
−2.80 ± 0.07
−1.99 ± 0.07
−1.78 ± 0.07
0.93 ± 0.07
1.41 ± 0.07
1.43 ± 0.09
a
6 (ICN,6)
4.82 ± 0.08
5.02 ± 0.08
5.74 ± 0.09
9.02 ± 0.08
10.44 ± 0.09
9.99 ± 0.11
a
7 (ICN,7)
7.29 ± 0.11
7.43 ± 0.11
8.61 ± 0.12
13.04 ± 0.11
17.48 ± 0.12
19.03 ± 0.14
a
8 (ICN,8)
−3.48 ± 0.14
−4.47 ± 0.13
−4.84 ± 0.14
0.50 ± 0.14
8.35 ± 0.14
12.49 ± 0.17
a
9 (ICN,9)
9.12 ± 0.16
9.37 ± 0.15
11.41 ± 0.17
11.72 ± 0.16
21.73 ± 0.17
31.74 ± 0.20
a
10 (ICN,10)
14.24 ± 0.15
15.75 ± 0.14
17.56 ± 0.16
16.60 ± 0.15
12.40 ± 0.16
13.85 ± 0.19
a
11 sqrt(Lon)
0.87 ± 0.17
−6.90 ± 0.16
−13.08 ± 0.18
−6.29 ± 0.17
−5.49 ± 0.18
−9.14 ± 0.21
a
12 sqrt(Lat)
−23.14 ± 0.32
−17.79 ± 0.30
−21.13 ± 0.34
−20.94 ± 0.32
−18.99 ± 0.33
−24.79 ± 0.40
a
13 sqrt(H)
−0.24 ± 0.00
−0.28 ± 0.00
−0.43 ± 0.00
−0.51 ± 0.00
−0.47 ± 0.00
−0.42 ± 0.01
λ
1
0.525
0.785
1.156
1.375
1.365
1.555
K
1
8.278
2.671
0.463
0.137
0.127
0.039
R
2 (ANOVA p)
0.88 (p<0.001)
0.89 (p<0.001)
0.90 (p<0.001)
0.93 (p<0.001)
0.94 (p<0.001)
0.93 (p<0.001)
Durbin-Watson (p)
1.99 (p=0.24)
2.00 (p=0.50)
2.00 (p=0.43)
2.00 (p=0.47)
1.99 (p=0.28)
2.00 (p=0.43)
Month
Jul
Aug
Sep
Oct
Nov
Dec
Regression coefficient
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
a
0 (intercept)
431.49 ± 3.74
381.54 ± 3.25
258.85 ± 2.06
210.73 ± 2.29
192.59 ± 2.08
163.20 ± 2.31
a
1 (ICN,1)
−9.68 ± 0.12
−9.61 ± 0.10
−7.70 ± 0.07
−9.92 ± 0.07
−7.88 ± 0.07
−8.42 ± 0.07
a
2 (ICN,2)
−24.04 ± 0.18
−22.25 ± 0.15
−16.37 ± 0.10
−16.90 ± 0.11
−13.49 ± 0.10
−13.02 ± 0.11
a
3 (ICN,3)
−40.47 ± 0.27
−36.64 ± 0.23
−26.21 ± 0.15
−24.19 ± 0.16
−18.10 ± 0.15
−16.08 ± 0.17
a
4 (ICN,4)
−58.04 ± 0.45
−52.24 ± 0.39
−35.99 ± 0.25
−29.00 ± 0.27
−19.62 ± 0.25
−17.01 ± 0.28
a
5 (ICN,5)
3.35 ± 0.12
1.79 ± 0.10
0.03 ± 0.06
−1.65 ± 0.07
−1.59 ± 0.06
−2.23 ± 0.07
a
6 (ICN,6)
13.52 ± 0.14
11.33 ± 0.12
8.33 ± 0.08
6.64 ± 0.09
5.79 ± 0.08
5.03 ± 0.09
a
7 (ICN,7)
25.19 ± 0.19
22.14 ± 0.16
15.58 ± 0.10
12.26 ± 0.11
9.10 ± 0.10
7.56 ± 0.12
a
8 (ICN,8)
20.64 ± 0.23
17.23 ± 0.20
7.98 ± 0.13
−1.21 ± 0.14
−3.54 ± 0.13
−4.17 ± 0.14
a
9 (ICN,9)
42.02 ± 0.27
37.20 ± 0.23
25.93 ± 0.15
18.76 ± 0.17
11.75 ± 0.15
9.71 ± 0.17
a
10 (ICN,10)
12.08 ± 0.26
13.06 ± 0.22
12.23 ± 0.14
18.90 ± 0.16
16.19 ± 0.14
16.63 ± 0.16
a
11 sqrt(Lon)
−8.93 ± 0.28
5.37 ± 0.25
12.92 ± 0.16
−2.16 ± 0.17
3.21 ± 0.16
3.45 ± 0.17
a
12 sqrt(Lat)
−18.17 ± 0.53
−22.56 ± 0.46
−29.67 ± 0.29
−20.50 ± 0.33
−24.50 ± 0.30
−21.50 ± 0.33
a
13 sqrt(H)
−0.41 ± 0.01
−0.45 ± 0.01
−0.44 ± 0.00
−0.41 ± 0.00
−0.32 ± 0.00
−0.29 ± 0.00
λ
1
0.663
0.620
0.895
1.211
0.783
0.656
K
1
8.675
11.186
1.828
0.337
2.767
4.578
R
2 (ANOVA p)
0.91 (p<0.001)
0.93 (p<0.001)
0.96 (p<0.001)
0.93 (p<0.001)
0.92 (p<0.001)
0.89 (p<0.001)
Durbin-Watson (p)
2.00 (p=0.41)
2.00 (p=0.39)
2.00 (p=0.38)
2.01 (p=0.15)
1.99 (p=0.29)
1.99 (p=0.25)Figure 5
Observed versus predicted values of the BoxCox transformed mean monthlyETo using the models of Table 5.The analysis of the general linear models verified the robustness of the cluster analysis in two ways (a) by the statistically significant difference between the clusters, which was derived by the ANOVA-LSD (Table4), and (b) by the high statistical significance of the CN parameter as independent variable in the models and their high square correlation coefficients. The contribution/effect of each cluster in the spatiotemporal variation of ETo in the Greek territory, as it is described by the sum term inside (2), is given in Figure 6. According to Figure 6, clusters 1, 2, 3, and 4 present a negative effect while clusters 6, 7, 9, 10, and 11 present a positive effect to all months. The effect of clusters 5 and 8 presents small positive values during the warm season and small negative values during the cold season. The contribution of each cluster (Figure 6) defines the relative behaviour of each cluster during each month using, as a reference, the mean monthly conditions of ETo in the Greek territory.Figure 6
The contribution/effect of each cluster in the spatiotemporal variation ofETo in the Greek territory as it is described by the sum term inside (2).


## 3.1. Spatial Variation of Mean AnnualETo

The 30 arc-sec resolution (~1 × 1 km) maps of elevationH and mean reference crop evapotranspiration ETo for the period 1950–2000 estimated by the ASCE standardized Penman-Monteith method are given in Figures 1(a) and 1(b), respectively, while the frequency distributions of pixel values are given, respectively, in Figures 2(a) and 2(b).Figure 1
Maps of (a) elevationH and (b) mean annual ETo for the period 1950–2000 estimated by the ASCE standardized Penman-Monteith method.
(a)
(b)Figure 2
Frequency analysis of the spatial variation of (a) elevationH and (b) mean annual ETo estimated by the ASCE standardized Penman-Monteith method for the period 1950–2000.
(a)
(b)TheETo map (Figure 1(b)) indicates that the lowland areas (below 250 m) with a distance less than 20–25 km from the shoreline and with latitude less than 39 degrees present annual values of ETo>1200 mm. Regions with extremely high mean annual ETo values >1400 mm occupy approximately 13.5% of the Greek territory and they are distributed mainly to lowlands along the coastlines of (a) Crete and Dodecanese islands (southern and eastern Aegean Sea), (b) eastern Attika, and (c) eastern Peloponnesus. These findings are in accordance to the study of [5] who used a 50 km resolution grid and the studies of [57, 58] who analyzed data from meteorological stations.


## 3.2. Spatial and Seasonal Variation of Mean MonthlyETo

The correlation matrix of the mean monthly values of reference evapotranspirationETo versus H, Lat, and Lon is quantified in Table 1. According to Table 1, the mean monthly ETo presents negative correlations with H and Lat while a rather week positive correlation between ETo and Lon exists. The absolute values of ETo-H correlations are maximized during the warm season while the ETo-Lat correlations are maximized during the winter. The absolute values of ETo-Lon correlations are generally low and they are maximized during the period of August–January.Table 1
Cross-correlation matrix of the mean monthlyETo (January to December) as well as the latitude, longitude, and elevation images. Correlations are performed between pairs of variables taking values from pixels of the same geographical position.
H
Lat
Lon
E
T
o
Jan
Feb
Mar
Apr
May
Jun
Jul
Aug
Sep
Oct
Nov
Dec
H
1
Lat
0.03
1
Lon
−0.30
−0.14
1
Jan
−0.55
−0.70
0.32
1
Feb
−0.63
−0.65
0.21
0.97
1
Mar
−0.70
−0.60
0.15
0.93
0.98
1
Apr
−0.77
−0.56
0.22
0.88
0.94
0.96
1
May
−0.76
−0.54
0.24
0.83
0.88
0.92
0.98
1
Jun
−0.72
−0.58
0.23
0.82
0.87
0.90
0.95
0.98
1
Jul
−0.70
−0.53
0.25
0.77
0.81
0.84
0.91
0.97
0.98
1
Aug
−0.72
−0.56
0.36
0.82
0.84
0.85
0.92
0.96
0.98
0.99
1
Sep
−0.72
−0.61
0.43
0.88
0.89
0.90
0.94
0.96
0.96
0.95
0.98
1
Oct
−0.72
−0.61
0.29
0.93
0.96
0.97
0.98
0.95
0.95
0.91
0.93
0.96
1
Nov
−0.64
−0.67
0.34
0.98
0.98
0.96
0.94
0.91
0.89
0.85
0.88
0.94
0.98
1
Dec
−0.59
−0.67
0.35
0.99
0.98
0.94
0.90
0.86
0.85
0.79
0.84
0.90
0.95
0.99
1PCA was implemented in order to interpret the monthly intercorrelations observed in Table1 and the corresponding eigenvalues and eigenvectors are presented in Table 2. According to the eigenvalues, the first two principal components (PC-1 and PC-2) account for the 97.9% of the total variance observed in the 12 monthly ETo images.Table 2
Eigenvalues and eigenvectors for the PCs of the mean monthlyETo.
Principal components (PCs)
Eigenvectors
1
2
3
4
5
6
7
8
9
10
11
12
Jan
0.28
−0.41
0.33
−0.19
−0.40
−0.03
0.00
0.04
−0.44
0.33
−0.34
−0.17
Feb
0.29
−0.32
−0.16
−0.36
0.06
0.06
0.16
−0.52
0.10
0.21
0.50
0.20
Mar
0.29
−0.20
−0.49
−0.19
0.11
0.61
−0.02
0.22
0.03
−0.28
−0.28
−0.14
Apr
0.29
0.05
−0.46
0.39
−0.14
−0.39
0.09
−0.43
−0.09
−0.12
−0.39
0.11
May
0.29
0.25
−0.25
0.19
−0.58
0.04
−0.26
0.26
−0.06
0.06
0.51
−0.16
Jun
0.29
0.31
−0.08
−0.30
0.07
−0.31
0.66
0.39
0.11
0.13
−0.03
−0.09
Jul
0.28
0.44
0.11
−0.46
0.00
−0.05
−0.47
−0.02
0.01
−0.04
−0.21
0.48
Aug
0.28
0.35
0.34
−0.04
0.11
0.12
−0.01
−0.44
0.02
−0.25
0.01
−0.63
Sep
0.29
0.17
0.34
0.46
0.02
0.49
0.33
−0.01
−0.06
0.10
0.00
0.44
Oct
0.30
−0.04
−0.08
0.21
0.67
−0.16
−0.28
0.18
−0.39
0.28
0.16
−0.10
Nov
0.29
−0.24
0.15
0.20
0.02
−0.10
−0.23
0.11
0.78
0.27
−0.16
−0.10
Dec
0.28
−0.36
0.28
0.04
0.01
−0.29
0.01
0.21
−0.03
−0.71
0.21
0.16
Eigenvalue
11.13
0.62
0.15
0.04
0.02
0.01
0.01
0.01
0.01
0.00
0.00
0.00
Variance %
92.71
5.17
1.25
0.32
0.18
0.12
0.08
0.06
0.04
0.04
0.02
0.02The PC-1 map (Figure3(a)) is composed of linear combinations of monthly ETo images which are not influenced by seasonality since the weights of eigenvectors slightly vary between 0.28 and 0.30 (Table 2). Taking into account the correlation coefficients of monthly ETo versus H and Lat from Table 1, it is easy to detect that PC-1 amplifies the effects of elevation and latitude. The effects of elevation and latitude indirectly express the effects of temperature, which is higher in the lowlands but also higher in lower latitudes. The PC-1 is spatially maximised in the lowlands of the southern continental Greece and in the islands.Figure 3
The first two principal components PC-1 and PC-2 that account for the 97.9% of theETo variance.
(a)
(b)For PC-2 map (Figure3(b)), the linear combinations of monthly ETo images amplify the difference between the periods of October–March and April–September which present negative and positive weights of eigenvectors, respectively (Table 2). The difference between these two periods is spatially maximised in central Greece, in central Macedonia (Thessaloniki plain and Chalkidiki regions), in eastern Thrace, in the islands of north-eastern Aegean sea and south Ionian sea. These patterns are in accordance with the spatiotemporal patterns given by [5, 8, 59].


## 3.3. Terrain Segmentation Based on the Spatial and Seasonal Variation ofETo

For the implementation of theK-means method, the following criteria were used: (a) small clusters with area extent (occurrence) less than 0.001% were eliminated by merging them with larger clusters that are closest to their centroids and the stopping criterion was defined as the percentage of the migrating pixels during a specific iteration (if it was less than 0.001% of the entire image pixels the clustering was terminated), (b) the maximum number of clusters K was defined by a trial and error procedure, and (c) the maximum number of iterations was set equal to 300. At the end of the analysis, the mean monthly values of ETo are given for each cluster (described as multipart polygons).Cluster analysis revealed 11 clusters (K=11) in which the mean monthly ETo observations approximate to representative mean values with similar temporal variability. The characteristics of their centroids (latitude, longitude), their percent area extent, mean elevation, and the variation of monthly ETo are given in Table 3. The spatial distribution of the 11 clusters is given in Figure 4.Table 3
Cluster centroids, percent area (%), meanH, mean monthly, and annual ETo.
CN
Centroid
Area %
H (m)
E
T
o (mm)
Lat
Lon
Jan
Feb
Mar
Apr
May
Jun
Jul
Aug
Sep
Oct
Nov
Dec
Annual
1
39.98
22.91
32.12
488
18
28
52
81
118
153
172
156
104
64
31
21
997
2
40.00
22.51
13.58
987
12
20
41
67
103
136
154
139
90
52
22
14
850
3
40.11
22.24
3.68
1427
9
15
32
53
88
119
136
122
77
41
16
9
718
4
40.05
22.57
0.85
1791
9
14
28
45
75
103
119
107
66
33
13
8
620
5
38.90
22.40
20.38
341
26
36
61
93
130
168
189
171
116
75
40
30
1134
6
37.88
22.67
13.54
234
38
46
71
103
141
179
202
184
128
86
51
41
1269
7
37.86
23.53
5.99
131
42
49
74
109
149
188
215
198
138
93
57
45
1356
8
38.68
24.70
2.83
413
25
31
54
90
134
175
203
186
126
74
38
28
1164
9
37.25
26.52
2.55
133
47
50
74
107
153
199
232
218
155
99
62
51
1447
10
35.90
24.70
3.33
218
57
60
83
113
145
185
202
191
140
101
68
61
1404
11
35.56
25.29
1.15
94
67
68
91
120
150
192
208
200
148
110
78
72
1504
Total Greece
25
34
57
87
125
161
182
165
112
71
38
28
1085Figure 4
The spatial distribution of the 11 clusters ofETo (on the bottom left corner of cluster figures are also given the CN values).In general,ETo is maximized during summer and minimized during winter where the 76.7% of annual ETo in Greece is observed during April to September (warm-dry period), while the 23.3% is observed during October to March (cold-wet period) (Table 3). Clusters of significant interest are the following:(i)
Clusters 9, 10, and 11 present mean annual values ofETo>1400 mm, cluster 9 presents the highest mean monthly ETo among all clusters during the period of May–September, and cluster 11 presents the highest mean monthly ETo among all clusters during the period of October–April.(ii)
Clusters 3 and 4 present the lower mean annual values ofETo and correspond to the higher elevation areas of Greece.(iii)
Clusters 1 and 5 cover the major part of the agricultural land in Greece.


## 3.4. Regression Models of Mean MonthlyETo

The eleven clusters (K=11) obtained by the K-means analysis were used to calibrate (2) for each month using the observations of monthly ETo, Lon, Lat, H, and CN values of the 29765 selected sampling positions. The sampling positions were distributed to each cluster based on the % cluster’s area extent (Table 3). Their number within each cluster is given in Table 4. The comparison of the mean monthly values of the dependent variable (the BoxCox transformed mean monthly ETo) of sampling positions belonging to each cluster based on ANOVA-LSD is also given in Table 4. This analysis verified that the sampling procedure retained the differences between clusters which were derived by K-means taking into account the fact that only 4 out of 660 combinations of pairs comparisons were found nonstatistically different (e.g., pairs of clusters 3 and 4 for January, clusters 9 and 10 in June, clusters 7 and 9 in August, and clusters 9 and 10 in October) (Table 4).Table 4
Number of sampling positions within each cluster and ANOVA-LSD comparisons of the mean values of the dependent variable between the 11 clusters for each month.
CN
Sampling positions
Jan
Feb
Mar
Apr
May
Jun
Jul
Aug
Sep
Oct
Nov
Dec
1
9628
c†
d
d
d
d
d
d
d
d
d
d
d
2
4071
b
c
c
c
c
c
c
c
c
c
c
c
3
1103
a
b
b
b
b
b
b
b
b
b
b
b
4
256
a
a
a
a
a
a
a
a
a
a
a
a
5
6096
e
f
f
f
f
f
f
f
f
f
f
f
6
4008
f
g
g
g
g
g
g
g
g
g
g
g
7
1737
g
h
h
h
h
h
h
h
h
h
h
h
8
847
d
e
e
e
e
e
e
e
e
e
e
e
9
741
h
i
i
i
i
i
i
h
i
i
i
i
10
961
i
j
j
j
j
i
j
i
j
i
j
j
11
317
j
k
k
k
k
j
k
j
k
j
k
k
†Different letters indicate statistically significant differences between clusters of each month at 95% confidence level (ANOVA-LSD).The variance inflation factors VIFs of the selected independent variables of (2) ranged between 1.26 and 3.3 (the maximum VIF value corresponds to the transformed value of elevation). The regression coefficients and the statistical tests of general multiple regression analysis are given in Table 5 (twelve models, one for each month) whereas the respective graphs of observed versus predicted values of the BoxCox transformed mean monthly ETo are given in Figure 5. For all models, type III sums of squares analysis showed that the independent variables CN, Lat, Lon, and H were statistically significant at 95% confidence level with p<0.001, while the Durbin-Watson test indicated no serial autocorrelation of the residuals at 95% confidence level. Square correlation coefficients R2 ranged between 0.88 and 0.96 indicating satisfactory accuracy of the twelve models.Table 5
Regression coefficients and statistical tests of the general linear model described by  (2) for each month.
Month
Jan
Feb
Mar
Apr
May
Jun
Regression coefficient
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
a
0 (intercept)
185.01 ± 2.24
193.20 ± 2.10
256.76 ± 2.37
238.60 ± 2.25
248.98 ± 2.32
315.45 ± 2.80
a
1 (ICN,1)
−7.82 ± 0.07
−7.34 ± 0.07
−7.24 ± 0.08
−7.02 ± 0.07
−6.87 ± 0.07
−9.63 ± 0.09
a
2 (ICN,2)
−11.67 ± 0.11
−12.27 ± 0.10
−13.52 ± 0.11
−15.10 ± 0.11
−16.37 ± 0.11
−20.87 ± 0.13
a
3 (ICN,3)
−13.85 ± 0.16
−16.10 ± 0.15
−19.25 ± 0.17
−23.78 ± 0.16
−27.12 ± 0.17
−32.94 ± 0.20
a
4 (ICN,4)
−14.18 ± 0.27
−16.68 ± 0.25
−20.39 ± 0.28
−27.84 ± 0.27
−35.58 ± 0.28
−43.84 ± 0.34
a
5 (ICN,5)
−2.80 ± 0.07
−1.99 ± 0.07
−1.78 ± 0.07
0.93 ± 0.07
1.41 ± 0.07
1.43 ± 0.09
a
6 (ICN,6)
4.82 ± 0.08
5.02 ± 0.08
5.74 ± 0.09
9.02 ± 0.08
10.44 ± 0.09
9.99 ± 0.11
a
7 (ICN,7)
7.29 ± 0.11
7.43 ± 0.11
8.61 ± 0.12
13.04 ± 0.11
17.48 ± 0.12
19.03 ± 0.14
a
8 (ICN,8)
−3.48 ± 0.14
−4.47 ± 0.13
−4.84 ± 0.14
0.50 ± 0.14
8.35 ± 0.14
12.49 ± 0.17
a
9 (ICN,9)
9.12 ± 0.16
9.37 ± 0.15
11.41 ± 0.17
11.72 ± 0.16
21.73 ± 0.17
31.74 ± 0.20
a
10 (ICN,10)
14.24 ± 0.15
15.75 ± 0.14
17.56 ± 0.16
16.60 ± 0.15
12.40 ± 0.16
13.85 ± 0.19
a
11 sqrt(Lon)
0.87 ± 0.17
−6.90 ± 0.16
−13.08 ± 0.18
−6.29 ± 0.17
−5.49 ± 0.18
−9.14 ± 0.21
a
12 sqrt(Lat)
−23.14 ± 0.32
−17.79 ± 0.30
−21.13 ± 0.34
−20.94 ± 0.32
−18.99 ± 0.33
−24.79 ± 0.40
a
13 sqrt(H)
−0.24 ± 0.00
−0.28 ± 0.00
−0.43 ± 0.00
−0.51 ± 0.00
−0.47 ± 0.00
−0.42 ± 0.01
λ
1
0.525
0.785
1.156
1.375
1.365
1.555
K
1
8.278
2.671
0.463
0.137
0.127
0.039
R
2 (ANOVA p)
0.88 (p<0.001)
0.89 (p<0.001)
0.90 (p<0.001)
0.93 (p<0.001)
0.94 (p<0.001)
0.93 (p<0.001)
Durbin-Watson (p)
1.99 (p=0.24)
2.00 (p=0.50)
2.00 (p=0.43)
2.00 (p=0.47)
1.99 (p=0.28)
2.00 (p=0.43)
Month
Jul
Aug
Sep
Oct
Nov
Dec
Regression coefficient
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
Mean ± S.E.
a
0 (intercept)
431.49 ± 3.74
381.54 ± 3.25
258.85 ± 2.06
210.73 ± 2.29
192.59 ± 2.08
163.20 ± 2.31
a
1 (ICN,1)
−9.68 ± 0.12
−9.61 ± 0.10
−7.70 ± 0.07
−9.92 ± 0.07
−7.88 ± 0.07
−8.42 ± 0.07
a
2 (ICN,2)
−24.04 ± 0.18
−22.25 ± 0.15
−16.37 ± 0.10
−16.90 ± 0.11
−13.49 ± 0.10
−13.02 ± 0.11
a
3 (ICN,3)
−40.47 ± 0.27
−36.64 ± 0.23
−26.21 ± 0.15
−24.19 ± 0.16
−18.10 ± 0.15
−16.08 ± 0.17
a
4 (ICN,4)
−58.04 ± 0.45
−52.24 ± 0.39
−35.99 ± 0.25
−29.00 ± 0.27
−19.62 ± 0.25
−17.01 ± 0.28
a
5 (ICN,5)
3.35 ± 0.12
1.79 ± 0.10
0.03 ± 0.06
−1.65 ± 0.07
−1.59 ± 0.06
−2.23 ± 0.07
a
6 (ICN,6)
13.52 ± 0.14
11.33 ± 0.12
8.33 ± 0.08
6.64 ± 0.09
5.79 ± 0.08
5.03 ± 0.09
a
7 (ICN,7)
25.19 ± 0.19
22.14 ± 0.16
15.58 ± 0.10
12.26 ± 0.11
9.10 ± 0.10
7.56 ± 0.12
a
8 (ICN,8)
20.64 ± 0.23
17.23 ± 0.20
7.98 ± 0.13
−1.21 ± 0.14
−3.54 ± 0.13
−4.17 ± 0.14
a
9 (ICN,9)
42.02 ± 0.27
37.20 ± 0.23
25.93 ± 0.15
18.76 ± 0.17
11.75 ± 0.15
9.71 ± 0.17
a
10 (ICN,10)
12.08 ± 0.26
13.06 ± 0.22
12.23 ± 0.14
18.90 ± 0.16
16.19 ± 0.14
16.63 ± 0.16
a
11 sqrt(Lon)
−8.93 ± 0.28
5.37 ± 0.25
12.92 ± 0.16
−2.16 ± 0.17
3.21 ± 0.16
3.45 ± 0.17
a
12 sqrt(Lat)
−18.17 ± 0.53
−22.56 ± 0.46
−29.67 ± 0.29
−20.50 ± 0.33
−24.50 ± 0.30
−21.50 ± 0.33
a
13 sqrt(H)
−0.41 ± 0.01
−0.45 ± 0.01
−0.44 ± 0.00
−0.41 ± 0.00
−0.32 ± 0.00
−0.29 ± 0.00
λ
1
0.663
0.620
0.895
1.211
0.783
0.656
K
1
8.675
11.186
1.828
0.337
2.767
4.578
R
2 (ANOVA p)
0.91 (p<0.001)
0.93 (p<0.001)
0.96 (p<0.001)
0.93 (p<0.001)
0.92 (p<0.001)
0.89 (p<0.001)
Durbin-Watson (p)
2.00 (p=0.41)
2.00 (p=0.39)
2.00 (p=0.38)
2.01 (p=0.15)
1.99 (p=0.29)
1.99 (p=0.25)Figure 5
Observed versus predicted values of the BoxCox transformed mean monthlyETo using the models of Table 5.The analysis of the general linear models verified the robustness of the cluster analysis in two ways (a) by the statistically significant difference between the clusters, which was derived by the ANOVA-LSD (Table4), and (b) by the high statistical significance of the CN parameter as independent variable in the models and their high square correlation coefficients. The contribution/effect of each cluster in the spatiotemporal variation of ETo in the Greek territory, as it is described by the sum term inside (2), is given in Figure 6. According to Figure 6, clusters 1, 2, 3, and 4 present a negative effect while clusters 6, 7, 9, 10, and 11 present a positive effect to all months. The effect of clusters 5 and 8 presents small positive values during the warm season and small negative values during the cold season. The contribution of each cluster (Figure 6) defines the relative behaviour of each cluster during each month using, as a reference, the mean monthly conditions of ETo in the Greek territory.Figure 6
The contribution/effect of each cluster in the spatiotemporal variation ofETo in the Greek territory as it is described by the sum term inside (2).


## 4. Conclusions

The derived high resolution dataset of mean monthlyETo and the statistical analysis captured the seasonal and spatial variation of ETo in Greece while the regional spatiotemporal trends were revealed and assigned to spatial clusters. More specifically, cross-correlation and principal component analysis quantified the temporal pattern of ETo and revealed season, latitude, longitude, and elevation dependencies, while K-means cluster analysis segmented the terrain to regions with distinct temporal variability of ETo. The derived spatiotemporal patterns are in accordance with previous research efforts.The parameterization of each cluster (revealing a unique temporal signature of monthlyETo) together with topography (elevation) and spatial arrangement (latitude and longitude) described satisfactorily the mean monthly ETo through generalized linear models verifying the robustness of the cluster analysis, providing also a tool for ETo estimations for agricultural purposes. The feasible development of high resolution grids of monthly ETo based on available meteorological data produced by climatic models and the terrain segmentation approach can be used as tools to design decision support and management schemes of water recourses.The procedures described in this study can be used as a general framework for the spatial and seasonal variation not only forETo but also for other climatic parameters while it can support future research efforts such as(a)
assessing any climate changes by integrating monthlyETo data sets or other climatic variables of the fore-coming years; the new data sets can be used to recompute the temporal pattern of the climatic variable within each cluster in an attempt to evaluate any trends in the centroid change between two periods;(b)
using the proposed framework at watershed scale in order to facilitate the surface integration of climatic variables obtained by meteorological stations based on clusters produced by long-term data such as the ones of this study.



---

*Source: 3092671-2016-01-11.xml*