冀东地区早侏罗世岩浆岩成因及其意义*

陈竟志 姜能 范文博 胡俊

1. 中国科学院地质与地球物理研究所，岩石圈演化国家重点实验室，北京 1000292. 中国科学院大学地球与行星科学学院，北京 1000493. 西北大学地质学系，大陆动力学国家重点实验室，西安 7100691.

作为地球上稳定的古老陆块之一，华北克拉通东部在中生代时期经历了重要的岩石圈减薄与克拉通破坏过程(朱日祥等, 2012)。中生代火成岩，特别是花岗岩在华北克拉通东部地区广泛分布，它们是华北克拉通破坏过程中的重要地质响应，记录了地球深部物质交换与区域地质演化的许多证据。已有的大量年代学研究资料显示，华北克拉通中生代岩浆岩主要形成于晚三叠世、侏罗纪和早白垩世三个时代。其中，晚三叠世、早白垩世岩浆岩受到广泛关注，它们被认为分别记录了华北克拉通破坏起始与峰期的时间(Zhangetal., 2014)。相比而言，侏罗纪的岩石在华北克拉通相对较少，已经识别出来的岩石主要集中在中-晚侏罗世时期(杨进辉等, 2007)。其中，早侏罗世岩石在华北普遍较为稀少，因此对这一时期的深部过程及其与克拉通破坏的联系，研究也十分有限。

冀东地区位于华北克拉通东部陆块，沿马兰峪背斜核部一带出露了大量的中生代侵入体，但对其年代学研究及不同时代岩体之间的成因联系，研究程度仍相对有限，对一些岩体的形成时代认识也很不足(马君, 2009; 杨付领等, 2015)。在对区域中生代岩浆岩开展系统研究的过程中，本文发现早侏罗世岩体在这一区域十分发育。这为我们探讨早侏罗世岩浆岩的成因以及与之相关的深部过程，提供了十分良好的地质素材。为此，本文以冀东青龙地区肖营子岩基内识别出的早侏罗世岩石为例，结合区域前人已有研究资料，对华北早侏罗世侵入岩进行系统研究，以揭示其成因，并探索其与克拉通破坏的关系。

1 区域地质背景

华北克拉通是世界上最古老的克拉通之一，由早太古代至早元古代的结晶基底和覆盖其上的沉积盖层组成。中生代时期，华北克拉通东部地区经历了广泛的构造-热事件，与之伴随的是岩石圈的减薄与克拉通的破坏，克拉通的破坏在早白垩世达到顶峰(吴福元等, 2008; 朱日祥等, 2012; Zhangetal., 2014)。

冀东青龙地区位于华北克拉通东部，区内中生代岩浆活动强烈(罗镇宽等, 2001; 陆继龙等, 2012; 叶浩等, 2014; 徐希阳等, 2016; Fanetal., 2017; 熊乐等, 2017)，形成了诸如都山、柏杖子、肖营子、青山口、罗文峪、贾家山、牛心山、峪耳崖、碾子峪等许许多多的侵入岩(图1)。肖营子岩基是区域上较大的一个侵入体，它位于冀东青龙县南部肖营子镇一带，呈不规则椭圆状，北东向展布，出露面积大约有230km2。肖营子岩基的围岩为太古宙变质岩系和中元古代高于庄组地层(图1)。前人曾将肖营子岩基视为同一时代的产物，并根据不同的测年结果解释将其置于侏罗纪的不同时代(马君, 2009; 杨付岭等, 2015; 尹业长等, 2018)。我们近年来的详细研究发现，它实际上至少由中侏罗世的白家店、早侏罗世的蛇盘兔两部分构成。本文将对研究区早侏罗世的蛇盘兔单元进行详细报道。蛇盘兔岩体位于岩基北东部位，出露面积约170km2，主体岩石为碱长花岗岩，其次为石英斑岩、正长花岗岩等(图2)。在其内部，还伴生出现了规模较小(～3km2)的龙须门闪长岩岩体(图1、图2)。蛇盘兔碱长花岗岩主要由石英(40%)、钾长石(35%)、斜长石(20%)组成，副矿物有磁铁矿、榍石等(图3)。花岗斑岩的矿物组合为黑云母(～5%)+石英(30%)+斜长石(20%)+钾长石(40%)，副矿物有磁铁矿、锆石、榍石等(图3)。龙须门闪长岩显示出半自形粒状结构，矿物组合为黑云母(～10%)、石英(～10%)、角闪石(20%～30%)、斜长石(50%～60%)，副矿物有磁铁矿、榍石、锆石和磷灰石，电子探针分析显示其斜长石类型为中长石(图3)。本文对早侏罗世蛇盘兔花岗岩与龙须门闪长岩进行了采样，采样位置详见图1；并开展了详细的年代学、岩石学和地球化学研究。

