岩石学报  2017, Vol. 33 Issue (8): 2509-2522   PDF    
藏东同普二叠纪高分异花岗岩的锆石U-Pb年龄和岩石成因
周士旭1, 朱弟成1,2, 张亮亮1, 谢锦程1, 王青1     
1. 地质过程与矿产资源国家重点实验室, 中国地质大学地球科学与资源学院, 北京 100083;
2. 中国科学院青藏高原地球科学卓越研究中心, 北京 100101
摘要: 对具有富U或高U锆石的高硅花岗岩进行准确定年并进行成因讨论是花岗岩研究中的难题。本文以青藏高原东部江达-维西构造带北部同普岩基边缘相的高硅花岗岩为对象,开展了锆石U-Pb年代学、锆石Hf同位素和全岩元素地球化学研究。采用大激光束斑(100μm)先剥蚀2~3个脉冲的预剥蚀方法,获得了这套含富U或高U锆石高硅花岗岩的可靠锆石U-Pb年龄(260±1Ma)。这套高硅花岗岩具有高的SiO2含量(74.92%~76.46%)、高的全碱含量(Na2O+K2O;7.61%~8.55%)和高的分异指数(92~96),明显的Eu(Eu/Eu*=0.17~0.41)异常和显著亏损Ba、Nb、Sr、P和Ti等特征,属于典型的高分异花岗岩。本文数据,结合文献数据,指示同普高硅花岗岩为高分异I型花岗岩。富集的锆石Hf同位素组成(εHft)=-9.9~-5.8)指示它们来源于古老地壳物质的部分熔融。定量模拟表明它们可能由同期花岗闪长质母岩浆经历斜长石、磷灰石和富钛矿物的分离结晶作用形成。
关键词: 藏东     同普岩基     高分异花岗岩     锆石U-Pb年龄     岩石成因    
Zircon U-Pb age and petrogenesis of the Permian highly fractionated granites in Tongpu, eastern Tibet
ZHOU ShiXu1, ZHU DiCheng1,2, ZHANG LiangLiang1, XIE JinCheng1, WANG Qing1     
1. State Key Laboratory of Geological Processes and Mineral Resources, and School of Earth Science and Resources, China University of Geosciences, Beijing 100083, China;
2. Center for Excellence in Tibetan Plateau Earth Sciences, Chinese Academy of Sciences, Beijing 100101, China
Abstract: How to precisely date the high-silica granites with U-rich or high-U zircons and how to effectively constrain their petrogenesis remain two puzzling questions in the research of granite. This paper reports zircon U-Pb ages and Hf isotope, whole-rock major and trace elemental data of the high-silica granites from the marginal facies of the Tongpu batholith in northern Jomda-Weixi belt, eastern Tibet. Reliable zircon U-Pb age data (260±1Ma) of the U-rich or high-U zircons from the high-silica granite were obtained by pre-ablation dating method, i.e., conducting 2~3 pulses ablation firstly with a large size beam spot (100μm) before normal ablation. These samples are characterized by high SiO2 (74.92%~76.46%), high Na2O+K2O (7.61%~8.55%), high differentiation index (92~96), and significant negative Eu anomalies (Eu/Eu*=0.17~0.41) and strong depletion of Ba, Nb, Sr, P, and Ti, similar to those of highly fractionated granites. These data, along with the data in the literature, indicate that the Tongpu high-silica granites can be assigned to be highly fractionated I-type granites. The negative zircon εHf(t) (-9.9~-5.8) values indicate that they are most likely derived from partial melting of ancient crustal material. Quantitative modeling shows that the Tongpu high-silica granites can be formed by fractional crystallization of plagioclase, apatite, and Ti-rich minerals from the coeval granodioritic magma.
Key words: Zircon U-Pb age     Petrogenesis     Highly fractionated granites     Tongpu batholith     Eastern Tibet    

在地球显生宙不同时代花岗岩体中,均发育数量不等的高硅花岗岩(SiO2 > 70%)。准确限定这些高硅花岗岩的侵位时代和岩浆演化过程,不但有助于了解其构造含义,而且也有利于揭示大陆地壳演化过程的元素传输过程和相关成矿作用机理(Lee and Morton, 2015Wu et al., 2017)。锆石U-Pb定年是限定花岗岩侵位时代的常见方法,但对高硅花岗岩而言,岩浆锆石常具有高的U、Th含量,易于发生蜕晶化(即形成所谓的“黑”锆石)而造成放射性成因Pb的丢失或获取,使得很难用常规的锆石U-Pb定年方法获得代表岩浆侵位的年龄(吴福元等,2015; Wu et al., 2017)。为此,一些研究者提出了采用物理磨蚀去掉锆石颗粒外层的“空蚀法”(Krogh,1982)和寻找其它合适的定年矿物(如独居石、磷灰石、榍石等)(李秋立,2016)来限定高硅花岗岩的侵位时代。另一方面,学术界对高硅花岗岩究竟是由母岩浆经历高程度分离结晶形成还是其本身就代表了母岩浆成分,也存在争议(Wu et al., 2017)。争议的原因在于:(1) 花岗质岩浆粘度大,不利于矿物的自由生长,也不利于先期结晶的矿物在岩浆中发生大尺度的移动和下沉,主要造岩矿物(例如斜长石)的密度与花岗质岩浆的密度相差无几,分离结晶作用很难进行(Pitcher, 1997张旗等,2007Wu et al., 2017);(2) 实验岩石学结果显示,不同类型沉积源岩(包括泥质岩、杂砂岩等)在不同物理化学条件下发生部分熔融(Scaillet et al., 1995Patiño Douce and Harris, 1998),也能形成高硅花岗岩。在这种情况下,有效探讨高硅花岗岩成因的关键,在于能否从地质上找到和高硅花岗岩具有成因联系的更镁铁质岩石序列(吴福元等,2015),并通过地球化学方法建立起高硅花岗岩和这些更镁铁质岩石之间的成因联系。

