An official website of the United States government
SUMMARY The recent developments in array-based surface-wave tomography have made it possible to directly measure apparent phase velocities through wave front tracking. While directionally dependent measurements have been used to infer intrinsic $$2\psi $$ azimuthal anisotropy (with a 180° periodicity), a few studies have also demonstrated strong but spurious $$1\psi $$ azimuthal anisotropy (360° periodicity) near major structure boundaries particularly for long period surface waves. In such observations, Rayleigh waves propagating in the direction perpendicular to the boundary from the slow to the fast side persistently show a higher apparent velocity compared to waves propagating in the opposite direction. In this study, we conduct numerical and theoretical investigations to explore the effect of scattering on the apparent Rayleigh-wave phase velocity measurement. Using 2-D spectral-element numerical wavefield simulations, we first reproduce the observation that waves propagating in opposite directions show different apparent phase velocities when passing through a major velocity contrast. Based on mode coupling theory and the locked mode approximation, we then investigate the effect of the scattered fundamental-mode Rayleigh wave and body waves interfering with the incident Rayleigh wave separately. We show that scattered fundamental-mode Rayleigh waves, while dominating the scattered wavefield, mostly cause short wavelength apparent phase velocity variations that could only be studied if the station spacing is less than about one tenth of the surface wave wavelength. Scattered body waves, on the other hand, cause longer wavelength velocity variations that correspond to the existing real data observations. Because of the sensitivity of the $$1\psi $$ apparent anisotropy to velocity contrasts, incorporating such measurements in surface wave tomography could improve the resolution and sharpen the structural boundaries of the inverted model.
more »
« less
Global reference seismological data sets: multimode surface wave dispersion
SUMMARY Global variations in the propagation of fundamental-mode and overtone surface waves provide unique constraints on the low-frequency source properties and structure of the Earth’s upper mantle, transition zone and mid mantle. We construct a reference data set of multimode dispersion measurements by reconciling large and diverse catalogues of Love-wave (49.65 million) and Rayleigh-wave dispersion (177.66 million) from eight groups worldwide. The reference data set summarizes measurements of dispersion of fundamental-mode surface waves and up to six overtone branches from 44 871 earthquakes recorded on 12 222 globally distributed seismographic stations. Dispersion curves are specified at a set of reference periods between 25 and 250 s to determine propagation-phase anomalies with respect to a reference Earth model. Our procedures for reconciling data sets include: (1) controlling quality and salvaging missing metadata; (2) identifying discrepant measurements and reasons for discrepancies; (3) equalizing geographic coverage by constructing summary rays for travel-time observations and (4) constructing phase velocity maps at various wavelengths with combination of data types to evaluate inter-dataset consistency. We retrieved missing station and earthquake metadata in several legacy compilations and codified scalable formats to facilitate reproducibility, easy storage and fast input/output on high-performance-computing systems. Outliers can be attributed to cycle skipping, station polarity issues or overtone interference at specific epicentral distances. By assessing inter-dataset consistency across similar paths, we empirically quantified uncertainties in traveltime measurements. More than 95 per cent measurements of fundamental-mode dispersion are internally consistent, but agreement deteriorates for overtones especially branches 5 and 6. Systematic discrepancies between raw phase anomalies from various techniques can be attributed to discrepant theoretical approximations, reference Earth models and processing schemes. Phase-velocity variations yielded by the inversion of the summary data set are highly correlated (R ≥ 0.8) with those from the quality-controlled contributing data sets. Long-wavelength variations in fundamental-mode dispersion (50–100 s) are largely independent of the measurement technique with high correlations extending up to degree ∼25. Agreement degrades with increasing branch number and period; highly correlated structure is found only up to degree ∼10 at longer periods (T > 150 s) and up to degree ∼8 for overtones. Only 2ζ azimuthal variations in phase velocity of fundamental-mode Rayleigh waves were required by the reference data set; maps of 2ζ azimuthal variations are highly consistent between catalogues ( R = 0.6–0.8). Reference data with uncertainties are useful for improving existing measurement techniques, validating models of interior structure, calculating teleseismic data corrections in local or multiscale investigations and developing a 3-D reference Earth model.
more »
« less
- Award ID(s):
- 1758198
- PAR ID:
- 10472031
- Publisher / Repository:
- Wiley
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 228
- Issue:
- 3
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- 1808 to 1849
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We use Eikonal tomography to derive phase and group velocities of surface waves for the plate boundary region in Southern California. Seismic noise data in the period range 2 and 20 s recorded in year 2014 by 346 stations with ~1‐ to 30‐km station spacing are analyzed. Rayleigh and Love wave phase travel times are measured using vertical‐vertical and transverse‐transverse noise cross correlations, and group travel times are derived from the phase measurements. Using the Eikonal equation for each location and period, isotropic phase and group velocities and 2‐psi azimuthal anisotropy are determined statistically with measurements from different virtual sources. Starting with the SCEC Community Velocity Model, the observed 2.5‐ to 16‐s isotropic phase and group dispersion curves are jointly inverted on a 0.05° × 0.05° grid to obtain local 1‐D piecewise shear wave velocity (Vs) models. Compared to the starting model, the final results have generally lowerVsin the shallow crust (top 3–10 km), particularly in areas such as basins and fault zones. The results also show clear velocity contrasts across the San Andreas, San Jacinto, Elsinore, and Garlock Faults and suggest that the San Andreas Fault southeast of San Gorgonio Pass is dipping to the northeast. Investigation of the nonuniqueness of the 1‐DVsinversion suggests that imaging the top 3‐kmVsstructure requires either shorter period (≤2 s) surface wave dispersion measurements or other types of data set such as Rayleigh wave ellipticity.more » « less
-
Summary For a weakly anisotropic medium, Rayleigh and Love wave phase speeds at angular frequency ω and propagation azimuth ψ are given approximately by V(ω, ψ) = A0 + A2ccos 2ψ + A2ssin 2ψ + A4ccos 4ψ + A4ssin 4ψ. Earlier theories of the propagation of surface waves in anisotropic media based on non-degenerate perturbation theory predict that the dominant components are expected to be 2ψ for Rayleigh waves and 4ψ for Love waves. This paper is motivated by recent observations of the the 2ψ component for Love waves and 4ψ for Rayleigh waves, referred to here as “unexpected anisotropy”. To explain these observations, we present a quasi-degenerate theory of Rayleigh-Love coupling in a weakly anisotropic medium based on Hamilton’s Principle in Cartesian coordinates, benchmarking this theory with numerical results based on SPECFEM3D. We show that unexpected anisotropy is expected to be present when Rayleigh-Love coupling is strong and recent observations of Rayleigh and Love wave 2ψ and 4ψ anisotropy can be fit successfully with physically plausible models of a depth-dependent tilted transversely isotropic (TTI) medium. In addition, when observations of the 2ψ and 4ψ components of Rayleigh and Love anisotropy are used in the inversion, the ellipticity parameter ηX, introduced here, is better constrained, we can constrain the absolute dip direction based on polarization measurements, and we provide evidence that the mantle should be modeled as a tilted orthorhombic medium rather than a TTI medium. Ignoring observations of unexpected anisotropy may bias the estimated seismic model significantly. We also provide information about the polarization of the quasi-Love waves and coupling between fundamental mode Love and overtone Rayleigh waves in both continental and oceanic settings. The theory of SV-SH coupling for horizontally propagating body waves is presented for comparison with the surface wave theory, with emphasis on results for a TTI medium.more » « less
-
Abstract We use surface wave measurements to reveal anisotropy as a function of depth within the Juan de Fuca and Gorda plate system. Using a two‐plane wave method, we measure phase velocity and azimuthal anisotropy of fundamental mode Rayleigh waves, solving for anisotropic shear velocity. These surface wave measurements are jointly inverted with constraints fromSKSsplitting studies using a Markov chain approach. We show that the two data sets are consistent and present inversions that offer new constraints on the vertical distribution of strain beneath the plates and the processes at spreading centers. Anisotropy of the Juan de Fuca plate interior is strongest (~2.4%) in the low‐velocity zone between ~40‐ to 90‐km depth, with ENE direction driven by relative shear between plate motion and mantle return flow from the Cascadia subduction zone. In disagreement withPnmeasurements, weak (~1.1%) lithospheric anisotropy in Juan de Fuca is highly oblique to the expected ridge‐perpendicular direction, perhaps connoting complex intralithospheric fabrics associated with melt or off‐axis downwelling. In the Gorda microplate, strong shallow anisotropy (~1.9%) is consistent withPninversions and aligned with spreading and may be enhanced by edge‐driven internal strain. Weak anisotropy with ambiguous orientation in the low‐velocity zone can be explained by Gorda's youth and modest motion relative to the Pacific. Deeper (≥90 km) fabric appears controlled by regional flow fields modulated by the Farallon slab edge: anisotropy is strong (~1.8%) beneath Gorda, but absent beneath the Juan de Fuca, which is in the strain shadow of the slab.more » « less
-
SUMMARY The sensitivity of Rayleigh wave amplitude to Earth structure has applications to seismic tomography, both in cases where amplitude information is used to supplement phase velocity data to improve images of elastic parameters, and to correct amplitudes for local Earth structure in attenuation tomography. We review the theoretical basis of the ray theoretical approximation, in which the wave amplitudes are controlled by a combination of geometrical spreading and local changes in energy density due to Earth structure. We focus mainly on the latter effect, which we term the constant energy flux approximation. We investigate the ray theoretical basis for this approximation, test it against a full waveform simulation that verifies its accuracy and show how it can be used to compute the sensitivity of amplitude to elastic moduli and density. We investigate how perturbing these parameters in a set of simple Earth models affects Rayleigh wave amplitudes, and demonstrate that a slow velocity heterogeneity can cause either increased or reduced amplitudes, depending upon the depth of the heterogeneity and the observation frequency. Consequently, amplitude sensitivity can be either positive or negative, and its magnitude can vary significantly with frequency. Although an added complication, the very different behaviour of phase velocity and amplitudes to changes in Earth structure implies that the two types of data are complementary and suggest the effectiveness of using both in Rayleigh wave tomography.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/gji/ggab418
