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.     

On crustal corrections in surface wave tomography

Mantle models from surface waves rely on good crustal corrections. We investigated how far ray theoretical and finite frequency approximations can predict crustal corrections for fundamental mode surface waves. Using a spectral element method, we calculated synthetic seismograms in transversely isotropic PREM and in the 3-D crustal model Crust2.0 on top of PREM, and measured the corresponding time-shifts as a function of period. We then applied phase corrections to the PREM seismograms using ray theory and finite frequency theory with exact local phase velocity perturbations from Crust2.0 and looked at the residual time-shifts. After crustal corrections, residuals fall within the uncertainty of measured phase velocities for periods longer than 60 and 80 s for Rayleigh and Love waves, respectively. Rayleigh and Love waves are affected in a highly non-linear way by the crustal type. Oceanic crust affects Love waves stronger, while Rayleigh waves change most in continental crust. As a consequence, we find that the imperfect crustal corrections could have a large impact on our inferences of radial anisotropy. If we want to map anisotropy correctly, we should invert simultaneously for mantle and crust. The latter can only be achieved by using perturbation theory from a good 3-D starting model, or implementing full non-linearity from a 1-D starting model.     

Extension of the Thomson-Haskell method to non-homogeneous spherical layers

Summary. Amplitude spectra of Rayleigh and Love waves in a layered non-gravitating spherical earth have been obtained using as a source, displacement and stress discontinuities. In each layer elastic parameters and density follow specified functions of radial distance and the solutions of the equations of motion are obtained in terms of exponential functions. The Thomson—Haskell method is extended to this case. The problem reduces to simple calculations as in a plane-layered medium. Numerical results of phase and group velocities up to periods of 300 s in various earth models when compared with earlier results (obtained by numerical integration) show that the present method can be used with sufficient accuracy. The differences in phase velocity, group velocity and amplitude (also surface ellipticity in the case of Rayleigh waves) between spherical- and flat-earth models have been investigated in the range 20–300–s period and expressed in polynomials in the period.     

Surface wave tomography of the western United States from ambient seismic noise: Rayleigh and Love wave phase velocity maps

We present the results of Rayleigh wave and Love wave phase velocity tomography in the western United States using ambient seismic noise observed at over 250 broad-band stations from the EarthScope/USArray Transportable Array and regional networks. All available three-component time-series for the 12-month span between 2005 November 1 and 2006 October 31 have been cross-correlated to yield estimated empirical Rayleigh and Love wave Green's functions. The Love wave signals were observed with higher average signal-to-noise ratio (SNR) than Rayleigh wave signals and hence cannot be fully explained by the scattering of Rayleigh waves. Phase velocity dispersion curves for both Rayleigh and Love waves between 5 and 40 speriod were measured for each interstation path by applying frequency–time analysis. The average uncertainty and systematic bias of the measurements are estimated using a method based on analysing thousands of nearly linearly aligned station-triplets. We find that empirical Green's functions can be estimated accurately from the negative time derivative of the symmetric component ambient noise cross-correlation without explicit knowledge of the source distribution. The average traveltime uncertainty is less than 1 s at periods shorter than 24 s. We present Rayleigh and Love wave phase speed maps at periods of 8, 12, 16,and 20 s. The maps show clear correlations with major geological structures and qualitative agreement with previous results based on Rayleigh wave group speeds.     

Global anisotropic phase velocity maps for higher mode Love and Rayleigh waves

