Wave Receiver Functions Research Articles

Translithospheric magma plumbing system fossilized in the Emeishan large igneous province

Lighting up the magma plumbing system beneath a large igneous province (LIP) is challenging because the complex magma migration paths are often covered by flood basalts and sediments. Here, we present a three-dimensional seismic image of the Permian Emeishan LIP in Southwest China, constructed by joint inversion of Rayleigh wave dispersions and receiver functions. The results outline a cylindrical, high-velocity anomaly extending to ~135 km depth below the inner zone of this continental LIP. The geometry and magnitude of the high-velocity anomaly suggest that it represents culminated crystallized materials of primary magmas, thereby mirroring a magma plumbing system fossilized in the lithosphere. Furthermore, our geodynamic models illustrate that the nearly vertical plumbing system was controlled by slow plate motion during the magma emplacement. The plume head beneath a nearly static plate has higher thermal buoyancy and thus promotes more intensive magma emplacement. This phenomenon may apply to other LIPs throughout Earth’s history

Fine Crust‐Mantle Structure of the Major Tectonic Boundaries Between the North China Craton and Central Asian Orogenic Belt Revealed From Rayleigh Wave Phase Velocities and Receiver Functions

AbstractThe Tanlu Fault Zone (TLFZ) and Chifeng‐Kaiyuan Fault (CKF) serve as tectonic boundaries between the North China Craton (NCC) and the Central Asian Orogenic Belt (CAOB). Clarifying the refined structure of these tectonic boundaries is crucial for understanding the relationships between the tectonic units and the heterogeneity in the destruction of the NCC. In this study, two linear seismic arrays were deployed across these tectonic boundaries. Based on the phase velocity dispersion and receiver functions extracted from the seismic arrays, the Hamiltonian Monte Carlo algorithm was employed for the joint inversion of the S‐wave velocity (Vs) in the crust and uppermost mantle. The Vs model was then used to correct the time differences in common conversion point (CCP) stacking. The CCP stacking results indicate that the boundary faults TLFZ and CKF are both whole‐crustal faults that separate the NCC and CAOB. The Vs structure showed a significant low‐velocity anomaly in the mantle beneath the NCC, with intense seismic activity within the crust. This suggests that the NCC was affected by the subduction of the Western Pacific, leading to crustal and mantle destruction. In contrast, the CAOB exhibited a clear high‐velocity anomaly with relatively stable crustal structures. We believe that the NCC and CAOB have undergone structural modification and destruction due to the closure of the Paleo‐Asian Ocean and the activities of the TLFZ since the Late Mesozoic. During the Cenozoic, the region east of the TLFZ experienced more significant destruction in the NCC than the other adjacent tectonic units.

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.

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.

Bayesian inversion for modelling 3-D density structures in the eastern margin of Bayan Har block and its tectonic implications

SUMMARY Recently, the eastern margin of the Tibetan plateau has experienced three moderate earthquakes. Deep crustal structure is a critical factor for understanding the seismotectonic environment and deformations in this region. Although several multidisciplinary geophysical approaches have been used to develop crustal structure models for this area, substantial inconsistencies still persist in the results due to the intricate structural properties of crust. In this study, we introduce a novel approach to construct a 3-D crustal model with assimilated density structure (MADS) at a resolution of 0.2° × 0.2° × 5 km. This model is built by integrating data from various source, including multisource gravity anomalies and diverse available crustal structure reference models. We use a model-assimilated gravity inversion process with Bayesian parameter optimization. First, we combine the terrestrial gravity profiles data set with published global gravity field models to create a data set rich in reliable high-frequency anomaly information. Secondly, we incorporate three reference models derived from different seismological methods as prior constraints for our model. These models encompass seismic tomography, surface wave dispersion and receiver function data. We optimize the hyper-parameters of these constraints using the Bayesian criterion. The results demonstrate that the MADS not only captures significant changes in the crustal density but also discerns subtle variations in the upper and middle crust, thereby providing detailed insights into the morphologies of major faults. For instance, the central section of the Longmenshan fault is revealed as a high-angle deep thrust feature, while the frontal section of the Longmenshan fault appears as a low-angle mid-deep thrust feature, and the Xianshuihe fault exhibits a vertical deep subduction feature. Additionally, our findings indicate a correlation between the locations of moderate-to-large earthquakes in this region and the high density-gradient zones or asperities with high density within MADS. We believe that the insights into density characteristics offered by the new MADS model can shed light on the study of asperities associated with recent moderate earthquakes and enhance our understanding of deformation in this region.

