Receiver Functions Research Articles

The Dynamic Crust of Northern Afar and Adjacent Rift Margins: New Evidence From Receiver Function Analysis in Eritrea and Ethiopia

AbstractAfar is undergoing the final stages of continental rifting and hosts the triple junction between the Red Sea, Gulf of Aden, and Main Ethiopian rifts. To better understand the nature of the crust and continental breakup in the region, we calculate teleseismic receiver functions across northeastern Afar and the Danakil microplate, using new data from a regional deployment in Eritrea. We estimate the Moho depth and bulk crustal VP/VS ratio using the H‐κ stacking method. The heterogeneity of our crustal thickness estimates (∼19–35 km) indicates that the Danakil microplate has undergone stretching and crustal thinning. By investigating the relationship between crustal thickness and topographic elevation, we estimate the regional crustal bulk density as ρc ≈ 2,850 ± 20 kg m−3, which is higher than expected, given the crustal thickness of the region. We show that topography is 1.5 ± 0.4 km higher than would be expected due to crustal isostasy alone. We propose that this topography is supported by the same hot mantle upwelling suggested to be responsible for the onset of rifting in East Africa. Uplift is generated due to the presence of a hot thermal anomaly beneath the plate and by thinning of the lithospheric mantle. Our results are consistent with a number of independent constraints on the thermal structure of the asthenospheric and lithospheric mantle. Evidence of melt within the crust is provided by anomalously high VP/VS ratios of >1.9, demonstrating that magma‐assisted extension continues to be important in the final stages of continental breakup.

Depth‐Dependent Anisotropy Along Northwest Segment of the North Anatolian Fault Zone: Evidence for Paleo‐Tectonic Structures Contributing to Overall Complexity

AbstractThe North Anatolian Fault Zone (NAFZ) is a prominent tectonic structure with a significant impact on the observed active deformation in Türkiye. Detailed knowledge of the seismic anisotropy in the crust and mantle along this nascent shear deformation zone provides insights into the kinematics associated with past and present tectonic events. We employed teleseismic earthquakes observed by the Dense Array North Anatolia seismic network to map 3‐ D variations in crustal and mantle anisotropy in/around the NW segment of the NAFZ. To achieve this, we first performed a harmonic decomposition analysis of P‐receiver functions. The results were then used as a priori information to conduct an anisotropic receiver function inversion with the Neighborhood Algorithm that enabled imaging of the actual orientation and geometry of anisotropic structures. SKS splitting measurements are further used to make a comparison between the anisotropic behavior of crustal and mantle structures. Crustal anisotropy parameters estimated in our analyses/models well identify the signature of deformation caused by accumulated strain in the earthquake cycle through the strike of shallow cutting faults in the brittle crust beneath the NAFZ. Diffuse intense anisotropic energy at lower crustal depths was attributed to lattice preferred orientation of crystals or partially molten lenses elongated along the shear direction. Strong harmonic energy variations beneath the northern part of the Istanbul Zone likely reflect imprints of LPO‐originated frozen fabric at shallow depths (0–20 km) associated with the palaeotectonic Odessa Shelf, Intra‐Pontide Suture Zones or remnants of the Tethys Ocean.

A 3D computational model for ground motion simulation in Peninsular India

