Receiver Functions Research Articles

A lower crust shear zone facilitates delamination and continental subduction under the Apennines

Physical properties and structure of the lithosphere are the first step to constrain the evolution of mountain belts. Here we show detailed shear wave velocity profiles of the lithosphere in the Apennines that clarify a controversial aspect of continental subduction: the intricate mechanism of crust delamination from the downgoing plate. From the analysis of complete and dense teleseismic Receiver Function data set, we find that the delamination of the continental lithosphere is favored by the development of a low seismic shear wave velocity zone in the middle-lower crust. We observe a double Moho below the external portions of the present mountain range, suggesting the progressive formation of the shallow interface. The delamination edge is located in the forearc, far eastward than expected, implying that the re-equilibration of the thermal unbalance, generated by the mantle substitution, may last 10-7 Myr.

Underthrusting of Tarim Lower Crust Beneath the Tibetan Plateau Revealed by Receiver Function Imaging

AbstractThe left‐lateral Altyn Tagh Fault (ATF) system is the northern boundary of the Tibetan Plateau resulted from the India–Eurasia continental collision. How intracontinental deformation across the central ATF responds to the distal collision remains elusive, primarily due to unclear crustal structure. We obtained detailed crustal structure across the central ATF using receiver functions recorded by ∼NW–SE oriented linear dense array. The images reveal the Tarim lower crust is underthrusting beneath the Tibetan Plateau and reaches to a maximum depth of ∼75 km and undergoing partial eclogitization. The two south‐dipping interfaces imaged beneath the Altyn Tagh Range (ATR) represent the thrusting Northern Altyn Fault and its branch fault. Oblique convergent forces extruded upper crustal materials along the thrust faults, creating the pop‐up structure of ATR, supported by low Vp/Vs ratios. Our balanced cross‐section for the Moho suggests intracontinental deformation in the ATR has accelerated since the late Miocene.

Low-strength shear zone in the western Makran subduction zone, southeastern Iran: insights from a receiver function analysis

SUMMARY To understand the seismic hazard of a subduction zone, it is necessary to know the geometry, location and mechanical characteristics of the interplate boundary below which an oceanic plate is thrust downward. By considering the azimuthal dependence of converted P-to-S (Ps) amplitudes in receiver functions, we have detected the interplate boundary in the Makran subduction zone, revealing significant seismic anisotropy at the base of the accretionary wedge above the slab before it bends down beneath the Jaz Murian basin. This anisotropic feature aligns with a zone of reduced seismic velocity and a high primary/secondary wave velocity ratio (Vp/Vs), as documented in previous studies. The presence of this low-velocity highly anisotropic layer at the base of the accretionary wedge, likely representing a low-strength shear zone, could possibly explain the unusually wide accretionary wedge in Makran. Additionally, it may impact the location and width of the locked zone along the interplate boundary.

Sensor Misorientation of CHINArray-II in Northeastern Tibetan Plateau from P- and Rayleigh-Wave Polarization Analysis

Abstract Accurate sensor orientation is important in providing reliable data used for seismological analysis such as P-wave receiver function analysis, shear-wave splitting, and ambient noise analysis. In this study, we used three distinct P-wave and Rayleigh-wave polarization analysis methods to estimate actual sensor orientation of 660 stations from CHINArray-II in northeastern margin of Tibetan plateau. We found that ∼42.12%–45.76% of the stations are well oriented with the absolute misorientation angle <3°, 42.42%–45.45% of the stations are fairly oriented with the absolute misorientation angle ranging from 3° to 10°, and 8.94%–11.82% of the stations are oriented with the absolute misorientation angle >10°. We further compared the results of some seismological analyses before and after sensor misorientation correction, such as Rayleigh-wave ellipticity (horizontal-to-vertical ratio) and P-wave receiver functions analysis. We found that when the sensor misorientation angle is large, it may lead to incorrect seismological results. With the same sensor misorientation, its influence on different seismological analysis is also different.

Geophysical Evidence of the Collisional Suture Zone in the Prydz Bay, East Antarctica

