Rayleigh Wave Phase Velocity Maps Research Articles

Long Period Rayleigh Wave Focal Spot Imaging Applied to USArray Data

AbstractWe demonstrate the effectiveness of seismic dense array surface wave focal spot imaging using USArray data from the western‐central United States. We study dispersion in the 60–310 s period range and assess the image quality of fundamental mode Rayleigh wave phase velocity maps. We apply isotropic spatial autocorrelation models to the time domain zero lag noise correlation wavefield data at distances of about one wavelength. Local estimates of the phase velocity, its uncertainty, and the regression quality imply overall better ZZ relative to ZR or RZ results. The extension of the depth resolution compared to passive surface wave tomography is demonstrated by the inversion of three clustered dispersion curves from different tectonic units. We observe anisotropic surface wave energy flux and the influence of body wave energy, but sensitivity tests at 60 s targeting the data range, correlation component, and processing choices show that the ZZ focal spots yield consistent high‐quality images compared to regional tomography results in the 60–150 s period range. In contrast, at 200–300 s the comparatively small scales of the imaged structures and the imperfect agreement with low‐resolution global tomography results highlight the persistent challenge to reconcile imaging results based on different data sources, theories, and techniques. Our study shows that surface wave focal spot imaging is an accurate, robust, local imaging approach. Better control over clean autocorrelation fields can further improve applications of this seismic imaging tool for increased resolution of the elastic structure below dense seismic arrays.

Mapping crustal structure across southern Australia using seismic ambient noise tomography

The rocks of southern Australia record over three billion years of Earth’s evolution, but the basement geology is veiled by sediment. Geophysical data are needed to unveil the geology. The 2018–2022 Lake Eyre Basin and 2020–2022 AusArray SA seismic arrays expand seismic coverage in South Australia. In conjunction with permanent and preceding temporary arrays, we extracted Rayleigh wave dispersion data from ambient noise recordings at a total of 501 seismic stations spanning the transition from Precambrian to Phanerozoic Australia. These data were used to develop Rayleigh wave phase velocity maps of southern Australia at periods 3–20s, and in turn, a shear wave velocity model to 20km depth. At upper-crustal depths, low velocity structure tracks Phanerozoic sedimentary accumulations. The Moyston Fault, regarded as the boundary between the Delamerian and Lachlan Orogens, has an intermittent expression in the shear velocity model: it has no obvious expression at depths shallower than ∼10km, but in the mid-crust is marked by a velocity contrast in southern Victoria and a velocity contrast tracing the southern edge of the Darling Basin in western New South Wales. An arcuate velocity contrast characterising the western edge of deep, sediment-filled troughs of the Darling Basin is a candidate for the transition from Precambrian to Phanerozoic crust. Fluid derived from neotectonic metamorphic devolatilization and/or remnant hydrated mantle is our preferred hypothesis for explaining seismicity and coincident seismic and conductivity anomalies in the mid-to-lower crust beneath the Ikara–Flinders Ranges. Our model suggests that the Olympic Dam and Carrapateena IOCG deposits reside above the margin of a mid-crustal low seismic velocity zone. We surmise that this might reflect past metalliferous fluid movement associated with the Olympic Cu-Au Province, akin to low reflectivity and low resistivity zones evident in 2D reflection seismic and magnetotelluric profiles, respectively.

Rayleigh Wave Phase Velocity Maps at Regional Scale Inferring from SPAC of Ambient Noise at a Dense Array: A Case Study in Northeastern Tibetan Plateau

Rayleigh Wave Phase Velocity Maps at Regional Scale Inferring from SPAC of Ambient Noise at a Dense Array: A Case Study in Northeastern Tibetan Plateau

Simultaneous Inversion for Surface Wave Phase Velocity and Earthquake Centroid Parameters: Methodology and Application

