Spectral-element simulations of global seismic wave propagation: II. Three-dimensional models, oceans, rotation and self-gravitation

SUMMARY We simulate global seismic wave propagation based upon a spectral-element method. We include the full complexity of 3-D Earth models, i.e. lateral variations in compressional-wave velocity, shear-wave velocity and density, a 3-D crustal model, ellipticity, as well as topography and bathymetry. We also include the effects of the oceans, rotation and self-gravitation in the context of the Cowling approximation. For the oceans we introduce a formulation based upon an equivalent load in which the oceans do not need to be meshed explicitly. Some of these effects, which are often considered negligible in global seismology, can in fact play a significant role for certain source‐receiver configurations. Anisotropy and attenuation, which were introduced and validated in a previous paper, are also incorporated in this study. The complex phenomena that are taken into account are introduced in such a way that we preserve the main advantages of the spectral-element method, which are an exactly diagonal mass matrix and very high computational efficiency on parallel computers. For self-gravitation and the oceans we benchmark spectral-element synthetic seismograms against normal-mode synthetics for the spherically symmetric reference model PREM. The two methods are in excellent agreement for all body- and surface-wave arrivals with periods greater than about 20 s in the case of self-gravitation and 25 s in the case of the oceans. At long periods the effect of gravity on multiorbit surface waves up to R4 is correctly reproduced. We subsequently present results of simulations for two real earthquakes in fully 3-D Earth models for which the fit to the data is significantly improved compared with classical normal-mode calculations based upon PREM. For example, we show that for trans-Pacific paths the Rayleigh wave can arrive more than a minute earlier than in PREM, and that the Love wave is much shorter in duration.

[1]  Jean-Pierre Vilotte,et al.  Solving elastodynamics in a fluid-solid heterogeneous sphere: a parallel spectral element approximation on non-conforming grids , 2003 .

[2]  D. Komatitsch,et al.  Spectral-element simulations of global seismic wave propagation—I. Validation , 2002 .

[3]  D. Komatitsch,et al.  Modeling of Seismic Wave Propagation at the Scale of the Earth on a Large Beowulf , 2001, ACM/IEEE SC 2001 Conference (SC'01).

[4]  Jack Dongarra,et al.  Top500 Supercomputer Sites , 1997 .

[5]  J. Woodhouse,et al.  Complex Shear Wave Velocity Structure Imaged Beneath Africa and Iceland. , 1999, Science.

[6]  J. Tromp,et al.  Theoretical Global Seismology , 1998 .

[7]  Gabi Laske,et al.  CRUST 5.1: A global crustal model at 5° × 5° , 1998 .

[8]  D. Kopriva.  A Conservative Staggered-Grid Chebyshev Multidomain Method for Compressible Flows. II. A Semi-Structured Method , 1996 .

[9]  P. Paolucci,et al.  The “Cubed Sphere” , 1996 .

[10]  Einar M. Rønquist.  Convection treatment using spectral elements of different order , 1996 .

[11]  Thomas Sterling,et al.  Enabling Technologies for Petaflops Computing , 1995 .

[12]  Bernard Valette,et al.  About the influence of pre‐stress upon adiabatic perturbations of the Earth , 1986 .

[13]  J. A. Deruntz,et al.  Finite Element Analysis of Shock-Induced Hull Cavitation , 1984 .

[14]  D. L. Anderson,et al.  Preliminary reference earth model , 1981 .

[15]  K. C. Park,et al.  Partitioned Transient Analysis Procedures for Coupled-Field Problems. , 1979 .

[16]  R. Sadourny.  Conservative Finite-Difference Approximations of the Primitive Equations on Quasi-Uniform Spherical Grids , 1972 .

[17]  T. G. Cowling,et al.  The non-radial oscillations of polytropic stars , 1941 .

[18]  Jean-Pierre Vilotte,et al.  Coupling the spectral element method with a modal solution for elastic wave propagation in global earth models , 2003 .

[19]  C. Bassin,et al.  The Current Limits of resolution for surface wave tomography in North America , 2000 .

[20]  Thomas Sterling,et al.  How to Build a Beowulf: A Guide to the Implementation and Application of PC Clusters 2nd Printing , 1999 .

[21]  M. Taylor.  The Spectral Element Method for the Shallow Water Equations on the Sphere , 1997 .

[22]  C. Bernardi,et al.  A New Nonconforming Approach to Domain Decomposition : The Mortar Element Method , 1994 .

[23]  Yvon Maday,et al.  Coupling finite element and spectral methods: first results , 1990 .

[24]  T. A. Zang,et al.  Spectral methods for fluid dynamics , 1987 .