Youqiang Yu

A uniform database of teleseismic shear wave splitting measurements for the western and central United States

We present a shear wave splitting (SWS) database for the western and central United States as part of a lasting effort to build a uniform SWS database for the entire North America. The SWS measurements were obtained by minimizing the energy on the transverse component of the PKS, SKKS, and SKS phases. Each of the individual measurements was visually checked to ensure quality. This version of the database contains 16,105 pairs of splitting parameters. The data used to generate the parameters were recorded by 1774 digital broadband seismic stations over the period of 1989–2012, and represented all the available data from both permanent and portable seismic networks archived at the Incorporated Research Institutions for Seismology Data Management Center in the area of 26.00°N to 50.00°N and 125.00°W to 90.00°W. About 10,000 pairs of the measurements were from the 1092 USArray Transportable Array stations. The results show that approximately 2/3 of the fast orientations are within 30° from the absolute plate motion (APM) direction of the North American plate, and most of the largest departures with the APM are located along the eastern boundary of the western US orogenic zone and in the central Great Basins. The splitting times observed in the western US are larger than, and those in the central US are comparable with the global average of 1.0 s. The uniform database has an unprecedented spatial coverage and can be used for various investigations of the structure and dynamics of the Earth.

Seismic anisotropy and mantle flow beneath the northern Great Plains of North America

A diverse set of tectonic features and the recent availability of high-quality broadband seismic data from the USArray and other stations on the northern Great Plains of North America provide a distinct opportunity to test different anisotropy-forming mechanisms. A total of 4138 pairs of well-defined splitting parameters observed at 445 stations show systematic spatial variations of anisotropic characteristics. Azimuthally invariant fast orientations subparallel to the absolute plate motion (APM) direction are observed at most of the stations on the Superior Craton and the southern Yavapai province, indicating that a single layer of anisotropy with a horizontal axis of symmetry is sufficient to explain the anisotropic structure. For areas with simple anisotropy, the application of a procedure for estimating the depth of anisotropy using spatial coherency of splitting parameters results in a depth of 200–250 km, suggesting that the observed anisotropy mostly resides in the upper asthenosphere. In the vicinity of the northern boundary of the Yavapai province and the Wyoming Craton, the splitting parameters can be adequately explained by a two-horizontal layer model. The lower layer has an APM-parallel fast orientation, and the upper layer has a fast orientation that is mostly consistent with the regional strike of the boundary. Based on the splitting measurements and previous results from seismic tomography and geodynamic modeling, we propose a model involving deflecting of asthenosphere flow by the bottom of the lithosphere and channeling of flow by a zone of thinned lithosphere approximately along the northern boundary of the Yavapai province.

Determining crustal structure beneath seismic stations overlying a low‐velocity sedimentary layer using receiver functions

The receiver function (RF) technique has been widely applied to investigate crustal and mantle layered structures using P-to-S converted (Ps) phases from velocity discontinuities. However, the presence of low-velocity (relative to that of the bedrock) sediments can give rise to strong reverberations in the resulting RFs, frequently masking the Ps phases from crustal and mantle boundaries. Such reverberations are caused by P-to-S conversions and their multiples associated with the strong impedance contrast across the bottom of the low-velocity sedimentary layer. Here we propose and test an approach to effectively remove the near-surface reverberations and decipher the Ps phases associated with the Moho discontinuity. Autocorrelation is first applied on the observed RFs to determine the strength and two-way traveltime of the reverberations, which are then used to construct a resonance removal filter in the frequency domain to remove or significantly reduce the reverberations. The filtered RFs are time corrected to eliminate the delay effects of the sedimentary layer and applied to estimate the subsediment crustal thickness and VP/VSusing a H-k stacking procedure. The resulting subsediment crustal parameters (thickness and VP/VS) are subsequently used to determine the thickness and VP/VS of the sedimentary layer, using a revised version of the H-k stacking procedure. Testing using both synthetic and real data suggests that this computationally inexpensive technique is efficient in resolving subsediment crustal properties beneath stations sitting on a low-velocity sedimentary layer and can also satisfactorily determine the thickness and VP/VS of the sedimentary layer.