AbstractThe location and origin of Neoproterozoic‐Cambrian sutures provide keys to understand the formation and evolution of the supercontinent Gondwana. The Larsemann Hills is located near a major Neoproterozoic‐Cambrian suture zone in the Prydz Belt, but has not been examined locally by comprehensive geophysical studies. In this study, we analyzed data collected from a one‐dimensional (1D) joint seismic‐MT array deployed during the 36th Chinese National Antarctic Research Expedition. We found that a sharp Moho discontinuity offset of 6–8 km shows up in the stacked image of teleseismic P‐wave receiver function analysis; coinciding with the abrupt Moho offset, a near‐vertical channel with (a) low resistivity extending to the uppermost mantle depths, and (b) high crustal Poisson's ratio in the crust is identified. These findings provide evidence for the determination of the location and collisional nature of the Prydz belt or a portion of it.

Possible layered lithospheric anisotropy around Longmenshan Faults by teleseismic S wave splitting and receiver functions

This study conducts an in-depth analysis of seismic anisotropy around the Longmenshan Faults. Utilizing a dataset of about 7710 earthquake catalogs from July 2007 to March 2023, we applied S wave splitting and receiver function methods to examine Pms and XKS waveforms collected from 12 fixed broadband stations across Gansu and Sichuan provinces. Our analysis revealed significant variations in seismic anisotropy between the crust and lithosphere, marked by distinct fast wave directions and delay times. These characteristics point to the possibility for layered anisotropy within the region. A two-layer anisotropy inversion analysis at key stations further delineated the differential anisotropic behaviors between the crust and the upper mantle, underscoring the impact of local geological structures and mantle dynamics. Crucially, our study posits the existence of layered anisotropy around the Longmenshan Fault Zone, a finding that significantly advances our comprehension of the region’s seismic anisotropy and adds a vital dimension to our understanding of its subsurface structural intricacies and tectonic evolution.

Crustal Imaging with Noisy Teleseismic Receiver Functions Using Sparse Radon Transforms

ABSTRACT The receiver function (RF) is a widely used crustal imaging technique. In principle, it assumes relatively noise-free traces that can be used to target receiver-side structures following source deconvolution. In practice, however, mode conversions and reflections may be severely degraded by noisy conditions, hampering robust estimation of crustal parameters. In this study, we use a sparsity-promoting Radon transform to decompose the observed RF traces into their wavefield contributions, that is, direct conversions, multiples, and incoherent noise. By applying a crustal mask on the Radon-transformed RF, we obtain noise-free RF traces with only Moho conversions and reflections. We demonstrate, using a synthetic experiment and a real-data example from the Sierra Nevada, that our approach can effectively denoise the RFs and extract the underlying Moho signals. This greatly improves the robustness of crustal structure recovery as exemplified by subsequent H−κ stacking. We further demonstrate, using a station sitting on loose sediments in the Upper Mississippi embayment, that a combination of our approach and frequency-domain filtering can significantly improve crustal imaging in reverberant settings. In the presence of complex crustal structures, for example, dipping Moho, intracrustal layers, and crustal anisotropy, we recommend caution when applying our proposed approach due to the difficulty of interpreting a possibly more complicated Radon image. We expect that our technique will enable high-resolution crustal imaging and inspire more applications of Radon transforms in seismic signal processing.

Seismic Evidence for a Transcrustal Magmatic Pathway Contributing to Critical Metal Deposits

AbstractThe Nanling metallogenic belt in South China hosts ∼50 percent of the world's tungsten (W), ∼20 percent tin (Sn), and other critical metals. However, the factors controlling ore distribution and types are still debated, due partly to the lack of high‐resolution crustal models. We apply an array‐based receiver function imaging method to resolve the fine‐scale crustal structure using a newly deployed dense nodal array across the belt. We find weak Moho converting phases in the northern belt, indicating a gradual crust‐mantle transition distinct from the sharp transitions observed elsewhere. The weak Moho spatially correlates with a major fault zone, A‐type granite distribution, and Sn ore deposits. Combining with geochemical and petrological evidence, we propose that the gradual crust‐mantle transition is attributed to, and likely facilitated, the mixing of mantle materials with the crust, thus contributing to Sn and other Mesozoic ore formations in the northern Nanling belt.

