西藏狮泉河蛇绿岩中侏罗世晚期(ca. 163Ma)OIB型辉绿岩及高镁闪长岩年代学及地球化学特征：早期洋壳俯冲产物？*

2019-04-04李志军李晨伟高一鸣曾敏

岩石学报 2019年3期

李志军李晨伟高一鸣曾敏, 3**

1. 成都理工大学地球科学学院，成都 6100592. 自然资源部成矿作用与资源评价重点实验室，中国地质科学院矿产资源研究所，北京 1000373. 油气藏地质及开发工程国家重点实验室，成都理工大学，成都 6100591.

板内洋岛玄武岩(OIB)的最初定义是远离构造板块边界的许多火山岛上发现的玄武岩，并认为这种构造板块内部的玄武岩喷发与板块构造作用没有直接关系(Niuetal., 2011)。自20世纪60年代首次提出地幔柱理论以来，OIB被认为与地幔柱有关(朱弟成等, 2006; 牛耀龄, 2010; 夏瑛等, 2012; Lietal., 2016)。但是近年来除被公认为典型板内地幔柱成因的洋岛型OIB(夏威夷岛链)和大火成岩省型OIB(峨眉山玄武岩)以外，地质学家在许多俯冲带周缘也发现了大量的具有和OIB地球化学性质极为相似的“富集型”基性岩。这些基性岩浆岩被解释为一些特定的俯冲过程的产物(如：洋中脊俯冲和板片断离)(Lietal., 2014b; Liuetal., 2014)。虽然，在上述俯冲过程中，OIB源区物质的来源仍有较大的争论(牛耀龄, 2010)，但是对其地球化学性质和成因的研究仍可以有效的帮助我们约束俯冲作用中的相关过程(Faccennaetal., 2010; Zhuetal., 2016)。

狮泉河-纳木错蛇绿岩带是青藏高原西部重要的蛇绿岩带之一，是研究特提斯洋东段演化及青藏高原地体增生和隆升的重要内容(Kappetal., 2003; 郑有业等, 2006; Zhuetal., 2013; Xuetal., 2014)。前人对于狮泉河-纳木错蛇绿岩带与北侧的班公湖-怒江蛇绿岩带之间的关系仍存在诸多争议(Zhangetal., 2012; Zhuetal., 2013., 2015)，而对于蛇绿岩本身岩石地球化学及年代学研究是解决这些争议的最为重要的手段。现有的研究多认为狮泉河蛇绿岩为俯冲带上(SSZ)型洋壳残余，其形成通常被认为与班公湖-怒江洋(后简称班-怒洋)壳南向或北向的俯冲作用相关(Kappetal., 2003, 2007; Liuetal., 2018)。然而，关于狮泉河蛇绿岩带所代表洋壳的时空演化过程仍需要进一步细化。

本文首次报道在狮泉河-纳木错蛇绿岩带中发现的侵位在地幔橄榄岩中的OIB型辉绿岩脉及高镁闪长岩脉。对该OIB型辉绿岩开展了锆石U-Pb年代学，对辉绿岩和闪长岩进行了岩石地球化学的研究工作，讨论其性质、构造环境和成因。研究数据将为理解狮泉河蛇绿岩所代表的洋壳演化历史提供重要依据，并为解决前人争议提供参考。

1 地质背景和样品描述

班公湖-怒江缝合带(后简称班-怒带)西起班公湖，向东经改则、安多、索县至丁青，再折而向南沿怒江直下，整体呈东西向展布的“S”型(Yin and Harrison, 2000；Zengetal., 2016a；Linetal., 2017；Liuetal., 2018；Songetal., 2018；李世民，2018)。该缝合带为南羌塘地块(北侧)和拉萨地块(南侧)的结合带，记录了班-怒洋的诞生、演变和消亡，是特提斯洋演化和盆山耦合的关键研究区域之一(Kappetal., 2003; 张玉修, 2007；Huetal., 2015; Zengetal., 2016a; Huangetal., 2017; 李世民，2018)。然而，对于洋盆开启的时间、洋盆的演化以及洋壳俯冲的极性等问题仍然存在着争议(Kappetal., 2003, 2005; Gehrelsetal., 2011; Zhuetal., 2013; Lietal., 2014a; Zengetal., 2016a;曾敏等, 2017)。

广义上的班-怒带包括位于羌塘地体和拉萨地体之间的由复理石及蛇绿混杂岩所组成的近东西展布的广泛区域(Pearce and Deng, 1988; Kappetal., 2003, 2005; Zhang and Zeng, 2018)，其垂向(南北向)宽度很大(局部>200km)，且有至少两条蛇绿岩带在其中断续出现(Girardeauetal., 1985; Kappetal., 2003; Zhuetal., 2013)，而许多研究工作所采用的班-怒带位置实际上为其北界(Kappetal., 2005; Zhangetal., 2014)。对于这两条蛇绿岩带成因的认识主要存在两种对立观点：1)形成在班-怒带内的蛇绿岩在持续构造挤压过程中向两侧发生逆冲推覆(尤其是向南向拉萨地体内部推覆)，多条蛇绿岩带均为源自班-怒带蛇绿岩的逆冲推覆席(Girardeauetal., 1985; Kappetal., 2003, 2005; 张玉修, 2007)；2)存在两条独立的蛇绿岩带，即北侧的班公湖-怒江蛇绿岩带和南侧的狮泉河-纳木错蛇绿岩带，分别代表不同性质的洋盆残余(Zhuetal., 2013; Xuetal., 2014)，二者之间主要为班-怒洋俯冲过程中改造或增生而成的新生地壳(“北拉萨地体”)(Zhuetal., 2013, 2016; Xuetal., 2014; Zengetal., 2018)。尽管存在上述争议，但几乎所有相关研究工作都把两条蛇绿岩带纳入班-怒洋演化的范畴(Kappetal., 2003, 2005; 张玉修, 2007; Zhangetal., 2012; Zhuetal., 2013; Xuetal., 2014)。对于班-怒洋的演化过程，大部分学者认为其经历了完整的威尔逊旋回，但是旋回各阶段的时限以及性质仍然存在争议(曾敏等, 2017)。因此，对班-怒带蛇绿岩及相关地质体开展进一步的深入研究是解决上述争议的关键，而此争议的解决是厘清班-怒带时空演化过程及蛇绿岩就位机制的关键，对理解与班-怒带相关的其它疑问也有重要意义。

