岩石学报  2020, Vol. 36 Issue (2): 426-442, doi: 10.18654/1000-0569/2020.02.06   PDF    
西藏尼玛县西拉木那勒地区晚白垩世埃达克岩岩石成因及构造指示
董宇超1, 王明1,2, 解超明1,2, 于云鹏1, 郝宇杰2     
1. 吉林大学地球科学学院, 长春 130061;
2. 自然资源部东北亚矿产资源评价重点实验室, 长春 130061
摘要: 青藏高原北拉萨板块上发育大规模的晚白垩世岩浆作用,长期以来,关于其岩石成因及构造属性一直存在较大争论。本次研究对产出于北拉萨板块尼玛县西约150km的拉木那勒晚白垩世岩体开展了锆石U-Pb年代学、地球化学等方面的工作。拉木那勒岩体主要为花岗闪长岩和二长花岗岩,具有较高的SiO2(68.12%~72.21%)、Al2O3(12.28%~17.02%)含量和较高的Sr/Y比值(26.7~145)、La/Yb比值(24.1~28.4),低含量的MgO(0.88%~1.86%)、Y(4.58×10-6~10.03×10-6)和Yb(0.33×10-6~0.73×10-6),具有埃达克岩的典型特征。本文共对4件年代学样品进行了锆石U-Pb定年,其中2件样品获得了88.4±0.5Ma和88.6±0.4Ma的谐和年龄,另外2件样品获得了89Ma和95Ma的下交点年龄。结合区域地质资料,我们倾向于认为拉木那勒岩体形成于班公湖-怒江洋闭合之后,是陆-陆碰撞加厚下地壳部分熔融的产物。此外,地壳增厚可能是北拉萨板块晚白垩世岩浆岩的一个重要的形成机制。
关键词: 青藏高原    晚白垩世    埃达克岩    地壳增厚    
Genesis and tectonic indication of the Late Cretaceous adakite rocks in the Lamunale area, Nima County, Tibet
DONG YuChao1, WANG Ming1,2, XIE ChaoMing1,2, YU YunPeng1, HAO YuJie2     
1. College of Earth Sciences, Jilin University, Changchun 130061, China;
2. MNR Key Laboratory of Mineral Resources Evaluation in Northeast Asia, Changchun 130061, China
Abstract: Large-scale Late Cretaceous magmatism occurs in the North Lhasa Block on the Tibet Plateau. For a long time, there have been large debates about its petrogenesis and tectonic properties. This time we carried out zircon U-Pb chronology, geochemistry and other aspects of the Lamunale Late Cretaceous intrusion located about 150km west of Nima County in the North Lhasa Block. The lithology of the Lamunale rock is composed of granodiorite and monzonitic granite with high content of SiO2 (68.12%~72.21%), Al2O3 (12.28%~17.02%), Sr/Y ratios (26.7~145), La/Yb ratios (24.1~28.4), low content of MgO (0.88%~1.86%), Y (4.58×10-6~10.03×10-6) and Yb (0.33×10-6~0.73×10-6), which is a typical characteristics of adakite. A total four samples were analyzed to determine their radiometric ages. Two samples achieved a concordant U-Pb age of 88.4±0.5Ma and 88.6±0.4Ma, whereas the other two samples provide the lower intersection ages of 89Ma and 95Ma. Combining regional geological data, we tend to think that the Lamunale rock formed after the closure of the Bangong-Nujiang Ocean Basin was a product of partial melting of the lower crust due to land-land collision thickening. Crustal thickening may be an important mechanism for the formation of the Late Cretaceous magmatic rocks in North Lhasa Block.
Key words: Tibet Plateau    Late Cretaceous    Adakite rock    Crust thickening    

拉萨板块夹持于羌塘板块与喜马拉雅板块之间,同时又被洛巴堆-米拉山断裂和狮泉河-纳木错蛇绿混杂岩带分为北、中、南三部分,构造背景十分复杂。拉萨板块上分布着大量从中生代到新生代的岩浆岩(朱弟成等, 2006, 2008; 马国林等, 2010; 吴珍汉等, 2014; 张硕等, 2014)。前人对这些岩浆岩的成因问题上存在较大的争论。在中、南拉萨板块上的中新生代岩浆岩普遍被认为与雅鲁藏布江新特提斯洋的北向俯冲有关,如平板俯冲、板片回返、板片断离、洋脊俯冲等(丁慧霞等, 2015; 王珍珍等, 2017);也有部分学者认为藏东大规模出露的白垩纪岩浆岩可能与澳大利亚同时期的岩浆岩具有亲缘性(朱弟成等, 2013);而学者们对于北拉萨板块上广泛分布的白垩纪岩浆岩的成因问题仍然具有较大的争论,主要争议在于这些岩浆岩究竟是与新特提斯洋板片北向俯冲有关还是与班公湖-怒江洋(后简称班-怒洋)南向俯冲作用有关,亦或是该期岩浆岩是狮泉河-永珠-嘉黎蛇绿混杂岩带闭合的产物(Wen et al., 2008)。本文对北拉萨板块尼玛县西约150km的拉木那勒地区新识别出的晚白垩世埃达克岩开展研究,通过锆石U-Pb同位素定年、全岩地球化学研究,并结合区域上已发表的数据,旨在对北拉萨板块的晚白垩世岩浆岩的成因及构造意义提供约束。

1 区域地质背景及样品情况

青藏高原位于特提斯构造域的东段,是世界上平均海拔最高的高原之一,具有全球罕见的地质形迹,是地质学研究的理想场所。研究表明,青藏高原是古生代以来冈瓦纳大陆北侧裂解的块体之间洋盆俯冲增生及相互聚敛碰撞形成的巨型拼贴体,后经历新生代(70~40Ma)印度大陆与欧亚大陆碰撞作用的产物,是研究陆-陆碰撞、造山隆升、海陆转换等地质作用的理想场所(Powell and Conaghan, 1973; Molnar and Tapponnier, 1975; 潘桂棠等, 1997; Yin and Harrison, 2000; 徐志琴等, 2013)。拉萨板块位于青藏高原的中南部,南北两侧分别与喜马拉雅板块和羌塘板块相邻,并分别以印度-雅鲁藏布江缝合带和班公湖-怒江缝合带(后简称班-怒带)为界线;拉萨板块从北到南分为北拉萨板块(NL)、中拉萨板块(CL)和南拉萨板块(SL),以洛巴堆-米拉山断裂(LMF)和狮泉河-纳木错蛇绿混杂岩带(SNMZ)为界线(图 1a)。拉萨板块之上广泛发育中、新生代的岩浆岩,南拉萨板块地壳大部分为新生地壳且仅有少部分的前寒武纪基底被保留下来(Mo et al., 2008; Ji et al., 2009; Zhu et al., 2013),以分布广泛的三叠-白垩纪冈底斯岩基为特征(Chu et al., 2006; Mo et al., 2008; Wen et al., 2008)。中拉萨板块主要以大范围出露晚侏罗-早白垩世的海陆交互相沉积地层、则弄群火山岩和早白垩世的侵入岩为特征(王珍珍等, 2017)(图 1a)。

