Effect of Climate Change on Spatio-Temporal Variability and Trends of Evapotranspiration , and Its Impact on Water Resources Management in The Kingdom of Saudi Arabia

© 2013 Elnesr and Alazba, licensee InTech. This is an open access chapter distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Effect of Climate Change on Spatio-Temporal Variability and Trends of Evapotranspiration, and Its Impact on Water Resources Management in The Kingdom of Saudi Arabia


Introduction
Recently, climate change is receiving much attention.Changes in the world's climate have significant effect on water resources which affect the livelihood of people especially in hyper arid regions such as the Kingdom of Saudi Arabia (KSA).The KSA suffers an enduring water shortage problem, despite the fact that the agricultural activities consume up to 90% of the water amount in the Kingdom.Reference Evapotranspiration (ETo) is an agro-climatic property that involves temperature, humidity, solar radiation, and wind speed.Identifying changes in ETo can also help in future planning of agriculture-water projects and identify lower and higher ETo zones for proper planning and management of agricultural projects in arid regions.

Water resources and climate change
Water shortage is a swelling problem in the arid and semi-arid regions.Affected by its geographic location and its climate, the Kingdom of Saudi Arabia (KSA) suffers a severe water deficit.Even rain, which is the only renewable water source, comes in flash short duration storms of high intensity and most of it vanishes to evaporation.Thus, almost all agriculture of the kingdom is irrigated.Irrigation water, though, consumes 80 to 88 % of the total water consumption (Abu-Ghobar, 2000;Abderrahman, 2001).In addition to these water scarcity conditions, but it seems getting scary by the effects of climate change on the hydrological cycle and water supply.The quantity of irrigation water is determined initially by identifying the reference evapotranspiration (ETo).Several researches was conducted to detect climate changes, trends and variability in various parts of the world using some climate parameters such as air temperature, rainfall depth, ETo, and pan evapotranspiration ETp (Shwartz and Randall, 2003;Garbrecht, et al., 2004;Hegerl, et al., 2007;Fu, et al., 2009;Hakan, et al., 2010;Elnesr and Alazba, 2010;Elnesr et al., 2010a;and Elnesr et al. 2010b).The ETo parameter has a special importance because it combines changes in many other climate parameters including temperature, radiation, humidity, and wind speed.It has, however, direct influence on hydrologic water balance, irrigation and drainage canal design, reservoir operation, potentials for rain-fed agricultural production, and crop water requirements (Dinpashoh, 2006).

