岩石学报  2019, Vol. 35 Issue (3): 665-686, doi: 10.18654/1000-0569/2019.03.04   PDF    
班公湖-怒江缝合带西段阿翁错复式岩体的岩浆混合成因:地球化学、年代学和暗色微粒包体证据
雷传扬1,2,3 , 唐菊兴4 , 李威2 , 袁华云2 , 张伟2 , 尹滔2 , 裴亚伦2 , 郎兴海1     
1. 成都理工大学地球科学学院, 成都 610059;
2. 四川省地质调查院, 成都 610081;
3. 四川省地矿局四○五地质队, 都江堰 611830;
4. 中国地质科学院矿产资源研究所, 自然资源部成矿作用与资源评价重点实验室, 北京 100037
摘要:阿翁错复式岩体位于班公湖-怒江缝合带西段,是班公湖-怒江特提斯洋俯冲消减,造山过程中岩浆响应的重要组成部分,以广泛发育暗色微粒包体和岩浆混合、不协调现象为特征。本文以阿翁错复式岩体为研究对象,对寄主岩和暗色微粒包体开展了系统的地质学、地球化学和锆石U-Pb年代学研究,探讨了阿翁错复式岩体的岩浆混合成因。暗色微粒包体塑性变形特征明显,与寄主岩呈截然或渐变接触,偶见反向脉发育,包体具细-中粗粒结构,含斜长石、钾长石、角闪石、暗色镶边石英等斑晶,偶见角闪石斑晶横跨包体和寄主岩,在包体及包体周围寄主岩中见长柱状斜长石、角闪石和针状磷灰石等结构特征,表明暗色微粒包体为岩浆混合作用的产物。寄主岩与包体均为准铝质、钙碱性-高钾钙碱性系列岩石,主要氧化物含量在Harker图解上具有良好的线性关系,稀土元素配分曲线图和微量元素蛛网图具有高度一致性,表明二者具有强烈的地球化学亲源关系,且经历了相似的岩浆演化过程。寄主岩和暗色微粒包体的成岩年龄分别为109.1±1.0Ma和107.4±0.7Ma,岩浆混合作用发生在早白垩世晚期,处于班公湖-怒江特提斯洋由弧-陆碰撞向陆陆碰撞的转换阶段即软碰撞阶段。研究表明,在班公湖-怒江特提斯洋向北向羌塘地块之下俯冲的背景下,洋壳脱水,引起上覆地幔楔发生部分熔融,形成镁铁质岩浆,镁铁质岩浆向上运移,并底侵于壳-幔边界,引发下地壳物质发生部分熔融,形成长英质岩浆,当镁铁质岩浆从底部注入长英质岩浆房时,镁铁质岩浆快速冷凝,形成部分色率高、粒度细,具冷凝边的包体,与寄主岩呈截然型接触,随着端元岩浆之间的温差逐渐降低,包体色率降低,粒度变大,与寄主岩呈渐变过渡。
关键词: 地球化学     年代学     暗色微粒包体     岩浆混合作用     阿翁错复式岩体     班公湖-怒江缝合带西段     西藏    
Genesis of magma mixing and mingling of the Awengcuo composite plutons in western segment of Bangongco-Nujiang Suture Zone: Evidence from geochemistry, geochronology and mafic microgranular enclaves (MMEs)
LEI ChuanYang1,2,3, TANG JuXing4, LI Wei2, YUAN HuanYun2, ZHANG Wei2, YIN Tao2, PEI YaLun2, LANG XingHai1     
1. College of Earth Sciences, Chengdu University of Technology, Chengdu 610059, China;
2. Geological Survey of Sichuan Province, Chengdu 610081, China;
3. 405 Geological Teem, Sichuan Bureau of Geology and Mineral Resources, Dujianyan 611830, China;
4. MNR Key Laboratory of Metallogeny and Mineral Assessment, Institute of Mineral Resources, CAGS, Beijing 100037, China
Abstract: Located at the western part of Bangong Co-Nujiang Suture Zone, the Awengcuo composite plutons are the important part result in the magmatic response during the subduction of the Bangongco-Nujiang Tethys Ocean, and is featured by widely developed mafic microgranular enclaves, the magma mingling and the magma incoordination phenomenon. Taking Awengcuo composite plutons as the research object, this paper conducts an in-depth study on the geological characteristics and geochronology of mafic microgranular enclave and its formation mechanism. Mafic microgranular enclave has obvious plastic deformation features, and has sharp or gradual contact with the host rock, presenting still a little reverse pulse development. With the structure of fine to medium-coarse, the enclave contains many phenocrysts such as plagioclase, potash feldspar, hornblende, and dark rimmed quartz. Besides, plagioclase phenocryst is occasionally found across the enclave and the host rock, and long columnar plagioclase, hornblende, acicular apatite are also existed in the enclave and the host rock. All of them show that mafic microgranular enclave is produced by magma mixing. Host rock and MMEs are transitional rock such as quasi-aluminum, high-calcium alkaline-potassium kaolin series, and the main oxide content in the Harker diagram has a good linear relationship. Rare earth element distribution curve and normalized spidergram for trace element primitive mantle are with high degree of consistency. It shows that there is a strong geophysical relationship between the two groups, and experience the similar magmatic evolution process. A good linear relationship is showed between the host rock and inclusions in the denominator covariant map of the principle and trace elements. And their relationship is closely related to the magmatic mixing. Both of them contain a lot of large ion lithophile elements such as Rb, K, U, Th, and lack high field strenghth elements Nb, Ta, P, Ti. Dark Microgranular is low-silicon, but rich with Mg, Fe, Ca, The ratio of Nb/Ta is 17.6, indicating that the original magma of the inclusions may be derived from the mantle wedge by fluid of subduction zone. While the host rock is relatively high in silicon and with lower Mg, Fe and Ca, Nb/Ta ratio is 9.85, indicating that the host rock is the product of mantle-derived magma mixed with the molten lower crust magma, mainly with molten lower crust magma. The LA-ICP-MS single-zircons U-Pb dating demonstrates that the magmatic age is 109.1±1.0Ma for the host rock and 107.4±0.7Ma for the enclave. The two ages are basically the same, suggesting that the mixing event took place during the period of Early Cretaceous, during the soft collision phrase which the Bangong Co-Nujiang Tethys Ocean is transitional from the ocean-continent subduction to continent-continent subduction. Research shows that during Bangongco-Nujiang Tethyan Oceanic plate north-forward subduction to the Qiangtang block, the dehydration of oceanic crust leads to the partial melting of the overlain mantle wedge, and produces mafic magma. The mafic magma transfers upward and eventually intrudes into crust-mantle boundary, which causes the partial melting of the lower crustal materials and forms felsic magma. When the mafic magma flows at the bottom into the felsic magma chamber, it will condensates rapidly and become some enclave with high color ratio, fine granularity and condensation rim. Such enclave is in sharp contact with the host rock, while with the gradually declining temperature between end-member magmas, the color ratio of the enclave is becoming lower and its granularity bigger, then the enclave has a gradational contact with the host rock.
Key words: Geochemistry     Geochronology     Mafic microgranular enclaves     Magmatic mixing     Awengcuo composite plutons     Bangongco-Nujiang Suture     Tibet    

岩浆混合作用与岩浆分异作用、分离结晶作用、同化混染作用一样,是岩浆演化的重要表现形式,也是导致岩浆岩多样性的重要原因(Barbarin, 2005; 王玉往等, 2012),岩浆混合作用以寄主岩中广泛发育各种暗色微粒包体和不协调现象为显著特征(张旗等, 2007; 王德滋和谢磊, 2008; 莫宣学, 2011)。花岗岩不仅仅是地壳物质再循环作用的产物,而且还有幔源岩浆的参与,幔源岩浆为地壳物质部分熔融提供热源,同时幔源物质注入长英质岩浆房,发生混合作用(Bonin, 2004; Castro, 2013)。根据岩浆混合作用程度将其分为低级和高级两个阶段,前者主要表现为机械混合(mingling),形成各种不协调现象及暗色微粒包体,后者主要表现为化学混合(mixing),形成均一的岩浆,二者之间是连续的,相互包容的(张旗等, 2007; 王玉往等, 2012),岩浆混合的能干性主要受岩浆黏度和温度影响(张旗等, 2007)。岩浆混合作用具有普遍性,不仅发生在Ⅰ型花岗岩中,S型和A型花岗岩中也存在岩浆混合作用(王涛, 2000; Bonin, 2007; Shellnutt et al., 2010; Champion and Bultitude, 2013),而花岗岩中暗色微粒包体的存在是岩浆混合最直接的证据,包体能提供上地幔与地壳的直接信息,是人们了解岩浆混合作用方式、成岩过程温压条件、端元岩浆等方面信息不可或缺的研究对象(莫宣学, 2011; 王玉往等, 2012)。

班公湖-怒江缝合带(下文简称班-怒带)为特提斯洋闭合后的残留,其西段出露大量中酸性岩浆岩,并发育一系列与中酸性火山岩-浅成岩组合有关的斑岩型-浅成低温热液型铜多金属矿床,以多龙矿集区为典型代表(唐菊兴等, 2014a, b, 2016, 2017; Lin et al., 2017a, b, 2019; Song et al., 2018),对其研究在解决基础地质问题和查明矿产分布规律等方面都有重要的意义。随着近年来各类地质调查工作的全面开展,公益性基础研究的不断深入,商业性勘查的投入,加之分析测试精度提高,班-怒带西段已取得了不少有价值的研究成果,以唐菊兴研究员和秦克章研究员为代表的研究团队对该缝合带西段出露的各类岩浆岩进行了系统的地质学、地球化学和年代学研究,初步查明了岩浆岩的成因、构造背景及深部动力学过程,为班公湖-怒江特提斯洋(下文简称班-怒洋)的形成演化提供了科学依据(李金祥等, 2006; 秦克章等, 2006; 李金祥, 2008; 雷传扬等, 2012, 2018, 2019; Li et al., 2013; 李光明等, 2015; 韦少港等, 2016, 2017a, b; 丁帅等, 2017; 林彬等, 2017; 宋扬等, 2017; 张志等, 2017; 王勤等, 2018)。尽管前人对缝合带西段岩浆岩开展了深入的研究,但一些关键问题并未涉及或未开展深入的研究工作,如岩体的岩浆混合成因。阿翁错复式岩体位于班-怒带西段,是班-怒洋俯冲消减,造山过程中岩浆响应的重要组成部分,以广泛发育暗色微粒包体为特征,是研究岩浆混合成因的理想对象,本文在对阿翁错复式岩体和暗色微粒包体开展详细的野外地质现象和特征研究基础上,分别采集了寄主岩和包体的锆石U-Pb年龄和地球化学分析样,系统总结了阿翁错复式岩体岩浆混合作用的岩相学证据,精确厘定了寄主岩和包体的成岩时代,探讨了岩体的岩浆混合成因、构造背景及动力学机制,以期为阿翁错复式岩体的岩浆混合成因及其动力学过程提供来自地质学、地球化学、年代学和暗色微粒包体的新证据。

① 秦克章, 李光明, 张旗, 李金祥, 缪宇, 肖波, 张天平, 多吉, 李金高, 陆彦. 2006.西藏浅成低温金-银矿的成矿条件与可能产出区分析——从斑岩-浅成低温铜金成矿系统的角度.见:中国科学院地质与地球物理研究所科技与成果转化处编.北京:中国科学院地质与地球物理研究所2006年论文摘要集, 666-670

1 区域地质及岩体地质特征

班-怒带夹持于羌塘地块和拉萨地块之间,其北界为班公湖-康托-兹格塘错-安多-丁青雪拉山断裂,南界为日土-改则-尼玛-丁青断裂(宋扬等, 2013),又以丁青、改则为界,把缝合带分为东、中、西三段(强巴扎西等, 2016),缝合带以产出北西西向不连续分布的蛇绿岩残片为特征。缝合带西段存在三条俯冲带,北面的日土、狮泉河-改则-洞错俯冲带,南面的拉果错俯冲带(曲晓明等, 2009)。阿翁错复式岩体位于班-怒带西段。研究区地层整体较简单,从侏罗系到第四系均有出露。区域性构造形迹主要表现为北西向的压性(兼具扭性)断裂,具有多期次、多阶段的活动特征。岩浆活动主要集中在燕山期(图 1)。

