Landslides in sensitive clays are catastrophic events that threaten life and property. When plastic strain occurs in sensitive clays, there is a significant reduction in the undrained shear strength. Consequently, even minor triggers can lead to multiple progressive landslides, causing extensive devastation. Investigating landslides in sensitive clay presents challenges, necessitating robust numerical methods to comprehend failure modes and collapse behaviors, particularly in the post-failure phase. Although several numerical methods have been utilized for simulations, most have primarily focused on upward progressive landslides. In this study, we adopted the nodal integration-based particle finite element method to simulate downward progressive landslides triggered by surcharge loads from an embankment near the crest in sensitive clay on a long slope. We also considered nonlinear strain softening in sensitive clays. The effects of the strain softening rate and strength gradient of the sensitive clay layer on the failure mechanisms and destruction extent were investigated. The simulation results reveal that a downward progressive landslide initiates with a localized flow slide surrounding the embankment and usually involves a combination of various failure modes.