岩石学报  2019, Vol. 35 Issue (10): 3065-3082, doi: 10.18654/1000-0569/2019.10.07   PDF    
西藏松多地区早侏罗世变质辉长岩的成因及其构造意义
解超明1,2, 段梦龙1, 于云鹏1, 王斌1, 宋宇航1, 张红雨3     
1. 吉林大学地球科学学院, 长春 130061;
2. 东北亚矿产资源评价自然资源部重点实验室, 长春 130061;
3. 中国地质大学, 北京 100083
摘要: 冈底斯岩浆岩带中东部松多地区新发现的变质辉长岩对于讨论松多古特提斯洋演化历史具有重要约束意义。在详细野外地质调查基础上,本文报道了松多变质辉长岩的LA-ICP-MS锆石U-Pb年龄、锆石Hf同位素成分分析和全岩地球化学组成,并探讨了岩石的成因及其对松多古特提斯洋构造演化的启示。结果表明,松多变质辉长岩原岩形成于195±6Ma~202±2Ma之间;SiO2含量为45.56%~53.69%,MgO含量较高为6.04%~11.75%,Mg#为66~77;高Cr(209×10-6~746×10-6)和Ni(80×10-6~234×10-6)含量。该岩石富集轻稀土元素,平滑右倾的稀土配分模式介于E-MORB和OIB之间;富集大离子亲石元素和Nb、Ta等高场强元素;具有较低的(La/Yb)N=2.06~4.81、(La/Sm)N=1.50~2.48、(Gd/Yb)N=1.21~1.67和(Ce/Y)N < 2。锆石εHft)值为-1.9~3.9,具相对古老Hf模式年龄(988~1361Ma)。综合对区域地质、年代学、岩石地球化学以及锆石Hf同位素等资料的全面分析,认为松多早侏罗世变质辉长岩原岩的岩浆可能起源于下地壳改造的亏损岩石圈地幔的部分熔融,形成于松多古特提斯洋壳板片断离事件所引发岩石圈减薄的伸展构造背景。
关键词: 青藏高原    松多地区    变质辉长岩    岩石成因    早侏罗世    
Genesis and geological significance of Early Jurassic metamorphic gabbro in the Sumdo area, Tibet
XIE ChaoMing1,2, DUAN MengLong1, YU YunPeng1, WANG Bin1, SONG YuHang1, ZHANG HongYu3     
1. College of Earth Science, Jilin University, Changchun 130061, China;
2. MNR Key Laboratory of Mineral Resources Evaluation in Northeast Asia, Changchun 130061, China;
3. China University of Geosciences, Beijing 100083, China
Abstract: The metamorphic gabbro newly discovered in the Sumdo area, the eastern part of the Gangdese magmatic belt, is of importance for discussing the evolution history of the Sumdo Paleo-Tethyan Ocean. In this paper, on the basis of detailed field geological investigation, we report LA-ICP-MS zircon U-Pb ages, zircon Hf isotopic and whole-rock major and trace element compositions from the metamorphic gabbro, with the aims of understanding their petrogenesis and the implications for Sumdo Paleo-Tethyan Ocean evolution. LA-ICP-MS zircon U-Pb isotopic dating indicates that the gabbros emplaced at 195±6Ma~202±2Ma. It is characterized by its contents of variable SiO2 (45.56%~53.69%), and high MgO (6.04%~11.75%), Mg# (66~77), Cr (209×10-6~746×10-6) and Ni (80×10-6~234×10-6). The REE distribution patterns are between E-MORB and OIB. They are characterized by smooth REE patterns curves inclining gently to right with LREE slightly enriched. The large ion lithophile elements (LILEs) and the high field strength element (HFSEs) (such as Nb, Ta) are riched. They have low (La/Yb)N=2.06~4.81, (La/Sm)N=1.50~2.48, (Gd/Yb)N=1.21~1.67 and (Ce/Y)N < 2. Lu-Hf isotope study demonstrates that its εHf(t) is between -1.9 and 3.9, and the value of ancient zircon Hf crustal modal ages changes greatly from 988Ma to 1361Ma. Combining with regional geology, geochronolagical, geochemical and zircon Hf isotopic data, we propose that the primary of Early Jurassic Sumdo Metamorphic gabbro was derived by partial melting of a depleted lithospheric mantle modified by ancient lower crustal material. They were formed in an extensional tectonic environment related to the lithospheric thinning by slab breakoff of Sumdo Paleo-Tethyan Ocean.
Key words: Tibetan Plateau    Sumdo    Metamorphic gabbro    Petrogenesis    Early Jurassic    

冈底斯构造-岩浆岩带作为青藏高原的重要组成部分,近东西向夹持于印度河-雅鲁藏布江结合带和班公湖-怒江结合带之间(图 1a),记录了特提斯洋俯冲消减、印度板块与欧亚板块碰撞、青藏高原隆升等一系列极其复杂的地质过程(肖序常和王军,1998潘桂棠等,2006)。冈底斯中部松多榴辉岩的发现(杨经绥等,2006),使松多地区的基础地质研究得到研究者的广泛关注(许志琴等,2013Cheng et al., 2015Weller et al., 2016Cao et al., 2017Yu et al., 2018Wang et al., 2019)。有学者认为拉萨地体中“松多-墨竹工卡”蛇绿岩和榴辉岩呈外来岩块群产出于松多岩组的浅变质岩系中,可能为异地推覆体(Liu et al., 2009Zhang and Tang, 2009);随着研究者在松多地区陆续发现与古大洋演化相关的蛇绿岩(王斌等,2017)、二叠纪洋岛(Wang et al., 2019)、榴辉岩及其变质作用(杨经绥等,2006Cheng et al., 2015李鹏等,2017),逐渐揭示在松多地区可能存在一条晚古生代的古特斯缝合带,其闭合时代为晚三叠世末期(李化启等,2008),但是关于其闭合后的构造-岩浆作用却未见报道,鲜有讨论。

图 1 冈底斯东段构造简图(a、b)及松多地区地质图(c) Fig. 1 Sketch map of the eastern Gangdise belt (a, b) and geological map of the Sumdo area of Tibet showing the localities of Early Jurassic metamorphic gabbro (c)

冈底斯地区以发育巨型的构造-岩浆活动为特色,也被称之为冈底斯岩浆岩带,发育的岩浆作用至少可以划分出4个阶段:220~152Ma、109~80Ma、65~41Ma、33~13Ma(Ji et al., 2012Yu et al., 2018),其中220~152Ma这一阶段的岩浆作用与古特提斯洋闭合、新特提斯洋开启、俯冲消减最密切相关,具有重要的科学意义(Chu et al., 2006; Zhu et al., 2008朱弟成等,2008Kang et al., 2014董汉文等,2016),但研究程度却相对较低。冈底斯晚三叠世-早侏罗世的岩浆活动主要集中分布于东段(图 1b),具有弧岩浆作用地球化学特征。黄丰等(2015)认为早侏罗世叶巴组和桑日群火山岩是雅鲁藏布江洋早期俯冲的岩浆产物。刘琦胜等(2006)根据冈底斯中部宁中早侏罗世白云母二长花岗岩的研究认为其属于冈底斯印支造山旋回晚期碰撞阶段的产物。也有学者推测具有侏罗纪岩浆作用可能形成于班公湖-怒江洋南向俯冲消减的结果(Geng et al., 2006)。更多学者倾向于冈底斯南缘在早侏罗世时期应处于雅鲁藏布江洋板片俯冲的构造背景(邱检生等,2015董汉文等,2016Yu et al., 2018)。

目前对冈底斯中部早侏罗世基性侵入岩体的研究相对较少,而花岗岩的地球化学数据解释常因其多解性难以限定和讨论成岩构造环境(张旗,2012),所以对该期弧岩浆事件成因及构造背景仍不清楚,该归属于北部的班公湖-怒江洋?还是南部的雅鲁藏布江洋?亦或中部的松多古特提斯洋?还存在不同的认识。本文以冈底斯中部松多地区新发现的早侏罗世变质辉长岩为例,通过详细的岩石地球化学、锆石U-Pb定年及锆石Hf同位素研究,探讨其岩石成因和构造背景,为冈底斯中部晚三叠世-早侏罗世构造-岩浆演化提供新的约束。

1 地质背景及样品描述

