Abstract: Estimating the spatiotemporal evolution and carbon sequestration potential of terrestrial vegetation under different emission scenarios in Heilongjiang Province from 2021 to 2060 is of great significance to achieving the "dual carbon" goals.Based on the Coupled Model Intercomparison Project Phase 6 (CMIP6) multi-model ensemble simulation data,this study drove the carbon cycle model (BEPS) to analyze the vegetation productivity of terrestrial ecosystems under three emission scenarios (SSP126,SSP245 and SSP585) in Heilongjiang Province from 2021 to 2060.Quantitative estimations of gross primary productivity (GPP) and net primary productivity (NPP) in Heilongjiang Province are generated at a spatial resolution of 1 km.The results showed that the corrected BEPS model driving data exhibited improved accuracy in Heilongjiang Province from 2021 to 2060,thereby enhancing the accuracy of vegetation productivity prediction.Under three emission scenarios,the vegetation NPP in most regions of Heilongjiang Province was higher than 500 gC·m^-2·a^-1 from 2021 to 2060.Compared with other regions,the forest areas in the central and southern part of Heilongjiang Province had the strongest vegetation productivity,and the vegetation NPP remaining above 700 gC·m^-2·a^-1.Under different emission scenarios,the overall vegetation productivity of Heilongjiang Province increased from 2021 to 2060,among which the Daxing'anling forest area and the eastern agricultural region showed the greatest increase.Compared with the SSP126 and SSP245 scenarios,the SSP585 scenario had the highest vegetation productivity in Heilongjiang Province,with the greatest increase and the largest interannual variability.