南海东部马尼拉俯冲带深部结构新认识*

赵明辉, 程锦辉, 高金尉, 孙龙涛, 徐亚, 张佳政, 杜峰

1. 中国科学院边缘海与大洋地质重点实验室, 中国科学院南海海洋研究所, 广东 广州 511458;

2. 南方海洋科学与工程广东省实验室, 广东 广州 511458;

3. 中国科学院深海科学与工程研究所, 海南 三亚 572000;

4. 中国科学院油气资源研究重点实验室, 中国科学院地质与地球物理研究所, 北京 100029;

5. 湖北省地震局, 地震预警湖北省重点实验室, 湖北 武汉 430071;

6. 中国科学院大学, 北京 100049

南海作为西太平洋最大的边缘海之一, 是我国挺进深海大洋的根据地, 也是我国系统开展深海基础科学研究的首选(汪品先, 2009)。国家基金委“南海深海过程演变”重大研究计划力争通过“解剖一只麻雀”, 揭示边缘海的演变规律及其对海底资源和宏观环境的影响。在众多科学家大力推动下, 我国成功在南海实施的国际大洋发现计划 IODP349/367/368/368X 航次, 为探讨南海深部地壳结构、大陆岩石圈破裂机制和南海扩张演化生命史提供了更多的数据支撑(Larsen et al, 2018; Sun et al, 2019; Ding et al, 2020), 也为边缘海的演变历史树立起系统研究的典范(汪品先, 2012, 2020), 使南海成为国际海洋地球科学研究的天然实验室。

南海东部马尼拉俯冲带是新生代以来才逐步形成、正在活动的年轻俯冲带, 地震与火山活动极为活跃, 代表着南海从陆缘张裂、洋脊扩张、俯冲消亡的完整威尔逊旋回的顶点 (Ludwig et al, 1979;李家彪 等, 2004; 尚继宏, 2008); 作为南海形成以来惟一的俯冲汇聚边界, 对于完整解剖南海这只“麻雀”的演化过程, 马尼拉俯冲带是至关重要的不可回避的环节, 是解决南海的形成与构造演化科学问题的关键之一(吴时国 等, 2004; 赵明辉 等,2016)。2014 年, 国家基金委重大研究计划“南海深海演变过程”重点支持项目“南海东部马尼拉俯冲带深部结构探测与研究”获得资助, 执行期间(2014—2018 年)在马尼拉俯冲带开展了多次综合地球物理探测, 采集了宝贵的第一手海上数据, 取得创新性进展。本文介绍了该项目在海上取得的实验数据,以方便更多的科学家申请与使用; 并总结了项目执行期间取得的4 项重要研究成果, 为重大研究计划所提出的“从海底扩张到洋壳俯冲的南海生命史的骨架”科学问题提供重要支撑。

1 重要科学实验与关键数据

在国家基金委共享航次的资助下, 先后开展了5 次综合地球物理探测, 投放海底地震仪(Ocean Bottom Seismometer, OBS)台站73 台次, 海底电磁仪(Ocean Bottom ElectroMagnetometers, OBEM)5 台次, 累积放炮达13,872 炮, 成功获得了60 台OBS和5 台OBEM 台站的宝贵数据(图1)。包括: 2015年二维主动源地震测线OBS2015-2; 2015 年7—8 月部署了7 台长期宽频带OBS 天然地震观测; 2016 年4—6 月完成的三维OBS 地震探测; 2016 年7—8 月完成的二维主动源地震测线OBS2016-2; 及2015—2017 年间三次共5 个OBEM 台次的海上实验。