松多地区处于冈底斯中东部,沙莫勒-米拉山断裂带东段,位于南冈底斯地块与中冈底斯地块的交界部位(图 1ab)。研究区出露的地质单元主要包括:石炭-二叠系松多岩组、侏罗-白垩系火山-沉积地层、晚古生代蛇绿混杂岩、松多榴辉岩以及大面积分布的中-新生代岩浆岩(图 1c)。其中石炭-二叠系松多岩组主体是一套小有序而大无序的构造地层,以陆源碎屑岩为主,夹少量碳酸盐岩的沉积岩系,主要岩性包括变质石英砂岩、白云母石英片岩、绢云石英片岩、二云(长石)石英片岩及少量黑云石英片岩组成。晚古生代蛇绿混杂岩岩石端元齐全,主要包括变质橄榄岩、堆晶辉长岩、辉长岩、玄武岩、斜长花岗岩等(王斌等,2017)。榴辉岩带近东西向展布,岩石呈包体形式产于松多岩组之中,主要岩性包括金红石榴辉岩、石英榴辉岩、角闪榴辉岩和多硅白云母榴辉岩等。现有研究表明,松多榴辉岩原岩为MORB型岩石(陈松永等,2007),原岩可能形成于石炭-二叠纪(徐向珍等,2007),变质年龄可能为晚二叠世(Cheng et al., 2012, 2015)。区内断裂构造十分发育,构造形迹以近东西向为主,具有多期次、脆性逆冲断裂构造发育的特点。

本文研究的变质辉长岩分布于松多地区加色村以北,变质辉长岩呈岩脉状近南北向产出,虽然受后期片理化改造较强,但仍可见变质辉长岩与石炭-二叠系松多岩组变质砂岩的侵入接触界线(图 2a)。岩石样品为灰绿色、灰黑色,块状构造(图 2b)。岩相学观察发现,辉长岩已经发生了明显的角闪岩相-低绿片岩相变质作用改造,具有中粗粒变余辉长结构,主要矿物由角闪石(15%~20%)、斜长石(35%~40%)、绿泥石(25%~35%)、绿帘石(10%~15%)等矿物组成(图 2cd)。辉石蚀变较强烈,仅局部可见辉石假象,已基本蚀变为角闪石,角闪石也部分绿泥石化;斜长石,板柱状,普遍绿泥石化和钠黝帘石化;绿泥石,淡绿色(图 2c),定向排列,具有一级灰干涉色(图 2d);黝帘石呈阴影状集合体;绿帘石,粒状,短柱状,粒度在0.1~0.6mm之间,正高突起,具有姜黄色异常干涉色(图 2d)。岩石中还含有少量石英颗粒,副矿物主要为磁铁矿。

图 2 西藏松多地区早侏罗世变质辉长岩野外露头及显微特征 (a)变质辉长岩与变质砂岩野外侵入接触界线露头;(b)变质辉长岩片理化现象露头近景照片;变质辉长岩单偏光下(c)和正交偏光下(d)显微照片. Hbl-角闪石;Pl-斜长石;Chl-绿泥石;Ep-绿帘石;Q-石英;Mag-磁铁矿 Fig. 2 Occurrence and micrographs of metamorphic gabbro in the Sumdo area, northern Tibet
2 分析方法 2.1 全岩地球化学测试

样品的无污染碎样在河北廊坊市(宇能)宇恒矿岩技术服务有限公司完成,将岩石样品用无污染碎样机精碎至小于5mm左右,选用无污染玛瑙球磨机磨至200目后送往地球化学分析实验室备用。全岩地球化学样品的主量元素、微量元素分析均在中国地质大学(北京)科学研究院实验中心完成。在超净实验室内使用电子天平上称取50mg样品,采用两酸(HNO3+HF)高压反应釜(Bomb)溶样方法进行样品的化学预处理,分析仪器为美国安捷伦公司生产Agilent7500a型等离子质谱仪,分析过程中使用美国地质调查局标样AGV2、W2、BHOV和中国地质测试中心岩石标样R1、R3进行分析质量检查和监控。

2.2 锆石U-Pb测年

锆石分选在在河北廊坊市(宇能)宇恒矿岩技术服务有限公司完成,采用常规粉碎、浮选和电磁选方法进行分选,最后在双目显微镜下挑纯。锆石阴极荧光图像分析(CL)在中国地质科学院电镜室的阴极荧光分析系统上完成。在中国地质大学(北京)地学实验中心进行了透射光和反射光显微照相及锆石LA-ICP-MS微区原位U-Pb同位素分析。分析时,测定部位的选取依据显微图像和阴极发光图像特点,选择无或者少包裹体的部位,尽可能避开裂纹部位;避免测定位置跨越不同世代的晶体区域,这些照片为数据解释提供依据。所使用的ICP-MS为美国Agilent科技有限公司的7500A ICP-MS,激光剥蚀系统为美国New Wave贸易有限公司,UP193SS型,深紫外(DUV)193nm、ArF准分子激光剥蚀系统。激光束斑直径为36μm,剥蚀采样时间为45s。实验中采用剥蚀物质载气为高纯度He,流速0.7L/min,用标准锆石91500为外标进行同位素比值校正(Wiedenbeck et al., 1995),用标准锆石TEM(417Ma)做监控盲样,元素含量以NIST612为外标进行标定,29Si作为内标,NIST612和NIST614做监控盲样。采用澳大利亚Glitter 4.4软件对同位素比值及元素含量计算进行处理(ver4.0,Macquarie University)(Ludwig,2003),采用ISOPLOT 3.0软件绘制U-Pb谐和图和计算加权平均值(Yuan et al., 2004)。

2.3 锆石Hf同位素测试

锆石Hf同位素分析在北京科萃测试技术有限公司实验室利用激光剥蚀多接收器电感耦合等离子体质谱仪完成。激光进样系统为NWR213nm固体激光器,分析系统为多接收等离子体质谱仪(NEPTUNE plus)。检测环境:温度18~22℃,相对湿度 < 65%。仪器运行条件及分析方法详见文献(Hu et al., 2012)。利用NWR213nm固体激光器对锆石进行剥蚀,激光剥蚀的斑束直径为40μm,能量密度为10~11J/cm2,频率为10Hz,激光剥蚀物质以高纯He为载气送入Neptune Plus (MC-ICPMS),接收器配置与溶液进样方式相同。分析数据处理(包括对样品和空白信号的选择、同位素质量分馏校正)采用软件ICP MSD ata Cal(Liu et al., 2010)完成。

3 分析结果 3.1 全岩地球化学特征

在主量元素分析表(表 1)中,变质辉长岩SiO2含量变化较大,具有较高的MgO、Fe2O3T、CaO含量。其中SiO2含量为45.56%~53.69%;MgO含量为6.04%~11.75%,Mg#为66~77,接近原生玄武质岩石(Mg#=70,Dupuy and Dostal, 1984);Fe2O3T含量为6.87%~9.53%;CaO含量为7.24%~11.76%。TiO2含量较低(0.77%~1.69%)。K2O+Na2O为4.52%~2.71%,K2O/Na2O < 0.5;样品的Al2O3含量变化较大,为12.68%~18.64%;低P2O5(0.08%~0.33%)含量。

表 1 松多变质辉长岩地球化学组成(主量元素:wt%;稀土和微量元素:×10-6) Table 1 Whole-rock geochemical compositions of the Sumdo metamorphic gabbro(Major elements: wt%; Teace elements: ×10-6)

在稀土微量元素分析表(表 1)中,岩石的稀土总量较低,∑REE介于32×10-6~91×10-6,均值为73×10-6,明显低于OIB(ΣREE=199×10-6),而略高于E-MORB(ΣREE=49×10-6)稀土元素总量(Sun and McDonough, 1989)。在球粒陨石标准化稀土元素配分图上,样品表现右倾平滑曲线(图 3a),显示轻稀土元素相对富集,重稀土元素元素比较平坦的分配模式。轻重稀土元素分异显著,(La/Yb)N=2.05~4.80,主要是轻稀土分馏((La/Sm)N=1.49~2.47),重稀土分馏较弱或基本无分馏((Gd/Yb)N=1.20~1.67)。其中4件样品出现了较为明显的Eu负异常(δEu=0.57~0.73),其他样品无Eu负异常。

图 3 松多地区变质辉长岩球粒陨石标准化稀土元素配分图(a)及原始地幔标准化微量元素蛛网图(b) 标准化值、OIB、N-MORB、E-MORB的数据据Sun and McDonough (1989). OIB-洋岛玄武岩;E-MORB-富集型洋脊玄武岩;N-MORB正常型洋脊玄武岩 Fig. 3 Chondrite-normalized REE patterns (a) and primitive mantle-normalized trace element spider diagrams (b) for the metamorphic gabbro in the Sumdo area