图1 工作区地质简图(据Liu et al., 2018; Zhu et al., 2013)JSSZ=金沙缝合带；LSSZ=龙木错-双湖缝合带；BNSZ=班公湖-怒江缝合带；SNMZ=狮泉河-纳木错蛇绿岩带；LMF=罗巴迪-米拉山断层；IYZSZ=雅鲁藏布江缝合带.红色五角星为本次狮泉河辉绿岩和闪长岩的采样地点;黑色五角星为研究区周缘可参考的前人年龄数据Fig.1 Simplified geologic map of the work area(after Liu et al., 2018; Zhu et al., 2013)JSSZ=Jinsha Suture Zone; LSSZ=Longmu Tso-Shuanghu Suture Zone; BNSZ=Bangong-Nujiang Suture Zone; SNMZ=Shiquanhe-Nam Tso Ophiolitic Mélange Zone; LMF=Luobadui-Milashan Fault; IYZSZ=Indus-Yarlung Zangbo Suture Zone. The red five-pointed star in the picture is the sampling location of the Shiquanhe diabase and diorite in this study; while the black five-pointed star is the predecessor age data for the study area

图2 狮泉河辉绿岩及闪长岩野外露头照片及镜下照片(a)侵位于橄榄岩中的透镜体化的辉绿岩；(b)辉绿岩脉与橄榄岩接触界面的冷凝边；(c)辉绿岩中的斜长石(Pl)斑晶；(d)辉绿岩中的辉绿结构；(e)侵位于橄榄岩中的闪长岩脉；(f)闪长岩中的辉石(Py)斑晶Fig.2 Shiquanhe diabase and diorite outcrop and microscopic photos(a) lenticular diabase dike in the peridotite; (b)chilled margin between the diabase and the peridotite; (c)plagioclase (Pl) phenocryst in the diabase; (d)diabaseic texture in the diabase; (e)diorite dike in the peridotite; (f)pyroxene (Py) phenocrysts in the diorite

表1狮泉河地区OIB型辉绿岩锆石U-Th-Pb同位素测定结果

Table 1 U-Th-Pb isotopic composition of zircons from OIB-type diabase in Shiquanhe area

测点号含量(×10-6)PbThUTh/U同位素比值年龄(Ma)207Pb/235U1sigma206Pb/238U1sigma207Pb/235U1sigma206Pb/238U1sigma协和度SQH-03-01175.3181237810.480.1670.0060.0260156.55.464164.51.92495%SQH-03-02153.4161030320.530.1760.0060.0260164.55.267165.51.90799%SQH-03-03677.6198939380.511.4970.0320.1570.002929.312.94948.38.85998%SQH-03-04172.8176836710.480.1670.0060.0260156.65.161163.21.73895%SQH-03-05154.1161030250.530.1770.0080.0260165.86.754164.32.11499%SQH-03-06154.5161930310.530.1760.0080.0260164.26.865162.42.17298%SQH-03-13172.9177536680.480.1640.0060.0250154.05.263160.61.76995%SQH-03-07375.4969.113690.711.1470.030.130.001785.714.17796.98.00798%SQH-03-08165.3169730650.550.1830.0070.0260170.45.784165.91.87997%SQH-03-09171.7173135620.490.1640.0060.0260153.85.574163.51.73793%SQH-03-10158.3165530250.550.1770.0060.0260165.45.455164.61.75499%SQH-03-11173.3175235950.490.1640.0070.0260154.15.745163.81.78893%SQH-03-12156.2164930350.540.1780.0070.0260166.05.695163.81.85598%SQH-03-14367.6480.5740.10.652.6840.0660.2190.002122318.11124912.3496%SQH-03-15153.6161430350.530.1730.0060.0260162.04.931163.31.78999%

本研究发现的辉绿岩脉位于西藏阿里地区狮泉河噶尔县县城北北西约10km处(图1；N32°34′42″、E80°03′42″)，从构造背景上来说属于狮泉河-纳木错蛇绿岩带(图1)。辉绿岩脉侵位于狮泉河蛇绿岩带的蛇纹石化橄榄岩中，且已变形成透镜状(图2a)，透镜体形态指示了南北向挤压缩短为主的构造变形。辉绿岩与橄榄岩的接触带可见明显的热接触冷凝边(图2b)，同样表明其侵入接触关系。辉绿岩镜下含少量斜长石斑晶(图2c)，基质主要为斜长石和单斜辉石，基质中斜长石较为自型，具有辉绿结构(图2d)。副矿物主要有钛铁矿和磁铁矿等。此外，在辉绿岩脉北北东向约200m处发现有同样侵位于地幔橄榄岩中的闪长岩脉(图1；采样位置坐标为N32°34′56″、E80°03′25″)(图2e)。闪长岩镜下具斑状结构，斑晶主要为角闪石，基质主要为斜长石、角闪石以及少量的黑云母(图2f)。

2 样品及测试分析方法

本次研究从狮泉河北部蛇绿岩内采集了约20kg辉绿岩样品用于进行了锆石U-Pb定年。为保证样品的可靠性，所采集的样品均去除表面杂质和次生风化薄壳以避免外部不相关锆石的混入。U-Pb锆石年代学研究样品首先采用重液分选和磁选将粉碎过后的样品进行初步分离的预处理。然后在双目显微镜下手工挑选出15颗锆石颗粒，所选锆石颗粒为自形到半自形，晶体长度为50～100μm，长宽比为2:1到3:1。挑选出的锆石颗粒被粘贴制成环氧树脂锆石靶并打磨抛光使其露出中心部位。最后，结合反射光、透射光和阴极发光图像观察锆石的内部结构并选择合适的U-Pb定年分析点，所选分析点(除3颗继承锆石)在透射光下无包裹体和裂隙，反射光下边缘无溶蚀，阴极发光图像下显示出明显的锆石震荡环带。锆石的分选、制靶和阴极发光采集均在北京锆年领航科技服务公司进行；锆石U-Pb同位素定年在中国科学院贵阳地球化学研究所利用Agilent 7500a LA-ICP-MS质谱仪分析完成，激光剥蚀束斑大小为24μm，采用91500标样作为锆石定年标样，NIST SRM 610标样作为锆石微量元素矫正标样。单个锆石测试数据的误差和206Pb/238U年龄的加权平均值误差均为1σ，由于所测试样品年龄较年轻故采用206Pb/238U年龄作为使用年龄，分析结果见表1。详细的仪器操作条件和数据处理方法参见Liuetal. ( 2008, 2010)。经过普通铅校正后的数据使用Isoplot3 (Ludwig, 2003；Andersen, 2002)进行年龄的拟合计算。

