Abstract

Droughts are projected to increase in severity and frequency on both regional and global scales. Despite the increasing occurrence and intensity of the 2009/2010 drought in southwestern China, the impacts of drought on vegetation in this region remain unclear. We examined the impacts of the 2009/2010 drought in southwestern China on vegetation by calculating the standardized anomalies of Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), Normalized Difference Water Index (NDWI), and Land Surface Temperature (LST). The standardized anomalies of NDVI, EVI, and NDWI exhibited positively skewed frequency distributions, while the standardized anomalies of LST exhibited a negatively skewed frequency distribution. These results implied that the NDVI, EVI, and NDWI declined, while LST increased in the 2009/2010 drought-stricken vegetated areas during the drought period. The responses of vegetation to the 2009/2010 drought differed substantially among biomes. Savannas, croplands, and mixed forests were more vulnerable to the 2009/2010 drought than deciduous forest and grasslands, while evergreen forest was resistant to the 2009/2010 drought in southwestern China. We concluded that the 2009/2010 drought had negative impacts on vegetation in southwestern China. The resulting assessment on the impacts of drought assists in evaluating and mitigating its adverse effects in southwestern China.

1. Introduction

Drought has increased in severity and frequency over the past decades [1]. Particularly, the 2009/2010 drought in southwestern China was reported as a “once-in-a-century” drought [2]. The occurrence and intensity of the 2009/2010 drought in southwestern China have been highlighted by [3], but the impacts of drought on vegetation in this region remain unclear. Hence, understanding the impacts of this drought is important for assisting in evaluation and prediction of droughts for reducing economic losses in southwestern China. Previous studies on the impacts of drought on vegetation have mainly focused on Amazon rainforests [46], North America [7, 8], Mongolia Plateau [9], and northern China [1012]. Little research has examined the responses of vegetation to drought in southwestern China.

Satellite-derived vegetation indices (e.g., NDVI and EVI) have been used widely to surrogate the vegetation canopy greenness [4, 13, 14]. Hence, it is vital to investigate the impacts of the 2009/2010 drought on vegetation in subtropical southwestern China using MODIS-derived NDVI and EVI in order to evaluate the responses of vegetation greenness to drought stress. NDWI reflects changes of spongy mesophyll and water content in vegetation canopies [15]. MODIS-derived NDWI has been widely used to detect the moisture conditions of vegetation canopies [16, 17]. MODIS-derived LST is a good indicator of the thermal conditions of the environment [18]. Therefore, the standardized anomalies of MODIS-derived NDWI and LST can reflect the changes in vegetation water content and thermal conditions of the environment, respectively.

In this study, the MODIS-derived NDVI, EVI, NDWI, and LST were proposed as variables for calculating the standardized anomalies to acquire specific information about vegetation greenness, vegetation water content, and land surface temperature, respectively. These variables are significant for examining the responses of vegetation status and thermal conditions of the environment to the 2009/2010 drought quantified by the global terrestrial drought severity index (DSI) [19]. The objective of this study is to assess the responses of vegetation to drought stress in areas of southwestern China that suffered from the 2009/2010 drought.

2. Study Area

The study area was located in southwestern China (21°N–35°N, 97°E–113°E), which consists of five provinces: Sichuan, Chongqing, Guizhou, Guangzhou, and Yunnan (Figure 1). According to the 2008 MODIS land cover data (MCD12Q1), these regions are mainly covered by mixed forests, deciduous forest, evergreen forest, savannas, grasslands, and croplands. The study area is in the subtropical monsoon climate zone [20]. The monthly mean air temperature and annual mean precipitation are approximately 10°C and 900 mm, respectively [21].

3. Data and Processing

