长白山新生代玄武岩中橄榄岩包体所揭示的岩石圈地幔特征*

徐青鹄 刘嘉麒 莫宣学 贺怀宇 张云辉 赵文斌

1. 中国地质大学(北京)地球科学与资源学院，北京 1000832. 中国科学院地质与地球物理研究所新生代地质与环境研究室，北京 1000293. 西南交通大学地球科学与环境工程学院，成都 611756

华北克拉通作为世界上最古老的克拉通之一(>3.8Ga)，在地质历史时期里经历了巨大的变迁(翟明国, 2019)。在古生代时期，华北克拉通为古老难熔的克拉通型岩石圈地幔；而在晚中生代-新生代时期，受到强烈的构造-热事件影响，华北克拉通被部分改造和破坏，岩石圈地幔厚度减薄并呈现出新生饱满的特征；近年来，华北克拉通岩石圈性质重大转变的研究一直是地球科学研究的焦点问题(郑建平, 1999; 冯亚洲等, 2020; 徐锡伟等, 2019; 朱日祥和徐义刚, 2019)。

已有研究表明，华北克拉通内部出露有大面积的古生代-新生代幔源火山岩，火山岩中含有丰富的橄榄岩、辉石岩等地幔包体(郑建平和戴宏坤, 2018)。地幔包体作为岩石圈地幔最直接的天然样品，可以揭示深部岩石圈地幔的部分熔融与地幔交代作用，是研究岩石圈地幔的性质、演化和岩石圈地幔减薄以及活化的重要媒介(Griffinetal., 2009)。许多学者对华北克拉通内部宽甸、汉诺坝和昌乐等幔源包体进行了详细的岩石学、地球化学及同位素研究，结合地球物理解译，发现华北克拉通岩石圈地幔属性和深部过程在空间上存在明显的不均一性，其东部岩石圈地幔的改造作用较为强烈，表现为新生饱满的特征；而中-西部则仍然整体保持着古老难熔的克拉通属性，仅局部岩石圈受到了改造作用(Yangetal., 2018; Dai and Zheng, 2019; Haoetal., 2019; Wangetal., 2019; Zouetal., 2020)。此外，制约华北克拉通不同部位改造作用的构造动力学机制也不尽相同：华北克拉通东部的岩石圈地幔可能受到了太平洋俯冲的影响(Lietal., 2016)，北部的岩石圈地幔可能受到了中亚造山带的影响(Daietal., 2019)，而其东南部则可能受到了华南克拉通与华北克拉通碰撞的影响(Yingetal., 2006)。

华北克拉通东北部出露了以长白山地区为代表的大范围新生代火山岩，主要包括天池和龙岗火山岩群。在天池和龙岗火山岩群均已发现丰富的地幔橄榄岩包体，为研究长白山地区岩石圈地幔的性质与演化提供了优越条件。前人对天池和龙岗火山岩群开展了大量的研究，在火山岩形成时代、成因机制和深部结构方面取得了丰硕的成果(樊祺诚, 2008; 郭文峰等, 2015; 陈棋福等, 2019)；但是，针对火山岩中橄榄岩包体的研究仍较为匮乏，成为制约学者们揭示长白山地区乃至华北克拉通东北部下覆岩石圈地幔特征的瓶颈。前人对天池和龙岗火山岩群中地幔橄榄岩包体进行了一定研究(Tangetal., 2012; Wangetal., 2013, 2017; Xuetal., 2019)；然而，由于样品数量有限和测试方法不足，长白山地区岩石圈地幔属性和深部过程(部分熔融和地幔交代)仍然未能得到准确限定，因此很有必要对长白山地区天池和龙岗火山岩群的橄榄岩包体开展详细、系统的地球化学研究。同时，华北克拉通东北部岩石圈地幔的转化过程及其构造动力学机制亦存在不同认识，诸如，转化过程包括：热力学-化学破坏作用(Xu, 2001)、古老地壳的拆沉作用(Gaoetal., 2008)和古大洋俯冲引起的熔体-橄榄岩相互作用(Linetal., 2019)；构造动力学机制可能为古亚洲洋的向南俯冲(Copeetal., 2005)或太平洋板块的向西俯冲(朱日祥和徐义刚, 2019)。地幔橄榄岩包体的地球化学证据将可以在解决这一争议中发挥重要作用。

因此，本文在前人研究基础上，广泛采集了长白山地区的天池和龙岗火山岩群中地幔橄榄岩包体样品，对橄榄岩包体的岩石学、全岩主微量元素、单矿物主微量元素、单斜辉石Sr-Nd-Hf同位素进行了详细、系统地分析，在此基础上查明长白山地区下覆岩石圈地幔属性，探讨岩石圈地幔深部的部分熔融和地幔交代作用，揭示制约华北克拉通东北部岩石圈地幔改造作用的构造动力学机制。

1 区域地质背景和橄榄岩包体特征