松多变质辉长岩富集大离子亲石元素(LILE)(表 1),如Cs、Rb、Ba等,相容元素Cr、Ni含量变化较大,分别介于209×10-6~746×10-6和80×10-6~234×10-6;Co含量变化范围较小,介于29×10-6~56×10-6之间,总体与正常洋中脊玄武岩(N-MORB)含量相当(Co=47×10-6Hofmann,1988);在剔除活动元素组成的原始地幔标准化微量元素蛛网图中(图 3b),富集高场强元素Nb、Ta等,并出现显著的Pb正异常。大部分样品的Nb(5.82×10-6~24.1×10-6)、Ta(0.34×10-6~1.47×10-6)、Zr(47.9×10-6~153×10-6)、Hf(1.23×10-6~3.56×10-6)等元素含量总体上略高于E-MORB而明显低于OIB含量(Sun and McDonough, 1989)。总体上,样品的稀土配分模式和微量元素配分曲线与N-MORB存在显著差异,而介于E-MORB和OIB之间。

3.2 锆石U-Pb定年

3件变质辉长岩(St1、St2、St3)的LA-ICP-MS锆石U-Pb定年结果见表 2,相对应的锆石测点稀土元素分析结果见表 3。锆石CL图像显示,锆石可以分为两类。一类锆石呈自形到半自形长柱状-板柱状,颗粒长径集中在150~200μm,长轴/短轴比值为2︰1~4︰1,所测锆石具有明显的生长韵律环带(图 4b, d, f);这些锆石测点的Th和U含量变化范围较大,分别为71×10-6~2950×10-6和114×10-6~1646×10-6,相应Th/U范围为0.54~1.79,均大于0.4。锆石的稀土元素球粒陨石标准化配分模式以HREE相对富集和具有明显Ce正异常和Eu负异常为主要特征(图 4b, d, f),结合锆石存在生长韵律环带、较高的Th/U比值,表明这些锆石为典型的基性岩浆锆石成因(Hoskin and Black, 2000)。另一类锆石数量较少,颗粒形态也偏小,主要呈现浑圆状,具有核幔边结构特点,具有岩浆生长环带,结合相对较高的Th/U值(0.35~1.45),表明这部分锆石属于具有岩浆成因的捕获锆石。

表 2 松多变质辉长岩LA-ICP-MS锆石U-Pb分析结果 Table 2 LA-ICP-MS U-Pb analytical results for zircons from the Sumdo metamorphic gabbro

表 3 松多变质辉长岩锆石稀土元素分析结果(×10-6) Table 3 The trace element compositions (×10-6) of the zircons from the Sumdo metamorphic gabbro

图 4 松多变质辉长岩锆石206Pb/238U-207Pb/235U年龄关系图、锆石球粒陨石标准化稀土元素配分模式和典型锆石阴极发光照片 实线圈为锆石U-Pb测试点位,虚线圈为锆石Lu-Hf同位素测试点位 Fig. 4 Zircon chondrite-normalized REE patterns, diagrams of 206Pb/238U-207Pb/235U relationship and representative CL images of the analyzed zircons in Sumdo metamorphic gabbro

样品(St1)锆石的16个测点中有13个测点的206Pb/238U-207Pb/235U谐和性较好,测点数据在谐和曲线上成群分布,计算的206Pb/238U加权平均年龄为200±1Ma(图 4a),代表变质辉长岩的原岩时代,其中捕获锆石的206Pb/238U年龄为2706±24Ma、2566±23Ma、588±6Ma。样品(St2)锆石的9个测点中存在3颗捕获锆石,206Pb/238U年龄为3434±30Ma、556±6Ma、206±2Ma;其余测点计算的206Pb/238U加权平均年龄为195±6Ma(图 4c)。样品(St3)锆石的7个测点中也存在3颗捕获锆石,206Pb/238U年龄为847±9Ma、250±6Ma、211±4Ma;其余测点计算的206Pb/238U加权平均年龄为202±2Ma(图 4e)。3件变质辉长岩样品的锆石U-Pb定年结果在误差范围内基本一致,暗示辉长岩的岩浆结晶年龄为早侏罗世早期。

3.3 锆石Hf同位素组成

样品锆石Hf同位素分析点位置和锆石U-Pb测点相对应,分析结果见表 4。锆石176Lu/177Hf介于0.0005~0.0014,说明这些锆石在其形成之后,没有明显的放射性成因Hf的积累,所以本文测定的176Lu/177Hf比值应代表了锆石形成时体系Hf同位素组成(吴福元等,2007)。每个测点的εHf(t)值和模式年龄根据其U-Pb年龄计算。εHf(t)介于-1.9~3.9之间,单阶段模式年龄tDM1介于774~927Ma之间,二阶段模式年龄tDM2介于988~1361Ma之间。fLu/Hf为-0.96~-0.98,明显小于镁铁质(下地壳)及硅铝质地壳(上地壳)的fLu/Hf值(分别为-0.34和-0.72),所以其二阶段模式年龄更能反映源区物质从亏损地幔被抽取的时间或源区物质在地壳的平均存留年龄(Amelin et al., 1999)。

表 4 松多变质辉长岩中锆石Lu-Hf同位素组成 Table 4 The Lu-Hf isotope compositions of the zircons from the Sumdo metamorphic gabbro
4 讨论 4.1 变质作用对全岩组成的影响

本文变质辉长岩经历了后期角闪石化、绿泥石化等角闪岩相-绿片岩相的变质作用改造,所以在利用元素讨论岩石成因和源区性质之前要首先解读元素的活动性(Xu et al., 2001)。在变质过程中,稀土元素与高场强元素(Th、Nb、Ta、Zr、Hf、Y)是相对稳定的(Rolliso,1993;Kerrich et al., 1999)。利用元素与最不活动性元素Zr在双变量图中的相关性也可以准确判别岩石遭受变质改造的条件(Polat et al., 2009, 2012),相关系数R < 0.75被认为是活动性元素,是判别元素活动性简便有效的方法(Polat and Hofmann, 2003刘平华等,2012)。松多变质辉长岩的元素Ti、REE、Th、Nd、Nb、Ta与Zr的相关性如图 5所示,TiO2与Zr具有非常好的相关性,相关系数R=0.91;稀土元素La、Gd、Yb、Sm、Nd与Zr均具有较高的相关性(R分别为0.91、0.91、0.78、0.95、0.97);高场强元素如Nb、Th、Ta与Zr也具有非常好的相关性(R分别为0.93、0.91、0.93)。其他微量元素Hf、U、Ni、Cr与Zr亦具有较高的相关性。上述元素与Zr的相关性分析表明,在变质过程中,这些元素基本保持稳定,无明显的迁移变化,代表变质辉长岩原岩的元素含量和特征;Cs、Ba和Pb与Zr也具有一定的相关性(R分别为0.70、0.71、0.64),遭受了轻微影响而发生迁移,基本也保持和继承了原岩特征,所以上述元素可以用来恢复和探讨变质辉长岩的岩浆系列、岩石成因和源区性质。与稀土元素和高场强元素相比较,而大离子亲石元素常常是活动元素(Rolliso,1993;Kerrich et al., 1999),大离子亲石元素(Sr、Rb)、主量元素(K、Na、Ca、Al)与Zr也没有明显的相关性(图 6),表明这些元素可能已经在变质作用过程中发生了迁移,所以这些强活动性的元素不宜用来讨论岩浆系列和岩石成因。

图 5 松多变质辉长岩Zr与主、微量元素(Ti、REE、Th、Nd、Nb、Ta)相关图 Fig. 5 Zr against selected major and trace elements (Ti, REE, Th, Nd, Nb, Ta) variation diagrams for the Sumdo metamorphic gabbro

图 6 松多变质辉长岩Zr与主、微量元素(Sr、Rb、Cs、K、Na、Ca、Al、Pb)相关图 Fig. 6 Zr against selected major and trace elements (Sr, Rb, Cs, K, Na, Ca, Al, Pb) variation diagrams for the Sumdo metamorphic gabbro
4.2 源区性质和岩石成因

由于变质作用的影响,K、Na等活动性元素已经发生了明显的迁移,因此本文选择没有发生显著变化的惰性组分对变质辉长岩进行原岩岩石系列的判别。在Nb/Y-Zr/TiO2(图 7a)与Nb/Y-SiO2(图 7b)分类图解中,松多变质辉长岩样品点落在亚碱性玄武岩和碱性玄武岩的交界区域。

图 7 松多变质辉长岩Nb/Y-Zr/TiO2 (a,据Pearce, 1996)与Nb/Y-SiO2 (b,据Winchester and Floyd, 1997)图解 Fig. 7 Nb/Y vs. Zr/TiO2 (a, after Pearce, 1996) and Nb/Y vs. SiO2 (b, after Winchester and Floyd, 1997) diagrams for the Sumdo metamorphic gabbro

