• <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国产免费在线观看| 欧美日韩综合久久久久久| 国产乱来视频区| 欧美精品国产亚洲| 亚洲欧美日韩无卡精品| 九九在线视频观看精品| 日韩,欧美,国产一区二区三区| 丝瓜视频免费看黄片| 成人综合一区亚洲| 美女高潮的动态| 秋霞在线观看毛片| 亚洲av免费高清在线观看| 国产精品av视频在线免费观看| 日本免费在线观看一区| 黄色日韩在线| 97超碰精品成人国产| 亚洲av在线观看美女高潮| 极品少妇高潮喷水抽搐| 亚洲真实伦在线观看| 麻豆成人午夜福利视频| 毛片女人毛片| 少妇的逼好多水| 99热网站在线观看| 国产高潮美女av| 边亲边吃奶的免费视频| 国产在线一区二区三区精| 国产成人a区在线观看| av卡一久久| 亚洲欧美日韩卡通动漫| 五月伊人婷婷丁香| 国产黄频视频在线观看| 夫妻性生交免费视频一级片| 精品熟女少妇av免费看| 伊人久久精品亚洲午夜| 18+在线观看网站| 毛片一级片免费看久久久久| 色播亚洲综合网| 国内精品美女久久久久久| 日韩欧美一区视频在线观看 | 岛国毛片在线播放| 欧美激情国产日韩精品一区| 亚洲美女搞黄在线观看| 内地一区二区视频在线| 亚洲精品中文字幕在线视频 | 国产成人a区在线观看| 最近2019中文字幕mv第一页| 久久久久性生活片| 午夜福利视频1000在线观看| 99久久精品国产国产毛片| 国产亚洲午夜精品一区二区久久 | 尾随美女入室| 久久久久九九精品影院| 嫩草影院精品99| 69人妻影院| 欧美日韩视频精品一区| 九九久久精品国产亚洲av麻豆| 欧美日韩精品成人综合77777| 精品久久久噜噜| 男人和女人高潮做爰伦理| av在线天堂中文字幕| 国产精品国产三级专区第一集| 国产黄频视频在线观看| 欧美区成人在线视频| 国产黄色免费在线视频| 成人亚洲欧美一区二区av| 欧美老熟妇乱子伦牲交| 亚洲国产日韩一区二区| 国产免费又黄又爽又色| 久久精品国产a三级三级三级| 人妻少妇偷人精品九色| 日本欧美国产在线视频| xxx大片免费视频| 国产精品精品国产色婷婷| 免费观看a级毛片全部| 亚洲欧美成人综合另类久久久| 国产精品国产三级国产av玫瑰| 国产精品成人在线| 免费av不卡在线播放| 免费黄网站久久成人精品| 18禁在线播放成人免费| 日日撸夜夜添| 五月玫瑰六月丁香| 日韩av免费高清视频| 神马国产精品三级电影在线观看| 波野结衣二区三区在线| 午夜精品一区二区三区免费看| 日日摸夜夜添夜夜添av毛片| 国产精品人妻久久久久久| 一个人看的www免费观看视频| av网站免费在线观看视频| 免费av不卡在线播放| 全区人妻精品视频| 欧美高清性xxxxhd video| 一二三四中文在线观看免费高清| 亚洲婷婷狠狠爱综合网| 一本一本综合久久| 欧美一区二区亚洲| 肉色欧美久久久久久久蜜桃 | 久久久久九九精品影院| 亚洲国产色片| 国产91av在线免费观看| 欧美丝袜亚洲另类| 美女cb高潮喷水在线观看| 色婷婷久久久亚洲欧美| 男男h啪啪无遮挡| 欧美日韩国产mv在线观看视频 | 嫩草影院新地址| 日韩中字成人| 国产精品一及| 在线观看三级黄色| 神马国产精品三级电影在线观看| 哪个播放器可以免费观看大片| 国产精品秋霞免费鲁丝片| 国产av码专区亚洲av| 欧美xxⅹ黑人| 中文资源天堂在线| kizo精华| h日本视频在线播放| 久久精品人妻少妇| 99视频精品全部免费 在线| 中文天堂在线官网| 欧美日韩在线观看h| 国产成年人精品一区二区| 国产一区亚洲一区在线观看| 亚洲国产av新网站| 欧美亚洲 丝袜 人妻 在线| 亚洲综合色惰| 精品午夜福利在线看| 欧美性感艳星| 少妇的逼好多水| 美女cb高潮喷水在线观看| 久久久久网色| av在线播放精品| 亚洲精品成人av观看孕妇| 下体分泌物呈黄色| 女的被弄到高潮叫床怎么办| 国产爽快片一区二区三区| 国产精品久久久久久久久免| 亚洲图色成人| 美女高潮的动态| 亚洲三级黄色毛片| 国产又色又爽无遮挡免| 另类亚洲欧美激情| 亚洲无线观看免费| 亚洲电影在线观看av| 国产精品.久久久| 国产黄频视频在线观看| 亚洲真实伦在线观看| 在线观看一区二区三区| 国产精品一及| 久久久精品94久久精品| 五月天丁香电影| 成年免费大片在线观看| 简卡轻食公司| 特大巨黑吊av在线直播| 日韩视频在线欧美| 精品亚洲乱码少妇综合久久| 色5月婷婷丁香| 99视频精品全部免费 在线| 少妇高潮的动态图| 亚洲精品视频女| 久久国产乱子免费精品| 久久ye,这里只有精品| a级毛片免费高清观看在线播放| 老司机影院成人| 国产大屁股一区二区在线视频| 国产成人精品一,二区| 欧美区成人在线视频| 纵有疾风起免费观看全集完整版| 少妇人妻 视频| 欧美亚洲 丝袜 人妻 在线| 久久精品综合一区二区三区| 只有这里有精品99| 男女边吃奶边做爰视频| 最近中文字幕高清免费大全6| 午夜老司机福利剧场| 人妻系列 视频| 亚洲av成人精品一区久久| 三级经典国产精品| 亚洲精品乱码久久久久久按摩| 一级a做视频免费观看| 日本av手机在线免费观看| 亚洲欧美中文字幕日韩二区| 久久久久久久大尺度免费视频| 国产大屁股一区二区在线视频| 国产精品.久久久| 国产白丝娇喘喷水9色精品| 欧美日韩精品成人综合77777| 综合色丁香网| 国国产精品蜜臀av免费| 国产男人的电影天堂91| 国产精品久久久久久av不卡| av女优亚洲男人天堂| 熟妇人妻不卡中文字幕| 国产精品99久久久久久久久| 可以在线观看毛片的网站| av天堂中文字幕网| 国产一区二区三区综合在线观看 | 两个人的视频大全免费| 国产精品秋霞免费鲁丝片| 国产免费视频播放在线视频| 中文字幕制服av| 国产探花在线观看一区二区| 亚洲在久久综合| 欧美日韩一区二区视频在线观看视频在线 | 尤物成人国产欧美一区二区三区| 国产精品不卡视频一区二区| 伦理电影大哥的女人| 国产一区二区在线观看日韩| 亚洲欧美成人综合另类久久久| 听说在线观看完整版免费高清| 在线观看免费高清a一片| 九九爱精品视频在线观看| 大又大粗又爽又黄少妇毛片口| 国产色爽女视频免费观看| h日本视频在线播放| 亚洲欧美精品自产自拍| 日韩制服骚丝袜av| 午夜激情久久久久久久| 午夜福利视频1000在线观看| 国产精品久久久久久精品古装| 国产av国产精品国产| 国产精品一二三区在线看| 久热久热在线精品观看| 精品久久久噜噜| 久久鲁丝午夜福利片| 3wmmmm亚洲av在线观看| 丰满人妻一区二区三区视频av| 亚洲精品国产av成人精品| 乱系列少妇在线播放| 天堂网av新在线| 在线天堂最新版资源| 亚洲av.av天堂| 男女下面进入的视频免费午夜| 久久精品久久精品一区二区三区| 国产人妻一区二区三区在| 好男人在线观看高清免费视频| 乱码一卡2卡4卡精品| 亚洲精品一区蜜桃| 国产熟女欧美一区二区| 亚洲真实伦在线观看| 国产伦精品一区二区三区视频9| 久久人人爽av亚洲精品天堂 | 99热6这里只有精品| 高清毛片免费看| 国产精品.久久久| 中文资源天堂在线| 人体艺术视频欧美日本| 亚洲欧美精品专区久久| 中文精品一卡2卡3卡4更新| 日本午夜av视频| 久久精品久久精品一区二区三区| 成人鲁丝片一二三区免费| 日本一本二区三区精品| 麻豆久久精品国产亚洲av| 97人妻精品一区二区三区麻豆| 国产爱豆传媒在线观看| 国国产精品蜜臀av免费| 在线观看美女被高潮喷水网站| 日韩亚洲欧美综合| 少妇猛男粗大的猛烈进出视频 | 男女下面进入的视频免费午夜| 亚洲欧美日韩无卡精品| 一个人观看的视频www高清免费观看| 性插视频无遮挡在线免费观看| 久久久a久久爽久久v久久| 国产免费一级a男人的天堂| 国产黄频视频在线观看| 国产成人精品婷婷| 精品一区二区三区视频在线| 欧美日韩亚洲高清精品| 国产毛片a区久久久久| 国产精品偷伦视频观看了| av国产精品久久久久影院| 国产一区二区在线观看日韩| 精品久久久久久久末码| 国产大屁股一区二区在线视频| 人妻一区二区av| 日韩精品有码人妻一区| 99热网站在线观看| videossex国产| 中文字幕亚洲精品专区| 真实男女啪啪啪动态图| 白带黄色成豆腐渣| 一个人看的www免费观看视频| 久久久久久伊人网av| 国产黄a三级三级三级人| 男插女下体视频免费在线播放| 亚洲精品国产av成人精品| 国产色爽女视频免费观看| 欧美区成人在线视频| 伦理电影大哥的女人| 亚洲一级一片aⅴ在线观看| videossex国产| 成人综合一区亚洲| 欧美+日韩+精品| 日韩国内少妇激情av| 免费电影在线观看免费观看| 2021少妇久久久久久久久久久| 最后的刺客免费高清国语| 成年女人在线观看亚洲视频 | 男女边吃奶边做爰视频| 建设人人有责人人尽责人人享有的 | 亚洲国产精品专区欧美| 亚洲精品久久午夜乱码| 久久久久久久国产电影| 国产在线男女| 日韩一区二区视频免费看| 国产精品一区www在线观看| 插阴视频在线观看视频| 日韩视频在线欧美| 天堂中文最新版在线下载 | 嫩草影院精品99| 熟女人妻精品中文字幕| 欧美潮喷喷水| av在线观看视频网站免费| 97热精品久久久久久| 赤兔流量卡办理| 国产精品福利在线免费观看| 日韩成人伦理影院| 日韩不卡一区二区三区视频在线| 欧美激情在线99| 亚洲人与动物交配视频| 岛国毛片在线播放| 精品久久久久久电影网| 欧美激情在线99| 成人一区二区视频在线观看| 国产爱豆传媒在线观看| 久久热精品热| 国产精品久久久久久精品电影| 国产美女午夜福利| 亚洲色图av天堂| 夜夜看夜夜爽夜夜摸| 噜噜噜噜噜久久久久久91| 丰满少妇做爰视频| 赤兔流量卡办理| 亚洲天堂av无毛| 午夜福利视频1000在线观看| 日韩在线高清观看一区二区三区| av在线蜜桃| 亚洲av成人精品一二三区| 精品少妇久久久久久888优播| 男人爽女人下面视频在线观看| 亚洲,欧美,日韩| 免费高清在线观看视频在线观看| 国产一区二区亚洲精品在线观看| 免费在线观看成人毛片| 成年版毛片免费区| 天天躁夜夜躁狠狠久久av| 毛片一级片免费看久久久久| 国产精品久久久久久精品古装| 亚洲av二区三区四区| 激情五月婷婷亚洲| 一级av片app| 最近中文字幕2019免费版| 日韩大片免费观看网站| 久久久久性生活片| 成人免费观看视频高清| 一边亲一边摸免费视频| 熟妇人妻不卡中文字幕| 国产女主播在线喷水免费视频网站| 久久久久久久精品精品| 真实男女啪啪啪动态图| 亚洲成人精品中文字幕电影| 亚洲国产精品成人综合色| 老女人水多毛片| 久久热精品热| 青青草视频在线视频观看| 内地一区二区视频在线| 精品人妻偷拍中文字幕| 国产亚洲最大av| 99久久人妻综合| 亚洲欧美精品自产自拍| 久久影院123| av线在线观看网站| 2021天堂中文幕一二区在线观| 大片免费播放器 马上看| 国产精品国产av在线观看| av免费在线看不卡| 国产精品伦人一区二区| 亚洲第一区二区三区不卡| 日韩成人av中文字幕在线观看| 欧美一级a爱片免费观看看| 五月伊人婷婷丁香| 97精品久久久久久久久久精品| 亚洲精华国产精华液的使用体验| 欧美bdsm另类| 精品人妻熟女av久视频| 网址你懂的国产日韩在线| 又粗又硬又长又爽又黄的视频| 蜜桃亚洲精品一区二区三区| 国产老妇女一区| 国产成人福利小说| 99热这里只有精品一区| 男女边摸边吃奶| 精品久久久久久久末码| 久久亚洲国产成人精品v| 自拍偷自拍亚洲精品老妇| 日本一二三区视频观看| 色综合色国产| 欧美另类一区| 亚洲在线观看片| 又大又黄又爽视频免费| 免费在线观看成人毛片| 大话2 男鬼变身卡| 亚洲,一卡二卡三卡| 波多野结衣巨乳人妻| 亚洲av.av天堂| 国产精品久久久久久精品古装| 欧美精品人与动牲交sv欧美| 亚洲最大成人av| 一级爰片在线观看| 日本午夜av视频| 最近的中文字幕免费完整| 在线a可以看的网站| 久久人人爽av亚洲精品天堂 | 特级一级黄色大片| 丰满人妻一区二区三区视频av| 免费看av在线观看网站| 网址你懂的国产日韩在线| 久久综合国产亚洲精品| 一级毛片aaaaaa免费看小| 亚洲自拍偷在线| 老司机影院成人| 五月开心婷婷网| 国产高清国产精品国产三级 | 亚洲国产精品999| 丰满乱子伦码专区| 高清在线视频一区二区三区| 精品国产乱码久久久久久小说| 99热全是精品| 永久网站在线| 男男h啪啪无遮挡| 伊人久久精品亚洲午夜| 国产视频首页在线观看| 精品视频人人做人人爽| 亚州av有码| 91精品国产九色| 日日啪夜夜撸| 婷婷色综合www| 在线a可以看的网站| 国产免费视频播放在线视频| 国产毛片在线视频| 亚洲av免费在线观看| 中文字幕制服av| 国产极品天堂在线| 美女cb高潮喷水在线观看| 亚洲成色77777| 久久6这里有精品| 国产伦理片在线播放av一区| 国国产精品蜜臀av免费| tube8黄色片| 亚洲av.av天堂| 国产精品蜜桃在线观看| 人妻夜夜爽99麻豆av| 最后的刺客免费高清国语| 成年版毛片免费区| 成人漫画全彩无遮挡| 自拍欧美九色日韩亚洲蝌蚪91 | 久久人人爽人人爽人人片va| 蜜桃久久精品国产亚洲av| 自拍偷自拍亚洲精品老妇| 亚洲精品日本国产第一区| 国产高清不卡午夜福利| 激情 狠狠 欧美| 男人添女人高潮全过程视频| 日韩三级伦理在线观看| 精品亚洲乱码少妇综合久久| 少妇人妻 视频| 国产精品久久久久久精品古装| 久久精品人妻少妇| 狂野欧美白嫩少妇大欣赏| 中文天堂在线官网| 精品一区二区三区视频在线| 制服丝袜香蕉在线| 免费观看的影片在线观看| 国产欧美亚洲国产| av播播在线观看一区| 精品久久久精品久久久| 国产精品一区二区在线观看99| 久久影院123| 亚洲在线观看片| 国产精品一及| 国产色爽女视频免费观看| 18禁裸乳无遮挡动漫免费视频 | 九九久久精品国产亚洲av麻豆| 内地一区二区视频在线| 精品国产三级普通话版| 亚洲久久久久久中文字幕| 美女主播在线视频| 国产v大片淫在线免费观看| 免费高清在线观看视频在线观看| 一级黄片播放器| 欧美xxⅹ黑人| 中文乱码字字幕精品一区二区三区| 你懂的网址亚洲精品在线观看| 欧美xxxx黑人xx丫x性爽| .国产精品久久| videossex国产| 国产有黄有色有爽视频| 精品久久久精品久久久| 国产有黄有色有爽视频| 亚洲最大成人手机在线| 99视频精品全部免费 在线| 婷婷色综合www| 亚洲av日韩在线播放| 久久久久精品性色| 亚洲电影在线观看av| 色婷婷久久久亚洲欧美| 纵有疾风起免费观看全集完整版| 禁无遮挡网站| 亚洲国产精品999| 天美传媒精品一区二区| 亚洲国产精品999| 校园人妻丝袜中文字幕| 成人高潮视频无遮挡免费网站| 日日摸夜夜添夜夜爱| .国产精品久久| 欧美成人一区二区免费高清观看| 国产男女超爽视频在线观看| 亚洲综合精品二区| 身体一侧抽搐| 亚洲va在线va天堂va国产| 亚洲国产高清在线一区二区三| 色婷婷久久久亚洲欧美| 三级国产精品欧美在线观看| 黄色视频在线播放观看不卡| 久久久午夜欧美精品| 免费在线观看成人毛片| 精品人妻视频免费看| 国产精品蜜桃在线观看| 国产伦精品一区二区三区视频9| 中文在线观看免费www的网站| 成人国产av品久久久| 亚洲丝袜综合中文字幕| 欧美少妇被猛烈插入视频| a级一级毛片免费在线观看| 久久久精品欧美日韩精品| 内射极品少妇av片p| av在线播放精品| 久久久久久久久久人人人人人人| 成人无遮挡网站| 你懂的网址亚洲精品在线观看| 成年人午夜在线观看视频| 婷婷色麻豆天堂久久| 日本黄色片子视频| 欧美bdsm另类| 又爽又黄a免费视频| 啦啦啦在线观看免费高清www| 人妻 亚洲 视频| 超碰av人人做人人爽久久| 久久精品熟女亚洲av麻豆精品| 国产成人91sexporn| 国产亚洲av片在线观看秒播厂| 欧美精品国产亚洲| 99久久精品一区二区三区| 91久久精品国产一区二区三区| 午夜免费男女啪啪视频观看| 国产精品国产三级国产av玫瑰| 亚洲不卡免费看| 丰满人妻一区二区三区视频av| 国产日韩欧美亚洲二区| 26uuu在线亚洲综合色| 欧美成人精品欧美一级黄| 国产伦精品一区二区三区视频9| 亚洲国产精品国产精品| 国产精品人妻久久久影院| 99久久精品一区二区三区| 视频中文字幕在线观看| 黄色一级大片看看| 午夜老司机福利剧场| 成人毛片a级毛片在线播放| 秋霞在线观看毛片| 亚洲自拍偷在线| 伊人久久国产一区二区| 国产亚洲91精品色在线| 欧美老熟妇乱子伦牲交| 国产免费视频播放在线视频| 国产成人一区二区在线| 大片免费播放器 马上看| 国产免费福利视频在线观看| 午夜激情福利司机影院| eeuss影院久久| 国产色爽女视频免费观看| 日韩强制内射视频| 亚洲色图av天堂| 亚洲精品日本国产第一区| 国产乱人视频| 久热久热在线精品观看| 菩萨蛮人人尽说江南好唐韦庄| 国产淫语在线视频| 久久影院123| 成人亚洲精品一区在线观看 | 美女xxoo啪啪120秒动态图| 国产成人精品福利久久| 日本欧美国产在线视频| 亚洲精品日韩在线中文字幕| 精品99又大又爽又粗少妇毛片| 在线观看三级黄色| av在线亚洲专区| 国产免费一区二区三区四区乱码| 亚洲av欧美aⅴ国产| 特大巨黑吊av在线直播| 亚洲精品色激情综合| 国产成人午夜福利电影在线观看| 一个人观看的视频www高清免费观看| 1000部很黄的大片| 麻豆精品久久久久久蜜桃| 免费播放大片免费观看视频在线观看| 日韩视频在线欧美| 久久久久久久精品精品| 欧美xxⅹ黑人| 亚洲精品中文字幕在线视频 | 真实男女啪啪啪动态图|