岩石学报  2018, Vol. 34 Issue (11): 3327-3340   PDF    
鲁西早白垩世岩石圈地幔的属性:大昆仑辉长岩和辉绿岩年代学、岩石地球化学和Sr-Nd-Pb-Hf同位素制约
杨浩田1 , 杨德彬1,2,3 , 师江朋1 , 许文良1,2 , 王枫1,2 , 谌瑛1     
1. 吉林大学地球科学学院, 长春 130061;
2. 国土资源部东北亚矿产资源评价重点实验室, 长春 130061;
3. 东北亚生物演化与环境教育部重点实验室, 长春 130026
摘要:大昆仑辉长岩和辉绿岩锆石LA-ICP-MS U-Pb年代学、锆石Hf同位素和全岩地球化学及全岩Sr-Nd-Pb同位素的研究为探讨鲁西地区早白垩世岩石圈地幔的属性提供了制约。结果表明,大昆仑辉长岩中少量锆石和辉绿岩中锆石发育条痕状吸收,具有相对高的Th/U比值(0.11~3.12),暗示它们均为岩浆成因,辉长岩中多数为岩浆成因的捕获锆石。辉长岩和辉绿岩中岩浆锆石的206Pb/238U加权平均年龄分别为132±3Ma和112±1Ma;辉长岩中捕获锆石的谐和年龄主要介于2469~1752Ma之间,暗示古老地壳物质的存在。辉长岩和辉绿岩整体具有低SiO2(49.08%~52.61%)、高MgO(4.52%~9.67%)、Na2O含量(2.58%~3.49%)以及高Mg#值(45~65)特征;富集轻稀土元素和大离子亲石元素(Rb、Ba、Th)、贫重稀土元素和高场强元素(Nb、Ta、Zr和Hf),其(La/Yb)N=8.98~10.2,Eu负异常不明显(δEu=0.83~0.93)。辉长岩的初始87Sr/86Sr比值变化于0.7046~0.7060之间,εNdt)值介于-16.71~-15.54之间,206Pb/204Pb、207Pb/204Pb和208Pb/204Pb初始比值分别为16.596~16.776、15.255~15.290和36.328~36.648。辉绿岩中岩浆锆石的εHft)值介于-12.0~-2.5之间。上述结果表明,大昆仑辉长岩和辉绿岩形成于早白垩世,与华北克拉通东部早白垩世岩石圈减薄的峰期时间一致;起源于受华北克拉通古老下地壳物质改造的富集岩石圈地幔的部分熔融,形成于岩石圈减薄的伸展构造背景。
关键词: 辉长岩     辉绿岩     岩石圈地幔     岩石成因     早白垩世     鲁西    
Nature of the Early Cretaceous lithospheric mantle in western Shangdong: Constraints from geochronology, geochemistry and Sr-Nd-Pb-Hf isotopic data of Dakunlun gabbros and diabases
YANG HaoTian1, YANG DeBin1,2,3, SHI JiangPeng1, XU WenLiang1,2, WANG Feng1,2, CHEN Ying1     
1. College of Earth Sciences, Jilin University, Changchun 130061, China;
2. Key Laboratory of Mineral Resources Evaluation in Northeast Asia, Ministry of Land and Resources, Changchun 130061, China;
3. Key Laboratory for Evolution of Past Life and Environment in Northeast Asia, Ministry of Education, Changchun 130026, China
Abstract: The LA-ICP-MS zircon U-Pb dating, Hf isotope, whole-rock geochemical and Sr-Nd-Pb data are reported for the Dakunlun gabbros and diabases in western Shandong to understand their nature of the Early Cretaceous lithospheric mantle. Few zircons from Dakunlun gabbros and the zircons from diabases all display striped absorption in cathodoluminescence images. It also indicates a magmatic origin that the high Th/U ratio (0.11~3.12). Moreover, most of zircons in the gabbros are captured and magmatic origin. The weighted mean 206Pb/238U ages of the magmatic zircons are 132±3Ma and 112±1Ma in the gabbros and diabases, respectively, implying an Early Cretaceous in age. The concordia ages of the captured zircons range from 2469Ma to 1752Ma, suggesting the existence of ancient crustal material. Geochemically, the Dakunlun gabbros and diabases are characterized by low SiO2 (49.08%~52.61%), high MgO (4.52%~9.67%), Na2O (2.58%~3.49%) and high Mg# (45~65). Moreover, they are enriched in light rare earth elements and large ion lithophile elements (e.g., Rb, Ba, Th), and are depleted in heavy rare earth elements and high field strength elements (e.g., Nb, Ta, Zr, and Hf). They have (La/Yb)N=8.98~10.2 and weak negative Eu anomaly (δEu=0.83~0.93). The (87Sr/86Sr)i ratios vary between 0.7046 and 0.7060 and the εNd(t) vary between -16.71 and -15.54 in the gabbros. Whole-rock initial 206Pb/204Pb, 207Pb/204Pb and 208Pb/204Pb ratios of the gabbros range from 16.596 to 16.776, 15.255 to 15.290 and 36.328 to 36.648, respectively. The εHf(t) values range from -12.0 to -2.5 in the magmatic zircons from diabases. The above data show that the Early Cretaceous Dakunlun gabbros and diabases are consistent with the period of lithospheric thinning beneath the eastern North China Craton. The primary magma was derived by partial melting of the enriched lithospheric mantle modified by the ancient lower crustal material. They were formed in an extensional tectonic environment related to the lithospheric thinning.
Key words: Gabbro     Diabase     Lithospheric mantle     Petrogenesis     Early Cretaceous     Western Shandong Province    

鲁西地区位于华北克拉通东部,东与郯庐断裂带和苏鲁造山带相邻,该区发育有丰富的晚中生代镁铁质岩浆岩及其携带的地幔橄榄岩包体,成为研究华北克拉通东部中生代岩石圈地幔的性质以及揭示华北克拉通岩石圈地幔减薄过程和机制的天然实验室(Zhang et al., 2004; Xu et al., 2004, 2008, 2013; 裴福萍等, 2004; 杨承海等, 2005; Yang et al., 2005, 2008, 2010, 2012; Huang et al., 2012; 钟军伟和黄小龙, 2012; Li and Li, 2016)。近年来,鲁西地区晚中生代辉长岩-辉长闪长岩-玄武岩的研究已经取得了丰富的成果,如Yang et al. (2012)基于对鲁西地区典型镁铁质岩石的全岩Sr-Nd-Pb同位素和锆石Hf同位素的空间变异研究认为,靠近郯庐断裂带一侧的鲁西岩石圈地幔(方城、费县玄武岩)主要受到了俯冲扬子板片部分熔融形成的熔体的改造,而远离郯庐断裂带一侧的鲁西西北部的岩石圈地幔(济南、邹平辉长岩)主要为华北克拉通古老下地壳物质的改造,进而提出扬子克拉通深俯冲对华北克拉通东南缘岩石圈地幔改造的空间范围约200km。此外,Zhang et al. (2004)认为鲁西地区中生代岩石圈地幔具高度不均一性,济南、邹平和莱芜等地区的岩石圈地幔表现出EMⅠ特征,继承了古老岩石圈地幔的属性并且很少受到富硅质流体的交代。而Huang et al. (2012)认为济南和邹平地区的岩石圈地幔不能用简单的地壳混染进行解释,鲁西西北部的岩石圈地幔在古元古代时期可能已经受到了拆沉陆壳物质的影响,鲁西地区的早白垩世基性岩浆起源于混合的地幔源区。由此可见,关于鲁西地区早白垩世岩石圈地幔的属性和改造过程仍存在较大争议。鉴于此,本文选择研究程度较低的鲁西大昆仑辉长岩和辉绿岩进行锆石LA-ICP-MS U-Pb年代学、岩石地球化学和全岩Sr-Nd-Pb同位素以及锆石Hf同位素研究,探讨它们的源区性质和岩石成因及其形成的构造背景,进而为揭示鲁西地区早白垩世岩石圈地幔的属性和深部作用过程提供依据。

