AbstractThe frequency of snowmelt‐induced soil slope instabilities is increasing in some seasonally cold regions because of climate change. Reliable hazard assessment and risk mitigation of snowmelt‐induced landslides require physically‐based prediction models. However, existing models either apply only at the slope scale or assume precipitation as the sole landslide trigger. In doing so, they neglect the complexity and coupled nature of the thermo‐hydro‐mechanical processes leading to slope instability in seasonally cold regions (such as snow accumulation and melting, infiltration and surface runoff, soil saturation, pore water pressure buildup and dissipation). Here, we present a spatially distributed and sequentially coupled numerical model to simulate snowmelt‐induced slope instabilities at the catchment scale. The model accounts for temperature‐dependent changes in the soil hydraulic behavior related to changes in water state by means of a routine implemented in a geographic information system. We verified the performance of the model using a case study of spring snowmelt‐induced soil slope failures that occurred after the 2004 Mid‐Niigata earthquake in Japan. Considering limitations and simplifications, the model was able to predict the triggering condition, magnitude, and spatial distribution of the snowmelt‐induced landslides with a satisfactory degree of accuracy. We believe that the robustness and simplicity of our numerical approach make it suitable for implementation in early warning systems.