长白山地区位于我国东北部的吉林省，区域构造上属华北克拉通东北缘，构造格架主要由NEE向与NNE向的断裂带构成，其次为EW向、NS向及NW向的断裂(图1)。长白山地区地处胶辽断块的隆起部位。胶辽断块的基底岩石从第四纪早期时开始缓慢隆起上升，此后从中更新世至今呈现出非持续性的上升，导致该区地貌出现严重变形，为玄武岩喷发创造了条件，最终在长白山地区形成了以天池和龙岗火山群为代表的大范围火山岩群(樊祺诚, 2008)。长白山地区火山岩类型以造盾期的玄武岩、造锥期粗面岩和晚期碱流岩为主，近代喷发伴有火山碎屑流的堆积(刘嘉麒等, 2015)。火山群中常见有富含地幔包体的新生代碱性玄武岩。

图1 华北克拉通东北缘区域地质图(a)和研究区玄武质岩石分布和采样点(b)F1-塔源-喜贵图断裂；F2-贺根山-嫩江断裂；F3-西拉木伦断裂；F4-开原-赤峰断裂；F5-伊兰-伊通断裂；F6-牡丹江断裂；F7-敦化-密山断裂Fig.1 Regional map of the northeastern North China Craton (a) and basaltic rock distribution and sampling locations in study area (b)F1-Tayuan-Xiguitu fault; F2-Hegenshan-Nenjiang fault; F3-Xilamulun fault; F4-Kaiyuan-Chifeng fault; F5-Yilan-Yitong fault; F6-Mudanjiang fault; F7-Dunhua-Mishan fault

天池火山群地处中国和朝鲜交界线附近，是中国现有的休眠活火山之一。天池火山喷发从新生代开始至今一直间断性发生，其喷发历史被划分为4个阶段：(1)渐新世到中新世，形成了大面积的碱性玄武岩；(2)上新世到更新世，形成了大型粗面玄武岩的熔岩盾；(3)约1Ma时期，主要为粗面质火山岩的喷发；(4)全新世，大规模火山强烈爆发活动形成了浮岩质火山碎屑堆积物和天池的破火山口(刘嘉麒等, 2015)。龙岗火山群坐落在吉林辉南和靖宇县辖区的龙岗山脉之中，其火山活动开始于新近纪早期，第四纪达到高峰，龙岗火山活动可分为3个阶段：(1)小浪山期；(2)中-晚更新世龙岗期，为产生火山渣和低平火山口的主要时期；(3)约一千五百年前至一千六百年前的金龙顶子期，产生火山碎屑(白志达等, 2006)。

长白山地区地幔橄榄岩包体采集于天池和龙岗火山群玄武岩中，共计23件。包体多呈棱角状至椭圆状，表面未见蛇纹石化特征，样品新鲜且呈翠绿-墨绿色，包体直径大小为0.5～20cm(大多约12cm)。包体岩性主要为尖晶石二辉橄榄岩，少量为方辉橄榄岩，具有原生粒状结构或等粒变晶结构(图2)，未发现金云母和石榴石。二辉橄榄岩主要矿物组合和含量为：橄榄石(55%～77%)+斜方辉石(9%～29%)+单斜辉石(3%～22%)+尖晶石(0.1%～3%)；方辉橄榄岩具有较高的橄榄石和较低的斜方辉石、单斜辉石含量，其矿物含量为：橄榄石(86%～89%)+斜方辉石(8%～12%)+单斜辉石(1%～2%)+尖晶石(～1%)(表1)。前人研究中曾在长白山地区橄榄岩包体中发现有角闪石(史兰斌等, 1999)，但是在本次研究未观察到。镜下二辉橄榄岩和方辉橄榄岩中的橄榄石和斜方辉石大小约0.5～3mm，为自形-半自形晶形。橄榄石发育有明显的裂纹；单斜辉石一般呈翠绿色，大小约0.1～1.5mm，为他形晶形；尖晶石多为褐色，经常出现在颗粒之间。此外，镜下观察样品发育有较多的熔体囊，可见比较明显的熔体交代特征。

表1 橄榄岩包体矿物组分和平衡温度表

2 测试方法

2.1 全岩主量和微量元素分析

橄榄岩样品的全岩主量和微量元素测试分析实验在中国科学院广州地球化学研究所完成。全岩主量元素测试工作在Rigaku ZSX-100e X荧光光谱仪上完成，分析精度优于±3%。全岩微量元素采用PE Elan 6000型电感耦合等离子体质谱议(ICPMS)完成，分析精度和准确度优于5%。

2.2 单矿物主量元素和微量元素分析

单矿物主量元素和微量元素分析在武汉上谱测试公司完成。主量元素测试所采用的仪器为JEOL-JXA-8100电子探针分析仪。分析矿物主要包括单斜辉石、斜方辉石、橄榄石和尖晶石，相同矿物在每个橄榄岩包体样品中选取5个分别进行测试，取其平均值作为分析数据。单斜辉石、斜方辉石和尖晶石的分析条件为加速电压15kV，探针电流20nA。橄榄石分析条件为加速电压为20kV，探针电流20nA，同时延长峰值时间，其中Mn、Ni、Cr、Ca和Al测试时间为100s。主量元素的分析精度优于±0.5%。

