华北克拉通辽北清原地体新太古代基性麻粒岩变质作用演化

崔润泽,魏春景

北京大学地球与空间科学学院,北京 100871

0 引言

太古宙克拉通是大陆地壳中的重要构造单元,许多太古宙克拉通具有 “穹隆-龙骨 (dome-andkeel) ” 构造,表现为大小不一的表壳岩呈块状/皮筏状产于TTG (tonalite-trondhjemite-granodiorite)穹隆内,或呈带状处在穹隆之间 (如: Collins et al.,1998;Hickman,2004;Lin and Beakhouse,2013;Anhaeusser,2014)。这些表壳岩经历了绿片岩相-麻粒岩相变质作用 (Condie,1981)。Brown(2007) 和Brown and Johnson (2018) 对17 个太古宙高级地体的总结表明,麻粒岩相表壳岩的P-T条件为0.65～1.40 GPa/800～1100 ℃,地热梯度很高,大于775 ℃/GPa。这些高温—超高温麻粒岩大多显示逆时针的P-T轨迹,但它们的构造体制仍存在争议。例如,Mezger et al.(1990) 认为苏必利尔克拉通 (Superior Craton) 的～2.7 Ga 皮奎托内 (Pikwitonei) 麻粒岩形成于大陆岩浆弧环境;Jayananda et al.(2012) 认为东达瓦尔克拉通(Eastern Dharwar Craton) 中部的～2.7 Ga 泥质麻粒岩可能与洋脊俯冲有关;然而,来自华北克拉通冀东地体的～2.5 Ga 麻粒岩被认为形成于垂直的沉落 (sagduction) 构造体制 (Duan et al.,2017;Liu and Wei,2018;Liu et al.,2020,2022a)。对太古宙麻粒岩相表壳岩变质演化的研究可以为研究地球早期的构造体制提供新的认识。

同世界上其他克拉通一样,华北克拉通具有漫长且复杂的演化历史 (Zhai and Santosh,2011;郑建平,2020;万渝生等,2022)。很多研究表明华北克拉通东部陆块太古宙表壳岩在新太古代晚期 (～2.50 Ga) 发生了角闪岩相-麻粒岩相变质作用,关于其形成的大地构造环境有三种不同的认识: 与微陆块拼贴时发生的俯冲-碰撞有关 (Zhai and Santosh,2013;翟明国,2019);与板块俯冲引起的碰撞造山有关 (Kusky et al.,2016);与地幔柱有关 (Geng et al.,2006;Kwan et al.,2016)。这些表壳岩记录了不同的峰期条件和P-T演化轨迹。其中角闪岩相表壳岩记录顺时针: (如Liu et al.,2020) 或逆时针型P-T轨迹 (如: Zhao,1995;Wu and Wei,2021);而麻粒岩相表壳岩一般记录逆时针型P-T轨迹 (如: Duan et al.,2017;Zhang et al.,2019;Lu and Wei,2020)。因此,需要对太古宙表壳岩进行更多的研究,以明晰其变质作用演化,并进一步探讨华北克拉通太古宙克拉通的构造体制。

清原地体位于华北克拉通北缘东段,其表壳岩在新太古代晚期经历了角闪岩相-麻粒岩相变质作用 (万渝生等,2005a,2005b;Peng et al.,2015;Wu et al.,2016;Li and Wei,2017;Wu and Wei,2021)。大部分表壳岩都经历了角闪岩相变质作用,记录了逆时针型P-T轨迹,峰期温压为0.7～1.0 GPa/750～890 ℃ (Wu et al.,2013;Wu and Wei,2021;Liu et al.,2022b)。只有东井沟、线金厂、唐力和平岭后地区出露着一些麻粒岩相表壳岩,但尚未对其变质演化进行研究。此外,清原地体新太古代晚期的构造模式也存在诸多争议,包括和地幔柱有关的垂向构造模型 (Zhai et al.,1985;Wu et al.,2013,2016;Li and Wei,2017),和洋壳俯冲有关的横向构造模型 (万渝生等,2005a;Peng et al.,2015;Wang et al.,2016a,2016b;王康等,2018;袁玲玲等,2020),以及太古宙特有的垂向/沉落构造体制 (Wu and Wei,2021)。沉落过程是指在花岗质岩浆底辟上涌的同时表壳岩下坠,最终形成穹隆-龙骨构造 (Collins et al.,1998;François et al.,2014)。为了进一步探讨清原地体新太古代晚期的构造体制,特别是在麻粒岩相变质演化的背景下,需要进行更多的研究。文章通过对清原地体东井沟地区基性麻粒岩的系统岩石学、矿物化学、相平衡模拟和锆石年代学研究,约束了其变质作用和构造演化。

1 地质背景

华北克拉通的构造演化是近年来研究的热点,主要有三种不同的构造模型。第一种模型认为,华北克拉通包括四个太古宙地块: 阴山、鄂尔多斯、龙岗和狼林地块 (图1a)。阴山地块和鄂尔多斯地块在～1.95 Ga 拼合形成了西部陆块和孔兹岩带,龙岗地块和朗林地块在～1.90 Ga 拼合形成了东部陆块和胶-辽-吉带,最后东、西陆块在～1.85 Ga 拼合形成了华北中部造山带 (TNCO) 和整个华北克拉通 (Zhao et al.,2005,2012)。第二种模型认为,在新太古代末期 (～2.5 Ga) 微陆块和岛弧拼贴形成了华北克拉通 (Zhai et al.,2000,2005;Zhai and Santosh,2011),随后在古元古代(约2.0～1.82 Ga) 经历了裂谷-俯冲-增生-碰撞等一系列构造过程,形成了丰镇带、晋豫带、胶-辽-吉带。第三种模型认为华北克拉通是由两个太古宙陆块 (东部、西部陆块) 在2.5～2.4 Ga 发生碰撞形成 (Kusky and Li,2003;Kusky et al.,2016),并且华北克拉通北缘在2.3～1.9 Ga 转变为安第斯型大陆边缘。随后,华北克拉通在1.90～1.85 Ga 与未知大陆发生碰撞,形成华北克拉通北缘的内蒙-冀北造山带 (Kusky et al.,2016)。

