Guided wave propagation in laterally varying media - I. Theoretical development

Summary. The propagation of surface waves in a laterally varying medium can be described by representing the wavetrain as a superposition of modal contributions for a reference structure. As the guided waves propagate through a heterogeneous zone the modal coefficients needed to describe the wavetrain vary with position, leading to interconversions between modes and reflection into backward travelling modes. The evolution of the modal terms may be described by a set of first-order differential equations which allow for coupling to both forward and backward travelling waves; the coefficients in these equations depend on the differences between the actual structure and the reference structure. This system is established using the orthogonality properties of the modal eigenfunctions and is valid for SH -waves, P - SV -waves and full anisotropy.
The reflected and transmitted wavefields for a region of heterogeneity can be related to the incident wave by introducing reflection and transmission matrices which connect the modal coefficients in these fields to those in the incident wavetrain. By considering a sequence of models with increasing width of heterogeneity we are able to derive a set of Ricatti equations for the reflection and transmission matrices which may be solved by initial value techniques. This avoids an awkward two-point boundary value problem for a large number of coupled equations. The method is demonstrated for 1 Hz Lg - and Sn -waves in a multilayered model for which there are 19 coupled modes.
The method is applicable to three-dimensional heterogeneity, and we are able to show that the interconversion between Love and Rayleigh waves, in the presence of gradients in seismic properties transverse to the propagation path, leads to a net rate of increase of the transverse components of the seismogram at the expense of the other components.     

Regional tomographic inversion of the amplitude and phase of Rayleigh waves with 2-D sensitivity kernels

In this study, we test the adequacy of 2-D sensitivity kernels for fundamental-mode Rayleigh waves based on the single-scattering (Born) approximation to account for the effects of heterogeneous structure on the wavefield in a regional surface wave study. The calculated phase and amplitude data using the 2-D sensitivity kernels are compared to phase and amplitude data obtained from seismic waveforms synthesized by the pseudo-spectral method for plane Rayleigh waves propagating through heterogeneous structure. We find that the kernels can accurately predict the perturbation of the wavefield even when the size of anomaly is larger than one wavelength. The only exception is a systematic bias in the amplitude within the anomaly itself due to a site response.
An inversion method of surface wave tomography based on the sensitivity kernels is developed and applied to synthesized data obtained from a numerical simulation modelling Rayleigh wave propagation over checkerboard structure. By comparing recovered images to input structure, we illustrate that the method can almost completely recover anomalies within an array of stations when the size of the anomalies is larger than or close to one wavelength of the surface waves. Surface wave amplitude contains important information about Earth structure and should be inverted together with phase data in surface wave tomography.     

Windowed f-k spectra of a three-dimensional wavefield excited by a point source in a two-dimensional multilayered elastic medium

When full 3-D modelling is too costly or cumbersome, computations of 3-D elastic wave propagation in laterally heterogeneous, multilayered 2-D geological structures may enhance considerably our ability to predict strong ground motion for seismological and engineering purposes. Towards this goal, we extend the method based on the combination of the thin-layer finite-element and boundary-element methods (TLFE-BEM) and calculate windowend f - k spectra of the 3-D wavefield. The windowed f - k spectra are spatially localized spectra from which the local properties of the wavefield can be extracted. The TLFE-BEM is particularly suited for calculating the complete wavefield where surface waves are dominant in multilayered media. The computations are performed in the frequency domain, providing the f - k spectra directly. From the results for the 3-D wavefield excited by a point source in a 2-D multilayered, sloped structure, it can be said that the phase velocity of the fundamental-mode Rayleigh wave in a laterally heterogeneous multilayered medium, estimated from the windowed f - k spectra, varies with the location of the point source. For the model calculated in this article, the phase velocity varies between the value for the flat layered structure of the thick-layer side and that for the structure just under the centre of the window. The exact subsurface structure just under the centre of an array in a laterally heterogeneous medium cannot be obtained if we use the f - k spectral analysis assuming a flat layered structure.     

A coupled mode examination of irregular waveguides including the continuum spectrum

