Mantle structure from inter-station Rayleigh wave dispersion and its tectonic implication in western China and neighboring regions

We processed more than 3000 inter-station great circle paths to determine the phase velocity for the fundamental mode of Rayleigh wave, and finally arrived at 110 paths of high quality dispersion data, which show good spatial coverage in western China and neighboring regions. Rayleigh wave phase velocity dispersion model WChina1D was obtained and compared with previous global and regional models. Phase velocity maps from 15 to 120 s were inverted and the maps of 20, 40, 80, and 120 s are presented in this paper. Checkerboard tests show the average lateral resolution in our area of interest is about 7°. Our tomographic results corroborate a prominent low-velocity anomaly lying mainly in the lower crust and uppermost mantle in the Chang Thang terrane. The apparent low-velocity anomaly also appears in the wide area of northeastern Tibet in the crust and upper mantle. The low-velocity area around southeastern Tibet may be created by the southeastern migration of the low-velocity mass of the Tibetan plateau. The eastern Tarim shows structure with higher velocities relative to that of central Tarim. A large-scale low-velocity anomaly is clearly seen in central and western Mongolia. Our high quality measurements were also used to evaluate the CUB global shear velocity model [Shapiro, N., Ritzwoller, M., 2002. Monte-Carlo inversion for a global shear-velocity model of the crust and upper mantle. Geophys. J. Int. 151, 88-105] of the crust and upper mantle. The 40 s Rayleigh phase velocity map predicted from CUB model shows an apparent discrepancy with our measurements in western China and western Mongolia, which implies a higher estimated (about +1-2%) phase velocity model in these regions, probably due to the Gaussian smoothing condition in their tomography inversion.     

基于背景噪声成像方法的新疆呼图壁储气库地区近地表速度结构研究

近年来,背景噪声成像方法在恢复高频面波信号及获取近地表速度结构方面得到了广泛的应用,本文将该方法应用于准噶尔盆地南缘的呼图壁背斜地区的呼图壁储气库.采用储气库及其周边区域22个台站记录的连续背景噪声数据的垂直分量,通过噪声互相关方法获得了台站对之间的瑞利面波经验格林函数,并进一步提取了0.5~1.5 s的基阶瑞利面波群速度频散曲线.首先根据区域平均频散曲线得到了该地区地下数百米的平均一维横波速度结构,然后利用基于面波射线路径追踪的面波频散直接成像方法得到该地区深度为500 m以上的三维横波速度结构.反演结果显示该地区沉积层较厚,整体横波速度值较小（0.4~0.9 km·s^-1）.储气库在地表投影区域的横波速度值较小,这可能是由于抽注水、气引起的沉积岩石裂隙所导致.储气库东北和东南方向均有明显的相对高速区,推测是区域地下水位和地形起伏综合作用的结果.本研究获得的近地表三维速度结构为呼图壁储气库地区的上覆地层物性研究、区域微震精定位、场地效应的评估和去除浅层影响的深部介质成像等研究提供了重要基础.     

Elastic-anelastic Regional Structures for the Iberian Peninsula Obtained from a Rayleigh Wave Tomography and a Causal Uncoupled Inversion

The elastic and anelastic structure of the lithosphere and asthenosphere of the Iberian Peninsula is derived by means of tomographic techniques applied to local phase and group velocities and local attenuation coefficients of Rayleigh wave fundamental mode. The database consists of surface wavetrains recorded at the broadband stations located in the Iberian Peninsula on the occasion of the ILIHA project. Path-averaged phase and group velocities and attenuation coefficients were previously obtained by standard filtering techniques of surface wavetrains and, subsequently, local dispersion curves were computed according to the Yanovskaya-Ditmar formulation. First, a principal component analysis (PCA) and the average linkage (AL) clustering algorithm are applied to these local values in order to classify the Iberian Peninsula in several rather homogeneous domains from the viewpoint of the similarity of the corresponding local dispersion curves, without previous seismotectonic constraints. Second, averaged phase and group velocities and attenuation coefficients representing each homogeneous region are used to derive the respective elastic and anelastic models of the lithosphere and asthenosphere. This purpose is achieved by using the uncoupled causal inversion of phase and group velocities and attenuation coefficients. The main features of the homogeneous regions are discussed by taking as reference the Hercynic, Alpine and Neogene domains of the Iberian Peninsula, and two questions affecting the reliability of the elastic-anelastic models are revised. First, the coherence of the shear-velocity and Q_β⁻¹ models obtained by causal uncoupled inversion for each region is analysed. Second, the influence of the causal phase and group velocities on the shear-velocity models is evaluated by comparing elastic and anelastic models derived from causal uncoupled inversion with those deduced from non-causal inversion.     

