研究论文

马达加斯加岛西南陆缘CMD01剖面重震联合反演与地壳结构分析

  • 张永祺 , 1, 2 ,
  • 唐勇 , 1, 2, 3, 4, * ,
  • 董崇志 1, 2 ,
  • 吴招才 1, 2 ,
  • 李赫 1, 2 ,
  • 郭楚枫 1, 2, 4 ,
  • 任建业 3 ,
  • 巢鹏 3 ,
  • 杨天一 1, 2, 3
展开
  • 1.自然资源部第二海洋研究所,浙江 杭州 310012
  • 2.自然资源部海底科学重点实验室,浙江 杭州 310012
  • 3.中国地质大学(武汉) 海洋学院,湖北 武汉 430074
  • 4.上海交通大学 海洋学院,上海 200030
*唐勇(1974—),男,研究员,主要从事海洋地球物理方面的研究,E-mail:

张永祺(1999—),男,福建省福州市人,主要从事海洋地质方面的研究,E-mail:

收稿日期: 2024-02-05

  修回日期: 2024-05-17

  网络出版日期: 2025-02-08

基金资助

国家重点研发计划(2023YFE0126100)

国家重点研发计划(2023YFC2808803)

浙江省万人计划科技创新领军人才计划(22618)

浙江省“高层次人才特殊支持计划”科技创新领军人才项目(2021R52058)

Joint inversion of gravity and seismic along the CMD01 profile on the southwestern margin of Madagascar and analysis of crustal structure

  • ZHANG Yongqi , 1, 2 ,
  • TANG Yong , 1, 2, 3, 4, * ,
  • DONG Chongzhi 1, 2 ,
  • WU Zhaocai 1, 2 ,
  • LI He 1, 2 ,
  • GUO Chufeng 1, 2, 4 ,
  • REN Jianye 3 ,
  • CHAO Peng 3 ,
  • YANG Tianyi 1, 2, 3
Expand
  • 1. Key Laboratory of Submarine Geosciences, MNR, Hangzhou 310012, China
  • 2. Second Institute of Oceanography, MNR, Hangzhou 310012,China
  • 3. College of Marine Science and Technology, China University of Geosciences(Wuhan), Wuhan 430074, China
  • 4. School of Oceanography, Shanghai Jiao Tong University, Shanghai 200030, China

Received date: 2024-02-05

  Revised date: 2024-05-17

  Online published: 2025-02-08

摘要

本文以2021年中马国际合作航次获得的实测重力、磁力和二维多道地震剖面资料为基础,对横穿马达加斯加岛西南陆缘的CMD01剖面进行解释分析和重震联合反演,研究马达加斯加岛西南陆缘的地壳结构、构造单元及其演化过程。通过建立马达加斯加岛西南陆缘地壳密度模型,发现该陆缘具有减薄的陆壳,部分区段存在岩浆底侵作用;由陆向洋可划分出减薄陆壳、洋陆过渡带和洋壳等构造单元。其中,由于戴维断裂带的走滑运动对陆缘的影响,洋陆过渡带呈现出典型的转换型陆缘的构造特征。结合地壳密度模型和板块重建研究,重建了马达加斯加岛西南陆缘的发育演化过程。研究成果对于深化冈瓦纳大陆的裂解分离、东非被动陆缘的发育演化和海域划界等均具有重要的理论意义和实际应用价值。

本文引用格式

张永祺 , 唐勇 , 董崇志 , 吴招才 , 李赫 , 郭楚枫 , 任建业 , 巢鹏 , 杨天一 . 马达加斯加岛西南陆缘CMD01剖面重震联合反演与地壳结构分析[J]. 海洋学研究, 2024 , 42(4) : 83 -99 . DOI: 10.3969/j.issn.1001-909X.2024.04.008

Abstract

This study was based on the actual gravity, magnetic and multichannel seismic data obtained during the international cooperative survey in 2021. Through joint inversion of gravity and seismic along the CMD01 profile traversing the southwestern margin of Madagascar, combined with interpretation and analysis of seismic profiles, the crustal structure, tectonic units, and evolution of the southwestern margin of Madagascar were studied. A crustal density model of the southwestern margin of Madagascar was established, revealing thinned continental crust and the presence of magmatic underplating in certain segments. Structural units such as thinned continental crust, continental-ocean transition zone and oceanic crust domain can be delineated from land to sea, among which, due to the influence of Davie Fracture Zone strike-slip tectonics on the margin, the continental-ocean transition zone exhibited typical characteristics of a transform margin. By combining plate reconstruction and crustal density model, the developmental evolution of the southwestern margin of Madagascar was reconstructed. The research results have important theoretical significance and practical application value for deepening the understanding of the breakup and separation of the Gondwana, the developmental evolution of the East African passive margin, and maritime delineation.

0 引言

马达加斯加岛位于非洲大陆东南的印度洋上,与非洲大陆隔莫桑比克海峡相望,是印度洋上的最大岛屿,面积约为62.7 万km2,原为冈瓦纳大陆的组成部分,随冈瓦纳大陆的裂解分离而漂移到现今的位置[1-4]。整体上,马达加斯加岛为向西倾斜而多山的微陆块。
前人曾对马达加斯加岛西南陆缘及其相邻的区域进行过研究,认为长期的裂解过程导致了马达加斯加岛周边陆缘结构和类型的复杂性。COFFIN等[5]根据野外露头观测、深海钻探计划和地震数据解释重建了冈瓦纳大陆早期分裂时马达加斯加岛和非洲板块之间的地质演化过程,确定了马达加斯加岛西部盆地的共轭陆缘为坦桑尼亚-莫桑比克北部陆缘。SINHA等[6]利用坦桑尼亚大陆边缘的二维地震剖面进行重力-地震反演(即重震联合反演),建立了地壳密度模型,认为坦桑尼亚陆缘洋壳与陆壳之间无明显的过渡性地壳,显示为转换陆缘结构。VORMANN等[7]基于莫桑比克北部陆缘的广角折射/反射地震资料,对莫桑比克北部陆缘和戴维断裂带(Davie Fracture Zone,DFZ)地壳结构进行了分析研究,通过重震联合反演,由陆到洋未发现向海倾斜反射层和由下地壳高速体构成的洋陆过渡带(ocean-continent transition,OCT),而是以DFZ截然过渡。虽然前人的研究认为马达加斯加岛西部陆缘与坦桑尼亚陆缘互为共轭陆缘,且DFZ为陆缘的洋陆边界,并在共轭陆缘中发现了转换型陆缘的特征,但是目前针对马达加斯加陆缘处研究较少,缺乏对马达加斯加岛西南陆缘结构特征和地质属性的明确界定,无法对以上推论进行验证,这与该区域重力、磁力异常数据、高品质的多道地震数据等关键资料稀少有关。
2021年,中国-马达加斯加国际合作航次实施,获得了多条实测重力、磁力和二维多道地震资料。本文选择其中横穿马达加斯加岛西南陆缘的CMD01剖面开展详细的地震剖面解释,通过地形资料和自由空间重力异常(free-air anomaly,FAA)(以下简称重力异常)资料分析以及重震联合反演,构建陆缘地壳密度模型,划分陆缘构造单元,阐明陆缘构造演化过程。现今,划定200海里以外大陆架的外部界限(简称外大陆架划界)事关国家海洋主权权益、资源利用和国家安全,而马达加斯加岛陆缘的地质属性和地壳结构是马达加斯加划界案能否通过的关键点。本研究可深化对冈瓦纳大陆裂解分离过程的认识,明确马达加斯加岛陆缘的地质属性,为马达加斯加大陆架划界提供数据支撑,加强中国-马达加斯加两国在海洋科技领域的合作与交流。

