Crustal and uppermost mantle structure of the northeastern Qinghai-Xizang Plateau from joint inversion of surface wave dispersions and receiver functions with P velocity constraints

Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a joint inversion of receiver functions and surface wave dispersions with P-wave velocity constraints using data from the Chin Array Ⅱ temporary stations deployed across the Qinghai-Xizang Plateau. Prior to joint inversion, we applied the H-κ-c method(Li JT et al., 2019) to the receiver function data in order to correct for the back-azimuthal variations in the arrival times of Ps phases and crustal multiples caused by crustal anisotropy and dipping interfaces. High-resolution images of vS, crustal thickness, and vP/vSstructures in the Qinghai-Xizang Plateau were simultaneously derived from the joint inversion. The seismic images reveal that crustal thickness decreases outward from the Qinghai-Xizang Plateau. The stable interiors of the Ordos and Alxa blocks exhibited higher velocities and lower crustal vP/vSratios. While, lower velocities and higher vP/vSratios were observed beneath the Qilian Orogen and Songpan-Ganzi terrane(SPGZ), which are geologically active and mechanically weak, especially in the mid-lower crust.Delamination or thermal erosion of the lithosphere triggered by hot asthenospheric flow contributes to the observed uppermost mantle low-velocity zones(LVZs) in the SPGZ. The crustal thickness, vS, and vP/vSratios suggest that whole lithospheric shortening is a plausible mechanism for crustal thickening in the Qinghai-Xizang Plateau, supporting the idea of coupled lithospheric-scale deformation in this region.

Crustal structure in the Anyuan Coal Mine and its adjacent areas of Jiangxi Province by P-wave receiver functions

We collected high-quality teleseismic events recorded by 12 broadband seismographs deployed in the Anyuan Coal Mine and its adjacent areas in Pingxiang City,Jiangxi Province for nearly two years.The H-κ-c stacking method was employed to obtain the crustal thickness and Poisson's ratio distribution,then the characteristics of crustal structure below the stations were obtained by using the time-domain linear inversion method.The crustal thickness in the Anyuan Coal Mine and its adjacent areas ranges from approximately 32~35 km,with an average thickness of 33 km,which is consistent with the crustal thickness results in South China from previous studies using the receiver function method.The average Poisson's ratio of the crustal bulk composition in the study area varies between 0.22 and 0.25,which is lower than the global value with a 0.27 average,indicating a predominantly intermediate-acidic or felsic crustal composition.There is a weak negative correlation between Poisson's ratio and crustal thickness estimates in the Anyuan Coal Mine and its adjacent areas,suggesting that the absence of mafic-ultramafic materials in the lower crust is associated with the process of crustal delamination.The velocity inversion results indicate that the crustal structure including three velocity discontinuity interfaces,with the first at a depth of approximately 1.5 km,the second at about 10~15 km,and the third being the Moho.The study also indicates that the results obtained by the H-κ-c stacking method are significantly better than those obtained by the H-κmethod,effectively reducing the standard deviation and dispersion of crustal thickness and vP/vSratio.

Constraining the crustal structure under the central and western Tian Shan based on teleseismic receiver functions and gravity anomalies

The Tian Shan is a vast range that spans several countries in Asia.Understanding its evolutionary history may provide valuable insights into intracontinental orogenic dynamics.In this study,we explored the crustal characteristics of the Tian Shan and their relationships to the tectonic evolution of the region.A new H-stacking method that combines the P receiver function and gravity anomalies was used to estimate the thickness and ratio of P-to S-wave velocities(Vp/Vs)for 91 broadband seismic stations in the central and western Tian Shan.Our results revealed significant lateral variations in crustal thickness and Vp/Vs.A—45-km-thick crust and an intermediate-high Vp/Vs(-1.74-1.84)were found in the Kazakh Shield and Tarim Basin,which we interpreted to indicate a mafic crystalline basement and lower crust.The central Tian Shan varied greatly in crustal thickness(40-64 km)and Vp/Vs ratio(1.65-2.00).which may be due to crustal shortening,mafic underplating,and crustal melting.In contrast,we observed a relatively thin crust(42-50 km)with an intermediate Vp/Vs ratio(-1.78)in the western Tian Shan.The differences in the crustal structures between the western and central Tian Shan imply that the Talas-Fergana Fault may be trans-lithospheric.

