大同第四纪玄武岩成因：主微量元素及Sr-Nd-Pb-Hf同位素研究*

孙嘉祥 李霓 张雯倩

1. 吉林长白山火山国家野外科学观测研究站，中国地震局地质研究所，北京 1000292. 中国地震局地震与火山灾害重点实验室，北京 100029

来源于地幔的玄武岩是探究深部地质过程的窗口，其地球化学特征可以反演源区组成，并为该区的地幔演化提供线索。前人应用玄武岩地球化学特征深入探讨了华北克拉通的地幔特征和演化历史，但是大多数成果集中在克拉通东部玄武岩及其地幔包体的研究上(Nakamuraetal., 1989; Fan and Menzies, 1992; 郑建平, 1999; Gaoetal., 2002; Zhangetal., 2002, 2009; Yangetal., 2010; Zengetal., 2010; 陈立辉等, 2012; Guoetal., 2013)，对于克拉通西部地区除了汉诺坝地区因为玄武岩规模较大并且携带地幔包体研究程度较高以外(Songetal., 1990; Zhietal., 1990; Basuetal., 1991; 解广轰和王俊文, 1992; 刘丛强和解广轰, 1996; Liuetal., 2001; Xu, 2002; Zhouetal., 2002; Rudnicketal., 2004; Choietal., 2008; Zhengetal., 2009; Qianetal., 2015)，重力梯度带西侧其它玄武岩的研究资料相对较少(Zhangetal., 2005; 叶蕾等, 2015)，对于大同玄武岩的研究自二十世纪八十年代末才开始起步。

陈孝德等(1997)采集了大同玄武岩中的幔源包体，并对包体中橄榄石Ca含量进行测定，采用地质温压计证明了大同火山群的岩浆主要来源于60km左右深处的上地幔，且第四纪以来大同地区地温梯度有增高的趋势。解广轰和王俊文(1989)报道了少量大同玄武岩Sr-Nd同位素特征，根据Sr-Nd的线性关系以及玄武岩εNd(t)>0的特征推测大同玄武岩来源于亏损的上地幔且未受到地壳的混染。马金龙和徐义刚(2004)对大同玄武岩进行了较为系统的研究，发现大同玄武岩具有类似洋岛玄武岩(OIB)的地球化学特征，并提出DM+EMⅠ的混合模型来解释大同玄武岩的成因。Xuetal. (2005)进一步指出大同玄武岩主要受控于软流圈物质，并加入了少量岩石圈成分。上述学者都得出大同玄武岩来源于上地幔，但对于玄武岩部分熔融程度、源区残留相、地幔端元成分以及岩石圈作用的认识仍然不够清晰。本文拟通过研究大同玄武岩地球化学特征来揭示该区玄武岩的形成机制，为华北克拉通西部地区破坏提供新的证据。

1 区域地质背景

华北克拉通是我国东部重要的构造单元，也是地球上克拉通破坏的典型地区之一。在华北克拉通周围发育一系列俯冲/碰撞带，包括北部的古生代古亚洲洋俯冲、南部的三叠纪扬子地块俯冲、东面的侏罗-白垩纪太平洋板块俯冲，形成了北临中亚造山带，南接秦岭-大别造山带，东临太平洋板块的大地构造(图1)(Jahnetal., 1987; Liuetal., 1992; Maruyamaetal., 1997; Zhaoetal., 2001; Li and Santosh, 2017)。自中生代以来的岩石圈减薄是去克拉通化的重要标志之一，虽然众多学者对于减薄的时间和空间有争论，但是大部分学者还是倾向于认为早白垩世是减薄的高峰期(Xuetal., 2001, 2009; 裴福萍, 2004; 吴福元等, 2008;朱日祥等, 2011; Kuskyetal., 2014)。空间上，重力梯度带东西两侧岩石圈厚度存在差异(娄辛辉等, 2017; 徐小兵等, 2018)，西部岩石圈自新生代逐渐减薄，东部岩石圈逐渐加厚(徐义刚, 2006; Xuetal., 2009)。总体来说，无论在时间或空间上，东西部岩石圈的“热状态”、厚度以及地球化学性质都有着明显的不同(Menziesetal., 1993; Niu, 2005; Xu, 2007; 吴福元等, 2008; Xuetal., 2009; Zhuetal., 2012; Guoetal., 2014)。