AbstractBoth the regional earthquake surface wave and seismic ambient noise provide important constraints on the Earth's structure, and yet no study satisfactorily combined them for the best imaging of subsurface seismic structure. In this study, we address this issue by developing a new method to simultaneously determine surface wave phase velocity and earthquake centroid parameters in three steps: (a) preliminary phase velocity inversion based on seismic ambient noise, (b) preliminary earthquake relocation based on earthquake surface wave data, and (c) simultaneous inversion for phase velocity and earthquake centroid parameters with constraints of interstation phase velocity measurements based on seismic ambient noise and event‐station phase velocity measurements based on earthquake surface wave data. Application of the method in the North South Seismic Belt region in China results in high‐resolution Rayleigh wave phase velocity maps and accurate earthquake centroid parameters. The additional earthquake data notably improve resolution of the inverted phase velocity model in west Yunnan and central Tibetan blocks, the regions with sparse seismic station coverage. The inverted phase velocity model exhibits high‐velocity anomalies in cratonic regions and the Emeishan Large Igneous Province, and low‐velocity anomalies in the interior and surrounding regions of the Tibetan Plateau. Relocation places earthquakes in shallow depths with geotherm above the crustal rock's brittle‐ductile transition temperature of , revealing thermal control on thickness of the seismogenic zone. With earthquake centroid parameters constrained, earthquake data are expected to provide further constraints on the deep seismic structure that is beyond the sampling limit of seismic ambient noise.

Surface-wave tomography using SeisLib: a Python package for multiscale seismic imaging

Summary To improve our understanding of the Earth’s interior, seismologists often have to deal with enormous amounts of data, requiring automatic tools for their analyses. It is the purpose of this study to present SeisLib, an open-source Python package for multiscale seismic imaging. At present, SeisLib includes routines for carrying out surface-wave tomography tasks based on seismic ambient noise and teleseismic earthquakes. We illustrate here these functionalities, both from the theoretical and algorithmic point of view and by application of our library to seismic data from North America. We first show how SeisLib retrieves surface-wave phase velocities from the ambient noise recorded at pairs of receivers, based on the zero crossings of their normalized cross-spectrum. We then present our implementation of the two-station method, to measure phase velocities from pairs of receivers approximately lying on the same great-circle path as the epicentre of distant earthquakes. We apply these methods to calculate dispersion curves across the conterminous United States, using continuous seismograms from the transportable component of USArray and earthquake recordings from the permanent networks. Overall, we measure 144 272 ambient-noise and 2055 earthquake-based dispersion curves, that we invert for Rayleigh-wave phase-velocity maps. To map the lateral variations in surface-wave velocity, SeisLib exploits a least-squares inversion algorithm based on ray theory. Our implementation supports both equal-area and adaptive parametrizations, with the latter allowing for a finer resolution in the areas characterized by high density of measurements. In the broad period range 4–100 s, the retrieved velocity maps of North America are highly correlated (on average, 96 per cent) and present very small average differences (0.14 ± 0.1 per cent) with those reported in the literature. This points to the robustness of our algorithms. We also produce a global phase-velocity map at the period of 40 s, combining our dispersion measurements with those collected at global scale in previous studies. This allows us to demonstrate the reliability and optimized computational speed of SeisLib, even in presence of very large seismic inverse problems and strong variability in the data coverage. The last part of the manuscript deals with the attenuation of Rayleigh waves, which can be estimated through SeisLib based on the seismic ambient noise recorded at dense arrays of receivers. We apply our algorithm to produce an attenuation map of the United States at the period of 4 s, which we find consistent with the relevant literature.

High‐Resolution 3D Shallow S Wave Velocity Structure of Tongzhou, Subcenter of Beijing, Inferred From Multimode Rayleigh Waves by Beamforming Seismic Noise at a Dense Array

