岩石学报  2020, Vol. 36 Issue (7): 2127-2148, doi: 10.18654/1000-0569/2020.07.14   PDF    
青藏高原东南缘腾冲后碰撞粗面安山岩形成的深部岩浆过程:来自辉长质物质的启示
成智慧1,2, 杨志军1,2, 赵文斌3, 张茂亮4, 雷鸣3, 马琳3, 李菊景3     
1. 中山大学地球科学与工程学院, 广州 510275;
2. 广东省地质过程与矿产资源探查重点实验室, 广州 510275;
3. 中国科学院地质与地球物理研究所, 新生代地质与环境重点实验室, 北京 100029;
4. 天津大学表层地球系统科学研究院, 天津 300072
摘要: 位于青藏高原东南缘的腾冲火山区发育大量的后碰撞高钾钙碱性岩浆岩(8Ma至今),形成了近连续的玄武岩-玄武安山岩-粗面安山岩-英安岩系列。在晚更新世(0.3~0.4Ma)粗面安山岩中发现了大量的辉长岩包体与辉长质矿物聚晶。辉长岩包体分为两类:I类为辉长苏长岩(直径2~4cm),主要由斜长石(50%~60%)、单斜辉石(20%~30%)和斜方辉石(5%~10%)组成,矿物间呈高角度接触的开放结构并包含少量的基质玻璃;Ⅱ类辉长岩包体(直径2~12cm)主要由斜长石(40%~50%)与单斜辉石(30%~40%)组成,含少量铁钛氧化物(5%~10%),矿物间见少量的基质玻璃与微晶斜方辉石(粒径 < 50μm)。矿物聚晶(粒径2~5mm)由斜长石(40%~60%)、单斜辉石(20%~30%)和斜方辉石(5%~10%)组成,矿物间呈高角度接触的开放结构,其间贯入基质玻璃。辉长岩包体、矿物聚晶与寄主粗面安山岩斑晶具有相同的矿物组成(斜长石+单斜辉石+斜方辉石);单斜辉石具有相似的微量元素组成特征;辉长岩包体全岩主量元素成分落在玄武岩-玄武安山岩-粗面安山岩-英安岩的演化序列中,它们与寄主粗面安山岩具有一致的Sr-Nd-Pb同位素组成。上述观测结果表明,辉长岩包体和矿物聚晶与寄主岩浆同源,它们均来自粗面安山质岩浆房。矿物温度与压力计算结果表明,辉长岩包体与矿物聚晶的结晶温度低于斑晶,形成深度位于粗面安山质岩浆房的中上部。这些具有开放结构的矿物聚晶与辉长岩包体可能代表粗面安山质岩浆喷发前形成晶粥的"碎片":矿物聚晶与I类辉长岩包体矿物间基质含量较高,矿物成分与寄主粗面安山质熔体(全岩成分)平衡,可能代表寄主粗面安山质岩浆形成的晶粥;Ⅱ类辉长岩包体矿物间紧密程度较高,矿物与较演化的岩浆平衡,可能代表早期较演化的岩浆形成的经历压实作用的晶粥。粗面安山质岩浆的快速上升将这些晶粥破碎并以辉长岩包体与矿物聚晶的方式运移至地表。岩浆房中基性岩浆的补给是导致安山质岩浆喷发的重要诱发机制。
关键词: 后碰撞粗面安山岩    辉长岩包体    矿物聚晶    温压估算    岩浆房系统    青藏高原东南缘    
Gabbroic xenoliths and glomerocrysts in the post-collisional trachyandestic rocks from Tengchong, SE Tibet: Implications for the magma chamber processes
CHENG ZhiHui1,2, YANG ZhiJun1,2, ZHAO WenBin3, ZHANG MaoLiang4, LEI Ming3, MA Lin3, LI JuJing3     
1. School of Earth Sciences and Engineering, Sun Yat-sen University, Guangzhou 510275, China;
2. Guangdong Provincial Key Lab of Geological Processes and Mineral Resource Survey, Guangzhou 510275, China;
3. Key Laboratory of Cenozoic Geology and Environment, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China;
4. Institute of Surface-Earth System Science, Tianjin University, Tianjin 300072, China
Abstract: Post-collisional high-K calc-alkaline magmatic activities (8Ma to present) with eruption of basalts,basaltic andesites,trachyandesites and dacites are extensively distributed in the Tengchong volcanic field,SE Tibet. Plenty of gabbroic xenoliths and glomerocrysts were entrained by the Late Pleistocene (0.3~0.4Ma) trachyandesites. There are two types of gabbroic xenoliths: Type I xenoliths (2~4cm in size) are in open structures and consist of plagioclase (50%~60%),clinopyroxene (20%~30%) and orthopyroxene (5%~10%) with interstitial glass; Type Ⅱ xenoliths (2~12cm in size) contain plagioclase (40%~50%),clinopyroxene (30%~40%) and Fe-Ti oxide (5%~10%) with interstitial orthopyroxene and glass. The glomerocrysts (2~5mm in size) are also in open structures,and they are composed of plagioclase (40%~60%),clinopyroxene (20%~30%) and orthopyroxene (5%~10%) with interstitial glass. The gabbroic xenoliths and glomerocrysts display similar mineral assemblage (plagioclase + clinopyroxene + orthopyroxene) with phenocrysts of the host lava. The clinopyroxenes from xenoliths,glomerocrysts and phenocrysts show similar trace element patterns. Furthermore,gabbroic xenoliths display similar whole-rock Sr,Nd and Pb isotopes to the host lava. We interpret these to indicate a genetic link between xenoliths (plus glomerocrysts) and the trachyandesitic magma,suggesting that the gabbroic xenoliths and glomerocrysts belong to the magmatic system in the Tengchong volcanic field. Results of P-T estimation indicate that the glomerocrysts and gabbroic xenoliths appear to have crystallized at the upper environment of the trachyandesitic magma reservoir. The glomerocrysts and gabbroic xenoliths with open structures contain residual interstitial glass,suggestting that they are fragments of partially solidified crystal mush: the glomerocrysts and Type I xenoliths which are in equilibrium with the melt appear to have crystallized from the host trachyandesitic melts; Type Ⅱ xenoliths may have crystallized from a more evolved magma based on the lower Mg# of the clinopyroxene. The crystal mush may have been disaggregated by the rapid ascending trachyandesitic magma and carried to the surface as glomerocrysts and gabbroic xenoliths. Recharge by mafic magmas in magma reservoir may be the main trigger for eruption.
Key words: Post-collisional trachyandesites    Gabbroic xenoliths    Glomerocrysts    P-T estimation    Magmatic plumbing system    SE Tibet    

作为大陆地壳的平均组成,安山岩的成因对于揭示陆壳形成与演化具有重要的意义(Taylor and McLennan, 1985; Rudnick and Gao, 2003; Kelemen and Behn, 2016)。安山质岩浆主要为玄武质熔体在地壳不同深度的岩浆房经历岩浆混合(Tepley et al., 2000; Reubi and Blundy, 2009; Kent et al., 2010; Laumonier et al., 2014; Fodor and Johnson, 2016)、分离结晶(FC)(Gill, 1981)或地壳混染与分离结晶(AFC)(Grove and Donnelly-Nolan, 1986; Grove et al., 2003; Annen et al., 2006; Price et al., 2016)等过程的产物。岩石学、矿物学和全岩地球化学等特征被广泛用于探讨安山岩的成因(Humphreys et al., 2006; Davidson et al., 2007; Price et al., 2016)。大量的研究成果表明,安山岩携带的超基性-基性岩包体为探究安山岩的成因提供了独特的窗口(Stamper et al., 2014; Price et al., 2016; Camejo-Harry et al., 2018; Cooper et al., 2019)。

