## Abbreviations \begin{table}
\begin{tabular}{c c} \hline \hline
**Plant Function Types** & \multicolumn{2}{c}{**Model Formula**} \\ \hline \multirow{4}{*}{Forest} & ENF & SOG = 684.40 \(-\) 115.13 \(\times\) mEVlsp – 165.15 \(\times\) sEVlsp – 70.16 \(\times\) EVI243 – 2.09 \(\times\) LST60 \\ \cline{2-3}  & DBF & SOC = 237.01 \(-\) 116.82 \(\times\) mEVlsp – 86.70 \(\times\) mEVlsu – 0.43 \(\times\) LST243 \\ \cline{2-3}  & MF & SOG = \(-\)73.94 \(-\) 94.49 \(\times\) EVI334 \(-\) 2.20 \(\times\) LST60 + 2.77 \(\times\) mLSTsu – 2362.37 \(\times\) cIStu \\ \hline \multirow{4}{*}{Non-forest} & CSH & \multirow{2}{*}{SOG = \(-\)51.57 \(-\) 243.87 \(\times\) mEVlsp – 117.05 \(\times\) EVI243 + 2.88 \(\times\) LST152 – 2.30 \(\times\) LST243 + 221.19 \(\times\) sEVlsu} \\ \cline{2-3}  & GRA & SOC = 180.18 – 180.55 \(\times\) mEVlsp – 207.39 \(\times\) sEVlsu – 3.25 \(\times\) sLStau \\ \cline{2-3}  & CRO & SOC = 88.28 – 114.86 \(\times\) mEVlsp – 141.83 \(\times\) maxEVlsu \\ \cline{2-3}  & WET & SOC = 521.47 – 1.40 \(\times\) mLSTsp – 352.64 \(\times\) sEVlsu – 62.34 \(\times\) mEVlsp \\ \hline \hline \end{tabular}
\end{table} 
Table 10: Model formulas for estimating the starting days of GPP (SOG). Evergreen needle-leaf forest (ENF), deciduous broadleaf forest (DBF), mixed forest (MF), closed shrublands (CSH), open shrublands (OSH), savannas (SAV), woody savannas (WSA), grasslands (GRA), wetlands (WET), and croplands (CRO). sp, spring; su, summer; au, autumn, m, mean; max, maximal; s, standard deviation; and c, coefficient of variation. For example, mEVlsp is the mean enhanced vegetation index (EVI) during spring. The number after EVI or LST means the day of the year (DOY); for example, EVI243 is the EVI value at DOY234. The full names for independent variables can be found in Table S2. * [PERSON] et al. (2012) [PERSON]; [PERSON]; [PERSON] [PERSON]; [PERSON]; [PERSON]; [PERSON] Interannual variability of net carbon exchange is related to the lag between the end-dates of net carbon uptake and photosynthesis: Evidence from long records at two contrasting forest stands. _Agric. For. Meteorol._**2012**, _164_, 29-38. [CrossRef]
* [PERSON] et al. (2016) [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; et al. Matching the phenology of Net Ecosystem Exchange and vegetation indices estimated with MODIS and FLUXNET in-situ observations. _Remote Sens. Environ._**2016**, _174_, 290-300. [CrossRef]
* [PERSON] et al. (2017) [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; et al. Land surface phenology derived from normalized difference vegetation index (NDVI) at global FLUXNET sites. _Agric. For. Meteorol._**2017**, _233_, 171-182. [CrossRef]
* [PERSON] et al. (2013) [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON] Climate change, phenology, and phenological control of vegetation feedbacks to the climate system. _Agric. For. Meteorol._**2013**, _169_, 156-173. [CrossRef]
* [PERSON] et al. (2015) [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON] Radiation contributed more than temperature to increased decadal autumn and annual carbon uptake of two eastern North America mature forests. _Agric. For. Meteorol._**2015**, _201_, 740-749. [CrossRef]
* [PERSON] et al. (2017) [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON] [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON] [PERSON]; et al. Climate controls over the net carbon uptake period and amplitude of net ecosystem production in temperate and boreal ecosystems. _Agric. For. Meteorol._**2017**, _243_, 9-18. [CrossRef]
* [PERSON] et al. (2011) [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON] A comparison of multiple phenology data sources for estimating seasonal transitions in deciduous forest carbon exchange. _Agric. For. Meteorol._**2011**, _151_, 1741-1752. [CrossRef]
* [PERSON] and [PERSON] (2016) [PERSON]; [PERSON] Improving vegetation phenological parameterization of, a land surface model. _Biogeosc. Discuss._**2016**. preprint.