It is well established that the Earth's uppermost mantle is anisotropic, but there are no clear observations of anisotropy in the deeper parts of the mantle. Surface waves are well suited to observe anisotropy since they carry information about both radial and azimuthal anisotropy. Fundamental mode surface waves, for commonly used periods up to 200 s, are sensitive to structure in the first few hundred kilometres, and therefore, do not provide information on anisotropy below. Higher mode surface waves have sensitivities that extend to and beyond the transition zone, and should thus give insight about azimuthal anisotropy at greater depths. We have measured higher mode Love and Rayleigh phase velocities using a model space search approach, which provides us with consistent relative uncertainties from measurement to measurement and from mode to mode. From these phase velocity measurements, we constructed global anisotropic phase velocity maps. Prior to inversion, we determine the optimum relative weighting for anisotropy. We present global azimuthal phase velocity maps for higher mode Rayleigh waves (up to the sixth higher mode) and Love waves (up to the fifth higher mode) with corresponding average model uncertainties. The anisotropy we derive is robust within the uncertainties for all modes. Given the ray theoretical sensitivity kernels of Rayleigh and Love wave modes, the source of anisotropy is complex, but mainly located in the asthenosphere and deeper. Our models show a good correspondence with other studies for the fundamental mode, but we have been able to achieve higher resolution.     

Measuring surface-wave overtone phase velocities using a mode-branch stripping technique

We present a method for the retrieval of the phase velocities of surface-wave overtones. The 'single-station' method is successful for several Love and Rayleigh overtone branches (up to at least four) in mode-specific period ranges between 40 and 200 s. It uses mode-branch cross-correlation functions and relies on adjusting the phase and amplitude of the mode branches one at a time. A standard statistical optimization technique is used. We discuss in detail the a priori information that is added to stabilize the retrieval procedure. In addition, we present a technique to estimate the reliability of individual phase and amplitude measurements. The retrieval method and the technique to estimate reliabilities can be used together in a highly automated way, making the methods especially suited for studying the large volume of digital data now available.
We include several applications to synthetic and recorded waveforms. We will discuss in detail an experiment with 90 waveforms that have travelled along very similar paths from Vanuatu to California. For this path, we will present average overtone phase velocities and an average 1-D velocity structure.     

Lithospheric structure of the continent–continent collision zone: eastern Turkey

We infer the lithospheric structure in eastern Turkey using teleseismic and regional events recorded by 29 broad-band stations from the Eastern Turkey Seismic Experiment (ETSE). We combine the surface wave group velocities (Rayleigh and Love) with telesesimic receiver functions to jointly invert for the S -wave velocity structure, Moho depth and mantle-lid (lithospheric mantle) thickness. We also estimated the transverse anisotropy due to Love and Rayleigh velocity discrepancies. We found anomalously low shear wave velocities underneath the Anatolian Plateau. Average crustal thickness is 36 km in the Arabian Plate, 44 km in Anatolian Block and 48 km in the Anatolian Plateau. We observe very low shear wave velocities at the crustal portion (30–38 km) of the northeastern part of the Anatolian Plateau. The lithospheric mantle thickness is either not thick enough to resolve it or it is completely removed underneath the Anatolian Plateau. The shear velocities and anisotropy down to 100 km depth suggest that the average lithosphere–asthenosphere boundary in the Arabian Plate is about 90 and 70 km in Anatolian block. Adding the surface waves to the receiver functions is necessary to constrain the trade-off between velocity and the thickness. We find slower velocities than with the receiver function data alone. The study reveals three different lithospheric structures in eastern Turkey: the Anatolian plateau (east of Karliova Triple Junction), the Anatolian block and the northernmost portion of the Arabian plate. The boundary of lithospheric structure differences coincides with the major tectonic boundaries.     

Surface wave dispersion and Earth structure in south-eastern China

Summary. A reconnaissance study of crust and mantle structure in southeastern China was made using surface waves confined to that region from recent earthquakes. Data from the WWSSN stations ANP arid SEO, and from the digital stations TATO and MAT, were used to measure fundamental-mode group velocities of Love and Rayleigh waves over nine paths in south-eastern China, an area which has been technically quiet since the early Cenozoic. Crustal structure in this region is typical of stable continents, but shear-wave velocities in the uppermost mantle are low for a continent, 4.45 km s⁻¹ or less. Other seismological data support this observation.     

Phase–velocity dispersion curves and small-scale geophysics using noise correlation slantstack technique

