Spatial and Temporal Differences in Alpine Meadow, Alpine Steppe and All Vegetation of the Qinghai-Tibetan Plateau and Their Responses to Climate Change

: Alpine meadow and alpine steppe are the two most widely distributed nonzonal vegetation types in the Qinghai-Tibet Plateau. In the context of global climate change, the differences in spatial-temporal variation trends and their responses to climate change are discussed. It is of great signiﬁcance to reveal the response of the Qinghai-Tibet Plateau to global climate change and the construction of ecological security barriers. This study takes alpine meadow, alpine steppe and the overall vegetation of the Qinghai-Tibet Plateau as the research objects. The normalized difference vegetation index (NDVI) data and meteorological data were used as the data sources between 2000 and 2018. By using the mean value method, threshold method, trend analysis method and correlation analysis method, the spatial and temporal variation trends in the alpine meadow, alpine steppe and the overall vegetation of the Qinghai-Tibet Plateau were compared and analyzed, and their differences in the responses to climate change were discussed. The results showed the following: (1) The growing season length of alpine meadow was 145~289 d, while that of alpine steppe and the overall vegetation of the Qinghai-Tibet Plateau was 161~273 d, and their growing season lengths were signiﬁcantly shorter than that of alpine meadow. (2) The annual variation trends of the growing season NDVI for the alpine meadow, alpine steppe and the overall vegetation of the Qinghai-Tibet Plateau increased obviously, but their ﬂuctuation range and change rate were signiﬁcantly different. (3) The overall vegetation improvement in the Qinghai-Tibet Plateau was primarily dominated by alpine steppe and alpine meadow, while the degradation was primarily dominated by alpine meadow. (4) The responses between the growing season NDVI and climatic factors in the alpine meadow, alpine steppe and the overall vegetation of the Qinghai-Tibet Plateau had great spatial heterogeneity in the Qinghai-Tibet Plateau. These ﬁndings provide evidence towards understanding the characteristics of the different vegetation types in the Qinghai-Tibet Plateau and their spatial differences in response to climate change.


Introduction
As an important part of the global terrestrial ecosystem, vegetation plays an important role in the global material and energy cycles, carbon balance regulation and maintenance of climate stability [1][2][3][4]. Additionally, vegetation is the natural link between soil, atmosphere and water, and is both affected by and positively responsive to climate change [5][6][7][8]. In recent years, in the context of global climate change, vegetation change has attracted more attention. Vegetation dynamics affect global environmental evolution, and the study of surface vegetation cover has become a popular field of global change research. The vegetation-climate response is of great practical significance to studying the relationship in vegetation variation and the change response on the QTP, and this information is of great significance to the sustainable development of grassland livestock husbandry and the improved function of grassland ecosystem services.
Based on the above research status, in this study, the NDVI data were synthesized every 16 d from MOD 13A2 from 2000 to 2018, and the observation data from the ground meteorological stations were selected as the data sources, the differences in growing seasons, the annual and seasonal variation trends of the NDVI in the AM, AS and the overall vegetation of the Qinghai-Tibetan Plateau (abbreviated as VQTP, NDVI > 0.1) were compared and analyzed, and the correlation between the growing season NDVI (AM, AS and VQTP) and climate factors (precipitation and temperature) were discussed from the spatial scale pixel by pixel. The results reveal the response of mixed vegetation and different alpine grassland types to climate change on the QTP to provide a scientific basis and decision-making services for monitoring the stability of the regional ecosystem, providing disaster warning, mitigating vegetation degradation and taking corresponding protection measures within the context of climate change.

Study Area
Located in Southwest China, the QTP is the highest physical geographic unit in the world, with an average altitude of over 4000 m. It is also an important ecological security barrier for China and Asia [59,60]. It ranges from 26 • 00 12"N~39 • 46 50"N and 73 • 18 52" E~104 • 46 59" E ( Figure 1). The south-north extent of the study area reaches from the southern foot of the Himalayas to the northern Kunlun and Qilian Mountains, and the eastwest extent reaches from the Hengduan Mountains to the Pamirs Plateau. The total area is approximately 2.57 × 10 6 km 2 , accounting for 26.8% of China's total land area [61,62]. The administrative divisions cover six provinces, including the Xizang (Tibet) Autonomous Region and Qinghai Province, as well as parts of Xinjiang Autonomous Region and Yunnan, Sichuan and Gansu Provinces ( Figure 1). Due to the fragile ecosystem of the QTP and its vulnerability to climate change, the behavior of the plateau vegetation system can predict global changes in a more timely and clear fashion than other regions, making it an ideal region to study global change [63].
ote Sens. 2021, 13, x FOR PEER REVIEW 3 of 24 heterogeneity in vegetation variation and the change response on the QTP, and this information is of great significance to the sustainable development of grassland livestock husbandry and the improved function of grassland ecosystem services. Based on the above research status, in this study, the NDVI data were synthesized every 16 d from MOD 13A2 from 2000 to 2018, and the observation data from the ground meteorological stations were selected as the data sources, the differences in growing seasons, the annual and seasonal variation trends of the NDVI in the AM, AS and the overall vegetation of the Qinghai-Tibetan Plateau (abbreviated as VQTP, NDVI > 0.1) were compared and analyzed, and the correlation between the growing season NDVI (AM, AS and VQTP) and climate factors (precipitation and temperature) were discussed from the spatial scale pixel by pixel. The results reveal the response of mixed vegetation and different alpine grassland types to climate change on the QTP to provide a scientific basis and decision-making services for monitoring the stability of the regional ecosystem, providing disaster warning, mitigating vegetation degradation and taking corresponding protection measures within the context of climate change.

