岩石学报  2017, Vol. 33 Issue (5): 1420-1444   PDF    
喜马拉雅碰撞造山带新生代地壳深熔作用与淡色花岗岩
曾令森, 高利娥     
中国地质科学院地质研究所, 北京 100037
摘要: 自从印度-欧亚大陆碰撞以来,伴随着构造演化和温度-压力-成分(P-T-X)的变化,喜马拉雅造山带中下地壳变质岩发生不同类型的部分熔融反应,形成性质各异的过铝质花岗岩。这些花岗岩在形成时代、矿物组成、全岩元素和放射性同位素地球化学特征上都表现出巨大的差异性。始新世构造岩浆作用形成高Sr/Y二云母花岗岩和演化程度较高的淡色花岗岩和淡色花岗玢岩,它们具有相似的Sr-Nd同位素组成,是碰撞早期增厚下地壳部分熔融的产物。渐新世淡色花岗岩主要为演化程度较高的淡色花岗岩,可能指示了喜马拉雅造山带的快速剥露作用起始于渐新世。早中新世以来的淡色花岗岩是喜马拉雅造山带淡色花岗岩的主体,是变泥质岩部分熔融的产物,包含两类部分熔融作用——水致白云母部分熔融作用(A类)和白云母脱水熔融作用(B类)。这两类部分熔融作用形成的花岗质熔体在元素和同位素地球化学特征上都表现出明显的差异性,主要受控于两类部分熔融作用过程中主要造岩矿物和副矿物的溶解行为。这些不同期次的地壳深熔作用都伴随着高分异淡色花岗岩,伴随着关键金属元素(Nb、Ta、Sn、Be等)的富集,是未来矿产勘探的重要靶区。新的观测结果表明:在碰撞造山带中,花岗岩岩石学和地球化学性质的变化是深部地壳物质对构造过程响应的结果,是深入理解碰撞造山带深部地壳物理和化学行为的重要岩石探针。
关键词: 淡色花岗岩     高分异花岗岩     地壳深熔作用     碰撞造山带     喜马拉雅造山带    
Cenozoic crustal anatexis and the leucogranites in the Himalayan collisional orogenic belt
ZENG LingSen, GAO LiE     
Institute of Geology, Chinese Academy of Geological Sciences, Beijing 100037, China
Abstract: Since the India-Eurasian continental collision, the Himalayan orogenic belt has experienced major tectonic transitions from earlier crustal compression and thickening to later extension and rapid exhumation, which induced pronounced changes in the pressure-temperature-composition (P-T-X) of high-grade metamorphic rocks. Consequently, the mid to lower crustal rocks has undergone correspondingly different partial melting processes and produced a wide spectrum of melts of leucogranitic compositions. Such granites show substantial differences in ages of crystallization, mineral compositions, whole-rock element as well as radiogenic isotope (Sr and Nd) compositions. The earliest anatexis is represented by the Eocene (~43Ma) high Sr/Y granites from melting of mafic rocks under thickened crustal conditions, followed by the melting of metasedimentary rocks in the Oligocene time, possibly induced by the initiation of rapid exhumation of the Himalayan high grade basements. A majority of leucogranites formed from ~25Ma to ~10Ma were derived either from fluxed melting of muscovite (A-type) or from muscovite dehydration melting (B-type) of metasediments. These two modes of crustal anatexis, possibly from vastly similar source rocks, generated granitic melts with substantial differences in major and trace element as well as Sr isotope geochemistry due to the coupled differences in the melting behavior of the major minerals (muscovite, feldspar) and accessory phases (zircon and monazite) during different modes of crustal anatexis. Each phases of leucogranite production are accompanied by various degrees of differentiation and the formation of highly fractionated leucogranites. Some of such leucogranites are highly enriched in key metal elements (Sn, Nb, Ta, and Be) and thus could be potential targets for future exploration of precious metals. Data summarized in this contribution suggests a strong coupling of granite chemical compositions with the changes in tectonic regimes, which in turn implies that leucogranites, by sorting out their original melt compositions, could serve as a valuable probe to investigate the physical and chemical behavior of deep crustal rocks in collisional orogenic belts worldwide.
Key words: Leucogranite     Highly fractionated granite     Crustal anatexis     Collisional orogenic belt     Himalaya    
1 引言

汇聚板块边缘是大陆地壳形成、活化和分异的重要场所。碰撞造山作用是板块汇聚作用的主要作用之一,可产生足够的热量,引发深部岩石的高级变质和深熔作用(Brown, 2007, 2010, 2013)。碰撞造山带通常具有顺时针P-T-t变质作用轨迹,在相应的地壳温度和压力条件下,变沉积岩(变杂砂岩和变泥质岩)是最容易发生部分熔融的岩石。但在有水或在特殊构造环境下(如裂谷带或快速剥露的高级变质地体),易熔组分(如基性岩和富含云母的花岗质岩石)也可发生地壳深熔作用。这些不同类型岩石的部分熔融都可以通过两种主要深熔作用方式来实现——含水部分熔融(也叫水致部分熔融)和脱水部分熔融。导致含水部分熔融的水的来源包括:(1) 高压变质作用的脱水作用;(2) 岩石网状裂隙中含水流体的汇聚,和(3) 熔体迁移过程中流体循环和结晶作用出溶的水(Thompson, 2001; Brown, 2010, 2013; Weinberg and Hasalová, 2015)。但是,“干”地壳在更高的温度条件下,发生脱水熔融反应,特征熔融产物指示温度达到800℃以上(Brown, 2010, 2013)。

喜马拉雅碰撞造山带,世界上最年轻且仍在活动的陆-陆碰撞型造山带,较完整地记录了与大陆深俯冲-快速折返-伸展作用等地质作用相关的变形、变质和岩浆作用的重要信息(Hodges, 2000; Yin and Harrison, 2000; Nabelek and Liu, 2004; Yin, 2006; Nábělek et al., 2009),是研究碰撞造山带深部物质物理和化学过程的野外基地,是检验各种造山作用理论和造山带演化模型的天然实验室。喜马拉雅造山带经历了早期(ca.50~30Ma)的地壳增厚作用和进变质作用(Vance and Mahar, 1998; Vance and Harris, 1999; Searle et al., 1999; Vannay and Grasemann, 2001; Gao et al., 2012; Kohn, 2014)和后期的快速折返和减压作用(Kohn, 2014)。喜马拉雅造山带发育三种类型的构造:近东西向展布的南向逆冲推覆构造(MCT、MBT和MFT等)和北向藏南拆离系(STDS)和近南北向展布、东西向伸展的藏南裂谷系(STRZ)(图 1a)。自~60Ma碰撞以来,该造山带经历了广泛的地壳水平缩短和垂向加厚作用(Yin, 2006; Aikman et al., 2008; Diedesch et al., 2016),形成了数条延伸上千千米以上的构造带,如大型逆冲断层-主中央逆冲断层(MCT)、主边界逆冲断层(MBT)和主前锋逆冲断层(MFT),和大型伸展构造-藏南拆离系(STDS)、北喜马拉雅片麻岩穹窿(NHGD)和南北向裂谷系(STRZ)(Zhang et al., 2012)。已有大量研究来限定藏南裂谷系E-W向伸展作用的启动时限(Coleman and Hodges, 1995; Blisniuk et al., 2001; Williams et al., 2001; Jolivet et al., 2003; Jessup and Cottle, 2010; Zhang et al., 2012; Mitsuishi et al., 2012),最近的研究表明藏南裂谷系最早启动可能开始于19Ma(Zhang et al., 2012; Mitsuishi et al., 2012)。

图 1 藏南喜马拉雅造山带地质简图(a)和雅拉香波穹隆地质简图(b) YTS-雅鲁藏布缝合带;STDS-藏南拆离系;MCT-主中央逆冲推覆带;MBT-主边界逆冲推覆带;LH-小喜马拉雅岩系.锆石U-Pb年龄数据包括Aikman et al.(2008, 2012b), Zeng et al.(2011, 2012, 2015)和Gao et al. (2012) Fig. 1 Simplified geologic map of the Himalayan orogenic Belt, southern Tibet (a) and simplified geologic map of the Yardoi area showing the Yardoi gneiss dome, Dala and Quedang plutons and the nearly east-west trending subvolcanic porphyritic leucogranite field from Qiaga to Longzi (b) YTS-Yarlung-Tsangpo suture; STDS-Southern Tibet Detachment System; MCT-Main Central Thrust; MBT-Main Boundary Thrust; LH-Lower Himalayan Crystalline Sequence. Zircon U/Pb age data from Aikman et al.(2008, 2012b), Zeng et al.(2011, 2012, 2015) and Gao et al. (2012)

在长度上千千米的喜马拉雅造山带中,广泛发育不同时代和不同地球化学性质的花岗岩类岩石(Le Fort et al., 1987; Harrison et al., 1997, 1999a, b; Prince et al., 2001; Zhang et al., 2004a, 2012; Ding et al., 2005; Aikman et al., 2008, 2012a, b; Searle et al., 2010; King et al., 2011; Zeng et al., 2011, 2012, 2015; Hou et al., 2012; Gao and Zeng, 2014; Wang et al., 2017; 王誉桦等, 2014; 吴福元等, 2015; 高利娥等, 2015, 2016)。花岗岩类岩石包括:(1) 原岩时代为古元古代(ca.1850~1800Ma)(Le Fort, 1981; Ahmad et al., 2000);(2) 原岩时代为新元古代(ca.850~800Ma)(DiPietro and Isachsen, 2001; Richards et al., 2006; Chu et al., 2011; Ahmad et al., 2013; Wang et al., 2017; 许志琴等, 2005);(3) 原岩时代为ca.520~480Ma的花岗岩(Schärer et al., 1986; Miller et al., 2001; Gehrels et al., 2006a, b; Quigley et al., 2008; Cawood et al., 2007; Wang et al., 2012; 许志琴等, 2005; 张泽明等, 2008; 高利娥等, 2012, 2015; 王誉桦等, 2014)和(4) 新生代( < 50Ma)淡色花岗岩。其中新生代淡色花岗岩的形成时间尺度从~45Ma一直到现今,且各主要时代的地球化学性质表现出较大的差异性。这些淡色花岗岩为了解中下地壳物质在碰撞造山过程的物理和化学行为提供了重要的岩石探针。厘定这些淡色花岗岩的产出状态、形成时限、岩石学、地球化学和地球物理性质及其成矿效应等是揭示碰撞造山带构造演化过程中,变质、岩浆和构造作用之间的互馈关系,浅表作用和深部地壳地质作用之间互动关系的关键。

在喜马拉雅造山带,按岩石的演化程度,淡色花岗岩包括二云母花岗岩、白云母花岗岩、含电气石或石榴子石淡色花岗岩和伟晶岩,多数高分异淡色花岗岩含电气石和石榴子石,少数含绿柱石、锡石和铌钽矿等(图 2)。自二十世纪八十年代以来,从不同侧面,大量学者开展了喜马拉雅淡色花岗岩的岩石学和地球化学(Vidal et al., 1982; Debon et al., 1986; Deniel et al., 1987; Harrison et al., 1997, 1999a, b)和地质年代学(Schärer et al., 1986; Harrison et al., 1995; Prince et al., 2001; Aoya et al., 2005; Lee and Whitehouse, 2007)的研究。

图 2 喜马拉雅造山带不同类型淡色花岗岩显微照片 (a)二云母花岗岩;(b)淡色花岗玢岩;(c)含黑云母淡色花岗岩;(d)含电气石淡色花岗岩;(e)含石榴子石淡色花岗岩; (f)含绿柱石淡色花岗岩. Qtz-石英; Pl-斜长石; Mus-白云母; Bi-黑云母; Grt-石榴子石; Tur-电气石; Be-绿柱石; Ap-磷灰石.比例尺为300μm Fig. 2 Microphotographs showing the texture and mineral assemblage (a) two mica granites; (b) porphyritic leucogranite; (c) biotite-bearing leucogranite; (d) tourmaline-bearing leucogranite; (e) garnet-bearing leucogranite; (f) beryl-bearing leucogranite from the Himalayan orogenic belt. Qtz-quartz; Pl-plagioclase; Mus-muscovite; Bi-biotite; Grt-garnet; Tur-tourmaline; Be-beryl; Ap-apatite. Scale bar=300μm

已有的结果表明,随着温压条件及原岩化学成分的变化,变泥质岩的部分熔融常常通过几个不同的部分熔融反应来实现,形成一系列地球化学特征不同的花岗岩熔体(Harris and Inger, 1992; Douce and Harris, 1998; Knesel and Davidson, 2002; Zeng et al., 2005a, b)。由于在部分熔融反应中,熔融相的化学计量上的差异,不同类型部分熔融反应形成的花岗质岩浆在元素和同位素地球化学特征上都可能表现出明显的差异性(Inger and Harris, 1993; Douce and Harris, 1998; Harrison et al., 1999a, b; Knesel and Davison, 2002; Zeng et al., 2005c; Gao and Zeng, 2014; Gao et al., 2017)。因此,在造山带构造演化过程中,随着P-T-X(温度-压力-含水量)条件的变化,地壳物质有可能发生不同类型的部分熔融反应,形成地球化学性质不同的花岗质岩浆,这些花岗质岩浆为探讨造山带深部物质如何响应造山过程提供了良好的岩石探针。

喜马拉雅造山带广泛发育淡色花岗岩,早期研究认为这些花岗岩主要形成于28~10Ma(Harris and Massey, 1994; Harris et al., 1995; Guillot and Le Fort; 1995; Harrison et al., 1997, 1999a, b; Simpson et al., 2000; Zhang et al., 2004b及其中的参考文献),是在喜马拉雅高级岩片快速折返过程中高级变泥质岩白云母脱水部分熔融的产物(Le Fort, 1981; Harrison et al., 1997; Douce and Harris, 1998)。但最近的研究揭示了在印度-欧亚大陆陆陆碰撞造山作用期间,喜马拉雅造山带经历了复杂多样的深熔作用,如角闪岩部分熔融作用(Zeng et al., 2009, 2011, 2015; Liu et al., 2014; 高利娥等, 2009)和变泥质岩水致白云母部分熔融作用(Prince et al., 2001; Zhang et al., 2004a; King et al., 2011; Guo and Wilson, 2012; Zeng et al., 2012)。同时,锆石U/Pb地质年代学研究也揭示上述深熔作用可追溯到始新世(~46Ma)(Ding et al., 2005; Aikman et al., 2008; Larson et al., 2010; Zeng et al., 2009, 2011, 2015; Hou et al., 2012; Liu et al., 2014, 2016; Carosi et al., 2014; Gao et al., 2016; 戚学祥等, 2008; 谢克家等, 2010; 胡古月等, 2011)。在本文中,综合近10年来有关喜马拉雅造山带新生代淡色花岗岩研究的进展,来刻画自50Ma以来,喜马拉雅造山带淡色花岗岩在不同时间段的元素和同位素(Sr和Nd)上的差异性,探讨(1) 淡色花岗岩源区的变化;(2) 地壳深熔作用的类型和(3) 与区域构造演化的关系。

2 主要地壳岩石部分熔融反应类型 2.1 变沉积岩部分熔融反应

在地壳岩石中,变泥质岩和变花岗岩的部分熔融反应类型基本上相似,包括较低温度下的水致部分熔融和较高温度下的云母脱水部分熔融。沉积岩的部分熔融反应包括含水和脱水部分熔融,在水饱和条件下,主要有以下类型部分熔融反应:

(1) 石英+斜长石+流体=熔体

(2) 石英+斜长石+白云母+水=熔体

(3) 石英+斜长石+黑云母+水=石榴子石+堇青石+熔体

在无水条件下,沉积岩的部分熔融主要通过含水矿物(白云母或黑云母)的脱水部分熔融来实现(Vielzeuf and Holloway, 1988; Le Breton and Thompson, 1988; Douce and Johnston, 1991; Vielzeuf and Montel, 1994; Gardien et al., 1995; Douce and Beard, 1995, 1996; Singh and Johannes, 1996; Stevens et al., 1997; Douce and Harris, 1998)。常见的部分熔融反应有:

石英+斜长石+白云母±黑云母=铝硅酸盐±石榴子石±钾长石+熔体

石英+斜长石+黑云母+铝硅酸盐=石榴子石±堇青石±钾长石+流体

石英+斜长石+黑云母=斜方辉石+石榴子石+钾长石+熔体

在上述部分熔融反应中,其中白云母脱水部分熔融(MS1) 和水致白云母部分熔融(MS2) 反应可能是最常见的部分熔融反应,具体为

白云母脱水部分熔融:22Ms+7Pl+8Qtz=25Melt+5Kfs+5Sil+2Bt (MS1)

水致白云母部分熔融:9Ms+15Pl+7Qtz+xH2O=31Melt (MS2)

其中MS=白云母;Pl=斜长石;Qtz=石英;Kfs=钾长石;Sil=矽线石;Bt=黑云母;Melt=熔体。

值得指出的是,在高压(P≥10kbar)条件下,钠长石相对于钙长石稳定,形成熔体应较富集Ca、Sr、Eu及其与Ca相似性质的元素(Douce and Harris, 1998)。在变泥质岩部分熔融过程中,水的存在可以稳定含水相。在高压和低温条件下,由于水的存在,白云母和黑云母趋于稳定,当P=10kbar,T=700~775℃,富钾沉积岩可以通过水致熔融产生富Na2O和CaO的熔体(Douce and Harris, 1998)。在相似的温压条件下,脱水熔融产生的熔体具有明显相反的特征,富K2O,相对亏损CaO和Na2O(Conrad et al., 1988; Douce, 1996; Douce and Harris, 1998)。在南迦巴瓦地区,赋存于变泥质岩的~25Ma淡色体具有与苦堆穹窿~28Ma二云母花岗岩相似的特征(Zhang et al., 2004b; King et al., 2011)。与较年轻淡色花岗岩相比,这些花岗岩含蓝晶石,具有较高的CaO和较低的结晶温度,代表了高压条件下变泥质岩水致部分熔融的产物(Zhang et al., 2004a; King et al., 2011; Zeng et al., 2012)。

2.2 基性岩部分熔融反应

在地壳温压条件下,基性岩石(玄武岩、辉长岩或角闪岩)的部分熔融通过角闪石脱水部分熔融来实现(Wolf and Wyllie, 1991, 1993, 1994; Wyllie and Wolf, 1993; Rapp et al., 1995; López and Castro, 2001),常见的部分熔融反应有:

角闪石+斜长石+石英+水=石榴子石+单斜辉石+榍石+熔体

角闪石+斜长石+石英=石榴子石+单斜辉石+熔体

上述两种反应都会形成由石榴子石和辉石组成的残留相,具有与榴辉岩类似的矿物组成和物理性质,密度比通常的上地幔岩石高0.2~0.3g/cm3。花岗岩岩基根部或增厚地壳的部分熔融主要是角闪岩脱水部分熔融,都可形成石榴辉石岩或榴辉岩(取决于角闪岩中角闪石和斜长石的组分)。富钙角闪岩部分熔融形成石榴辉石岩;富钠角闪岩部分熔融可能形成严格意义上的榴辉岩残留体。在高压和含水条件下,黑云母部分熔融可在温度低于800℃的条件下发生,形成石榴子石+角闪石+熔体的组合(Gardien et al., 2000), 也形成高密度的石榴角闪岩残留体;这种高密度残留相导致重力失稳,是引发造山带或大陆岛弧大规模拆沉作用的主要因素(Ducea and Saleeby, 1996; Saleeby et al., 2003; 曾令森等, 2006)。值得指出的是,地壳岩石常常含两种或两种以上的含水矿物,如角闪岩中,常出现角闪石和黑云母共存的情况,这种岩石的部分熔融行为较复杂,取决于温度、压力、是否有自由水及含水矿物的化学组成,有可能同时发生角闪石和黑云母脱水部分熔融的情况,形成具有混合特征的花岗质熔体。

3 地壳岩石部分熔融过程中矿物反应

在地壳深熔作用中,与造岩矿物直接参与部分熔融不同,副矿物(锆石、磷灰石、独居石、榍石等)通常通过溶解作用,间接调整花岗质熔体的微量元素和同位素地球化学特征。要正确解译花岗质熔体的关键微量元素(如REE)和放射性同位素(Sr、Nd和Hf)的地球化学特征,需要了解不同副矿物在花岗质熔体中的溶解行为。

3.1 磷灰石在岩浆中的溶解度

磷是磷灰石的必要结构组分之一,与磷灰石保持平衡关系的岩浆的磷的含量受磷灰石的溶解度的缓冲。实验结果表明了在熔体中,磷灰石的饱和度是熔体温度和成分的函数,在水含量为0~10%范围内,与熔体的含水量无关(Harrison and Watson, 1984; Pichavant et al., 1992; Wolf and London, 1994; London et al., 1999)。磷灰石的溶解度可以用磷灰石和熔体之间磷的配分系数来描述,磷的配分系数越小,越多的磷灰石要溶解到熔体中,直到熔体达到磷的饱和浓度。Harrison and Watson (1984)给出了下述公式来描述磷在熔体中的溶解度:

其中SiO2是熔体中SiO2的重量百分比,T是绝对温度。该方程适用于SiO2含量为45%到75%,含水量小于10%的熔体,及在地壳压力条件下。对过铝酸性岩浆,Pichavant et al. (1992)进一步的实验研究对Harrison-Watson模型进行了进一步的修正,认为:

其中P2O5HW按公式(1) 来计算,P2O5Per按下述公式来计算:

其中,A/CNK是Al/(2Ca+K+Na)的摩尔比值,PMR、HW和Per分别代表Pichavant-Montel-Richard、Harrison-Watson和过铝质。修正的模型可以解释过铝质花岗岩的较高的P2O5

其中CApP=43%。

3.2 独居石溶解模型

轻稀土元素(从La到Gd)及磷都是独居石基本组构组分(essential structural constituent)。实验研究表明:(1) 岩浆中的磷的浓度控制着含磷矿物的稳定性,如果磷足够多,独居石在控制岩浆的LREE的浓度中起重要作用;(2) 岩浆中的LREE饱和浓度是熔融体温度、成分及含水量的函数。在变泥质岩和高铝花岗岩中,独居石是轻稀土元素和生热元素U和Th的主要赋存矿物之一。独居石的稳定性影响着地壳深熔作用导致的地壳中包括轻稀土、U和Th在内的微量元素的重新分布。因此了解独居石在地壳深熔作用中的地球化学行为是认识某些花岗岩地球化学性质的关键。

在部分熔融过程中,独居石是否参与熔融作用(即独居石溶解到岩浆中)或在岩浆的演化后期是否结晶独居石取决于熔体的磷及LREE的浓度。独居石在贫钙的酸性岩浆中的溶解度可以通过以下公式来计算(Montel, 1993)

其中 Na、K、Li、Ca、Al和Si按原子百分比,H2O是重量百分比来计算。T是绝对温度。REE包含La到Gd,但排除Eu。

在部分熔融稳定低于800℃的条件下,取决于源区的初始LREE的浓度,深熔作用产生过铝花岗质岩浆,相应熔融区残留部分则富集独居石。但是当熔融温度接近于850℃时,地壳深熔作用导致的下地壳富集独居石不太可能,因为在高温情况下产生的岩浆具有非常高的独居石饱和溶解度(Rapp and Watson, 1986)。在地壳深熔作用中,独居石的稳定程度是缓冲岩浆P2O5的浓度的主要因素之一。在磷灰石和独居石共存的情况下,根据Harrison-Waston的磷灰石溶解模型,磷灰石的溶解将抑制独居石在岩浆中的溶解程度,在致使岩浆的P2O5含量是否达到饱和上,独居石比磷灰石所起的作用弱些。

3.3 锆石溶解模型

锆元素是锆石的必要结构组分。前人对不同成分岩石部分熔融所产生的熔体的锆溶解度实验结果表明在不同成分的岩浆(SiO2含量为54.2%~69.5%,M=1.01~2.09) 中,锆的饱和程度取决于岩浆温度和岩浆成分(Watson and Harrison, 1983; Harrison and Watson, 1984; Watson and Harrison, 1984)。锆石在岩浆中的溶解度满足下述方程:

其中DZrZircon/Melt是Zr在锆石和岩浆中浓度的比值,T绝对温度,和M是阳离子的含量比值,由公式(Na+K+2Ca)/(Al×Si)来计算。值得指出的是,在这些实验中,过碱性岩浆具有异常高的Zr溶解度(Watson, 1979; Watson and Harrison, 1983),如Watson and Harrison (1983)实验中的第21次实验,熔体具有9.8%的Na2O+K2O,为过碱性,Zr含量为860×10-6,熔体初始温度为860℃。这种现象的原因是,Zr4+与熔体中不与Al结合的自由碱离子(K+或Na+)形成络合物,Zr在过碱性岩浆中的溶解度增加。但对于非过碱性岩浆,Zr4+和碱离子不形成络合物。

最近的系统研究表明了上述公式可以用来估算岩浆的温度,即锆石饱和温度计(Miller et al., 2003)。他们的研究结果表明:如果岩浆中Zr不饱和,锆石饱和温度计所给出的锆石饱和温度(TZr)为初始岩浆温度的下限;而如果岩浆中Zr达到饱和状态,所给出的锆石饱和温度为初始岩浆温度的上限。如果已知某一岩浆的主要成分及锆的含量(Zrmelt),可以根据下述公式来计算锆石饱和温度(Miller et al., 2003):

3.4 含钛矿物的溶解模型

Ryerson and Watson (1987)实验研究了压力为8~30kbar,温度在1000~1300℃,在含水-CO2饱和及不含水条件下,金红石在硅酸盐岩浆中的溶解度。硅酸盐岩浆的TiO2饱和溶解度由下述公式给出:

其中CmeltCrutile分别是熔体和金红石中TiO2的浓度;T为温度,单位为K;P为压力,单位为kbar;FM是岩浆的成分参数,由下述公式计算:

其中金属阳离子为摩尔数。压力的大小轻微影响TiO2在硅酸盐岩浆中的溶解度,但含水量对其几乎没有影响(Green and Pearson, 1987; Ryerson and Watson, 1987)。

4 喜马拉雅造山带地壳深熔作用的表现方式和产物

地壳深熔作用表现在不同的尺度下,(1) 野外露头尺度(米-几十米级);(2) 手标本尺度(厘米级);(3) 显微尺度(亚毫米级)。高喜马拉雅结晶岩系由石榴黑云母片麻岩、眼球状片麻岩、变泥质岩(含变沙质岩块体或层)、石榴斜长角闪岩、含透辉石大理岩及石英岩组成。其中角闪岩、泥质岩和黑云母片麻岩都不同程度地混合岩化,包含厚度不等(亚毫米至几个厘米)的浅色体。浅色体类型包括:(1) 含石榴子石和角闪石的浅色花岗岩,以近平行于面理或呈梭形透镜体的形式,存在于斜长角闪岩中;(2) 含石榴子石淡色体;(3) 含电气石淡色体。同构造花岗岩脉呈两种产状,平行于面理和垂直/斜交面理。

喜马拉雅碰撞造山带发育多种典型的新生代S-型花岗岩。这些花岗岩以透镜体、岩脉或大型岩体的形式分布于北喜马拉雅片麻岩穹窿和高喜马拉雅结晶岩系中(图 1)。众多喜马拉雅淡色花岗岩体都存在多幕淡色花岗岩,在短短的几个百万年内,形成年龄和地球化学特征存在明显差异的淡色花岗岩。如玛纳斯鲁(Manaslu)岩体包括至少两阶段淡色花岗岩,早期形成于22~24Ma和后期形成于19.0~19.6Ma(Harrison et al., 1999b);佩古错岩体包含28.2±0.5Ma电气石淡色花岗岩和19.8±0.5Ma二云母花岗岩(Gao et al., 2013);在夏如穹隆中,除了原岩形成于古生代的淡色花岗质片麻岩外(Liu et al., 2016),还包括分别形成于~35Ma和~28Ma的含电气石淡色花岗岩(Liu et al., 2016; Gao et al., 2016)。因此,要深入了解各个淡色花岗岩体的形成过程和成因机理,需要细致确定它们的形成年代和全岩地球化学特征。

5 喜马拉雅造山带各主要阶段淡色花岗岩的地球化学特征

在利用花岗岩全岩元素地球化学组成来反演地壳深熔作用的源区和部分熔融反应类型的过程中,需要甄别它们是否代表较原始的熔体。在花岗质熔体的形成、抽取、汇聚、融合和最终的侵位过程中,影响花岗岩最终组分的因素包括:源区特征、部分熔融反应类型、岩浆形成后的后期岩浆作用(包括分离结晶作用和同化混染作用)。在分析大量数据的基础上,考虑它们的矿物组合、元素和同位素(Sr和Nd)组成及其副矿物在花岗质熔体的溶解行为(Miller and Stoddard, 1981; Miller and Mittlefehldt, 1982, 1984; Mittlefehldt and Miller, 1983; Miller, 1985; Shearer and Papike, 1987; Watt and Harley, 1993; Bea, 1996; Sylvester, 1998; Zeng et al., 2005a, b; Claiborne et al., 2006; Clarke, 2007; Stevens et al., 2007; Ballouard et al., 2016; Liu et al., 2016; Gao et al., 2017; 吴福元等, 2015; 郭春丽等, 2017),可利用以下判别标志(表 1)来甄别较原始岩浆, 其中较原始的花岗质岩浆是指(1) 对于某类原岩,全岩花岗岩的矿物组成和主量元素地球化学特征较接近实验岩石学产生的熔体组分的那些花岗岩;或(2) 对某一岩体,全岩地球化学(主量元素和微量元素)组成演化程度最低的样品。从始新世到中新世,按其形成时代、元素和同位素地球化学特征,喜马拉雅淡色花岗岩可分为始新世(包括二云母花岗岩、淡色花岗岩和淡色花岗玢岩),渐新世淡色花岗岩和中新世(25~10Ma)淡色花岗岩。其中始新世淡色花岗岩包括二云母花岗岩、淡色花岗岩和淡色花岗玢岩(图 2a, b),渐新世主要为分异程度较高的淡色花岗岩,中新世包括A、B、C和F四类淡色花岗岩,其中A和B类代表较原始的岩浆。各类花岗岩的总体元素地球化学特征见表 2

表 1 高分异或含堆晶花岗岩的典型地球化学特征 Table 1 Major geochemical characteristics for the highly fractionated granites or granites subjected to accumulation of plagioclase or K-feldspar