图 1 北拉萨板块拉木那勒地区地质简图及采样位置 (a)拉萨板块中-新生代岩浆岩分布简图(据Zhu et al., 2009);(b)研究区地质简图. BNSZ-班公湖-怒江缝合带;SNMZ-狮泉河-纳木错蛇绿混杂岩带;LMF-洛巴堆-米拉山断裂带;IYSSZ-印度-雅鲁藏布缝合带 Fig. 1 Simplified geological map of Lamunale in North Lhasa Block and the location of sampling (a) simplified Mesozoic-Cenozoic magmatic activity distribution of Lhasa Terrane (after Zhu et al., 2009); (b) simplified geological map of research area. BNSZ-Bangong-Nujiang suture zone; SNMZ-Shiquanhe-Nam Tso Mélange Zone; LMF-Luobadui-Milashan Fault; IYSSZ-Indus-Yarlung Zangbo Suture Zone

研究区位于北拉萨板块尼玛县西约150km的拉木那勒地区(图 1a, b),区域内出露的地层较为简单,主要由下白垩统朗山组生物碎屑灰岩、含生物碎屑泥晶灰岩和上白垩统马莫勒组砂砾岩组成。区内断裂构造发育,在研究区中部发育两条近东西延伸的逆断层和一条近南北向延伸的断层,分割了晚白垩世的侵入体和朗山组二段生物碎屑灰岩地层。拉木那勒岩体位于研究区的中部,整体出露面积约3km2(图 1b),与朗山组灰岩和复理石片岩呈侵入接触关系,在接触边界可见矽卡岩化和角岩化现象(图 2a),岩石受后期改造作用明显,风化破碎严重,部分露头可见基岩出露(图 2b-d)。

图 2 拉木那勒埃达克岩及竟柱山组砾岩野外照片 (a)花岗闪长岩(N17T1)与朗山组灰岩的侵入接触界线;(b)二长花岗岩(N17T44)近景照片;(c)花岗闪长岩(N17T1)近景照片;竟柱山组砾岩远景照片(d)和近景照片(e) Fig. 2 The field photos of the Lamunale adakites and Jingzhushan conglomerates (a) intrusive boundary between the limestone of Langshan Formation and granodiorite (N17T1); (b) the close-up photo of sample (N17T44); (c) the close-up photo of sample (N17T1); the distant phoot (d) and the close up photo (e) of Jingzhushan Fortmation conglomerate

本次研究我们重点采集了拉木那勒岩体中的花岗闪长岩样品(N17T1、N17T8、N17T9)和二长花岗岩样品(N17T44)。花岗闪长岩较为新鲜,表面呈灰白色或灰黄色,花岗结构,块状构造(图 2c),矿物组成为石英(35%~45%)、钾长石(15%~25%)、斜长石(30%~40%)以及少量的角闪石(5%~10%)。石英多呈他形粒状(图 3a-c),粒度在0.1~0.2mm;斜长石多呈自形-半自形板状,粒度在0.4~0.6mm,大多数的斜长石发育典型的聚片双晶结构,部分样品(N17T8)的斜长石发育良好的环带结构(图 3b);角闪石具有一定的蚀变,自形程度较低,多呈他形板状(图 3a-c),粒度在0.3~0.5mm。二长花岗岩整体为灰白色带黄褐色调,岩石表面有一定的蚀变现象(图 2b),但岩石内部较为新鲜,矿物组成为石英(30%~40%)、钾长石(25%~35%)、斜长石(25%~35%)以及少量的角闪石(< 5%)。钾长石成自形-半自形板状或粒状形态,表面有轻微的高岭土化(图 3d),粒度在0.2~0.4mm;斜长石呈灰白色,聚片双晶不明显,粒度在0.3~0.5mm之间。野外共采集4件年龄样品(样品号及采样坐标分别为N17T1:N31°45′42″、E86°4′10″;N17T8:N31°45′12″、E86°4′35″;N17T9:N31°45′44″、E86°6′15″;N17T44:N31°44′48″、E86°6′48″)和18件地球化学样品。

图 3 拉木那勒埃达克岩镜下照片 Pl-斜长石;Hb-角闪石;Q-石英;Or-正长石;ZP-环带的斜长石 Fig. 3 Photomicrographs of the Lamunale adakites Pl-plagioclase; Hb-hornblende; Q-quartz; Or-orthoclase; ZP-zoned plagioclase
2 分析方法 2.1 锆石U-Pb测年

样品锆石的分选和收集在河北省地质调查研究院实验室完成,在双目镜下挑选干净透明、无裂隙、无包裹体的锆石颗粒至环氧树胶中,后进行打磨剖光,使锆石的中心位置暴露出来,锆石的阴极荧光图像(CL)在中国地质科学院完成,锆石的LA-ICP-MS原位分析在中国地质大学(北京)地学实验中心进行,ICP-MS采用美国Agilent科技公司的7500a型,剥蚀系统为美国New Wave公司的UP193SS型。激光束斑直径36μm,剥蚀时间为45s,载气为氦气,流速为0.7L/min,锆石91500和QH为外标样进行同位素比值校正,Si为内标计算,更详细的分析步骤参照Yuan et al. (2004),普通Pb元素的校正采用Anderson (2002)。最后使用GLITTER 4.4软件对同位素数据处理,用Isoplot/Ex(version 4.15)软件绘制U-Pb谐和图与计算年龄加权平均值(Ludwig, 2003)。

2.2 全岩地球化学分析

样品在河北省地质调查研究院实验室无污染碎至200目,之后在中国地质大学(北京)地学实验中心完成全岩主量及微量元素分析。烧失量(LOI)是将1g样品至于熔炉中以980℃加热4~6个小时,随后在干燥器中冷却至常温重新称重计算得来。全岩主量元素分析采用美国Agilent科技公司的7500a型等离子体质谱仪(ICP-MS)测试,测试精度采用国际岩石参考物质AGV-2(USGS)与GSR-3。元素分析仪器为PS-950等离子提光谱仪,实测数据与推荐值具有相对偏差,分析精度通过测量值和推荐值之间的相对差异来估算,对于大部分样品,分析精度优于5%。更详细的实验步骤见参考文献(Hu et al., 2013; Zhai et al., 2013)。

3 分析结果 3.1 锆石U-Pb同位素测年结果

本文对3件花岗闪长岩样品(N17T1、N17T8、N17T9)和1件二长花岗岩样品(N17T44)进行了锆石U-Pb定年,分析结果见表 1

表 1 拉木那勒埃达克岩锆石U-Pb年龄数据 Table 1 Zircon U-Pb isotopic data of the Lamunale adakites

样品N17T1和N17T9的锆石自形程度高,整体呈灰白色或浅灰色,表面干净,包体较少,在个别颗粒内部保存有古老的继承核,锆石形态多呈短柱状或长柱状,锆石颗粒较大(75~150μm)且发育良好的岩浆震荡环带(图 4a, c),锆石的Th/U比值在0.28~1.10之间,微量元素的稀土配分曲线表现出明显的Ce的正异常(Ce*=1.01~119)和Eu的负异常(Eu*=0.07~0.86),以及富集的重稀土等特征(图 5a),我们认为二者为典型的岩浆成因锆石(Barbey et al., 1995; Belousova et al., 2002; Hoskin, 2005)。N17T1共进行了20个锆石颗粒的U-Pb同位素年代学分析,所有分析点的年龄结果十分集中,在锆石年龄谐和曲线中表现出良好的谐和性,206Pb/238U年龄在88~89Ma之间,加权平均年龄为88.6±0.4Ma(MSWD=0.26)(图 4a);N17T9共进行20个测试点的分析,其中6个分析点(01, 03, 08, 10, 12, 17)偏离谐和线,另外14个测点表现出较好的谐和性,年龄结果分布非常集中,206Pb/238U年龄在88~89Ma,加权平均年龄为88.4±0.5Ma(MSWD=0.24)(图 4c)。

