孫賽軍 張麗鵬 丁興 孫衛(wèi)東 張兆榮

1.中國科學院廣州地球化學研究所，礦物學與成礦學重點實驗室，廣州 510640

2.中國科學院大學，北京 100049

3.中國科學院廣州地球化學研究所，同位素地球化學國家重點實驗室，廣州 510640

4.中國科學院青藏高原地球科學卓越創(chuàng)新中心，北京 100101

5.福建省閩西地質(zhì)大隊，廈門 361004

1 引言

拉薩地塊一般指位于北部的班公湖-怒江縫合帶(BNSZ)和南側(cè)的印度河-雅魯藏布江縫合帶(IYZS)之間的近東西向的狹長地域，長約2500km，南北寬150～300km，面積達4.5 ×105km2的巨型構(gòu)造巖漿帶。班公湖-怒江縫合帶是拉薩地塊與羌塘-保山地塊的匯聚縫合線，兩側(cè)發(fā)育了大量的變質(zhì)巖和巖漿巖(Ding et al.，2003;Kapp et al.，2003;Zhu et al.，2011;Sui et al.，2013;常青松等，2011;黃玉等，2012;吳浩等，2013)。這些巖石記錄了這兩個板塊在青藏高原形成過程中由新特提斯洋俯沖到陸-陸碰撞、拼合的歷史。通過這些巖石的成因研究，將有助于更深入了解青藏高原造山機制和深部動力學過程。

白堊世時期，是現(xiàn)今歐亞大陸匯聚拼合的重要時期，也是青藏高原各微陸塊匯聚拼合的關(guān)鍵時期(Zhu et al.，2013;Ding et al.，2014)。通過對拉薩地塊班公湖-怒江縫合帶南部(岡底斯中北部)的早白堊世巖漿巖的研究，前人多認為它們的形成與新特提斯洋殼北向俯沖有關(guān)(Xu et al.，1985;Coulon et al.，1986;Pearce and Mei，1988;Ding et al.，2003);但是近年更深入廣泛的研究表明，它們可能與班公湖-怒江特提斯洋殼南向俯沖所引起的板片斷離有關(guān)(Zhu et al.，2009b，2011，2013;莫宣學等，2005;潘桂棠等，2006;朱弟成等，2008;陳越等，2010;張亮亮等，2010，2011;張曉倩等，2010;黃玉等，2012;吳浩等，2013)?？梢姡P(guān)于這些早白堊世巖漿巖成因的更清晰的動力學過程仍有待更多研究的解析。本文通過研究拉薩地塊中北部那曲地區(qū)的早白堊世火山巖的鋯石U-Pb 年齡、Hf 同位素和其地球化學特征，來討論那曲地區(qū)火山巖的巖漿源區(qū)、巖石成因及構(gòu)造背景意義，以期為拉薩地塊班公湖-怒江縫合帶南部巖漿巖的成因模型提供制約，繼而為更深入了解青藏高原造山機制和深部動力學過程提供新的資料。

表1 那曲研究區(qū)巖石地層簡表Table 1 The stratigraphic profile of the Nagqu study area

圖1 拉薩地塊在青藏高原中的地理位置(a，據(jù)Zhu et al.，2011)和那曲地質(zhì)簡圖(b、c)Qhal、Qhpl、Qhf、Qpl-全新世沉積物;E2γδ-始新世細粒角閃黑云花崗閃長巖;K1d1、K1d2-多尼組上段、下段;K2ξγab-中粒黑云母花崗巖;K2δο-中粒角閃黑云石英閃長巖;K2ηγ-中細粒黑云二長花崗巖;J2-3 l2-拉貢塘組上段;J2-3 l1-拉貢塘組下段;J2s-桑卡拉傭組;J2m2-馬里組上段;Oph(T)-橄欖巖、輝長巖;JSSZ-金沙江縫合帶;BNSZ-班公湖-怒江縫合帶;YZSZ-雅魯藏布縫合帶.黑色圓點標注為采樣位置Fig.1 Location of Lhasa Terrane and simplified geological map of the studied area (after Zhu et al.，2011)

2 區(qū)域地質(zhì)背景與樣品描述

那曲地區(qū)位于青藏高原腹地，西藏自治區(qū)北部(圖1)，地處兩個板片結(jié)合地帶，南為岡底斯-念青唐古拉板片北緣(桑雄-麥地卡陸緣巖漿弧帶，為岡底斯-念青唐古拉板片的一部分)，北為班公湖-怒江結(jié)合帶中段，班-怒結(jié)合帶明顯收斂變窄的部位。測區(qū)地層屬岡底斯-騰沖地層區(qū)班戈-八宿地層分區(qū)之桑雄-麥地卡地層小區(qū)，出露地層主要為中生代地層，詳見表1。本文采集的中酸性火山巖樣品來自于多尼組上段，位于那曲南東方向約95km 處，屬于哈爾麥火山巖范疇。哈爾麥火山巖發(fā)育于哈爾麥鄉(xiāng)附近，被認為夾于拉貢塘組砂板巖中，此套火山巖巖石類型以安山巖、安山質(zhì)火山角礫巖和流紋質(zhì)火山角礫凝灰?guī)r為主(李奮其等，2010)。從野外觀察來看，本次采集火山巖樣品風化較嚴重、露頭出露一般，產(chǎn)狀不清，與拉貢塘組之間應(yīng)為不整合接觸關(guān)系(彭智敏等，2011)。

安山巖 巖石呈淺灰-淺灰綠色，殘余斑狀結(jié)構(gòu)，塊狀構(gòu)造。斑晶主要為斜長石，更長石占25%，粒徑0.5～2.5mm;基質(zhì)主要為斜長石，其中微粒更長石占50%。也見絹云母呈鱗片狀不均勻分布其中，可見少量石英。副礦物為磷灰石、褐鐵礦等。

流紋質(zhì)火山角礫凝灰?guī)r 巖石呈灰綠色，巖屑晶屑火山角礫凝灰結(jié)構(gòu)，塊狀構(gòu)造?；鹕浇堑[占20%，主要由酸性火山巖，晶屑凝灰?guī)r等組成，粒徑0.5～1.5cm，多呈熔蝕狀;晶屑成份有鈉更長石、石英、正長石及少量金屬礦物等，呈次棱角狀，熔蝕狀;巖屑成分有酸性霏細狀長英質(zhì)，呈棱角狀或長橢圓狀。巖石中的膠結(jié)物多為次生的綠泥石。

3 分析方法