CRUSTAL THICKNESS BENEATH UTAKO ABUJA USING RECEIVER FUNCTION ANALYSIS

Understanding the Mohorovicc discontinuity (Moho depth) is crucial in unraveling the Earth’s subsurface structure. In this research, we present the report of the crustal thickness and the ratio of the primary wave to the shear wave velocities (p-to-s) beneath the seismic station located at the Nigerian Geological Survey Agency (NGSA) Abuja using Receiver Function Analysis (RFA). The primary objective of this study is to determine the depth of the Moho boundary and gain insights into the geological properties of the Earth’s crust and upper mantle beneath the receiver station. This analysis focused on teleseismic earthquakes that were recorded at the station between 2020 and 2021, with epicentere distances between 30o and 90o and magnitudes 5.5 or greater. The P wave Receiver Function (PRF) was modeled from iterative time domain deconvolution using the converted phases (Ps) at the Moho interface.

Differences in Scattering Properties of the Shallow Crusts of Earth, Mars, and the Moon Revealed by P‐Wave Receiver Functions

AbstractThe scattering properties of terrestrial planetary bodies can provide valuable insights into their shallow seismic structure, meteoritic impact history, and geological activity. Scattering properties of the shallow crusts of Earth, Mars, and the Moon are investigated by constructing P‐wave receiver functions (PRFs) from teleseismic waveforms with high signal‐to‐noise ratios. The authors’ analysis reveals that strong coda waves lead to significant variations in the PRF waveforms calculated using different time windows, and the stability of the PRF is primarily influenced by the fractional velocity fluctuation. Synthetic PRFs for various scattering media confirm these observations. Comparing the observed and synthetic PRFs, it is found that the fractional velocity fluctuation in the shallow crust is greater than ∼0.2 for the Moon but less than ∼0.2 for Earth and Mars. The authors further discuss possible mechanisms that could have affected the fractional velocity fluctuation and suggest that the distinct fractional velocity fluctuation between the Moon and Earth/Mars is mainly due to differences in the water content of the crustal rocks of the three planetary bodies.

Crustal structure beneath and surrounding the Appalachian Basin region of Eastern North America imaged by joint inversion of P wave receiver functions and surface wave dispersion measurements

SUMMARY To evaluate the plate flexure model for the formation of the Appalachian Basin, we investigate the extent to which crustal structure beneath and surrounding the basin was modified by the Palaeozoic orogenic events that created the basin. We jointly invert receiver functions and surface wave dispersion measurements to obtain 1-D crustal Vs profiles for 261 seismic stations located within and around the basin. The average crustal thickness for the region is 44 km, and the crust gradually thins to the east, consistent with previous studies. Four areas of anomalous crust are identified with respect to the eastward thinning of the crust. An area of thick crust is found along the Grenville Front on the western side of the Appalachian Basin where the crust thickens by ∼5–10 km. Moho depths of up to 54 km in this region likely result from suture-thickening. The crust is thinner beneath the Neoproterozoic Scranton rift by ∼5–7 km, coincident with a ∼40 mGal Bouguer gravity high. Across the Neoproterozoic Rome Trough, the crust thins by ∼4–5 km, coincident with a ∼10 mGal Bouguer gravity high. Density models for these rifts show that the rift-related crustal thinning is sufficient to explain the gravity anomalies. The Vs models obtained for stations in the rifts indicate little, if any, mafic layering in the mid-crust and only a modest amount of mafic layering in the lower crust. In the northwestern portion of the Appalachian Basin in northeastern Ohio and northwestern Pennsylvania within the Elzevir block, another area of anomalously thick crust (50–52 km) is found. This region is not associated with any known tectonic structures or boundaries or a gravity anomaly. The lower ∼5–10 km of the crust in this region is characterized by high (>3.9 km s−1) shear wave velocities and thus appears to be mafic. The origin of anomalous crustal structure in all four areas is best attributed to Precambrian tectonic events that predate the formation of the Appalachian Basin, indicating that the crystalline crust beneath and surrounding the basin was not significantly affected by the Palaeozoic basin-forming orogenic events, a finding which supports the use of plate flexure models for understanding basin formation.

Lithosphere structure beneath the eastern Dharwar craton kimberlite field, India, inferred from joint inversion of surface wave dispersion and receiver function data

