Quasi-isotropic approximation of ray theory for anisotropic media

Wave propagation in weakly anisotropic inhomogeneous media is studied by the quasi-isotropic approximation of ray theory. The approach is based on the ray-tracing and dynamic ray-tracing differential equations for an isotropic background medium. In addition, it requires the integration of a system of two complex coupled differential equations along the isotropic ray.
The interference of the qS waves is described by traveltime and polarization corrections of interacting isotropic S waves. For qP waves the approach leads to a correction of the traveltime of the P wave in the isotropic background medium.
Seismograms and particle-motion diagrams obtained from numerical computations are presented for models with different strengths of anisotropy.
The equivalence of the quasi-isotropic approximation and the quasi-shear-wave coupling theory is demonstrated. The quasi-isotropic approximation allows for a consideration of the limit from weak anisotropy to isotropy, especially in the case of qS waves, where the usual ray theory for anisotropic media fails.     

Computation of qP-wave rays, traveltimes and slowness vectors in orthorhombic media

The eikonal equation is the equation of the phase slowness surface for isotropic and anisotropic media. In general anisotropic media, there is no simple explicit expression for the phase slowness surface. An approximate expression of the eikonal equation may be obtained in weakly anisotropic media. In orthorhombic media, the approximate eikonal equation of the qP wave is the sum of an ellipsoidal form and a more complicated term. The ellipsoidal form corresponds to what we call ellipsoidal anisotropy. Ray equations written in the Hamiltonian formulation are characteristics of the eikonal equation. Ray perturbation theory may be used to compute changes in ray paths and physical attributes (traveltime, polarization, amplitude) due to changes in the medium with respect to a reference medium. Examples obtained in homogeneous orthorhombic media show that a reference medium with ellipsoidal anisotropy is a better choice to develop the perturbation approach than an isotropic reference medium. Models with strong anisotropy can be considered. The comparison with results obtained by an exact ray program shows a relative traveltime error of less than 0.5 per cent for a model with relatively strong anisotropy. We propose a finite element approach in which the medium is divided into a set of elements with polynomial elastic parameter distributions. Inside each element, using a perturbation approach, analytical expressions for rays and traveltimes are obtained Ray tracing reduces to connecting these analytical solutions at the vertices of the cells.     

First-order P-wave ray synthetic seismograms in inhomogeneous weakly anisotropic media

We propose approximate equations for P -wave ray theory Green's function for smooth inhomogeneous weakly anisotropic media. Equations are based on perturbation theory, in which deviations of anisotropy from isotropy are considered to be the first-order quantities. For evaluation of the approximate Green's function, earlier derived first-order ray tracing equations and in this paper derived first-order dynamic ray tracing equations are used.
The first-order ray theory P -wave Green's function for inhomogeneous, weakly anisotropic media of arbitrary symmetry depends, at most, on 15 weak-anisotropy parameters. For anisotropic media of higher-symmetry than monoclinic, all equations involved differ only slightly from the corresponding equations for isotropic media. For vanishing anisotropy, the equations reduce to equations for computation of standard ray theory Green's function for isotropic media. These properties make the proposed approximate Green's function an easy and natural substitute of traditional Green's function for isotropic media.
Numerical tests for configuration and models used in seismic prospecting indicate negligible dependence of accuracy of the approximate Green's function on inhomogeneity of the medium. Accuracy depends more strongly on strength of anisotropy in general and on angular variation of phase velocity due to anisotropy in particular. For example, for anisotropy of about 8 per cent, considered in the examples presented, the relative errors of the geometrical spreading are usually under 1 per cent; for anisotropy of about 20 per cent, however, they may locally reach as much as 20 per cent.     

Effects of slight anisotropy on surface waves

We present a complete ray theory for the calculation of surface-wave observables from anisotropic phase-velocity maps. Starting with the surface-wave dispersion relation in an anisotropic earth model, we derive practical dynamical ray-tracing equations. These equations allow calculation of the observables phase, arrival-angle and amplitude in a ray theoretical framework. Using perturbation theory, we also obtain approximate expressions for these observables. We assess the accuracy of the first-order approximations by using both theories to make predictions on a sample anisotropic phase-velocity map. A comparison of the two methods illustrates the size and type of errors which are introduced by perturbation theory. Perturbation theory phase and arrival-angle predictions agree well with the exact calculation, but amplitude predictions are poor. Many previous studies have modelled surface-wave propagation using only isotropic structure, not allowing for anisotropy. We present hypothetical examples to simulate isotropic modelling of surface waves which pass through anisotropic material. Synthetic data sets of phase and arrival angle are produced by ray tracing with exact ray theory on anisotropic phase-velocity maps. The isotropic models obtained by inverting synthetic anisotropic phase data sets produce deceptively high variance reductions because the effects of anisotropy are mapped into short-wavelength isotropic structure. Inversion of synthetic arrival-angle data sets for isotropic models results in poor variance reductions and poor recovery of the isotropic part of the anisotropic input map. Therefore, successful anisotropic phase-velocity inversions of real data require the inclusion of both phase and arrival-angle measurements.     

