Spectral element modeling of three dimensional wave propagation in a self-gravitating Earth with an arbitrarily stratified outer core

Abstract

This paper deals with the spectral element modeling of seismic wave propagation at the global scale. Two aspects relevant to low-frequency studies are particularly emphasized. First, the method is generalized beyond the Cowling approximation in order to fully account for the effects of self-gravitation. In particular, the perturbation of the gravity field outside the Earth is handled by a projection of the spectral element solution onto the basis of spherical harmonics. Second, we propose a new formulation inside the fluid which allows to account for an arbitrary density stratification. It is based upon a decomposition of the displacement into two scalar potentials, and results in a fully explicit fluid-solid coupling strategy. The implementation of the method is carefully detailed and its accuracy is demonstrated through a series of benchmark tests.

Sent to Geophysical Journal International on July 29, 2003

Cited by in corpus (6)

• An axisymmetric time-domain spectral-element method for full-wave simulations: Application to ocean acoustics
• A Discontinuous Galerkin method with a modified penalty flux for the propagation and scattering of acousto-elastic waves
• A spectral-infinite-element solution of Poisson's equation: an application to self gravity
• A weight-adjusted discontinuous Galerkin method for wave propagation in coupled elastic-acoustic media
• Acoustic and inertial modes in planetary-like rotating ellipsoids
• Geophysical flows over topography, a playground for laboratory experiments