Drought Trends and Temperature Influence in Zhanghe River Basin, China

Our study area is one of the semiarid region of the China with under water stress condition that causes economic damage. The main objective of this study is to apply standardized precipitation evapotranspiration index (SPEI) and to use linear regression to calculate drought conditions in the study area. For this purpose, data from 1980 to 2010 was analyzed at different (1, 6, 12, and 24 months) time scales. Results depicted both dry and wet periods in the study area; occurrence of dry span with different frequency andmagnitudewas increased over last decades (2000–2010) atmost of the stations. Statistical results demonstrated that temperature was decreased in the 1st decade inmost of stations but in two decades from 1990 to 2000 and 2001 to 2010, temperaturewas increased except in Changzhi station. These results could be a future reference for developing information programs about monitoring and early drought information, planning of existing reservoirs, and management of water resources under climate conditions.


Introduction
All over the world, effects of climate change have been observed and effects of droughts are very common [1,2].Also, it is expected that as changes in the climate continue especially in regard to global warming, the potential for extreme periods of droughts will also be on the increase [3,4].Effects of climatic variations particularly drought in the agriculture can be measured by using various models and by analyzing drought indices.Such a kind of studies relies upon some variables like precipitation, temperature, solar radiation, wind velocity, and evapotranspiration.Mohammed H. I. Dore [5] studied global precipitation trends and found spatial and temporal changes in the precipitation.The basic reason of the drought is unbalanced relationship between precipitation and the required amount of water that creates worst effects on the agriculture such as decrease in yield production, economic losses due to forests loss, unbalance in ecosystem, and environmental loss such as soil erosion and desertification.In case of agricultural drought the availability of soil moisture for plant growth is a basic component, if evapotranspiration increases and water is insufficient in the soil and in these conditions water stress takes place and this stress over longer duration leads to droughts.Evapotranspiration is a main factor to assess water stress; several methods and techniques have developed to calculate evapotranspiration [6][7][8].Different scientists worked on different models such as [9][10][11][12][13][14][15][16][17][18][19][20] to improve these methods and models to measure evapotranspiration.Basically evapotranspiration is not a water stress indicator but should be considered as crop water stress index [21][22][23] as it is influenced by a combination of other meteorological variables like wind speed and solar radiation and is mostly used in agriculture for irrigation, soil water requirements, and forecasting of production loss due to shortage of water.Evapotranspiration deficit index and evaporative stress index should be used to estimate drought conditions [24,25].According to [26], standardized precipitation evapotranspiration is one of the comparable indexes to identify drought characters in arid regions due to warming.Increasing rates in evaporation leads to droughts [27].The occurrence, severity and duration of droughts are always varying that causes the greatest loss amongst all the natural disasters in the world.Considering the effect of precipitation on the occurrence of droughts, it is also important that evapotranspiration also be considered due to its rather direct effect on the amount of water that is available in the soil for the growing season [28].Zhanghe River basin is a semiarid region due to water stress conditions recently drought appears in that area more frequently.
However a comprehensive analysis on effects of global warming and variability in precipitation still needs more research.
The main objectives of this study are:

Study Area and Data
Collection.This study is based on nine weather stations located at upstream of Zhanghe River basin of China.Data spans over 30-year (1980 to 2010) period.The river has a catchment area of 18,200 km 2 .The watershed falls approximately between 112 ∘ and 115 ∘ east and 34 ∘ and 39 ∘ north.Figure 1 shows the locations of meteorological stations at Zhanghe River basin.

Methodology.
The SPEI is a very simple index to calculate and is based on the original SPI calculation procedure.The SPI is calculated using monthly precipitation as the input data.The SPEI uses the monthly (or weekly) difference between precipitation and PET.This represents a simple climatic water balance which is calculated at different time scales to obtain the SPEI.

