This paper introduces an improved modal pushover analysis (IMPA) which can effectively evaluate the seismic response of multi-span continuous bridge structures on the basis of modal pushover analysis (MPA). Differently from previous modal pushover analyses which cause the numerical unstability because of the occurrence of reversed relation between the pushover load and displacement, the proposed method eliminates this numerical instability and, in advance the coupling effects induced from the direct application of modal decomposition by introducing an identical stiffness ratio for each dynamic mode at the post-yielding stage together with an approximate elastic deformation. In addition to these two introductions, the use of an effective seismic load, calculated from the modal spatial force and applied as the distributed load, makes it possible to predict the dynamic responses of all bridge structures through a simpler analysis procedure than those in conventional modal pushover analyses. Finally, in order to establish validity and applicability of the proposed method, correlation studies between a rigorous nonlinear time history analysis and the proposed method were conducted for multi-span continuous bridges.