图1 华北克拉通及研究区示意图(a)、冀东地区侵入体分布图(b)及冀东青龙地区区域地质简图(c) (据Fan et al., 2017修改)

图2 冀东地区蛇盘兔与龙须门岩体野外照片

图3 冀东地区蛇盘兔与龙须门岩石显微照片

2 分析测试方法

2.1 锆石U-Pb年龄

锆石单矿物在河北廊坊物化探实验室分选完成，锆石制靶及阴极发光图像拍摄在中国科学院地质与地球物理研究所进行。锆石U-Pb年龄分析在中国科学院地质与地球物理研究所离子探针实验室完成，所用仪器为CAMECA IMS-1280二次离子质谱仪。

U-Pb年龄详细分析方法见Lietal. (2009)。U-Th-Pb同位素比值用标准锆石Plésovice校正获得，U含量采用标准锆石91500校正获得，以长期监测标准样品获得的标准偏差和单点测试内部精度共同传递得到样品单点误差，以标准样品Qinghu作为未知样监测数据的精确度。普通Pb校正采用实测204Pb值。单点分析的同位素比值及年龄误差均为1σ；加权平均值误差为2σ，U-Pb的加权平均年龄、谐和年龄误差在95%的置信区间内。数据结果处理采用Isoplot(ver.3.0)软件完成(Ludwig, 2003)。

2.2 锆石Hf同位素分析

锆石原位Hf同位素测试在中国科学院地质与地球物理研究所多接收等离子体质谱仪实验室完成。分析时，所用仪器为装配有193nm激光剥蚀系统的Neptune多接收电感耦合等离子体质谱仪。采用He作为剥蚀物质载气，激光束斑直径为63μm，激光脉冲频率为8～10Hz，激光束的脉冲能量为100J。每个分析点的激光剥蚀时间26s，气体背景采集时间为30s，信号采集时间为30s。仪器的运行条件，详细的分析过程和数据处理见Wuetal. (2006)。分析过程中，用国际标样91500锆石为外标。数据结果计算中，衰减常量λLu取1.867×10-11y-1(Söderlundetal., 2004)，普通球粒陨石的176Hf/177Hf与176Lu/177Hf的比值以及现今的亏损地幔估计值分别为0.282772、0.0332和0.0384(Blichert-Toft and Albarède, 1997)。

2.3 全岩主微量元素分析

主量元素分析在西北大学大陆动力学国家重点实验室完成。首先称取0.5g样品粉末放入陶瓷坩埚中，在1100℃高温下灼烧，恒重后测定烧失量。然后将样品转移至玛瑙研钵中，加入5g Li2B4O7试剂，混合均匀后完全转移至铂金坩埚，滴入4滴发泡剂(NH4Br)，利用M4自动熔样机高温熔融，制成玻璃片上机测试。测试工作通过X-射线荧光光谱仪(AXIOS和RIX2100)完成。分析方法采用标准曲线法，岩石标样为GSR-1和GSR-3。样品的分析精度为～1%(含量>10%)和～10%(含量<1.0%)。微量元素分析在中国地质大学(武汉)地质过程与矿产资源国家重点实验室完成，所用仪器为Agilent 7500a ICP-MS，样品制备采用混合酸溶样法完成。数据分析的精确度优于10%，其中稀土元素的分析进度优于5%。

3 分析测试结果

3.1 锆石微区分析

图4 冀东地区蛇盘兔与龙须门岩体锆石CL照片