Gaussian beams in elastic 2-D laterally varying layered structures

Summary. In a paper by Červený & Pšenčik, high-frequency Gaussian beams in elastic 2-D, laterally inhomogeneous, smooth media were investigated as asymptotic high-frequency solutions of elastodynamic equations, concentrated close to rays of P - and S -waves. This paper generalizes the above results for 2-D, laterally inhomogeneous, layered structures. Gaussian beams concentrated close to any multiply-reflected, possibly converted, ray are investigated. Gaussian beams are regular everywhere, including caustic regions. The paraxial ray approximation, which allows the wavefield in the zero-order ray approximation to be evaluated not only directly on the ray, but also in its vicinity, is derived as a limiting case of the Gaussian beams.     

Dynamic ray tracing on Lagrangian manifolds

Summary. We show that Maslov's extension of the WKBJ method allows an extension of the dynamic ray tracing to wavefields involving caustics of arbitrary form. If the receiver lies off the caustics, then the synthetic seismogram can be obtained by integrating the DRT system along a single ray joining the receiver to the source which may touch caustics. If the receiver-lies in the vicinity of a caustic then DRT has to be carried out along a bunch of rays covering a neighbourhood of the receiver. Our approach encompasses pre-stressed and/or anisotropic media. Initial boundary conditions for a point source embedded in an anisotropic elastic medium are also presented.     

Computation of Large Anisotropic Seismic Heterogeneities (CLASH)

A general tomographic technique is designed in order (i) to operate in anisotropic media; (ii) to account for the uneven seismic sampling and (iii) to handle massive data sets in a reasonable computing time. One modus operandi to compute a 3-D body wave velocity model relies on surface wave phase velocity measurements. An intermediate step, shared by other approaches, consists in translating, for each period of a given mode branch, the phase velocities integrated along ray paths into local velocity perturbations. To this end, we develop a method, which accounts for the azimuthal anisotropy in its comprehensive form. The weakly non-linear forward problem allows to use a conjugate gradient optimization. The Earth's surface is regularly discretized and the partial derivatives are assigned to the individual grid points. Possible lack of lateral resolution, due to the inescapable uneven ray path coverage, is taken into account through the a priori covariances on parameters with laterally variable correlation lengths. This method allows to efficiently separate the 2ψ and the 4ψ anisotropic effects from the isotropic perturbations. Fundamental mode and overtone phase velocity maps, derived with real Rayleigh wave data sets, are presented and compared with previous maps. The isotropic models concur well with the results of Trampert & Woodhouse. Large 4ψ heterogeneities are located in the tectonically active regions and over the continental lithospheres such as North America, Antarctica or Australia. At various periods, a significant 4ψ signature is correlated with the Hawaii hotspot track. Finally, concurring with the conclusions of Trampert & Woodhouse, our phase velocity maps show that Rayleigh wave data sets do need both 2ψ and 4ψ anisotropic terms.     

Wavefield smoothing and the effect of rough velocity perturbations on arrival times and amplitudes

Geometric ray theory is an extremely efficient tool for modelling wave propagation through heterogeneous media. Its use is, however, only justified when the inhomogeneity satisfies certain smoothness criteria. These criteria are often not satisfied, for example in wave propagation through turbulent media. In this paper, the effect of velocity perturbations on the phase and amplitude of transient wavefields is investigated for the situation that the velocity perturbation is not necessarily smooth enough to justify the use of ray theory. It is shown that the phase and amplitude perturbations of transient arrivals can to first order be written as weighted averages of the velocity perturbation over the first Fresnel zone. The resulting averaging integrals are derived for a homogeneous reference medium as well as for inhomogeneous reference media where the equations of dynamic ray tracing need to be invoked. The use of the averaging integrals is illustrated with a numerical example. This example also shows that the derived averaging integrals form a useful starting point for further approximations. The fact that the delay time due to the velocity perturbation can be expressed as a weighted average over the first Fresnel zone explains the success of tomographic inversions schemes that are based on ray theory in situations where ray theory is strictly not justified; in that situation one merely collapses the true sensitivity function over the first Fresnel zone to a line integral along a geometric ray.     

Computation of wave fields in inhomogeneous media — Gaussian beam approach

