青藏高原纳木错西缘新元古代中期岩浆事件:对北拉萨地块起源的约束*

胡培远 翟庆国 赵国春 唐跃 王海涛 朱志才 王伟 吴昊

1. 自然资源部深地动力学重点实验室，中国地质科学院地质研究所，北京 1000372. 香港大学地球科学系，香港3. 大陆动力学国家重点实验室，西北大学地质学系，西安 7100691.

图1 青藏高原中部构造划分简图 NGSZ-北冈底斯缝合带；BNSZ-班公湖-怒江缝合带；IYZSZ-雅鲁藏布江缝合带Fig.1 Simplified tectonic map of the central Tibetan Plateau

青藏高原是地球上最大最高的高原，其演化历史涉及到多个特提斯洋盆的打开、俯冲和消亡以及多个微陆块的裂离、漂移和拼贴过程(Yin and Harrison, 2000; Metcalfe, 2013; Zhuetal., 2013)，并且至今仍然处于活动状态，因此在构造地质学研究中处于特殊的地位。拉萨地块位于青藏高原中部，大地构造位置上夹持于班公湖-怒江和雅鲁藏布江缝合带之间，是组成青藏高原的主体地块之一(Yin and Harrison, 2000; Zhangetal., 2012a; Yangetal., 2009; Chenetal., 2009)。传统观点认为，拉萨地块起源于印度大陆东北缘(Yin and Harrison, 2000; Gehrelsetal., 2011; Guoetal., 2017)或者澳大利亚大陆北缘(Audley-Charles, 1983; Allègre, 1984; Zhuetal., 2011)。然而，新近在拉萨地块中部识别出了一条古特提斯板块缝合带——北冈底斯板块缝合带(图1)，其岩石组合包括蛇绿混杂岩、岛弧岩浆岩和榴辉岩等(Yangetal., 2009; Chenetal., 2009; Chengetal., 2012; Caoetal., 2017)。因此，拉萨地块可以被进一步划分为北拉萨和南拉萨两个地块，它们的起源和早期构造演化历史有必要重新考虑。

近年来，有学者对南拉萨地块的前寒武纪演化历史开展了较为系统的研究，结果表明该陆块起源于印度大陆东北缘(Linetal., 2013; Xuetal., 2013; Guoetal., 2017)。然而，北拉萨地块的起源依然存有争议。目前主要有以下两种不同观点：(1)南、北拉萨地块具有类似的前寒武纪演化历史，均起源于印度大陆东北缘(Yin and Harrison, 2000; Gehrelsetal., 2011; Guoetal., 2017)；(2)南、北拉萨地块具有不同的前寒武纪演化历史，其中北拉萨地块可能起源于澳大利亚大陆北缘(Ferrarietal., 2008; Zhuetal., 2011)或者东非造山带北段(Zhangetal., 2012a; Huetal., 2018a, b, c)。前寒武纪岩浆和变质事件的区域对比是解决这一分歧的有效方法之一。

图2 纳木错西缘地区区域地质简图 前人年龄(单位：Ma)资料引自张修政等(2013); Hu et al. (2005, 2018b, c); Dong et al. (2011); Zhang et al. (2012a)Fig.2 Geological map of the western Nam Tso area

念青唐古拉岩群自提出以来一直被认为是北拉萨地块古老变质基底的典型代表(李璞, 1955)。从东部的工布江达、到中部的纳木错西缘和念青唐古拉山，再到西部的日土地区，念青唐古拉岩群大致覆盖了整个北拉萨地块(图1)。同时，念青唐古拉岩群中包含较多的基性岩(Huetal., 2005; 胡培远等, 2016)，这为探讨其形成的构造环境提供了较好的素材。本文试图通过对纳木错西缘地区念青唐古拉岩群中新元古代中期变质辉长岩和花岗片麻岩的全岩地球化学、同位素年代学以及锆石Hf同位素的研究，揭示其岩浆源区特征以及构造背景，从而为探索拉萨地块的起源提供新资料。

1 地质概况

前人研究表明，在青藏高原中部地区存在3个地块，分别为南羌塘、北拉萨和南拉萨地块(Yin and Harrison, 2000; Yangetal., 2009; Zhangetal., 2012a)。这些地块被班公湖-怒江和北冈底斯板块缝合带分隔。同时，在班公湖-怒江缝合带内部还有一个相对较小的陆块——安多微陆块，夹持于蛇绿岩带内(Zhangetal., 2012b; Xieetal., 2013; 解超明等, 2014;王明等, 2012; 辜平阳等, 2012)。安多片麻岩(Guynnetal., 2012)、念青唐古拉岩群(胡培远等, 2016)和林芝岩群(Linetal., 2013)分别代表了安多微陆块、北拉萨地块和南拉萨地块的前寒武纪基底(图1)。