1 地质背景与样品描述

华北克拉通由东部陆块、西部陆块和中部造山带组成(Zhao et al., 2001)。鲁西地区位于华北克拉通东部陆块的东部,郯庐断裂带以西,东南与苏鲁造山带相邻(图 1a)。该区地层发育良好,属典型华北地层区,基底主要由新太古代泰山群组成,其上依次为元古界济宁群和土门群,古生界寒武系、中-下奥陶统、石炭系-二叠系(山东省地质矿产局, 2003)。中生界主要为断陷盆地沉积,发育三叠系刘家沟组、侏罗系坊子组和三台组以及白垩系汶南组,新生界主要为河湖相沉积(Yang et al., 2013)。鲁西地区除广泛分布有前寒武纪岩浆岩外,中生代侵入岩和火山岩也较为发育(Yang et al., 2012)。本文研究的大昆仑辉长岩以及侵入辉长岩中的辉绿岩是其典型代表之一。

图 1 华北克拉通构造简图(a)和鲁西大昆仑辉长岩地质图及采样位置(b) Fig. 1 Tectonic sketch map of North China Craton (a) and geological map showing the locations of the samples from the Dakunlun gabbros in western Shandong (b)

大昆仑辉长岩位于淄博市大昆仑镇宋家坊村附近(图 1b),出露面积约3km2,岩体侵入三叠系刘家沟组和侏罗系坊子组及三台组砂岩中,与围岩接触处多发育冷凝边(图 2a),辉长岩和侏罗系砂岩中发育较多北西向展布的晚期辉绿岩脉。大昆仑辉长岩(样号ZB-K),风化面呈土黄色,新鲜面呈灰黑色,辉绿辉长结构,块状构造。主要矿物为斜长石(~60%)、普通辉石(~38%)和少量的普通角闪石(~2%),粒度介于0.5~2mm之间(图 2b),斜长石牌号为55~58,属于拉长石,副矿物主要为磁铁矿。斜长石发育聚片双晶,少量绢云母化,普通辉石具二级中部干涉色,普通角闪石发生少量的绿泥石化,岩石中不含橄榄石。大昆仑辉绿岩(样号SK1)侵入大昆仑辉长岩中,新鲜面呈灰黑色,辉绿结构,块状构造。主要由斜长石(~60%)、普通辉石(~38%)和少量的石英(~2%)组成,粒度为0.5~1.5mm(图 2c)。斜长石呈自形长板状,斜长石牌号为53~55,属于拉长石,粒间被他形普通辉石充填构成辉绿结构,部分斜长石发育绢云母化,普通辉石存在少量的绿泥石化,副矿物主要为磁铁矿。

图 2 大昆仑辉长岩野外产状(a)和辉长岩(b)、辉绿岩(c)的矿物组成 Aug-普通辉石;Pl-斜长石;Hb-普通角闪石;Q-石英 Fig. 2 Field photo (a) for the Dakunlun gabbro and micrographs for the Dakunlun gabbro (b) and diabase (c)
2 分析方法

锆石分选在河北省廊坊市区域地质矿产调查研究所实验室完成。野外采集样品10kg左右,首先将岩石进行粉碎,然后采用淘选与电磁选方法进行锆石分选。在双目镜下挑选出没有包体和裂隙,透明度与晶形较好的锆石颗粒,将其用双面胶粘在载玻片上并套上靶环。注入环氧树脂后,冷却和打磨并露出锆石表面,在此基础上对锆石进行透射光、反射光和阴极发光(CL)图像的采集。锆石U-Pb定年在中国地质大学(武汉)地质过程与矿产资源国家重点实验室利用GeoLas 2005型激光剥蚀系统和Agilent 7500a型电感耦合等离子体质谱仪(ICP-MS)分析完成。激光剥蚀过程中采用氦气作为载气,氩气作为补偿气以调节灵敏度。U-Pb定年中采用国际锆石标样91500作为外标进行同位素分馏校正,锆石微量元素采用NIST 610作为外标、Si作为内标的方法进行定量计算。样品的U-Th-Pb同位素比值和年龄以及元素含量处理均采用ICPMSDataCal软件完成。详细的实验条件和数据处理方法见Liu et al. (2010)

主量元素、微量元素以及Sr-Nd-Pb同位素分析均在中国地质大学(武汉)地质过程与矿产资源国家重点实验室完成。主量元素采用X射线荧光光谱法(XRF)分析,在XRF上运用外标法测定样品中氧化物的含量,分析过程中运用BCR-2标样进行校正,分析精确度和准确度优于5%。微量元素采用Agilent 7500a型ICP-MS完成,分析精确度和准确度优于10%。Sr-Nd同位素分析采用HF+HNO3混合酸液解,用阳离子交换技术进行分离,分析精度优于0.002%。国际标样La Jolla和BCR-2的143Nd/144Nd测定平均值分别为0.511859±6(2σ, n=20)和0.512637±3(2σ, n=12),NBS-987给出的87Sr/86Sr值为0.710245±15(2σ, n=9)。Pb同位素采用AG50W-X12,200-400mesh离子交换树脂在超净实验室内进行化学分离,在测定样品之前,用NBS 981标准校定仪器,测试过程中采用NBS 997 Tl作为内标进行质量分馏校正,随时测定NBS981来监测仪器。分析过程中NBS 981标准测定值的平均值为:206Pb/204Pb=16.940、207Pb/204Pb=15.486和208Pb/204Pb=36.722。锆石Hf同位素在中国科学院地质与地球物理研究所配备193nm激光取样系统的Neptune MC-ICP-MS上进行。测定时锆石标样91500作为外标,所用的激光脉冲速率为6~8Hz,能量为100mJ。仪器的运行条件和分析流程见Wu et al. (2006)。对176Yb的干扰校正采用TIMS测定值176Yb/172Yb=0.5886,而对每个分析点的βYbβHf则利用对该分析点实测得出的平均值进行校正,锆石标样91500的176Hf/177Hf=0.2822934±24(2σ, n=35)。

3 分析结果 3.1 锆石U-Pb年代学

