2. 成都理工大学地球探测与信息技术教育部重点实验室, 成都 610059;
3. 成都理工大学地质灾害防治与地质环境保护国家重点实验室, 成都 610059;
4. 国土资源部构造成矿成藏重点实验室, 成都 610059;
5. 四川省地震局, 成都 610041
2. Key Lab. of Earth Exploration and Information Technique of Ministry of Education(Chengdu University of Technology), Chengdu 610059, China;
3. State Key Lab. of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology), Chengdu 610059, China;
4. Key Laboratory of Tectonic Controls on Mineralization and Hydrocarbon Accumulation, Ministry of Land and Resources(Chengdu University of Technology), Chengdu 610059, China;
5. Sichuan Earthquake Agency, Chengdu 610041, China
青藏高原是全球造山带研究的热点地区,关于青藏高原的隆升及变形机制,已成为研究的重要主题.近年来,国内外学者对该问题的讨论提出了许多模型(Tapponnier et al., 1982;Houseman and England, 1993;Royden et al., 1997;Xu et al., 1999;Clark and Royden, 2000;Beaumont et al., 2001;Burchfiel et al., 2008;Liu et al., 2014;Zhao et al., 2013),包括上地壳的纯剪切增厚和中下地壳的黏性物质管道流.前者认为脆性上地壳发生增厚,主要为上地壳的断层和褶皱作用,而黏性的中下地壳发生缩短(Hubbard and Shaw, 2009);后者认为黏性中下地壳使上地壳从上地幔解耦,地表不存在显著的水平缩短(Clark et al., 2005).在这个争论中,青藏高原东缘被认为是高原内部物质向东运移途径的地区之一.因此,青藏高原东缘的地壳结构,对认识青藏高原的隆升变形和动力学演化具有重要意义.
此外,青藏高原也是中国大陆强震频繁发生的地区,是中国大陆7级地震活动的主要区域(邓起东等, 2014).自1995年以来,青藏高原已进入又一次地震丛集高潮期(邓起东等, 2014).在青藏高原东缘,一系列规模不等、力学性质不同的活动断裂带控制着该区地震的不均匀空间分布(徐锡伟等, 2003),近年来即已密集发生了2008年汶川M8.0地震(朱介寿, 2008;易桂喜等, 2012)、2013年芦山M7.0地震(杨宜海等, 2015;易桂喜等, 2016)和2017年九寨沟M7.0地震(杨宜海等, 2017;易桂喜等, 2017a).因此,研究该地区的地壳精细结构,可以帮助提升对该地区地震成因机制的认识.
背景噪声层析成像摆脱了传统成像方法对地震事件的依赖性,利用相同时间段不同台站记录的连续背景噪声,采用互相关获取台站间的经验格林函数(Weaver and Lobkis, 2001;Wapenaar, 2004;Roux et al., 2005;Yao et al., 2006;Bensen et al., 2007),目前已有多个研究组在青藏高原东缘及邻区开展噪声成像研究(如Yao et al., 2008, 2010;Li et al., 2009;Yang et al., 2012;Liu et al., 2014).然而,已有的研究大部分采用线性共轭梯度法进行三维速度反演(Mitchell and Herrmann, 1979),该方法的主要优点是反演速度快,但对初始模型的依赖较大,当反演区域仅使用单一的初始模型时,区域内各点的值均向初始模型的局部极优解收敛,同时会导致反演的三维模型平滑程度过大.现有研究表明,青藏高原东缘的深部结构存在显著的横向变化(如Liu et al., 2014;王苏等,2015;郑晨等,2016),因此很难给出一个单一合理的初始模型来约束反演.此外,该方法一般仅用于反演S波速度变化,很少用来反演其他物性特征(如波速比、地壳各向异性),而这些特征在综合解释成因机制及构造变形方式中同样具有重要作用.
因此,本文采用Liang和Langston(2009)应用的模拟退火(Simulated Annealing)程序进行三维结构反演.该方法对初始模型的要求较小,在给定的区域模型中搜寻全局最优解,可以更真实的反演研究区横向速度变化的同时反演波速比(Vp/Vs)等物性特征,结合其他地球物理资料,共同对青藏高原东缘的侧向生长机制及区域地震成因进行探讨.
1 数据与方法 1.1 数据本研究共搜集了127个宽频带固定地震台站记录的连续一年的波形数据,分别来自四川地震台网、云南地震台网全部宽频带固定地震台站及部分中国地震台网宽频带固定地震台站(图 1).在提取互相关函数前,对垂直分量波形数据以天为单位进行预处理,包括降采样到0.2 s、去均势、去平均、去仪器响应以及带通滤波(0.01~0.5 Hz),最后对预处理数据进行挑选,去除数据有大量缺失的天数据.
通常,背景噪声层析成像研究有四个主要步骤:
(1) 通过互相关计算提取两台站间的互相关函数作为经验格林函数(Bensen et al., 2007;Weaver and Lobkis, 2001).在时间域中,以天为单位将预处理数据进行台站对的互相关计算,去掉信噪比较低的互相关函数后,将以天为单位的互相关函数叠加,得到每个台站对一年期的经验格林函数,获得更为清晰的面波信号.图 2给出与MAL台间距在100~800 km的台站之间的互相关函数,按台间距排列后,可以看出MAL台来自各方位的互相关函数显示出均匀且明显的位移信号.
(2) 利用FTAN(Frequency Time Analysis)方法从经验格林函数中提取群速度频散曲线(Dziewonski et al., 1969;Herrmann, 1973).很多学者基于此方法原理开发了不同程序,本研究所用的程序来自Liang和Langston(2008),该程序加入了相位均衡和相位拟合技术来提高信噪比,同时设置理论时频曲线用于对比真实速度值,确保群速度提取的准确性.理论时频曲线(图 3a)参考了AK135全球模型及杨宜海等(2015)利用接收函数方法建立的区域速度模型.通过人工选取,用于群速度成像及三维S波速度反演的频散曲线数量分布如图 3b所示,最大数量频散曲线主要分布在10 s到30 s周期范围内.
(3) 群速度成像的获取本文参考Barmin等(2001),将研究区划分成0.5°×0.5°的网格,频散曲线的速度在每个网格中的分布基于最小二乘法计算,得到每个网格点确定的群速度值,不同周期下,每个网格点的群速度经过插值计算便可组成群速度结构分布.为了更好的检测反演的可靠性,同时利用该方法对所选频散曲线的分辨率和稳定性进行了测试.
(4) 通过提取不同周期下各网格点的群速度,可以得到该网格点独立的频散曲线,进一步反演频散曲线可以获取相应的一维S波速度,而所有网格点一维S波速度即构成了三维S波速度结构.由于各网格点的频散曲线与一维速度模型是非线性关系,所以利用频散曲线反演三维S波速度一般采用两种方法,一种是线性共轭梯度法(Herrmann and Ammon, 2004),另一种为非线性类反演方法(Panza, 1981; Zhang et al., 2014).前者定义的分层模型一般较为精细并假设相邻层的速度变化比较连续,且容易以局部极值代替最大值,因此要获取较准确合理的结果,对初始模型的要求较高;而非线性类反演法中的源于统计热力物理学的模拟退火(SA)法(Kirkpatrick et al., 1983; Rothman, 1985)则较好的规避了线性反演的弊端.SA是基于统计力学的模拟热平衡退火过程的一种Monte Carlo法(姚姚,1995).该方法在不断的统计实验中进行组合优化,在解空间中寻找拟合度最佳的模型.Rothman(1985)首次将该方法引入到地震勘探资料中,以解决剩余静校正问题.随后,Sen和Stoffa(1991)再次将其应用于一维地震波形反演和一维电测深反演中.自此,产生了大量基于模拟退火法的地震资料研究,Shapiro和Ritzwoller(2002)在对全球壳幔S波速度结构研究中,将此方法改进并应用于面波频散曲线反演.本文延用Liang和Langston(2009)介绍的基于Boltzmann采样技术(Ingber, 1989)的模拟退火法反演程序对三维S波速度和Vp/Vs进行反演.研究中,初始模型相邻层的速度变化可以是任意范围,能够较好体现速度跃变,降低了最优解模型对初始模型的依赖性,同时,除了S波速度的反演,该方法还可设定其他物性参数(如:地壳各向异性、密度、Vp/Vs等)进行共同反演.
2 结果 2.1 检测板测试结果检测板测试利用的射线路径和反演方法与真实群速度成像相同.将研究区划分为0.5°×0.5°的网格,每个网格采用±0.3 km·s-1的速度扰动,反演出较为清晰的异常速度变化模式的网格大小可达1.0°×1.0°.图 4展示了30 s的射线分布图和检测板初始棋盘格模型以及10 s、20 s、30 s和45 s的反演结果.从结果可以看出本文选取的射线可清晰的反演出研究区域的速度异常结构.
为了测试检测板模型和真实模型之间的差异,本文参考Zelt(1998),采用分辨系数R作为衡量参数,R的变化范围为0~1,R值越大,表明检测板模型与真实模型越接近.Zelt(1998)定义R>0.7的区域为可分辨区域.图 4中10 s、20 s、30 s和45 s检测板反演结果图像中用白色曲线圈出的范围为R>0.7的区域,也就是分辨率达到1.0°的区域.可见研究区绝大部分区域反演结果的R值都大于0.7,说明检测板测试结果可以很好地代表真实反演,且本文结果具有较高的分辨率和可靠性.
2.2 群速度反演结果图 5选取了不同周期的群速度反演结果.周期为10 s的结果代表了近地表的面波群速度分布.10 s周期图显示,青藏高原东缘近地表的平均速度特征为高速异常,其速度值基本在3.0 km·s-1以上;而四川盆地(SCB)由于受近地表沉积层的影响(Xu et al., 2008),表现出明显的低速特征(< 2.8 km·s-1),这与该地区的接收函数及面波成像研究结果基本一致(如Zhu et al., 2012; 范莉苹等, 2015;朱介寿等, 2017).
15~35 s范围的群速度主要对10~40 km深度的S波速度较为敏感,而36~45 s群速度可以反映深至60 km的S波速度特征.青藏高原东缘在30 s的群速度显示出大范围的低速异常,与10 s的结果形成明显差异;在40 s低速异常范围有所缩小,主要分布在松潘—甘孜地块(SGB)的西部和川滇地块(CDB)大部分地区,其中川滇地块北部(NCDB)比南部(SCDB)具有更明显的低速特征,而滇南地块(SYB)、滇东地块(EYB)和SCB的群速度则为高速异常.SCB的群速度在10~40 s周期与青藏高原东缘形成明显差异,表现出典型的克拉通构造特征.
2.3 三维S波速度结构及泊松比分布除了S波速度,地壳的泊松比变化也密切的反映了壳内结构特征,在一定程度上与地壳的物质组成及熔融程度有关.胡家富等(2003)及Hu等(2005)在对云南的泊松比分布研究中指出岩石的泊松比变化范围大致在0.20~0.35,铁镁质含量越高泊松比越大,长英质含量越高泊松比越小.同时,泊松比值与壳内二氧化硅含量也紧密相关,与地壳熔融程度成正比,当泊松比大于0.3时地壳可能存在部分熔融(张培震, 2008).因此,研究地壳泊松比为探测地壳岩石类型和组成提供了一个有效方法.
由此,本文根据三维反演得到的Vp/Vs推导出泊松比,从而进行泊松比成像与分析.反演应用的初始模型参数设置如表 1所示,该初始模型参数参考了全球模型AK135以及杨宜海等(2015)通过接收函数反演得到的区域速度模型.对于模型第一层,给定沉积层的速度,反演沉积层的厚度及Vp/Vs,以下各层给定地层厚度,反演每层的速度及Vp/Vs,通过Vp/Vs与泊松比的关系(王兴臣等,2017)进一步推导每一层的泊松比值.
图 6给出了10~40 km的S波速度分布.在上地壳,10 km深度SCB表现为相对低的速度特征(~3.0 km·s-1),而这一深度下青藏高原东缘普遍为高速特征.20 km的低速异常主要分布在SYB和EYB,而SGB和CDB依然显示高速特征.30 km,在青藏高原东缘出现大范围的低速异常,具体而言,SGB的速度特征大致以龙日坝断裂为界,与鲜水河断裂带、龙门山断裂带及岷山隆起所围限区域呈现明显的低速异常,而龙日坝断裂以西为高速异常;CDB基本全部为低速异常;SCB、SYB和EYB均显示明显的高速异常,推测已进入高速的下地壳(>3.6 km·s-1).40 km的低速异常范围明显减小,其中龙日坝断裂以东地区从低速异常变为高速异常,而以西则变为低速;CDB大范围的低速异常也基本消失,S波速度大致在3.8~4.0 km·s-1,仅NCDB仍有少量低速异常.SCB、SYB和EYB的S波速度可达4.0 km·s-1以上,表明这些地区已到下地壳底部或岩石圈上地幔.
图 7对应给出了与S波速度成像相同深度的泊松比分布.结果显示,上地壳高泊松比与地壳S波低速异常分布较为吻合,10 km深度的高泊松比(~0.30)主要分布在SCB边缘,龙门山断裂带至龙日坝断裂以及CDB内部.到20 km高泊松比范围逐渐缩小到主要断裂带附近,如东昆仑断裂带、龙门山断裂带、鲜水河断裂带、丽江断裂等.30 km深度的泊松比分布延续了20 km的泊松比特征,东昆仑断裂带、金沙江断裂带、安宁河断裂带附近都出现了大范围的高泊松比分布;在CDB内,高泊松比值区域同样有所增大,但主要集中在SCDB内,NCDB基本不存在高泊松比,这与该深度下S波速度特征完全相反.40 km的泊松比普遍在0.25左右,只有极少地区具有高泊松比结构,这一特征与该深度下S波速度相一致.
青藏高原是由不同时期的多个地块拼合而成,尤其在新生代发生了显著的地壳增厚和地表抬升(Yin and Harrison, 2000).松潘—甘孜地块与扬子克拉通的相互碰撞导致了龙门山在早新生代即开始抬升,其抬升速率在青藏高原的隆升过程中显著增大(Xu et al., 2008).毫无疑问,青藏高原与扬子克拉通的碰撞,以及随后龙门山的加速抬升和高原在东南方向的逃逸,对青藏高原东缘,主要是SGB和CDB,以及龙门山的地壳结构及变形都有巨大的影响.本文在不同构造单元获取的地壳结构剖面(图 8),也揭示了不同构造单元之间的地壳结构及强度存在着明显差异.
SGB是古老的三叠纪造山带,N-S向宽约1000 km,E-W向长约2000 km,以龙门山造山带为其东部边界,在平面上呈三角形展布(Xu et al., 2008).在Clark和Royden(2000)提出的下地壳流模型中,SGB是地壳流主要途经地区之一,而本文研究区覆盖了SGB的东部地区,因此提供了一个探讨下地壳流模型的窗口.此外,在SGB东部边界,近50年来密集发生了一系列强烈地震,包括1976年松潘—平武双震(M7.2),2008年汶川M8.0地震,2013年芦山M7.0地震及2017年九寨沟M7.0地震,所以根据成像结果,本文还将对该区域的孕震环境进行探讨.
根据三维反演结果,本文选取了2条横跨龙门山断裂带的S波速度及泊松比剖面(图 8中P1和P2).剖面成像显示,P1剖面在SGB的中下地壳存在显著的壳内低速异常,其低速层厚度自西向东发生减薄,穿过构造边界带(TB),截止于龙门山断裂带下方;而泊松比剖面显示SGB中下地壳泊松比值基本在0.25左右,不具有显著的高泊松比特征,与SGB巨厚的壳内低速层特征不吻合;P2剖面在SGB也存在明显的中下地壳(~30 km)低速层,低速层同样没有延伸到龙门山下方,与Liu等(2014)采用噪声与接收函数联合反演结果相吻合.龙门山下方的低速软弱物质只发育在15 km以上的地壳中,而汶川和芦山震源区的地震矩心深度大部分分布在10~20 km,也就是高低速的跃变区,并未有深源地震的发生,因此难以用下地壳流模型解释龙门山的抬升.
徐锡伟等(2008)利用GPS资料研究认为龙门山断裂带到龙日坝断裂带是巴颜喀拉地块SE向运移过程中受到华南地块的强烈阻挡而形成的后展式推覆构造,本文在SGB获取的中下地壳低速异常分布与其后展式推覆构造的范围基本吻合.因此本文推测,SGB中下地壳的低速物质可能是青藏高原与扬子块体长期相互作用产生的塑性低速滑脱层;上地壳脆性物质在板块作用下沿中地壳低速滑脱层顶面发生逆冲增厚,造成龙门山的持续抬升和地形起伏,并在TB形成了应变积累和应力集中;而龙门山断裂带的上地壳低速软弱物质为地壳发生破裂提供了有利条件,从而在某种程度上促进了汶川地震和芦山地震的发生.
P1剖面不仅横跨龙门山造山带,同时也横跨岷山隆起.岷山隆起及以南的龙门山中南段第四纪活动比较活跃,被认为共同构成了川西高原的最新活动边界(邓起东等,1994).本文在TB及岷山隆起一带获取了显著的中地壳(~20 km)低速层及呈“凸起”形态的下地壳高泊松比,与接收函数CCP成像结果(Liu et al., 2017;Wang et al., 2018)基本一致,此前的研究显示该地区具有较高的热流及岩石快速抬升等现象(Kirby and Harkins, 2013;姜光政等, 2016),因此本文推测岷山隆起附近大范围且较厚的壳内高泊松比物质可能是由于岩石圈拆沉导致的地幔热物质上涌所致.
3.2 川滇菱形块体的变形机制探讨已有的研究中已揭示了川滇菱形块体壳内低速层发育的复杂性(Yao et al., 2010;Yang et al., 2012;Chen et al., 2013;Xie et al., 2013;Bao et al., 2015;王苏等,2015).关于该区域低速层的形成机制也有很多不一样的解释:其中一种认为川滇菱形块体的低速层来自青藏高原中部,其中部软弱物质在横向压力梯度差的作用下向东南缘挤出,从而形成川滇地区中下地壳大面积的低速区域(如Li et al., 2009;Liu et al., 2014);另一种解释认为川西高原存在的中下地壳低速层可能为地壳的部分熔融,其高泊松比、低Q值和高导低阻层的物性特征,以及贡嘎山出露的浅色花岗岩都为部分熔融现象提供了有力证据(如Xu et al., 2007;王椿镛等, 2008).为了更清晰地显示川滇菱形块体下方的速度结构,本文给出了2条横跨川滇菱形块体的剖面成像结果.
P3剖面显示了2条不同形态的中下地壳低速层:其西侧的低速层从QTB延伸到NCDB,主要分布在~40 km深的下地壳,而东侧的低速层分布在NCDB东部至SGB约30 km的下地壳,两条低速层在理塘断裂附近交汇增厚.然而,在P4剖面中SCDB已不存在明显的壳内低速层,仅在构造边界及两侧的SYB及EYB存在少量中地壳低速异常,表明SCDB的地壳结构相对稳定,其相对高速的中下地壳部分可能和峨眉山大火成岩省有关.泊松比结果显示P3和P4剖面都仅在中地壳存在一条厚度不足5 km的高泊松比异常带,因此难以支持下地壳流模型.
本文在CDB获取的剖面成像结果与Liu等(2014)利用联合反演得到的成像结果存在差异,但是与Chen等(2013)根据接收函数方法计算的Vp/Vs结果相吻合;前者的结果显示低速层在整个川滇菱形块体内具有很好的延续性,而后者根据获取的Vp/Vs推测SCDB不存在明显的壳内低速层,该地区深部地壳变形较小,主要受到红河断裂、鲜水河—小江断裂等活动边界的控制,所以丽江断裂以北的低速层也就无法达到“层流”的条件.GPS观测也发现SCDB的地表形变速率明显低于NCDB(Shen et al., 2005).此外,通过Lg波能量的高衰减来探明地壳高温熔融区域的研究表明,低QLg值(~120)截止在NCDB,以丽江断裂为界,与SCDB的高QLg值(~285)形成强烈对比,可能预示着NCDB与SCDB不同的壳内结构(Zhao et al., 2013).利用SKS分裂测量获取的地幔各向异性结果表明,26°以南的快波方向发生明显变化,同样说明NCDB和SCDB具有不同的深部变形特征(Lev et al., 2006).Chen等(2013)综合Pms和SKS分裂测量结果认为,NCDB存在壳内解耦,可能为下地壳流变形模式,且岩石圈地幔变形主要受软流圈作用,而SCDB主要受到活动边界断裂的控制.地震震源机制揭示了上地壳构造应力场的特征:Luo等(2016)发现JSJF存在大量的正断型地震,表明CDB的西侧可能主要通过上地壳剪切增厚发生构造变形,同时,该研究还发现在SCDB的构造应力场为低倾角近N-S向挤压作用,而NCDB的构造应力场为近E-W向高倾角挤压作用.此外,在易桂喜等(2017b)对四川理塘M4.9和M5.1地震发震构造分析中同样得到NCDB上地壳近E-W向高倾角挤压作用的构造应力场特征.
因此,综合以上研究,本文认为CDB在软流圈地幔的驱动下,其岩石圈沿NW-SE向发生变形,壳内低速物质在NCDB向理塘地块的浅部汇聚,同时形成了左旋走滑的理塘断裂(4.0±1.0 mm·a-1;徐锡伟等, 2005),吸收了高原在SE方向生长的部分形变.作为NCDB和SCDB边界的丽江断裂再次吸收了剩余的大部分形变量,因此SCDB受高原侧向生长的作用相对较弱;结合前人在该地区的地震各向异性研究(Chen et al., 2013),本文推测SCDB变形主要受活动边界断裂控制.
4 结论本文搜集了青藏高原东缘及邻区布设的127个宽频带固定地震台站记录的连续波形资料,反演获取了研究区的面波群速度、三维S波速度及泊松比分布.结合GPS观测、构造应力场和各向异性等研究结果,对青藏高原东缘的深部结构及变形过程进行了探讨,取得以下主要结论:
(1) SGB的中下地壳低速异常主要分布在龙日坝断裂带、鲜水河断裂带、龙门山断裂带和岷山隆起所围限的区域,而该区域的中下地壳仅具有中等泊松比值,推测SGB中下地壳低速物质可能是青藏高原与扬子块体长期相互作用产生的塑性低速滑脱层;上地壳脆性物质在板块作用下沿中地壳低速滑脱层顶面发生逆冲增厚,造成龙门山的持续抬升和地形起伏,并在TB形成了应变积累和应力集中;而龙门山断裂带的上地壳低速软弱物质为地壳发生破裂提供了有利条件,从而在某种程度上促进了汶川地震和芦山地震的发生.
(2) 岷山隆起一带中下地壳的高泊松比异常呈“凸起”形态,结合前人研究发现的较高热流和岩石快速抬升现象,推测岷山隆起一带可能存在岩石圈的拆沉,导致地幔热物质上涌而形成下地壳高泊松比物质.
(3) NCDB和SCDB具有不同的S波速度和泊松比分布特征.30 km深的NCDB具有明显的低速异常,而该深度下并不具有明显的高泊松比值特征.此外剖面成像结果也显示CDB内的低速异常与高泊松比的分布不一致,因此与下地壳流模型不符.综合其它地震学证据,本文认为CDB的变形模式为上地壳纯剪切增厚,块体变形主要受块体内部的走滑断裂及活动边界断裂控制.
致谢 感谢中国地震台网中心、四川省地震局和云南省地震局为本研究提供原始波形数据.感谢两位审稿专家提出的宝贵修改意见和建议.文中大部分图件采用GMT(Wessel and Smith, 1995)绘制.
Bao X W, Sun X X, Xu M J, et al. 2015. Two crustal low-velocity channels beneath SE Tibet revealed by joint inversion of Rayleigh wave dispersion and receiver function. Earth and Planetary Science Letters, 415: 16-24. DOI:10.1016/j.epsl.2015.01.020 |
Barmin M P, Ritzwoller M H, Levshin A. 2001. A fast and reliable method for surface wave tomography. Pure and Applied Geophysics, 158(8): 1351-1375. DOI:10.1007/PL00001225 |
Beaumont C, Jamieson R A, Nguyen M H, et al. 2001. Himalayan tectonics explained by extrusion of a low-viscosity crustal channel coupled to focused surface denudation. Nature, 414(6865): 738-742. DOI:10.1038/414738a |
Bensen G B, Ritzwoller M H, Barmin M P, et al. 2007. Processing seismic ambient noise data to obtain reliable broad-band surface wave dispersion measurements. Geophysical Journal International, 169(3): 1239-1260. DOI:10.1111/j.1365-246X.2007.03374.x |
Burchfiel B C, Royden L H, Van Der Hilst R D, et al. 2008. A geological and geophysical context for the Wenchuan earthquake of 12 May 2008, Sichuan, People's Republic of China. GSA Today, 18(7): 4-11. DOI:10.1130/GSATG18A.1 |
Chen Y, Zhang Z J, Sun C Q, et al. 2013. Crustal anisotropy from Moho converted Ps wave splitting analysis and geodynamic implications beneath the eastern margin of Tibet and surrounding regions. Gondwana Research, 24(3-4): 946-957. DOI:10.1016/j.gr.2012.04.003 |
Clark M K, Royden L H. 2000. Topographic ooze:building the eastern margin of Tibet by lower crustal flow. Geology, 28(8): 703-706. |
Clark M K, Bush J W M, Royden L H. 2005. Dynamic topography produced by lower crustal flow against rheological strength heterogeneities bordering the Tibetan Plateau. Geophysical Journal International, 162(2): 575-590. DOI:10.1111/j.1365-246X.2005.02580.x |
Deng Q D, Chen S F, Zhao X L. 1994. Tectonics, scismisity and dynamics of Longmenshan Mountains and its adjacent regions. Seismology and Geology (in Chinese), 16(4): 389-403. |
Deng Q D, Cheng S P, Ma J, et al. 2014. Seismic activities and earthquake potential in the Tibetan Planteau. Chinese Journal of Geophysics (in Chinese), 57(7): 2025-2042. DOI:10.6038/cjg20140701 |
Dziewonski A, Bloch S, Landisman M. 1969. A technique for the analysis of transient seismic signals. Bulletin of the Seismological Society of America, 59(1): 427-444. |
Fan L P, Wu J P, Fang L H, et al. 2015. The characteristic of Rayleigh wave group velocities in the southeastern margin of the Tibetan Plateau and its tectonic implication. Chinese Journal of Geophysics (in Chinese), 58(5): 1555-1567. DOI:10.6038/cjg20150509 |
Herrmann R B. 1973. Some aspects of band-pass filtering of surface waves. Bulletin of the Seismological Society of America, 63(2): 663-671. |
Herrmann R B, Ammon C J. 2004. Surface waves, receiver functions and crustal structure.//Computer Programs in Seismology, Version 3.30. Saint Louis: Saint Louis University.
|
Houseman G, England P. 1993. Crustal thickening versus lateral expulsion in the Indian-Asian continental collision. Journal of Geophysical Research:Solid Earth, 98(B7): 12233-12249. DOI:10.1029/93JB00443 |
Hu J F, Su Y J, Zhu X G, et al. 2005. S-wave velocity and Passion's ratio structure of crust in Yunnan and its implication. Science in China Series D:Earth Sciences, 48(2): 210-218. DOI:10.1360/03yd0062 |
Hubbard J, Shaw J H. 2009. Uplift of the Longmen shan and Tibetan plateau, and the 2008 Wenchuan (M=7.9) earthquake. Nature, 458(7235): 194-197. DOI:10.1038/nature07837 |
Ingber L. 1989. Very fast simulated re-annealing. Mathematical and Computer Modelling, 12(8): 967-973. DOI:10.1016/0895-7177(89)90202-1 |
Jiang G Z, Gao P, Rao S, et al. 2016. Compilation of heat flow data in the continental area of China (4th edition). Chinese Journal of Geophysics (in Chinese), 59(8): 2892-2910. DOI:10.6038/cjg20160815 |
Kirby E, Harkins N. 2013. Distributed deformation around the eastern tip of the Kunlun fault. International Journal of Earth Sciences, 102(7): 1759-1772. DOI:10.1007/s00531-013-0872-x |
Kirkpatrick S, Gelatt C D Jr, Vecchi M P. 1983. Optimization by simulated annealing. Science, 220(4598): 671-680. DOI:10.1126/science.220.4598.671 |
Lev E, Long D M, Van Der Hilst R D. 2006. Seismic anisotropy in Eastern Tibet from shear wave splitting reveals changes in lithospheric deformation. Earth and Planetary Science Letters, 251(3-4): 293-304. DOI:10.1016/j.epsl.2006.09.018 |
Li H Y, Su W, Wang C Y, et al. 2009. Ambient noise Rayleigh wave tomography in western Sichuan and eastern Tibet. Earth and Planetary Science Letters, 282(1-4): 201-211. DOI:10.1016/j.epsl.2009.03.021 |
Liang C T, Langston C A. 2008. Ambient seismic noise tomography and structure of eastern North America. Journal of Geophysical Research:Solid Earth, 113(B3): B03309. DOI:10.1029/2007JB005350 |
Liang C T, Langston C A. 2009. Three-dimensional crustal structure of eastern North America extracted from ambient noise. Journal of Geophysical Research:Solid Earth, 114(B3): B03310. DOI:10.1029/2008JB005919 |
Liu Q Y, Van Der Hilst R D, Li Y, et al. 2014. Eastward expansion of the Tibetan Plateau by crustal flow and strain partitioning across faults. Nature Geoscience, 7(5): 361-365. |
Liu Z, Tian X B, Gao R, et al. 2017. New images of the crustal structure beneath eastern Tibet from a high-density seismic array. Earth and Planetary Science Letters, 480: 33-41. DOI:10.1016/j.epsl.2017.09.048 |
Luo J, Zhao C P, Lü J, et al. 2016. Characteristics of focal mechanisms and the stress field in the southeastern margin of the Tibetan Plateau. Pure and Applied Geophysics, 173(8): 2687-2710. DOI:10.1007/s00024-016-1350-8 |
Panza G F. 1981. The resolving power of seismic surface waves with respect to crust and upper mantle structural models.//Cassinis R ed. The Solution of the Inverse Problem in Geophysical Interpretation. Boston, MA: Springer, 39-77, doi: 10.1007/978-1-4684-3962-5_3.
|
Rothman D H. 1985. Nonlinear inversion, statistical mechanics, and residual statics estimation. Geophysics, 50(12): 2784-2796. DOI:10.1190/1.1441899 |
Roux P, Sabra K G, Kuperman W A, et al. 2005. Ambient noise cross-correlation in free space:theoretical approach. The Journal of the Acoustical Society of America, 117(1): 79-84. DOI:10.1121/1.1830673 |
Royden L H, Burchfiel B C, King R W, et al. 1997. Surface deformation and lower crustal flow in eastern Tibet. Science, 276(5313): 788-790. DOI:10.1126/science.276.5313.788 |
Sen M K, Stoffa P L. 1991. Nonlinear one-dimensional seismic waveform inversion using simulated annealing. Geophysics, 56(10): 1624-1638. DOI:10.1190/1.1442973 |
Shapiro N M, Ritzwoller M H. 2002. Monte-Carlo inversion for a global shear-velocity model of the crust and upper mantle. Geophysical Journal International, 151(1): 88-105. DOI:10.1046/j.1365-246X.2002.01742.x |
Shen Z K, Lü J N, Wang M, et al. 2005. Contemporary crustal deformation around the southeast borderland of the Tibetan Plateau. Journal of Geophysical Research:Solid Earth, 110(B11): B11409. DOI:10.1029/2004JB003421 |
Tapponnier P, Peltzer G, Le Dain A Y, et al. 1982. Propagating extrusion tectonics in Asia:new insights from simple experiments with plasticine. Geology, 10(12): 611-616. DOI:10.1130/0091-7613(1982)10<611:PETIAN>2.0.CO;2 |
Wang C Y, Lou H, Lü Z Y, et al. 2008. S-wave crustal and upper mantle's velocity structure in the eastern Tibetan Plateau-Deep environment of lower crustal flow. Science in China Series D:Earth Sciences, 51(2): 263-274. DOI:10.1007/s11430-008-0008-5 |
Wang S, Xu X Y, Hu J F. 2015. Review on the study of crustal structure and geodynamic models for the southeast margin of the Tibetan Plateau. Chinese Journal of Geophysics (in Chinese), 58(11): 4235-4253. DOI:10.6038/cjg20151129 |
Wang X, Chen L, Ai Y S, et al. 2018. Crustal structure and deformation beneath eastern and northeastern Tibet revealed by P-wave receiver functions. Earth and Planetary Science Letters, 497: 69-79. DOI:10.1016/j.epsl.2018.06.007 |
Wang X C, Ding Z F, Wu Y, et al. 2017. Crustal thicknesses and Poisson's ratios beneath the northern section of the north-south seismic belt and surrounding areas in China. Chinese Journal of Geophysics (in Chinese), 60(6): 2080-2090. DOI:10.6038/cjg20170605 |
Wapenaar K. 2004. Retrieving the elastodynamic Green's function of an arbitrary inhomogeneous medium by cross correlation. Physical Review Letters, 93(25): 254301. DOI:10.1103/PhysRevLett.93.254301 |
Weaver R L, Lobkis O I. 2001. Ultrasonics without a source:thermal fluctuation correlations at MHz frequencies. Physical Review Letters, 87(13): 134301. DOI:10.1103/PhysRevLett.87.134301 |
Wessel P, Smith W H F. 1995. New version of the Generic Mapping Tools. Eos, 76(33): 329. |
Xie J Y, Ritzwoller M H, Shen W S, et al. 2013. Crustal radial anisotropy across eastern Tibet and the western Yangtze craton. Journal of Geophysical Research:Solid Earth, 118(8): 4226-4252. DOI:10.1002/jgrb.50296 |
Xu L L, Rondenay S, Van Der Hilst R D. 2007. Structure of the crust beneath the southeastern Tibetan Plateau from teleseismic receiver functions. Physics of the Earth and Planetary Interiors, 165(3-4): 176-193. DOI:10.1016/j.pepi.2007.09.002 |
Xu X W, Wen X Z, Zheng R Z, et al. 2003. Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region, China. Science in China Series D:Earth Sciences, 46(2): 210-226. |
Xu X W, Wen X Z, Yu G H, et al. 2005. Average slip rate, earthquake rupturing segmentation and recurrence behavior on the Litang fault zone, western Sichuan Province, China. Science in China Series D:Earth Sciences, 48(8): 1183-1196. DOI:10.1360/04yd0072 |
Xu X W, Wen X Z, Chen G H, et al. 2008. Discovery of the Longriba fault zone in eastern Bayan Har block, China and its tectonic implication. Science in China Series D:Earth Sciences, 51(9): 1209-1223. DOI:10.1007/s11430-008-0097-1 |
Xu Z Q, Jiang M, Yang J S, et al. 1999. Mantle diapir and inward intracontinental subduction:A discussion on the mechanism of uplift of the Qinghai-Tibet Plateau. Geological Society of America Special Papers, 328: 19-31. |
Xu Z Q, Ji S C, Li H B, et al. 2008. Uplift of the Longmen Shan range and the Wenchuan earthquake. Episodes, 31(3): 291-301. DOI:10.18814/epiiugs/2008/v31i3/002 |
Yang Y H, Liang C T, Su J R. 2015. Focal mechanism inversion based on regional model inverted from receiver function and its application to the Lushan earthquake sequence. Chinese Journal of Geophysics (in Chinese), 58(10): 3583-3600. DOI:10.6038/cjg20151013 |
Yang Y H, Fan J, Hua Q, et al. 2017. Inversion for the focal mechanisms of the 2017 Jiuzhaigou M7.0 earthquake sequence using near-field full waveforms. Chinese Journal of Geophysics (in Chinese), 60(10): 4098-4101. DOI:10.6038/cjg20171034 |
Yang Y H, Liang C T, Fang L H, et al. 2018. A comprehensive analysis on the stress field and seismic anisotropy in Eastern Tibet. Tectonics, 37(6): 1648-1657. DOI:10.1029/2018TC005011 |
Yang Y J, Ritzwoller M H, Zheng Y, et al. 2012. A synoptic view of the distribution and connectivity of the mid-crustal low velocity zone beneath Tibet. Journal of Geophysical Research:Solid Earth, 117(B4): B04303. DOI:10.1029/2011JB008810 |
Yao H J, van Der Hilst R D, De Hoop M V. 2006. Surface-wave array tomography in SE Tibet from ambient seismic noise and two-station analysis-I. Phase velocity maps. Geophysical Journal International, 166(2): 732-744. DOI:10.1111/j.1365-246X.2006.03028.x |
Yao H J, Beghein C, Van Der Hilst R D. 2008. Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis-Ⅱ. Crustal and upper-mantle. Geophysical Journal International, 173(1): 205-219. DOI:10.1111/j.1365-246X.2007.03696.x |
Yao H J, Van Der Hilst R D, Montagner J P. 2010. Heterogeneity and anisotropy of the lithosphere of SE Tibet from surface wave array tomography. Journal of Geophysical Research, 115(B12): B12307. DOI:10.1029/2009JB007142 |
Yao Y. 1995. Improvement on nonlinear geophysical inversion simulated annealing. Acta Geophysica Sinica (in Chinese), 38(5): 643-650. |
Yi G X, Long F, Zhang Z W. 2012. Spatial and temporal variation of focal mechanisms for aftershocks of the 2008 MS8.0 Wenchuan earthquake. Chinese Journal of Geophysics (in Chinese), 55(4): 1213-1227. DOI:10.6038/j.issn.0001-5733.2012.04.017 |
Yi G X, Long F, Vallage A, et al. 2016. Focal mechanism and tectonic deformation in the seismogenic area of the 2013 Lushan earthquake sequence, southwestern China. Chinese Journal of Geophysics (in Chinese), 59(10): 3711-3731. DOI:10.6038/cjg20161017 |
Yi G X, Long F, Liang M J, et al. 2017a. Focal mechanism solutions and seismogenic structure of the 8 August 2017 M7.0 Jiuzhaigou earthquake and its aftershocks, northern Sichuan. Chinese Journal of Geophysics (in Chinese), 60(10): 4083-4097. DOI:10.6038/cjg20171033 |
Yi G X, Long F, Liang M J, et al. 2017b. Seismogenic structure of the M4.9 and M5.1 Litang earthquakes on 23 September 2016 in southwestern China. Seismology and Geology (in Chinese), 39(5): 949-963. |
Yin A, Harrison T M. 2000. Geologic evolution of the Himalayan-Tibetan orogen. Annual Review of Earth and Planetary Sciences, 28: 211-280. DOI:10.1146/annurev.earth.28.1.211 |
Zelt C A. 1998. Lateral velocity resolution from three-D seismic refraction data. Geophysical Journal International, 135(3): 1101-1112. DOI:10.1046/j.1365-246X.1998.00695.x |
Zhang X M, Teng J W, Sun R M, et al. 2014. Structural model of the lithosphere-asthenosphere system beneath the Qinghai-Tibet Plateau and its adjacent areas. Tectonophysics, 634: 208-226. DOI:10.1016/j.tecto.2014.08017 |
Zhao L F, Xie X B, He J K, et al. 2013. Crustal flow pattern beneath the Tibetan Plateau constrained by regional Lg-wave Q tomography. Earth and Planetary Science Letters, 383: 113-122. DOI:10.1016/j.epsl.2013.09.038 |
Zheng C, Ding Z F, Song X D. 2016. Joint inversion of surface wave dispersion and receiver functions for crustal and uppermost mantle structure in Southeast Tibetan Plateau. Chinese Journal of Geophysics (in Chinese), 59(9): 3223-3236. DOI:10.6038/cjg20160908 |
Zhu J S. 2008. The Wenchuan earthquake occurrence background in deep structure and dynamics of lithosphere. Journal of Chengdu University of Technology (Science & Technology Edition) (in Chinese), 35(4): 348-356. DOI:10.3969/j.issn.1671-9727.2008.04.002 |
Zhu J S, Zhao J M, Jiang X T, et al. 2012. Crustal flow beneath the eastern margin of the Tibetan plateau. Earthquake Science, 25(5-6): 469-483. DOI:10.1007/s11589-012-0871-1 |
Zhu J S, Wang X B, Yang Y H, et al. 2017. The crustal flow beneath the eastern margin of the Tibetan Plateau and its process of dynamic. Chinese Journal of Geophysics (in Chinese), 60(6): 2038-2057. DOI:10.6038/cjg20170602 |
邓起东, 陈社发, 赵小麟. 1994. 龙门山及其邻区的构造和地震活动及动力学. 地震地质, 16(4): 389-403. |
邓起东, 程绍平, 马冀, 等. 2014. 青藏高原地震活动特征及当前地震活动形势. 地球物理学报, 57(7): 2025-2042. DOI:10.6038/cjg20140701 |
范莉苹, 吴建平, 房立华, 等. 2015. 青藏高原东南缘瑞利波群速度分布特征及其构造意义探讨. 地球物理学报, 58(5): 1555-1567. DOI:10.6038/cjg20150509 |
姜光政, 高鹏, 饶松, 等. 2016. 中国大陆地区大地热流数据汇编(第四版). 地球物理学报, 59(8): 2892-2910. DOI:10.6038/cjg20160815 |
胡家富, 苏有锦, 朱雄关, 等. 2003. 云南的地壳S波速度与泊松比结构及其意义. 中国科学(D辑), 33(8): 714-722. |
王椿镛, 楼海, 吕智勇, 等. 2008. 青藏高原东部地壳上地幔S波速度结构——下地壳流的深部环境. 中国科学D辑:地球科学, 38(1): 22-32. |
王苏, 徐晓雅, 胡家富. 2015. 青藏高原东南缘的地壳结构与动力学模式研究综述. 地球物理学报, 58(11): 4235-4253. DOI:10.6038/cjg20151129 |
王兴臣, 丁志峰, 武岩, 等. 2017. 中国南北地震带北段及其周缘地壳厚度与泊松比研究. 地球物理学报, 60(6): 2080-2090. DOI:10.6038/cjg20170605 |
徐锡伟, 闻学泽, 郑荣章, 等. 2003. 川滇地区活动块体最新构造变动样式及其动力来源. 中国科学(D辑), 33(S1): 151-162. |
徐锡伟, 闻学泽, 于贵华, 等. 2005. 川西理塘断裂带平均滑动速率、地震破裂分段与复发特征. 中国科学D辑:地球科学, 35(6): 540-551. |
徐锡伟, 闻学泽, 陈桂华, 等. 2008. 巴颜喀拉地块东部龙日坝断裂带的发现及其大地构造意义. 中国科学D辑:地球科学, 38(5): 529-542. |
杨宜海, 梁春涛, 苏金蓉. 2015. 用接收函数建立区域模型的震源机制反演及其在芦山地震序列研究中的应用. 地球物理学报, 58(10): 3583-3600. DOI:10.6038/cjg20151013 |
杨宜海, 范军, 花茜, 等. 2017. 近震全波形反演2017年九寨沟M7.0地震序列震源机制解. 地球物理学报, 60(10): 4098-4104. DOI:10.6038/cjg20171034 |
姚姚. 1995. 地球物理非线性反演模拟退火法的改进. 地球物理学报, 38(5): 643-650. DOI:10.3321/j.issn:0001-5733.1995.05.011 |
易桂喜, 龙锋, 张致伟. 2012. 汶川MS8.0地震余震震源机制时空分布特征. 地球物理学报, 55(4): 1213-1227. DOI:10.6038/j.issn.0001-5733.2012.04.017 |
易桂喜, 龙锋, Vallage A, 等. 2016. 2013年芦山地震序列震源机制与震源区构造变形特征分析. 地球物理学报, 59(10): 3711-3731. DOI:10.6038/cjg20161017 |
易桂喜, 龙锋, 梁明剑, 等. 2017a. 2017年8月8日九寨沟M7.0地震及余震震源机制解与发震构造分析. 地球物理学报, 60(10): 4083-4097. DOI:10.6038/cjg20171033 |
易桂喜, 龙锋, 梁明剑, 等. 2017b. 2016年9月23日四川理塘M4.9和M5.1地震发震构造分析. 地震地质, 39(5): 949-963. |
张培震. 2008. 青藏高原东缘川西地区的现今构造变形、应变分配与深部动力过程. 中国科学D辑:地球科学, 38(9): 1041-1056. |
郑晨, 丁志峰, 宋晓东. 2016. 利用面波频散与接收函数联合反演青藏高原东南缘地壳上地幔速度结构. 地球物理学报, 59(9): 3223-3236. DOI:10.6038/cjg20160908 |
朱介寿. 2008. 汶川地震的岩石圈深部结构与动力学背景. 成都理工大学学报(自然科学版), 35(4): 348-356. DOI:10.3969/j.issn.1671-9727.2008.04.002 |
朱介寿, 王绪本, 杨宜海, 等. 2017. 青藏高原东缘的地壳流及动力过程. 地球物理学报, 60(6): 2038-2057. DOI:10.6038/cjg20170602 |