岩石学报  2019, Vol. 35 Issue (4): 989-1014, doi: 10.18654/1000-0569/2019.04.02   PDF    
中生代华北克拉通岩石圈减薄的证据——以河北武安固镇杂岩体为例
霍延安 , 苏尚国 , 杨誉博 , 顾大鹏     
中国地质大学(北京)地球科学与资源学院, 北京 100083
摘要:本文对位于华北克拉通中部的固镇侵入杂岩体(闪长岩、二长岩、石英二长岩和正长岩)进行了野外地质调查、岩相学、矿物学、地球化学、锆石U-Pb定年、Sr-Nd同位素和Lu-Hf同位素研究,进而探讨其形成时代、源区和成因,以期得到关于华北克拉通中部中生代岩石圈减薄研究的新启示。研究区岩体形成时代和岩石结构有一定关联性,总体来看,等粒结构的岩体形成时代早于似斑状结构岩体,研究区岩浆活动可分为两个期次:第一期次岩浆活动形成等粒结构的深成岩体,包括中细粒闪长岩、二长岩、石英二长岩和正长岩,形成时代约为130~135Ma;第二期次岩浆活动形成浅部的似斑状岩体,包括似斑状闪长岩和石英二长岩,形成时代约为127Ma。结合野外观察、岩相学、矿物学、地球化学等特征,讨论分析了固镇杂岩体的成因:固镇闪长岩富镁,高Mg#值,Cr、Ni含量较高,富集LILE和LREE,亏损HFSE,铕异常不明显,具有高镁安山岩(HMA)性质,是与其相对应的侵入岩。固镇高镁闪长岩是被俯冲板片熔体和流体交代后的地幔橄榄岩部分熔融的产物,岩浆在上升途中与地壳发生了较弱的同化混染。固镇正长岩形成于加厚陆壳的部分熔融,地壳厚度约为40~67km。固镇二长岩是由区内正长岩和高镁闪长岩混合形成。研究区石英二长岩是区内高镁闪长岩与另一种更富硅的岩浆混合形成。俯冲板片物质(流体+熔体)对固镇高镁闪长岩的形成有重要贡献,显示出古太平洋板块俯冲作用的影响也涉及到了NCC中部,但受影响的程度小于NCC东部;正长岩形成于加厚陆壳部分熔融,显示出拆沉作用在中生代NCC中部岩石圈减薄过程中发挥了重要作用。结合角闪石温压计的计算结果和前人研究成果可以推测,NCC中部大规模拆沉作用发生的时间节点为127Ma。
关键词: 固镇侵入杂岩体     岩石成因     华北克拉通岩石圈减薄     古太平洋板块俯冲     拆沉作用    
The evidence for lithospheric thinning of Mesozoic North China Craton: Taking Guzhen intrusive complex as an example
HUO YanAn, SU ShangGuo, YANG YuBo, GU DaPeng     
School of Geosciences and Resources, China University of Geosciences, Beijing 100083, China
Abstract: This paper takes the Guzhen intrusive complex which is located in Hebei Province in the central part of the North China Craton (NCC) as a research object. The author conducted detailed field geological survey, combined with indoor research, including petrography, mineralogy, geochemistry, Sr-Nd isotope, Lu-Hf isotope and U-Pb zircon chronology analysis, so as to explore the formation age, source area, petrogenesis of Guzhen intrusive complex, in order to draw new inspiration for the Mesozoic lithospheric thinning in the central NCC. In this paper, it is believed that the formation age of the Guzhen complex is related to the structure of the rock. The formation age of the rock with equal grain structure is earlier than that of the rock with porphyritic structure. Magmatic activity can be divided into two phases:the first stage (130~135Ma) magmatic activity forms plutonic bodies with equal grain structure, including medium-fine diorite, monzonite, quartz monzonite and syenite; the second stage (127Ma) magmatism forms hypabyssal bodies, including porphyritic diorite and porphyritic quartz monzonite. Combined with field observation, petrography, mineralogy and geochemical characteristics, the petrogenesis of the Guzhen complex is analyzed. The Guzheng diorite has HMA properties. Guzhen diorite is a product of partial melting of the mantle peridotite metasomatized by the subduction slab melt and fluid, and the magma has weak assimilation and contamination with the crustal material on the ascending route. The Guzhen syenite is a product of partial melting of the thickened continental crust with the depth of 40~67km. The Guzhen monzonite is formed by the mixing of syenite and diorite in the study area. The Guzhen quartz monzonite is formed by mixing of Guzhen diorite and another silicon-rich magma. The subduction material (fluid+melt) has an important contribution to the formation of Guzheng high-Mg diorite, which shows that the influence of the subduction of the Paleo-Pacific plate also involves the central NCC, but the degree of the impact is less than that of the eastern NCC. The Guzhen syenite is formed by partial melting of the thickening continental crust, which shows that the delamination plays an important role in the thinning process of the NCC lithosphere in the Middle Mesozoic. Combined with the calculation results of the amphibole thermobarometer and the previous research results, it can be inferred that the time node of the large-scale delaminatiom of the central NCC is 127Ma.
Key words: Guzhen intrusive complex     Petrogenesis     Lithosphere thinning of the NCC     Paleo-Pacific plate subduction     Delamination    

克拉通是地球上长期稳定的构造单元,约占陆地面积的50%。典型的克拉通岩石圈厚约为200km,具有低密度、低热流值和高刚性的特点,所以后期的地质改造对其影响很小(Pearson, 1999; 朱日祥等, 2011)。华北克拉通(North China Craton)于18亿年左右形成统一的结晶基底,至早中生代一直持续稳定的克拉通演化阶段(Gao et al., 2002; 赵国春等, 2002)。但自中生代以来,NCC发生了大规模的岩浆活动以及构造变形,并且丢失了约120km的岩石圈地幔(Menzies et al., 1993; Fan et al., 2000)。目前,关于NCC岩石圈减薄主导机制存在着不同的观点(吴福元等,2008),主要为:岩石圈拆沉(Gao et al., 2004; 邓晋福等, 2006; Deng et al., 2007; Menzies et al., 2007; Liu et al., 2008; 许文良等, 2009)、热机械-化学侵蚀(Xu et al., 2004; Zhang et al., 2012; 李洪颜等, 2013)、地幔对流失稳(朱日祥等, 2011)和岩石圈地幔水化作用(Niu et al., 2005, 2015)。岩石圈减薄的时代为晚三叠世-早白垩世(Zhang et al., 2009Yang et al., 2007; Wu et al., 2005; 朱日祥等, 2011; 李洪颜等, 2013; Niu et al., 2015; )。对于NCC岩石圈减薄的动力学背景也存在着不同的观点,但是古太平洋板块俯冲的重要影响已经被广泛接受。研究区河北省武安固镇位于NCC中部,就NCC中部来说:前人研究认为其岩石圈减薄机制主要为拆沉作用(邓晋福等, 2006许文良等,2009王玉,2018霍腾飞等,2016)和热侵蚀作用(Xu et al., 2004; Zhang et al., 2012; 李洪颜等, 2013);关于古太平洋板块俯冲的影响是否涉及到了中部,也存在着不同的见解(罗照华等, 1997; 陈斌等, 2005a, b; Niu et al., 2015; 孙熠, 2016)。那么导致NCC中部中生代岩石圈减薄的主导机制到底是拆沉还是热侵蚀?古太平洋板块俯冲的影响是否涉及到了中部?本文对固镇侵入杂岩体进行了野外地质调查、岩石学、矿物学、锆石U-Pb年代学、全岩地球化学和Lu-Hf同位素研究,详细讨论了各岩性的形成时代、源区及成因,以期得到关于NCC中部中生代岩石圈减薄的新启示。

1 区域地质背景

研究区武安位于华北克拉通中部的太行山南段,处于赞皇隆起带以南,东临太行山山前断裂,西接西柏峪深大断裂,是邯邢地区的重要成矿地带(图 1)。该地区基底地层为太古宇赞皇群,上覆地层由老至新分别为:中元古代长城系,古生代寒武系、奥陶系、石炭系、二叠系,中生代三叠系、白垩系,新生代第三系、第四系,与基底地层呈角度不整合接触。研究区内广泛发育燕山期岩浆岩,出露面积约为350km2,岩体受北北东向深大断裂控制分为三个岩浆岩带,呈串珠状分布,分别为:东部岩浆岩带的新城岩体、洪山岩体、鼓山岩体,中部岩浆岩带的綦村岩体、矿山村岩体、武安岩体、固镇岩体,以及西部岩浆岩带的符山岩体。东部岩浆岩带主要岩性为正长岩,中部岩浆岩带为闪长岩-二长岩系列,西部岩浆岩带则为辉长岩-闪长岩系列,岩石自东向西基性程度逐渐增加(罗照华等, 1997陈斌等,2005b)。

图 1 华北克拉通区域地质图(a)和研究区区域地质图(b) Fig. 1 Regional geological map of the North China Craton (a) and regional geological map of study area (b)
2 岩体地质和岩相学特征

固镇侵入杂岩体位于武安岩体南部,岩体附近地层为奥陶系马家沟组灰岩及第四系沉积物,侵入岩可分为闪长岩、二长岩、石英二长岩和正长岩。结合野外穿插关系和测年数据可以总结出侵入关系由新至老为似斑状闪长岩→中细粒闪长岩→中细粒二长岩→中细粒正长岩。似斑状石英二长岩侵入到中细粒石英二长岩中。岩体间常见渐变过渡,部分岩体间接触界线截然。隐爆角砾岩呈条带状、不规则状分布在杂岩体边缘(图 2)。

图 2 固镇侵入杂岩体地质图 Fig. 2 Geological maps of Guzhen intrusive complex

野外观察闪长岩体中流体脉发育(图 3a),部分地区分布密集,形成“流体脉墙”,暗示形成闪长岩的岩浆为一个富水的体系。等粒结构的闪长岩为灰白色,中细粒结构,块状构造,主要矿物为角闪石(35%~45%)和斜长石(45%~55%),云母、绿帘石、磁铁矿及其他矿物含量约为5%~10%左右。似斑状闪长岩的斑晶为角闪石,多呈长柱状和六边形状,粒径约为3~7mm,含量约为5%~8%;基质主要为斜长石、角闪石、云母、碱性长石及少量其他矿物。镜下观察角闪石斑晶多呈长柱状和六边形状(图 4i),部分角闪石周围有磁铁矿圈。基质中的斜长石环带发育较少。