Due to the gradual and constant accumulation of seismic energy, Peninsular India (PI) is typically considered seismically stable with low to moderate seismicity. The seismic studies in Peninsular India always resorted to synthetic ground motion simulations, because of the limited instrumentation and hence lack of recorded data. In the absence of a well-defined medium model for PI, the usual practice is to use simple site proxies or one-dimensional velocity structures for ground motion simulations. However, the region consists of multi-scale geometric complexities, significant topography, and sedimentary basins and is surrounded by deep oceans. Thus, the radiated seismic wave field in the region is influenced by the medium properties and in the absence of a well-defined tomography model the reliable estimation of seismic hazard is a challenging problem in PI. Therefore, the seismic wave propagation in PI can be investigated using numerical simulation with reliable 3D computational model for PI, incorporating the knowledge of the underlying Earth structure. Hence, the present study attempts to develop a sophisticated three-dimensional (3D) medium model of Peninsular India for physics-based ground motion simulations for regional earthquakes. This is aided by the availability of one-dimensional (1D) velocity models and the crustal structure from the receiver function analysis which provides valuable insight into the variation of material properties in the region. In the present study, >100 s of 1D velocity profiles are collected from various literature, which is then grouped under 23 different geological regions identified in PI (as per GSI (2000)). The averaged material properties are assigned per each geological region and the information on sediment depths, basin geometry, topography, and bathymetry are incorporated. We use the spectral element method (SEM) to calibrate our 3D computational model by simulating synthetic seismograms and comparing them to recorded ground motions for two past earthquakes: the 2001 Mw 7.6 Bhuj earthquake and the 1997 Mw 5.8 Jabalpur earthquake. Further, the seismic waveforms at the near field of 2001 Mw 7.6 Bhuj event are simulated using a refined regional model. The spatial variability of associated seismic intensities and peak ground velocity (PGV) amplification are investigated. In addition, a study of the impact of model depth truncation and sphericity on ground motion is also conducted. The implemented medium model is the first of its kind for Peninsular India and can reliably be used in seismic wave propagation studies in the region. The simulated outcomes from the model are of engineering importance as these results can be used for seismic hazard assessment of the region.

Understanding the structure of crust and shallow upper mantle beneath western Tibet through the joint inversion of Rayleigh wave group velocity dispersion with interpolated receiver functions

We present the 3-D shear wave velocity variationin the crust and upper mantle structure of western Tibet using a joint inversion of Rayleigh wave group velocity dispersion and interpolated receiver function. The lateral sampling of the receiver function and the surface wave dispersion is equalized by the interpolation scheme. This new method of spatial interpolation eliminates difficulties caused by back azimuthal variation in the receiver function. We have jointly inverted these interpolated receiver functions with the Rayleigh wave group velocity dispersions obtained from the earthquake and ambient noise tomography to map the tectonic structure of western Tibet. The study reveals remarkable variation of the seismic characteristics in the crust and shallow upper mantle of western Tibet. The mapped Moho depth is ~74 km in the northern part of the plateau, whereas it is ~69 km in the southern part. A thick mid-crustal slow-velocity zone from ~10km up to ~40 km depth has been observed, which could be attributed to partial melt. A noteworthy finding is the variation in slower crustal velocity between the south and north of the Lhasa Block (LB). We interpret this as the Indian crust that has been underthrusted beneath the Tibetan plateau, and the northern limit of the Indian lithosphere extends to the southern part of the LB.

Deep Structure of the Baikal Rift Zone and Central Mongolia

The upper mantle and the transition zone of the Baikal rift zone (BRZ) are studied. The observations are analyzed using P-wave receiver functions. It is found that in the BRZ central and northeastern part, the P410s converted seismic phase is preceded by a precursory wave with negative polarity which is formed in the low S-wave velocity layer at a depth of 350–410 km. A similar precursory wave with low S-wave velocity and negative polarity is formed at a depth of 600–660 km. The low-velocity layers are interpreted as resulting from the hydration of wadsleyite and ringwoodite during the subduction of the Pacific lithosphere. A similar study of the mantle in Central Mongolia found no expected signs of hydration. Modeling of the lithosphere–asthenosphere system in Central Mongolia by joint inversion of the body wave receiver functions and surface wave dispersion curves reveals a very thin lithospheric lid beneath Khangai and a thick layered asthenosphere to a depth of 200 km with a lithospheric inclusion between low-velocity layers.

Crustal thickness variation of Kalimantan and Sulawesi Region from Teleseismic receiver function

Crustal thickness variation of Kalimantan and Sulawesi Region from Teleseismic receiver function