Summary. An asymptotic procedure for the computation of wave fields in two-dimensional laterally inhomogeneous media is proposed. It is based on the simulation of the wave field by a system of Gaussian beams. Each beam is continued independently through an arbitrary inhomogeneous structure. The complete wave field at a receiver is then obtained as an integral superposition of all Gaussian beams arriving in some neighbourhood of the receiver. The corresponding integral formula is valid even in various singular regions where the ray method fails (the vicinity of caustic, critical point, etc.). Numerical examples are given.     

Computation of high-frequency seismic wavefields in 3-D laterally inhomogeneous anisotropic media

Summary. An algorithm for the computation of travel times, ray amplitudes and ray synthetic seismograms in 3-D laterally inhomogeneous media composed of isotropic and anisotropic layers is described. All 21 independent elastic parameters may vary within the anisotropic layers. Rays and travel times are evaluated by numerical solution of the ray tracing equations. Ray amplitudes are determined by evaluating reflection/ transmission coefficients and the geometrical spreading along individual rays. The geometrical spreading is computed approximately by numerical measurement of the cross-sectional area of the ray tube formed by three neighbouring rays. A similar approximate procedure is used for the determination of the coefficients of the paraxial ray approximation. The ray paraxial approximation makes computation of synthetic seismograms on the surface of the model very efficient. Examples of ray synthetic seismograms computed with a program package based on the described algorithm are presented.     

Synthetic body wave seismograms for three-dimensional laterally varying media

Summary. Two methods of computing body wave synthetic seismograms in three-dimensional laterally varying media are discussed. Both these methods are based on the summation of Gaussian beams. In the first, the initial beam parameters are chosen at the source, in the second at the beam endpoints. Both these variants eliminate the ray method singularities. The expansion of the wavefield into plane waves may be considered as the limiting case of the first approach and the Chapman–Maslov method as the limiting case of the second approach. Computer algorithms are briefly described and numerical examples presented. In the first numerical example, the comparisons of the two approaches, based on summing Gaussian beams, with the reflectivity method indicate that the computed synthetic seismograms are satisfactorily accurate even in the caustic region. The next example suggests that the two methods discussed can be simply and effectively applied to 3-D laterally inhomogeneous structures.     

Finite-frequency sensitivity of body waves to anisotropy based upon adjoint methods

We investigate the sensitivity of finite-frequency body-wave observables to mantle anisotropy based upon kernels calculated by combining adjoint methods and spectral-element modelling of seismic wave propagation. Anisotropy is described by 21 density-normalized elastic parameters naturally involved in asymptotic wave propagation in weakly anisotropic media. In a 1-D reference model, body-wave sensitivity to anisotropy is characterized by 'banana–doughnut' kernels which exhibit large, path-dependent variations and even sign changes. P -wave traveltimes appear much more sensitive to certain azimuthally anisotropic parameters than to the usual isotropic parameters, suggesting that isotropic P -wave tomography could be significantly biased by coherent anisotropic structures, such as slabs. Because of shear-wave splitting, the common cross-correlation traveltime anomaly is not an appropriate observable for S waves propagating in anisotropic media. We propose two new observables for shear waves. The first observable is a generalized cross-correlation traveltime anomaly, and the second a generalized 'splitting intensity'. Like P waves, S waves analysed based upon these observables are generally sensitive to a large number of the 21 anisotropic parameters and show significant path-dependent variations. The specific path-geometry of SKS waves results in favourable properties for imaging based upon the splitting intensity, because it is sensitive to a smaller number of anisotropic parameters, and the region which is sampled is mainly limited to the upper mantle beneath the receiver.     

Rayleigh-wave dispersion function for a transversely isotropic layered spherical earth using the Thomson-Haskell matrix method

We consider the two coupled differential equations of the two radial functions appearing in the displacement components of spheroidal oscillations for a transversely isotropic (TI) medium in spherical coordinates. Elements of the layer matrix have been explicitly written—perhaps for the first time—to extend the use of the Thomson-Haskell matrix method to the derivation of the dispersion function of Rayleigh waves in a transversely isotropic spherical layered earth. Furthermore, an earth-flattening transformation (EFT) is found and effectively used for spheroidal oscillations. The exponential function solutions obtained for each layer give the dispersion function for TI spherical media the same form as that on a flat earth. This has been achieved by assuming that the five elastic parameters involved vary as r ^p and that the density varies as r ^p-2, where p is an arbitrary constant and r is the radial distance. A numerical illustration with p = - 2 shows that, in spite of the inhomogeneity assumed within layers, the results for spherical harmonic degree n , versus time period T , obtained here for the Primary Reference Earth Model (PREM), agree well with those obtained earlier by other authors using numerical integration or variational methods. The results for isotropic media derived here are also in agreement with previous results. The effect of transverse isotropy on phase velocity for the first two modes of Rayleigh waves in the period range 20 to 240 s is calculated and discussed for continental and oceanic models.     

