岩石学报  2017, Vol. 33 Issue (12): 3729-3740   PDF    
喜马拉雅造山带片麻岩中石榴子石的多期生长
高利娥1 , 曾令森1 , 赵令浩2 , 胡明月2 , 王莉3 , 高家昊1 , 尚振1     
1. 中国地质科学院地质研究所, 北京 100037;
2. 中国地质科学院国家地质实验测试中心, 北京 100037;
3. 河南省有色金属地质勘查总院, 郑州 450052
摘要:石榴子石是高级变质作用的重要矿物之一,能够保留成分环带和不同变质阶段的矿物或流体包裹体,为示踪寄主岩石经历的变质演化历史提供重要信息。本文对采自雅拉香波片麻岩穹窿内高级变质岩中石榴子石进行了详细的包裹体成分、主量元素环带和微量元素特征的研究,揭示出石榴子石黑云母片麻岩至少记录了五期岩浆或变质热事件。第Ⅰ期石榴子石为来源于源区的岩浆型石榴子石。第Ⅱ期、第Ⅲ期和第Ⅳ期石榴子石为变质型石榴子石,但不同期次变质作用的温压条件和生长介质、矿物组合不同。第Ⅴ期石榴子石为转熔型石榴子石,是黑云母脱水熔融形成,记录了喜马拉雅造山过程早期加厚地壳条件下的深熔作用。喜马拉雅造山带变质岩中石榴子石具有复杂的成因机制和演化历史,在应用石榴子石进行变质作用研究时,需要仔细甄别,否则会得到错误的结论。
关键词: 喜马拉雅造山带     片麻岩     石榴子石     多阶段生长    
Multiphase growth of garnet in gneisses from the Himalayan orogenic belt
GAO LiE1, ZENG LingSen1, ZHAO LingHao2, HU MingYue2, WANG Li3, GAO JiaHao1, SHANG Zhen1     
1. Institute of Geology, Chinese Academy of Geological Sciences, Beijing 100037, China;
2. National Research Center for Geoanalysis, Chinese Academy of Geological Sciences, Beijing 100037, China;
3. Henan Nonferrous Metals Geological Exploration Institute, Zhengzhou 450052, China
Abstract: Garnet is one of the important minerals in the metamorphic rocks, can well preserved the major and trace element zonings and mineral or fluid inclusions formed during different metamorphic processes. It provides significant information to trace metamorphic evolution of hosted rocks. Base on carefully checking composition assemblages of inclusions, geochemistry feature of major and trace element of garnets in high grade metamorphic rocks from the Yardoi gneiss dome, we observed at least five generations of garnet growth. Garnet-Ⅰ is magmatic origin and derived from source rock. Garnet-Ⅱ, Ⅲ, Ⅳ are metamorphic origin, however, the minerals assemblages and temperature-pressure condition during three periods of metamorphism are different. Garnet-Ⅴ is peritectic origin and formed during dehydration of biotite occurred in the thickened crustal conditions in the Himalaya time. The evolution process of garnets in the Himalayan high grade metamorphic rocks are complicated; therefore, we need to carefully check the origins of garnets when discussing the metamorphic processes of hosted rocks.
Key words: The Himalayan orogenic belt     Gneiss     Garnet     Multiphase growth    

自然岩石包含三种成因类型的石榴子石:岩浆型、变质型和转熔型(夏琼霞等, 2013; 焦淑娟等, 2013)。不同类型石榴子石常常包含不同矿物包裹体,呈现不同主量和微量元素地球化学特征(高利娥等, 2012)。在中高级变质岩中,石榴子石是重要造岩矿物之一,易保存,为示踪石榴子石寄主岩石经历的变质演化历史提供重要记录(夏琼霞等, 2013),是反演造山带(如:大别-苏鲁,喜马拉雅)构造演化历史的重要探针(Lee et al., 2000; Liu et al., 2006; Guillot et al., 2008; Zhou et al., 2011; Xu et al., 2013; Xia et al., 2012, 2016)。在变质作用的研究中,不管是利用传统温压计,还是热门的相平衡模拟,了解石榴子石的组分特征都起着关键作用。但在一些变质花岗岩中,存在岩浆成因石榴子石的残留核或转熔成因石榴子石的增生边(Habler et al., 2007; Xia et al., 2016; 夏琼霞等, 2013)。与变质成因石榴子石相比,岩浆成因或转熔成因石榴子石具有明显不同的元素成分环带,如岩浆成因石榴石具有较高的锰铝榴石含量,主量元素呈“振荡环带”,Eu为明显的负异常(高利娥等, 2012; 夏琼霞等, 2013)。因此,正确甄别石榴子石的成因机制对于变质岩的研究非常重要,否则简单应用石榴子石的元素地球化学特征来反演寄主岩石的变质演化过程,可能导致错误的结论。

喜马拉雅造山带是世界上最年轻且仍在活动的陆-陆碰撞造山带,较完整地记录了与大陆俯冲-碰撞-折返相关的变质作用,是研究造山带演化历史的天然实验室。位于造山带核部的高喜马拉雅结晶岩系记录了印度-欧亚碰撞过程中高级变质作用,是揭示喜马拉雅造山带演化历史的关键(England and Thompson, 1984; Guillot, 1999; Lee et al., 2000)。自20世纪80年代以来,高喜马拉雅结晶岩系变质作用研究,尤其是倒转变质作用,一直是青藏高原研究的热点(许荣华等, 1986; 任留东和陈炳蔚, 2004; Larson et al., 2013; Wang et al., 2015),并提出了不同类型的造山模型,如挤出模式(Burchfiel and Royden, 1985; Chemenda et al., 2000)、隧道流模型(Beaumont et al., 2001, 2004),但对造山演化历史并没有形成统一的认识。石榴子石作为喜马拉雅高级变质岩中重要的矿物,携带了关于演化历史的关键信息。石榴子石的化学成分环带与温度、压力以及与之平衡的共生矿物组合密切相关。但主量元素易受到后期热扩散的影响,发生重平衡作用,而微量元素环带对岩石形成的化学动力学和矿物组合的变化更为敏感, 可以很好地指示与石榴子石的生长或重结晶相关的矿物化学反应(Konrad-Schmolke et al., 2008)。另外, 石榴子石内的各种包裹体,是揭示石榴子石生长时矿物组合的直接证据。所以可以通过对喜马拉雅变泥质岩中石榴子石的晶格常数(金成伟和王英兰, 1981)、包裹体的成分、主量元素特征(Harris et al., 2004; 李旺超等, 2015)和微量元素特征(高利娥等, 2012)进行综合详细研究,反演喜马拉雅造山带变质演化历史。