SPEI Definition and Theory.
A number of equations exist to model PET based on available data for instance the Thornthwaite equation, the Penman-Monteith equation, the Hargreaves equation, and so forth and the SPEI is not linked to any particular one.We used the FAO-56 Penman-Monteith equation [29].
By default the original parameterization of [30] is used, corresponding to a short reference crop of 0.12.The method requires data on the incoming solar radiation and this is estimated from data on the bright sunshine duration  sun,  height.
Similarly, if data on the saturation water pressure are not available, it is possible to estimate it from the dewpoint temperature  dew , from the relative humidity RH or even from the minimum temperature  min (sorted from least to most uncertain method).
Similarly, the atmospheric surface pressure required for computing the psychometric constant can be calculated from the atmospheric pressure at sea level  and the elevation , or else it will be assumed to be constant (101.3 kPa).
With a value for potential evapotranspiration (PET), the difference between the precipitation () and PET for the month  is calculated: This (  ) provides a simple measure of the water surplus or deficit for the analyzed month.The calculated   values are aggregated at different time scales, following the same procedure as for the SPI.

Standardization of the Variable. Selection of the most
suitable statistical distribution to model the  series was difficult, given the similarity among the four distributions (Pearson III, lognormal, log-logistic, and general extreme value).The selection was based on the behavior at the most extreme values.Log-logistic distribution showed a gradual decrease in the curve for low values and coherent probabilities were obtained for very low values of , corresponding to 1 occurrence in 200 to 500 years.Additionally, no values were found below the origin parameter of the distribution.These results were used for the log-logistic distribution for standardizing  series to calculate the SPEI.
Parameters of the log-logistic distribution can be obtained following different procedures.Among them, the -moment procedure is the most robust and easy approach [31].
Calculate -moment's ratios ( skewness  3 and  kurtosis  4 ) as where  2 ,  3 , and  4 are the  moments of  series, calculated from the probability-weighted moments (PWMs) by using the following equations: The probability-weighted moments (PWMs) are obtained by means of the unbiased estimator given by [32].The unbiased PWMs are obtained according to where   is a frequency estimator and computed by using [33] where  is the range of observations and is arranged in increasing order and  is the number of data points.The probability density function of a three-parameter log-logistic distributed variable is computed as The probability density function of a three-parameter loglogistic distributed variable is expressed as follows: , , and  are scale, shape, and origin parameters, respectively, for  values in the range ( >  < ∞).When -moments are calculated, the parameters of the Pearson III distribution can be obtained in the following [34]: where Γ() is the gamma function of .
The probability distribution function of  according to the log-logistic distribution is then given by With () the SPEI can easily be obtained as the standardized values of ().For example, consider the classical approximation of [35] for  ≤ 0.5,  being the probability of exceeding a determined  value,  = 1−().If  > 0.5,  is replaced by 1− and the sign of the resultant SPEI is reversed.The constants are  0 = 2.515517,  1 = 0.802853,  2 = 0.010328,  1 = 1.432788,  2 = 0.189269, and  3 = 0.001308.The average value of the SPEI is 0, and the standard deviation is 1.The SPEI is a standardized variable, and it can therefore be compared with other SPEI values over time and space.An SPEI of 0 indicates a value corresponding to 50% of the cumulative probability of , according to a log-logistic distribution.
SPEI may be calculated for any month in a series of monthly precipitation data for a given location for the previous  months where  = 1, 2, 3, . . ., 12, . . ., 48, depending on the time series of interest.A 3-month SPEI makes use of the previous 2 months and the month being considered.A 12-month SPEI in the same way uses the last 11 months and the month being considered and as such, the SPEI is time dependent.
In this study, SPEI is determined for 6 time scales, 1, 3, 6, 9, 12, and 24 months.As such, short term through intermediate to long term seasonal analyses become possible.The SPI was determined using the SPEI R package version 1.6 by Santiago Beguerí.

Linear Regression Method.
Linear regression method is one of the methods which are used to compute slope.The slope indicates the mean temporal changes of the analyzed variable.Positive values of the slope show increasing trends, while negative values of the slope indicate decreasing trends.
A linear regression line has an equation of the form where  = the explanatory variable,  = independent variable,  = slope of the line, and  = intercept.Results also depict that Xingtai, Xinxiang andYangcheng stations faced more frequent droughts and two drought periods were observed in these stations, 1980-1995 and 2005-2010, 1980-1985 and 2000-2010, and 1982-1983 and 2002-2003, respectively.From 1990 to 2010 more droughts with high magnitude was recorded in the Yushe station.This result shows that relationship between climatic conditions in which changes in the temperature leads mainly to the variability in precipitation.Reference [39] reported that effects of temporal and spatial evolution of plateau monsoon caused unusual droughts and floods in China and East Asia.Figures 2(a)-2(b) show an interesting finding that index responds slowly as the time scale increased.Drought index is best in both situations: one when there are more droughts but with less frequency and the other when there are few droughts with longer duration to interpret different water resources.From Figure 2(c), it is clear that as results move to 6-month and   (1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000) and temperature was increased in all stations but maximum temperature was increased in Shijiazhuang and Yushe stations at the rate of 0.25 ∘ C/10 years and 0.19 ∘ C/10 years, respectively.The trend of precipitation and temperature may be the result of global warming and due to variation in climatic variables [40,41].In the third decade (2001-2010), temperature was also increased in all stations except Changzhi station; at this station temperature was decreased at rate of −0.01 ∘ C/10 years.On the basis of these results, it is concluded that temperature resulted in evapotranspiration and precipitation decreased in drought periods.References [42,43] found that climatic factors effect the rate of potential evapotranspiration.