Receiver function study of the crustal structure of Northeast China:Seismic evidence for a mantle upwelling beneath the eastern flank of the Songliao Basin and the Changbaishan region

Northeast China is a unique place to study intra-plate volcanism. We analyzed P-wave receiver function data recorded by 111 permanent broadband seismic stations in Northeast China. The results show that the crustal thickness varies from 27.9 km beneath the eastern flank of the Songliao Basin to 40.7 km beneath the Great Xing＇an Range region. The large depth variations of the Moho can be largely but not completely explained by surface topography. The residual Moho depth calculated based on the Airy＇s isostasy model indicates that the Moho is dynamically uplifted by 〈3 km beneath the eastern flank of the Songliao Basin and the Changbaishan region. We suggest that a mantle upwelling, which has been proposed by several recent seismic studies, might have caused the uplift.

Stepwise joint inversion of surface wave dispersion,Rayleigh wave ZH ratio,and receiver function data for 1D crustal shear wave velocity structure

Accurate determination of seismic velocity of the crust is important for understanding regional tectonics and crustal evolution of the Earth. We propose a stepwise joint linearized inversion method using surface wave dispersion, Rayleigh wave ZH ratio （i.e., ellipticity）, and receiver function data to better resolve 1D crustal shear wave velocity （Vs） structure. Surface wave dispersion and Rayleigh wave ZH ratio data are more sensitive to absolute variations of shear wave speed at depths, but their sensi- tivity kernels to shear wave speeds are different and complimentary. However, receiver function data are more sensitive to sharp velocity contrast （e.g., due to the existence of crustal interfaces） and Vp/Vs ratios. The stepwise inversion method takes advantages of the complementary sensitivities of each dataset to better constrain the Vs model in the crust. We firstly invert surface wave dispersion and ZH ratio data to obtain a 1D smooth absolute vs model and then incorporate receiver function data in the joint inver- sion to obtain a finer Vs model with better constraints on interface structures. Through synthetic tests, Monte Carlo error analyses, and application to real data, we demonstrate that the proposed joint inversion method can resolve robust crustal Vs structures and with little initial model dependency.

Imaging the mantle transition zone beneath eastern and central China with CEArray receiver functions

We analyzed a total of 37 427 receiver-function data recorded by national and regional broadband seismic networks of the China Earthquake Administration to study the mantle transition zone beneath eastern and central China. Significant topography on both the 410-km and 660-km discontinuities was clearly imaged in the 3D volume of CCP （common-conversion-point） stacked images that cover an area of 102.5°F-122.5°E and 22.0°N-42.0°N. 3D crustal and mantle velocity models were used in computing the Ps time moveouts to better image the absolute depths of the two discontinuities. We found that the 660-km discontinuity is depressed up to 25 km along the east coast of China, indicating the presence of the subducted Pacific slab in the region. More interestingly, a double 660-km was observed beneath the Huanghai Sea, which was not seen from the 410-km discontinuity. We found a strong and localized anomaly beneath the Quaternary Datong volcano located at the northeastern edge of the Ordos Plateau in North China. The 410-km is depressed by as much as 15 km beneath the volcano. If the amount of the depression is caused by unmodeled low velocity anomaly in the upper mantle, it requires either a strong （4%-8%） localized low velocity anomaly near the lithosphere-asthenosphere boundary or -1.5% low velocity column across the entire upper mantle above the 410-km. The latter could suggest a deep origin of the Datong volcano, although it was usually attributed to edge-driven convection by most of the literatures.

S-wave velocity structure beneath Changbaishan volcano inferred from receiver function

The S wave velocity structure in Changbaishan volcanic region was obtained from teleseismic receiver function modeling. The results show that there exist distinct low velocity layers in crust in volcano area. Beneath WQD station near to the Tianchi caldera the low velocity layer at 8 km depth is 20 km thick with the lowest S-wave velocity about 2.2 km/s At EDO station located 50 km north of Tianchi caldera, no obvious crustal low velocity layer is detected. In the volcanic region, the thickness of crustal low velocity layer is greater and the lowest velocity is more obvious with the distance shorter to the caldera. It indicates the existence of the high temperature material or magma reservoir in crust near the Tianchi caldera. The receiver functions and inversion result from different back azimuths at CBS permanent seismic station show that the thickness of near surface low velocity layer and Moho depth change with directions. The near surface low velocity layer is obviously thicker in south direction. The Moho depth shows slight uplifting in the direction of the caldera located. We con- sider that the special near surface velocity structure is the main cause of relatively lower prominent frequency of volcanic earthquake waveforms recorded by CBS station. The slight uplifting of Moho beneath Tianchi caldera indicates there is a material exchanging channel between upper mantle and magma reservoir in crust.