单矿物原位微区微量元素含量在武汉上谱分析科技有限责任公司利用LA-ICP-MS完成。GeolasPro激光剥蚀系统由COMPexPro 102ArF 193nm准分子激光器和MicroLas光学系统组成，ICP-MS型号为Agilent 7700e。激光剥蚀过程中采用氦气作载气、氩气为补偿气以调节灵敏度，二者在进入ICP之前通过一个T型接头混合，激光剥蚀系统配置有信号平滑装置(Huetal., 2015)。单矿物微量元素含量处理中采用玻璃标准物质BHVO-2G，BCR-2G和BIR-1G进行多外标无内标校正。每个时间分辨分析数据包括大约20～30s空白信号和50s样品信号。详细的仪器参数和分析流程见(Zongetal., 2017)。对分析数据的离线处理采用软件ICPMSDataCal(Liuetal., 2008)完成。每个橄榄岩包体样品的单斜辉石选取5个点进行测试，取其平均值开展分析。

2.3 单斜辉石Sr-Nd-Hf同位素分析

单斜辉石Sr-Nd-Hf同位素的测定是在桂林理工大学地球科学学院MC-ICP-MS同位素实验室完成的。从橄榄岩样品中挑选出足量单斜辉石颗粒，研磨至粉末以备分析。采用Neptune Plus多接受电感耦合等离子质谱仪(MC-ICP-MS)测试Sr-Nd-Hf同位素。Sr同位素和Nd同位素的测量监控采用不同的标准，分别为国际标样NBS987、实验室标准Sr-GIG和国际标准Jndi-1、实验室标准Nd-GIG，数据的标准化值分别采用：87Sr/86Sr=0.1194、143Nd/144Nd=0.7219。标样NBS987的87Sr/86Sr比值和Jndi-1的143Nd/144Nd比值分别为0.710288±28(σ)和0.512109±12(2σ)。Sr-Nd同位素分析精度好于0.002%。利用Neptune型MC-ICP-MS仪器分别对Lu、Hf进行质谱测试，Yb用于校正出Lu的分馏，179Hf/177Hf=0.7325和Alfa Hf标准的176Hf/177Hf=0.282185可以对Hf的分馏校正。详细流程见文献(韦刚健等, 2002; 梁细荣等, 2003)。

P86兼有辅助捕收、起泡、消泡等多项特性，合理使用P86对减少捕收剂TL-1的用量、调整泡沫现象和改善操作有较大关系，合理的P86用量会对锡石浮选产生良好的效果。

3 测试结果

3.1 全岩主量元素和微量元素特征

长白山地区橄榄岩包体的全岩主量和微量元素测试结果列于表2。长白山大部分样品由于Fe2O3T含量高，使得烧失量(LOI)的结果为负值。与方辉橄榄岩相比，二辉橄榄岩含有更高的Al2O3(2.32%～3.72%)和更低的MgO(38.2%～41.9%)。橄榄岩包体的主量元素含量低于原始上地幔，二辉橄榄岩具有较低的MgO含量，表现出新生饱满的特征，少数的方辉橄榄岩表现出新生饱满向古老难熔过渡的特征(图3)。

图3 橄榄岩包体全岩主量元素与微量元素二元图解原始上地幔(PUM)引自Hart and Zindler (1986)；亏损地幔和大洋橄榄岩趋势线引自Boyd (1989)和Griffin et al. (1998b)；熔体抽取趋势线引自Frey et al. (1985)和Xu et al. (2008)Fig.3 Binary diagrams of major element and trace element compositions of peridotite xenolithsPrimitive upper mantle from Hart and Zindler (1986); depleted mantle and oceanic trend line from the Boyd (1989) and Griffin et al. (1998b); melting extraction trend from Frey et al. (1985) and Xu et al. (2008)

表2 橄榄岩包体全岩主量元素(wt%)和微量元素(×10-6)分析测试结果

3.2 单矿物主量元素特征

本次研究对长白山地区橄榄岩包体中橄榄石、斜方辉石、单斜辉石和尖晶石的主量元素含量开展电子探针分析，测试结果如表3所示。主微量元素二元图解表明长白山地区方辉橄榄岩具有过渡型岩石圈地幔特征，二辉橄榄岩表现出新生饱满的岩石圈地幔特征(图4)。长白山地区二辉橄榄岩包体中橄榄石MgO含量为45.32%～49.44%，Fo值(Fo=100×Mg2+/(Mg2++Fe2+))为85.6～91.0，FeO含量为8.84%～13.74%。与此对比，方辉橄榄岩包体中橄榄石具有更高的MgO(49.01%～49.70%)、Fo值(90.8～91.3)和更低的FeO(8.47%～8.90%)。二辉橄榄岩包体中斜方辉石的Al2O3含量为2.95%～4.77%，Mg#值为86.2～91.4；方辉橄榄岩包体中斜方辉石的Al2O3含量为2.94%～3.31%，Mg#值在91.5～91.9之间。方辉橄榄岩包体比二辉橄榄岩具有更高的斜方辉石Cr#(Cr#=100×molar Cr/(Cr+Al))值(10.8>4.1～8.8)，和Mg#值呈现较好的正相关关系。二辉橄榄岩包体的单斜辉石Al2O3分别为5.03%～6.37%，Mg#值变化较大为85.9～92.1，Cr#值在5.7～10.9之间。方辉橄榄岩包体的单斜辉石Al2O3、Mg#值和Cr#值含量分别为3.85%～4.00%、92.1～92.4和15.8。二辉橄榄岩包体的尖晶石Cr#值为8.9～18.8，Mg#值为66.9～76.4，TiO2含量为0.01%～0.24%；方辉橄榄岩包体的尖晶石具有更高的Cr#值(33.2～34.8)，Mg#值(71.0～71.7)，TiO2含量为0.09%～0.16%。

