Abstract:Soil organic carbon (SOC) plays a crucial role in the global carbon cycle, and with the influence of global climate change and human activities, soil organic carbon density is constantly changing. A soil organic carbon density (SOCD) estimation method was proposed based on climate regionalization and a random forest model. It also developed a SOCD product with a long time series from the 1980s to 2020s and a spatial resolution of 1 km. The spatial heterogeneity and evolution patterns of SOCD in China from the 1980s to the 2020s were analyzed. Using Landsat series satellite images, elevation data, meteorological data, and measured SOCD data, a digital soil mapping method based on the random forest model was constructed to estimate the spatio-temporal distribution of 0~20 cm surface SOCD in China. The results showed that the prediction accuracy of the model considering climate zoning (R2=0.55, RMSE was 2.19 kg/m2) was better than that of the global model (R2=0.46, RMSE was 2.36 kg/m2). Meteorological factors had a significant impact on SOCD. Increasing temperature would accelerate the metabolic rate of microorganisms, promote the decomposition of soil organic matter, and lead to the increase of soil organic carbon release. Precipitation had a direct effect on soil water status, and suitable soil water content was conducive to SOC accumulation. At the same time, through verification with the measured data of the Heihe River basin, a high consistency was achieved between the model estimation results and the measured data (R2=0.69, RMSE was 2.01 kg/m2). The research result can provide a scientific basis for the accurate estimation and analysis of SOCD in China and it had important guiding significance for optimizing agricultural practice, improving soil carbon sink function, and realizing the national “double carbon” goal, which was conducive to promoting sustainable agricultural development and ecological environmental protection.