青藏高原东部三江地区广泛发育不同时代的各种花岗岩类(莫宣学等,1993邓军等,2011王舫等,2013Deng et al., 2014a, b, 2015, 2017),其中江达-维西构造带同普花岗岩基(图 1a, b)是代表性的二叠纪花岗岩基之一。已有研究显示,同普岩基的岩石类型较为复杂,包括石英闪长岩、花岗闪长岩、黑云母花岗岩和斑状花岗岩等(吴涛等,2013)。重要的是,这些黑云母花岗岩和斑状花岗岩发育在同普花岗岩基的边缘(图 1c),SiO2含量高(>75%),这种岩石组合和地质特征提供了研究高硅花岗岩侵位时代和成因的难得机会。因此,本文目的是,以同普岩基边缘相的高硅花岗岩为对象,采用“空蚀法”(Krogh,1982),限定这些高硅花岗岩的侵位时代,利用本文获得的新的全岩地球化学和锆石Hf同位素数据,并结合同普岩基报道的其它数据,讨论同普岩基边缘相高硅花岗岩的成因。

图 1 研究区地质简图和样品野外照片 (a)青藏高原-三江地区构造单元划分(据Zhu et al., 2013Deng et al., 2014a, b);(b)三江地区地质简图(据Deng et al., 2014a, b;谢尧武等,2007改绘);(c)同普岩基地质简图(据吴涛等,2013;谢尧武等,2007);(d)同普岩基和下三叠统普水桥组接触界面;(e)样品14TP03-1野外露头照片;(f)样品14TP03-1照片;(g)样品14TP03-1显微照片 Fig. 1 Simplified geological map of the studied area and petrographical photos of samples (a) tectonic subdivision of the Tibetan Plateau-Sanjiang Region (modified after Zhu et al., 2013; Deng et al., 2014a, b); (b) simplified geological map of Sanjiang area (modified after Deng et al., 2014a, b); (c) simplified geological map of Tongpu Batholith (modified after Wu et al., 2013) (d) contact interface between Tongpu Batholith and Pushuiqiao Formation (Lower Triassic); (e) outcrops of sample 14TP03-1; (f) photo of sample 14TP03-1; (g) petrographical photo of sample 14TP03-1

① 谢尧武,彭兴阶,陈应明,陈德泉等. 2007.中华人民共和国1:25万囊谦县、昌都县、江达县幅区域地质调查报告

1 地质背景和样品描述

青藏高原东部的三江构造域,物质组成非常复杂,不同混杂岩带或蛇绿岩带分隔了众多微陆块或地体。例如,在三江构造域北段,主要由西羌塘地体、昌都地体和义敦岛弧等地体组成,其间分别间夹着北澜沧江缝合带、金沙江缝合带等(图 1a)(王立全等,2000Hou et al., 2007Deng et al., 2014a, b)。本文所研究的江达-维西构造带,位于昌都地体东缘,紧邻金沙江缝合带,南北延伸约500km(图 1a)。带内最古老的地层是以片麻岩和变粒岩为代表的古-中元古代宁多群,泥盆系-侏罗系地层广泛发育(谢尧武等,2007)。江达-维西构造带内岩浆活动较为发育,主要为早二叠世-晚三叠世,并存在3个明显的岩浆活动峰值(约265Ma,约248Ma和约235Ma)(莫宣学等,2001Zi et al., 2013Wang et al., 2014a)。早二叠世晚期-晚二叠世岩浆岩分布较零散,主要由基性-中酸性岩石组成(李兴振等,1998王立全等,2000)。早-中三叠世岩浆活动以中酸性岩石为主,带内广泛发育的晚三叠世岩浆活动具有明显的双峰式岩浆岩特征(罗建宁,1992莫宣学等,1993王立全等,2000)。

江达-维西构造带记录了金沙江洋复杂的演化历史。前人研究认为金沙江洋于早石炭世打开,自晚石炭世开始向昌都地体之下俯冲,形成陆缘弧(何龙清,1998王立全等,2000彭勇民等,2000董树义等, 2006, 2010),中三叠世末昌都地体与德格-中甸地体碰撞,江达-维西构造带成为三江复合造山带的一部分(刘肇昌等,1995汪啸风等,1999王立全等,2003邓军等,2013)。

同普岩基位于江达-维西构造带北段(图 1b),呈北西-南东向展布,出露面积约400km2,为藏东地区最大的二叠纪侵入杂岩体(图 1c)。同普岩基主要侵位于上泥盆统东拉组(D3d)砂岩、板岩、中基性火山岩中,其上被下三叠统普水桥组(T1p)砂岩和砾岩覆盖(图 1d-g)。本文样品采自同普岩基北侧,位于同普乡八一桥南约5km处(图 1d)。岩体经历了较强烈的风化蚀变,风化面呈淡黄色,新鲜面呈灰白色。岩石类型主要为正长花岗岩,浅肉红色中粒半自形结构,块状构造。矿物主要包括:钾长石(约为45%),石英(约为25%),斜长石(约为20%),黑云母和白云母(约为5%)。副矿物主要有锆石、磷灰石和榍石(图 1g)。

2 分析方法

本文对同普岩基2件正长花岗岩样品进行了激光锆石U-Pb定年并对其中1件正长花岗岩样品进行锆石原位Hf同位素分析。样品碎样和锆石分选是在河北廊坊市宇能岩石矿物分选技术服务有限公司采用传统的浮选和电磁选方法处理的。锆石阴极发光(CL)显微照相在中国科学院地质与地球物理研究所扫描电镜室完成的,仪器型号为聚焦离子束-扫描电子显微镜(FIB/SEM)双束系统-Zeiss Auriga Compact。锆石LA-ICP-MS U-Pb定年测试工作在中国地质大学(北京)矿物激光微区分析实验室(Milma Lab)完成,测试仪器分别为NewWave 193UC型准分子激光器和Agilent 7900四级杆型等离子质谱仪(ICP-MS)。测试过程中,激光剥蚀斑束直径为35μm,剥蚀频率为8Hz,激光能量为6J/cm2。锆石年龄计算采用国际标准锆石91500作为外标,元素含量采用美国国家标准物质局人工合成硅酸盐玻璃NIST610作为外标,29Si作为内标元素进行校正。数据校正采用ICPMSDataCal软件来处理样品的同位素比值和元素含量数据(Liu et al., 2008a, b, 2010),普通铅校正是在Andersen(2002)程序中完成的,最终的年龄计算和谐和图绘制采用ISOPLOT(2.49版;Ludwig, 2000)宏程序完成。

