Wave propagation in spherical sections by the Chebyshev spectral method
Igel, H. (1999),
Wave propagation in spherical sections by the Chebyshev spectral method,
Geophys. J. Int., 136(3), 559566.
 Abstract
 The elastic wave equation in spherical coordinates is solved by a Chebyshev spectral method. In the algorithm presented the singularities in the governing equations are avoided by centring the physical domain around the equator. The highly accurate pseudospectral (PS) derivative operators reduce the required grid size compared to finite difference (FD) algorithms. The nonstaggered grid scheme allows easy extension to general material anisotropy without additional interpolations being required as in staggered FD schemes. The boundary conditions previously derived for curvilinear coordinate systems can be applied directly to the velocity vector and stress tensor in the spherical basis. The algorithm is applied to the problem of a doublecouple source located in a highvelocity region at the top of the mantle (slab). The synthetic seismograms show azimuthdependent traveltime and waveform effects which are likely to be observable in regions where subduction takes place. Such techniques are important in modelling the fullwave characteristics of the Earth's 3D structure and in providing accurate reference solutions for 3D global models.
 BibTeX

@article{id323,
author = {H. Igel},
journal = {Geophys. J. Int.},
number = {3},
pages = {559566},
title = {{Wave propagation in spherical sections by the Chebyshev spectral method}},
volume = {136},
year = {1999},
language = {en},
}
 EndNote

%0 Journal Article
%A Igel, H.
%D 1999
%N 3
%V 136
%J Geophys. J. Int.
%P 559566
%T Wave propagation in spherical sections by the Chebyshev spectral method