AbstractThe 3D S wave velocity of shallow structures, especially the quaternary sediments at 0–1 km near the surface, is an important issue of concern in urban planning and construction for seismic hazard assessment and disaster mitigation. Due to the facility and less dependence on the site environment, noise‐based technique is an ideal way to acquire the fine structure of urban sedimentary basin. Based on the dense array composed of more than 900 stations deployed in Tongzhou at a local scale of 20 × 40 km2, we prove that the lateral variation of the phase velocity of multimode surface waves can be estimated directly with sufficient accuracy by beamforming seismic noise with moving subarrays. Rayleigh wave phase velocity maps, at frequencies between 0.3 and 2.5 Hz for the fundamental mode as well as 0.8 and 3.0 Hz for the first overtone, are obtained. The 3D S wave velocity model at 0–1‐km depth with lateral resolution of 1 km is then established by inverting phase velocity maps of the two modes. The sediment thickness is delineated by the impedance interface given by microtremor H/V (horizontal‐to‐vertical) spectral ratio. The proposed model is in good agreement with the distribution of tectonic units. The sediment thickness of Daxing high and two sags around Gantang and Xiadian are, respectively, 100–400 and 400–600 m, which correlates well with the 1 km/s isosurface of S wave velocity. The model presents some evidence on the extension of Daxing fault along NE direction.

Distributed Extension Across the Ethiopian Rift and Plateau Illuminated by Joint Inversion of Surface Waves and Scattered Body Waves

AbstractThe East African Rift System provides a rare location in which to observe a wide scope of rifting states. Well‐defined active narrow rifting in the Main Ethiopian Rift (MER) transitions to incipient extension and eventually pre‐rifted lithosphere through the northwestern flank of the Ethiopian Plateau (EP). Although the MER is well studied, the off‐axis region has received less attention. We develop Rayleigh wave phase velocity maps, Ps receiver functions, and H‐κ stack surfaces, and jointly invert these data using a trans‐dimensional, hierarchical Bayesian inversion algorithm to create shear velocity profiles across the MER and EP. All shear velocities observed are slower than the PREM global average, a reflection of the elevated temperatures that persist from plume impingement. In the EP, we find a shallow mantle slow shear velocity lineament parallel to the MER axis, amidst otherwise faster shear velocities. The crust is shallow in the MER, and also in the northwestern‐most EP flank. Thicker crust found elsewhere throughout the plateau is caused by crustal underplating and flood basalt emplacement. Shear velocities more reduced than the already low regional average, in concert with surficial volcanic features, geodetic observations, and slow P‐ and S‐wave anomalies, support off‐axis extension in the Ethiopian plateau, requiring reevaluation of the localization of continental breakup in the narrow MER.

Eikonal Tomography Using Coherent Surface Waves Extracted From Ambient Noise by Iterative Matched Filtering—Application to the Large‐N Maupasacq Array

AbstractStandard ambient noise tomography relies on cross‐correlation of noise records between pairs of sensors to estimate empirical Green's functions. This approach is challenging if the distribution of noise sources is heterogeneous and can get computationally intensive for large‐N seismic arrays. Here, we propose an iterative matched filtering method to isolate and extract coherent wave fronts that travel across a dense array of seismic sensors. The method can separate interfering wave trains coming from different directions, to provide amplitude and travel time fields for each detected wave front. We use the eikonal equation to derive phase velocity maps from the gradient of these travel time fields. Artifacts originating from scattered waves are removed by azimuthal averaging and spatial smoothing. The method is validated on a synthetic test and then applied to the data of the Maupasacq experiment. Rayleigh wave phase velocity maps are obtained for periods between 2 and 9 s. These maps correlate with surface geology at short period (T<3 s) and reveal the deep architecture of the Arzacq and Mauleon basins at longer periods (T>4 s).

Surface-Wave Tomography of Eastern and Central Tibet from Two-Plane-Wave Inversion: Rayleigh-Wave and Love-Wave Phase Velocity Maps

