<p>Based on the limit analysis theorem, the expression for the stability coefficient of multi-stage loess slopes considering the amplification effect of ground shaking is derived using the pseudo-dynamic method combined with the geometrical characteristics of the slopes. The cyclic dynamics of the slope stability and the scale of landslides are revealed. Through the parameter analysis of the three-stage loess slopes, the results show that the slope-stability coefficients in the seismic time interval obtained by the pseudo-dynamic method show simple harmonic variation and the calculation results are slightly higher than those of the pseudo-static method. The geometric parameters have a significant influence on the stability of multi-stage loess slopes. At the same height, slowing down the slope rate, increasing the platform width, or increasing the number of slope stages can improve the overall stability of the slopes; with the increase in seismic coefficients and magnification factors, the seismic performance of the slopes is significantly reduced.</p>