Crustal Structures From Receiver Functions and Gravity Modeling in Central Mongolia

Abstract3D forward gravity modeling combined with receiver function (RF) analysis characterizes the crustal structures of the southern part of the Mongolian collage. The seismic signals of the 48 stations of the MOBAL2003 and the IRIS‐PASSCAL experiments were analyzed to get the RFs. This analysis revealed a significant difference between the crustal structures of the Hangay dome and the tectonic zones in the south. In addition, seismic stations south of the Hangay dome display significant signals related to the occurrence of a low‐velocity zone at lower crustal level confirmed by the gravity anomalies. Finally, these seismic analysis inputs, the boundaries, the lithologies, and the density values from rock samples of the different tectonic zones constitute the starting points from the 3D forward gravity modeling. The resulting crustal density model indicates: (a) the likely absence of a Precambrian basement block beneath the Hangay dome, (b) an alternation of two low‐velocity/low‐density zones (LVLDZs) with high‐density zones in the Baydrag microcontinent interpreted as fragments of early Tonian plutons, (c) the occurrence of an LVLDZ at the lower crustal level beneath the Lake zone, the Mongol‐Altai Accretionary Wedge, and the Trans‐Altai Zone. Therefore, the combination of the seismic RF with gravity analysis and modeling reveals new crustal structures of the Mongolian collage and enhances the occurrence and the extent of an LVLDZ at lower crustal level. These LVLDZ may demonstrate the existence of the relamination of a hydrous material in southern Mongolian collage.

Characterization of anisotropy in basin-scale subsurface using teleseismic receiver function analysis

Teleseismic receiver function (RF) analysis offers a passive-source analogue to detect impedance boundaries using the converted body waves generated by earthquakes. While the technique traditionally has targeted deep earth structures such as the Moho and transition zones, there is growing interest in assessing its applicability in basin-scale seismic characterization, ultimately aimed for onshore commercial integration as a cost-effective complement to existing active-source seismic surveys. Specifically, the conventional broadband seismometers used in global observational seismic experiments are not only logistically simple in terms of data acquisition, but they also record ground motions in three mutually orthogonal time series, enabling effective detection of shear waves and directional variations of observed signals. Here, we perform teleseismic RF analysis to detect shear-wave anisotropy and related symmetry axes orientations in a basin setting, using open-source seismic data recorded at 55 closely spaced seismic stations in the LaBarge Passive Seismic Experiment deployed in Wyoming between November 2008 and June 2009. We find that the strengths and geometry of the observed anisotropy are variable along the array. Significantly, not only can anisotropy effectively delineate subsurface fault interfaces, it can also substantiate and reveal additional interpretable signals that are otherwise disregarded. The estimated fast axes orientations compare favorably with the complex fracture systems documented in the region. Finally, we show that P-to-S amplitude variations with P incidence are systematic and modelable using existing computational tools, offering an opportunity to develop an analysis technique similar to amplitude variation with offset with the products of RF analysis.

Africa's Crustal Architecture Inferred From Probabilistic and Perturbational Inversion of Ambient Noise: ADAMA

