The inversion of Megaseism M 712 in coastal area of Quanzhou in 1604 is performed by means of deterministic parsing of seismic hazards

to attain the matrix horizontal peak acceleration （191gal） as well as the response spectrum of matrix horizontal acceleration at the pagoda site caused by the earthquake. Then combined with the site geotechnical characteristic determined in the test

one-dimensional equivalent linearized fluctuation is applied to calculate the earthquake responses of the soil layer in order to attain the horizontal acceleration time-history in the site ground. Based on the MATLAB platform together with the findings of seismic hazard analysis

the program Wilson θ-method is used to conduct the elastic time-history dynamic analysis of structures of East and West Towers. According to the destruction conditions of East and West Towers in the earthquake from historic records

the limited value 1/532 of the elastic interlayer displacement in East and West Towers is obtained. There is no great difference between the analytical results in this paper and those of Reference 4

but the former seems to be more theoretical reasonable. This analysis can be used as a reference for investigation of the aseismic deformability of ancient masonry pagodas and prediction of seismic damages.