Study Area
Located in Southwest China, the QTP is the highest physical geographic unit in the world, with an average altitude of over 4000 m. It is also an important ecological security barrier for China and Asia [59,60]. It ranges from 26°00′12″N ~ 39°46′50″N and 73°18′52″ E ~ 104°46′59" E ( Figure 1). The south-north extent of the study area reaches from the southern foot of the Himalayas to the northern Kunlun and Qilian Mountains, and the east-west extent reaches from the Hengduan Mountains to the Pamirs Plateau. The total area is approximately 2.57 × 10 6 km 2 , accounting for 26.8% of China's total land area [61,62]. The administrative divisions cover six provinces, including the Xizang (Tibet) Autonomous Region and Qinghai Province, as well as parts of Xinjiang Autonomous Region and Yunnan, Sichuan and Gansu Provinces ( Figure 1). Due to the fragile ecosystem of the QTP and its vulnerability to climate change, the behavior of the plateau vegetation system can predict global changes in a more timely and clear fashion than other regions, making it an ideal region to study global change [63]. Due to its special terrain and geographical position, the QTP has formed a unique nonzonal plateau climate [64]. Its climate is characterized by strong solar radiation, low Due to its special terrain and geographical position, the QTP has formed a unique nonzonal plateau climate [64]. Its climate is characterized by strong solar radiation, low air temperatures throughout the year, large regional differences and diurnal temperature Remote Sens. 2021, 13, 669 4 of 24 differences, distinct dry and wet conditions [65], less precipitation and an uneven distribution, with more in the southeast and less in the northwest [48]. Under the influence of hydrothermal conditions and altitude, the climate of the QTP changes from humid and warm in the southeast to cold and dry in the northwest. The hydrothermal, dry and wet conditions of the QTP indirectly affect the growth process of alpine grasslands [66]. There are various vegetation types in the QTP, including mountain forest, mountain shrub, AM, AS and alpine desert from the southeast to northwest [67], among which grassland resources are the most widely distributed. Therefore, the study of vegetation change in alpine grasslands (i.e., AM and AS) on the QTP is of great value for the rational use of grassland resources and the protection of alpine ecosystems.

Data Sources and Preprocessing
The data used in this study mainly include remote sensing data set (NDVI), meteorological data, vegetation type data and Digital Elevation Model data of the study area, and their detailed introduction and preprocessing are as follows:

