Interpretable machine learning guided by physical mechanisms reveals drivers of runoff under dynamic land use changes

Human activities continuously impact water balances and cycling in watersheds, making it essential to accurately identify the responses of runoff to dynamic changes in land use types. Although machine learning models demonstrate promise in capturing the intricate interplay between hydrological factors, their "black box" nature makes it challenging to identify the dynamic drivers of runoff. To overcome this challenge, we employed an interpretable machine learning method to inversely deduce the dynamic determinants within hydrological processes. In this study, we analyzed land use changes in the Ningxia section of the middle Yellow River across four periods, laying the foundation for revealing how these changes affect runoff. The sub-watershed attributes and meteorological characteristics generated by the Soil and Water Assessment Tool (SWAT) model were used as input variables of the Extreme Gradient Boosting (XGBoost) model to simulate substantial sub-watershed rainfall runoff in the region. The XGBoost was interpreted using the SHapley Additive exPlanations (SHAP) to identify the dynamic responses of runoff to the land use changes over different periods. The results revealed increasingly frequent interchanges between the land use types in the study area. The XGBoost effectively captured the characteristics of the hydrological processes in the SWAT-derived sub-watersheds. The SHAP analysis results demonstrated that the promoting effect of agricultural land (AGRL) on runoff gradually weakens, while forests (FRST) continuously strengthen their restraining effect on runoff. Relevant land use policies provide empirical support for these findings. Furthermore, the interaction between meteorological variables and land use impacts the runoff generation mechanism and exhibits a threshold effect, with the thresholds for relative humidity (RH), maximum temperature (MaxT), and minimum temperature (MinT) determined to be 0.8, 25 °C, and 15 °C, respectively. This reverse deduction method can reveal hydrological patterns and the mechanisms of interaction between variables, helping to effectively addressing constantly changing human activities and meteorological conditions.