AbstractAfrica's continental crust hosts a variety of geologic terrains and is crucial for understanding the evolution of its longest‐lived cratons. However, few of its seismological models are yet to incorporate the largest continent‐wide noise dispersion data sets. Here, we report on new insights into Africa's crustal architecture obtained using a new data set and model assessment product, ADAMA, which comprises a large ensemble of short‐period surface wave dispersion measurements: 5–40 s. We construct a continent‐wide model of Africa's Crust Evaluated with ADAMA's Rayleigh Phase maps (ACE‐ADAMA‐RP). Dispersion maps, and uncertainties, are obtained with a probabilistic approach. This model update, and a crustal taxonomy derived from unsupervised machine learning, reveals that the architecture of Africa's crust can be classified into two main types: primitive (C1: faster velocities with little gradients) and modified (C2–C4: slower velocities in the shallow crust with more pronounced gradients). The Archean shields are “primitive,” showing little variation or secular evolution. The basins, orogens, and continental margins are “modified” and retain imprints of surface deformation. The crustal taxonomy is obtained without a‐priori geological information and differs from previous classification schemes. While most of our reported features are robust, probabilistic modeling suggests caution in the quantitative interpretations where illumination is compromised by low‐quality measurements, sparse coverage or both. Future extension of our approach to other complementary seismological and geophysical data sets—for example, multimode earthquake dispersion, receiver functions, gravity, and mineral physics, will enable continent‐wide lithospheric modeling that extends resolution to the upper mantle.

Investigation of the quality of functioning of a quadrature signals receiver of a continuous automatic cab signaling on records of real signals and interference

Ensuring a safe transportation process is impossible without the use of complex systems for interval regulation of train traffic. In particular, such a system is a continuous automatic cab signaling (CACS) where signals are affected by various kinds of interference, leading to system failures. To reduce the number of system failures, the authors previously proposed a quadrature signal receiver of CACS (QSR of CACS). The article presents the methodology and results of an experimental study of the quality of the functioning of the QSR of CACS on records of real signals and interference. A simulation model of the test facility has been developed in the Simulink simulation environment. The study showed that the QSR of CACS provides reliable reception of CACS code signals with a low level of interference. However, it needs to apply additional measures in it to reduce the influence of impulse noise.

Jrfapp: A Python Package for Joint Inversion of Apparent S-Wave Velocity and Receiver Function Time Series

jrfapp: A Python Package for Joint Inversion of Apparent S-Wave Velocity and Receiver Function Time Series

Geological Implications of Crustal Thickness, Vp/Vs, and Poisson's Ratio Beneath the Western Mesopotamian Plain Edge, Iraq

To know the crustal structure, composition, deformation, and evolution, it is crucial to appreciate the ratio between the velocity of the primary wave and the velocity of the secondary wave, Poisson's ratio, and crustal thickness. The receiver function inversion method was used to determine the (h), the primary wave and velocity of the secondary wave, and the Poisson's ratio. The chosen earthquakes' receiver functions were examined after being recorded at the two seismic stations in Anbar ANB1 and Karbala KAR2. According to the research, the westernmost part of the Mesopotamian plain has a crustal thickness between 44 and 46 km. The calculated crustal thickness value agrees with data from platform regions. The study area's crustal primary wave and velocity of the secondary wave is 1.79. For the platform region and andesitic rocks, this value is typical. The crustal Poisson's ratio was reported to have an intermediate value (0.27). It demonstrates that there are not many differences between the felsic and mafic compositions of the crust.

Sedimentary structures of the western part of the Indo-Gangetic Plain and Siwalik Himalaya inferred from receiver function inversion

SUMMARY We investigated sedimentary thickness and shear wave velocity structure in the western part of the Indo-Gangetic Plain (Punjab and Haryana Plain) and adjoining Siwalik Himalaya with the help of receiver function inversion at 20 broad-band seismological stations. This region is one of the most seismically vulnerable zones of the world due to the presence of thick surface sediments in the foreland basin that can amplify seismic waves and cause huge damage due to the earthquakes of the Himalaya. The study reveals a progressive thickening of sediments from southwest to northeast. The basement depth varies from ∼1.5 to 1.7 km in the Central Alluvium Plain, ∼1.8 to 2.8 km in the Zone of Terminal Fans, and attains a maximum of ∼3.8 km near the Himalayan Frontal Thrust. The inverted models show the presence of soft alluvial with extremely low Vs (< 0.5 km s−1) and high Vp/Vs (∼2.5–3.0) at the top ∼400–700 m of the surface at most of the stations. A comparatively higher velocity of surface sediments observed at northern stations suggests the presence of compact sediments at the surface. The layered sedimentary structure revealed by the S-wave velocity models supports the previous geophysical investigations using borehole data. The velocity–depth structure obtained in this study is important for evaluating the seismic hazard of the densely populated urban areas spread over this region.

