Spatial and Temporal Characteristics of Meteorological Drought in Shandong Province , China , from 1961 to 2008

Shandong province is located in the northern part of China and tends to be a drought-prone region. This study is dedicated to making a comprehensive and quantitative analysis of the spatial patterns of drought frequency and its climate trend coefficient, drought grades, and temporal characteristics of drought coverage area, drought duration, and drought intensity from 1961 to 2008 by using the meteorological drought composite index (CI). The results indicated that the occurrence frequency of meteorological drought in Shandong province was generally high and some part of this region such as Jiaodong peninsula had suffered drought pressure with an evident ascending trend. The drought extent and influence were very severe in 1980s and 1990s but very slight in 1964; large-area drought mainly occurred after the 1970s and the yearly and seasonally interdecadal drought duration both showed a fluctuation of 10a periodic cycle approximately. Furthermore, the slight drought mainly appeared in the northwest and southwest while the other grades of drought exhibited much significant spatial and temporal variability. Besides, drought in spring was more serious than that in winter. This study is anticipated to support the mitigation of drought hazards and to improve the management practices of environment system in Shandong province.


Introduction
According to the Intergovernmental Panel on Climate Change (IPCC) report [1], the land surface temperatures have risen globally over the past century and a half, which have shown significant regional variations also.And global warming intensifies the global hydrological cycle in terms of the averaged precipitation, evaporation, and runoff; thus it further compounds the drought occurrence around the world [2].In particular in recent decades, drought occurred frequently, and large-scale intensive drought has been observed in many continents, such as Africa [3,4], Asia [5,6], Europe [7,8], and America [9,10].Being located in East Asia, China has been suffering large-scale and long-lasting severe droughts since the early 1950s, which had negative impacts on environment, ecosystem, social-economic development, and people's living conditions [11,12].For example, the frequent occurrence of severe droughts in 1997, 1999, and 2002 in Shandong province, Hebei province, Gansu province, Inner Mongolia Autonomous Region, and so on had caused very serious economic and societal losses [13].In 1997, the very severe drought in northern China resulted in a period of 226 days with no streamflow in the Yellow River, which is the longest drying-up duration on record [2].Thus it is of great importance to carry on more thorough drought research at local and regional as well as national scales, especially to those areas with insufficient water resources.However, due to the complexity of drought dynamics and its diverse influencing factors such as temperature, relative humidity, high winds, monsoons, rains, and latitudes, the accurate monitoring and comprehensive evaluating of drought are still a big challenge.
The American Meteorological Society classifies drought into four categories, including meteorological, hydrological, agricultural, and socioeconomic drought [14].The meteorological drought, which is mainly used for monitoring and assessing the water deficit degree aroused by climate 2 Advances in Meteorology anomalies in a region within a certain period of time, is the fundamental reason for the occurrence of other droughts.Therefore, it contributes to further researches of other types of drought to explore the meteorological drought's causes and laws and to carry out effective monitoring, early warning, and historical assessment of meteorological drought hazards.As a primary method for assessing the effect of drought, the drought index is a numerical expression of drought parameters, such as intensity, severity, duration, and spatial extent.Many indices have been developed so far to quantify a drought among which some meteorological drought indices have been implemented widely at present, such as the standardized precipitation index (SPI) [15][16][17], Palmer drought severity index (PDSI) [18,19], and vegetation condition index (VCI) [20,21].Recently, the meteorological drought composite index (CI), which was proposed by China Meteorological Administration, has been used in many regions in China.Compared to meteorological indices previously reported such as SPI and PDSI as well as VCI, CI has the advantage of taking water-heat balance into consideration for drought stress characterization because it combines both precipitation and temperature information; it has the ability to identify the drought onset and termination, severity, and duration as well as drought-hit area; therefore it is suitable for real-time meteorological drought monitoring and assessment of the historical meteorological drought [22].For instance, Zhao et al. [23] made a prediction of meteorological drought risk with CI in Yunnan province, supporting the quantitative analysis of drought influencing factors such as LUCC in this region.Based on the calculation of CI, Zou et al. [24] studied the drought areas variations for ten major Chinese river basins and China as a whole from 1951 to 2008, revealing the drought areas trends and pointing out the significant drought areas in China.What is more, by employing CI, several other meteorological drought researches had been taken out in eastern [25] and northwest [26] and some other regions in China.
In this study, we made a thorough and quantitative analysis of the spatial-temporal dynamics of the meteorological drought that occurred in Shandong province from 1961 to 2008 by adopting CI which has not been employed in this region.Specifically, we analyzed the spatial patterns of drought frequency and its climate trend coefficient, drought grades, and temporal characteristics of drought coverage area, drought duration, and drought intensity in the expectation of providing reasonable reference and scientific support for alleviating and mitigating the impacts of meteorological drought., respectively.Being located in the eastern coast area of China, Bohai Sea borders its north and Yellow Sea borders its northeast and southeast (Figure 1).Being affected by the typical warm temperature monsoon climate, Shandong province is characterized by cold and less rain in spring, being hot and rainy in summer, moderate precipitation and temperature in autumn, and being cold and dry in winter, respectively.The average annual temperature is 11∼14 ∘ C, increasing form northeast to southwest.The average annual precipitation is generally between 550 and 950 mm, decreasing from southeast to northwest.The seasonal distribution of rainfall is uneven because nearly 60∼70% precipitation concentrate in summer, tending to induce or form floods.However, there are so less rainfall in autumn, winter and spring, and the droughts arouses easily in these seasons.Furthermore, with a large population and widespread farmland but limited water resources, Shandong province is facing severe water scarcity due to an increasing demand of water resource by population's livelihood, agriculture, energy, and industry sectors.Therefore, Shandong province becomes a drought-prone region, which threats regional social-economic development and living standard of local dwellers.However, as for researches that have been reported on drought hazards in Shandong province, it mainly focused on drought causes from the aspect of climatology background in terms of atmospheric circulation and monsoon [27], or drought characters by assessing statistical data such as temperature, precipitation, and disaster losses [28,29]; there is a shortage of quantitative exploration of meteorological drought from the perspective of spatialtemporal characteristics [30,31] and different severity levels with a long time series in this region.