鲁西地区大昆仑辉长岩和辉绿岩中典型锆石CL图像如图 3所示。辉长岩中锆石主要分为两种类型(图 3a),一种为自形-半自形晶,数量少且粒径小,大小为30~50μm,不发育核边结构,CL图像呈均匀的条痕状吸收,与基性岩中锆石的特征相一致,暗示其为岩浆锆石(Koschek, 1993),这与它们高的Th/U比值(0.73~1.64) (表 1)吻合;另一种锆石呈浑圆状,粒径为80~100μm,发育核边结构,具有岩浆生长环带,结合相对高的Th/U比值(0.11~1.99)(表 1),表明其为具岩浆成因的捕获锆石。辉绿岩中锆石较单一,均为自形-半自形晶(图 3b),内部结构均匀,发育条痕状吸收,粒径50~100μm,其Th/U比值较高为0.83~3.12,暗示它们为岩浆锆石。

图 3 大昆仑辉长岩(a)和辉绿岩(b)中典型锆石阴极发光图像 Fig. 3 Cathodoluminescence images of selected zircons from the Dakunlun gabbro (a) and diabase (b)

表 1 大昆仑辉长岩和辉绿岩锆石LA-ICP-MS U-Pb定年结果 Table 1 LA-ICP-MS zircon U-Pb dating data from the Dakunlun gabbros and diabases

大昆仑辉长岩和辉绿岩锆石LA-ICP-MS U-Pb定年结果见表 1图 4。对辉长岩中锆石U-Pb同位素分析结果,较老的锆石(>1000Ma)采用207Pb/206Pb年龄,而对于年轻锆石(<1000Ma)采用206Pb/238U年龄,结果显示,岩浆锆石的年龄介于129~137Ma之间,206Pb/238U加权平均年龄为132±3Ma(MSWD=0.6, n=5),该年龄代表了大昆仑辉长岩的形成时代(图 4a)。捕获锆石的年龄主要介于2669~1432 Ma之间,少量锆石发生Pb丢失,但它们均沿不一致线分布,下交点均位于岩体形成时代附近,说明最晚期岩浆事件可能是造成Pb丢失的主要原因,其中谐和年龄可分为4个主要年龄段:2469~2386Ma、2263~2232Ma、2147~1997Ma和1908~1752Ma,其加权平均年龄分别为2411Ma(n=7)、2252Ma(n=4)、2039Ma(n=6)和1844Ma(n=10)。辉绿岩中岩浆锆石的年龄比较均一,不存在捕获锆石,其年龄介于109~114Ma之间,加权平均年龄为112±1Ma(MSWD=1.7, n=18),代表了辉绿岩的形成时代(图 4b)。

图 4 大昆仑辉长岩(a)和辉绿岩(b)锆石LA-ICP-MS U-Pb年龄谐和图 Fig. 4 Concordia diagrams showing LA-ICP-MS zircon U-Pb dating for the Dakunlun gabbro (a) and diabase (b)
3.2 岩石地球化学 3.2.1 主量元素

大昆仑辉长岩3个样品和辉绿岩2个样品的主量元素和微量元素分析结果见表 2。辉长岩和辉绿岩具有类似的主量元素组成,整体具低硅和富镁特征,它们的SiO2介于49.08%~52.61%之间,TiO2变化于0.92%~1.65%之间,Al2O3为12.90%~14.83%,MgO为4.52%~9.67%,Mg#[Mg#=Mg/(Mg+Fe2+)×100]值为45~65,CaO为7.59%~9.27%,Na2O为2.58%~3.49%,K2O为0.85%~2.11%。在TAS图上(图 5),辉长岩样品投影在辉长岩和二长辉长岩区域,辉绿岩样品分布在辉长闪长岩及其与辉长岩分界线附近,除ZB-K1样品外,其它样品均位于亚碱性系列范围内。在主量元素组成上,大昆仑辉长岩和辉绿岩与济南辉长岩相类似(Yang et al., 2012),而与同时代的鲁西地区辉绿岩墙相比(Liu et al., 2016),具有相对低的碱含量。在Harker图解上(图 6),大昆仑辉长岩和辉绿岩的SiO2总体上与Fe2O3T、MgO、P2O5和CaO呈负相关性,而与Na2O呈正相关性,Al2O3含量的变化趋势不明显,其演化趋势整体与济南辉长岩相类似(Yang et al., 2012)。

表 2 大昆仑辉长岩和辉绿岩主量元素(wt%)和微量元素(×10-6)分析结果 Table 2 Major (wt%) and trace (×10-6) element data of the Dakunlun gabbros and diabases

图 5 大昆仑辉长岩和辉绿岩TAS图解 数据来源:铁铜沟辉石闪长岩据Yang et al. (2005);上峪辉长闪长岩据Yang et al. (2008);济南辉长岩据Yang et al. (2012);鲁西辉绿岩墙据Liu et al. (2016)图 6-图 9图 12 Fig. 5 Total alkali vs. SiO2 variation diagram for Dakunlun gabbros and diabases

图 6 大昆仑辉长岩和辉绿岩Harker图解 Fig. 6 Harker diagrams for Dakunlun gabbros and diabases

图 7 大昆仑辉长岩和辉绿岩球粒陨石标准化稀土元素配分型式(a,标准化值据Boynton, 1984)和原始地幔标准化微量元素蛛网图(b, 标准化值据Sun and McDonough, 1989) 数据来源:方城-费县玄武岩据Zhang et al. (2002)Yang et al. (2012), 图 8 Fig. 7 Chondrite-normalized REE patterns (a, normalization values after Boynton, 1984) and primitive mantle-normalized trace element spider diagrams (b, normalization values after Sun and McDonough, 1989) for the Dakunlun gabbros and diabases

图 8 大昆仑辉长岩(87Sr/86Sr)i-εNd(t)变异图 华北克拉通下地壳和大陆上地壳据Jahn et al. (1999);金岭黑云母闪长岩据Yang et al. (2005)图 9 Fig. 8 Initial (87Sr/86Sr)i vs. εNd(t) diagram for the Dakunlun gabbros

图 9 大昆仑辉长岩Pb同位素组成(据Yang et al., 2012修改) 数据来源:华北克拉通中生代镁铁质岩据Li and Yang (2003);扬子克拉通中生代铁镁质据Li and Yang (2003)闫峻等(2003);北半球参考线据Hart (1984);地球等时线据Zindler and Hart (1986);徐淮高镁闪长岩据杨德彬等(2008) Fig. 9 Pb isotopic compositions of the Dakunlun gabbros (modified after Yang et al., 2012)

图 12 大昆仑辉长岩和辉绿岩的Ba/Rb-Rb/Sr (a)和K/Yb-Dy/Yb (b)图解(据Ma et al., 2014修改) Fig. 12 Ba/Rb vs. Rb/Sr (a) and K/Yb vs. Dy/Yb (b) for Dakunlun gabbros and diabases (modified after Ma et al., 2014)
3.2.2 微量元素

大昆仑辉长岩和辉绿岩具有类似的球粒陨石标准化稀土元素配分型式和原始地幔标准化微量元素蛛网图(图 7)。它们整体具轻稀土元素(LREE)中等程度富集、重稀土元素(HREE)亏损的特征(图 7a),轻、重稀土元素分异明显[LREE/HREE=7.52~8.72,(La/Yb)N=8.98~10.2],Eu负异常不明显,δEu值介于0.83~0.93之间。类似于同时代的鲁西地区辉绿岩墙的REE丰度(Liu et al., 2016),而与鲁西地区方城-费县玄武岩相比(Zhang et al., 2002; Yang et al., 2012),大昆仑辉长岩和辉绿岩具有相对偏低的LREE丰度,但它们整体高于济南辉长岩的REE丰度(Yang et al., 2012)。