Climate change effect on evapotranspiration worldwide
Several studies conducted in North America have shown that some climate parameters are on the rise including ETo (Fehrman, 2007;Garbrecht et al., 2004;Szilagyi, 2001).Fehrman, 2007 found an increasing trend in ETo over the Mississippi area and that most of ETo increase can be attributed to the increase in July.He also found that the rate of ETo increase was 0.29 mm/years when his study period extended from 1940 to 1999 compared to 0.88 mm/year when the study period was limited to 1950 to 1999 records.The accelerated ET over North America is presumed to be due to a rise in temperature over the past century (Myeni et al. 1997, Milly andDunne 2001).In the contrary ETo and pan evaporation has shown to decrease in China (Thomas, 2000, Liu et al., 2004) and at a rate of 1.19 mm/year (Song et al., 2010) despite the rise in maximum daily temperature.In the Tibetan Plateau ETo decreases as well at a rate of 1.31 mm/year or 2.0% of the annual total evapotranspiration (Shenbin et al., 2006).The decrease in ETo has been attributed to the decrease in wind speed and net radiation.In another study Gao et al., (2007) found that the actual evapotranspiration had a decreasing trend in most of the eastern part of china and there was an increasing trend in the western and the northern parts of northeast China and that the change in precipitation played a key role for the change of estimated actual evapotranspiration.Similar negative trends in pan evaporation were found in 24 out of 27 observation stations in a 19-year study in Thailand (Tebakari et al., 2005).In India, a significant decreasing trend was found in ETo all over the Indian plateau during the past 40 years, which was mainly caused by a significant increase in the relative humidity and a consistent significant decrease in the wind speed throughout the country (Bandyopadhyay et al., 2009).In Australia, Roderick and Farquhar (2004) found a decreasing trend in pan evaporation and conclude that Australia is becoming less arid.However, there is enough evidence now that a decrease in pan evaporation is an indicator to an increase in actual evaporation.This is what known now as the evaporation paradox (Hobbins et al., 2004).Some researchers developed a hypothetical scenario to study the effect of possible increase on temperature over the ETo and subsequently on water supply.A study conducted by Abderrahman et al. (1991) concluded that in the KSA, a 1℃ increase in temperature would increase ETo from 1 to 4.5%.In another study, that includes selected cities in KSA, United Arab Emirates and Kuwait, Abderahman and Al-Harazin (2003) concluded that an increase in temperature by 1℃ would increase ETo over these area by a maximum of 20%.In general, studies involving ETo calculation seemed to be more limited worldwide compared to other climate parameters.In the other hand, regarding other climatic parameters, Hakan et al. (2010) reported an increasing trend in temperature and ETo in most of stations they analyzed in Turkey using Mann-Kendall analysis.Cohen and Stanhill, (1996) studied rainfall changes in the Jordan Valley/Jordan and found a tangible but insignificant decrease at a rate of -0.47 and -0.16 mm/year for two different stations.Similar conclusions were observed by Al-Ansari et al (1999) who observed a general decrease in rainfall intensity.Smadi (2006), and Smadi and Zghoul (2006) found a prompt shift in rainfall and temperature in Jordan.ElNesr et al (2010b) concluded that the Saudi Arabia and the Arabian Peninsula are suffering from a considerable warming trend form year 1980 to 2008.Still, Elnesr et al. (2010a) concluded that the percentage land area with annual ETo>4000 mm increased from about 20% to 40% in the period they studied.On the other hand, lower ETo values, less than 3600 mm, contracted from about 30% to 12%.

Objective of the study
This study aims to trace the ETo values over time throughout all the area of the Saudi Arabia, then to quantify the future of water demand according to the ETo trends

Geography of the Saudi Arabia
Saudi Arabia is the largest country of the Arabian Peninsula; it occupies about 80% of its area (Wynbrandt, 2004).The country lies between latitudes 16°21'58''N, and 32°9'57''N, and longitudes 34°33'48''E and 55°41'29''E, as illustrated in Fig. 1.Saudi Arabia has a desert dry climate with high temperatures in most of the country.However, the country falls in the tropical and subtropical desert region.Winds reaching the country are generally dry, and almost all the area is arid.Because of the aridity and the relatively cloudless skies, there are great extremes in temperature, but there are also wide variations between the seasons and regions (AQUASTAT, 2008).

Evapotranspiration calculation
Evapotranspiration was calculated using Food and Agricultural Organization (FAO) Penman-Monteith (PM) procedure, FAO 56 method, presented by Allen et al. (1998).In this method, ETo is expressed as follows: where ETo is the daily reference evapotranspiration  The measured meteorological data available were Ta, Relative humidity (RH) and U2 whereas soil heat flux (G) was taken equal to zero, (Allen et al, 2005).The slope of the saturation vapor pressure curve (Δ) is computed by the following equation as in Murray (1967): where o a e T   is calculated according to (Tetens, 1930): The net radiation Rn was estimated as the difference between the net short wave incoming radiation Rns and the net long wave outgoing radiation Rnl.The calculation of Rns, and Rnl, followed the procedures outlined in Allen et al. (1998) and Doorenbos and Pruitt (1977).All radiation was computed in daily energy flux units (MJ m −2 day −1 ).where z: altitude [m].The latent heat λ depends on the average temperature, Eqn(7) , while it can be taken as an approximate value of 2.45 as reported by Harrison (1963) for Ta=20 ℃.
In the current study, we chose to calculate the latent heat using Eqn(7) .
2.5 0.00236 a T   The saturation vapour pressure, es, and actual vapour pressure, ea, are calculated according to Allen et al (2005) as: where RHx, RHn: maximum and minimum relative humidity [%] respectively.
The average daily ETo in a specific month was calculated by taking the arithmetic average of the daily values in that month.The summation of all ETo daily values in a year for a station will give the total annual ETo for that station.