图 4 拉木那勒埃达克岩锆石U-Pb年龄谐和图及典型锆石阴极发光图像 Fig. 4 Zircon U-Pb concordia diagrams and representative CL images of Lamunale adakites

图 5 拉木那勒埃达克岩锆石球粒陨石标准化稀土元素配分曲线(a, 标准化值据Boynton, 1984)与U/Yb-Y锆石成因判别图解(b, 底图据Grimes et al., 2007) Fig. 5 Chondrite-normalized REE patterns (a, normalization values after Boynton, 1984) and U/Yb vs. Y genesis discrimination diagram (b, base map after Grimes et al., 2007) of the zircon from Lamunale adakites

样品N17T8和N17T44锆石形态类似,均呈半自形-自形的短柱状或菱形形态,锆石颗粒较N17T1和N17T9略小,长径粒度在60~100μm之间,长宽比在1:1~2:1之间。与另外2件测年样品不同的是,在N17T8和N17T44部分锆石中保留着较大的古老的核部(图 4b, d),表现出类似变质成因锆石的核边结构,但与变质锆石不同的是,锆石普遍自形程度高,且Th/U在0.16~1.42之间,锆石轻重稀土分异明显,不具备重稀土亏损的特征(图 5a),这都表明了二者为典型的岩浆成因锆石而非变质成因锆石(Hoskin and Schaltegger, 2003; Chen et al., 2010)。这2件样品在锆石U-Pb年龄谐和图上未能形成有效的谐和年龄,一方面由于古老继承核的影响,另一方面也可能是由于Pb丢失的原因导致。可以通过构筑不一致线与谐和线的交点来分析其年龄,就所研究的2件样品而言,由于它们可能形成于相对年轻的时代,因此通过下交点的年龄来判断岩体的形成时代是确定分析样品年龄的有效手段,N17T8的下交点年龄为95Ma,N17T44的下交点年龄为89Ma,所以,通过前文的分析,我们认为这2件样品的结晶年龄也应在晚白垩世时期。

3.2 全岩地球化学分析结果

18件地球化学样品整体较为新鲜,烧失量(LOI)在1.01%~2.06%之间,详细分析结果见表 2

表 2 拉木那勒埃达克岩主量(wt%)和微量(×10-6)元素成分 Table 2 Whole-rock major (wt%) and trace (×10-6) elements compositions of the Lamunale adakite

拉木那勒花岗闪长岩包括样品整体具有高含量的SiO2(68.56%~72.21%)、Al2O3(12.88%~15.00%)和Na2O(3.17%~4.67%),低含量的MgO(1.13%~1.86%)和K2O(1.61%~2.65 %)。在R1-R2岩石定名图解中,样品全部落入花岗闪长岩区域(图 6a)。微量元素方面,样品表现出高的Sr含量(247.2×10-6~522.0×10-6)、高的Sr/Y比值(26.7~83.8)以及低的Y含量(6.14×10-6~10.03×10-6)和Yb含量(0.46×10-6~0.73×10-6),在YbN-(La/Yb)N和Y-Sr/Y地球化学判别图解中,样品全部落入埃达克岩区域(图 7a, b)。另外,样品具有低含量的Co(4.10×10-6~8.91×10-6)和Th(1.94×10-6~5.06×10-6),所有样品都落入了钙碱性系列(图 6b);在原始地幔标准化微量元素蛛网图上,拉木那勒花岗闪长岩富集大离子亲石元素(LILES)K、Rb、Ba等,亏损高场强元素Nb、Ta、Ti等(图 8a);在球粒陨石标准化稀土配分曲线中,样品整体轻重稀土分异作用较弱((La/Yb)N=17.3~20.6),表现轻微的轻稀土(LREE)富集和重稀土(LREE)亏损以及弱的铕异常(Eu*=0.27~0.30)(图 8b)。

图 6 拉木那勒埃达克岩岩石类型和系列划分图解 (a) R1-R2分类图解(据LeBas et al., 1986);(b) K2O-SiO2图解(据Hastie et al., 2007).区域数据据余红霞等, 2011; 张硕等, 2014; Sun et al., 2015; 图 7-图 10数据来源同此图 Fig. 6 Classification and series diagrams of the Lamunale adakites (a) R1 vs. R2 diagram (after LeBas et al., 1986); (b) K2O vs. SiO2 (after Hastie et al., 2007). Regional data according to Yu et al., 2011; Zhang et al., 2014; Sun et al., 2015; also in Fig. 7-Fig. 10

图 7 埃达克岩(La/Yb)N-YbN (a)和Sr/Y-Y (b)判别图解(据Defant and Drummond, 1990) Fig. 7 (La/Yb)N vs. YbN (a) and Sr/Y vs. Y (b) discrimination diagrams of adakites (after Defant and Drummond, 1990)

图 8 拉木那勒埃达克岩原始地幔标准化微量元素蛛网图(a, 标准化值据Sun and McDonough, 1989)与球粒陨石标准化稀土元素配分曲线图(b, 标准化值据Boynton, 1984) Fig. 8 Primitive mantle-normalized trace-element spidergrams (a, normalization values after Sun and McDonough, 1989) and chondrite-normalized REE patterns (b, normalization values after Boynton, 1984) of the Lamunale adakite

拉木那勒二长花岗岩地球化学特征与花岗闪长岩类似,整体表现为高含量的SiO2(68.12%~70.67%)、Al2O3(15.16%~17.02%)、Na2O(6.41%~6.94%)以及极高的Sr(613.2×10-6~683.0×10-6)和Sr/Y比值(134~145),相对低的MgO(0.88%~1.00%)、K2O(0.87%~1.03 %)以及Y(4.58×10-6~4.77×10-6)和Yb(0.33×10-6~0.34×10-6)含量。在R1-R2岩石定名图解中,样品全部落入二长花岗岩区域(图 6a)。样品具有极高的Sr含量,高的Al2O3含量和Sr/Y比,低的Y、Yb等地球化学特征与典型的埃达克岩十分相符,并且样品在YbN-(La/Yb)N和Y-Sr/Y地球化学判别图解中也全部落入了埃达克岩的区域(图 7a, b)。与拉木那勒花岗闪长岩类似,二长花岗岩在球粒陨石标准化稀土配分曲线中表现较弱的轻重稀土的分异((La/Yb)N=17.6~18.6)以及轻微的铕异常(Eu*=0.35~0.38)(图 8b)。另外,在原始地幔标准化微量元素蛛网图上,拉木那勒二长花岗岩同样富集大离子亲石元素(LILES)K、Rb、Ba等,亏损高场强元素Nb、Ta、Ti等(图 8a)。

4 讨论 4.1 岩石类型

