The accumulation and ablation processes of seasonal snow significantly affect the land surface phenology in a mountainous ecosystem. However, the ability of snow to regulate the alpine land surface phenology in the arid regions is not well described in the context of climate change. The impact of snowpack changes on land surface phenology and its driving factors were investigated in the Tianshan Mountains using the land surface phenology metrics derived from satellited products and a snow dataset from downscaled regional climate model simulations covering the period from 1983 to 2015. The results demonstrated that the annual mean start of growing season (SOS) and length of growing season (LOS) experienced a significant (p < 0.05) decrease and increase with a rate of -2.45 days/decade and 2.98 days/decade, respectively. The significantly advanced SOS and increased LOS were mainly seen in the Western Tianshan Mountains and Ili Valley regions with elevations from 2500 to 3500 m a.s.l and below 3000 m a.s.l, respectively. During the early spring, the significant decline in snow cover fraction (SCF) could advance the SOS. In contrast, snowmelt amount and annual maximum snow water equivalent (SWE) have an almost equally substantial positive correlation with annual maximum vegetation greenness. In particular, the SOS of grassland was the most sensitive to variations of snow cover fraction during early spring than that of other vegetation types, and their strong relationship was mainly located at elevations from 1500 to 2500 m a.s.l. Its greenness was significantly controlled by the annual maximum snow water equivalent in all elevation bands. Both decreased SCF and increased temperature in the early spring caused a significant advance of the SOS, consequently prolonging the LOS. Meanwhile, more SWE and snowmelt amount could significantly promote vegetation greenness by regulating the soil moisture. The results can improve the understanding of the snow ecosystem services in the alpine regions under climate change.

Recent efforts have been made to monitor the seasonal metrics of plant canopy variations globally from space, using optical remote sensing. However, phenological estimations based on vegetation indices (VIs) in high-latitude regions such as the pan-Arctic remain challenging and are rarely validated. Nevertheless, pan-Arctic ecosystems are vulnerable and also crucial in the context of climate change. We reported the limitations and challenges of using MODerate-resolution Imaging Spectroradiometer (MODIS) measurements, a widely exploited set of satellite measurements, to estimate phenological transition dates in pan-Arctic regions. Four indices including normalized vegetation difference index (NDVI), enhanced vegetation index (EVI), phenology index (PI), plant phenological index (PPI) and a MODIS Land Cover Dynamics Product MCD12Q2, were evaluated and compared against eddy covariance (EC) estimates at 11 flux sites of 102 site-years during the period from 2000 to 2014. All the indices were influenced by snow cover and soil moisture during the transition dates. While relationships existed between VI-based and EC-estimated phenological transition dates, the R-2 values were generally low (0.01-0.68). Among the VIs, PPI-estimated metrics showed an inter-annual pattern that was mostly closely related to the EC-based estimations. Thus, further studies are needed to develop region-specific indices to provide more reliable estimates of phenological transition dates.

Alpine vegetation plays a crucial role in global carbon cycle. Snow cover is an essential component of alpine land cover and shows high sensitivity to climate change. The Tibetan Plateau (TP) has a typical alpine vegetation ecosystem and is rich of snow resources. With global warming, the snow of the TP has undergone significant changes that will inevitably affect the growth of alpine vegetation, but observed evidence of such interaction is limited. In particular, a comprehensive understanding of the responses of alpine vegetation growth to snow cover variability is still not well characterized on TP region. To investigate this, we calculated three indicators, the start (SOS) and length (LOS) of growing season, and the maximum of normalized difference vegetation index (NOVImax) as proxies of vegetation growth dynamics from the Moderate Resolution Imaging Spectroradiometer (MODIS) data for 2000-2015. Snow cover duration (SCD) and melt (SCM) dates were also extracted during the same time frame from the combination of MODIS and the Interactive Multi-sensor Snow and Ice Mapping System (IMS) data. We found that the snow cover phenology had a strong control on alpine vegetation growth dynamics. Furthermore, the responses of SOS, LOS and NDVImax to snow cover phenology varied among different biomes, eco-geographical zones, and temperature and precipitation gradients. The alpine steppes showed a much stronger negative correlation between SOS and SCD, and also a more evidently positive relationship between LOS and SCD than other types, indicating a longer SCD would lead to an earlier SOS and longer LOS. Most areas showed positive correlation between SOS and SCM, while a contrary response was also found in the warm but drier areas. Both SCD and SCM showed positive correlations with NDVImax, but the relationship became weaker with the increase of precipitation. Our findings provided strong evidence between vegetation growth and snow cover phenology, and changes in snow cover should be also considered when analyzing alpine vegetation growth dynamics in future,