此外，在辉绿岩脉的不同位置选取了5个样品，在细晶闪长岩脉中采集了2个样品以进行全岩地球化学分析。为保证测试分析数据的可靠性，所采样品均为无蚀变的新鲜样品，同时采样过程尽量避开了岩石内部的后期微小细脉。样品的全岩地球化学分析在西南冶金地质测试中心完成。样品粉碎前先用无水乙醇清洁其表面，然后将样品粉碎成200目粉末并加热至920℃持续2h来计算灼烧量(LOI)。主量元素用Rigaku RIX2100型X荧光光谱仪分析完成，分析精度优于1%；微量元素通过ICP-MS测定，仪器型号为Agilent 7500a质谱仪，详细的样品消解处理过程、分析精密度和准确度同Liuetal. (2008)，分析结果见表2。

表2狮泉河辉绿岩和闪长岩主量元素(wt%)和微量元素(×10-6)数据

Table 2 Whole-rock major (wt%) and trace (×10-6) element data of the diabase and diorite samples in Shiquanhe

样品号SQH-03-S2SQH-03-S3SQH-03-S4SQH-03-S5SQH-03-S6SQH-06-S1SQH-06-S2岩性辉绿岩闪长岩SiO247.8340.9344.5749.8849.4555.3956.81Al2O316.6520.3517.7417.8818.5015.6916.13Fe2O31.890.601.861.170.491.622.00FeO7.9611.039.358.859.414.654.06MgO8.197.836.314.755.148.827.07CaO5.696.558.136.166.324.083.91Na2O4.416.874.745.986.395.236.08K2O3.291.812.481.010.712.231.94TiO21.962.272.262.021.981.000.94MnO0.180.220.210.190.190.110.10P2O50.510.630.590.520.530.230.21LOI0.930.890.910.930.940.960.98FeOT9.6611.5811.039.909.866.115.86Mg#60.454.950.846.348.472.268.5FeOT/MgO1.181.481.752.091.920.690.83MgO/SiO20.170.190.140.100.100.160.12Nb/Zr0.130.130.130.130.130.070.06La/Yb 12.113.612.912.914.012.414.1Nb/Th7.808.5710.09.399.703.112.86Nb/Y1.221.251.171.271.320.730.66La/Sm3.833.873.783.863.834.675.02Nb/Zr0.130.130.130.130.130.070.06La23.1129.7830.1725.1722.9513.8816.39Ce47.4262.1963.9953.5947.9726.4529.78Pr6.808.698.927.336.703.523.87Nd31.9639.8941.6434.7331.3315.7117.02Sm6.037.697.976.515.992.973.27Eu1.832.232.331.981.790.920.98Gd5.617.037.166.105.572.963.17Tb0.780.981.030.860.760.420.46Dy4.145.125.284.453.942.272.42Ho0.800.991.010.840.760.450.48Er2.052.512.592.131.941.181.26Tm0.340.400.410.350.290.190.20Yb1.912.192.341.961.641.121.16Lu0.270.310.320.270.210.160.17Y19.8924.4324.6820.4318.8311.7612.45Ba657232697883757218210Co25.0931.9027.6824.9624.4526.2522.26Cr13312838.7527.9124.52164154Hf5.276.166.025.635.173.733.71Ni71.6381.7112.1510.038.12145.090.51Pb27.4119.8826.4124.7224.4829.7016.70Rb11065.1083.5729.8318.2143.0431.97Sr392230388465483337335Ta1.732.162.071.981.710.680.63Th3.113.572.892.762.562.772.86U0.830.870.900.770.690.710.70V222270259235238131126Zr188230218197192131133

3 分析结果

3.1 锆石U-Pb年代学

通常较高的Th/U比值(>0.4)和震荡环带被认为是岩浆成因锆石的特征，而无环带或弱环带和低的Th/U比值(<0.1)被认为是变质成因的锆石特征(Hoskin and Ireland, 2000; Wangetal., 2016)。本次所分析的锆石样品具有较高的Th/U比值(0.48～0.71)(表1)，以及部分明显的震荡环带(图3a)。在锆石的球粒陨石标准化稀土元素配分模式图上(图3b)，锆石的轻稀土元素(LREE)亏损，重稀土元素(HREE)富集，并具有明显的Ce正异常(Ce/Ce*平均值=65.68)和Eu负异常(Eu/Eu*平均值=0.22)，显示出岩浆锆石的成分特征(Hoskin and Schaltegger, 2003)。前人研究表明，岩浆锆石的微量元素含量从超基性岩，经过基性岩到中酸性岩石,其含量总体上具有增长的趋势(∑REE=50×10-6～10000×10-6)，其中基性岩中锆石REE总量约为2000×10-6(Hoskin, 2005; 赵志丹等，2018)。因此，狮泉河辉绿岩中的锆石极可能为基性岩浆成因的锆石(∑REE平均值=2327×10-6)。此外，对锆石部分微量元素及REE比值的判别((Sm/La)N=96～915；Ce/Ce*=12～215)也显示出了本次研究所涉及的锆石大部分为同岩浆期的岩浆锆石(Hoskin, 2005; 赵志丹等，2018)。