图4 矿物主量元素关系图(a)橄榄石Mg#值-CaO含量，(b)斜方辉石Cr#值-Mg#值，(c)单斜辉石Mg#值-Al2O3含量，(d)尖晶石Cr#值-Mg#值.山旺数据来源于Zheng et al. (1998)；鹤壁数据来源于Zheng et al. (2001)Fig.4 Binary diagrams of major element and trace element compositions of minerals in peridotite xenoliths(a) olivine Mg# values vs. CaO contents, (b) orthopyroxene Cr# values vs. Mg# values, (c) clinopyroxene Mg# values vs. Al2O3 contents, (d) spinel Cr# values vs. Mg# values. Data of Shanwang area from Zheng et al. (1998); Data of Hebi area from Zheng et al. (2001)

3.3 单斜辉石微量元素特征

长白山地区橄榄岩包体的单斜辉石微量元素含量如表4所示。本次研究根据球粒陨石标准化稀土元素配分曲线特征，进一步将单斜辉石划分为以下四个类型(图5)：(1)LREE严重亏损型：单斜辉石的轻重稀土元素分异较为明显，(La/Yb)N值最低，为0.01～0.02，总体呈现LREE严重亏损和HREE平坦的特征，反映了二辉橄榄岩可能在尖晶石稳定区经历了不同程度熔融抽取后的残留，熔融比例在1.5%；(2)LREE亏损型：单斜辉石的轻重稀土元素分异略明显，(La/Yb)N值较低，为0.07～0.33，总体呈现LREE亏损和HREE平坦的特征。(3)向右微倾型：单斜辉石轻稀土元素相对重稀土元素富集的特征，(La/Yb)N值为1.04～1.14，整体表现为平缓向右倾斜的直线；(4)“勺型”：单斜辉石LREE相对于MREE略微富集，(La/Yb)N值为0.95～1.25。橄榄岩包体中单斜辉石的微量元素表现出相近的特征，均富集大离子亲石元素(Rb、Ba、Th、U)，亏损高场强元素(Nb、Ta、Zr、Hf、Ti)和Pb。

图5 长白山地区橄榄岩包体中单斜辉石球粒陨石标准化稀土元素配分图和原始地幔标准化微量元素蛛网图解(标准化值引自Sun and McDonough, 1989)Fig.5 Chondrite-normalized REE patterns and spider diagrams of clinopyroxenes in peridotite xenoliths from the Changbaishan (normalization values after Sun and McDonough, 1989)

3.4 平衡温度

本次研究的橄榄岩包体主要为尖晶石相二辉橄榄岩，不含斜长石和石榴子石矿物。尖晶石相橄榄岩通常稳定存在于0.8～2.5GPa环境中(Nickel, 1986; Webb and Wood, 1986)，其形成时的准确压力值难以用压力计计算，本次研究选取前人在华北克拉通东北部岩石圈地幔研究所用的压力值1.5GPa为参考值(Linetal., 2019)。为了避免使用单个地质温度计造成的系统误差，本次研究采用二辉石温度计(TWells)(Wells, 1977)和斜方辉石Ca温度计(TBK)(Brey and Köhler, 1990)共同估算长白山和龙岗火山岩中橄榄岩包体的平衡温度。当Mg#值比较高(>80)的橄榄岩包体，平衡温度范围在900～1400℃时，选用二辉石温度计能够给出更为准确的温度。因此，在本次研究的计算结果中，优先考虑二辉石温度计的平衡温度，与斜方辉石Ca温度计的结果进行对比，来确保结果的可靠性。两种地质温度计给出的结果见图6和表1。二辉石温度计和斜方辉石Ca温度计的计算结果较为一致， 分别为928～1101℃和947～1139℃，说明长白山地区橄榄岩包体的平衡温度范围大约为900～1100℃。

图6 橄榄岩包体平衡温度频率直方图(a)二辉石温度计(TWells)；(b)斜方辉石Ca温度计(TBK)Fig.6 Frequency histograms of equilibrium temperatures(a) TWells, Na-in-orthopyroxene-clinopyroxene; (b) TBK, Ca-in-orthopyroxene geothermometer

3.5 单斜辉石Sr-Nd-Hf同位素特征

长白山地区橄榄岩包体中单斜辉石的Sr、Nd和Hf同位素分析结果如表5所示。87Sr/86Sr=0.702749～0.707276、143Nd/144Nd=0.512886～0.51333(εNd=+4.8～+13.6)、176Hf/177Hf=0.283273～0.284179(εHf=+17.7～+49.8)。由图7可知，研究区二辉橄榄岩包体中单斜辉石的87Sr/86Sr和143Nd/144Nd比值变化范围较大。整体上可以看出，长白山地幔橄榄岩样品中单斜辉石表现出整体Sr同位素亏损，部分样品Sr同位素富集的特征，样品Nd-Hf同位素整体亏损的特征，整体基本与中国东部新生岩石圈地幔的Sr-Nd-Hf同位素特征相似。部分橄榄岩样品沿地幔演化线发生了明显的偏离，表现出不同程度的Nd-Hf同位素解耦。