挑選較新鮮無污染樣品碎至200 目，主微量元素分析在福建省地質(zhì)礦產(chǎn)局三明實驗室完成。主量元素采用原子吸收分光光度法在GGX-9 儀器上分析，微量元素在熱電ICP 6300 型電感耦合等離子體發(fā)射光譜儀上測定，精度均優(yōu)于5%。

鋯石分選是在河北省區(qū)域地質(zhì)調(diào)查研究院完成。鋯石制靶在中國科學院廣州地球化學研究所完成，首先在雙目鏡下挑選出晶形完整、無包裹體和裂隙的鋯石，和標樣TEMORA 粘于雙面膠上，然后用環(huán)氧樹脂固化制靶，鋯石靶打磨、拋光至將近一半進行鋯石陰極發(fā)光(CL)顯微照相。鋯石U-Pb 定年在中科院廣州地球化學研究所利用LA-ICPMS 測定。儀器采用美國Resonetics 公司的RESOlution M-5激光剝蝕系統(tǒng)和Agilent 7500a 型的ICP-MS。儀器工作參數(shù):激光剝蝕斑束為33μm，頻率為8Hz，剝蝕時間為40s，背景積分時間25s(涂湘林等，2011;Li et al.，2012)。鋯石年齡計算采用國際標準鋯石TEMORA 作為外標，元素含量采用美國國家標準物質(zhì)局人工合成硅酸鹽玻璃NIST SRM610 作為外標，29Si 作為內(nèi)標元素進行校正。樣品的同位素比值和元素含量數(shù)據(jù)處理采用軟件ICPMSDataCal8.3 (Liu et al.，2008)。鋯石諧和年齡繪圖采用軟件Isoplot (Ludwig，2003)。

鋯石Hf 同位素分析是基于U-Pb 定年同點或同等效位置在中科院廣州地球化學研究所同位素地球化學國家重點實驗室LA-MC-ICPMS 上完成。實驗運行條件斑束為45μm，頻率8Hz，推薦176Hf/177Hf =0.282906 ±0.0000010 (2σ)的蓬萊鋯石作為標樣(Li et al.，2010)。具體分析流程見Wu et al.(2006)。

圖2 那曲中酸性火山巖的鋯石CL 圖像、U-Pb 定年點(紅色實線圓圈)與Hf 同位素測試點(黑色實線圓圈)Fig.2 Cathodoluminescence (CL)images of zircon from Nagqu intermediate-felsic volcanic rocks Red circles indicate locations of U-Pb dating，while black circles indicate Hf analyses positions

4 分析結(jié)果

4.1 鋯石U-Pb 年齡

西藏岡底斯那曲地區(qū)火山巖2 件樣品的鋯石分析結(jié)果見表2。

流紋質(zhì)英安巖PM201-2 樣品鋯石多為長、短柱狀，半自形-自形晶形，多數(shù)鋯石具有明顯的巖漿振蕩環(huán)帶，部分鋯石顯示核邊結(jié)構(gòu)(圖2)。鋯石Th/U 比值為0.38～0.92(＞0.1)，均值為0.49，顯示典型的巖漿鋯石特征(Hoskin and Black，2000)。除去繼承核年齡點和由于包體或打穿引起的諧和度不好的數(shù)據(jù)點外，其余年齡點均落在一致曲線上或附近(圖3)，顯示較弱的Pb 丟失，它們的加權(quán)平均年齡為109.2 ±3.5Ma(MSWD =3.5)，代表英安巖漿活動的時代。而位于核部的年齡點PM201-2-16，其Th/U 比值與巖漿鋯石相似，207Pb/206Pb 年齡為1439Ma，該繼承年齡的出現(xiàn)說明它有可能來源于源區(qū)古老殘留物或者是英安巖漿上侵過程中捕獲了元古代的物質(zhì)。

安山巖樣品PM205-1 中的鋯石較之PM201-2 普遍較大，短柱狀，自形-半自形晶為主，具明顯的巖漿振蕩環(huán)帶(圖2)。鋯石Th/U 比值介于0.56～1.78 之間，均值為0.95，屬巖漿成因鋯石。排除由于包體或鉛丟失引起的諧和度不好的數(shù)據(jù)點外，其206Pb/238U 加權(quán)平均年齡為111.9 ±1.2Ma(MSWD=1.4)，代表了安山巖的成巖時代。

4.2 全巖地球化學

12 件那曲多尼組火山巖樣品主微量元素分析結(jié)果見表3，PM201-1～PM201-5 及PM205-1 樣品均發(fā)生了不同程度的熱液蝕變，燒失量較大(2.60%～8.86%)，故如Ba、K、Na、Rb、Sr、U 等活動性元素不宜用來進行相關(guān)解釋和判別。在蝕變過程中，一般認為高場強元素(Nb、Ta、Zr、Hf 等)、相容元素(Cr、Ni)和稀土元素受到的影響較小，可用來討論蝕變巖石的類型和成因(Winchester and Floyd，1977;Hastie et al.，2007)。

表2 那曲中酸性火山巖的鋯石LA-ICP-MS U-Pb 定年結(jié)果Table 2 Zircon LA-ICP-MS U-Pb data for Nagqu intermediate-felsic volcanic rocks

圖3 那曲中酸性火山巖的鋯石U-Pb 年齡諧和圖Fig.3 The zircon U-Pb concordia diagrams for Nagqu intermediate-felsic volcanic rocks

表3 那曲中酸性火山巖的全巖主量(wt%)與微量元素(×10 -6)組成Table 3 Whole rock major element (wt%)and trace elements (×10 -6)data for Nagqu intermediate-felsic volcanic rocks

續(xù)表3Continued Table 3

圖4 那曲中酸性火山巖巖石類型劃分圖解(據(jù)Winchester and Floyd，1977)Fig.4 Diagrams of classification of Nagqu intermediatefelsic volcanic rocks (after Winchester and Floyd，1977)

12 件那曲火山巖樣品的SiO2含量為51.88%～74.62%，平均為66.03%;K2O + Na2O 含量為4.89%～10.19%，平均為7.65%。在SiO2-Zr/TiO2圖解中(圖4)，PM201-1～PM201-5 和PM205-1～PM205-7 樣品大致分為中性安山巖和酸性流紋英安巖，大多落于鈣堿性系列范圍內(nèi)。同時A/CNK 比值為0.90～1.48，平均為1.17，并且標準礦物中幾乎都出現(xiàn)剛玉分子，少見透輝石(Di)出現(xiàn)，屬過鋁質(zhì)巖石類型。在主量元素的Harker 圖解(圖5)中，隨著SiO2含量增加，Al2O3、CaO、MgO、TiO2的含量與其呈現(xiàn)一定的負相關(guān)性，暗示PM201 和PM205 之間存在密切的成因聯(lián)系，樣品也經(jīng)歷過分離結(jié)晶作用。