Conclusions
The main objective of this study was to evaluate the presence of temporal trends of drought in Zhanghe River basin over three decades by using Standardized precipitation index (SPEI).Results exhibited extreme dry periods in all stations with different frequency.Data also reveals that there is increasing trend towards drought frequency, particularly at Anyang, Changzhi, Shijiazhuang, and Zhengzhou stations from 1980 to 1995 whereas Taiyuan, Xingtai, Xinxiang, and Yangcheng stations suffered due to drought conditions between 1996 and 2010.These results show two periods of drought in the study area from 1980 to 1995 and 1996 to 2010.Statistical results also depicted that during three decades temperature was increased particularly over the last two decades (1990-2010) over the period under study.Our results of standardized precipitation index (SPEI) at different time scales could be helpful for the planning of efficient use of water resources, hydroelectric projects, and agricultural production and designing drought information systems as well as mitigating drought impacts.

( 1 )
to research variability in precipitation at 1, 6, 12 and 24 months' time scales by applying Standardized Precipitation Evapotranspiration index (SPEI) and linear regression (2) to consider the impact of global warming and calculate change in temperature per decade (3) to investigate the drought in Zhanghe River basin during 1980-2010.

Figure 1 :
Figure 1: Locations of meteorological stations at Zhanghe River basin.
[36][37][38]om 12 to 24 months it was ranged between +2 and −2, the negative trend in SPEI values indicating a drought conditions in the study area.From Figures2(a)-2(d), it is clear that drought frequency changes with time scale changes, while as time scales increase drought becomes less frequent but persists for longer time.During 1-month time scale drought conditions appeared in all stations but in Shijiazhuang maximum drought appeared from 1980 to 1990 but drought trend was decreased as time scale was increased.In Xinxiang and Yushe stations drought conditions increased from 2000 to 2010; thus as shown in Table2, temperature was increased in these stations during these time periods.Global warming creates many problems in the hydrological cycle phenomenon like changes in the frequency and intensity of precipitation, soil moisture, flow in rivers, evapotranspiration, and ground water recharge capacity[36][37][38]. Figure2(a) also shows that drought appearance is inversely proportional to its duration; in other words, drought frequency increases but duration decreases.Another more important thing is that SPEI reciprocates quickly to increasing and decreasing trends of droughts.Figure2(b)shows the 6-month SPEI results of the study stations.Data shows that at shortest time scales like 1-month and 6-month timescales drought appurtenance was high for short time, while at the longest time scales 12month and 24-month time scales the drought was persistent for longer duration with low frequency.Thus short time scales are related to soil water content and river discharge in head water areas, while medium time scales are concerned with reservoir storages and discharge in the head water areas, and rivers and longtime scales are related to variations in ground water storage.Therefore, Standardized Precipitation Evapotranspiration index (SPEI) at different time scales are useful for monitoring drought conditions in different hydrological sub systems.At 6-month time scales all stations suffered due to droughts with different drought intensity.It is clear that temporal distribution of droughts and its frequencies are different in all regions.Dry spans are more persistent at 6month time scale as compared to 1-month scale.Figures also show that in Changzhi three droughts peaks were observed in 1985, 1990, and in 2000, while in Shijiazhuang maximum droughts appeared from 1980 to 1995.Taiyuan station comes under the stress of drought from 2000 to 2010 in last ten years.
3.3.Summary of Statistical Parameters.Statistical resultsof mean temperature at nine weather stations during the period of 1980-2010 are summarized in Table2.The mean