念青唐古拉岩群由李璞(1955)所称的念青唐古拉片麻岩和那更拉片岩系演变而来，主要由三部分岩石组成：(1)副片麻岩，以片岩为主，主要岩石类型为白云母石英片岩、白云母石英岩、长石石英岩、黑云母石英片岩等；(2)正片麻岩，原岩以基性岩和酸性岩为主，包含少量的中性岩，主要岩石类型为变质辉长岩(斜长角闪岩)和花岗片麻岩；(3)后期侵入脉体，广泛分布于副片麻岩和正片麻岩中，以花岗质岩石为主，也可见少量的辉长岩，变质变形程度较弱(胡道功等, 2003; Huetal., 2005; 胡培远等, 2016)。

虽然念青唐古拉岩群广泛分布于工布江达、纳木错西缘、念青唐古拉山、日土等地区，但是不同地区的研究程度差别较大。在工布江达和日土地区，念青唐古拉岩群主要出露于高山峡谷区，开展研究工作十分困难，研究程度很低；同时，这两个地区位于拉萨地块的边缘位置，能否代表拉萨地块的基底性质尚未可知。在念青唐古拉山一带，虽然出露有较大规模的片麻岩系，但是曾有学者在这一地区的正片麻岩中获得了古新世的岩浆结晶年龄(64～59Ma)(胡道功等, 2003)，从而对该地区念青唐古拉岩群的时代提出了质疑。目前，念青唐古拉岩群中确切的新元古代构造岩浆记录主要集中于纳木错西缘地区。Huetal. (2005)曾在纳木错西缘念青唐古拉岩群中识别出了新元古代拉斑玄武岩和花岗岩，其中拉斑玄武岩的继承锆石给出1766～988Ma的中元古代年龄，而花岗岩则获得了748Ma和787Ma的谐和年龄；张修政等(2013)报道该地区的念青唐古拉岩群中存在E-MORB型斜长角闪岩(约758Ma)和近同时代的弧岩浆岩(约742～731Ma)；Dongetal. (2011)报道念青唐古拉岩群中斜长角闪岩的原岩形成于850Ma左右；Huetal. (2018b, c)在邻近地区报道了822～760Ma的弧后盆地岩浆岩；Zhangetal. (2012a)在纳木错西缘念青唐古拉岩群中发现了基性高压麻粒岩，认为其原岩形成时代在897～886Ma左右，并且具有类似于N-MORB的地球化学特征；最近，这一900Ma左右的N-MORB型基性岩浆事件得到了进一步的确认(胡培远等, 2016; Huetal., 2018d; Zengetal., 2018)。

图4 纳木错西缘地区变质辉长岩(a)和花岗片麻岩(b)中典型锆石的阴极荧光图像(年龄单位：Ma)Fig.4 CL images of typical zircons of the dated meta-gabbro (a) and granitic gneiss (b) samples from the western Nam Tso area (age unit: Ma)

本文研究的变质辉长岩和花岗片麻岩侵入于念青唐古拉岩群副片麻岩中，出露于仁错东南约20km，纳木错以西约45km(图2)。岩相学观察发现，辉长岩发生了较为明显的角闪岩相变质作用，辉石部分已经变质成为角闪石，但是仍可见残余的辉长结构(图3a)；岩石矿物组合以斜长石(40%～45%)、角闪石(30%～35%)和辉石(25%～30%)为主(图3b)，此外还有少量绿帘石、钛铁矿、榍石、磁铁矿等，未见石榴石；矿物粒度为0.1～0.6mm。花岗片麻岩总体呈灰白色(图3c)，具细粒花岗结构，片麻状构造。经镜下鉴定，主要矿物组成为石英(40%～50%)和斜长石(35%～55%)，钾长石含量较少(<5%)；斜长石主要为钠长石，普遍发育聚片双晶；副矿物有锆石、磷灰石、褐帘石和榍石等(图3d)。

2 样品测试方法

