-
板块俯冲系统对地球板块运动起着至关重要的作用,它在关闭大洋的同时也带来了岩石圈内的岩浆活动并导致上覆板片的地壳增生和下插板片的回落,在威尔逊循环中起到了承前启后的作用,是地球构造演化不可或缺的一环。其俯冲过程重要的一个阶段—俯冲作用的起始过程对应着构造应力环境的转变,但由于这种地质现象不易观测到,俯冲起始过程因此成为板块构造理论中一个需要解决的问题[1]。针对起始过程的研究包括了俯冲起始机制的模拟和俯冲带形成的构造环境的控制因素两方面的内容。目前关于板块俯冲初始机制的大部分模式都来自数值和物理模拟的理论研究,缺乏直接的观测证据。其原因不仅在于俯冲带需要海陆联合观测,也因处于俯冲初始阶段的俯冲带数目非常有限[2]。Robert J Stern在回顾了全球针对俯冲带机制的研究后[3],根据俯冲带的起始机制将俯冲带分为两类[4-5],诱发俯冲型(Induced Subduction Initiation,简称ISI)和自发俯冲型(Spontaneous Subduction Initiation,简称SSI)。将可以自我维持进行的板块俯冲作为起始俯冲作用完成的标志的话[6],那么ISI型需要在长时间压应力作用下,使已经存在的薄弱结构在主导走向上逐渐收敛,使俯冲板片俯冲到一定深度直至其板块运动可以自我维持进行。而另一种SSI型,通常发育在大洋转换断层或者被动大陆边缘颈缩地区,一般是在走滑断裂导致或经岩石圈伸展的薄弱位置处发生。由于该地区两侧不同性质的板片的重力失稳进而塌陷,导致一侧板块的下插俯冲。这两种初始机制的分类是依据早期对由菲律宾海转换断层发育的伊豆-小笠原-马里亚纳(Izu-Bonin-Marianas)弧前火山弧和阿曼地幔蛇绿岩(Ophiolite of Oman)剥露的研究证据所建立的,但由于两个研究区俯冲作用至少已经分别持续了50和90 百万年[7-10],所以,仅依靠来自两个地区的岩石学证据和模拟实验结果来检验俯冲起始机制模式仍然比较困难。特别是由于针对俯冲初始过程所带来深部地幔脱水、深部岩浆运移、上覆地壳形变和增生范围等地质现象的观测还远远不足,因此导致对地幔深部作用机理研究较为薄弱。究其原因就是缺乏足够的、针对处于俯冲初始阶段的俯冲带深部地震观测和探测数据,影响了俯冲初始机制和地幔深部作用机理的研究。
印度尼西亚东部海域(苏拉威西海、北班达海)有着世界上最丰富的板块俯冲系统,其形成年代较晚,现今活动俯冲带主要发育在早中新世(约20 Ma)以后[11]。该地区俯冲带并不像世界上其他处于起始阶段的俯冲带—发育自洋中脊的转换断层地区(如南大西洋南三明治海沟和地中海直布罗陀弧),而是发育自被动大陆边缘—俯冲带不是由两个构造板块在压应力环境下发育,而是由伸展应力环境转变发育而成[11]。位于班达海的苏拉深渊(Sula Deep)和Tolo海槽,苏拉威西海的科托巴托海沟(Cotabato trench)和北苏拉威西海沟(North Sulawesi trench),它们的俯冲分别起始于洋中脊的转换断层或被动陆缘的洋陆过渡带(Ocean-Continent Transition)地区,属于典型的自发型俯冲。更为重要的是,以上4个俯冲带都是独立发育,并处于初始俯冲状态下的不同阶段,因此具有研究俯冲起始机制的独特优势。我们选择北苏拉威西俯冲带开展深部地震学观测,不仅可以获得其俯冲带岩石圈结构和浅层孕震区的地震活动性,还可以利用初始俯冲阶段的岩石圈结构开展动力学模拟,从而讨论自发型俯冲过程的构造控制因素。
-
印度尼西亚群岛被环太平洋俯冲带所围绕,是世界上火山活动最频繁的国家。通常认为,火山活动代表着板块已经俯冲到岩石圈300 km以下,同时也标志着俯冲系统进入了稳定的威尔逊循环。印尼东部地区位于欧亚大陆板块、印度-澳洲板块和菲律宾板块的汇聚边缘之间,不仅存在俯冲作用所导致的陆壳隆升、海沟后撤、增生杂岩等典型地质现象,而且在苏拉维西海还存在着俯冲转换传递边缘[12]、对向俯冲系统等独特的地质构造(图1)。
苏拉维西海在20世纪80—90年代经过海洋地球物理探测和海底钻探实验后(ODP leg 124)[14-19]基本确定了其边缘海的属性,它在海底扩张时期属于大洋的一部分,与西菲律宾海一样同为俯冲作用形成弧后盆地[20]。Robert Hall综合了前人的构造地质及地球物理证据,对东南亚区域开展了古板块重建研究[21]:苏拉威西海自晚始新世开始扩张(磁条带20,约47 Ma,位置见图1),持续了大约8 Ma,其洋壳与另一个正在扩张的大洋(西菲律宾海)通过一条转换断裂带相连。中始新世(~40 Ma)苏拉威西海扩张停止,而西菲律宾海持续扩张至早渐新世(~32 Ma),两个海洋板块之间转换断层因为不同的扩张速度和持续时间,断裂逐渐发育。晚渐新世(~24 Ma)随着澳洲板块持续向菲律宾板块俯冲,连接苏拉威西海与菲律宾海的转换断裂逐渐发育为新的俯冲带—摩鹿加海俯冲。摩鹿加板块俯冲到苏拉威西海东部之下,形成了桑义赫岛弧/海沟。在同一时期,位于苏拉威西岛南支(位置参考图1)与澳洲板块之间转换断层(Palu-Koro断裂的前身)也停止了活动,苏拉威西岛东南支随着俯冲作用从澳洲板块漂移,逐渐与南支拼接,同时苏拉威西岛的东支与北支开始缝合,形成了现今苏拉威西岛的原型。晚始新世也同时是南海扩张的活跃期,南海南部陆缘(现今南沙地块)随着南海的扩张,以南海西端为轴顺时针旋转,其陆缘与婆罗洲开始碰撞;这一时期巽他古陆南部(包括苏门答腊岛、巽他陆架、婆罗洲、爪哇岛等)也开始以古陆中心为轴逆时针旋转;巽他海沟范围也开始扩大,爪哇俯冲带向东延伸,小巽他群岛火山弧发育。中中新世(~15 Ma)随着南海关闭,南海的南部陆缘与西北婆罗洲(沙捞越地区)的碰撞也逐渐停止,这时由苏拉威西海向苏禄群岛的西北向俯冲,苏禄海开始了弧后海底扩张的过程。晚中新世(~12 Ma)随着爪哇俯冲带继续向东延伸,另一个弧后边缘海(北班达海)开始扩张。晚中新世至早上新世(~12~5 Ma)由于同时受到了南沙地块的板块碰撞和苏拉威西海的板块俯冲的双重作用,沙巴地区开始了造山运动,在此期间,苏拉威西-苏禄俯冲带边俯冲边后撤,爪哇俯冲带在后撤的同时也一直向东生长;早上新世(约5 Ma),苏拉威西-苏禄俯冲带停止了俯冲,苏禄海也停止其海底扩张的过程。难以理解的是在这一时期苏拉威西海另一侧陆缘开始了俯冲作用,俯冲方向调转为南向,苏拉威西海板块向苏拉威西岛北支下插入,而位于苏拉威西岛南支、北支和东支之间的古近纪断裂带重新活动。自上新世以来,随着苏拉威西海的南向俯冲作用,北苏拉威西海沟持续向北后撤(图2),导致了位于俯冲带上覆板块的岛北支和未受到俯冲作用的岛南支之间的古近纪断裂发展为Palu-Koro左旋走滑断裂,岛北支也同时以其东端为轴,开始顺时针方向的旋转(参见图2)且旋转角速度大约为4°/Ma,因此北苏拉威西俯冲带的板片俯冲速度是西快东慢。此外沿着海沟走向,长度大于500 km的增生楔前缘由于不同的俯冲速度,增生楔(或推覆岩席)的范围和厚度都不相同[22]。
经过5百万年的俯冲作用,虽然苏拉威西海古洋脊已经俯冲至苏拉威西岛北支之下,但是由于其俯冲深度有限[23],该俯冲并未在岛北支形成对应的火山弧(而位于岛北支东端的火山群是由于摩鹿加海俯冲作用形成的火山弧)。此外位于滑脱面的增生楔范围有限且保留了较为完整的变形序列[13],俯冲带从被动陆缘发展成主动陆缘这一过程仅仅持续了不到5百万年,且该俯冲带东西两侧板片俯冲程度也不同—由于岛北支的旋转,其俯冲速度自东向西线性递增。针对改俯冲带开展三维地震观测研究,可以在同一俯冲带下观测不同俯冲速度作用下的岩石圈结构,为研究自发俯冲型的初始机制提供第一手地震观测资料。
-
为了能够较为完整地观测北苏拉威西俯冲系统,地震台网的设计需要涵盖未受俯冲带影响的苏拉威西海洋盆、俯冲增生楔和弧前盆地地区。但是由于北苏拉威西俯冲带年龄仅有5百万年,尚未形成弧前盆地[13],为此,可将海底地震观测台网布设于增生楔后的陆坡上,同时还要考虑俯冲带转换边缘的影响,需要在岛北支西侧布设地震台以控制俯冲转换边缘(Palu-Koro断裂带)的地质结构。据此设计理念,中国科学院地质与地球物理研究所和英国剑桥大学地球科学系的科研人员共同合作,设计出北苏拉威西海俯冲带三维地震观测台网(图3)。2019年8月15—25日,中国科学院、英国剑桥大学、印尼万隆科技学院的科研人员使用印尼国家海洋地质研究所的科学调查船“海洋地质三号”和中国科学院地质与地球物理研究所研发的长周期宽频带海底地震仪(OBS),完成了OBS投放航次,共投放27台OBS,站位间隔为 50~70 km(图3),三维地震观测区域为苏拉威西海南部及望加锡海峡北部海区(俯冲转换边缘)。三维OBS台网预计观测时间为2019年9月至2020年6月,从而获得第一手的天然地震数据观测资料。
-
在苏拉威西岛北支及其周边海域开展关于北苏拉威西俯冲系统的三维地震观测研究的主要目的,是为了揭示该俯冲系统在有限的俯冲深度和角度作用下的上地幔活动过程,认识板块转换陆缘与海沟后撤之间的联系。在这个地区进行海底地震观测并开展天然地震学研究,可以利用陆地永久地震台数据和前人的地震探测数据(图4),开展不同尺度的地震学成像研究,获得北苏拉威西俯冲带的岩石圈结构、自东向西由于海沟后撤速度不同导致的增生楔/褶皱变形带的构造变化,并通过对比俯冲转换边缘两侧的岩石圈三维结构与属性,分析俯冲转换边缘在初始俯冲系统中扮演的角色,并讨论该俯冲带与古俯冲带形成的对向俯冲系统与该区域唯一的火山una-una活动之间的联系。以该区的地震学俯冲带深部结构研究为基础,利用有限元数值模拟分析技术,分析不同俯冲速度下岩石圈的形变程度,还原北苏拉威西海俯冲带的演化过程,为进一步认识自发俯冲型俯冲带的初始机制提供新的地球物理证据。
本项研究工作需要在国际合作基础上开展,主要以中国科学院地质与地球物理研究所和英国剑桥大学地球科学系、自然资源部第二海洋研究所等机构于2018年开始的针对研究区科学问题的国际合作为基础,同时得到了国家自然科学基金委员会重大研究计划“西太平洋地球系统多圈层相互作用”下属项目“苏拉威西海与古南海对向俯冲系统的三维地震观测与板块活动机制”的支持。
-
以揭示对向俯冲板片系统构造及地幔活动过程为目标,通过对海陆天然地震联合数据的处理,获得北苏拉威西俯冲带的岩石圈结构、自东向西由于海沟后撤速度不同所导致的增生楔/褶皱变形带的构造变化,并以此为基础开展有限元模拟,还原苏拉威西海俯冲演化过程,为研究自发型俯冲系统的初始机制提供重要的地球物理证据。因此,针对研究区的大规模的海底地震观测是实现以上目标的基础与关键,国家自然科学基金重大研究计划的设立与实施为实现这一目标提供了良好的契机与强有力的支撑。
致谢:感谢国家自然科学基金重大研究计划“西太平洋地球系统多圈层相互作用”项目的支持和国家重点研发计划资助项目国家质量基础的共性技术研究与应用“战略性新兴产业关键国际标准研究课题(2016YFF0202704)的资助。本项研究同时也得到了自然资源部第二海洋研究所李家彪院士的学术指导和中国科学院地质与地球物理研究所游庆瑜研究员及其OBS研发团队的大力支持,OBS投放航次是由中科院地球物理研究所郝天珧研究员团队与剑桥大学地球科学系Nick Rawlinson教授团队合作完成,航次得到了印尼合作方万隆科技学院和印尼海洋地质研究所的支持,特别是印尼科学调查船“海洋地质三号”的全体船员为投放航次的完成提供了有力保障,在此一并致谢!
Introduction to the 3-D seismic observation of the north Sulawesi subduction zone and the study of initiation mechanism
More Information-
摘要: 苏拉威西海具有独特的演化过程—位于北婆罗洲及苏禄岛弧下的苏拉威西海板块停止俯冲后,其俯冲板块的另一侧开始向苏拉威西岛俯冲,形成了诸如俯冲转换边缘、对向俯冲系统等独特的地质现象。针对这种同一板块一侧俯冲停止后在板块另一侧发育新俯冲的过程,科学界提出了众多的板块初始俯冲机制模式,如前陆盆地碰撞后发育新俯冲、大陆/岛弧边缘重力失稳后解耦后俯冲、岩石圈浅部小范围的地幔对流活动导致俯冲等14种初始模式。由于难以获得处于俯冲初始阶段俯冲带的深部岩石圈的地震观测资料,造成这些板块活动模式仍止步于推测。苏拉威西海区的天然地震观测资料不仅是研究苏拉威西海俯冲系统岩石圈深部构造及演化的基础,而且也是认识俯冲过程初始机制的关键一环。本文重点介绍了国家自然科学基金委员会重大研究计划“西太平洋地球系统多圈层相互作用”所属重点项目“苏拉威西海与古南海对向俯冲系统的三维地震观测与板块活动机制”在苏拉威西海区针对北苏拉威西海俯冲板片结构所开展的三维地震观测研究。该项研究实验以揭示俯冲系统深部结构与板块活动机制为目标,通过多边国际合作,于2019年8月15—25日在苏拉威西海区成功投放了27台国产宽频带海底地震仪(BBOBS),拟采集10个月的天然地震数据。同时计划利用印尼合作方的陆区地震台站数据,与OBS数据一同开展研究区岩石圈三维结构的地震学研究与数值模拟,认识俯冲下插板片的变形形态与特征、上覆地壳增厚程度及与俯冲系统持续作用过程之间的关系,讨论该地区独特俯冲系统的初始机制。Abstract: The Celebes sea has a unique evolutionary history. When the oceanic crust of Celebes below the north Borneo and Sulu arcs stopped subducting, the plate on the other side began to subduct towards the northern arm of the Sulawesi Island. Several unique geological structures, such as the subduction transform edge propagator, double inward subduction were formed. This dynamic tectonic process, from subduction cessation to initiation in a relatively short period of time is a distinctive feature of the western Pacific subduction system. To explain its physical mechanisms of the process, many models to explain the subduction initiation mechanism have been proposed, such as the post-collision development from foreland arcs, gravity instability initialization at the arc/continental margins and sub-lithospheric mantle convection etc.. Due to deficiency of high resolution seismological data of the incipient subduction zone structure, these mantle dynamic models are still remaining in debate. Acquiring earthquake observation data from the incipient subduction system is not only the basis for studying the deep structure and evolution of its lithosphere, but also a key element in understanding the initial mechanism of the subduction process. This paper introduces the 3-D offshore earthquake observation experiment focusing on the structure of the incipient subduction zone of the North Sulawesi, which is supported by the National Science Foundation of China key project " 3-D seismic observation and plate movement mechanism of the inward subduction system in Celebes Sea and proto-South China Sea”. The objective of this experiment is to reveal the deep structure of the incipient subduction system and the mechanism of plate movement, 27 broadband ocean bottom seismometers (OBS) were deployed from August 15 to 25, 2019 through international collaboration to record seismic data for 10 months. It plans to study the 3-D structure of the lithosphere of the north Sulaweis subduction zone and carries out numerical simulation analysis of its plate movement. Using the OBS and the land seismic stations data together, we can study the deformation pattern and characteristics of subducting plates, overlying crust and its relationship with the continuous shorting of double subduction, to reveal the mechanism of this incipient subduction system.
-
图 1 苏拉威西海及苏拉威西岛构造简图
主要俯冲带、磁条带及断裂带位置根据Kopp等. 1999[13])。古俯冲带为灰色虚线加三角,磁条带为蓝色虚线,俯冲带为黑色实线加三角,黑色实线为主要的断裂带
Figure 1. Tectonic sketch of Celebes Sea and Sulawesi Island
Major subduction zones, magnetic isochrons and fault zones referred to Kopp et al. 1999[13]). The paleo-subduction zone is shown by gray dashed line with triangle, the magnetic strip shown by dotted line, the subduction zone shown by black solid line with triangle, and the black solid line with arrows indicate the main fault zone
图 3 三维海底地震观测及OBS投放位置图
红色五角星为实际OBS投放点;白色三角为当地永久地震观测台站,灰色三角为剑桥大学临时地震观测台站,紫色三角为火山
Figure 3. Deployment locations of OBS and 3-D seismic survey
(The red dots are the originally designed location of OBS, and the orange stars are the actual deployment location. The white triangles are the local permanent seismic stations. The grey triangles are the temporary seismic stations of University of Cambridge. The purple triangles are volcanos)
-
[1] Mckenzie D P. The initiation of trenches: a finite amplitude instability[M]//Talwani M, Pitman W C. Island Arcs, Deep Sea Trenches, and Back-Arc Basins. Washington, D. C.: AGU, 1977, 105: 57-61. [2] Gerya T. Introduction to Numerical Geodynamic Modelling[M]. 2nd ed. Cambridge: Cambridge University Press, 2019: 369-405. [3] Stern R J. Subduction zones [J]. Review of Geophysics, 2002, 40(4): 1012. doi: 10.1029/2001RG000108 [4] Stern R J. Subduction initiation: spontaneous and induced [J]. Earth and Planetary Science Letters, 2004, 226(3-4): 275-292. doi: 10.1016/S0012-821X(04)00498-4 [5] Stern R J, Gerya T. Subduction initiation in nature and models: a review [J]. Tectonophysics, 2018, 746: 173-198. doi: 10.1016/j.tecto.2017.10.014 [6] Lévy F, Jaupart C. The initiation of subduction by crustal extension at a continental margin [J]. Geophysical Journal International, 2012, 188(3): 779-797. doi: 10.1111/j.1365-246X.2011.05303.x [7] Uyeda S, Ben-Avraham Z. Origin and development of the Philippine Sea [J]. Nature Physical Science, 1972, 240(104): 176-178. doi: 10.1038/physci240176a0 [8] Casey J F, Dewey J F. Initiation of subduction zones along transform and accreting plate boundaries, triple-junction evolution, and forearc spreading centres–implications for ophiolitic geology and obduction[M]//Gass I G, Lippard S J, Shelton A W. Ophiolites and Ocean Lithosphere. Geological Society, London, Special Publications, 1984, 13(1): 269-290. [9] Coleman R G. Tectonic setting for ophiolite obduction in Oman [J]. Journal of Geophysical Research: Solid Earth, 1981, 86(B4): 2497-2508. doi: 10.1029/JB086iB04p02497 [10] Boudier F, Ceuleneer G, Nicolas A. Shear zones, thrusts and related magmatism in the Oman ophiolite: initiation of thrusting on an oceanic ridge [J]. Tectonophysics, 1988, 151(1-4): 275-296. doi: 10.1016/0040-1951(88)90249-1 [11] Hall R. The subduction initiation stage of the wilson cycle [J]. Geological Society, London, Special Publications, 2018, 470(1): 415-437. [12] Govers R, Wortel M J R. Lithosphere tearing at STEP faults: response to edges of subduction zones [J]. Earth and Planetary Science Letters, 2005, 236(1-2): 505-523. doi: 10.1016/j.jpgl.2005.03.022 [13] Kopp C, Flueh E R, Neben S. Rupture and accretion of the celebes sea crustrelated to the North–Sulawesi subduction: combined interpretation of reflection and refraction seismic measurements [J]. Journal of Geodynamics, 1999, 27(3): 309-325. doi: 10.1016/S0264-3707(98)00004-0 [14] Karig D E. Origin and development of marginal basins in the western Pacific [J]. Journal of Geophysical Research, 1971, 76(11): 2542-2561. doi: 10.1029/JB076i011p02542 [15] Weissel J K. Evidence for Eocene oceanic crust in the Celebes basin[M]//Hayes D E. The Tectonic and Geologic Evolution of Southeast Asian Seas and Islands. Washington, D. C.: AGU, 1980: 37-47. [16] Silver E A. A new tectonic map of the Molucca Sea and East Sulawesi, Indonesia, with implications for hydrocarbon potential and metallogenesis[M]//Barber A J, Wiryosujono S. The Geology and Tectonics of Eastern Indonesia. Bandung: Geological Research and Development Centre, Special Publications, 1981, 2: 343-347. [17] Silver E A, McCaffrey R, Smith R B. Collision, Rotation, and the Initiation of Subduction in the Evolution of Sulawesi, Indonesia [J]. Journal of Geophysical Research: Solid Earth, 1983, 88(B11): 9407-9418. doi: 10.1029/JB088iB11p09407 [18] Hinz K, Fritsch J. Geoscientific Investigations in the Celebes Sea and on the adjacent Subduction Zones (GIGICS)-SO98, Report on the SONNE cruise SO98/2, Geophysical Survey inthe Indonesian Part of the Celebes Sea, 1 November 1994-28 November 1994[R]. Unpubl. report BGRno. 113.095, 1995. [19] Rangin C, Silver E A, Von Breymann M T, et al. Summary of shipboard results[C]//Proceedings of the Ocean Drilling Program, Initial Repo. College Station, TX, 1990, 124: 415-419. [20] Rangin C, Silver E A, Von Breymann M T, et al. General introduction[C]//Proceedings of the Ocean Drilling Program, Initial Repo. College Station, TX, 1990. [21] Hall R. Late Jurassic–Cenozoic Reconstructions of the Indonesian Region and the Indian Ocean [J]. Tectonophysics, 2012, 570-571: 1-41. doi: 10.1016/j.tecto.2012.04.021 [22] Otofuji Y I, Sasajima S, Nishimura S, et al. Paleomagnetic evidence for clockwise rotation of the northern arm of Sulawesi, Indonesia [J]. Earth and Planetary Science Letters, 1981, 54(2): 272-280. doi: 10.1016/0012-821X(81)90010-8 [23] Surmont J, Laj C, Kissel C, et al. New paleomagnetic constraints on the Cenozoic tectonic evolution of the North Arm of Sulawesi, Indonesia [J]. Earth and Planetary Science Letters, 1994, 121(3-4): 629-638. doi: 10.1016/0012-821X(94)90096-5 -