A high-resolution seismic velocity model for East Asia using full-waveform tomography: Constraints on India-Asia collisional tectonics

East Asia's complex tectonic history makes it a unique place for plate-tectonic studies. While the crust and mantle structures beneath this region have been extensively investigated by receiver function and many different kinds of tomography studies, discrepancies among seismic velocity models have complicated our understanding of East Asia's tectonic evolution. Here, we use full-waveform adjoint tomography to study the seismic P- and S-wave velocity structures of the crust and mantle beneath East Asia. We use a large waveform dataset from 141 earthquakes recorded by 4640 broadband seismic stations. Our inversion optimizes the normalized correlation coefficient between synthetic and observed seismograms within individual time windows of different seismic phases, from P arrivals to the slowest surface waves. This approach allows us to fit regional multipathed waveforms and provides very well-resolved seismic P- and S-wave velocity structures from the crust to depths of about 800 km. Furthermore, we are able to obtain this resolution over most of East Asia allowing us to compare the amplitude and shapes of anomalies across the entire region. The observed structures in our new full-waveform tomography model (FWEA23) correlate well with tectonic units, revealing sharp contrasts across tectonic boundaries. Our model reveals narrow and strong linear fast anomalies associated with subduction zones in the western Pacific and Southeast Asia. Compared to these oceanic slabs, we observe different and more complicated mantle structures beneath the India-Asia continental collision zone. While our model spans a wide region in East Asia, we focus our interpretations on the seismic velocity anomalies beneath Tibet and surrounding regions, as our model provides significant implications for India-Asia continental collision tectonics. Our analysis suggests that cratonic Indian lithosphere collided with Eurasian plate around 25–20 Ma, subsequently underthrusting beneath Tibet horizontally with little deformation. This resulted in thickening, destabilization, and delamination of the rheologically weaker and less buoyant mantle lithosphere beneath north-central Tibet. Moreover, our model reveals isolated lithospheric fragments within the mantle transition zone and lower mantle beneath southern Tibet and northern India, suggesting that non-cratonic Greater India continental lithosphere has broken off and sunk into the deeper mantle. The lack of continuity among these lithospheric fragments suggests a convective destabilization process, distinct from typical oceanic plate subduction. Our results emphasize the important role of rheological properties of continental lithospheres in shaping the dynamics of continental collision.

Examination of shallow and deep S-wave velocity structures from microtremor array measurements and receiver function analysis at strong-motion stations in Kathmandu basin, Nepal

The Himalayan collision zone, where the Indian Plate subducts beneath the Eurasian Plate at a low angle, has caused many devastating earthquakes. The Kathmandu basin, situated in this region, is surrounded by mountains on all sides and is filled with distinct soft lake sediments with a highly undulating bedrock topography. The basin has been experiencing rapid urbanization, and the growing population in its major cities has increased the vulnerability to seismic risk during future earthquakes. Several strong-motion stations have recently been deployed in the Kathmandu basin. It is expected that the data captured by this strong-motion station array will further enhance our understanding of site amplification in sedimentary basins. Clear P-to-S converted waves have been observed in the strong-motion records. In this study, we investigate the medium boundary that generated these converted waves. First, we estimate the shallow velocity structures, which correspond to the topographic slopes or surface geology, beneath the strong-motion stations. We then apply a receiver function analysis to the strong-motion records. The receiver function indicates that the interface between the soft sediment and seismic bedrock serves as a boundary that generates converted waves. The obtained results can be used for tuning three-dimensional velocity structures.Graphical

Evidence of a Low‐Velocity Zone in the Upper Mantle Beneath Cumbre Vieja Volcano (Canary Islands) Through Receiver Functions Analysis

