Nonlinear Threshold Responses and Spatial Heterogeneity of Soil Organic Carbon Under Contrasting Pedoclimatic Regimes

Jin Cui¹Yuzhi Xu¹Mengqi Wang¹Aiju Liu¹Lin Sun¹Xinyu Feng¹Qingrun Yang¹Shang Wang^2*Hongqiang Liu³Yujuan Lv⁴Kai Liu^1*

• ¹Shandong University of Technology, Zibo, China
• ²UFZ-Helmholtz Centre for Environmental Research, Halle, Germany
• ³Central South Academy of Inventory and Planning of National Forestry and Grassland Administration, Changsha, China
• ⁴MEE Nanjing Institute of Environmental Sciences, Nanjing, China

Soil organic carbon (SOC) exhibits distinct spatial heterogeneity across different pedoclimatic regions, yet the underlying regulatory mechanisms and their threshold responses remain poorly understood. In this study, the spatial patterns and underlying region-specific regulatory factors controlling SOC dynamics were investigated across a pedoclimatic gradient represented by the Jiaodong Peninsula (maritime monsoon climate) and Southwest Shandong (continental climate) in Shandong Province, China. Geostatistical analysis coupled with sequential Gaussian simulations provided probabilistic assessment of SOC spatial patterns, while machine learning algorithms (Linear Regression, Random Forest, XGBoost and Support Vector Machine) integrated with SHAP analysis enabled quantification of nonlinear threshold responses and identification of dominant factors governing SOC dynamics. The results showed that SOC in Jiaodong exhibited a west-high-east-low gradient (9.64-21.04 g·kg⁻¹) characterized by local-scale structure (range = 42.6 km, nugget = 39.14%), whereas Southwest Shandong showed higher SOC contents (11.32-27.15 g·kg⁻¹) dominated by macro-scale gradients (range = 68 km, nugget = 1.12%). The Random Forest model, which demonstrated robust predictive capacity (Adj R² > 0.98, RMSE < 0.021 g·kg⁻¹ across both regions), identified distinct regulatory mechanisms in Jiaodong, where NO₃⁻-N and extractable Fe exhibited a dual-threshold domain (NO₃⁻-N = 10.0 mg·kg⁻¹, Fe = 12.0 mg·kg⁻¹), with the marginal effect of Fe on SOC shifting from negative (SHAP = -0.82) to positive when NO₃⁻-N exceeded its threshold concentration. In Southwest Shandong, total nitrogen (TN) was revealed as the dominant predictor (67.12% of model variance), with a critical threshold at 3.25 g·kg⁻¹ above which SOC increased by 2.0 g·kg⁻¹, while NO₃⁻-N showed negative effects above 27 mg·kg⁻¹. This study demonstrates that the combination of interpretable machine learning and geostatistical approaches can effectively elucidate region-specific threshold mechanisms and nonlinear controls governing SOC dynamics. This approach is critical for developing spatially-explicit soil carbon management strategies under varying pedoclimatic conditions.