Annual gross primary productivity (AGPP) is the basis for grain production and terrestrial carbon sequestration. Mapping regional AGPP from site measurements provides methodological support for analysing AGPP spatiotemporal variations thereby ensures regional food security and mitigates climate change. Based on 641 site-year eddy covariance measuring AGPP from China, we built an AGPP mapping scheme based on its formation and selected the optimal mapping way, which was conducted through analysing the predicting performances of divergent mapping tools, variable combinations, and mapping approaches in predicting observed AGPP variations. The reasonability of the selected optimal scheme was confirmed by assessing the consistency between its generating AGPP and previous products in spatiotemporal variations and total amount. Random forest regression tree explained 85 % of observed AGPP variations, outperforming other machine learning algorithms and classical statistical methods. Variable combinations containing climate, soil, and biological factors showed superior performance to other variable combinations. Mapping AGPP through predicting AGPP per leaf area (PAGPP) explained 86 % of AGPP variations, which was superior to other approaches. The optimal scheme was thus using a random forest regression tree, combining climate, soil, and biological variables, and predicting PAGPP. The optimal scheme generating AGPP of Chinese terrestrial ecosystems decreased from southeast to northwest, which was highly consistent with previous products. The interannual trend and interannual variation of our generating AGPP showed a decreasing trend from east to west and from southeast to northwest, respectively, which was consistent with data-oriented products. The mean total amount of generated AGPP was 7.03 ± 0.45 PgC yr