It has been demonstrated both theoretically and experimentally that the Green's function between two receivers can be retrieved from the cross-correlation of isotropic noise records. Since surface waves dominate noise records in geophysics, tomographic inversion using noise correlation techniques have been performed from Rayleigh waves so far. However, very few numerical studies implying surface waves have been conducted to confirm the extraction of the true dispersion curves from noise correlation in a complicated soil structure. In this paper, synthetic noise has been generated in a small-scale (<1 km) numerical realistic environment and classical processing techniques are applied to retrieve the phase velocity dispersion curves, first step toward an inversion. We compare results obtained from spatial autocorrelation method (SPAC), high-resolution frequency-wavenumber method (HRFK) and noise correlation slantstack techniques on a 10-sensor array. Two cases are presented in the (1–20 Hz) frequency band that corresponds to an isotropic or a directional noise wavefield. Results show that noise correlation slantstack provides very accurate phase velocity estimates of Rayleigh waves within a wider frequency band than classical techniques and is also suitable for accurately retrieving Love waves dispersion curves.     

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.     

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.     

Off‐great‐circle propagation of intermediate‐period surface waves observed on a dense array in the French Alps

Array analysis is performed on surface waves recorded in the French Alps using a small‐aperture (25 km) temporary array of six broad‐band stations. The analysis shows that both Rayleigh and Love waves deviate relative to the great‐circle path. The deviations are particularly strong, up to 30°, between 20 and 40 s period. To interpret these observations, we first study the effect of large‐scale structures using ray tracing in a smooth, laterally heterogeneous model of the Earth. Second, we evaluate the local effect by considering a model for the French Alps including strong lateral heterogeneities around the array that were not taken into account in the ray tracing. By combining these two possible causes of the observed deviations, we propose an explanation for the general trend in the observed deviations. Finally, we show that by taking into account azimuthal deviations, phase velocities measured at a regional scale can be significantly improved.     

Transitional continental–oceanic structure beneath the Norwegian Sea from inversion of surface wave group velocity data

We have analysed the fundamental mode of Love and Rayleigh waves generated by 12 earthquakes located in the mid-Atlantic ridge and Jan Mayen fracture zone. Using the multiple filter analysis technique, we isolated the Rayleigh and Love wave group velocities for periods between 10 and 50  s. The surface wave propagation paths were divided into five groups, and average group velocities calculated for each group. The average group velocities were inverted and produced shear wave velocity models that correspond to a quasi-continental oceanic structure in the Greenland–Norwegian Sea region. Although resolution is poor at shallow depth, we obtained crustal thickness values of about 18  km in the Norwegian Sea area and 9  km in the region between Svalbard and Iceland. The abnormally thick crust in the Norwegian Sea area is ascribed to magmatic underplating and the thermal blanketing effect of sedimentary layers. Maximum crustal shear velocities vary between 3.5 and 3.9  km  s⁻¹ for most paths. An average lithospheric thickness of 60  km was observed, which is lower than expected for oceanic-type structure of similar age. We also observed low shear wave velocities in the lower crust and upper mantle. We suggest that high heat flow extending to depths of about 30  km beneath the surface can account for the thin lithosphere and observed low velocities. Anisotropy coefficients of 1–5 per cent in the shallow layers and >7 per cent in the upper mantle point to the existence of polarization anisotropy in the region.     

Possible mode conversion between Love and Rayleigh waves at a continental margin

Summary. Mode conversion at a continental margin between normal modes of surface waves is investigated by theoretical calculations for oblique incidence for periods longer than 15 s. It is suggested that significant conversion takes place between the various modes of Love waves in the period interval between 15 and 40 s, while there is negligible mode conversion for longer periods. The largest mode conversion involves the lowest modes. In addition the calculations have revealed a small but significant conversion between Love and Rayleigh fundamental modes around 20-s period. Reflections of Love waves are found to be significant only for the continental fundamental mode.     

A parabolic approximation for surface waves