RS Data
The vegetation index used in this study was from the Land Processes Distributed Active Archive Center (LP DAAC) of the United States (http://lp-daac.usgs.gov/main.asp). The MOD13A2 Version 6 product, which was extracted from the Terra Moderate Resolution Imaging Spectroradiometer (MODIS), has a temporal resolution of 16 d, a spatial resolution of 1 km, and a time span of 2000-2018. The track numbers of MODIS products in the study area were h23v05, h24v05, h25v04, h25v05, h25v06, h26v05 and h26v06. The data format is HDF-EOS, and the projection is sinusoidal.
The MODIS Re-Projection Tool (MRT) was used for image mosaic, format conversion, resampling and projection transformation. The original HDF format was converted into Geotiff, and the sinusoidal projection was converted into the Albers conical equal area with the datum D_WGS_1984. The resampling method was the nearest neighbor method with a resolution of 1 km. The maximum value composite (MVC) method was then used to obtain the monthly NDVI from 2000 to 2018, which effectively removed the influence of the clouds, atmosphere and solar elevation angle [68]. Lastly, the mean value method was used to obtain the NDVI during the growing season and different seasons to eliminate the influence of climatic anomalies in extreme years on vegetation growth. The vector boundary of the QTP was used to clip the NDVI data in ArcGIS to obtain the NDVI time-series data set of the study area. The value range of the NDVI is −1.0~1.0. According to a large number of research results in the QTP [13,69,70], it is generally believed that an average NDVI value in the growing season above 0.1 indicates vegetation cover, while a value below 0.1 indicates no vegetation cover on the surface, such as bare soil, desert, the Gobi, a water body, snow, ice and clouds. Because the object of this study was vegetation, to calculate the results more accurately, only pixels with an average NDVI value during the growing season of greater than 0.1 were calculated in this study. The vegetation cover area after the mask was employed as the targeted region of this paper.

Meteorological Data
The meteorological data used in this study were provided by the Resource and Environment Science and Data Center (RESDC) and the China Meteorological Data Service Center (CMDC). The spatial interpolation data sets of annual precipitation and annual average temperature from 2000 to 2015 were derived from the RESDC, with a spatial resolution of 1 km. The spatial interpolation data sets of annual precipitation and annual average temperatures from 2016 to 2018 were obtained by using ANUSPLIN software based on the monitoring data from meteorological stations, and the spatial resolution of the data sets was the same as that of the MODIS-NDVI remote sensing image. In this study, a total of 180 meteorological stations in the QTP and in its surrounding provinces were selected for spatial interpolation ( Figure 1). The data sources, interpolation methods and spatial resolutions used during both periods were all the same, to establish the spatial interpolation data set of temperature and precipitation on the QTP from 2000 to 2018.

Vegetation and Digital Elevation Model Data
The vegetation type data were obtained from a digitized 1:1,000,000 vegetation map of China compiled by the Chinese Academy of Sciences (CAS), which was provided by the RESDC (http://www.resdc.cn). The vector boundary of the QTP was used to clip out the vegetation type data from the QTP, and the spatial distribution range of AM and AS on the QTP was classified by summarizing and extracting the same vegetation type (Figure 2) [71]. study, a total of 180 meteorological stations in the QTP and in its surrounding provinces were selected for spatial interpolation ( Figure 1). The data sources, interpolation methods and spatial resolutions used during both periods were all the same, to establish the spatial interpolation data set of temperature and precipitation on the QTP from 2000 to 2018.

Vegetation and Digital Elevation Model Data
The vegetation type data were obtained from a digitized 1:1,000,000 vegetation map of China compiled by the Chinese Academy of Sciences (CAS), which was provided by the RESDC (http://www.resdc.cn). The vector boundary of the QTP was used to clip out the vegetation type data from the QTP, and the spatial distribution range of AM and AS on the QTP was classified by summarizing and extracting the same vegetation type (Figure 2) [71].
The Vegetation and Digital Elevation Model (DEM) data were taken from radar terrain mapping SRTM (Shuttle Radar Topography Mission, SRTM) data from the U.S. space shuttle Endeavour. The data set was based on the latest SRTM V4.1 data and downloaded from the RESDC with a spatial resolution of 1 km, which was primarily used as the covariate of meteorological data interpolation and topographic map production.

Determination of the Growing Season
In this paper, three phenological parameters, namely, the start of the growing season (SGS), the end of the growing season (EOS) and the length of the growing season (LGS), were selected as indicators for monitoring the vegetation growth season on the QTP to determine the spatial-temporal differences in the SGS, EOS and LGS of AM, AS and VQTP. Based on the MODIS-NDVI time-series data set, the maximum slope method was used to determine the start and end dates of the vegetation growing season. According to the vegetation growth process, the maximum slope method defines the time corresponding to the point with the maximum rate of change in the NDVI fitting curve as the The Vegetation and Digital Elevation Model (DEM) data were taken from radar terrain mapping SRTM (Shuttle Radar Topography Mission, SRTM) data from the U.S. space shuttle Endeavour. The data set was based on the latest SRTM V4.1 data and downloaded from the RESDC with a spatial resolution of 1 km, which was primarily used as the covariate of meteorological data interpolation and topographic map production.

Determination of the Growing Season
In this paper, three phenological parameters, namely, the start of the growing season (SGS), the end of the growing season (EOS) and the length of the growing season (LGS), were selected as indicators for monitoring the vegetation growth season on the QTP to determine the spatial-temporal differences in the SGS, EOS and LGS of AM, AS and VQTP. Based on the MODIS-NDVI time-series data set, the maximum slope method was used to determine the start and end dates of the vegetation growing season. According to the vegetation growth process, the maximum slope method defines the time corresponding to the point with the maximum rate of change in the NDVI fitting curve as the corresponding Remote Sens. 2021, 13, 669 6 of 24 key phenological period for the vegetation [2,72]. The slope of the fitting curve is calculated as follows: where NDVI slope(t) is the slope of the NDVI fitting curve at time t, and NDVI (t+1) and NDVI (t) are the NDVI values at times t + 1 and t, respectively. The date corresponding to the maximum change rate of the NDVI slope(t) on the left side of the growth curve (i.e., the inflection point when rising) is defined as the beginning of the growing season, and the date corresponding to the maximum rate of change on the right side of the growth curve (i.e., the inflection point when falling) is defined as the end of the growing season. The length between the beginning and end of the growing season is defined as the growing season length.

Trend Analysis Method
The dynamic variation trend of the NDVI in each pixel was simulated by using the linear regression method, and the spatial and temporal pattern evolution of vegetation during the growing season was comprehensively reflected by the pixel scale. Slope _NDVI is defined as the slope of the NDVI for the annual growing season as fitted using the least square method within a given time range, which can better reflect the variation characteristics of aboveground biomass and vegetation cover, thus revealing the spatial pattern and variation trend in vegetation on the QTP over the last 20 years. The Slope _NDVI is calculated as follows [73,74]: where Slope _NDVI is the trend of the growing season NDVI, n is the range of the time series (n = 19), i is the order of years from 1 to n, and NDVI i is the growing season NDVI for the ith year. Slope _NDVI > 0 indicates that the growing season NDVI shows an increasing trend during the study period and vice versa.

Correlational Analysis Method
The correlation analysis method was used to study the relationship between vegetation and climatic factors. The correlation coefficient between the annual growing season NDVI and the annual precipitation and annual mean temperature of the corresponding year on the QTP from 2000 to 2018 was calculated, and a significance test was conducted to reflect the degree of correlation between different climatic factors and the NDVI. Its formula is as follows [75,76]: where R is the correlation coefficient between the variables x and y, which represents the correlation degree between the two factors, and its value ranges from −1 to 1. n is the cumulative number of years during the monitoring period. x i and y i are the values of independent variables x and y in the ith year. x and y are the means of the values of independent variables x and y for all years. According to the correlation coefficient test table, the correlation degree is divided into five grades: extremely significant negative correlation (R < 0, p < 0.01), significantly negative correlation (R < 0, 0.01 < p < 0.05), no significant correlation (p > 0.05), significantly positive correlation (R > 0, 0.01 < p < 0.05), and extremely significant positive correlation (R > 0, p < 0.01).

Mann-Kendall (MK) Test
The MK test, which was developed by Mann and Kendall [77], is a nonparametric test method used for detecting statistically significant trends in the NDVI and the continuous environmental variable time series [78]. It does not require the samples to be independent Remote Sens. 2021, 13, 669 7 of 24 and normally distributed and has been widely used in vegetation dynamics [79,80]. The MK test statistic Z is computed as follows: in which n is the length of the time series, x j and x i are the sequential data values in time series i and j (j > i), and sgn(x j − x i ) is the symbolic function as follows: The variance is computed as follows: in which Var(S) represents the variance of S, m is the number of tied groups and t i is the number of ties for extent i. A tied group is a set of sample data that have the same value.
Positive Z values indicate increasing trends, while negative Z values show decreasing trends in the time series. The significance test is estimated at the given α significance level. When |Z| > Z 1−α/2 , the null hypothesis is rejected, and a significant trend is present in the time series. In this study, significance levels of α = 0.05 and α = 0.01 were used, and their Z 1−α/2 values were 1.96 and 2.576, respectively.

Phenological Characteristic Differences of AM, AS and QTP
The three phenological parameters of AM, AS and VQTP, namely, the SGS, EGS and LGS, were extracted from the NDVI time-series curves for 16 d from 2000 to 2018 using the maximum slope method. As shown in Figure 3, the NDVI time-series curves of AM, AS and VQTP all show parabolic shapes that are low on both sides and high in the middle, reflecting the changing rules of each NDVI curve over different periods of a year. Through the comparative analysis of various phenological parameters, it was found that the SGS of the AM occurred in late May (145 d), and the SGS of the AS and VQTP was in early June (161 d). The SGS of the AM occurred approximately half a month earlier than that of AM and VQTP. The EGS of the AM was in the middle of October (289 d), and the EGS of the AS and VQTP was in late September (273 d), half a month earlier than that of the AM, which was consistent with the field observation results. The SGS of the AM lasted for 145 days, from 145 d to 289 d, and the SGS of the AS and VQTP started from 161 d to 273 d and lasted for 112 days. The LGSs of the AS and VQTP were significantly shorter than that of the AM. It can be observed from the variation curves of the NDVI in the AM, AS and VQTP that the NDVI of VQTP was greater than that of AM in the non-growing season (that was before the SGS and after the EGS) ( Figure 3), while the NDVI of AM was greater than that of VQTP in the growing season. This indicated that the NDVI of AM was lower as a whole because it was dormant in the non-growing season. However, due to the influence of forests and evergreen vegetation in the southeastern part of the QTP, the overall NDVI of VQTP during this period was relatively high. When the growing season came, AM began to grow and gradually reached the growth peak period, while the NDVI Remote Sens. 2021, 13, 669 8 of 24 of VQTP was slightly lower than that of AM in the growing season due to the influence of large-area low vegetation cover regions such as AS and desert grassland. The maximum NDVI values of the AM and VQTP appeared at the end of July (209 d), while the maximum value of the AS appeared in the middle of August (225 d), which was slightly delayed. The primary reasons for the differences in phenological characteristics were as follows: the AM was mainly distributed over the central and eastern parts of the plateau, while the AS was mainly distributed over the northwestern part of the plateau. There was a significant regional difference in the surface temperature. The precipitation was largely controlled by the southwest monsoon and gradually decreased from the southeast to northwest. The hydrothermal conditions in the distribution area of the AM were better than those in the distribution area of the AS. AM began to grow and gradually reached the growth peak period, while the NDVI of VQTP was slightly lower than that of AM in the growing season due to the influence of large-area low vegetation cover regions such as AS and desert grassland. The maximum NDVI values of the AM and VQTP appeared at the end of July (209 d), while the maximum value of the AS appeared in the middle of August (225 d), which was slightly delayed. The primary reasons for the differences in phenological characteristics were as follows: the AM was mainly distributed over the central and eastern parts of the plateau, while the AS was mainly distributed over the northwestern part of the plateau. There was a significant regional difference in the surface temperature. The precipitation was largely controlled by the southwest monsoon and gradually decreased from the southeast to northwest. The hydrothermal conditions in the distribution area of the AM were better than those in the distribution area of the AS. LGS, length of the growing season; AM, alpine meadow; AS, alpine steppe; VQTP, overall vegetation of the Qinghai-Tibetan Plateau.

Interannual Variation Characteristics of the NDVI during the Growing Seasons of AM, AS and VQTP
The mean value method was used to obtain the mean value of the NDVI during the growing seasons of the AM, AS and VQTP, and then the annual variation trend of the three was analyzed. As shown in the change curves in Figure 4

Interannual Variation Characteristics of the NDVI during the Growing Seasons of AM, AS and VQTP
The mean value method was used to obtain the mean value of the NDVI during the growing seasons of the AM, AS and VQTP, and then the annual variation trend of the three was analyzed. As shown in the change curves in Figure 4 in the growing season, AS is different from them in all periods. The NDVI of the three types has obvious differences in the value range, fluctuation amplitude and change rate. and VQTP in the growing season, AS is different from them in all periods. The NDVI of the three types has obvious differences in the value range, fluctuation amplitude and change rate.

Interannual Variation Characteristics of the AM, AS and VQTP during Different Seasons
Through the statistics on the annual average NDVI values for the AM, AS and VQTP during the growing season, spring, summer and autumn, the annual NDVI variation curves were obtained for different seasons from 2000 to 2018 ( Figure 5). Due to the influence of snow cover, the winter vegetation was not considered. The average seasonal NDVI of the AM, AS and VQTP showed increasing trends from 2000 to 2018, but the variation trends were different to some extent. Among them, the NDVI of the AM and VQTP increased at the fastest rate in the spring, with rates of 0.0014 yr −1 and 0.0015 yr −1 , respectively. The NDVI of the AS increased at the fastest rate, at 0.0011 yr -1 , during the growing season. The VQTP and AM curves primarily occurred between 0.30 and 0.40 during the growing season, and the AS curve was between 0.15 and 0.20. In the spring, the VQTP and AM curves primarily ranged between 0.15 and 0.20, the value of the VQTP curve was slightly larger than that of AM, and the AS curve primarily fluctuated at approximately 0.10. The NDVI values of all the types during the spring were significantly lower than those of the growing season. In the summer, the AM curve primarily fluctuated at approximately 0.40, which was significantly higher than that during the growing season and in the spring, while the VQTP and AS curves were slightly higher than that during the growing season. In the autumn, the AM, AS and VQTP curves were all significantly higher than those in the spring but lower than those for the summer and the growing season. Therefore, from a seasonal point of view, the AM, AS and VQTP flourishing periods all occurred in the summer, followed by autumn and lastly spring. The mean value of the growing

Interannual Variation Characteristics of the AM, AS and VQTP during Different Seasons
Through the statistics on the annual average NDVI values for the AM, AS and VQTP during the growing season, spring, summer and autumn, the annual NDVI variation curves were obtained for different seasons from 2000 to 2018 ( Figure 5). Due to the influence of snow cover, the winter vegetation was not considered. The average seasonal NDVI of the AM, AS and VQTP showed increasing trends from 2000 to 2018, but the variation trends were different to some extent. Among them, the NDVI of the AM and VQTP increased at the fastest rate in the spring, with rates of 0.0014 yr −1 and 0.0015 yr −1 , respectively. The NDVI of the AS increased at the fastest rate, at 0.0011 yr -1 , during the growing season. The VQTP and AM curves primarily occurred between 0.30 and 0.40 during the growing season, and the AS curve was between 0.15 and 0.20. In the spring, the VQTP and AM curves primarily ranged between 0.15 and 0.20, the value of the VQTP curve was slightly larger than that of AM, and the AS curve primarily fluctuated at approximately 0.10. The NDVI values of all the types during the spring were significantly lower than those of the growing season. In the summer, the AM curve primarily fluctuated at approximately 0.40, which was significantly higher than that during the growing season and in the spring, while the VQTP and AS curves were slightly higher than that during the growing season. In the autumn, the AM, AS and VQTP curves were all significantly higher than those in the spring but lower than those for the summer and the growing season. Therefore, from a seasonal point of view, the AM, AS and VQTP flourishing periods all occurred in the summer, followed by autumn and lastly spring. The mean value of the growing season was used as a reference, and the mean value of the NDVI was only less than that of the summer.
highly consistent, and the two were significantly correlated at the 0.01 level (Table 1). This result indicated that the AM made a high contribution to the overall vegetation change on the QTP. The fluctuation trends of the AS and VQTP during different seasons were also relatively consistent. The correlation between AS and VQTP was slightly weaker in the growing season (p < 0.05) but not significant in summer (p > 0.05) ( Table 1). The variation trends in the AM and AS were highly consistent in the growing season, spring and autumn (p < 0.01) but different in the summer (p > 0.05). The difference of AS in summer may be due to the overall influence of the QTP by other vegetation types or the sensitive response of the AS to climate change. In general, AM and AS are the two primary vegetation types on the QTP, and their changes play key and decisive roles in the overall vegetation change trend for the QTP.   According to the interannual variation curve of the NDVI during each season ( Figure 5), the interannual fluctuation trends of the AM and VQTP during different seasons were highly consistent, and the two were significantly correlated at the 0.01 level (Table 1). This result indicated that the AM made a high contribution to the overall vegetation change on the QTP. The fluctuation trends of the AS and VQTP during different seasons were also relatively consistent. The correlation between AS and VQTP was slightly weaker in the growing season (p < 0.05) but not significant in summer (p > 0.05) ( Table 1). The variation trends in the AM and AS were highly consistent in the growing season, spring and autumn (p < 0.01) but different in the summer (p > 0.05). The difference of AS in summer may be due to the overall influence of the QTP by other vegetation types or the sensitive response of the AS to climate change. In general, AM and AS are the two primary vegetation types on the QTP, and their changes play key and decisive roles in the overall vegetation change trend for the QTP.

Spatiotemporal Change Trends in Growing Season NDVI (AM, AS and VQTP) and Significance Test
Based on the trend analysis method, the variation trend in the NDVI during the growing season was estimated from 2000 to 2018 in the study area, and the corresponding significance level was tested by the MK method. According to the trend analysis and significance tests over the past 19 years, the spatial distribution of the slope of the NDVI during the growing season was divided into six categories: extremely significant increase (p < 0.01 and slope > 0), significant increase (0.01 < p < 0.05 and slope > 0), nonsignificant increase (0.05 < p and slope > 0), nonsignificant decrease (0.05 < p and slope < 0), significant decrease (0.01 < p < 0.05 and slope < 0), and extremely significant decrease (p < 0.01 and slope < 0) ( Figure 6). As shown in Figure 6 and Table 2 Table 2). The significantly increased areas were largely distributed throughout the north-central and northeastern parts of the QTP, especially in the Qilian Mountains, and this result was closely related to the ecological and environmental protection project and policy implementation in the Qilian Mountain region in recent years. The area showing a decreasing trend accounted for 30.56% of the total monitoring area (Figure 6a, Table 2), of which only 1.90% passed the significance test ( Figure 6b, Table 2). The areas with significant declines were mainly distributed within the surrounding areas of Namtso and Siling Co in Tibet and the source region of the Yangtze River in Qinghai. These two regions have clearly been disturbed by human activities and were the primary distribution areas of vegetation degradation throughout the QTP.
Using the AM and AS vector boundaries to extract the spatial and temporal variation trends of the growing season NDVI and its test results for the QTP, the spatial distribution of the variation trends and significance test results for the AM and AS were obtained, respectively (Figure 7). The results showed that the increasing region of the AM was mainly distributed over the central and eastern parts of the QTP, and its area accounted for 29.25% of the vegetation increase in the monitored area, of which the region with a significant increase (p < 0.05) accounted for 22.63% ( Figure 7a 1 and a 2 Table 3). In addition to the mountainous region in the southeastern QTP, the vegetation degradation trend in the AM distribution region was clear. The region with a decreasing trend in the AM accounted for 51.49% of the vegetation degradation in the monitored area, of which the significantly decreased area (p < 0.05) accounted for 39.78% of the monitored area with a significant decrease (p < 0.05), which was mainly distributed within the source region of the Yangtze River, northwest of Lhasa and along National Highway 109 (within Lhasa) (Figure 7a 1 and a 2 , Table 3). The above analysis showed that more than half of the vegetation degradation on the QTP occurred in the AM, and the significant degradation region accounted for a higher proportion of the significant degradation in the monitored area. Using the AM and AS vector boundaries to extract the spatial and temporal variation trends of the growing season NDVI and its test results for the QTP, the spatial distribution of the variation trends and significance test results for the AM and AS were obtained, respectively (Figure 7). The results showed that the increasing region of the AM was mainly distributed over the central and eastern parts of the QTP, and its area accounted for 29.25% of the vegetation increase in the monitored area, of which the region with a significant increase (p < 0.05) accounted for 22.63% (Figure 7a1 and a2, Table 3). In addition to the mountainous region in the southeastern QTP, the vegetation degradation trend in Figure 6. Spatial distribution of the NDVI slope (a) and significance level (b) in the Qinghai-Tibet Plateau (QTP) during the growing season from 2000 to 2018. p(decreasing) and p(increasing) are the p values of the decrease and increase for the growing season NDVI, respectively, which are divided into three levels: p < 0.01, 0.01 < p < 0.05, and p > 0.05. the AM distribution region was clear. The region with a decreasing trend in the AM accounted for 51.49% of the vegetation degradation in the monitored area, of which the significantly decreased area (p < 0.05) accounted for 39.78% of the monitored area with a significant decrease (p < 0.05), which was mainly distributed within the source region of the Yangtze River, northwest of Lhasa and along National Highway 109 (within Lhasa) (Figure 7a1 and a2, Table 3). The above analysis showed that more than half of the vegetation degradation on the QTP occurred in the AM, and the significant degradation region accounted for a higher proportion of the significant degradation in the monitored area. The results of the AS change trend showed (Figure 7b1 and b2, Table 3) that the region where the growing season NDVI of AS showed an increasing trend was much larger than the region showing a decreasing trend. The increasing region of the AS was largely distributed over the central and western parts of the QTP, the Qinghai Lake Basin and Qilian Mountains in the northeast, and its area accounted for 33.17% of the vegetation increase in the monitored area, among which the significantly increased region (p < 0.05) accounted for 45.39% of the monitored area that had a significant increase (p < 0.05). The region where the AS showed a downward trend accounted for 18.22% of the vegetation decline in the monitored area, and the region with a significant decline (p < 0.05) accounted for 26.05% of the monitored area that had a significant decrease (p < 0.05). The overall improvement in vegetation on the QTP was primarily based on alpine meadow and alpine steppe data. They were also the main distribution regions of vegetation degradation on the QTP. Although they accounted for a relatively high proportion, the absolute area was relatively small, with a value of only 2.35 × 10 4 km 2 . Figure 7. Spatial distribution of the AM_NDVI slope (a 1 ) and its significance level (SL) (a 2 ) and the AS_NDVI slope (b 1 ) and its significance level (SL) (b 2 ) during the growing season from 2000 to 2018. p(decreasing) and p(increasing) are the p values of the decrease and increase for the growing season NDVI, respectively, which are divided into three levels: p < 0.01, 0.01 < p < 0.05, and p > 0.05. The results of the AS change trend showed ( Figure 7b 1 and b 2 , Table 3) that the region where the growing season NDVI of AS showed an increasing trend was much larger than the region showing a decreasing trend. The increasing region of the AS was largely distributed over the central and western parts of the QTP, the Qinghai Lake Basin and Qilian Mountains in the northeast, and its area accounted for 33.17% of the vegetation increase in the monitored area, among which the significantly increased region (p < 0.05) accounted for 45.39% of the monitored area that had a significant increase (p < 0.05). The region where the AS showed a downward trend accounted for 18.22% of the vegetation decline in the monitored area, and the region with a significant decline (p < 0.05) accounted for 26.05% of the monitored area that had a significant decrease (p < 0.05). The overall improvement in vegetation on the QTP was primarily based on alpine meadow and alpine steppe data. They were also the main distribution regions of vegetation degradation on the QTP. Although they accounted for a relatively high proportion, the absolute area was relatively small, with a value of only 2.35 × 10 4 km 2 .

Interannual Variation Characteristics of Climatic Factors
The average value of the meteorological interpolation data was calculated in the study area from 2000 to 2018, and the spatial distribution of the annual precipitation and annual mean temperature was obtained for the past 19 years (Figure 8a,b). As shown in Figure 8a, the average multiyear precipitation on the QTP showed an obvious step-like distribution, gradually decreasing from the southeast to northwest, with a range of 50.81~2450.21 mm, which represented a large span. The regions with high precipitation were primarily concentrated in the western Sichuan Plateau and the northwest mountainous Yunnan Province, and the average annual precipitation was over 1500 mm. The regions with low precipitation were largely distributed in the Qaidam Basin and the central and western parts of the Ali region of Tibet. The average annual precipitation was below 300 mm. As shown in Figure 8b, the spatial distribution range of multi-year mean temperature is from −19.50 • C to 23.99 • C. Except for the higher temperature in the Qaidam Basin, the distribution of the regional average annual temperatures was similar to that of the precipitation, which gradually decreased from the southeast to northwest as the altitude increased. Due to the effect of the altitude and terrain, the regions with high temperatures were mainly distributed in the southern and southeastern valleys of Tibet at lower altitudes. The regions with low temperatures were mostly distributed in the Kunlun Mountains and the Qiangtang Plateau region. The average altitude of this region was above 4000 m, and it was the main distribution area of AS and alpine desert.
perature rise trend in the AS was extremely significant, which passed the significance test of 0.01. Combined with the spatial variation trend in the annual mean temperature (Figure  8f), except for some regions in the north and south QTP, the annual mean temperature in most other regions showed an increasing trend, among which the most obvious increasing trend occurred in the Ali region, Nagqu region and Zhongba County in the central and western regions of the plateau. This result indicates that under the influence of global warming, the warming effect is relatively obvious in areas where the AS is widely distributed.

Correlation between Growing Season NDVI (AM, AS and VQTP) and Precipitation
The spatial distribution of the correlation between the growing season NDVI and the climatic factors was determined, the significance level of the correlation coefficient between them was tested, and the correlation coefficient was graded according to the significance level (Figure 9). The spatial distribution of the significance level in the correlation coefficient between the growing season NDVI and annual precipitation (Figure 9a) shows that the regions with a positive correlation between the growing season NDVI and the Through a statistical analysis of the spatial interpolation data set on the annual mean temperature and annual precipitation on the QTP, the interannual variation in the annual mean temperature and annual precipitation was obtained for the study area for the last 19 years (Figure 8c,d). As shown in Figure 8c, the annual precipitation in the AM, AS and VQTP distribution regions showed fluctuating and decreasing trends from 2000 to 2018, and the variation trends of the three were similar. Among them, the change trends of the three were relatively stable from 2000 to 2014, showing slight increasing trends. Due to the impact of extreme drought years in 2015, the annual precipitation declined sharply after 2014, reached the lowest point in 2015, and gradually increased thereafter. The annual precipitation of the AM, AS and VQTP experienced a process of slow increase-sharp decline-gradual increase. Viewed over the entire study period, the annual precipitation of the VQTP and AS showed significant downward trends, and both passed the 0.05 significance test (Figure 8c). Viewed from the spatial variation trend in annual precipitation (Figure 8e), the annual precipitation in most regions of the QTP showed a decreasing trend, among which the decline was more obvious in the southeastern and southwestern parts of the plateau; only the southeastern part of Qinghai Province and the northwestern area of Sichuan Province increased, in addition to parts of southern Tibet. The annual variation trend in the average temperature (Figure 8d) indicates that the annual average temperatures for the AM, AS and VQTP all showed fluctuating upward trends, with rates of increase of 0.26 • C/10 yr −1 , 0.46 • C/10 yr −1 and 0.17 • C/10 yr −1 , respectively, of which the warming rate of the AS was the fastest, followed by the AM. Between them, the temperature change curves of the QTP and AM were extremely close, with a fluctuation range between −1.0 • C and 1.0 • C. The average annual temperature of the AS was below −1 • C, and the change curve was between −3.0 • C and −1.0 • C. The annual mean temperature rising trends in the QTP and AM were not obvious, while the annual mean temperature rise trend in the AS was extremely significant, which passed the significance test of 0.01. Combined with the spatial variation trend in the annual mean temperature (Figure 8f), except for some regions in the north and south QTP, the annual mean temperature in most other regions showed an increasing trend, among which the most obvious increasing trend occurred in the Ali region, Nagqu region and Zhongba County in the central and western regions of the plateau. This result indicates that under the influence of global warming, the warming effect is relatively obvious in areas where the AS is widely distributed.

Correlation between Growing Season NDVI (AM, AS and VQTP) and Precipitation
The spatial distribution of the correlation between the growing season NDVI and the climatic factors was determined, the significance level of the correlation coefficient between them was tested, and the correlation coefficient was graded according to the significance level ( Figure 9). The spatial distribution of the significance level in the correlation coefficient between the growing season NDVI and annual precipitation (Figure 9a) shows that the regions with a positive correlation between the growing season NDVI and the annual precipitation were mainly distributed in the central region of the QTP, along the northeast to the southwest direction, accounting for 51.11% of the total monitored area. The significant positive correlation areas were mainly distributed within the source region of the Yangtze and Yellow Rivers in southeastern Qinghai Province and in the central part of the Tibet Autonomous Region, accounting for 12.11% of the total monitored area, indicating that the positive feedback effect of precipitation on vegetation changes in these areas was relatively obvious. The regions with a negative correlation between the growing season NDVI and annual precipitation were mainly distributed within the western Sichuan Plateau, southeastern Tibet, northwestern Tibet and the Qilian Mountains, accounting for 48.89% of the total monitored area, among which the regions with a significant negative correlation accounted for 5.37%. Based on the analysis of the vegetation change trend, it is clear that except for the increase in annual precipitation in the western Sichuan Plateau, the annual precipitation in other negatively correlated regions decreased. However, the growing season NDVI in these regions tended to increase significantly during the growing season, especially in the Qilian Mountains. This result shows that the improvement in vegetation in these areas was not significantly affected by the precipitation. In general, the regions with a significant positive correlation between the growing season NDVI and annual precipitation on the QTP were significantly larger than the regions with a significant negative correlation.
3.5.3. Correlation between Growing Season NDVI (AM, AS and VQTP) and Temperature The correlation coefficient between the growing season NDVI and the annual mean temperature was calculated to describe the closeness between the two quantitatively, and Figure 9. The correlation and spatial distribution of the significance level between the growing season NDVI and annual precipitation, annual mean temperature: (a) VQTP_NDVI and precipitation, (a 1 ) AM_NDVI and precipitation, (a 2 ) AS_NDVI and precipitation, (b) VQTP_NDVI and temperature, (b 1 ) AM_NDVI and temperature, (b 2 ) AS_NDVI and temperature.
To understand the differences in the AM and AS responses to precipitation change on the QTP, the spatial distribution characteristics of the correlation between the AM, AS and annual precipitation were discussed in this study (Figure 9a 1 ,a 2 ). The results showed that the regions with a significant positive correlation between the growing season NDVI and the annual precipitation of the AM were mainly distributed between the 400 mm and 600 mm isohyets along the "northeast to southwest direction," accounting for 46.02% of the regions with a significant positive correlation in the study area. The significant negative correlation area accounted for 23.88% of the regions with a significant positive correlation in the study area. The regions with a significant positive correlation between the growing season NDVI and the annual precipitation of the AS were largely distributed around the 400 mm isohyet in central and eastern Qinghai Province and central Tibet, accounting for 35.85% of the regions with a significant positive correlation in the study area. The significant negative correlation area accounted for 28.32% of the regions with a significant negative correlation in the study area. Given the above analysis, it was concluded that the significant positive feedback region of AM to annual precipitation was mainly in the isohyet region of 400~600 mm, while the positive feedback region of AS to annual precipitation was mainly distributed around the 400 mm isohyet. The total proportion of regions with significant positive correlations between the AM and AS and annual precipitation in the study area was as high as 81.87%, indicating that the regions with significant positive correlations between the growing season NDVI and annual precipitation on the QTP were primarily AM and AS.

Correlation between Growing Season NDVI (AM, AS and VQTP) and Temperature
The correlation coefficient between the growing season NDVI and the annual mean temperature was calculated to describe the closeness between the two quantitatively, and the spatial distribution of the correlation coefficient was classified according to the correlation between the growing season NDVI and the annual mean temperature and its significance level (Figure 9b). The spatial distribution map showing the significance level of the correlation coefficient indicates that the growing season NDVI and annual mean temperature in most regions of the QTP were positively correlated, and its proportion accounted for 59.06% of the total monitored area, of which the significant positive correlation area accounted for 5.73% of the total monitored area, which was mainly distributed in the Ali region of northwest Tibet and in southeastern Qinghai Province. The NDVI in these regions increased with the increasing temperature. The increase in temperature promoted the growth of vegetation in this area, especially in the Ali region of Tibet. The regions with a negative correlation between the growing season NDVI and the annual mean temperature were also widely distributed. The regions with a negative correlation accounted for 40.94% of the total monitored area, and the significantly negatively correlated regions accounted for only 3.68% of the total monitored area, primarily in the western Qilian Mountains and the eastern margin of the Qaidam Basin. As indicated by the temperature variation trend in Figure 8f, the temperature in this region showed a downward trend, while the vegetation displayed an obvious increasing trend, indicating that the vegetation increased with the decreasing temperature. In other words, once the temperature increases in this region, vegetation degradation will occur, which is the total opposite to the situation in the Ali region.
To understand the differences in the responses of the AM and AS to temperature changes on the QTP, the spatial distribution characteristics of the correlation between the AM and AS and the annual mean temperature were discussed (Figure 9b 1 ,b 2 ). The results showed that the regions with a significant positive correlation between the growing season NDVI and the annual mean temperature of the AM were mainly distributed in the source region of the Yellow River in the southeastern Qinghai Province, between the −5 • C and 1 • C isotherms, of which the significant positive correlation area accounted for 28.79% of the total significant positive correlation in the study area. The regions with a significant positive correlation between the AM and AS and the annual mean temperature were all located between the −5 • C and 1 • C isotherms, accounting for 78.62% of the regions with a significant positive correlation and 58.97% of the regions with a significant negative correlation in the study area, indicating that the significant influence of the temperature on the vegetation of the QTP was mainly dominated by alpine grassland.

Discussion
At present, it is very common to use the vegetation index to monitor the dynamic changes in vegetation. The existing research includes many related studies using the annual average NDVI, annual maximum NDVI, summer NDVI and growing season NDVI. Among them, the growing season NDVI has been widely recognized and has achieved better results. However, there is no strict theoretical basis for the definition of the growing season in different vegetation monitoring studies, most of which use April to October (according to empirical values) or March to November (according to the season) as the growing season. The selection of the growing season is arbitrary and unscientific, which is likely to cause deviations in the monitoring results on vegetation changes in the same area based on the same data source. Additionally, due to the influence of comprehensive factors such as the topography, climate, soil and vegetation types, the changes in different vegetation types in different regions, the start and end of the growing season and its length are different. If the average value of the regional vegetation index is used to reflect the vegetation change trend of the entire study area and its response to climate change, the inherent details of regional vegetation change and their sensitivity differences to climate change are ignored. Based on the above problems, this study, from the perspective of remote sensing phenology and considering the differences in vegetation types, first determined the start, end and length of the growing seasons of the AM, AS and VQTP. Using the VQTP monitoring results as a reference, the phenological characteristics of the AM, AS and VQTP, the difference in the NDVI in different seasons, the spatial-temporal trends of the growing season NDVI and their correlations with climatic factors were compared. Compared with traditional vegetation change monitoring methods, this study defined the vegetation growing season more strictly and accurately and classified different types of vegetation in a targeted way, which more intuitively and comprehensively reflected the spatial and temporal differentiation characteristics of regional vegetation and its spatial heterogeneity in response to climate change.
Several studies [81][82][83] have shown that, whether in the 1980s or after 2000, the growing season NDVI of the QTP showed an overall increasing trend, and grassland dominated, with some regions being degraded, showing an "overall improvement and local degradation" pattern [84], which was consistent with the results of this study. The climate change trend showed that before 2014, the annual precipitation and annual mean temperature of the QTP both increased, and the regional climate was "warming and humid." However, after 2014, an extreme drought year, the annual precipitation in 2015 dropped rapidly, and this result has also been confirmed by relevant research [85]. Moreover, this study found that the annual precipitation on the QTP began to decline slowly from 2000 to 2018, while the annual mean temperature displayed an upward trend. The results of this study on vegetation change and climate change trends are consistent with previous studies. In comparing the spatial correlation between the growing season NDVI and climate factors, it was observed that the response degrees and spatial distributions of the growing season NDVI to climate change were different from the results of previous studies. This difference was primarily due to the data sources, research periods, definitions of the growing season, and selection of the number of weather stations. All these factors may cause a given deviation in the spatial distribution of the NDVI data, meteorological data and their correlation coefficients.
There are significant differences in the spatial distribution of vegetation NDVI and climatic factors of the QTP. It is not scientific to use only the statistical data of the research area or some site data to represent the vegetation change trend of the entire region and its response to climate change. Therefore, this study used spatial interpolation technology to interpolate temperature and precipitation into raster data with the same spatial resolution as that of the NDVI and to analyze the NDVI change trend and its correlation with meteorological factors at the pixel scale. This approach works around the deficiency in traditional results in which the surface is represented by the results on the point so that the spatial heterogeneity of the monitoring results can be displayed more intuitively in space. However, due to the vast territory of the QTP, there are relatively few ground meteorological stations, most of which are mainly located in the central and eastern regions, and the number of stations in the northwestern region is small. This characteristic will inevitably affect the accuracy of the spatial interpolation of meteorological data. In addition, due to the influence of data access permissions, this study analyzed only the interannual response relationship between the growing season NDVI and climatic factors on the temporal scale and did not consider the time lag in climate factor impacts on the growth season NDVI at different periods, such as months and seasons. Therefore, more appropriate interpolation methods will be developed in studies to follow, while focusing on improving the interpolation accuracy of meteorological factors in regions without meteorological stations and obtaining more comprehensive meteorological observation data to reveal the impact of climate factors on vegetation change in the QTP.
The changes in the NDVI of different vegetation types have different responses to climate change, and the spatial-temporal evolution of the vegetation NDVI is the combined result of natural and human activities such as climate change, soil and water conservation, and land use. This paper discusses only the influences of the two primary meteorological factors, temperature and precipitation, on vegetation change in the QTP, without considering other climatic factors, such as sunshine, evaporation, humidity and CO 2 . Additionally, since it is difficult to quantify the relevant indicators of human activities in space and because we cannot match and establish a corresponding relationship with the NDVI on the spatial scale, it was not considered in this paper. In future research, the quantitative exploration of human factors should be strengthened.

Conclusions
Based on the MOD 13A2 NDVI time-series data and the observational data from surface meteorological stations on the Qinghai-Tibet Plateau from 2000 to 2018, this study extracted the related phenological parameters from the AM, AS and VQTP and determined the phenological characteristics and differences in the three. The trend analysis method was used to analyze the spatiotemporal dynamic change trend of the growing season NDVI in the AM, AS and VQTP from the perspectives of time and space. Lastly, a correlation analysis was used to discuss the correlation between the growing season NDVI and climatic factors in the AM, AS and VQTP and their response differences. The primary conclusions are as follows: (1) Due to different vegetation types and regional differences in temperature and precipitation, the SGS, EGS and LGS of the AM, AS and VQTP were different. The SGS of the AS and VQTP were approximately half a month later than that of the AM, while the EGS were approximately half a month earlier than that of the AM. Their LGS values were clearly shorter than that of AM. (2) The growing season NDVI of the AM, AS and VQTP all showed significant upward trends from 2000 to 2018, and there were obvious differences in the fluctuation range and change rate between them, among which the AS had the fastest change rate. (3) From 2000 to 2018, the NDVI of the AM, AS and VQTP all showed upward trends in each season. The NDVI of the AM and VQTP increased at the fastest rate in the spring, while the NDVI of the AS increased at the fastest rate during the growing season. The peak growth periods of the AM, AS and VQTP were all in the summer. The fluctuation trends of the NDVI for the AM, AS and VQTP were relatively consistent during different seasons, indicating that the AM and AS played key and decisive roles in the overall vegetation change trend on the Qinghai-Tibet Plateau.
(4) The growing season NDVI showed an increasing trend in most regions of the Qinghai-Tibet Plateau from 2000 to 2018, among which the significantly increased regions were mainly distributed along the midwestern and northeastern regions of the Qinghai-Tibet Plateau, with the Qilian Mountains being the most typical region. The areas of significant decline were mainly distributed among the surrounding areas of Namtso and Siling Co in Tibet and the source region of the Yangtze River in Qinghai Province. The overall improvement in vegetation on the Qinghai-Tibet Plateau was mainly dominated by AM and AS, while vegetation degradation primarily occurred in the AM, and the seriously degraded regions were largely concentrated in areas where AM was distributed. (5) The significant positive feedback region of the AM on annual precipitation was primarily between the 400 mm and 600 mm isohyet areas, while the positive feedback region of the AS on annual precipitation was mostly distributed around the 400 mm isohyet area. The regions with a significant positive correlation between AM and AS and annual mean temperatures were all located between the −5 • C and 1 • C isotherm areas, and the significant impact of the temperature on the vegetation of the Qinghai-Tibet Plateau was mostly based on alpine grassland. The correlation between the NDVI and climatic factors in different regions and different vegetation types had great spatial heterogeneity.