Coupling the spectral element method with a modal solution for elastic wave propagation in global earth models

Abstract

International audienceWe present a new method for wave propagation in global earth models based upon the coupling between the spectral element method and a modal solution method. The Earth is decomposed into two parts, an outer shell with 3-D lateral heterogeneities and an inner sphere with only spherically symmetric heterogeneities. Depending on the problem, the outer heterogeneous shell can be mapped as the whole mantle or restricted only to the upper mantle or the crust. In the outer shell, the solution is sought in terms of the spectral element method, which stem from a high order variational formulation in space and a second-order explicit scheme in time. In the inner sphere, the solution is sought in terms of a modal solution in frequency after expansion on the spherical harmonics basis. The spectral element method combines the geometrical flexibility of finite element methods with the exponential convergence rate of spectral methods. It avoids the pole problems and allows for local mesh refinement, using a non-conforming discretization, for the resolution of sharp variations and topography along interfaces. The modal solution allows for an accurate isotropic representation in the inner sphere. The coupling is introduced within the spectral element method via a Dirichlet-to-Neumann (DtN) operator. The operator is explicitly constructed in frequency and in generalized spherical harmonics. The inverse transform in space and time requires special attention and an asymptotic regularization. The coupled method allows a significant speed-up in the simulation of the wave propagation in earth models. For spherically symmetric earth model, the method is shown to have the accuracy of spectral transform methods and allow the resolution of wavefield propagation, in 3-D laterally heterogeneous models, without any perturbation hypothesis

An efficient coupled acoustic-elastic local solver applied to phase inversion

“…Only in the relatively small region of in-terest, they are coupled with a 3D heterogeneous solver. Capdeville et al (2003) use a similar approach, where the earth is again separated into two parts. This approach again limits the expensive solver to a small region, using a simplified solver for the remainder.…”

Regional wave propagation using the discontinuous Galerkin method

“…A particular FE scheme is the spectral-element (SE) method, which is state-of-the-art in regional and global earthquake simulations at the moment (e.g., Priolo et al, 1994;Komatitsch, 1997;Seriani, 1998;Chaljub, 2000;Capdeville et al, 2003;Nissen-Meyer et al, 2008;Fichtner et al, 2009;Cupillard et al, 2012). The SE method supports computational grids of deformable hexahedral elements, to simulate fully 3-D wavefields with high-order spatial accuracy, and explicit time extrapolation schemes.…”

Spectral Element Simulation of Rupture Dynamics

“…Spectral element method (SEM) locates at the center of this framework, combining the geometrical flexibility of finite elements with the convergence rate of the spectral methods and it already demonstrated a powerful tool in modeling the seismic wave propagation at regional and global scales (Komatitcsh and Vilotte [5], Capdeville et al [6]). Taking advantage from the variational approximation, the rupture modeling along non planar flaults through friction laws can be introduced in SEM as an interface condition, in constrast with finite differences, which generally deal with velocity-stress staggered grids, where mixed boundary conditions require necessarily some interpolation.…”