图1 华北克拉通构造简图(a,据Yang et al., 2010修改)及大同第四纪玄武岩分布简图(b,据马金龙和徐义刚, 2004)DTGL-大兴安岭-太行山重力梯度带；TLFZ-郯庐断裂带；①汉诺坝；②集宁；③大同；④繁峙Fig.1 Geological map of the North China Craton (a, modified after Yang et al., 2010) and the Quaternary basalt in Datong (b, after Ma and Xu, 2004)DTGL-Daxing’anling-Taihangshan Gravity Lineament; TLFZ-Tanlu Fault Zone; ①Hannuoba; ②Jining; ③Datong; ④Fanshi

大同火山区分布在华北克拉通太行山-大兴安岭重力梯度带西侧，位于山陕地堑最北端盆地，处于中国东部NNE向剪切拉张断陷盆地的构造地貌格局之中，北界是阳高-天镇断裂，南至六棱山山前断裂，西临口泉断裂(王乃樑等, 1996; 岑敏等, 2015)。大同第四纪玄武岩出露范围大致位于39°55′22″～40°10′35″N，113°34′17″～113°56′45″E。前人以陈庄-许堡为界将大同火山区大致分为东西两个部分，西区玄武岩多以爆发和溢流相产出，形成大小形状不同的30余座火山锥，由于西区玄武岩来源深度较深，所以形成低硅高碱的碱性玄武岩；东区火山爆发系数较小，常沿着断裂位置(桑干河、六棱山断裂)溢出，形成溢流玄武岩，由于东区玄武岩来源较浅，形成高硅低碱的拉斑玄武岩(韩军青, 1992; 李文宣等, 1994; 陈孝德等, 2001; 安卫平和苏宗正, 2008)。前人通过释光法、古地磁以及K-Ar测年等方法对大同玄武岩喷发年龄进行了研究(裴静娴, 1981; 李虎侯和孙建中, 1984; 朱日祥等, 1985, 1990; 陈文寄等, 1992; 赵华等, 2012)，虽然测定结果不完全一致，但是总体上大同玄武岩喷发可分为三个阶段：第一阶段是以东区册田水库溢流玄武岩为代表距今0.74～0.54Ma的喷发；第二阶段西区爆破相和东区溢流相玄武岩均在0.4～0.3Ma左右大规模产出；第三阶段的大规模喷发以秋林村山前断裂溢流玄武岩为代表，年龄在0.2Ma左右，即大同玄武岩的喷发时代发生在早更新世晚期-晚更新世早期(李虎侯和孙建中, 1984; 陈文寄等, 1992)。本文所采集的样品包括西区碱性玄武岩和东区拉斑玄武岩。

2 样品的采集与分析

大同玄武岩呈灰-灰黑色，块状构造，少数玄武岩含有排列规则的气孔，气孔中充填白色碳酸盐矿物呈杏仁状构造。本文采样点涵盖东区和西区，分析所用样品为新鲜致密块状灰黑色玄武岩。与汉诺坝玄武岩含有丰富的地幔橄榄岩包体相比，仅在大同金山、狼窝山处发现少量尖晶石二辉橄榄岩包体。大同玄武岩呈块状构造，显微斑状结构。西区玄武岩斑晶含量约为20%～30%，基质含量约为40%～50%。斑晶为橄榄石(ol)和单斜辉石(cpx)，橄榄石常呈半自形或他形粒状，单斜辉石多为半自形粒状。基质由橄榄石、单斜辉石、斜长石(pl)、铁钛氧化物以及玻璃质组成。东区玄武岩主要矿物成分与碱性玄武岩基本一致，差别在于斑晶组成上，橄榄石含量下降，辉石含量上升；基质中斜长石微晶明显增多，其组成的三角格架中充填他形粒状辉石、橄榄石和磁铁矿的细小颗粒(图2)。

图2 大同东区溢流玄武岩野外照片(a)和镜下显微照片(b)Fig.2 Field (a) and microscopic (b) photographs of the eastern Datong basalts

