A further discussion on geophysical characteristics of Mazhan magma pocket in Tengchong volcano-tectonic zone
-
摘要:
云南腾冲火山盆地深部存在岩浆囊是人们关注的问题,多方面进行了大量的研究且存在不同意见。作者此前依据低电阻异常推断出位于小空山、大空山、黑空山等火山口的下部,深13~30 km范围内东西向25 km,南北向30km的低电阻率的岩浆囊,其东侧为大盈江断裂。文章利用反射地震、重力、航空磁测资料进一步探讨了岩浆囊的其他地球物理特征。由于火山喷发,不仅在火山口附近形成了局部物质亏损,而且,岩浆囊本身也亏损了岩浆物质,在重力上延图中滤去了表层的火山局部负重力异常后,呈现出大空山北侧的与低阻体相对应的负重力异常,这与岩浆囊高温熔融状态含水低密度物质有密切关系。整体的重力负异常印证了低阻岩浆囊的存在。腾冲马站岩浆囊的温度为397~651℃,平均温度为524℃,其主体中心部位不可能存在磁性体。对航磁资料的处理,消除表层火山岩影响后,仍然发现有较强磁异常,其平面投影范围与低阻体相近,推断磁异常是岩浆囊顶层,在15 km深度以上范围内,随着岩浆囊体温度逐步下降,当低于480℃以后会形成新的铁磁性矿物,因此,在囊体上层出现了范围与囊体相近的较强的磁性体。反射地震仅仅在浅部较清楚地观测到较连续的反射波界面,是浅部火山岩、含水层、以及花岗岩顶界面的反映。向深部反射界面很不连续推测是由于多次岩浆上涌,其喷发时间有先后,成分有差异,故岩浆囊物质的不均匀性,虽然由于含水、矿化、熔融体构成了低电阻的共性,可温度的不均匀,却又显示了不同期次岩浆成分有变化,由岩浆囊向地表的通道上,东部花岗岩下为高速,相对于较低温度的物质,在小空山、大空山和黑空山深部沿大盈江断裂带有东西两个大的隐伏花岗岩体之间是最新的火山喷发通道,具有相对低速的通道。
Abstract:The existence of magma pocket in the depth of Tengchong volcanic basin of Yunnan Province is a topic which has aroused much attention. Extensive studies have been conducted from many aspects but there still exist different opinions. It was previously inferred by the authors based on low-resistivity anomalies that, underneath the craters of Xiaokong Mountain, Dakong Mountain and Heikong Mountain, there is a low-resistivity magma pocket with a depth range of 13-30 km and sizes of 25 km (from east to west) by 30 km (from south to north), on the east of which is Dayingjiang Fault. Based on seismic reflection survey, gravity and aeromagnetic survey data, this paper further discusses other geophysical characteristics of the magma pocket. Because of volcanic eruption, not only local material depletion is present near the craters, but the magma pocket itself depletes magma material as well; after the filtering of local negative gravity anomaly of surficial volcano, the gravity upward continuation map shows a negative gravity anomaly corresponding to the low-resistivity body on the north of Dakong Mountain, which is closely related to the depletion of aqueous low-density material from magma pocket in molten state at high temperatures. The entire negative gravity anomaly corroborates the existence of low-resistivity magma pocket. The Tengchong Mazhan magma pocket has temperatures of 397-651℃, with a mean value of 527℃, so magnetic body is impossibly present at the main central place of the pocket. After processing of the aeromagnetic data to eliminate the influence of surficial volcanic rock, large strong magnetic anomaly was still found, close to low-resistivity body in size, and it is presumed that the magnetic anomaly is at the top of the magma pocket at a depth above 15 km; as magma pocket temperature decreases gradually, new magnetic mineral will be formed below 480℃; as a result, there appears relatively strong magnetic body at the top of the pocket with a similar range to that of the pocket. Only seismic reflection at shallow depths reveals clearly that the relatively continuous reflected wave interface is a reflection of the interface among shallow volcanic rock, aquifer and granite. The reflected wave velocity map exhibits a number of velocity zones relatively continuous in vertical direction, the reflection interface towards deep formation is very discontinuous, presumably because of multiple times of magma upwelling, variable eruption times and different magma components, so the inhomogeneity of magma pocket material gives a novel geophysical model, indicating that magma pocket components are uneven; although aqueous mineralized melt constitutes a generality of low resistivity, uneven temperature shows variation in components of magma at different periods; in the pathway of magma pocket going towards ground surface, the eastern granite is underlain by high-velocity low-temperature material; and between two large concealed granite plutons in the east and west of the Dayingjiang River fault zone deeply underneath Xiaokong Mountain, Dakong Mountain and Heikong Mountain, is the newest volcanic eruption pathway, which is a low-velocity pathway.
-
1. 引言
蛇绿岩是分布在板块缝合带的一套特殊类型的镁铁质-超镁铁质岩石组合,被认为是仰冲就位于陆地的大洋岩石圈板块残片,是板块汇聚最重要的标志之一。由于蛇绿岩的岩石单元可以和现代大洋岩石圈对比,记录了大洋的形成、扩张、消减到碰撞等一系列过程,因此蛇绿岩被认为蕴含了地球内部动力学和板块运动过程的丰富信息,是研究大洋岩石圈演化的最佳载体(Dilek and Furnes, 2011;Rampone and Hofmann, 2012)。
西藏雅鲁藏布构造带是印度板块和欧亚板块的分界线。狭窄的缝合带不连续分布有一系列蛇绿岩体,西段以大规模出露新鲜地幔橄榄岩和非常少的镁铁质岩石为特征,中段以普遍蛇纹石化的地幔橄榄岩伴生相对完整的镁铁质岩石单元为特征,东段以较少出露镁铁质岩石且地幔橄榄岩中赋存大规模铬铁矿床为特征。然而,雅鲁藏布构造带蛇绿岩体的形成环境目前仍然存在较大的争议(Bao et al., 2013)。一些学者认为雅鲁藏布蛇绿岩形于洋中脊慢速扩张中心(Nicolas et al., 1981;Allegre et al., 1984;Girardeau and Mercier, 1988;吴福元等,2014;Liu et al., 2014;Zhang et al., 2016;Zhang et al., 2016, 2017;Liu et al., 2016)。然而,很多学者认为雅鲁藏布蛇绿岩形成于俯冲带环境(张旗等2003;Dubois-Cote et al., 2005;钟立峰等,2006;Bedard et al., 2009;Hebert et al., 2012;Dai et al., 2012, 2013;李文霞等,2012)。本研究展示新的地质和地球化学证据进一步探讨雅鲁藏布构造带中段日喀则蛇绿岩的形成构造环境。
2. 地质背景
青藏高原由不同的地体拼接而成,由南向北包括拉萨地体、羌唐地体、松潘—甘孜地体、昆仑—祁连地体(Yin and Harrison, 2000)。雅鲁藏布构造带位于拉萨地体和印度大陆之间,一系列断续出露的蛇绿岩体东西向展布超过2500 km(图 1a)(Gansser, 1980)。前人对这些蛇绿岩体中镁铁质岩石开展了大量高精度U-Pb锆石年代学工作,定年结果表明这些镁铁质岩石形成时代在118~149 Ma范围内(图 1a)。根据地理分布和岩体特征,雅鲁藏布蛇绿岩带从西向东可以分为3段:西段、中段和东段(图 1a)。
图 1 雅鲁藏布构造带及蛇绿岩简图Figure 1. Simplified geological map of Yarlung -Zangbo suture zone and ophiolitic massifsa-Distribution of the Yarlung Zangbo ophiolites and their zircon U-Pb ages from mafic rocks (modified from An et al., 2014; Wu Fuyuan et al., 2014); b-Distribution of the Xigaze ophiolitic massifs and sampling locations of this study (modified from Wang Xibin et al., 1987)雅鲁藏布西段蛇绿岩主要以东波、普兰、休古嘎布、仲巴、萨嘎、桑桑等岩体为代表,以出露大量地幔橄榄岩及少量镁铁质岩石为特征(徐德明等,2007;Bedard et al., 2009;刘钊等,2011;Bezard et al., 熊发挥等,2011;刘飞等,2013;2011; Liu et al., 2014; Liu et al., 2015)。东波、普兰、休古嘎布、仲巴岩体的地幔橄榄岩主要以新鲜的方辉橄榄岩为主,含少量二辉橄榄岩和纯橄岩,局部出露小规模豆荚状铬铁矿化。而萨嘎和桑桑岩体的地幔橄榄岩强烈蛇纹石化(Bedard et al., 2009)。这些岩体的地幔橄榄岩中普遍被不同规模的辉长岩和辉绿岩脉侵入,局部出露少量玄武岩。
雅鲁藏布中段蛇绿岩分布于昂仁到仁布地区,从西向东包括昂仁、吉定、路曲、下鲁、冲堆、群让、得几、白朗、大竹曲和仁布等岩体,统称为日喀则蛇绿岩(图 1a)。日喀则蛇绿岩发育最为典型,研究程度也最高(吴福元等,2014)。日喀则蛇绿岩体以广泛出露地幔橄榄岩和相对完整的镁铁质岩石层序为特征(Nicolas et al., 1981; 王希斌等,1987)。蛇绿岩层序从南向北,大致以地幔橄榄岩-辉长岩-辉绿岩-玄武岩-硅质岩的岩性变化。新鲜的地幔橄榄岩主要分布在路曲和大竹曲岩体,主要由方辉橄榄岩和少量二辉橄榄岩和纯橄岩组成,部分橄榄岩中含豆荚状铬铁矿化,而其他岩体的地幔橄榄岩多已不同程度的蛇纹石化。日喀则蛇绿岩的吉定、白朗、大竹曲等岩体发育洋壳,主要有层状辉长岩、块状辉长岩、辉绿岩席、枕状熔岩组成。而昂仁、路曲、下鲁、冲堆、群让、得几、仁布等岩体则缺失辉长岩,仅仅包含相对较薄的辉绿岩席和枕状熔岩。在壳幔过渡带可见一定规模的辉绿岩和辉长岩脉侵入到下覆的地幔橄榄岩中(陈根文等, 2003, 2015;牛晓露等,2006;李文霞等,2012;Liu et al., 2016;Zhang et al., 2016;Zhang et al., 2017)。
雅鲁藏布东段蛇绿岩主要包括泽当、罗布莎等岩体,以大量地幔橄榄岩和发育纯橄岩和铬铁矿床为特征。泽当岩体由地幔橄榄岩及一套镁铁质杂岩组成,但是与日喀则蛇绿岩不同的是,泽当岩体的地幔橄榄岩部分位于北侧,南侧出现少量辉绿岩及硅质岩。罗布莎岩体主要由地幔橄榄岩、辉长岩-辉绿岩、玄武岩组成。辉绿岩和辉长岩规模延伸较大,位于地幔橄榄岩体的北侧,玄武岩不发育,少量玄武岩位于地幔橄榄岩体的北侧(韦栋梁等,2004;钟立峰等,2006;Zhang et al., 2016;Xiong et al., 2016;刘良志等,2018)。
3. 野外地质和岩相学特征
3.1 路曲岩体
路曲岩体位于日喀则市区向南10 km左右(图 1b)。岩体北侧与古近—新近纪日喀则组磨拉石建造和第四纪堆积物呈构造接触,岩体南侧与厚层的三叠纪嘎学群复理石建造呈断层接触。路曲岩体最底部为蛇绿混杂岩,主要为强烈蛇纹石化的地幔橄榄岩和少量异剥钙榴岩脉组成(图 2a)。新鲜的橄榄岩主要位于地幔橄榄岩段中部,岩性为方辉橄榄岩和少量二辉橄榄岩。壳幔过渡带橄榄岩强烈蛇纹石化,其中被少量不同规模的辉绿岩脉侵入(图 2b)。路曲岩体镁铁质岩石主要为辉绿岩和少量枕状玄武岩(图 2f),缺失辉长岩。辉绿岩位于地幔橄榄岩和玄武岩之间,与上覆玄武岩过渡接触,与下覆地幔橄榄岩呈侵入接触关系。玄武岩呈枕状和气孔-杏仁状构造,主要位于岩体北侧,与顶部硅质岩整合接触。
图 2 日喀则蛇绿岩野外照片a—路曲岩体蛇绿混杂岩;b—路曲岩体辉绿岩脉侵入地幔橄榄岩(蛇纹岩);c—大竹曲岩体层状辉长岩;d—大竹曲岩体块状辉长岩和辉绿岩;e—大竹曲岩体含斜长石斑晶的辉绿岩;f—路曲岩体枕状玄武岩Figure 2. Field photos of the Xigaze ophiolitea-Ophiolitic mélange of the Luqu massif; b-Diabase dike that intruded into mantle peridotite (serpentinite) of the Luqu massif; c-Layered gabbro of the Dazhuqu massif; d-Massive gabbro and diabase of the Dazhuqu massif; e-Plagioclase phenocrysts in diabase of the Dazhuqu massif; f-Pillow basalts of the Luqu massif3.2 大竹曲岩体
大竹曲岩体位于日喀则蛇绿岩东部,紧邻仁布岩体,距日喀则市区80 km左右(图 1b)。岩体北侧为古近—新近纪日喀则组,南侧为侏罗—白垩系深海沉积物,为断层接触关系。大竹曲岩体的地幔橄榄岩部分主要由方辉橄榄岩和和少量纯橄岩组成。壳幔过渡带的地幔橄榄岩强烈蛇纹石化且被辉绿岩和辉长岩脉侵入。而洋壳序列底部为镁铁质-超镁铁质堆积杂岩,包括橄长岩、辉石岩、层状辉长岩、块状辉长岩,但是难以观察到这些岩石的界限和接触关系(王希斌等,1987)。层状辉长岩具典型的暗示辉长岩和浅色辉长岩交替韵律层结构(图 2c)。暗色辉长岩一般含有60%~70%单斜辉石和30%~40%斜长石,而浅色辉长岩含50%~60%斜长石和40%~50%单斜辉石。斜长石多呈自形—半自形晶,中-粗粒结构,多已蚀变,单斜辉石多呈半自形粒状,并与斜长石组成典型的辉长结构(图 3a)。
图 3 日喀则蛇绿岩镁铁质岩石显微照片a—辉长岩中半自形单斜辉石和自形-半自形斜长石构成辉长结构;b—辉长岩中粗粒单斜辉石中包裹自形—半自形的斜长石晶体;c—辉绿岩中细粒辉石和自形斜长石组成典型的辉绿结构;d—辉绿岩中自形斜长石斑晶;e—玄武岩中自形斜长石和他形单斜辉石;f—玄武岩中自形斜长石和他形单斜辉石斑晶Figure 3. Petrographic microphotos of the mafic rocks from the Xigaze ophiolitea-Gabbro texture consisting of subhedral clinopyroxene and subhedral to euhedral plagioclase in gabbro; b-Subhedral-euhedral plagioclase enclosed by coarse-grain clinopyroxene in gabbro; c-Diabasic texture consisting of fine grained clinopyroxene and euhedral plagioclase in diabase; d-Euhedral plagioclase phenocrysts in diabase; e-Fine grained plagioclase and clinopyroxene in basalt; f-Euhedral plagioclase and anhedral clinopyroxene phenocrysts in basalt常见粗大的单斜辉石晶体包裹自形—半自形的斜长石晶体(图 3b)。块状辉长岩常显示块状构造,中粒结构,岩石常被辉绿岩顺层侵入切割(图 2d)。块状辉长岩主要组成为半自形的单斜辉石和已蚀变的自形-半自形的斜长石。大竹曲辉绿岩多为顺层侵入的辉绿岩席,且这些辉绿岩席也普遍侵入到下部的地幔橄榄岩中。部分侵入地幔橄榄岩的辉绿岩发生了异剥钙榴岩化。这些辉绿岩一般粒度较细,呈典型的辉绿结构(图 3c),辉绿岩中常出现自形—半自形的斜长石斑晶(图 2e,图 3d)。在辉绿岩席的上部为一定厚度的玄武岩,这些玄武质喷出岩粒度变化较大(图 3e~f),部分玄武岩出现斑晶,以自形—半自形斜长石为主,单斜辉石较少且自形程度较低(图 3f)。
4. 分析方法
4.1 全岩主量和微量元素分析
全岩主量元素分析在澳实分析检测(广州)有限公司完成,分析仪器为荷兰帕纳科公司生产的X荧光光谱仪,型号为PANalytical MagXFast。通过XRF射线光谱仪对制成四硼酸锂玻璃进行分析,分析精度优于5%。分析方法为:首先将已扣除烧失量的准确称量的粉末样品和Li2B4O7-LiBO2助熔剂,混合均匀,转移至铂金坩埚中,在1200℃下使样品完全熔融,再将熔体倒入铂金模具,冷却后制成玻璃再上机测试。全岩氧化亚铁(FeO)含量分析在澳实分析检测(广州)有限公司完成,根据强酸消解、重铬酸钾滴定测定样品中Fe2+含量。
全岩微量元素测试分析在中国科学院地球化学研究所矿床地球化学国家重点实验室完成,测试仪器为Perkin-Elmer Sciex ELAN DRC-e电感耦合等离子体质谱仪(ICP-MS),分析精度优于10%。准确称取50 mg粉末(200目)样品放入Teflon杯中,加入1 mL亚沸HF,并在电热板上蒸干以去除大部分SiO2,再加入1 mL亚沸HF和0.5 mL亚沸HNO3,然后用不锈钢钢套密封后放入烘箱中,温度185℃加热36 h。冷却至室温后在电热板上150℃蒸干去除HF,然后再加入1 mL亚沸HNO3再次蒸干,并重复一次全部除去HF。最后加入2 mL HNO3、5 mL去离子水和500 ng Rh内标溶液,重新放入不锈钢钢套密封放入烘箱,温度140℃加热8 h,冷却至室温后吸取其中400 μL溶液移到10 mL离心管,加入去离子水定容至10 mL,最后上机测试。
4.2 全岩Pb同位素分析
Pb同位素分析在中国地质科学院地质研究所国土资源部同位素地质研究重点实验室完成。化学流程如下:吸取样品溶液0.2 mL转移至AG1-X8阴离子交换树脂,并用1 mol/L HBr溶液淋洗5次,每次1 mL,然后用2mol/L HCl溶液0.5 mL淋洗1次,然后再用6 mol/L HCl溶液1 mL淋洗1次并接Pb,最后继续用6 mol/L HCl溶液0.5 mL淋洗1次并接Pb。Pb同位素比值用多接收器等离子体质谱法(MC-ICPMS)测定,所用仪器为英国Nu Plasma HR,仪器的质量分馏以Tl同位素外标校正。样品中Tl的加入量约为铅含量的1/2。国际标样NBS 981长期测定的统计结果:208Pb/206Pb =2.16736 ± 0.00066 (± 2σ),207Pb/206Pb =0.91488 ± 0.00028 (± 2σ),206Pb/204Pb=16.9386 ± 0.0131 (± 2σ),207Pb/204Pb=15.4968±0.0107 (±2σ),208Pb/204Pb= 36.7119±0.0331 (± 2σ)。详细的化学分离方法和质谱测定方法见何学贤等(2007)。
5. 日喀则蛇绿岩镁铁质岩石地球化学特征
5.1 全岩主量元素特征
日喀则蛇绿岩路曲和大竹曲岩石样品主微量元素分析结果见表 1。路曲和大竹曲岩体的玄武岩、辉绿岩、辉长岩均有一定程度的蚀变,烧失量分别为2.28%~3.20%、2.32%~3.42 %、3.06%~3.24%。路曲镁铁质岩石的SiO2含量在49.38%~50.62%,Na2O含量为2.73% ~4.35%,K2O含量为0.25% ~1.14%。TAS投图显示路曲样品均落在玄武岩/辉长岩区域(图 4)。大竹曲岩体镁铁质岩石的SiO2含量为48.89%~51.17%,Na2O含量为2.38%~4.71%,K2O含量为0.07%~1.32 %。TAS投图显示大竹曲样品也全部落在玄武岩/辉长岩范围之内(图 4)。
表 1 日喀则蛇绿岩路曲和大竹曲岩体镁铁质岩石全岩主量元素(%)和微量元素(10-6)成分Table 1. The compositions of whole-rock major (%) and trace (10-6) elements of the mafic rocks from the Luqu and Dazhuqu massifs of the Xigaze ophiolite5.2 全岩微量元素特征
路曲和大竹曲岩石样品的球粒陨石标准化REE配分模式整体上显示LREE轻微亏损和HREE相对平坦,其特征与N-MORB非常相似(图 5a、5c)。在球粒陨石标准化REE配分模式中,路曲所有样品均无明显的Eu异常(图 5a),大竹曲岩体除了辉长岩样品显示轻微Eu正异常外,玄武岩和辉绿岩样品均没有明显的Eu异常(图 5c)。由于Rb、Ba、Cs、U、K、Sr、P等元素容易受到后期海水和热液的影响发生迁移,而高场强元素Nb、Ta、Zr、Hf、Ti和REE元素在后期热液活动过程中相对不容易发生迁移,它们主要受地幔源区和岩浆演化过程控制(Pearce et al., 2014)。因此,本次研究不选取这些容易受后期影响的元素进行标准化作图。在NMORB标准化微量元素配分模式中,所有路曲样品具有相似的微量元素特征且均显示弱Nb-Ta负异常,而Th、Zr、Hf、Ti、REE等元素特征与N-MORB一致(图 5b)。在大竹曲岩体中,大部分样品的微量元素特征保持一致,均显示弱Nb-Ta负异常(图 5d)。大竹曲玄武岩和辉绿岩的配分模式总体上相似于N-MORB,辉长岩相对于N-MORB显示略低的微量元素含量,但配分图形状接近于N-MORB。
图 5 日喀则蛇绿岩路曲和大竹曲岩石球粒陨石标准化REE配分模式(a,c)及N-MORB标准化微量元素配分模式(b,d)(球粒陨石和N-MORB标准化值据Sun and McDonough, 1989; 用于对比的OIB、E-MORB、N-MORB数值据Sun and McDonough, 1989)Figure 5. Chlondrite-normalized REE patterns (a, c) and N-MORB normalized trace elements patterns (b, d)of Xigaze ophioliteNormalized values, and N-MORB, OIB, E-MORB values after Sun and McDonough (1989)5.3 全岩Pb同位素特征
路曲和大竹曲样品全岩Pb同位素比值结果见表 2。路曲和大竹曲岩体样品的Pb同位素比值变化范围相对较小,其206Pb/204Pb的变化范围为17.54~18.05,207Pb/204Pb变化范围为15.42~15.56,208Pb/204Pb变化范围为37.67~38.17。根据样品中Th、U、Pb含量对Pb同位素比值的时间校正计算(t=140 Ma),校正结果显示(206Pb/204Pb)t的变化范围为17.46~17.94,(207Pb/204Pb)t变化范围为15.42~15.56,(208Pb/204Pb)t变化范围为37.35~38.13。在(208Pb/204Pb)t-(206Pb/204Pb)t和(207Pb/204Pb)t-(206Pb/204Pb)t图解上(图 6a、6b),路曲和大竹曲样品Pb同位素比值接近于雅鲁藏布蛇绿岩带其他岩体镁铁质岩石,且完全落入印度洋MORB区域。这些特征表明这些镁铁质岩石的地幔源区具有印度洋MORB型同位素组成特征。而雅鲁藏布蛇绿岩镁铁质岩石的Pb同位素比值则明显区别于现代大洋沉积物(Ben Othman et al., 1989)。
表 2 日喀则蛇绿岩路曲和大竹曲镁铁质岩石Pb同位素成分Table 2. The compositions of Pb isotope of mafic rocks from the Luqu and Dazhuqu massifs of the Xigaze ophiolite图 6 全岩Pb同位素比值图解a-(208Pb/204Pb)t-(206Pb/204Pb)t; b-(207Pb/204Pb)t-(206Pb/204Pb)t;矫正年龄为t=140 Ma。太平洋和北大西洋MORB和印度洋MORB范围据Mahoney et al.1998和Zhang et al. 2005。现代大洋沉积物(GLOSS, 据Ben Othman et al., 1989, 北半球参考线(NHRL)据Hart, 1984。用于对比的雅鲁藏布蛇绿岩的Pb同位素数据源自文献(牛晓露等, 2006; 钟立峰等, 2006; 徐德明等, 2007; 刘飞等, 2013; Liu et al., 2015)Figure 6. Plots of Pb isotope ratiosa-(208Pb/204Pb)tversus (206Pb/204Pb)t; b-(207Pb/204Pb)t versus (206Pb/204Pb)t at 140 Ma. Fields of Pacific and North Atlantic MORB, and India MORB after Mahoney et al. (1998) and Zhang et al. (2005). Field of GLOSS after Ben Othman et al., 1898 and NHRL after Hart (1984). The data of Pb isotope of the Yarlung-Zangbo ophiolite after Niu et al. (2006), Zhong et al. (2006), Xu et al. (2007), Liu et al.(2013, 2015)6. 讨论
6.1 地幔源区性质
由于单斜辉石相对富集HREE,因此单斜辉石分离结晶将会导致残余岩浆亏损HREE而相对富集LREE(Hauri et al., 1994)。然而,路曲和大竹曲岩石样品的球粒陨石标准化REE配分模式均显示轻微LREE亏损(图 5a,5c),说明岩浆就位之前在深部没有经历重要的单斜辉石分离结晶。此外,路曲和大竹曲镁铁质岩石的球粒陨石标准化REE配分模式没有明显Eu负异常,而Eu相容于斜长石(Aigner-Torres et al., 2007),这说明形成这些镁铁质岩石的岩浆在深部没有经历重要的斜长石分离结晶(5a,5c)。这些特征说明经地幔部分熔融产生而聚集的岩浆在上升过程中没有经历明显的结晶分异作用,因此它们的微量元素能够反演地幔源区性质和部分熔融程度。
路曲和大竹曲岩体的玄武岩、辉绿岩、辉长岩均显示相似的球粒陨石标准化REE配分模式和NMORB标准化微量元素配分模式(图 5),说明形成这些岩石的岩浆来自同一地幔源区。此外,这些岩石的REE与N-MORB的REE配分模式非常接近(图 5a,5c),且N-MORB标准化微量元素配分模式也非常平坦且接近于1(图 5b,5d),表明形成这些岩石的岩浆来自亏损地幔源区。路曲和大竹曲岩石样品的Pb同位素均落在印度洋MORB范围之内(图 6),一方面说明这些岩石起源于亏损地幔,另一方面暗示新特提斯地幔源区与印度洋地幔域有相似的地球化学特征(Mahoney et al., 1998;Xu and Castillo, 2004;Zhang et al., 2005;牛晓露等,2006)。尖晶石二辉橄榄岩部分熔融过程中熔融残留物和熔体之间不会发生中稀土(如Sm)和重稀土(如Yb)分异,但轻稀土(如La)和中稀土(如Sm)之间会随着地幔部分熔融程度的增高而发生明显分异(Aldanmaz et al., 2000)。据图 7所示,大竹曲和路曲岩体样品均落在尖晶石二辉橄榄岩熔融趋势线上,其部分熔融程度约在10%。
图 7 全岩La/Sm-Sm/Yb图解(据Aldanmaz et al., 2000)PM、N-MORB、E-MORB数据源自Sun and McDonough (1989)Figure 7. Diagram of La/Sm versus Sm/Yb (after Aldanmaz et al., 2000)PM, N-MORB, E-MORB data after Sun and McDonough (1989)6.2 形成构造环境
作为雅鲁藏布蛇绿岩带上重要的组成部分,日喀则蛇绿岩的形成环境一直是众多学者关注的焦点(吴福元等,2014)。从野外地质关系来看,路曲岩体洋壳部分缺失堆晶岩,表明路曲岩体洋壳形成过程中可能不存在岩浆房。而路曲辉绿岩大多以岩席的形式顺层侵入,甚至少量辉绿岩脉侵入到下部的地幔橄榄岩中(图 2b)。这些特征与慢速扩张脊特征一致(Purdy and Detrick, 1986)。大竹曲岩体的辉绿岩和玄武岩野外特征与路曲岩体类似,但是发育规模较小的层状辉长岩和块状辉长岩(图 2c,2d),表明大竹曲岩体洋壳形成过程中存在一定规模的岩浆房。层状辉长岩显示浅色辉长岩层和暗色辉长岩层的韵律交替结构(图 2c),表明岩浆房发生一定程度的分异。辉长岩中常见粗粒单斜辉石包裹自形斜长石晶体(图 3b),说明斜长石结晶明显早于单斜辉石。此外,大竹曲辉绿岩和玄武岩中常出现自形—半自形的斜长石斑晶(图 2d,2f),这些特征都说明斜长石较早结晶并被上升岩浆捕获到浅部。实验岩石学研究表明玄武质岩浆中水增加将降低斜长石结晶温度,从而推迟斜长石结晶时间,然而在无水环境下,斜长石则较早从玄武质岩浆中结晶(Sisson and Grove, 1993; Botcharnikov et al., 2008)。因此,可以推断形成大竹曲镁铁质岩石的岩浆是一种起源于大洋中脊扩张中心的无水玄武质岩浆。
传统上判断蛇绿岩形成于俯冲带(SSZ)环境主要是根据蛇绿岩石序列中镁铁质岩石地球化学特征而没有其他直接地质证据(Miyashiro et al., 1973; Pearce et al., 1984, 2014)。然而,从Zr-Zr/Y图解来看(图 8a),日喀则蛇绿岩路曲和大竹曲岩体大部分岩石覆盖于洋中脊玄武岩范围,明显区别于弧前玄武岩,但是不能与弧后盆地玄武岩区别。在(Nb/Yb)N-(Th/Yb)N相关图中(图 8b),路曲和大竹曲岩体岩石分布范围和趋势与洋中脊玄武岩一致,而弧后盆地玄武岩则显示明显更高(Nb/Yb)N和(Th/Yb)N比值。在(Nb/Th)N-(Th/La)N相关图中(图 8c),路曲和大竹曲样品明显区别于弧前和弧后盆地玄武岩,完全落入洋中脊玄武岩区域。这些特征再次表明路曲和大竹曲镁铁质岩体岩石形成于洋中脊环境。此外,路曲和大竹曲岩体岩石的球粒陨石标准化REE配分模式与N-MORB配分模式完全一致(图 5a,5c),也暗示这些岩石形成于洋中脊环境。在NMORB标准化微量元素配分模式图中(图 5b,5d),路曲和大竹曲镁铁质岩石微量元素特征除了相对于Th和LREE显示轻微Nb-Ta负异常,而Th、LREE、Zr、Hf、Ti总体上与N-MORB保持一致。
图 8 微量元素构造环境判别图a-Zr/Y-Zr(据Pearce, 2000);b-Nb/Yb-Th/Yb(据Pearce, 2008);c-Nb/Th-Th/La(据Godard et al., 2006);弧后盆地、洋中脊扩张中心、弧前玄武岩数据源自于PETDB数据库(http://www.petdb.org)Figure 8. Discriminant diagrams of tectonic settinga-Zr/Y versus Zr (after Pearce, 2000), b-Nb/Yb versus Th/Yb(after Pearce, 2008); c-Nb/Th versus Th/La(after Godard et al., 2006). The data of back-arc basalt, mid-oceanic ridge basalt, fore-arc basalt after PETDB database虽然路曲和大竹曲样品普遍存在轻微Nb-Ta负异常,而Nb-Ta负异常也被众多学者认为是SSZ成因的主要印记(Hebert et al., 2012;李文霞等,2012)。然而,值得注意的是路曲和大竹曲岩石的Pb同位素结果均落入印度洋MORB地幔域,而这些岩石的Nb-Ta负异常可能与印度洋MORB地幔域的地球化学不均一性有关(Moores et al., 2000)。先前的研究表明,青藏高原晚震旦世以来的蛇绿岩镁铁质岩石的同位素特征与现代印度洋MORB的同位素地球化学特征一致,明显区别于太平洋和北大西洋洋中脊玄武岩,表明现代印度洋MORB地幔域继承了特提斯洋地幔原来的空间位置和地球化学特征(Mahoney et al., 1998; Xu and Castillo, 2004; Zhang et al., 2005)。然而,在地质演化期间,特提斯洋岩石圈地幔域曾发生过大陆地壳、古老俯冲带物质混入(Castillo, 1988;Klein et al., 1988;Le Roex et al., 1989;Storey et al., 1989;Rehkamper and Hofmann, 1997)。由于在俯冲过程中Th容易迁移而Nb-Ta相对不易迁移,因此大多数与岛弧有关的岩石显示较高的Th/Nb和Th/Ta比值。此外,大陆地壳也显示较高的Th/Nb比值(Pearce et al., 2008),因此若地幔源区交代古老俯冲带物质或大陆地壳物质,则交代地幔形成的岩浆可能显示Nb-Ta负异常。而路曲和大竹曲镁铁质岩石的Pb同位素结果均落入印度洋MORB范围之内,暗示这些镁铁质岩石的地幔源区可能属于与现代印度洋MORB地幔域相似的异常地幔。因此,路曲和大竹曲镁铁质岩石中轻微的Nb-Ta负异常可能是继承而来。
7. 结论
(1)日喀则蛇绿岩路曲和大竹曲岩体镁铁质岩石的母岩浆起源亏损地幔源区,其地幔源区部分熔融程度在10%左右。这些镁铁质岩石的地幔源区的地球化学特征与印度洋MORB地幔域相似。
(2)路曲和大竹曲岩体镁铁质岩石的野外地质特征、岩相学特征和地球化学特征均支持这些岩石起源于大洋中脊扩张中心。微量元素比值也接近洋中脊玄武岩且明显区别于弧后盆地玄武岩和弧前玄武岩。N-MORB标准化微量元素配分模式显示弱Nb-Ta负异常可能是由于其地幔源区交代了古老俯冲物质。
致谢: 作为大陆科学钻探选址及钻探试验项目的组成部分, 该项工作得到了项目和课题负责人中国地质科学院地质研究所杨经绥研究员和许志琴院士, 中国科学院地质与地球物理研究所刘嘉麒院士的大力帮助和支持, 在工作过程中得到了中国地质调查局发展研究中心袁学诚教授, 中国地质大学(北京)的李金铭教授、余钦范教授、王光锷教授等,中国地质科学院地质研究所戚学祥研究员以及中国科学院地质与地球物理研究所郭正府研究员的指导和帮助, 审稿专家对论文提出了宝贵修改意见, 在此表示衷心感谢。中国地质大学(北京)的许多同学参加了野外工作和资料处理, 同时予以感谢。 -
图 1 云南腾冲火山构造区地球物理剖面位置图
图中1 线、2线、3 线(绿色)剖面为可控源音频大地电磁剖面; 1线、3 线同时为反射地震剖面; 5 线、9 线(蓝色)剖面为重磁、大地电磁剖面,AA’剖面为岩浆囊纵向剖面位置,5 线为横向剖面岩浆囊位置; 左上方为穿越小空山、大空山、黑空山的重磁剖面位置图; 背景地质图来自文献[4]
Figure 1. Location of geophysics profiles in Tengchong volcano-tectonic zone,Yunnan
where lines 1,2,and 3 (green) indicate CSAMT profiles; lines 1 and 3 indicate seismic reflection profiles,too; lines 5 and 9 (blue) indicate GM and MT profiles,respectively; profile AA’is the longitudinal profile of magma pocket,and line 5 indicates the transverse profile of magma pocket; the top left inset is the location map of GM profile traversing Xiaokong Mountain,Dakong Mountain,and Heikong Mountain; and the background geological map is after reference [4]
图 2 云南腾冲1∶20 万航磁异常平面图(a),1∶20 万航磁化极异常平面图(b),航磁化极上延2000 m异常平面图(c)和航磁化极上延5000 m异常平面图(d) (等值线距10 nT)1—ΔT磁测剖面; 2—小空山火山口磁测剖面
Figure 2. 1∶200000 aeromagnetic anomaly map (a),1∶200000RTP aeromagnetic anomaly map (b),RTP and upwardcontinued aeromagnetic anomaly map (to 2000 m) (c),andRTP and upward continued aeromagnetic anomaly map (to5000 m) (d) (contour interval: 10 nT) of Tengchong,Yunnan1-ΔT aeromagnetic survey profile; 2-Aeromagnetic survey profileacross the crater of Xiaokong Mountain
图 3 a—腾冲火山构造区小空山—黑空山布格重力异常平面图; b—腾冲重力上延500 m异常平面图; c—小空山、大空山和黑空山地面磁测异常、重力布格异常与地形高程剖面图,红线—重力异常; 棕色线—磁异常; 蓝线—沿测线测点的GPS地形高程; d—向上延拓的航空磁测异常及反演结果,计算磁性体的深度坐标考虑了上延2000 m的高度
Figure 3. a-Xiaokong Mountain—Heikong Mountain Bouguer gravity anomaly map of Tengchong volcano-tectonic area; b-Upwardcontinued gravity anomaly map of Tengchong (to 500 m); c-Ground magnetic anomaly,Bouguer gravity anomaly and terrainelevation profiles of Xiaokong Mountain,Dakong Mountain and Heikong Mountain: Red lines-Gravity anomalies,brown lines-Magnetic anomalies,blue lines-GPS terrain elevations along survey points on survey lines; d-Results of upward continuation andinversion of aeromagnetic anomalies,where upward continuation to 2000 m was considered in calculation of depth coordinates ofmagnetic objects
图 4 腾冲北线地震波速度与反射地震、大地电磁法剖面图
a—反射地震速度剖面; b—同等比例尺的反射地震叠加剖面和大地电磁法电阻率剖面叠合图
Figure 4. Maps of seismic wave velocity profile,seismic reflection profiles,and MT profile for North Tengchong line
a-Seismic reflection velocity profile; b-Superposed map of stacked seismic reflection profile and MT resistivity profile at the same scale
图 5 云南腾冲或山区马站岩浆囊横剖面地球物理模式
同等比例尺的反射地震速度剖面和大地电磁法电阻率剖面叠合图,航磁反演结果见图 3-d
Figure 5. Geophysical model for transverse profile of the magma pocket in Tengchong or mountainous Mazhan,Yunnan
See the superposed map of seismic reflection velocity profile and MT resistivity profile at the same scale,and result of aeromagnetic inversion in Fig. 3-d
-
[1] 刘嘉麒.中国火山[M].北京:科学出版社,1999:219. Liu Jiaqi.Volcanos in China[M].Beijing:Science Press,1999:219(in Chinese).
[2] 罗照华,刘嘉琦,赵慈平,等.深部流体与岩浆活动:兼论腾冲火山群的深部过程[J].岩石学报,2011,27(10):2855-2862. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201110006.htm Luo Zhaohua,Liu Jiaqi,Zhao Ciping,et al.Deep fluids and magmatism:The deep processes beneath the Tengchong volcano group[J].Acta Petrologica Sinica,2011,27(10):2855-2862(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201110006.htm
[3] 罗照华,莫宣学,侯增谦,等.青藏高原新生代形成演化的整合模型——来自火成岩的约束[J].地学前缘,2006,13(4):196-211. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200604020.htm Luo Zhaohua,Mo Xuanxue,Hou Zengqian,et al.An integrated model for the Cenozoic evolution of the Tibetan plateau:constraints from igneous rocks[J].Earth Science Frontiers,2006,13(4):196-211(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200604020.htm
[4] 戚学祥,赵宇浩,朱路华,等.青藏高原东南缘哀牢山构造带泥质高压麻粒岩的发现及其构造意义[J].岩石学报,2012,28(6):1846-1856. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201206012.htm Qi Xuexiang,Zhao Yuhao,Zhu Luhua,et al.Discovery of highpressure pelitic granulite in Ailaoshan orogenic belt,southeastern Tibet,and its tectonic implications[J].Acta Petrologica Sinica,2012,28(6):1846-1856(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201206012.htm
[5] Nicolasa,Ildefonse B.Flow mechanism and viscosity in basaltic magma chambers[J].Geophysical Research Letters,1996,23(16):2013-2016. doi: 10.1029/96GL02073
[6] 穆治国,佟伟,Garniss H.Curtis.腾冲火山活动的时代和岩浆来源问题[J].地球物理学报,1987,30(3):261-270. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX198703004.htm Mu Zhiguo,Tong Wei,GarnissH.Curtis.Times of volcanic activity and origin of magma in Tengchong geothermal area,west Yunnan province[J].Acta Geophysica Sinica,1987,30(3):261-270(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX198703004.htm
[7] 白登海,廖志杰,赵国泽,等.从MT探测结果推论腾冲热海热田的岩浆热源[J].科学通报,1994,39(4):344-347. http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB199404017.htm Bai Denghai,Liao Zhijie,Zhao Guoze,et al.The inference of magmatic heat source beneath the Rehai (Hot Sea) field of Tengchong from the result of magnetotelluric sounding[J].Chinese Science Bulletin,1994,39(4):344-347(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB199404017.htm
[8] 孙洁,徐常芳,江钊,等.滇西地区地壳上地幔电性结构与地壳构造活动的关系[J].地震地质,1989,11(1):35-45. http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ198901004.htm Sun Jie,Xu Changfang,Jiang Zhao,et al.The electrical structure of the crust and upper mantle in the west part of Yunnan province and its relation to crustal tectonics[J].Seismology and Geology,1989,11(1):35-45(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ198901004.htm
[9] 王椿镛,楼海,吴建平,等.腾冲火山地热区地壳结构的地震学研究[J].地震学报,2002,24(3):231-242. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXB200203002.htm Wang Chunyong,Lou Hai,Wu Jianping,et al.Seismological study on the crustal structure of Tengchong volcano-geothermal area[J].Acta Seismologica Sinica,2002,24(3):231-242(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-DZXB200203002.htm
[10] Nercessian A,Hirn A,Tarantola A.Three-dimensional seismic transmission prospecting of the Mont Dore volcano,France[J].Geophysical Journal of the Royal Astronomical Society,1984,76(2):307-315. doi: 10.1111/j.1365-246X.1984.tb05048.x
[11] 上官志冠,都吉夔,臧伟,等.郯庐断裂及胶辽断块区现代地热流体地球化学[J].中国科学(地球科学),1998,28(1):23-29. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK199801004.htm Shang Guan Zhiguan,Du Jikui,Zang Wei,et al.1998.Modern hot spring geochemistry at the Tanlu fault and Jiaoliao block in eastern China[J].Science in China (Series D:Earth Sciences),1998,28(1):23-29(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK199801004.htm
[12] 楼海,王椿镛,皇甫岗,等.云南腾冲火山区上部地壳三维地震速度层析成像[J].地震学报,2002,24(3):243-251. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXB200203003.htm Lou Hai,Wang Chunyong,Huang fuhang,et al.Three-Demensional Seismic Velocity Tomography of the upper crust in Tengchong volcanic area,Yunnan province[J].Acta Seismologica Sinica,2002,24(3):243-251(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-DZXB200203003.htm
[13] 谭捍东,姜枚,林昌洪,等.云南腾冲火山构造区的电性结构特征及其地质意义[J].中国地质,2013,40(3):800-805. http://geochina.cgs.gov.cn/ch/reader/view_abstract.aspx?file_no=20130312&flag=1 Tan Handong,Jiang Mei,Lin Changhong,et al.Characteristics of electrical structure of Tengchong volcano-tectonic belt in Yunnan Province[J].Geology in China,2013,40(3):800-805(in Chinese with English abstract). http://geochina.cgs.gov.cn/ch/reader/view_abstract.aspx?file_no=20130312&flag=1
[14] 郭正府,张茂亮,成智慧,等.火山"熔岩流气泡古高度计"及其在云南腾冲火山区的应用[J].岩石学报,2011,27(10):2863-2872. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201110007.htm Guo Zhengfu,Zhang Maoliang,Gheng Zhihui,et al.A link of measurements of lava flows to Palaeoelevation estimations and its application in Tengchong volcanic eruptive field in Yunnan Province (SW China)[J].Acta Petrologica Sinica,2011,27(10):2863-2872(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201110007.htm
[15] 姜朝松,周真恒,赵慈平.腾冲火山区地壳及上地幔结构特征[J].地震研究,2004,27(S1):1-6. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYJ2004S1000.htm Jiang Chaosong,Zhou Zhenheng,Zhao Ciping.The structural characteristics of the crust and upper mantle in the area of Tengchong volcano[J].Journal of Seismological Research,2004,27(S1):1-6(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-DZYJ2004S1000.htm
[16] 李成波,施行觉,刘苏苏,等.腾冲火山区的GPS形变特征[J].地球物理学进展,2007,22(3):765-770. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ200703015.htm Li Chengbo,Shi Xingjue,Liu Susu,et al.Research on the character of the GPS deformationin the Tengchong volcano area[J].Progress in Geophysics,2007,22(3):765-770(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ200703015.htm
[17] Jiang Mei,Tan Handong,Zhang Yuwen,et al.Geophysical mode of Mazhan-Gudong magma chamber in Tengchong volcanotectonic area[J].Acta Geoscientica Sinica,2012,33(5):731-739(in Chinese with English abstract). [18] 赵慈平,冉华,陈坤华.由相对地热梯度推断的腾冲火山区现存岩浆囊[J].岩石学报,2006,22(6):1517-1528. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200606010.htm Zhao Ciping,Ran Hua,Chen Kunhua.Present-day magma chambers in Tengchong volcano area inferred from relative geothermal gradient[J].Acta Petrolovica Sinica,2006,22(6):1517-1528(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200606010.htm
[19] 汪缉安,徐青,张文仁.云南大地热流及地热地质问题[J].地震地质,1990,12(4):367-377. http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ199004012.htm Wang Jian,Xu Qing,Zhang Wenren.Heat flow data and some geologic-geothermal problems in Yunnan Province[J].Seismology and Geology.1990,12(4):367-377(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ199004012.htm
[20] Clark D A.Magnetic petrophysics and magnet icpetrophy:aids to geological interpret at ion of magnetic surveys[J].Journal of Australian Geology and Geophysics,1997,17(2):83-103.
[21] Schmidt P W and Clark D A.The calculaion of magnetic components and moments from TMI:A case study from the Tuckers igneous complex,Queensland[J].Exploration Geophysics,1998,29(4):609-614. doi: 10.1071/EG998609
[22] 永田武.岩石磁学[M].北京:地质出版社.1959. Takesi Nagata.Rock-Magnetism[M].Beijing:Geological Publishing House,1959(in Chinese).
[23] 姜枚,谭捍东,等.沿雅鲁藏布江含铬蛇绿岩体、金川铜镍矿区、腾冲火山盆地深部结构的地球物理研究[M].北京:地质出版社,2015,45-60. Jiang Mei,Tan Handong,et al.Geophysics Study of the Deep Structure of Chromium-Bearing Serpentinite Masses,along the Yarlung Zangbo River Suture Zone,the Jinchuan Cu-Ni Mine Area,and the Tengchong Volcanic Basin[M].Beijing:Geological Publishing House,2015,45-60(in Chinese).
[24] Hetland E A,Wu F T,Song J L.Crustal structure in the Changbaishan volcanic area,China,determined by modeling receiver functions[J].Tectonophysics,2004,386(3/4):157-175. http://cn.bing.com/academic/profile?id=2043847102&encoded=0&v=paper_preview&mkt=zh-cn
[25] Fountain D M.Composition of the continental crust and upper mantle:A review[J].Geol.Soc.Amer.Bull.,1989,172:711-742.
[26] Mainprice D.Modelling the anisotropic seismic properties of partially molten rocks found at mid-ocean ridges[J].Tectonophysics,1997,279(1/4):161-179. http://cn.bing.com/academic/profile?id=1970709691&encoded=0&v=paper_preview&mkt=zh-cn
[27] Collier J S,Singh S C.Detailed structure of the top of the melt body beneath the East Pacific Rise at 9° 40'N from waveform inversion of seismic reflection data[J].Journal of Geophysical Research,1997,102(B9):20287-20304. doi: 10.1029/97JB01514
[28] 刘若新,魏海泉,李继泰,等.长白山天池火山近代喷发[M].北京:科学出版社,1998:108-122. Liu Ruoxin,Wei Haiquan,Li Jitai,et al.Modern Eruption of Changbaishan,Tianchi Volcano[M].Beijing:Sciece Press,1998:108-122(in Chinese).
[29] 郭履灿,马石庄,张禹慎.应用地震CT技术研究长白山火山的岩浆囊[J]. http://www.cnki.com.cn/Article/CJFDTOTAL-CTLL199601014.htm CT理论与应用研究,1996,(1):47-52.Guo Lucan,Ma Shizhuang,Zhang Yushen.Research on "Magma Chamber" of Changbai Mountain volcanoes by means of seismic tomography[J].CT Theory and Applications.1996(1):47-52(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-CTLL199601014.htm
[30] 胡亚轩,王雄.应用形变和重力资料分析腾冲火山区岩浆的活动特征[J].地震地质,2009,31(4):655-663. http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ200904010.htm Hu Yaxuan,Wang Xiong.Anaysis of the characteristics of magma activity based on the deformation and gravity measurements in Tengchong volcano region[J].Seismology and Geology,2009,31(4):655-663(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ200904010.htm
[31] 樊祺诚,刘若新,魏海泉,等.腾冲活火山的岩浆演化[J].地质论评,1999,45(S1):895-904. http://www.cnki.com.cn/Article/CJFDTOTAL-DZLP1999S1133.htm Fan Qicheng,Liu Ruoxin,Wei Haiquan,et al.The magmatic evolution of the active volcano in the Tengchong area[J].Geological Review,1999,45(S1):895-904(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-DZLP1999S1133.htm
[32] 杨经绥,许志琴,李天福,等.青藏高原拉萨地块中的大洋俯冲型榴辉岩:古特提斯洋盆的残留?[J].地质通报,2007,26(10):1277-1287. http://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD200710008.htm Yang Jingsui,Xu Zhiqin,Li Tianfu,et al.Oceanic subduction-type eclogite in the Lhasa block,Tibet,China:Remains of the Paleo-Tethys ocean basin?[J].Geological Bulletin of China,2007,26(10):1277-1287(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD200710008.htm
[33] 徐向珍,杨经绥,李天福,等.青藏高原拉萨地块松多榴辉岩的锆石SHRIMPU-Pb年龄及锆石中的包裹体[J].地质通报,2007,26(10):1340-1355. http://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD200710014.htm Xu Xiangzhen,Yang Jingsui,Li Tianfu,et al.SHRIMP U-Pb ages and inclusions of zircons from the Sumdo eclogite in the Lhasa block,Tibet,China[J].Geological Bulletin of China,2007,26(10):1340-1355(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD200710014.htm
[34] 陈松永,杨经绥,徐向珍,等.西藏拉萨地块松多榴辉岩的锆石Lu/Hf同位素研究及La-Icpms U-Pb定年[J].岩石学报,2008,24(7):1528-1538. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200807011.htm Chen Songyong,Yang Jingsui,Xu Xiangzhen,et al.Study of Lu-Hf geochemical tracing and LA-ICPMS U-Pb isotopic dating of the Sumdo eclogite from the Lhasa block,Tibet[J].Acta Petrologica Sinica,2008,24(7):1528-1538(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200807011.htm
[35] 李恒忠,杨存宝.腾冲热海地下流体观测研究[J].地震研究,2000,23(2):231-238. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYJ200002017.htm Li Hengzhong,Yang Cunbao.Underground fluid study in Tengchong hot sea[J].Journal of Seismological Research,2000,1-193(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-DZYJ200002017.htm
[36] 徐青,李翠华,汪缉安,等.云南地热资源——以腾冲地区为重点进行解剖[J].地质地球化学,1997,4:77-84. http://www.cnki.com.cn/Article/CJFDTOTAL-DZDQ199704011.htm Xu Qing,Li Cuihua,Wang Jian,et al.Geothermal resources in Tengchong region[J].Geology Geochemistry.1997,4:77-84(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-DZDQ199704011.htm
[37] 上官志冠.腾冲热海地热田热储结构与岩浆热源的温度[J].岩石学报,2000,16(1):83-90. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200001008.htm Shangguan Zhiguan.Structure of geothermal reservoirs and temperature of mantle-derived magma hot source in the Rehai area,Tengchong[J].Acta Petrologica Sinica,2000,16(1):83-90(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200001008.htm
[38] 佟伟,章铭陶.腾冲地热.北京:科学出版社,1989,1-262. Tong Wei,Zhang Mingtao.Geothermics in Tengchong[M].Beijing:Science Press,1989,1-262(in Chinese with English abstract).