用面波方法研究上扬子克拉通壳幔速度结构

本文研究采用单台法和双台法提取了穿越上扬子的基阶面波相速度和群速度频散;通过对提取的面波群速度和相速度频散进行联合反演,得到的1-D SV速度模型显示上扬子块体下地壳S波速度与典型克拉通区域相当,其上地幔顶部80~170 km深处存在高速的岩石圈盖层,较AK135模型要快2%~3%,其岩石圈厚度约为180 km.在上扬子地区,径向各向异性集中分布在300 km以浅的岩石圈与软流圈部分,其中岩石圈部分SH波比SV波波速要快2%~4%,软流圈部分SH波比SV波波速要快3%~5%;Rayleigh波相速度方位各向异性分析结果显示,上扬子块体周期为25~45 s（大致相当于30~70 km深度范围内）的Rayleigh波相速度存在1.8%~2.7%不等的方位各向异性,其快波方向介于147°~174°.我们认为上扬子块体径向各向异性集中分布在岩石圈、软流圈部分,且各向异性随深度变化, 其岩石圈部分各向异性为大陆克拉通化的遗迹,软流圈部分各向异性与现今板块运动相关.     

Shallow Seismic Velocity Structure of the Betic Cordillera (Southern Spain) from Modelling of Rayleigh Wave Dispersion

A detailed dispersion analysis of Rayleigh waves generated by local earthquakes and occasionally by blasts that occurred in southern Spain, was undertaken to obtain the shear-wave velocity structure of the region at shallow depth. Our database includes seismograms generated by 35 seismic events that were recorded by 15 single-component short-period stations from 1990 to 1995. All these events have focal depths less than 10 km and body-wave magnitudes between 3.0 and 4.0, and they were all recorded at distances between 40 and 300 km from the epicentre. We analysed a total of 90 source-station Rayleigh-wave paths. The collected data were processed by standard digital filtering techniques to obtain Rayleigh-wave group-velocity dispersion measurements. The path-averaged group velocities vary from 1.12 to 2.25 km/s within the 1.0-6.0 s period interval. Then, using a stochastic inversion approach we obtained 1-D shear-wave velocity–depth models across the study area, which were resolved to a depth of circa 5 km. The inverted shear-wave velocities range approximately between 1.0 and 3.8 km/s with a standard deviation range of 0.05–0.16 km/s, and show significant variations from region to region. These results were combined to produce 3-D images via volumetric modelling and data visualization. We present images that show different shear velocity patterns for the Betic Cordillera. Looking at the velocity distribution at various depths and at vertical sections, we discuss of the study area in terms of subsurface structure and S-wave velocity distribution (low velocity channels, basement depth, etc.) at very shallow depths (0–5 km). Our results characterize the region sufficiently and lead to a correlation of shear-wave velocity with the different geological units features.     

蒙古中南部地区面波相速度层析成像

本文采用小波变换频时分析技术提取了1893条蒙古中南部地区双台间基阶Rayleigh波相速度频散曲线,通过对提取到的相速度频散进行二维反演,首次重构了蒙古高原中南部10~80s周期内分辨率为0.5°×0.5°的相速度分布图.结果表明,短周期相速度频散(10~20s)明显受地表地形控制,表现为杭爱—肯特山盆地为相对高速异常,而南部的戈壁带表现为相对低速异常;在30s到60s周期内,南部的戈壁带和北部的杭爱—肯特山盆地都显示出稳定的高速异常,而中部戈壁带则表现为低速异常.该低速异常区与新生代火山岩出露位置的一致性,暗示该低速异常可能与新生代火山活动有关.穿过蒙古高原中部的两条主要断裂带与研究区内的地震波速度分布具有很好的一致性,可能暗示两条断裂带一直延伸到整个岩石圈.此外,蒙古中南部地区的壳幔速度低于全球平均大陆值,且其相速度频散曲线与遭受破坏的克拉通,尤其大陆裂谷很相似,可能暗示着研究区具有薄的、活跃的岩石圈构造.     

