An Interpretable Machine Learning Framework for Unraveling the Dynamics of Surface Soil Moisture Drivers

Understanding the impacts of the spatial non-stationarity of environmental factors on surface soil moisture (SSM) in different seasons is crucial for effective environmental management. Yet, our knowledge of this phenomenon remains limited. This study introduces an interpretable machine learning framework that combines the SHapley Additive exPlanations (SHAP) method with two-step clustering to unravel the spatial drivers of SSM across Iran. Due to the limited availability of in situ SSM data, the performance of three global SSM datasets—SMAP, MERRA-2, and CFSv2—from 2015 to 2023 was evaluated using agrometeorological stations. SMAP outperformed the others, showing the highest median correlation and the lowest Root Mean Square Error (RMSE). Using SMAP, we estimated SSM across 609 catchments employing the Random Forest (RF) algorithm. The RF model yielded R2 values of 0.89, 0.83, 0.70, and 0.75 for winter, spring, summer, and autumn, respectively, with corresponding RMSE values of 0.076, 0.081, 0.098, and 0.061 m3/m3. SHAP analysis revealed that climatic factors primarily drive SSM in winter and autumn, while vegetation and soil characteristics are more influential in spring and summer. The clustering results showed that Iran’s catchments can be grouped into five categories based on the SHAP method coefficients, highlighting regional differences in SSM controls.