西藏拉萨地块北部晚白垩世晚期基性岩墙的成因:来自锆石U-Pb年代学及地球化学的制约*
2020-04-01于云鹏王明解超明范建军罗安波曾孝文
于云鹏 王明, 2** 解超明, 2 范建军 罗安波 曾孝文
1. 吉林大学地球科学学院,长春 1300612. 东北亚矿产资源评价自然资源部重点实验室,长春 130061
班公湖-怒江缝合带(后简称班-怒带)横贯青藏高原中部,其北部为羌塘地块,南部为拉萨地块。作为青藏高原众多缝合带中重要的一条,班-怒带的性质、大洋俯冲极性、演化过程与闭合机制过去一直受到众多学者的研究(Yin and Harrison, 2000; Kappetal., 2003; 潘桂棠等, 2006; 张玉修, 2007)。目前主流的观点根据沉积及岩浆的记录,认为羌塘地块与拉萨地块于早白垩世期间沿班-怒带发生碰撞(Kappetal., 2007; Suietal., 2013; Zhuetal., 2016),随后的持续碰撞导致缝合带附近的地壳在晚白垩世发生了明显的增厚缩短(Guynnetal., 2006; Kappetal., 2007)。近年来在大规模区域地质调查的基础上,大量岩石年代学、岩石地球化学以及同位素数据的应用,使得拉萨地块北部的晚白垩世岩浆岩演化过程得到明确。拉萨地块中北部94~79Ma的埃达克岩及花岗质岩石为碰撞后加厚下地壳部分熔融及拆沉的产物(Zhaoetal., 2008; 余红霞等, 2011; 王保弟等, 2013; 李华亮等, 2014; 雷鸣等, 2015; 雷鸣, 2016; 刘涵等, 2015; Sunetal., 2015; Yietal., 2018; Liuetal., 2019),80Ma左右的花岗质岩石则被解释为后碰撞环境下的岩浆作用产物(高顺宝等, 2011; 马蓁, 2013),而在班公湖地区发现的88Ma辉绿岩脉、79Ma花岗斑岩脉以及措勤地区87Ma双峰式火山岩被认为是伸展环境下的产物(辛洪波和曲晓明, 2006; 江军华等, 2011; 张硕等, 2014)。这些晚白垩世岩浆活动不仅记录了班公湖-怒江洋(后简称班-怒洋)闭合后的碰撞过程:碰撞-地壳加厚-拆沉(伴随伸展),同时也是认识碰撞后地壳的生长与组成的重要突破口。但由于以往对碰撞后岩浆岩的研究主要集中在晚白垩世早期,而对于晚白垩世晚期(<80Ma)的基性岩浆事件鲜有报道(Zhuetal., 2016)。因此本文在以往研究的基础上,结合拉萨地块北部晚白垩世晚期基性岩墙的锆石U-Pb年代学及岩石地球化学数据,以期进一步确定拉萨地块北部碰撞后的地壳演化过程。
1 地质概况及样品描述
拉萨地块位于班公湖-怒江缝合带和印度-雅鲁藏布江缝合带之间,并被狮泉河-纳木错蛇绿混杂岩和洛巴堆-米拉山断裂带由南至北分为北、中、南三个块体(图1a)。冈底斯岩浆带位于拉萨地块内部,其主体为白垩纪和古近纪的岩浆活动产物(莫宣学等, 2005; Harrisonetal., 2000; Jietal., 2012; Zhuetal., 2011; Zhangetal., 2010; Suietal., 2013)。研究区位于北拉萨地块尼玛县阿索乡一带,在阿索乡以南50km的目思地区出露有一系列近NE-SW走向的晚白垩世基性岩墙(图1b)。
目思地区出露的地层主要为晚古生代地层,包括中二叠统昂杰组及上石炭统-下二叠统拉嘎组(图1b)。其中昂杰组主要以杂砂岩、变质砂岩为主,拉嘎组以含生物碎屑灰岩、泥晶灰岩为主,反映地层的形成环境为滨浅海相。为了进一步明确该区的晚白垩世岩浆事件,并探究基性岩墙的成因,对上述岩体进行样品采集并对其中一个岩体进行剖面测制。剖面的详细位置见图1b,剖面起始坐标为北纬31°39′36″、东经85°57′00″。
在目思地区测制的基性岩墙剖面可反映岩墙的野外产状及围岩性质(图2)。岩墙侵出产状近直立,宽约6~15m,长度约130m,其围岩主要为上石炭统-下二叠统拉嘎组(图3a),包括二段变质砂岩、变质长石石英砂岩、变质复成分砾岩及含砾板岩。岩墙岩性为典型的辉绿玢岩,风化程度较强烈,部分岩石的斑晶被剥蚀形成孔洞且辉石矿物破碎严重(图3b, d)。岩石组构为似斑状结构,块状构造。斑晶的主要成分为斜长石,其次为单斜辉石。镜下明显可见岩石发生弱蚀变,板柱状斜长石斑晶发育聚片双晶,长度可达5mm(图3c)。岩石基质由大量斜长石、少量单斜辉石及微量的角闪石、黑云母微晶组成(图3d),粒径在0.1~0.7mm之间。副矿物为少量锆石及磷灰石。
2 分析方法
在目思基性岩墙中共采集9件辉绿玢岩地球化学样品用于全岩地球化学测试分析(表1)。为了避免弱蚀变对样品全岩地球化学分析的影响即确保实验结果的可信程度,样品在野外采集过程中均挑选最新鲜岩石并去除风化面。之后用<5% HNO3和去离子水清洗样品,待样品干燥后将样品放入无污染玛瑙球磨机粉碎至200目,以上实验过程均在河北省廊坊市区域地质矿产调查研究院无污染实验室完成。所有粉碎的样品送到中国地质大学(北京)地学实验中心用于全岩地球化学主微量元素分析。全岩地球化学主量元素测试在等离子体发射光谱仪(ICP-OEC)实验中进行,测试仪器为PS-950等离子体光谱仪;微量元素的化学预处理采用两酸(HNO3+HF)高压反应釜溶样法,测试方法为ICP-MS法,测试仪器为Agilent-7500a电感耦合等离子质谱仪,实验过程中采用国际标准参考样品AGV-2和GSR-3进行校对(Govindaraju, 1994)。
图1 西藏南部冈底斯岩浆带白垩世火山岩分布(a,据Zhu et al., 2011修改)和目思地区地质简图(b)年龄数据来自:Zhao et al., 2008; 余红霞等, 2011; 刘涵等, 2015; 雷鸣, 2016; Sun et al., 2015; Yi et al., 2018; Liu et al., 2019Fig.1 Distribution of the Cretaceous volcanic rocks on Gangdese belt, South Tibet (a, modified after Zhu et al., 2011) and simplified geological map of the Musi area (b)
图2 西藏尼玛县阿索乡北部目思基性岩墙剖面1-变质砂岩;2-变质长石石英砂岩;3-变质复成分砾岩;4-含砾板岩;5-辉绿玢岩. C2P1l2-上石炭统-下二叠统拉嘎组二段;K2βμ-晚白垩世辉绿玢岩;N17T53-基性岩墙年龄及化学样品采样点Fig.2 Section of Musi mafic dyke in north of Asa Village, Nyima County, Tibet
此外在目思基性岩墙中共采集3件年龄样品,采样坐标见图1b及表2。定年样品的U、Th和Pb同位素元素分析在中国地质大学(北京)的实验室完成,分析仪器为193nm激光剥蚀进样系统UP 193SS和Agilent 7500a型四级杆等离子体质谱仪构成的激光等离子体质谱仪(LA-ICP-MS)。激光束直径为36μm,以He作为剥蚀物质载气。国际标样91500作为外标校对同位素分馏,U、Th含量采用NIST610中的29Si作为内标校正仪器引起的信号漂移。详细的方法和分析流程见Yuanetal. (2004)。锆石定年数据处理采用Isoplot (version 3.0)和Glitter (version 4.4)软件(Ludwig, 2003)。
3 分析结果
3.1 全岩地球化学
表1目思基性岩墙主量(wt%)和微量(×10-6)元素分析结果
Table 1 Analytic results of major (wt%) and trace (×10-6) elements of the Musi mafic dyke
样品号T51H1T51H2T51H3T53H1T53H2T53H3T56H1T56H2T56H3SiO253.5352.4753.5552.8551.9053.1852.5053.0652.34TiO21.651.631.681.571.531.451.611.621.74Al2O314.4414.7814.5615.6516.2316.1015.8014.8615.26Fe2OT39.199.148.818.878.828.348.438.509.28MnO0.150.150.150.120.120.120.120.130.15MgO4.574.754.304.824.974.854.694.343.98CaO7.848.198.107.427.637.297.697.087.60Na2O2.732.882.943.063.113.133.113.303.34K2O0.950.980.971.131.011.001.141.191.25P2O50.500.500.510.520.500.460.540.540.58LOl4.254.434.263.844.133.914.295.164.38Total99.7999.9199.8399.8699.9599.8399.9099.7899.90Mg#53.6654.7553.1955.8856.7557.5256.4654.3349.99Li40.1044.7939.0641.0736.1135.3438.5528.4628.25Sc22.8824.3623.3422.6321.5221.0221.9022.5023.48V191.3190.9198.2183.3174.0174.4182.3189.8200.2Cr58.3658.6457.3668.7879.4875.6269.8651.3055.00Co41.7435.4629.8634.3632.8631.5634.1432.1434.84Ni69.6468.4066.6465.7374.1770.2669.7055.9458.54Cu24.6225.3224.6622.9124.2823.1222.4221.5222.52Zn79.0880.0086.9677.5678.3675.4088.5087.8489.78Ga20.3520.6920.2520.4720.8120.5220.8820.6921.47Rb18.7219.4721.8227.5924.2224.7426.5631.7234.12Sr629.2643.4610.2605.3580.6585.6561.2552.6563.2Y33.3433.6633.6431.4131.1430.6832.0233.3635.32Zr290.4287.2295.6267.5266.8260.0273.4284.0293.6Nb16.6016.5216.4515.5915.2614.9215.8716.4416.97Cs2.562.614.664.454.725.024.666.807.62Ba319.8328.0314.6345.6344.2333.2349.8387.4465.0La41.4440.9240.7439.7438.0437.5039.4640.9642.66Ce90.5489.1093.6486.8682.3281.4285.6288.1294.24Pr11.4511.6011.5310.9510.5010.3510.9211.2911.94Nd41.5841.9641.8439.3238.1837.5239.5241.1243.72Sm8.108.238.177.707.517.417.848.128.61Eu2.082.092.071.951.951.901.982.052.17Gd7.177.197.296.736.666.576.887.217.62Tb1.081.081.091.011.000.991.031.081.14Dy6.256.176.225.735.715.675.956.166.51Ho1.261.241.261.161.171.141.201.241.31Er3.453.383.453.143.183.153.273.393.55Tm0.500.490.500.460.460.460.480.500.52Yb3.082.993.042.772.802.782.892.993.12Lu0.470.450.460.420.420.420.440.450.46Hf6.986.646.796.096.316.266.536.776.88Ta1.020.980.970.920.930.950.910.991.00Pb6.215.445.987.498.368.149.069.139.39Th5.865.505.715.305.625.585.976.126.11U0.810.750.780.950.960.951.071.091.07Eu∗0.830.830.820.830.840.830.830.820.82(La/Yb)N9.649.829.6110.309.749.679.799.819.82(Ce/Yb)N8.158.288.568.728.168.138.238.188.40∑REE218.5216.9221.3207.9199.9197.3207.5214.7227.6
注:Mg#=100×Mg/(Mg+Fe);Eu*=EuN/(SmN×GdN)1/2;EuN、SmN、GdN、(La/Yb)N、(Ce/Yb)N均为球粒陨石标准化后的比值
表2目思基性岩墙锆石U-Pb同位素测试结果
Table 2 Zircon U-Pb isotope analytic results of the Musi mafic dyke
测点号含量(×10-6)ThUTh/UU-Pb同位素比值年龄(Ma)207Pb/206Pb1σ207Pb/235U1σ206Pb/238U1σ207Pb/235U1σ206Pb/238U1σ谐和度N17T51:E85°57′35″、N31°40′16″01416.49423.410.980.04750.00130.07460.00210.01140.0001732730.90.00%0257.35109.890.520.04760.00320.07480.00490.01140.00027357310.00%0625.4435.130.720.04730.00950.07800.01560.01200.00037615772-1.30 %07100.7114.230.880.05350.00340.08750.00550.01190.000285576111.84%10162.87220.670.730.04730.00230.07710.00360.01180.0002753761-1.32%1240.5259.060.690.04610.00420.07920.00720.01250.0002777801-3.75%15114.85163.390.700.04750.00250.07150.00380.01090.000270469.910.14%16269.61281.060.960.05000.00190.07760.00300.01130.000176372.10.95.41%20158.62219.390.720.04800.00220.07590.00360.01150.000274373.510.68%N17T53:E85°57′35″、N31°40′14″0137.6744.360.850.04730.00710.07890.01180.01210.00027711781-1.28%02164.74176.920.930.04750.00240.07610.00380.01160.00027447410.00%0326.8246.850.570.04650.00630.07950.01070.01240.00027810801-2.50%0453.8271.320.750.04730.00490.07830.00810.01200.00027787710.00%05526.31363.651.450.04690.00150.07570.00240.01170.000274274.91-1.20%0637.1946.900.790.04710.00720.07690.01180.01180.00027511761-1.32%0799.41129.760.770.04750.00270.07530.00430.01150.00027447410.00%0979.3989.370.890.04860.00350.08170.00580.01210.00028057812.56%1019.4228.040.690.04710.01190.07860.01980.01210.00037719782-1.28%12191.16196.160.970.04750.00230.07550.00370.01150.00027437410.00%1351.6363.470.810.04650.00550.07600.00890.01190.0002748761-2.63%1434.0043.350.780.04760.00750.08110.01270.01230.000279127920.00%1566.0793.140.710.04760.00360.07990.00590.01220.00027867810.00%1624.9233.760.740.04750.00900.07870.01490.01200.000377147720.00%1784.16112.890.750.04750.00350.07820.00570.01190.00027657610.00%18154.35158.890.970.04750.00280.07430.00440.01130.00027347310.00%19151.65184.390.820.04760.00260.07590.00410.01160.00027447410.00%20286.18247.991.150.05220.00190.08170.00290.01140.000280372.819.89%N17T56:E85°57′35″、N31°40′13″01117.24115.431.020.04790.00290.08380.00490.01270.00028258111.23%02140.03152.190.920.04780.00290.07100.00430.01080.000270469.111.30%0631.5253.200.590.05110.00690.08230.01090.01170.000280107526.67%0739.9450.870.790.04670.00670.07990.01140.01240.00027811791-1.27%0855.8264.790.860.05310.00720.08860.01190.01210.0002861178110.26%1175.2889.340.840.04750.00380.08070.00650.01230.00027967910.00%1266.5281.050.820.04650.00350.07780.00570.01210.0002765781-2.56%14122.32136.600.900.05010.00290.07990.00460.01160.00027847415.41%
图3 目思基性岩墙野外产状及岩相学特征(a)目思基性岩墙;(b)岩石近景照片;(c、d)目思基性岩墙岩相学特征(样品N17T51、N17T56,正交偏光). Pl-斜长石;Cpx-单斜辉石Fig.3 Occurrences and petrographic features of the Musi mafic dyke
图4 目思基性岩墙Nb/Y-Zr/TiO2×0.0001图解(a,底图据Winchester and Floyd, 1977)和AFM判别图解(b,底图据Irvine and Baragar, 1971)Fig.4 Nb/Y vs. Zr/TiO2×0.0001 (a, base map after Winchester and Floyd, 1977) and AFM (b, base map after Irvine and Baragar, 1971) discriminant diagrams for the Musi mafic dyke
在球粒陨石化稀土元素配分曲线图中,目思基性岩墙呈现出轻稀土元素富集的右倾曲线(图5a),(La/Yb)N(9.60~10.46)、(Ce/Yb)N(7.85~8.79)同样呈现出较高的比值。在原始地幔标准化微量元素蛛网图中(图5b),目思基性岩墙呈现出明显的大离子亲石元素富集(如Ba、Th、U和K),亏损Nb、Ta、Ti等高场强元素,并具有Pb的正异常。同时,Eu具有明显的负异常(Eu*=0.81~0.84)。
3.2 锆石U-Pb年龄
目思基性岩墙样品中获取的锆石的形态并不规则,在阴极发光图像中(图6),锆石多呈半自形长柱状,部分呈他形粒状,粒径范围100~300μm,长宽比1:1~3:1。锆石晶面较光滑,弱环带结构,Th/U比值为0.52~1.45,均大于0.3,表明锆石为岩浆成因。
35个测点206Pb/238U年龄值变化在69~81Ma,锆石年龄较为分散,存在72~76Ma和76~80Ma两个峰值年龄(图7a)。对两组年龄分别进行投图,在207Pb/235U-206Pb/238U谐和曲线图中,两组年龄样品的测点均分布于谐和线上(图7b, c),加权平均年龄分别为73.5±0.6Ma(MSWD=0.62, N=13)和77.3±0.5Ma(MSWD=1.1, N=17)。其中较为年轻的年龄73.5±0.6Ma最有可能代表岩体的形成时代,指示目思基性岩墙的岩浆结晶年龄为晚白垩世。
图5 目思基性岩墙球粒陨石标准化稀土元素配分曲线图(a)和原始地幔标准化微量元素蛛网图(b)(标准化值据Sun and McDonough, 1989)Fig.5 Chondrite-normalized REE patterns (a) and primitive mantle normalized muti-element diagrams (b) for the Musi mafic dyke (normalization values after Sun and McDonough, 1989)
图6 目思基性岩墙代表性锆石阴极发光(CL)图像Fig.6 Representative cathodoluminescence (CL) images of zircons from the Musi mafic dyke
4 讨论
4.1 岩石成因
在讨论岩浆来源之前,应首先考虑岩浆在上升过程中受地壳混染程度。相较于活跃的大离子亲石元素,稳定富集的高场强元素可以用于讨论岩浆是否与地壳混染(Winchester and Floyd, 1977)。Nb/Yb-Th/Yb图解和(La/Yb)PM-(Th/Nb)PM图解能够有效判别地壳物质的贡献(Xuetal., 2002)。目思岩墙的岩石样品全部落入岩浆混染区域中(图8a),而在(La/Yb)PM-(Th/Nb)PM图解中同样呈现出受大陆地壳混染的趋势(图8b)。岩石中的La/Ta和La/Sm受岩浆混合或地壳混染的影响较大,受地壳混染后的玄武岩La/Ta会在25以上,La/Sm会在5以上(王明等, 2010)。目思基性岩墙La/Ta介于39.46~43.40之间,La/Sm介于4.95~5.16之间。结合Nb/Yb-Th/Yb图解和(La/Yb)PM-(Th/Nb)PM图解的投图结果,认为岩浆在上升过程中受到了地壳物质的混染。
图7 目思基性岩墙锆石U-Pb谐和图Fig.7 U-Pb concordia diagrams of zircons from the Musi mafic dyke
现有研究表明,玄武质岩浆一般来源于地幔橄榄岩的部分熔融。由尖晶石二辉橄榄岩部分熔融的玄武质岩浆通常具有平缓的稀土元素分布模式,而由石榴子石二辉橄榄岩低程度熔融生成的岩浆,通常重稀土元素显著分馏(Hart and Dunn, 1993)。目思基性岩墙具有轻微的重稀土元素分馏(图5a),同时La/Sm-Sm/Yb的投图结果显示,初始岩浆来自于尖晶石+石榴石二辉橄榄岩地幔源区1%~5%的部分熔融,尖晶石与石榴石比值接近1:1(图9a)。在Ta/Yb-Th/Yb图解(Pearce, 2008)中,目思基性岩墙落入钙碱性玄武岩区域中,同样落入地幔序列以外(图9b)。结合目思基性岩墙中富集Ba、Th、U、K等大离子亲石元素,亏损Nb、Ta、Ti等高场强元素的地球化学特征,推断目思基性岩墙来源于富集地幔部分熔融。
目思基性岩墙中的Mg#值以及Cr、Ni质量分数(Mg#=50.0~57.5、Cr=51.30×10-6~79.48×10-6、Ni=55.94×10-6~74.168×10-6)均低于原生玄武质岩浆(Mg#=68~75、Cr=300×10-6~500×10-6、Ni=300×10-6~400×10-6; Freyetal., 1978; 朱弟成等, 2005),表明其经历了一定程度的结晶分异作用。在La/Yb-La和La-La/Sm中同样呈现出分离结晶的趋势(图10a, b)。而在Hacker图解中,Mg#与Cr、Ni呈现出明显的正相关(图10c, d),表明岩浆在演化过程中发生了橄榄石、单斜辉石及角闪石等镁铁质矿物的结晶分异作用。同时目思基性岩墙具有Eu负异常(Eu*=0.81~0.84)的特征表明岩石还发生了斜长石的分离结晶。
图8 目思基性岩墙Nb/Yb-Th/Yb图解(a,底图据Pearce, 2008)和 (La/Nb)PM-(Th/Nb)PM图解(b,底图据Frey et al., 2002)原始地幔标准化数据引自Sun and McDonough (1989),下地壳与中地壳数据引自Rudnick and Gao (2003)Fig.8 Nb/Yb vs. Th/Yb (a, base map after Pearce, 2008) and (La/Nb)PM vs. (Th/Nb)PM (b, base map after Frey et al., 2002) plots of Musi mafic dyke (base map after Xu et al., 2002)
图9 目思基性岩墙La/Sm-Sm/Yb图解(a,底图据Aldanmaz et al., 2000)和Ta/Yb-Th/Yb图解(b,底图据Pearce, 2008)Fig.9 La/Sm vs. Sm/Yb (a, base map after Aldanmaz et al., 2000) and Ta/Yb vs. Th/Yb (b, base map after Pearce, 2008) plots of Musi mafic dyke
图10 目思基性岩墙La/Yb-La图解(a,底图据Chung et al., 2009)、La-La/Sm图解(b,底图据陈玲等, 2013)、Mg#与Cr、Ni元素相关图(c、d)Fig.10 La/Yb vs. La plot (a, base map after Chung et al., 2009), La vs. La/Sm plot (b, base map after Chen et al., 2013), correlation diagram of Mg# and Cr, Ni (c, d) of Musi mafic dyke
综上所述,目思基性岩墙的源区为地幔尖晶石+石榴石二辉橄榄岩,在岩浆上升过程中遭受地壳物质混染。在岩浆演化的过程中,发生了明显的橄榄石、单斜辉石及斜长石矿物的结晶分异。
4.2 构造环境
以往的研究认为羌塘地块与拉萨地块于早白垩世沿闭合的班-怒带发生碰撞(Kappetal., 2007; Suietal., 2013; Zhuetal., 2016)。伴随着班-怒洋的洋壳陆下断离(113Ma),拉萨地块中北部与羌塘-拉萨碰撞带进入陆内碰撞环境(Zhuetal., 2009; Zhuetal., 2011)。板片的断离引起软流圈地幔的上涌,于莫霍面附近形成广泛分布的镁铁质岩浆,从而增加了地壳的厚度(Zhuetal., 2017)。Kappetal. (2007)根据尼玛盆地中地层的沉积特征,认为拉萨地块中北部于125~95Ma发生了明显的地壳增厚缩短,竟柱山组磨拉石建造与下覆海相地层的不整合进一步表明地壳最有可能在113~93Ma发生南北向的增厚缩短。而碰撞后末期的岩石圈增厚正是发生岩石圈拆沉作用的必要条件(Bonin, 2004)。随着拉萨中北部地壳的两次增厚,下地壳的深度与压力逐渐增大,成分逐渐由角闪岩相向榴辉岩相转变,密度逐渐增大(Wangetal., 2014)。张旗等(2006)认为该阶段的榴辉岩下地壳的密度小于下伏地幔,拆沉作用的发生仍需要大规模的岩浆活动使低密度的中酸性岩浆岩移出下地壳。而区域内大规模的同期(94~79Ma)埃达克质岩浆活动为拆沉模式提供了强有力的支持(Zhaoetal., 2008; 余红霞等, 2011; 王保弟等, 2013; 李华亮等, 2014; 雷鸣等, 2015; 刘涵等, 2015; Sunetal., 2015; Chenetal., 2015; Yietal., 2018; Liuetal., 2019)。
考虑到下地壳的拆沉作用受到空间及密度的制约,软流圈只能经受有限的岩浆上涌与部分熔融,大部分的拆沉模式并非下地壳+岩石圈地幔的大规模整体拆沉,而是下地壳如同冰块融化般的小规模逐渐拆沉(Houseman and Molnar, 1997; 张旗等, 2006)。区域内94~79Ma小规模岩浆活动埃达克质和富镁的地球化学特征全部与岩石圈拆沉有关(Zhaoetal., 2008; 余红霞等, 2011; 雷鸣等, 2015; 雷鸣, 2016; 刘涵等, 2015; Sunetal., 2015; Yietal., 2018; Liuetal., 2019),这些独特的地球化学特征恰好反映了小规模岩石圈拆沉的结果。随着下地壳拆沉作用的进行,地壳逐渐减薄,部分地区由于地壳过薄发生垮塌,形成伸展环境,于中北拉萨地块内部形成了与拆沉环境有关岩浆活动(94~79Ma,图1a)同时期的伸展背景下的岩浆活动(88~79Ma; 辛洪波和曲晓明, 2006; 江军华等, 2011; 张硕等, 2014)。本次研究于拉萨地块北部发现的基性岩墙具有74Ma的加权平均年龄,推测拉萨中北部的小规模拆沉作用持续到晚白垩世晚期。同时,由于小规模的拆沉模式导致基性岩浆的不发育(DeCellesetal., 2015),镁铁质岩浆在产生过程中被同期大量的花岗质岩浆所消耗(Wangetal., 2014)。本次所研究的对象岩性更偏中性,同样是因为岩石圈地幔产生的镁铁质岩浆在上升过程中与下地壳物质发生混染(图8)。
目思基性岩墙的构造成因可以总结如下:班-怒洋闭合后的后碰撞作用使拉萨地块增厚缩短,致使下地壳压力增大,形成榴辉岩相下地壳。下地壳发生岩浆上涌,拉萨地块中北部形成大规模埃达克质的岩浆活动,下地壳密度进一步加大。由于下地壳榴辉岩密度大于下伏橄榄岩地幔,重力的不稳定性引起下地壳发生拆沉。拆沉作用导致岩石圈地幔上涌,相关的地球化学特征显示幔源岩浆于源区内发生分离结晶,上升过程中与下地壳物质发生混染。与区域内同期拆沉成因岩浆岩的对比表明拆沉模式并非传统认为的下地壳+地幔整体拆沉,而是小规模的部分逐渐拆沉,该期拆沉作用沿羌塘-拉萨碰撞带持续至晚白垩世晚期。
5 结论
基于拉萨地块北部阿索地区目思基性岩墙的岩石学、全岩地球化学以及U-Pb年代学,可以得到以下几点认识:
(1)目思基性岩墙中最年轻的一组锆石U-Pb谐和年龄为74Ma,表明其形成时代为晚白垩世晚期;
(2)目思基性岩墙具有较高的SiO2、Mg#、Cr和Ni含量。这些地球化学特征表明岩浆来源于部分熔融的岩石圈地幔,并在上升过程中与拆沉的下地壳发生混染。
(3)目思基性岩墙形成于班-怒洋闭合后的后碰撞拆沉环境中,结合目思基性岩墙的锆石U-Pb年龄表明羌塘-拉萨碰撞带的拆沉作用一直持续到晚白垩世晚期。
致谢分析测试工作得到中国地质大学(北京)地学实验中心苏犁教授与张红雨老师的指导;两位审稿人提出了建设性的意见;在此一并致谢。