图3 狮泉河OIB型辉绿岩锆石阴极发光及测试点位置图(a)、锆石的球粒陨石标准化REE模式图(b，标准化值据Sun and McDonough, 1989)及锆石U-Pb年龄协和图(c)Fig.3 CL images of zircons from OIB-type diabases in Shiquanhe (a), chondrite-normalized REE patterns for dated zircons (b, normalization values after Sun and McDonough, 1989) and U-Pb zircon concordia diagram (c)

本次研究所选取的15颗锆石中有3颗锆石有“核心-环带结构”(图3a)，并且在锆石核心部位分别产生了1249±12.3Ma(协和度=96%)、948±8.8Ma(协和度=98%)和796±8.8Ma(协和度=98%)的206Pb/238U年龄，推测为岩浆上涌过程中捕获锆石(Wangetal., 2016)。同时，12颗岩浆锆石测点均投影于谐和线上或谐和线附近(图3c)，计算出了163.7±0.54Ma(n=12，MSWD=3.8)的206Pb/238U年龄加权平均值，该年龄基本可以代表矿物结晶年龄，即辉绿岩脉的侵入年龄。

3.2 全岩地球化学

本次分析的7件全岩地球化学样品的烧失量均小于1%，结合野外调查和镜下观察结果(图2)，认为所有样品仅发生了轻微程度的蚀变。因此，对于所有样品，主微量元素均重新归一化为无水的质量百分比，并且仅使用不易迁移的微量元素(如，高场强元素(HFSE)和稀土元素(REE))讨论上述样品的岩石学成因和演化过程(Wangetal., 2016)。由于岩浆来源相同的次火山岩和火山岩有相似的地球化学组成和性质(Daietal., 2012)，同时前人已有大量研究表明全岩地球化学性质相同的高镁安山岩与高镁闪长岩以及OIB型玄武岩和有着OIB化学性质的辉绿岩拥有相同的岩浆起源(Kameietal., 2004; Zhouetal., 2009; 邓晋福等, 2010; Lietal., 2016)。此外，这些化学性质相同的次火山岩和火山岩常同时用于讨论岩石成因和构造背景(Kameietal., 2004; Lietal., 2016)。因此，本次研究选用了一些火山岩的地球化学图解用于讨论上述次火山岩性质的细晶闪长岩和辉绿岩样品。

3.2.1 闪长岩

本次所采集的细晶闪长岩样品的SiO2含量为55%～57%，K2O含量为1.9%～2.2%，Al2O3含量为15.7%～16.1%，MgO含量为7.1%～8.8%，FeOT(=FeO+ 0.9×Fe2O3)含量为5.86%～6.11%，Mg#为68.5～72.2，TiO2变化于0.94%～1.0%。在SiO2-Zr/TiO2图解中(图4a)，上述样品被分类为安山岩/闪长岩; 在FeOT/MgO-SiO2图解中显示为钙碱性的特征(图4b); 在K2O-SiO2图解中所有样品均显示出了高K的特征，显示为高钾钙碱性系列(图4c); 在SiO2-Mg#图解上(图4d)，所有闪长岩样品均投点到高镁安山岩/闪长岩(HMA)区域。

图4 狮泉河闪长岩SiO2-Zr/TiO2图(a, 据Winchester and Floyd, 1977)、FeOT/MgO-SiO2图(b, 据Miyashiro, 1974)、 K2O-SiO2图(c, 底图据Liu et al., 2018)和Mg#-SiO2图(d, 底图据Liu et al., 2018)蓝色圆形为本次狮泉河闪长岩样品;橙色菱形为赞岐岩型高镁安山岩样品(数据据Tatsumi, 2006)，蓝色方形为Liu et al. (2018)所报道的狮泉河高镁安山岩样品Fig.4 Diagrams of SiO2 vs. Zr/TiO2(a, after Winchester and Floyd, 1977), FeOT/MgO vs. SiO2(b, after Miyashiro, 1974), K2O vs. SiO2(c, after Liu et al., 2018)and Mg# vs. SiO2(d, 据Liu et al., 2018) for diorite in the Shiquanhe areaThe blue circle is the diorite sample in the Shiquanhe area in this study; the orange diamond is the sample of the sanukite-type high-magnesium andesite (data from Tatsumi, 2006); the blue diamond is the lion spring high magnesium reported by Liu et al. (2018)

此外，它们的稀土元素分配曲线为轻稀土富集型((La/Yb)N=8.9～10.1)，基本无铕异常(Eu*=0.92～0.93)(图5a)。在微量元素蛛网图上(图5b)，岩石富集Rb、Ba、Th等大离子亲石元素(LILE)，并显示出Nb负异常的典型大陆地壳微量元素配分模式。同时，岩石具有有较高的Cr含量(153×10-6～163×10-6)和Ni含量(90.5×10-6～145×10-6)。

图5 狮泉河样品的球粒陨石标准化稀土元素配分模式图(a)和N-MORB标准化微量元素蛛网图(b)数据来源：球粒陨石、N-MORB、全球平均OIB(Sun and McDonough, 1989);洋脊俯冲成因的玄武岩(Li et al., 2016);赞岐岩型高镁安山岩(Liu et al., 2018; Tatsumi, 2006);初始俯冲弧前玄武岩(Stern and Gerya, 2018) Fig.5 Chondrite-normalized REE pattern (a) and primitive mantle-normalized trace-element spidergram (b) for the samples in the Shiquanhe areaData source: chondrite, N-MORB, global average OIB (Sun and McDonough, 1989); basalts formed by oceanic ridge subduction (Li et al., 2016); high-magnesia andesite associated with the sanukite rock (Liu et al., 2018; Tatsumi, 2006); fore-arc basalts associated with subduction initiation (Stern and Gerya, 2018)

上述这些地球化学特征与日本西南部Setouchi火山带中的新生代高镁安山岩(Miyashiro, 1974; Shimodaetal., 1998; Tatsumi and Hanyu, 2003)，青藏高原北部的晚三叠世Hohxil高镁安山岩(Wangetal., 2011)，青藏高原BNSZ中段的中侏罗世Darutso高镁安山岩(Zengetal., 2016b)以及西藏狮泉河北部的高镁安山岩(Liuetal., 2018)类似，符合高镁安山岩/闪长岩(HMA)的性质。

