• <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
    a级一级毛片免费在线观看| 99久久精品热视频| 精品一区在线观看国产| 女人十人毛片免费观看3o分钟| 97热精品久久久久久| 国产片特级美女逼逼视频| 啦啦啦啦在线视频资源| 极品少妇高潮喷水抽搐| 岛国毛片在线播放| 国产黄频视频在线观看| 国产黄a三级三级三级人| 91在线精品国自产拍蜜月| 美女cb高潮喷水在线观看| 国产精品伦人一区二区| 久久99热6这里只有精品| 99久久精品国产国产毛片| 搡老妇女老女人老熟妇| 在线观看一区二区三区| 22中文网久久字幕| 日本色播在线视频| 偷拍熟女少妇极品色| 18禁裸乳无遮挡免费网站照片| 日韩一区二区三区影片| 日日啪夜夜爽| 黄色配什么色好看| 中文欧美无线码| 爱豆传媒免费全集在线观看| 女人十人毛片免费观看3o分钟| 成人性生交大片免费视频hd| 久久久久久久亚洲中文字幕| 亚洲av中文字字幕乱码综合| 丝袜喷水一区| 亚洲欧美日韩东京热| 精品一区二区三区人妻视频| 成人国产麻豆网| 舔av片在线| 午夜老司机福利剧场| a级毛色黄片| 夜夜爽夜夜爽视频| 国产黄片美女视频| 久久精品国产自在天天线| 床上黄色一级片| av卡一久久| 三级国产精品片| 亚洲av成人精品一区久久| 少妇熟女aⅴ在线视频| 99视频精品全部免费 在线| 国产精品久久久久久久久免| 亚洲美女视频黄频| 在线观看一区二区三区| 成人国产麻豆网| 久久国产乱子免费精品| 久久99蜜桃精品久久| 91午夜精品亚洲一区二区三区| 永久免费av网站大全| 亚洲欧美精品自产自拍| 久久久久久久久久久免费av| 黄色配什么色好看| 成人亚洲精品一区在线观看 | videos熟女内射| 在线观看人妻少妇| 日本免费在线观看一区| 亚洲美女视频黄频| 国产成人a∨麻豆精品| 久久久久九九精品影院| 99久久人妻综合| 日韩国内少妇激情av| 国产精品综合久久久久久久免费| 丰满人妻一区二区三区视频av| 日韩成人av中文字幕在线观看| 免费观看的影片在线观看| 久久午夜福利片| 免费看日本二区| 国产午夜精品一二区理论片| 夜夜看夜夜爽夜夜摸| 国产 一区精品| 欧美精品国产亚洲| 久久国内精品自在自线图片| 丝瓜视频免费看黄片| 尾随美女入室| 肉色欧美久久久久久久蜜桃 | 又大又黄又爽视频免费| 一个人免费在线观看电影| 久久久久精品性色| 成人综合一区亚洲| 永久免费av网站大全| 色吧在线观看| 日日摸夜夜添夜夜爱| 免费在线观看成人毛片| 欧美高清成人免费视频www| 国产精品一区二区性色av| 热99在线观看视频| 毛片一级片免费看久久久久| 午夜福利成人在线免费观看| 黄色配什么色好看| 午夜日本视频在线| 听说在线观看完整版免费高清| 亚洲精品日本国产第一区| 国产在视频线在精品| 久久久久国产网址| 亚洲精品日韩av片在线观看| 搡老乐熟女国产| 亚洲最大成人手机在线| 丝袜喷水一区| 身体一侧抽搐| 国产精品国产三级国产av玫瑰| 精品久久久久久久末码| 免费黄网站久久成人精品| 国产伦一二天堂av在线观看| 国产精品一区二区性色av| 日韩制服骚丝袜av| 一级毛片我不卡| 精品久久久久久成人av| 成人亚洲欧美一区二区av| av免费观看日本| 亚洲av一区综合| 在线免费观看不下载黄p国产| 99热这里只有是精品50| 美女脱内裤让男人舔精品视频| 日本与韩国留学比较| 精品久久久久久久久av| 男人和女人高潮做爰伦理| 伊人久久国产一区二区| 熟女电影av网| 国产成人精品婷婷| 日韩 亚洲 欧美在线| 国产淫语在线视频| 国产永久视频网站| 国产一区二区亚洲精品在线观看| 午夜激情欧美在线| 免费无遮挡裸体视频| 国产不卡一卡二| 欧美激情在线99| 嘟嘟电影网在线观看| 国产伦在线观看视频一区| 只有这里有精品99| 青春草亚洲视频在线观看| 51国产日韩欧美| 秋霞伦理黄片| 午夜日本视频在线| 国产精品女同一区二区软件| 女的被弄到高潮叫床怎么办| 夫妻性生交免费视频一级片| av又黄又爽大尺度在线免费看| 久久久久九九精品影院| 日韩欧美一区视频在线观看 | 免费看光身美女| 伊人久久精品亚洲午夜| 最近最新中文字幕免费大全7| 国产在视频线精品| 夜夜爽夜夜爽视频| 国产精品国产三级专区第一集| 久久综合国产亚洲精品| 蜜桃久久精品国产亚洲av| 久久久久国产网址| 少妇丰满av| 日日摸夜夜添夜夜爱| 亚洲精品国产av蜜桃| 国产精品一区二区三区四区免费观看| 99热网站在线观看| 精品一区二区三区视频在线| 91精品一卡2卡3卡4卡| 2022亚洲国产成人精品| 一级毛片我不卡| 国产日韩欧美在线精品| 亚洲av中文av极速乱| 狂野欧美激情性xxxx在线观看| 少妇人妻精品综合一区二区| 99久久九九国产精品国产免费| 婷婷六月久久综合丁香| 亚洲欧美清纯卡通| 国产成人精品一,二区| 亚洲欧美精品自产自拍| 国产一区二区亚洲精品在线观看| 高清在线视频一区二区三区| 国产高潮美女av| 精品久久国产蜜桃| 日本黄色片子视频| 人人妻人人澡人人爽人人夜夜 | 国产亚洲最大av| 国产黄色视频一区二区在线观看| av在线亚洲专区| 美女大奶头视频| 亚洲国产精品成人综合色| 欧美高清性xxxxhd video| 大又大粗又爽又黄少妇毛片口| 80岁老熟妇乱子伦牲交| 久久久久久久久久人人人人人人| 国产亚洲91精品色在线| 亚洲国产精品成人久久小说| 91精品一卡2卡3卡4卡| 日产精品乱码卡一卡2卡三| 亚洲aⅴ乱码一区二区在线播放| 亚洲不卡免费看| 最近中文字幕2019免费版| 婷婷六月久久综合丁香| 99久国产av精品国产电影| 亚洲丝袜综合中文字幕| 亚洲精品日韩在线中文字幕| 69av精品久久久久久| 九九爱精品视频在线观看| 国产老妇女一区| 日韩三级伦理在线观看| 国产成人aa在线观看| 亚洲国产日韩欧美精品在线观看| 日韩制服骚丝袜av| 亚洲成人久久爱视频| 免费看不卡的av| 亚洲成人精品中文字幕电影| 禁无遮挡网站| 中文字幕av在线有码专区| 在线观看av片永久免费下载| 成人漫画全彩无遮挡| 国产大屁股一区二区在线视频| 国产精品一区二区性色av| 国产高清国产精品国产三级 | 91久久精品国产一区二区三区| a级毛色黄片| 成年版毛片免费区| 国产精品日韩av在线免费观看| 午夜福利视频1000在线观看| 日韩一区二区视频免费看| av免费观看日本| 亚洲综合精品二区| 亚洲人成网站在线观看播放| 男人舔女人下体高潮全视频| 日本一二三区视频观看| 久久久久久久久中文| 国产黄频视频在线观看| 免费观看无遮挡的男女| 99热这里只有精品一区| 午夜福利高清视频| 爱豆传媒免费全集在线观看| 色播亚洲综合网| 三级毛片av免费| 精品99又大又爽又粗少妇毛片| 亚洲欧美成人精品一区二区| 男女下面进入的视频免费午夜| 一级毛片我不卡| 精品久久久久久久人妻蜜臀av| 深夜a级毛片| 大片免费播放器 马上看| 一级爰片在线观看| 边亲边吃奶的免费视频| 在线观看免费高清a一片| 国产色婷婷99| 国产三级在线视频| 美女cb高潮喷水在线观看| 赤兔流量卡办理| av免费观看日本| 插阴视频在线观看视频| 免费av观看视频| 日本午夜av视频| 亚洲va在线va天堂va国产| 国产欧美日韩精品一区二区| 国产伦精品一区二区三区视频9| 久久久色成人| 日本熟妇午夜| 一边亲一边摸免费视频| 成年免费大片在线观看| 看免费成人av毛片| 国产一区二区亚洲精品在线观看| 亚洲久久久久久中文字幕| 三级国产精品欧美在线观看| 欧美性感艳星| 天堂影院成人在线观看| 中文字幕av成人在线电影| 日本爱情动作片www.在线观看| 乱码一卡2卡4卡精品| 男女啪啪激烈高潮av片| 女人十人毛片免费观看3o分钟| 精品久久久噜噜| 只有这里有精品99| 精品欧美国产一区二区三| 亚洲欧美成人精品一区二区| 久久久久久久久久黄片| 成人亚洲精品av一区二区| 久久久亚洲精品成人影院| 日韩av在线大香蕉| a级一级毛片免费在线观看| 又大又黄又爽视频免费| a级毛色黄片| 九九在线视频观看精品| 欧美日韩视频高清一区二区三区二| 一级毛片aaaaaa免费看小| 日本色播在线视频| 午夜爱爱视频在线播放| 亚洲成人av在线免费| 久久久久久国产a免费观看| 乱系列少妇在线播放| 欧美日韩综合久久久久久| 搡老妇女老女人老熟妇| 久久精品国产鲁丝片午夜精品| 自拍偷自拍亚洲精品老妇| 91在线精品国自产拍蜜月| 国产免费又黄又爽又色| 精品一区二区免费观看| 99热这里只有是精品在线观看| 亚洲第一区二区三区不卡| 搡老妇女老女人老熟妇| 亚洲三级黄色毛片| 亚洲av男天堂| 女人被狂操c到高潮| 一级毛片aaaaaa免费看小| 哪个播放器可以免费观看大片| 2021少妇久久久久久久久久久| 尾随美女入室| 久久精品夜夜夜夜夜久久蜜豆| 午夜免费激情av| 国产成人午夜福利电影在线观看| 美女国产视频在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲美女搞黄在线观看| 成人美女网站在线观看视频| 少妇熟女aⅴ在线视频| 国产黄色免费在线视频| av又黄又爽大尺度在线免费看| 日本猛色少妇xxxxx猛交久久| 在线观看免费高清a一片| 亚洲美女搞黄在线观看| 久久精品国产亚洲av涩爱| 日韩一本色道免费dvd| 亚洲国产欧美在线一区| 午夜精品在线福利| 国产高清三级在线| 大陆偷拍与自拍| 精华霜和精华液先用哪个| 嫩草影院新地址| 欧美极品一区二区三区四区| 白带黄色成豆腐渣| 亚洲精品成人久久久久久| 日韩成人av中文字幕在线观看| 国产色婷婷99| 国产伦精品一区二区三区视频9| 麻豆精品久久久久久蜜桃| 国产精品1区2区在线观看.| 久久6这里有精品| 国产精品人妻久久久久久| 亚洲av免费在线观看| 久久草成人影院| 永久免费av网站大全| 亚洲婷婷狠狠爱综合网| 亚洲人成网站在线观看播放| 亚洲av男天堂| 国产麻豆成人av免费视频| 直男gayav资源| 51国产日韩欧美| 欧美日韩亚洲高清精品| 内射极品少妇av片p| 成人国产麻豆网| 亚洲国产色片| 大香蕉久久网| 国产精品三级大全| 99热这里只有精品一区| 亚洲欧美精品自产自拍| 99九九线精品视频在线观看视频| 日韩欧美国产在线观看| 自拍偷自拍亚洲精品老妇| 日本欧美国产在线视频| 麻豆成人av视频| 久久久久久伊人网av| 国产精品伦人一区二区| 亚洲第一区二区三区不卡| 18禁在线播放成人免费| 韩国av在线不卡| 欧美激情久久久久久爽电影| 神马国产精品三级电影在线观看| 久久久久久国产a免费观看| 欧美3d第一页| 欧美人与善性xxx| 国产精品福利在线免费观看| 国产黄色小视频在线观看| 91精品国产九色| 熟妇人妻久久中文字幕3abv| 亚洲欧美日韩无卡精品| 少妇裸体淫交视频免费看高清| 亚洲精品第二区| 欧美高清成人免费视频www| 天堂√8在线中文| 精品久久久久久久久久久久久| 麻豆久久精品国产亚洲av| 日韩av在线大香蕉| 久久久久久久久久久免费av| 超碰av人人做人人爽久久| 国产淫语在线视频| 亚洲国产欧美在线一区| 一级毛片aaaaaa免费看小| 日韩大片免费观看网站| 中文在线观看免费www的网站| 91aial.com中文字幕在线观看| 简卡轻食公司| 日本免费在线观看一区| 一个人观看的视频www高清免费观看| 国产伦理片在线播放av一区| 听说在线观看完整版免费高清| 老女人水多毛片| 国产精品人妻久久久影院| 亚洲天堂国产精品一区在线| 亚洲综合色惰| 国产在线一区二区三区精| 久久99蜜桃精品久久| 亚洲va在线va天堂va国产| av在线老鸭窝| 国产伦精品一区二区三区视频9| 如何舔出高潮| 亚洲乱码一区二区免费版| 国产片特级美女逼逼视频| 欧美区成人在线视频| 久久韩国三级中文字幕| 熟妇人妻久久中文字幕3abv| 日韩强制内射视频| 免费人成在线观看视频色| 日本wwww免费看| 亚洲国产成人一精品久久久| 蜜臀久久99精品久久宅男| 五月天丁香电影| 国产老妇女一区| 亚洲av电影在线观看一区二区三区 | 亚洲欧美成人精品一区二区| 青春草亚洲视频在线观看| 中文资源天堂在线| 91在线精品国自产拍蜜月| 看免费成人av毛片| 午夜老司机福利剧场| 国产成人精品久久久久久| 亚洲精品日本国产第一区| 精品人妻偷拍中文字幕| 韩国高清视频一区二区三区| 91av网一区二区| 亚洲精品乱码久久久v下载方式| 欧美bdsm另类| 真实男女啪啪啪动态图| 国产国拍精品亚洲av在线观看| 中文字幕av在线有码专区| 欧美一区二区亚洲| 久久精品国产亚洲av涩爱| 一区二区三区四区激情视频| 老司机影院毛片| 波多野结衣巨乳人妻| 最近视频中文字幕2019在线8| 麻豆成人av视频| 日韩成人av中文字幕在线观看| 国产精品一区二区三区四区免费观看| 亚洲综合精品二区| 天天躁夜夜躁狠狠久久av| 水蜜桃什么品种好| 少妇人妻一区二区三区视频| 亚洲av电影不卡..在线观看| 色综合亚洲欧美另类图片| 日本一本二区三区精品| 国产一区二区亚洲精品在线观看| 国产精品99久久久久久久久| 国产精品久久久久久av不卡| 日本黄大片高清| 51国产日韩欧美| 精品人妻视频免费看| 亚洲精品国产av蜜桃| 高清日韩中文字幕在线| 丝袜美腿在线中文| 伦理电影大哥的女人| 久久韩国三级中文字幕| 身体一侧抽搐| 久久精品国产亚洲网站| 成人国产麻豆网| 观看免费一级毛片| 亚洲最大成人av| 日韩,欧美,国产一区二区三区| 男女视频在线观看网站免费| 大话2 男鬼变身卡| 欧美丝袜亚洲另类| 亚洲成人中文字幕在线播放| 欧美高清性xxxxhd video| 亚洲自偷自拍三级| 80岁老熟妇乱子伦牲交| 国精品久久久久久国模美| 老司机影院毛片| 精品一区二区三区视频在线| 国产精品福利在线免费观看| 亚洲精品乱码久久久v下载方式| 国产欧美另类精品又又久久亚洲欧美| 亚洲精华国产精华液的使用体验| av在线播放精品| 亚洲成人av在线免费| 能在线免费观看的黄片| 精品人妻视频免费看| av黄色大香蕉| 国产黄a三级三级三级人| 亚洲性久久影院| 嫩草影院入口| 国产高潮美女av| 女人被狂操c到高潮| 国产精品久久久久久精品电影| 国产淫语在线视频| 国产黄片视频在线免费观看| 九九在线视频观看精品| 51国产日韩欧美| 黄片wwwwww| 亚洲人成网站高清观看| 深爱激情五月婷婷| av福利片在线观看| 国产麻豆成人av免费视频| 五月伊人婷婷丁香| 老师上课跳d突然被开到最大视频| 97在线视频观看| 欧美激情久久久久久爽电影| 五月天丁香电影| 亚洲熟女精品中文字幕| 久久久久久久亚洲中文字幕| 欧美性猛交╳xxx乱大交人| 亚洲av二区三区四区| 成人性生交大片免费视频hd| 哪个播放器可以免费观看大片| 国产大屁股一区二区在线视频| 淫秽高清视频在线观看| 波野结衣二区三区在线| 最近的中文字幕免费完整| 成人国产麻豆网| av专区在线播放| 一二三四中文在线观看免费高清| 一区二区三区高清视频在线| 成人亚洲精品av一区二区| 天堂√8在线中文| 日韩,欧美,国产一区二区三区| 午夜精品在线福利| av免费观看日本| 国产免费又黄又爽又色| 免费观看性生交大片5| 夜夜看夜夜爽夜夜摸| 三级国产精品片| 国产成人91sexporn| 少妇的逼水好多| 欧美成人a在线观看| 91午夜精品亚洲一区二区三区| 精品99又大又爽又粗少妇毛片| 国产国拍精品亚洲av在线观看| 不卡视频在线观看欧美| 精品一区二区三区视频在线| 日韩中字成人| 亚洲三级黄色毛片| 日韩av在线大香蕉| 夜夜看夜夜爽夜夜摸| 亚洲欧美成人精品一区二区| 日韩av免费高清视频| 国产毛片a区久久久久| 婷婷色综合www| 国产午夜精品久久久久久一区二区三区| 免费黄频网站在线观看国产| 在线观看美女被高潮喷水网站| 亚洲第一区二区三区不卡| 一级av片app| 亚洲av成人av| 成人高潮视频无遮挡免费网站| 亚洲国产色片| 草草在线视频免费看| 欧美高清成人免费视频www| 日韩欧美精品免费久久| 街头女战士在线观看网站| 99热这里只有是精品50| 两个人视频免费观看高清| 91久久精品电影网| 久久精品国产亚洲av天美| 丰满乱子伦码专区| 欧美区成人在线视频| 啦啦啦啦在线视频资源| 又粗又硬又长又爽又黄的视频| 国产成人精品福利久久| 午夜激情欧美在线| 最近最新中文字幕免费大全7| 欧美+日韩+精品| 亚洲精品日本国产第一区| 日产精品乱码卡一卡2卡三| 国产v大片淫在线免费观看| 久久精品久久久久久噜噜老黄| 狠狠精品人妻久久久久久综合| 99热网站在线观看| 久久精品熟女亚洲av麻豆精品 | 中文乱码字字幕精品一区二区三区 | 午夜激情久久久久久久| 美女被艹到高潮喷水动态| 乱人视频在线观看| 人妻制服诱惑在线中文字幕| 国产精品久久久久久av不卡| 成人亚洲欧美一区二区av| 大香蕉97超碰在线| 亚洲精品乱码久久久久久按摩| 国产亚洲一区二区精品| 亚洲欧美成人综合另类久久久| 亚洲av在线观看美女高潮| 亚洲av免费在线观看| 波多野结衣巨乳人妻| 嫩草影院新地址| 人体艺术视频欧美日本| 午夜精品在线福利| 亚洲精品日韩在线中文字幕| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| 成人综合一区亚洲| 欧美性感艳星| 黄色日韩在线| av女优亚洲男人天堂| 国产极品天堂在线| 亚洲美女搞黄在线观看| 国产男人的电影天堂91| 赤兔流量卡办理| 最近中文字幕高清免费大全6| 国产精品久久久久久久电影| 久99久视频精品免费| 精品国产三级普通话版| 人妻一区二区av| 久久久久网色| 婷婷六月久久综合丁香| 特级一级黄色大片| 51国产日韩欧美| 亚州av有码| 一个人看的www免费观看视频|