5 期 范文渊等 : 青藏高原东部和周边地区地壳速度结构的背景噪声层析成像 1569 meltinglayerintheanomalouslythickcrustiscriticaltounderstandthedynamicsandevolution historyoftheformationan

Transcription

1 第 58 卷第 5 期 2015 年 5 月 地球物理学报 CHINESE JOURNAL OF GEOPHYSICS Vol.58,No.5 May,2015 范文渊, 陈永顺, 唐有彩等 青藏高原东部和周边地区地壳速度结构的背景噪声层析成像. 地球物理学报,58(5): ,doi: /cjg Fan W Y,ChenY S,Tang Y C,etal.2015.CrustanduppermantlevelocitystructureoftheeasternTibetanPlateauand adjacentregionsfromambientnoisetomography. 犆犺犻狀犲狊犲犑. 犌犲狅狆犺狔狊.(inChinese),58(5): ,doi: /cjg 青藏高原东部和周边地区地壳速度结构的背景噪声层析成像 范文渊, 陈永顺, 唐有彩, 周仕勇, 冯永革, 岳汉, 王海洋, 金戈, 魏松峤, 王彦宾, 盖增喜, 宁杰远 北京大学地球与空间科学学院地球物理研究所, 北京 摘要利用连续地震背景噪声记录和互相关技术获得瑞利面波格林函数, 进而反演获得了青藏高原东部和周边地区的地壳三维速度结构. 地震数据源于北京大学宽频带流动观测地震台阵, 国家数字测震台网数据备分中心提供的部分固定台站的连续记录及 INDEPTHIV 宽频带流动观测地震台阵. 首先对观测数据进行处理和分析取得所有可能台站对的面波经验格林函数和瑞利波相速度频散曲线, 反演得到了观测台阵下方周期从 6~60s 的瑞利波相速度异常分布图像. 并且进一步反演获得研究区域三维剪切波速度结构和莫霍面深度分布. 短周期 (6~14s) 相速度异常分布与地表地质构造特征吻合较好, 在青藏高原和四川盆地之间存在一个明显的南北向转换带. 而本文最重要的结果是周期大于 25s 的相速度异常分布图像显示, 以昆仑断裂带为界, 柴达木盆地和祁连山脉地区呈现与青藏高原截然不同的中地壳速度结构, 反而与青藏高原东缘地区和川滇菱形块体速度结构相似. 反演获得的剪切波速度在 27.5~45km 深度的切片也明显地揭示 : 青藏高原的松潘 甘孜地块和羌塘地块呈现均一的低速层 ; 然而, 柴达木盆地和祁连山脉地区则呈现较强的横向不均匀性, 尤其是柴达木盆地的高速异常和四川盆地的高速异常相对应. 这些结果为前人提出的青藏高原东北向台阶式增长模式提供了重要的地震学观测证据. 与全球一维平均速度模型 (AK135) 相比较发现, 本文测量和反演获得的研究区域内平均相速度和剪切波速度都比 AK135 模型慢很多, 尤其是青藏高原的中地壳 (25~40km) 剪切波速度显著低于全球平均速度模型. 进一步的层析成像反演证实松潘 甘孜和羌塘地块中地壳 (27.5~45km) 呈现大范围均一的低速层, 为青藏高原可能存在大规模中下地壳 层流 提供地震学观测证据. 在祁连山脉的 27.5~45km 深度观测到的明显低速异常体可能对应于该造山带下地幔岩浆活动导致的底侵作用, 表明引起该地区地壳增厚的主要机制可能是来自地幔岩浆的底侵作用. 关键词背景噪声 ; 瑞利面波 ; 层析成像 ; 剪切波速度结构 ; 青藏高原 ; 柴达木盆地和祁连山脉地区 doi: /cjg 中图分类号 P315 收稿日期 , 收修定稿 犆狉狌狊狋犪狀犱狌狆狆犲狉犿犪狀狋犾犲狏犲犾狅犮犻狋狔狊狋狉狌犮狋狌狉犲狅犳狋犺犲犲犪狊狋犲狉狀犜犻犫犲狋犪狀犘犾犪狋犲犪狌犪狀犱犪犱犼犪犮犲狀狋狉犲犵犻狅狀狊犳狉狅犿犪犿犫犻犲狀狋狀狅犻狊犲狋狅犿狅犵狉犪狆犺狔 FAN Wen Yuan,CHEN Yong Shun,TANG You Cai,ZHOUShi Yong,FENG Yong Ge,YUE Han, WANG Hai Yang,JIN Ge,WEISong Qiao,WANG Yan Bin,GEZeng Xi,NINGJie Yuan 犐狀狊狋犻狋狌狋犲狅犳犜犺犲狅狉犲狋犻犮犪犾犪狀犱犃狆狆犾犻犲犱犌犲狅狆犺狔狊犻犮狊, 犛犮犺狅狅犾狅犳犈犪狉狋犺犪狀犱犛狆犪犮犲犛犮犻犲狀犮犲狊, 犘犲犽犻狀犵犝狀犻狏犲狉狊犻狋狔, 犅犲犻犼犻狀犵 , 犆犺犻狀犪 犃犫狊狋狉犪犮狋 ThecruststructureoftheTibetPlateauisstilcontroversial.Whetherthereisapartial 基金项目国家自然科学基金委项目 ( ; ) 资助. 作者简介范文渊, 男, 硕士研究生, 主要从事面波和地震震源研究. 通讯作者 陈永顺, 男, 教授, 从事地球动力学研究.E mail:johnyc@pku.edu.cn

2 5 期 范文渊等 : 青藏高原东部和周边地区地壳速度结构的背景噪声层析成像 1569 meltinglayerintheanomalouslythickcrustiscriticaltounderstandthedynamicsandevolution historyoftheformationand on going process. Weaddresstheissue withrobustlyinverted surfacewaveandshearwavevelocitymodelsfromthewestchinadensearraydata. Ambientnoisetomographywasappliedtoobtainthecrustanduppermantlevelocitymodels ofthe Tibetan Plateau.The methodisbased onextracting Rayleigh wavephaseand group velocitydispersioncurvefromcontinuouswaveformcross correlation.ourdatacomesfromthe Peking University seismic arrays,indepth IV passive seismic array,and four Chinese Provincialseismicnetworks.Theusabledatawascolectedfrom05,2007to12,2009.Thelong timespanassureshighsignaltonoiseratiowhenthecross correlationisperformed.weextract theempiricalgreen sfunctionsfromthecross correlationandperformasurfacewavetomography inversionwiththecross correlationfunctions.withthephaseandgroupvelocitymodelsfromthe previousstep,a3dshearwavevelocitymodelisobtained. Oursurfacephase/group velocity modelsandinvertedshearvelocity modelshow strong heterogeneityandcomplexity overthe wholestudiedregion.surface geologicalfeaturesare accuratelydepictedbytheshortperiodsurfacewaves(6sto14s,bothphaseandgroupvelocity model).theqaidambasinandsichuanbasinarewelcorrelatedtolowvelocityanomaliesdueto thethicksedimentsovertheseregions.mountains,likeqilian,kunlun,andlongmenshanare markedbyhighvelocityanomalies.wealsoobserveacontinuoustransitionzoneinbetweenthe TibetplateauandSichuanBasinatshortperiod.ForRayleigh wavetomography modelslonger than25s,therearesignificantdiferencesbetweenqaidam Qilianregionandtheplateausouthof Kunlunfaultzone.Thenewlyformedplateauarea,i.e.Qaidam Qilianregionsharessimilar velocitystructuresignatureswiththeeasternboundaryofthetibetanplateauandsichuanbasin. Bothourphasevelocity measurementsandinvertedshear wavevelocity modelsshowthatthe Tibetanplateauhasalarge scaleuniformlowvelocitylayerinthemiddlecrust (25~40kmin depth).thevelocityatthisdepthoftheregionismuchslowerthantheglobalaveragevelocityof AK135.ItisalsoshowsdistinctdeviationfromtheeasternboundaryoftheTibetanPlateau.The low velocity layer supports the channel flow model that explains the eastward material transportationoftibet.fromtheintermediateperiodvelocity maps,weobservealowvelocity zoneatqilian Mountainregionanditmightberelatedtotheelevatedtemperatureinthemiddle andlowercrustofthisregion.thevariationofthetemperatureatthisdepth mightcorrelateto theunderplatingofmagmafromthelocalmantleupweling.theseobservationssuggestthatthe on goingcrustalthickeningatqilian Mountainareaispossiblyduetotheunderplatingprocess. OurdetailedstructureofTibetanPlateaucrustiscriticaltounderstandtheon goingplateau formationprocess.therayleigh wavegroup/phasevelocity modelsand3dshearwavevelocity modelindicatethereisalarge scaleuniformlow velocity withintheplateau.thelow velocity layercouldexplaintheeastward materialtransportationoftibet.thenewlyformed Qaidam Qilianregion,asanewplateaushowsdistinctvelocitysignaturesinourvelocitymodels.Andthe low velocity zone below Qilian Mountain is possibly an indication ofthe ongoing magma underplatingprocessinthatregion. 犓犲狔狑狅狉犱狊 Ambientnoise;Rayleigh wave;phasevelocity Tomography;Shear wavevelocity structure;tibetanplateau;qaidambasinandqilianmountain

