Trends of Climate Change in Saudi Arabia : Implications on Water Resources

Climate change is an important factor for sustainable water resource management in the arid and semi-arid countries. In this study, future trends of temperature and rainfall were assessed for several regions in Saudi Arabia. The linear and Mann–Kendall analyses showed an increase of temperature in all regions and decrease of rainfall in many regions. Following trend analysis, the outputs of the NCAR Community Climate System Model were obtained for three emission scenarios (high: representative concentration pathways RCP8.5; high medium: RCP6; and low: RCP2.6) for the assessment periods of 2025–2044, 2045–2064 and 2065–2084 respectively, and compared with the average values from the reference period (1986–2005). In all emission scenarios, temperature showed an increase from 1986 to 2005 in all regions. For RCP8.5, increase of temperature are in the ranges of 0.8–1.6 ◦C, 0.9–2.7 ◦C and 0.7–4.1 ◦C during 2025–2044, 2045–2064 and 2065–2084 respectively. However, rainfall showed variable patterns with respect to emission scenarios and assessment periods. In most regions, the RCP6 showed decrease in rainfall from the reference period while the RCP8.5 and RCP2.6 showed variable patterns. The increase of temperature and variable pattern of rainfall may increase uncertainty in developing sustainable water resource management strategies.


Introduction
Climate change is an essential component for strategic water resource management in arid and semi-arid countries, including Saudi Arabia.Using regional models, Hasanean and Almazroui [1] showed the decreasing trends of rainfall in most regions during 1978-2009 while no trend was reported for the southwest region.Chowdhury and Al-Zahrani [2] predicted an increase of rainfall by 15-25 mm/year in the central, western and eastern regions by 2050 while Al-Zawad [3] predicted the rainfall increase by 26-35 mm/year during 2070-2100.In the southwest region, Chowdhury and Al-Zahrani [2] predicted the rainfall increase of 109.7-130.4mm/year while Al-Zawad [3] predicted an increase of 96.7 mm/year.
Using the linear and Mann-Kendall (M-K) trends, Rehman and Al-Hadrami [4] reported an extreme temperature trend in the west coast of Saudi Arabia.This study showed an increase in summer temperature and the number of hot days per year, due possibly to the combined effects of higher temperature and lower rainfall in this period.The mean annual maximum temperature could increase at a higher rate than the mean annual minimum temperature [4].Almazroui [5] reported the increase of temperature by 0.65 • C per decade while the central region to the coast of the Red Sea would be affected by the increased extreme rainfall events.This study applied the regional climate model (PRECIS) for the predictions.Chowdhury and Al-Zahrani [2] showed an increase of temperature by 1.8-4.1 • C in different regions of Saudi Arabia, which was consistent to the global positive trends [6].Addisu et al. [7] reported the extremely variable temperature and rainfall in 17 (out of 109) stations in Ethiopia where the coefficients of variation (CV) were in the range of 20-89%.The multi-model ensemble combination showed an enhanced predictive capacity for climatic parameters in comparison to a single model [8].
Several studies have reported the effects of climate change on the availability and quality of water resources, increased morbidity and mortality due to thermal extremes, epidemics, malnutrition, geographical and seasonal spread of infectious vector, rodent borne diseases, and cardiovascular and respiratory illness [9][10][11][12].In the past three decades, climate change related diseases were responsible for approximately 150,000 deaths and five million disability-adjusted life years (DALYs) [10,11].The Intergovernmental Panel on Climate Change (IPCC) has explained the vulnerability of climate change in terms of exposure to hazards and the capacity to adapt.Understanding the interaction of climatic parameters and environment is likely to reduce vulnerability, such as growing new crops suitable to climate change and efficient management of water resources [13,14].Increase in temperature might decrease dissolved oxygen (DO) and increase dissolved organic matter (DOM) and salinity in source water [2,15].Past studies have reported an increase of agricultural water demand by 5-15% during 2050, which was attributed to increased evapotranspiration [2].An increase of temperature by 1 • C might change the thermal limits of a crop by 10-30%, which could push the crop over the brink and affect crop yields [16].In addition, an increase in temperature by 1 • C might increase the capacity of air to hold water vapor by seven percent, resulting in reduced precipitation [17].As such, the water resources and soil-water balance might be affected further.Global warming will affect rainfall and temperature, which can have effects on groundwater reserves and soil moisture [18].Many countries have incorporated the effects of climate change for strategic water resource management.
Saudi Arabia is characterized by a semi-arid environment with high temperature variability, low annual rainfall, no natural perennial flow and limited groundwater reserves [2].The policies on agriculture, industry and water resources are greatly affected by the climatic condition.Despite the availability of several studies on climate forecasting, comprehensive information on climatic variability and its implications is limited.This study aims to use the historical data on climatic parameters from 1985-2013 to predict the current trends for temperature and rainfall in the central, north and southwest regions of Saudi Arabia.Following trend analysis, the outputs from 14 weather stations in Saudi Arabia were obtained from the National Center for Atmospheric Research (NCAR) Community Climate System Model (CCSM4) for the periods of 2025-2044, 2045-2064 and 2065-2084 under three standard emission scenarios, known as the representative concentration pathways (RCP).The RCP were defined as: high (RCP8.5),high-medium (RCP6) and low (RCP2.6)emission scenarios.The outputs were compared with temperature and rainfall data from the reference period (1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005), and the change in temperature and rainfall were predicted for the periods of 2025-2044, 2045-2064 and 2065-2084 respectively.The possible implications of climatic variability on water resources were discussed.

