Additional file 1: Figure S1. The spatial patterns of the end of growing season (EGS) in subtropical vegetation in China from 2000 to 2018: NIRv dataset. a) derivative method and b) dynamic threshold method. Inset plots (the bottom-left of each figure) display the frequency distribution of EGS. Figure S2. Comparison of the end of growing season (EGS) retrieved from NIRv and flux data. Two ecological monitoring stations were listed, including Qianyanzhou station (QYZ) and Dinghushan station (DHS). Both methods were listed, including derivative and dynamic threshold methods. Figure S3. Linear trends of the end of growing season (EGS) across China’s Subtropical biomes from 2000 to 2018. Five main vegetation types in study area were listed, including evergreen coniferous forest (ECF), evergreen broadleaved forest (EBF), deciduous broadleaved forest (DBF), shrub, and grassland. a) derivative and b) dynamic threshold method. A negative value indicates an advance, and a positive value indicates a delay. * indicates statistically significant trends at the 90% (P