青藏高原东部和周边地区地壳速度结构的背景噪声层析成像

利用连续地震背景噪声记录和互相关技术获得瑞利面波格林函数,进而反演获得了青藏高原东部和周边地区的地壳三维速度结构.地震数据源于北京大学宽频带流动观测地震台阵,国家数字测震台网数据备分中心提供的部分固定台站的连续记录及INDEPTH IV宽频带流动观测地震台阵.首先对观测数据进行处理和分析取得所有可能台站对的面波经验格林函数和瑞利波相速度频散曲线,反演得到了观测台阵下方周期从6~60s的瑞利波相速度异常分布图像.并且进一步反演获得研究区域三维剪切波速度结构和莫霍面深度分布.短周期(6~14s)相速度异常分布与地表地质构造特征吻合较好,在青藏高原和四川盆地之间存在一个明显的南北向转换带.而本文最重要的结果是周期大于25s的相速度异常分布图像显示,以昆仑断裂带为界,柴达木盆地和祁连山脉地区呈现与青藏高原截然不同的中地壳速度结构,反而与青藏高原东缘地区和川滇菱形块体速度结构相似.反演获得的剪切波速度在27.5~45km深度的切片也明显地揭示:青藏高原的松潘—甘孜地块和羌塘地块呈现均一的低速层;然而,柴达木盆地和祁连山脉地区则呈现较强的横向不均匀性,尤其是柴达木盆地的高速异常和四川盆地的高速异常相对应.这些结果为前人提出的青藏高原东北向台阶式增长模式提供了重要的地震学观测证据.与全球一维平均速度模型(AK135)相比较发现,本文测量和反演获得的研究区域内平均相速度和剪切波速度都比AK135模型慢很多,尤其是青藏高原的中地壳(25~40km)剪切波速度显著低于全球平均速度模型.进一步的层析成像反演证实松潘—甘孜和羌塘地块中地壳(27.5~45km)呈现大范围均一的低速层,为青藏高原可能存在大规模中下地壳"层流"提供地震学观测证据.在祁连山脉的27.5~45km深度观测到的明显低速异常体可能对应于该造山带下地幔岩浆活动导致的底侵作用,表明引起该地区地壳增厚的主要机制可能是来自地幔岩浆的底侵作用.     

利用Rayleigh波相速度和群速度联合反演青藏高原东北缘S波速度结构

利用青藏高原东北缘地区固定和流动地震台网2007年8月到2012年1月期间记录的远震波形,运用小波变换频时分析方法分别测定了1216和653条周期从15到140 s的台站间基阶Rayleigh相速度和群速度频散曲线.通过对上述频散进行反演,重构了青藏高原东北缘分辨率高达0.5°×0.5°的2-D相速度和群速度分布图.然后通过对所提取到的每个格网点Rayleigh波相速度和群速度频散进行联合反演,得到了研究区下方一维S波速度结构.最后通过线性插值,得到了青藏高原东北缘下方地壳上地幔三维S波结构.结果表明,印度板块向北俯冲已经达到班公-怒江缝合带附近;在柴达木盆地北部祁连山下面我们发现了亚洲板块,且其没有表现出明显的向南俯冲的迹象;在两大板块中间,我们观测到延伸到250 km深度的低速异常,该低速异常可能是地幔物质底辟上涌现象造成的.     

声波测井中的相速度与群速度讨论(英文)

声波测井过程中获取的速度到底是相速度还是群速度,目前仍存在一些争议,本文从理论分析和数值模拟的角度,使用三种模型对这一问题进行了研究。首先,构造一个相速度与群速度可调的稳态声波传播模型——不同声速的两个平面波叠加模型,利用慢度时间相关(STC)方法提取声波波速,数值模拟结果表明,无论相速度较大或是群速度较大,STC方法提取出来的波速都是相速度;其次,通过频散分析和割线积分得到刚性壁圆柱流体模型中的频散曲线与分波波形,使用STC方法得到的速度与相速度的频散曲线吻合较好,而直接读取波至获得的速度与群速度的频散曲线趋势一致;最后,利用频散分析和实轴积分方法,获得偶极子在慢地层中激发的模式及全波波形,得到的结果再次验证了刚性壁圆柱流体模型中的结论。     