表1冀东地区早侏罗世岩浆岩SIMS锆石U-Pb年龄

Table 1 SIMS zircon U-Pb ages of the Early Jurassic igneous rocks from eastern Hebei Province

测点号UTh(×10-6)Th/U同位素比值年龄(Ma)207Pb/235U±σ206Pb/238U±σ207Pb/235U±σ206Pb/238U±σJN09114蛇盘兔花岗岩17659481.2400.214261.860.03081.501973.341952.8933062280.7450.201692.670.03021.531874.561922.91484311141.3220.212701.830.03111.501963.261972.9261131060.9380.208163.240.02981.531925.691892.8582482270.9130.206822.730.03091.521914.761962.9392712360.8730.197152.600.02981.511834.351902.83105266091.1580.205982.580.02981.551904.481892.89112572310.8960.184027.480.03001.6817211.871913.15128668230.9500.199012.240.02881.511843.791832.72146226601.0620.196584.840.02931.531828.111862.80166226721.0810.214222.040.03121.501973.671982.93184694220.9000.206452.070.03081.521913.601962.931994610691.1310.204722.480.03041.501894.281932.86205455240.9610.207182.070.03051.501913.611932.86211952031.0390.211893.010.03051.621955.361933.10JN09116龙须门闪长岩167400.5870.216755.370.02991.551999.761902.9121822931.6120.198633.120.02911.521845.271852.7631451861.2820.217102.950.03031.571995.361932.98582520.6390.208784.170.03001.661937.351903.12652290.5460.206364.540.02991.901907.911903.56769400.5790.213663.970.02941.591977.121872.9381021161.1340.205925.050.03001.501908.791902.82977490.6330.210213.760.03051.501946.651942.871075450.5980.206635.360.02991.551919.361902.901168590.8680.202364.670.02981.501878.011892.80121451761.2110.195642.990.03041.551814.981932.94131241461.1750.218093.100.03011.532005.661912.881484520.6240.213244.250.02981.541967.611902.881555310.5540.187554.460.03031.561757.191932.96162263581.5850.203803.380.02991.651885.831903.101733240.7320.219955.330.03051.512029.811942.88181061441.3590.207113.730.03011.521916.521912.8619951271.3380.200774.340.03001.601867.391913.012061731.1940.223865.840.02991.5220510.911902.84

图5 蛇盘兔花岗岩与龙须门闪长岩锆石U-Pb年龄谐和图

表2冀东地区早侏罗世岩浆岩锆石原位Hf同位素组成

Table 2 The in-situ Hf isotopic compositions of zircons for the Early Jurassic igneous rocks from eastern Hebei Province

测点号176Lu/177Hf2σ176Hf/177Hf2σ(176Hf/177Hf)iεHf(t)tCDM (Ma)JN09114蛇盘兔花岗岩10.003790.0000460.2824830.0000180.282470-6.5163230.002060.0000090.2824470.0000180.282439-7.6171640.003950.0000380.2824100.0000210.282395-9.1180160.001530.0000160.2824120.0000190.282407-8.7179480.001400.0000400.2824100.0000190.282405-8.8179990.001610.0000240.2824410.0000170.282435-7.71729100.003670.0000460.2824490.0000190.282436-7.71711110.004170.0000630.2824590.0000220.282444-7.41688120.002910.0000140.2824630.0000180.282452-7.11679140.003410.0001720.2824370.0000190.282425-8.11738160.005100.0001000.2824540.0000150.282435-7.71700180.003990.0000530.2824250.0000170.282410-8.61766190.004060.0000750.2824210.0000150.282407-8.71774200.001330.0000180.2824290.0000180.282424-8.11756210.006260.0000280.2824590.0000170.282437-7.71687JN09116龙须门闪长岩10.000730.0000090.2824920.0000160.282490-5.8161320.004590.0001000.2825500.0000220.282533-4.3148130.002480.0000670.2825050.0000190.282496-5.6158350.001360.0000790.2824870.0000190.282482-6.1162460.000740.0000030.2825200.0000180.282518-4.8154870.000890.0000230.2825570.0000190.282554-3.5146380.001220.0000100.2824990.0000170.282495-5.6159790.000820.0000100.2825140.0000180.282511-5.01562100.000690.0000070.2824990.0000160.282496-5.61597110.001370.0000110.2825190.0000180.282514-4.91551120.001920.0000210.2824440.0000190.282437-7.71723130.001150.0000130.2824920.0000180.282488-5.91613140.000880.0000020.2824820.0000160.282479-6.21636150.000690.0000190.2824710.0000170.282468-6.61662160.003380.0000410.2825650.0000220.282553-3.61445170.001500.0000570.2824830.0000200.282477-6.31635180.002420.0000690.2824850.0000230.282477-6.31628190.004110.0000500.2825190.0000230.282504-5.31552200.001410.0000630.2825050.0000180.282500-5.51584

