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.     

Shear-wave polarizations near the North Anatolian Fault – I. Evidence for anisotropy-induced shear-wave splitting

Summary. The Turkish Dilatancy Projects (TDP1 in 1979 and TDP2 in 1980) recorded small earthquakes near the North Anatolian Fault with closely-spaced networks of three-component seismometers in order to investigate the possibility of diagnosing dilatancy from its effects of shear-wave propagation. This paper examines the polarizations of shear wavetrains recorded in the shear-wave window immediately above the earthquake foci. Abrupt changes in the orientation and/or ellipticity of the shear-wave polarizations are almost always observed during the first few cycles following the initial shear-wave arrival on each seismogram. The horizontal projections of the polarizations of the first shear-wave arrivals at any given station show nearly parallel alignments with approximately the same orientations at each of the recording sites (with one exception). It is difficult to explain this uniform alignment over a wide area in terms of scattering at the irregular surface topography or by earthquake focal mechanisms. We demonstrate that the shear-wave splitting is likely to be the result of anisotropy in the region above the earthquake foci, which could produce polarizations displaying the observed alignments. The temporal change of the azimuth of alignment, observed at one locality between 1979 and 1980, may be due to the release of a local stress anomaly by a very near earthquake.     

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.     

Effects of the free surface on shear wavetrains

Summary. The behaviour of shear-waves is of great importance in identifying and investigating seismic anisotropy in the Earth. However, shear wavetrains recorded at the Earth's surface do not always reflect the motion at depth, introducing practical problems of interpretation. Shear wavetrains incident on the surface of an isotropic half-space at angles less than critical (about 35°) are broadly preserved, but at greater angles substantial distortions can occur. For stations situated close to the source, as in local earthquake studies, the local SP phase, a radially polarized precursor to S , may occur. The behaviour at the surface of an anisotropic half-space is further complicated by the divergence of phase and energy propagation vectors. All of these complications suggest that detailed seismogram modelling is essential to any study of shear wave propagation in the Earth, and in particular to investigations of anisotropy-induced shear-wave splitting.     

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.     

Modelling the compliance of crustal rock—I. Response of shear-wave splitting to differential stress

We show that seismic shear waves may be used to monitor the in situ stress state of deep inaccessible rocks in the crust. The most widespread manifestation of the stress-related behaviour of seismic waves is the shear-wave splitting (shear-wave birefringence) observed in almost all rocks, where the polarizations of the leading split shear waves are usually subparallel to the direction of the local maximum horizontal stress. It has been recognized that such shear-wave splitting is typically the result of propagation through distributions of stress-aligned fluid-filled microcracks and pores, known as extensive-dilatancy anisotropy or EDA. This paper provides a quantitative basis for the EDA hypothesis. We model the evolution of anisotropic distributions of microcracks in triaxial differential stress, where the driving mechanism is fluid migration along pressure gradients between neighbouring microcracks and pores at different orientations to the stress field. This leads to a non-linear anisotropic poroelasticity (APE) model for the stress-sensitive behaviour of fluid-saturated microcracked rocks. A companion paper shows that APE modelling matches a range of observed phenomena and is a good approximation to the equation of state of a stressed fluid-saturated rock mass.     

Existence of transverse waves in anisotropic poroelastic media

Out of the four waves in an anisotropic poroelastic medium, two are termed as quasi-transverse waves. The prefix 'quasi' refers to their polarizations being nearly, but not exactly, perpendicular to direction of propagation. In this composite medium, unlike perfectly elastic medium, the propagation of a longitudinal wave along a phase direction may not be accompanied by transverse waves. The existence of a transverse wave in anisotropic poroelastic media is ensured by the two equations restricting the choice of elastic coefficients of porous aggregate as well as fluid–solid coupling. Necessary and sufficient conditions for the existence of transverse waves along the coordinate axes and in the coordinate planes for general anisotropy are discussed. The discussion is extended to the case of orthotropic materials and existence for few specific phase directions is also explored. The conditions for the transverse waves decided on the basis of their apparent polarizations, that is, particle motion being perpendicular to ray direction, are also discussed. For a particular numerical model, the existence of these apparent transverse waves is solved numerically for phase directions in coordinate planes. For general directions of phase propagation, the existence of these transverse waves is checked graphically for the chosen numerical model.     

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.     

Experimental study of the anisotropy of longitudinal and transverse waves from local earthquake records