Climatic data source and description
Basic climatic data were taken from the Presidency of Meteorology and Environment in KSA, the official climate agency in the country.The data set is the most accurate one in KSA and used by all other governmental and academic agencies for climate research and prediction.Weather stations are equipped with up-to-date monitoring devices and subjected to regular inspection and replacement for defected devices (personal communication with the Presidency of Meteorology and Environment).Data represents 29 meteorological stations as shown in Fig. 1.These stations represent all the 13 districts of the KSA.The data covers 29 years of daily meteorological records for 20 stations, 24 years for 6 stations, and 3 stations with less than 20 years as shown in

Data grouping and contouring
After correction the data sets, daily ETo values were calculated for each station, then aggregated to annual and monthly values.Annual ETo value (mm/year) for each station was calculated by summation of the daily ETo for the entire year.On the other hand, the monthly average ETo value was calculated by taking the average of the daily ETo values during each month.
Evapotranspiration data were graphically represented by contour maps irrespective of stations altitude.Analysis of ET variations with stations' altitude for each of the 30 years under study revealed no trends.Other researchers have also found no correlation between ET and station altitude in China (Thomas, 2000).Contour maps present clearly zones of common ET values as well as clarify vividly ET differences between zones and viability a long months or years.This approach has also been adopted by other researchers to study ET variability in China (Thomas, 2000;Shenbin et al., 2006).
Data was arranged in three columns format namely, longitude, latitude, and ETo.Each set of data was gridded separately using the ordinary point-Kriging method which estimates the values of the points at the grid nodes (Abramowitz andStegun, 1972, andIsaaks andSrivastava, 1989).This procedure is used by SURFER™ Software which has been used in our calculations.The resulted grid was blanked outside the political borders of the KSA.The political borders' information of the KSA was grabbed from electronic map of NIMA ( 2003).
The electronic map was digitized and converted to DMS geographic coordinate system.The blanked grid was plotted as a contour map using Surfer™ 8.0 software (Surfer, 2002).Sample plots for the average daily ETo during a month, June in this case and the ETo, in a year, 1991 in this case, is shown in  All of the data are daily values, the obtained climatic data records were carefully inspected for missing and erroneous reading.Very few errors were found, (median value of 00.45%).
Errors were classified into four categories: Errors because of mistaken extreme values such as a relative humidity exceeds 100% or below 0%.Illogical errors such as the recorded maximum daily temperature (Tx) was less than the minimum daily temperature (Tn) in the same day, or if Tx = Tn.Missing values; i.e.Tx is present but Tn is missing.Recording an error-indication number (like 999 or 777) if the sensor is not functioning.On analyzing the data record, any value contains one or more errors was considered missing record unless the missing record could be predicted with minimal error, i.e. if the average temperature (Ta) is missing while Tx and Tn are logged with no errors; in this case Ta = (Tx + Tn)/2.However, the amount of missing data in the recorded period could be considered negligible in most of the stations, where the average amount of missing data is 0.78%.