自20世纪70年代以来,地质学界掀起了对花岗质岩石的分类热潮,约有20余种对花岗岩的分类方案被提出(Chappell and White, 1974; Chappell, 1999),而“埃达克岩”的分类方案就是在这种大背景下由Defant and Drummond (1990)最早提出。“埃达克岩”起初被认为是一种由俯冲的年轻且热的洋壳在榴辉岩相条件下熔融形成的一类具特殊地球化学特征的中酸性火成岩(Defant and Drummond, 1990; 王强等, 2008)。这些岩石通常具有以下特征:SiO2>56%,Al2O3>15%,MgO < 3%,Sr/Y比值(>20),La/Yb比值(>20),亏损Y(< 18×10-6)和重稀土元素(HREE)(Yb < 1.9×10-6),高的Sr(>400×10-6或300×10-6),无或正Eu、Sr异常,贫高场强元素(HFSE)等(Martin, 1999; Clemens, 2003; Castillo, 2006, 2012; 王强等, 2008)。拉木那勒侵入体(花岗闪长岩和二长花岗岩)整体上具有高含量的SiO2(68.12%~72.21%)和Al2O3(12.28%~17.02%)和高的Sr/Y比值(26.7~145)、La/Yb比值(24.1~28.4)、低含量MgO(0.88%~1.86%)、Y(4.58×10-6~10.03×10-6)和Yb(0.33×10-6~0.73×10-6),地球化学组成上十分符合埃达克岩的特征,另外在YbN-(La/Yb)N和Y-Sr/Y判别图解中,所有样品都落入埃达克岩的区域(图 7a, b),以上特征都表明了拉木那勒花岗闪长岩和二长花岗岩属于埃达克岩。样品N17T8具有更低的Sr含量(247.2×10-6~274.6×10-6),镜下观察样品发现其发育有环带明显的斜长石斑晶,由于斜长石在结晶过程会吸收一定量的Sr,可能会导致样品Sr含量的降低,但与之相悖的是斜长石斑晶的大量结晶作用又会使岩体里的Eu出现负异常,而该样品并未表现出明显的Eu的负异常。因此,我们认为N17T8样品中Sr含量的偏低不是由斜长石结晶作用所致,而是源区的Sr含量偏低导致(Zhu et al., 2009)。

4.2 岩石成因

埃达克岩提出伊始被认为是形成于年轻的(≤25Ma)与俯冲大洋岩石圈有关的岛弧环境中的火山岩或侵入岩(Defant and Drummond, 1990)。其不仅具有某些特殊的地球化学特征,还被认定与特定的构造背景相关(洋壳的部分熔融)。随后,越来越多的研究证明,具有与埃达克岩相似地球化学特征的岩石不仅可以出现在岛弧环境,而且在很多陆陆碰撞造山带、板内伸展等构造背景中都发现了埃达克岩的存在(王强等, 2008)。因此,越来越多的地质学者将“埃达克岩”的岩石分类方案仅仅定义为具有某些特定地球化学特征的岩石,而与岩石形成的构造背景无关。最近的研究表明,埃达克岩的成因主要分为以下5种:(1)洋壳部分熔融有或无地幔楔组分的加入(Defant and Drummond, 1990; Rapp et al., 1999; Martin et al., 2005; Zhang et al., 2010; Eyuboglu et al., 2011);(2)在高压下原始玄武质岩浆的分离结晶与同化混染(Defant and Drummond, 1990; Castillo et al., 1999; 康志强等, 2009);(3)酸性与基性岩浆的混合作用(Guo et al., 2007; Streck et al., 2007);(4)增厚下地壳或拆沉下地壳高温条件下的部分熔融(Wang et al., 2005; 王强等, 2008; Lai and Qin, 2013);(5)俯冲陆壳的部分熔融(王强等, 2008; Lai and Qin, 2013)。

拉木那勒埃达克岩具有高的SiO2含量(68.12%~72.21%),在稀土配分曲线中无明显的铕异常(0.27~0.35),表明了其不是原始玄武质岩浆高压下持续分离结晶与同化混染的产物。由于玄武质岩浆不可能直接形成酸性岩浆岩,地幔岩浆必须经历显著的分离结晶作用才能形成花岗质岩石,因此,这个过程必定会在区域上形成大量的基性岩浆岩(Fan et al., 2016; Wu et al., 2016),然而,区域研究结果表明在研究区附近的基性岩浆岩十分匮乏。另外,在Rb/Nd-Rb、La/Sm-Sm、La-La/Yb和Th/Nb-Th图解中(图 9),也表明了拉木那勒侵入体并不是原始基性岩浆分离结晶的产物;拉木那勒侵入体中含有大量古老的继承锆石,研究表明,在俯冲碰撞的构造环境下,一般只有陆壳物质的熔融才会形成大量的继承锆石,而洋壳直接熔融形成的岩浆岩中继承锆石往往也是由于陆壳物质的熔融所致。此外,锆石的微量元素分析表明样品锆石具有相对高的U/Yb比值和低的Y含量,在U/Yb-Y图解中,研究锆石全部落入大陆成因锆石的区域,这也进一步说明了拉木那勒侵入体并不是洋壳直接熔融的产物(Grimes et al., 2007)(图 5b)。拉木那勒岩体普遍具有高的SiO2含量,而且在岩石中没有基性岩捕掳体的存在,表明其并非酸性岩浆与基性玄武质岩浆混合形成的产物,这也与区域上基性岩浆岩整体出露十分匮乏相一致。前人研究表明,由于在班-怒洋两侧块体普遍缺失140~130Ma的岩浆活动,班-怒洋在此期间可能已经闭合,北拉萨板块与南羌塘板块发生弧-弧软碰撞(Zhu et al., 2016),拉木那勒侵入体侵位年龄在晚白垩世时期,此时北拉萨板块与南羌塘板块碰撞已过去40~50Ma,即使存在陆壳俯冲过程,陆壳也应该早已发生断离或熔融入地幔楔组分,无法在晚白垩世由熔融作用形成岩浆岩。

图 9 拉木那勒埃达克岩岩石成因判别图解 (a) Rb/Nd-Rb(据Schiano et al., 2010);(b) La/Sm-La(据Wu et al., 2015);(c) La-La/Yb(据Sun et al., 2015);(d) Th/Nb-Th(据Sun et al., 2015) Fig. 9 Discrimination diagrams of petrogenesis of Lamunale adakites (a) Rb/Nd vs. Rb (after Schiano et al., 2010); (b) La/Sm vs. La (after Wu et al., 2015); (c) La vs. La/Yb (after Sun et al., 2015) (d) Th/Nb vs. Th (after Sun et al., 2015)

综上所述,我们认为拉木那勒埃达克岩可能是增厚或拆沉下地壳部分熔融的产物,在主量元素与微量元素的哈克图解中,研究样品表现弱的相关性并且全部落入增厚下地壳的区域(图 10)。一般来说,由拆沉地壳部分熔融形成的埃达克岩,在岩浆上涌过程中受地幔组分的作用,往往具有高的MgO、Cr和Ni含量(Smithies, 2000; Prouteau et al., 2001; Martin et al., 2005; 余红霞等, 2011),拉木那勒埃达克质岩具有较低的MgO(0.88%~1.86%)、Ni(4.89×10-6~24.58×10-6)和Cr(6.19×10-6~47.16×10-6)含量,表明其不太可能是拆沉下地壳部分熔融的产物,而是与那些增厚下地壳部分熔融形成的埃达克岩类似,在哈克图解中拉木那勒埃达克岩大部分落入增厚下地壳成因的区域(图 10)。另外,前人在藏北尼玛县尼则地区的岩体中发现了89~89Ma的埃达克岩,具有较高的La/Yb比值(12~26.2),被认为可能是下地壳增厚超过30km部分熔融的产物(Liu et al., 2019),相比之下,拉木那勒岩体具有更高的La/Yb值(24.1~28.7),因此,我们认为拉木那勒侵入体可能是地壳增厚超过40km部分熔融的产物(图 11)。

