南海上地幔各向异性结构及动力学含义

薛梅, 李琳, 杨挺, 刘晨光, 华清峰, 夏少红, 黄海波, 黎伯孟, 霍达, 潘谟晗

薛梅,李琳,杨挺,刘晨光,华清峰,夏少红,黄海波,黎伯孟,霍达,潘谟晗. 2023. 南海上地幔各向异性结构及动力学含义. 地震学报,45(3):494−520. DOI: 10.11939/jass.20230054
引用本文: 薛梅,李琳,杨挺,刘晨光,华清峰,夏少红,黄海波,黎伯孟,霍达,潘谟晗. 2023. 南海上地幔各向异性结构及动力学含义. 地震学报,45(3):494−520. DOI: 10.11939/jass.20230054
Xue M,Li L,Yang T,Liu C G,Hua Q F,Xia S H,Huang H B,Le B M,Huo D,Pan M H. 2023. Anisotropic structure and dynamic implications of the upper mantle in the South China Sea. Acta Seismologica Sinica45(3):494−520. DOI: 10.11939/jass.20230054
Citation: Xue M,Li L,Yang T,Liu C G,Hua Q F,Xia S H,Huang H B,Le B M,Huo D,Pan M H. 2023. Anisotropic structure and dynamic implications of the upper mantle in the South China Sea. Acta Seismologica Sinica45(3):494−520. DOI: 10.11939/jass.20230054

南海上地幔各向异性结构及动力学含义

基金项目: 国家自然科学基金面上项目(42076064)和上海佘山地球物理国家野外科学观测研究站开放基金(2020K03)共同资助
详细信息
    通讯作者:

    薛梅,博士,教授,主要从事海洋地球物理与地球动力学方面的研究,e-mail:meixue@tongji.edu.cn

  • 中图分类号: P315.2

