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.     

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.     

Processing seismic ambient noise data to obtain reliable broad-band surface wave dispersion measurements

Ambient noise tomography is a rapidly emerging field of seismological research. This paper presents the current status of ambient noise data processing as it has developed over the past several years and is intended to explain and justify this development through salient examples. The ambient noise data processing procedure divides into four principal phases: (1) single station data preparation, (2) cross-correlation and temporal stacking, (3) measurement of dispersion curves (performed with frequency–time analysis for both group and phase speeds) and (4) quality control, including error analysis and selection of the acceptable measurements. The procedures that are described herein have been designed not only to deliver reliable measurements, but to be flexible, applicable to a wide variety of observational settings, as well as being fully automated. For an automated data processing procedure, data quality control measures are particularly important to identify and reject bad measurements and compute quality assurance statistics for the accepted measurements. The principal metric on which to base a judgment of quality is stability, the robustness of the measurement to perturbations in the conditions under which it is obtained. Temporal repeatability, in particular, is a significant indicator of reliability and is elevated to a high position in our assessment, as we equate seasonal repeatability with measurement uncertainty. Proxy curves relating observed signal-to-noise ratios to average measurement uncertainties show promise to provide useful expected measurement error estimates in the absence of the long time-series needed for temporal subsetting.     

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.