Data Statistics
In this study, the monthly data on temperature and rainfall in different regions of Saudi Arabia were obtained for the period of 1985-2013 from the Presidency of Meteorology and Environment (PME), Saudi Arabia [19].The data from Riyadh, Tabouk and Abha weather stations, representing the central, north and southwest regions of Saudi Arabia, respectively, were selected for analysis [19].These data were grouped into different seasons and specific months representing the cold and hot periods.The variability of data was discussed and the extreme values (e.g., outliers) were identified using the box plots.The coefficients of variation (CV) were calculated as CV = std.dev/mean.These data were used for trend analysis using the linear and Mann-Kendall methods.

Trend Analysis
Among many methods, the Mann-Kendall (M-K) test is widely used for trend analysis [19].The M-K test is advantageous due to its non-sensitivity to data normality (e.g., non-parametric) and the test statistics are not affected by the non-homogeneous time series.The M-K test assumes the hypothesis, H 0 , indicating no trend and that the observations (x i ) are randomly variable with time.The alternative hypothesis H 1 , indicating an increasing or decreasing monotonic trend.The test assumes time series of n records and x j and x k as two subsets of data where k = 1, 2, 3, . . ., n − 1; and j = k +1, k + 2, k + 3, . . ., n.The data are evaluated as a time series in sequence.Each value in the dataset is compared with all of the subsequent data points.If the value in the later period is higher than the earlier period, the M-K statistics, S is decremented by 1.If it is lower, S will be decremented by −1.Following the M-K test, final value of S can be computed as: where x j and x k are used as the annual or monthly values and/or j and k, j > k r respectively.A high positive value of S indicates an increasing trend while a low negative value indicates a decreasing trend.If n ≥ 10, S is considered normally distributed with variance as: where n = number of data points; q = number of tied groups; and t p = number data values in group p.Consequently, the standardized Z statistics follows the normal standard distribution as: The trend is computed using the Z value.The positive and negative values of Z indicate the increasing and decreasing trends respectively.For the non-parametric test, the Kendall's tau(τ) assess the statistical associations based on the ranks of the data.The Kendall's tau for two sets of data can be determined as: The value of τ varies between −1 and +1 representing a strong negative association to strong positive association between the data sets.
The linear and sequential version of M-K trends (U(t) and Û(t)) were investigated.The U(t) and Û(t) demonstrate the change of trend with time.This test considers the relative values of all terms in a time series (x 1 , x 2 , . . ., x n ).The magnitudes of x j annual mean time series (j = 1, 2, 3, . . ., n) are compared with x k (k = 1, 2, . . ., j − 1).At each comparison, the number of cases x j > x k is counted and denoted by n j .The test statistics (t) are computed as: The mean and variance of t is calculated as: 7) The sequential values of the statistics (Û) are calculated as Similarly, values of Û(t) are computed backward starting from the end of the series.The sequential version of the Mann-Kendall could be considered as an effective way of locating the beginning year(s) of a trend.The M-K statistics, U(t) is a value that indicates direction (or sign) and statistical magnitude of a trend in a series.When U(t) is significant at the 5% level, it can be decided whether it is an increasing or a decreasing trend depending on whether U(t) > 0.0 or U(t) < 0.0.

