The use of equivalent current systems in the interpretation of Geomagnetic Deep Sounding data

Summary. Many geomagnetic variation anomalies are probably caused by the channelling, through small-scale bodies, of electric currents induced in much larger conductors elsewhere. Consequently, the direct interpretation of anomalous magnetic fields by modelling the electromagnetic response of conductive structures may give misleading results. It is suggested that, rather than attempting to proceed directly from the electromagnetic fields to conductivity models, we should instead take the intermediate step of determining the distribution of anomalous current flow.
Maps of the anomalous fields over a conductive structure can be generated from inter-station transfer functions. If it is assumed that the internal currents are concentrated in a thin sheet at a specified depth, the equivalent current system in the sheet can be computed directly from the vertical magnetic field. The most straightforward method of performing this calculation is to compute the Fast Fourier Transform of the magnetic field data, and then to apply a wavenumber filter.
The presence of any vertical currents invalidates the thin sheet model. However, if the spatial distribution of a horizontal component of the anomalous magnetic field is also known, the presence of any vertical currents can be detected directly, and their position determined. The value of the methods is illustrated by applying them to the interpretation of a Geomagnetic Deep Sounding survey of the Kenya rift valley.     

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.     

Inverse problems for torsional modes

Summary. We consider a spherically symmétric, non-rotating earth consisting of an isotropic, perfect elastic material where the density and the S -wave velocity may have one or two discontinuities in the upper mantle. We show that given the velocity throughout the mantle and the crust and given the density in the lower mantle, then the frequencies of the torsional oscillations of one angular order (one torsional spectrum), determine the density in the upper mantle and in the crust uniquely. If the velocity is known only in the lower mantle, then the frequencies of the torsional oscillations of two angular orders uniquely determine both the density and the velocity in the upper mantle and in the crust. In particular, the position and size of the discontinuities in the density and velocity are uniquely determined by two torsional spectra.     

Anisotropy of heterogeneity scale lengths in the lower mantle from PKIKP precursors

The anisotropy of heterogeneity scale lengths in the lower mantle is investigated by modelling its effect on the high-frequency precursors of PKIKP scattered by the heterogeneities. Although models having either an isotropic or an anisotropic distribution of scale lengths can fit the observed coda shapes of short-period precursors, the frequency content of broad-band PKIKP precursors favours a dominantly isotropic distribution of scale lengths. Precursor coda shapes are consistent with 1 per cent fluctuations in P velocity in the wavenumber band 0.05–0.5  km⁻¹ extending to 1000  km above the core–mantle boundary, and with a D" region open to circulation throughout the lower mantle. The level of excitation of PKIKP precursors observed in the frequency band 0.02–2  Hz requires a power spectrum of heterogeneity that is nearly white or slowly increasing with wavenumber. Anisotropy of scale lengths may exist in a D" layer having larger horizontal than vertical scale lengths and produce little or no detectable effects on PKIKP precursors for P -velocity perturbations as high as 3 per cent when averaged over a vertical scale of several kilometres, and much higher when averaged over scales of hundreds of metres or less.     

The significance of isotropic inversion of anisotropic surface-wave dispersion

Summary. The limitations of isotropic modelling in the inversion of anisotropic surface-wave phase velocities are examined. Inversion of synthetic dispersion data for some model ocean-basin structures is used to demonstrate that isotropic inversions can give inaccurate and misleading estimates of upper-mantle properties when anisotropy is present.     

Inversion of P-wave velocity anisotropy

Summary. We examine the way in which measurements of velocity anisotropy can add to our understanding of upper mantle structure. Measurements of P -wave velocity anisotropy in a single plane contain very little direct information about the anisotropic structure. A promising technique is to fit the observed velocity variation with a mixture of an assumed anisotropic constituent and a proportion of isotropic material. Using this technique, mixtures of orthorhombic and transversely isotropic olivine are obtained, which are in excellent agreement with observed velocity variations in the Pacific.     

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.     

Partial derivatives of surface wave phase velocities for flat anisotropic models

