Alpine Grassland Reviving Response to Seasonal Snow Cover on the Tibetan Plateau

: Season snow cover plays an important role in vegetation growth in alpine regions. In this study, we analyzed the spatial and temporal variations in seasonal snow cover and the start of the growing season (SOS) of alpine grasslands and preliminarily studied the mechanism by which snow cover affects SOS changes by modifying the soil temperature (ST) and soil moisture (SM) in spring. The results showed that signiﬁcant interannual trends in the SOS, snow end date (SED), snow cover days (SCD), ST, and SM existed over the Tibetan Plateau (TP) in China from 2000 to 2020. The SOS advanced by 2.0 d/10 a over the TP over this period. Moreover, the SOS showed advancing trends in the eastern and central parts of the TP and a delayed trend in the west. The SED and SCD exhibited an advancing trend and a decreasing trend in high-elevation areas, respectively, and the opposite trends in low-elevation areas. The ST showed a decreasing trend in low-elevation areas and an increasing trend in high-elevation areas. The SM tended to increase in most areas. The effects of the seasonal snow cover on the ST and SM indirectly inﬂuenced the SOS of alpine grasslands. The delayed SEDs and more SCD observed herein could provide increasingly wet soil conditions optimal for the advancement of the SOS, while less snow and shorter snow seasons could delay the SOS of alpine grasslands on the TP. SED and SCD were both signiﬁcantly negatively correlated with the ST, suggesting that fewer SCD and earlier SEDs contributed to the increased STs, while more SCD and later SEDs inhibited soil warming. The SED was negatively correlated with the SM, and SCD was positively correlated with the SM, respectively. This result indicated that the more SCD there were, the earlier the SED occurred, and the more snowmelt water replenished the soil in spring, thus increasing the SM. Negative correlations were found between SOS and ST and between SOS and SM, indicating that relatively high STs and SMs promoted SOS advancement, and vice versa.


Introduction
Vegetation phenology refers to the annual recurrence of vegetation growth during the life cycle of vegetation; this pattern is influenced by changes in the seasonal climate conditions and other environmental factors [1,2]. As a macroscopic and integrated representation of seasonal vegetation turnover, vegetation phenology can characterize the combined effects of environmental changes on plant life processes over time and is a key parameter used to monitor the response mechanisms of biological processes to global climate change [3]. Vegetation phenology affects the seasonal changes in biogeochemical cycling and biophysical properties through factors such as the reflectance, surface roughness, canopy conductivity, moisture, energy, CO 2 , and biogenic volatile organic matter [4][5][6].
Currently, global warming is an unignorably environmental issue that will continue in the coming decades [7]. Since 1860, the global average temperature has reportedly increased by nearly 1.0 • C [8,9], and this trend is more pronounced in high-latitude and high-elevation areas. The Tibetan Plateau (TP) warmed by approximately 0.4 • C/10 a ( • C/10 a means a temperature change over 10 years) from 1981 to 2011, almost twice the average rate of global warming [10]. Vegetation phenology is regarded as the most intuitive and sensitive integrated indicator of global climate change and is also an important observable indicator of the responses of ecosystems to climate change, which, in turn, determine the behavior of land surface carbon fluxes [11]. The sensitivity of vegetation phenology to climate change is more pronounced and observable in alpine regions than in low-elevation regions [12,13]. Dynamic variations in vegetation phenology affect the

Study Area
The TP is located in southwestern China and spans from 76°-104°E and from 25°-40°N [33], with a length of 2800 km from east to west and a width of approximately 300-1500 km from north to south, covering an overall area of approximately 2.6 × 10 6 km 2 [34]. Due to the large elevation gradient and complex topography on the TP, the plateau is rich in vegetation types and contains a wide range of plant species. The TP vegetation exhibits obvious horizontal and vertical zonal distribution patterns from southeast to northwest, varying from forests, scrublands, meadows, and grasslands to deserts [35]. Grasslands compose the main vegetation type on the TP and account for 62.2% of its total area ( Figure  1b). The TP is rich in snow, and snowmelt plays an important role in regulating river runoff and in ensuring the water supply of western China; because of this, the TP is known as the "Asian Water Tower" [33,35]. The TP is one of the most sensitive ecosystems to global changes; because of its unique geographical features and low level of human interference, it is an ideal area for studying the response mechanisms of spatial and temporal alpine grassland phenology heterogeneities to seasonal snow cover changes [36,37].