3 1570 地球物理学报 (ChineseJ.Geophys.) 58 卷 1 引言 平均海拔 5000m 的青藏高原是现今地球上最大的陆地高原之一. 它形成于约 50 Ma 之前的印度板块与欧亚板块的陆 陆碰撞 (Harrisonetal., 2000;Tapponieretal.,2001). 自北向南, 青藏高原可以大体上分为柴达木块体, 松潘 甘孜块体, 羌塘块体和拉萨块体 ; 昆仑断裂带, 金沙缝合带和班公怒江缝合带依次将这些块体分开. 时至今日, 印度板块仍以每年 6cm 的速度向北运动 (Bilhametal., 1997). 由陆 陆碰撞引起的压缩变形由南向北延伸至整个青藏高原, 导致其地壳增厚 (60~70km). 位于青藏高原向北新生长的柴达木块体正在经历着地形的隆升和地壳增厚 (Meyeretal.,1998). 然而, 青藏高原巨厚的地壳结构仍无定论. 国际上用于解释青藏高原的隆升和地壳增厚的两种主要的学说是在 20 世纪 80 年代提出的大陆逃逸模型 (Tapponnier etal.,1982) 和连续形变模型 (EnglandandSearle, 1986). 地震层析成像是用于研究地球内部地震波速度结构不均匀分布的最有效手段之一. 在过去的 30 年中, 地球物理学家利用宽频带流动观测地震台阵记录, 运用各种地震层析成像技术来探测青藏高原及其周边地区的地壳结构. 前人通过瑞利面波的频散特性大尺度地研究了整个欧亚大陆 (Curtisetal., 1998;FengandTeng,1983;RitzwolerandLevshin, 1998;Ritzwoleretal.,1998;WuandLevshin,1994), 也曾利用更多数据, 采用更小网格精细地获得区域面波层析成像结果 (BourjotandRomanowicz,1992; Brandonand Romanowicz,1986;Griotetal.,1998; Huangetal.,2003;Romanowicz,1982). 被认为青藏高原东部物质溢出出口的川滇菱形块体也被很多学者所研究 (Huangetal.,2002;Yaoetal.,2008). 同时, 青藏高原内有诸多体波研究层析成像结果. 若干研究 (HuangandZhao,2006,Xuetal.,2002) 使用区域和远震地震数据来反演青藏高原及其周边区域的三维 P 波结构. 还有学者 (BaronandPriestley, 2009;Hearnetal.,2004;Liangand Song,2006; Liangetal.,2004) 采用 Pn 波和 Sn 波来研究青藏高原地壳底部的横向不均匀性和各向异性. 然而, 已有的研究结果对于全面了解青藏高原及其周边地区的地壳速度结构的横向变化仍存在较大的局限性.(1) 限于台站覆盖的区域和空间分辨 率, 现有结果对所感兴趣区域浅部结构分辨率有限 ; 而且大多数流动台站分布于藏南地区, 柴达木块体及青藏高原北缘的流动台站观测稀缺.(2) 基于传统方法的面波研究因为短周期信号衰减较快, 对地壳浅部结构分辨率有限 ; 而采用天然地震的体波成像结果对纵向分辨率略逊于面波成像.(3) 除少数反射折射地震学方法所得的二维地壳成像结果外, 整个青藏高原及其周边地区的高分辨率地壳三维成像结果稀缺. 所幸近几年布设结束的 INDEPTHIV 的流动台阵覆盖了大部分青藏高原东部地区和柴达木块体. 本文所用数据源自多个流动地震台阵和固定台网, 覆盖了青藏高原东部和周边地区 ( 图 1). 采用背景噪声层析成像方法对连续地震记录数据处理, 得到了研究区域的地壳瑞利面波相速度异常分布图像. 研究目的是利用连续记录的背景噪声数据和相应的相速度层析成像技术, 反演获得青藏高原东部地区和东缘的三维地壳结构. 进一步, 通过瑞利面波相速度分布结果反演得到研究区域的三维剪切波速度结构和莫霍面深度分布. 本研究所得青藏高原东部和周边地区的地壳三维地震波速度结构, 为进一步深入讨论青藏高原形成和生长的动力学机制提供了重要地震观测约束. 2 数据及处理流程 本文所用的连续记录地震资料来源于 : 在过去 5 年里完成的 3 个北京大学宽频带流动观测地震台阵, 包括北京大学银川台阵, 于田台阵和东昆仑台阵 ; 国家数字测震台网数据备分中心提供的甘肃 青海 四川及西藏台网的固定台站的连续记录及 INDEPTHIV 宽频带流动观测地震台阵 (Zhenget al.,2010). 图 1 为地震台阵分布图. 表 1 为所用地震台阵及台网的仪器信息. 由表 1 可以看出, 不同地震台阵 台网之间至少有 10 个月的记录时间重叠, 得以保障最终所得各台站对之间的互相关函数信噪比 (SNR) 较高. 表 1 还显示, 不同的台阵 台网所采用的仪器不同, 所以在做互相关之前, 需要去除仪器响应, 将记录转换为位移记录. 由于研究目标是获取瑞利面波信号, 所以仅采用垂向记录. 同时, 将记录时间全部换算至格林威治时间 (GMT). 这些数据的预处理保证了之后获取的瑞利面波群速度和相速度频散曲线的一致性和准确性.