Meteorological Data.
The meteorological data records with excellent quality through a strict quality control were selected from the Chinese Terrestrial Climate Information of Daily Datasets.There are 32 representative meteorological stations in these datasets in Shandong province; however, due to the limitation of inconsistency observation records in some stations, we finally choose 18 meteorological stations' daily climate data including measured average temperature and precipitation data from 1961 to 2008.In order to satisfy the requirements of CI calculation, the meteorological records in nearly 50 years were readily preprocessed according to the data processing document.

Methodology
3.1.Meteorological Drought Composite Index.In this study, CI was employed to analyze the spatial and temporal characteristics of meteorological drought in Shandong province at time scales of seasonal and interannual.CI is calculated based on standardized precipitation index (SPI) and relative moisture index (MI) as follows [22]: where Z 30 and Z 90 refer to the SPI of 30 days and 90 days, respectively, M 30 represents the MI of 30 days, and a, b, and c are coefficients.SPI was developed for the purpose of defining and assessing drought by McKee et al. [15], and it was computed following Zhang et al. [22]   balance between precipitation and evaporation during a period of time, and it was calculated as follows: where  refers to the total amount of precipitation in the recent 30 days (unit : mm) and PE  is the potential evapotranspiration in the recent 30 days (unit : mm).The potential evapotranspiration was calculated based on the Thornthwaite function according to Zhang et al. [22].The meteorological data including the daily average temperature and daily precipitation were employed to calculate the CI according to formula (1), and the CI was classified based on Table 1 [23].

Determination of the Meteorological Drought
Process.We adopted CI to determine the meteorological drought process.A meteorological drought process is initiated when CI in 10 consecutive days is above slight drought level (CI 10 ≤ −0.6), and this process is terminated until another 10 consecutive nondrought days appear (CI 10 > −0.6).The drought duration refers to the period between the first date and last date of the whole process.A drought event for a specific time scale (month, season, and year) appears when there is at least one occurrence of a drought process and the length of cumulative drought duration exceeds a quarter of the timescale [22].In addition, the drought intensity is referred to the sum of all CI values above slight drought level during a drought process in which smaller values indicate stronger drought intensity.
In this paper, the annual and seasonal drought duration in Shandong province were obtained by calculating the arithmetic mean of the total number of days of all drought processes for total meteorological stations within the corresponding periods.And the annual drought intensity was defined as the arithmetic mean of the drought intensity for total meteorological stations for each year.The definition of season is as follows: winter covers January, February, and December of last year, spring is from March to May, summer is from June to August, and autumn covers from September to November.

Drought Frequency and Its Climate Trend Coefficient.
The drought frequency is calculated as follows: where  is the drought frequency,  refers to the number of drought years, and  represents the total number of years.In this paper, there were 48 years from 1961 to 2008; thus  took the value of 47 because CI was calculated by rolling back.
The climate trend coefficient is used to describe the temporal characteristics of climate trend change intensity quantitatively [32]; it is calculated as follows: where   refers to the correlation coefficient between factor sequence and natural sequence in  years,   is the factor value in the  year,  represents the sample mean, and  = ( + 1)/2.