Detailed crustal deformation in the Ailaoshan Orogenic Belt revealed by receiver functions from a dense array

The Ailaoshan Orogenic Belt (AOB), located at the southeastern boundary of the Tibetan Plateau, is an ideal place for investigating the mechanisms of lateral growth of Tibet. Using the data recorded by a dense seismic array across the Ailaoshan belt, we investigate the detailed lateral variations of crustal anisotropy on the basis of Pms phase of receiver functions. Remarkable crustal anisotropy is observed throughout this study region with a mean delay time of 0.33 ± 0.19 s, indicating the anisotropy primarily originates in the middle-lower crust. The fast directions beneath the AOB including the Ailaoshan-Red River shear zone (ARRSZ) and its western low-grade metamorphic unit generally align with the NW-SE strike of ARRSZ. The weak anisotropy in the South China Block (SCB) argues that the block is relatively stable, with limited internal deformation. Meanwhile, the anisotropy beneath the western boundary of the SCB is strong, and the N-S oriented fast direction is influenced by both the crustal stress and Xiaojiang Fault. Combining the high Vp/Vs and significant lateral variations of crustal anisotropy parameters, we suggest that the strike-slip motion along the ARRSZ induces the partial melting and pronounced anisotropy in the middle-lower crust of AOB, without the presence of crustal flow. The differences between crustal and mantle anisotropy indicate crust-mantle decoupling deformation of the AOB, supporting the block extrusion model occurring only in the crustal scale as the primary deformation pattern.

Estimating seismometer component orientation of the Brazilian seismographic network using teleseismic P-wave particle motion analysis and directional statistics

Seismometer component orientation is one of the most critical parameters for seismology, especially for methods that take wave field polarization effects into account, such as receiver functions, seismic anisotropy, surface wave dispersion and seismic tomography. Using P-wave particle motion combined with directional statistics, we estimated the orientation for 90 stations of the Brazilian Seismographic Network (RSBR), corresponding to 105 seismometers. The difference between the number of stations (90) and seismometers (105) is due to instrument replacement during station operation. We observed that 81 seismometers have an orientation of ±5° and 19 seismometers greater than ±5°. Only one seismometer has an orientation greater than 90° (BOAV station). ON and BR are the subnetworks of RSBR with the largest number of seismometers with orientation errors bigger than ±10° (6 for each one) and NB is the only subnetwork with no misoriented station. All the estimated orientations were close to the measured in the field, which is evidence that our approach is more accurate to detect misoriented stations than previous studies that analyzed the orientation of RSBR stations. The Watson's U2 statistical test indicated that the back azimuth deviations for Precambrian Basement and Phanerozoic Cover were significantly different, suggesting that the regional geology could affect back azimuth estimates.

New constraints on the shear wave velocity structure of the Ivrea geophysical body from seismic ambient noise tomography (Ivrea-Verbano Zone, Alps)