图 3 固镇侵入杂岩体野外特征 (a)闪长岩中的热液脉;(b)二长岩中的暗色微粒包体;(c)石英二长岩中的石英眼斑;(d)二长岩中发育反向脉的暗色微粒包体;(e)石英二长岩中的暗色微粒包体;(f)石英二长岩中的正长质“团块” Fig. 3 Field characteristics of Guzhen intrusive complex

图 4 固镇侵入杂岩体镜下特征 (a)石英二长岩中的暗色微粒包体(MME);(b)石英二长岩MME中的短柱状磷灰石与针状磷灰石共生;(c)斜长石周围发育细粒角闪石;(d)二长岩中的斜长石环带;(e)二长岩中辉石、石英和云母共生;(f)二长岩中的MME;(g)二长岩中发育针状磷灰石的斜长石核部;(h)石英二长岩中的暗色矿物聚集生长;(i)似斑状闪长岩中的角闪石斑晶 Fig. 4 Microscope photos of Guzhen intrusive complex

二长岩体中常见有不均匀现象,正长质岩浆呈团块状、阴影状、米级条带状存在于二长岩体中;暗色微粒包体也是十分常见的(图 3b),部分包体存在反向脉的切割(图 3d),即寄主二长岩质细脉穿切进暗色微粒包体的现象。二长岩,灰白色,中细粒结构,块状构造,主要矿物为斜长石、钾长石和角闪石。斜长石粒度约为2~3mm,含量约为40%~45%;角闪石粒度约为3~4mm,含量约为5%~10%;钾长石粒度约为2~3mm,含量约为35%~40%;石英含量约4%;副矿物为绿泥石、磁铁矿、榍石等,含量约为5%~7%。斜长石普遍发育聚片双晶和斜长石环带(图 4d),另发现部分斜长石核部发育针状磷灰石(图 4g),后被不具磷灰石的外圈包裹。二长岩中可见石英与辉石共生(图 4e)、角闪石包裹长石等现象。暗色微粒包体(图 4f)主要由粒度极细的角闪石和云母组成。

研究区石英二长岩暗色微粒包体发育,部分暗色微粒包体被寄主石英二长岩所晶出的长石切穿(图 3e)。不均匀现象也十分常见,例如石英二长岩体中存在不规则状的正长岩团块(图 3f)。石英二长岩中还可见石英眼斑,即石英斑晶周围镶有一圈暗色矿物微粒(图 3c)。似斑状石英二长岩斑晶为角闪石,自形长柱状,粒度约为3~5mm,含量约为5%~8%,基质主要为钾长石(35%~40%)、斜长石(35%~40%)、石英(5%~10%)和少量绿帘石、磁铁矿等副矿物(5%~7%)。中细粒石英二长岩,灰白色,块状构造,主要矿物为斜长石、钾长石和角闪石。斜长石粒度约为2~3mm,含量约为40%~45%;角闪石粒度约为3~4mm,含量约为5%~10%;钾长石粒度约为2~3mm,含量约为35%~40%;石英含量约5%~10%;副矿物为磁铁矿、磷灰石等,含量约为5%~7%。石英二长岩中存在暗色微粒包体(图 4a),并且有短柱状磷灰石和针状磷灰石在其中共存(图 4b);斜长石环带发育,部分斜长石具有钠化边,呈现“膜长石”特征;部分斜长石晶形较为浑圆,周围有细粒的暗色矿物(角闪石)围绕生长(图 4c)。

正长岩,浅肉红色,中细粒结构,块状构造,主要矿物为钾长石、斜长石和角闪石。钾长石粒径约为0.5~2.5mm,含量约为70%;斜长石粒度约为0.5~2.5mm,含量约为10%;角闪石粒径约为0.5~2.5mm,含量约为9%;石英含量约5%;其他矿物约占7%。镜下观察发现部分碱性长石具有绢云母化和高岭土化现象,斜长石可见显著的聚片双晶现象。

3 样品分析测试方法

矿物主量元素分析在中国冶金地质总局山东测试中心完成,测试仪器型号为:日本电子[JEOL]JXA-8230探针仪。分析条件:15kV电压,20nA电流,40°检出角,分析束斑:1~5μm,校正方法为ZAF法校正。

全岩主微量分析在河北廊坊市更新地质有限公司完成。主量元素检测标准为GB/T14506.28-2010,使用检测仪器为AxiosmaxX射线荧光光谱仪;FeO含量为滴定法测得;其灼失量为DZG20-1,H2O-为GB/T14506.1-2010,H2O+仪器为GB/T14506.2-2010,检测仪器为P1245电子分析天平。全岩微量元素分析采用ICP-MS溶液法分析。详细样品分析过程以及ICP测试条件见Qi et al.(2000)

本文部分锆石U-Pb测年在北京离子探针中心SHRIMP Ⅱ上完成。测试时一次流O2-强度为3~5nA,测试束斑直径为25~30μm。标样为M257 (U=840×10-6)和TEM年龄为417Ma(Black et al., 2003), 分别用于锆石U含量和年龄校正。另一部分锆石采用了LA-MC-ICP-MS锆石U-Pb同位素定年分析,在内蒙古自治区地质调查院分析测试中心完成,测试所用仪器为多接收电感耦合等离子体质谱仪(Neptune Plus)和激光剥蚀系统(GeoLasHD193nm)。GJ-1、Plešovice标准锆石为本次实验的锆石标样。

全岩Sr-Nd同位素分析在中国地质大学(北京)科学研究院同位素实验室完成。首先称量0.1~0.2g全岩粉末样品,使用混合酸(HF+HNO3+HClO4)溶样,待样品完全溶解后蒸干,加入6mol/L的HCl溶液,将样品转化为氯化物蒸干。然后用0.5mol/L的HCl溶液溶解,使用Rb-Sr、Sm-Nd交换柱进行元素分离提纯和元素提取。最后用ISOPROBE-T热电离质谱仪完成元素比值测定。

锆石原位Lu-Hf同位素测试分析工作在内蒙古自治区地质调查院分析测试中心完成。测试仪器为为多接收电感耦合等离子体质谱仪(Neptune Plus)和激光剥蚀系统(GeoLasHD193nm)。实验中激光剥蚀物质的载气是He,剥蚀直径为32μm,标准锆石为GJ-1和Plešovice。

4 锆石年代学

本文针对固镇杂岩体的不同岩性,包括闪长岩、二长岩、石英二长岩和正长岩,进行了锆石U-Pb定年。测试结果详见表 1表 2

表 1 固镇杂岩体SHRIMP锆石U-Pb定年数据 Table 1 SHRIMP zircon U-Pb data of Guzhen intrusive complex

表 2 固镇杂岩体激光拉曼锆石U-Pb定年数据 Table 2 LA-MC-ICP-MS zircon U-Pb data of Guzhen intrusive complex

闪长岩样品(B6368-2和B7909)的锆石自形程度良好,以短柱状为主,表面光滑(图 5)。其中B7909为中细粒闪长岩,Th/U值为0.30~0.67(Hermann et al., 2001),结合锆石特征,显示为岩浆锆石,计算的206Pb/238U加权平均年龄为131.7±1.2Ma(图 6)。B6368-2为似斑状闪长岩,Th/U值为1.20~3.19,显示出岩浆锆石特征(Hermann et al., 2001),计算的206Pb/238U加权平均年龄为126.7±1.1Ma (图 6)。

图 5 锆石阴极发光照片 Fig. 5 CL images of zircons

中细粒二长岩样品(B7909-1)的锆石形态以短柱状为主,并且发育有较为显著的震荡环带(图 5),Th/U值较高为0.81~6.94,计算的206Pb/238U加权平均年龄为133.7±3.2Ma(图 6)。

图 6 固镇侵入杂岩体U-Pb年龄协和图 Fig. 6 U-Pb age concordia diagram of Guzhen intrusive complex

石英二长岩样品(B7824和HY1-3A)的锆石自形程度良好,以短柱状为主,发育震荡环带(图 5)。HY1-3A为中细粒石英二长岩,Th/U值为0.30~0.83,计算的206Pb/238U加权平均年龄为129.5±1.1Ma(图 6);B7824为似斑状石英二长岩,Th/U值较高为1.28~3.41,计算的206Pb/238U加权平均年龄为126.6±1.3Ma(图 6)。

正长岩样品(ZZ0726-2和PM14-5b1)的锆石自形程度较差(图 5),Th/U比值低,显示继承锆石特征(Hermann et al., 2001),年龄集中在2300~2700Ma,可能为岩浆源区的继承或捕获锆石。以上2个样品测试结果协和度不高,未生成有效协和图。有研究者对区内的正长岩进行测年,测得其形成时代为135.6±1.1 Ma(王玉,2018),本文采用此年龄来进行后续讨论。

5 地球化学特征 5.1 全岩主量元素

固镇侵入杂岩体的主量元素测试结果见表 3

表 3 固镇杂岩体的主量元素(wt%)和稀土微量元素(×10-6)组成 Table 3 Major (wt%) and trace elements (×10-6) compositions for Guzhen intrusive complex

闪长岩SiO2含量为52.9%~54.0%,Al2O3含量为14.3%~15.9%,K2O+Na2O含量为6.3%~6.4%,CaO含量为6.2%~7.0%,MgO含量为6.0%~8.4%。在TAS图解中(图 7b)数据点均落入二长闪长岩区。在SiO2-K2O图解中(图 7a),样品点均落入高钾钙碱性区域。

图 7 固镇侵入杂岩体SiO2-K2O图解(a)和TAS图解(b) Fig. 7 SiO2 vs. K2O diagram (a) and TAS diagram (b) of Guzhen intrusive complex

二长岩SiO2含量为56.7%~59.8%,Al2O3含量较高为14.9%~16.9%。K2O+Na2O含量为6.8%~7.9%。CaO含量为4.7%~5.6%。MgO含量为2.7%~5.5%。A/CNK值平均为1.0~1.1,表明研究区二长岩为准铝质岩石。在TAS图解(图 7b)中落入二长岩区域。在SiO2-K2O图解中(图 7a),样品主要分布在高钾钙碱性系列,1个样品落入中钾钙碱性系列。

石英二长岩SiO2含量为62.7%~67.7%,Al2O3含量较高为15.9%~17.2%。K2O+Na2O含量为7.5%~9.5%。CaO含量为2.4%~4.2%。MgO含量为0.8%~1.6%。A/CNK值介于1.0~1.2,表明研究区石英二长岩为准铝质岩石。在TAS图解(图 7b)中落入石英二长岩区域。在SiO2-K2O图解中(图 7a),样品分布在高钾钙碱性系列。

