Abstract

Drought is a natural disaster caused by long-term water deficit. Because the growth of crops and vegetation is closely related to soil moisture environment, it is of great significance to study the soil drought and vegetation response. In this paper, the soil moisture availability index (SMAI) was developed for quantifying soil drought conditions. The effectiveness and the ability of SMAI to recognize drought events were analyzed, while the vegetation condition index (VCI) was used to characterize the vegetation status. Temporal and spatial variations of soil drought and vegetation condition as well as the impacts of drought on vegetation in North China during 2001–2015 were comprehensively examined. We firstly concluded that SMAI related well with standardized precipitation evapotranspiration index (SPEI), and drought events can be detected by SMAI. Next, the mean value of SMAI in North China showed a decreasing trend in recent 15 years. Finally, the SMAI positively correlated with VCI in most areas of North China, and the response of four types of vegetation to SMAI differed over time. The results of SMAI on vegetation would assist drought research and application in North China.

1. Introduction

Drought is a natural phenomenon caused by imbalance in water budget, thus leading to water shortage, a sharp drop in groundwater reserves, crop reduction, and other issues. Unlike other forms of natural disasters, droughts often develop slowly and its effect varies from region to region. In the background of global warming, evaporation is getting more intense, so drought is usually becoming more and more serious. The area and severity of drought-affected regions have increased in recent years in China [1, 2], especially in the northern China where water resources are relatively scarce. Drought is a long and sustained process; the impact of the initial drought is not obvious, but with the increase of the drought duration, it will cause a series of problems. Vegetation is an important component of land surface coverage, and the response of vegetation to climate change is an important part of global or regional change research issue. Vegetation is the indicator of an ecological environment, including the feedback of drought. Therefore, the study of regional drought and vegetation response is helpful to understand the evolution and regularity of drought and the impact on terrestrial ecological environment, which can provide a reference for the monitoring and evaluation of the regional ecological environment.

Su et al. have summarized drought monitoring and assessment using remote sensing data. Many indices have been proposed for drought assessment, and each of these indices addresses a specific aspect of drought, including meteorological, soil moisture, hydrological, and socioeconomic drought [3]. The selection of drought indicators is crucial to drought evaluation, while the lack of precipitation and the high temperature are often the beginning of meteorological drought. Therefore, the previous studies of drought are mostly based on the precipitation from the meteorological perspective, mainly by the Z-index of available rainfall [4], standardized precipitation index (SPI) [5], palmer drought severity index (PDSI) [6], and standardized precipitation evapotranspiration index (SPEI) [7]. The correlation between vegetation parameters and meteorological parameters (such as precipitation and temperature) has been studied [810], as well as the response of vegetation to the meteorological drought index [11, 12]. As the reflection of precipitation and evapotranspiration, soil moisture deficit is a comprehensive feedback of meteorological drought in a certain period of time [13, 14]. Furthermore, the growth status of vegetation is closely related to the soil water; therefore, it is also necessary to analyze the response of vegetation to soil drought conditions. Considering the seasonal and cumulative effects of drought, the soil moisture deficit index (SMDI) [15] is chosen as the evaluation index of drought.

With the characteristics of the complex and diverse land surface belonging to semihumid areas, North China is sensitive to global climate changes [14]. The seasonal variation of vegetation is significant and especially susceptible to drought conditions because of the low annual precipitation in North China. Considering the impact of drought on vegetation varied in different regions and vegetation types [16], as well as water demand of vegetation differed over seasons [17], it is necessary to study the response of vegetation conditions to soil drought under different conditions from a more detailed spatial and temporal scale [18].

2. Study Area

The study area is located in North China (31°23′–42°37′N, 110°15′–122°43′E). It includes the province of Hebei, Henan, Shanxi, and Shandong and the municipalities of Beijing and Tianjin (Figure 1), with an area of 6.929 × 105 km2. The study area is in the temperate semihumid continental monsoon climate zone. The annual precipitation is between 400 and 800 mm (decreasing from southeast to northwest), mainly affected by the Pacific monsoon intensity. Due to the precipitation, evaporation and temperature varied greatly at interannual and intra-annual timescales, and drought and flood hazards occur mainly in spring and summer, respectively. North China is an important agricultural and livestock breeding base in China. The DEM in Figure 1 shows the terrain of the study area. The northern and western regions of the study area are mountainous, and the main vegetation types are woodland and grassland. While the central and eastern regions are a large area of plain, and the main crops are wheat, corn, peanuts, cotton, and so on.