NCAR Community Climate System Model (CCSM4) Data
The Community Climate Model (CCM) was developed by the National Center for Atmospheric Research (NCAR) in 1983.The CCM has been improved to the Community Climate System Model (CCSM), which is a coupled climate model system for simulating the Earth's climate.In CCSM, the earth's atmosphere, ocean, land surface and sea-ice, and one central coupler component are simulated concurrently, which assists in explaining the earth's past, present and future climate.It allows the downscaling of climatic information following the dynamic and statistical approaches [20,21].The details on downscaling can be found in [21].
In 2008, a total of 20 modeling groups and the World Climate Research Programs (WCRP) Working Group on Coupled Modelling agreed to promote a new set of coordinated climate model experiments with the inputs from the Analysis, Integration and Modeling of the Earth System (AIMES) project of the International Geosphere-Biosphere Programme (IGBP).These experiments represented the fifth phase of the Couple Model Inter-comparison Project (CMIP5).The CMIP5 was involved in the modeling of long-term (century scale) integration and decadal prediction, or near-term integration (10-30 years) [22][23][24].In CMIP5, the long-term simulation was coupled with the biogeochemical components that accounted for the fluxes of carbon between the ocean, atmosphere and terrestrial biosphere carbon reservoirs.The spatial resolution for the atmospheric and ocean components were in the ranges of 0.5-4.0• and 0.2-2.0• respectively.The CCSM4 model was simulated for the RCP8.5, RCP6 and RCP2.6 emission scenarios.The output data on temperature and rainfall were recorded.The data on temperature and rainfall in the central, north, southwest, eastern and western regions of Saudi Arabia were analyzed for the periods of 2025-2044, 2045-2064 and 2065-2084 respectively.These data were compared with the reference period (1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005).The weather stations were located in five regions: central (Riyadh, Unaizah); north (Hail, Al-Jouf and Tabouk); southwest (Abha, Al-Baha and Jazan); eastern (Dammam and Al-Hofuf); and western (Jeddah, Taif, Makkah and Al-Madinah) respectively [19].Further details on the model and outputs can be found in Figure 1 and from the NCAR Community Climate System Model (CCSM4) web portal.

Central Region
The seasonal variability of temperature in Riyadh, Tabouk and Abha representing the central, north and southwest regions are shown in Table 1.The data are presented for mean temperature in  four seasons: winter (December, January and February), spring (March, April and May), summer (June, July and August) and autumn (September, October and November).In addition, the data in January and June are summarized to explain the cold and hot months respectively (Table 1).The mean temperature in summer was 2.4 times the mean temperature in winter (35.5 • C and 15 • C respectively), which is characterized as continental climate.The mean temperature in spring and autumn were similar (26.2 • C and 26.6 • C respectively).The mean and median were similar in each season while skewness and kurtosis were small (Table 1).Skewness measures of the lack of symmetry of data while kurtosis indicates whether the data are heavy-tailed or light-tailed relative to a normal distribution.The negative and positive skewness indicated that the data were skewed toward left and right respectively (Table 1).In addition, CV were low (range: 0.020-0.107).The similar mean, median, low skewness, kurtosis and CV indicated the normal distributions (Table 1).The box plots (Figure 2) show the percentile distributions of data where the bottom and top of the boxes represent the first quartile (25-percentile) and the third quartile (75-percentile) respectively.The lower and upper fringes represent the lowest and highest limits of the data respectively.Any data beyond this limit is an outlier, such as 12.6 • C for winter (Figure 2a), representing an extreme low mean temperature.

