Abstract. Top-down inversions of China's terrestrial carbon sink are known to be uncertain because of errors related to the relatively coarse resolution of global transport models and the sparseness of in situ observations. Taking advantage of regional chemistry transport models for mesoscale simulation and spaceborne sensors for spatial coverage, the Greenhouse Gases Observing Satellite (GOSAT) retrievals of column-mean dry mole fraction of carbon dioxide (XCO2) were introduced in the Models-3 (a flexible software framework) Community Multi-scale Air Quality (CMAQ) and ensemble Kalman smoother (EnKS)-based regional inversion system to constrain China's biosphere sink at a spatiotemporal resolution of 64 km and 1 h. In general, the annual, monthly, and daily variation in biosphere flux was reliably delivered, attributable to the novel flux forecast model, reasonable CMAQ background simulation, well-designed observational operator, and Joint Data Assimilation Scheme (JDAS) of CO2 concentrations and natural fluxes. The size of the assimilated biosphere sink in China was −0.47 Pg C yr−1, which was comparable with most global estimates (i.e., −0.27 to −0.68 Pg C yr−1). Furthermore, the seasonal patterns were recalibrated well, with a growing season that shifted earlier in the year over central and south China. Moreover, the provincial-scale biosphere flux was re-estimated, and the difference between the a posteriori and a priori flux ranged from −7.03 Tg C yr−1 in Heilongjiang to 2.95 Tg C yr−1 in Shandong. Additionally, better performance of the a posteriori flux in contrast to the a priori flux was statistically detectable when the simulation was fitted to independent observations, indicating sufficient to robustly constrained state variables and improved fluxes estimation. This study serves as a basis for future fine-scale top-down carbon assimilation.