3.2.2 辉绿岩

本次选取的辉绿岩样品的SiO2含量较低且变化较大(40.93%～49.88%)，并伴随有较高含量的MgO(4.75%～8.19%)和Mg#(=Mg/(Mg+Fe)×100；46.32～60.42)，TiO2含量变化于1.96%～2.27%，P2O5含量变化于0.51%～0.63%，FeOT含量变化于9.66%～11.58%，Al2O3含量变化于16.65%～20.35%，Na2O含量变化于4.41%～6.87%，K2O含量变化于0.70%～3.29%。同时，所测样品含有较高含量的Cr(25×10-6～133×10-6)和Ni(8.1×10-6～81.7×10-6)。在Zr/TiO2-Nb/Y图解中(图6)，所有样品均投点到碱性玄武岩中。

图6 狮泉河辉绿岩Nb/Y-Zr/TiO2图解(据Winchester and Floyd, 1977)红色三角型为本次辉绿岩样品；绿色圆形为本次闪长岩样品；蓝色菱形为夏威夷玄武岩样品Fig.6 Nb /Y-Zr / TiO2 diagram of the diabase dike in Shiquanhe area (after Winchester and Floyd, 1977)The red triangle is the diabase sample; the green circle is the diorite sample; the blue diamond is the Hawaiian basalt sample

在球粒陨石标准化的稀土元素配分模式图(图5a)和N-MORB标准化的微量元素蛛网图(图5b)中，所有样品均富集轻稀土元素((La/Yb)N=8.67～10.1)和高场强元素(HFSEs)，同时伴随有轻微的Nb负异常和强烈的Pb正异常，无明显Eu异常(Eu/Eu*=0.91～0.95)。上述特征与已证实为洋脊俯冲成因的的班公湖怒江缝合带中段的OIB型基性岩相似(Lietal., 2016)。

4 讨论

4.1 形成时代及其地质意义

从狮泉河北部辉绿岩脉中挑出的12颗岩浆锆石的加权平均年龄为163.7±0.54Ma(n=12，MSWD=3.8)。指示了该辉绿岩在俯冲带增生楔中的就位时间。该年龄与近年来狮泉河-纳木错蛇绿岩带附近报道的俯冲带相关的中-晚侏罗世的各类岩浆岩年龄较为接近(图1)，如：165.8±1.7Ma(闪长岩；郑有业等, 2006)、163.35±0.75Ma(辉长岩-闪长岩；郑有业等, 2006)、160.8±2.3Ma(高镁安山岩；Liuetal., 2018)。上述时代接近的不同岩浆岩被解释为不同的成因类型，其中Liuetal. (2018)所报道的160.8±2.3Ma的高镁安山岩被解释为俯冲沉积物的部分熔融和年轻洋壳的部分熔融作为物源的赞岐岩型高镁安山岩。

4.2 岩石形成的P-T条件

4.2.1 闪长岩

在岩浆岩源区结晶分异的过程中，许多微量元素的浓度和比值取决于残余矿物相的性质(Johnson and Plank, 2000; Spandleretal., 2008; Hermann and Rubatto, 2009; Skora and Blundy, 2010; Liuetal., 2018)。例如，在锆石中由于Zr元素的分配系数高于Hf元素(Blundy and Wood, 2003)，所以，球粒陨石标准化后Zr/Hf比值则能指示了部分熔融后岩浆岩中锆石的相对含量(Blundy and Wood, 2003; Nebeletal., 2011)。锆石由于存在较高的饱和结晶温度(731～846℃; Watson and Harrison, 1983)，所以本次研究的闪长岩中较高(Zr/Hf)N(0.97～0.99；球粒陨石标准化值，Sun and McDonough, 1989)反映了闪长岩形成于较高温的环境下。

稀土元素在角闪石中具有高度的相容性，有着角闪石残余相的部分熔融的俯冲洋壳可以产生较高的Zr/Sm比值和较低的Sm/Hf比值(Green, 1994; Sisson, 1994)。本次研究的闪长岩与全球平均OIB相比(Sm/Hf =1.42，Zr/Sm=22.49；Plank and Langmuir, 1998)具有的较低的Nb/Ta比值(12.8～13.0)，Hf/Sm比值(1.13～1.25)和较高的Zr/Sm比值(40.6～44.0)，指示了岩浆形成于有残余角闪石相的低压条件下(Hanyuetal., 2006)。因此，本文认为本次研究的闪长岩形成于高温-低压环境。

图7 赞岐岩与埃达克岩、玻安岩和bajaite 的对比(据Kamei et al., 2004)Fig.7 Comparison of sanukitic HMA to adakitic HMA, boninitic HMA and bajaitic HMA(after Kamei et al., 2004)

4.2.2 辉绿岩

通常，地幔衍生的初级熔体具有较高的Ni含量(>400×10-6)、Cr含量(>1000×10-6)和Mg#(73～81)(Wilson, 2007)，而本次研究的OIB型辉绿岩显示出了变化较大的Mg#值(46.3～60.4)，并伴随有变化较大的相容元素含量：Cr=25×10-6～133×10-6)，Ni=8.1×10-6～81.7×10-6)，表明橄榄石已经发生了显着的分馏。同时，在辉绿岩中观察到的主要斑晶是斜长石斑晶(图2c)，指示了低压分馏阶段(Tangetal., 2012)。

4.3 岩石成因

4.3.1 HMA型闪长岩