出露于北喜马拉雅片麻岩穹窿核部的高级变质岩,经历了和高喜马拉雅变质岩相似的变质和深熔作用(Lee et al., 2004; Zhang et al., 2004; Lee and Whitehouse, 2007; Zeng et al., 2011; Gao et al., 2012a),被认为是高喜马拉雅结晶岩系的北延(Harris and Massey, 1994)。本文对采自雅拉香波片麻岩穹窿的高级变质岩中石榴子石进行了详细的包裹体成分,主量元素环带和微量元素特征的对比研究,来揭示石榴子石生长时的物理化学条件和矿物组合,从而识别石榴子石多阶段生长的成因机制,进而制约寄主岩石演化的P-T-t轨迹。

1 地质背景

喜马拉雅造山带呈E-W向弧形展布(图 1a),自北向南由四个构造单元组成,依次为特提斯喜马拉雅带(内部发育北喜马拉雅片麻岩穹窿,NHGD)、高喜马拉雅带结晶岩系(HHCS)、低喜马拉雅岩系(LHS)和次喜马拉雅岩系(SHS)。它们之间的界限分别为藏南拆离系(STDS)、主中央逆冲断层(MCT)、主边界逆冲断层(MBT)(图 1a)。雅拉香波穹窿位于北喜马拉雅穹窿带的最东端,自核部向边部依次由高级变质岩系、中级变质岩系和特提斯沉积岩系及侵入其中的花岗岩体组成,各岩石单元之间为韧性或脆韧性拆离断层环绕穹隆分布(图 1b)。锆石U-Pb年龄和地球化学特征显示本地区发育三期花岗岩:(1)~44Ma高Sr/Y比值二云母花岗岩,淡色花岗岩和淡色花岗玢岩(Zeng et al., 2015);(2)~35Ma淡色花岗岩(Zeng et al., 2009; 高利娥等, 2009);(3)17~20Ma高Sr/Y比值花岗斑岩和淡色花岗岩(Zeng et al., 2011)。除了这些岩浆岩外,在打拉南侧,~273Ma辉绿岩体侵入到页岩和石英砂岩中(曾令森等, 2012)。高级变质岩由眼球状花岗片麻岩、石榴子石黑云母花岗片麻岩,含石榴子石的泥质片麻岩、石榴角闪岩、角闪岩等组成。石榴角闪岩和花岗质片麻岩记录了~45Ma的变质作用(Zeng et al., 2011Gao et al., 2012a; Ding et al., 2016b),同时花岗质片麻岩的原岩形成于~518.4Ma(Gao et al., 2012a)。T0392为石榴子石黑云母片麻岩,主要由石英、斜长石、钾长石、黑云母、白云母、石榴子石,以及少的磷灰石、蓝晶石、十字石、榍石、电气石和锆石组成(图 2)。岩石中含有大量的黑云母和白云母(图 2a, b),原岩应该是沉积岩,T0392为副变质岩。石榴子石的形态、结构和所含包裹体的成分相差较大(图 3图 4),可能反映了这些片麻岩经历了多期叠加变质作用。为了详细全面地解译片麻岩的演化历史,我们对此岩体进行了系统的采样(T0392-1、T0392-2、T0392-3、T0392-4)。在这些片麻岩中,石榴子石大致有三类:(1)含有大量的石英包裹体;(2)非常干净,无包裹体;和(3)含有大量定向排列的钛铁矿包裹体。我们从T0392系列样品中挑选了两粒结构完整,粒度较大的石榴子石进行研究,来揭示石榴子石的生长环带和成分环带之间的关系,并且甄别不同环带的成因机制。

图 1 藏南喜马拉雅造山带地质简图(a, 据Zeng et al., 2009修改)和雅拉香波穹隆地质简图,示样品T0392位置(b) 数据来自Zeng et al., 2011, 2015; 高利娥等, 2012. YTS-雅鲁藏布缝合带; STDS-藏南拆离系; MCT-主中央逆冲推覆带; MBT-主边界逆冲推覆带; LH-小喜马拉雅岩系 Fig. 1 Simplified geologic map of the Himalayan orogenic belt, southern Tibet (a, after Zeng et al., 2009) and simplified geologic map of the Yardoi area showing the Yardoi gneiss dome (YGD), as well as the locations of sample T0392 Age source data from Zeng et al., 2011, 2015; Gao et al., 2012. YTS-Yarlung-Tsangpo suture; STDS-Southern Tibet Detachment System; MCT-Main Central Thrust; MBT-Main Boundary Thrust; LH-Lower Himalayan Crystalline Sequence

图 2 石榴子石黑云母片麻岩T0392的显微照片 Bt-黑云母;Gt-石榴子石; Mus-白云母; Pl-斜长石; Qtz-石英 Fig. 2 Photomicrographs showing the texture and mineral assemblage of the garnet-biotite gneiss T0392 Bt-biotite; Gt-garnet; Mus-muscovite; Pl-plagioclase; Qtz-quartz

图 3 片麻岩T0392-1中多期石榴子石的背散射图像 Bt-黑云母;IIm-钛铁矿;Qtz-石英;Tour-电气石 Fig. 3 Backscatter images (BSE) show the texture and composition of inclusions for multiphase growth of garnet in the gneiss T0392-1 Bt-biotite; IIm-ilmenite; Qtz-quartz; Tour-tourmaline

图 4 片麻岩中T0392-4多期石榴子石的背散射图像 Ky-蓝晶石;IIm-钛铁矿;Qtz-石英;St-十字石 Fig. 4 Backscatter images (BSE) show the texture and composition of inclusions for multiphase growth of garnet in the gneiss T0392-4 Ky-kyanite; IIm-ilmenite; Qtz-quartz; St-staurolite
2 分析方法

石榴子石主量元素的成分测试在北京大学地球与空间科学学院造山带与成矿作用重点实验室的JEOL-JXA-8800型电子探针仪上进行,加速电压为15kV; 电子束1×10-8A,电子束斑1μm。石榴子石单矿物电子探针分析结果见电子版附表 1

附表 1 片麻岩T0392中石榴子石的主量元素地球化学特征(wt%) AppendixTable1 Chemical composition of garnets in gneiss T0392(wt%)

扫描电镜分析在中国地质科学院地质研究所大陆动力学实验室进行, 所使用的仪器型号为日本电子JEOLJSM-5610LV,工作时的加速电压为20kV, 聚焦距离为20mm, X光能谱仪为英国牛津公司INCA系统。用纯Co金属优化能谱系统, 用背散射电子图像(BSE)进行矿物观察和EDS进行矿物化学成分分析,BSE图像可以很好地显示石榴子石的内部结构特征。石榴子石成分面扫描图见图 5

图 5 片麻岩T0392-1 (a)和T0392-4 (b)中多期石榴子石的Si、Mg、Fe、Mn、Ca元素的X光成分图 Fig. 5 X-ray mapping of Si, Mg, Fe, Mn and Ca element for multiphase growth of garnets in the gneiss T0392-1 (a) and T0392-4 (b)