表5 橄榄岩包体单斜辉石Sr-Nd-Hf同位素特征

图7 橄榄岩包体中单斜辉石143Nd/144Nd-87Sr/86Sr(a)和176Hf/177Hf-143Nd/144Nd(b)同位素比值图解OIB-洋岛玄武岩；MORB-洋中脊玄武岩；BSE-整体硅酸岩地球. 洋中脊玄武岩数据和地幔演化线引自Zindler and Hart (1986)；洋岛玄武岩数据引自Chauvel and Blichert-Toft (2001)和Salters and White (1998)；整体硅酸岩地球数据引自Bouvier et al. (2008)和Blichert-Toft et al. (1999)；长白山玄武岩数据引自Basu et al. (1991)、Guo et al. (2014)、Kuritani et al. (2009)和Zhang et al. (2014)；三个交代作用的混合端元包括：长白山玄武岩(Nd=35.9×10-6、Hf=5.1×10-6, Kuritani et al., 2009)、大洋岛弧玄武岩(Nd=10.1×10-6、Hf=1.7×10-6，Kelemen et al., 2005)、富水的熔/流体(Nd=11.0×10-6、Hf=1.0×10-6, Kessel et al., 2005)Fig.7 143Nd/144Nd vs. 87Sr/86Sr (a) and 176Hf/177Hf vs. 143Nd/144Nd (b) isotopic ratios diagrams of peridotite xenolithsOIB-oceanic-island basalt; MORB-mid-ocean-ridge basalt; BSE-bulk silicate earth. Data for MORB and mantle array from Zindler and Hart (1986); data for OIB from Chauvel and Blichert-Toft (2001), Salters and White (1998); Data for BSE are from Bouvier et al. (2008) and Blichert-Toft et al. (1999); basaltic rocks from Changbaishan from Basu et al. (1991), Guo et al. (2014), Kuritani et al. (2009), and Zhang et al. (2014); three mixing end members include: Changbaishan basalt (Nd=35.9×10-6、Hf=5.1×10-6, Kuritani et al., 2009), oceanic arc basalts (Nd=10.1×10-6、Hf=1.7×10-6，Kelemen et al., 2005), aqueous fluids (Nd=11.0×10-6、Hf=1.0×10-6, Kessel et al., 2005)

4 讨论

4.1 岩石圈地幔的部分熔融过程

长白山地区橄榄岩包体岩石类型以二辉橄榄岩为主，含少量方辉橄榄岩，表明长白山地区的岩石圈地幔组成主要为富集和中等亏损型地幔。长白山地区包体橄榄岩的MgO含量和其它主微量元素含量表现出明显负相关的线性关系，与世界其它地区地幔橄榄岩类似，说明长白山地区橄榄岩包体是经历不同程度的熔融抽取后的结果(Freyetal., 1985; Xuetal., 2008)。地幔橄榄岩中橄榄石的Fo值、尖晶石的Cr#值和TiO2含量可以有效指示地幔岩石熔体抽取的程度，随着部分熔融程度的增加，橄榄石的Fo值和尖晶石中的Cr#值都会增大(Hellebrandetal., 2001)。长白山地区橄榄岩包体中橄榄石的Fo值(85.6～91.3)和尖晶石的Cr#(8.9～34.8)变化较大，较低的橄榄石Fo值可排除典型的古老难熔克拉通型的地幔的可能性(Fo>92)。由图8可知，长白山-龙岗二辉橄榄岩经历的部分熔融程度均小于10%。其中，2件方辉橄榄岩包体样品比二辉橄榄岩包体样品具有较高的Cr#值，可能代表了部分熔融程度更高的地幔。此外，尖晶石沿着拟合熔融曲线有不同程度的偏移，说明除了部分熔融作用以外可能还发生了其它的地质过程。方辉橄榄岩包体中的尖晶石在同样Cr#值的情况下更富集TiO2，这个偏离说明方辉橄榄岩在熔融残留的后期可能与更富集Ti的熔体发生了交代作用；橄榄岩包体成分点落入深海橄榄岩区域并指向更为富集的趋势，Cr2O3含量和Cr#值的变化较大也可能是由于后期交代作用所导致的。

图8 橄榄岩包体中尖晶石Cr#-TiO2图(a,据Park et al., 2017)和橄榄石Fo值-尖晶石Cr#值图(b)(a)点线为深海橄榄岩中尖晶石成分区域(Dick and Bullen, 1984)，灰色箭头代表了橄榄岩与似MORB熔体交代反应的变化趋势，图中虚线是熔融趋势线，引自Choi et al. (2008)；(b)地幔演化线和橄榄石-尖晶石熔融趋势引自Arai (1994)；超俯冲型橄榄岩和深海橄榄岩数据引自Ishii et al. (1992), Pearce et al. (2000)Fig.8 Diagrams of Spinel Cr#-TiO2 (a, modified after Park et al., 2017) and spinel Cr#-Fo contents in the coexisting olivine (b) of peridotite xenoliths(a) the dotted curve is a field of abyssal peridotite spinel (Dick and Bullen, 1984). Gray arrow represents an effect of metasomatic reaction with MORB-like melt. Melting trend, shown in dashed curve, is from Choi et al. (2008). (b) Arai (1994) was the source of mantle array and melting trend values for olivine-spinel. Data on supra-subduction zone (SSZ) peridotites and abyssal peridotites sourced from Pearce et al. (2000) and Ishii et al. (1992)