通过全岩地球化学分析结果认为本次发现的侵位于超基性岩内的闪长岩脉为高镁安山岩/闪长岩(HMA)性质，且与最近前人所报道的狮泉河北部高镁安山岩(Liuetal., 2018)性质较为接近。高镁安山岩的岩石成因目前有几种模型：1)俯冲板部分熔融后与上覆增生楔的相互作用(Hanyuetal., 2003, 2006; Wangetal., 2011, Zengetal., 2016b)；2)拆沉后的下地壳在软流圈中熔融形成的熔体在上升过程中与地幔橄榄岩的相互作用(Xuetal., 2002; Gaoetal., 2004)；3)含水地幔橄榄岩的部分熔融(Wood and Turner, 2009; Straubetal., 2011)；4)不同性质的岩浆的混合作用(Qian and Hermann, 2010)；5)部分熔融的俯冲海沟沉积物与地幔橄榄岩之间相互作用(Zengetal., 2016b; Liuetal., 2018)；6)较为年轻的热的洋壳的部分熔融(Yogodzinskietal., 1995; Furukawa and Tatsumi, 1999; Kameietal., 2004)。这几种成因模式主要产生了4种高镁安山岩的岩石类型：赞岐岩型、埃达克型、巴哈岩型和玻安岩型(Kay, 1978; Rogersetal., 1985; Bloomer and Hawkins, 1987; Tatsumi, 2006;唐功建和王强, 2010)。与赞岐岩型相比：埃达克型和巴哈岩型具有较高的Sr、Sr/Y和(La/Yb)N值(Rogersetal., 1985; Martin, 1999)，玻安岩型具有较低的TiO2、HFSE、REE和大离子亲石元素(LILE)(Bloomer and Hawkins, 1987)，通过对上述不同地球化学特征的对比(图7)，本次研究的狮泉河北部的高镁HMA具有较高的TiO2含量，Y含量和Yb含量以及较低的Sr/Y比值和La/Yb比值(图7)，结合微量元素蛛网图的对比(图5)，本文认为狮泉河高镁安山岩的特征与出露于日本中新世Setouchi火山岩带的赞岐岩极为相似。

赞岐岩已被公认为是熔融俯冲板片沉积物熔体与上覆地幔楔相互作用的产物(Shimodaetal., 1998; Wangetal., 2011)。本文认为狮泉河赞岐岩型闪长岩也是这种教科书模式的成因，该推断同时基于以下证据：1)Th元素和Ta元素大量富集于海相沉积物中，同时U元素和Nb元素也与海相沉积物密切相关(Plank and Langmuir, 1998)。因此，本次研究的HMA中较低的U/Th(0.24～0.25)和Nb/Ta(12.7～13.0)具有极其类似于全球俯冲沉积物(GLOSS)的特点(GLOSS：U/Th=0.24，Nb/Ta=14.19；Plank and Langmuir，1998)；2)根据Labaniehetal. (2018)对小安地列斯群岛Martinique岛弧的研究结果，La/Sm比值可以指示混入到地幔楔中熔融的海沟沉积物的相对含量。因此，本次研究较高的La/Sm比值指示了海沟沉积物部分熔融为HMA提供了部分物源；3)本次研究的HMA的MgO、Cr和Ni的含量均较高，表明其源区有地幔的部分熔融(Stern and Hanson, 1991; Smithies and Champion, 2000)，而不可能为单一的洋壳沉积物熔融体。

图8 狮泉河辉绿岩和闪长岩的Harker图解红色三角形为本次研究的辉绿岩，蓝色圆形为闪长岩Fig.8 Harker diagrams of the diabase and the diorite in shiquanhe The red triangle is the diabase of this study, and the blue circle is diorite

通过上述分析，笔者认为狮泉河赞岐岩型的闪长岩为俯冲沉积物与地幔熔体相互作用的结果。俯冲带正常的地热梯度通常只会在俯冲板块(包括沉积物)开始融化之前使洋壳脱水(Tatsumi, 2006)。因此，通过沉积物衍生的熔体与年轻上覆地幔楔的相互作用形成的高镁安山岩(赞岐岩)只能存在于异常的地热梯度区域，如洋脊俯冲的板片窗环境(Wangetal., 2011)和俯冲早期阶段的弧前热异常环境(Tatsumi, 2006；Zengetal., 2016b)。但洋脊俯冲模式下的高镁安山岩会有埃达克岩的地球化学性质(Tangetal., 2010)，该特征与本次研究的闪长岩不符。因此我们在此提出，狮泉河赞岐岩型闪长岩形成于俯冲早期阶段的弧前热异常环境，并且很可能是狮泉河蛇绿岩带的初始俯冲阶段。在较高的地热梯度下，板片的熔融将会发生在俯冲的较浅部地区，并伴随含白云母矿物的熔融产生贫Sr的安山质岩浆(Weinberg and Hasalová, 2015)，该特征符合本次研究的闪长岩的特征(图5)。此外，狮泉河闪长岩较高的(Zr/Hf)N(0.97～0.99)、Zr/Sm比值(40.6～44.0)和较低的Hf/Sm比值推测出的高温-低压成岩环境也与初始俯冲模型中的温压条件相符合。

4.3.2 OIB型辉绿岩

OIB型玄武岩浆的形成和演化，涉及到来自地幔柱、软流圈、岩石圈地幔和地壳等不同端元组分的贡献(徐义刚, 2002)。因此，要分析其地幔源区的性质，首先要考察地壳物质对岩浆源区的影响(朱弟成等, 2006)。在N-MORB标准化的微量元素蛛网图上可以很明显的看出，本次OIB型的玄武岩与全球平均的OIB型的玄武岩相比有明显的Nb负异常、轻微的Ti负异常、强烈Pb正异常以及变化的Sr含量与Sr负异常(图5b)。该特征反映了岩浆作用过程中的地壳混染以及岩石圈拉张过程中的岩石圈富集组分部分熔融(Huangetal., 2013;杨帆等, 2018)。

图9 狮泉河辉绿岩不活动微量元素构造环境判别图解Fig.9 Tectonic discrimination diagrams based on immobile trace elements for the diabase in the Shiquanhe area(a) Zr/Y vs. Zr diagram (Pearce and Norry, 1979), A=volcanic arc basalt; B=mid-ocean ridge basalt; C=within-plate basalt; D=mid-ocean ridge basalt and volcanic arc basalt; E=mid-ocean ridge basalt and within-plate basalt. (b) TiO2/Yb vs. Nb/Yb diagram (Pearce et al., 1992). (c) 2×Nb-Zr/4-Y diagram (Meschede, 1986), A1=within-plate alkaline basalt; A2=within-plate alkaline basalt and intraplate tholeiite; B=enriched mid-ocean ridge basalt; C=intraplate tholeiite and volcanic arc basalt; D=volcanic arc basalt and normal-mid-ocean ridge basalt. (d) Ti/100-Zr-Y×3 diagram (Pearce and Cann, 1973), A=island arc tholeiite; B=mid-ocean ridge basalt; C=calc-alkaline basalt; D=within-plate basalt