Summary. A coupled mode theory is used to examine surface wave propagation in a laterally inhomogeneous acoustic waveguide. The theory is developed from the equations of motion for the pressure and velocity fields. The presence of lateral inhomogeneities in the form of varying layer thickness causes coupling among the discrete modes of the waveguide and radiation to the continuum. Expressions for the coupling coefficients among all mode types including coupling to the continuum spectrum are derived. The coupling coefficients are proportional to the horizontal derivative of the function describing the interface between layers of constant material properties but varying thickness. The coupled mode equations are solved in approximation for the case of a sinusoidal boundary and a sloping boundary. The results for radiation losses due to interaction with the irregular boundary of the waveguide are presented in analytical form, which clearly show the primary physical effects on the wavefield of the interaction. The far field amplitude of the scattered modes, excited by the interaction of some incident signal with a weak boundary irregularity, is modulated by the spatial Fourier transform of the irregularity.     

Scattering of surface waves modelled by the integral equation method

The integral equation method is used to model the propagation of surface waves in 3-D structures. The wavefield is represented by the Fredholm integral equation, and the scattered surface waves are calculated by solving the integral equation numerically. The integration of the Green's function elements is given analytically by treating the singularity of the Hankel function at   R = 0  , based on the proper expression of the Green's function and the addition theorem of the Hankel function. No far-field and Born approximation is made. We investigate the scattering of surface waves propagating in layered reference models imbedding a heterogeneity with different density, as well as Lamé constant contrasts, both in frequency and time domains, for incident plane waves and point sources.     

Reciprocity theorems for one-way wavefields

Acoustic reciprocity theorems have proved their usefulness in the study of forward and inverse scattering problems. The reciprocity theorems in the literature apply to the two-way (i.e. total) wavefield, and are thus not compatible with one-way wave theory, which is often applied in seismic exploration. By transforming the two-way wave equation into a coupled system of one-way wave equations for downgoing and upgoing waves it appears to be possible to derive 'one-way reciprocity theorems" along the same lines as the usual derivation of the 'two-way reciprocity theorems'. However, for the one-way reciprocity theorems it is not directly obvious that the 'contrast term' vanishes when the medium parameters in the two different states are identical. By introducing a modal expansion of the Helraholtz operator, its square root can be derived, which appears to have a symmetric kernel. This symmetry property appears to be sufficient to let the contrast term vanish in the above-mentioned situation.
The one-way reciprocity theorem of the convolution type is exact, whereas the one-way reciprocity theorem of the correlation type ignores evanescent wave modes. The extension to the elastodynamic situation is not trivial, but it can be shown relatively easily that similar reciprocity theorems apply if the (non-unique) decomposition of the elastodynamic two-way operator is done in such a way that the elastodynamic one-way operators satisfy similar symmetry properties to the acoustic one-way operators.     

Evidence for anisotropy in the upper mantle beneath Eurasia from the polarization of higher mode seismic surface waves

Summary. Analysis of NORSAR records and a number of Soviet microfilms reveals second-mode surface Caves propagating along paths covering a large part of Eurasia. These second modes in the 6–15-s period band are frequently disturbed by other surface-wave modes and by body-wave arrivals. However, in all cases, where the modes appear to be undisturbed and show normal dispersion, the Second Rayleigh modes have a slowly varying phase difference with the Second Love modes. This coupling has the particle motion of Inclined Rayleigh waves characteristic of surface-wave propagation in anisotropic media, where the anisotropy possesses a horizontal plane of symmetry. Numerical examination of surface wave propagating in Earth models, with an anisotropic layer in the upper mantle, demonstrate that comparatively small thicknesses of material with weak velocity anisotropy can produce large deviations in the polarizations of Inclined Rayleigh Second modes. In many structures, these inclinations are very sensitive to small changes in anisotropic orientation and to small changes in the surrounding isotropic structure. It is suggested that examination of second mode inclination anomalies of second mode surface waves may be a powerful technique for examining the detailed anisotropic structure of the upper mantle.     

A review of the effects of anisotropic layering on the propagation of seismic waves