Mann-kendall test
The Mann-Kendall test is a non-parametric test used for identifying trends in time series data.The test compares the relative magnitudes of sample data rather than the data values themselves Both Kendall tau coefficient (τ) and Mann-Kendall coefficient (s) are nonparametric statistics used to find rank correlation.Kendall (τ) is a ratio between the actual rating score of correlation, to the maximum possible score.To obtain the rating score for a time series, the dataset is sorted in ascending order according to time, and then the following formula is applied: where s: the rating score (also called the Mann-Kendall sum); x: the data value; i and j: counters; n: number of data values in the series; Sign is a function having values of +1, 0, or -1 if (xi-xj) is positive, zero, or negative, respectively.According this formula, the maximum value of s is: Hence, the Kendall (τ) is calculated as: A positive value of s or  is an indicator of an increasing trend, and a negative value indicates a decreasing trend.However, it is necessary to compute the probability associated with s or  and the sample size, n, to quantify the significance of the trend statistically.Kendall and Gibbons (1990) introduced a normal-approximation test that could be applied on datasets of more than ten values with s variance (σ 2 ): where CFR: repetition correction factor, to fix the effect of tied groups of data (when some of the data values appear more than one time in the dataset, this group of values are called a tied group); g: number of tied groups; k: a counter; m: number of data values in each tied group.Then normal distribution parameter (called the Mann-Kendall statistic, Z) is calculated as follows: The last step is to find the minimum probability level at which the parameter Z is significant, this could be found using two-tailed t statistical Tables or as mentioned by Abramowitz and Stegun (1972): where αmin: Minimum level of significance; q: counter; bx: constants: b0= 0.3989, b1= 0.3194, b2= -0.3566, b3= 1.7814, b4= -1.8213, b5= 1.3303, b6= 0.2316, ABS(Z): the absolute value of Z. Kendall tau is considered significant when alpha min is less than a specified alpha value, i.e 0.05.

Sen-slope estimator test
Sen's statistic is the median slope of each point-pair slope in a dataset (Sen, 1968).To perform the complete Sen's test, several rules and conditions should be satisfied; the time series should be equally spaced, i.e. the interval between data points should be equal.However, Sen's method considers missing data.The data should be sorted ascending according to time, and then apply the following formula to calculate Sen's slope estimator (Q) as the median of Sen's matrix members.Median Its sign reflect the trend's direction, while its value reflects how steep the trend is.To determine whether the median slope is statistically different than zero, the variance is calculated using Eqn.( 13), to obtain the confidence interval of Q at a specific probability level, e.g 95%.The area (Z) under two-tailed normal distribution curve is calculated at the level (1-α/2), where α=1-confidence level.For example, for a confidence level of 95%, Z should be evaluated at 0.975, hence Z= 1.96.Next, the parameter Cα is calculated as follows: The upper and lower confidence boundaries for Q are then calculated as follows: where int() represents the integer value; Mu and Ml are the upper and lower boundaries for Q at 1-α probability level; nq is the number of Sen's matrix members calculated from Equation ( 17), equal to nq=n(n-1).The median slope is then defined as statistically different from zero for the selected confidence interval if the zero does not lie between the upper and lower confidence limits.

Results and discussion
The daily ETo data for each of the studied stations were calculated for the study period.Then we summarize the data on monthly basis to find the average, maximum and minimum values per month for the whole country, Figure 3.The lowest values of ETo occurred in the winter season, December and January whereas the highest values occurred during summer months; June, July and August.The results showed high variation in ETo from about 5 mm/day to 15 mm/day in July whereas the absolute minimum and maximum ETo showed even higher variation and ranged from 3.9 mm/day in January to as high as 18.5 mm/day in July.However, Figure 2 indicates the high variability of climate conditions over Saudi