Crust and upper mantle structures beneath Northeast China from receiver function studies

P-wave and S-wave receiver function analyses have been performed along a profile consisted of 27 broadband seismic stations to image the crustal and upper mantle discontinuities across Northeast China. The results show that the average Moho depth varies from about 37 km beneath the Daxing＇anling orogenic belt in the west to about 33 km beneath the Songliao Basin, and to about 35 km beneath the Changbai mountain region in the east. Our results reveal that the Moho is generally flat beneath the Daxing＇anling region and a remarkable Moho offset （about 4 km） exists beneath the basin-mountain boundary, the Daxing＇anling-Taihang Gravity Line. Beneath the Tanlu faults zone, which seperates the Songliao Basin and Changbai region, the Moho is uplift and the crustal thick- ness changes rapidly. We interpret this feature as that the Tanlu faults might deeply penetrate into the upper mantle, and facilitate the mantle upwelling along the faults during the Cenozoic era. The average depth of the lithosphereasthenosphere boundary （LAB） is - 80 km along the profile which is thinner than an average thickness of a continental lithosphere. The LAB shows an arc-like shape in the basin, with the shallowest part approximately beneath the center of the basin. The uplift LAB beneath the basin might be related to the extensive lithospheric stretching in the Mesozoic. In the mantle transition zone, a structurally complicated 660 km discontinuity with a maximum 35 km depression beneath the Changbai region is observed. The 35 km depression is roughly coincident with the location of the stagnant western pacific slab on top of the 660 km discontinuity revealed by the recent P wave tomography.

Crustal structure in the Binchuan basin of Yunnan constrained from receiver functions on a 2-D seismic dense array

The Binchuan region is located in a seismically active area in northwestern Yunnan,China.The detailed crustal structure is important to understand the tectonic evolution and to assess the seismic hazard in the study area.With a 2-D dense array deployed in this region,we use teleseismic receiver function traditional imaging methods,including the H-κand common-conversion-point stacking methods,to derive high-resolution crustal thickness and vP/vS ratio maps.Our results indicate that the crustal thickness increases from~40 km to~46 km in the south-north direction,and the average crustal thickness beneath the Binchuan basin is~42 km.Our results agree with previous results but have higher resolution due to dense interstation spacing.

Crustal structure beneath Qiangtang and Lhasa terrane from receiver function

To determine the crustal structure in central Tibet, we used teleseismic waveform data recorded by 18 stations in the INDEPTH-Ⅲ seismic array across the central Tibet from the central Lhasa terrane to the central Qiangtang terrane. The S-wave velocity structures beneath stations are determined by inverting the stacked radial receiver function using the GA method. The first order features in the receiver function are modeled. Our results show that the Moho in Qiangtang is about 8 km shallower than that in Lhasa terrane along the INDEPTH-Ⅲ profile. It maybe suggests the northward subduction of the Lhasa mantle lid beneath the Qiangtang terrane is affected by the India-Asia collision. We conclude that there exist low velocity zone in the middle crust across the northern Lhasa and Qiangtang terrane, which can be related to the high temperature upper mantle beneath that.

A study on deep structure using teleseismic receiver function in Western Yunnan

s Western Yunnan is located at the boundary of collision or underthrusting zone of Eurasian plate and is influenced by many times tectonic movements. With very complex geological environment and tectonic background, it is one of the seismically active areas. In the paper, the teleseismic records were selected from 16 national, local and mo-bile stations, including 4 very-wide-band mobile stations of PASSCAL. And nearly 2 000 receiver functions were extracted. Two measuring lines are 650 km and 450 km, respectively and across some major tectonic units in Western Yunnan. It is indicated that Nujiang might be a seam characterized by underthrusting. The western and eastern boundaries of Sichuan-Yunnan rhombus block, i.e., Honghe and Xiaojiang faults, might be an erection seam or collision belt. Panxi tectonic zone still has the characteristics of continental rift valley, that is, the surface is hollow and the upper mantle is upwarping. The tectonic situation in Western Yunnan is of certain regulation with the interlacing distribution of orogenic zone and seam. The crustal thickness decreases gradually from the north to the south and the S wave velocity is globally lower here.