松多变质辉长岩以高MgO(6.04%~11.75%)、高Fe2O3T(6.87%~9.53%)、高Mg#值(66~77)和富Na2O(2.31%~4.15%)为特征,并具有较高的Sc(28.0×10-6~46.7×10-6)、Co(28.6×10-6~55.6×10-6)、Ni(80.3×10-6~234.0×10-6)含量,表明其具有地幔的源区性质(Rapp and Watson, 1995)。那么,是起源于岩石圈地幔还是软流圈地幔源区(Sklyarov et al., 2003Zhao and Zhou, 2007)?松多变质辉长岩的稀土总量较低,具有轻稀土元素相对富集,重稀土元素元素比较平坦的分配模式,同时富集大离子亲石元素(LILE,如Cs、Rb、Ba),富集高场强元素Nb、Ta等,这与正常洋中脊玄武岩(N-MORB)具有亏损LREE、LILE特征的软流圈地幔存在显著差异(Sun and McDonough, 1989)。样品的Nb/La(1.14~1.61)>1,Nb/U比值平均为53.9,与OIB和E-MORB的值相一致(Hofmann,1988)。蛛网图呈现介于OIB和E-MORB之间的平滑配分曲线(图 3b)。在2Nb-Zr/4-Y图解(图 8a)中,样品主要投在AII和富集型洋中脊玄武岩区域内;在Nb/Yb-Th/Yb图解(图 8b)和Nb/Yb-TiO2/Yb图解(图 8c)中,样品点投在地幔演化趋势中E-MORB和OIB过度区域。但考虑到辉长岩形成过程中很可能存在堆晶作用,其原始岩浆的成分应该更加接近OIB。上述特征暗示松多变质辉长岩原岩应起源于亏损岩石圈地幔的部分熔融。

图 8 松多变质辉长岩的地球化学判别图解 (a) 2Nb-Zr/4-Y图解(Meschede,1986);(b) Nb/Yb-Th/Yb图解(Pearce,2008);(c) Nb/Yb-TiO2/Yb图解(Pearce,2008).AI-板内碱性玄武岩;AⅡ-板内碱性玄武岩和板内拉斑玄武岩;B-富集型洋中脊玄武岩;C-板内拉斑玄武岩和火山弧玄武岩;D-正常型洋中脊玄武岩;OIB-洋岛玄武岩;N-MORB-正常型洋中脊玄武岩;E-MORB-富集型洋中脊玄武岩 Fig. 8 Geochemical discrimination diagrams for the Sumdo metamorphic gabbro

锆石Hf同位素组成因其很少受到后期地质作用的影响,近年来成为探讨岩浆源区属性的重要约束(Griffin et al., 2000)。松多变质辉长岩中岩浆锆石的εHf(t)介于-1.9~3.9之间,显示出富集的锆石Hf同位素组成,进一步揭示其岩浆起源于亏损岩石圈地幔的部分熔融。松多变质辉长岩出现显著的Pb正异常,并在锆石中出现较多较老具有岩浆成因的捕获锆石,表明岩浆在上升侵位过程中可能发生了地壳物质的混染。由于在分离结晶和部分熔融过程中,总分配系数相同或相近的元素比值基本不受影响,所以常常用来检验地壳混染是否存在及其程度。通常认为高的(La/Sm)N值(>4.5)指示了地壳物质的混染(张永明等,2019),松多变质辉长岩中(La/Sm)N值较低,并具有较高Mg#值(66~77),而Nb、Ta也并未出现亏损,反而相对富集,这都表明岩浆演化过程中同化混染作用是有限的。

对于被下地壳物质改造的地幔源区部分熔融形成的岩石存在两种主要机制,即拆沉陆壳物质对岩浆源区的直接改造和岩浆上升侵位过程中被陆壳物质的同化混染(杨浩田等,2018)。这两种岩浆演化过程的本质区别在于,前者对岩浆源区的改造是地幔物质与地壳物质发生熔融并最终形成均一的元素-同位素地球化学体系,而后者是幔源岩浆上升侵位到地壳深度过程中同化混染了地壳物质,并导致岩浆的元素-同位素组成发生明显的变化(Yang et al., 2008杨浩田等,2018)。而松多变质辉长岩中相容元素含量变化较大可能与岩浆结晶分异作用有关,如Cr和Ni含量表现出较宽的变化范围(209×10-6~746×10-6和80×10-6~234×10-6),表明岩浆在演化过程中存在不同程度的镁铁质矿物分离结晶。在MgO与Cr和Ni的Harker图解中,Cr和Ni与MgO具有强烈的正相关性(图 9),说明该区变质辉长岩原岩形成过程中存在明显的橄榄石、尖晶石和单斜辉石的分离结晶。个别样品中还出现了Eu的负异常,暗示在其源区也可能存在斜长石的残留。综上所述,改造松多变质辉长岩原岩亏损岩石圈地幔源区的壳源熔体可能主要来自于拆沉的古老下地壳物质的部分熔融,下地壳物质对幔源源区的改造应主要发生在岩浆源区。

图 9 松多变质辉长岩MgO-Cr (a)与MgO-Ni (b)图解 Fig. 9 MgO vs. Cr (a) and MgO vs. Ni (b) variation diagrams for the Sumdo metamorphic gabbro

在地幔橄榄岩的熔融过程中稀土元素属于中等不相容元素(Johnson,1994),一般其浓度不会受地幔亏损和流体混入的明显影响(Pearce and Peate, 1995Münker,2000)。如Sm、Yb和Ce、Y在尖晶石二辉橄榄岩体系中分配系数相近(McKenzie and O’Nions,1991Deniel,1998),当石榴石二辉橄榄岩发生低程度熔融时,常常出现HREE显著的分馏,由于溶体中Yb和Y的含量由于受控于残留的石榴石(Johnson,1994),分配系数分别高于Sm和Ce,所以具有较低的Sm/Yb比值(1.26~1.91)和Ce/Y比值,以及相对于尖晶石二辉橄榄岩部分熔融较陡的演化趋势(图 10Deniel,1998);而由尖晶石二辉橄榄岩部分熔融形成的岩浆往往具有HREE相对平坦的稀土配分特点(蓝江波等,2007)。因此松多变质辉长岩HREE相对平坦的稀土配分曲线(图 3b)、低的(La/Yb)N值(2.06~4.81)、(La/Sm)N值(1.50~2.48)、(Gd/Yb)N值(1.21~1.67)和较低的(Ce/Y)N值(小于2)均表明其岩浆源区可能与角闪石的残留有关,而不存在石榴石。岩石中低(Tb/Yb)N比值(1.16~1.48,),均值1.31,也表明其主要来源于含尖晶石二辉橄榄岩的部分熔融(Wang et al., 2002)。Deniel(1998)提出了以不活动高场强元素判别玄武质岩浆源区矿物组合的有效图解(图 10),本文样品均位于原始尖晶石相二辉橄榄岩和亏损尖晶石相二辉橄榄岩熔融源区之间,暗示其源区深度可能在70km左右(Lambert and Wyllie, 1968)。

图 10 松多变质辉长岩Zr/Nb-Ce/Y图解 Fig. 10 Variation of Zr/Nb vs. Ce/Y for the Sumdo metamorphic gabbro
4.3 大地构造背景

冈底斯中东部早侏罗世存在较大规模火山-岩浆作用,是中生代冈底斯陆缘岩浆弧形成的重要时期(Chu et al., 2006Zhu et al., 2008Ji et al., 2012Kang et al., 2014董汉文等,2016Yu et al., 2018)。除了本文报道的早侏罗世变质辉长岩外,在中冈底斯东部已经识别出多处早侏罗世的中酸性火山-岩浆作用(图 1a),但是对具有弧岩浆作用地球化学特征的成因及大地构造背景还存在不同的认识,该归属于“哪个洋”的弧?是北部的班公湖-怒江洋?还是南部的雅鲁藏布江洋?有学者推测叶巴组火山岩可能形成于班公湖-怒江洋南向俯冲消减的结果(Geng et al., 2006耿全如等,2006);班公湖-怒江大洋板片俯冲到南冈底斯以下,在板片回转作用下形成弧后拉张,从而使地幔和下地壳部分熔融,导致叶巴组及相关早侏罗世岩浆作用的出现(Zhu et al., 2008)。如果冈底斯南缘的侏罗系桑日群和叶巴组形成于班公湖-怒江洋南向俯冲而导致的弧后盆地环境,那么冈底斯南缘应该存在一些古老地壳,但是冈底斯南缘却以新生地壳为主(Zhu et al., 2011);叶巴组分布范围现今距离班公湖-怒江缝合带至少在300km以上,而且在早白垩世冈底斯南部还曾经出现大规模的地壳缩短(England and Houseman, 1986Murphy et al., 1997),所以在早侏罗世叶巴组产出的位置距离班公湖-怒江洋的位置更远,受到班公湖-怒江洋南向俯冲的影响更是不太可能(黄丰等,2015)。董彦辉等(2006)认为雅鲁藏布江洋开始俯冲的时代应在早中侏罗世或更早;雅鲁藏布江洋可能在早三叠世之前已经打开(王玉净等,2002Chen et al., 2019)。在叶巴组火山岩南部还发育着大规模早侏罗世桑日群火山岩,最新研究表明叶巴组火山岩表现为典型的大陆边缘弧特征,而桑日群具有洋内弧火山岩,可能分别代表了雅鲁藏布江洋北向俯冲形成的陆缘弧和洋内弧(黄丰等,2015)。