4 5 期 范文渊等 : 青藏高原东部和周边地区地壳速度结构的背景噪声层析成像 1571 图 1 台站分布及研究区域地质块体边界 Fig.1 Stationmapofstudyingregion 台阵 北京大学银川台阵 北京大学于田台阵 北京大学东昆仑台阵 INDEPTHIV 台阵 家数字测震台网数据备分中心台网 表 1 地震台阵及仪器信息 犜犪犫犾犲 1 犛狋犪狋犻狅狀犪狀犱犻狀狊狋狉狌犿犲狀狋犾犻狊狋 记录日期 台站数目 地震仪类型 CMG 3ESP CMG 3ESPC CMG 3ESPC CMG 3ESP STS 2ESPCMG CMG 3ESPC BBVS 60CTS 1JCZ 1CTS 1EKS 2000M 背景噪声层析成像 (AmbientNoiseTomography) 近年来在地震大地构造学研究中得到广泛的应用, 从地震尾波的互相关 (CampiloandPaul,2003; Pauletal.,2005)(SeismicCoda) 或背景噪声的互相关 (Shapiroetal.,2005;ShapiroandCampilo, 2004)(AmbientNoise) 函数中, 可提取出面波经验格林函数. 近几年来被广泛地应用于区域层析成像研究 (Shapiroetal.,2005;KangandShin2006) 和大陆尺度的层析成像研究 (Yangetal.,2007;Zheng etal.,2008). 相关文献 (ShapiroandCampilo,2004; Campilo,2006;Lobkisand Weaver,2001;Snieder, 2004;Wapenaar,2004) 就背景噪声层析成像的理论和方法做过详细讨论. 这些研究成果是本文的理论依据. 本文数据处理主要依据 Bensen 等 (Bensenet al.,2007) 发展的由背景噪声数据测量瑞利面波相 速度频散曲线及相速度层析成像的方法及处理流程. 首先, 对一个台站对的垂向记录做互相关运算, 之后从互相关函数里获取该台站对之间的经验格林函数. 并且对所获取的台站对经验格林函数做带通滤波, 以保留高信噪比记录. 本文采用 Bensen(Bensen etal.,2007) 所提出对信噪比的定义, 对噪声即 在一给定时窗内的面波信号与之后的噪声信号的均方根之比, 对面波为最大振幅之比. 图 2 为叠加正负部分后的台站 A01 与其余台站之间的互相关函数. 其次, 采用时频分析方法 (FTAN) 获取所有台站对之间的周期从 6~60s 的瑞利面波群速度及相速度频散曲线. 台站对 A01 和 XHA( 图 1) 的群速度频散曲线如图 3 所示. 该台站对跨越青藏高原东缘 ( 图 1 中的浅蓝色虚线 ( 张培震等,2003)), 由图 3 可以看出, 在 12s 周期附近群速度出现局部极值, 这与其特殊路径相关. 基于台站对距离需大于三个波长的挑选准则, 因此对于本文使用的地震台站分布 ( 图 1), 大于 60s 周期的群速度无法可靠的从互相关函数所得的经验格林函数中获取. 最后, 在得到了所有台站对之间的瑞利面波群速度和相速度的频散曲线之后, 采用 Barmin(Barminet al.,2001) 提出的面波层析成像反演方法来获取研究区域的瑞利面波群速度和相速度的速度异常分布. 计算格点为 , 所用平滑阻尼参数依据所采用反演格点而定. 瑞利面波相速度的变化对介质密度和 P 波速度变化不敏感, 而对 S 波速度的变化非常敏感 (Takeuchietal.,1964). 基于所得的瑞利面波相速

5 １５７２ 地 球 物 理 学 报 Ｃｈ ｎｅ ｅｊ Ｇｅ ｏｐｈｙ ５８ 卷 图 ２ Ａ０１ 及其台站对之间的互相关函数 ５ １００ Ｆ ２ Ｃ ｏ ｃ ｏ ｅ ａ ｏｎｏ ｆａ０１ａｎｄ ａ ｏｎｐａ ｇ 度层析成像 反演得 到 研 究 区 域 的 三 层 地 壳 剪 切 波 利面波频散曲线反演获得该区域内的一维剪切波平 速度模型及莫霍深 度 分 布 剪 切 波 速 度 模 型 反 演 对 均速度模型 之后以 获 得 的 平 均 速 度 模 型 为 初 始 参 初始参考模型敏感 而 本 文 的 研 究 区 域 涵 盖 面 积 较 考模型 在各子区域 内 计 算 每 一 个 格 点 的 一 维 剪 切 大 包括青 藏 高 原 东 部 和 周 边 地 区 其 地 质 构 造 复 波速度及地壳厚度 最 终 将 各 格 点 反 演 获 得 的 一 维 杂 因此我们将研究 区 域 以 昆 仑 断 裂 带 ３６ Ｎ 和 龙 速度结果 拼 合 可 得 研 究 区 域 的 三 维 剪 切 波 速 度 模 门山 断 裂 带 １０２ Ｅ 为 界 划 分 为 ３ 个 子 区 域 图 型 在 剪 切 波 速 度 反 演 过 程 中 利 用 Ｓａ ｏ 的程序 ４ 青 藏 高 原 新 造 高 原 即 柴 达 木 地 块 青 藏 高 原 东缘地区 在 ３ 个子区域内 分别依据区域内平均瑞 Ｓａ ｏ １９８８ 计 算 理 论 相 速 度 和 其 对 速 度 的 偏 微 分 在线性反演过程中 用一系列水平层来描述地球

6 5 期 范文渊等 : 青藏高原东部和周边地区地壳速度结构的背景噪声层析成像 1573 在 ±2.5% 之间变化. 反演格点为 , 输入模型和检测结果如图 5 所示. 从检测结果可以看出, 在周期为 6s 和 60s 范围内, 青藏高原的东缘和北缘可以较好的得以重建, 但其余地区分辨率略低. 在 8~55s 周期范围, 本文的地震台站分布和反演方法能够在研究区内得到比较可靠而且稳定的结果, 包括祁连山 柴达木盆地 松潘 甘孜地块 羌塘地块 拉萨地块的北部和青藏高原东缘地区, 速度异常均可被精确的分辨重建. 在台阵内部, 或射线覆盖好的区域, 分辨率小于 200km, 不仅能够探测到速度异常体的形态, 而且能较好地反映异常体的速度差异. 而在射线覆盖不好的区域, 分辨率随着周期的增加逐渐降低. 3.2 群速度和相速度层析成像图 6 为观测到的研究区内瑞利面波平均群速度和相速度频散曲线. 与全球一维平均速度模型 AK135 的频散曲线相比较可以得出一个很重要的结论 : 本文测量的研究区域 ( 包括青藏高原 柴达木地块和青藏高原东缘 ) 内平均相速度都比 AK135 模型的相速度慢很多. 这可以部分归因于研究区内增厚的地壳 ( 主要影响在 30~60s 范围 ), 但是对应于 6~20s 范围则可能是由于研究区内物质不同所致 ( 如存在比平均大陆厚很多的沉积层 ). 这一点在后面展示反演获得的一维剪切波速度模型时还要进行 图 3 不同频段经验格林函数 (a) 及台站对 A01 XHA 群速度频散曲线 (b) Fig.3 (a) Empirical Green s functions in diferent frequencybandsand (b)groupvelocitydispersioncurse ofstationpaira01 XHA 介质, 每层厚度 10km 到 30km 不等. 假定这些水平层厚度不变, 拟合频散曲线来反演每一层的剪切波速度. 另外, 为了得到莫霍面的深度, 莫霍面上下两层的厚度也设为待反演的参数, 即在反演过程中可以在一定范围内自由变化. 3 结果 3.1 检测板测试为了检验本文层析成像结果的分辨率和可靠性, 按照惯例做了检测板测试. 输入检测板的速度异常体大小约为 200km(2 2 ), 输入模型为相速度, 详细讨论. 图 7 和图 8 的各个周期图像中的二维速度异常分布均相对于图 6 中的平均速度频散曲线而言. 不同周期的群速度和相速度的敏感核对应于不同深度, 同时不同周期的群速度和相速度的敏感核又有部分重叠. 长周期的速度敏感核较为宽泛, 因而所得的不同周期的瑞利波相速度异常变化图像随周期加大而趋于平缓. 经验而言, 基阶瑞利波群速度和相速度对约 1/3 波长深度附近介质的 S 波速度结构最为敏感 ; 面波成像反应平均结果, 而非特定深度结果. 图 7 和图 8 分别展示了瑞利面波群速度异常和相速度异常分布图像. 图 7 中的群速度异常与图 8 中的大多数相速度异常一致, 些许不同源于群速度和相速度不同的敏感核. 短周期相速度异常分布图像中 ( 图 8,6~14s), 高速异常与低速异常同地表地质构造相关性较高. 例如, 柴达木盆地和四川盆地的范围被低速异常精确地勾绘出来. 短周期图像反映了上地壳的地震波速度特征, 因此盆地的低速异常反映了这两个盆地内大范围的沉积层. 柴达木盆地和四川盆地沉积层

