Linking Vegetation Phenology to Net Ecosystem Productivity: Climate Change Impacts in the Northern Hemisphere Using Satellite Data

With global climate change, linking vegetation phenology with net ecosystem productivity (NEP) is crucial for assessing vegetation carbon storage capacity and predicting terrestrial ecosystem changes. However, there have been few studies investigating the relationship between vegetation phenology and NEP in the middle and high latitudes of the Northern Hemisphere. This study comprehensively analyzed vegetation phenological changes and their climate drivers using satellite data. It also investigated the spatial distribution and climate drivers of NEP and further analyzed the sensitivity of NEP to vegetation phenology. The results indicated that the average land surface phenology (LSP) was dominated by a monotonic trend in the study area. LSP derived from different satellite products and retrieval methods exhibited relatively consistent responses to climate. The average SOS and POS for different retrieval methods showed a higher negative correlation with nighttime temperatures compared to daytime temperatures. The average EOS exhibited a higher negative correlation with daytime temperatures than a positive correlation. The correlations between VPD and the average SOS, POS, and EOS showed that the proportion of negative correlations was higher than that of positive correlations. The average annual NEP ranged from 0 to 1000 gC·m−2. The cumulative trends of NEP were mainly monotonically increasing, accounting for 61.04%, followed by monotonically decreasing trends, which accounted for 17.95%. In high-latitude regions, the proportion of positive correlation between VPD and NEP was predominant, while the proportion of negative correlation was predominant in middle-latitude regions. The positive and negative correlations between soil moisture and NEP (48.08% vs. 51.92%) were basically consistent in the study area. The correlation between SOS and POS with NEP was predominantly negative. The correlation between EOS and NEP was overall characterized by a greater proportion of negative correlations than positive correlations. The correlation between LOS and NEP exhibited a positive relationship in most areas. The sensitivity of NEP to vegetation phenological parameters (SOS, POS, and EOS) was negative, while the sensitivity of NEP to LOS was positive (0.75 gC·m−2/d for EVI vs. 0.63 gC·m−2/d for LAI vs. 0.30 gC·m−2/d for SIF). This study provides new insights and a theoretical basis for exploring the relationship between vegetation phenology and NEP under global climate change.