AbstractOceanic volcanic islands have a complex internal structure due to their geological evolution. La Palma is the second youngest in the Canary Islands. Historical eruptions occurred in the Cumbre Vieja (CV) volcanic system, including the last 2021 Tajogaite eruption. The seismicity before, during, and after the Tajogaite eruption provides valuable information for a better understanding of deep volcanic processes. We applied the receiver function technique to image the crustal and upper mantle structures beneath the CV volcanic complex. We investigated the first 50 km of the lithosphere, identifying various seismic discontinuities. After a thin layer of recent volcanic rocks, we found a layer corresponding to the old oceanic crust. Beneath the Moho, we found a significant decrease in the seismic velocities until about 37 km of depth. This low‐velocity layer corresponds to a zone of partial melting with interconnected magmatic chambers that feed the volcanic activity in CV.

Ancient Craton‐Wide Mid‐Lithosphere Discontinuity Controlled by Pargasite Channels

AbstractThe mechanisms governing a commonly observed seismic velocity drop in the cratonic lithosphere, referred to as the mid‐lithospheric discontinuity (MLD), have been widely debated. To identify the composition and seismic structure of MLDs, we have analyzed Sp receiver functions (SRF) and mantle xenocrysts for six regions across Australia. We utilize locations where seismic stations and kimberlite‐hosted mantle xenocrysts are both available, allowing for comparison between seismological and petrological constraints. Our results show negative SRF phases indicative of the MLD coincide with clinopyroxene‐depleted zones at 60–140 km depth. Clinopyroxenes with different chemical compositions across the MLD define a litho‐chemical discontinuity. Modeling and experimental data show that MLDs may be explained by modified lherzolite with 10%–20% modal pargasite. Pargasite MLDs may form when rising H2O‐bearing melts cross the amphibole dehydration curve and react with clinopyroxene in lherzolite. Because the amphibole dehydration curve is isobaric at 80–120 km, pargasite will be precipitated as horizontal channels.

Joint Inversion of Surface-Wave Dispersions and Receiver Functions Based on Deep Learning

Abstract We proposed a deep learning (DL) method to derive VS models from joint inversion of Rayleigh-wave dispersions and receiver functions, which is based on multilabel convolutional neural network and recurrent neural network. We used a spline-based approach to generate synthetic models instead of directly using existing models to build the training data set, which improves the generalization of the method. Unlike the traditional methods, which usually set a fixed VP/VS ratio, our method makes full use of the powerful data mining ability of DL to invert the VS models assuming different VP/VS ratios. A loss function is specially designed that focuses on key features of the model space, for example, the shape and depth of Moho. Synthetic tests demonstrate that the proposed method is accurate and fast. Application to the southeast margin of the Tibetan Plateau shows results consistent with the previous joint inversion with P constraints, indicating the proposed method is reliable and robust.

Waveform Fitting of Receiver Functions for Enhanced Retrieval of Crustal Structure in the Presence of Sediments

AbstractThe receiver function technique is widely used to image crustal structure using P‐to‐S converted phases at the Moho discontinuity. However, the presence of sedimentary layer generates additional P‐to‐S conversions and reverberations, which can overprint the Moho phases and pose problems in imaging crustal structure with standard receiver function techniques. We introduce a robust two‐step method that uses H‐κ stacking to determine average thickness and Vp/Vs of the sedimentary layer, followed by waveform‐fitting of the observed receiver function to constrain the average crustal thickness and sub‐sediment Vp/Vs. We tested the method using both synthetic data and real‐data from stations located on sedimentary layers in the Netherlands and USA. We show that the new method outperforms other common approaches in retrieving accurate Moho depth and sub‐sediment Vp/Vs estimates, even in cases where the Moho phases are completely overprinted by large‐amplitude phases related to sedimentary layers.

A Taxonomy of Upper‐Mantle Stratification in the US