图1 华北克拉通构造分区图及清原地体地质简图a—华北克拉通构造分区 (据Zhao et al.,2005 修改);b—清原地体地质简图 (据万渝生等,2005b;Peng et al.,2015;Li and Wei,2017修改)Fig.1 Tectonic subdivision of the North China Craton and geological sketch map of the Qingyuan Terrane(a) Tectonic subdivision of the North China Craton (modified from Zhao et al.,2005);(b) Geological sketch map of the Qingyuan Terrane (modified from Wan et al.,2005b;Peng et al.,2015;Li and Wei,2017)

辽北清原地体地处龙岗陆块的东北部 (图1a),具有典型的 “穹隆-龙骨” 构造。清原地体主要由TTG 质 (英云闪长质-奥长花岗质-花岗闪长质) -花岗质片麻岩和少量表壳岩组成,TTG质-花岗质片麻岩呈穹隆状,而表壳岩则以皮筏状或带状处于穹隆内或穹隆之间 (图1b)。清原地体TTG 质-花岗质片麻岩中岩浆锆石的年龄显示了其侵位时间为2.57～2.49 Ga (万渝生等,2005b;白翔等,2014;Peng et al.,2015;Wang et al.,2016a,2016b;Wu et al.,2016;王康等,2018;Li et al.,2020;袁玲玲等,2020),并且它们和表壳岩一样经历了角闪岩相-麻粒岩相变质作用,变质年龄为2.55～2.42 Ga (白翔等,2014;Wang et al.,2016a,2016b;Li et al.,2020)。清原地体的表壳岩由下至上可划分为石棚子组、红透山组和南天门组 (Zhai et al.,1985)。石棚子组以麻粒岩、斜长角闪岩、辉石角闪岩、黑云母斜长片麻岩、石榴方辉石岩等为主,可能的原岩有科马提岩、拉斑玄武岩、英安岩和富铁沉积物;红透山组主要由斜长角闪岩、黑云母斜长片麻岩等组成,含少量夕线石或蓝晶石片麻岩,其原岩可能为拉斑岩、英安岩、泥质-杂砂质沉积岩;南天门组以黑云斜长片麻岩、二云母片岩、斜长角闪岩、磁铁石英岩和大理岩为主,其原岩可能为凝灰岩、泥质-杂砂质沉积岩、基性火山岩、条带状铁质建造以及碳酸盐岩。锆石年代学显示,表壳岩的变质年龄约为2.55～2.43 Ga (万渝生等,2005b;白翔等,2014;Peng et al.,2015;Wang et al.,2016a,2016b;Li and Wei,2017;Wu and Wei,2021)。Li and Wei (2017) 将红透山组的玄武质岩石按照稀土元素特征划分为富集型和亏损型两个系列,并推测其成因与地幔柱活动有关。其中富集型玄武质岩石稀土元素总量高、富集轻稀土元素,具有岛弧型地球化学特征;而亏损型玄武质岩石稀土元素总量较低。

清原地体还有很多古元古代变质基性岩墙(图2),它们侵入了新太古代TTG 质片麻岩和表壳岩。Duan et al.(2019) 和 Wu and Wei (2021)报道称变质基性岩墙经历了角闪岩相-高压麻粒岩相变质作用,具有顺时针型P-T轨迹,峰值条件为0.9～1.1 GPa/720～780 ℃。岩墙侵位时间为～2.12 Ga,变质时间为～1.85 Ga,其变质作用可能与地壳增厚造山事件有关。

图2 显示岩性分布和样品采集点的清原地区地质图 (据Duan et al.,2019 修改)Fig.2 Geological sketch map of the Qingyuan area showing lithological distribution and sample localities (modified from Duan et al.,2019)

文中研究的基性麻粒岩样品 (19DJ07) 采自东井沟村公路旁的采石坑 (124°57′55.523″N,42°12′14.957″E;图2)。清原地体高级表壳岩包含三种岩性: 基性麻粒岩、中酸性麻粒岩以及变富铁沉积岩,三种表壳岩互层产出;TTG 质片麻岩侵入表壳岩,其内可见大小不一的表壳岩捕虏体。基性麻粒岩与变富铁沉积岩互层;基性麻粒岩部分熔融程度较高,露头上常见宽度不等的浅色脉体 (图3a)。部分基性麻粒岩的浅色体中包含石榴石和斜方辉石。基性麻粒岩的部分区域含有石榴石 (图3a、3b)。

