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

    海南島新生代玄武巖成因的PGE與Os同位素地球化學(xué)制約

    2022-07-08 07:13:50王國慶
    大地構(gòu)造與成礦學(xué) 2022年3期
    關(guān)鍵詞:洋殼源區(qū)海南島

    王國慶, 郭 鋒, 趙 亮

    海南島新生代玄武巖成因的PGE與Os同位素地球化學(xué)制約

    王國慶1, 2, 郭 鋒1, 3*, 趙 亮1, 3

    (1. 中國科學(xué)院 廣州地球化學(xué)研究所, 同位素地球化學(xué)國家重點實驗室, 廣東 廣州 510640; 2. 中國科學(xué)院大學(xué) 地球與行星科學(xué)學(xué)院, 北京 100049; 3. 中國科學(xué)院深地科學(xué)卓越研究中心, 廣東 廣州 510640)

    本文對海南島新生代玄武巖開展了PGE和Re-Os同位素體系的系統(tǒng)研究, 結(jié)合全巖主量、微量元素和Sr-Nd-Pb同位素數(shù)據(jù), 取得了如下認(rèn)識: ①海南島新生代玄武巖根據(jù)PGE含量和配分模式可以劃分為兩類, 第一類PGE總量為1.88×10?9~4.54×10?9, 配分模式為相對平坦型(Pd/Ir<3.5), 高Os含量和低187Os/188Os值, 代表了來自虧損軟流圈地幔的熔體; 第二類PGE總量為0.35×10?9~2.67×10?9, 配分模式為向右上傾型(Pd/Ir>5.0), 低Os含量和高187Os/188Os值, 反映其熔融地幔源區(qū)包含了大量的再循環(huán)洋殼; ②海南島新生代玄武巖PGE地球化學(xué)特征受控于巖漿演化過程中硫化物的熔離作用; ③Os-Pb同位素模擬結(jié)果暗示海南島新生代玄武巖再循環(huán)洋殼的年齡在100~300 Ma之間, 為近期/同期俯沖的大洋板片。本次研究為進(jìn)一步理解海南島新生代玄武巖源區(qū)的再循環(huán)組分來源、硫化物在巖漿起源與演化的作用以及地幔柱演化動力學(xué)方面提供了新的依據(jù)。

    再循環(huán)洋殼; PGE; Re-Os同位素; 新生代玄武巖; 海南島

    0 引 言

    海南島及周邊地區(qū)包括雷州半島、南海海盆、北部灣地區(qū)、印支半島等地廣泛分布新生代玄武巖(圖1b), 前人對這些玄武巖進(jìn)行了大量研究。Hoang and Flower (1998)認(rèn)為在海南島周緣可能存在一個分散的大火成巖省, 其成因與地幔柱作用相關(guān)(Wang et al., 2012, 2013; Xia et al., 2016)。已有的年代學(xué)和巖石地球化學(xué)研究表明, 這些玄武巖主要噴發(fā)于晚新生代(Ho et al., 2000; 韓江偉等, 2009; Wang et al., 2012), 巖石類型包括拉斑玄武巖和堿性玄武巖, 微量元素蛛網(wǎng)圖特征與典型洋島型玄武巖(OIB)一致, 全巖Sr-Nd-Pb-Hf同位素結(jié)果指示玄武巖源區(qū)具有兩端元混合的特征, 即DMM和EM2, 且具有Dupal異常(Tu et al., 1992; Zou and Fan, 2010; Wang et al., 2012, 2013; Liu et al., 2015; 梅盛旺和任鐘元, 2019)。然而, 關(guān)于EM2的來源一直存在爭議, 有人認(rèn)為其可能來源于陸下巖石圈地幔(Tu et al., 1992; Hoang et al., 1996),也有學(xué)者認(rèn)為源于再循環(huán)到核?幔邊界的太平洋洋殼(Wang et al., 2013; Li et al., 2020)。

    Re和Os具有強(qiáng)烈的親鐵和親銅性質(zhì), 在地球早期形成和分異過程中富集于地核, 而在地幔中含量較低。187Re可通過衰變形成187Os。在地幔熔融過程中, Re和Os表現(xiàn)出不同的地球化學(xué)行為, Os是相容元素, 易殘留在地幔中; Re是中等不相容元素, 傾向進(jìn)入熔體及流體相中?;赗e和Os的地球化學(xué)性質(zhì), Re-Os同位素體系在大陸巖石圈地幔定年、示蹤殼幔演化和超基性?基性巖成因研究等方面都具有重要意義(Shirey and Walker, 1998; Gao et al., 2002)。

    鉑族元素(PGE)包括鋨(Os)、銥(Ir)、釕(Ru)、銠(Rh)、鉑(Pt)和鈀(Pd)6個元素, 具高度親硫性(Barnes et al., 1985; Barnes, 1993)。根據(jù)耐熔性和化學(xué)性質(zhì), 鉑族元素可分為兩個亞組, 即相對易熔的PPGE(Rh、Pt、Pd)和難熔的IPGE(Os、Ir、Ru)(Pearson et al., 2004; Becker et al., 2006)。由于IPGE元素相容性遠(yuǎn)大于PPGE, 在地幔分異過程中, Os、Ir、Ru傾向于保留在地幔殘留相, 而Rh、Pt、Pd傾向進(jìn)入到熔體相。因此, 經(jīng)歷了高度熔體提取的難熔克拉通型地幔具有PPGE相對IPGE虧損的特征, 而新生的飽滿型巖石圈地幔具有相對平坦的PGE配分模式(Lorand et al., 2013); 熔體相(玄武巖)則表現(xiàn)出IPGE相對PPGE的虧損配分模式。由于鉑族元素的親硫親銅性, 它們?nèi)菀走M(jìn)入到硫化物中, 因此可以用PGE配分模式來判斷熔體中硫的飽和程度?;阢K族元素具有不同于其他微量元素的特殊化學(xué)性質(zhì),運用PGE地球化學(xué)行為來示蹤鎂鐵質(zhì)巖漿起源和演化過程中硫化物的貢獻(xiàn), 為理解深部動力學(xué)過程提供依據(jù)。

    圖b據(jù)Yan et al., 2018修改, 年齡數(shù)據(jù)來自Ho et al., 2000; 圖c據(jù)Wang et al., 2012修改。

    本次對海南島新生代玄武巖開展PGE和Re-Os同位素體系綜合研究, 結(jié)合全巖的主量、微量元素和Sr-Nd-Pb同位素組成, 探討硫化物在巖漿起源和演化過程中的作用, 厘定再循環(huán)地殼組分的性質(zhì)和可能滯留時間, 在此基礎(chǔ)上探討海南島地幔柱演化的深部動力學(xué)過程。

    1 地質(zhì)背景與樣品描述

    海南島位于華南大陸的南端, 南海的北緣, 構(gòu)造位置上位于歐亞板塊、太平洋(菲律賓)板塊和印度?澳大利亞板塊交匯點的前緣, 具復(fù)雜的構(gòu)造背景和演化歷史(圖1a; Tu et al., 1992)。Ho et al. (2000)通過對雷瓊地區(qū)新生代玄武巖的年代學(xué)研究發(fā)現(xiàn), 該區(qū)巖漿活動始于上新世(3.8~5.6 Ma), 廣泛分布于早更新世(2.0~8.0 Ma), 并持續(xù)到晚中更新世(0.1~0.7 Ma)。Wang et al. (2012)對海南島玄武巖開展了Ar-Ar年代學(xué)研究, 獲得了海南島玄武巖噴發(fā)年齡為0.7~ 12.9 Ma。火山活動主要沿著深大斷裂分布, 早期的火山巖主要為石英拉斑玄武巖和橄欖石拉斑玄武巖; 后期則以堿性橄欖石玄武巖為主。

    本次研究對海南島新生代玄武巖開展系統(tǒng)的野外地質(zhì)調(diào)查, 并采集了25件巖石樣品, 具體采樣位置見圖1c。在野外, 可觀察到玄武巖柱狀節(jié)理(圖2a), 所采集樣品均較為新鮮, 呈灰黑色?黑色, 無明顯蝕變現(xiàn)象。巖石為致密塊狀構(gòu)造, 斑狀結(jié)構(gòu)(圖2b), 斑晶礦物主要為橄欖石(圖2b、c), 顯示出自形?半自形粒狀晶, 斑晶未見蝕變, 含量在5%~8%之間。其中標(biāo)號為YX樣品的橄欖石斑晶含量高達(dá)18%, 在橄欖石斑晶中可以觀察到尖晶石和硫化物包裹體(圖2c、d); 單斜輝石與斜長石斑晶少見; 基質(zhì)主要由斜長石、橄欖石、輝石微晶以及玻璃和不透明礦物組成。標(biāo)號HSL和GT的橄欖玄武巖樣品中可見地幔巖捕虜體, 捕虜體主要由單斜輝石、斜方輝石和橄欖石等礦物組成的尖晶石二輝橄欖巖,大者可達(dá)20 cm,小者約2~3 cm, 常呈橢圓狀或具菱角的四邊形。

    (a) 海南島玄武巖柱狀節(jié)理; (b) 橄欖玄武巖的手標(biāo)本照片; (c) 橄欖玄武巖中橄欖石為柱狀自形晶體, 包含有尖晶石; (d) 橄欖玄武巖中橄欖石內(nèi)的硫化物包裹體。礦物代號: Ol. 橄欖石; Sp. 尖晶石; Py. 黃鐵礦。

    2 分析方法

    2.1 全巖主量、微量元素及Sr-Nd-Pb同位素

    全巖的主量、微量元素分析測試在澳實分析檢測(廣州)有限公司完成。主量元素通過X射線熒光光譜儀(XRF)進(jìn)行測定, 分析誤差通常<2%。微量元素含量采用電感耦合等離子體質(zhì)譜(ICP-MS)分析獲得, 分析誤差在5%~10%之間(Zhou et al., 2014)。

    全巖Sr-Nd-Pb同位素含量在中國科學(xué)院廣州地球化學(xué)研究所同位素地球化學(xué)國家重點實驗室完成。稱取巖石粉末樣品100 mg, 并置于干凈的Teflon杯中。加入1 mL 1∶1高純濃HNO3以及2 mL高純濃HF, 輕輕搖勻并將其封蓋。隨后, 將它們置于120 ℃電熱板上保溫7天, 使樣品充分溶解, 然后按照實驗室詳細(xì)的操作流程及規(guī)則依次進(jìn)行全巖Sr-Nd-Pb同位素的分離工作。

    通過AG1X8陽離子交換樹脂以及稀釋的HBr將Pb從樣品溶液中分離出來, 化學(xué)流程中空白樣品Pb<50 pg, 使用國際標(biāo)準(zhǔn)溶液SRM981作為Pb同位素的監(jiān)測標(biāo)樣(Chen et al., 2002), 本次分析過程中的SRM 981分析結(jié)果為:206Pb/204Pb=16.9422±0.0006,207Pb/204Pb=15.4857±0.0006,208Pb/204Pb=36.6968±0.0015。通過陽離子交換樹脂將Sr與REE從溶液中分離, Sr同位素詳細(xì)分離流程見韋剛健等(2002)?;瘜W(xué)流程中空白樣品Sr<200 pg, 分析中使用Sr同位素國際標(biāo)準(zhǔn)樣品NBSSRM987(87Sr/86Sr=0.710238± 0.000010)進(jìn)行數(shù)據(jù)監(jiān)測, 采用86Sr/88Sr=0.1194進(jìn)行Sr同位素數(shù)據(jù)校正。Nd與其他稀土元素使用HDEHP柱進(jìn)行分離, 具體流程見梁細(xì)榮等(2003)?;瘜W(xué)流程中空白樣品Nd<30 pg, 分析過程中使用Nd同位素國際標(biāo)準(zhǔn)樣品為Shin Etsu JNdi-1(143Nd/144Nd=0.512114± 0.000006)進(jìn)行數(shù)據(jù)監(jiān)測, 采用146Nd/144Nd=0.7219進(jìn)行Nd同位素數(shù)據(jù)校正。

    2.2 Re-Os同位素

    為了避免樣品破碎過程中的污染, 用于Re-Os同位素和PGE分析的樣品, 在加工過程中用粗布或紙包裹樣品進(jìn)行破碎。中碎和細(xì)碎分別在剛玉和瑪瑙內(nèi)膽裝置中進(jìn)行。將約2 g粉末樣品與190Os和185Re稀釋劑一起溶解在王水中, 在Carius管中冷凍并密封, 然后在240 ℃溶解24小時(Shirey and Walker, 1998)。打開Carius管后, Os被提取出來, 用CCl4并反萃取到HBr中(Pearson and Woodland, 2000)。將剩余的含Re溶液干燥并重新溶解在HCl中, 在陰離子交換柱(AG1X8樹脂, 100~200目)中分離和純化Re。最后, 通過同位素稀釋?電感耦合等離子體質(zhì)譜法(ID-ICP-MS, Thermo-Scientific XSERIES-2)測量Re含量。進(jìn)一步蒸餾純化的Os溶液在Thermo-Finnigan Triton上測量OsO3和ReO4的同位素豐度質(zhì)譜儀處于負(fù)離子檢測模式(Creaser et al., 1991)。分別使用靜態(tài)多重法拉第收集器和脈沖計數(shù)電子倍增器測量Re和Os同位素組成。通過將測得的192Os/188Os比率歸一化為3.08271, 對Os的儀器質(zhì)量分餾進(jìn)行了校正。Re和Os均針對空白進(jìn)行了校正, Re和Os總空白分別為5.2±1.3和0.60±0.14 pg, 空白187Os/188Os值為0.26572±0.03895,因此空白對樣品Os和Re含量的貢獻(xiàn)很小(<1%)。測試詳細(xì)流程見Li et al. (2010, 2014)。

    2.3 PGE分析

    全巖PGE分析測試在中國科學(xué)院廣州地球化學(xué)研究所同位素地球化學(xué)國家重點實驗室利用硫鎳火試金?碲共沉淀方法完成。稱樣10~20 g, 加入硫磺、羥基鎳、硼酸鈉、碳酸鈉、硼砂、面粉均勻混合后置于坩堝中, 再加Os同位素稀釋劑。在1150 ℃條件下加熱熔融1.5小時后將熔體注入鐵模, 冷卻后取出鎳扣。將鎳扣粉碎后用濃鹽酸溶解, 加入碲共沉淀后, 將不溶物濾出, 加入王水在100 ℃加熱溶解1小時, 轉(zhuǎn)入比色管中定容, 然后利用ICP-MS (TJAPQ-EXCELL)進(jìn)行元素含量分析。分析測試過程中Os、Ir、Ru和Pd空白值小于0.03 pg, Pt和Rh空白值小于0.06 pg。詳細(xì)流程見孫亞莉等(1997)和李杰等(2011), 分析誤差小于10%。

    3 分析結(jié)果

    根據(jù)海南島玄武巖PGE和Re-Os同位素組成特征, 本次研究將海南島玄武巖分為高Os玄武巖(>40×10?12)和低Os玄武巖(<40×10?12)兩種類型進(jìn)行討論。

    3.1 主量、微量元素組成

    海南島新生代玄武巖全巖主量、微量元素分析結(jié)果見表1。本次研究的海南島新生代玄武巖具有低的LOI(<2.0%), MgO范圍為6.03%~14.70%, NaO+K2O為3.33%~5.12%, 在TAS圖解中(圖3)落在玄武巖、玄武安山巖和粗面玄武巖范圍內(nèi), 與前人的分析結(jié)果相似。

    在原始地幔標(biāo)準(zhǔn)化微量元素蛛網(wǎng)圖中, 海南島新生代玄武巖呈Nb-Ta正異常, HREE相對虧損, 總體上表現(xiàn)出與OIB相似特征(圖4a、c)。在球粒隕石標(biāo)準(zhǔn)化稀土元素配分模式圖中, 樣品呈右傾分布模式, 富集LREE, 虧損HREE, 具輕微Eu正異常(圖4b、d)。

    前人數(shù)據(jù)源: Wang et al., 2012; Liu et al., 2015; 梅盛旺和任鐘元, 2019。

    3.2 Sr-Nd-Pb-Os同位素組成

    由于海南島玄武巖時代年輕(多小于5.0 Ma), 因此所測得同位素數(shù)據(jù)在誤差范圍內(nèi)可認(rèn)為是其初始同位素值, 因此不再進(jìn)行同位素初始值校正。本次研究的海南島玄武巖Sr-Nd-Pb同位素分析結(jié)果見表2。玄武巖同位素組成變化較小,87Sr/86Sr=0.703220~ 0.704211,143Nd/144Nd=0.512803~0.512963,Nd=+3.22~ +6.34,206Pb/204Pb=18.5315~18.7571,207Pb/204Pb= 15.5166~ 15.6252,208Pb/204Pb=38.6847~38.9566, 與前人的結(jié)果基本一致。結(jié)合前人的相關(guān)數(shù)據(jù), 在Sr-Nd同位素圖解中, 海南島新生代玄武巖表現(xiàn)出DMM和EM2兩端元混合的特征(圖5a), Pb同位素組成位于北半球參考線(NHRL)之上(圖5b、c)。

    本次研究選取了Liu et al. (2015)、梅盛旺和任鐘元(2019)以及本次采集的高M(jìn)gO樣品進(jìn)行了Re-Os同位素分析, 結(jié)果見表3。高Os玄武巖的Re含量為36.1×10?12~131×10?12, Os含量為54.4×10?12~ 303×10?12,187Re/188Os值為0.85~11.5,187Os/188Os值為0.1251~0.1554, 變化相對較小, 虧損模式年齡MA介于32~347 Ma之間, 與Wang et al. (2013)報道的海南島Re-Os同位素數(shù)據(jù)相似。低Os玄武巖的Re含量為38.1×10?12~152×10?12, Os含量為2.9×10?12~ 38.4×10?12,187Re/188Os值為9.79~176.53之間,187Os/188Os值為0.1469~0.3018, 明顯高于軟流圈地幔值, 虧損模式年齡MA介于49~534 Ma之間。187Os/188Os-Os圖解顯示, 當(dāng)Os>40×10?12,187Os/188Os值變化很小; 而當(dāng)Os<40×10?12時, 隨著Os含量降低,187Os/188Os值顯著增加(圖5d); 且187Os/188Os<0.20低Os玄武巖樣品中, Os含量與187Os/188Os值之間顯示出一定的正相關(guān)關(guān)系。

    圖4 海南島玄武巖原始地幔標(biāo)準(zhǔn)化微量元素蛛網(wǎng)圖(a、c)和球粒隕石標(biāo)準(zhǔn)化稀土元素配分模式(b、d)(粒隕石和原始地幔數(shù)據(jù)來自Sun and McDonough, 1989)

    表2 海南島新生代玄武巖Sr-Nd-Pb同位素組成

    前人數(shù)據(jù)來源: Wang et al., 2012, 2013; Liu et al., 2015; 梅盛旺和任鐘元, 2019。

    為了理解巖石圈地幔對玄武巖成因的影響, 我們對海南島玄武巖攜帶的地幔橄欖巖包體(HN-1、2、3)進(jìn)行Re-Os同位素體系分析測試, 其結(jié)果見表3。這些尖晶石二輝橄欖巖的Re含量為41.8×10?12~ 207×10?12, Os含量為2390×10?12~5944×10?12,187Re/188Os值為0.08~0.17,187Os/188Os值為0.1186~0.1227, 明顯低于軟流圈地幔值(0.129; Shirey and Walker, 1998)和海南島玄武巖的187Os/188Os值。其虧損模式年齡MA介于1150~2443 Ma之間, 顯示出古老巖石圈地幔的同位素特點。

    3.3 PGE地球化學(xué)特征

    海南島新生代玄武巖全巖PGE分析結(jié)果見表4。原始地幔標(biāo)準(zhǔn)化鉑族元素圖(圖6)顯示, 玄武巖樣品PGE含量明顯低于地幔橄欖巖包體, 其ΣPGE含量為0.35×10?9~4.54×10?9, 富集PPGE, 分布模式總體呈向右上傾趨勢, 與地幔柱來源的夏威夷玄武巖特征相似(Tatsumi et al., 1999; Bennett et al., 2000; Jamais et al., 2008), 但是明顯不同于大洋中脊玄武巖(MORB)強(qiáng)烈虧損Ir的特征(Barnes et al., 2015)。根據(jù)其配分模式, 把海南玄武巖分為兩類: 高Os玄武巖配分曲線近水平, Pd/Ir值較低, 總體小于3.5; 低Os玄武巖配分曲線為向右陡上傾, Pd/Ir值較高, 總體大于5.0。相對于低Os玄武巖, 高Os玄武巖具有較低Pd/Ir值和較高La/Sm值(圖7a)。在187Os/188Os-Pd/Ir圖解中, 高Os玄武巖Pd/Ir值與Os同位素組成之間相關(guān)關(guān)系不明顯, 而低Os玄武巖中則存在明顯的正相關(guān)關(guān)系(圖7b)。與兩類玄武巖明顯不同, 地幔包體的PGE配分曲線向右下傾, Pd/Ir值變化在0.4~0.8之間, 表現(xiàn)出巖漿抽提之后殘留地幔的特點。

    表3 海南島新生代玄武巖Re-Os同位素組成

    表4 鉑族元素分析結(jié)果(×10?9)

    夏威夷玄武巖的PGE數(shù)據(jù)來源: Tatsumi et al., 1999; Bennett et al., 2000; Jamais et al., 2008。原始地幔值來自Becker et al., 2006。

    4 討 論

    海南島新生代玄武巖在噴發(fā)地表過程中要經(jīng)過大陸巖石圈地幔和地殼, 因此在討論其成因和地幔源區(qū)之前很有必要討論其低溫巖漿過程, 如低溫蝕變, 地殼混染, 分離結(jié)晶等過程的影響。

    4.1 低溫蝕變與地殼混染

    海南島玄武巖樣品新鮮, 斑晶礦物表面蝕變較弱, 總體上具有低燒失量, 強(qiáng)不相容元素Ba、Sr、La、U和Nb(蝕變過程中穩(wěn)定元素)呈明顯的正相關(guān)關(guān)系, 因此可以排除蝕變作用對全巖元素地球化學(xué)特征的影響。前人研究表明海南島地區(qū)新生代玄武巖受到地殼混染的影響可忽略不計(Wang et al., 2013; Liu et al., 2015; 梅盛旺和任鐘元, 2019)。本次研究的玄武巖樣品具有較高M(jìn)gO(>8%)含量, 微量元素蛛網(wǎng)圖中具有正Nb-Ta異常, 與大陸地殼強(qiáng)烈虧損Nb-Ta的特征截然相反。Ce/Pb和Nb/U值分別為16.3~26.3與34.7~52.7, 與全球的MORB和OIB的一致(Ce/Pb=25±5, Nb/U=47±7; Hofmann et al., 1986),也說明地殼混染影響很小。玄武巖全巖Nd和Sr同位素組成較為均一, 隨著MgO含量降低并沒有表現(xiàn)出類似于地殼混染的趨勢, 因此其受到地殼混染的影響可以忽略。

    通常認(rèn)為, 玄武巖中Os>40×10?12的樣品遭受到地殼混染的概率較小(Molzahn et al., 1996; Widom, 1997)。在187Os/188Os-Os圖解中(圖5d), Os>40×10?12的樣品, 其187Os/188Os值變化較小(0.1251~0.1554), Os同位素基本沒有受到地殼混染的作用, 總體反映了源區(qū)的Os同位素組成。Os<40×10?12的樣品, 隨著Os含量降低,187Os/188Os值快速上升, 顯示出類似地殼混染作用的特征。但是187Os/188Os<0.20時, 低Os玄武巖樣品顯示出Os含量與187Os/188Os值的正相關(guān)關(guān)系, 顯然這種變化趨勢不適合采用地殼同化混染作用來解釋。對于187Os/188Os>0.20低Os玄武巖樣品, 隨著Os含量降低,187Os/188Os值大幅增加, 顯示為地殼混染/同化混染?結(jié)晶分異過程(AFC)的結(jié)果。這與鏡下觀察到橄欖石斑晶中的硫化物包裹體相一致。AFC模擬結(jié)果顯示, 地殼物質(zhì)混染速率()變化在0.3~0.7之間時, 殘余巖漿比例為95%~97%, 可以較好地擬合部分低Os玄武巖的Os同位素變化(圖5d)。這么小比例的地殼同化混染作用對于全巖的主量、微量元素和Sr-Nd-Pb同位素的影響較弱。

    4.2 部分熔融和分離結(jié)晶作用

    海南島地區(qū)新生代玄武巖的巖相學(xué)特征和巖石地球化學(xué)特征均反映這些巖石經(jīng)歷了不同程度的分離結(jié)晶作用。其中MgO和Ni、Co之間呈正相關(guān)關(guān)系, 指示了橄欖石的分離結(jié)晶作用; 而MgO和CaO/Al2O3之間無明顯的相關(guān)關(guān)系說明單斜輝石不是最主要的分離結(jié)晶相, 這與鏡下觀察到僅有少量的單斜輝石斑晶一致。本次研究中斜長石的分離結(jié)晶作用也很少, 在稀土元素配分模式上甚至出現(xiàn)了Eu的正異常, 這與鏡下所觀察到的橄欖石斑晶的廣泛分布, 少見有單斜輝石和斜長石斑晶的巖相學(xué)特征一致。

    部分熔融程度可以導(dǎo)致PGE內(nèi)部發(fā)生分餾作用(Barnes et al., 1985, 2015; Barnes, 1993; Shirey and Walker, 1998; Momme et al., 2003)。高部分熔融程度(>20%)的熔體, 例如科馬提巖, 通常為S不飽和巖漿, 熔融過程中, 源區(qū)的PGE能夠完全進(jìn)入巖漿中, 這種熔體往往具有相對平坦的球粒隕石標(biāo)準(zhǔn)化PGE 模式(Pd/Ir值接近球粒隕石)和較高的PGE含量(Puchtel and Humayun, 2000; Puchtel et al., 2004)。相反, 較低程度的部分熔融熔體(<20%), 例如大多數(shù)玄武巖, 往往是S飽和的巖漿(Barnes et al., 2015), 導(dǎo)致地幔中硫化物的不完全溶解和部分殘留在源區(qū)。由于IPGE在單體硫化物(MSS)中是相容元素, 如果熔融過程中MSS進(jìn)入熔體, 那么巖漿就會富集IPGE; 相反地, 如果部分MSS殘留在源區(qū), 巖漿則顯示出IPGE的虧損, 如MORB具有非常低的Ir和Ru含量, 以及大部分OIB都顯示出IPGE的虧損(Barnes et al., 2015)。粒間富Cu硫化物(ISS)則不會產(chǎn)生IPGE與PPGE之間的分餾作用。地幔熔融過程中, ISS通常在MSS之前發(fā)生熔融, 由于源區(qū)存在MSS殘留, 導(dǎo)致衍生巖漿中PPGE與IPGE的分餾作用。海南島低Os玄武巖具有較低IPGE含量和較高的Pd/Ir值, PPGE與IPGE分餾作用較為明顯。如果低Os玄武巖的高Pd/Ir值反映了較低的地幔部分熔融程度, 其應(yīng)該相對于高Os玄武巖具有更高的La/Sm值(La/Sm值與地幔熔融程度具有負(fù)相關(guān)關(guān)系)。然而, 在La/Sm-Pd/Ir圖解中(圖7a), 并沒有觀察到低程度部分熔融對應(yīng)高的La/Sm值, 與此相反, 低Os玄武巖相對高Os玄武巖具有更低的La/Sm值, 這說明部分熔融程度不是導(dǎo)致兩類玄武巖具有差異性PGE配分模式的主要原因。因此兩類玄武巖不同的PGE配分模式可能與源區(qū)不均一性有關(guān)。

    圖7 海南島玄武巖的La/Sm-Pd/Ir(a)和187Os/188Os-Pd/Ir(b)圖解

    由于PGE在硫化物中具有極高的分配系數(shù), 故巖漿中微量的硫化物熔離作用或者源區(qū)中極少量硫化物的殘留都會造成巖漿中PGE含量的明顯降低。在S不飽和巖漿中, Cu 和Pd 均為不相容元素, 因此原始巖漿的Cu/Pd值應(yīng)接近地幔值(Cu/Pd≈10000)。如果部分硫化物殘留于地幔, 由于PGE在硫化物中的分配系數(shù)遠(yuǎn)高于Cu, 故Cu/Pd值將高于地幔值(Fleet et al., 1996)。本區(qū)所有樣品的Cu/Pd值均高于地幔值, 且其PGE含量相對原始地幔明顯虧損, 說明巖漿上升過程中可能存在硫化物的熔離作用。實驗巖石學(xué)已證明PGE元素在硫化物和硅酸鹽中的分配系數(shù)差異很大(Bezmen et al., 1994), 因此海南島玄武巖中低Os玄武巖的PGE陡向上右傾斜的原始地幔標(biāo)準(zhǔn)化配分模式不能用硫化物的熔離作用來解釋。前人多認(rèn)為這種差異反映巖漿演化歷史中, 分離礦物相的差異。比如, 鉻鐵礦的分餾可能導(dǎo)致IPGE的分餾, 因為鉻鐵礦可形成包裹IPGE的合金(Ballhaus et al., 2006)。不混溶硫化物中, 單硫化物的熔離也會導(dǎo)致IPGE與PPGE的分餾。因此我們推測低Os玄武巖的PGE配分模式是由于巖漿演化過程中富IPGE相的熔離(單硫化物或是賦存在鉻鐵礦中的IPGE合金)。而IPGE合金相在地幔源區(qū)的殘留可能是形成高Os玄武巖PGE配分模式的主要原因。

    為了定量的評估硫化物對PGE的影響, 我們利用Cu/Pd-Pd圖解來模擬S飽和巖漿的熔融和分異過程。在該模型中, 用MORB地幔(DMM)的PGE含量作為初始值, 計算不同比例硫化物熔離作用或者殘留源區(qū)時Pd和Cu/Pd值的變化情況。模擬結(jié)果顯示,高Os玄武巖的Pd與Cu/Pd值的變化關(guān)系可以通過0.002%~0.001%的硫化物殘留在源區(qū)的熔融過程進(jìn)行解釋; 低Os玄武巖極高Cu/Pd值則很可能是巖漿上升過程中硫化物分離結(jié)晶作用的結(jié)果(圖8)。

    4.3 海南島新生代玄武巖源區(qū)性質(zhì)

    4.3.1 巖石圈地幔或軟流圈地幔?

    海南島新生代玄武巖同位素特征反映該區(qū)玄武巖源區(qū)至少包含了虧損與富集兩種端元組分(圖5), 韓江偉等(2009)把海南島玄武巖富集不相容元素和同位素的地球化學(xué)特征歸結(jié)于陸下巖石圈的貢獻(xiàn)。范蔚茗等(1993)對該區(qū)地幔捕虜體的研究結(jié)果顯示, 其具有比寄主玄武巖具有更虧損的Sr-Nd同位素組成, 顯然說明了該區(qū)巖石圈的部分熔融不能產(chǎn)生富集的同位素特征, 也不能作為EM2的源區(qū)。

    黑色曲線代表含殘余硫化物的橄欖巖的簡單分批熔化曲線。橄欖石、斜輝石、斜輝石和石榴石的 Pd 分配系數(shù)來自Chazey and Neal (2005)。橄欖石、斜輝石、斜輝石和石榴石對 Cu 的分配系數(shù)分別來自Yurimoto and Ohtani (1992); Hart and Dunn (1993); Klemme et al. (2006)。Pd和Cu在硫化物中的分配系數(shù)來自Mungall and Brenan (2014)。橄欖巖熔融過程中殘余相的比例假定為62%的橄欖石, 14.9%~14.9999%的單斜輝石, 15%的斜方輝石, 8%的石榴子石以及0.1%~0.001%的硫化物。DMM的數(shù)據(jù)來自Salters and Stracke (2004)。

    本次研究的三個地幔包體Os同位素組成較為虧損,187Os/188Os值在0.1186~0.1227之間(表3), 如果此類物質(zhì)的加入將會造成玄武巖Os含量的增加以及Os同位素比值的減少, 但在Os-187Os/188Os圖中并未見到這樣的趨勢, 同樣說明了大陸巖石圈地幔對其源區(qū)的影響很小。由于地幔包體相對于海南島高Os與低Os兩類玄武巖具有更高的Nd同位素組成(范蔚茗等, 1993), 巖石圈地幔的加入會同時導(dǎo)致187Os/188Os值降低、PGE含量和143Nd/144Nd值的增加, 然而高Os玄武巖盡管相對低Os玄武巖具有低187Os/188Os值和高PGE含量, 但是其143Nd/144Nd值總體較低Os玄武巖略低(表2), 顯然也不支持大陸巖石圈地幔的貢獻(xiàn)。因此, 高Os玄武巖相對平坦的分布模式(Pd/Ir值小)和變化較小的Os同位素組成反映其主要來源于軟流圈地幔。

    近年來海南島地區(qū)開展詳細(xì)的地球物理和巖石地球化學(xué)研究, 普遍認(rèn)為海南島新生代玄武巖為地幔柱作用的結(jié)果(Huang and Zhao, 2006; Zhao, 2007; Wang et al., 2012, 2013; Liu et al., 2015; Xia et al., 2016)。海南島玄武巖表現(xiàn)出類似夏威夷地幔柱起源玄武巖的PGE配分模式(圖6), 也支持其地幔柱成因。

    4.3.2 再循環(huán)洋殼的貢獻(xiàn)

    微量元素蛛網(wǎng)圖中, 研究區(qū)所有玄武巖樣品均顯示Nb、Ta富集和Th、U、LREE虧損(圖4c、d)。Nb、Ta是流體不流動元素, 在俯沖帶脫水過程中保留在洋殼中, 而 K、Li、Pb 和 LILE 是流體流動元素, 俯沖過程中轉(zhuǎn)移到上覆地幔楔(Elliott et al., 1997; Kessel et al., 2005)。因此, 大洋板塊經(jīng)過俯沖脫水變成榴輝巖, Pb隨著流體被遷出, 而Nb和Ta則殘留在板片中, 這些富含Nb、Ta的板片在下沉過程發(fā)生部分熔融, 形成的熔體交代周圍地幔, 形成Nb-Ta正異常的OIB型地幔源區(qū)(Wu et al., 2020; Guo et al., 2021)。這些特征與海南島玄武巖的OIB型微量元素地球化學(xué)特征相吻合。此外, 巖石Sr和Eu正異常表明了源區(qū)中富含斜長石。一般來說, 地幔中有兩個潛在的富含斜長石的區(qū)域: ①輝長巖下洋殼(現(xiàn)在榴輝巖或輝石巖; Hofmann and Jochum, 1996); ②大陸下地殼(Rudnick and Fountain, 1995), 例如鎂鐵質(zhì)麻粒巖。但是大陸下地殼的參與, 會導(dǎo)致Pb同位素比值降低, 與該區(qū)的Pb同位素組成不符。因此, 我們推斷, 作為俯沖地殼的一部分進(jìn)入源區(qū)的再循環(huán)大洋輝長巖的存在, 可以較好地解釋Sr和Eu正異常特征(Wang et al., 2012; Liu et al., 2015; 梅盛旺和任鐘元, 2019)。

    前人對夏威夷等地玄武巖Os同位素研究揭示, 放射性 Os 和Pb同位素組成為地幔中循環(huán)洋殼的存在提供了強(qiáng)有力的支持(Escrig et al., 2005; Dale et al., 2007; Day et al., 2009)。研究區(qū)沒有遭受地殼混染的玄武巖樣品, 其187Os/188Os值為0.1225~0.1554, 明顯高于原始地幔的值說明了源區(qū)中有循環(huán)洋殼的參與。

    4.3.3 再循環(huán)洋殼年齡

    關(guān)于海南島地幔柱再循環(huán)洋殼的年齡, Wang et al. (2013)通過Os-Pb同位素模擬結(jié)果認(rèn)為地幔熔融源區(qū)中存在年輕洋殼(200~500 Ma)組分; 梅盛旺和任仲元(2019)通過對海南島新生代玄武巖的Pb-Pb同位素組成估算再循環(huán)洋殼年齡在1000 Ma以內(nèi), 并以此解釋為何海南島玄武巖沒有表現(xiàn)出洋島玄武巖的HIMU特征。Li et al. (2020)通過對海南島及東南沿海新生代玄武巖的Nd-Pb同位素組成進(jìn)行計算, 認(rèn)為海南島玄武巖的組分中包括了古老的再循環(huán)洋殼和太平洋地殼物質(zhì)以及滯留的俯沖板片, 而滯留的年輕板片來自太平洋板塊。

    研究區(qū)玄武巖Re-Os同位素數(shù)據(jù)顯示, Os模式年齡(MA)為32~534 Ma, 明顯低于上覆的大陸巖石圈地幔, 暗示其源巖為較為年輕的地幔源區(qū)。為了進(jìn)一步限定再循環(huán)洋殼的年齡, 我們根據(jù)海南島全巖的Pb和Os同位素組成的變化特征進(jìn)行了模擬計算。通過設(shè)置初始地幔、再循環(huán)洋殼/輝長巖洋殼的Pb含量、Pb同位素比值和Os含量以及Os同位素比值, 可以大體限制再循環(huán)洋殼的年齡。研究結(jié)果顯示再循環(huán)洋殼的年齡在100~300 Ma之間(圖9), 主要來源于周緣俯沖再循環(huán)洋殼。

    通過對海南島新生代玄武巖以及巖石圈地幔包體的PGE和Re-Os同位素的綜合研究, 發(fā)現(xiàn)其地幔源區(qū)中缺乏古老巖石圈地幔的貢獻(xiàn)。新生代, 海南島經(jīng)歷了多期次的大洋板塊俯沖作用, 包括太平洋、菲律賓、澳大利亞和古南海等板塊, 地幔柱物質(zhì)在上升過程中被軟流圈強(qiáng)烈的對流作用驅(qū)散(Li et al., 2020), 或者被來自周圍的俯沖板片壓制而難以影響到巖石圈地幔(Wang et al., 2013)。另外, 海南島玄武巖中水含量研究顯示海南島地幔柱相對貧水(Gu et al., 2019), 因此其上升過程中源區(qū)熔融程度相對較低, 黏度較大, 對上覆巖石圈的影響較弱。Li et al. (2020)對海南帶及東南沿海地區(qū)新生代玄武巖的綜合研究也顯示海南島地幔柱的影響可能覆蓋了整個中國東南部, 但是地幔熱柱作用的表現(xiàn)為沿主要的巖石圈斷裂帶發(fā)育了玄武巖, 認(rèn)為海南島地幔柱可能類似古生代的塔里木地幔柱(徐義剛等, 2013), 地幔柱的能量較低, 很可能屬于潛伏模式。

    前人數(shù)據(jù)引自Wang et al. (2013); 計算方法據(jù)Wu et al. (2020); 循環(huán)洋殼和輝長巖的Os含量及組成來自Peucker-Ehrenbrink et al. (2012); 用Stracke et al. (2003)的方法計算不同循環(huán)年齡的Pb同位素組成。

    5 結(jié) 論

    (1) 海南島新生代玄武巖總體上具有相對原始地幔低的PGE含量, 可以劃分為兩種類型: PGE相對平坦型(低Pd/Ir值)和向右陡上傾型(高Pd/Ir值), 前者對應(yīng)了較高Os含量和低187Os/188Os同位素比值, 代表了主要來自軟流圈地幔的熔體; 后者對應(yīng)了低Os含量和高187Os/188Os同位素比值, 反映了源區(qū)中含有大量的再循環(huán)洋殼組分。根據(jù)Pb-Os同位素協(xié)變關(guān)系, 估計再循環(huán)洋殼的年齡在100~300 Ma之間, 主要來源于周緣俯沖再循環(huán)洋殼。

    (2) 海南島玄武巖的PGE特征主要反映了巖漿在上升過程中發(fā)生了硫化物熔離作用, 熔融過程中可能伴隨有少量硫化物殘留源區(qū)。

    (3) 海南島新生代玄武巖源區(qū)中缺乏上覆古老巖石圈地幔的貢獻(xiàn), 暗示地幔柱物質(zhì)在上升過程中被對流軟流圈或俯沖板片壓制, 從而導(dǎo)致潛伏的地幔柱側(cè)向流動到周緣的中南半島、南海和東南沿海地區(qū), 形成彌散型的大火成巖省。

    致謝:研究過程中得到了中國科學(xué)院廣州地球化學(xué)研究所任鐘元研究員的指導(dǎo), 野外工作得到了中國科學(xué)院廣州地球化學(xué)研究所宋茂雙研究員和胡琪偉同學(xué)的大力協(xié)助, 中國科學(xué)院廣州地球化學(xué)研究所任鐘元研究員和云南大學(xué)王選策教授為本文提出中肯評審意見和建議, 在此一并致以誠摯的謝意。

    范蔚茗, Menzies M A, 尹漢輝, 郝新才, 陳欣, 鄒和平. 1993. 中國東南沿海深部巖石圈的性質(zhì)和深部作用過程初探. 大地構(gòu)造與成礦學(xué), 17(1): 23–30.

    韓江偉, 熊小林, 朱照宇. 2009. 雷瓊地區(qū)晚新生代玄武巖地球化學(xué): EM2成分來源及大陸巖石圈地幔的貢獻(xiàn). 巖石學(xué)報, 25(12): 3208–3220.

    李杰, 鐘立峰, 涂湘林, 胡光黔, 孫彥敏, 梁細(xì)榮, 許繼峰. 2011. 利用同一化學(xué)流程分析地質(zhì)樣品中的鉑族元素含量和錸?鋨同位素組成. 地球化學(xué), 40(4): 372– 380.

    梁細(xì)榮, 韋剛健, 李獻(xiàn)華, 劉穎. 2003. 利用MC-ICPMS精確測定143Nd/144Nd和Sm/Nd比值. 地球化學(xué), 32(1): 91–96.

    梅盛旺, 任鐘元. 2019. 海南島新生代玄武質(zhì)熔巖源區(qū)中再循環(huán)物質(zhì)及其年齡的限定: 來自Hf-Sr-Nd-Pb同位素的制約. 大地構(gòu)造與成礦學(xué), 43(5): 1036–1051.

    孫亞莉, 管希云, 杜安道. 1997. 锍試金富集貴金屬元素: Ⅰ等離子體質(zhì)譜法測定地質(zhì)樣品中痕量鉑族元素. 巖礦測試, 16(1): 12–17.

    韋剛健, 梁細(xì)榮, 李獻(xiàn)華, 劉穎. 2002. MC-ICPMS方法精確測定液體和固體樣品的Sr同位素組成. 地球化學(xué), 31(3): 295–299.

    徐義剛, 何斌, 羅震宇, 劉海泉. 2013. 我國大火成巖省和地幔柱研究進(jìn)展與展望. 礦物巖石地球化學(xué)通報, 32(1): 25–39.

    Ballhaus C, Bockrath C, Wohlgemuth-Ueberwasser C, Laurenz V, Berndt J. 2006. Fractionation of the noble metals by physical processes., 152(6): 667–684.

    Barnes S J. 1993. Partitioning of the platinum group elements and gold between silicate and sulphide magmas in the Munni Munni Complex, Western Australia., 57(6): 1277–1290.

    Barnes S J, Mungall J E, Maier W D. 2015. Platinum group elements in mantle melts and mantle samples., 232: 395–417.

    Barnes S J, Naldrett A J, Gorton M P. 1985. The origin of the fractionation of platinum-group elements in terrestrial magmas., 53(3–4): 303–323.

    Becker H, Horan M F, Walker R J, Gao S, Lorand J P, Rudnick R L. 2006. Highly siderophile element composition of the Earth’s primitive upper mantle: Constraints from new data on peridotite massifs and xenoliths., 70(17): 4528–4550.

    Bennett V C, Norman M D, Garcia M O. 2000. Rhenium and platinum group element abundances correlated with mantlesource components in Hawaiian picrites: Sulphides in the plume., 183(3–4): 513–521.

    Bezmen N I, Asif M, Brugmann G E, Romanenko I M, Naldrett A J. 1994. Distribution of Pd, Rh, Ru, Ir, Os, and Au between sulfide and silicate metals., 58(4): 1251–1260.

    Chazey W J, Neal C R. 2005. Platinum-group element constraints on source composition and magma evolution of the Kerguelen Plateau using basalts from ODP Leg 183., 69(19): 4685–4701.

    Chen F K, Satir M, Ji J Q, Zhang D. 2002. Nd-Sr-Pb isotopes of Tengchong Cenozoic volcanic rocks from western Yunnan, China: Evidence for an enriched mantle source., 21(1): 39–45.

    Creaser R A, Papanastassiou D A, Wasserburg G J. 1991. Negative thermal ion mass spectrometry of osmium, rhenium and iridium., 55(1): 397–401.

    Dale C W, Gannoun A, Burton K W, Argles T W, Parkinson I J. 2007. Rhenium-osmium isotope and elemental behaviour during subduction of oceanic crust and the implications for mantle recycling., 253(1–2): 211–225.

    Day J M D, Pearson D G, Macpherson C G, Lowry D, Carracedo J C. 2009. Pyroxenite-rich mantle formed by recycled oceanic lithosphere: Oxygen-osmium isotope evidence from Canary Island lavas., 37(6): 555–558.

    Elliott T, Plank T, Zindler A, White W, Bourdon B. 1997. Element transport from slab to volcanic front at the Mariana arc.:, 102(B7): 14991–15019.

    Escrig S, Doucelance R, Moreira M, Allègre C J. 2005. Os isotope systematics in Fogo Island: Evidence for lower continental crust fragments under the Cape Verde Southern Islands., 219(1–4): 93–113.

    Fleet M E, Crocket J H, Stone W E. 1996. Partitioning of platinum-group elements (Os, Ir, Ru, Pt, Pd) and gold between sulfide liquid and basalt melt., 60(13): 2397–2412.

    Gao S, Rudnick R L, Carlson R W, McDonough W F, Liu Y S. 2002. Re-Os evidence for replacement of ancient mantle lithosphere beneath the North China craton., 198(3–4): 307–322.

    Gu X Y, Wang P Y, Kuritani T, Hanski E, Xia Q K, Wang Q Y. 2019. Low water content in the mantle source of the Hainan plume as a factor inhibiting the formation of a large igneous province., 515: 221–230.

    Guo F, Wu Y M, Zhang B, Zhang X B, Zhao L, Liao J. 2021. Magmatic responses to Cretaceous subduction and tearing of the paleo-Pacific Plate in SE China: An overview., 212(1), 103448.

    Hart S R, Dunn T. 1993. Experimental cpx/melt partitioning of 24 trace elements., 113(1): 1–8.

    Ho K S, Chen J C, Juang W S. 2000. Geochronology and geochemistry of late Cenozoic basalts from the Leiqiong area, southern China., 18(3): 307–324.

    Hoang N, Flower M F J. 1998. Petrogenesis of Cenozoic basalts from Vietnam: Implication for origins of a ‘diffuseigneous province’., 39(3): 369–395.

    Hoang N, Flower M F J, Carlson R W. 1996. Major, trace element, and isotopic compositions of Vietnamese basalts: Interaction of hydrous EM1-rich asthenosphere with thinned Eurasian lithosphere., 60(22): 4329–4351.

    Hofmann A W, Jochum K P. 1996. Source characteristics derived from very incompatible trace elements in Mauna Loa and Mauna Kea basalts, Hawaii scientific drilling project.:, 101(B5): 11831–11839.

    Hofmann A W, Jochum K P, Seufert M, White W M. 1986. Nb and Pb in oceanic basalts: New constraints on mantleevolution., 79(1–2): 33–45.

    Huang J L, Zhao D P. 2006. High-resolution mantle tomography of China and surrounding regions.:, 111, B09305.

    Jamais M, Lassiter J C, Brügmann G. 2008. PGE and Os-isotopic variations in lavas from Kohala Volcano, Hawaii: Constraints on PGE behavior and melt/crust interaction., 250(1–4): 16–28.

    Kessel R, Schmidt M W, Ulmer P, Pettke T. 2005. Trace element signature of subduction-zone fluids, melts and supercritical liquids at 120–180?km depth., 437(7059): 724–727.

    Klemme S, Günther D, Hametner K, Prowatke S, Zack T. 2006. The partitioning of trace elements between ilmenite, ulvospinel, armalcolite and silicate melts with implicationsfor the early differentiation of the moon., 234(3–4): 251–263.

    Li J, Jiang X Y, Xu J F, Zhong L F, Wang X C, Wang G Q, Zhao P P. 2014. Determination of platinum-group elements and Re-Os isotopes using ID-ICP-MS and N-TIMS from a single digestion after two-stage column separation., 38(1): 37–50.

    Li J, Xu J F, Suzuki K, He B, Xu Y G, Ren Z Y. 2010. Os, Nd and Sr isotope and trace element geochemistry of the Mulipicrites: Insights into the mantle source of the Emeishan Large Igneous Province., 119(1–2): 108–122.

    Li Y Q, Kitagawa H, Nakamura E, Ma C Q, Hu X Y, Kobayashi K, Sakaguchi C. 2020. Various ages of recycled material in the source of Cenozoic Basalts in SE China: Implications for the role of the Hainan Plume., 61(6), egaa060.

    Liu J Q, Ren Z Y, Nichols A R L, Song M S, Qian S P, Zhang Y, Zhao P P. 2015. Petrogenesis of Late Cenozoic basalts from North Hainan Island: Constraints from melt inclusions and their host olivines., 152: 89–121.

    Lorand J P, Luguet A, Alard O. 2013. Platinum-group element systematics and petrogenetic processing of the continental upper mantle: A review., 164–167: 2–21.

    Molzahn M, Reisberg L, W?rner G. 1996. Os, Sr, Nd, Pb, O isotope and trace element data from the Ferrar flood basalts, Antarctica: Evidence for an enriched subcontinentallithospheric source., 144(3): 529–545.

    Momme P, óskarsson N, Keays R R. 2003. Platinum-group elements in the Icelandic rift system: Melting processes and mantle sources beneath Iceland., 196(1–4): 209–234.

    Mungall J E, Brenan J M. 2014. Partitioning of platinum- group elements and Au between sulfide liquid and basalt and the origins of mantle-crust fractionation of the chalcophile elements., 125: 265–289.

    Pearson D G, Irvine G J, Ionov D A, Boyd F R, Dreibus G E. 2004. Re-Os isotope systematics and platinum group element fractionation during mantle melt extraction: A study of massif and xenolith peridotite suites., 208(1–4): 29–59.

    Pearson D G, Woodland S J. 2000. Solvent extraction/anion exchange separation and determination of PGEs Os, Ir, Pt, Pd, Ru and Re-Os isotopes in geological samples by isotope dilution ICP-MS., 165(1–2): 87–107.

    Peucker-Ehrenbrink B, Hanghoj K, Atwood T, Kelemen P B. 2012. Rhenium-osmium isotope systematics and platinum group element concentrations in oceanic crust., 40(3): 199–202.

    Puchtel I S, Humayun M. 2000. Platinum group elements in Kostomuksha komatiites and basalts: Implications for oceanic crust recycling and core-mantle interaction., 64(24): 4227–4242.

    Puchtel I S, Humayun M, Campbell A J, Sproule R A, Lesher C M. 2004. Platinum group element geochemistry of komatiites from the Alexo and Pyke Hill areas, Ontario, Canada., 68(6): 1361–1383.

    Rudnick R L, Fountain D M. 1995. Nature and composition of the continental crust: A lower crustal perspective., 33(3): 267–309.

    Salters V J M, Stracke A. 2004. Composition of the depleted mantle.,,, 5, Q05B07.

    Shirey S B, Walker R J. 1998. The Re-Os isotope system in cosmochemistry and high-temperature geochemistry., 26: 423–500.

    Stracke A, Bizimis M, Salters V J M. 2003. Recycling oceaniccrust: Quantitative constraints., 4(3), 8003.

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

    Tatsumi Y, Oguri K, Shimoda G. 1999. The behaviour of platinum-group elements during magmatic differentiation in Hawaiian tholeiites., 33(4): 237–247.

    Tu K, Flower M F J, Carlson R W, Xie G H, Chen C Y, Zhang M. 1992. Magmatism in the South China Basin: 1. Isotopic and trace-element evidence for an endogenous Dupal mantle component., 97(1–2): 47–63.

    Wang X C, Li Z X, Li X H, Li J, Liu Y, Long W G, Zhou J B, Wang F. 2012. Temperature, pressure, and composition of the mantle source region of Late Cenozoic basalts in Hainan Island, SE Asia: A consequence of a young thermal mantle plume close to subduction zones?, 53(1): 177–233.

    Wang X C, Li Z X, Li X H, Li J, Xu Y G, Li W X. 2013. Identification of an ancient mantle reservoir and young recycled materials in the source region of a young mantle plume: Implications for potential linkages between plume and plate tectonics., 377–378: 248–259.

    Widom E. 1997. Sources of ocean island basalts: A review of theosmium isotope evidence.:, 244(1): 484–496.

    Wu Y M, Guo F, Wang X C, Zhang B, Zhang X B, Alemayehu M, Wang G Q. 2020. Generation of late Cretaceous Ji’an basalts through asthenosphere-slab interaction in South China., 132(5–6): 1316–1332.

    Xia S H, Zhao D P, Sun J L, Huang H B. 2016. Teleseismic imaging of the mantle beneath southernmost China: New insights into the Hainan plume., 36: 46–56.

    Yan Q S, Shi X F, Metcalfe I, Liu S F, Xu T Y, Kornkanitnan N, Sirichaiseth T, Yuan L, Zhang Y, Zhang H. 2018. Hainan mantle plume produced late Cenozoic basaltic rocks in Thailand, Southeast Asia., 8(1): 2640.

    Yurimoto H, Ohtani E. 1992. Element partitioning between majorite and liquid: A secondary ion mass spectrometric study., 19(1): 17–20.

    Zhao D P. 2007. Seismic images under 60 hotspots: Search for mantle plumes., 12(4): 335–355.

    Zhou L G, Xia X Q, Zheng Y F, Hu Z C. 2014. Polyphase growth of garnet in eclogite from the Hong’an orogen: Constraints from garnet zoning and phase equilibrium., 206–207: 79–99.

    Zou H B, Fan Q C. 2010. U-Th isotopes in Hainan basalts: Implications for sub-asthenospheric origin of EM2 mantle endmember and the dynamics of melting beneath Hainan Island., 116(1–2): 145–152.

    Origin of Cenozoic Basalts in Hainan Island: Constraints from Re-Os Isotops and PGE Geochemistry

    WANG Guoqing1, 2, GUO Feng1, 3*, ZHAO Liang1, 3

    (1. State Key Laboratory of Isotope Geochemistry, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences,Guangzhou 510640, Guangdong, China; 2. College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China; 3. CAS Center for Excellence in Deep Earth Science, Guangzhou 510640, Guangdong, China)

    In this paper, we report PGE and Re-Os isotope compositions of the Cenozoic basalts in Hainan Island. In combination with the whole-rock major, trace element and Sr-Nd-Pb isotopic results, the following conclusions can be reached. (1) The Hainan basalts can be subdivided into two groups on the basis of chondrite-normalized PGE patterns. The high-Os (Os>40×10?12) basalts have PGE concentrations varying from 1.88×10?9to 4.54×10?9and show a relatively flat PGE pattern with Pd/Ir<3.5; whereas the low-Os (Os<40×10?12) basalts have PGE contents varying from 0.35×10?9to 2.67×10?9and demonstrate an upright-inclined PGE patterns with Pd/Ir>5.0. The high-Os basalts also have low187Os/188Os ratios (0.1251–0.1554), representative of melt from the depleted asthenospheric mantle. By contrast, the low-Os basalts have higher187Os/188Os ratios (0.1469–0.3018), suggestive of significant involvement of recycled crust either in the mantle source or during the magmatic differentiation. The samples with187Os/188Os<0.20 show a positive correlation between Os and187Os/188Os, which cannot be produced by crustal assimilation. The radiogenic Os compositions of the low-Os basalts might have primarily been derived from the recycled oceanic crust in the melting source. (2) The PGE geochemistry of the Cenozoic Hainan basalts might have been controlled by sulfide segregation during the magma differentiation. (3) The Os-Pb isotope modelling results imply that the age of recycled oceanic crust for the Cenozoic Hainan basalts varies between 100 and 300 Ma, which represents young or contemporaneous subducted oceanic slabs. The new geochemical results of Re-Os isotopes and PGE concentrations provide further insights into the recycled source mantle components and the role of sulfides in the origin and evolution of the primitive magma of the Cenozoic Hainan basalts, and the geodynamics of mantle plume activity.

    recycled oceanic crust; PGE geochemistry; Re-Os isotopes; Cenozoic basalts; Hainan Island

    10.16539/j.ddgzyckx.2022.03.006

    2021-10-12;

    2022-01-21

    國家基金委NSFC–廣東省聯(lián)合基金項目(U1701641)和國家自然科學(xué)基金項目(41525006、42073032、42021002)聯(lián)合資助。

    王國慶(1994–), 女, 碩士研究生, 地球化學(xué)專業(yè)。E-mail: wanggq@gig.ac.cn

    郭鋒(1971–), 男, 研究員, 主要從事巖石學(xué)和大地構(gòu)造學(xué)研究。E-mail: guofengt@263.net; fengguo@gig.ac.cn

    P581; P595; P597

    A

    1001-1552(2022)03-0483-018

    猜你喜歡
    洋殼源區(qū)海南島
    洋殼活動斷層特征和動力作用
    冬小麥蒸散源區(qū)代表性分析
    印度洋、大西洋洋殼流運動對地形地貌的影響
    巖漿與構(gòu)造作用對洋殼厚度的影響
    ——以西北印度洋為例
    太平洋洋殼流運動對地形地貌的影響
    渭河源區(qū)徑流量變化特征及趨勢分析
    是海南省還是海南島?
    關(guān)于海南島戰(zhàn)役作戰(zhàn)方針的考察
    軍事歷史(1994年6期)1994-08-15 08:56:38
    回憶海南島渡海登陸戰(zhàn)役
    軍事歷史(1987年2期)1987-08-20 06:10:48
    海南島戰(zhàn)役中的“土艦隊”
    軍事歷史(1981年1期)1981-08-21 03:38:54
    国产真人三级小视频在线观看| 美女福利国产在线| 妹子高潮喷水视频| 日本av手机在线免费观看| 黄片播放在线免费| 国产高清国产精品国产三级| 精品一区在线观看国产| 欧美国产精品一级二级三级| 日本vs欧美在线观看视频| 99精品久久久久人妻精品| 最黄视频免费看| 日韩制服骚丝袜av| 2018国产大陆天天弄谢| 精品少妇久久久久久888优播| 99久久综合免费| 久久国产精品影院| 老熟女久久久| 蜜桃在线观看..| 岛国毛片在线播放| 美女脱内裤让男人舔精品视频| 精品福利永久在线观看| 中国美女看黄片| 最黄视频免费看| 久久精品久久久久久噜噜老黄| 欧美日韩亚洲国产一区二区在线观看 | www.999成人在线观看| 日韩 欧美 亚洲 中文字幕| 热99久久久久精品小说推荐| 久久久欧美国产精品| 成人午夜精彩视频在线观看| 成年动漫av网址| 777久久人妻少妇嫩草av网站| 男男h啪啪无遮挡| 亚洲国产最新在线播放| 女性被躁到高潮视频| 一本综合久久免费| 日日夜夜操网爽| 91老司机精品| 丁香六月欧美| 热99国产精品久久久久久7| 9热在线视频观看99| 国产精品国产三级国产专区5o| 悠悠久久av| 又紧又爽又黄一区二区| 成人国产av品久久久| 高潮久久久久久久久久久不卡| 91麻豆av在线| 这个男人来自地球电影免费观看| 麻豆国产av国片精品| 亚洲黑人精品在线| 天天操日日干夜夜撸| 高清黄色对白视频在线免费看| 高清视频免费观看一区二区| 久久人妻熟女aⅴ| 久久久久精品人妻al黑| 看十八女毛片水多多多| 久久亚洲精品不卡| 亚洲av日韩精品久久久久久密 | 国产成人精品久久久久久| av视频免费观看在线观看| 日韩熟女老妇一区二区性免费视频| 午夜福利视频在线观看免费| 国产91精品成人一区二区三区 | 久久久久久亚洲精品国产蜜桃av| 免费av中文字幕在线| av又黄又爽大尺度在线免费看| 黑人猛操日本美女一级片| 十分钟在线观看高清视频www| 久久影院123| 青青草视频在线视频观看| 久久久精品国产亚洲av高清涩受| 久久热在线av| 伦理电影免费视频| 欧美日韩国产mv在线观看视频| 夜夜骑夜夜射夜夜干| 9热在线视频观看99| 国产老妇伦熟女老妇高清| 无限看片的www在线观看| 又紧又爽又黄一区二区| 欧美精品亚洲一区二区| 欧美日韩亚洲高清精品| 一级黄片播放器| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成人国产一区在线观看 | 丰满人妻熟妇乱又伦精品不卡| 亚洲精品久久午夜乱码| 成人18禁高潮啪啪吃奶动态图| 18禁黄网站禁片午夜丰满| 国产亚洲午夜精品一区二区久久| 久久久精品国产亚洲av高清涩受| svipshipincom国产片| 亚洲国产中文字幕在线视频| 国精品久久久久久国模美| 国产又爽黄色视频| 精品久久久精品久久久| 国产老妇伦熟女老妇高清| 777久久人妻少妇嫩草av网站| 久热爱精品视频在线9| 亚洲精品成人av观看孕妇| 国产无遮挡羞羞视频在线观看| 欧美av亚洲av综合av国产av| tube8黄色片| 国产av一区二区精品久久| 欧美+亚洲+日韩+国产| 欧美变态另类bdsm刘玥| 高清视频免费观看一区二区| 国产一卡二卡三卡精品| 久久精品国产亚洲av涩爱| 老汉色av国产亚洲站长工具| 久久久精品区二区三区| 国产一区有黄有色的免费视频| 日韩,欧美,国产一区二区三区| 国产主播在线观看一区二区 | 91精品三级在线观看| 天天操日日干夜夜撸| 一级,二级,三级黄色视频| 日本色播在线视频| 久久精品亚洲av国产电影网| 国产一区二区在线观看av| 午夜激情久久久久久久| 黄片播放在线免费| 最近手机中文字幕大全| √禁漫天堂资源中文www| 好男人视频免费观看在线| 五月天丁香电影| 青青草视频在线视频观看| 精品福利观看| 我要看黄色一级片免费的| 国产精品一区二区在线不卡| 最近中文字幕2019免费版| 在线观看免费午夜福利视频| 精品国产乱码久久久久久小说| 制服诱惑二区| 午夜久久久在线观看| 欧美乱码精品一区二区三区| a 毛片基地| 免费日韩欧美在线观看| 又紧又爽又黄一区二区| 中文字幕亚洲精品专区| 一本久久精品| 午夜日韩欧美国产| 亚洲精品av麻豆狂野| 一边摸一边做爽爽视频免费| 在线观看www视频免费| 一个人免费看片子| 国产成人系列免费观看| 日韩熟女老妇一区二区性免费视频| 久久精品熟女亚洲av麻豆精品| 国产一区二区三区综合在线观看| 国产av一区二区精品久久| 大型av网站在线播放| 操美女的视频在线观看| 国产成人一区二区三区免费视频网站 | 一区二区三区精品91| 高清av免费在线| 首页视频小说图片口味搜索 | 精品人妻1区二区| 少妇精品久久久久久久| 久久精品人人爽人人爽视色| 亚洲av电影在线进入| 国产精品成人在线| 午夜av观看不卡| 两个人看的免费小视频| 老司机影院成人| 亚洲一区中文字幕在线| www.自偷自拍.com| 久久亚洲国产成人精品v| 欧美在线黄色| 成在线人永久免费视频| 亚洲国产av新网站| √禁漫天堂资源中文www| 亚洲欧美中文字幕日韩二区| 亚洲av成人不卡在线观看播放网 | 老司机靠b影院| 男女无遮挡免费网站观看| 国产成人精品久久二区二区91| 99热国产这里只有精品6| 亚洲欧美日韩高清在线视频 | 男女下面插进去视频免费观看| 精品人妻熟女毛片av久久网站| 男女边摸边吃奶| 天堂中文最新版在线下载| 亚洲欧美激情在线| 精品久久蜜臀av无| 麻豆av在线久日| 国产精品人妻久久久影院| 国产精品三级大全| 看免费av毛片| 欧美亚洲日本最大视频资源| tube8黄色片| 51午夜福利影视在线观看| 美女高潮到喷水免费观看| 满18在线观看网站| 国产伦人伦偷精品视频| 嫁个100分男人电影在线观看 | 妹子高潮喷水视频| 99re6热这里在线精品视频| 亚洲中文av在线| 99热国产这里只有精品6| 青青草视频在线视频观看| 一二三四在线观看免费中文在| 亚洲国产最新在线播放| 极品少妇高潮喷水抽搐| 精品卡一卡二卡四卡免费| 亚洲专区中文字幕在线| 18禁国产床啪视频网站| 男女免费视频国产| 18禁裸乳无遮挡动漫免费视频| 久久久久久久久免费视频了| 999精品在线视频| 亚洲精品国产av蜜桃| 日韩熟女老妇一区二区性免费视频| 纵有疾风起免费观看全集完整版| 又粗又硬又长又爽又黄的视频| 女人爽到高潮嗷嗷叫在线视频| 一边摸一边做爽爽视频免费| 两个人免费观看高清视频| avwww免费| 久久免费观看电影| 亚洲一区二区三区欧美精品| 成年女人毛片免费观看观看9 | 日日摸夜夜添夜夜爱| 热re99久久精品国产66热6| www.999成人在线观看| 免费看十八禁软件| 欧美精品av麻豆av| 蜜桃国产av成人99| 欧美在线黄色| 久久免费观看电影| 亚洲精品乱久久久久久| 欧美性长视频在线观看| 欧美97在线视频| 99国产综合亚洲精品| 久久久久久久久久久久大奶| 男女国产视频网站| 老司机在亚洲福利影院| 一区二区日韩欧美中文字幕| 久久天躁狠狠躁夜夜2o2o | 亚洲国产看品久久| 捣出白浆h1v1| 欧美日韩福利视频一区二区| 黑人欧美特级aaaaaa片| 亚洲精品国产区一区二| 亚洲,欧美精品.| 午夜福利在线免费观看网站| 国产一区亚洲一区在线观看| 久久久久国产一级毛片高清牌| 老司机亚洲免费影院| 国产在线一区二区三区精| 国产男女超爽视频在线观看| av网站免费在线观看视频| www.熟女人妻精品国产| 不卡av一区二区三区| av欧美777| 国产精品成人在线| 黄色毛片三级朝国网站| 亚洲一区中文字幕在线| 80岁老熟妇乱子伦牲交| 亚洲久久久国产精品| 色婷婷av一区二区三区视频| 亚洲精品国产色婷婷电影| 日韩中文字幕欧美一区二区 | 日本wwww免费看| 色播在线永久视频| 国产野战对白在线观看| 1024香蕉在线观看| 亚洲天堂av无毛| 国产一区亚洲一区在线观看| 精品福利观看| 少妇被粗大的猛进出69影院| 亚洲综合色网址| 久久久久久久精品精品| 18禁观看日本| 你懂的网址亚洲精品在线观看| 别揉我奶头~嗯~啊~动态视频 | 精品少妇内射三级| 丝袜喷水一区| 宅男免费午夜| 69精品国产乱码久久久| 国产成人一区二区三区免费视频网站 | 欧美精品一区二区免费开放| 色视频在线一区二区三区| 久久精品国产a三级三级三级| 日本a在线网址| 在线av久久热| 九色亚洲精品在线播放| 久久亚洲精品不卡| 这个男人来自地球电影免费观看| 成人18禁高潮啪啪吃奶动态图| 精品国产超薄肉色丝袜足j| 90打野战视频偷拍视频| 日韩人妻精品一区2区三区| 免费在线观看日本一区| 精品人妻一区二区三区麻豆| 国产精品国产三级国产专区5o| 精品亚洲成国产av| 国产女主播在线喷水免费视频网站| 极品人妻少妇av视频| 亚洲国产最新在线播放| 亚洲午夜精品一区,二区,三区| 两个人免费观看高清视频| 国产精品一区二区在线观看99| 国产精品偷伦视频观看了| 国产成人免费观看mmmm| 精品一区二区三区四区五区乱码 | 伊人久久大香线蕉亚洲五| 久久99一区二区三区| 9色porny在线观看| 狠狠精品人妻久久久久久综合| 久久久欧美国产精品| 赤兔流量卡办理| 婷婷色综合大香蕉| 国产在视频线精品| 18禁国产床啪视频网站| 两人在一起打扑克的视频| 777久久人妻少妇嫩草av网站| 国产精品九九99| 可以免费在线观看a视频的电影网站| 一边摸一边抽搐一进一出视频| 亚洲成av片中文字幕在线观看| 精品久久蜜臀av无| 久久久久久久大尺度免费视频| 国产高清不卡午夜福利| 国产精品亚洲av一区麻豆| 日日爽夜夜爽网站| 午夜福利一区二区在线看| 最黄视频免费看| 成年人午夜在线观看视频| www.av在线官网国产| 久久久国产一区二区| 91字幕亚洲| 一级片免费观看大全| 午夜免费鲁丝| avwww免费| 欧美日韩一级在线毛片| 国语对白做爰xxxⅹ性视频网站| 1024香蕉在线观看| 首页视频小说图片口味搜索 | 女人久久www免费人成看片| 中文精品一卡2卡3卡4更新| 999精品在线视频| 色94色欧美一区二区| 国产精品人妻久久久影院| 一本综合久久免费| 青春草视频在线免费观看| av网站免费在线观看视频| 男女边吃奶边做爰视频| 黄网站色视频无遮挡免费观看| 一本一本久久a久久精品综合妖精| 精品熟女少妇八av免费久了| 亚洲 欧美一区二区三区| 日本五十路高清| 亚洲七黄色美女视频| 在线观看国产h片| 老司机在亚洲福利影院| 高清视频免费观看一区二区| 一区在线观看完整版| 久久精品国产综合久久久| 精品人妻在线不人妻| 国产成人一区二区三区免费视频网站 | 欧美精品av麻豆av| 亚洲图色成人| 婷婷色麻豆天堂久久| 午夜91福利影院| 另类亚洲欧美激情| 嫁个100分男人电影在线观看 | 亚洲成人国产一区在线观看 | 中文字幕制服av| 国产精品一区二区在线观看99| 91麻豆av在线| 婷婷色综合www| 99re6热这里在线精品视频| 国产男人的电影天堂91| 啦啦啦啦在线视频资源| 亚洲人成77777在线视频| 婷婷丁香在线五月| 侵犯人妻中文字幕一二三四区| 国产精品二区激情视频| 亚洲精品乱久久久久久| 精品国产超薄肉色丝袜足j| 国产亚洲欧美在线一区二区| 亚洲免费av在线视频| 欧美人与性动交α欧美软件| 精品人妻在线不人妻| 日韩伦理黄色片| 国产精品免费视频内射| 色94色欧美一区二区| 伦理电影免费视频| 91字幕亚洲| 国产精品国产三级专区第一集| 亚洲国产最新在线播放| 国产又色又爽无遮挡免| 亚洲欧美成人综合另类久久久| a 毛片基地| 免费高清在线观看日韩| 国产成人av激情在线播放| 欧美av亚洲av综合av国产av| 欧美日韩一级在线毛片| 又紧又爽又黄一区二区| 大陆偷拍与自拍| 大码成人一级视频| av不卡在线播放| 纯流量卡能插随身wifi吗| 亚洲欧美精品自产自拍| 亚洲欧美日韩高清在线视频 | 亚洲中文日韩欧美视频| 熟女av电影| 多毛熟女@视频| 精品人妻一区二区三区麻豆| bbb黄色大片| 妹子高潮喷水视频| 最新的欧美精品一区二区| 高清不卡的av网站| 亚洲,欧美,日韩| 大片电影免费在线观看免费| 午夜视频精品福利| 国产黄色免费在线视频| 日本欧美国产在线视频| 高清视频免费观看一区二区| 亚洲欧洲日产国产| 极品少妇高潮喷水抽搐| 午夜av观看不卡| 久久人人爽人人片av| 大话2 男鬼变身卡| 日日爽夜夜爽网站| 国产精品国产av在线观看| 一本色道久久久久久精品综合| 亚洲欧美一区二区三区国产| 90打野战视频偷拍视频| 久久影院123| 久久久国产欧美日韩av| 香蕉国产在线看| www日本在线高清视频| 美女高潮到喷水免费观看| 欧美av亚洲av综合av国产av| 国产欧美日韩精品亚洲av| 青草久久国产| 黄色片一级片一级黄色片| 一本综合久久免费| 欧美在线黄色| 国产麻豆69| 飞空精品影院首页| 国产高清不卡午夜福利| 大香蕉久久成人网| 不卡av一区二区三区| 两个人免费观看高清视频| 99热网站在线观看| 欧美亚洲日本最大视频资源| 夫妻性生交免费视频一级片| 日本av手机在线免费观看| 老司机靠b影院| 制服诱惑二区| 又紧又爽又黄一区二区| 波野结衣二区三区在线| 91精品国产国语对白视频| 亚洲人成电影观看| 99香蕉大伊视频| 久久久久网色| 国产三级黄色录像| 一级片'在线观看视频| 亚洲欧美激情在线| 老汉色av国产亚洲站长工具| 日韩中文字幕视频在线看片| 女人爽到高潮嗷嗷叫在线视频| 性高湖久久久久久久久免费观看| 天堂8中文在线网| 亚洲精品国产区一区二| 欧美日韩视频精品一区| www.精华液| 精品国产一区二区三区四区第35| 七月丁香在线播放| 久久青草综合色| 免费在线观看完整版高清| 大香蕉久久成人网| 国产精品国产三级国产专区5o| av天堂在线播放| 久久这里只有精品19| 热re99久久精品国产66热6| 欧美激情 高清一区二区三区| 国产成人啪精品午夜网站| 91字幕亚洲| 亚洲人成77777在线视频| 久久国产亚洲av麻豆专区| 久久亚洲国产成人精品v| 99热国产这里只有精品6| 欧美 日韩 精品 国产| 国产免费现黄频在线看| 中文字幕av电影在线播放| 亚洲国产日韩一区二区| 亚洲天堂av无毛| 国产精品久久久久久精品古装| 国产伦人伦偷精品视频| 伦理电影免费视频| 考比视频在线观看| 久久久久久亚洲精品国产蜜桃av| 午夜免费成人在线视频| 日韩人妻精品一区2区三区| 老司机在亚洲福利影院| 99国产综合亚洲精品| 婷婷丁香在线五月| 波多野结衣av一区二区av| 国产精品偷伦视频观看了| 精品国产超薄肉色丝袜足j| 国产精品成人在线| 国产精品99久久99久久久不卡| 在线精品无人区一区二区三| 自拍欧美九色日韩亚洲蝌蚪91| 国产xxxxx性猛交| 最新的欧美精品一区二区| 国产亚洲av高清不卡| 大型av网站在线播放| 亚洲国产精品成人久久小说| 高清视频免费观看一区二区| 久久99热这里只频精品6学生| 亚洲第一青青草原| 亚洲人成电影免费在线| 婷婷成人精品国产| 久久人妻福利社区极品人妻图片 | 悠悠久久av| 超色免费av| 国产精品久久久久久精品古装| 亚洲五月婷婷丁香| 色婷婷久久久亚洲欧美| 亚洲国产看品久久| 国产精品免费大片| 欧美日韩亚洲高清精品| 免费高清在线观看视频在线观看| 亚洲精品成人av观看孕妇| 日本欧美国产在线视频| 午夜福利在线免费观看网站| 狂野欧美激情性xxxx| 大片免费播放器 马上看| 十八禁高潮呻吟视频| 亚洲成国产人片在线观看| 国产又爽黄色视频| 日韩av免费高清视频| 精品人妻一区二区三区麻豆| 亚洲国产欧美在线一区| 只有这里有精品99| 国产片内射在线| 人人妻,人人澡人人爽秒播 | 大型av网站在线播放| 国产在视频线精品| 人人妻人人澡人人爽人人夜夜| 真人做人爱边吃奶动态| 五月开心婷婷网| 亚洲av片天天在线观看| 国产日韩欧美在线精品| av网站在线播放免费| 国产高清视频在线播放一区 | 日本wwww免费看| 国产免费一区二区三区四区乱码| 成人亚洲欧美一区二区av| av在线播放精品| 看十八女毛片水多多多| 国产在视频线精品| 看十八女毛片水多多多| 国产精品久久久久成人av| 亚洲av欧美aⅴ国产| 精品少妇久久久久久888优播| 亚洲中文av在线| 国产成人免费无遮挡视频| 成在线人永久免费视频| 欧美日韩视频高清一区二区三区二| 两个人免费观看高清视频| 国产一区二区 视频在线| videos熟女内射| av网站在线播放免费| av视频免费观看在线观看| 亚洲欧美日韩高清在线视频 | 高潮久久久久久久久久久不卡| 精品国产乱码久久久久久男人| 亚洲精品久久午夜乱码| 看免费成人av毛片| 亚洲,一卡二卡三卡| 久久精品aⅴ一区二区三区四区| 国产视频首页在线观看| 最近手机中文字幕大全| 交换朋友夫妻互换小说| 丝袜美足系列| 国产成人精品在线电影| 不卡av一区二区三区| 天堂中文最新版在线下载| 欧美国产精品一级二级三级| 一二三四社区在线视频社区8| 91精品国产国语对白视频| 考比视频在线观看| 亚洲,一卡二卡三卡| 国产欧美日韩精品亚洲av| 咕卡用的链子| cao死你这个sao货| 久久精品久久久久久久性| 日本av免费视频播放| 亚洲伊人色综图| 亚洲国产中文字幕在线视频| 亚洲午夜精品一区,二区,三区| 久久精品亚洲av国产电影网| 老熟女久久久| 久久精品成人免费网站| 少妇粗大呻吟视频| 久久久久视频综合| 亚洲国产成人一精品久久久| 成人国语在线视频| 国产精品久久久人人做人人爽| 久久午夜综合久久蜜桃| 亚洲成人免费av在线播放| 国产成人啪精品午夜网站| 国产在线视频一区二区| 青春草亚洲视频在线观看| 岛国毛片在线播放| 老司机亚洲免费影院| 亚洲人成电影观看| 建设人人有责人人尽责人人享有的| 汤姆久久久久久久影院中文字幕| 免费在线观看视频国产中文字幕亚洲 |