Next Article in Journal
Spatiotemporal Evolution and Driving Forces of Production-Living-Ecological Space in Arid Ecological Transition Zone Based on Functional and Structural Perspectives: A Case Study of the Hexi Corridor
Previous Article in Journal
Sustainable Agriculture and Agri-Food
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Altitude-Shifted Climate Variables Dominate the Drought Effects on Alpine Grasslands over the Qinghai–Tibetan Plateau

1
School of Life Sciences, Guizhou Normal University, Guiyang 550025, China
2
Qiangtang Alpine Grassland Ecosystem Research Station (Jointly Built with Lanzhou University), Tibet Agricultural and Animal Husbandry University, Nyingchi 860100, China
3
Lhasa Plateau Ecosystem Research Station, Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
4
College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100190, China
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(15), 6697; https://doi.org/10.3390/su16156697
Submission received: 20 May 2024 / Revised: 28 July 2024 / Accepted: 31 July 2024 / Published: 5 August 2024

Abstract

:
Drought has broad and deep influences on ecosystem dynamics and functions, particularly considering the lagged and cumulative effects of drought. Yet the individual role of climate variables in mediating such drought effects on vegetation remains largely unknown. Based on the Normalized Difference Vegetation Index (NDVI) and the standard precipitation evapotranspiration index (SPEI), here, we investigated the patterns and mechanisms of drought effects on alpine grasslands in the Qinghai–Tibetan Plateau (QTP) from 1982 to 2015. Drought imposed widespread lagged and cumulative impacts on alpine grasslands with notable spatial heterogeneity, showing that the southwestern and northeastern parts of the plateau were more sensitive and responded quickly to drought. Further, drought effects showed an evident elevation dependence across different grassland types, which could be explained by altitudinal shifts in climatic factors, including temperature and precipitation. Precipitation was the dominant factor in drought effects on alpine meadows, while temperature dominated the drought impacts on the alpine steppes. Such a divergent dominant factor implied that there would be different vegetation responses to future climate change among diverse types of alpine grasslands. To maintain the sustainability of alpine grassland, more effort should be applied to alpine steppes regarding pasture management, particularly in response to extreme drought due to warmer climates in the future.

1. Introduction

In the context of ongoing climate change, climate extremes have caused severe impacts on the structure and function of terrestrial ecosystems [1,2,3,4]. As one of the most devastating extreme weather events, drought imposes a far-reaching influence on ecosystem productivity [5,6,7]. Drought can significantly reduce vegetation productivity by constraining vegetation growth [8,9], causing massive plant mortality [10,11,12]. Additionally, the drought-related risks of reducing vegetation productivity are expected to increase under future warming [13,14]. Vegetation provides a lot of ecosystem services, such as food production and carbon uptake, and plays a crucial role in maintaining ecosystem sustainability. An in-depth understanding of the vegetation response to drought is thus critical to evaluating the current sustainability of ecosystems in response to climate extremes and predicting it under future climate changes.
Generally, drought chronically affects vegetation growth [15,16]. For example, current vegetation dynamics might be driven by the impacts of previous water stress due to the lagged effects of climate variables [17,18]. Moreover, drought can also impose a cumulative effect on vegetation performance since it spans a period of time, leading to persistent water constraints in plants [19,20,21]. Combining such lagged and cumulative effects enables a comprehensive assessment of drought impacts on vegetation. Previous studies have demonstrated that the magnitude of lagged and cumulative effects vary among the water conditions typically represented by the annual standard precipitation-evapotranspiration index (SPEI) [22,23,24,25,26]. Yet the individual role of temperature in mediating drought effects has received less attention. Temperature is a crucial factor in vegetation growth and regulates the vegetation response to external disturbance [27]. Growing evidence has also documented how the response of ecosystems to drought could change a lot under warmer environments [28,29]. Therefore, quantifying the potential contributions of temperature to the variability of drought effects is also imperative to comprehend the underlying mechanisms in the vegetation response to drought.
As a key ecological security barrier, the Qinghai–Tibetan Plateau (QTP) is crucial for regional climate and global energy–water cycles [30,31,32]. At the same time, the QTP is a vulnerable region in regard to ongoing climate change due to a typical highland climate [33,34]. Recently, large-scale and rapid warming (roughly twice the global average) has been observed on the QTP [35,36,37,38]. The evident warming has extended the vegetation growth season length and enhanced vegetation productivity over the QTP [39,40]. Nevertheless, the increasing temperature has also accelerated evaporation, exacerbated water deficiency, and thereby exerted damaging impacts on vegetation growth [41,42]. For instance, several studies have shown that warming triggers the substantial browning of vegetation by reducing water availability in the middle and southwestern QTP [43,44]. Hence, considering the rapidly changing climate, drought is projected to trigger more negative influences on vegetation productivity on the QTP against a warming background [45,46].
Covering roughly two-thirds of the total plateau area, alpine grasslands are the major vegetation types on the QTP and maintain fundamental ecosystem services (such as food production) to local human society [47,48]. Alpine grasslands on the QTP can be classified into diverse types according to their dominant species. Among them, the most widely distributed grassland types are alpine meadows and steppe [49,50]. Alpine meadows are characterized by relatively abundant water availability, while the alpine steppe has suffered from limited precipitation [51,52]. Meanwhile, climate variables, mainly including temperature and precipitation, display obvious elevation-dependent characteristics on the QTP [53,54,55]. Further, the alpine grasslands’ sensitivity to climate change also varies with elevation gradients [56,57,58]. Here, the two widespread grassland types and huge elevational ranges on the QTP provide an opportunity to explore the potential drivers of drought effects on vegetation growth.
In this study, leveraged by the Normalized Difference Vegetation Index (NDVI) and the corresponding standard precipitation evapotranspiration index (SPEI), we calculated drought effects, including lagged and cumulative effects, on QTP grasslands and further investigated their variations along the gradients of elevation and climate variables. Our overall objectives are as follows: (1) to detect the spatial distributions of drought effects on alpine grasslands; (2) to reveal how these effects change along elevation gradients and, further, to identify the predominant climatic contributor.

2. Materials and Methods

2.1. Study Area

The QTP is located in southwest China and has a unique environment with low temperatures, intensified solar radiation, limited precipitation, and complex terrain [59,60]. The elevation of alpine grasslands ranges from less than 2500 m to more than 5500 m (Figure 1a,b). The mean temperature of the growing season from 1982 to 2015 in the grasslands increased from less than 0 °C to greater than 12 °C (Figure 1c). The total precipitation of the growing season ranges from 100 mm to 500 mm (Figure 1d). The growing season is limited to May to September of each year.

2.2. Data

2.2.1. NDVI Dataset

To capture the vegetation dynamics, we used the NDVI as the indicator of vegetation growth. The GIMMS3g NDVI for 1982–2015 was derived at https://poles.tpdc.ac.cn/en/data/9775f2b4-7370-4e5e-a537-3482c9a83d88/, accessed on 15 May 2024 [61]. This biweekly NDVI dataset is one of the longest time series of satellite-based vegetation indicators with an 8 km × 8 km grid resolution. This NDVI dataset minimized the influence of sensor degradation, cloud cover, and volcanic aerosols and showed good performance in temporal consistency. To eliminate the impacts of sparsely vegetated areas, grid cells with a mean NDVI for growing season less than 0.1 were excluded from the analysis. To match the resolution of the climate dataset, the NDVI dataset was resampled to 0.1° grid resolution by adopting the nearest-neighbor algorithm and aggregated into a monthly timescale via maximum value composite.

2.2.2. Climate Data