Cratons are the oldest parts of continents with cold and thick, high velocity, depleted lithosphere, which remained stable for billions of years. However, many of them are susceptible to deformation and erosion under various geological settings, such as near subduction zones and above mantle plumes. The role of mantle plumes in eroding the overlying lithosphere remains elusive, particularly for the fast-moving cratons. In this study, we examine the plume-lithosphere interaction in the Dharwar Craton located in South India by comparing the lithospheric property derived from the diamondiferous kimberlite xenoliths with its present-day shear velocity structure computed from collocated seismological data. The lithosphere-asthenosphere system beneath a 200 km long profile is explored at 13 locations through joint inversion of receiver function and surface wave phase dispersion data. The velocity model suggests 165 ± 15 km thick present-day lithosphere, marginally lower than 180 ± 20 km derived from the kimberlite xenoliths which erupted around 1100 Myr. The velocity increases from the Moho to ∼100 km depth, possibly due to the spinel-garnet peridotite transformation. The cold cratonic lithosphere of Eastern Dharwar craton (EDC) is manifested in very high shear velocities, up to >4.7 km/s indicating harzburgite-dunite composition, correlating well with the composition inferred from analysis of xenolith data. Our result suggests no significant thermal alteration to the lithosphere of the eastern Dharwar Craton despite being affected by three mantle plumes in the last 90 Myr.

The crust-mantle transition beneath Northeast China from P–wave receiver functions

Introduction: The continental lithosphere of the Northeast China Plain (NECP) is probed using P–wave receiver function analysis of passive seismic data. The NECP is much–discussed as it includes different geological provinces of varying tecto–magmatic origin and dates back from Archean to Holocene. Quantifying the tectonic and magmatic influences on the structure and composition of the lithosphere puts important constraints on evolution of NECP. For this, we explore 75 sites across NECP using receiver function analysis.Methods: A recently developed technique of inverting for 1–D S–wave velocity profiles beneath seismic stations that is based on the principles of Bayesian statistics (hierarchical transdimensional Bayesian Inversion; HTBI) is applied to receiver functions from the NECP. In addition, an improved crustal thickness–compressional to shear wave velocity ratio (H–κ) analysis was conducted to retrieve the crustal thickness and Vp/Vs ratio of the region. These estimated point measurements are integrated and systematically studied for a regional view of the current crustal architecture.Results and Discussion: We observe a laterally varying and highly complex lithosphere beneath the NECP. A shallower crust–mantle transition (≤32 km) characterises the Precambrian North China craton and Late Mesozoic–Cenozoic Songliao Basin from the adjacent Central Asian Orogenic Belt and the Changbaishan Volcanic field (35–40 km). Beneath the latter, low Vp/Vs ratios (~1.65) are obtained, whereas all other regions feature ratios in excess of 1.75. Multiple velocity gradients are observed at crustal depths within the craton, in contrast to the adjacent orogen, which indicates a higher degree of crustal complexity of the former. The width of the crust–mantle transition across the NECP is found to be mainly intermediate (2–7 km) and occasionally sharp (≤2 km). From our observations, we infer that there is a substantial difference between the eastern North China Craton’s lithospheric architecture and the rest of the NECP, with most of the NECP exhibiting more complexity than previously reported.

Ocean bottom seismometer (OBS) noise reduction from horizontal and vertical components using harmonic–percussive separation algorithms

Abstract. Records from ocean bottom seismometers (OBSs) are highly contaminated by noise, which is much stronger compared to data from most land stations, especially on the horizontal components. As a consequence, the high energy of the oceanic noise at frequencies below 1 Hz considerably complicates the analysis of the teleseismic earthquake signals recorded by OBSs. Previous studies suggested different approaches to remove low-frequency noises from OBS recordings but mainly focused on the vertical component. The records of horizontal components, which are crucial for the application of many methods in passive seismological analysis of body and surface waves, could not be much improved in the teleseismic frequency band. Here we introduce a noise reduction method, which is derived from the harmonic–percussive separation algorithms used in Zali et al. (2021), in order to separate long-lasting narrowband signals from broadband transients in the OBS signal. This leads to significant noise reduction of OBS records on both the vertical and horizontal components and increases the earthquake signal-to-noise ratio (SNR) without distortion of the broadband earthquake waveforms. This is demonstrated through tests with synthetic data. Both SNR and cross-correlation coefficients showed significant improvements for different realistic noise realizations. The application of denoised signals in surface wave analysis and receiver functions is discussed through tests with synthetic and real data.