那曲火山巖樣品的稀土總量為148.5 ×10-6～421.1 ×10-6，具有明顯的Eu 負異常(δEu=0.15～0.68)的輕稀土富集型((La/Yb)N=8.79～14.85)配分模式(圖6a)?？傮w上，那曲火山巖稀土配分模式與上地殼及前人研究(黃玉等，2012)類似，但稀土總含量明顯偏高(圖6a)。原始地幔標準化微量元素蛛網(wǎng)圖(圖6b)中，樣品相對富集不相容元素Th、Hf、Zr、Rb 等元素，而虧損P、Ba、Sr、Nb、Ta、Ti 等元素，與上地殼有著相似的變化趨勢。P、Nb、Ta、Ti 的虧損常認為與俯沖作用有關(guān)的火成巖的特征(Condie，2001)。

表4 那曲中酸性火山巖鋯石Ti 溫度和Ce(Ⅳ)/Ce(Ⅲ)比值Table 4 Ti-in-zircon thermometers and Ce (Ⅳ)/Ce (Ⅲ)rations of Nagqu intermediate-felsic volcanic rocks

圖5 那曲中酸性火山巖地球化學圖解(a，據(jù)Zorpi et al.，1991)Fig.5 The geochemical diagrams of Nagqu intermediate-felsic volcanic rocks (a，after Zorpi et al.，1991)

4.3 成巖溫度和氧逸度

巖漿鋯石常有Ce 的正異常，并且Ce 的異常程度取決于巖漿的氧逸度。Ce 有+3 價和+4 價，并且鋯石中的Zr4+的半徑(離子半徑0.97?)與Ce4+離的半徑(離子半徑0.84?)較相近，因此可以通過Ce4+/Ce3+比值來反映巖漿氧逸度的變化(Zhang et al.，2013)。本文根據(jù)Ballard et al.(2002)提出的Ce 氧逸度公式得出那曲火山巖樣品PM201-2 和PM205-1 的Ce(Ⅳ)/Ce(Ⅲ)具有類似的特征，即在相似的年齡范圍內(nèi)Ce 比值普遍小于300，并且大多均小于100(表4)，與西藏東緣玉龍非成礦巖體相似(＜120)，指示巖體具有較低的氧逸度(梁華英等，2006)。我們用鋯石Ti 溫度計計算得出(Watson et al.，2006;Ferry and Watson.，2007)那曲中酸性火山巖樣品PM201-2 約～720℃，PM205-1 約～710℃，其中PM205-1-22 由于具有異常高的Ti 含量，其對應(yīng)的溫度為1930℃(表4)?？傮w上指示那曲地區(qū)早白堊世火山巖的巖漿具有低的溫度。

圖6 那曲中酸性火山巖球粒隕石標準化稀土元素配分圖(a)和原始地幔標準化微量元素蛛網(wǎng)圖(b)(標準化值據(jù)Sun and McDonough，1989)數(shù)據(jù)來源:灰色陰影(黃玉等，2012);上地殼(Rudnick and Gao，2003)Fig.6 Chondrite-normalized REE patterns (a)and primitive-mantle-normalized trace element patterns (b)for Nagqu intermediate-felsic volcanic rocks (normalization values after Sun and McDonough，1989)Data sources:The shades of gray (after Huang et al.，2012);Upper continental crust (after Rudnick and Gao，2003)

圖7 那曲中酸性火山巖εHf(t)-Age 圖解Fig.7 εHf(t)-Age diagram of Nagqu intermediate-felsic volcanic rocks

4.4 鋯石Hf 同位素

那曲中酸性火山巖樣品(PM201-2、PM205-1)鋯石Hf 同位素分析結(jié)果見表5。樣品中鋯石的176Yb/177Hf 和176Lu/177Hf 比值范圍分別為:0.026231～0.139581 和0.000682～0.003431。176Lu/177Hf 比值大多小于0.002，表明這些鋯石在形成后基本上沒有明顯的放射性成因Hf 的積累，故所測定的176Hf/177Hf 比值可代表其形成時的Hf 同位素組成(吳福元等，2007)。樣品鋯石176Hf/177Hf 比值為0.281861～0.282730，εHf(t)值基本為負值(-13.6～0.5)，其中PM201-2 鋯石樣品中核部兩點(PM201-2-03 和PM201-2-09)更為偏負(-30.4 和-16.6)(圖7)，所對應(yīng)的Hf 同位素模式年齡也較老(2481Ma 和1785Ma)。樣品Hf 模式年齡大多集中于918～1336Ma，平均為1128Ma。與前面所得的鋯石繼承核年齡1439Ma 較為接近，我們推測鋯石繼承核可能來源于古老源區(qū)殘留物，而后期經(jīng)歷了～110Ma 變質(zhì)巖漿事件。

5 討論

5.1 巖漿源區(qū)與巖石成因

本次研究區(qū)多尼組火山巖主要以中酸性火山巖為主，玄武巖幾乎露頭不可見。一般認為中酸性火山巖的成因有二種可能:(1)來自于玄武質(zhì)巖漿中廣泛的分離結(jié)晶作用及同化混染作用演化而產(chǎn)生，故與玄武巖具有相似的幔源地球化學特征(Bacon and Druitt，1988;Bonin，2004);(2)與玄武巖不同同位素特征的幔源基性巖漿引起的地殼物質(zhì)部分熔融形成(Tepper et al.，1993;Guffanti et al.，1996)。因為玄武質(zhì)巖漿只有經(jīng)過強烈的結(jié)晶分離作用(＞90%)才能產(chǎn)生酸性熔體(Wilson，1993)，所以一般情況下第一種過程產(chǎn)生的大多是小規(guī)模的中酸性熔體(Shinjo and Kato，2000)。鑒于目前在研究區(qū)及附近幾乎不可見玄武質(zhì)火山巖及缺乏相應(yīng)的研究，作者認為第二種過程更能產(chǎn)生以中酸性火山巖為主的多尼組火山巖系列?？抵緩姷?2009)研究的拉薩北部的白堊紀多尼組火山巖也證明了這一觀點。

