<p>The stability analysis of jointed rock slopes under seismic activity is a crucial and challenging task. To assess the various potential damage modes of rock slopes featuring tension fracture planes, it is necessary to take into consideration the nonlinear properties of shear strength and the dynamic features of seismic activity. A modified pseudo-dynamic technique has been used to represent the slope safety coefficient under seismic action by utilizing the nonlinear Barton-Bandis failure criterion. Furthermore, a system reliability analysis of slope stability is carried out using a Monte Carlo simulation method. The findings indicate that the slope safety coefficient exhibits periodic variations over time when considering seismic action in dynamic conditions. Additionally, the interaction force between the blocks significantly influences the factor of safety. The slope system stability is notably affected by the horizontal seismic coefficient, while the vertical seismic coefficient has a less effect. Among the various potential failure modes, failure modes 2 and 4 are identified as the primary failure modes. The likelihood of failure mode 1 occurrence gradually increases with an increase in the nodal strength parameter.</p>