Seismic body waves in anisotropic media: reflection and refraction at a plane interface

Summary. A formulation is derived for calculating the energy division among waves generated by plane waves incident on a boundary between generally anisotropic media. A comprehensive account is presented for P, SV and SH waves incident from an isotropic half-space on an orthorhombic olivine half-space, where the interface is parallel to a plane of elastic symmetry. For comparison, a less anisotropic medium having transverse isotropy with a horizontal axis of symmetry is also considered. The particle motion polarizations of waves in anisotropic medium differ greatly from the polarizations in isotropic media, and are an important diagnostic of the presence of anisotropy. Incident P and SV waves generate quasi- SH waves, and incident SH waves generate quasi- P and quasi- SV waves, often of considerable relative magnitude. The direction of energy transport diverges from the propagation direction.     

Seismic body waves in anisotropic media: synthetic seismograms

Summary. Synthetic seismograms and particle motion diagrams are computed for simple, layered Earth models containing an anisotropic layer. The presence of anisotropy couples the P, SV and SH wave motion so that P waves incident on the anisotropic layer from below produce P, SV and small-amplitude SH waves at the surface both the P velocity and the amplitudes of the converted phases vary with azimuth. Significant SH amplitudes may be generated even when the wavelength of the P wave is much greater than the thickness of the anisotropic layer. Incident SV or SH waves may each generate large amplitudes of both SV and SH motion. This strong coupling is largely independent of the degree of velocity anisotropy of the medium. The arrivals from short-period S waves exhibit S-wave splitting, but arrivals from longer period S waves superpose into a modified waveform. This strong coupling does not allow the arrival of separate phases with pure SV and SH polarization except along directions of symmetry where the motion decouples.     

Elastodynamic and elastostatic Green tensors for homogeneous weak transversely isotropic media

An explicit analytical formula for the complete elastodynamic Green tensor for homogeneous unbounded weak transversely isotropic media is presented. The formula was derived by analytical calculations of higher-order approximations of the ray series. The ray series is finite and consists of seven non-zero terms. The formula for the Green tensor is complete and correct for the whole frequency range, thus it describes correctly the wavefield at all distances and at all directions including the shear-wave singularity direction. The Green tensor consists of P, SV and SH far-field waves and four coupling waves. Three of them couple P and SV waves, and the fourth wave couples the SV and SH waves. The P-SV coupling waves behave similarly to the near-field waves in isotropy. However, the SV-SH coupling wave, which is called 'shear-wave coupling', behaves exceptionally and it has no analogy in the Green tensor for isotropy. The formula for the elastostatic Green tensor is also derived.     

Energy flux in viscoelastic anisotropic media

We study properties of the energy-flux vector and other related energy quantities of homogeneous and inhomogeneous time-harmonic P and S plane waves, propagating in unbounded viscoelastic anisotropic media, both analytically and numerically. We propose an algorithm for the computation of the energy-flux vector, which can be used for media of unrestricted anisotropy and viscoelasticity, and for arbitrary homogeneous or inhomogeneous plane waves. Basic part of the algorithm is determination of the slowness vector of a homogeneous or inhomogeneous wave, which satisfies certain constraints following from the equation of motion. Approaches for determination of a slowness vector commonly used in viscoelastic isotropic media are usually difficult to use in viscoelastic anisotropic media. Sometimes they may even lead to non-physical solutions. To avoid these problems, we use the so-called mixed specification of the slowness vector, which requires, in a general case, solution of a complex-valued algebraic equation of the sixth degree. For simpler cases, as for SH waves propagating in symmetry planes, the algorithm yields simple analytic solutions. Once the slowness vector is known, determination of energy flux and of other energy quantities is easy. We present numerical examples illustrating the behaviour of the energy-flux vector and other energy quantities, for homogeneous and inhomogeneous plane P , SV and SH waves.     

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.     

Anisotropy of the mantle inferred from observations of P to S converted waves

