Abstract: Material Point Method(MPM)belongs to one of the meshless numerical analysis methods. It is able to integrate a variety of constitutive models and consider the fluid-solid coupling behavior within soils, so as to effectively simulate the large and extra-large deformation of landslides. This paper introduces the fundamental principles of MPM, the governing equations and the solving formats. Based on the 5 types of porous media characterization models of MPM, we briefly review the state-of-the-art of run-out process simulation of soil landslides using MPM.Subsequently, we simulate the run-out process of Shenzhen "12·20" landfill landslide using MPM with the single-layer and sing-phase material point model. Linear loading is applied to identify the initial stress field. Elasto-palstic constitutive models using Drucker-Prager yield criterion, GIMP algorithm and MUSL solving format are also used, respectively. The simulation results show that the maximal sliding distance reaches 510m, after the slope losing its stability. Moreover, the longitudinal distance of the landslide scope is 1050m and the minimal landslide angle is 5.95°, both agreeing well with the observed results. The evolution trends of the equivalent plastic strain zone within the soil exhibit that the landslide mechanism is progressive failure. Specifically, shear failure of the soil occurs at the slope toe initially, and tensile fracture occurs at the slope top subsequently, with both plastic zones extending toward the interior of the slope along the bedrock. The extra-large deformation occurs once the sliding surface is cut-through to accelerate the sliding. The MPM simulation for the run-out process of landslides can shed light on the hazardous behaviors of landslides and provide assistants for landslide prevention.