7 1574 地球物理学报 (ChineseJ.Geophys.) 58 卷 图 4 研究子区域划分 Fig.4 Sub regionsdivisionofstudyingarea 厚度的差异也为两块低速异常体所描述, 相较于四川盆地的低速异常消失于 12s 周期之后的分布图像, 柴达木较厚沉积层的影响, 使该地区低速异常持续至更长周期. 山脉为高速异常体所勾绘, 例如四川盆地西部的岷山和龙门山山脉均显示为条带状高速异常. 昆仑山和祁连山对短周期相速度异常图像有同样的影响, 但其高速异常被低速异常间或分隔为块体状而非条带状. 在祁连山脉和昆仑山脉交汇的地方, 共和盆地的沉积层为 8~12s 异常分布图像中的低速异常所描述. 基于本文数据分辨率, 青藏高原主要显现低速异常. 同时, 不同地块之间存在差异, 相较于羌塘和拉萨地块的较为均一的低速异常体, 松潘 甘孜的低速异常间或被小的高速异常分割. 使用较大阻尼系数加强平滑反演后仍得到相似结果. 尤其值得注意的是, 在青藏高原东缘并无一明显边界来界定青藏高原和四川盆地, 而是在两者之间存在若干小的速度异常体. 在图 8,10s 周期的相速度异常分布图像中, 可以看到 A 和 B 两块速度异常体.A 位于尼玛 玉树 鲜水河断裂带和金沙 红河断裂带的交汇点, 走向位于喜马拉雅东构造结的北北东延长线上, 在 A 和青藏高原东缘地区之间为川滇菱形块体. 川滇菱形块体在短周期图像中, 呈现为高速异常, 但随周期加大, 异常逐渐转换为低速异常.B 位于青藏高原东缘地区, 鄂尔多斯地块和柴达木地块的交汇处, 为低速异常. 中周期相速度异常分布图像 ( 图 8,16~35s) 描 绘了更深的速度结构.16s 及更长周期的图像与较短周期的图像差异很大. 这是由于 16~30s 这段周期的相速度对 30~40km 深的速度结构最为敏感, 而此深度也正好是平均大陆莫霍面的深度. 青藏高原的莫霍面远深于这个深度, 因此, 在青藏高原与青藏高原东缘地区之间应该存在明显速度异常差异. 如图 8, 中周期相速度异常分布图像所示, 所得结果与预期一致, 高速区对应于莫霍面较浅的地区 ( 如青藏高原东缘地区 ), 而低速异常对应于莫霍面较深的地区 ( 青藏高原 ). 短周期图像中存在的青藏高原东部与青藏高原东缘地区之间的若干破碎状速度异常体逐渐融合. 在短周期图像中松潘 甘孜地块间或出现的高速异常逐渐消失, 整个块体呈现均一低速异常 ; 但仍与羌塘和拉萨地块存在差异. 整个青藏高原呈现低速异常. 青藏高原东缘地区呈现均一高速异常, 这是因为下地壳速度结构影响了该地区的中周期的速度异常分布. 相似的低速异常转为高速异常的特点也出现在柴达木盆地和祁连山下. 但其低速异常的消失较四川盆地缓慢, 低速异常一直影响至 30s 周期的相速度异常分布图像. 这可能是源于柴达木盆地更厚的沉积层. 从 30s 周期图像开始, 明显的青藏高原东缘和北缘边界出现. 速度异常的对比勾勒出青藏高原与周边地体的边界. 前面提到的高速异常 A 消失于 30s 周期的相速度异常图像之后. 同时, 川滇菱形块体呈现低速异常. 但这一异常以鲜水河断裂和金沙红河断裂为界并未与青藏高原的低速异常融合. 从图 8,20s 的相

8 ５期 范文渊等 青藏高原东部和周边地区地壳速度结构的背景噪声层析成像 １５７５ 图 ５ 相速度层析成像检测板实验 Ｆ ５ Ｃｈｅ ｃｋｅ ｂｏａ ｄ ｅ ｆ ｏ ｐｈａ ｅｖｅ ｏ ｃ ｏｍｏｇ ａｐｈｙ ｇ ｙ 速度异常分布图像 中 可 以 看 出 在 柴 达 木 盆 地 以 北 的祁连山脉中出现 东 西 走 向 低 速 异 常 条 带 标 记 这 长周期相速度异常分布图像中 图 ８ ４ ０ ６ ０ 青藏高原东缘地区和柴达木地块的速度异常反映了 一低速异常体为 Ｃ 它 的 出 现 可 能 与 祁 连 山 脉 下 物 上地幔结构信息 而青藏高原则为下地壳信息 祁连 质流变性质有关 也可能与过深的莫霍深度相关 山脉中的低速异常 Ｃ 逐 渐 消 失 祁 连 山 脉 整 体 高 速

9 1576 地球物理学报 (ChineseJ.Geophys.) 58 卷 居中. 反演所得子区域的一维剪切波速度模型特征与其对应的相速度频散曲线特征一致, 青藏高原除在最浅层剪切波速度高于其他两个子区域外, 在其余层中 ( 中下地壳 ) 速度均为最低. 尤其注意到青藏高原在 25~40km 深度范围内存在显著低速层, 这是否对应于青藏高原中地壳的 层流 (Channel flow)? 该结果与藏南瑞利面波所得结果一致 (Jiangetal.,2011). 由图 9 中还可看出青藏高原的平均地壳厚度为 70km, 比其他两个子区域厚约 图 6 平均相速度群速度频散曲线 Fig.6 Averagedispersioncurvesofgroupand phasevelocity 异常强度减弱, 趋于均值. 青藏高原东缘地区呈现均一高速异常. 川滇菱形块体下的低速异常消失, 出现高速异常, 并与青藏高原东缘地区的高速异常融合. 青藏高原仍然呈现低速异常, 但与其东缘和北缘地区速度异常相比强度减弱, 高速异常与低速异常数值均减小趋于平均值. 青藏高原东部和西部速度结构出现差异, 西部速度异常强度大于东部异常. 在 40~60s 周期范围内的速度异常图像中, 速度异常结构主要为莫霍面深度控制, 因此青藏高原的东西速度差异可能源于东西不同的地壳厚度. 同时, 青藏高原下的印度板块俯冲深度不同 (Lietal.,2008), 也可能引起青藏高原南 中部地区地壳底部物质成分及温度的东西差异. 这些原因均可能引起青藏高原在长周期的相速度异常分布图像中出现东西差异, 需要更多观测证据约束. 3.3 剪切波速度结构和莫霍面深度根据瑞利面波相速度分布结果进一步反演获得剪切波速度分布, 可以更直接地认识地下速度结构. 研究区域内青藏高原东缘地区, 柴达木地块和青藏高原的地壳厚度差异很大, 而不同的初始速度模型在我们采用的反演剪切波速度结构及莫霍面深度的算法中对结果影响较大. 因而, 将研究区域划分成 3 个子区域 ( 图 4). 在不同子区域内分别进行剪切波速度及莫霍面深度反演, 最终将结果拼合便得研究区域内的三维剪切波速度模型. 图 9 为 3 个子区域的平均相速度频散曲线及反演所获得的一维平均剪切波速度结构. 由图 9 可以看出 3 个子区域存在明显差异 : 除最短周期时青藏高原的平均相速度最高外, 其余周期的相速度青藏高原均为最低, 青藏高原东缘地区最高, 柴达木地块 10km. 图 10 展示了剪切波速度分布的 3 个深度切片及莫霍面深度分布图像. 射线密度决定反演所获得的相速度层析成像的精度, 而相速度频散的准确性限定了反演所得的剪切波速度结构的准确性. 为防止误差放大, 我们仅反演计算了射线密度大于 500 条区域的剪切波速度模型. 如图 10 所示, 所得结果与图 8 所示的瑞利面波相速度层析成像结果特征比较一致. 在图 10 深度为 7.5km 的切片上可以看出, 剪切波速约为 3.0km s -1 的剪切波速约为 3.2km s -1 的剪切波速将四川盆地和柴达木盆地与青藏高原分离. 在深度为 27.5km 的切片中可知, 其速度分布与图 8 相速度异常分布的中周期部分对应. 前面提到的高速异常 A 也在该切片显现. 在此深度下, 剪切波波速在青藏高原中为 3.2km s -1, 与相速度异常层析成像中的低速异常相对应. 其中若干速度异常分块与青藏高原内部的横向不均匀性相对应. 从深度为 45km 的切片中可知, 青藏高原内部速度分布均一, 但绝对剪切波速度显示其仍为中地壳速度值, 而青藏高原东缘地区和柴达木块体则已经达到了下地壳的速度值. 限于我们采用的计算莫霍面方法, 所得莫霍面深度结果较接收函数等方法而言精度较差. 但是莫霍面深度分布模式还是与其他研究结果有较好吻合. 在松潘 甘孜地块内, 青藏高原的莫霍面深度达到最大值, 约 80km; 柴达木盆地的地壳厚度为 50km 左右, 较青藏高原东缘地区深. 川滇菱形块体的中部莫霍面深度较深. 祁连山地壳厚度较深, 为 70km 左右, 这可能是孕育出现图 8 中低速特征 C 的原因之一. 4 讨论 4.1 青藏高原 : 松潘 甘孜地块, 羌塘地块和拉萨地块青藏高原主要由 3 个地块构成 : 松潘 甘孜地