石榴子石微区原位微量元素分析在国土资源部国家地质测试分析研究中心进行。所用仪器为Finnigan ELEMENT 2 ICP-MS。激光剥蚀系统为新浪UP-213nm激光器。分析采用的激光剥蚀孔径为44μm, 激光脉冲为10Hz, 能量为90mJ。详细的分析方法和流程同胡明月等(2008)。数据处理过程中选择Ca作为内标元素, 国际标准玻璃NIST610用于外标校正。绝大部分微量元素分析的准确度优于10%, Be、P、Ni、Zn、Ga和Yb的准确度优于20%。石榴子石单矿物LA-ICP-MS分析结果见电子版附表 2

附表 2 片麻岩T0392中石榴子石的微量元素(×10-6)地球化学特征 AppendixTable2 trace elements data (×10-6) of garnet in gneiss T0392

为了查明石榴子石黑云母片麻岩记录的变质作用年代学特征,从样品T0392中分选出锆石,经过挑选、制靶和抛光,然后进行阴极发光和扫描电子显微镜成像观察,揭示锆石的内部结构,最后选点进行了SHRIMP锆石U-Pb地质年代学测试。阴极发光成像观察在北京离子探针中心进行。在中国地质科学院地质研究所,利用扫描电子显微镜进行了BSE图像和锆石内部包裹体的成分测试。通过阴极发光和BSE图像来查明锆石内部生长层的分布和结构,选取测试点。锆石U-Pb同位素定年测试在北京离子探针中心进行,所用仪器为高分辨率高灵敏度离子探针SHRIMP Ⅱ。分析时所用标样为TEMORA锆石,每测定3个未知点,插入1次标样,以便及时校正,保障测试精度,测试结果见表 1

表 1 片麻岩T0392的锆石SHRIMP定年数据 Table 1 SHRIMP zircon U-Pb data for gneiss T0392
3 数据及结果 3.1 包裹体的成分

从显微照片和BSE可以看出,单个石榴子石结构比较复杂,从核部到边部,包裹体的成分不同(图 2图 3图 4),包裹体的方向斜切外侧片麻理方向,外侧包裹体斜切内侧包裹体。在仔细对比石榴子石的形态,包裹体矿物组合,包裹关系和X光成分扫描图的基础上,结合石榴石主量元素和微量元素含量,从核部到边部把石榴子石分成五期。

第Ⅰ期石榴子石含小颗粒的包裹体,主要为圆形-椭圆形的石英、黑云母、白云母和钛铁矿和电气石。电气石包裹体,由于硬度大而粒小,在抛光过程中大部分已经剥落,只剩下个别颗粒和大量的空洞。第Ⅱ期石榴子石含有大量的包裹体,呈筛状变晶结构,包裹体主要由石英和少量云母和钛铁矿组成。石英形态不规则,部分为港湾状,粒度显著增大。这期石榴子石内的包裹体与外侧片麻理或斜交或垂直(图 3a),表明该期石榴子石与片麻理不属于同一期次。第Ⅲ期和第Ⅳ期石榴子石很干净,不含有包裹体。第Ⅴ期石榴子石内部含多组定向排列的钛铁矿和金红石,以及黑云母、白云母、石英、蓝晶石、十字石、锆石和磷灰石(图 4b, c),这期包裹体斜切第Ⅱ期包裹体和外侧片麻理(图 4a)。

3.2 主量元素特征

为了查明不同期次石榴子石的成因,选取了粒径为5~6mm的2粒石榴子石为研究对象。样品T0392-1包含四期石榴子石的生长域(图 3图 5a),穿越这些域总共测试了32个点,其中1~6为第Ⅳ期,7~9为第Ⅲ期,10~17为第Ⅱ期,18~31为第Ⅰ期。样品T0392-4也包含四期石榴子石的生长域(图 4图 5b),穿越这些域总共测试了19个点,其中1~5为第Ⅴ期,6~7为第Ⅲ期,8~12和15~19为第Ⅱ期,13~14为第Ⅰ期。首先对这两颗石榴子石进行成分剖面测试,确定石榴子石中Si、Mg、Fe、Mn和Ca的分布特征(图 5)。从X光成分图像来看:T0392-1中第Ⅳ期石榴子石比前三期石榴子石相对富Ca, 贫Mn, 其它成分没明显变化(图 5a);在样品T0392-4中,第Ⅴ期石榴子石也比第Ⅰ、Ⅱ、Ⅲ期的石榴子石富Ca, 贫Mn(图 5b)。从主量元素测试结果来看,这些石榴子石都相对富FeO,贫MgO、CaO和MnO,但不同期次的石榴子石主量元素具有明显的变化(表 2)。T0392-1中(图 6a),总体来看,从第Ⅰ期到第Ⅳ期XMn逐渐降低,XCa呈锯齿状变化,先降低,后升高又降低,XFe和XMg变化不大,石榴子石具体成分变化见表 2。在第Ⅰ期石榴子石中,这四种端元成分锯齿状变化Fe/(Fe+Mg)比值变化较复杂,总体来看,第Ⅰ期和第Ⅱ期Fe/(Fe+Mg)比值一致,第Ⅲ期和第Ⅳ期Fe/(Fe+Mg)比值一致,从第Ⅰ期和第Ⅱ期到第Ⅲ期和第Ⅳ期Fe/(Fe+Mg)比值略升高,第Ⅰ期、第Ⅲ期和第Ⅳ期边部Fe/(Fe+Mg)比值略升高。在样品T0392-4中(图 6b),总体来看,从第Ⅰ期到第Ⅴ期XMn和Fe/(Fe+Mg)比值逐渐降低,XCa、XFe和XMg锯齿状变化, 具体成分变化见表 2

表 2 片麻岩T0392中五期石榴子石的成分差别 Table 2 The difference of composition for five stages growth of garnet in gneiss T0392

图 6 片麻岩T0392-1和T0392-4中多期石榴子石的端元成分剖面图 Fig. 6 Composition cross-section for multiphase growth of garnets in the gneiss T0392-1 and T0392-4
3.3 微量元素特征

为了查明石榴子石在不同期次生长过程中,微量元素的响应,利用LA-ICP-MS,对典型石榴子石进行了微量元素成分剖面的测试(附表 2)。通过对照显微照片和BSE图像,尽量地避开含包裹体区域进行原位LA-ICP-MS测试,但个别点还是受微米级包裹体的影响,因此,对最终数据进行了细致的分析,鉴别那些明显影响石榴子石的稀土元素、Zr、Hf、Y等元素的包裹体,来获取石榴子石的有效微量元素地球化学特征。样品T0392-1包含四期石榴子石的生长域(图 3图 5a),穿越这些域总共测试了33个点,其中1~7为第Ⅳ期,8~10为第Ⅲ期,11~19为第Ⅱ期,20~33为第Ⅰ期。样品T0392-4包含四期石榴子石的生长域(图 4图 5b),穿越这些域总共测试了48个点,其中14~22和36~48为第Ⅴ期,12~13和31~35为第Ⅲ期,1~6、9~11和23~30为第Ⅱ期,7~8为第Ⅰ期。测试结果表明:与主量元素不同,在不同期次石榴子石中微量元素变化明显。

