The Red Sea-Gulf of Suez-Cairo-Alexandria Clysmic-Trend in northern Egypt is the main earthquake zone in the country,with a moderate-to-high seismic hazard and a history of significant earthquakes caused by rifting an...The Red Sea-Gulf of Suez-Cairo-Alexandria Clysmic-Trend in northern Egypt is the main earthquake zone in the country,with a moderate-to-high seismic hazard and a history of significant earthquakes caused by rifting and active faulting.To improve our understanding of the tectonic and seismic processes in this area,more comprehensive imaging of the crustal structure is required.This can be achieved by increasing the number of receiver functions(RFs)recorded by the seismic stations in northern Egypt and the southeastern Mediterranean.Data handling and processing should also be automated to increase process efficiency.In this study,we developed a capsule neural network for automated selection of RFs.The model was trained on a dataset containing RFs(both selected and unselected)from five broadband stations in northern Egypt.Stations SLM,SIWA,KOT,NBNS,and NKL are located in the unstable shelf region of Egypt,where limited knowledge of the deep crustal structure is available.The proposed capsule neural network achieved an average precision of 80%on the test set.The automated selection of RFs using a capsule neural network has the potential to significantly improve the efficiency and accuracy of RF analysis,as demonstrated by the stacking test.This could lead to a better understanding of crustal structure and tectonic processes in northern Egypt and the southeastern Mediterranean.展开更多
Based on the observational data from 60 short-period stations deployed in the Jishishan M6.2 earthquake epicenter and adjacent regions(Gansu Province,2023),this study inverted the near-surface S-wave velocity structur...Based on the observational data from 60 short-period stations deployed in the Jishishan M6.2 earthquake epicenter and adjacent regions(Gansu Province,2023),this study inverted the near-surface S-wave velocity structure through teleseismic receiver function analysis by using the amplitude of direct P-wave.The results reveal that the epicentral area(Liugou Township and surroundings)exhibits markedly low S-wave velocities of 400-600 m/s,with a mean value of(500±50)m/s.In contrast,intermountain basins-Guanting Basin and Dahejia Basin-demonstrate significantly elevated velocities,exceeding the epicentral zone by 100-300 m/s,with values concentrated at 600-900 m/s.Notably,localized areas such as Jintian Village and Caotan Village maintain stable S-wave velocities of(700±30)m/s.The western margin tectonic belt of Jishishan displays distinctive velocity differentiation:A pronounced velocity gradient zone along the 35.8°N latitude boundary separates northern areas(<550 m/s)from southern regions(>750 m/s).These findings demonstrate significant spatial heterogeneity in shallow S-wave velocity structures,primarily controlled by three factors:(1)topographic-geomorphic units,(2)stratigraphic lithological contrasts,and(3)anthropogenic modifications.The persistent low-velocity anomalies(<600 m/s)in the epicentral zone and northern Yellow River T2 terrace likely correlate with Quaternary unconsolidated sediments,enhanced groundwater circulation,and bedrock weathering.These results provide critical geophysical constraints for understanding both the seismogenic environment of the Jishishan earthquake and its damage distribution patterns.Furthermore,they establish a foundational framework for regional seismic intensity evaluation,site amplification analysis,and secondary hazard risk assessment.展开更多
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 j...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.展开更多
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 me...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.展开更多
The Sumatra subduction zone is located in the southwest of the Suna plate,between the Euro-Asia Plate and Indo-Austrilian Plates.With the obliquely subducting of the Indo-Austrilian Plate toward the the Euro-Asia Plat...The Sumatra subduction zone is located in the southwest of the Suna plate,between the Euro-Asia Plate and Indo-Austrilian Plates.With the obliquely subducting of the Indo-Austrilian Plate toward the the Euro-Asia Plate,complex tectonics,strong earthquakes and volcanoes have been observed in this area which has become a well experimental field used to study the subduction zone.In this work,we employed receiver function method to evaluate the S-wave velocity structure beneath 5 broadband seismic stations along the Sumatra subduction zone.We selected 332 receiver function waveforms with intelligent software and manual picking methods,including 130,34,42,29 and 97 receiver function waveforms corresponding to BKNI,GSI,LHMI,MNAI,and PMBI sta-tions,respectively.These stacked receiver function waveforms were applied to inversion to estimate Swave structure beneath each station based on a Neighborhood Algorithm(NA).Our results indicate that the sediment layers for GSI,LHMI and MNAI stations are more than 3 km thick,two stations of which are thicker than 6 km(e.g.GSI and LHMI).The difference of receiver function waveforms for NE,SW and W orientation at station GSI where is accompanied with strong thrust earthquakes suggests that there is a complicated structure beneath this station.Station BKNI and PMBI are located on the eastern side of the Sumatra fault and the thickness of their sediment layers is only~1 km.The crustal thickness for back-arc basin is within 30-36 km.However,the crustal thickness of forearc area varies from~26 km of the forearc ridge to 26-30 km of the forearc basin toward continent and its,which suggests that the down dip limit(slab-Moho intersection)of seismogenic zone is within 29-36 km in forearc and explains why the shallow-focus earthquakes play a dominant role in this area.The stable state for the inner wedge of forearc within a seismogenic circle provides a favorable environment for storing stress.Meanwhile,these faults caused by the subducting of Indo-Austrilian Plate constructed a condition(e.g.cracking of intact rocks and frictional sliding)in which it would trigger shallow-focus seismic activities(releasing stress).展开更多
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 functio...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.展开更多
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 th...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.展开更多
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 ...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.展开更多
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 ...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.展开更多
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 ...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.展开更多
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 sho...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.展开更多
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 a...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.展开更多
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 cha...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.展开更多
Yunnan Province is located on the southeastern margin of Tibet and represents an important marker in understanding the tectonic evolution of Tibetan Plateau. In this study, we calculated teleseismic P-wave receiver fu...Yunnan Province is located on the southeastern margin of Tibet and represents an important marker in understanding the tectonic evolution of Tibetan Plateau. In this study, we calculated teleseismic P-wave receiver functions at 49 permanent broadband seismic stations in Yunnan Province and estimated crustal thickness and the bulk crust ratios of P-wave to S-wave velocities using the H-κ method together with more detailed crustal structural profiles from the common conversion point stacking method. There is a significant transition of Moho interface and lower crustal composition along latitude 26°N in northwestern Yunnan. Decrease of crustal thickness with a concomitant increase of Poisson’s ratio occurs at station CUX. An interesting phenomenon is that a step-like Moho fashion is observed at several stations, which might correspond to local thermal activities, such as partial melt/lower crust delamination. Our results show changes in crustal properties appear to be associated with varieties in upper mantle structure and compositions. We propose the controlling factor of the dynamic processes below 26°N is the result of eastern forward subduction of the Indian Plate; the northern part is controlled by the redirected material flow from the SE Tibet.展开更多
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 Qia...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.展开更多
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, i...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.展开更多
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&...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.展开更多
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 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 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 amp...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.展开更多
We performed a receiver function analysis on teleseismic data recorded along two dense seismic profiles and from 4 broadband regional seismic stations across the northeastern Tibetan plateau.The crustal thickness and ...We performed a receiver function analysis on teleseismic data recorded along two dense seismic profiles and from 4 broadband regional seismic stations across the northeastern Tibetan plateau.The crustal thickness and vP/vS ratio were measured by the H-κdomain search algorithm.The Moho discontinuity across the Haiyuan arc fault zone was also revealed by common conversion point(CCP)imaging.Our study results show that the crustal thickness and the vP/vS ratio were 42–56 km and 1.60–1.88,respectively.The crustal thickening on the northeastern margin indicates that the crust is shortening or that there was a superimposition of crusts during the collision of the Tibetan plateau with Eurasian block.Our results suggest that Haiyuan fault likely resulted from the interactions of high temperature and pressure conditions during the collision of the Indian and Asian continents.The Moho beneath the Haiyuan tectonic region exhibits an obvious offset and a vague discontinuity according to CCP imaging.This study suggests that the Haiyuan arc fault zone is a trans-crustal fault that cuts through the Moho in the northeastern Tibetan Plateau.Moreover,there are indications of strong deformation in the intensive crustal extrusion from the interior of the Tibetan Plateau to its northeastern margin.展开更多
文摘The Red Sea-Gulf of Suez-Cairo-Alexandria Clysmic-Trend in northern Egypt is the main earthquake zone in the country,with a moderate-to-high seismic hazard and a history of significant earthquakes caused by rifting and active faulting.To improve our understanding of the tectonic and seismic processes in this area,more comprehensive imaging of the crustal structure is required.This can be achieved by increasing the number of receiver functions(RFs)recorded by the seismic stations in northern Egypt and the southeastern Mediterranean.Data handling and processing should also be automated to increase process efficiency.In this study,we developed a capsule neural network for automated selection of RFs.The model was trained on a dataset containing RFs(both selected and unselected)from five broadband stations in northern Egypt.Stations SLM,SIWA,KOT,NBNS,and NKL are located in the unstable shelf region of Egypt,where limited knowledge of the deep crustal structure is available.The proposed capsule neural network achieved an average precision of 80%on the test set.The automated selection of RFs using a capsule neural network has the potential to significantly improve the efficiency and accuracy of RF analysis,as demonstrated by the stacking test.This could lead to a better understanding of crustal structure and tectonic processes in northern Egypt and the southeastern Mediterranean.
基金project is supported in part by Broadband Seismic 3D Array Detection(PhaseⅠ),Deep Earth Probe and Mineral Resources Exploration-National Science and Technology Major Project(Grant No.2024ZD1000300)National Natural Science Foundation of China(42204061)Gansu Jishishan 6.2 magnitude earthquake scientific investigation(DQJB23Y45)program。
文摘Based on the observational data from 60 short-period stations deployed in the Jishishan M6.2 earthquake epicenter and adjacent regions(Gansu Province,2023),this study inverted the near-surface S-wave velocity structure through teleseismic receiver function analysis by using the amplitude of direct P-wave.The results reveal that the epicentral area(Liugou Township and surroundings)exhibits markedly low S-wave velocities of 400-600 m/s,with a mean value of(500±50)m/s.In contrast,intermountain basins-Guanting Basin and Dahejia Basin-demonstrate significantly elevated velocities,exceeding the epicentral zone by 100-300 m/s,with values concentrated at 600-900 m/s.Notably,localized areas such as Jintian Village and Caotan Village maintain stable S-wave velocities of(700±30)m/s.The western margin tectonic belt of Jishishan displays distinctive velocity differentiation:A pronounced velocity gradient zone along the 35.8°N latitude boundary separates northern areas(<550 m/s)from southern regions(>750 m/s).These findings demonstrate significant spatial heterogeneity in shallow S-wave velocity structures,primarily controlled by three factors:(1)topographic-geomorphic units,(2)stratigraphic lithological contrasts,and(3)anthropogenic modifications.The persistent low-velocity anomalies(<600 m/s)in the epicentral zone and northern Yellow River T2 terrace likely correlate with Quaternary unconsolidated sediments,enhanced groundwater circulation,and bedrock weathering.These results provide critical geophysical constraints for understanding both the seismogenic environment of the Jishishan earthquake and its damage distribution patterns.Furthermore,they establish a foundational framework for regional seismic intensity evaluation,site amplification analysis,and secondary hazard risk assessment.
基金supported by the Natural Science Basic Research Program of Shaanxi(No.2023-JC-QN-0306)the Special Fund of the Institute of Geophysics,China Earthquake Administration(No.DQJB21B32)the National Natural Science Foundation of China(No.42174069).
文摘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.
基金supported by the Spark Program of Earthquake Technology of CEA,China(XH20032)open Research Project from the State Key Laboratory of Nuclear Resources and Environment,East China University of Technology(2022NRE17)+1 种基金open Research Project from the State Key Laboratory of Geological Processes,Mineral Resources,China University of Geosciences(GPMR202114)the Graduate Student Innovation Fund of East China University of Technology(YC2022-s628)。
文摘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.
基金supported by NSFC(42306084)Hunan Provincial Natural Science Foundation of China(2023JJ40223)Scientific Research Fund of Hunan Provincial Education Department(21B0566).
文摘The Sumatra subduction zone is located in the southwest of the Suna plate,between the Euro-Asia Plate and Indo-Austrilian Plates.With the obliquely subducting of the Indo-Austrilian Plate toward the the Euro-Asia Plate,complex tectonics,strong earthquakes and volcanoes have been observed in this area which has become a well experimental field used to study the subduction zone.In this work,we employed receiver function method to evaluate the S-wave velocity structure beneath 5 broadband seismic stations along the Sumatra subduction zone.We selected 332 receiver function waveforms with intelligent software and manual picking methods,including 130,34,42,29 and 97 receiver function waveforms corresponding to BKNI,GSI,LHMI,MNAI,and PMBI sta-tions,respectively.These stacked receiver function waveforms were applied to inversion to estimate Swave structure beneath each station based on a Neighborhood Algorithm(NA).Our results indicate that the sediment layers for GSI,LHMI and MNAI stations are more than 3 km thick,two stations of which are thicker than 6 km(e.g.GSI and LHMI).The difference of receiver function waveforms for NE,SW and W orientation at station GSI where is accompanied with strong thrust earthquakes suggests that there is a complicated structure beneath this station.Station BKNI and PMBI are located on the eastern side of the Sumatra fault and the thickness of their sediment layers is only~1 km.The crustal thickness for back-arc basin is within 30-36 km.However,the crustal thickness of forearc area varies from~26 km of the forearc ridge to 26-30 km of the forearc basin toward continent and its,which suggests that the down dip limit(slab-Moho intersection)of seismogenic zone is within 29-36 km in forearc and explains why the shallow-focus earthquakes play a dominant role in this area.The stable state for the inner wedge of forearc within a seismogenic circle provides a favorable environment for storing stress.Meanwhile,these faults caused by the subducting of Indo-Austrilian Plate constructed a condition(e.g.cracking of intact rocks and frictional sliding)in which it would trigger shallow-focus seismic activities(releasing stress).
基金State Key Basic Research Development and Programming Project (G199804070201) State Natural Science Foundation (40074008).
文摘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.
文摘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.
基金supported by the National Earthquake Science Experiment in Sichuan and Yunnan Provinces of China(#2016 CESE 0201)National Natural Science Foundation of China(#41574034)China National Special Fund for Earthquake Scientific Research in Public Interest(#201508008)
文摘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.
文摘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.
基金supported by National Key Project of Scientific and Technical Supporting Programs Funded by Ministry of Science & Technology of China (2006BAC01B04)Joint Seismological Science Foundation of China (106023)Contribution No. is 09FE3006 of Institute of Geophysics,China Earthquake Administration
文摘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.
文摘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.
基金supported by the National Key R&D Program of China(No.2018YFC1503400)China Earthquake Science Experiment Project,CEA(Nos.2018CSES0101,2018CSES0102 and 2019CSES0107)+2 种基金HKSAR Research Grant Council GRF(No.14305617)CUHK Direct Grant from Faculty of Science,National Natural Scientific Foundation of China(Nos.41974069 and 41790463)Chen Yong Academician Workstation of Yunnan Province in China(No.2014IC007)。
文摘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.
基金supported by the National Science Foundation of China (Nos. U1839210, 41874097, 41474072, and 41874108)
文摘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.
基金supported by the 973 Project of China(No.2013CB733303)the National Natural Science Foundation of China(No.41474093)
文摘Yunnan Province is located on the southeastern margin of Tibet and represents an important marker in understanding the tectonic evolution of Tibetan Plateau. In this study, we calculated teleseismic P-wave receiver functions at 49 permanent broadband seismic stations in Yunnan Province and estimated crustal thickness and the bulk crust ratios of P-wave to S-wave velocities using the H-κ method together with more detailed crustal structural profiles from the common conversion point stacking method. There is a significant transition of Moho interface and lower crustal composition along latitude 26°N in northwestern Yunnan. Decrease of crustal thickness with a concomitant increase of Poisson’s ratio occurs at station CUX. An interesting phenomenon is that a step-like Moho fashion is observed at several stations, which might correspond to local thermal activities, such as partial melt/lower crust delamination. Our results show changes in crustal properties appear to be associated with varieties in upper mantle structure and compositions. We propose the controlling factor of the dynamic processes below 26°N is the result of eastern forward subduction of the Indian Plate; the northern part is controlled by the redirected material flow from the SE Tibet.
基金National Natural Science Foundation of China (40274029) and Joint Seismological Science Foundation of China (105096).
文摘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.
文摘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.
基金supported by the National Natural Science Foundation of China grants 40674028 and 40474022
文摘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.
基金supported by the basic research and development fund from Institute of Earthquake Science,China Earthquake Administration(grant No.2011IESLZ05)the National Natural Science Foundation of China(grant No.40904014)
文摘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.
基金financially supported by the National Natural Science Foundation of China(Grant 91755214).
文摘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.
基金supported by National Natural Science foundation of China(Nos.41474032 and U1839209)the Basic Research Project of the Institute of Earthquake Forecasting,China Earthquake Administration(Nos.2019CSES0102 and 2012IES010101).
文摘We performed a receiver function analysis on teleseismic data recorded along two dense seismic profiles and from 4 broadband regional seismic stations across the northeastern Tibetan plateau.The crustal thickness and vP/vS ratio were measured by the H-κdomain search algorithm.The Moho discontinuity across the Haiyuan arc fault zone was also revealed by common conversion point(CCP)imaging.Our study results show that the crustal thickness and the vP/vS ratio were 42–56 km and 1.60–1.88,respectively.The crustal thickening on the northeastern margin indicates that the crust is shortening or that there was a superimposition of crusts during the collision of the Tibetan plateau with Eurasian block.Our results suggest that Haiyuan fault likely resulted from the interactions of high temperature and pressure conditions during the collision of the Indian and Asian continents.The Moho beneath the Haiyuan tectonic region exhibits an obvious offset and a vague discontinuity according to CCP imaging.This study suggests that the Haiyuan arc fault zone is a trans-crustal fault that cuts through the Moho in the northeastern Tibetan Plateau.Moreover,there are indications of strong deformation in the intensive crustal extrusion from the interior of the Tibetan Plateau to its northeastern margin.