锆石U-Pb测年样品中锆石的分选在河北省区域地质调查院完成，采用常规的重液和磁选方法进行分选，最后在双目显微镜下挑纯。样品靶的制备在中国地质科学院地质研究所完成，制成的样品靶直径为25mm。锆石的阴极荧光图像分析在中国地质科学院地质研究所的阴极荧光分析系统(HITACH S-3000N型场发射环境扫描电镜和Gatan公司Chroma阴极荧光谱仪)上完成。样品的锆石U-Th-Pb分析在北京科荟测试技术有限公司完成，分析仪器为美国ESI公司生产的NWR 193nm激光剥蚀进样系统和德国AnlyitikJena公司生产的PQMS Elite型四级杆等离子体质谱仪联合构成的激光等离子体质谱仪(LA-ICP-MS)。本次分析193nm激光器工作频率为10Hz；测试点束斑直径为25μm，剥蚀采样时间为45s，具体分析流程见侯可军等(2009)。锆石GJ-1(Jacksonetal., 2004)作为外部标准来校正分析过程中的同位素分馏；NIST610作为外部标准来获得分析点的Th和U的含量。锆石U-Pb年龄用ICPMSDataCal数据处理软件(Liuetal., 2010)计算获得，加权平均年龄的计算和谐和图的绘制采用ISOPLOT 3.0程序(Ludwig, 2003)。锆石Hf同位素分析在中国科学院地质与地球物理研究所Neptune多接收电感耦合等离子质谱仪(MC-ICPMS)和193nm激光取样系统上进行， 仪器的运行条件及详细的分析过程参见Wuetal. (2006)。采用单点剥蚀模式，斑束固定为44μm。实验测定过程中，MUD标准锆石的176Hf/177Hf的测定结果是0.282505±21，与前人获得的结果一致(Wuetal., 2006)。地球化学样品的主量元素、微量元素和稀土元素的分析均在国家地质实验测试中心。主量元素采用X-射线荧光光谱仪(PW4400)分析。微量元素和稀土元素的分析仪器为X-series等离子质谱仪，实验室分析详细方法见相关参考文献(Zhaietal., 2016)。

表1纳木错西缘新元古代中期岩浆岩的锆石LA-ICP-MSU-Pb-Th分析结果

Table 1 LA-ICP-MS U-Pb-Th data for zircons from the Middle Neoproterozoic magmatic rocks in the western Nam Tso area

测点号含量(×10-6)Pb232Th238UTh/U同位素比值年龄(Ma)207Pb/206Pb207Pb/235U206Pb/238U207Pb/206Pb207Pb/235U206Pb/238U比值1σ比值1σ比值1σ年龄1σ年龄1σ年龄1σ16T160 变质辉长岩 30°49.128′N、90°01.946′E1.1271091610.680.06430.00101.05400.02420.11890.00247503373112724142.122761400.540.06340.00091.03290.02470.11790.00247242772012719143.1301051880.560.06400.00081.04380.02410.11830.00267432772612721154.1341292120.610.06390.00071.03730.02170.11790.00247391972311718145.1341402130.660.06750.00171.08200.04050.11580.00248545474520707146.1341302040.640.06290.00081.02660.01810.11870.0021703287179723127.1381342410.560.06500.00091.05150.02200.11730.002177622773011715128.121711320.540.06490.00101.06560.02900.11900.00297723373714725169.1261141530.750.06590.00111.07000.02180.11800.0018120037739117191110.128761760.430.06500.00111.06000.02630.11930.003177234734137271811.1281161760.660.06400.00091.04530.02380.11910.002773931727127251512.119691180.580.06600.00111.06970.02670.11750.002380634739137161313.1271271680.760.06600.00091.07180.02410.11810.002580630740127191414.1351362070.660.06560.00101.05990.02120.11750.002379436734107161315.1241001430.700.06390.00101.04460.02070.11860.001973935726107231116.140592700.220.06640.00081.07610.02350.11760.002582031742117171417.123961370.700.06490.00121.05790.02600.11800.002177239733137191218.1291111760.630.06400.00091.04500.02140.11790.001874325726117191119.124741480.500.06350.00101.03950.02250.11850.002072436724117221120.1401362460.550.06460.00101.05910.02200.11860.002076131733117231121.1371352280.590.06430.00091.04730.02710.11890.003975034728137242316T097 花岗片麻岩 30°48.297′N、90°0.404′E1.11093226970.460.06450.00071.07180.02300.12040.00237671874011733132.117443711210.390.06500.00091.07240.03200.11930.00287762974016726163.1872925300.550.06520.00101.07060.02060.11920.00217833173910726124.11555299280.570.06370.00061.06100.02040.12110.00257311473410737145.1641873900.480.06470.00071.06700.01680.11950.001576524737872896.1491452970.490.06420.00071.06840.02200.12070.00257502073811735157.124731470.50.06470.00101.07000.02800.11970.00247653373914729148.118261911040.560.06410.00071.06740.02490.12100.00287432373712736169.1762564710.540.06470.00081.06710.02370.12060.003376526737127341910.11466288720.720.06350.00071.06260.02450.12140.002672424735127391511.11266847320.930.06420.00071.06620.02450.12080.003074622737127351712.116350310270.490.06460.00091.06420.02180.12120.003976128736117382213.11184137380.560.06470.00061.07140.02370.12020.002776520739127321614.11074086340.640.06490.00051.07530.01930.12010.00227721774197311215.11304358020.540.06540.00081.06600.02960.11990.005178726737157302916.11114266690.640.06470.00071.07650.02240.12040.002376522742117331317.11214826960.690.06470.00071.07420.02110.11990.002076523741107301118.11384378210.530.06550.00241.06890.03040.12040.004179177738157332419.11194847100.680.06470.00081.06590.02560.12070.004776527737137352720.11012026550.310.06460.00091.07060.04080.12000.0044761297392073025

图5 纳木错西缘地区变质辉长岩(a)和花岗片麻岩(b)的锆石的U-Pb谐和图Fig.5 U-Pb zircon concordia diagrams of the dated meta-gabbro (a) and granitic gneiss (b) samples from the western Nam Tso area

图6 纳木错西缘地区变质辉长岩和花岗片麻岩的SiO2-Zr/Ti (a)、Th-Co (b)、FeOT/MgO-SiO2 (c)和Nb/Th-Nb/La (d)图解(底图据Winchester and Floyd, 1977; Hastie et al., 2007; Miyashiro, 1974; Li et al., 2006) 塞浦路斯斜长花岗岩的数据引自Freund et al. (2014)Fig.6 SiO2 vs. Zr/Ti (a), Th vs. Co (b), FeOT/MgO vs. SiO2 (c) and Nb/Th vs. Nb/La (d) diagrams of the meta-gabbro and granitic gneiss samples from the western Nam Tso area (base map after Winchester and Floyd, 1977; Hastie et al., 2007; Miyashiro, 1974; Li et al., 2006)

表2纳木错西缘新元古代中期岩浆岩的锆石Hf同位素组成

Table 2 Hf isotopic compositions for zircons of the Middle Neoproterozoic magmatic rocks in the western Nam Tso area

测点号年龄(Ma)176Yb177Hf2σ176Lu177Hf2σ176Hf177Hf2σ176Hf177HfiεHf(0)εHf(t)2σtDM (Ma)tCDM (Ma)fLu/Hf16T1601.17240.0061110.0001740.0003030.0000080.2825990.0000180.282595-6.19.70.69081023-0.992.17190.0038570.0000090.0001950.0000000.2825560.0000190.282553-7.78.10.79661121-0.993.17210.0033930.0000260.0001710.0000010.2825410.0000160.282539-8.27.70.69851152-0.994.17180.0037850.0000300.0001950.0000010.2825770.0000170.282574-6.98.80.69371074-0.995.17070.0038820.0000180.0001950.0000010.2825840.0000210.282581-6.78.80.79271066-0.996.17230.0043340.0000470.0002240.0000030.2825350.0000160.282532-8.47.50.69951165-0.997.17150.0050520.0000110.0002530.0000010.2825800.0000160.282576-6.88.90.69341071-0.998.17250.0028620.0000560.0001490.0000030.2825180.0000150.282516-9.06.90.510171201-1.009.17190.0047890.0000670.0002360.0000030.2825570.0000150.282554-7.68.20.59641118-0.9910.17270.0595990.0008750.0025760.0000270.2825010.0000210.282466-9.65.20.811101312-0.9211.17250.0416380.0006940.0016440.0000260.2825800.0000140.282558-6.88.40.59681106-0.9512.17160.0023640.0000110.0001200.0000000.2825600.0000160.282558-7.58.20.69581110-1.0013.17190.0041310.0000610.0002100.0000030.2825610.0000160.282558-7.58.30.69591109-0.9914.17160.0056830.0000380.0002810.0000020.2825400.0000200.282536-8.27.50.79901160-0.9915.17230.0101910.0008130.0005570.0000490.2825580.0000180.282551-7.68.10.79711123-0.9816T0971.17330.0366620.0003520.0015560.0000140.2823150.0000180.282294-16.2-0.70.613431691-0.953.17260.0355540.0005860.0015210.0000240.2822790.0000160.282258-17.4-2.20.613931775-0.954.17370.0462160.0005710.0019540.0000250.2822910.0000150.282264-17.0-1.70.513921754-0.945.17280.0344660.0005290.0014530.0000230.2823110.0000170.282291-16.3-1.00.613461701-0.967.17290.0357950.0009820.0014880.0000420.2821950.0000280.282307-15.9-0.40.713191665-0.968.17360.0279200.0001440.0011800.0000050.2823230.0000190.282219-18.8-3.31.114491856-0.959.17340.0381080.0002980.0015210.0000130.2822400.0000320.282278-16.6-1.30.613691725-0.9510.17390.0412480.0003800.0017690.0000150.2823030.0000160.282315-15.20.10.513201641-0.9411.17350.0482960.0001010.0020000.0000050.2823430.0000150.282271-17.1-1.50.513711741-0.9612.17380.0439360.0005980.0018490.0000220.2823140.0000140.282303-15.7-0.30.513351668-0.9413.17320.0304900.0001790.0012870.0000070.2822890.0000130.282324-15.10.30.512991624-0.9514.17310.0439350.0006200.0018400.0000260.2823290.0000140.282286-16.3-1.10.613571710-0.9515.17300.0361090.0003690.0015030.0000170.2823450.0000140.282319-15.20.10.613121638-0.9516.17330.0399000.0006430.0017220.0000290.2823100.0000160.282289-16.2-0.90.513561703-0.94

3 分析结果

3.1 锆石LA-ICP-MS U-Pb年代学

阴极荧光照片显示，变质辉长岩样品(16T160)中锆石粒度为60～120μm，多具有弱环带或补丁状结构(图4a)。21个LA-ICP-MS U-Pb分析点的Th含量为59×10-6～140×10-6，U含量为118×10-6～270×10-6，Th/U比值为0.22～0.76 (表1)，与典型岩浆成因锆石相似(吴元保和郑永飞, 2004)。所有分析点的U-Pb同位素年龄在误差范围内谐和，206Pb/238U年龄的变化范围为707～727Ma(表1)，206Pb/238U年龄加权平均值为720±6Ma(2σ；MSWD=0.10)(图5a)，相当于新元古代中期，代表了变质辉长岩的岩浆结晶年龄。

花岗片麻岩样品中锆石具有较典型的岩浆振荡环带结构(图4b)，且晶形比较完整，呈自形晶-半自形晶，长约50～150μm，长宽比为1:2～1:3，显示出岩浆锆石的特点。锆石测点的Th含量为73×10-6～684×10-6，U为147×10-6～1121×10-6，Th、U含量呈现出较好的正相关关系，Th/U比值介于0.31～0.93之间(表1)，为典型的岩浆锆石(吴元保和郑永飞, 2004)。所有测点在U-Pb谐和图(图5b)上集中落在谐和线上或其附近，获得732±7Ma(2σ；MSWD=0.06)的206Pb/238U加权平均年龄，代表花岗片麻岩的原岩形成年龄，即新元古代中期。

表3纳木错西缘新元古代中期岩浆岩的全岩地球化学分析结果(主量元素:wt%；稀土和微量元素:×10-6)

Table 3 Concentrations of major elements (wt%) and trace elements (×10-6) of the Middle Neoproterozoic magmatic rocks in the western Nam Tso area

样品号16T16017T08217T08317T08417T08516T09716T09816T09916T10016T10116T102岩性变质辉长岩花岗片麻岩SiO251.2152.7049.7949.1751.8275.3976.7476.2176.7176.1875.71TiO20.500.570.430.350.630.190.170.170.170.180.18Al2O316.4217.1115.4014.0914.2914.0613.7413.7313.9914.0613.95Fe2O3T7.235.826.446.6210.691.481.100.970.881.001.12MnO0.100.080.100.110.150.020.010.010.010.010.02MgO8.117.489.8510.807.470.540.460.840.110.140.54CaO11.7610.8113.8513.619.780.490.370.380.340.480.44Na2O2.953.912.101.413.666.995.946.406.675.986.75K2O0.600.630.521.420.590.410.770.530.490.890.36P2O50.050.060.040.030.050.060.030.050.040.060.06LOL1.221.111.312.040.810.761.200.890.931.030.93Total100.15100.2799.8299.6599.95100.39100.53100.18100.34100.01100.06Li4.784.016.578.434.301.811.284.021.903.230.14Be0.480.970.620.090.251.681.591.080.842.140.91Sc35.539.037.652.447.62.572.462.732.202.042.21V18018217317324810.911.810.610.711.614.9Cr50.083.828921233.31.752.532.041.791.801.78Co33.528.037.142.046.01.971.351.421.001.571.88Ni62.984.817515555.91.232.481.651.572.522.02Ga14.914.612.311.214.813.114.313.713.813.711.4Rb23.825.024.710818.719.838.333.119.841.818.1Sr15116212310414148.247.850.637.856.534.7Zr34.344.124.115.323.0151145130134159141Nb2.772.441.380.831.567.048.256.879.058.877.56Cs0.390.350.882.420.250.240.280.240.170.260.15Ba36.555.634.394.649.660.489.248.744.396.725.3Ta0.210.210.130.090.140.921.100.890.671.860.95Pb9.397.502.894.383.053.162.752.042.475.523.15Th0.920.850.370.250.3624.926.724.218.930.016.1U0.570.280.130.080.113.032.071.552.012.981.97Hf1.041.410.820.540.844.484.243.843.944.754.00Y11.010.69.048.1212.49.447.576.5910.106.877.58La3.143.782.301.452.1821.019.831.126.329.521.3Ce6.997.725.303.435.6242.039.646.843.649.130.9Pr0.940.980.710.500.843.813.394.714.855.113.32Nd4.404.433.402.574.2812.110.717.515.116.711.7Sm1.261.261.060.851.401.801.732.162.532.781.81Eu0.510.540.430.360.590.480.310.430.370.440.37Gd1.731.611.451.201.841.561.601.522.052.301.49Tb0.310.300.270.230.350.270.220.220.330.330.25Dy2.001.931.701.472.321.671.441.291.761.591.37Ho0.420.410.370.320.500.350.270.240.340.270.25Er1.281.231.120.931.491.060.830.691.030.740.75Tm0.200.180.150.130.210.180.130.110.160.100.11Yb1.311.201.010.881.441.120.850.710.950.660.65Lu0.200.180.150.130.220.180.150.120.160.120.13

图7 纳木错西缘地区变质辉长岩和花岗片麻岩的哈克图解Fig.7 Harker diagrams of the meta-gabbro and granitic gneiss samples from the western Nam Tso area

图8 纳木错西缘地区变质辉长岩(a、b)和花岗片麻岩(c、d)的球粒陨石标准化稀土元素配分图和原始地幔标准化微量元素蛛网图 N-MORB=正常洋中脊型玄武岩；E-MORB=富集型洋中脊玄武岩；IAB=岛弧玄武岩. 标准化值、N-MORB和E-MORB数据引自Sun and McDonough (1989)；IAB数据引自Pearce and Peate (1995)Fig.8 Chondrite-normalized rare-earth element patterns and primitive mantle-normalized spider diagrams for the meta-gabbro (a, b) and granitic gneiss (c, d) samples from the western Nam Tso area

图9 纳木错西缘地区变质辉长岩和花岗片麻岩的V-Ti (a)、Hf-Th-Ta (b)、Th/Yb-Nb/Yb (c)和Nb-Y (d)构造环境判别图解(底图据Li et al., 2015; Pearce and Peate, 1995)Fig.9 V vs. Ti (a), Hf-Th-Ta (b), Th/Yb vs. Nb/Yb (c) and Nb vs. Y (d) diagrams of the meta-gabbro and granitic gneiss samples from the western Nam Tso area (base map after Li et al., 2015; Pearce and Peate, 1995)

图10 北拉萨地块及其邻区的元古代-早古生代主要构造岩浆-变质事件的时空分布(据Hu et al., 2018c)Fig.10 Condensed timetable of major Proterozoic-Early Paleozoic geological events in the North Lhasa terrane and its adjacent regions in the Gondwana (after Hu et al., 2018c)

3.2 锆石Hf同位素

3.3 全岩地球化学

变质辉长岩和花岗片麻岩样品的主量和微量元素的分析结果见表3。镜下和野外观察结果表明变质辉长岩受到了后期变质作用的影响，样品中大离子亲石元素(例如K、Rb、Sr和Ba等)的含量可能会发生迁移。因此，本次研究主要依据高场强元素(例如Ti、Nb、Ta、Zr和Hf等)、稀土元素和过渡元素(例如V、Fe和Mn等)的含量来对样品进行岩石学分类和成因讨论。

变质辉长岩样品具有较均一的SiO2含量，变化范围集中在50.38%～53.15%之间(去烧失量归一化处理之后，下同)。在SiO2-Zr/Ti(图6a)和FeOT/MgO-SiO2(图6c)岩石分类图上，样品投点落入亚碱性钙碱性玄武岩区域。样品的稀土元素(REE)总体含量在14.5×10-6～25.8×10-6之间。样品具有变化的Mg#值(62～79)；Fe2O3T(图7b)、TiO2(图7c)和V(图7f)与Mg#呈现负相关关系，而Ni(图7d)和Cr(图7e)与Mg#具有正相关关系。在粒陨石标准化稀土元素球配分图(图8a)和原始地幔标准化微量元素蛛网图(图8b)上，所有样品的曲线一致性较好，均表现为平坦的曲线，兼具Nb、Ta和Ti负异常，未见明显的Eu异常(Eu/Eu*=1.06～1.16)。

所有花岗片麻岩样品均表现出高硅的特征(SiO2含量为75.67%～77.26%)。在Th-Co图解上，样品投点落在高钾钙碱性系列区(图6b)。样品的ΣREE总体含量较低，在223×10-6～283×10-6之间，轻稀土元素(LREE)相对重稀土元素(HREE)富集，LREE/HREE比值为1.73～2.66。在球粒陨石标准化稀土元素配分图(图8c)上，所有样品的曲线一致性较好，均表现为右倾的海鸥型，(La/Yb)N比值为13.4～32.1，同时具有明显负Eu异常，Eu/Eu*比值为0.50～0.88。在原始地幔标准化微量元素蛛网图(图8d)上，样品亏损Nb、Ta和Ti元素，富集Zr、Hf和Th元素。

4 讨论

4.1 岩石成因

4.1.1 变质辉长岩

变质辉长岩样品的成分类似于亚碱性钙碱性玄武岩(图6c)，并且亏损Nb、Ta和Ti元素(图8b)，表现出了岛弧亲缘性(Pearce and Peate, 1995)。在洋壳消减过程中，Th表现出活动性，而Nb表现出不活动性，因此，Th/Nb比值可以用来验证俯冲消减事件是否存在(Pearce, 2014)。在Th/Yb-Nb/Yb图解上(图9c)，变质辉长岩具有相对较高的Th/Nb比值，高于MORB-OIB地幔演化线；这一地球化学特征与岛弧玄武岩类似(Pearce, 2014)。虽然地壳混染过程也可以形成类似的地球化学特征，但是这些样品并没有在Nb/Th-Nb/La图解(图6d)上表现出相关的特征。在V-Ti(图9a)和Hf-Th-Ta(图9b)构造环境判别图解上，变质辉长岩样品投点均落在岛弧玄武岩上，进一步确定了这一岛弧亲缘性。考虑到变质辉长岩中锆石具有较高的εHf(t)值(+5.2～+9.7)，这些岩石应当是源自俯冲环境下相对亏损的地幔楔(Wuetal., 2006; Pearce and Peate, 1995)。

变质辉长岩样品具有变化的Mg#值(62～79)，表明它们可能经历一定程度的结晶分异作用。Fe2O3T(图7b)、TiO2(图7c)和V(图7f)与Mg#呈现负相关关系，指示岩浆演化晚期存在含Fe、Ti矿物的结晶分异。Ni(图7d)和Cr(图7e)与Mg#具有正相关关系，与橄榄石和/或辉石的分异一致。值得注意的是，Al2O3与Mg#没有表现出线性关系，表明长石的结晶分异作用不明显。

4.1.2 花岗片麻岩

花岗质岩石通常被分为I型、S型、M型和A型，前三种主要根据其源岩性质划分，A型花岗岩则是一类具有特殊的地球化学特征以及特定构造背景的花岗岩。M型花岗岩是洋壳的组成部分，一般具有低Th的特点，例如塞浦路斯蛇绿岩中的斜长花岗岩(Freundetal., 2014)，与本文研究的花岗片麻岩明显不同(图6b)。A型花岗岩普遍具有较高的10000Ga/Al比值(>2.7)和高场强元素含量(Zr+Nb+Ce+Y>250×10-6)，也与本文研究的花岗片麻岩(10000Ga/Al=1.53～1.95；Zr+Nb+Ce+Y=187×10-6～224×10-6)具有明显差别。实验研究表明，在准铝质到弱过铝质的I型花岗岩浆中，磷灰石的溶解度很低，磷灰石容易达到过饱和而结晶分离，在岩浆分异过程中随SiO2的增加而降低；而在强过铝质S型花岗岩浆中，磷灰石溶解度变化趋势与此相反，常具有高的P2O5含量(>0.26%)，并且随铝饱和指数的增加而增大(Wolf and London, 1994)。本文的数据显示，花岗片麻岩具有较低的P2O5含量(0.03×10-6～0.06×10-6)，并且P2O5含量随着SiO2的增加而降低(图7i)，所以应当属于I型花岗质岩石。

4.2 对北拉萨地块起源的约束

前人工作已在北拉萨地块上识别出了新元古代弧岩浆记录。张修政等(2013)报道北拉萨地块永珠地区的念青唐古拉岩群中存在近同时代的弧岩浆岩(约742 Ma)；Huetal. (2018b, c)在邻近地区报道了822～760Ma的弧后盆地岩浆岩。本文研究的变质辉长岩形成于新元古代中期(约720Ma)洋壳俯冲消减过程。同时，花岗片麻岩的原岩形成于732Ma左右，为I型花岗质岩石，并且在Nb-Y构造环境判别图解中投点落入同碰撞和火山弧花岗岩区(图9d)，与变质辉长岩的形成时代和构造环境相似。这些变质辉长岩和花岗片麻岩代表了北拉萨地块上最年轻的新元古代弧岩浆记录，从而为探索北拉萨地块的起源提供了新的约束。

