As a key biotic factor, phenology exerts fundamental influences on ecosystem carbon sequestration. However, whether spring phenology affects the subsequent seasonal ecosystem productivity and the underlying resource limitation mechanism remains unclear for the alpine grasslands of the Tibetan Plateau (TP). In this study, we investigated the direct and lagged seasonal responses of net primary productivity (NPP) to the beginning of growing season (BGS) along a precipitation gradient by integrating field observations, remote sensing monitoring and ecosystem model simulations. The results revealed distinct response patterns of seasonal NPP to BGS. Specifically, the BGS showed a significant and negative correlation with spring NPP (R = −0.73, p < 0.01), as evidenced by the direct boosting effects of earlier BGS on spring NPP. Moreover, spring NPP was more responsive to BGS in areas with more annual precipitation. The boosting effects of earlier BGS on NPP tended to weaken in summer compared with that in spring. Sequentially, BGS exhibited stronger positive correlation with autumn NPP in areas with less annual precipitation, which suggested the enhanced lagged suppressing effects of earlier spring phenology on ecosystem carbon assimilation during the later growing season under aggravated water stress. Overall, the strengthened NPP in spring was offset by its decrement in autumn, resulting in no obvious relationship between BGS and annual NPP (R = −0.34, p > 0.05) for the entire grasslands on the TP. The findings of this study imply that the lagged effects of phenology on the ecosystem productivity during the subsequent seasons should not be neglected in the future studies.

Climate change variation on a small scale may alter the underlying processes determining a pattern operating at large scale and vice versa. Plant response to climate change on individual plant levels on a fine scale tends to change population structure, community composition and ecosystem processes and functioning. Therefore, we reviewed the literature on plant response and resilience to climate change in space and time at different scales on the Tibetan Plateau. We report that spatiotemporal variation in temperature and precipitation dynamics drives the vegetation and ecosystem function on the Tibetan Plateau (TP), following the water–energy dynamics hypothesis. Increasing temperature with respect to time increased the net primary productivity (NPP) on most parts of the Tibetan Plateau, but the productivity dynamics on some parts were constrained by 0.3 °C decade⁻¹ rising temperature. Moreover, we report that accelerating studies on plant community assemblage and their contribution to ecosystem functioning may help to identify the community response and resilience to climate extremes. Furthermore, records on species losses help to build the sustainable management plan for the entire Tibetan Plateau. We recommend that incorporating long-term temporal data with multiple factor analyses will be helpful to formulate the appropriate measures for a healthy ecosystem on the Tibetan Plateau.

Vegetation phenology is considered a sensitive indicator of climate change, which controls carbon, nitrogen, and water cycles within terrestrial ecosystems. The Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) is an important moderate resolution remote sensing data for monitoring vegetation phenology. However, Terra MODIS Collection 5 (C5) vegetation index products were identified to be affected by sensor degradation, which has been addressed in the recently released MODIS Collection 6 (C6) vegetation index products. In order to compare the difference between MODIS C5 and C6 NDVI in monitoring vegetation phenology, the start and end of growing season (SOS and EOS) of the alpine grassland on the Tibetan Plateau (TP) were extracted using four common methods. Then, the C5 and C6 NDVI-derived SOS (SOS_C5 and SOS_C6) and EOS (EOS_C5 and EOS_C6) were compared with ground-observed phenology data. Results showed that the multi-year average growing season NDVIs of C6 were lower than those of C5 in most areas, while the inter-annual variation patterns of regional average SOS_C5 and SOS_C6 (EOS_C5 and EOS_C6) were consistent. However, large spatial differences in phenological trends were found between C5 and C6 NDVI products. From C5 to C6, pixels with a SOS (EOS) trend shifting from significant to insignificant or from insignificant to significant accounted for at least 14.58% (9.07%) of the total pixels. SOS_C5 was more consistent than SOS_C6 with the ground-observed green-up dates. C5 NDVI may be more appropriate for monitoring SOS than C6 NDVI in the study region, but more ground-observed phenology records are needed to confirm it due to only four observational sites in this study. However, large differences and poor correlations existed between EOS_C5 (EOS_C6) and the ground-observed beginning of leaf coloring. To further evaluate the uncertainty of MODIS C5 and C6 NDVI in monitoring vegetation phenology, higher resolution near-surface remote sensing data and corresponding validation methods should be applied.