North Region
The summer mean temperature (31 • C) was almost 2.5 times the winter mean temperature (12.2 • C), indicating the continental climate (Table 1).The skewness were small (range: −0.38-0.54)while standard deviation (range: 0.69-1.6)and CV (range: 0.023-0.147)were low (Table 1), indicating the normal distributions.Figure 2e-h shows that the fourth quartile temperature in winter were in the range of 13-13.2• C (e.g., highest) while the first quartile temperature were in the range of 9.8-11.6• C (Figure 2e).Two outliers for the highest (15.2 • C) and lowest (9 fourth quartile temperature was in the range of 31.5-32.2• C with the maximum of 33 • C (Figure 2g).In autumn, approximately 50% data were in the range of 23-23.9• C (Figure 2h).

North Region
The summer mean temperature (31 °C) was almost 2.5 times the winter mean temperature (12.2 °C), indicating the continental climate (Table 1).The skewness were small (range: −0.38-0.54)while standard deviation (range: 0.69-1.6)and CV (range: 0.023-0.147)were low (Table 1), indicating the normal distributions.Figure 2e-h shows that the fourth quartile temperature in winter were in the

Southwest Region
The southwest region was different from the central and north regions.In the southwest region, mean values in summer and winter were 24.3 • C and 19.4 • C (Table 1) respectively.The intra-seasonal difference was less than 5 • C while in the central and north regions, the differences were up to 20 • C (Table 1).The wet summer was likely to reduce the temperature difference in the southwest region while the other regions had very low rainfall in summer.The data had low skewness and low CV (range: 0.062-0.119),and the mean and median were similar in most seasons (except spring) (Table 1), indicating the normal distributions of data.In summer, the range of the first to third quartiles (50% data) were 22.9 • C to 25.4 • C (Figure 2k) while in winter, the range was 12.3-13.9• C (Figure 2i).The winter temperature in the southwest was higher than the central and north regions, due mainly to the effects of the moist air mass from the Red Sea in the equatorial side.The central and north regions were affected by the dry continental air mass from the Siberian high system in winter.Overall temperature was lower in the southwest region due to the effects of elevation, while the other regions were affected by tropical air mass in spring, summer and autumn.Further details on temperature distributions can be found in Figure 2i-l.

Rainfall
The rainfall showed spatial and temporal variability, which was irregular, torrential and often associated with thunderstorms [1].In the southwest, rainfall was often affected by the wet summer monsoon.The mean rainfall in the central, north and southwest regions were 106.9, 31.7 and 206 mm/year respectively and the ranges were 32-309.4,1.1-89.2and 87.4-568.2mm/year respectively (Table 2).The lower rainfall in the central and north regions indicated the aridity in these regions (Table 2).The skewness were in the range of 0.84-1.6indicating the distributions with right tails.The higher values of skewness and CV indicated non-normal distributions of rainfall (Table 2).Figure 3 shows that 50% data in the central region were in the range of 65-111 mm/year.The fourth quartile data were in the range of 111-168 mm/year and the outlier was 309 mm/year.In the north region, 50% data had the range of 11-45 mm/year and the fourth quartile data were in the range of 45-89 mm/year.In the southwest, 50% data had the range of 125-252 mm/year while the fourth quartile data were in the range of 252-413 mm/year.The extreme value (568 mm/year) was an outlier (Figure 3).2e).Two outliers for the highest (15.2 °C) and lowest (9 °C) were noted.In summer, the fourth quartile temperature was in the range of 31.5-32.2°C with the maximum of 33 °C (Figure 2g).In autumn, approximately 50% data were in the range of 23-23.9°C (Figure 2h).

Southwest Region
The southwest region was different from the central and north regions.In the southwest region, mean values in summer and winter were 24.3 °C and 19.4 °C (Table 1) respectively.The intra-seasonal difference was less than 5 °C while in the central and north regions, the differences were up to 20 °C (Table 1).The wet summer was likely to reduce the temperature difference in the southwest region while the other regions had very low rainfall in summer.The data had low skewness and low CV (range: 0.062-0.119),and the mean and median were similar in most seasons (except spring) (Table 1), indicating the normal distributions of data.In summer, the range of the first to third quartiles (50% data) were 22.9 °C to 25.4 °C (Figure 2k) while in winter, the range was 12.3-13.9°C (Figure 2i).The winter temperature in the southwest was higher than the central and north regions, due mainly to the effects of the moist air mass from the Red Sea in the equatorial side.The central and north regions were affected by the dry continental air mass from the Siberian high system in winter.Overall temperature was lower in the southwest region due to the effects of elevation, while the other regions were affected by tropical air mass in spring, summer and autumn.Further details on temperature distributions can be found in Figure 2i-l.