Summary. Seismic anisotropy has been previously studied at depths usually not exceeding 100 or 150 km. In this paper we present a method of analysis of seismic records which is very sensitive to azimuthal anisotropy and is applicable at almost any depth range. The idea of the method is to detect and analyse the SH -component of the waves, converted from P to S in the mantle. The procedure of record processing includes frequency filtering, axis rotation, transformation of the record to a standard form, stacking the standardized SH -component records of many seismic events, and the harmonic analysis of amplitude as a function of the direction of wave propagation. When applied to the long-period records of NORSAR the procedure detected a converted wave with the properties implying the possibility of its propagation in a transversely isotropic medium with a horizontal axis of symmetry . Our preferred model postulates anisotropy of ∼ 1 per cent in a layer 50 km thick at the base of the upper mantle.     

Shear-wave polarization anisotropy in the Pacific Basin

Summary. Inversion of 295. Love- and Rayleigh-wave phase travel times across the Pacific Basin has yielded a structure which has a channel that is anisotropic with respect to the polarization of shear waves. The velocity of SH waves is approximately 4.24 km/s, and the velocity of SV waves is approximately 4.10 km/s in the low-velocity channel. The lid to the channel is isotropic with respect to the polarization of S waves and the velocity is approximately 4.60 km/s. The lid to the low-velocity channel increases in thickness with lithospheric age at the expense of the channel, and its thickness is apparently still increasing at a sea-floor age of 150 Myr.
These results can be explained in terms of a model with both randomly-and preferentially-oriented, liquid-filled cracks in the channel. In the model, it is assumed that the liquid-filled cracks are due to partial melting in the channel, and that any preferred orientation is caused by a shear-flow gradient resulting from differential motion between the lid and the deeper parts of the mantle.     

Numerical simulation of the propagation of P waves in fractured media

We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q ⁻¹ of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q ⁻¹ of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) ⁻¹ in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.     

Full reconstruction of a layered elastic medium from P-SV slant-stack data

Summary. We develop a méthod of reconstructing the elastic paraméters as functions of depth, for a horizontally stratified, isotropic elastic half-space. Unlike previous schemes, which have been able to retrieve the shear wave speed and density from SH seismograms slant stacked at two angles, our méthod makes use of P - SV data at a single stacking paraméter to obtain all three elastic constants. The data required are the elements of the full reflection matrix at the surface, corresponding to measurements of two separate components of the response to two independent sources, one explosive, the other generating shear waves.
In developing this inverse scheme fundamental differences emerge between the acoustic or SH problem, and the coupled P - SV case, the most important being in the nature of the interfacial scattering matrix. We show that it is not possible to make use of the downward reflection data for an interface to determine directly the remaining reflection and transmission coefficients, but that the scattering data may be completed by applying a simple iterative procedure at each interface.
We show the result of applying our inverse scheme to seismograms generated for a six-layered model, including a low-velocity layer. We are able to reconstruct both wave speeds and the density as functions of depth, all quantities being in close agreement with the original model.     

The Backus–Gilbert approach to the 3-D structure in the upper mantle – II. SH and SV velocity

Summary. Phase velocity variations obtained in the previous paper are inverted by the Backus–Gilbert method for the velocity structure of the upper mantle. Spheroidal modes and toroidal modes in the period range of 125–260 s are used in the inversion. The data cannot constrain all six parameters in a transversely isotropic medium and we chose to perturb only two parameters, SH and SV velocities. SV velocities are resolved between the depths of about 200 and 400 km and SH velocities between 0 and 200 km. Resolution kernels have half-peak widths of about 200–300 km in depth, becoming broader for deeper target depths. SV velocity kernels show secondary peaks near the surface of the Earth, with widths varying from 50 to 100 km. The deeper the target depths, the wider the secondary peaks near the surface. SH velocity kernels do not possess such secondary peaks. The trade-off between SV and SH velocities is small. SV velocity is essentially determined by spheroidal modes and SH velocity by toroidal modes. Because of the broad width of the resolution kernels, the structure in the resolved region is difficult to detect from our data set; for example the differences in SV velocity structure between 250 and 350 km or the differences in SH velocity between 100 and 200 km are difficult to distinguish. Considering the horizontal resolution of about 2000 km, obtained in the previous paper, averaging kernels for 3-D structure are quite elongated in the horizontal dimension.     

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.     