Crustal velocity structures beneath North China revealed by ambient noise tomography

We collected continuous noise waveform data from January 2007 to February 2008 recorded by 190 broadband and 10 very broadband stations of the North China Seismic Array.The study region is divided into grid with interval 0.25°×0.25°,and group velocity distribution maps between 4 s and 30 s are obtained using ambient noise tomography method.The lateral resolution is estimated to be 20-50 km for most of the study area.We construct a 3-D S wave velocity model by inverting the pure path dispersion curve at each grid using a genetic algorithm with smoothing constraint.The crustal structure observed in the model includes sedimentary basins such as North China basin,Yanqing-Huailai basin and Datong basin.A well-defined low velocity zone is observed in the Beijing-Tianjin-Tangshan region in 22-30 km depth range,which may be related to the upwelling of hot mantle material.The high velocity zone near Datong,Shuozhou and Qingshuihe within the depth range of 1-23 km reveals stable characteristics of Ordos block.The Taihangshan front fault extends to 12 km depth at least.     

辽宁及邻区背景噪声相速度结构研究

收集辽宁及其周边地区（吉林、河北、山东、内蒙）70个宽频带地震仪2012年连续背景噪声波形数据,基于地震背景噪声层析成像方法,得到研究区面波群速度及相速度图像。利用台站对互相关方法,提取瑞利面波格林函数,采用时频分析法（FTAN）获取2 416条相速度频散曲线,从中筛选1 661条信噪比较高的频散曲线。将研究区以0.25°×0.25°进行网格化,采用Ditmar等提出的层析成像反演方法,得到周期10—40 s的瑞利面波群速度及相速度结构分布图。与群速度结果相比,分辨率更高,研究区大部可达0.5°×0.5°（局部可达0.25°×0.25°）。结果表明,辽宁地区地壳及上地幔面波相速度结构存在显著的横向不均匀性。在周期10—15 s的群速度图中,浅层及中上地壳速度分布与研究区地形地貌及主要地质构造单元具有较好的对应关系,盆地及沉积层低速,山区隆起高速,且在高低速转换带多为地震孕震区;在周期20—30 s相速度结构图中,下地壳至上地幔顶部深度范围内,相速度速度结构主要受地壳厚度及渤海湾内巨厚沉积层的影响,在海城至大连区域内出现的低速异常推测为地下热物质上涌;随着深度的增加,在周期30—40 s的相速度图中,速度分布逐渐受控于莫霍面起伏,明显变化出现在辽东半岛,由高速变为低速。     

3-D Velocity Structure of the Middle and Upper Crust in the Yunnan Region,China

Using short-period (1–18 s) surface wave data recorded by 23 stations of the Yunnan Digital Seismic Network of China we determined phase velocities of the fundamental Rayleigh wave along 209 paths by the two-station narrowband filtering and cross-correlation method, followed by an inversion for phase velocity distributions at various periods using the Ditmar-Yanovskaya method. We then obtained a 3-D S-wave velocity structure of the middle and upper crust in the Yunnan region using the genetic algorithm. The results show strong lateral variation of phase velocity in the region. The short-period phase velocity variation is closely related to thickness variation of sedimentary layer in the shallow crust. Within the depth range of 26–30 km, the S-wave velocity in the Sichuan-Yunnan rhombic block is lower than in the surrounding areas. Most large earthquakes of M > 6.0 in Yunnan occurred in the transition zones between low and high velocities.     

中国东北地区噪声层析成像