AbstractThe investigation of upper mantle structure beneath the US has revealed a growing diversity of discontinuities within, across, and underneath the sub‐continental lithosphere. As the complexity and variability of these detected discontinuities increase—for example, velocity increase/decrease, number of layers and depth—it is hard to judge which constraints are robust and which explanatory models generalize to the largest set of constraints. Much work has been done to image discontinuities of interest using S‐waves that convert to P‐waves (or top‐side reflected SS waves). A higher resolution method using P‐to‐S scattered waves is preferred but often obscured by multiply reflected waves trapped in a shallower layer, limiting the visibility of deeper boundaries. Here, we address the interference problem and re‐evaluate upper mantle stratification using filtered P‐to‐S receiver functions (Ps‐RFs) interpreted using unsupervised machine‐learning. Robust insight into upper mantle layering is facilitated with CRISP‐RF: Clean Receiver‐Function Imaging using Sparse Radon Filters. Subsequent sequencing and clustering organizes the polarity‐filtered Ps‐RFs into distinct depth‐based clusters. We find three types of upper mantle stratification beneath the old and stable continental US: (a) intra‐lithosphere discontinuities (paired or single boundary), (b) transitional discontinuities (single boundary or with a top layer), and (c) sub‐lithosphere discontinuities. Our findings contribute a more nuanced understanding of mantle discontinuities, offering new perspectives on the nature of upper mantle layering beneath continents.

Unveiling the Distinct Structure of the Upper Mantle Beneath the Canary and Madeira Hotspots, as Depicted by the 660, 410, and X Discontinuities

AbstractThe Canary and Madeira Islands are two distinct hotspots in the Central‐East Atlantic that are close to each other. Their volcanism is generally attributed to underlying mantle plumes, but the detailed structure of these plumes is still not well understood. The thermal and compositional structure of the plume introduces complexities in the phase transitions of the mantle, which impact the depth and magnitude of seismic discontinuities. We use 1,268 high‐quality receiver functions from stations located at the two hotspots to detect P‐to‐s converted phases through a common‐conversion point stacking approach and conduct a detailed analysis of mantle seismic discontinuities. The results show that both hotspots are characterized by a thin mantle transition zone (MTZ), with sharp 410 and 660 discontinuities at depths of 428–421 km and 647–664 km, beneath the Canaries and Madeira respectively. The results indicate that the Canary plume crosses the MTZ, whereas the Madeira plume mainly influences the upper portion of the MTZ. Furthermore, we find reliable detections of a sharp X discontinuity beneath the Canaries at 287 km. Its presence suggests the accumulation of silica‐rich recycled eclogite at these depths. We also use the amplitudes of P410s and PXs to derive velocity jumps at corresponding discontinuities. Based on these measurements, we estimate that the basalt proportion is 60%–80%, with accumulation being more significant in the Canaries than in Madeira. The MTZ thickness, the presence of the X discontinuity, and the high basalt proportion provide compelling evidence for a deep‐rooted thermochemical plume beneath the study area.

Optimized wireless communication for 6G signal processing using memory‐augmented deep unfolding network

SummaryIn this manuscript, optimized wireless communication for 6G signals processing using memory‐augmented deep unfolding network (WC‐6G‐SP‐MADUN) is proposed. The aim of this work is to give a memory‐augmented deep unfolding network approach with a focus on belief propagation decoding of error correction codes and multiple‐input multiple‐output (MIMO) wireless systems. In contrast, compares and tabulates the performance of memory‐augmented deep unfolding methodologies with classic principled procedures, which could potentially enable 6G. Here, the wireless communication for 6G signal processing methods was used for a variety of receiver functions, including self‐interference cancelation, signal estimation and detection using memory‐augmented deep unfolding Network (MADUN), and advanced error correction by clouded leopard optimization algorithm (CLOA). Furthermore, outline the general framework of memory‐augmented deep unfolding method for any kind of signal processing application and talk about cutting edge research paths that will make it possible to build communication networks in the future. The proposed approach is implemented in MATLAB; the performance is analyzed with performance metrics such as reliability, complexity, SNR gain, convergence rate, hardware efficiency, and quantization bit width. The proposed approach attains higher SNR gain 17.96%, 24.75%, and 30.09%; lower complexity 17.92%, 23.41%, and 30.13%; and higher convergence 18.01%, 25.41%, and 31.39% when analyzed with existing techniques such as redefined wireless communication for 6G: signal processing meets deep learning and deep unfolding (RWC‐6G‐SP‐DU), lower complication deep unfolded neural network receiver for MIMO under probability data association detector (LC‐DUNN‐MIMO‐PDAD), and deep learning in physical layer communications: evolution with prospects on 5G and 6G networks (DL‐PLC‐6G), respectively.

