Landslides represent a large-deformation process influenced by various factors of uncertainty, such as types of ground motion (GM), randomness of GMs, and spatial variability of soils. The behavior of landslides is profoundly affected, making their assessment and measurement challenging. This paper proposes a computational approach for simulating the large-deformation process of landslides based on three-dimensional (3D) non-ordinary state-based peridynamics (NOSBPD). The analysis results of NOSBPD indicate that the mean runout distance triggered by pulse-like ground motions (PLGMs) is 16% greater than that induced by non-pulse ground motions (NPGMs), suggesting that PLGMs exhibit higher destructiveness. Besides, this paper establishes a runout distance assessment framework by considering the stochastic nature of PLGMs. This framework allows for the evaluation of the specific risk probability of landslides caused by stochastic PLGMs that match the target spectrum specified by codes. By introducing the theory of random fields and implementing a coupled procedure in peridynamics, we conducted a detailed analysis of the impact of spatial heterogeneity on the evolution process and consequences of landslides. Additionally, compared to two-dimensional (2D) analysis, the mean runout distance obtained from 3D analysis increased by 27.5%. This suggests that 2D analysis may underestimate the consequences of landslides. The findings of this study can serve as a scientific foundation for predicting the extent and scope of landslides triggered by near-fault earthquakes.