Abstract: An efficient model reduction-based topology optimization method is developed for solving the design problem of continuous structures under seismic acceleration. A topology optimization model of minimizing the maximum dynamic compliance under material volume constraint is established. To ensure the differentiability of the objective function, the generalized mean function is used to represent the maximum dynamic compliance in the time domain. Topology optimization is carried out using the enhanced parametric level set band method, which allows gradient-based optimization algorithms to solve optimization problems. The derivatives of maximum dynamic compliance and volume constraint with respect to design variables are detailed by combining the direct method and adjoint vector method. Subsequently, the method of moving asymptotes (MMA) is employed to update design variables. To enhance the efficiency of optimization calculation, a model reduction technique based on Quasi-Static Ritz Vector (QSRV) is incorporated into the Newmark method to achieve efficient solution of transient dynamics equations and adjoint equations in sensitivity analysis. Finally, three typical examples under artificial earthquake are given to illustrate the effectiveness of the proposed method. At the same time, the optimization results based on the QSRV reduction method and the Newmark total analysis method are compared. The results indicate that the topology optimization method based on QSRV reduction technology can significantly reduce the calculation time of transient dynamic optimization iteration and improve the calculation efficiency.