石榴子石是Zr、Hf、Y和HREE的重要载体,不同期次石榴子石在这些元素含量上表现了极大的差异性。稀土元素特征来看,在T0392-1中(图 7a),第Ⅰ期石榴子石重稀土平坦,即(Gd/Yb)N=0.5~5.5;第Ⅱ期石榴子石重稀土也显示平坦的特征,即(Gd/Yb)N=0.4~2.3;第Ⅲ期石榴子石重稀土富集,即(Gd/Yb)N < 0.1;第Ⅳ期石榴子石重稀土亏损,即(Gd/Yb)N > 11.1。在T0392-4中(图 7b),第Ⅰ期石榴子石重稀土平坦,即(Gd/Yb)N=0.7~1.0;第Ⅱ期石榴子石重稀土也平坦,即(Gd/Yb)N=0.1~3.5;第Ⅲ期石榴子石重稀土富集,即(Gd/Yb)N < 0.1;第Ⅴ期石榴子石重稀土平坦,即(Gd/Yb)N=0.3~5.0。上面数据表明这两粒石榴子石同一期显示了相同的稀土元素特征。另外,在两颗石榴子石中五期区域内,Y/Yb比值显示了和(Gd/Yb)N相似的变化,(Gd/Yb)N和Y/Yb比值为正相关(图 8a)。第Ⅳ期Y/Yb比值最高,Y/Yb > 58.8,第Ⅲ期最低,Y/Yb < 8.5,其它三期次,Y/Yb=4.5~37.7。Eu/Eu*在不同期次石榴子石中也有区别(图 8b),第Ⅴ期石榴子石Eu/Eu*高度变化,从弱的负异常到正异常,Eu/Eu*从0.7到3.8,第Ⅰ期石榴子石Eu为明显的负异常,Eu/Eu* < 0.4,其他期次石榴子石一致,显示弱负异常的特征,Eu/Eu*=0.3~0.8。另外,Zr和Hf的含量在不同期次石榴子石中不同(图 8c),第Ⅰ期和第Ⅱ期石榴子石Zr和Hf含量较高且高度变化,Zr和Hf呈正相关关系,Zr分别为7.44×10-6~155.8×10-6和1.66×10-6~541×10-6,Hf分别为0.09×10-6~3.94×10-6和0.01×10-6~15.14×10-6。其它期次石榴子石具有较低的Zr(<7.23×10-6)和Hf含量(<0.28×10-6),Zr和Hf无相关性。

图 7 片麻岩T0392-1和T0392-4中多期石榴子石的球粒陨石标准化稀土元素配分图(标准化值据Sun and McDonough, 1989) Fig. 7 Chondrite-normalized REE distribution patterns for multiphase growth of garnets in the gneiss T0392-1 and T0392-4(normalized values after Sun and McDonough, 1989)

图 8 片麻岩T0392中多期石榴子石的微量元素差异 Fig. 8 The difference of trace elements for multiphase growth of garnets in the gneiss T0392
3.4 锆石年代学特征

样品T0392中大部分锆石为半自形-他形、长柱状,粒度在100~200μm之间,长宽比一般为2:1。锆石磨圆度较好,表明T0392的原岩为沉积岩,经过了一定距离的搬运。大部分锆石具有核-幔-边结构(图 9a),即白色含有包裹体的核部,具韵律生长环带的干净幔部和较窄(~30μm)的暗灰色边部。个别锆石结构不完整。边部锆石207Pb/206Pb年龄为115~539Ma。这些年龄远离谐和线,显示均一化的灰白色或微弱环带,可能是受到后期喜马拉雅造山事件的影响。具有明显韵律环带区域的207Pb/206Pb年龄为1311~1394Ma,5个点的加权平均年龄为1364±31Ma (图 9c), 记录了一期岩浆事件的年龄。微弱环带或均一化的锆石域207Pb/206Pb年龄跨度很大,在805~1546Ma之间(图 9b),相对集中于993±36Ma和1153±15Ma (图 9c), 与格林威尔造山作用的时间相近(张泽明等, 2008; 王誉桦等, 2014)。以上的年龄数据表明:石榴子石黑云母片麻岩T0392中碎屑锆石来源于1364Ma的岩浆事件和格林威尔期的变质事件,并受到喜马拉雅造山事件的影响。

图 9 片麻岩T0392中锆石的CL图像(a)和SHRIMP U-Pb定年谐和图(b) Fig. 9 Cathodoluminescence (CL) images (a) showing the texture, spot and respective age of SHRIMP zircon U-Pb dating and U-Pb concordia diagram (b) for the gneiss T0392
4 讨论 4.1 石榴石的多阶段生长

石榴子石是高级变质岩中常见的矿物,具有难熔和较低溶解度的特点,通常可以保存复杂的成分环带,以及矿物或流体包裹体,为解释石榴子石寄主岩石经历的变质演化历史提供了重要信息。已有研究发现,当岩石经历角闪岩相和麻粒岩相的变质作用时, 温度超过650℃, 并且持续时间足够长时, 由于Mn在石榴子石中的扩散速率较Fe、Mg和Ca高,石榴子石的Mn环带可能经历扩散再平衡过程,出现逐渐平坦的分布特征(Carlson,2002; Gao et al., 2012b)。但微量元素具有较大的离子半径, 热扩散作用对它们的影响不大, 常可以保存很好的成分环带。所以, 石榴子石内微量元素环带对岩石形成的化学动力学和矿物组合的变化更为敏感, 可以较好地反演石榴子石的生长历史(Konrad-Schmolke et al., 2008)。另外,石榴子石内的各种包裹体,包括矿物、熔体或流体包裹体,是揭示石榴子石生长时的物理化学以及矿物组合的直接证据,可以反映石榴子石的多期生长历史(Cutts et al., 2010)。综上所述,可以通过系统分析石榴子石的主量元素、微量元素以及矿物包裹体组合特征,确定石榴子石的生长阶段,揭示寄主岩石的演化历史。