由于单斜辉石中不相容的重稀土元素含量不易被后期地幔交代作用改变，因此可以利用单斜辉石中的微量元素Yb和Y含量模拟岩石圈地幔的部分熔融程度(Norman, 1998)。本次研究采用分离部分熔融模型进行估算，结果表明长白山地区的二辉橄榄岩部分熔融程度介于3%～5%之间，方辉橄榄岩的部分熔融程度较大，为10%～12%(图9)。Hellebrandetal. (2001)认为当Cr#值在10～60范围时，可利用公式F=10ln(Cr#)+24，计算橄榄岩的部分熔融程度F。计算结果表明，长白山二辉橄榄岩的部分熔融程度(0.8%～7.3%)低于方辉橄榄岩的部分熔融程度(13.0%～13.4%)，这也与图8、图9中的部分熔融估算结果一致。

图9 橄榄岩分离熔融模拟图(分离熔融模拟曲线据Norman, 1998;标准化值据McDonough and Sun, 1995)Fig.9 Fractional partial melting modeling for the Changbaishan and Longgang peridotite (fractional partial melting model after Norman, 1998; normalized to primitive mantle after McDonough and Sun, 1995)

4.2 岩石圈地幔的交代富集作用

岩石圈地幔在受到部分熔融抽取后往往还会受到后期熔体/流体交代作用的改造，交代作用通常分为两种形式：显性交代和隐性交代。当发生显性交代时，会形成新的矿物如角闪石、磷灰石或金云母(Gordon Medarisetal., 2015)，隐性交代通常会导致不相容元素的富集而不会形成其它矿物相(Uenver-Thieleetal., 2017)。长白山地区橄榄岩包体中缺乏角闪石、磷灰石或金云母等显性交代的矿物，在全岩和单斜辉石微量元素上具有明显的富集轻稀土元素的特征，并且表现出Nd-Hf同位素解耦现象，说明长白山地区橄榄岩在形成时发生了隐性交代作用。

已有研究表明，岩石圈地幔的交代介质主要来源有碳酸盐熔体(Rudnicketal., 1993)、富水硅酸盐熔体(Zanganaetal., 1999)和富水流体(Xuetal., 2003)。在长白山地区橄榄岩中第一、四组单斜辉石具有亏损Nb的特征，第二、三组的单斜辉石Nb含量相对高一些，暗示在地幔交代中可能有Nb的加入(图5)。由于富水流体并不能传输Nb(Eggler, 1987)，因此排除了交代介质为富水流体的可能。组一和组二的样品中轻重稀土元素分异明显，并且略微富集大离子亲石元素，大多伴有Zr的亏损，说明长白山地区橄榄岩包体的形成受到熔体交代影响(Blusztajn and Shimizu, 1994)。进而，需要进一步辨别交代介质是碳酸盐还是硅酸盐熔体。交代反应中通常用单斜辉石(La/Yb)N和Ti/Eu的关系图来指示熔体的性质(Coltortietal., 1999)。高的(La/Yb)N和低的Ti/Eu比值多体现碳酸盐熔体的交代作用，一般与稳定大陆的岩浆作用有关；而低的(La/Yb)N和高的Ti/Eu比值多与硅酸盐熔体交代有关(郑建平, 2009)。长白山地区橄榄岩包体中的单斜辉石Ti/Eu比值较高且变化较大，而(La/Yb)N比值相对很低，可以明显反映出硅酸盐熔体的交代作用(图10)。此外，碳酸盐熔体形成的橄榄岩通常富集轻稀土元素、Sr、Ba和Rb元素，亏损高场强元素，具有较高的 (La/Yb)N值(Suetal., 2012)，与长白山地区橄榄岩包体的微量元素含量明显不同；橄榄岩包体中还缺乏碳酸盐交代的典型矿物，如方解石、磷灰石等(Yaxleyetal., 1998)，排除了交代介质来源于碳酸盐熔体的可能性。

图10 橄榄岩中单斜辉石的(La/Yb)N-Ti/Eu关系图(据Coltorti et al., 1999)粒陨石标准化值据Sun and McDonough (1989)Fig.10 Plot of (La/Yb)N vs. Ti/Eu of clinopyroxene from peridotite (after Coltorti et al., 1999)Chondrite-normalized values after Sun and McDonough (1989)