辉长岩包体是安山岩中较为常见的一类包体,在俯冲带和洋岛等火山活动伴随的玄武安山-安山质喷发物中均发现了大量的辉长岩包体。前者如南安第斯火山带(Costa et al., 2002)、小安德烈斯火山带(Arculus and Wills, 1980; Stamper et al., 2014; Camejo-Harry et al., 2018; Cooper et al., 2019)、新西兰北岛火山(Price et al., 2016)等,后者如夏威夷火山(Fodor et al., 1993; Shamberger and Hammer, 2006)和基拉韦尔火山(Fodor and Moore, 1994)等。这些辉长岩包体可能代表岩浆上升过程中捕获的地壳组分,也可能是岩浆系统早期或同期堆晶或分离结晶的产物。地壳辉长岩包体通常具有高压麻粒岩的粒状变晶结构,矿物经常有明显的变形或破碎(Costa et al., 2002; Stamper et al., 2014; Price et al., 2016; Camejo-Harry et al., 2018)。来自岩浆房系统的辉长岩通常与寄主岩浆具有同系列的主微量元素组成和较一致的Sr-Nd-Pb同位素组成(Yang et al., 2012; Stamper et al., 2014; Price et al., 2016)。其中,具有嵌晶结构的辉长岩包体通常为岩浆分离结晶形成的堆晶体(Arculus and Wills, 1980; Fodor and Moore, 1994; Costa et al., 2002; Stamper et al., 2014),它们具有较寄主岩浆更低的SiO2和更高的CaO含量,堆晶矿物也具有更基性的成分(例如,单斜辉石较高的Mg#值, 斜长石较高的An含量)(Stamper et al., 2014; Camejo-Harry et al., 2018; Cooper et al., 2019)。(中粗)粒状结构的辉长岩包体自形程度较高,矿物间呈包含火山玻璃的开放结构(Fodor and Moore, 1994; Costa et al., 2002; Cashman et al., 2017),除辉长岩包体以外,寄主岩浆还携带大量辉长质的矿物聚晶(glomerocryst)与循环晶(antecryst),这些辉长质的物质通常为岩浆上升过程中捕获的岩浆房中的晶粥(Costa et al., 2002; Jerram and Martin, 2008; Cashman and Sparks, 2013; Edmonds et al., 2019; Lissenberg et al., 2019)。由此可见,安山岩及其携带的辉长岩包体的岩石学、矿物学、矿物化学和岩石地球化学特征为揭示安山岩的成因及深部演化过程提供了独特的信息。

受印度-欧亚大陆碰撞(65~60Ma; Wu et al., 2014; Hu et al., 2015)及印度板块持续向北俯冲的影响,青藏高原东南缘的腾冲地区发育了大量的后碰撞高钾钙碱性火山岩,形成了近连续的玄武岩-玄武安山岩-粗面安山岩-英安岩系列(Zhou et al., 2012; Guo et al., 2015; Cheng et al., 2020; Li et al., 2020)。目前,关于腾冲安山岩的成因仍具有较大的争议。岩石学与地球化学研究成果表明,腾冲安山岩主要为玄武质岩浆在地壳岩浆房经历分离结晶(FC)(樊祺诚等, 1999;赵勇伟和樊祺诚, 2010)、地壳混染和分离结晶(AFC)(Zou et al., 2017; Duan et al., 2019; Cheng et al., 2020)或岩浆混合(Gao et al., 2015)等过程的产物。腾冲粗面安山岩中包含较丰富的岩石包体,它们的岩石学、矿物学及地球化学特征为揭示地壳深部岩浆房结构、探讨安山质岩浆的成因与演化过程提供了良好的契机。有学者对粗面安山岩中的辉长岩(直径 < 2.5cm)和麻粒岩(直径3~5cm)包体进行了较全面的岩石学与矿物学研究(于红梅等, 2010;林木森等, 2014),但岩石地球化学研究工作还较少。众所周知,包体与寄主火山岩的主微量元素尤其是Sr-Nd-Pb同位素特征对限定包体的来源与深部岩浆过程具有重要意义。本次研究在腾冲粗面安山岩中同时发现了大量的辉长岩包体(直径2~13cm)与矿物聚晶(粒径2~5mm),首次对辉长岩包体与粗面安山岩开展了的全岩主微量元素和Sr-Nd-Pb同位素的对比研究,以及各类单斜辉石(辉长岩包体、聚晶和斑晶)原位微量元素特征的对比研究,结合矿物温压计算结果,试图限定腾冲安山质岩浆房的结构,进而探讨辉长岩包体与粗面安山岩的成因。

1 地质背景

腾冲地块位于青藏高原的东南缘,其东部以高黎贡韧性剪切带为界,西部止于密支那缝合带。在大地构造位置上,腾冲地块处于密支那与怒江缝合带之间,为青藏高原南部拉萨地块的东南延伸(Yin and Harrison, 2000; Searle et al., 2007;戚学祥等, 2011; Qi et al., 2019)(图 1)。因处于印度-欧亚板块碰撞的最前缘,这里发育强烈的构造变形、广泛的岩浆活动和强烈的水热活动。腾冲地区的断裂主要为NS或NNE走向,包括西部的梁河断裂、中部的大盈江断裂和东部的固东-腾冲断裂等(皇甫岗和姜朝松,2000;戚学祥等,2011)。区内出露的地层由老至新分别为:元古界高黎贡群变质岩,岩性主要为千枚岩、片岩与片麻岩,厚度大于1000m;石炭系勐洪群泥岩、板岩和白云岩,厚度约3500m;之后,地层缺失,直至古近系南林组陆相碎屑沉积不整合于元古界变质岩之上(陈廷方等,2003)。由于经历了中特提斯洋(怒江洋,Yin and Harrison, 2000;莫宣学和潘桂棠,2006)与新特提斯洋(密支那洋,莫宣学和潘桂棠,2006)的闭合和块体之间的碰撞,腾冲地块内中、新生代爆发了燕山期和喜马拉雅期两次大规模的岩浆活动(杨启军等,2006戚学祥等,2011)。新生代以来,区内以广泛发育高钾钙碱性火山活动为特征,在南北狭长的盆地中形成了70多座火山,熔岩出露面积达790km2(姜朝松,1998樊祺诚等,1999)。腾冲地区水热活动(喷沸泉、沸泉和温泉等)强烈,地震活动频繁(佟伟和章铭陶,1989上官志冠等,2004成智慧等, 2012, 2014; Zhang et al., 2016)。

图 1 腾冲火山区后碰撞高钾钙碱性火山岩分布图(据Mo et al., 2006; Guo et al., 2015; Cheng et al., 2020修改) 五角星代表粗面安山岩及辉长岩包体的采样位置 Fig. 1 Simplified geological map showing the distribution of the post-collisional high-K calc-alkaline (HKCA) magmatic rocks in the Tengchong volcanic field (modified after Mo et al., 2006, Guo et al., 2015; Cheng et al., 2020) The star denotes sampling location of trachyandesites with gabbroic xenoliths

腾冲后碰撞火山活动始于上新世(8Ma至今),形成了一系列的高钾钙碱性火山岩(Zhu et al., 1983;皇甫岗和姜朝松, 2000;李大明等, 2000; Zhou et al., 2012; Guo et al., 2015),最近的一次火山活动发生于公元1609年(徐弘祖, 1639)。根据火山岩的年代学、岩石学及地球化学特征,腾冲火山活动可被划分为四个期次:上新世(8~3Ma)玄武岩-粗面玄武岩-玄武质粗面安山岩系列,早更新世(2.7~0.8Ma)粗面安山岩-英安岩系列,晚更新世(0.8~0.2Ma)玄武岩-粗面玄武岩系列和全新世(< 0.2Ma)粗面玄武岩-粗面安山岩系列(Cheng et al., 2020)。关于原始玄武质岩浆成因的研究较为成熟,认为其形成于富集地幔的低程度部分熔融,地幔富集过程可能受到了新特提斯洋壳(从柏林等, 1994; Tian et al., 2018; Duan et al., 2019)、俯冲的印度陆壳(Zou et al., 2014; Guo et al., 2015; Xu et al., 2018b; Cheng et al., 2020)或俯冲的印度大洋岩石圈(90°E海岭)(Zhou et al., 2012; Liu et al., 2018)的影响。关于中酸性岩浆成因的研究较少,认为其为玄武质岩浆上升至地表的过程中经历分离结晶(樊祺诚等, 1999;赵勇伟和樊祺诚, 2010)、地壳混染和分离结晶(AFC)(Zou et al., 2017; Duan et al., 2019; Cheng et al., 2020)或岩浆混合(Gao et al., 2015)等过程的产物。地震与地球物理资料显示,腾冲地壳厚度约38~42km(周真恒等, 2001;王椿镛等, 2002;张龙等, 2015; Xu et al., 2018a),在地壳不同深度(9~30km)存在多个部分熔融体,可能为壳内岩浆房(秦嘉政等, 2000; Bai et al., 2001;姜枚等, 2012;谭捍东等, 2013; Xu et al., 2018b; Hua et al., 2019)。

2 岩相学特征

所有的样品均采自黑空山火山以北5km出露面积约500m2的露天采石场(25°16′N、98°30′E)(图 1图 2)。共采集了6件粗面安山岩与12件包含辉长岩包体的粗面安山岩样品,其中,包含辉长岩包体的火山岩样品均采自该采石场北部厚度约3.5m的熔岩流剖面(图 2a)。

图 2 辉长岩包体与寄主火山岩野外照片 (a、b)辉长岩包体野外出露点,厚度约3.5m的安山质熔岩流;(c)包含辉长岩包体的粗面安山岩;(d)粗面安山岩中发育大量白色的矿物聚晶;(e、f)分别为I类和II类辉长岩包体,两类包体与火山岩均具有清晰的接触线 Fig. 2 Macroscopic overview of representative gabbroic xenoliths entrained by the trachyandesites in the Tengchong volcanic field (a, b) trachyandesite lava flow with 3.5m thick which entrains gabbroic xenoliths; (c) trachyandesites which entrain gabbroic xenoliths; (d) trachyandesites with plenty of glomerocrysts; (e, f) gabbroic xenoliths of Type I and Type II, respectively, displaying clear boundary with trachyandesites

火山岩样品新鲜,气孔状构造,斑状结构(图 2d)。斑晶以斜长石(5%~10%)和单斜辉石(5%~8%)为主,含少量斜方辉石(3%~5%)(图 3g)。斜长石呈自形-半自形板状,韵律结构发育,少数核部发育筛状结构(图 3g)。单斜辉石与斜方辉石为自形-半自形粒状结构,简单双晶发育,少量发育环带结构(图 3g)。粗面安山岩中发育大量的矿物聚晶(15% ~20%),粒径2~5mm,主要由自形-半自形的斜长石(40%~60%)、单斜辉石(20%~30%)和斜方辉石(5%~10%)构成,含少量Fe-Ti氧化物(< 5%)(图 3g)。斜长石颗粒之间呈较高的角度接触并构成格架,单斜辉石与斜方辉石颗粒充填其中,在聚晶矿物间隙贯入少量的玻璃质与气孔(图 3h)。

图 3 辉长岩包体和矿物聚晶的显微结构与背散射图像照片 (a)I类辉长岩包体,粒状结构,与粗面安山岩平整接触,无反应边结构(单偏光下);(b)II类辉长岩包体,粒状结构,与粗面安山岩平整接触(正交偏光下);(c)和(e)I类辉长岩包体正交偏光镜下照片与BSE图像,主要矿物组成为斜长石、单斜辉石与斜方辉石,含少量Fe-Ti氧化物,矿物间贯入玻璃质;(d)和(f)II类辉长岩包体正交偏光镜下照片与BSE图像,主要由斜长石与单斜辉石组成,含少量Fe-Ti氧化物,矿物间充填少量玻璃质;(g)和(h)矿物聚晶单偏光镜下照片与BSE图像,斜长石、单斜辉石与斜方辉石构成开放结构,矿物间贯入熔体.矿物缩写:Ap-磷灰石;Pl-斜长石;Cpx-单斜辉石;Opx-斜方辉石;Fe-Ti-Fe-Ti氧化物 Fig. 3 Representative microphotographs of gabbroic xenoliths and glomerocrysts in trachyandesites (a) Type I gabbroic xenoliths with granular texture, which display clear boundary with trachyandesites; (b) Type II gabbroic xenoliths with granular texture, which display clear boundary with trachyandesites; (c) and (e) Type I gabbroic xenoliths composed of plagioclase (Pl), clinopyroxene (Cpx) and orthopyroxene (Opx) with interstitial glass; (d) and (f) Type II gabbroic xenolith which is composed of Pl and Cpx with interstitial orthopyroxene and glass; (g) and (h) glomerocrysts composed of Pl, Cpx and Opx with open structure. Mineral abbreviations: Ap-apatite; Fe-Ti-Fe-Ti oxide

辉长岩包体野外呈灰白色,直径2~12cm,包体新鲜,呈等粒状或椭圆状(图 2e, f)。包体与寄主粗面安山岩具有较清晰的界限(图 2e, f图 3a, b)。根据矿物组成与结构特征,辉长岩包体分为两类:I类辉长苏长岩(样品数,3个)(图 2e图 3a)和II类辉长岩(样品数,9个)(图 2f图 3b)。I类辉长岩包体粒径较小(直径2~4cm),主要的矿物组成为斜长石(50%~60%)+单斜辉石(20%~30%)+斜方辉石(5%~10%)+Fe-Ti氧化物(< 5%)(图 3c, e)。其中,斜长石呈半自形粒状(0.5~5mm)结构,复合双晶发育;单斜辉石与斜方辉石呈半自形粒状(0.2~2mm)结构;矿物间隙贯入少量的玻璃质(图 3c, e)。II类辉长岩包体较大(直径2~12cm),主要的矿物组成为斜长石(40%~50%)+单斜辉石(30%~40%)+ Fe-Ti氧化物(5%~10%)(图 3d, f)。其中,斜长石呈自形至半自形板状,粒径2~3mm,少数边部发育海绵结构;单斜辉石为半自形至他形粒状结构,粒径1~3mm;斜长石与单斜辉石矿物间隙贯入玻璃质,同时发育少量微晶斜方辉石(粒径≤50μm)(图 3d, f)。

3 分析方法

本次研究对粗面安山岩和辉长岩包体开展了全岩主微量元素和Sr-Nd-Pb同位素分析,并对包含辉长岩包体的粗面安山岩进行了K-Ar定年工作。同时,详细的测试了辉长岩、矿物聚晶和粗面安山岩斑晶矿物(斜长石、单斜辉石与斜方辉石)的主量元素组成与单斜辉石的原位微区微量元素组成。

利用薄片切割法将辉长岩包体与粗面安山岩分离,分别粉碎至200目,用于全岩地球化学成分测试工作。全岩主微量元素和Sr-Nd-Pb同位素测试均在中国科学院地质与地球物理研究所完成。全岩主量元素成分采用X-射线荧光光谱仪XRF-1500测试,分析精度为~1%(含量>10%)和~5%(含量 < 1.0%)。全岩微量元素分析利用Finnigan Element型ICP-MS(电感耦合等离子体质谱仪)完成,精度优于10%。Rb-Sr和Sm-Nd同位素比值测试采用同位素稀释法完成,测量仪器为德国Finnigan公司MAT-262型热电离质谱计。分别采用146Nd/144Nd=0.7219和86Sr/88Sr=0.1194对Nd和Sr同位素比值进行质量分馏校正。详细测试方法参见Guo et al. (2006)。Pb同位素分析利用德国Finnigan公司MAT-262型热电离质谱仪完成,分析精度优于0.005%,利用标样NBS981进行校正,详细的测试方法参见Wang et al. (2018)。将包含辉长岩包体的粗面安山岩粉碎至40~60目,利用重选、磁选等选矿方法初步将基质与斑晶分离,然后在显微镜下将单晶矿物尤其是斜长石从基质中分离,挑选基质3~5g,利用K-Ar法对粗面安山岩开展定年工作,该项测试在中国地震局地质研究所利用高精度MM-1200质谱完成,详细的测试方法参见Aslan et al. (2014)

矿物的主量元素分析在中国科学院地质与地球物理研究所电子探针实验室利用JXA8100电子探针完成,加速电压为15kV,电子束电流10nA,束斑直径3μm。单斜辉石原位微量元素分析在中国科学院地质与地球物理研究所利用LA-ICP-MS完成,采用的激光剥蚀孔径为40μm,激光脉冲为10Hz,数据处理过程中选择Ca作为内标元素,国际标样NIST610用于外标校正。详细的分析方法和流程同杨岳衡等(2009)

4 分析结果 4.1 K-Ar定年结果

基质K-Ar定年结果显示,粗面安山质岩浆喷发于晚更新世(0.29~0.38Ma)(表 1),属于腾冲新生代火山区第三期岩浆活动(0.8~0.2Ma; Cheng et al., 2020)。本期岩浆活动主要分布于腾冲火山区的中部及中北部地区,形成了以玄武岩-玄武质安山岩-粗面安山岩为主的火山岩系列(图 1)。

表 1 粗面安山岩基质K-Ar定年结果 Table 1 Whole-rock K-Ar dating results for glass matrix of the trachyandesitic rocks
4.2 全岩地球化学特征

寄主火山岩与II类辉长岩包体全岩的主量、微量元素和Sr-Nd-Pb同位素成分见表 2。寄主火山岩具有较高的SiO2(58.28%~59.29%),Al2O3(16.54%~17.20%),Fe2O3(5.91% ~6.44%),CaO (5.48%~5.84%)和MgO(3.33%~3.87%)含量(表 2图 4)。它们具有较高的总碱(Na2O + K2O)(6.63%~7.02%)含量(表 2),根据火山岩TAS分类图将其定名为粗面安山岩。这些火山岩具有较高的K2O(3.22%~3.44%)含量,属于高钾钙碱性系列(表 2图 4)。II类辉长岩包体具有较低的SiO2(50.41%~51.55%)含量,较高的Al2O3(18.45%~20.26%),Fe2O3(6.39%~9.53%),CaO(10.59%~10.64%)和MgO(6.22%~6.86%)含量(表 2)。在主量元素(SiO2、CaO、Al2O3和FeO)与MgO的相关图解中(图 4),II类辉长岩包体均落在腾冲近连续的玄武岩-玄武安山岩-粗面安山岩-英安岩的演化序列内。

表 2 粗面安山岩与辉长岩包体全岩主量(wt%)、微量元素(×10-6)和Sr-Nd-Pb同位素组成 Table 2 Table 2 Whole-rock major element (wt%), trace element (×10-6) and Sr-Nd-Pb isotopic compositions for trachyandesitic rocks and gabbroic xenoliths

图 4 腾冲后碰撞高钾钙碱性火山岩与II类辉长岩包体哈克图解 实心与空心数据点分别代表本文与前人发表的数据.前人数据来源:樊祺诚等, 1999;赵勇伟和樊祺诚, 2010; Li et al., 2012; Zhang et al., 2012; Zhou et al., 2012; Guo et al., 2015; Tian et al., 2018; Cheng et al., 2020.图 5图 6数据来源同此图 Fig. 4 Harker diagrams for the post-collisional HKCA magmas and Type II gabbroic xenoliths in the Tengchong volcanic field Filled and open symbols represent data from this study and previous publications, respectively. Sources of published data: Fan et al., 1999; Zhao and Fan, 2010; Li et al., 2012; Zhang et al., 2012; Zhou et al., 2012; Guo et al., 2015; Tian et al., 2018; Cheng et al., 2020. Sources of published data in Fig. 5 and Fig. 6 are the same as in this figure

全岩微量元素测试结果显示,粗面安山岩具有较高的稀土元素含量(∑REE=345×10-6~426×10-6)。在球粒陨石标准化稀土元素配分图中(图 5a),粗面安山岩具有右倾的配分曲线((La/Yb)N =21~23)与较强的负铕异常(δEu=0.63~0.65;δEu=2×EuN /(SmN + GdN)),与前人报道的腾冲粗面安山岩稀土元素配分型式一致。II类辉长岩包体稀土元素总量较低(∑REE=100×10-6~120×10-6),它们也呈现右倾的稀土配分曲线((La/Yb)N =11~16),但铕呈较弱的正异常(δEu=1.01~1.29)(图 5a)。在原始地幔标准化微量元素蛛网图中(图 5b),粗面安山岩具有较高的微量元素含量,相对于高场强元素,它们富集大离子亲石元素(例如,Rb、Ba和K等),同时具有Nb-Ta-Ti的负异常与Pb的正异常(图 5b)。II类辉长岩包体与粗面安山岩具有相似的微量元素配分曲线,即富集大离子亲石元素,亏损高场强元素,同时具有Nb-Ta-Ti的负异常与Pb的正异常(图 5b)。与粗面安山岩相比,II类辉长岩的微量元素总量较低,同时具有更显著的Nb-Ta负异常(图 5b)。

图 5 粗面安山岩与II类辉长岩包体全岩的球粒陨石标准化稀土配分图(a)和原始地幔标准化的微量元素蛛网图(b)(标准化值据Sun and McDonough, 1989) Fig. 5 Chondrite-normalized rare earth element diagrams (a) and primitive mantle-normalized trace element patterns (b) for the trachyandesitic magmas and the entrained gabbroic xenoliths (normalization values from Sun and McDonough, 1989)

同位素测试结果显示,粗面安山岩具有较高87Sr/86Sr(0.7084~0.7085)与较低的143Nd/144Nd(0.51217~0.51220)同位素比值,较高的206Pb/204Pb(18.02~18.07)、207Pb/204Pb(15.66~15.69)与208Pb/204Pb(39.11~39.20)同位素比值(表 2图 6)。与前人已发表的腾冲第三期玄武岩-粗面安山岩-安山岩成分对比,这些粗面安山岩具有较高的87Sr/86Sr与较低的143Nd/144Nd同位素比值(图 6b)。II类辉长岩包体与粗面安山岩具有较一致的Sr-Nd-Pb同位素组成,即较高的87Sr/86Sr(0.7068~0.7079)比值与较低的143Nd/144Nd(0.51217~0.51223)比值,较高的206Pb/204Pb(18.00~18.09)、207Pb/204Pb(15.65~15.67)与208Pb/204Pb(39.09~39.11)比值(表 2图 6)。

图 6 腾冲后碰撞高钾钙碱性火山岩与辉长岩包体Sr-Nd-Pb同位素相关图解 (a)辉长岩包体与粗面安山岩(143Nd/144Nd)i-(87Sr/86Sr)i相关图;(b)放大版的图a;(c)(207Pb/204Pb)i-(206Pb/204Pb)i相关图;(d)(208Pb/204Pb)i-(206Pb/204Pb)i相关图.腾冲地区花岗岩与片麻岩成分据Chen et al. (2002)杨启军等(2009);藏南上地壳与下地壳成分据Miller et al. (1999); MORB(Workman and Hart, 2005)、EMI(Hofmann, 1997)、EMII(Hofmann, 1997)和东印度大陆边缘沉积物(EIS; 潘桂棠等, 2004; Najman, 2006)成分范围以箭头或轮廓线显示 Fig. 6 Sr-Nd-Pb isotopic compositions for the post-collisional HKCA magmas and Type II gabbroic xenoliths in the Tengchong volcanic field (a) (143Nd/144Nd)i vs. (87Sr/86Sr)i of Type II gabbroic xenoliths and trachyandesitic magmas; (b) Enlargement of (143Nd/144Nd)i vs. (87Sr/86Sr)i in Fig. 6a; (c) (207Pb/204Pb)i vs. (206Pb/204Pb)i; (d) (208Pb/204Pb)i vs. (206Pb/204Pb). Data of granites and gneisses in the Tengchong terrane from Chen et al. (2002) and Yang et al. (2009). Lower (LC) and upper (UC) continental crust of southern Tibet from Miller et al. (1999). The compositions of DMM (Workman and Hart, 2005), EMI (Hofmann, 1997), EMII (Hofmann, 1997) and eastern Indian continental margin sediments (EIS) (Pan et al., 2004; Najman, 2006 and references therein) are also shown as arrows or outline fields
4.3 矿物地球化学特征 4.3.1 矿物主量元素特征

粗面安山岩斑晶的主量元素成分见表 3。斜长石斑晶以拉长石为主,具有较高的An(45 ~68)值(图 7a)。Na-Ca交换平衡计算结果显示,斜长石与熔体(全岩)成分平衡(Kd=0.14 ~0.28; Kd(An-Ab)Pl-liq=[XAbPlXAlO1.5liqXCaOliq/XAnPlXNaO0.5liqXSiO2liq]=0.27±0.11, Putirka, 2008)(图 8a)。少数他形的斜长石颗粒具有较高的An(69~73)值,它们与熔体不平衡(Kd=0.12~0.14),可能为早期基性岩浆的循环晶(图 8a表 3)。单斜辉石为普通辉石(En42-48Wo40-43Fs11-19),具有较高的Mg#(69~83,多数>78)值(图 7b表 3)。Fe-Mg交换平衡计算结果显示,多数单斜辉石(Mg#=78~83)与熔体平衡(Kd=0.26~0.38; KdCpx-liq=[XFeCpxXMgliq/XMgCpxXFeliq]=0.28±0.08, Putirka, 2008)(图 8b)。具反环带的单斜辉石自核部向边缘Mg#值升高,TiO2、Cr2O3和MnO含量也陡然升高(表 3)。斜方辉石为斜顽辉石(En72-79Wo3-4Fs18-25),具有较高的Mg#(70~82,多数>78)值(表 3)。Mg#(78~82)值较高的斜方辉石与熔体平衡(Kd=0.28~0.37; KOpx-liq d=[XFeOpxXMgliq/XMgOpxXFeliq]=0.29±0.06,Putirka, 2008)。与单斜辉石相似,发育反环带结构的斜方辉石自核向边Mg#值升高,TiO2、Cr2O3和MnO含量也陡然升高(表 3)。

表 3 粗面安山岩斑晶矿物核部主量元素组成(wt%) Table 3 Major element composition (wt%) of phenocrysts from the trachyandesitic rocks

图 7 斑晶、聚晶及辉长岩包体中斜长石(a)与辉石(b)分类定名图 Fig. 7 Denomination of plagioclase on the ternary feldspar diagram (a) and denomination of clinopyroxene and orthopyroxene based on the enstatite-ferrosilite-dedenbergite quadrilateral (b) from phenocrysts, glomerocrysts and gabbroic xenoliths

聚晶斜长石、单斜辉石与斜方辉石具有较均一的化学成分(表 4图 7)。斜长石以拉长石(An58-64)为主,与熔体Na-Ca交换平衡(Kd=0.17~0.22)(表 4图 8a)。单斜辉石为普通辉石(En45-49Wo39-42Fs11-13),具有较高的Mg#(78~82)值(表 4图 7b)。Fe-Mg交换平衡计算结果显示,这些单斜辉石与全岩(Kd=0.23~0.29)平衡(图 8b)。斜方辉石以斜顽辉石(En75-80Wo3-4Fs16-21)为主,具有较高的Mg#(78~83)值(表 4),它们与全岩Fe-Mg交换平衡(Kd=0.26~0.37)。整体而言,聚晶矿物成分与斑晶相似,均与寄主火山岩平衡。

表 4 聚晶矿物核部的主量元素组成(wt%) Table 4 Major element composition (wt%) of minerals from the glomerocrysts

图 8 斜长石An含量与Kd(a)和单斜辉石Mg#值与Kd(b)相关图解 Fig. 8 Diagrams of Kd vs. An contents of plagioclase (a) and Kd vs. Mg# of clinopyroxene (b)

辉长岩包体主要矿物的代表性成分见表 5。包体斜长石成分与斑晶斜长石较一致,其中,I类辉长岩包体斜长石为拉长石(An57-64),具有较高的SiO2(52.53%~54.45%)、CaO(11.43% ~12.71%)和Al2O3(28.92%~30.44%)含量(表 5图 7a);II类辉长岩包体斜长石也以拉长石(An54-61)为主,少量发育海绵结构的斜长石An(44~51)值较低,但FeO与MgO含量无明显的变化,可能形成于斜长石的降压分解。辉长岩包体中单斜辉石成分具有较大的差异(表 5图 7b):I类辉长岩包体中的单斜辉石(En41-46Wo41-45Fs13-14)与斑晶单斜辉石成分相似,具有较高的Mg#(77~79)值,较高的CaO(11.42%~12.69%)和FeO(7.51%~8.24%)含量(表 5图 7b);II类辉长岩包体单斜辉石也为普通辉石(En40-47Wo41-46Fs14-17),但Mg#(71~76)值较低,CaO(19.62%~21.68%)和FeO(8.22%~10.29%)含量较高(表 5图 7b)。斜方辉石主要发育于I类辉长岩包体中,成分与斑晶斜方辉石较一致,以斜顽辉石(En74-79Wo3-4Fs13-14)为主并具有较高的Mg#(76~82)值(表 5图 7b)。II类辉长岩包体中仅发育少量微晶斜方辉石,它们具有较低的Mg#值(64~70)和较高的FeO(18.42%~22.01%)含量(表 5),其成分与粗面安山岩基质中的斜方辉石微晶(Mg#=67~72; FeO=16.56%~28.03%)较一致,可能为寄主岩浆晚期快速冷却结晶的产物。矿物与熔体交换平衡计算结果表明:I类辉长岩包体中的斜长石、单斜辉石和斜方辉石均与寄主岩浆(Mg#=0.61)平衡(图 8);II类辉长岩斜长石An值较低,单斜辉石的Mg#值也较低,可能结晶自演化程度更高的岩浆。以该地区演化程度较高的粗面安山岩全岩(TC1311: Mg#=0.56)代表熔体,II类辉长岩包体中的斜长石与单斜辉石的成分与之平衡(图 8)。

表 5 辉长岩包体矿物核部的主量元素组成(wt%) Table 5 Major element composition (wt%) of minerals from the gabbroic xenoliths
4.3.2 单斜辉石微量元素特征

单斜辉石原位微区微量元素测试结果表明,不同岩石类型(包括I类辉长岩、II类辉长岩、聚晶和斑晶)中单斜辉石具有相似的微量元素分布型式(表 6图 9)。在球粒陨石标准化稀土元素配分图中(图 9a),单斜辉石均具有右倾的稀土配分型式,即相对富集轻稀土元素、亏损重稀土元素,同时具有微弱的负铕异常。其中,I类辉长岩包体、聚晶与斑晶单斜辉石具有相似稀土配分型式;II类辉长岩包体的单斜辉石具有更高的轻稀土元素含量和更强烈的负铕异常(图 9a)。在原始地幔标准化微量元素蛛网图中(图 9b),单斜辉石具有相似的微量元素配分型式:相对于原始地幔,它们更富集微量元素,同时具有Th、Sr和Ti的负异常;I类辉长岩包体、聚晶与斑晶单斜辉石的微量元素配分型式较相似,II类辉长岩单斜辉石具有更高的微量元素含量(图 9b)。

表 6 斑晶、矿物聚晶与辉长岩包体中单斜辉石原位微区微量元素(×10-6)组成 Table 6 Trace element compositions (×10-6) of clinopyroxene from phenocrysts, glomerocrysts and gabbroic xenoliths

图 9 斑晶、聚晶与辉长岩包体单斜辉石的球粒陨石标准化稀土配分图(a)及原始地幔标准化微量元素蛛网图(b)(标准化值据Sun and McDonough, 1989) Fig. 9 Chondrite-normalized rare earth element diagrams (a) and primitive mantle-normalized trace element patterns (b) for clinopyroxene from phenocrysts, glomerocrysts and gabbroic xenoliths (normalization values from Sun and McDonough, 1989)
5 讨论 5.1 矿物聚晶的来源

腾冲后碰撞高钾钙碱性岩浆作用形成了近连续的玄武岩-玄武安山岩-粗面安山岩-英安岩系列(图 4)。岩石学、矿物学、矿物温压计算与EC-AFC模拟结果显示,粗面安山岩为玄武质岩浆在中下地壳岩浆房经历地壳混染与分离结晶(AFC)的产物(Cheng et al., 2020)。近年来,大量的岩石学、地球化学、数值模拟和地球物理学研究结果均表明,地壳岩浆房系统由相互连接的以熔体为主的区域和以晶体为主的晶粥区域构成(Nakagawa et al., 1999; Cashman et al., 2017; Jackson et al., 2018; Luo et al., 2018; Bennett et al., 2019; Edmonds et al., 2019; Caracciolo et al., 2020)。由于温度较低、晶体含量较高和粘度较高,晶粥不能直接喷出地表,但可以被快速上升的岩浆分解并以单个晶体、聚晶(glomerocryst)、堆晶团块(cumulate nodules)或包体的方式带至地表(Nakagawa et al., 1999; Cashman et al., 2017; Bennett et al., 2019)。

本研究报道的粗面安山岩斑晶含量较高,同时包含大量的矿物聚晶和辉长岩包体(图 2图 3)。这些矿物聚晶具有如下特点:(1)具有自形-半自形中粗粒结构,晶形完整,无折断或变形迹象;(2)与粗面安山岩斑晶矿物组成相同,主要由斜长石、单斜辉石与斜方辉石构成;(3)聚晶矿物之间(尤其是斜长石)呈高角度接触的开放结构;(4)矿物间充填大量的基质玻璃或气孔。火山岩中类似的聚晶结构通常是岩浆捕获的岩浆房的半固结部分,代表岩浆喷发前形成的晶粥(Crystal mush)(Nakagawa et al., 1999; Holness et al., 2005; Cashman et al., 2017; Edmond et al., 2019; Lissenberg et al., 2019)。例如,东太平洋海隆(Pan and Batiza, 2003; Moore et al., 2014)与北冰洋加克利(Bennett et al., 2019)火山岩中的矿物聚晶均为熔体捕获的岩浆房内部晶粥的碎片。这些聚晶矿物的成分可能与寄主熔体平衡或不平衡,分别代表寄主熔体或早期熔体形成的晶粥(Pan and Baitza, 2003; Moore et al., 2013, 2014)。本研究中,聚晶矿物(包括拉长石、普通辉石和斜顽辉石)与寄主粗面安山岩斑晶具有一致的主量元素组成(图 7图 8),它们均与粗面安山岩的全岩成分平衡(图 8)。同时,聚晶单斜辉石与斑晶具有相似的稀土和微量元素组成(图 9)。矿物组成与矿物化学成分上的相似性暗示这些聚晶矿物为寄主粗面安山质岩浆结晶作用形成的晶粥。

5.2 辉长岩包体的成因

辉长岩包体样品新鲜,与寄主火山岩接触清晰,无熔蚀、烘烤、冷凝等反应结构(图 2)。包体具有自形-半自形粒状结构,主要由拉长石、普通辉石与斜顽辉石组成,矿物晶形完整且未发生明显的变形(图 3)。前人在腾冲晚更新世粗面安山岩中发现了一些长英质麻粒岩和二辉质麻粒岩包体,它们具有粒状变晶结构并以包含紫苏辉石为特征(林木森等,2014)。本研究发现的辉长岩包体的结构、矿物组成和化学成分均与前人报道的麻粒岩包体具有较大的差异,因此,排除这些辉长岩包体为岩浆上升过程中捕获的地壳组分。

I类辉长岩包体的结构、矿物组成特征和化学成分均与矿物聚晶具有较大的相似性(图 3)。例如,I类辉长岩包体也主要由自形-半自形的拉长石、普通辉石和斜顽辉石组成,并且矿物间呈高角度接触,间隙还充填少量的基质玻璃与气孔(图 3);I类辉长岩包体中的斜长石、单斜辉石及斜方辉石与聚晶矿物具有一致的主量元素组成,它们均与粗面安山岩达到了成分上的平衡(图 8);尤为重要的是,I类辉长岩包体的单斜辉石与聚晶单斜辉石呈现较一致的稀土和微量元素分配型式(图 5表 5)。基于以上岩石结构、矿物组成与矿物化学成分上的相似性,我们认为I类辉长岩包体与矿物聚晶相同,为岩浆喷发前形成晶粥的“碎片”。与矿物聚晶相比,I类辉长岩包体粒径较大,矿物间紧密程度较高,玻璃基质含量较低,它们可能代表岩浆房晶体含量更高的晶粥。

与I类辉长岩包体和矿物聚晶相比,II类辉长岩包体在岩石结构与矿物组成上均具有差异,主要表现为II类辉长岩包体的矿物间接触更为紧密,矿物间仅发现少量的玻璃基质及斜方辉石微晶(图 3)。为了更好的限定II类辉长岩包体的成因,我们对比研究了辉长岩包体与寄主粗面安山岩全岩的主微量元素与同位素地球化学特征(图 4图 5图 6)。在全岩主量元素(SiO2、CaO、FeO和K2O)和微量元素(Sr和Ba)与MgO含量的相关图解中(图 4),II类辉长岩包体的成分均落在腾冲玄武岩-玄武安山岩-粗面安山岩-英安岩的演化序列中。全岩微量元素测试结果显示,II类辉长岩包体与寄主粗面安山岩具有相似的稀土元素与微量元素配分型式(图 5)。尤为重要的是,辉长岩包体与寄主粗面安山岩具有一致的Sr-Nd-Pb同位素组成(图 6)。这些地球化学特征表明II类辉长岩包体与寄主粗面安山质岩浆是同源的,它们均来自腾冲地壳安山质岩浆房。II类辉长岩包体与斑晶单斜辉石矿物相似的微量元素组成进一步证实了这些辉长岩包体与寄主粗面安山岩的同源性(图 9)。手标本与显微镜下观察结果均显示辉长岩包体具有粒状结构,未见明显的矿物分带与定向分布;同时,单斜辉石Mg#值较低, 斜长石An含量较低,全岩成分也未表现出异常高的CaO含量。因而我们推断,II类辉长岩包体未经历显著的堆晶过程,是岩浆房结晶作用的产物。II类辉长岩包体与矿物聚晶和I类辉长岩包体同时被粗面安山质岩浆携带至地表,并且II辉长岩包体的矿物间也发现了少量的基质玻璃(图 2),我们认为II类辉长岩包体的来源与矿物聚晶及I类辉长岩包体相似,也为岩浆房中晶体含量较高的晶粥。与I类辉长岩相比,II类辉长岩包体矿物间紧密程度较高,可能代表岩浆房较早期形成的经过部分压实的晶粥(Bennett et al., 2019)。II类辉长岩包体单斜辉石Mg#值较低,斜长石An含量也较低(图 8),指示其结晶于演化程度较高的岩浆(Stamper et al., 2014),这也解释了II类辉长岩包体单斜辉石更高的微量元素含量与更强烈的轻重稀土元素分异(图 5)。

5.3 矿物平衡温度与压力估算 5.3.1 单斜辉石-熔体温压计

大量的实验岩石学研究结果表明,单斜辉石硬玉成分(NaAlSi2O6)的含量对压力较敏感,因而可以较好的限定单斜辉石结晶的压力条件(Putirka et al., 2003; Armienti et al., 2007; Putirka, 2008; Mollo et al., 2013; Giacomoni et al., 2016; Neave and Putirka, 2017; Ubide and Kamber, 2018)。单斜辉石在粗面安山岩与辉长岩包体中均较好的发育(图 3),为我们限定其形成的温度和压力条件提供了良好的基础。Neave and Putirka (2017)提出的单斜辉石-熔体温压计很好地提高了压力估算的精度(±1.4kbar, 1σ),本研究利用该方法对单斜辉石结晶条件进行了估算。温度与压力估算时首先需判断单斜辉石与熔体是否平衡,其次还需对熔体的水含量进行估算。早期多利用单斜辉石与熔体间Fe-Mg交换系数(Kd=0.28±0.08)判断二者是否平衡,近年来大量的研究结果表明,该交换系数并非常数,而是与熔体的成分、温度和氧逸度等条件有关(Putirka, 2008; Neave and Putirka, 2017)。Mollo et al. (2013)提出根据回归分析的方法获得与熔体平衡的单斜辉石成分,再将该成分与实际测出的单斜辉石成分进行对比,可以更好地判断单斜辉石与熔体是否平衡。该方法近年来得到了较广泛的应用(Haddadi et al., 2017; Neave and Putirka, 2017; Eskandari et al., 2018; Ubide et al., 2019; Di et al., 2020)。因此,在Fe-Mg交换平衡计算的基础上,我们进一步利用该方法验证单斜辉石与熔体的平衡(图 8b图 10a)。虽然II类辉长岩包体与寄主火山岩成分不平衡,但以研究区Mg#值较低的粗面安山岩(TC1311)全岩代表熔体,II类辉长岩包体中的单斜辉石成分与之达到平衡(图 8b),因此我们以该全岩成分代表II类辉长岩包体的熔体。根据Waters and Lange (2015)提出的斜长石-熔体湿度计,初步估算粗面安山质岩浆的水含量为2.5~3.0±0.4%。

图 10 单斜辉石-熔体间平衡判断(a,据Mollo et al., 2013)以及聚晶、辉长岩包体与粗面安山岩单斜辉石-熔体(Neave and Putirka, 2017)和单斜辉石-斜方辉石(Putirka, 2008)温压计算结果(b) Fig. 10 Equilibrium testing of clinopyroxene based on a comparison of observed and predicted diopside-hedenbergite (DiHd) components (a, based on Mollo et al., 2013) and P-T calculation results of the clinopyroxene-liquid thermobarometer (Neave and Putirka, 2017) and the clinopyroxene-orthopyroxene (Putirka, 2008) for glomerocrysts, gabbroic xenoliths and trachyandesites (b)

选取与熔体成分平衡的单斜辉石进行温度与压力计算。结果显示,单斜辉石斑晶的结晶温度为1088~1137±25℃,压力为5.3~9.2±1.4kbar(深度20~31km)(图 10b),与前人单斜辉石斑晶的温压计算结果较一致(Hu et al., 2018; Cheng et al., 2020)。聚晶单斜辉石的结晶温度(1083~1106±25℃)与压力(4.1~6.3±1.4kbar;17~23km)略低于单斜辉石斑晶(图 10b)。I类辉长岩包体单斜辉石的结晶温度(1098~1110±25℃)与聚晶单斜辉石较一致,但结晶压力(5.4~6.7±1.4kbar;21~25km)略高(图 10b)。与I类辉长岩包体及矿物聚晶对比,II类辉长岩包体的结晶压力(4.8~6.9±1.4kbar;深度18~25km)与前者一致,但结晶温度(1046~1082±25℃)略低(图 10b)。整体而言,两类辉长岩包体与聚晶矿物的结晶温度与压力条件均低于火山岩斑晶(图 10b)。

5.3.2 单斜辉石-斜方辉石温压计

单斜辉石-斜方辉石温压计也可以较好的估算矿物结晶的温度和压力条件(Brey and Khler, 1990; Putirka, 2008)。显微镜下观察结果表明,I类辉长岩包体与矿物聚晶的单斜辉石与斜方辉石较自形,无环带或变形现象,二者在结构上平衡。电子探针分析结果显示,同一单斜辉石或斜方辉石颗粒的不同部位成分差异很小。Fe-Mg交换平衡计算结果(KdFe-Mg =0.94~1.22; 1.09±0.16, Putirka, 2008)也证实单斜辉石与斜方辉石成分上平衡。因此,单斜辉石-斜方辉石温压计可以为I类辉长岩包体和矿物聚晶提供较可靠的温压条件的限定。计算结果显示,聚晶矿物结晶的温度(997~1041±50℃)与压力(3.0~5.1±2.8kbar)条件略低于I类辉长包体(1024~1045±50℃; 4.0~4.4±2.8kbar)。利用单斜辉石-斜方辉石矿物对获得的结晶压力与单斜辉石-熔体压力结果较一致,但温度略低(图 9b)。考虑到二辉石温度计较大的误差范围(±50℃;Putirka, 2008),我们认为两种方法计算的温度结果基本一致。

5.4 岩浆房系统与深部岩浆过程

腾冲新生代火山区上新世以来发育近连续的玄武岩-玄武安山岩-粗面安山质-英安岩系列(图 4)。岩石学、矿物学、温压估算与EC-AFC模拟结果均显示,腾冲上地壳(9~14km)与中下地壳(13~28km)均存在岩浆房(Cheng et al., 2020),这与该地区大地电磁低阻异常(12~30km;姜枚等,2012谭捍东等,2013)与地震S波低速异常(9~16km;Xu et al., 2018b)特征指示的岩浆房分布相吻合。腾冲粗面安山质岩浆为玄武质岩浆在中下地壳岩浆房(13~28km)经历分离结晶与地壳混染(AFC)过程的产物(Cheng et al., 2020)。本研究中,寄主粗面安山岩斑晶的压力计算结果(5.3~9.2±1.4kbar;深度20~31km)也显示单斜辉石斑晶来自腾冲深部地壳岩浆房(图 10b图 11)。这些粗面安山岩包含大量的辉长岩包体与矿物聚晶,岩石学、矿物学、矿物化学和全岩地球化学特征均表明,它们均来自粗面安山质岩浆房(图 3-图 6图 10)。矿物聚晶与辉长岩包体均发育粒状结构,矿物间呈高角度接触并贯入基质玻璃,呈现典型的“开放”的结构特征(图 3)。与东太平洋海隆(Pan and Batiza, 2003; Moore et al., 2014)、北冰洋加克利(Bennett et al., 2019)和新西兰鲁阿佩胡(Nakagawa et al., 1999, 2002)火山岩中包含的矿物聚晶与辉长岩包体相似,它们可能代表岩浆房中上部以晶体为主的晶粥(图 11)。矿物聚晶和I类辉长岩包体与粗面安山岩的斑晶呈现较一致的矿物组成和化学成分,矿物间贯入较高的基质玻璃,指示它们可能代表寄主粗面安山质岩浆形成的晶粥,这些晶粥未经历明显的压实作用。II类辉长岩包体矿物具有演化的成分(单斜辉石较低的Mg#值和斜长石较低的An含量)和具有较低的结晶温度,它们可能代表演化程度更高的粗面安山质岩浆形成的晶粥,可能形成于岩浆房系统更靠近围岩的部位(图 11)。II类辉长岩包体更紧密的岩石结构指示它们可能经历了部分压实作用(图 11)。矿物聚晶与辉长岩包体的形成深度(17~25km)均处于腾冲地壳深部岩浆房的中上部(图 11)。在岩浆房中上部,斜长石等密度较小的矿物易于聚集,这与矿物聚晶和辉长岩包体中斜长石含量较高一致(图 3)。同时,岩浆房中上部挥发分或气相组分的存在也利于矿物聚晶的形成(Cashman et al., 2017),这也与粗面安山岩大量发育的气孔构造一致(图 2d图 3g)。粗面安山质岩浆快速上升并喷出地表的过程中将位于岩浆房中上部的晶粥破碎,并以辉长岩包体与聚晶的形式携带至地表(Nakagawa et al., 1999, 2002; Yang et al., 2012; Price et al., 2017; Cashman et al., 2017; Maclennan, 2019; Edmonds et al., 2019; Caracciolo et al., 2020),这也解释了腾冲晚更新世粗面安山岩中同时包含了大量的矿物聚晶与辉长岩包体(图 2c-f)。腾冲粗面安山岩中单斜辉石与斜方辉石斑晶的反环带结构记录了岩浆房中基性熔体的补给过程(Hu et al., 2018; Cheng et al., 2020),可能是粗面安山质岩浆喷发的主要诱导因素。

图 11 腾冲粗面安山岩岩浆房系统示意图 ① II类辉长岩包体代表粗面安山质岩浆形成的较紧密的“晶粥”;② I类辉长岩包体代表寄主粗面安山质岩浆形成的“晶粥”;③矿物聚晶代表寄主粗面安山质岩浆形成的“晶粥”碎片 Fig. 11 Schematic diagram to illustrate the plumbing system of trachyandesitic magmas in the Tengchong volcanic field Numbers ①, ② and ③ denote the distribution of Type II gabbroic xenoliths, Type I gabbroic xenoliths and glomerocrysts, respectively
6 结论

(1) 腾冲新生代火山区晚更新世粗面安山岩中包含大量的辉长岩包体和矿物聚晶。辉长岩包体、聚晶与火山岩斑晶的岩石学、矿物学特征,尤其是单斜辉石的微量元素与全岩Sr-Nd-Pb同位素组成的对比研究结果表明,矿物聚晶和辉长岩包体与寄主粗面安山质岩浆同源,均来自腾冲深部地壳粗面安山质岩浆房。

(2) 结合辉长岩包体与聚晶的结构特征与矿物温压计算结果,认为聚晶与辉长岩包体均代表粗面安山质岩浆房中上部的晶粥,快速上升的安山质岩浆将这些晶粥破碎并以聚晶和辉长岩包体的形式带至地表,基性岩浆的补给是导致粗面安山质岩浆喷发的重要诱发机制。

致谢      中国地震局地质研究所樊祺诚研究员、中山大学地质科学与工程学院刘志超副教授与中国地质大学(北京)黄丰副教授给予了有益的建议;K-Ar年龄由中国地震局地质研究所李大明研究员完成;中国科学院地质与地球物理研究所高金亮与丁磊磊博士在野外考察与采样工作中给予了大量帮助;中国科学院地质与地球物理研究所张丽红博士在全岩地球化学成分测试中给予了帮助;匿名审稿专家提出宝贵的修改意见;在此一并致以衷心的感谢!

参考文献
Annen C, Blundy JD and Sparks RSJ. 2006. The genesis of intermediate and silicic magmas in deep crustal hot zones. Journal of Petrology, 47(3): 505-539 DOI:10.1093/petrology/egi084
Arculus RJ and Wills KJA. 1980. The petrology of plutonic blocks and inclusions from the Lesser Antilles island arc. Journal of Petrology, 21(4): 743-799 DOI:10.1093/petrology/21.4.743
Armienti P, Tonarini S, Innocenti F and D'Orazio M. 2007. Mount Etna pyroxene as tracer of petrogenetic processes and dynamics of the feeding system. In: Beccaluva L, Bianchini G and Wilson M (eds.). Cenozoic Volcanism in the Mediterranean Area. Geological Society of America Special Paper, 418: 265-276
Aslan Z, Arslan M, Temizel İ and Kaygusuz A. 2014. K-Ar dating, whole-rock and Sr-Nd isotope geochemistry of calc-alkaline volcanic rocks around the Gümüşhane area:Implications for post-collisional volcanism in the Eastern Pontides, Northeast Turkey. Mineralogy and Petrology, 108(2): 245-267
Bai DH, Meju MA and Liao ZJ. 2001. Magnetotelluric images of deep crustal structure of the Rehai geothermal field near Tengchong, southern China. Geophysical Journal International, 147(3): 677-687 DOI:10.1046/j.0956-540x.2001.01568.x
Bennett EN, Lissenberg CJ and Cashman KV. 2019. The significance of plagioclase textures in mid-ocean ridge basalt (Gakkel Ridge, Arctic Ocean). Contributions to Mineralogy and Petrology, 174(6): 49 DOI:10.1007/s00410-019-1587-1
Brey GP and Köhler T. 1990. Geothermobarometry in four-phase lherzolites:Ⅱ. New thermobarometers, and practical assessment of existing thermobarometers. Journal of Petrology, 31(6): 1353-1378
Camejo-Harry M, Melekhova E, Blundy J, Attridge W, Robertson R and Christopher T. 2018. Magma evolution beneath Bequia, Lesser Antilles, deduced from petrology of lavas and plutonic xenoliths. Contributions to Mineralogy and Petrology, 173(10): 77 DOI:10.1007/s00410-018-1504-z
Caracciolo A, Bali E, Guefinnsson GH, Kahl M, Halldórsson SA, Hartley ME and Gunnarsson H. 2020. Temporal evolution of magma and crystal mush storage conditions in the Bárðarbunga-Veiðivötn volcanic system, Iceland. Lithos, 352-353: 105234 DOI:10.1016/j.lithos.2019.105234
Cashman KV and Sparks RSJ. 2013. How volcanoes work:A 25 year perspective. Geological Society of America Bulletin, 125(5-6): 664-690 DOI:10.1130/B30720.1
Cashman KV, Sparks RSJ and Blundy JD. 2017. Vertically extensive and unstable magmatic systems:A unified view of igneous processes. Science, 355(6331): eaag3055 DOI:10.1126/science.aag3055
Chen F, Satir, M, Ji J and Zhong D. 2002. Nd-Sr-Pb isotopes of Tengchong Cenozoic volcanic rocks from western Yunnan, China:Evidence for an enriched-mantle source. Journal of Asian Earth Sciences, 21(1): 39-45
Chen TF. 2003. The petrology of the volcanic rocks in Tengchong, Yunnan. Sedimentary Geology and Tethyan Geology, 23(4): 56-61 (in Chinese with English abstract)
Cheng ZH, Guo ZF, Zhang ML and Zhang LH. 2012. CO2 flux estimations of hot springs in the Tengchong Cenozoic volcanic field, Yunnan Province, SW China. Acta Petrologica Sinica, 28(4): 1217-1224 (in Chinese with English abstract)
Cheng ZH, Guo ZF, Zhang ML and Zhang LH. 2014. Carbon dioxide emissions from Tengchong Cenozoic volcanic field, Yunnan Province, SW China. Acta Petrologica Sinica, 30(12): 3657-3670 (in Chinese with English abstract)
Cheng ZH, Guo ZF, Dingwell DB, Li XH, Zhang ML, Liu JQ, Zhao WB and Lei M. 2020. Geochemistry and petrogenesis of the post-collisional high-K calc-alkaline magmatic rocks in Tengchong, SE Tibet. Journal of Asian Earth Sciences, 193: 104309 DOI:10.1016/j.jseaes.2020.104309
Cong BL, Chen QY, Zhang RY, Wu GY and Xu P. 1994. Petrogenesis of Cenozoic volcanic rocks in Tengchong region of western Yunnan Province, China. Science in China (Series B), 24(4): 441-448 (in Chinese)
Cooper GF, Blundy JD, Macpherson CG, Humphreys MCS and Davidson JP. 2019. Evidence from plutonic xenoliths for magma differentiation, mixing and storage in a volatile-rich crystal mush beneath St. Eustatius, Lesser Antilles. Contributions to Mineralogy and Petrology, 174(5): 39 DOI:10.1007/s00410-019-1576-4
Costa F, Dungan MA and Singer BS. 2002. Hornblende- and phlogopite-bearing gabbroic xenoliths from Volcán San Pedro (36°S), Chilean Andes:Evidence for melt and fluid migration and reactions in subduction-related plutons. Journal of Petrology, 43(2): 219-241 DOI:10.1093/petrology/43.2.219
Davidson JP, Morgan DJ, Charlier BLA, Harlou R and Hora JM. 2007. Microsampling and isotopic analysis of igneous rocks:Implications for the study of magmatic systems. Annual Review of Earth and Planetary Sciences, 35(1): 273-311 DOI:10.1146/annurev.earth.35.031306.140211
Di YK, Tian W, Chen MM, Li ZF, Chu ZY and Liang J. 2020. A method to estimate the pre-eruptive water content of basalts:Application to the Wudalianchi-Erkeshan-Keluo volcanic field, northeastern China. American Mineralogist, 105(2): 149-161
Duan XZ, Fan HR, Zhang HF, Yaxley G, Santosh M, Tian HC, Tan KX, Tang ZP, Xie YS, Xiao YL, Hou ZH and Guo HF. 2019. Melt inclusions in phenocrysts track enriched upper mantle source for Cenozoic Tengchong volcanic field, Yunnan Province, SW China. Lithos, 324-325: 180-201 DOI:10.1016/j.lithos.2018.10.034
Edmonds M, Cashman KV, Holness M and Jackson M. 2019. Architecture and dynamics of magma reservoirs. Philosophical Transactions of the Royal Society A:Mathematical, Physical and Engineering Sciences, 377(2139): 20180298 DOI:10.1098/rsta.2018.0298
Eskandari A, Amini S, De Rosa R and Donato P. 2018. Nature of the magma storage system beneath the Damavand volcano (N. Iran):An integrated study. Lithos, 300-301: 154-176 DOI:10.1016/j.lithos.2017.12.002
Fan QC, Liu RX, Wei HQ, Shi LB and Sui JL. 1999. The magmatic evolution of the active volcano in the Tengchong area. Geological Review, 45(Suppl.1): 895-904 (in Chinese with English abstract)
Fodor RV, Rudek EA and Bauer GR. 1993. Hawaiian magma-reservoir processes as inferred from the petrology of gabbro xenoliths in basalt, Kahoolawe Island. Bulletin of Volcanology, 55(3): 204-218 DOI:10.1007/BF00301517
Fodor RV and Moore RB. 1994. Petrology of gabbroic xenoliths in 1960 Kilauea basalt:Crystalline remnants of prior (1955) magmatism. Bulletin of Volcanology, 56(1): 62-74 DOI:10.1007/BF00279729
Fodor RV and Johnson KG. 2016. Origin of Miocene andesite and dacite in the Goldfield-Superstition volcanic province, central Arizona:Hybrids of mafic and silicic magma mixing. Geochimica et Cosmochimica Acta, 185: 394-417 DOI:10.1016/j.gca.2016.04.001
Gao JF, Zhou MF, Robinson PT, Wang CY, Zhao JH and Malpas J. 2015. Magma mixing recorded by Sr isotopes of plagioclase from dacites of the Quaternary Tengchong volcanic field, SE Tibetan Plateau. Journal of Asian Earth Sciences, 98: 1-17 DOI:10.1016/j.jseaes.2014.10.036
Giacomoni PP, Coltorti M, Bryce JG, Fahnestock MF and Guitreau M. 2016. Mt. Etna plumbing system revealed by combined textural, compositional, and thermobarometric studies in clinopyroxenes. Contributions to Mineralogy and Petrology, 171(4): 34
Gill JB. 1981. Orogenic Andesites and Plate Tectonics. Berlin: Springer, 1-390
Grove TL and Donnelly-Nolan JM. 1986. The evolution of young silicic lavas at Medicine Lake Volcano, California:Implications for the origin of compositional gaps in calc-alkaline series lavas. Contributions to Mineralogy and Petrology, 92(3): 281-302 DOI:10.1007/BF00572157
Grove TL, Elkins-Tanton LT, Parman SW, Chatterjee N, Müntener O and Gaetani GA. 2003. Fractional crystallization and mantle-melting controls on calc-alkaline differentiation trends. Contributions to Mineralogy and Petrology, 145(5): 515-533 DOI:10.1007/s00410-003-0448-z
Guo ZF, Wilson M, Liu JQ and Mao Q. 2006. Post-collisional, potassic and ultrapotassic magmatism of the northern Tibetan Plateau:Constraints on characteristics of the mantle source, geodynamic setting and uplift mechanisms. Journal of Petrology, 47(6): 1177-1220 DOI:10.1093/petrology/egl007
Guo ZF, Cheng ZH, Zhang ML, Zhang LH, Li XH and Liu JQ. 2015. Post-collisional high-K calc-alkaline volcanism in Tengchong volcanic field, SE Tibet:Constraints on Indian eastward subduction and slab detachment. Journal of the Geological Society, 172(5): 624-640 DOI:10.1144/jgs2014-078
Haddadi B, Sigmarsson O and Larsen G. 2017. Magma storage beneath Grímsvötn volcano, Iceland, constrained by clinopyroxene-melt thermobarometry and volatiles in melt inclusions and groundmass glass. Journal of Geophysical Research:Solid Earth, 122(9): 6984-6997 DOI:10.1002/2017JB014067
Hofmann AW. 1997. Mantle geochemistry:The message from oceanic volcanism. Nature, 385(6613): 219-229 DOI:10.1038/385219a0
Holness MB, Cheadle MJ and McKenzie D. 2005. On the use of changes in dihedral angle to decode late-stage textural evolution in cumulates. Journal of Petrology, 46(8): 1565-1583 DOI:10.1093/petrology/egi026
Hu JH, Song XY, He HL, Zheng WQ, Yu SY, Chen LM and Lai CK. 2018. Constraints of texture and composition of clinopyroxene phenocrysts of Holocene volcanic rocks on a magmatic plumbing system beneath Tengchong, SW China. Journal of Asian Earth Sciences, 154: 342-353 DOI:10.1016/j.jseaes.2017.12.029
Hu X, Garzanti E, Moore T and Raffi I. 2015. Direct stratigraphic dating of India-Asia collision onset at the Selandian (Middle Paleocene, 59±1Ma). Geology, 43(10): 859-862 DOI:10.1130/G36872.1
Hua YJ, Zhang SX, Li MK, Wu TF, Zou CY and Liu L. 2019. Magma system beneath Tengchong volcanic zone inferred from local earthquake seismic tomography. Journal of Volcanology and Geothermal Research, 377: 1-16 DOI:10.1016/j.jvolgeores.2019.04.002
Huangfu G and Jiang CS. 2000. Study on Tengchong Volcanic Activity. Kunming: Yunnan Scientific and Technology Press, 1-418 (in Chinese)
Humphreys MCS, Blundy JD and Sparks RSJ. 2006. Magma evolution and open-system processes at Shiveluch volcano:Insights from phenocryst zoning. Journal of Petrology, 47(12): 2303-2334 DOI:10.1093/petrology/egl045
Jackson MD, Blundy J and Sparks RSJ. 2018. Chemical differentiation, cold storage and remobilization of magma in the Earth's crust. Nature, 564(7736): 405-409 DOI:10.1038/s41586-018-0746-2
Jerram DA and Martin VM. 2008. Understanding crystal populations and their significance through the magma plumbing system. Geological Society, London, Special Publications, 304(1): 133-148 DOI:10.1144/SP304.7
Jiang CS. 1998. Period division of volcano activities in the Cenozoic era of Tengchong. Journal of Seismological Research, 21(4): 320-329 (in Chinese with English abstract)
Jiang M, Tan HD, Zhang JW, Peng M, Li QQ, Zhang LS, Xu LH and Wang W. 2012. Geophysical mode of Mazhan-Gudong magma chamber in Tengchong volcano-tectonic area. Acta Geoscientica Sinica, 33(5): 731-739 (in Chinese with English abstract)
Kelemen PB and Behn MD. 2016. Formation of lower continental crust by relamination of buoyant arc lavas and plutons. Nature Geoscience, 9(3): 197-205 DOI:10.1038/ngeo2662
Kent AJR, Darr C, Koleszar AM, Salisbury MJ and Cooper KM. 2010. Preferential eruption of andesitic magmas through recharge filtering. Nature Geoscience, 3(9): 631-636 DOI:10.1038/ngeo924
Laumonier M, Scaillet B, Pichavant M, Champallier R, Andujar J and Arbaret L. 2014. On the conditions of magma mixing and its bearing on andesite production in the crust. Nature Communications, 5: 5607 DOI:10.1038/ncomms6607
Li DM, Li Q and Chen WJ. 2000. Volcanic activities in the Tengchong volcano area since Pliocene. Acta Petrologica Sinica, 16(3): 362-370 (in Chinese with English abstract)
Li DP, Luo ZH, Liu JQ, Chen YL and Jin Y. 2012. Magma origin and evolution of Tengchong Cenozoic volcanic rocks from West Yunnan, China:Evidence from whole rock geochemistry and Nd-Sr-Pb isotopes. Acta Geologica Sinica, 86(4): 867-878 DOI:10.1111/j.1755-6724.2012.00712.x
Li N, Zhao YW, Zhang LY and Wang JL. 2020. The Quaternary eruptive sequence of the Tengchong volcanic group, southwestern China. Lithos, 354-355: 105173 DOI:10.1016/j.lithos.2019.105173
Lin MS, Peng SB, Qiao WT, Li CN and Xu HJ. 2014. Characteristics and genetic significance of high temperature granulite xenoliths in Cenozoic volcanic rocks, Tengchong, western Yunnan Province, China. Earth Science (Journal of China University of Geosciences), 39(7): 807-819 (in Chinese with English abstract) DOI:10.3799/dqkx.2014.076
Lissenberg CJ, Macleod CJ and Bennett EN. 2019. Consequences of a crystal mush-dominated magma plumbing system:A mid-ocean ridge perspective. Philosophical Transactions of the Royal Society A:Mathematical, Physical and Engineering Sciences, 377(2139): 20180014 DOI:10.1098/rsta.2018.0014
Liu Y, Shi ZT and Zhang WX. 2018. Age of a Late Quaternary (13.6 cal kyr BP) eruption and associated marker tephra of the Tengchong volcanoes, Yunnan Province, southwestern China. Palaeogeography, Palaeoclimatology, Palaeoecology, 511: 588-596 DOI:10.1016/j.palaeo.2018.09.024
Luo BJ, Zhang HF, Xu WC, Yang H, Zhao JH, Guo L, Zhang LQ, Tao L, Pan FB and Gao Z. 2018. The magmatic plumbing system for Mesozoic high-Mg andesites, garnet-bearing dacites and porphyries, rhyolites and leucogranites from West Qinling, Central China. Journal of Petrology, 59(3): 447-482 DOI:10.1093/petrology/egy035
Luo ZH, Yang ZF, Dai G, Cheng LL and Zhou JL. 2013. Crystal populations of igneous rocks and their implications in genetic mineralogy. Geology in China, 40(1): 176-181 (in Chinese with English abstract)
Maclennan J. 2019. Mafic tiers and transient mushes:Evidence from Iceland. Philosophical Transactions of the Royal Society A:Mathematical, Physical and Engineering Sciences, 377(2139): 20180021 DOI:10.1098/rsta.2018.0021
Miller C, Schuster R, Klötzli U, Frank W and Purtscheller F. 1999. Post-collisional potassic and ultrapotassic magmatism in SW Tibet:Geochemical and Sr-Nd-Pb-O isotopic constraints for mantle source characteristics and petrogenesis. Journal of Petrology, 40(9): 1399-1424 DOI:10.1093/petroj/40.9.1399
Mo XX and Pan GT. 2006. From the Tethys to the formation of the Qinghai-Tibet Plateau:Constrained by tectono-magmatic events. Earth Science Frontiers, 13(6): 43-51 (in Chinese with English abstract)
Mo XX, Zhao ZD, Deng JF, Flower M, Yu XH, Luo ZH, Li YG, Zhou S, Dong GC, Zhu DC and Wang LL. 2006. Petrology and geochemistry of postcollisional volcanic rocks from the Tibetan Plateau: Implications for lithosphere heterogeneity and collision-induced asthenospheric mantle flow. In: Postcollisional Tectonics and Magmatism in the Mediterranean Region and Asia. Geological Society of America Special Papers, 409: 507-530
Mollo S, Putirka K, Misiti V, Soligo M and Scarlato P. 2013. A new test for equilibrium based on clinopyroxene-melt pairs:Clues on the solidification temperatures of Etnean alkaline melts at post-eruptive conditions. Chemical Geology, 352: 92-100 DOI:10.1016/j.chemgeo.2013.05.026
Moore A, Coogan LA, Costa F and Perfit MR. 2014. Primitive melt replenishment and crystal-mush disaggregation in the weeks preceding the 2005~2006 eruption 9°50'N, EPR. Earth and Planetary Science Letters, 403: 15-26 DOI:10.1016/j.epsl.2014.06.015
Najman Y. 2006. The detrital record of orogenesis:A review of approaches and techniques used in the Himalayan sedimentary basins. Earth-Science Reviews, 74(1-2): 1-72
Nakagawa M, Wada K, Thordarson T, Wood CP and Gamble JA. 1999. Petrologic investigations of the 1995 and 1996 eruptions of Ruapehu volcano, New Zealand:Formation of discrete and small magma pockets and their intermittent discharge. Bulletin of Volcanology, 61(1): 15-31
Nakagawa M, Wada K and Wood CP. 2002. Mixed magmas, mush chambers and eruption triggers:Evidence from zoned clinopyroxene phenocrysts in andesitic scoria from the 1995 eruptions of Ruapehu volcano, New Zealand. Journal of Petrology, 43(12): 2279-2303 DOI:10.1093/petrology/43.12.2279
Neave DA and Putirka KD. 2017. A new clinopyroxene-liquid barometer, and implications for magma storage pressures under Icelandic rift zones. American Mineralogist, 102(4): 777-794
Pan GT, Ding J, Yao DS and Wang LQ. 2004. An Introduction of Geological Map of the Tibetan Plateau and Neighbouring Area. Chengdu: Chengdu Map Press (in Chinese)
Pan YC and Batiza R. 2003. Magmatic processes under mid-ocean ridge:A detailed mineralogic study of lavas from East Pacific Rise 9°30'N, 10°30'N, and 11°20'N. Geochemistry, Geophysics, Geosystems, 4(1): 1-33
Price RC, Smith IEM, Stewart RB, Gamble JA, Gruender K and Maas R. 2016. High-K andesite petrogenesis and crustal evolution:Evidence from mafic and ultramafic xenoliths, Egmont Volcano (Mt. Taranaki) and comparisons with Ruapehu volcano, North Island, New Zealand. Geochimica et Cosmochimica Acta, 185: 328-357 DOI:10.1016/j.gca.2015.12.009
Price RC, Smith IEM, Maas R and Stewart RB. 2017. Sub-volcanic xenoliths from Raoul Volcano:A window on petrogenetic processes in an intraoceanic subduction system. New Zealand Journal of Geology and Geophysics, 60(2): 89-111 DOI:10.1080/00288306.2016.1275708
Putirka KD, Mikaelian H, Ryerson F and Shaw H. 2003. New clinopyroxene-liquid thermobarometers for mafic, evolved, and volatile-bearing lava compositions, with applications to lavas from Tibet and the Snake River Plain, Idaho. American Mineralogist, 88(10): 1542-1554 DOI:10.2138/am-2003-1017
Putirka KD. 2008. Thermometers and barometers for volcanic systems. Reviews in Mineralogy and Geochemistry, 69(1): 61-120 DOI:10.2138/rmg.2008.69.3
Qi XX, Zhu LH, Hu ZC and Li ZQ. 2011. Zircon SHRIMP U-Pb dating and Lu-Hf isotopic composition for Early Cretaceous plutonic rocks in Tengchong block, southeastern Tibet, and its tectonic implications. Acta Petrologica Sinica, 27(11): 3409-3421 (in Chinese with English abstract)
Qi XX, Wei C, Zhang C, Zhang SQ, Hu ZC and Ji FB. 2019. Southward extension of the Bangonghu-Nujiang Suture:Evidence from Early Cretaceous intermediate and felsic magmatism in the Gaoligong orogen, China. Journal of Asian Earth Sciences, 175: 1-25 DOI:10.1016/j.jseaes.2018.09.007
Qin JZ, Huangfu G, Li Q, Qian XD, Su YJ and Cai MJ. 2000. 3-D chromatography of velocity structure in Tengchong volcano areas and nearby. Journal of Seismological Research, 23(2): 157-164 (in Chinese with English abstract)
Reubi O and Blundy J. 2009. A dearth of intermediate melts at subduction zone volcanoes and the petrogenesis of arc andesites. Nature, 461(7268): 1269-1273 DOI:10.1038/nature08510
Rudnick RL and Gao S. 2003. Composition of the continental crust. In: Holland HD and Turekian KK (eds.). Treatise on Geochemistry. Oxford: Elsevier, 1-64
Searle MP, Noble SR, Cottle JM, Waters DJ, Mitchell AHG, Hlaing T and Horstwood MSA. 2007. Tectonic evolution of the Mogok metamorphic belt, Burma (Myanmar) constrained by U-Th-Pb dating of metamorphic and magmatic rocks. Tectonics, 26(3): TC3014
Shamberger PJ and Hammer JE. 2006. Leucocratic and gabbroic xenoliths from Hualalai volcano, Hawai'i. Journal of Petrology, 47(9): 1785-1808 DOI:10.1093/petrology/egl027
Shangguan ZG, Zhao CP, Li HZ, Gao QW and Sun ML. 2004. Evolutionary characteristics of recent hydrothermal explosions in the Tengchong Rehai volcanic geothermal region. Bulletin of Mineralogy, Petrology and Geochemistry, 23(2): 124-128 (in Chinese with English abstract)
Stamper CC, Blundy JD, Arculus RJ and Melekhova E. 2014. Petrology of plutonic xenoliths and volcanic rocks from Grenada, Lesser Antilles. Journal of Petrology, 55(7): 1353-1387 DOI:10.1093/petrology/egu027
Sun SS and McDonough WF. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. In: Saunders AD and Norry MJ (eds.). Magmatism in the Ocean Basins. Geological Society, London, Special Publications, 42(1): 313-345
Tan HD, Jiang M, Lin CH, Peng M, Wu LS, Ma H, Wang ZK, Li JJ and Tan JY. 2013. Characteristics of electrical structure of Tengchong volcano-tectonic belt in Yunnan Province. Geology in China, 40(3): 800-806 (in Chinese with English abstract)
Taylor SR and McLennan SM. 1995. The geochemical evolution of the continental crust. Reviews of Geophysics, 33(2): 241-265
Tepley IFJ, Davidson JP, Tilling RI and Arth JG. 2000. Magma mixing, recharge and eruption histories recorded in plagioclase phenocrysts from El Chichón volcano, Mexico. Journal of Petrology, 41(9): 1397-1411 DOI:10.1093/petrology/41.9.1397
Tian HC, Yang W, Li SG, Ke S and Duan XZ. 2018. Low δ26Mg volcanic rocks of Tengchong in southwestern China:A deep carbon cycle induced by supercritical liquids. Geochimica et Cosmochimica Acta, 240: 191-219 DOI:10.1016/j.gca.2018.08.032
Tong W and Zhang MT. 1989. Geothermics in Tengchong. Beijing: Science Press, 1-262 (in Chinese)
Ubide T and Kamber BS. 2018. Volcanic crystals as time capsules of eruption history. Nature Communications, 9: 326 DOI:10.1038/s41467-017-02274-w
Ubide T, Mollo S, Zhao JX, Nazzari M and Scarlato P. 2019. Sector-zoned clinopyroxene as a recorder of magma history, eruption triggers, and ascent rates. Geochimica et Cosmochimica Acta, 251: 265-283 DOI:10.1016/j.gca.2019.02.021
Wang CY, Lou H, Wu JP, Bai ZM, Huangfu G and Qin JZ. 2002. Seismological study on the crustal structure of Tengchong volcano-geothermal area. Acta Seismologica Sinica, 24(3): 231-242 (in Chinese with English abstract)
Wang Y, Tang JX, Wang LQ, Duan JL, Danzhen WX, Li S and Li Z. 2018. Petrogenesis of Jurassic granitoids in the west central Lhasa subterrane, Tibet, China:The Geji example. International Geology Review, 60(9): 1155-1171 DOI:10.1080/00206814.2017.1375438
Waters LE and Lange RA. 2015. An updated calibration of the plagioclase-liquid hygrometer-thermometer applicable to basalts through rhyolites. American Mineralogist, 100(10): 2172-2184 DOI:10.2138/am-2015-5232
Workman RK and Hart SR. 2005. Major and trace element composition of the depleted MORB mantle (DMM). Earth and Planetary Science Letters, 231(1-2): 53-72 DOI:10.1016/j.epsl.2004.12.005
Wu FY, Ji WQ, Wang JG, Liu CZ, Chung SL and Clift PD. 2014. Zircon U-Pb and Hf isotopic constraints on the onset time of India-Asia collision. American Journal of Science, 314(2): 548-579
Xu HZ. 1639. Diary of the Travels of Xuxiake. In: Zhu ST and Wu Y (arranged). Shanghai: Shanghai Ancient Books Press (in Chinese)
Xu MJ, Huang H, Huang ZC, Wang P, Wang LS, Xu MJ, Mi N, Li H, Yu DY and Yuan XH. 2018a. Insight into the subducted Indian slab and origin of the Tengchong volcano in SE Tibet from receiver function analysis. Earth and Planetary Science Letters, 482: 567-579 DOI:10.1016/j.epsl.2017.11.048
Xu Y, Li XL and Wang S. 2018b. Seismic structure beneath the Tengchong volcanic area (Southwest China) from receiver function analysis. Journal of Volcanology and Geothermal Research, 357: 339-348
Yang K, Arai S, Yu JE, Yun SH, Kim JS and Hwang JY. 2012. Gabbroic xenoliths and megacrysts in the Pleisto-Holocene alkali basalts from Jeju Island, South Korea:The implications for metasomatism of the lower continental crust. Lithos, 142-143: 201-215 DOI:10.1016/j.lithos.2012.03.006
Yang QJ, Xu YG, Huang XL and Luo ZY. 2006. Geochronology and geochemistry of granites in the Gaoligong tectonic belt, western Yunnan:Tectonic implications. Acta Petrologica Sinica, 22(4): 817-834 (in Chinese with English abstract)
Yang QJ, Xu YG, Huang XL, Luo ZY and Shi YR. 2009. Geochronology and geochemistry of granites in the Tengliang area, western Yunnan:Tectonic implication. Acta Petrologica Sinica, 25(5): 1092-1104 (in Chinese with English abstract)
Yang YH, Wu FY, Xie LW, Yang JH and Zhang YB. 2009. In-situ Sr isotopic measurement of natural geological samples by LA-MC-ICP-MS. Acta Petrologica Sinica, 25(12): 3431-3441 (in Chinese with English abstract)
Yin A and Harrison TM. 2000. Geologic evolution of the Himalayan-Tibetan orogen. Annual Review of Earth and Planetary Sciences, 28(1): 211-280 DOI:10.1146/annurev.earth.28.1.211
Yu HM, Lin CY, Shi LB, Xu JD and Chen XD. 2010. Characteristics and origin of mafic and ultramafic xenoliths in trachyandesite lavas from Heikongshan volcano, Tengchong, Yunnan Province, China. Science China (Earth Sciences), 53(9): 1295-1306 DOI:10.1007/s11430-010-4044-6
Zhang L, Hu YL, Qin M, Duan Y, Duan YZ, Peng HC and Zhao H. 2015. Study on crustal and lithosphere thicknesses of Tengchong volcanic area in Yunnan. Chinese Journal of Geophysics, 58(5): 1622-1633 (in Chinese with English abstract)
Zhang ML, Guo ZF, Sano YJ, Zhang LH, Sun YT, Cheng ZH and Yang TF. 2016. Magma-derived CO2 emissions in the Tengchong volcanic field, SE Tibet:Implications for deep carbon cycle at intra-continent subduction zone. Journal of Asian Earth Sciences, 127: 76-90 DOI:10.1016/j.jseaes.2016.06.009
Zhang YT, Liu JQ and Meng FC. 2012. Geochemistry of Cenozoic volcanic rocks in Tengchong, SW China:Relationship with the uplift of the Tibetan Plateau. Island Arc, 21(4): 255-269 DOI:10.1111/j.1440-1738.2012.00819.x
Zhao YW and Fan QC. 2010. Magma origin and evolution of Maanshan volcano, Dayingshan volcano and Heikongshan volcano in Tengchong area. Acta Petrologica Sinica, 26(4): 1133-1140 (in Chinese with English abstract)
Zhou MF, Robinson PT, Wang CY, Zhao JH, Yan DP, Gao JF and Malpas J. 2012. Heterogeneous mantle source and magma differentiation of Quaternary arc-like volcanic rocks from Tengchong, SE margin of the Tibetan Plateau. Contributions to Mineralogy and Petrology, 163(5): 841-860 DOI:10.1007/s00410-011-0702-8
Zhou ZH, Xiang CY and Yang KJ. 2001. Petrological structure of the crust and upper matle in Yunnan, China. Seismology and Geology, 23(1): 69-78 (in Chinese with English abstract)
Zhu BQ, Mao CX, Lugmair GW and Macdougall JD. 1983. Isotopic and geochemical evidence for the origin of Plio-Pleistocene volcanic rocks near the Indo-Eurasian collisional margin at Tengchong, China. Earth and Planetary Science Letters, 65(2): 263-275
Zou HB, Shen CC, Fan QC and Lin K. 2014. U-series disequilibrium in young Tengchong volcanics:Recycling of mature clay sediments or mudstones into the SE Tibetan mantle. Lithos, 192-195: 132-141 DOI:10.1016/j.lithos.2014.01.017
Zou HB, Ma MJ, Fan QC, Xu B, Li SQ, Zhao YW and King DT Jr. 2017. Genesis and open-system evolution of Quaternary magmas beneath southeastern margin of Tibet:Constraints from Sr-Nd-Pb-Hf isotope systematics. Lithos, 272-273: 278-290 DOI:10.1016/j.lithos.2016.12.012
陈廷方. 2003. 云南腾冲火山岩岩石学特征. 沉积与特提斯地质, 23(4): 56-61. DOI:10.3969/j.issn.1009-3850.2003.04.009
成智慧, 郭正府, 张茂亮, 张丽红. 2012. 腾冲新生代火山区温泉CO2气体排放通量研究. 岩石学报, 28(4): 1217-1224.
成智慧, 郭正府, 张茂亮, 张丽红. 2014. 腾冲新生代火山区CO2气体释放通量及其成因. 岩石学报, 30(12): 3657-3670.
从柏林, 陈秋媛, 张儒瑗, 吴根耀, 徐平. 1994. 中国滇西腾冲新生代火山岩的成因. 中国科学(B辑), 24(4): 441-448.
樊祺诚, 刘若新, 魏海泉, 史兰斌, 隋建立. 1999. 腾冲活火山的岩浆演化. 地质论评, 45(增1): 895-904.
皇甫岗, 姜朝松. 2000. 腾冲火山研究. 昆明: 云南科技出版社, 1-418.
姜朝松. 1998. 腾冲地区新生代火山活动分期. 地震研究, 21(4): 320-329.
姜枚, 谭捍东, 张聿文, 彭淼, 李庆庆, 张立树, 许乐红, 王伟. 2012. 云南腾冲火山构造区马站-固东岩浆囊的地球物理模式. 地球学报, 33(5): 731-739.
李大明, 李齐, 陈文寄. 2000. 腾冲火山区上新世以来的火山活动. 岩石学报, 16(3): 362-370.
林木森, 彭松柏, 乔卫涛, 李昌年, 徐海军. 2014. 滇西腾冲地块新生代火山岩中高温麻粒岩包体的发现及成因. 地球科学(中国地质大学学报), 39(7): 807-819.
莫宣学, 潘桂棠. 2006. 从特提斯到青藏高原形成:构造-岩浆事件的约束. 地学前缘, 13(6): 43-51. DOI:10.3321/j.issn:1005-2321.2006.06.007
潘桂棠, 丁俊, 姚冬生, 王立全. 2004. 青藏高原及邻区地质图(1:1500000). 成都: 成都制图出版社.
戚学祥, 朱路华, 胡兆初, 李志群. 2011. 青藏高原东南缘腾冲早白垩世岩浆岩锆石SHRIMP U-Pb定年和Lu-Hf同位素组成及其构造意义. 岩石学报, 27(11): 3409-3421.
秦嘉政, 皇甫岗, 李强, 钱晓东, 苏有锦, 蔡明军. 2000. 腾冲火山及邻区速度结构的三维层析成象. 地震研究, 23(2): 157-164. DOI:10.3969/j.issn.1000-0666.2000.02.010
上官志冠, 赵慈平, 李恒忠, 高清武, 孙明良. 2004. 腾冲热海火山地热区近期水热爆炸的阶段性演化特征. 矿物岩石地球化学通报, 23(2): 124-128. DOI:10.3969/j.issn.1007-2802.2004.02.007
谭捍东, 姜枚, 林昌洪, 彭淼, 吴良士, 马欢, 王忠凯, 李晋津, 谭嘉言. 2013. 云南腾冲火山构造区的电性结构特征及其地质意义. 中国地质, 40(3): 800-806. DOI:10.3969/j.issn.1000-3657.2013.03.012
佟伟, 章铭陶. 1989. 腾冲地热. 北京: 科学出版社, 1-262.
王椿镛, 楼海, 吴建平, 白志明, 皇甫岗, 秦嘉政. 2002. 腾冲火山地热区地壳结构的地震学研究. 地震学报, 24(3): 231-242. DOI:10.3321/j.issn:0253-3782.2002.03.002
徐弘祖. 1639.徐霞客游记.绪绍唐, 吴应寿整理.上海: 上海古籍出版社
杨启军, 徐义刚, 黄小龙, 罗震宇. 2006. 高黎贡构造带花岗岩的年代学和地球化学及其构造意义. 岩石学报, 22(4): 817-834.
杨启军, 徐义刚, 黄小龙, 罗震宇, 石玉若. 2009. 滇西腾冲-梁河地区花岗岩的年代学、地球化学及其构造意义. 岩石学报, 25(5): 1092-1104.
杨岳衡, 吴福元, 谢烈文, 杨进辉, 张艳斌. 2009. 地质样品Sr同位素激光原位等离子体质谱(LA-MC-ICP-MS)测定. 岩石学报, 25(12): 3431-3441.
于红梅, 林传勇, 史兰斌, 许建东, 陈孝德. 2010. 腾冲黑空山粗安岩中镁铁质-超镁铁质包体的特征及成因初探. 中国科学(地球科学), 40(12): 1697-1709.
张龙, 胡毅力, 秦敏, 段毅, 段元泽, 彭恒初, 赵宏. 2015. 云南腾冲火山区地壳及岩石圈厚度研究. 地球物理学报, 58(5): 1622-1633.
赵勇伟, 樊祺诚. 2010. 腾冲马鞍山、打鹰山、黑空山火山岩浆来源与演化. 岩石学报, 26(4): 1133-1140.
周真恒, 向才英, 杨昆杰. 2001. 云南地壳和上地幔的岩石学结构. 地震地质, 23(1): 69-78. DOI:10.3969/j.issn.0253-4967.2001.01.009