Reflection—refraction of general P-and type-I S-waves in elastic and anelastic solids

Summary. The reflection and refraction of general (homogeneous or inhomo-geneous) plane P and type-I S ( SV ) body waves incident on plane boundaries are considered for general linear viscoelastic solids. Reflection—refraction laws, physical characteristics of the waves, and the nature of critical angles are examined in detail at welded boundaries and a free surface. General visco-elasticity with no low-loss approximations predicts that contrasts in intrinsic absorption at boundaries give rise to inhomogeneous reflected and refracted waves with elliptical particle motions, velocities and maximum attenuations that vary with frequency and angle of incidence, energy propagation at speeds and directions different from phase propagation, phase propagation that in general is parallel to the boundary for at most one angle of incidence, and reflection—transmission coefficients dependent on energy flow due to wave interaction. None of these physical characteristics are predicted for waves incident on boundaries that respond instantaneously.     

Flexure of the ocean lithosphere from island uplift, bathymetry and geoid height observations: the Society Islands

Summary. The response of a stratified elastic medium can be conveniently characterized by the Green's tensor for the medium. For coupled seismic wave propagation ( P—SV or fully anisotropic), the Green's tensor may be constructed directly from two matrices of linearly independent displacement solutions. Rather simple forms for the Green's tensor can be found if each displacement matrix satisfies one of the boundary conditions on the seismic field. This approach relates directly to 'reflection matrix' representations of the seismic field.
For a stratified elastic half space the Green's tensor is used to give a spectral representation for coupled seismic waves. By means of a contour integration a general completeness relation is obtained for the 'body wave' and 'surface wave' parts of the seismic field. This relation is appropriate for SH and P–SV waves in an isotropic medium and also for full anisotropy.     

Seismic anisotropy in the core–mantle transition zone

Split S waves observed at Hockley, Texas from events in the Tonga–Fiji region of the southwest Pacific show predominantly vertically polarized shear-wave ( SV  ) energy arriving earlier than horizontally polarized ( SH ) energy for rays propagating horizontally through D" . After corrections are made for the effects of upper-mantle anisotropy beneath Hockley, a time lag of 1.5 to 2.0  s remains for the furthest events (93.9°–100.6° ), while the time lags of the nearer observations (90.5°–92.9° ) nearly disappear. At closer distances, the S waves from these same events do not penetrate as deeply into the lower mantle, and are not split. These observations suggest that a patch of D" beneath the central Pacific is anisotropic, while the mantle immediately above the patch is isotropic. The thickness of the anisotropic zone appears to be of the order of 100–200  km.
  Observations of shear-wave splitting have previously been made for paths that traverse D" under the Caribbean and under Alaska. SH leads SV , the reverse of the Hockley observations, but in these areas the fact that SV  leads SH in the HKT data shown here suggests a different sort of anisotropy under the central Pacific from that under Alaska and the Caribbean. The case of SH travelling faster than SV  is consistent with transverse isotropy with a vertical axis of symmetry (VTI) and does not require variations with azimuth. The case of SV  leading SH is consistent with transverse isotropy with a horizontal axis of symmetry (HTI), an azimuthally anisotropic medium, and with a VTI medium formed by a hexagonal crystal. Given that (Mg,Fe)SiO₃ perovskite appears unlikely to form anisotropic fabrics on a large scale, the presence of anisotropy may point to chemical heterogeneity in the lowermost mantle, possibly due to mantle–core interactions.     

Measuring anisotropy in rocks using laser-generated ultrasound

