Open AccessGaussian beams for surface waves in laterally slowly‐varying media
Author(s) -
Yomogida Kiyoshi
Publication year - 1985
Publication title -
geophysical journal of the royal astronomical society
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0016-8009
DOI - 10.1111/j.1365-246x.1985.tb05149.x
Subject(s) - eigenfunction , ray tracing (physics) , gaussian curvature , gaussian , rayleigh wave , surface wave , perpendicular , wave propagation , physics , surface (topology) , geometry , internal wave , curvature , mathematical analysis , classical mechanics , eigenvalues and eigenvectors , optics , mathematics , quantum mechanics
Summary. Asymptotic ray theory is applied to surface waves in a medium where the lateral variations of structure are very smooth. Using ray‐centred coordinates, parabolic equations are obtained for lateral variations while vertical structural variations at a given point are specified by eigenfunctions of normal mode theory as for the laterally homogeneous case. Final results on wavefields close to a ray can be expressed by formulations similar to those for elastic body waves in 2‐D laterally heterogeneous media, except that the vertical dependence is described by eigenfunctions of ‘local’ Love or Rayleigh waves. The transport equation is written in terms of geometrical‐ray spreading, group velocity and an energy integral. For the horizontal components there are both principal and additional components to describe the curvature of rays along the surface, as in the case of elastic body waves. The vertical component is decoupled from the horizontal components. With complex parameters the solutions for the dynamic ray tracing system correspond to Gaussian beams: the amplitude distribution is bell‐shaped along the direction perpendicular to the ray and the solution is regular everywhere, even at caustics. Most of the characteristics of Gaussian beams for 2‐D elastic body waves are also applicable to the surface wave case. At each frequency the solution may be regarded as a set of eigenfunctions propagating over a 2‐D surface according to the phase velocity mapping.