Summary. This paper reviews recent work, much of it unpublished, on the effects of anisotropy on seismic waves, and lays the theoretical background for some of the other papers in this number of the Geophysical Journal .
The propagation of both body and surface waves in anisotropic media is fundamentally different from their propagation in isotropic media, although the differences in behaviour may be comparatively subtle and difficult to observe. One of the most diagnostic of these anomalies, which has been observed on some surface-wave trains, and should be evident in body-wave arrivals, is generalized, three-dimensional polarization, where the Rayleigh motion is coupled to the Love, and the P and SV motion is coupled to the SH . This coupling introduces polarization anomalies which may be used to investigate anisotropy within the Earth.     

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.     

Ray asymptotic propagation of the spectra of elastic surface waves in a waveguide with laterally smooth inhomogeneities

Summary. A technique based on ray asymptotics has been developed to propagate complex spectra of elastic normal mode surface waves in a waveguide with material and geometrical properties varying smoothly in the lateral directions. In the technique, the original problem defined in the unstretched coordinates has been transformed into an eiconal equation as well as into a certain number of transport equations defined in stretched coordinates.
The solution of the eiconal equation is equal to the solution of the eigenproblem of the eiconal operator A₀. Due to the self-adjointness of A₀, in each of the relevant local inner product spaces, LIPS, the solution of the eigenproblem, A₀ψ= v ψ results in the set { v _t} of real local eigenvalues and in the orthonormal system {ψ_t} of local eigenvectors.
As the Hamiltonian function of an initial value problem, each eigenvalues gives birth to a bicharacteristic curve as well as to the related ray. The introduction of the rays induces connections between the vertical cross-sections of the waveguide.
Finally, for each asymptotic order j , the LIPS-valued transport equations are reduced to a set of matricial propagation equations in the local spectral amplitude vectors, LSAVs. Consequently, a knowledge of the initial conditions at a vertical cross-section makes it possible to propagate the LSAVs along the rays of the relevant modes. However, to complete the propagation one needs, in addition to the initial values, information about certain additional quantities, non-diagonal terms of order j , diagonal terms of orders lower than j and the auxiliary boundary terms of orders from 1 to j . The treatment has been completed by the propagation of the modal phases along the relevant rays.     

Exact synthetic seismograms for an inhomogeneity in a layered elastic half-space

Summary. The propagation of a pulsed elastic wave in the following geometry is considered. An elastic half-space has a surface layer of a different material and the layer furthermore contains a bounded 3-D inhomogeneity. The exciting source is an explosion, modelled as an isotropic pressure point source with Gaussian behaviour in time.
The time-harmonic problem is solved using the null field approach (the T matrix method), and a frequency integral then gives the time-domain response. The main tools of the null field approach are integral representations containing the free space Green's dyadic, expansions in plane and spherical vector wave functions, and transformations between plane and spherical vector wave functions. It should be noted that the null field approach gives the solution to the full elastodynamic equations with, in principle, an arbitrarily high accuracy. Thus no ray approximations or the like are used. The main numerical limitation is that only low and intermediate frequencies, in the sense that the diameter of the inhomogeneity can only be a few wavelengths, can be considered.
The numerical examples show synthetic seismograms consisting of data from 15 observation points at increasing distances from the source. The normal component of the velocity field is computed and the anomalous field due to the inhomogeneity is sometimes shown separately. The shape of the inhomogeneity, the location and depth of the source, and the material parameters are all varied to illustrate the relative importance of the various parameters. Several specific wave types can be identified in the seismograms: Rayleigh waves, direct and reflected P -waves, and head waves.     

Time-lapse traveltime change of singly scattered acoustic waves

We present a technique based on the single-scattering approximation that relates time-lapse localized changes in the propagation velocity to changes in the traveltime of singly scattered waves. We describe wave propagation in a random medium with homogeneous statistical properties as a single-scattering process where the fluctuations of the velocity with respect to the background velocity are assumed to be weak. This corresponds to one of two end-member regimes of wave propagation in a random medium, the first being single scattering, and the second multiple scattering. We present a formulation that relates the change in the traveltime of the scattered waves to a localized change in the propagation velocity by means of the Born approximation for the scattered wavefield. We validate the methodology with synthetic seismograms calculated with finite differences for 2-D acoustic waves. Potential applications of this technique include non-destructive evaluation of heterogeneous materials and time-lapse monitoring of heterogeneous reservoirs.     