10 ５期 范文渊等 青藏高原东部和周边地区地壳速度结构的背景噪声层析成像 １５７７ 图 ７ 瑞利波群速度速度异常层析成像 Ｆ ７ Ｇ ｏｕｐｖｅ ｏ ｃ ｏｍｏｇ ａｐｈｙｏ ｆｒａｙ ｅ ｇ ｙ ｇｈｗａｖｅ 块 羌塘 地 块 和 拉 萨 地 块 ＤｅＣｅ ｅ ｅ ａ ２００２ 纪时期自冈瓦纳古 陆 漂 移 至 此 后 于 晚 侏 罗 纪 中 白 松潘 甘孜地块地壳主要由晚三叠纪早侏罗纪闭合 垩纪时期与羌塘地 块 碰 撞 拉 萨 地 块 地 壳 主 要 由 中 大洋残留海盆的 浊 流 岩 构 成 羌 塘 地 块 位 于 松 潘 级及高级变质沉积岩 和古 生代 沉积 岩构 成 ＤｅＣｅ ｅ 甘孜地块南部 于晚 古 生 代 时 期 自 冈 瓦 纳 古 陆 漂 移 至 北 后 于 晚 三 叠 纪 早 侏 罗 纪 时 期 与 松 潘 甘 孜 地 ｅ ａ ２００２ 本文研究结果表明青藏高原的 ３ 个块体东西向 块碰撞 羌塘地块地 壳 主 要 由 变 质 岩 和 古 生 代 及 中 速度结构均有差异 由 图 ８ 中 短 周 期 相 速 度 异 常 分 生代的沉积岩构成 拉萨地块位于最南部 于晚三叠 布 图 像 可 知 松 潘 甘 孜 地 块 上 地 壳 呈 现 强 烈 东 西

11 １５７８ 地 球 物 理 学 报 Ｃｈ ｎｅ ｅｊ Ｇｅ ｏｐｈｙ ５８ 卷 图 ８ 瑞利波相速度速度异常层析成像 Ｆ ８ Ｐｈａ ｅｖｅ ｏ ｃ ｏｍｏｇ ａｐｈｙｏ ｆｒａｙ ｅ ｇ ｙ ｇｈｗａｖｅ 向不均匀性 若干 低 速 异 常 体 被 高 速 体 分 隔 同 时 渡带内不仅存在东 西 向 速 度 结 构 差 异 同 时 存 在 南 居于青藏高原东缘地区和青藏高原之间的川滇菱形 北向速度结构差异 过 渡 带 内 的 东 西 向 速 度 异 常 差 块体与周边块体速度结构差异很大 因此 我们认为 异可能是青藏高原 物 质 向 东 非 均 匀 溢 出 的 结 果 同 在青藏高原和青藏 高 原 东 缘 地 区 的 上 地 壳 部 分 二 时 过 渡 带 内 南 北 向 的 速 度 异 常 差 异 可 能 是 源 于 ３ 者并没有一个明显的边界 而是存在一个过渡带 如 个块体的不同物质 差 异 在 图 ８ 的 中 长 周 期 相 速 度 在图 ８ １２ 周 期 速 度 异 常 分 布 图 像 中 Ｌ ｎｅ Ｒ 和 异常图像中可以看出 这一过渡带逐渐消失 最终融 Ｌ ｎｅ Ｌ 所 圈 定 的 区 域 这 一 过 渡 带 宽 约 ５００ｋｍ 过 入青藏高原的低速异常 因此 青藏高原中下地壳较

12 5 期 范文渊等 : 青藏高原东部和周边地区地壳速度结构的背景噪声层析成像 km 2, 北边界为祁连山南山山脉, 盆地沉积层厚度为 8~15km(GuandDi,1989). 图 8 周期为 8s 的相速度异常分布图像中, 柴达木盆地形状精确的被低速异常所描绘. 从图 8 中的周期为 8~ 20s 的图像中可以看出, 盆地西部低速异常消失略快于东部低速异常, 相对应于盆地东 西部的不同沉积层厚度. 柴达木盆地东部的共和盆地为一年轻的盆地 (Paresetal.,2003), 它对图 8, 相速度异常图像的影响主要体现在 8~12s 的范围. 将两个盆地分开的高速异常体与地表出露的古生代和中生代的火成岩相关 (Burchfieletal,1991). 若干研究 (Burchfiel etal.,1991;tapponnieretal.,1990;métivieret al.,1998) 指出, 共和盆地在晚上新世有过明显的抬升运动. 需要指出的是, 两个地区上地壳本身存在差异, 地壳厚度是可能造成所观测到差别的一种解释. 图 8 中的周期大于 25s 的图像显示, 以昆仑断裂带为界, 被称为 新增高原 的柴达木盆地和祁连山脉地区呈现与青藏高原截然不同的速度结构, 反而与青藏高原东缘地区和川滇菱形块体速度结构相似. 反演获得的剪切波速度在 45km 的切片 ( 图 10) 也明显地揭示 : 青藏高原的松潘 甘孜和羌塘地块都呈现均一的低速 (3.4km s -1 ) 结构 ( 图 9); 然而, 柴达木盆地和祁连山脉地区则呈现较强的横向不均匀性, 尤其是柴达木盆地的高速异常和四川盆地的高速异常相对应. 本文反演结果所显示的柴达木盆地和祁连山脉地区呈现与青藏高原截然不同的速度结构, 支持前 图 9 子区域平均相速度频散曲线 (a) 及一维剪切波模型 (b) Fig.9 Averagephasedispersioncurve ofthreesub regions 上地壳要均一. 最后, 根据瑞利波相速度异常图 ( 图 8) 可以看出羌塘地块和拉萨地块速度结构差异很小, 但与松潘 甘孜地块有较大差异. 剪切波反演结果 ( 图 10) 显示松潘 甘孜和羌塘地块中地壳 (27.5~45km) 呈现均一的低速结构, 可能对应于大范围的中地壳低速层 ( 图 9), 为青藏高原中下地壳可能存在大规模 层流 提供地震学观测证据. 类似地震学观测研究也在藏南推测存在中下地壳层流 (Jiangetal.,2011). 4.2 新增高原 : 柴达木盆地和祁连山脉位于昆仑山山脉以北的柴达木盆地占地面积 人提出的青藏高原东北向台阶式增长模式 (Tapponnier etal.,1990;métivieretal.,1998), 即青藏高原东北部边界从昆仑山山脉跨越到北祁连山脉, 使得柴达木盆地和祁连山脉地区成为新增高原. 按照这个模式, 现今柴达木盆地和祁连山脉地区的构造运动和隆升过程也就是青藏高原早期的形成模式. 根据柴达木盆地和祁连山脉地区的相速度频散曲线更接近于青藏高原东缘地区 ( 图 9), 我们推测该新增高原还仅仅处于初始阶段, 这也和柴达木盆地平均海拔高度仅为 2500m 相吻合. 图 9 显示祁连山脉在 27.5~45km 深度范围内呈现明显低速异常, 对应于相速度异常分布 ( 图 8) 从 20~45s 的低速异常体 C. 野外地质调查显示祁连山脉主要由古生代的变质岩构成, 该古老的缝合带现在受青藏高原南北向的挤压而重新活跃 隆升造山 ( 万景林 etal.,2010; 宋述光 etal.,2004). 推