1 区域地质背景

本文研究区为马达加斯加岛西南陆缘,研究区向洋一侧为南北向扩张形成的莫桑比克海盆,该海盆在早侏罗世时期开始扩张,与索马里海盆皆为冈瓦纳大陆裂解以及南北向扩张过程中形成的最古老的深海海盆;研究区南部为脊状地形的马达加斯加海台,玄武岩和拉斑玄武岩证据表明其形成与马里昂热点有关,但其地壳属性目前仍存在争议[8-9];北部与戴维断裂带相接,该断裂带为东非海岸洋壳与陆壳的转换边界,在13°S—23°S 构造特征较为明显,表现为近南北走向,是连接东非海岸与马达加斯加的构造带;东侧发育相邻的穆伦达瓦盆地与马任加盆地,两者均是中生代以来随着冈瓦纳大陆裂解、印度洋形成而产生的沉积盆地,具有明显的双层结构,下部裂谷层系包括上石炭统、二叠系、三叠系及侏罗系,以陆相沉积充填为主,上部坳陷层系为白垩系及新生界,全部为海相沉积充填(图1)。
图1 马达加斯加岛西南部陆缘周边区域构造图

Fig.1 The tectonic map of the region around the southwestern margin of Madagascar Island

马达加斯加岛地质背景复杂,受冈瓦纳大陆裂离影响,主要经历了三个阶段:陆内裂谷阶段(晚石炭世至早侏罗世)、马达加斯加漂移阶段(中侏罗世至早白垩世)、被动陆缘阶段(晚白垩世至今)。
晚石炭世至早侏罗世即陆内裂谷时期,卡鲁地幔柱剧烈活动引起裂谷作用,在冈瓦纳大陆东部产生大量可容空间,充填了卡鲁群地层[5,10-12]
中侏罗世至早白垩世,东、西冈瓦纳大陆开始分离[5],马达加斯加岛西部及西北部陆缘开始发育。现今在莫桑比克海盆内识别出的磁力异常条带以及坦桑尼亚边缘附近发现的一系列拉张盆地,证明初始洋壳形成的主要构造作用为索马里海盆与南部莫桑比克海盆的斜向扩张。直至距今153 Ma,东、西冈瓦纳之间的伸展运动逐渐被转换运动取代[2],东冈瓦纳大陆(南极洲、印度、澳大利亚、马达加斯加)沿戴维断裂带向东南方向斜向漂移。
晚白垩世,索马里海盆扩张停止,马达加斯加岛不再漂移,其西部陆缘基本定型,东非大陆边缘进入被动陆缘时期[4]。晚白垩世末,马达加斯加岛受马里昂(Marion)地幔柱的影响,与印度板块之间发生海底扩张,形成马斯克林海盆[13]。古新世早期,塞舌尔陆块与印度陆块分离[14],随后印度陆块向东北方向移动[15]
马达加斯加岛陆缘的演化过程较为复杂,包括裂谷初期的斜向伸展、索马里海盆和莫桑比克海盆引起的渐进式海底扩张以及沿大型转换断层(戴维断裂带)由北至南的漂移[16]。多种运动的叠加在马达加斯加岛形成了不同成因、类型和特征的大陆边缘[17-18],主要包括发育较厚卡鲁群地层的西部陆缘、长而笔直的东部陆缘以及与海台相接的南部陆缘。马达加斯加岛西部陆缘及其共轭陆缘坦桑尼亚-莫桑比克北部[2,19-22],是马达加斯加岛演化过程中重要的应力转变位置。

2 数据及方法

本研究采用地震剖面构造地层解释和重震联合反演的方法研究马达加斯加岛西南部陆缘地壳结构。基于地震反射、折射以及沿测线的实测重力异常数据,应用重震联合反演软件进行人机交互物性反演,以地震反射、折射剖面作为约束条件,确定关键界面的埋藏深度。通过拟合迭代反演物性结构,使得模型计算的重力异常与实测重力异常数据之间的误差最小,最终得到研究区的地壳密度模型。

2.1 数据来源

