EARTHQUAKE GROUND MOTION SIMULATION THROUGH THE 2-D SPECTRAL ELEMENT METHOD

Abstract

The application of the 2-D Chebyshev spectral element method (SPEM) to engineering seismology problems is reviewed in this paper. The SPEM is a high-order finite element technique which solves the variational formulation of the seismic wave propagation equations. The computational domain is discretised into an unstructured grid composed by irregular quadrilateral elements. This property makes the SPEM particularly suitable to compute numerically accurate solutions of the full wave equations in complex media. The earthquake is simulated following an approach that can be considered "global", that is all the factors influencing the wave propagation — source, crustal heterogeneity, fine details of the near-surface structure, and topography — are taken into account and solved simultaneously. The basic earthquake source is represented by a 2-D point double couple model. Ruptures propagating along fault segments placed on the model plane are simulated as a finite summation of elementary point sources. After a general introduction, the paper first gives an overview of the method; then it concentrates on some methodological topics of interest for practical applications, such as quadrangular mesh generation, source definition and scaling, numerical accuracy and computational efficiency. Limitations and advantages of using a 2-D approach, although sophisticated such as the SPEM, are addressed, as well. The effectiveness of the method is illustrated through two case histories, i.e. the ground shaking prediction in Catania (Sicily, Italy) for a catastrophic earthquake, and the analysis of the ground motion in the presence of a massive structure.