Anisotropic structure and dynamic implications of the upper mantle in the South China Sea

  • 摘要: 南海处于欧亚板块、太平洋板块和印度—澳大利亚板块的交会区,是西北太平洋一系列边缘海中最大的边缘海。关于南海的打开以往研究提出了如印度板块与欧亚板块碰撞驱动挤出以及古南海俯冲拖拽等诸多模型。本文力图通过南海海盆及周边各向异性结构来约束南海演化机制。基于同济大学2012和2014年在南海中央海盆进行的两次被动源宽频带海底地震观测试验回收的10台OBS记录仪近1年的地震数据,本文采用三种不同的横波分裂方法,获取了中央海盆针对两次远震的XKS分裂结果以及南海周边20次区域地震提供的S震相分裂结果。SKS分裂结果显示,南海中央海盆下方存在快轴方向为NE-SW向的各向异性,其成因可能与海底扩张时期沿洋脊方向的地幔流以及南海海洋板块俯冲拖拽的地幔流有关。南海及其周边上地幔存在强各向异性,且不同方位观测到的各向异性不同,快轴方向与前人SKS横波分裂结果、GPS和板块运动一致,较好地对应了区域构造运动或者地幔对流模型。各向异性结果与印度—欧亚板块碰撞驱动挤出模型以及古南海俯冲板块拖拽模型预期结果一致,与理想的地幔柱上涌驱动模型不一致。由于海盆各向异性观测特别有限,各向异性结果不能证实亦不能证伪“大西洋型”海底扩张模型、弧后扩张模型和板缘破裂模型,后续还需要更多的观测结果来证实或证伪上述模型。
    Abstract: The South China Sea (SCS) is located at the intersection of the Eurasian, Pacific, and India-Australia plates. It is the largest marginal sea in a series of marginal seas in the Northwest Pacific Ocean. Many models have been proposed for the opening of the SCS, such as the extrusion model driven by the collision of the Indian plate and the Eurasian plate, and the slab pull model related to the subduction of the proto-SCS. This study aims to constrain the models of opening the SCS through the anisotropic structure of the central basin of the SCS and its surroundings. Based on the seismic data recorded by ten ocean bottom seismometers recovered from two passive seismic experiments conducted by Tongji University in the central basin of the SCS in 2012 and 2014, three different shear wave splitting methods are used to obtain the XKS splitting results of the central basin for two global earthquakes and the S phase splitting results provided by 20 regional earthquakes surrounding the SCS. The SKS splitting results demonstrate the presence of strong anisotropy with the NE fast direction in the central basin of the SCS, which may be related to mantle flow along the ocean ridge during seafloor expansion and the mantle flow dragged by the subduction of the proto-SCS plate. Strong anisotropy is also observed in the upper mantle surrounding SCS, and the anisotropy observed in different azimuths is different. The fast directions obtained are consistent with previous SKS-splitting results, GPS, and plate motions, and importantly correspond well to the regional tectonics or mantle convection models. The anisotropic results are consistent with the expected results of the extrusion model driven by the collision of the Indian-Eurasian plate and the slab pull of proto-SCS. The anisotropy results are inconsistent with the ideal upwelling driven model of the mantle plume. Unfortunately, due to the limited splitting observations in the central basin, the anisotropic results cannot confirm or falsify the “Atlantic-type” seafloor spreading model, the backarc spreading model, or the plate-edge rifting model. To verify the above models, further observations are needed.
  • 在均匀水平层状介质中,面波具有多模特征,即在相同频率下存在多个传播速度。传统上将这些传播速度从低到高排列,速度最低的模式称为基阶模式,也称为零阶模式,向上依次称为一阶、二阶等高阶模式(夏江海,2015)。对于速度递增的均匀水平层状介质,通常基阶模式能量占据主导地位,观测频散曲线通常对应基阶模式(Socco et al,2010)。在反演时,基于理论预测模型计算出频散曲线后,可以直接取速度最小的模式与观测频散曲线进行拟合。如果基于观测数据得到分离的高阶模式的频散曲线,也可以在反演时同时拟合高阶模式的频散曲线。只要将同一频率计算的速度由低到高排列,依次与基阶、一阶和二阶等模式频散曲线进行拟合即可。反演时考虑高阶模式瑞雷面波具有两个明显的优势:① 相同频率下,高阶模式面波比基阶模式具有更深的穿透深度,可以约束更深层的S波速度(Gabriels et al,1987Socco et al,2010);② 同时考虑基阶和高阶模式频散曲线可以提高反演的稳定性(Xu et al,2006)和准确性(Maraschini et al,2010)。

    然而在浅地表地球物理研究领域,由于沉积环境不同或经过人为改造或其它自然条件,会形成存在低速夹层的介质。根据频散曲线的形态和分布,可以将低速层存在时的频散曲线简单地分为两类。其中一类为含低速层模型的频散曲线,在低频段,基阶模式瑞雷面波在某个频率处表现出明显的低速特征;在高频段,频散曲线存在视觉上的交叉现象,观测到的频散曲线在各个模式之间跳跃,或许是不同模式共同作用的结果。Foti (2000)称此时的频散曲线为有效频散曲线(effective dispersion curve),凡友华等(2009)将其称之为能量最大模。此时,很难精确地将频散曲线与具体的模式对应。如果按照速度递增正频散介质的反演方法,无法根据相速度的大小选择合适的模式进行拟合,即存在模式误判(mode-misidentified)现象(Zhang,Chan,2003Gao et al,2014)。对于此类情形,可以采用Foti (2000)的方式,对预测的理论模型计算时域观测记录,采用与实际记录相同的处理方式,获得有效的相速度并与观测的有效相速度进行拟合,或者直接拟合观测与理论波场(Forbriger,2003ab)。通过理论和实验模拟,多名研究人员提出一种考虑高阶模式频散曲线的反演方法(鲁来玉等,2006;Lu,Zhang ,2006;Lu et al2007),该方法通过计算各个模式能量的大小,选取各频段能量最大的模式,与观测的频散曲线进行拟合。上述所列方法需要计算预测模型的波场记录,比计算频散曲线需要更多的计算时间,但相对于反演精度的提高,这种代价是值得的(Maraschini,Foti,2010Soccoet al,2010Ikeda et al,2012)。

    很多研究者开展了基于高阶瑞雷面波模式的反演研究(Xia et al,2003鲁来玉等,2006Luo et al,2007罗银河等,2008李建平,2018Matsuzawa,Yoshizawa,2019Wang et al,2019)。本文主要讨论另一类含低速层介质模型的高阶模式面波反演问题。对于这类模型,频散曲线在低频段与前述的低速层模型类似,基阶模式在某些频段表现出明显的低速特征;但在实际能提取到的频段范围,即我们感兴趣的范围内,各个模式的频散曲线仍然是明显分离的,视觉上不存在交叉现象。限于激发源和复杂介质模型的原因,实际观测到的基阶模式频散曲线并不包含明显低速特征的频段,在有效的频段内,观测频散曲线与速度递增介质的频散曲线的介质类似,正如下文第二节中所讲的实例。即便如此,通过考虑高阶模式的反演,依然可以基于观测频散曲线重建介质的低速特征,而仅考虑基阶模式反演时则无法重建介质的低速层特征。鉴于此,本文拟通过实例分析和数值模拟,研究高阶模式频散信息对此类模型低速层的约束作用。在案例分析中,我们选择浅层地震反射勘探中被视作干扰的面波信号作为处理对象,联合基阶和高阶面波频散曲线反演浅地表S波速度,为浅层地震反射资料解释提供额外的约束。因此,本文的研究目的是双重的:一是讨论高阶模式瑞雷面波的频散对低速层的约束;二是通过案例分析证明,浅层地震反射数据中的面波可以被重新利用,在不增加经济成本的情况下,通过多模式面波反演为反射资料解释提供独立的信息约束。

    本文采用的实例分析数据来自“天津活断层探测与地震危险性评价”项目实施的浅层地震反射勘探的数据(王辉,丁志峰,2006)。如图1所示,该地震测线(L02)近垂直于天津北断裂布设。沿测线布设TN4,TN5,TN6,TN7共四个钻孔,钻孔间距分别为50,65,50 m,孔深分别为200,100,100,200 m。地震勘探数据采用可控源车激发能量,扫描频率为20—120 Hz,扫描时间长度为8 s,探测方向由北西至南东,数据采集系统采用德国DMT公司生产的SUMMIT高分辨率遥测地震仪。测线总长度为2.0 km,共采集103个点,炮点每次沿测线移动18 m;1—92号炮点由120道检波器接收,93号炮点接收道数为116道,此后每移动一次炮点记录减少6道,到103号炮点时记录为56道;设置偏移距为24 m,道间距为3 m,采样间隔取0.5 ms,记录长度为1 024 ms,采用三次震源叠加以增加信噪比。

    图  1  跨天津北断裂的地震勘探测线及钻孔位置示意图(陈宇坤等,2013
    Figure  1.  Schematic diagram of seismic survey line and borehole location across the north fault of Tianjin (Chen et al,2013

    图2给出了沿测线第41个炮点的时域记录,可以看出:体波信号主要集中在前200 ms,200 ms之后可以看到较强的面波信号;由于面波传播速度较慢,在第64道记录之后,面波能量的主要到时在1 024 ms之后,超过了当初主要研究体波反射设置的记录长度。对各炮点的多道记录,采用二维傅里叶变换(Mcmechan,Yedlin,1981Gabriels et al,1987Alleyne,Cawley,1991Forchap,Schmid,1998Lu,Zhang,2004),将其变换到频率-波数域内。如图3a所示,在频率-波数域中,体波能量在整个扫频范围20—120 Hz内均有分布;面波能量主要集中在20—60 Hz之间,表现为斜率不同的近似直线,斜率的大小取决于波的传播速度。将体波剔除后的结果如图3b所示,可以明显观察到不同模式的瑞雷面波相速度,其中基阶、一阶及二阶频散能量均较为明显。

    图  2  跨天津北某活动断层浅层地震勘探数据的时间序列
    Figure  2.  Time series of the shallow seismic exploration data across the north fault of Tianjin
    图  3  图2所示时间序列在频率-波数域和频率-相速度域中的能量分布
    (a) 未剔除体波的结果;(b) 剔除体波后的结果;(c) 提取的三个模式的频散点
    Figure  3.  The energy distribution in the frequency-wavenumber and frequency-phase velocity domains for the time series shown in Fig.2
    (a) The result before removing body waves;(b) The result after removing body waves; (c) The extraction of surface-wave dispersion curves for three modes

    根据相速度vR、频率f和波数k之间的关系,可将频率-波数域变换到频率-相速度域。如图3c所示,在频率-相速度域中,面波能量的极大值点对应面波频散方程的解,即面波的特征值。根据能量极大值对应的频率和相速度的连线,可以提取观测的不同模式的相速度频散。如图3c中的圆点所示,这些频散曲线根据相速度由低到高依次称为基阶模式、一阶模式和二阶模式。可以看出,基阶、一阶和二阶模式的能量较强,三阶模式之后能量逐渐减弱。因此,我们在后面的反演中将主要考虑前三阶模式。

    图3c显示,基阶、一阶和二阶模式的频散曲线是相互分离的,不存在实际勘探中复杂地质条件下可能会出现的模式跳跃或“之”字形回折现象。因此,在反演时不考虑激发强度,只考虑振型相速度的大小。在相同频率下,通过相速度的大小判断该速度属于基阶、一阶或二阶等模式,从而与相应模式的观测频散曲线拟合。

    反演采用地震学计算机程序(computer programs in seismology,缩写为CPS)中的surf96 (Herrmann,2013),由于该程序的反演方法基于线性最小二乘原理,对初始模型的依赖性较高。因此我们根据早期面波勘探解释频散曲线时相速度与深度的对应关系经验公式:

    $$ {v}_{{\rm{S}}}=1.1{v}_{{\rm{R}}} \qquad {\textit{z}}=\frac{\lambda}{3} \text{,} $$ (1)

    式中,vS为深度z处的S波速度,vR为波长λ对应的相速度(Foti,2000)。根据频率、相速度、波长的关系,计算出提取的基阶模式频散点在此经验公式中对应深度的S波速度;然后依据钻孔岩土介质性质给出了此文反演的初始模型。由于瑞雷面波相速度对密度和纵波速度的敏感度远远不及S波速度(Xia et al,1999鲁来玉等,2011),本文仅反演S波速度。密度和P波速度根据场地条件给出,完整的初始模型参数如表1所示。

    表  1  初始模型参数和各模式频散曲线联合反演的S波速度
    Table  1.  The parameters of the initial model and the S-wave velocities inverted from multi-mode dispersion curves
    层序层厚/m深度/mρ/(g·cm−3vP/(m·s−1初始vS/(m·s−1 反演vS/(m·s−1
    11.01.01.60300152135.4
    21.02.01.73400156145.3
    31.03.01.77470163200.8
    41.04.01.79550165221.5
    51.05.01.79570170195.4
    61.06.01.79590180169.6
    71.07.01.79600185155.3
    81.08.01.79700190102.4
    92.010.01.79800200152.8
    10 3.013.01.79900210254.6
    11 4.017.01.801 100220259.3
    半空间$ \infty $$ \infty $1.801 200230250.1
    下载: 导出CSV 
    | 显示表格

    首先仅考虑基阶面波模式的反演,结果如图4a所示,S波速度模型随迭代次数的变化如图4a左所示,将30次迭代之后的结果作为最终模型。可观察到迭代进行到15次左右时已经基本收敛,反演得到的模型是一个速度递增的S波速度剖面。图4a右中反演模型的理论频散曲线很好地拟合了观测的基阶模式频散曲线(图中红色误差棒),而反演模型的理论频散曲线与观测的一阶和二阶模式频散曲线(灰色误差棒)没有拟合。

    图  4  基于多模式瑞雷面波频散曲线的反演(左)及其最终模型理论频散曲线与观测频散曲线的拟合(右)
    (a) 仅用基阶模式;(b) 基阶、一阶模式联合;(c) 基阶、一阶和二阶模式联合
    Figure  4.  The inversion based on multi-mode Rayleigh-wave dispersion curves (left) and the fitting between the theoretical and observed dispersion curves (right)
    (a) Using only the fundamental mode;(b) Using the fundamental and the first modes;(c) Using the fundamental,the first and the second modes

    图4b左是使用基阶和一阶模式频散曲线联合反演的结果,为了增加收敛速度,反演时将基阶模式反演的结果作为初始模型。与图4a相比,加入一阶模式之后,反演的S波速度剖面出现明显的低速层,深度在4—10 m之间。此时,如图4b右所示,该模型的频散曲线很好地拟合了观测的基阶(红色误差棒)和一阶模式(绿色误差棒)频散曲线;二阶观测频散曲线拟合情况较之前稍好,但仍有偏差。

    图4c左是使用三个模式的观测频散曲线联合反演的结果,使用的初始模型是基阶和一阶模式联合反演的结果,即图4b左中的反演结果。可以看出,基阶、一阶和二阶模式的联合反演结果与基阶、一阶模式的联合反演结果相近,但拟合效果更佳。如图4c右所示,反演模型的理论频散曲线很好地拟合了各阶模式的观测频散曲线。

    上述过程体现了高阶模式面波反演对反演多解性的约束。在只考虑基阶模式的情况下,得到的地下S波速度是递增的,加入高阶模式之后,在4—10 m左右出现明显的低速层。这一方面表明,此频带范围的基阶模式可能对5 m以下深度的约束不够;另一方面说明,在实际应用时,由于频散数据的限制,仅仅考虑基阶模式可能无法获得真实的S波速度。

    通常,对于介质存在低速夹层的情形,其频散曲线形态可简单分为两类。第一类为含低速层模型的频散曲线,基阶模式的瑞雷面波在某个反映低速层深度的频率附近表现出明显的低速特征;在高频段,频散曲线存在视觉上的交叉现象,观测的频散曲线,在各个模式之间跳跃或呈现“之”字形的回折。另一类低速层介质,即如图4c右模型的频散曲线所示,各个模式的频散曲线是明显分离的,视觉上不存在交叉现象。图4c右中,观测的多模式频散曲线呈现分离的状态,在有效的频带范围(23—40 Hz)内,并没有明显的低速指示(比如类似10 Hz附近基阶模式的相速度极小值点),其形态和速度递增的模型频散曲线类似。对于这样的频散曲线形态,图4c左的结果表明,考虑三个模式的最终反演结果具有明显的低速层。由于每次考虑更高一阶的模式反演时,都是利用前一次的反演结果作为初始模型,最终的低速层结果是否可信?反演结果是否与初始模型的选择有关?为了回答这些问题,并验证反演结果的可信度,我们将通过速度递增和含有低速层的两个模型进行数值模拟,以研究高阶模式对低速层的约束。

    反演结果会受到两类因素影响:一类是比较稳定的因素,如反演方法本身带来的影响等;另一类是不稳定的因素,如实际数据提取产生的误差、频带范围限制和初始模型等先验信息的影响。为了判定前述反演结果的可靠性,我们选取两个模型进行数值模拟研究。模型参数如表2所示,其中:模型1是由图4a中仅采用基阶模式反演得到的模型再按速度大小构建的正频散介质模型;模型2是含有低速层的模型,即图4c中的最终反演模型。两个模型的理论频散曲线如图5所示。将图5中的理论频散曲线作为观测数据,采用与第二节不同的初始模型、不同的频带范围并考虑不同模式的频散曲线,讨论这些因素对反演结果的影响。在实际应用中,频散曲线的提取误差总是存在的,这是面波频散曲线反演的一般性问题。这里我们不考虑对频散曲线进行加噪处理,原因在于:① 频散曲线的提取误差即使在实际应用中,通常是利用频率-速度域的能量分布给出的,也不是统计意义上的误差;② 在反演中,如果要考虑频散曲线提取误差的影响,可以通过多次测量的标准差作为误差,直接作为输入。因此,本文将聚焦于初始模型和观测频段以及高阶模式面波对低速层的约束。本节将分作三个部分,前两小节分别介绍模型1和模型2的情况,最后进行分析讨论。

    表  2  模型1和模型2的模型参数
    Table  2.  The parameters of models 1 and 2
    层序层厚/m深度/m模型1模型2
    $\mathrm{\rho }/ ( \mathrm{g}\cdot {\mathrm{c}\mathrm{m} }^{-3} ) $${v}_{{\rm{P}}}/ ( {\rm{m}}\cdot {{\rm{s}}}^{-1} ) $${{v} }_{{{\rm{S}}} }/ ( \mathrm{m}\cdot {\mathrm{s} }^{-1} ) $$\mathrm{\rho }/ ( \mathrm{g}\cdot {\mathrm{c}\mathrm{m} }^{-3} ) $${v}_{{\rm{P}}}/ ( {\rm{m}}\cdot {{\rm{s}}}^{-1} ) $${{v} }_{{{\rm{S}}} }/ ( \mathrm{m}\cdot {\mathrm{s} }^{-1} ) $
    1 1.0 1.0 1.581 7 249.1 126.5 1.605 3 266.2 135.4
    2 1.0 2.0 1.654 6 399.7 156.0 1.629 7 373.2 145.3
    3 1.0 3.0 1.728 3 556.2 192.9 1.724 1 580.1 200.8
    4 1.0 4.0 1.735 2 659.8 196.8 1.776 2 738.1 221.5
    5 1.0 5.0 1.739 5 646.0 199.2 1.732 7 654.9 195.4
    6 1.0 6.0 1.739 5 664.8 199.3 1.683 6 554.8 169.6
    7 1.0 7.0 1.739 7 653.0 199.4 1.653 0 505.1 155.3
    8 1.0 8.0 1.742 1 739.8 200.8 1.508 2 376.8 102.4
    9 2.0 10.0 1.754 8 832.9 208.2 1.647 2 611.3 152.8
    10 3.0 13.0 1.767 5 926.3 216.0 1.824 7 1 089.9 254.6
    11 4.0 17.0 1.780 3 1 120.5 224.1 1.831 0 1 297.7 259.3
    半空间 $ \infty $ $ \infty $ 1.792 9 1 210.7 232.3 1.818 5 1 303.7 250.1
    下载: 导出CSV 
    | 显示表格
    图  5  表2中两个模型的理论频散曲线
    红点表示的频率范围与本文实际数据采集得到的观测频散曲线的频带范围相同,如图3c所示(a) 模型1:速度递增模型;(b) 模型2:含有低速层的模型
    Figure  5.  The theoretical dispersion curves calculated from two models 1 and 2 listed in Table 2
    The red dots indicate the theoretical dispersion curves with the same frequency range as the observed dispersion curves shown in Fig. 3c (a) Model 1:A normal layered model with S-wave velocity increasing with depth; (b) Model 2:A model with a low-velocity layer

    图6是利用较宽频带范围的频散数据进行反演的结果,观测到的频散点是模型1正演计算得到的1—60 Hz (间隔为1 Hz)的数据。依旧采用表1中的初始模型参数,但区别于第二节在加入新的一阶进行反演时用前一次反演的最终结果更新初始模型以加快收敛速度的做法,本节在反演时均使用同一个初始模型。图6左列分别表示仅使用基阶(图6a)、基阶和一阶(图6b)以及基阶、一阶和二阶(图6c)模式反演的结果。灰色实线表示初始模型,红色到蓝色表示不同迭代次数的结果,蓝色实线是30次迭代最终得到的结果,黑色虚线表示理论模型;右列对应左侧同一行反演所得最终模型计算的理论频散曲线与所使用的观测频散点的拟合情况。

    图  6  基于模型1多模式瑞雷面波频散曲线的反演(左)及其最终模型理论频散曲线与观测频散曲线的拟合(右)
    (a) 仅用基阶模式;(b) 基阶和一阶模式联合;(c) 基阶、一阶和二阶模式联合
    Figure  6.  The inversion based on multi-mode Rayleigh-wave dispersion curves of model 1 (left) and the fitting between the theoretical and observed dispersion curves (right)
    (a) Using only the fundamental mode;(b) Using the fundamental and the first modes;(c) Using the fundamental,the first and the second modes

    对比图6中利用不同模式频散曲线的反演结果,可以看出,对于正频散的介质模型(模型1),如果观测频散曲线的频段较宽,仅利用基阶模式的信息就能很好地拟合理论频散曲线,而且仅利用基阶模式反演得到的模型,对高阶模式的观测频散曲线也有很好的拟合。与加入高阶信息反演的结果相对比,唯一能分辨的区别可能在于在反演迭代过程中收敛速度上的微小差异。对比图6a6b6c中的反演迭代次数可知,增加高阶模式频散信息联合反演迭代过程的收敛速度相对于仅使用基阶反演时稍快。

    图7是基于图5a中红色圆点所示的23—40 Hz之间的频散点进行反演的结果,其中基阶、一阶和二阶模式的频散曲线与实例分析中的观测频散曲线频带范围一致。图7中左列分别表示仅使用基阶(图7a)、基阶和一阶(图7b)以及基阶、一阶和二阶(图7c)模式反演的结果。灰色实线表示初始模型,红色到蓝色表示不同迭代次数的结果,蓝色实线是30次迭代最终得到的结果,黑色虚线表示的是理论模型;右列分别对应左列最终模型计算的频散曲线与所使用的频散点的拟合情况。

    图  7  反演中仅考虑较窄频段(23—40 Hz)的数据,其它同图6
    Figure  7.  The same as Fig. 6,but only narrow frequency band (23−40 Hz) data are considered in the inversion

    观察图7可以发现,即使观测频带范围较窄,对于速度递增的正频散介质(模型1),仅用基阶反演与加入高阶信息反演的结果都能很好地收敛于理论模型。但是,随着高阶模式信息的加入,反演的收敛速度明显加快。此外,与图6对比发现,仅使用基阶频散信息反演时,频带范围也会影响反演的收敛速度。图6a显示迭代小于5次就基本收敛于理论模型,图7a经过了更多次的迭代才达到与其相当的收敛效果。此时,加入高阶模式反演能缩小甚至消除这种差距(如图7b图7c)。上述结果表明,仅使用基阶频散数据进行反演时,提取的数据的频段范围会影响反演的收敛速度,加入高阶信息可以减小对反演收敛速度的影响。因此,在实际运用中,加入高阶模式信息可以加快反演收敛速度,使反演更加稳定。

    在第二节的实例分析中,为了加快收敛速度,我们通常采用低阶模式反演的最终模型作为考虑更高阶模式反演的初始模型。在图6图7中,仅考虑基阶模式的反演结果已相当接近理论的模型。三种反演情况的初始模型是相同的,考虑不同模式不同频带范围的反演结果都未出现低速层。因此这个结果说明,在实例分析中,加入高阶模式之后,反演得到的低速特征并非初始模型选择的原因,而是实际介质的真实反映,这一点从下文对含有低速层的模型2的模拟结果可得到进一步印证。

    图8图9图10给出了包含低速层的模型2的数值模拟结果。图8是基于较宽频带的数据反演的结果,左侧分别是利用基阶模式(图8a左)、基阶和一阶(图8b左)以及基阶、一阶和二阶(图8c左)模式反演的模型随迭代次数的变化,图中灰色实线表示初始模型,红色到蓝色实线代表不同迭代次数的结果,蓝色实线是30次迭代的最终结果,黑色虚线表示理论模型;图8c图8d的区别在于后者反演使用的初始模型为基阶和一阶的反演结果。图8右侧分别对应左侧最终模型计算的理论频散曲线与所使用的观测频散点的拟合情况。图9将二阶观测模式频散曲线截掉几个高频点并采用三种模式共同反演,图9b图9a的区别是前者将基阶和一阶的反演结果作为初始模型。图10图8类似,差别是图10各模式的频带范围较窄,在23—40 Hz之间。

    图  8  基于模型2多模式瑞雷面波频散曲线的反演(左)及其最终模型理论频散曲线与观测频散曲线的拟合 (右)
    (a) 仅用基阶模式;(b) 基阶和一阶模式联合
    Figure  8.  The inversion based on multi-mode Rayleigh-wave dispersion curves of model 2 (left) and the fitting between the theoretical and observed dispersion curves (right)
    (a) Using only the fundamental mode;(b) Using the fundamental and the first modes
    图  8  基于模型2多模式瑞雷面波频散曲线的反演(左)及其最终模型理论频散曲线与观测频散曲线的拟合 (右)
    (c) 基阶、一阶和二阶模式联合;(d) 基阶、一阶和二阶模式联合,其中初始模型为图(c)的最终模型
    Figure  8.  The inversion based on multimode Rayleigh-wave dispersion curves of model 2 (left) and the fitting between the theoretical and observed dispersion curves (right)
    (c) Using the fundamental,the first and the second modes;(d) Using the fundamental, the first and the second modes,and the initial model is the final model of Fig.c;
    图  9  基于模型2多模式瑞雷面波频散曲线的反演 (左) 及其最终模型理论频散曲线与观测频散曲线的拟合 (右)
    (a) 基阶、一阶和二阶模式联合(仅使用部分二阶频散数据);(b) 基阶、一阶和二阶模式联合,其中初始模型为图8b的最终模型(仅使用部分二阶频散数据)
    Figure  9.  The inversion based on multi-mode Rayleigh-wave dispersion curves of model 2 (left) and the fitting between the theoretical and observed dispersion curves (right)
    (a) Using the fundamental,the first higher and the second higher modes (only part of the second higher mode dispersion data are used);(b) Using the fundamental,the first and the second modes,and the initial model is the final model of Fig.8b (only part of the second mode dispersion data are used)
    图  10  反演中仅考虑较窄频段(23—40 Hz)的数据,其它同图8
    Figure  10.  The same as Fig. 8,but only narrow frequency band (23−40 Hz) data are considered in the inversion

    由于模型2的反演结果差异较大,因此只讨论结果好坏,不再讨论收敛速度。从图8中可观察到:① 仅利用基阶反演的结果在5 m内是很理想的,但在5—13 m与理论模型有所偏差,特别是在6—8 m和10—13 m处(图8a左),此时,理论频散曲线与基阶观测频散点拟合很好,但并没有拟合上(图8a右);② 加入一阶进行反演的结果明显好很多,所得结果与理论模型相当接近,理论频散曲线与基阶、一阶拟合很好,与二阶在相对低频部分拟合很好,高频部分仅有微小偏差(图8b左);③ 图8c左在加入二阶信息后,结果却与理论模型有了很大偏差;观测频散点与理论频散曲线拟合也不好(图8c右);④ 图8d左是利用基阶和一阶的反演结果作为初始模型,加入二阶信息反演的结果,对比前面三种情况,此时得到的结果是最好的,其频散点拟合情况也很好。

    图9,截去二阶54—60 Hz范围内的七个观测频散点,然后利用基阶、一阶和二阶联合反演的结果(图9a)相较图8b左的结果在2—10 m深度更好;而以基阶、一阶反演结果作为初始模型反演(图9b左)的结果较此还要更好一些,与图8d左的结果基本一致。

    图10a仅考虑基阶模式的反演结果显示,4—10 m的低速层未能很好地恢复,反演模型的频散曲线在10 Hz左右的极小值也不够明显,最终结果与真实模型尽管相差较大,但所表现出的高速层和低速层深度基本一致。随着加入高阶信息反演,得到的反演模型的低速层逐渐接近真实模型,且反演模型的基阶模式频散曲线中10 Hz处的极小值也越来越明显。

    图8a中的基阶模式频散曲线具有较宽的频带范围,在整个频段上,频散曲线的变化特征就指示了低速层的存在,比如图8b中基阶模式频散曲线在10 Hz左右具有速度的极小值。因此,仅利用基阶模式的反演结果,虽然与真实模型有一定的差距,但反演模型表现出明显的低速特征。正因为如此,在实际应用中,如果观测的频散曲线不含明显的低速特征,经验会使人们以为反演的模型不应具有明显的低速特征,这种经验并不总是有用的。如图10所示,如果观测的频散曲线只是整个频段的一部分,且不具有明显的低速特征,按照经验,其反演模型应该不具有明显的低速特征。但图10的结果表明,虽然仅考虑基阶模式时,反演模型频散曲线的低速层特征并不明显;但加入模式之后,反演模型频散曲线的低速特征相当接近真实模型,如图10b10c所示。这里对于仅考虑基阶和模式的反演,采用了同样的初始模型,且给定的初始模型不具有明显的低速特征,因此可以排除低速层的存在是初始模型造成的。对比图10b图10c,可以看出基阶、一阶联合反演结果与基阶、一阶及二阶的联合反演结果相比并不能确切地说哪个更好。具体来说,前者对5 m以上的层约束更好,后者对5 m以下的约束更好;且二者的结果均比仅使用基阶反演的结果要好。在图10d中,采用基阶和一阶反演的结果作为考虑更多模式反演时的初始模型,其反演结果与理论模型基本一致。图10的模拟结果表明,图4中反演得到的低速层特征,并不是初始模型选择的原因。处理实际数据时,用前一次的反演结果作为考虑更反演时的初始模型也是合理的。

    Liang等(2008)的研究认为低速层以下的相速度具有很低的敏感度。对于相速度对S波速度敏感度很低的地层,其反演结果具有较强的不确定性,需慎重对待(尹晓菲等,2020)。基于这个观点,我们计算了频率为30,35和40 Hz在不同深度的瑞雷面波相速度对S波速度的敏感度(图11)。整体而言,两个模型的基阶模式敏感度主要集中在浅地表4 m之内,4 m以下敏感度迅速衰减为0,这佐证了第二节我们推测实际提取的基阶频散数据对5 mm以下的深度约束不够的观点。对于速度递增的模型1,随着阶数的增加,相同频率下,阶数越高,其敏感度反映的深度越深,如图11所示,模型1的一阶模式敏感度主要集中在0—12 m以内,二阶模式在0—17 m均有敏感度。对含有低速层的模型2,基阶模式的敏感度特征与速度递增的模型1类似,但是高阶模式的敏感度表现出明显不同的特征,其敏感度的最大值集中在4—10 m,即模型2的低速层之内,这是典型的能陷模式的导波特征(张碧星,鲁来玉,2002),与包含低速层的多层介质高阶模式面波产生的机理相关,一阶和二阶模式对应在低速层内来回反射,能量陷于低速层内的能陷模式,因此一阶和二阶模式对低速层更为敏感,即使观测频散曲线不包含明显指示低速特征的频段,也可以通过一阶和二阶模式的加入反演约束模型的低速层。这也证明,图4的实例分析中反演得到的低速层,是实际介质模型的反映,与初始模型的选择无关。

    图  11  模型1 (a)和模型2 (b)在不同频率和不同模式下瑞雷波相速度的S波速度深度敏感核
    Figure  11.  The sensitivity analysis of the multimode Rayleigh-wave phase velocities at different frequencies for models 1 (a) and 2 (b)

    如Liang等(2008)所说,模型2的1阶高阶灵敏度在低速层4—10 m下确实逐渐衰减为0,但二阶在16 m深度以上仍有一定敏感度,且两种模式在低速层均有较高的敏感度。在4 m以内,基阶的灵敏度高出一阶和二阶很多。因此,在加入二阶反演时出现无法收敛的原因可能是由于反演具有多解性,基阶和高阶都对此范围有约束,反而影响了反演结果。加入二阶数据反演时出现的情况说明反演算法的局限性限制了数据利用的最大化,因此反演并非利用越多的观测数据越好,应根据实际情况使用数据。实际应用中,最可靠的标准应该是各阶观测频散数据的拟合情况,拟合好的情况下,加入高阶信息进行反演的结果明显更准确,即加入高阶信息可以提高反演准确性;而由于反演的多解性及反演算法对初始模型的依赖性,在观测数据拟合不好的情况下,不能盲目加入高阶信息,可以尝试对高阶信息先筛选后使用,或将未加入更高阶信息时拟合情况较好的反演结果作为初始模型,再进行反演从而得到更准确的结果。

    为了更直观地观察考虑不同模式的反演结果,图12给出了模型1和模型2考虑不同模式、基于不同频段频散曲线的反演结果。图12a为速度递增的模型1的反演结果,左侧为较窄频段的频散曲线反演结果,右侧为较宽频段的频散曲线反演的结果。整体而言,仅用基阶模式与加入高阶模式的反演结果差异较小,仅用基阶模式反演也能较好地恢复真实模型。

    图  12  模型1 (a)和模型2 (b)中采用不同模式基于较窄(左)和较宽(右)频段的频散曲线的反演结果
    Figure  12.  Inversion results of dispersion curves based on narrower (left) and wider (right) frequency bands using different modes for models 1 (a) and 2 (b)

    图12b为包含低速层的模型2的反演结果,左侧为窄频段的频散曲线反演结果,右侧为较宽频段的频散曲线反演的结果。观察图12b可以发现,仅使用基阶及使用基阶和一阶反演时,宽频带范围的结果明显更好;但加入二阶后,窄、宽频段的反演结果相差无几。这一点并不奇怪,在图12b右中,反演所用的频带包含了高频和低频段数据,相应地提供了更多地约束。因此,如果能获取感兴趣的深度范围对应频率范围完整的基阶频散数据,那此时仅用基阶反演得到的结果相对来说也是比较可靠的。另外,与模型1相比,此时高阶模式对低速层的约束更为显著,仅考虑基阶模式时,如果频带较窄,反演模型低速特征并不明显(图12b左灰色实线);如果频带范围较宽,反演的模型包含低速特征(图12b右灰色实线),但与真实模型也相差较远;对于窄频段情况,仅考虑基阶模式则难以恢复模型的低速特征(图12b左灰色实线)。这与之前众多学者研究多模式瑞雷面波频散反演得出的观点—高阶模式对地下介质分析有着重要作用,尤其是存在低速层时作用更为显著相契合(如Tokimatsu et al,1992)。

    通过以上数值模拟,我们论证了:① 对于本文所考虑的这类包含低速层的介质模型,仅从观测的基阶模式无法准确地重建介质的低速特征,高阶模式瑞雷面波对低速层具有重要的约束作用;② 相对于基阶模式,包含低速层的高阶模式对应着时域多径在低速夹层中来回反射、能量陷于低速层内部的能陷导波,因此其灵敏度主要集中在低速层区域;③ 在第二节中,由实际观测数据反演得到的低速层并不是初始模型选取的原因,而是实际介质的真实反应。

    通过理论分析,我们论证了基于实际观测的多模面波反演的可靠性。将第二节中单个阵列的反演结果作为阵列中心位置下方的速度结构,然后移动阵列,对各个炮点的数据进行类似处理,再将不同阵列的一维S波速度剖面拼合得到拟二维S波速度剖面。

    按照上述方法,我们对测线上103个炮点的排列都进行处理,得到了各阵列下方的一维S波速度结构,再通过插值构建出沿测线的拟二维S波速度剖面,如图13所示。由于缺乏该区的钻孔波速资料,我们只能结合该区的第四纪地层特征对所得结果进行定性分析。前人利用基于天津地区钻孔资料开展的古地磁与磁性地层的研究所得结果,并结合钻孔资料,给出的天津地区的第四纪地层剖面(如陈宇坤等,2008)显示,天津地区的第四纪地层全新统(Qh)底板埋深约为18.6 m。该深度范围正是我们反演的有效深度范围,顶部多有人工填土,底部常含有泥炭层和咸、淡水化合物,为海侵前的海相沉积。沿测线的钻孔TN7,TN6,TN5及TN4的结果显示,第一海相层的深度约在2—18 m之间(陈宇坤等,2013),其图层性质如图13右所示(陈宇坤等,2013)。反演的整个剖面显示,在浅地表2 m之内,有一个低速层,对应地表的杂填土;在2—5 m范围,对应第一海相层的顶层;其土层性质为黏土,相对地表的杂填土和下面的海相沉积,表现为高速特征,速度范围为200—300 m/s;在5—10 m范围出现了一个低速层,其速度范围在100—200 m/s,推测是粉砂和细粉砂海相沉积层的反映。在此低速层之后速度逐渐增加至12 m深度的200—300 m/s。

    图  13  S波速度剖面 (a)及钻孔土层 (b)
    Figure  13.  S-wave velocity profile (a) and borehole soil layer (b)

    在实际应用中,存在一类包含低速层的介质模型,在感兴趣的频率范围内,瑞雷面波各模式频散曲线并无视觉上的交叉现象,各个模式的观测频散曲线也无明显的低速特征指示。针对这类包含低速层的介质模型,本文通过对天津某地一条测线的浅层地震反射勘探资料进行二次处理提取出其中的面波信息,结合数值模拟,研究了高阶瑞雷面波对低速层的约束作用。研究结果表明:① 对于该类包含低速层的介质模型,仅仅由频散点数较少的基阶模式的频散曲线无法反演重建介质的低速特征,想要约束模型的低速层需要加入高阶模式的频散曲线;② 这类含低速层的模型,其观测频散曲线在感兴趣的频段可能无明显的低速特征指示,如果同时考虑高阶模式反演,反演模型的低速特征可能是真实模型的体现,并不是由初始模型的选择或其它原因引起;③ 这类模型的存在表明,观测的基阶模式无低速特征指示,并不意味着真实模型不含低速层,因此在实际应用中,仅由观测曲线的速度变化特征定性判定模型是否具有低速特征,具有一定局限性,尤其是观测频段较窄时;④ 验证了传统高阶模式反演的两个优势,即相同频率下,高阶模式面波比基阶模式具有更深的穿透深度,以及同时考虑基阶和高阶模式频散曲线可以加快收敛速度,增加反演的稳定性;⑤ 浅层地震反射中的勘探数据,经过一定处理,可以对之前被视为干扰的面波信号进行二次利用;且由于此类数据能提取到的频散数据频率覆盖范围较窄,此时高阶模式反演很适用。就本文考虑的实例而言,利用较少模式反演的结果作为增加更多模式反演时的初始模型是合理的选择,反演模型中低速层的出现并非出于初始模型选择。

    天津市地震局高武平高级工程师在资料处理和论文撰写过程中提供了帮助,两位审稿专家为文章的完善提出了宝贵意见,作者在此一并表示感谢!

  • 图  2   (a) 布设的海底地震仪及提供了有用的区域S波震相分裂解的20次区域地震的分布;(b) 中央海盆的地震仪分布;(c) 提供了有用的SKS震相分裂解的两次远震的分布

    Figure  2.   (a) The deployed ocean bottom seismometers and the distribution of regional earthquakes which provide useful regional S-wave splitting results;(b) The enlarged version of the distribution of OBSs in central basin;(c) The distribution of the two global earthquakes which provids useful SKS splitting results

    图  1   南海区域构造地质、地震台站及历史地震分布图

    南海的海陆过渡带位置参考Li和Song(2012),马尼拉海沟位置来自Jiang和Wang(2021);其中中南断裂和残留洋中脊的位置来自Zhao 等(2019),红河断裂带来自Yang 等 (2015)。箭头表示板块的绝对运动方向(Gripp,Gordon,2002

    Figure  1.   Tectonics of South China Sea and the distributions of seismic stations as well as history seismic events

    Continental Ocean Boundaryof South China Sea refers to Li and Song (2012),and the location of Manila trench is from Jiang and Wang (2021). The locations of Zhongnan fault and the remnant mid-ocean ridges are from Zhao et al2019) and the location of Red River Fault is from Yang et al2015). Black arrows indicate the absolute plate motions (Gripp,Gordon,2002

    图  3   (a) 利用Taup软件并采用IASP91速度模型得到的SKS,S,pS,sS和PS震相射线路径示意图;(b)文中提供有用横波分裂解的S震相的入射角度统计图

    Figure  3.   (a) Schematic raypaths of SKS,S,pS,sS,and PS phases obtained using the Taup package and the IASP91 velocity model (Kennett,Engdahl,1991Crotwell et al,1999);(b) Statistical histogram of the incident angles of the S phases used in shear-wave splitting

    图  4   采用切向能量最小化(左)和最小特征值最小化(右)方法对HY16台站记录的2012年11月7日危地马拉沿海地震激发的PKS震相进行横波分裂分析

    (a)—(c)横波分裂校正前后PKS震相滤波后的地震波形、投影到快/慢轴的波形及质点振动图(由原始的椭圆状线性化)

    Figure  4.   Shear wave splitting analysis on the PKS phase excited by the Guatemalan coastal earthquake on November 7,2012 recorded at station HY16 using the tangential energy minimization (left) and the minimum eigenvalue minimization (right) methods

    (a)−(c) The filtered seismic waveforms,waveforms projected onto the fast /slow axis and the particle motions which arelinearized from the original elliptical shape of the PKS phase before and after correcting anisotropy

    图  4   采用切向能量最小化(左)和最小特征值最小化(右)方法对HY16台站记录的2012年11月7日危地马拉沿海地震激发的PKS震相进行横波分裂分析

    (d)能量等值线图,能量被校正到最小的那个点即为最优解(星号)

    Figure  4.   Shear wave splitting analysis on the PKS phase excited by the Guatemalan coastal earthquake on November 7,2012 recorded at station HY16 using the tangential energy minimization (left) and the minimum eigenvalue minimization (right) methods

    (d) Contours of energy and the point where the energy corrected to the minimum is the optimal solution (asterisk)

    图  6   采用最小特征值最小化(左)和波形旋转互相关(右)方法分析HY08台站记录的2012年8月14日日本本州东部区域M7.7地震激发的区域S震相的横波分裂结果

    (a)—(c)为横波分裂校正前后PKS震相的地震波形、投影到快慢轴的波形及质点振动图(原始的椭圆状线性化),其中右侧图采用的是质点的位移而非质点的振动速度,下同

    Figure  6.   The S-phases are triggered by the M7.7 earthquake in the eastern region of Honshu,Japan on August 14,2012 recorded at station HY08,and a comparison is also made between the shear wave splitting methods of the minimum eigenvalue minimization (left) and the cross-correlation (right)

    (a)−(c) The filtered seismic waveforms,waveforms projected onto the fast /slow axis and the particle motions which are linearized from the original elliptical shape of the PKS phase before and after correcting anisotropy,figures in the right column use displacements rather than vibration velocities of the ground,the same below

    图  6   采用最小特征值最小化(左)和波形旋转互相关(右)方法分析HY08台站记录的2012年8月14日日本本州东部区域M7.7地震激发的区域S震相的横波分裂结果

    (d) 等值线图,能量和相关系数分别被校正到最小和最大的那个点即为最优解

    Figure  6.   The S-phases are triggered by the M7.7 earthquake in the eastern region of Honshu,Japan on August 14,2012 recorded at station HY08,and a comparison is also made between the shear wave splitting methods of the minimum eigenvalue minimization (left) and the cross-correlation (right)

    (d) Contour maps. The point where the energy and correlation coefficient are corrected to the minimum and maximum,respectively,is the optimal solution (asterisk)

    图  5   采用最小特征值最小化(a-b)和互相关横波分裂(c-e)方法分析HY17台站记录到的2012年11月7日危地马拉沿海地震激发的SKS震相得到的空解特征

    (a)表示SKS震相投影到最大及最小特征值方向的波形图;(b)为最小特征值方向能量的等值线图;(c)代表N-S和E-W向地震波位移记录;(d)表示质点运动图;(e)表示互相关方法得到的快慢波相关系数等值线图,两种方法的空解特征均为等值线图沿两个候选的快轴方向呈条带状分布

    Figure  5.   Examples of null results obtained through shear wave splitting analysis of the SKS phase excited by the November 7,2012 Guatemalan coastal earthquake recorded at HY17 station,and a comparison of the null results,obtained using two shear wave splitting methods of the minimum eigenvalue minimization (a-b) and the cross-correlation methods (c-e)

    (a) The SKS phases projected onto the directions of the maximum and minimum eigenvector;(b) The contours of the energy in the direction of the minimum eigenvector;(c) Seismic wave displacements in the N-S and E-W directions;(d) The particle motion diagram;(e) The contours of the cross-correlation coefficients between candidate fast and slow waves. Note that the common features of both methods for null results are that the contours are distributed in strips along the two candidate fast directions

    图  7   HY08地震−台站对采用最小特征值最小化方法的横波分裂结果进行聚类分析实例

    (a-b)为图6所示的在0.05—0.1 Hz滤波频段针对不同横波分裂时窗得到的快轴方向和分裂时间;(c-d)为针对不同滤波频段获得的最优快轴方向和分裂时间;(e-f)为0.05-0.1 Hz和不同滤波频段的聚类分析结果,其中相同颜色的解表示通过聚类分析后划分为同一类的解

    Figure  7.   An example of clustering analysis on shear wave splitting results obtained using the minimum eigenvalue minimization method for the earthquake-station HY08 pair

    (a-b) The fast directions and splitting times obtained at different splitting windows for the frequency band of 0.05−0.1 Hz in fig. 6;(c-d) The optimum fast directions and splitting times obtained at different filtering frequency bands;(e-f) The corresponding clustering analysis results for the frequency band of 0.05−0.1 Hz and for different frequency bands,where solutions with the same color represent solutions that are classified into the same class after clustering analysis

    图  8   南海中央海盆SKS分裂结果

    Figure  8.   The SKS splitting results in the eastern sub-basin of the South China Sea

    图  10   南海海盆及周边横波分裂及P波各向异性结果比较

    黄色双箭头表示本研究得到的SKS分裂解;蓝色和蓝绿色双箭头分别为采用互相关方法和最小特征值最小化方法针对区域地震激发的S震相得到的横波分裂解;针对SKS震相在不同台站获得单个分裂解以及台站叠加解来自文献(Bai et al,2009Hammond et al,2010Huang et al,2011Xue et al,2013Yu et al,2018Song et al,2021);P波各向异性结果来自Huang等(2015),其走向和长短分别代表快波的极化方向和各向异性的强度。台站东南方位马鲁古海区域根据横波分裂结果预期的该区域的地幔对流模型:黑色渐变粗实线代表Singihe板块下方与海沟平行的地幔流,板块两侧则是与板缘有关的偏转流,红色渐变粗实线则代表地幔楔里的角流,箭头指向流动方向

    Figure  10.   Comparisons of shear wave splitting results from this study with previous ones as well as the P-wave anisotropy results in the South China Sea Basin and its surrounding areas

    The yellow double arrows represent the SKS splits obtained in this study;the blue and blue-green double arrows represent the splits obtained using the cross correlation method and the minimum eigenvalue minimization method for the S-phase excited by regional earthquakes,respectively;the individual SKS splits and stacked SKS splits are from previous studies (Bai et al,2009Hammond et al,2010Huang et al,2011Xue et al,2013Yu et al,2018Song et al,2021);the P-wave anisotropy at different depths are from Huang et al (2015),and the direction and length of the short solid lines represent fast directions and anisotropic intensities respectively. The mantle convection model of the Molucca Sea region in the southeast of the OBSs consistent with most shear wave splitting results obtained in this study: the black gradient thick solid lines represent the trench-parallel flow under the Sangihe plate and the deflection flow related to the plate edge is on both sides of the plate,and the red gradient thick solid line represents the corner flow in the mantle wedge,and the arrows indicate the flow directions

    图  9   南海周边区域地震的S震相分裂结果

    蓝色短实线和蓝绿色短实线分别代表采用互相关方法和最小特征值最小化方法得到的横波分裂解

    Figure  9.   The S-wave splitting results for regional earthquakes surrounding the South China Sea

    Blue solid lines and light blue solid lines represent splitting results obtained by the cross-correlation method and the smallest eigenvalue minimization method,respectively. The lengths of the lines are proportional to their splitting time,as shown in the lower left corner of the figure

    表  1   提供有用数据的OBS信息

    Table  1   The information of OBSs used

    台站东经/°北纬/°水深/m记录起止日期钟漂/s水平分量方位角/°倾斜/°倾斜方向/°
    (年-月-日—年-月-日)
    HY02116.306 116.203 33 7502012-04-21—2012-12-0386.53203.5160
    HY08117.796 513.800 44 104 2012-04-22—2012-08-1902180.7285
    HY10116.998 314.598 94 276 2012-04-23—2012-12-09193.21660.8215
    HY15117.537 016.503 33 753 2012-04-25—2012-12-190.7415416.3 5
    HY16118.213 416.451 33 920 2012-04-25—2012-12-13−0.51260.3 91
    HY17118.803 716.201 03 870 2012-04-25—2012-12-10−10.543302.0131
    B04117.011 014.024 34 246 2014-07-11—2014-11-2301856.5280
    B19116.499 014.500 64 301 2014-07-11—2014-11-2303506.5110
    B32117.999 014.396 53 8592014-07-13—2014-11-2401100N/A
    注:台站钟漂的误差在2—3 s之间,水平分量方位角误差约为8°;倾斜指OBS的垂直分量偏离垂向的角度(刘晨光等,2014Le et al,2018Hung et al,2021
    下载: 导出CSV

    表  2   挑选后可用于横波分裂的地震事件信息

    Table  2   The information of selected events used for shear wave splitting

    事件序号年-月-日儒略日发震时刻纬度经度震源深度/kmMW注释
    12012-08-27240043719.4312.139°N88.590°W28.07.3远震
    22012-11-07312163546.6913.963°N91.854°W24.07.4远震
    32012-04-29120015751.903.059°S136.109°E18.45.2区域震
    42012-05-23144150225.3141.335°N142.082°E46.05.9区域震
    52012-06-01153065620.240.720°S133.269°E25.05.8区域震
    62012-06-29181210733.8643.433°N84.700°E18.06.3区域震
    72012-07-16198163310.081.296°S137.053°E13.15.6区域震
    82012-07-25207002745.262.707°N96.045°E22.06.4区域震
    92012-08-12225104706.4535.661°N82.518°E13.06.2区域震
    102012-08-14227025938.4649.800°N145.064°E583.27.7区域震
    112012-09-03247182305.2310.708°S113.931°E14.06.1区域震
    122012-09-07251031942.5327.575°N103.983°E10.05.5区域震
    132012-10-08282114331.424.472°S129.129°E10.06.1区域震
    142012-10-17291044230.404.232°N124.520°E32.66.0区域震
    152012-11-27332025906.522.952°S129.219°E11.25.7区域震
    16 2014-07-11192192200.8237.005°N142.452°E20.06.5区域震
    17 2014-07-20201183247.7944.642°N148.784°E61.06.2区域震
    18 2014-08-03215083013.5727.189°N103.409°E12.06.2区域震
    19 2014-08-06218114522.687.274°S128.036°E10.06.2区域震
    20 2014-09-10253051653.2118.400°N125.125°E30.05.9区域震
    21 2014-09-17260061445.4113.764°N144.429°E130.06.7区域震
    222014-11-18322044716.631.869°N126.475°E30.05.8区域震
    下载: 导出CSV

    表  3   远震的横波分裂结果

    Table  3   The shear-wave splitting results for teleseismic earthquakes

    台站事件后方位角/°震中距/°震相ϕδt/s方法
    HY02146.2142.6SKKS46.2;−43.8空解切向能量
    46.2;−43.8空解最小特征值
    46.2;−43.8空解互相关
    HY10246.5139.7PKS46.5;−43.5空解切向能量
    46.5;−43.5空解最小特征值
    46.5;−43.5空解互相关
    HY15245.5138.1SKKS45.5;−44.5空解切向能量
    45.5;−44.5空解最小特征值
    45.5;−44.5空解互相关
    PKS45.5;−44.5空解切向能量
    45.5;−44.5空解最小特征值
    45.5;−44.5空解互相关
    HY16246.2137.6PKS35±1.52.6±0.35切向能量
    33±6.52.25±0.55最小特征值
    46.2;−43.8空解互相关
    SKKS46.2;−43.8空解切向能量
    46.2;−43.8空解最小特征值
    46.2;−43.8空解互相关
    HY17247.1137.3SKS47.1;−42.9空解切向能量
    空解最小特征值
    空解互相关
    空解切向能量
    PKS47.1;−42.9空解最小特征值
    空解互相关
    下载: 导出CSV

    表  4   区域性地震的横波分裂结果

    Table  4   The shear-wave splitting results obtained for regional earthquakes

    台站事件后方位角/°震中距/°震相φ±σφ)/°δt±σδt)/s方法
    HY026321.938.2S,pS,sS−42±17.09.75±0.96互相关
    9309.235.8S,pS,sS,PS28±4.52.8±0.4最小特征值
    11185.226.8S,pS,sS−47±13.00.75±0.23最小特征值
    −48±24.00.75±0.60互相关
    12316.816.1S,sS,Sn,SS−1±8.51.8±0.375最小特征值
    −36±25.00.75±0.42互相关
    HY083131.424.7S,pS,sS,PS−47±52.72±0.5最小特征值
    −11±211.2±0.50互相关
    433.034.6S50±72.8±0.575最小特征值
    11±221.15±0.39互相关
    5132.021.1S,pS,sS−34±71.0±0.23最小特征值
    −48;42空解互相关
    6322.641.0S,pS,sS,PS−46±5.09.6±0.375最小特征值
    −49±8.09.85 ±0.45互相关
    7126.724.3S,pS,sS,PS27±15.52.65±0.83最小特征值
    59±24.01.5±0.55互相关
    1026.242.3S4±6.52.6±0.2最小特征值
    8±72.35±0.17互相关
    HY156320.838.7S,pS,sS,PS−38±4.59.6±0.28最小特征值
    −46±7.09.75 ±0.41互相关
    8239.325.2S,pS,sS,PS69±9.56.6±1.35最小特征值
    58±6.07.4±0.56s互相关
    1027.840.0S36±4.07.65±0.35最小特征值
    −42;48空解 互相关
    HY166320.439.1S,pS,sS−64±4.08.85±0.25最小特征值
    −63±3.09.0 ±0.14互相关
    11189.227.3S,pS,sS,PS−79±4.52.85±0.7最小特征值
    −69±24.65±0.24互相关
    13151.623.4S,pS,sS,PS20±9.51.9±0.275最小特征值
    17±181.75±0.27互相关
    14152.313.6S,Sn14±1.52.2±0.2最小特征值
    −67;23空解互相关
    15149.622.1S,pS,sS,PS42±5.03.25±0.85 最小特征值
    54±3.04.9±0.43互相关
    HY176320.239.8S,pS,sS,PS−44±5.08.80±0.23最小特征值
    −45±5.08.95 ±0.22互相关
    B041640.232.2S,pS,sS,PS−67±9.03.55±0.775最小特征值
    −80±114.25±0.65互相关
    22141.415.3S,Sn,sS−79±202.05±0.61互相关
    B191735.840.6S,pS,sS,PS−34±5.53.1±0.15最小特征值
    18318.017.6S,Sn,sS,SS−7±2.03.65±0.3最小特征值
    19151.424.5S,pS,sS,PS−41±1.58.95±0.25最小特征值
    −41±29.1±0.18互相关
    20148.916.9S,Sn,sS,SS−36±38.05±1.0最小特征值
    −41;49空解互相关
    2188.127.1S,sS5±3.54.1±0.325最小特征值
    −12±9.04.85±0.35互相关
    B321734.739.9S,sS−26±2.03.65±0.2最小特征值
    19154.523.7S,pS,sS,PS−75±2.53.7±0.4最小特征值
    −75;15空解
    2188.125.6S,sS7±52.95±0.2最小特征值
    2±102.95±0.4
    下载: 导出CSV
  • 陈立,薛梅,Phon L K,杨挺. 2012. 南海瑞雷面波群速度层析成像及其地球动力学意义[J]. 地震学报,34(6):754–772. doi: 10.3969/j.issn.0253-3782.2012.06.003

    Chen L,Xue M,Phon L K,Yang T. 2012. Group velocity tomography of Rayleigh waves in South China Sea and its geodynamic implications[J]. Acta Seismologica Sinica,34(6):754–772 (in Chinese).

    李家彪,金翔龙,高金耀. 2002. 南海东部海盆晚期扩张的构造地貌研究[J]. 中国科学:D辑,32(3):239–249.

    Li J B,Jin X L,Gao J Y. 2002. A study on the tectonic geomorphology of the late expansion of the Eastern sub-basin of South China Sea[J]. Science in China:Series D,32(3):239–249 (in Chinese).

    李家彪,金翔龙,阮爱国,吴世敏,吴自银,刘建华. 2004. 马尼拉海沟增生楔中段的挤入构造[J]. 科学通报,49(10):1000–1008.

    Li J B,Jin X L,Ruan A G,Wu S M,Wu Z Y,Liu J H. 2004. Intrusion structures in the middle of the accretionary wedge of the Manila Trench[J]. Chinese Science Bulletin,49(10):1000–1008 (in Chinese).

    李琳. 2016. 利用OBS数据反演南海海盆的各向异性结构[D]. 上海: 同济大学: 68.

    Li L. 2016. The Anisotropic Structure of South China Sea Constrained by OBS Data Geophysics[D]. Shanghai: Tongji University: 68 (in Chinese).

    李思田,林畅松,张启明,杨士恭,吴培康. 1998. 南海北部大陆边缘盆地幕式裂陷的动力过程及10Ma以来的构造事件[J]. 科学通报,43(8):797–810.

    Li S T,Lin C S,Zhang Q M,Yang S G,Wu P K. 1998. Dynamic process of episodic rifting in the northern continental margin basin of the South China Sea and tectonic events since 10 Ma[J]. Science Bulletin,43(8):797–810 (in Chinese).

    刘晨光,华清峰,裴彦良,杨挺,夏少红,薛梅,黎伯孟,霍达,刘芳,黄海波. 2014. 南海海底天然地震台阵观测实验及其数据质量分析[J]. 科学通报,59(16):1542–1552.

    Liu C G,Hua Q F,Pei Y L,Yang T,Xia S H,Xue M,Le B M,Huo D,Liu F,Huang H B. 2014. Passive-source ocean bottom seismograph (OBS) array experiment in South China Sea and data quality analyses[J]. Chinese Science Bull,59(16):1542–1552 (in Chinese).

    任建业,李思田. 2000. 西太平洋边缘海盆地的扩张过程和动力学背景[J]. 地学前缘,7(3):203–213.

    Ren J Y,Li S T. 2000. Spreading and dynamic setting of marginal basins of the western Pacific[J]. Earth Science Frontiers,7(3):203–213 (in Chinese).

    谢建华,夏斌,张宴华,王喜臣. 2005. 南海形成演化探究[J]. 海洋科学进展,23(2):212–218.

    Xie J H,Xia B,Zhang Y H,Wang X C. 2005. Study on formation and evolution of the South China Sea[J]. Advances in Marine Science,23(2):212–218 (in Chinese).

    胥颐,李志伟,郝天珧,刘建华,刘劲松. 2007. 南海东北部及其邻近地区的Pn波速度结构与各向异性[J]. 地球物理学报,50(5):1473–1479.

    Xu Y,Li Z W,Hao T Y,Liu J H,Liu J S. 2007. Pn wave velocity and anisotropy in the northeastern South China Sea and adjacent region[J]. Chinese Journal of Geophysics,50(5):1473–1479 (in Chinese).

    薛梅. 2008. 横波分裂方法研究[C]//中国地球物理学会第二十四届年会. 北京: 中国大地出版社: 390.

    Xue M. 2008. Study on shear-wave splitting [C]//The 24th Annual Meeting of the Chinese Geophysical Society. Beijing: China Land Press: 390 (in Chinese).

    鄢全树,石学法. 2007. 海南地幔柱与南海形成演化[J]. 高校地质学报,13(2):311–322.

    Yan Q S,Shi X F. 2007. Hainan mantle plume and the formation and evolution of the South China Sea[J]. Geological Journal of China Universities,13(2):311–322 (in Chinese).

    张健,熊亮萍,汪集旸. 2001. 南海深部地球动力学特征及其演化机制[J]. 地球物理学报,44(5):602–610.

    Zhang J,Xiong L P,Wang J Y. 2001. Characteristics and mechanism of geodynamic evolution of the South China Sea[J]. Chinese Journal of Geophysics,44(5):602–610 (in Chinese).

    张亮. 2012. 南海构造演化模式及其数值模拟[D]. 青岛: 中国科学院研究生院(海洋研究所):1−185.

    Zhang L. 2012. Tectonic Evolution of South China Sea and A Numerical Modeling[D].Qingdao: Graduate University of Chinese Academy of Sciences (Institute of Oceanography):1−185 (in Chinese).

    钟广见,林珍,高红芳,金华峰. 2006. 南海南北缘的构造特征对比[J]. 南海地质研究,(1):30–40.

    Zhong G J,Lin Z,Gao H F,Jin H F. 2006. Comparison of the tectonic characteristic between north margin and south margin of the South China Sea[J]. Geological Research of South China Sea,(1):30–40 (in Chinese).

    周蒂,陈汉宗,吴世敏,俞何兴. 2002. 南海的右行陆缘裂解成因[J]. 地质学报,76(2):180–190.

    Zhou D,Chen H Z,Wu S M,Yu H X. 2002. Opening of the South China Sea by dextral splitting of the East Asian continental margin[J]. Acta Geologica Sinica,76(2):180–190 (in Chinese).

    Bai L,Iidaka T,Kawakatsu H,Morita Y,Dzung N Q. 2009. Upper mantle anisotropy beneath Indochina block and adjacent regions from shear-wave splitting analysis of Vietnam broadband seismograph array data[J]. Phys Earth Planet Inter,176(1/2):33–43.

    Bell S W,Forsyth D W,Ruan Y Y. 2015. Removing noise from the vertical component records of ocean-bottom seismometers:Results from year one of the cascadia initiative[J]. Bull Seismol Soc Am,105(1):300–313. doi: 10.1785/0120140054

    Ben-Avraham Z,Uyeda S. 1973. The evolution of the China Basin and the mesozoic paleogeography of Borneo[J]. Earth Planet Sci Lett,18(2):365–376. doi: 10.1016/0012-821X(73)90077-0

    Bowman J R,Ando M. 1987. Shear-wave splitting in the upper-mantle wedge above the Tonga subduction zone[J]. Geophys J Roy Astron Soc,88(1):25–41. doi: 10.1111/j.1365-246X.1987.tb01367.x

    Briais A,Patriat P,Tapponnier P. 1993. Updated interpretation of magnetic anomalies and seafloor spreading stages in the South China Sea:Implications for the tertiary tectonics of Southeast Asia[J]. J Geophys Res:Solid Earth,98(B4):6299–6328. doi: 10.1029/92JB02280

    Cao L M,He X B,Zhao L,Lü C C,Hao T Y,Zhao M H,Qiu X L. 2021. Mantle flow patterns beneath the junction of multiple subduction systems between the Pacific and Tethys Domains,SE Asia:Constraints from SKS-wave splitting measurements[J]. Geochem Geophys Geosyst,22(9):e2021GC009700.

    Chevrot S. 2000. Multichannel analysis of shear wave splitting[J]. J Geophys Res:Solid Earth,105(B9):21579–21590. doi: 10.1029/2000JB900199

    Clift P,Lee G H,Anh Duc N,Barckhausen U,Van Long H,Zhen S. 2008. Seismic reflection evidence for a Dangerous Grounds miniplate:No extrusion origin for the South China Sea[J]. Tectonics,27(3):TC3008.

    Crampin S,Gao Y. 2006. A review of techniques for measuring shear-wave splitting above small earthquakes[J]. Phys Earth Planet Inter,159(1/2):1–14.

    Crotwell H P,Owens T J,Ritsema J. 1999. The Taup ToolKit:Flexible seismic travel-time and ray-path utilities[J]. Seismo Res Lett,70(2):154–160. doi: 10.1785/gssrl.70.2.154

    Di Leo J F,Wookey J,Hammond J O S,Kendall J M,Kaneshima S,Inoue H,Yamashina T,Harjadi P. 2012a. Mantle flow in regions of complex tectonics:Insights from Indonesia[J]. Geochem Geophys Geosyst,13(12):Q12008.

    Di Leo J F,Wookey J,Hammond J O S,Kendall J M,Kaneshima S,Inoue H,Yamashina T,Harjadi P. 2012b. Deformation and mantle flow beneath the Sangihe subduction zone from seismic anisotropy[J]. Phys Earth Planet Inter,194-195:38–54. doi: 10.1016/j.pepi.2012.01.008

    Dziewonski A M,Chou T A,Woodhouse J H. 1981. Determination of earthquake source parameters from waveform data for studies of global and regional seismicity[J]. Journal of Geophysical Research:Solid Earth,86(B4):2825–2852. doi: 10.1029/JB086iB04p02825

    Eakin C M,Long M D,Scire A,Beck S L,Wagner L S,Zandt G,Tavera H. 2016. Internal deformation of the subducted Nazca slab inferred from seismic anisotropy[J]. Nat Geosci,9(1):56–59. doi: 10.1038/ngeo2592

    England P,Houseman G. 1986. Finite strain calculations of continental deformation 2:Comparison with the India-Asia collision zone[J]. J Geophys Res:Solid Earth,91(B3):3664–3676. doi: 10.1029/JB091iB03p03664

    England P,Molnar P. 1990. Right-lateral shear and rotation as the explanation for strike-slip faulting in eastern Tibet[J]. Nature,344(6262):140–142. doi: 10.1038/344140a0

    Gripp A E,Gordon R G. 2002. Young tracks of hotspots and current plate velocities[J]. Geophys J Int,150(2):321–361. doi: 10.1046/j.1365-246X.2002.01627.x

    Hall R,Spakman W. 2015. Mantle structure and tectonic history of SE Asia[J]. Tectonophysics,658:14–45. doi: 10.1016/j.tecto.2015.07.003

    Hammond J O S,Wookey J,Kaneshima S,Inoue H,Yamashina T,Harjadi P. 2010. Systematic variation in anisotropy beneath the mantle wedge in the Java-Sumatra subduction system from shear-wave splitting[J]. Phys Earth Planet Inter,178(3/4):189–201.

    Holtzman B K,Kohlstedt D L,Zimmerman M E,Heidelbach F,Hiraga T,Hustoft J. 2003. Melt segregation and strain partitioning:Implications for seismic anisotropy and mantle flow[J]. Science,301(5637):1227–1230. doi: 10.1126/science.1087132

    Huang Z C,Wang L S,Zhao D P,Mi N,Xu M J. 2011. Seismic anisotropy and mantle dynamics beneath China[J]. Earth Planet Sci Lett,306(1/2):105–117.

    Huang Z C,Zhao D P,Wang L S. 2015. P wave tomography and anisotropy beneath Southeast Asia:Insight into mantle dynamics[J]. J Geophys Res:Solid Earth,120(7):5154–5174. doi: 10.1002/2015JB012098

    Huang C Y,Wang P X,Yu M M,You C F,Liu C S,Zhao X X,Shao L,Zhong G F,Yumul G P Jr. 2019. Potential role of strike-slip faults in opening up the South China Sea[J]. Natl Sci Rev,6(5):891–901. doi: 10.1093/nsr/nwz119

    Hung T D,Yang T,Le B M,Yu Y Q. 2019. Effects of failure of the ocean-bottom seismograph leveling system on receiver function analysis[J]. Seismol Res Lett,90(3):1191–1199. doi: 10.1785/0220180276

    Hung T D,Yang T,Le B M,Yu Y Q,Xue M,Liu B H,Liu C G,Wang J,Pan M H,Huong P T,Liu F,Morgan J P. 2021. Crustal structure across the extinct mid-ocean ridge in South China Sea from OBS receiver functions:Insights into the spreading rate and magma supply prior to the ridge cessation[J]. Geophys Res Lett,48(3):e2020GL089755.

    Jiang X,Wang Z C. 2021. A critical review of existing models for the origin of the South China Sea and a new proposed model[J]. J Asian Earth Sci:X,6:100065.

    Karig D E. 1971. Origin and development of marginal basins in the Western Pacific[J]. J Geophys Res,76(11):2542–2561. doi: 10.1029/JB076i011p02542

    Kennett B L N,Engdahl E R. 1991. Traveltimes for global earthquake location and phase identification[J]. Geophys J Int,105(2):429–465. doi: 10.1111/j.1365-246X.1991.tb06724.x

    Kuo B Y,Chen C C,Shin T C. 1994. Split S waveforms observed in northern Taiwan:Implications for crustal anisotropy[J]. Geophys Res Lett,21(14):1491–1494. doi: 10.1029/94GL01254

    Le B M,Yang T,Chen Y J,Yao H J. 2018. Correction of OBS clock errors using Scholte waves retrieved from cross-correlating hydrophone recordings[J]. Geophys J Int,212(2):891–899. doi: 10.1093/gji/ggx449

    Li C F,Song T R. 2012. Magnetic recording of the Cenozoic oceanic crustal accretion and evolution of the South China Sea basin[J]. China Science Bulletin,57(24):3165–3181. doi: 10.1007/s11434-012-5063-9

    Li L, Xue M, Yang T, Liu C, Hua Q, Le B M, Xia S, Huang H, Pan M, Huo D. 2016. A comparison of various methods of determining the horizontal orientation of Ocean Bottom Seismometers (OBS): Applying to South China Sea OBS data[C]//AGU Fall Meeting. San Francisco: American Geophysical Union: 2016AGUFMOS51C2090L

    Long M D,Silver P G. 2008. The subduction zone flow field from seismic anisotropy:A global view[J]. Science,319(5861):315–318. doi: 10.1126/science.1150809

    Long M D,Becker T W. 2010. Mantle dynamics and seismic anisotropy[J]. Earth Planet Sci Lett,297(3-4):341–354. doi: 10.1016/j.jpgl.2010.06.036

    Mainprice D,Tommasi A,Couvy H,Cordier P,Frost D J. 2005. Pressure sensitivity of olivine slip systems and seismic anisotropy of Earth’s upper mantle[J]. Nature,433(7027):731–733. doi: 10.1038/nature03266

    Miyashiro A. 1986. Hot regions and the origin of marginal basins in the western Pacific[J]. Tectonophysics,122(3/4):195–216.

    Morley C K. 2002. A tectonic model for the Tertiary evolution of strike–slip faults and rift basins in SE Asia[J]. Tectonophysics,347(4):189–215. doi: 10.1016/S0040-1951(02)00061-6

    Rangin C,Klein M,Roques D,Le Pichon X,Van Trong L. 1995. The Red River fault system in the Tonkin Gulf,Vietnam[J]. Tectonophysics,243(3/4):209–222.

    Savage M K. 1999. Seismic anisotropy and mantle deformation:What have we learned from shear wave splitting?[J]. Rev Geophys,37(1):65–106. doi: 10.1029/98RG02075

    Shen Z K,Zhao C K,Yin A,Li Y X,Jackson D D,Fang P,Dong D N. 2000. Contemporary crustal deformation in East Asia constrained by global positioning system measurements[J]. J Geophys Res:Solid Earth,105(B3):5721–5734. doi: 10.1029/1999JB900391

    Silver P G,Chan W W. 1988. Implications for continental structure and evolution from seismic anisotropy[J]. Nature,335(6185):34–39. doi: 10.1038/335034a0

    Silver P G,Chan W W. 1991. Shear wave splitting and subcontinental mantle deformation[J]. J Geophys Res:Solid Earth,96(B10):16429–16454. doi: 10.1029/91JB00899

    Silver P G. 1996. Seismic anisotropy beneath the continents:Probing the depths of geology[J]. Annu Rev Earth Planet Sci,24:385–432. doi: 10.1146/annurev.earth.24.1.385

    Song W K,Yu Y Q,Gao S S,Liu K H,Fu Y F. 2021. Seismic anisotropy and mantle deformation beneath the central Sunda Plate[J]. J Geophys Res:Solid Earth,126(3):e2020JB021259.

    Sun Z,Lin J,Qiu N,Jian Z M,Wang P X,Pang X,Zheng J Y,Zhu B D. 2019. The role of magmatism in the thinning and breakup of the South China Sea continental margin[J]. Natl Sci Rev,6(5):871–876. doi: 10.1093/nsr/nwz116

    Tapponnier P,Molnar P. 1976. Slip-line field theory and large-scale continental tectonics[J]. Nature,264(5584):319–324. doi: 10.1038/264319a0

    Tapponnier P,Peltzer G,Le Dain A Y,Armijo R,Cobbold P. 1982. Propagating extrusion tectonics in Asia:New insights from simple experiments with plasticine[J]. Geology,10(12):611–616. doi: 10.1130/0091-7613(1982)10<611:PETIAN>2.0.CO;2

    Tapponnier P, Peltzer G, Armijo R. 1986. On the mechanics of the collision between India and Asia[C]//Collision Tectonics. London: Geological Society, Special Publications, 19: 113–157.

    Tapponnier P,Lacassin R,Leloup P H,Schärer U,Zhong D L,Wu H W,Liu X H,Ji S H,Zhang L S,Zhong J Y. 1990. The Ailao Shan/Red River metamorphic belt:Tertiary left-lateral shear between Indochina and South China[J]. Nature,343(6257):431–437. doi: 10.1038/343431a0

    Taylor B, Hayes D E. 1983. Origin and history of the South China Sea basin[C]//The Tectonic and Geologic Evolution of Southeast Asian Seas and Islands: Part 2. Washington D C: AGU: 27: 23–56.

    Teanby N A,Kendall J M,van der Baan M. 2004. Automation of shear-wave splitting measurements using cluster analysis[J]. Bull Seismol Soc Am,94(2):453–463. doi: 10.1785/0120030123

    Wang C Y,Yang J M,Zhu W L,Zheng M L. 1995. Some problems in understanding basin evolution[J]. Earth Sci Front,2(3/4):29–44.

    Wang L L,He X B. 2020. Seismic anisotropy in the Java-Banda and philippine subduction zones and its implications for the mantle flow system beneath the Sunda Plate[J]. Geochem Geophys Geosyst,21(4):e2019GC008658.

    Wang P X,Huang C Y,Lin J,Jian Z M,Sun Z,Zhao M H. 2019. The South China Sea is not a mini-Atlantic:Plate-edge rifting vs intra-plate rifting[J]. Natl Sci Rev,6(5):902–913. doi: 10.1093/nsr/nwz135

    Wu J,Suppe J. 2018. Proto-South China Sea plate tectonics using subducted slab constraints from tomography[J]. J Earth Sci,29(6):1304–1318. doi: 10.1007/s12583-017-0813-x

    Xue M,Allen R. 2005. Asthenospheric channeling of the Icelandic upwelling:Evidence from seismic anisotropy[J]. Earth Planet Sci Lett,235(1/2):167–182.

    Xue M,Le K P,Yang T. 2013. Seismic anisotropy surrounding South China Sea and its geodynamic implications[J]. Mar Geophys Res,34(3/4):407–429.

    Yang T,Liu F,Harmon N,Le K P,Gu S Y,Xue M. 2015. Lithospheric structure beneath Indochina block from Rayleigh wave phase velocity tomography[J]. Geophys J Int,200(3):1582–1595. doi: 10.1093/gji/ggu488

    Yu Y Q,Gao S S,Liu K H,Yang T,Xue M,Le K P,Gao J Y. 2018. Characteristics of the mantle flow system beneath the Indochina Peninsula revealed by teleseismic shear wave splitting analysis[J]. Geochem Geophys Geosyst,19(5):1519–1532. doi: 10.1029/2018GC007474

    Zhao M H,Sibuet J C,Wu J. 2019. Intermingled fates of the South China Sea and Philippine Sea plate[J]. Natl Sci Rev,6(5):886–890. doi: 10.1093/nsr/nwz107

图(12)  /  表(4)
计量
  • 文章访问数:  400
  • HTML全文浏览量:  136
  • PDF下载量:  157
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-04-10
  • 修回日期:  2023-04-30
  • 网络出版日期:  2023-05-29
  • 发布日期:  2023-05-14

目录

/

返回文章
返回