Rainfall
The rainfall showed spatial and temporal variability, which was irregular, torrential and often associated with thunderstorms [1].In the southwest, rainfall was often affected by the wet summer monsoon.The mean rainfall in the central, north and southwest regions were 106.9, 31.7 and 206 mm/year respectively and the ranges were 32-309.4,1.1-89.2and 87.4-568.2mm/year respectively (Table 2).The lower rainfall in the central and north regions indicated the aridity in these regions (Table 2).The skewness were in the range of 0.84-1.6indicating the distributions with right tails.The higher values of skewness and CV indicated non-normal distributions of rainfall (Table 2).Figure 3 shows that 50% data in the central region were in the range of 65-111 mm/year.The fourth quartile data were in the range of 111-168 mm/year and the outlier was 309 mm/year.In the north region, 50% data had the range of 11-45 mm/year and the fourth quartile data were in the range of 45-89 mm/year.In the southwest, 50% data had the range of 125-252 mm/year while the fourth quartile data were in the range of 252-413 mm/year.The extreme value (568 mm/year) was an outlier (Figure 3).The M-K trends were analyzed for annual mean temperature, and minimum and maximum temperature in January and June for Riyadh, Tabouk and Abha representing the central, north and southwest regions respectively.The parameters are shown in Table 3.The p-value ≤ 0.05 indicated the support toward rejection of H 0 and acceptance of H 1 , indicating the presence of trend.The positive values of Kendall's tau(τ) showed an increasing trend while the negative values showed a decreasing trend.The high positive value of S statistic indicated an increasing trend while the low negative value indicated a decreasing trend.In Riyadh, temperature showed increasing trends in most scenarios except the maximum temperature in January (Table 3).In Tabouk, the M-K test showed increasing trends for annual mean temperature and minimum temperature in June, while no trend was observed for the other data.In Abha, the M-K test showed increasing trends in all cases except the maximum temperature in January (Table 3).The linear trends for the winter and summer mean temperature are shown Figure 4.In Riyadh, Tabouk and Abha, the data showed increasing trends in both winter and summer seasons (Figure 4).In all regions, the rates of increase of temperature were higher in summer than winter.In Riyadh, averages of temperature in summer and winter during 1985-1989 were 33.54 • C and 12.96 • C respectively.During 2009-2013, the averages were 35.4 • C and 14.6 • C respectively.During this period, increase of temperature in summer and winter were 1.86 • C and 1.64 • C respectively.The linear rates of increase during summer and winter were 0.0676 • C/year and 0.0561 • C/year respectively (Figure 4).In Tabouk, averages of temperature in summer and winter during 1984-1988 were 30.6 • C and 12.6 • C respectively.During 2009-2013, the averages were 31.9• C and 13.4 • C respectively.The linear rates of increase of temperature during summer and winter were 0.0583 • C/year and 0.0427 • C/year respectively (Figure 4).In Abha, averages of temperature in summer and winter during 1975-1979 were 21.9 • C and 13.2 • C respectively.During 2005-2009, the averages were 25.5 • C and 14.3 • C respectively.The linear rates of increase of temperature during summer and winter were 0.125 • C/year and 0.0662 • C/year respectively.Overall, the rates of increase of temperature in summer and winter were in the ranges of 0.0583-0.125 • C/year and 0.0427-0.0662• C/year respectively.The higher rates of increase of summer temperature were consistent to past study [4].Almazroui [5] reported that the average rate of increase of temperature was 0.65 • C per decade, which was consistent to the current findings.The higher rates of temperature increase during summer deserve further attention.
period, increase of temperature in summer and winter were 1.86 °C and 1.64 °C respectively.The linear rates of increase during summer and winter were 0.0676 °C/year and 0.0561 °C/year respectively (Figure 4).In Tabouk, averages of temperature in summer and winter during 1984-1988 were 30.6 °C and 12.6 °C respectively.During 2009-2013, the averages were 31.9°C and 13.4 °C respectively.The linear rates of increase of temperature during summer and winter were 0.0583 °C/year and 0.0427 °C/year respectively (Figure 4).In Abha, averages of temperature in summer and winter during 1975-1979 were 21.9 °C and 13.2 °C respectively.During 2005-2009, the averages were 25.5 °C and 14.3 °C respectively.The linear rates of increase of temperature during summer and winter were 0.125 °C/year and 0.0662 °C/year respectively.Overall, the rates of increase of temperature in summer and winter were in the ranges of 0.0583-0.125 °C/year and 0.0427-0.0662°C/year respectively.The higher rates of increase of summer temperature were consistent to past study [4].Almazroui [5] reported that the average rate of increase of temperature was 0.65 °C per decade, which was consistent to the current findings.The higher rates of temperature increase during summer deserve further attention.