中国东北地区是中国唯一的深震孕育区,获取该区的壳幔结构,对于研究板块俯冲、深震以及火山活动等有重要的科学意义.本文利用该区159个固定台站2011年1月至2012年6月和27个流动台站2011年1月至2011年6月间的垂向波形连续记录,计算了台站间的预估格林函数,并采用基于连续小波变换的时频分析方法,测量了双台路径上瑞雷波的群速度和相速度频散曲线.通过质量控制和筛选,最终得到了2204条路径上周期5~40 s的群/相速度频散曲线.检测板测试表明,反演结果可以达到2°×2°的分辨.利用Ditmar & Yanovskaya反演方法,我们得到了研究区（105°E-135°E,39°N-52°N）周期8~30 s的瑞雷波的群速度和相速度分布图.不同周期的群/相速度分布图,反映了不同深度S波速度的横向变化情况.研究结果显示：中国东北地区的地壳上地幔S波速度结构存在横向非均匀性.短周期（如8 s）的群/相速度分布与地表地质构造具有明显的相关性,具体来说,山区显示为高速,沉积盆地显示为低速;随着周期的增大（如15 s,22 s）,地形的控制作用相对减弱;较长周期（如30 s）的群/相速度分布与地壳厚度密切相关.     

Great-circle Rayleigh wave attenuation and group velocity,Part III: Inversion of global average group velocities and attenuation coefficients

Average shear-velocity models for the upper mantle have been derived by controlled Monte Carlo inversion of global average Rayleigh wave group velocity (GAGV) data for periods between 50 and 300 seconds. GAGV data have been corrected for attenuative dispersion using a method based on the theory of Liu, Anderson and Kanamori. Two types of model bounds have been used with one- or two-layer low-velocity zones beginning at depths of 70 and 100 km. All models fitting GAGV data within one standard deviation have low-velocity zones in the 100–200 km depth range. Models with low-velocity zones beginning at 70 km, as well as 100 km, fit GAGV data within one standard deviation, so the average thickness of the lithosphere (taken as the depth to the top of the low-velocity zone) cannot be determined with precision.Global average models for shear-wave attenuation (Q^?1_β) have been derived from global average Rayleigh wave attenuation coefficients for periods between 50 and 300 s and average shear-velocity models. Zones of high Q^?1_β coincide with the low-velocity zones of all shear-velocity models, however, models with low-velocity zones beginning at a depth of 70 km have the highest-attenuation layer in the lower half of the low-velocity zone. Resolution kernels for these attenuation models show that parameters for layers shallower than the lower part of the low-velocity-high-attenuation zone are strongly coupled but are distinct from the lower part of this zone. This suggests that the deeper part of the low-velocity-high-attenuation zone is the most mobile part of the zone or that on the average, the top of the zone is deeper than 70 km.The average Q_β of the lithosphere, low-velocity zone, and sub-low-velocity layer (asthenosphere) are approximately 200, 85–110 and 170–200, respectively.     

Crust and upper mantle structure of East Asia from ambient noise and earthquake surface wave tomography

The complex tectonic background of East Asia makes it an ideal region for investigating the evolution of the continental lithosphere, for which high-resolution lithospheric structural models are essential. In this study, we measured Rayleigh-wave phase-velocity dispersion curves at periods of 10–120 s and group velocity dispersion curves at periods of 10–140 s using event records from more than 1,000 seismic stations in and around China. By jointly inverting new and previously published dispersion data from ambient noise and earthquakes, we developed a high-resolution shear-wave velocity model down to a depth of ～300 km beneath East Asia. Our model revealed heterogeneous lithospheric structures beneath East Asia, and allowed us to investigate the velocity structure of the entire lithosphere. We also derived crustal and lithospheric thickness models from the three-dimensional (3D) shear-wave model, revealing strong spatial heterogeneity and a general thinning trend of lithospheric thickness from west to east across the study region. Overall, our models reveal important lithospheric features beneath East Asia and provide a valuable baseline dataset for understanding continental-scale dynamics and evolution.     

深地震测深揭示的华南地区地壳结构及其动力学意义

20世纪70年代以来,我国实施与华南地区有关的深地震测深剖面达57条.本文收集该区的深地震测深研究成果,利用保真能力强的三维克里金插值技术构建了100°E～125°E,18°N～34°N区域内的三维地壳速度模型.基于三维地壳结构模型,分别探讨了华南地区不同构造单元地壳厚度空间变化特征、地壳属性、上地幔顶部地震波速变化以...     

基于密集流动地震台阵的青藏高原东北缘及邻区Rayleigh波相速度层析成像