正长岩SiO2含量介于59.9%~61.2%之间,Al2O3含量为16.6%~17.4%,K2O+Na2O含量为8.9%~9.9%,CaO含量为2.2%~3.4%,MgO含量为1.7%~2.9%。A/CNK值为0.9~1.1,表明研究区正长岩为准铝质岩石。在TAS图解中(图 7b)数据点均落入正长岩区域。在SiO2-K2O图解中(图 7a),样品主要分布于高钾钙碱性系列,1个样品位于中钾钙碱性系列。

5.2 全岩微量元素和稀土元素

微量元素与稀土元素分析测试结果见表 3

固镇杂岩体各岩性球粒陨石标准化稀土元素配分曲线皆呈现“右倾”,且样品均无明显Eu异常或具有微弱Eu正异常(图 8a-d)。原始地幔标准化微量元素蛛网图(图 8e-h)显示,固镇杂岩体各岩性均富集Rb、Ba、Sr等大离子亲石元素,亏损Nb、Ta、Ti等高场强元素,有较深的Nb-Ta槽。各岩性中,闪长岩的Cr、Ni含量远高于其他岩性。

图 8 固镇侵入杂岩体球粒陨石标准化稀土元素配分图(a-d, 标准化值据McDonough and Sun, 1995)和原始地幔标准化微量元素蛛网图(e-h, 标准化值据Sun and McDonough, 1989) Fig. 8 Chondrite-normalized REE patterns (a-d, normalization values after Mcdonough and Sun, 1995) and primitive mantle-normalized trace element patterns (e-h, normalization values after Sun and McDonough, 1989) of Guzhen intrusive complex
5.3 全岩Sr-Nd同位素

Sr-Nd同位素测试结果见表 4

表 4 固镇侵入杂岩体全岩Sr-Nd分析数据 Table 4 Whole-rock Sr and Nd isotopic data of Guzhen intrusive complex

对研究区固镇杂岩体的闪长岩、二长岩、石英二长岩和正长岩进行了Sr-Nd同位素测试,结果如下:闪长岩ISr=0.70580,εNd(t)=-16.3,tDM=2.20Ga;二长岩ISr=0.70602,εNd(t)=-17.0,tDM=2.19Ga;石英二长岩ISr=0.70570,εNd(t)=-13.9,tDM=1.95Ga;正长岩ISr=0.70575~0.70610,εNd(t)=-19.7~-17.6,tDM=2.03~2.14Ga。

5.4 锆石Lu-Hf同位素

在对固镇侵入杂岩体闪长岩、二长岩和石英二长岩锆石进行Lu-Hf同位素的分析。固镇侵入杂岩体锆石Lu-Hf同位素组成测试结果详见表 5

表 5 固镇侵入杂岩体锆石Lu-Hf同位素测试结果 Table 5 Zircon Lu and Hf isotopic data of Guzhen intrusive complex

固镇闪长岩(B7909)176Hf/177Hf值为0.282783~0.283013,εHf(t)值为2.84~10.70,tDM1值为425.3~748.3Ma,tDM2值为535~1109.9Ma。闪长岩中未发现老的继承锆石。

固镇二长岩(B7907)根据锆石年龄可分为两组:第一组年龄为130.4~147.0Ma,176Hf/177Hf值介于0.282561~0.282679之间。εHf(t)值介于-4.90~-0.98之间。tDM1值为948.2~1111.7Ma,tDM2值为1367.6~1661.5Ma。第二组年龄为1893.0~2524.0Ma,176Hf/177Hf值介于0.281271~0.281758之间。εHf(t)值介于1.45~13.65之间,tDM1值为2123.9~2710.3Ma,tDM2值为2084.9~2939.5Ma。

固镇石英二长岩(HY16和HY1-3): 176Hf/177Hf值介于0.282173~0.28284之间,εHf(t)值介于-18.58~4.41之间,tDM1值为729.6~1583.3Ma,tDM2值为964.1~2673.4Ma。

有研究者对本区正长岩进行了Lu-Hf同位素组成测试,176Hf/177Hf值介于0.28217~0.28245,εHf(t)值介于-18.4~-14.0,tDM1值为1318.1~1580.4Ma,tDM2值为1930.2~2671.8Ma(蒋俊毅,未刊数据)。

5.5 矿物化学特征

角闪石和斜长石电子探针分析结果见表 6表 7

表 6 固镇侵入杂岩体角闪石探针分析结果(wt%) Table 6 Major elements compositions (wt%) of amphiboles for Guzhen intrusive complex

表 7 固镇侵入杂岩体斜长石探针数据分析结果(wt%) Table 7 Major elements compositions (wt%) of feldspars for Guzhen intrusive complex

本文根据Ridolfi et al. (2010)提出的钙质角闪石温压计计算公式,对研究区闪长岩和石英二长岩中的角闪石进行了一系列计算(表 6),所得结果如下:中细粒石英二长岩中角闪石结晶压力为129.3~180.4MPa,对应的结晶深度为5~6km;角闪石结晶时的熔体含水量为3.94%~4.87%。似斑状石英二长岩中角闪石斑晶核部结晶压力为323.0~324.8MPa,对应的结晶深度为12.20~12.27km,熔体含水量为5.05%~5.77%;角闪石斑晶边部和基质角闪石的结晶压力为54.58~81.13MPa,对应的结晶深度为2.06~3.06km,熔体含水量为3.68%~3.83%。似斑状闪长岩中角闪石斑晶核部结晶压力为419.4~449.6MPa,对应的结晶深度为15.84~17.01km,熔体含水量为6.12%~56.27%;角闪石斑晶边部和基质角闪石的结晶压力为65.4~81.9MPa,对应的结晶深度为2.47~3.10km,熔体含水量为4.52%~4.92%。

研究区二长岩和石英二长岩中均发育斜长石的异常环带(表 7)。二长岩中发育反环带的斜长石:核部An牌号介于35~37之间,边部突然增长至47~60>45。石英二长岩中部分斜长石环带为正环带,但是从核部到边部An牌号呈陡然降低的趋势,核部An牌号介于50~53之间,边部陡然降低至33~36;另有部分斜长石为反环带,核部An=23~24,边部An牌号升至47~52。

6 讨论 6.1 岩体的形成时代

固镇杂岩体中细粒闪长岩和似斑状闪长岩的206Pb/238U加权平均年龄分别为131.7±1.2Ma(图 6)和126.7±1.1Ma (图 6);中细粒二长岩的206Pb/238U加权平均年龄为133.7±3.2Ma;中细粒石英二长岩和似斑状石英二长岩的206Pb/238U加权平均年龄分别为129.5±1.1Ma(图 6)和126.6±1.3Ma(图 6);前人对区内的中细粒正长岩进行测年,测得其形成时代为135.6±1.1Ma(王玉,2018)。

根据以上测试结果,本文认为研究区岩体形成时代和岩石结构有一定关联性,总体看,等粒结构岩体的形成时代早于似斑状结构岩体,岩浆活动可分为两个期次:第一期次岩浆活动形成等粒结构的深成岩体,包括中细粒闪长岩、中细粒二长岩、中细粒石英二长岩和中细粒正长岩,形成时代约为130~135Ma;第二期次岩浆活动形成浅成的似斑状岩体,包括似斑状闪长岩、似斑状石英二长岩,形成时代约为127Ma。

6.2 岩石成因及源区

Sr-Nd同位素通常被认为可以指示岩石的源区。本文中所有样品在ISr-εNd(t)图解中(图 9),均落在壳幔混合区域,且εNd(t)值相差不大。Leisher(1990)通过实验证明,当有不同物质组分的端元加入岩浆时,同位素的交换明显的快于元素的交换和混合;不相容元素的交换快于主量元素。Sr-Nd同位素的交换是在瞬间完成的。是以本文认为,研究区所有样品εNd(t)特征所能明确的仅仅是固镇杂岩体各岩性的Sr、Nd同位素已然达到了混染后十分均一的状态,表明固镇杂岩体经历了壳幔相互作用,对于原始源区的性质,已经不具备指示意义。而Hf同位素则不同,锆石的Hf同位素体系封闭温度很高,即使在高级变质条件下仍可保持最初的Hf同位素组成,使Hf同位素示踪具有一定的优越性,能更好的显示原始的源区性质(吴福元等,2007)。故而本文采用Hf同位素来进行有关岩石源区的讨论。

图 9 固镇侵入杂岩体εNd(t)-Isr图解(底图据张旗等,2008) Fig. 9 εNd(t) vs. Isr diagram for Guzhen intrusive complex (base map after Zhang et al., 2008)

La-La/Sm图解(图 10)可以指示岩浆过程中是否存在分离结晶作用,因为La为超岩浆元素,固-液相分配系数很低(与0.2~0.5比较可以忽略不计),在分离结晶过程中浓度变化缓慢,在部分熔融过程中浓度变化大;而Sm为亲岩浆元素,固-液相分配系数大于超岩浆元素,但与1比较可仍忽略不计,所以在分离结晶过程中同样浓度变化缓慢,而在部分熔融过程中浓度的变化小于La (韩吟文等, 2003)。如果本文岩浆经历了分离结晶过程,La-La/Sm图解中的轨迹应为图示水平虚线。在La-La/Sm图解中(图 10),研究区不同岩性的样品La与La/Sm为正相关关系,说明它们之间不具结晶分异关系。

图 10 固镇侵入杂岩体La-La/Sm图解 Fig. 10 La vs. La/Sm diagram for Guzhen intrusive complex

通过Sr-Nd同位素特征明确了固镇杂岩体各岩性都经历了壳幔相互作用,通过La-La/Sm图解明确了研究区不同岩性之间不具结晶分异关系,下面本文将分别讨论固镇杂岩体不同岩性的源区和成因。

6.2.1 闪长岩