锆石Hf同位素分析是在中国地质大学(武汉)地质过程与矿产资源重点实验室完成的,分析仪器为Neptune plus MC-ICP-MS和GeoLas 2005激光剥蚀系统。锆石Hf同位素测试在锆石U-Pb定年相同的分析位置进行。测试过程中,激光束直径为44μm,激光脉冲频率为6Hz。整个测试以91500、GJ-1和Monastery作为标样监测数据质量。标样91500的176Hf/177Hf测试结果为0.282291±35(2σ,n=89),GJ-1的176Hf/177Hf测试结果为0.282019±4(2σ,n = 87),均在误差范围内与推荐值一致(Woodhead and Hergt, 2005)。详细的仪器分析条件和数据获取方法见文献(Liu et al., 2010)。

全岩主量元素分析在中国科学院地质与地球物理研究所采用XRF法完成,分析精度优于5%。微量元素是在武汉上谱分析科技有限责任公司分析测试中心利用电感耦合等离子体质谱法(ICP-MS)完成,分析精度为5%~10%。测试过程中,全部样品在Agilent 7700e ICP-MS上测定,分析过程中采用标样BHVO-2、BCR-2和GSP进行质量监控(Liu et al., 2008a, b)。

3 分析结果

本文对2件同普花岗岩样品(14TP02-1和14TP03-1) 进行了锆石U-Pb定年(表 1)。这两件样品的锆石颗粒大多呈自形-半自形长柱或短柱状(长轴方向80~150μm),部分颗粒发育清晰的震荡环带,多数锆石显示黑色CL图像(图 2),具富U或高U锆石特征。本文先后采用了两种方式对这些锆石进行激光U-Pb年龄测定:(1) 按照常规激光参数直接对锆石进行剥蚀,获得的Th和U含量分别为517×10-6和1117×10-6(图 3a),18个测点206Pb/238U年龄数据变化范围大(200~280Ma)(图 3b);(2) 采用激光预剥蚀的方法进行了第二次分析,即先采用大的激光束斑(100μm)对锆石表面剥蚀2~3个脉冲,去除因放射性损伤导致的具有隐伏裂隙的锆石颗粒,然后再对无明显裂隙的锆石颗粒进行常规激光剥蚀,获得了非常好的数据结果(图 3c)。22个测点的Th和U含量分别为366×10-6和792×10-6(图 3a),206Pb/238U年龄变化于259~262Ma,加权平均年龄为260±1Ma(n=22)。采用预剥蚀方法对样品14TP03-1进行了测试,由于锆石颗粒较小,24个测点中17个打穿未获得有效年龄,剩下7个测点提供了较为一致的206Pb/238U加权平均年龄(260±2Ma,图 3d)。

表 1 藏东江达同普岩基正长花岗岩的锆石U-Pb年龄数据 Table 1 Zircon U-Pb age data of the syengranite from the Tongpu batholith, Jomda in eastern Tibet

图 2 同普岩基样品锆石预剥蚀前后反射图像对比图和锆石CL图像特征 Fig. 2 Comparison of reflection images of zircons before and after pre-ablation and cathodoluminescence (CL) images of zircons from samples from the Tongpu batholith

图 3 同普岩基高U锆石两次分析结果 (a)样品14TP02-1两次测试结果的U含量和年龄关系图;(b)样品14TP02-1第一次锆石U-Pb年龄谐和图;(c)样品14TP02-1第二次锆石U-Pb年龄谐和图;(d)样品14TP03-1第二次锆石U-Pb年龄谐和图 Fig. 3 Comparison of two results of high-U zircon and results of zircon U-Pb data of Tongpu batholith (a) U vs. age diagram of two-round running; (b) concordia diagram of zircons from sample 14TP02-1 during the first running; (c) concordia diagram of zircons from sample 14TP02-1 during the second running; (d) concordia diagram of zircons from sample 14TP03-1 during the second running

对样品14TP03-1共进行了20个测试点的锆石的Hf同位素测试(表 2),获得的176Lu/177Hf值为0.001273~0.005926,176Hf/177Hf值为0.282466~0.282357,对应的εHf(t)值为-9.9~-5.8,地壳模式年龄(tDMC)为2.0~1.8Ga。

表 2 藏东江达同普岩基正长花岗岩的锆石Hf同位素分析结果 Table 2 ircon Hf-isotope data of the syengranite from the Tongpu batholith, Jomda in eastern Tibet

本文获得的同普岩基样品,均为正长花岗岩(表 3图 4a),以高硅(SiO2=75.92%~77.42%)、高全碱(Na2O+K2O=7.61%~8.55%)和较高K2O/Na2O(1.05~1.43) 比值为特征。这些样品为高钾钙碱性系列(图 4b),并具有非常高的分异指数(DI=92~96)(图 4c)。这些样品显示强烈的负Eu异常(Eu/Eu*=0.17~0.41)(图 5a)和显著的负Ba、Nb、Sr、P和Ti异常(图 5b)。

表 3 藏东江达同普岩基正长花岗岩的主量元素(wt%)和微量元素(×10-6)分析结果 Table 3 Whole-rock major (wt%) and trace element (×10-6) data of the syengranite from the Tongpu batholith, Jomda in eastern Tibet