Mars Seismology

For the first time, from early 2019 to the end of 2022, Mars’ shallow and deep interiors have been explored by seismology with the InSight mission. Thanks to the performances of its seismometers and the quality of their robotic installation on the ground, 1,319 seismic events have been detected, including about 90 marsquakes at teleseismic distances, with Mw from 2.5 to 4.7 and at least 6 impacts, the largest ones with craters larger than 130 m. A large fraction of these marsquakes occur in Cerberus Fossae, demonstrating active regional tectonics. Records of pressure-induced seismic noise and signals from the penetration of a heat flow probe have provided subsurface models below the lander. Deeper direct and secondary body wave phase travel time, receiver function, and surface wave analysis have provided the first interior models of Mars, including crustal thickness and crustal layering, mantle structure, thermal lithospheric thickness, and core radius and state. ▪ With InSight's SEIS (Seismic Experiment for Interior Structure of Mars) experiment and for the first time in planetary exploration, Mars’ internal structure and seismicity are constrained. ▪ More than 1,300 seismic events and seismic noise records enable the first comparative seismology studies together with Earth and lunar seismic data. ▪ Inversion of seismic travel times and waveforms provided the first interior model of another terrestrial planet, down to the core. ▪ Several impacts were also seismically recorded with their craters imaged from orbit, providing the first data on impact dynamic on Mars.

3D velocity and anisotropy of the southeastern Tibetan plateau extracted by joint inversion of wave gradiometry, ambient noise, and receiver function

With the deformation pattern of the southeastern Tibetan Plateau remaining controversial, its anisotropic structure and the three-dimensional (3D) velocity model are powerful tools for understanding this issue. The surface wave dispersion curve (DC) and receiver function (RF) are sensitive to absolute velocity and elasticity contrasts, respectively. In this work, we propose a new method for the joint inversion of 3D shear velocity (Vs) and anisotropy using the datasets of RFs and DCs, which are obtained by combining the wave gradiometry method and ambient noise tomography. Because DCs and RFs are divided into back-azimuth bins, and both datasets in each bin are jointly inverted to obtain velocity in the corresponding direction, the method is called the azimuth-dependent joint inversion method. By applying this new method to the Temporary West Sichuan Array, we obtained a high-resolution 3D Vs model and anisotropic structure of the southeastern Tibetan Plateau. In our model, the vertically consistent fast propagation direction in the Songpan-Ganzi block (SGB) and the north Chuandian block (NCDB) suggest that the deformation in this region is largely coupled vertically. Our model provides anisotropic evidence of the Yangtze crystalline basement extending to the SGB. The low-velocity zone (LVZ) was much thicker and more prominent in the NCDB than in the SGD. Anisotropy in the middle lower crust of the NCDB does not show any sign of an LVZ flowing from the SGB or the Qiangtang block to the NCDB. These observations are contrary to those of the crust flow model, suggesting that the LVZ was generated locally. The mushroom-shaped high-velocity structure beneath the southern Chuandian block and boundary fault zone may be a relic of the Emeishan paleo-mantle plume.

Moho depth variation and shear wave velocity structure in northern Algeria from joint inversion of P-wave receiver functions and Rayleigh wave dispersion data

SUMMARY In this work, the Moho depth and the velocity structure of the crust and upper mantle beneath broad-band seismic stations of the Algerian broad-band seismic network are investigated. Teleseismic P-wave receiver functions jointly inverted with Rayleigh wave dispersion curves obtained from local earthquakes have been used. The seismic stations are located in different geological settings including the Tell Atlas, High Plateaus and the Saharan Atlas. The crustal thickness and the Vp/Vs ratio are first derived by the H–κ stacking method of receiver functions. The inversion results show the variation in Moho depth in the different geological contexts. The shallowest depths of the Moho (∼20–30 km) are estimated along the Algerian continental margin and Tell Atlas. In the High Plateaus region, the Moho depths vary from 30–36 km, whereas the deepest Moho depths are found in the Saharan Atlas (36–44 km). Two-layer crust is observed in the whole study area. In the upper crust, ∼8–14 km thick, the average shear wave velocity is ∼3.0 km s−1. The lower crust of about 12–30 km thick has an average shear wave velocity that ranges between 3.4 and 3.8 km s−1. The lower crust is thicker than the upper crust particularly in the Saharan Atlas. The upper mantle shear wave velocity varies from 4.1 to 4.5 km s−1 maximum and is stable, generally, below ∼60 km depth. Two low-velocity zones are clearly observed particularly in the eastern part of the Tell Atlas and the High Plateaus. The first one about 10 km thick is in the lower part of the lower crust and the other one is in the upper mantle between 40 and 60 km depth. The obtained results are in accordance with the previous results found in the region, particularly those using land gravity and seismic data. As the first estimate of the Moho depth from earthquake data in northern Algeria, using the receiver function method, this study sheds new light on the crustal structure and the Moho depth in this region of the world.