那曲中酸性火山巖富集Th、Zr、Hf，虧損Nb、Ta、Ti、P 等高場強元素，具有類似于島弧火山巖的特征(圖6b)。樣品的La/Nb 比值高于陸殼均值(2.5)，而Nb/U 比值普遍低于陸殼均值(6.2)(Rudnick and Gao，2003)，簡單的陸殼混染難以解釋。MgO 與Nb/La、Hf/Sm 之間無相關(guān)關(guān)系(圖8)，說明高場強元素的異常是巖漿的原始特征(李海勇等，2008)。同時樣品具有較低的MgO 和Mg#＜40，也暗示了其形成可能與下地殼鎂鐵質(zhì)巖石部分熔融有關(guān)(Rapp et al.，1999)。從Yb/Ta-Y/Nb 圖(圖9)中可以看出安山質(zhì)火山巖和幔源組分高的流紋質(zhì)火山巖都在下地殼和虧損地?；旌暇€上，只是大部分流紋質(zhì)樣品偏離，而偏離的原因是Y/Nb 比值降低。由于在中酸性巖漿中Y 更易進入角閃石、黑云母、磷灰石、鋯石中，而隨著這些礦物的分離結(jié)晶，殘余熔體Y/Nb 比值變小。可見，簡單的陸殼與虧損地?；旌喜荒芡昝澜忉屗袠悠?，樣品的分離結(jié)晶作用或許是一種可行的解釋，同時我們發(fā)現(xiàn)地殼的參與(＞80%)在那曲早白堊世火山巖形成過程中起著更重要的作用(圖9)。

表5 那曲中酸性火山巖的鋯石Hf 同位素組成Table 5 Zircon Hf isotope data for Nagqu intermediate-felsic volcanic rocks

圖8 那曲中酸性火山巖MgO 與Nb/La (a)和Hf/Sm (b)的關(guān)系圖解Fig.8 Nb/La-MgO (a)and Hf/Sm-MgO (b)diagrams of Nagqu intermediate-felsic volcanic rocks

圖9 那曲中酸性火山巖Yb/Ta-Y/Nb 圖解數(shù)據(jù)來源:BCC-平均大陸地殼(Rudnick and Gao，2003);LCC-大陸下地殼(Rudnick and Gao，2003);DMM-虧損地幔(Salters and Stracke，2004)Fig.9 Yb/Ta-Y/Nb discrimination diagram of Nagqu intermediate-felsic volcanic rocksData source:BCC-Bulk Continental Crust (Rudnick and Gao，2003);LCC-Lower Continental Crust (Rudnick and Gao，2003);DMM-Depleted MORB Mantle (Salters and Stracke，2004)

樣品SiO2與Al2O3、P2O5、TiO2等呈現(xiàn)負相關(guān)關(guān)系，暗示了在形成中酸性火山巖過程中經(jīng)歷了結(jié)晶分異作用。研究區(qū)內(nèi)流紋英安巖和流紋巖樣品具有較高的FeOT/MgO 值和分異指數(shù)，是明顯經(jīng)歷了強結(jié)晶分異作用形成。高分異流紋質(zhì)火山巖微量元素蛛網(wǎng)圖曲線與安山巖曲線變化趨勢相近，與安山巖相比，高分異流紋質(zhì)火山巖的Eu、Sr、Ba、P、Ti 更加虧損，在酸性巖漿體系中，由于Sr、Eu 及Ba 在斜長石和鉀長石中為強相容元素，那么它們的虧損表明發(fā)生了斜長石和鉀長石的分離結(jié)晶作用。而Ti、P 的虧損很可能與含Ti 礦物(如鈦鐵礦和金紅石)及磷灰石的分離有關(guān)。本文和前人數(shù)據(jù)(黃玉等，2012)的微量元素及稀土元素配分曲線與上地殼具有近一致的曲線特征，但是區(qū)內(nèi)少有玄武巖出露，難以用基性巖漿的結(jié)晶分異作用來解釋大面積的中酸性巖漿活動，筆者認為可能與地殼物質(zhì)部分熔融有關(guān)。由于研究區(qū)內(nèi)安山巖和流紋質(zhì)火山巖在地球化學方面有著高度的相似性，同時在圖10 中我們可以看出安山質(zhì)火山巖和流紋質(zhì)火山巖經(jīng)歷了角閃石和長石的不同程度的分離結(jié)晶作用(Ding et al.，2014)。因此筆者初步認為，二者是同源巖漿受不同程度結(jié)晶分異作用的產(chǎn)物。

圖10 那曲中酸性火山巖Rb/Sr-Sr 圖解(據(jù)Ding et al.，2014)Fig.10 Rb/Sr-Sr diagram of Nagqu intermediate-felsic volcanic rocks (after Ding et al.，2014)

研究區(qū)樣品具有與西藏東緣玉龍非成礦巖體相似(＜120)的Ce 比值，指示巖體具有較低的氧逸度(梁華英等，2006)，同時在該區(qū)附近也未見礦體出露。安山巖樣品鋯石εHf(t)值為- 5.5～0.5，對應(yīng)的Hf 模式年齡為918～1226Ma。流紋巖鋯石εHf(t)大多為為-2.8～-13.6，對應(yīng)的Hf 模式年齡集中在1084～1636Ma 之間，其中鋯石核部區(qū)PM201-2-03 (107.6Ma)和PM201-2-09 (105.6Ma)顯示更負的εHf(t)分別為-30.4 和-16.6，相應(yīng)的Hf 模式年齡為2481Ma 和1785Ma。樣品不均一的Hf 同位素成分，需要一個開放體系來引起熔體中176Hf/177Hf 比值的變化，不均一的εHf(t)可能來源于虧損地幔與地殼的混合(Kemp et al.，2007)。并且這些火山巖的虧損地幔組分最多不超過20%(圖9)，可能是微陸塊的古老巖石圈地幔。樣品較低的成巖溫度(～715℃)與幔源物質(zhì)貢獻量不超過20%相符，不同于藏東緣察隅較高的成巖溫度(789～812℃)(Zhu et al.，2009a)。

那么古老巖石圈地幔是中拉薩地塊下伏的古老巖石圈地幔還是北邊具寒武紀結(jié)晶基底的安多微陸塊?前人研究表明，110Ma 時具寒武紀結(jié)晶基底的安多微陸塊與北拉薩地塊北緣已拼接在一起(Zhu et al.，2011)，考慮到晚白堊至古新世期間拉薩地殼的縮短，約110Ma 時那曲地區(qū)與具有寒武紀結(jié)晶基底的中拉薩地塊距離應(yīng)該大于100km，并且在桑巴地區(qū)酸性侵入巖鋯石εHf(t)=-6.0～5.7，說明了那曲與中拉薩之間的地塊至少有新生地殼的生成。因此古老巖石圈可能是具寒武紀結(jié)晶基底的安多微陸塊(Kapp et al.，2003;Zhu et al.，2011)。