原始地幔标准化微量元素蛛网图显示,大昆仑辉长岩与辉绿岩富集大离子亲石元素(LILE,如Cs、Rb、Ba),明显亏损高场强元素(HFSE,如Nb、Ta、Zr、Hf),具有明显的Ti负异常和Pb正异常(图 7b)。整体与鲁西地区同时代的济南辉长岩特征相类似(Yang et al., 2012),而与方城和费县玄武岩以及鲁西辉绿岩墙相比(Zhang et al., 2002; Yang et al., 2012; Liu et al., 2016),大昆仑辉长岩与辉绿岩具有相对低的Th、U和Pb丰度,它们的分布模式与大陆弧钙碱性玄武岩相类似,而不同于碱性洋岛玄武岩和板内玄武岩的特征。此外,大昆仑辉长岩与辉绿岩的微量元素组成与MORB相比(Pearce, 1982),具有相对高的Rb、Ba、Th含量(19.6×10-6~53.7×10-6、519×10-6~1129×10-6、1.12×10-6~2.30×10-6)和相对低的Nb和Ta含量(4.29×10-6~7.74×10-6和0.19×10-6~0.35×10-6)。

3.3 全岩Sr-Nd-Pb同位素

大昆仑辉长岩的Sr-Nd-Pb同位素组成列于表 3。从表 3图 8中可以看出,大昆仑辉长岩3个样品整体具有相对低的初始87Sr/86Sr比值和εNd(t)值,它们分别变化于0.7046~0.7060之间和-16.71~-15.54之间, 其中样品ZB-K2具有相对较低的初始87Sr/86Sr比值,而样品ZB-K1和ZB-K3的Sr-Nd同位素组成相接近。与亏损地幔相比(Sun and McDonough, 1989),大昆仑辉长岩表现出相对较高的初始87Sr/86Sr比值和较低的εNd(t)值。与同时代的鲁西地区济南辉长岩和金岭黑云母闪长岩及铁铜沟辉石闪长岩相比(Yang et al., 2008, 2012),大昆仑辉长岩的εNd(t)值较低,而与方城和费县玄武岩、鲁西辉绿岩墙以及上峪辉长闪长岩相比(Zhang et al., 2002; Yang et al., 2008, 2012; Liu et al., 2016),大昆仑辉长岩的初始87Sr/86Sr比值明显偏低,而εNd(t)值相近(图 8)。

表 3 大昆仑辉长岩Sr-Nd-Pb同位素组成 Table 3 Sr-Nd-Pb isotopic composition of the Dakunlun gabbros

大昆仑辉长岩的(206Pb/204Pb)t、(207Pb/204Pb)t和(208Pb/204Pb)t比值分别变化于16.596~16.776、15.255~15.290和36.328~36.648之间。在Pb同位素组成图解上(图 9),大昆仑辉长岩整体具有较低的Pb同位素组成,均位于北半球参考线之上,与鲁西地区金岭黑云母闪长岩和济南辉长岩(Yang et al., 2012)相类似,但低于铁铜沟辉石闪长岩和上峪辉长闪长岩(Yang et al., 2008, 2012)。鲁西地区多数样品的初始Pb同位素组成投影于华北克拉通中生代镁铁质岩浆岩Pb同位素组成范围内(图 9),而明显低于徐淮地区高镁闪长岩的Pb同位素组成(杨德彬等, 2008),后者均投影在扬子克拉通中生代镁铁质岩浆岩Pb同位素组成范围内(图 9)。

3.4 锆石原位Hf同位素

在锆石U-Pb年代学分析的基础上,对大昆仑辉长岩和辉绿岩进行了锆石原位Hf同位素分析,因辉长岩中代表岩浆形成时代的岩浆锆石粒度较小,在锆石定年过程中已全部消耗掉,因此仅对辉长岩中少量捕获锆石进行了Hf同位素分析,结果见表 4图 10。从表 4中可以看出,辉长岩中1752Ma、1755Ma锆石的176Hf/177Hf值是0.281324、0.281341,其εHf(t)值分别为-13.1、-14.1,而2147Ma、2248Ma锆石的176Hf/177Hf值是0.281261和0.281271,其εHf(t)值分别为-6.9、-3.6。辉绿岩中岩浆锆石的176Hf/177Hf值介于0.282369~0.282644之间,其εHf(t)值介于-12.0~-2.5之间。在εHf(t)对U-Pb年龄图中(图 10),辉绿岩锆石的Hf同位素组成均投影于球粒陨石之下,而明显低于亏损地幔组成。

表 4 大昆仑辉长岩和辉绿岩锆石Hf同位素分析结果 Table 4 Zircon Hf isotopic data for the Dakunlun gabbros and diabases

图 10 大昆仑辉长岩和辉绿岩锆石Hf同位素组成 Fig. 10 Zircon Hf isotopic compositions of the Dakunlun gabbros and diabases
4 讨论 4.1 大昆仑辉长岩和辉绿岩的形成时代

山东省区域地质志(山东省地质矿产局, 2003)基于岩体与围岩的接触关系,即大昆仑辉长岩侵入三叠系刘家沟组和侏罗系三台组(图 2),将该辉长岩的时代定为早白垩世,但长期以来缺乏精确同位素年代学的制约。大昆仑辉长岩中少量锆石为自形-半自形晶,呈条痕状均匀吸收,与典型基性岩中锆石的特征相一致,结合较高的Th/U比值(0.73~1.64),表明它们为岩浆锆石,其206Pb/238U年龄介于129~137Ma之间,加权平均年龄为132±3Ma,这表明辉长岩的结晶年龄为早白垩世,这与辉长岩侵入侏罗系三台组的地质事实相一致(山东省地质矿产局, 2003)。此外,对侵入辉长岩中辉绿岩脉的锆石U-Pb定年结果表明,辉绿岩中锆石均呈条痕状吸收,暗示它们为岩浆成因。岩浆锆石的定年结果介于109~114Ma之间,其206Pb/238U加权平均年龄为112±1Ma,即形成于早白垩世晚期,稍晚于辉长岩的结晶时代,并与野外的侵入接触关系相吻合(图 1)。上述定年结果表明,大昆仑辉长岩和辉绿岩均形成于早白垩世,这与济南辉长岩的时代相一致(127~131Ma; 杨承海等, 2005),同时与鲁西地区辉绿岩墙的锆石U-Pb定年结果相类似(122~124Ma; Liu et al., 2016),也与鲁西地区金岭黑云母闪长岩和铁铜沟辉石闪长岩的锆石U-Pb定年结果(131~135Ma; Yang et al., 2005)和上峪辉长闪长岩锆石U-Pb定年结果相类似(129~134Ma; Yang et al., 2008),并与华北克拉通东部(蚌埠、徐淮、胶东和辽东地区)早白垩世侵入岩(辉长岩和花岗质岩石)中锆石LA-ICP-MS和SHRIMP U-Pb定年结果相一致(Xu et al., 2004; Yang et al., 2004, 2010; Wu et al., 2005)。鲁西地区大昆仑辉长岩和辉绿岩以及鲁西地区乃至华北克拉通东部辉长质和花岗质岩石的形成时代主要为早白垩世,这与华北克拉通最强烈的一次岩浆活动相吻合(Xu et al., 2004; Yang et al., 2004, 2005, 2008, 2010, 2017; Wu et al., 2005; Li et al., 2017)。