图 10 拉木那勒埃达克岩哈克图解 不同埃达克岩成因区域划分据Wang et al. (2005); Lai and Qin (2013) Fig. 10 Harker diagrams for the Lamunale adakites Fields representing different genesis of adakite according to Wang et al. (2005); Lai and Qin (2013)

图 11 北拉萨板块白垩世时期(120~88Ma)地球动力学演化示意图(a, 据Zhu et al., 2016) (a)早白垩世时期地球动力学模型;(b)晚白垩世时期地球动力学模型 Fig. 11 Schematic illustrations showing the geodynamic evolution of the North Lhasa Block during the Cretaceous (120~88Ma) (a, after Zhu et al., 2016) (a) geodynamic model during the Early Cretaceous; (b) geodynamic model during the Late Cretaceous
4.3 地球动力学意义

北拉萨板块分布着大量的晚白垩世岩浆岩(100~79Ma),整体呈近东西条带状分布,延伸超过1000km,如此大规模的岩浆活动很可能是区域上某些构造运动作用的结果,然而有关北拉萨板块晚白垩世地球动力学背景的研究一直没有定论,主要可以归纳为3点:(1)狮泉河-纳木错蛇绿混杂岩带闭合的产物;(2)班-怒洋闭合的产物;(3)雅鲁藏布江新特提斯洋北向平板俯冲的产物。

拉木那勒埃达克岩出露位置夹持于班-怒带与狮泉河-永珠-嘉黎蛇绿混杂岩带之间(图 1a),研究表明在晚白垩世时期,两个大洋均已闭合,拉萨板块和南羌塘板块发生碰撞(Chen et al., 2014, 2015; Wang et al., 2014),同时期竟柱山组复成分砾岩(96~88Ma)角度不整合覆盖在朗山组灰岩(113~96Ma)之上,被认为是海陆相转化快速造山的重要证据(Sun et al., 2015; 李华亮等, 2016)。根据报道,在狮泉河-永珠-嘉黎蛇绿混杂岩带出露大量的板内玄武岩(Chen et al., 2014; Sui et al., 2013)、A型花岗岩和双峰式火山岩,这表明其明显处于伸展的构造背景,狮泉河-纳木错蛇绿混杂岩带可能是班-怒洋的弧后洋盆,很难形成拉木那勒埃达克岩。Wen et al. (2008)认为冈底斯岩浆带上的晚白垩世埃达克岩是雅鲁藏布江洋北向平板俯冲的产物,而平板俯冲由于在俯冲过程中往往会形成隔热层,一般不会形成岩浆作用(Gutscher et al., 2000),而在拉萨板块南缘发现众多白垩世岛弧成因的岩浆岩(Ji et al., 2009),这与平板俯冲的模式相悖。另外,对中、北拉萨板块地壳厚度的研究表明,中冈底斯板块具有比北拉萨板块更大的地壳厚度,反映拉萨板块由北向南陆壳成熟度逐渐变高的特征,暗示岩浆活动南向俯冲的极性(朱弟成等, 2006)。

近些年的研究表明,班-怒洋可能存在双向俯冲过程(Zhu et al., 2016),在羌塘板块南缘和中北拉萨板块分布的大量早白垩世岩浆岩被认为是双向俯冲的重要证据(康志强等, 2009, 2010; Zhang et al., 2017)。Zhu et al. (2016)认为班-怒洋洋壳向南俯冲过程中在120~110Ma发生板片断离,在之后进入陆-陆碰撞阶段,拉萨与羌塘板块持续的拼贴导致了北拉萨板块地壳的加厚与区域显著构造缩减(Kapp et al., 2007)。另外,在研究区广泛发育一套晚白垩世的竟柱山组沉积建造,竟柱山组属于典型的磨拉石建造,厚度巨大,前人认为,竟柱山组是班-怒洋闭合后快速造山的产物,是古大洋闭合后转向陆相沉积环境的重要标志(李华亮等, 2016)。最新的研究结果表明,在研究区的竟柱山组砾岩中发现了88Ma的最小锆石U-Pb谐和年龄(数据未发表)(图 2e),这与拉木那勒岩体的侵位时代一致,这也进一步佐证了在该时期研究区发育显著的地壳增厚与抬升。因此,综上所述,我们认为拉木那勒岩体是班-怒洋闭合后,陆-陆碰撞加厚下地壳部分熔融的产物(图 11),而地壳增厚可能是解释北拉萨板块晚白垩世岩浆岩成因的一个重要形成机制。

5 结论

(1) 西藏尼玛县拉木那勒地区埃达克岩的侵位时代在88Ma左右,可能是北拉萨板块晚白垩世岩浆活动的产物。

(2) 拉木那勒埃达克岩具有高SiO2、Al2O3、Sr、Sr/Y比值,低Y、Yb,轻稀土富集,轻重稀土分异明显,无明显负Eu异常等特征,具有埃达克岩的地球化学特征。

(3) 拉木那勒埃达克岩可能是班-怒洋壳南向俯冲闭合后,陆-陆碰撞加厚下地壳部分熔融的产物,地壳增厚可能是解释北拉萨板块晚白垩世岩浆岩成因的一个重要机制。

致谢      衷心的感谢《岩石学报》副主编俞良军对本文的耐心指导和帮助。感谢中国地质科学院唐菊兴研究员、宋扬副研究员以及中国科学院地质与地球物理研究所兰中伍副研究员对本文的细心指导及提出的宝贵意见。感谢吉林大学西藏科研队的老师及同学们为本文提供的无私帮助。