Regression analysis between
ETo and time has confirmed the positive trend and showed that the slope of the line was 0.020 with R 2 = 0.50.However, this relationship was not statistically significant at 95% probability level.Longer period of data analysis is needed to confirm this result.Nevertheless, the present analysis indicates clearly that climate variability is indeed affecting the country and the evapotranspiration demand is increasing with time.
Due to some restrictions in Man-Kendall and Sen's methods, two stations out of the 29 stations were omitted from calculations due to the small number of years they had (less than 10 years); those were stations #3 (Guraiat) and # 13 (Dammam).Mann-Kendall and Sen Slope statistics were performed on the rest 27 stations on monthly basis to confirm trends direction and test its significance.Two parameters were calculated namely Kendall  and Sen Slope Q and their confidence limits at 95% and 99% probability level as described in Materials and Methods.A group of selected results is shown in Figure 6 where the parameters of Mann-Kendall and Sen Slope and their significant tests are presented.The Figure represents ET trends in January for four stations, Tabuk, Sharurrah Yenbo, and Hail, showing possible combinations of Mann-Kendall (MK) and Sen Statistic, in addition to their significance under increasing or decreasing ETo conditions.That is, Figure 6a showing a downtrend with MK and Sen significant at 95% and 99%. Figure 6b showing a downtrend with only MK is significant at 95%. Figure 6c showing an uptrend with all statistics was significant.Figure 6d showing an uptrend with only MK is significant at 95%.The two tests gave similar results in all of these cases but Sen Slope test were found to be more conservative.A positive sign in τ or Q indicates an increasing trend, Figure 6 C and D while a negative value indicates a decreasing trend, Figure 6 A and B. The significance of τ was tested by comparing the calculated min with 0.05 or 0.01 for 95% and 99% confidence level, respectively; min < 0.05 or 0.01.The corresponding significant test for Q was carried out by calculating its confidence intervals at 95% and 99% indicated by respectively.If the two limits have similar sign, then the calculated Sen Slope Q value is confirmed not to be zero and therefore the slope is significantly different from zero, indicating a positive or negative trend for +Q or -Q, respectively.Sharurrah, both Q (-0.020, -0.028) and  (-0.365, -0.281), respectively, were negative indicating a decreasing trend for ETo.However, this decrease is significant for Tabuk and not significant for Sharurrah at 99% probability level as confirmed by both statistical methods.In the case of Tabuk , the value of calculated min = 0.0029 is less than 0.01 indicating a significant trend according to Mann-Kendall test; for Sen slope test the 99% minimum and maximum values of the confidence interval are both negative indicating that the slope Q is not zero and therefore a negative trend is confirmed.However, for Sharurrah stations min = 0.03 is less than 0.05 but larger than 0.01 indicating that the positive trend is significant at 95% but not at 99% level according to Mann-kendall test.A slightly different result was found in Sharurrah with Sen slope test.The upper limits for the Q confidence interval at 95% and 99% were both larger than zero, 95% max 0.004, that the median slope of the ETo series, Q, can actually be zero and therefore the negative trend is not significant at both levels.This results show that Sen Slope test can be more conservative than Mann-Kendall test.
Both Yenbo and Hail have increasing ETo trends, as shown in Figure 6C and D, because  and Q values were both positive (Q= 0.022, = 0.515 for Yenbo; while Q= 0.01 and  = 0.24 for Hail).Both tests showed that ETo in Yenbo is increasing significantly at 95% as well as at 99% level since τ and Q were positive and min = 0.000 < 0.05 and 95% min Q > 0.0.However, in Hail this uptrend was significant according to Mann-Kendall (min = 0.000 < 0.05) at 95% level but not significant according to Sen Slope since 95% min 0.001 Q  indicating the possibility of the slope being zero.These results indicated the validity of these two statistical methods to detect trends in a time series data.
The previous analyses shown in Figure 6 were carried out for the 27 stations out of 29 under study.Man-Kendall, and Sen's methods' can deal with data series with 10 or more data points.However, Gurrayat and Dammam have less than 10 years of data and they were excluded from trend analysis.Average ETo time series were analyzed for each month and the resultant Sen slope Q and Kendall  are shown in Table 2 and Table 3, respectively followed by up or down arrows to indicate their significance.Up-arrows in light or dark black indicate significance up trend while similar light and dark black down arrows indicate a significant down trends at 95% and 99% probability level, respectively.Numbers without arrows are not statistically significant.The total number of stations with a decreasing or increasing trend in each month were calculated and shown at the bottom of the Table .whereas the number of months at which stations showed a decreasing or increasing trends were shown for each station at the right side of the Table .Numbers between brackets indicates the number of months or stations with the corresponding significant trend.
The tests were carried out for maximum, minimum and average monthly ETo but only the average ETo is shown in the Tables, because extreme ETo showed similar behavior to that of average ETo.Fourteen stations have a positive Q and , for at least 10 months in a year therefore an uptrend in ETo namely; Turaif, Arar, Al jouf, Hafr Al-Baten, Hail, Gassim, Dhahran, Riyadh (old), Yenbo, Taif, W-dawaser, Bisha, Abha, and khamis Mushait.Another six stations showed a negative or zero Q during the whole year, therefore a downtrend in EToincluding; Tabuk, Wejh, Makkah, Nejran, Sharurrah, and Gizan.The other seven stations showed a mix of increasing and decreasing trends during the year and those are Rafha, Qaisumah, Ahsa, Madina, Riyadh (new), Jeddah, and Baha.
However, the up or down trends or downtrends in ETo in the first mentioned group were not always significant as indicated by the upward arrows and summed in the last two columns of Tables 2 and 3.Only Yenbo had a confirmed significant trend at 95% level during the entire year.Other stations showed a significant up trends in ETo for several months during the year including Hail and Taif , 10 months; Turaif and Arar, 9 months; Al jouf and Khamis Mushait, 8 months.The other stations among the uptrend group had significant uptrend in 3 months to 7 months in a year as shown in Table 2.
The number of stations with a decreasing trend is far less than those with increasing ETo.Few stations showed a decreasing trend in ETo for 9 months or higher including, Qaisumah, Tabuk, Wejh, Ahsa, Makkah, Nijran and Sharurrah.However, only Ahsa station had a significant decreasing trend for 7 months followed by Qaisumah and Sharurrah, 4 months, and Tabuk with only 3 months of declining in ETo.The rest of stations, Wejh, Makkah, Nejran, Gizan, had a decreasing trend but this trend is not significant at 95% probability level.
The numbers of stations with increasing or decreasing trend are shown in the last two rows of Table 2 and 3. Figure 7, as well, illustrates the number of stations with significant/nonsignificant increasing/decreasing trend of ET throughout the studied areas.At least 15 stations or higher showed an increasing trend for the entire year except in January at which 14 stations showed a decreasing trend.March, April and June showed the highest number of stations with increasing ETo.However, the significant increase in ETo were confirmed for about 10 stations and for 9 months, February to October.During the months of October to January about 10 stations showed a decreasing trend but this decrease was significant for only 4 stations in September and October, one station in November and December and 3 stations in January.
Further inspection on the location of stations with increasing trends in ETo revealed that most of these stations are located in the northern part of the Arabian peninsula north to the latitude line of 22 degrees.However, some other stations were located southern of this line at the southern west corner of Saudi Arabia.It seemed that stations located along the longitudinal line of 45 degrees showed an increasing trend.Actually the wind direction over the Arabian peninsula seemed to follow this line from south west to north in rainy seasons and from north to south west in the dry seasons.
To have an aerial graph for the regions with a decreasing or increasing trend in ETo a contour map were plotted for each months and the results are shown in Figure 8. Certainly, regions with increasing trends are concentrated in the northern part of Saudi Arabia and extended to the south along the 45-degree longitudinal line.Significant and increasing regions, indicated by black and grey regions (P>95%) are prevail for most of the year except in January and July to some extent.In January most of SA areas have decreasing ETo as shown in Figure 8 and also in Table 3 where 14 stations have a decreasing trends; although this decrease is not significant at 95% level except for 3 stations.The southeastern parts seemed to have decreasing trends most of the years, but also this trend is not significant at 95% except in July and October as indicated by the darker dotted regions.