近年来,松多古特提斯洋的发现和研究,为认识上述争论提供了新的思路。松多地区已经陆续发现了与古大洋演化相关的中二叠世蛇绿岩(王斌等,2017)、中二叠世洋岛(Wang et al., 2019),证实了松多古特提斯洋的存在。大量有关榴辉岩及其变质作用的深入研究(杨经绥等,2006Cheng et al., 2015李鹏等,2017),也进一步揭示在松多地区存在以松多古特提斯洋北向俯冲闭合产物组成的缝合带,其闭合时代为晚三叠世末期(李化启等,2008)。李化启等(2008)通过对松多地区贡布帕拉韧性剪切带的构造解析研究,获得含白云母绿片岩40Ar/39Ar得年龄为230~220Ma,提出拉萨地块松多岩组经历了印支期造山事件的认识,并进一步指出这是松多古特提斯洋深俯冲之后南、北拉萨地体碰撞的结果。李化启等(2011)对松多地区的区域构造变形的深入研究,又获得糜棱质绿片岩和白云母石英片岩的白云母40Ar/39Ar同位素年龄为~220Ma,进一步证实该期地壳浅部构造变形事件的存在;与此同时,获得榴辉岩及退变榴辉岩中白云母和角闪石40Ar/39Ar同位素年龄为240~220Ma,代表榴辉岩的折返时代,不仅进一步表明冈底斯中部于220~240Ma经历过碰撞造山事件,同时还揭示出松多古特提斯洋在220Ma以后已经消失,南、北拉萨地体已经发生碰撞。碰撞后必然是地壳厚度减薄,处于伸展的构造背景(张旗,2012)。有研究表明,在陆陆碰撞之后10~15Ma俯冲消减到大陆之下的残余洋壳发生板片断离而逐渐下沉(van de Zedde and Wortel, 2001)。如果松多地区在~220Ma发生了陆陆碰撞,那么持续俯冲的松多古特提斯洋壳便可能会在~200Ma左右发生板片断离和板片回返而引起弧后伸展,引发软流圈地幔上涌并发生减压熔融,在同时有俯冲板片提供流体或者俯冲组分的情况下,引发岩石圈地幔的部分熔融,从而有可能形成本文研究的变质辉长岩的岩浆源区。正是岩石圈地幔的部分熔融,由此产生的幔源岩浆能提供足够的热引起上覆古老的浅部地壳物质重熔,形成大规模的中酸性岩浆岩(张旗,2012)。已有大量研究表明,冈底斯中部确实存在晚三叠世-早侏罗世的花岗岩带,包括尼木县城以南的尼木花岗岩(205Ma、202Ma)(张宏飞等,2007)、南木林县折无巨斑状花岗岩(217Ma)(李才等,2003)、当雄南宁中地区白云母-二长花岗岩(190Ma)(刘琦胜等,2006)、松多地区的花岗闪长岩和二长花岗岩(215~190Ma)(和钟铧等,2006Yu et al., 2018),均为冈底斯南部晚三叠世-早侏罗世(217~190Ma)中酸性岩浆事件,该期恰好晚于松多地区记录的碰撞造山事件和榴辉岩折返退变时代约10~15Ma。综合本文变质辉长岩岩石成因和冈底斯地区晚三叠世-早侏罗世中酸性岩体的展布特征,这些花岗岩可能是和冈底斯中部曾发生过的印支期碰撞造山事件有关的晚碰撞或后碰撞花岗岩(李化启等,2011)。

5 初步结论

综合以上的讨论,本文得出如下几点初步认识:

(1) 松多变质辉长岩具高MgO和Mg#值以及Cr、Ni含量,富集大离子亲石元素和稀土元素,Nb、Ta等高场强元素富集,锆石εHf(t)介于-1.9~3.9之间,显示其岩浆起源于在下地壳改造的亏损岩石圈地幔的部分熔融。

(2) 松多地区变质辉长岩原岩形成于早侏罗世早期(~200Ma),反应了松多古特提斯洋闭合后,晚三叠世-早侏罗世冈底斯中部有可能存在松多古特提斯洋壳板片断离事件,引发岩石圈地幔的部分熔融。

致谢      样品锆石LA-ICP-MS原位微区U-Pb同位素测定得到中国地质大学(北京)苏犁教授的热情帮助,在此致以诚挚谢意。

恰逢肖序常先生90华诞,谨以此文表达对他的衷心祝愿,祝福先生永远健康。

