A Lagrangian numerical approach for the simulation of rapid landslide runouts is presented and discussed. The simulation approach is based on the so-called Particle Finite Element Method. The moving soil mass is assumed to obey a rigid-viscoplastic, non-dilatant Drucker-Prager constitutive law, which is cast in the form of a regularized, pressure-sensitive Bingham model. Unlike in classical formulations of computational fluid mechanics, where no-slip boundary conditions are assumed, basal slip boundary conditions are introduced to account for the specific nature of the landslide-basal surface interface. The basal slip conditions are formulated in the form of modified Navier boundary conditions, with a pressure-sensitive threshold. A special mixed Eulerian-Lagrangian formulation is used for the elements on the basal interface to accommodate the new slip conditions into the Particle Finite Element Method framework. To avoid inconsistencies in the presence of complex shapes of the basal surface, the no-flux condition through the basal surface is relaxed using a penalty approach. The proposed model is validated by simulating both laboratory tests and a real large-scale problem, and the critical role of the basal slip is elucidated. Copyright (c) 2016 John Wiley & Sons, Ltd.