Crust and upper mantle structure beneath Bohai Sea inferred from receiver function study

We analyzed teleseismic waveforms recorded by 36 stations near Bohai Sea region and obtained 2 248 high quality receiver functions.The crustal thickness (H) and average crustal vP/vS ratio (κ) as well as the Poisson's ratios beneath 34 stations were estimated using the H-κ stacking method.The results indicate that crustal thicknesses near the Liaoning province range from 30.0 to 35.5 km,and the corresponding vP/vS ratios vary from 1.72 to 1.89 which corresponds to Poisson's ratio with a range from 0.243 to 0.305.We also apply a common conversion point (CCP) stacking method of receiver function (RF) to image the upper mantle discontinuity structure beneath Bohai Sea region.Both the 410-km and the 660-km discontinuities (hereafter called the 410 and the 660) show clearly in the study region.The transition zone (TZ) thickness shows a different picture from the west to the east of the study region,which is a little bit thicker than that of the global average in the west of longitude 122°E,however,thinner in the east of longitude 122°E.We suggested that the dehydration of sinking slab into the lower mantle or a small-scale mantle plume from the lower mantle generated hot upwelling beneath this region.

Preliminary analysis of teleseismic receiver functions of the Ningxia and its adjacent area

The teleseismic receiver functions of digital seismic network of Ningxia and its adjacent area are calculated with two different Gauss filter factors. The accuracy and stability of the receiver functions are discussed. The h-k stacking method is applied to estimate the crustal thickness and velocity ratio beneath seismic stations. The results indicate that there are sharp changes of crustal thickness and velocity ratio in the studied region. This region is located in the northeastern margin of Tibet, and also a junction of several first-grade blocks. The large contrast of crustal structure in this region is considered to be resulted from the interaction of these blocks. Our results are helpful to construct the completed model of the formation and evolution of the Tibet. Some local structures are also discussed combining with the geological faults.

A new approach for inversion of receiver function for crustal structure in the depth domain

A method for reconstructing crustal velocity structure using the optimization of stacking receiver function amplitude in the depth domain,named common conversion amplitude(CCA)inversion,is presented.The conversion amplitude in the depth domain,which represents the impedance change in the medium,is obtained by assigning the receiver function amplitude to the corresponding conversion position where the P-to-S conversion occurred.Utilizing the conversion amplitude variation with depth as an optimization objective,imposing reliable prior constraints on the structural model frame and velocity range,and adopting a stepwise search inversion technique,this method efficiently weakens the tendency of easily falling into the local extremum in conventional receiver function inversion.Synthetic tests show that the CCA inversion can reconstruct complex crustal velocity structures well and is especially suitable for revealing crustal evolution by estimating diverse velocity distributions.Its performance in reconstructing crustal structure is superior to that of the conventional receiver function imaging method.

Receiver Function Estimated by Wiener Filtering

Wiener filtering is used to estimate receiver function in a time_domain. With the vertical component of 3_component teleseismic P waveform as the input of a Wiener filter, receiver function as the filter response, and radial and tangential components as the expected output, receiver function is estimated by minimizing the error between expected and actual outputs. Receiver function can be obtained by solving the Toeplitz equation using the Levinson algorithm. The non_singularity of the Toeplitz equation ensures the stability of Wiener Deconvolution. Both synthetic and observational seismogram checks show that Wiener Deconvolution is an effective time_domain method to estimate receiver function from teleseismic P waveform.

Crustal structure beneath the central and western North China from receiver function analysis