Rainfall
The future trends were assessed for the annual rainfall in Riyadh, Tabouk and Abha.The data did not show any trend in the central and north regions (Table 3) and were supported by the low Kendall's tau (−0.075 and −0.076 respectively).The southwest region showed decreasing trend (pvalue = 0.005), which was also supported by the large negative S (S = −157) and negative Kendall's tau (−0.361) (Table 3).The linear and sequential version of M-K trends (U and Û) for Riyadh, Tabouk and Abha are shown in Figure 5a-f.In Riyadh, no consistent trend was observed (Figure 5a

Rainfall
The future trends were assessed for the annual rainfall in Riyadh, Tabouk and Abha.The data did not show any trend in the central and north regions (Table 3) and were supported by the low Kendall's tau (−0.075 and −0.076 respectively).The southwest region showed decreasing trend (p-value = 0.005), which was also supported by the large negative S (S = −157) and negative Kendall's tau (−0.361) (Table 3).The linear and sequential version of M-K trends (U and Û) for Riyadh, Tabouk and Abha are shown in Figure 5a-f.In Riyadh, no consistent trend was observed (Figure 5a) with increasing trend during 1992-1997 and 2006-2010, and decreasing trends during 1985-1992 and 2010-2015 (Figure 5a,b).In Tabouk, no consistent trend was also observed with decreasing trend during 1985-2003 and increasing trend during 2004-2015 (Figure 5c,d).The rainfall in Abha (southwest region) showed an overall decreasing trend (Figure 5e).The sequential plot (Figure 5f) showed an increasing trend during 1988-1997and 2004-2005. However, during 1998-2003, 2006-2009and 2010-2012, rainfall showed decreasing trends.Overall, multiple trends were observed in all regions during the data periods.However, the overall rainfall in Abha has been reduced in the recent years and was also evidenced by the reduction in the natural vegetation [25].
Climate 2018, 6, 8 region) showed an overall decreasing trend (Figure 4e).The sequential plot (Figure 5f) showed an increasing trend during 1988-1997and 2004-2005. However, during 1998-2003, 2006-2009and 2010-2012, rainfall showed decreasing trends.Overall, multiple trends were observed in all regions during the data periods.However, the overall rainfall in Abha has been reduced in the recent years and was also evidenced by the reduction in the natural vegetation [25].
Overall, the highest increase of rainfall was predicted for Jazan (southwest) by 152.4 mm/year during 2065-2084 under RCP8.5.The RCP6 showed the decrease in rainfall in most stations during 2025-2044 and 2045-2064, and increase in 2065-2084 (Table 4).In RCP2.6, all stations showed increase in rainfall from the reference period except Abha (southwest) and Taif (western) during 2065-2084.Rainfall in Abha and Taif are likely to be reduced by 7.65 mm/year and 2.5 mm/year respectively.The data showed significant variability among the regions, which needs better understanding for water resource management.It is to be noted that the CCSM4 model was simulated based on the single valued parameters for the RCP8.5, RCP6 and RCP2.6 emission scenarios while the parameters were likely to have uncertainty.While the study presented the averages of temperature and rainfall under the low, medium and high emission scenarios, tempo-spatial variability of input parameters could impart significant uncertainty to the output data.As such, the forecasted temperature and rainfall are likely to have a range instead of single values.Upon availability of more reliable output data in future, the comparison with the reference period can be updated in future.