Advances in Meteorology
In this study, the climate trend coefficient of drought frequency was calculated based on the occurrence frequency of drought in five decades for each meteorological station at time scales of year and season.The positive and negative values of   reflect linear increase and decrease trends of drought frequency in  years, respectively.And the spatial distribution of drought frequency and its climate trend coefficient as well as different drought grades were obtained by employing the inverse distance weighting (IDW) interpolation method in ArcGIS 9.3 platform.

Drought Coverage Area.
The probability of drought coverage area (PDCA) is defined as the area ratio of drought events that occur for a given timescale in the study region and is calculated as follows: where  is the probability of drought coverage area,  is the number of meteorological stations in which drought events that occur for a given timescale, and  refers to the total meteorological stations.The PDCA of more than 90% is defined as the large-area drought in this research [25].Furthermore, the temporal trend analysis of drought frequency showed that the meteorological drought in Shandong province tended to be aggravating from 1962 to 2008 (Figure 3(a)).The spatial distribution of drought frequency tendency in spring (Figure 3(b)) was similar to that in autumn (Figure 3(d)) and winter (Figure 3(e)), presenting an ascent trend in eastern and northern area of Shandong province, especially in Jiaodong peninsula.In summer, the increasing trend of drought frequency was observed in finite areas, around Weifang city and Jiaodong peninsula mainly (Figure 3(c)).In general, it demonstrates that regions of significant aggravating drought during 1962 to 2008 mainly concentrated in the northern and southern part and Jiaodong peninsula, showing a certain degree of meteorological drought aggravation, and need to be focused on for mitigating the drought stresses in the future.4(c)), respectively.Spatially, the days of slight drought decreased from west to east, and the days of moderate drought showed a zonal distribution, with south being more than north, which was conversely with the pattern of severe and extreme drought.Overall, the distributions of high-value centers for every drought grade have no common features whereas low-value centers were all in the Taishan area.

Spatial-Temporal Characteristics of Meteorological Drought
In spring, the annual days of slight drought, moderate drought and severe and extreme drought were 16∼25 d (Figure 4(d)), 11∼16 d (Figure 4(e)), and 4∼9 d (Figure 4(f)), respectively.And the high value centers of slight drought were located in the western and northern regions, with drought days exceeding more than 22 d; on the contrary, the Taishan area and Weihai-Shidao in Jiaodong peninsula were not prone to slight drought.The moderate drought was mostly observed in the west to Weifang-Juxian areas.Moreover, the severe and extreme drought was the most severity in four seasons, exhibiting a high-value in the eastern part of Jiaodong peninsula and low-value in the Taishan area.
In summer, the annual days of slight drought, moderate drought, and severe and extreme drought were 10∼17 d (Figure 4(g)), 5∼11 d (Figure 4(h)), and 3∼7 d (Figure 4(i)), respectively.And the high incidence areas for each drought grades mainly located in the north, northwest, southwest, and the Jiaodong peninsula; whereas the low incidence areas were all in and around the Taishan area.
In autumn, the annual days for the three drought grades were 14∼22 d (Figure 4(j)), 9∼15 d (Figure 4(k)), and 3∼ 6 d (Figure 4(l)), respectively.The slight drought was likely to occur in the northwest, southwest, and the east area of Jiaodong peninsula.And the moderate drought days were 11∼13 d in most areas of the province, except the Yellow River delta and Rizhao city (high-value centers of moderate drought) and Taishan area (low-value center of moderate drought).In addition, the northwest and southwest were vulnerable to severe and extreme drought in autumn, with drought days more than 4 d.In winter, the annual drought days for different grades were 14∼22 d (Figure 4(m)), 9∼15 d (Figure 4(n)), and 3∼6 d (Figure 4(o)), respectively.The northwest, southwest, and eastern coastal areas were prone to slight drought in winter; however, the moderate drought days distributed evenly in the study area, ranging between 7 and 9 d generally.Besides, the severe and extreme drought in winter was the least severity in four seasons, having only 0 to 2 drought days in most parts of the province.
To sum up, each grade of drought that lasts a relatively long time has occurred in Shandong province from 1962 to 2008.And the slight drought was generally observed in the northwest and southwest both for annual and seasonal, but the other grades of drought were observed in different areas for different timescales.