13 １５８０ 地 球 物 理 学 报 Ｃｈ ｎｅ ｅｊ Ｇｅ ｏｐｈｙ ５８ 卷 图 １０ 三维剪切波模型及莫霍面深度分布图 Ｆ １０ Ｓｈｅ ａ ｗａｖｅ ｏｍｏｇ ａｐｈｙａｎｄ Ｍｏｈｏｄｅｐ ｈｍａｐ ｇ 测在祁连山脉 地 区 ２７ ５ ４５ｋｍ 深 度 范 围 内 发 现 的显著低速异常体对应于该造山带下可能的地幔岩 Ｃｈｅｎｅ ａ １９９４ Ｃｈｅｎａｎｄ Ｗ ｏｎ １９９６ 提 出 中下地壳逆冲形成 的 推 覆 体 来 解 释 龙 门 山 岷 山 的 浆活动造成的底侵作用 ｕｎｄｅ ｐ ａ ｎｇ 反演 同时 获 得的莫霍面深 度 分 布 图 １０ 也 显 示 祁 连 山 脉 地 区 隆升机制 龙门山逆 冲 带 的 变 形 历 史 可 以 追 溯 至 早 莫霍面显著变深 说 明 该 地 区 的 地 壳 增 厚 主 要 是 由 Ｗａｎｇｅ ａ １９８９ 岷 山 位 于 松 潘 甘 孜 褶 皱 带 内 长约 １８０ｋｍ 宽 ３０ ５０ｋｍ 从 新生 代起 至第 四 底侵作用引起的 同 一 地 区 的 接 收 函 数 研 究 也 揭 示 祁连山脉地区莫霍面变深 Ｙｕｅｅ ａ ２０１２ ４ ３ 四川盆地 川滇菱形块体和龙门山山脉 从柴达木盆地向东南看去 松潘 甘孜褶皱带 龙门山逆 冲 带 和 西 四 川 盆 地 构 建 起 青 藏 高 原 的 东 缘 龙门山的前缘恰 好 为 青 藏 高 原 和 华 南 扬 子 克 拉 白垩 后 印 欧 陆 陆 碰 撞 之 后 在 中 新 世 又 重 新 活 跃 纪为一上升 带 Ｌｏｎｇ １９８３ 四 川 盆 地 内 的 沉 积 岩 大多为第三纪的沉积岩 比柴达木盆地内沉积岩老 四川 盆 地 的 沉 积 层 厚 度 约 为 １０ ｋｍ Ｘｕｅ ａ ２００７ Ｘｕｅ ａ ２００７ 从 图 ８ 中 的 瑞 利 波 短 周 期 相速度速度异 常 分 布 图 像 和 图 １０ 中 的 ７ ５ｋｍ 深 通的碰撞边界 其长约 ５００ｋｍ 为四川盆 地西缘 龙 门山 岷山地区是世 界 上 最 陡 峭 的 大 陆 断 崖 陡 坡 带 的剪切波速度模型切片中均可看出四川盆地的沉积 之一 在水平 跨 度 仅 为 ５０ｋｍ 左 右 的 范 围 内 从 四 这一差异很好地勾绘了四川盆地的西边界 然而 从 川盆地海 拔 约 ４００ ｍ 起 陡 升 至 龙 门 山 及 松 潘 甘 图 ８ 中可知 四川盆 地 与 龙 门 山 脉 及 岷 山 山 脉 的 速 孜地 块 的 ５０００ ｍ 相 关 文 献 Ｃｈｅｎｅ ａ １９９５ 度异常仅存在 于 ８ 至 １４ 的速 度异 常图 相 中 中 岩与周围龙门山 岷山山脉的火成岩存在明显差异

14 5 期 范文渊等 : 青藏高原东部和周边地区地壳速度结构的背景噪声层析成像 1581 长周期的速度异常图像中, 四川盆地与龙门山脉及岷山山脉的速度结构均一稳定. 同样使用背景噪声层析成像的方法, 但是基于不同的地震数据, 李昱在研究这一区域时, 得到相同的结果 (Lietal., 2009). 因此, 根据图 10 中的剪切波速度结构与莫霍面深度图像, 我们认为四川盆地与龙门山山脉的主要差异在于上地壳. 在中下地壳, 青藏高原东缘地区为一稳定均一地质块体, 与青藏高原速度结构差异很大. 由鲜水河小江断裂带 红河断裂带和大理断裂带围起的川滇菱形块体为青藏高原的东南缘主要物质溢出出口. 研究 (Clarketal.,2005;Roydenet al.,1997) 发现该块体内无显著变形 且地表高程变化缓慢, 与四川盆地和龙门山前高程的激烈变化相反.GPS 观测显示该块体连贯一致地绕喜马拉雅东构造结运动 (Zhangetal.,2004). 大地电磁测深成像 (Baietal.,2010) 和面波成像 (Yaoetal., 2008) 研究发现该区域中地壳 (20~40km) 存在层流, 且层流主要存在于该地区主要断裂带之间. 本文的瑞利面波相速度异常 ( 图 8) 结果显示在短周期, 即上地壳, 这一区域呈现高速异常 ; 而在中周期, 即中地壳, 这一区域则呈现低速异常, 并且与青藏高原大范围的低速异常的连接被红河断裂带所割断, 进而从地震学观测的角度验证了近期大地电磁测深成像 (Baietal.,2010) 的发现. 5 结论 利用地震背景噪声记录和互相关技术, 反演获得了青藏高原东部和周边地区的地壳三维地震波速度结构. 瑞利波相速度分布图像结果表明短周期 (6 ~14s) 相速度异常分布与地表地质构造特征吻合较好, 在青藏高原和四川盆地之间存在一个明显的南北向转换带. 中周期 (16~35s) 相速度分布揭示新造高原 柴达木盆地和祁连山脉地区呈现与青藏高原截然不同的中地壳速度结构, 反而与青藏高原东缘地区和川滇菱形块体速度结构相似. 反演获得的剪切波速度在 27.5~45km 深度的切片也明显地揭示 : 青藏高原的松潘 甘孜地块和羌塘地块呈现均一的低速层 ; 然而, 柴达木盆地和祁连山脉地区则呈现较强的横向不均匀性, 尤其是柴达木盆地的高速异常和四川盆地的高速异常相对应. 这些地震学观测证据支持前人提出的青藏高原东北向台阶式增长模式. 与全球一维平均速度模型 (AK135) 相比较发现, 不考虑地壳厚度的显著不同, 本文测量和反演获得的研究区域内平均相速度和剪切波速度都比 AK135 模型慢很多, 尤其是青藏高原的中地壳 (25~ 40km) 剪切波速度显著低于全球平均速度模型, 进一步的层析成像反演证实松潘 甘孜和羌塘地块中地壳 (27.5~45km) 呈现大范围均一的低速层, 为青藏高原可能存在大规模中下地壳 层流 (Channel flow) 提供地震学观测证据. 最后, 在祁连山脉的 27.5~45km 深度观测到的明显低速异常体可能对应于该造山带下地幔岩浆活动导致的底侵作用 (Underplating), 表明引起该地区地壳增厚的主要机制是来自地幔岩浆的底侵作用. 根据这些地震学观测证据, 我们推测现今柴达木盆地和祁连山脉地区的构造运动和隆升过程也就是青藏高原早期的形成模式. 因此深入开展该地区的研究对于认识青藏高原早期形成和隆升过程具有重要科学意义. 致谢 作者向参与北京大学流动地震台阵的架设 数据采集工作的北京大学 地震大地构造学研究小组 的同学和老师表示衷心的感谢. 感谢中国地震局地球物理研究所 国家数字测震台网数据备份中心 为本研究提供地震波形数据. 犚犲犳犲狉犲狀犮犲狊 Bai D H,Unsworth M J, Meju M A,et al Crustal deformation of the eastern Tibetan plateau revealed by magnetoteluricimaging. 犖犪狋狌狉犲犌犲狅狊犮犻犲狀犮犲,3(5): Barmin M P,Ritzwoler M H,Levshin A L.2001.A fastand reliable method for surface wave tomography. 犘狌狉犲犪狀犱犃狆狆犾犻犲犱犌犲狅狆犺狔狊犻犮狊,158(8): BarronJ,Priestley K.2009.Observationsoffrequency dependent SnpropagationinNorthernTibet. 犌犲狅狆犺狔狊犻犮犪犾犑狅狌狉狀犪犾犐狀狋犲狉狀犪狋犻狅狀犪犾, 179(1): BensenG D,RitzwolerM H,BarminM P,etal.2007.Processing seismicambientnoisedatatoobtainreliablebroad bandsurface wavedispersion measurements. 犌犲狅狆犺狔狊犻犮犪犾犑狅狌狉狀犪犾犐狀狋犲狉狀犪狋犻狅狀犪犾, 169(3): Bilham R,LarsonK,FreymuelerJ,etal.1997.GPSmeasurementsof present dayconvergenceacrossthe Nepal Himalaya. 犖犪狋狌狉犲, 386(6620): BourjotL,RomanowiczB.1992.Crustanduppermantletomographyin tibetusingsurface waves. 犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺犔犲狋犲狉狊,19 (9): BrandonC,RomanowiczB.1986.A no lid zoneinthecentral Chang ThangplatformofTibet:Evidencefrompurepathphase