图6 冀东地区早侏罗世侵入岩地球化学分类图解

3.2 全岩主微量元素

蛇盘兔花岗岩和龙须门闪长岩的稀土与微量元素组成明显不同(表4、图8、图9)。稀土元素配分曲线显示，蛇盘兔花岗岩显示不同程度的中稀土元素(MREE)亏损，富集K、Rb但亏损Ba、Sr、P和Ti，具有不同程度的负铕异常(Eu/Eu*=0.18～0.58)，REE含量也有较大的变化。龙须门闪长岩显示出右倾的平滑稀土配分型式，富集LILEs(Rb、K、Sr和Ba)，亏损HFSE(Nb、Ta、Th和U)，无或仅显示轻微的铕正异常(Eu/Eu*=1.0～1.1)。

4 讨论

4.1 早侏罗世岩石的识别

冀东地区中生代岩浆活动强烈，形成了晚三叠世、侏罗纪、早白垩世等不同时代的大量中酸性侵入岩。其中，已经识别出的早侏罗世岩体包括青山口花岗岩(199±2Ma; Miaoetal., 2008)、罗文峪花岗岩(197±7.0Ma; 陆继龙等, 2012)、高家店岩体(199±1Ma; 尹业长等, 2018)。另外，在肖营子镇以西的四拨子-六拨子矿区，也发现有少量的早侏罗世小型侵入单元，例如老商家和五拨子花岗斑岩的锆石年龄分别为190±0.7Ma和196±0.8Ma(李强等, 2012)。本文进一步获得青龙地区的蛇盘兔花岗岩和龙须门闪长岩的锆石年龄分别为192.0±2.4Ma和190.5±1.3Ma，表明冀东地区早侏罗世是一重要的岩浆作用时期。

同时，蛇盘兔花岗岩和龙须门闪长岩的早侏罗世锆石年龄与它们之南白家店花岗岩的中侏罗世年龄(约170Ma)明显不同(Fanetal., 2017)，它们原来被当作统一的燕山期肖营子岩体，因此，肖营子岩体应为两个时代形成的复式岩体。

4.2 早侏罗世岩石成因

表3冀东地区早侏罗世岩浆岩全岩主量元素(wt%)分析结果

Table 3 Whole-rock major elements (wt%) of Early Jurassic igneous rocks from eastern Hebei Province

样品号JN09114QL24QL25QL1401QL1405QL1406JD050∗JD052∗JN09116QL35QL1403岩性蛇盘兔花岗岩龙须门闪长岩SiO278.1477.7672.7372.5872.0776.2969.7563.154.1451.7950.69TiO20.090.10.30.320.380.270.600.641.491.921.90Al2O312.1312.2713.9914.8514.6712.9214.316.1217.3717.2817.43Fe2OT30.70.321.81.652.31.644.575.828.569.669.93MnO0.040.010.050.030.050.030.080.070.110.130.12MgO0.050.060.40.150.450.150.730.483.543.423.39CaO0.250.250.930.971.020.491.220.655.335.836.23Na2O3.923.94.4654.855.264.663.85.024.62.76K2O4.634.574.964.384.423.394.555.152.322.864.32P2O50.020.010.090.100.110.060.140.250.951.181.30LOI0.410.30.30.90.280.580.590.860.70.88Total100.3899.55100.01100.93100.60101.08101.1996.0899.6999.3798.95Mg#12.427.130.615.327.915.324.014.045.041.240.3A/CNK1.021.040.971.011.000.980.971.240.80.80.8

