Structure of lithospheric discontinuities beneath the southeast margin of Tibetan Plateau revealed by seismic daylight imaging and its dynamic significance
-
摘要: 为了进一步认识青藏高原东南缘的构造演化等动力学问题以及该区域的深部孕震机理，本文使用位于该区域内的中国地震科学台阵探测项目的台站所记录到的远震P波波形数据，采用地震光照成像法获取了岩石圈间断面的结构，并讨论了该方法的准确性和稳定性。研究结果显示，青藏高原东南缘的岩石圈西薄东厚，其中：滇缅泰地块腾冲火山附近最薄，约为60 km，其较薄的岩石圈可能是软流圈地幔物质上涌造成的；扬子地块岩石圈厚度从四川盆地向南逐渐减薄，特别是四川盆地下方最厚，可达190 km左右；滇缅泰地块腾冲火山下方150 km深度左右探测到明显的间断面，该间断面可能是腾冲火山原始岩浆源的位置即岩浆源。本研究所得结果 “印支地块与滇缅泰地块结构的连续性” 进一步为印度板块的推挤作用造成腾冲火山低速物质向东溢出的结论提供了地震学证据。此外，研究区域最北端的剖面显示，峨眉山大火成岩省的内带在50—250 km深度范围及其上方地壳内存在明显的局部高速异常，其不均匀分布特征可能与二叠纪火山喷发过程中岩浆底侵及中新生代以来多期次构造活动有关。Abstract: In order to further understand the dynamic problems such as tectonic evolution and deep seismogenic mechanism of the southeastern margin of the Tibetan Plateau, the structure of the lithospheric discontinuity was obtained by seismic daylight imaging based on teleseismic P-wave waveform data recorded by stations of China Seismological Science Array Exploration Project located in the region, and the accuracy and stability of the method were also discussed. The results show that the lithospheric thickness in the southeastern margin of the Tibetan Plateau is thinner in the west and thicker in the east. The thinnest lithospheric thickness near Tengchong volcano in Yunnan-Myanmar-Thailand block is about 60 km. The thinner lithosphere may be caused by the upwelling of asthenosphere mantle material. The lithospheric thickness of the Yangtze block is gradually thinning from Sichuan basin to the south, especially the thickest beneath the Sichuan basin, whose lithospheric thickness can reach about 190 km. Meanwhile, obvious discontinuities have been detected at a depth of about 150 km beneath Tengchong volcano, which may be the original magma source of volcano. The continuity of the structure of the Indo-China block and the Yunnan-Myanmar-Thailand block further provides seismological evidence for the eastward spillover of low-velocity materials from Tengchong volcano caused by the pushing of the Indian Plate. The northernmost section of the studied area shows that there are obvious local high-velocity variations in the depth range of 50−250 km within the Emeishan Large Igneous Province and the crust above it. The uneven distribution characteristics may be related to the magma bottom intrusion during the Permian volcanic eruption and the multiple tectonic activities since the Mesozoic and Cenozoic.
-
引言
印度−欧亚板块的陆−陆碰撞造山运动造就了世界上海拔最高、面积最大的高原—青藏高原，喜马拉雅造山带作为印度板块与欧亚板块相互碰撞的前沿地带，其造山作用一直为地球科学家所关注（Tapponnier et al，1982；Yin，Harrison，2000；许志琴等，2008）。喜马拉雅东构造结处于喜马拉雅造山带东侧（图1），是地壳大规模缩短和构造运动方向发生转变的轴心区域，山脉和水系的走向发生了90°大转弯，地球物理观测结果也显示出相似的剧烈变形，例如GPS测量的地壳运动速度场和上地幔各向异性远震SKS波分裂快波方向都围绕喜马拉雅东构造结发生了顺时针90°旋转（Gan et al，2007；Sol et al，2007；Chang et al，2015）。该地区主要发育有NNE−NE向和WNW−NW向两组活动断裂，例如NNE−NE的米林断裂和墨脱断裂，WNW−NW向的嘉黎断裂和阿帕龙断裂，且地震活动与断裂构造密切相关（邵翠茹等，2008；王晓楠等，2018）。喜马拉雅东构造结的地形变化剧烈，深部结构复杂，强震频发，曾在1950年发生过察隅M8.6大地震，是开展地球动力学研究的天然实验场所。
![]()
图 1 喜马拉雅东构造结区域构造（邓起东等，2002）及台站和2016年近震的分布F1：雅鲁藏布江断裂；F2：墨脱断裂；F3：嘉黎断裂；F4：迫龙—旁辛断裂；F5：阿帕龙断裂；F6：主边界断裂；F7：墨竹工卡断裂；F8：班公错—怒江断裂。左上角为GPS计算的青藏高原地壳运动速度场（Gan et al，2007），其中黑框为研究区Figure 1. The regional geologic setting （Deng et al，2003） and distribution of local earthquakes in 2016 andstations in the eastern Himalayan syntaxisF1：Yarlung Zangbo River fault；F2：Motuo fault；F3：Jiali fault；F4：Polong-Pangxin fault；F5：Apalon fault；F6：Main boundary fault；F7：Mozhugongka fault；F8：Bangongcuo-Nujiang fault. Inset shows the surface velocity field determined from GPS observations in the Tibetan Plateau （Gan et al，2007），and black frame is the studied region地震波在各向异性介质中传播时，其传播速度和偏振方向及其它特性随波的传播方向的变化而变化，这种现象称为地震各向异性，可以用于探测地壳变形。横波分裂是地震波在各向异性介质中最明显的表现形式，有两个分裂参数，即快波偏振方向和快、慢波的到时延迟，分别反映了介质的变形方向和强度（Silver，1996）。横波分裂被广泛应用于大陆壳幔各向异性变形特征研究中，较深的下地壳和地幔的各向异性被认为是由于矿物晶格优势排列所引起（Mainprice，Nicolas，1989；Silver，1996）；浅部上地壳的各向异性一般被认为是由大量微裂隙在应力作用下定向排列而形成（Crampin，Atkinson，1985）。深部各向异性特征分析一般采用远震横波分裂，在喜马拉雅东构造结及周边地区，Sol等（2007）和Chang等（2015）利用远震SKS波分裂开展了上地幔各向异性特征分析，其结果显示快波方向呈现出绕研究区顺时针旋转的特征，Chang等（2015）基于SKS波分裂、GPS和断裂第四纪滑动速率联合分析了岩石圈变形特征，动力学模拟结果显示喜马拉雅东构造结周边地区的岩石圈变形特征为垂直连贯变形模式，即力学上的壳幔耦合模式。浅部上地壳各向异性特征一般采用近震横波分裂，近震横波分裂特征揭示出研究区的断裂构造和应力场控制的微裂隙结构，经常被用来评估区域应力场分布和断裂的活动特性（Boness，Zoback，2006；高原等，2018）。近震横波分裂研究曾在国内多地开展，例如卢龙地区（姚陈等，1992）、唐山地区（高原等，1999）、汶川地震震源区（丁志峰等，2008）和青藏高原东北缘（郭桂红等，2015）等地区。尽管喜马拉雅东构造结地区的地震活动性较强，但是因为近震横波分裂工作需要近震分布区内有台站分布，而研究区野外工作条件恶劣，开展地震观测相对其它地区困难，因此，此项工作一直未开展。本文拟利用2016年雅鲁藏布江下游布设的地震台阵波形记录开展近震横波分裂研究，分析喜马拉雅东构造结上地壳各向异性特征，以期得到该特征与活动断裂和应力场特征的相关性，以及喜马拉雅东构造结与板块碰撞相关的动力学特征。
1. 资料和方法
雅鲁藏布江下游台阵由16个宽频带流动地震台站组成（图1），这些流动地震台站均统一应用REFTEK数据采集器和CMG-3ESP地震计（频带范围为50 Hz−60 s）。台站代码分别是BAX，BEB，DAM，DEX，DOJ，LAD，NYG，QID，SAM，SCY，TOM，WOL，YIG，ZIB，ZOB和ZOS。从图1可以看出，研究所用的地震台站较好地覆盖了地震分布区域，为开展近震横波分裂提供了较为丰富的地震观测资料。本文分析采用记录较为完备的2016年近震波形数据。
在计算地震各向异性参数时，要求地震记录必须清晰记录在横波窗内，横波窗口即S波到达台站时，以sin−1（vP/vS）为顶角的入射角锥形体。在横波入射到地表时，如果入射角大于临界角就会发生全反射现象，造成波形畸变，为避免锥形体内的S波记录被S－P转换波干扰，一般情况下，介质泊松比为0.25时，横波入射角窗口大约为35°。在实际分析中，由于浅部低速层的存在所造成的射线弯曲，可以适度将横波窗的入射角扩至45°—50° （Crampin，Peacock，2005）。因此，此次选取入射角≤45°的波形资料，可使所有波形记录满足窗口要求。
横波分裂参数测量有可视化测量（Chen et al，1987；Peacock et al，1988；Liu et al，1997）和计算机程序自动测定（Shih et al，1989）两种方法。可视化测量方法的步骤为：首先判断选取地震射线能否在横波窗口内截取合适的横波波形；其次，根据质点运动轨迹判定横波是否发生分裂和分析快波偏振方向；最后，测量快、慢波的到时延迟以及检验等工作。该方法每一步均经过人工检测，相较于计算机程序自动测定，该方法更准确可靠，且所提取的地震波形数据更全面，因此本文采用该方法进行测量。
以BAX台记录的一次近震为例（图2），介绍在三分量地震波形记录中进行横波分析分裂的具体过程（丁志峰等，2008；常利军等，2010）：
![]()
图 2 BAX台站的一次近震波形记录的横波分裂分析示例（a） 原始横波的三分量波形记录；（b） 图（a）中P波初动窗口放大后的波形记录；（c） 两水平分量的波形记录；（d） 两水平分量的横波质点运动；（e） 旋转到快、慢波方向的水平分量波形；（f） 去除各向异性校正后快、慢波水平分量；（g） 校正后的快、慢波水平分量的质点运动图；（h）用于横波分裂分析的近震震源深度分布Figure 2. An example of shear wave splitting analysis for a recording at station BAX（a） Three-component records of original seismic waveform；（b） Amplified three-component records of P wave in Fig.（a）；（c） Two horizontal components of seismograms；（d） Shear wave particle motion in horizontal；（e） Seismograms rotated to the fast and slow shear wave directions；（f） Corrected seismograms in fast-slow coordinate system；（g） Partical motion of the corrected phases；（h） The focal depth distribution of the earthquakes used in shear wave splitting1） 选取清晰可靠的初动信号，测量P波初动的水平（EW和NS）分量和垂直（UD）分量，并计算其入射角和入射方向；同时，要求P波初动的垂直分量应大于水平分量两倍以上，以确保所有地震波都被记录在横波窗内（图2a，b）。
2） 根据两个水平（EW和NS）分量的波形情况选取合适的横波分裂窗口，一般时间间隔在0.3—0.5 s之间，需选取高信噪比地震波形记录，较高的信噪比对后续步骤的顺利进行和结果的准确性影响很大（图2c）。
3） 根据两水平分量的质点运动图是否为椭圆来判断是否存在横波分裂现象。如图2d所示，质点运动图为一个近似椭圆，且从S波初动的质点运动轨迹可以清晰地看出快波偏振方向为近NS方向，慢波方向为与其垂直的近EW方向，依此可判断其为一个合适的近震横波分裂例子；适当旋转角度，可测得快波偏振方向。
4） 将两个水平分量的波形旋转至快、慢波方向（图2e），能够清晰地观察到两个到时不同但波形相似的记录，就此直接测量慢波延迟时间，即快、慢波的到时差，且各向异性校正后的快、慢波分量的到时差消失（图2f），其质点运动轨迹变为近似直线（图2h）。然后根据台站到震源的射线路径距离到时延迟进行归一化处理。
2. 结果
本文对喜马拉雅东构造结区域内16个地震台站的近震波形记录开展横波分裂测量，采用可视化测量方法，从大量的近震横波记录中选取符合横波分裂窗口的记录数据，严格按照上一节所述计算步骤和检验方法，最后得到了这16个台站的横波分裂测量参数，慢波延迟时间要经过归一化处理，在研究区各台站总计得到了369个有效分裂结果。研究区各台站基本参数和横波分裂参数列于表1。图3为研究区各台站快波偏振方向的等面积投影玫瑰图分布。对照表1和图3可以看出，一般在地震密集分布区域内的台站可用于近震横波分裂的有效记录个数较多，比如台站BAX 31个，DAM 12个，DEX 21个，DOJ 54个，LAD 27个，QID 26个，SAM 59个，TOM 27个，ZIB 60个和ZOS 29个等，近震横波分裂的有效个数除了DAM台12个外，其余台站都在20个以上；一般距离地震密集分布区较远的台站可用于近震横波分裂的有效记录个数较少，例如台站BEB 2个，NYG 3个，SCY 7个，WOL 3个，YIG 4个和ZOB 4个等，近震横波分裂有效记录个数都在10个以内。从图3所示的各台站快波偏振优势方向来看，总体上各台站的快波偏振优势方向明显，均呈现了一个主要的优势方向，表1中各台站快波偏振方向的平均值基本上与图3中相应台站的快波偏振优势方向一致，各台站得到的慢波到时延迟平均值为0.99—2.26 ms/km。
表 1 各台站基本参数与横波分裂参数Table 1. Station parameters and results of shear-wave splitting台站代码 东经/° 北纬/° 快波偏振
方向/°快波偏振方向
标准差/°到时延迟
/（ms·km−1）到时延迟标准差
/（ms·km−1）有效记
录条数BAX 95.4 29.6 177 17 2.04 0.65 31 BEB 95.2 29.2 41 4 1.35 0.60 2 DAM 95.5 29.5 165 29 1.83 0.50 12 DEX 95.3 29.3 26 33 1.91 0.92 21 DOJ 94.8 30.0 55 25 1.57 0.68 54 LAD 93.1 29.0 1 12 2.26 0.49 27 NYG 94.2 29.1 60 9 1.07 0.39 3 QID 95.6 30.1 108 9 1.67 0.47 26 SAM 97.0 28.4 165 22 1.87 0.47 59 SCY 96.7 28.8 156 32 0.99 0.25 7 TOM 95.1 30.1 18 13 1.36 0.40 27 WOL 93.7 29.1 96 10 1.48 0.13 3 YIG 94.8 30.2 12 17 1.68 0.66 4 ZIB 94.9 29.6 38 16 0.97 0.28 60 ZOB 96.3 29.6 122 13 2.23 0.46 4 ZOS 93.3 29.8 89 21 1.61 0.32 29 ![]()
图 3 喜马拉雅东构造结各台站快波偏振方向等面积投影玫瑰图分布Figure 3. The distribution of the homolographic projection rose diagrams of the fast directions at the stations in the eastern Himalayan syntaxis3. 分析和讨论
程成等（2017）和Wang等（2019）利用P波接收函数方法研究了喜马拉雅东构造结的地壳结构，其结果显示地壳厚度从SW向NE方向逐渐增厚的趋势，大约从SW侧的55 km过渡到NE侧的70 km，其中Wang等（2019）综合布格重力异常剖面分析得到上地壳埋深为20—25 km。该区以往地震活动性研究显示地震主要发生在浅部上地壳（杨建亚等，2017；Bai et al，2017；王晓楠等，2018）。如图2h所示，本文用于横波分裂的近震事件震源深度主要分布在8—22 km范围内，因此，本文所得到的近震横波分裂特征主要反映了研究区上地壳各向异性的变形特征。
上地壳介质中广泛分布的微裂隙在应力的作用下会沿一定优势方向排列，这是导致上地壳各向异性的主要原因（Crampin，Atkinson，1985）。通常横波在上地壳各向异性介质中产生的分裂特征可以反映区域应力场特征和活动断裂特征（Crampin，Peacock，2005；高原等，2018）。Boness和Zoback （2006）归纳出上地壳各向异性主要为应力作用和断裂结构所致，即应力和结构控制的各向异性。上地壳近震横波分裂的快波偏振优势方向与微裂隙的定向排列一致，与原地最大主压应力方向一致；靠近或位于活动断裂带上的台站，其快波偏振优势方向通常与活动断裂的走向密切相关，慢波延迟时间对应于上地壳微裂隙在构造应力作用下的几何形态，反映了构造应力的强度和效应，能够较好地反映当前断裂的活动特征（丁志峰等，2008；吴晶等，2009；高原等，2018；Wu et al，2019）。
研究区各台站得到的快波偏振优势方向自西向东，由LAD台的近NS方向到ZOS台和WOL台的近EW方向，再到台站NYG，ZIB，DOJ，BEB和DEX的NE方向，然后是台站YIG，TOM，BAX和DAM的近NS方向或NNE方向，最后到研究区东部台站QID，ZOB，SCY和SAM的NW方向（图3）。大部分位于或靠近活动断裂的台站，其快波偏振优势方向与断裂的走向相一致：位于墨竹工卡断裂上的ZOS台，快波偏振优势方向与近EW的断裂走向一致；位于雅鲁藏布江断裂西段上的WOL台，其快波偏振优势方向与近EW方向的断裂走向一致；位于雅鲁藏布江断裂东段的台站NYG，ZIB和DOJ和墨脱断裂上的台站BEB和DEX的快波偏振优势方向与NE向的断裂走向一致；位于研究区东部NW走向的嘉黎断裂的台站QID，ZOB，SCY和SAM快波偏振优势方向为NW向，与断裂走向一致；位于迫龙—旁辛断裂的台站BAX和DAM的快波偏振优势方向与此断裂的近NS方向的走向一致。距离雅鲁藏布江断裂西段和东段有一定距离的台站LAD和YIG，以及位于雅鲁藏布江断裂东段与嘉黎断裂交会处的TOM台的快波偏振优势方向与该两条断裂走向存在一定角度，但是其与印度板块持续向欧亚板块的推挤作用下在喜马拉雅东构造结及周边地区形成的NNE向主压应力场方向一致（许忠淮，2001），而且也与由GPS和断裂第四纪滑动速率数据分析得到的NNE向压应变方向一致。由以上分析可以看出，喜马拉雅东构造结各台站的快波偏振优势方向基本上可以归纳为受活动断裂控制和应力控制的两组，具体为13个靠近或位于墨竹工卡断裂、雅鲁藏布江断裂、墨脱断裂和嘉黎断裂的台站的上地壳各向异性主要受断裂结构影响，而另外2个距活动断裂有一定距离的台站LAD和YIG，以及处于断裂交会处的TOM台，其上地壳各向异性主要受应力控制。
上地壳近震横波分裂的参数对区域应力场的变化比较敏感，可以通过横波分裂参数随时间的变化特征来分析区域应力场的变化特征（Gao，Crampin，2003）。前人的研究显示，位于地震震源区的台站其近震横波分裂参数在主震前后变化明显，快波偏振方向在主震后一段时间离散度增大，到时延迟在主震后一段时间会增大，随后快波偏振方向和到时延迟会逐渐恢复到正常水平，例如Crampin等（1999）关于冰岛一次M5地震的前震和余震序列的分析显示其横波分裂参数具有随时间变化的规律特征，国内在汶川MS8.0地震和玉树MS7.1地震的震源区横波分裂参数也表现出随时间的变换规律（丁志峰等，2008；常利军等，2010）。为了分析喜马拉雅东构造结横波分裂参数和区域应力场的特征，绘制了2016年各台站横波分裂参数结果随时间的变化图（图4），该图显示快波偏振方向和到时延迟并未表现出随时间的规律变化，整体上比较稳定。根据中国地震台网中心（2016）目录，研究区最大地震为3月21日0点54分发生于工布江达县的ML3.8地震，震中位置为（29.79°N，93.16°E），震源深度为7 km。分析显示，2016年研究区各台站的横波分裂参数未表现出随时间的规律变化特征，主要是由于2016年研究区的地震强度较低，尽管小震定位研究区发生了数百次小震（图1），但未发生过M＞5.0较强地震，研究区由小震导致的应力调整还不足以引起横波分裂参数的规律变化。尽管快波偏振方向和到时延迟未表现出随时间的规律变化，而快波偏振优势方向比较突出，但是各台站快波偏振方向和到时延迟时间的离散度还是比较大（图4和表1）。16个台站中，9个台站的快波偏振方向标准差超过15°，15个台站的到时延迟时间标准差超过0.25 ms/km，而且台站间的横波分裂参数变化也较大，快波偏振优势方向有近EW，近NS，NE和NW等各个方向，到时延迟时间的平均值范围为0.99—2.26 ms/km。这也体现了在印度板块持续向欧亚板块碰撞过程中形成的复杂应力场导致了喜马拉雅东构造结复杂的构造和剧烈的变形特征。
![]()
图 4 各台站快波偏振方向和快慢波到时延迟随时间的变化图（a—d）Figure 4. Temporal variations of the fast polarization directions and the normalized delay times at the stations （a−d）![]()
图 4 各台站快波偏振方向和快慢波到时延迟随时间的变化图（e—l）Figure 4. Temporal variations of the fast polarization directions and the normalized delay times at the stations （e−l）![]()
图 4 各台站快波偏振方向和快慢波到时延迟随时间的变化图（m—p）Figure 4. Temporal variations of the fast polarization directions and the normalized delay times at the stations （m−p）4. 结论
本文基于雅鲁藏布江下游台阵16个台站2016年的近震记录，使用可视化横波分裂分析方法计算了各台站横波分裂参数，获得了喜马拉雅东构造结的上地壳各向异性特征。整体上，研究区位于或靠近活动断裂带上的台站，其快波偏振优势方向与断裂的走向基本一致，距离活动断裂有一定距离或位于断裂交会处的台站，其快波偏振优势方向与主压应力场方向一致，反映了结构控制和应力控制作用下的上地壳各向异性特征。研究区各台站之间的横波分裂结果呈现较大差异，自身离散度较大，这些特征反映了东构造结复杂的结构和剧烈的变形。由于研究区2016年的地震活动强度不大，横波分裂参数未表现出随时间的规律变化。此项研究需要对其它震例进一步深入研究，如2017年11月17日的米林MS6.9地震，可能会得到随时间的变化特征。
-
![]()
图 1 青藏高原东南缘的构造背景及台站、剖面分布
F1：怒江断裂；F2：红河断裂；F3：小江断裂，下同
Figure 1. Tectonic settings and distribution of stations and profiles in the southeast margin of TibetanPlateau
F1：Nujiang fault；F2：Honghe fault；F3：Xiaojiang fault，the same below
![]()
图 3 多层介质模型（a）及其利用地震光照成像法构建的深度域地震反射波（b）
图（b）中stk道集为入射慢度为0.04—0.08 s/km时对应的远震事件地震反射波叠加的结果
Figure 3. A multi-layer medium model （a） and the reflection wave in depth domain constructed by seismic daylight imaging method （b）
In Fig. （b） seismic gathers labeled stk is the result of stacked reflected waves from teleseismic events with slowness 0.04−0.08 s/km
![]()
图 6 台站53080 （a），53107 （b），53034 （c）和4501 （d）上地震光照成像法稳定性分析结果
每幅图中左侧黑色实线为叠加自相关道集，灰色实线为波形包络，灰色阴影区为波形包络的绝对值
Figure 6. The stability analyses on the seismic daylight imaging for the stations 53080 （a），53107 （b），53034 （c） and 4501 （d）
The black line on the left side of each image is the stacked autocorrelograms gather，the grey line is the seismic envelope，and the grey shadow area is the absolute value of the seismic envelope
![]()
图 4 台站53085 （a）和53162 （b）的莫霍面结果对比
左侧黑色实线为叠加自相关道集，浅灰色实线为波形包络，灰色阴影区为波形包络的绝对值；深红色实线为根据 SDI方法得到的莫霍面位置，绿色虚线为Sun等（2014）的莫霍面结果，右侧的速度曲线为其反演得到的S波速度模型曲线
Figure 4. Comparison of Moho results for the stations 53085 （a） and 53162 （b）
The black solid line on the left is stacked autocorrelograms gather，the light grey solid line is seismic envelope，the grey shadow area is absolute value of seismic envelope；the deep red solid line is Moho position obtained by SDI method，and the green dashed line is Moho surface result from Sun et al （2014）. The velocity curve on the right side is S-wave velocity model curve obtained by inversion
![]()
图 5 剖面A （a），B （b）和C （c）的间断面结构图
黑色实线为叠加自相关道集的波形图，灰色阴影为波形包络，蓝色实线为探测到的较深主间断面的位置，绿色实线为确定的位于莫霍面与主间断面之间的次级间断面的位置，白色实线为地震光照成像法确定的台站下方莫霍面的深度；图中50 km左右的黑色实线为Sun等（2014）的莫霍面的结果。为了凸显50 km以下区域的速度对比，更好地讨论间断面的位置，图中50 km以上的区域为粉红色并且不包含在色标中，4.4 km/s以下的速度也未包含在色标中
Figure 5. The structure of the discontinuities for the profiles A （a），B （b）和C （c）
The black solid line is the waveform of the superimposed autocorrelation gather，the grey shadow is the seismic envelope，the blue solid line is the location of the deep main discontinuity detected，the green solid line is location of the secondary discontinuity located between Moho surface and the main discontinuity，and white solid lines are the Moho depth beneath the stations. The black solid line near 50 km is the Moho result of Sun et al （2014）. In order to highlight the velocity contrast of the area deeper than 50 km，and to better discuss the position of the discontinuities，the area above the depth of 50 km in the figs. is pink and not included in the color mark，and the velocity smaller than 4.4 km/s is not included in the color mark here
![]()
图 7 剖面A （a），B （b）和C （c）的岩石圈-软流圈边界（LAB）结果图
黑色实线为叠加自相关道集的波形图，灰色阴影为波形包络，绿色实线为确定的可能的LAB的位置；图中50 km左右的黑色实线为Sun等（2014）的莫霍面结果。为了凸显50 km以下区域的速度对比，更好地讨论间断面的位置，图中50 km以上的区域为粉红色并且不包含在色标中，4.4 km/s以下的速度未包含于色标中
Figure 7. The results of the LAB for the profiles A （a），B （b） and C （c）
The black line is the waveform of the stacked autocorrelograms，the grey shadow is the seismic envelope，and the green line is the location of the possible LAB. The black solid line near 50 km depth is the Moho result of Sun et al （2014）. In order to highlight the velocity contrast of the area below 50 km，and to better discuss the position of the discontinuity，the area above 50 km in the fig. is pink and not included in the color mark，and the velocity below 4.4 km/s is not included in the color mark
表 1 剖面B台站下方间断面的具体判别准则
Table 1 Specific criteria for discrimination of discontinuities beneath the stations of section B
台站名 间断面1 间断面2 间断面3 间断面4 53080 地震波组振幅、频率的变化
及明显的反射信号地震波组频率的变化 地震波组频率及波形包络 53076 地震波组频率的变化及明
显的反射信号地震波组频率的变化及明显
的反射信号地震波组频率的变化及地震波形包络 53073 地震波组振幅、频率的变化 地震波组振幅、频率的变化 地震波组振幅、频率的变化 53089 地震波组振幅、频率的变化 地震波组振幅、频率的变化 地震波组振幅、频率的变化及波形包络 53151 地震波组频率的变化 地震波组频率的变化 地震波组频率的变化及地震波形包络 53149 地震波组频率的变化 地震波组振幅、频率的变化
及明显的反射信号地震波组振幅、频率的变化 地震波组振幅、
频率的变化53160 地震波组频率的变化 地震波组频率的变化 地震波组振幅、频率的变化及波形包络 53158 地震波组频率的变化及
明显的反射信号地震波组振幅、频率的变化
及明显的反射信号地震波组振幅、频率的变化及波形包络 地震波组频率
的变化53226 地震波组振幅、频率的变化 地震波组频率的变化及明显
的反射信号地震波组频率的变化及明显的反射信号 53218 地震波组频率的变化 地震波组频率的变化及明显
的反射信号地震波组振幅、频率的变化 52036 地震波组振幅、频率的变化 地震波组振幅、频率的变化 地震波组振幅、频率的变化 地震波组振幅、
频率的变化52037 地震波组振幅、频率的变化
及明显的反射信号地震波组频率的变化及明显
的反射信号地震波组振幅、频率的变化 52039 地震波组振幅、频率的变化
及明显的反射信号地震波组振幅、频率的变化 地震波组振幅、频率的变化及明显的反射信号 
下载: 导出CSV
-
高名修. 1996. 青藏高原东南缘现今地球动力学研究[J]. 地震地质,18(2):129–142. Gao M X. 1996. A study on recent geodynamics at the southeastern margin of the Qinghai-Xizang (Tibet) Plateau[J]. Seismology and Geology,18(2):129–142 (in Chinese).
阚荣举,赵晋明,阚丹. 1996. 腾冲火山地热区的构造演化与火山喷发[J]. 地震地磁观测与研究,17(4):28–33. Kan R J,Zhao J M,Kan D. 1996. The tectonic evolution and volcanic eruption in Tengchong volcanogeothermic region[J]. Seismological and Geomagnetic Observation and Research,17(4):28–33 (in Chinese).
李永华,吴庆举,田小波,张瑞青,潘佳铁,曾融生. 2009. 用接收函数方法研究云南及其邻区地壳上地幔结构[J]. 地球物理学报,52(1):67–80. Li Y H,Wu Q J,Tian X B,Zhang R Q,Pan J T,Zeng R S. 2009. Crustal structure in the Yunnan region determined by modeling receiver functions[J]. Chinese Journal of Geophysics,52(1):67–80 (in Chinese).
刘福田,曲克信,吴华,李强,刘建华,胡戈. 1989. 中国大陆及其邻近地区的地震层析成象[J]. 地球物理学报,32(3):281–291. doi: 10.3321/j.issn:0001-5733.1989.03.005 Liu F T,Qu K X,Wu H,Li Q,Liu J H,Hu G. 1989. Seismic tomography of the Chinese continent and adjacent region[J]. Acta Geophysica Sinica,32(3):281–291 (in Chinese).
穆治国,佟伟,Curtis G H. 1987. 腾冲火山活动的时代和岩浆来源问题[J]. 地球物理学报,30(3):261–270. doi: 10.3321/j.issn:0001-5733.1987.03.005 Mu Z G,Tong W,Curtis G H. 1987. Times of volcanic activity and origin of magma in Tengchong geothermal area,west Yunnan Province[J]. Acta Geophysica Sinica,30(3):261–270 (in Chinese).
潘佳铁,李永华,吴庆举,丁志峰. 2015. 青藏高原东南部地区瑞雷波相速度层析成像[J]. 地球物理学报,58(11):3993–4006. Pan J T,Li Y H,Wu Q J,Ding Z F. 2015. Phase velocity maps of Rayleigh waves in the southeast Tibetan Plateau[J]. Chinese Journal of Geophysics,58(11):3993–4006 (in Chinese).
乔学军,王琪,杜瑞林. 2004. 川滇地区活动地块现今地壳形变特征[J]. 地球物理学报,47(5):805–811. doi: 10.3321/j.issn:0001-5733.2004.05.011 Qiao X J,Wang Q,Du R L. 2004. Characteristics of current crustal deformation of active blocks in the Sichuan-Yunnan region[J]. Chinese Journal of Geophysics,47(5):805–811 (in Chinese).
孙若昧,刘福田,刘建华. 1991. 四川地区的地震层析成像[J]. 地球物理学报,34(6):708–719. doi: 10.3321/j.issn:0001-5733.1991.06.005 Sun R W,Liu F T,Liu J H. 1991. Seismic tomography of Sichuan[J]. Acta Geophysica Sinica,34(6):708–719 (in Chinese).
孙伟家,符力耘,魏伟,唐清雅. 2019. 探测岩石圈间断面精细结构的地震光照成像新方法[J]. 中国科学: 地球科学,49(3):521–536. Sun W J,Fu L Y,Wei W,Tang Q Y. 2019. A new seismic daylight imaging method for determining the structure of lithospheric discontinuity[J]. Science China Earth Sciences,49(3):521–636 (in Chinese).
王夫运,潘素珍,刘兰,刘宝峰,张建狮,邓晓果,马策军,张彩军. 2014. 玉溪—临沧剖面宽角地震探测:红河断裂带及滇南地壳结构研究[J]. 地球物理学报,57(10):3247–3258. doi: 10.6038/cjg20141013 Wang F Y,Pan S Z,Liu L,Liu B F,Zhang J S,Deng X G,Ma C J,Zhang C J. 2014. Wide angle seismic exploration of Yuxi-Lincang profile:The research of crustal structure of the Red River fault zone and southern Yunnan[J]. Chinese Journal of Geophysics,57(10):3247–3258 (in Chinese).
王雅明,佟殿君,任建业. 2007. 印支地块的逃逸旋转与莺歌海盆地的发育演化[J]. 断块油气田,14(2):33–35. doi: 10.3969/j.issn.1005-8907.2007.02.011 Wang Y M,Tong D J,Ren J Y. 2007. Extrusion and rotation of Indo-China block and development and evolutional process of Yinggehai Basin,South China Sea[J]. Fault-Block Oil and Gas Field,14(2):33–35 (in Chinese).
吴建平,杨婷,王未来,明跃红,张天中. 2013. 小江断裂带周边地区三维P波速度结构及其构造意义[J]. 地球物理学报,56(7):2257–2267. doi: 10.6038/cjg20130713 Wu J P,Yang T,Wang W L,Ming Y H,Zhang T Z. 2013. Three dimensional P-wave velocity structure around Xiaojiang fault system and its tectonic implications[J]. Chinese Journal of Geophysics,56(7):2257–2267 (in Chinese).
吴庆举,曾融生. 1998. 用宽频带远震接收函数研究青藏高原的地壳结构[J]. 地球物理学报,41(5):669–679. doi: 10.3321/j.issn:0001-5733.1998.05.010 Wu Q J,Zeng R S. 1998. The crustal structure of Qinghai-Xizang Plateau inferred from broadband teleseismic waveform[J]. Acta Geophysica Sinica,41(5):669–679 (in Chinese).
徐强,赵俊猛,崔仲雄,刘明乾. 2009. 利用接收函数研究青藏高原东南缘的地壳上地幔结构[J]. 地球物理学报,52(12):3001–3008. doi: 10.3969/j.issn.0001-5733.2009.12.009 Xu Q,Zhao J M,Cui Z X,Liu M Q. 2009. Structure of the crust and upper mantle beneath the southeastern Tibetan Plateau by P and S receiver functions[J]. Chinese Journal of Geophysics,52(12):3001–3008 (in Chinese).
徐涛,张忠杰,刘宝峰,陈赟,张明辉,田小波,徐义刚,滕吉文. 2015. 峨眉山大火成岩省地壳速度结构与古地幔柱活动遗迹: 来自丽江—清镇宽角地震资料的约束[J]. 中国科学: 地球科学,45(5):561–576. Xu T,Zhang Z J,Liu B F,Chen Y,Zhang M H,Tian X B,Xu Y G,Teng J W. 2015. Crustal velocity structure in the Emeishan Large Igneous Province and evidence of the Permian mantle plume activity[J]. Science China Earth Sciences,58(7):1133–1147. doi: 10.1007/s11430-015-5094-6
徐义刚,钟孙霖. 2001. 峨眉山大火成岩省: 地幔柱活动的证据及其熔融条件[J]. 地球化学,30(1):1–9. doi: 10.3321/j.issn:0379-1726.2001.01.002 Xu Y G,Zhong S L. 2001. The Emeishan Large Igneous Province:Evidence for mantle plume activity and melting conditions[J]. Geochimica,30(1):1–9 (in Chinese).
杨晓涛,胥颐,刘建华,李志伟. 2011. 腾冲火山区的地震层析成像及其构造意义[J]. 地球物理学报,54(8):2050–2059. doi: 10.3969/j.issn.0001-5733.2011.08.012 Yang X T,Xu Y,Liu J H,Li Z W. 2011. Seismic tomography in the Tengchong volcanic area and its tectonic implication[J]. Chinese Journal of Geophysics,54(8):2050–2059 (in Chinese).
张风雪,吴庆举,丁志峰. 2018. 青藏高原东部P波速度结构及其对高原隆升的启示[J]. 科学通报,63(19):1949–1961. Zhang F X,Wu Q J,Ding Z F. 2018. A P-wave velocity study beneath the eastern region of Tibetan Plateau and its implication for plateau growth[J]. Chinese Science Bulletin,63(19):1949–1961 (in Chinese). doi: 10.1360/N972018-00337
张忠杰,陈赟,田小波. 2009. 青藏高原东缘地壳上地幔结构及其动力学意义[J]. 地质科学,44(4):1136–1150. doi: 10.3321/j.issn:0563-5020.2009.04.009 Zhang Z J,Chen Y,Tian X B. 2009. Crust-upper mantle structure on the eastern margin of Tibet Plateau and its geodynamic implications[J]. Chinese Journal of Geology,44(4):1136–1150 (in Chinese).
郑晨,丁志峰,宋晓东. 2016. 利用面波频散与接收函数联合反演青藏高原东南缘地壳上地幔速度结构[J]. 地球物理学报,59(9):3223–3236. doi: 10.6038/cjg20160908 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[J]. Chinese Journal of Geophysics,59(9):3223–3236 (in Chinese).
Bai D H,Unsworth M J,Meju M A,Ma X B,Teng J W,Kong X R,Sun Y,Sun J,Wang L F,Jiang C S,Zhao C P,Xiao P F,Liu M. 2010. Crustal deformation of the eastern Tibetan Plateau revealed by magnetotelluric imaging[J]. Nat Geosci,3(5):358–362. doi: 10.1038/ngeo830
Bao X W,Sun X X,Xu M J,Eaton D W,Song X D,Wang L S,Ding Z F,Mi N,Li H,Yu D Y,Huang Z C,Wang P. 2015. Two crustal low-velocity channels beneath SE Tibet revealed by joint inversion of Rayleigh wave dispersion and receiver functions[J]. Earth Planet Sci Lett,415:16–24. doi: 10.1016/j.jpgl.2015.01.020
Copley A. 2008. Kinematics and dynamics of the southeastern margin of the Tibetan Plateau[J]. Geophys J Int,174(3):1081–1100. doi: 10.1111/j.1365-246X.2008.03853.x
Courtillot V,Jaupart C,Manighetti I,Tapponnier P,Besse J. 1999. On causal links between flood basalts and continental breakup[J]. Earth Planet Sci Lett,166(3/4):177–195.
Farra V,Vinnik L. 2000. Upper mantle stratification by P and S receiver functions[J].Geophys J Int,141(3):699–712. doi: 10.1046/j.1365-246x.2000.00118.x
Ford H A,Fischer K M,Abt D L,Rychert C A,Elkins-Tanton L T. 2010. The lithosphere-asthenosphere boundary and cratonic lithospheric layering beneath Australia from Sp wave imaging[J]. Earth Planet Sci Lett,300(3/4):299–310.
Fu Y Y,Gao Y,Li A B,Li L,Chen A G. 2017. Lithospheric structure of the southeastern margin of the Tibetan Plateau from Rayleigh wave tomography[J]. J Geophys Res,122(6):4631–4644. doi: 10.1002/2016JB013096
Houston H. 2001. Influence of depth,focal mechanism,and tectonic setting on the shape and duration of earthquake source time functions[J]. J Geophys Res,106(B6):11137–11150. doi: 10.1029/2000JB900468
Hu J F,Yang H Y,Xu X Q,Wen L M,Li G Q. 2012. Lithospheric structure and crust-mantle decoupling in the southeast edge of the Tibetan Plateau[J]. Gondwana Res,22(3/4):1060–1067.
Huang Z C,Wang P,Xu M J,Wang L S,Ding Z F,Wu Y,Xu M J,Mi N,Yu D Y,Li H. 2015. Mantle structure and dynamics beneath SE Tibet revealed by new seismic images[J]. Earth Planet Sci Lett,411:100–111. doi: 10.1016/j.jpgl.2014.11.040
Kennett B L N. 2015. Lithosphere-asthenosphere P-wave reflectivity across Australia[J]. Earth Planet Sci Lett,431:225–235. doi: 10.1016/j.jpgl.2015.09.039
Kennett B L N,Engdahl E R,Buland R. 1995. Constraints on seismic velocities in the Earth from traveltimes[J]. Geophys J Int,122(1):108–124. doi: 10.1111/j.1365-246X.1995.tb03540.x
Kumar M R,Bostock M G. 2006. Transmission to reflection transformation of teleseismic wavefields[J]. J Geophys Res,111(B8):B08306.
Langston C A. 1979. Structure under Mount Rainier,Washington,inferred from teleseismic body waves[J]. J Geophys Res,84(B9):4749–4762. doi: 10.1029/JB084iB09p04749
Sun W J,Kennett B L N. 2016. Receiver structure from teleseisms: Autocorrelation and cross correlation[J]. Geophys Res Lett,43(12):6234–6242. doi: 10.1002/2016GL069564
Sun W J,Kennett B L N. 2017. Mid-lithosphere discontinuities beneath the western and central North China Craton[J]. Geophys Res Lett,44:1302–1310.
Sun W J,Fu L Y,Saygin E,Zhao L. 2018. Insights into layering in the cratonic lithosphere beneath western Australia[J]. J Geophys Res,123(2):1405–1418. doi: 10.1002/2017JB014904
Sun X X,Bao X W,Xu M J,Eaton D W,Song X D,Wang L S,Ding Z F,Mi N,Yu D Y,Li H. 2014. Crustal structure beneath SE Tibet from joint analysis of receiver functions and Rayleigh wave dispersion[J]. Geophys Res Lett,41(5):1479–1484. doi: 10.1002/2014GL059269
Sun Y,Niu F L,Liu H F,Chen Y L,Liu J X. 2012. Crustal structure and deformation of the SE Tibetan Plateau revealed by receiver function data[J]. Earth Planet Sci Lett,349/350:186–197. doi: 10.1016/j.jpgl.2012.07.007
Vinnik L P. 1977. Detection of waves converted from P to SV in the mantle[J]. Phys Earth Planet Inter,15(1):39–45. doi: 10.1016/0031-9201(77)90008-5
Vinnik L P,Avetisjan R A,Mikhailova N G. 1983. Heterogeneities in the mantle transition zone from observations of P-to-SV converted waves[J]. Phys Earth Planet Inter,33(3):149–163. doi: 10.1016/0031-9201(83)90117-6
Wang Y. 2001. Heat flow pattern and lateral variations of lithosphere strength in China mainland:Constraints on active deformation[J]. Phys Earth Planet Inter,126(3/4):121–146.
Wu T F,Zhang S X,Li M K,Qin W B,Zhang C Y. 2016. Two crustal flowing channels and volcanic magma migration underneath the SE margin of the Tibetan Plateau as revealed by surface wave tomography[J]. J Asian Earth Sci,132:25–39. doi: 10.1016/j.jseaes.2016.09.017
Xu L L,Rondenay S,van der Hilst R D. 2007. Structure of the crust beneath the southeastern Tibetan Plateau from teleseismic receiver functions[J]. Phys Earth Planet Inter,165(3/4):176–193.
Xu Y G,Chung S L,Jahn B M,Wu G Y. 2001. Petrologic and geochemical constraints on the petrogenesis of Permian-Triassic Emeishan flood basalts in southwestern China[J]. Lithos,58(3/4):145–168.
Yang H Y,Peng H C,Hu J F. 2017. The lithospheric structure beneath southeast Tibet revealed by P and S receiver functions[J]. J Asian Earth Sci,138:62–71. doi: 10.1016/j.jseaes.2017.02.001
Zhang F X,Wu Q J,Li Y H,Zhang R Q,Sun L,Pan J T,Ding Z F. 2018. Seismic tomography of eastern Tibet: Implications for the Tibetan Plateau growth[J]. Tectonics,37(9):2833–2847. doi: 10.1029/2018TC004977
-
期刊类型引用(2)
1. 周进举，张伟. 青藏高原南缘中部的地壳结构：远震P波尾波自相关的新约束. 地球物理学报. 2023(09): 3662-3673 . 
百度学术
2. 韩光洁，邓文泽. 中国西北、西南地区地震台站地震计方位角偏差估计. 地震地磁观测与研究. 2023(04): 1-8 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 1322
- HTML全文浏览量: 737
- PDF下载量: 142
- 被引次数: 2