橄榄岩中单斜辉石的Nd-Hf同位素也可以约束交代熔体的性质。在Nd-Hf同位素关系图中(图7b)，研究区的橄榄岩包体落在Nd-Hf同位素拟合的地幔演化线上。这样的Nd-Hf同位素解耦在地幔橄榄岩中经常观察到(Choi and Mukasa, 2012; Liuetal., 2017)。为查明交代介质的性质，我们用亏损地幔(143Nd/144Nd=0.5153、176Hf/177Hf=0.2853、Nd=1×10-6、Hf=0.4×10-6)和三个交代介质分别模拟出混合曲线，包括长白山玄武岩(Nd=35.9×10-6、Hf=5.1×10-6, Kuritanietal., 2009)、岛弧玄武岩(Nd=10.1×10-6、Hf=1.7×10-6, Kelemenetal., 2005)、富水的熔/流体(Nd=11.0×10-6、Hf=1.0×10-6, Kesseletal., 2005)。由图7b可知，长白山和龙岗二辉橄榄岩样品主要落在亏损地幔与富水流体端元和岛弧玄武岩端元混合的曲线上。综上所述，长白山地区二辉橄榄岩样品经历的交代作用中交代介质为来源于与俯冲作用相关的富水硅酸盐熔体。

4.3 岩石圈地幔的形成时代

长白山地区地幔包体均为不含石榴子石的尖晶石橄榄岩，表明其来自低于80km深度的地幔，全岩和矿物主量元素显示形成尖晶石橄榄岩的地幔具有饱满和高温的特征。Fo值大于91的橄榄石被认为是太古代-元古代岩石圈地幔碎片的标志，Fo值小于90的橄榄岩被认为是显生宙新增生的饱满岩石圈地幔(Griffinetal., 1998a; Zhengetal., 2007)。长白山地区地幔橄榄岩包体中橄榄石的Fo值变化范围较大，显示出新增生的岩石圈地幔和过渡型岩石圈地幔的特征。在图11中可见，长白山地区二辉橄榄岩包体大多落在显生宙区域，而方辉橄榄岩包体则落在元古宙区域附近，样品整体沿着大洋演化趋势线分布，说明在长白山地区下部存在不同时代的岩石圈地幔。

图11 橄榄岩中橄榄石含量与橄榄石Fo值关系图解大洋趋势演化线来自Boyd (1989).太古代、元古代和显生宙区域范围来自Griffin et al. (1998a)Fig.11 Plot of modal olivine in peridotite versus Fo in olivineOceanic trend derived from Boyd (1989). Fields for the Archean, Proterozoic, and Phanerozoic from Griffin et al. (1998a)

在图12中，长白山地区橄榄岩包体中橄榄石的Fo值与平衡温度之间不存在相关性，说明长白山地区存在新生饱满和过渡型岩石圈地幔，但是二者在垂向上并未表现出明显的分层效应。Xuetal. (2019)对研究区橄榄岩含水量的研究发现，全岩水含量超过50×10-6(最高达132×10-6)的样品代表新增生的软流圈来源的岩石圈地幔，而全岩水含量小于50×10-6(低至9×10-6)的样品代表着减薄的古老的岩石圈地幔，这也与前人研究的华北克拉通地区新生代玄武岩所携带的地幔橄榄岩水含量相一致(Xiaetal., 2010)。

图12 地幔橄榄岩橄榄石Fo值与平衡温度关系图解Fig.12 Equilibration temperature versus Fo in olivine for the peridotite

长白山地区橄榄岩包体的单斜辉石表现出整体Sr同位素亏损，部分样品Sr同位素富集的特征，以及整体Nd-Hf同位素亏损的特征。在图13中，长白山地区橄榄岩包体的Sm-Nd和Lu-Hf同位素特征表现出一定的正相关关系，在橄榄岩包体中选取不受熔体交代及受熔体交代作用影响弱的样品，进行等时线拟合分析，分别得到Sm-Nd和Lu-Hf同位素模式年龄127～510Ma和260～853Ma，二者的年龄范围具有较大的差异，Nd-Hf同位素发生了一定的解耦现象(图7)，说明橄榄岩包体的Nd-Hf同位素可能受到对流地幔分离后岩石圈晚期交代作用改造(Bizimisetal., 2007)。研究表明Nd比Hf在含水溶液或者熔体中更易迁移(Ferrandoetal., 2019)。因此，相比于Lu-Hf同位素体系，熔/流体交代作用更容易重置Sm-Nd同位素体系。同时，Lu-Hf同位素体系具有较高的封闭温度，在岩石圈地幔形成后不容易受到后期扰动，可以得到较为可信的年龄。此外，我们挑选出不受熔体交代作用以及受熔体交代弱的橄榄岩样品。这些不受熔体交代作用或者交代作用不明显的样品的单斜辉石Hf同位素的模式年龄能够提供一个比较可靠的岩石圈地幔年龄，这些样品的单斜辉石Hf同位素计算出模式年龄的结果为260～853Ma，橄榄岩包体的形成时间范围与华北克拉通的演化历史相符(翟明国, 2019)。以上矿物学、地球化学和同位素证据表明，长白山地区下覆岩石圈地幔以新增生的饱满地幔为主，夹有少量古老难熔的岩石圈地幔碎片。