Temporal Characteristics of Drought Coverage Area.
The temporal sequence of PDCA in Shandong province from 1962 to 2008 is shown in Figure 5.The annual average PDCA was 75%, and the annual large-area drought was observed in 14 a which was evenly distributed in the 1980s and 1990s.Relatively, the probabilities of drought extent are zero in 1964 and 2003, revealing that there were no drought events at time scale of year in these two years in Shandong province.For each season, the annual probability of an area affected by a drought in spring, summer, autumn, and winter was 72%, 61%, 65%, and 58%, respectively.There were 21 a of large-area drought which was mainly distributed in the 1980s, 1990s, and early 21 s century in spring, and 9 a, 12 a, and 11 a of large-area drought in summer, autumn, and winter whereas almost all of these years concentrated in the 1980s and 1990s.Meanwhile, the nondrought years for every season were relatively few over the past nearly 50 years, only 4 a in spring, 3 a in summer, and 2 a both in autumn and winter.In summary, the largearea drought at annual and seasonal time scales in Shandong province mostly occurred in the 1980s and 1990s.
Moreover, time-continuous large-area drought events have been observed in a few years in Shandong province; in specific, transseasonal large-area droughts that covers spring, summer, and autumn occurred in 1981, 1988, and 2001, respectively.Autumn-winter continuous large-area droughts 1964 1968 1972 1976 1980 1984 1988 1992 1996 2000  happened in 1990 and 2002, respectively.Serious springsummer large-area drought occurred in 1992 and summerautumn large-area drought occurred in 1997, respectively.It can be seen that time-continuous large-area droughts mainly took place after the 1970s in Shandong province, revealing that time-continuous droughts were steadily increasing with time during the study period.However, it is also found that the linear changing tendencies were not obvious in the annual and seasonal PDCA in Shandong province by carrying out the  significance test with the confidence level of  = 0.01 and  = 0.05.

Temporal Characteristics of Drought Duration.
As is shown in Figure 6, the annual and seasonal drought durations in Shandong province were both fluctuated significantly from 1962 to 2008.The annual drought duration was 153 d, with the longest drought duration of 291 d in 1999 and the shortest drought duration of 8 d in 1964.There were 15 years that the annual drought duration of which was more than 183 days (half of the total days of a year), and furthermore, these years were mainly in the 1980s and 1990s.On seasonal, the annual drought durations in summer and winter were both 34 d and in spring and autumn were 45 d and 39 d, respectively.Over the past nearly 50 years, there were 24 a, 12 a, 16  It reveals that the periodic variation of the annual and seasonal interdecadal drought duration was roughly similar, showing a cycle of 10 a approximately (Table 2).But the linear changing tendencies of drought time duration in Shandong province from 1962 to 2008 were insignificant through executing the  significance test with the confidence level of  = 0.01 and  = 0.05.

Temporal Characteristics of Drought Intensity.
Over the past nearly 50 years, the annual drought intensity, which is referred to the arithmetic mean of the sum of all the days of CI value that are all above slight drought level for total meteorological stations for each year in Shandong province, showed a fluctuation tendency (Figure 7).The minimum CI value which indicates the strongest drought intensity was −369.93 and appeared in 1988, followed by −333.56 in 1999, and the maximum CI value conversely implying the lightest drought intensity was −15.22 in 1964.There are 19 a that annual drought intensity was above the mean intensity for total years from 1962 to 2008, and those years were mainly in the 1980s, 1990s, and early 21 s century.However, the linear tendency of annual drought intensity has not passed the  significance test with the confidence level of  = 0.01 and  = 0.05, revealing that the variation tendency of annual drought intensity was much complicated in Shandong province.
The linear tendencies of PDCA and drought duration for annual and seasonal and annual drought intensity were all insignificant; however, it is still obvious that the temporal drought parameters were salient in the 1980s and 1990s than in any other decades, which gives an evidence of the highlighting drought extent and influence in these two decades in Shandong province.In addition, the three temporal drought parameters were all the palest in 1964, demonstrating that the drought was lightest in this year in Shandong province.Moreover, taking into account the annual drought duration and intensity, it is found that the drought duration in 1999 1964 1968 1972 1976 1980 1984 1988 1992 1996 2000    In general, the three temporal parameters of drought, namely, PDCA, drought duration, and drought intensity, effectively depict the drought from different aspects and all of them are of great importance to the evaluation of the meteorological drought.

