• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    西藏謝通門縣切瓊地區(qū)典中組流紋巖鋯石U-Pb年齡、Hf同位素及地球化學(xué)特征

    2022-03-04 01:41:20閆茂強(qiáng)沈志遠(yuǎn)魏俊浩毛國(guó)正鄧永明
    大地構(gòu)造與成礦學(xué) 2022年1期

    曾 成, 閆茂強(qiáng), 沈志遠(yuǎn), 高 強(qiáng), 魏俊浩*, 毛國(guó)正, 鄧永明

    西藏謝通門縣切瓊地區(qū)典中組流紋巖鋯石U-Pb年齡、Hf同位素及地球化學(xué)特征

    曾 成1, 閆茂強(qiáng)1, 沈志遠(yuǎn)1, 高 強(qiáng)1, 魏俊浩1*, 毛國(guó)正2, 鄧永明2

    (1. 中國(guó)地質(zhì)大學(xué)(武漢) 資源學(xué)院, 湖北 武漢 430074; 2. 西藏工程勘查施工(集團(tuán))有限公司, 西藏 拉薩 850000)

    西藏中拉薩地塊切瓊地區(qū)廣泛發(fā)育的林子宗群火山巖被認(rèn)為是印度大陸與亞洲大陸碰撞火山作用的響應(yīng), 對(duì)于揭示大陸碰撞的時(shí)限和過(guò)程具有重要的意義。本文對(duì)區(qū)內(nèi)流紋巖進(jìn)行了LA-ICP-MS鋯石U-Pb年齡、Hf同位素及全巖元素地球化學(xué)研究, 結(jié)果顯示: 研究區(qū)流紋巖年齡為67.1±0.1 Ma, 符合林子宗群典中組火山巖年齡范圍。該流紋巖具有富Si、高K、低Ti、過(guò)鋁質(zhì)特征; 富集輕稀土元素, 虧損重稀土元素, (Gd/Yb)N=0.57~1.16, 普遍出現(xiàn)負(fù)Eu異常(δEu=0.05~0.08), 相對(duì)富集大離子親石元素(Pb、Rb、Th、U), 虧損高場(chǎng)強(qiáng)元素(P、Ti), 并且明顯虧損 Ba和Sr, 屬于高鉀鈣堿性過(guò)鋁質(zhì)系列, 具有陸源弧火山巖的特點(diǎn); 鋯石176Hf/177Hf值在0.282673~0.282804之間,Hf()在?2.1~+2.3之間, 平均為?0.6,2DM變化于875~1117 Ma之間。綜合分析表明, 切瓊地區(qū)林子宗群典中組流紋巖形成于俯沖構(gòu)造環(huán)境, 源于古老地殼物質(zhì)再循環(huán), 且混染了少量幔源物質(zhì)。對(duì)比區(qū)域上相關(guān)礦床, 結(jié)合切瓊地區(qū)已發(fā)現(xiàn)礦點(diǎn), 表明切瓊地區(qū)典中組火山巖具有一定尋找淺成低溫?zé)嵋盒汀⑽◣r型礦床的潛力。

    林子宗群; 典中組; 流紋巖; 晚白堊世; 中拉薩地塊

    0 引 言

    印度?亞洲大陸碰撞作為地球上重大的地質(zhì)事件之一, 一直是國(guó)內(nèi)外地學(xué)工作者關(guān)注的焦點(diǎn)。近30年來(lái), 眾多地質(zhì)學(xué)家在沉積地層(DeCelles et al., 2004, 2014; Zhu et al., 2005; Najman et al., 2010; Cai et al., 2011; Roddaz et al., 2011; Wang et al., 2011; Hu et al., 2012a, 2016; Orme et al., 2015)、古地磁(Dupont-Nivet et al., 2010; Yi et al., 2011; Meng et al., 2012; Lippert et al., 2014; Huang et al., 2015)和超高壓變質(zhì)(O’brien et al., 2001; Leech et al., 2005;Donaldson et al., 2013)等方面研究對(duì)印度?亞洲大陸碰撞的認(rèn)識(shí)取得了巨大的進(jìn)展(Wang et al., 2014)。印度?亞洲大陸碰撞經(jīng)歷了從弧?陸碰撞到陸?陸碰撞復(fù)雜的演化過(guò)程(潘桂棠等, 2006; 侯增謙等, 2006a, 2006b; 莫宣學(xué)和潘桂棠, 2006; 朱弟成等, 2008; Zhu et al., 2013; 丁林等, 2017; Liu et al., 2018; Colleps et al., 2020), 從洋殼的俯沖到后續(xù)碰撞造山, 并伴隨著大量的巖漿作用, 關(guān)于巖石成因和形成機(jī)制的研究對(duì)認(rèn)識(shí)和了解造山作用、深部動(dòng)力學(xué)過(guò)程, 以及限定初始碰撞時(shí)限具有重要意義。林子宗群火山巖作為中、新生代巖漿作用重要的產(chǎn)物, 廣泛分布于拉薩地塊南部, 西起獅泉河, 東至拉薩東, 構(gòu)成一條長(zhǎng)達(dá)1000 km的火山巖帶, 總厚約5000 m。自下而上, 林子宗群分為典中組、年波組和帕那組(莫宣學(xué)等, 2003; 董國(guó)臣等, 2005; 侯增謙等, 2006b), 具有由鈣堿性的陸緣弧火山巖經(jīng)鈣性?鈣堿性的陸緣弧?陸內(nèi)火山巖到高鉀鈣堿性?鉀玄巖系列的板內(nèi)火山巖過(guò)渡的地球化學(xué)特征(Mo et al., 2007, 2008; Lee et al., 2009), 一定程度上記錄了新特提斯洋俯沖消減過(guò)渡到印度?亞洲大陸碰撞的過(guò)程(莫宣學(xué)等, 2003; Ding, 2003)。處于林子宗群底部的典中組火山巖, 與下伏地層為區(qū)域性角度不整合接觸, 該不整合關(guān)系被認(rèn)為可能代表了印度大陸與亞洲大陸之間的碰撞事件。因此對(duì)典中組火山巖的研究, 可以更好地揭示該地區(qū)碰撞開始的時(shí)限與過(guò)程(莫宣學(xué)等, 2003, 2007, 2009; 莫宣學(xué)和潘桂棠, 2006; Lee et al., 2009, 2012)。

    近年來(lái), 林子宗群典中組火山巖的研究區(qū)域, 由早期的拉薩地塊東部林周地區(qū)(90°E以東)(李皓揚(yáng)等, 2007; He et al., 2007; Lee et al., 2012; 楊輝等, 2013; Huang et al., 2015; Zhu et al., 2015; 陳貝貝等, 2016; 唐攀等, 2018; Yan et al., 2019; 韋乃韶等, 2019), 逐漸擴(kuò)大到拉薩地塊的中部(84°E~90°E)(胡新偉等, 2007; 李再會(huì)等, 2008; 梁銀平等, 2010; 謝冰晶等, 2013; 張小強(qiáng), 2013; 丁帥等, 2017; 付燕剛, 2017; 李勇等, 2018; 周鵬等, 2019)和西部(84°E以西)(王喬林, 2011; Jiang et al., 2018; 姜軍勝, 2018; 曹延等, 2019), 進(jìn)一步完善了林子宗群典中組火山巖系的年代學(xué)框架、巖石成因及構(gòu)造背景體系。另一方面, 隨著找礦工作的進(jìn)行以及成巖成礦研究的不斷深入, 在林子宗群典中組火山巖中發(fā)現(xiàn)了一批多金屬礦床, 比如斯弄多淺成低溫?zé)嵋盒虯g-Pb-Zn礦床(唐菊興等, 2016; 付燕剛等, 2017)、仁堆Pb-Zn礦床(張小強(qiáng), 2013)、諾倉(cāng)Pb-Zn(Cu)礦床、北納Pb-Zn(Cu)礦床、打加錯(cuò)Ag-Pb-Zn礦床等(姜軍勝, 2018), 呈現(xiàn)出林子宗群典中組火山巖較好的找礦潛力, 拓展了對(duì)其研究的意義。

    本次對(duì)研究程度較低的拉薩地塊中段切瓊地區(qū)林子宗群典中組火山巖進(jìn)行了詳細(xì)的剖面觀察, 并針對(duì)區(qū)內(nèi)廣泛出露的流紋巖進(jìn)行巖相學(xué)觀察、LA-ICP-MS鋯石U-Pb定年、鋯石Hf同位素及巖石地球化學(xué)研究, 同時(shí)與拉薩地塊東段、中段和西段出露的典型典中組火山巖進(jìn)行對(duì)比, 探討其形成時(shí)代、巖石成因及構(gòu)造環(huán)境, 豐富和完善林子宗群典中組火山巖的研究, 為討論典中組火山巖區(qū)域性變化特征與成礦潛力提供基礎(chǔ)資料。

    1 地質(zhì)背景與巖石特征

    拉薩地塊是指位于青藏高原南部的印度河?雅魯藏布江縫合帶(IYZSZ)與北部的班公湖?怒江縫合帶(BNSZ)之間近東西向的狹長(zhǎng)地域。根據(jù)基底巖石和沉積蓋層的差異, 以沙莫勒?麥拉?洛巴堆?米拉山斷裂(SMLMF)和獅泉河?永珠?納木錯(cuò)?嘉黎蛇綠混雜巖帶(SYNJF)為界, 由南向北將拉薩地塊劃分為南拉薩地塊、中拉薩地塊和北拉薩地塊(圖1b; Zhu et al., 2011, 2013)。

    切瓊地區(qū)大地構(gòu)造位置處于中拉薩地塊中段南部(圖1b), 措勤?納木錯(cuò)初始弧間裂谷盆地的西段。林子宗群火山巖在切瓊地區(qū)出露面積約893.3 km2, 不整合覆蓋于石炭系永珠組(C1-2)、上石炭統(tǒng)?下二疊統(tǒng)拉嘎組(C2~P1)、二疊系昂杰組(P1)之上, 其上被烏郁群(N2)及第四系(Q)覆蓋。該火山巖大部分為典中組火山巖, 僅在西南角出露少量年波組火山巖(圖2)。通過(guò)詳細(xì)的野外剖面觀察, 切瓊地區(qū)典中旋回由18個(gè)韻律組成, 巖性及厚度如圖3所示, 典中組火山巖系在該區(qū)出露厚達(dá)766 m, 上覆年波組為含角礫晶屑凝灰?guī)r, 下伏昂杰組為石英砂巖, 因剖面位置昂杰組出露較少, 在圖2中未體現(xiàn)。典中組火山活動(dòng)多以噴溢?爆發(fā)為主, 也有噴(爆)發(fā)?噴溢?沉積方式。巖石成分總體上具有從早到晚由中性向中酸性(酸性)演化的特征?;鹕交顒?dòng)較強(qiáng)烈, 爆發(fā)指數(shù)為15%~23%, 早期爆發(fā)較強(qiáng)烈, 中晚期以噴溢為主, 屬陸相噴發(fā)產(chǎn)物。整個(gè)切瓊地區(qū), 空間上, 火山巖東西方向上巖性及厚度變化不大, 南北方向上變化較大, 以爆發(fā)相的酸性巖石為主, 溢流相的巖石次之, 并以沉積夾層的相對(duì)發(fā)育為特征; 時(shí)間上, 其下部以爆發(fā)相和溢流相互層為特征, 主要巖石類型有英安質(zhì)晶屑凝灰?guī)r, 流紋質(zhì)含角礫晶屑凝灰?guī)r、晶屑凝灰?guī)r、安山巖、英安巖、流紋巖; 中部以發(fā)育爆發(fā)相和沉積相為特征, 其巖石類型有火山角礫巖、英安巖或流紋質(zhì)含角礫晶屑凝灰?guī)r、英安質(zhì)玻屑凝灰熔巖、含礫雜砂巖、雜砂巖、細(xì)砂巖; 上部以爆發(fā)相和溢流相為主, 分為上、下兩段, 下段以火山碎屑巖為主, 其巖石類型有英安質(zhì)火山角礫巖、英安質(zhì)晶屑玻屑凝灰?guī)r、凝灰?guī)r和流紋質(zhì)晶屑凝灰熔巖, 上段以熔巖為主, 次為碎屑巖和沉積巖, 其巖石類型有似球粒流紋巖、流紋質(zhì)晶屑凝灰?guī)r、硅質(zhì)巖, 自下而上火山活動(dòng)具有弱→強(qiáng)→弱的演化特征。野外宏觀觀察結(jié)果表明, 典中組火山巖多呈層狀、似層狀產(chǎn)出, 沉積夾層較少, 表現(xiàn)為火山物質(zhì)的正常堆積, 地層結(jié)構(gòu)主要為熔巖與火山碎屑巖互層特征。根據(jù)礦調(diào)報(bào)告(西藏自治區(qū)地質(zhì)調(diào)查院, 2016), 在切瓊地區(qū)林子宗群火山巖中發(fā)現(xiàn)鉛鋅礦、銅多金屬礦、銀礦等多個(gè)礦點(diǎn), 展現(xiàn)出較好的找礦潛力。其中, 查嘎勒姐礦點(diǎn)賦存在典中組火山角礫巖裂隙中, 色拉鉛鋅礦化點(diǎn)賦存在殘留于典中組火山角礫巖中的泥質(zhì)粉砂巖內(nèi), 均與典中組火山巖密切相關(guān)。

    圖1 西藏拉薩地塊林子宗群火山巖及部分典中組年齡分布簡(jiǎn)圖(據(jù)Zhu et al., 2011修改)

    1. 第四系; 2. 烏郁群碎屑巖段; 3. 年波組火山巖; 4. 典中組火山巖; 5. 敵布錯(cuò)組; 6. 二疊系; 7. 拉嘎組; 8. 永珠組; 9. 始新世二長(zhǎng)花崗巖; 10. 晚三疊世二長(zhǎng)花崗巖; 11. 斷層; 12. 地質(zhì)界線; 13. 不整合界線; 14. 磁鐵礦點(diǎn); 15. 鉛鋅礦點(diǎn); 16. 銅多金屬礦點(diǎn); 17. 銀礦點(diǎn); 18. 實(shí)測(cè)剖面位置; 19. 采樣點(diǎn)。

    圖3 切瓊地區(qū)典中組火山巖柱狀圖

    用于鋯石U-Pb定年、Hf同位素及巖石地球化學(xué)分析的樣品均為流紋巖, 采樣位置見圖2。流紋巖, 主要分布于研究區(qū)東南部, 巖石呈淺灰、淺黃色(圖4b), 球粒結(jié)構(gòu)(圖4c), 塊狀構(gòu)造、流紋構(gòu)造(圖4d), 主要由變余長(zhǎng)英質(zhì)球粒(70%)、石英(10%)、蝕變長(zhǎng)石(15%)、褐鐵礦(2%)和少量蝕變鉀長(zhǎng)石斑晶等組成。變余長(zhǎng)英質(zhì)球粒呈球粒狀, 粒度多為0.2 mm左右, 具長(zhǎng)石泥霧狀蝕變, 部分球粒變質(zhì)重結(jié)晶成顯微文象長(zhǎng)石狀, 但仍顯示放射狀結(jié)構(gòu)。石英呈不規(guī)則粒狀, 多數(shù)粒度小于0.1 mm, 充填于球粒之間, 少數(shù)石英粒度稍大, 達(dá)0.2 mm×0.4 mm。蝕變長(zhǎng)石呈不規(guī)則粒狀, 粒度多小于0.1 mm, 具泥霧狀蝕變, 偶見聚片雙晶, 既有斜長(zhǎng)石, 也有鉀長(zhǎng)石。褐鐵礦呈紅褐色, 褐黑色, 不規(guī)則團(tuán)塊狀, 團(tuán)塊粒度達(dá)0.6 mm,為原巖磁鐵礦風(fēng)化而成。蝕變鉀長(zhǎng)石斑晶(圖4d)呈自形?半自形板狀, 粒度大者達(dá)1.6 mm×2.4 mm, 發(fā)育泥狀高嶺石化、綠簾石化。

    2 分析方法

    2.1 LA-ICP-MS鋯石U-Pb定年和Hf同位素測(cè)試

    本文用于鋯石U-Pb定年的樣品為流紋巖(-Zr), 巖石新鮮且表面經(jīng)過(guò)打磨清理。鋯石分選由廊坊巖拓地質(zhì)服務(wù)有限公司完成, 挑選晶形完好、未蝕變的鋯石顆粒制成樣品靶。鋯石陰極發(fā)光圖像拍攝在武漢上譜分析科技有限責(zé)任公司完成, 儀器為高真空掃描電子顯微鏡(JSM-IT100), 配備有GATAN MINICL系統(tǒng)。

    鋯石U-Pb同位素定年和微量元素含量在武漢上譜分析科技有限責(zé)任公司利用LA-ICP-MS同時(shí)分析完成。詳細(xì)的儀器參數(shù)和分析流程見Zonget al., (2010)。本次分析的激光束斑為32 μm, 頻率為5 Hz。U-Pb同位素定年和微量元素含量處理中采用鋯石標(biāo)準(zhǔn)91500和玻璃標(biāo)準(zhǔn)物質(zhì)NIST610作外標(biāo)分別進(jìn)行同位素和微量元素分餾校正。數(shù)據(jù)的離線處理采用軟件ICPMSDataCal(Liuet al., 2008, 2010)完成。鋯石U-Pb年齡諧和圖繪制和加權(quán)平均年齡計(jì)算采用Isoplot/Ex_ver3(Ludwig, 2003)軟件完成。

    鋯石原位微區(qū)Hf同位素測(cè)試在武漢上譜分析科技有限責(zé)任公司利用Neptune型多接收等離子體質(zhì)譜儀和Geolas Pro型激光剝蝕系統(tǒng)聯(lián)用的方法完成, 詳細(xì)的分析流程見Hu et al. (2012b)。Hf同位素測(cè)試點(diǎn)位置與U-Pb測(cè)試點(diǎn)位置一致, 測(cè)試束斑直徑為32 μm。

    (a) 巖石露頭; (b) 流紋巖手標(biāo)本照片; (c) 流紋巖球粒結(jié)構(gòu)(正交偏光); (d) 流紋構(gòu)造、蝕變鉀長(zhǎng)石斑晶(正交偏光)。礦物代號(hào): Qtz. 石英; Kfs. 鉀長(zhǎng)石。

    2.2 巖石地球化學(xué)測(cè)試

    樣品全巖主量、微量和稀土元素測(cè)試均在澳實(shí)分析檢測(cè)(廣州)有限公司完成。主量元素分析釆用ME-XRF26d X射線熒光光譜儀熔融法完成, 微量和稀土元素分析采用M61-MS81電感耦合等離子體質(zhì)譜法(ICP-MS)完成, 主量元素分析精度和準(zhǔn)確度優(yōu)于5%, 微量和稀土元素分析精度和準(zhǔn)確度優(yōu)于10%。

    3 分析結(jié)果

    3.1 LA-ICP-MS鋯石U-Pb年齡

    切瓊地區(qū)典中組流紋巖樣品(-Zr)LA-ICP-MS鋯石U-Pb年齡測(cè)試結(jié)果見表1。鋯石均為無(wú)色透明, 顆粒自形程度較好, 多為自形?半自形晶體, 粒徑在80~200 μm之間, 長(zhǎng)寬比為1∶1~3∶1。從CL圖像(圖5)可見, 所有鋯石均發(fā)育典型的巖漿振蕩環(huán)帶(Wildeet al., 2001)。數(shù)據(jù)顯示Th/U值為0.56~2.70, 均值為1.12, 具典型的巖漿鋯石特征(吳元保和鄭永飛, 2004)。25個(gè)測(cè)試點(diǎn)中除2個(gè)不諧和點(diǎn)(點(diǎn)4和 23)外, 其余23個(gè)測(cè)試點(diǎn)206Pb/238U年齡分布在67.0~67.6 Ma之間, 在諧和圖上測(cè)試點(diǎn)均位于諧和線上或諧和線附近(圖6a),206Pb/238U加權(quán)平均年齡為67.1±0.1 Ma (MSWD=0.20)。因此, 該流紋巖的形成年齡為67.1± 0.1 Ma, 為晚白堊世巖漿作用的產(chǎn)物。

    3.2 全巖地球化學(xué)特征

    為探討切瓊地區(qū)林子宗群火山巖的地球化學(xué)特征, 對(duì)6件流紋巖樣品進(jìn)行了主量、微量和稀土元素分析, 分析結(jié)果見表2。本次測(cè)試樣品均采自地表新鮮露頭且不含裂隙與細(xì)脈, 全巖燒失量均小于3%, 表明測(cè)試的巖石樣品不受后期風(fēng)化作用, 或受后期風(fēng)化作用以及流體作用的影響較小。

    3.2.1 主量元素

    流紋巖樣品的SiO2含量普遍較高, 在75.93%~ 79.75%之間, 平均含量為78.18%; 全堿K2O+Na2O含量為5.82%~7.99%, 平均含量為7.16%; 其中K2O含量為3.23%~5.10%, 平均含量為4.18%; Na2O含量為2.18%~3.52%, 平均含量為2.98%, K2O/Na2O值為0.96~2.01, 平均值為1.45。在火山巖TAS圖中, 6件樣品全部投在流紋巖區(qū)域, 且處于亞堿性范圍內(nèi)(圖7a);在K2O-SiO2圖中, 6件樣品中主要投在高鉀鈣堿性系列區(qū)域, 部分處于高鉀鈣堿性系列與鈣堿性系列界線附近(圖7c)。Al2O3含量相對(duì)較高, 在11.58%~ 13.30%之間, 平均含量為12.41%; A/CNK值為1.19~ 1.61, 均大于1, 平均值為1.32, 在鋁飽和指數(shù)圖解中, 6件樣品全部投在過(guò)鋁質(zhì)區(qū)域(圖7d)。因此, 該火山巖樣品為過(guò)鋁質(zhì)鈣堿性?高鉀鈣堿性流紋巖, 與已報(bào)道的拉薩地塊中段典中組火山巖特征基本一致。

    表1 切瓊地區(qū)流紋巖LA-ICP-MS鋯石U-Pb同位素分析結(jié)果

    圖5 切瓊地區(qū)流紋巖鋯石陰極發(fā)光圖像(括號(hào)內(nèi)為εHf(t)值, 年齡單位: Ma)

    圖6 切瓊地區(qū)流紋巖鋯石U-Pb年齡諧和圖

    表2 切瓊地區(qū)典中組流紋巖全巖地球化學(xué)分析結(jié)果

    續(xù)表2:

    注: 主量元素(%), 微量和稀土元素(×10?6), LOI為燒失量, A/NK=Al2O3/(Na2O+K2O)(摩爾比), A/CNK=Al2O3/(CaO+Na2O+K2O)(摩爾比), 戈蒂里指數(shù)=(Al2O3?Na2O)/TiO2, 里特曼指數(shù)=(Na2O+K2O)2/(SiO2?43), δEu=2×EuN/(SmN+GdN), 下標(biāo)N表示球粒隕石標(biāo)準(zhǔn)化后數(shù)值, 標(biāo)準(zhǔn)化值據(jù)Sun and McDonough, 1989。

    3.2.2 微量元素

    原始地幔標(biāo)準(zhǔn)化微量元素蛛網(wǎng)圖呈現(xiàn)多峰多谷的右傾模式(圖8a), 可以看出, 該流紋巖相對(duì)富集Pb、Rb、Th、K、U等元素, 虧損Ba、Ti、P、Sr等元素。Zr/TiO2-Nb/Y圖解也進(jìn)一步確定該火山巖為流紋巖(圖7b)。

    3.2.3 稀土元素

    該流紋巖樣品稀土元素總量為73.5×10?6~ 172×10?6, 整體偏高且變化較大。其中輕稀土元素(LREE)含量為47.3×10?6~149×10?6, 重稀土元素(HREE)含量為19.6×10?6~26.2×10?6, LREE/HREE值介于1.80~6.90之間, (La/Yb)N為1.31~5.51, 具明顯的輕重稀土元素分異, 相對(duì)富集輕稀土元素, 虧損重稀土元素; (Gd/Yb)N為0.574~0.960, 表明重稀土元素未發(fā)生分餾。球粒隕石標(biāo)準(zhǔn)化稀土元素分布型勢(shì)圖(圖8b)呈明顯的右傾模式, 具有顯著的負(fù)Eu異常(δEu=0.05~0.07)等特點(diǎn)。

    3.3 鋯石Hf同位素

    切瓊地區(qū)典中組流紋巖鋯石Hf同位素測(cè)試結(jié)果見表3。樣品中23顆鋯石176Hf/177Hf值在0.282673~ 0.282804之間, 平均為0.282716;Hf()在?2.1~2.3之間, 平均為?0.6, 大部分小于0(圖9b、c); 對(duì)應(yīng)的二階段模式年齡(2DM)變化于875~1117 Ma之間, 平均為1037 Ma(圖9a)。

    4 討 論

    4.1 火山巖形成時(shí)代

    切瓊地區(qū)典中組流紋巖加權(quán)平均年齡為67.1± 0.1 Ma, 其代表了該火山巖的形成時(shí)代, 與拉薩地塊中段典中組火山巖的形成年齡(70.89~57.42 Ma)相符。

    在林子宗群火山巖研究早期, 一批學(xué)者對(duì)林周盆地林子宗群火山巖進(jìn)行系統(tǒng)的40Ar-39Ar同位素測(cè)年, 給出了林子宗群典中組火山巖與下伏地層間不整合形成最晚時(shí)限為65 Ma左右(莫宣學(xué)等, 2003; 周肅等, 2004; 董國(guó)臣等, 2005)。目前已有數(shù)據(jù)表明, 整個(gè)拉薩地塊林子宗群典中組火山巖年齡為70.89~ 57.42 Ma(圖10), 其中拉薩地塊東段典中組火山巖年齡為70.89~58.30 Ma(李皓揚(yáng)等, 2007; He et al., 2007; 楊輝等, 2013; Huang et al., 2015; Zhu et al., 2015; 陳貝貝等, 2016; 唐攀等, 2018; Yan et al., 2019; 韋乃韶等, 2019; 周鵬等, 2019), 拉薩地塊中段典中組火山巖年齡為70.70~57.42 Ma(梁銀平等, 2010; 謝冰晶等, 2013; 張小強(qiáng), 2013; 付燕剛, 2017; 丁帥等, 2017; 李勇等, 2018; 陳蘭樸等, 2019), 拉薩地塊西段典中組火山巖年齡為70.80~64.25 Ma(王喬林, 2011; 姜軍勝, 2018; 曹延等, 2019)。綜合前人資料表明, 拉薩地塊林子宗群典中組火山巖的活動(dòng)時(shí)間為晚白堊世至始新世, 且具有東段、中段比西段持續(xù)時(shí)間長(zhǎng)的特點(diǎn)。

    4.2 巖石成因

    酸性火山巖是巖漿演化最后階段的產(chǎn)物, 是活動(dòng)大陸邊緣、大陸裂谷以及大火山巖省等大陸火成作用的重要組成部分(McCulloch et al., 1994; Feeley and Dungan, 1996; Pankhurst et al., 2000; Riley et al., 2001; Miller and Harris, 2007), 其成因研究對(duì)認(rèn)識(shí)大陸地殼演化過(guò)程具有重要意義(Annen et al., 2006)。研究區(qū)典中組流紋巖為過(guò)鋁質(zhì)鈣堿性?高鉀鈣堿性流紋巖, SiO2含量普遍較高(75.93%~79.75%, 平均值為78.18%), 為典型的酸性巖石。

    從巖石地球化學(xué)特征可以看出, 該樣品過(guò)鋁質(zhì)特征明顯, 相對(duì)富集大離子親石元素Rb、Th、U和K, 虧損高場(chǎng)強(qiáng)元素Ti、P和大離子親石元素Ba、Sr。輕重稀土元素分異明顯, 相對(duì)富集輕稀土元素, 虧損重稀土元素, 負(fù)Eu異常特征顯著。通過(guò)與拉薩地塊中段已報(bào)道的謝通門(李再會(huì)等, 2008)、朱諾(梁銀平等, 2010)、廳宮(周鵬等, 2019)、孔隆?丁仁勒(謝冰晶等, 2013)、斯弄多(丁帥等, 2017)等地區(qū)典中組火山巖的巖石地球化學(xué)特征對(duì)比分析, 切瓊地區(qū)流紋巖與拉薩地塊中段其他地區(qū)典中組火山巖具有一致的陸緣弧火山巖的特點(diǎn)。

    數(shù)據(jù)來(lái)源: 謝通門地區(qū)據(jù)李再會(huì)等, 2008; 朱諾地區(qū)據(jù)梁銀平等, 2010; 孔隆至丁仁勒地區(qū)據(jù)謝冰晶等, 2013; 南木林地區(qū)據(jù)丁帥等, 2017; 廳宮地區(qū)據(jù)周鵬等, 2019。底圖a圖據(jù)Middlemost, 1994; b圖據(jù)Winchester and Floyd, 1977; c圖據(jù)Peccerillo and Taylor, 1976; d圖據(jù)Maniar and Piccoli, 1989。

    拉薩地塊中段數(shù)據(jù)來(lái)源同圖7, 球粒隕石和原始地幔標(biāo)準(zhǔn)化數(shù)據(jù)Sun and McDonough, 1989。

    表3 切瓊地區(qū)典中組流紋巖的鋯石Hf同位素分析結(jié)果

    注:Lu/Hf=(176Lu/177Hf)s/(176Lu/177Hf)CHUR?1,Hf()=10000×{[(176Hf/177Hf)s?(176Lu/177Hf)s×(e?1)]/[(176Hf/177Hf)CHUR, 0?(176Lu/177Hf)CHUR×(e?1)]?1},DM= 1/λ×ln{1+[(176Hf/177Hf)s?(176Hf/177Hf)DM]/[(176Lu/177Hf)s?(176Lu/177Hf)DM,DM2=DM(Hf)?(DM(Hf)?)×[(cc?s)/(cc?DM)], (176Lu/177Hf)CHUR=0.0332, (176Hf/177Hf)CHUR, 0=0.28272(Blichert-Toft et al., 1997), (176Lu/177Hf)DM=0.0332, (176Hf/177Hf)DM=0.282772(Griffin et al., 2000),cc=0.015,DM= ?0.548,=1.867×10?11a?1(S?derlund et al., 2004),為對(duì)應(yīng)鋯石的年齡。

    數(shù)據(jù)來(lái)源: 林周地區(qū)據(jù)李皓揚(yáng)等, 2007; 新嘎果地區(qū)據(jù)唐攀等, 2018; 南木林地區(qū)據(jù)付燕剛, 2017; 拉薩基底據(jù)Zhu et al., 2011。

    數(shù)據(jù)來(lái)源: (1) 拉薩地塊東段據(jù)李皓揚(yáng)等, 2007; He et al., 2007; 楊輝等, 2013; Huang et al., 2015; Zhu et al., 2015; 陳貝貝等, 2016; 唐攀等, 2018; Yan et al., 2019; 韋乃韶等, 2019; 周鵬等, 2019; (2) 拉薩地塊中段據(jù)梁銀平等, 2010; 謝冰晶等, 2013; 張小強(qiáng), 2013; 付燕剛, 2017; 丁帥等, 2017; 李勇等, 2018; 陳蘭樸等, 2019; (3) 拉薩地塊西段據(jù)王喬林, 2011; 姜軍勝, 2018; 曹延等, 2019。

    流紋巖巖漿源區(qū)方面, 該樣品具有較低Cr(8.00×10?6~17.0×10?6)和Ni(0.40×10?6~0.80×10?6)含量, 其更接近地殼來(lái)源(地殼Cr和Ni含量分別約為5.00×10?6~33.0×10?6和3.00×10?6~39.0×10?6; Lara et al., 2017)。大離子親石元素Ba和Sr出現(xiàn)明顯的負(fù)異常, 也指示巖漿源區(qū)中有地殼物質(zhì)參與(趙志丹等, 2001; 丁楓等, 2017; 唐攀等, 2018)。另一方面, 該火山巖的Th/U值為4.23~8.26(平均值為6.13); Nb/Ta值為11.57~13.14(平均值為12.69), 顯示殼幔混合的特點(diǎn)(下地殼Th/U值為3.80~6.00, 原始地幔Th/U值為4.00, 幔源巖漿Nb/Ta值為17.5±2, 殼源巖漿Nb/Ta值為11.0~12.0; Gao et al., 2004)。

    鑒于鋯石Hf同位素在地質(zhì)演化示蹤上的優(yōu)越性, 其被廣泛應(yīng)用到巖漿源區(qū)研究中(吳福元等, 2007)。在拉薩地塊東段, 李皓揚(yáng)等(2007)、唐攀等(2018)分別對(duì)林周盆地典中組最底部英安質(zhì)火山角礫巖和新嘎果地區(qū)典中組流紋質(zhì)巖屑晶屑凝灰?guī)r進(jìn)行鋯石Hf同位素研究, 得到了相似的Hf()值和2DM, 其Hf()均為正值(Hf()=+0.5~+9),2DM在556~1104 Ma之間, 表明拉薩地塊東段典中組火山巖主要來(lái)源于地幔組分, 但混染了古老地殼物質(zhì)。在拉薩地塊中段, 付燕剛(2017)對(duì)南木林地區(qū)典中組火山巖進(jìn)行鋯石Hf同位素原位分析, 兩套樣品Hf()在正負(fù)值范圍內(nèi)波動(dòng)(?4.7~+3.5與?5.2~+2), 均值小于0,2DM在933~1548 Ma之間, 表明了該火山巖為殼幔互相作用的產(chǎn)物。在拉薩地塊西段, Jiang et al. (2018)對(duì)諾倉(cāng)地區(qū)典中組玄武質(zhì)?安山質(zhì)晶屑凝灰?guī)r和同期的花崗斑巖分別進(jìn)行Sr-Nd同位素研究(Nd()= ?5.9~?4.3)與鋯石Hf同位素研究(Hf()=?22~?6, 均值為?13.3), 結(jié)合地球化學(xué)特征, 認(rèn)為該地區(qū)典中組中酸性火山巖主要來(lái)源于古拉薩地塊基底物質(zhì)的重熔, 并混入了少量幔源物質(zhì)。

    研究區(qū)典中組流紋巖鋯石Hf同位素研究結(jié)果顯示,Hf()值在?2.1~+2.3之間, 平均值為?0.6, 大部分小于0(圖9b、c), 對(duì)應(yīng)的二階段模式年齡2DM變化于875~1117 Ma之間, 平均為1037 Ma(圖9a), 表明源區(qū)有少量幔源巖漿加入, 但大部分為先存古老地殼組分的改造、再循環(huán)。對(duì)比前人研究成果, 進(jìn)一步證實(shí)了于楓等(2010)通過(guò)對(duì)比措麥地區(qū)與林周地區(qū)林子宗群火山巖Sr-Nd-Pb同位素組成, 得出拉薩地塊地殼基底組成不均一性的結(jié)論。拉薩地塊東段火山巖主要來(lái)源于地幔組分, 并在巖漿演化過(guò)程中混染了古老地殼(Mo et al., 2007, 2008; 韋乃韶等, 2019); 拉薩地塊西段火山巖可能源于幔源巖漿與殼源物質(zhì)間的相互作用(王喬林, 2011), 或者主要來(lái)源于古拉薩地塊基底物質(zhì)的重熔, 并混入了少量幔源物質(zhì)(Jiang et al., 2018); 拉薩地塊中段林子宗群中酸性火山巖為地幔源區(qū)基性分異巖漿和陸殼重熔的酸性巖漿按不同比例混合而成(梁銀平等, 2010), 或主要是古老地殼物質(zhì)再循環(huán), 并混有少量幔源巖漿??傮w看來(lái), 典中組火山巖巖漿源區(qū)從東向西幔源組分不斷減少, 而古老地殼物質(zhì)增多的變化趨勢(shì)愈加明顯, 預(yù)示拉薩地塊中段和西段存在古老結(jié)晶基底的可能。

    另一方面, 拉薩地塊中段林子宗群典中組火山巖顯示了較大范圍的Sr/Nd和Ba/Th變化值以及較小Th/Yb和Th/Nb變化值(圖11a、b), 這一特征與經(jīng)歷了相關(guān)流體富集的地幔楔產(chǎn)生的巖漿一致, 反映了典中組火山巖初始巖漿源區(qū)流體富集; 巖石相對(duì)虧損Ti、P等高場(chǎng)強(qiáng)元素以及明顯的陸緣弧火山巖特征, 表明源區(qū)經(jīng)歷過(guò)由俯沖板片脫水產(chǎn)生的流體交代作用(Rapp et al., 2003; 侯增謙等, 2006a); 從La/Sm-La圖解(圖12a)可以看出, 切瓊地區(qū)流紋巖與拉薩地塊中段典中組火山巖均表現(xiàn)為主要受部分熔融作用的控制, 暗示俯沖的新特提斯洋殼板片脫水產(chǎn)生的流體交代地幔楔使其發(fā)生部分熔融。Inger and Harris (1993)提出Rb/Sr值較高的熔體一般是熔融過(guò)程中流體較少的結(jié)果, 而Rb/Sr值升高是由于較低程度的部分熔融所致。切瓊地區(qū)流紋巖Rb/Sr值(20.59~51.81)相對(duì)較高, 且Rb/Sr值與Ba含量呈負(fù)相關(guān)關(guān)系, 表明該成巖過(guò)程中流體參與較少, 可能來(lái)源于白云母的脫水熔融(圖12b), 進(jìn)一步說(shuō)明該流紋巖巖漿源區(qū)也經(jīng)歷了新特提斯洋板片回卷過(guò)程中的古地殼的部分熔融。SiO2含量與(La/Sm)N值呈正相關(guān)關(guān)系, 與(Dy/Yb)N值呈負(fù)相關(guān)關(guān)系(圖12c、d), 反映巖漿演化過(guò)程中存在一定的角閃石分異結(jié)晶(Macpherson et al., 2006)。同時(shí), Ti、P呈負(fù)異常, 表明了金紅石、磷灰石的殘余, Sr和Ba的虧損以及明顯的負(fù)Eu異常, 說(shuō)明斜長(zhǎng)石和鉀長(zhǎng)石也經(jīng)歷了一定的結(jié)晶分異作用。

    綜上所述, 切瓊地區(qū)流紋巖具有陸緣弧火山巖特征, 巖漿源區(qū)主要是古老地殼物質(zhì)再循環(huán), 且混有少量幔源巖漿, 其巖漿演化過(guò)程主要受部分熔融控制, 存在一定的結(jié)晶分異, 總體顯示了具有新特提斯俯沖洋殼、大洋巖石圈地幔與大陸地殼的混合趨勢(shì)。

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

    研究區(qū)典中組火山巖是一套以爆發(fā)相、溢流相為主的中性、中酸性和酸性火山巖系, 主要由中性、中酸性熔巖、火山碎屑巖類組成, 是一套從中性火山巖向中酸性(酸性)火山巖, 由鈣堿性系列向高鉀鈣堿性系列連續(xù)演化的火山巖組合。地球化學(xué)特征表明, 研究區(qū)流紋巖屬高鉀鈣堿性系列、少量為鈣堿性系列, 反映其發(fā)育于擠壓造山帶(島弧、活動(dòng)大陸邊緣、陸?陸與陸內(nèi)造山帶)。在lg-lg圖解中, 研究區(qū)流紋巖與拉薩地塊中段典中組樣品全部落入造山帶(島弧及活動(dòng)大陸邊緣)環(huán)境區(qū)域(圖13a), 也進(jìn)一步表明了研究區(qū)流紋巖形成于造山帶環(huán)境。該流紋巖高Al2O3、低TiO2以及較低的Ni、Cr含量(Ni=0.40×10?6~0.80×10?6, 平均值為0.62×10?6; Cr=8.00×10?6~17.0×10?6, 平均值為12.7×10?6)等特征與島弧或活動(dòng)大陸邊緣弧環(huán)境相吻合; 另一方面, 樣品虧損Ti、P等高場(chǎng)強(qiáng)元素, 以及呈輕稀土元素富集、重稀土元素虧損的右傾稀土元素配分模式, 均指示其具有島弧或活動(dòng)大陸邊緣火山巖的構(gòu)造屬性, 也表明了該流紋巖是在特提斯洋北向俯沖于亞洲大陸之下的構(gòu)造背景下形成。在Rb/30-Hf-3Ta圖解中, 流紋巖樣品與拉薩地塊中段典中組火山巖主要落在火山弧區(qū)域, 并有向同碰撞?碰撞后區(qū)域演化的趨勢(shì)(圖13b), 暗示典中組火山巖的構(gòu)造環(huán)境可能經(jīng)歷了從俯沖到碰撞的變化過(guò)程。巖石成因方面, 該流紋巖主要為古老地殼物質(zhì)部分熔融產(chǎn)生的酸性巖漿的產(chǎn)物。

    拉薩地塊中段數(shù)據(jù)來(lái)源同圖7, 底圖據(jù)Jiang et al., 2018。

    拉薩地塊中段數(shù)據(jù)來(lái)源同圖7, 圖b據(jù)Inger and Harris, 1993。

    前人研究認(rèn)為, 在印度?亞洲大陸碰撞之前, 新特提斯洋板片以低角度或平向俯沖于亞洲大陸之下(Yin and Harrison, 2000; Ding et al., 2003; Chung et al., 2005; Wen et al., 2008; Lee et al., 2009)。70~40 Ma之間, 新特提斯洋板片以較陡的角度向北俯沖于亞洲大陸之下, 而后經(jīng)歷了新特提斯洋板片的回卷和斷裂, 引發(fā)了一系列巖漿爆發(fā)事件(Wen et al., 2008; Ji et al., 2009; Lee et al., 2012; Zhu et al., 2015)。其中, 新特提斯洋板片陡向俯沖于亞洲大陸之下, 降低了地幔楔初熔溫度, 隨后俯沖流體交代地幔楔使其熔融(Lee et al., 2009, 2012), 形成基性巖漿, 并在陸殼底部聚集。隨著俯沖的不斷加深, 新特提斯洋板片發(fā)生回卷, 擾動(dòng)了巖石圈地幔使其產(chǎn)生對(duì)流, 產(chǎn)生的熱量導(dǎo)致古老地殼發(fā)生部分熔融, 進(jìn)而產(chǎn)生混入少量幔源物質(zhì)的酸性巖漿。之后, 酸性巖漿與上侵的基性巖漿混合, 形成了典中組火山巖(Chung et al., 2005; 侯增謙等, 2006a; 張小強(qiáng), 2013)。這一過(guò)程常常用來(lái)推導(dǎo)隨后的俯沖末期和印度?亞洲大陸碰撞早期的新特提斯洋板片斷裂時(shí)限問(wèn)題, 即印度?亞洲大陸初始碰撞時(shí)限(Lee et al., 2009, 2012)。具有陸緣弧火山巖特征的林子宗群典中組火山巖應(yīng)是新特提斯洋巖石圈北向俯沖的結(jié)果, 其年齡應(yīng)該早于印度?亞洲大陸碰撞的初始年齡, 其底部年齡可能代表了俯沖作用年齡。現(xiàn)已報(bào)道的典中組底部年齡集中在70 Ma左右(李皓揚(yáng)等, 2007; 王喬林, 2011; 楊輝等, 2013; Zhu et al., 2015; 陳貝貝等, 2016; 唐攀等, 2018; 周鵬等, 2019; 陳蘭樸等, 2019), 表明新特提斯洋此時(shí)已經(jīng)陡沖; 而典中組頂部年齡集中在60 Ma(楊輝等, 2013; Zhu et al., 2015; 陳貝貝等, 2016; 陳蘭樸等, 2019), 表明印度?亞洲大陸初始碰撞發(fā)生在60 Ma之后。

    因此, 根據(jù)前人對(duì)于典中組火山巖的構(gòu)造演化的研究, 結(jié)合其地球化學(xué)特征及構(gòu)造環(huán)境判別圖解, 切瓊地區(qū)流紋巖最有可能形成于俯沖的島弧環(huán)境, 為新特提斯洋板片回卷引發(fā)的古地殼熔融的產(chǎn)物。

    4.4 成礦潛力分析

    在切瓊地區(qū), 由于印度?亞洲陸塊的碰撞, 出現(xiàn)大規(guī)模的巖漿活動(dòng), 使得以高鉀鈣堿性火山巖為主的典中組火山巖廣泛出露, 接觸變質(zhì)、熱液蝕變等現(xiàn)象較為普遍, 斷裂構(gòu)造十分發(fā)育, 為內(nèi)生礦產(chǎn)的富集提供了良好的成礦條件, 形成了研究區(qū)鉛、鋅、銀、鐵等多金屬為主的礦產(chǎn)。例如, 查嘎勒姐鉛鋅礦點(diǎn), 其含礦巖性為典中組灰黑色火山角礫巖, 金屬礦物呈膠結(jié)物充填于火山碎屑之間或節(jié)理裂隙中; 礦石礦物為方鉛礦、磁鐵礦, 脈石礦物為白云母、綠簾石、綠泥石等; 圍巖蝕變主要有綠簾石化、綠泥石化和黃鐵礦化等; 該礦點(diǎn)類型為火山熱液型。色拉鉛鋅礦(化)點(diǎn), 含礦巖性為典中組火山角礫巖中殘留的泥質(zhì)粉砂巖, 礦化分布于巖石顆??紫吨谢蚬?jié)理裂隙面上, 礦化蝕變類型有磁鐵礦化、黃鐵礦化、角巖化和鉛鋅礦化等; 礦石礦物主要為鉛、鋅、褐鐵礦等。土壤地球化學(xué)顯示, Pb、Zn、Cu等成礦元素峰值區(qū)與地表礦化體套合較好, 該礦化蝕變與巖漿熱液和火山熱液先后疊加作用有關(guān)。

    區(qū)域上, 林子宗群典中組火山巖年齡集中在70~60 Ma, 與同一構(gòu)造背景下的查個(gè)勒Pb-Zn礦床(高順寶等, 2012)、龍根Ag-Zn礦床(段志明等, 2014)等成巖成礦年齡一致。作為重要的賦礦層位, 典中組火山巖區(qū)內(nèi)發(fā)現(xiàn)如納如松多Pb-Zn礦床(紀(jì)現(xiàn)華等, 2014)、斯弄多Ag-Pb-Zn多金屬礦床(唐菊興等, 2016; 丁帥等, 2017)等淺成低溫?zé)嵋盒偷V床, 打加錯(cuò)、諾倉(cāng)等矽卡巖型Ag-Pb-Zn礦床, 北納Pb-Zn礦床等熱液脈型礦床(姜軍勝, 2018)?;鹕綆r源區(qū)方面, 岡底斯成礦帶中西段林子宗群典中組火山巖受控于拉薩地塊地殼基底, 是俯沖和再循環(huán)的大陸地殼物質(zhì)與幔源巖漿混合作用的產(chǎn)物, 有利于鉛鋅銀礦床的形成(黃瀚霄等, 2019)。

    拉薩地塊中段數(shù)據(jù)來(lái)源同圖7; a圖據(jù)Rittmann, 1971; b圖據(jù)Harris et al., 1986。

    綜上所述, 切瓊地區(qū)大面積出露的林子宗群典中組火山巖, 具有一定的尋找淺成低溫?zé)嵋盒?、矽卡巖型礦床的潛力。

    5 結(jié) 論

    (1) 切瓊地區(qū)流紋巖形成于晚白堊世, 鋯石U-Pb年齡為67.1±0.1 Ma, 屬林子宗群典中組火山巖。拉薩地塊林子宗群典中組火山巖噴發(fā)時(shí)間上具有東段和中段比西段持續(xù)時(shí)間長(zhǎng)的特點(diǎn)。

    (2) 切瓊地區(qū)流紋巖屬過(guò)鋁質(zhì)鈣堿性?高鉀鈣堿性系列, 富集大離子親石元素Rb、Th、U和K, 虧損高場(chǎng)強(qiáng)元素Ti、P和大離子親石元素Ba、Sr, 富集輕稀土元素, 虧損重稀土元素, 且具有明顯的負(fù)Eu異常, 總體顯示陸緣弧火山巖特征。鋯石Hf同位素特征(Hf()=?2.1~+2.3, 平均為?0.6)顯示該火山巖巖漿源區(qū)主要是古老地殼物質(zhì)再循環(huán), 并且混有少量幔源成分, 巖漿演化過(guò)程中存在一定的結(jié)晶分異。與拉薩地塊其他地區(qū)典中組火山巖的對(duì)比表明, 拉薩地塊典中組火山巖巖漿源區(qū)存在不均一性, 從東向西幔源組分減少而古老地殼物質(zhì)逐漸增多。

    (3) 切瓊地區(qū)典中組火山巖形成于俯沖的島弧環(huán)境, 具有大陸邊緣弧火山巖特征, 為新特提斯巖石圈板片回卷引發(fā)的古地殼熔融的主要產(chǎn)物, 具一定的尋找淺成低溫?zé)嵋盒?、矽卡巖型礦床潛力。

    致謝:中國(guó)地質(zhì)大學(xué)(武漢)資源學(xué)院閆茂強(qiáng)同學(xué)和高強(qiáng)同學(xué)在野外工作及論文撰寫過(guò)程中提供了幫助與指導(dǎo), 中國(guó)地質(zhì)科學(xué)院礦產(chǎn)資源研究所唐菊興研究員和一位匿名審稿人提出了建設(shè)性的修改意見, 在此一并致以誠(chéng)摯的感謝。

    曹延, 康志強(qiáng), 許繼峰, 李強(qiáng), 楊鋒, 韋天偉, 韋乃韶, 王睿. 2019. 拉薩地塊西部獅泉河地區(qū)典中組火山巖年代學(xué)、地球化學(xué)特征及其構(gòu)造意義. 地球科學(xué), 45(5): 1573–1592.

    陳貝貝, 丁林, 許強(qiáng), 岳雅慧, 謝靜. 2016. 西藏林周盆地林子宗群火山巖的精細(xì)年代框架. 第四紀(jì)研究, 36(5): 1037–1054.

    陳蘭樸, 黃澤森, 江巴多吉, 達(dá)瓦次仁, 塔爾杰. 2019. 西藏神公地區(qū)典中組火成巖鋯石U-Pb年齡及地球化學(xué)特征. 地質(zhì)通報(bào), 38(7): 1127–1135.

    丁楓, 徐忠彪, 劉壽航, 李躍, 高建國(guó), 范宇航, 李青. 2017. 西藏措勤赤馬地區(qū)古近系年波組火山巖LA-ICP-MS鋯石U-Pb年齡、地球化學(xué)特征及其成因. 地質(zhì)論評(píng), 63(4): 1102–1116.

    丁林, Satybaev M, 蔡福龍, 王厚起, 宋培平, 紀(jì)偉強(qiáng), 許強(qiáng), 張利云, Qasim M, Baral U. 2017. 印度與歐亞大陸初始碰撞時(shí)限、封閉方式和過(guò)程. 中國(guó)科學(xué): 地球科學(xué), 47(3): 293–309.

    丁帥, 陳毓川, 唐菊興, 謝富偉, 胡古月, 楊宗耀, 施碩, 李于海, 楊洪鈺. 2017. 林子宗群火山巖與成礦關(guān)系: 以斯弄多淺成低溫?zé)嵋盒偷V床為例. 礦床地質(zhì), 36(5): 1074–1092.

    董國(guó)臣, 莫宣學(xué), 趙志丹, 王亮, 周肅. 2005. 拉薩北部林周盆地林子宗火山巖層序新議. 地質(zhì)通報(bào), 24(6): 549–557.

    段志明, 李光明, 李應(yīng)栩, 王保弟, 張暉, 黃勇, 郭琳, 段瑤瑤. 2014. 中岡底斯成礦帶龍根鉛鋅礦床含礦斑巖年代學(xué)與地球化學(xué)特征. 礦床地質(zhì), 33(3): 625–638.

    付燕剛. 2017. 西藏拉薩地塊南部新特提斯洋俯沖成巖成礦作用. 北京: 中國(guó)地質(zhì)科學(xué)院博士學(xué)位論文: 64–67.

    付燕剛, 胡古月, 唐菊興, 劉秋平, 王豪. 2017. 西藏斯弄多低硫化型淺成低溫?zé)嵋篈g-Pb-Zn礦床: Si-H-O同位素的示蹤應(yīng)用. 地質(zhì)學(xué)報(bào), 91(4): 836–848.

    高順寶, 鄭有業(yè), 田立明, 張眾, 屈文俊, 劉敏院, 鄭海濤, 鄭磊, 朱繼華. 2012. 西藏查個(gè)勒銅鉛鋅礦成巖成礦時(shí)代及意義. 地球科學(xué), 37(3): 507–514.

    侯增謙, 莫宣學(xué), 高永豐, 楊志明, 董國(guó)臣, 丁林. 2006a. 印度大陸與亞洲大陸早期碰撞過(guò)程與動(dòng)力學(xué)模型——來(lái)自西藏岡底斯新生代火成巖證據(jù). 地質(zhì)學(xué)報(bào), 80(9): 1233–1248.

    侯增謙, 趙志丹, 高永豐, 楊志明, 江萬(wàn). 2006b. 印度大陸板片前緣撕裂與分段俯沖: 來(lái)自岡底斯新生代火山?巖漿作用證據(jù). 巖石學(xué)報(bào), 22(4): 761–774.

    胡新偉, 馬潤(rùn)則, 陶曉風(fēng), 劉登忠, 石和. 2007. 西藏措勤地區(qū)典中組火山巖地球化學(xué)特征及構(gòu)造背景. 成都理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 34(1): 15–22.

    黃瀚霄, 張林奎, 劉洪, 李光明, 黃勇, 蘭雙雙, 呂夢(mèng)鴻. 2019. 西藏岡底斯成礦帶西段礦床類型、成礦作用和找礦方向. 地球科學(xué), 44(6): 1876–1887.

    紀(jì)現(xiàn)華, 孟祥金, 楊竹森, 張乾, 田世洪, 李振清, 劉英超, 于玉帥. 2014. 西藏納如松多隱爆角礫巖型鉛鋅礦床絹云母Ar-Ar定年及其地質(zhì)意義. 地質(zhì)與勘探, 50(2): 281–290.

    姜軍勝. 2018. 岡底斯西段林子宗群火山巖區(qū)多金屬礦床成因及找礦潛力. 武漢: 中國(guó)地質(zhì)大學(xué)(武漢)博士學(xué)位論文: 48–65.

    李皓揚(yáng), 鐘孫霖, 王彥斌, 朱弟成, 楊進(jìn)輝, 宋彪, 劉敦一, 吳福元. 2007. 藏南林周盆地林子宗火山巖的時(shí)代、成因及其地質(zhì)意義: 鋯石U-Pb年齡和Hf同位素證據(jù). 巖石學(xué)報(bào), 23(2): 493–500.

    李勇, 張士貞, 李?yuàn)^其, 秦雅東, 鞏小棟. 2018. 拉薩地塊中段查孜地區(qū)典中組火山巖鋯石U-Pb年齡及地質(zhì)意義. 地球科學(xué), 43(8): 2755–2766.

    李再會(huì), 鄭來(lái)林, 李軍敏, 夏祥標(biāo). 2008. 岡底斯中段林子宗火山巖巖石地球化學(xué)特征. 礦物巖石地球化學(xué)通報(bào), 27(1): 20–27.

    梁銀平, 朱杰, 次邛, 何衛(wèi)紅, 張克信. 2010. 青藏高原岡底斯帶中部朱諾地區(qū)林子宗群火山巖鋯石U-Pb年齡和地球化學(xué)特征. 地球科學(xué), 35(2): 211–223.

    莫宣學(xué), 潘桂棠. 2006. 從特提斯到青藏高原形成: 構(gòu)造?巖漿事件的約束. 地學(xué)前緣, 13(6): 43–51.

    莫宣學(xué), 趙志丹, 鄧晉福, 董國(guó)臣, 周肅, 郭鐵鷹, 張雙全, 王亮亮. 2003. 印度?亞洲大陸主碰撞過(guò)程的火山作用響應(yīng). 地學(xué)前緣, 10(3): 135–148.

    莫宣學(xué), 趙志丹, 周肅, 董國(guó)臣, 廖忠禮. 2007. 印度?亞洲大陸碰撞的時(shí)限. 地質(zhì)通報(bào), 26(10): 1240–1244.

    莫宣學(xué), 趙志丹, 朱弟成, 喻學(xué)惠, 董國(guó)臣, 周肅. 2009. 西藏南部印度?亞洲碰撞帶巖石圈: 巖石學(xué)?地球化學(xué)約束. 地球科學(xué), 34(1): 17–27.

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

    唐菊興, 丁帥, 孟展, 胡古月, 高一鳴, 謝富偉, 李壯, 袁梅, 楊宗耀, 陳國(guó)榮, 李于海, 楊洪鈺, 付燕剛. 2016. 西藏林子宗群火山巖中首次發(fā)現(xiàn)低硫化型淺成低溫?zé)嵋盒偷V床——以斯弄多銀多金屬礦為例. 地球?qū)W報(bào), 37(4): 461–470.

    唐攀, 唐菊興, 鄭文寶, 冷秋鋒, 林彬, 唐曉倩, 王豪. 2018. 西藏新嘎果地區(qū)典中組火山巖年代學(xué)、Hf同位素及地球化學(xué)特征. 巖石礦物學(xué)雜志, 37(1): 47–60.

    王喬林. 2011. 岡底斯西段林子宗群火山巖的地球化學(xué)特征及鋯石年代學(xué)研究. 北京: 中國(guó)地質(zhì)大學(xué)(北京)碩士學(xué)位論文: 45–50.

    韋乃韶, 康志強(qiáng), 楊鋒, 冉孟蘭, 李強(qiáng), 韋天偉, 劉迪, 曹延, 陳歡, 李岱鮮. 2019. 西藏拉薩地塊東部南緣昌果地區(qū)典中組火山巖的年代學(xué)、地球化學(xué)及成因. 地球化學(xué), 48(1): 30–42.

    吳福元, 李獻(xiàn)華, 鄭永飛, 高山. 2007. Lu-Hf同位素體系及其巖石學(xué)應(yīng)用. 巖石學(xué)報(bào), 23(2): 185–220.

    吳元保, 鄭永飛. 2004. 鋯石成因礦物學(xué)研究及其對(duì)U-Pb年齡解釋的制約. 科學(xué)通報(bào), 49(16): 1589–1604.

    西藏自治區(qū)地質(zhì)調(diào)查院. 2016. 西藏切窮地區(qū)礦產(chǎn)遠(yuǎn)景調(diào)查報(bào)告.

    謝冰晶, 周肅, 謝國(guó)剛, 田明中, 廖忠禮. 2013. 西藏岡底斯中段孔隆至丁仁勒地區(qū)林子宗群火山巖鋯石SHRIMP年齡和地球化學(xué)特征的區(qū)域?qū)Ρ? 巖石學(xué)報(bào), 29(11): 3803–3814.

    楊輝, 向樹元, 王欣, 拉巴頓珠, 王勇. 2013. 西藏馬鄉(xiāng)地區(qū)典中組年齡厘定及其構(gòu)造背景. 地質(zhì)科技情報(bào), 32(4): 89–96.

    于楓, 李志國(guó), 趙志丹, 謝國(guó)剛, 董國(guó)臣, 周肅, 朱弟成, 莫宣學(xué). 2010. 西藏岡底斯帶中西部措麥地區(qū)林子宗火山巖地球化學(xué)特征及意義. 巖石學(xué)報(bào), 26(7): 2217– 2225.

    張小強(qiáng). 2013. 西藏則學(xué)地區(qū)林子宗群火山巖巖石成因及其與成礦關(guān)系的研究. 武漢: 中國(guó)地質(zhì)大學(xué)(武漢)碩士學(xué)位論文: 33–68.

    趙志丹, 莫宣學(xué), 張雙全, 郭鐵鷹, 周肅, 董國(guó)臣, 王勇. 2001. 西藏中部烏郁盆地碰撞后巖漿作用——特提斯洋殼俯沖再循環(huán)的證據(jù). 中國(guó)科學(xué)(D輯), 31(S1): 20–26.

    周鵬, 劉恭喜, 李強(qiáng), 余波. 2019. 岡底斯中段廳宮地區(qū)林子宗火山巖鋯石U-Pb年齡和地球化學(xué)特征. 礦物巖石地球化學(xué)通報(bào), 38(2): 352–360.

    周肅, 莫宣學(xué), 董國(guó)臣, 趙志丹, 邱瑞照, 王亮亮, 郭鐵鷹. 2004. 西藏林周盆地林子宗火山巖40Ar/39Ar年代格架. 科學(xué)通報(bào), 49(20): 2095–2103.

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

    Annen C, Blundy J D, Sparks R S J. 2006. The genesis of intermediate and silicic magmas in deep crustal hot zones., 47(3): 505–539.

    Blichert-Toft J, Chauvel C, Albarède F. 1997. Separation of Hf and Lu for high-precision isotope analysis of rock samples by magnetic sector-multiple collector ICP-MS., 127(3): 248–260.

    Cai F L, Ding L, Yue Y H. 2011. Provenance analysis of upper Cretaceous strata in the Tethys Himalaya, southern Tibet: Implications for timing of India-Asia collision., 305(1): 195–206.

    Chung S L, Chu M F, Zhang Y Q, Xie Y W, Lo C H, Lee T Y, Lan C Y, Li X H, Zhang Q, Wang Y Z. 2005. Tibetan tectonic evolution inferred from spatial and temporal variations in post-collisional magmatism., 68(3–4): 173–196.

    Colleps C L, Mckenzie N R, Horton B K, Webb A A G, Ng Y W, Singh B P. 2020. Sediment provenance of pre- and post-collisional Cretaceous-Paleogene strata from the frontal Himalaya of northwest India., 534: 1–15.

    Decelles P G, Gehrels G E, Najman Y, Martin A J, Carter A, Garzanti E. 2004. Detrital geochronology and geochemistryof Cretaceous-Early Miocene strata of Nepal: Implicationsfor timing and diachroneity of initial Himalayan orogenesis., 227(3–4): 313–330.

    Decelles P G, Kapp P, Gehrels G E, Ding L. 2014. Paleocene- Eocene foreland basin evolution in the Himalaya of southern Tibet and Nepal: Implications for the age of initial India-Asia collision., 33(5): 824–849.

    Ding L, Kapp P, Zhong D L, Deng W M. 2003. Cenozoic volcanism in Tibet: Evidence for a transition from oceanic to continental subduction., 44(10): 1833–1865.

    Donaldson D G, Webb A A G, Menold C A, Kylander-Clark A R C, Hacker B R. 2013. Petrochronology of Himalayan ultrahigh-pressure eclogite.,41(8): 835–838.

    Dupont-Nivet G, Lippert P C, Van Hinsbergen D J J, Meijers M J M, Kapp P. 2010. Palaeolatitude and age of the Indo-Asia collision: Palaeomagnetic constraints., 182(3): 1189–1198.

    Feeley T C, Dungan M A. 1996. Compositional and dynamic controls on mafic-silicic magma interactions at continental arc volcanoes: Evidence from Cordón El Guadal, Tatara-San Pedro Complex, Chile., 37(6): 1547–1577.

    Gao S, Rudnick R L, Yuan H L, Liu X M, Liu Y S, Xu W L, Ling W L, John A, Wang X C, Wang Q H. 2004. Recycling lower continental crust in the North China craton., 432(7019): 892–897.

    Griffin W L, Belousova E A, Shee S R, Pearson N J, O’Reilly S Y. 2004. Archean crustal evolution in the northern Yilgarn Craton: U-Pb and Hf-isotope evidence from detrital zircons., 131(3–4): 231–282.

    Harris N B W, Pearce J A, Tindle A G. 1986. Geochemicalcharacteristics of collision-zone magmatism.,,, 19(1): 67–81.

    He S D, Kapp P, Decelles P G, Gehrels G E, Heizler M. 2007. Cretaceous-Tertiary geology of the Gangdese Arc in the Linzhou area, southern Tibet., 433(1–4): 15–37.

    Hu X M, Sinclair H D, Wang J G, Jiang H H, Wu F Y. 2012a. Late Cretaceous-Palaeogene stratigraphic and basin evolution in the Zhepure Mountain of southern Tibet: Implications for the timing of India-Asia initial collision., 24(5): 520–543.

    Hu X M, Wang J G, Boudagher-Fadel M, Garzanti E, An W. 2016. New insights into the timing of the India-Asia collision from the Paleogene Quxia and Jialazi Formations of the Xigaze forearc basin, South Tibet., 32: 76–92.

    Hu Z C, Liu Y S, Gao S, Liu W G, Zhang W, Tong X R, Lin L, Zong K Q, Li M, Chen H H. 2012b. ImprovedHf isotope ratio analysis of zircon using newly designed X skimmer cone and jet sample cone in combination with the addition of nitrogen by laser ablationmultiple collector ICP-MS., 27(9): 1391–1399.

    Huang W T, Dupont-Nivet G, Lippert P C, van Hinsbergen D J J, Dekkers M J, Waldrip R, Ganer?d M, Li X C, Guo Z J, Kapp P. 2015. What was the Paleogene latitude of the Lhasa terrane? A reassessment of the geochronologyand paleomagnetism of Linzizong volcanic rocks (Linzhou basin, Tibet)., 34(3): 594–622.

    Inger S, Harris N. 1993. Geochemical constraints on leucogranite magmatism in the Langtang Valley, Nepal Himalaya., 34(2): 345–368.

    Ji W Q, Wu F Y, Chung S L, Li J X, Liu C Z. 2009. Zircon U-Pb geochronology and Hf isotopic constraints on petrogenesis of the Gangdese batholith, southern Tibet., 262(3–4): 229–245.

    Jiang J S, Zheng Y Y, Gao S B, Zhang Y C, Huang J, Liu J, Wu S, Xu J, Huang L L. 2018. The newly-discovered Late Cretaceous igneous rocks in the Nuocang district: Products of ancient crust melting trigged by Neo-Tethyan slab rollback in the western Gangdese., 308–309: 294–315.

    Lara P, Oyhant?abal P, Dadd K. 2017. Post-collisional, Late Neoproterozoic, high-Ba-Sr granitic magmatism from the Dom Feliciano Belt and its cratonic foreland, Uruguay: Petrography, geochemistry, geochronology, and tectonic implications., 277: 178–198.

    Lee H Y, Chung S L, Ji J Q, Qian Q, Gallet S, Lo C H, Lee T Y, Zhang Q. 2012. Geochemical and Sr-Nd isotopic constraints on the genesis of the Cenozoic Linzizong volcanic successions, southern Tibet., 53: 96–114.

    Lee H Y, Chung S L, Lo C H, Ji J Q, Lee T Y, Qian Q, Zhang Q. 2009. Eocene Neotethyan slab breakoff in southern Tibet inferred from the Linzizong volcanic record., 477(1–2): 20–35.

    Leech M, Singh S, Jain A, Klemperer S, Manickavasagam R. 2005. The onset of India-Asia continental collision: Early, steep subduction required by the timing of UHP metamorphism in the western Himalaya., 234(1–2): 83–97.

    Lippert P C, van Hinsbergen D J J, Dupont-Nivet G. 2014. Early Cretaceous to present latitude of the central proto-Tibetan Plateau: A paleomagnetic synthesis with implications for Cenozoic tectonics, paleogeography, and climate of Asia., (507): 1–21.

    Liu A L, Wang Q, Zhu D C, Zhao Z D, Liu S A, Wang R, Dai J G, Zheng Y C, Zhang L L. 2018. Origin of the. 50 Ma Linzizong shoshonitic volcanic rocks in the eastern Gangdese arc, southern Tibet., 304–307: 374–387.

    Liu Y S, Gao S, Hu Z C, Gao C G, Zong K Q, Wang D B. 2010. Continental and oceanic crust recycling-induced melt-peridotite interactions in the Trans-North China Orogen: U-Pb dating, Hf isotopes and trace elements in Zircons from mantle xenoliths., 51: 392–399.

    Liu Y S, Hu Z C, Gao S, Günther D, Xu J, Gao C G, Chen H H. 2008.analysis of major and trace elements of anhydrous minerals by LA-ICP-MS without applying an internal standard., 257(1–2): 34–43.

    Ludwig K R. 2003. ISOPLOT 3.0: A Geochronological Toolkit for Microsoft Excel. Berkeley Geochronology Center, Special Publication, 3: 1–70.

    Macpherson C G, Dreher S T, Thirlwall M F. 2006. Adakites without slab melting: High pressure differentiation of island arc magma, Mindanao, the Philippines., 243(3–4): 581–593.

    Maniar P D, Piccoli P M. 1989. Tectonic discrimination of granitoids., 101(5): 635–643.

    Mcculloch M T, Kyser T K, Woodhead J D, Kinsley L. 1994. Pb-Sr-Nd-O isotopic constraints on the origin of rhyolites from the Taupo Volcanic Zone of New Zealand: Evidence for assimilation followed by fractionation from basalt., 115(3): 303–312.

    Meng J, Wang C S, Zhao X X, Coe R, Li Y L, Finn D. 2012. India-Asia collision was at 24°N and 50 Ma: Palaeomagnetic proof from southernmost Asia., 2(1): 1–11.

    Middlemost E A K. 1994. Naming materials in the magma/ igneous rock system., 37(3): 215–224.

    Miller J A, Harris C. 2007. Petrogenesis of the Swaziland and northern natal rhyolites of the Lebombo Rifted volcanic margin, south east Africa., 48(1): 185–218.

    Mo X X, Hou Z Q, Niu Y Q, Dong G C, Qu X M, Zhao Z D, Yang Z M. 2007. Mantle contributions to crustal thickening during continental collision: Evidence from Cenozoic igneous rocks in southern Tibet., 96(1–2): 225–242.

    Mo X X, Niu Y L, Dong G C, Zhao Z D, Hou Z Q, Zhou S, Ke S. 2008. Contribution of syncollisional felsic magmatism to continental crust growth: A case study of the Paleogene Linzizong volcanic Succession in southern Tibet., 250(1–4): 49–67.

    Najman Y, Appel E, Boudagher-Fadel M, Bown P, Carter A, Garzanti E, Godin L, Han J, Liebke U, Oliver G, Parrish R, Vezzoli G. 2010. Timing of India-Asia collision: Geological, biostratigraphic, and palaeomagnetic constraints., 115, B12416.

    O’Brien P J, Zotov N, Law R, Khari M A, Jan M Q. 2001. Coesite in eclogite from the Upper Kaghan Valley, Pakistan: A first record and implications., 29: 435–438.

    Orme D A, Carrapa B, Kapp P. 2015. Sedimentology, provenance and geochronology of the upper Cretaceous- lower Eocene western Xigaze forearc basin, southern Tibet., 27(4): 387–411.

    Pankhurst R J, Riley T R, Fanning C M, Kelley S P. 2000. Episodic silicic volcanism in Patagonia and the Antarctic Peninsula: Chronology of magmatism associated with the break-up of Gondwana., 41(5): 605–625.

    Peccerillo A, Taylor S R. 1976. Geochemistry of Eocene calc-alkaline volcanic rocks from the Kastamonu area, Northern Turkey., 58(1): 63–81.

    Rapp R P, Shimizu N, Norman M D. 2003. Growth of early continental crust by partial melting of eclogite., 425(6958): 605–609.

    Riley T R, Leat P T, Pankhurst R J, Harris C. 2001. Origins of large volume rhyolitic volcanism in the Antarctic Peninsula and Patagonia by crustal melting., 42(6): 1043–1065.

    Rittmann A. 1971. The probable origin of high-alumina basalts., 34(2): 414–420.

    Roddaz M, Said A, Guillot S, Antoine P, Montel J, Martin F, Darrozes J. 2011. Provenance of Cenozoic sedimentary rocks from the Sulaiman fold and thrust belt, Pakistan: Implications for the palaeogeography of the Indus drainage system., 168(2): 499–516.

    S?derlund U, Patchett P J, Vervoort J D, Isachsen C E. 2004. The176Lu decay constant determined by Lu-Hf and U-Pb isotope systematics of Precambrian mafic intrusions., 219(3–4): 311–324.

    Sun S S, Mcdonough W F. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes.,,, 42(1): 313–345.

    Wang C S, Dai J G, Zhao X X, Li Y L, Graham S A, He D F, Ran B, Meng J. 2014. Outward-growth of the Tibetan Plateau during the Cenozoic: A review., 621: 1–43.

    Wang J G, Hu X M, Jansa L, Huang Z C. 2011. Provenanceof the Upper Cretaceous-Eocene Deep-Water sandstones in Sangdanlin, Southern Tibet: Constraints on the timingof initial India-Asia collision., (119): 293–309.

    Wen D R, Liu D Y, Chung S L, Chu M F, Ji J Q, Zhang Q, Song B, Lee T Y, Yeh M W, Lo C H. 2008. Zircon SHRIMP U-Pb ages of the Gangdese Batholith and implications for Neotethyan subduction in southern Tibet., 252(3–4): 191–201.

    Wilde S A, Valley J W, Peck W H, Graham C M. 2001. Evidence from detrital zircons for the existence of continental crust and oceans on the Earth 4.4 Gyr ago., 409(6817): 175–178.

    Winchester J A, Floyd P A. 1977. Geochemical discriminationof different magma series and their differentiation products using immobile elements., 20(4): 325–343.

    Yan H Y, Long X P, Li J, Wang Q, Zhao B S, Shu C T, Gou L L, Zuo R. 2019. Arc andesitic rocks derived from partial melts of mélange diapir in subduction zones: Evidence from whole-rock geochemistry and Sr-Nd-Mo isotopes of the Paleogene Linzizong Volcanic succession in Southern Tibet.:, 124(1): 456–475.

    Yi Z Y, Huang B C, Chen J S, Chen L W, Wang H L. 2011. Paleomagnetism of early Paleogene marine sediments in southern Tibet, China: Implications to onset of the India-Asia collision and size of Greater India., 309: 153–165.

    Yin A, Harrison T M. 2000. Geologic evolution of the Himalayan Tibetan Orogen., 28(1): 211–280.

    Zhu B, Kidd W S F, Rowley D B, Currie B S, Shafique N. 2005. Age of initiation of the India-Asia collision in the East-Central Himalaya., (113): 265–285.

    Zhu D C, Wang Q, Zhao Z D, Chung S L, Cawood P A, Niu Y L, Liu S A, Wu F Y, Mo X X. 2015. Magmatic record of India-Asia collision., 5(1): 1–8.

    Zhu D C, Zhao Z D, Niu Y L, Dilek Y, Hou Z Q, Mo X X. 2013. The origin and pre-Cenozoic evolution of the Tibetan Plateau., 23(4): 1429–1454.

    Zhu D C, Zhao Z D, Niu Y L, Mo X X, Chung S L, Hou Z Q, Wang L Q, Wu F Y. 2011. The Lhasa Terrane: Record of a microcontinent and its histories of drift and growth., 301(1–2): 241– 255.

    Zong K, Liu Y S, Gao C G, Hu Z C, Gao S, Gong H J. 2010.U-Pb dating and trace element analysis of zircons in thin sections of eclogite: Refining constraints on the ultra high-pressure metamorphism of the Sulu terrane, China, 269(3): 237–251.

    Zircon U-Pb Ages, Hf Isotopes and Geochemistry of the Rhyolite in the Dianzhong Formation from the Qieqiong Area, Xietongmen County, Tibet

    ZENG Cheng1, YAN Maoqiang1, SHEN Zhiyuan1, GAO Qiang1, WEI Junhao1*, MAO Guozheng2, DENG Yongming2

    (1. Faculty of Earth Resources, China University of Geosciences (Wuhan), Wuhan 430074, Hubei, China; 2. Tibet Engineering Exploration and Construction (Group) Co., Ltd., Lhasa 850000, Tibet, China)

    TheQieqiong area is located in the middle Lhasa block, Tibet, where volcanic rocks of the Linzizong Group are widespread. The Linzizong volcanic succession is regarded as the magmatism response to the continent collision between India and Asia, which is essentially important in discovering the initial collisional time and processes. This paper reports the whole-rock geochemical, Hf isotope and zircon U-Pb geochronological results of the rhyolite in this area. Tthe dating results show that the age of the rhyolite in the study area is 67.1±0.1 Ma and belongs to the volcanic rock of the Dianzhong Formation of the Linzizong Group. Chemical analyses show that the rhyolites in this area have high Si, high K, and low Ti, with peraluminous nature; they are characterized by LREE enrichment and obvious REE differentiation (Gd/Yb)N=0.57–1.16, with significant negative Eu anomalies (δEu=0.05–0.08), they are enriched in LILEs (Pb, Rb, Th and U), depleted in HFSEs (P and Ti), and obviously depleted in Ba and Sr, thus belong to calc-alkalic and peraluminous series with characteristics of continental margin arc volcanic rocks. Zircon grains from the rhyolite have176Hf/177Hf values in the range of 0.282673 to 0.282804,Hf() values in the range of ?2.1 to +2.3, and two-stage model ages (2DM) of 875 Ma to 1117 Ma. Overall, the rhyolite in the Dianzhong Formation from the Qieqiong area was formed in a subduction tectonic setting, derived mainly from recycled crustal material, and mixed with a small amount of mantle-derived magma. Considering the related deposits in the area and the discovered ore spots in the Qieqiong area, we suggest that the volcanic rocks of the Dianzhong Formation in the Qieqiong area have promising potential for epithermal and skarn deposits prospecting.

    Linzizong Group; Dianzhong Formation; rhyolite; Late Cretaceous; middle Lhasa block

    2020-05-22;

    2020-11-06

    國(guó)家自然科學(xué)基金項(xiàng)目(41772071)資助。

    曾成(1994–), 男, 碩士研究生, 礦產(chǎn)普查與勘探專業(yè)。E-mail: 13075979289@163.com

    魏俊浩(1961–), 男, 教授, 博士生導(dǎo)師, 從事礦產(chǎn)勘查研究工作。E-mail: junhaow@163.com

    P595; P597; P581

    A

    1001-1552(2022)01-0154-021

    国内精品久久久久精免费| 国产蜜桃级精品一区二区三区| 久久久水蜜桃国产精品网| 亚洲一区二区三区不卡视频| 99riav亚洲国产免费| 久久精品综合一区二区三区| 99久久无色码亚洲精品果冻| 99riav亚洲国产免费| 给我免费播放毛片高清在线观看| 久久久国产精品麻豆| 日日干狠狠操夜夜爽| 亚洲精品久久国产高清桃花| 美女cb高潮喷水在线观看 | 亚洲精品中文字幕一二三四区| 国产精品香港三级国产av潘金莲| 动漫黄色视频在线观看| 不卡一级毛片| 国产精品亚洲美女久久久| 国产亚洲精品久久久com| 精品一区二区三区视频在线 | 91av网站免费观看| 久久久久国内视频| 男女那种视频在线观看| 亚洲成人久久性| 深夜精品福利| 成人性生交大片免费视频hd| 99国产精品99久久久久| 两个人的视频大全免费| 久久久久久大精品| 国产亚洲欧美98| 波多野结衣巨乳人妻| av在线蜜桃| 在线免费观看的www视频| 18美女黄网站色大片免费观看| 亚洲av电影不卡..在线观看| av福利片在线观看| 久久婷婷人人爽人人干人人爱| 91在线精品国自产拍蜜月 | 人人妻人人看人人澡| 成年女人看的毛片在线观看| 19禁男女啪啪无遮挡网站| 这个男人来自地球电影免费观看| 亚洲男人的天堂狠狠| 亚洲最大成人中文| 色综合亚洲欧美另类图片| 精品免费久久久久久久清纯| 色av中文字幕| 欧美精品啪啪一区二区三区| 一级a爱片免费观看的视频| 亚洲熟妇熟女久久| 最近最新免费中文字幕在线| 久久热在线av| 色哟哟哟哟哟哟| 国产精品免费一区二区三区在线| 国产亚洲av嫩草精品影院| 精品午夜福利视频在线观看一区| 啦啦啦观看免费观看视频高清| 亚洲熟女毛片儿| 亚洲成人久久性| 亚洲18禁久久av| 久久久色成人| 啦啦啦观看免费观看视频高清| 人妻夜夜爽99麻豆av| 天堂动漫精品| 后天国语完整版免费观看| 天堂网av新在线| 熟女电影av网| 久久精品国产清高在天天线| 亚洲欧美日韩高清专用| 国产成人精品久久二区二区免费| 在线视频色国产色| 一边摸一边抽搐一进一小说| 怎么达到女性高潮| 国产精品久久久久久亚洲av鲁大| 国产成人精品无人区| 成熟少妇高潮喷水视频| 国产亚洲av高清不卡| 免费av毛片视频| 十八禁人妻一区二区| 黑人欧美特级aaaaaa片| 国产三级中文精品| 叶爱在线成人免费视频播放| 国产成人啪精品午夜网站| 欧美不卡视频在线免费观看| 97超视频在线观看视频| 久久九九热精品免费| 亚洲av成人精品一区久久| 精品国产超薄肉色丝袜足j| 波多野结衣高清作品| 成人性生交大片免费视频hd| 中文资源天堂在线| 精品国产三级普通话版| 久久午夜亚洲精品久久| 99久久99久久久精品蜜桃| 18禁观看日本| 丰满的人妻完整版| 日本五十路高清| 99久久成人亚洲精品观看| 精品国内亚洲2022精品成人| 亚洲欧美日韩东京热| 色老头精品视频在线观看| 日韩中文字幕欧美一区二区| 一区二区三区高清视频在线| 一区二区三区高清视频在线| 美女大奶头视频| 婷婷六月久久综合丁香| 精品乱码久久久久久99久播| 深夜精品福利| 午夜亚洲福利在线播放| 最近视频中文字幕2019在线8| 亚洲av日韩精品久久久久久密| 听说在线观看完整版免费高清| 91在线精品国自产拍蜜月 | 嫩草影院入口| 日日摸夜夜添夜夜添小说| 午夜免费成人在线视频| 一边摸一边抽搐一进一小说| 91老司机精品| 最好的美女福利视频网| 中国美女看黄片| 国内精品久久久久精免费| 日韩欧美国产一区二区入口| 老鸭窝网址在线观看| 国产又黄又爽又无遮挡在线| av天堂在线播放| 搡老妇女老女人老熟妇| 成人18禁在线播放| 午夜视频精品福利| 91久久精品国产一区二区成人 | 国产激情偷乱视频一区二区| 日韩欧美国产一区二区入口| 日本a在线网址| 成人18禁在线播放| 久久久久亚洲av毛片大全| 午夜亚洲福利在线播放| 国产精品永久免费网站| 国产精品影院久久| 性色avwww在线观看| 91老司机精品| 一级a爱片免费观看的视频| a级毛片在线看网站| 国产精品一区二区免费欧美| 久久人人精品亚洲av| cao死你这个sao货| 看片在线看免费视频| 天天躁狠狠躁夜夜躁狠狠躁| tocl精华| 亚洲国产精品合色在线| 99在线人妻在线中文字幕| 欧美色欧美亚洲另类二区| 白带黄色成豆腐渣| 亚洲无线在线观看| 99久久成人亚洲精品观看| 国产高潮美女av| 偷拍熟女少妇极品色| 亚洲国产中文字幕在线视频| 久久久久久人人人人人| 久久亚洲真实| 在线国产一区二区在线| 他把我摸到了高潮在线观看| 在线观看日韩欧美| 亚洲专区国产一区二区| 亚洲无线在线观看| 国产 一区 欧美 日韩| 欧美+亚洲+日韩+国产| 国产黄色小视频在线观看| 国产人伦9x9x在线观看| 日本五十路高清| 中亚洲国语对白在线视频| 99热精品在线国产| 国产伦精品一区二区三区视频9 | 级片在线观看| 久久人人精品亚洲av| 免费看十八禁软件| 91在线观看av| 五月玫瑰六月丁香| 国产激情久久老熟女| av在线蜜桃| 黄频高清免费视频| 天堂av国产一区二区熟女人妻| www日本黄色视频网| 免费观看的影片在线观看| 国产亚洲精品综合一区在线观看| 国产高潮美女av| 一区福利在线观看| 一级毛片高清免费大全| 亚洲欧美精品综合一区二区三区| 日韩欧美精品v在线| 丰满人妻一区二区三区视频av | 亚洲精品一卡2卡三卡4卡5卡| 成年女人看的毛片在线观看| 国产又色又爽无遮挡免费看| 19禁男女啪啪无遮挡网站| 香蕉丝袜av| 桃红色精品国产亚洲av| 国产成人aa在线观看| 亚洲天堂国产精品一区在线| 好看av亚洲va欧美ⅴa在| 两个人的视频大全免费| 国产欧美日韩一区二区三| 亚洲狠狠婷婷综合久久图片| 亚洲成人久久爱视频| 丰满人妻熟妇乱又伦精品不卡| 日韩精品中文字幕看吧| 最新中文字幕久久久久 | 99久久综合精品五月天人人| 一级黄色大片毛片| 我要搜黄色片| 在线播放国产精品三级| 成人特级黄色片久久久久久久| 中文在线观看免费www的网站| 亚洲国产中文字幕在线视频| 国产av麻豆久久久久久久| 我的老师免费观看完整版| 国产又色又爽无遮挡免费看| 精品久久久久久久毛片微露脸| 天天躁日日操中文字幕| 老司机午夜十八禁免费视频| 亚洲欧洲精品一区二区精品久久久| 久久精品亚洲精品国产色婷小说| av福利片在线观看| 欧美在线黄色| 亚洲av成人av| 精品久久久久久久末码| www国产在线视频色| 麻豆成人av在线观看| 麻豆国产97在线/欧美| 国产欧美日韩一区二区精品| 美女被艹到高潮喷水动态| 老司机在亚洲福利影院| 久久久国产成人免费| 又爽又黄无遮挡网站| 搞女人的毛片| www.自偷自拍.com| 他把我摸到了高潮在线观看| 一个人免费在线观看电影 | 国产精品乱码一区二三区的特点| 亚洲成人久久性| 日韩欧美三级三区| 成人无遮挡网站| 特大巨黑吊av在线直播| 一个人看视频在线观看www免费 | av国产免费在线观看| 国产精品av久久久久免费| 午夜精品久久久久久毛片777| 啦啦啦韩国在线观看视频| 国产精品久久久av美女十八| 熟女人妻精品中文字幕| 国产99白浆流出| а√天堂www在线а√下载| av片东京热男人的天堂| 亚洲欧美日韩东京热| 别揉我奶头~嗯~啊~动态视频| 国内久久婷婷六月综合欲色啪| 国产精品99久久久久久久久| 啦啦啦免费观看视频1| 亚洲五月婷婷丁香| 制服人妻中文乱码| 岛国在线观看网站| 国产精品,欧美在线| 真人一进一出gif抽搐免费| 一级毛片女人18水好多| 小说图片视频综合网站| 成年女人永久免费观看视频| 丝袜人妻中文字幕| x7x7x7水蜜桃| 91九色精品人成在线观看| 亚洲人成网站在线播放欧美日韩| 国产不卡一卡二| 欧美一区二区国产精品久久精品| 亚洲国产色片| 色尼玛亚洲综合影院| 欧美激情久久久久久爽电影| 精品人妻1区二区| 琪琪午夜伦伦电影理论片6080| 国模一区二区三区四区视频 | 女同久久另类99精品国产91| 久久中文看片网| 成人性生交大片免费视频hd| 婷婷亚洲欧美| 国产精品精品国产色婷婷| 久久久久国产一级毛片高清牌| 亚洲精品中文字幕一二三四区| 国产精品久久久久久人妻精品电影| 老司机午夜福利在线观看视频| 亚洲无线观看免费| 国产探花在线观看一区二区| 亚洲中文字幕一区二区三区有码在线看 | 在线十欧美十亚洲十日本专区| 国产精品久久久久久亚洲av鲁大| 五月玫瑰六月丁香| 国产美女午夜福利| 超碰成人久久| 亚洲激情在线av| 在线观看免费午夜福利视频| 亚洲国产日韩欧美精品在线观看 | 国内少妇人妻偷人精品xxx网站 | 国产1区2区3区精品| 亚洲av成人不卡在线观看播放网| 观看美女的网站| 久久久久性生活片| 久久久久国产一级毛片高清牌| bbb黄色大片| 欧美最黄视频在线播放免费| av天堂在线播放| 免费看美女性在线毛片视频| 美女 人体艺术 gogo| 色综合婷婷激情| 日韩成人在线观看一区二区三区| 国产亚洲欧美98| 久久久久久九九精品二区国产| 国产成人影院久久av| 免费无遮挡裸体视频| 男女之事视频高清在线观看| 日本熟妇午夜| 51午夜福利影视在线观看| netflix在线观看网站| 日韩三级视频一区二区三区| 怎么达到女性高潮| 99在线人妻在线中文字幕| xxx96com| 国产高清视频在线观看网站| 午夜精品一区二区三区免费看| 国产精品一及| 天堂影院成人在线观看| 丰满的人妻完整版| 一本精品99久久精品77| 丰满人妻一区二区三区视频av | 在线观看美女被高潮喷水网站 | 在线播放国产精品三级| 一个人看视频在线观看www免费 | 一a级毛片在线观看| 最近最新中文字幕大全免费视频| 99久久综合精品五月天人人| av天堂中文字幕网| 哪里可以看免费的av片| 欧美另类亚洲清纯唯美| 精品国产三级普通话版| 高清在线国产一区| 精品国产三级普通话版| 黄色丝袜av网址大全| 国产主播在线观看一区二区| 日本在线视频免费播放| 久久天堂一区二区三区四区| 久久国产精品人妻蜜桃| 亚洲av免费在线观看| 久久午夜综合久久蜜桃| 日韩欧美三级三区| 国内精品久久久久精免费| 日本精品一区二区三区蜜桃| 人妻丰满熟妇av一区二区三区| 欧美午夜高清在线| 午夜激情欧美在线| 国产av麻豆久久久久久久| 又爽又黄无遮挡网站| 亚洲专区国产一区二区| 亚洲熟妇中文字幕五十中出| www国产在线视频色| 三级男女做爰猛烈吃奶摸视频| 亚洲熟妇熟女久久| 亚洲av电影不卡..在线观看| 成人特级黄色片久久久久久久| 国产三级黄色录像| 久久草成人影院| 成年版毛片免费区| 精品久久久久久久毛片微露脸| 国产v大片淫在线免费观看| 亚洲电影在线观看av| 一个人免费在线观看的高清视频| 成人18禁在线播放| 国产男靠女视频免费网站| 欧美成狂野欧美在线观看| 午夜日韩欧美国产| 日韩 欧美 亚洲 中文字幕| 成人午夜高清在线视频| 国产综合懂色| 2021天堂中文幕一二区在线观| 最近视频中文字幕2019在线8| 久久久国产成人免费| 国产成人av激情在线播放| 2021天堂中文幕一二区在线观| 夜夜爽天天搞| xxx96com| avwww免费| 欧美色欧美亚洲另类二区| 色老头精品视频在线观看| 亚洲性夜色夜夜综合| 久久午夜综合久久蜜桃| 亚洲国产精品sss在线观看| 午夜福利在线观看吧| 91九色精品人成在线观看| 国产精品一区二区三区四区久久| 99精品久久久久人妻精品| 变态另类成人亚洲欧美熟女| 99在线人妻在线中文字幕| 美女扒开内裤让男人捅视频| 黄色日韩在线| 亚洲,欧美精品.| 日本一二三区视频观看| 99国产精品99久久久久| 午夜免费成人在线视频| 听说在线观看完整版免费高清| bbb黄色大片| 欧美激情久久久久久爽电影| 亚洲美女黄片视频| 久久性视频一级片| 日本一本二区三区精品| 亚洲av中文av极速乱| av视频在线观看入口| 少妇的逼水好多| 国产精品久久久久久久电影| 自拍偷自拍亚洲精品老妇| 美女cb高潮喷水在线观看| 不卡视频在线观看欧美| 国产精品三级大全| 99热这里只有精品一区| 波多野结衣高清无吗| 精品久久久久久久末码| 国产又色又爽无遮挡免| 在线免费观看的www视频| 国产精品久久久久久精品电影小说 | 男女那种视频在线观看| 亚洲美女视频黄频| 久久久久久久亚洲中文字幕| 久久久久久大精品| 久久鲁丝午夜福利片| av卡一久久| 你懂的网址亚洲精品在线观看 | 一区二区三区高清视频在线| 一本久久精品| 两性午夜刺激爽爽歪歪视频在线观看| 欧美一区二区国产精品久久精品| 在线播放国产精品三级| 国产一区二区三区av在线| 国产麻豆成人av免费视频| 久久精品91蜜桃| 国产免费视频播放在线视频 | 久久久久久久久中文| 精品午夜福利在线看| 国产一级毛片在线| 91久久精品电影网| 精品人妻视频免费看| 国产成年人精品一区二区| 色综合站精品国产| 国产单亲对白刺激| 日韩人妻高清精品专区| 国产亚洲av嫩草精品影院| 黄色一级大片看看| 亚洲精品国产成人久久av| 天堂网av新在线| 99久久精品热视频| 中文字幕制服av| 久久国产乱子免费精品| 韩国av在线不卡| 欧美一区二区精品小视频在线| 大香蕉97超碰在线| 国产老妇女一区| 欧美成人a在线观看| 免费看光身美女| 久久久精品欧美日韩精品| 亚洲丝袜综合中文字幕| 亚洲一级一片aⅴ在线观看| 特大巨黑吊av在线直播| 久久欧美精品欧美久久欧美| 三级国产精品欧美在线观看| 69av精品久久久久久| 久久亚洲国产成人精品v| 成人特级av手机在线观看| 熟女人妻精品中文字幕| 男女那种视频在线观看| 草草在线视频免费看| 欧美3d第一页| 成人毛片60女人毛片免费| 国产精华一区二区三区| 小说图片视频综合网站| 久久精品久久精品一区二区三区| 日日摸夜夜添夜夜添av毛片| 黄色日韩在线| 日本午夜av视频| 久久久国产成人免费| 欧美日韩一区二区视频在线观看视频在线 | av免费在线看不卡| 亚洲国产精品久久男人天堂| 婷婷色综合大香蕉| 婷婷色av中文字幕| 亚洲成人精品中文字幕电影| 免费在线观看成人毛片| 国产精品综合久久久久久久免费| 免费观看性生交大片5| av在线播放精品| av专区在线播放| 亚洲一级一片aⅴ在线观看| 亚洲av电影不卡..在线观看| 神马国产精品三级电影在线观看| 亚洲欧美一区二区三区国产| 午夜福利高清视频| 国产淫片久久久久久久久| 欧美区成人在线视频| 美女内射精品一级片tv| 国产白丝娇喘喷水9色精品| 亚洲av电影在线观看一区二区三区 | 男女视频在线观看网站免费| 纵有疾风起免费观看全集完整版 | 国产精品伦人一区二区| 亚洲国产精品专区欧美| 看黄色毛片网站| 国产精品久久久久久av不卡| 天天躁夜夜躁狠狠久久av| 老司机福利观看| 在线免费观看的www视频| 亚洲国产日韩欧美精品在线观看| 欧美成人一区二区免费高清观看| 日韩一本色道免费dvd| 久久久久久久国产电影| 高清av免费在线| 少妇熟女aⅴ在线视频| 亚洲av免费高清在线观看| 91aial.com中文字幕在线观看| 国产亚洲av片在线观看秒播厂 | 日日干狠狠操夜夜爽| 亚洲中文字幕一区二区三区有码在线看| 麻豆av噜噜一区二区三区| 婷婷六月久久综合丁香| 看黄色毛片网站| 美女内射精品一级片tv| 日本免费一区二区三区高清不卡| 成人亚洲精品av一区二区| 国产av不卡久久| 看免费成人av毛片| 美女内射精品一级片tv| 午夜日本视频在线| 久久精品国产鲁丝片午夜精品| 男女国产视频网站| 国产高清不卡午夜福利| 久久久久性生活片| 丰满人妻一区二区三区视频av| 久久久久免费精品人妻一区二区| 成人无遮挡网站| 日本五十路高清| 国产综合懂色| 观看美女的网站| 看十八女毛片水多多多| 国模一区二区三区四区视频| 人人妻人人澡欧美一区二区| 日本一本二区三区精品| 少妇裸体淫交视频免费看高清| 人体艺术视频欧美日本| 国产精品国产三级国产专区5o | 久久久久久久久中文| 欧美成人一区二区免费高清观看| 国产av在哪里看| 国产大屁股一区二区在线视频| 国产精品女同一区二区软件| 国产免费福利视频在线观看| 午夜福利成人在线免费观看| 久久久久久久午夜电影| 人人妻人人澡人人爽人人夜夜 | 亚洲人成网站在线播| 久久久久久国产a免费观看| 波野结衣二区三区在线| 国产视频内射| 国产欧美另类精品又又久久亚洲欧美| 免费看光身美女| 成人性生交大片免费视频hd| 亚洲怡红院男人天堂| 小蜜桃在线观看免费完整版高清| 最后的刺客免费高清国语| 国产午夜精品久久久久久一区二区三区| 一级爰片在线观看| 69av精品久久久久久| 又粗又爽又猛毛片免费看| 国产精品永久免费网站| 最近视频中文字幕2019在线8| 亚洲国产高清在线一区二区三| 免费看日本二区| 天堂√8在线中文| h日本视频在线播放| 人妻系列 视频| 五月玫瑰六月丁香| 白带黄色成豆腐渣| 久久精品国产亚洲网站| 亚洲18禁久久av| 亚洲国产精品合色在线| 中文字幕精品亚洲无线码一区| 丰满人妻一区二区三区视频av| 国产成人freesex在线| 在线免费十八禁| 国产成人91sexporn| 搞女人的毛片| 国产成人aa在线观看| 国产精品一及| 久久久久国产网址| 国产黄色视频一区二区在线观看 | 日韩在线高清观看一区二区三区| 国产亚洲av嫩草精品影院| 在线观看一区二区三区| 精品一区二区三区人妻视频| 日韩欧美三级三区| 久久精品国产99精品国产亚洲性色| 91午夜精品亚洲一区二区三区| 国产午夜精品久久久久久一区二区三区| 国产中年淑女户外野战色| 免费观看a级毛片全部| 天堂影院成人在线观看| 国产亚洲5aaaaa淫片| 国产精品女同一区二区软件| 久久精品国产亚洲av涩爱| 久久99热这里只频精品6学生 | 欧美日韩综合久久久久久| 成人性生交大片免费视频hd| 亚洲欧美精品专区久久| 国产在线一区二区三区精 | 一二三四中文在线观看免费高清| 夫妻性生交免费视频一级片| 亚洲欧美中文字幕日韩二区|