15 1582 地球物理学报 (ChineseJ.Geophys.) 58 卷 velocitymeasurementsoflongperiodrayleigh waves. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 犛狅犾犻犱犈犪狉狋犺犪狀犱犘犾犪狀犲狋狊,91(B6): BurchfielBC,ZhangPZ,WangYP,etal.1991.Geologyofthe Haiyuanfaultzone,Ningxia HuiAutonomousRegion,China, anditsrelationtotheevolutionofthenortheastern marginof thetibetanplateau. 犜犲犮狋狅狀犻犮狊,10(6): CampiloM,PaulA.2003.Long rangecorrelationsinthedifuse seismiccoda. 犛犮犻犲狀犮犲,299(5606): Campilo M.2006.Phaseandcorrelationin Random seismicfields andthereconstructionofthegreenfunction. 犘狌狉犲犪狀犱犃狆狆犾犻犲犱犌犲狅狆犺狔狊犻犮狊,163(2 3): ChenSF,WilsonCJL,DengQD,etal.1994.Activefaultingand block movementassociated withlargeearthquakesinthe Min ShanandLongmen Mountains,northeastern TibetanPlateau. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺,99(B12): ChenSF,WilsonCJL,WorleyB A.1995.Tectonictransition fromthesongpan GarzêfoldbelttotheSichuanBasin,south westernchina. 犅犪狊犻狀犚犲狊犲犪狉犮犺,7(3): ChenSF,WilsonCJL.1996.EmplacementoftheLongmenShan Thrust Nappebeltalongthe eastern margin ofthe Tibetan plateau. 犑狅狌狉狀犪犾狅犳犛狋狉狌犮狋狌狉犪犾犌犲狅犾狅犵狔,18(4): Clark M K,BushJW M,RoydenLH.2005.Dynamictopography produced bylowercrustalflow againstrheologicalstrength heterogeneities bordering the Tibetan Plateau. 犌犲狅狆犺狔狊犻犮犪犾犑狅狌狉狀犪犾犐狀狋犲狉狀犪狋犻狅狀犪犾,162(2): CurtisA,TrampertJ,SniederR,etal.1998.Eurasianfundamentalmode surfacewavephasevelocitiesandtheirrelationshipwithtectonic structures. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 犛狅犾犻犱犈犪狉狋犺,103 (B11): DeCelesP G,Robinson D M,Zandt G.2002.Implicationsof shorteninginthe Himalayanfold thrustbeltforupliftofthe TibetanPlateau. 犜犲犮狋狅狀犻犮狊,21(6): EnglandP,SearleM.1986.Thecretaceous tertiarydeformationof thelhasablockanditsimplicationsforcrustalthickeningin Tibet. 犜犲犮狋狅狀犻犮狊,5(1):1 14. FengC C,Teng T L.1983.Three dimensionalcrustandupper mantlestructureoftheeurasiancontinent. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺,88(B3): GriotD A,MontagnerJP,TapponnierP.1998.Phasevelocity structurefrom Rayleigh and Love wavesin Tibet and its neighboringregions. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 犛狅犾犻犱犈犪狉狋犺,103(B9): GuSS,DiHS.1989.MechanismofformationoftheQaidambasin anditscontrolonpetroleum.//zhuxed.chinesesedimentary Basins.Amsterdam:Elsevier, HarrisonT M,YinA,GroveM,etal.2000.TheZedongWindow: Arecordofsuperposed Tertiaryconvergenceinsoutheastern Tibet. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 犛狅犾犻犱犈犪狉狋犺,105 (B8): HearnT M,WangS Y,NiJF,etal.2004.Uppermostmantle velocitiesbeneath Chinaandsurroundingregions. 犑狅狌狉狀犪犾狅犳 犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 犛狅犾犻犱犈犪狉狋犺,109(B11),doi: / 2003JB HuangZX,SuW,PengYJ,etal.2003.Rayleighwavetomographyof Chinaandadjacentregions. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 犛狅犾犻犱犈犪狉狋犺,108(B2),doi: /2001JB HuangJL,ZhaoDP,ZhengSH.2002.Lithosphericstructureand itsrelationshiptoseismicand volcanicactivityin southwest China. 犑. 犌犲狅狆犺狔狊. 犚犲狊.,107(B10):ESE13 1 ESE HuangJL,ZhaoDP.2006.High resolutionmantletomographyof China and surrounding regions. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 犛狅犾犻犱犈犪狉狋犺,111(B9),doi: /2005JB Jiang M M,ZhouS Y,SandvolE,etal Dlithospheric structurebeneathsouthern Tibetfrom Rayleigh wavetomography witha2 Dseismicaray. 犌犲狅狆犺狔狊犻犮犪犾犑狅狌狉狀犪犾犐狀狋犲狉狀犪狋犻狅狀犪犾,185(2): KangTS,ShinJS.2006.Surface wavetomographyfromambient seismicnoise ofaccelerograph networksin southern Korea. 犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺犔犲狋犲狉狊,33(17),doi: /2006GL LiC,vanderHilstR D,MeltzerAS,etal.2008.Subductionof theindianlithospherebeneaththetibetanplateauandburma. 犈犪狉狋犺犪狀犱犘犾犪狀犲狋犪狉狔犛犮犻犲狀犮犲犔犲狋犲狉狊,274(1 2): LiH Y,Su W,WangC Y,etal.2009.AmbientnoiseRayleigh wavetomographyinwesternsichuanandeasterntibet. 犈犪狉狋犺犪狀犱犘犾犪狀犲狋犪狉狔犛犮犻犲狀犮犲犔犲狋犲狉狊,282(1 4): LiangCT,SongXD,HuangJL.2004.Tomographicinversionof Pntraveltimesin China. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 犛狅犾犻犱犈犪狉狋犺,109(B11),doi: /2003JB LiangCT,SongX D.2006.Alowvelocitybeltbeneathnorthern andeasterntibetanplateaufrom Pntomography. 犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺犔犲狋犲狉狊,33(22),doi: /2006GL LobkisOI,WeaverR L.2001.OntheemergenceoftheGreen's functioninthecorrelationsofadifusefield. 犑狅狌狉狀犪犾狅犳狋犺犲犃犮狅狌狊狋犻犮犪犾犛狅犮犻犲狋狔狅犳犃犿犲狉犻犮犪,110(6): Long D X.1983.A preliminary study on the songpan pingwu earthquakesandtheirstressfield. 犑狅狌狉狀犪犾狅犳犛犲犻狊犿狅犾狅犵犻犮犪犾犚犲狊犲犪狉犮犺,(S1). MétivierF,Gaudemer Y,Tapponnier P,etal.1998.Northeastward growthofthetibetplateaudeducedfrombalancedreconstructionof twodepositionalareas:theqaidamandhexicoridorbasins,china. 犜犲犮狋狅狀犻犮狊,17(6): MeyerB,TapponnierP,BourjotL,etal.1998.Crustalthickening ingansu Qinghai,lithosphericmantlesubduction,andoblique, strike slipcontroledgrowthofthetibetplateau. 犌犲狅狆犺狔狊犻犮犪犾犑狅狌狉狀犪犾犐狀狋犲狉狀犪狋犻狅狀犪犾,135(1):1 47. Pares J M, Van der Voo R, Downs W R, et al Northeastward growth and uplift of the Tibetan Plateau: MagnetostratigraphicinsightsfromtheGuideBasin. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 : 犛狅犾犻犱犈犪狉狋犺 ( ),108(B1): EPM1 1 EPM1 11. PaulA,CampiloM,MargerinL,etal.2005.Empiricalsynthesis oftime asymmetricalgreenfunctionsfrom thecorrelation of codawaves. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 犛狅犾犻犱犈犪狉狋犺,110

