Infiltration and groundwater have been widely considered as the main factors that cause shallow landslides; however, the effect of runoff has received less attention. In this study, an in-house physical-process-based shallow landslide model is developed to demonstrate the influence of runoff. The model is controlled by coupling the shallow water equation (dynamic) and Richards' equation. An infinite slope stability analysis is applied to evaluate the possibility of regional landslides. A real, small catchment topography is adopted as a demonstration example. The simulation illustrates the variations of runoff and the factor of safety (FS) during a storm. The results indicate that, after the surface becomes saturated, the FS may keep varying due to the increasing pressure head, which is caused by increasing surface water depth. This phenomenon most likely occurs downstream where the slopes easily accumulate water. The depth of the surface water may also be a factor of slope failure. Therefore, it is essential to increase the accuracy of calculating the runoff depth when assessing regional shallow landslides. stimulated by the rise of the groundwater table. Tsai and Wang [10] demonstrated the effect of different patterns of rainfall on shallow landslides in unsaturated soils, stimulated by dissipation through matric suction. Chen et al. [11] further explored the effects of rainfall duration, rainfall amount, and lateral flow-induced slope failures, using a vertical two-dimensional (2D) numerical landslide model. The relationships between rainfall infiltration and slope failures have been discussed extensively.Recently, some studies have connected landslides with runoff depth. For example, Chan et al.[12] combined a hydrological model and a landslide susceptibility model to establish a landslide analysis procedure. The results of their study indicated that the runoff flow depth may be selected as an analysis factor instead of the rainfall depth or maximum rainfall intensity. However, studies have seldom directly confirmed the physical mechanism of runoff on shallow landslides; although surface runoff is considered to increase the water pressure and affect the boundary conditions of infiltration, this behavior has not been demonstrated clearly.Furthermore, most landslide models were simplified for modeling the surface water flow. Some models focusing on the simulation of the infiltration mechanism did not consider the water depth [13,14]. Some assumed a kinematic or diffusion wave to model surface water flow [15][16][17]. However, ignoring the runoff or simplifying the momentum equations may fail to evaluate the effect of surface runoff. For the purpose of accurately evaluating the influence of runoff, dynamic equations may be more suitable.This study develops an in-house physical-process-based landslide model to observe the influence of runoff on shallow landslides in unsaturated soil. The model is governed by the coupling of the 2D shallow water equation (dynamic) and one-dimensional (1D) Richards' equation to simu...