The monthly mean temperature and total precipitation were obtained from the China meteorological forcing dataset (CMFD) (http://data.tpdc.ac.cn/en/, accessed on 15 May 2024), which was developed by fusing the climate variables of satellite-based observations, reanalysis datasets, and meteorological station data [62]. With a grid resolution of 0.1°, the CMFD is generally better at evaluating accuracy than other climatic products in China.

2.2.3. Vegetation Cover and Elevation Dataset

The grassland types on the QTP were derived from the 1:1,000,000 vegetation atlas of China (https://poles.tpdc.ac.cn/zh-hans/data/eac4f2cf-d527-4140-a35d-79992957f043/, accessed on 15 May 2024), which was completed by more than 200 scientists after more than 30 years of work and has been widely used in previous ecological studies. The vegetation map shows the distribution of 868 basic vegetation taxonomic units in China. It shows, in detail, the distribution and geographical pattern of vegetation, including horizontal distribution, vertical distribution, and its relationship with climatic factors and ground environmental factors. The elevation (spatial resolution of ~1 km) of the Global 30 Arc-Second Elevation Dataset (GTOPO30, https://www.usgs.gov/centers, accessed on 15 May 2024) was downloaded from the Google Earth Engine (GEE) platform. GTOPO30, completed in late 1996, is a global digital elevation model (DEM) derived from several raster and vector sources of topographic information. To match the spatial resolution of the climate data, both the vegetation cover and elevation dataset were resampled to a 0.1° grid resolution by adopting bilinear interpolation.

2.3. Method

2.3.1. Drought Index

The standard precipitation evapotranspiration index (SPEI), derived from the deviations from average climatic water balance (Equation (1)), is widely used in the assessment of the intensity and duration of drought [63]. The SPEI synthesizes the effects of temperature variability into drought evaluation and has the capacity to identify droughts at different scales [19,64]. The calculation method for climatic water balance is as follows:
B i = P i P E T i
where i is the month, and B i , P i , and P E T i are the monthly water balance (mm), the monthly total precipitation (mm), and the monthly potential evapotranspiration (mm), respectively.
To obtain the SPEI values, we first calculated the gridded potential evapotranspiration following the method of Thornthwaite [65] as follows:
P E T i = 16 × 10 T i H A × N 12 × D 30
H = i = 1 12 T i 5 1.154
A = 0.49239 + 1792 × 10 5 H 771 × 10 7 H 2 + 675 × 10 9 H 3
where T i , A, N, D, and H are the mean temperature, constant, maximum sunny hours, number of days per month, and heat index, respectively. It should be noted that the N is estimated by the latitude.
Then, we normalized the monthly aggregated water balance from Equation (1) at multiple timescales (from 1 to 24 months) to obtain the SPEI values (http://spei.csic.es/home.html, accessed on 15 May 2024). The values of the SPEI indicate climatic water availability, with larger positive values suggesting more water surplus and vice versa. Moreover, the SPEI can represent the previous cumulated water availability with multiple timescales. For instance, the values of the 1-month and 3-month SPEIs indicate the water conditions of the current month and the previous three months, respectively.

2.3.2. Lagged and Cumulative Effects of Drought on Vegetation

The lagged and cumulative effects were calculated based on the maximum Pearson’s correlation coefficient between NDVI and SPEI [19,21,22]. Specifically, for the lagged effect, the correlation coefficients between the NDVI and 1-month SPEI among the preceding 24 months were first calculated. For the cumulative effect, the correlation coefficients were derived from the NDVI and the corresponding different timescales (i.e., from 1 to 24 months) of the SPEI. Then, the lagged month-of-drought effect was determined, in which the preceding month where the maximum correlation coefficient between the NDVI and the 1-month SPEI was observed. The accumulated month-of-drought effect was determined, in which the timescale of the SPEI where the maximum correlation coefficient between the NDVI and the SPEI is shown. Meanwhile, the maximum correlation coefficients indicated the strength of the lagged and cumulative effects. Generally, the shorter drought timescales (i.e., lagged and accumulated times) imply that vegetation reacts more quickly to drought, indicating higher sensitivity [66,67].
R i = c o r r N D V I , S P E I i 1 i 24
R m a x _ l a g = max R i 1 i 24
where R i represents the correlation coefficient at the lagged time of the i month, and S P E I i and R m a x _ l a g denote the 1-month timescale SPEI at previous i month and the maximum value of R i , respectively.
R j = c o r r N D V I , S P E I j 1 j 24
R m a x _ c m l = max R j 1 j 24
where R j represents the correlation coefficient between the NDVI and SPEI at a j-month timescale. S P E I j and R m a x _ c m l denote the j-month timescale SPEI and the maximum value of R j , repectively.

2.3.3. Elevation and Climate Gradients

To investigate the elevation and climate gradients of drought effects in grassland vegetation, we calculated the mean lagged and cumulative effects in bins of elevation at an interval of 100 m. Then, the average temperature and total precipitation were summarized in each elevation bin. When calculating drought effects along elevation, only the grid cells with significant correlation coefficients between the NDVI and the SPEI were included, and elevation bins with a sample size of less than 100 were removed.

2.4. Statistical Analysis

Linear or quadratic regression was adopted to estimate the relationships between drought effects and elevation, as well as climatic variables. It should be noted that we chose the above two regression methods due to their explicit ecological implications. For example, the linear regression between temperature and drought effects implies that drought effects increase or decrease along with temperature, while quadratic regression suggests that there is an optimum temperature. Then, the Akaike information criterion (AIC) was used to determine the best regression between linear and quadratic fitting, with the lower values suggesting better performance. Finally, we employed the random forest algorithm (RF) [68] to evaluate the relative contributions of elevation and climatic variables (temperature and precipitation) to the effects of variation in drought effects on grassland vegetation. We used 70% of the data for model training and 30% for model validation. Finally, we built an RF model with 500 regression trees. In the RF, the mean square error (%IncMSE) represents the importance of each predictor (i.e., higher values indicate higher variable importance) (Ho 1998).
All calculations and statistical tests were conducted in R software version 4.3.0 [69]. The “SPEI” package was used to calculate the SPEI value, the “Hmisc” package to calculate Pearson’s correlation, the “stats” package to conduct linear and quadratic regression, and the “randomForest” package to create a random forest model.

3. Results

3.1. Spatial Patterns of Drought Effects

For the lagged effects, about 99.6% of grasslands presented positive correlations, and 54.52% of those were statistically significant (p < 0.05), primarily distributed in the northeastern and southwestern parts of the plateau (Figure 2a). Among the pixels showing lagged effects, about 40.1% of them exhibited shorter lagged times (range of 1 to 4 months), mainly located in the western and northeastern parts of the plateau (Figure 2b). Conversely, approximately 12.2% of the areas showed longer lagged times (range of 21 to 24 months) and were concentrated on the central plateau.
For the cumulative effects, most grassland areas (80.28%) in the QTP were dominated by positive correlation, while 19.71% of areas showed negative correlations (Figure 3a). In addition, 32.01% of grasslands, mostly located in the southwestern and northeastern parts of the plateau, exhibited significantly positive correlations (p < 0.05). Negative correlations, despite almost none meeting the significance criteria, were primarily observed in the northwestern and southeastern parts of the plateau. Considering the accumulated timescale, the short-term (1–4 months) accumulated months principally occupied the southwest and northeast plateau (covering 12.20% of the total plateau areas), whereas medium (13–16 months) and long-term (21–24 months) accumulated months mostly occurred in the central plateau (amounting to 23.30% and 30.27% of total plateau areas, respectively) (Figure 3b).

3.2. Variation in Drought Effects along the Gradients of Elevation, Temperature, and Precipitation

For alpine meadows, a parabolic-shaped relationship was observed between lagged months and elevation (p = 0.001) and a U-shaped relationship between Rmax_lag and elevation (p = 0.001) (Figure 4a; Supplementary Figure S1a). The lagged months first increased at altitudes below 4300 m and then decreased at altitudes greater than 4300 m. However, there are no relationships between accumulated months/Rmax_cml and elevation (Figure 4c, Supplementary Figure S1c). For the alpine steppe, the lagged (p < 0.001) and accumulated (p < 0.001) months were significantly positively correlated with elevation (Figure 4b,d). Additionally, Rmax_lag (R2 = 0.72, p < 0.001) and Rmax_cml (R2 = 0.72, p < 0.001) decreased significantly along the increasing elevation. These results suggest that the alpine steppes are more susceptible and respond more quickly to drought at lower altitudes (Supplementary Figure S1b,d).
For alpine meadows, a parabolic-shaped relationship was observed between lagged months and temperature (p = 0.001) and a U-shaped relationship for Rmax_lag and temperature (p = 0.001) (Figure 5a; Supplementary Figure S2a). The lagged months first increased at temperatures below 6.1 °C and then decreased at altitudes greater than 6.1 °C. However, the accumulated months and Rmax_cml did not show significant linear or quadratic changes along the temperature gradients (Figure 5c; Supplementary Figure S2c). For the alpine steppe, the lagged (p < 0.001) and accumulated (p < 0.001) months were significantly negatively correlated with temperature (Figure 5b,d). Meanwhile, Rmax_lag (p < 0.001) and Rmax_cml (p < 0.001) significantly increased along with increasing temperatures. Both of these findings indicate that alpine steppes are more sensitive and respond more rapidly to drought in hotter areas (Supplementary Figure S2b,d).
For the alpine meadow, the lagged months were significantly positively correlated with precipitation during the growing season (p < 0.001), while Rmax_lag was significantly negatively correlated (p < 0.001) (Figure 6a; Supplementary Figure S3a). Meanwhile, a parabolic-shaped relationship was exhibited between the accumulated months and precipitation (p < 0.01), while there was a U-shaped relationship for Rmax_cml and precipitation (p < 0.05) (Figure 6c, Supplementary Figure S3c). The accumulated months first increased at temperatures below 385 mm and then decreased at altitudes greater than 385 mm. For the alpine steppe, the lagged months, accumulated months, Rmax_lag, and Rmax_cml showed no significant relationships with precipitation (Figure 6b,d; Supplementary Figure S3b,d).

3.3. Divergent Contributions of Elevation, Temperature, and Precipitation to Drought Effects

Overall, the random forest model showed good performance in simulating the observed variables regarding drought effects on the alpine meadow and steppe (Supplementary Figures S4–S6). As for alpine meadows, elevation contributed the most variations in lagged effects and accumulated months (Figure 7a,c). Among the climatic predictors, precipitation contributed more variations in accumulated months than temperature. For the alpine steppe, elevation contributed the most variations in lagged and accumulated months (Figure 7b,d). Among the climatic predictors, temperature was a more important variable than precipitation. In addition, the results of Rmax_lag and Rmax_cml were similar to the findings for the lagged and accumulated months (Supplementary Figures S5 and S6).

4. Discussion

4.1. The Spatial Patterns of Drought Effects on Alpine Grassland

Overall, alpine grasslands were dominated by positive correlation coefficients in regard to the lagged and cumulative effects, indicating vegetation growth on grasslands is substantially influenced by previous climatic water conditions (Figure 2a and Figure 3a). The underlying mechanisms may be related to the unique climatic environment of alpine grasslands and the functional traits of herbaceous plants. First, alpine grasslands mainly belong to arid and semi-arid regions [70]. As they are characterized by infrequent and highly variable precipitation, water availability is a critical factor in affecting vegetation growth in alpine grasslands [71,72]. Second, the grasses merely capture the water from the surface of soil layers due to their relatively short and shallow roots [73,74]. Such root characteristics prevent the herbs from absorbing the deep soil water to escape drought stress, leading to a strong relationship between vegetation growth and water availability in alpine grasslands.
Additionally, we found that the drought effects were spatially heterogeneous. On the one hand, the significant positive correlation coefficients were primarily distributed in the southwestern and northeastern parts of the plateau (Figure 2a and Figure 3a). On the other hand, the short-term (1–4 months) lagged/accumulated months principally occupied the southwestern and northeastern parts of the plateau, whereas medium (13–16 months) and long-term (21–24 months) lagged/accumulated months mostly occurred in the central plateau (Figure 2b and Figure 3b). These results indicated that the southwestern and northeastern parts of the plateau are more susceptible and respond more quickly to climatic water deficits. This phenomenon may be attributable to the distribution patterns of different grassland types. Spatially, alpine steppes are largely aggregated in the western and northern parts of the QTP, while alpine meadows are mainly distributed in the central and eastern parts of the plateau [70]. Alpine steppes generally exhibit higher vulnerability to drought stress than alpine meadows [29].

4.2. Divergent Variation in Drought Effects along Elevation and Climate Gradients among Different Grassland Types

Generally, the drought effects on alpine grasslands varied with elevation, which can be attributed to altitudinal shifts in climate variables (e.g., temperature and precipitation). For alpine meadows, the magnitude of lagged months first increased at lower altitudes (below 4300 m) and then decreased at higher altitudes (greater than 4300 m) (Figure 4a). This phenomenon was consistent with the finding that the drought effect generally increased with increasing elevation but showed a decline when the altitude ranged from 3000 to 5000 m in the global mountain areas [21]. Meanwhile, the variation in lagged months exhibited a parabolic-shaped relationship with temperature gradients and a positive relationship with precipitation (Figure 5a and Figure 6a). Furthermore, the random forest results imply that water availability is the dominant factor of lagged months. These results suggest that alpine meadows are sensitive to drought in more arid conditions (Figure 6a). Sufficient precipitation could serve as a buffer by increasing the soil moisture content for vegetation growth when drought occurs [20,75]. Meanwhile, the coupling of vegetation dynamics and water availability could be strengthened when ecosystems experience declined precipitation [66,76]. Therefore, the magnitude of drought effects on alpine meadows varies along with the water conditions, which was also corroborated by previous studies [23,26]. For instance, previous studies have revealed that the drought effects on global grasslands were stronger in arid regions than in humid areas. Meanwhile, there is also numerical evidence that global vegetation tends to be more tolerant to drought in humid ecosystems.
For the alpine steppe, drought effects showed linear altitude dependencies, as demonstrated by lagged and accumulated months rapidly increasing along elevation gradients (Figure 4b,d). These results suggested that alpine steppes are more susceptible to water deficits and respond more quickly to short-term drought at lower altitudes, which is in line with previous evidence showing that drought effects decrease with increasing elevation [29,77]. At the same time, positive relationships were found in the variation in lagged and accumulated months along with increasing temperature, indicating a more pronounced drought effect on alpine steppes under warmer conditions (Figure 5b,d). This result was in contrast to previous studies that typically revealed that water availability, not temperature, is the dominant factor affecting the vegetation response to drought. Our study object, that is, the alpine steppe, is the primary reason for this contrasting result. Compared with the other ecosystems investigated in previous studies, the alpine steppe is more likely to suffer from the limitations of energy and temperature. Thus, temperature plays a critical role in mediating vegetation productivity and, therefore, the vegetation response to drought. A recent study found that the magnitude of vegetation growth responses to drought is mainly influenced by temperature in the QTP, particularly at 3000–4500 m intervals [29,78]. Moreover, it has also been documented that drought accompanied by warmer temperatures causes more damage to plants compared with simple water deficits [79,80,81]. Additionally, there were non-significant relationships between drought effects and precipitation (Figure 6b,d). Largely located in hyper-arid and arid areas, alpine steppes suffer from chronic water stress [82]. We speculate that the physiological and functional characteristics of plants in alpine steppes may evolve to adapt to the arid environment, leading to high resilience and resistance to precipitation variations [83,84]. Finally, the results of the random forest model also showed that temperature is a more important variable than precipitation for drought effects (Figure 7b,d). Therefore, against the background of rapid warming, the alpine steppe is expected to be more vulnerable to water deficiency since increasing temperatures will shorten its response time to drought.

4.3. Limitations and Implications for Future Studies

Here, we detected the spatial patterns of drought effects on alpine grassland and further analyzed how elevation and climatic variables affect such drought effects among various grassland types. Specifically, we found that precipitation was the dominant factor in drought effects on alpine meadows, whereas temperature is the main driving factor for the responses of alpine steppe vegetation to drought. Such a divergent dominant factor implies that there will be different vegetation responses to future climate change among the two types of alpine grassland. Considering that the QTP is projected to become warmer and wetter, drought effects on vegetation could be exacerbated in the alpine steppe due to elevated temperatures and be alleviated in alpine meadows because of increased precipitation. Thus, our study highlighted that we should pay more attention and apply more effort to the alpine steppe in response to future droughts. Nevertheless, our study is still encumbered by some limitations and uncertainties, which hamper our comprehensive understanding of the vegetation response to drought. First, the spatial heterogeneity of the vegetation growth response to drought in different elevation intervals needs further exploration, as there are varying hydrothermal conditions at various elevations. Multiple drought indexes should be incorporated to enhance the robustness of the results, for example, the normalized Reconnaissance Drought Index (RDI). It was shown that there was no significant influence on RDI results detected by four popular PET methods (Hargreaves, Thornthwaite, Blaney–Criddle, and FAO Penman–Monteith (only T)), regardless of the reference period analysis. Second, we argue that the apparent spatial heterogeneity in drought effects on the alpine grasslands was associated with different hydrothermal conditions and diverse functional strategies for different grassland types. Other environmental factors, such as soil properties and biodiversity [85,86,87], may also alter the vegetation response to drought. For instance, the high species diversity and complex community structure could enhance resistance to drought stresses, as illustrated by forests with higher species diversity that are less affected by changing water conditions [88,89,90]. Additionally, we also did not take the impact of ecological disturbances on vegetation growth into consideration. Ecological disturbances, including land use changes and human activities [35,59], may pronouncedly affect the growth of vegetation and thereby reshape the influence of drought on the alpine grasslands. Thus, future studies may account for more environmental indicators to identify potential control factors in shaping drought effects. Third, due to the limited large-scale control experiments conducted on the QTP, we employed statistical methods (including regression analysis and random forest arithmetic) to explore the relationship between drought effects and their possible driving forces. However, more field drought experiments are urgently needed in the future, which would help us to reveal the physiological and ecological mechanisms of alpine vegetation in response to drought.

5. Conclusions

Quantifying the role of climate variables in shaping the vegetation response to drought is important for a comprehensive understanding of drought impacts. Using the NDVI and corresponding SPEI, we detected the spatial patterns of drought effects on alpine grasslands and analyzed their variations along the gradients of elevation and climate variables across different grassland types in the QTP. In general, drought imposed widespread lagged and cumulative effects on alpine grasslands. Meanwhile, the stronger sensitivity and short-term response time of vegetation were primarily distributed in the southwestern and northeastern parts of the plateau. Further, drought effects showed evident elevation dependence across different grassland types, which could be explained by altitudinal shifts in climate factors, including temperature and precipitation. Precipitation dominated the response of alpine meadows to drought, whereas temperature was the dominant driving factor of the vegetation response to drought on the alpine steppe. Our findings highlighted the divergent effects of hydrothermal conditions on different alpine grasslands in response to drought, which further implied a divergent change when facing future climate change. Additionally, more drought-related experiments are urgently needed to reveal the physiological and ecological mechanisms of alpine vegetation under drought.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/su16156697/s1, Table S1. The Akaike information criterion (AIC) of for different regression methods. Figure S1. Relationships between the maximum correlation coefficients (Rmax_lag and Rmax_cml) and elevation among different grassland types. Relationships between Rmax_lag and elevation in alpine meadows (a) and the alpine steppe (b); relationships between Rmax_cml and elevation in alpine meadows (c) and the alpine steppe (d). Dots represent the average values of the corresponding elevation bins (at 100 m intervals). Rmax_lag represents the maximum correlation coefficients between the NDVI and the 1-month SPEI in the preceding 1 to 24 months. Rmax_cml represents the maximum correlation coefficients between the 1-to-24-month SPEI and the NDVI. The error bars and shade denote the standard errors and 95% confidence interval of regression fitting. Figure S2. Relationships between the maximum correlation coefficients (Rmax_lag and Rmax_cml) and temperature among different grassland types. Relationships between Rmax_lag and temperature in alpine meadows (a) and the alpine steppe (b); relationships between Rmax_cml and temperature in alpine meadows (c) and on the alpine steppe (d). Dots represent the average values of the corresponding elevation bins (at 100 m intervals). The error bars and shade denote the standard errors and 95% confidence interval of regression fitting. Figure S3. Relationships between the maximum correlation coefficients (Rmax_lag and Rmax_cml) and precipitation among different grassland types. Relationships between Rmax_lag and precipitation in alpine meadows (a) and the alpine steppe (b); relationships between Rmax_cml and precipitation in alpine meadows (c) and the alpine steppe (d). Dots represent the average values of the corresponding elevation bins (at 100 m intervals). The error bars and shade denote the standard errors and 95% confidence interval of regression fitting. Figure S4. The performance of the random forest model. The relationships between the observed lagged months and simulated lagged months in alpine meadows (a) and the alpine steppe (b). The relationships between the observed accumulated months and simulated accumulated months in alpine meadows (c) and the alpine steppe (d). Figure S5. The performance of the random forest model. The relationships between observed Rmax_lag and simulated Rmax_lag in alpine meadows (a) and the alpine steppe (c). Relative contributions of elevation, temperature, and precipitation to Rmax_lag in alpine meadows (b) and the alpine steppe (d). Figure S6. The performance of the random forest model. The relationships between observed Rmax_cml and simulated Rmax_cml in alpine meadows (a) and the alpine steppe (c). Relative contributions of elevation, temperature, and precipitation to Rmax_cml in alpine meadows (b) and the alpine steppe (d). Figure S7. Scatterplot of residuals of different regressions for alpine meadows. Figure S8. Scatterplot of residuals of different regressions for the alpine steppe.

Author Contributions

Conceptualization, X.W.; Writing—original draft, X.W.; Writing—review & editing, Z.H., Z.Z., J.T. and B.N.; Funding acquisition, X.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was jointly supported by National Natural Sciences Foundation of China (42161012), National Key R&D Program of China (2023YFF1304300), and the Tibet Autonomous Region Science and Technology Planning Project (XZ202201ZY0016G).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in [GIMMS3g NDVI] [https://poles.tpdc.ac.cn/en/data/9775f2b4-7370-4e5e-a537-3482c9a83d88/, accessed on 15 May 2024.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Frank, D.; Reichstein, M.; Bahn, M.; Thonicke, K.; Frank, D.; Mahecha, M.D.; Smith, P.; van der Velde, M.; Vicca, S.; Babst, F.; et al. Effects of climate extremes on the terrestrial carbon cycle: Concepts, processes and potential future impacts. Glob. Chang. Biol. 2015, 21, 2861–2880. [Google Scholar] [CrossRef] [PubMed]
  2. Masson-Delmotte, V.; Zhai, P.; Pirani, A.; Connors, S.L.; Péan, C.; Berger, S.; Caud, N.; Chen, Y.; Goldfarb, L.; Gomis, M. Climate Change 2021: The Physical Science Basis; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2021; Volume 2. [Google Scholar]
  3. Piao, S.; Zhang, X.; Chen, A.; Liu, Q.; Lian, X.; Wang, X.; Peng, S.; Wu, X. The impacts of climate extremes on the terrestrial carbon cycle: A review. Sci. China Earth Sci. 2019, 62, 1551–1563. [Google Scholar] [CrossRef]
  4. Reichstein, M.; Bahn, M.; Ciais, P.; Frank, D.; Mahecha, M.D.; Seneviratne, S.I.; Zscheischler, J.; Beer, C.; Buchmann, N.; Frank, D.C.; et al. Climate extremes and the carbon cycle. Nature 2013, 500, 287–295. [Google Scholar] [CrossRef] [PubMed]
  5. Dai, A.G. Increasing drought under global warming in observations and models. Nat. Clim. Chang. 2013, 3, 52–58. [Google Scholar] [CrossRef]
  6. Jiang, D.; Wang, X. A brief interpretation of drought change from IPCC Sixth Assessment Report. Trans. Atmos. Sci. 2021, 44, 650–653. [Google Scholar]
  7. Trenberth, K.E.; Dai, A.; van der Schrier, G.; Jones, P.D.; Barichivich, J.; Briffa, K.R.; Sheffield, J. Global warming and changes in drought. Nat. Clim. Chang. 2013, 4, 17–22. [Google Scholar] [CrossRef]
  8. Ciais, P.; Reichstein, M.; Viovy, N.; Granier, A.; Ogee, J.; Allard, V.; Aubinet, M.; Buchmann, N.; Bernhofer, C.; Carrara, A.; et al. Europe-wide reduction in primary productivity caused by the heat and drought in 2003. Nature 2005, 437, 529–533. [Google Scholar] [CrossRef] [PubMed]
  9. Esquivel-Muelbert, A.; Baker, T.R.; Dexter, K.G.; Lewis, S.L.; Brienen, R.J.W.; Feldpausch, T.R.; Lloyd, J.; Monteagudo-Mendoza, A.; Arroyo, L.; Alvarez-Davila, E.; et al. Compositional response of Amazon forests to climate change. Glob. Chang. Biol. 2019, 25, 39–56. [Google Scholar] [CrossRef] [PubMed]
  10. McDowell, N.; Pockman, W.T.; Allen, C.D.; Breshears, D.D.; Cobb, N.; Kolb, T.; Plaut, J. Mechanisms of plant survival and mortality during drought: Why do some plants survive while others succumb to drought? New Phytol. 2008, 178, 719–739. [Google Scholar] [CrossRef]
  11. McDowell, N.G.; Sapes, G.; Pivovaroff, A.; Adams, H.D.; Allen, C.D.; Anderegg, W.R.L.; Arend, M.; Breshears, D.D.; Brodribb, T.; Choat, B.; et al. Mechanisms of woody-plant mortality under rising drought, CO2 and vapour pressure deficit. Nat. Rev. Earth Environ. 2022, 3, 294–308. [Google Scholar] [CrossRef]
  12. van der Molen, M.K.; Dolman, A.J.; Ciais, P.; Eglin, T.; Gobron, N.; Law, B.E.; Meir, P.; Peters, W.; Phillips, O.L.; Reichstein, M.; et al. Drought and ecosystem carbon cycling. Agric. For. Meteorol. 2011, 151, 765–773. [Google Scholar] [CrossRef]
  13. Xu, C.; McDowell, N.G.; Fisher, R.A.; Wei, L.; Sevanto, S.; Christoffersen, B.O.; Weng, E.; Middleton, R.S. Increasing impacts of extreme droughts on vegetation productivity under climate change. Nat. Clim. Chang. 2019, 9, 948–953. [Google Scholar] [CrossRef]
  14. Zeng, Z.; Wu, W.; Li, Y.; Huang, C.; Zhang, X.; Peñuelas, J.; Zhang, Y.; Gentine, P.; Li, Z.; Wang, X.; et al. Increasing meteorological drought under climate change reduces terrestrial ecosystem productivity and carbon storage. One Earth 2023, 6, 1326–1339. [Google Scholar] [CrossRef]
  15. Vicente-Serrano, S.M.; Quiring, S.M.; Pena-Gallardo, M.; Yuan, S.; Dominguez-Castro, F.J.E.-S.R. A review of environmental droughts: Increased risk under global warming? Earth-Sci. Rev. 2020, 201, 102953. [Google Scholar] [CrossRef]
  16. Wu, D.; Zhao, X.; Liang, S.; Zhou, T.; Huang, K.; Tang, B.; Zhao, W. Time-lag effects of global vegetation responses to climate change. Glob. Chang. Biol. 2015, 21, 3520–3531. [Google Scholar] [CrossRef] [PubMed]
  17. Kannenberg, S.A.; Schwalm, C.R.; Anderegg, W.R.L. Ghosts of the past: How drought legacy effects shape forest functioning and carbon cycling. Ecol. Lett. 2020, 23, 891–901. [Google Scholar] [CrossRef] [PubMed]
  18. Muller, L.M.; Bahn, M. Drought legacies and ecosystem responses to subsequent drought. Glob. Chang. Biol. 2022, 28, 5086–5103. [Google Scholar] [CrossRef] [PubMed]
  19. Vicente-Serrano, S.M.; Gouveia, C.; Camarero, J.J.; Begueria, S.; Trigo, R.; Lopez-Moreno, J.I.; Azorin-Molina, C.; Pasho, E.; Lorenzo-Lacruz, J.; Revuelto, J.; et al. Response of vegetation to drought time-scales across global land biomes. Proc. Natl. Acad. Sci. USA 2013, 110, 52–57. [Google Scholar] [CrossRef] [PubMed]
  20. Xu, H.-J.; Wang, X.-P.; Zhao, C.-Y.; Yang, X.-M. Diverse responses of vegetation growth to meteorological drought across climate zones and land biomes in northern China from 1981 to 2014. Agric. For. Meteorol. 2018, 262, 1–13. [Google Scholar] [CrossRef]
  21. Zhang, Z.; Ju, W.; Zhou, Y.; Li, X. Revisiting the cumulative effects of drought on global gross primary productivity based on new long-term series data (1982–2018). Glob. Chang. Biol. 2022, 28, 3620–3635. [Google Scholar] [CrossRef]
  22. Peng, J.; Wu, C.; Zhang, X.; Wang, X.; Gonsamo, A. Satellite detection of cumulative and lagged effects of drought on autumn leaf senescence over the Northern Hemisphere. Glob. Chang. Biol. 2019, 25, 2174–2188. [Google Scholar] [CrossRef] [PubMed]
  23. Wei, X.; He, W.; Zhou, Y.; Ju, W.; Xiao, J.; Li, X.; Liu, Y.; Xu, S.; Bi, W.; Zhang, X.; et al. Global assessment of lagged and cumulative effects of drought on grassland gross primary production. Ecol. Indic. 2022, 136, 108646. [Google Scholar] [CrossRef]
  24. Zhan, C.; Liang, C.; Zhao, L.; Jiang, S.; Niu, K.; Zhang, Y. Drought-related cumulative and time-lag effects on vegetation dynamics across the Yellow River Basin, China. Ecol. Indic. 2022, 143, 109409. [Google Scholar] [CrossRef]
  25. Zhang, Q.; Kong, D.; Singh, V.P.; Shi, P. Response of vegetation to different time-scales drought across China: Spatiotemporal patterns, causes and implications. Glob. Planet. Chang. 2017, 152, 1–11. [Google Scholar] [CrossRef]
  26. Zhao, A.; Yu, Q.; Feng, L.; Zhang, A.; Pei, T. Evaluating the cumulative and time-lag effects of drought on grassland vegetation: A case study in the Chinese Loess Plateau. J. Environ. Manag. 2020, 261, 110214. [Google Scholar] [CrossRef] [PubMed]
  27. Huang, M.T.; Piao, S.L.; Ciais, P.; Penuelas, J.; Wang, X.H.; Keenan, T.F.; Peng, S.S.; Berry, J.A.; Wang, K.; Mao, J.F.; et al. Air temperature optima of vegetation productivity across global biomes. Nat. Ecol. Evol. 2019, 3, 772–779. [Google Scholar] [CrossRef] [PubMed]
  28. Liu, Z.; Zhu, J.; Xia, J.; Huang, K. Declining resistance of vegetation productivity to droughts across global biomes. Agric. For. Meteorol. 2023, 340, 109602. [Google Scholar] [CrossRef]
  29. Wang, Y.; Fu, B.; Liu, Y.; Li, Y.; Feng, X.; Wang, S. Response of vegetation to drought in the Tibetan Plateau: Elevation differentiation and the dominant factors. Agric. For. Meteorol. 2021, 306, 108468. [Google Scholar] [CrossRef]
  30. Immerzeel, W.W.; Lutz, A.F.; Andrade, M.; Bahl, A.; Biemans, H.; Bolch, T.; Hyde, S.; Brumby, S.; Davies, B.J.; Elmore, A.C.; et al. Importance and vulnerability of the world’s water towers. Nature 2020, 577, 364–369. [Google Scholar] [CrossRef]
  31. Qiu, J. The third pole. Nature 2008, 454, 393–396. [Google Scholar] [CrossRef]
  32. Yao, T.; Bolch, T.; Chen, D.; Gao, J.; Immerzeel, W.; Piao, S.; Su, F.; Thompson, L.; Wada, Y.; Wang, L.; et al. The imbalance of the Asian water tower. Nat. Rev. Earth Environ. 2022, 3, 618–632. [Google Scholar] [CrossRef]
  33. Piao, S.; Zhang, X.; Wang, T.; Liang, E.; Wang, S.; Zhu, J.; Niu, B. Responses and feedback of the Tibetan Plateau’s alpine ecosystem to climate change. Chin. Sci. Bull. 2019, 64, 2842–2855. [Google Scholar] [CrossRef]
  34. Zhang, X.; Yang, Y.; Piao, S.; Bao, W.; Wang, S.; Wang, G.; Sun, H.; Luo, T.; Zhang, Y.; Shi, P.; et al. Ecological change on the Tibetan Plateau. Chin. Sci. Bull. 2015, 60, 3048–3056. [Google Scholar]
  35. Chen, H.; Zhu, Q.A.; Peng, C.H.; Wu, N.; Wang, Y.F.; Fang, X.Q.; Gao, Y.H.; Zhu, D.; Yang, G.; Tian, J.Q.; et al. The impacts of climate change and human activities on biogeochemical cycles on the Qinghai-Tibetan Plateau. Glob. Chang. Biol. 2013, 19, 2940–2955. [Google Scholar] [CrossRef] [PubMed]
  36. Kang, S.C.; Xu, Y.W.; You, Q.L.; Flugel, W.A.; Pepin, N.; Yao, T.D. Review of climate and cryospheric change in the Tibetan Plateau. Environ. Res. Lett. 2010, 5, 8. [Google Scholar] [CrossRef]
  37. Kuang, X.X.; Jiao, J.J. Review on climate change on the Tibetan Plateau during the last half century. J. Geophys. Res.-Atmos. 2016, 121, 3979–4007. [Google Scholar] [CrossRef]
  38. Yao, T.D.; Xue, Y.K.; Chen, D.L.; Chen, F.H.; Thompson, L.; Cui, P.; Koike, T.; Lau, W.K.M.; Lettenmaier, D.; Mosbrugger, V.; et al. Recent Third Pole’s Rapid Warming Accompanies Cryospheric Melt and Water Cycle Intensification and Interactions between Monsoon and Environment: Multidisciplinary Approach with Observations, Modeling, and Analysis. Bull. Am. Meteorol. Soc. 2019, 100, 423–444. [Google Scholar] [CrossRef]
  39. Shen, M.G.; Wang, S.P.; Jiang, N.; Sun, J.P.; Cao, R.Y.; Ling, X.F.; Fang, B.; Zhang, L.; Zhang, L.H.; Xu, X.Y.; et al. Plant phenology changes and drivers on the Qinghai-Tibetan Plateau. Nat. Rev. Earth Environ. 2022, 3, 717. [Google Scholar] [CrossRef]
  40. Zhang, G.L.; Zhang, Y.J.; Dong, J.W.; Xiao, X.M. Green-up dates in the Tibetan Plateau have continuously advanced from 1982 to 2011. Proc. Natl. Acad. Sci. USA 2013, 110, 4309–4314. [Google Scholar] [CrossRef] [PubMed]
  41. Fu, G.; Shen, Z.X.; Zhang, X.Z. Increased precipitation has stronger effects on plant production of an alpine meadow than does experimental warming in the Northern Tibetan Plateau. Agric. For. Meteorol. 2018, 249, 11–21. [Google Scholar] [CrossRef]
  42. Ganjurjav, H.; Hu, G.Z.; Zhang, Y.; Gornish, E.S.; Yu, T.Q.; Gao, Q.Z. Warming tends to decrease ecosystem carbon and water use efficiency in dissimilar ways in an alpine meadow and a cultivated grassland in the Tibetan Plateau. Agric. For. Meteorol. 2022, 323, 9. [Google Scholar] [CrossRef]
  43. Li, P.L.; Hu, Z.M.; Liu, Y.W. Shift in the trend of browning in Southwestern Tibetan Plateau in the past two decades. Agric. For. Meteorol. 2020, 287, 9. [Google Scholar] [CrossRef]
  44. Liu, Y.; Li, Z.; Chen, Y. Continuous warming shift greening towards browning in the Southeast and Northwest High Mountain Asia. Sci. Rep. 2021, 11, 17920. [Google Scholar] [CrossRef]
  45. Liu, D.; Wang, T.; Yang, T.; Yan, Z.J.; Liu, Y.W.; Zhao, Y.T.; Piao, S.L. Deciphering impacts of climate extremes on Tibetan grasslands in the last fifteen years. Sci. Bull. 2019, 64, 446–454. [Google Scholar] [CrossRef]
  46. Liu, L.; Wang, Y.; Wang, Z.; Li, D.; Zhang, Y.; Qin, D.; Li, S. Elevation-dependent decline in vegetation greening rate driven by increasing dryness based on three satellite NDVI datasets on the Tibetan Plateau. Ecol. Indic. 2019, 107, 105569. [Google Scholar] [CrossRef]
  47. Cui, X.F.; Graf, H.F. Recent land cover changes on the Tibetan Plateau: A review. Clim. Chang. 2009, 94, 47–61. [Google Scholar] [CrossRef]
  48. Harris, R.B. Rangeland degradation on the Qinghai-Tibetan plateau: A review of the evidence of its magnitude and causes. J. Arid. Environ. 2010, 74, 1–12. [Google Scholar] [CrossRef]
  49. Wang, Z.P.; Wu, J.S.; Niu, B.; He, Y.T.; Zu, J.X.; Li, M.; Zhang, X.Z. Vegetation Expansion on the Tibetan Plateau and Its Relationship with Climate Change. Remote Sens. 2020, 12, 4150. [Google Scholar] [CrossRef]
  50. Wang, Z.P.; Zhang, X.Z.; Niu, B.; Zheng, Y.P.; He, Y.T.; Cao, Y.A.; Feng, Y.F.; Wu, J.S. Divergent Climate Sensitivities of the Alpine Grasslands to Early Growing Season Precipitation on the Tibetan Plateau. Remote Sens. 2022, 14, 2484. [Google Scholar] [CrossRef]
  51. Li, M.; Wu, J.; He, Y.; Wu, L.; Niu, B.; Song, M.; Zhang, X. Dimensionality of grassland stability shifts along with altitudes on the Tibetan Plateau. Agric. For. Meteorol. 2020, 291, 108080. [Google Scholar] [CrossRef]
  52. Wu, J.S.; Wurst, S.; Zhang, X.Z. Plant functional trait diversity regulates the nonlinear response of productivity to regional climate change in Tibetan alpine grasslands. Sci. Rep. 2016, 6, 35649. [Google Scholar] [CrossRef] [PubMed]
  53. Li, X.; Wang, L.; Guo, X.; Chen, D. Does summer precipitation trend over and around the Tibetan Plateau depend on elevation? Int. J. Climatol. 2017, 37, 1278–1284. [Google Scholar] [CrossRef]
  54. Yao, J.; Yang, Q.; Mao, W.; Zhao, Y.; Xu, X. Precipitation trend–Elevation relationship in arid regions of the China. Glob. Planet. Chang. 2016, 143, 1–9. [Google Scholar] [CrossRef]
  55. You, Q.; Chen, D.; Wu, F.; Pepin, N.; Cai, Z.; Ahrens, B.; Jiang, Z.; Wu, Z.; Kang, S.; AghaKouchak, A. Elevation dependent warming over the Tibetan Plateau: Patterns, mechanisms and perspectives. Earth-Sci. Rev. 2020, 210, 103349. [Google Scholar] [CrossRef]
  56. Li, L.; Zhang, Y.; Wu, J.; Li, S.; Zhang, B.; Zu, J.; Zhang, H.; Ding, M.; Paudel, B. Increasing sensitivity of alpine grasslands to climate variability along an elevational gradient on the Qinghai-Tibet Plateau. Sci. Total Environ. 2019, 678, 21–29. [Google Scholar] [CrossRef] [PubMed]
  57. Piao, S.; Cui, M.; Chen, A.; Wang, X.; Ciais, P.; Liu, J.; Tang, Y. Altitude and temperature dependence of change in the spring vegetation green-up date from 1982 to 2006 in the Qinghai-Xizang Plateau. Agric. For. Meteorol. 2011, 151, 1599–1608. [Google Scholar] [CrossRef]
  58. Shen, M.; Zhang, G.; Cong, N.; Wang, S.; Kong, W.; Piao, S. Increasing altitudinal gradient of spring vegetation phenology during the last decade on the Qinghai–Tibetan Plateau. Agric. For. Meteorol. 2014, 189–190, 71–80. [Google Scholar] [CrossRef]
  59. Chen, B.; Zhang, X.; Tao, J.; Wu, J.; Wang, J.; Shi, P.; Zhang, Y.; Yu, C. The impact of climate change and anthropogenic activities on alpine grassland over the Qinghai-Tibet Plateau. Agric. For. Meteorol. 2014, 189–190, 11–18. [Google Scholar] [CrossRef]
  60. Yao, T.; Wu, G.; Xu, B.; Wang, W.; Gao, J.; An, B. Asian Water Tower Change and Its Impacts. Bull. Chin. Acad. Sci. 2019, 34, 1203–1209. [Google Scholar]
  61. Pinzon, J.E.; Tucker, C.J. A Non-Stationary 1981-2012 AVHRR NDVI3g Time Series. Remote Sens. 2014, 6, 6929–6960. [Google Scholar] [CrossRef]
  62. He, J.; Yang, K.; Tang, W.J.; Lu, H.; Qin, J.; Chen, Y.Y.; Li, X. The first high-resolution meteorological forcing dataset for land process studies over China. Sci. Data 2020, 7, 25. [Google Scholar] [CrossRef] [PubMed]
  63. Vicente-Serrano, S.M.; Begueria, S.; Lopez-Moreno, J.I. A Multiscalar Drought Index Sensitive to Global Warming: The Standardized Precipitation Evapotranspiration Index. J. Clim. 2010, 23, 1696–1718. [Google Scholar] [CrossRef]
  64. Begueria, S.; Vicente-Serrano, S.M.; Reig, F.; Latorre, B. Standardized precipitation evapotranspiration index (SPEI) revisited: Parameter fitting, evapotranspiration models, tools, datasets and drought monitoring. Int. J. Climatol. 2014, 34, 3001–3023. [Google Scholar] [CrossRef]
  65. Thornthwaite, C.W. An approach toward a rational classification of climate. Geogr. Rev. 1948, 38, 55–94. [Google Scholar] [CrossRef]
  66. Jiao, W.; Wang, L.; Smith, W.K.; Chang, Q.; Wang, H.; D’Odorico, P. Observed increasing water constraint on vegetation growth over the last three decades. Nat. Commun. 2021, 12, 3777. [Google Scholar] [CrossRef] [PubMed]
  67. Li, D.; An, L.; Zhong, S.; Shen, L.; Wu, S. Declining coupling between vegetation and drought over the past three decades. Glob. Chang. Biol. 2024, 30, e17141. [Google Scholar] [CrossRef] [PubMed]
  68. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  69. R Core TeamR: A Language and Environment for Statistical Computing. 2023. Available online: https://www.R-project.org/ (accessed on 20 January 2024).
  70. Zheng, D. The system of physico-geographical regions of the Qinghai-Xizang (Tibet) plateau. Sci. China Ser. D-Earth Sci. 1996, 39, 410–417. [Google Scholar]
  71. Zha, X.J.; Niu, B.; Li, M.; Duan, C. Increasing Impact of Precipitation on Alpine-Grassland Productivity over Last Two Decades on the Tibetan Plateau. Remote Sens. 2022, 14, 3430. [Google Scholar] [CrossRef]
  72. Zhang, T.; Zhang, Y.J.; Xu, M.J.; Zhu, J.T.; Chen, N.; Jiang, Y.B.; Huang, K.; Zu, J.X.; Liu, Y.J.; Yu, G.R. Water availability is more important than temperature in driving the carbon fluxes of an alpine meadow on the Tibetan Plateau. Agric. For. Meteorol. 2018, 256, 22–31. [Google Scholar] [CrossRef]
  73. Knapp, A.K.; Carroll, C.J.W.; Denton, E.M.; La Pierre, K.J.; Collins, S.L.; Smith, M.D. Differential sensitivity to regional-scale drought in six central US grasslands. Oecologia 2015, 177, 949–957. [Google Scholar] [CrossRef] [PubMed]
  74. Ma, Z.Q.; Guo, D.L.; Xu, X.L.; Lu, M.Z.; Bardgett, R.D.; Eissenstat, D.M.; McCormack, M.L.; Hedin, L.O. Evolutionary history resolves global organization of root functional traits. Nature 2018, 555, 94–97. [Google Scholar] [CrossRef] [PubMed]
  75. Fan, Y.; Miguez-Macho, G.; Jobbagy, E.G.; Jackson, R.B.; Otero-Casal, C. Hydrologic regulation of plant rooting depth. Proc. Natl. Acad. Sci. USA 2017, 114, 10572–10577. [Google Scholar] [CrossRef] [PubMed]
  76. Li, W.; Migliavacca, M.; Forkel, M.; Denissen, J.M.C.; Reichstein, M.; Yang, H.; Duveiller, G.; Weber, U.; Orth, R. Widespread increasing vegetation sensitivity to soil moisture. Nat. Commun. 2022, 13, 3959. [Google Scholar] [CrossRef] [PubMed]
  77. Rita, A.; Camarero, J.J.; Nole, A.; Borghetti, M.; Brunetti, M.; Pergola, N.; Serio, C.; Vicente-Serrano, S.M.; Tramutoli, V.; Ripullone, F. The impact of drought spells on forests depends on site conditions: The case of 2017 summer heat wave in southern Europe. Glob. Chang. Biol. 2020, 26, 851–863. [Google Scholar] [CrossRef] [PubMed]
  78. Yao, Y.; Liu, Y.; Fu, B.; Wang, Y.; Wang, Y.; Chen, P.; Zhan, T. A warmer winter followed by a colder summer contributed to a longer recovery time in the high latitudes of Northeast China. Agric. For. Meteorol. 2022, 321, 108979. [Google Scholar] [CrossRef]
  79. De Boeck, H.J.; Bassin, S.; Verlinden, M.; Zeiter, M.; Hiltbrunner, E.J.N.P. Simulated heat waves affected alpine grassland only in combination with drought. New Phytol. 2016, 209, 531–541. [Google Scholar] [CrossRef] [PubMed]
  80. Teuling, A.J. A hot future for European droughts. Nat. Clim. Chang. 2018, 8, 364–365. [Google Scholar] [CrossRef]
  81. Zhang, Y.; Xiao, X.; Zhou, S.; Ciais, P.; McCarthy, H.; Luo, Y. Canopy and physiological controls of GPP during drought and heat wave. Geophys. Res. Lett. 2016, 43, 3325–3333. [Google Scholar] [CrossRef]
  82. Wu, J.S.; Shen, Z.X.; Zhang, X.Z. Precipitation and species composition primarily determine the diversity-productivity relationship of alpine grasslands on the Northern Tibetan Plateau. Alp. Bot. 2014, 124, 13–25. [Google Scholar] [CrossRef]
  83. Sun, J.; Du, W.P. Effects of precipitation and temperature on net primary productivity and precipitation use efficiency across China’s grasslands. Gisci. Remote Sens. 2017, 54, 881–897. [Google Scholar] [CrossRef]
  84. Wu, J.S.; Zhang, X.Z.; Shen, Z.X.; Shi, P.L.; Xu, X.L.; Li, X.J. Grazing-Exclusion Effects on Aboveground Biomass and Water-Use Efficiency of Alpine Grasslands on the Northern Tibetan Plateau. Rangel. Ecol. Manag. 2013, 66, 454–461. [Google Scholar] [CrossRef]
  85. Li, C.; Cao, Z.; Chang, J.; Zhang, Y.; Zhu, G.; Zong, N.; He, Y.; Zhang, J.; He, N. Elevational gradient affect functional fractions of soil organic carbon and aggregates stability in a Tibetan alpine meadow. Catena 2017, 156, 139–148. [Google Scholar] [CrossRef]
  86. Tashi, S.; Singh, B.; Keitel, C.; Adams, M. Soil carbon and nitrogen stocks in forests along an altitudinal gradient in the eastern Himalayas and a meta-analysis of global data. Glob. Chang. Biol. 2016, 22, 2255–2268. [Google Scholar] [CrossRef] [PubMed]
  87. Thakur, D.; Chawla, A. Functional diversity along elevational gradients in the high altitude vegetation of the western Himalaya. Biodivers. Conserv. 2019, 28, 1977–1996. [Google Scholar] [CrossRef]
  88. Anderegg, W.R.L.; Konings, A.G.; Trugman, A.T.; Yu, K.L.; Bowling, D.R.; Gabbitas, R.; Karp, D.S.; Pacala, S.; Sperry, J.S.; Sulman, B.N.; et al. Hydraulic diversity of forests regulates ecosystem resilience during drought. Nature 2018, 561, 538–541. [Google Scholar] [CrossRef] [PubMed]
  89. Liu, D.; Wang, T.; Peñuelas, J.; Piao, S. Drought resistance enhanced by tree species diversity in global forests. Nat. Geosci. 2022, 15, 800–804. [Google Scholar] [CrossRef]
  90. Ma, Q.; Su, Y.; Niu, C.; Ma, Q.; Hu, T.; Luo, X.; Tai, X.; Qiu, T.; Zhang, Y.; Bales, R.C.; et al. Tree mortality during long-term droughts is lower in structurally complex forest stands. Nat. Commun. 2023, 14, 7467. [Google Scholar] [CrossRef]
Figure 1. Spatial distributions of vegetation types, elevation, and climate variables on the alpine grasslands of the QTP: (a) alpine grassland types; (b) elevation; (c,d) average air temperature and total precipitation of growing season from 1982 to 2015.
Figure 1. Spatial distributions of vegetation types, elevation, and climate variables on the alpine grasslands of the QTP: (a) alpine grassland types; (b) elevation; (c,d) average air temperature and total precipitation of growing season from 1982 to 2015.
Sustainability 16 06697 g001
Figure 2. Spatial patterns of the maximum correlation coefficients ((a) Rmax_lag) and corresponding months ((b) lagged month) between the NDVI and the 1-month SPEI and the preceding 1 to 24 months from 1982 to 2015. Inset in (a) represents the pixels with significant correlation coefficients with p < 0.05 (red: positive; blue: negative).
Figure 2. Spatial patterns of the maximum correlation coefficients ((a) Rmax_lag) and corresponding months ((b) lagged month) between the NDVI and the 1-month SPEI and the preceding 1 to 24 months from 1982 to 2015. Inset in (a) represents the pixels with significant correlation coefficients with p < 0.05 (red: positive; blue: negative).
Sustainability 16 06697 g002
Figure 3. Spatial distributions of the maximum correlation coefficients ((a) Rmax_cml) and corresponding months ((b) accumulated month) between the 1-to-24-month SPEI and the NDVI from 1982 to 2015. Inset in (a) represents the pixels with significant correlation coefficients with p < 0.05 (red: positive; blue: negative).
Figure 3. Spatial distributions of the maximum correlation coefficients ((a) Rmax_cml) and corresponding months ((b) accumulated month) between the 1-to-24-month SPEI and the NDVI from 1982 to 2015. Inset in (a) represents the pixels with significant correlation coefficients with p < 0.05 (red: positive; blue: negative).
Sustainability 16 06697 g003
Figure 4. Relationships between timescales (lagged and accumulated months) of drought effects and elevation during the growing season. Relationships between lagged months and elevation on alpine meadows (a) and the alpine steppe (b); relationships between accumulated months and elevation in alpine meadows (c) and the alpine steppe (d). The error bars and shade denote the standard errors and 95% confidence interval of regression fitting.
Figure 4. Relationships between timescales (lagged and accumulated months) of drought effects and elevation during the growing season. Relationships between lagged months and elevation on alpine meadows (a) and the alpine steppe (b); relationships between accumulated months and elevation in alpine meadows (c) and the alpine steppe (d). The error bars and shade denote the standard errors and 95% confidence interval of regression fitting.
Sustainability 16 06697 g004
Figure 5. Relationships between timescales (lagged and accumulated months) of drought effects and temperature during the growing season. Relationships between lagged months and temperature in alpine meadows (a) and the alpine steppe (b); relationships between accumulated months and temperature in alpine meadows (c) and the alpine steppe (d). The error bars and shade denote the standard errors and 95% confidence interval of regression fitting.
Figure 5. Relationships between timescales (lagged and accumulated months) of drought effects and temperature during the growing season. Relationships between lagged months and temperature in alpine meadows (a) and the alpine steppe (b); relationships between accumulated months and temperature in alpine meadows (c) and the alpine steppe (d). The error bars and shade denote the standard errors and 95% confidence interval of regression fitting.
Sustainability 16 06697 g005
Figure 6. Relationships between timescales (lagged and accumulated months) of drought effects and precipitation during the growing season. Relationships between lagged months and precipitation in alpine meadows (a) and the alpine steppe (b); relationships between accumulated months and precipitation in alpine meadows (c) and the alpine steppe (d). The error bars and shade denote the standard errors and 95% confidence interval of regression fitting.
Figure 6. Relationships between timescales (lagged and accumulated months) of drought effects and precipitation during the growing season. Relationships between lagged months and precipitation in alpine meadows (a) and the alpine steppe (b); relationships between accumulated months and precipitation in alpine meadows (c) and the alpine steppe (d). The error bars and shade denote the standard errors and 95% confidence interval of regression fitting.
Sustainability 16 06697 g006
Figure 7. Relative contributions of elevation, temperature, and precipitation to timescales of drought effects on grassland NDVI. Relative contributions of predicting variables to lagged months in alpine meadows (a) and the alpine steppe (b). Relative contributions of prediction variables to accumulated months in alpine meadows (c) and the alpine steppe (d).
Figure 7. Relative contributions of elevation, temperature, and precipitation to timescales of drought effects on grassland NDVI. Relative contributions of predicting variables to lagged months in alpine meadows (a) and the alpine steppe (b). Relative contributions of prediction variables to accumulated months in alpine meadows (c) and the alpine steppe (d).
Sustainability 16 06697 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, X.; Hu, Z.; Zhang, Z.; Tang, J.; Niu, B. Altitude-Shifted Climate Variables Dominate the Drought Effects on Alpine Grasslands over the Qinghai–Tibetan Plateau. Sustainability 2024, 16, 6697. https://doi.org/10.3390/su16156697

AMA Style

Wang X, Hu Z, Zhang Z, Tang J, Niu B. Altitude-Shifted Climate Variables Dominate the Drought Effects on Alpine Grasslands over the Qinghai–Tibetan Plateau. Sustainability. 2024; 16(15):6697. https://doi.org/10.3390/su16156697

Chicago/Turabian Style

Wang, Xiangtao, Zhigang Hu, Ziwei Zhang, Jiwang Tang, and Ben Niu. 2024. "Altitude-Shifted Climate Variables Dominate the Drought Effects on Alpine Grasslands over the Qinghai–Tibetan Plateau" Sustainability 16, no. 15: 6697. https://doi.org/10.3390/su16156697

APA Style

Wang, X., Hu, Z., Zhang, Z., Tang, J., & Niu, B. (2024). Altitude-Shifted Climate Variables Dominate the Drought Effects on Alpine Grasslands over the Qinghai–Tibetan Plateau. Sustainability, 16(15), 6697. https://doi.org/10.3390/su16156697

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop