甘保平 第五春荣王伯隆 王庭怡
西北大学大陆动力学国家重点实验室,西北大学地质学系,西安 710069
华北克拉通是世界上最古老的克拉通块体之一,出露了大量的前寒武纪的岩石,且保留了3.85Ga的大陆地壳岩石(Liuetal., 1992; Songetal., 1996;Wanetal., 2005),这些前寒武纪基底岩石在早前寒武纪时期(太古宙-古元古代)经历了广泛的变质变形-岩浆构造热事件(Zhaoetal., 2001, 2005,2006,2010;Santosh, 2010;耿元生等,2010;翟明国,2011)。因此,华北克拉通已成为研究早前寒武纪大陆地壳增生和演化的关键地区之一。按照Zhaoetal. (2001, 2005)的划分方案,华北克拉通主要由西部陆块、东部陆块和二者之间近南北向的中部造山带组成(图1a),并将西部陆块进一步划分为北部的阴山地块,南部的鄂尔多斯地块及夹在二者之间近东西向的孔兹岩带(图1b)。其中,孔兹岩带主要由古元古代孔兹岩系、少量的TTG片麻岩、紫苏花岗岩和S型花岗岩组成,被认为是由鄂尔多斯地块和阴山地块大约在~1.95Ga碰撞-拼贴-聚合形成的(Zhaoetal., 2005, 2010, 2012)。
图1 华北克拉通基底构造单元划分图(a, 据Zhao et al., 2005)、西部陆块构造位置略图(b,据Zhao et al., 2005)及贺兰山北段宗别立地区地质简图(c,据周喜文等,2010;Yin et al., 2011修改)Fig.1 Basement tectonic subdivisions of North China Craton (a, modified after Zhao et al., 2005); tectonic subdivision of the Western Block (b, modified after Zhao et al., 2005)and geologic map of the Zongbieli area in the northern Helanshan (c, modified after Zhou et al., 2010;Yin et al., 2011)
近年来,诸多学者对华北克拉通西部孔兹岩带中的孔兹岩系岩石开展了研究工作,取得了大量的研究进展,认为贺兰山地区孔兹岩的原岩沉积时代可能在2.0~1.95Ga,且经历了大约两次变质事件(~1.95Ga和~1.85Ga),而东部乌拉山-千里山孔兹岩的变质年龄峰值也主要集中在这两个时期(Zhaoetal., 2005, 2006, 2010; Yinetal., 2009, 2011, 2015; Danetal., 2012; Guoetal., 2012; Dongetal., 2014; Pengetal., 2014; 蔡佳等,2015)。在华北克拉通西部贺兰山-千里山地区也有大量年龄为2.0~1.85Ga的古元古代S型花岗岩出露(耿元生等,2009;Yinetal., 2011; Danetal., 2012),其中一些研究者对孔兹岩带西部贺兰山地区~1.95Ga的S型花岗岩研究表明其形成与阴山地块和鄂尔多斯地块碰撞-拼合有关(Danetal., 2012; 刘金科等,2016;Lietal., 2017;Zhangetal., 2017),认为是孔兹岩系部分熔融的产物(Yinetal., 2011;张道涵,2017;余淳梅等,2019)。而年龄为~2.0Ga的S型花岗岩在华北克拉通西部贺兰山地区亦有出露,主要分布在贺兰山北段(宗别立地区)(耿元生等,2009;Lietal., 2017)和中段(黄旗口、赵池沟地区)(杨华本,2013;Lietal., 2018;余淳梅等,2019)。相对而言,对贺兰山地区出露古元古代中期(~2.0Ga)的花岗质岩石研究相对较少,并且S型花岗质岩石作为华北克拉通早前寒武纪地质体中最主要的岩石之一,对其进行系统的岩石学,年代学和地球化学分析,可为深入探讨华北克拉通西部孔兹岩带古元古代的构造演化提供更多的依据。
本文以华北克拉通西部贺兰山地区的古元古代S型花岗岩为研究对象,通过详细的岩石学、岩石地球化学、锆石U-Pb年代学和Hf同位素研究,厘定S型花岗岩的侵位年龄,进而探讨其岩石成因机制及构造意义。结合已有的区域地质构造热事件,约束其华北克拉通西部阴山地块和鄂尔多斯地块俯冲碰撞的时限,为理解孔兹岩带早前寒武纪构造演化的过程提供重要的地质依据。
孔兹岩带位于华北克拉通西部陆块的中部,其北侧为阴山地块,南侧为鄂尔多斯地块,呈近东西向线状展布,宽约200km(图1a),主要由东部地区的集宁-凉城杂岩、中部地区的大青山和乌拉山杂岩以及西部地区千里山和贺兰山变质杂岩组成(图1b)。孔兹岩带的岩石主要由一套含石榴子石石英岩、含石墨富铝片麻岩、长英质副片麻岩、大理岩、钙硅酸盐岩和少量的斜长角闪岩组成。传统的观点认为这些岩石的原岩通常形成于稳定大陆边缘的构造环境(卢良兆等,1996),其变质程度多达高角闪岩相,局部可达到高压麻粒岩相(Zhouetal., 2010; Yinetal., 2014,2015)。
研究区位于贺兰山北段宗别立地区,区内孔兹岩系与变形花岗岩呈渐变或复杂的构造接触,所以把这些变质表壳岩和花岗岩类统称为贺兰山变质杂岩。区内的变质表壳岩主要以贺兰山岩群(富铝岩系)为主,出露面积约1000km2,是贺兰山杂岩的主要组成之一。贺兰山岩群主要由富铝片麻岩、麻粒岩、变粒岩、夹大理岩和钙硅酸盐岩等变质岩组成,其中变粒岩类和片麻岩是主体,约占变质地层的90%以上(胡能高和杨家喜,1995;卢良兆等,1996)。根据野外岩石组合特征,将其划分为上部的秃鲁根变粒岩-大理岩岩组,中部的阿楞呼都格变粒岩岩组及下部的柳树根富铝片麻岩岩组(图1c)(胡能高,1994)。这套杂岩普遍经历了中高级变质作用,特征变质矿物主要有石榴子石、夕线石和堇青石等,局部地区可见紫苏辉石。其中,中部和下部岩层中的片麻岩整体变质程度已达到了麻粒岩相(卢良兆等,1996)。此外,贺兰山岩群中还出露一些古元古代辉绿岩和辉绿玢岩,走向以近东西向和北北东向为主,风化面呈深灰色-深灰绿色,主要由普通辉石和斜长石组成,出露宽约4~15m不等,局部有少量的基性麻粒岩透镜体出露(胡能高,1994)。
研究区内花岗质岩石出露较为广泛,岩石类型主要有含石榴黑云母花岗岩、石榴石花岗岩和似斑状黑云母花岗岩等,其中以石榴石花岗岩出露最为广泛。此外,在宗别立地区北部有变形花岗岩出露,形成年龄为2053±58Ma (耿元生等,2009)。本次研究样品主要采自于贺兰山北段宗别立地区(图1c),样品采集GPS点位见后文。
片麻状含石榴黑云母花岗岩出露面积较小,以蒙果特东北部分布居多,岩体东西向分布,岩石风化严重,野外与围岩接触关系不明显。新鲜面呈浅灰色-灰色,片麻状构造,斑状变晶结构(图2c)。主要矿物为斜长石(~25%),钾长石(~30%),石英(~35%),黑云母(~5%)和石榴子石(<5%)。其中斜长石多发生绢云母化,石英呈他形粒状,石榴子石呈大斑晶,粒径多数为2~5mm。副矿物主要有锆石、磷灰石和磁铁矿等(图2e, f)。
似斑状石榴黑云母花岗岩主要分布在巴音恩德勒,扣笨和浩尧尔地区,出露规模较大,局部与片麻岩呈侵入接触关系(图2a)。新鲜面呈浅肉红色,块状构造,花岗变晶结构(图2b)。主要矿物组成为石英(~25%),钾长石(~20%),斜长石(~25%),微斜长石(~15%),黑云母(~8%),石榴子石(~7%)。副矿物主要由锆石,榍石和磷灰石等。其中斑晶为钾长石,多数呈半自形-他形,大多数斑晶粒度在0.2×0.6cm左右,可见简单双晶;斜长石呈板状,半自形-他形,发育聚片双晶;微斜长石具有明显的格子双晶;石榴子石颗粒呈大斑晶产出,一般粒径约0.5~1cm,最大粒径可达4~5cm左右(图2d, g, h)。
锆石单矿物分选在河北省区域地质矿产调查研究所实验室完成。首先采用常规重力和磁选方法从岩石样品中分选出锆石颗粒,最后在双目镜下手工挑纯。本次研究中的锆石制靶、CL图像拍摄、锆石U-Pb、Lu-Hf同位素分析、主量和微量元素测定均在西北大学大陆动力学国家重点实验室完成。
首先将双目镜下挑选出来的代表性锆石颗粒粘贴于环氧树脂中,其次将其固化抛光,再利用超声波在纯净水中将其表面的污染物除尽,洗净,最后进行透、反射光及CL图像的拍摄。锆石U-Pb同位素组成分析在激光剥蚀电感耦合等离子体质谱(LA-ICP-MS)仪上完成。激光剥蚀系统为193nm ArF-excimer激光器(Geolas 2005),ICP-MS为安捷伦公司的Agilent 7500a ICP-MS,分析时采用激光剥蚀斑束为32μm,激光脉冲为6Hz,能量为10mJ/cm2。年龄计算以标准锆石91500为外标进行同位素比值进行分馏校正,而元素浓度计算采用NIST610为外标,29Si为内标,具体方法可见Liuetal. (2007)。
锆石微区原位Lu-Hf同位素测定在英国Nu公司生产的型号为Nu Plasma HR多接收电感耦合等离子体质谱仪(MC-ICP-MS)上完成,激光能量密度为6mJ/cm2, 频率为5Hz,斑束为43μm,载气为高纯氦气,流速为280mL/min。Lu-Hf同位素分馏校正采用指数法则计算,采用176Lu/175Lu=0.02656和176Yb/173Yb=0.78696比值扣除176Lu和176Yb对176Hf的干扰,最终获得准确的176Hf信号值。详细的分析方法和仪器参数详见Yuanetal. (2008)。
首先选择新鲜的、均匀的代表性样品,在小型颚式破碎机击上将其碎成粒度约5~10mm的颗粒,再将其用玛瑙研钵托盘在振动式碎样机中碎至约200目,将细碎后的粉末样品用来做化学成分分析测试。
全岩主量元素分析在XRF(Rugaku RIX 2100)仪器上测定。 微量元素分析在Agilent 7500a ICP-MS上完成,样品测试中用国际标样BCR-2、BHVO-1和AGV-1作为外标进行监控及校正。所有元素的分析精度和准确度优于5%,详细的分析测试流程见刘晔等(2007)。
图2 贺兰山北段宗别立地区花岗岩的野外照片(a-d)及显微镜下照片(e-h) (a)花岗岩与围岩(含榴片麻岩)的野外接触关系;(b)似斑状石榴黑云母花岗岩;(c)片麻状含石榴黑云母花岗岩;(d)似斑状石榴黑云母花岗岩中的石榴子石; (e、f)片麻状含石榴黑云母花岗岩(其中e为单偏光照片); (g、h)似斑状石榴黑云母花岗岩. 矿物缩写据沈其韩(2009):Kfs-钾长石;Pl-斜长石;Mc-微斜长石;Bt-黑云母;Qtz-石英;Grt-石榴子石Fig.2 Field (a-d) and microscopic (e-h) photos of the granites from the Zongbieli area in the northern Helanshan (a) field photograph of the contact relationship between granite and wall rock (garnet-bearing gneiss); (b) porphyritic garnet biotite granite; (c) gneissic bearing-garnet biotite granite; (d) garnet photographs in porphyritic garnet biotite granite; (e, f) microscopic photos of gneissic bearing-garnet biotite granite (e under plane-polarized light); (g, h) microscopic photos of porphyritic garnet biotite granite. Mineral abbreviation after Shen (2009): Kfs-K-feldspar; Pl-plagioclase; Mc-picrocline; Bi-biotite; Qtz-quartz; Grt-garnet
对贺兰山杂岩中的片麻状含石榴黑云母花岗岩(08HL93和08HL100)和似斑状石榴黑云母花岗岩(08HL95)进行锆石U-Pb年龄测试分析,其结果列于表1,锆石CL图像如图3所示。
片麻状含石榴黑云母花岗岩(08HL93)中锆石大部分为无色透明,多数呈半自形-自形长柱状,锆石粒度约50~150μm,长宽比约为2:1~3:1,锆石的CL图像中显示其阴极发光性相对较差,韵律环带不明显(图3a)。对16粒锆石进行测试分析,其中8、9和15号测点的207Pb/206Pb年龄值约为~2.06Ga, 通过观察CL图像,这3粒锆石可能为继承锆石。其余13颗锆石的Th、U含量分别为10.7×10-6~375×10-6、171×10-6~983×10-6,Th/U比值除了13和14号测点小于0.1外,其余测点均大于0.1(0.11~0.84),表明这些锆石为岩浆成因(Belousovaetal., 2002; Williamsetal., 2009)。这13粒锆石的测试结果构成了一条线性较好的不一致线,表明这些锆石在后期的地质演化过程中经历了不同程度Pb丢失,其207Pb/206Pb上交点年龄值为2032±39Ma (MSWD=2.7,n=13;图3b),该年龄应解释为此片麻状含石榴黑云母花岗岩的侵位年龄。
表1 贺兰山北段古元古代花岗岩LA-ICP-MS 锆石U-Pb定年分析结果
Table 1 LA-ICP-MS zircon U-Pb data of the Paleoproterozoic granites in the northern Helanshan
测点号ThU(×10-6)Th/U 同位素比值年龄(Ma)207Pb206Pb1σ207Pb235U1σ206Pb238U1σ207Pb206Pb1σ207Pb235U1σ206Pb238U1σ08HL93 片麻状含石榴黑云母花岗岩-012515510.460.119280.003084.308890.060460.262070.00340194645169512150017 -0266.63420.190.120920.003165.169420.075610.310150.00406197046184812174220 -031079830.110.101090.002592.306960.031550.165550.0021316444712141098812 -0465.24680.140.122130.003174.842120.069530.287590.00374198846179212163019 -0595.95450.180.114480.002933.762080.051730.238370.00307187246158511137816 -061563570.440.116400.002993.588560.049670.223620.00288190245154711130115 -071815120.350.116070.002954.253890.057050.265810.00341189745168511152017 -0818.735.60.530.127730.003666.372700.119550.361840.00509206750202916199124 -0935.463.70.560.127220.003505.695570.096420.324680.00441206048193115181321 -101742630.660.115920.003094.298810.066890.268930.00355189447169313153518 -113755060.740.117200.003004.177850.057260.258510.00332191445167011148217 -121441710.840.125000.003236.030720.085880.349870.00454202945198012193422 -1310.77390.010.113130.002883.975190.053780.254780.00326185045162911146317 -1416.68230.020.105540.002722.489920.035060.171070.00220172447126910101812 -151671770.940.127470.003366.036780.091690.343390.00452206346198113190322 -162797730.360.115370.002933.509900.047690.220540.0028118864515291112851508HL100 片麻状含石榴黑云母花岗岩-012831232.310.123080.003396.050790.112150.356210.00515200148198316196425 -021642760.600.121450.003055.661730.081420.337850.00445197844192612187621 -0365.81560.420.118860.003045.713100.086510.348400.00466193945193313192722 -042443800.640.120840.003084.197140.062090.251790.00332196945167412144817 -0584.41800.470.121640.003125.987470.090160.356890.00477198045197413196723 -061902710.700.111950.002883.045040.045510.197240.00259183146141911116114 -0791.41150.790.120850.003255.714640.096000.342990.00475196947193415190123 -081782920.610.120500.003185.758840.091570.346660.00471196446194014191923 -091512020.750.117350.003235.650340.100290.349270.00493191649192415193124 -1055.61940.290.123050.003365.009910.086620.295360.00412200148182115166820 -112293320.690.124840.003335.441160.087920.316230.00431202747189114177121 -1258.63040.190.115500.003073.812490.060080.239500.00321188847159513138417 -136807180.950.112060.003003.347200.053900.216730.00291183348149213126515 -141072320.460.121420.003135.863460.084400.350400.00461197745195612193722 -153241841.760.123050.003445.595060.101380.329960.00469200149191516183823 -1682.13080.270.121560.003154.929400.071510.294260.0038719794518071216631908HL95 似斑状石榴黑云母花岗岩-011213130.390.122820.003794.078290.065460.240910.00360199854165013139119 -0296.82920.330.125200.003864.502630.072540.260920.00391203254173213149520 -0363.23740.170.122740.003784.189720.066770.247660.00370199654167213142619 -0492.52000.460.117060.003644.494410.074510.278560.00421191255173014158421 -0573.13790.190.127330.003924.767790.076180.271670.00408206153177913154921 -0645.181.30.550.122930.004105.944550.121650.350870.00567199958196818193927 -0778.21690.460.120570.003906.450920.121580.388210.00613196557203917211528 -0875.22450.310.125150.003874.681940.075870.271430.00411203154176414154821 -0961.694.80.650.124620.004066.107720.117470.355600.00566202357199117196127 -1070.95600.130.123460.003913.961480.069220.232810.00359200755162614134919 -1136.62620.140.119570.003785.999490.104970.364050.00564195055197615200127 -121493360.440.123580.003885.442180.092980.319510.00492200955189215178724 -1346.04010.110.119550.003764.937090.084450.299630.00461195055180914169023 -1432.72650.120.119990.003755.062880.085240.306150.00470195655183014172223 -1543.82230.200.120720.003866.090010.110670.366020.00575196756198916201127
续表1
Continued Table 1
测点号ThU(×10-6)Th/U 同位素比值年龄(Ma)207Pb206Pb1σ207Pb235U1σ206Pb238U1σ207Pb206Pb1σ207Pb235U1σ206Pb238U1σ-1672.03330.220.123070.003855.229590.087900.308300.00475200154185814173223 -1746.04780.100.120470.003834.438040.078590.267270.00417196356172015152721 -1828.63090.090.120370.003744.821020.079450.290580.00446196254178914164422 -1972.12830.250.118300.003775.894820.106120.361520.00568193156196116198927 -2082.73970.210.125780.003915.092220.083720.293730.00451204054183514166022 -2140.86130.070.117760.003704.417460.075410.272130.00422192355171614155221 -2242.73900.110.122640.003855.572060.094820.329600.00513199555191215183625 -231211690.710.120650.003896.169610.114390.370950.00591196656200016203428 -2460.43070.200.121040.003836.257510.109450.375040.00589197255201315205328 -2538.13430.110.121480.003815.620510.095990.335630.00523197855191915186625
表2 贺兰山北段古元古代花岗岩锆石Lu-Hf同位素分析结果
Table 2 Zircon Lu-Hf isotope compositions of the Paleoproterozoic granites in the northern Helanshan
测点号Age(Ma)176Yb177Hfσ176Lu177Hfσ176Hf177HfσεHf(0)εHf(t)tDM1(Ma)tDM2(Ma)fLu/Hf08HL100 片麻状含石榴黑云母花岗岩-01-02-03-05-07-08-09-10-11-14-1519910.0140160.00034600.0004250.0000120.2815100.000016-44.63-0.724002651-0.990.0156000.00050200.0004630.0000150.2816030.000016-41.342.622762452-0.990.0313990.00065500.0008870.0000190.2815150.000020-44.45-1.124222679-0.970.0197460.00005100.0006790.0000050.2815650.000017-42.681.023412552-0.980.0296210.00071400.0009270.0000250.2815510.000020-43.180.123752603-0.970.0422010.00016900.0012270.0000040.2815650.000016-42.680.223742598-0.960.0376020.00068400.0010170.0000200.2815620.000016-42.790.423662587-0.970.0259190.00015400.0008980.0000070.2815820.000018-42.081.323312533-0.970.0228160.00024100.0007060.0000060.2815740.000016-42.371.223302535-0.980.0196420.00012100.0006190.0000020.2815800.000023-42.151.623172515-0.980.0226210.00004900.0006040.0000010.2815220.000018-44.21-0.523942640-0.9808HL95 似斑状石榴黑云母花岗岩-06-07-19-11-12-14-15-16-19-20-22-23-24-2519730.0233780.00015600.0006650.0000040.2815360.000018-43.71-0.523792625-0.980.0462130.00114000.0012850.0000300.2816320.000018-40.322.122862467-0.960.0167710.00012700.0004660.0000030.2816160.000014-40.882.622592434-0.990.0274100.00029400.0007760.0000090.2816230.000015-40.632.522682445-0.980.0655820.00190000.0016800.0000450.2815280.000018-43.99-2.124552726-0.950.0217040.00005400.0006350.0000020.2816200.000016-40.742.622632440-0.980.0226780.00006600.0006590.0000030.2815850.000014-41.981.323122518-0.980.0309310.00042800.0008710.0000120.2815880.000016-41.871.123212529-0.970.0266500.00011900.0007340.0000040.2816190.000015-40.772.422712450-0.980.0299250.00008400.0008010.0000030.2815770.000016-42.260.823322547-0.980.0408050.00023400.0011500.0000070.2816160.000017-40.881.722992491-0.970.0267800.00032900.0007280.0000080.2815060.000017-44.77-1.624242696-0.980.0246470.00003000.0006860.0000010.2816570.000018-39.433.822162363-0.980.0392880.00017400.0010900.0000070.2815770.000017-42.260.423492571-0.97
注:176Lu的衰变常数为1.867×10-11yr-1(Albarèdeetal., 2006); (176Lu/177Hf)CHUR=0.0332±0.0002, (176Hf/177Hf)CHUR=0.0282772±0.000029 (Bouvieretal., 2008); (176Lu/177Hf)DM=0.0384, (176Hf/177Hf)DM=0.28325(Griffinetal., 2000);εHf(0)=((176Hf/177Hf)s/(176Hf/177Hf)CHUR,0-1)×10000;εHf(t)=((176Hf/177Hf)s-(176Lu/177Hf)s×(eλt-1))/((176Hf/177Hf)CHUR,0-(176Lu/177Hf)CHUR×(eλt-1))-1)×10000;tDM1=1/λ×ln[1+((176Hf/177Hf)s-(176Hf/177Hf)DM)/((176Lu/177Hf)s-(176Lu/177Hf)DM)];tDM2=tDM1-(tDM1-t)((fcc-fs)/(fcc-fDM));fLu/Hf=(176Lu/177Hf)s/(176Lu/177Hf)CHUR-1; 其中,两阶段Hf模式年龄(tDM2)计算时,采用地壳平均值176Lu /177Hf=0.015 (Rudnick and Gao, 2003)
图3 贺兰山北段花岗岩锆石CL图像(a、c、e)和U-Pb年龄谐和图(b、d、f)Fig.3 Cathodoluminescence (CL) images of representative zircons (a, c, e) and zircon U-Pb concordia diagrams (b, d, f) for the granites from the northern Helanshan area
图4 贺兰山北段花岗岩(Na2O+K2O)-SiO2图解(a,据 Middlemost, 1994)、 A/CNK-A/NK图解(b,据 Maniar and Piccoli, 1989)、 SiO2-K2O图解(c,据 Peccerillo and Taylor, 1976)及 (Al2O3-(K2O+Na2O))-CaO-(FeOT+MgO) 图解(d,据 Chappell and White,1974)Fig.4 Na2O+K2O vs. SiO2 classification diagram(a, after Middlemost, 1994), A/CNK vs. A/NK diagram(b, after Maniar and Piccoli, 1989), SiO2 vs. K2O diagram(c, after Peccerillo and Taylor, 1976)and (Al2O3-(K2O+Na2O))-CaO-(FeOT+MgO) diagram(d, after Chappell and White,1974)for the granites in the northern Helanshan
片麻状含石榴黑云母花岗岩(08HL100)中的锆石呈短柱状,锆石粒度约30~100μm,长宽比约为1:1~3:1。CL图像显示锆石阴极发光性较差,大部分锆石显示模糊的韵律环带(图3c)。锆石的Th、U含量分别为55.6×10-6~680×10-6、123×10-6~718×10-6,Th/U比值主要在0.27~2.31之间,均大于0.1,具有典型的岩浆成因锆石的特征。选取16粒锆石进行测试分析(表1),获得其207Pb/206Pb上交点年龄为1991±29Ma(MSWD=3.5,n=16;图3d),该结果解释为此片麻状含石榴黑云母花岗岩的侵位年龄。
似斑状石榴黑云母花岗岩(08HL95)中的锆石呈半透明,半自形-自形长柱状,磨圆度相对较差,锆石粒度约为50~200μm,长宽比约为2:1~4:1。CL图像显示锆石呈暗灰色,成分环带模糊,部分锆石颗粒显示较好的韵律环带(图3e)。锆石的Th、U含量分别为28.6×10-6~149×10-6、81.3×10-6~613×10-6,Th/U比值除了18和21号测点小于0.1之外,其余测点的Th/U比值均大于0.1(0.10~0.71),结合锆石CL图像,判断为岩浆成因的锆石(Belousovaetal., 2002; Williamsetal., 2009)。选取25粒锆石进行测试分析,获得其207Pb/206Pb上交点年龄为1973±23Ma(MSWD=3.4,n=25;图3f),该年龄应代表了似斑状石榴黑云母花岗岩的结晶年龄。
片麻状含石榴黑云母花岗岩(08HL100)和似斑状石榴黑云母花岗岩(08HL95)中锆石原位Hf同位素分析结果见表2。
对片麻状含石榴黑云母花岗岩样品中完成锆石U-Pb测试分析的11颗锆石进行原位Hf同位素测试分析。结果显示:锆石的176Lu/177Hf的比值为0.000425~0.001227,176Hf/177Hf的值为0.281510~0.281603,利用上交点年龄(t=1991Ma)计算获得其εHf(t)值为-1.1~+2.6(平均值为+0.6),单阶段Hf模式年龄为(tDM1)2.28~2.42Ga,两阶段模式年龄(tDM2)为2.45~2.68Ga。
图5 贺兰山北段花岗岩球粒陨石标准化稀土元素配分图(a、c)和原始地幔标准化微量元素蛛网图(b、d)(标准化值据Sun and McDonough, 1989)Fig.5 Chondrite-normalized REE pattern diagrams (a, c) and primitive mantle-normalized trace element spidergrams (b, d) for the granites in the northern Helanshan area (normalization values after Sun and McDonough, 1989)
对似斑状石榴黑云母花岗岩样品中完成锆石U-Pb年龄分析的14颗锆石进行原位Hf同位素测试分析。获得其176Lu/177Hf的比值为0.000466~0.001685,176Hf/177Hf的值为0.281506~0.281657,利用上交点年龄(t=1973Ma)计算获得其εHf(t)值为-2.1~+3.8(平均值为+1.2),单阶段Hf模式年龄为(tDM1)2.22~2.46Ga,两阶段模式年龄(tDM2)为2.36~2.73Ga。
3.3.1 主量元素特征
选取6件代表性的样品进行主量、微量元素分析测试,其结果见表3。片麻状含石榴黑云母花岗岩具有低SiO2(60.88%~65.42 %), 高Al2O3(16.84%~19.35 %)、 Fe2O3T(5.94%~6.14%)和K2O(3.25%~5.10%),K2O/Na2O=1.33~4.47,相对较高的Mg#(47~51)。而似斑状石榴黑云母花岗岩具有相对较高的SiO2(71.93%~74.98%)、K2O(4.58%~6.98%)、Al2O3(13.25%~14.79%)和NaO2(3.13%~3.94%), K2O/Na2O=1.16~2.23,相对富钾,中等Mg#值(17~43)。在(Na2O+K2O)-SiO2图解中(图4a),片麻状含榴黑云母花岗岩样品点落入闪长岩-花岗闪长岩区域内,其余样品点落入花岗岩区域内,贺兰山花岗岩与凉城花岗岩具有相类似的硅碱含量特征。在A/CNK-A/NK图解(图4b)贺兰山花岗岩的铝饱和指数A/CNK=1.08~2.34,刚玉含量(C)相对较大(1.11~12.5),表明样品属于过铝质。里特曼指数σ=1.45~3.53,平均为2.44,在SiO2-K2O图解中(图4c),样品位于高钾钙碱性-钾玄质系列区域。在(Al2O3-(K2O+Na2O))-CaO-(FeOT+MgO)图解中(图4d),贺兰山花岗岩样品落入S型花岗岩区域内。以上特征显示,本文研究的贺兰山两类花岗岩属于过铝质S型花岗岩。
3.3.2 稀土及微量元素特征
研究区内的片麻状含石榴黑云母花岗岩具有较高的稀土元素含量(∑REE=262×10-6~290×10-6)。在球粒陨石标准化稀土元素配分模式图上(图5a),样品显示右倾型的模式,轻、重稀土分异较为明显((La/Yb)N=11.6~14.5),且Eu负异常明显(δEu=0.30~0.55),表明源区斜长石残留或在岩浆演化中经历了较为显著的斜长石结晶分异。似斑状石榴黑云母花岗岩具有相对较低的稀土元素总量(∑REE=35.3×10-6~127×10-6),中等轻、重稀土分异程度((La/Yb)N=1.79~13.6),且具有明显的Eu负异常或正异常(δEu=0.07~3.58)(图5c)。在原始地幔标准化微量元素模式图上(图5b,d),贺兰山早期花岗岩和晚期花岗岩均显示明显富集大离子亲石元素(如Rb、K、U等),亏损高场强元素 (如Nb、 Ta等) 的特征,但早期片麻状含石榴黑云母花岗岩中具有相对高的Zr、Hf元素。
表3 贺兰山北段花岗岩样品位置及其主量(wt%)、微量(×10-6)元素分析结果
Table 3 Major (wt%) and trace (×10-6) elements of the granites in the northern Helanshan and sample locations
样品号08HL9308HL10008HL8808HL9108HL9608HL101岩性片麻状含石榴黑云母花岗岩似斑状石榴黑云母花岗岩GPS39°13'23.75″N105°58'47.52″E39°20'25″N106°15'54″E39°15'30.76″N106°02'19″E39°15'25.86″N106°05'14.79″E39°16'04.42″N106°03'03.50″E39°21'07″N106°15'20″ESiO260.8865.4271.9373.2974.9872.78TiO20.810.650.020.070.040.07Al2O319.3516.8414.7914.4713.2514.37Fe2O3T5.946.140.810.691.141.39MnO0.040.080.010.030.010.03MgO2.632.350.260.160.100.30CaO0.491.140.560.700.460.69Na2O1.142.453.133.943.423.16K2O5.103.256.984.585.386.08P2O50.090.080.220.170.050.12LOI3.261.530.831.740.680.58Total99.7399.9399.5499.8499.5199.57Mg#514743351734A/CNK2.341.751.081.141.081.10C12.58.031.632.221.111.66σ2.181.453.532.402.422.87K2O/Na2O4.471.332.231.161.571.92Li51.547.55.976.976.059.36Be2.7916.50.631.5811.50.98Sc9.7910.82.722.489.096.13V94.079.32.552.981.968.47Cr89.591.51.131.971.5013.0Co28.854.345.046.739.654.8Ni27.836.82.281.691.847.25Cu3.318.452.371.802.432.39Zn77.281.56.568.6916.710.8Ga27.522.713.116.120.214.1Ge1.551.631.362.002.491.33Rb175115179145451156Sr56.517310153.312.0159Y29.937.24.4810.318.417.0Zr40630611.433.951.2116Nb23.616.30.313.4641.91.89Cs3.562.271.250.6730.00.63Ba59976251312334.8673La59.355.09.698.107.1733.0Ce12511415.817.817.555.7Pr15.413.51.441.991.905.61Nd55.646.84.357.155.8318.1Sm10.17.800.661.502.172.81Eu0.901.300.730.240.051.13Gd8.446.750.591.321.972.37Tb1.121.040.090.240.460.36Dy5.806.490.671.582.942.51Ho1.021.300.150.330.580.59Er2.903.750.491.081.912.03Tm0.440.540.080.190.360.35Yb2.933.390.511.382.882.46Lu0.470.520.070.210.470.41
续表3
Continued Table 3
样品号08HL9308HL10008HL8808HL9108HL9608HL101岩性片麻状含石榴黑云母花岗岩似斑状石榴黑云母花岗岩GPS39°13'23.75″N105°58'47.52″E39°20'25″N106°15'54″E39°15'30.76″N106°02'19″E39°15'25.86″N106°05'14.79″E39°16'04.42″N106°03'03.50″E39°21'07″N106°15'20″EHf11.48.370.391.442.723.39Ta1.590.960.150.306.200.26Pb9.2922.432.825.733.835.1Th22.416.30.751.783.894.74U3.842.230.642.107.451.10δEu0.300.553.580.510.071.33(La/Yb)N14.511.613.64.211.799.62∑REE29026235.343.146.2127
图6 贺兰山地区古元古代岩浆年龄(a)和变质年龄(b)频率分布图(数据来自表 4) Fig.6 Histogram diagrams of Paleoproterozoic magmatic (a) and metamorphic (b) ages from the Helanshan area (data from Table 4)
近年来,很多学者对华北克拉通西部的孔兹岩系进行了高精度锆石U-Pb年代学的研究,结果表明孔兹岩系原岩中的碎屑物质主要以2.2~2.0Ga的物质为主,源区物质的峰期年龄为2.03Ga,其中最小的碎屑锆石年龄约为~2.00Ga,变质时代约为1.95~1.92Ga,推测其孔兹岩系原岩沉积时间在2.0~1.95Ga之间,属于古元古代的产物(Zhaoetal., 2005, 2010; Xiaetal., 2006; 董春艳等,2007;Wanetal., 2009; Yinetal., 2009,2011; Danetal., 2012;李黎明等,2014;Qiaoetal., 2016),而不是之前所认为的新太古代(沈其韩等, 1990;胡能高和杨家喜,1995;Qian and Li, 1999)。锆石U-Pb年代学研究表明,孔兹岩带东段乌拉山-大青山孔兹岩系变质年龄峰值主要集中在~1.95Ga 和~1.86Ga(Zhaoetal., 2010; Dongetal., 2014; 蔡佳等,2015)。孔兹岩带西段贺兰山群峰期变质年龄大约为~1.95Ga,此外,还记录了1.92~1.90Ga的变质年龄,此年龄解释为变质峰期后降压变质作用的时代(Yinetal., 2009, 2011;Qiaoetal., 2016)。
耿元生等(2009)将华北克拉通北缘古元古代晚期花岗岩带划分为三个主要时段:大于2.0Ga、2.0~1.87Ga和1.85~1.80Ga。已有的研究表明,贺兰山地区的岩浆活动主要集中在~2.05Ga、~1.95Ga和~1.85Ga(李俊建等,2004; 耿元生等,2009; Yinetal., 2011; Danetal., 2012; Lietal., 2017, 2018; Zhangetal., 2017; 张道涵,2017),而1.85~1.80Ga的岩浆活动在贺兰山-千里山地区记录相对较少(Yinetal., 2009,2011),反而在孔兹岩带东部较为发育 (王惠初等,2005; Zhaoetal., 2008; 马铭株等,2012; Liuetal., 2013; Wanetal., 2013; 陈佩嘉等,2017)。
图7 贺兰山北段花岗岩源区判别图解(a) (Na2O+K2O+MgO+TiO2+FeO)-(Na2O+K2O)/(MgO+TiO2+FeO) (据 Patio Douce,1999);(b) Al2O3/TiO2-CaO/Na2O;(c)Rb/Sr-Rb/Ba(b, c 据 Sylvester,1998);图例同图4Fig.7 Discrimination diagrams of granites source area in the northern Helanshan area
研究区内出露的S型花岗岩在空间上与大面积出露的孔兹岩系呈复杂构造接触关系,并且一些研究者认为该区古元古代时期的花岗岩(~1.95Ga)是孔兹岩系部分熔融的结果(Yinetal., 2009;Danetal., 2014;Zhangetal., 2017)。贺兰山孔兹岩具有较高的∑REE含量(162×10-6~337×10-6),Eu负异常明显(δEu=0.30~0.65),轻、重稀土强烈分异,相对富集K、Rb等大离子亲石元素,具有较低的Ti、Cr、Sr、Y等元素(胡能高和杨家喜,1995)。如图5所示,本文所研究的S型花岗岩稀土和微量元素配分模式演化特征与上地壳演化特征相类似,而与贺兰山孔兹岩的稀土、微量元素特征明显不同。需要指出的是:贺兰山花岗岩稀土元素含量与孔兹岩类似,甚至更低,特别是晚期花岗岩具有较低的稀土含量(图5c)。总体来讲,贺兰山两期花岗岩相对富集轻稀土元素和大离子亲石元素(如:K、Rb、Sr等),亏损高场强元素(如:Nb、Ta等),但是早期片麻状黑云母花岗岩相对富集Zr和Hf元素。Zr和Hf均属于高场强元素,它们具有相类似的离子半径和相同的离子价态,并且被认为具有相类似的地球化学特征。研究表明由于基性-碱性岩浆中具有较高的Zr溶解度和熔融温度,所以在这些岩浆岩中几乎不含继承锆石(赵子福等,2013; Mengetal., 2019)。但是在非碱性、花岗质岩浆中,Zr 的溶解度较低(Watson and Harrison, 1983),尤其是对于粒度较大的锆石晶体(>50~100μm),在地壳物质受部分熔融过程时可能得以保存并最终进入到花岗质岩浆中(Watson, 1996),特别是在S型花岗岩中,继承锆石较为常见(赵子福等,2013),所以贺兰山早期花岗岩微量元素中Zr和Hf轻微富集可能是由于锆石中存在一些继承锆石所致。
岩石学和地球化学特征表明,S型花岗岩通常是由地壳物质部分熔融所形成的(Barbarin, 1999; Chappell, 1999),其源岩主要为变泥质岩和杂砂岩(Sylvester, 1998; Clemens, 2003)。实验岩石学表明,过铝质的S型花岗岩CaO/Na2O比值可能与斜长石/粘土矿物的比值有关,若CaO/Na2O>0.3,则源区为杂砂岩的部分熔融;反之,则源区为泥质岩石的部分熔融(Sylvester,1998)。在花岗岩源区判别图解中(图7a),片麻状含石榴黑云母花岗岩样品点靠近变质杂砂岩区域,似斑状石榴黑云母花岗岩样品点主要落入长英质泥质岩区域内。在Al2O3/TiO2-CaO/Na2O图解中(图7b),片麻状含石榴黑云母花岗岩样品点落入砂砾岩源区,而似斑状石榴黑云母花岗岩样品点落入泥质岩源区。本文中的花岗岩样品除了08HL96样品具有高的Rb/Sr和Rb/Ba比值外,其余样品均具有低Rb/Sr(0.67~3.18)和Rb/Ba(0.15~0.30),其源区主要为粘土质源区(Sylvester,1998;图7c)。以上地球化学特征表明贺兰山花岗岩原岩是由地壳泥质岩夹部分杂砂岩部分熔融形成的。
图8 εHf(t)-锆石U-Pb年龄图解Fig.8 U-Pb age vs. εHf(t) diagram
贺兰山两期花岗岩中所测锆石的176Lu/177Hf比值均小于0.002,表明由176Lu衰变形成的放射性成因177Hf同位素的积累很小,可忽略不计。所以实际测得176Hf/177Hf的比值代表了锆石形成时的初始Hf同位素组成。贺兰山两期S型花岗岩具有较低的176Hf/177Hf比值(0.281506~0.281657),其εHf(t)值分别为-1.1~+2.6和-2.1~+3.8,Hf两阶段模式年龄分别为2.45~2.68Ga和2.36~2.73Ga,明显老于岩石的形成年龄(图8;表2)。在εHf(t)-t图解中(图8),贺兰山两期花岗岩中锆石的εHf(t)值明显低于同时期亏损地幔的值,主要分布在球粒陨石演化线附近,表明这些花岗岩的源区主要是以古老地壳物质部分熔融为主,可能有少量的新生地壳物质加入。由于贺兰山副变质岩碎屑锆石的εHf(t)变化较大(-2.9~+8.9),锆石大多数具有正的εHf(t)值,指示其源区性质可能以亏损地幔物质为主。此外,这些碎屑锆石的δ18O 的峰值为6.6‰和8.2‰,进一步表明孔兹岩的碎屑物质源区以新生的地壳物质为主,并有少量古老地壳物质加入,指示其可能形成于活动大陆边缘的构造环境(Yinetal., 2009, 2011; 但卫,2011;Danetal., 2012)。同时,如上所述贺兰山花岗岩稀土元素含量与孔兹岩类似,甚至更低,由此推断贺兰山北部宗别立地区的片麻状含石榴黑云母花岗岩和似斑状石榴黑云母花岗岩可能并不是由孔兹岩部分熔融的而成,而是古老变沉积岩部分熔融的产物。
近年来,高精度年代学工作研究表明华北克拉通西部陆块孔兹岩带中的孔兹岩原岩沉积时间为2.0~1.95Ga(Xiaetal., 2009; Yinetal., 2011;校培喜等,2011;Qiaoetal., 2016; Lietal., 2017),变质作用时限约为1.95~1.92Ga(周喜文和耿元生,2009;Yinetal., 2009,2011,2014,2015;Zhaoetal., 2010;但卫,2011; Danetal., 2012;Guoetal., 2012;Jiaoetal., 2013)。一些学者对华北克拉通西部孔兹岩带东部的千里山群,乌拉山群和集宁群中的孔兹岩及贺兰山地区高压泥质麻粒岩进行了变质演化过程的研究,结果显示都具有相类似的顺时针PT轨迹演化的特征,指示了一个陆-陆碰撞,之后又快速折返并抬升到地壳的动力学过程(周喜文等,2010;Wangetal., 2011;Yinetal., 2014, 2015;Caietal., 2015),这从岩石学的角度进一步表明了华北克拉通西部孔兹岩带是由阴山地块和鄂尔多斯地块在~1.95Ga期间碰撞所形成的(赵国春,2009;周喜文等,2010;Yinetal., 2011)。
贺兰山地区古元古代晚期的岩浆活动以发育S型花岗岩为主。可是,诸多研究者对这些不同时期花岗岩形成环境的解释不尽相同。一些研究者认为贺兰山北段~2.0Ga的S型花岗岩(张道涵,2017)和中段小口子花岗岩(Lietal., 2018)的地球化学特征可能反映了一个与俯冲相关的陆缘弧环境。另外一些研究者发现贺兰山北部~1.95Ga的S型花岗岩具有较高的锆石δ18O 值(8.7‰~10.6‰),变化范围较小的εHf(t)值(-1.0~+0.8),相对应的二阶段Hf模式年龄(tDM2)为2.41~2.66Ga,这些S型花岗岩的成因模式可用俯冲板片-断离的构造机制来解释(Danetal., 2014)。还有一些研究者基于孔兹岩带中记录了~1.92Ga的超高温变质作用事件和发现了同时期的辉长苏长岩(Pengetal., 2010; Guoetal., 2012),认为该区~1.92Ga的地质事件可能是由碰撞造山后拉伸所引起的地幔岩浆上涌底侵所形成的(赵国春,2009;Zhaoetal., 2012;Zhangetal., 2017)。综合已有的区域地质资料,推测华北克拉通西部贺兰山地区在古元古代晚期可能处于活动大陆边缘的构造环境(Wanetal., 2009; Yinetal., 2009; Danetal., 2012; Zhangetal., 2015, 2017; Lietal., 2017, 2018)。贺兰山北部宗别立地区的片麻状含石榴黑云母花岗岩和似斑状石榴黑云母花岗岩是古元古代晚期碰撞造山相关的两期岩浆作用的产物,是由变沉积岩部分熔融形成的,且片麻状含石榴黑云母花岗岩的侵位年龄(~2.0Ga)明显早于孔兹岩的峰期变质时间(~1.95Ga),其岩浆活动形成于陆陆碰撞之前,推测这次岩浆事件可能是由洋壳俯冲所引起的;而晚期似斑状石榴黑云母花岗岩(~1.97Ga)的侵位时间与区内孔兹岩峰期变质时间大致相同或稍早,可能为同碰撞阶段岩浆作用的产物,反映两个微陆块碰撞时限可能稍早于~1.95Ga。总之,本文所研究的两期花岗岩记录了华北克拉通阴山地块和鄂尔多斯地块从初始俯冲-碰撞到同碰撞的重要信息,这为理解华北克拉通西部陆块的拼合碰撞演化时限提供了新的约束。
(1)华北克拉通西部贺兰山地区出露的片麻状含石榴黑云母花岗岩和似斑状石榴黑云母花岗岩侵位年龄分别为~2.0Ga 和~1.97Ga,为贺兰山地区古元古代不同阶段岩浆作用的产物。
(2)地球化学特征显示贺兰山两期花岗岩属于过铝质S型花岗岩。Hf同位素结果表明它们的原岩主要来自于古老变沉积岩,并非为孔兹岩部分熔融产物。
(3)贺兰山地区在古元古代晚期可能为陆缘弧的构造环境,早期片麻状含石榴黑云母花岗岩和晚期似斑状石榴黑云母花岗岩分别形成于阴山地块和鄂尔多斯地块俯冲-碰撞前和同碰撞的构造环境。
致谢 耿元生研究员、杜利林研究员、王洛娟博士、张华锋副教授及匿名审稿人对本文进行了详细的评审和指导,并提出了宝贵的修改建议,在此表示衷心感谢!