Local earthquake tomography of the Aegean crust: Implications for active deformation, large earthquakes, and arc volcanism

Three-dimensional velocity models can enable accurate earthquake location, improved seismic hazard assessment, and can enhance our understanding of geodynamic processes. This is particularly true for areas such as the Aegean, where the crust is marked by active volcanoes as well as shear and rift zones, all of which make this region highly heterogeneous. This work describes the application of local earthquake tomography in order to elucidate the crustal structure of the Aegean. We utilized a dataset of 2,135,625 P arrivals and 1,095,515 S arrivals from 178,676 events recorded across 324 stations. Results indicate that the boundary separating the slow Vp from the fast anomalies at 20–40 km depth is well correlated with the Aegean Moho obtained from receiver functions. A slow Vp (−5%) zone at depths >30 km along the Hellenic arc is likely associated with basal underplating that forms the base of the crust. On the other hand, fast Vp and moderate to high Vp/Vs (1.77–1.92) in the upper-middle crust may indicate evolving metamorphic core complexes. Steeply dipping clusters of earthquakes could highlight pathways of fluid migration from the base to the middle/upper crust, albeit more detailed seismological studies are needed to confirm this interpretation. The 3D velocity model can be also utilized in order to investigate the amount of melt fraction beneath active volcanoes and the influence of fluids on the rupture zone of large earthquakes. Our results show a melt fraction between 4 and 10% beneath active volcanoes, with the largest volume of melt present beneath Santorini caldera. The ruptures zones of 10 large earthquakes (Mw ≥ 6.0), show slow Vp anomalies and moderate to high Vp/Vs (> 1.77), suggesting that these events were likely triggered by weakening of their source zones by fluid activity.

Three‐Dimensional Teleseismic Elastic Reverse‐Time Migration With Deconvolution Imaging Condition and Its Application to Southwest Japan

AbstractWe have developed a novel deconvolution‐based reverse time migration method to image lithospheric structures using teleseismic data recorded by dense seismic arrays. Unlike traditional approaches that rely on the retrieval of Green's functions (or receiver functions), the new method directly utilizes the recorded three‐component (3‐C) seismic waveforms to reconstruct subsurface wavefields, which has the advantage of avoiding the estimation and removal of source time function. Importantly, it enables the asymptotic estimation of P‐to‐S transmission conversion coefficients at the elastic discontinuities and enhances resolving power for the fine‐scale heterogeneities. Taking these improvements, we have obtained a full three‐dimensional (3‐D) high‐resolution image of the subduction zone beneath southwest Japan. This image provides a comprehensive configuration of the severely deformed subducting plate beneath the Kii Channel and Kinki Peninsula.

Significant variations in the structure and composition of the crust beneath northwestern India: Imprints of magmatism