注：*引自尹业长等(2018)

图7 冀东地区早侏罗世侵入岩Harker图解

图8 冀东地区蛇盘兔花岗岩与龙须门闪长岩球粒陨石标准化稀土元素配分图(标准化值据Sun and McDonough, 1989)

图9 冀东地区蛇盘兔花岗岩与龙须门闪长岩原始地幔标准化微量元素蛛网图(标准化值据Sun and McDonough, 1989)

表4冀东地区早侏罗世岩浆岩全岩微量元素(×10-6)分析结果

Table 4 Whole-rock trace elements (×10-6) of Early Jurassic igneous rocks from eastern Hebei Province

样品号JN09114QL24QL25JD050∗JN09116QL35岩性蛇盘兔花岗岩龙须门闪长岩V3.691.6813.6116125Cr0.620.501.400.981.00Co0.290.081.7117.717.3Ni0.390.080.830.610.47Cu1.558.741.232.962.87Zn10.55.1030.284.2103Ga15.516.016.720.121.7Rb13514799.286.040.653.7Sr19.517.017119715831369Y17.710.624.127.620.826.2Zr10776.3202208114Nb25.023.419.931.811.716.2Ba36.225.664770311611573La25.86.3945.076.128.849.8Ce55.77.3489.315662103Pr4.971.019.7716.97.8811.8Nd14.82.9433.559.633.842.1Sm2.550.485.659.316.496.69Eu0.140.060.921.402.132.02Gd2.130.444.208.195.625.65Tb0.370.10.681.160.750.76Dy2.441.063.956.064.164.65Ho0.550.280.771.130.770.90Er1.801.262.443.272.062.40Tm0.290.220.370.520.280.33Yb2.161.762.383.381.742.00Lu0.370.270.360.490.270.28Hf4.183.265.062.324.582.41Ta2.171.751.462.830.630.70Pb24.421.0521.410.010.8Th11.66.869.1113.81.621.73U1.824.721.5325.60.380.48

注：*引自尹业长等(2018)

图10 冀东地区早侏罗世岩石锆石Hf同位素组成与该地区其它中生代岩石的对比图

然而，与蛇盘兔花岗岩相比，龙须门闪长岩在锆石Hf同位素上具有一定的差别。龙须门闪长岩锆石εHf(t)值在-7.7～-3.5，略高于蛇盘兔花岗岩锆石εHf(t)的-9.1～-6.5，表明随着岩浆的结晶分异，可能存在少量围岩的同化混染。

在讨论早侏罗世岩浆岩的岩浆来源时，需要与该地区不同时期岩浆岩进行对比。在青龙地区的晚三叠世都山花岗岩和中侏罗世白家店花岗岩中都存在暗色包体(叶浩等, 2014; Fanetal., 2017)。中侏罗世白家店花岗岩中的暗色包体(SiO2=52%)与寄主花岗岩之间在锆石年龄上完全相同，但是在锆石Hf同位素上存在明显差别，暗色包体比花岗岩具有更高的锆石εHf(t)和更大的εHf(t)变化范围(图10)。暗色包体的锆石εHf(t)为-3.7～+3.5，正的εHf(t)表明亏损地幔物质的加入，而其大的εHf(t)变化范围则表明壳幔岩浆的混合，而寄主花岗岩的锆石εHf(t)多数落入华北古老下地壳的范围，表明主要来源于古老下地壳的部分熔融。晚三叠世都山岩体主体为花岗岩(SiO2=70%～73%)，其锆石εHf(t)为-13.7～-9.6，而1个采自该岩基边部的闪长岩(SiO2=54%)样品在锆石年龄上和花岗岩完全一致，均为221Ma，但是其锆石εHf(t)为-5.2～+3.9(叶浩等, 2014)，与白家店暗色包体相似，其正的εHf(t)表明亏损地幔物质的加入，而其大的εHf(t)变化范围也表明壳幔岩浆的混合。