ABSTRACT An understanding of mantle dynamics occurring beneath the Tibetan plateau requires a detailed image of its seismic velocity and anisotropic structure. Surface waves at long periods (>50 s) could provide such critical information. Though Rayleigh-wave phase velocity maps have been constructed in the Tibetan regions using ambient-noise tomography (ANT) and regional earthquake surface-wave tomography, Love-wave phase velocity maps, especially those at longer periods (>50 s), are rare. In this study, two-plane-wave teleseismic surface-wave tomography is applied to develop 2D Rayleigh-wave and Love-wave phase velocity maps at periods between 20 and 143 s across eastern and central Tibet and its surroundings using four temporary broadband seismic experiments. These phase velocity maps share similar patterns and show high consistency with those previously obtained from ANT at overlapping periods (20–50 s), whereas our phase velocity maps carry useful information at longer periods (50–143 s). Prominent slow velocity is imaged at periods of 20–143 s beneath the interior of the Tibetan plateau (i.e., the Songpan–Ganzi terrane, the Qiangtang terrane, and the Lhasa terrane), implying the existence of thick Tibetan crust along with warm and weak Tibetan lithosphere. In contrast, the dispersal of fast velocity anomalies coincides with mechanically strong, cold tectonic blocks, such as the Sichuan basin and the Qaidam basin. These phase velocity maps could be used to construct 3D shear-wave velocity and radial seismic anisotropy models of the crust and upper mantle down to 250 km across the eastern and central Tibetan plateau.

3-D crustal structure of the Iran plateau using phase velocity ambient noise tomography

SUMMARY More accurate crustal structure models will help us to better understand the tectonic convergence between Arabian and Eurasian plates in the Iran plateau. In this study, the crustal and uppermost mantle velocity structure of the Iran plateau is investigated using ambient noise tomography. Three years of continuous data are correlated to retrieve Rayleigh wave empirical Green's functions, and phase velocity dispersion curves are extracted using the spectral method. High-resolution Rayleigh wave phase velocity maps are presented at periods of 8–60 s. The tomographic maps show a clear consistency with geological structures such as sedimentary basins and seismotectonic zones, especially at short periods. A quasi-3-D shear wave velocity model is determined from the surface down to 100 km beneath the Iran plateau. A transect of the shear wave velocity model has been considered along with a profile extending across the southern Zagros, the Sanandaj-Sirjan Zone (SSZ), the Urumieh-Dokhtar Magmatic Arc (UDMA) and Central Iran and Kopeh-Dagh (KD). Obvious crustal thinning and thickening are observable along the transect of the shear wave velocity model beneath Central Iran and the SSZ, respectively. The observed shear wave velocities beneath the Iran plateau, specifically Central Iran, support the slab break-off idea in which low density asthenospheric materials drive towards the upper layers, replacing materials in the subcrustal lithosphere.

Amphibious surface-wave phase-velocity measurements of the Cascadia subduction zone

SUMMARY A new amphibious seismic data set from the Cascadia subduction zone is used to characterize the lithosphere structure from the Juan de Fuca ridge to the Cascades backarc. These seismic data are allowing the imaging of an entire tectonic plate from its creation at the ridge through the onset of the subduction to beyond the volcanic arc, along the entire strike of the Cascadia subduction zone. We develop a tilt and compliance correction procedure for ocean-bottom seismometers that employs automated quality control to calculate robust station noise properties. To elucidate crust and upper-mantle structure, we present shoreline-crossing Rayleigh-wave phase-velocity maps for the Cascadia subduction zone, calculated from earthquake data from 20 to 160 s period and from ambient-noise correlations from 9 to 20 s period. We interpret the phase-velocity maps in terms of the tectonics associated with the Juan de Fuca plate history and the Cascadia subduction system. We find that thermal oceanic plate cooling models cannot explain velocity anomalies observed beneath the Juan de Fuca plate. Instead, they may be explained by a ≤1 per cent partial melt region beneath the ridge and are spatially collocated with patches of hydration and increased faulting in the crust and upper mantle near the deformation front. In the forearc, slow velocities appear to be more prevalent in areas that experienced high slip in past Cascadia megathrust earthquakes and generally occur updip of the highest-density tremor regions and locations of intraplate earthquakes. Beneath the volcanic arc, the slowest phase velocities correlate with regions of highest magma production volume.

