Spatiotemporal Variation of Vegetation on the Qinghai–Tibet Plateau and the Inﬂuence of Climatic Factors and Human Activities on Vegetation Trend (2000–2019)

: Vegetation is the terrestrial ecosystem component most sensitive to climate change. The Qinghai–Tibet Plateau (QTP), characterized by a cold climate and vulnerable ecosystems, has experienced signiﬁcant warming in previous decades. Identifying the variation in vegetation coverage and elucidating its main driving factors are critical for ecological protection on the QTP. In this study, MOD13A2 Normalized Di ﬀ erence Vegetation Index (NDVI) data in the growing season (May to September) was used to represent QTP vegetation coverage during 2000–2019. The univariate linear regression, partial correlation analysis, residual analysis, and the Hurst exponent were used to detect the vegetation spatiotemporal dynamic, analyze the relationship between the vegetation and main driving factors, and predict the future vegetation dynamic. The growing season NDVI (GNDVI) of the QTP showed an extremely signiﬁcant rate of increase (0.0011 / a) during the study period, and 79.29% of the vegetated areas showed a greening trend. Over the past 20 years, the northeast, mid-east, and western edges of the plateau have been cooling and wetting, while the southwest, mid-west, and southeast have been warming and drying. Di ﬀ erent climatic conditions lead to spatial di ﬀ erences in the response of plateau vegetation to climatic factors with generally 1–4 months lag time. The vegetation in the north of the plateau was mainly positively correlated with moisture, and negatively correlated with temperature, while the southern part showed positive correlation with temperature and negative correlation with moisture. Due to the enhancement of cooling and wetting trend in the last decade (2010–2019), especially in the south of the plateau, the greening trend of the plateau vegetation slowed down appreciably and even degraded in some areas. Human activities were mainly concentrated in the eastern part of the plateau—and its positive e ﬀ ect on vegetation was gradually increasing in most areas during study period, especially in the northeastern part. However, vegetation degradation caused by human activities in the southeast of the plateau should not be ignored. The future vegetation dynamic based on the Hurst exponent showed that the plateau faces a higher risk of vegetation degradation, which deserves more attention. This study explored the e ﬀ ect of climatic factors and human activities on vegetation of the QTP, thereby providing some guidance for the study of vegetation dynamic in the alpine areas.


