This study examined long-term changes in soil carbon stock dynamics 11 and 19 years after fire under different severities at 0–5 and 0–25 cm depths with a digital soil mapping approach. Linear (MLR) and non-linear models (RF, SVR, XGBoost) combined with feature selection methods (r < 0.8, FFS, Boruta) were used to predict bulk density (BD), total C, and C stock. Distributional biases were evaluated with Kolmogorov–Smirnov statistics and corrected by Quantile Mapping (QM). RF-FFS performed best for BD and total C at 0–5, while RF-SVR outperformed for C stock and all properties at 0–25. Total C was 49% higher at 0–5, whereas C stock was 7.57 times greater at 0–25. Both models underestimated variability, especially for C stock. At 0–25, bulk density decreased after fire, particularly under conditions of medium severity, while total C increased following the same tendency. The results showed that fire’s legacy is still present in the ecosystem after one and two decades. This is particularly evident at greater depths, where long-term C stock is lower.
Benhalima et al. (Tue,) studied this question.