图 1 研究区大地构造位置(a)及地质简图(b) ①阿尼玛卿-昆仑断裂;②羊湖-金沙江断裂;③双湖-澜沧江断裂;④班公湖-怒江断裂;⑤雅鲁藏布江断裂;1-第四系;2-渐新统丁青湖组上段;3-渐新统丁青湖组下段;4-古新统-始新统牛堡组上段;5-古新统-始新统牛堡组下段;6-上侏罗统-下白垩统沙木罗组下段;7-早白垩世似斑状黑云母二长花岗岩;8-早白垩世中粒黑云母二长花岗岩;9-早白垩世似斑状黑云母钾长花岗岩;10-早白垩世中粒黑云母钾长花岗岩;11-早白垩世似斑状钾长花岗岩;12-早白垩世中粒钾长花岗岩;13-早白垩世中粒黑云母花岗闪长岩;14-早白垩世石英闪长岩;15-闪长岩脉/辉绿岩脉;16-地质界线;17-角度不整合接触;18-渐变过渡;19-断层;20-地质剖面;21-包体长短轴比值;22薄片编号及采样位置;23-照片编号及拍摄位置;24-中基性岩体出露位置 Fig. 1 Tectonic position (a) and sketch geological map (b) of the study area

阿翁错复式岩体以岩基形式呈北西-南东向产出,主要出露在阿翁错、丁字杂沼、热帮乡一带,地表出露长度约20km,宽约5~8km,面积约100km2。复式岩体侵入到沙木罗组下段(J3K1s1)地层中,在岩体与地层的接触带能见到明显的角岩化、硅化现象,由于同化作用不彻底,在接触带及附近岩体中能见到少量沙木罗组地层的捕虏体,面积较小,多呈不规则形态产出。根据地质特征及岩石学特征将阿翁错复式岩体划分为以下8种类型:细粒石英闪长岩、中粒黑云母花岗闪长岩、中粒黑云母钾长花岗岩、似斑状黑云母钾长花岗岩、中粒钾长花岗岩、似斑状钾长花岗岩、中粒黑云母二长花岗岩和似斑状黑云母二长花岗岩。其中石英闪长岩出露面积较小,其余岩石出露面积相当。不同岩性之间呈侵入接触,同一种岩性不同结构之间呈渐变过渡接触(图 2)。在黑云母钾长花岗岩中见中基性岩体出露(图 3a),宽约10~15m,走向上延伸275m,面积较小,空间上与黑云母钾长花岗岩相伴生,呈截然型接触,在接触部位见大量暗色微粒包体成群分布(图 3b),接触带往外包体数量逐渐减少。

图 2 阿翁错复式岩体剖面图 1-第四系;2-含砾砂岩;3-板岩;4-似斑状黑云母二长花岗岩;5-中粒黑云母二长花岗岩;6-似斑状黑云母钾长花岗岩;7-中粒黑云母钾长花岗岩;8-似斑状钾长花岗岩;9-中粒钾长花岗岩;10-中粒黑云母花岗闪长岩;11-细粒石英闪长岩;12-暗色微粒包体;13-侵入接触界线/渐变过渡界线;14-地质代号.剖面位置见图 1 Fig. 2 Sections of the Awengcuo composite plutons

图 3 阿翁错复式岩体中暗色微粒包体的野外露头特征 (a)黑云母钾长花岗岩中出露的中基性岩体;(b)密集型暗色微粒包体;(c)云雾状暗色微粒包体;(d)不规则形态暗色微粒包体;(e)同一露头暗色微粒包体具不同方向;(f)细长条状定向包体.照片拍摄地点详见图 1,下同 Fig. 3 Field outcrop characteristics of mafic microgranular enclaves in the Awengcuo composite plutons
2 暗色微粒包体地质特征 2.1 野外地质特征

在阿翁错复式岩体各类型岩石中几乎都有暗色微粒包体分布,从石英闪长岩到黑云母二长花岗岩包体含量逐渐增加,石英闪长岩中仅偶见包体,含量约占0.5%,黑云母花岗闪长岩中包体含量约占1%~2%,黑云母钾长花岗岩和黑云母二长花岗岩中包体含量约占5%~8%,在黑云母钾长花岗岩与中基性岩体接触部位包体成群、成带分布,含量约占50%~60%,局部可高达90%。以上显示从早期到晚期,从中性岩石到酸性岩石包体的含量逐渐增加,但野外观察发现同一类型岩石中包体的分布是随机的,未见规律性。

暗色微粒包体直径悬殊明显,可分为四类:小型(< 10cm)、中型(10~30cm)、大型(30~50cm)、巨型(> 50cm)。据野外观测统计,整体以小型为主,但在黑云母钾长花岗岩与中基性岩体接触部位以巨型为主。包体形态各异,见细长条状、椭圆状、透镜状、云雾状以及不规则形态,以椭圆状和细长条状居多,其中细长条状包体定向性明显,其它类型包体未见明显定向(图 3b-f)。

采用长短轴法对阿翁错复式岩体中6个包体分布较多且定向性明显的露头进行了暗色微粒包体测量统计(表 1),为了测量结果更加可靠,每个露头统计数大于30(陈国超等, 2016)。从表中可以看出6个露头包体的长短轴比值为4.14~5.32,整体上较为接近,长轴方向以北西向为主,其中D5378点长轴方向为255°,D5209点长轴方向为10°,与主体方向稍有偏差。

表 1 阿翁错复试岩体中包体轴比和长轴方向 Table 1 Axial ratios and long-axis directions of the mafic microgranular enclaves developed in the Awengcuo composite plutons
2.2 岩相学特征

暗色微粒包体呈浅灰色-黑色,色率明显高于寄主岩,以细粒为主,见少量粗粒。野外观察发现暗色微粒包体粒度越细,其色率越高,色率与暗色矿物含量呈正相关,色率越高,表明岩浆混合作用程度越低(陈国超等, 2016)。根据色率将暗色微粒包体划分为浅灰色、深灰色、黑色三类,浅灰色包体呈中粒结构,深灰色包体呈中细粒结构,黑色包体呈细粒结构,以黑色细粒包体为主,但有的包体边部矿物粒度较细,往中部矿物粒度逐渐变大(图 4a),在包体中可以见到寄主岩的矿物成分(图 4b-d)。据镜下观察将包体分为黑云闪长岩包体,具细粒半自形粒状结构,矿物组成为斜长石(61%)、普通角闪石(29%)、黑云母(6%)及少量石英(1%),副矿物见榍石(2%)和少量磷灰石(图 4e);石英闪长岩包体,具细粒半自形粒状结构,矿物组成为斜长石(52%)、角闪石(32%)、石英(5%)、碱性长石(4%)、黑云母(4%),副矿物见磷灰石(图 4f);黑云石英二长闪长岩包体,具中细粒半自形粒状结构,矿物组成为斜长石(63%)、钾长石(10%)、石英(12%)、黑云母(11%)、角闪石(3%)副矿物见少量磷灰石和榍石(图 4g)。

图 4 暗色微粒包体及寄主岩的岩石学特征 (a)包体与寄主岩接触部位矿物粒度较细,往包体中部矿物粒度逐渐变大;(b)包体中见来自寄主岩的钾长石斑晶;(c)同一露头不同类型的暗色微粒包体,包体中见来自寄主岩的暗色环边石英板岩和长柱状角闪石斑晶;(d)包体中见来自寄主岩的斜长石斑晶;(e)黑云闪长岩包体的细粒半自形粒状结构(+);(f)石英闪长岩包体的细粒半自形粒状结构(+);(g)黑云石英二长闪长岩包体的中细粒半自形粒状结构(+);(h)寄主黑云母二长花岗岩的中粒似斑状结构(+);(i)寄主黑云母钾长花岗岩中粒似斑状结构(+);(j)寄主岩黑云母花岗闪长岩的中粒花岗结构(+);(k)寄主石英闪长岩的细粒半自形粒状结构,见针状磷灰石(-);(l)包体中的长柱状斜长石斑晶(+).矿物代号:Kf-钾长石;Pl-斜长石;Am-角闪石;Bi-黑云母;Qz-石英;Ap-磷灰石;Mt-磁铁矿;Spn-榍石.薄片样品采样位置详见图 1,下同 Fig. 4 Petrological characteristics of mafic microgranular enclaves and their host rocks

寄主岩与暗色微粒包体的矿物组成一致,差异主要体现在包体中暗色矿物(角闪石)的含量明显比寄主岩高,而寄主岩中石英含量明显高于包体,此外包体粒度明显小于寄主岩,且晶型较差,多呈半自形-他形,在二者接触处上述差异十分明显。黑云母二长花岗岩呈灰白色,具中粒似斑状结构-中粒花岗结构,矿物组成为斜长石(35%),钾长石(24%),石英(33%),黑云母(7%),角闪石(1%),副矿物见榍石和磷灰石(图 4h);黑云母钾长花岗岩呈肉红色-浅肉红色,具中粒似斑状结构-中粒花岗结构,矿物组成为钾长石(48%),斜长石(24%),石英(20%),黑云母(6%),角闪石(0.5%),榍石(1%)以及少量磷灰石(图 4i);黑云母花岗闪长岩呈灰白色,具中粒结构,矿物组成为斜长石(63%)、石英(20%)、钾长石(5%)、黑云母(6%)、角闪石(2%),副矿物见少量磁铁矿(图 4j);石英闪长岩呈浅灰色,具细粒半自形粒状结构,矿物组成为斜长石(66%)、角闪石(26%)、黑云母(8%)、石英(2%),副矿物见榍石及少量磁铁矿(0.5%)(图 4k)。

浅灰色中粒暗色微粒包体广泛发育斑晶,主要见钾长石、斜长石、角闪石和石英斑晶(图 4b-d),而深灰色中细粒暗色微粒包体和黑色细粒暗色微粒包体所含斑晶矿物明显减少(图 4l),或者不含斑晶。与寄主岩矿物组成对比发现,暗色微粒包体中的斑晶矿物主要来源于寄主岩,表明两种岩浆的温差不大,彼此之间存在成分交换。

2.3 包体与寄主岩接触关系

根据野外观察发现,暗色微粒包体与寄主岩的接触关系主要为截然型和渐变过渡型。一部分包体与寄主岩的接触关系为截然,该类型接触关系的最大特征是接触带暗色微粒包体的色率、粒度与寄主岩差异明显(图 5a),有时在暗色微粒包体一侧还能见到冷凝边(图 5b);另一部分暗色微粒包体与寄主岩呈渐变过渡,该类型接触关系的特征是暗色微粒包体和寄主岩之间色率和粒度差异较小,在暗色微粒包体中能见到寄主岩斑晶,偶见斑晶矿物横跨寄主岩和暗色微粒包体(图 5c),两者之间偶见明显的过渡带(图 5d);此外,还发现少量的反向脉(图 5e),在暗色微粒包体中见到黑云母钾长花岗岩细脉,脉体呈网脉状,与包体呈截然或渐变过渡关系,脉体色率明显低于包体,但粒度比包体大。其中截然型表明两种端元岩浆的温差较大,由高温的镁铁质岩浆注入温度相对偏低的长英质岩浆中,快速冷凝形成,多分布在岩浆混合边界及附近,而渐变型表明两种端元岩浆之间温差较小,相互之间进行了充分的机械混合和化学混合(康磊等, 2009; 顾枫华等, 2015)。

图 5 暗色微粒包体与寄主岩的接触关系 (a)包体与寄主岩之间为截然型接触关系;(b)包体与寄主岩接触带发育冷凝边;(c)角闪石斑晶横跨寄主岩和包体;(d)包体与寄主岩之后存在一条不规则的过渡带;(e)包体中发育的方向脉;(f)包体与寄主岩之间相互包裹;(g)包体包裹寄主岩;(h)包体之间相互包裹;(i)同一露头出露两种类型的包体 Fig. 5 Contact relationship of mafic microgranular enclaves and their host rocks

主要见寄主岩包裹暗色微粒包体,但偶尔可见包体与寄主岩相互包裹、穿插(图 5f),或包体包裹寄主岩(图 5g),以及包体之间相互包裹(图 5h)。在同一个露头上可见到两种矿物组成(或色率)、结构完全不同的暗色微粒包体(图 5i),偶见寄主岩(黑云母钾长花岗岩)和暗色微粒包体共同被后期的细粒钾长花岗岩脉穿切。

3 分析方法