Summary. Measurements using standard contacting piezoelectric transducers and non-contacting laser sources and detectors, have been investigated for the study of ultrasonic anisotropy in rocks. An ultrasonic polariscope has been constructed in order to obtain reproducible travel-time and amplitude measurements. Three case studies are described to demonstrate the apparatus, namely isotropic halite, anisotropic calcite and transversely anisotropic mudstone. A novel technique has been developed in order to construct pseudo-particle motion diagrams, to highlight shear-wave birefringence in rock samples using 2.25 MHz transducers. A pulsed laser has been used to generate compressional and shear waves for comparison with piezoelectric transducer results. The pulses generated by laser irradiation have many advantages for the study of velocity and attenuation anisotropy because of their known characteristics, broad bandwidth and high level of reproducibility. The use of a non-contacting laser source and detector eliminates the need for elaborate coupling agents, stress bonding or immersion techniques. Point-source and line-focusing of the laser beam provides an indirect method of studying shear-wave polarization phenomena. Results from rotation of the line-focused laser beam and rotation of piezoelectric shear-wave transducers with respect to anisotropy, are compared for both velocity and amplitude phenomena in an anisotropic rock sample.     

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.     

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.     

ELASTIC WAVES IN A CONTINUOUSLY STRATIFIED MEDIUM

Summary. Approximate methods are applied to problems of internal reflection in a continuous medium when the changes of properties in a wave-length are small. It is found that total internal reflection gives a phase shift of 1/2π in harmonic waves and hence a blunting of a pulse and the introduction of some oscillation.
In the absence of a free surface P and SV movements are substantially independent even if total reflection takes place. SH and SV can travel in an internal layer of low velocity, but will not produce much motion at a free surface unless the total variation of properties is small. P waves in such a channel will suffer attenuation by conversion to SV at the free surface.     

The nature of particle motion in regional seismograms and its utilization for phase identification

The particle motion of regional arrivals is frequently treated in automatic phase-recognition schemes as that appropriate to simple P or S waves incident on an elastic, laterally homogeneous half-space. This model implies that the motion in ' P -type' phases can be described in terms of a single, generalized signal process and ' S -type' phases in terms of two independent processes ( SV and SH ) and thus, all regional arrivals could be fully characterized by three components of motion. In this paper, we present anlyses of the particle-motion patterns of various regional arrivals recorded at the ARCESS array from closely spaced events in the Kola Peninsula. We have found that only Pn -particle motion, described in terms of two independent signal processes, can be reliably characterized by three-component recordings. On the other hand, the various regional arrivals following Pn , such as Pg and Sn and Lg , can only be poorly characterized on the basis of three-component recordings alone. The reason is that these arrivals must be described in terms of more than two independent generalized signal processes, at least three for Pg and Sn , and possibly up to five for Lg . Recognition of these phases will thus require the use of more sensors than signal processes in the observing sensor configuration, such as three-component sensors combined with a small tripartite array. We have investigated the feasibility of adaptive, automatic recognition of regional arrivials by a wavefield extrapolation scheme utilizing such a mini-array. The process, which appear to be promising, adaptively learns the particle-motion patterns of individual arrivals, including complex site-response functions, from examples of closely located regional events.     

Seismic vertical array analysis for phase decomposition

We propose a vertical array analysis method that decomposes complex seismograms into body and surface wave time histories by using a velocity structure at the vertical array site. We assume that the vertical array records are the sum of vertically incident plane P and S waves, and laterally incident Love and Rayleigh waves. Each phase at the surface is related to that at a certain depth by the transfer function in the frequency domain; the transfer function is obtained by Haskell's matrix method, assuming a 1-D velocity structure. Decomposed P , S and surface waves at the surface are estimated from the vertical array records and the transfer functions by using a least-squares method in the frequency domain; their time histories are obtained by the inverse Fourier transform. We carried out numerical tests of this method based on synthetic vertical array records consisting of vertically incident plane P and S waves and laterally incident plane Love and Rayleigh waves. Perfect results of the decomposed P , S , Love and Rayleigh waves were obtained for synthetic records without noise. A test of the synthetic records in which a small amount of white noise was added yielded a reasonable result for the decomposed P , S and surface waves. We applied this method to real vertical array records from the Ashigara valley, a moderate-sized sedimentary valley. The array records from two earthquakes occurring at depths of 123 and 148 km near the array (epicentral distance of about 31 km) exhibited long-duration later phases. The analysis showed that duration of the decomposed S waves was a few seconds and that the decomposed surface waves appeared a few seconds after the direct S -wave arrival and had very long duration. This result indicated that the long-duration later phases were generated not by multireflected S waves, but by basin-induced surface waves.