样品的主量、微量和同位素分析测试均在武汉上谱分析科技有限责任公司完成。挑选新鲜样品，切去风化面，研磨至5mm， 剔除杂质， 最后磨至200目粉末。全岩主量元素在波长色散X射线荧光光谱仪(ZSXPrimusⅡ)上完成，标样采用GBW07105(标准值)来保证测试精度，分析精度为1%～5%。微量元素测试在电感耦合等离子体质谱仪(Agilent 7700e ICP-MS)上采用GB/T14506.30—2010硅酸盐岩石化学分析方法完成，分析精度好于10%。具体流程：(1)将200目样品置于105℃烘箱中烘干12小时；(2)准确称取粉末样品50mg置于Teflon溶样弹中；(3)依次缓慢加入1mL高纯硝酸和1mL高纯氢氟酸；(4)将Teflon溶样弹放入钢套，拧紧后置于190℃烘箱中加热24小时以上；(5)待溶样弹冷却，开盖后置于140℃电热板上蒸干，然后加入1mL硝酸并再次蒸干(确保溶样弹壁无液体)；(6)加入1mL高纯硝酸、1mL MQ水和1mL内标In(浓度为1×10-6)，再次将Teflon溶样弹放入钢套，拧紧后置于190℃烘箱中加热12小时以上；(7)将溶液转入聚乙烯塑料瓶中，并用2%硝酸稀释至100g以备ICP-MS测试。岩石主微量分析结果见表1。Sr-Nd-Pb-Hf同位素测试采用BCR-2、RGM-2两个标样在MC-ICP-MS上完成，分析结果见表2。同位素分析利用美国Thermo Fisher Scientific公司的MC-ICP-MS(Neptune Plus)完成。

表1 大同玄武岩主量元素(wt%)和微量元素(×10-6)含量

续表1

图3 大同玄武岩TAS图解(底图据Le Bas et al., 1986；碱性玄武岩和拉斑玄武岩分界虚线据Irvine and Baragar, 1971)实心图标数据引自马金龙和徐义刚(2004)；图9同Fig.3 TAS diagram of Datong basalts (after Le Bas et al., 1986, the boundary line between alkaline basalt and tholeiitic basalt after Irvine and Baragar, 1971)The data of solid triangles and squares from Ma and Xu, 2004

3 分析结果

3.1 主量元素

表1经标准化计算得知，大同西区玄武岩SiO2含量为45.02%～49.07%，K2O+Na2O含量为5.23%～6.53%，TiO2含量为2.46%～3.05%；东区玄武岩SiO2含量为49.76%～53.30%，K2O+Na2O含量为3.60%～4.56%，TiO2含量为1.81%～2.28%。相较于东区玄武岩，西区玄武岩整体呈低硅高碱高钛的特征。在TAS投图上(图3)，大同西区玄武岩属于碱性玄武岩系列，东区玄武岩属于拉斑玄武岩系列。主量元素Fe2O3、TiO2、Na2O、K2O与MgO之间均无明显协变关系，这和汉诺坝、繁峙、集宁玄武岩类似，但可以看出相对于东区拉斑玄武岩，西区碱性玄武岩具有明显高的主量元素含量(图4)。SiO2-MgO图解表明，随着SiO2含量增加，MgO含量降低(图4a)；CaO/Al2O3-MgO图解显示简单的正相关(图4b)；大同碱性玄武岩和拉斑玄武岩MgO含量虽然不同，但两种玄武岩Mg#并没有明显的演化关系。

图4 大同第四纪玄武岩主量元素氧化物与MgO协变图汉诺坝，繁峙，集宁数据来自Zhi et al., 1990; Zhang et al., 2005; 叶蕾等, 2015Fig.4 Diagrams of major element oxides vs. MgO for Quaternary Datong basaltsData sources of Hannuoba, Fanshi, Jining after Zhi et al., 1990; Zhang et al., 2005; Ye et al., 2015

3.2 微量元素