3D shallow structures in the Baogutu area, Karamay, determined by eikonal tomography of short-period ambient noise surface waves

Eikonal tomography based on ambient noise data is one of the most effective methods to reveal shallow earth structures. By tracking surface wave phase fronts, constructing travel time surfaces, and computing the gradients of travel time surfaces to generate phase velocity maps, eikonal tomography avoids the ray tracing and matrix construction and inversion in the traditional surface wave tomography methods.In this study, we collect continuous ambient noise data recorded by a dense seismic array in Karamay, Xinjiang to construct a 3D model of shallow structures using eikonal tomography. The seismic array consists of 35 stations with shortest interstation distance close to 1km. 890 empirical surface wave Green's functions (EGFs) between each station pair are retrieved by cross-correlating one or two months of continuous ambient noise data. From these EGFs, surface wave travel times in the frequency range of 1.8 to 4.0 Hz are measured by a frequency–time analysis technique (FTAN). Then, eikonal tomography is adopted to construct Rayleigh wave phase velocity maps and estimate the phase velocity uncertainties. Finally, we invert the obtained phase velocity dispersion curves for 1D shear velocity profiles and then assemble these 1D profiles to construct a 3D shear velocity model. Major velocity features of our 3D model are correlated well with the known geological features. A shallow east–west velocity discontinuity is observed, which clearly reflects the lithological change between Baogutu formation (C1b) and Xibeikulasi formation (C1x) of lower Carboniferous system. Low shear velocities are observed beneath the location of porphyry copper deposit (V), possibly related to stockwork fracture and hydrothermal brecciation developed during the intrusion of deep magma in forming the deposit.

Shear-wave velocity structure of the crust and uppermost mantle in the Shanxi rift zone

The Shanxi rift zone is one of the largest and active Cenozoic grabens in the world, studying the velocity structure of the crust and upper mantle in this region may help us to understand the mechanisms of rift processes and the seismogenic environment of active seismicity in continental rifts. In this work, using the broadband seismic data of Shanxi, Hebei, Henan, Shaanxi provinces, and the Inner Mongolia Autonomous Region from February 2009 to November 2011, we have picked out 350 high-quality phase velocity dispersion curves of fundamental mode Rayleigh waves at periods from 8 to 75 s, and Rayleigh wave phase velocity maps have been constructed from 8 to 75 s period with horizontal resolution ranging from 40 to 50 km by two-station surface-wave tomography. Then, using a genetic algorithm, a 3D shear-wave speed model of the crust and uppermost mantle have been derived from these maps with a spatial resolution of 0.4° × 0.4°. Four characteristics can be outlined from the results: (1) Except in the Datong volcanic zone, in the depth range of 11–30 km, the location of a transition zone between the high- and low-velocity regions is in agreement with the seismicity pattern in the study region, and the earthquakes are mostly concentrated near this transition zone; (2) In the depth range of 31–40 km, shear-wave velocities are higher to the south of the Taiyuan Basin and lower to the north, which is similar to the distribution pattern of Moho depth variations in the Shanxi region; (3) The shear-wave velocity pattern of higher velocities to the south of 38°N and lower velocities to the north is found to be consistent with that from the upper crustal levels to depth of 70 km. At the deeper depths, the spatial scale of the low-velocity anomalies zone in the north is gradually shrinking with depth increasing, the low-velocity anomalies are gradually disappearing beneath the Datong volcanic zone at the depth of 151–200 km. We proposed that the root of the Datong volcano may reach to a depth around 150 km; (4) Along the N–S vertical profile at 112.8°E, the 38°N latitude is the boundary between high and low velocities, arguing the tectonic difference between the Shanxi rift zone and its flanks, in the rift zone the seismic velocity is dominated by low-velocity anomalies while in the flanks it is high.