表6 多尼組中基性火山巖與其他地區(qū)中基性火山巖在微量元素(×10 -6)的對比Table 6 Comparison between trace elements (×10 -6)of basic volcanic rocks from Duoni Formation and from other areas

綜合以上證據(jù)，我們初步認為研究區(qū)那曲早白堊世火山巖很可能來源于:(1)班公-怒江洋南向俯沖過程中受流體改造的安多微陸塊古老巖石圈地幔(貢獻小部分)與(2)羌塘和北拉薩地塊陸-陸碰撞地殼加厚過程中下地殼部分熔融(貢獻大部分)及后期經(jīng)歷不同程度的分離結(jié)晶作用形成的。

5.2 構(gòu)造環(huán)境

拉薩地塊北部早白堊世火山巖從西到東依次為班-怒帶內(nèi)及南緣斷續(xù)出露的去申拉組，烏木壟鉛波巖組，多尼組及臥榮溝組火山巖。前人在詳細的地球化學及年代學基礎(chǔ)上，提出拉薩地塊中北部在113Ma 左右發(fā)生的大規(guī)模的巖漿活動，可能是班公湖-怒江洋巖石圈南向俯沖過程中板片斷離導致的(Zhu et al.，2011;Sui et al.，2013;Chen et al.，2014)。由于Yb 為不活潑元素，與不相容元素行為類似，那么在部分熔融和分離結(jié)晶過程中，Th/Yb、Nb/Yb 比值將保持不變(Pearce，1983;Rollinson，1993)，在該圖解中(圖11)樣品均落在大陸邊緣弧火山巖范圍內(nèi)，暗示那曲早白堊世中酸性火山巖成因可能與班公湖-怒江洋的南向俯沖有關(guān)。并且樣品富集Th、U 和輕稀土元素，虧損Nb、Ta、Ti 和P 等高場強元素，暗示源區(qū)曾受到與俯沖相關(guān)流體或熔體的改造。中生代拉薩地塊處于班公-怒江南向俯沖和雅魯藏布新特提斯北向雙向俯沖下，導致了拉薩地塊的縮短、陸殼加厚并產(chǎn)生了大量的弧火山巖、弧盆體系(潘桂棠等，2006)。康志強等(2009)將多尼組中基性火山巖與新西蘭Egmont 火山巖、阿根廷西北部火山巖及則弄群火山巖進行對比發(fā)現(xiàn)(表6)，它們具有相似的微量元素分布特征(Zr、Hf、Y、Yb 等)，并且本次研究的安山巖樣品具有同樣相似的特征。由于La/Yb 可以指示地殼厚度并反應(yīng)巖漿來源的深度(Kay et al.，1991)，他們發(fā)現(xiàn)多尼組火山巖有著與則弄群類似的巖漿形成深度，暗示了當時北拉薩地塊已經(jīng)加厚。班公怒江洋南向俯沖開始于中二疊紀(～263Ma)，終止于晚白堊紀(～113Ma)，之后轉(zhuǎn)入羌塘地塊與拉薩地塊陸-陸碰撞環(huán)境(Zhu et al.，2011)。而本文早白堊世中酸性火山巖加厚地殼的構(gòu)造背景很可能是班公湖-怒江洋俯沖結(jié)束進入陸-陸碰撞過程中引起的。

圖11 那曲中酸性火山巖Th/Yb-Nb/Yb 圖解Fig.11 Th/Yb-Nb/Yb discrimination diagram for Nagqu intermediate-felsic volcanic rocks

6 結(jié)論

北拉薩地塊那曲地區(qū)哈爾麥鄉(xiāng)中酸性火山巖漿活動侵位于110Ma 左右，是北部拉薩地塊～113Ma 巖漿大爆發(fā)事件的產(chǎn)物之一;

那曲早白堊世中酸性火山巖富集Th、Zr、Hf，虧損Nb，Ta，Ti，P 等高場強元素，顯示為負的εHf(t)值，暗示其形成來源于班公-怒江南向俯沖過程中受流體改造的安多微陸塊古老巖石圈地幔(貢獻小部分)與羌塘和北拉薩地塊陸陸碰撞地殼加厚過程中下地殼部分熔融(貢獻大部分)，并經(jīng)歷了不同程度的分離結(jié)晶作用形成的。

致謝 主微量元素測試得到福建省地質(zhì)礦產(chǎn)局三明實驗室工作人員的幫助;LA-ICP-MS 鋯石U-Pb 定年得到廣州地化所李聰穎博士的幫助;鋯石Hf 同位素測試過程中得到過廣州地球化學研究所張樂的幫助;另外在本文修改過程中得到胡永斌、張潺蟬、祝紅麗、吳凱、劉芳等的幫助;在此表示誠摯的感謝。

Bacon CR and Druitt TH.1988.Compositional evolution of the zoned calcalkaline magma chamber of Mount Mazama，Crater Lake，Oregon.Contributions to Mineralogy and Petrology，98(2):224-256

Ballard JR，Palin JM and Campbell IH.2002.Relative oxidation states of magmas inferred from Ce(IV)/Ce(III)in zircon:Application to porphyry copper deposits of northern Chile.Contributions to Mineralogy and Petrology，144(3):347-364

Bonin B.2004.Do coeval mafic and felsic magmas in post-collisional to within-plate regimes necessarily imply two contrasting，mantle and crustal，sources?A review.Lithos，78(1-2):1-24

Bouvier A，Vervoort JD and Patchett PJ.2008.The Lu-Hf and Sm-Nd isotopic composition of CHUR:Constraints from unequilibrated chondrites and implications for the bulk composition of terrestrial planets.Earth and Planetary Science Letters，273(1-2):48-57

Chang QS，Zhu ZD，Zhao ZD，Dong GC，Mo XX，Liu YS and Hu ZC.2011.Zircon U-Pb geochonology and Hf isotopes of the Early Cretaceous Rena-Co rhyolites from southern margin of Qiangtang，Tibet，and their implications.Acta Petrologica Sinica，27(7):2034-2044 (in Chinese with English abstract)

