Objectives This study aimed to investigate calibration methods for improving the accuracy of the RothC model in paddy fields of Northeast China, and to evaluate the spatiotemporal dynamics of soil organic carbon (SOC) under different straw return scenarios in typical rice-planting areas of Northeast China, thereby providing a basis for enhancing the soil fertility improvement effect of straw return.

Methods The study area was the 850 Farm located in Qiqihar City, Heilongjiang Province. Soil organic matter data for this farm were obtained from HWSD2, land use data from the 2020 CLUD dataset, and soil type information from the 1:4000000 Soil Map of China compiled by the Institute of Soil Science, Chinese Academy of Sciences. The RothC model was used to predict the dynamic changes of the soil organic carbon pool in paddy fields at both the experimental site and regional scales. The site-scale model was calibrated and validated using measured data collected from the field and long-term positioning observation data. Finally, by integrating rice planting area, soil properties, and meteorological data of Qiqihar City, the spatiotemporal changes of SOC under different straw return scenarios at the regional scale were simulated.

Results The tested models included the RothC model and its two modified versions, RothC_p and RothC_0.6. Each model was calibrated using three methods (M1, M2, M3), resulting in nine models. All models and calibration methods overestimated the observed trends to varying degrees. The RothC_0.6 model, which represented the anaerobic environment of paddy fields by fixing the moisture factor, aligned more closely with the actual straw decomposition processes in cold-region paddy fields. Consequently, it performed best in simulating local paddy soil organic carbon dynamics, especially under the treatment with only inorganic fertilizer and no straw return (NPK), where simulation accuracy was extremely high (normalized root mean square error, nRMSE<1%). For the treatment with inorganic fertilizer plus straw return (NPKS), the initial simulations overestimated SOC (nRMSE: 11.21%−12.12%). However, after semi-quantitatively adjusting the carbon input based on actual field management practices, model performance was substantially improved (nRMSE decreased to 1.42%−2.34%). Therefore, the RothC_0.6 model calibrated by the transfer function method was selected for subsequent regional-scale simulations. Twenty-year regional scenario simulations in Qiqihar showed that remaining root stubble (CK), with limited carbon input, could not compensate for SOC mineralization losses, ‌leading to a dramatic decrease of 10.49 Tg in SOC stock‌ relative to the baseline. Remaining root stubbles and shallow rotary tillage in Spring (CQX) sustained SOC stock with only a 0.7 Tg reduction, as carbon input roughly balanced decomposition. In contrast, full straw return significantly increased SOC stock. Among the straw return practices, all stubbles were cut and buried into soil and watering in Autumn (QSQM), with its semi-quantified carbon input closer to the actual amount and a more favorable decomposition environment during the fallow period, demonstrated significantly better carbon sequestration than all the straw and root stubbles remained in field and deep ploughing in Autumn (QSF), increasing SOC stock by 3.93 Tg (QSF) and 13.83 Tg (QSQM) relative to the baseline. The spatial distribution of SOC density in Qiqihar paddy fields was generally lower in the southwest and higher in the northeast, and straw return further enhanced the spatial heterogeneity of SOC density.

Conclusions  After being calibrated by the transfer function method (M3), the RothC_0.6 model exhibits high accuracy in simulating soil organic carbon (SOC) dynamics in single-cropping paddy fields of the cold region in Northeast China, and is suitable for regional-scale SOC variation research. The simulation results based on the calibrated model show that the SOC density of paddy fields in Qiqihar generally presents a spatial distribution pattern of “lower in the southwest and higher in the northeast”. Different straw return practices have significant effects on SOC stock. Specifically, remaining root stubble (CK) and remaining root stubbles and shallow rotary tillage in Spring (CQX) are insufficient to maintain the existing SOC stock, while full straw return significantly increases SOC stock by enhancing the absolute content of the humus carbon pool. Further regional-scale simulation indicates that the paddy fields in the northwest of Qiqihar have higher carbon sequestration potential, and straw return enhances the spatial heterogeneity of SOC distribution. These results can provide scientific basis and decision support for straw resource utilization management and black soil carbon sequestration technology optimization in cold-region paddy fields of Northeast China.