Introduction
The world has experienced significant warming over the past few decades, and this warming has produced a significant impact on vegetation [1][2][3]. Vegetation plays an irreplaceable role in a negative impact on vegetation [19,[54][55][56][57]. On the other hand, some human activities also have obvious positive effects on vegetation, such as environmental and ecological construction projects, have been implemented [58,59]. Restricted by natural conditions, the human activities intensity of the QTP was low overall; however, it has shown an increasing trend, which is significantly higher than the global level in the past few decades [60]. Therefore, in the context of climate warming, the impact of continuously increasing human activities on the QTP, a fragile (but important) region, has become a hot issue of global change.
The trend of vegetation dynamics not only reflects the conditions of vegetation changes during the study period, but also indicates its possible future directions [61]. The latter is more important because it allows us to adopt specific measures to cope with predicted changes in vegetation [62]. The Hurst exponent, which was an effective method to quantitatively describe the long-range dependence of time series data and to predict possible future change trends, has been widely used in hydrology, meteorology, finance, and other fields [63]. Several previous studies have used the Hurst exponent to analyze future vegetation dynamics [61,[64][65][66]. Peng et al. [61] used the Hurst exponent method based on an AVHRR GIMMS NDVI dataset obtained during 1982-2003 to predict future vegetation change on the QTP. The result showed the Hurst exponent in most areas of the QTP was >0.5, which represented the consistency of the trends of vegetation dynamics after the study period. However, predictions of future vegetation variation based on recent trends of vegetation change on the QTP have not been reported.
As a simple and effective vegetation index, NDVI has been widely used in the study of vegetation dynamics, and numerous studies have used NDVI as a representative of vegetation activities [2,33,42,67]. This study used MOD13A2 NDVI data with 16-d temporal resolution and 1-km spatial resolution, acquired during 2000-2019, to analyze recent vegetation dynamics in the growing season on the QTP. Temperature data were obtained from MOD11A2 land surface temperature. In this paper, considering the influence of potential evapotranspiration, the moist coefficient was used instead of precipitation to more accurately reflect the relationship between water and vegetation dynamics. The moist coefficient was calculated from GPM precipitation data and MOD16A2 potential evapotranspiration data. We analyzed the climate change during the study period and calculated the correlation between vegetation and climatic factors on the QTP. On this basis, we used residual analysis to explore the impact of human activities on vegetation. Finally, based on the Hurst exponent, we predict the future vegetation trend. The purpose of this paper is (1) to understand the changing pattern of vegetation on the QTP in the past 20 years; (2) to understand the main driving mechanism of vegetation change from the perspective of both climatic and human factors.

Study Area
The QTP is located in southwestern China, which covers an area of >2.5 million km 2 . It encompasses most of Tibet and Qinghai Province, and parts of Xinjiang, Gansu, Sichuan, and Yunnan provinces ( Figure 1). The mean elevation of the QTP is more than 4000 m [7], which makes the plateau colder than any other regions at the same latitude [68]. In most of the QTP, the mean annual temperature is <0 • C, and the annual accumulated temperature >10 • C is <2000 • C [69]. The mean annual precipitation across the plateau is < 400 mm, which decreases from the southeast to the northwest with notable spatial difference [6]. Areas of forest and shrub are distributed mostly in the southeast of the plateau, while the desert is found primarily in northwestern areas of the QTP. As the two principal vegetation types, meadow and steppe occupy more than half the area of the plateau.

NDVI Data
MOD13A2 NDVI data acquired from the National Aeronautics and Space Administration (https://ladsweb.modaps.eosdis.nasa.gov/) were used in this research. This data comprised the 16day synthesized L3 data product with 1-km spatial resolution. Image splicing and projection transformation were applied using the MODIS Reprojection Tools (MRT). The Albers equal-area conic projection model was adopted with the Krasovsky datum. To further eliminate the contamination of clouds, we used a Savitzky-Golay approach to smooth the annual NDVI cycle, and then replaced the low-quality pixels based on the reliability data of MOD13A2 according to the following conditions: if the pixel value in the reliability layer is equal to 0, the original NDVI value is retained, while the pixel value is not equal to 0, then the reconstructed NDVI value is used to replace it. Through this process, the new annual NDVI data series were reconstructed [52,70]. The QTP has a short growing season, which generally concentrates from May to September; therefore, the mean annual growing season NDVI (GNDVI) value was averaged based on NDVI data from May-September each year. We calculated the 20-year mean GNDVI, and then masked out the pixels with a perennial mean of GNDVI less than 0.1 to exclude sparse and non-vegetated areas [2,10,52].

Temperature Data
The temperature data were acquired from MOD11A2 data. This data comprised the 8-day synthesized L3 data product with 1-km spatial resolution, and included daytime and nighttime landsurface temperature. The temperature data were processed using the same method as NDVI data above, and the difference is that the temperature data quality control was through the QC layer in MOD11A2. The daily temperature were obtained by averaging the daytime and nighttime landsurface temperature. The mean temperature of the growing season was obtained by averaging the daily temperature from May to September.

NDVI Data
MOD13A2 NDVI data acquired from the National Aeronautics and Space Administration (https://ladsweb.modaps.eosdis.nasa.gov/) were used in this research. This data comprised the 16-day synthesized L3 data product with 1-km spatial resolution. Image splicing and projection transformation were applied using the MODIS Reprojection Tools (MRT). The Albers equal-area conic projection model was adopted with the Krasovsky datum. To further eliminate the contamination of clouds, we used a Savitzky-Golay approach to smooth the annual NDVI cycle, and then replaced the low-quality pixels based on the reliability data of MOD13A2 according to the following conditions: if the pixel value in the reliability layer is equal to 0, the original NDVI value is retained, while the pixel value is not equal to 0, then the reconstructed NDVI value is used to replace it. Through this process, the new annual NDVI data series were reconstructed [52,70]. The QTP has a short growing season, which generally concentrates from May to September; therefore, the mean annual growing season NDVI (GNDVI) value was averaged based on NDVI data from May-September each year. We calculated the 20-year mean GNDVI, and then masked out the pixels with a perennial mean of GNDVI less than 0.1 to exclude sparse and non-vegetated areas [2,10,52].

Temperature Data
The temperature data were acquired from MOD11A2 data. This data comprised the 8-day synthesized L3 data product with 1-km spatial resolution, and included daytime and nighttime land-surface temperature. The temperature data were processed using the same method as NDVI data above, and the difference is that the temperature data quality control was through the QC layer in MOD11A2. The daily temperature were obtained by averaging the daytime and nighttime land-surface temperature. The mean temperature of the growing season was obtained by averaging the daily temperature from May to September.

Moisture Data
Considering the effect of evapotranspiration on available water, we used moist coefficient as an indicator of water availability: Moist coe f f icient = Precipitation Potential evapotranspiration (1) The greater the moist coefficient is, the better the moisture conditions for vegetation growth. The precipitation data were obtained from the Global Precipitation Measurement (GPM) data, which comprised daily precipitation data with 0.1 • spatial resolution. The potential evapotranspiration data were obtained from MOD16A2. This data is the summation of 8-day total potential water with 500 m spatial resolution. After preprocessing, both data were resampled to the same spatial resolution as NDVI data (1 km), and then counted total precipitation and the total potential evapotranspiration in the growing season (May to September), respectively, the ratio of total precipitation to total potential evapotranspiration was the moist coefficient in the growing season. We acknowledge that, when precipitation and potential evapotranspiration were resampled to the same spatial resolution as NDVI data, it was impossible to eliminate uncertainties in spatial variation because of the change of spatial resolution, especially for precipitation data with coarse resolution (0.1 • × 0.1 • ).

Linear Trend Analysis
The univariate linear regression was applied to calculate the trends in the GNDVI and meteorological data of the QTP during 2000-2019. The specific calculation formula is as follows: where n is the length of the research period (n = 20 in this study), D i represents the GNDVI and meteorological data value in the i-th year, and slope is the change rate of the regression equation.
If slope > 0, the trend of the GNDVI and meteorological data during the 20-year period was increasing; if slope < 0, the trend was decreasing. In addition, the Mann-Kendall (MK) test was used to test slope significance.

Partial Correlation Analysis
The relationships between the GNDVI, temperature, and moist coefficient were evaluated by partial correlation coefficients [67]. The partial correlation coefficients were calculated as follows: (3) where r 12 , r 13 , and r 23 are the correlation coefficients between the GNDVI and temperature, GNDVI and moist coefficient, and temperature and moist coefficient, respectively. Here, r 12•3 represents the partial correlation coefficient of the GNDVI and temperature, and r 13•2 represents the partial correlation coefficient between the GNDVI and moist coefficient.

Residual analysis
The residual analysis method was used to analyze the impact of human activities on vegetation change [71,72]. Based on 20 years of the GNDVI, temperature, and moist coefficient data, we establish the multiple linear regression model of vegetation and climatic factors through the least square method at pixel scale:ŷ whereŷ is the calculated value of the GNDVI, x 1 is temperature, x 2 is the moist coefficient, ε is the random error, a and b are regression coefficients. The calculated GND VI value was obtained through the temperature and moist coefficient of each year. Then, the annual GNDVI residual was obtained by calculating the interpolation between the calculated GNDVI value and the actual value, and then linear regression was performed with the year to obtain the GNDVI residual trend. Positive residuals represent the promotion effect of human activities on the vegetation, while negative residuals represent inhibition effects.

Hurst Exponent Method
The Hurst exponent based on the annual GNDVI values during 2000-2019 was used to predict vegetation change trends and to assess the sustainability of the vegetation time series. The most commonly used Rescaled Range Analysis (R/S) analysis method, was adopted to calculate the Hurst exponent [73]: (a) Divide the long time series x i (I = 1, 2, 3 . . . , N) into A subseries with length of n, in which A × n = N. The subseries is represented by I a (a = 1, 2, 3 . . . , A), and the element in I a could be represented by N (k, m), k = 1, 2, 3 . . . , n, and m = 1, 2, 3 . . . , A.
(b) Calculate the mean value (E a ) of I a: (c) Create the cumulative deviation time series X k,a of I a : k = 1, 2, 3..., n (d) Calculate the range R a of I a : k = 1, 2, 3..., n; (e) Create the standard deviation S a of I a : (f) Create the rescale range R a /S a (g) For A subseries with length of n, A rescale ranges can be calculated, and their mean value can be calculated to get the rescale range when the length is n: (h) Starting from n = 3, repeat the above operations for all integers n that satisfy A × n = N to obtain the sequence of the double range of [R/S] n , n = 1, 2, 3 . . . , N.
(i) Linear regression was performed with Log(R/S) n as factor and Log(n) as dependent variable: where c is a constant. H is the Hurst exponent with the range of 0-1. When 0 < H < 0.5, the future trend of the series is opposite to that of the past, and the closer the value of H to 0, the more obvious the reversal of the trend. When 0.5 < H < 1, the future trend of the series is consistent with that of the past, and the closer the value of H to 0, the more obvious the consistency of the trend. When H = 0.5, the time series information is a mutually independent random sequence for which the future change trend is uncertain. The vegetation trends in the future was obtained from superimposing the GNDVI slope and Hurst exponent as the pattern given in Table 1. In addition, during a relatively long study period, the dynamics of vegetation in different time periods are likely to show different trends. Therefore, in the process of the GNDVI dynamic analysis and correlation analysis, we divided the study period into three time periods so that we can understand the dynamics of vegetation in different periods: The whole study period (2000-2019), the first decade (2000-2009) and the second decade (2010-2019). During 2000-2019, the mean GNDVI in the whole QTP was 0.3201, showing a decreasing trend from east to west ( Figure 2). The sparse and non-vegetated areas (<0.1) occupied 30.74% of the whole plateau area, mainly including the Qaidam Basin in the north and the Kunlun Mountains in the northwest. For vegetated areas, the areal proportions of different GNDVI levels were 2.22% (>0.7), 20.36% (0.5-0.7), 22.72% (0.3-0.5), and 54.70% (0.1-0.3), respectively. During 2000-2019, the change rate of the GNDVI across the entire QTP was 0.0011/a, showing an extremely significant increasing trend (p < 0.01) ( Figure 3). The increase rate of the GNDVI in the latter decade (2010-2019) was significantly lower than that in the previous decade (2000-2009), and the increasing trend also became non-significant (p > 0.05). The GNDVI dropped significantly from 2011 to 2014, and then continued to rise, forming an obvious U-shaped curve. This difference in GNDVI trends led to the non-significant GNDVI variation in the latter decade on the QTP. During 2000-2019, the change rate of the GNDVI across the entire QTP was 0.0011/a, showing an extremely significant increasing trend (p < 0.01) ( Figure 3). The increase rate of the GNDVI in the latter decade (2010-2019) was significantly lower than that in the previous decade (2000-2009), and the increasing trend also became non-significant (p > 0.05). The GNDVI dropped significantly from 2011 to 2014, and then continued to rise, forming an obvious U-shaped curve. This difference in GNDVI trends led to the non-significant GNDVI variation in the latter decade on the QTP.  During 2000-2019, the GNDVI in 79.29% of the vegetated areas showed an increasing trend (Table 2), and 31.77% of the vegetated area reached a significant level (p < 0.05), which were mainly concentrated in the north of the QTP, especially the northeastern parts ( Figure 4). The strongest increasing trend areas (slope > 0.004) were mainly located around Qinghai Lake and the Yellow River Basin in the northeast of the plateau, and the increasing trend in the east of the plateau was stronger than that in the west. On the other hand, the GNDVI decreasing areas accounted for 20.71% of the vegetated areas, but only 1.23% of areas reached a significant level (p < 0.05). The GNDVI decreasing areas mainly concentrated in the Qumalai prefecture in the central plateau, the Nagqu prefecture in north Tibet, the alpine valleys of the Yarlung Zangbo River and the Nu River in the southern parts and the eastern edge of the plateau. During 2000-2019, the GNDVI in 79.29% of the vegetated areas showed an increasing trend (Table 2), and 31.77% of the vegetated area reached a significant level (p < 0.05), which were mainly concentrated in the north of the QTP, especially the northeastern parts ( Figure 4). The strongest increasing trend areas (slope > 0.004) were mainly located around Qinghai Lake and the Yellow River Basin in the northeast of the plateau, and the increasing trend in the east of the plateau was stronger Remote Sens. 2020, 12, 3150 9 of 22 than that in the west. On the other hand, the GNDVI decreasing areas accounted for 20.71% of the vegetated areas, but only 1.23% of areas reached a significant level (p < 0.05). The GNDVI decreasing areas mainly concentrated in the Qumalai prefecture in the central plateau, the Nagqu prefecture in north Tibet, the alpine valleys of the Yarlung Zangbo River and the Nu River in the southern parts and the eastern edge of the plateau.  During 2000-2019, the GNDVI in 79.29% of the vegetated areas showed an increasing trend (Table 2), and 31.77% of the vegetated area reached a significant level (p < 0.05), which were mainly concentrated in the north of the QTP, especially the northeastern parts ( Figure 4). The strongest increasing trend areas (slope > 0.004) were mainly located around Qinghai Lake and the Yellow River Basin in the northeast of the plateau, and the increasing trend in the east of the plateau was stronger than that in the west. On the other hand, the GNDVI decreasing areas accounted for 20.71% of the vegetated areas, but only 1.23% of areas reached a significant level (p < 0.05). The GNDVI decreasing areas mainly concentrated in the Qumalai prefecture in the central plateau, the Nagqu prefecture in north Tibet, the alpine valleys of the Yarlung Zangbo River and the Nu River in the southern parts and the eastern edge of the plateau.   Table 2), indicating that the increasing trend of the GNDVI on the QTP has been weakened in the past decade. The areas where the GNDVI trend has changed significantly included the southwestern, mid-eastern and southeastern parts of the QTP. The trend change of the southwestern part was mainly from negative to positive, while that of the mid-eastern and southeastern parts plateau were mainly from positive to negative ( Figure 4).

Vegetation Responses to Climate Change on the QTP
In Figure 5a, the temperature distribution of QTP displayed a pattern of high in northeast and southwest, low in central and southern, with the temperature in most areas (80.18%) between 5−15°C. The spatial distribution of moisture on the QTP decreased from southeast to northwest, and about half (47.27%) of the areas were in arid and semi-arid climate (Figure 5b).

Vegetation Responses to Climate Change on the QTP
In Figure 5a, the temperature distribution of QTP displayed a pattern of high in northeast and southwest, low in central and southern, with the temperature in most areas (80.18%) between 5−15 ℃. The spatial distribution of moisture on the QTP decreased from southeast to northwest, and about half (47.27%) of the areas were in arid and semi-arid climate (Figure 5b). In the past 20 years, the mean temperature in the growing season across the whole QTP showed a significant increasing trend (p < 0.05) with a rate of 0.0229 ℃/a, among which 65.65% of vegetated areas showed an increasing trend in temperature, while 34.35% showed a decreasing trend ( Table 3). The cooling areas were mainly distributed in the northeast, mid-east, and western edge of the plateau, and the most cooling area was in the south of the Qinghai Lake. In contrast, the warming areas included the mid-west, southwest, and southeast of the plateau, especially the Nagqu Prefecture in northern Tibet ( Figure 6). The moist coefficient in the growing season across the whole QTP showed a slightly increasing trend (p > 0.05) with a rate of 0.0038/a, indicating that the water supply was basically stable during the study period. The areas of increasing moist coefficient were mainly in the northeast of the plateau, accounting for 59.61% of the vegetated areas, while the decreasing areas were mainly in the southwest and southeast of the plateau, accounting for 40.39%. On the whole, western and southeastern parts of the QTP were warming and drying, while most other areas were warming and wetting during the study period. In the past 20 years, the mean temperature in the growing season across the whole QTP showed a significant increasing trend (p < 0.05) with a rate of 0.0229°C/a, among which 65.65% of vegetated areas showed an increasing trend in temperature, while 34.35% showed a decreasing trend ( Table 3). The cooling areas were mainly distributed in the northeast, mid-east, and western edge of the plateau, and the most cooling area was in the south of the Qinghai Lake. In contrast, the warming areas included the mid-west, southwest, and southeast of the plateau, especially the Nagqu Prefecture in northern Tibet ( Figure 6). The moist coefficient in the growing season across the whole QTP showed a slightly increasing trend (p > 0.05) with a rate of 0.0038/a, indicating that the water supply was basically stable during the study period. The areas of increasing moist coefficient were mainly in the northeast of the plateau, accounting for 59.61% of the vegetated areas, while the decreasing areas were mainly in the southwest and southeast of the plateau, accounting for 40.39%. On the whole, western and southeastern parts of the QTP were warming and drying, while most other areas were warming and wetting during the study period.   (Table 3). These changes indicated that the plateau has undergone a transition from a warming and drying trend to a cooling and wetting trend as a whole, and this transition in the climate change trend was more obvious in the western and southern part of the plateau.  Figure 7 showed the difference in response of vegetation to temperature and moisture under different environmental conditions. In cold areas (mean temperature < 5 °C), the temperature had a high positive correlation with the GNDVI, while the moisture showed a significantly negative correlation with GNVDI. With the warming of the environment, the correlation coefficient between the temperature and GNDVI decreased, while the correlation coefficient between water and the GNDVI gradually increased. Under relatively warm conditions (mean temperature > 10 °C), the effect of moisture on vegetation exceeded that of temperature. With the further increase of temperature, the correlation coefficient between the temperature and GNDVI becomes negative, which means that However, climate change trends on the QTP were not consistent in different decades. The temperature change rate was 0.0609°C/a in the first decade (2000-2009), and then dropped to −0.0176°C/a in the latter decade (2010-2019), while the moist coefficient change rate was −0.0397/a in the first decade, and then rose to 0.0155/a in the latter decade. The areal proportion of the warming area also decreased, while that of the wetting area increased in different periods (Table 3). These changes indicated that the plateau has undergone a transition from a warming and drying trend to a cooling and wetting trend as a whole, and this transition in the climate change trend was more obvious in the western and southern part of the plateau. Figure 7 showed the difference in response of vegetation to temperature and moisture under different environmental conditions. In cold areas (mean temperature < 5 • C), the temperature had a high positive correlation with the GNDVI, while the moisture showed a significantly negative correlation with GNVDI. With the warming of the environment, the correlation coefficient between the temperature and GNDVI decreased, while the correlation coefficient between water and the GNDVI gradually increased. Under relatively warm conditions (mean temperature > 10 • C), the effect of moisture on vegetation exceeded that of temperature. With the further increase of temperature, the correlation coefficient between the temperature and GNDVI becomes negative, which means that in the high-temperature areas, the warming will limit the vegetation growth. The result showed that low temperature was the main factor limiting vegetation growth in the low-temperature zone, and the increase of water had a negative inhibitory effect on vegetation growth. With the continuous rise of the ambient temperature, the roles of temperature and water are gradually exchanged. Finally, when the temperature reached a certain degree, the negative effect of drought brought by further warming would inhibit vegetation growth, and water would become the main factor affecting vegetation growth. The correlation coefficient between the moisture and GNDVI also experienced similar changes to that of temperature with the change of different moisture conditions. The above results showed that the variation of environmental conditions changed the limiting factors of vegetation growth and made the responses of vegetation to climate factors change.
when the temperature reached a certain degree, the negative effect of drought brought by further warming would inhibit vegetation growth, and water would become the main factor affecting vegetation growth. The correlation coefficient between the moisture and GNDVI also experienced similar changes to that of temperature with the change of different moisture conditions. The above results showed that the variation of environmental conditions changed the limiting factors of vegetation growth and made the responses of vegetation to climate factors change. Due to the great spatial heterogeneity of the climatic conditions on the QTP, the response of vegetation to climatic factors also had a strong spatial difference. During 2000-2019, the partial correlation coefficient between the GNDVI and temperature was negative for 46.41% pixels (8.16% significantly with p < 0.05), and positive for 53.59% pixels (11.60% significantly with p < 0.05) ( Table  4). The positive correlation pixels were mainly distributed in the northeast and southwest of the plateau, while the negative correlation pixels were mainly distributed in the central and southeast of the plateau (Figure 8). The spatial pattern of the partial correlation coefficient between the GNDVI and moist coefficient showed a significant north-south difference. The positive correlation pixels were mainly distributed in the north of the plateau, accounting for 58.33% of total pixels (13.77% significantly with p < 0.05), while that negative correlation pixels were mainly distributed in the south, accounting for 41.67% of total pixels (4.95% significantly with p < 0.05). In different decades, the spatial pattern of the partial correlation coefficient between temperature and moist coefficient was basically consistent, except for some parts of the North Tibet Plateau, and the areal proportion of different partial correlation levels between the GNDVI and temperature and moist coefficient also did not change significantly over different decades (Table 4). Both phenomena suggested that the response mechanism of vegetation to climate change has not changed significantly in most areas of the QTP. In the North Tibet Plateau, the positive correlation between the temperature and GNDVI decreased, while the positive correlation between the moisture and GNDVI increased, which meant that available water gradually replaces temperature as the main factor affecting vegetation growth. This transformation might be associated with the lack of water, causing by the continuous climate warming and drying in this region. Due to the great spatial heterogeneity of the climatic conditions on the QTP, the response of vegetation to climatic factors also had a strong spatial difference. During 2000-2019, the partial correlation coefficient between the GNDVI and temperature was negative for 46.41% pixels (8.16% significantly with p < 0.05), and positive for 53.59% pixels (11.60% significantly with p < 0.05) ( Table 4). The positive correlation pixels were mainly distributed in the northeast and southwest of the plateau, while the negative correlation pixels were mainly distributed in the central and southeast of the plateau (Figure 8). The spatial pattern of the partial correlation coefficient between the GNDVI and moist coefficient showed a significant north-south difference. The positive correlation pixels were mainly distributed in the north of the plateau, accounting for 58.33% of total pixels (13.77% significantly with p < 0.05), while that negative correlation pixels were mainly distributed in the south, accounting for 41.67% of total pixels (4.95% significantly with p < 0.05). In different decades, the spatial pattern of the partial correlation coefficient between temperature and moist coefficient was basically consistent, except for some parts of the North Tibet Plateau, and the areal proportion of different partial correlation levels between the GNDVI and temperature and moist coefficient also did not change significantly over different decades (Table 4). Both phenomena suggested that the response mechanism of vegetation to climate change has not changed significantly in most areas of the QTP. In the North Tibet Plateau, the positive correlation between the temperature and GNDVI decreased, while the positive correlation between the moisture and GNDVI increased, which meant that available water gradually replaces temperature as the main factor affecting vegetation growth. This transformation might be associated with the lack of water, causing by the continuous climate warming and drying in this region.   Previous studies have highlighted that the relationship between the NDVI and various climatic factors has a general time lag [2,42,64], the shorter the lag time is, the more immediate the response of vegetation to climate change is. We studied the lag response between the GNDVI and climatic factors at an interval of eight days, and took the date corresponding to the maximum partial correlation coefficient as the hysteresis time. In order to prevent the negative correlation coefficient from being concealed, we converted the partial correlation coefficient to the absolute value in advance. The lag time between the GNDVI and temperature was 1-month for 36.45% pixels, 2-month for 20.93% pixels, 3-month for 19.17% pixels, and 4-month for 23.45% pixels, while the lag time for moist coefficient was 1-month for 24.87% pixels, 2-month for 19.60% pixels, 3-month for 25.83% pixels, and 4-month for 29.69% pixels, respectively. This result indicated that the hysteretic response of vegetation to climatic factors is common on the QTP. In terms of spatial distribution, the vegetation with a short lag response to temperature was mainly in the southeast, northwest, and parts of the northeast, while the vegetation with long lag response was mainly in the central and southwest of the plateau (Figure 9a). The vegetation with a short lag response to moisture was mainly in the northwest and parts of the southeast, while the vegetation with long lag response was mainly in the east and west of the plateau. Previous studies have highlighted that the relationship between the NDVI and various climatic factors has a general time lag [2,42,64], the shorter the lag time is, the more immediate the response of vegetation to climate change is. We studied the lag response between the GNDVI and climatic factors at an interval of eight days, and took the date corresponding to the maximum partial correlation coefficient as the hysteresis time. In order to prevent the negative correlation coefficient from being concealed, we converted the partial correlation coefficient to the absolute value in advance. The lag time between the GNDVI and temperature was 1-month for 36.45% pixels, 2-month for 20.93% pixels, 3-month for 19.17% pixels, and 4-month for 23.45% pixels, while the lag time for moist coefficient was 1-month for 24.87% pixels, 2-month for 19.60% pixels, 3-month for 25.83% pixels, and 4-month for 29.69% pixels, respectively. This result indicated that the hysteretic response of vegetation to climatic factors is common on the QTP. In terms of spatial distribution, the vegetation with a short lag response to temperature was mainly in the southeast, northwest, and parts of the northeast, while the vegetation with long lag response was mainly in the central and southwest of the plateau (Figure 9a). The vegetation with a short lag response to moisture was mainly in the northwest and parts of the southeast, while the vegetation with long lag response was mainly in the east and west of the plateau.

The Influence of Human Factors on Vegetation Change on the QTP
In the whole QTP, the GNDVI residuals were generally negative in the first decade (2000-2009), indicating that human activities had a negative effect on vegetation change at this stage (Figure 10a). Since 2010, the GNDVI residual turned positive, indicating that human activities began to play a promoting role in vegetation change. During the whole study period (2000-2019), the GNDVI showed an extremely significant upward trend with a rate of 0.0007/a (p < 0.001), indicating that human activities had an increasingly significant positive impact on vegetation on the QTP. In terms of the spatial pattern, the impact of human activities was low in the western plateau because of the sparsely populated, so the variation trend of the GNDVI residual was weak (Figure 10b). In contrast,

The Influence of Human Factors on Vegetation Change on the QTP
In the whole QTP, the GNDVI residuals were generally negative in the first decade (2000-2009), indicating that human activities had a negative effect on vegetation change at this stage (Figure 10a). Since 2010, the GNDVI residual turned positive, indicating that human activities began to play a promoting role in vegetation change. During the whole study period (2000-2019), the GNDVI showed an extremely significant upward trend with a rate of 0.0007/a (p < 0.001), indicating that human activities had an increasingly significant positive impact on vegetation on the QTP. In terms of the spatial pattern, the impact of human activities was low in the western plateau because of the sparsely populated, so the variation trend of the GNDVI residual was weak (Figure 10b). In contrast, the impact of human activities on vegetation in the eastern part of the plateau was relatively strong, especially in the northeastern part of the plateau, where the vast majority of the plateau's population was located. The GNDVI residual in the northeast and mid-east of the plateau showed a relatively high increasing trend, represented that the human activities in these areas have strong positive promotion effects on vegetation. The negative trend of the GNDVI residual was mainly distributed in three regions, including the three-river source region in the south of Qinghai Province, the Lhasa region in the southwest of the plateau, and the southeastern edge region of the plateau.
correlation coefficient between the GNDVI and (a) temperature and (b) moist coefficient on the QTP during 2000-2019.

The Influence of Human Factors on Vegetation Change on the QTP
In the whole QTP, the GNDVI residuals were generally negative in the first decade (2000-2009), indicating that human activities had a negative effect on vegetation change at this stage (Figure 10a). Since 2010, the GNDVI residual turned positive, indicating that human activities began to play a promoting role in vegetation change. During the whole study period (2000-2019), the GNDVI showed an extremely significant upward trend with a rate of 0.0007/a (p < 0.001), indicating that human activities had an increasingly significant positive impact on vegetation on the QTP. In terms of the spatial pattern, the impact of human activities was low in the western plateau because of the sparsely populated, so the variation trend of the GNDVI residual was weak (Figure 10b). In contrast, the impact of human activities on vegetation in the eastern part of the plateau was relatively strong, especially in the northeastern part of the plateau, where the vast majority of the plateau's population was located. The GNDVI residual in the northeast and mid-east of the plateau showed a relatively high increasing trend, represented that the human activities in these areas have strong positive promotion effects on vegetation. The negative trend of the GNDVI residual was mainly distributed in three regions, including the three-river source region in the south of Qinghai Province, the Lhasa region in the southwest of the plateau, and the southeastern edge region of the plateau.

Prediction of Future Trends of Vegetation Trends
The Hurst exponent (H) of the GNDVI of the whole QTP was in the range 0.06-0.96 (mean value: 0.44). Pixels with H > 0.5 accounted for 25.21% of the total pixels, while pixels with H < 0.5 accounted for 74.79% (Figure 11a), indicating that the reversal of the current GNDVI trends would occur in most areas of the plateau in the future, and this reverse trend was probably related to the slowing of the rising trend of GNDVI in the second decade of the study period. However, the H values of 60.09% pixels were between 0.4 and 0.6, meaning that their future GNDVI trends are more random than a

Prediction of Future Trends of Vegetation Trends
The Hurst exponent (H) of the GNDVI of the whole QTP was in the range 0.06-0.96 (mean value: 0.44). Pixels with H > 0.5 accounted for 25.21% of the total pixels, while pixels with H < 0.5 accounted for 74.79% (Figure 11a), indicating that the reversal of the current GNDVI trends would occur in most areas of the plateau in the future, and this reverse trend was probably related to the slowing of the rising trend of GNDVI in the second decade of the study period. However, the H values of 60.09% pixels were between 0.4 and 0.6, meaning that their future GNDVI trends are more random than a definite direction. The pixels with strong, consistent trends (H > 0.6) and strong reverse trends (H < 0.4) accounted for 5.14% and 34.75% of the total pixels, respectively.
Remote Sens. 2018, 10, x FOR PEER REVIEW 14 of 21 definite direction. The pixels with strong, consistent trends (H > 0.6) and strong reverse trends (H < 0.4) accounted for 5.14% and 34.75% of the total pixels, respectively. Based on the Hurst exponent and current GNDVI trend, we predicted the spatial pattern of the GNDVI trend in the future on the QTP (Figure 11b). For most parts of the vegetated areas, the GNDVI future trends are not clear in future, as the H value ranges from 0.4 to 0.6. 17.69% of the plateau vegetated areas will remain stable because of the low current GNDVI change rate, and these areas are distributed primarily in the western parts of the plateau. The areas with a decreasing trend, which account for 18.43% of the vegetated area, are largely scattered in the eastern parts of the plateau. The areas with increasing trend only account for 3.82%, which are mainly distributed in the central and southern parts of the plateau. In terms of the whole plateau, the future vegetation degradation area is larger than the improvement area, which represents a certain risk of vegetation degradation on the QTP in the future, and this risk is most obvious in the northeast of the plateau.

Discussion
Over the past 20 years, the whole QTP showed an increasing trend in GNDVI, which was consistent with the trend identified in an earlier period   [61], indicating the whole QTP has undergone a process of continuous vegetation greening since 1982. Due to the spatial heterogeneity of climate change, the GNDVI trends of the plateau showed significant spatial difference [74,75].
Low temperature is one of the main factors limiting the growth of high-altitude vegetation [52], and climate warming can significantly promote vegetation growth on the QTP [2,76]. Our results showed that the positive response to temperature appeared in the southeast and northwest, while that negative mainly appeared in the southeast and northwest, which was similar to a previous study [77]. In the southeastern and northwestern parts, which were dominated by cold weather (Figure 5a), the increased temperature would directly promote photosynthesis of vegetation, and break the restriction of low temperature on vegetation growth [78,79]. The negative correlation with temperature appears in the northeast and southwest of the plateau, indicating that the high temperature would inhibit vegetation growth, and the likely cause is increased respiration caused by the higher temperature [80]. In addition, the increased temperature could also thin the active layer of the permafrost [33] and increase the soil evapotranspiration [36,81], both could reduce the amount of water available in the topsoil, and further lead to loss of soil organic matter and release of soil carbon and nitrogen, which in turn would limit vegetation growth [82,83].
The northern and southwestern parts of the QTP were dominated by the arid and semi-arid climate ( Figure 5b) and we found that vegetation in these areas is positively correlated with moisture Based on the Hurst exponent and current GNDVI trend, we predicted the spatial pattern of the GNDVI trend in the future on the QTP (Figure 11b). For most parts of the vegetated areas, the GNDVI future trends are not clear in future, as the H value ranges from 0.4 to 0.6. 17.69% of the plateau vegetated areas will remain stable because of the low current GNDVI change rate, and these areas are distributed primarily in the western parts of the plateau. The areas with a decreasing trend, which account for 18.43% of the vegetated area, are largely scattered in the eastern parts of the plateau. The areas with increasing trend only account for 3.82%, which are mainly distributed in the central and southern parts of the plateau. In terms of the whole plateau, the future vegetation degradation area is larger than the improvement area, which represents a certain risk of vegetation degradation on the QTP in the future, and this risk is most obvious in the northeast of the plateau.

Discussion
Over the past 20 years, the whole QTP showed an increasing trend in GNDVI, which was consistent with the trend identified in an earlier period   [61], indicating the whole QTP has undergone a process of continuous vegetation greening since 1982. Due to the spatial heterogeneity of climate change, the GNDVI trends of the plateau showed significant spatial difference [74,75].
Low temperature is one of the main factors limiting the growth of high-altitude vegetation [52], and climate warming can significantly promote vegetation growth on the QTP [2,76]. Our results showed that the positive response to temperature appeared in the southeast and northwest, while that negative mainly appeared in the southeast and northwest, which was similar to a previous study [77]. In the southeastern and northwestern parts, which were dominated by cold weather (Figure 5a), the increased temperature would directly promote photosynthesis of vegetation, and break the restriction of low temperature on vegetation growth [78,79]. The negative correlation with temperature appears in the northeast and southwest of the plateau, indicating that the high temperature would inhibit vegetation growth, and the likely cause is increased respiration caused by the higher temperature [80]. In addition, the increased temperature could also thin the active layer of the permafrost [33] and increase the soil evapotranspiration [36,81], both could reduce the amount of water available in the topsoil, and further lead to loss of soil organic matter and release of soil carbon and nitrogen, which in turn would limit vegetation growth [82,83].
The northern and southwestern parts of the QTP were dominated by the arid and semi-arid climate ( Figure 5b) and we found that vegetation in these areas is positively correlated with moisture ( Figure 8), which was consistent with previous studies [2,36]. The relatively dry environment limits the supply of water to vegetation growth, making the vegetation in these areas very sensitive to available water, especially under the background of continuous climate warming [84,85]. In contrast, the climate in the mid-east and southeast of the plateau is relatively humid with abundant water resources to offset the drought in soils caused by rising temperatures; increased moisture could reduce the plant photosynthesis and temperature, due to decreased solar radiation [42]. In addition, the increasing precipitation would restrict vegetation growth by producing soil erosion [43] and temperature decreasing [44], so the vegetation in these areas showed a negative correlation with moisture, and the influence of temperature on vegetation becomes more positive [37].
Although the GNDVI of the QTP showed an upward trend throughout the whole study period, the increasing trend of GNDVI slowed markedly in the last decade (Figure 3), and the most obvious change areas were in the central and south of the plateau (Figure 4). Compared the climate trends in different decades (Figure 6), we found that the causes of the GNDVI trend decline are different. In the central area, a warming and drying climate caused a decrease in available water, which led to a downward trend in GNDVI, while in the south, especially in the southeast areas, the downward trend in GNDVI was mainly caused by a cooling and wetting climate. Interestingly, the same cooling and wetting climate change in the southwest of the plateau has led to an upward trend in GNDVI, which may be due to the difference in vegetation types between the two areas. These changes in GNDVI and climatic factors well explains the different mechanisms of vegetation response to climate change in different areas of the QTP.
The vegetation in the QTP has shown hysteresis response to climate change generally with a lag of 1-4 months, which was similar to the findings of Piao's study on the national scale [42]. The lag time is most likely due to the growth of vegetation that need time to accumulate the soil moisture, heat, and nutrient [86]. In addition, the low-temperature environment on the QTP made the rate of organic nitrogen and organic phosphate mineralization of soil slow, and the soil lacked available nitrogen and organic phosphorus that could be absorbed and utilized by vegetation, which might lead to the lag effect of vegetation growth in summer on spring precipitation [87]. Our results showed that the higher the correlation coefficient between vegetation and a certain climate factor is, the shorter the hysteresis response time is, which reflected the timely response of vegetation to the change of growth restriction factors [64,88].
The unique geographical conditions made the most of the QTP not suitable for human habitation, especially in the west, and most human activities were concentrated in the eastern part of the plateau [54,60,89]; the GNVDI residual analysis in this paper showed the similar result. The overall residual trend of the plateau was mainly positive, especially in the northeast of the plateau, indicating that a series of ecological projects which have been implemented has produced a positive influence on the vegetation of the QTP [36,43,90,91]. However, the GNDVI residual trend remained negative in some areas of the plateau, including the mid-eastern, southeastern and southwestern areas around Lhasa, where urban expansion and infrastructure construction might be the main factors affecting vegetation change [54,60]. Therefore, appropriate measures should be taken to prevent the further degradation of vegetation in these areas under human influence.
The Hurst exponent based on the GNDVI data (2000-2019) of the entire QTP is 0.44, indicating a very random future change of vegetation. This randomness was likely to be the significant decline in GNDVI between 2011 and 2014, which interrupted the increasing trend of the GNDVI in the latter decade, further undermining the consistency of the increasing trend of the GNDVI throughout the whole study period. Hurst exponent is a calculation based on the time series GNDVI value and does not involve other factors; but since the GNDVI is affected by climatic factors, climate change also indirectly affects Hurst exponent. The change of climate trends in different decades made the GNDVI trend more fluctuant, which further affected the Hurst exponent. Therefore, the relationship between Hurst exponent and climate change deserves further study. In addition, although Hurst exponent can be used to predict the trend of future vegetation dynamic, it cannot provide a clear duration, which makes vegetation prediction more like a possible risk indication. Thus, determining the likely duration is a more important issue that needs further study in the future [61].
Considering the important role of heat and water in vegetation growth, we used land surface temperature and moisture coefficient to quantify the two main environmental factors, and analyzed the relationship between them and vegetation dynamics. However, due to the large spatial heterogeneity on the QTP, other factors may also play important roles affecting vegetation growth, such as solar radiation [2,52], winter snow [37], and permafrost soil [75,92,93]. The combination of many factors may lead to errors in our assessment of the relationship between vegetation and major climatic factors. Therefore, the influence of other factors on vegetation on the QTP should be considered comprehensively in future studies.

Conclusions
The GNDVI of whole QTP showed a significant increasing trend in the past 20-years with a rate of 0.0011/a, with the exception of parts of the southwest, 79.29% of the plateau displayed greening trends. During the whole study period, both temperature and moist coefficient of the whole plateau dispalyed an upward trend, with the rising rates of 0.0229 • C/a and 0.0038/a, respectively. Spatially, there was a strong spatial heterogeneity of climate trends on the plateau: The climate trends in the mid-east, northeast, and west edge of the plateau were mainly cooling and wetting, while that in the mid-west, southwest, and southeast of the plateau mainly showed warming and drying trend. Because of the spatially heterogeneous of climatic conditions, the response of vegetation on the plateau to climatic factors has also formed obvious spatial differences: The northeast and southwest vegetation was negatively correlated with temperature and positively correlated with moisture. In contrast, the southeast vegetation was negatively correlated with moisture and positively correlated with temperature; while the northwest vegetation was positively correlated with both climatic factors. In the second decade of the study period (2010-2019), the cooling and wetting climate trend strengthened, which weakened the greening trend of the plateau, especially in the southern parts. The GNDVI residual of the whole plateau showed an extremely significant upward trend with the rate of 0.0007/a in the past 20 years, indicating the positive effect of human activities on vegetation, especially in the northeastern part of the plateau. However, the destruction of vegetation by human beings in the southeast of the plateau should not be ignored and more attention should be given to prevent further disturbance of human activities. The future vegetation dynamic based on the Hurst exponent showed that the area of degradation is likely to become larger than that of improvement, indicating a higher risk of vegetation degradation. To reduce the risk of ecological degradation on the plateau, we need to take a series of measures to protect the plateau's ecological environment: Further progressing the ecological protection projects have been carried out; establishing natural ecological protection regions; enhancing local people's awareness of ecological and environmental protection; carrying out afforestation and the construction of artificial grassland and prohibiting excessive development. In addition, rational resource planning and allocation will help us better cope with the ecological security challenges on the QTP.