Summary. Surface wave behaviour in flat anisotropic structures is first illustrated by performing an exact computation on a simple two-layer model. The variational procedure of Smith & Dahlen is then used to compute the partial derivatives of surface wave phase velocities with respect to the elastic parameters in more realistic earth models. Linear relationships between the partial derivatives for a general anisotropic structure and those for a transversely isotropic structure are derived. When considering waves propagating in a fixed direction, there are only four independent derivatives for Rayleigh waves, and two for Love waves. To avoid the lack of resolution in an inverse method, we propose to use physically constrained models. These results are illustrated by using a model with hexagonal symmetry and a symmetry axis oriented either vertically or horizontally. Quasi-Love- and quasi-Rayleigh-wave partial derivatives are computed for both axis orientations. Modes up to the second overtone and periods ranging between 45 and 130 s have been considered. Finally, anomalies of phase velocity are computed in an oceanic model made of 1/6 oriented olivine crystals with horizontal or vertical preferred orientations of the a -axis.     

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.     

Length scales, patterns and origin of azimuthal seismic anisotropy in the upper mantle as mapped by Rayleigh waves

We measure the degree of consistency between published models of azimuthal seismic anisotropy from surface waves, focusing on Rayleigh wave phase-velocity models. Some models agree up to wavelengths of ∼2000 km, albeit at small values of linear correlation coefficients. Others are, however, not well correlated at all, also with regard to isotropic structure. This points to differences in the underlying data sets and inversion strategies, particularly the relative 'damping' of mapped isotropic versus anisotropic anomalies. Yet, there is more agreement between published models than commonly held, encouraging further analysis. Employing a generalized spherical harmonic representation, we analyse power spectra of orientational (2Ψ) anisotropic heterogeneity from seismology. We find that the anisotropic component of some models is characterized by stronger short-wavelength power than the associated isotropic structure. This spectral signal is consistent with predictions from new geodynamic models, based on olivine texturing in mantle flow. The flow models are also successful in predicting some of the seismologically mapped patterns. We substantiate earlier findings that flow computations significantly outperform models of fast azimuths based on absolute plate velocities. Moreover, further evidence for the importance of active upwellings and downwellings as inferred from seismic tomography is presented. Deterministic estimates of expected anisotropic structure based on mantle flow computations such as ours can help guide future seismologic inversions, particularly in oceanic plate regions. We propose to consider such a priori information when addressing open questions about the averaging properties and resolution of surface and body wave based estimates of anisotropy.     

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.     

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.     

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.     

P-wave traveltime and polarization tomography of VSP data

The presence of anisotropy requires that tomographic methods be generalized to account for anisotropy. This generalization allows geological structure to be correctly imaged and allows the anisotropic parameters to be estimated. Use of isotropic inversion for imaging anisotropic structures gives systematic trends in the traveltime and polarization residuals. However, due to the limited directional coverage, the traveltimes along may not be sufficient to study the anisotropic properties of the structure. Polarizations can provide independent information on the structure. Traveltime and polarization inversion are applied to synthetic examples simulating VSP experiments. Transverse isotropy and 1-D structure are assumed. Plots of traveltime and polarization residuals are an important tool to detect the anomalies due to the presence of anisotropy. For receivers located in anisotropic layers, polarization residuals display consistent anomalies of several degrees. The synthetic examples show that even the simple 1-D problem is difficult, when using direct arrivals only. Large a posteriori errors in anisotropic parameters are obtained by traveltime inversion in layers where available incidence angles are less than 45°. Resolution of the tomographic image of VSP data is greatly improved by a combination of traveltime and polarization information. In order to obtain accurate inversion results, the measurement error of polarization data should be kept to within a few degrees.     

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.     

Formulation of the multiple anisotropic scattering process in two dimensions for anisotropic source radiation