主动源二维地震测线 OBS2015-2(图 1a), 以20km 间隔共投放了9 台OBS, 回收9 台, 回收率100%, 剖面长度200km, 共激发有效炮1994 炮, 放炮期间同步进行了近水面单道电缆接收, 获得单道反射地震剖面200km(刘思青 等, 2017); 主动源二维地震测线OBS2016-2, 是沿着OBS2015-2 测线东西向延长(图1a), 共投放9 台OBS, 回收9 台, 回收率100%, 共放炮1078 炮, 长130km; 以及沿着OBS测线的多道反射地震调查, 圆满完成了横跨马尼拉俯冲带主剖面的深地震调查任务。在巴士海峡区域(图1a)实施三维主动源OBS 地震探测; 此次实验共投放48 台国产便携式OBS, 回收41 台, 放炮测线长2300km, 有效放炮10800 次, 其中5 条北东—南西向测线(Lx), 11 条北西—南东向测线(Ly)(Du et al,2018; 张佳政 等, 2018)。

图1 项目实施的重要科学实验与关键数据a 中圆圈与三角形分别代表主动源和被动源OBS 台站; b 中红色圆点为OBEM 台站, L1 测线对应于图4 中多道反射剖面的位置Fig. 1 Significant scientific experiments and key data acquired during this project. Circles and a triangle represent Ocean Bottom Seismometers (OBS) stations recording signals caused by active airgun-shooting and passive earthquakes, respectively,in (a). Circles in (b) indicate the Ocean Bottom ElectroMagnetometers (OBEM) stations; and the seismic line L1 corresponds to the multichannel reflection profile in Fig. 4

2015 年7—8 月, 通过搭载国家基金委共享航次, 在研究区部署了7 台长期宽频带OBS, 于2016年2 月仅回收1 台(POBS04)(图1a), 有效记录时段为6 个月, 记录了马尼拉海沟附近的天然地震数据,有效约束了马尼拉俯冲板片形态(任昱 等, 2020)。

2015 年7—8 月和2017 年5—6 月通过基金委共享航次, 进行了3 次OBEM 测试和采集工作。第1 次海试选择在南海北部进行海上采集实验, 主要检验OBEM 在深水环境工作情况, 验证OBEM 样机的工作状态及对深水环境探测能力, 为OBEM 投入实际应用积累实践经验。第2 次海试重点验证了OBEM 投放与回收操作性能。由于OBEM 电极臂较长, 投放与回收比较困难。经过在实验2 号的两次投放与回收实验对比, 采用尾部投放、侧位回收的方案较为理想。第3次海试, 共获得5 个OBEM 站位的有效数据(图1b)。其中OBEM-01 和OBEM-02 站点位于陆坡处, 海底地形起伏较大, OBEM-03 站点位于海盆边缘,OBEM-04、OBEM-05 站点位于中沙群岛附近, 进一步数据处理与分析正在进行之中(Xu et al, 2017)。

在项目成员与基金委南海北部地球物理共享航次通力合作和共同努力, 完成了各项地球物理数据的采集任务。为了更好地利用新采集的数据资料, 充分发挥共享航次的作用, 下船后主导科学家和有数据需求的科学家共同签署《数据共享协议》, 保证主导科学家项目的完成和数据优先使用权, 同时实现数据资料共享, 所有数据交给国家基金委资料共享中心, 共同推进南海深地震探测研究快速发展。

2 南海东北部地壳属性及洋陆边界的确定

2015 年7—8 月期间采集的广角反射/折射深地震测线OBS2015-2 位于马尼拉俯冲带前缘(图1a)。基于2015 年采集的多道地震剖面建立初始正演模型, 通过时深转换, 获取海底、沉积层、基底深度信息, 模型长189.6km。9 台OBS 记录的地震数据, 震相清晰。通过走时拟合与射线追踪, 发现来自深部的P 波震相信息, 包括PsP、Pg、PmP 和Pn, 并发现了下地壳高速层折射震相(PhP)与反射震相(Ph); 8个台站共拾取到7977 个走时(OBS02 台站无有效数据)(刘思青 等, 2017)。接着采用二维射线追踪程序RayInvr(Zelt et al, 1992)与走时层析成像方法Tomo2D(Korenaga et al, 2000)获取OBS2015-2 测线下方的速度结构(图2)。速度结构模型及OBS 台站下方一维速度结构揭示南海东北部地壳厚度为12～15km, 发育有下地壳底侵高速层, 前人认为这是东西向海底扩张形成的磁异常条带(C15—C17)(Hsu et al, 2004), 最新的研究认为是岩浆活动形成的具有高磁特征的火山岩脊; 在本质上, 南海东北部的地壳属于受到张裂期后岩浆活动影响的减薄陆壳(Liu et al, 2018; Zhao et al, 2018)。基于最新的IODP349 和IODP367/368 航次数据, 确定了南海北部陆缘的洋陆边界COB(图3)(Liu et al, 2018), 为南海停止扩张时洋壳范围的推断提供了证据。