4.2 源区性质和岩石成因

鲁西早白垩世大昆仑辉长岩和辉绿岩具有低硅(SiO2=49.08%~52.61%)、高镁(MgO=4.52%~9.67%)和高的Mg#值(45~65)以及富钠(Na2O=2.58%~3.49%)的特征,结合相对高的Sc、Co、Ni含量(25.3×10-6~32.5×10-6、30.0×10-6~47.8×10-6和14.5×10-6~91.6×10-6),表明它们具有地幔的源区性质(Rapp and Watson, 1995)。大昆仑辉长岩相对高的87Sr/86Sr初始比值(0.7046~0.7060)、低的εNd(t)值(-16.71~-15.54),暗示大昆仑辉长岩应起源于富集岩石圈地幔的部分熔融,而非软流圈地幔,因为后者通常具有亏损的Nd同位素组成(εNd(t)值>0)和低的87Sr/86Sr初始比值。此外,大昆仑辉长岩富集LREE和LILE,与类似MORB具有亏损LREE、LILE特征的软流圈地幔明显不同也证明了这一点(Sun and McDonough, 1989)。锆石Hf同位素组成很少受到后期地质事件的影响,是探讨岩浆源区属性的重要工具(Griffin et al., 2002)。早白垩世大昆仑辉绿岩中岩浆锆石的εHf(t)值介于-12.0~-2.5之间,这与大昆仑辉长岩全岩的εNd(t)值(-16.71~-15.54)相类似,总体上也显示出富集的锆石Hf同位素组成。结合大昆仑辉长岩和辉绿岩具有相同的产出位置和相类似的主量和微量元素组成,暗示它们应具有相类似的源区性质,即起源于富集岩石圈地幔的部分熔融。大昆仑辉长岩和辉绿岩明显亏损高场强元素Nb、Ta、Zr和Hf以及具有Pb的正异常,暗示它们的地幔源区存在陆壳物质的改造,这与辉长岩中古元古代(2421~1837Ma)捕获锆石的存在相一致(图 4a)。

那么,改造富集岩石圈地幔源区属性的壳源物质是来自于大陆上地壳还是下部大陆地壳。首先,大昆仑辉长岩的87Sr/86Sr初始比值介于0.7046~0.7060之间,整体表现出EMⅠ的演化特征,与受下地壳物质改造的济南辉长岩的Sr-Nd同位素组成相接近(Yang et al., 2012),而与受上地壳物质影响的方城和费县玄武岩的同位素组成明显不同(图 8);其次,大昆仑辉长岩全岩Pb同位素初始比值[(206Pb/204Pb)t=16.596~16.776、(207Pb/204Pb)t=15.255~15.290]与全球麻粒岩相下地壳岩石Pb同位素组成相接近(西格陵兰和北美拉布拉多地区206Pb/204Pb=12~18、207Pb/204Pb=13.5~15.0;苏格兰Lewisian麻粒岩206Pb/204Pb=13.3~16.8、207Pb/204Pb=14.5~15.3;南印度麻粒岩206Pb/204Pb=14.2~16.0、207Pb/204Pb=15.2~15.6),并且介于华北麻粒岩相下地壳物质Pb同位素组成范围内(206Pb/204Pb=14.3~17.8、207Pb/204Pb=14.8~15.5)(涂湘林等, 1993);最后,大昆仑辉长岩中存在大量的古元古代捕获锆石(图 3),也暗示了古老下部大陆地壳物质的参与。综上所述,我们认为改造大昆仑辉长岩和辉绿岩地幔源区的陆壳物质主要来自于下部大陆地壳。然而,对于地幔源区部分熔融形成的岩石,被下地壳物质改造的机制主要有两种,即拆沉陆壳物质对岩浆源区的直接改造和岩浆上升侵位过程中被陆壳物质的同化混染。这两种岩浆演化过程的本质区别在于,前者对岩浆源区的改造是地幔物质与地壳物质发生熔融并最终形成均一的元素-同位素地球化学体系,而后者是幔源岩浆上升侵位到地壳深度过程中同化混染了地壳物质,并导致岩浆的元素-同位素组成发生明显的变化(Yang et al., 2008)。在εNd(t)值-SiO2和(87Sr/86Sr)i-SiO2图中(图 11),大昆仑辉长岩样品没有显示明显的线性关系,而与分离结晶的演化趋势相一致,表明岩浆演化过程中同化混染作用不明显,因此,下地壳物质对幔源源区的改造应主要发生在岩浆源区,而辉长岩中古元古代锆石被捕获的过程应主要发生在岩浆房内,而非岩浆上升过程中捕获浅部陆壳物质的。鲁西地区铁铜沟高镁闪长岩中地幔橄榄岩包体(Xu et al., 2008)和费县玄武岩中幔源捕掳体和捕虏晶的岩相学、矿物化学和矿物原位Sr-O同位素研究表明(裴福萍等, 2004; Xu et al., 2013),鲁西地区早白垩世岩石圈地幔受到了拆沉大陆下地壳物质的强烈改造。因此,改造鲁西地区大昆仑辉长岩和辉绿岩富集岩石圈地幔源区的壳源熔体主要来自于拆沉的古老下地壳物质的部分熔融。

图 11 大昆仑辉长岩的(87Sr/86Sr)i-SiO2 (a)和εNd(t)-SiO2 (b)变异图 Fig. 11 (87Sr/86Sr)i vs. SiO2 (a) and εNd(t) vs. SiO2 (b) for the Dakunlun gabbros

那么,改造富集岩石圈地幔源区属性的大陆下地壳物质是拆沉的华北古老下地壳还是俯冲断离的扬子陆壳物质,这可以通过Pb同位素的示踪得到回答(Yang et al., 2012)。中国东部几个重要构造块体之间Pb同位素组成具有明显的差异,Li and Yang (2003)根据中国东部中生代镁铁质岩浆岩的Pb同位素研究表明,华北克拉通的镁铁质岩浆岩具有较低的放射性Pb同位素初始比值,通常206Pb/204Pb<17.80、207Pb/204Pb<15.50和208Pb/204Pb<38.00,而扬子克拉通镁铁质岩浆岩则以高的放射成因Pb同位素初始比值为特征,其中206Pb/204Pb>17.80、207Pb/204Pb>15.50和208Pb/204Pb>38.00,并以此作为区别扬子克拉通基底和华北克拉通基底物质的重要标志(杨德彬等, 2008; Yang et al., 2010, 2012, 2016)。鲁西早白垩世大昆仑辉长岩具有低的全岩Pb同位素初始比值[(206Pb/204Pb)t=16.596~16.776、(207Pb/204Pb)t=15.255~15.290、(208Pb/204Pb)t=36.328~36.648],介于华北克拉通麻粒岩相带Pb同位素组成范围内(206Pb/204Pb=14.3~17.8、207Pb/204Pb=14.8~15.5、208Pb/204Pb=33.4~37.8,涂湘林等, 1993),它们均投影在华北克拉通中生代镁铁质岩的Pb同位素组成范围内(图 9),与济南辉长岩的Pb同位素组成相类似,但低于具有扬子和华北克拉通基底混合源区特征的上峪和铁铜沟辉长闪长质岩石的Pb同位素组成(Yang et al., 2012),而与徐淮地区具有高放射性Pb同位素组成的早白垩世高镁闪长岩明显不同,后者具有扬子克拉通中生代镁铁质岩浆的属性(杨德彬等, 2008)。因此,上述结果进一步表明,改造鲁西地区早白垩世大昆仑辉长岩和辉绿岩岩浆源区的壳源熔体主要来源于拆沉的华北克拉通下地壳物质的部分熔融,这与岩石中大量古元古代捕获锆石的存在相吻合。