The North China Craton(NCC)is one of the oldest cratons on earth.Several important tectonic transformations of MesozoicCenozoic tectonic regime led to the destruction of the North China craton.The knowledge of crustal structure can provide important constraints for the formation and evolution of cratons.New maps of sediment thickness,crustal thickness(H)and v_(P)/v_(S)(κ)in the central and western NCC were obtained using sequential H-κstacking.P-wave receiver functions are calculated using teleseismic waveform data recorded by 405 stations from Chin Array project.Benefiting from the densely distribution of temporary seismic stations,our results reveal details of the crustal structure in the study area.The thickness of sedimentary layer in North China ranges from 0–6.4 km,and the thickest sedimentary layer is in Ordos block and its surroundings(about 2.8–6 km);The thickness of sedimentary layer in the Mongolia fold belt and Yinshan orogenic belt is relatively thin(less than 1 km).The crustal thickness of the study area varies between 27–48 km,of which the crust of the North China Plain is about 30–33 km,the central NCC is about 33–40 km,and the Ordos block is 40–48 km thick.The average v_(P)/v_(S)ratios in the study area is mostly between 1.66 and 1.90,and that in the Yanshan-Taihang mountain fold belt is between 1.70 and 1.85,and that in the Ordos block is between 1.65 and 1.90,with an average value of 1.77,indicating the absence of a thick basaltic lower crust.The obvious negative correlation between crustal thickness and average v_(P)/v_(S)ratio within Ordos and Central Asia orogenic belt may be related to magmatic underplating during the crustal formation.There is no significant correlation between the crustal thickness and the v_(P)/v_(S)ratio in the Lüliang-Taihang mountain fold belt,which may be related to the multiple geological processes such as underplating and crustal extension and thinning in this area.The lack of correlation between crust thickness and topography in the central orogenic belt and the North China Basin indicates the topography of these areas are controlled not only by crustal isostatic adjustment but also by the lithospheric mantle processes.

Receiver function estimated by maximum entropy deconvolution

Maximum entropy deconvolution is presented to estimate receiver function, with the maximum entropy as the rule to determine auto-correlation and cross-correlation functions. The Toeplitz equation and Levinson algorithm are used to calculate the iterative formula of error-predicting filter, and receiver function is then estimated. During extrapolation, reflective coefficient is always less than 1, which keeps maximum entropy deconvolution stable. The maximum entropy of the data outside window increases the resolution of receiver function. Both synthetic and real seismograms show that maximum entropy deconvolution is an effective method to measure receiver function in time-domain.

Mapping crustal S-wave velocity structure with SV-component receiver function method

In this article, we analyze the characters of SV-component receiver function of teleseismic body waves and its advantages in mapping the S-wave velocity structure of crust in detail. Similar to radial receiver function, SV-component receiver function can be obtained by directly deconvolving the P-component from the SV-component of teleseismic recordings. Our analyses indicate that the change of amplitude of SV-component receiver function against the change of epicentral distance is less than that of radial receiver function. Moreover, the waveform of SV-component receiver function is simpler than the radial receiver function and gives prominence to the PS converted phases that are the most sensitive to the shear wave velocity structure in the inversion. The synthetic tests show that the convergence of SV-component receiver function inversion is faster than that of the radial receiver function inversion. As an example, we investigate the S-wave velocity structure beneath HIA sta-tion by using the SV-component receiver function inversion method.

S-wave velocity structure inferred from receiver function inversion in Tengchong volcanic area

Tengchong volcanic area is located near the impinging and underthrust margin of India and Eurasia plates. The volcanic activity is closely related to the tectonic environment. The deep structure characteristics are inferred from the receiver function inversion with the teleseismic records in the paper. The results show that the low velocity zone is influenced by the NE-trending Dayingjiang fault. The S-wave low velocity structure occurs obviously in the southern part of the fault, but unobviously in its northern part. There are low velocity zones in the shallow po-sition, which coincides with the seismicity. It also demonstrates that the low velocity zone is directly related to the thermal activity in the volcanic area. Therefore, we consider that the volcano may be alive again.

Wavefield features of teleseismic receiver function in laterally inhomogeneous media

In this study we systematically investigate the effects of laterally inhomogeneous crustal structure with planar dipping interfaces on the teleseismic receiver function. For this purpose the synthetic receiver functions in the 3-D laterally inhomogeneous media are computed using the Maslov asymptotic ray method (CHEN, LIU, 1999). Our results demonstrate that the symmetric and asymmetric variations above the sources azimuths appear on the radial and transverse components of the receiver function, respectively. The dipping direction and dipping angle of the interface beneath the station can be estimated from these features. But, when several different dipping interfaces exist within the crust, only an apparent dipping direction can be inferred from the waveform variations of the receiver function over different azimuths at a single station. By analyzing our observational data it is shown that the planar dipping interface is a well first-order approximation for investigating the crustal lateral inhomogeneity, when only data recorded at a single station are used. And our results given in this study are consistent with those given by seismic array study.