本研究使用的多道地震剖面(CMD01)位于马达加斯加岛西南部陆缘(位置见图1),2021年由大洋号采集,剖面长度约为502 km。多道地震电缆长度为6 km,共480道,气枪总体积为0.1 m3,炮间距为37.5 m,采样率为2 ms,记录长度为14 s。通过叠前去噪、地表一致性振幅补偿、直达波压制、多次波组合压制、偏移速度分析和叠前时间偏移等处理方式,获得叠前时间偏移剖面。使用ZL11-A型海洋重力仪同步采集重力数据,使用G880磁力仪同步采集磁力数据。盆地边界数据来源于Energy Portal by S&P Global(https://energyportal.ci.spglobal.com),海底地形和海洋卫星测高重力异常数据来源于全球卫星重力异常数据库(https://topex.ucsd.edu/marine_grav/mar_grav.html),数据的网格间距约为1'×1',局部地区重力异常数据的精度达到1.0 mGal。

2.2 地层界面识别与地壳密度模型建立

2.2.1 关键地质界面识别

在叠前时间偏移剖面的基础上,根据剖面的地震反射特征,参照前人在东非陆缘和马达加斯加岛西部陆缘盆地划分的地层格架[9,16,23-25],开展了构造地层的解释和分析。在CMD01地震剖面中识别出6个主要地层界面与5个地震反射层组,分别标记为T0、T1、T2、T3、T4和Tb以及U1、U2、U3、U4和U5,同时还识别出了地壳底界面——莫霍面(表1,图2,图3)。
表1 马达加斯加岛西南陆缘层序地震反射特征

Tab.1 Seismic reflection characteristics of sequence in the southwestern margin of Madagascar Island

层序 内部反射特征 频率 振幅 连续性 地震相
U5(T0~T1) 平行-亚平行 高频,变化较小 低-中振幅,局部空白 好-中等
U4(T1~Tb) 亚平行,局部嘈杂状 中频 低-中振幅,局部空白 中等-差
U3(T1~T2) 平行-亚平行,局部波状 高-中频,局部有变化 中振幅,局部空白 好-中等
U2(T2~T3) 亚平行,局部嘈杂状 中频 中振幅
U1(T3~T4) 亚平行,局部嘈杂状 变化较大 高-中振幅 好-中等
图2 马达加斯加岛西南陆缘地层格架图

(T4形成于晚石炭世陆内裂谷开始初期,未在图中标注。图件根据文献[16]修改。)

Fig.2 Stratigraphic framework of the southwestern margin of Madagascar Island (T4 was formed at the beginning of the Late Carboniferous intracontinental rift period and is not noted in the figure. Figure is modified from reference [16].)

图3 CMD01二维地震剖面

(U表示地层,其中U1为T3与T4之间的地层,U2为T2与T3之间的地层,U3为T1和T2之间的地层,U4为T1和Tb之间的地层,U5为T0和T1之间的地层。)

Fig.3 CMD01 2D seismic profile

(U represents stratum, U1 is the stratum between T3 and T4, U2 is the stratum between T2 and T3, U3 is the stratum between T1 and T2, U4 is the stratum between T1 and Tb, U5 is the stratum between T0 and T1.)

图3所示,在CMD01剖面左端起始一直到268 km,可见清晰的声波基底Tb处于双程走时7.1~8.7 s深度。在剖面0~50 km,双程走时10.2 s左右的位置,可见比较强的断续反射,结合前人研究发现的这个区域比较显著的线性磁力异常特征[18-19],我们解释该界面为洋壳的莫霍面。Tb至莫霍面之间的双程走时间隔大致为2 s,与6~7 km的标准洋壳厚度的双程走时2 s相同。从50 km向右,该剖面上的断续反射不再明显,主要根据“10 s规则”来推断莫霍面的位置。WARNER[26]研究了全球洋壳的莫霍面,提出在地震反射走时剖面上,莫霍面的表观符合所谓的“10 s规则”,即地壳处于均衡状态时,莫霍面会相当稳定地出现在双程走时为10 s的深度,并表现出相当的平坦度。莫桑比克海盆属于中生代的洋壳,已达到均衡状态,因此基于上述识别标准和规则,推断CMD01剖面50~268 km区段的莫霍面在双程走时为10~10.2 s的位置。在这个区段,尽管缺乏连续的强反射,但是断续的较强反射仍然模糊可见。在0~268 km区段内,除了局部的洋壳厚度小于2 s双程走时(如50 km处为1.6 s)外,其余洋壳的厚度均在标准洋壳的厚度范围内。从268 km向右,代表莫霍面的断续反射逐渐加深到双程走时10.5 s左右,推断可能与这个区段后来的地幔柱活动有关。
如上所述,Tb为莫桑比克海盆洋壳顶界面,而在莫桑比克海盆中识别出的磁力异常条带M0r,测线洋壳年龄为距今120 Ma左右[18-19]。总体上看,Tb界面在0~268 km区段基本水平延伸,上方存在明显不整合面。在268~500 km区段内,双程走时由7.7 s逐渐向上抬升至2.2 s附近,界面的反射强度不大,中等连续,中间有断层切割而错移,该界面标注为T3。T3界面之下为楔状地层U1,表现为连续性较好的亚平行反射特征,振幅高-中,局部为嘈杂状。T3受同生断层控制,其上沉积层的厚度变化不大,并且还具有对U1的削截反射,由此认为T3为同裂陷地层顶界面和裂后沉积地层底界面,属于马达加斯加西南被动陆缘发育的裂后不整合界面,代表晚侏罗世陆缘张裂活动结束、索马里海盆与莫桑比克海盆扩张开始、洋壳形成以及马达加斯加岛进入漂移阶段。该界面在东非大陆边缘广泛分布,是冈瓦纳大陆破裂形成的不整合面[24-25]
在剖面306 km处,断层右侧出现T2界面,位于T3界面上方,界面反射强度不大,连续性一般,被断层切割错移。由306 km向右(在306~500 km区段内),T2界面双程走时由6.5 s逐渐向上抬升至2.0 s附近。T2与T3界面之间的地层U2(由陆向洋方向),厚度基本没有变化,内部表现为亚平行状、较好连续性的反射,振幅中等,在306 km处受断层影响尖灭。T2之上为披覆地层。由早白垩世时期马达加斯加岛停止漂移的时间与该区域洋壳形成时间相近[20,27]可知,剖面中马达加斯加陆缘漂移期沉积应只分布于陆壳区,U2符合该沉积层特征,由此确定T2为漂移期沉积的顶界面。
Tb界面从268 km处向右分开的下支界面被标定为T4,其反射强度和连续性均不及上述T2、T3界面。从地震剖面可见,此T4界面可从268 km处向右直至剖面末端,即268~500 km处均可追踪。值得注意的是,在306~440 km区段,T4界面双程走时逐渐加深,直至8.2 s附近,从440 km向右,双程走时由于断层作用而抬升至4 s附近。T4实际上为声波基底,其下为地壳基底,上为陆缘沉积充填的起始。考虑到冈瓦纳基底的背景,结合相邻陆缘盆地油气勘探揭示的盆地基底,本文将T4界面确定为晚石炭世冈瓦纳陆内裂谷时期的开始。长期的中生代陆内裂谷作用,使T4与上覆侏罗世界面之间沉积了巨厚的地层[23-24]
在马达加斯加陆缘盆地中,T1界面为最显著的一条强反射界面,与相邻盆地油气勘探资料[3,28-30]的对比分析表明,该界面是早、晚白垩世的分界面。地球化学与岩石学数据分析表明,此时马里昂热点产生了剧烈的岩浆活动[8-9],并形成大量熔岩流,使沉积环境发生剧烈变化,导致了该界面在地震剖面上的强反射和显著发育。同时,由于该界面的屏蔽效应,使得下伏地层的地震反射明显弱化。T1界面在整个研究区都能够连续识别,表现为强振幅、连续的地震反射特征。在剖面的0~268 km,T1双程走时在7~8 s附近;在268~500 km,双程走时不断抬升,最高至1.2 s。受戴维断裂带走滑运动的影响,剖面中洋壳部分的年龄大致为距今120 Ma,洋盆中不存在T2界面,因此界面T1下伏U3与U4两组不同地层。其中,U3以界面T2为底界面,存在于陆缘部分,为后漂移沉积,地震反射表现为平行-亚平行的局部波状反射,连续性较好,中振幅,存在部分空白区域;U4以界面Tb为底界面,存在于洋盆部分,厚度薄于U3,内部表现为亚平行状、局部嘈杂状反射,连续性中等-差,振幅低-中等,局部区域存在空白反射。
T0为现今海底面,是沉积地层和海水的分界面。在二维地震剖面上(由陆向洋方向),T0界面不断加深,在洋盆区域变化不明显,坡度较小;在洋陆边界附近抬升明显,表明陆架存在明显陡坡。T0下伏地层U5,以T1为底界面,地震反射特征表现为平行-亚平行,连续性较好,振幅较低,受T1界面的强振幅屏蔽,部分区域呈现空白反射。

2.2.2 重震联合反演

晚古生代以来,马达加斯加岛西南部陆缘经历了多期构造活动。基于地震剖面的构造地层解释,以东西冈瓦纳大陆裂解、马达加斯加岛停止漂移、被动陆缘稳定发育为时间节点,将马达加斯加岛西南部陆缘演化过程划分为晚石炭世-早侏罗世、中侏罗世-早白垩世及晚白垩世以来三个阶段[4,23,31-32]。根据测线地层界面,对应地可划分4个沉积层,其中沉积层2与沉积层1皆为被动陆缘阶段沉积层。由于缺乏马达加斯加岛西南部陆缘深部钻井岩芯资料,不能获得裂陷期和更古老地层准确的岩石密度。因此,本文通过对前人研究资料[6,28-29,33]进行合理推测来设置研究区重震联合反演所需要的物性参数,具体的参数见表2
表2 马达加斯加岛西南部陆缘地壳密度与速度模型[7,28-29]

Tab.2 Crustal density model and velocity model of the southwestern margin of Madagascar Island[7,28-29]

地层 密度/(g·cm-3) 层速度/(m·s-1)
海水 1.03 1 500
沉积层1 2.20 2 000
沉积层2 2.30~2.45 3 300
沉积层3 2.52~2.56 3 500
沉积层4 2.57~2.64 4 000
洋壳 2.93 6 500
上地壳 2.67~2.74 5 700~6 300
下地壳 2.90~2.94 6 400~7 000
下地壳高密度体 3.10~3.16 7 000~7 500
地幔 3.30 8 000~8 100
为实现重震联合反演,需要对地震反射时间剖面进行时深转换,形成深度剖面。由于本区域缺少层位速度模型,因此以多道地震剖面中识别的7个界面(莫霍面、Tb、T0、T1、T2 、T3、T4)及其年龄为基础,采用前人在莫桑比克北部与中部陆缘深部地震勘探得到的速度模型[7,28-29],推断本研究剖面中各个地层的速度(表2)。

3 结果

3.1 区域地形和重力特征

图4展示了马达加斯加岛西南部的卫星测高重力异常和地形。自西向东,从莫桑比克海盆向马达加斯加岛和马达加斯加海台方向,重力异常值逐渐增加,整体呈现东高西低的特征(图4a)。马达加斯加岛内、马达加斯加海台和西侧的戴维断裂带表现为重力异常正高值,莫桑比克海盆区表现为重力异常负低值。最为显著的是在莫桑比克海盆与马达加斯加岛之间存在一条重力异常负低值带,以25°S—27°S之间的椭圆形高重力异常区(图4中的“H”)为界分为南北两个区段,北侧为NNW向延伸的北部区段,南侧为NNE向延伸区段,条带内部重力异常幅值整体小于-100 mGal。条带东侧为重力异常高值区,西侧为重力异常低值区或过渡区。
图4 马达加斯加岛西南陆缘自由空间重力异常(a)与地形图(b)

Fig.4 The free-space gravity anomaly (a) and topographic map (b) of the southwestern margin of Madagascar Island

图4可见,地形与重力异常在大部分区域有较好的对应关系,但在莫桑比克海盆与马达加斯加岛之间的重力异常负低值带内不存在坳陷和隆起地形,因此可以大致排除其受海沟地形的影响,同时椭圆形的高重力异常区H也并非海山的反应。前人资料[29]显示,莫桑比克海盆区域发育明显的NE向条带磁力异常,这些条带磁力异常终止于上述重力异常负低值带的西侧边界,且内部不存在明显的条带状磁力异常。因此,该重力异常显著负低值带的地壳性质不明,可能是地壳的密度变化或深部存在异常体而导致的。

3.2 CMD01地震反射特征

图5为测线CMD01多道地震剖面和对应的实测重力、磁力异常。图6图5中洋陆过渡区地震剖面的局部放大。根据前文所述,在268 km这个位置Tb 被截断,向陆方向分支为两个不整合面,分别为向上抬升的T3和向下延伸的T4两条分支界面。268 km处还发育了近直立的断层,该断层可能为转换断层戴维断裂带的马尾状延伸。基于上述特征,以该断层作为洋陆边界,左侧为洋盆部分,右侧为马达加斯加陆缘部分。向陆方向的重力异常先升高后逐渐下降,洋盆内重力异常整体无大幅度变化,受深部构造影响存在较小起伏。图7显示,剖面的26~68 km段,T1呈现下陷趋势,重力异常值由东向西减小,位置与下陷区域一致,判断该位置为拗陷构造;剖面75~87 km,重力异常曲线隆起,并且二维地震剖面显示有明显的岩浆活动痕迹,在深部形成的侵入体刺穿T1并破坏Tb的连续性。图8显示,剖面290~315 km可识别出多条走滑断层,断层右侧(向陆一侧)地形与地层出现抬升,此时重力异常急剧升高,磁力异常迅速降低后迅速升高至>400 nT。由于该区段内无明显岩浆活动痕迹,因此认为影响重力、磁力异常曲线变化的物质可能位于地壳深部。在306 km处,漂移期沉积U2被截断,T2界面尖灭。从320 km向右,重力异常不增反降直至425 km,425~500 km,重力异常曲线上升(图5c),根据二维地震剖面推断是受到了地形和基底抬升的影响。其中,在470~485 km的地形抬升处可识别出代表走滑断层的花状结构,且断层两侧地层倾向相反,表明该位置疑似存在早期压扭走滑,晚期转为张扭的过程,使地形发生改变,并生成地堑结构(图9)。
图5 CMD01的二维地震剖面(a)、磁力异常曲线(b)和自由空间重力异常曲线(c)

Fig.5 CMD01 2D seismic profile (a), magnetic anomaly curve (b) and gravity anomaly curve (c)

图6 测线CMD01洋陆边界的地震剖面

Fig.6 Continental-ocean boundary in CMD01 seismic profile

图7 测线CMD01海盆内拗陷区与岩浆侵入的地震剖面

Fig.7 Rift area and magmatic intrusions within the basin in CMD01 seismic profile

图8 测线CMD01陆缘区地层抬升的地震剖面

Fig.8 Uplift stratigraphic in CMD01 seismic profile

图9 测线CMD01陆缘区地堑构造的地震剖面

Fig.9 Graben structure in CMD01 seismic profile

3.3 重震联合反演

对基于时深转换得到的深度剖面进行重震联合反演,得到地壳密度模型(以下简称模型),并对其反演的莫霍面深度进行评估。将反演莫霍面(绿色实线)与时深转换后的莫霍面(蓝色虚线)进行对比,发现两条线在起伏趋势上基本相似,但在300~400 km存在最大约3 km的深度差距(图10c),而在该区域,重力异常曲线(图10b)也表现出与两侧明显不同,存在明显的升高趋势,因此,猜想在该区域深部可能存在异常体。为了验证该猜想,在不改变原模型其他反演层位深度以及地层密度的前提下,在300~400 km之间添加一层下地壳高密度体作为新模型,新模型反演的莫霍面深度(图10e)与基于时深转换的莫霍面深度趋于一致。图10d显示,重力异常实测值(黑色实线)为-170 ~-113 mGal,由陆向洋方向上变化明显,在洋壳区域和陆壳区域呈现明显的分带特征。新模型的拟合误差为0.727 mGal。结合CMD01地壳密度模型和自由空间重力异常,马达加斯加西南陆缘可以划分为减薄陆壳区、洋陆过渡带、洋壳区三个陆缘构造单元。
图10 CMD01剖面地壳密度模型界面反演结果

(D表示密度。)

Fig.10 The inversion results of the crustal density model interface for CMD01

(D represents density.)

1)减薄陆壳区
减薄陆壳区位于模型303~500 km之间,重力异常变化强烈,呈两段式变化(图10d)。并且,图10e的地壳结构中显示,在425~500 km之间不存在地壳高密度体,而在303~425km之间存在下地壳高密度体。因此根据重力变化和地壳结构的不同将减薄陆壳区域划分为两部分进行讨论。
在模型的425~500 km之间,重力异常值的变化范围在-112 ~-170 mGal之间,形成了由陆向洋方向的第一个重力异常峰值区,对应莫霍面埋深变化明显,由25 km抬升至16 km。由于结晶基底向洋方向缓慢下降,地壳厚度由22 km减薄至8 km,属于典型的减薄型陆壳。各地层厚度在该区域内的变化不同,沉积层1在470~500 km之间被完全剥蚀;沉积层2厚度变化不明显;沉积层3无明显变化;沉积层4逐渐增厚。
在模型的303~425 km之间,重力异常值在-115~-170 mGal之间变化,是第二个重力异常峰值区,曲线变化坡度相比第一个峰值区较缓。通过模型拟合发现,导致重力异常抬升的原因是在该区域下地壳中存在高密度体,该高密度体位于模型300~400 km,厚度最大为2 km。在该区域的磁力异常曲线上也有较明显的变化(图10a)。莫霍面深度由16 km抬升至14 km,地壳厚度由8 km减薄至6 km,沉积层3和沉积4厚度逐渐变小。
2)洋陆过渡带
洋陆过渡带位于模型268~303 km,在洋壳和减薄陆壳之间,重力异常曲线先减小后增加形成了碗状结构。相较于减薄陆壳区,莫霍面变化不明显,地壳厚度大致为7 km,其与洋壳以断层为界,与减薄陆壳以一个大型走滑断层带为界。
3)洋壳区
洋壳区位于模型左侧25~268 km之间,重力异常的变化范围在-150 ~-130 mGal之间,变化幅度与陆壳相比明显较缓,其中75~85 km之间隆起,该区域对应二维地震剖面上有明显的岩浆活动痕迹(图7)。洋壳区磁力异常较减薄陆壳区平缓,表现为明显的洋壳特征,洋壳厚度在6~8 km之间。