实验岩石学研究证明含水矿物金云母和角闪岩只能稳定存在于岩石圈地幔(Olafsson and Eggler, 1983)。Rb和Ba元素在金云母中为不相容元素,而Rb、Sr和Ba在角闪石中为中等相容性(LaTourrette et al., 1995),因此,含金云母源区的熔体具有低的Ba和Ba/Rb比值(<20),而含角闪石源区的熔体具有较低的Rb/Sr(<0.1)和高的Ba/Rb比值(>20)(Furman and Graham, 1999)。大昆仑辉长岩和辉绿岩的Rb/Sr比值为0.03~0.07,Ba/Rb比值为18.6~32.5,在图 12a中可以看出,大昆仑辉长岩和辉绿岩样品均位于源区为含交代成因角闪石的岩石圈地幔区域,角闪石为难熔的古老岩石圈地幔提供了绝大多数不相容元素,这与鲁西地区上峪辉石闪长岩和济南辉长岩的源区性质相似(Yang et al., 2008, 2012)。此外,K/Yb和Dy/Yb比值可用于区分含金云母石榴石相或含角闪石尖晶石相二辉橄榄岩部分熔融的源区(Duggen et al., 2005)。大昆仑辉长岩和辉绿岩轻稀土元素中等程度富集[LREE/HREE=7.52~8.72; (La/Yb)N=8.98~10.2],表明它们经历了较高程度的部分熔融,大昆仑辉长岩和辉绿岩的K/Yb×1000=4.37~6.92,Dy/Yb=2.11~2.37,在图 12b中投影于源区深度在尖晶石和石榴石转换带附近,暗示它们的源区深度可能在60~80km左右(Lambert and Wyllie, 1968)。大昆仑辉长岩和辉绿岩中具有Ti的负异常,表明岩浆演化过程中存在钛铁矿的分离结晶,而弱的Eu负异常暗示,斜长石不是源区主要的残留相矿物。在Harker图解上(图 6),大昆仑辉长岩和辉绿岩的SiO2含量总体上与Fe2O3T、MgO、P2O5和CaO含量呈负相关性(图 6),也说明了存在钛铁矿和橄榄石以及辉石等矿物的分离结晶。

4.3 构造背景

大昆仑辉长岩和辉绿岩形成于早白垩世(112~132Ma),这与华北克拉通东部广泛发育的早白垩世岩浆活动相一致,暗示它们具有一个共同的形成环境。鲁西地区早白垩世岩浆活动广泛而强烈,形成了大量的侵入岩和火山岩,这包括济南辉长岩(杨承海等, 2005)、鲁西辉绿岩墙(Liu et al., 2016)和金岭黑云母闪长岩、铁铜沟辉石闪长岩(Yang et al., 2005)以及上峪辉长闪长岩(Yang et al., 2008)和早白垩世方城和费县玄武岩等(Zhang et al., 2002; 裴福萍等, 2004)。此外,形成于早白垩世的蚌埠隆起区花岗岩(Yang et al., 2010)、徐淮高镁闪长岩(Xu et al., 2004)、胶东花岗岩(Yang et al., 2004)以及辽东地区大面积的早白垩世花岗岩(Wu et al., 2005)均为中国东部早白垩世时期最重要的一次岩浆事件。早白垩世的岩浆事件在华北克拉通东部乃至中国东部主要表现为双峰式的火成岩组合(许文良等, 2004),表明它们形成于伸展的构造背景,这与华北克拉通东部变质核杂岩和大型伸展构造以及拉分盆地的形成时间相吻合(Wu et al., 2007),同时与华北克拉通东部岩石圈强烈减薄的峰期时间相对应(许文良等, 2004)。综上所述,鲁西地区大昆仑辉长岩和辉绿岩形成于华北克拉通东部早白垩世时期岩石圈强烈伸展的构造背景,而这一时期与古太平洋板块向欧亚大陆东缘之下的斜向俯冲时间相一致,包括大昆仑辉长岩和辉绿岩在内的鲁西地区大量早白垩世辉长质和花岗质岩石均应形成于这一共同的构造背景之下。

5 结论

(1) 鲁西大昆仑辉长岩和辉绿岩中岩浆锆石的206Pb/238U加权平均年龄分别为132Ma和112Ma,即它们均形成于早白垩世。

(2) 大昆仑辉长岩和辉绿岩整体具有低SiO2、高MgO、Na2O含量以及高Mg#值特征,富集轻稀土元素和大离子亲石元素、贫重稀土元素和高场强元素,具有相对高的初始87Sr/86Sr比值和低的εNd(t)值、εHf(t)值,暗示它们起源于含交代成因角闪石的富集岩石圈地幔源区的部分熔融,同时岩浆源区受到了拆沉华北克拉通古老下地壳物质的强烈改造。

(3) 大昆仑辉长岩和辉绿岩形成于华北克拉通东部早白垩世时期岩石圈强烈伸展的构造背景。

致谢      样品分析得到了中国地质大学(武汉)地质过程与矿产资源国家重点实验室和中国科学院地质与地球物理研究所岩石圈演化国家重点实验室的大力支持;两位匿名审稿专家和编辑提供了建设性修改意见;谨此一并表示感谢。