Northwestern India was affected by the Neoproterozoic Malani and Erinpura magmatism, Cretaceous rifting, and magmatism associated with the Deccan precursors. Receiver function analysis was done using 1704 RFs from 18 stations to image and decipher the imprints of the magmatic events and comprehend the crustal modifications in terms of variations in structure and composition. The H-κ grid search and neighbourhood inversion techniques are used to retrieve the crustal structure and Vp/Vs ratios (κ). The study reveals significant variations in the crustal thickness (H = 34 km to 43 km), composition (κ=1.74 to 1.92) and shear wave velocity structure across northwestern India. The terrains encompassing the Erinpura granite and the Malani igneous suite in the vicinity of the Barmer rift are characterized by a thick (H ≥ 41 km) crust with felsic to intermediate composition (κ ≤ 1.81). The crust beneath the southern part of the Marwar basin is ≈ 38–40 km thick with mafic composition (κ> 1.81). The region around the Barmer rift has a thin crust (H = 34–36 km) with intermediate to mafic composition. The region hosting the Early Cretaceous to Paleogene alkaline complexes exhibits a high Vp/Vs ratio (κ= 1.91) that may be associated with mafic cumulates emplaced by magmatic events that overprint the signatures of the Malani event. The crust is heterogeneous with low/high velocity intracrustal layers that reflect the fractionation of magma at different depths. The mafic residue, together with magmatic intrusions, results in a mafic lowermost crust with high shear velocities of 3.9–4.0 km/s beneath most stations. Overall, northwestern India is characterized by a thick crust with intermediate crustal composition and intracrustal layering resulting from large scale magmatic events linked to the Neoproterozoic reorganization of plates and younger magmatic events.

Imaging the Sedimentary and Crustal Structure of the Luoyang Basin, Central China, Using a Dense Nodal Seismic Array

Abstract The Luoyang basin lies in the southern margin of the North China block, separating the trans-North China orogen to the north and the Qinling-Dabie orogen to the south. Determining how the basin formed is important for understanding the history of the North China block and its evolution during the Mesozoic and Cenozoic times. Based on the teleseismic data recorded by a dense nodal seismic array, we used the receiver function method to image the sedimentary and crustal structures in the Luoyang basin. Common conversion point stacking images show that the Moho is at a depth of ∼35 km on the south and west sides and slightly uplifted to ∼30 km below the northeastern basin. Two sets of P-to-S conversions are imaged in the shallow crust, separating the near-surface sediments into consolidated, semiconsolidated, and unconsolidated layers. The top of the consolidated sedimentary layer is close to the surface at the southern basin and present at a depth of ∼2 km beneath the central basin, then deepens to a depth of ∼3 km below the northern basin. The discontinuous interface in the sediments indicates that the sedimentary layer was truncated by some blind north-dipping normal faults. The northeastward thinning crust, thickening sedimentary layers, and dipping normal faults together indicate that the Luoyang basin evolved in association with the deep crustal extension response to the lithospheric thinning of the North China block. By superimposing the deep crustal extension, we propose that the present-day landform of the Luoyang basin was also shaped by fluvial erosion at the surface, which was accompanied by the expansion of Yihe and Luohe riverbeds during the Quaternary.

The crust-mantle velocity structure beneath the Wudalianchi-Erkeshan-Keluo volcanic belt by joint inversion of receiver functions and ambient noise

The Wudalianchi-Erkeshan-Keluo (WEK) volcanic belt is a significant component of intraplate volcanism in Northeast China and is composed of the Wudalianchi, Erkeshan, and Keluo volcanic clusters. Using joint inversion of receiver functions and ambient noise, we construct a high-resolution 3-D S-wave velocity model of the WEK volcanic belt and its adjacent region, taking advantage of a deployed dense seismic array around this volcanic belt. There is a prominent low-velocity anomaly at 8–15 km depth beneath the Wudalianchi volcanic cluster, suggesting the presence of a crustal magma chamber. Low-velocity anomalies are also observed at 30–35 km depth beneath the Erkeshan volcanic cluster and 30–40 km depth beneath the Keluo volcanic cluster, resulting in discontinuous velocity structures at the Moho discontinuity. We further find a distinct low-velocity anomaly in the uppermost mantle beneath the WEK volcanic belt. Combined with previous geophysical and geochemistry studies, we propose a magma system scenario for the WEK volcanic belt. The upwelling molten material from the asthenosphere accumulated in the uppermost mantle and the magma chamber was formed, which provided the same uppermost mantle magma sources for the WEK volcanic belt.