Conclusion
Water scarcity problem can be solved by proper management of water usage.Most of the depleted water in KSA is consumed through agriculture.Identifying the ETo trend and knowing the zones having the least ETo values can help in determining the future plans of agricultural and water extensions.Historical analysis of daily ETo in Saudi Arabia was carried out using Penman Monteith equation (FAO-56) for 29 meteorological stations distributed all over Saudi Arabia for the period 1980 to 2008.The long time average daily ETo varied from about 5 mm/d in Jan to 15 mm/day in July which is one of the hottest months in the country.ETo time series analysis using Mann-Kendall and Sen slope statistics revealed that ETo has been increasing steadily during the study period.The average minimum and maximum daily ETo increased steadily and ETo average increased from about 9.6 to about 10.4 mm/day in 2008.Trend analysis revealed that about 14 of the weather stations showed a significant increasing trend in ETo during the year for more than 7 months.Only 4 stations showed decreasing trends in three months, September, October and January.
Increasing ETo trends prevail in the northern and south-west areas along the longitudinal line of 45 degrees while decreasing trends prevail in the north western spot along the red sea and south eastern parts along the Arabian Gulf.This demonstrates that ETo fluctuation is increasing with time that can be considered a significant sign for climate variability in the Arabian peninsula.This increase in ETo seemed to be mainly affected by the global warming or the increase in temperature in the Arabian peninsula which was confirmed by several studies mentioned in this paper.Analyses of longer historic data are needed to confirm these findings.This demonstrates that ETo fluctuation is increasing with time that can be considered a significant sign for climate change.Though, the findings of this research suggest the needs to consider ETo changes in the planning for agricultural and water resources projects.Thus to rank the areas with fixed and decreasing ETo trend as highly recommended zones for future agricultural projects, and to do the opposite with the increasing ETo trends' zones.Finally, if the low ranked zones are essential due to other circumstances, then the water management policy should consider the increment rate in ETo and its effect on water consumption.
[mm day-1], Rn is the net radiation at the crop surface [MJ m-2 day-1], G is the soil heat flux density [MJ m-2 day-1], Ta is the mean daily air temperature at 2 m height [°C], U2 is the wind speed at 2 m height [m s-1], es: saturation vapor pressure [kPa], ea: actual vapor pressure [kPa], ∆ is the slope of vapor pressure curve [kPa °C-1 ], and γ is the psychometric constant [kPa °C-1 ].