Stepwise Lithospheric Delamination Leads to Pulsed Cenozoic Uplifts of Central Tien Shan

Abstract The Tien Shan provides an ideal site to study mechanism of intracontinental orogeny due to distant effect of Indo-Asian collision. We investigate lithospheric structures, in particular the lithosphere-asthenosphere boundary (LAB), of Central Tien Shan (CTS) using S wave receiver functions. The results show distinct structures across the orogen. Under the southern CTS, the LAB is shallower than that of the Tarim Basin; a 50 km vertical offset implies that part of the lithosphere has been delaminated. Under the middle CTS, two phases of negative velocity gradient are obtained, which may indicate a new LAB and an ongoing delamination underneath. Under the northern CTS and Kazakh Shield northward, the lithosphere is stable although the LAB inclines southward slightly. The two periods of lithospheric delamination under the southern and middle CTS account well for pulsed uplifts of the Tien Shan at ~11-8 Ma and ~5-0 Ma, respectively.

Joint Inversion of Body Wave Arrival Times, Surface Wave Dispersion Data and Receiver Functions: Method and Application to South China

AbstractWe have developed a new joint inversion method that incorporates body wave arrival times, surface wave dispersion and receiver functions to simultaneously update earthquake locations and constrain three‐dimensional P‐wave (Vp) and S‐wave velocity (Vs) models. Due to complementary sensitivities of the three types of data, the proposed joint inversion algorithm can reduce the intrinsic non‐uniqueness of inversions using fewer types of data and better determine smooth velocity variations and velocity discontinuities. Synthetic tests demonstrate the advantages of this new joint inversion algorithm in resolving velocity structures, especially in constraining velocity gradients across the Moho interface. We have applied the proposed joint inversion algorithm to image the lithosphere velocity structure of south China. The inverted Vp and Vs models fit body wave arrival times, surface wave dispersion and receiver functions well. We further analyzed the distribution of Vs gradients across the Moho interface in detail, which helps us better understand tectonics in south China.

Test of P-wave receiver functions for a seismic velocity and gravity model across the Baikal Rift Zone

SUMMARY The seismic receiver function (RF) technique is widely used as an economic method to image earth's deep interior in a large number of seismic experiments. P-wave receiver functions (RFs) constrain crustal thickness and average Vp/Vs in the crust by analysis of the Ps phase and multiples (reflected/converted waves) from the Moho. Regional studies often show significant differences between the Moho depth constrained by RF and by reflection/refraction methods. We compare the results from RF and controlled source seismology for the Baikal Rift Zone by calculating 1480 synthetic RFs for a seismic refraction/reflection velocity model and processing them with two common RF techniques [H–κ and Common Conversion Point (CCP) stacking]. We compare the resulting synthetic RF structure with the velocity model, a density model (derived from gravity and the velocity model), and with observed RFs. Our results demonstrate that the use of different frequency filters, the presence of complex phases from sediments and gradual changes in the properties of crustal layers can lead to erroneous interpretation of RFs and incorrect geological interpretations. We suggest that the interpretation of RFs should be combined with other geophysical methods, in particular in complex tectonic regions and that the long-wavelength Bouguer gravity anomaly signal may provide effective calibration for the determination of the correct Moho depth from RF results. We propose and validate a new automated, efficient method for this calibration.

Decoding the Emeishan Permian mantle plume in the southeastern margin of Tibet from the seismic signature of the local lithosphere