参考文献
Anderson T. 2002. Correction of common lead in U-Pb analyses that do not report 204Pb. Chemical Geology, 192: 9-79
Barbey P, Allé P, Brouand M and Albarède F. 1995. Rare-earth patterns in zircons from the Manaslu granite and Tibetan Slab migmatites (Himalaya):Insights in the origin and evolution of a crustally-derived granite magma. Chemical Geology, 125(1-2): 1-17 DOI:10.1016/0009-2541(95)00068-W
Belousova E, Griffin W, O'Reilly SY and Fisher N. 2002. Igneous zircon:Trace element composition as an indicator of source rock type. Contributions to Mineralogy and Petrology, 143(5): 602-622 DOI:10.1007/s00410-002-0364-7
Boynton WV. 1984. Cosmochemistry of the rare earth elements:Meteorite studies. Developments in Geochemistry, 2: 63-114 DOI:10.1016/B978-0-444-42148-7.50008-3
Castillo PR, Janney PE and Solidum RU. 1999. Petrology and geochemistry of Camiguin island, southern Philippines:Insights to the source of adakites and other lavas in a complex arc setting. Contributions to Mineralogy and Petrology, 134(1): 33-51 DOI:10.1007/s004100050467
Castillo PR. 2006. An overview of adakite petrogenesis. Chinese Science Bulletin, 51(3): 257-268 DOI:10.1007/s11434-006-0257-7
Castillo PR. 2012. Adakite petrogenesis. Lithos, 134-135: 304-316
Chappell BW and White AJR. 1974. Two contrasting granite types. Pacific Geology, 8: 173-174
Chappell BW. 1999. Aluminium saturation in I- and S-type granites and the characterization of fractionated haplogranites. Lithos, 46(3): 535-551 DOI:10.1016/S0024-4937(98)00086-3
Chen JL, Xu JF, Yu HX, Wang BD, Wu JB and Feng YX. 2015. Late Cretaceous high-Mg# granitoids in southern Tibet:Implications for the early crustal thickening and tectonic evolution of the Tibetan Plateau?. Lithos, 232: 12-22 DOI:10.1016/j.lithos.2015.06.020
Chen RX, Zheng YF and Xie LW. 2010. Metamorphic growth and recrystallization of zircon:Distinction by simultaneous in-situ analyses of trace elements, U-Th-Pb and Lu-Hf isotopes in zircons from eclogite-facies rocks in the Sulu orogen. Lithos, 114(1-2): 132-154 DOI:10.1016/j.lithos.2009.08.006
Chen Y, Zhu DC, Zhao ZD, Meng FY, Wang Q, Santosh M, Wang LQ, Dong GC and Mo XX. 2014. Slab breakoff triggered ca.113Ma magmatism around Xainza area of the Lhasa Terrane, Tibet. Gondwana Research, 26(2): 449-463 DOI:10.1016/j.gr.2013.06.005
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
Clemens JD. 2003. S-type granitic magmas:Petrogenetic issues, models and evidence. Earth-Science Reviews, 61(1-2): 1-18 DOI:10.1016/S0012-8252(02)00107-1
Defant MJ and Drummond MS. 1990. Derivation of some modern arc magmas by melting of young subducted lithosphere. Nature, 347(6294): 662-665 DOI:10.1038/347662a0
Ding HX, Zhang ZM, Xiang H, Qi M, Gou ZB and Lei HC. 2015. The petrogenesis and tectonic significance of the Early Cretaceous volcanics from the northern Lhasa terrane, Tibet. Acta Petrologica Sinica, 31(5): 1247-1267 (in Chinese with English abstract)
Eyuboglu Y, Chung SL, Santosh M, Dudas FO and Akaryalı E. 2011. Transition from shoshonitic to adakitic magmatism in the eastern Pontides, NE Turkey:Implications for slab window melting. Gondwana Research, 19(2): 413-429 DOI:10.1016/j.gr.2010.07.006
Fan JJ, Li C, Wu H, Zhang TY, Wang M, Chen JW and Xu JX. 2016. Late Jurassic Adakitic granodiorite in the Dong Co area, northern Tibet:Implications for subduction of the Bangong-Nujiang oceanic lithosphere and related accretion of the Southern Qiangtang terrane. Tectonophysics, 691: 345-361 DOI:10.1016/j.tecto.2016.10.026
Grimes CB, John BE, Kelemen PB, Mazdab FK, Wooden J L, Cheadle MJ, Hanghøj K and Schwartz JJ. 2007. Trace element chemistry of zircons from oceanic crust:A method for distinguishing detrital zircon provenance. Geology, 35(7): 643-646
Guo F, Nakamuru E, Fan WM, Kobayoshi K and Li CW. 2007. Generation of Palaeocene adakitic andesites by magma mixing in Yanji Area, NE China. Journal of Petrology, 48(4): 661-692 DOI:10.1093/petrology/egl077
Gutscher MA, Maury R, Eissen JP and Bourdon E. 2000. Can slab melting be caused by flat subduction?. Geology, 28(6): 535-538 DOI:10.1130/0091-7613(2000)28<535:CSMBCB>2.0.CO;2
Hastie AR, Kerr AC, Pearce JA and Mitchell SF. 2007. Classification of altered volcanic island arc rocks using immobile trace elements:Development of the Th-Co discrimination diagram. Journal of Petrology, 48(12): 2341-2357 DOI:10.1093/petrology/egm062
Hoskin PWO and Schaltegger U. 2003. The composition of zircon and igneous and metamorphic petrogenesis. Reviews in Mineralogy and Geochemistry, 53(1): 27-62 DOI:10.2113/0530027
Hoskin PWO. 2005. Trace-element composition of hydrothermal zircon and the alteration of Hadean zircon from the Jack Hills, Australia. Geochimica et Cosmochimica Acta, 69(3): 637-648 DOI:10.1016/j.gca.2004.07.006
Hu PY, Li C, Wang M, Xie CM and Wu YW. 2013. Cambrian volcanism in the Lhasa terrane, southern Tibet:Record of an Early Paleozoic Andean-type magmatic arc along the Gondwana proto-Tethyan margin. Journal of Asian Earth Sciences, 77: 91-107 DOI:10.1016/j.jseaes.2013.08.015
Ji WQ, Wu FY, Chung SL, Li JX and Liu CZ. 2009. Zircon U-Pb geochronology and Hf isotopic constraints on petrogenesis of the Gangdese batholith, southern Tibet. Chemical Geology, 262(3-4): 229-245 DOI:10.1016/j.chemgeo.2009.01.020
Kang ZQ, Xu JF, Wang BD, Dong YH, Wang SQ and Chen JL. 2009. Geochemistry of Cretaceous volcanic rocks of Duoni Formation in northern Lhasa Block:Discussion of tectonic setting. Earth Science (Journal of China University of Geosciences), 34(1): 89-104 (in Chinese with English abstract) DOI:10.3799/dqkx.2009.009
Kang ZQ, Xu JF, Wang BD and Chen JL. 2010. Qushenla Formation volcanic rocks in north Lhasa block:Products of Bangong Co-Nujiang Tethy's southward subduction. Acta Petrologica Sinica, 26(10): 3106-3116 (in Chinese with English abstract)
Kapp P, DeCelles PG, Gehrels GE, Heizler M and Ding L. 2007. Geological records of the Lhasa-Qiangtang and Indo-Asian collisions in the Nima area of Central Tibet. Geological Society of America Bulletin, 119(7-8): 917-932 DOI:10.1130/B26033.1
Lai SC and Qin JF. 2013. Adakitic rocks derived from the partial melting of subducted continental crust:Evidence from the Eocene volcanic rocks in the northern Qiangtang block. Gondwana Research, 23(2): 812-824 DOI:10.1016/j.gr.2012.06.003
Le Bas MJ, Le Maitre RW, Streckeisen A and Zanettin B. 1986. A chemical classification of volcanic-rocks based on the total alkali silica diagram. Journal of Petrology, 27: 745-750 DOI:10.1093/petrology/27.3.745
Li HL, Gao C, Li ZH, Zhang Z, Peng ZM and Guan JL. 2016. Age and tectonic Significance of Jingzhushan Formation in Bangong Lake area, Tibet. Geotectonica et Metallogenia, 40(4): 663-673 (in Chinese with English abstract)
Liu YM, Wang M, Li C, Li SZ, Xie CM, Zeng XW, Dong YC and Liu JH. 2019. Late Cretaceous tectono-magmatic activity in the Nize region, central Tibet:Evidence for lithospheric delamination beneath the Qiangtang-Lhasa collision zone. International Geology Review, 61(5): 562-583 DOI:10.1080/00206814.2018.1437789
Ludwig KJ. 2003. ISOPLOT 3.0. Vol.4. Berkeley Geochronology Center Special Publication, 4: 1-70
Ma GL and Yue YH. 2010. Cretaceous volcanic rocks in northern Lhasa Block:Constraints on the tectonic evolution of the Gangdise Arc. Acta Petrologica et Mineralogica, 29(5): 525-538 (in Chinese with English abstract)
Martin H. 1999. Adakitic magmas:Modern analogues of Archaean granitoids. Lithos, 46(3): 411-429 DOI:10.1016/S0024-4937(98)00076-0
Martin H, Smithies RH, Rapp R, Moyen JF and Champion D. 2005. An overview of adakite, tonalite-trondhjemite-granodiorite (TTG), and sanukitoid:Relationships and some implications for crustal evolution. Lithos, 79(1-2): 1-24 DOI:10.1016/j.lithos.2004.04.048
Mo XX, Niu YL, Dong GC, Zhao ZD, Hou ZQ, Zhou S and Ke S. 2008. Contribution of syncollisional felsic magmatism to continental crust growth:A case study of the Paleogene Linzizong volcanic succession in southern Tibet. Chemical Geology, 250(1-4): 49-67 DOI:10.1016/j.chemgeo.2008.02.003
Molnar P and Tapponnier P. 1975. Cenozoic tectonics of Asia:Effects of a continental collision:Features of recent continental tectonics in Asia can be interpreted as results of the India-Eurasia collision. Science, 189(4201): 419-426 DOI:10.1126/science.189.4201.419
Pan GT, Chen ZL, Li XZ, Yan YJ, Xu XS, Xu Q, Jiang XS, Wu YL, Luo JN, Zhu TX and Peng YM. 1997. Geological-tectonic Evolution in the Eastern Tethys. Beijing: Geological Publishing House, 1-218 (in Chinese)
Powell CM and Conaghan PJ. 1973. Plate tectonics and the Himalayas. Earth and Planetary Science Letters, 20(1): 1-12 DOI:10.1016/0012-821X(73)90134-9
Prouteau G, Scaillet B, Pichavant M and Maury R. 2001. Evidence for mantle metasomatism by hydrous silicic melts derived from subducted oceanic crust. Nature, 410(6825): 197-200 DOI:10.1038/35065583
Rapp RP, Shimizu N, Norman MD and Applegate GS. 1999. Reaction between slab-derived melts and peridotite in the mantle wedge:Experimental constraints at 3.8GPa. Chemical Geology, 160(4): 335-356 DOI:10.1016/S0009-2541(99)00106-0
Schiano P, Monzier M, Eissen JP, Martin H and Koga KT. 2010. Simple mixing as the major control of the evolution of volcanic suites in the Ecuadorian Andes. Contributions to Mineralogy and Petrology, 160(2): 297-312 DOI:10.1007/s00410-009-0478-2
Smithies RH. 2000. The Archaean tonalite-trondhjemite-granodiorite (TTG) series is not an analogue of Cenozoic adakite. Earth and Planetary Science Letters, 182(1): 115-125 DOI:10.1016/S0012-821X(00)00236-3
Streck MJ, Leeman WP and Chesley J. 2007. High-magnesian andesite from Mount Shasta:A product of magma mixing and contamination, not a primitive mantle melt. Geology, 35(4): 351-354
Sui QL, Wang Q, Zhu DC, Zhao ZD, Chen Y, Santosh M, Hu ZC, Yuan HL and Mo XX. 2013. Compositional diversity of ca. 110Ma magmatism in the northern Lhasa Terrane, Tibet:Implications for the magmatic origin and crustal growth in a continent-continent collision zone. Lithos, 168-169: 144-159 DOI:10.1016/j.lithos.2013.01.012
Sun GY, Hu XM, Zhu DC, Hong WT, Wang JG and Wang Q. 2015. Thickened juvenile lower crust-derived~90Ma adakitic rocks in the central Lhasa terrane, Tibet. Lithos, 224-225: 225-239 DOI:10.1016/j.lithos.2015.03.010
Sun SS and McDonough WF. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. In: Sanders AD and Norry MJ (eds.). Magmatism in Ocean Basins. Geological Society of London Special Publications, 42(1): 313-345
Wang Q, McDermott F, Xu JF, Bellon H and Zhu YT. 2005. Cenozoic K-rich adakitic volcanic rocks in the Hohxil area, northern Tibet:Lower-crustal melting in an intracontinental setting. Geology, 33(6): 465-468 DOI:10.1130/G21522.1
Wang Q, Xu JF, Zhao ZH, Zi F, Tang GJ, Jia XH and Jiang ZQ. 2008. Tectonic setting and associated rock suites of adakitic rocks. Bulletin of Mineralogy, Petrology and Geochemistry, 27(4): 344-350 (in Chinese with English abstract)
Wang Q, Zhu DC, Zhao ZD, Liu SA, Chung SL, Li SM, Liu D, Dai JG, Wang LQ and Mo XX. 2014. Origin of the ca. 90Ma magnesia-rich volcanic rocks in SE Nyima, central Tibet:Products of lithospheric delamination beneath the Lhasa-Qiangtang collision zone. Lithos, 198-199: 24-37 DOI:10.1016/j.lithos.2014.03.019
Wang ZZ, Liu D, Zhao ZD, Yan JJ, Shi QS and Mo XX. 2017. The Sangri highly fractionated I-type granites in southern Gangdese:Petrogenesis and dynamic implication. Acta Petrologica Sinica, 33(8): 2479-2493 (in Chinese with English abstract)
Wen DR, Chung SL, Song B, Iizuka Y, Yang HJ, Ji JQ, Liu DY and Gallet S. 2008. Late Cretaceous Gangdese intrusions of adakitic geochemical characteristics, SE Tibet:Petrogenesis and tectonic implications. Lithos, 105(1-2): 1-11 DOI:10.1016/j.lithos.2008.02.005
Wu H, Li C, Xu MJ and Li XK. 2015. Early cretaceous adakitic magmatism in the Dachagou Area, northern Lhasa terrane, Tibet:Implications for slab roll-back and subsequent slab break-off of the lithosphere of the Bangong-Nujiang Ocean. Journal of Asian Earth Sciences, 97: 51-66 DOI:10.1016/j.jseaes.2014.10.014
Wu H, Xie CM, Li C, Wang M, Fan JJ and Xu WL. 2016. Tectonic shortening and crustal thickening in subduction zones:Evidence from Middle-Late Jurassic magmatism in southern Qiangtang, China. Gondwana Research, 39: 1-13 DOI:10.1016/j.gr.2016.06.009
Wu ZH, Wu XW, Zhao Z, Lu L, Ye PS and Zhang YL. 2014. Shrimp U-Pb isotopic dating of the Late Cretaceous volcanic rocks and its chronological constraint on the red-beds in southern Qiangtang Block. Acta Geoscientica Sinica, 35(5): 567-572 (in Chinese with English abstract)
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)
Yin A and Harrison TM. 2000. Geologic evolution of the Himalayan-Tibetan Orogen. Annual Review of Earth and Planetary Sciences, 28: 211-280 DOI:10.1146/annurev.earth.28.1.211
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 Geoanalytical Research, 28: 3-370
Yu HX, Chen JL, Xu JF, Wang BD, Wu JB and Liang YH. 2011. Geochemistry and origin of Late Cretaceous (~90Ma) ore-bearing porphyry of Balazha in mid-northern Lhasa terrane, Tibet. Acta Petrologica Sinica, 27(7): 2011-2022 (in Chinese with English abstract)
Zhai QG, Jahn BM, Wang J, Su L, Mo XX, Wang KL, Tang SH and Lee HY. 2013. The carboniferous ophiolite in the middle of the Qiangtang terrane, Northern Tibet:SHRIMP U-Pb dating, geochemical and Sr-Nd-Hf isotopic characteristics. Lithos, 168-169: 186-199 DOI:10.1016/j.lithos.2013.02.005
Zhang S, Shi HF, Hao HJ, Li DW, Lin Y and Feng MX. 2014. Geochronology, geochemistry and tectonic significance of Late Cretaceous adakites in Bangong Lake, Tibet. Earth Science (Journal of China University of Geosciences), 39(5): 509-524 (in Chinese with English abstract)
Zhang YX, Li ZW, Yang WG, Zhu LD, Jin X, Zhou XY, Tao G and Zhang KJ. 2017. Late Jurassic-Early Cretaceous episodic development of the Bangong Meso-Tethyan subduction:Evidence from elemental and Sr-Nd isotopic geochemistry of arc magmatic rocks, Gaize region, central Tibet, China. Journal of Asian Earth Sciences, 135: 212-242 DOI:10.1016/j.jseaes.2016.12.043
Zhang ZM, Zhao GC, Santosh M, Wang JL, Dong X and Shen K. 2010. Late Cretaceous charnockite with adakitic affinities from the Gangdese batholith, southeastern Tibet:Evidence for Neo-Tethyan mid-ocean ridge subduction?. Gondwana Research, 17(4): 615-631 DOI:10.1016/j.gr.2009.10.007
Zhu DC, Pan GT, Mo XX, Wang LQ, Liao ZL, Zhao ZD, Dong GC and Zhou CY. 2006. Late Jurassic-Early Cretaceous geodynamic setting in middle-northern Gangdese:New insights from volcanic rocks. Acta Petrologica Sinica, 22(3): 534-546 (in Chinese with English abstract)
Zhu DC, Pan GT, Wang LQ, Mo XX, Zhao ZD, Zhou CY, Liao ZL, Dong GC and Yuan SH. 2008. Tempo-spatial variations of Mesozoic magmatic rocks in the Gangdise belt, Tibet, China, with a discussion of geodynamic setting-related issues. Geological Bulletin of China, 27(9): 1535-1550 (in Chinese with English abstract)
Zhu DC, Pan GT, Zhao ZD, Pan GT, Lee HY, Kang ZQ, Liao ZL, Wang LQ, Li GM, Dong GC and Liu B. 2009. Early Cretaceous subduction-related adakite-like rocks of the Gangdese Belt, southern Tibet:Products of slab melting and subsequent melt-peridotite interaction?. Journal of Asian Earth Sciences, 34(3): 298-309 DOI:10.1016/j.jseaes.2008.05.003
Zhu DC, Zhao ZD, Niu YL, Dilek Y, Hou ZQ and Mo XX. 2013. The origin and pre-Cenozoic evolution of the Tibetan Plateau. Gondwana Research, 23(4): 1429-1454 DOI:10.1016/j.gr.2012.02.002
Zhu DC, Xia Y, Qiu BB, Wang Q and Zhao ZD. 2013. Why do we need to propose the Early Cretaceous Comei large igneous province in southeastern Tibet?. Acta Petrologica Sinica, 29(11): 3659-3670 (in Chinese with English abstract)
Zhu DC, Li SM, Cawood PA, Wang Q, Zhao ZD, Liu SA and Wang LQ. 2016. Assembly of the Lhasa and Qiangtang terranes in central Tibet by divergent double subduction. Lithos, 245: 7-17
丁慧霞, 张泽明, 向华, 祁敏, 苟正彬, 雷恒聪. 2015. 青藏高原拉萨地体北部早白垩世火山岩的成因及意义. 岩石学报, 31(5): 1247-1267.
康志强, 许继峰, 王保弟, 董彦辉, 王树庆, 陈建林. 2009. 拉萨地块北部白垩纪多尼组火山岩的地球化学:形成的构造环境. 地球科学(中国地质大学学报), 34(1): 89-104.
康志强, 许继峰, 王保弟, 陈建林. 2010. 拉萨地块北部去申拉组火山岩:班公湖-怒江特提斯洋南向俯冲的产物?. 岩石学报, 26(10): 3106-3116.
李华亮, 高成, 李正汉, 张璋, 彭智敏, 关俊雷. 2016. 西藏班公湖地区竟柱山组时代及其构造意义. 大地构造与成矿学, 40(4): 663-673.
马国林, 岳雅慧. 2010. 西藏拉萨地块北部白垩纪火山岩及其对冈底斯岛弧构造演化的制约. 岩石矿物学杂志, 29(5): 525-538. DOI:10.3969/j.issn.1000-6524.2010.05.008
潘桂棠, 陈智梁, 李兴振, 颜仰基, 许效松, 徐强, 江新胜, 吴应林, 罗建宁, 朱同兴, 彭勇民. 1997. 东特提斯地质构造形成演化. 北京: 地质出版社, 1-218.
王强, 许继峰, 赵振华, 资锋, 唐功建, 贾小辉, 姜子琦. 2008. 埃达克质岩的构造背景与岩石组合. 矿物岩石地球化学通报, 27(4): 344-350. DOI:10.3969/j.issn.1007-2802.2008.04.003
王珍珍, 刘栋, 赵志丹, 闫晶晶, 石卿尚, 莫宣学. 2017. 冈底斯带南部桑日高分异Ⅰ型花岗岩的岩石成因及其动力学意义. 岩石学报, 33(8): 2479-2493.
吴珍汉, 吴学文, 赵珍, 陆露, 叶培盛, 张耀玲. 2014. 羌塘地块南部晚白垩世火山岩离子探针测年及其对红层时代的约束. 地球学报, 35(5): 567-572.
许志琴, 杨经绥, 李文昌, 李化启, 蔡志慧, 闫臻, 马昌前. 2013. 青藏高原中的古特提斯体制与增生造山作用. 岩石学报, 29(6): 1847-1860.
余红霞, 陈建林, 许继峰, 王宝弟, 邬建斌, 梁华英. 2011. 拉萨地块中北部晚白垩世(约90Ma)拔拉扎含矿斑岩地球化学特征及其成因. 岩石学报, 27(7): 2011-2022.
张硕, 史洪峰, 郝海健, 李德威, 吝岩, 冯旻譞. 2014. 青藏高原班公湖地区晚白垩世埃达克岩年代学、地球化学及构造意义. 地球科学(中国地质大学学报), 39(5): 509-524.
朱弟成, 潘桂棠, 莫宣学, 王立全, 廖忠礼, 赵志丹, 董国臣, 周长勇. 2006. 冈底斯中北部晚侏罗世-早白垩世地球动力学环境:火山岩约束. 岩石学报, 22(3): 534-546.
朱弟成, 潘桂棠, 王立全, 莫宣学, 赵志丹, 周长勇, 廖忠礼, 董国臣, 袁四化. 2008. 西藏冈底斯带中生代岩浆岩的时空分布和相关问题的讨论. 地质通报, 27(9): 1535-1550. DOI:10.3969/j.issn.1671-2552.2008.09.013
朱弟成, 夏瑛, 裘碧波, 王青, 赵志丹. 2013. 为什么要提出西藏东南部早白垩世措美大火成岩省. 岩石学报, 29(11): 3659-3670.