Three-Dimensional Modeling of Soil Organic Carbon Stocks in Forest Ecosystems of Northeastern China Under Future Climate Warming Scenarios

Understanding the detailed spatiotemporal variations in soil organic carbon (SOC) stocks is essential for assessing soil carbon sequestration potential. However, most existing studies predominantly focus on topsoil SOC stocks, leaving significant knowledge gaps regarding critical zones, depth-dependent variations, and key influencing factors associated with deeper SOC stock dynamics. This study adopted a comprehensive methodology that integrates random forest modeling, equal-area soil profile analysis, and space-for-time substitution to predict depth-specific SOC stock dynamics under climate warming in Northeast China’s forest ecosystems. By combining these techniques, the approach effectively addresses existing research limitations and provides robust projections of soil carbon changes across various depth intervals. The analysis utilized 63 comprehensive soil profiles and 12 environmental predictors encompassing climatic, topographic, biological, and soil property variables. The model’s predictive accuracy was assessed using 10-fold cross-validation with four evaluation metrics: MAE, RMSE, R2, and LCCC, ensuring comprehensive performance evaluation. Validation results demonstrated the model’s robust predictive capability across all soil layers, achieving high accuracy with minimized MAE and RMSE values while maintaining elevated R2 and LCCC scores. Three-dimensional spatial projections revealed distinct SOC distribution patterns, with higher stocks concentrated in central regions and lower stocks prevalent in northern areas. Under simulated warming conditions (1.5 °C, 2 °C, and 4 °C increases), both topsoil (0–30 cm) and deep-layer (100 cm) SOC stocks exhibited consistent declining trends, with the most pronounced reductions observed under the 4 °C warming scenario. Additionally, the study identified mean annual temperature (MAT) and normalized difference vegetation index (NDVI) as dominant environmental drivers controlling three-dimensional SOC spatial variability. These findings underscore the importance of depth-resolved SOC stock assessments and suggest that precise three-dimensional mapping of SOC distribution under various climate change projections can inform more effective land management strategies, ultimately enhancing regional soil carbon storage capacity in forest ecosystems.