Implications on Water Resources and Quality
Three emission scenarios were investigated for temperature and rainfall change during 2025-2044, 2045-2064 and 2065-2084.The increase in temperature can vary in the ranges of 0.8-1.9• C, 0.9-3.7 • C and 0.7-4.1 • C respectively.The changes in rainfall can be in the ranges of −33-38.1,−38.1-111.8 and −7.6-152.4mm/year respectively.The increase in temperature and variable patterns of rainfall can have consequences on water resources, water quality, agricultural productions and agricultural water demand.
Increase in temperature is likely to increase evaporation from the open reservoirs.In Saudi Arabia, there are 449 dams to collect, store and recharge runoff and to control flash floods [26].These dams are often used as potential sources for drinking water.The southern part has the valley basins, which often dry up due to higher evaporation, excessive use of water and low precipitation [26][27][28].Past studies reported an approximately 50% increase in evaporation between 2011 and 2050 [2].Few studies have reported that the free surface evaporation from wadies in Saudi Arabia could be in the range of 5-80% of water on annual basis [25,27].For a dam reservoir in the southwestern region, loss through surface evaporation was reported to be in the range of 4.7-6.0m per year [28].Missimer et al. [28] reported that up to 80% of water could be lost in old wadies through evaporation.Increase in temperature is likely to increase the evaporative loss from the dam reservoirs and affect the water sources.Future study should investigate the extent of such effects.In addition, an increase in temperature by 1 • C and 5 • C may reduce surface runoff by 115-184 million m 3 (MCM) and 600-960 MCM per year respectively [29,30].The recharge of water to all aquifers was estimated to reduce by 91.4 MCM and 475 MCM per year due to 1 • C and 5 • C increase in temperature respectively.The overall reduction in water sources was estimated to be 241-1435 MCM per year [29,30].Other studies reported the possible increase in agricultural water demand by 5-15% in 2050 for the unchanged level of productions [2,[31][32][33].
Increase in temperature is likely to increase the domestic water demand.In Phoenix, Arizona, increase in temperature by 1 • C was reported to increase domestic water demand by 6.7% [34].Using a similar rate, domestic water demand in different regions were estimated to increase by 5.4-12.7%,6.0-24.8%and 4.7-27.5% during 2025-2044, 2045-2064 and 2065-2084 respectively.Upon availability of the effects of temperature on domestic water demand in Saudi Arabia, the assessment can be updated further.Future study should investigate the effects of temperature increase in Saudi Arabia.In addition, quality of water in the drinking water sources are likely to be degraded with the increase of temperature [2,30,31].
The outliers in Figure 3 demonstrate occasional intense rainfall events in the central, north and southwest regions.These rainfall events are likely to generate a significant amount of surface runoff, which can recharge the groundwater aquifers [32,33].In addition, the runoff can be collected and used for other purposes (e.g., domestic, agricultural uses).In Saudi Arabia, occasional intense rainfall events have caused several flash floods in the past, which were responsible for the loss of lives and properties [35][36][37].There is a need to better protect lives and properties during flash floods from the extreme rainfall events.Future study should assess the flood hazard, runoff collection potential and aquifer recharge from such events.The increase in temperature and variable pattern of rainfall are likely to increase uncertainty on the available water resources [37][38][39].Through a study in Scotland, Afzal et al. [39] demonstrated the decrease in the time-based and volumetric reliability of reservoirs under the assumption of an unchanging demand.The variability of rainfall had the greatest effect on reservoir reliability.In addition, increase in water demand could reduce reservoir reliability [39].In Saudi Arabia, domestic water demands are likely to increase in future while the demands are highly variable in Makkah and Al-Madinah, which are very high during four months (Arabic months: Ramadan through Dhu Al-Hijjah).The decline in reliability of water resources necessitates careful planning for water demand satisfaction under the highly variable demand characteristics in Saudi Arabia.

Conclusions
This study presented the trends of temperature and rainfall in the central, north and southwestern regions of Saudi Arabia based on the historical data.The change in temperature and rainfall were predicted for the periods of 2025-2044, 2045-2064 and 2065-2084 using the data from the National Center for Atmospheric Research (NCAR) Community Climate System Model (CCSM4).Temperature was predicted to increase in all regions throughout 2025-2084 with the highest increase in the central and north regions while rainfall showed variable patterns with significant increase in the southwest and western regions.The regional variability of temperature increase and rainfall are significant, indicating the need for region-specific action plans to minimize the effects of climate change.In the context of water resource management, inclusion of such variability is important.Further, consumption of water in different sectors (e.g., domestic, industrial and agriculture) are likely to be affected by climate change.The quality of water in the sources (e.g., groundwater, seawater used for desalination, etc.) can be affected, which may increase treatment costs.Future study is essential to explain the extent of such effects on groundwater resources, groundwater and seawater quality, and cost of water treatment and desalination.The implications of climate change on water treatment, supply cost and human health risks deserve further attention.Despite several limitations, this study sheds light on the possible changes of temperature and rainfall in different regions of Saudi Arabia and the implications on water resource management.

Figure 1 .
Figure 1.Locations of 14 weather stations in Saudi Arabia included in the current study.The arrows in red and blue colors show the wind directions in summer and winter respectively.Central region (Riyadh, Unaizah); North region (Hail, Al-Jouf and Tabouk); Southwest region (Abha, Al-Baha and Jazan); Eastern region (Dammam and Al-Hofuf); and Western region (Jeddah, Taif, Makkah and Al-Madinah).

Figure 1 .
Figure 1.Locations of 14 weather stations in Saudi Arabia included in the current study.The arrows in red and blue colors show the wind directions in summer and winter respectively.Central region (Riyadh, Unaizah); North region (Hail, Al-Jouf and Tabouk); Southwest region (Abha, Al-Baha and Jazan); Eastern region (Dammam and Al-Hofuf); and Western region (Jeddah, Taif, Makkah and Al-Madinah).

Figure 2 .
Figure 2. Box plots of mean temperature for different seasons in Riyadh, Tabouk and Abha representing the central, north and southwest regions of Saudi Arabia.(a-d): winter, spring, summer and autumn in Riyadh respectively.(e-h): winter, spring, summer and autumn in Tabouk respectively.(i-l): winter, spring, summer and autumn in Abha respectively.Bottom and top of the boxes represent the first quartile (Q1) and the third quartile (Q3) respectively.The lower whisker extends up to (Q1 − 1.5 × (Q3 − Q1)) and upper whisker extends up to (Q3 + 1.5 × (Q3 − Q1)).The horizontal lines within the box represent the medians of the data.Values larger than Q3 or smaller than Q1 by at least 1.5 times the interquartile range (IQR) are outliers.The (+) sign indicates the mean values of the data.

Figure 2 .
Figure 2. Box plots of mean temperature for different seasons in Riyadh, Tabouk and Abha representing the central, north and southwest regions of Saudi Arabia.(a-d) winter, spring, summer and autumn in Riyadh respectively.(e-h) winter, spring, summer and autumn in Tabouk respectively.(i-l) winter, spring, summer and autumn in Abha respectively.Bottom and top of the boxes represent the first quartile (Q1) and the third quartile (Q3) respectively.The lower whisker extends up to (Q1 − 1.5 × (Q3 − Q1)) and upper whisker extends up to (Q3 + 1.5 × (Q3 − Q1)).The horizontal lines within the box represent the medians of the data.Values larger than Q3 or smaller than Q1 by at least 1.5 times the interquartile range (IQR) are outliers.The (+) sign indicates the mean values of the data.

Figure 3 .Figure 3 .
Figure 3. Box plots of rainfall for Riyadh (a), Tabouk (b) and Abha (c) representing the central, north and southwest regions.Bottom and top of the boxes represent the first quartile (Q1) and the third quartile (Q3) respectively.The lower whisker extends up to (Q1 − 1.5 × (Q3 − Q1)) and upper whisker extends up to (Q3 + 1.5 × (Q3 − Q1)).The horizontal lines within the box represent the medians of the

Figure 5 .
Figure 5. Linear trends (left) and the sequential version of the Mann-Kendall trends (right) of the mean annual rainfall for the central (a,b), north (c,d) and southwest (e,f) regions (1985-2013).U(t) and Û(t) are the sequential M-K statistics, computed through forward and backward starting of the time series data.

Figure 5 .
Figure 5. Linear trends (left) and the sequential version of the Mann-Kendall trends (right) of the mean annual rainfall for the central (a,b), north (c,d) and southwest (e,f) regions (1985-2013).U(t) and Û(t) are the sequential M-K statistics, computed through forward and backward starting of the time series data.

Figure 6
Figure 6 shows the temperature increase from 1986 to 2005.The central region is represented by two stations: Riyadh and Unaizah (Figure 6a).Increase in temperature during 2025-2044, 2045-2064 and 2065-2084 are in the ranges of 1.2-1.6 • C, 1.3-2.6 • C and 1.2-4.1 • C respectively (Figure 6a).In RCP8.5, increase in temperature during 2065-2084 are in the range of 3.9-4.1 • C. In 2025-2044 and 2045-2064, these are 1.6 • C and 2.6 • C respectively (Figure6a).In RCP8.5 and RCP6, increase in temperature was minimum in 2025-2044 and maximum in 2065-2084 respectively.However, in RCP2.6, increase in temperature was maximum during 2045-2064.The variability under RCP2.6 may be due to the interaction effects of rainfall, temperature and the flow of cold/warm wind.

Table 1 .
The data are presented for mean temperature in

Table 1 .
Seasonal variability of mean temperature in the central (Riyadh), north (Tabouk) and southwest (Abha) regions of Saudi Arabia.
CV: Coefficient of variation.

Table 2 .
Summary of rainfall in the central (Riyadh), north (Tabouk) and southwest (Abha) regions of Saudi Arabia.°C (e.g., highest) while the first quartile temperature were in the range of 9.8-11.6°C (Figure

Table 3 .
Mann-Kendall trend analysis of temperature and rainfall for Riyadh, Tabouk and Abha representing the central, north and southwest regions of Saudi Arabia.