西太平洋俯冲板块正断层系统与动力学特征

王悠昆，周志远，林间,，张帆

1. 南方科技大学海洋科学与工程系， 深圳 518055

2. 中国科学院边缘海与大洋地质重点实验室，中国科学院南海海洋研究所，南海生态环境工程创新研究院，广州 511458

3. 中国-巴基斯坦地球科学研究中心，中国科学院-巴基斯坦高等教育委员会，巴基斯坦伊斯兰堡 45320

俯冲带是地球循环系统的重要组成部分。在俯冲过程中，俯冲板块沉入地幔深处，并因其板块弯曲，在板块上部产生拉伸应力，在板块下部产生挤压应力[1-4]。大量研究表明，俯冲带中由板块弯曲引起的正断层在上地幔蛇纹石化、板内地震、板块中的流体活动以及外缘隆起带断层引起的海啸中发挥着重要作用[5-6]。西太平洋处于4 个板块交界之处，北为北美板块，西为欧亚板块，西南为菲律宾海板块，东南为太平洋板块。其中，北美板块和欧亚板块为大陆板块，菲律宾海板块和太平洋板块为大洋板块，两个大洋板块向两个大陆板块俯冲汇聚。西太平洋俯冲带是板块年龄最老、海沟最深和板块挠曲程度最大的俯冲系统。其洋盆发育了众多海山和海底高原，对海沟的几何形态产生了较大影响[7]。板块边界地震、火山活动活跃，西太平洋俯冲板块中的正断层地震通常发生在板块边界附近，因此可能构成巨大的海啸威胁。例如，1933 年，日本海沟发生了8.4 级地震，这是有记录以来最大的俯冲板块外缘隆起带（Outer rise）正断层地震，在日本三陆沿海地区引发了海啸[8]。2009 年9 月29 日，汤加海沟发生8.1 级外缘隆起带正断层地震，并引发了毁灭性的海啸[9]。因此，研究西太平洋俯冲板块正断层的动力学机制具有重要意义。

已观察到正断层在外缘隆起区到海沟轴部之间普遍存在。观测显示，与板块弯曲有关的正断层可能是板块俯冲过程中产生的新断层，或是在大洋中脊形成的重新激活的深海丘陵断层[10]。前期地球动力学模拟研究表明，正断层开始在距离海沟轴部一定距离处形成，然后向海沟轴部方向生长[11-12]。海沟附近发育了丰富的挠曲正断层，这些断层除了会诱发板内地震之外，也是流体进入板块内部和俯冲带深部的主要通道[13]。西太平洋海沟的观测与模拟研究相对丰富，包括汤加、日本、伊豆-小笠原和马里亚纳海沟等。这些俯冲带的构造特征变化很大，包括海沟深度、俯冲倾角、板块内和板块间地震活动，这使它们成为研究海沟动力学和板块相互作用的理想场所。本文综述了西太平洋俯冲板块弯曲与正断层的观测，并总结分析了正断层模拟研究揭示的正断层形成过程，这对进一步揭示俯冲带动力学机制有着十分重要的意义。

1 西太平洋俯冲板块正断层特征

汤加、日本、伊豆-小笠原和马里亚纳海沟均位于西太平洋（图1），且都是板块年龄相对较老（均超过100 Ma）的俯冲板块，因而远端的板块厚度可能相对较大，并且在4 个海沟之间变化不大，因此可以直接比较各海沟的正断层特征，并揭示它们的共同特征。高分辨率海底多波束测深数据是从NOAA 美国环境信息中心的多波束测深数据库（MBBDB）和海洋地球科学数据系统（MGDS）的全球多分辨率地形合成（图1）而成，网格平均分辨率约为100 m[14]。基于此数据，前人分别提取了汤加、日本、伊豆-小笠原和马里亚纳海沟的14、9、15 和15 个垂直于海沟的剖面，以计算平均断层走向和密度，得到4 个海沟的正断层特征（图2）。高分辨率多波束测深数据显示，正断层在4 个海沟的俯冲板块上普遍存在。