On the nature of regional seismic phases–II. On the influence of structural barriers

The blockage of the L _g wave by crustal barriers such as continental margins and graben structures has long been recognized as providing a very useful tool for mapping large-scale lateral crustal variations along the propagation path. Numerical simulation of L _g -wave propagation in complex anelastic media using the pseudospectral method provides insight into the nature of the propagation process using both snapshots of the wavefield and synthetic seismograms. A variety of 2-D structures have been investigated, including the influence of sediments, crustal thickness and attenuation.
Thick sedimentary basins covering a graben structure can have a major influence, since they remove L _g energy by generating P conversion and scattering–the principal mechanisms for strong L _g attenuation across a graben. The reduction of the L _g energy is reinforced by anelastic attenuation in the sediments as well as the influence of the gradually thinning crustal waveguide associated with an elevated Moho.
The extinction of L _g in a sequence of explosions fired across the central graben of the North Sea can be simulated by numerical calculations for the structure derived from refraction experiments.     

Polarization anomaly of Love waves caused by lateral heterogeneity

We calculate surface waves propagating in a laterally heterogeneous structure beneath the Kuril trench, where significant Love-wave polarization anomalies, called quasi-Love waves, are generated. Since 3-D wave propagation in the two-dimensionally heterogeneous structure can be assumed, we apply the 2.5-D finite difference method to the surface-wave calculations. The calculations show that a velocity contrast of 7 per cent at depths of less than 210 km beneath the Kuril trench cannot generate quasi-Love waves, and that an unlikely contrast of 20 per cent is required to generate clear quasi-Love waves. The possible cause of the quasi-Love waves inferred from previous studies on coupled free oscillations is a lateral variation in azimuthal anisotropy. The lateral variation in azimuthal anisotropy beneath the Kuril trench suggests a change in the mantle flow induced by the subducting slab.     

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.     

Effect of 3-D viscoelastic structure on post-seismic relaxation from the 2004 M= 9.2 Sumatra earthquake

The 2004 M = 9.2 Sumatra–Andaman earthquake profoundly altered the state of stress in a large volume surrounding the ∼1400 km long rupture. Induced mantle flow fields and coupled surface deformation are sensitive to the 3-D rheology structure. To predict the post-seismic motions from this earthquake, relaxation of a 3-D spherical viscoelastic earth model is simulated using the theory of coupled normal modes. The quasi-static deformation basis set and solution on the 3-D model is constructed using: a spherically stratified viscoelastic earth model with a linear stress–strain relation; an aspherical perturbation in viscoelastic structure; a 'static' mode basis set consisting of Earth's spheroidal and toroidal free oscillations; a "viscoelastic" mode basis set; and interaction kernels that describe the coupling among viscoelastic and static modes. Application to the 2004 Sumatra–Andaman earthquake illustrates the profound modification of the post-seismic flow field at depth by a slab structure and similarly large effects on the near-field post-seismic deformation field at Earth's surface. Comparison with post-seismic GPS observations illustrates the extent to which viscoelastic relaxation contributes to the regional post-seismic deformation.     

Tectonic modelling of the Middle Jurassic synrift stratigraphy in the Oseberg–Brage area, northern Viking Graben