Figure 2 .
Figure 2. Sample contour map of daily and annual evapotranspiration in the KSA.(a): average daily evapotranspiration for the month of June over 30 years period (mm/d).(b): annual evapotranspiration of the year 1991 (mm/y)

Figure 3 .
Figure 3.The monthly average ETo variations during the year in KSA, vertical bars indicate standard deviation values.

Figure 4 .
Figure 4. Relationship between the average annual ETo (mm/d) as related to station altitudes, latitudes, and longitudes.

Figure 5 .
Figure 5. Temporal change of the average, maximum, and minimum ETo through the study period.

Figure 6 .
Figure 6.Analysis of monthly average ETo trends using Mann-Kendall and Sen slope estimator and their significant tests parameters: A: Tabuk area, B: Sharurrah area, C: Yenbo area, and D: Hail area.

Figure
Figure 6A-D represents four possible cases of Q and τ and their significance.In Tabuk and

Figure 7 .
Figure 7. Number of stations with significant/non-significant increasing/decreasing trends of ET in the studied areas.

Figure 8 .
Figure 8. Contour plots showing the distribution ETo trends' direction and confidence level over the area of KSA.

Table 1 .
All of the data ends in 2008 and started at 1980 and 1985 for the 29 and 24 years logging.Effect of Climate Change on Spatio-Temporal Variability and Trends of Evapotranspiration, and Its Impact on Water Resources Management in The Kingdom of Saudi Arabia 279 *: Stations having less than 10 years of data.**:averageerrorratio in data recording is less than 0.7% including missing records if any.Table1.Geographical information of the meteorological stations included in this study.

Table 2 .
Monthly trends in the average ETo as estimated from Sen's slope statistics (Q) for various meteorological Station.Numbers in parenthesis indicate number of significant values.

Table 3 .
Monthly trends in the Average ETo as estimated from Mann-kendall statistics () for various meteorological Station.Numbers in parenthesis indicate number of significant values.