表 2 喜马拉雅造山带新生代淡色花岗岩全岩地球化学特征(主量元素:wt%; 微量元素:×10-6 Table 2 Summary of the element geochemistry of the Himalayan Cenozoic leucogranites (major oxides in wt%; trace elements in ×10-6)
5.1 始新世花岗岩

沿喜马拉雅造山带,始新世(ca.41~46Ma,锆石U-Pb年龄)地壳深熔作用表现为(1) 形成较原始的二云母花岗岩(TMG)(Zeng et al., 2011, 2015; Hou et al., 2012; Liu et al., 2014; 戚学祥等, 2008; 高利娥等, 2009; 谢克家等, 2010)及其演化程度较高的淡色花岗岩(LG)和淡色花岗玢岩(LGP)(Zeng et al., 2015; 胡古月等, 2011);(2) 石榴子石包裹的纳米花岗岩(Carosi et al., 2014)。另外除了雅拉香波-隆子地区(图 1b)和然巴穹窿外(图 1),在萨嘎的纽库地区也发育年龄相近的二云母花岗岩(Ding et al., 2005),但未见全岩地球化学数据。

5.1.1 始新世二云母花岗岩

始新世二云母花岗岩具有以下特征:(1) 较高的SiO2(68.3%~72.5%)和Al2O3(15.1%~17.0%);(2) 较低的FeO( < 2.0%)和MgO( < 1.5%)(图 3);(3) Na2O/K2O(0.94~1.42) 和A/CNK(1.04~1.23) 比值都较高(图 4),表明它们是富钠过铝质花岗岩。

图 3 喜马拉雅淡色花岗岩TiO2 (a)、Al2O3(b)、FeO (c)、MgO (d)、CaO (e)和P2O5(f)与SiO2的关系图 Fig. 3 Diagrams showing the relationships of selected major oxides of TiO2 (a), Al2O3(b), FeO (c), MgO (d), CaO (e) and P2O5(f) plotted against SiO2 for the Himalayan Cenozoic leucogranites Data from Zhang et al., 2004a; Inger and Harris, 1993; Searle et al., 1999; Aoya et al., 2005; King et al., 2011; Zeng et al., 2011, 2015; Liu et al., 2014, 2016; Gao and Zeng, 2014; Gao et al., 2016, 2017. Data sources in Fig. 4-Fig. 8 and Fig. 10 are same as in this figure

图 4 喜马拉雅新生代淡色花岗岩Na2O/K2O-SiO2 (a)和P2O5-A/CNK (b)关系图 Fig. 4 Na2O/K2O vs. SiO2(a) and P2O5 vs. A/CNK (b) diagrams for the Himalayan Cenozoic leucogranites

二云母花岗岩都富集LREE((La/Yb)N>20.0,高达56.8),亏损HREE且HREE(从Ho到Lu)较平直,(Ho/Yb)N比值为0.88到1.50,具弱或无Eu异常(Eu/Eu*=0.79~1.10)(图 5a)。在微量元素地球化学特征上,二云母花岗岩具有(1) 较高的Sr(252×10-6~416×10-6)但较低的Y( < 8.5×10-6),和较高的Sr/Y比(32.8~96.4)(图 6);(2) 与同期的淡色花岗岩和淡色花岗玢岩相比,在相似的Rb含量下,这些二云母花岗岩的Sr含量高2倍以上(图 7);和(3) 较低的Na/Ta比值( < 13.4) 和变化范围较大的Zr/Hf(28.6~36.5) 比值(图 8)。二云母花岗岩的87Sr/86Sr(t)和εNd(t)值都较低,为0.7112~0.7193和-14.9~-9.2(Zeng et al., 2011, 2015; Liu et al., 2014)(图 9)。

图 5 喜马拉雅新生代淡色花岗岩球粒陨石标准化稀土元素分布模式图(标准化值据Sun and McDonough, 1989) Fig. 5 Chondrite-normalized rare earth element distribution patterns for the Himalayan Cenozoic leucogranites (normalization values after Sun and McDonough, 1989)

图 6 喜马拉雅新生代淡色花岗岩Sr/Y-Y (a)和Sr/Y-Sr (b)关系图 Fig. 6 Sr/Y vs. Y (a) and Sr/Y vs. Sr (b) diagrams for the Himalayan Cenozoic leucogranites

图 7 喜马拉雅新生代淡色花岗岩Rb/Sr-Ba (a)、Rb/Sr-Sr (b)、Sr-Rb (c)和Sr-Eu/Eu* (d)关系图 A=水致白云母部分熔融;B=白云母脱水部分熔融 Fig. 7 Rb/Sr vs. Ba (a), Rb/Sr vs. Sr (b), Sr vs. Rb (c) and Sr vs. Eu/Eu*(d) diagrams for the Himalayan Cenozoic leucogranites A=fluxed muscovite melting; B=muscovite dehydration melting

图 8 喜马拉雅新生代淡色花岗岩Nb/Ta-Nb (a)和Zr/Hf-Zr (b)关系图 A=正常岩浆分异作用;B=熔体结构变化导致的Zr-Hf分馏 Fig. 8 Nb/Ta vs. Nb (a) and Zr/Hf vs. Zr (b) diagrams for the Himalayan Cenozoic leucogranites A=normal magmatic differentiation; B=melt's structure changes induced fractionation of Zr over Hf

图 9 喜马拉雅新生代淡色花岗岩εNd-87Sr/86Sr (a)、εNd-Nd (b)、87Sr/86Sr-Sr (c)和εNd-Sm/Nd (d)关系图 基性岩数据来自Zeng et al., 2011; Liu et al., 2014;变泥质岩数据来源于Zhang et al., 2004a; King et al., 2011; Zeng et al., 2009, 2011, 2012, 其它数据来源与图 3相同 Fig. 9 εNd vs. 87Sr/86Sr (a), εNd vs. Nd (b), 87Sr/86Sr vs. Sr (c) and εNd vs. Sm/Nd (d) diagrams for the Himalayan Cenozoic leucogranites Field for the mafic rocks from Zeng et al., 2011; Liu et al., 2014; Metapelite field from Zhang et al., 2004a; King et al., 2011; Zeng et al., 2009, 2011, 2012; the other data sources same as in the Fig. 3
5.1.2 始新世淡色花岗岩

在主量元素组成上,该套淡色花岗岩的特征如下:(1) SiO2和Al2O3含量较高,分别为72.8%~75.5%和14.4%~17.0%;(2) TiO2( < 0.1%)、FeO( < 0.8%)、MgO( < 0.3%)和MnO( < 0.1%)含量较低(图 3);(3) 具有较高的Na/K比值(Na2O/K2O=1.10~2.29) 和A/CNK(1.05~1.63)(图 4),显示过铝质的特征,表明这些淡色花岗岩为过铝质富钠花岗岩。

在稀土元素上,该套花岗岩表现出较复杂的稀土配分模式,具有以下特征:(1) LREE总体上表现出较一致,但从MREE到HREE,含量变化较大;(2) 在11件样品中,2件样品具有强分馏的稀土元素分配模式,富集LREE,亏损HREE(图 5);(3) 与二云母花岗岩不同的是,除了具有强烈负Eu异常(Eu/Eu*=0.15~0.55) 的样品外,个别样品具有强烈的正Eu异常或无Eu异常(Eu/Eu*=1.03~1.46)(图 5b)。和二云母花岗岩相比,该套淡色花岗岩的Rb(92.5×10-6~206×10-6),Sr( < 172×10-6)、Ba( < 370×10-6)(图 6图 7)和Zr( < 45×10-6)(图 8)都较低,Na/Ta(3.7~9.7) 和Zr/Hf(9.4~20.0) 比值较低(图 8),明显低于球粒陨石,但Rb/Sr明显较高(0.98~2.68)。淡色花岗岩的87Sr/86Sr(t)和εNd(t)值都较低,为0.7112~0.7153和-12.5~-9.8(Zeng et al., 2011, 2015),与二云母花岗岩的类似(图 9)。

5.1.3 始新世淡色花岗玢岩

到目前为止,以花岗玢岩状态产出的淡色花岗质岩石仅在隆子县扎西岗-隆子县城一带出露(图 1b)。在主量元素组成上,该套岩石的特征如下:(1) SiO2(73.9%~77.3%)和Al2O3(13.8%~16.0%)含量较高;(2) TiO2( < 0.05%)、FeO( < 0.8%)、MgO( < 0.3%)和MnO( < 0.03%)含量较低(图 3);(3) 具有较高的Na/K比值(Na2O/K2O=1.09~2.12) 和A/CNK(1.14~1.71)(图 4),表明这些淡色花岗岩为过铝质富钠花岗岩。

在稀土元素上,与二云母和淡色花岗岩相比,该套花岗岩具有以下特征:(1) LREE含量较低,HREE相似;(2) 多数样品表现出四分组REE分配模式(胡古月等, 2011);(3) 明显的负Eu异常(Eu/Eu*=0.45~0.65)(图 5c),与高分异淡色花岗岩的特征类似。和二云母花岗岩相比,该套淡色花岗岩的Rb(134×10-6~275×10-6),Sr( < 92×10-6)、Zr( < 30×10-6)和Ba( < 135×10-6)都较低(图 6图 7);(4) Na/Ta(4.8~7.0) 和Zr/Hf(14.9~18.6) 比值较低(图 8),明显低于球粒陨石,但Rb/Sr(0.98~2.68) 明显较高和Sr/Y(8.1~27.4) 比较低(图 6图 7)。与同期的二云母花岗岩和淡色花岗岩相似,87Sr/86Sr(t)和εNd(t)值都较低,为0.7112~0.7185和-13.8~-13.2(图 9)(Zeng et al., 2011, 2015)。

5.2 渐新世淡色花岗岩

在喜马拉雅造山带,已报道的具有年龄和地球化学限定的渐新世淡色花岗岩都是含石榴子石/电气石淡色花岗岩,包括夏如穹窿(Gao et al., 2016),苦堆穹窿(Zhang et al., 2004b; King et al., 2011)、康巴穹窿(Liu et al., 2016)和佩古错穹窿(Aoya et al., 2005; Gao et al., 2013, 2016; 张进江等, 2011; 王晓先等, 2015, 2016)。另外昌果穹窿也发育年龄~31Ma的淡色花岗岩(Larson et al., 2010),但未报道地球化学数据。

在主量元素组成上,该期淡色花岗岩表现为:(1) SiO2(72.3%~78.0%),Al2O3(12.9%~15.1%)和FeO(0.1%~2.1%)含量变化都较大;(2) TiO2( < 0.3%)、MgO( < 0.6%)和MnO( < 0.05%)含量较低(图 3);(3) Na/K比值(Na2O/K2O=0.37~0.91) 较低和A/CNK(1.00~1.14) 较高(图 4),为过铝质富钾淡色花岗岩。

在稀土元素上,该期花岗岩具有以下特征:(1) LREE含量较低,LREE和HREE含量相似,呈现出“海鸥状”的稀土分配模式;(3) 明显的负Eu异常(Eu/Eu*=0.08~0.81)(图 5d),与高分异花岗岩的特征相似(Gao et al., 2016; Liu et al., 2016; 吴福元等, 2015)。

该期花岗岩具有(1) 较高的Rb(281×10-6~712×10-6),但Sr(10.8×10-6~171×10-6),Ba(17.9×10-6~519×10-6)和Zr( < 19.9×10-6~133.0×10-6)变化较大,多数样品的Sr含量小于100×10-6,Zr含量小于80×10-6(图 6图 7);(2) Nb/Ta(1.5~5.8) 和Zr/Hf(12.3~27.9) 比值较低,明显低于球粒陨石(Nb/Ta=19.9;Zr/Hf=36.0)(图 8);(3) 较高的Rb/Sr(1.64~45.19),多数样品的Rb/Sr比值大于6.0,但Sr/Y(0.7~17.4) 比较低(图 6图 7)。

该期淡色花岗岩的87Sr/86Sr(t)较高且变化较大,为0.7537~1.0648,但εNd(t)值较低,和-13.8~-9.4(图 9)(Zhang et al., 2004a; King et al., 2011; Gao and Zeng, 2014; Gao et al., 2016)。虽然在Sr同位素组成上,与始新世花岗岩相比,该套淡色花岗岩明显不同,但两者具有相似的Nd同位素组成。值得指出的是,夏如渐新世淡色花岗岩具有异常高的初始Sr同位素组成(87Sr/86Sr(t)>0.8400)(Gao et al., 2016),可能反映了在喜马拉雅造山带折返早期,以黑云母脱水部分熔融为主的部分熔融作用。

在已经报道的渐新世淡色花岗岩中,除了个别样品保存较原始熔体特征外,多数样品都经历了不同程度的岩浆分异作用,表现为较低的Sr和Ba(图 6),较高的Rb和Rb/Sr比值(图 7),“海鸥状”的REE分配模式(图 5d)。

5.3 早中新世淡色花岗岩

沿喜马拉雅造山带,无论是高喜马拉雅带还是特提斯喜马拉雅带,该期淡色花岗岩都广泛发育,是喜马拉雅新生代花岗岩的主体(Le Fort, 1981; Schärer et al., 1986; Debon et al., 1986; Harris and Massey, 1994; Harris et al., 1995; Harrison et al., 1997, 1999a, b; Searle et al., 1999, 2010; Simpson et al., 2000; Daniel et al., 2003; Zhang et al., 2004a; Aoya et al., 2005; King et al., 2011; Aikman et al., 2012a, b; Guo and Wilson, 2012; Zeng et al., 2012; Gao and Zeng, 2014; Gao et al., 2017; 于俊杰等, 2011; 吴福元等, 2015; 胡古月等, 2016; 王晓先等, 2016; 及其中的参考文献)。几乎在每一个淡色花岗岩岩体,都发育含黑云母淡色花岗岩,淡色花岗岩和含石榴子石或电气石的淡色花岗岩,个别岩体还发育含绿柱石粗粒淡色花岗岩和伟晶岩。如位于错那北部的拿日雍措穹隆,从西到东呈现为一向西倾斜的淡色花岗岩岩体,从西部的二云母花岗岩,过渡到淡色花岗岩,再到含石榴子石/电气石淡色花岗岩,最后到最东边的含绿柱石伟晶花岗岩,同时局部含堆晶作用形成的电气石团块(电气石岩),主体花岗岩形成于20.0±0.1Ma,侵入到十字石石墨片岩的淡色花岗岩脉形成于21.8±0.3Ma(锆石U-Pb年龄)(高利娥等, 2017)。其他淡色花岗岩体也表现出相似特征。为探讨在喜马拉雅造山带构造演化过程中,中下地壳岩石部分熔融的类型和作用方式,需要甄别较原始的淡色花岗质熔体。与基性岩相比,甄别较原始的花岗质熔体特征难度较大,至今还未有统一的标准。在本文中,较原始的花岗质岩浆是指(1) 对于某类原岩,全岩花岗岩的矿物组成和主量元素地球化学特征较接近实验岩石学产生的熔体组分的那些花岗岩;或(2) 对某一岩体,全岩地球化学(主量元素和微量元素)组成演化程度最低的样品。根据各类淡色花岗岩的矿物组合和结构特征,结合各自的元素和同位素地球化学特征,依据表 1的判别原则,把这些淡色花岗岩分为四类:A类、B类、C类和F类淡色花岗岩。其中A和B类都代表较原始的岩浆组分,C类代表经历明显堆晶作用的淡色花岗岩,在矿物组成上,斜长石或钾长石含量明显较高,同时在手标本尺度上,常常含富集长石的条带或团块,类似现象在淡色花岗岩体中普遍存在。F类为经历明显分离结晶作用的淡色花岗岩,表现为缺乏云母矿物,斜长石的钠长石组分明显升高甚至以钠长石为主。各类花岗岩的元素和同位素地球化学特征(表 2)分述如下。值得指出的是,即使是相同源区和相同的部分熔融反应,由于部分熔融程度的差异,形成的淡色花岗质熔体都可以在主量和微量元素地球化学特征上表现出一定程度的内在变化,但总体上,不同类型的部分熔融反应所形成的熔体在主量和微量元素地球化学特征应表现出明显的差别。

5.3.1 A类淡色花岗岩:水致白云母部分熔融

在A类淡色花岗岩中,SiO2(69.5%~77.1%)、Al2O3(12.9%~16.2%)和CaO(0.6%~2.0%)都变化较大,但多数花岗岩的CaO大于1.0%;FeO( < 2.0%)、MgO( < 0.7%)、MnO( < 0.05%)和TiO2( < 0.31%)都较低(图 3)。除个别样品外,Na2O/K2O比值较高外,多数比值小于1.0,且A/CNK比值大于1.0,高达1.3(图 4),为过铝质花岗岩。

这些淡色花岗岩富集轻稀土(LREE),亏损重稀土(HREE)((La/Yb)N>12.5);显示较平直的重稀土分布特征((Gd/Yb)N=2.9~5.0) 和几乎无或弱负Eu异常(Eu/Eu*=0.7~0.9)(图 5e)。该类淡色花岗岩具有较高的Ba(>305×10-6)和Sr(>105×10-6),较低的Rb( < 270×10-6)和Rb/Sr比值( < 2.2)(图 6图 7)。同时,随着Ba浓度的增加,Rb/Sr比值几乎保持不变(图 6),与变质岩水致白云母部分熔融的特征一致(Inger and Harris, 1993; Zeng et al., 2005b, c; Gao and Zeng, 2014; Gao et al., 2017)。Nb/Ta和Zr/Hf比值都小于球粒陨石,分别为5.3~10.2和24.1~40.6(图 8)。较均一的放射性同位素(Sr和Nd)组成(图 9),其中87Sr/86Sr(t)=0.7391~0.7523,εNd(t)=-14.4~-13.7(Gao and Zeng, 2014; Gao et al., 2017; 王晓先等, 2015)。

5.3.2 B类淡色花岗岩:白云母脱水部分熔融

与A类淡色花岗岩相比,B类淡色花岗岩的SiO2(72.4%~74.2%)、Al2O3(14.2%~15.4%)都变化范围都较小,FeO、MgO、MnO和TiO2也都较低(图 3)。A/CNK比值同样大于1.0,高达1.28。除个别样品外,多数样品的Na2O/K2O比值小于1.0和CaO小于0.9%(图 3图 4)。

与A类淡色花岗岩不同的是,此类花岗岩的稀土元素地球化学特征较一致(图 5f),表现为(1) 轻重稀土分馏明显,富集轻稀土(LREE),亏损重稀土(HREE),(La/Yb)N=4.56~19.8;(2) 较显著的负Eu异常(Eu/Eu*=0.37~0.71)(图 5f)。B类淡色花岗岩具有较高的Rb(>325×10-6)(图 6),但较低的Ba( < 265×10-6)和Sr( < 76.0×10-6)(图 7)和Rb/Sr比值较高(>4.6)。Ba和Rb/Sr比值呈负相关关系,随Rb/Sr比值的升高,Ba浓度降低(图 7),反映了变泥质岩白云母脱水部分熔融的特征。Nb/Ta和Zr/Hf比值都低于球粒陨石且变化较大(图 8),分别为3.2~10.3和17.0~33.5。较均一的放射性同位素(Sr和Nd)组成(图 9),其中87Sr/86Sr(t)=0.7441~0.7604,εNd(t)=-14.7~-12.3(Gao and Zeng, 2014; Gao et al., 2017)。

5.3.3 C类淡色花岗岩:堆晶淡色花岗岩

除了与典型淡色花岗岩相似的低FeO、MgO和TiO2,高A/CNK等特征外(图 3图 4),该类花岗岩最显著地特征是显著的正Eu异常(图 5g)和多数样品具有高CaO(>1.1%)(图 3)或较高的K2O(>5.0%),代表斜长石堆晶或部分熔融过程中捕获源区转熔型或残留钾长石的效应。无论在主量元素,还是REE或微量元素地球化学特征上,该类花岗岩都呈现高度变化特征,如(1) SiO2(70.8%~76.3%)、Al2O3(13.6%~16.4%)和CaO(0.5%~2.7%)(图 3);(2) 较复杂的REE分配模式,从轻重稀土高度分馏到不分馏((La/Yb)N=1.0~21.1)(图 5g);(3) 高度变化的Ba(14×10-6~2467×10-6)、Sr(22×10-6~439×10-6),Rb(43×10-6~375×10-6)和Rb/Sr比值(0.2~15.9)(图 6图 7);(4) 高度变化的Nb/Ta(1.4~11.9) 和Zr/Hf(13.5~35.6) 比值(图 8)。

5.3.4 F类淡色花岗岩:高度分异淡色花岗岩

与前三类淡色花岗岩相比,该类花岗岩含石榴子石、电气石或绿柱石(图 2),与渐新世淡色花岗岩相似。该类花岗岩除了较高的SiO2外,其它主量元素和其它淡色花岗岩相似(图 3)。虽然在Rb/Sr-Sr和Rb/Sr-Ba的系统关系上与白云母脱水部分熔融形成的熔体相似(图 7),这些花岗岩最显著地特征是具有“海鸥型”稀土分布模式和强烈的负Eu异常(Eu/Eu* < 0.5)(图 5h)。除了个别样品外,多数样品具有较低的Sr( < 80×10-6)和Ba( < 153×10-6),较高Rb(>150×10-6)和Rb/Sr比值(>3.0)(图 6图 7),但较低的Nb/Ta( < 5.0) 和Zr/Hf( < 25.0) 比值(图 8)。

6 讨论 6.1 缩短增厚下地壳部分熔融

在喜马拉雅造山带的构造演化过程中,变泥质岩和角闪岩是两种易发生部分熔融作用的地壳物质。变泥质岩具有高Rb,低Sr和Rb/Sr比值的特征,角闪岩正好相反。因此,变泥质岩和角闪岩部分熔融产生的熔体具有明显不同的Rb-Sr系统关系(Zeng et al., 2011)。高喜马拉雅和NHGD其他穹窿地区的淡色花岗岩一般有:(1) 高Rb、低Sr以及较高的Rb/Sr比(>1.0);(2) 高初始Sr(87Sr/86Sr(t)>0.7300) 和低初始Nd(εNd(t) < -12) 同位素组成(Vidal et al., 1982; Debon et al., 1986; Deniel et al., 1987; Ferrara et al., 1991; Searle et al., 1997; 杨晓松和金振民, 2001; Zhang et al., 2004a; King et al., 2011);(3) 形成于小于30Ma的年龄(Noble and Searle, 1995; Harrison et al., 1997; Zhang et al., 2004a, b; King et al., 2011)等特征。雅拉香波和然巴始新世二云母花岗岩具有较低的Sr同位素组成(87Sr/86Sr(t) < 0.7120)(图 9),较老的形成年龄(>42.0Ma),要求与年轻( < 30Ma)淡色花岗岩不同的源区和部分熔融反应。

与大多数喜马拉雅淡色花岗岩不同,始新世二云母花岗岩具有高Sr/Y和La/Yb比值的特征,Sr和LREE含量较高,Y( < 10×10-6)和Yb( < 1.0×10-6)都较低的典型特征(Zeng et al., 2011; Liu et al., 2014; 高利娥等, 2009; 谢克家等, 2010)。结合喜马拉雅造山带碰撞早期经历了明显的缩短变形等构造过程(Aikman et al., 2008; Zhang et al., 2012),元素和Sr-Nd同位素地球化学特征表明:要形成类似喜马拉雅二云母花岗岩的岩浆,需要以角闪岩脱水部分熔融为主的部分熔融过程,同时石榴子石或角闪石呈残留相,但斜长石为熔融相。形成该类花岗岩的最有利条件是增厚下地壳构造背景。在雅拉香波穹隆(图 1b),混合岩化石榴角闪岩的部分熔融时代为43.5±1.3Ma(锆石U-Pb年龄)(Zeng et al., 2011),与雅拉香波的二云母花岗岩的形成时间相似,同时混合岩化石榴角闪岩具有与二云母花岗岩类似的Sr和Nd同位素组成,也支持变基性岩在增厚地壳条件下部分熔融形成始新世高Sr/Y比二云母花岗岩的推论。

6.2 A和B类淡色花岗岩的成因机制

实验岩石学(Douce and Harris, 1998; Knesel and Davison, 2002)、理论计算(Zeng et al., 2005a, b)和野外观测(Zhang et al., 2004a; King et al., 2011; Zeng et al., 2012)结果都表明:随着温压条件和含水量的变化,变泥质岩在无水或含水条件下发生从白云母到黑云母的递进部分熔融作用,形成性质各异、具有复杂Rb-Sr-Ba关系的花岗质岩浆(Harris and Inger, 1992; Knesel and Davison, 2002; Zeng et al., 2005b)。变泥质岩水致白云母熔融反应(MS2) 为7Qtz+15Pl+9Mus+xH2O→31M,白云母脱水部分熔融反应(MS1) 为:22Ms+7Pl+8Qtz=25Melt+5Kfs+5Sil+2Bt(Inger and Harris, 1993; Douce and Harris, 1998; Knesel and Davison, 2002)。与白云母脱水熔融形成的熔体相比,水致白云母熔融产生的熔体具有明显不同的Rb-Sr-Ba地球化学特征。在变泥质岩中,长石和云母分别是Sr和Rb的主要寄主矿物,与长石和全岩相比,云母具有较高的Rb/Sr比值和87Sr/86Sr同位素组成。Ca为斜长石的必要组成元素,Sr、Ba和Eu可以通过类质同象的方式替换Ca。因此,在变泥质岩水致白云母熔融过程中,与白云母脱水熔融相比,由于较多斜长石的加入,熔体具有较高的Sr、Ba和Ca,较低的Rb和Sr同位素比值(Inger and Harris, 1993; Knesel and Davidson, 2002; Zeng et al., 2005b)。白云母脱水部分熔融形成的熔体具有相反的特征。因此,花岗质岩浆的Rb-Sr-Ba关系是鉴定地壳岩石发生部分熔融作用类型的探针,可从地壳熔体的Rb-Sr-Ba系统关系来反推部分熔融类型和源岩性质(Inger and Harris, 1993; Knesel and Davidson, 2002; Zeng et al., 2005b, 2011)。

从Rb/Sr-Ba关系来看(图 7),A类花岗岩具有较高的Sr,较低的Rb和Rb/Sr比值。在Rb/Sr-Ba系统关系中,随着Ba浓度的增加,A类花岗岩Rb/Sr比值保持不变(图 7),与水致白云母部分熔融作用的特征一致(Inger and Harris, 1993),是变泥质岩水致白云母熔融的产物。如果年龄小于25Ma的喜马拉雅淡色花岗岩都源自相似的变泥质岩,B类淡色花岗岩具有较高的Rb(>325×10-6),但较低的Ba( < 265×10-6)和Sr( < 76.0×10-6)(图 7)和Rb/Sr比值较高(>4.6)。Ba和Rb/Sr比值呈负相关关系,随Rb/Sr比值的升高,Ba浓度降低(图 7),反应了变泥质岩白云母脱水部分熔融的特征。

上述两种部分熔融反应形成的熔体在主量、微量和同位素地球化学特征上都表现出明显的差异性。类似的现象在其它淡色花岗岩也有报道,如高喜马拉雅的玛纳斯鲁岩体。玛那斯鲁(Manaslu)淡色花岗岩体是喜马拉雅造山带中典型岩体之一,该岩体的淡色花岗岩存在两期,早期形成于22~24Ma,具有较高Rb(>300×10-6),较低的Sr(87×10-6~41×10-6)和较高的Rb/Sr(3.6~10.5) 和初始87Sr/86Sr(0.7523~0.7738) 比值;相反地,晚期淡色花岗岩形成于19.0~19.6Ma,具有较低的Rb( < 320×10-6),较高的Sr(83×10-6~114×10-6)和较低Rb/Sr(1.0~3.3) 和初始87Sr/86Sr(0.7446~0.7490) 比值(Vidal et al., 1982; Deniel et al., 1987; Guillot and Le Fort, 1995; Harrison et al., 1995, 1999a)。实验岩石学结果(Knesel and Davidson, 2002)和地球化学分析(Guillot and Le Fort, 1995)都表明:早期部分熔融为白云母脱水部分熔融,后期为含水白云母部分熔融的结果。

对比磷灰石和独居石在花岗质熔体的溶解模型(Harrison and Watson, 1984; Montel, 1993),影响这两种主要含磷和LREE载体矿物的溶解度的因素包括熔体温度,CaO活度和含水量。高温和干条件有利于磷灰石的溶解作用;湿条件有利于独居石的溶解作用(Zeng et al., 2005a, b)。与B类淡色花岗岩相比,A类淡色花岗岩具有较高的CaO,但较低的P2O5(图 4)和较低的Sm/Nd比值(图 9),反映了在含水部分熔融过程中,斜长石参与部分熔融比例较高,熔体的CaO含量和活度较高,抑制磷灰石溶解作用,但有利于独居石溶解,导致熔体的Sm/Nd和εNd值都较低(图 9)。

在获取了大量的地质年代学、全岩元素和同位素和矿物化学数据,甄别出两类地壳深熔作用—变泥质岩的白云母脱水熔融作用和水致白云母部分熔融作用。在大陆地壳的温度压力条件下,这两类地壳深熔作用是最常见、最普遍的深部地质过程,产生的熔体不仅在主量元素和微量元素,而且在同位素(Sr、Nd和Hf)和高场强元素(HFSE)上都表现出明显的差异性(Gao et al., 2017)。在这两类变泥质岩部分熔融作用中,主要造岩矿物和副矿物的溶解行为的差异性是调控熔体地球化学特征差异性的决定性因素。

7 未来研究方向和关键课题 7.1 花岗质片麻岩在地壳深熔作用中的作用

喜马拉雅造山带发育大量花岗质片麻岩,原岩主要形成于寒武世-早奥陶世(Miller et al., 2001; Gehrels et al., 2006a, b; Cawood et al., 2007; Wang et al., 2012; 许志琴等, 2005; 张泽明等, 2008; 王誉桦等, 2014; 高利娥等, 2015)和新元古代(Singh et al., 2002; Ahmad et al., 2013; Wang et al., 2017)。许多花岗片麻岩经历了广泛的混合岩化作用,发育含石榴子石的淡色体(Pognante et al., 1990),同时还记录了20~18Ma的变质和部分熔融年龄(许志琴等, 2005),表明花岗质片麻岩也可能发生深熔作用,形成花岗质岩浆。但是,花岗质片麻岩在淡色花岗岩形成中的作用并不清楚。在日玛那穹窿和麻珈穹窿中,中新世淡色花岗岩具有特殊的Sr-Nd同位素系统特征(Zhang et al., 2004a; 于俊杰等, 2011; 胡古月等, 2016)。于俊杰等(2011)利用Sr-Nd同位素组成二元混合计算得出,花岗片麻岩参与部分熔融可以形成类似于淡色花岗岩成分的岩浆,定结淡色花岗岩的源岩可能是有花岗片麻岩和变泥质岩混合形成。虽然目前大多数学者认为高喜马拉雅变沉积岩是淡色花岗岩的主要源岩,但有限数据已揭示在特定阶段或构造背景下,花岗质片麻岩可能发生部分熔融,形成具有异常Sr和Nd同位素组成的淡色花岗岩(Zhang et al., 2004a; Gao et al., 2016; 于俊杰等, 2011)。但与变泥质岩相比,多数花岗片麻岩含云母和斜长石组分较少,较难熔融,只有在较高温或富流体的条件下,发生较高程度的部分熔融,形成淡色花岗质岩浆。

因此,有必要加强花岗质片麻岩的元素和同位素地球化学特征的研究,确定它们的变质和部分熔融反应的结构和时限,限定它们的Sr和Nd同位素组成,深入了解它们在喜马拉雅新生代地壳深熔作用中扮演的角色,甄别来自花岗质片麻岩的贡献程度,有可能在识别喜马拉雅造山带和其它碰撞造山带可能经历过的特殊构造过程提供重要的限定。

7.2 副矿物在地壳深熔作用中的作用

和地幔部分熔融不同,地壳深熔作用往往通过一系列复杂的部分熔融反应来完成。富水原岩的部分熔融反应一般在较低温度下发生,熔融反应主要以石英+长石±白云母+水(水致部分熔融)进行。缺水或干燥原岩的熔融反应主要以含水矿物(如白云母、黑云母及角闪石)的脱水熔融反应进行。在变泥质岩、变杂砂岩和变基性岩中,上述两种反应是产生花岗岩质岩浆的重要机制。由于常见的副矿物(锆石、磷灰石、独居石等)是REE和关键放射性同位素(Nd、Hf)的主要赋存矿物,在这些反应中,副矿物是否参与部分熔融强烈影响着花岗岩质岩浆的关键微量元素和同位素(Sr、Nd、Hf和Pb)地球化学特征具有控制作用(Watson and Harrison, 1984; Ayres and Harris, 1997; Zeng et al., 2005a, b, 2012; Perini et al., 2009; Gao et al., 2017)。例如,在地壳深熔作用中,锆石影响花岗质熔体的Pb及Hf,石榴子石影响Hf,磷灰石及独居石影响Nd的同位素地球化学特征。通常条件下,由于不同副矿物的差异溶解行为,地壳物质部分熔融形成的熔体较源岩的LREE和HREE都低,但和源岩相比,熔体的LREE斜率((La/Gd)N)比源岩的陡, 导致残留体的(La/Gd)N小于源岩,但LREE含量增高。

野外、元素地球化学及放射性同位素地球化学观测揭示了许多具有独特地球化学特征的花岗岩岩体,如喜马拉雅的玛纳斯鲁(Malasulu)淡色花岗岩体、美国黑山(Black Hills)的Harney Peak花岗岩体(Shearer and Papike, 1987; Nabelek and Glasock, 1995; Harrison et al., 1999a, b)。相对于典型的大陆上地壳,这些花岗岩具有:(1) 强烈亏损轻稀土;(2) 高度的Sm-Nd的分馏;(3) 与邻近同时代花岗岩差别明显的初始143Nd/144Nd比值。这种怪异的地球化学特征往往归结于:岩浆后期副矿物分离结晶作用或在较深部部分熔融区,副矿物的差异溶解作用。甄别这些不同的岩浆作用,需要综合考虑副矿物的溶解动力学模型。

除了副矿物本身在深熔熔体中的溶解动力学行为外,副矿物在源岩中的赋存方式也是影响深熔熔体地球化学特征的因素之一。在变质岩中,主要造岩矿物(白云母、黑云母)通常都包含大量的副矿物(Watson et al., 1989; Watt and Harley, 1993),在部分熔融过程中,这些造岩矿物是否分解和参与部分熔融,释放包裹的副矿物并参与部分熔融是导致熔体较复杂的元素和同位素地球化学特征的原因。例如,在高级变质岩中,大部分独居石包裹在黑云母或石榴子石中。由于白云母脱水熔融反应温度普遍地低于黑云母脱水熔融反应,在白云母脱水部分熔融过程中,由于黑云母和石榴子石是非熔融相,多数独居石并不释放到熔体中,不参与部分熔融反应,导致岩浆与独居石的元素和同位素不平衡,产生的岩浆具有强烈亏损LREE的特征。相反地,在较高温情况下,黑云母的分解及熔融反应施放包含的独居石,独居石参与熔融反应,岩浆具有较高的LREE和Th(Nabelek and Glasock, 1995)。

在过去几十年来,已有一系列的实验研究了磷灰石、独居石、锆石及含Nb-Ta矿物在花岗质熔体的溶解行为(Watson and Harrison, 1983, 1984; Harrison and Watson, 1984; Rapp and Watson, 1986; Montel, 1993; Wolf and London, 1994, 1995; Linnen and Keppler, 1997, 2002; Linnen, 1998, 2005; Chevychelov et al., 2010; Bartels et al., 2010; Fiege et al., 2011; Linnen et al., 2012; Boehnke et al., 2013),为建立符合地质情况的地壳深熔作用的理论模型提供了关键的限定。如地壳深熔作用中,Pb、Nd和Hf同位素地球化学行为(Hogan and Sinha, 1991; Watson, 1996; Zeng et al., 2005a; Tang et al., 2014),但是至今还缺乏石榴子石和榍石在花岗质熔体的溶解行为的实验研究。

石榴子石是众多岩石类型中最常见的造岩矿物之一,在较宽的温度和压力条件下,都稳定存在。石榴子石不仅记录了生长时的时间和物理条件,而且即使在后期热扰动下,还常常保留了这些信息。和其它矿物相比,石榴子石具有以下特性:(1) 较高的Sm/Nd比值和Nd同位素比值;(2) 常具有化学成分环带,从核部到边部,无论是主量元素,还是微量元素(Sc、Y、HREE),甚至Nd和Hf放射性同位素常表现出环带特征;(3) Sm、Nd、Lu和Hf在石榴子石的扩散速率较低(Burton et al., 1995);在部分熔融中,无论是以固态残留相还是以熔融相存在,除了影响熔体中与石榴子石相容的元素(HREE、Y、Sc)的地球化学特征外,还可能导致Lu-Hf和Sm-Nd同位素系统的分馏和熔体的Hf和Nd同位素不平衡(Zeng et al., 2012)。在喜马拉雅造山带中,许多变质岩都含石榴子石,如果石榴子石形成时间远早于部分熔融时间,石榴子石将具有较高的Sm/Nd和143Nd/144Nd比值。在部分熔融过程中,如果石榴子石不参与反应,会降低熔体的143Nd/144Nd比值。如果熔体抽离,远离残留型石榴子石,那么熔体的Nd同位素就会出现同位素不平衡。但是,如果这些石榴子石被带入熔体,对淡色体的微量元素和放射性同位素产生影响。如南迦巴瓦~25Ma的高钙淡色花岗岩,与变泥质岩相比,含石榴子石的淡色花岗岩明显富集石榴子石相容元素(Sc、Y和HREE)和较高的εNd值,是变泥质岩部分熔融过程中,花岗质熔体携带源区石榴子石的结果(Zeng et al., 2012)。

除了石榴子石外,榍石是基性-花岗闪长质变质岩中的常见副矿物,富含Sr、REE、Nb和Ta。在这些岩石的部分熔融过程中,榍石的溶解行为及其如何影响花岗质熔体的Sr和Nd同位素组成及Nb-Ta的系统关系也需要野外实测和实验岩石学的研究。

7.3 高分异淡色花岗岩元素地球化学行为和成矿效应

在花岗质岩浆从正常岩浆向高硅岩浆(通常SiO2>72.0%)演化过程中,岩浆不仅SiO2升高,同时挥发分(H2O、Cl或F)也随之升高,熔体结构发生实质性的变化,导致矿物组成和某些关键元素的地球化学行为发生变异,包括:(1) 主要造岩矿物相溶解度的改变(Ren et al., 2012);(2) 关键微量元素(Rb、Sr、Ba、Cs、REE等)分配系数的变化(Icenhower and London, 1995; Glazner et al., 2008; Ren et al., 2012);(3) 副矿物(如锆石、独居石、磷灰石、铪石等)溶解行为的变化(Ayers and Watson, 1991; Wolf and London, 1994, 1995; Bea, 1996; Linnen and Keppler, 1997, 2002);(4) 结晶石榴子石、电气石或绿柱石(Liu et al., 2016; 高利娥等, 2012)和(5) 金属元素(Nb、Ta、Sn等)的富集(Taylor and Wall, 1992)。

多数淡色花岗岩具有较低的Zr/Hf比值,导致高分异淡色花岗岩Zr/Hf比值降低的原因包括:(1) 锆石的分离结晶作用(Dostal and Chatterjee, 2000; Claiborne et al., 2006)或(2) 锆石(ZrSiO4)和铪石(HfSiO4)的溶解行为的变化(Linnen and Keppler, 2002)。Zr和Hf虽然具有相似的地球化学性质,但是随着分异作用的不断增强,花岗质岩浆的结构将逐渐发生变化,在达到熔体-流体相互作用阶段时,锆石(ZrSiO4)和铪石(HfSiO4)的溶解行为随之也会发生变化,虽然Zr和Hf在高分异花岗质岩浆中的溶解度都升高,熔体中的Zr和Hf的浓度都升高,但Hf的溶解度将明显高于Zr的溶解度(Linnen and Keppler, 2002),导致熔体的Zr/Hf比值逐渐变低。类似的情形也适用于Na和Ta孪生元素对(Linnen and Keppler, 1997; Bartels et al., 2010; Fiege et al., 2011)。但在确定高分异花岗岩的挥发分组成和熔体结构变化特征等方面的研究明显不足。

在上述高分异淡色花岗岩中,部分样品具有明显升高的B、Be、Sn,Nb和Ta含量(图 10),这些淡色花岗岩含电气石、绿柱石、锡石、钶铁矿和钽铁矿等,显示具有Be、B、Sn、Nb和Ta等元素富集和矿化特征,需要进一步的地质调查和研究,来确定喜马拉雅淡色花岗岩,尤其是高分异淡色花岗岩的关键元素,如B-Be、Nb-Ta和Zr-Hf等元素对的分异和富集行为及其控制因素,厘定它们的成矿潜力。

图 10 喜马拉雅新生代淡色花岗岩Sn (a)、W (b)、Nb (c)和Ta (d)和Rb/Sr比值关系图 Fig. 10 Sn (a), W (b), Nb (c) and Ta (d) versus Rb/Sr diagrams for the Himalayan Cenozoic leucogranites
8 结论

在众多碰撞造山带中,壳源花岗岩的岩石学和地球化学特征取决于温度-压力-成分(P-T-X),P-T条件与造山带物质所经历的构造相紧密相关。虽然在细节上可能存在差异,多数碰撞型造山带都经历了相似的顺时针P-T演化历史。早期为同碰撞缩短增厚作用,表现为高压-低温的P-T条件,在外来热的作用下,下地壳物质可在石榴子石或角闪石稳定域内发生部分熔融,形成高CaO,高Na/K和Sr/Y比的花岗岩;后期经历伸展垮塌作用,为低压-高温阶段,取决于含水条件,可能发生水致部分熔融和脱水部分熔融。与其他碰撞造山带相似,喜马拉雅造山带也经历了顺时针P-T-t构造演化和从早期缩短增厚向晚期伸展构造体制的转换,伴随相应的深熔作用,具体表现为(1) 在地壳缩短阶段(>35Ma),发生高角闪岩相-榴辉岩相的变质作用(55~46Ma, De Sigoyer et al., 2000; Kaneko et al., 2003; Parrish et al., 2006; Gao et al., 2012),增厚基性下地壳发生部分熔融作用, 形成高Na/K和Sr/Y比的过铝质花岗岩(46~35 Ma, Ding et al., 2005; Aikman et al., 2008; Zeng et al., 2011, 2015; Hou et al., 2012; Liu et al., 2014; 戚学祥等, 2008; 高利娥等, 2009, 2010);(2) 壳内伸展变形阶段,高喜马拉雅结晶岩系在快速折返过程中发生变泥质岩减压脱水部分熔融作用,形成高Rb/Sr、低Sr/Y比的过铝质花岗岩(28~10Ma; Le Fort, 1981; Deniel et al., 1987; Harris and Massey, 1994; Harris et al., 1995; Harrison et al., 1997, 1999a; Douce and Harris, 1998; Zhang et al., 2004b; Aoya et al., 2005; King et al., 2011; Aikman et al., 2012a, b);(3) 在部分地区,变泥质岩发生水致部分熔融作用(Harrison et al., 1999b; Prince et al., 2001; Zhang et al., 2004a; King et al., 2011; Zeng et al., 2012; Gao and Zeng, 2014; Gao et al., 2017),可能与藏南裂谷系的开启有关。始新世-早渐新世为碰撞缩短作用,中-下地壳物质发生部分熔融作用,这些熔体有效地改变了深部岩石的构造物理性质,促使地壳变形从缩短增厚向伸展减薄转换,诱发STDS的启动(Lee and Whitehouse, 2007; Yang et al., 2009; Zeng et al., 2011; Zhang et al., 2012)。伴随着STDS进一步伸展,深部地壳物质快速剥露,导致变泥质岩发生大规模减压部分熔融作用,形成了大量中新世典型的喜马拉雅淡色花岗岩(Schärer et al., 1986; Harrison et al., 1999a; Zhang et al., 2004a, 2012; Aoya et al., 2005; Lee and Whitehouse, 2007; King et al., 2011; Wang et al., 2012; Gao et al., 2013; Gao and Zeng, 2014; 于俊杰等, 2011)。由此可见,淡色花岗岩的形成与构造转换密切相关,地壳物质的深熔作用引起了构造变形的转换,同时构造环境控制了淡色花岗岩的地球化学特征。因此,准确厘定淡色花岗岩的地球化学特征和形成时代以及产出的构造背景是探讨喜马拉雅造山带及其他碰撞造山带构造演化过程的重要途径之一(Lee and Whitehouse, 2007; Zeng et al., 2011; Hou et al., 2012; Gao et al., 2017)。

致谢 感谢翟明国院士热情约稿和张旗研究员提出众多建设性修改意见。
参考文献
[] Ahmad I, Khan S, Lapen T, Burke K, Jehan N. 2013. Isotopic ages for alkaline igneous rocks, including a 26Ma ignimbrite, from the Peshawar plain of northern Pakistan and their tectonic implications. Journal of Asian Earth Sciences, 62: 414–424. DOI:10.1016/j.jseaes.2012.10.025
[] Ahmad T, Harris N, Bickle M, Chapman H, Bunbury J, Prince C. 2000. Isotopic constraints on the structural relationships between the Lesser Himalayan Series and the High Himalayan Crystalline Series, Garhwal Himalaya. Geological Society of America Bulletin, 112(3): 467–477. DOI:10.1130/0016-7606(2000)112<467:ICOTSR>2.0.CO;2
[] Aikman AB, Harrison TM, Ding L. 2008. Evidence for early (>44Ma) Himalayan crustal thickening, Tethyan Himalaya, southeastern Tibet. Earth and Planetary Science Letters, 274(1-2): 14–23. DOI:10.1016/j.epsl.2008.06.038
[] Aikman AB, Harrison TM, Hermann J. 2012a. The origin of Eo-and Neo-Himalayan granitoids, eastern Tibet. Journal of Asian Earth Sciences, 58: 143–157. DOI:10.1016/j.jseaes.2012.05.018
[] Aikman AB, Harrison TM, Hermann J. 2012b. Age and thermal history of Eo-and Neohimalayan granitoids, eastern Himalaya. Journal of Asian Earth Sciences, 51: 85–97. DOI:10.1016/j.jseaes.2012.01.011
[] Aoya M, Wallis SR, Terada K, Lee J, Kawakami T, Wang Y, Heizler M. 2005. North-south extension in the Tibetan crust triggered by granite emplacement. Geology, 33(11): 853–856. DOI:10.1130/G21806.1
[] Ayers JC, Watson EB. 1991. Solubility of apatite, monazite, zircon, and rutile in supercritical aqueous fluids with implications for subduction zone geochemistry. Philosophical Transactions: Physical Sciences and Engineering, 355(1638): 365–375.
[] Ayres M, Harris N. 1997. REE fractionation and Nd-isotope disequilibrium during crustal anatexis: Constraints from Himalayan leucogranites. Chemical Geology, 139(1-4): 249–269. DOI:10.1016/S0009-2541(97)00038-7
[] Ballouard C, Poujol M, Boulvais P, Branquet Y, Tartèse R, Vigneresse JL. 2016. Nb-Ta fractionation in peraluminous granites: A marker of the magmatic-hydrothermal transition. Geology, 44(3): 231–234. DOI:10.1130/G37475.1
[] Bartels A, Holtz F, Linnen RL. 2010. Solubility of manganotantalite and manganocolumbite in pegmatitic melts. American Mineralogist, 95(4): 537–544. DOI:10.2138/am.2010.3157
[] Bea F. 1996. Controls on the trace element composition of crustal melts. Transactions of the Royal Society of Edinburg. Earth Sciences, 87(1-2): 33–41. DOI:10.1017/S0263593300006453
[] Blisniuk MP, Hacker RB, Glodny J, Ratschbacher L, Bi SW, Wu ZH, McWilliams MO, Calvert A. 2001. Normal faulting in central Tibet since at least 13. 5Myr ago. Nature, 412(6847): 628–632. DOI:10.1038/35088045
[] Boehnke P, Watson EB, Trail D, Harrison TM, Schmitt AK. 2013. Zircon saturation re-revisited. Chemical Geology, 351: 324–334. DOI:10.1016/j.chemgeo.2013.05.028
[] Brown M. 2007. Crustal melting and melt extraction, ascent and emplacement in orogens: Mechanisms and consequences. Journal of the Geological Society, 164(4): 709–730. DOI:10.1144/0016-76492006-171
[] Brown M. 2010. Melting of the continental crust during orogenesis: The thermal, rheological, and compositional consequences of melt transport from lower to upper continental crust. Canadian Journal of Earth Sciences, 47(5): 655–694. DOI:10.1139/E09-057
[] Brown M. 2013. Granite: From genesis to emplacement. Geological Society of America Bulletin, 125(7-8): 1079–1113. DOI:10.1130/B30877.1
[] Burton KW, Kohn MJ, Cohen AS, O'Nions K. 1995. The relative diffusion of Pb, Nd, Sr and O in garnet. Earth and Planetary Science Letters, 133(1-2): 199–211. DOI:10.1016/0012-821X(95)00067-M
[] Carosi R, Montomoli C, Langone A, Turina A, Cesare B, Iaccarino S, Fascioli L, Visonà D, Ronchi AR, Rai SM. 2014. Eocene partial melting recorded in peritectic garnets from kyanite-gneiss, greater Himalayan sequence, central Nepal. In: Gómez-Tuena A, Straub SM and Zellmer GF (eds.). Orogenic Andesites and Crustal Growth. Geological Society, London, Special Publications, 412(1): 111–129.
[] Cawood PA, Johnson MRW, Nemchin AA. 2007. Early Palaeozoic orogenesis along the Indian margin of Gondwana: Tectonic response to Gondwana assembly. Earth and Planetary Science Letters, 255(1-2): 70–84. DOI:10.1016/j.epsl.2006.12.006
[] Chevychelov VY, Borodulin GP, Zaraisky GP. 2010. Solubility of columbite, (Mn, Fe)(Nb, Ta)2O6, in granitoid and alkaline melts at 650~850℃ and 30~400MPa: An experimental investigation. Geochemistry International, 48(5): 456–464. DOI:10.1134/S0016702910050034
[] Chu MF, Chung SL, O'Reilly SY, Pearson NJ, Wu FY, Li XH, Liu DY, Ji JQ, Chu CH, Lee HY. 2011. India's hidden inputs to Tibetan orogeny revealed by Hf isotopes of Transhimalayan zircons and host rocks. Earth and Planetary Science Letters, 307(3-4): 479–486. DOI:10.1016/j.epsl.2011.05.020
[] Claiborne LL, Miller CF, Walker BA, Wooden JL, Mazdab FK, Bea F. 2006. Tracking magmatic processes through Zr/Hf ratios in rocks and Hf and Ti zoning in zircons: An example from the Spirit Mountain batholith, Nevada. Mineralogical Magazine, 70(5): 517–543. DOI:10.1180/0026461067050348
[] Clarke DB. 2007. Assimilation of xenocrysts in granitic magmas: Principles, processes, proxies, and problems. Canadian Mineralogist, 45(1): 5–30. DOI:10.2113/gscanmin.45.1.5
[] Coleman M, Hodges K. 1995. Evidence for Tibetan Plateau uplift before 14Myr ago from a new minimum age for east-west extension. Nature, 374(6517): 49–52. DOI:10.1038/374049a0
[] Conrad WK, Nicholls IA, Wall VJ. 1988. Water-saturated and-undersaturated melting of metaluminous and peraluminous crustal compositions at 10kb: Evidence for the origin of silicic magmas in the Taupo volcanic zone, New Zealand, and other occurrences. Journal of Petrology, 29(4): 765–803. DOI:10.1093/petrology/29.4.765
[] Daniel CG, Hollister LS, Parrish RR, Grujic D. 2003. Exhumation of the main central thrust from lower crustal depths, eastern Bhutan Himalaya. Journal of Metamorphic Geology, 21(4): 317–334. DOI:10.1046/j.1525-1314.2003.00445.x
[] De Sigoyer J, Chavagnac V, Blichert-Toft J, Villa IM, Luais B, Guillot S, Cosca M, Mascle G. 2000. Dating the Indian continental subduction and collisional thickening in the northwest Himalaya: Multichronology of the Tso Morari eclogites. Geology, 28(6): 487–490. DOI:10.1130/0091-7613(2000)28<487:DTICSA>2.0.CO;2
[] Debon F, Le Fort P, Sheppard SMF, Sonet J. 1986. The four plutonic belts of the Transhimalaya-Himalaya: A chemical, mineralogical, isotopic, and chronological synthesis along a Tibet-Nepal Section. Journal of Petrology, 27(1): 219–250. DOI:10.1093/petrology/27.1.219
[] Deniel C, Vidal P, Fernandez A, Le Fort P, Peucat JJ. 1987. Isotopic study of the Manaslu granite (Himalaya, Nepal): Inferences on the age and source of Himalayan leucogranites. Contributions to Mineralogy and Petrology, 96(1): 78–92. DOI:10.1007/BF00375529
[] Diedesch TF, Jessup MJ, Cottle JM, Zeng LS. 2016. Tectonic evolution of the middle crust in southern Tibet from structural and kinematic studies in the Lhagoi Kangri gneiss dome. Lithosphere, 8(5): 480–504. DOI:10.1130/L506.1
[] Ding L, Kapp P, Wan XQ. 2005. Paleocene-Eocene record of ophiolite obduction and initial India-Asia collision, south central Tibet. Tectonics, 24(3): 1–18.
[] DiPietro JA, Isachsen CE. 2001. U-Pb zircon ages from the Indian plate in Northwest Pakistan and their significance to Himalayan and pre-Himalayan geologic history. Tectonics, 20(4): 510–525. DOI:10.1029/2000TC001193
[] Dostal J, Chatterjee AK. 2000. Contrasting behaviour of Nb/Ta and Zr/Hf ratios in a peraluminous granitic pluton (Nova Scotia, Canada). Chemical Geology, 163(1-4): 207–218. DOI:10.1016/S0009-2541(99)00113-8
[] Douce AEP, 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
[] Douce AEP, Beard JS. 1995. Dehydration-melting of biotite gneiss and quartz amphibolite from 3 to 15 kbar. Journal of Petrology, 36(3): 707–738. DOI:10.1093/petrology/36.3.707
[] Douce AEP. 1996. Effects of pressure and H2O content on the compositions of primary crustal melts. Transactions of the Royal Society of Edinburgh: Earth Sciences, 87(1-2): 11–21. DOI:10.1017/S026359330000643X
[] Douce AEP, Beard JS. 1996. Effects of P, f(O2) and Mg/Fe ratio on dehydration melting of model metagreywackes. Journal of Petrology, 37(5): 999–1024. DOI:10.1093/petrology/37.5.999
[] Douce AEP, Harris N. 1998. Experimental constraints on Himalayan Anatexis. Journal of Petrology, 39(4): 689–710. DOI:10.1093/petroj/39.4.689
[] Ducea MN, Saleeby JB. 1996. Buoyancy sources for a large unrooted mountain range, the Sierra Nevada, California: Evidence from xenolith thermobarometry. Journal of Geophysical Research: Solid Earth, 101: 8229–8241. DOI:10.1029/95JB03452
[] Ferrara G, Lombardo B, Tonarini S, Turi B. 1991. Sr, Nd and O isotopic characterization of the Gophu La and Gumburanjon leucogranites (High Himalaya). Schweizer Mineralogische and Petrographische Mitteilungen, 71: 35–51.
[] Fiege A, Kirchner C, Holtz F, Linnen RL, Dziony W. 2011. Influence of fluorine on the solubility of manganotantalite (MnTa2O6) and manganocolumbite (MnNb2O6) in granitic melts: An experimental study. Lithos, 122(3-4): 165–174. DOI:10.1016/j.lithos.2010.12.012
[] Gao LE, Zeng LS, Liu J, 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, Hu GY. 2010. High Sr/Y two-mica granite from Quedang area, southern Tibet: Formation mechanism and tectonic implication. Geological Bulletin of China, 29(2-3): 214–226.
[] Gao LE, Zeng LS, Xie KJ. 2012. Eocene high grade metamorphism and crustal anatexis in the North Himalaya Gneiss Domes, Southern Tibet. Chinese Science Bulletin, 57(6): 639–650. DOI:10.1007/s11434-011-4805-4
[] Gao LE, Zeng LS, Shi WG, Chen ZY, Hu MY, 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 LE, Zeng LS, Hou KJ, Guo CL, Tang SH, Xie KJ, Hu GY, Wang L. 2013. Episodic crustal anatexis and the formation of Paiku composite leucogranitic pluton in the Malashan Gneiss Dome, Southern Tibet. Chinese Science Bulletin, 58(28-29): 3546–3563. DOI:10.1007/s11434-013-5792-4
[] Gao LE, Zeng LS. 2014. Fluxed melting of metapelite and the formation of Miocene high-CaO two-mica granites in the Malashan gneiss dome, southern Tibet. Geochimica et Cosmochimica Acta, 130: 136–155. DOI:10.1016/j.gca.2014.01.003
[] Gao LE, Zeng LS, Xu ZQ, Wang L. 2015. Himalaya in the Caledonia Time: A record from the Malashan-Gyirong area, southern Tibet. Acta Petrologica Sinica, 31(5): 1200–1218.
[] Gao LE, Zeng LS, Gao JH, Shang Z, Hou KJ, Wang Q. 2016. Oligocene crustal anatexis in the Tethyan Himalaya, southern Tibet. Lithos, 264: 201–209. DOI:10.1016/j.lithos.2016.08.038
[] Gao LE, Zeng LS, Wang L, Hou KJ, Gao JH, Shang Z. 2016. Timing of different crustal partial melting in the Himalayan orogenic belt and its tectonic implications. Acta Geologica Sinica, 90(11): 3039–3059.
[] Gao LE, Zeng LS, Asimow PD. 2017. Contrasting geochemical signatures of fluid-absent versus fluid-fluxed melting of muscovite in metasedimentary sources: The Himalayan leucogranites. Geology, 45(1): 39–42. DOI:10.1130/G38336.1
[] Gao LE, Zeng LS, Gao JH, Zhao LH, Hou KJ, Tang SH. 2017. The Miocene leucogranite in the Nariyongcuo Gneiss Domes, southern Tibet: Products from melting metapelite and fractional crystallization. Acta Petrologica Sinica, 33(7).
[] Gardien V, Thompson AB, Grujic D, Ulmer P. 1995. Experimental melting of biotite+plagioclase+quartz±muscovite assemblages and implications for crustal melting. Journal of Geophysical Research-Solid Earth, 100(B8): 15581–15591. DOI:10.1029/95JB00916
[] Gardien V, Thompson AB, Ulmer P. 2000. Melting of biotite+plagioclase+quartz gneisses: The role of H2O in the stability of amphibole. Journal of Petrology, 41(5): 651–666. DOI:10.1093/petrology/41.5.651
[] Gehrels GE, Decelles PG, Ojha TP, Upreti BN. 2006a. Geologic and U-Th-Pb geochronologic evidence for early Paleozoic tectonism in the Kathmandu thrust sheet, central Nepal Himalaya. Geological Society of America Bulletin, 118(1-2): 185–198. DOI:10.1130/B25753.1
[] Gehrels GE, DeCelles PG, Ojha TP, Upreti BN. 2006b. Geologic and U-Pb geochronologic evidence for Early Paleozoic tectonism in the Dadeldhura thrust sheet, far-west Nepal Himalaya. Journal of Asian Earth Sciences, 28(4-6): 385–408. DOI:10.1016/j.jseaes.2005.09.012
[] Glazner AF, Coleman DS, Bartley JM. 2008. The tenuous connection between high-silica rhyolites and granodiorite plutons. Geology, 36(2): 183–186. DOI:10.1130/G24496A.1
[] Green TH, Pearson NJ. 1987. An experimental study of Nb and Ta partitioning between Ti-rich minerals and silicate liquids at high pressure and temperature. Geochimica et Cosmochimica Acta, 51(1): 55–62. DOI:10.1016/0016-7037(87)90006-8
[] Guillot S, Le Fort P. 1995. Geochemical constraints on the bimodal origin of High Himalayan leucogranites. Lithos, 35(3-4): 221–234. DOI:10.1016/0024-4937(94)00052-4
[] Guo CL, Zeng LS, Gao LE, Su HZ, Ma XH and Yin B. 2017. Mineral and whole-rock geochemistry as mineralization prospecting indicators for highly fractionated granites: Taking comparison of Hetian and Xihuashan granites as an example. Acta Geologica Sinica, (8): in press (in Chinese with English abstract)
[] Guo ZF, Wilson M. 2012. The Himalayan leucogranites: Constraints on the nature of their crustal source region and geodynamic setting. Gondwana Research, 22(2): 360–376. DOI:10.1016/j.gr.2011.07.027
[] Harris N, Massey J. 1994. Decompression and anatexis of Himalayan metapelites. Tectonics, 13(6): 1537–1546. DOI:10.1029/94TC01611
[] Harris N, Ayres M, Massey J. 1995. Geochemistry of granitic melts produced during the incongruent melting of muscovite: Implications for the extraction of Himalayan leucogranite magmas. Journal of Geophysical Research: Solid Earth, 100(B8): 15767–15777. DOI:10.1029/94JB02623
[] Harris NBW, Inger S. 1992. Trace element modelling of pelite-derived granites. Contributions to Mineralogy and Petrology, 110(1): 46–56. DOI:10.1007/BF00310881
[] Harrison TM, Watson EB. 1984. The behavior of apatite during crustal anatexis: Equilibrium and kinetic considerations. Geochimica et Cosmochimica Acta, 48(7): 1467–1478. DOI:10.1016/0016-7037(84)90403-4
[] Harrison TM, Copeland P, Kidd WSF, Lovera OM. 1995. Activation of the Nyainqentanghla shear zone: Implications for uplift of the southern Tibetan Plateau. Tectonics, 14(3): 658–676. DOI:10.1029/95TC00608
[] Harrison TM, Lovera OM, Grove M. 1997. New insights into the origin of two contrasting Himalayan granite belts. Geology, 25(10): 899–902. DOI:10.1130/0091-7613(1997)025<0899:NIITOO>2.3.CO;2
[] Harrison TM, Grove M, Lovera OM, Catlos EJ, D'Andrea J. 1999a. The origin of Himalayan anatexis and inverted metamorphism: Models and constraints. Journal of Asian Earth Sciences, 17(5-6): 755–772. DOI:10.1016/S1367-9120(99)00018-8
[] Harrison TM, Grove M, Mckeegan KD, Coath CD, Lovera O, Le Fort P. 1999b. Origin and episodic emplacement of the Manaslu Intrusive Complex, Central Himalaya. Journal of Petrology, 40(1): 3–19. DOI:10.1093/petroj/40.1.3
[] Hodges KV. 2000. Tectonics of the Himalaya and southern Tibet from two perspectives. Geological Society of America Bulletin, 112(3): 324–350. DOI:10.1130/0016-7606(2000)112<324:TOTHAS>2.0.CO;2
[] Hogan JP, Sinha AK. 1991. The effect of accessory minerals on the redistribution of lead isotopes during crustal anatexis: A model. Geochimica et Cosmochimica Acta, 55(1): 335–348. DOI:10.1016/0016-7037(91)90422-2
[] Hou ZQ, Zheng YC, Zeng LS, Gao LE, Huang KX, Li W, Li QY, Fu Q, Liang W, Sun QZ. 2012. Eocene-Oligocene granitoids in southern Tibet: Constraints on crustal anatexis and tectonic evolution of the Himalayan orogen. Earth and Planetary Science Letters, 349-350: 38–52. DOI:10.1016/j.epsl.2012.06.030
[] Hu GY, Zeng LS, Qi XX, Hou KJ, Gao LE. 2011. The Mid-Eocene subvolcanic field in the Lhunze-Qiaga area, Tethyan Himalaya, southern Tibet: A high-level magmatic suite related to the Yardio two-mica granite. Acta Petrologica Sinica, 27(11): 3308–3318.
[] Hu GY, Zeng LS, Chen H, Gao LE. 2016. Assembling mechanism for the Late-Miocene Mabja leucogranite pluton, Himalayan orogen. Acta Geologica Sinica, 90(8): 1737–1754.
[] Icenhower J, London D. 1995. An experimental study of element partitioning among biotite, muscovite, and coexisting peraluminous silicic melt at 200MPa (H2O). American Mineralogist, 80(10-11): 1229–1251.
[] Inger S, Harris N. 1993. Geochemical constraints on leucogranite magmatism in the Langtang Valley, Nepal Himalaya. Journal of Petrology, 34(2): 345–368. DOI:10.1093/petrology/34.2.345
[] Jessup MJ, Cottle JM. 2010. Progression from south-directed extrusion to Orogen-parallel extension in the southern margin of the Tibetan Plateau, Mount Everest Region, Tibet. The Journal of Geology, 118(5): 467–486. DOI:10.1086/655011
[] Jolivet M, Brunel M, Seward D, Xu Z, Yang J, Malavieille J, Roger F, Leyreloup A, Aenand N, Wu C. 2003. Neogene extension and volcanism in the Kunlun fault zone, northern Tibet: New constraints on the age of the Kunlun fault. Tectonics, 22(5): 1052.
[] Kaneko Y, Katayama I, Yamamoto H, Misawa K, Ishikawa M, Rehman HU, Kausar AB, Shiraishi K. 2003. Timing of Himalayan ultrahigh-pressure metamorphism: Sinking rate and subduction angle of the Indian continental crust beneath Asia. Journal of Metamorphic Geology, 21(6): 589–599. DOI:10.1046/j.1525-1314.2003.00466.x
[] King J, Harris N, Argles T, Parrish R, Zhang HF. 2011. Contribution of crustal anatexis to the tectonic evolution of Indian crust beneath southern Tibet. Geological Society of America Bulletin, 123(1-2): 218–239. DOI:10.1130/B30085.1
[] Knesel KM, Davidson JP. 2002. Insights into collisional magmatism from isotopic fingerprints of melting reactions. Science, 296(5576): 2206–2208. DOI:10.1126/science.1070622
[] Kohn MJ. 2014. Himalayan metamorphism and its tectonic implications. Annual Review of Earth and Planetary Sciences, 42(1): 381–419. DOI:10.1146/annurev-earth-060313-055005
[] Larson KP, Godin L, Davis WJ, Davis DW. 2010. Out-of-sequence deformation and expansion of the Himalayan orogenic wedge: Insight from the Changgo culmination, south central Tibet. Tectonics, 29(4): TC4013.
[] Le Breton N, Thompson AB. 1988. Fluid-absent (dehydration) melting of biotite in metapelites in the early stages of crustal anatexis. Contributions to Mineralogy and Petrology, 99(2): 226–237. DOI:10.1007/BF00371463
[] Le Fort P. 1981. Manaslu leucogranite: A collision signature of the Himalaya: A model for its genesis and emplacement. Journal of Geophysical Research: Solid Earth, 86(B11): 10545–10568. DOI:10.1029/JB086iB11p10545
[] Le Fort P, Cuney M, Deniel C, France-Lanord C, Sheppard SMF, Upreti BN, Vidal P. 1987. Crustal generation of the Himalayan leucogranites. Tectonophysics, 134(1-3): 39–57. DOI:10.1016/0040-1951(87)90248-4
[] Lee J, 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
[] Linnen RL, Keppler H. 1997. Columbite solubility in granitic melts: Consequences for the enrichment and fractionation of Nb and Ta in the Earth's crust. Contributions to Mineralogy and Petrology, 128(2): 213–227.
[] Linnen RL. 1998. The solubility of Nb-Ta-Zr-Hf-W in granitic melts with Li and Li+F: Constraints for mineralization in rare metal granites and pegmatites. Economic Geology, 93(7): 1013–1025. DOI:10.2113/gsecongeo.93.7.1013
[] Linnen RL, Keppler H. 2002. Melt composition control of Zr/Hf fractionation in magmatic processes. Geochimica et Cosmochimica Acta, 66(18): 3293–3301. DOI:10.1016/S0016-7037(02)00924-9
[] Linnen RL. 2005. The effect of water on accessory phase solubility in subaluminous and peralkaline granitic melts. Lithos, 80(1-4): 267–280. DOI:10.1016/j.lithos.2004.04.060
[] Linnen RL, Van Lichtervelde M, Černy P. 2012. Granitic pegmatites as sources of strategic metals. Elements, 8(4): 275–280. DOI:10.2113/gselements.8.4.275
[] Liu ZC, Wu FY, Ji WQ, Wang JG, Liu CZ. 2014. Petrogenesis of the Ramba leucogranite in the Tethyan Himalaya and constraints on the channel flow model. Lithos, 208-209: 118–136. DOI:10.1016/j.lithos.2014.08.022
[] Liu ZC, Wu FY, Ding L, Liu XC, Wang JG, Ji WQ. 2016. Highly fractionated Late Eocene (~35Ma) leucogranite in the Xiaru Dome, Tethyan Himalaya, South Tibet. Lithos, 240-243: 337–354. DOI:10.1016/j.lithos.2015.11.026
[] London D, Wolf MB, Morgan Ⅵ GB, Garrido MG. 1999. Experimental silicate-phosphate equilibria in peraluminous granitic magmas, with a case study of the Alburquerque Batholith at Tres Arroyos, Badajoz, Spain. Journal of Petrology, 40(1): 215–240. DOI:10.1093/petroj/40.1.215
[] López S, Castro A. 2001. Determination of the fluid-absent solidus and supersolidus phase relationships of MORB-derived amphibolites in the range 4~14kbar. American Mineralogist, 86(11-12): 1396–1403. DOI:10.2138/am-2001-11-1208
[] Miller C, Thöni M, Frank W, Grasemann B, Klötzli U, Guntli P, Draganits E. 2001. The Early Palaeozoic magmatic event in the Northwest Himalaya, India: Source, tectonic setting and age of emplacement. Geological Magazine, 138(3): 237–251.
[] Miller CF, Stoddard EF. 1981. The role of manganese in the paragenesis of magmatic garnet: An example from the Old Woman-Piute Range, California. The Journal of Geology, 89(2): 233–246. DOI:10.1086/628582
[] Miller CF, Mittlefehldt DW. 1982. Depletion of light rare-earth elements in felsic magmas. Geology, 10(3): 129–133. DOI:10.1130/0091-7613(1982)10<129:DOLREI>2.0.CO;2
[] Miller CF, Mittlefehldt DW. 1984. Extreme fractionation in felsic magma chambers: A product of liquid-state diffusion or fractional crystallization?. Earth and Planetary Science Letters, 68(1): 151–158. DOI:10.1016/0012-821X(84)90147-X
[] Miller CF. 1985. Are strongly peraluminous magmas derived from pelitic sedimentary sources?. The Journal of Geolog, 93(6): 673–689. DOI:10.1086/628995
[] Miller CF, McDowell SM, Mapes RW. 2003. Hot and cold granites? Implications of zircon saturation temperatures and preservation of inheritance. Geology, 31(6): 529–532. DOI:10.1130/0091-7613(2003)031<0529:HACGIO>2.0.CO;2
[] Mitsuishi M, Wallis SR, Aoya M, Lee J, Wang Y. 2012. E-W extension at 19Ma in the Kung Co area, S. Tibet: Evidence for contemporaneous E-W and N-S extension in the Himalayan orogen. Earth and Planetary Science Letters, 325-326: 10–20.
[] Mittlefehldt DW, Miller CF. 1983. Geochemistry of the Sweetwater Wash Pluton, California: Implications for "anomalous" trace element behavior during differentiation of felsic magmas. Geochimica et Cosmochimica Acta, 47(1): 109–124. DOI:10.1016/0016-7037(83)90095-9
[] Montel JM. 1993. A model for monazite/melt equilibrium and application to the generation of granitic magmas. Chemical Geology, 110(1-3): 127–146. DOI:10.1016/0009-2541(93)90250-M
[] Nábělek J, Hetényi G, Vergne J, Sapkota S, Kafle B, Jiang M, Su HP, Chen J, Huang BS, Team HC. 2009. Underplating in the Himalaya-Tibet collision zone revealed by the Hi-CLIMB experiment. Science, 325(5946): 1371–1374. DOI:10.1126/science.1167719
[] Nabelek PI, Glascock MD. 1995. REE-depleted leucogranites, Black Hills, South Dakota: A consequence of disequilibrium melting of monazite-bearing schists. Journal of Petrology, 36(4): 1055–1071. DOI:10.1093/petrology/36.4.1055
[] Nabelek PI, Liu M. 2004. Petrologic and thermal constraints on the origin of leucogranites in collisional orogens. Earth and Environmental Science Transactions of the Royal Society of Edinburgh, 95(1-2): 73–85. DOI:10.1017/S0263593300000936
[] Noble SR, Searle MP. 1995. Age of crustal melting and leucogranite formation from U-Pb zircon and monazite dating in the western Himalaya, Zanskar, India. Geology, 23(12): 1135–1138. DOI:10.1130/0091-7613(1995)023<1135:AOCMAL>2.3.CO;2
[] Perini G, Cesare B, Gömez-Pugnaire MT, Ghezzi L, Tommasini S. 2009. Armouring effect on Sr-Nd isotopes during disequilibrium crustal melting: The case study of frozen migmatites from El Hoyazo and Mazarrón, SE Spain. European Journal of Mineralogy, 21(1): 117–131. DOI:10.1127/0935-1221/2009/0021-1882
[] Pichavant M, Montel JM, Richard LR. 1992. Apatite solubility in peraluminous liquids: Experimental data and an extension of the Harrison-Watson model. Geochimica et Cosmochimica Acta, 56(10): 3855–3861. DOI:10.1016/0016-7037(92)90178-L
[] Pognante U, Castelli D, Benna P, Genovese G, Oberli F, Meier M, Tonarini S. 1990. The crystalline units of the High Himalayas in the Lahul-Zanskar region (Northwest India): Metamorphic-tectonic history and geochronology of the collided and imbricated Indian plate. Geological Magazine, 127(2): 101–116. DOI:10.1017/S0016756800013807
[] Prince C, Harris N, Vance D. 2001. Fluid-enhanced melting during prograde metamorphism. Journal of the Geological Society, 158(2): 233–241. DOI:10.1144/jgs.158.2.233
[] Qi XX, Zeng LS, Meng XJ, Xu ZQ, Li TF. 2008. Zircon SHRIMP U-Pb dating for Dala granite in the Tethyan Himalaya and its geochemical implication. Acta Petrologica Sinica, 24(7): 1501–1508.
[] Quigley MC, Yu LJ, Gregory C, Corvino A, Sandiford M, Wilson CJL, Liu XH. 2008. U-Pb SHRIMP zircon geochronology and T-t-d history of the Kampa Dome, southern Tibet. Tectonophysics, 446(1-4): 97–113. DOI:10.1016/j.tecto.2007.11.004
[] Rapp RP, Watson EB. 1986. Monazite solubility and dissolution kinetics: Implications for the thorium and light rare earth chemistry of felsic magmas. Contributions to Mineralogy and Petrology, 94(3): 304–316. DOI:10.1007/BF00371439
[] Rapp RP, 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 MH. 2012. Partitioning of Sr, Ba, Rb, Y, and LREE between alkali feldspar and peraluminous silicic magma. American Mineralogist, 88(8-9): 1290–1303.
[] Richards A, Parrish R, Harris N, Argles T, Zhang L. 2006. Correlation of lithotectonic units across the eastern Himalaya, Bhutan. Geology, 34(5): 341–344. DOI:10.1130/G22169.1
[] Ryerson FJ, Watson EB. 1987. Rutile saturation in magmas: Implications for Ti-Nb-Ta depletion in island-arc basalts. Earth and Planetary Science Letters, 86(2-4): 225–239. DOI:10.1016/0012-821X(87)90223-8
[] Saleeby J, Ducea M, Clemens-Knott D. 2003. Production and loss of high-density batholithic root, southern Sierra Nevada, California. Tectonics, 22(6): 1064.
[] Schärer U, Xu RH, Allègre C. 1986. U-(Th)-Pb systematics and ages of Himalayan leucogranites, south Tibet. Earth and Planetary Science Letters, 77(1): 35–48. DOI:10.1016/0012-821X(86)90130-5
[] Searle MP, Parrish RR, Hodges KV, Hurford A, Ayres MW, Whitehouse MJ. 1997. Shisha Pangma leucogranite, South Tibetan Himalaya: Field relations, geochemistry, age, origin, and emplacement. The Journal of Geology, 105(3): 295–318. DOI:10.1086/515924
[] Searle MP, Noble SR, Hurford AJ, Rex DC. 1999. Age of crustal melting, emplacement and exhumation history of the Shivling leucogranite, Garhwal Himalaya. Geological Magazine, 136(5): 513–525. DOI:10.1017/S0016756899002885
[] Searle MP, Parrish RR, Thow AV, Noble S, Phillips RJ, Waters DJ. 2010. Anatomy, age and evolution of a collisional mountain belt: The baltoro granite batholith and Karakoram metamorphic complex, Pakistani Karakoram. Journal of the Geological Society, 167(1): 183–202. DOI:10.1144/0016-76492009-043
[] Shearer CK, Papike JJ. 1987. Mineralogical and chemical evolution of a rare-element granite-pegmatite system: Harney peak granite, black hills, south Dakota. Geochimica et Cosmochimica Acta, 51(3): 473–486. DOI:10.1016/0016-7037(87)90062-7
[] Simpson RL, Parrish RR, Searle MP, Waters DJ. 2000. Two episodes of monazite crystallization during metamorphism and crustal melting in the Everest region of the Nepalese Himalaya. Geology, 28(5): 403–406. DOI:10.1130/0091-7613(2000)28<403:TEOMCD>2.0.CO;2
[] Singh J, Johannes W. 1996. Dehydration melting of tonalites. Part Ⅱ. Composition of melts and solids. Contributions to Mineralogy and Petrology, 125(1): 26–44. DOI:10.1007/s004100050204
[] Singh S, Barley ME, Brown SJ, Brown SJ, Jain AK, Manickavasagam RM. 2002. SHRIMP U-Pb in zircon geochronology of the Chor granitoid: Evidence for Neoproterozoic magmatism in the Lesser Himalayan granite belt of NW India. Precambrian Research, 18(3-4): 285–292.
[] Stevens G, Clemens JD, Droop GTR. 1997. Melt production during granulite-facies anatexis: Experimental data from "primitive" metasedimentary protoliths. Contributions to Mineralogy and Petrology, 128(4): 352–370. DOI:10.1007/s004100050314
[] Stevens G, Villaros A, Moyen JF. 2007. Selective peritectic garnet entrainment as the origin of geochemical diversity in S-type granites. Geology, 35(1): 9–12. DOI:10.1130/G22959A.1
[] Sun SS, 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. DOI:10.1144/GSL.SP.1989.042.01.19
[] Sylvester PJ. 1998. Post-collisional strongly peraluminous granites. Lithos, 45(1-4): 29–44. DOI:10.1016/S0024-4937(98)00024-3
[] Tang M, Wang XL, Shu XJ, Wang D, Yang T, Gopon P. 2014. Hafnium isotopic heterogeneity in zircons from granitic rocks: Geochemical evaluation and modeling of "zircon effect" in crustal anatexis. Earth and Planetary Science Letters, 389: 188–199. DOI:10.1016/j.epsl.2013.12.036
[] Taylor JR, Wall VJ. 1992. The behavior of tin in granitoid magmas. Economic Geology, 87(2): 403–420. DOI:10.2113/gsecongeo.87.2.403
[] Thompson AB. 2001. Clockwise P-T paths for crustal melting and H2O recycling in granite source regions and migmatite terrains. Lithos, 56(1): 33–45. DOI:10.1016/S0024-4937(00)00058-X
[] Vance D, Mahar E. 1998. Pressure-temperature paths from P-T pseudosections and zoned garnets: Potential, limitations and examples from the Zanskar Himalaya, NW India. Contributions to Mineralogy and Petrology, 132(3): 225–245. DOI:10.1007/s004100050419
[] Vance D, Harris N. 1999. Timing of prograde metamorphism in the Zanskar Himalaya. Geology, 27(5): 395–398. DOI:10.1130/0091-7613(1999)027<0395:TOPMIT>2.3.CO;2
[] Vannay JC, Grasemann B. 2001. Himalayan inverted metamorphism and syn-convergence extension as a consequence of a general shear extrusion. Geological Magazine, 138(3): 253–276.
[] Vidal P, Cocherie A, Le Fort P. 1982. Geochemical investigations of the origin of the Manaslu leucogranite (Himalaya, Nepal). Geochimica et Cosmochimica Acta, 46(11): 2279–2292. DOI:10.1016/0016-7037(82)90201-0
[] Vielzeuf D, Holloway JR. 1988. Experimental determination of the fluid-absent melting relations in the pelitic system. Contributions to Mineralogy and Petrology, 98(3): 257–276. DOI:10.1007/BF00375178
[] Vielzeuf D, Montel JM. 1994. Partial melting of metagreywackes. Part Ⅰ. Fluid-absent experiments and phase relationships. Contributions to Mineralogy and Petrology, 117(4): 375–393.
[] Wang XX, Zhang JJ, Santosh M, Liu J, Yan SY, Guo L. 2012. Andean-type orogeny in the Himalayas of South Tibet: Implications for early Paleozoic tectonics along the Indian margin of Gondwana. Lithos, 154: 248–262. DOI:10.1016/j.lithos.2012.07.011
[] Wang XX, Zhang JJ, Yan SY, Liu J, Wang JM, Guo L. 2015. Age and formation mechanism of the Qiamangba two-mica granite in northern Himalaya. Chinese Journal of Geology, 50(3): 708–727.
[] Wang XX, Zhang JJ, Yan SY, Liu J. 2016. Age and geochemistry of the Cuona leucogranite in southern Tibet and its geological implications. Geological Bulletin of China, 35(1): 91–103.
[] Wang YH, Zeng LS, Gao LE, Zhang LF, 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.
[] Wang YH, Zeng LS, Gao LE, Hou KJ, Guo CL, Zhang LF, Wang W, Sun HY. 2017. Neoproterozoic magmatism in eastern Himalayan Terrane. Science Bulletin, 62(6): 415–424. DOI:10.1016/j.scib.2017.02.003
[] Watson EB. 1979. Zircon saturation in felsic liquids: Experimental results and applications to trace element geochemistry. Contributions to Mineralogy and Petrology, 70(4): 407–419. DOI:10.1007/BF00371047
[] Watson EB, Harrison TM. 1983. Zircon saturation revisited: Temperature and composition effects in a variety of crustal magma types. Earth and Planetary Science Letters, 64(2): 295–304. DOI:10.1016/0012-821X(83)90211-X
[] Watson EB, Harrison TM. 1984. Accessory minerals and the geochemical evolution of crustal magmatic systems: A summary and prospectus of experimental approaches. Physics of the Earth and Planetary Interiors, 35(1-3): 19–30. DOI:10.1016/0031-9201(84)90031-1
[] Watson EB, Vicenzi EP, Rapp RP. 1989. Inclusion/host relations involving accessory minerals in high-grade metamorphic and anatectic rocks. Contributions to Mineralogy and Petrology, 101(2): 220–231. DOI:10.1007/BF00375308
[] Watson EB. 1996. Dissolution, growth and survival of zircons during crustal fusion: Kinetic principals, geological models and implications for isotopic inheritance. Earth and Environmental Science Transactions of the Royal Society of Edinburgh, 87(1-2): 43–56. DOI:10.1017/S0263593300006465
[] Watt GR, Harley SL. 1993. Accessory phase controls on the geochemistry of crustal melts and restites produced during water-undersaturated partial melting. Contributions to Mineralogy and Petrology, 114(4): 550–566. DOI:10.1007/BF00321759
[] Weinberg RF, Hasalová P. 2015. Water-fluxed melting of the continental crust: A review. Lithos, 212-215: 158–188. DOI:10.1016/j.lithos.2014.08.021
[] Williams H, Turner S, Kelley S, Harris N. 2001. Age and composition of dikes in Southern Tibet: New constraints on the timing of east-west extension and its relationship to postcollisional volcanism. Geology, 29(4): 339–342. DOI:10.1130/0091-7613(2001)029<0339:AACODI>2.0.CO;2
[] Wolf MB, Wyllie PJ. 1991. Dehydration-melting of solid amphibolite at 10kbar: Textural development, liquid interconnectivity and applications to the segregation of magmas. Contributions to Mineralogy and Petrology, 44(3-4): 151–179. DOI:10.1007/BF01166961
[] Wolf MB, Wyllie PJ. 1993. Garnet growth during amphibolite anatexis: Implications of a garnetiferous restite. Journal of Geology, 101(3): 357–373. DOI:10.1086/648229
[] Wolf MB, London D. 1994. Apatite dissolution into peraluminous haplogranitic melts: An experimental study of solubilities and mechanisms. Geochimica et Cosmochimica Acta, 58(19): 4127–4145. DOI:10.1016/0016-7037(94)90269-0
[] Wolf MB, Wyllie PJ. 1994. Dehydration-melting of amphibolite at 10kbar: The effects of temperature and time. Contributions to Mineralogy and Petrology, 115(4): 369–383. DOI:10.1007/BF00320972
[] Wolf MB, London D. 1995. Incongruent dissolution of REE-and Sr-rich apatite in peraluminous granitic liquids: Differential apatite, monazite, and xenotime solubilities during anatexis. American Mineralogist, 80(7-8): 765–775. DOI:10.2138/am-1995-7-814
[] Wu FY, Liu ZC, Liu XC, Ji WQ. 2015. Himalayan leucogranite: Petrogenesis and implications to orogenesis and plateau uplift. Acta Petrologica Sinica, 31(1): 1–36.
[] Wyllie PJ, Wolf MB. 1993. Amphibolite dehydration-melting: Sorting out the solidus. In: Prichard HM, Alabaster T, Harris NBW and Neary CR (eds.). Magmatic Processes and Plate Tectonics. Geological Society, London, Special Publications, 76(1): 405–416.
[] Xie KJ, Zeng LS, Liu J, Gao LE. 2010. Late-Eocene Dala adakitic granite, southern Tibet and geological implications. Acta Petrologica Sinica, 26(4): 1016–1026.
[] Xu ZQ, Yang JS, Liang FH, Qi XX, Liu FL, Zeng LS, Liu DY, Li HB, Wu CL, Shi RD, Chen SY. 2005. Pan-African and early Paleozoic orogenic events in the Himalaya terrane: Inference from SHRIMP U-Pb zircon ages. Acta Petrologica Sinica, 21(1): 1–12.
[] Yang XS, Jin ZM. 2001. Studies on Rb-Sr and Sm-Nd isotope of Yadong leucogranite in Tibet: Constraint on its age and source material. Geological Review, 47(3): 294–300.
[] Yang XY, Zhang JJ, Qi GW, Wang DC, Guo L, Li PY, Liu J. 2009. Structure and deformation around the Gyirong basin, North Himalaya, and onset of the South Tibetan detachment system. Science in China (Series D), 52(8): 1046–1058. DOI:10.1007/s11430-009-0111-2
[] Yin A, 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
[] Yin A. 2006. Cenozoic tectonic evolution of the Himalayan orogen as constrained by along-strike variation of structural geometry, exhumation history, and foreland sedimentation. Earth Science Reviews, 76(1): 1–131.
[] Yu JJ, Zeng LS, Liu J, Gao LE, Xie KJ. 2011. Early Miocene leucogranites in Dinggye area, southern Tibet: Formation mechanism and tectonic implications. Acta Petrologica Sinica, 27(7): 1961–1972.
[] Zeng LS, Asimow PD, Saleeby JB. 2005a. Coupling of anatectic reactions and dissolution of accessory phases and the Sr and Nd isotope systematics of anatectic melts from a metasedimentary source. Geochimica et Cosmochimica Acta, 69(14): 3671–3682. DOI:10.1016/j.gca.2005.02.035
[] Zeng LS, Saleeby JB, Asimow P. 2005b. Nd isotope disequilibrium during crustal anatexis: A record from the Goat Ranch migmatite complex, southern Sierra Nevada batholith, California. Geology, 33(1): 53–56. DOI:10.1130/G20831.1
[] Zeng LS, Saleeby JB, Ducea M. 2005c. Geochemical characteristics of crustal anatexis during the formation of migmatite at the Southern Sierra Nevada, California. Contributions to Mineralogy and Petrology, 150(4): 386–402. DOI:10.1007/s00410-005-0010-2
[] Zeng LS, Liu J, Saleeby JB. 2006. Root zone processes for the formation and evolution of a large-scale granitic batholith: Drip structures, potassic volcanism and Earth surface processes. Geological Bulletin of China, 25(11): 1257–1273.
[] Zeng LS, Liu J, Gao LE, Xie KJ, 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, Zeng-Liu J. 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, Dong CY, Tang SH. 2012. High-pressure melting of metapelite and the formation of Ca-rich granitic melts in the Namche Barwa Massif, Southern Tibet. Gondwana Research, 21(1): 138–151. DOI:10.1016/j.gr.2011.07.023
[] Zeng LS, Gao LE, Tang SH, Hou KJ, Guo CL, Hu GY. 2015. Eocene magmatism in the Tethyan Himalaya, southern Tibet. In: Jenkin GRT and Lusty PAJ (eds.). Ore Deposits in an Evolving Earth. Geological Society, London, Special Publications, 412(1): 287–316.
[] Zhang HF, Harris N, Parrish R, Kelley S, Zhang L, Rogers N, Argles T, King J. 2004a. 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 HF, Harris N, Parrish R, Zhang L, Zhao ZD. 2004b. U-Pb ages of Kude and Sajia leucogranites in Sajia dome from North Himalaya and their geological implications. Chinese Science Bulletin, 49(19): 2087–2092. DOI:10.1360/04wd0198
[] Zhang JJ, Yang XY, Qi GW, Wang DC. 2011. Geochronology of the Malashan dome and its application in formation of the Southern Tibet detachment system (STDS) and Northern Himalayan gneiss domes (NHGD). Acta Petrologica Sinica, 27(12): 3535–3544.
[] Zhang JJ, Santosh M, Wang XX, Guo L, Yang XY, Zhang B. 2012. Tectonics of the northern Himalaya since the India-Asia collision. Gondwana Research, 21(4): 939–960. DOI:10.1016/j.gr.2011.11.004
[] Zhang ZM, Wang JL, Shen K, Shi C. 2008. Paleozoic circum-Gondwana orogens: Petrology and geochronology of the Namche Barwa Complex in the eastern Himalayan syntaxis, Tibet. Acta Petrologica Sinica, 24(7): 1627–1637.
[] 高利娥, 曾令森, 刘静, 谢克家. 2009. 藏南也拉香波早渐新世富钠过铝质淡色花岗岩的成因机制及其构造动力学意义. 岩石学报, 25(9): 2289–2302.
[] 高利娥, 曾令森, 胡古月. 2010. 藏南确当地区高Sr/Y比值二云母花岗岩的形成机制及其构造动力学意义. 地质通报, 29(2-3): 214–226.
[] 高利娥, 曾令森, 石卫刚, 陈振宇, 胡明月, 孙东阳. 2012. 喜马拉雅造山带新生代花岗岩中两类石榴石的地球化学特征及其在地壳深熔作用中的意义. 岩石学报, 28(9): 2963–2980.
[] 高利娥, 曾令森, 许志琴, 王莉. 2015. 喜马拉雅造山带加里东期构造作用:以马拉山-吉隆构造带为例. 岩石学报, 31(5): 1200–1218.
[] 高利娥, 曾令森, 王莉, 侯可军, 高家昊, 尚振. 2016. 喜马拉雅碰撞造山带不同类型部分熔融作用的时限及其构造动力学意义. 地质学报, 90(11): 3039–3059. DOI:10.3969/j.issn.0001-5717.2016.11.006
[] 高利娥, 曾令森, 高家昊, 赵令浩, 侯可军, 唐索寒. 2017. 藏南拿日雍措片穹窿中新世淡色花岗岩的形成过程: 变泥质岩部分熔融与分离结晶作用. 岩石学报, 33(7): 出版中
[] 郭春丽, 曾令森, 高利娥, 苏红中, 马星华, 尹冰. 2017. 高分异花岗岩中矿物和全岩地球化学找矿标志——以河田与西华山花岗岩对比为例. 地质学报, (8): 出版中
[] 胡古月, 曾令森, 戚学祥, 侯可军, 高利娥. 2011. 藏南特提斯喜马拉雅带始新世隆子-恰嘎次火山岩区:雅拉香波二云母花岗岩的高位岩浆体系. 岩石学报, 27(11): 3308–3318.
[] 胡古月, 曾令森, 陈翰, 高利娥. 2016. 喜马拉雅造山带晚中新世麻迦淡色花岗岩的构建机制. 地质学报, 90(8): 1737–1754.
[] 戚学祥, 曾令森, 孟祥金, 许志琴, 李天福. 2008. 特提斯喜马拉雅打拉花岗岩的锆石SHRIMP U-Pb定年及其地质意义. 岩石学报, 24(7): 1501–1508.
[] 王晓先, 张进江, 闫淑玉, 刘江, 王佳敏, 郭磊. 2015. 北喜马拉雅恰芒巴二云母花岗岩的年龄及形成机制. 地质科学, 50(3): 708–727.
[] 王晓先, 张进江, 闫淑玉, 刘江. 2016. 藏南错那淡色花岗岩LA-MC-ICP-MS锆石U-Pb年龄、岩石地球化学及其地质意义. 地质通报, 35(1): 91–103.
[] 王誉桦, 曾令森, 高利娥, 张立飞, 侯可军. 2014. 喜马拉雅造山带东构造结拉布拉多期和格林威尔期造山作用的记录. 岩石学报, 30(8): 2241–2252.
[] 吴福元, 刘志超, 刘小驰, 纪伟强. 2015. 喜马拉雅淡色花岗岩. 岩石学报, 31(1): 1–36.
[] 谢克家, 曾令森, 刘静, 高利娥. 2010. 西藏南部晚始新世打拉埃达克质花岗岩及其构造动力学意义. 岩石学报, 26(4): 1016–1026.
[] 许志琴, 杨经绥, 梁凤华, 戚学祥, 刘福来, 曾令森, 刘敦一, 李海兵, 吴才来, 史仁灯, 陈松永. 2005. 喜马拉雅地体的泛非-早古生代造山事件年龄记录. 岩石学报, 21(1): 1–12.
[] 杨晓松, 金振民. 2001. 西藏亚东淡色花岗岩Rb-Sr和Sm-Nd同位素研究——关于其年龄和源岩的证据. 地质论评, 47(3): 294–300.
[] 于俊杰, 曾令森, 刘静, 高利娥, 谢克家. 2011. 藏南定结地区早中新世淡色花岗岩的形成机制及其构造动力学意义. 岩石学报, 27(7): 1961–1972.
[] 曾令森, 刘静, SaleebyJB. 2006. 大型花岗岩基形成和演化的深部动力学过程:滴水构造、钾质火山作用与地表地质过程. 地质通报, 25(11): 1257–1273. DOI:10.3969/j.issn.1671-2552.2006.11.003
[] 张进江, 杨雄英, 戚国伟, 王德朝. 2011. 马拉山穹窿的活动时限及其在藏南拆离系——北喜马拉雅片麻岩穹窿形成机制的应用. 岩石学报, 27(12): 3535–3544.
[] 张泽明, 王金丽, 沈昆, 石超. 2008. 环东冈瓦纳大陆周缘的古生代造山作用:东喜马拉雅构造结南迦巴瓦岩群的岩石学和年代学证据. 岩石学报, 24(7): 1627–1637.