所有样品均采自阿翁错复式岩体,选择新鲜具有代表性的样品,室内处理样品首先用切割机将样品表面部分切割掉,然后将样品粉碎至200目。主量、微量和稀土元素分析在澳实矿物实验室(广州)完成。主量元素的分析,首先将试样加入包含硝酸锂在内的助熔剂,充分混合后,高温熔融,熔融物倒入铂金模子形成扁平玻璃片后,再用产自荷兰的Philips PW2404型X射线荧光光谱仪(XRF)分析完成,主量元素分析所用标样为NCSDC47009、SARM-2、SARM-3、SARM-4,其中FeO分析所用标样为SY-4和YSBC-11701;微量元素和稀土元素的分析,首先将试样加入到偏硼酸锂/四硼酸锂熔剂中,混合均匀,在1025℃以上的熔炉中熔化,待熔液冷却后,用硝酸、盐酸和氢氟酸定容,再用产自美国的Perkin Elmer Elan 9000型电感耦合等离子体质谱仪(ICP-MS)分析完成,微量元素和稀土元素分析所用标样为OGGeo08和OREAS-45e。主量元素分析误差优于5%,稀土、微量元素的分析误差优于10%。

锆石样品的挑选由河北省区域地质矿产调查研究所实验室完成,岩石样品经人工破碎至80~100目,按照常规方法分选出锆石,然后在双目镜下挑选透明,晶型完好,无裂痕和包裹体且具代表性的锆石颗粒粘贴于环氧树脂的表面,固化后打磨抛光至露出一个光洁表面(宋彪等, 2002)。锆石的透、反射照相,阴极发光(CL)成像以及LA-ICP-MS锆石U-Pb定年测试分析在武汉上谱分析科技有限责任公司实验室完成,锆石定年分析所用仪器为安捷伦电感耦合等离子体质谱仪(Agilent 7700)及与之配套的相干193nm准分子激光剥蚀系统(GeoLasPro),激光能量80mJ,频率5Hz,激光束斑直径32μm。锆石U-Pb同位素比值校正采用91500标准锆石作为外标,微量元素含量测定采用NIST610作为外标,29Si作为内标,同位素比值监控标准样品为GJ-1,此次分析的分析流程为:NIST610→2个91500→2个GJ-1→5个样品(U-Pb年龄)→2个91500→……→2个GJ-1→2个91500→NIST610。具体分析条件及流程详见文献Liu el al. (2008, 2010)。采用ICPMSDataCal10.7程序处理数据。

4 分析结果 4.1 暗色包体与寄主岩的锆石U-Pb年龄

本次研究工作采集了寄主岩(黑云母钾长花岗岩),以及出露在黑云母钾长花岗岩中的暗色微粒包体进行锆石U-Pb定年。

寄主岩的锆石呈自形长柱状,长径120~200μm,长宽比约为1.5:1~3:1,锆石具有清晰的内部结构和典型的岩浆成因震荡环带(吴元保和郑永飞, 2004),晶型完好,裂纹不发育(图 6a)。暗色微粒包体中的锆石多呈短柱状,它形柱状,长径100~200μm,长宽比约1:1~2:1,部分锆石具有明显的震荡环带,晶型较好,整体上裂纹不发育(图 6b)。同位素分析结果(表 2)显示寄主岩和暗色微粒包体锆石的Th,U含量分别为233×10-6~871×10-6,263×10-6~694×10-6和135×10-6~2494×10-6,134×10-6~1385×10-6,二者呈现出良好的正相关,Th/U比值分别为0.89~1.29(均值1.02)和0.47~2.42(均值1.34),具有岩浆锆石的特征(Belousova et al., 2002)。

图 6 寄主岩(a)与暗色包体(b)锆石阴极发光图像 Fig. 6 CL images of zircon of the MMEs (b) and their host rocks (a)

表 2 阿翁错复式岩体及其中暗色微粒包体的LA-ICP-MS锆石U-Pb分析结果 Table 2 LA-ICP-MS U-Pb isotope data for zircons from the Awengcuo composite plutons and entrained MMEs

本次采集的寄主岩样品PM18N1共测定了20颗锆石20个测点,其中18个测点的206Pb/238U表面年龄为106.2~113.5Ma,且均位于谐和线上(图 7a),表明这些锆石形成后其U-Pb同位素体系保持封闭状态。PM18N1-7测点的年龄为99.3Ma,明显偏小,可能是所测定的锆石发生了少量放射成因铅丢失的结果,PM18N1-16测点的207Pb/235U和206Pb/238U年龄不谐和,为了数据的质量,剔除这两个测点的年龄,剩下18个测点的206Pb/238U年龄加权平均值为109.1±1.0Ma(MSWD=1.4,n=18)(图 7a)。暗色微粒包体样品PM40N1共测定了25颗锆石25个测点,PM40N1-2、PM40N1-20两个测点的207Pb/235U和206Pb/238U年龄明显不谐和(图 7b),根据表 2图 7b,暗色微粒包体样品的25个测点中,除了PM40N1-2、PM40N1-20两个测点明显不谐和外,其余23个测点的206Pb/238U表面年龄可以分为两组,其中一组14个测点的206Pb/238U表面年龄为104.9~110.6Ma,加权平均值为107.4±0.7Ma(MSWD=1.1,n=14);另一组9个测点的206Pb/238U表面年龄为98.4~101.8Ma,明显偏小(图 7b浅色线表示),可能是所测定的锆石发生了少量放射成因铅丢失的结果,为了数据的质量,剔除这11个测点的年龄,以107.4±0.7Ma作为暗色微粒包体样品PM40N1的成岩年龄。上述测年结果表明阿翁错复式岩体寄主岩和包体的成岩年龄分别为109.1±1.0Ma和107.4±0.7Ma,在误差范围内基本上一致。

图 7 寄主岩(a)与暗色包体(b)锆石U-Pb谐和图 Fig. 7 U-Pb concordia diagram of age data for zircons from MMEs (b) and their host rocks (a)
4.2 暗色包体与寄主岩的岩石化学特征

暗色包体与寄主岩的主微量元素分析结果见表 3。寄主岩的SiO2含量为67.10%~70.02%(平均68.00%),为酸性岩类,暗色微粒包体SiO2含量为52.08%~57.07%(平均54.07%),为中性岩类。寄主岩的全碱含量为8.05%~8.39%,平均8.25%,包体的全碱含量为3.95%~6.92%,平均5.05%,其中寄主岩的K2O/Na2O>1,而包体的K2O/Na2O < 1,表明寄主岩比包体更富钾。在K2O-SiO2关系图中,寄主岩落在高钾钙碱性系列区域,而暗色包体落在钙碱性系列和高钾钙碱性系列区域(图 8a)。寄主岩的A/CNK=0.95~1.00,平均0.98,而包体的A/CNK=0.76~0.95,平均0.85,整体为准铝质(图 8b)。寄主岩的分异指数(DI)=79.02~84.01,平均80.44,包体的分异指数=36.84~54.30,平均45.04,表明寄主岩分异程度明显高于包体。

表 3 阿翁错复式岩体及暗色微粒包体全岩地球化学数据(主量元素:wt%; 稀土和微量元素:×10-6) Table 3 Geochemical data of Awengcuo composite plutons and their MMEs (major elements: wt%; trace elements: ×10-6)

图 8 阿翁错复式岩体的SiO2与K2O关系图解(a, 据Rickwood, 1989)和A/CNK-A/NK图解(b, 据肖庆辉等, 2002) Fig. 8 Diagrams of SiO2 vs. K2O (a, after Rickwood, 1989) and A/CNK vs. A/NK (b, after Xiao et al., 2002) for Awengcuo composite plutons and their MMEs

寄主岩和暗色微粒包体的稀土元素含量分别为166.5×10-6~231.6×10-6(平均213.3×10-6)和76.57×10-6~228.9×10-6(平均148.0×10-6),前者高于后者(表 2);LREE/HREE分别为9.96~16.8(平均12.2)和3.84~9.32(平均6.30),(La/Yb)N分别为9.35~16.8(平均12.1)和2.68~10.5(平均6.08),(Gd/Yb)N分别为1.16~1.54(平均1.32)和1.29~1.92(平均1.54),表明二者均为LREE强烈富集,分馏明显,但寄主岩的分异程度明显高于包体,而HREE相对亏损,分布较缓,分馏较明显;δEu分别为0.60~0.75(平均0.69)和0.70~0.94(平均0.84),寄主岩负铕异常中等,而暗色包体负铕异常相对较弱,表明岩浆演化过程中可能存在斜长石的分离结晶,从包体(中性岩)到寄主岩(酸性岩)负铕异常逐渐增强。稀土元素配分曲线图(图 9a, b)呈向右陡倾斜下凹浅的“Ⅴ”字型,轻、重稀土分馏明显,二者配分样式基本上一致。

图 9 暗色包体与寄主岩球粒陨石标准化稀土元素配分曲线(a、b, 标准化值据Boynton, 1984)和原始地幔标准化微量元素蛛网图(c、d, 标准化值据McDonough et al., 1992) Fig. 9 Chondrite-normalized REE patterns (a, b, normalization values after Boynton, 1984) and primitive mantle-normalized trace elements spidergram (c, d, normalization values after McDonough et al., 1992) of MMEs and their host rocks

寄主岩和暗色微粒包体具有一致的痕量元素特征和标准化蛛网图样式,表现为向右缓倾的“多峰多谷”形态(图 9c, d),二者均富集大离子亲石元素Rb、K和不相容元素U、Th,相对亏损高场强元素Nb、Ta、P、Ti和大离子亲石元素Ba、Sr等。

5 讨论 5.1 岩体的岩浆混合成因

已有研究认为暗色微粒包体具有以下几种成因类型:(1)寄主岩母岩浆早期结晶的凝块、基性岩墙经破碎而成(同源捕虏体)(张德全, 1981);(2)寄主岩源区熔融后的残留体(Chappell and White, 1991);(3)寄主岩侵位过程中捕获的围岩捕虏体(异源捕掳体,根据来源不同分为就地和深部两种)(Maas et al., 1997),是同化混染作用不彻底的产物;(4)寄主岩早期的堆晶矿物,即析离体或异离体(Didier, 1973);(5)镁铁质岩浆与长英质岩浆混合作用不彻底的产物(MME)(Perugini et al., 2003; Feeley et al., 2008; 陈广俊等, 2014)。

阿翁错复式岩体的围岩是上侏罗-下白垩统沙木罗组下段,主要为一套细砂质板岩、粉砂质板岩、泥质板岩、变中-细粒钙质岩屑石英砂岩、变钙质石英粉砂岩、变细粒钙质凝灰质长石杂砂岩组合,在成分和结构上与暗色微粒包体存在明显差异,野外调查在阿翁错复式岩体中发现了沙木罗组下段地层中的板岩捕掳体,捕掳体呈不规则棱角状,见角岩化、硅化等蚀变,具变余层理构造,而暗色微粒包体未见角岩化、硅化等蚀变现象,具典型的火成结构;此外,暗色微粒包体塑性变形特征明显,呈椭圆形、细长条形、透镜状以及不规则形态,与围岩捕掳体的不规则棱角状明显不同,表明暗色微粒包体不是围岩捕掳体。镜下观察发现暗色微粒包体具有典型的岩浆岩矿物组合和结构特征,未见石榴石、红柱石、堇青石等特征蚀变矿物,表明暗色微粒包体不是源区熔融后的残留体(陈广俊等, 2014)。暗色微粒包体为寄主岩早期的堆晶矿物的最大特征是暗色微粒包体具堆晶结构,矿物粒度与寄主岩相近,接触关系多为渐变(Noyes et al., 1983)。阿翁错复式岩体中的暗色微粒包体未发现堆晶结构,大部分包体的粒度明显小于寄主岩粒度,接触关系见截然型和渐变型,表明暗色微粒包体不是寄主岩早期的堆晶矿物(肖庆辉等, 2002; 陈国超等, 2016)。少部分暗色微粒包体与寄主岩接触带见宽约几毫米的冷凝边(图 5b),偶见反向脉(图 5e),而冷凝边和反向脉的出现是镁铁质岩浆与长英质岩浆发生混合作用的可靠证据(王德滋和谢磊, 2008)。

以上分析表明镁铁质岩浆注入长英质岩浆发生混合作用的观点可以合理解释阿翁错复式岩体中暗色微粒包体的成因。