图2 OBS2015-2 测线下方的速度结构a 、b 分别为正演和反演速度结构模型; c、d 是分别从a 和b 速度模型中提取的OBS 台站下方一维速度结构Fig. 2 Velocity structure along the profile OBS2015-2. (a) and (b) showing the forward and inverted velocity structures,respectively. (c) and (d) are one-dimensional velocity structures beneath each OBS station, which was extracted from the forward (a) and inverse velocity structures (b), respectively

图3 南海东北部陆缘洋陆边界图a 为南海东北部地形水深图; b 为自由空间重力异常图; 其中图中灰色线段为磁异常条带, 黑色交叉线条为张裂期后火山岩脊形成的构造高地(Briais et al, 1993; Hsu et al, 2004; Sibuet et al, 2016)Fig. 3 Continent-Ocean Boundaries in the northeastern continental margin of the South China Sea. (a) The bathymetry map in the northeastern SCS. (b) The free-air gravity map in the northeastern SCS. Grey lines indicate magnetic anomalies, and black lines with crosses in (b) indicate structural highs (elongated volcanic intrusions) referenced by Briais et al (1993), Hsu et al (2004), and Sibuet et al (2016)

3 马尼拉俯冲带增生楔前缘的精细结构

基于高分辨率多道地震反射剖面, 对马尼拉俯冲带增生楔前缘的构造变形、岩浆活动和深部反射等特征进行了分析和刻画(高金尉 等, 2018)。增生楔下陆坡部分由盲冲断层、构造楔和叠瓦逆冲断层构成, 逆冲断层归并于一条位于下中新统的滑脱面上, 滑脱面向海方向的展布明显受到增生楔之下埋藏海山和基底隆起的影响; 靠近海沟的海盆区域以及增生楔前缘岩浆活动开始于晚中新世末期并持续至第四纪。重震联合模拟结果表明, 俯冲带北段19°N—21°N 区域是南海东北部陆缘向海盆的延伸, 是陆壳高度减薄的区域, 并受到岩浆活动的强烈改造, 从而导致其浮力较大, 产生的平缓俯冲作用造成海沟北段的几何形态明显地向东凹进(图4)。

4 南海在停止扩张时(约15Ma)洋壳范围

南海在约15Ma 停止扩张时具有多大洋壳范围？马尼拉海沟何时开始俯冲？将沿着马尼拉海沟已经俯冲下去的南海板片恢复出来是解决南海形成演化动力学的关键。基于天然地震层析成像方法生成的P 波速度波动(dVP)图, 展示了南海俯冲板片已经向东俯冲到达海沟以下450km 深度(Wu et al,2016)。将南海东边已经俯冲到马尼拉海沟之下的洋壳范围展平到地球表面上, 俯冲板片沿海沟向东恢复到400～500km, 板片南北两侧dVP负值区域对应着陆壳, dVP正值区域通常对应着洋壳, 较小的负值区(-0.1%dVP)对应着减薄陆壳; 根据dVP正负差异, COB 在恢复的俯冲板片上向北延伸至400km,位于西琉球俯冲带以下, 在南北向上欧亚板块与碰撞前的吕宋岛弧平行(图5), 锆石裂变径迹研究(Lee et al, 2015)显示, 在5—6Ma, 吕宋岛弧在约400km长的欧亚板块边界上同时碰撞。上述研究结果很好地为南海地球动力学机制提供了新的约束。

5 建立南海与菲律宾海板块的构造演化模式