Vegetation Index Data
Moderate-resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) data were obtained from the MOD13Q1 V006 dataset released by the United States Geological Survey (USGS; https://lpdaac.usgs.gov/products/mod13q1v006/, accessed on 2 April 2021) at a temporal resolution of 16 d, a spatial resolution of 500 m, and a temporal range of 28 February 2000 to 31 December 2020. These data were preprocessed through stitching and reprojection steps and used to extract the SOS of alpine grasslands on the TP.

Snow Cover Data
The utilized MODIS daily cloud-free snow cover dataset was obtained from the National Cryosphere Desert Data Center (http://www.ncdc.ac.cn, accessed on 8 April 2021). This dataset was based on two MODIS reflectance products, MOD/MYD09GA, and a

Vegetation Index Data
Moderate-resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) data were obtained from the MOD13Q1 V006 dataset released by the United States Geological Survey (USGS; https://lpdaac.usgs.gov/products/mod13q1v006/, accessed on 2 April 2021) at a temporal resolution of 16 d, a spatial resolution of 500 m, and a temporal range of 28 February 2000 to 31 December 2020. These data were preprocessed through stitching and reprojection steps and used to extract the SOS of alpine grasslands on the TP.

Snow Cover Data
The utilized MODIS daily cloud-free snow cover dataset was obtained from the National Cryosphere Desert Data Center (http://www.ncdc.ac.cn, accessed on 8 April 2021). This dataset was based on two MODIS reflectance products, MOD/MYD09GA, and a multi-index snow cover mapping algorithm was developed in consideration of different land cover types to improve the identification accuracy of snow-covered areas in forests and mountainous regions. Moreover, the hidden Markov random field (HMRF) algorithm and multisource data fusion method were used to remove clouds from the daily snow cover products. The reported overall accuracy of this dataset was above 93% [38], confirming that this dataset can be used in snow cover change analyses with high credibility.

Soil Temperature and Moisture Data
Daily soil temperature data were obtained from the National Tibetan Plateau Data Center (TPDC; http://data.tpdc.ac.cn/en, accessed on 17 October 2021) at a temporal resolution of twice daily, divided into daytime and nighttime data, and a spatial resolution of 1 km from 2000-2020; this dataset covered the spatial range from 72 • E-104 • E and from 20 • N-45 • N and contained values in units of 0.01 K. A novel surface temperature time decomposition model was used to produce the daily soil temperature products, with input data such as Aqua MODIS ST products, GLDAS, and auxiliary data. The dataset is spatially seamless with good image quality and accuracy. Using MODIS ST as a reference, the mean bias error and the standard deviation of this dataset are 0.08-0.16 K and 1.12-1.46 K, respectively [39].
Monthly SM data were also obtained from the TPDC at a monthly temporal resolution, a 0.05 • spatial resolution, a spatial range of 72 • E-142 • E and 16 • N-56 • N, and a temporal range of July 2002 to December 2018 (these data were in m 3 /m 3 ). The SM dataset is a fusion of products and composited the Microwave Scanning Radiometer-Earth Observing System (AMSR-E), the Advanced Microwave Scanning Radiometer 2 (AMSR2), and the Soil Moisture and Ocean Salinity (SMOS) mission soil moisture products. To overcome the deficiencies of the passive microwave SM products with low resolution, the dataset is constructed using a spatially weighted decomposition algorithm to generate the 0.05 • spatial resolution product. The SM products were consistent with the in situ measurements (R > 0.78, RMSE < 0.05 m 3 /m 3 ), confirming the good accuracy throughout the time series [40].
For this study, the units of the SM data were first converted to K; then, the STs were calculated by taking the average value between the daytime and nighttime data. Then, the springtime ST and SM values were calculated (from March to May of the Julian year) over the TP.

Auxiliary Data
The 3572 samples utilized herein were all collected from field sampling sites in TP grasslands from 2003 to 2019, of which 1338 samples were from Xia et al. [41]. The locations of these sampling sites on the TP are shown in Figure 1. The samples were used in this study to ensure the reliability of the correlation analysis results derived between snow cover and grassland phenology to avoid uncertainties caused by land cover type map misclassifications.

Start of Growing Season (SOS)
Although MOD13Q1 is already a time series product after the noise removal step and 16 d synthesis process; residual fluctuations and noise still exist in these data due to the presence of clouds and atmospheric, sensor, and surface reflectivity. Therefore, in this study, smooth NDVI variation curves were obtained using the Savitzky-Golay filtering algorithm (S-G) with a filter window size of 4 [42,43]. Then, the dynamic threshold method [44,45] with a threshold value of 0.2 was used to extract the SOS of alpine grasslands on the TP [46,47].

Snow Melting Period Extraction
The hydrological year spans from September 1 of a given year to August 31 of the following year; in this study, the analyzed time range is 2000 to 2020, covering a total of 20 snow seasons. Based on daily cloud-free snow cover data, the snow end date (SED) and snow-covered days (SCD) were calculated according to the algorithm proposed by [48]. The utilized formulas are as follows: where Fd is the fixed date of March 1, SCD Fd denotes the SCD after March 1, n is the number of days, and S i is 0 or 1 corresponding to a nonsnow covered area or a snow covered area, respectively.

Tendency and Impact Analyses
Currently, slope trend analysis is the most commonly used method for visually exploring the temporal and spatial trends of SOS, SED, SCD, ST, and SM. The relevant formula is as follows: where Slope is the trend of the variable; y is the SOS, SED, SCD, ST, or SM value in year i; i is the annual variable varying from 1 to n, and n is the number of years.
To further determine the degree of interannual variation of each variable, F tests were performed on the slope variation trends to determine whether these trends were significant at the 95% confidence level. The corresponding formula is as follows: where i is the number of years, y i is the SOS, SED, SCD ST, or SM value in year i, y i is the fitted value, and y is the average SOS, SED, SCD, ST, or SM value. Pearson correlation analyses are usually used to measure linear correlations between two variables. In this study, this method was used to analyze the correlations between snow cover changes and the ST and between snow cover changes and the SM as well as the relationship between the ST/SM and alpine grassland phenology on the TP. The formula is as follows: where r is the correlation coefficient, x i is the SOS, SED, or SCD value in the corresponding year i, y i indicates the ST/SM/snow accumulation variable, x and y are the means of the corresponding variables in the corresponding year, and n is the number of years. The Pearson correlation coefficient varies from −1 ≤ r ≤ 1, and r values > 0 indicate positive correlations between variables, while r values < 0 indicate negative correlations between variables. The larger the absolute value of r is, the stronger the correlation is, while the closer to 0 the value is, the weaker the correlation is. An r value equal to 0 indicates that no linear correlation exists between the analyzed variables. Figure 2 shows the intra-annual variations in NDVI, snow cover extent (SCE), ST, and SM on the TP from 2000 to 2020. The annual mean NDVI had a minimum value on February 18 and then rose rapidly after April, reaching a maximum value of 0.6 on July 27 in the alpine grasslands; then, it began to decline gradually. The lowest NDVI value appeared during the period corresponding to the maximum SCE in mid-February. With the gradual decrease in SCE, the NDVI started to increase. The NDVI reached a maximum in July at the peak grass-growing stage, when SCE was the smallest. The snow cover began Figure 2 shows the intra-annual variations in NDVI, snow cover extent (SC and SM on the TP from 2000 to 2020. The annual mean NDVI had a minimum v February 18 and then rose rapidly after April, reaching a maximum value of 0.6 27 in the alpine grasslands; then, it began to decline gradually. The lowest NDV appeared during the period corresponding to the maximum SCE in mid-Februar the gradual decrease in SCE, the NDVI started to increase. The NDVI reached a ma in July at the peak grass-growing stage, when SCE was the smallest. The snow cove to accumulate gradually in September. With increasing SCE, the NDVI began to rapidly, and the end of the growing season was recorded around mid-October. The minimum ST of 266.7 K appeared on January 11 and then began to rise gra on March 9, the soil temperature rose above the freezing point, and on July 26, th mum soil temperature reached 286.7 K. Then, the ST began to decline gradually below freezing on November 14. The annual average SM in January was 0.05 m 3 /m which it increased slowly. Then, the SM increased slowly, possibly due to the in soil water content caused by melting snow. The SM increased rapidly in May, main plemented by liquid precipitation and thawing permafrost, and the grasslands gr turned green, thus enhancing the soil water holding capacity [49]. The SM reached imum value of 0.2 m 3 /m 3 in July and then gradually decreased. After September, decreased rapidly, likely due to the increased SCE. Then, the ST gradually droppe freezing point, and the SM declined sharply to the minimum value of 0.05 m 3 /m cember.

Intra-Annual Variations in Snow Cover Extent (SCE), NDVI, ST, and SM
In summary, there was good consistency among the changes in snow cover, SM, as well as the intra-annual changes in alpine grassland. In general, as snow the reduction in SCE promoted the absorption of solar radiation energy by the la face, the ST started to rise, the soil thawed, and the SM started to increase, thus pr suitable hydrothermal conditions for grassland growth. The minimum ST of 266.7 K appeared on January 11 and then began to rise gradually; on March 9, the soil temperature rose above the freezing point, and on July 26, the maximum soil temperature reached 286.7 K. Then, the ST began to decline gradually, falling below freezing on November 14. The annual average SM in January was 0.05 m 3 /m 3 , after which it increased slowly. Then, the SM increased slowly, possibly due to the increased soil water content caused by melting snow. The SM increased rapidly in May, mainly supplemented by liquid precipitation and thawing permafrost, and the grasslands gradually turned green, thus enhancing the soil water holding capacity [49]. The SM reached a maximum value of 0.2 m 3 /m 3 in July and then gradually decreased. After September, the SM decreased rapidly, likely due to the increased SCE. Then, the ST gradually dropped to the freezing point, and the SM declined sharply to the minimum value of 0.05 m 3 /m 3 in December.
In summary, there was good consistency among the changes in snow cover, ST, and SM, as well as the intra-annual changes in alpine grassland. In general, as snow melted, the reduction in SCE promoted the absorption of solar radiation energy by the land surface, the ST started to rise, the soil thawed, and the SM started to increase, thus providing suitable hydrothermal conditions for grassland growth.

Snow Cover
Snow cover affects the SOS of alpine grasslands through the influences of snowmelt amount and timing variations. In this study, the SED and SCD were used to analyze the effects of snow cover on the SOS of alpine grasslands. Figure 3 shows that the pattern and trends of SED and SCD exhibited obvious spatial heterogeneity on the TP from 2000 to 2020. The annual SED ranged from March 1 to August 31, with an average SED of March 27. The SED showed a clear vertical zonal distribution and was gradually delayed as the elevation gradient rose. In addition, there was no obvious snowmelt date in high-mountain areas, such as the Karakoram, Kunlun, Himalayas, and Nyingchi Tanggula, as these areas have permanent snow cover. The SED was delayed by 1.2 d/10 a on average over 21 years on the TP. Almost 61.0% of the study regions showed trends of delayed SEDs, and significant delays accounted for 8.3% of these, most of which were distributed on the Pamir Plateau, in the Qilian Mountains, in the Bayankara Mountains, and in the Himalayas and Gangdis Mountains. The areas with advancing SED trends accounted for 38.8% of the region, mainly on the northwestern plateau in Ngari Prefecture, in the Arjinshan Mountains, Nyingchi Tanggula, and in the southern part of the Hengduan Mountains. Areas with significant SED advancement accounted for only 0.9% of the study region, and these locations were mainly in the Karakorum Mountains and in the southern part of the Hengduan Mountains.
trends of SED and SCD exhibited obvious spatial heterogeneity on the TP from 2 2020. The annual SED ranged from March 1 to August 31, with an average SED of M 27. The SED showed a clear vertical zonal distribution and was gradually delayed elevation gradient rose. In addition, there was no obvious snowmelt date in high-m tain areas, such as the Karakoram, Kunlun, Himalayas, and Nyingchi Tanggula, as areas have permanent snow cover. The SED was delayed by 1.2 d/10 a on average o years on the TP. Almost 61.0% of the study regions showed trends of delayed SED significant delays accounted for 8.3% of these, most of which were distributed on t mir Plateau, in the Qilian Mountains, in the Bayankara Mountains, and in the Him and Gangdis Mountains. The areas with advancing SED trends accounted for 38.8% region, mainly on the northwestern plateau in Ngari Prefecture, in the Arjinshan M tains, Nyingchi Tanggula, and in the southern part of the Hengduan Mountains. with significant SED advancement accounted for only 0.9% of the study region, and locations were mainly in the Karakorum Mountains and in the southern part of the duan Mountains. Combined with the spatial distribution of elevation (Figure 1a), the SCD show increasing trend with increasing elevation on TP (Figure 3c). The average SCD on t from September to April in the hydrological year ranged from 0 to 242 d, with an av value of 54.8 d. The Karakorum, Kunlun, Qilian, Himalaya, and Hengduan Mou were the areas with the highest SCDs, and some high-mountain areas were covered snow all year round. Relatively small SCD values were distributed in the hinterland the northern Qaidam Basin due to the low water vapor inputs in these regions caus high-mountain obstructions. SCD increased by 2.2 d/10 a from 2000 to 2020 on the total of 39.3% of the TP showed a decreasing SCD trend, with 1.38% decreasing s cantly, mainly in the southern part of the Hengduan Mountains and the Karak Combined with the spatial distribution of elevation (Figure 1a), the SCD showed an increasing trend with increasing elevation on TP (Figure 3c). The average SCD on the TP from September to April in the hydrological year ranged from 0 to 242 d, with an average value of 54.8 d. The Karakorum, Kunlun, Qilian, Himalaya, and Hengduan Mountains were the areas with the highest SCDs, and some high-mountain areas were covered with snow all year round. Relatively small SCD values were distributed in the hinterlands and the northern Qaidam Basin due to the low water vapor inputs in these regions caused by high-mountain obstructions. SCD increased by 2.2 d/10 a from 2000 to 2020 on the TP. A total of 39.3% of the TP showed a decreasing SCD trend, with 1.38% decreasing significantly, mainly in the southern part of the Hengduan Mountains and the Karakorum Mountains. The SCD in 59.0% of the study area showed an increasing trend, of which 4.7% increased significantly. The regions with decreasing SCDs largely overlapped with the regions with advancing SEDs on the plateau, while areas with increasing SCDs largely overlapped with those with delayed SEDs. The results showed that snow cover changes on the TP exhibited spatial distribution patterns with two scenarios corresponding to increased or decreased snow cover, the spatial distributions of these patterns were basically stable, and the proportions of both were also the same. Overall, the SED in the central and southern parts of the plateau had been delayed over the past 20 years, while SCD has been increasing. In the other regions, the SED has advanced, while the SCD has been reduced.

Start of the Growing Season (SOS)
The spatial distribution characteristics of the SOS of alpine grasslands on the TP were consistent with the elevation and hydrothermal conditions described above and showed a gradual delayed trend from southeast to northwest (Figure 4). The annual average SOS ranged from the 89th day to the 192nd day, with an average value of the 127th day (early May) of the Julian year from 2000 to 2020. A total of 78.1% of the regions entered the greening period before June. The earliest SOS dates were observed in the eastern valley region, especially in Aba Prefecture of Sichuan Province, where the average SOS occurred on approximately April 10. The latest SOS dates were observed after June on the southwestern TP, especially along the northern Himalayas. The interannual SOS variation showed significant spatial differences on the TP over 21 years. In general, the early SOSs observed in the northeast were more extensive than those in the south, while delayed SOSs were found mainly in the relatively high-elevation areas of the southwestern region. Among the 65.4% of TP that exhibited advancing SOS trends, approximately 4.9% showed significant advances (p < 0.05), mainly in the Yellow River source area, with an average advance of 4.9 d/10 a; the maximum advancing rate reached 5.2 d/10 a (p < 0.05). Approximately 34.6% of the TP showed a delayed SOSs, with extremely small areas exhibiting significant delays. Overall, the average SOS over the TP advanced by 2.1 d/10 a from 2000 to 2020. Moreover, compared to the two different scenarios of SED and SCD, the regions with delayed SEDs and increased SCDs on TP, the SOS basically showed advanced trend. In contrast, SOS tended to delay in areas where the SEDs advanced and SCD decreased. stable, and the proportions of both were also the same. Overall, the SED in the central and southern parts of the plateau had been delayed over the past 20 years, while SCD has been increasing. In the other regions, the SED has advanced, while the SCD has been reduced.

Start of the Growing Season (SOS)
The spatial distribution characteristics of the SOS of alpine grasslands on the TP were consistent with the elevation and hydrothermal conditions described above and showed a gradual delayed trend from southeast to northwest (Figure 4). The annual average SOS ranged from the 89th day to the 192nd day, with an average value of the 127th day (early May) of the Julian year from 2000 to 2020. A total of 78.1% of the regions entered the greening period before June. The earliest SOS dates were observed in the eastern valley region, especially in Aba Prefecture of Sichuan Province, where the average SOS occurred on approximately April 10. The latest SOS dates were observed after June on the southwestern TP, especially along the northern Himalayas. The interannual SOS variation showed significant spatial differences on the TP over 21 years. In general, the early SOSs observed in the northeast were more extensive than those in the south, while delayed SOSs were found mainly in the relatively high-elevation areas of the southwestern region. Among the 65.4% of TP that exhibited advancing SOS trends, approximately 4.9% showed significant advances (p < 0.05), mainly in the Yellow River source area, with an average advance of 4.9 d/10 a; the maximum advancing rate reached 5.2 d/10 a (p < 0.05). Approximately 34.6% of the TP showed a delayed SOSs, with extremely small areas exhibiting significant delays. Overall, the average SOS over the TP advanced by 2.1 d/10 a from 2000 to 2020. Moreover, compared to the two different scenarios of SED and SCD, the regions with delayed SEDs and increased SCDs on TP, the SOS basically showed advanced trend. In contrast, SOS tended to delay in areas where the SEDs advanced and SCD decreased.   (Figure 5a). The ST increased by 0.9 K/10 a on the TP from 2000 to 2020. The areas with increasing ST trends accounted for 60.9% of the whole plateau, of which the ST increased significantly in 0.6% (p < 0.05), mainly distributed in the northern Karakorum Mountains and the southern foothills of the Himalayas with increasing trends of >0.5 K/10 a. On the northern TP, in the Himalayas, the Gangdis Mountains, and the Yellow River source area, areas with decreasing ST trends accounted for 38.2% of the study area. However, only part of southeastern Tibet exhibited significant decreases, accounting for approximately 0.3%. Comparing the spatial and temporal changes in SED and SCD, regions with delayed SEDs and increased SCDs on the TP were found to have had lower STs (Figure 5b). In contrast, STs tended to increase in areas where the SEDs occurred earlier and the SCD decreased. This finding indicated that the amount of snow cover directly affected the ST. The ST was relatively low in areas with more and longer-lasting snow cover, and vice versa.
creasing SCDs. This result indicated that SED advancement directly leads to regional S increases. The SM increases were not significant in areas with increased SCD and delay SEDs due to the lack of replenished SM over time. However, this finding may have a been caused by the thawing of permafrost caused by rising surface temperatures, wh would further increase the SM content [49]. The SM tended to increase in 90.3% of the TP in spring, of which 58.6% of regions exhibited significant increases (p < 0.05). Areas with decreasing SM trends were mainly distributed around watersheds, in the Arjinshan, and in the Qaidam Basin, accounting for 8.4% of the study area (Figure 5d). Most areas with significantly increasing SM were located in regions with advancing SEDs and decreasing SCDs. This result indicated that SED advancement directly leads to regional SM increases. The SM increases were not significant in areas with increased SCD and delayed SEDs due to the lack of replenished SM over time. However, this finding may have also been caused by the thawing of permafrost caused by rising surface temperatures, which would further increase the SM content [49].

SOS Response to Snow Cover, Soil Moisture, and Temperature
Regions with increased SCD were associated with decreased ST and increased SM, and delayed SED was associated with decreased STs and increased SM, respectively. The SOS in these regions showed overall advancement. In contrast, in regions with reduced SCD and earlier SED, the ST increased and the SM increased significantly, while the SOS basically showed delayed tendencies in spatial. This result suggested that SCD and SED variations might lead to changes in regional ST and SM, thus inducing SOS alteration in alpine grasslands. The interaction of several effects need further validation. Figure 6 shows that the Pearson correlation coefficients derived among the SOS, SED, SCD, ST, and SM were all highly significant (p < 0.001). These results indicated that the SED and SCD were both significantly negatively correlated with the ST, suggesting that fewer SCD and earlier SEDs contributed to the increased STs, while more SCD and later SEDs inhibited soil warming. The SED was negatively correlated with the SM, and SCD was positively correlated with the SM, respectively. This result indicated that the more SCD there were, the earlier the SED occurred, and the more snowmelt water replenished the soil in spring, thus increasing the SM. Negative correlations were found between SOS and ST and between SOS and SM, indicating that relatively high STs and SMs promoted SOS advancement, and vice versa.
variations might lead to changes in regional ST and SM, thus inducing SOS alter alpine grasslands. The interaction of several effects need further validation. Figure 6 shows that the Pearson correlation coefficients derived among the SO SCD, ST, and SM were all highly significant (p < 0.001). These results indicated SED and SCD were both significantly negatively correlated with the ST, suggest fewer SCD and earlier SEDs contributed to the increased STs, while more SCD a SEDs inhibited soil warming. The SED was negatively correlated with the SM, a was positively correlated with the SM, respectively. This result indicated that th SCD there were, the earlier the SED occurred, and the more snowmelt water repl the soil in spring, thus increasing the SM. Negative correlations were found betw and ST and between SOS and SM, indicating that relatively high STs and SMs pr SOS advancement, and vice versa. In summary, more snow and the delay of melting time will inhibit the rise o can provide enough SM to promote the alpine grassland to turn green in adva reduction of snow cover and the advance of melting will lead to the delay of SOS, although less snow cover will promote the rise of soil temperature, the impact o relatively complex, and the result of the impact is completely opposite. This pro needs to supplement the precipitation and temperature data in these areas for furt ification. These findings again confirmed that the amount of snow and the melting were the main factors causing ST and SM changes and subsequently leading to in the SOS of alpine grasslands. Figure 6. Pearson correlation analysis results obtained for the start of the growing season, snow end date, snow-covered days, soil temperature, and moisture. *** denote an extremely significant correlation (p < 0.001).
In summary, more snow and the delay of melting time will inhibit the rise of ST but can provide enough SM to promote the alpine grassland to turn green in advance; the reduction of snow cover and the advance of melting will lead to the delay of SOS, because although less snow cover will promote the rise of soil temperature, the impact on SM is relatively complex, and the result of the impact is completely opposite. This process still needs to supplement the precipitation and temperature data in these areas for further verification. These findings again confirmed that the amount of snow and the melting timing were the main factors causing ST and SM changes and subsequently leading to changes in the SOS of alpine grasslands.

Discussion
The spatial distribution of snow cover showed vertical zonal distribution characteristics on the TP. More and longer-term snow cover was found in high-elevation mountainous areas where the SEDs advanced and the SCD decreased from 2000 to 2020. In low-elevation areas, the snow cover ranged from accumulation to complete melting throughout the course of a hydrological year, and alpine grasslands are widely distributed in these regions. Previous studies revealed that the increased SCD and delayed SEDs contributed to earlier SOSs, while the decreased SCD on the western TP promoted the delay of the SOS [50][51][52]. Those conclusions are consistent with the founding of this paper.
However, past studies mainly focused on the effects of snow on the SOS on the TP but did not reveal the mechanism by which snow affects other environmental factors, such as the ST or SM, thus indirectly influencing the SOS. It is well known that the presence of snow protects newborn seedlings in the spring when the vegetation conditions are suitable for sprouting and growth [53]. Snowmelt plays an important role in affecting the SM [54], and the bare ground exposed following snowmelt is warmed by the absorption of solar radiation [55]. Thus, the hydrothermal conditions during the SOS period are important factors affecting vegetation emergence and growth [56]. Thus, the purpose of this study, which was based on the premise that snow cover plays an important role in affecting the SOS of alpine grasslands, was to further reveal the effects that snow cover has on the SOS by influencing soil hydrothermal conditions.
In this study, we found that SCD reductions and advanced SEDs provided a more suitable soil hydrothermal environment for forage germination and growth. This was because the amount of snow cover affected the absorption of solar radiation energy at the land surface. In areas with relatively low snow cover and significant decreasing snow cover trends, the land surface could absorb more energy, thus prompting the ST to rise and, in turn, leading to advanced permafrost melting, increased SM, and, finally, an advanced SOS [49]. SOS trends may vary depending on the local plant community and the ecogeographic zone as the TP has complex topography, a large elevation gradient, and a wide range of plant species [57]. Yi et al. [58] found that the interaction mechanism between multiyear permafrost and alpine grassland changes varied greatly from region to region. In this variable climate system, the causes of SOS changes on the TP are complex [59]. Therefore, the generalizability of this effect of snow cover on the SOS by altering the ST and SM needs further study and verification in other regions on the TP.
We found that a single factor influencing an SOS change could also act in conjunction with other factors or by influencing other environmental factors. In this study, we found that snow affected the SOS by affecting the ST and SM conditions in spring. There are many other ways in which snow affects SOS through the ST and SM. For example, the effect of cold-season snow on the subsurface ST is one way by which snow affects SOS changes via the soil environment. Yu et al. [24] found that the winter and spring snow depths combined to affect vegetation growth by regulating the ST. In addition, the effect of snow cover through ST and SM changes on the SOS may also be influenced by climate warming. In other words, snow cover and climate factors interact to influence the SOS. Hoye et al. [60] found that the average temperature during the snowmelt period together with the average temperature during the period following snowmelt, but before the SOS, influenced the dynamics of the SOS in the Arctic. Furthermore, increased SMs resulting from snowmelt are accompanied by corresponding increases in evaporation, which in turn might cause precipitation to increase and the air temperature to decrease [61,62]. This persistent effect of snow cover on precipitation and air temperature conditions moderated through the soil may also affect the SOS.
In this study, we found that the snow cover influenced SOS changes by increasing the SM, and this process was complemented by snow melt and the thawing of permafrost. In addition, the occurrence of spring rain/snow events also critically impacts SM changes. Brandt et al. [63] found that rain-snow events penetrated the soil more quickly than snowmelt water because of the closer connections among spring rain, snow, and soil. In this study, the average SED in the sampled area was March 10. Most of the study area experienced snow season termination before April. Therefore, we speculate that the primary source of SM is snowmelt water, while the effect of overall spring precipitation may also contribute to the SM.

Conclusions
In this paper, the dynamics of vegetation revival, snow cover, SM, and ST in alpine grasslands in TP were examined from 2000 to 2020. By analyzing the effect of the seasonal snow cover on the soil hydrothermal conditions in spring, we clarified that there is a certain relationship between seasonal snow cover and SOS of alpine grassland, and the snow cover affects the SOS by affecting the soil temperature and moisture in spring. The main conclusions are summarized as follows: (1) There was both a significant negative correlation between the SOS with ST and SM in TP (p < 0.001), indicating that alpine grassland would undergo revival and growth earlier in warmer and wetter soil environment in spring.
(2) In TP, the SCD was significantly positively correlated with the SM, and significantly negatively correlated with the ST, suggesting that increased SCD can lead to spring soil temperature decrease and moisture increase in spring. (3) The SED was both significantly negatively correlated with the ST and SM. It indicated that the end of the snow season is earlier, and the spring soil would be warmer and wetter in TP. (4) Both the delayed SED and the increased SCD contributed to the advanced SOS, and the reduction of snow cover and the advance of melting will lead to the delay of SOS.
There is no doubt that the raised ST and increased SM in spring will lead to an earlier SOS. However, the spatial trends showed that the SM increased in the whole TP. It indicated that although the decreased SCD and delayed SED limited the increased SM, it is still replenished, probably by the permafrost thawing or liquid precipitation. In addition, the snow cover in TP is not continuous and melts very rapidly, which increases the complexity of analyzing the effect mechanism of how the changes of SOS of alpine grassland are affected by snow cover. In order to clarify the contribution of SED to the SOS of alpine grasslands, further analysis should be made in combination with spring precipitation and temperature changes in alpine grassland over the Tibetan Plateau.  Data Availability Statement: The NDVI datasets are available at https://lpdaac.usgs.gov/products/ mod13q1v006/ (accessed on 2 April 2021). The snow cover datasets are available at http://www. ncdc.ac.cn (accessed on 8 April 2021). The daily soil temperature and moisture datasets are available at http://data.tpdc.ac.cn (accessed on 17 October 2021).

Conflicts of Interest:
The authors declare no conflict of interest.