Surface wave tomography: global membrane waves and adjoint methods

We implement the wave equation on a spherical membrane, with a finite-difference algorithm that accounts for finite-frequency effects in the smooth-Earth approximation, and use the resulting 'membrane waves' as an analogue for surface wave propagation in the Earth. In this formulation, we derive fully numerical 2-D sensitivity kernels for phase anomaly measurements, and employ them in a preliminary tomographic application. To speed up the computation of kernels, so that it is practical to formulate the inverse problem also with respect to a laterally heterogeneous starting model, we calculate them via the adjoint method, based on backpropagation, and parallelize our software on a Linux cluster. Our method is a step forward from ray theory, as it surpasses the inherent infinite-frequency approximation. It differs from analytical Born theory in that it does not involve a far-field approximation, and accounts, in principle, for non-linear effects like multiple scattering and wave front healing. It is much cheaper than the more accurate, fully 3-D numerical solution of the Earth's equations of motion, which has not yet been applied to large-scale tomography. Our tomographic results and trade-off analysis are compatible with those found in the ray- and analytical-Born-theory approaches.     

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.