图 4 同普岩基样品的地球化学图解 (a)R1-R2分类图解(De la Roche et al., 1980);(b)SiO2-K2O图解(Rickwood, 1989);(c)分异指数(DI)-A/CNK图解;花岗闪长岩样品、石英闪长岩样品、黑云母花岗岩和斑状花岗岩样品均引自(吴涛等,2013) Fig. 4 Geochemical diagram of samples from the Tongpu batholith (a) R2 vs. R1 diagram (De la Roche et al., 1980) for classification; (b) K2O vs. SiO2 (Rickwood, 1989); (c) A/CNK vs. differentiation index (DI) plots showing sample compositional variation. A/CNK = Al2O3/(CaO+Na2O+K2O) (molar ratio), DI = Quartz (Qtz)+Orthoclase (Or)+Albite (Ab)+Nepheline (Ne)+Leucite (Lc)+K-feldspar (Kfs), from CIPW calculating values; Data sources: granodiorite, quartz diorite, biotite granite, and porphyritic granite (Wu et al., 2013)

图 5 同普岩基样品的球粒陨石标准化稀土元素配分模式图(a)和原始地幔标准化微量元素配分图(b)(标准化值据Sun and McDonough, 1989) Fig. 5 Chondrite-normalized REE pattern (a) and primitive mantle-normalized trace-element spidergram (b) of samples from the Tongpu batholith (normalization values after Sun and McDonough, 1989)
4 讨论 4.1 同普岩基高硅花岗岩的年龄

黑色的锆石CL图像(图 2)和变化的U含量(图 3a)均表明,本文正长花岗岩样品中的锆石为富U或高U锆石。高U含量的锆石常出现在特定岩石类型中,如喜马拉雅淡色花岗岩和华南含W-Zn矿花岗岩(李秋立,2016),这些花岗岩被解释为高分异花岗岩(李献华等,2007吴福元等,2015)。这种特征的锆石很难用来开展U-Pb定年工作,这是因为这种锆石在形成时,高含量的U、Th元素发生放射性衰变,当这种放射性损伤超过临界点时,将破坏锆石颗粒的内部结构,且随着时间积累,破坏会加剧(Ewing et al., 2003; Farnan et al., 2003),形成裂隙。裂隙的存在使得放射性成因铅的活动(Pb丢失或得到)增强(Wang et al., 2014b),导致所获得的U-Pb表观年龄会随着U含量的升高而变低(李秋立,2016)。

已有研究提出了两种方法来解决这类富U或高U锆石的定年问题,一是长期使用的“空蚀法”,即通过物理磨蚀去掉锆石颗粒的外层,再测定剩余的中间部分以获得更为谐和的年龄(Krogh,1982),二是采用其它不会造成严重放射性损伤而发生显著Pb丢失的矿物(如斜锆石、独居石、磷灰石、榍石等)(李秋立,2016)。受第一种方法启示,为有效避免无法利用反射光、透射光图像和CL图像识别锆石颗粒中隐伏裂隙和放射性损伤部位的缺陷,本文在实验过程中,先用大的激光束(100μm)对样品14TP02-1中选定的锆石颗粒剥蚀2~3个脉冲(即预剥蚀),结果获得了可靠的LA-ICP-MS锆石U-Pb谐和年龄(260±1Ma)(图 3a, b)。这一年龄与已有研究在同普黑云母花岗岩(263.7±1.6Ma)、花岗闪长岩(263.9±1.9Ma)和石英闪长岩(262.8±1.5Ma)中获得的锆石U-Pb年龄(吴涛等,2013)在误差范围内一致。这些年龄数据表明,同普岩基中的高硅花岗岩,与其它岩性同期侵位于约260Ma。

4.2 同普岩基高硅花岗岩的成因类型

本文获得的同普岩基边缘相正长花岗岩样品,主要由长石和石英组成,见少量黑云母,未见角闪石。这些样品具有高的SiO2含量(75.92%~77.42%)、高的分异指数(DI=92~96) 和低的P2O5含量(0.01%~0.04%),表明它们经历了高程度结晶分异。样品同时还具有强烈的Eu负异常(Eu/Eu*=0.17~0.41)(图 5a)和显著的Ba、Nb、Sr、P和Ti的负异常(图 5b),在Zr+Nb+Ce+Y-(Na2O+K2O)/CaO和Zr+Nb+Ce+Y-FeOT/MgO图解中位于高分异花岗岩区域(图 6a, b)。显然,地球化学特征显示同普岩基边缘相的高硅花岗岩为高分异花岗岩。

图 6 同普岩基样品的Zr+Nb+Ce+Y-(Na2O+K2O)/CaO (a)和Zr+Nb+Ce+Y-FeOT/MgO (b)(据Whalen et al., 1987) 花岗闪长岩样品、黑云母花岗岩和斑状花岗岩样品均引自(吴涛等,2013) Fig. 6 Zr+Nb+Ce+Y vs. (Na2O+K2O)/CaO (a) and Zr+Nb+Ce+Y vs. FeOT/MgO (b) plots of samples from the Tongpu batholith(after Whalen et al., 1987) Data sources: granodiorite, biotite granite and porphyritic granite (Wu et al., 2013)