Chen Y，Zhu DC，Zhao ZD，Zhang LL，Liu M，Yu F，Guan Q and Mo XX.2010.Geochronology，geochemistry and petrogenesis of the Bamco andesites from the northern Gangdese， Tibet.Acta Petrologica Sinica，26(7):2193-2206 (in Chinese with English abstract)

Chen Y，Zhu DC，Zhao ZD，Meng FY，Wang Q，Santosh M，Wang LQ，Dong GC and Mo XX.2014.Slab breakoff triggered ca.113Ma magmatism around Xainza area of the Lhasa Terrane，Tibet.Gondwana Research，26(2):449-463

Condie KC.2001.Mantle Plume and Their Record in Earth History.London:Cambridge University Press

Coulon C，Maluski H，Bollinger C and Wang S.1986.Mesozoic and Cenozoic volcanic rocks from central and southern Tibet:39Ar-40Ar dating，petrological characteristics and geodynamical significance.Earth and Planetary Science Letters，79(3-4):281-302

Ding L，Kapp P，Zhong DL and Deng WM.2003.Cenozoic volcanism in Tibet: Evidence for a transition from oceanic to continental subduction.Journal of Petrology，44(10):1833-1865

Ding X，Sun WD，Chen WF，Chen PR，Sun T，Sun SJ，Lin CT and Chen FK.2014.Multiple Mesozoic magma processes formed the 240～185Ma composite Weishan pluton，South China:Evidence from geochronology，geochemistry，and Sr-Nd isotopes.International Geology Review，doi:10.1080/00206814.2014.905997

Ferry JM and Watson EB.2007.New thermodynamic models and revised calibrations for the Ti-in-zircon and Zr-in-rutile thermometers.Contributions to Mineralogy and Petrology，154(4):429-437

Griffin WL，Pearson NJ，Belousova E，Jackson SE，Van Achterbergh E，O’Reilly SY and Shee SR.2000.The Hf isotope composition of cratonic mantle:LAM-MC-ICPMS analysis of zircon megacrysts in kimberlites.Geochimica et Cosmochim Acta，64(1):133-147

Guffanti M，Clynne MA and Muffler L.1996.Thermal and mass implications of magmatic evolution in the Lassen volcanic region，California，and minimum constraints on basalt influx to the lower crust.Journal of Geophysical Research:Solid Earth (1978～2012)，101(B2):3003-3013

Hastie AR，Kerr AC，Pearce JA and Mitchell SF.2007.Classification of altered volcanic island arc rocks using immobile trace elements:Development of the Th-Co discrimination diagram.Journal of Petrology，48(12):2341-2357

Hoskin PWO and Black LP.2000.Metamorphic zircon formation by solid-state recrystallization of protolith igneous zircon.Journal of Metamorphic Geology，18(4):423-439

Huang Y，Zhu DC，Zhao ZD，Zhang LL，DePaolo D，Hu ZC，Yuan HL and Mo XX.2012.Petrogenesis and implication of the andesites at～113Ma in the Nagqu region in the northern Lhasa subterrane.Acta Petrologica Sinica，28(5):1603-1614 (in Chinese with English abstract)

Kang ZQ，Xu JF，Dong YH and Wang BD.2008.Cretaceous volcanic rocks of Zenong Group in north-middle Lhasa block:Products of southward subduction of the Slainajap ocean?Acta Petrologica Sinica，24(2):303-314 (in Chinese with English abstract)

Kang ZQ，Xu JF，Wang BD，Dong YH，Wang SQ and Chen JL.2009.Geochemistry of Cretaceous volcanic rocks of Duoni Formation in Northern Lhasa block:Discussion of tectonic setting.Earth Science，34(1):89-104 (in Chinese with English abstract)

Kapp P，Murphy MA，Yin A，Harrison TM，Ding L and Guo JH.2003.Mesozoic and Cenozoic tectonic evolution of the Shiquanhe area of western Tibet.Tectonics，22(4)，doi:10.1029/2001TC001332

Kay SM，Mpodozis C，Ramos VA and Munizaga F.1991.Magma source variations for Mid-Late Tertiary magmatic rocks associated with a shallowing subduction zone and a thickening crust in the central Andes (28° to 33° S).Geological Society of America Special Papers，265:113-138

Kemp AIS，Hawkesworth CJ，F(xiàn)oster GL，Paterson BA，Woodhead JD，Hergt JM，Gray CM and Whitehouse MJ.2007.Magmatic and crustal differentiation history of granitic rocks from Hf-O isotopes in zircon.Science，315(5814):980-983

Liang HY，Sun WD，Yu HX，Xie YW，Mo JH and Zhang YQ.2006.The age of ore-bearing porphyries in the Yulong copper belt and the formation process of porphyry Cu (Au)deposits.Mineral Deposits，25(Suppl.):415-418 (in Chinese with English abstract)

Li FQ，Liu W and Geng QR.2010.Zircon LA-ICP-MS U-Pb ages of the Mesozoic volcanic rocks in Nagqu area of Gangdise in Tibet and their geological significance.Acta Geoscientica Sinica，31(6):781-790(in Chinese with English abstract)

Li H，Ling MX，Li CY，Zhang H，Ding X，Yang XY，F(xiàn)an WM，Li YL and Sun WD.2012.A-type granite belts of two chemical subgroups in central eastern China:Indication of ridge subduction.Lithos，150:26-36

Li HY，Xu ZW，Lu XC，Chen LH，Liu SM，Yang XN，Zhang J，Li HC and Chen W.2008.Evolution of Mesozoic volcanic rocks in the Zouping basin，western Shandong Province:Constraints for mantle sources.Acta Petrologica Sinica，24(11):2537-2547 (in Chinese with English abstract)

Li XH，Long WG，Li QL，Liu Y，Zheng YF，Yang YH，Chamberlain KR，Wan DF，Guo CH，Wang XC and Tao H.2010.Penglai zircon megacrysts: A potential new working reference material for microbeam determination of Hf-O isotopes and U-Pb age.Geostandards and Geoanalytical Research，34(2):117-134

Liu YS，Hu ZC，Gao S，Gunther D，Xu J，Gao CG and Chen HH.2008.In situ analysis of major and trace elements of anhydrous minerals by LA-ICP-MS without applying an internal standard.Chemical Geology，257(1-2):34-43

Ludwig KR.2003.User’s Manual for ISOPLOT 3.00: A Geochronological Toolkit for Microsoft Excel.Berkeley:Berkeley Geochronology Center，Special Publication，4:1-70

