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

    山東昌樂地區(qū)王家大山玄武巖成因研究

    2016-06-23 06:11:06張銀慧任鐘元洪路兵黃小龍陳林麗
    地球化學(xué) 2016年1期
    關(guān)鍵詞:橄欖石輝石源區(qū)

    張銀慧, 任鐘元, 洪路兵, 張?樂, 黃小龍, 陳林麗

    ?

    山東昌樂地區(qū)王家大山玄武巖成因研究

    張銀慧1,2, 任鐘元1*, 洪路兵1,3, 張?樂1, 黃小龍1, 陳林麗1

    (1. 中國科學(xué)院 廣州地球化學(xué)研究所 同位素地球化學(xué)國家重點(diǎn)實(shí)驗(yàn)室, 廣東 廣州?510640, 2. 中國科學(xué)院大學(xué), 北京?100049, 3. 中國科學(xué)院 南海海洋研究所 邊緣海地質(zhì)重點(diǎn)實(shí)驗(yàn)室, 廣東 廣州?510301)

    報(bào)道了山東昌樂地區(qū)新生代王家大山玄武巖的全巖主元素、微量元素, 橄欖石斑晶及被其捕獲的熔體包裹體主元素和鉛同位素組成。王家大山玄武巖為弱堿性-拉斑玄武巖。全巖主元素成分和熔體包裹體的SiO2、TiO2、Al2O3、CaO、Na2O、K2O和P2O5含量均與MgO含量呈負(fù)相關(guān)關(guān)系, CaO/Al2O3值基本不變; 在微量元素原始地幔標(biāo)準(zhǔn)化蛛網(wǎng)圖上與EMI洋島玄武巖(OIB)相似, 表現(xiàn)為明顯的Ba、K和Sr正異常, Th和Pb負(fù)異常, 無Nb-Ta和Ti異常。熔體包裹體鉛同位素組成為208Pb/206Pb=2.091~2.166和207Pb/206Pb= 0.846~0.899, 其變化范圍覆蓋了山東地區(qū)新生代強(qiáng)堿性玄武巖到弱堿性-拉斑玄武巖的鉛同位素范圍; 根據(jù)鉛同位素組成, 熔體包裹體可以分為低鉛同位素組(207Pb/206Pb<0.86)和高Pb同位素組(207Pb/206Pb>0.87)。橄欖石斑晶Ni、Fe/Mn和Ca分別為1312~2417 μg/g、71~92和745~2068 μg/g, 與典型輝石巖熔體結(jié)晶的橄欖石斑晶成分相當(dāng), 但是, 比橄欖巖熔體結(jié)晶的橄欖石成分高Ni和Fe/Mn, 低Ca。全巖和熔體包裹體主元素的成分變化, 結(jié)合巖相學(xué)和MELTS軟件模擬結(jié)果, 指示王家大山玄武巖主要經(jīng)歷了橄欖石結(jié)晶分離作用。橄欖石斑晶高Ni和Fe/Mn, 低Ca, 結(jié)合全巖高Fe/MnO值和低CaO, 說明其源區(qū)巖石為輝石巖。全巖微量元素和熔體包裹體鉛同位素指示源區(qū)有EMI組分, 可能和再循環(huán)洋殼輝長巖有關(guān)。熔體包裹體兩組鉛同位素組成指示王家大山玄武巖源區(qū)是不均一的, 暗示山東地區(qū)新生代玄武巖的源區(qū)是高度不均一的。

    玄武巖; 橄欖石; 輝石巖; 熔體包裹體; 地幔不均一性; 山東省

    0?引?言

    分布于山東地區(qū)的新生代玄武巖, 是研究該地區(qū)深部地幔物質(zhì)組成, 熱狀態(tài)和巖石圈演化的重要巖石探針, 多年來一直備受關(guān)注[1?14]。已有的研究表明, 山東地區(qū)新生代玄武巖的源區(qū)巖石主要為輝石巖, 源區(qū)中有EMI組分[10,11,15]。然而, 關(guān)于源區(qū)EMI的成因, 卻存在著爭議:部分學(xué)者根據(jù)山東及其鄰區(qū)新生代玄武巖的Sr-Nd-Pb同位素具有大陸下地殼的屬性, 認(rèn)為再循環(huán)大陸下地殼是源區(qū)EMI組分的來源[10,16,17]; 另外一部分學(xué)者則對山東及其鄰區(qū)新生代玄武巖的橄欖石斑晶進(jìn)行了氧同位素研究, 發(fā)現(xiàn)異常低18O的橄欖石斑晶(18O<5.0‰; 正常地幔18O=(5.2±0.2)‰[18]), 結(jié)合總體虧損的Sr-Nd-Pb同位素和微量元素?zé)oNb-Ta等高場強(qiáng)元素虧損, 認(rèn)為源區(qū)的富集組分與再循環(huán)洋殼及其沉積物有關(guān), 其源區(qū)輝石巖是來自再循環(huán)洋殼及其沉積物的熔體交代巖石圈地幔的結(jié)果[11,19,20]。另外, 前人對山東地區(qū)新生代玄武巖的研究還表明山東地區(qū)新生代玄武巖的源區(qū)在大尺度上是非常不均一的: 強(qiáng)堿性玄武巖顯示HIMU-洋島玄武巖(OIB)的微量元素特征和總體虧損Sr-Nd-Pb同位素, 源區(qū)含有碳酸鹽; 弱堿性-拉斑玄武巖則具有EMI-OIB的微量元素特征和總體更富集Sr-Nd-Pb同位素[9?11,15]。然而, 這種現(xiàn)象是否也出現(xiàn)在單一的玄武質(zhì)巖漿的源區(qū)中, 還不是很清楚。值得提出的是, 目前國內(nèi)外對火山巖成因方面的研究主要是利用火山巖的全巖地球化學(xué)分析手段, 火山巖的全巖成分反映的不僅是原生巖漿的成分, 還反映了巖漿形成到噴發(fā)過程中的一系列巖漿后期作用的影響, 如不同來源巖漿的混合、結(jié)晶分異及同化混染等作用對巖漿成分的影響, 致使原生巖漿成分及其演化、源區(qū)母巖方面的重要信息等變得模糊或消失, 這給巖石學(xué)成因方面的研究帶來了較大的困難。熔體包裹體是巖漿結(jié)晶礦物捕獲的熔體液滴, 記錄捕獲時(shí)的熔體組分[21?22]。因此, 與全巖成分不同, 巖漿中較早結(jié)晶的橄欖石中的熔融包裹體能較好地記錄原生巖漿的成分及巖漿來源方面的重要信息, 而晚期結(jié)晶的橄欖石中的熔體包裹體則可以記錄巖漿演化過程的重要信息[23?27]。

    本研究擬以山東昌樂地區(qū)的王家大山玄武巖為研究目標(biāo)(圖1), 利用橄欖石及其被捕獲的熔體包裹體, 結(jié)合全巖組成, 探討其源區(qū)巖性, 并利用熔體包裹體成分揭示其地幔源區(qū)不均一性。

    1?地質(zhì)背景及樣品特征

    研究區(qū)位于山東省郯廬斷裂帶中段, 其北部被大面積的第四系沉積物覆蓋, 南部出露的地層有太古代泰山群片麻巖、混合巖和侵入巖基底; 震旦系石英巖和頁巖, 古生代寒武系饅頭組、張夏組頁巖和徐莊組的砂巖, 頁巖泥灰?guī)r和灰?guī)r等; 中生代侏羅系汶南組礫巖、石英砂巖, 白堊系青山組的雜色安山巖及王氏組砂頁巖、砂巖; 新生代古近系五圖組砂巖、頁巖和礫巖。研究區(qū)新生代巖漿活動(dòng)比較廣泛, 在沂水、臨朐以及昌樂地區(qū)均有大面積的中新世玄武巖出露, 它們與太古代, 古生代和中生代地層呈不整合接觸產(chǎn)出(圖1b)。王家大山玄武巖位于山東省昌樂市郊區(qū), 為早-中中新世玄武巖(圖1b)。玄武巖為黑色, 斑狀結(jié)構(gòu), 塊狀構(gòu)造。斑晶主要為橄欖石, 體積約占6%~8%, 偶見輝石斑晶, 未發(fā)現(xiàn)斜長石斑晶; 基質(zhì)為間粒結(jié)構(gòu), 主要是輝石、斜長石和橄欖石, 可見磁鐵礦富礦物。橄欖石斑晶多呈自形到半自形, 邊部輕微蝕變, 大小以0.5~2 mm為主。大部分橄欖石斑晶中含有熔體包裹體(圖2), 20~150 μm; 有些橄欖石有尖晶石包裹體。

    2?分析方法

    本研究分析所涉及的前處理工作, 全巖主元素、微量元素分析, 橄欖石斑晶和熔體包裹體電子探針(EPMA)分析以及熔體包裹體鉛同位素分析都在中國科學(xué)院廣州地球化學(xué)研究所同位素地球化學(xué)國家重點(diǎn)實(shí)驗(yàn)室完成。

    2.1?全巖樣品分析方法

    將巖石碎成小塊, 挑出最新鮮的, 用去離子水沖洗, 并超聲震蕩20 min, 反復(fù)3次后烘干, 磨成200目粉末, 用于主元素、微量元素分析。碎樣的剩余部分, 用于挑選橄欖石單礦物。

    稱取約1.3 g 200目的巖石粉末樣品, 在920 ℃下灼燒3 h之后, 制成玻璃餅, 然后上機(jī)測定, 詳細(xì)方法見文獻(xiàn)[29]。主元素分析的儀器是Rigaku ZSX-100e型X熒光光譜儀(XRF), 所有主元素的分析精度優(yōu)于5%。微量元素分析采用酸溶法, 分析儀器是PE Elan 6000型電感耦合等離子體質(zhì)譜儀(ICP-MS), 過程詳見[30], 大部分微量元素的分析精度優(yōu)于5%。

    圖1?王家大山地質(zhì)圖(底圖引自文獻(xiàn)[28])

    圖2?王家大山玄武巖橄欖石中的熔體包裹體照片

    (a) 單偏光下, 巖石薄片中未均一化的橄欖石熔體包裹體, 其表面比較渾濁, 可能是由于包裹體在捕獲后, 子礦物的結(jié)晶造成的。(b) 反光鏡下, 環(huán)氧樹脂靶中均一化之后的包裹體, 該包裹體由均一的玻璃和氣泡組成

    (a) A melt inclusion hosted by an olivine in thin sections without homogenization, plane light. Note that the melt inclusion is not clean, suggesting the presence of daughter minerals crystallized after Being trapped. (b) A melt inclusion hosted in an olivine after homogenization in epoxy resin target, with a rounded shape and a shrinkage void, reflected light

    2.2?橄欖石及熔體包裹體EPMA分析方法

    將挑選出來的橄欖石在通有混合氣體H2-CO2的加熱爐中加熱, 使熔體包裹體均一化, 氧逸度和加熱溫度分別為石英-鐵橄欖石-磁鐵礦(QFM)和1250 ℃, 詳細(xì)過程見[25]。將加熱過的橄欖石制靶, 拋磨至熔體包裹體露出表面。

    橄欖石和熔體包裹體的化學(xué)組成用JEOL JXA-8100 EPMA測定。橄欖石的分析條件為: 電壓20 kV, 束斑直徑2 μm, 電流強(qiáng)度3.0×10?7A, 分析方法見[31]。分析過程中, 選用漢諾壩橄欖巖包體的橄欖石作為監(jiān)控標(biāo)樣, 監(jiān)測儀器的穩(wěn)定性, 監(jiān)測標(biāo)樣的標(biāo)準(zhǔn)偏差(2,=30): SiO2、MgO、FeO、CaO、NiO和MnO分析精度分別優(yōu)于0.1%、0.24%、0.13%、0.004%、0.017%和0.005%。熔體包裹體的分析條件: 電壓15 kV, 束斑直徑3 μm, 電流2.0×10?8A, 分析方法與[32]相似, 分析過程中, 選用JB-2玄武巖玻璃標(biāo)樣作為監(jiān)控標(biāo)樣, 監(jiān)測儀器的穩(wěn)定性, 監(jiān)測標(biāo)樣的標(biāo)準(zhǔn)偏差(2,=30): SiO2、TiO2、Al2O3、FeO、MgO、MnO、CaO、Na2O、K2O和P2O5分析精度分別優(yōu)于0.1%、0.24%、0.13%、0.004% (Ca: 27 μg/g), 0.017% (Ni: 137 μg/g)和0.005% (Mn: 39 μg/g)。

    2.3?熔體包裹體鉛同位素分析方法

    用激光剝蝕-多接收等離子體質(zhì)譜儀(LA-MC- ICPMS)測定熔體包裹體的原位鉛同位素。分析條件: 能量80 mJ, 能量衰減值為25%, 束斑大小45 μm, 激光頻率3 Hz, 積分時(shí)間0.262 s, 混合氣體為 800 mL/min He和2 mL/min N2。測定過程中, 選用錐組合JET樣品錐(sample cone)+X截取錐(skimmer cone)增強(qiáng)208Pb的信號強(qiáng)度, 提高鉛同位素組成的精度, 具體方法見[33]。玻璃標(biāo)樣(BHVO-2G)的分析準(zhǔn)確度(2,=75)分別為:208Pb/206Pb優(yōu)于0.15%,207Pb/206Pb優(yōu)于–0.20%; 相對分析精度(2,=75):208Pb/206Pb優(yōu)于0.25%,207Pb/206Pb優(yōu)于0.29%。

    3?結(jié)?果

    3.1?全巖地球化學(xué)組成

    王家大山玄武巖的全巖主元素和微量元素化學(xué)組成見表1。SiO2和全堿(Na2O+K2O)含量分別為46.8%~47.7%和3.4%~4.6%, 在TAS圖解中顯示弱堿性-拉斑玄武巖(圖3)。王家大山玄武巖的MgO為7.7%~10.5%(表1), 隨著MgO降低, 其他氧化物(SiO2、TiO2、Al2O3、CaO、Na2O、K2O和P2O5)均增加, 而CaO/Al2O3值基本不變(圖4)。在微量元素蛛網(wǎng)圖上, 王家大山玄武巖的微量元素組成與典型的EMI-OIB(Gough)相似, 富集大離子親石元素, 具有明顯的Ba、K和Sr正異常, Th和Pb負(fù)異常, 無高場強(qiáng)元素(Nb-Ta、Zr-Hf)異常(圖5)。這些特征與下地殼存在明顯差別(圖5)。

    圖3?王家大山玄武巖的全巖和熔體包裹體TAS圖解(底圖引自文獻(xiàn)[34])

    圖中, 4、10和11分別代表樣品DSW09-4、DSW09-10和DSW09-11

    4, 10, and 11 represent samples DSW09-4, DSW09-10, and DSW09-11, respectively

    3.2?橄欖石斑晶成分

    本次工作分析了108顆橄欖石的化學(xué)組成, 分析結(jié)果見表2和圖6。橄欖石的Fo為80~85。橄欖石的Ca含量745~2068 μg/g, 高于典型的地幔捕擄晶(CaO<0.1%[40?43]), 說明它們均是從巖漿中結(jié)晶的; 橄欖石的Ni含量和Fe/Mn值分別為1312~2417 μg/g和71.2~91.8。王家大山玄武巖的橄欖石斑晶成分落在典型輝石巖熔體(如Koolau和EI Herrio)[38]結(jié)晶的橄欖石斑晶成分演化趨勢上; 在相同橄欖石Fo時(shí), 它們比典型的橄欖巖熔體(如MORB和古老的科馬提巖)[38]結(jié)晶的橄欖石斑晶高Ni、Fe/Mn和低Ca。

    3.3?熔體包裹體主元素

    熔體包裹體在巖漿房或巖漿通道中被橄欖石捕獲后, 可能經(jīng)歷結(jié)晶作用以及與寄主橄欖石發(fā)生Fe-Mg交換再平衡過程[44], 因此, 需要對熔體包裹體成分進(jìn)行恢復(fù)?;謴?fù)熔體包裹體組成可以用軟件PROGRAM FEO_EQ2 V.3.2[44]實(shí)現(xiàn); 恢復(fù)熔體包裹體組成的參數(shù)為: Ford.[45]的熔體與橄欖石平衡模型, 氧逸度為Fe2+/FeT=0.9 (FeT為全Fe), 以及FeOT=11.4%(高M(jìn)gO玄武巖(MgO>8)的全巖平均FeOT含量)。

    本次工作分析了3個(gè)王家大山玄武巖的112個(gè)熔體包裹體的主元素, 恢復(fù)的熔體包裹體組成見表3。熔體包裹體的MgO為6.9%~10%, 總體上比全巖略低, 但是比全巖具有更寬的范圍, 這可能與全巖中橄欖石的堆晶作用有關(guān), 或者是沒有找到被更高Fo橄欖石熔體包裹體。在TAS圖解中(圖3), 熔體包裹體的成分落在堿性-拉斑玄武巖的分界線上, 指示弱堿性-拉斑玄武質(zhì)成分; 與全巖相比, 熔體包裹體的全堿含量總體上偏高, 這與熔體包裹體普遍具有更高的Na2O有關(guān)(圖4f), 可能是全巖在地表發(fā)生輕微蝕變作用的結(jié)果。圖4中, 熔體包裹體的成分變化趨勢與全巖基本相似, MgO與其他氧化物(SiO2、TiO2、Al2O3、CaO、Na2O、K2O和P2O5)呈負(fù)相關(guān)關(guān)系, 而CaO/Al2O3值與MgO的關(guān)系基本不變。

    圖4?王家大山玄武巖全巖和熔體包裹體MgO與其他主元素相關(guān)圖

    圖中, 4、10和11分別代表樣品DSW09-4、DSW09-10和DSW09-11

    4, 10, and 11 represent samples DSW09-4, DSW09-10, and DSW09-11, respectively

    表1?王家大山玄武巖主元素(%)和微量元素(μg/g)分析結(jié)果

    (續(xù)表1)

    圖5?王家大山玄武巖微量元素原始地幔標(biāo)準(zhǔn)化蛛網(wǎng)圖

    圖中, 4、10和11分別代表樣品DSW09-4、DSW09-10和DSW09-11。原始地幔數(shù)據(jù)引自文獻(xiàn)[35]; EMI(Gough)數(shù)據(jù)引自文獻(xiàn)http://georoc.mpch-mainz.gwdg.de/georoc/Start.asp; 下地殼數(shù)據(jù)引自文獻(xiàn)[36]; 輝石巖數(shù)據(jù)引自文獻(xiàn)[37]。其中, 輝石巖包體的原巖為再循環(huán)下洋殼輝長巖[37]

    4, 10, and 11 represent samples DSW09-4, DSW09-10, and DSW09-11, respectively. Data source: Primitive mantle compositions?[35]; EMI (Gough)?http://georoc.mpch-mainz.gwdg.de/georoc/Start.asp; Lower crust?[36]; Pyroxenite?[37]. The protolith of the pyroxenite xenolith is the gabbro of recycled lower oceanic crust[37]

    3.4?熔體包裹體的原位鉛同位素成分

    本次工作共測定了56個(gè)熔體包裹體的原位鉛同位素組成, 結(jié)果見表4。熔體包裹體的208Pb/206Pb和207Pb/206Pb值變化分別為2.091~2.166和0.846~ 0.889(表4、圖7a), 它們落在EMI、EMII、DMM之間, 部分樣品的鉛同位素組成與EMI相當(dāng)(圖7a)。與山東其他地區(qū)新生代玄武巖的鉛同位素相比, 王家大山玄武巖的熔體包裹體鉛同位素組成覆蓋了強(qiáng)堿性和弱堿性-拉斑玄武巖的鉛同位素范圍, 但是比魯西地區(qū)周村玄武巖的鉛同位素組成低。在207Pb/206Pb-SiO2(MgO)圖解中, 熔體包裹體的鉛同位素組成可以劃分為兩組: 低207Pb/206Pb(<0.86)熔體包裹體和高207Pb/206Pb(>0.87)熔體包裹體, 落在中間的樣品很少(圖7b、圖7c)。這兩組熔體包裹體鉛同位素組成(207Pb/206Pb)不隨SiO2(MgO)變化而變化。

    表2?王家大山玄武巖的橄欖石電子探針分析結(jié)果

    (續(xù)表2)

    (續(xù)表2)

    圖6?王家大山玄武巖橄欖石斑晶的成分圖(底圖引自文獻(xiàn)[38])

    Koolau、Loihi、Indian MORB、EPR MORB和古老的科馬提巖(2.7Ga)的數(shù)據(jù)引自文獻(xiàn)[31]); Mangaia的數(shù)據(jù)引自文獻(xiàn)[38]; EI Hierro的數(shù)據(jù)引自文獻(xiàn)[39]

    The data source: Koolau, Loihi, Indian MORB, EPR MORB and old Komatiites (2.7Ga)?[31]; Mangaia?[38]; EI Hierro?[39]

    4?討?論

    4.1?地殼混染和結(jié)晶分離

    由于巖漿在通過巖漿通道上升過程或在巖漿房中, 可能會(huì)經(jīng)歷分離結(jié)晶和地殼混染, 因此, 需要對這些過程進(jìn)行評估。

    在微量元素蛛網(wǎng)圖上, 王家大山玄武巖顯示無Nb-Ta異常和Pb負(fù)異常(圖5), 這和大陸地殼普遍存在的Nb-Ta負(fù)異常和Pb正異常有明顯區(qū)別[36], 暗示王家大山玄武巖可能沒有經(jīng)歷明顯的地殼混染。王家大山玄武巖Ba/Nb值為11.3~12.4, 遠(yuǎn)低于地殼(圖8a); Nb/La和Nb/U值分別為1.27~1.32和39.7~47.3, 均落在EMI-OIB范圍[48]內(nèi), 且顯著高于地殼(圖8b、圖8c); 這些比值與MgO沒有明顯相關(guān)關(guān)系, 也支持王家大山玄武巖無明顯地殼混染的觀點(diǎn)。然而, 上述觀點(diǎn)似乎和Ce/Pb值相矛盾: Ce/Pb值為15.1~25.7, 盡管也落在EMI-OIB范圍[48]內(nèi), 但是, Ce/Pb與MgO弱的負(fù)相關(guān)關(guān)系(圖8d), 暗示部分樣品(低Ce/Pb值樣品)可能經(jīng)歷地殼混染; 這也和簡單混合模擬計(jì)算結(jié)果相吻合(低Ce/Pb值樣品剛好落在與下地殼的混合線上, 圖9)。由于Pb在地殼中的含量非常高, Ce/Pb值示蹤地殼混染比其他元素比值更敏感, 上述看似矛盾的現(xiàn)象, 可以解釋為: 王家大山玄武巖經(jīng)歷了輕微地殼混染, 但是, 由于樣品數(shù)量有限, 該過程除Ce/Pb外, 沒有在其他元素比值上表現(xiàn)出來。然而, 地殼混染模擬結(jié)果需要高達(dá)30%下地殼組分加入, 而不是少量地殼組分, 這可能與下地殼組分非常復(fù)雜, 而我們選取的是全球平均下地殼組分作為端元有關(guān)。為了減少地殼混染的影響, 在接下來的討論中, 將不再討論樣品DSW09-4、DSW09-10和DSW09-11。

    表3?王家大山玄武巖熔體包裹體電子探針分析結(jié)果(%)

    (續(xù)表3)

    (續(xù)表3)

    (續(xù)表3)

    注: Fo_Host代表熔體包裹體寄主橄欖石的Fo

    表4?王家大山玄武巖的熔體包裹體原位鉛同位素組成

    注: Stderr為標(biāo)準(zhǔn)誤差

    隨著SiO2增加(MgO降低), 王家大山玄武巖的熔體包裹體鉛同位素組成基本不發(fā)生變化(圖7b、圖7c), 說明被橄欖石捕獲的熔體包裹體經(jīng)歷的地殼混染程度不明顯。考慮到王家大山玄武巖部分樣品經(jīng)歷了輕微的地殼混染, 我們認(rèn)為可能是由于橄欖石結(jié)晶較早, 捕獲的熔體包裹體代表巖漿較早期的成分, 而地殼混染發(fā)生在巖漿晚期, 因而沒有被熔體包裹體所記錄。因此, 為了盡可能減低地殼混染的影響, 在接下來的討論中, 也將不再討論樣品DSW09-10的熔體包裹體。

    圖7?王家大山玄武巖熔體包裹體鉛同位素組成圖

    無棣大山玄武巖的數(shù)據(jù)引自文獻(xiàn)[2]; 強(qiáng)堿性和弱堿性-拉斑玄武巖數(shù)據(jù)分別引自文獻(xiàn)[2,11]和[11]; EMI、EMII、DMM和HIMU端元組分引自文獻(xiàn)[46]; 周村玄武巖熔體包裹體數(shù)據(jù)來自[47]。二元混合計(jì)算時(shí), 強(qiáng)堿性玄武巖端元來自無棣大山玄武巖(207Pb/206Pb=0.843和208Pb/206Pb=2.067[2]; 平均Pb=5.59 μg/g[9]); 弱堿性-拉斑玄武巖端元來自昌樂玄武巖(207Pb/206Pb=0.890和208Pb/206Pb=2.174; Pb= 2.3 μg/g)[11]

    Data source: Wudi – [2]; The strongly alkaline basalts (light dark area) and the weakly alkaline-theloiitic (dark grey area) basalts–[2, 11] and [11], respectively; mantle endmembers EMI, EMII, DMM and HIMU– [46]; The data of Pb isotope of melt inclusion for Zhoucun basalts–[47]. Note in mixing calculation we use data for basalts from Wudidashan (207Pb/206Pb=0.843,208Pb/206Pb=2.067[2]and average Pb=5.59 μg/g[9]) and from Changle (207Pb/206Pb=0.890 and208Pb/206Pb=2.174, Pb= 2.3 μg/g)[11]as the endmembers for strongly alkaline basalts and weakly alkaline-theloiitic basalts, respectively

    玄武巖在巖漿房或者巖漿通道中結(jié)晶的礦物主要為橄欖石、輝石和斜長石。薄片中, 王家大山玄武巖的斑晶主要為橄欖石, 偶見單斜輝石, 幾乎沒有斜長石斑晶, 初步說明王家大山玄武巖主要經(jīng)歷了橄欖石的分離結(jié)晶。玄武質(zhì)巖漿的主元素成分變化可以示蹤巖漿的結(jié)晶分離過程, 一般地, 輝石結(jié)晶導(dǎo)致殘留熔體SiO2、CaO和CaO/Al2O3值突然下降; 斜長石結(jié)晶導(dǎo)致殘留熔體Al2O3下降和CaO/Al2O3值快速增加; 而橄欖石結(jié)晶導(dǎo)致除MgO和FeO外的其他各氧化物緩慢增加, CaO/Al2O3值基本不變。在MgO與其他氧化物相關(guān)性圖解(圖4)中, 全巖和熔體包裹體的各氧化物均與MgO呈反相關(guān)關(guān)系, 而CaO/Al2O3不隨MgO變化(圖4), 說明王家大山玄武巖的主元素成分主要受控于橄欖石結(jié)晶分離作用。最后, 我們用MELTS軟件[49]模擬了王家大山玄武巖的結(jié)晶分離過程, 選取MgO含量最高的熔體包裹體DSW09-7(2)-5(MgO約10.0%, 表3)作為母巖漿, 氧逸度為QFM, 壓力0.1~0.5 GPa, 溫度為1400~1000 ℃, H2O=0% (嘗試了含水量0%~0.5%范圍的模擬, 由于水對熔體中只有橄欖石分離結(jié)晶期間的影響不大, 文中沒有給出)。模擬結(jié)果(圖10)表明, 王家大山玄武巖的成分和只經(jīng)歷橄欖石結(jié)晶分離趨勢最吻合。因此, 本文認(rèn)為王家大山玄武巖主要經(jīng)歷了橄欖石結(jié)晶分離過程。

    4.2?源區(qū)輝石巖及其形成機(jī)制

    最近的研究表明, 地幔中的橄欖巖和輝石巖, 都可以是玄武巖的源區(qū)巖石[25,31,50?54], 因此, 需要對玄武巖的源區(qū)巖性進(jìn)行甄別。到目前為止, 有效識(shí)別玄武巖源區(qū)巖性是橄欖巖還是輝石巖, 主要有三種指標(biāo): 全巖FeO/MnO值、全巖CaO和橄欖石組成。這是因?yàn)榈蒯V械牡V物主要為橄欖石、輝石和石榴子石, 其中, 橄欖石和輝石分別是Ni和Ca的主要載體, 同時(shí), 橄欖石的Fe-Mn分配系數(shù)(Fe/Mn)>1, 而輝石和石榴子石的Fe/Mn均小于1[55]。因此, 貧橄欖石的輝石巖, 與橄欖巖相比, 具有較低的Ni和Fe/Mn, 較高的Ca, 相應(yīng)地, 部分熔融形成的輝石巖熔體具有高Ni和高FeO/MnO值, 低CaO[31,53,54]。作為玄武巖巖漿中結(jié)晶最早的礦物, 橄欖石會(huì)記錄巖漿的這些特征, 即輝石巖熔體結(jié)晶的橄欖石, 與橄欖巖熔體相比, 具有更高的Ni和FeO/MnO值, 更低的CaO[31,52,54]。王家大山玄武巖全巖和熔體包裹體的CaO較低, 在CaO-MgO圖中, 落在輝石巖熔體的范圍內(nèi)(圖11a), 暗示其源區(qū)巖性為輝石巖。同時(shí), 王家大山玄武巖的FeO/MnO=69~74, 高于正常的橄欖巖熔體FeO/MnO值(小于60; 圖11b)[56], 進(jìn)一步指示王家大山玄武巖的源區(qū)可能為輝石巖。最后, 相對于典型的橄欖巖熔體結(jié)晶的橄欖石(如科馬提巖熔巖、MORB)[38], 王家大山玄武巖的橄欖石斑晶具有高Ni、Fe/Mn值, 低Ca含量, 它們與典型輝石巖熔體結(jié)晶的橄欖石成分相當(dāng)(圖6a、圖6b、圖6c)。盡管被硅過飽和熔體交代形成的再飽滿橄欖巖(refertilized peridotite)熔體, 如Mangaia和Loihi (圖6d、圖6e、圖6f), 也可以結(jié)晶高Ni橄欖石[38], 但是, 這些橄欖石一般具有高Ca和較低的Fe/Mn值 (圖6d、圖6e、圖6f), 它們和典型輝石巖熔體結(jié)晶的橄欖石有明顯差別[38]。因此, 本文認(rèn)為王家大山玄武巖的源區(qū)巖石為輝石巖。

    圖8?王家大山玄武巖的Ba/Nb、Nb/La和Nb/U與MgO關(guān)系圖解

    圖中, 4, 10和11分別指樣品DSW09-4, DSW09-10和DSW09-11。上地殼、下地殼和地殼的成分引自文獻(xiàn)[36], EMI-OIB范圍引自文獻(xiàn)[48]

    4, 10, and 11 represent samples DSW09-4, DSW09-10, and DSW09-11, respectively. The compositions of upper crust, lower crust and bulk crust are from [36], The ranges of EMI-OIBs are after [48]

    圖9?王家大山玄武巖地殼混染模擬圖解

    圖中, 4, 10, 和11分別指樣品DSW09-4, DSW09-10, 和DSW09-11。地殼數(shù)據(jù)(上地殼, 下地殼和全地殼)引自文獻(xiàn)[36]

    4, 10, and 11 represent samples DSW09-4, DSW09-10, and DSW09-11, respectively.Crustdataforthelower,upperandbulkcrustare from [36]

    滿足玄武質(zhì)熔體結(jié)晶高Ni橄欖石的源區(qū)輝石巖必然是由橄欖巖經(jīng)歷富硅熔體交代反應(yīng)形成的[31,52,54,58]。這類輝石巖有兩種方式:在島弧地區(qū), 俯沖板片釋放富SiO2熔體/流體, 它們與地幔楔橄欖巖反應(yīng), 形成富斜方輝石(Opx)輝石巖[58]。由于俯沖板片釋放的富SiO2熔體/流體虧損Nb-Ta等高場強(qiáng)元素, 反應(yīng)形成的富Opx輝石巖以及由它為源區(qū), 熔融形成的玄武質(zhì)熔體均應(yīng)當(dāng)虧損Nb-Ta等高場強(qiáng)元素[58]。然而, 王家大山玄武巖的全巖微量元素并未顯示Nb-Ta等高場強(qiáng)元素虧損特征(圖5), 說明該方式不適合王家大山玄武巖的源區(qū)輝石巖, 可以排除。我們傾向于第二種方式:對流地幔中的再循環(huán)地殼組分, 它們可能是經(jīng)歷俯沖脫水的洋殼(再循環(huán)洋殼), 也可能是在洋殼俯沖過程中帶入地?;蛘叽箨憥r石圈拆沉進(jìn)入地幔的大陸下地殼(再循環(huán)大陸下地殼), 被上升地幔流(如地幔柱)攜帶至地幔淺部的過程中, 由于它們具有比橄欖巖更低的固相線, 優(yōu)先發(fā)生熔融形成SiO2過飽和熔體, 并快速與周圍橄欖巖反應(yīng), 生成輝石巖[50?52]。

    圖10?MELTS分離結(jié)晶模擬結(jié)果與王家大山玄武巖和熔體包裹體成分對比圖

    粗虛線和粗實(shí)線分別代表斜長石(Pl)和單斜輝石(Cpx)開始結(jié)晶

    Note that the thick dotted and solid lines represent the beginning of crystallization of plagioclase (Pl) and clinopyroxene (Cpx), respectively

    圖11?王家大山玄武巖全巖和熔體包裹體的Ca和FeO/MnO與MgO關(guān)系圖解

    圖中FeO為全鐵。圖11a底圖引自文獻(xiàn)[56]。圖1s1b橄欖巖部分熔融熔體的成分引自文獻(xiàn)[54]。橄欖巖熔體和輝石巖熔體的分界線FeO/MnO = 60(虛線)引自文獻(xiàn)[57]

    FeO data is normalized total Fe. (a) The range of CaO contents for the peridotite partial melts is after [56]. (b) FeO/MnO ratios of the peridotite partial melts are after [54]. The boundary of FeO/MnO=60 between peridotite- and pyroxenite- partial melts (the dashed line) is after [57]

    4.3?源區(qū)EMI組分

    山東地區(qū)新生代玄武巖的微量元素和同位素組成顯示源區(qū)存在EMI組分[10?11]。然而, 關(guān)于EMI組分來自再循環(huán)大陸下地殼還是再循環(huán)洋殼, 卻存在爭議[10?11]。王家大山玄武巖的全巖微量元素和熔體包裹體鉛同位素組成均顯示源區(qū)中有EMI組分的貢獻(xiàn)(圖5、圖7a), 因而可以為山東地區(qū)新生代玄武巖的源區(qū)EMI組分提供約束??偟膩碚f, 地幔EMI組分有三種可能成因: 再循環(huán)大陸下地殼, 再循環(huán)遠(yuǎn)洋沉積物和再循環(huán)洋殼輝長巖[59]。盡管再循環(huán)大陸下地殼也具有Ba、K和Sr正異常, 但是, 其明顯的Nb-Ta負(fù)異常和Pb正異常有別于王家大山玄武巖(圖5), 可以排除。橄欖石斑晶的氧同位素可以很好區(qū)分EMI組分是再循環(huán)遠(yuǎn)洋沉積物還是洋殼輝石巖。一般地, 遠(yuǎn)洋沉積物由于在洋殼表面經(jīng)歷低溫?zé)嵋何g變作用, 往往具有比正常地幔高的18O值, 而洋殼輝石巖則在洋殼深部經(jīng)歷高溫?zé)嵋何g變作用, 具有比正常地幔低的18O值[18,60]。盡管本文沒有分析王家大山玄武巖橄欖石斑晶的氧同位素組成, 但是, Xu.[11]報(bào)道了研究區(qū)內(nèi)同期的昌樂玄武巖的橄欖石斑晶具有異常低的氧同位素組成(18O低至4.09‰); 另外, 王家大山玄武巖全巖的微量元素組成蛛網(wǎng)圖分布模式基本和典型的再循環(huán)洋殼輝長巖非常相似(圖5), 這些特征都和洋殼輝長巖吻合。因此, 本文認(rèn)為王家大山玄武巖的源區(qū)EMI組分和含輝長巖的再循環(huán)洋殼有關(guān)。

    4.4?源區(qū)不均一性

    已有的研究表明, 山東地區(qū)新生代玄武巖的源區(qū)是非常不均一的[9?11,15], 這種現(xiàn)象也出現(xiàn)在王家大山玄武巖中。王家大山玄武巖的熔體包裹體鉛同位素組成從強(qiáng)堿性玄武巖區(qū)域逐漸變化到弱堿性-拉斑玄武巖區(qū)域(圖7a); 同時(shí), 熔體包裹體鉛同位素組成顯示王家大山玄武巖存在兩組熔體包裹體(低207Pb/206Pb和高207Pb/206Pb, 圖7b、圖7c), 這些特征說明王家大山玄武巖的源區(qū)是非常不均一的, 存在堿性玄武巖的源區(qū)端元和弱堿性-拉斑玄武巖源區(qū)端元。鉛同位素的二元混合計(jì)算結(jié)果表明 (圖7a), 王家大山玄武巖高鉛同位素熔體包裹體以弱堿性-拉斑玄武巖源區(qū)端元為主, 強(qiáng)堿性玄武巖源區(qū)源區(qū)端元小于20%, 而低鉛同位素熔體包裹體含更高比例的強(qiáng)堿性玄武巖源區(qū)端元, 強(qiáng)堿性玄武巖源區(qū)端元可達(dá)60%。

    然而, 必須注意到王家大山玄武巖的熔體包裹體組分均為弱堿性-拉斑玄武質(zhì), 沒有強(qiáng)堿性玄武質(zhì)成分, 這似乎和源區(qū)中有強(qiáng)堿性玄武巖源區(qū)端元相矛盾。這種現(xiàn)象可以用Ol-CaTs-Qz相圖來解釋。含CO2橄欖巖和輝石巖的熔融實(shí)驗(yàn)結(jié)果表明: 低程度熔融形成強(qiáng)堿性熔體, 落在Ol-CaTs的左側(cè); 隨著熔融程度增加, 熔體堿性減弱, 在相圖上逐漸往右側(cè)移動(dòng); 最終落在Ol-CaTs的右側(cè)(圖12a)。山東地區(qū)強(qiáng)堿性玄武巖的源區(qū)有CO2[9,15], 它們落在Ol-CaTs的左側(cè)(圖12b), 說明可能經(jīng)歷低程度的部分熔融; 所有王家大山玄武巖的熔體包裹體組分落在Ol-CaTs右側(cè)(圖12b), 可能經(jīng)歷了更高程度的部分熔融。因此, 本文認(rèn)為王家大山玄武巖中強(qiáng)堿性玄武巖源區(qū)產(chǎn)生弱堿性-拉斑玄武質(zhì)熔體可能是高程度部分熔融的結(jié)果。

    圖12?從透輝石[Di]投影的橄欖石(Ol)-Ca-Tschermaks(CaTs)-石英(Qz)相圖(底圖引自文獻(xiàn)[53])

    (a)實(shí)驗(yàn)熔體; (b) 山東地區(qū)新生代玄武巖和王家大山玄武巖熔體包裹體(本研究)。山東地區(qū)新生代玄武巖均為只經(jīng)歷橄欖石結(jié)晶分離作用, 并通過加減平衡橄欖石校正到與Fo=90橄欖石平衡。實(shí)驗(yàn)熔體數(shù)據(jù)來源:橄欖巖+CO2引自文獻(xiàn)[61?62]; 硅不飽和輝石巖+CO2引自文獻(xiàn)[63]; MORB-like輝石巖+CO2引自文獻(xiàn)[64?65]。山東地區(qū)新生代玄武巖強(qiáng)堿性玄武巖和弱堿性-拉斑玄武巖的數(shù)據(jù)分別引自文獻(xiàn)[9,57,66,67]和[10,11,66]

    (a): high pressure experiments; (b) the compositions of Cenozoic basalts in Shandong Province and melt inclusions of Wangjiadashan basalts (this study). All Cenozoic basalts in Shandong Province are olivine-controlled, and corrected to in equilibration witholivine Fo=90 by adding/subtracting the equilibrated olivine. Data source for experimental partial melts: peridotite+CO2-[61?62]; SiO2-undersaturated pyroxenite+CO2-[63]; and SiO2-oversaturated pyroxenite +CO2-[64?65]. Data source for the strongly alkaline basalts and weakly alkaline-theloiitic basalts of Cenozoic basalts in Shandong is from [9,57,66,67] and [10,11,66], respectively

    5?結(jié)?論

    (1) 山東昌樂地區(qū)王家大山玄武巖為弱堿性-拉斑玄武巖, 主要經(jīng)歷橄欖石結(jié)晶分離作用。

    (2) 根據(jù)橄欖石斑晶成分, 結(jié)合全巖FeO/MnO值和CaO, 指示源區(qū)巖性為輝石巖。

    (3) 王家大山玄武巖的全巖微量元素和熔體包裹體鉛同位素組成, 指示源區(qū)有EMI組分, 它可能和再循環(huán)洋殼輝長巖有關(guān)。

    (4) 王家大山玄武巖的熔體包裹體鉛同位素組成指示其源區(qū)是不均一的。

    感謝中國科學(xué)院廣州地球?qū)W研究所劉穎工程師、涂湘林工程師和胡光黔工程師以及吳蕾研究員助理在全巖主元素以及熔體包裹體樣品制作和電子探針分析工作中的幫助。感謝兩位匿名審稿人提出的寶貴意見。

    [1] Zhou X H, Armstrong R L. Cenozoic volcanic rocks of eastern China-secular and geographic trends in chemistry and strontium isotopic composition[J]. Earth Planet Sci Lett, 1982, 58(3): 301–329.

    [2] Peng Z C, Zartman R E, Futa K, Chen D G. Pb-, Sr- and Nd-isotopic systematics and chemical characteristics of Cenozoic basalts, Eastern China[J]. Chem Geol Isot Geosci, 1986, 59(86): 3–33.

    [3] Zhi X C, Song Y, Frey F A, Fen J L, Zhai M Z. Geochemistry of Hannuoba basalts, eastern China: Constraints on the origin of continental alkalic and tholeiitic basalt[J]. Chem Geol, 1990, 88(1/2): 1–33.

    [4] Basu A R, Wang J W, Huang W K, Xie G H, Tatsumoto M. Major element, REE, and Pb, Nd and Sr isotopic geochemistry of Cenozoic volcanic rocks of eastern China: Implications for their origin from suboceanic-type mantle reservoirs[J]. Earth Planet Sci Lett, 1991, 105(1–3): 149–169.

    [5] Fan Q C, Hooper P R. The Cenozoic basaltic rocks of eastern China: Petrology and chemical composition[J]. J Petrol, 1991, 32(4): 765–810.

    [6] Liu C Q, Masuda A, Xie G H. Major element compositions and rare-earth element abundances of Cenozoic basalts in Eastern China: Implications for a pressure control over LREE/HREE fractionation in continental basalt [J]. Chinese J Geochem, 1992, 11(4): 289–313.

    [7] Liu C Q, Masuda A, Xie G H. Major-and trace-element compositions of Cenozoic basalts in eastern China: Petrogenesis and mantle source[J]. Chem Geol, 1994, 114(94): 19–42.

    [8] Chung S L. Trace element and isotope characteristics of Cenozoic basalts around the Tanlu fault with implications for the eastern plate boundary between North and South China[J]. J Geol, 1999, 107(3): 301–312.

    [9] Zeng G, Chen L H, Xu X S, Jiang S Y, Hofmann A W. Carbonated mantle sources for Cenozoic intra-plate alkaline basalts in Shandong, North China[J]. Chem Geol, 2010, 273(1/2): 35–45.

    [10] Zeng G, Chen L H, Hofmann A W, Jiang S Y, Xu X S. Crust recycling in the sources of two parallel volcanic chains in Shandong, North China[J]. Earth Planet Sci Lett, 2011, 302(3): 359–368.

    [11] Xu Z, Zhao Z F, Zheng Y F. Slab-mantle interaction for thinning of cratonic lithospheric mantle in North China: Geochemical evidence from Cenozoic continental basalts in central Shandong[J]. Lithos, 2012, 146–147: 202–217.

    [12] Xu Y G. Recycled oceanic crust in the source of 90–40 Ma basalts in North and Northeast China: Evidence, provenance and signi?cance[J]. Geochim Cosmochim Acta, 2014, 143: 49–67.

    [13] 王慧芬, 朱炳泉, 張前鋒, 范彩云, 董麗敏. 山東臨朐地區(qū)新生代玄武巖同位素鉀-氬年齡研究[J].地球化學(xué), 1981, 10(4): 321–328. Wang Hui-fen, Zhu Bing-quan, Zhang Qian-feng, Fan Cai-yun, Dong Li-min. A study on K-Ar isotopic ages of Cenozoic basalts from Linqu area, Shandong Province[J]. Geochimica, 1981, 10(4): 321–328 (in Chinese with English abstract).

    [14] 陳道公, 彭子成. 山東新生代火山巖K-Ar年齡和Pb-Sr同位素特征[J]. 地球化學(xué), 1985, 14(4): 293–303 Chen Dao-gong, Peng Zi-cheng. K-Ar ages and Pb, Sr isotopic characteristics of Cenozoic volcanic rocks in Shandong, China[J].Geochimica, 1985, 14(4): 293–303 (in Chinese with English abstract).

    [15] Sakuyama T, Tian W, Kimura J I, Fukao Y, Hirahara Y, Takahashi T, Senda R, Chang Q, Miyazaki T, Obayashi M, Kawabata H, Tatsumi Y. Melting of dehydrated oceanic crust from the stagnant slab and of the hydrated mantle transition zone: Constraints from Cenozoic alkaline basalts in eastern China[J]. Chem Geol, 2013, 359: 32–48.

    [16] Chen L H, Zeng G, Jiang S Y, Hofmann A W, Xu X S. Sources of Anfengshan basalts: Subducted lower crust in the Sulu UHP belt, China[J]. Earth Planet Sci Lett, 2009, 286(3/4): 426–435.

    [17] 陳立輝, 曾罡, 胡森林, 俞恂, 陳霞玉. 地殼再循環(huán)與大陸堿性玄武巖的成因: 以山東新生代堿性玄武巖為例[J]. 高校地質(zhì)學(xué)報(bào), 2012, 18(1): 16–27.Chen Li-hui, Zeng Gang, Hu Sen-lin, Yu Xun, Chen Xia-yu. Crustal recycling and genesis of continental alkaline basalts: Case study of the Cenozoic alkaline basalts from Shandong Province, eastern China[J]. Geol J China Univ, 2012, 18(1): 16–27 (in Chinese with English abstract).

    [18] Eiler J M, Farley K A, Valley J W, Hofmann A W, Stolper E M. Oxygen isotope constraints on the sources of Hawaiian volcanism[J]. Earth Planet Sci Lett, 1996, 144(3/4): 453–468.

    [19] Zhang J J, Zheng Y F, Zhao Z F. Geochemical evidence for interaction between oceanic crust and lithospheric mantle in the origin of Cenozoic continental basalts in east-central China[J]. Lithos, 2009, 110: 305–326.

    [20] Wang Y, Zhao Z F, Zheng Y F, Zhang J J. Geochemical constraints on the nature of mantle source for Cenozoic continental basalts in east-central China[J]. Lithos, 2011, 125: 940–955.

    [21] Roedder E. Fluid Inclusions, Reviews in Mineralogy and Geochemistry[M]. Chantilly: Mineralogical Society of America, 1984: 644p.

    [22] Frezzotti M L. Silicate-melt inclusions in magmatic rocks: Applications to petrology[J]. Lithos, 2001, 55(1–4): 273–299.

    [23] Sobolev A V. Melt inclusions in minerals as a source of principle petrological information[J]. Petrology, 1996, 4(3): 209–220.

    [24] Sobolev A V, Hofmann A W, Nikogosian I K. Recycled oceanic crust observed in ‘ghost plagioclase’ within the source of Mauna Loa lavas[J]. Nature, 2000, 404(6781): 986–989.

    [25] Ren Z Y, Ingle S, Takahashi E, Hirano N, Hirata T. The chemical structure of the Hawaiian mantle plume[J]. Nature, 2005, 436(7052): 837–840.

    [26] Hong L B, Zhang Y H, Qian S P, Liu J Q, Ren Z Y, Xu Y G. Constraints from melt inclusions and their host olivines on the petrogenesis of Oligocene-Early Miocene Xindian basalts, Chifeng area, North China Craton[J]. Contrib Mineral Petrol, 2013, 165(2): 305–326.

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

    [28] 山東省地質(zhì)礦產(chǎn)局. 山東省區(qū)域地質(zhì)志(中華人民共和國地質(zhì)礦產(chǎn)部地質(zhì)專報(bào)—區(qū)域地質(zhì)第26號)[M]. 北京: 地質(zhì)出版社, 1991: 1–595. Bureau of Geology and Mineral Resources of Shandong Province. Regional Geology of Shandong Province[M]. Beijing: Geological Publishing House, 1991: 1–595 (in Chinese with English abstract).

    [29] Goto A, Tatsumi Y. Quantitative analysis of rock samples by an X-ray fluorescence spectrometer (II)[J]. The Rigaku Journal, 1996, 13(2): 20–39.

    [30] 劉穎, 劉海臣, 李獻(xiàn)華. 用ICP-MS準(zhǔn)確測定巖石樣品中的40余種微量元素[J]. 地球化學(xué), 1996, 25(6): 552–558.Liu Ying, Liu Hai-chen, Li Xian-hua. Simultaneous and precise determination of 40 trace elements in rock samples using ICP-MS[J]. Geochimica, 1996, 25(6): 552–558 (in Chinese with English abstract).

    [31] Sobolev A V, Hofmann A W, Kuzmin D V, Yaxley G M, Arndt N T, Chung S L, Danyushevsky L V, Elliott E, Frey F A, Garcia M O, Gurenko A S, Kamenetsky V S, Kerr A C, Krivolutskaya N A, Matvienkov V V, Nikogosian I K, Rocholl A, Sigurdsson I A, Sushchevskaya N M, Teklay M. The amount of recycled crust in sources of mantle-derived melts[J]. Science, 2007, 316(5823): 412–417.

    [32] Wang Z, Gaetani G A. Partitioning of Ni between olivine and siliceous eclogite partial melt: Experimental constraints on the mantle source of Hawaiian basalts[J]. Contrib Mineral Petrol, 2008, 156: 661–678.

    [33] Zhang L, Ren Z Y, Nichols A R L, Zhang Y H, Zhang Y, Qian SP, Liu J Q. Lead isotope analysis of melt inclusions by LA-MC-ICP-MS[J]. J Anal At Spectrom, 2014, 29(8): 1393–1405.

    [34] Le Bas M J, Le Maitre R W, Strekeisen A L, Zanettin B A. A chemical classification of volcanic rocks based on the total alkali-silica diagram[J]. J Petrol, 1986, 27(3): 745–750.

    [35] McDonough W F, Sun S-s. The composition of the Earth[J]. Chem Geol, 1995, 120(3/4): 223–253.

    [36] Rudnick R L, Gao S. The composition of the continental crust[M]//Holland H D, Turekian K K. Treatise on Geochemistry. Oxford: Elsevier-Pergamon, 2003, 3: 1–64.

    [37] Yu S Y, Xu Y G, Ma J L, Zheng Y F, Kuang Y S, Hong L B, Ge W C, Tong L X. Remnants of oceanic lower crust in the subcontinental lithospheric mantle: Trace element and Sr-Nd-O isotope evidence from aluminous garnet pyroxenite xenoliths from Jiaohe, Northeast China[J].Earth Planet Sci Lett, 2010, 297(3/4): 413–422.

    [38] Herzberg C, Cabral R A, Jackson M G, Vidito C, Day J M D, Hauri E H. Phantom Archean crust in mangaia hotspot lavas and the meaning of heterogenous mantle[J]. Earth Planet Sci Lett, 2014, 396: 97–106.

    [39] Gurenko A A, Sobolev A V, Hoernle K A, Hauff F, Schmincke H-U. Enriched, HIMU-type peridotite and depleted recycled pyroxenite in the Canary plume: A mixed-up mantle[J]. Earth Planet Sci Lett, 2009, 277(3/4): 514–524.

    [40] Ramsay W R H, Crawford A J, Foden J D. Field setting, mineralogy, chemistry, and genesis of arc picrites, New-Georgia, Solomon-Islands[J]. Contrib Mineral Petrol, 1984, 88(4): 386–402.

    [41] Cameron W E. Petrology and origin of primitive lavas from the troodos ophiolite, Cyprus[J]. Contrib Mineral Petrol, 1985, 89(2): 239–255.

    [42] Boudier F. Olivine xenocrysts in picritic magmas-an experimental and microstructural study[J]. Contrib Mineral Petrology, 1991, 109(1): 114–123.

    [43] Rohrbach A, Schuth S, Ballhaus C, Münker C, Matveev S, Qopoto C. Petrological constraints on the origin of arc picrites, New Georgia Group, Solomon Islands[J]. Contrib Mineral Petrol, 2005, 149(6): 685–698.

    [44] Danyushevsky L V, Della-Pasqua F N, Sokolov S. Re-equilibration of melt inclusions trapped by magnesian olivine phenocrysts from subduction-related magmas: Petrological implication[J]. Contrib Mineral Petrol, 2000, 138(1): 68–83.

    [45] Ford C E, Russel D G, Craven J A, Fisk M R. Olivine-liquid equilibria: temperature, pressure and composition dependence of the crystal/liquid cation partition cofficients for Mg, Fe2+, Ca and Mn[J]. J Petrol, 1983, 24(3): 256–265.

    [46] Yurimoto H, Kogiso T, Abe K, Barsczus H G, Utsunomiya A, Maryuama S. Lead isotopic compositions in olivine-hosted melt inclusions from HIMU basalts and possible link to sulfide components[J]. Phys Earth Planet In, 2004, 146(1/2): 231–242.

    [47] 張銀慧, 任鐘元, 洪路兵, 張樂, 錢生平, 吳亞東, 丁相禮, 徐義剛, 陳林麗. 魯西地區(qū)新生代周村玄武巖的成因研究[J]. 大地構(gòu)造與成礦學(xué), 2015. (已接收)Zhang Yin-hui, Ren Zhong-yuan, Hong Lu-bing, Zhang Le, Qian Sheng-ping, Wu Ya-dong, Ding Xiang-li, Xu Yi-gang, Chen Lin-li. Petrogenesis of Cenozoic Zhoucun basalts, Western Shandong[J]. Geotectont Metallogen, 2015. (in press)

    [48] Willbold M, Stracke A. Trace element composition of mantle end-members: Implications for recycling of oceanic and upper and lower continental crust[J]. Geochem Geophys Geosyst, 2006, 7. doi: 10.1029/2005GC001005.

    [49] Ghiorso M S, Sack R O. Chemical mass transfer in magmatic processes. IV. A revised and internally consistent thermodynamic model for the interpolation and extrapolation of liquid-solid equilibria in magmatic systems at elevated and pressures[J]. Contrib Mineral Petrol, 1995, 119: 197–212.

    [50] Ren Z Y, Shibata T, Yoshikawa M, Johnson K T M, Takahashi E. Isotope compositions of the submarine Hana Ridge lavas, Haleakala volcano, Hawaii: Implications for source compositions, melting process and the structure of Hawaiian plume[J]. J Petrol, 2006, 47(2): 255–275.

    [51] Ren Z Y, Hanyu T, Miyazaki T, Change Q, Kawabata H, Takahashi T, Hirahara Y, Nichols A R L, Tatsumi Y. Geochemical differences of the Hawaiian shield lavas: Implications for melting process in the heterogeneous Hawaiian plume[J]. J Petrol, 2009, 50(8): 1553–1573.

    [52] Sobolev A V, Hofmann A W, Sobolev S V, Nikogosian I K. An olivine-free mantle source of Hawaiian shield basalts[J]. Nature, 2005, 434(7033): 590–597.

    [53] Herzberg C. Petrology and thermal structure of the Hawaiian plume from Mauna Kea volcano[J]. Nature, 2006, 444(7119): 605–609.

    [54] Herzberg C. Identification of source lithology in the Hawaiian and Canary Islands: Implications for origins[J]. J Petrol, 2011, 52(1): 113–146.

    [55] Davis F A, Humayun M, Hirschmann M M, Cooper R S. Experimentally determined mineral/melt partitioning of first-row transition elements (FRTE) during partial melting of peridotite at 3 GPa[J]. Geochim Cosmochim Acta, 2013, 104: 232–260.

    [56] Herzberg C, Asimow P D. Petrology of some oceanic island basalts: PRIMELT2.XLS software for primary magma calculation[J]. Geochem Geophys Geosyst, 2008, 9. doi: 10.1029/2008GC002057.

    [57] Liu Y S, Gao S, Kelemen P B, Xu W L. Recycled crust controls contrasting source compositions of Mesozoic and Cenozoic basalts in the North China Craton[J]. Geochim Cosmochim Acta, 2008, 72(9): 2349–2376.

    [58] Straub S M, LaGatta A B, Pozzo A L M D, Langmuir C H. Evidence from high-Ni olivines for a hybridized peridotite/pyroxenite source for orogenic andesites from the central Mexican Volcanic Belt[J]. Geochem Geophys Geosyst, 2008, 9(3). doi:10.1029/2007GC001583

    [59] Hofmann A W. Sampling mantle heterogeneity through oceanic basalts: Isotopes and trace elements[M]//Holland H D, Turekian K K. Treatise on Geochemistry. Oxford: Elsevier-Pergamon, 2003, 2: 61–101.

    [60] Eiler J M. Oxygen Isotope Variations of basatic lavas and upper manlt rocks[J]. Rev Mineral Geochem, 2001, 43(1): 319–364.

    [61] Hirose K. Partial melt compositions of carbonated peridotite at 3 GPa and role of CO2in alkali-basalt magma generation[J]. Geophys Res Lett, 1997, 24(22): 2837–2840.

    [62] Dasgupta R, Hirschmann M M, Smith N D. Partial melting experiments of peridotite+CO2at 3 GPa and genesis of alkalic ocean island basalts[J]. J Petrol, 2007, 48(11): 2093–2124.

    [63] Dasgupta R, Hirschmann M M, Stalker K. Immiscible transition from carbonate-rich to silicate-rich melts in the 3 Gpa melting interval of eclogite+CO2and genesis of silica-undersaturaed ocean island lavas[J]. J Petrol, 2006, 47(4): 647–671.

    [64] Gerbode C, Dasgupta R. Carbonate-fluxed melting of MORB-like pyroxenite at 2.9 GPa and genesis of HIMU ocean island basalts[J]. J Petrol, 2010, 51(10): 2067–2088.

    [65] Kiseeva E S, Yaxley G M, Hermann J, Litasov K D, Rosenthal A, Kamentesky V S. An experimental study of carbonated eclogite at 3.5–5.5 GPa: Implications for silicate and carbonate metasomatism in the cratonic mantle[J]. J Petrol, 2012, 53(4): 727–759.

    [66] Chen Y, Zhang Y, Graham D, Su S, Deng J. Geochemistry of Cenozoic basalts and mantle xenoliths in Northeast China[J]. Lithos, 2007, 96(1/2): 108–126.

    [67] 羅丹, 陳立輝, 曾罡. 陸內(nèi)強(qiáng)堿性火山巖的成因: 以山東無棣大山霞石巖為例[J]. 巖石學(xué)報(bào), 2009, 25(2): 311–319.Luo Dan, Chen Li-hui, Zeng Gang. Gensis of intra-continental strongly alkaline volcanic rocks: A case study of Dashan nephelinites in Wudi, Shandong Province, North China[J]. Acta Petrol Sinica, 2009, 25(2): 311–319 (in Chinese with English abstract).

    Petrogenesis of Cenozoic Wangjiadashan basalts in Changle area, Shandong Province

    ZHANG Yin-hui1,2, REN Zhong-yuan1*, HONG Lu-bing1,3, ZHANG Le1,HUANG Xiao-long1and CHEN Lin-li1

    1. State Key Laboratory of Isotope Geochemistry, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou?510640, China; 2. University of Chinese Academy of Sciences, Beijing?100049, China; 3. Key Laboratory of Marginal Sea Geology, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou?510301, China

    This study reported major and trace element compositions of bulk rocks, major element andPb isotope compositions of melt inclusions and chemistry of their host olivine in Wangjiadashan basalts, Changle area, Shandong Province to constrain their petrogenesis. The results show that the Wangjiadashan basalts are weakly alkaline to tholeiitic. The major element compositions of the bulk rocks of the Wangjiadashan basalts are similar to those of the melt inclusions there, both showing negative correlation between MgO and other oxides, such as SiO2, TiO2, Al2O3, CaO, Na2O, K2O and P2O5, and relative constant CaO/Al2O3ratios. However, the trace element compositions exhibit EMI-oceanic island basalts’ flavor, characterized by significant positive anomalies of Ba, K and Sr, negative anomalies of Th and Pb, but no anomalies of Nb-Ta and Ti. Pb isotope compositions of the melt inclusions in the Wangjiadashan basalts range from 2.091 ? 2.166 for208Pb/206Pb and 0.846 ? 0.899 for207Pb/206Pb, overlapping the Pb isotope range for the strongly alkaline basalts and weakly alkaline-tholeiitic basalts in Shandong Province. On the basis of Pb isotope compositions, melt inclusions can be classified into two groups: low and high Pb isotope group (207Pb/206Pb<0.86 and207Pb/206Pb>0.87, respectively). Olivine phenocrysts in the Wangjiadashan basalts show Ni 1312 ? 2417 μg/g, Fe/Mn ratios 71 ? 92, Ca 745 ? 2068 μg/g, which are comparable to those from the typical pyroxenitic melts. In contrast to the peridotitic melts-derived olivines, they generally exhibit higher Ni and Fe/Mn, lower Ca at given Fo values. The major element compositions of the bulk rocks and the melt inclusions, together with the petrography and MELTS modeling results, suggest that the Wangjiadashan basalts mainly underwent olivine-controlled fractionation. High Ni and Fe/Mn, low Ca of olivine compositions, combined with low CaO and high FeO/MnO of the bulk rocks, all argue for pyroxenite as the source lithology for the Wangjiadashan basalts. Both the trace element compositions of bulk rocks and the Pb isotope compositions of melt inclusions point to the recycled oceanic gabbro-associated EMI component present in the source. Low and high Pb isotope groups of the melt inclusions indicate the heterogeneous nature of sources for the Wangjiadashan basalts, thus highlighting the highly mantle heterogeneity in Shandong Province.

    basalt; olivine; pyroxenite; melt inclusion; mantle heterogeneity; Shandong Province

    P581; P597

    A

    0379-1726(2016)01-0001-23

    2015-02-10;

    2015-04-03;

    2015-04-24

    國家自然科學(xué)基金(91214202)

    張銀慧(1985–), 女, 博士, 巖石地球化學(xué)專業(yè)。E-mail: jiangno.1@163.com

    REN Zhong-yuan, E-mail: zyren@gig.ac.cn; Tel: +86-20-85292969

    猜你喜歡
    橄欖石輝石源區(qū)
    冬小麥蒸散源區(qū)代表性分析
    新疆和巴西薔薇輝石玉的對比研究
    綠色之星橄欖石
    化石(2021年1期)2021-03-16 01:20:50
    粉色薔薇輝石的寶石學(xué)及礦物學(xué)特征
    中國寶玉石(2020年6期)2021-01-15 10:25:14
    不同溫度、壓強(qiáng)、氧逸度條件下斜方輝石含水性的實(shí)驗(yàn)研究
    橄欖石項(xiàng)鏈
    渭河源區(qū)徑流量變化特征及趨勢分析
    天然橄欖石單晶的壓縮性*
    綠色小精靈橄欖石
    西部資源(2014年1期)2014-04-29 00:44:03
    江西廣豐早白堊世堿性橄欖玄武巖中單斜輝石礦物化學(xué)及其地質(zhì)意義
    国产乱来视频区| 欧美激情久久久久久爽电影| 人妻系列 视频| 有码 亚洲区| 国产一区二区在线观看日韩| 成人亚洲欧美一区二区av| 午夜免费观看性视频| 国产伦精品一区二区三区四那| 97超碰精品成人国产| 国产精品久久久久久精品古装| 下体分泌物呈黄色| 亚洲av日韩在线播放| 三级国产精品欧美在线观看| 日韩成人伦理影院| 啦啦啦啦在线视频资源| 国产一区有黄有色的免费视频| 久久精品熟女亚洲av麻豆精品| 日韩国内少妇激情av| 熟女电影av网| 最新中文字幕久久久久| 噜噜噜噜噜久久久久久91| 亚洲欧美日韩卡通动漫| 精品一区二区免费观看| 偷拍熟女少妇极品色| 丰满少妇做爰视频| 又爽又黄无遮挡网站| 国产欧美日韩一区二区三区在线 | 91久久精品国产一区二区成人| 成人无遮挡网站| 国产黄色免费在线视频| 国产男女超爽视频在线观看| 国产 一区 欧美 日韩| 成年免费大片在线观看| 欧美成人a在线观看| 亚洲成人中文字幕在线播放| 国产永久视频网站| 中文在线观看免费www的网站| 天堂俺去俺来也www色官网| 男人和女人高潮做爰伦理| 久久精品国产亚洲网站| 国产久久久一区二区三区| 中文字幕久久专区| 欧美+日韩+精品| 波野结衣二区三区在线| 国产午夜精品久久久久久一区二区三区| 亚洲一区二区三区欧美精品 | 欧美成人一区二区免费高清观看| videossex国产| 一区二区三区精品91| 午夜福利在线在线| 伦理电影大哥的女人| 国产色婷婷99| 国产色婷婷99| 国产乱人偷精品视频| 亚洲av免费高清在线观看| 日韩三级伦理在线观看| 日韩成人av中文字幕在线观看| 视频中文字幕在线观看| 一级黄片播放器| 欧美成人午夜免费资源| 久久久精品94久久精品| 国产亚洲av嫩草精品影院| 韩国高清视频一区二区三区| 直男gayav资源| 男女无遮挡免费网站观看| 国产老妇女一区| 国产亚洲5aaaaa淫片| 日韩制服骚丝袜av| 午夜免费鲁丝| 亚洲人与动物交配视频| 成人综合一区亚洲| 18+在线观看网站| 高清视频免费观看一区二区| 在线a可以看的网站| freevideosex欧美| av一本久久久久| 亚洲精品影视一区二区三区av| 卡戴珊不雅视频在线播放| 2021天堂中文幕一二区在线观| 国产有黄有色有爽视频| 久久久久久久久久久丰满| 中文乱码字字幕精品一区二区三区| 18禁在线无遮挡免费观看视频| 一级a做视频免费观看| 丰满乱子伦码专区| 日韩成人av中文字幕在线观看| 日本欧美国产在线视频| 国产v大片淫在线免费观看| 亚洲最大成人手机在线| 又粗又硬又长又爽又黄的视频| 日产精品乱码卡一卡2卡三| www.av在线官网国产| 91狼人影院| 五月开心婷婷网| 国产乱人视频| 一级片'在线观看视频| 一区二区三区四区激情视频| 精品人妻熟女av久视频| videos熟女内射| 欧美潮喷喷水| 欧美xxⅹ黑人| 男女国产视频网站| 国产v大片淫在线免费观看| 看免费成人av毛片| 亚洲自拍偷在线| 中文精品一卡2卡3卡4更新| 91久久精品国产一区二区成人| 久久精品夜色国产| 97在线人人人人妻| www.av在线官网国产| 亚洲精品日韩在线中文字幕| 欧美日韩在线观看h| 七月丁香在线播放| 免费观看无遮挡的男女| 精品人妻熟女av久视频| 国产色婷婷99| h日本视频在线播放| 亚洲av中文av极速乱| 中国三级夫妇交换| 欧美一区二区亚洲| 最近的中文字幕免费完整| 久久久久久久久久久丰满| 尾随美女入室| 成年女人看的毛片在线观看| 少妇人妻久久综合中文| 欧美激情在线99| 午夜爱爱视频在线播放| 日本一本二区三区精品| 日韩欧美一区视频在线观看 | freevideosex欧美| 精品人妻视频免费看| 国国产精品蜜臀av免费| 国语对白做爰xxxⅹ性视频网站| 香蕉精品网在线| 80岁老熟妇乱子伦牲交| 91久久精品电影网| 日韩三级伦理在线观看| 精品一区二区三区视频在线| 午夜免费观看性视频| 身体一侧抽搐| 国产精品人妻久久久影院| 六月丁香七月| 亚洲精品影视一区二区三区av| 18禁在线播放成人免费| 汤姆久久久久久久影院中文字幕| 在线看a的网站| 久久久久国产网址| 亚洲国产精品999| 欧美+日韩+精品| 午夜福利在线观看免费完整高清在| 久久久午夜欧美精品| 在线观看三级黄色| 亚洲综合精品二区| 亚洲婷婷狠狠爱综合网| 国产高清有码在线观看视频| 成人特级av手机在线观看| 亚洲av日韩在线播放| 午夜免费男女啪啪视频观看| 看黄色毛片网站| 国产免费福利视频在线观看| 久久久久久伊人网av| 69av精品久久久久久| 国产亚洲一区二区精品| 亚洲不卡免费看| 人人妻人人看人人澡| 欧美精品一区二区大全| 国产亚洲91精品色在线| 亚洲国产最新在线播放| 乱码一卡2卡4卡精品| 波野结衣二区三区在线| 一区二区av电影网| 国产精品一区二区在线观看99| 日韩欧美精品v在线| 爱豆传媒免费全集在线观看| videossex国产| 久久精品国产鲁丝片午夜精品| 国产探花极品一区二区| 国产精品国产三级专区第一集| 国产亚洲91精品色在线| 韩国av在线不卡| 成人一区二区视频在线观看| 久久久亚洲精品成人影院| eeuss影院久久| 久久综合国产亚洲精品| 女的被弄到高潮叫床怎么办| 免费大片黄手机在线观看| 国产精品成人在线| 亚洲av电影在线观看一区二区三区 | a级毛片免费高清观看在线播放| 高清欧美精品videossex| 精品99又大又爽又粗少妇毛片| 天天一区二区日本电影三级| 亚洲精品成人久久久久久| 欧美性感艳星| 久久久久性生活片| 高清视频免费观看一区二区| 少妇被粗大猛烈的视频| 亚洲高清免费不卡视频| 涩涩av久久男人的天堂| 欧美精品人与动牲交sv欧美| 夫妻午夜视频| 午夜福利视频精品| 丝瓜视频免费看黄片| 人人妻人人澡人人爽人人夜夜| 精品久久久久久久久亚洲| 久久精品久久精品一区二区三区| 亚洲不卡免费看| 国产精品一区二区三区四区免费观看| 亚洲美女视频黄频| 一二三四中文在线观看免费高清| 亚洲欧美日韩卡通动漫| av在线播放精品| 国产亚洲最大av| 国产精品久久久久久久电影| 久久精品熟女亚洲av麻豆精品| 欧美 日韩 精品 国产| 在线亚洲精品国产二区图片欧美 | 久久亚洲国产成人精品v| 日日啪夜夜撸| 久久精品综合一区二区三区| 亚洲国产精品成人久久小说| 成人美女网站在线观看视频| 日韩一区二区视频免费看| 欧美一级a爱片免费观看看| 免费不卡的大黄色大毛片视频在线观看| 日韩不卡一区二区三区视频在线| 日日撸夜夜添| 舔av片在线| 免费观看性生交大片5| 卡戴珊不雅视频在线播放| 18禁动态无遮挡网站| 国产精品一区二区三区四区免费观看| 亚洲一区二区三区欧美精品 | 国产精品99久久99久久久不卡 | 女人被狂操c到高潮| 看黄色毛片网站| 在线观看av片永久免费下载| 国产国拍精品亚洲av在线观看| 一区二区三区乱码不卡18| av黄色大香蕉| 九九在线视频观看精品| 少妇的逼好多水| 亚洲欧美中文字幕日韩二区| av在线亚洲专区| 亚洲欧美精品专区久久| 视频中文字幕在线观看| 国产精品伦人一区二区| 看十八女毛片水多多多| 草草在线视频免费看| 久久韩国三级中文字幕| 亚洲国产精品国产精品| 精品熟女少妇av免费看| 日本猛色少妇xxxxx猛交久久| 久久久久精品性色| 久久久欧美国产精品| 国产精品成人在线| 亚洲aⅴ乱码一区二区在线播放| 精品99又大又爽又粗少妇毛片| 午夜视频国产福利| 亚洲欧美精品专区久久| 日韩免费高清中文字幕av| 国产免费又黄又爽又色| 99久久九九国产精品国产免费| 国产伦精品一区二区三区四那| 日韩视频在线欧美| 九九爱精品视频在线观看| 成年女人在线观看亚洲视频 | 免费观看av网站的网址| 在线看a的网站| 自拍偷自拍亚洲精品老妇| 精华霜和精华液先用哪个| 九九久久精品国产亚洲av麻豆| 自拍欧美九色日韩亚洲蝌蚪91 | 永久免费av网站大全| 99精国产麻豆久久婷婷| 国产亚洲5aaaaa淫片| 亚洲综合色惰| 国产v大片淫在线免费观看| 成人亚洲精品一区在线观看 | 亚洲av中文av极速乱| 2018国产大陆天天弄谢| 国产在线一区二区三区精| 欧美xxⅹ黑人| 91久久精品国产一区二区三区| 日日啪夜夜爽| 精品人妻偷拍中文字幕| 91精品伊人久久大香线蕉| 国产一区有黄有色的免费视频| 久久精品人妻少妇| 国精品久久久久久国模美| 丝瓜视频免费看黄片| 在线观看av片永久免费下载| 午夜福利视频精品| 可以在线观看毛片的网站| 免费高清在线观看视频在线观看| 免费观看性生交大片5| 成年av动漫网址| 麻豆成人av视频| 777米奇影视久久| 国产精品国产三级国产专区5o| 国产精品.久久久| 日日撸夜夜添| 99热这里只有精品一区| 亚洲欧美精品自产自拍| 观看美女的网站| 在线亚洲精品国产二区图片欧美 | 在线a可以看的网站| 日本爱情动作片www.在线观看| 亚洲av在线观看美女高潮| 边亲边吃奶的免费视频| 久久韩国三级中文字幕| 精品一区二区三卡| 亚洲国产av新网站| 国产视频内射| 美女脱内裤让男人舔精品视频| 国国产精品蜜臀av免费| 亚洲天堂国产精品一区在线| 黄片无遮挡物在线观看| 在现免费观看毛片| 欧美高清成人免费视频www| 黄片wwwwww| 久久国内精品自在自线图片| 亚洲欧美日韩东京热| 爱豆传媒免费全集在线观看| 少妇熟女欧美另类| 色5月婷婷丁香| 亚洲国产欧美人成| 欧美 日韩 精品 国产| 午夜福利在线观看免费完整高清在| 视频中文字幕在线观看| 日韩电影二区| 中文欧美无线码| 97精品久久久久久久久久精品| 白带黄色成豆腐渣| 国产成人福利小说| 欧美日韩亚洲高清精品| 国产精品国产av在线观看| 插逼视频在线观看| 国产一区二区在线观看日韩| 最近中文字幕高清免费大全6| 99久国产av精品国产电影| 亚洲精品456在线播放app| 中国美白少妇内射xxxbb| 国产 一区精品| 午夜福利视频1000在线观看| 80岁老熟妇乱子伦牲交| 午夜福利在线在线| 王馨瑶露胸无遮挡在线观看| 80岁老熟妇乱子伦牲交| eeuss影院久久| 亚洲国产精品成人久久小说| 欧美+日韩+精品| 亚洲自拍偷在线| 日日摸夜夜添夜夜爱| av专区在线播放| 免费播放大片免费观看视频在线观看| 在线免费观看不下载黄p国产| av天堂中文字幕网| 观看美女的网站| 老师上课跳d突然被开到最大视频| 亚洲欧美精品专区久久| 热99国产精品久久久久久7| 免费看光身美女| 亚洲国产欧美在线一区| 久久久久久久大尺度免费视频| 国产免费一级a男人的天堂| 久久久久久久精品精品| 18+在线观看网站| 18禁裸乳无遮挡动漫免费视频 | 高清在线视频一区二区三区| 国产真实伦视频高清在线观看| 男插女下体视频免费在线播放| 免费看不卡的av| 亚洲av中文av极速乱| 午夜老司机福利剧场| 欧美高清成人免费视频www| 少妇人妻久久综合中文| 欧美3d第一页| av网站免费在线观看视频| 亚洲精品久久午夜乱码| 久久这里有精品视频免费| 国产精品av视频在线免费观看| 欧美成人一区二区免费高清观看| 国产精品99久久久久久久久| 国产亚洲91精品色在线| 亚洲美女搞黄在线观看| 2021天堂中文幕一二区在线观| 一本久久精品| 特级一级黄色大片| 色哟哟·www| 亚洲国产精品成人久久小说| 欧美区成人在线视频| 欧美另类一区| 涩涩av久久男人的天堂| 久久精品国产自在天天线| 麻豆成人午夜福利视频| 色网站视频免费| 午夜免费鲁丝| 成人一区二区视频在线观看| 日韩不卡一区二区三区视频在线| 97精品久久久久久久久久精品| 一个人看的www免费观看视频| 中国国产av一级| 国产一区有黄有色的免费视频| 黄片wwwwww| 亚洲欧美中文字幕日韩二区| 99精国产麻豆久久婷婷| 亚洲aⅴ乱码一区二区在线播放| 国产又色又爽无遮挡免| 亚洲av.av天堂| 99热这里只有是精品在线观看| 国产高清不卡午夜福利| 久久99蜜桃精品久久| 听说在线观看完整版免费高清| 在线免费观看不下载黄p国产| 国产色婷婷99| 国产熟女欧美一区二区| 又大又黄又爽视频免费| 久久综合国产亚洲精品| 春色校园在线视频观看| 精品一区二区三区视频在线| 国产精品一区二区三区四区免费观看| 一区二区av电影网| 亚洲精品中文字幕在线视频 | 国产白丝娇喘喷水9色精品| 少妇人妻 视频| 国产免费视频播放在线视频| 中文精品一卡2卡3卡4更新| 国产在视频线精品| 亚洲不卡免费看| 亚洲第一区二区三区不卡| 久久久久久久精品精品| 午夜免费男女啪啪视频观看| 免费不卡的大黄色大毛片视频在线观看| 99热全是精品| 插逼视频在线观看| 免费观看性生交大片5| 蜜桃久久精品国产亚洲av| 婷婷色麻豆天堂久久| 丰满少妇做爰视频| 亚洲国产日韩一区二区| 国产精品秋霞免费鲁丝片| 99热这里只有是精品在线观看| 麻豆久久精品国产亚洲av| 免费观看性生交大片5| 国产淫语在线视频| 免费播放大片免费观看视频在线观看| 久久99热6这里只有精品| 国产高潮美女av| 精品久久久久久久末码| 成年版毛片免费区| 欧美一区二区亚洲| 国产精品国产av在线观看| 大码成人一级视频| 午夜精品国产一区二区电影 | 国产精品人妻久久久久久| 亚州av有码| 在线免费观看不下载黄p国产| 国产午夜福利久久久久久| 国产精品人妻久久久影院| 别揉我奶头 嗯啊视频| 亚洲熟女精品中文字幕| 男人爽女人下面视频在线观看| 天天一区二区日本电影三级| 狂野欧美激情性xxxx在线观看| 夜夜爽夜夜爽视频| 中国三级夫妇交换| 99热全是精品| 青春草国产在线视频| 成人欧美大片| 少妇丰满av| 干丝袜人妻中文字幕| av在线观看视频网站免费| 欧美3d第一页| 九九久久精品国产亚洲av麻豆| av在线app专区| 精品国产三级普通话版| videos熟女内射| 国产成人免费观看mmmm| 久久久a久久爽久久v久久| 又粗又硬又长又爽又黄的视频| 成年女人在线观看亚洲视频 | 日韩三级伦理在线观看| 一级毛片aaaaaa免费看小| 色哟哟·www| 午夜视频国产福利| 久久久久久久精品精品| 久久女婷五月综合色啪小说 | 久久鲁丝午夜福利片| 久久久久久久精品精品| 在线a可以看的网站| 精品久久久久久久末码| 人妻夜夜爽99麻豆av| 国产精品一二三区在线看| 国产精品福利在线免费观看| 特大巨黑吊av在线直播| 国产成人免费观看mmmm| 人妻 亚洲 视频| 色视频在线一区二区三区| 亚洲精华国产精华液的使用体验| 亚洲色图综合在线观看| 欧美日韩综合久久久久久| 国产精品嫩草影院av在线观看| 亚洲av欧美aⅴ国产| 啦啦啦啦在线视频资源| 99久久中文字幕三级久久日本| 中文字幕久久专区| 99九九线精品视频在线观看视频| 国产毛片a区久久久久| 最后的刺客免费高清国语| 日本黄大片高清| 久久精品国产亚洲网站| 少妇被粗大猛烈的视频| 涩涩av久久男人的天堂| 夫妻午夜视频| 看免费成人av毛片| 禁无遮挡网站| 久久久午夜欧美精品| 国产91av在线免费观看| 亚洲自偷自拍三级| 久久亚洲国产成人精品v| 色播亚洲综合网| 夜夜看夜夜爽夜夜摸| 国产黄色视频一区二区在线观看| 日韩伦理黄色片| 国产一区有黄有色的免费视频| 男女边吃奶边做爰视频| 日韩国内少妇激情av| 一区二区三区免费毛片| 22中文网久久字幕| 欧美 日韩 精品 国产| 七月丁香在线播放| 欧美日本视频| 国产成人a∨麻豆精品| 国产精品一区二区性色av| 国产乱来视频区| 精品国产三级普通话版| 精品久久久久久久人妻蜜臀av| 一区二区三区四区激情视频| 美女脱内裤让男人舔精品视频| 最近2019中文字幕mv第一页| 国产一区二区在线观看日韩| 亚洲在线观看片| 听说在线观看完整版免费高清| 白带黄色成豆腐渣| 黄色视频在线播放观看不卡| 少妇被粗大猛烈的视频| 男女无遮挡免费网站观看| 午夜亚洲福利在线播放| 日本-黄色视频高清免费观看| 久久99热6这里只有精品| 少妇人妻 视频| 欧美潮喷喷水| 久久韩国三级中文字幕| 交换朋友夫妻互换小说| 女的被弄到高潮叫床怎么办| 国产av不卡久久| 久久精品国产亚洲av涩爱| 高清欧美精品videossex| 亚洲av国产av综合av卡| 久久精品国产a三级三级三级| 男女边吃奶边做爰视频| 久久女婷五月综合色啪小说 | 韩国高清视频一区二区三区| 身体一侧抽搐| 免费观看性生交大片5| 国产成人精品一,二区| 国产精品精品国产色婷婷| 十八禁网站网址无遮挡 | 少妇人妻久久综合中文| 在线观看人妻少妇| 国产一区亚洲一区在线观看| 国产精品人妻久久久久久| 六月丁香七月| 2022亚洲国产成人精品| 国产日韩欧美亚洲二区| 舔av片在线| 精品国产露脸久久av麻豆| 99热这里只有是精品在线观看| 国产精品秋霞免费鲁丝片| 赤兔流量卡办理| 18禁在线无遮挡免费观看视频| 日产精品乱码卡一卡2卡三| 久久精品国产鲁丝片午夜精品| 成人特级av手机在线观看| 久久久久国产网址| 国产美女午夜福利| 久久久成人免费电影| 中国美白少妇内射xxxbb| 一区二区三区免费毛片| 国产黄色免费在线视频| 久久精品国产鲁丝片午夜精品| 欧美性感艳星| 国产成人精品久久久久久| 成人二区视频| 精品久久国产蜜桃| 成年版毛片免费区| 黄色视频在线播放观看不卡| 又大又黄又爽视频免费| 午夜视频国产福利| 一区二区三区乱码不卡18| 汤姆久久久久久久影院中文字幕| 亚洲av在线观看美女高潮| 超碰av人人做人人爽久久| 精品熟女少妇av免费看| 欧美3d第一页| 国产成人精品婷婷| 黄色日韩在线| 精品一区二区三区视频在线| 又大又黄又爽视频免费| 80岁老熟妇乱子伦牲交| 亚洲精品影视一区二区三区av| 国产精品国产三级专区第一集| 97在线人人人人妻| 在线观看国产h片|