黑龙江虎林杂岩变形样式与时代:对中国东北东部早白垩世古太平洋板块俯冲的启示*
2020-04-01梁琛岳刘永江宋志伟张骞杨岩
梁琛岳 刘永江 宋志伟 张骞 杨岩
1. 吉林大学地球科学学院,长春 1300612. 吉林大学东北亚矿产资源评价自然资源部重点实验室,长春 1300613. 中国海洋大学海底科学与探测技术教育部重点实验室,海洋高等研究院,海洋地球科学学院,青岛 2661004. 青岛海洋科学与技术国家实验室海洋矿产资源评价与探测技术功能实验室,青岛 2662371.
黑龙江虎林地区位于兴凯地块北部,出露一套呈NE走向展布的花岗质-变质岩体,前人称之为虎林杂岩(图2; Zhouetal., 2010; Wildeetal., 2010; 胥嘉等, 2012; 何松, 2017),主要由绿片岩-蓝片岩相变质的“下元古界黑龙江群”,少量元古代混合花岗岩,以及呈不规则岩基状侵入二者之中的中生代花岗质岩体共同组成,岩体内部局部可见上新统玄武岩不整合覆盖其上(图2、图3; 黑龙江省地质矿产局, 1993; Wuetal., 2007; Zhouetal., 2010; 胥嘉等, 2012; 何松, 2017)。部分学者对虎林杂岩的年代学和地球化学属性进行了研究(Wildeetal., 2010; Zhouetal., 2010; 胥嘉等, 2012; 何松, 2017),初步认为该套杂岩由一套副变质岩云母片岩(~250Ma, Zhouetal., 2010)和中生代侵入岩体(112~128Ma, Wildeetal., 2010; 胥嘉等, 2012; 何松, 2017)组成,地球化学数据表明,杂岩内早白垩世岩体形成与古太平洋板块俯冲诱发的岩石圈减薄拆沉引发的伸展背景相关(胥嘉等, 2012; 何松, 2017)。但虎林杂岩内部发育多期次花岗质岩脉侵入,并且岩石强烈韧性变形,具有较为统一的片麻理,迄今为止,尚未有学者对虎林杂岩的变形样式及其动力学背景进行探讨,并且后期(同)变形花岗质岩脉的时代和性质也尚不明确。
图1 中国东北东部大地构造图(a,据Liu et al., 2017a修改)及虎林地区地质简图(b,据Zhou et al., 2010修改)
图2 虎林地区地质图(a, 据Zhou et al., 2010; 何松, 2017修改)及构造剖面(b)和构造要素赤平投影图(c, 下半球投影图)
图3 宏-微观构造变形特征
本文聚焦于虎林盆地中央隆起区出露的虎林杂岩(图1、图2),对其变形样式和构造特征进行解析,同时结合岩石学和岩相学研究,开展变形岩脉的LA-ICP-MS锆石U-Pb同位素年代学、变形岩石的石英EBSD组构测试等多手段综合分析,结合前人年代学研究成果以及与虎林盆地的耦合关系,探讨其变形机制和动力学背景,为进一步研究兴凯地块乃至中国东北东部中生代构造属性和古太平洋构造体制下构造演化过程提供理论支持。
1 区域地质背景
兴凯地块主体位于俄罗斯境内,少量出露于我国境内,其北部与完达山增生地体相邻,西部以敦化-密山断裂为界与佳木斯地块相接,西南以牡丹江-依兰缝合带为界与松嫩地块相邻(Liuetal., 2017a, b; 周建波等, 2018; Zhouetal., 2018; 刘永江等, 2019; 图1),是兴蒙造山带的重要组成部分(Zhouetal., 2018; 许文良等, 2019; 刘永江等, 2019)。兴蒙造山带(古亚洲洋东段)在古生代期间,额尔古纳地块、兴安地块、松嫩地块与华北板块相继拼合,并于二叠纪末期-中三叠世最终闭合(Wuetal., 2011; Liuetal., 2017a; 刘永江等, 2019; Zhangetal., 2019; Liangetal., 2019b)。之后晚三叠世开始,进入古太平洋的演化阶段,持续性俯冲导致佳木斯-兴凯地块和松嫩地块拼合(Wildeetal., 2000, 2010; Zhouetal., 2010; Wuetal., 2007, 2011; Liuetal., 2017a, b; 周建波等, 2016, 2018; 张夺等, 2018; 刘永江等, 2019),以及包括跃进山杂岩、那丹哈达增生杂岩以及饶河杂岩在内的增生杂岩系统以及大规模的岩浆事件和成矿作用(黑龙江省地质矿产局, 1993; Wuetal., 2011; 孙明道, 2016; Bietal., 2016, 2017; Wangetal., 2017; 李三忠等, 2017; 周建波等, 2016, 2018; 张夺等, 2018; Lietal., 2019; 许文良等, 2019)。因此,理清楚兴凯地块内部虎林地区出露的虎林杂岩的构造变形样式和大地构造背景,对于理解古太平板块的俯冲过程具有重要的意义。
兴凯地块的变质基底主要沿虎林-虎头一线展布(黑龙江省地质矿产局, 1993; Zhouetal., 2010; Wildeetal., 2010)。虎头地区分布的虎头杂岩为一套绿片岩-蓝片岩相变质杂岩,主体为一套副片麻岩以石榴花岗片麻岩为主,记录了490~510Ma的早古生代岩浆和变质事件(Wildeetal., 2010; Zhouetal., 2010),与佳木斯地基底块麻山群年龄较为一致,同属中国东北~500Ma泛非期孔兹岩带的一部分,指示兴凯地块与佳木斯地块具有相同的泛非期变质基底(周建波等, 2011; Wildeetal., 2010; Zhouetal., 2010)。虎林盆地中央隆起带内出露的虎林杂岩主要由云母片岩和侵入其中的花岗质岩体组成,黑龙江省地质矿产局(1993)将该套云母片岩划分为下元古界黑龙江群,岩体遭受强烈糜棱岩化,片理发育,Zhouetal. (2010)在云母片岩中获得233~270Ma的碎屑锆石年龄,指示其原岩沉积时代不早于233Ma(图2a)。侵入其中的中生代岩体局部变形强烈,部分学者的年代学工作证实,该时期花岗质岩体/脉侵位时间集中在112~128Ma(Wildeetal., 2010; 胥嘉等, 2012; 何松, 2017; 图2a),可能与古太平洋板块俯冲诱发的岩石圈减薄拆沉引发的伸展背景相关(胥嘉等, 2012; 何松, 2017)。
值得注意的是,虎林杂岩内云母片岩遭受强烈的变质-变形改造,具有强烈的糜棱面理,而部分侵位于其中的花岗质岩石也具有统一的片理和矿物拉伸线理,这可能代表了晚期岩浆侵位之后二者共同遭受一期变质-变形作用改造。对这一期次变形事件的厘定,对于理解白垩纪早期古太平洋板块的俯冲过程、虎林杂岩的隆升以及虎林盆地沉积过程具有重要的意义。本文正是从虎林杂岩的构造变形入手,细致分析其构造样式,结合变形岩石锆石U-Pb年代学证据,探讨虎林杂岩的隆升时代及变形样式,进而探讨古太平洋板块俯冲过程。
2 宏-微观变形样式
研究区位于黑龙江省虎林市西北(图2a; 采样位置:N45°46′46.4″、E132°56′17.3″),岩体呈北东向展布,主体为黑龙江群,整体被后期早白垩世花岗质岩石侵入和新生代辉绿岩脉不整合覆盖(黑龙江省地质矿产局, 1993; Wildeetal., 2010; Zhouetal., 2010; 胥嘉等, 2012; 何松, 2017)。变质杂岩明显塑性变形,片麻理倾向NW(339°~352°),走向NE-SW,倾角较缓(倾角10°~25°),片理面上普遍发育矿物拉伸线理,主要由暗色黑云母和石英颗粒相间排列组成(图3a),主体表现为NE向左行走滑特征(图2c、图3)。主体岩石为云母片岩,主要组成矿物为石英、黑云母、白云母、较少量长石和副矿物(图3c-g, i)。侵位其中的闪长玢岩脉宽约1m,强烈变形,不规则侵位界线经变形改造,近似平行于片麻理,野外露头可清晰识别出钠长石残斑与暗色角闪石定向排列(图3b),在显微镜下暗色角闪石和钠长石斑晶显示类似的定向排列特征(图3h),钠长石双晶与纤维状角闪石均近平行于片麻理(图3h),基质由长石、石英组成,并且总体呈定向排列,组成弱S-C组构,共同指示左行剪切特征(图3h),表明左行走滑事件发生在闪长玢岩侵位之后。本次工作选取其中的强变形云母片岩和侵位其中的变形闪长玢岩作为研究对象,采样位置见(图2b)。
云母片岩或云母石英片岩均具有片麻状构造,由较为粗粒的石英和黑/白云母相间定向排列而成,各种矿物颗粒或集合体不同程度拉长定向(图3c, d, f-i)。细粒石英塑性拉长成条带状而形成糜棱线理,石英颗粒细粒化,波状消光,颗粒边界产生膨凸式动态重结晶(图3c, d, f-i)。片状黑云母定向排列,大量云母平行定向排列形成叶理面,呈丝带状,由于糜棱岩化而被韧性拉长,部分扭折形成云母鱼和S-C组构,指示左行韧性走滑特征(图3c, d, f)。弱变形闪长玢岩中,钠长石和角闪石斑晶弱定向排列,钠长石颗粒边界相对模糊,可能是受韧性变形改造的结果,颗粒内部脆性破裂明显,与基质内定向排列的细粒长石、石英组成的片麻理面平行,局部角闪石被扭折拉长呈“S-C”状,共同指示左行走滑特征(图3e)。这些显微组构和矿物组合显示岩石变形发生在浅部地壳的绿片岩相环境下(Stippetal., 2002, 2010; 向必伟等, 2007; 胡玲等, 2009; Liangetal., 2015, 2019b; 表1、图3),共同指示了约350~450℃的变形温度。这些中低温组构特征与该地区所表现的走滑相匹配。
3 岩石EBSD组构特征
EBSD(电子背散射衍射)技术是通过分析晶体背散射衍射图像来确定晶轴方向,进而确定晶体颗粒排列的取向性,确定晶体内发育的活动滑移系,进而可以估算矿物的变形温度和滑移系(刘俊来等, 2008; 许志琴等, 2009; 夏浩然和刘俊来, 2011; Liangetal., 2015)。本次EBSD石英C轴组构分析在吉林大学东北亚矿产资源评价自然资源部重点实验室完成,具体流程及判别方法参见Liangetal. (2019b)。为保证测试精度,测试颗粒多选颗粒细小且糜棱质发育部位,采用机器自动识别标定,做图采用下半球投影,极点等密度线以所占百分数表示,测试结果见图4。
图4 典型变形岩石石英C轴组构图
虎林杂岩的石英
表1典型变形岩石显微变形特征
Table 1 Microscopic deformation characteristics of typical deformed rocks
样品号岩性结构构造主要矿物组合微观变形特征18JX1-1云母片岩变晶结构、片状构造Qtz+Bt+Ms+Fel矿物韧性变形,细粒石英塑性拉长成条带状且细粒化,具不连续波状消光,局部石英颗粒产生膨凸式动态重结晶;云母片状定向排列,扭折变形呈云母鱼,指示左行走滑18JX1-2云母片岩变晶结构、片状构造Bt+Qtz+Ms+Fel矿物韧性变形,长石脆性或韧-脆性变形,呈棒状;细粒石英塑性拉长成条带状,部分细粒化,波状消光,部分石英颗粒产生膨凸式动态重结晶.片状黑/白云母定向排列,部分扭折呈S-C组构特征,指示左行韧性走滑特征18JX1-3变形闪长玢岩变晶/斑状结构、弱片麻状构造斑晶:Ab+Am;基质:Qtz+Fel钠长石斑晶表现为脆性破裂,边部模糊细粒化,弱定向排列;角闪石蚀变严重,但整体定向排列,组成弱S-C组成,指示左行走滑特征;基质为细粒长石和石英组成,粒度较小18JX1-4云母片岩变晶结构、片状构造Qtz+Bt+Fel矿物韧性变形,石英细粒化,条带状,波状消光,部分颗粒膨凸式动态重结晶;黑云母定向排列,与石英颗粒组成片麻状构造,指示左行走滑18JX1-5云母石英片岩变晶结构、片状构造Qtz+Ms+Bt+Fel矿物韧性变形,石英颗粒塑性拉长,动态重结晶,波状消光,细粒化明显,产生膨凸式动态重结晶,颗粒边界呈现为镶嵌式接触或锯齿状;片状白云母定向排列18JX1-6云母石英片岩变晶结构、片状构造Qtz+Ms+Fel石英颗粒塑性拉长,波状消光,细粒化,产生膨凸式动态重结晶;白云母丝带状定向排列18JX1-7云母石英片岩变晶结构、片状构造Qtz+Ms+Fel石英颗粒塑性拉长,波状消光,细粒化,产生膨凸式动态重结晶;白云母丝带状定向排列
注:Ab-钠长石;Am-角闪石;Qtz-石英;Bt-黑云母;Ms-白云母;Fel-长石
表2变形闪长玢岩(样品18JX1-3)锆石U-Pb同位素数据
Table 2 Zircon U-Pb isotopic data of deformed diorite porphyrite (Sample 18JX1-3)
测点号Th/U同位素比值年龄(Ma)207Pb206Pb1σ207Pb235U1σ206Pb238U1σ208Pb232Th1σ207Pb206Pb1σ207Pb235U1σ206Pb238U1σ208Pb232Th1σ-010.040.04880.00300.10850.00660.01610.00050.00500.000213788105610331014-020.060.04900.00280.11480.00650.01700.00050.00520.000214581110610931054-030.050.16660.004111.00280.29680.47930.01190.12480.0036252320252325252452237765-040.120.04950.00880.16190.02790.02370.00120.00750.000617327215224151815113-050.020.04860.00150.10240.00330.01530.00040.00640.0002128369939821285-060.130.04830.00290.11330.00680.01700.00050.00600.000211686109610931215-070.140.04840.00410.10880.00900.01630.00050.00530.0003118126105810431065-080.150.04900.00410.11100.00910.01640.00050.00550.0003146126107810531115-090.090.04920.00300.15210.00920.02240.00060.00840.000415987144814341688-100.140.04910.00350.11540.00800.01710.00050.00540.0003153105111710931095-110.060.04930.00350.15800.01090.02330.00070.01000.000616110314910148420112-120.210.04860.00630.12450.01580.01860.00080.00610.00041272081191411951228-130.210.05160.00320.28510.01740.04010.00110.01280.00062698925514253725811-140.290.04900.00520.12740.01310.01890.00070.00680.00041451651221212141368-150.310.04850.00390.11180.00870.01670.00050.00550.0003121117108810731105-160.030.04830.00140.11470.00350.01720.00040.00660.000311433110311031335-170.160.05010.00320.19160.01210.02780.00080.00970.0005197931781017751959-180.360.04850.00290.11390.00670.01700.00050.00560.000212384109610931124-190.400.04870.00390.11170.00870.01670.00050.00540.0003131118107810631085-200.250.04860.00310.10970.00700.01640.00050.00550.000313193106610531115-210.370.05160.00450.12650.01080.01780.00060.00600.00042661351211011441207-220.500.04840.00380.11040.00840.01650.00050.00510.0002120115106810631024-230.390.04860.00360.10820.00780.01610.00050.00580.0003131108104710331165-240.280.04880.00360.11210.00800.01670.00050.00520.0003138107108710731045-250.470.04880.00290.11490.00670.01710.00050.00540.000214083110610931084-260.470.04780.00330.10330.00710.01570.00050.00480.0002889910071003974-270.480.04860.00330.10990.00730.01640.00050.00480.00021279610671053974-280.450.04870.00400.10830.00870.01610.00050.00470.000313312210481033955-290.530.04810.00260.11700.00620.01760.00050.00540.000210673112611331084-300.580.04840.00440.10970.00980.01640.00050.00570.0003121138106910531156-310.680.04830.00390.10900.00870.01640.00050.00520.0002112121105810531055-320.570.04870.00360.12050.00860.01790.00050.00560.0003134107116811531125-330.440.04940.00420.11340.00940.01670.00050.00550.0003164129109910731116-340.530.04850.00410.11800.00980.01770.00060.00520.0003121127113911341046-350.480.04910.00360.11420.00830.01690.00050.00550.0003151110110810831116-360.040.05610.00270.56290.02750.07270.00190.05790.00764586245318453121137145-370.740.05190.00870.11280.01860.01580.00050.00500.00022813421091710131004-380.620.04820.00370.11450.00860.01730.00050.00570.0003107112110811031156-390.220.04900.00220.11350.00510.01680.00040.00530.000214559109510831075-400.790.04890.00590.12830.01510.01900.00070.00620.00051451911231412151249-410.780.04870.00500.10820.01090.01610.00060.00530.00031311601041010341076-420.690.04810.00400.11890.00960.01790.00060.00580.0003104122114911541166-430.370.05580.00180.55620.01830.07230.00180.02510.000744535449124501150114
4 锆石U-Pb年代学
4.1 测试样品及方法
为进一步限定变形时代,选取变形的闪长玢岩脉(18JX1-3)进行LA-ICP-MS锆石U-Pb年代学测试。锆石分选在河北廊坊地质调查院完成,选取5~10kg不等的样品经过破碎、浮选和电磁选方法进行分选,之后在显微镜下观察锆石的标型特征,挑选出各个样品中具有不同特征的锆石单矿物。样品靶在北京锆年领航科技有限公司制备,所选锆石约150颗,置于环氧树脂中,待固结后将锆石抛磨掉大约1/2,使锆石内部结构充分暴露,然后进行锆石的光学、阴极发光图像分析,来选取用于测试微区及测年锆石的部位。LA-ICP-MS锆石U-Pb年代学测试在吉林大学东北亚矿产资源评价自然资源部重点实验室完成。激光剥蚀使用德国相干公司(Coherent)COMPExPro型ArF准分子激光器,质谱仪为美国安捷伦公司7500A型四极杆等离子质谱。激光条件为:激光束斑直径32μm,激光能量密度10J/cm2,剥蚀频率8Hz。剥蚀样品前首先采集30s的空白,随后进行30s的样品剥蚀,剥蚀完成后进行2min的样品池冲洗。载气使用高纯度He气,气流量为600mL/min;辅助气为Ar气,气流量为1.15L/min。对于不用同位素的采集时间,204Pb、206Pb、207Pb和208Pb为20ms,232Th、238U为15ms,49Ti为20ms,其余元素为6ms。使用标准锆石91500(1062Ma)作为外标进行同位素比值校正,标准锆石PLE/GJ-1/Qing Hu为监控盲样。元素含量以国际标样NIST610为外标,Si为内标元素进行计算,NIST612和NIST614为监控盲样。使用Glitter软件进行同位素比值及元素含量的计算。谐和年龄计算及图像绘制采用国际标准程序Isoplot(ver 3.23; Ludwig, 2003)。普通铅校正使用Andersen (2002)给出的程序计算。分析数据及锆石U-Pb谐和图给出误差为1σ,表示95%的置信度(表2)。
4.2 测试结果
对变形闪长玢岩(18JX1-3)进行锆石U-Pb测年,以期限定其变形年龄。锆石CL图像见图5a,定年结果见图5b, c、表2。CL图像显示选取的锆石大多结晶较好,自形到半自形,晶型较完整,多为灰白色长柱状晶形,颗粒长宽比2:1~3:1,粒径多分布在80~150μm之间。部分锆石晶棱较圆滑,呈深灰色,振荡环带变较弱,局部由核到边发生模糊、退化的现象,且核部宽颜色发暗、边部窄而亮,呈现浑圆粒状、椭圆粒状及粒状,表现出继承性锆石的特点(图5a、表2)。但多数锆石内部具有清晰的振荡环带结构,锆石Th/U比值为0.02~0.79,均大于0.01,显示为岩浆成因(图5a、表2)。说明该样品中锆石的主体是岩浆锆石,但存在早期受不同程度变质作用改造的捕获锆石,部分捕获锆石具有明显核边结构(测点13,43;图5a)。总体而言,主体锆石显示岩浆锆石特点,此次共测试43颗锆石,可识别出6期峰值(图5b, c),其中8个测点年龄较老,~146.5Ma(n=3)、~175.5Ma(n=1)、~255.5Ma(n=1)、~453.5Ma(n=2)和~2520Ma(n=1),CL图像显示这些较老年龄锆石以灰白色短柱状为主,长度100~150μm,长宽比为1:1~2:1,锆石颗粒内部呈现无分带、弱分带、面状分带和溶蚀结构等(图5a),揭示这些锆石可能为继承锆石、捕获锆石或变质锆石,所代表的年龄可能分别对应于寄主岩石的源区残留锆石、捕获锆石及变质事件(一期和/或多期)的年龄,当然,也不排除一些测点年龄值为混合年龄的可能。35颗锆石年龄显示(图5c),206Pb/238U加权平均年龄为107.0±1.7Ma(MSWD=2.7,n=35),结合锆石的同位素特征及内部结构,确定变形闪长玢岩的岩浆结晶年龄或岩脉就位年龄为107.0±1.7Ma,形成于早白垩世晚期,因此,可以限定晚期NE向左行走滑事件晚于闪长玢岩就位时间。
5 讨论
5.1 虎林杂岩伸展变形样式与变形机制
虎林杂岩塑性变形特征明显,片麻理倾向NW,具有统一的NE向低角度矿物拉伸线理,云母鱼以及S-C组构均指示一期左行走滑事件。细粒石英波状消光和不同程度细粒化,指示石英变形以晶内塑性变形、位错滑移和膨凸式动态重结晶为主要的变形机制。云母被韧性拉长,部分扭折形成云母鱼、S-C组构,也指示一种中低温变形样式。变形闪长玢岩中钠长石斑晶颗粒边界模糊,可能指示长石由脆性破裂向韧性变形的转变。结合石英EBSD组构图中可识别出两期变形的叠加,早期稍高温度变形,石英以菱面滑移系为主,晚期底面滑移系启动,但结合矿物组合和矿物微观变形样式,岩石总体变形为中低温变形,变形温度不超过高绿片岩相(350~450℃),石英变形机制以位错滑移和膨凸重结晶为主。结合宏观变形样式可以得知,早期变形以NW倾向片麻理为标志,指示NW-SE向伸展作用,晚期变形以NE走向低角度矿物拉伸线理为标志,指示NE-SW向左行走滑事件。
图5 典型锆石CL图像(a)与锆石U-Pb谐和图(b、c)
5.2 兴凯地块早白垩世岩浆-构造事件与虎林杂岩隆升关系
兴凯地块主体位于俄罗斯境内,少量出露于我国境内,是兴蒙造山带的重要组成部分(李旭平等, 2010; Xuetal., 2013a; 徐备等, 2014; Liuetal., 2017a; Yangetal., 2017a, b; 周建波等, 2018; 许文良等, 2019; 刘永江等, 2019; 肖文交等, 2019)。兴凯地块广泛出露古生代-中生代花岗质岩体,为进一步限定虎林杂岩乃至整个兴凯地块的岩浆序列,将兴凯地块已有锆石年龄数据进行统计(图6),主体年龄可划分为前寒武纪、早古生代、晚古生代晚期和白垩纪。
图6 兴凯地块花岗质岩石年龄分布图
少量前寒武纪结晶基底的年龄,如~2500Ma、~1500Ma、~700Ma的变质年龄,来源于基底岩石锆石的捕获,记录了各微板块的基底年龄,乃至从Rodinia到Gondwana大陆的聚合与离散记录(Zhaoetal., 2000; Zhou and Wilde, 2013; 孙立新等, 2013; 吕长禄等, 2014; 周建波等, 2009, 2011, 2016; Sorokinetal., 2016; Zhouetal., 2011, 2018; 刘永江等, 2019)。大量430~550Ma的锆石年龄,记录了兴凯地块早古生代岩浆和变质事件(颉颃强等, 2008a, b; 温泉波等, 2008; Wildeetal., 2010; Zhouetal., 2010; Khanchuk, 2010),与佳木斯地块麻山群的年龄一致,表明二者具有相同的泛非期变质基底(Wildeetal., 2010; Zhouetal., 2010; 周建波等, 2009, 2011; Yangetal., 2014, 2017a, b; Bietal., 2015)。
晚古生代晚期~256.2Ma的峰值年龄与Zhouetal. (2010)在虎林云母片岩内获取的257±3Ma结晶年龄一致,可能与佳木斯-兴凯地块裂解事件相关(Zhouetal., 2010; Khanchuketal., 2010),之后早中侏罗世(~200Ma和170~180Ma)的捕获锆石年龄可能与佳木斯-兴凯地块与松嫩地块拼合以及拼合后区域伸展环境导致的岩浆活动有关(Wuetal., 2007; 黄映聪等, 2008; Zhouetal., 2010; Wuetal., 2011; 张磊等, 2013; Wangetal., 2015; 丛智超等, 2016; Zhuetal., 2015, 2017; 曾振等, 2017a, b; 周建波等, 2018),当然也不排除与古亚洲洋在晚二叠-早三叠世闭合后的伸展造山环境有关,或佳木斯地块东侧古洋板块的西向俯冲作用有关(毕君辉等, 2014; 于介江等, 2013)。三叠纪末-侏罗纪初,古亚洲洋构造域结束演化并进入古太平洋构造域演化阶段(Zhouetal., 2010; Khanchuketal., 2010; Wuetal., 2011; Guoetal., 2015; Liuetal., 2017a, b, c)。
兴凯地块虎林杂岩内早白垩世侵入体也明显显示3个峰期(图6),~104Ma、~115Ma和~128Ma(图6),但总体较为连续,也可以认为是一期岩浆事件的不同演化阶段。早白垩世晚期(128~104Ma)期间,佳木斯-兴凯地块已与松嫩地块拼合完成,完达山-锡霍特-阿林造山带主体也已拼贴完成,中国东北东部整体进入古太平洋板块持续俯冲阶段(Wildeetal., 2010; Wuetal., 2007, 2011; 张兴洲等, 2011; 胥嘉等, 2012; 许文良等, 2013; Xuetal., 2013b; Guoetal., 2015; Liuetal., 2017a, b; 周建波等, 2010, 2018),该期岩浆事件主体形成于大陆弧后伸展环境,而相应的广泛分布于东北东部乃至俄罗斯远东的早白垩晚期钙碱性火山岩也多形成于活动大陆边缘的构造背景,与大兴安岭和松辽盆地内部双峰式火山岩一致,均与陆内伸展有关(Dengetal., 2019; Xuetal., 2013b; 孙明道, 2016),是古太平洋板块俯冲后撤相关产物(胥嘉等, 2012; 何松, 2017; Kirillova, 2005; 程瑞玉等, 2006; 孙景贵等, 2008, 2009; Wuetal., 2002, 2011; Xuetal., 2013b; Ouyangetal., 2013; Yangetal., 2014; Sunetal., 2015a, b)。
图7 中国东北东部早白垩世构造演化图(据朱光等, 2018修改)
虎林杂岩紧邻敦密断裂(图1),前人对敦密断裂研究限定其在早白垩世100~130Ma期间为强烈伸展期,伴随着大规模的韧性变形和火山喷发作用(孙晓猛等, 2016),同时,虽然可以将早白垩世岩浆划分为三期,但很明显的是~115Ma岩浆活动最为强烈,之后趋于平静,仅有的地球化学分析表明,早白垩世花岗质岩石为I型花岗岩,形成于大陆弧引张环境,源区为壳幔混合源,主要来自加厚下地壳基性岩石部分熔融,并受俯冲带流体交代作用形成的富集地幔部分熔融作用的影响(胥嘉等, 2012; 何松, 2017),与中国东北东部同时期的伸展环境相一致。
虎林杂岩主体由原属中元古界“黑龙江群”(黑龙江省地质矿产局, 1993)的云母/石英片岩组成,Zhouetal. (2010)在云母片岩中获得233~270Ma的年龄,指示其原岩形成时代为233Ma(中三叠统)之后,并非前人认为的中元古界。云母片岩被变形的闪长玢岩脉(107.0±1.7Ma,本文)所侵入,指示主期左行走滑事件晚于岩脉侵位时间。详细的野外观察发现,虎林杂岩强烈韧性变形,早期变形以NW倾向片麻理为标志,指示NW-SE向伸展作用,该期伸展作用与区域岩浆事件相吻合。晚期变形以NE向低角度矿物拉伸线理为标志,指示NE-SW向左行走滑事件。闪长玢岩脉在侵位过程中,伴随着虎林杂岩整体隆升,之后在虎林片岩继续向地壳浅部抬升过程中,在挤压背景下,形成NE-SW向左行走滑事件,产生低角度矿物拉伸线理。
图8 虎林盆地地层表(a, 据韩国卿等, 2008修改)及地震反射剖面解释图(b, 据刘志宏等, 2014修改; 位置见图1a)
综上所述,早白垩世期间,古太平洋-伊泽奈崎板块由向NWW运动转变为向NW运动,从而导致东北东部乃至中国东部东部伸展方向由WNW-ESE向转变为NW-SE向(Wuetal., 2011; Zhuetal., 2010, 2012; Wangetal., 2018; 朱日祥等, 2012; Kuskyetal., 2014; 周丽云等, 2015; Zhou and Li, 2017; 郑建平和戴宏坤, 2018; 图7a),对应着古太平洋-伊泽奈崎俯冲板片的后撤,在弧后拉张的动力学背景下,诱发岩石圈大规模减薄拆沉,致使软流圈地幔和地壳直接接触,岩浆底垫作用和深部地壳的高温变质及部分熔融作用,导致包括虎林花岗质岩石在内的中国东北东部地区早白垩世(128~104Ma)发生剧烈火山-岩浆活动和突发性成矿作用(毛景文等, 2005; 程瑞玉等, 2006; 孙景贵等, 2008, 2009; Wuetal., 2002, 2011; Xuetal., 2013b),该期伸展活动同时造就了中国东北东部地区大量伸展断陷盆地的形成,伸展方向为NW-SE向(周建波, 2009; 张兴洲和马志红, 2010; 林伟等, 2013; 张兴洲等, 2015; 孙晓猛等, 2016; Wangetal., 2017, 2018; 李三忠等, 2017; Daietal., 2018; Lietal., 2019),这也与中国东部岩石圈减薄高潮期(130~120Ma)的观点基本一致(吴福元等, 2003; 董树文等, 2007, 2008; Wuetal., 2002, 2011; 朱日祥等, 2012; Dongetal., 2018; Wangetal., 2018),一直持续到早白垩世末期。该期伸展事件也导致了虎林杂岩隆升、系列犁式正断层的发育和虎林盆地的大范围坳沉(具体见下节)。
在早白垩世末(104Ma之后),古太平洋-伊泽奈崎板块向NNW向运动可能再次加速(Maruyamaetal., 1997; 刘志宏等, 2014; 朱光等, 2018),由NNW转向N,由低角度俯冲替代高角度俯冲,使得中国东部再次遭受区域性挤压,导致虎林杂岩发育大规模NE向左行走滑变形和持续性隆升,此时敦-密断裂、佳-伊断裂乃至整个郯庐断裂带显示为左行汇聚型走滑,同时,岩石圈减薄拆沉告一段落,大规模岩浆事件较少出现(图7b)。同样,该期挤压事件在整个中国东部也可以发现,如整个郯庐断裂带的伸展活动减弱,盆地反转、隆升且范围缩小,岩浆活动明显减弱(朱光等, 2018),是华北克拉通破坏结束的标志之一(Wuetal., 2002, 2011; 朱日祥等, 2012; 朱光等, 2018; Daietal., 2018),也是燕山运动最后一期挤压事件导致的不整合的表现(董树文等, 2007, 2008; 朱光等, 2018; Dongetal., 2018; Wangetal., 2018)。类似的左行走滑构造在锡霍特-阿林地区能够很清晰地识别,早期形成的浊积岩盆地、增生楔和火山岛弧,均沿着东亚大陆边缘向北东移动,如佳木斯-兴凯地块在这一时期(Albian期末期)的运移、拼贴,同时导致侏罗纪形成的与古太平洋俯冲相关的褶皱构造再一次遭受挤压,轴面陡倾变为紧闭褶皱(Kemkinetal., 2016; 刘凯等, 2016),在虎林地区识别出的大量紧闭褶皱,可能与该期构造有关。而晚白垩世-古新世期间,古太平洋-伊泽奈崎板块又转为斜向俯冲,区域性左行走滑作用减弱,域内岩浆活动频发,形成东锡霍特-阿林陆源火山弧(杂岩带)和西库页前陆盆地等地质单元(Dengetal., 2019; Xuetal., 2013b; 邵济安和唐克东, 2015; 孙明道, 2016; Kemkinetal., 2016; Wangetal., 2019)。
5.3 虎林杂岩隆升与虎林盆地沉积的耦合关系
虎林盆地是佳木斯隆起周缘中新生代盆地之一,与敦-密断裂带毗邻,总体上表现为北东向隆坳相间分布,其内的北部坳陷区和南部坳陷区被NEE向展布的中央隆起所分隔(曹成润等, 2001; 张凤旭等, 2010)。中央隆起区主要由虎林杂岩和少量白垩纪花岗岩侵入体组成,上部几乎没有中-新生代地层覆盖。南北凹陷区主要发育有白垩系、渐新统和中新统,自下而上可见裴德组(K1p)、七虎林组(K1q)、云山组(K1y)、珠山组(K1z)、东山组(K1d)、虎林组(E2-3h)、富锦组(N1f)和第四纪系(图8a)。但重磁异常显示盆地存在明显的高磁高重力异常,说明盆地下部分布着统一的较老的结晶基底(曹成润等, 2001; 刘银萍和孟令顺, 2007; 张凤旭等, 2010)。
在裴德期-东山期,虎林盆地处于强烈伸展构造体制当中,受呈NE-NNE向展布的主要伸展控陷断层影响,形成一系列相互独立的伸展断陷,在断陷中不仅有正常沉积的碎屑岩类,而且伴随广泛分布的早白垩世钙碱性火山碎屑岩建造主要为岛弧系列岩浆岩(王伟涛等, 2007)。早白垩世碎屑岩建造主要为海陆交互相沉积,并见有植物化石和大量海相化石,如菊石、双壳类、腹足类、介形虫和沟鞭藻化石,揭示盆地的形成时代为Hauterivian期(132.9Ma)至Albian期(100.5Ma)(何承全等, 1999)。同时,伴随的广泛分布的早白垩世晚期东山组火山岩,其形成时代为110~117Ma(张兴洲和马志红, 2010),二者共同指示虎林盆地早白垩世沉积时代应该为132.9~100.5Ma,与区域岩浆事件相吻合(图6)。值得关注的是,盆地南部穆棱河凹陷内的白垩系以变质岩岩屑为主,与中央隆起出露的基底变质岩组成一致(张兴洲和马志红, 2010),暗示盆地基底岩石是其物源之一。在虎林盆地北缘跨过敦-密断裂的地震反射剖面中可以发现,在控陷断层上盘,离断裂较近的部位,下白垩统的沉积厚度更大一些,远离断裂的部位,下白垩统的沉积厚度逐渐减小,反映出虎林盆地在早白垩世受到NW-SE向的持续伸展作用(刘志宏等, 2014; 图8b),可能与古太平洋-伊泽纳奇板块NW向俯冲和板片后撤有关,形成了系列的NE向铲式正断层控制了普遍存在的NE-NEE向伸展断陷盆地群的形成、广泛的沉积作用和盆地间的断垒构造(图7a、图8b; 刘志宏等, 2014),虎林盆地中央隆起和周缘地层区为其主要的沉积物源(王伟涛等, 2007; 韩国卿等, 2008; 张兴洲和马志红, 2010)。
早白垩晚期区域持续隆升机制与古太平洋-伊泽奈崎板块高速向欧亚大陆俯冲由NNW转向N有关,造成区域性的挤压隆升(图7b; 韩国卿等, 2008; 朱光等, 2018),敦-密断裂大规模的左行走滑,同时,虎林盆地由早白垩世的伸展应力体制转换为挤压应力体制(图7b),出现系列的构造反转,使早期形成的伸展控陷正断层转换为NE向压扭性的左行走滑(刘志宏等, 2014),伴随虎林杂岩的强烈韧性变形,形成NE向近水平矿物拉伸线理,同时导致早白垩世形成的伸展型虎林盆地大幅度快速隆升。该时期的挤压隆升导致盆地在东山组沉积期后整体处于剥蚀状态,而且隆升剥蚀速率大于沉积速率,持续为周缘盆地提供物源,使得盆内晚白垩世及古近纪初期地层保留甚少(王伟涛等, 2007; 韩国卿等, 2008; 刘志宏等, 2014; 图8a),直至始新世开始接受沉积,沉积以粗碎屑岩建造为主的虎林组(韩国卿等, 2008; 图8a),盆地转换为挤压拗陷型盆地。
6 结论
(1)虎林杂岩塑性变形特征明显,发育片麻理倾向NW,低角度矿物拉伸线理走向NE-SW。微观构造特征和石英EBSD组构分析显示,岩石总体变形为中低温变形,变形温度350~450℃,石英变形机制以位错滑移和膨凸重结晶为主。
(2)结合宏观变形样式可以得知,虎林杂岩早期变形以NW倾向片麻理为标志,指示NW-SE向伸展作用,晚期变形以NE向低角度矿物拉伸线理为标志,指示NE-SW向左行走滑事件。变形闪长玢岩的锆石206Pb/238U加权平均年龄为107.0±1.7Ma(MSWD=2.7, n=35),指示岩脉形成于早白垩世晚期,限定NE向左行走滑事件晚于闪长玢岩就位时间。
(3)早白垩世期间,伊泽奈崎板块NW向俯冲和板片后撤,导致东北东部乃至中国东部东部NW-SE向伸展,造成大陆岩石圈减薄、岩浆活动和包括虎林盆地在内的伸展断陷盆地的形成。在早白垩世末(104Ma之后),伊泽奈崎板块高速向欧亚大陆俯冲由NNW转向N,由低角度俯冲替代高角度俯冲,使得中国东部再次遭受区域性挤压,导致虎林杂岩发育大规模NE向压扭性左行走滑变形和盆地的持续隆升。
致谢对吉林大学东北亚矿产资源评价自然资源部重点实验室在锆石U-Pb测试和EBSD组构测试方面给予的支持表示感谢。感谢吉林大学郑常青教授和周建波教授在野外工作中的帮助与指导,感谢刘志宏教授在论文修改过程中提供的帮助。同时感谢两位审稿人提出的宝贵而细致的修改意见。