SUMMARY The Emeishan large igneous province (ELIP), which is located along the southeastern margin of Tibet, is widely regarded as the erosional remnants of voluminous Permian mafic volcanic successions that were associated with mantle plume activity. In addition, evidence for the eruption of the mantle plume and its signature of lithospheric modification should have been retained in the lithosphere after cooling of the plume. Here we use calculated compressional- and shear wave receiver functions (PRFs and SRFs, respectively) from >750 broad-band stations, together with heat flow and Bouguer gravity anomaly data, to investigate the seismic signature within the lithosphere beneath the southeastern margin of Tibet. Our results from the inner zone (INZ) of the ELIP indicate that the ∼60 km depth of the Moho interpreted from depth-migrated PRF images is consistent with that inferred from SRFs and is deeper than that in the region surrounding the INZ, which shows a continuous positive converted phase at ∼35–40 km depth. The crustal geometry and physical properties show a 15–20-km-thick mafic underplated layer that extends for 150–200 km in the E–W direction and ∼400 km in the N–S direction along the base of the crust in the northeastern INZ and its surrounding region. We interpret the underplated layer as the concealed Emeishan plume hotspot track and suggest that the head of the postulated mantle plume is centred approximately beneath Panzhihua City and its surroundings. We propose a new geodynamic model for the formation of the ELIP, in which the lower–middle crust has been strengthened by voluminous mafic intrusions and underplating, whereas the Yangtze Craton in the plume region has been destroyed by removing the base of the lithosphere.

A distinct contrast in the lithospheric structure and limited crustal flow across the northeastern Tibetan Plateau: Evidence from Vs and Vp/Vs imaging

The uplift mechanism and geodynamic model of the northeastern (NE) Tibetan Plateau remain controversial. Two competing models have been proposed for the uplift of the NE Tibetan Plateau, including crustal shortening through folding or thrust faulting and the middle-lower crustal flow model. Here, we applied the joint inversion of surface wave dispersions and receiver functions with P-wave velocity constraints on a dense linear broadband seismic array to obtain the Vs and Vp/Vs profiles of the crust and uppermost mantle across northeastern Tibet. The inversion yields robust structural images that show stark Vs and Vp/Vs jumps across the West Qinling fault from the Songpan-Ganzi and Kunlun-West Qinling terranes (south) to the Qilian and North China terranes (north). The results reveal that regional low-Vs/high-Vp/Vs anomalies present in the middle-to-lower crust, which may be the signature of viscous flow, are distinctly limited south of the West Qinling fault. In contrast, the imaged physical properties of the crust and mantle lithosphere north of this fault manifest a rigid pattern, which is inferred to be originated from the North China craton underthrust below the northern margin of the Tibetan plateau. Our seismic imaging indicates that the NE corner of the Tibetan Plateau does not currently serve as the flow channel of the ductile material from the central and northern Tibet.

Nature of the lithosphere-asthenosphere boundary beneath the Eastern Dharwar Craton of the Indian Shield

• Thickness and nature of lithosphere-asthenosphere boundary has been discussed. • Partial melts may explain the shallow and sharp lithosphere-asthenosphere boundary. • Plume induced lithospheric thinning argued beneath the Eastern Dharwar Craton . • The Eastern Dharwar Craton lithosphere has been substantially modified. The lithosphere-asthenosphere boundary (LAB) is a fundamental element of the plate tectonic hypothesis that accommodates the differential motion of rigid lithosphere over the weaker asthenosphere . In recent times, various usages have been used to define the LAB, depending on the nature of their measurements. Here, we investigate the lithospheric structure beneath the Eastern Dharwar Craton (EDC) of the Indian Shield using geochemical, thermal and seismological data sets. We analysed S-receiver functions from the stations deployed in the EDC along with the surface wave dispersion tomography. We also added thermal measurements from 5 different locations and geochemical data from 34 Kimberlite/Lamproite xenolith samples to constrain the nature of the LAB. The seismological measurements using Rayleigh wave dispersion and receiver function analysis indicate the lithospheric thickness of 98–118 and 94–118 km respectively, with sharp transition across the LAB. The P-T results from xenoliths are interpreted in concurrence with the heat-flow measurements suggesting a thick thermal lithosphere of ∼ 200 km for the normal mantle solidus with cold geotherm. To reconcile our observations, we invoke partial melts or melts enriched in volatiles, which significantly lowers the viscosity of mantle rocks inducing a zone of weakness between the rigid lithosphere (∼125 km) and the convective asthenosphere. Further, we favour the view that the thick lithosphere beneath the Indian plate has been thinned by a plume during the Gondwanaland breakup at ∼ 130 Ma. The presence of younger kimberlites from the Indian shield support that it is further degenerated by the delamination leading to an uneven topography in the LAB.