上述数据分析已揭示出石榴子石黑云母片麻岩T0392中石榴子石比较复杂,从核部到边部石榴子石具有以下特征:1)完全不同的包裹体成分;2)后期包裹体切割前期包裹体;3)不同的主量元素和微量元素含量,这些现象表明石榴子石经历了多阶段的生长。从主量元素来看,第Ⅰ期到第Ⅴ期,石榴子石的Spes逐渐降低(图 6),这两粒石榴子石具有保存完好的正生长环带(Spear, 1993),表明此片麻岩没有经历较长时间的高温变质作用。实验岩石学研究表明,在封闭系统中,岩石内石榴子石的Ca元素含量随着压力升高而升高(Hermann and Green, 2001)。第Ⅲ期、Ⅳ期和Ⅴ期石榴子石中,钙铝榴石含量急剧上升,指示了这三期石榴子石形成于较高压的环境。

自然界存在三种类型的石榴子石,即变质型、岩浆型和转熔型,它们在主量元素和微量元素含量上表现出较明显的差异(高利娥等, 2012; 夏琼霞等, 2013; 焦淑娟等, 2013)。从微量元素来看,第Ⅰ期石榴子石显示Eu为明显的负异常,Eu/Eu* < 0.4(图 7),HREE和Y含量较高(图 7)。同时,这期石榴子石具有较高的锰铝榴石(5~14), 主量元素环带呈“振荡型”(图 6)。这些特征与喜马拉雅淡色花岗岩中岩浆型石榴子石(高利娥等, 2012)或大别-苏鲁地区中温/超高压花岗片麻岩中岩浆型石榴石核部的地球化学特征相似(夏琼霞等, 2013),表明第Ⅰ期为残留的岩浆型石榴子石。在岩浆的结晶过程中,由于:1)斜长石的结晶作用摄取了大量的Eu元素;2)缺乏富Mn元素的矿物;3)熔体中含有较高浓度的稀土元素。因此,从岩浆中结晶形成的石榴子石含有较高的Spes,富集HREE,和Eu为负异常。同时,从岩浆型石榴子石的核部到边部,由于瑞利分馏作用,HREE逐渐降低(高利娥等, 2012)。相同的情况,即石榴子石变斑晶内部残留岩浆成因核部,已被前人报道过(Habler et al., 2007)。由于岩浆成因和变质成因石榴子石具有不同的元素成分,因此,在经历多期变质的花岗质岩石中,识别是否含有岩浆型石榴子石核部对于确定变质岩的P-T-t演化轨迹中非常重要。在花岗质岩石中,岩浆型石榴子石多出现在伟晶岩或高度演化的花岗岩中(Dahlquist et al., 2007; 高利娥等, 2012)。由于赋存岩浆型石榴子石的花岗质岩石通常具有富铝和高A/CNK比值,特殊的稀土配分模式(海鸥型)和高Rb/Sr比值等特征,可以应用这些特征,通过对比含石榴子石花岗质片麻岩的全岩元素地球化学特征,来预测这类片麻岩是否残留岩浆型石榴子石,更准确地反演变质历史。

在较高的温度和压力条件下,黑云母脱水部分熔融作用和角闪石部分熔融作用可以形成转熔型石榴子石(Wyllie and Wolf, 1993; Rapp and Watson, 1995; Patiño Douce and Harris, 1998)。黑云母常见的部分熔融反应有:(1)石英+斜长石+白云母±黑云母→铝硅酸盐±石榴石±钾长石+熔体;(2)黑云母+斜长石+铝硅酸盐+石英→石榴子石+钾长石+熔体+含钛相;(3)石英+斜长石+黑云母→斜方辉石+石榴子石+钾长石+熔体。角闪石常见的部分熔融反应有:(4)角闪石+斜长石+石英+水→石榴子石+单斜辉石+榍石+熔体;(5)角闪石+斜长石+石英→石榴子石+单斜辉石+熔体。第Ⅴ期石榴子石为长条状、内部含多组定向排列的含钛矿物相及少量的石英和锆石(图 4)。这种组合与许多实验得出的黑云母脱水部分熔融反应相吻合(黑云母+斜长石+铝硅酸盐+石英→石榴子石+钾长石+熔体+含钛相, Patiño Douce and Johnston, 1991; Patiño Douce and Harris, 1998)。在温度约为900℃条件下,石榴石、石英和钛铁矿/金红石组合构成典型的含Ti黑云母脱水熔融作用的残留相。从上面的分析可以推测出第Ⅴ期石榴子石可能是黑云母脱水熔融形成的转熔型石榴子石,包裹体定向排列的方向为反应物黑云母的解理方向, 石榴石继承了黑云母的形态。蓝晶石包裹体的出现暗示了第Ⅴ期石榴子石生长压力较高。相同特征的转熔型石榴石在雅拉香波穹窿内同一构造域内被报道过,形成于T > 800℃,P > 8kbar条件下,反映了雅拉香波高级变质岩经历过麻粒岩相变质作用和伴随的部分熔融作用(曾令森等,2008)。黑云母脱水部分熔融过程中,斜长石作为主要的熔融相,分解释放Ca和Eu元素到熔体中,被与熔体平衡共生的转熔型石榴子石摄取,使转熔石榴子石表现出含有较高的钙铝榴石端元和Eu正异常(Eu/Eu*=0.7~3.8, 图 7图 8, Xu et al., 2013)。在高级变质岩中,原有的变质石榴石边部增生新生成的转熔石榴石已被前人报道过(Xia et al., 2016)。由于转熔成因和变质成因石榴子石是寄主岩石不同变质阶段的产物,因此,正确的识别转熔石榴石有助于准确反演P-T-t演化轨迹。

第Ⅱ期、第Ⅲ期和第Ⅳ期石榴子石,锰铝榴石含量较低(Spes < 11,图 6),Eu为弱负异常(Eu/Eu*=0.3~0.8,图 8),为变质型石榴子石的特征(Ding et al., 2016a; 李旺超等, 2015)。包裹体的成分和微量元素含量略不同(图 3图 4图 6图 7),第Ⅱ期石榴子石含有大量的以不规则石英颗粒为主的包裹体和HREE含量中等,第Ⅲ期石榴子石不含包裹体和富集HREE,第Ⅳ期石榴子石不含包裹体和亏损HREE。这三期石榴子石没有包裹高压-超高压的矿物,Ca含量中等,具有特征的Mn元素生长环带,表明形成于中-低压变质作用过程中。HREE含量的差异是由于:1)生长介质的地球化学特征;2)同期矿物组合;或3)外来流体/熔体不同导致的。

4.2 石榴子石黑云母片麻岩的演化过程