Mo XX，Dong GC，Zhao ZD，Zhou S，Wang LL，Qiu RZ and Zhang FQ.2005.Spatial and temporal distribution and characteristics of granitoids in the Grangdese，Tibet and implication for crustal growth and evolution.Geological Journal of China Universities，11(3):281-290 (in Chinese with English abstract)

Pan GT，Mo XX，Hou ZQ，Zhu DC，Wang LQ，Li GM，Zhao ZD，Geng QR and Liao ZL.2006.Spatial-temporal framework of the Gangdese Orogenic Belt and its evolution.Acta Petrologica Sinica，22(3):521-533 (in Chinese with English abstract)

Pearce JA.1983.Role of the sub-continental lithosphere in magma genesis at active continental margins.In:Hawkesworth CJ and Morry MJ (eds.).Continental Basalts and Mantle Xenoliths.Nantwich，Cheshire:Shiva Publications，230-249

Pearce JA and Mei HJ.1988.Volcanic rocks of the 1985 Tibet Geotraverse:Lhasa to Golmud.Philosophical Transactions of the Royal Society of London，Series A，Mathematical and Physical Sciences，327(1594):169-201

Peng ZM，Geng QR，Zhang Z and Cong F.2011.Zircon LA-ICP-MS UPb ages and geochemical characteristics of rhyolite in Nagqu area，Tibet.Geological Bulletin of China，30(7):1050- 1059 (in Chinese with English abstract)

Price RC， Stewart RB， Woodhead JD and Smith IEM.1999.Petrogenesis of high-K arc magmas:Evidence from Egmont volcano，North Island，New Zealand.Journal of Petrology，40(1):167-197

Rapp RP，Shimizu N，Norman MD and Applegate GS.1999.Reaction between slab-derived melts and peridotite in the mantle wedge:Experimental constraints at 3.8 GPa.Chemical Geology，160(4):335-356

Richards JP and Villeneuve M.2002.Characteristics of late Cenozoic volcanism along the Archibarca lineament from Cerro Llullaillaco to Corrida de Cori，northwest Argentina.Journal of Volcanology and Geothermal Research，116(3-4):161-200

Rollinson HR.1993.Using Geochemical Data: Evaluation，Presentation，interpretation.New York:Longan Scientific and Technical，1-352

Rudnick R and Gao S.2003.Composition of the continental crust.Treatise on Geochemistry，3:1-64

Salters V and Stracke A.2004.Composition of the depleted mantle.Geochemistry， Geophysics， Geosystems， 5 (5 )， doi: 10.1029/2003GC000597

Shinjo R and Kato Y.2000.Geochemical constraints on the origin of bimodal magmatism at the Okinawa Trough，an incipient back-arc basin.Lithos，54(3-4):117-137

S?derlund U，Patchett PJ，Vervoort JD and Isachsen CE.2004.The176Lu decay constant determined by Lu-Hf and U-Pb isotope systematics of Precambrian mafic intrusions.Earth and Planet Science Letters，219(3-4):311-324

Sui QL，Wang Q，Zhu DC，Zhao ZD，Chen Y，Santosh M，Hu ZC，Yuan HL and Mo XX.2013.Compositional diversity of ca.110Ma magmatism in the northern Lhasa Terrane，Tibet:Implications for the magmatic origin and crustal growth in a continent-continent collision zone.Lithos，168-169:144-159

Sun SS and McDonough W.1989.Chemical and isotopic systematics of oceanic basalts:Implications for mantle composition and processes.In:Saunders AD and Norry MJ (eds.).Magmatism in the Ocean Basins.Geological Society，London，Special Publications，42(1):313-345

Tepper JH，Nelson BK，Bergantz GW and Irving AJ.1993.Petrology of the Chilliwack batholith，North Cascades，Washington:Generation of calc-alkaline granitoids by melting of mafic lower crust with variable water fugacity.Contributions to Mineralogy and Petrology，113(3):333-351

Tu XL，Zhang H，Deng WF，Ling MX，Liang HY，Liu Y and Sun WD.2011.Application of RESOlution in-situ laser ablation ICP-MS in trace element analyses.Geochimica，40(1):83-98 (in Chinese with English abstract)

Watson EB， Wark DA and Thomas JB.2006.Crystallization thermometers for zircon and rutile.Contributions to Mineralogy and Petrology，151(4):413-433

Wilson M.1993.Magmatism and the geodynamics of basin formation.Sediment Geology，86(1-2):5-29

Winchester JA and Floyd PA.1977.Geochemical discrimination of different magma series and their differentiation products using immobile elements.Chemical Geology，20:325-343

Wu H，Li C，Hu PY，F(xiàn)an JJ，Zhang HY and Li J.2013.The discovery of Qushenla volcanic rocks in Tasepule area of Nyima County，Tibet，and its geological significance.Geological Bulletin of China，32(7):1014-1026 (in Chinese with English abstract)

Wu FY，Yang YH，Xie LW，Yang JH and Xu P.2006.Hf isotopic compositions of the standard zircons and baddeleyites used in U-Pb geochronology.Chemical Geology，234(1-2):105-126

Wu FY，Li XH，Yang JH and Zheng YF.2007.Discussions on the petrogenesis of granites.Acta Petrologica Sinica，23(6):1217-1238 (in Chinese with English abstract)

Xu RH，Sch?rer U and Allègre CJ.1985.Magmatism and metamorphism in the Lhasa Block (Tibet):A geochronological study.Journal of Petrology，93(1):41-57

Zhang H，Ling MX，Liu YL，Tu XL，Wang FY，Li CY，Liang HY，Yang XY，Arndt NT and Sun WD.2013.High oxygen fugacity and slab melting linked to Cu mineralization:Evidence from Dexing porphyry copper deposits，southeastern China.Journal of Geology，121(3):289-305

Zhang LL，Zhu DC，Zhao ZD，Dong GC，Mo XX，Guan Q，Liu M and Liu MH.2010.Petrogenesis of magmatism in the Baerda region of northern Gangdese， Tibet: Constraints from geochemistry，geochronology and Sr-Nd-Hf isotopes.Acta Petrologica Sinica，26(6):1871-1888 (in Chinese with English abstract)

Zhang LL，Zhu DC，Zhao ZD，Liao ZL，Wang LQ and Mo XX.2011.Early Cretaceous granitoids in Xainza，Tibet:Evidence of slab break-off.Acta Petrologica Sinica，27 (7):1938- 1948 (in Chinese with English abstract)