某些特定元素和特定指数之间的相关性可以指示了辉绿岩的结晶分异过程(Daietal., 2012)。在以Mg#为横轴，Cr和Ni分别为纵轴的二元图解上，Cr和Ni均表现出和Mg#较强的相关性(图8a, b)，指示了样品经历了镁铁质矿物(如橄榄石、尖晶石和辉石)的分离结晶作用。在以Mg#为横轴，Al2O3和Sr分别为纵轴的二元图解上，Al2O3和Sr与Mg#的负相关性指示了斜长石没有发生明显的分异过程(图8c, d)，这与REE配分图中缺失Eu负异常的特征所反映出的结论一致(图5a)。

狮泉河辉绿岩具有明显的富集不相容元素和REE(图5a, b)，同时其球粒陨石标准化之后的REE配分图和N-MORB标准化之后的微量元素蛛网图显示出和全球平均洋岛玄武岩(OIB)较强的相似性(Fig.5a, b，Sun and McDonough, 1989)。上述特征均指示了该辉绿岩和典型的洋岛玄武岩(OIB)有相同的物质组成(OIB; Hofmann, 1997)。此外，所研究的样品在Zr/Y-Zr图解上显示出较高的Zr/Y比值(图9a)，在TiO2/Yb-Nb/Yb图解上显示出较高的TiO2/Yb比值和Nb/Yb比值(图9b)，结合在2×Nb-Zr/4-Y三元图解和Ti/100-Zr-Y×3三元图解中的投点结果(图9c, d)，本文认为该基性岩具有OIB地球化学特征，即其来源于富集的地幔软流圈(朱弟成等, 2006; Daietal., 2012)。

图10 狮泉河辉绿岩Nb/Th-Ti/Yb 图解数据来源：N-MORB, E-MORB(Sun and McDonough, 1989)；上地壳和下地壳(Rudnick and Fountain, 1995)；峨眉山高Ti玄武岩(Xu et al., 2001)Fig.10 Nb/Th vs. Ti/Yb diagram of the Shiquanhe diabaseData sources: N-MORB, E-MORB(Sun and McDonough, 1989); upper crust and lower crust (Rudnick and Fountain, 1995); Emei high-Ti basalt(Xu et al., 2001)

图11 狮泉河OIB-型镁铁质岩的锆石Ti温度-(Gd/Yb)N散点图(据夏瑛等, 2012; 高晓英和郑永飞, 2011)Fig.11 Scatter plot of Ti-in-zircon temperature vs. (Gd/Yb)N of the OIB-type mafic rocks in the Shiquanhe area(after Xia et al., 2012; Gao and Zheng, 2011)

图12 狮泉河镁铁质岩石的Sm/Yb-La/Yb图(据Zhao and Zhou, 2009)Fig.12 Sm/Yb vs. La/Yb diagrams for the Shiquanhe mafic rocks (melting curves is after Zhao and Zhou，2009)

板内洋岛玄武岩(OIB)通常被认为与“热点”或“地幔柱”有关(牛耀龄, 2010)，其典型的类型有夏威夷洋岛碱性玄武岩和峨眉山高Ti型陆内玄武岩(徐义刚，2002)。此外，在现代俯冲带中也大量发现了具有OIB型地球化学特征的碱性基性岩(Rogersetal., 1985; Reagan and Gill, 1989; Gorring and Kay, 2001)。这些基性岩指示了一些特定的俯冲过程，如：洋中脊俯冲(Abratis and Wörner, 2001)，俯冲终止导致的板片断离(蓝江波等, 2007；Lietal., 2014b, Liuetal., 2014, Zhuetal., 2016)，俯初始俯冲之后弧后拉张成因(Faccennaetal., 2010)。在上述类型中，夏威夷玄武岩由于没有陆壳岩石圈，故其成因上上没有任何岩石圈混染(朱弟成，2006)，峨眉山大火成岩省玄武岩由于存在特殊的地幔柱-岩石圈的相互作用而显示出了高Ti(TiO2=4.25%～4.63%)的地球化学特征(图10，Xuetal., 2001)。上述特征不符合狮泉河OIB型辉绿岩存在地壳混染(Nb负异常、轻微的Ti负异常、强烈Pb正异常)和低含量的Ti含量(TiO2=1.96%～2.27%)的特征。而就位于俯冲带上的基性岩通常与板片流体和板片熔融体交代过的地幔楔中的橄榄岩的部分熔融有关(Wilson, 1989)。这些基性岩以钙碱性系列(成熟的岛弧和陆缘弧)为主，并伴随少量的拉斑玄武岩系列(原始岛弧)(Stern, 2002)，通常显示出HFSE(如Nb和Ta)的亏损和LILE的富集(Wilson, 1989; Stern, 2002)，该特征与本次研究的辉绿岩一致，此外，通过岩浆锆石反映出的岩石结晶的大致温度与俯冲带成因的铁镁质岩的结晶温度区间相似(图11)，故本次研究的辉绿岩的成因应该与俯冲作用密切相关。