A range of MODIS datasets was used in this study (Table 1). They consist of MODIS vegetation indices (NDVI and EVI) products (MOD13A3) (Collection 5), MODIS surface reflectance product (MOD09A1) (Collection 5), MODIS land surface temperature product (MOD11A2) (Collection 5), and MODIS land cover product (MCD12Q1) (Collection 5) acquired from the US National Aeronautics and Space Administration (NASA) Land Processes Distributed Active Archive Center (LP DAAC) (https://lpdaac.usgs.gov/).

The MOD13A3 (monthly, 1 km), MOD09A1 (8-day, 0.5 km), and MOD11A2 (8-day, 1 km) products for the period 2000–2011 were used to calculate the standardized anomalies of vegetation greenness (NDVI and EVI), vegetation water content (NDWI), and land surface temperature (LST) to examine the effects of the 2009/2010 drought on vegetation in southwestern China, respectively. These MODIS products were aggregated into a monthly interval by calculating the average value over the corresponding period and then resampled at a 0.05° spatial resolution to maintain the same temporal and spatial resolutions for different types of datasets.

The 2008 MODIS land cover product (MCD12Q1) was used to locate different vegetation types in the study area: deciduous forest, evergreen forest, mixed forests, croplands, grasslands, and savannas. Other land cover classes, such as water, urban and built-up, barren or sparsely vegetated, were aggregated as nonvegetation.

The global terrestrial drought severity index (DSI) products were generated by [19] to monitor the global terrestrial drought for the period 2000–2011 by integrating the MODIS-derived NDVI and Evapotranspiration/Potential Evapotranspiration (ET/PET) data. The reasons of choosing DSI as a standard to evaluate the results from different RS-based indices are the following: () DSI has advantages in high spatial resolution with 0.05 degrees comparing to meteorological-based drought index such as Palmer Drought Severity Index (PDSI) with 0.25 degrees [22]; () DSI could be more sensitive to drought stress comparing to other remotely sensed drought indexes such as Vegetation Condition Index (VCI), which is driven only by NDVI [23, 24], as DSI is designed to detect the anomalies of photosynthetic activity of vegetation conditions, evapotranspiration and potential evapotranspiration of soil and vegetation activities under the stress of drought [19]; () this study is a further work following the previous one by Zhang and Yamaguchi [3]. So we used the same drought severity index (DSI) to evaluate the impacts of a serious drought event on vegetation over southwestern China.

DSI data was acquired from the Numerical Terradynamical Simulation Group (NTSG) at the University of Montana (http://files.ntsg.umt.edu/data/NTSG_Products/DSI/). In this study, the DSI products at 8-day intervals with 0.05° spatial resolution from 2000 to 2011 were composited to a monthly interval to characterize the extent, duration, and severity of the 2009/2010 drought in southwestern China. The DSI values less than or equal to −0.3 were categorized as drought [19]. Monthly DSI values for the drought period from November 2009 to March 2010 were averaged to represent the dry period mean of the 2009/2010 drought in southwestern China.

Our study area (21°N–35°N, 97°E–113°E) was covered by five tiles of MODIS images. The MODIS tile numbers are h26v05, h26v06, h28v06, h27v05, and h27v06. Five tiles of original MODIS images (MOD13A3, MOD09A1, and MOD11A2) were mosaicked and cut using MODIS Reprojection Tool (MRT) developed by NASA.

4. Methodology

4.1. Calculating the Standardized Anomaly (SA)

The standardized anomalies (SAs) of NDVI, EVI, NDWI, and LST are calculated aswhere SA is a given variable (e.g., NDVI, EVI, NDWI, and LST) calculated from its value for the 2009/2010 drought period () relative to long-term mean () and standard deviation () over a reference period. The 2009/2010 drought period for southwestern China was characterized by [3] as “November 2009 to March 2010.” Thus, in this study, the monthly MODIS composites (NDVI, EVI, NDWI, and LST) were averaged for November, December, January, February, and March to represent the dry period mean for each year over the reference period 2000–2011, but excluding the 2009/2010 drought period.

In this study, only vegetated areas that suffered from the 2009/2010 drought were selected for calculating the standardized anomalies. Regions labeled “no drought” (DSI equal to or higher than −0.29) and “no data” (see Figure 2) were excluded.

4.2. Evaluating the Impacts of 2009/2010 Drought on Vegetation in Southwestern China

The impacts of the 2009/2010 drought on vegetation in southwestern China were evaluated and analyzed using monthly averaged standardized anomalies of NDVI, EVI, NDWI, and LST. They were first calculated pixel by pixel for vegetated areas affected by the 2009/2010 drought quantified by DSI during the drought; then frequency histograms of the standardized anomalies of NDVI, EVI, NDWI, and LST were generated; finally, nonlinear curves of the frequency histograms of the standardized anomalies for NDVI, EVI, NDWI, and LST were fitted by a Gaussian function to show the skewed frequency distributions.

The NDVI, EVI, NDWI, and LST were used as RS-based indices to compare with DSI in this study. Pearson’s coefficient () and value were used to evaluate the relationships between DSI and standardized anomalies of NDVI, EVI, NDWI, and LST for different biomes. Noticeably, NDVI is one of driving parameters for DSI, which is a little different from that used for calculating the standardized anomaly of NDVI in this study. As the former one is only for growing season, the latter one is based on whole season.

5. Results and Discussion

5.1. The 2009/2010 Drought Quantified by the MODIS-Derived DSI

Distribution of the drought in southwestern China as detected by DSI is shown in Figure 2. Spatially, western Guizhou, central and eastern Yunnan, and western Guangxi experienced moderate to extreme drought, showing higher drought stress in the drought-stricken regions (Figure 2). Nondrought areas were mainly distributed in the horns of northwestern Sichuan, northeastern Sichuan, and Chongqing, as well as south and east of Guangxi (Figure 2).

As shown in Figure 2, 74.4% of the study area experienced drought whose intensities ranged from incipient to extreme. Percentages of areas that experienced incipient, mild, moderate, severe, and extreme drought were 14.6%, 11.7%, 14.9%, 16.3%, and 16.9%, respectively (Figure 3).

5.2. Spatial Distribution of Standardized Anomalies for NDVI, EVI, NDWI, and LST

Figure 4 shows the spatial distributions of standardized anomalies of NDVI, EV, NDWI, and LST in vegetated areas that suffered from the 2009/2010 drought. The spatial patterns of standardized anomalies varied among NDVI, EVI, NDWI, and LST, as well as among different drought-stricken regions.

For both NDVI and EVI during the drought period, areas experiencing negative anomalies had similar spatial distributions in central and eastern Yunnan, parts of southern Sichuan and western Guizhou, and northern Sichuan, which suffered from moderate to extreme drought (Figures 4(a) and 4(b)). It was implied that the vegetation greenness (NDVI and EVI) declined in those drought-stricken regions during the drought period. We found that the spatial distributions of standardized anomalies of NDVI corresponded favorably with the distributions of drought, in comparison with the EVI for overall distribution trends, implying that standardized anomalies of NDVI were more sensitive to the 2009/2010 drought than standardized anomalies of EVI in southwestern China (Figures 2, 4(a), and 4(b)).

The areas experiencing negative anomalies of NDWI were smaller than those of NDVI and EVI, which were mainly distributed in the eastern Yunnan, western Guizhou, and central Sichuan, which suffered from moderate to extreme drought (Figure 4(c)). In comparison, areas experiencing positive anomalies of NDWI were distributed more widely than those of NDVI and EVI, indicating that the greenness of vegetation (NDVI and EVI) was more easily affected by the 2009/2010 drought than vegetation water content (NDWI), especially in central and western Yunnan (Figures 2, 4(a), 4(b), and 4(c)). It is notable that the results of NDWI were contradictory with NDVI, EVI, and LST in some areas of eastern Sichuan (croplands) and southern Yunnan (deciduous forests). The reasons concerning in this phenomenon were unknown (Figures 4(a), 4(b), and 4(c)).

The areas experiencing positive anomalies for LST included most areas of Yunnan, western Guangzhou, and southern Sichuan, which suffered from mild to extreme drought (Figure 4(d)). It was implied that the LST increased in those drought-stricken areas due to the 2009/2010 drought. Particularly, there were extreme droughts in areas bordering Yunnan and Guizhou, where positive anomalies were greatest for LST (Figures 2 and 4(d)).

5.3. Frequency Distributions of Standardized Anomalies for NDVI, EVI, NDWI, and LST

To examine the impacts of the 2009/2010 drought on vegetation, we generated frequency histograms of standardized anomalies of NDVI, EVI, NDWI, and LST (Figure 5). We found that the standardized anomalies of NDVI, EVI, and NDWI exhibited positively skewed frequency distributions with the peak values between −1 and 0 (Figures 5(a), 5(b), and 5(c)), while the standardized anomalies of LST exhibited a negatively skewed frequency distribution with the peak value between 0 and 1 (Figure 5(d)). Both the vegetation greenness (NDVI and EVI) and vegetation water content (NDWI) declined, while the land surface temperature (LST) increased in the 2009/2010 drought-stricken vegetated areas, implying that the 2009/2010 drought had negative impacts on vegetation in southwestern China.

5.4. Relationships between DSI and Standardized Anomalies of NDVI, EVI, NDWI, and LST

Figure 6 shows the relationships between DSI and standardized anomalies of NDVI, EVI, NDWI, and LST during the drought period. We found that the standardized anomalies of NDVI (, ) (Figure 6(a)), EVI (, ) (Figure 6(b)), and NDWI (, ) (Figure 6(c)) had positive correlations with DSI, while LST (, ) (Figure 6(d)) had a negative correlation with DSI during the drought period in southwestern China, implying that the 2009/2010 drought reduced NDVI, EVI, and NDWI but increased LST in southwestern China. It is notable that the plots when DSI is equal to or higher than −0.29 were masked in order to calculate the standardized anomalies of RS-based indices in vegetated areas, which suffered from the 2009/2010 drought over southwestern China. In other words, only the pixels when DSI is equal to or lower than −0.30 (drought) were included in this study; others were masked.

5.5. Reponses of Vegetation Types to the 2009/2010 Drought

Our study regions were mainly covered by mixed forests, deciduous forest, evergreen forest, savannas, grasslands, and croplands (Figure 1). We intended to analyze the effects of the 2009/2010 drought for different types of vegetation in southwestern China. We found that responses of vegetation to the 2009/2010 drought varied among different biomes.

Forest ecosystems showed complex patterns in their responses to the 2009/2010 drought. Among the mixed forests, deciduous forest and evergreen forest, negative anomalies of NDVI and EVI were the greatest for the mixed forests in central and western Yunnan that suffered from severe to extreme drought (Figures 1, 2, and 4). It was indicated that the 2009/2010 drought-induced greenness of mixed forests declined substantially. The standardized anomalies of NDVI (, ), EVI (, ), and NDWI (, ) had positive correlations with DSI, while LST (, ) had a negative correlation with DSI for mixed forests (Figure 7). More than half of the areas with deciduous forest did not experience the 2009/2010 drought (Figures 1 and 2); the standardized anomalies of NDVI (, ), EVI (, ), and NDWI (, ) had positive correlations with DSI, while no relationship was identified between LST and DSI (Figure 7). In comparison with mixed forests and deciduous forest, evergreen forest exhibited different response patterns to the 2009/2010 drought, and the standardized anomalies of NDWI (, ) had a positive correlation with DSI, while NDVI (, ), EVI (, ), and LST (, ) had negative correlations with DSI (Figure 7). It was suggested that the 2009/2010 drought reduced NDWI but increased NDVI, EVI, and LST for evergreen forest.

Grasslands were mainly distributed in north and west Sichuan, which suffered from incipient to mild drought, except for the horn of northwestern Sichuan (no drought) (Figures 1, 2, and 4). The standardized anomalies of NDVI (, ), EVI (, ), and NDWI (, ) had positive correlations with DSI, while LST (, ) had a negative correlation with DSI for grasslands (Figure 7).

Areas experiencing the greatest negative anomalies of NDWI and positive anomalies of LST were covered by mixes of croplands, savannas, and mixed forests and found bordering Yunnan, Guizhou, and Guangxi, which suffered from mild to extreme drought (Figures 1, 2, and 4). The 2009/2010 drought-induced vegetation water content of croplands, savannas, shrublands, and mixed forests declined, while LST increased in those regions. Standardized anomalies of NDVI (, ), EVI (, ), and NDWI (, ) had positive correlations with DSI, while LST (, ) had a negative correlation with DSI for savannas, implying that the 2009/2010 drought reduced NDVI, EVI, and NDWI but increased LST for savannas (Figure 7). For croplands, standardized anomalies of NDVI (, ), EVI (, ), and NDWI (, ) had positive correlations with DSI, while LST (, ) had a negative correlation with DSI, implying that the water content of croplands was reduced largely due to the 2009/2010 drought in comparison with other biomes (Figure 7).

For different biomes in our study area, the savannas were most affected by the 2009/2010 drought quantified by DSI, which had the highest positive correlations with standardized anomalies of NDVI and EVI, but the highest negative correlation with LST. The largest proportions of areas for savannas experienced declines in NDVI and EVI and an increase in LST. DSI had positive correlations with standardized anomalies of NDVI and EVI, but a lower negative correlation with LST, indicating that evergreen forest was less affected by the 2009/2010 drought. Hence, evergreen forest was resistant to the 2009/2010 drought in southwestern China compared with other biomes. The largest positive correlation between DSI and NDWI was seen for croplands, implying that the largest proportions of areas for croplands experienced a decline in NDWI during the 2009/2010 drought.

Our results suggested that the responses of vegetation to the 2009/2010 drought differed substantially among different biomes. Savannas, croplands, and mixed forests were more vulnerable to the 2009/2010 drought than deciduous forest and grasslands, while evergreen forest was resistant to the 2009/2010 drought in southwestern China. The possible reasons of different drought resistances for vegetation types lie in the following: forests were more resistant to this extreme drought than croplands and savannas because trees were generally able to have great access to ground water with deeper roots.

6. Conclusion

We assessed the severity of the 2009/2010 drought using MODIS-derived DSI data and examined the vegetation response to this extreme event using standardized anomalies of NDVI, EVI, NDWI, and LST derived from MODIS data for the period 2000–2011. We found that () NDVI, EVI, and NDWI decreased, while LST increased in the 2009/2010 drought-stricken vegetated areas in southwestern China, as shown by the skewed frequency distributions of standardized anomalies of NDVI, EVI, NDWI, and LST; () the standardized anomalies of NDVI, EVI, and NDWI had positive correlations with DSI, while LST had a negative correlation with DSI; () responses of vegetation to the 2009/2010 drought varied among biomes: savannas, croplands, and mixed forests were more vulnerable to the 2009/2010 drought than deciduous forest and grasslands, while evergreen forest was resistant to the 2009/2010 drought. We concluded that the 2009/2010 drought had negative impacts on vegetation in southwestern China.

There is a need to further investigate the response of croplands to the 2009/2010 drought using other drought indices, such as the vegetation water thermal stress index (VWTCI) [25], as there is a lack of DSI data for cropland cover in southeastern Sichuan and southwestern Chongqing. Further study should focus on comparing the responses of vegetation indices and LST to drought and soil moisture in dry period and wet period [2632].

Competing Interests

The authors declare no conflict of interests.

Acknowledgments

This study was financially supported by the China Scholarship Council (CSC) and also sponsored by the National Natural Science Foundation of China (nos. 41471460, 41130750), the “One Hundred Talents Program” of the Chinese Academy of Sciences (no. Y5BR011001), and the Natural Science Foundation of Jiangsu Province (BK20151060). The authors thank the Land Processes Distributed Active Archive Center (LP DAAC) and Numerical Terradynamical Simulation Group (NTSG) at the University of Montana for providing the MODIS datasets.