南海与菲律宾海板块在构造演化上关系密切。南海陆缘张裂主要分为两期; 第一次裂谷发生在始新世早期(56—40Ma), 以多裂谷盆地呈30°～45°排列为特征(Zhou et al, 1995; Lei et al, 2019); 第二次裂谷发生在始新世晚期(40—33Ma), 以新的45°～70°张裂方向, 形成的裂谷盆地中充满了黑色页岩, 夹砂岩层(Zhou et al, 1995; Lei et al, 2019)。第一次裂谷阶段仅发生在陆缘上部, IODP 钻孔U1501 的北部, 其南部没有发生(Larsen et al, 2018; Lei et al,2019)。而在菲律宾海域, 从磁异常分布可以分为两个主要扩张阶段(Doo et al, 2015): 扩张第一阶段始于Chron 24(53.2Ma), 可能更早于Chron 25(56.9Ma)或 26(58.6Ma); 由于Oki Daito 洋脊的后期形成过程中, 破坏了原来的构造, 导致这两个时间得不到很好的确定。第一阶段扩张结束发生在 Chron 20(42.5Ma)(Wu et al, 2016); 第二阶段开启始于Chron 19(40.4Ma), 止于Chron 13(33.7Ma)。

图4 马尼拉俯冲带北段增生楔多道地震反射剖面和解释[据高金尉等 (2018)修改]剖面位置见图1bFig. 4 Multi-channel seismic reflection profile (upper) and interpretation (lower) of accretionary prism of northern Manila subduction (after Gao et al, 2018). The location of this profile is shown in Fig. 1b

图5 南海在停止扩张时(15Ma)洋壳范围据Zhao 等(2019)修改。将俯冲下去的南海板片(速度波动图)展平到地球表面上, 沿海沟向东伸展了400～500km。绿色实线圈闭的是现今南海洋壳范围(Sibuet et al, 2016), 东部为已俯冲下去的洋壳范围。灰色阴影区域(宽度＜100km)为具有边界效应的不确实区域。粉色、黄色和浅蓝色分别代表南海不同扩张时期产生的洋壳, 扩张方向分别为N055°、N075°、N085°。洋壳内部的红色双虚线代表南海残余扩张脊(He et al, 2016)。红色粗线为COB 边界, 实心圆为IODP 钻探井位Fig. 5 The whole oceanic crust domain when the SCS seafloor spreading stopped (15Ma); after Zhao et al (2019). The subducted SCS plate (dVP) is unfolded and restored on the surface of the Earth, and the plate extends 400-500 km east along the Manila trench. The green line is the current SCS domain (Sibuet et al, 2016), and the east part is the domain of subducted SCS oceanic crust. The grey shaded area (＜100 km width) is probable artifacts. The SCS oceanic domain, characterized by N055°, N075° and N085° seafloor spreading directions, is delimited by pink, yellow and light blue areas, respectively. The red double dashed line indicates the extinct spreading ridge of the SCS (He et al, 2016), the red line indicates COB, and the circles indicate the location of IODP drilling sites

南海陆缘张裂的第一阶段(56—40Ma)刚好对应菲律宾海扩张的第一阶段, 从 53.2Ma(可能在Chrons 25 或26 之前几Ma)到42.5Ma(可能在Chrons 20 和19 之间)。南海陆缘张裂的第二阶段(40—33Ma)刚好对应于菲律宾海扩张的第二阶段(42.5—33.7Ma)。因此, 在始新世, 南海大陆边缘和菲律宾海板块(Philippine sea plate, PSP)均处于张性构造阶段(Zhao et al, 2019)。由于PSP 最初位于南海南侧约2200km 处, 有一个由南向北的运动形态, 在南海大陆边缘有伸展过程和海底扩张时, 这两个区域由左旋剪切板块边界连接。在15Ma(图6a), 将加瓜海脊和早白垩纪花东海盆外推到琉球海沟的当前位置(Hall, 2002)。截面AA′(图6c)显示了从南海东部盆地到加瓜海脊西部的PSP。随着PSP 不断向北移动,马尼拉俯冲带开始呈现左旋剪切运动, 大部分剪切运动持续到今天(Zhao et al, 2019)。吕宋弧在约12Ma 时开始形成, 大约6Ma 吕宋弧的北部开始沿南北走向、长400km 的部分, 与欧亚板块发生碰撞(图6b), 对应的截面BB’(图6c)展示了这一运动过程。