图13 长白山地区橄榄岩包体单斜辉石147Sm/144Nd-143Nd/144Nd关系图解(a)和176Lu/177Hf-176Hf/177Hf关系图解(b)亏损地幔(DM)和球粒陨石均一储库(CHUR)数据来自Blichert-Toft and Albarède (1997), Pearson and Nowell (2004)和Bouvier et al. (2008). 计算Lu-Hf等时线时衰变常数λ=1.865×10-11y-1(Scherer et al., 2001)Fig.13 Plot of 147Sm/144Nd vs. 143Nd/144Nd (a) and plot of 176Lu/177Hf vs. 176Hf/177Hf (b) for clinopyroxenes from Changbaishan peridotitesData for depleted mantle (DM) and chondritic uniform reservoir (CHUR) from Blichert-Toft and Albarède (1997), Pearson and Nowell (2004) and Bouvier et al. (2008). Calculation of the Lu-Hf errorchron was based on the decay constant λ=1.865×10-11y-1 (Scherer et al., 2001)

4.4 岩石圈地幔转化的构造动力学机制

在华北克拉通内，古生代金伯利岩的包体指示厚、古老和难熔的前寒武系岩石圈地幔，新生代玄武岩中橄榄岩包体指示薄、新生和饱满的新生代岩石圈地幔(郑建平, 2009)。本次研究发现，长白山地区以新生饱满岩石圈地幔为主，含少量古老难熔岩石圈地幔，说明在长白山地区所处的华北克拉通东北部发生了明显的岩石圈地幔的转化。华北克拉通岩石圈地幔的转化过程一直是学者研究的焦点，其构造动力学机制仍存在不同认识，包括热化学侵蚀(Xu, 2001)、古老地壳的拆沉作用(Gaoetal., 2008)和板片俯冲作用(Linetal., 2019)。

研究表明，华北克拉通东部岩石圈地幔具有“东薄西厚”的特征，是太平洋板块向西俯冲的结果(Zhengetal., 2009)。长白山地区位于华北克拉通东北部，地球物理(Zhao and Ohtani, 2009)和Os同位素(Zhangetal., 2017)证据表明在太平样板片形成的“地幔楔”越过长白山地区，向东俯冲至华北克拉通东部的郯庐断裂带下部。构造地貌学和构造应力场研究发现，在华北克拉通东部的盆地、断裂和变质核杂岩都表现出NE-NNE的线性展布规律(Zhuetal., 2012)，也与太平洋板块向西俯冲的构造动力学机制一致。此外，华北克拉通东部侏罗-白垩纪岩浆活动也证明是太平洋板块向西的俯冲挤压和伸展拉伸阶段的产物(Liangetal., 2020)。综上所述，太平洋板块的俯冲作用对华北克拉通东北部岩石圈地幔造成的影响有：华北克拉通东北部岩石圈的拉张和伸展的时空变化与俯冲的太平洋板块的方向和角度相一致；俯冲太平洋板块为中-新生代岩浆提供了拉张伸展的环境且提供了源区物质。

地球物理的证据表明东北地区地幔过渡带深度存在着滞留的太平洋板块(Zhao and Ohtani, 2009)研究表明地幔过渡带是地幔内部一个较大的含水储库(Pearsonetal., 2014)。地幔过渡带中水的长期循环流动过程中会导致局部的瑞利-泰勒不稳定，进而产生上涌的水柱(Richard and Bercovici, 2009)。层析图像的结果也表明在长白山地区从上地幔至410km(地幔过渡带的最上端)存在着一个显著的低速异常(Huang and Zhao, 2006; Zhaoetal., 2009)。长白山地区上地幔深度存在的低速异常很可能是由于来源于地幔过渡带的水柱所致(Zhaoetal., 2009)。这些水柱在上升过程中，将先前存在地幔过渡带的太平洋板块的物质(洋壳或和沉积物)运送到上地幔。这些再循环的太平洋板块的物质在上地幔的某一深度发生部分熔融产生熔体，这些熔体在上升过程中与长白山地区古老难熔的岩石圈地幔发生相互作用，导致岩石圈地幔由古老难熔向新生饱满发生转变。在转变过程中，少量古老难熔的成分会残留到岩石圈地幔，最终在长白山地区形成以新生饱满为主，含少量古老难熔的岩石圈地幔。因此，太平洋板片俯冲不仅为新生代长白山地区提供了构造动力学背景，还为其岩石圈形成与破坏提供所需的物质。

5 结论

本文基于长白山地区地幔橄榄岩包体的岩相学、全岩主微量元素、单矿物原位主微量元素、单斜辉石Sr-Nd-Hf同位素特征研究，得到以下结论：

(1)长白山地区的橄榄岩包体主要是尖晶石相二辉橄榄岩和少量的方辉橄榄岩，平衡温度为900～1100℃，经历了的不同程度的部分熔融(二辉橄榄岩<10%；方辉橄榄岩>10%)，后期受到地幔交代作用影响，交代介质为富水硅酸盐熔体。

(2)长白山地区晚新生代岩石圈地幔表现出不均一的特征，以新生饱满地幔为主，夹有少量古老难熔的岩石圈地幔碎片。

(3)太平洋板块向西俯冲的过程使得长白山地区所处的华北克拉通东北部下覆软流圈上涌，携带大量壳源物质进入地幔深部，富水硅酸盐熔体上升过程中与岩石圈地幔发生橄榄岩-熔体反应，导致岩石圈地幔由古老难熔向新生饱满发生转变。

致谢感谢中国科学院地质与地球物理所雷鸣博士、成智慧及张茂亮博士后等在野外工作和数据分析中的大力支持和帮助。感谢匿名审稿人和期刊主编对论文提出的宝贵意见。