Anisotropic Rayleigh-wave phase velocities beneath northern Vietnam

We explore the Rayleigh-wave phase-velocity structure beneath northern Vietnam over a broad period range of 5 to 250 s. We use the two-stations technique to derive the dispersion curves from the waveforms of 798 teleseismic events recoded by a set of 23 broadband seismic stations deployed in northern Vietnam. These dispersion curves are then inverted for both isotropic and azimuthally anisotropic Rayleigh-wave phase-velocity maps in the frequency range of 10 to 50 s. Main findings include a crustal expression of the Red River Shear Zone and the Song Ma Fault. Northern Vietnam displays a northeast/southwest dichotomy in the lithosphere with fast velocities beneath the South China Block and slow velocities beneath the Simao Block and between the Red River Fault and the Song Da Fault. The anisotropy in the region is relatively simple, with a high amplitude and fast directions parallel to the Red River Shear Zone in the western part. In the eastern part, the amplitudes are generally smaller and the fast axis displays more variations with periods.

Seismic anisotropy changes across upper mantle phase transitions

The mantle transition zone is believed to play an important role in the thermochemical evolution of our planet and in its deep water cycle. Constraining mantle flow at these depths can help elucidate its nature and better understand mantle dynamics and the history of plate tectonics. Seismic anisotropy, i.e., the directional dependence of seismic wave velocity, provides us with the most direct constraints on mantle deformation. Its detection below ∼250km depth is challenging, and it is often assumed that the deep upper mantle is seismically isotropic due to a change in mantle deformation mechanism. Here, we present a global model of azimuthal anisotropy for the top 1000km of the mantle. We used a dataset composed of fundamental and higher mode Rayleigh wave phase velocity maps, which provides resolution of azimuthal anisotropy to much greater depths than in previous studies. Our model unravels the presence of significant anisotropy in the transition zone, challenging common views of mantle deformation mechanisms, and reveals a striking correlation between changes in anisotropy amplitudes and in the fast direction of wave propagation where upper mantle phase transitions occur. The newly found relation between anisotropy changes and phase transformations gives new insight on the nature of the MTZ and suggests that the anisotropy originates from the lattice preferred orientation of anisotropic material. While the interpretation of our results in terms of mantle deformation is not straightforward due to many possible scenarios, possible explanations for our findings include recrystallization during phase transformations, or changes in the slip system across the MTZ boundaries, which in turn could be explained by changes in the water content of mantle material, consistent with the idea that the transition zone acts as a water filter.

A synoptic view of the distribution and connectivity of the mid‐crustal low velocity zone beneath Tibet

Based on 1–2 years of continuous observations of seismic ambient noise data obtained at more than 600 stations in and around Tibet, Rayleigh wave phase velocity maps are constructed from 10 s to 60 s period. A 3‐D Vsv model of the crust and uppermost mantle is derived from these maps. The 3‐D model exhibits significant apparently inter‐connected low shear velocity features across most of the Tibetan middle crust at depths between 20 and 40 km. These low velocity zones (LVZs) do not conform to surface faults and, significantly, are most prominent near the periphery of Tibet. The observations support the internal deformation model in which strain is dispersed in the deeper crust into broad ductile shear zones, rather than being localized horizontally near the edges of rigid blocks. The prominent LVZs are coincident with strong mid‐crustal radial anisotropy in western and central Tibet and probably result at least partially from anisotropic minerals aligned by deformation, which mitigates the need to invoke partial melt to explain the observations. Irrespective of their cause in partial melt or mineral alignment, mid‐crustal LVZs reflect deformation and their amplification near the periphery of Tibet provides new information about the mode of deformation across Tibet.

Crustal structure beneath the Dabie orogenic belt from ambient noise tomography