高镁安山岩(High-Mg andesites)是指比典型的岛弧安山岩具有更高的MgO(>5%)和Mg#(>45),更低的Al2O3(<16%)和CaO(<10%)的安山岩,高镁闪长岩是与之对应的侵入岩(Tatsumi, 2001)。从地球化学特征来看,固镇闪长岩富镁(6.0%~8.4%),高Mg#值(56~67),Al2O3含量为14.3%~15.9%,CaO含量为6.2%~7.0%。从形成的构造环境来看,太行山造山带被认为是中生代环太平洋岩浆弧的组成部分(罗照华等,1997),138Ma左右古太平洋板块俯冲的弧岩浆前锋移至太行山(陈斌等,2005a),在微量元素构造判别图解(图 11)中所有闪长岩样品均落入火山弧花岗岩区域。以上特征说明固镇闪长岩属于高镁闪长岩(Tatsumi, 2001)。高镁闪长岩有可能是非高镁闪长岩质岩浆形成的堆晶岩,经观察,固镇闪长岩中并不存在堆晶结构,判断其不是堆晶岩,而是岩浆直接结晶产生的,可以用全岩化学成分对本区的高镁闪长岩成因展开进一步的讨论。

图 11 固镇闪长岩微量元素构造判别图解(据Pearce et al., 1996) Fig. 11 Tectonic setting discrimination diagram with trace elemental plots of Guzhen diorites(after Pearce et al., 1996)

关于高镁闪长岩的成因的观点,目前并不统一,有观点认为它的形成是俯冲板片洋壳部分熔融的熔体和/或脱水产生的流体与地幔楔交代的结果(Tatsumi et al., 1982; Kelemen et al., 1995 2003; Tatsumi, 1995);也有研究认为,富硅熔体与地幔橄榄岩相互作用也可以形成高镁闪长岩,富硅质熔体可能来自于地壳,也可能来自于俯冲板片的部分熔融(Xu et al., 2002, 2008; Gao et al., 2004)。相同的是,两种假说的主要源区都是交代地幔(Yin et al., 2010)。

固镇高镁闪长岩Mg#值较高(56~67),而地壳部分熔融形成的岩石Mg#值不会大于42(Wolf et al., 1992; Rapp et al., 1995)。εHf(t)值大于0(图 12),176Hf/177Hf介于0.28310~0.28297之间,与Schmidberger et al. (2002)所测得的加拿大Somerset岛橄榄岩176Hf/177Hf值相近也表明其幔源特征。同时固镇闪长又表现出富Na以及富含Cr、Ni、Co等过渡族元素的特征,Zr/Nb的值介于17~31之间,远远大于OIB的Zr/Nb值(5.83),以上特征均表明初始岩浆并非来自于软流圈地幔,而是源自岩石圈地幔。那么交代岩石圈地幔的物质是俯冲板片脱水流体?来自洋壳部分熔融产生的熔体?还是地壳部分熔融产生的富硅质熔体?须进行进一步讨论解决。

图 12 固镇侵入杂岩体U-Pb年龄与εHf(t)图解 符山辉长岩-闪长岩样品引自Chen et al., 2008;正长岩样品引自蒋俊毅,未刊数据 Fig. 12 U-Pb ages vs. εHf(t) diagram for Guzhen intrusive complex

华北克拉通已发现较多源自SCLM (Subcontinental Lithospheric mantle)的岩石,可以作为了解其化学成分的窗口。部分岩石起源于低金云母/单斜辉石比率的SCLM最上部(Schimidt et al., 1999; Guo et al., 2003),具有低的87Sr/86Sr(0.7039~0.7058)和低的143Nd/144Nd(0.5114~0.5123)(Guo et al., 2001),岩性主要为辉长岩和玄武岩,分布于太行地区、NCC北缘以及山东西部等地;部分岩石起源于受扬子克拉通俯冲改造的地幔橄榄岩(Zhang et al., 2002; Ying et al., 2004),具有较高的87Sr/86Sr(0.7061~0.7118)和低的143Nd/144Nd(0.5114~0.5121)(Zhang et al., 2002; Ying et al., 2004; Yang et al., 2004),岩性主要为玄武岩、安山岩碳酸岩等,分布在山东东部、山东西部等地。而本文中以及部分NCC其他地区(如鲁西等地)的高镁闪长岩与以上两种均有不同,具有中等的87Sr/86Sr(0.7053~0.7074)和较高的143Nd/144Nd(0.5117~0.51247)(Xu et al., 2004),表明源区地幔经历了不同的改造过程。

固镇高镁闪长岩中有大量的黑云母和角闪石,暗示岩浆源区存在较高的水逸度。本文利用Ridolfi et al. (2010)总结出的钙质角闪石温压计对固镇似斑状高镁闪长岩中的角闪石进行了相关计算,结果显示角闪石斑晶在16~17km深度结晶时,岩浆含水量高达6%左右,即使在浅部定位后,含水量也大于4%。而研究区西侧地幔源区受拆沉地壳影响的高镁闪长岩(许文良等,2009),用同样方法测得的岩浆含水量仅为2%左右,部分小于1%(孙熠,2016)。微量元素的比值(La/Yb、Th/Yb、Ba/La、La/Nb等)被广泛应用于识别含水流体和俯冲带沉积物(Hawkesworth et al., 1991; Pearce and Peate, 1995)。Ba比LREE更易溶于流体,另一方面Th和Yb在含水流体中不易迁移,因此当俯冲板片的流体加入到地幔楔源区中时,Ba/La值应升高,Th/Yb不变。本文中高镁闪长岩样品具有较低且集中的Th/Yb值和升高的Ba/La值,样品分布呈现出俯冲板片流体加入的趋势。以上特征均表明固镇高镁闪长岩的形成过程中可能有俯冲板片脱水产生的流体加入。且实验岩石学已证明,地幔橄榄岩的局部熔融只能产生玄武岩,然而在含水条件下,顽火辉石发生不一致熔融可以产生镁橄榄石和SiO2,从而形成富Si富Mg的高镁闪长岩(邓晋福等,2004),与本文讨论结果相一致。

研究区西侧的符山高镁闪长岩的地幔源区经历了拆沉下地壳的强烈改造(许文良等,2009),致使其εHf(t)值小于0,在Hf同位素图解中主要分布在-10~-20区间(图 12),与之相比,固镇高镁闪长岩的εHf(t)值则大于0。另外,固镇闪长岩具有更高的Mg#值(56~67)和过渡族元素(Cr、Ni、Co等)含量。同时值得注意的是,符山高镁闪长岩中存在大量的下地壳继承锆石,而研究区闪长岩中并未发现。以上特征说明固镇高镁闪长岩的地幔源区并没有经历地壳物质的强烈改造,其受地壳物质改造的程度远远小于符山高镁闪长岩。那么造成其明显富集LILE亏损HFSE的应当另有其他原因。

地幔橄榄岩中,各残留相的KdYbmin/liq远远小于1,使Yb表现为极不相容元素,典型的岛弧英安岩中Yb含量平均为4.4×10-6,而源区为地幔橄榄岩的固镇高镁闪长岩Yb含量较低(1.6×10-6~2.2×10-6),推测其亏损HREE是由于源区残留相中存在KdYbmin/liq>1的矿物,而石榴子石是地幔中能使岩浆富集LILE亏损HREE的矿物。已有研究表明,被俯冲板片熔体交代的地幔橄榄岩能大规模的转化为辉石+石榴子石+角闪石+金云母组合(Prouteau et al., 2001),此矿物组合在后续的熔融过程中,以3:1的比例产生石榴子石和单斜辉石(Francis and Ludden, 1995),这个过程使石榴子石得以在更大的熔融程度下保持稳定,进而产生像固镇高镁闪长岩一样的富集LILE亏损HREE的岩浆。同时,由于源区石榴子石在较大区间保持稳定,大量的Lu存留在石榴石中,而Hf则进入熔体相,导致熔体相具相对较高的176Hf/177Hf和较低的176Lu/177Hf,也可以解释研究区εHf(t)值较高的现象(Patchett et al., 1981; Bizimis et al., 2003; Schmitz et al., 2004)。另外,即使固镇高镁闪长岩显示出明显的Nb、Ta亏损,但其Nb的浓度(4.41×10-6~6.2×10-6)大于地幔Nb含量(通常<1×10-6),Nb在流体中不易迁移(Tatsumi et al., 1986),却能够跟随熔体转移,推测一定量的Nb元素跟随俯冲板片熔体进入了地幔橄榄岩。综上所述,固镇高镁闪长岩源区地幔橄榄岩可能经历了俯冲板片富硅熔体的交代。

根据上述讨论得出,固镇高镁闪长岩是被俯冲板片熔体和流体交代后的地幔橄榄岩部分熔融的产物,且高镁闪长质岩浆在上升途中与地壳发生了较弱的同化混染。

6.2.2 正长岩

邯邢地区正长岩出露将对较少,主要出露在研究区东侧的永安县洪山地区,武安地区也有少量出露,目前研究最多的为洪山正长岩体,认为其源区为EMI型富集地幔(陈斌等,2005b)。而对武安杂地壳的正长岩则缺少研究。那么武安的固镇正长岩具有怎样的成因呢?

研究区正长岩εHf(t)值介于-18.4~-14.0之间,显示壳源特征,不同于源区为EM1型富集地幔的洪山正长岩(εHf(t)值介于-7.0~-10.0之间, 陈斌等,2005b)。对固镇杂岩体样品进行SHRIMP定年,发现正长岩样品所测定的大部分锆石均为介于2300~2700Ma的继承锆石(表 1),也说明了固镇正长岩的源区可能为华北地台古老下地壳,幔源的洪山岩体则从未发现过大规模的继承锆石。实验岩石学已证明,地壳中玄武质岩石部分熔融形成花岗岩类Mg#值不大于42(Wolf et al., 1992; Rapp et al., 1995),而研究区部分正长岩样品Mg#较高(35.5~49.6),证明了幔源物质的贡献。推测是由于幔源岩浆的底侵作用(以固镇闪长岩为代表)导致了幔源物质的混入。

固镇正长岩稀土元素球粒陨石标准化图解中(图 8c)可见Eu几乎无异常或呈Eu正异常,表明源区无斜长石残留,重稀土元素的亏损,表明源区有石榴子石残留,根据残留相推断研究区高Sr低Yb型正长岩形成深度可能大于40km(张旗等, 2006)。即使有幔源组分的加入,本文中正长岩样品的Nb/La比值仍很低(平均0.2),等同于华北下地壳岩石Nb/La的下限,推测下地壳部分熔融可能发生于金红石稳定的条件,指示研究区正长岩起源于增厚地壳。另,Condie (1973)提出了安山质岩石SiO2为60%时K2O的质量分数K60与地壳厚度的关系,而后周珣若和吴克隆(1994)将其扩大到成分近似安山质的侵入岩中,证明具有较好的可移植性。基于K60与地壳厚度的依赖关系(邓晋福等,2004),可以从岩石系列得到地壳厚度的总体框架:低钾钙碱性系列地壳厚度<17km,中钾钙碱性系列地壳厚度为17~40km,高钾钙碱性系列地壳厚度为40~67km。在SiO2-K2O图解中(图 7a),固镇正长岩样品点主要落入高钾钙碱性系列,推测地壳厚度在40~67km之间。仅一个样品落入中钾钙碱性系列,此样品在野外和镜下均显示经历了强烈的钠化作用,钠化作用导致了K2O含量的降低。综上所述,固镇正长岩形成于加厚陆壳的部分熔融,地壳厚度约为40~67km。