图2 汤加海沟、日本海沟、伊豆-小笠原海沟、马里亚纳海沟海底地形图a-d 分别为汤加海沟、日本海沟、伊豆-小笠原海沟和马里亚纳海沟，白色线段为选取的断层剖面位置[15]。Fig.2 Seafloor bathymetry of the Tonga, Japan, Izu-Bonin, and Mariana Trenchesa-d: the Tonga Trench, Japan Trench, Izu-Bonin Trench, and Mariana Trench, respectively. White lines depict the deployment of transaction profiling[15].

通过对比实际观测与弹塑性变形模型，前人研究了沿着汤加、日本、伊豆-小笠原、马里亚纳海沟的板块挠曲与正断层特征（图2）。观测表明，平均海沟挠曲量在日本海沟最小（3 km），在马里亚纳海沟最大（4.9 km），而平均正断层垂直断距在日本海沟最小（113 m），汤加海沟最大（284 m）。而后模拟了俯冲板块在3 种构造加载的作用下发生弯曲变形并产生正断层的过程，3 种构造加载分别为：垂向加载（V0）、弯矩（M0）和水平拉张力（F0）。在板块挠曲与正断层特征的双重约束下，反演得到了4 个海沟的最佳模型解。

汤加海沟的平均断层垂直断距最大，最大值为420 m，平均值为284 m（图3a）。伊豆-小笠原海沟和马里亚纳海沟的断层垂直断距相似，最大值均为320 m，平均值分别为238 m 和148 m。在汤加海沟、日本海沟、伊豆-小笠原海沟和马里亚纳海沟，可识别的断层起始点距离海沟轴线分别为85、80、100 和115 km。

图3 汤加海沟、日本海沟、伊豆-小笠原海沟和马里亚纳海沟跨海沟剖面上每 5 km 的平均断层错距（a）和断层密度（b）[16]Fig.3 The average fault offset (a) and fault density (b) per 5 kilometers on the cross-trench profiles of the Tonga Trench, Japan Trench, Izu-Bonin Trench, and Mariana Trench[16]

据观测，马里亚纳海沟正断层密度最大，伊豆-小笠原海沟的正断层密度最小（图3b）。在马里亚纳海沟，从距海沟轴线 80 km 处开始，断层密度开始明显增加，而其他海沟的断层密度则在距海沟轴部近 50 km 处开始增加，表明在马里亚纳海沟，正断层带最宽，断层密度最大。

2 西太平洋俯冲板块变形机制

观测到的海底地形受到各种组成部分的影响，包括沉积物厚度、板块冷却引起的沉降和艾里均衡补偿地形[16-18]，去除这些影响后，可得到非均衡地形，能够最大程度地反映板块弯曲变形的情况。以非均衡地形作为板块变形程度的观测，结合薄板弯曲理论模型，前人反演了4 个海沟的最佳构造载荷。图4 为西太平洋4 个海沟俯冲板块的平滑弯曲形态。对于每个海沟，黑色细虚线表示单个剖面，红色粗曲线表示海沟的平均剖面。每个海沟截取的多条剖面上的最大变形量（W0）平均值用蓝点标记。远场参考海底深度用灰色线标记。红色箭头标记表示弯曲曲率降低到可忽略值（0.1×10-6m-1）的特征距离（Xc）。

图4 汤加海沟、日本海沟、伊豆-小笠原海沟和马里亚纳海沟南北部的观测挠曲量及其平均值黑色虚线细曲线显示单个剖面，红色粗曲线显示沟槽的平均剖面，远场参考海底深度用灰色线标记。蓝点所在位置为多个海沟轴部最大变形量W0 的平均值，红色箭头所指位置Xc 为弯曲曲率可忽略值[19]。Fig.4 Observed flexures and their average values of the Tonga Trench, Japan Trench, Izu-Bonin Trench, and Mariana TrenchThe thin black dashed curves: the individual profiles; the thick red curves: the average profile of the trenches; the grey line: the far-field reference seafloor depth. W0: the average of the maximum deformation in the axial part of several trenches. Xc: the negligible value of the bending curvature[19].