The Qinling–Dabie–Sulu orogenic belt in east-central China is the largest high and ultrahigh pressure (HP and UHP) metamorphic zone in the world. The Dabie Mountains are the central segment of this orogenic belt between the North China and Yangtze cratons. This work studies the nature of the crustal structure beneath the Dabie orogenic belt to better understand the orogeny. To do that, we apply ambient noise tomography to the Dabie orogenic belt using ambient noise data from 40 stations of the China National Seismic Network (CNSN) between January 2008 and December 2009. We retrieve high signal noise ratio (SNR) Rayleigh waves by cross-correlating ambient noise data between most of the station pairs and then extract phase velocity dispersion measurements from those cross-correlations using a spectral method. Taking those dispersion measurements, we obtain high-resolution phase velocity maps at 8–35 second periods. By inverting Rayleigh wave phase velocity maps, we construct a high-resolution 3D shear velocity model of the crust in the Dabie orogenic belt. The resulting 3D model reveals interesting crustal features related to the orogeny. High shear wave velocities are imaged beneath the HP/UHP metaphoric zones at depths shallower than 9 km, suggesting that HP/UHP metaphoric rocks are primarily concentrated in the upper crust. Underlying the high velocity HP/UHP metamorphic zones, low shear velocities are observed in the middle crust, probably representing ductile shear zones and/or brittle fracture zones developed during the exhumation of the HP/UHP metamorphic rocks. Strong high velocities are present beneath the Northern Dabie complex unit in the middle crust, possibly related to cooling and crystallization of intrusive igneous rocks in the middle crust resulting from the post-collisional lithosphere delamination and subsequent magmatism. A north-dipping Moho is revealed in the eastern Dabie with the deepest Moho appearing beneath the Northern Dabie complex unit, consistent with the model of Triassic northward subduction of the Yangtze Craton beneath the North China Craton.

Rayleigh wave phase velocity maps of Tibet and the surrounding regions from ambient seismic noise tomography

Ambient noise tomography is applied to the significant data resources now available across Tibet and surrounding regions to produce Rayleigh wave phase speed maps at periods between 6 and 50 s. Data resources include the permanent Federation of Digital Seismographic Networks, five temporary U.S. Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) experiments in and around Tibet, and Chinese provincial networks surrounding Tibet from 2003 to 2009, totaling ∼600 stations and ∼150,000 interstation paths. With such a heterogeneous data set, data quality control is of utmost importance. We apply conservative data quality control criteria to accept between ∼5000 and ∼45,000 measurements as a function of period, which produce a lateral resolution between 100 and 200 km across most of the Tibetan Plateau and adjacent regions to the east. Misfits to the accepted measurements among PASSCAL stations and among Chinese stations are similar, with a standard deviation of ∼1.7 s, which indicates that the final dispersion measurements from Chinese and PASSCAL stations are of similar quality. Phase velocities across the Tibetan Plateau are lower, on average, than those in the surrounding nonbasin regions. Phase velocities in northern Tibet are lower than those in southern Tibet, perhaps implying different spatial and temporal variations in the way the high elevations of the plateau are created and maintained. At short periods (<20 s), very low phase velocities are imaged in the major basins, including the Tarim, Qaidam, Junggar, and Sichuan basins, and in the Ordos Block. At intermediate and long periods (>20 s), very high velocities are imaged in the Tarim Basin, the Ordos Block, and the Sichuan Basin. These phase velocity dispersion maps provide information needed to construct a 3‐D shear velocity model of the crust across the Tibetan Plateau and surrounding regions.

Rayleigh-wave phase-velocity maps and three-dimensional shear velocity structure of the western US from local non-plane surface wave tomography

