Abstract:The Huang-Huai-Hai Plain is an important agricultural production base in China. Accurately quantifying the carbon budget in agricultural production could provide a scientific basis for the agriculture development in China. The Biome-BGCMuSo model was used to simulate the gross primary productivity (GPP), ecosystem respiration (RECO), and net ecosystem productivity (NEP)in the winter wheat summer maize rotation systems of the Huang Huai Hai Plain in period of 1981—2100, and their spatiotemporal change trends were analyzed as well. The results showed that after parameter optimization through the differential evolution algorithms, the simulation accuracy of the Biome BGCMuSo model was high. The determination coefficients (R2 )between the simulated and observed daily values of GPP, RECO, and NEP were 0. 80, 0. 64, and 0. 68, respectively. And the R2 values were 0. 92, 0. 73, and 0. 83 at monthly scale, respectively. In historical period of 1981—2020, the average annual values of GPP, RECO, and NEP were 1 682. 40 gC / (m2 ·a ), 1 293. 12 gC / (m2 ·a ), and 389. 28 gC / (m2·a), respectively. Among them, GPP and RECO were increased significantly, but NEP was increased significantly before 2000 and then changed insignificantly. In future period of 2020—2100, the average annual values of GPP, RECO, and NEP were 2 015. 75 gC / (m2·a), 1 613. 06 gC / (m2·a), and 402. 69 gC / (m2 · a )under the SSP245 scenario, while they were 2 144. 67 gC / (m2 · a ), 1 729. 52 gC / (m2·a), and 415. 15 gC / (m2·a)under the SSP585 scenario, respectively. GPP and RECO showed obvious increasing trends in the future, while NEP showed a relatively stable trend first and then decreased significantly after 2051. At spatial scale, GPP and RECO showed a distribution pattern of high in the middle and low in the northern and southern parts of the Huang-Huai-Hai Plain, while NEP was increased from southeast to northwest. Under future climate change, the weakening carbon sink function of the winter wheat summer maize rotation system will be mainly because the increasing rate of ecosystem respiration exceeded the increasing rate of crop gross primary productivity in the future. Since the quantification of agricultural carbon budget was of great significance for the ‘ carbon peak and carbon neutrality ’ goals and food security of China, more studies were needed to further quantify the carbon budget and explain the underlying change mechanism in the Huang-Huai-Hai Plain in the future.