skip to main content


Title: Love Wave Tomography of the United States
Abstract

Love wave phase velocity maps provide essential constraints on radial anisotropy and deformation in the crust and upper mantle. However, the phenomenon of overtone interference causes scatter and systematic bias in the velocity measurements and impedes efforts to image small‐scale anisotropic variations. We develop an approach for identifying Love wave measurements that are biased by overtone interference, demonstrate its efficacy with EarthScope USArray data, and determine the first earthquake‐derived Love wave phase velocity maps for the entire conterminous U.S. in the period range 35–75 s. We show that radial anisotropy in parts of the crust and most of the lithospheric mantle is necessary to reconcile these maps with Rayleigh wave phase velocities. Our results convey the impact and geographic variability of overtone interference, offer an easy‐to‐implement method to ameliorate this impact, and present high‐resolution constraints on radial anisotropy beneath North America.

 
more » « less
NSF-PAR ID:
10393790
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
49
Issue:
24
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Seismic anisotropy provides essential information for characterizing the orientation of deformation and flow in the crust and mantle. The isotropic structure of the Antarctic crust and upper mantle has been determined by previous studies, but the azimuthal anisotropy structure has only been constrained by mantle core phase (SKS) splitting observations. This study determines the azimuthal anisotropic structure of the crust and mantle beneath the central and West Antarctica based on 8—55 s Rayleigh wave phase velocities from ambient noise cross‐correlation. An anisotropic Rayleigh wave phase velocity map was created using a ray—based tomography method. These data are inverted using a Bayesian Monte Carlo method to obtain an azimuthal anisotropy model with uncertainties. The azimuthal anisotropy structure in most of the study region can be fit by a two‐layer structure, with one layer at depths of 0–15 km in the shallow crust and the other layer in the uppermost mantle. The azimuthal anisotropic layer in the shallow crust of West Antarctica, where it coincides with strong positive radial anisotropy quantified by the previous study, shows a fast direction that is subparallel to the inferred extension direction of the West Antarctic Rift System. Fast directions of upper mantle azimuthal anisotropy generally align with teleseismic shear wave splitting fast directions, suggesting a thin lithosphere or similar lithosphere‐asthenosphere deformation. However, inconsistencies in this exist in Marie Byrd Land, indicating differing ancient deformation patterns in the shallow mantle lithosphere sampled by the surface waves and deformation in the deeper mantle and asthenosphere sampled more strongly by splitting measurements.

     
    more » « less
  2. 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
  3. Abstract

    Lithospheric seismic anisotropy illuminates mid‐ocean ridge dynamics and the thermal evolution of oceanic plates. We utilize short‐period (5–7.5 s) ambient‐noise surface waves and 15‐ to 150‐s Rayleigh waves measured across the NoMelt ocean‐bottom array to invert for the complete radial and azimuthal anisotropy in the upper ∼35 km of ∼70‐Ma Pacific lithospheric mantle, and azimuthal anisotropy through the underlying asthenosphere. Strong azimuthal variations in Rayleigh‐ and Love‐wave velocity are observed, including the first clearly measured Love‐wave 2θand 4θvariations. Inversion of averaged dispersion requires radial anisotropy in the shallow mantle (2‐3%) and the lower crust (4‐5%), with horizontal velocities (VSH) faster than vertical velocities (VSV). Azimuthal anisotropy is strong in the mantle, with 4.5–6% 2θvariation inVSVwith fast propagation parallel to the fossil‐spreading direction (FSD), and 2–2.5% 4θvariation inVSHwith a fast direction 45° from FSD. The relative behavior of 2θ, 4θ, and radial anisotropy in the mantle are consistent with ophiolite petrofabrics, linking outcrop and surface‐wave length scales.VSVremains fast parallel to FSD to ∼80 km depth where the direction changes, suggesting spreading‐dominated deformation at the ridge. The transition at ∼80 km perhaps marks the dehydration boundary and base of the lithosphere. Azimuthal anisotropy strength increases from the Moho to ∼30 km depth, consistent with flow models of passive upwelling at the ridge. Strong azimuthal anisotropy suggests extremely coherent olivine fabric. Weaker radial anisotropy implies slightly nonhorizontal fabric or the presence of alternative (so‐called E‐type) peridotite fabric. Presence of radial anisotropy in the crust suggests subhorizontal layering and/or shearing during crustal accretion.

     
    more » « less
  4. Abstract

    Comprehensive observations of surface wave anisotropy across Alaska and the Aleutian subduction zone would help to improve understanding of its tectonics, mantle dynamics, and earthquake risk. We produce such observations, using stations from the USArray Transportable Array, regional networks across Alaska, and the Alaska Amphibious Community Seismic Experiment in the Alaska‐Aleutian subduction zone both onshore and offshore. Our data include Rayleigh and Love wave phase dispersion from earthquakes (28–85 s) and ambient noise two‐ and three‐station interferometry (8–50 s). Compared with using two‐station interferometry alone, three‐station interferometry significantly improves the signal‐to‐noise ratio and approximately doubles the number of measurements retained. Average differences between both isotropic and anisotropic tomographic maps constructed from different methods lie within their uncertainties, which is justification for combining the measurements. The composite tomographic maps include Rayleigh wave isotropy and azimuthal anisotropy from 8 to 85 s both onshore and offshore, and onshore Love wave isotropy from 8 to 80 s. In the Alaska‐Aleutian subduction zone, Rayleigh wave fast directions vary from trench parallel to perpendicular and back to parallel with increasing periods, apparently reflecting the effect of the subducted Pacific Plate. The tomographic maps provide a basis for inferring the 3‐D anisotropic crustal and uppermost mantle structure across Alaska and the Aleutian subduction zone.

     
    more » « less
  5. null (Ed.)
    Two types of surface wave anisotropy are observed regularly by seismologists but are only rarely interpreted jointly: apparent radial anisotropy, which is the difference in propagation speed between horizontally and vertically polarized waves inferred from Love and Rayleigh waves, and apparent azimuthal anisotropy, which is the directional dependence of surface wave speeds (usually Rayleigh waves). We show that a new data set of Love and Rayleigh wave isotropic phase speeds and Rayleigh wave azimuthal anisotropy observed within and surrounding eastern Tibet can be explained simultaneously by modeling the crust as a depth-dependent tilted hexagonally symmetric (THS) medium. We specify the THS medium with depth-dependent hexagonally symmetric elastic tensors tilted and rotated through dip and strike angles and estimate these quantities using a Bayesian Monte Carlo inversion to produce a 3-D model of the crust and uppermost mantle on a 0.5° × 0.5° spatial grid. In the interior of eastern Tibet and in the Yunnan-Guizhou plateau, we infer a steeply dipping THS upper crustal medium overlying a shallowly dipping THS medium in the middle-to-lower crust. Such vertical stratification of anisotropy may reflect a brittle to ductile transition in which shallow fractures and faults control upper crustal anisotropy and the crystal-preferred orientation of anisotropic (perhaps micaceous) minerals governs the anisotropy of the deeper crust. In contrast, near the periphery of the Tibetan Plateau the anisotropic medium is steeply dipping throughout the entire crust, which may be caused by the reorientation of the symmetry axes of deeper crustal anisotropic minerals as crustal flows are rotated near the borders of Tibet. 
    more » « less