本文利用喜马拉雅二期台阵674个流动地震台一年的远震垂直向观测资料,获取了18491条独立路径上的Rayleigh波相速度频散曲线,并反演得到了周期10~80 s的Rayleigh波的相速度分布图.通过对比,本文与已有成像结果具有较一致的高低速分布特征,表明了本文结果的可靠性.结果显示,研究区的相速度分布存在明显的横向非均匀性.短周期(如10~15 s)的相速度分布与地表地形密切相关,中等周期(如20~40 s)的相速度分布受地壳厚度的影响较大.在长周期(如60~80 s),鄂尔多斯块体的高速比阿拉善块体更显著、完整,表明同属于华北板块的阿拉善地块,其上地幔结构并没有鄂尔多斯稳定.从短周期至长周期,与周缘地块相比,青藏高原始终表现出较明显的低速异常,可能暗示其具有较活跃的地壳上地幔结构.松潘一甘孜和北祁连块体的中上地壳均存在低速层.全球参考模型Crust1.0和Lithos1.0均不能很好地解释我们的观测频散,基于本文获得的相速度结果可在很大程度上对Crust1.0和Lithos1.0模型进行补充和完善.     

中国大陆及邻区Rayleigh面波相速度分布特征

本文根据102个数字化台站记录的长周期垂直向面波资料,利用双台互相关方法测量了538条独立路径的基阶Rayleigh面波相速度频散资料,反演获得了中国大陆及边邻地区（70°E～140°E,18°N～55°N）20～120 s（周期间隔为5 s） 共21个周期的Rayleigh波相速度空间分布图像. 检测板测试结果显示中国大陆中东部地区横向分辨率可达3°,而西部及边邻地区大约5°. 研究表明,中国大陆地区的Rayleigh波相速度分布横向差异显著,大致以104°E为界,可分成具有不同速度结构特征的东、西两部分. 一般较短周期（20～35 s）的相速度分布受地形和地壳厚度的影响较大,总体表现为东部速度高,西部速度低;塔里木盆地、青藏地块及其东缘的松潘-甘孜地块形成整个研究区内最为突出的低速异常体,蒙古西部低速特征也较清晰;东部的四川盆地、扬子地块、华南地块、松辽盆地、日本海及蒙古东部高速特征明显. 随着周期的增大,青藏地块中部的低速异常体横向尺度逐渐缩小,而喜马拉雅冲断带、塔里木盆地相速度不断升高,意味着青藏低速区受到南、西北、东三个方向的高速区夹击,可能导致高原中部软弱的低速物质向东南方向迁移;同时,东部地区由高速逐渐转变为大面积的低速分布,反映东部地区岩石圈较薄而软流圈发育. 随着青藏地块低速特征的减弱,印支地块北部及相邻海域、东海、东北吉林深震区、日本海、中-朝地块至蒙古东部成为120 s周期上突出的低速异常体,而上扬子地块包括四川盆地高速特征依然明显,显示出稳定的古板块特征. 南北地震带始终呈现出相对较低的速度特征,并成为划分中国大陆具有不同岩石圈相速度特征的东部与西部的天然分界.     

福建-台湾地区一维地壳速度结构的初始模型

文中收集整理了利用人工地震测深、 天然地震到时以及布格重力异常、 地震地质资料反演福建-台湾地区地壳结构的不同结果。 对这些不同的初始模型进行了分析和比较, 给出了该区域的综合一维平均地壳速度结构模型, 为进一步开展该区域二维或三维地壳精细结构反演提供一个合适的一维初始地壳模型。     

The statistical laws of strong seismic motion phase characteristics

ＴｈｅｓｔａｔｉｓｔｉｃａｌｌａｗｓｏｆｓｔｒｏｎｇｓｅｉｓｍｉｃｍｏｔｉｏｎｐｈａｓｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓＸｉｎｇＪＩＮ（金星）ａｎｄＺｈｅｎ－ＰｅｎｇＬＩＡＯ（廖振鹏）（ＩｎｓｔｉｔｕｔｅｏｆＥｎｇｉｎｅｅｒｉｎｇＭｅｃｈａｎｉｃｓ，Ｓ...