4 讨论

4.1 下地壳高密度体成因

地壳密度模型,由浅到深分别为海水层、沉积层、地壳层与地幔层,模型中的侵入体为侵入岩与下地壳高密度体。
下地壳中高密度异常体有许多成因。考虑到本研究区处于被动大陆边缘,多强调两个成因,即:海水沿贯穿薄化地壳的深大断裂向下渗透,使上地幔发生蛇纹石化作用[34-35]和地幔上涌后,岩浆侵入下地壳[36-37]
CMD01剖面地壳密度模型显示,地壳的厚度最薄仅为7 km,而一般正常的地壳厚度为30 km左右,已有的研究认为在这种强烈的减薄陆壳中一般存在拆离断层作用[38-39]。马达加斯加岛陆缘在发生漂移运动前经历过斜向拉张,而该异常体大致发育在减薄陆壳区且靠近走滑断层,因此,走滑、拉张背景下发育的拆离断层可能将上地幔剥露,导致海水与地幔的橄榄岩发生化学作用形成蛇纹石化橄榄岩。并且,走滑断层作为深大断裂,也有可能使海水沿该断裂下渗至上地幔中,形成蛇纹石化上地幔[40]。然而,蛇纹石化作用形成的异常地幔一般位于地壳之下,而反演的结果表明该高密度体侵入于正常地壳之间,与实际情况不符。
所以,我们更倾向于认为CMD01剖面地壳密度模型中的异常体为岩浆底侵形成的下地壳高密度体,主要的原因分析如下。侏罗纪时期,冈瓦纳大陆受到地幔柱的显著影响开始裂离,同时伴随着强烈的岩浆活动,形成了典型的火山型陆缘,如莫桑比克中部陆缘,其洋陆过渡带中发育有向海倾斜的反射层和下地壳高速体[29,41];根据马达加斯加岛西部的穆伦达瓦盆地的深部地震资料,发现在洋陆过渡带和陆壳区域内存在下地壳高速体,最大厚度可达4.1 km[28];共轭陆缘(坦桑尼亚-北莫桑比克陆缘)二维地震剖面深部的反射结构与重震联合反演的地壳密度模型都表明其洋陆过渡带深部发生了岩浆底侵增生[6]。目前发现的大面积岩浆活动主要集中于莫桑比克中部与南部,而坦桑尼亚陆缘洋壳下方的底侵表明马达加斯加岛陆缘在裂离过程中同样存在岩浆活动,只不过由于该处岩浆活动的强度不及莫桑比克中部陆缘而未形成火成岩省或向海倾斜的反射层。