16 5 期 范文渊等 : 青藏高原东部和周边地区地壳速度结构的背景噪声层析成像 1583 (B8),doi: /2004JB RitzwolerM H,LevshinAL.1998.Eurasiansurfacewavetomography: Groupvelocities. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 犛狅犾犻犱犈犪狉狋犺,103 (B3): Ritzwoler M H,Levshin A L,Ratnikova L I,etal Intermediate period group velocity mapsacross Central Asia, western China and parts ofthe Middle East. 犌犲狅狆犺狔狊犻犮犪犾犑狅狌狉狀犪犾犐狀狋犲狉狀犪狋犻狅狀犪犾,134(2): RomanowiczB A.1982.ConstraintsonthestructureoftheTibet Plateaufrom purepathphasevelocitiesofloveand Rayleigh waves. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺,87(B8): RoydenL H,BurchfielB C,King R W,etal.1997.Surface deformationandlowercrustalflowineastern Tibet. 犛犮犻犲狀犮犲, 276(5313): Saito M.1988.DISPER80:Asubroutinepackageforthecalculation ofseismicnormalmodesolutions.//doornbosdjed.seismological Algorithms:Computational Methodsand ComputerPrograms.San Diego:AcademicPres, ShapiroN M,Campilo M,StehlyL,etal.2005.High resolution surface wavetomographyfromambientseismicnoise. 犛犮犻犲狀犮犲, 307(5715): ShapiroN M,CampiloM.2004.EmergenceofbroadbandRayleigh wavesfrom corelationsoftheambientseismicnoise. 犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺犔犲狋犲狉狊,31(7),doi: /2004GL Snieder R Extracting the Green s function from the correlation ofcoda waves:a derivation based onstationary phase. 犘犺狔狊犻犮犪犾犚犲狏犻犲狑犈,69(4): SongS G,Zhang L F,Niu Y,etal.2004.Northern Tibetan Paleozoic bloc evolution and continent subduction. 犌犲狅犾狅犵狔犅狉犻犲犳 (inchinese),23(009): Takeuchi H,Dorman J,Saito M.1964.Partialderivatives of surfacewavephasevelocity withrespecttophysicalparameter changeswithintheearth. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺,69 (16): TapponnierP,MeyerB,AvouacJP,etal.1990.ActiveThrusting andfoldingintheqilianshan,anddecouplingbetweenupper crustand mantlein Northeastern Tibet. 犈犪狉狋犺犪狀犱犘犾犪狀犲狋犪狉狔犛犮犻犲狀犮犲犔犲狋犲狉狊,97(3 4): , TapponnierP,PeltzerG,Ledain A Y,etal.1982.Propagating extrusion tectonics in Asia: New insights from simple experimentswithplasticine. 犌犲狅犾狅犵狔,10(12): TapponnierP,XuZQ,RogerF,etal.2001.Obliquestepwiserise andgrowthofthe Tibetplateau. 犛犮犻犲狀犮犲,294(5547): WanJL,Zheng W J,ZhengD W,etal.2010.Lowtemperature thermalevidenceofnorthern Qilian mountaintectonicsatthe latecenozoicera. 犌犲狅犮犺犲犿犻狊狋狉狔 (inchinese),39(5): WangJ,BaoC,LouZ,etal.1989.Formationanddevelopmentof thesichuan Basin.//ZhuJed.ChineseSedimentary Basins. Amsterdam:Elsevier,1: WapenaarK.2004.RetrievingtheelastodynamicGreen'sfunctionof an arbitrary inhomogeneous medium by cross correlation. 犘犺狔狊犻犮犪犾犚犲狏犻犲狑犔犲狋犲狉狊,93(25): WuFT,LevshinA.1994.Surface wavegroupvelocitytomography ofeastasia. 犘犺狔狊犻犮狊狅犳狋犺犲犈犪狉狋犺犪狀犱犘犾犪狀犲狋犪狉狔犐狀狋犲狉犻狅狉狊,84 (1 4): XuG M,Yao H J,Zhu L B,etal.2007.Shear wavevelocity structureofthecrustanduppermantleinwesternchinaandits adjacentarea. 犆犺犻狀犲狊犲犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮狊,50(1): , doi: /cjg XuLL,RondenayS,vanderHilstR D.2007.Structureofthe crustbeneaththesoutheasterntibetanplateaufromteleseismic receiver functions. 犘犺狔狊犻犮狊狅犳狋犺犲犈犪狉狋犺犪狀犱犘犾犪狀犲狋犪狉狔犐狀狋犲狉犻狅狉狊,165(3 4): XuY,Liu F T,LiuJ H,etal.2002.Crustandupper mantle structure beneath western China from P wave traveltime tomography. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 犛狅犾犻犱犈犪狉狋犺, 107(B10):ESE4 1 ESE4 15. YangYJ,Ritzwoler M H,Levshin A L,etal.2007.Ambient noiserayleigh wavetomographyacross Europe. 犌犲狅狆犺狔狊犻犮犪犾犑狅狌狉狀犪犾犐狀狋犲狉狀犪狋犻狅狀犪犾,168(1): YaoHJ,BegheinC,vanderHilstR D.2008.Surface wavearray tomographyinse Tibetfrom ambientseismicnoiseandtwo stationanalysis I.Crustalandupper mantlestructure. 犌犲狅狆犺狔狊犻犮犪犾犑狅狌狉狀犪犾犐狀狋犲狉狀犪狋犻狅狀犪犾,173(1): YueH,ChenYJ,SandvolE,etal.2012.Lithosphericandupper mantlestructureofthenortheasterntibetanplateau. 犑狅狌狉狀犪犾狅犳犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺 : 犛狅犾犻犱犈犪狉狋犺 ( ),117 (B5),doi: /2011JB ZhangPZ,Deng Q D,Zhang G M,etal.2003.Seismicityand seismicactiveblocksofchina. 犆犺犻狀犪犛犮犻犲狀犮犲 ( 犇 )(inchinese), 33(B04): ZhangPZ,ShenZK,WangM,etal.2004.Continuousdeformationof thetibetanplateaufromglobalpositioningsystemdata. 犌犲狅犾狅犵狔,32 (9): ZhengS H,Sun X L,Song X D,etal.2008.Surface wave tomographyofchinafrom ambientseismicnoisecorrelation. 犌犲狅犮犺犲犿犻狊狋狉狔犌犲狅狆犺狔狊犻犮狊犌犲狅狊狔狊狋犲犿狊,9(5),doi: / 2008GC ZhengXF,YaoZX,LiangJH,etal.2010.Theroleplayedand opportunitiesprovidedbyigp DMCofChinanationalseismic networkinwenchuanearthquakedisasterreliefandresearches. 犅狌犾犲狋犻狀狅犳狋犺犲犛犲犻狊犿狅犾狅犵犻犮犪犾犛狅犮犻犲狋狔狅犳犃犿犲狉犻犮犪,100(5B): 附中文参考文献 宋述光, 张立飞,NiuY, 等 青藏高原北缘早古生代板块构造演化和大陆深俯冲. 地质通报,23(9): 万景林, 郑文俊, 郑德文等 祁连山北缘晚新生代构造活动的低温热年代学证据. 地球化学,39(5): 张培震, 邓起东, 张国民等 中国大陆的强震活动与活动地块. 中国科学 :D 辑,33(B04): ( 本文编辑汪海英 )