6.2.3 二长岩和石英二长岩

二长岩以及石英二长岩在本文固镇研究区广泛出露。在野外地质勘查过程中,就发现了固镇二长岩和石英二长岩具有极其显著的岩浆混合成因的野外宏观特征,具体描述如下:

暗色微粒包体(Mafic microgranular enclave)是指分布于花岗质岩石之中,暗色、细粒的,且比寄主岩石更富铁镁矿物的岩石包体,是岩浆混合最重要的特征(王德滋等,2008)。研究区二长岩、石英二长岩中普遍发育闪长质暗色微粒包体(MME),呈水滴状、椭圆状、不规则状等形态(图 3b, d, e图 4a, f),与寄主岩体接触边界多不截然,表明包体进入寄主岩浆中时,二者均为塑性的或处于流动状态。部分包体发育反向脉,这是由于淬冷导致MME突然冷凝收缩形成了不规则裂缝,使寄主岩浆得以注入到包体之中(图 3d)。二长岩和石英二长岩体中混合导致的成分不均匀现象也十分常见(Barharin, 1999; 王涛, 2000),如:正长岩呈团块状包卷在石英二长岩体中;正长岩呈阴影状、不规则状赋存与二长岩中(图 3f);石英二长岩岩体中有米级正长岩条带等,体现出不同端元岩浆共存时的塑性状态,甚至是紊流状态。区内二长岩、石英二长岩内亦可见石英斑晶外圈被细小的暗色矿物镶边(图 3c),这种现象通常称为石英眼斑,也是一种典型的岩浆混合现象,当偏基性岩浆中已经结晶的石英由于岩浆混合而进入偏基性端元时,作为捕掳晶的石英会发生融解变得浑圆,并吸收岩浆中的结晶潜热,使其晶体周围形成一个局部过冷的环境,致使暗色矿物围绕其结晶。以上现象的普遍存在,使我们不得不考虑,研究区的二长岩以及石英二长岩是岩浆混合的产物。

研究区的二长岩以及石英二长岩的岩浆混合成因也得到了地球化学和同位素的支持。在Harker图解中(图 13),固镇二长岩与区内高镁闪长岩和正长岩呈现出良好的成分过渡关系,在Hf同位素图解中(图 12),二长岩样品分布在闪长岩和正长岩之间,结合野外特征可以推测二长岩是由区内正长岩和闪长岩混合形成的,二者分别为酸性端元和基性端元。与固镇二长岩不同,研究区石英二长岩的SiO2含量大于正长岩,很显然不可能是正长岩与闪长岩混合的产物。同时对石英二长岩中的基性包体做了地球化学测试,结果显示其与区内高镁闪长岩具有相似的地球化学性质(图 8a, e),故而推测石英二长岩的基性端元与二长岩相同,均为区内闪长岩,其酸性端元可能是一种比区内正长岩更富Si的岩浆,在混合过程中已经耗尽或者未出露地表。在Hf同位素图解中(图 12),石英二长岩的样品分布在较大的区间,也可以用岩浆混合作用来解释(吴福元等,2007)。

图 13 固镇侵入杂岩Harker图解 Fig. 13 Harker diagrams of Guzhen intrusive complex

除野外宏观特征和地球化学特征以外,矿物的显微构造和化学成分的变化能更好的为岩浆混合作用的发生提供有力证据。镜下观察二长岩和石英二长岩中普遍发育具有岩浆混合特征的显微构造,如:暗色微粒包体以及其中共存的短柱状磷灰石和针状磷灰石(图 4a, b)(Wyllie et al., 1962);辉石、石英和云母共生的不平衡矿物共生现象(图 4e);斜长石核部发育针状磷灰石后被无磷灰石的边部包裹(图 4g)(Wyllie et al., 1962);暗色矿物聚集生长(图 4h)等。另外,斜长石的异常环带也普遍存在于二长岩和石英二长岩之中,斜长石存在两种异常环带:第一种就是我们通常说得“反环带”——斜长石核部An牌号低为,二长岩集中在35~37,石英二长岩集中在23~33,外圈被富Ca的长石包裹,An牌号增至>45;第二种异常环带从核部到边部整体看来是An牌号降低的正环带,但是这种牌号的变化并不是渐变的,而是突变的,通常是核部的牌号为50左右,到边部陡然降低到<40,最低可达<30。一般认为,当结晶环境较为稳定和封闭时,会形成斜长石的正环带,或者岩浆房内的对流将斜长石带到岩浆房的不同部位导致结晶条件的微弱变化会导致斜长石具有震荡环带(Viccaro et al., 2010),但总体来说,斜长石An牌号的变化幅度较小,不会超过10。本文中二长岩和石英二长岩的斜长石呈现的An牌号异常,暗示晶体形成于开放的岩浆系统,有基性岩浆的注入(Tepley et al., 2000)。

6.3 对华北克拉通中部中生代岩石圈减薄的指示意义

部分研究认为,古太平洋板块俯冲的影响没有涉及到华北克拉通中部,因为华北克拉通中部距海岸线距离过远,超过了1000km。并且有研究表明,NCC中部的下部为贫水地幔(夏群科等,2017),另外,在山东铁铜沟发育有直接源自俯冲板片熔融的花岗岩(Martin et al., 2005),而研究区则未见出露。然而有地球物理数据显示,克拉通下滞留在地幔过渡带的古太平洋板块属于水平俯冲(陈斌等,2005a),板块水平俯冲的影响可以远至2000km的距离,如果将中生代以来中国东部巨大的伸展量以及新生代形成的日本海的距离去除,太行山当时距离古太平洋俯冲带的距离应该小于2000km(陈斌等,2005b)。本文研究的固镇高镁闪长岩是被俯冲板片熔体和流体交代后的地幔橄榄岩部分熔融的产物,是古太平洋板块俯冲影响涉及到了NCC中部的印证,本文未发育类似山东铁铜沟直接由俯冲板片熔融形成的花岗岩(Martin et al., 2005),可能是由于体量较少,在交代岩石圈地幔时已经耗尽。综上所述,本文认为古太平洋板块俯冲作用的影响也涉及到了NCC中部,但受影响的程度小于NCC东部。

Kay et al.(1993)提出,当地壳厚度>50km时,下地壳可向榴辉岩相转变,使下地壳密度大于地幔,进而导致下地壳的拆沉作用。经上文讨论,固镇正长岩形成于加厚陆壳的部分熔融,形成深度约为40~67km,显然研究区的地壳已具备拆沉作用发生的厚度。同时,下地壳熔出正长质岩浆后,留下比重较大的榴辉岩相,会导致重力失稳而促进下地壳的拆沉。不仅本文研究区,华北克拉通中部的其他地区也存在加厚陆壳部分熔融所形成的岩体,例如,武安地区西石门正长岩(135Ma)(王玉, 2018)、山西省运城的蚕坊岩体和孤峰岩体(130Ma)(齐玥等, 2016)等。以上区域加厚地壳熔出富含易熔组分的岩浆后,都会导致下地壳向榴辉岩相的转变,进而很可能产生拆沉作用的发生。并且在华北克拉通中部也确实有多处成因与拆沉作用相关的岩体,且形成时代集中在126~130Ma之间,正是既加厚陆壳熔融成因的岩体形成时代之后,例如位于华北克拉通南部临汾市万荣岩体、万峰山岩体、狐偃山岩体和黄龙垴岩体等,均形成于127~130Ma之间(霍腾飞等, 2016);太行山南段的符山高镁闪长岩(126Ma)中包含有被地壳物质改造的幔源包体,被认为是地壳拆沉导致的结果(许文良等, 2009)。

那么是否可以对NCC中部拆沉作用发生的时间做出制约呢?研究区固镇杂岩体存在中细粒结构岩体与似斑状结构岩体密切共生的现象,对固镇杂岩体似斑状闪长岩、似斑状石英二长岩和中细粒石英二长岩中的角闪石进行全铝温压计计算(Ridolfi et al., 2010),从计算结果(表 6)可知:中细粒石英二长岩中角闪石结晶深度,即等粒结构岩体的就位深度为5~6km。似斑状闪长岩中角闪石核部的结晶深度为16~17km,角闪石斑晶边部和基质角闪石的结晶深度为3km(似斑状闪长岩定位深度);似斑状石英二长岩角闪石斑晶核部的结晶深度为12km,角闪石斑晶边部和基质角闪石的结晶深度为2~3km(似斑状石英二长岩定位深度)。综合研究区锆石U-Pb测年数据结果(表 1表 2):等粒结构岩体的形成时代集中在130~135Ma,形成最晚的等粒结构的岩体(中细粒石英二长岩)年龄约为130Ma;随后,似斑状结构的岩体集中形成于127Ma。也就是说,在大约3Myr之间,地壳的隆升剥蚀达2~4km。这说明研究区的此次的地壳隆升剥蚀是快速的、剧烈的、突发的,这很有可能与拆沉作用有关。拆沉作用具有“灾变性”(邓晋福等, 2006),引起地壳减薄和隆升的过程是突然的,而不是缓慢的,与本文研究区所体现的隆升过程一致。研究区似斑状岩体的快速上升侵位均发生在127Ma,暗示发生拆沉作用的时间节点在127Ma左右。总结前人的研究资料发现:NCC中部中生代岩浆活动主要集中在138~124Ma,源自加厚地壳部分熔融或富集地幔的岩体形成时代区间为138~130Ma,源自加厚陆壳拆沉或经历了快速上升侵位的岩体形成时代集中在127Ma左右,个别为130Ma(蔡剑辉等,2003许文良等,2004陈斌等, 2005a, b苏尚国等,2006刘玲等,2009霍腾飞等,2016齐玥等,2016;),与本文所得的127Ma时间节点相一致。推测NCC中部在127Ma爆发了大规模的拆沉作用,130Ma可能为此次大规模拆沉的前奏。