4.2 陆缘结构单元的地壳特征

根据实测重力数据和反演结果(图10)可知,减薄陆壳区域位于303~500 km之间,约为197 km,范围较宽,地壳底部岩浆侵入形成厚度约为2 km左右的下地壳高密度体,重力异常变化幅度大,呈现多个峰值,莫霍面埋深与地壳厚度的变化明显;洋陆过渡带位于268~303 km内,宽度约为35 km,向洋方向重力异常曲线逐渐平缓,莫霍面埋深无明显变化,与洋壳和陆壳之间以断层为界;洋壳区位于25~268 km,重力异常变化小,地壳厚度变化约为2 km,重力、磁力异常呈现明显的洋壳特征。
以上陆缘结构单元的地壳特征表明,马达加斯加陆缘的洋陆过渡带较窄,没有出现贫岩浆型被动大陆边缘宽阔的地幔剥露[34,42-43],或富岩浆型被动大陆边缘向海倾斜的火山岩层、岩浆底侵形成的下地壳高速体和大量火成岩体[44-45]。受控于走滑断裂带,紧邻转换断层的洋壳厚度相对于正常洋壳明显较薄,因此属于转换陆缘,其发育受戴维断裂带走滑运动的影响,是在马达加斯加岛与非洲大陆分离的过程中,在从拉张运动转变为走滑运动的时期发育形成的。陆壳区薄且范围较宽,推断可能是受拉张运动的影响,原因如下:1)在多道地震剖面中减薄陆壳部分能够识别出的走滑断层多为小型断层且较为模糊,缺少确切的大型走滑断层,不符合走滑运动主导的减薄陆壳特征;2)陆缘演化过程表明,在拉张运动转变为走滑运动时,大部分陆缘已发育完全,因此戴维断裂带的走滑运动主要破坏洋壳与洋陆过渡带,对陆壳的影响较小[18,46]