大同玄武岩微量元素分析结果见表1，其中碱性玄武岩的稀土元素总量为192.1×10-6～254.6×10-6，明显高于拉斑玄武岩的69.25×10-6～103.7×10-6含量范围。在稀土元素配分图上，无论西区碱性玄武岩还是东区拉斑玄武岩均显示富集轻稀土(LREE)、亏损重稀土(HREE)的配分模式(图5a)，但两种玄武岩的轻重稀土分馏程度不同，碱性玄武岩(La/Yb)N=19.74～31.62，(La/Sm)PM=2.90～3.40；拉斑玄武岩(La/Yb)N=5.82～12.27，(La/Sm)PM=1.39～2.37，可见碱性玄武岩的轻重稀土分馏程度明显更高。此外，两种玄武岩均未显示Eu异常。原始地幔标准化蛛网图中显示两种玄武岩类似洋岛玄武岩(OIB)的特征，富集Ba、Sr、K等大离子亲石元素(LILE)，Nb、Ta两种高场强元素(HFSE)轻微富集，在图上基本显示平滑曲线，且不亏损Zr的特征(图5b)。

3.3 Sr-Nd-Pb-Hf同位素

大同玄武岩Sr-Nd-Hf同位素组成分别是87Sr/86Sr为0.703302～0.705102，143Nd/144Nd为0.512561～0.512963，以及176Hf/177Hf为0.282922～0.283072(表2)。Sr-Nd二元同位素图解显示大同玄武岩整体具有负相关特征，并落在洋岛玄武岩(OIB)范围内，接近全硅酸盐地球值(BSE)(图6a)。εNd(t)为-1.5～6.3，εHf(t)为5.3～10.6，在εNd(t)-εHf(t)图解(图6b)上，大同玄武岩样品沿着地幔演化趋势线呈正相关分布，并且也落在OIB范围内。相较于大同拉斑玄武岩，大同碱性玄武岩具有较高Nd和较低Sr同位素特征，显示了同位素组成与岩性的相关关系。大同玄武岩Pb同位素比值如下：206Pb/204Pb为17.131～18.266，207Pb/204Pb为15.382～15.531，208Pb/204Pb为37.591～38.381，在208Pb/204Pb-206Pb/204Pb图解(图7a)显示大同碱性玄武岩落在I-MORB(印度洋MORB)范围内，而拉斑玄武岩则更接近五大连池钾质玄武岩，这点在207Pb/204Pb-206Pb/204Pb图解(图7b)得到了验证。所有研究区样品都在NHRL (北半球参考线)之上，且相较于富集端元EMⅡ具有更低的放射性Pb同位素组成，同时可以看出大同碱性玄武岩Pb同位素比值要高于拉斑玄武岩，这些特征与重力梯度带附近繁峙、汉诺坝、集宁等地玄武岩类似(解广轰等, 1989; Xu, 2002; 马金龙和徐义刚, 2004; Zhangetal., 2005; Chenetal., 2007; Yan and Zhao, 2008; Yangetal., 2010; 叶蕾等, 2015)。

表2 大同玄武岩Sr-Nd-Pb-Hf同位素组成

图5 大同玄武岩球粒陨石标准化稀土元素配分曲线(a)和原始地幔标准化不相容元素曲线(b)球粒陨石、原始地幔、OIB、E-MORB和N-MORB参考值据Sun and McDonough, 1989Fig.5 Chondrite-normalized REE patterns (a) and primitive-mantle normalized incompatible element curves (b) for Datong basaltsChondrite, primitive mantle, OIB, E-MORB and N-MORB reference values after Sun and McDonough, 1989

图6 大同及周边玄武岩87Sr/86Sr-143Nd/144Nd图解(a)和εNd(t)-εHf(t)图解(b)(a)数据来源：PREMA、MORB、OIB、EMⅠ和EMⅡ (Zindle and Hart, 1986)、古老岩石圈地幔(Zhang et al., 2002)、集宁玄武岩(Zhang et al., 2005)、汉诺坝玄武岩(Song et al., 1990)、繁峙玄武岩(叶蕾等，2015)；(b)数据来源：汉诺坝(Choi et al., 2008)、集宁玄武岩(Zhang et al., 2005; Ho et al., 2011; Zhao et al., 2013)、方城玄武岩(Zhang et al., 2002)，繁峙玄武岩数据同图6aFig.6 Diagrams of 87Sr/86Sr vs. 143Nd/144Nd (a) and εNd(t) vs. εHf(t) (b) for basalts at Datong and adjacent areasData sources in Fig.6a: MORB, OIB, EMⅠ, PREMA and EMⅡ (Zindle and Hart, 1986), ancient lithosphere mantle (Zhang et al., 2002), Jining basalt (Zhang et al., 2005), Hannuoba basalt (Song et al., 1990), Fanshi basalt (Ye et al., 2015). Fig.6b data source: Hannuoba (Choi et al., 2008), Jining basalt (Zhang et al., 2005; Ho et al., 2011; Zhao et al., 2013)，Fangcheng basalt (Zhang et al., 2002), Fanshi basalt data are the same as those of Fig.6a