Conclusion and Discussion
Owing to various factors such as global warming, abnormity of East Asia monsoon wind, and local precipitation decreasing, the drought hazard in Shandong province is becoming more and more serious.Therefore, on the basis of drought analysis by adopting and calculating the CI for 18 meteorological stations, we comprehensively evaluated the spatial-temporal characteristics of meteorological drought in Shandong province from 1962 to 2008.Main conclusions are as follows.
The frequency of meteorological drought was generally high in Shandong province, and high-frequency regions are mainly concentrated in the southwest and north areas whereas low-frequency regions are generally centered in the central and southeast.The spatial pattern of drought frequency is closely related to the distribution of precipitation which decreases from southeast to northwest and temperature which increases from northeast to southwest.
The seasonal drought frequency followed such an order of spring > autumn > summer > winter.The precipitation is rare in spring but the evaporation is high due to the land surface warming rapidly, which aggravates the spring drought.In contrast, abundant rain falls in summer and low evaporation in winter owing to low temperature result in relatively low drought frequency in these two seasons.
The slight drought, moderate drought, and severe and extreme drought were all observed in Shandong province from 1962 to 2008.And the slight drought generally occurred in the northwest and southwest, but drought of other grades was observed in different areas at different timescales.Additionally, various drought grades were longest in spring but shortest in winter which was consistent with the statistics of the calculated drought parameters and, therefore, further revealed that the spring drought was much more serious than that in the other three seasons.
With respect to the temporal trend of drought frequency, it demonstrated that regions including north, south, and Jiaodong peninsula were all suffering an ascend drought over the past nearly 50 years, but the upward trend in other regions was much evident.Besides, temporal trends of multiple drought parameters including PDCA, drought duration, and drought intensity all indicated that the drought extent and influence were very severe in the 1980s and 1990s but much mild in 1964; the time-continuous largearea droughts mainly occurred after the 1970s, and the interdecadal drought duration at annual and seasonal scales both showed a fluctuation tendency of 10 a periodic cycle approximately.
In brief, this study proves that Shandong province is a sensitive area to global climate change from the perspective of meteorological drought occurrence.What is more, the much complicated spatial-temporal characteristics of meteorological drought in Shandong province detected by this study do confirm that the regional response to climate change exhibits very significant spatial-temporal differences.The findings summed up above will in turn give us essential knowledge and guidelines on adaptation and mitigation options to address climate change.In specific, in Shandong province, critical works to be carried out are as follows: (1) to conduct more wide and more extensive monitoring and observations in terms of climate change and drought hazards, and to carry out more intensive scientific researches that focus on the dynamics of meteorological drought; (2) to strengthen drought risk management practices in the public sector, especially to choose and enact defense-oriented strategies, and to develop drought emergency plans and measures at the same time; (3) to maintain and make use of the ecosystem functions and services in terms of drought risk defense and water conversation by protecting and restoring ecosystems in ecocritical areas and drought prone areas; (4) to investigate, choose, develop, and promote effective measures of drought defense and prevention among agricultural, industrial, and energy sectors and people's lives, such as dry farming and water-saving agriculture, green-water utilization technologies, industrial water-saving technologies, and domestic water conservation technologies; (5) to promote public education and to motivate public participation with respect to drought defense and prevention, and so on.

Figure 1 :
Figure 1: Location map of Shandong province and the meteorological stations.

4. 2 .
Spatial Patterns of Different Drought Grades.The spatial patterns of different drought grades from 1962 to 2008 in Shandong province are shown in Figure 4.In the respect of annual drought, over the past nearly 50 years, the number of slight drought, moderate drought, and severe and extreme drought days that occurred in Shandong province were between 57∼81 d (Figure 4(a)), 32∼49 d (Figure 4(b)), and 11∼23 d (Figure

Figure 5 :
Figure 5: Annual, spring, summer, autumn, and winter probability of drought coverage area in Shandong province.

Figure 7 :
Figure 7: The annual drought intensity in Shandong province.
2.1.Study Area.Shandong province is situated in the downstream area of the Yellow River; it lies between E 114 ∘ 36  ∼E 122 ∘ 43  and N 34 ∘ 25  ∼N 38 ∘ 23 in this paper.MI indicates the

Table 1 :
Classification standard of the meteorological drought composite index (CI).
S l i g h t d r o u g h t −1.2 < CI ≤ −0.6 3 M o d e r a t e d r o u g h t −1.8 < CI ≤ −1.2 4 Heavy drought −2.4 < CI ≤ −1.8 5 Extreme drought CI ≤ −2.4 a, and 12 a that annual drought duration was more than 46 d (half of the total days of one season) in spring, summer, autumn, and winter, respectively.The longest annual drought duration for spring, summer, autumn, and winter was 84 d in 1988, 76 d in 1992, 78 d in 2002, and 89 d in 1999, respectively; correspondingly, the shortest annual drought duration in spring and summer was 0 d in 1990, while in autumn and winter was 0 d in 1964 and 2003 and 1 d in 1964, respectively.