4.3 马达加斯加岛西南部陆缘演化

根据板块重建可知,在冈瓦纳大陆裂离之前(距今约182 Ma),马达加斯加岛西南部陆缘与莫桑比克北部陆缘是一对共轭陆缘(图11a)[2,20-22,27]
图11 马达加斯加岛西南部陆缘演化示意图

(图件据文献[19]和[53]修改。)

Fig.11 Evolution map of the southwestern margin of Madagascar Island (Figure is modified from references [19] and [53].)

在距今164 Ma,马达加斯加西部陆缘首先受到东、西冈瓦纳大陆张裂控制,形成了张裂型陆缘[2,47-49],本研究中的减薄陆壳区也在此张裂过程中形成。大量的研究表明沿戴维断裂带的走滑运动是形成转换型陆缘的重要机制[6-7,19,28,50-51]。REEVES 等[2]在坦桑尼亚陆缘南部的二维地震剖面中发现戴维断裂带的东西两侧分别存在不同时期的洋壳,认为东、西冈瓦纳大陆最早的张裂运动方向是NW—SE向,伸展方向大致垂直于断裂带,在发生马达加斯加岛漂移运动前,先行的裂谷运动已经形成了扩张中心与洋壳。SINHA 等[6]通过多道地震剖面在坦桑尼亚近海找到了疑似裂谷时期的扩张中心,佐证了该观点(图11b)。
在距今约153 Ma,坦桑尼亚陆缘由斜向伸展过渡为转换运动,戴维断裂带由陆内的走滑断层转变为转换断层[2,6],影响了原本正常发育的张裂型洋陆过渡带,随之形成了本研究中张裂型与转换型并存的洋陆过渡带(图11c)。VORMANN等[7,28]根据深部地震探测数据解释分析,认为该陆缘由陆内岩石圈斜向伸展后沿着转换断层的走滑运动形成。中侏罗世时期,马达加斯加岛与非洲大陆以NW—SE向裂离,并且在海底扩张形成的索马里海盆中识别出与最老洋壳对应的磁力异常条带,为NE—SW走向,证实了北部陆缘在受走滑运动影响之前初始裂离方向为NW—SE向[6-7,28,52]
最终,于距今120 Ma,索马里海盆停止扩张,马达加斯加岛停止漂移,其西部陆缘基本形成(图11d),洋壳也在该阶段形成[5-6,48,53]。本研究中的地壳密度模型存在着较宽的减薄陆壳与较窄的洋陆过渡带,表明该陆缘前期经历伸展运动,后期受走滑断裂影响。

5 结论

通过对马达加斯加岛西南部陆缘的实测重力、地震剖面构造地层解释和在此基础上的重震联合反演研究,结合磁力异常特征综合分析,得到以下几点结论。
1)剖面CMD01的地壳密度模型内存在下地壳高密度体,与重力异常剖面中的椭圆高重力异常区对应,底侵的岩浆来源于冈瓦纳大陆裂解初期的岩浆活动,不同陆缘的岩浆活动程度不同。
2)CMD01剖面根据重力异常曲线、地层和地壳厚度等信息,由洋向陆划分为洋壳、洋陆过渡带、减薄陆壳三个区域。其中,洋陆过渡带位于两者之间,宽度大致为35 km,地壳厚度薄至7 km,推测是受戴维断裂带走滑运动影响,具备部分转换型洋陆过渡带的特征。
3)地壳密度模型内减薄陆壳宽度约为197 km,洋陆过渡带宽度约为35 km,较窄的洋陆过渡带与较广的减薄陆壳表明,马达加斯加岛西南部陆缘在演化过程中存在拉张运动向走滑运动的转变。
[1]
张光亚, 刘小兵, 赵健, 等. 东非被动大陆边缘盆地演化及大气田形成主控因素:以鲁武马盆地为例[J]. 地学前缘, 2018, 25(2):24-32.