图7 大同及其它玄武岩208Pb/204Pb-206Pb/204Pb(a)和207Pb/204Pb-206Pb/204Pb(b)图解数据来源：太平洋、大西洋MORB和印度洋MORB (Barry and Kent, 1998; Zou et al., 2000; Chauvel and Blichert-Toft, 2001)、汉诺坝玄武岩(Song et al., 1990; Basu et al., 1991)、五大连池玄武岩(Zhang et al., 1998; Zou et al., 2003)、NHRL (Hart, 1984)，其他图例同图6Fig.7 Plots of 208Pb/204Pb vs. 206Pb/204Pb (a) and 207Pb/204Pb vs. 206Pb/204Pb (b) for Datong and other basaltsData sources: Pacific, Atlantic and Indian MORB (Barry and Kent, 1998; Zou et al., 2000; Chauvel and Blichert-Toft, 2001), Hannuoba basalt (Song et al., 1990; Basu et al., 1991), Wudalianchi basalt (Zhang et al., 1998; Zou et al., 2003), NHRL (Hart, 1984), other legends same as in Fig.6

图8 大同玄武岩La/Nb-Ba/Nb图解(a)和Nb/U-Nb图解(b)数据源：OIB、N-MORB (Sun and McDonough, 1989)，全大陆地壳BCC (Rudnick and Gao, 2003)Fig.8 Plot of La/Nb vs. Ba/Nb (a) and Nb/U vs. Nb (b) for Datong basaltsData sources: OIB, N-MORB (Sun and McDonough, 1989), BCC (Rudnick and Gao, 2003)

图9 大同玄武岩MgO-Ni(a)和MgO-Cr(b)图解Fig.9 Diagrams of MgO vs. Ni (a) and MgO vs. Cr (b) for Datong basalts

图10 大同玄武岩La/Yb-Sm/Yb (a)和Zr-Zr/Y (b)图解曲线旁边数字代表熔融比例，源区矿物组分、矿物熔融比例和元素分配系数据Johnson et al., 1990; McKenzie and O’Nions, 1991Fig.10 La/Yb vs. Sm/Yb (a) and Zr vs. Zr/Y (b) diagrams of Datong basaltThe numbers next to the curve represent the melting ratio. Mineral composition, melting ratio, and elemental distribution data of the source area after Johnson et al., 1990; McKenzie and O’Nions, 1991

4 讨论

4.1 地壳混染