Summary. A parabolic approximation to the equation of motion of elastic waves as a sum of surface modes and discovering a parabolic approximation be applied directly to surface waves. The approximation depends on the material properties varying slowly within a wavelength, whereas surface waves may travel in a surface wave guide whose depth is of the same order of magnitude as a wavelength. This difficulty is overcome by representing the waves as a sum of surface modes and discovering a parabolic approximation for the amplitudes as a function of position on the surface. The theory is applicable to the propagation of Love or Rayleigh waves in a structure which is vertically stratified in an arbitrary way, but varies slowly in any horizontal direction.     

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.     

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.     

A simple method of representing azimuthal anisotropy on a sphere

We describe a method of expressing azimuthally anisotropic surface wave velocities on the Earth using a local and smooth spherical-spline parametrization. Anisotropy in the Earth leads to azimuthally varying Love and Rayleigh wave velocities that can be expressed as (cos 2ζ, sin 2ζ) and (cos 4ζ, sin 4ζ) perturbations to the isotropic velocities, where ζ is the direction of surface-wave propagation. The strength of the perturbations varies laterally, and a current goal of seismic tomography is the detailed global mapping of these variations. Several parametrizations have previously been used to describe azimuthally varying velocities. The representation proposed here uses spherical splines and is designed to describe smooth variations in both the strength and geometry of azimuthal anisotropy. The method builds on a simple geometrical approximation for the local azimuth of propagation expressed at the defining spline knot points. It avoids the singularities at the poles that result when azimuthal variations are parametrized using traditional scalar spherical harmonics. Compared with a generalized spherical-harmonic expansion of the tensor fields that represent 2ζ and 4ζ azimuthal variations smoothly on a sphere, the new method offers the advantages of local geographical support and simplicity of implementation.     

Gaussian beams for surface waves in laterally slowly-varying media

Summary. Asymptotic ray theory is applied to surface waves in a medium where the lateral variations of structure are very smooth. Using ray-centred coordinates, parabolic equations are obtained for lateral variations while vertical structural variations at a given point are specified by eigenfunctions of normal mode theory as for the laterally homogeneous case. Final results on wavefields close to a ray can be expressed by formulations similar to those for elastic body waves in 2-D laterally heterogeneous media, except that the vertical dependence is described by eigenfunctions of 'local' Love or Rayleigh waves. The transport equation is written in terms of geometrical-ray spreading, group velocity and an energy integral. For the horizontal components there are both principal and additional components to describe the curvature of rays along the surface, as in the case of elastic body waves. The vertical component is decoupled from the horizontal components. With complex parameters the solutions for the dynamic ray tracing system correspond to Gaussian beams: the amplitude distribution is bell-shaped along the direction perpendicular to the ray and the solution is regular everywhere, even at caustics. Most of the characteristics of Gaussian beams for 2-D elastic body waves are also applicable to the surface wave case. At each frequency the solution may be regarded as a set of eigenfunctions propagating over a 2-D surface according to the phase velocity mapping.     

The utility of regional Chinese seismograms for source and path studies in central Asia

Summary. Broadband seismograms from the National Seismic Network of the People's Republic of China (PRC) have recently become available through a data exchange programme between NOAA and the State Seismological Bureau of the PRC. In this study, regional surface waves recorded at the Urumchi station located about 700 km north of the Tibetan Plateau in the Sinkiang Province are used to study East Kazakh explosions and wave propagation in central Asia. The data consist of broadband (flat to displacement between 0.1 and 10 Hz), photographic records from an SK Kirnos galvanometric system. Simultaneous inversion of Rayleigh wave phase and group velocities for the path from East Kazakh through the Dzhungarian Basin yields a crustal model dominated by the presence of very low velocities and a strong positive velocity gradient above 15 km depth. Velocities below 15 km depth are not significantly different from other continental structures underlain by Palaeozoic or Precambrian basement. Seismic moments were estimated for seven East Kazakh explosions using models of explosion sources with associated tectonic strain release. The largest explosion studied occurred on 1980 September 14 and had an m_b of 6.2 and a seismic moment of 2.7 × 10²³ dyn cm. The observed amplitude spectra of Rayleigh waves are richer in high frequencies than spectra calculated from our models. This could be caused by a path effect involving seismic wave focusing by the Dzhungarian Basin, although source medium effects cannot be ruled out.