前人系统性地研究了西太平洋汤加-克马德克、马尼拉、菲律宾、日本以及马里亚纳海沟的俯冲板块变形研究并进行对比。图5a 为海沟附近板块弯曲示意模型，显示预期的拉伸屈服变形区（红色网格纹）和压缩屈服变形区（蓝色网格纹）。图5b中X0是板块的宽度，也是垂直变形W=0 的位置。在板块弯曲过程中，正断层发育的最大深度主要是由轴向垂直载荷（V0）控制，而正断层最大深度离海沟的距离是由轴向弯矩（M0）所控制的。西太平洋4 个海沟板块变形对比研究表明，最大变形量（W0）和板块宽度（X0）主要由靠近海沟的由于断层作用降低的有效弹性厚度控制，并且几乎不受远端俯冲板块的初始有效弹性厚度的影响。板块有效弹性厚度的降低导致了海沟的显著加深和变窄，而板块的有效弹性厚度的变化与其年龄又存在着密不可分的联系[19]。通过汤加-克马德克海沟与马尼拉、日本、马里亚纳海沟等俯冲板片弯曲的分析对比（图5b），发现对于较年轻或较老板块的情况，无论海沟处加载量如何变化，板块的年龄都可能是控制海沟弯曲形状的主要因素（图5）。

图5 海沟板块附近弯曲的示意图模型(a)和汤加-克马德克、马尼拉、菲律宾、日本以及马里亚纳海沟俯冲板块年龄与板块弯曲参数的关系(b)X0 为板块宽度，W0 为最大变形量，V0 和M0 是弯曲参数[19]。Fig.5 Schematic modelling of buckling near the trench plate (a) and the relationship between subducting plate age and plate buckling parameters in the Tonga-Kermadec, Manila, Philippine, Japan, and Mariana Trench subducting plates (b)X0 is the slab width, W0 is the maximum deformation, and V0 and M0 are the bending parameters [19].

3 西太平洋俯冲带板块正断层模拟与含水量估算

西太平洋的俯冲带板块正断层的模拟研究中，俯冲板块在板块弯曲演变过程中的水平偏差应力能直接反映断层形态。根据最佳拟合模型得到4 个海沟的俯冲板块的正断层模式（图6a-d）。根据计算，正断层发育在上塑性屈服破裂带内。

图6 汤加、日本、伊豆-小笠原、马里亚纳海沟俯冲板块正断层特征以及板块有效弹性厚度变化、计算的水平偏应力和有效弹性板厚度a-d 中黑色虚线表示伸展屈服带的最大深度，带误差条的黑圈显示了研究区域内可重新定位的外升正断层地震；e-h 中是4 个沟槽计算的Te（黑色曲线）和计算的面积SΔTe（白色区域）[15]。Fig.6 Normal faulting characteristics in the subducting plates of the Tonga, Japan, Izu-Bonin, and Mariana Trench and variations in the effective plate elastic thicknessThe black dashed lines in a-d indicate the maximum depth of the stretching yield zone, and the black circles with error bars show the relocatable outgoing uplift normal fault earthquakes in the study area; e-h are the calculated Te (black curves) and calculated area SΔTe (white areas) for the four trenches [15].

结果显示，在日本海沟、伊豆-小笠原海沟和马里亚纳海沟中，大多数正断层都是向海沟方向倾斜的。但汤加海沟既有向海沟倾斜的断层，也有向海洋倾斜的断层。经计算，日本海沟和汤加海沟的正断层比伊豆-小笠原海沟和马里亚纳海沟的正断层浅。研究区域内现有的重定位正断层地震均位于计算出的拉伸屈服破裂带内（图6a-d）。黑色虚线曲线表示拉伸屈服带的最大深度。带误差条的黑圈显示了Emry 和Wiens 研究区域内可重新定位的正断层地震[20-24]。

在构造加载的作用下，俯冲板块在距离海沟100 km 左右处的外缘隆起区开始产生正断层，逐步向海沟轴部发育，随着断层横向发育的过程中断层深度也逐渐增大，直至断层形态趋于稳定[16]。模型结果显示，日本海沟的水平张力分别比马里亚纳、汤加和伊豆-小笠原海沟小33%、50%和60%。汤加、日本、伊豆-小笠原、马里亚纳海沟的正断层最深可达海底以下29、23、32 和32 km（图6），这与重新定位后的日本与伊豆-小笠原地震深度一致。此外，反演得到的水平张拉力与观测到的平均垂直断距呈一定正相关性，而计算得到的有效弹性厚度减少量与观测到的海沟挠曲量也相关。这些结果表明，水平张拉力在正断层发展过程中起着关键控制作用，板块弱化可导致板块挠曲量的显著增加。

