Detecting the Spatiotemporal Variation of Vegetation Phenology in Northeastern China Based on MODIS NDVI and Solar-Induced Chlorophyll Fluorescence Dataset

Vegetation phenology is a crucial biological indicator for monitoring changes in terrestrial ecosystems and global climate. Currently, there are limitations in using traditional vegetation indices for phenology monitoring (e.g., greenness saturation in high-density vegetation areas). Solar-induced chlorophyll fluorescence (SIF), a novel remote sensing product, has great potential in depicting seasonal vegetation dynamics across various regions with different vegetation covers and latitudes. In this study, based on the GOSIF and MODIS NDVI data from 2001 to 2020, we extracted vegetation phenological parameters in Northeastern China by using Double Logistic (D-L) fitting function and the dynamic threshold method. Then, we analyzed the discrepancy in phenological period and temporal trend derived from SIF and NDVI data at multiple spatiotemporal scales. Furthermore, we explored the response of vegetation phenology to climate change and the persistence of phenological trends (Hurst exponent) in Northeastern China. Generally, there is a significant difference in trends between SIF and NDVI, but with similar spatial patterns of phenology. However, the dates of key phenological parameters are distinct based on SIF and MODIS NDVI data. Specifically, the start of season (SOS) of SIF started later (about 10 days), and the end of season (EOS) ended earlier (about 36 days on average). In contrast, the fall attenuation of SIF showed a lag process compared to NDVI. This implies that the actual period of photosynthesis, that is, length of season (LOS), was shorter (by 46 days on average) than the greenness index. The position of peak (POP) is almost the same between them. The great difference in results from SIF and NDVI products indicated that the vegetation indexes seem to overestimate the time of vegetation photosynthesis in Northeastern China. The Hurst exponent identified that the future trend of SOS, EOS, and POP is dominated by weak inverse sustainability, indicating that the future trend may be opposite to the past. The future trend of LOSSIF and LOSNDVI are opposite; the former is dominated by weak inverse sustainability, and the latter is mainly weak positive sustainability. In addition, we speculate that the difference between SIF and NDVI phenology is closely related to their different responses to climate. The vegetation phenology estimated by SIF is mainly controlled by temperature, while NDVI is mainly controlled by precipitation and relative humidity. Different phenological periods based on SIF and NDVI showed inconsistent responses to pre-season climate. This may be the cause of the difference in the phenology of SIF and NDVI extraction. Our results imply that canopy structure-based vegetation indices overestimate the photosynthetic cycle, and the SIF product can better track the phenological changes. We conclude that the two data products provide a reference for monitoring the phenology of photosynthesis and vegetation greenness, and the results also have a certain significance for the response of plants to climate change.