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

    西藏南木林晚白堊世輝綠巖與花崗質(zhì)脈巖成因及其揭示的伸展背景*

    2015-07-21 08:51:38葉麗娟趙志丹劉棟朱弟成董國臣莫宣學(xué)胡兆初劉勇勝
    巖石學(xué)報 2015年5期
    關(guān)鍵詞:輝綠巖白堊鋯石

    葉麗娟 趙志丹 劉棟 朱弟成 董國臣 莫宣學(xué) 胡兆初 劉勇勝

    1.中國地質(zhì)大學(xué)地質(zhì)過程與礦產(chǎn)資源國家重點實驗室,地球科學(xué)與資源學(xué)院,北京 100083

    2.天津地質(zhì)礦產(chǎn)研究所,天津 300170

    3.中國地質(zhì)大學(xué)地質(zhì)過程與礦產(chǎn)資源國家重點實驗室,地球科學(xué)學(xué)院,武漢 430074

    1 引言

    青藏高原雅魯藏布縫合帶(IYZSZ)是一條板塊俯沖消減帶,它標(biāo)志著印度與亞洲大陸碰撞之前曾經(jīng)存在新特提斯洋(常承法,1984),特提斯洋的形成演化及其有關(guān)的構(gòu)造-巖漿作用過程近年來一直受到國內(nèi)外地學(xué)研究者的關(guān)注,被認(rèn)為是更好揭示印度與亞洲大陸碰撞以來青藏高原構(gòu)造演化的重要基礎(chǔ)(Zhu et al.,2013)。早至195Ma 以來,特提斯洋就開始向北俯沖削減到拉薩地塊南部,形成了島弧巖漿作用(Kang et al.,2014),并持續(xù)北向俯沖直到特提斯洋閉合、印度與亞洲大陸的對接和碰撞(Xu et al.,1985;Coulon et al.,1986;Harris et al.,1988;Yin and Harrison,2000;Mo et al.,2007,2008;Zhu et al.,2009;Zhao et al.,2009)。

    新一輪地質(zhì)大調(diào)查獲得的1∶25 萬區(qū)域地質(zhì)資料表明,中生代晚白堊世巖漿作用出露比較廣泛,呈大面積分布于南部拉薩地塊(朱弟成等,2008),目前已經(jīng)發(fā)現(xiàn)的侵入巖包括輝長巖、閃長巖、英云閃長巖、花崗閃長巖和花崗巖等各種巖石類型,構(gòu)成了平行于雅魯藏布縫合帶的大陸邊緣弧巖漿作用帶(Ji et al.,2009)。關(guān)于印度-亞洲大陸碰撞前的巖漿作用及深部動力學(xué)過程目前還存在頗多爭議,存在多種觀點和成因模型,例如新特提斯洋殼低角度北向俯沖(Wen et al.,2008a,b)、高角度北向俯沖(王莉等,2013)和新特提斯洋脊俯沖(管琪等,2010;Zhang et al.,2010a;Zhu et al.,2011)等。

    作為拉薩地塊南部晚白堊世巖漿作用研究的一部分,本文在南木林縣南部發(fā)現(xiàn)了一套輝綠巖體以及呈巖脈狀侵入的花崗質(zhì)脈巖,開展了這套基性巖和酸性脈巖的鋯石U-Pb年代學(xué)和地球化學(xué)研究,結(jié)合拉薩地塊南緣新近報道的同時期巖漿巖成果,探討了拉薩地塊南部晚白堊世巖漿巖的成因及其構(gòu)造意義。

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

    以金沙江縫合帶(JSSZ)、班公湖-怒江縫合帶(BNSZ)和雅魯藏布江縫合帶(YZSZ)為界將青藏高原從北向南分為松潘-甘孜復(fù)理石雜巖帶、羌塘、岡底斯和喜馬拉雅帶(Yin and Harrison,2000)。拉薩地塊位于班公湖-怒江縫合帶(BNSZ)和雅魯藏布縫合帶(YZSZ)之間,以獅泉河-拉果錯-永珠-納木錯-嘉黎蛇綠混雜巖帶(SLYNJOMZ)和沙莫勒-麥拉-洛巴堆-米拉山斷裂帶(SMLMF)為界,由北向南分為北部拉薩地塊、中部拉薩地塊和南部拉薩地塊(圖1a)(Zhu et al.,2013)。拉薩地塊南緣,沿著雅魯藏布縫合帶北側(cè)近東西向展布一條狹長花崗巖帶為岡底斯巖基(Chu et al.,2006;Ji et al.,2009;Wen et al.,2008a,b),其西端起始于岡仁波齊峰地區(qū),以喀喇昆侖斷裂為界與科希斯坦-拉達(dá)克巖基斷開,東側(cè)延伸到林芝地區(qū),是一條主要由I 型巖漿巖組成的侵入雜巖帶(Debon et al.,1986;Harris et al.,1988),與其并行展布的還有中生代、新生代火山巖和沉積地層,如中侏羅統(tǒng)-下侏羅統(tǒng)葉巴組(董彥輝等,2006;朱弟成等,2008;王偉等,2013;Kang et al.,2014)、上侏羅統(tǒng)-下白堊統(tǒng)桑日群(姚鵬等,2006;Zhu et al.,2009)、古近紀(jì)林子宗群(劉鴻飛,1993;莫宣學(xué)等,2003;周肅等,2004;李皓揚等,2007;Mo et al.,2008;Lee et al.,2009)。目前已獲得的南部拉薩地塊晚白堊世數(shù)據(jù)大部分均采自南木林至林芝之間(圖1b)。

    圖1 研究區(qū)地質(zhì)圖(a)青藏高原及拉薩地塊構(gòu)造單元劃分簡圖(據(jù)Zhu et al.,2011 修改),JSSZ-金沙江縫合帶;BNSZ-班公湖-怒江縫合帶;SNMZ-獅泉河-納木錯蛇綠混雜巖帶;LMF-洛巴堆-米拉山斷裂帶;YZSZ-雅魯藏布縫合帶;(b)南部拉薩地塊南木林至林芝地區(qū)地質(zhì)簡圖(據(jù)朱弟成等,2008 修改),文獻數(shù)據(jù)引自黃玉等,2010;管琪等,2010,2011;Wen et al.,2008a,b;Zhang et al.,2010;Zhu et al.,2011;(c)研究區(qū)輝綠巖和花崗質(zhì)脈巖的野外露頭;(d)南木林南部地區(qū)地質(zhì)簡圖Fig.1 Geological map of the study area(a)simplified geological map of Lhasa terrane and main units in the Tibetan Plateau (modified after Zhu et al.,2011),JSSZ:Jinsha Suture Zone;BNSZ:Bangong-Nujiang Suture Zone;SNMZ:Shiquan River-Nam Tso Ophiolitic Melange Zone;LMF:Luobadui-Mila Mountain Fault;YZSZ:Yarlung Zangbo Suture Zone;(b)geological map of the Namling to Linzhi area in southern Lhasa terrane (modified after Zhu et al.,2008),literature data:Huang et al.,2010;Guan et al.,2010,2011;Wen et al.,2008a,b;Zhang et al.,2010;Zhu et al.,2011;(c)field outcrop of the diabase and intruded granitic dikes in Namling region;(d)geological map of the Namling region

    本文研究區(qū)位于拉薩地塊南部南木林縣以南約25km 的棉將雜巖體內(nèi)(圖1c,d),該巖體出露面積約15km2,呈巖瘤狀產(chǎn)出,巖體主體巖性為輝綠巖,巖體內(nèi)部發(fā)育了多條花崗巖脈,巖脈多沿節(jié)理和裂隙貫入,脈寬1~2m,長幾十米,沿北東10°~30°延伸。野外共采集了4 件樣品,包括輝綠巖1件、花崗巖脈3 件。采樣點GPS 坐標(biāo)為N29°27.235’,E89°05.939’。樣品NM1104 為輝綠巖,巖石具典型的輝綠結(jié)構(gòu)(圖2a),主要礦物為斜長石(>50%)、輝石(20%)、角閃石(10%)、黑云母(10%)等。斜長石呈自形板狀晶體,之間成近三角形,聚片雙晶發(fā)育,輝石及角閃石等發(fā)育兩組解理,暗色礦物部分發(fā)生蝕變。樣品NM1101、NM1102、NM1103 均為花崗巖脈,巖石呈淺肉紅色,主要礦物為石英(40%)、鉀長石(表面高嶺土化)(30%)、斜長石(15%)及少量黑云母(5%)和白云母(5%)(圖2b)。

    3 分析方法

    鋯石單礦物分選在河北省區(qū)域地質(zhì)礦產(chǎn)調(diào)查研究所實驗室完成。鋯石陰極發(fā)光(CL)顯微照相在中國科學(xué)院地質(zhì)與地球物理研究所電子探針室完成。鋯石微量元素含量和U-Pb 同位素定年在中國地質(zhì)大學(xué)(武漢)地質(zhì)過程與礦產(chǎn)資源國家重點實驗室(GPMR)利用LA-ICP-MS 同時分析完成。激光剝蝕系統(tǒng)為GeoLas 2005,ICP-MS 為Agilent 7500a。對分析數(shù)據(jù)的離線處理(包括對樣品和空白信號的選擇、儀器靈敏度漂移校正、元素含量及U-Th-Pb 同位素比值和年齡計算)采用軟件ICPMSDataCal(Liu et al.,2008,2010a)完成。詳細(xì)的儀器操作條件和數(shù)據(jù)處理方法同Liu et al.(2008,2010a,b)。

    圖2 樣品顯微結(jié)構(gòu)照片(a)NM1104(輝綠巖),呈輝綠結(jié)構(gòu),正交偏光;(b)NM1101(正長花崗巖),花崗結(jié)構(gòu).視域左半部分為正交偏光,右半部分為單偏光Fig.2 Photo of micrographs(a)sample NM1104,diabase,diabasic texture,under cross-polarized light;(b)sample NM1101,granites,granular texture.Left half under crosspolarized light,right half under polarized light

    鋯石Hf 同位素測試是在中國地質(zhì)科學(xué)院礦產(chǎn)資源研究所國土資源部成礦作用與資源評價重點實驗室Neptune 多接收等離子質(zhì)譜和Newwave UP213 紫外激光剝蝕系統(tǒng)(LAMC-ICP-MS)上進行的,實驗過程中采用He 作為剝蝕物質(zhì)載氣,剝蝕直徑采用55μm,測定時使用鋯石國際標(biāo)樣GJ-1 和Plesovice 作為參考物質(zhì),分析點與U-Pb 定年分析點為同一位置。相關(guān)儀器運行條件及詳細(xì)分析流程見(侯可軍等,2007)。分析過程中鋯石標(biāo)準(zhǔn)GJ-1 的176Hf/177Hf 測試加權(quán)平均值分別為0.282007 ±0.000007(2σ,n =36),與文獻報道值(侯可軍等,2007;Morel et al.,2008)在誤差范圍內(nèi)完全一致。

    以上定年和鋯石微量元素、Hf 同位素測試結(jié)果分別列于表1、表2 和表3 中。

    主量和微量元素在中國地質(zhì)大學(xué)(武漢)地質(zhì)過程與礦產(chǎn)資源國家重點實驗室測定。主量元素采用X 射線熒光方法,實驗儀器為XRF-1800。微量元素采用Agilent 7500a 等離子體質(zhì)譜儀(ICP-MS)測定,分析精度優(yōu)于5%~10%。詳細(xì)測試方法和分析流程見Gao et al.(2002)。主、微量元素測試結(jié)果列于表4 中。

    4 分析結(jié)果

    4.1 鋯石U-Pb 年齡

    本文分別對1 件輝綠巖(NM1104)和1 件花崗巖脈(NM1103)進行了鋯石U-Pb 定年。鋯石陰極發(fā)光(CL)圖像和年齡諧和圖見圖3。輝綠巖(NM1104)樣品中的鋯石具有寬板狀韻律環(huán)帶,顆粒粒度較大(粒徑一般>100μm),鋯石的Th/U 比值均大于0.1(0.49~0.97),為典型的巖漿成因鋯石(Hoskin and Schaltegger,2003),剔除18 個測點中5 個不諧和點,其余13 個測點的206Pb/238U 加權(quán)平均年齡為93.8±0.8Ma(MSWD=0.17),為輝綠巖的形成時代。

    花崗巖脈(NM1103)樣品中鋯石大多數(shù)為自形,粒徑80~100μm,多具典型巖漿震蕩環(huán)帶,鋯石的Th/U 比值為0.35~0.68,剔除18 個測點中6 個不諧和點,剩余12 個點的206Pb/238U 加權(quán)平均年齡為92.4 ±1.0Ma(MSWD=0.078)。

    從鋯石U-Pb 年齡看,正長花崗巖脈(92.4 ±1.0Ma)以略晚的年齡,穿插到輝長-輝綠巖(93.8 ±0.8Ma)中。在測試誤差范圍內(nèi),兩者年齡相似,它們應(yīng)屬于同期巖漿活動的產(chǎn)物。

    4.2 鋯石微量元素

    花崗巖脈NM1103 鋯石Th 含量變化于156 ×10-6~503×10-6,Pb 含量幾乎都小于39 ×10-6,并具有高的Th/Pb 比值和變化范圍大的Nb 含量(圖4b)。LuN值變化于1289~2474,處于Hoskin and Ireland (2000)報道的鋯石的LuN值是球粒隕石的1000 到20000 倍范圍內(nèi),為典型的巖漿鋯石;鋯石的Y 含量為667 × 10-6~1464 × 10-6,與Poller et al.(2001)報道的花崗巖類鋯石的Y 豐度500 ×10-6~4534 ×10-6相當(dāng)。這些I 型花崗巖類鋯石以中等程度的負(fù)Eu 異常(Eu/Eu*= 0.36~0.60)、較高的(Nb/Pb)N比值(0.46~1.27)為特征,其富集HREE,具有明顯的正Ce 異常(圖4a),總體與巖漿鋯石的微量元素特征類似(Hoskin and Schaltegger,2003)。

    輝綠巖樣品NM1104 的稀土元素含量整體偏高,ΣREE=800 ×10-6~2546 ×10-6,并且與全巖稀土元素含量整體偏高相一致,而不符合鋯石稀土元素豐度隨巖性由超基性巖到酸性巖逐漸升高的趨勢,因此輝綠巖和花崗巖脈之間在成分上不是同一來源的巖漿的分異演化關(guān)系,而是來自于不同的源區(qū)。在球粒隕石標(biāo)準(zhǔn)化稀土元素圖解上(圖4a),樣品顯示為LREE 虧損、HREE 富集的特點,且具明顯的Ce 正異常和較強烈的Eu 負(fù)異常(Eu/Eu*=0.16~0.26)。

    表1 西藏南木林巖漿巖鋯石U-Pb 定年數(shù)據(jù)Table 1 Zircons U-Pb data of the magmatic rocks from Namling region

    Wang et al.(2012)的研究表明,雖然I 型和S 型花崗巖類鋯石均具有類似的火成巖類鋯石稀土元素配分型式(如正Ce 異常、負(fù)Eu 異常、富集HREE 等)(Hoskin and Schaltegger,2003),但在一些元素豐度(如Pb、Th 元素)和比值(如Th/Pb、(Nb/Pb)N和Eu/Eu*)等方面,卻存在明顯差別(圖4b),即I 型花崗巖類以低Pb 含量、高Th/Pb 和(Nb/Pb)N比值為特征,S 型花崗巖類以高Pb 含量、顯著負(fù)Eu 異常、低Th/Pb和(Nb/Pb)N比值為特征,而A 型花崗巖類這些參數(shù)的變化范圍介于I 型和S 型花崗巖類之間。

    源區(qū)物質(zhì)成分的差異,決定了其部分熔融產(chǎn)物(即I 型或S 型花崗質(zhì)熔體)的獨特地球化學(xué)特征。作為從花崗質(zhì)熔體中結(jié)晶出來的鋯石,必然繼承了源區(qū)物質(zhì)成分特點,從而表現(xiàn)出S 型花崗巖鋯石的Pb 含量比值高于I 型花崗巖,而Nb 含量和(Nb/Pb)N比值卻低于I 型花崗巖。本文的花崗巖脈與I 型相同或者相似,相比Th 低Pb 高,但仍不同于A 型和S 型,全巖數(shù)據(jù)也表示出了相似的特征,可能是鋯石中隨時間由Th 和U 衰變而積累了Pb 含量。因此鋯石和全巖的微量元素都表明本文的花崗巖脈成分屬于I 型花崗巖。

    圖3 南木林輝綠巖和花崗巖的鋯石U-Pb 年齡諧和圖(a、c)和陰極發(fā)光(CL)圖像(b、d)Fig.3 Zircon U-Pb concordia diagrams (a,c)and CL images (b,d)of representative zircons from diabase and granite from Namling region

    圖4 鋯石稀土和微量元素圖(a)稀土元素球粒隕石標(biāo)準(zhǔn)化圖(球粒隕石據(jù)Boynton,1984);(b)微量元素鋯石類型區(qū)分圖解,I 型、S 型花崗巖類的鋯石微量元素含量和比值代表性散點圖,內(nèi)置插圖為全巖地球化學(xué)數(shù)據(jù)投圖,I 和S 型花崗巖的成分區(qū)域根據(jù)Wang et al.(2012)Fig.4 Zircon REE and trace elements diagrams(a)chondrite-normalized REE patterns for zircons (chondrite values after Boynton,1984);(b)selected plots of trace element abundances and ratios in zircons from I-,S-type granitoids.The whole-rock geochemical data are also plotted (see insets)for comparison (Wang et al.,2012)

    表2 南木林巖漿巖鋯石微量元素數(shù)據(jù)(×10 -6)Table 2 Zircons trace element data of the magmatic rocks from Namling region (×10 -6)

    表3 南木林巖漿巖鋯石Hf 同位素數(shù)據(jù)Table 3 Zircon Hf isotopic data of the magmatic rocks from Namling region

    4.3 鋯石Hf 同位素

    對南木林花崗巖脈和輝綠巖的25 顆鋯石的原位Hf 同位素組成進行了分析(表3)表明,正長花崗巖脈樣品中鋯石的176Yb/177Hf 和176Lu/177Hf 比值范圍分別為0.046091~0.129083 和0.000751~0.003441。εHf(t)范圍為+ 7.7~+13.0,Hf 同位素地殼模式年齡tDMC=301~637Ma。輝綠巖中鋯石的176Yb/177Hf 和176Lu/177Hf 比值范圍分別為0.021627~0.098605 和0.000325~0.001670。εHf(t)范圍為+8.5~+14.2,Hf 同位素地殼模式年齡tDMC=219~592Ma。這表明,輝綠巖和其中的花崗巖脈都是來自于強烈虧損的源區(qū),顯示了新生幔源物質(zhì)的性質(zhì)。

    4.4 全巖地球化學(xué)

    花崗巖脈巖的主量元素特征(表4)表現(xiàn)為SiO2含量比較均勻(74.76%~76.02%),按照主量元素多參數(shù)分類結(jié)果(圖5a),花崗巖可詳細(xì)命名為正長花崗巖。巖石具有較高的K2O 含量(4.75%~6.18%)和全堿(K2O + Na2O)含量(7.57%~8.44%),低的TiO2(0.08%~0.11%)和P2O5(0.01%~0.02%)含量;巖石的鋁飽和指數(shù)A/CNK 為1.10~1.18,為具有過鋁質(zhì)特征的正長花崗巖,屬于高鉀鈣堿性-鉀玄質(zhì)系列的巖石(圖5b)。而輝綠巖樣品為典型的基性輝長質(zhì)成分(圖5a),具有較高的Fe2O3T(9.37%)和低MgO 含量(4.32%),全堿(K2O+Na2O)含量(5.03%)也較高,屬于高鉀鈣堿性系列巖石(圖5b)。

    圖5 南木林輝綠巖和花崗巖的巖石類型和系列劃分圖(a)巖石類型R1-R2圖解(據(jù)De La Roche et al.,1980);(b)巖石系列K2O-SiO2劃分圖解(據(jù)Rollinson,1993).文獻數(shù)據(jù)引自黃玉等,2010;管琪等,2010,2011;Wen et al.,2008a,b;Zhang et al.,2010;Zhu et al.,2011Fig.5 Discrimination diagrams of rock types and series from the Namling rocks(a)R1-R2 diagram(after De La Roche et al.,1980);(b)K2O-SiO2 diagram(after Rollinson,1993).Literature data:Huang et al.,2010;Guan et al.,2010,2011;Wen et al.,2008a,b;Zhang et al.,2010;Zhu et al.,2011

    圖6 南木林南部巖漿巖球粒隕石標(biāo)準(zhǔn)化稀土元素配分曲線圖(a,球粒隕石據(jù)Boynton,1984)和原始地幔標(biāo)準(zhǔn)化微量元素配分曲線圖(b,原始地幔數(shù)據(jù)Sun and McDonough,1989)拉薩地塊南部晚白堊世巖類數(shù)據(jù)自黃玉等,2010;管琪等,2010,2011;Wen et al.,2008;Zhang et al.,2010;Zhu et al.,2011Fig.6 Chondrite-normalized REE patterns (a,after Boynton,1984)and primitive mantle-normalized trace element patterns (b,after Sun and McDonough,1989)Late Cretaceous magmatic rocks from the southern Lhasa Terrane are from literatures (Huang et al.,2010;Guan et al.,2010,2011;Wen et al.,2008a,b;Zhang et al.,2010;Zhu et al.,2011)

    從微量元素特征看,輝綠巖顯示相對平坦的配分形式(圖6a),輕重稀土元素分異不明顯,具弱的Eu 異常(δEu =0.87),與拉薩地塊南緣報道的含紫蘇輝石的閃長巖和花崗閃長巖均類似(管琪等,2011;Zhu et al.,2011;Zhang et al.,2010;Ma et al.,2013a,b,c)。輝綠巖富集K、Rb、Sr、Ba 等大離子親石元素,虧損Nb、Ta、Ti 等高場強元素,顯示了楔形地幔源區(qū)部分熔融的島弧巖漿的特征,應(yīng)形成于俯沖帶環(huán)境;Pb 的含量高說明基性巖漿受到了地殼物質(zhì)的影響(圖6b)。正長花崗巖脈輕重稀土分餾較大,(La/Yb)N=11~64,具有微弱的Eu 異常(δEu=0.79~0.90),Nb、P、Ti 虧損,較中基性巖類具有輕稀土含量相似,重稀土含量較低的特征,而該地區(qū)報道的近于同時代的花崗閃長巖類有四個樣品具有Eu 正異常(管琪等,2010;黃玉等,2010;Zhu et al.,2011;Wen et al.,2008a,b),可能是斜長石堆晶作用所致。在微量元素方面的另一重要特征是,花崗巖具有明顯的高Sr/Y(41~84)、低Y (2.4 ×10-6~4.6 ×10-6)、低Yb (0.2×10-6~0.4 ×10-6)的特征,結(jié)合(La/Yb)N等特征看,巖石成分處于典型埃達(dá)克巖成分區(qū)域的邊緣,具有類似埃達(dá)克巖的特征。

    此外,輝綠巖全巖稀土元素含量明顯高于3 個花崗巖樣品,與前述的鋯石微量元素類似,不符合同一巖漿體系稀土元素豐度越到酸性巖越高的演化趨勢,表明輝綠巖和花崗巖脈之間不存在分異演化關(guān)系,即它們應(yīng)來自于不同源區(qū)。

    5 討論

    5.1 拉薩地塊南部晚白堊世巖漿作用時空分布

    根據(jù)區(qū)域地質(zhì)填圖資料,朱弟成等(2008)發(fā)現(xiàn)拉薩地塊南部存在晚白堊世巖漿活動,大致平行于雅魯藏布江縫合帶呈帶狀分布。近年來在拉薩地塊南部的日喀則到米林一線,越來越多的研究(Wen et al.,2008a,b;Ji et al.,2009,2014;Zhang et al.,2010;黃玉等,2010;管琪等,2010,2011;Zhu et al.,2011;Ma et al.,2013a,b,c)進一步揭示了晚白堊世研究活動的大規(guī)模分布,其年齡主要集中于78~100Ma,被稱為晚白堊世的“巖漿大爆發(fā)”(Ma et al.,2013a)。本文獲得的正長花崗巖脈和輝綠巖形成年齡分別為92.4Ma 和93.8Ma,也屬于晚白堊世早期巖漿作用。這些研究表明,在岡底斯帶南緣,確實明顯存在一條晚白堊世的巖漿巖帶,其巖石成分從基性的輝長-輝綠巖、到角閃石巖、閃長巖,以及大量的花崗巖類。

    晚白堊世的巖石類型也存在多樣性,出現(xiàn)了紫蘇花崗巖(Zhang et al.,2010;Guan et al.,2012;Ma et al.,2013c)、蘇長巖和角閃輝長巖(Ma et al.,2013a)、以及采自朗縣到臥龍鎮(zhèn)(管琪等,2010;Guan et al.,2012;Zheng et al.,2014)、里龍與米林一帶(Zhang et al.,2010a;Ma et al.,2013c)報道了87~90Ma 的埃達(dá)克質(zhì)巖漿活動埃達(dá)克質(zhì)巖石,反映了新特提斯洋俯沖消減過程中既有俯沖洋殼、同時也有洋脊俯沖等多種構(gòu)造環(huán)境的巖漿作用記錄。

    從巖石出露的空間特征看,78~100Ma 的巖石主要出露在拉薩地塊南緣,疊加在岡底斯巖基的南緣(Ji et al.,2009,2014;Zhang et al.,2010;黃玉等,2010;管琪等,2010,2011;Zhu et al.,2011;Guan et al.,2012;Ma et al.,2013a,b,c),尤其是在臥龍到米林一帶分布很多(Zhang et al.,2010;管琪等,2010,2011;Zhu et al.,2011;Guan et al.,2012;Ma et al.,2013a,b,c;Zheng et al.,2014)。在中部拉薩地塊(也稱岡底斯弧背斷隆帶)的門巴地區(qū)也出現(xiàn)了88Ma 的埃達(dá)克質(zhì)巖石(孟繁一等,2010;Meng et al.,2014)。本文位于南木林的研究區(qū)與上述研究的區(qū)別是,巖體出露疊加在岡底斯巖基的北緣,靠近中部拉薩地塊,處于岡底斯島弧的北緣或者弧后地區(qū)。

    5.2 拉薩地塊南部晚白堊世巖石成因與構(gòu)造環(huán)境

    本文研究的南木林輝綠巖和花崗巖脈,結(jié)合區(qū)域同時代巖漿巖研究結(jié)果,具有如下特征:

    圖7 鋯石εHf(t)-年齡圖解(a)和地殼模式年齡直方圖(b)文獻數(shù)據(jù)引自黃玉等,2010;管琪等,2010,2011;Zhu et al.,2011Fig.7 Zircon εHf(t)-age (Ma)diagram (a)and histogram of tCM values (b)Literature data from Huang et al.,2010;Guan et al.,2010,2011;Zhu et al.,2011

    圖8 拉薩地塊晚白堊世巖石構(gòu)造判別圖(a)花崗巖類Rb-Y+Nb 判別圖(Pearce et al.,1984);(b)鎂鐵質(zhì)巖類Th/Hf-Ta/Hf 圖解(汪云亮等,2001).文獻數(shù)據(jù)引自黃玉等,2010;管琪等,2010,2011;Wen et al.,2008a,b;Zhang et al.,2010;Zhu et al.,2011Fig.8 Geochemical discrimination plots of the magmatic rocks in southern Lhasa terrane(a)Rb-Y+Nb for granites (after Pearce et al.,1984 );(b)Th/Hf-Ta/Hf diagram for mafic rocks (after Wang et al.,2001).Literature data:Huang et al.,2010;Guan et al.,2010,2011;Wen et al.,2008a,b;Zhang et al.,2010;Zhu et al.,2011

    (1)巖石來自于虧損的源區(qū),且具有源區(qū)組成的不均一性。南木林輝綠巖(εHf(t)= +8.5~+14.2)和正長花崗巖(εHf(t)= +7.7~+13.0)具有虧損程度較大的鋯石Hf 同位素組成(圖7a),該特征與岡底斯帶南緣的晚白堊世巖漿巖完全一致(Ji et al.,2014),同時輝綠巖微量元素顯示了明顯的輕稀土和LILE 的富集、以及HFSE 的虧損(圖6),表明輝綠巖應(yīng)來自于遭受了俯沖交代的地幔楔部分熔融,以及來自于深部的強烈虧損的軟流圈地幔物質(zhì)的加入,因此該巖漿作用可能代表了巖石圈-軟流圈的相互作用過程,在拉薩地塊南緣的米林地區(qū)的93Ma 的蘇長巖和角閃輝長巖也代表了巖石圈-軟流圈的相互作用的產(chǎn)物(Ma et al.,2013a)。而花崗巖具有虧損的Hf 同位素特征,源區(qū)物質(zhì)的鋯石Hf 模式年齡峰值為300~400Ma(圖7b),表明其源于島弧基底的地幔注入到下地殼的新生地殼的部分熔融。同時,上述鋯石Hf 同位素的變化范圍均大于5 個ε 單位,顯示了源區(qū)組成的不均一性。

    (2)花崗巖脈屬于高分異的I 型花崗巖。南木林花崗巖脈的鋯石和全巖的微量元素(圖4)特征表明巖石屬于I 型花崗巖。同時,巖石又具有富SiO2(74.76%~76.02%)、富K2O(4.75%~6.18%)、富堿(K2O + Na2O = 7.57%~8.44%),以及過鋁質(zhì)特征(A/CNK=1.10~1.18),表明巖石是高分異的I 型花崗巖,在巖石中出現(xiàn)了白云母,可以作為過鋁質(zhì)的巖石學(xué)證據(jù)。那么,略晚期侵入到基性巖體中的花崗巖脈,是否可以代表基性巖漿分異作用的產(chǎn)物呢?前文述及輝綠巖的鋯石和全巖稀土元素含量都分別明顯整體高于花崗巖脈的稀土含量(圖4、圖6a),與同一巖漿體系稀土元素豐度越到酸性巖越高的演化趨勢矛盾,表明花崗巖脈不是基性巖漿分異演化的產(chǎn)物,它們具有不同的巖漿源區(qū)。

    (3)巖石可能代表了弧后伸展作用的產(chǎn)物。從研究區(qū)的巖石組合看,出現(xiàn)了基性巖石與后期的花崗巖脈共存的特征,顯示了類似“雙峰式”巖漿作用的侵入巖組合;脈體的產(chǎn)出也意味著該時期具有伸展的構(gòu)造背景。

    已有的研究結(jié)果表明,采用比較謹(jǐn)慎的方法,還是可以從花崗巖類的巖石類型和地球化學(xué)特征中獲得一些構(gòu)造環(huán)境方面的重要信息,如:島弧巖漿組合中英云閃長巖和花崗閃長巖占優(yōu)勢(也包括部分花崗巖或輝長巖),后期形成鉀玄巖系列巖石(圖5b),地球化學(xué)性質(zhì)主要屬鈣堿性偏鋁質(zhì)到弱過鋁質(zhì)的I 型花崗巖(也包括部分S 型花崗巖)(朱弟成等,2008)。南部拉薩地塊晚白堊世已有的所有數(shù)據(jù)特征完全符合以上島弧巖漿組合的特征。將南部拉薩地塊花崗巖類投在Pearce et al.(1984)提出的花崗巖構(gòu)造背景判別圖上,樣品落入火山弧花崗巖的范圍內(nèi)(圖8a),按照花崗巖構(gòu)造判別圖及微量元素特征,包括南木林花崗巖脈的南部拉薩地塊花崗巖類應(yīng)屬于島弧花崗巖,其總體構(gòu)造環(huán)境應(yīng)是新特提斯洋板塊一直向著拉薩地塊持續(xù)的俯沖和擠壓。

    同時拉薩地塊南緣出現(xiàn)的基性和酸性脈巖出現(xiàn),則表明在擠壓同時局部存在伸展作用。南木林地區(qū)晚白堊世輝綠巖具有較高的Zr 含量(84.93 ×10-6),在Th/Hf-Ta/Hf 圖解(圖8b)中落入大陸伸展帶/初始裂谷玄武巖區(qū)域,顯示出板內(nèi)鈣堿性玄武巖親緣性。

    結(jié)合前述的研究區(qū)處于岡底斯巖基帶的北緣、產(chǎn)出基性和酸性的雙峰式巖石組合、出現(xiàn)脈巖侵入等特征,本文認(rèn)為南木林地區(qū),甚至整個拉薩地塊南部,晚白堊世時期在總體為特提斯洋向北俯沖的大構(gòu)造背景控制下,也局部存在弧后伸展的構(gòu)造背景,向北在中部拉薩地塊南緣的門巴地區(qū)近于同時也出現(xiàn)了伸展性質(zhì)的巖漿作用(孟繁一等,2010;Meng et al.,2014)。

    5.3 區(qū)域構(gòu)造-巖漿演化意義

    新特提斯洋巖石圈的北向俯沖消減于拉薩地塊南緣,是解釋拉薩地塊南緣晚白堊世巖漿作用成因的主要模型。具體到不同研究地區(qū),則存在幾種不同觀點:(1)正常安第斯型島弧巖漿作用(Ji et al.,2009;Chu et al.,2006;李曉雄等,2015);(2)雅魯藏布新特提斯洋殼低角度或平板俯沖(Coulon et al.,1986;Wen et al.,2008a);(3)東部朗縣-米林地區(qū)的新特提斯洋脊俯沖(70~95Ma,Zhang et al.,2010;管琪等,2010;Guan et al.,2012);(4)俯沖的特提斯洋板片的回轉(zhuǎn)(Ma et al.,2014a)。

    從本文的南木林地區(qū)產(chǎn)出的基性巖和酸性脈巖的產(chǎn)出結(jié)果以及上述討論看,該巖漿作用顯示了靠近弧后的伸展構(gòu)造背景、軟流圈和巖石圈的共同參與、上地幔和下地殼同時產(chǎn)生巖漿作用等過程,結(jié)合目前認(rèn)為的78~100Ma 巖漿作用大量出現(xiàn)(Ma et al.,2013a),本文認(rèn)為除了特殊巖石類型需要應(yīng)用洋脊俯沖模式來解釋之外,板片斷離、板片回轉(zhuǎn)等都可以解釋上述晚白堊世的大規(guī)模巖漿作用過程。如果應(yīng)用板片回轉(zhuǎn)模式(Ma et al.,2013a),在先期總體俯沖作用的背景下,拉薩地塊可能在晚白堊時期已經(jīng)存在地殼縮短、抬升、變形等事件(Murphy et al.,1997;Yin and Harrison,2000;Kapp et al.,2007;Leier et al.,2007);大約在100~90Ma 期間,之前低角度北向俯沖的新特提斯洋板片由于重力原因,發(fā)生板塊回轉(zhuǎn)(roll-back)(Chung et al.,2005;Ma et al.,2013a),伴隨著上涌對流軟流圈提供熱量與物質(zhì),使得受俯沖交代過的上覆地幔楔發(fā)生部分熔融,產(chǎn)生了基性的巖漿;基性巖漿上侵到地殼后,誘發(fā)了島弧基底物質(zhì)的部分熔融,形成了中酸性的具有新生地殼性質(zhì)的巖漿。上述巖漿疊加侵位到岡底斯巖基中。本文的南木林巖漿巖位于巖基北緣,可能代表了弧后伸展的構(gòu)造背景下的產(chǎn)物,而俯沖帶深部的板片斷離和板片回轉(zhuǎn)過程,都可以導(dǎo)致地殼淺部形成伸展的構(gòu)造環(huán)境。

    6 結(jié)論

    (1)獲得了南部拉薩地塊南木林南部輝綠巖年齡為93.8 ± 0.8Ma,侵入其中的正長花崗巖脈年齡為92.4±1.0Ma。

    (2)輝綠巖和花崗巖的組合可能代表了弧后伸展背景的巖漿作用產(chǎn)物,輝綠巖可能是巖石圈和軟流圈共同作用的產(chǎn)物;花崗巖是新生的島弧基底的部分熔融,為高分異的I 型花崗巖。

    (3)南部拉薩地塊晚白堊世巖漿巖可能是在新特提斯洋殼北向俯沖過程中,發(fā)生板片斷離或者板片回轉(zhuǎn)過程中,交代的地幔楔與軟流圈作用熔融形成基性巖漿,基性巖漿底侵下地殼誘發(fā)下地殼部分熔融形成酸性巖漿,形成了后期侵入的脈巖。

    致謝 實驗中得到秦虹、陳海紅和劉碩的幫助;曾令森和戴緊根審閱論文并提出了很好的修改意見;在此一并表示感謝。

    Blichert-Toft J and Albarède F.1997.The Lu-Hf isotope geochemistry of chondrites and the evolution of the mantle-crust system.Earth and Planetary Science Letters,148(1-2):243-258

    Boynton WV.1984.Geochemistry of the rare earth elements:Meteorite studies.In: Henderson P (ed.).Rare Earth Element Geochemistry.Amsterdam:Elsevier,63-114

    Chang CF.1984.Geological structural features and evolution of the Yarlung Zangbo Suture Zone.In:Sino-French Investigation Report on Himalayan.Beijing:Science Press,327-339 (in Chinese)

    Chu MF,Chung SL,Song B,Liu DY,O’Reilly SY,Pearson N,Ji J and Wen DJ.2006.Zircon U-Pb and Hf isotope constraints on the Mesozoic tectonics and crustal evolution of southern Tibet.Geology,34(9):745-748

    Chung SL,Chu MF,Zhang YQ,Xie YW,Lo CH,Lee TY,Lan CY,Li XH,Zheng Q and Wang YZ.2005.Tibetan tectonic evolution inferred from spatial and temporal variations in post-collisional magmatism.Earth-Science Reviews,68(3-4):173-196

    Coulon C,Maluski H,Bollinger C and Wand 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

    De La Roche H,Leterrier J,Grandclaude P and Marchal M.1980.Classification of volcanic and plutonic rocks using R1-R2 diagrams and major element analysis-its relationships and current nomenclature.Chemical Geology,29(1-4):183-210

    Debon F,Le Fort P,Sheppard SM,Simon MF and Sonet J.1986.The four plutonic belts of the Transhimalaya-Himalaya:A chemical,mineralogical,isotopic,and chronological synthesis along a Tibet-Nepal section.Journal of Petrology,27(1):219-250

    Dong YH,Xu JF,Zeng QG,Wang Q,Mao GZ and Li J.2006.Is there a Neo-Tethys’subduction record earlier than arc volcanic rocks in the Sangri Group?Acta Petrologica Sinica,22(3):661-668 (in Chinese with English abstract)

    Gao S,Liu XM,Yuan HL,Hattendorf B,Günther D,Chen L and Hu SH.2002.Determ ination of forty two major and trace elements in USGS and NIST SRM glasses by laser ablation inductively coupled plasmamass spectrometry.Geostandards Newsletter-Journal of Geostandards and Geoanalysis,26(2):191-196

    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 Cosmochimica Acta,64(1):133-147

    Guan Q,Zhu DC,Zhao ZD,Zhang LL,Liu M,Li XW,Yu F and Mo XX.2010.Late Cretaceous adakites from the eastern segment of the Gangdese Belt,southern Tibet:Products of Neo-Tetheyan mid-ocean ridge subduction.Acta Petrologica Sinica,26(7):2165-2179 (in Chinese with English abstract)

    Guan Q,Zhu DC,Zhao ZD,Dong GC,Mo XX,Liu YS,Hu ZC and Yuan HL.2011.Zircon U-Pb chronology,geochemistry of the Late Cretaceous mafic magnatism in the southern Lhasa Terrane and its implications.Acta Petrologica Sinica,27(7):2083- 2094 (in Chinese with English abstract)Guan Q,Zhu DC,Zhao ZD,Dong GC,Zhang LL,Li XW,Liu M,Mo

    XX,Liu YS and Yuan HL.2012.Crustal thickening prior to 38Ma in southern Tibet:Evidence from lower crust-derived adakitic magmatism in the Gangdese Batholith.Gondwana Research,21(1):88-99

    Harris NBW,Xu RH,Lewis CL and Jin CW.1988.Plutonic 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):145-146

    Hoskin PWO and Ireland TR.2000.Rare earth element chemistry of zircon and its use as a provenance indicator.Geology,28(7):627-630

    Hoskin PWO and Schaltegger U.2003.The composition of zircon and igneous and metamorphic petrogenesis.Reviews in Mineralogy and Geochemistry,53(1):27-62

    Hou KJ,Li YH,Zou TR,Qu XM,Shi YR and Xie GQ.2007.Laser ablation-MC-ICP-MS technique for Hf isotope microanalysis of zircon and its geological applications.Acta Petrologica Sinica,23(10):2595-2604 (in Chinese with English abstract)

    Huang Y,Zhao ZD,Zhang FQ,Zhu DC,Dong GC,Zhou S and Mo XX.2010.Geochemistry and implication of the Gangdese batholiths from Renbu and Lhasa areas in southern Gangdese,Tibet.Acta Petrologica Sinica,26(10):3131-3142 (in Chinese with English abstract)

    Ji WQ,Wu FY,Chung SL,Li JX and Liu CZ.2009.Zircon U-Pb geochronology and Hf isotopic constraints on petrogenesis of the Gangdese batholith,southern Tibet.Chemical Geology,262(3-4):229-245

    Ji WQ,Wu FY,Chung SL and Liu CZ.2014.The Gangdese magmatic constraints on a latest Cretaceous lithospheric delamination of the Lhasa terrane,southern Tibet.Lithos,210-211:168-180

    Kang ZQ,Xu JF,Wilde SA,F(xiàn)eng ZH,Chen JL,Wang BD,F(xiàn)u WC and Pan HB.2014.Geochronology and geochemistry of the Sangri Group volcanic rocks,Southern Lhasa Terrane:Implications for the early subduction history of the Neo-Tethys and Gangdese magmatic arc.Lithos,200-201:157-168

    Kapp P,DeCelles PG,Leier AL,F(xiàn)abijanic JM,He SD,Pullen A and Gehrels GE.2007.The Gangdese retroarc thurst belt revealed.GSA Today,17(7):4-9

    Lee HY,Chung SL,Wang YB,Zhu DC,Yang JH,Song B,Liu DY and Wu FY.2007.Age,petrogenesis and geological significance of Linzizong volcanic successions in the Linzhou basin,southern Tibet:Evidence from zircon U-Pb dates and Hf isotopes.Acta Petrologica Sinica,23(2):493-500 (in Chinese with English abstract)

    Lee HY,Chung SL,Lo CH,Ji JQ,Lee TY,Qian Q and Zhang Q.2009.Eocene Neotethyan slab breakoff in southern Tibet inferred from the Linzizong volcanic record.Tectonophysics,477(1-2):20-35

    Leier AL,DeCelles PG,Kapp P and Ding L.2007.The Takena Formation of the Lhasa terrane,southern Tibet:The record of a Late Cretaceous retroarc foreland basin.Geological Society of America Bulletin,119(1-2):31-48

    Li XX,Jiang W,Liang JH,Zhao ZD,Liu D and Mo XX.2014.The geochemical characteristics and significance of the basalt from Shexing Formation in Linzhou basin, southern Tibet.Acta Petrologica Sinica,31(5):1285-1297 (in Chinese with English abstract)

    Liu HF.1993.Division of Linzizong volcanic rock system and belong to time in Lhasa area.Tibet Geology,(2):59-69 (in Chinese)

    Liu YS,Hu ZC,Gao S,Günther 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

    Liu YS,Gao S,Hu ZC,Gao CG,Zong KQ and Wang DB.2010a.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 of mantle xenoliths.Journal of Petrology,51(1-2):537-571

    Liu YS,Hu ZC,Zong KQ,Gao CG,Gao S,Xu J and Chen HH.2010b.Reappraisement and refinement of zircon U-Pb isotope and trace element analyses by LA-ICP-MS.Chinese Science Bulletin,55(15):1535-1546

    Ma L,Wang Q,Li ZX,Wyman DA,Jiang ZQ,Yang JH,Gou GN and Guo HF.2013a.Early Late Cretaceous (ca.93Ma)norites and hornblendites in the Milin area,eastern Gangdese:Lithosphereasthenosphere interaction during slab rollback and an insight into early Late Cretaceous (ca.100~80Ma)magmatic“flare-up”in southern Lhasa (Tibet).Lithos,172-173:17-30 Ma L,Wang Q,Wyman DA,Jiang ZQ,Yang JH,Li QL,Gou GN and

    Guo HF.2013b.Late Cretaceous crustal growth in the Gangdese area,southern Tibet:Petrological and Sr-Nd-Hf-O isotopic evidence from Zhengga diorite-gabbro.Chemical Geology,349-350:54-70 Ma L,Wang Q,Wyman DA,Li ZX,Jiang ZQ,Yang JH,Gou GN and Guo HF.2013c.Late Cretaceous (100~89Ma) magnesian charnockites with adakitic affinities in the Milin area,eastern Gangdese: Partial melting of subducted oceanic crust and implications for crustal growth in southern Tibet.Lithos,175-176:315-332

    Meng FY,Zhao ZD,Zhu DC,Zhang LL,Guan Q,Liu M,Yu F,Liu MH and Mo XX.2010.The petrogenesis of Late Cretaceous adakitelike rocks in Mamba from the eastern Gangdese,Tibet.Acta Petrologica Sinica,26(7):2180-2192 (in Chinese with English abstract)

    Meng FY,Zhao Z,Zhu DC,Mo X,Guan Q,Huang Y,Dong GC,Zhou S,DePaolo DJ,Harrison TM,Zhang ZC,Liu JL,Liu YS,Hu ZC and Yuan HL.2014.Late Cretaceous magmatism in Mamba area,central Lhasa subterrane:Products of back-arc extension of Neo-Tethyan Ocean?Gondwana Research,26(2):505-520

    Mo XX,Zhao ZD,Deng JF,Dong GC,Zhou S,Guo TY,Zhang SQ and Wang LL.2003.Response of volcanism to the India-Asia collision.Earth Science Frontiers,10(3):135- 148 (in Chinese with English abstract)

    Mo XX,Hou ZQ,Niu Y,Dong GC,Qu XM,Zhao ZD and Yang ZM.2007.Mantle contributions to crustal thickening during continental collision:Evidence from Cenozoic igneous rocks in southern Tibet.Lithos,96(1-2):225-242

    Mo XX,Niu Y,Dong GC,Zhao ZD,Hou ZQ,Zhou S and Ke S.2008.Contribution of syncollisional felsic magmatism to continental crust growth:A case study of the Paleogene Linzizong volcanic succession in southern Tibet.Chemical Geology,250(1-4):49-67

    Morel MLA,Nebel O,Nebel-Jacobsen YJ,Miller JS and Vroon PZ.2008.Hafnium isotope characterization of the GJ-1 zircon reference material by solution and laser-ablation MC-ICPMS.Chemical Geology,255(1-2):231-235

    Murphy MA,Harrison TM,Durr SB,Chen Z,Ryerson FJ,Kidd WSF,Wang X and Zhou X.1997.Significant crustal shortening in southcentral Tibet prior to the Indo-Asian collision.Geology,25:719-722

    Pearce JA, Harris NBW and Tindle AG.1984.Trace element discrimination diagrams for the tectonic interpretation of granitic rocks.Journal of Petrology,25(4):956-983

    Poller U,Huth J,Hoppe P and Williams IS.2001.REE,U,Th,and Hf distribution in zircon from western Carpathian Variscan granitoids:A combined cathodoluminescence and ion microprobe study.American Journal of Science,301(10):858-876

    Rollinson HR.1993.Using Geochemical Data: Evaluation,Presentation,Interpretation.New York:Longman Group UK Ltd.,1-352

    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 Planetary Science Letters,219(3-4):311-324

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

    Wang L,Zeng LS,Gao LE and Chen ZY.2013.Early Cretaceous high Mg#and high Sr/Y clinopyroxene-bearing diorite in the southeast Gangdese batholith,Southern Tibet.Acta Petrologica Sinica,29(6):1977-1994 (in Chinese with English abstract)

    Wang Q,Zhu DC,Zhao ZD,Guan Q,Zhang XQ,Sui QL,Hu ZC and Mo XX.2012.Magmatic zircons from I-,S-and A-type granitoids in Tibet:Trace element characteristics and their application to detrital zircon provenance study.Journal of Asian Earth Sciences,53:59-66

    Wang W,Zeng LS,Liu J,Xiao P and Gao LE.2013.Geochemical characteristics of the Late-Cretaceous dacite in Coqen,western Gangdese batholith,Tibet.Chinese Journal of Geology,48(2):484-500 (in Chinese with English abstract)

    Wang YL,Zhang CJ and Xiu SZ.2001.Th/Hf-Ta/Hf identification of tectonic setting of basalts.Acta Petrologica Sinica,17(3):413-421 (in Chinese with English abstract)

    Wen DR,Chung SL,Song B,Iizuka Y,Yang HJ,Ji J,Liu D and Gallet S.2008a.Late Cretaceous Gangdese intrusions of adakitic geochemical characteristics,SE Tibet:Petrogenesis and tectonic implication.Lithos,105(1-2):1-11

    Wen DR,Liu DY,Chung SL,Chu MF,Ji JQ,Zhang Q,Song B,Lee TY,Yeh MW and Lo CH.2008b.Zircon SHRIMP U-Pb ages of the Gangdese batholith and implications for Neotethyan subduction in southern Tibet.Chemical Geology,252(3-4):191-201

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

    Yao P,Li JG,Wang QH,Gu XX,Tang JX and Hui L.2006.Discovery and geological significance of the adakite in Gangdise island arc belt,Xizang (Tibet).Acta Petrologica Sinica,22(3):612-620 (in Chinese with English abstract)

    Yin A and Harrison TM.2000.Geologic evolution of the Himalayan-Tibetan Orogen.Annual Review of Earth and Planetary Sciences,28:211-280

    Zhang ZM,Zhao GC,Santosh M,Wang JL,Dong X and Shen K.2010.Late Cretaceous charnockite with adakitic affinities from the Gangdese batholith,southeastern Tibet:Evidence for Neo-Tethyan mid-ocean ridge subduction?Gondwana Research,17 (4):615-631

    Zhao ZD,Mo XX,Dilek Y,Niu Y,DePaolo DJ,Robinson P,Zhu DC,Sun CG,Dong GC,Zhou S,Luo ZH and Hou ZQ.2009.Geochemical and Sr-Nd-Pb-O isotopic compositions of the postcollisional ultrapotassic magmatism in SW Tibet:Petrogenesis and implications for India intra-continental subduction beneath southern Tibet.Lithos,113(1-2):190-212

    Zheng YC,Hou ZQ,Gong YL,Liang W,Sun QZ,Zhang S,F(xiàn)u Q,Huang KX,Li QY and Li W.2014.Petrogenesis of Cretaceous adakite-like intrusions of the Gangdese Plutonic Belt,southern Tibet:Implications for mid-ocean ridge subduction and crustal growth.Lithos,190-191:240-263

    Zhou S,Mo XX,Dong GC,Zhao ZD,Qiu RZ,Guo TY and Wang LL.2004.40Ar-39Ar geochronology of Cenozoic Linzizong volcanic rocks from Linzhou Basin,Tibet,China,and their geological implications.Chinese Science Bulletin,49(18):1970-1979

    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,27(9):1535-1550 (in Chinese with English abstract)

    Zhu DC,Mo XX,Niu Y,Zhao ZD,Wang LQ,Liu YS and Wu FY.2009.Geochemical investigation of Early Cretaceous igneous rocks along an east-west traverse throughout the central Lhasa Terrane,Tibet.Chemical Geology,268(3-4):298-312

    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 and Planetary 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

    附中文參考文獻

    常承法.1984.雅魯藏布江縫合帶地質(zhì)構(gòu)造特征及其演化.見:中法喜馬拉雅考察報告.北京:科學(xué)出版社,327-339

    董彥輝,許繼峰,曾慶高,王強,毛國政,李杰.2006.存在比桑日群弧火山巖更早的新特提斯洋俯沖記錄么?巖石學(xué)報,22(3):661-668

    管琪,朱弟成,趙志丹,張亮亮,劉敏,李小偉,于楓,莫宣學(xué).2010.西藏南部岡底斯帶東段晚白堊世埃達(dá)克巖:新特提斯洋脊俯沖的產(chǎn)物?巖石學(xué)報,26(7):2165-2179

    管琪,朱弟成,趙忠丹,董國臣,莫宣學(xué),劉勇勝,胡兆初,袁洪林.2011.西藏拉薩地塊南緣晚白堊世鎂鐵質(zhì)巖漿作用的年代學(xué)、地球化學(xué)及意義.巖石學(xué)報,27(7):2083-2094

    侯可軍,李延河,鄒天人,曲曉明,石玉若,謝桂青.2007.LA-MCICP-MS 鋯石Hf 同位素的分析方法及地質(zhì)應(yīng)用.巖石學(xué)報,23(10):2595-2604

    黃玉,趙志丹,張鳳琴,朱弟成,董國臣,周肅,莫宣學(xué).2010.西藏岡底斯仁布-拉薩一帶花崗巖基的地球化學(xué)及其意義.巖石學(xué)報,26(10):3131-3142

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

    李曉雄,江萬,梁錦海,趙志丹,劉棟,莫宣學(xué).2015.西藏林周盆地設(shè)興組玄武巖地球化學(xué)特征及意義.巖石學(xué)報,31(5):1285-1297

    劉鴻飛.1993.拉薩地區(qū)林子宗火山巖系的劃分和時代歸屬.西藏地質(zhì),(2):59-69

    孟繁一,趙志丹,朱弟成,張亮亮,管琪,劉敏,于楓,劉美華,莫宣學(xué).2010.西藏岡底斯東部門巴地區(qū)晚白堊世埃達(dá)克質(zhì)巖的巖石成因.巖石學(xué)報,26(7):2180-2192

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

    王莉,曾令森,高利娥,陳振宇.2013.藏南岡底斯巖基東南緣早白堊世高鎂-高Sr/Y 含單斜輝石閃長巖.巖石學(xué)報,29(6):1977-1994

    王偉,曾令森,劉靜,肖萍,高利娥.2013.西藏措勤晚白堊世英安巖的厘定與地球化學(xué)特征.地質(zhì)科學(xué),48(2):484-500

    汪云亮,張成江,修淑芝.2001.玄武巖類形成的大地構(gòu)造環(huán)境的Th/Hf-Ta/Hf 圖解判別.巖石學(xué)報,17(3):413-421

    姚鵬,李金高,王全海,顧雪祥,唐菊興,惠蘭.2006.西藏岡底斯南緣火山-巖漿弧帶中桑日群adakite 的發(fā)現(xiàn)及其意義.巖石學(xué)報,22(3):612-620

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

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

    猜你喜歡
    輝綠巖白堊鋯石
    鋯石成因礦物學(xué)及Lu-Hf同位素的應(yīng)用
    SN 地區(qū)輝綠巖下弱信號寬方位采集與成像
    阜平輝綠巖飾面石材礦廢石利用簡析
    微量元素在鶴崗盆地早白堊世沉積環(huán)境分析中的應(yīng)用
    柴北緣牦牛山地區(qū)輝綠巖體鋯石U-Pb年齡及地質(zhì)意義
    俄成功試射“鋯石”高超音速巡航導(dǎo)彈
    軍事文摘(2020年24期)2020-02-06 05:56:36
    新疆博格達(dá)地區(qū)發(fā)現(xiàn)類似稀土四分組效應(yīng)的輝綠巖
    鋯石微區(qū)原位U-Pb定年的測定位置選擇方法
    治白堊病良方
    鉆石與鋯石的區(qū)別知多少?
    精品亚洲乱码少妇综合久久| 国产黄频视频在线观看| 日韩不卡一区二区三区视频在线| 国产高清有码在线观看视频| 亚洲精品中文字幕在线视频 | 国产精品成人在线| 亚洲精品国产成人久久av| av不卡在线播放| 伊人久久国产一区二区| 久久国产精品大桥未久av | 观看免费一级毛片| 老司机影院毛片| 欧美bdsm另类| 国产精品久久久久久av不卡| 丰满饥渴人妻一区二区三| 国产免费又黄又爽又色| 人体艺术视频欧美日本| 18禁裸乳无遮挡动漫免费视频| 在线观看免费高清a一片| 亚洲成人一二三区av| 国产又色又爽无遮挡免| 久久久亚洲精品成人影院| 久久久久视频综合| 在线 av 中文字幕| 99热这里只有是精品50| 久久99一区二区三区| 亚洲欧美成人综合另类久久久| 亚洲精品久久午夜乱码| 国产精品蜜桃在线观看| 成人免费观看视频高清| 桃花免费在线播放| av不卡在线播放| 观看免费一级毛片| 永久免费av网站大全| 日韩三级伦理在线观看| 曰老女人黄片| kizo精华| 自拍偷自拍亚洲精品老妇| 亚洲欧洲精品一区二区精品久久久 | 国产一区有黄有色的免费视频| 五月开心婷婷网| 国产欧美日韩一区二区三区在线 | 亚洲精品自拍成人| 久久影院123| 男人爽女人下面视频在线观看| 久久久国产欧美日韩av| 亚洲国产精品一区二区三区在线| 免费久久久久久久精品成人欧美视频 | 永久网站在线| √禁漫天堂资源中文www| 五月天丁香电影| 欧美精品高潮呻吟av久久| 99视频精品全部免费 在线| 日本wwww免费看| 国产免费视频播放在线视频| 欧美成人午夜免费资源| 精品国产一区二区三区久久久樱花| 日本wwww免费看| 国产精品久久久久久av不卡| 精品亚洲成国产av| 天天操日日干夜夜撸| 桃花免费在线播放| 亚洲av免费高清在线观看| 精品人妻一区二区三区麻豆| 天天操日日干夜夜撸| 国产男女内射视频| 亚洲av免费高清在线观看| 成人美女网站在线观看视频| 亚洲国产精品一区二区三区在线| 精品午夜福利在线看| 欧美 亚洲 国产 日韩一| 亚洲av在线观看美女高潮| 久久女婷五月综合色啪小说| 视频区图区小说| 一二三四中文在线观看免费高清| 亚洲av在线观看美女高潮| 亚洲成人av在线免费| 亚洲国产精品专区欧美| 99九九在线精品视频 | 色哟哟·www| 亚洲欧美中文字幕日韩二区| 少妇的逼水好多| 男女啪啪激烈高潮av片| 亚洲精品久久午夜乱码| 精品亚洲成国产av| 亚洲精品视频女| av天堂久久9| 国产69精品久久久久777片| 国产伦理片在线播放av一区| 国产伦精品一区二区三区视频9| 99久久综合免费| 国产欧美日韩一区二区三区在线 | 免费看日本二区| 久久久久久久久大av| 亚洲精品自拍成人| 日韩不卡一区二区三区视频在线| 亚洲精品国产成人久久av| 日韩成人av中文字幕在线观看| 久久精品久久精品一区二区三区| 黄色视频在线播放观看不卡| 黑人猛操日本美女一级片| 成人国产av品久久久| 啦啦啦中文免费视频观看日本| 少妇猛男粗大的猛烈进出视频| 99久国产av精品国产电影| 51国产日韩欧美| 人妻系列 视频| 丝瓜视频免费看黄片| 国产精品人妻久久久久久| 国产高清不卡午夜福利| 99热网站在线观看| 能在线免费看毛片的网站| 丁香六月天网| 一级,二级,三级黄色视频| 国产高清国产精品国产三级| 久热久热在线精品观看| 永久网站在线| 日韩成人av中文字幕在线观看| 熟妇人妻不卡中文字幕| 两个人的视频大全免费| 久久青草综合色| 亚洲美女搞黄在线观看| 午夜福利网站1000一区二区三区| 色视频在线一区二区三区| 女的被弄到高潮叫床怎么办| 欧美最新免费一区二区三区| 日韩中文字幕视频在线看片| 久久精品久久久久久久性| 久久久久国产网址| 久久人人爽人人爽人人片va| 免费看av在线观看网站| 精品国产露脸久久av麻豆| 一本—道久久a久久精品蜜桃钙片| av不卡在线播放| 久久毛片免费看一区二区三区| 午夜av观看不卡| 一二三四中文在线观看免费高清| 精品视频人人做人人爽| 国产 一区精品| 国产免费视频播放在线视频| 成人午夜精彩视频在线观看| 日韩伦理黄色片| 午夜91福利影院| 黄色欧美视频在线观看| 亚洲精品,欧美精品| 91精品国产国语对白视频| 精品亚洲成a人片在线观看| 成人毛片60女人毛片免费| 中国国产av一级| 婷婷色av中文字幕| 黑丝袜美女国产一区| 中文精品一卡2卡3卡4更新| 国产片特级美女逼逼视频| 成人午夜精彩视频在线观看| 精品久久久久久久久亚洲| 欧美日韩国产mv在线观看视频| 成人综合一区亚洲| 国产一区二区三区av在线| 嫩草影院入口| 一级毛片aaaaaa免费看小| 久久99精品国语久久久| 黄色毛片三级朝国网站 | 麻豆精品久久久久久蜜桃| 亚洲在久久综合| 七月丁香在线播放| 丰满饥渴人妻一区二区三| 人人妻人人澡人人爽人人夜夜| 亚洲精品色激情综合| 久久精品夜色国产| 美女脱内裤让男人舔精品视频| 精品久久久久久久久av| 成年av动漫网址| 亚洲欧美精品自产自拍| 嫩草影院新地址| 久久精品国产a三级三级三级| 久久人妻熟女aⅴ| 色哟哟·www| 亚洲av综合色区一区| 夫妻性生交免费视频一级片| 一级二级三级毛片免费看| 黑人巨大精品欧美一区二区蜜桃 | 国产在线免费精品| av一本久久久久| 国产精品久久久久成人av| a级片在线免费高清观看视频| 美女主播在线视频| 少妇 在线观看| 成人免费观看视频高清| 亚洲国产精品一区三区| 国产成人精品一,二区| 国产极品天堂在线| 久久6这里有精品| 亚洲第一av免费看| 搡女人真爽免费视频火全软件| 国产老妇伦熟女老妇高清| 麻豆成人av视频| av免费在线看不卡| 亚洲国产欧美日韩在线播放 | 狂野欧美白嫩少妇大欣赏| 少妇高潮的动态图| 噜噜噜噜噜久久久久久91| 亚洲av在线观看美女高潮| 一级,二级,三级黄色视频| 在线观看免费高清a一片| 在线观看免费视频网站a站| 久久精品国产自在天天线| 国产精品嫩草影院av在线观看| 亚洲久久久国产精品| 最近最新中文字幕免费大全7| 亚洲欧美日韩另类电影网站| 高清午夜精品一区二区三区| 人人澡人人妻人| 丝袜喷水一区| 免费播放大片免费观看视频在线观看| 国产亚洲一区二区精品| 夜夜看夜夜爽夜夜摸| 亚洲av电影在线观看一区二区三区| 欧美 日韩 精品 国产| 亚洲av日韩在线播放| 天美传媒精品一区二区| 国产av一区二区精品久久| 成人国产麻豆网| 国产成人aa在线观看| 大香蕉久久网| 国产男人的电影天堂91| 伊人久久国产一区二区| 3wmmmm亚洲av在线观看| 建设人人有责人人尽责人人享有的| 亚洲成人av在线免费| 久久影院123| 一个人免费看片子| 国产欧美日韩综合在线一区二区 | 大陆偷拍与自拍| 91久久精品国产一区二区三区| 丝瓜视频免费看黄片| 午夜日本视频在线| 在线 av 中文字幕| 国产精品久久久久久久电影| 91久久精品国产一区二区成人| 成人综合一区亚洲| 极品人妻少妇av视频| 亚洲国产色片| 国产黄片美女视频| 欧美另类一区| 男女无遮挡免费网站观看| 26uuu在线亚洲综合色| 久久精品夜色国产| 午夜福利网站1000一区二区三区| 一本久久精品| 成人免费观看视频高清| av天堂久久9| 亚洲国产毛片av蜜桃av| 亚洲欧洲国产日韩| 晚上一个人看的免费电影| 日韩人妻高清精品专区| 2021少妇久久久久久久久久久| 日日摸夜夜添夜夜添av毛片| 国国产精品蜜臀av免费| 大又大粗又爽又黄少妇毛片口| 精品国产乱码久久久久久小说| 久久久久国产精品人妻一区二区| 男女免费视频国产| 大片免费播放器 马上看| 色5月婷婷丁香| 亚洲国产欧美在线一区| 午夜精品国产一区二区电影| 汤姆久久久久久久影院中文字幕| 老司机影院成人| av有码第一页| 狂野欧美激情性xxxx在线观看| av在线老鸭窝| 精品一区二区免费观看| 亚洲av综合色区一区| 永久网站在线| 亚洲综合精品二区| 国产精品一区二区在线观看99| 青青草视频在线视频观看| 久久久久久久亚洲中文字幕| 婷婷色综合大香蕉| 久久精品久久久久久久性| 天天操日日干夜夜撸| 国产免费一区二区三区四区乱码| 九草在线视频观看| 丰满乱子伦码专区| 一个人免费看片子| 成人国产av品久久久| 又爽又黄a免费视频| 汤姆久久久久久久影院中文字幕| 下体分泌物呈黄色| 人妻夜夜爽99麻豆av| av在线播放精品| 午夜免费男女啪啪视频观看| 日日撸夜夜添| 亚洲av在线观看美女高潮| 欧美精品一区二区免费开放| 国产亚洲一区二区精品| 男人舔奶头视频| 亚洲四区av| 国产精品一区二区三区四区免费观看| 黄色怎么调成土黄色| 久久这里有精品视频免费| 日韩精品有码人妻一区| 亚洲av不卡在线观看| 日韩一本色道免费dvd| 亚洲精品国产av成人精品| 成年人免费黄色播放视频 | 国产深夜福利视频在线观看| 街头女战士在线观看网站| 亚洲国产精品成人久久小说| 两个人的视频大全免费| 国产亚洲精品久久久com| 欧美最新免费一区二区三区| 国产av国产精品国产| 亚洲无线观看免费| 高清av免费在线| 亚洲图色成人| 国产一区二区在线观看av| 美女视频免费永久观看网站| 欧美日韩在线观看h| 秋霞在线观看毛片| 天天躁夜夜躁狠狠久久av| 色婷婷久久久亚洲欧美| 一级片'在线观看视频| 国产爽快片一区二区三区| 欧美日韩视频精品一区| 韩国av在线不卡| 欧美精品一区二区大全| 在线 av 中文字幕| 边亲边吃奶的免费视频| 人妻人人澡人人爽人人| 黄色日韩在线| 热re99久久国产66热| 亚洲精品一二三| 国产男人的电影天堂91| 精品久久久噜噜| 免费播放大片免费观看视频在线观看| 一级av片app| 男女边摸边吃奶| 久久久久网色| 成年美女黄网站色视频大全免费 | 丝瓜视频免费看黄片| 亚洲熟女精品中文字幕| 久久久久久久大尺度免费视频| 91久久精品电影网| 免费观看的影片在线观看| 最近中文字幕2019免费版| 亚洲欧美日韩卡通动漫| 国产欧美亚洲国产| 卡戴珊不雅视频在线播放| 久久精品国产亚洲av涩爱| 亚洲情色 制服丝袜| 久久久欧美国产精品| 亚洲内射少妇av| 一级爰片在线观看| 亚洲精品乱久久久久久| 欧美97在线视频| 日韩大片免费观看网站| 久久久久久久亚洲中文字幕| 亚洲自偷自拍三级| 免费少妇av软件| 18禁裸乳无遮挡动漫免费视频| 99久久人妻综合| 亚洲国产av新网站| 99久久人妻综合| 久久精品国产自在天天线| 欧美日韩亚洲高清精品| 亚洲精品,欧美精品| a级毛色黄片| 天天躁夜夜躁狠狠久久av| 麻豆成人午夜福利视频| 亚洲精品自拍成人| 麻豆乱淫一区二区| 国产 一区精品| 欧美丝袜亚洲另类| 如日韩欧美国产精品一区二区三区 | 久久99热这里只频精品6学生| 亚洲精品一二三| 亚洲欧洲精品一区二区精品久久久 | 内地一区二区视频在线| 日本色播在线视频| 美女视频免费永久观看网站| 亚洲欧洲国产日韩| 亚洲精品一区蜜桃| www.色视频.com| 国产深夜福利视频在线观看| 亚洲欧美成人综合另类久久久| 纯流量卡能插随身wifi吗| 少妇被粗大猛烈的视频| 午夜精品国产一区二区电影| 亚洲美女视频黄频| 老女人水多毛片| 亚洲综合色惰| 国产精品麻豆人妻色哟哟久久| 国产精品国产三级国产av玫瑰| 国产精品麻豆人妻色哟哟久久| 亚洲一级一片aⅴ在线观看| 中文在线观看免费www的网站| 婷婷色麻豆天堂久久| 大香蕉97超碰在线| 尾随美女入室| 免费不卡的大黄色大毛片视频在线观看| 男人狂女人下面高潮的视频| 日本黄大片高清| 免费观看性生交大片5| 丁香六月天网| av在线观看视频网站免费| 久久韩国三级中文字幕| 中文字幕久久专区| 精华霜和精华液先用哪个| 99久久精品国产国产毛片| 免费看光身美女| 中文精品一卡2卡3卡4更新| 妹子高潮喷水视频| 在线天堂最新版资源| 搡老乐熟女国产| av黄色大香蕉| 国产精品国产三级专区第一集| 国内精品宾馆在线| 色5月婷婷丁香| 精品午夜福利在线看| 美女视频免费永久观看网站| 两个人的视频大全免费| 国产淫片久久久久久久久| 高清黄色对白视频在线免费看 | 只有这里有精品99| 国产一级毛片在线| 少妇丰满av| 又爽又黄a免费视频| 中国三级夫妇交换| 久久久久久久亚洲中文字幕| 国国产精品蜜臀av免费| 午夜福利,免费看| 有码 亚洲区| 国产午夜精品一二区理论片| 9色porny在线观看| 日本欧美视频一区| 亚洲综合色惰| 国产精品久久久久久久久免| 日韩av在线免费看完整版不卡| 老司机影院成人| 精品国产露脸久久av麻豆| 日韩成人av中文字幕在线观看| 成人综合一区亚洲| 中国三级夫妇交换| 久久精品国产亚洲av天美| 69精品国产乱码久久久| 久久久久国产精品人妻一区二区| 99视频精品全部免费 在线| 亚洲精品国产成人久久av| 久久久国产欧美日韩av| 又大又黄又爽视频免费| 麻豆乱淫一区二区| 欧美xxⅹ黑人| 性色avwww在线观看| 亚洲在久久综合| 国产高清国产精品国产三级| 午夜激情久久久久久久| 一级毛片我不卡| 国产成人午夜福利电影在线观看| 丰满人妻一区二区三区视频av| 中文字幕亚洲精品专区| 亚洲激情五月婷婷啪啪| 亚洲欧美清纯卡通| 大香蕉久久网| 天堂8中文在线网| 秋霞伦理黄片| av又黄又爽大尺度在线免费看| 内地一区二区视频在线| 亚洲精品日韩av片在线观看| 99热6这里只有精品| 精品熟女少妇av免费看| 精品国产国语对白av| 久久鲁丝午夜福利片| 久久影院123| 免费看不卡的av| a级毛片免费高清观看在线播放| 蜜桃在线观看..| 日韩成人av中文字幕在线观看| 国产白丝娇喘喷水9色精品| 亚洲av免费高清在线观看| 一级,二级,三级黄色视频| 国产日韩欧美视频二区| 3wmmmm亚洲av在线观看| 久久久国产精品麻豆| 老司机亚洲免费影院| 我的老师免费观看完整版| 人妻系列 视频| 日本色播在线视频| 十八禁高潮呻吟视频 | 午夜免费观看性视频| 狂野欧美激情性xxxx在线观看| 一个人看视频在线观看www免费| 亚洲成人手机| 欧美激情极品国产一区二区三区 | 国产视频首页在线观看| 久久久久国产网址| 丰满少妇做爰视频| 日韩一本色道免费dvd| 欧美3d第一页| 少妇人妻久久综合中文| 欧美变态另类bdsm刘玥| 视频中文字幕在线观看| 99精国产麻豆久久婷婷| 午夜激情福利司机影院| 日本午夜av视频| 一级爰片在线观看| 中国国产av一级| 女性生殖器流出的白浆| 亚洲精品日本国产第一区| 男女啪啪激烈高潮av片| 18禁动态无遮挡网站| 五月玫瑰六月丁香| 18+在线观看网站| 男男h啪啪无遮挡| 99精国产麻豆久久婷婷| 欧美xxxx性猛交bbbb| 免费播放大片免费观看视频在线观看| 国产日韩欧美在线精品| 久久午夜综合久久蜜桃| 大陆偷拍与自拍| 丝袜脚勾引网站| 久久国产乱子免费精品| 成人美女网站在线观看视频| 搡老乐熟女国产| 老司机影院毛片| 精品一区二区三区视频在线| 欧美日韩在线观看h| av在线老鸭窝| 成人影院久久| 久久国产精品大桥未久av | 国产 一区精品| 日本91视频免费播放| 少妇的逼好多水| 国产成人freesex在线| 大香蕉97超碰在线| 国产极品天堂在线| 中文在线观看免费www的网站| 亚洲国产欧美在线一区| 久久久久久久久久成人| 亚洲欧美清纯卡通| 欧美亚洲 丝袜 人妻 在线| 欧美成人精品欧美一级黄| 精华霜和精华液先用哪个| 男人狂女人下面高潮的视频| 日日撸夜夜添| 熟女av电影| 人人妻人人添人人爽欧美一区卜| 亚洲性久久影院| 亚洲成人手机| 国产极品天堂在线| av线在线观看网站| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产精品国产精品| 国产精品三级大全| 日本黄色日本黄色录像| 男人爽女人下面视频在线观看| 亚洲天堂av无毛| 精品久久久噜噜| 日日爽夜夜爽网站| 亚洲欧美精品自产自拍| 特大巨黑吊av在线直播| 七月丁香在线播放| 熟女av电影| 国产精品久久久久成人av| 免费高清在线观看视频在线观看| 精品人妻一区二区三区麻豆| 国产精品女同一区二区软件| 人体艺术视频欧美日本| 成人二区视频| 各种免费的搞黄视频| 伦精品一区二区三区| av不卡在线播放| 中文欧美无线码| 精品一区二区三卡| 久久99热这里只频精品6学生| 男女免费视频国产| 国产爽快片一区二区三区| 国产国拍精品亚洲av在线观看| 国产欧美亚洲国产| 精品人妻熟女av久视频| 久久精品久久久久久噜噜老黄| 国产欧美日韩精品一区二区| 美女福利国产在线| 久久久亚洲精品成人影院| 建设人人有责人人尽责人人享有的| 国产黄频视频在线观看| 亚洲成人一二三区av| 一级爰片在线观看| 精品一区二区免费观看| 国产日韩欧美在线精品| av一本久久久久| 国产精品.久久久| 亚洲精品国产av成人精品| 亚洲av日韩在线播放| 丝袜脚勾引网站| 91aial.com中文字幕在线观看| 纯流量卡能插随身wifi吗| 国产成人免费观看mmmm| 男女无遮挡免费网站观看| 美女xxoo啪啪120秒动态图| 欧美人与善性xxx| 国产精品久久久久久av不卡| 国产有黄有色有爽视频| 亚洲精品国产av成人精品| 边亲边吃奶的免费视频| 国产69精品久久久久777片| 日本午夜av视频| 久久精品国产a三级三级三级| 亚洲综合色惰| 国产色婷婷99| 成人特级av手机在线观看| 欧美变态另类bdsm刘玥| 久久午夜综合久久蜜桃| 九九在线视频观看精品| 色视频www国产|