锆石U-Pb年代学显示(图 9),石榴子石黑云母片麻岩中碎屑锆石来源于1364±31Ma的岩浆事件,1153±15Ma和993±36Ma的变质事件。第Ⅰ期石榴子石为岩浆型石榴石,第Ⅱ期、第Ⅲ期和第Ⅳ期石榴子石为变质型石榴子石,第Ⅴ期石榴子石为转熔型石榴子石。虽然寄主岩石没有直接记录与喜马拉雅造山作用相关的岩浆或者变质热事件,但是115~539Ma的边部年龄表明受到了喜马拉雅造山事件的影响。第Ⅴ期石榴子石中Zr和Hf的含量较低(Zr<7.16×10-6,Hf<0.28×10-6),且无相关性(图 8c),表明在这期地质热事件中流体的成分不利于锆石的结晶。另外,雅拉香波穹窿内同一构造位置的角闪岩和花岗质片麻岩记录了48.2~44.8Ma的变质作用(图 1b, Gao et al., 2012a; Ding et al., 2016b),表明穹窿核部的高级变质岩经历了喜马拉雅造山作用早期的麻粒岩相变质作用,部分岩石可能由于流体不足,锆石没有得到充分的生长。在46~35Ma期间,喜马拉雅造山带处于地壳缩短加厚阶段,以角闪岩为主和变泥质岩为辅的下地壳物质发生脱水部分熔融作用(Zeng et al., 2011),形成高Sr/Y比值的二云母花岗岩,同时在角闪岩和较富沃的高级变质岩中形成转熔型石榴子石。结合花岗质片麻岩中石榴石的地球化学特征,锆石年代学特征和所处的构造背景,可以推测出第Ⅴ期石榴子石为喜马拉雅造山作用的产物。第Ⅰ期岩浆型石榴石与锆石类似,可能来源于源区。与第Ⅴ期石榴子石相似,第Ⅲ期和第Ⅳ期石榴子石具有较低的Zr和Hf含量,Zr/Hf比值无相关性(图 8c),表明在这两期地质热事件中流体的成分也不利于锆石的结晶,锆石没有参与第Ⅲ期和第Ⅳ期石榴子石的形成过程。所以,第Ⅱ期、第Ⅲ期和第Ⅳ期变质石榴子石是否与喜马拉雅造山作用相关,目前没法推测,还需要进一步对不同期次石榴子石中独居石包裹体进行原位U-Pb定年。

5 结论

高喜马拉雅结晶岩系是喜马拉雅造山带重要的岩石组合,是揭示喜马拉雅造山带演化历史的关键。出露于北喜马拉雅片麻岩穹窿核部的高级变质岩,经历了和高喜马拉雅变质岩相似的变质和深熔作用,被认为是高喜马拉雅结晶岩系的北延。通过对采自雅拉香波片麻岩穹窿内高级变质岩中石榴子石进行详细的包裹体成分,主量元素环带和微量元素特征的对比研究,揭示出石榴子石黑云母片麻岩至少记录了五期岩浆或变质热事件。第Ⅰ期石榴子石具有较高的锰铝榴石,端元成分呈“振荡型环带”,Eu为明显的负异常,HREE和Y含量较高,为岩浆型石榴子石。第Ⅱ期、第Ⅲ期和第Ⅳ期石榴子石中,锰铝榴石含量较低,Eu为弱的负异常,为变质型石榴子石,记录了中低压条件下的变质作用,但不同期次变质作用的温压条件和生长介质,矿物组合不同。第Ⅴ期石榴子石内部含多组定向排列的钛铁矿和金红石、黑云母、白云母、石英、蓝晶石、十字石、锆石和磷灰石包裹体,Eu从弱的负异常到正异常,HREE含量中等,是黑云母脱水熔融形成,为转熔型石榴子石,记录了喜马拉雅造山过程早期加厚地壳条件下的深熔作用。喜马拉雅造山带变质岩中石榴子石具有复杂的结构和多样的成因机制(岩浆型、变质型和转熔型)。在应用石榴子石的元素地球化学特征来解释寄主岩石的变质历史时需要正确甄别石榴子石的成因机制。