Zhang XQ，Zhu DC，Zhao ZD，Wang LQ，Huang JC and Mo XX.2010.Petrogenesis of the Nixiong pluton in Coqen，Tibet and its potential significance for the Nixiong Fe-rich mineralization.Acta Petrologica Sinica，26(6):1793-1804 (in Chinese with English abstract)

Zhu DC，Pan GT，Wang LQ，Mo XX，Zhao ZD，Zhou CY，Liao ZL，Dong GC and Yuan SH.2008.Tempo-spatial variations of Mesozoic magmatic rocks in the Gangdese belt，Tibet，China，with a discussion of geodynamic setting-related issues.Geological Bulletin of China，27(9):1535-1550 (in Chinese with English abstract)

Zhu DC，Mo XX，Wang LQ，Zhao ZD，Niu YL，Zhou CY and Yang YH.2009a.Petrogenesis of highly fractionated I-type granites in the Zayu area of eastern Gangdese，Tibet:Constraints from zircon U-Pb geochronology，geochemistry and Sr-Nd-Hf isotopes.Science in China (Series D)，52(9):1223-1239

Zhu DC，Zhao ZD，Pan GT，Lee HY，Kang ZQ，Liao ZL，Wang LQ，Li GM，Dong GC and Liu B.2009b.Early Cretaceous subductionrelated adakite-like rocks of the Gangdese Belt，southern Tibet:Products of slab melting and subsequent melt-peridotite interaction?Journal of Asian Earth Sciences，34(3):298-309

Zhu DC，Zhao ZD，Niu YL，Mo XX，Chung SL，Hou ZQ，Wang LQ and Wu FY.2011.The Lhasa Terrane:Record of a microcontinent and its histories of drift and growth.Earth Planet Science Letters，301(1-2):241-255

Zhu DC，Zhao ZD，Niu YL，Dilek Y，Hou ZQ and Mo XX.2013.The origin and pre-Cenozoic evolution of the Tibetan Plateau.Gondwana Research，23(4):1429-1454

Zorpi MJ，Coulon C and Orsini JB.1991.Hybridization between felsic and mafic magmas in calc-alkaline granitoids:A case study in northern Sardinia，Italy.Chemical Geology，92(1-3):45-86

附中文參考文獻

常青松，朱弟成，趙志丹，董國臣，莫宣學，劉勇勝，胡兆初.2011.西藏羌塘南緣熱那錯早白堊世流紋巖鋯石U-Pb 年代學和Hf 同位素及其意義.巖石學報，7(7):2034-2044

陳越，朱弟成，趙志丹，張亮亮，劉敏，于楓，管琪，莫宣學.2010.西藏北岡底斯巴木錯安山巖的年代學、地球化學及巖石成因.巖石學報，26(7):2193-2206

黃玉，朱弟成，趙志丹，張亮亮，DePaolo D，胡兆初，袁洪林，莫宣學.2012.西藏北部拉薩地塊那曲地區(qū)約113Ma 安山巖巖石成因與意義.巖石學報，28(5):1603-1614

康志強，許繼峰，董彥輝，王保弟.2008.拉薩地塊中北部白堊紀則弄群火山巖:Slainajap 洋南向俯沖的產(chǎn)物?巖石學報，24(2):303-314

康志強，許繼峰，王保弟，董彥輝，王樹慶，陳建林.2009.拉薩地塊北部白堊紀多尼組火山巖的地球化學:形成的構(gòu)造環(huán)境.地球科學，34(1):89-104

梁華英，孫衛(wèi)東，喻亨祥，謝應(yīng)雯，莫濟海，張玉泉.2006.西藏東緣玉龍斑巖銅礦帶含礦巖體時代及斑巖銅金礦床形成研究.礦床地質(zhì)，25(增刊):415-418

李奮其，劉偉，耿全如.2010.西藏岡底斯帶那曲地區(qū)中生代火山巖的LA-ICP-MS 鋯石U-Pb 年齡和地質(zhì)意義.地球?qū)W報，31(6):781-790

李海勇，徐兆文，陸現(xiàn)彩，陳立輝，劉蘇明，楊小男，張軍，李紅超，陳偉.2008.魯西鄒平盆地中生代火山巖的演化:對地幔源區(qū)的約束.巖石學報，24(11):2537-2547

莫宣學，董國臣，趙志丹，周肅，王亮亮，邱瑞照，張風琴.2005.西藏岡底斯帶花崗巖的時空分布特征及地殼生長演化信息.高校地質(zhì)學報，11(3):281-290

潘桂棠，莫宣學，侯增謙，朱弟成，王立全，李光明，趙志丹，耿全如，廖忠禮.2006.岡底斯造山帶的時空結(jié)構(gòu)及演化.巖石學報，22(3):521-533

彭智敏，耿全如，張璋，叢峰.2011.西藏那曲地區(qū)流紋巖LA-ICPMS 鋯石U-Pb 年齡和地球化學特征.地質(zhì)通報，30(7):1050-1059

涂湘林，張紅，鄧文峰，凌明星，梁華英，劉穎，孫衛(wèi)東.2011.RESOlution 激光剝蝕系統(tǒng)在微量元素原位微區(qū)分析中的應(yīng)用.地球化學，40(1):83-98

吳浩，李才，胡培遠，范建軍，張紅雨，李嬌.2013.西藏尼瑪縣塔色普勒地區(qū)去申拉組火山巖的發(fā)現(xiàn)及其地質(zhì)意義.地質(zhì)通報，32(7):1014-1026

吳福元，李獻華，楊進輝，鄭永飛.2007.花崗巖成因研究的若干問題.巖石學報，23(6):1217-1238

張亮亮，朱弟成，趙志丹，董國臣，莫宣學，管琪，劉敏，劉美華.2010.西藏北岡底斯巴爾達地區(qū)巖漿作用的成因:地球化學、年代學及Sr-Nd-Hf 同位素約束.巖石學報，26(6):1871-1888

張亮亮，朱弟成，趙志丹，廖忠禮，王立全，莫宣學.2011.西藏申扎早白堊世花崗巖類:板片斷離的證據(jù).巖石學報，27(7):1938-1948

張曉倩，朱弟成，趙志丹，王立全，黃建村，莫宣學.2010.西藏措勤尼雄巖體的巖石成因及其對富Fe 成礦作用的潛在意義.巖石學報，26(6):1793-1804

朱弟成，潘桂棠，王立全，莫宣學，趙志丹，周長勇，廖忠禮，董國臣，袁四化.2008.西藏岡底斯帶中生代巖漿巖的時空分布和相關(guān)問題的討論.地質(zhì)通報，27(9):1535-1550