由于高分异花岗岩的矿物组成和化学成分都趋近于低共结的花岗岩,从而为鉴定其成因类型(如I型、A型、S型或M型花岗岩)带来了困难(吴福元等,2007a)。已有研究表明,可采用两种方法鉴别高分异花岗岩的成因类型(I型或S型花岗岩),一是寻找与高分异花岗岩相关的偏镁铁质岩石,因为它们较多地保留了母岩浆的矿物学和地球化学特征(吴福元等,2007bWu et al., 2017),钙碱性原始玄武质岩浆在温度较低时,通过角闪石、斜长石和Fe-Ti氧化物在跨度达280℃温度范围内的连续分离结晶,可导致熔体SiO2从53%连续上升到78%(Nandedkar et al., 2014);二是利用一些元素地球化学散点图(如SiO2与P2O5、La、Y和Rb-Th、Rb-Y等),这是因为在准铝质和过铝质岩浆中,这些元素具有不同的地球化学行为(Chappell, 1999; 李献华等, 2007)。本文样品与文献报道的同普岩基花岗闪长岩、黑云母花岗岩和斑状花岗岩一起,在SiO2-P2O5散点图上,共同定义了良好的与I型花岗岩特征类似的负相关趋势(图 7a)。对I型花岗岩而言,随着SiO2的增高,Y作为强不相容元素,其含量将急剧升高,在SiO2-Y图解上形成明显的“翘尾”现象(Chappell, 1999),这一特征与同普岩基中SiO2 > 65%样品所定义的趋势相符(图 7b)。同普岩基富硅花岗岩所定义的I型花岗岩趋势,在Rb-Y图解上也得到了很好体现(图 7c)。因为富Y矿物不会在准铝质I型岩浆演化的早期阶段结晶,从而引起分异的I型花岗岩的Y含量高,并与Rb含量呈正相关关系(李献华等,2007朱弟成等,2009)。因此,同普岩基边缘相的高硅花岗岩样品,很可能来源于I型花岗质母岩浆的结晶分异,属于高分异I型花岗岩。

图 7 同普岩基样品的哈克图解 Fig. 7 Harker plots of samples from the Tongpu batholith
4.3 同普岩基高分异花岗岩的成因

一般认为,I型花岗岩主要来源于壳内变质中基性火成岩的部分熔融(Roberts et al., 1993程顺波等,2016)。本文同普岩基边缘相高分异花岗岩负的锆石εHf(t)(-9.9~-5.8) 和古老的地壳模式年龄(1.8~2.0Ga),与同普岩基斑状花岗岩和花岗闪长岩的全岩εNd(t)(-10.3~-8.6) 和二阶段地壳模式年龄(1.6~1.5Ga)(吴涛等,2013)接近,表明这些高分异花岗岩的母岩浆主要来源于古老下地壳的部分熔融。由于这些富硅花岗岩具有与同普岩基中石英闪长岩(εNd(t)=1.37) 不同的同位素组成,表明它们不太可能由闪长质岩石分离结晶而来。考虑到这些边缘相的高分异岩石具有与同普花岗闪长岩类似的同位素组成,指示可能来自共同源区,因此如果假定花岗闪长岩中分异程度最低的样品(JD-35H)可用来代表本文样品的母岩浆,那么定量模拟表明,这些高分异岩石的Ba/Sr比值和Ba、Sr丰度,可由母岩浆经历10%~15%斜长石分离结晶作用来解释(图 8a, b)。由于Sr、Eu对斜长石是相容元素,Nb、Ta和Ti在富钛矿物相(如钛铁矿和/或金红石)中是相容元素,P在磷灰石中是相容元素,因此同普岩基中高硅花岗岩强烈亏损Sr、Eu、Nb、P和Ti等特征(图 5)指示了这些矿物相的分离。

图 8 同普岩基样品的Sr-Ba/Sr (a)和Sr-Ba (b)图解 分配系数引自Rollinson (1993), Green and Pearson (1986), Mahood and Hildreth (1983), Yurimoto et al. (1990) Fig. 8 Ba/Sr vs. Sr (a) and Ba-Sr (b) plots of samples from the Tongpu batholith Partition coefficients are from Rollinson (1993), Green and Pearson (1986), Mahood and Hildreth (1983), Yurimoto et al. (1990)
5 结论

(1) 通过预剥蚀,解决了高硅花岗岩高U或富U锆石的定年难题,并在同普岩基高硅花岗岩中获得了260±1Ma的高质量锆石U-Pb年代学数据。

(2) 同普岩基高硅花岗岩为高分异I型花岗岩,来源于古老地壳物质的部分熔融。

(3) 同普岩基高分异花岗岩可来源于同期花岗闪长质母岩浆经历斜长石、磷灰石和富钛矿物相的分离结晶而成。