综上所述,古太平洋板块俯冲作用的影响也涉及到了NCC中部,但受影响的程度小于NCC东部;拆沉作用在NCC中部岩石圈减薄过程中发挥重要作用,大规模拆沉作用发生的时间节点为127Ma。

7 结论

(1) 研究区岩体形成时代和岩石结构有一定关联性,总体来看,等粒结构的岩体形成时代早于似斑状结构岩体,研究区岩浆活动可分为两个期次:第一期次岩浆活动形成等粒结构的深成岩体,包括中细粒闪长岩、二长岩、石英二长岩和正长岩,形成时代约为130~135Ma;第二期次岩浆活动形成浅部的似斑状岩体,包括似斑状闪长岩和石英二长岩,形成时代约为127Ma。

(2) 结合野外观察、岩相学、矿物学、地球化学等特征,讨论分析了固镇杂岩体的成因:固镇闪长为高镁闪长岩,是被俯冲板片熔体和流体交代后的地幔橄榄岩部分熔融的产物,岩浆在上升途中与地壳发生了较弱的同化混染。固镇正长岩形成于加厚陆壳的部分熔融,地壳厚度约为40~67km。固镇二长岩是由区内正长岩和闪长岩混合形成。研究区石英二长岩是区内闪长岩与另一种更富硅的岩浆混合形成。

(3) 俯冲板片物质(流体+熔体)对固镇高镁闪长岩的形成有重要贡献,显示出古太平洋板块俯冲作用的影响也涉及到了NCC中部,但受影响的程度小于NCC东部;正长岩形成于加厚陆壳部分熔融,显示出拆沉作用在NCC中部岩石圈减薄过程中发挥重要作用。结合角闪石温压计计算结果和前人研究成果可以推测,NCC中部大规模拆沉作用发生的时间节点为127Ma。

致谢      感谢邓晋福老师、罗照华老师在成文过程中以及野外工作时的教导与帮助;同时感谢师门的同学们在野外工作过程中的帮助!