图3 表壳岩的露头照片、基性麻粒岩19DJ07 的野外露头及显微照片Grt—石榴石;Cpx—单斜辉石;Opx—斜方辉石;Bt—黑云母;Amp—角闪石;Ilm—钛铁矿;Pl—斜长石;Qz—石英下标数字: 0—包裹体状态的第一期矿物;1—第一期矿物;2—第二期矿物a、b 标记的线条对应图4 中的矿物成分剖面a—基性麻粒岩 (含石榴石域和无石榴石域) 和富铁沉积物互层,露头见部分熔融产生的浅色脉体;b—含石榴石细带的样品19DJ07 露头照片;c—样品薄片的扫描照片,岩石分为石榴石域 (19DJ07-GD) 和无石榴石域 (19DJ07-NGD);d、e—两期石榴石显微照片对比,其中Grt2 和Qz 构成冠状后成合晶,发育两期斜方辉石 (Opx1、Opx2)、三期角闪石 (包裹体Amp0、岩石主体矿物之一Amp1 呈变晶结构、Amp2 粒度细且生长在辉石周围;f—岩石局部叠加变质组合显微特征,发育Grt2+Qz+Cpx2 的后成合晶和两期斜长石 (Pl1、Pl2);g—两期单斜辉石 (Cpx1、Cpx2) 的显微特征;h—岩石近粒变晶结构显微特征,发育被辉石包裹的Amp0 和两期斜长石 (Pl1、Pl2)Fig.3 Outcrop photo of supracrustal rock,and field outcrop photo and microphotographs of mafic granulite 19DJ07(a) Interbedding of mafic granulite (containing garnet-bearing and non-garnet-bearing domains) with iron-rich metasediments,with outcrops displaying light-colored veins formed by partial melting;(b) Outcrop photo of sample 19DJ07,which contains fine bands of garnet;(c) Scanned photos of sample thin-section,categorizing the rocks into garnet-bearing domain (19DJ07-GD) and non garnetbearing domain (19DJ07-NGD);(d and e) Comparative microphotographs of two generations of garnet,where Grt2 and Qz form coronal symplectites,and two generations of orthopyroxene (Opx1/2) and three types of amphibole (Amp0/1/2) are developed;Amp0 occurs as inclusions,Amp1 as one of the main minerals of rock exhibits blastic texture,and Amp2 is fine-grained and grows around pyroxene;(f) Microscopic characteristics of locally superimposed metamorphic combinations in rocks,showing the development of symplectites of Grt2 +Qz+Cpx2 and two generations of plagioclase (Pl1/2);(g) Microscopic characteristics of two generations of clinopyroxene (Cpx1/2);(h) Homeoblast texture in 19DJ07-NGD,with Amp0 enclosed by pyroxene and two generations of plagioclase(Pl1/2)The lines labeled with a and b corresponds to the composition zoning profile in Fig.4;mineral abbreviations are according to Whitney and Evans,2010;Warr (2021): Grt-garnet;Cpx-clinopyroxene;Opx-orthopyroxene;Bt-biotite;Amp-hornblende;Ilm-ilmenite;Pl-plagioclase;Qz-quartz

2 岩石学分析

2.1 全岩成分分析

研究样品的全岩主量元素测试分析于中国地质大学 (北京) 的相关实验室完成,使用仪器为电感耦合等离子体发射光谱仪 (ICP-OES),测试结果见表1。

表1 清原地区基性麻粒岩的全岩成分以及相平衡模拟用的有效全岩成分Table 1 Bulk-rock compositions and effective bulk-rock compositions of samples from Qingyuan terrane

2.2 岩相学和矿物化学

矿物化学分析使用了北京大学地球与空间科学学院的JXA-8100 型号的电子探针,测试条件为:加速电压15 kV,电流10 nA,束斑直径为2 μm(测试黑云母时为5 μm),使用SPI 公司的53 种矿物作为标准样品。

基性麻粒岩样品19DJ07 中含石榴石的区域呈条带状分布,条带宽度0.4～1.0 cm (图3b、3c),因此,可以划分为含石榴石域 (19DJ07-GD) 和无石榴石域 (19DJ07-NGD) 两部分,分别为石榴二辉麻粒岩和二辉麻粒岩组合。两个结构域的代表性矿物成分见表2、表3。

表3 清原地体基性麻粒岩样品19DJ07-NGD 中代表性矿物的探针分析Table 3 Selected microprobe analyses for mafic granulite 19DJ07-NGD from the Qingyuan terrane

(1) 19DJ07-GD 的矿物组合各成分体积分数φ分别为石榴石 (5%)、单斜辉石 (15%)、斜方辉石 (6%)、角闪石 (28%)、黑云母 (1%)、斜长石 (40%)、石英 (3%) 及钛铁矿+磁铁矿(2%)。具有典型等粒粒状变晶结构 (图3d—3g)。

石榴石分为Grt1和Grt2两类。Grt1为他形—半自形,粒度0.1～0.5 mm,部分颗粒含有极少的包裹体 (石英或黑云母) (图3d)。Grt1的成分中xGrs=Ca/ (Mg+Fe+Ca+Mn),xAlm(铁铝榴石摩尔分数)、xSps(锰铝榴石摩尔分数)、xPy(镁铝榴石摩尔分数) 依此类推,从核部到边部逐渐降低(0.26→0.21);xPy具有与xGrs相似的环带 (0.13→0.10);xAlm表现出与xPy相反的环带,从核部到边部升高 (0.59 →0.65);xSps只有0.02～0.03,整体较为均匀 (图4a)。Grt2作为0.1～0.2 mm 宽的冠状环带生长在其他矿物周围 (图3d、3g);Grt2呈筛状,常与细粒石英包裹体及Cpx2构成后成合晶结构 (图3f)。相比于Grt1,Grt2的xGrs更低,为0.19～0.24,但是xAlm更高,为0.60～0.68;xPy为0.09～0.13,xSps为0.02～0.03。

图4 样品19DJ07 中代表性矿物的成分特征图a—19DJ07-GD 的Grt1 的成分剖面 (xAlm=Fe2+/ (Fe2++Mn+Mg+Ca),xSps、xPy、xGrs 依此类推);b、c—典型角闪石的Ti 成分环带 (b—19DJ07-GD 的Amp1,c—19DJ07-NGD 的Amp1);d、e—典型斜长石的的xAn 成分剖面图, xAn=Ca/ (Ca+Na+K) (d—19DJ07-GD 的Pl1,e—19DJ07-NGD 的Pl1)Fig.4 Composition diagrams of representative minerals in Sample 19DJ07(a) Grt1 chemical zoning profiles in 19DJ07-GD (xAlm=Fe2+/(Fe2+ +Mn+Mg+Ca);xSps,xPy,and xGrs follow this pattern,respectively);(b and c) Typical composition zoning of Ti for Amp1 in 19DJ07-GD (b) and 19DJ07-NGD (c);(d) Typical composition zoning of xAn(=Ca/(Ca+Na+K)) for Pl1 in 19DJ07-GD (d) and 19DJ07-NGD (e)

单斜辉石包含Cpx1和Cpx2两类。Cpx1作为主体矿物之一,粒度0.2～0.3 mm,他形粒状 (图3d、3f);而Cpx2则是0.01～0.02 mm 的他形的颗粒或冠状体,生长于辉石或角闪石周围 (图3f、3g)。Cpx1和Cpx2具有相似的透辉石成分:xWo=0.46～0.51、xEn=0.29～0.36 (Morimoto,1988)。

斜方辉石包含Opx1和Opx2两类。Opx1是粒度0.1～0.3 mm 的他形—半自形柱状颗粒 (图3d、3e);Opx2常生长于叠加变质结构域的二辉石周围,粒度～0.01 mm (图3e、3g)。两类斜方辉石都具有斜铁辉石的成分,xEn(=Mg/(Mg+Fe+Ca))=0.39～0.46、xWo(=Ca/(Mg +Fe +Ca))=0.01(Morimoto,1988)。

根据岩石结构,角闪石可以分为Amp0、Amp1和Amp2三类。Amp0粒度约为0.01～0.10 mm,作为斜方辉石或斜长石的包裹体存在 (图3e、3g);Amp1作为主体矿物之一,粒度为0.1～0.5 mm,他形—半自形 (图3d—3g);Amp2为他形,粒度为0.01～0.05 mm,生长于其他矿物边部 (图3d—3f)。根据Hawthorne et al.(2012) 的分类,三类角闪石都属于韭闪石 ((Na+K)A=0.59～0.75,(Al+Fe3++2Ti)C=0.87～1.20)。Amp0成分中Ti 含量为0.19～0.27 p.f.u.(per formula unit)、xMg=0.41～0.47。Amp1的Ti 含量为0.17～0.25 p.f.u.,且具有从核部到幔部升高再到边部降低的成分环带 (图4b);xMg=0.39～0.47。Amp2的Ti 含量最低,为0.14～0.21 p.f.u.,xMg=0.44～0.46。

黑云母也分Bt1和Bt2两类,Bt1为他形或半自形,粒度为0.01～0.1 mm;Bt2为他形,粒度为～0.01 mm (图3d、3f);两类黑云母成分相似,xMg=0.45～0.46,Ti 含量为0.27～0.28 p.f.u.。

斜长石可以分为Pl1和Pl2两类。Pl1为他形—半自形,粒度为0.1～5.0 mm (图3d、3e、3g)。Pl1的xAn从核部到幔部升高再到内边降低 (0.42→0.45→0.40),部分Pl1颗粒还生长了xAn更低的外边部 (0.26～0.36;图4d),可能为熔体局部汇聚结晶形成 (Liu et al.,2020;Zhang et al.,2020,2021)。部分Pl1颗粒的核部xAn较低,低至0.36。Pl2一般生长在叠加变质结构域,其xAn更低,且从核部到边部降低 (0.32→0.22),粒度也较Pl1更细 (0.01～0.05 mm)。

根据矿物的结构和化学成分,可以划分出三个阶段的矿物组合,即第一期变质峰期阶段、第一期变质固相线阶段以及叠加变质阶段。第一期变质峰期阶段的矿物组合为Grt1+Cpx1+Opx1+Pl1+Ilm;第一期变质固相线的矿物组合与峰期阶段相比的差别是峰后阶段伴随熔体结晶生长的角闪石、黑云母和斜长石等;叠加变质的矿物组合为Grt2+Qz±Cpx2的后成合晶以及Opx2、Pl2和Amp2。

(2) 19DJ07-NGD 的矿物组合中各成分体积分数φ分别为单斜辉石 (10%)、斜方辉石 (7%)、角闪石 (37%)、黑云母 (1%)、斜长石 (43%)、石英(1%) 及钛铁矿+磁铁矿(1%) (图3g、3h)。岩石具有麻粒岩典型的等粒变晶结构、三联点结构。

单斜辉石和斜方辉石为他形,粒度0.1～0.5 mm,内部常含角闪石包裹体 (图3g、3h)。单斜辉石中xEn=0.28～0.32、xWo=0.47～0.51,斜方辉石中xEn=0.39～0.41、xWo=0.01,根 据 Morimoto(1988) 的分类,它们分别为透辉石和斜铁辉石。

角闪石可以分为三类 (图3h): Amp0、Amp1和Amp2。Amp0是斜方辉石或斜长石中的包裹体,Ti 含量为0.18～0.24 p.f.u.,xMg=0.38～0.44。Amp1粒度为0.1～0.5 mm,他形—半自形,Ti 含量为0.17～0.29 p.f.u.,且呈现从核部到幔部升高再到边部降低的成分环带 (图4c);xMg=0.37～0.47。Amp2是粒度～0.01 mm 的他形颗粒,生长在辉石周围,Ti 含量为0.14～0.24 p.f.u.,xMg=0.42～0.48。三类角闪 石的 (Na+K)A=0.40～0.75、(Al+Fe3++2Ti)C=0.49～1.32,根据Hawthorne et al.(2012) 的分类,大部分颗粒属于韭闪石,少部分颗粒属于透闪石、镁质闪石及钙镁闪石。黑云母为0.05～0.1 mm 片状,xMg=0.42,Ti 含量为0.27 p.f.u.。

斜长石有Pl1、Pl2两类。Pl1为他形,粒度0.1～0.5 mm (图3g、3h)。Pl1的xAn从核部到幔部升高再在到边部降低 (图4e)。岩石中Pl1核部-幔部的xAn一般集中在0.42～0.45,内边部xAn集中在0.40 左右;但是部分颗粒的核部xAn可以低至0.35。大部分Pl1颗粒外边部的xAn更低(0.25～0.33),可能代表着局部的熔体汇聚、结晶(如: Liu et al.,2020;Zhang et al.,2020,2021)。Pl2为他形颗粒,粒度～0.1 mm,成分XAn从核部到边部降低 (0.35→0.24)。

根据观察到的岩石结构及矿物成分,可以划分出三个阶段的矿物组合: 第一期变质峰期阶段矿物组合Opx+Cpx+Pl1+Ilm±Amp、以峰后冷却阶段生长角闪石、斜长石为特征的固相线组合Opx+Cpx+Pl1+Amp1+Bt+Ilm+Qz 以及叠加变质矿物Pl2。包裹在二辉石中的Amp0也许代表着进入变质阶段。

3 相平衡模拟

基于上述样品的矿物组合及矿物化学成分,在NCKFMASHTO (Na2O-CaO-K2O-FeO-MgO-Al2O3-SiO2-H2O-TiO2-Fe2O3) 体系中,利用 THERMOCALC 3.45 (Holland and Powell,1998) 和热力学数据库ds62 (Holland and Powell,2011) 进行了视剖面图计算。所用矿物相活度-成分 (a-x) 模型包括:石榴石 (White et al.,2007,2014),角闪石、单斜辉石、熔体 (Green et al.,2016),斜方辉石、黑云母 (White et al.,2014),斜长石 (Holland and Powell,2003),钛铁矿 (White et al.,2000),石英为纯端元组分。模拟中所使用的H2O 和O 值是依据T-M (H2O/O) 图解来确定的,以确保模拟的固相线组合与实际观察的一致 (Korhonen et al.,2013),并且这样确定的O 值可避免样品制备过程中潜在的污染或氧化以及在抬升出露过程中可能的流体作用导致的氧逸度的变化 (如: Dos Santos et al.,2011;Cao et al.,2011;Korhonen et al.,2012)。

为了使视剖面图更好地反映样品的实际情况,采用有效全岩成分的摩尔分数进行相平衡模拟。有效全岩是据样品的矿物组合、矿物的体积分数及化学成分进行整合计算得到的。所有样品的有效全岩均列在了表1 中。

19DJ07-GD 的P-T视剖面图 (图5) 中,P-T范围是0.2～1.2 GPa/700～1000 ℃。在相关的矿物组合区域绘制出了石榴石的xGrs等值线、斜长石的xAn等值线以及角闪石的Ti 含量等值线。固相线组合Grt1+Opx1+Cpx1+Amp1+Bt1+Pl1+Qz±Ilm 稳定的温压范围是0.75～0.92 GPa/780～810 ℃,推测的峰期组合Grt1+Opx1+Cpx1+Pl1+Ilm+Qz 可以给出一个更大的温压范围0.75～0.9 GPa/860～980 ℃。Amp0中最高的Ti 含量 (0.27 p.f.u.) 对应着大约900 ℃的温度条件;将实测Pl1的xAn(核部-幔部-边部: 0.42→0.45→0.36) 投到图中,再结合推测的峰期矿物组合可以确定峰期条件为0.8～0.9 GPa/900～940 ℃。依据Pl1的xAn从核部到幔部再到内边部先升高后降低的成分剖面 (图4d)可以给出一条逆时针的变质轨迹,且最小值0.36正好和固相线组合预测的等值线数值一致。xAn=0.25～0.33 的外边部则代表着亚固相线条件下的更低温压条件,可能由局部熔体结晶形成。Amp1从核部到幔部升高再到边部降低的Ti 成分环带(0.17～0.25 p.f.u.;图4b) 同样对应着逆时针的演化轨迹。图中角闪石最高Ti 含量投影超过了角闪石的最大稳定域,可能原因是角闪石稳定域被低估,如果模拟中考虑F、Cl 等组分,角闪石的稳定域会扩大 (Sajeev et al.,2009)。所以19DJ07-GD 的峰期温度一定在900 ℃以上。或许由于受到熔体在某些区域汇聚而产生局部平衡的影响,实测的石榴石的成分并不能投到图5 中,图中预测的xGrs最大值 (0.20) 小于实测的xGrs最小值(0.21)。在逆时针演化的峰后降温降压过程中,斜长石、角闪石、黑云母以及石英都会生长。

图5 样品19DJ07-GD 在NCKFMASHTO 体系下的P-T 视剖面图及P-T 轨迹计算P-T 视剖面图所使用的有效全岩成分见表1Grt—石榴石;Cpx—单斜辉石;Opx—斜方辉石;Bt—黑云母;Amp—角闪石;Ilm—钛铁矿;Pl—斜长石;Qz—石英;L—熔体Fig.5 P-T pseudosection with proposed P-T vectors for 19DJ07-GD in the NCKFMASHTO systemThe effective bulk-rock compositions used to calculate the P-T pseudosection are shown in Table 1.Grt-garnet;Cpx-clinopyroxene;Opx-orthopyroxene;Bt-biotite;Amp-hornblende;Ilm-ilmenite;Pl-plagioclase;Qz-quartz;L-melt

19DJ07-NGD 的P-T视剖面图 (图6) 中,在相关的矿物组合区域绘制出了斜长石的xAn等值线以及角闪石的Ti 等值线。观察到的固相线组合Cpx+Opx+Amp1+Pl1+Ilm+Bt+Qz 在视剖面图中被预测稳定在一个狭窄的P-T区域0.3～0.9 GPa/740～830 ℃,推测的峰期矿物组合Opx+Cpx+Pl1+Ilm±Amp 稳定在<0.9 GPa/>840 ℃的宽泛区域。依据Pl1幔部的最大xAn(0.45) 和Amp1最高Ti 值以确定峰期温压条件为～0.9 GPa/>950 ℃,并位于含石榴石的组合中。Pl1的xAn从幔部到外边部降低(0.45→0.40) 对应着峰后降温降压的演化阶段,岩石中实测的Pl1最低xAn=0.35 和固相线区域预测的数值一致。在峰后降温及轻微降压的过程中,石榴石会通过Grt+Cpx+Ilm+L=Opx+Amp+Pl 的变质反应被消耗殆尽 (图6)。沿着这条降压冷却轨迹,相图中的角闪石的Ti 值、斜长石的xAn降低。实测中从核部到幔部升高再到内边部降低的Pl1的xAn成分环带 (图4e) 以及从核部到幔部升高再到边部降低的Amp1的Ti 成分环带 (图4c) 都对应着一条逆时针的演化轨迹 (图6)。在图6 中XAn<0.35 的等值线将投到亚固相线的条件下 (Pl1的外边部xAn=0.24～0.35),反映其可能生长于由熔体汇聚形成的局部区域,并且记录了比图中固相线更低的温度 (如: Liu et al.,2020;Zhang et al.,2021)。

图6 样品19DJ07-NGD 在NCKFMASHTO 体系下的P-T 视剖面图及P-T 轨迹计算P-T 视剖面图所使用的有效全岩成分见表1Grt—石榴石;Cpx—单斜辉石;Opx—斜方辉石;Bt—黑云母;Amp—角闪石;Ilm—钛铁矿;Pl—斜长石;Qz—石英;Mt—磁铁矿;L—熔体Fig.6 P-T pseudosection with proposed P-T vectors for 19DJ07-NGD in the NCKFMASHTO systemThe effective bulk-rock compositions used to calculate the P-T pseudosection are shown in Table 1.Grt-garnet;Cpx-clinopyroxene;Opx-orthopyroxene;Bt-biotite;Amp-hornblende;Ilm-ilmenite;Pl-plagioclase;Qz-quartz;Mtmagnetite;L-melt

4 锆石年代学

19DJ07 的锆石U-Pb 年代学测试及微量元素分析在北京大学造山带与地壳演化教育部重点实验室完成。测试的16 个标准锆石91500 的207Pb/206Pb加权平均年龄为1066 ± 24 Ma (2SE;MSWD=0.19),和参考的TIMS 年龄 (1065.4±0.3 Ma;2SE) 几乎一致 (Wiedenbeck et al.,1995),误差在2σ以内。测试的9 个Plešovice 锆石的206Pb/238U加权平均年龄为336.7±2.3 Ma (2SE;MSWD=0.76),和参考 的TIMS 年龄 (337.1 ± 0.4 Ma;2SE) 几乎一致 (Sláma et al.,2008)。所有19DJ07锆石的测试结果详细数据见表4。

表4 清原地体基性麻粒岩样品19DJ07 的锆石U-Pb 同位素数据Table 4 U-Pb isotopic data for zircon of mafic granulite 19DJ07 from the Qingyuan terrane

锆石颗粒大部分为长柱状,少数为浑圆状,直径10～100 μm,长宽比为1 ∶1～5 ∶1。锆石的CL 图像显示了其内部结构为补丁状、冷杉叶状,少数具有震荡环带或较均匀,少数锆石具有核-边结构。锆石的颜色有亮白色、浅灰色、深灰色或黑色 (图7a)。57 个U-Pb 同位素分析点位显示207Pb/206Pb 表观年龄范围为2508±18 Ma 到2497±10 Ma,加权平均年龄为2494.7±4.4 Ma (MSWD=1.7),谐和图上交点年龄为2498±6.9 Ma (MSWD=0.39) (图7b)。锆石Th/U=0.11～0.46,(Lu/Gd)N=41.4～111.9;大部分锆石具有负Eu 异常,极个别为正Eu 异常 (Eu/Eu∗=0.29～1.3);重稀土元素呈现上升的配分模式 (图7c)。上述锆石结构和化学特征表明,锆石可能生长于麻粒岩相变质过程 (Vavra et al.,1999;Rubatto,2002;Corfu et al.,2003)。

图7 石榴二辉麻粒岩样品19DJ07 的锆石分析结果a—代表性锆石颗粒的阴极发光图像,显示其内部结构和分析位置,编号详情见表4;b—通过LA-ICP-MS 分析结果计算的U-Pb 同位素谐和图、加权平均年龄图和U-Pb 年龄分布图;c—球粒陨石标准化的稀土配分模式图 (根据Sun and McDonough,1989 进行标准化)Fig.7 Analytical results of zircons from mafic granulite 19DJ07(a) Cathodoluminescence images of selected zircon grains showing the inner structures and analyzed locations with the identification numbers as in Table 4;(b) U-Pb isotopic concordia,weighted mean age diagram,and U-Pb age distribution diagram showing the LAICP-MS analytical results;(c) Chondrite-normalized REE patterns of zircons (normalized according to Sun and McDonough,1989)

5 讨论

5.1 变质演化过程

基于清原地体基性麻粒岩的岩相学特征和相平衡模拟,划分了三个变质阶段,包括峰前升温升压至峰期阶段、峰后减压冷却阶段和后期叠加变质阶段。

峰期前的升温升压阶段主要依据矿物成分的变化趋势确定。如图4b、4c 所示,无论19DJ07-GD 还是-NGD,Amp1的Ti 成分环带都呈现出从核部到幔部升高再到边部降低的成分环带。两个区域Amp1幔部及Amp0核部的最高Ti 值 (0.27～0.29) 指示峰期温度为900～950 ℃,从核部到幔部升高的Ti 值预示着峰期前的升温过程。同样,如图4d、4e 所示,Pl1从核部到幔部升高的xAn成分环带也对应着升温升压的进变质过程 (图5,图6)。依据所推测的岩石峰期矿物组合 (-GD 和-NGD 均为Grt+Cpx+Opx+Pl+Ilm ± Amp) 及相应的矿物等值线推测岩石峰期的压力为0.8～0.9 GPa。

19DJ07-GD 和-NGD 峰后减压冷却过程均以Grt+Cpx+L=Opx+Pl+Amp 的反应为主。在此过程中石榴石被消耗而斜长石、角闪石生长;斜长石幔部到边部的xAn降低 (图4d、4e),角闪石幔部到边部的Ti 含量降低 (图4b、4c)。图5、图6 的模拟结果表明,类似的演化过程可以使-NGD 的石榴石被完全消耗,而-GD 的石榴石存留下来,因为后者比前者具有更低的Mg#值,更有利于石榴石的稳定。如图5、图6 所示,降温降压过程终止于固相线,含有黑云母的固相线组合给出了0.70～0.78 GPa/～800 ℃的P-T条件。

19DJ07-GD 和-NGD 都发育了叠加变质组合,尤其以-GD 更为丰富。叠加变质组合以Grt2+Qz±Cpx2的后成合晶/冠状体及第二期矿物Opx2/Pl2为特征 (图3f—3h),和岩石主体矿物的细粒变晶结构有明显差别。虽然19DJ07 没有古元古代的年龄记录,但是Grt2+Qz ± Cpx2后成合晶和区内经历了古元古代晚期 (～1.85 Ga) 高压麻粒岩相变质的基性岩墙中石榴石+石英± 单斜辉石后成合晶一致。Duan et al.(2019) 认为该后成合晶组合记录了顺时针型P-T轨迹,可能与古元古代华北克拉通北缘的地壳增厚造山事件有关。在冀东地体诸多记录了新太古代—古元古代两期变质年龄的表壳岩中都发现了类似的后成合晶叠加结构 (如:Yang and Wei,2017;Lu and Wei,2020)。Yang and Wei (2017) 通过石榴子石-全岩Lu-Hf 等时线方法得到了基性麻粒岩中和Qz 构成冠状后成合晶的Grt2年龄为1.79～1.77 Ga。在清原地体发育着丰富石榴石+石英后成合晶结构的中酸性麻粒岩中也发现了两期的年龄记录。于是认为基性麻粒岩19DJ07 同样经历了古元古代高压麻粒岩相变质的叠加,并且保留了形成于缺流体的升温升压阶段的冠状体/后成合晶结构。

5.2 年代学启示

基性麻粒岩19DJ07 的锆石记录了新太古代末期变质年龄位2508±18～2497±10 Ma,加权平均年龄2494.7±4.4 Ma (MSWD=1.7) (图7)。变质锆石的生长主要受熔体/流体行为控制。Zr 在熔体中有较高的溶解度 (Watson and Harrison,1984),所以锆石在含有大量熔体的进变质过程中不会生长(Roberts and Finger,1997),而会在峰后降温冷却过程中随着熔体结晶生长。所以高级变质岩中的变质锆石一般记录峰后冷却过程的时间 (Kelsey and Powell,2011;Zhang et al.,2013;Wei et al.,2014;Yakymchuk et al.,2014,2017)。如图8 所示,统计已有研究在该地区的工作,TTG 质-花岗质片麻岩中岩浆锆石的年龄为2.57～2.49 Ga,代表了其侵位时间;角闪岩相表壳岩的变质年龄为2.55～2.43 Ga,并且在2.51 Ga、2.49 Ga、2.47 Ga出现峰期;麻粒岩相表壳岩的变质年龄为2.50～2.46 Ga,并且在2.49 Ga、2.47 Ga 出现峰期。所以19DJ07 获得的新太古代末期的变质年龄与清原地体其它角闪岩相-麻粒岩相表壳岩变质年龄基本一致,并与区内TTG 质-花岗质岩浆活动晚期脉冲年龄相同 (图8)。至于岩石并未记录古元古代变质年龄的原因有以下三种: 古元古代叠加变质事件发生于缺流体的条件下,不利于锆石生长;第二期叠加变质事件峰期温度较第一期变质事件更低,因此无法改造第一期变质事件中形成的锆石;“Oswald 熟化” 表明大颗粒的表面能更低,因此第一期较高温度的变质事件形成的粗粒锆石很难被改造而形成小颗粒锆石 (Nemchin et al.,2001)。

图8 清原地区各类岩石的年龄汇总图图中展示了其他角闪岩相-麻粒岩相表壳岩的年龄 (万渝生等,2005b;白翔等,2014;Peng et al.,2015;Wang et al.,2016a;Wu et al.,2016;Li and Wei,2017;Wu and Wei,2021)、TTG 质-花岗质片麻岩的结晶年龄 (万渝生等,2005b;Bai et al.,2014;Peng et al.,2015;Wang et al.,2016a,2016b;Wu et al.,2016;王康等,2018;Li et al.,2020;袁玲玲等,2020) 以及研究样品19DJ07的年龄Fig.8 Age summary chart of different rock types from the Qingyuan TerraneThe ages of other amphibolite-granulite facies supracrustal rocks(Wan et al.,2005b;Bai et al.,2014;Peng et al.,2015;Wang et al.,2016;Wu et al.,2016;Li and Wei,2017;Wu and Wei,2021) and the crystallization ages of TTG-granitic gneiss (Wan et al.,2005b;Bai et al.,2014;Peng et al.,2015;Wang et al.,2016a,2016b;Wu et al.,2016;Wang et al.,2018;Li et al.,2020;Yuan et al.,2020) are shown in the diagram.The age data for the investigated samples 19DJ07 are presented for comparison with other research data

5.3 大地构造意义

基性麻粒岩19DJ07 记录了逆时针的P-T轨迹,包括峰前升温升压至峰期和峰后降温降压至固相线的变质阶段 (图9),峰期变质条件为0.8～0.9 GPa/900～950 ℃,对应地温梯度30 ℃/km。清原地体基性麻粒岩的P-T轨迹与冀东新太古代麻粒岩相表壳岩几乎一致,只是峰期温度略低 (图9;Liu and Wei,2018,2020)。清原地体具有典型的太古宙 “穹隆-龙骨构造” (图1;Collins et al.,1998;Hickman,2004;Lin and Beakhouse,2013;Anhaeusser,2014),并且表壳岩变质时间与TTG质-花岗质岩浆活动的晚期脉冲基本相同 (图8),由此推断麻粒岩相表壳岩变质作用演化与太古宙特有的垂向/沉落构造体制有关 (如: Duan et al.,2017;Liu and Wei,2018,2020;Liu et al.,2020;Wu and Wei,2021)。沉落构造体制包括以下过程:首先地壳浅部表壳岩被侵入的高温TTG 质-花岗质岩浆加热;然后受热的表壳岩碎裂成大小不一的块体,坠入巨大 “TTG 岩浆海” 的深部,对应峰前升压过程;最后TTG 岩浆/岩石穹窿抬升并冷却,裹挟其中的表壳岩块体发生峰后降压降温的变质演化 (Duan et al.,2017;Liu and Wei,2018)。

图9 清原地区和冀东地区新太古代表壳岩的P-T-t 轨迹总结清原地区: 19DJ07 为文中研究样品;Wu21a 为Wu and Wei (2021) 报道的新太古代石榴角闪岩。冀东地区: LW18 和LW20 为Liu and Wei(2018,2020) 报道的新太古代麻粒岩相表壳岩;Liu20 为Liu et al.(2020) 报道的新太古代角闪岩相表壳岩。基性岩的湿固相线 (WS) 来自Lambert and Wyllie (1972)。角闪石的消失线: Amp-out-a 来自Wyllie and Wolf (1993),Amp-out-b 来自Sen and Dunn (1994);石榴石生成线: Grt-in-a 来自Winther and Newton (1991),Grt-in-b 来自Liu et al.(1996)。变质相界线来自魏春景等 (2017),特别指出麻粒岩相包括正常麻粒岩亚相 (<900 ℃) 和超高温麻粒岩亚相 (>900 ℃) (Brown,2007)Fig.9 A summary of the P-T-t paths for Neoarchean supracrustal rocks from the Qingyuan Terrane with comparison to those from the Eastern Hebei TerraneFor the Qingyuan terrane: 19DJ07,the investigated sample of this study;Wu21a,Neoarchean garnet amphibolite in Wu and Wei(2021);LW18 and LW20,granulite facies supracrustal rocks in Liu and Wei (2018,2020);Liu20,amphibolite facies supracrustal rock in Liu et al.(2020).The wet solidus of mafic rocks (WS) is after Lambert and Wyllie (1972);the amphibole-out curves:Amp-out-a is from Wyllie and Wolf (1993),and Amp-out-b is from Sen and Dunn (1994);the garnet-in curves: Gtr-in-a is from Winther and Newton (1991),and Gtr-in-b from Liu et al.(1996).The distributions of metamorphic facies are from Wei et al.(2017) It is pointed out that the granulite facies include normal granulite subfacies (<900 ℃) and UHT granulite subfacies(>900 ℃)

5 条新太古代晚期表壳岩的变质演化P-T轨迹(19DJ07、LW18、LW20、Wu21a 和Liu20,图9)显示出虽然其峰前升压加热过程以及峰期温度不同,但峰期压力却很相似,都对应大约30～35 km地壳深度。这说明表壳岩在峰期前经历了不同的加热过程,但块体均沉落到了下地壳深度。事实上,华北克拉通东部陆块的角闪岩相-麻粒岩相表壳岩记录着十分复杂的P-T轨迹。以冀东地区为例,太古宙表壳岩-TTG 片麻岩同样呈 “穹隆-龙骨” 构造,其中新太古代角闪岩相表壳岩记录了顺时针型P-T轨迹 (Liu et al.,2020),而新太古代麻粒岩相表壳岩则记录着逆时针型P-T轨迹(Liu and Wei,2018,2020)。针对华北克拉通东部陆块新太古代表壳岩具有时间上共存的不同P-T-t轨迹的现象,Yu et al.(2022) 进行了一系列数值模拟研究,认为这很可能是在和地幔柱相关的地球动力学机制下发展起来的。

6 结论

通过对清原地区基性麻粒岩进行的详细的岩石学、地球化学、相平衡模拟以及锆石年代学的研究,主要得出以下几点认识。

(1) 清原地体基性麻粒岩记录了逆时针的PT轨迹,峰期温压条件为0.8～0.9 GPa/900～950 ℃,峰前以升温升压为主,峰后为降温降压演化。

(2) 锆石定年结果表明清原地体新太古代晚期高温—超高温麻粒岩相变质作用发生在～2.50 Ga,与区域上TTG 质-花岗质岩浆活动的晚期脉冲时间一致;结合区域上的穹隆-龙骨构造、基性麻粒岩逆时针的P-T轨迹,推测清原新太古宙麻粒岩相变质作用形成于太古宙特有的垂向构造/沉落构造体制。

(3) 清原地体的表壳岩都经历了古元古代高压麻粒岩相变质作用的叠加,特征为岩石局部发育后成合晶和冠状体结构。

致谢:感谢秦红、李小犁、马芳、武现伟对实验分析的帮助,感谢段站站、吴定、董杰、杨子珍对野外工作的参与帮助。