致谢 感谢宋婉婧、田园、刘颖和金鹭在野外工作和室内测试中的帮助!感谢审稿老师的评审意见。
参考文献
[] Andersen T. 2002. Correction of common lead in U-Pb analyses that do not report 204Pb. Chemical Geology, 192(1-2): 59–79. DOI:10.1016/S0009-2541(02)00195-X
[] 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
[] Cheng SB, Fu JM, Ma LY, Lu YY, Kou XH, Zhang LG, Huang HL. 2016. Origin of the Yuechengling Caledonian Granitic Batholith, northeastern Guangxi:Constraint from zircon U-Pb geochronology, geochemistry and Nd-Hf isotopes. Geotectonica et Metallogenia, 40(4): 853–872.
[] De la Roche H, Leterrier J, Grande Claude P, Marchal M. 1980. A classification of volcanic and plutonic rocks using R1-R2 diagrams and major element analyses:Its relationship and current nomenclature. Chemical Geology, 29: 183–210. DOI:10.1016/0009-2541(80)90020-0
[] Deng J, Yang LQ, Wang CM. 2011. Research advances of superimposed orogenesis and metallogenesis in the Sanjiang Tethys. Acta Petrologica Sinica, 27(9): 2501–2509.
[] Deng J, Ge LS, Yang LQ. 2013. Tectonic dynamic system and compound orogeny:Additionally discussing the temporal-spatial evolution of Sanjiang orogeny, southwest China. Acta Petrologica Sinica, 29(4): 1099–1114.
[] Deng J, Wang QF, Li GJ, Li CS, Wang CM. 2014a. Tethys tectonic evolution and its bearing on the distribution of important mineral deposits in the Sanjiang region, SW China. Gondwana Research, 26(2): 419–437. DOI:10.1016/j.gr.2013.08.002
[] Deng J, Wang QF, Li G, Santosh M. 2014b. Cenozoic tectono-magmatic and metallogenic processes in the Sanjiang region, southwestern China. Earth-Science Reviews, 138: 268–299. DOI:10.1016/j.earscirev.2014.05.015
[] Deng J, Wang QF, Li GJ, Hou ZQ, Jiang CZ, Danyushevsky L. 2015. Geology and genesis of the giant Beiya porphyry-skarn gold deposit, northwestern Yangtze Block, China. Ore Geology Reviews, 70: 457–485. DOI:10.1016/j.oregeorev.2015.02.015
[] Deng J, Wang QF and Li GJ. 2017. Tectonic evolution, superimposed orogeny, and composite metallogenic system in China. Gondwana Research, doi:10.1016/j.gr.2017.02.005
[] Dong SY, Tang JX, Zhong KH, Gu XX. 2006. Metallogenic subsystems of Jomda tectonic belt in eastern Tibet Plateau. Mineral Deposits(Suppl.1): 4.3–406.
[] Dong SY, Zhong KH, Tang JX, Gu XX, Feng J. 2010. Intracontinental rift evolution and mineralization of the Jiangda tectonic belt on the eastern margin of the Tibetan Plateau, SW China. Journal of Chengdu University of Technology (Science & Technology Edition), 37(4): 475–480.
[] Ewing RC, Meldrum A, Wang L, Weber WJ, Corrales LR. 2003. Radiation effects in zircon. Reviews in Mineralogy and Geochemistry, 53(1): 387–425. DOI:10.2113/0530387
[] Farnan I, Balan E, Pickard CJ, Mauri F. 2003. The effect of radiation damage on local structure in the crystalline fraction of ZrSiO4:Investigating the 29Si NMR response to pressure in zircon and reidite. American Mineralogist, 88(11-12): 1663–1667. DOI:10.2138/am-2003-11-1205
[] Green TH, Pearson NJ. 1986. Rare-earth element partitioning between titanite and coexisting silicate liquid at high pressure and temperature. Chemical Geology, 55: 105–119. DOI:10.1016/0009-2541(86)90131-2
[] He LQ. 1998. Tectonic setting and evolutionary pattern of the Jinshajiang orogenic belt. Geoscience, 12(2): 185–191.
[] Hou ZQ, Zaw K, Pan GT, Mo XX, Xu Q, Hu YZ, Li XC. 2007. Sanjiang tethyan metallogenesis in S. W. China:Tectonic setting, metallogenic epochs and deposit types. Ore Geology Reviews, 31(1-4): 48–87.
[] Krogh TE. 1982. Improved accuracy of U-Pb zircon ages by the creation of more concordant systems using an air abrasion technique. Geochimica et Cosmochimica Acta, 46(4): 637–649. DOI:10.1016/0016-7037(82)90165-X
[] Lee CTA, Morton DM. 2015. High silica granites:Terminal porosity and crystal settling in shallow magma chambers. Earth & Planetary Science Letters, 409: 23–31.
[] Li QL. 2016. "High-U effect" during SIMS zircon U-Pb dating. Bulletin of Mineralogy, Petrology and Geochemistry, 35(3): 405–412.
[] Li XH, Li WX, Li ZX. 2007. On the genetic classification and tectonic implications of the Early Yanshanian granitoids in the Nanling Range, South China. Chinese Science Bulletin, 52(14): 1873–1885. DOI:10.1007/s11434-007-0259-0
[] Li XZ, Du DX, Wang YZ. 1998. Basin-mountain conversion and its mineralization:A study of Qamdo-Simao Basin and Jinsha River-Ailao Mountain. Sedimentary Geology and Tethyan Geology, 22(A01): 5–20.
[] Liu YS, Zong KQ, Kelemen PB, Gao S. 2008a. Geochemistry and magmatic history of eclogites and ultramafic rocks from the Chinese continental scientific drill hole:Subduction and ultrahigh-pressure metamorphism of lower crustal cumulates. Chemical Geology, 247(1-2): 133–153. DOI:10.1016/j.chemgeo.2007.10.016
[] Liu YS, Hu ZC, Gao S, Günther D, Xu J, Gao CG, Chen HL. 2008b. In situ analysis of major and trace elements of anhydrous minerals by LA-ICP-MS without applying an internal standard. Chemical Geology, 257(1-2): 34–43. DOI:10.1016/j.chemgeo.2008.08.004
[] Liu YS, Hu ZC, Zong KQ, Gao CG, Gao S, Xu J, Chen HH. 2010. Reappraisement and refinement of zircon U-Pb isotope and trace element analyses by LA-ICP-MS. Chinese Science Bulletin, 55(15): 1535–1546. DOI:10.1007/s11434-010-3052-4
[] Liu ZC, Li FY, Zhong KH, Li W. 1995. The rift (depression) and metallogenisis of the western margin of the Yangtze platform and its neighbouring region. Geological Exploration for Non-Ferrous Metals, 4(2): 70–76.
[] Ludwig KR. 2000. Users Manual for Isoplot/EX, Vresion 2:A Geochronological Toolkit Microsoft Excel. Berkeley:Berkeley Geochronology Centre Special Publication, No. 1a:53
[] Luo JN. 1992. Sedimentary Geology and Mineralization of Sanjiang Tethys.:People's Republis of China Ministry of Geology and Mineral Resources Geological Memoirs. Beijing: Geological Publishing House: 1-231.
[] Mahood G, Hildreth W. 1983. Large partition coefficients for trace elements in highsilica rhyolites. Geochimica et Cosmochimica Acta, 47: 11–30. DOI:10.1016/0016-7037(83)90087-X
[] Mo XX, Lu FX, Shen SY. 1993. Sanjiang Tethys Volcanism and Related Mineralization. Beijing:Geological Publishing House: 1–269.
[] Mo XX, Deng JF, Dong FL, Yu XH, Wang Y, Yang S, Yang WG. 2001. Volcanic petrotectonic assemblages in Sanjiang orogenic belt, SW China and implication for tectonics. Geological Journal of China University, 7(2): 21–138.
[] Nandedkar RH, Ulmer P, Müntener O. 2014. Fractional crystallization of primitive, hydrous arc magmas:An experimental study at 0. 7Gpa. Contributions to Mineralogy and Petrology, 167(6): 1015. DOI:10.1007/s00410-014-1015-5
[] Patiño Douce AE, Harris N. 1998. Experimental constraints on Himalayan anatexis. Journal of Petrology, 39(4): 689–710. DOI:10.1093/petroj/39.4.689
[] Peng YM, Pan GT, Luo JN. 2000. Arc-basin framework and geologic-tectonic evolution in middle-northern sector of the Jinsha-Nujiang-Lancang rivers area. Acta Geologica Sichuan, 20(3): 176–182.
[] Pitcher WS. 1997. The Nature and Origin of Granite. 2nd Edition. Dordrecht, Netherlands:Springer
[] Rickwood PC. 1989. Boundary lines within petrologic diagrams which use oxides of major and minor elements. Lithos, 22: 247–263. DOI:10.1016/0024-4937(89)90028-5
[] Roberts MP, Clemens JD. 1993. Origin of high-potassium, talc-alkaline, I-type granitoids. Geology, 21(9): 825–828. DOI:10.1130/0091-7613(1993)021<0825:OOHPTA>2.3.CO;2
[] Rollinson HR.1993.Using Geochemical Data:Evaluation, Presentation, Interpretation. Longman Scientific & Technical, London, 1-352
[] Scaillet B, Pichavamt M, Roux J. 1995. Experimental crystallization of leucogranite magmas. Journal of Petrology, 36(3): 663–705. DOI:10.1093/petrology/36.3.663
[] Sun SS, McDonough WF. 1989. Chemical and isotope systematics of oceanic basalts:Implications for mantle composition and processes. In:Sanders AD and Norry MJ (eds.). Magmatism in Ocean Basins. Geological Society Publication, 42(1): 313–345.
[] Wang BD, Wang LQ, Chen JL, Yin FG, Wang DB, Zhang WP, Chen LK, Liu H. 2014a. Triassic three-stage collision in the Paleo-Tethys:Constraints from magmatism in the Jiangda-Deqen-Weixi continental margin arc, SW China. Gondwana Research, 26(2): 475–491. DOI:10.1016/j.gr.2013.07.023
[] Wang F, Liu FL, Liu PH. 2013. Multiple granitic magma events in Sanjiang complex belt, Yunnan Province:Implications for tectonic evolution. Acta Petrologica Sinica, 29(6): 2141–2160.
[] Wang LQ, Pan GT, Li DM, Xu TR. 2000. The evolution and mineralization of the Jomda-Weixi continental marginal volcanic arc, southwestern China. Sedimentary Geology and Tethyan Geology, 20(2): 1–17.
[] Wang LQ, Hou ZQ, Li DM, Pan GT, Guan SP and Xu TR. 2003. The tectonic implications of the post-collisional volcanic rocks of the Late Triassic in the middle segment of the Mid-Jinshajiang tectonic zone. Proceedings of the Symposium on the Tibetan Plateau and Its Adjacent Area Geological Environment and Resources, 77-78 (in Chinese)
[] Wang XF, Metcalfe I, Jian P, He LQ, Wang CS. 2000. The Jinshajiang suture zone:Tectono-stratigraphic subdivision and revision of age. Science in China (Series D), 43(1): 10–22. DOI:10.1007/BF02877827
[] Wang XL, Coble MA, Valley JW, Shu XJ, Kitajima K, Spicuzza MJ, Sun T. 2014b. Influence of radiation damage on Late Jurassic zircon from southern China:Evidence from in situ measurements of oxygen isotopes, laser Raman, U-Pb ages, and trace elements. Chemical Geology, 389: 122–136. DOI:10.1016/j.chemgeo.2014.09.013
[] Whalen JB, Currie KL, Chappell BW. 1987. A-type granites:Geochemical characteristics, discrimination and petrogenesis. Contributions to Mineralogy and Petrology, 95: 407–419. DOI:10.1007/BF00402202
[] Woodhead JD, Hergt JM. 2005. A preliminary appraisal of seven natural zircon reference materials for in situ Hf isotope determination. Geostandards and Geoanalytical Research, 29(2): 183–195. DOI:10.1111/ggr.2005.29.issue-2
[] Wu FY, Li XY, Yang JH, Zheng YF. 2007a. Discussions on the petrogenesis of granites. Acta Geologcia Sinica, 23(6): 1217–1238.
[] Wu FY, Li XY, Zheng YF, Gao S. 2007b. Lu-Hf isotopic systematics and their applications in petrology. Acta Geologcia Sinica, 23(2): 185–220.
[] Wu FY, Liu ZC, Liu XC, Ji WQ. 2015. Himalayan leucogranite:Petrogenesis and implications to orogenesis and plateau uplift. Acta Geologcia Sinica, 31(1): 1–36. DOI:10.1007/s10114-015-4321-6
[] Wu FY, Liu XC, Ji WQ, Wang JM, Yang L. 2017. Highly fractionated granites:Recognition and research. Science China (Earth Sciences), 60(7): 1201–1219. DOI:10.1007/s11430-016-5139-1
[] Wu T, Xiao L, Ma CQ, Huang W. 2013. The geochronological, geochemical and Sr-Nd isotopic characteristics of Tongpu intrusive complex and its implications. Acta Geologcia Sinica, 29(10): 3567–3580.
[] Yurimoto H, Duke EF, Papike JJ, Shearer CK. 1990. Are discontinuous chondritenormalized REE patterns in pegmatitic granite systems the results of monazite fractionation?. Geochimica et Cosmochimica Acta, 54: 2141–2145. DOI:10.1016/0016-7037(90)90277-R
[] Zhang Q, Pan GQ, Li CD, Jin WJ, Jia XQ. 2007. Does fractional crystallization occur in granitic magma? Some crucial questions on granite study (2). Acta Geologcia Sinica, 23(6): 1239–1251.
[] Zhu DC, Mo XX, Wang LQ, Zhao ZD, Niu YL, Zhou CY, Yang YH. 2009. Petrogenesis of highly fractionated I-type granites in the Zayu area of eastern Gangdese, Tibet:Constraints from zircon U-Pb geochronology, geochemistry and Sr-Nd-Hf isotopes. Science in China (Earth Sciences), 52(9): 1223–1239. DOI:10.1007/s11430-009-0132-x
[] Zhu DC, Zhao ZD, Niu YL, Dilek Y, Hou ZQ, Mo XX. 2013. The origin and pre-Cenozoic evolution of the Tibetan Plateau. Gondwana Research, 23: 1429–1454. DOI:10.1016/j.gr.2012.02.002
[] Zi JW, Cawood PA, Fan WM, Tohver E, Wang YJ, Mccuaig TC. 2013. Late Permian-Triassic magmatic evolution in the Jinshajiang orogenic belt, SW China and implications for orogenic processes following closure of the Paleo-Tethys. American Journal of Science, 313(2): 81–112. DOI:10.2475/02.2013.02
[] 程顺波, 付建明, 马丽艳, 卢友月, 寇晓虎, 张利国, 黄惠兰. 2016. 桂东北越城岭岩体加里东期成岩作用:锆石U-Pb年代学、地球化学和Nd-Hf同位素制约. 大地构造与成矿学, 40(4): 853–872.
[] 邓军, 杨立强, 王长明. 2011. 三江特提斯复合造山与成矿作用研究进展——兼论三江复合造山带时空演化. 岩石学报, 27(9): 2501–2509.
[] 邓军, 葛良胜, 杨立强. 2013. 构造动力体制与复合造山作用. 岩石学报, 29(4): 1099–1114.
[] 董树义, 唐菊兴, 钟康惠, 顾雪祥. 2006. 西藏高原东部江达构造带成矿体系. 矿床地质(增1): 403–406.
[] 董树义, 钟康惠, 唐菊兴, 顾雪祥, 冯军. 2010. 西藏高原东部江达构造带陆内裂谷演化与成矿作用. 成都理工大学学报(自然科学版), 37(4): 475–480.
[] 何龙清. 1998. 金沙江造山带的大地构造环境及演化模式. 现代地质, 12(2): 185–191.
[] 李秋立. 2016. 离子探针锆石U-Pb定年的"高U效应". 矿物岩石地球化学通报, 35(3): 405–412.
[] 李献华, 李武显, 李正祥. 2007. 再论南岭燕山早期花岗岩的成因类型与构造意义. 科学通报, 52(9): 81–99.
[] 李兴振, 杜德勋, 王义昭. 1998. 盆山转换及其成矿作用——以昌都-思茅盆地和金沙江-哀牢山带为例. 沉积与特提斯地质, 22(A01): 5–20.
[] 刘肇昌, 李凡友, 钟康惠, 李伟. 1995. 扬子地台西缘及邻区裂谷(陷)构造与金属成矿. 有色金属矿产与勘查, 4(2): 70–76.
[] 罗建宁. 1992. 中华人民共和国地质矿产部地质专报.三, 岩石、矿物、地球化学.第17号.三江特提斯沉积地质与成矿. 北京: 地质出版社: 1-231.
[] 莫宣学, 路凤香, 沈上越. 1993. 三江特提斯火山作用与成矿. 北京: 地质出版社: 1-269.
[] 莫宣学, 邓晋福, 董方浏, 喻学惠, 王勇, 杨肃, 杨伟光. 2001. 西南三江造山带火山岩-构造组合及其意义. 高校地质学报, 7(2): 21–138.
[] 彭勇民, 潘桂棠, 罗建宁. 2000. 三江中北段弧-盆格架与地质构造演化. 四川地质学报, 20(3): 176–182.
[] 王舫, 刘福来, 刘平华. 2013. 云南"三江"变质杂岩带多期花岗质岩浆事件及其构造意义. 岩石学报, 29(6): 2141–2160.
[] 王立全, 潘桂棠, 李定谋, 须同瑞. 2000. 江达-维西陆缘火山弧的形成演化及成矿作用. 沉积于特提斯地质, 20(2): 1–17.
[] 王立全, 侯增谦, 李定谋, 潘桂棠, 管士平, 须同瑞. 2003. 金沙江构造带中段晚三叠世碰撞后火山岩及其构造意义. 见: 青藏高原及邻区地质与资源环境学术讨论会摘要汇编, 77-78 http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDJ200304001117.htm
[] 汪啸风, MetcalfeI, 简平, 何龙清, 王传尚. 1999. 金沙江缝合带构造地层划分及时代厘定. 中国科学(D辑), 29(4): 289–297.
[] 吴福元, 李献华, 杨进辉, 郑永飞. 2007a. 花岗岩成因研究的若干问题. 岩石学报, 23(6): 1217–1238.
[] 吴福元, 李献华, 郑永飞, 高山. 2007b. Lu-Hf同位素体系及其岩石学应用. 岩石学报, 23(2): 185–220.
[] 吴福元, 刘志超, 刘小驰, 纪伟强. 2015. 喜马拉雅淡色花岗岩. 岩石学报, 31(1): 1–36.
[] 吴涛, 肖龙, 马昌前, 黄婉. 2013. 藏东同普杂岩体年代学、地球化学、Sr-Nd同位素特征及大地构造意义. 岩石学报, 29(10): 3567–3580.
[] 张旗, 潘国强, 李承东, 金惟俊, 贾秀琴. 2007. 花岗岩结晶分离作用问题——关于花岗岩研究的思考之二. 岩石学报, 23(6): 1239–1251.
[] 朱弟成, 莫宣学, 王立全, 赵志丹, 牛耀龄, 周长勇, 杨岳衡. 2009. 西藏冈底斯东部察隅高分异I型花岗岩的成因锆石U-Pb年代学、地球化学和Sr-Nd-Hf同位素约束. 中国科学(地球科学), 39(7): 833–848.