Wave propagation in transversely isotropic and periodically layered isotropic media

Summary. A set of stable algorithms for computing synthetic seismograms in attenuating transversely isotropic media is presented. The structures of these algorithms for anisotropic media are formally equivalent to their counterparts for isotropic media. The seismic responses of a periodically layered isotropic medium are compared with those of its long-wave equivalent transversely isotropic medium. The synthetics for the two media show observable differences in the range of frequencies considered. The differences are small in the P -waves, but partly large in later arrivals.     

On the inversion of Love- and Rayleigh-wave dispersion and implications for Earth structure and anisotropy

Summary. Various factors can make it difficult to explain observations of Love- and Rayleigh-wave dispersion with the same relatively simple isotropic model. These factors include systematic errors which might occur in determinations of observed group and phase velocities, lateral variations in structure along the path of travel, and the attempt to explain observations with a model comprised of only a small number of thick layers. The last of these factors is illustrated by an inversion of dispersion data in the central United States where shear-wave anisotropy had previously been invoked as one way to explain incompatible Love- and Rayleigh-wave velocities. It is shown that the data can be satisfied equally well by an isotropic model consisting of several thin layers.
In cases where the incompatibility of Love- and Rayleigh-wave data might be produced by intrinsic anisotropy, it is necessary to invert those data using an anisotropic theory rather than by separate isotropic inversions of Love and Rayleigh waves. Inversions of fundamental-mode data for a region of the Pacific, assuming anisotropic media in which the layers are transversely isotropic with a vertical axis of symmetry, lead to models which are highly non-unique. Even if the inversions solve only for shear velocities in the litho-sphere and asthenosphere it is not possible, without supplementary information, to ascertain the depth interval over which anisotropy occurs or to determine the thickness of the lithosphere or asthenosphere with much precision.     

Transverse isotropy of thinly layered media

Summary. Three problems of seismic anisotropy in thinly layered media (TPM) are discussed: (1) A dependence is established for the character of the ray velocity of longitudinal low-frequency waves on the ratio of P - and S -wave velocities in thin layers. (2) Conditions are specified for cusps on SV -wave surfaces. Nomograms are suggested for quick estimation of these conditions. (3) A comparison is made between TPM anisotropy and other types of transversely isotropic media.     

Caustics and focusing produced by sedimentary basins: applications of catastrophe theory to earthquake seismology

Summary. As high-frequency elastic waves propagate through real media, it is common for caustics and focusing to occur. Typically, rays may envelop a caustic surface in space, or exceptionally they may all coalesce at a focal point. In strong motion seismology, the observed large fluctuations in peak acceleration and intensity of ground shaking may just be a consequence of focusing and caustics created by waves propagating through irregularly shaped sedimentary basins. These basins, acting as deformed optical lenses, are capable of producing a complex network of patches and seemingly isolated pockets of intensified damage or high intensity shaking where the caustic intersects the ground surface. We adapt methods from optics and catastrophe theory to study the properties of caustics induced by typical sedimentary basins. Several hypothetical examples are shown that reflect the fact that these properties are useful to assess quantitatively the degree of wavefield amplification to be expected. A good correlation is found between actual damage patterns and caustic locations computed for the Caracas, Venezuela earthquake of 1967.     

Investigation of 3D wavefields of reflected shear-waves and converted waves: mathematical modelling and reflection data processing

Summary. Seismic investigations using shear-wave and converted wave techniques show that very often reflected PS - and SS -waves have anomalous polarizations ( accessory components ). This phenomenon cannot be explained in terms of isotropic models with dipping boundaries. Computations of synthetic seismograms of reflected PS - and SS -waves were made for different models of transversely isotropic media with dipping anisotropic symmetry axes not normal to the boundaries. Synthetic seismograms were computed by ray techniques using an optimization algorithm to construct all rays arriving at a given receiver. These computations indicate that accessory components arise when the medium above the boundary is anisotropic, where they are caused by the constructive interference of qSV - and qSH -waves. If a low-velocity layer is present, displacement vectors of both waves have horizontal projections which are approximately orthogonal. The algorithm for wave separation is presented and some results of its use are given.     

Anisotropic tomography of P-wave traveltimes

The ray path of a P -wave is specified in terms of the ray parameter and three Euler angles. the P -wave traveltime depends only on the ray parameter for a spherically symmetric earth. If we introduce an aspherical perturbation, including general ani-sotropy, the dependence on Euler angles can be expanded in terms of the rotation matrix for a fixed ray parameter. If the perturbation is isotropic, the expansion coefficients satisfy certain relations which may be used to obtain definite evidence for anisotropy rather than isotropic lateral heterogeneity.