阿翁错复式岩体发生岩浆混合作用具有诸多识别标志。在复式岩体中广泛分布闪长质暗色微粒包体,包体直径悬殊明显,以小型包体为主,具明显的塑性变形特征,与寄主岩接触关系以截然为主,见少量渐变过渡型;包体与寄主岩之间偶见宽约几厘米到几十厘米不等的过渡带,过渡带沿着包体与寄主岩的接触带呈环状或脉状产出,该带以色率介于包体和寄主岩石之间,粒度大于寄主岩而小于包体为显著特征(图 5d);有时可见包体与寄主岩之间相互穿插、包裹(图 5f, g);在包体与寄主岩接触带包体一侧偶见宽约1~2mm的冷凝边(图 5b),包体中偶见反向脉发育(图 5e);在包体中可见到大量从寄主岩捕获的斑晶(Troll and Schmincke, 2002; 陈斌等, 2006; 陈国超等, 2016),如发育暗色镶边的石英和钾长石斑晶(图 4b, c),镶边是由于基性岩浆的温度高于酸性岩浆约300℃(王德滋和谢磊, 2008),当两者之间发生混合作用时,吸热效应在石英、钾长石边缘的一圈熔体中形成局部过冷的条件,导致细粒他形角闪石、黑云母等暗色矿物围绕石英、钾长石捕掳晶晶出(Vernon, 2014; 陈广俊等, 2014; 陈国超等, 2016),此外,在包体中还见到斜长石斑晶,这些斑晶以发育溶蚀状和港湾状外形为特征(图 4d);在包体及包体周围寄主岩中见大量不平衡结构:包体中的长柱状斜长石(图 10a)、角闪石(图 10b, c)和针状磷灰石等结构特征(图 4k)。

图 10 暗色微粒包体中的不平衡结构 (a)包体中见来自寄主岩的长柱状斜长石斑晶;(b)包体中的长柱状角闪石;(c)镜下包体中的长柱状角闪石(-) Fig. 10 Unbalanced structure in mafic microgranular enclaves
5.2 岩浆混合的程度

岩浆混合的能干性主要受岩浆的温度和黏性影响(张旗等, 2007; 王玉往等, 2012),基岩岩浆的温度恒高于酸性岩浆(王德滋和谢磊, 2008),且基岩岩浆的黏性小于酸性岩浆(张旗等, 2007),故基性岩浆混合的能干性远远强于酸性岩浆,在酸性岩石中常见基性暗色微粒包体,而在基岩岩石中很少见酸性包体。阿翁错复式岩体亦如此。在暗色微粒包体与寄主岩的接触带见冷凝边(图 5b),且从接触带到中心矿物粒度逐渐变大(图 4a),表明镁铁质岩浆注入长英质岩浆的早期阶段由于两者之间具有较大温差,包体边部温度快速降低,矿物快速结晶形成颗粒细小的冷凝边,由于未经历充分的结晶过程,矿物晶形差,粒度偏小,多呈他形-半自形细粒结构;晚期阶段随着暗色微粒包体与寄主岩之间温差逐渐降低,包体内部残余的岩浆经历了相对充分的结晶过程,形成的矿物晶型较好,粒度偏大。此外包体中见长柱状斜长石(图 6a)、角闪石(图 6b, c)和针状磷灰石(图 4h)也是温度较高的镁铁质岩浆注入温度相对较低的长英质岩浆中快速冷却的标志(Baxter and Feely, 2002)。部分暗色微粒包体与寄主岩之间呈渐变过渡,偶见两者之间存在明显的过渡带(图 5d),表明两者之间温差较小。

根据岩浆混合的程度可将混合作用划分为低级和高级两个阶段,低级阶段以机械混合(mixling)为主,而高级阶段主要为化学混合(mingling),两个阶段是连续的,宏观和微观下见到的各种不协调、不均一现象主要出现在低级混合阶段,而高级阶段混合作用进行得彻底,岩石在组构上具均一性,见不到端元组分或混合残余物(李昌年, 1999; 张旗等, 2007; 王玉往等, 2012)。阿翁错复式岩体中广泛分布暗色微粒包体,能见到各种尺度的不协调现象,如包体与寄主岩之间相互包裹、包体中见反向脉、包体矿物晶形较差、斑晶具港湾状结构(石英斑晶发育暗色镶边)、针状磷灰石以及包体中矿物异常共生等,这些现象表明阿翁错复式岩体是岩浆混合作用不彻底的产物。

阿翁错复式岩体各类型岩石中均见暗色微粒包体分布,但分布极不均匀,从石英闪长岩-花岗闪长岩-花岗岩,暗色微粒包体数量呈增加趋势,尤其是在基性岩墙与花岗岩的接触带及附近暗色微粒包体成群、成带分布。这一现象表明从中性岩石到酸性岩石,岩浆混合程度在逐渐降低,这种趋势与中性岩浆的黏度低于酸性岩浆,而温度高于酸性岩浆有关。

5.3 岩浆起源及构造背景

阿翁错复式岩体的岩浆混合成因意味着暗色微粒包体和寄主岩相应端元之间存在成分的交换,但是包体明显比寄主岩偏中性,且包体与寄主岩之间的接触关系以截然为主,表明岩浆混合作用不彻底,存在物理混合(mingling),那么寄主岩和包体继承了相应端元的部分岩浆源区特征,为本文通过地球化学数据讨论阿翁错复式岩体的岩浆源区以及构造环境、地球动力学提供了大量信息。在以SiO2含量为横坐标的Harker图解中(图 11),寄主岩和暗色微粒包体的FeOT、Al2O3、MgO、CaO、P2O5、TiO2含量与之呈较好的负相关性,表明在岩浆演化过程中二者经历了辉石、斜长石、磷灰石及钛铁氧化物的分离结晶作用(李献华等, 2000; 付强等, 2011),二者的稀土元素配分曲线图和微量元素原始地幔标准化蛛网图具有高度一致性,表明二者具有强烈的地球化学亲源关系,且经历了相似的岩浆演化过程(陈广俊等, 2014)。

图 11 暗色微粒包体与寄主岩的Harker图解 Fig. 11 Harker diagram of MMEs and their host rocks

暗色微粒包体低硅(52.08%~57.07%),富Mg、Fe、Ca,高Mg#(45~58,平均51)的地球化学特征,暗示包体的原始岩浆起源于地幔的部分熔融(Zhou et al., 2005; 于介江等, 2013),但是包体富集大离子亲石元素Rb、K,亏损高场强元素Nb、Ta、P、Ti,暗示岩浆演化过程中受到地壳物质的混染或岩浆源区经历了俯冲带流体的交代,低硅又说明地壳物质混染不显著。包体的Nb/Ta比值为12.2~25.8(平均17.6)接近原始地幔值(17.5±2.0, Jochum et al., 1989)。由此认为,包体的原始岩浆可能来源于受俯冲带流体交代的地幔楔(许继峰等, 2001; 郝百武, 2012; 陈广俊等, 2014)。而寄主岩相对高硅(67.10%~70.02%),贫Mg、Fe、Ca,具有较低的Mg#(36~38,平均37),同样富集大离子亲石元素Rb、K,亏损高场强元素Nb、Ta、P、Ti,其Nb/Ta比值为7.17~11.6(平均9.85)介于下地壳(8.3, Rudnick and Gao, 2003)和原始地幔值(17.5±2.0, Jochum et al., 1989)之间,更倾向于下地壳,以上特征表明寄主岩为幔源岩浆与熔融下地壳岩浆混合的产物,以熔融下地壳岩浆为主。

寄主岩和暗色微粒包体的CaO/Na2O比值分别为0.59~0.74(平均0.67)和1.34~2.96(平均2.10),已有研究表明该比值介于0.3~1.5之间的中酸性岩石由杂砂岩或火成岩熔融形成,而由角闪岩熔融形成的中酸性岩石该比值变化较大,以0.1~3为主,可高达10(Jung and Pfänder, 2007)。此外二者的A/CNK值分别为0.95~1.01(平均0.98)和0.76~0.95(平均0.85),均小于1.1,根据已有研究认为二者不可能有硬砂岩和泥质岩石部分熔融形成(Rushmer, 1991; Montel and Vielzeuf, 1997; 马乐天等, 2010),但是角闪岩在水不饱和的条件下,随着水压的降低,脱水熔融形成岩浆的A/CNK值小于1(Rushmer, 1991; 马乐天等, 2010),结合二者均富Na、Ca,具有较低的K2O/Na2O比值(0.63~1.43,平均0.96),推断二者的物质源区可能以角闪岩为主(赵永久等, 2007),但是寄主岩物质源区有向杂砂岩过渡的趋势(图 12)。

图 12 暗色微粒包体与寄主岩的源区组成判别图(据Altherr and Siebel, 2002) Fig. 12 Discrimination diagrams of magma source for MMEs and their host rocks (after Altherr and Siebel, 2002)

目前对于班-怒洋的认识尚存诸多争议。不同地区的研究表明班-怒洋盘的打开具有穿时性,在印度板块地区打开的时间为中-晚二叠世(Stampfli and Borel, 2002),Spring et al. (1993)在冈瓦纳大陆北缘采集的花岗岩年龄为281±1Ma、270Ma和268±5Ma,该地区的岩浆活动集中在早二叠世晚期-晚二叠世(Zhu et al., 2010),黄启帅等(2012)Shi et al. (2012)分别对班公湖MOR型蛇绿岩和东巧蛇绿岩的地幔橄榄岩进行测年,获得了254±28Ma和251±65Ma两个年龄,可以代表班-怒洋的初始裂解时间,该时间比曲晓明等(2009)在日土县东面由洋脊型镁铁质岩墙测得的洋盆打开时间186.6~181.9Ma要早很多,林彬等(2017)对在多龙矿集区尕尔勤矿区发现的枕状玄武岩进行测年,获得其年龄为298±11Ma,将班-怒洋的初始裂解时限提前至早二叠世,Zhang et al. (2016)获得改则地区的榴辉岩年龄为260±5Ma,强巴扎西等(2009)对藏东丁青地区蛇绿岩中堆晶辉长岩进行测年,获得其年龄217±1.6Ma,这些年龄数据进一步证实班-怒洋盘的打开具有穿时性。班-怒洋盆存在三条俯冲带(杜德道等, 2011; 曲晓明等, 2009):北面的日土和狮泉河-改则-洞错俯冲带,南面的拉果错俯冲带,前者开始向北向羌塘地块俯冲的时间是165.5±1.9Ma~166.4±2.0Ma(曲晓明等, 2009),但羌塘南缘广泛发育早侏罗世(190~180Ma)SSZ型蛇绿岩指示北向俯冲作用此时已开始(Wang et al., 2016),后者开始向南向拉萨地块俯冲的时间为134Ma(杜德道等, 2011),到侏罗纪末-白垩纪初洋盆闭合(曲晓明等, 2009),随后进入陆内俯冲,该过程一直持续到约95Ma(Kapp et al., 2005),然而李华亮等(2016)根据上白垩统竞柱山组(K2j)为陆相磨拉石沉积,且角度不整合覆盖在蛇绿岩及特提斯残留海相地层之上,认为班公湖地区在96Ma左右全面完成了由洋到陆的转换,进入陆内环境,约88.8Ma完成了从陆内俯冲向局部伸展的构造环境转变(雷鸣等, 2015),李发桥等(2018)对班-怒带中段玛日埃错地区花岗斑岩进行研究,认为其形成于后碰撞伸展环境,表明晚白垩世晚期(78.3Ma)班怒带中段地区已经进入后碰撞伸展阶段。

阿翁错复式岩体中寄主岩和暗色微粒包体的Na2O含量平均值为2.80%~4.83%(平均均3.62%),小于4.5%,表明区内花岗岩类形成深度接近50km(马乐天等, 2010; Winther, 1996),显示此时的地壳已经加厚,证实了狮泉河-改则-洞错俯冲带向北向羌塘地块之下俯冲的观点(杜德道等, 2011)。岩体的成岩年龄是109.1±1.0Ma,此时正处于班-怒洋由弧-陆碰撞向陆陆碰撞的转换阶段即软碰撞阶段(任纪舜等, 1999)。寄主岩和暗色微粒包体均富集不相容元素Rb、K、U、Th,相对亏损高场强元素Nb、Ta、P、Ti,具有岛弧型岩浆岩的特征(杜德道等, 2011),在Pearce (1996)的构造环境判别图上所有样品主要落在火山弧构造环境(图 13),有个别样品落在板内构造环境,表明其形成与狮泉河-改则-洞错俯冲带向北向羌塘地块之下俯冲关系密切。

图 13 暗色微粒包体与寄主岩的Rb-Nb+Y图(a)和Rb-Yb+Ta图(b)(据Pearce, 1996) Fig. 13 Rb vs. Nb+Y diagram (a) and Rb vs. Yb+Ta diagram (b) for MMEs and their host rocks (after Pearce, 1996)

