2021, Vol. 64
2. 成都理工大学地球物理学院, 成都 610059
2. College of Geophysics, Chengdu University of Technology, Chengdu 610059, China
在新生代印度与欧亚板块碰撞汇聚过程中，青藏高原发生南北向缩短和地壳增厚.大约在渐新世至中新世，青藏高原由地面隆升转为侧向挤出(Yin and Harrison, 2000).青藏高原东南缘是高原物质侧向挤出和构造逃逸的重要出口(Yin, 2010; Tapponnier et al., 2001; Royden et al., 2008).有关青藏高原东南缘的地壳变形机制主要包括：地壳沿大型走滑断裂向东逃逸(Tapponnier et al., 1982, 2001); 地壳连续变形、增厚(England and Houseman, 1986, 1989)和黏滞性地壳物质流动(Royden et al., 1997, 2008; Royden, 1996; Clark and Royden, 2000; Clark et al., 2003, 2005; Schoenbohm et al., 2006; Klemperer, 2006)等.这些模型都在不同程度上揭示了青藏高原东南缘地壳上地幔的物质运动和构造变形过程.
图 1是青藏高原东南缘地形及构造分区略图.研究区域内有一系列共轭交叉的断裂带(张培震，2008)，如北西向的鲜水河断裂带(XSHF)，北东向的小金河断裂带(XJHF)和龙门山断裂带(LMSF)，南北向的小江断裂带(XJF)等.这些断裂带将该区划分为不同地块，由北向南主要包括松潘甘孜地块(SGB)，川西地块(WSB)，滇中地块(CYB)，滇西南地块(SWYB)等.川滇块体向南东方向的不均衡挤出作用使木里盐源地区由原来的挤压状态转变为拉张状态，形成构造逸出盆地(李勇等，2001).
|
图 1 青藏高原东南缘及相邻区域的地形及构造简图 WQLF—西秦岭断裂带; LMSF—龙门山断裂带; ；HYSF—华蓥山断裂带; XSHF—鲜水河断裂带; JSJF—金沙江断裂带；LCJF—澜沧江断裂带; NJF—怒江断裂带; LTF—理塘断裂带；XJHF—小金河断裂带; ZMHF—则木河断裂带; XJF—小江断裂带; RRF—红河断裂带; LMSB—龙门山造山带；SGB—松潘甘孜地块; WSB—川西地块; MYB—木里盐源盆地；CYB—滇中地块；SCB—四川盆地；YZB—扬子地块；SCHB—华南地块；SWYB—滇西南地块; EHS—喜马拉雅东构造结. Fig. 1 Brief topographic and the tectonic map of the southeastern margin of the Tibetan Plateau WQLF—West Qinling fault; LMSF—Longmenshan fault; ; HYSF—Huayingshan fault; XSHF—Xianshuihe fault; JSJF—Jingshajiang fault; LCJF—Lancangjiang fault; NJF—Nujiang fault; LTF—Litang fault; XJHF—Xiaojinhe fault; ZMHF—Zemuhe fault; XJF—Xiaojiang fault; RRF—Red River fault; LMSB—Longmenshan belt; SGB—Songpan-Ganzi block; WSB—Western Sichuan block; MYB—Muli-Yanyuan basin; CYB—Central Yunnan block; SCB—Sichuan basin; YZB—Yangtze block; SCHB—South China block; SWYB—Southwest Yunnan block; EHS—East Himalayan syntaxis. |
地球内部结构及变形所引起的地震各向异性是地球介质的一个重要性质，它与地壳岩石中的片理、节理及裂缝等微构造排列方向，以及矿物结晶优势方向有密切的关系.青藏高原东南缘的地震各向异性分布记录了地壳上地幔现今结构特征及其历史演化过程.观测该区地震各向异性对研究该区构造演化及动力学问题有着重要的意义.
近20多年来，已发表了一系列青藏高原东南缘地震各向异性研究成果.高原等(2020)综合前人研究的结果，发现青藏高原东南缘下地壳和上地幔顶部存在南北分区特征，各向异性方向在26°N以北近似为NS方向，在南部近似为EW方向，推测在26°N北侧上地幔有较厚的高速体，其南侧呈现近似直立构造.王琼等(2015)利用背景噪声等资料，得到青藏高原东南缘上地壳各向异性特征，壳幔各向异性具有不同的特征与形成机制.地壳各向异性与地表现代运动相符，大体上与主压应力的方向和断裂分布有关(太龄雪等，2015；范莉萍等，2015).Yao等(2008, 2009, 2010)利用环境噪声和远震面波的瑞利波相速度频散曲线反演青藏高原东缘的岩石圈结构，揭示了地壳深部区域的剪切波低速区域以及各向异性随深度的变化，上地壳方位各向异性方向围绕喜马拉雅东构造结作旋转运动，其快速方向通常与主走滑断层平行.Huang等(2015, 2018)、Wei等(2013)、常利军等(2015)利用远震P波和S波的各向异性研究，以及利用在Moho面P波转换为S波的Pms波形时差反演地壳各向异性(Zheng et al., 2018; Kong et al., 2016; Cai et al., 2016; Chen et al., 2012; Chen et al., 2013)，都发现了不同深度地壳各向异性的变化，指出高度变形的青藏高原与稳定的扬子板块地震各向异性差异，中下地壳地震快波方向与地壳低速层方向一致.这些结果对青藏高原东南缘地壳上地幔内部结构及变形动力机制提供了重要依据.
随着观测资料的增多和研究方法的改进，地震各向异性的纵向分辨率有望进一步提高.例如原有的Pms震相只能反应整体地壳内部的平均各向异性，对地壳内不同深度的各向异性不敏感.与前人的研究相比, 本文收集密集的宽频带地震固定台站和流动台阵的观测资料，具有更加完备的数据覆盖，利用环境噪声数据提取Rayleigh波频散曲线，采用多角度频散曲线反演方法，获得了地壳上地幔高分辨率的S波速度与各向异性图像，研究地壳内部物质结构及变形，为青藏高原东南缘地壳上地幔构造演化和动力过程提供了新的证据.
1 数据和方法本文使用的数据来源于以下三个部分：(1)四川与云南地震台网120个固定台站记录的2015年1月至2017年12月的连续波形数据；(2)中国地震台阵I期(中国地震科学台阵, 2006)240个流动台站记录的2011年4月至2013年6月的连续波形数据；(3)川西台阵(刘启元等，2008)200个流动台站记录的2007年6月至2008年8月的连续波形数据.(图 2).
|
图 2 青藏高原东南缘台站分布图 黑色三角代表中国地震台阵I期, 绿色三角代表四川与云南国家数字地震台网，红色三角代表川西台阵. Fig. 2 Distribution map of stations on the eastern margin of the Tibetan Plateau The black triangle represents phase I of the China seismograph array, and the green triangle represents the national digital seismic network of Sichuan and Yunnan province. The red triangle represents the western Sichuan seismic array. |
本文采用以下步骤进行数据处理(Liang and Langston, 2008; Bensen et al., 2007)：
(1) 对原始数据按照采样率为10 Hz进行重采样；(2)去倾斜，去均值，消除仪器响应等；(3)在0.02至0.2 Hz的频带范围内进行滤波；(4)波形互相关提取Rayleigh波波形.
将台站对的垂直分量(Z-Z)按24 h的时段进行互相关计算，再对全年数据进行叠加，获得各个台站对之间的瑞利波经验格林函数(图 3a).对30000多条互相关波形按信噪比(SNR)大于5的标准筛选出21473条频散曲线(图 3b)，使用频率-时间分析(FTAN)算法提取频散曲线的群速度(Levshin et al., 1972).将研究区域划分为0.25°×0.25°的网格.利用棋盘格测试试验来评估研究区域的横向分辨率，设置检测板的网格大小为0.5°×0.5°，速度扰动为各周期平均速度的±0.3 km·s-1，得到检测板如图 4所示.
|
图 3 左图表示部分台站对之间的互相关波形，横、纵坐标分别表示时间和台站之间的距离; 右图表示不同周期瑞利波频散曲线数 Fig. 3 The figure on the left shows the cross-correlation waveform with individual station-pair. The horizontal and vertical axis represent the time and distance between stations, respectively. The figure on the right shows the number of Rayleigh wave dispersion curves in different periods |
|
图 4 周期为10 s，20 s，30 s，40 s的群速度层析成像检测板结果 Fig. 4 Checkboard test results of group velocity tomography with a period of 10 s, 20 s, 30 s, 40 s |
在计算过程中加入平滑参数可以减少反演的误差，但可能降低分辨率(Liang and Langston, 2009).利用不同的平滑参数，随机选取数据迭代计算500次，统计500次反演结果，可以得到速度与各向异性的平均值和标准差.通过选取参数，使得各参数误差较小(图 5)也能确保有较好的分辨率.
|
图 5 左、中、右图依次为周期为20 s，40 s各向异性方向(da)、各向异性强度(dm)、群速度(dv)的误差分布 Fig. 5 The left, middle and right images show the error of anisotropy direction (da), anisotropy magnitude (dm) and group velocity (dv) with a period of 20 s and 40 s |
在各向异性介质中，使用多角度频散曲线反演方法(Liang et al., 2020)，在每个块体中对某一周期不同方位的群速度通过拟合来确定相应的各向异性常数ak与bk，再使用线性反演方法(Hermann, 2013)得到VS速度随深度变化的一维速度模型.最后，将所有块体模型进行组合得到三维速度模型.对各向异性介质而言，网格中的走时差与平均慢度和方位各向异性有关.
|
(1) |
式中Δtij是第i个台站到第j个台站的走时差；dijk是第i个台站到第j个台站射线在第k个网格中的长度；sk是第k个网格的慢度.ak和bk是第k个网格中的各向异性参数，ϕijk是第i个台站到第j个台站射线在第k个网格中的方位角，N为网格总数.
将公式(1)化简得到第k个网格中群速度与慢度和各向异性的参数之间的关系：
|
(2) |
第k个网格中各向同性速度vk0，各向异性强度γk和快波方向ϕk分别为
|
(3) |
对于每个网格，得到一个群速度参数和两个各向异性参数(ak和bk)，利用反演的ak和bk参数计算出各向异性强度与快波方向.
2 结果 2.1 不同周期的群速度和各向异性图根据高密度台网的资料，反演得到周期范围为5至50 s的群速度图像.图 6显示了10 s、20 s、30 s和40 s的群速度与各向异性分布图.
|
图 6 周期为10 s, 20 s, 30 s，40 s瑞利面波群速度与各向异性分布 Fig. 6 Images of the group velocities and anisotropies of the Rayleigh waves at various periods of 10, 20, 30 and 40 s, respectively |
由图 6可见，10 s周期群速度四川盆地明显低于其他地区，反映了盆地地壳上部为中新生代的沉积层.青藏高原东南缘与滇东南，黔西南群速度较高.青藏高原东南缘的北部快波方向显示为北西向，中南部逐渐转为南西向.20 s周期与10 s周期群速度大体相似.30 s周期，四川盆地与华南地块群速度较高，松潘甘孜地块出现明显低速异常.40 s周期，松潘甘孜地块仍显示低速分布，川西地块的低速异常范围明显减小，各向异性方向转为南西.四川盆地、扬子地块及华南地块都显示为高速分布.
由于研究区域地壳厚度有较大的变化(Wang et al., 2017；朱介寿等, 2017, Huang et al., 2018)，如图 7可所示，在地下同一深度，各地块所反映的地壳上地幔的部位有较大的差异.在讨论地块某一块体各向异性和速度分布时，必须根据该图来确定其地壳及上地幔所在深度范围.
|
图 7 青藏高原东南缘及相邻地区Moho界面深度图(根据Wang et al., 2017; 朱介寿等, 2017等修改) Fig. 7 Depth map of the Moho interface on the eastern margin of the Tibet Plateau and adjacent areas (modified by Wang et al., 2017; Zhu et al., 2017) |
基于多角度频散曲线反演方法(Liang et al., 2020)，获得了不同深度的速度和各向异性图像(图 8).
|
图 8 (a)、(b)、(c)、(d)表示青藏高原东南缘及相邻区域上地壳顶部(2~5 km), 上地壳(10~20 km), 中下地壳(30~40 km)，下地壳底部以及上地幔顶部(50~60 km)的S波速度和各向异性分布图 红色玫瑰图为各地块的各向异性方向统计结果, 紫色箭头显示青藏高原东南缘上地壳及下地壳底部及上地幔顶部的各向异性方向，蓝色箭头表示中下地壳低速层分布区域各向异性方位，(d)中TC表示腾冲火山活动带. Fig. 8 (a)、(b)、(c)、(d) show the S-wave velocity and anisotropy maps of the top upper crust (2~5 km); upper crust (10~20 km), mid-lower crust (30~40 km), bottom of lower crust to top of upper mantle (50~60 km) in the southeastern margin of the Tibetan Plateau and adjacent areas The red rose diagrams show the statistical results of the anisotropy direction of each block, the purple arrows. Show the anisotropy direction of the southeastern of the Tibetan Plateau, the blue arrows show the anisotropy direction of the low velocity zone in mid-lower crust. TC represents the Tengchong volcanic activity zone in (d). |
在2~5 km深度(图 8a)，青藏高原东南缘，扬子地块等区域，都位于上地壳顶部.松潘甘孜地块各向异性方位角为NW向，与鲜水河断裂带走向一致.川西地块各向异性由NW转为SW向.木里盐源盆地和滇中地块都为NE-SW向，与小金河断裂带走向一致，滇西南地块亦为SW向.滇中地块由NNW逐渐转为NS向，与则木河断裂带和小江断裂带等平行.青藏高原东南缘各地块显示的各向异性方向，整体上是围绕喜马拉雅东构造结作顺时针旋转运动的, 如图10a中的紫色箭头所示.四川盆地及扬子地块上地壳2~5 km的各向异性特征如下.龙门山造山带为NE向，与龙门山断裂带走向基本一致，四川盆地为NE向，与盆地内部的主要断裂带，如华蓥山断裂带的走向一致.扬子地块及华南地块各向异性特征优势方向为SE向.
在10~20 km深度(图 8b)，青藏高原东南缘为上地壳底部, 四川盆地，扬子地块，华南地块及滇西南地块已为中地壳.各地块的各向异性方向(图 8b)整体上仍围绕东构造结(EHS)作顺时针旋转运动(如图中紫色箭头所示).四川盆地，扬子地块以及华南地块的各向异性方位角大体从NS向再转为SE向.
在30~40 km深度(图 8c)，松潘甘孜地块及川西地块在这一深度位于中下地壳，滇中及滇西南位于下地壳与上地幔顶部.四川盆地，扬子地块为地壳底层，华南地块为地壳底部或上地幔顶部.青藏高原东南缘中下地壳方位角相对上地壳发生了很大变化，松潘甘孜地块方向仍为NNW，川西地块仍为NE向，但是木里盐源盆地由上地壳的NE向逆转为NW向，到达滇中地块也由NE向逆转为NW向.
低速层由松潘甘孜地块向东，向南折向川西地块转为NS及NE向，在木里盐源盆地转为NW向，向东转进滇中地块.也有部分低速层物质进入扬子地块，使扬子地块西缘局部地区也成为低速带，如图 8c中蓝色箭头所示.
毗邻的扬子地块及华南地块在这一深度已接近地壳下部.龙门山造山带方位角仍为NE向，与龙门山断裂带走向一致.四川盆地为NS或NNW走向，扬子地块为NNE向，华南地块仍为SE向.中下地壳各向异性与上地壳分布大体相同，仅局部有一些变化.
在50~60 km深度(图 8d), 青藏高原东南缘及相邻区域为下地壳底部及上地幔顶部.这一深度地震各向异性分布与上地壳相比，除川西地块稍有偏离外，其余如松潘甘孜地块，木里盐源盆地、滇中地块等与上地壳分布基本相同.各地块总体趋势仍然是围绕东构造结(EHS)作顺时针环绕运动, 如图中紫色箭头所示.其中西南缘的低速分布带可能由腾冲(TC)火山活动引起.
四川盆地及扬子地块这一深度已位于上地幔顶部，从龙门山造山带至四川盆地扬子板块和华南板块，各地块的方位角排列总体上仍然与上地壳和中下地壳相似.四川盆地及华南地块均位于上地幔顶部，各向异性强度(MOA)较地壳有所增大.
3 讨论为了解中下地壳低速层的分布，图 9显示了不同方位的速度剖面.地壳低速层的上界面位于地下18~20 km深度，下界面位于30~42 km深度，地壳低速层的厚度变化范围为10~25 km.松潘甘孜地块与川西地块，木里盐源盆地、滇中地块，厚度可达20~25 km左右.低速层向东及东南逐渐减薄，到龙门山断裂带和小江断裂带，其厚度减小到10~15 km左右.
|
图 9 A、B、C、D、E为五条地壳速度剖面 红色线条表示中下地壳低速层上下界面，红色箭头表示低速层物质移动方向，黑色线条表示莫霍界面. Fig. 9 A、B、C、D、E represent five crustal velocity profiles Red line represents the upper and lower interface of the low velocity layer in the mid-lower crust (crustal flow), the red arrows show the moving direction of the low velocity layer matter, and black line indicate the Moho boundary. |
地震各向异性是地球介质的一个重要特征.由于大陆地壳具有漫长的历史，地壳上地幔不同深度的地震各向异性，记录了它的现今构造特征及其历史演变过程, 为研究构造演化及动力学过程提供重要依据.
在脆性变形的上地壳中，各向异性可以由岩体中微破裂及微裂纹的排列、岩石层理或熔体组织显示优势方向引起, 这称为形状优势方向(SPO)(Crampin and Chastin, 2003, 2008).在许多情况下，SPO受应力场的控制，快波方向平行于岩体的最大剪切或拉张应变轴，与构造活动区的断裂带的走向或其边界一致(Dreiling et al., 2018, Huang et al., 2015).
在塑性变形的中下地壳和上地幔中，各向异性可由岩石中的矿物结构特性产生.如中下地壳中的云母、角闪石和上地幔中的橄榄石等矿物的晶格优势方向(LPO)及晶体优势方向(CPO)(Mainprice, 2007；Brownlee et al., 2017)与地震快波方向一致.
地球动力学问题大都与挤压或拉张应力有关，由于应力方向重新定向，将影响到地壳中矿物的LPO和CPO, 使各向异性发生变化.对于软流层中的地幔流, 各向异性的快速方向一般与地幔流方向相同(Karato et al., 2008).
青藏高原东南缘上地壳各地块环绕东构造结(ESH)作顺时针旋转运动, 与青藏高原东南缘相对于扬子地块地面GPS水平运动方向是一致的(Gan et al., 2007；王阎昭等，2008; Shen et al., 2005；Pan and Shen, 2017；Wang and Shen, 2020).断层最大走滑量沿鲜水河断裂带转为则木河断裂带(ZMHF)至小江断裂带.这些走滑断裂带也与青藏高原东南缘的最大剪切应变带相符合(Pan and Shen, 2017；Gan et al., 2007；Huang et al., 2018).
青藏高原东南缘地壳低速层的物质来源及性质，以及动力过程已有大量研究成果.其中包括地壳流导致地壳加厚和地表隆升机制的研究(Royden, 1996; Royden et al., 1997；Schoenbohm et al., 2006；Clark et al., 2003, 2005; Clark and Royden, 2000；Burchfiel et al., 2008；张培震，2008；Yi et al., 2016)；基于天然地震对中下地壳低速层及地壳流的研究(Liu et al., 2014, 2018；Yang et al., 2012；王椿镛等，2008；朱介寿等，2017；李昱等，2010；Yao et al., 2008, 2009；范莉苹等，2015；郑晨等，2016；Bao et al., 2015).通过对地壳温度及深部岩石包体分析，认为青藏高原地壳内20~40 km深观测到的低速层是在700 ℃~850 ℃温度下，脱水含云母的硅酸岩部份熔融的黏滞性流体，其晶体优势方向(CPO)与中地壳快波方向一致(Hacker et al., 2014, Xie et al., 2013, 2017).
青藏高原东南缘地壳中部的黏滞性流体在向东及东南方向运移时，除了带动地表各地块以不同速度移动外，它的动力学效应还表现在地表隆升及地壳增厚.青藏高原东南缘海拔1500 m至4000 m以上的高原地区，地壳内普遍存在黏滞性低速层，而在1500 m以下区域，几乎没有低速层.
青藏高原东南缘的上地壳与下地壳各向异性特征相一致，表明它们在较早的历史演化时期(如古生代至中生代)可能是耦合在一起的.在新生代时期，青藏高原隆升和向东挤出，低速的黏滞性物质流入中地壳，使原有上下地壳分离解耦，但它们的相对位置及各向异性特征仍然不变.
稳定的扬子地块及华南地块为前寒武纪克拉通地区.扬子地块的核心区(四川盆地)基本没有变形. 在漫长的地质历史中，从地壳到上地幔顶部，各向异性方位角基本保持不变，因此，这里显示的可能是古构造活动遗留下来的各向异性(Huang et al., 2015).
远震P波对青藏高原东南缘及相邻区域各向异性成像(Wei et al., 2013)发现，在中下地壳中存在一个显著的厚约20 km的低速地壳流，该地壳流在向东流动时受刚性的四川盆地阻挡，地震各向异性也发生变化.在上地壳发生的强烈地震都分布在塑性流体层的边界，表明地震活动受深部动力过程控制.这一研究与本文得到的中下低速地壳流分布范围与动力过程是相同的.用P波研究青藏高原东南缘的成果(Huang et al., 2018)指出，青藏高原东南缘围绕喜马拉雅东构造结(EHS)顺时针旋转同时，扬子地台西缘部分地区受到低速物质侵入影响，成为构造活动区.地壳上地幔各向异性显示，地震快波方向与地面主要断裂带走向一致.
利用Moho界面转换波PmS的正弦波时差反演青藏高原东南缘各向异性(Zheng et al., 2018) 发现松潘甘孜、川西、滇中等地块各向异性方向与地面主干断裂带平行.中下地壳各向异性快波方向与低速地壳流流向一致.利用Moho界面及地壳内的速度界面的P波转换为S波的波形时差反演的各向异性成果指出(Kong et al., 2016)，快波方向与地表主干断裂平行，而且主要集中在地块的边界部份，并显示中下地壳存在塑性地壳流.基于程函方程的面波层析成像，对青藏高原东南缘各向异性及速度成像，发现20 km以上地震各向异性与地面断裂带走向及最大主压应力方向一致.在川滇菱形块体北部，存在明显的低速地壳流分布，地壳流流向与块体向南挤出方向一致(王怀富等，2020).
由上可见，本文所得到的青藏高原东南缘地震各向异性的研究结果，与近年来发表的相关结论基本是一致的.由于我们的反演具有更高的纵向分布率，可以获得地壳及上地幔不同部位的速度及各向异性的细节特征，能对其物质结构及动力过程进行更深入研究.
4 结论对青藏高原东南缘及相邻地区中国地震台阵I期，四川与云南国家数字地震台网，川西台阵观测资料进行处理, 用环境噪声成像提取瑞利面波, 采用面波多方位角反演方法，获得青藏高原东南缘及相邻地区高分辨率的S波各向异性及速度分布图像.
在上地壳(从地表到20 km深度), 青藏高原东南缘各地块的快波方向与其邻近的走滑断裂带走向一致.各向异性形状优势方向(SPO)与相对于华南地区的地面GPS水平速度场方向一致，也与青藏高原东南缘的最大剪切应变带相符合，都是围绕喜马拉雅东构造结作顺时针旋转运动.
在中下地壳(30~40 km深度)，地震各向异性相对上地壳发生了很大变化.地震各向异性方位角从北向南在滇中地块转为北东向.各向异性快波方向由中下地壳中厚约10~25 km的低速流体物质流向所控制.
下地壳底部至上地幔顶部(50~60 km深度)，青藏高原东南缘的各向异性分布仍与上地壳大体相同，总体上仍围绕喜马拉雅东构造结(EHS)作顺时针旋转运动.青藏高原东南缘的上地壳与下地壳各向异性特征相一致，表明它们在较早的历史演化时期(如古生代至中生代)可能是耦合在一起的.
根据热年代学测定结果，发现青藏高原东南缘在新生代至少有两期隆升事件，分别发生在渐新世30~25 Ma及中新世15~10 Ma, 黏滞性流体进入中下地壳最可能是后一时期(Wang et al., 2012).对青藏高原新生代隆升及挤出的模拟，也得出中新世15 Ma以来为青藏高原东南缘隆升及低速地壳流活动时期的结论(Yang et al., 2009).
在新生代时期，随着青藏高原隆升向东挤出，特别在中新世以来低速的黏滞性物质流入青藏高原东南缘的中下地壳，改变了各向异性方向，但原有上地壳与下地壳底部的相对位置及各向异性特征仍保持不变.青藏高原东南缘地壳整体沿喜马拉雅山东构造结(EHS)顺时针旋转，上地壳与下地壳底部作同步运动，仍保持了原有各向异性特征.但流入中下地壳的黏滞性流体使各向异性方向和速度分布发生颠转变化.新生代以来青藏高原东南缘地壳结构的变化及外部物质的渗入，是地壳变形和构造活动的根本原因.
四川盆地、扬子地块及华南地块从地壳到上地幔顶部，各向异性特征基本保持不变.它们的构造特征及演化过程与青藏高原东南缘是不同的.扬子地块及华南地块为稳定的前寒武纪克拉通地区.在漫长的地质历史中，显示出化石各向异性特征.
致谢 感谢“中国地震科学探测台阵数据中心”和“川西台阵”，及四川省地震局、云南地震局为本文提供的地震数据资料；感谢两位匿名审稿人对本文认真严谨的审阅，并提出了宝贵的修改意见和建议.本文的图件均由GMT，Surfer，Matlab绘制.
Bai D H, Unsworth M J, Meju M A, et al. 2010. Crustal deformation of the eastern Tibetan plateau revealed by magnetotelluric imaging. Nature Geoscience, 3(5): 358-362. DOI:10.1038/ngeo830 |
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 functions. Earth and Planetary Science Letters, 415: 16-24. DOI:10.1016/j.epsl.2015.01.020 |
Bensen G D, 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 of the Royal Astronomical Society, 169(3): 1239-1260. DOI:10.1111/j.1365-246X.2007.03374.x |
Brownlee S J, Schulte-Pelkum V, Raju A, et al. 2017. Characteristics of deep crustal seismic anisotropy from a compilation of rock elasticity tensors and their expression in receiver functions. Tectonics, 36(9): 1835-1857. DOI:10.1002/2017TC004625 |
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 |
Cai Y, Wu J P, Fang L H, et al. 2016. Crustal anisotropy and deformation of the southeastern margin of the Tibetan Plateau revealed by Pms splitting. Journal of Asian Earth Sciences, 121: 120-126. DOI:10.1016/j.jseaes.2016.02.005 |
Chang L J, DING Z F, Wang C Y. 2015. Upper mantle anisotropy beneath the southern segment of North-South tectonic belt, China. Chinese Journal of Geophysics (in Chinese), 58(11): 4052-4067. DOI:10.6038/cjg20151114 |
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 |
ChinArray. 2006. China Seismic Array waveform data (in Chinese). China Earthquake Adminstration. DOI:10.12001/ChinArray.Data
|
Clark M K, Royden L H. 2000. Topographic ooze: building the eastern margin of Tibet by lower crustal flow. Geology, 28(8): 703-706. DOI:10.1130/0091-7613(2000)28<703:TOBTEM>2.0.CO;2 |
Clark M K, Royden L, Burchfiel B C, et al. 2003. Late Cenozoic uplift of southeastern Tibet: implications for lower crustal flow. Geophysical Research Abstracts, 5: 12969. |
Clark M K, House M A, Royden L H, et al. 2005. Late Cenozoic uplift of southeastern Tibet. Geology, 33(6): 525-528. DOI:10.1130/G21265.1 |
Crampin S, Chastin S. 2003. A review of shear wave splitting in the crack-critical crust. Geophysical Journal International, 155(1): 221-240. DOI:10.1046/j.1365-246X.2003.02037.x |
Crampin S, Peacock S. 2008. A review of the current understanding of seismic shear-wave splitting in the Earth's crust and common fallacies in interpretation. Wave Motion, 45(6): 675-722. DOI:10.1016/j.wavemoti.2008.01.003 |
Dreiling J, Tilmann F, Yuan X H, et al. 2018. Crustal radial anisotropy and linkage to geodynamic processes: A study based on seismic ambient noise in southern Madagascar. Journal of Geophysical Research: Solid Geophysical Research: Solid Earth, 123(6): 5130-5146. DOI:10.1029/2017JB015273 |
England P, Houseman G. 1986. Finite strain calculations of continental deformation: 2. Comparison with the India-Asia collision zone. Journal of Geophysical Research: Solid Earth, 91(B3): 3664-3676. DOI:10.1029/JB091iB03p03664 |
England P, Houseman G. 1989. Extension during continental convergence, with application to the Tibetan Plateau. Journal of Geophysical Research: Solid Earth, 94(B12): 17561-17579. DOI:10.1029/JB094iB12p17561 |
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 implications. Chinese Journal of Geophysics (in Chinese), 58(5): 1555-1567. DOI:10.6038/cjg20150509 |
Gan W J, Zhang P Z, Shen Z K, et al. 2007. Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements. Journal of Geophysical Research: Solid Earth, 112(B8): B08416. DOI:10.1029/2005JB004120 |
Gao Y, Shi Y T, Wang Q. 2020. Seismic anisotropy in the southeastern margin of the Tibetan Plateau and its deep tectonic significances. Chinese Journal of Geophysics (in Chinese), 63(3): 802-816. DOI:10.6038/cjg2020O0003 |
Hacker B R, Ritzwoller M H, Xie J. 2014. Partially melted, mica-bearing crust in Central Tibet. Tectonics, 33(7): 1408-1424. DOI:10.1002/2014TC003545 |
Herrmann R B. 2013. Computer programs in seismology: An evolving tool for instruction and research. Seismological Research Letters, 84(6): 1081-1088. DOI:10.1785/0220110096 |
Huang Z C, Wang L S, Xu M J, et al. 2015. Teleseismic shear-wave splitting in SE Tibet: Insight into complex crust and upper-mantle deformation. Earth and Planetary Science Letters, 432: 354-362. DOI:10.1016/j.epsl.2015.10.027 |
Huang Z C, Wang L S, Xu M J, et al. 2018. P wave anisotropic tomography of the SE Tibetan Plateau: Evidence for the crustal and upper-mantle deformations. Journal of Geophysical Research: Solid Earth, 123(10): 8957-8978. DOI:10.1029/2018JB016048 |
Karato S I, Jung H, Katayama I, et al. 2008. Geodynamic significance of seismic anisotropy of the upper mantle: new insights from laboratory studies. Annual Review of Earth and Planetary Sciences, 36: 59-95. DOI:10.1146/annurev.earth.36.031207.124120 |
Klemperer S L. 2006. Crustal flow in Tibet: geophysical evidence for the physical state of Tibetan lithosphere, and inferred patterns of active flow. Geological Society, London, Special Publications, 268(1): 39-70.
|
Kong F S, Wu J, Liu K H, et al. 2016. Crustal anisotropy and ductile flow beneath the eastern Tibetan Plateau and adjacent areas. Earth and Planetary Science Letters, 442: 72-79. DOI:10.1016/j.epsl.2016.03.003 |
Levshin A L, Pisarenko V F, Pogrebinsky G A. 1972. On a frequency-time analysis of oscillations. Annales Geophysicae, 28(2): 211-218. |
Li Y, Hou Z J, Si G Y, et al. 2001. Sedimentary characteristics and tectonic controls in Neogene Yanyuan tectonic escape basin in Southeastern Qinghai-Tibet Plateau. Journal of Mineralogy and Petrology (in Chinese), 21(3): 34-43. DOI:10.3969/j.issn.1001-6872.2001.03.006 |
Li Y, Yao H J, Liu Q Y, et al. 2010. Phase velocity array tomography of Rayleigh waves in western Sichuan from ambient seismic noise. Chinese Journal of Geophysics (in Chinese), 53(4): 842-852. DOI:10.3969/j.issn.0001-5733.2010.04.009 |
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 |
Liang C T, Liu Z Q, Hua Q, et al. 2020. The 3D seismic azimuthal anisotropies and velocities in the eastern Tibetan Plateau extracted by an azimuth-dependent dispersion curve inversion method. Tectonics, 39(4): e2019TC005747. DOI:10.1029/2019TC005747 |
Liu Q, 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. DOI:10.1038/ngeo2130 |
Liu Q Y, Chen J H, Li S C, et al. 2008. The MS8.0 Wenchuan earthquake: preliminary results from the Western Sichuan mobile seismic array observations. Seismology and Geology (in Chinese), 30(3): 584-596. DOI:10.3969/j.issn.0253-4967.2008.03.002 |
Liu Q Y, Li Y, Chen J H, et al. 2009. Wenchuan MS8.0 earthquake: preliminary study of the S-wave velocity structure of the crust and upper mantle. Chinese Journal of Geophysics (in Chinese), 52(2): 309-319. |
MAINPRICE D. 2007. Seismic Anisotropy of the Deep Earth from a Mineral and Rock Physics Perspective. Treatise on Geophysics: 437-491. DOI:10.1016/b978-044452748-6/00045-6 |
Pan Y J, Shen W B. 2017. Contemporary crustal movement of southeastern Tibet: Constraints from dense GPS measurements. Scientific Reports, 7: 45348. DOI:10.1038/srep45348 |
Royden L H. 1996. Coupling and decoupling of crust and mantle in convergent orogens: implications for strain partitioning in the crust. Journal of Geophysical Research: Solid Earth, 101(B8): 17679-17705. DOI:10.1029/96JB00951 |
Royden L H, Burchfiel B C, King W R, et al. 1997. Surface deformation and lower crustal flow in eastern Tibet. Science, 276(5313): 788-790. DOI:10.1126/Science.276.5313.788 |
Royden L H, Burchfiel B C, Van Der Hilst R D. 2008. The geological evolution of the Tibetan Plateau. Science, 321(5892): 1054-1058. DOI:10.1126/science.1155371 |
Schoenbohm I M, Burchfiel B C, Chen L Z. 2006. Propagation of surface uplift, lower crustal flow, and Cenozoic tectonics of the southeast margin of the Tibetan plateau. Geology, 34(10): 813-816. DOI:10.1130/G22679.1 |
Shen Z K, Lu 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 |
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 |
Tai L X, Gao Y, Liu G, et al. 2015. Crustal seismic anisotropy in the southeastern margin of Tibetan plateau by China Array data: shear-wave splitting from temporary observation of the first phase. Chinese Journal of Geophysics (in Chinese), 58(11): 4079-4091. DOI:10.6038/cjg20151116 |
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 |
Tapponnier P, Xu Z Q, Roger F, et al. 2001. Oblique stepwise rise and growth of the Tibet Plateau. Science, 294(5547): 1671-1677. DOI:10.1126/science.105978 |
Wang C Y, Lou H, Lv 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 E, Kirby E, Furlong K P, et al. 2012. Two-phase growth of high topography in eastern Tibet during the Cenozoic. Nature Geoscience, 5(9): 640-645. DOI:10.1038/ngeo1538 |
Wang H F, Wu J P, Zhou S Y, et al. 2020. Rayleigh wave azimuthal anisotropy in the southeastern Tibetan plateau from Eikonal tomography. Chinese Journal of Geophysics (in Chinese), 63(3): 1070-1084. DOI:10.6038/cjg2020N0104 |
Wang M, Shen Z K. 2020. Present-day crustal deformation of continental China derived from GPS and its tectonic implications. Journal of Geophysical Research: Solid Earth, 125(2): e2019JB018774. DOI:10.1029/2019JB018774 |
Wang Q, Gao Y, Shi Y T. 2015. Rayleigh wave azimuthal anisotropy on the southeastern front of the Tibetan plateau from seismic ambient noise. Chinese Journal of Geophysics (in Chinese), 58(11): 4068-4078. DOI:10.6038/cjg20151115 |
Wang W L, Wu J P, Fang L H, et al. 2017. Crustal thickness and Poisson's ratio in Southwest China based on data from dense seismic arrays. Journal of Geophysical Research: Solid Earth, 122(9): 7219-7235. DOI:10.1002/2017JB013978 |
Wang X B, Yu N, Gao S, et al. 2017. Research progress in research on electrical structure of crust and upper mantle beneath the eastern margin of the Tibetan plateau. Chinese Journal of Geophysics (in Chinese), 60(6): 2350-2370. DOI:10.6038/cjg20170626 |
Wang Y Z, Wang E N, Shen Z K, et al. 2008. GPS-constrained inversion of present-day slip rates along major faults of the Sichuan-Yunnan region, China. Science in China Series D: Earth Sciences, 51(9): 1267-1283. DOI:10.1007/s11430-008-0106-4 |
Wei W, Zhao D P, Xu J D. 2013. P-wave anisotropic tomography in Southeast Tibet: New insight into the lower crustal flow and seismotectonics. Physics of the Earth and Planetary Interiors, 222: 47-57. DOI:10.1016/j.pepi.2013.07.002 |
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 |
Xie J Y, Ritzwoller M H, Shen W S, et al. 2017. Crustal anisotropy across eastern Tibet and surroundings modeled as a depth-dependent tilted hexagonally symmetric medium. Geophysical Journal International, 209(1): 466-491. DOI:10.1093/gji/ggx004 |
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 |
Yang Y Q, Liu M. 2009. Crustal thickening and lateral extrusion during the Indo-Asian collision: A 3D viscous flow model. Tectonophysics, 465(1-4): 128-135. DOI:10.1016/j.tecto.2008.11.002 |
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 structure.. Geophysical Journal International, 173(1): 205-219. DOI:10.1111/j.1365-246X.2007.03696.x |
Yao H J, Robert D, Van Der Hilst E D. 2009. Analysis of ambient noise energy distribution and phase velocity bias in ambient noise tomography, with application to SE Tibet. Geophysical Journal International, 179(2): 1113-1132. DOI:10.1111/j.1365-246X.2009.04329.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: Solid Earth, 115(B12): B12307. DOI:10.1029/2009JB007142 |
Yi S, Freymueller J T, Sun W K. 2016. How fast is the middle-lower crust flowing in eastern Tibet? A constraint from geodetic observations. Journal of Geophysical Research: Solid Earth, 121(9): 6903-6915. DOI:10.1002/2016JB013151 |
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 |
Yin A. 2010. Cenozoic tectonic evolution of Asia: A preliminary synthesis. Tectonophysics, 488(1-4): 293-325. DOI:10.1016/j.tecto.2009.06.002 |
Zhang P Z. 2008. A study on the present tectonic deformation, strain partitioning and deep dynamic process of west Sichuan region on eastern margin of Qinghai-Tibet Plateau. Science in China (Series D), 38(9): 1041-1056. |
Zhao G Z, Chen X B, Wang L F, et al. 2008. Evidence of crustal 'channel flow' in the eastern margin of Tibetan Plateau from MT measurements. Chinese Science Bulletin (in Chinese), 53(12): 1887-1893. |
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 |
Zheng C, Ding Z F, Song X D. 2018. Joint inversion of surface wave dispersion and receiver functions for crustal and uppermost mantle structure beneath the northern north-south seismic zone. Chinese Journal of Geophysics (in Chinese), 61(4): 1211-1224. DOI:10.6038/cjg2018L0443 |
Zheng T, Ding Z F, Ning J Y, et al. 2018. Crustal azimuthal anisotropy beneath the southeastern Tibetan Plateau and its geodynamic implications. Journal of Geophysical Research: Solid Earth, 123(11): 9733-9749. DOI:10.1029/2018JB015995 |
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 dynamics. Chinese Journal of Geophysics (in Chinese), 60(6): 2038-2057. DOI:10.6038/cjg20170602 |
常利军, 丁志峰, 王椿镛. 2015. 南北构造带南段上地幔各向异性特征. 地球物理学报, 58(11): 4052-4067. DOI:10.6038/cjg20151114 |
范莉苹, 吴建平, 房立华, 等. 2015. 青藏高原东南缘瑞利波群速度分布特征及其构造意义探讨. 地球物理学报, 58(5): 1555-1567. DOI:10.6038/cjg20150509 |
高原, 石玉涛, 王琼. 2020. 青藏高原东南缘地震各向异性及其深部构造意义. 地球物理学报, 63(3): 802-816. |
李勇, 侯中健, 司光影, 等. 2001. 青藏高原东南缘晚第三纪盐源构造逸出盆地的沉积特征与构造控制. 矿物岩石, 21(3): 34-43. DOI:10.3969/j.issn.1001-6872.2001.03.006 |
李昱, 姚华建, 刘启元, 等. 2010. 川西地区台阵环境噪声瑞利波相速度层析成像. 地球物理学报, 53(4): 842-852. DOI:10.3969/j.issn.0001-5733.2010.04.009 |
刘启元, 陈九辉, 李顺成, 等. 2008. 汶川MS8.0地震: 川西流动地震台阵观测数据的初步分析. 地震地质, 30(3): 584-596. DOI:10.3969/j.issn.0253-4967.2008.03.002 |
太龄雪, 高原, 刘庚, 等. 2015. 利用中国地震科学台阵研究青藏高原东南缘地壳各向异性: 第一期观测资料的剪切波分裂特征. 地球物理学报, 58(11): 4079-4091. DOI:10.6038/cjg20151116 |
王椿镛, 楼海, 吕智勇, 等. 2008. 青藏高原东部地壳上地幔S波速度结构-下地壳流的深部环境. 中国科学D辑: 地球科学, 38(1): 22-32. DOI:10.3321/j.issn:1006-9267.2008.01.003 |
王怀富, 吴建平, 周仕勇, 等. 2020. 青藏高原东南缘基于程函方程的面波方位各向异性研究. 地球物理学报, 63(3): 1070-1084. DOI:10.6038/cjg2020N0104 |
王琼, 高原, 石玉涛. 2015. 青藏高原东南缘基于背景噪声的Rayleigh面波方位各向异性研究. 地球物理学报, 58(11): 4068-4078. DOI:10.6038/cjg20151115 |
王阎昭, 王恩宁, 沈正康, 等. 2008. 基于GPS资料约束反演川滇地区主要断裂现今活动速率. 中国科学D辑: 地球科学, 38(5): 582-597. DOI:10.3321/j.issn:1006-9267.2008.05.006 |
张培震. 2008. 青藏高原东缘川西地区的现今构造变形、应变分配与深部动力过程. 中国科学D辑: 地球科学, 38(9): 1041-1056. DOI:10.3321/j.issn:1006-9267.2008.09.001 |
郑晨, 丁志峰, 宋晓东. 2016. 利用面波频散与接收函数联合反演青藏高原东南缘地壳上地幔速度结构. 地球物理学报, 59(9): 3223-3236. DOI:10.6038/cjg20160908 |
中国地震科学台阵. 2006. 中国地震科学探测台阵波形数据. 中国地震局. DOI:10.12001/ChinArray.Data
|
朱介寿, 王绪本, 杨宜海, 等. 2017. 青藏高原东缘的地壳流及动力过程. 地球物理学报, 60(6): 2038-2057. DOI:10.6038/cjg20170602 |