参考文献
Barbarin B. 1999. A review of the relationships between granitoid types, their origins and their geodynamic environments. Lithos, 46: 605-626. DOI:10.1016/S0024-4937(98)00085-1
Black LP, Kamo SL, Allen CM, Aleinikoff JN, Davis DW, Korsch RJ and Foudoulis C. 2003. TEMORA 1:A new zircon standard for Phanerozoic U-Pb geochronology. Chemical Geology, 200(1-2): 155-170. DOI:10.1016/S0009-2541(03)00165-7
Cai JH, Yan GH, Chang ZS, Wang XF, Shao HX and Chu ZY. 2003. Petrological and geochemical characteristic of the Wanganzhen complex and discussion on its genesis. Acta Petrologica Sinica, 19(1): 81-92.
Chen B, Tian W, Zhai MG and Arakawa Y. 2005a. Zircon U-Pb geochronology and geochemistry of the Mesozoic magmatism in the Taihang Mountains and other places of the North China craton, with implications for petrogenesis and geodynamic setting. Acta Petrologica Sinica, 21(1): 13-24.
Chen B, Zhai MG, Tian W and Jiang BM. 2005b. Petrogenesis of the Mesozoic intrusive complexes from the southern Taihang Orogen, North China Craton:Elemental and Sr-Nd-Pb isotopic constraints. Bulletin of Mineralogy, Petrology and Geochemistry, 24(2): 93-102.
Chen LH and Zhou XH. 2005. Subduction-related metasomatism in the thinning Lithosphere:Evidence from a composite dunite-orthopyroxenite xenolith entrained in Mesozoic Laiwu high-Mg diorite, North China Craton. Geochemistry, Geophysics, Geosystems, 6(6): 1-20.
Condie KC. 1973. Archean magmatism and crustal thickening. Geological Society of America Bulletin, 84: 2981-1992. DOI:10.1130/0016-7606(1973)84<2981:AMACT>2.0.CO;2
Deng JF, Luo ZH, Su SG, Mo XX, Yu BS, Lai XY and Chen HW. 2004. Rock Genesis, Tectonic Setting and Mineralization. Beijing: Geological Publishing House, 26-28 (in Chinese)
Deng JF, Su SG, Liu C, Zhao GC, Zhao XG, Zhou S and Wu ZX. 2006. Discussion on the Lithospheric thinning of the North China craton:Delamination? or thermal erosion and chemical metasomatism?. Earth Science Frontiers, 13(2): 105-119.
Deng JF, Su SG, Niu YL, Liu C, Zhao GC, Zhao XG, Zhou S and Wu ZX. 2007. A possible model for the lithospheric thinning of North China Craton:Evidence from the Yanshanian (Jura-Cretaceous) magmatism and tectonic deformation. Lithos, 96: 22-35. DOI:10.1016/j.lithos.2006.09.009
Fan WM, Zhang HF, Baker J, Jarvis KE, Mason PRD and Menzies MA. 2000. On and off the North China Craton:Where is the Archaean keel. Journal of Petrology, 41: 933-950. DOI:10.1093/petrology/41.7.933
Francis D and Ludden J. 1995. The signature of amphibole in mafic alkaline lavas, a study in the Northern Canadian Cordillera. Journal of Petrology, 36: 1171-1191. DOI:10.1093/petrology/36.5.1171
Gao S, Rudnick RL, Carlson RW, Mcdonough WF and Liu YS. 2002. Re-Os evidence for replacement of ancient mantle Lithosphere beneath the North China Craton. Earth and Planetary Science Letters, 198: 307-322. DOI:10.1016/S0012-821X(02)00489-2
Gao S, Rudnick RL, Yuan HL, Liu XM, Liu YS, Xu WL, Ling WL, Ayers J, Wang XC and Wang QH. 2004. Recycling lower continental crust in the North China Craton. Nature, 432: 892-897. DOI:10.1038/nature03162
Guo F, Fan WM, Wang YJ and Lin G. 2001. Late Mesozoic mafic intrusive complexes in North China Block:Constraints on the nature of subcontinental lithospheric mantle. Physics and Chemistry of the Earth Part A Solid Earth and Geodesy, 26(9-10): 759-771. DOI:10.1016/S1464-1895(01)00125-9
Guo F, Fan W, Wang YJ and Lin G. 2003. Geochemistry of Late Mesozoic mafic magmatism in West Shandong Province, eastern China:Characterizing the lost lithospheric mantle beneath the North China Block. Geochemical Journal, 37: 63-77. DOI:10.2343/geochemj.37.63
Han YW and Ma ZD. 2003. Geochemistry. Beijing: Geological Publishing House, 203-204 (in Chinese)
Hawkesworth CJ, Hergt JM, Ellam RM and McDermott F. 1991. Element fluxes associated with subduction related magmatism. Philosophical Transactions of the Royal Society of London. Series A:Physical and Engineering Sciences, 335: 393-405. DOI:10.1098/rsta.1991.0054
Hermann J, Rubatto D, Korsakov A and Shatsky VS. 2001. Multiple zircon growth during fast exhumation of diamondiferous, deeply subducted continental crust (Kokchetav massif, Kazakhstan). Contributions to Mineralogy and Petrology, 141(1): 66-82. DOI:10.1007/s004100000218
Huo TF, Yang DB, Shi JP, Xu WL and Yang HT. 2016. Petrogenesis of the Early Cretaceous alkali-rich intrusive rocks in the central North China Block:Constraints from zircon U-Pb chonology and Sr-Nd-Hf isotopes. Acta Petrologica Sinica, 32(3): 697-712.
Kay RW and Kay SM. 1993. Delamination and delamination magmatism. Tectonophys, 219: 177-189. DOI:10.1016/0040-1951(93)90295-U
Kelemen PB, Yogodzzinski GM and Scholl DW. 2003. Along-strike variation in the Aleutian island arc:Genesis of high Mg# andesite and implications for continental crustal inside the subduction factory. American Geophysical Union Geophysical Monograph, 138: 223-276.
Kelemen PB. 1995. Genesis of high Mg-number andesites and the continental-crust. Contributions to Mineralogy and Petrology, 120(1): 1-19. DOI:10.1007/BF00311004
Lesher CE. 1990. Decoupling of chemical and isotopic exchange during magma mixing. Nature, 344: 235-237. DOI:10.1038/344235a0
Li CN and Liao QA. 2006. Magmatic mixing and mingling of Precambrian Gangbian complex in northeastern Jiangxi Province and their geological significance. Acta Petrologica et Minerlogica, 25(5): 357-376.
Li HY. 2013. Destruction of North China Craton:Insights from temporal and spatial evolution of the proto-basins and magmatism. Scientia Sinica (Terrae), 43(9): 1396-1409.
Liu L, Chen B and Liu AK. 2009. Petrogenesis of the Zijingguan mafic pluton, northern Taihang Orogen:Constraints from petrology and geochemistry. Earth Science, 34(1): 165-178.
Liu SC and Xia QK. 2014. Water content in the Early Cretaceous lithospheric mantle beneath the south-central Taihang Mountains:Implications for the destruction of the North China Craton. Chinese Science Bulletin, 59: 1280-1283. DOI:10.1007/s11434-014-0136-6
Liu Y, Gao S, Kelemen PB and Xu W. 2008. Recycled crust controls contrasting source compositions of Mesozoic and Cenozoic basalts in the North China Craton. Geochimica et Cosmochimica Acta, 72: 2349-2376. DOI:10.1016/j.gca.2008.02.018
Luo ZH, Deng JF, Zhap GC and Cao YQ. 1997. Magmatic activity characteristics of Taihangshan orogenic belt and its inversion of orogenic process. Earth Science, 22(3): 279-284.
Ma CQ, Wang RJ and Qiu JX. 1992. Signs of granitic magma origin and multiple magma mixing:Inclusions:Taking the Zhoukoudian rock mass in Beijing as an example. Geological Review, 38(2): 109-119.
Martin H, Smithies RH, Rapp R, Moyen JF and Champion D. 2005. An overview of adakite, tonalite-trondhjemite-granodiorite (TTG), and sanukitoid:Relationships and some implications for crustal evolution. Lithos, 79(1-2): 1-24. DOI:10.1016/j.lithos.2004.04.048
Mcdonough WF and Sun SS. 1995. The composition of the Earth. Chemical Geology, 120(3-4): 223-253. DOI:10.1016/0009-2541(94)00140-4
Menzies MA, Fan WM and Zhang M. 1993. Paleozoic and Cenozoic lithoprobe and the loss of >20km of Archean lithosphere, Sino-Korean Craton, China. Geological Society, London, Special Publication, 76: 71-78. DOI:10.1144/GSL.SP.1993.076.01.04
Menzies MA, Xu YG, Zhang H and Fan W. 2007. Integration of geology, geophysics and geochemistry:A key to understanding the North China Craton. Lithos, 96: 1-21. DOI:10.1016/j.lithos.2006.09.008
Niu YL. 2005. Generation and evolution of basaltic magmas:Some basic concepts and a hypothesis for the origin of the Mesozoic-Cenozoic volcanism in eastern China. Geological Journal of China Universities, 11: 9-46.
Niu YL, Liu Y, Xue QQ, Shao FL, Chen S, Duan M, Guo PY, Gong HM, Hu Y, Hu ZH, Kong JJ, Li JY, Liu JJ, Sun P, Sun WL, Ye L, Xiao YY and Zhang Y. 2015. Exotic origin of the Chinese continental shelf:New insights into the tectonic evolution of the western Pacific and eastern China since the Mesozoic. Earth Sciences, 60(18): 1598-1616.
Pearce JA and Peate DW. 1995. Tectonic implications of the composition of volcanic arc magmas. Annual Review of Earth and Planetary Sciences, 23: 251-285. DOI:10.1146/annurev.ea.23.050195.001343
Pearce JA. 1996. Sources and setting of granitic rocks. Episodes, 19(4): 120-125.
Pearson DG. 1999. The age of continental roots. Lithos, 48: 171-194. DOI:10.1016/S0024-4937(99)00026-2
Prouteau G, Scaillet B, Pichavant M and Maury R. 2001. Evidence for mantle metasomatism by hydrous silicic melts derived from subducted oceanic crust. Nature, 410: 197-200. DOI:10.1038/35065583
Qi L, Hu J and Gregoire DC. 2000. Determination of trace elements in granites by inductively coupled plasma-mass spectrometry. Talanta, 51(3): 507-513. DOI:10.1016/S0039-9140(99)00318-5
Qi Y, Luo JH, Wu JD and Chen FK. 2016. Geochemical and Sr-Nd-Pb isotopic composition of the Canfang and Gufengshan granodiorite plutons in central-southern North China. Acta Petrologica Sinica, 32(7): 2015-2028.
Rapp RP and Wastson EB. 1995. Dehydration melting of metabasalt at 8~32kbar:Implications for continental growth and crust-mantle recycling. Journal of Petrology, 36: 891-931. DOI:10.1093/petrology/36.4.891
Ridolfi F, Renzulli A and Puerini M. 2010. Stability and chemical equilibrium of amphibole in calc-alkaline magmas:An overview, new thermobarometric formulations and application to subduction-related volcanoes. Contributions to Mineralogy and Petrology, 160: 45-66. DOI:10.1007/s00410-009-0465-7
Schmidberger SS, Simonetti A, Francis D and Gariépy C. 2002. Probing Archean Lithosphere using the Lu-Hf isotope systematics of peridotite xenoliths from Somerest Island kimberlites, Canada. Earth and Planetary Science Letters, 197: 245-259. DOI:10.1016/S0012-821X(02)00491-0
Schmidt, Bottazzi P, Vannucci R and Mengel K. 1999. Trace element partitioning between phlogopite, clinopyroxene and leucite lamproite melt. Earth and Planetary Science Letters, 168: 287-299. DOI:10.1016/S0012-821X(99)00056-4
Su SG, Deng JF, Zhao GC, Zhao XG and Liu C. 2006. Characteristics, origin and resource property of Xuejiashiliang complex, Beijing area and their relationship with the way of lithospheric thining. Earth Science Frontiers, 13(2): 148-157.
Sun SS and McDonough WF. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. In: Sanders AD and Norry MJ (eds.). Magmatism in Ocean Basins. Geological Society, London, Special Publication, 42(1): 313-345
Sun Y. 2016. Petrogenesis and Mineralization of the Early Cretaceous intrusive rocks in Handan-Xingtai area, Hebei Province, China. Ph. D. Dissertation. Beijing: China University of Geosciences, 1-134 (in Chinese with English summary)
Tatsumi Y and Ishizaka K. 1982. Origin of high-magnesian andesites in the Setochi volcanic belt, Southwest Japan, Petrographical and chemical characteristics. Earth and Planetary Science Letters, 60(2): 293-304. DOI:10.1016/0012-821X(82)90008-5
Tatsumi Y and Nakamura N. 1986. Composition of aquaeous fluid from serpentinite in the subducted lithosphere. Geochemical Journal, 20: 191-196. DOI:10.2343/geochemj.20.191
Tatsumi Y. 1995. Subduction Zone Magmatism. Boston: Blackwell Publishing, 1-244
Tatsumi Y. 2001. Geochemical modeling of partial melting of subducting sediments and subsequent melt-mantle interaction:Generation of high-Mg andesites in the Setouchi volcanic belt, Southwest Japan. Geology, 29(4): 323-326. DOI:10.1130/0091-7613(2001)029<0323:GMOPMO>2.0.CO;2
Taylor SR and McLennan SM. 1985. The Continental Crust: Its Composition and Evolution. Oxford: Blackwell Scientific Publications, 1-312
Tepley FJ, Davidson JP, Tilling RI and Arth JG. 2000. Magma mixing, recharge and eruption histories recorded in plagioclase phenocrysts from EI Chichon Volcano, Mexico. Journal of Petrology, 41(9): 1397-1411. DOI:10.1093/petrology/41.9.1397
Viccaro M, Giacomoni PP, Ferlito C and Cristofolini R. 2010. Dynamics of magma supply at Mt. Etna volcano (Southern Italy) as revealed by textural and compositional features of plagioclase phenocrysts. Lithos, 116(1-2): 77-91.
Wang DZ and Xie L. 2008. Magma Mingling:Evidence from Enclaves. Geological Journal of China Universities, 14(1): 16-21.
Wang T. 2000. Origin of hybrid granitoids and the implications for continental dynamics. Acta Petrologica Sinica, 16(2): 161-168.
Wang XX, Wang T, Lu XX and Xiao QH. 2002. Petrographic evidence of the magmatic mixing of Laojunshan and Qinling Lianghuan granitoids in the North Qinling Mountains and its significance. Geological Bulletin of China, 24(8-9): 523-529.
Wang Y. 2018. Origin and formation mechanism of pyroxene syenite in Hongshan, Xishimen, Hebei Province. Master Degree Thesis. Beijing: China University of Geosciences, 1-83 (in Chinese with English summary)
Weaver BL and Tarney J. 1984. Empirical approach to estimating the composition of the continental crust. Nature, 310(5978): 575-577. DOI:10.1038/310575a0
Wolf MB and Wyllie PJ. 1992. The formation of tonalitic liquids during the vapor-absent partial melting of amphibolite at 10kbar. Eos, 70: 506-518.
Wu FY, Lin JQ, Wilde SA, Zhang XO and Yang JH. 2005. Nature and significance of the Early Cretaceous giant igneous event in eastern China. Earth and Planetary Science Letters, 233: 103-119. DOI:10.1016/j.epsl.2005.02.019
Wu FY, Li XH, Zheng YF and Gao S. 2007. Lu-Hf isotopic systematics and their applications in petrology. Acta Petroloica Sinica, 23(2): 185-220.
Wu FY, Xu YG, Gao S and Zheng JP. 2008. Lithospheric thining and destruction of North China Craton. Acta Petroloica Sinica, 24(6): 1145-1174.
Wyllie PJ, Cox KG and Biggar GM. 1962. The habit of apatite in synthetic systems and igneous rocks. Journal of Petrology, 3(2): 238-243. DOI:10.1093/petrology/3.2.238
Xia QK and Hao YT. 2013. The distribution of water in the continental Lithospheric mantle and its implications for the stability of continents. Chinese Science Bulletin, 58: 3879-3889. DOI:10.1007/s11434-013-5949-1
Xia QK, Cheng H and Liu J. 2017. The distribution of the Early Cretaceous hydrous lithophere mantla in the North China Ctaton:Constraints from water content in peridotites of Tietonggou. Earth Science, 42(6): 853-861.
Xu JF, Shinjo R, Defant MJ, Wang Q and Rapp RP. 2002. Origin of Mesozoic adakitic intrusive rocks in the Ningzhen area of East China:Partial melting of delaminated lower continental crust. Geology, 30(12): 1111-1114. DOI:10.1130/0091-7613(2002)030<1111:OOMAIR>2.0.CO;2
Xu W, Hergt JM, Gao S, Pei F, Wang W and Yang D. 2008. Interaction of adakitic melt-peridotite:Implications for the high-Mg# signature of Mesozoic adakitic rocks in the eastern China Craton. Earth and Planetary Science Letters, 265(1-2): 123-137. DOI:10.1016/j.epsl.2007.09.041
Xu WL, Wang QH, Wang DY, Pei FP and Gao S. 2004. Processes and mechanism of Mesozoic lithospheric thinning in eastern North China Craton:Evidence from Mesozoic igneous rocks and deep-seated xenoliths. Earth Science Frontiers, 11(3): 310-317.
Xu WL, Yang DB, Pei FP and Yu Y. 2009. Petrogenesis of Fushan high-Mg# diorites from the southern Taihang Mts. in the central North China Craton:Resulting from interaction of peridotite-melt derived from partial melting of delamination lower continental crust. Acta Petrologica Sinica, 25(8): 1947-1961.
Xu YG, Huang XL, Ma JL, Wang YB, Lizuka Y, Xu JF, Wang Q and Wu XY. 2004. Crust-mantle interaction during the tectono-thermal reactivation of the North China craton:Constraints from SHRIMP zircon U-Pb chronology and geochemistry of Mesozoic plutons from western Shandong. Contributions to Mineralogy and Petrology, 147: 750-767. DOI:10.1007/s00410-004-0594-y
Yang CH, Xu WL, Yang DB, Wang W, Wang WD and Liu JM. 2008. Petrogenesis of Shangyu gabbro-diorites in the west of Shandong:Evidence of chronology and geochemistry. Science in China (Series D), 38(1): 44-55.
Yang JH, Chung SL, Zhai MG and Zhou XH. 2004. Geochemical and Sr-Nd-Pb isotopic compositions of mafic dikes from the Jiaodong Peninsula, China:Evidence for vein-plus-peridotite melting in the lithospheric mantle. Lithos, 73: 145-160. DOI:10.1016/j.lithos.2003.12.003
Yang JH, Sun JF and Chen FK. 2007. Sources and petrogenesis of Late Triassic dolerite dikes in the Liaodong Peninsula:Implications for post-collisional Lithosphere thinning of eastern North China craton. Journal of Petrology, 48: 1973-1997. DOI:10.1093/petrology/egm046
Yin JY, Yuan C, Sun M, Long XP, Zhao GC, Wong KP, Geng HY and Cai KD. 2010. Late Carboniferous high-Mg dioritic dikes in western Junggar, NW China:Geochemical features, petrogenesis and tectonic implications. Gondwana Research, 17: 145-152. DOI:10.1016/j.gr.2009.05.011
Ying JF, Zhou XH and Zhang HF. 2004. Geochemical and isotopic investigation of the Laiwu-Zibo carbonatites from western Shandong Province, China, and implications for their petrogenesis and enriched mantle source. Lithos, 75: 413-426. DOI:10.1016/j.lithos.2004.04.037
Zhang HF, Sun M, Zhou XH, Zhou XH, Fan WM, Zhai MG and Yin JF. 2002. Mesozoic Lithosphere destruction beneath the North China Craton:Evidence from major, trace element, and Sr-Nd-Pb isotope studies of Fangcheng basalts. Contributions to Mineralogy and Petrology, 144: 241-253. DOI:10.1007/s00410-002-0395-0
Zhang HF, Sun YL, Tang YJ, Xiao Y, Zhang WH, Zhao XM, Santosh M and Menzies MA. 2012. Melt peridotite interaction in the Pre-Cambrian mantle beneath the western North China Craton:Petrology, geochemistry and Sr, Nd and Re isotopes. Lithos, 149: 100-114. DOI:10.1016/j.lithos.2012.01.027
Zhang Q, Wang Y, Li CD, Wang YL, Jin WJ and Jia XQ. 2006. Granite classification on the basis of Sr and Yb contents and its implications. Acta Petrologica Sinica, 22(9): 2249-2269.
Zhang Q, Wang Y, Xiong XL and Li CD. 2008. Adakite and Granite Challenge and Opportunity. Beijing: China Land Press: 107-129.
Zhang SH, Zhao Y, Liu XC, Liu DY, Chen FK, Xie LW and Chen HH. 2009. Late Paleozoic to Early Mesozoic mafic-ultramafic complexes from the northern North China block:Constraints on the composition and evolution of the lithospheric mantle. Lithos, 110: 229-246. DOI:10.1016/j.lithos.2009.01.008
Zhao GC, Sun M and Wilde SA. 2002. Characteristics of the basement tectonic units of the North China Craton and the Early Proterozoic colliding. Science in China (Series D), 32: 538-549.
Zhao Y, Chen B, Zhang SH, Liu JM, Hu JM, Liu J and Pei JL. 2010. Major geological events in the northern margin of the North China Craton and the adjacent area in Early Yanshan Period. Geology in China, 37(4): 900-915.
Zheng JP, Lu FX, O'Reilly SY, Griffin WL and Zhang M. 1999. Characteristics and evolution of the Paleozoic and Cenozoic lithospheric mantles in the eastern North China Platform. Acta Geologica Sinica, 73(1): 47-56. DOI:10.1111/acgs.1999.73.issue-1
Zhou R and Wu KL. 1994. Zhangzhou I-A type granite. Beijing: Science Press: 1-152.
Zhu RX, Chen L, Wu FY and Liu JL. 2011. Timing, scale and mechanism of the destruction of the North China Craton. Scientia Sinica (Terrae), 41(5): 583-592.
蔡剑辉, 闫国翰, 常兆山, 王晓芳, 邵宏翔, 储著银. 2003. 王安镇岩体岩石地球化学特征及成因探讨. 岩石学报, 19(1): 81-92.
陈斌, 田伟, 翟明国, 荒川洋二. 2005a. 太行山和华北其它地区中生代岩浆作用的锆石U-Pb年代学和地球化学特征及其岩浆成因和地球动力学意义. 岩石学报, 21(1): 13-24.
陈斌, 翟明国, 田伟, 江博明. 2005b. 太行山南段中生代杂岩体的岩石成因:元素和Nd-Sr-Pb同位素地球化学证据. 矿物岩石地球化学通报, 24(2): 93-102.
邓晋福, 罗照华, 苏尚国, 莫宣学, 于炳松, 赖兴运, 谌宏伟. 2004. 岩石成因、构造环境与成矿作用. 北京: 地质出版社: 26-28.
邓晋福, 苏尚国, 刘翠, 赵国春, 赵兴国, 周肃, 吴宗絮. 2006. 关于华北克拉通燕山期岩石圈减薄的机制与过程的讨论:是拆沉, 还是热侵蚀和化学交代?. 地学前缘, (2): 105-119. DOI:10.3321/j.issn:1005-2321.2006.02.009
韩吟文, 马振东. 2003. 地球化学. 北京: 地质出版社: 203-204.
霍腾飞, 杨德彬, 师江朋, 许文良, 杨浩田. 2016. 华北地块中部早白垩世富碱侵入岩的成因:锆石U-Pb年代学和Sr-Nd-Hf同位素制约. 岩石学报, 32(3): 697-712.
李昌年, 廖群安. 2006. 赣东北前寒武纪港边杂岩体的岩浆混合(和)作用及其地质意义. 岩石矿物学杂志, 25(5): 357-376. DOI:10.3969/j.issn.1000-6524.2006.05.001
李洪颜. 2013. 华北克拉通原型盆地及岩浆活动时空演化对克拉通破坏的制约. 中国科学(地球科学), 43(9): 1396-1409.
刘玲, 陈斌, 刘安坤. 2009. 北太行紫荆关基性岩体的成因:岩石学和地球化学的证据. 地球科学, 34(1): 165-178. DOI:10.3321/j.issn:1000-2383.2009.01.015
刘少辰, 夏群科. 2014. 太行山中南段早白垩世岩石圈地幔的水含量:对华北克拉通破坏的动力因素制约. 科学通报, 59: 1280-1283.
罗照华, 邓晋福, 赵国春, 曹永清. 1997. 太行山造山带岩浆活动特征及其造山过程反演. 地球科学, 22(3): 279-284.
马昌前, 王人镜, 邱家骧. 1992. 花岗质岩浆起源和多次岩浆混合的标志:包体——以北京周口店岩体为例. 地质论评, 38(2): 109-119. DOI:10.3321/j.issn:0371-5736.1992.02.002
齐玥, 罗金海, 巫嘉德, 陈福坤. 2016. 华北中南部蚕坊和孤峰山花岗闪长岩岩体的地球化学特征和Sr-Nd-Pb同位素组成. 岩石学报, 32(7): 2015-2028.
苏尚国, 邓晋福, 赵国春, 赵兴国, 刘翠. 2006. 北京燕山地区薛家石梁杂岩体特征、成因、源区性质及岩石圈减薄方式. 地学前缘, 13(2): 148-157. DOI:10.3321/j.issn:1005-2321.2006.02.013
孙熠. 2016.河北邯邢地区早白垩世岩浆作用及成矿意义.博士学位论文.北京: 中国地质大学, 1-134 http://cdmd.cnki.com.cn/Article/CDMD-10491-1016312065.htm
王德滋, 谢磊. 2008. 岩浆混合作用:来自岩石包体的证据. 高校地质学报, 14(1): 16-21. DOI:10.3969/j.issn.1006-7493.2008.01.002
王涛. 2000. 花岗岩混合成因研究及大陆动力学意义. 岩石学报, 16(2): 161-168.
王晓霞, 王涛, 卢欣祥, 肖庆辉. 2002. 北秦岭老君山、秦岭梁环斑结构花岗岩岩浆混合的岩相学证据及其意义. 地质通报, 24(8-9): 523-529.
王玉. 2018.河北武安西石门洪山辉石正长岩的成因和形成机制.硕士学位论文.北京: 中国地质大学, 1-83
吴福元, 李献华, 郑永飞, 高山. 2007. Lu-Hf同位素体系及其岩石学应用. 岩石学报, 23(2): 185-220.
吴福元, 徐义刚, 高山, 郑建平. 2008. 华北岩石圈减薄与克拉通破坏研究的主要学术争论. 岩石学报, 24(6): 1145-1174.
谢群科, 程徽, 刘佳. 2017. 山东铁铜沟橄榄岩的水含量:华北克拉通早白垩世富水岩石圈的分布. 地球科学, 42(6): 853-861.
许文良, 王清海, 王冬艳, 裴福萍, 高山. 2004. 华北克拉通东部中生代岩石圈减薄的过程与机制:中生代火成岩和深源捕虏体证据. 地学前缘, 11(3): 310-317.
许文良, 杨德彬, 裴福萍, 于洋. 2009. 太行山南段符山高镁闪长岩的成因——拆沉陆壳物质熔融的熔体与地幔橄榄岩反应的结果. 岩石学报, 25(8): 1947-1961.
杨承海, 许文良, 杨德彬, 王伟, 王伟德, 刘金民. 2008. 鲁西上峪辉长-闪长岩的成因:年代学与岩石地球化学证据. 中国科学(D辑), 38(1): 44-55. DOI:10.3321/j.issn:1006-9267.2008.01.005
张旗, 王焰, 李承东, 王元龙, 金惟俊, 贾秀勤. 2006. 花岗岩的Sr-Yb分类及其地质意义. 岩石学报, 22(9): 2249-2269.
张旗, 王焰, 熊小林, 李承东. 2008. 埃达克岩和花岗岩:挑战和机遇. 北京: 中国大地出版社: 107-129.
赵国春, 孙敏, Wilde SA. 2002. 华北克拉通基底构造单元特征及早元古代拼合. 中国科学(D辑), 32: 538-549.
赵越, 陈斌, 张拴宏, 刘建民, 胡健民, 刘健, 裴军令. 2010. 华北克拉通北缘及邻区前燕山期主要地质事件. 中国地质, 37(4): 900-915. DOI:10.3969/j.issn.1000-3657.2010.04.007
郑建平, 路凤祥, Reilly SYO, Griffin WL, 张明. 1999. 华北地台东部古生代与新生代岩石圈地幔特征及其演化. 地质学报, 73(1): 47-56. DOI:10.3321/j.issn:0001-5717.1999.01.006
周珣若, 吴克隆. 1994. 漳州I-A型花岗岩. 北京: 科学出版社: 1-152.
朱日祥, 陈凌, 吴福元, 刘俊来. 2011. 华北克拉通破坏的时间、范围与机制. 中国科学(地球科学), 41(5): 583-592.