DOI

ZHANG G Y, LIU X B, ZHAO J, et al. Passive continental margin basin evolution of East Africa and the main controlling factors of giant gas fields: An example from the Rovuma Basin[J]. Earth Science Frontiers, 2018, 25(2): 24-32.

DOI

[2]
REEVES C V, TEASDALE J P, MAHANJANE E S. Insight into the Eastern Margin of Africa from a new tectonic model of the Indian Ocean[J]. Geological Society, London, Special Publications, 2016, 431(1): 299-322.

[3]
温志新, 王兆明, 宋成鹏, 等. 东非被动大陆边缘盆地结构构造差异与油气勘探[J]. 石油勘探与开发, 2015, 42(5):671-680.

WEN Z X, WANG Z M, SONG C P, et al. Structural architecture difference and petroleum exploration of passive continental margin basins in East Africa[J]. Petroleum Exploration and Development, 2015, 42(5): 671-680.

[4]
许志刚, 韩文明, 孙玉梅. 东非大陆边缘构造演化过程与油气勘探潜力[J]. 中国地质, 2014, 41(3):961-969.

XU Z G, HAN W M, SUN Y M. Tectonic evolution and petroleum exploration prospect of East Africa[J]. Geology in China, 2014, 41(3): 961-969.

[5]
COFFIN M F, RABINOWITZ P D. Reconstruction of Mada-gascar and Africa: Evidence from the Davie fracture zone and western Somali Basin[J]. Journal of Geophysical Research, 1987, 92(B9): 9385-9406.

[6]
SINHA S T, SAHA S, LONGACRE M, et al. Crustal archi-tecture and nature of continental breakup along a transform margin: New insights from Tanzania-Mozambique margin[J]. Tectonics, 2019, 38(4): 1273-1291.

[7]
VORMANN M, JOKAT W. The crustal structure of the Kerimbas Basin across the offshore branch of the East African Rift System[J]. Geophysical Journal International, 2021, 226(3): 2073-2102.

[8]
MEERT J G, TAMRAT E. Paleomagnetic evidence for a stationary Marion hotspot: Additional paleomagnetic data from Madagascar[J]. Gondwana Research, 2006, 10(3/4): 340-348.

[9]
STOREY M, MAHONEY J J, SAUNDERS A D, et al. Timing of hot spot: Related volcanism and the breakup of Madagascar and India[J]. Science, 1995, 267(5199): 852-855.

[10]
JOURDAN F, BERTRAND H, FÉRAUD G, et al. Litho-spheric mantle evolution monitored by overlapping large igneous provinces: Case study in southern Africa[J]. Lithos, 2009, 107(3/4): 257-268.

[11]
CATUNEANU O, WOPFNER H, ERIKSSON P G, et al. The Karoo Basins of south-central Africa[J]. Journal of African Earth Sciences, 2005, 43(1/2/3): 211-253.

[12]
COX K G. Karoo igneous activity, and the early stages of the break-up of Gondwanaland[J]. Journal of the Geological Society, 1992, 68(1): 137-148.

[13]
EAGLES G, HOANG H H. Cretaceous to present kinematics of the Indian, African and Seychelles plates[J]. Geophysical Journal International, 2014, 196(1): 1-14.

[14]
COLLIER J S, SANSOM V, ISHIZUKA O, et al. Age of Seychelles-India break-up[J]. Earth and Planetary Science Letters, 2008, 272(1/2): 264-277.

[15]
GEIGER M, CLARK D N, METTE W. Reappraisal of the timing of the breakup of Gondwana based on sedimentological and seismic evidence from the Morondava Basin, Madagascar[J]. Journal of African Earth Sciences, 2004, 38(4): 363-381.

[16]
覃阳亮, 何幼斌, 蔡俊, 等. 东非海岸Davie构造带的构造演化特征及其成因机制[J]. 岩性油气藏, 2021, 33(2):104-115.

QIN Y L, HE Y B, CAI J, et al. Tectonic evolution and formation mechanism of Davie Fracture Zone in East Africa coast[J]. Lithologic Reservoirs, 2021, 33(2): 104-115.

[17]
KLIMKE J, FRANKE D, MAHANJANE E S, et al. Tie points for Gondwana reconstructions from a structural inter-pretation of the Mozambique Basin, East Africa and the Riiser-Larsen Sea, Antarctica[J]. Solid Earth, 2018, 9(1): 25-37.

[18]
REEVES C. The position of Madagascar within Gondwana and its movements during Gondwana dispersal[J]. Journal of African Earth Sciences, 2014, 94: 45-57.

[19]
KLIMKE J, FRANKE D. Gondwana breakup: No evidence for a Davie Fracture Zone offshore northern Mozambique, Tanzania and Kenya[J]. Terra Nova, 2016, 28(4): 233-244.

[20]
MAHANJANE E S. A geotectonic history of the northern Mozambique Basin including the Beira High—A contribution for the understanding of its development[J]. Marine and Petroleum Geology, 2012, 36(1): 1-12.

[21]
LEINWEBER V T, JOKAT W. The Jurassic history of the Africa-Antarctica corridor—New constraints from magnetic data on the conjugate continental margins[J]. Tectonophysics, 2012, 530/531: 87-101.

[22]
KÖNIG M, JOKAT W. The Mesozoic breakup of the Weddell Sea[J]. Journal of Geophysical Research: Solid Earth, 2006, 111: B12102.

[23]
SCHANDELMEIER H, BREMER F, HOLL H G. Kinematic evolution of the Morondava rift basin of SW Madagascar—From wrench tectonics to normal extension[J]. Journal of African Earth Sciences, 2004, 38(4): 321-330.

[24]
SALMAN G, ABDULA I. Development of the Mozambique and Ruvuma sedimentary basins, offshore Mozambique[J]. Sedimentary Geology, 1995, 96(1/2): 7-41.

[25]
NAIRN A E M, LERCHE I, ILIFFE J E. Geology, basin analysis, and hydrocarbon potential of Mozambique and the Mozambique Channel[J]. Earth-Science Reviews, 1991, 30(1/2): 81-123.

[26]
WARNER M R. Seismic reflections from the Moho—The effect of isostasy[J]. Geophysical Journal International, 1987, 88(2): 425-435.