近年来，大量碎屑锆石数据表明，南拉萨地块起源于印度大陆东北缘(Guoetal., 2017)。虽然有学者认为南、北拉萨地块在前寒武纪是连接在一起的(Yin and Harrison, 2000; Gehrelsetal., 2011; Guoetal., 2017)，但是南拉萨地块的前寒武纪林芝岩群与北拉萨地块的念青唐古拉岩群具有完全不一致的物质组成(图10)。林芝岩群中包含1780Ma、1343～1276Ma和1250Ma的花岗质岩石，分别与哥伦比亚超大陆边缘岩浆弧、哥伦比亚超大陆裂解和罗迪尼亚超大陆聚合相关(Xuetal., 2013; Linetal., 2013)。这些花岗质岩石在1117Ma和618Ma左右发生了变质，可能与印度-澳大利亚格林威尔期造山事件和新元古代晚期Kuunga造山带相关(Linetal., 2013)。与此形成鲜明对比的是，念青唐古拉岩群中没有中元古代的岩浆或者变质记录。念青唐古拉岩群主要包含新元古代早期N-MORB型基性岩(925～850Ma; 胡培远等, 2016; Zhangetal., 2012a; Dongetal., 2011)、新元古代中期岛弧岩浆岩和变质岩(822～660Ma; 张修政等, 2013; Dongetal., 2011; Huetal., 2018b, c)和新元古代中期碰撞型高压麻粒岩(650 Ma; Zhangetal., 2012a)。此外，北拉萨地块上存在大规模的埃迪卡拉纪-早古生代安第斯型岩浆弧(Zhuetal., 2012; Huetal., 2013, 2018a; Dingetal., 2015)，而南拉萨地块上并没有可以与之对比的岩浆记录。因此，南、北拉萨地块应该具有不同的前寒武纪演化历史。

如果北拉萨地块不是起源于印度大陆东北缘，目前有另外两个假说：(1)澳大利亚大陆北缘(Ferrarietal., 2008; Zhuetal., 2011)和(2)东非造山带北段(Zhangetal., 2012a; Huetal., 2018a)。东非造山带记录了莫桑比克洋的漫长演化历史(从1080Ma以前到600Ma左右; Merdithetal., 2017; Moleetal., 2018)。北拉萨地块上的新元古代早期N-MORB型基性岩与莫桑比克洋的洋壳时代一致(Fritzetal., 2013; Huetal., 2018d)。新元古代中期岛弧岩浆岩和变质岩(张修政等, 2013; Dongetal., 2011; Huetal., 2018b, c, 2019; 本文报道)也与东非造山带中的岩石可以对比(图10)。然而，澳大利亚大陆缺少这些岩浆和变质记录，所以北拉萨地块应当起源于东非造山带北段。

虽然安多微陆块位于班公湖-怒江缝合带内部，但是传统观点仍然将其视作北拉萨地块的一部分(Xuetal., 1985)。近年来，有学者对这一观点提出了质疑。Guynnetal. (2013)指出安多微陆块可能与南羌塘地块更具有亲缘性。解超明等(2014)提出安多微陆块可能来源于扬子地块。安多片麻岩的主要组成岩石为新元古代早-中期(920～800Ma; Guynnetal., 2012; Zhangetal., 2012b; 解超明等, 2014;辜平阳等, 2012)和早古生代(532～488Ma; Guynnetal., 2012; Xieetal., 2013)的岩浆岩。这些早古生代岩浆岩的成因争议较小，通常认为其与冈瓦纳大陆北缘的早古生代安第斯型岩浆弧相关(Guynnetal., 2012; Xieetal., 2013)。新元古代早-中期岩浆岩的成因争议较大。解超明等(2014)和辜平阳等(2012)认为安多微陆块上的新元古代早-中期玄武质岩石和A型花岗岩形成于大陆裂谷环境，与罗迪尼亚超大陆的裂解相关。然而，Zhangetal. (2012b)在安多微陆块上识别出了同时代的弧岩浆岩。Huetal. (2018b)指出弧后盆地环境可以解释同时代岛弧型和裂解型岩浆岩的共存现象。如图10所示，安多微陆块的前寒武纪-早古生代岩浆记录与南羌塘和扬子地块区别明显，但是与北拉萨地块可以对比。因此，本文倾向于认为安多微陆块与北拉萨地块具有类似的前寒武纪演化历史。

5 结论

综合上述分析讨论，得出以下初步结论：

(1)北拉萨地块念青唐古拉岩群中变质辉长岩和花岗片麻岩的原岩形成时代分别为720±6 Ma和732±7Ma，相当于新元古代中期。

(2)变质辉长岩为亚碱性钙碱性系列，具有Nb、Ta和Ti负异常，与岛弧玄武岩类似。变质辉长岩中锆石具有较高的εHf(t)值(+5.2～+9.7)，应当是源自俯冲环境下相对亏损的地幔楔。花岗片麻岩为I型花岗质岩石，并且具有较为均一的锆石εHf(t)值(-3.3～+0.3)，可能形成于地壳内古元古代变质火成岩的部分熔融作用。

(3)念青唐古拉岩群中的前寒武纪岩浆和变质记录与东非造山带的活动时限较为一致，因而北拉萨地块可能与东非造山带具有亲缘性。

致谢感谢张拴宏研究员和另外一位审稿人对本文提出的中肯的、建设性的修改意见。

谨以此文祝贺肖序常院士九十华诞，感谢肖先生多年来对我们的教导、关心和大力支持，敬祝肖先生生日快乐、健康长寿！