An infinite slope stability numerical model driven by the comprehensive physically-based integrated hydrology model(InHM) is presented.In this approach,the failure plane is assumed to be parallel to the hydraulic gradient instead of the slope surface.The method helps with irregularities in complex terrain since depressions and flat areas are allowed in the model.The present model has been tested for two synthetic single slopes and a small catchment in the Mettman Ridge study area in Oregon,United States,to estimate the shallow landslide susceptibility.The results show that the present approach can reduce the simulation error of hydrological factors caused by the rolling topography and depressions,and is capable of estimating spatial-temporal variations for landslide susceptibilities at simple slopes as well as at catchment scale,providing a valuable tool for the prediction of shallow landslides. An infinite slope stability numerical model driven by the comprehensive physically-based integrated hydrology model (InHM) is presented. In this approach, the failure plane is assumed to be parallel to the hydraulic gradient instead of the slope surface. The method helps with irregularities in complex terrain since depressions and flat areas are allowed in the model the present model has been tested for two synthetic single slopes and a small catchment in the Mettman Ridge study area in Oregon, United States, to estimate the shallow landslide susceptibility.The results show that the present approach can reduce the simulation error of hydrological factors caused by the rolling topography and depressions, and is capable of estimating spatial-temporal variations for landslide susceptibilities at simple slopes as well as at catchment scale, providing a valuable tool for the prediction of shallow landslides