“…We used six long-term and three short-term satellite remote sensing datasets in this study, covering four types of vegetation indicators -NDVI, LAI, SIF and VOD (Table 1). These remote sensing products include: (a) VIP15 NDVI product (1981, which has 0.05°and 15-day resolutions and is developed by harmonizing the observations of Advanced Very High Resolution Radiometer (AVHRR) from 1981 to 1999 and Moderate Resolution Imaging Spectroradiometer (MODIS) C5 from 2000 to 2014 (Didan et al, 2015); (b) GIMMS NDVI3g , which has 15-day and 1/12°resolutions, and was produced by aggregating daily AVHRR surface reflectance (Pinzon & Tucker, 2014); (c) GIMMS LAI3g product , which was further produced by GIMMS NDVI3g product using a neural network algorithm (Zhu et al, 2013); (d) PKU GIMMS NDVI (1982NDVI ( -2020, which is a new version of GIMMS NDVI product produced by a machine learning model incorporating Landsat images (Li et al, 2023); (e) GLASS LAI product (1981, which has 8-day temporal resolution and 0.05°spatial resolution, and was reconstructed by combing AVHRR LAI from 1981 to 1999 and MODIS LAI from 2000 to 2018 using a bidirectional long short-term memory (Bi-LSTM) model (Ma & Liang, 2022); (f) GLOBMAP LAI (1982LAI ( -2019 dataset at a spatial resolution of ∼0.07°, covering the period from 1982 to 2019, has half-month (1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000) and 8-day (2001-2019) temporal resolutions, and was produced by establishing a pixel-level AVHRR Simple Ratio (SR)-MODIS LAI relationship (Liu et al, 2012); (g) MOD13C1 NDVI (2000-2020) product (C61), which has a 0.05°spatial resolution and a 16-day temporal resolution (Didan & Munoz, 2019); (h) OCO2 SIF product (2000 at resolutions of 0.05°and 4 days, which was generated from MODIS surface reflectance and OCO2 SIF data using a neural network approach…”