参考文献
Amelin Y, Lee DC and Halliday AN. 2000. Early-Middle Archaean crustal evolution deduced from Lu-Hf and U-Pb isotopic studies of single zircon grains. Geochimica et Cosmochimica Acta, 64(24): 4205-4225. DOI:10.1016/S0016-7037(00)00493-2
Bizzarro M, Simonetti A, Stevenson RK and David J. 2002. Hf isotope evidence for a hidden mantle reservoir. Geology, 30(9): 771-774. DOI:10.1130/0091-7613(2002)030<0771:HIEFAH>2.0.CO;2
Boynton WV. 1984. Geochemistry of the rare earth elements: Meteorite studies. In: Henderson P (ed.). Rare Earth Element Geochemistry. Amsterdam: Elsevier, 63-114
Bureau of Geology and Mineral Resources of Shandong Province. 2003. Regional Geology of Shandong Province. Ji'nan: Shandong Mapping Press: 1-350.
Duggen S, Hoernle K, Van Den Bogaard P and Garbe-Schönberg D. 2005. Post-collisional transition from subduction-to intraplate-type magmatism in the westernmost Mediterranean:Evidence for continental-edge delamination of subcontinental lithosphere. Journal of Petrology, 46(6): 1155-1201. DOI:10.1093/petrology/egi013
Furman T and Graham D. 1999. Erosion of lithospheric mantle beneath the East African Rift system:Geochemical evidence from the Kivu volcanic province. Developments in Geotectonics, 24: 237-262. DOI:10.1016/S0419-0254(99)80014-7
Griffin WL, Pearson NJ, Belousova E, Jackson SE, Van Achterbergh E, O'Reilly SY and Shee SR. 2000. The Hf isotope composition of cratonic mantle:LAM-MC-ICPMS analysis of zircon megacrysts in kimberlites. Geochimica et Cosmochimica Acta, 64(1): 133-147. DOI:10.1016/S0016-7037(99)00343-9
Griffin WL, Wang X, Jackson SE, Pearson NJ, O'Reilly SY, Xu XS and Zhou XM. 2002. Zircon chemistry and magma mixing, SE China:In-situ analysis of Hf isotopes, Tonglu and Pingtan igneous complexes. Lithos, 61(3-4): 237-269. DOI:10.1016/S0024-4937(02)00082-8
Hart SR. 1984. A large-scale isotope anomaly in the Southern Hemisphere mantle. Nature, 309: 753-757. DOI:10.1038/309753a0
Huang XL, Zhong JW and Xu YG. 2012. Two tales of the continental lithospheric mantle prior to the destruction of the North China Craton:Insights from Early Cretaceous mafic intrusions in western Shandong, East China. Geochimica et Cosmochimica Acta, 96: 193-214. DOI:10.1016/j.gca.2012.08.014
Jahn BM, Wu FY, Lo CH and Tsai CH. 1999. Crust-mantle interaction induced by deep subduction of the continental crust:Geochemical and Sr-Nd isotopic evidence from post-collisional mafic-ultramafic intrusions of the northern Dabie complex, central China. Chemical Geology, 157(1-2): 119-146. DOI:10.1016/S0009-2541(98)00197-1
Koschek G. 1993. Origin and significance of the SEM cathodoluminescence from zircon. Journal of Microscopy, 171(3): 223-232. DOI:10.1111/jmi.1993.171.issue-3
Lambert IB and Wyllie PJ. 1968. Stability of hornblende and a model for the low velocity zone. Nature, 219(5160): 1240-1241. DOI:10.1038/2191240a0
LaTourrette T, Hervig RL and Holloway JR. 1995. Trace element partitioning between amphibole, phlogopite, and basanite melt. Earth and Planetary Science Letters, 135(1-4): 13-30. DOI:10.1016/0012-821X(95)00146-4
Li FL and Li WS. 2016. Controlling factors for dawsonite diagenesis:A case study of the Binnan region in Dongying Sag, Bohai bay Basin, China. Australian Journal of Earth Sciences, 63(2): 217-233. DOI:10.1080/08120099.2016.1173096
Li FL, Li WS, Yu ZC, Liu N, Yang HD and Liu L. 2017. Dawsonite occurrences related to the age and origin of CO2 influx in sandstone reservoirs:A case study in the Songliao Basin, NE China. Geochemistry, Geophysics, Geosystems, 18(1): 346-368. DOI:10.1002/2016GC006555
Li SG and Yang W. 2003. Decoupling of surface and subsurface sutures in the Dabie orogen and a continent-collisional lithospheric-wedging model:Sr-Nd-Pb isotopic evidences of Mesozoic igneous rocks in eastern China. Chinese Science Bulletin, 48(8): 831-838. DOI:10.1007/BF03187062
Liu S, Feng CX, Zhai MG, Hu RZ, Lai SC, Chen JJ and Yan J. 2016. Zircon U-Pb age, geochemical, and Sr-Nd-Hf isotopic constraints on the origin of Early Cretaceous mafic dykes from western Shandong Province, eastern North China Craton, China. Acta Petrologica Sinica, 32(3): 629-645.
Liu YS, Gao S, Hu ZC, Gao CG, Zong KQ and Wang DB. 2010. Continental and oceanic crust recycling-induced melt-peridotite interactions in the trans-north China Orogen:U-Pb dating, Hf isotopes and trace elements in zircons from mantle xenoliths. Journal of Petrology, 51(1-2): 537-571. DOI:10.1093/petrology/egp082
Ma L, Jiang SY, Hofmann AW, Dai BZ, Hou ML, Zhao KD, Chen LH, Li JW and Jiang YH. 2014. Lithospheric and asthenospheric sources of lamprophyres in the Jiaodong Peninsula:A consequence of rapid lithospheric thinning beneath the North China Craton?. Geochimica et Cosmochimica Acta, 124: 250-271. DOI:10.1016/j.gca.2013.09.035
Olafsson M and Eggler DH. 1983. Phase relations of amphibole, amphibole-carbonate, and phlogopite-carbonate peridotite:Petrologic constraints on the asthenosphere. Earth and Planetary Science Letters, 64(2): 305-315. DOI:10.1016/0012-821X(83)90212-1
Pearce JA. 1982. Trace element characteristics of lavas from destructive plate boundaries. In: Thorpe RS (ed.). Andesites: Orogenic Andesites and Related Rocks. New York: John Wiley and Sons, 525-548
Pei FP, Xu WL, Wang QH, Wang DY and Lin JQ. 2004. Mesozoic basalt and mineral chemistry of the mantle-derived Xenocrysts in Feixian, western Shandong, China:Constraints on nature of Mesozoic lithospheric mantle. Geological Journal of China Universities, 10(1): 88-97.
Rapp RP and Watson EB. 1995. Dehydration melting of metabasalt at 8~32kbar:Implications for continental growth and crust-mantle recycling. Journal of Petrology, 36(4): 891-931. DOI:10.1093/petrology/36.4.891
Söderlund U, Patchett PJ, Vervoort JD and Isachsen CE. 2004. The 176Lu decay constant determined by Lu-Hf and U-Pb isotope systematics of Precambrian mafic intrusions. Earth and Planetary Science Letters, 219(3-4): 311-324. DOI:10.1016/S0012-821X(04)00012-3
Sun SS and McDonough WF. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. In: Saunders AD and Norry MJ (eds.). Magmatism in the Ocean Basins. Geological Society, London, Special Publications, 42(1): 313-345
Tu XL, Zhu BQ, Fan SK and Hu AQ. 1993. Pb isotopic compositions of granulite in North China and global comparison. Science in China (Series B), 23(5): 537-544.
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(1-2): 103-119. DOI:10.1016/j.epsl.2005.02.019
Wu FY, Yang YH, Xie LW, Yang JH and Xu P. 2006. Hf isotopic compositions of the standard zircons and baddeleyites used in U-Pb geochronology. Chemical Geology, 234(1-2): 105-126. DOI:10.1016/j.chemgeo.2006.05.003
Wu FY, Han RH, Yang JH, Wilde SA, Zhai MG and Park SC. 2007. Initial constraints on the timing of granitic magmatism in North Korea using U-Pb zircon geochronology. Chemical Geology, 238(3-4): 232-248. DOI:10.1016/j.chemgeo.2006.11.012
Xu WL, Wang QH, Liu XC, Wang DY and Guo JH. 2004. Chronology and sources of mesozoic intrusive complexes in the Xuzhou-Huainan region, central China:Constraints from SHRIMP zircon U-Pb dating. Acta Geologica Sinica, 78(1): 96-106.
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): 309-317.
Xu WL, Hergt JA, Gao S, Pei FP, Wang W and Yang DB. 2008. Interaction of adakitic melt-peridotite:Implications for the high-Mg# signature of Mesozoic adakitic rocks in the eastern North China Craton. Earth and Planetary Science Letters, 265(1-2): 123-137. DOI:10.1016/j.epsl.2007.09.041
Xu WL, Zhou QJ, Pei FP, Yang DB, Gao S, Li QL and Yang YH. 2013. Destruction of the North China Craton:Delamination or thermal/chemical erosion? Mineral chemistry and oxygen isotope insights from websterite xenoliths. Gondwana Research, 23(1): 119-129. DOI:10.1016/j.gr.2012.02.008
Yan J, Chen JF, Yu G, Qian H and Zhou TX. 2003. Pb isotopic characteristics of Late Mesozoic mafic rocks from the Lower Yangtze Region:Evidence for enriched mantle. Geological Journal of China Universities, 9(2): 195-206.
Yang CH, Xu WL, Yang DB, Liu CC, Liu XM and Hu ZC. 2005. Petrogenesis of Mesozoic high-mg diorites in western Shandong:Evidence from chronology and petro-geochemistry. Journal of China University of Geosciences, 16(4): 297-308.
Yang CH, Xu WL, Yang DB, Liu CC, Liu XM and Hu ZC. 2005. Chronology of the Jinan gabbro in western Shandong:Evidence from LA-ICP-MS zircon U-Pb dating. Acta Geoscientica Sinica, 26(4): 321-325.
Yang CH, Xu WL, Yang DB, Wang W, Wang WD and Liu JM. 2008. Petrogenesis of Shangyu gabbro-diorites in western Shandong:Geochronological and geochemical evidence. Science in China (Series D), 51(4): 481-492. DOI:10.1007/s11430-008-0029-0
Yang DB, Xu WL, Pei FP, Wang QH and Gao S. 2008. Chronology and Pb isotope compositions of Early Cretaceous adakitic rocks in Xuzhou-Huaibei area, central China:Constraints on magma sources and tectonic evolution in the eastern North China Craton. Acta Petrologica Sinica, 24(8): 1745-1758.
Yang DB, Xu WL, Wang QH and Pei FP. 2010. Chronology and geochemistry of Mesozoic granitoids in the Bengbu area, central China:Constraints on the tectonic evolution of the eastern North China Craton. Lithos, 114(1-2): 200-216. DOI:10.1016/j.lithos.2009.08.009
Yang DB, Xu WL, Pei FP, Yang CH and Wang QH. 2012. Spatial extent of the influence of the deeply subducted South China Block on the southeastern North China Block:Constraints from Sr-Nd-Pb isotopes in Mesozoic mafic igneous rocks. Lithos, 136-139: 246-260. DOI:10.1016/j.lithos.2011.06.004
Yang DB, Xu WL, Xu YG, Pei FP and Wang F. 2013. Provenance of sediments from Mesozoic basins in western Shandong:Implications for the evolution of the eastern North China Block. Journal of Asian Earth Sciences, 76: 12-29. DOI:10.1016/j.jseaes.2013.07.027
Yang DB, Xu WL, Zhao GC, Huo TF, Shi JP and Yang HT. 2016. Tectonic implications of Early Cretaceous low-Mg adakitic rocks generated by partial melting of thickened lower continental crust at the southern margin of the central North China Craton. Gondwana Research, 38: 220-237. DOI:10.1016/j.gr.2015.11.013
Yang DB, Yang HT, Shi JP, Xu WL and Wang F. 2017. Sedimentary response to the paleogeographic and tectonic evolution of the southern North China Craton during the Late Paleozoic and Mesozoic. Gondwana Research, 49: 278-295. DOI:10.1016/j.gr.2017.06.009
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(3-4): 145-160. DOI:10.1016/j.lithos.2003.12.003
Zhang HF, Sun M, 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(2): 241-254. DOI:10.1007/s00410-002-0395-0
Zhang HF, Sun M, Zhou MF, Fan WM, Zhou XH and Zhai MG. 2004. Highly heterogeneous Late Mesozoic lithospheric mantle beneath the North China Craton:Evidence from Sr-Nd-Pb isotopic systematics of mafic igneous rocks. Geological Magazine, 141(1): 55-62. DOI:10.1017/S0016756803008331
Zhao GC, Wilde SA, Cawood PA and Sun M. 2001. Archean blocks and their boundaries in the North China Craton:Lithological, geochemical, structural and P-T path constraints and tectonic evolution. Precambrian Research, 107(1-2): 45-73. DOI:10.1016/S0301-9268(00)00154-6
Zhong JW and Huang XL. 2012. Spatial variation of zircon Hf isotopes for the Early Cretaceous mafic intrusions in western Shandong and its genesis. Geotectonica et Metallogenia, 36(4): 572-580.
Zindler A and Hart S. 1986. Chemical geodynamics. Annual Review of Earth and Planetary Sciences, 14: 493-571. DOI:10.1146/annurev.ea.14.050186.002425
刘燊, 冯彩霞, 翟明国, 胡瑞忠, 赖绍聪, 陈俊瑾, 颜军. 2016. 华北克拉通东部鲁西早白垩世基性岩墙成因:来自锆石年龄、地球化学和Sr-Nd-Hf同位素的证据. 岩石学报, 32(3): 629-645.
裴福萍, 许文良, 王清海, 王冬艳, 林景仟. 2004. 鲁西费县中生代玄武岩及幔源捕掳晶的矿物化学:对岩石圈地幔性质的制约. 高校地质学报, 10(1): 88-97. DOI:10.3969/j.issn.1006-7493.2004.01.008
山东省地质矿产局. 2003. 山东省区域地质志. 济南: 山东省地图出版社: 1-350.
涂湘林, 朱炳泉, 范嗣昆, 胡霭琴. 1993. 华北麻粒岩相带铅同位素组成特征及其与全球对比. 中国科学(B辑), 23(5): 537-544.
许文良, 王清海, 王冬艳, 裴福萍, 高山. 2004. 华北克拉通东部中生代岩石圈减薄的过程与机制:中生代火成岩和深源捕虏体证据. 地学前缘, 11(3): 309-317. DOI:10.3321/j.issn:1005-2321.2004.03.029
闫峻, 陈江峰, 喻钢, 钱卉, 周泰禧. 2003. 长江中下游晚中生代中基性岩的铅同位素特征:富集地幔的证据. 高校地质学报, 9(2): 195-206. DOI:10.3969/j.issn.1006-7493.2003.02.005
杨承海, 许文良, 杨德彬, 刘长春, 柳小明, 胡兆初. 2005. 鲁西济南辉长岩的形成时代:锆石LA-ICP-MS U-Pb定年证据. 地球学报, 26(4): 321-325. DOI:10.3321/j.issn:1006-3021.2005.04.005
杨德彬, 许文良, 裴福萍, 王清海, 高山. 2008. 徐淮地区早白垩世adakitic岩石的年代学和Pb同位素组成:对岩浆源区与华北克拉通东部构造演化的制约. 岩石学报, 24(8): 1745-1758.
钟军伟, 黄小龙. 2012. 鲁西早白垩世基性侵入岩的锆石Hf同位素组成变化及其成因. 大地构造与成矿学, 36(4): 572-580. DOI:10.3969/j.issn.1001-1552.2012.04.011