克拉通背景下的大陆玄武质岩浆在上升至地表过程中会经历较厚的大陆地壳，存在着被地壳混染的可能性(Carteretal., 1978; Watson, 1982; Huppertetal., 1985; Reinersetal., 1995)。大同玄武岩是华北克拉通西部玄武岩的重要组成部分，部分拉斑玄武岩在微量元素蛛网图上显示Pb正异常，暗示岩浆同样可能被大陆地壳混染，但是基于以下事实可以排除地壳混染的可能性：(1)狼窝山和金山所见的玄武岩含有橄榄岩包体，这一般被认为是玄武质岩浆快速上升的标志，所以与地壳混染的可能性不大。(2)大陆地壳相较于地幔来源的物质具有低Nb、富集LREE、更低的Nb/U和更高的La/Nb比值(Hofmann, 1988)，如果大同玄武岩经历了地壳混染，则Nb/U-Nb以及La/Nb-Ba/Nb应该有正相关关系，而在图8a和图8b中并未发现大同碱性玄武岩样品具有这种特征，且投点落在OIB范围内，具有比地壳更低的La/Nb比值，更高的Nb含量和Nb/U比值，远离地壳投点。(3)在微量元素蛛网图中并未观察到Nb、Ti、Ta亏损的特征，相反绝大部分样品(Nb/La)N>1，并且显示Ti、Ta的轻微富集，所以基本可以排除地壳混染的可能性。大同拉斑玄武岩La/Nb-Ba/Nb正相关，Nb、Ti、Ta元素相较于碱性玄武岩含量更低，微量元素蛛网图显示Pb正异常，可能暗示地壳混染，但是Ce/Pb-Nb/U图解显示拉斑玄武岩仍然远离地壳投点(图略)，且如果经历了较大程度的地壳混染，则εNd值应该与MgO含量存在相关性，但是拉斑玄武岩并没有显示这样的特征，这与AFC预测的混染模型不一致，总体而言微量元素仍然以类似OIB特征为主，所以拉斑玄武岩经历了轻微的地壳混染，但根据上述讨论可以排除地壳混染对于玄武岩成因的影响。基于这个前提，大同玄武岩样品的地球化学特征可以更好地反映地幔源区的特征，这与前人对重力梯度带附近的新生代玄武岩区地壳混染程度不高的评估结果一致(Zhietal., 1990; Xuetal., 2005; Zhangetal., 2005; 刘金菊等, 2014; 叶蕾等, 2015)。

4.2 分离结晶作用

大同玄武岩MgO变化范围为6.63%～8.04%，低于中国东部新生代原始岩浆参考值(MgO 10%～12%)，样品的Ni(130×10-6～175×10-6)、Cr(165×10-6～265×10-6)含量又远远低于地幔原始岩浆的含量(Ni>235×10-6, Cr>400×10-6, Sato, 1977; Wilson, 1989)，镜下观察到大同玄武岩含有橄榄石、单斜辉石斑晶和斜长石微晶，这些都是大同玄武岩发生分离结晶作用的重要证据。由于Ni在橄榄石中的分配系数为5.9～29，Cr在单斜辉石中的分配系数约为34(Arth, 1976)，图9中显示随着MgO的降低，Ni、Cr含量降低，所以暗示岩浆过程中发生了橄榄石和单斜辉石的分离结晶作用。同时观察到碱性玄武岩MgO与Cr、Ni的正相关趋势明显，而拉斑玄武岩的正相关趋势则稍弱，表明碱性玄武岩分离结晶程度更高；另外CaO/Al2O3-MgO协变图(图4b)的良好正相关关系也是大同碱性玄武岩和拉斑玄武岩经历单斜辉石分离结晶的证据。微量元素蛛网图(图5b)并未见Sr、Eu的负异常，说明斜长石的分离结晶作用不明显。总的来说，虽然两种玄武岩均经历了橄榄石、单斜辉石的分离结晶作用，未见明显斜长石结晶，但大同碱性玄武岩分离结晶程度要高于拉斑玄武岩，这与重力梯度带附近汉诺坝等地玄武岩的结晶演化作用基本一致(Song and Frey, 1989; Zhietal., 1990; Liuetal., 2001; Xuetal., 2005; 叶蕾等, 2015)。

4.3 岩浆源区

中国东部新生代玄武岩被认为是源区两端元不同混合比例形成的(Pengetal., 1986; Basuetal., 1991; Fan and Hooper, 1991; Tatsumotoetal., 1992; Liuetal., 1994)。大同附近的汉诺坝、集宁和繁峙等重力梯度带周边的新生代玄武岩源区同样被认为是来自二端元的混合(解广轰等, 1989; 刘丛强和解广轰, 1996; Coulonetal., 1996; Tangetal., 2006; Zhangetal., 2012; Guoetal., 2014)。那么大同玄武岩的岩浆源区是否同样如此呢，这可从所采集样品的主微量和同位素的地球化学特征来推断其源区成因。

微量元素蛛网图和稀土配分图(图5b)均显示大同玄武岩有类似OIB的元素特征，和MORB相比非常富集Rb、Ba、Th、U、La等强不相容元素，碱性玄武岩的上述元素含量要高于拉斑玄武岩，这些强不相容元素的含量一般和源区性质有关，暗示两种岩性的岩浆源区有差异。这种类似OIB富集不相容元素的特征和重力梯度带西侧集宁、汉诺坝、繁峙等地玄武岩极为类似(解广轰等, 1989; Zhietal., 1990; 马金龙和徐义刚, 2004; Xuetal., 2005; Zhangetal., 2005; 叶蕾等, 2015)，后者被认为来自软流圈地幔(Song and Frey, 1989; Basuetal., 1991; Fan and Hooper, 1991)，结合大同玄武岩右倾的稀土元素配分模式要求源区残留石榴石，据此推测大同玄武岩也可能来自软流圈地幔。