参考文献
Amelin Y, Lee DC, Halliday AN and Pidqeon RT. 1999. Nature of the Earth's earliest crust from hafnium isotopes in single detrital zircons. Nature, 399(6733): 252-255 DOI:10.1038/20426
Cao DD, Cheng H, Zhang LM and Wang K. 2017. Post-peak metamorphic evolution of the Sumdo eclogite from the Lhasa terrane of Southeast Tibet. Journal of Asian Earth Sciences, 143: 156-170 DOI:10.1016/j.jseaes.2017.04.020
Chen DS, Luo H, Wang XH, Xu B and Matsuoka A. 2019. Late Anisian radiolarian assemblages from the Yarlung-Tsangpo Suture Zone in the Jinlu area, Zedong, southern Tibet:Implications for the evolution of Neotethys. Island Arc, 28(4): e12302
Chen SY, Yang JS, Luo LQ, Li ZL, Xu XZ, Li TF, Ren YF and Li HQ. 2007. MORB-type eclogites in the Lhasa block, Tibet, China:Petrochemical evidence. Geological Bulletin of China, 26(10): 1327-1339 (in Chinese with English abstract)
Cheng H, Zhang C, Vervoort JD, Lu HH, Wang C and Cao DD. 2012. Zircon U-Pb and garnet Lu-Hf geochronology of eclogites from the Lhasa Block, Tibet. Lithos, 155: 341-359 DOI:10.1016/j.lithos.2012.09.011
Cheng H, Liu YM, Vervoort JD and Lu HH. 2015. Combined U-Pb, Lu-Hf, Sm-Nd and Ar-Ar multichronometric dating on the Bailang eclogite constrains the closure timing of the Paleo-Tethys Ocean in the Lhasa terrane, Tibet. Gondwana Research, 28(4): 1482-1499 DOI:10.1016/j.gr.2014.09.017
Chu MF, Chung SL, Song B, Liu DY, O'Reilly SY, Pearson NJ, Ji JQ and Wen DJ. 2006. Zircon U-Pb and Hf isotope constraints on the Mesozoic tectonics and crustal evolution of southern Tibet. Geology, 34(9): 745-748 DOI:10.1130/G22725.1
Deniel C. 1998. Geochemical and isotopic (Sr, Nd, Pb) evidence for plume-lithosphere interactions in the genesis of Grande Comore magmas (Indian Ocean). Chemical Geology, 144(3-4): 281-303 DOI:10.1016/S0009-2541(97)00139-3
Dong HW, Xu ZQ, Li Y and Liu Z. 2016. Genesis and geological significance of Early Jurassic gabbro in Medog area, the eastern part of the Gangdese magmatic belt. Acta Petrologica Sinica, 32(12): 3624-3634 (in Chinese with English abstract)
Dong YH, Xu JF, Zeng QG, Wang Q, Mao GZ and Li J. 2006. Is there a Neo-Tethys' subduction record earlier than arc volcanic rocks in the Sangri Group?. Acta Petrologica Sinica, 22(3): 661-668 (in Chinese with English abstract)
Dupuy C and Dostal J. 1984. Trace element geochemistry of some continental tholeiites. Earth and Planetary Science Letters, 67(1): 61-69 DOI:10.1016/0012-821X(84)90038-4
England P and Houseman G. 1986. Finite strain calculations of continental deformation:2. Comparison with the India-Asia collision zone. Journal of Geophysical Research:Solid Earth, 91(B3): 3664-3676
Geng QR, Pan GT, Zheng LL, Chen ZL, Fisher RD, Sun ZM, Ou CS, Dong H, Wang XW, Li S, Lou XY and Fu H. 2006. The eastern Himalayan syntaxis:Major tectonic domains, ophiolitic mélanges and geologic evolution. Journal of Asian Earth Sciences, 27(3): 265-285 DOI:10.1016/j.jseaes.2005.03.009
Geng QR, Pan GT, Wang LQ, Zhu DC and Liao ZL. 2006. Isotopic geochronology of the volcanic rocks from the Yeba Formation in the Gangdise zone, Xizang. Sedimentary Geology and Tethyan Geology, 26(1): 1-7 (in Chinese with English abstract)
Griffin WL, Pearson NJ, Belousova E, Jackson SE, Van Achterbergh E, O'Reilly SY and Shee SR. 2000. The Hf isotope composition of cratonic mantle:LAM-MC-ICPMS analysis of zircon megacrysts in kimberlites. Geochimica et Cosmochimica Acta, 64(1): 133-147 DOI:10.1016/S0016-7037(99)00343-9
He ZH, Yang DM, Zheng CQ and Wang TW. 2006. Isotopic dating of the mamba granitoid in the gangdise tectonic belt and its constraint on the subduction time of the Neotethys. Geological Review, 52(1): 100-106 (in Chinese with English abstract)
Hofmann AW. 1988. Chemical differentiation of the earth:The relationship between mantle, continental crust, and oceanic crust. Earth and Planetary Science Letters, 90(3): 297-314 DOI:10.1016/0012-821X(88)90132-X
Hoskin PWO and Black LP. 2000. Metamorphic zircon formation by solid-state recrystallization of protolith igneous zircon. Journal of Metamorphic Geology, 18(4): 423-439
Hu ZC, Liu YS, Gao S, Liu WG, Zhang W, Tong XR, Lin L, Zong KQ, Li M, Chen HH, Zhou L and Yang L. 2012. Improved in situ Hf isotope ratio analysis of zircon using newly designed X skimmer cone and Jet sample cone in combination with the addition of nitrogen by laser ablation multiple collector ICP-MS. Journal of Analytical Atomic Spectrometry, 27(9): 1391-1399 DOI:10.1039/c2ja30078h
Huang F, Xu JF, Chen JL, Kang ZQ and Dong YH. 2015. Early Jurassic volcanic rocks from the Yeba Formation and Sangri Group:Products of continental marginal arc and intra-oceanic arc during the subduction of Neo-Tethys Ocean?. Acta Petrologica Sinica, 31(7): 2089-2100 (in Chinese with English abstract)
Ji WQ, Wu FY, Chung SL and Liu CZ. 2012. Identification of Early Carboniferous granitoids from southern Tibet and implications for terrane assembly related to the Paleo-Tethyan evolution. Journal of Geology, 120(5): 531-541 DOI:10.1086/666742
Johnson KTM. 1994. Experimental cpx/and garnet/melt partitioning of REE and other trace elements at high pressures:Petrogenetic implications. Mineralogical Magazine, 58A(1): 454-455 DOI:10.1180/minmag.1994.58A.1.236
Kang ZQ, Xu JF, Wilde SA, Feng ZH, Chen JL, Wang BD, Fu WC and Pan HB. 2014. Geochronology and geochemistry of the Sangri Group volcanic rocks, southern Lhasa terrane:Implications for the early subduction history of the Neo-Tethys and Gangdese magmatic arc. Lithos, 200-201: 157-168 DOI:10.1016/j.lithos.2014.04.019
Kerrich R, Polat A, Wyman D and Hollings P. 1999. Trace element systematics of Mg-, to Fe-tholeiitic basalt suites of the Superior Province:Implications for Archean mantle reservoirs and greenstone belt genesis. Lithos, 46(1): 163-187
Lambert IB and Wyllie PJ. 1968. Stability of hornblende and a model for the low velocity zone. Nature, 219(5160): 1240-1241 DOI:10.1038/2191240a0
Lan JB, Xu YG, Yang QJ and Huang XL. 2007. ~40Ma OIB-type mafic magmatism in the Gaoligong belt:Results of break-off between subducting Tethyan slab and Indian plate?. Acta Petrologica Sinica, 23(6): 1334-1346 (in Chinese with English abstract)
Li C, Wang TW, Li HM and Zeng QG. 2003. Discovery of Indosinian megaporphyritic granodiorite in the Gangdise area:Evidence for the existence of Paleo-Gangdise. Geological Bulletin of China, 22(5): 364-366 (in Chinese with English abstract)
Li HQ, Cai ZH, Chen SY, Tang ZM and Yang M. 2008. The Indosinian orogenesis occurred in Lhasa terrain and the evidence from muscovite 40Ar-39Ar geochronology. Acta Petrologica Sinica, 24(7): 1595-1604 (in Chinese with English abstract)
Li HQ, Xu ZQ, Yang JS, Tang ZM and Yang M. 2011. Sys-collisional exhumation of Sumdo eclogite in the Lhasa Terrane, Tibet:Evidences from structural deformation and 40Ar/39Ar geochronology. Earth Science Frontiers, 18(3): 66-78 (in Chinese with English abstract)
Li P, Zhang C, Liu XY, Shen TT, Qiu T and Yang JS. 2017. The metamorphic processes of the Xindaduo eclogite in Tibet and its constrain on the evolutionary of the Paleo-Tethys subduction zone. Acta Petrologica Sinica, 33(12): 3753-3765 (in Chinese with English abstract)
Liu PH, Liu FL, Wang F, Liu JH, Yang H and Shi JR. 2012. Geochemical characteristics and genesis of the high-pressure mafic granulite in the Jiaobei high-grade metamorphic basement, eastern Shandong, China. Acta Petrologica Sinica, 28(9): 2705-2720 (in Chinese with English abstract)
Liu QS, Jiang W, Jian P, Ye PS, Wu ZH and Hu DG. 2006. Zircon SHRIMP U-Pb age and petrochemical and geochemical features of Mesozoic muscovite monzonitic granite at Ningzhong, Tibet. Acta Petrologica Sinica, 22(3): 643-652 (in Chinese with English abstract)
Liu Y, Liu HF, Theye T and Massonne HJ. 2009. Evidence for oceanic subduction at the NE Gondwana margin during Permo-Triassic times. Terra Nova, 21(3): 195-202 DOI:10.1111/j.1365-3121.2009.00874.x
Liu YS, Gao S, Hu ZC, Gao CG, Zong KQ and Wang DB. 2010. Continental and oceanic crust recycling-induced melt-peridotite interactions in the Trans-North China Orogen:U-Pb dating, Hf isotopes and trace elements in zircons from mantle xenoliths. Journal of Petrology, 51(1-2): 537-571 DOI:10.1093/petrology/egp082
Ludwig KR. 2003. ISOPLOT 3.0:A Geochronological Toolkit for Microsoft Excel. Berkeley Geochronology Center Special Publication
McKenzie D and O'Nions RK. 1991. Partial melt distributions from inversion of rare earth element concentrations. Journal of Petrology, 32(5): 1021-1091 DOI:10.1093/petrology/32.5.1021
Meschede M. 1986. A method of discriminating between different types of mid-ocean ridge basalts and continental tholeiites with the Nb-Zr-Y diagram. Chemical Geology, 56(3-4): 207-218 DOI:10.1016/0009-2541(86)90004-5
Münker C. 2000. The isotope and trace element budget of the Cambrian Devil River arc system, New Zealand:Identification of four source components. Journal of Petrology, 41(6): 759-788 DOI:10.1093/petrology/41.6.759
Murphy MA, Yin A, Harrison TM, Dürr BB, Chen Z, Ryerson FJ, Kidd WSF, Wang X and Zhou X. 1997. Did the Indo-Asian collision alone create the Tibetan Plateau?. Geology, 25(8): 719-722 DOI:10.1130/0091-7613(1997)025<0719:DTIACA>2.3.CO;2
Pan GT, Mo XX, Hou ZQ, Zhu DC, Wang LQ, Li GM, Zhao ZD, Geng QR and Liao ZL. 2006. Spatial-temporal framework of the Gangdese Orogenic Belt and its evolution. Acta Petrologica Sinica, 22(3): 521-533 (in Chinese with English abstract)
Pearce JA and Peate DW. 1995. Tectonic implications of the composition of volcanic ARC magmas. Annual Review of Earth and Planetary Sciences, 23: 251-285 DOI:10.1146/annurev.ea.23.050195.001343
Pearce JA. 1996. A user's guide to basalt discrimination diagrams. In: Wyman DA (ed.). Trace Element Geochemistry of Volcanic Rocks: Applications for Massive Sulphide Exploration. Geological Association of Canada, Short Course Notes, 12: 79-113
Pearce JA. 2008. Geochemical fingerprinting of oceanic basalts with applications to ophiolite classification and the search for Archean oceanic crust. Lithos, 100(1-4): 14-48 DOI:10.1016/j.lithos.2007.06.016
Polat A and Hofmann AW. 2003. Alteration and geochemical patterns in the 3.7~3.8Ga Isua greenstone belt, West Greenland. Precambrian Research, 126(3-4): 197-218 DOI:10.1016/S0301-9268(03)00095-0
Polat A, Appel PWU, Fryer B, Windley B, Frei R, Samson IM and Huang H. 2009. Trace element systematics of the Neoarchean Fiskenaesset anorthosite complex and associated meta-volcanic rocks, SW Greenland: Evidence for a magmatic arc origin. Precambrian Research, 175(1-4): 87-115
Polat A, Fryer BJ, Samson IM, Weisener C, Appel PWU, Frei R and Windley BF. 2012. Geochemistry of ultramafic rocks and hornblendite veins in the Fiskenaesset layered anorthosite complex, SW Greenland:Evidence for hydrous upper mantle in the Archean. Precambrian Research, 214-215: 124-153 DOI:10.1016/j.precamres.2011.11.013
Qiu JS, Wang RQ, Zhao JL and Yu SB. 2015. Petrogenesis of the Early Jurassic gabbro-granite complex in the middle segment of the Gangdese belt and its implications for tectonic evolution of Neo-Tethys:A case study of the Dongga pluton in Xi'gaze. Acta Petrologica Sinica, 31(12): 3569-3580 (in Chinese with English abstract)
Rapp RP and Watson EB. 1995. Dehydration melting of metabasalt at 8~32kbar:Implications for continental growth and crust-mantle recycling. Journal of Petrology, 36(4): 891-931 DOI:10.1093/petrology/36.4.891
Rollison HR. 1993. Using Geochemical Data:Evaluation, Presentation, Interpretation. London: Longman Group, 1-352
Sklyarov EV, Gladkochub DP, Mazukabzov AM, Menshagin YV, Watanabe T and Pisarevsky SA. 2003. Neoproterozoic mafic dike swarms of the Sharyzhalgai metamorphic massif, southern Siberian Craton. Precambrian Research, 122(1-4): 359-376 DOI:10.1016/S0301-9268(02)00219-X
Sun SS and McDonough WF. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. In: Saunders AD and Norry MJ (eds.). Magmatism in the Ocean Basins. Geological Society, London, Special Publications, 42(1): 313-345
van de Zedde DMA and Wortel MJR. 2001. Shallow slab detachment as a transient source of heat at midlithospheric depths. Tectonics, 20(6): 868-882 DOI:10.1029/2001TC900018
Wang B, Xie CM, Li C, Liu YM and Yu YP. 2017. The discovery of Wenmulang ophiolite in Songduo area of the Tibetan Plateau and its geological significance. Geological Bulletin of China, 36(11): 2076-2081 (in Chinese with English abstract)
Wang B, Xie CM, Fan JJ, Wang M, Yu YP, Dong YC and Hao YJ. 2019. Genesis and tectonic setting of Middle Permian OIB-type mafic rocks in the Sumdo area, southern Lhasa terrane. Lithos, 324-325: 429-438 DOI:10.1016/j.lithos.2018.11.015
Wang K, Plant T, Walker JD and Smith EI. 2002. A mantle melting profile across the Basin and Range, SW USA. Journal of Geophysical Research:Solid Earth, 107(B1): ECV 5-1-ECV 5-21 DOI:10.1029/2001JB000209
Wang YJ, Yang Q, Matsuoka A, Kobayashi K, Nagahashi T and Zeng QG. 2002. Triassic radiolarians from the Yarlung Zangbo Suture Zone in the Jinlu area, Zetang County, southern Tibet. Acta Micropalaeontologica Sinica, 19(3): 215-227 (in Chinese with English abstract)
Weller OM, St-Onge MR, Rayner N, Waters DJ, Searle MP and Palin RM. 2016. U-Pb zircon geochronology and phase equilibria modelling of a mafic eclogite from the Sumdo complex of south-east Tibet:Insights into prograde zircon growth and the assembly of the Tibetan plateau. Lithos, 262: 729-741 DOI:10.1016/j.lithos.2016.06.005
Wiedenbeck M, Allé P, Corfu F, Griffin WL, Meier M, Oberli F, von Quadt A, Roddick JC and Spiegel W. 1995. Three natural zircon standards for U-Th-Pb, Lu-Hf, trace element and REE analyses. Geostandards and Geoanalytical Research, 19(1): 1-23 DOI:10.1111/j.1751-908X.1995.tb00147.x
Winchester J and Floyd PA. 1977. Geochemical discrimination of different magma series and their differentiation products using immobile elements. Chemical Geology, 20: 325-343 DOI:10.1016/0009-2541(77)90057-2
Wu FY, Li XH, Zheng YF and Gao S. 2007. Lu-Hf isotopic systematics and their applications in petrology. Acta Petrologica Sinica, 23(2): 185-220 (in Chinese with English abstract)
Xiao XC and Wang J. 1998. A brief review of tectonic evolution and uplift of the Qinghai-Tibet. Geological Review, 44(4): 372-381 (in Chinese with English abstract)
Xu XZ, Yang JS, Li TF, Chen SY, Ren YF, Li ZL and Shi YR. 2007. SHRIMP U-Pb ages and inclusions of zircons from the Sumdo eclogite in the Lhasa block, Tibet, China. Geological Bulletin of China, 26(10): 1340-1355 (in Chinese with English abstract)
Xu YG, Menzies MA, Thirlwall MF and Xie GH. 2001. Exotic lithosphere mantle beneath the western Yangtze craton:Petrogenetic links to Tibet using highly magnesian ultrapotassic rocks. Geology, 29(9): 863-866 DOI:10.1130/0091-7613(2001)029<0863:ELMBTW>2.0.CO;2
Xu ZQ, Yang JS, Li WC, Li HQ, Cai ZH, Yan Z and Ma CQ. 2013. Paleo-Tethys system and accretionary orogen in the Tibet Plateau. Acta Petrologica Sinica, 29(6): 1847-1860 (in Chinese with English abstract)
Yang CH, Xu WL, Yang DB, Wang W, Wang WD and Liu JM. 2008. Petrogenesis of Shangyu gabbro-diorites in western Shandong:Geochronological and geochemical evidence. Science in China (Series D), 51(4): 481-492 DOI:10.1007/s11430-008-0029-0
Yang HT, Yang DB, Shi JP, Xu WL, Wang F and Chen Y. 2018. Nature of the Early Cretaceous lithospheric mantle in western Shangdong:Constraints from geochronology, geochemistry and Sr-Nd-Pb-Hf isotopic data of Dakunlun gabbros and diabases. Acta Petrologica Sinica, 34(11): 3327-3340 (in Chinese with English abstract)
Yang JS, Xu ZQ, Geng QR, Li ZL, Xu XZ, Li TF, Ren YF, Li HQ, Cai ZH, Liang FH and Chen SY. 2006. A possible new HP/UHP (?) metamorphic belt in China:Discovery of eclogite in the Lasha Terrane, Tibet. Acta Geologica Sinica, 80(12): 1787-1792 (in Chinese with English abstract)
Yu YP, Xie CM, Fan JJ, Wang M, Dong YC, Wang B and Hao YJ. 2018. Zircon U-Pb geochronology and geochemistry of Early Jurassic granodiorites in Sumdo area, Tibet:Constraints on petrogenesis and the evolution of the Neo-Tethyan Ocean. Lithos, 320-321: 134-143 DOI:10.1016/j.lithos.2018.09.006
Yuan HL, Gao S, Liu XM, Li HM, Günther D and Wu FY. 2004. Accurate U-Pb age and trace element determinations of zircon by Laser Ablation-Inductively Coupled Plasma-Mass Spectrometry. Geostandards and Geoanlytical Research, 28(3): 353-370 DOI:10.1111/j.1751-908X.2004.tb00755.x
Zhang HF, Xu WC, Guo JQ, Zong KQ, Cai HM and Yuan HL. 2007. Zircon U-Pb and Hf isotopic composition of deformed granite in the southern margin of the Gangdese belt, Tibet:Evidence for Early Jurassic subduction of Neo-Tethyan oceanic slab. Acta Petrologica Sinica, 23(6): 1347-1353 (in Chinese with English abstract)
Zhang KJ and Tang XC. 2009. Eclogites in the interior of the Tibetan Plateau and their geodynamic implications. Chinese Science Bulletin, 54(15): 2556-2567
Zhang Q. 2012. Collision and granite:Collision is a tectonic event, not a tectonic environment. Acta Petrologica et Mineralogica, 31(5): 745-749 (in Chinese with English abstract)
Zhang YM, Pei XZ, Li ZC, Li RB, Liu CJ, Pei L, Chen YX and Wang M. 2019. Zircon U-Pb geochronogy, geochemistry and its geological implication of the Early Indosinian basic complex in the Qinghai Nanshan tectonic belt. Earth Science, 44(7): 2461-2481 (in Chinese with English abstract)
Zhao JH and Zhou MF. 2007. Geochemistry of Neoproterozoic mafic intrusions in the Panzhihua district (Sichuan Province, SW China):Implications for subduction-related metasomatism in the upper mantle. Precambrian Research, 152(1-2): 27-47 DOI:10.1016/j.precamres.2006.09.002
Zhu DC, Pan GT, Chung SL, Liao ZL, Wang LQ and Li GM. 2008. SHRIMP zircon age and geochemical constraints on the origin of Lower Jurassic volcanic rocks from the Yeba Formation, southern Gangdese, South Tibet. International Geology Review, 50(5): 442-471 DOI:10.2747/0020-6814.50.5.442
Zhu DC, Pan GT, Wang LQ, Mo XX, Zhao ZD, Zhou CY, Liao ZL, Dong GC and Yuan SH. 2008. Spatial-temporal distribution and tectonic setting of Jurassic magmatism in the Gangdise belt, Tibet, China. Geological Bulletin of China, 27(4): 458-468 (in Chinese with English abstract)
Zhu DC, Zhao ZD, Niu YL, Mo XX, Chung SL, Hou ZQ, Wang Q and Wu FY. 2011. The Lhasa Terrane:Record of a microcontinent and its histories of drift and growth. Earth and Planetary Science Letters, 301(1-2): 241-255 DOI:10.1016/j.epsl.2010.11.005
陈松永, 杨经绥, 罗立强, 李兆丽, 徐向珍, 李天福, 任玉峰, 李化启. 2007. 西藏拉萨地块MORB型榴辉岩的岩石地球化学特征. 地质通报, 26(10): 1327-1339. DOI:10.3969/j.issn.1671-2552.2007.10.011
董汉文, 许志琴, 李源, 刘钊. 2016. 冈底斯带东段墨脱地区早侏罗世辉长岩的成因及其构造意义. 岩石学报, 32(12): 3624-3634.
董彦辉, 许继峰, 曾庆高, 王强, 毛国政, 李杰. 2006. 存在比桑日群弧火山岩更早的新特提斯洋俯冲记录么?. 岩石学报, 22(3): 661-668.
耿全如, 潘桂棠, 王立全, 朱弟成, 廖忠礼. 2006. 西藏冈底斯带叶巴组火山岩同位素地质年代. 沉积与特提斯地质, 26(1): 1-7. DOI:10.3969/j.issn.1009-3850.2006.01.001
和钟铧, 杨德明, 郑常青, 王天武. 2006. 冈底斯带门巴花岗岩同位素测年及其对新特提斯洋俯冲时代的约束. 地质论评, 52(1): 100-106. DOI:10.3321/j.issn:0371-5736.2006.01.013
黄丰, 许继峰, 陈建林, 康志强, 董彦辉. 2015. 早侏罗世叶巴组与桑日群火山岩:特提斯洋俯冲过程中的陆缘弧与洋内弧?. 岩石学报, 31(7): 2089-2100.
蓝江波, 徐义刚, 杨启军, 黄小龙. 2007. 滇西高黎贡带~40Ma OIB型基性岩浆活动:消减特提斯洋片与印度板块断离的产物?. 岩石学报, 23(6): 1334-1346. DOI:10.3969/j.issn.1000-0569.2007.06.010
李才, 王天武, 李惠民, 曾庆高. 2003. 冈底斯地区发现印支期巨斑花岗闪长岩——古冈底斯造山的存在证据. 地质通报, 22(5): 364-366. DOI:10.3969/j.issn.1671-2552.2003.05.011
李化启, 蔡志慧, 陈松永, 唐哲民, 杨梅. 2008. 拉萨地体中的印支造山事件及年代学证据. 岩石学报, 24(7): 1595-1604.
李化启, 许志琴, 杨经绥, 唐哲民, 杨梅. 2011. 拉萨地体内松多榴辉岩的同碰撞折返:来自构造变形和40Ar-39Ar年代学的证据. 地学前缘, 18(3): 66-78.
李鹏, 张聪, 刘晓瑜, 申婷婷, 邱添, 杨经绥. 2017. 西藏新达多地区榴辉岩的变质过程研究及其对古特提斯俯冲带演化的限定意义. 岩石学报, 33(12): 3753-3765.
刘平华, 刘福来, 王舫, 刘建辉, 杨红, 施建荣. 2012. 胶北高级变质基底中高压基性麻粒岩的地球化学特征及其成因. 岩石学报, 28(9): 2705-2720.
刘琦胜, 江万, 简平, 叶培盛, 吴珍汉, 胡道功. 2006. 宁中白云母二长花岗岩SHRIMP锆石U-Pb年龄及岩石地球化学特征. 岩石学报, 22(3): 643-652.
潘桂棠, 莫宣学, 侯增谦, 朱弟成, 王立全, 李光明, 赵志丹, 耿全如, 廖忠礼. 2006. 冈底斯造山带的时空结构及演化. 岩石学报, 22(3): 521-533.
邱检生, 王睿强, 赵姣龙, 喻思斌. 2015. 冈底斯中段早侏罗世辉长岩-花岗岩杂岩体成因及其对新特提斯构造演化的启示:以日喀则东嘎岩体为例. 岩石学报, 31(12): 3569-3580.
王斌, 解超明, 李才, 刘一鸣, 于云鹏. 2017. 青藏高原松多地区温木朗蛇绿岩的发现及其地质意义. 地质通报, 36(11): 2076-2081. DOI:10.3969/j.issn.1671-2552.2017.11.017
王玉净, 杨群, 松冈笃, 小林健太, 武井雅彦, 曾庆高. 2002. 藏南泽当雅鲁藏布缝合带中的三叠纪放射虫. 微体古生物学报, 19(3): 215-227. DOI:10.3969/j.issn.1000-0674.2002.03.001
吴福元, 李献华, 郑永飞, 高山. 2007. Lu-Hf同位素体系及其岩石学应用. 岩石学报, 23(2): 185-220.
肖序常, 王军. 1998. 青藏高原构造演化及隆升的简要评述. 地质论评, 44(4): 372-381. DOI:10.3321/j.issn:0371-5736.1998.04.006
徐向珍, 杨经绥, 李天福, 陈松永, 任玉峰, 李兆丽, 石玉若. 2007. 青藏高原拉萨地块松多榴辉岩的锆石SHRIMP U-Pb年龄及锆石中的包裹体. 地质通报, 26(10): 1340-1355. DOI:10.3969/j.issn.1671-2552.2007.10.012
许志琴, 杨经绥, 李文昌, 李化启, 蔡志慧, 闫臻, 马昌前. 2013. 青藏高原中的古特提斯体制与增生造山作用. 岩石学报, 29(6): 1847-1860.
杨浩田, 杨德彬, 师江朋, 许文良, 王枫, 谌瑛. 2018. 鲁西早白垩世岩石圈地幔的属性:大昆仑辉长岩和辉绿岩年代学、岩石地球化学和Sr-Nd-Pb-Hf同位素制约. 岩石学报, 34(11): 3327-3340.
杨经绥, 许志琴, 耿全如, 李兆丽, 徐向珍, 李天福, 任玉峰, 李化启, 蔡志慧, 梁凤华, 陈松永. 2006. 中国境内可能存在一条新的高压/超高压(?)变质带——青藏高原拉萨地体中发现榴辉岩带. 地质学报, 80(12): 1787-1792. DOI:10.3321/j.issn:0001-5717.2006.12.001
张宏飞, 徐旺春, 郭建秋, 宗克清, 蔡宏明, 袁洪林. 2007. 冈底斯南缘变形花岗岩锆石U-Pb年龄和Hf同位素组成:新特提斯洋早侏罗世俯冲作用的证据. 岩石学报, 23(6): 1347-1353. DOI:10.3969/j.issn.1000-0569.2007.06.011
张旗. 2012. 碰撞与花岗岩——碰撞是构造事件, 不是构造环境. 岩石矿物学杂志, 31(5): 745-749. DOI:10.3969/j.issn.1000-6524.2012.05.012
张永明, 裴先治, 李佐臣, 李瑞保, 刘成军, 裴磊, 陈有炘, 王盟. 2019. 青海南山构造带印支早期基性杂岩体年代学、地球化学特征及地质意义. 地球科学, 44(7): 2461-2481.
朱弟成, 潘桂棠, 王立全, 莫宣学, 赵志丹, 周长勇, 廖忠礼, 董国臣, 袁四化. 2008. 西藏冈底斯带侏罗纪岩浆作用的时空分布及构造环境. 地质通报, 27(4): 458-468. DOI:10.3969/j.issn.1671-2552.2008.04.003