关于岩浆混合成因的机制有以下几种观点:(1)层状岩浆房对流观点。岩浆房中存在分层明显的两层岩浆,上部为长英质岩浆,下部为镁铁质岩浆,由于密度、黏度、温度等的不均衡性,两种岩浆发生大规模的对流作用,下部镁铁质岩浆的对流作用将长英质岩浆拖入,镁铁质岩浆混合能干性强,二者之间发生化学混合,形成均一岩浆;上部长英质岩浆的对流作用将镁铁质岩浆带入,而长英质岩浆混合能干性相对较差,二者之间仅发生机械混合作用,形成各种形态的暗色微粒包体(Vernon, 1983; 马昌前等, 1994; 王德滋等, 1992; 张旗等, 2007; 王德滋和谢磊, 2008)。(2)暗色微粒包体来源于富气的镁铁质岩浆上浮到长英质岩浆中结晶形成(Eichelberger, 1980; 王德滋等, 1992; 王德滋和谢磊, 2008)。(3)镁铁质岩浆以喷泉或射流的形式从岩浆房的一侧喷入到长英质岩浆(Frost and Mahood, 1987)。(4)深部岩浆房的镁铁质岩浆经伸展构造注入到浅部长英质岩浆房,发生岩浆混合作用,形成暗色微粒包体或基性岩墙(王德滋等, 1992; 王德滋和谢磊, 2008; Kumar and Pieru, 2010)。(5)岩浆混合作用主要发生在深部源区,在岩浆上升或侵位阶段发生混合作用的可能性很小(张旗等, 2007)

根据野外观察统计,按照拉伸程度将研究区暗色微粒包体分为两类。第一类包体呈细长条状,定向明显,长短轴比值介于4.14~5.32之间(表 1),包体分布的位置距离中基岩岩墙与寄主岩的接触部位较远,包体和寄主岩均未见固态变形,其拉伸程度与塑性状态下岩浆的流动有关,其长轴方向指示岩浆的流动方向(马昌前等, 1994; 陈国超等, 2016);第二类包体呈椭圆状、透镜状、云雾状以及不规则形态,直径明显大于第一类包体,未见明显拉伸,多成群、成带分布在中基性岩墙与寄主岩接触部位附近,包体与寄主岩结构差异明显,包体常发育冷凝边,这类包体未见明显定向,或在同一露头可见几种不同的方向,为岩浆局部对流、搅动的结果(Paterson et al., 2004; 陈国超等, 2016)。在岩浆混合边界,端元岩浆的温度、黏度差异明显(Kumar and Rino, 2006),岩浆运移速度较小,包体拉伸程度较低,形成第二类包体,随着远离混合边界,岩浆运移速度增大,包体拉伸程度也随之变大,形成第一类包体(陈国超等, 2016)。

岩浆混合作用的发生与构造环境关系密切,俯冲造山带、大陆热点、地壳伸展、深大断裂等构造环境有利于岩浆混合作用的发生(Shamberger and Garcia, 2007; 齐有强等, 2008)。阿翁错复式岩体中寄主岩和暗色微粒包体的成岩年龄分别为109.1±1.0Ma和107.4±0.7Ma,此时正处于班-怒洋北向俯冲羌塘地块之下的弧-陆碰撞造山阶段(Kapp et al., 2005),当洋壳俯冲到一定深度,随着温度升高、压力增大,洋壳开始脱水,产生流体,引起上覆地幔楔发生部分熔融,形成镁铁质岩浆(金性春和于开平, 2003; 董国臣等, 2006; 邱检生等, 2015),在密度差的作用下镁铁质岩浆向上运移,并底侵于壳-幔边界,为下地壳带来高温热能(温差300~500℃)和富含挥发分的流体,引发下地壳物质发生部分熔融,形成长英质岩浆(Huppert and Sparks, 1988; 董国臣等, 2006; 王德滋和谢磊, 2008),当镁铁质岩浆从底部注入到已经部分结晶的长英质岩浆中时,由于二者之间温度、黏度差异明显,不能发生彻底的化学混合(mixing)形成均一的岩浆,而以发生机械混合(mingling)为主(张旗等, 2007; 邱检生等, 2015),表现为各种岩浆混合现象以及不协调现象(张旗等, 2007)。

6 结论

(1) 岩浆混合作用以机械混合为主,存在化学混合,且岩浆混合作用具有多期次性,岩浆混合作用的影响因素主要是端元岩浆的温度和黏度,其中温差导致暗色微粒包体和寄主岩之间存在不同的组构。

(2) 暗色微粒包体的拉伸程度与塑性状态下岩浆的流动有关,具有明显定向性包体的长轴方向指示岩浆的流动方向,而包体未见明显定向,或在同一露头可见几种不同的方向,为岩浆局部对流、搅动的结果。

(3) 阿翁错复式岩体中寄主岩和暗色微粒包体的年龄分别为109.1±1.0Ma和107.4±0.7Ma,基本上一致,岩浆混合作用发生在早白垩世,此时正处于班-怒洋由弧-陆碰撞向陆陆碰撞的转换阶段即软碰撞阶段。

(4) 在班-怒洋北向俯冲羌塘地块之下的背景下,洋壳脱水,引起上覆地幔楔发生部分熔融,形成镁铁质岩浆,镁铁质岩浆向上运移,并底侵于壳-幔边界,引发下地壳物质发生部分熔融,形成长英质岩浆,当镁铁质岩浆从底部注入长英质岩浆房时,镁铁质岩浆快速冷凝,形成部分色率高、粒度细,具冷凝边的包体,与寄主岩呈截然型接触,随着端元岩浆之间的温差逐渐降低,包体色率降低,粒度变大,与寄主岩呈渐变过渡。

致谢      四川省地质调查院阚泽忠教授级高工在成文过程中给予了帮助;两位匿名审稿专家在审稿过程中提出了宝贵的意见和建议,使文章得到了进一步提升和完善;宋扬副研究员、林彬博士及编辑部老师为本文的顺利刊出付出了努力;在此一并表示衷心的感谢。

