魯洪濤,黃維平,劉振宇,鄭 明,李海峰
1. 中國地質(zhì)科學院礦產(chǎn)資源研究所 成礦作用與資源評價重點實驗室,北京 100037;2. 吉林省地質(zhì)調(diào)查院,吉林 長春 130102;3. 中國地質(zhì)大學(北京)地球科學與資源學院,北京 100083
班公湖—怒江縫合帶(以下簡稱班怒帶)是拉薩地體與羌塘地體的分界線(圖 1a),班怒帶兩側(cè)分布的大量中生代巖漿巖常被用來研究還原班怒帶構(gòu)造演化過程[1-2]。雪如巖體地理位置位于西藏班戈縣以南20 km,大地構(gòu)造位置位于班怒帶中段南側(cè)班戈巖漿弧中部(圖 1a)。目前對雪如巖體的研究工作存在兩點不足:一是對雪如巖體的巖石成因及構(gòu)造環(huán)境的判別尚存在爭議,二是對于雪如巖體的成巖條件(溫壓、氧逸度)缺乏研究。本文對雪如巖體開展了地質(zhì)年代學、鋯石微量元素及全巖地球化學的研究工作,旨在查明雪如巖體的形成時代、巖石成因、成巖條件及構(gòu)造背景,進而約束班怒帶中段的演化過程。
圖1 西藏班戈地區(qū)雪如巖體地質(zhì)簡圖(圖a據(jù)[3])Fig.1 Geological sketch map of the Xueru Rocks in Bange area,Tibet
班公湖—怒江縫合帶位于青藏高原中部,是班公湖—怒江洋洋盆閉合過程中北拉薩地體與南羌塘地體碰撞拼合的產(chǎn)物[3]??p合帶中地層主要出露有侏羅系木嘎崗日群(JM)、上侏羅統(tǒng)—下白堊統(tǒng)沙木羅組(J3K1s)、上白堊統(tǒng)竟柱山組(K2j)等。木嘎崗日群巖性主要為玄武巖、硅質(zhì)巖、灰?guī)r等,沙木羅組巖性主要為灰?guī)r,竟柱山組巖性主要為粗碎屑巖??p合帶內(nèi)出露的巖漿巖主要有晚二疊世至早三疊世的蛇綠混雜巖、早白堊世OIB型玄武巖和去申拉組火山巖以及白堊紀的花崗質(zhì)侵入巖。
研究區(qū)出露地層主要有上白堊統(tǒng)竟柱山組(K2j),下白堊統(tǒng)朗山組(K1l)和多尼組(K1d)。竟柱山組巖性主要為礫巖、砂巖,朗山組為碳酸鹽巖,多尼組巖性主要為砂巖、礫巖。侵入巖發(fā)育有晚白堊世二長花崗巖、早白堊世中晚期二云二長花崗巖及早白堊世中期花崗閃長巖。雪如巖體呈橢圓狀近東西向展布,南北兩側(cè)與下白堊統(tǒng)朗山組呈侵入接觸,北側(cè)與上白堊統(tǒng)竟柱山組呈角度不整合接觸(圖 1b)。雪如巖體主要由二長花崗巖組成。二長花崗巖呈中細?!至=Y(jié)構(gòu)和似斑狀結(jié)構(gòu),塊狀構(gòu)造。斑晶主要由堿性長石(約10%)和石英(約5%)等組成,粒徑一般5~15 mm?;|(zhì)中主要礦物為石英、斜長石、堿性長石、黑云母等。
本文樣品采集位置見圖1b。測試工作包括鋯石U-Pb定年、鋯石微量元素分析、全巖地球化學分析。鋯石測年及鋯石微量元素分析樣品為似斑狀二長花崗巖(17BGA04-3)。全巖地球化學分析樣品為雪如巖體中中細粒二長花崗巖(17BGA18-2)、粗粒二長花崗巖(17BGA18-7)及似斑狀二長花崗巖(17BGA04-2、17BGA19-5、17BGA20-2、17BGA28-2)。樣品粉碎、鋯石挑選、制靶照相工作由北京鋯年領(lǐng)航科技有限公司完成,鋯石U-Pb測年和鋯石原位微量元素測試均在中國地質(zhì)大學(北京)進行。全巖地球化學分析在澳實分析檢測(廣州)有限公司完成。
二長花崗巖中的鋯石陰極發(fā)光CL圖像見圖2a,鋯石顆粒晶型比較完整,多為灰黑色,大多呈長柱狀,長軸80~150 μm,短軸50~90 μm,大部分發(fā)育明顯的震蕩環(huán)帶結(jié)構(gòu),少數(shù)發(fā)育核幔結(jié)構(gòu)。鋯石U-Pb定年結(jié)果顯示(表1),鋯石中Th質(zhì)量分數(shù)為270×10-6~3 872×10-6,U質(zhì)量分數(shù)351×10-6~2 315×10-6,Th/U平均比值為0.48(小于1)。綜上表明鋯石為巖漿成因[4]。18個測點中,因測點04、09諧和度較差故排除該兩點,有效測點16個,最終206Pb/238U加權(quán)平均年齡為(76±1)Ma(MSWD=0.81)(圖 2b),表明其結(jié)晶年齡為晚白堊世晚期。
圖2 雪如巖體似斑狀二長花崗巖17BGA04-3的鋯石陰極發(fā)光圖像(a)及U-Pb年齡諧和圖和加權(quán)平均年齡(b)Fig.2 Zircon CL photos(a),zircon U-Pb dating concordia diagram and weighted average ages diagram(b) from 17BGA04-3 of monzonitic granite in Xueru rocks
鋯石稀土元素總量(∑REE)為436×10-6~2 307×10-6,均值為772×10-6。稀土配分模式圖(圖 3a)表現(xiàn)為虧損LREE、富集HREE的明顯左傾型。LREE/HREE=0.05~0.59(均值為0.13),LaN/YbN=0.000 1~0.156(均值為0.015 4),輕重稀土元素分餾明顯。鋯石微量元素分析結(jié)果顯示見表2。Th、U、Hf元素富集,而Ti、La、Ce虧損,鋯石微量元素蛛網(wǎng)圖見圖3b。
巖石樣品全巖地球化學分析結(jié)果 (表 3),SiO2質(zhì)量分數(shù)為72.25%~76.15%,平均為74.05%;全堿(Na2O+K2O)質(zhì)量分數(shù)介于7.70%~8.51%之間,平均為8.11%。TAS圖解中,所有樣品均落在花崗巖范圍內(nèi)(圖 4a);里特曼指數(shù)σ=2.01~2.26(平均2.09),為鈣堿性花崗巖,根據(jù)SiO2-K2O圖解(圖 4b)進一步確認其高鉀鈣堿性的特征。鋁飽和指數(shù)A/CNK=0.96~1.04,A/CNK-A/NK圖解(圖 4c)顯示其為準鋁質(zhì)-弱過鋁質(zhì)花崗巖。二長花崗巖稀土元素總量(ΣREE)為198.39×10-6~256.34×10-6,均值為200.39×10-6,稀土配分模式圖表現(xiàn)出明顯的右傾型(圖 5a),且LREE/HREE=12.28~20.63,LaN/YbN=10.04~30.70,為輕稀土富集型花崗巖,Eu表現(xiàn)出中等負異常(δEu=0.37~0.58)。微量元素原始地幔標準化蛛網(wǎng)圖(圖 5b)顯示,大離子親石元素(Rb、Ba、Th、U、K)相對富集,高場強元素(Nb、Ta、P、Ti)相對虧損。
圖3 二長花崗巖鋯石稀土元素球粒隕石標準化配分圖解(a)和微量元素原始地幔標準化蛛網(wǎng)圖(b)Fig.3 Chondrite-normalized REE distribution patterns(a) and primitive-mantle-normalized trace elements spider diagrams(b) of zircon from monzonitic granite
表2 西藏班戈縣雪如巖體二長花崗巖(17BGA04-3)LA-ICP-MS鋯石微量元素分析結(jié)果
Table 2 Zircon trace elements (10-6) from the Xueru monzonitic granite (17BGA04-3) in Bange County,Tibet
編號12345.06789101112131415Ti5.425.752.515.3211.90.322.0313.03.003.911.519.151.671.322.15Sr0.270.590.101.150.10.090.100.20.240.200.266.267.740.600.03Y975631592890325 2697542113 5791521921126 5461893573Nb5.665.084.046.8210.34.264.042.45.822.595.055.413.117.954.07La0.274.100.044.870.10.480.630.01.590.680.661.2241.904.190.00Ce55.453.640.768.4153.639.631.226.241.736.368.762.5115.351.043.4Pr0.180.800.051.441.00.190.210.10.510.150.270.4311.960.980.03Nd2.393.521.077.9314.51.241.793.12.832.222.394.7451.035.251.28Sm3.362.251.703.8625.42.631.606.02.922.583.995.529.982.451.54Eu0.900.550.421.156.50.490.451.60.530.791.401.751.130.610.46Gd19.310.029.1815.88103.212.478.0824.513.169.7420.629.213.7214.198.82Tb6.633.383.315.4229.34.193.108.54.783.316.819.123.164.583.23Dy78.541.942.167.0322.354.536.6100.757.340.478.810635.761.341.0Ho31.018.1717.8226.61111.522.1915.6037.123.9816.530.541.512.9326.3916.73Er1459590137471.210781166.5120801391897013585Tm32.223.7821.8932.6993.024.2020.6035.228.3017.830.440.717.0833.2120.86Yb317265.8236.1352.5823.1249.0226.2329.8303.9183301.9390192.8356.0230.2Lu65.465.2055.2679.75152.954.9155.9066.170.7842.371.383.446.1684.0654.89Hf9 76110 19010 42510 2317 25710 26110 3197 944.110 6849 0739 6377 73610 18310 83310 016Ta2.551.841.842.482.72.001.520.82.021.102.011.871.342.851.39Pb9.3211.6712.0112.2112.59.1410.513.214.143.5713.3015.308.6315.8511.31Th453446504478106360341178.65191335871044294599465U560723769772594.7559667184.19102138126165451052712ΣREE757588519805230 8573483805672436756965623779508LREE62654488201453637504377762316447HREE695523475717210 6528447768622394679889392715461LREE/HREE0.090.120.090.120.100.090.080.050.080.110.110.090.590.090.10(La/Yb)N0.000.010.000.010.000.000.00?0.000.000.000.000.160.01?
注:*表示無法取值(La含量為0),微量元素為×10-6
表3 西藏班戈縣雪如巖體二長花崗巖主量元素(%)、稀土元素(×10-6)和微量元素(×10-6)分析結(jié)果
Table 3 Analysis results of major elements(%),REE(×10-6) and trace elements(×10-6) of the Xueru monzonitic granite in Bange County,Tibet
樣品17BGA042 17BGA19517BGA20217BGA28217BGA18217BGA187樣品17BGA04217BGA19517BGA20217BGA28217BGA18217BGA187SiO273.6172.2573.4774.1676.1574.68Lu0.300.360.260.210.300.36TiO20.330.390.310.260.170.28Y17.1021.4016.4012.7017.1021.40Al2O313.1813.2813.3013.2612.5113.09ΣREE256.34243.43231.30198.39131.34210.62FeO1.291.361.051.140.780.74LREE244.49229.42220.02189.45121.45197.49Fe2O30.871.010.990.880.450.44HREE11.8514.0111.288.949.8913.13MnO0.030.060.050.050.040.04LREE/HREE20.6316.3819.5121.1912.2815.04MgO0.630.810.600.440.270.43LaN/YbN25.7719.0426.1030.7010.0416.52CaO1.441.711.421.240.971.43YbN10.5312.769.126.9412.3512.82Na2O3.233.193.143.113.553.20δEu0.580.590.600.620.380.49K2O4.634.514.715.174.965.26Y/Yb9.609.9010.6010.807.809.90P2O50.120.140.110.100.040.08Li104.5075.3087.4081.7027.3036.70LOl0.840.530.650.470.460.41P523.75611.04480.11436.46174.58349.17Total100.2099.2099.80100.30100.40100.10Sc5.807.305.505.405.906.40N/K0.700.710.670.600.720.61Ti1 978.352 338.051 858.451 558.701 019.151 678.60A/CNK1.021.001.041.020.960.97V28.0033.0024.0018.009.0018.00Mg#35.1938.9235.5228.8628.8540.35Cr13.0013.0013.0010.008.0013.00R12 6072 5512 6152 5652 6142 569Co2.003.502.902.101.301.70R2444484443415363431Ni3.003.202.902.401.601.80σ2.022.032.022.202.182.26Cu1.501.700.8019.201.305.20DI87.2185.0587.3389.0893.7291.04Zn23.0031.0030.0020.0022.0017.00La64.3057.6056.4050.5064.3057.60Ga18.5517.7017.7017.2016.6017.80Ce121.50112.00109.5092.10121.50112.00Rb368.00309.00343.00360.00402.00347.00Pr12.4512.0011.059.5812.4512.00Sr197.00221.00226.00180.0070.40166.00Nd39.4040.6036.5031.6039.4040.60Zr182.00219.00150.00147.00123.00180.00Sm5.936.215.654.875.936.21Y17.1021.4016.4012.7016.3021.50Eu0.911.010.920.800.911.01Nb17.7020.0015.6018.1021.7022.20Gd3.874.443.903.173.874.44Cs15.9016.9512.8511.8516.106.18Tb0.550.620.520.420.550.62Ba252.00328.00259.00262.0067.00243.00Dy2.883.442.742.132.883.44Hf5.406.004.304.504.405.30Ho0.560.690.520.400.560.69Ta1.802.502.001.602.602.40Er1.651.971.561.261.651.97Pb35.0027.4032.6037.2034.2027.60Tm0.250.320.230.170.250.32Th42.7038.3038.8035.6041.9037.20Yb1.792.171.551.181.792.17U8.095.037.219.423.785.28
注:N/K=Na2O/K2O(wt%);A/CNK=Al2O3/(CaO+Na2O+K2O)(mol);Mg#=100×Mg/(Mg+Fe)(mol);
σ=(Na2O+K2O)2/(SiO2-43)(wt%);DI=Qz+Or+Ab+Ne+Lc+Kp。
圖4 二長花崗巖TAS圖解(a)、SiO2-K2O解(b)、A/CNK-A/NK圖解(c)、YbN-(La/Yb)N圖解(d)(據(jù)[5])Fig.4 TAS diagram(a),SiO2- K2O diagram(b), A/CNK-A/NK diagram(c), YbN-(La/Yb)N diagram(d) of monzonitic granite
圖5 二長花崗巖稀土元素球粒隕石標準化配分圖解(a)和微量元素原始地幔標準化蛛網(wǎng)圖(b)Fig.5 Chondrite-normalized REE distribution patterns(a) and primitive-mantle-normalized trace elements spider diagrams(b) of monzonitic granite
雪如巖體根據(jù)巖性可分為三種花崗巖:中細粒二長花崗巖、中粗粒斑狀二長花崗巖和似斑狀二長花崗巖。前兩者巖石分別成巖于(79.25±1)Ma和(79.72±0.51)Ma[6],而本文通過鋯石U-Pb定年,結(jié)果為(76±1)Ma,其結(jié)晶年齡較之稍晚。
本文二長花崗巖A/CNK=0.96~1.04(小于1.1),說明二長花崗巖為I型花崗巖[7]。而ACF圖解(圖 6a)和Th-Rb圖解(圖 6b),進一步確認其為I型花崗巖類型。綜上,二長花崗巖為I型花崗巖。
巖石微量元素比值顯示,二長花崗巖Zr/Hf比值為27.95~36.50(平均33.28),接近于中地殼比值(33.9)[8],Rb/Cs比值為18.23~56.15(平均29.93),接近于中地殼比值(30)[8],Th/U比值為3.78~11.05(平均6.70),高于下地殼的Th/U比值(6.00)[8]。以上微量元素比值表明雪如巖體來源于中地殼。La /Ta 值大于25代表存在地幔物質(zhì)的混染[9],本文La /Ta為23~35. 7(均值27.9),說明源區(qū)受到少量的幔源物質(zhì)混染。
實驗熔融曲線(La/Yb)N-YbN圖解(圖 4d)所示,源區(qū)殘留相為斜長角閃巖相。微量元素蛛網(wǎng)圖(圖5b)顯示Ti元素虧損,指示巖漿源區(qū)可能富含流體,流體可能是早期俯沖到深部的洋殼板片釋放的流體[10]。在Rb/Sr-Sr圖解(圖6c)和Ba-Sr圖解(圖6d)中,二長花崗巖分別沿著角閃石、斜長石的分離結(jié)晶趨勢分布,說明巖漿結(jié)晶過程中發(fā)生過角閃石、斜長石的分離結(jié)晶作用。綜上,中地殼物質(zhì)在早期俯沖洋殼釋放的流體作用下,發(fā)生部分熔融,源區(qū)殘留相為斜長角閃巖相,巖漿結(jié)晶過程中發(fā)生斜長石、角閃石的分離結(jié)晶作用,最終形成雪如二長花崗巖。
圖6 二長花崗巖ACF圖解(a)和Th-Rb圖解(b)(據(jù)[11])、Ba-Sr圖解(c)和Rb/Sr-Sr圖解(d)(據(jù)[12])Fig.6 ACF diagram(a) and Th-Rb diagram(b), Ba-Sr diagram(c) and Rb/Sr-Sr diagram(d)
巖漿巖成巖壓強利用CIPW標準礦物Ab-Or-Q圖解投圖判定。結(jié)果顯示(圖7),樣品落入1 GPa左右,表明雪如巖體形成壓強為1 GPa。壓力梯度取地殼平均值0.03 GPa/km計算[13],巖漿形成深度為33.3 km左右,介于雪如地區(qū)上地殼厚度(22.5 km)和中上地殼總厚度(49.3 km)之間[14],表明雪如巖體形成于中地殼。
圖7 標準礦物Ab-Or-Q圖解(據(jù)[15])Fig.7 Ab-Or-Q diagram for standard mineral
巖漿巖溫度的計算有兩種:鋯石飽和溫度計、鋯石Ti溫度計。鋯石飽和溫度計公式為:TZr(℃)= 12 900/(2.95+0.85M+lnDZr)-273.15[16],其中,Si、Al、Fe、Mg、Ca、Na、K、P進行原子數(shù)歸一化計算,M=(2Ca+K+Na)/(Si×Al)(均為原子百分比)。DZr為純鋯石中的Zr含量(496 000×10-6)與全巖中Zr含量的比值。計算結(jié)果顯示(表4),二長花崗巖溫度為760~810 ℃(均值786 ℃)。鋯石 Ti 溫度計公式:TTi(℃) =(5 080 ± 30)/(6.01±0.03-lg(Ti)) - 273.15[17],Ti 為鋯石中 Ti 的質(zhì)量分數(shù)(×10-6)。因為15個 測點中有9個測點的Ti含量低于檢測限,故只有6個測點的數(shù)據(jù)有效。Ti溫度計計算結(jié)果顯示(表4):二長花崗巖溫度為688~764 ℃(均值721 ℃)。鋯石Ti溫度計結(jié)果主要受到Ti含量和壓力的影響[18]。上文所求得二長花崗巖形成壓力為1 GPa,屬于低壓情況,壓力效應(yīng)為50 ℃ /GPa[18]。Ti 溫度計結(jié)果校正后應(yīng)為738~814 ℃,均值為771 ℃,與鋯石飽和溫度計接近,屬于中高溫花崗巖。
表4 二長花崗巖鋯石飽和溫度計和鋯石Ti 溫度計計算結(jié)果
Table 4 Calculated zircon saturation temperature and zircon Ti temperature of monzonitic granite
樣品TZr/℃鋯石測點TTi/℃17BGA042 797 17BGA0430169017BGA19581017BGA0430269417BGA202781 17BGA0430468817BGA282779 17BGA0430575717BGA182760 17BGA0430876417BGA187791 17BGA04312733
研究表明,與俯沖帶相關(guān)的巖漿巖的鋯石氧逸度比值非常高[19],而鋯石Ce4+/Ce3+及δEu 能夠反映鋯石氧逸度。鋯石Ce4+/Ce3+計算公式:Ce4+/Ce3+=(Ce熔體-Ce鋯石/DCe3+)/(Ce鋯石/DCe4+-Ce熔體),其中,Ce鋯石、Ce熔體分別為鋯石中和全巖中的 Ce 含量;DCe3+、DCe4+分別為Ce3+、Ce4+在鋯石礦物相和熔體相的分配系數(shù),計算方式參考文獻[20]。Ce4+/Ce3+及δEu的結(jié)算結(jié)果如表5所示,二長花崗巖鋯石Ce4+/Ce3+=1.38~939.5,均值為170.15;鋯石δEu=0.26~0.48,均值為0.37,其較低的氧逸度暗示該巖體的形成可能與洋殼的俯沖無關(guān)。
前人對班怒洋的演化過程做過大量的研究??p合帶中MOR型蛇綠巖年齡介于254~217 Ma,表明主洋盆可能在晚二疊世至早三疊世時期發(fā)生裂解[21];班公湖SSZ型蛇綠巖年齡介于177~162 Ma[22],說明班怒洋盆至少在中侏羅世開始俯沖消減;對于班怒帶俯沖極性的認識主要存在三種說法:北向俯沖[23-25],南向俯沖[26-28],雙向俯沖[27-28]。與蛇綠混雜巖呈角度不整合接觸的上白堊統(tǒng)竟柱山組地層年齡測定為96 Ma[29],說明班怒洋閉合于晚白堊世之前。
本文通過全巖地球化學數(shù)據(jù)投圖所示,R1-R2圖解(圖 8a)中顯示雪如巖體樣品主要位于同碰撞花崗巖區(qū)域,Rb/10-Hf-Ta×3圖解(圖 8b)顯示為碰撞大地構(gòu)造背景下的花崗巖。Nb-Y圖解(圖 8c)顯示,雪如巖體二長花崗巖均落在火山弧或同碰撞區(qū)域。綜上,雪如巖體可能形成于同碰撞階段。但Rb-(Y+Nb)圖解(圖 8d)顯示,巖石樣品均落在同碰撞與后碰撞花崗巖疊加區(qū)域,因此在討論構(gòu)造環(huán)境時還需結(jié)合區(qū)域構(gòu)造演化過程。
通過統(tǒng)計班怒帶內(nèi)成巖于80 Ma左右的花崗質(zhì)巖體的形成年代和構(gòu)造環(huán)境(表 6),筆者發(fā)現(xiàn)這些花崗質(zhì)巖體多形成于后碰撞伸展階段,成巖年齡自97.4 Ma至76.1 Ma,故判斷雪如巖體可能也形成于后碰撞伸展環(huán)境。結(jié)合前人和本文研究成果,可總結(jié)出班怒帶的構(gòu)造演化過程如下:班怒洋于晚二疊世至早三疊世(254~217 Ma)開始裂解,中侏羅世由擴張轉(zhuǎn)換為俯沖(北向俯沖),早白堊世晚期至晚白堊世早期(116~96 Ma)班怒洋洋盆閉合,晚白堊世(100~76 Ma)兩側(cè)地體發(fā)生碰撞并進入后碰撞伸展階段。
表5 鋯石Ce4+/Ce3+和δEu計算結(jié)果
圖8 二長花崗巖R1-R2圖解(a)(據(jù)[30])、Rb/10-Hf-3×Ta圖解(b)(據(jù)[31])、Y-Nb(c)(據(jù)[32]) 和Rb-(Y+Nb)(d)(據(jù)[33])Fig.8 Digrams of R1-R2 (a), Rb/10-Hf-3×Ta (b), Y-Nb (c), Rb-(Y+Nb) (d)
表6 班公湖—怒江成礦帶80 Ma左右花崗質(zhì)侵入巖體特征對比
班戈縣雪如巖體二長花崗巖成巖年齡為(76±1)Ma,形成于晚白堊世晚期,形成壓強為1 GPa,形成于深度約33.3 km的中地殼。鋯石飽和溫度計和Ti溫度計計算結(jié)果顯示其為中高溫花崗巖,較低的氧逸度說明其形成于還原環(huán)境,表明本期巖漿活動與俯沖事件無關(guān)。二長花崗巖為I型花崗巖,物質(zhì)來源為中地殼,受到少量地幔物質(zhì)的混染。在早期俯沖洋殼釋放的流體的作用下,中地殼物質(zhì)發(fā)生部分熔融,源區(qū)主要殘留斜長角閃巖相,在巖漿結(jié)晶過程中發(fā)生角閃石和斜長石的分離結(jié)晶作用。二長花崗巖形成于后碰撞伸展環(huán)境,表明晚白堊世晚期(76 Ma)班怒帶中段地區(qū)已經(jīng)進入后碰撞伸展階段。
致謝感謝中國地質(zhì)大學(北京)相鵬老師在鋯石U-Pb測年實驗上的大力協(xié)助,感謝李發(fā)橋同學在撰寫論文過程中提供的幫助。