Submarine landslides have attracted widespread attention, with the continuous development of ocean engineering. Due to the recent developments of in-situ investigation and modelling techniques of submarine landslides, significant improvements were achieved in the evolution studies on submarine landslides. The general characteristics of typical submarine landslides in the world are analyzed. Based on this, three stages of submarine landslide disaster evolution are proposed, namely, the submarine slope instability evolution stage, the large deformation landslide movement stage, and the stage of submarine landslide deposition. Given these three stages, the evolution process of submarine landslide disaster is revealed from the perspectives of in-situ investigation techniques, physical simulation, and numerical simulation methods, respectively. For long-term investigation of submarine landslides, an in-situ monitoring system with long-term service and multi-parameter collaborative observation deserves to be developed. The mechanism of submarine landslide evolution and the early warning factors need to be further studied by physical modelling experiments. The whole process of the numerical simulation of submarine landslides, from seabed instability to large deformation sliding to the impact on marine structures, and economizing the computational costs of models by advanced techniques such as parallel processing and GPU-accelerators, are the key development directions in numerical simulation. The current research deficiencies and future development directions in the subject of submarine landslides are proposed to provide a useful reference for the prediction and early warning of submarine landslide disasters.

The reliable numerical simulation of the landslide process contributes to the establishment of evidence-based disaster mitigation measures in seismically active zones. To achieve this goal, a simple and unified state-dependent shear strength model of discontinuities is presented to describe the shear strength degradation in a seismic landslide process. The proposed model establishes a relationship between the shear strength parameters and the global safety factor of the slope by assuming that the slope instability (or landslide initiation) is accompanied by an instantaneous shear strength degradation of discontinuities. To realize the model numerically, the algorithms for the computation of global safety factor and the modification of shear strength parameters were incorporated into the discontinuous deformation analysis (DDA). Subsequently, the kinematic accuracy of the improved DDA method was validated by comparisons with theoretical solutions for the dynamic sliding of a block on an inclined plane. Numerical simulations of the Daguangbao landslide triggered by the Wenchuan earthquake were performed using the improved DDA method. The results illustrate that the shear strength degradation of discontinuities affect the evolution process, travel distance, and post-failure shape of the seismic landslide significantly.