参考文献
Altherr R and Siebel W. 2002. Ⅰ-type plutonism in a continental back-arc setting:Miocene granitoids and monzonites from the central Aegean Sea, Greece. Contributions to Mineralogy and Petrology, 143(4): 397-415. DOI:10.1007/s00410-002-0352-y
Barbarin B. 2005. Mafic magmatic enclaves and mafic rocks associated with some granitoids of the central Sierra Nevada batholith, California:Nature, origin, and relations with the hosts. Lithos, 80(1-4): 155-177. DOI:10.1016/j.lithos.2004.05.010
Baxter S and Feely M. 2002. Magma mixing and mingling textures in granitoids:Examples from the Galway granite, Connemara, Ireland. Mineralogy and Petrology, 76(1-2): 63-74.
Belousova EA, Griffin WL, O'reilly SY and Fisher NI. 2002. Igneous zircon:Trace element composition as an indicator of source rock type. Contributions to Mineralogy and Petrology, 143: 602-622. DOI:10.1007/s00410-002-0364-7
Bonin B. 2004. Do coeval mafic and felsic magmas in post-collisional to within-plate regimes necessarily imply two contrasting, mantle and crustal, sources? A review. Lithos, 78(1-2): 1-24. DOI:10.1016/j.lithos.2004.04.042
Bonin B. 2007. A-type granites and related rocks:Evolution of a concept, problems and prospects. Lithos, 97(1-2): 1-29. DOI:10.1016/j.lithos.2006.12.007
Boynton WV. 1984. Geochemistry of the rare earth elements: Meteorite studies. In: Henderson P (ed.). Rare Earth Element Geochemistry. Amsterdam: Elsevier, 63-114
Castro A. 2013. Tonalite-granodiorite suites as cotectic systems:A review of experimental studies with applications to granitoid petrogenesis. Earth-Science Reviews, 124: 68-95. DOI:10.1016/j.earscirev.2013.05.006
Champion DC and Bultitude RJ. 2013. The geochemical and Sr-Nd isotopic characteristics of Paleozoic fractionated S-types granites of North Queensland:Implications for S-type granite petrogenesis. Lithos, 162-163: 37-56. DOI:10.1016/j.lithos.2012.11.022
Chappell BW and White AJR. 1991. Restite enclaves and the restite model. In: Didier J and Barbarin B (eds.). Enclaves and Granite Petrology. Amsterdam: Elsevier, 375-381
Chen B, Liu CQ and Tian W. 2006. Magma-mixing between mantle-and crustal-derived melts in the process of Mesozoic magmatism, Taihangshan:Constraints from petrology and geochemistry. Earth Science Frontiers, 13(2): 140-147.
Chen GJ, Sun FY, Li BL, Wang G, Qian Y, Xu QL and Liu GC. 2014. Geochemistry and petrogenesis of Gouli mafic enclaves and their host rocks in eastern Kunlun. Journal of Jilin University (Earth Science Edition), 44(3): 892-904.
Chen GC, Pei XZ, Li RB, Li ZC, Liu CJ, Chen YX, Xu T and Zhang YM. 2016. Genesis of magma mixing and mingling of Xiangjiananshan granite batholith in the eastern section of East Kunlun orogen:Evidence from mafic microgranular enclaves (MMEs). Earth Science Frontiers, 23(4): 226-240.
Didier J. 1973. Granites and Their Enclaves: The Bearing of Enclaves on the Origin of Granites. Amsterdam: Elsevier, 1-393
Ding S, Tang JX, Zheng WB, Yang C, Zhang Z, Wang Q and Wang YY. 2017. Geochronology and geochemistry of Naruo porphyry Cu (Au) deposit in Duolong ore-concentrated area, Tibet, and their geological significance. Earth Science, 42(1): 1-23.
Dong GC, Mo XX, Zhao ZD, Zhu DC, Wang LL, Chen T and Li B. 2006. Magma mixing in middle part of Gangdise magma belt:Evidences from granitoid complex. Acta Petrologica Sinica, 22(4): 835-844.
Du DD, Qu XM, Wang GH, Xin HB and Liu ZB. 2011. Bidirectional subduction of the Middle Tethys oceanic basin in the west segment of Bangonghu-Nujiang suture, Tibet:Evidence from zircon U-Pb LAICPMS dating and petrogeochemistry of arc granites. Acta Petrologica Sinica, 27(7): 1993-2002.
Eichelberger JC. 1980. Vesiculation of mafic magma during replenishment of silicic magma reservoirs. Nature, 288(5790): 446-450. DOI:10.1038/288446a0
Feeley TC, Wilson LF and Underwood SJ. 2008. Distribution and compositions of magmatic inclusions in the Mount Helen dome, Lassen Volcanic Center, California:Insights into magma chamber processes. Lithos, 106(1-2): 173-189. DOI:10.1016/j.lithos.2008.07.010
Frost TP and Mahood GA. 1987. Field, chemical, and physical constraints on mafic-felsic magma interaction in the Lamarck granodiorite, Sierra Nevada, California. Geological Society of America Bulletin, 99(2): 272-291. DOI:10.1130/0016-7606(1987)99<272:FCAPCO>2.0.CO;2
Fu Q, Ge WS, Wen CS, Cai KQ, Li SF, Zhang ZW and Li XF. 2011. Geochemistry and genesis of Michang granites and their dark microgranular enclaves in Guangxi. Acta Geoscientica Sinica, 32(3): 293-303.
Gu FH, Zhang YM, Liu RP, Zheng L and Sun X. 2015. Magma mixing and mingling of the Shadegai granite in Inner Mongolia:Evidence from petrography, mineral chemistry and geochronology. Acta Petrologica Sinica, 31(5): 1374-1390.
Hao BW. 2012. Discovery, genesis and tectonic significance of the Late-Paleozoic miarolitic K-feldspar granite in southern Xianghuangqi, Inner Mongolia. Journal of Jilin University (Earth Science Edition), 42(2): 269-284.
Huang QS, Shi RD, Ding BH, Liu DL, Zhang XR, Fan SQ and Zhi XC. 2012. Re-Os isotopic evidence of MOR-type ophiolite from the Bangong Co for the opening of Bangong-Nujiang Tethys Ocean. Acta Petrologica et Mineralogica, 31(4): 465-478.
Huppert EH and Sparks RSJ. 1988. The generation of granitic magmas by intrusion of basalt into continental crust. Journal of Petrology, 29(3): 599-624. DOI:10.1093/petrology/29.3.599
Jin XC and Yu KP. 2003. Subduction factory and subduction recycling of continental material. Advance in Earth Sciences, 18(5): 737-744.
Jochum KP, McDonough WF, Palme H and Spettel B. 1989. Compositional constraints on the continental lithospheric mantle from trace elements in spinel peridotite xenoliths. Nature, 340(6234): 548-550. DOI:10.1038/340548a0
Jung S and Pfänder JA. 2007. Source composition and melting temperatures of orogenic granitoids:Constraints from CaO/Na2O, Al2O3/TiO2 and accessory mineral saturation thermometry. European Journal of Mineralogy, 19(6): 859-870. DOI:10.1127/0935-1221/2007/0019-1774
Kang L, Li YJ, Zhang B, Zhang HW and Wang JN. 2009. Petrographic evidence for magma mixing of Xiaerpu granite in West Junggar, Xinjiang. Acta Petrologica et Mineralogica, 28(5): 423-432.
Kapp P, Yin A, Harrison TM and Lin D. 2005. Cretaceous-Tertiary shortening, basin development, and volcanism in central Tibet. Geological Society of America Bulletin, 117(7-8): 865-878.
Kumar S and Rino V. 2006. Mineralogy and geochemistry of microgranular enclaves in Palaeoproterozoic Malanjkhand granitoids, central India:Evidence of magma mixing, mingling, and chemical equilibration. Contributions to Mineralogy and Petrology, 152(5): 591-609. DOI:10.1007/s00410-006-0122-3
Kumar S and Pieru T. 2010. Petrography and major elements geochemistry of microgranular enclaves and Neoproterozoic granitoids of South Khasi, Meghalaya:Evidence of magma mixing and alkali diffusion. Journal of the Geological Society of India, 76(4): 345-360. DOI:10.1007/s12594-010-0106-9
Lei CY, Li ZJ, Zhang Z, Hu ZH, Wang HX and Song JL. 2012. Geochemical characteristics and geodynamic significance of the granites in the Ga'erqiong Cu-Au deposit, Tibet. Acta Geoscientica Sinica, 33(4): 601-612.
Lei CY, Wu JL, Yin XK, Liu W, Wang B, Li W, Yuan HY, Zhang W, Yin T and Pei YL. 2018. New discovery of the diorite porphyry dyke in the Shamuluo formation in western segment of the Bangongco-Nujiang suture zone and its geological significance. Bulletin of Mineralogy, Petrology and Geochemistry, 37(2): 250-259.
Lei CY, Wu JL, Yin XK, Liu W, Wang B, Li W, Yuan HY, Zhang W, Yin T and Pei YL. 2019. Petrogenesis of intermediate-acid intrusive rocks in Awengcuo area in the west segment of the Bangongco-Nujiang Suture, Tibet, and its instructions to the regional tectonic evolution. Geotectonica et Metallogenia, 43(1): 168-185.
Lei M, Chen JL, Xu JF and Zeng YC. 2015. Geochemistry of early Late Cretaceous Gaerqiong high-Mg# diorite porphyry in midnorthern Lhasa terrane:Partial melting of delaminated lower continental crust?. Geological Bulletin of China, 34(2-3): 337-346.
Li CN. 1999. Discovery for Precambrian magmatic mixed zone from northeastern Jiangxi Province and their geological significance. Earth Science Frontiers, 6(4): 331-337.
Li FQ, Liu ZB, Tang JX, Song Y, Gao K, Li HF, Zhao ZQ and Teng L. 2018. Petrogenesis of granite porphyry in mariaicuo area, Shuanghu County, Tibet, and constraints on the evolution in the middle section of Bangonghu-Nujiang suture zone. Earth Science, 43(4): 1051-1069.
Li GM, Zhang XN, Qin KZ, Sun XG, Zhao JX, Yin XB, Li JX and Yuan HS. 2015. The telescoped porphyry-high sulfidation epithermal Cu (-Au) mineralization of Rongna deposit in Duolong ore cluster at the southern margin of Qiangtang Terane, Central Tibet:Integrated evidence from geology, hydrothermal alteration and sulfide assemblages. Acta Petrologica Sinica, 31(8): 2307-2324.
Li HL, Gao C, Li ZH, Zhang Z, Peng ZM and Guan JL. 2016. Age and tectonic significance of Jingzhushan Formation in Bangong Lake area, Tibet. Geotectonica et Metallogenia, 40(4): 663-673.
Li JX, Qin KZ and Li GM. 2006. Basic characteristics of gold-rich porphyry copper deposits and their ore sources and evolving processes of high oxidation magma and ore-forming fluid. Acta Petrologica Sinica, 22(3): 678-688.
Li JX. 2008. Geochronology, petrology and metallogeneses of high oxidized magma-hydrothermal fluid of Duobuza gold-rich porphyry copper deposit in Bangonghu belt, Northern Tibet. Ph. D. Dissertation. Beijing: Institute of Geology and Geophysics, Chinese Academy of Sciences 1-224 (in Chinese with English summary)
Li JX, Qin KZ, Li GM, Xiao B, Zhao JX, Cao MJ and Chen L. 2013. Petrogenesis of ore-bearing porphyries from the Duolong porphyry Cu-Au deposit, central Tibet:Evidence from U-Pb geochronology, petrochemistry and Sr-Nd-Hf-O isotope characteristics. Lithos, 160-161: 216-227. DOI:10.1016/j.lithos.2012.12.015
Li XH, Zhou HW, Liu Y, Lee CY, Chen ZH, Yu JS and Gui XT. 2000. Mesozoic shoshonitic intrusives in the Yangchun Basin, western Guangdong, and their tectonic significance:Ⅰ. Petrology and isotope geochronology. Geochimica, 29(6): 513-520.
Lin B, Chen YC, Tang JX, Wang Q, Song Y, Yang C, Wang WL, He W and Zhang LJ. 2017a. 40Ar/39Ar and Rb-Sr ages of the Tiegelongnan porphyry Cu-(Au) deposit in the Bangong Co-Nujiang metallogenic belt of Tibet, China:Implication for generation of super-large deposit. Acta Geologica Sinica, 91(2): 602-616. DOI:10.1111/acgs.2017.91.issue-2
Lin B, Tang JX, Chen YC, Song Y, Hall G, Wang Q, Yang C, Fang X, Duan JL, Yang HH, Liu ZB, Wang YY and Feng J. 2017b. Geochronology and genesis of the Tiegelongnan porphyry Cu(Au) deposit in Tibet:Evidence from U-Pb, Re-Os dating and Hf, S, and H-O isotopes. Resource Geology, 67(1): 1-21. DOI:10.1111/rge.2017.67.issue-1
Lin B, Tang JX, Song Y, Wang Q, Cao HW, Baker M, Zhang LJ and Li YB. 2017. Petrogeochemistry and SHRIMP dating of Ga'erqin pillow basalt in Duolong, Northern Tibet. Acta Geoscientica Sinica, 38(5): 702-710.
Lin B, Tang JX, Chen YC, Baker M, Song Y, Yang HH, Wang Q, He W and Liu ZB. 2019. Geology and geochronology of Naruo large porphyry-breccia Cu deposit in the Duolong district, Tibet. Gondwana Research, 66: 168-182. DOI:10.1016/j.gr.2018.07.009
Liu YS, Hu ZC, Gao S, Günther D, Xu J, Gao CG and Chen HH. 2008. In situ analysis of major and trace elements of anhydrous minerals by LA-ICP-MS without applying an internal standard. Chemical Geology, 257(1-2): 34-43. DOI:10.1016/j.chemgeo.2008.08.004
Liu YS, Hu ZC, Zong KQ, Gao CG, Gao S, Xu J and Chen HH. 2010. Reappraisement and refinement of zircon U-Pb isotope and trace element analyses by LA-ICP-MS. Chinese Science Bulletin, 55(15): 1535-1546. DOI:10.1007/s11434-010-3052-4
Ma CQ, Yang KG, Tang ZH and Li ZT. 1994. Magma-Dynamics of Granitoids Theory, Method and A Case Study of the Eastern Hubei Granitoids. Wuhan: China University of Geosciences Press.
Ma LT, Zhang ZC, Dong SY, Zhang S, Zhang DY and Huang H. 2010. Geology and geochemistry of the Yingmailai granitic intrusion in the Southern Tianshan and its implications. Earth Science (Journal of China University of Geosciences), 35(6): 908-920. DOI:10.3799/dqkx.2010.106
Maas R, Nicholls IA and Legg C. 1997. Igneous and metamorphic enclaves in the S-type Deddick granodiorite, Lachlan Fold Belt, SE Australia:Petrographic, geochemical and Nd-Sr isotopic evidence for crustal melting and magma mixing. Journal of Petrology, 38(7): 815-841. DOI:10.1093/petroj/38.7.815
McDonough WF, Sun SS, Ringwood AE, Jagoutz E and Hofmann AW. 1992. Potassium, rubidium, and cesium in the Earth and Moon and the evolution of the mantle of the Earth. Geochimica et Cosmochimica Acta, 56(3): 1001-1012. DOI:10.1016/0016-7037(92)90043-I
Mo XX. 2011. Magma and magmatic/igneous rocks:A lithoprobe into the deep earth and records of the earth's evolution. Chinese Journal of Nature, 33(5): 255-259.
Montel JM and Vielzeuf D. 1997. Partial melting of metagreywackes, Part Ⅱ. Compositions of minerals and melts. Contributions to Mineralogy and Petrology, 128(2-3): 176-196. DOI:10.1007/s004100050302
Noyes HJ, Frey FA and Wones DR. 1983. A tale of two plutons:Geochemical evidence bearing on the origin and differentiation of the Red Lake and Eagle Peak plutons, Central Sierra Nevada, California. The Journal of Geology, 91(5): 487-509. DOI:10.1086/628801
Paterson SR, Pignotta GS and Vernon RH. 2004. The significance of microgranitoid enclave shapes and orientations. Journal of Structural Geology, 26(8): 1465-1481. DOI:10.1016/j.jsg.2003.08.013
Pearce JA. 1996. Sources and settings of granitic rocks. Episodes, 19(4): 120-125.
Perugini D, Poli G, Christofides G and Eleftheriadis G. 2003. Magma mixing in the Sithonia plutonic complex, Greece:Evidence from mafic microgranular enclaves. Mineralogy and Petrology, 78(3-4): 173-200. DOI:10.1007/s00710-002-0225-0
Qi YQ, Hu RZ, Liu S and Tian JJ. 2008. Review on magma mixing and mingling. Bulletin of Mineralogy, Petrology and Geochemistry, 27(4): 409-416.
Qiangba ZX, Xie YW, Wu YW, Xie CM, Li QL and Qiu JQ. 2009. Zircon SIMS U-Pb dating and its significance of cumulate gabbro from Dengqen ophiolite, eastern Tibet, China. Geological Bulletin of China, 28(9): 1253-1258.
Qiangba ZX, Wu H, Gesang WD, Ciren OZ, Basang DZ, Qiong D and Niu DW. 2016. Early Cretaceous magmatism in Dongqiao, Tibet:Implications for the evolution of the Bangong-Nujiang Ocean and crustal growth in a continent-continent collision zone. Geological Bulletin of China, 35(5): 648-666.
Qiu JS, Wang RQ, Zhao JL and Yu SB. 2015. Petrogenesis of the Early Jurassic gabbro-granite complex in the middle segment of the Gangdese belt and its implications for tectonic evolution of Neo-Tethys:A case study of the Dongga pluton in Xi'gaze. Acta Petrologica Sinica, 31(22): 3569-3580.
Qu XM, Wang RJ, Xin HB, Zhao YY and Fan XT. 2009. Geochronology and geochemistry of igneous rocks related to the subduction of the Tethys oceanic plate along the Bangong Lake arc zone, the western Tibetan Plateau. Geochimica, 38(6): 523-535.
Ren JS, Niu BG and Liu ZG. 1999. Soft collision, superposition orogeny and polycyclic suturing. Earth Science Frontiers, 6(3): 85-93.
Rickwood PC. 1989. Boundary lines within petrologic diagrams which use oxides of major and minor elements. Lithos, 22(4): 247-263. DOI:10.1016/0024-4937(89)90028-5
Rudnick RL and Gao S. 2003. Composition of the continental crust. In: Rudniek RL (ed.). Treatise on Geochemistry. Amsterdam: Elsevier, 1-64 http://www.sciencedirect.com/science/article/pii/0016703795000382
Rushmer T. 1991. Partial melting of two amphibolites:Contrasting experimental results under fluid-absent conditions. Contributions to Mineralogy and Petrology, 107(1): 41-59. DOI:10.1007/BF00311184
Shamberger PJ and Garcia MO. 2007. Geochemical modeling of magma mixing and magma reservoir volumes during early episodes of Kīlauea Volcano's Pu'u 'Ō'ō eruption. Bulletin of Volcanology, 69(4): 345-352. DOI:10.1007/s00445-006-0074-5
Shellnutt JG, Jahn BM and Dostal J. 2010. Elemental and Sr-Nd isotope geochemistry of microgranular enclaves from peralkaline A-type granitic plutons of the Emeishan large igneous province, SW China. Lithos, 119(1): 34-46.
Shi RD, Griffin WL, O'Reilly SY, Huang QS, Zhang XR, Liu DL, Zhi XC, Xia QX and Ding L. 2012. Melt/mantle mixing produces podiform chromite deposits in ophiolites:Implications of Re-Os systematics in the Dongqiao Neo-Tethyan ophiolite, northern Tibet. Gondwana Research, 21(1): 194-206. DOI:10.1016/j.gr.2011.05.011
Song B, Zhang YH, Wan YS and Jian P. 2002. Mount making and procedure of the SHRIMP dating. Geological Review, 48(Suppl.): 26-30.
Song Y, Qu XM, Xin HB, Wang ZY and Du DD. 2013. The crustal carbon cycling in the Late Cretaceous extension of Bangong Lake arc zone in the Xizang (Tibetan) Plateau. Geological Review, 59(2): 225-233.
Song Y, Yang HH, Lin B, Liu ZB, Wang Q, Gao K, Yang C and Fang X. 2017. The preservation system of epithermal deposits in South Qiangtang terrane of central Tibetan Plateau and its significance:A case study of the Tiegelongnan superlarge deposit. Acta Geoscientica Sinica, 38(5): 659-669.
Song Y, Yang C, Wei SG, Yang HH, Fang X and Lu HT. 2018. Tectonic control, reconstruction and preservation of the Tiegelongnan porphyry and epithermal overprinting Cu (Au) deposit, Central Tibet, China. Minerals, 8(9): 398. DOI:10.3390/min8090398
Spring L, Bussy F, Vannay J C, Huon S and Cosca MA. 1993. Early Permian granitic dykes of alkaline affinity in the Indian High Himalaya of upper Lahul and SE Zanskar:Geochemical characterization and geotectonic implications. Geological Society, London, Special Publications, 74(1): 251-264. DOI:10.1144/GSL.SP.1993.074.01.18
Stampfli GM and Borel GD. 2002. A plate tectonic model for the Paleozoic and Mesozoic constrained by dynamic plate boundaries and restored synthetic oceanic isochrons. Earth and Planetary Science Letters, 196(1-2): 17-33. DOI:10.1016/S0012-821X(01)00588-X
Tang JX, Sun XG, Ding S, Wang Q, Wang YY, Yang C, Chen HQ, Li YB, Li YB, Wei LJ, Zhang Z, Song JL, Yang HH, Duan JL, Gao K, Fang X and Tan JY. 2014a. Discovery of the epithermal deposit of Cu (Au-Ag) in the Duolong ore concentrating area, Tibet. Acta Geoscientica Sinica, 35(1): 6-10.
Tang JX, Wang Q, Yang C, Ding S, Lang XH, Liu HF, Huang Y, Zheng WB, Wang LQ, Gao YM, Feng J, Duan JL, Song Y, Wang YY, Lin B, Fang X, Zhang Z and Yang HH. 2014b. Two porphyry-epithermal deposit metallogenic subseries in Tibetan Plateau:Practice of "absence prospecting" deposit metallogenic series. Mineral Deposits, 33(6): 1151-1170.
Tang JX, Song Y, Wang Q, Lin B, Yang C, Guo N, Fang X, Yang HH, Wang YY, Gao K, Ding S, Zhang Z, Duan JL, Chen HQ, Su DK, Feng J, Liu ZB, Wei SG, He W, Song JL, Li YB and Wei LJ. 2016. Geological characteristics and exploration model of the Tiegelongnan Cu (Au-Ag) deposit:The first ten million tons metal resources of a porphyry-epithermal deposit in Tibet. Acta Geoscientica Sinica, 37(6): 663-690.
Tang JX, Wang Q, Yang HH, Gao X, Zhang ZB and Zou B. 2017. Mineralization, exploration and resource potential of porphyry-skarn-epithermal copper polymetallic deposits in Tibet. Acta Geoscientica Sinica, 38(5): 571-613.
Troll VR and Schmincke HU. 2002. Magma mixing and crustal recycling recorded in ternary feldspar from compositionally zoned peralkaline ignimbrite 'A', Gran Canaria, Canary Islands. Journal of Petrology, 43(2): 243-270. DOI:10.1093/petrology/43.2.243
Vernon RH. 1983. Restite, xenoliths and microgranitoid enclaves in granites. Journal and Proceedings, Royal Society of New South Wales, 116(3): 77-103.
Vernon RH. 2014. Microstructures of microgranitoid enclaves and the origin of S-type granitoids. Australian Journal of Earth Sciences, 61(2): 227-239. DOI:10.1080/08120099.2014.886623
Wang BD, Wang LQ, Chung SL, Chen JL, Yin FG, Liu H, Li XB and Chen LK. 2016. Evolution of the Bangong-Nujiang Tethyan Ocean:Insights from the geochronology and geochemistry of mafic rocks within ophiolites. Lithos, 245: 18-33. DOI:10.1016/j.lithos.2015.07.016
Wang DZ, Zhou XM, Xu XS and Yao YP. 1992. On geneses of microgranitoid enclaves. Journal of Guilin College of Geology, 12(3): 235-241.
Wang DZ and Xie L. 2008. Magma mingling:Evidence from enclaves. Geological Journal of China Universities, 14(1): 16-21.
Wang Q, Lin B, Tang JX, Song Y, Li YB, Hou JF, Li YB and Wei LJ. 2018. Diagenesis, lithogenesis and geodynamic setting of intrusions in Senadong area, Duolong district, Tibet. Earth Science, 43(4): 1125-1141.
Wang T. 2000. Origin of hybrid granitoids and the implications for continental dynamics. Acta Petrologica Sinica, 16(2): 161-168.
Wang YW, Wang JB, Long LL, Zou T, Tang PZ and Wang LJ. 2012. Type, indicator, mechanism, model and relationship with mineralization of magma mixing:A case study in North Xinjiang. Acta Petrologica Sinica, 28(8): 2317-2330.
Wei SG, Song Y, Tang JX, Gao K, Feng J, Li YB and Hou L. 2016. Geochronology, geochemistry and petrogenesis of quartz diorite porphyrite from the Sena copper (gold) deposit, Tibet. Geology in China, 43(6): 1894-1912.
Wei SG, Tang JX, Song Y, Liu ZB, Wang Q, Lin B, He W and Feng J. 2017a. Petrogenesis, zircon U-Pb geochronology and Sr-Nd-Hf isotopes of the intermediate-felsic volcanic rocks from the Duolong deposit in the Bangonghu-Nujiang suture zone, Tibet, and its tectonic significance. Acta Geologica Sinica, 91(1): 132-150.
Wei SG, Tang JX, Song Y, Liu ZB, Wang Q, Lin B, Wang YY, Fang X, Yang HH and Feng J. 2017b. Zircons LA-MC-ICP-MS U-Pb ages, petrochemical, petrological and its significance of the Potassic monzonitic granite porphyry from the Duolong ore concentrated district, Gaize County, Xizang (Tibet). Geological Review, 63(1): 189-206.
Winther KT. 1996. An experimentally based model for the origin of tonalitic and trondhjemitic melts. Chemical Geology, 127(1-3): 43-59. DOI:10.1016/0009-2541(95)00087-9
Wu YB and Zheng YF. 2004. Genesis of zircon and its constraints on interpretation of U-Pb age. Chinese Science Bulletin, 49(15): 1554-1569. DOI:10.1007/BF03184122
Xiao QH, Deng JF, Ma DQ, Hong DW, Mo XX, Lu XX, Li ZC, Wang XW, Ma CQ, Wu FY, Luo ZH and Wang T. 2002. The Ways of Investigation on Granitoids. Beijing: Geological Publishing House.
Xu JF, Mei HJ, Yu XY, Bai ZH, Niu HC, Chen FR, Zheng ZP and Wang Q. 2001. Adakites related to subduction in the northern margin of Junggar arc for the Late Paleozoic:Products of slab melting. Chinese Science Bulletin, 46: 684-688.
Yu JJ, Hou XG, Ge WC, Zhang YL and Liu JC. 2013. Magma mixing genesis of the Early Permian Liulian pluton at the northeastern margin of the Jiamusi massif in NE China:Evidences from petrography, geochronology and geochemistry. Acta Petrologica Sinica, 29(9): 2971-2986.
Zhang C, Liu CZ, Wu FY, Zhang LL and Ji WQ. 2016. Geochemistry and geochronology of mafic rocks from the Luobusa ophiolite, South Tibet. Lithos, 245: 93-108. DOI:10.1016/j.lithos.2015.06.031
Zhang DQ. 1981. The characteristics and origin of the dioritic inclusions in the Fangshan intrusion, Beijing. Geological Review, 27(2): 143-149.
Zhang Q, Pan GQ, Li CD, Jin WJ and Jia XQ. 2007. Granitic magma mixing versus basaltic magma mixing:New viewpoints on granitic magma mixing process:Some crucial questions on granite study (1). Acta Petrologica Sinica, 23(5): 1141-1152.
Zhang Z, Fang X, Tang JX, Wang Q, Yang C, Wang YY, Ding S and Yang HH. 2017. Chronology, geochemical characteristics of the Gaerqin porphyry copper deposit in the Duolong ore concentration area in Tibet and discussion about the identification of the lithoscaps and the possible epithermal deposit. Acta Petrologica Sinica, 33(2): 476-494.
Zhao YJ, Yuan C, Zhou MF, Yan DP, Long XP and Li JG. 2007. Geochemistry and petrogenesis of Laojungou and Mengtonggou granites in western Sichuan, China:Constraints on the nature of Songpan-Ganzi basement. Acta Petrologica Sinica, 23(5): 995-1006.
Zhou WG, Xie HS, Liu YG, Zheng XG, Zhao ZD and Zhou H. 2005. Dehydration melting of solid amphibolite at 2.0GPa:Effects of time and temperature. Science in China (Series D), 48(8): 1120-1133. DOI:10.1360/01yd0534
Zhu DC, Mo XX, Zhao ZD, Niu Y, Wang LQ, Chu QH, Pan GT, Xu JF and Zhou CY. 2010. Presence of Permian extension-and arc-type magmatism in southern Tibet:Paleogeographic implications. Geological Society of America Bulletin, 122(7-8): 979-993. DOI:10.1130/B30062.1
陈斌, 刘超群, 田伟. 2006. 太行山中生代岩浆作用过程中的壳幔岩浆混合作用:岩石学和地球化学证据. 地学前缘, 13(2): 140-147. DOI:10.3321/j.issn:1005-2321.2006.02.012
陈广俊, 孙丰月, 李碧乐, 王冠, 钱烨, 许庆林, 刘国才. 2014. 东昆仑沟里地区暗色包体及其寄主岩石地球化学特征及成因. 吉林大学学报(地球科学版), 44(3): 892-904.
陈国超, 裴先治, 李瑞保, 李佐臣, 刘成军, 陈有炘, 徐通, 张永明. 2016. 东昆仑东段香加南山花岗岩基的岩浆混合成因:来自镁铁质微粒包体的证据. 地学前缘, 23(4): 226-240.
丁帅, 唐菊兴, 郑文宝, 杨超, 张志, 王勤, 王艺云. 2017. 西藏拿若斑岩型铜(金)矿含矿岩体年代学、地球化学及地质意义. 地球科学, 42(1): 1-23. DOI:10.3969/j.issn.1672-6561.2017.01.001
董国臣, 莫宣学, 赵志丹, 朱弟成, 王亮亮, 陈涛, 李冰. 2006. 冈底斯岩浆带中段岩浆混合作用:来自花岗杂岩的证据. 岩石学报, 22(4): 835-844.
杜德道, 曲晓明, 王根厚, 辛洪波, 刘治博. 2011. 西藏班公湖-怒江缝合带西段中特提斯洋盆的双向俯冲:来自岛弧型花岗岩锆石U-Pb年龄和元素地球化学的证据. 岩石学报, 27(7): 1993-2002.
付强, 葛文胜, 温长顺, 蔡克勤, 李世富, 张志伟, 李小飞. 2011. 广西米场花岗岩及其暗色微粒包体的地球化学特征和成因分析. 地球学报, 32(3): 293-303. DOI:10.3975/cagsb.2011.03.04
顾枫华, 章永梅, 刘瑞萍, 郑硌, 孙玄. 2015. 内蒙古沙德盖花岗岩岩浆混合作用:岩相学、矿物化学和年代学证据. 岩石学报, 31(5): 1374-1390.
郝百武. 2012. 内蒙古镶黄旗南晚古生代晶洞钾长花岗岩杂岩的发现、岩石成因及其构造意义. 吉林大学学报(地球科学版), 42(2): 269-284.
黄启帅, 史仁灯, 丁炳华, 刘德亮, 张晓冉, 樊帅权, 支霞臣. 2012. 班公湖MOR型蛇绿岩Re-Os同位素特征对班公湖-怒江特提斯洋裂解时间的制约. 岩石矿物学杂志, 31(4): 465-478. DOI:10.3969/j.issn.1000-6524.2012.04.001
金性春, 于开平. 2003. 俯冲工厂和大陆物质的俯冲再循环研究. 地球科学进展, 18(5): 737-744. DOI:10.3321/j.issn:1001-8166.2003.05.014
康磊, 李永军, 张兵, 张洪伟, 王军年. 2009. 新疆西准噶尔夏尔莆岩体岩浆混合的岩相学证据. 岩石矿物学杂志, 28(5): 423-432. DOI:10.3969/j.issn.1000-6524.2009.05.002
雷传扬, 李志军, 张志, 胡正华, 王红星, 宋俊龙. 2012. 西藏尕尔穷铜金矿床花岗岩类的地球化学特征及其地球动力学意义. 地球学报, 33(4): 601-612. DOI:10.3975/cagsb.2012.04.20
雷传扬, 吴建亮, 尹显科, 刘文, 王波, 李威, 袁华云, 张伟, 尹滔, 裴亚伦. 2018. 班公湖-怒江缝合带西段沙木罗组地层中闪长玢岩脉的发现及其地质意义. 矿物岩石地球化学通报, 37(2): 250-259.
雷传扬, 吴建亮, 尹显科, 刘文, 王波, 李威, 袁华云, 张伟, 尹滔, 裴亚伦. 2019. 班公湖-怒江缝合带西段阿翁错地区中酸性侵入岩的成因及其对区域构造演化的指示. 大地构造与成矿学, 43(1): 168-185.
雷鸣, 陈建林, 许继峰, 曾云川. 2015. 拉萨地体中北部尕尔穷晚白垩世早期高镁闪长玢岩地球化学特征指示:加厚下地壳的拆沉?. 地质通报, 34(2-3): 337-346.
李昌年. 1999. 赣东北前寒武纪岩浆混合作用岩带的发现及其地质意义. 地学前缘, 6(4): 331-337. DOI:10.3321/j.issn:1005-2321.1999.04.016
李发桥, 刘治博, 唐菊兴, 宋扬, 高轲, 李海峰, 赵志强, 滕磊. 2018. 西藏玛日埃错地区花岗斑岩岩石成因及其对班公湖-怒江缝合带中段演化的制约. 地球科学, 43(4): 1051-1069.
李光明, 张夏楠, 秦克章, 孙兴国, 赵俊兴, 印贤波, 李金祥, 袁华山. 2015. 羌塘南缘多龙矿集区荣那斑岩-高硫型浅成低温热液Cu-(Au)套合成矿:综合地质、热液蚀变及金属矿物组合证据. 岩石学报, 31(8): 2307-2324.
李华亮, 高成, 李正汉, 张璋, 彭智敏, 关俊雷. 2016. 西藏班公湖地区竟柱山组时代及其构造意义. 大地构造与成矿学, 40(4): 663-673.
李金祥, 秦克章, 李光明. 2006. 富金斑岩型铜矿床的基本特征、成矿物质来源与成矿高氧化岩浆-流体演化. 岩石学报, 22(3): 678-688.
李金祥. 2008.班公湖带多不杂超大型富金斑岩铜矿床的成岩成矿年代学、岩石学及高氧化岩浆-流体-成矿作用.博士学位论文.北京: 中国科学院地质与地球物理研究所, 1-224 http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y1625435
李献华, 周汉文, 刘颖, 李寄嵎, 陈正宏, 于津生, 桂训唐. 2000. 粤西阳春中生代钾玄质侵入岩及其构造意义:Ⅰ.岩石学和同位素地质年代学. 地球化学, 29(6): 513-520. DOI:10.3321/j.issn:0379-1726.2000.06.001
林彬, 唐菊兴, 宋扬, 王勤, 曹华文, Baker M, 张乐骏, 李玉彬. 2017. 藏北多龙矿集区尕尔勤枕状玄武岩地球化学及SHRIMP测年. 地球学报, 38(5): 702-710.
马昌前, 杨坤光, 唐仲华, 李增田. 1994. 花岗岩类岩浆动力学:理论方法及鄂东花岗岩类例析. 武汉: 中国地质大学出版社.
马乐天, 张招崇, 董书云, 张舒, 张东阳, 黄河. 2010. 南天山英买来花岗岩的地质、地球化学特征及其地质意义. 地球科学-中国地质大学学报, 35(6): 908-920.
莫宣学. 2011. 岩浆与岩浆岩:地球深部"探针"与演化记录. 自然杂志, 33(5): 255-259.
齐有强, 胡瑞忠, 刘燊, 田建吉. 2008. 岩浆混合作用研究综述. 矿物岩石地球化学通报, 27(4): 409-416. DOI:10.3969/j.issn.1007-2802.2008.04.012
强巴扎西, 谢尧武, 吴彦旺, 解超明, 李秋立, 邱军强. 2009. 藏东丁青蛇绿岩中堆晶辉长岩锆石SIMS U-Pb定年及其意义. 地质通报, 28(9): 1253-1258. DOI:10.3969/j.issn.1671-2552.2009.09.013
强巴扎西, 吴浩, 格桑旺堆, 次仁欧珠, 巴桑顿珠, 琼达, 女达娃. 2016. 班公湖-怒江缝合带中段东巧地区早白垩世岩浆作用——对大洋演化和地壳增厚的指示. 地质通报, 35(5): 648-666. DOI:10.3969/j.issn.1671-2552.2016.05.003
邱检生, 王睿强, 赵姣龙, 喻思斌. 2015. 冈底斯中段早侏罗世辉长岩-花岗岩杂岩体成因及其对新特提斯构造演化的启示:以日喀则东嘎岩体为例. 岩石学报, 31(12): 3569-3580.
曲晓明, 王瑞江, 辛洪波, 赵元艺, 樊兴涛. 2009. 西藏西部与班公湖特提斯洋盆俯冲相关的火成岩年代学和地球化学. 地球化学, 38(6): 523-525. DOI:10.3321/j.issn:0379-1726.2009.06.002
任纪舜, 牛宝贵, 刘志刚. 1999. 软碰撞、叠覆造山和多旋回缝合作用. 地学前缘, 6(3): 85-93. DOI:10.3321/j.issn:1005-2321.1999.03.008
宋彪, 张玉海, 万渝生, 简平. 2002. 锆石SHRIMP样品靶制作、年龄测定及有关现象讨论. 地质论评, 48(增): 26-30.
宋扬, 曲晓明, 辛洪波, 王志云, 杜德道. 2013. 西藏班公湖岛弧带晚白垩世地壳伸展期间碳的壳内循环. 地质论评, 59(2): 225-233. DOI:10.3969/j.issn.0371-5736.2013.02.005
宋扬, 杨欢欢, 林彬, 刘治博, 王勤, 高轲, 杨超, 方向. 2017. 青藏高原羌塘地体南缘浅成低温热液成矿系统的保存机制及其重要意义——以铁格隆南超大型矿床为例. 地球学报, 38(5): 659-669.
唐菊兴, 孙兴国, 丁帅, 王勤, 王艺云, 杨超, 陈红旗, 李彦波, 李玉彬, 卫鲁杰, 张志, 宋俊龙, 杨欢欢, 段吉琳, 高轲, 方向, 谭江云. 2014a. 西藏多龙矿集区发现浅成低温热液型铜(金银)矿床. 地球学报, 35(1): 6-10.
唐菊兴, 王勤, 杨超, 丁帅, 郎兴海, 刘鸿飞, 黄勇, 郑文宝, 王立强, 高一鸣, 冯军, 段吉琳, 宋杨, 王艺云, 林彬, 方向, 张志, 杨欢欢. 2014b. 青藏高原两个斑岩-浅成低温热液矿床成矿亚系列及其"缺位找矿"之实践. 矿床地质, 33(6): 1151-1170.
唐菊兴, 宋扬, 王勤, 林彬, 杨超, 郭娜, 方向, 杨欢欢, 王艺云, 高轲, 丁帅, 张志, 段吉琳, 陈红旗, 粟登逵, 冯军, 刘治博, 韦少港, 贺文, 宋俊龙, 李彦波, 卫鲁杰. 2016. 西藏铁格隆南铜(金银)矿床地质特征及勘查模型——西藏首例千万吨级斑岩-浅成低温热液型矿床. 地球学报, 37(6): 663-690.
唐菊兴, 王勤, 杨欢欢, 高昕, 张泽斌, 邹兵. 2017. 西藏斑岩-矽卡岩-浅成低温热液铜多金属矿成矿作用、勘查方向与资源潜力. 地球学报, 38(5): 571-613.
王德滋, 周新民, 徐夕生, 姚玉鹏. 1992. 微粒花岗岩类包体的成因. 桂林冶金地质学院学报, 12(3): 235-241.
王德滋, 谢磊. 2008. 岩浆混合作用:来自岩石包体的证据. 高校地质学报, 14(1): 16-21. DOI:10.3969/j.issn.1006-7493.2008.01.002
王勤, 林彬, 唐菊兴, 宋扬, 李彦波, 侯俊富, 李玉彬, 卫鲁杰. 2018. 多龙矿集区色那东岩体年龄、成因与动力学背景. 地球科学, 43(4): 1125-1141.
王涛. 2000. 花岗岩混合成因研究及大陆动力学意义. 岩石学报, 16(2): 161-168.
王玉往, 王京彬, 龙灵利, 邹滔, 唐萍芝, 王莉娟. 2012. 岩浆混合作用的类型、标志、机制、模式及其与成矿的关系——以新疆北部为例. 岩石学报, 28(8): 2317-2330.
韦少港, 宋扬, 唐菊兴, 高轲, 冯军, 李彦波, 侯淋. 2016. 西藏色那铜(金)矿床石英闪长玢岩年代学、地球化学与岩石成因. 中国地质, 43(6): 1894-1912.
韦少港, 唐菊兴, 宋扬, 刘治博, 王勤, 林彬, 贺文, 冯军. 2017a. 西藏班公湖-怒江缝合带美日切错组中酸性火山岩锆石U-Pb年龄、Sr-Nd-Hf同位素、岩石成因及其构造意义. 地质学报, 91(1): 132-150.
韦少港, 唐菊兴, 宋扬, 刘治博, 王勤, 林彬, 王艺云, 方向, 杨欢欢, 冯军. 2017b. 西藏改则多龙矿集区地堡那木岗矿床钾玄质二长花岗斑岩锆石LA-MC-ICP-MS U-Pb年龄、地球化学特征及其地质意义. 地质论评, 63(1): 189-206.
吴元保, 郑永飞. 2004. 锆石成因矿物学研究及其对U-Pb年龄解释的制约. 科学通报, 49(16): 1589-1604. DOI:10.3321/j.issn:0023-074X.2004.16.002
肖庆辉, 邓晋福, 马大铨, 洪大卫, 莫宣学, 卢欣祥, 李志昌, 汪雄武, 马昌前, 吴福元, 罗照华, 王涛. 2002. 花岗岩研究思维与方法. 北京: 地质出版社.
许继峰, 梅厚钧, 于学元, 白正华, 牛贺才, 陈繁荣, 郑作平, 王强. 2001. 准噶尔北缘晚古生代岛弧中与俯冲作用有关的adakite火山岩:消减板片部分熔融的产物. 科学通报, 46(8): 684-688. DOI:10.3321/j.issn:0023-074X.2001.08.016
于介江, 侯雪刚, 葛文春, 张彦龙, 柳佳成. 2013. 佳木斯地块东北缘早二叠世六连岩体的岩浆混合成因:岩相学、年代学和地球化学证据. 岩石学报, 29(9): 2971-2986.
张德全. 1981. 北京房山侵入体闪长质包体的特征及其成因探讨. 地质论评, 27(2): 143-149. DOI:10.3321/j.issn:0371-5736.1981.02.008
张旗, 潘国强, 李承东, 金惟俊, 贾秀勤. 2007. 花岗岩混合问题:与玄武岩对比的启示——关于花岗岩研究的思考之一. 岩石学报, 23(5): 1141-1152.
张志, 方向, 唐菊兴, 王勤, 杨超, 王艺云, 丁帅, 杨欢欢. 2017. 西藏多龙矿集区尕尔勤斑岩铜矿床年代学及地球化学——兼论硅帽的识别与可能的浅成低温热液矿床. 岩石学报, 33(2): 476-494.
赵永久, 袁超, 周美夫, 颜丹平, 龙晓平, 李继亮. 2007. 川西老君沟和孟通沟花岗岩的地球化学特征、成因机制及对松潘-甘孜地体基底性质的制约. 岩石学报, 23(5): 995-1006.