3. Data and Methodology

3.1. Drought Monitoring Index

Soil condition is closely related to agriculture and vegetation status, and soil moisture is often regarded as one of the indices of agricultural drought. SMDI [14] was selected as the drought monitoring index in this study. Because SMDI increases with increasing soil moisture, we renamed it as “soil moisture availability index” (SMAI). The SMAI is calculated as follows:where is the actual monthly soil moisture at month and , , and are the long-term maximum, minimum, and median soil moisture at month , respectively. is defined as the monthly soil moisture availability with the range of −100 to 100. SMAI considers both the soil moisture condition at the current month and that of the previous month. The initial value of SMAI was calculated only considering the SA with the range of −2 to 2. With the iterative computations, the range of SMAI values is −4 to 4 [14].

Considering that the soil moisture observation data are limited and the geographical distribution is uneven [19], while with the advance of remote sensing, soil moisture data can be extracted in a short time for large areas. With the continuous improvement and development of global soil water monitoring, the global land surface data products by the Global Land Data Assimilation System (GLDAS) have been widely used in related research. The applicability of these data in North China has also been documented [20]. Therefore, we used GLDAS data as input data for soil drought index calculations in this study. The monthly soil moisture data (unit: kg/m2) of GLDAS_NOAH025 (http://disc.sci.gsfc.nasa.gov/hydrology/data-holdings) were used in this study, and the data version is 2.1, which consists of 4 layers of soil water (0-0.1 m, 0.1–0.4 m, 0.4–1.0 m, and 1.0-2.0 m underground, resp.) with 0.25° × 0.25° spatial resolution, about 25 km × 25 km.

As an indicator of meteorological drought monitoring, SPEI has multiple timescale features while considering the precipitation and potential evapotranspiration simultaneously. Studies have shown that the SPEI is in good agreement with the actual drought condition both in China and in North China [21, 22]. Therefore, TRMM and MODIS global evapotranspiration product data are used to calculate the SPEI in this study, and the relationship between SPEI and SMAI will be analyzed for the applicability validation of SMAI in our study area.

3.2. Vegetation Monitoring Index

The normalized difference vegetation index (NDVI) is the most widely used vegetation index. It has a strong correlation with vegetation coverage, vegetation growth status, leaf area, biomass, chlorophyll, and other physiological parameters [23]. The vegetation condition index (VCI) obtained by normalizing NDVI for many years can be used as an indicator of the degree of environmental stress. It is widely used in agricultural drought monitoring and other fields. The formula is as follows.where , , and are monthly NDVI, multiyear maximum NDVI, and multiyear minimum NDVI, respectively, for each grid cell. VCI changes from 0 to 100 corresponding to changes in vegetation condition from extremely unfavorable to optimal. Because the effect of drought on vegetation growth was mainly in the growing season and the growing season of vegetation and crops in the North China is usually in April to October, the monthly 1 km MOD13A3 NDVI product (https://lpdaac.usgs.gov/dataset_discovery/modis/modis_products_table/mod13a3) of growing season was used in this study.

3.3. Trend Analysis

In this paper, the temporal evolutions of drought and vegetation index were estimated using the slope of the time series of data. Taking the monthly variation of SMAI as an example, the formula is as follows:where is the total number of months and is the SMAI value in months.

3.4. Correlation Analysis

In order to analyze the response of vegetation growth to drought, the Pearson correlation analysis was done between SMAI and VCI. Because the variables are normal continuous variables, their correlation degree can be measured by the Pearson correlation coefficient.

3.5. Drought Event Identification and Its Relationship with Meteorological Conditions

It is difficult to capture both the duration and the location of the drought simultaneously based on the temporal and spatial characteristics, respectively. Therefore, the combination of the above two characteristics can identify the duration, influence area, and severity of drought. In the past, the drought index in a three-month scale, SPI-3, RDI-3 (reconnaissance drought index), and SPEI-3 were selected for seasonal drought study.

With the SMAI calculated for each 0.25° × 0.25° grid from 2001 to 2015, drought events are extracted using the proposed space-time continuum identification method. According to the study area and the data spatial resolution, we will identify drought events that are over three months and cover an area of more than 10000 km2. A threshold value of SMAI-1 should be set to identify the drought condition.

For each month during 2001–2015, each successive arid region is marked as a uniquely coded drought patch. When the area of overlapping regions of two consecutive months of dry patches is greater than 10000 km2, it is marked as the same drought event. So, the area where the drought occurs, its duration, and the severity of the drought can be tracked [24].

Drought are closely related with precipitation and temperature, so we also collected the average precipitation and temperature in each month from 2001 to 2015 (data from China Meteorological Data Network: www.data.cma.cn) and calculated the precipitation and temperature anomalies during the drought period.

4. Results and Discussion

4.1. Relationship between SPEI and SMAI

Different drought indices are concerned with different perspectives of drought. The SPEI focuses on the meteorological drought, and the SMAI focuses on the soil drought. The various links of the drought are interrelated. In order to verify the consistency and accuracy of the drought index selected for this study, the correlation coefficient is calculated, including four-layer soil drought conditions reflected by SMAI and five timescale meteorological drought characteristics reflected by SPEI. The monthly SMAI was calculated by the soil moisture data of the GLDAS in the four soil layers (SMAI-1: 0.0 m-0.1 m, SMAI-2: 0.1 m–0.4 m, SMAI-3: 0.4 m–1.0 m, and SMAI-4: 1.0 m-2.0 m) in 2001–2015. The SPEI has an advantage of the multiple timescales. SPEI-1, SPEI-3, SPEI-6, SPEI-9, and SPEI-12 are calculated, respectively, on 1-, 3-, 6-, 9-, and 12-month scales. For example, SPEI-1 indicates the degree of meteorological drought in the current month, and SPEI-3 indicates the comprehensive meteorological drought situation including the current month and the previous two months.

The results (in Table 1) showed that the correlation coefficient of SMAI-1 and SPEI-3 was the highest and reached a significance level of α = 0.01. Overall, SMAI responded with SPEI, with the increase in timescale and the decrease in soil depth, and the correlation coefficient between SMAI and SPEI increased.

4.2. Temporal and Spatial Variations of Drought in North China

Time series of the monthly mean SMAI of the four soil layers in the 180 months are shown in Figure 2 so as to reflect the change trend of soil drought in North China.

The results showed that the variation trends of SMAI in the four layers were almost the same in the past 15 years in the North China region. However, SMAI-2, SMAI-3, and SMAI-4 have a time lag of 1 to 7 months compared with SMAI-1. The overall change trends of SMAI decreased; however, the change trends of SMAI in the top three soil layers were not significant, and only the trend of SMAI-4 was significant and it reached a significance level of α = 0.01. The results indicated that the degree of soil drought increased.

The values and variation among SMAI in the first three layers are similar, while SMAI-4 has the smallest fluctuation. Taking SMAI-1 as an example, the soil moisture condition was the best during 2003 to 2005, and there was a significant long-term drought period before and after this time. From the end of 2003 to the end of 2006, the average SMAI value dropped most significantly from 3.30 to −2.29. After that, SMAI fluctuated up and down around 0, declining sharply in 2007, 2009, and 2011, and soil drought indices reached low values. The drought in 2014 was also severe; on average, SMAI was negative for ten months in 2014.

Considering that the values of SMAI in the top three layers are similar, the property of surface soil is closely related to vegetation growth. The annual SMAI-1 was calculated by the mean of monthly SMAI-1. The spatial distributions of annual SMAI-1 during 2001–2015 are shown in Figure 3. The results indicated that the annual SMAI-1 was relatively high in 2003, 2004, and 2005 in the whole study area, and there is no obvious drought area. On the contrary, the annual SMAI-1 was relatively low in 2001, 2002, 2006, 2007, 2011, 2014, and 2015. In order to facilitate the comparison of similar studies, the zero baseline for SMAI was adopted to analyze soil drought/moistness. The area where the SMAI is less than 0 is considered to be drought, and the drought area accounted for 97%, 79%, 90%, 68%, 94%, 83%, and 81%, respectively, in the above year.

The monthly slope of SMAI-1 was calculated by formula (3) for each pixel, and the result is shown in Figure 4. The negative value indicates that the area tends to be drought, and with the decrease of the slope value, the degree of drought increases. The results in Figure 4 showed that the trend of drying was the most obvious in the southern part of the study area, and the trend of wetting was obvious in the north. The area where the slope of SMAI is less than 0 is considered to be drying areas, which account for 72% of the study area, and 50% of the drying areas reached statistical significance. The degree of drying is more serious in the southern part of the study area, including Henan Province, Shandong Province, and part of the mountainous area in Shanxi Province. The overall situation of soil wetness in the Beijing-Tianjin-Hebei region has improved, especially in the northeastern part of Hebei Province. Overall, SMAI trend analysis reflected that Henan Province tended to be drying significantly during 2001–2015.

4.3. Persistent Drought Events in North China

We tested the threshold of SMAI-1 values with −1, −1.5, and −2.0 to detect drought events, and the drought events detected by a value of −1.5 matched well with the records from China flood and drought disaster bulletin (2006–2015). Therefore, the SMAI threshold was set as −1.5 in this study. A total of 19 persistent drought events were detected during 2001–2015 in North China (Table 2), of which the duration is the number of months that the drought event lasts, the area is the maximum area affected by drought during the drought period, and the integrated drought severity index (IDSI) is the sum of the monthly SMAI values in the drought area during the drought period.

The time period of 19 persistent drought events in North China during 2001–2015 was drawn on the monthly time series of average SMAI in Figure 5. The persistent drought events also occurred when the whole average area was in drought conditions in North China. Because the IDSI can reflect the duration, area, and soil moisture of drought events, we sort the intensity of the 19 drought events by IDSI. The most comprehensive drought event from November 2010 to September 2011 was as long as 11 months, with an IDSI of −14835, extending 97% of the study area, followed by the drought from April 2001 to January 2002, which lasted up to 10 months, and almost the entire study area was affected by the above drought events. However, in other time periods, although the overall drought situation in North China was not obvious, drought also occurred in some areas. For example, in Figure 5, the overall drought index in North China in 2004 was above the average, but there were still persistent droughts in southern Shanxi.

In order to verify the accuracy of the drought event identified in the study, the drought duration and distribution reflected by the SMAI were compared with the China flood and drought disaster bulletin records (2006–2015). The results showed in Table 2 that the drought conditions reflected by SMAI are basically consistent with the data records.

Taking one of the most serious drought events as an example, the status of precipitation and temperature during the drought period is analyzed. The drought event began in November 2010 and ended in September 2011, affecting most of the study area. The monthly precipitation anomaly, the temperature anomaly, and the SMAI during the drought event are shown in Figure 6. It can be seen from November 2010 to June 2011 that the precipitation in the study area has been at a low level. However, no obvious high temperature appeared during this period, and the drought was mainly caused by the less rainfall. The calculated SMAI of the same period also shows a similar rule to the change in precipitation.

4.4. Temporal and Spatial Variations of Vegetation in North China

Due to the phenology of vegetation, the vegetation condition varied over time. VCI normalizes NDVI and allows for a comparison of different ecosystems. The results showed that the VCI value was the lowest in April, followed by October, May, and June in North China. The monthly mean values of VCI in these above months are below 50, indicating that the vegetation status in these months was poor. In order to explore the spatial distribution of the vegetation index in North China, we calculated the mean value of growing season VCI in North China (in Figure 7). In 2001 and 2002, the VCI value was much lower than that in other years. Then, the vegetation status has a stable improvement in whole regions until 2011. Since 2012, the vegetation condition had improved more significantly, especially in the northern and western regions of the study area.

The trend analysis method was also used to study the spatial-temporal variation characteristics of vegetation. The spatial distribution of the annual VCI linear trend is shown in Figure 8. The slope value above or below 0 indicates that the vegetation condition is improved or degenerated.

The trend of VCI data indicated that 88% of the areas were improved in the study area, mainly distributed in the high-altitude area. In contrast, the area of vegetation degradation was distributed in the southeast part with low altitude. It is concentrated in the southern part of Henan Province, the central and eastern coastal areas of Hebei Province, and the small part of coastal areas in Shandong. For the temporal change of VCI, 70.41%, 65.00%, and 52.28% areas reached the positive significance level of α = 0.1, α = 0.05, and α = 0.01, respectively, while 3.78%, 2.90%, and 1.57% areas reached the negative significance level of α = 0.1, α = 0.05, and α = 0.01, respectively.

4.5. Vegetation Response to Drought in North China
4.5.1. Response of Annual Vegetation Condition to Drought

To analyze the vegetation response to the drought, we calculated the correlation coefficient between SMAI and the corresponding vegetation index at the pixel scale during 2001–2015. The annual values of SMAI-1 were calculated from the average of monthly data of SMAI-1. The SMAI data were resampled to 1 km by 1 km spatial resolution. MODIS Land Cover data in MCD12Q1 in 2013 were used to analyze whether the differences of the drought effects were related to vegetation types. According to the MCD12Q1’s IGBP global land cover classification scheme and primary vegetation types in the study area, the land cover data were reclassified into eight categories, including woodland, grassland, shrub, farmland, city, bareland, waterbody, and wetland, and the first four were related with vegetation.

The correlation coefficient between SMAI and VCI is shown in Figure 9(a). There was 76% of the total area with positive value, with a mean value of 0.30, and 19%, 13%, and 5% areas reached the positive significance level of α = 0.1, α = 0.05, and α = 0.01, respectively. The high positive correlation between SMAI and VCI was mainly distributed in the northeast of Hebei Province, the junction of Hebei Province and Shanxi Province, and the central and eastern parts of Shandong Province. The mean values of the positive correlation coefficient for shrub, woodland, grassland, and farmland were 0.42, 0.39, 0.32, and 0.28, respectively, and none reached the significance level of α = 0.1.

4.5.2. Response of Monthly Vegetation Condition to Drought

To explore the response of vegetation to soil drought at different phenological periods, we analyzed the relationship between vegetation condition and drought at the monthly scale in the vegetation growing season (from April to October). Through normalization, the VCI eliminates the difference to make vegetation changes comparable on the monthly scale. The response of monthly vegetation condition to drought was analyzed between SMAI and VCI. The results are shown in Figures 10 and 11.

Generally, the results showed that the positive coefficient between SMAI-1 and VCI accounted for almost 91% of the total area, and about 37% of them reached the significance level of α = 0.01, which indicated that surface soil drought affected vegetation significantly in the growing season. Vegetation water demand varies with the vegetation type and seasonal variation. From April to October, the mean values of SMAI-VCI correlation coefficients in the whole study area were 0.161, 0.218, 0.222, 0.223, 0.259, 0.158, and 0.171, respectively, and the positive correlation proportions were 76.9%, 78.5%, 80.8%, 80.2%, 81.9%, 71.4%, and 73.4%, respectively, with mean values of 0.259, 0.322, 0.310, 0.317, 0.354, 0.293, and 0.293.

From the temporal evolution analysis, the correlation coefficient between SMAI and VCI increased from April to August gradually and then decreased significantly from August to September overall, which indicates the vegetation response to soil drought most significantly in August. In terms of spatial distribution, the correlation coefficient was high in the north part of the study area in the whole growing season, while as for the middle part of the study area, the correlation coefficient was high during May to August.

The monthly mean values of SMAI-VCI correlation coefficients of the four vegetation types in which the significance level reached α = 0.1 are shown in Figure 12. There are two peaks of the correlation coefficient between SMAI and VCI in the growing season for the four types of vegetations, one occurred in May and the other in August or September. The maximum correlation coefficient of woodland and shrub occurred in September, while for grassland and farmland, it occurred in August. The maximum correlation coefficient between SMAI and VCI for the four types of vegetation decreased by shrub, woodland, grassland, and farmland, which is reasonable, besides that the soil drought, soil fertility, irrigation, and other human interventions have a great influence on crop growth.

5. Conclusions

This paper investigated the variation of soil drought and vegetation indices and their correlations during 2001–2015 in North China. The following important conclusions can be drawn from this study:(1)The SMAI related well with SPEI which is widely used in meteorological drought assessment, and drought events detected by SMAI matched well with the “China flood and drought disaster bulletin (2006–2015)”.(2)The monthly average drought degree by the time-series SMAI has an increasing tendency generally in the study area, and the trend of drought is the most significant in the deepest soil. The drought degree varied spatially. The most severe drought area was located in the southern part of the study area, mainly in Henan Province. While in the north part of the study area, there is a trend of wetting.(3)The annual trend analysis of VCI showed that the vegetation condition was generally improved in North China.(4)The correlation coefficient between the vegetation index and drought index was positive both on the annual scale and monthly scale in most of the study area. The response of VCI to SMAI in the grassland, shrub, woodland, and farmland was significantly different. The correlation coefficient of farmland was much lower than that of the other vegetation types because human factors or natural factors including land management, soil fertility, crop harvesting, and other effects will also affect the crop growth.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This paper was jointly supported by the National Key R&D Program of China (2017YFC0406004). The authors also would like to thank Professor Richard P. Greene and Professor Bob Su who provided help.