2. 中国科学院地质与地球物理研究所, 岩石圈演化国家重点实验室, 北京 100029
2. State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
核幔边界D″层是指从核幔边界(Core-Mantle Boundary, CMB)到其上方250 km左右范围内的下地幔底部区域(Bullen, 1949),是一个强烈的热-化学边界层.在全地幔对流模型中这一区域可能是地幔柱起源以及板块俯冲终结的位置,其热结构对地幔以及外核对流都有十分重要的影响(Lay and Garnero, 2004);同时下地幔底部一些区域可能堆积了与上部地幔化学组分存在显著差异的物质,反映了地幔物质和地球内部热状态的演化过程(Labrosse et al., 2007; Nomura et al., 2011).
地震学方法是研究地球内部结构的强有力手段,已经开展的地震学研究揭示了核幔边界D″层的多尺度复杂结构.全球层析成像结果显示位于非洲和太平洋下方的核幔边界区域分别存在大尺度的剪切波低速异常区,而在环太平洋区域则存在高速异常区(Grand, 2002; Panning and Romanowicz, 2006).波形拟合研究进一步显示非洲地区下地幔底部的剪切波低速异常扰动约为-2%~-12%(Wang and Wen, 2007a);太平洋地区下地幔底部的剪切波低速异常扰动约为-3%~-5%,并且在低速区局部边缘存在小尺度的超低速区(Ultra Low Velocity Zone, ULVZ),剪切波速度扰动可达-10%(He and Wen, 2009).对于环太平洋剪切波高速异常区,波形拟合研究显示局部区域D″层上部存在速度间断面,其剪切波速度扰动约为+2%~+3%(Lay and Helmberger, 1983).
除了剪切波速度的不均匀分布,矿物物理实验和地震学研究揭示D″层还具有地震各向异性特征,显著区别于下地幔的各向同性性质(Meade et al., 1995),表明D″层具有产生地震各向异性的独特物质结构/化学成分以及地球动力学条件.早期地震学研究中发现核幔边界反射震相ScS的径向分量和切向分量存在明显的走时差.Cormier(1986)提出D″层的具有垂直对称轴的横向各向同性模型(Vertical Transversely Isotropy, VTI)来解释这一现象,这成为后来学者利用走时分析方法研究D″层地震各向异性的基本假设.对核幔边界衍射震相Sdiff的波形拟合研究表明太平洋东部局部区域D″层可能存在方位各向异性(Maupin, 1994),Garnero等(2004)则进一步用D″层的具有倾斜对称轴的横向各向同性模型(Tilted Transverse Isotropy, TTI)来解释穿过加勒比海地区下方核幔边界区域的Sdiff波形异常.非洲下方下地幔底部低速异常区边界存在复杂的方位各向异性,可能与地幔流动状态和强烈的剪切变形作用有关(Wang and Wen, 2007b; Lynner and Long, 2014; Ford et al., 2015).剪切波分裂研究显示,太平洋北部和东部以及西伯利亚地区的D″层剪切波高速异常区域均存在TTI结构的各向异性,可能与俯冲板片物质及地幔流动状态有关(Wookey et al., 2005; Wookey and Kendall, 2008; Long, 2009).对更多的D″层剪切波高速异常区域的地震各向异性进行研究,将会有助于我们进一步理解各向异性的来源,认清各向异性的形成机制.全球三维径向各向异性模型显示东南亚下方核幔边界D″层是剪切波高速异常区,并且横向偏振的SH波速要大于垂直偏振的SV波速(Panning and Romanowicz, 2006).前人的研究结果进一步表明这个区域的剪切波高速异常区被一个非异常区分为东西两部分,推测和不同的古俯冲板块相关(He and Wen, 2011).对这个区域核幔边界D″层的地震各向异性分布特征的研究有助于进一步理解这一区域的地球动力学演化过程,具有重要的科学意义.
剪切波分裂参数分析方法是获取各向异性信息的有效手段.国内外学者通过分析多种震相资料的分裂参数可以约束地球内部地壳和上地幔各向异性,深刻反映了过去或者现今大规模的构造运动和上地幔变形特征(Chen et al., 2005; Li and Chen, 2006; Levin et al., 2008; 江丽君等, 2010; Li et al., 2010, 2011; Huang et al., 2011; Nagaya et al., 2011; Zhao et al., 2011; 冯强强等,2012; Chen et al., 2015; 于勇等,2016).由于地震各向异性广泛分布于地壳和上地幔区域,因此穿过D″层震相的剪切波分裂参数通常无法直接反映D″层的地震各向异性.但通过选择特定震相对进行剪切波分裂参数分析,消除地壳和上地幔各向异性对单个震相的影响,也能够有效地获得D″层的各向异性信息(Wookey et al., 2005; Wookey and Kendall, 2008; Long, 2009).例如直达S和ScS这对震相(图 1)的剪切波分裂参数,可以对核幔边界D″各向异性进行很好的约束.直达S波携带的是震源一侧上地幔以及接收台站下方地壳和上地幔的各向异性信息;而ScS波是核幔边界的反射波,在一定震中距范围内与直达S波在地壳和上地幔区域的射线路径非常相近,不同之处在于其经过核幔边界D″层,从而携带有D″层的各向异性信息.因此直达S和ScS的剪切波分裂参数若存在差异,则可能反映的是D″层地震各向异性.
在本研究中,我们使用来自深震的S和ScS震相对,通过把S震相的分裂参数作为上地幔各向异性校正量,对相应的ScS波形数据进行上地幔各向异性校正,获得了东南亚下方下地幔底部的各向异性的异常分布.结合前人速度结构的研究结果,我们推测这里的地震各向异性变化可能与古俯冲板块相关.
2 数据资料与处理结果本研究收集了中国国家数字地震台网宽频带台站记录的2001年至2015年,震中距55°~85°,震级大于5.8的地震资料(郑秀芬等,2009),同时利用国际数字地震台网FDSN的台站来增加研究区域的数据覆盖(图 2).这些台站主要记录了西太平洋以及东南亚地区的地震事件,ScS在CMB的反射点对于中国大陆和东南亚下方的核幔边界D″层有较好的覆盖.
对采集的地震数据去中值和趋势,同时考虑到S和ScS震相的主频范围约为0.05~0.1 Hz,对所有数据进行0.01~0.3 Hz的双通道双极点Butterworth带通滤波以压制噪声,提高信噪比,然后挑选S和ScS波形清晰的高信噪比地震记录进行后续数据处理.为了减少震源一侧上地幔区域各向异性的影响,选择震源深度大于400 km的事件,同时深震的数据具有较高的信噪比,保证随后剪切波分裂参数测量结果的可靠性.实际本研究开始时,并未对研究所使用事件的震源深度做出限制,因为S和ScS差异分裂研究方法同样可以应用于浅震,结合SK(K)S分裂结果估计震源一侧的各向异性,用于校正ScS震相波形.但随着研究展开,高质量的S和ScS震相对分裂参数测量结果均来自深震数据,因此这些测量结果不包含震源一侧地震各向异性信息.
在进行剪切波分裂参数分析之前要保证台站仪器三分量极性的可靠性,因此我们采用Niu和Li (2011)的研究结果对定向存在偏差的部分CDSN台站进行极性校正,然后利用SplitLab (Wüstefeld et al., 2008)程序,同时采用旋转互相关法(Fukao, 1984)和特征值法(Silver and Chan, 1991)对S和ScS震相进行分裂参数分析(图 3).为了确保测量结果的可靠性,我们只保留旋转互相关法和特征值法测量结果一致的分裂参数,并进一步根据各向异性校正后的质点运动轨迹筛选线性偏振程度高的结果,最终获得了24对高质量的S和ScS分裂参数,24对数据的地理分布如图 2b所示.表 1给出了24对S和ScS震相的事件列表,图 3给出了分裂参数测量的一个实例.所有结果及相关信息将汇总在表 2中.
至此本研究获得了24对S和ScS的各向异性分裂参数,同一对S和ScS分裂参数并不完全一致,有些存在明显差异(特别是分裂时间).这种差异可能包含了来自核幔边界D″层的地震各向异性信息,需要做进一步分析从数据中提取出来.对核幔边界D″层地震各向异性研究需要对上地幔各向异性有较好的认识,因为所有能够反映D″层地震各向异性的震相都不可避免的受到上地幔各向异性的影响,包括此处的ScS震相.因此需要对ScS进行上地幔各向异性校正,得到ScS剩余分裂参数做进一步分析.
2.1 上地幔各向异性前人在台站所在区域利用SKS、SKKS和深震S等多种震相资料进行剪切波分裂参数观测研究,取得了上地幔地震各向异性的可靠结果(Chen et al., 2005; Huang et al, 2011; 冯强强等, 2012).将本研究的深震S震相分裂参数ΓS与前人上地幔各向异性研究结果ΓUM进行对比,发现两者具有高度的一致性(表 2),表明本研究中深震S波的分裂参数主要反映了台站下方的上地幔各向异性.
以深震S震相的分裂参数作为本研究的上地幔各向异性校正量,对波形数据进行上地幔各向异性校正,结果显示校正后深震S震相的分裂参数测量结果均为空值(Null),即无分裂状态,并且质点运动轨迹呈现很好的线性特征;同时我们采用Harvard Global CMT的震源机制解(www.globalcmt.org),利用位错点源模型计算远场S波的理论偏振方向,并与校正后的深震S波质点运动进行对比,发现质点运动与理论偏振方向两者吻合得很好(图 4).以上结果表明了上地幔各向异性校正的可靠性,同时避免在校正过程中人为引入新的各向异性.
对地震数据进行上地幔各向异性校正后,S震相各向异性效应已被消除,质点振动呈线性(图 4中Ⅱ列);对ScS震相质点振动分析结果显示,19个波形数据的ScS震相质点呈现出显著的椭圆运动轨迹(如图 4中的ZSU和HYS),表明这些数据中仍然存在显著的ScS剩余各向异性,进一步对ScS剩余各向异性进行测量,分裂参数ΓScS-S见表 2;其余5个ScS震相质点呈现较为明显的线性运动轨迹(如图 4中的HEF和JHE),表明上地幔各向异性校正后的ScS剩余各向异性为Null.
由于下地幔区域除了核幔边界D″层有显著的地震各向异性,其余部分主要表现为各向同性(Meade et al., 1995);同一对ScS与S的差异分裂结果,即ScS剩余各向异性分裂参数,反映的是ScS在核幔边界区域路径上的D″层各向异性信息.为了便于分析分裂参数与核幔边界反射点以及射线路径的几何关系,我们采用将ScS剩余各向异性分裂参数从地理坐标系变换到射线坐标系中:(Φ, δt)→(Φ*, δt), Φ*=Baz-Φ,Φ*表示在ScS核幔边界反射点位置沿着近水平传播的射线方向,观察到的快轴方向与垂直方向的夹角(Wookey and Kendall, 2008).图 5展示了本研究最终获得的射线坐标系下ScS剩余分裂参数的分布情况.
本文利用CDSN和FDSN台网地震数据,采用S和ScS差异分裂方法对东南亚下方的核幔边界D″层各向异性进行研究,最终得到了24个射线坐标系下高质量的ScS剩余各向异性分裂参数(图 5),可能反映了研究区域核幔边界D″层各向异性分布信息.从图 5中红色线段的分布可以看出,19个ScS剩余分裂参数之间差异不大,平均分裂时间δt=1.08±0.35 s,平均快轴偏振方向Φ*=71.4°±30.2°,比较接近SH水平方向(Φ*=±90°),这与Thomas等(2007)的研究结果比较一致.因此我们推测研究区域内核幔边界D″层各向异性结构整体上为对称面接近水平的TTI结构,与VTI结构差异不大,具有VSH>VSV的特征,这也符合全球径向各向异性层析成像模型(Panning and Romanowicz, 2006)的结果,即核幔边界剪切波高速区具有一阶VTI对称结构.
另外5个ScS剩余分裂参数结果为Null,其中4个分布于在经度107°~112°范围内,表明该核幔边界区域为各向同性,或者ScS入射波偏振方向与各向异性对称轴垂直或平行.但以表 2中事件3分裂结果为例,可以看到SCH,WSU和JHE三个台站的剩余分裂为Null,FUY和HYS台站的剩余分裂观测值分别为(28°, 0.9 s)和(41°, 1.1 s),而这5个台站的入射波偏振方向几乎是一致的,表明该经度范围内的核幔边界各向异性显著区别于两侧区域.同时该区域的另一个Null分裂来自事件5的HEF台站,其入射波偏振方向又与事件3中各台站不同.因此我们推断经度107°~112°范围内的Null分裂应该不是由于特定的入射波偏振方向引起的,其结果最可能代表的是局部区域的各向同性.
4 结论与讨论综合上述分析结果,24个ScS剩余分裂参数的分布具有区域特征:在经度107°~112°范围内,剩余各向异性分裂参数均为Null(图 5黄色圆点所示);而在该经度范围东西两侧区域,各向异性分裂结果基本均为非Null(西侧一个Null结果除外).剩余分裂参数区域分布特征表明研究区域核幔边界D″层各向异性存在区域变化.He和Wen(2011)对欧亚大陆东部下方核幔边界区域开展的S-ScS走时残差分析和波形拟合研究揭示,在我们研究区域附近的核幔边界D″层存在具有不同速度结构和不同间断面性质的东西两块高速区域,可能与古特提斯和古太平洋俯冲板块相关(van der Meer et al., 2010; He and Wen, 2011).McNamara等(2002)的地球动力学模拟结果显示,核幔边界区域俯冲板块存在强烈的剪切应变.矿物物理实验已经发现在D″层温压条件下可能存在钙钛矿(pv)到过钙钛矿(ppv)的相变,导致速度结构的异常.在强烈的剪切应变条件下,过钙钛矿的晶格优势排列(LPO,Lattice Preferred Orientation)则可能是造成剪切波高速区D″层存在显著各向异性的微观机制(Murakami et al., 2004; Nowacki et al., 2013).结合前人研究结果,我们推测本研究中各向异性的区域变化可能同样与古特提斯和/或古太平洋俯冲板块相关(图 6).
前人研究结果表明约束TTI对称轴的几何参数需要采用多路径多方位角分布的各向异性分裂参数(Wookey et al., 2008; Nowacki et al., 2010).受地震事件和台站分布限制,我们获取的高质量分裂参数均来自同一后方位角分布(125°~144°)的地震数据,因此无法利用获得的剩余分裂参数对该区域核幔边界D″层TTI对称结构做出进一步的约束.另外本研究中高质量的各向异性参数数量仍然有限,目前还无法约束不同各向异性区域的分布范围.进一步的研究需要结合其他震相对的研究结果,如SKS和SKKS差异分裂(Niu and Perez, 2004; Long, 2009; He and Long, 2011;Ford et al., 2015)等,来共同约束可能的各向异性对称结构,提高对各向异性精细结构的分辨能力,从而进一步探讨核幔边界D″层各向异性产生的机制及相应的地球动力学过程.
致谢感谢中国地震局地球物理研究所国家数字测震台网数据备份中心为本研究提供CDSN台网地震波形数据,感谢IRIS为本研究提供FDSN台网的地震波形数据.研究中使用Wüstefeld的SplitLab软件对S和ScS进行剪切波分裂参数的分析.
Bullen K E. 1949. Compressibility-pressure hypothesis and the Earth's interior. Geophysical Journal International , 5 (9) : 355-368. | |
Chen Y P, Wang L S, Mi N, et al. 2005. Shear wave splitting observations in the Chinese Tianshan orogenic belt. Geophysical Research Letters , 32 (7) : 303-341. | |
Chen Y, Li W, Yuan X H, et al. 2015. Tearing of the Indian lithospheric slab beneath southern Tibet revealed by SKS-wave splitting measurements. Earth & Planetary Science Letters , 413 : 13-24. | |
Cormier V F. 1986. Synthesis of body waves in transversely isotropic Earth models. Bulletin of the Seismological Society of America , 76 (1) : 231-240. | |
Dziewonski A M, Anderson D L. 1981. Preliminary reference Earth model. Physics of the Earth & Planetary Interiors , 25 (4) : 297-356. | |
Feng Q Q, Wu Q J, Li Y H, et al. 2012. Shear wave splitting in Xinjiang region. Acta Seismologica Sinica , 34 (3) : 296-307. | |
Ford H A, Long M D, He X B, et al. 2015. Lowermost mantle flow at the eastern edge of the African Large Low Shear Velocity Province. Earth & Planetary Science Letters , 420 : 12-22. | |
Fukao Y. 1984. Evidence from core-reflected shear waves for anisotropy in the Earth's mantle. Nature , 309 (5970) : 695-698. DOI:10.1038/309695a0 | |
Garnero E J, Maupin V, Lay T, et al. 2004. Variable azimuthal anisotropy in Earth's lowermost mantle. Science , 306 (5694) : 259-261. DOI:10.1126/science.1103411 | |
Grand S P. 2002. Mantle shear-wave tomography and the fate of subducted slabs. Philosophical Transactions of the Royal Society A:Mathematical, Physical & Engineering Sciences , 360 (1800) : 2475-2491. | |
He X B, Long M D. 2011. Lowermost mantle anisotropy beneath the northwestern Pacific:Evidence from PcS, ScS, SKS, and SKKS phases. Geochemistry, Geophysics, Geosystems , 12 (12) . DOI:10.1029/2011GC003779 | |
He Y M, Wen L X. 2009. Structural features and shear-velocity structure of the "Pacific Anomaly". Journal of Geophysical Research , 114 (B2) . DOI:10.1029/2008JB005814 | |
He Y M, Wen L X. 2011. Seismic velocity structures and detailed features of the D" discontinuity near the core-mantle boundary beneath eastern Eurasia. Physics of the Earth & Planetary Interiors , 189 (3-4) : 176-184. | |
Huang Z C, Wang L S, Zhao D P, et al. 2011. Seismic anisotropy and mantle dynamics beneath China. Earth & Planetary Science Letters , 306 (1-2) : 105-117. | |
Jiang L J, Li Y H, Wu Q J. 2010. The shear wave splitting of Central Tien Shan and its implications. Chinese J. Geophys. , 53 (6) : 1399-1408. DOI:10.3969/j.issn.0001-5733.2010.06.018 | |
Labrosse S, Hernlund J W, Coltice N. 2007. A crystallizing dense magma ocean at the base of the Earth's mantle. Nature , 450 (7171) : 866-869. DOI:10.1038/nature06355 | |
Lay T, Helmberger D V. 1983. A lower mantle S-wave triplication and the shear velocity structure of D". Geophysical Journal International , 75 (3) : 799-837. DOI:10.1111/j.1365-246X.1983.tb05010.x | |
Lay T, Garnero E J. Core-mantle boundary structures and processes.//Sparks R S J, Hawkesworth C J eds. The State of the Planet:Frontiers and Challenges in Geophysics. Washington, DC: AGU, 2004 . | |
Levin V, Roecker S, Graham P, et al. 2008. Seismic anisotropy indicators in Western Tibet:Shear wave splitting and receiver function analysis. Tectonophysics , 462 (1-4) : 99-108. DOI:10.1016/j.tecto.2008.03.019 | |
Li A B, Chen C Z. 2006. Shear wave splitting beneath the central Tien Shan and tectonic implications. Geophysical Research Letters , 33 (22) : 217-234. | |
Li Y H, Wu Q J, Jiang L J, et al. 2010. Complex seismic anisotropic structure beneath the central Tien Shan revealed by shear wave splitting analyses. Geophysical Journal International , 181 (3) : 1678-1686. | |
Li Y H, Wu Q J, Zhang F X, et al. 2011. Seismic anisotropy of the Northeastern Tibetan Plateau from shear wave splitting analysis. Earth & Planetary Science Letters , 304 (1-2) : 147-157. | |
Long M D. 2009. Complex anisotropy in D" beneath the eastern Pacific from SKS-SKKS splitting discrepancies. Earth & Planetary Science Letters , 283 (1-4) : 181-189. | |
Lynner C, Long M D. 2014. Lowermost mantle anisotropy and deformation along the boundary of the African LLSVP. Geophysical Research Letters , 41 (10) : 3447-3454. DOI:10.1002/2014GL059875 | |
Maupin V. 1994. On the possibility of anisotropy in the D" layer as inferred from the polarization of diffracted S-waves. Physics of the Earth & Planetary Interiors , 87 (1-2) : 1-32. | |
McNamara A K, van Keken P E, Karato S I. 2002. Development of anisotropic structure in the Earth's lower mantle by solid-state convection. Nature , 416 (6878) : 310-314. DOI:10.1038/416310a | |
Meade C P, Silver P G, Kaneshima S. 1995. Laboratory and seismological observations of lower mantle isotropy. Geophysical Research Letters , 22 (10) : 1293-1296. DOI:10.1029/95GL01091 | |
Murakami M, Hirose K, Kawamura K, et al. 2004. Post-perovskite phase transition in MgSiO3. Science , 304 (5672) : 855-858. DOI:10.1126/science.1095932 | |
Nagaya M, Oda H, Kamimoto T. 2011. Regional variation in shear-wave polarization anisotropy of the crust in southwest Japan as estimated by splitting analysis of Ps-converted waves on receiver functions. Physics of the Earth & Planetary Interiors , 187 (1-2) : 56-65. | |
Niu F L, Perez A M. 2004. Seismic anisotropy in the lower mantle:A comparison of waveform splitting of SKS and SKKS. Geophysical Research Letter , 31 (24) : 357-370. | |
Niu F L, Li J. 2011. Component azimuths of the CEArray stations estimated from P-wave particle motion. Earthquake Science , 24 (1) : 3-13. DOI:10.1007/s11589-011-0764-8 | |
Nomura R, Ozawa H, Tateno S, et al. 2011. Spin crossover and iron-rich silicate melt in the Earth's deep mantle. Nature , 473 (7346) : 199-203. DOI:10.1038/nature09940 | |
Nowacki A, Wookey J, Kendall J M. 2010. Deformation of the lowermost mantle from seismic anisotropy. Nature , 467 (7319) : 1091-1094. DOI:10.1038/nature09507 | |
Nowacki A, Walker A M, Wookey J, et al. 2013. Evaluating post-perovskite as a cause of D" anisotropy in regions of palaeosubduction. Geophysical Journal International , 192 (3) : 1085-1090. DOI:10.1093/gji/ggs068 | |
Panning M, Romanowicz B. 2006. A three-dimensional radially anisotropic model of shear velocity in the whole mantle. Geophysical Journal International , 167 (1) : 361-379. DOI:10.1111/gji.2006.167.issue-1 | |
Silver P G, Chan W W. 1991. Shear wave splitting and subcontinental mantle deformation. Journal of Geophysical Research:Solid Earth , 96 (B10) : 16429-16454. DOI:10.1029/91JB00899 | |
Simmons N A, Forte A M, Boschi L, et al. 2010. GyPSuM:A joint tomographic model of mantle density and seismic wave speeds. Journal of Geophysical Research , 115 (B12) : 424-440. | |
Thomas C, Wookey J, Simpson M. 2007. D" anisotropy beneath Southeast Asia. Geophysical Research Letters , 34 (4) . DOI:10.1029/2006GL028965 | |
van der Meer D G, Spakman W, van Hinsbergen D, et al. 2010. Towards absolute plate motions constrained by lower-mantle slab remnants. Nature Geoscience , 3 (1) : 36-40. DOI:10.1038/ngeo708 | |
Wang Y, Wen L X. 2007a. Geometry and P and S velocity structure of the "African Anomaly". Journal of Geophysical Research , 112 (B5) . DOI:10.1029/2006JB004483 | |
Wang Y, Wen L X. 2007b. Complex seismic anisotropy at the border of a very low velocity province at the base of the Earth's mantle. Journal of Geophysical Research , 112 (B9) . DOI:10.1029/2006JB004719 | |
Wookey J, Kendall J M. 2008. Constraints on lowermost mantle mineralogy and fabric beneath Siberia from seismic anisotropy. Earth & Planetary Science Letters , 275 : 32-42. | |
Wookey J, Kendall J M, Rumpker G. 2005. Lowermost mantle anisotropy beneath the north Pacific from differential S-ScS splitting. Geophysical Journal International , 161 (3) : 829-838. DOI:10.1111/gji.2005.161.issue-3 | |
Wüstefeld A, Bokelmann G, Zaroli C, et al. 2008. SplitLab:A shear-wave splitting environment in Matlab. Computers & Geosciences , 34 (5) : 515-528. | |
Yu Y, Chen Y S, Jian H C, et al. 2016. SKS wave splitting study of the transition zone at the central portion of the North China Craton. Chinese J. Geophys. , 59 (1) : 141-151. DOI:10.6038/cjg20160111 | |
Zhao L, Zheng T Y, Lü G, et al. 2011. No direct correlation of mantle flow beneath the North China Craton to the India-Eurasia collision:Constraints from new SKS wave splitting measurements. Geophysical Journal International , 187 (2) : 1027-1037. DOI:10.1111/gji.2011.187.issue-2 | |
Zheng X F, Ouyang B, Zhang D N, et al. 2009. Technical system construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake. Chinese J. Geophys. , 52 (5) : 1412-1417. DOI:10.3969/j.issn.0001-5733.2009.05.031 | |
冯强强, 吴庆举, 李永华, 等. 2012. 新疆地区S波分裂研究. 地震学报 , 34 (3) : 296–307. | |
江丽君, 李永华, 吴庆举. 2010. 中天山及邻区S波分裂研究及其动力学意义. 地球物理学报 , 53 (6) : 1399–1408. DOI:10.3969/j.issn.0001-5733.2010.06.018 | |
于勇, 陈永顺, 菅汉超, 等. 2016. 华北克拉通中部过渡带SKS波分裂研究:鄂尔多斯东南角的局部软流圈绕流. 地球物理学报 , 59 (1) : 141–151. DOI:10.6038/cjg20160111 | |
郑秀芬, 欧阳飚, 张东宁, 等. 2009. "国家数字测震台网数据备份中心"技术系统建设及其对汶川大地震研究的数据支撑. 地球物理学报 , 52 (5) : 1412–1417. DOI:10.3969/j.issn.0001-5733.2009.05.031 | |