This study developed a discrete element simulation model based on the 2D particle flow code (PFC2D). of which the mesoscopic parameters were calibrated by the indoor experiments. to investigate the rainfall erosion damage of residual soil slope in the intermittently frozen area. It is to be noted that the runoff scouring action was simulated according to the equivalent rainfall method. https://www.drcristinosuarez.com/