SUMMARY We utilize two-and-three-quarter years of vertical-component recordings made by the Transportable Array (TA) component of Earthscope to constrain three-dimensional (3-D) seismic shear wave velocity structure in the upper 200 km of the western United States. Single-taper spectral estimation is used to compile measurements of complex spectral amplitudes from 44 317 seismograms generated by 123 teleseismic events. In the first step employed to determine the Rayleigh-wave phase-velocity structure, we implement a new tomographic method, which is simpler and more robust than scattering-based methods (e.g. multi-plane surface wave tomography). The TA is effectively implemented as a large number of local arrays by defining a horizontal Gaussian smoothing distance that weights observations near a given target point. The complex spectral-amplitude measurements are interpreted with the spherical Helmholtz equation using local observations about a succession of target points, resulting in Rayleigh-wave phase-velocity maps at periods over the range of 18–125 s. The derived maps depend on the form of local fits to the Helmholtz equation, which generally involve the non-plane-wave solutions of Friederich et al. In a second step, the phase-velocity maps are used to derive 3-D shear velocity structure. The 3-D velocity images confirm details witnessed in prior body-wave and surface-wave studies and reveal new structures, including a deep (>100 km deep) high-velocity lineament, of width ∼200 km, stretching from the southern Great Valley to northern Utah that may be a relic of plate subduction or, alternatively, either a remnant of the Mojave Precambrian Province or a mantle downwelling. Mantle seismic velocity is highly correlated with heat flow, Holocene volcanism, elastic plate thickness and seismicity. This suggests that shallow mantle structure provides the heat source for associated magmatism, as well as thinning of the thermal lithosphere, leading to relatively high stress concentration. Our images also confirm the presence of high-velocity mantle at ≳100 km depth beneath areas of suspected mantle delamination (southern Sierra Nevada; Grande Ronde uplift), low velocity mantle underlying active rift zones, and high velocity mantle associated with the subducting Juan de Fuca plate. Structure established during the Proterozoic appears to exert a lasting influence on subsequent volcanism and tectonism up to the Present.

Depth constraints on azimuthal anisotropy in the Great Basin from Rayleigh-wave phase velocity maps

We present fundamental-mode Rayleigh-wave azimuthally anisotropic phase velocity maps obtained for the Great Basin region at periods between 16 s and 102 s. These maps offer the first depth constraints on the origin of the semi-circular shear-wave splitting pattern observed in central Nevada, around a weak azimuthal anisotropy zone. A variety of explanations have been proposed to explain this signal, including an upwelling, toroidal mantle flow around a slab, lithospheric drip, and a megadetachment, but no consensus has been reached. Our phase velocity study helps constrain the three-dimensional anisotropic structure of the upper mantle in this region and contributes to a better understanding of the deformation mechanisms taking place beneath the western United States. The dispersion measurements were made using data from the USArray Transportable Array. At periods of 16 s and 18 s, which mostly sample the crust, we find a region of low anisotropy in central Nevada coinciding with locally reduced phase velocities, and surrounded by a semi-circular pattern of fast seismic directions. Away from central Nevada the fast directions are ∼ N–S in the eastern Great Basin, NW–SE in the Walker Lane region, and they transition from E–W to N–S in the northwestern Great Basin. Our short-period phase velocity maps, combined with recent crustal receiver function results, are consistent with the presence of a semi-circular anisotropy signal in the lithosphere in the vicinity of a locally thick crust. At longer periods (28–102 s), which sample the uppermost mantle, isotropic phase velocities are significantly reduced across the study region, and fast directions are more uniform with an ∼ E–W fast axis. The transition in phase velocities and anisotropy can be attributed to the lithosphere–asthenosphere boundary at depths of ∼ 60 km. We interpret the fast seismic directions observed at longer periods in terms of present-day asthenospheric flow-driven deformation, possibly related to a combination of Juan de Fuca slab rollback and eastward-driven mantle flow from the Pacific asthenosphere. Our results also provide context to regional SKS splitting observations. We find that our short-period phase velocity anisotropy can only explain ∼ 30% of the SKS splitting times, despite similar patterns in fast directions. This implies that the origin of the regional shear-wave splitting signal is complex and must also have a significant sublithospheric component.