根据计算的水平偏应力，按照前人的方法[25-27]，计算了由于板块弯曲和正断层作用而产生的有效弹性板块厚度Te变化（图6e-h）。结果表明，有效板块厚度向海沟轴线逐渐减小，汤加、日本、伊豆-小笠原和马里亚纳海沟的最大Te减少量分别为25、24、22 和26 km。然后，通过对跨轴距离上的Te还原进行积分来计算Te还原的面积SΔTe。计算得出SΔTe的变化值中马里亚纳海沟最大，日本海沟最小。同时，计算得出Te还原区的宽度中日本海沟最大，其他3 个海沟几乎相同。因此，日本海沟的平均Te减少量最小，马里亚纳海沟的减少量最大。

大多数利用薄板弯曲理论的研究已经认识到垂直载荷和弯矩的重要性，然而很少有研究调查水平张力（F0）的重要性。基于先前的研究表明，F0对于解释板块弯曲和断层垂直断距至关重要，特别是在控制最大断层断距离海沟轴部距离方面[28-30]。计算的F0与观测到的4 个海沟的平均断层垂直断距呈正相关。这一结果表明，较大的F0有利于俯冲板块中较大正断层的发育。基于以上结果，可以推断F0在控制正断层模式中起着关键作用。由于研究存在局限性，即只研究了相对古老板块的俯冲带，因此想要证实这个结论还需要进一步的研究。

西太平洋俯冲板块因其年龄老、具有较高的刚度，形成的断裂分布广且断距大，进一步促进了流体进入地幔，并引起地幔蛇纹石化[31-32]。前人通过研究表明进入板块内地幔水化的范围和程度可以被用来估计带入俯冲带的水量[33-35]。汤加、日本、伊豆-小笠原和马里亚纳海沟的累计断层长度（即单个可识别断层的总和）分别为240、260、360 和450 km。并且通过对Cascadia 海沟的地震反射研究，估算单个断层周围的透水断层带宽度为75～600 km，从而估计透水断层带的体积为18.0～144.0、19.5～156.0、27.0～216.0 和34.0～270.0 km3，汤加、日本、伊豆-小笠原和马里亚纳海沟的地幔蛇纹石化百分比分别为0.4%～3.4%、0.4%～3.1%、0.6 %～5.1%和1.4%～10.8%。这一结果表明，马里亚纳海沟的地幔蛇纹石化程度可能明显大于其他3 条海沟，相当于汤加和日本海沟的350%，伊豆-小笠原海沟的230%。

4 结论与展望

（1）观测表明马里亚纳海沟、日本海沟、伊豆-小笠原海沟和汤加海沟都有显著的正断层特征，且汤加海沟的平均断层错距最大，马里亚纳海沟的正断层平均密度最大；（2）板块有效弹性厚度的降低导致了海沟的显著加深和变窄，并且无论海沟处的加载如何变化，板块的年龄都可能是控制海沟弯曲形状的主要因素；（3）屈服带模型揭示马里亚纳海沟的有效弹性厚度变化最多，导致其正断层特征更为明显，这也符合对正断层的观测。这些发现对于理解俯冲带的动力学过程具有重要意义。

当前的研究仍存在一定的局限性，包括：（1）地球动力学模型可能无法完全考虑实际地质过程中的所有复杂因素；（2）当前的俯冲板块弯曲动力学模型基本上都是二维模型，而实际的海沟走向并非直线而全部为曲线形态，因而亟需三维地球动力学模拟方法来解释观测到的板片弯曲和正断层形态；（3）很多海沟仍缺乏实测的高精度海底地形数据，限制了不同区域的对比研究。今后的研究应朝着以上方向去探索，以提高对西太平洋俯冲板块动力机制的更深理解。