Ecosystems play an important role in global carbon cycle and climate change (Nemani et al., 2003;S. Q. Wang et al., 2015). Terrestrial gross primary productivity (GPP) is the primary driver of the global carbon cycle (Yuan et al., 2014;Zhang, Song, et al., 2016). Moderate fluctuations in terrestrial production can result in substantial variabilities for global warming (Cai et al., 2014). Accurate quantification of GPP can promote the understanding of the feedback between terrestrial biosphere and climate systems (Xiao et al., 2004).As a traditional method, eddy covariance technique derives site-level GPP as the difference between measured net ecosystem exchange (NEE) and estimated respiration (Reichstein et al., 2005). These GPP measurements are useful for exploring regional carbon budget and validating vegetation models (Mercado et al., 2009;Moore et al., 2018). However, the limited spatiotemporal coverage of site-level observations restricts the assessment of carbon fluxes at large scales and/or over long-term periods (Ran et al., 2016). Satellite remote sensing provides continuous retrieval of ecosystem productivity on a much wider scale than the flux tower data (Running et al., 2000). However, such products are not direct measurements, and as a result may introduce biases in the estimate of global GPP (Li & Xiao, 2019). Therefore, dynamic global vegetation models (DGVMs) have been developed for estimating regional to global carbon budget, in combination with site-level and remote sensing data (He et al., 2013;.The DGVMs help estimate GPP through the modeling of complex biophysical and biogeochemical processes but with uncertainties from multiple sources, including parameterization schemes (Zaehle et al., 2005), model framework (Song & Zeng, 2014), and/or input data (Wu et al., 2017). The first two types of uncertainties are related to model inherent characteristics, and can be improved step by step with the developed understanding of physical processes or better data constraints (Dietze et al., 2014;Raczka et al., 2018). The latter one, however, is related to the external forcing and can be improved independent of model development. Most of DGVMs are forced with climate reanalyzes, which were derived by assimilating numerical