A finite difference model, allowing for episodic movements along different faults, is used to examine the effect of tectonics on the stratigraphic signature in the Oseberg–Brage area in the northern Viking Graben. Constraints are provided by local exploration and production well data and 3-D seismic coverage, and a regional depth-converted seismic line.
In the modelling, we focus on the influence of varying rates of fault movement on stratigraphic signatures such as upflank unconformities and changes in layer thickness. We couple the basinwide features of the northern Viking Graben with the fault-block-scale features of the Oseberg–Brage area by using parameter constraints derived by large-scale modelling as input for the local-scale model. In addition, subsidence patterns resulting from the basinwide model were used as background subsidence for the fault block model of the Oseberg–Brage area.
The model results indicate that the alternating activation of different faults with varying extension rates can cause stratigraphic features such as unconformities, condensation and onlap/offlap patterns. Onlap occurs during periods of low extension rates. An increase in extension rate along a fault causes footwall uplift, resulting in condensation or upflank erosion yielding unconformities. This influence can also affect sub-basins further away from the fault. Downdip layer thickening reflects the local tilting of fault blocks.
The coupling of the local and regional scales turns out to be essential in explaining the stratigraphy of the Oseberg–Brage area: basinward and, notably, updip layer thickening as observed on some of the fault blocks can only be explained by activity of the boundary fault on the opposing, western margin of the northern Viking Graben.
     

The Acoustics of Pipe Arrays

Summary A theoretical analysis is given for the acoustical behaviour of the pipe-microbarograph systems used to detect acoustic gravity waves and other modes of infrasound. It is shown how to compute the response of the microbarograph to a fluctuating pressure at any one inlet port of the pipe and how the results of such computations may be used to calculate the response to a plane sound wave traversing the system.
The analysis is illustrated by numerical examples obtained by means of a computer program. These examples confirm that the tapered tube modelled after Daniels' line microphone has very good characteristics, but that good results may also be obtained using pipes of uniform bore. The work leans heavily on Benade's calculations of sound propagation in a circular conduit.     

Waveform inversion of mantle Love waves:the Born seismogram approach

Summary. Normal mode theory, extended to the slightly laterally heterogeneous earth by the first-order Born approximation, is applied to the waveform inversion of mantle Love wave (200–500 s) for the Earth's lateral heterogeneity at l = 2 and a spherically symmétric anelasticity ( Q _μ) structure. The data are from the Global Digital Seismograph Network (GDSN). The l =2 pattern is very similar to the results of other studies that used either different méthods, such as phase velocity measurements and multiplet location measurements, or a different data set, such as mantle Rayleigh waves from different instruments. The results are carefully analysed for variance reduction and are most naturally explained by heterogeneity in the upper 420 km. Because of the poor resolution of the data set for the deep interior, however, a fairly large heterogeneity in the transition zones, of the order of up to 3.5 per cent in shear wave velocity, is allowed. It is noteworthy that Love waves of this period range cannot constrain the structure below 420 km and thus any model presented by similar studies below this depth are likely to be constrained by Rayleigh waves (spheroidal modes) only.
The calculated modal Q values for the obtained Q _μ model fall within the error bars of the observations. The result demonstrates the discrepancy of Rayleigh wave Q and Love wave Q and indicates that care must be taken when both Rayleigh and Love wave data, including amplitude information, are inverted simultaneously.
Anomalous amplitude inversions of G2 and G3, for example, are observed for some source-receiver pairs. This is due to multipathing effects. One example near the epicentral region, which is modelled by the obtained l = 2 heterogeneity, is shown.     

Diffraction of P, S and Rayleigh waves by three-dimensional topographies

The diffraction of P, S and Rayleigh waves by 3-D topographies in an elastic half-space is studied using a simplified indirect boundary element method (IBEM). This technique is based on the integral representation of the diffracted elastic fields in terms of single-layer boundary sources. It can be seen as a numerical realization of Huygens principle because diffracted waves are constructed at the boundaries from where they are radiated by means of boundary sources. A Fredholm integral equation of the second kind for such sources is obtained from the stress-free boundary conditions. A simplified discretization scheme for the numerical and analytical integration of the exact Green's functions, which employs circles of various sizes to cover most of the boundary surface, is used.
The incidence of elastic waves on 3-D topographical profiles is studied. We analyse the displacement amplitudes in the frequency, space and time domains. The results show that the vertical walls of a cylindrical cavity are strong diffractors producing emission of energy in all directions. In the case of a mountain and incident P, SV and SH waves the results show a great variability of the surface ground motion. These spatial variations are due to the interference between locally generated diffracted waves. A polarization analysis of the surface displacement at different locations shows that the diffracted waves are mostly surface and creeping waves.