致谢 感谢吴春明教授和戚学祥研究员仔细审阅稿件,提出众多建设性修改意见。感谢张建新研究员、张泽明研究员和董昕副研究员在写作过程中给予的有益指导。
参考文献
Beaumont C, Jamieson RA, Nguyen MH and Lee B. 2001. Himalayan tectonics explained by extrusion of a low-viscosity crustal channel coupled to focused surface denudation. Nature, 414(6865): 738-742. DOI:10.1038/414738a
Beaumont C, Jamieson RA, Mai HN and Medvedev S. 2004. Crustal channel flows:1.Numerical models with applications to the tectonics of the Himalayan-Tibetan orogen. Journal of Geophysical Research, 109(B6): B06406. DOI:10.1029/2003JB002809
Burchfiel BC and Royden LH. 1985. North-south extension within the convergent Himalayan region. Geology, 13(10): 679-682. DOI:10.1130/0091-7613(1985)13<679:NEWTCH>2.0.CO;2
Carlson WD. 2002. Scales of disequilibrium and rates of equilibration during metamorphism. American Mineralogist, 87(2-3): 185-204. DOI:10.2138/am-2002-2-301
Chemenda AI, Burg JP and Mattauer M. 2000. Evolutionary model of the Himalaya Tibet system:Geopoem:Based on new modelling, geological and geophysical data. Earth and Planetary Science Letters, 174(3-4): 397-409. DOI:10.1016/S0012-821X(99)00277-0
Cutts KA, Kinny PD, Strachan RA, Hand M, Kelsey DE, Emery M, Frend CRL and Leslie AG. 2010. Three metamorphic events recorded in a single garnet:Integrated phase modelling, in situ LA-ICPMS and SIMS geochronology from the Moine Supergroup, NW Scotland. Journal of Metamorphic Geology, 28(3): 249-267. DOI:10.1111/jmg.2010.28.issue-3
Dahlquist JA, Galindo C, Pankhurst RJ, Rapela CW, Alasino PH, Saavedra J and Fanning CM. 2007. Magmatic evolution of the Peñón Rosado granite:Petrogenesis of garnet-bearing granitoids. Lithos, 95(3-4): 177-207. DOI:10.1016/j.lithos.2006.07.010
Ding HX, Zhang ZM, Hu KM, Dong X, Xiang H and Mu HC. 2016a. P-T-t-D paths of the North Himalayan metamorphic rocks:Implications for the Himalayan orogeny. Tectonophysics, 683: 393-404. DOI:10.1016/j.tecto.2016.06.035
Ding HX, Zhang ZM, Dong X, Tian ZL, Xiang H, Mu HC, Gou ZB, Shui XF, Li WC and Mao LJ. 2016b. Early Eocene (ca. 50Ma) collision of the Indian and Asian continents:Constraints from the North Himalayan metamorphic rocks, southeastern Tibet. Earth and Planetary Science Letters, 435: 64-73.
England PC and Thompson AB. 1984. Pressure-temperature-time paths of regional metamorphism 1.Heat transfer during the evolution of regions of thickened continental crust. Journal of Petrology, 25(4): 894-928. DOI:10.1093/petrology/25.4.894
Gao LE, Zeng LS, Liu J and Xie KJ. 2009. Early Oligocene Na-rich peraluminous leucogranites in the Yardoi gneiss dome, southern Tibet:Formation mechanism and tectonic implications. Acta Petrologica Sinica, 25(9): 2289-2302.
Gao LE, Zeng LS and Xie KJ. 2012a. Eocene high grade metamorphism and crustal anatexis in the North Himalaya Gneiss Domes, Southern Tibet. Chinese Science Bulletin, 56: 639-650.
Gao LE, Zeng LS, Shi WG, Chen ZY, Hu MY and Sun DY. 2012. Two types of garnets in the Cenozoic granites from the Himayalan Orogenic Belt:Geochemical characteristics and implications for crustal anatexis. Acta Petrologica Sinica, 28(9): 2963-2980.
Gao XY, Zheng YF and Chen YX. 2012b. Dehydration melting of ultrahigh-pressure eclogite in the Dabie orogen:Evidence from multiphase solid inclusions in garnet. Journal of Metamorphic Geology, 30(2): 193-212. DOI:10.1111/jmg.2012.30.issue-2
Guillot S. 1999. An overview of the metamorphic evolution in central Nepal. Journal of Asian Earth Sciences, 17(5-6): 713-725. DOI:10.1016/S1367-9120(99)00045-0
Guillot S, Mahéo G, de Sigoyer J, Hattori KH and Pêcher A. 2008. Tethyan and Indian subduction viewed from the Himalayan high-to ultrahigh-pressure metamorphic rocks. Tectonophysics, 451(1-4): 225-241. DOI:10.1016/j.tecto.2007.11.059
Habler G, Thöni M and Miller C. 2007. Major and trace element chemistry and Sm-Nd age correlation of magmatic pegmatite garnet overprinted by eclogite-facies metamorphism. Chemical Geology, 241(1-2): 4-22. DOI:10.1016/j.chemgeo.2007.01.026
Harris N and Massey J. 1994. Decompression and anatexis of Himalayan metapelites. Tectonics, 13(6): 1537-1546. DOI:10.1029/94TC01611
Harris NBW, Caddick M, Kosler J, Goswami S, Vance D and Tindle AG. 2004. The pressure-temperature-time path of migmatites from the Sikkim Himalaya. Journal of Metamorphic Geology, 22(3): 249-264. DOI:10.1111/jmg.2004.22.issue-3
Hermann J and Gree DH. 2001. Experimental constraints on high pressure melting in subducted crust. Earth and Planetary Science Letters, 188(1-2): 149-186. DOI:10.1016/S0012-821X(01)00321-1
Hu MY, He HL, Zhan XC, Fan XT, Wang G and Jia ZR. 2008. Matrix normalization for in-situ multi-element quantitative analysis of zircon in laser ablation-inductively coupled plasma mass spectrometry. Chinese Journal of Analytical Chemistry, 36(7): 947-953.
Jiao SJ, Guo JH and Peng SB. 2013. Petrogenesis of garnet in the Darongshan-Shiwandashan granitic suite of the South China Block and the metamorphism of the granulite enclave. Acta Petrologica Sinica, 29(5): 1740-1758.
Jin CW and Wang YL. 1981. Garnets of crystalline rocks in the Himalayas. Scientia Geologica Sinica, (1): 12-20.
Konrad-Schmolke M, Zack T, O'Brien PJ and Jacob DE. 2008. Combined thermodynamic and rare earth element modelling of garnet growth during subduction:Examples from ultrahigh-pressure eclogite of the Western Gneiss Region, Norway. Earth and Planetary Science Letters, 272(1-2): 488-498. DOI:10.1016/j.epsl.2008.05.018
Larson KP, Gervais F and Kellett DA. 2013. A P-T-t-D discontinuity in east-central Nepal:Implications for the evolution of the Himalayan mid-crust. Lithos, 179: 275-292. DOI:10.1016/j.lithos.2013.08.012
Lee J, Hacker BR, Dinklage WS, Wang Y, Gans P, Calcert A, Wan JL, Chen WJ, Blythe AE and McClelland W. 2000. Evolution of the Kangmar Dome, southern Tibet:Structural, petrologic, and thermochronologic constraints. Tectonics, 19(5): 872-895. DOI:10.1029/1999TC001147
Lee J, Hacker B and Wang Y. 2004. Evolution of North Himalayan gneiss domes:Structural and metamorphic studies in Mabja Dome, southern Tibet. Journal of Structural Geology, 26(12): 2297-2316. DOI:10.1016/j.jsg.2004.02.013
Lee J and Whitehouse MJ. 2007. Onset of mid-crustal extensional flow in southern Tibet:Evidence from U/Pb zircon ages. Geology, 35(1): 45-48. DOI:10.1130/G22842A.1
Li WC, Zhang ZM, Xiang H, Guo ZB and Ding HX. 2015. Metamorphism and anatexis of the Himalayan orogen:Petrology and geochronology of HP pelitic granulites from the Yadong area, Southern Tibet. Acta Petrologica Sinica, 31(5): 1219-1234.
Liu FL, Gerdes A, Liou JG, Xue HM and Liang FH. 2006. SHRIMP U-Pb zircon dating from Sulu-Dabie dolomitic marble, eastern China:Constraints on prograde, ultrahigh-pressure and retrograde metamorphic ages. Journal of Metamorphic Geology, 24(7): 569-589. DOI:10.1111/jmg.2006.24.issue-7
Patiño Douce AE and Harris N. 1998. Experimental constraints on Himalayan anatexis. Journal of Petrology, 39(4): 689-710. DOI:10.1093/petroj/39.4.689
Patiño Douce AE and Johnston AD. 1991. Phase equilibria and melt productivity in the pelitic system:Implications for the origin of peraluminous granitoids and aluminous granulites. Contributions to Mineralogy and Petrology, 107(2): 202-218. DOI:10.1007/BF00310707
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
Ren LD and Chen BW. 2004. The implication of the "silimanite zone" in the high Himalayan crystalline rocks. Advances in Earth Science, 19(5): 715-721.
Spear FS. 1993. Metamorphic Phase Equilibria and Pressure-Temperature-Time Paths. Washington, D.C.:Mineralogical Society of America
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
Wang JM, Rubatto D and Zhang JJ. 2015. Timing of partial melting and cooling across the Greater Himalayan Crystalline Complex (Nyalam, Central Himalaya):In-sequence thrusting and its implications. Journal of Petrology, 56(9): 1677-1702. DOI:10.1093/petrology/egv050
Wang YH, Zeng LS, Gao LE, Zhang LF and Hou KJ. 2014. Labradonian and Greenvillian orogenic events in the Namche Barwa Massif of the Himalayan orogenic belt. Acta Petrologica Sinica, 30(8): 2241-2252.
Wyllie PJ and Wolf MB. 1993. Amphibolite dehydration-melting:Sorting out the solidus. Geological Society, London, Special Publication, 76(1): 405-416. DOI:10.1144/GSL.SP.1993.076.01.20
Xia QX, Zheng YF, Lu XN, Hu ZC and Xu HJ. 2012. Formation of metamorphic and metamorphosed garnets in the low-T/UHP metagranite during continental collision in the Dabie orogen. Lithos, 136-139: 73-92. DOI:10.1016/j.lithos.2011.10.004
Xia QX, Zhou LG and Zheng YF. 2013. Zonation and multiphase growth of garnet in UHP metamorphic rocks of continental subduciton zone. Chinese Science Bulletin, 58(22): 2138-2144.
Xia QX, Wang HZ, Zhou LG, Gao XY, Zheng YF, Van Orman AJ, Xu HJ and Hu ZC. 2016. Growth of metamorphic and peritectic garnets in ultrahigh-pressure metagranite during continental subduction and exhumation in the Dabie orogen. Lithos, 266-267: 158-181. DOI:10.1016/j.lithos.2016.08.043
Xu LJ, Xiao YL, Wu F, Li SG, Simon K and Wörner G. 2013. Anatomy of garnets in a Jurassic granite from the south-eastern margin of the North China Craton:Magma sources and tectonic implications. Journal of Asian Earth Sciences, 78: 198-221. DOI:10.1016/j.jseaes.2012.11.026
Xu RH, Cheng ZL, Gui XT and Wang JW. 1986. A study of the main metamorphic age of Nyalam Group (Tibet). Acta Petrologica Sinica, 2(2): 13-22.
Zeng LS, Liang FH, Xu ZQ and Qi XX. 2008. Metapelites in the Himalayan orogenic belt and their protoliths. Acta Petrologica Sinica, 24(7): 1517-1527.
Zeng LS, Liu J, Gao LE, Xie KJ and Wen L. 2009. Early Oligocene anatexis in the Yardoi gneiss dome, southern Tibet and geological implications. Chinese Science Bulletin, 54(1): 104-112. DOI:10.1007/s11434-008-0362-x
Zeng LS, Gao LE, Xie KJ and Jiang LZ. 2011. Mid-Eocene high Sr/Y granites in the Northern Himalayan Gneiss Domes:Melting thickened lower continental crust. Earth and Planetary Science Letters, 303(3-4): 251-266. DOI:10.1016/j.epsl.2011.01.005
Zeng LS, Gao LE, Hou KJ and Tang SH. 2012. Late Permian mafic magmatism along the Tethyan Himalayan Belt, southern Tibet and tectonic implications. Acta Petrologica Sinica, 28(6): 1731-1740.
Zeng LS, Gao LE, Tang SH, Hou KJ, Guo CL and Hu GY. 2015. Eocene magmatism in the Tethyan Himalaya, southern Tibet. Geological Society, London, Special Publications, 412(1): 287-316. DOI:10.1144/SP412.8
Zhang HF, Harris N, Parrish R, Kelley S, Zhang L, Rogers N, Argles T and King J. 2004. Causes and consequences of protracted melting of the mid-crust exposed in the North Himalayan antiform. Earth and Planetary Science Letters, 228(1-2): 195-212. DOI:10.1016/j.epsl.2004.09.031
Zhang ZM, Wang JL, Zhao GC and Shi C. 2008. Geochronology and Precambrian tectonic evolution of the Namche Barwa complex from the eastern Himalayan syntaxis, Tibet. Acta Petrologica Sinica, 24(7): 1477-1487.
Zhou LG, Xia QX, Zheng YF and Chen RX. 2011. Multistage growth of garnet in ultrahigh-pressure eclogite during continental collision in the Dabie orogen:Constrained by trace elements and U-Pb ages. Lithos, 127(1-2): 101-127. DOI:10.1016/j.lithos.2011.08.010
高利娥, 曾令森, 刘静, 谢克家. 2009. 藏南也拉香波早渐新世富钠过铝质淡色花岗岩的成因机制及其构造动力学意义. 岩石学报, 25(9): 2289-2302.
高利娥, 曾令森, 石卫刚, 陈振宇, 胡明月, 孙东阳. 2012. 喜马拉雅造山带新生代花岗岩中两类石榴石的地球化学特征及其在地壳深熔作用中的意义. 岩石学报, 28(9): 2963-2980.
胡明月, 何红蓼, 詹秀春, 樊兴涛, 王广, 贾泽荣. 2008. 基体归一定量技术在激光烧蚀-等离子体质谱法锆石原位多元素分析中的应用. 分析化学, 36(7): 947-953.
焦淑娟, 郭敬辉, 彭松柏. 2013. 华南大容山-十万大山花岗岩体中石榴石成因以及麻粒岩包体变质作用研究. 岩石学报, 29(5): 1740-1758.
金成伟, 王英兰. 1981. 喜马拉雅结晶岩系中的石榴子石. 地质科学, (1): 12-20.
李旺超, 张泽明, 向华, 苟正彬, 丁慧霞. 2015. 喜马拉雅造山带核部的变质作用与部分熔融:亚东地区高压泥质麻粒岩的岩石学与年代学研究. 岩石学报, 31(5): 1219-1234.
任留东, 陈炳蔚. 2004. 高喜马拉雅变质岩"矽线石带"的地质意义. 地球科学进展, 19(5): 715-721.
王誉桦, 曾令森, 高利娥, 张立飞, 侯可军. 2014. 喜马拉雅造山带东构造结拉布拉多期和格林威尔期造山作用的记录. 岩石学报, 30(8): 2241-2252.
夏琼霞, 周李岗, 郑永飞. 2013. 大陆俯冲带超高压变质岩石中石榴石的环带和多阶段生长. 科学通报, 58(22): 2138-2144.
许荣华, 成忠礼, 桂训唐, 王俊文. 1986. 西藏聂拉木群主变质时代的讨论. 岩石学报, 2(2): 13-22.
曾令森, 梁凤华, 许志琴, 戚学祥. 2008. 喜马拉雅造山带变泥质岩系及其地球化学特征. 岩石学报, 24(7): 1517-1527.
曾令森, 高利娥, 侯可军, 唐索寒. 2012. 藏南特提斯喜马拉雅带晚二叠纪基性岩浆作用及其构造地质意义. 岩石学报, 28(6): 1731-1740.
张泽明, 王金丽, 赵国春, 石超. 2008. 喜马拉雅造山带东构造结南迦巴瓦岩群地质年代学和前寒武纪构造演化. 岩石学报, 24(7): 1477-1487.