Summary. The paper gives the results of a study of the anisotropy of seismic wave velocities within the Ashkhabad test field in Central Asia. The anisotropy was studied by analysing variations in the values of apparent velocities of first arrivals for epicentral distances ranging from 30 to 130 km and by analysing the delays (Δ t_s1-_s2 ) between the arrival times of shear waves with different polarizations.
The velocities of P -waves vary with azimuth from 5.3 to 6.27 km s^-1 and the velocities of S -waves vary from 3.15 to 3.5 km s^-1.
The delay times Δ t_S1 - _S2 depend on the direction of the propagation. The character of the variation of the propagation velocity of the longitudinal wave, the presence of two differently polarized shear waves S 1 and S 2 propagating at different velocities, and the character of the distribution of Δ t_S1 - _S2 on the stereogram suggest that the symmetry of the anisotropic medium is close to hexagonal with a nearly horizontal symmetry axis coinciding with the direction of maximal velocity. The azimuth of the symmetry axis of the medium is 140° and coincides with the direction of geological faults.     

Seismic-wave propagation through a cracked solid: polarization as a possible dilatancy diagnostic

summary . A new technique is presented for modelling the elastic constants of cracked structures with application to systems with weak concentrations of parallel cracks, and of simple biplanar and triplanar cracks. The velocities and V_p/V_s ratios of these anisotropic structures are used to provide quantitative models for some earthquake precursors. These results indicate the great importance of crack geometry to the behaviour of precursors. The velocities of saturated cracks appear to favour the dilatancy-diffusion model of precursory phenomena. Synthetic seismograms are calculated for propagation through possible dilatancy zones. The seismograms show some characteristic features which may be useful for the investigation of earthquake dilatancy.     

Geometrical structure of shear wave surfaces near singularity directions in anisotropic media

There are three types of surfaces which are used for studying wave propagation in anisotropic media: normal surfaces, slowness surfaces and wave surfaces. Normal surfaces and slowness surfaces have been researched in detail. Wave surfaces are the most complicated and comparatively poorly known compared with the other two. Areas of complicated geometrical structure of the wave surfaces are located in the vicinity of conical acoustic axes. There is an elliptical hole on the quick shear wave surface and complicated folds and cusps on the slow shear wave surface. Decomposition of the slow shear wave surface into smooth sheets is used for the study of its geometrical structure. Complexity of shear wave surfaces can be expressed by the number of waves corresponding to a fixed ray. An original approach to the calculation of wave normals depending on ray direction is presented.     

Guided waves in three-dimensional structures

A new formulation for the propagation of surface waves in three-dimensionally varying media is developed in terms of modal interactions. A variety of assumptions can be made about the nature of the modal field: a single set of reference modes, a set of local modes for the structure beneath a point, or a set of local modes for a laterally varying reference structure. Each modal contribution is represented locally as a spectrum of plane waves propagating in different directions in the horizontal plane. The influence of 3-D structure is included by allowing coupling between different modal branches and propagation directions. For anisotropic models, with allowance for attenuation, the treatment leads to a set of coupled 2-D partial differential equations for the weight functions for different modal orders.
The representation of the guided wavefield requires the inclusion of a full set of modes, so that, even for isotropic models, both Love and Rayleigh modes appear as different polarization states of the modal spectrum. The coupling equations describe the interaction between the different polarizations induced by the presence of the 3-D structure.
The level of lateral variation within the 3-D model is not required to be small. Horizontal refraction or reflection of the surface wavefield can be included by allowing for transfer between modes travelling in different directions. Approximate forms of the coupled equation system can be employed when the level of heterogeneity is small, for example the coupling between the fundamental mode and higher modes can often be neglected, or forward propagation can be emphasized by restricting the interaction to a limited band of plane waves covering the expected direction of propagation.     

Seismic anisotropy — the state of the art: II

Summary. The theory, causes, observations, and possible applications of seismic anisotropy in the Earth have developed considerably since the previous state of the art paper was published in 1977. The behaviour of waves in layered anisotropic media is now much better understood and the evidence for seismic anisotropy indicates that anisotropy is likely to be present throughout much of the crust and upper mantle. The top few hundred kilometres of the mantle appears to be anisotropic with the orientations aligned by the present or palaeo stress-field. The upper part of the crust is frequently anisotropic, probably due to cracks differentially aligned by the non-lithostatic stresses. The possibility of being able to monitor crack geometry by seismic techniques opens a wide range of applications in currently important activities.     

Shear-wave polarizations near the North Anatolian Fault – II. Interpretation in terms of crack-induced anisotropy

