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.展开更多
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.展开更多
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 algorith...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.展开更多
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 in...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.展开更多
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 functio...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.展开更多
文摘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.
基金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.
基金State Natural Science Foundation of China (49974021).
文摘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.
文摘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.
基金State National Natural Science Foundation (49674213).
文摘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.