[27]
KLIMKE J, FRANKE D, GAEDICKE C, et al. How to identify oceanic crust—Evidence for a complex break-up in the Mozambique Channel, off East Africa[J]. Tectono-physics, 2016, 693: 436-452.

[28]
VORMANN M, FRANKE D, JOKAT W. The crustal structure of the southern Davie Ridge offshore northern Mozambique—A wide-angle seismic and potential field study[J]. Tectonophysics, 2020, 778: 228370.

[29]
MUELLER C O, JOKAT W, SCHRECKENBERGER B. The crustal structure of Beira High, central Mozambique—Combined investigation of wide-angle seismic and potential field data[J]. Tectonophysics, 2016, 683: 233-254.

[30]
于璇, 侯贵廷, 代双河, 等. 东非大陆边缘构造演化与油气成藏模式探析[J]. 地质科技情报, 2015, 34(6):147-154,158.

YU X, HOU G T, DAI S H, et al. Tectonic evolution and hydrocarbon pooling patterns analysis in East Africa continental margin[J]. Bulletin of Geological Science and Technology, 2015, 34(6): 147-154, 158.

[31]
黄国平, 胡清乐, 陈冬明, 等. 马达加斯加地质矿产概况[J]. 资源环境与工程, 2014, 28(5):626-632.

HUANG G P, HU Q L, CHEN D M, et al. General situation of geology and mineral resources in Madagascar[J]. Resources Environment & Engineering, 2014, 28(5): 626-632.

[32]
COLLINSA S, WINDLEY B F. The tectonic evolution of central and northern Madagascar and its place in the final assembly of Gondwana[J]. The Journal of Geology, 2002, 110(3): 325-339.

[33]
RYZHOVA D A, KOSNYREVA M V, DUBININ E P, et al. The structure tectonosphere of the Mozambique and Mada-gascar ridges by geophysical data[J]. Moscow University Geology Bulletin, 2022, 77(1): 18-27.

[34]
FUNCK T, HOPPER J R, LARSEN H C, et al. Crustal structure of the ocean-continent transition at Flemish Cap: Seismic refractionresults[J]. Journal of Geophysical Research, 2003, 108: 2531.

[35]
BOILLOT G, GRIMAUD S, MAUFFRET A, et al. Ocean-continent boundary off the Iberian margin: Aserpentinite diapir west of the Galicia Bank[J]. Earth and Planetary Science Letters, 1980, 48(1): 23-34.

[36]
HIRSCH K K, BAUER K, SCHECK-WENDEROTH M. Deep structure of the western South African passive margin—Results of a combined approach of seismic, gravity and isostatic investigations[J]. Tectonophysics, 2009, 470(1/2): 57-70.

[37]
MINSHULL T A, LANE C I, COLLIER J S, et al. The relationship between rifting and magmatism in the north-eastern Arabian Sea[J]. Nature Geoscience, 2008, 1: 463-467.

[38]
MASINI E, MANATSCHAL G, MOHN G. The Alpine Tethys rifted margins: Reconciling old and new ideas to understand the stratigraphic architecture of magma-poor rifted margins[J]. Sedimentology, 2013, 60(1): 174-196.

[39]
PÉRON-PINVIDIC G, MANATSCHAL G. From microconti-nents to extensional allochthons: Witnesses of how continents rift and break apart?[J]. Petroleum Geoscience, 2010, 16(3): 189-197.

[40]
MAIA M, SICHEL S, BRIAIS A, et al. Extreme mantle uplift and exhumation along a transpressive transform fault[J]. Nature Geoscience, 2016, 9: 619-623.

[41]
VOLKER T L, FRAUKE K, SÖNKE N, et al. The crustal structure of the Central Mozambique continental margin—Wide-angle seismic, gravity and magnetic study in the Mozambique Channel, Eastern Africa[J]. Tectonophysics, 2013, 599: 170-196.

[42]
HUISMANS R, BEAUMONT C. Depth-dependent extension, two-stage breakup and cratonic underplating at rifted margins[J]. Nature, 2011, 473: 74-78.

[43]
AUTIN J, LEROY S, BESLIER M O, et al. Continental break-up history of a deep magma-poor margin based on seismic reflection data (northeastern Gulf of Aden margin, offshore Oman)[J]. Geophysical Journal International, 2010, 180(2): 501-519.

[44]
FRANKE D. Rifting, lithosphere breakup and volcanism: Comparison of magma-poor and volcanic rifted margins[J]. Marine and Petroleum Geology, 2013, 43: 63-87.

[45]
GEOFFROY L. Volcanic passivemargins[J]. Comptes Rendus Géoscience, 2005, 337(16): 1395-1408.

[46]
SINHA S T, NEMČOK M, CHOUDHURI M, et al. The role of break-up localization in microcontinent separation along a strike-slip margin: The East India-Elan Bank case study[J]. Geological Society, London, Special Publications, 2016, 431(1): 95-123.

[47]
PHETHEAN J J J, KALNINS L M, VAN HUNEN J, et al. Madagascar’s escape from Africa: A high-resolution plate reconstruction for the Western Somali Basin and implications for supercontinent dispersal[J]. Geochemistry, Geophysics, Geosystems, 2016, 17(12): 5036-5055.

[48]
DAVISON I, FAULL T, GREENHALGH J, et al. Trans-pressional structures and hydrocarbon potential along the Romanche Fracture Zone: A review[J]. Geological Society, London, Special Publications, 2016, 431(1): 235-248.

[49]
MUELLER C O, JOKAT W. The initial Gondwana break-up: A synthesis based on new potential field data of the Africa-Antarctica Corridor[J]. Tectonophysics, 2019, 750: 301-328.

[50]
FRANKE D, JOKAT W, LADAGE S, et al. The offshore East African Rift System: Structural framework at the toe of a juvenile rift[J]. Tectonics, 2015, 34(10): 2086-2104.

[51]
SAUTER D, UNTERNEHR P, MANATSCHAL G, et al. Evidence for magma entrapment below oceanic crust from deep seismic reflections in the Western Somali Basin[J]. Geology, 2016, 44(6): 407-410.

[52]
GAINA C, TORSVIK T H, VAN HINSBERGEN D J J, et al. The African Plate: A history of oceanic crust accretion and subduction since the Jurassic[J]. Tectonophysics, 2013, 604: 4-25.

[53]
VORMANN M, JOKAT W. Crustal variability along the rifted/sheared East African margin: A review[J]. Geo-Marine Letters, 2021, 41(2): 19.

文章导航

/