Summary. The polarizations of shear waves recorded by networks of digital three-component seismometers immediately above small earthquakes near the North Anatolian Fault in Turkey display shear-wave splitting on almost all shear-wave seismograms recorded within the shear-wave window. This splitting is incompatible with source radiation-patterns propagating through simple isotropic structures but is compatible with effective anisotropy of the internal structure of the rock along the ray paths. This paper interprets the phenomena in terms of widespread crack-induced anisotropy. Distributions of stress-induced cracks model many features of the observations, and synthetic polarization diagrams calculated for propagation through simulated cracked rock are similar to the observed patterns. This evidence for widespread crack-induced anisotropy lends strong support to the hypothesis of extensive-dilatancy anisotropy (EDA) suggested by laboratory experiments in subcritical crack-growth. The crucial evidence confirming some form of EDA would be observations of temporal changes in shear-wave splitting as the stress field alters the crack density and crack geometry. There is some weak evidence for such temporal changes at one site, but further analysis of suitable digital three-component seismometer networks in seismic areas is required to confirm EDA.     

The polarization of P-waves in anisotropic media

Summary. The paper examines P -wave propagation in anisotropic solids, and demonstrates the effect of anisotropy on the polarizations of quasi P -waves. The deviation of the polarization of the quasi P -wave from the propagation vector may be significant, but is in almost the same direction as the deviation of the group-velocity vector. Since the group-velocity, or energy propagation, vector follows seismic ray paths, the apparent deviation is the difference of the polarization and group-velocity deviations and is small. Consequently, it may be easily overlooked, hidden by noise, or attributed to the effects of inhomogeneity.     

Surface-wave propagation in an ocean basin with an anisotropic upper mantle: observations of polarization anomalies

Summary. Observations of surface waves crossing ocean basins indicate that, at least in some cases, a better interpretation of what was previously considered to be the simultaneous arrival of fundamental Love and second Rayleigh modes, with similar group velocities, is the coupled motion of the second-generalized mode surface wave propagating in an anisotropic upper-mantle. This mode displays the polarization anomalies expected from anisotropic alignments caused by syntectonic recrystallization of olivine in a zone of horizontal shear. The polarity of the pattern of anomalies, found in waves crossing the Pacific Basin, suggests that the lithosphere is dragging the asthenosphere.     

Effective anisotropic elastic constants for wave propagation through cracked solids

Summary. Theoretical developments of Hudson demonstrate how to calculate the variations of velocity and attenuation of seismic waves propagating through solids containing aligned cracks. The analysis can handle a wide variety of crack configurations and crack geometries. Hudson associates the velocity variations with effective elastic constants. In this paper we associate the variation of attenuation with the imaginary parts of complex effective elastic constants. These complex elastic constants permit the simulation of wave propagation through two-phase materials by the calculation of wave propagation through homogeneous anisotropic solids.     

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.     

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.     

P-SH conversions in a flat-layered medium with anisotropy of arbitrary orientation

P-SH conversion is commonly observed in teleseismic P waves, and is often attributed to dipping interfaces beneath the receiver. Our modelling suggests an alternative explanation in terms of flat-layered anisotropy. We use reflectivity techniques to compute three-component synthetic seismograms in a 1-D anisotropic layered medium. For each layer of the medium, we prescribe values of seismic velocities and hexagonally symmetric anisotropy about a common symmetry axis of arbitrary orientation. A compressional wave in an anisotropic velocity structure suffers conversion to both SV -and SH -polarized shear waves, unless the axis of symmetry is everywhere vertical or the wave travels parallel to all symmetry axes. The P-SV conversion forms the basis of the widely used 'receiver function' technique. The P-SH conversion occurs at interfaces where one or both layers are anisotropic. A tilted axis of symmetry and a dipping interface in isotropic media produce similar amplitudes of both direct ( P ) and converted ( Ps ) phases, leaving the backazimuth variation of the P-Ps delay as the main discriminant. Seismic anisotropy with a tilted symmetry axis leads to complex synthetic seismograms in velocity models composed of just a few flat homogeneous layers. It is possible therefore to model observations of P coda with prominent transverse components with relatively simple 1-D velocity structures. Successful retrieval of salient model characteristics appears possible using multiple realizations of a genetic-algorithm (GA) inversion of P coda from several backazimuths. Using GA inversion, we determine that six P coda recorded at station ARU in central Russia are consistent with models that possess strong (> 10 per cent) anisotropy in the top 5 km and between 30 and 43 km depth. The symmetry axes are tilted, and appear aligned with the seismic anisotropy orientation in the mantle under ARU suggested by SKS splitting.