对比可以发现，尽管龙须门闪长岩与上述都山岩体和白家店岩体中的闪长岩三者在锆石年龄上明显不同，但它们在稀土配分模式以及微量元素配分模式上非常相似(图8、图9)。在锆石εHf(t)上，龙须门闪长岩明显低于都山和白家店闪长岩的锆石εHf(t)(图10)，但高于都山花岗岩和白家店花岗岩的锆石εHf(t)，都山花岗岩和白家店花岗岩的锆石εHf(t)主体上落在华北古老下地壳的范围之内，反映它们主要来源于古老下地壳部分熔融。尽管龙须门闪长岩在主、微量元素上与晚三叠世都山及中侏罗世白家店的闪长岩相似，但龙须门闪长岩的锆石Hf同位素变化范围很小，难以用壳幔岩浆混合来解释。另外，龙须门闪长岩具有极低的Cr、Ni含量(均≤1×10-6)，也排除具有高Cr、Ni含量幔源岩浆的加入。

蛇盘兔花岗岩在锆石εHf(t)上只是略低于龙须门闪长岩，两者的锆石εHf(t)显著高于华北古老下地壳的范围(Jiangetal., 2013)，说明它们可能并不是主要由古老下地壳重熔产生。考虑到该地区在三叠纪都山花岗岩体重的闪长岩中记录有亏损地幔物质的加入，我们认为在该地区在三叠纪可能存在亏损地幔物质的底侵作用，底侵作用形成的新生地壳与古老地壳在早侏罗世发生部分熔融，形成的岩浆在侵位前就已混合，具有比较均一的、高于古老下地壳部分熔融而低于新生地壳部分熔融岩浆的εHf(t)特征，蛇盘兔花岗岩比龙须门闪长岩稍低的锆石εHf(t)，表明随着岩浆的结晶分异，前者可能存在少量围岩的同化混染。

4.3 早侏罗世岩浆岩与区域演化的联系

华北克拉通在显生宙经历了岩石圈减薄和克拉通破坏，通常认为破坏的峰期在白垩纪(吴福元等, 2008)。但是一些研究显示，华北克拉通破坏可能在晚三叠世就已开始(Zhangetal., 2014)，并且始于周缘，这主要是基于华北北缘和辽东的一些晚三叠世岩浆岩中显示出软流圈地幔物质的加入(Yangetal., 2008; Zhangetal., 2014)得出的认识。

在冀东地区，晚三叠世都山岩体中存在亏损地幔物质的加入，进一步支持华北克拉通破坏可能在晚三叠世就已开始。我们的研究表明，青龙地区的早侏罗世岩浆岩可能来源于三叠纪底侵作用形成的新生地壳与古老地壳的部分熔融，反映了早侏罗世岩浆作用与华北克拉通破坏息息相关，早期的破坏为早侏罗世岩浆岩提供了所需的部分物质基础。另外，在冀东地区的中侏罗世岩浆岩中也存在亏损地幔物质的加入，同样反映了克拉通的破坏。这样直到华北克拉通破坏在早白垩世达到峰期，冀东地区的中生代岩浆作用表明克拉通的破坏是一个持续的长期过程。

5 结论

(1)锆石U-Pb定年结果显示，冀东青龙地区的蛇盘兔花岗岩和龙须门闪长岩分别形成于192.0±2.4Ma和190.5±1.3Ma，结合区内前人报道的多个同期岩体，表明冀东地区早侏罗世存在一期十分明显的岩浆作用，这一期岩浆活动在华北其它地区鲜有报道。

(2)岩石学、地球化学证据显示，冀东早侏罗世岩石主要由经历了不同程度结晶分异的岩浆岩组成，其岩浆可能来源于底侵作用形成的新生地壳与古老地壳的部分熔融。

(3)冀东地区中生代不同时期的岩浆岩，都与华北克拉通破坏有关，反映了华北克拉通破坏是一个持续的长期过程。

致谢作者感谢中国科学院地质与地球物理研究所杨岳衡、谢烈文、刘宇、唐国强、闫欣等老师在样品分析测试过程中提供的帮助与指导。感谢审稿人提出的宝贵的修改意见和建议。