143Nd/144Nd-87Sr/86Sr、207Pb/204Pb-206Pb/204Pb图解(图6a、图7b)可以看出大同玄武岩主要来自PREMA与EMⅠ端元的混合，并有少量EMⅡ端元的影响，故PREMA是大同玄武岩的亏损端元。PREMA是由国外学者根据全球不同地方玄武岩的同位素特征划分出来的地幔端元中的一个，它的主要特征有(1)与洋岛玄武岩极为相似的同位素特征，富集地壳亏损的元素；(2)可以代表很多大陆玄武岩的亏损端元；(3)具有PREMA同位素特征的玄武岩来源深度可能大于100km(Zindleetal., 1982; Zindle and Hart, 1986)。大同玄武岩在143Nd/144Nd-87Sr/86Sr、εNd(t)-εHf(t)图解(图6)均落在OIB范围内，根据现有地球物理资料显示大同部分区域岩石圈厚度要超过100km(娄辛辉等, 2017; 徐小兵等, 2018)，并且大同附近的汉诺坝和集宁玄武岩Pb同位素特征也显示了PREMA作为亏损端元参与了新生代玄武岩的形成(解广轰和王俊文, 1992; Zhangetal., 2005)，所以认为来自软流圈的PREMA是大同玄武岩的亏损端元。

Sr-Nd二元图解显示大同拉斑玄武岩可能受到少量富集组分的影响(图6a)，但是不能确定富集端元EMⅠ和EMⅡ哪一个参与得更多。首先，EMⅡ组分相对大同玄武岩具有较高的206Pb/204Pb和87Sr/86Sr比值，而且它一般起源于大陆下地壳或陆源沉积物(Carlson, 1984; Zindle and Hart, 1986)，前文中的讨论以及微量元素蛛网图上也显示大同玄武岩并未有明显的地壳混染特征，因此可以排除EMⅡ作为主要端元参与了大同玄武岩的形成，但不排除有少部分EMⅡ的参与。图7显示大同玄武岩均有良好的从PREMA到EMⅠ的混合趋势，所以推断EMⅠ是主要的富集端元。关于EMⅠ端元的形成，国内外学者看法不一致，部分学者推测EMⅠ可能来自古洋壳及其携带的沉积物俯冲至深部交代富集软流圈地幔形成(McDonough, 1990; Weaver, 1991; 石林等, 1998)，另一部分学者认为EMⅠ来源于古老克拉通之下的大陆岩石圈地幔下部(Tatsumotoetal., 1992)。大同玄武岩的富集端元EMⅠ可能是第二种成因，即来自于古老岩石圈地幔。证据如下：(1)同为重力梯度带西侧且临近大同的汉诺坝、集宁等地的大陆岩石圈地幔下部具有和EMⅠ一样低的143Nd/144Nd和高的87Sr/86Sr同位素特征(Tatsumotoetal., 1992; Tangetal., 2006)；(2)中国东部包括上述区域新生代玄武岩EMⅠ型特征也被广泛认为是大陆岩石圈地幔(SCLM)下部的贡献(Zhietal., 1990; Basuetal., 1991; Fan and Hooper, 1991; Tatsumotoetal., 1992; 马金龙和徐义刚, 2004; Zhangetal., 2005; Hoetal., 2011; 叶蕾等, 2015)；(3)208Pb/204Pb-206Pb/204Pb图解上(图7a)，大同拉斑玄武岩更靠近五大连池玄武岩(板块内部SCLM来源)，因此本文认为来自古老岩石圈地幔的EMⅠ是大同玄武岩的富集端元。