6 结论与展望

图6 15Ma 以来南海与菲律宾海板块演化过程示意图据Zhao 等(2019 )修改。a. 约15Ma 时, 南海洋壳沿着左侧剪切边界开始俯冲, 至菲律宾海板块下面; b. 菲律宾断裂带、花东海盆的增生部分和吕宋岛弧西部的当前位置; c. 截面AA′和BB′简化示意图, 截面位置在a 和b 中黑线。深绿色: 早白垩世花东海盆; 浅绿色:现今的菲律宾海板块; 中绿色: 浅灰色的菲律宾海板块部分在花东海盆的延伸, 随后俯冲消亡; 黄色: 海底扩张结束时南海的东部延伸; 橙色: Wu 等(2016)定义的东亚板块的一部分; 浅粉色: 苏禄海东北延伸Fig. 6 The kinematic evolution of the SCS and PSP since 15 Ma (after Zhao et al, 2019). (a) About 15 Ma, the oceanic crust of the SCS started to subduct along the left shear plate boundary to the bottom of the PSP. (b)The current location of the emplacement of the Philippine belt, the accretion of the western part of the Huatung Basin and the western Luzon arc. (c) A simplified diagram of profile AA’ and BB’ located in (a) and (b), respectively. Dark green: early Cretaceous HB; light green:current PSP; extension of the HB within the portion of PSP in light grey, which will subduct to be disappeared; yellow area:eastern extension of the SCS at the end of seafloor spreading; orange domain: part of East Asian Plate defined by Wu et al(2016); and pink domain: the northeastern extension of the Sulu Sea

为期8 年的南海深部计划研究积累了丰富的探索深海的基础资料, 取得了一系列创新研究成果(汪品先, 2020)。其中 “南海东部马尼拉海沟俯冲带深部结构探测与研究”是重点支持项目之一, 基于人工地震探测及天然地震层析成像结果, 确定了南海东北部地壳属于受到张裂后期岩浆活动影响的减薄陆壳, 同时联合国际大洋钻探成果, 划分了南海北部陆缘洋陆边界COB 位置, 刻画了马尼拉俯冲带北段增生楔前缘的精细结构, 进而圈定了南海停止扩张时的洋壳范围, 并初步构建了南海与菲律宾海板块的构造演化模型。该项目对于揭示南海形成演化及其动力学过程具有重要意义, 而南海也逐步成为深海科学研究的天然实验室及边缘海系统研究的典范。

随着基础数据的不断充实, 以及新技术新手段的不断引入, 南海研究也不断涌现新的科学问题,例如: 南海扩张后期岩浆活动的特征与机制; 马尼拉俯冲带地貌、构造和南海俯冲板块形态由北向南的变化及其空间耦合关系; 南海张裂-破裂机制与西太平洋边缘海系统的成因联系; 南海由东向西打开的动力来源; 此外, 位于欧亚板块与菲律宾海板块交汇处的花东海盆, 是现今西太平洋俯冲体系中唯一遗留的早白垩纪洋盆, 也成为研究西太平洋整体构造演化与动力学机制的关键区域(李春峰 等,2007; 黄奇瑜, 2017; 林间 等, 2017)。国家基金委“西太平洋地球系统多圈层相互作用”重大研究计划(2018—2025)立足科学前沿和国家需求, 其重点支持项目“花东海盆及其东/西部板块边界的跨圈层综合地震探测(2020—2023)”旨在加强花东海盆多学科的构造属性及形成演化研究, 为构建西太平洋俯冲构造体系跨圈层的构造演化模式提供重要的地质和地球物理证据(赵明辉 等, 2020), 推动板块构造理论发展与完善; 同时, 该项目也将为中国未来在花东海盆的大洋钻探计划(黄奇瑜, 2017)提供重要的基础数据, 引领科学前沿。