The radiative transfer theory (RTT) describes the energy transport through a random heterogeneous medium, neglecting phase information. It provides an adequate framework for modelling high-frequency seismogram envelopes. For isotropic scattering and sources, the radiative transfer equation (RTE) has been formulated analytically and numerically simulated using Monte Carlo methods for acoustic and elastic media. Here, we derive an exact analytical solution of the RTE in 2-D space for the acoustic case, including anisotropic scattering for a anisotropic point-like impulsive source. For this purpose, we generalize the path integral method, which has been used before in the isotropic case, to take into account the anisotropy of both the source radiation pattern and scattering processes, simultaneously. Then we obtain a general solution, which is written in a closed form in the Fourier space. To illustrate the theoretical results, we compute the full space and time evolution of the specific intensity for an arbitrary case. We also compare the time traces computed from our general solution with cases in which the source and/or the scattering process are isotropic. The importance of taking into account both anisotropies simultaneously becomes obvious in our examples. We also show that at long lapse time, our example approaches the solution of the diffusion equation.     

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.     

Surface-wave polarization data and global anisotropic structure

In the past few years, seismic tomography has begun to provide detailed images of seismic velocity in the Earth's interior which, for the first time, give direct observational constraints on the mechanisms of heat and mass transfer. The study of surface waves has led to quite detailed maps of upper-mantle structure, and the current global models agree reasonably well down to wavelengths of approximately 2000 km. Usually, the models contain only elastic isotropic structure, which provides an excellent fit to the data in most cases. For example, the variance reduction for minor and major arc phase data in the frequency range 7–15 mHz is typically 65–92 per cent and the data are fit to within 1–2 standard deviations. The fit to great-circle phase data, which are not subject to bias from unknown source or instrument effects, is even better. However, there is clear evidence for seismic anisotropy in various places on the globe. This study demonstrates how much (or little) the fit to the data is improved by including anisotropy in the modelling process. It also illuminates some of the trade-offs between isotropic and anisotropic structure and gives an estimate of how much bias is introduced by neglecting anisotropy. Finally, we show that the addition of polarization data has the potential for improving recovery of anisotropic structure by diminishing the trade-offs between isotropic and anisotropic effects.     

A new method for computing synthetic seismograms

Summary. The computation of theoretical seismograms for models in which the elastic parameters and density vary only with depth (in a plane, cylindrical or spherical geometry) reduces to the solution of an ordinary differential equation plus the evaluation of inverse transformations. In principle, the problem is straightforward. In practice, many techniques and approximations can be used at each stage and many combinations and variants are possible. In this paper, we discuss a new method of evaluating the inverse transforms. Any method can be used to solve the differential equation and we only discuss a few analytic approximations to illustrate the new method. The inverse transformations are a frequency and wavenumber integral. Essentially four techniques can be used to evaluate these depending on the order of integration and whether the wavenumber integral is distorted from the real axis. Three of these have been widely used, but the technique of evaluating the frequency integral first and keeping the wavenumber real is new. In this paper, we discuss some of the advantages of this combination.     

Upper mantle anisotropy: evidence from free oscillations

Summary Isotropic earth models are unable to provide uniform fits to the gross Earth normal mode data set or, in many cases, to regional Love-and Rayleigh-wave data. Anisotropic inversion provides a good fit to the data and indicates that the upper 200km of the mantle is anisotropic. The nature and magnitude of the required anisotropy, moreover, is similar to that found in body wave studies and in studies of ultramafic samples from the upper mantle. Pronounced upper mantle low-velocity zones are characteristic of models resulting from isotropic inversion of global or regional data sets. Anisotropic models have more nearly constant velocities in the upper mantle.
Normal mode partial (Frediét) derivatives are calculated for a transversely isotropic earth model with a radial axis of symmetry. For this type of anisotropy there are five elastic constant. The two shear-type moduli can be determined from the toroidal modes. Spheroidal and Rayleigh modes are sensitive to all five elastic constants but are mainly controlled by the two compressional-type moduli, one of the shear-type moduli and the remaining, mixed-mode, modulus. The lack of sensitivity of Rayleigh waves to compressional wave velocities is a characteristic only of the isotropic case. The partial derivatives of the horizontal and vertical components of the compressional velocity are nearly equal and opposite in the region of the mantle where the shear velocity sensitivity is the greatest. The net compressional wave partial derivative, at depth, is therefore very small for isotropic perturbations. Compressional wave anisotropy, however, has a significant effect on Rayleigh-wave dispersion. Once it has been established that transverse anisotropy is important it is necessary to invert for all five elastic constants. If the azimuthal effect has not been averaged out a more general anisotropy may have to be allowed for.