同位素特征揭示了大同玄武岩是软流圈(PREMA)+大陆古老岩石圈地幔下部(EMⅠ)形成的。因为古老岩石圈地幔一般比较难熔，但是玄武岩的形成需要饱满岩石熔融(Griffinetal., 1999)，Sr-Nd组分混合模拟曲线显示大同拉斑玄武岩有不超过20%、大同碱性玄武岩不超过10%的富集组分的参与(图6a)，所以大同第四纪玄武岩主要来源于软流圈的部分熔融，并熔融了少量古老岩石圈地幔。同位素图解上碱性玄武岩更靠近亏损端元(PREMA)，拉斑玄武岩更靠近富集端元(EMⅠ)，说明拉斑玄武岩中岩石圈组分比例较碱性玄武岩大。

4.4 部分熔融程度和深度

玄武岩部分熔融的深度和硅的饱和程度有较大关系，通常硅越饱和玄武岩形成深度越浅，所以碱性玄武岩形成深度要大于拉斑玄武岩(Hirose and Kushiro, 1993; Kushiro, 2001)。前人研究表明华北碱性玄武岩形成深度>80km，拉斑玄武岩的形成深度为50～60km(Nohdaetal., 1991)。大同西区主要为碱性玄武岩，形成深度应该大于80km。主量上碱性玄武岩SiO2含量大于拉斑玄武岩，在硅碱投图上处于硅不饱和区域，拉斑玄武岩则相对硅饱和。微量上(La/Sm)碱性>(La/Sm)拉斑、(La/Yb)碱性>(La/Yb)拉斑，显示拉斑玄武岩轻重稀土分馏程度低于碱性玄武岩，这种轻重稀土的分馏程度受控于岩石圈厚度(源区深度)，岩石圈厚度越大，发生部分熔融程度越低，轻重稀土的分馏程度就越明显，即源区来源越深的玄武岩，其轻重稀土越分馏(Ellam, 1992; Niu, 2005, 2016; Guoetal., 2020)，所以东区拉斑玄武岩形成深度应浅于西区碱性玄武岩。

对于石榴子石，Yb是相容元素，而La、Sm为不相容元素，石榴子石相橄榄岩部分熔融的程度越低，分异程度越明显；而在尖晶石相橄榄岩部分熔融作用中，La/Yb变化较小，Sm/Yb基本不变，因此La/Yb-Sm/Yb图常用于区分来自石榴子石相橄榄岩和尖晶石相橄榄岩的玄武岩(Xuetal., 2005)。从图10中可以看出大同玄武岩落在石榴子石二辉橄榄岩熔融曲线上，碱性玄武岩熔融程度较小为1.5%～3%，拉斑玄武岩熔融程度明显比碱性玄武岩大，约为4%～8%。不相容元素的含量和比值通常也是判断部分熔融程度的重要方法，Zr/Y比值不受分离结晶影响而受控于部分熔融程度，这是因为Zr的不相容性更高，随着熔融程度的增加，Zr/Y比值下降，图10b显示拉斑玄武岩更小的Zr/Y比值，这与图10a的模拟结果一致。

5 结论

(1)大同第四纪火山岩以陈庄-许堡断裂为界分为东西两区，西区是碱性玄武岩，东区为拉斑玄武岩。

(2)大同玄武岩含有大量橄榄石和单斜辉石斑晶，结合Ni、Cr与MgO的负相关特征，推断这两种矿物是主要的分离结晶矿物。

(3)La/Yb-Sm/Yb图解模拟计算得出大同玄武岩是石榴子石二辉橄榄岩低程度部分熔融的结果，碱性玄武岩约为1.5%～3%，拉斑玄武岩为4%～8%，部分熔融程度不同控制着两种玄武岩主、微量元素的差异。

(4)大同碱性玄武岩和拉斑玄武岩相对富集LREE、LILE和HFSE，具有类似洋岛玄武岩(OIB)的特征，同位素图解表明它们都是以PREMA端元为代表的软流圈地幔低程度部分熔融，并有较少比例以EMⅠ端元为代表的古老岩石圈地幔的加入(碱性玄武岩不超过10%，拉斑玄武岩不超过20%)，其同位素特征的差异则是岩石圈物质参与比例不同造成的。

致谢感谢三位评审专家及编辑部俞良军博士对本文的建设性修改意见，使本文的论述及表达更为严谨清晰。