Seismic Anisotropy and Subduction-Induced Mantle Fabrics beneath the Arabian and Nubian Plates Adjacent to the Red Sea

For most continental areas, the mechanisms leading to mantle fabrics responsible for the observed anisotropy remain ambiguous, partially due to the lack of sufficient spatial coverage of reliable seismological observations. Here we report the first joint analysis of shear-wave splitting measurements obtained at stations on the Arabian and Nubian Plates adjacent to the Red Sea. More than 1100 pairs of high-quality splitting parameters show dominantly N-S fast orientations at all 47 stations and larger-than-normal splitting times beneath the Afro-Arabian Dome (AAD). The uniformly N-S fast orientations and large splitting times up to 1.5 s are inconsistent with significant contributions from the lithosphere, which is about 50-80 km thick beneath the AAD and even thinner beneath the Red Sea. The results can best be explained by simple shear between the lithosphere and the asthenosphere associated with northward subduction of the African/Arabian Plates over the past 150 Ma.

Seismic Arrays to Study African Rift Initiation

Rifting of stable continents is a key element of plate tectonic cycles. In spite of numerous studies, the mechanism responsible for the initiation and evolution of rift valleys such as the East African Rift System (EARS) is still poorly understood, partly because most previous investigations focused on rift segments that were in the mature stage. Geodynamic modeling [Huismans et al., 2001] suggests that upwelling of the asthenosphere ubiquitously observed beneath mature rifts can either originate from thermal or dynamic anomalies in the deep mantle (active rifting) or be induced by thinning of the lithosphere from far-field stresses (passive rifting) [Sengor and Burke, 1978].

Azimuthal anisotropy beneath north central Africa from shear wave splitting analyses

This study represents the first multistation investigation of azimuthal anisotropy beneath the interior of north central Africa, including Libya and adjacent regions, using shear wave splitting (SWS) analysis. Data used in the study include recently available broadband seismic data obtained from 15 stations managed by the Libyan Center for Remote Sensing and Space Science, and those from five other stations at which data are publicly accessible. A total of 583 pairs of high-quality SWS measurements utilizing the PKS, SKKS, and SKS phases demonstrate primarily N-S fast orientations with an average splitting delay time of approximately 1.2 s. An absence of periodic azimuthal variation of the observed splitting parameters indicates the presence of simple anisotropy, and lack of correlation between surficial features and the splitting parameters suggests that the origin of the observed anisotropy is primarily asthenospheric. This conclusion is enhanced by nonperiodic azimuthal variation of the splitting parameters observed at one of the stations located near the boundary of areas with different anisotropic properties. We interpret the observed anisotropy to be the consequence of northward movement of the African plate relative to the asthenosphere toward the Hellenic and Calabrian subduction zones. Local variance in fast orientations may be attributable to flow deflection by the northern edge of the African continental root. The observations provide critical and previously lacking constraints on mantle dynamic models in the vicinity of the convergent boundary between the African and Eurasian plates.

Seismic anisotropy and mantle dynamics beneath the Malawi Rift Zone, East Africa†

SKS, SKKS, and PKS splitting parameters measured at 34 seismic stations that we deployed in the vicinity of the Cenozoic Malawi Rift Zone (MRZ) of the East African Rift System demonstrate systematic spatial variations with an average splitting time of 1.0±0.3 s. The overall NE-SW fast orientations are consistent with absolute plate motion (APM) models of the African Plate constructed under the assumption of no-net rotation of the global lithosphere, and are inconsistent with predicted APM directions from models employing a fixed hotspot reference frame. They also depart considerably from the trend of most of the major tectonic features. These observations, together with the results of anisotropy depth estimation using the spatial coherency of the splitting parameters, suggest a mostly asthenospheric origin of the observed azimuthal anisotropy. The single-layered anisotropy observed at 30 and two-layered anisotropy observed at four of the 34 stations can be explained by APM-related simple shear within the rheologically transitional layer between the lithosphere and asthenosphere, as well as by the horizontal deflection of asthenospheric flow along the southern and western edges of a continental block with relatively thick lithosphere revealed by previous seismic tomography and receiver function investigations. This first regional-scale shear wave splitting investigation of the MRZ suggests the absence of rifting-related active mantle upwelling or small-scale mantle convection, and supports a passive-rifting process for the MRZ.