俯冲带中就位的OIB型岩浆的起源目前有两种公认的成岩模式：1)被俯冲板片的部分熔融物(埃达克质熔体)交代过的地幔楔源的部分熔融(Defantetal., 1992, Sajonaetal., 1996)；2)洋脊俯冲成因的OIB(Reagan and Gill, 1989; Thorkelsonetal., 2011)。由于俯冲板片熔融物在交代地幔橄榄岩时，会形成富集HFSF的矿物相(Kepezhinskasetal., 1996)，因此，板片熔体交代过的地幔源区派生的玄武岩具有富Nb的特点(Sajonaetal., 1996; Wangetal., 2008)，而本次研究的OIB型基性岩不具有富Nb的特征。所以，俯冲板片对地幔的交代产生OIB的成因不符合本次研究的辉绿岩的特征。此外，由于HREE在石榴子石中具有较强的相容性，所以Sm/Yb比值可以反映地幔橄榄岩源区中石榴子石和尖晶石含量(Rollinson，1993)，同时实验表明La/Yb比值随着部分熔化程度的增加而降低(Aldanmazetal., 2000; Zhuetal., 2008)。从含尖晶石的地幔中分异出的基性岩浆比含石榴子石的地幔中分异出的岩浆有更低的Sm/Yb比值和La/Yb比值，所以这两个比值可以用于判断基性岩源区性质以及地幔部分熔融的深度(Aldanmazetal., 2000)。在Sm/Yb-La/Yb图中(图12)，狮泉河OIB型辉绿岩与洋脊俯冲成因的OIB相比，有更深的源区熔融深度并且显示出了更接近含石榴子石橄榄岩的源区性质(图12)。综上所述，现有的俯冲带相关的OIB成因理论并不能很好的解释狮泉河OIB型辉绿岩的就位机制。

图13 狮泉河俯冲带早期洋壳俯冲演化的地球动力学模型 Fig.13 Geodynamic model for the evolution of the northward subduction of oceanic crust along the Shiquanhe suture at the early stage

狮泉河辉绿岩与上述证明为初始俯冲成因的赞岐岩均侵入于狮泉河蛇绿岩带上的橄榄岩内(图2)，这种就位产状与最近Liuetal. (2018)报道的狮泉河北东侧的160.8±2.3Ma的洋内初始俯冲成因的高镁安山岩的就位产状一致。鉴于两处高镁安山岩处于同一大地构造背景，在野外的分布距离较近并且有着相似的成因机制，推测两者的侵位时间应该均为中侏罗世晚期。因此，同样为中侏罗世晚期侵位的狮泉河OIB型辉绿岩应该与狮泉河高镁安山岩有着相似的初始俯冲成因机制。然而，即使在研究程度较高Izu-Bonin-Mariana (IBM)的初始俯冲相关的弧前序列中也没有相关的OIB型玄武岩(Stern, 2002; Reaganetal., 2013)。但是，同时也注意到在初始俯冲之后的一段时间内在俯冲带的弧后区域可以形成地幔柱相关的OIB(Reaganetal., 2013; Liuetal., 2016)。在这种特殊的OIB成因模式中，俯冲作用与地幔柱相互之间的诱发机制尚未研究清楚。但是Faccennaetal.(2010)通过模拟实验尝试提出了一种解释：俯冲过程开始之后不久由于“回返作用(rollback)”，会在弧后区域形成一个拉张环境(Faccennaetal., 2005)。同时，下行的俯冲板片边缘显示出有活跃的上升物质，这种机制增强了深源地幔和上部物质之间的混合，产生与之相关的碱性OIB型岩浆。这种有着典型地幔亲缘性的OIB岩浆在弧后张性环境中伴随有上部地壳的混染特征而上涌，产生了与俯冲成因的岛弧不同的“异常”岩浆作用(OIB型基性岩)。这种模式的提出解释了现代很多俯冲带附近的与经典岛弧岩浆岩地球化学性质完全不同的OIB型岩浆岩，如地中海直布罗陀俯冲带中的碱性OIB型玄武岩(Duggenetal., 2003; Faccennaetal., 2004)。笔者认为狮泉河OIB型辉绿岩可能为这种成因机制下的产物，理由如下：1)典型的OIB地球化学性质(极度富集的LREE模式以及较高的Zr/Y比值、TiO2/Yb比值和Nb/Yb比值)指示了地幔来源；2)Nb、Ti的负异常和强烈Pb正异表现出了岩浆作用过程中岩石圈的地壳混染；3)狮泉河OIB型辉绿岩的物源形成于中高温-低压环境，与地幔柱(高温)成因以及初始俯冲拉张(低压)的环境一致。因此，上述OIB型辉绿岩可能指示了狮泉河蛇绿岩带在洋内初始俯冲之后一段时间内存在俯冲相关的地幔柱。

4.4 岩石形成的演化过程

本次在狮泉河北约10km发现了具有OIB型的辉绿岩以及高镁安山岩型的闪长岩。闪长岩具有与典型的赞岐岩相似的较高的TiO2含量，Y含量和Yb含量以及较低的Sr/Y比值和La/Yb比值。其U/Th、Nb/Ta和La/Sm指示了源区有俯冲沉积物的贡献，较高的MgO、Cr和Ni含量指示了有地幔物质的参与。其成因需要较高的地温梯度，同时排除洋脊俯冲的成因可能。因此，认为狮泉河高镁安山岩型闪长岩为俯冲早期阶段的弧前热异常环境下熔融的俯冲板片沉积物与上覆地幔楔相互作用的产物。OIB型的辉绿岩的REE配分图和微量元素蛛网图显示出和全球平均洋岛玄武岩(OIB)较强的相似性。Nb负异常、Pb正异常、以及变化的Sr异常指示了岩浆作用过程中的地壳混染以及富集组分部分熔融。Zr/TiO2比值、Nb/Y比值、MgO/SiO2比值、Nb/Zr比值、Nb/Th比值、Ti/Yb比值、La/Sm比值和Sm/Yb比值指示了与现有的俯冲带相关的OIB成因理论不一致的一个俯冲相关的地幔柱成因。结合OIB型辉绿岩中获得的协和度较高的163.7±0.54Ma的岩浆年龄，我们提出了一个中侏罗世晚期班-怒洋在狮泉河地区的早期俯冲模式(图13)：中侏罗世晚期(163.7±0.54Ma)班-怒洋在狮泉河蛇绿岩带的北向初始俯冲过程开始之后不久，俯冲沉积物在浅部发生部分熔融，熔体与上覆地幔楔反应，在弧前区域产生了和日本中新世Setouchi火山岩带的赞岐岩较为相似的高镁安山质岩浆。随后在增生楔下由于弧后拉张作用(或其他机制)产生了与俯冲相关的地幔柱成因的OIB型岩浆。