SUMMARY We performed seismic ambient noise tomography to investigate the shallow crustal structure around the Ivrea geophysical body (IGB) in the Ivrea-Verbano Zone (IVZ). We achieved higher resolution with respect to previous tomographic works covering the Western Alps, by processing seismic data collected by both permanent and temporary seismic networks (61 broad-band seismic stations in total). This included IvreaArray, a temporary, passive seismic experiment designed to investigate the IVZ crustal structure. Starting from continuous seismic ambient noise recordings, we measured and inverted the dispersion of the group velocity of surface Rayleigh waves (fundamental mode) in the period range 4–25 s. We obtained a new, 3-D vS model of the IVZ crust via the stochastic neighbourhood algorithm (NA), with the highest resolution between 3 to 40 km depth. The fast and shallow shear wave velocity anomaly associated with the IGB presents velocities of 3.6 km s−1 directly at the surface, in remarkable agreement with the location of the exposed lower-to-middle crustal and mantle outcrops. This suggests a continuity between the surface geological observations and the subsurface geophysical anomalies. The fast IGB structure reaches vS of 4 km s−1 at 20–25 km depth, at the boundary between the European and Adriatic tectonic plates, and in correspondence with the earlier identified Moho jump in the same area. The interpretation of a very shallow reaching IGB is further supported by the comparison of our new results with recent geophysical investigations, based on receiver functions and gravity anomaly data. By combining the new geophysical constraints and the geological observations at the surface, we provide a new structural interpretation of the IGB, which features lower crustal and mantle rocks at upper crustal depths. The comparison of the obtained vS values with the physical properties from laboratory analysis of local rock samples suggests that the bulk of the IGB consists of a combination of mantle peridotite, ultramafic and lower crustal rocks, bound in a heterogeneous structure. These new findings, based on vS tomography, corroborate the recent interpretation for which the Balmuccia peridotite outcrops are continuously linked to the IGB structure beneath. The new outcomes contribute to a multidisciplinary framework for the interpretation of the forthcoming results of the scientific drilling project DIVE. DIVE aims at probing the lower continental crust and its transition to the mantle, with two ongoing and one future boreholes (down to 4 km depth) in the IVZ area, providing new, complementary information on rock structure and composition across scales. In this framework, we constrain the upper crustal IGB geometries and lithology based on new evidence for vS, connecting prior crustal knowledge to recent active seismic investigations.

Crustal properties investigation in the Meratus Complex Zone, Southeast Borneo, using H-κ stacking of receiver function analysis

We study crustal properties using teleseismic receiver functions extracted from two broadband seismic stations around the Meratus Complex region, Southeast Borneo. The two seismic stations, BBKI and KBKI, are distributed in the Meratus Mountains and Pulau Laut. We then employ the H-κ stacking technique to measure the crustal thickness (H) and the crustal average Vp/Vs ratio (κ). The method works quite well for receiver functions recorded on seismic station KBKI. The thickness of the crust and Vp/Vs ratio underneath this station are around 38 km and 1.64, respectively. On the contrary, the standard H-κ method fails in estimating the crustal properties beneath station BBKI due to significant sedimentary contamination on the receiver functions. Thus, we apply the sequential H-κ stacking analysis by removing the sedimentary effect and increasing the Moho Ps amplitude. The result shows the thickness of the crust is around 29 km, although the Vp/Vs is not well determined. The research findings from both stations are comparable to those derived from the previous receiver function study using inversion modelling. The crustal properties around this area may reflect the continental material part of the Sunda continental block.

Insights From Layered Anisotropy Beneath Southern New England: From Ancient Tectonism to Present‐Day Mantle Flow

AbstractSeismic anisotropy beneath eastern North America, as expressed in shear wave splitting observations, has been attributed to plate motion‐parallel shear in the asthenosphere, resulting in fast axes aligned with the plate motion. However, deviations of fast axes from plate motion directions are observed near major tectonic boundaries of the Appalachians, indicating contributions from lithospheric anisotropy associated with past tectonic processes. In this study, we conduct anisotropic receiver function (RF) analysis using data from a dense seismic array traversing the New England Appalachians in Connecticut to examine anisotropic layers in the crust and upper mantle and correlate them with past tectonic processes as well as present‐day mantle flow. We use the harmonic decomposition method to separate directionally‐dependent variations of RFs and focus on features with the same harmonic signals observed across multiple stations. Within the crust, there are multiple features that may be correlated with stratification in the Hartford Basin, faults in the Taconic thrust belt, shear zones formed during Salinic/Acadian terrane accretion events, and orogen‐parallel crustal flow in the Acadian orogenic plateau. We apply a Bayesian inversion method to obtain quantitative constraints on the direction and strength of intra‐crustal anisotropy beneath the Hartford Basin. In the upper mantle, we identify a fossil shear zone possibly formed during oblique subduction of Rheic Ocean lithosphere. We also find evidence for a plate motion‐parallel flow zone in the asthenosphere that is likely disturbed by mantle upwelling near the southern margin of the Northern Appalachian Anomaly in the eastern part of the study area.