Passive Rifting of Thick Lithosphere in the Southern East African Rift: Evidence from Mantle Transition Zone Discontinuity Topography

To investigate the mechanisms responsible for the initiation and early-stage evolution of the non-volcanic southernmost segments of the East African Rift System (EARS), we installed and operated 35 broadband seismic stations across the Malawi and Luangwa rift zones over a two-year period from mid-2012 to mid-2014. Stacking of over 1,900 high-quality receiver functions provides the first regional-scale image of the 410 and 660-km seismic discontinuities bounding the mantle transition zone (MTZ) within the vicinity of the rift zones. When a 1-D standard Earth model is used for time-depth conversion, a normal MTZ thickness of 250 km is found beneath most of the study area. In addition, the apparent depths of both discontinuities are shallower than normal with a maximum apparent uplift of 20 km, suggesting widespread upper-mantle high-velocity anomalies. These findings suggest that it is unlikely for a low-velocity province to reside within the upper mantle or MTZ beneath the non-volcanic southern EARS. They also support the existence of relatively thick and strong lithosphere corresponding to the widest section of the Malawi rift zone, an observation that is consistent with strain-localization models and fault polarity and geometry observations. We postulate that the Malawi Rift is driven primarily by passive extension within the lithosphere attributed to the divergent rotation of the Rovuma microplate relative to the Nubian plate, and that contributions of thermal upwelling from the lower mantle are insignificant in the initiation and early-stage development of rift zones in southern Africa.

Mantle Structure beneath the Incipient Okavango Rift Zone in Southern Africa

Numerous investigations of the mature segments of the East African rift system (EARS) have significantly improved our understanding of the structure and processes associated with well-developed continental rifts. In contrast, knowledge of rifting processes at their early stage is still significantly limited. Here we present results from a teleseismic P-wave tomography investigation of the incipient Okavango rift zone (ORZ), which is located at the southwestern terminus of the EARS. P-wave relative travel-time residuals recorded by 17 recently deployed portable seismic stations were manually picked and inverted for three-dimensional upper-mantle and mantle transition-zone tomographic images beneath the ORZ and its adjacent areas. High-velocity anomalies probably representing cratonic lithosphere are visible under the Congo and Kalahari cratons, extending to depths of ∼250–350 km. The tectonic boundary of the Congo craton is observed along the western edge of the ORZ. A localized low-velocity anomaly of about –1% in magnitude is revealed in the upper asthenosphere beneath the ORZ, which is interpreted to represent decompression melting induced by lithospheric thinning. The results support the notion that the initiation and early-stage development of the ORZ are mostly due to lithospheric stretching resulted from the relative motion between the Archean Congo and Kalahari cratons along preexisting ancient orogenic zones.

Characteristics of the Mantle Flow System Beneath the Indochina Peninsula Revealed by Teleseismic Shear Wave Splitting Analysis: INDOCHINA SEISMIC ANISOTROPY

State Key Laboratory of Marine Geology, Tongji University, Shanghai, China, Geology and Geophysics Program, Missouri University of Science and Technology, Rolla, Missouri, USA, School of Oceanography, Southern University of Science and Technology, Shenzhen, China, Faculty of Oil and Gas, Hanoi University of Mining and Geology, Hanoi, Vietnam, Second Institute of Oceanography, State Oceanic Administration, Hangzhou, China