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

    華北地幔捕虜體Fe-Mg-Ca同位素研究進(jìn)展*

    2022-12-29 11:14:52趙新苗賈喜良易春霞張兆峰朱祥坤
    巖石學(xué)報(bào) 2022年12期
    關(guān)鍵詞:橄欖巖巖石圈輝石

    趙新苗 賈喜良, 2 易春霞, 2 張兆峰 朱祥坤

    克拉通是地球上古老的大陸,它最基本的特征是穩(wěn)定,而克拉通能夠長(zhǎng)期穩(wěn)定存在是因?yàn)樗哂泄爬喜⑶揖藓竦膸r石圈地幔。我國(guó)華北是典型的克拉通,而越來(lái)越多的研究表明,華北克拉通巖石圈地幔發(fā)生了破壞,失去了穩(wěn)定性,其主要原因是它的巖石圈地幔的組成和性質(zhì)發(fā)生了根本性的轉(zhuǎn)變,是自顯生宙以來(lái)全球古老克拉通破壞最激烈的地區(qū),因而是研究大陸形成與演化的最佳場(chǎng)所(Menziesetal., 1993; Griffinetal., 1998; Fanetal., 2000; Carlsonetal., 2005)。研究巖石圈地幔組成和演化的理想樣品是地幔橄欖巖捕虜體和捕虜晶,這是因?yàn)樗鼈兪莵?lái)源于巖石圈地幔最直接的樣品,能夠反映巖石圈地幔組成和演化的直接信息。學(xué)者們通過(guò)對(duì)華北克拉通地幔橄欖巖捕虜體和捕虜晶的研究發(fā)現(xiàn),華北巖石圈地幔組成和性質(zhì)自顯生宙以來(lái)發(fā)生了根本性改變,在古生代的時(shí)候它還具有古老而且非常厚的克拉通型巖石圈地幔,然而到了新生代,它變成了比較薄而且年輕的大洋型巖石圈地幔(高山等,2009;徐義剛等,2009;張宏福,2009;鄭建平,2009;鄭永飛和吳福元,2009;朱日祥等,2011, 2012, 2020及其參考文獻(xiàn))。而大量研究表明,橄欖巖-熔體反應(yīng)是華北巖石圈地幔性質(zhì)和組成發(fā)生轉(zhuǎn)變的重要途徑,為認(rèn)清克拉通破壞本質(zhì)奠定了基礎(chǔ)(張宏福,2009;Tangetal., 2013, 2021及其參考文獻(xiàn))。但是如何認(rèn)識(shí)華北巖石圈地幔組成和屬性的轉(zhuǎn)變過(guò)程、尤其是地幔改造過(guò)程中熔體的性質(zhì)和來(lái)源,是鑒別不同克拉通破壞機(jī)制的關(guān)鍵。由于華北幔源巖石的研究已開(kāi)展了近四十年,已有的研究表明僅依靠傳統(tǒng)的研究手段未能很好地解決上述問(wèn)題,亟需創(chuàng)新思路和通過(guò)新途徑對(duì)華北克拉通破壞關(guān)鍵地區(qū)和時(shí)期的地幔捕虜體做進(jìn)一步的研究。

    非傳統(tǒng)穩(wěn)定同位素(如Fe、Mg和Ca等)地球化學(xué)是當(dāng)前同位素地球化學(xué)研究領(lǐng)域最具活力的方向,為研究大陸巖石圈演化和地幔物質(zhì)再循環(huán)等地幔地球化學(xué)帶來(lái)了全新的視野和手段(張宏福等,2007;周新華等,2013;朱祥坤等,2013;趙新苗,2019;張宏福和陳立輝,2021)。Fe、Mg和Ca是大陸巖石圈地幔的主要組成元素,在與熔體/流體有關(guān)的許多地質(zhì)過(guò)程中會(huì)發(fā)生活化遷移和再分配,從而可能引起同位素組成的變化,而這一變化直接反映了巖石圈地幔組成的變化,因此對(duì)熔體的性質(zhì)和來(lái)源判斷更加直接。本文對(duì)目前已報(bào)道的華北克拉通不同地區(qū)、不同類(lèi)型地幔捕虜體Fe、Mg和Ca同位素組成進(jìn)行了總結(jié),揭示華北克拉通大陸巖石圈地幔Fe、Mg和Ca同位素組成的空間差異性及其可能原因,呈現(xiàn)Fe、Mg和Ca同位素體系在示蹤大陸巖石圈地幔演化過(guò)程中的重要作用和發(fā)展前景。

    1 Fe、Mg和Ca同位素體系概述

    1.1 Fe同位素體系簡(jiǎn)介

    鐵的原子序數(shù)是26,為第4周期、第VIII族過(guò)渡金屬元素。Fe元素在在自然界中廣泛分布,它不僅是地球和類(lèi)地行星的主要元素之一,還是地球上豐度最高的變價(jià)(Fe0、Fe2+、Fe3+)金屬元素,并且與生命活動(dòng)密切相關(guān)。因此,F(xiàn)e的地球化學(xué)循環(huán)與成礦作用、地幔深部過(guò)程、氧化還原條件過(guò)程以及環(huán)境與生命演化等重大科學(xué)問(wèn)題有密切關(guān)系。Fe在自然界中存在四個(gè)穩(wěn)定同位素:54Fe、56Fe、57Fe和58Fe,其相對(duì)豐度分別為5.85%、91.75%、2.12%、0.28%(Berglund and Wieser 2011)。近十幾年來(lái),隨著同位素分析技術(shù)的提高,尤其是多接收電感耦合等離子體質(zhì)譜儀(multi-collector inductively coupled plasma mass spectrometry)的開(kāi)發(fā)和應(yīng)用,F(xiàn)e同位素分析測(cè)試精度得到顯著提高,普遍可達(dá)0.05‰(2SD)或更高(Zhuetal., 2000; Weyer and Schwieters 2003; Poitrasson and Freydier 2005; Dauphasetal., 2009; Craddock and Dauphas, 2011; Heetal., 2015)。在此基礎(chǔ)上,近年來(lái)學(xué)者們開(kāi)展了大量的Fe同位素研究工作,大體上查明了地球各個(gè)主要儲(chǔ)庫(kù)的Fe同位素組成特征,隨后大大開(kāi)拓了其在地球化學(xué)、天體化學(xué)、環(huán)境科學(xué)以及生命科學(xué)等方面的應(yīng)用(Dauphasetal., 2017; Johnsonetal., 2020及其參考文獻(xiàn))。Fe同位素組成的分析結(jié)果通常用相對(duì)于國(guó)際標(biāo)準(zhǔn)物質(zhì)IRMM-014的千分偏差δ值表示:

    δ56/54Fe(‰)=[(56Fe/54Fe)樣品/(56Fe/54Fe)標(biāo)準(zhǔn)-1]×103

    δ57/54Fe(‰)=[(57Fe/54Fe)樣品/(57Fe/54Fe)標(biāo)準(zhǔn)-1]×103

    自然界中的鐵同位素符合質(zhì)量分餾:δ56/54Fe≈0.678×δ57/54Fe(Zhuetal., 2000)。

    1.2 Mg同位素體系簡(jiǎn)介

    鎂的原子系數(shù)是12,為第3周期、第IIA族堿土金屬元素。Mg是地球中豐度僅次于氧和硅的主要造巖元素,主要富集于地幔中,參與殼幔幾乎所有的地球化學(xué)過(guò)程,同時(shí)廣泛分布于水圈和生物圈中。Mg在自然界中有3個(gè)穩(wěn)定同位素:24Mg、25Mg和26Mg,相對(duì)豐度分別為78.99%、10.00%和11.01%(Berglund and Wieser, 2011)。24Mg和26Mg之間較大的相對(duì)質(zhì)量差(~8%)使得Mg同位素在自然界中產(chǎn)生顯著的同位素分餾(Teng, 2017及其參考文獻(xiàn))。Mg同位素組成的分析結(jié)果通常用相對(duì)于國(guó)際標(biāo)準(zhǔn)物質(zhì)DSM3(Galyetal., 2003)的千分偏差δ值表示:

    δ26/24Mg(‰)=[(26Mg/24Mg)樣品/(26Mg/24Mg)標(biāo)準(zhǔn)-1]×103

    δ25/24Mg(‰)=[(25Mg/24Mg)樣品/(25Mg/24Mg)標(biāo)準(zhǔn)-1]×103

    1.3 Ca同位素體系簡(jiǎn)介

    鈣的原子系數(shù)是20,為第四周期、第IIA族堿土金屬元素。Ca廣泛分布在地球的各個(gè)圈層中,它不僅是主要的造巖元素之一,也是生物必須的元素之一。自然界Ca有40Ca、42Ca、43Ca、44Ca、46Ca和48Ca六種同位素,其相對(duì)豐度為96.941%、0.647%、0.135%、2.086%、0.004%、0.187%(Coplen, 2011)。Ca具有最大的相對(duì)質(zhì)量差(48Ca和40Ca相差20%),僅次于氫和氦元素,這一大的相對(duì)質(zhì)量差使得很多地質(zhì)過(guò)程中Ca同位素能夠發(fā)生顯著的質(zhì)量分餾。在自然界中已發(fā)現(xiàn)δ44/40Ca有~6‰的變化,這樣的特性決定了Ca同位素體系有助于了解地球和各圈層的演化與相互作用,也有助于對(duì)太陽(yáng)系和行星演化的認(rèn)識(shí),同時(shí)也有可能對(duì)示蹤深部碳循環(huán)、重建古氣候變化、反演古環(huán)境及其限定礦床成因等重大理論和應(yīng)用問(wèn)題提供新的認(rèn)識(shí)(DePaolo, 2004)。Ca同位素組成的分析結(jié)果通常用相對(duì)于國(guó)際標(biāo)準(zhǔn)物質(zhì)NIST SRM915a的千分偏差δ值表示:δ44/40Ca(‰)=[(44Ca/40Ca)樣品/(44Ca/40Ca)標(biāo)準(zhǔn)-1] ×103。

    2 華北地幔捕擄體Fe-Mg-Ca同位素組成

    2.1 華北地幔捕虜體的Fe同位素組成特征

    鐵不僅是巖石圈地幔直接樣品——橄欖巖捕虜體的主要組成元素,還是地幔最主要的變價(jià)(Fe2+和Fe3+)元素,在與熔體/流體有關(guān)的許多地質(zhì)過(guò)程中會(huì)發(fā)生活化遷移和再分配。因此,對(duì)地幔捕虜Fe同位素組成的研究可以為深入認(rèn)識(shí)大陸巖石圈地幔的演化過(guò)程,尤其是地幔改造過(guò)程中熔體的性質(zhì)和來(lái)源帶來(lái)新的切入點(diǎn)。已有的研究結(jié)果表明,地幔橄欖巖和輝石巖捕虜體及其中單礦物的Fe同位素組成具有非常大的變化范圍(圖1、圖2;Johnsonetal., 2020及其參考文獻(xiàn)),尤其是我國(guó)華北地區(qū)地幔捕虜體的Fe同位素組成有~2‰的變化(δ57/54Fe=-1.04‰~1.00‰;圖1、圖2;Zhaoetal., 2010, 2012, 2015, 2017a, b)。中奧陶世金伯利巖的噴發(fā)是華北克拉通化后首次出現(xiàn)的巖漿活動(dòng)。因此,這些金伯利巖及其攜帶的橄欖巖捕虜體可以很好的限定克拉通破壞前的巖石圈狀態(tài)。筆者對(duì)華北遼寧復(fù)縣和山東蒙陰地區(qū)石榴子石二輝橄欖巖和尖晶石二輝橄欖巖開(kāi)展了Fe同位素研究,結(jié)果發(fā)現(xiàn)這些地幔橄欖巖Fe同位素組成變化范圍較大(δ57/54Fe=-0.09‰~0.55‰,平均值為0.17±0.12‰,2SE;n=12;趙新苗等,未發(fā)表數(shù)據(jù)),并且其Fe同位素組成隨其燒失量、CaO含量及其Mg#[=100×Mg2+/(Mg2++Fe2+)]的升高而降低,說(shuō)明蛇紋石化和碳酸鹽化過(guò)程會(huì)引起橄欖巖Fe同位素組成的改變,無(wú)法準(zhǔn)確反映原巖的Fe同位素組成信息。為此,筆者選擇了我國(guó)華南地區(qū)古生代大洪山金伯利巖中攜帶的新鮮的,沒(méi)有遭受蝕變的石榴子石二輝橄欖巖進(jìn)行了Fe同位素研究,結(jié)果表明大洪山地區(qū)的石榴子石二輝橄欖巖Fe同位素組成比較均一(δ57/54Fe=0.02‰~0.03‰;Zhaoetal., 2012),其平均值(δ57/54Fe=0.03±0.01‰,2SE;n=2)與球粒隕石的平均值接近(δ57/54Fe=-0.01±0.01‰,2SE;n=146;Johnsonetal., 2020),并且與Weyer and Ionov(2007)根據(jù)橄欖巖Mg#=89.5估算的飽滿的上地幔(Fertile upper mantle)的Fe同位素組成接近(δ57/54Fe=0.03±0.03‰)。中生代是華北克拉通巖石圈減薄的關(guān)鍵時(shí)期。所以,華北中生代玄武質(zhì)巖石攜帶的地幔橄欖巖捕虜體是探討華北克拉通型巖石圈之所以能夠被破壞和減薄的最直接的途徑。筆者對(duì)華北大西莊地區(qū)白堊紀(jì)玄武巖石攜帶的地幔橄欖巖及其單礦物進(jìn)行了詳細(xì)的Fe同位素研究,發(fā)現(xiàn)中生代玄武巖攜帶的地幔橄欖巖Fe同位素組成變化范圍比較小(δ57/54Fe=-0.36‰~0.12‰;趙新苗等,未發(fā)表數(shù)據(jù)),其平均值(δ57/54Fe=0.00±0.07‰,2SE;n=12)與球粒隕石的平均值(δ57/54Fe=-0.01±0.01‰,2SE;n=146;Johnsonetal., 2020)和飽滿的上地幔Fe同位素組成接近(δ57/54Fe=0.03±0.03‰,2SD; Weyer and Ionov, 2007)。而華北克拉通典型地區(qū)如北緣(漢諾壩、三義堂)、內(nèi)部(陽(yáng)原、鶴壁)和東部郯廬斷裂帶內(nèi)部(昌樂(lè))的新生代玄武巖攜帶的地幔捕虜體及其單礦物具有很大的Fe同位素變化范圍(δ57/54Fe=-1.04‰~ 1.00‰;圖1、圖2;Zhaoetal., 2010, 2012, 2015, 2017a, b),不同類(lèi)型地幔捕虜體及其單礦物的Fe同位素組成不同(圖1、圖2;Zhaoetal., 2010, 2012, 2015, 2017a, b)。例如華北北緣漢諾壩、三義堂、陽(yáng)原、南北重力梯度帶鶴壁地區(qū)以及郯廬斷裂帶昌樂(lè)地區(qū)的鎂橄欖巖(Mg#≥90)的Fe同位素組成變化范圍(δ57/54Fe=-0.28‰~0.19‰;平均值δ57/54Fe=-0.01± 0.05‰,2SE;n=28)與華北中生代玄武質(zhì)巖石攜帶的地幔橄欖巖捕虜體(δ57/54Fe=-0.36‰~0.12‰;趙新苗等,未發(fā)表數(shù)據(jù))、華南地區(qū)鎂橄欖巖(δ57/54Fe=-0.29‰~0.28‰;平均值δ57/54Fe=0.02±0.04‰,2SE;n=30;Huangetal., 2011; Zhaoetal., 2021)以及世界其它地區(qū)的鎂橄欖巖的Fe同位素組成接近(δ57/54Fe=-0.22‰~0.22‰;平均值δ57/54Fe=0.01±0.02‰,2SE;n=79;圖1)。同時(shí)這類(lèi)橄欖巖的不同共生礦物之間存在Fe同位素分餾,并且大多數(shù)高溫礦物之間具有δ57/54Fe橄欖石<δ57/54Fe輝石<δ57/54Fe尖晶石的規(guī)律(圖2;Zhaoetal., 2010, 2012, 2015, 2017a)。與華北鎂橄欖巖不同,華北陽(yáng)原和郯廬斷裂帶昌樂(lè)地區(qū)的富Fe橄欖巖(Mg#<89)的Fe同位素組成(δ57/54Fe=-1.04‰~-0.33‰;圖1)顯著低于華北地區(qū)鎂橄欖巖的Fe同位素組成,尤其是在華北巖石圈減薄強(qiáng)烈的地區(qū)(如郯廬斷裂帶內(nèi)昌樂(lè))地幔橄欖巖捕虜體的Fe同位素組成存在~1.3‰的巨大變化,并具有明顯低于全球其他地區(qū)輕的Fe同位素組成(δ57/54Fe低至-1.04‰;圖1)。這些富Fe橄欖巖組成礦物間Fe同位素分餾較大(圖2),礦物間Fe同位素組成沒(méi)有完全分布在平衡分餾線上(Zhaoetal., 2010, 2012, 2015, 2017a),說(shuō)明這些單礦物間的Fe同位素組成尚沒(méi)有達(dá)到平衡。與鎂橄欖巖和富Fe橄欖巖不同,華北昌樂(lè)地區(qū)的異剝橄欖巖(Mg#<88)明顯富集重Fe同位素(δ57/54Fe=0.17‰~0.26‰,平均值為0.21±0.09‰,2SE;n=2;圖2),其平均值比世界其它地區(qū)的異剝橄欖巖Fe同位素平均值偏重(δ57/54Fe=0.05‰~0.17‰,平均值為0.09±0.04‰,2SE;n=5;圖2),而與華北新生代玄武巖的Fe同位素平均值接近(δ57/54Fe=0.20±0.02‰,2SE;n=9;Zhaoetal., 2012)。同時(shí)異剝橄欖巖的共生礦物之間δ57/54Fe橄欖石<δ57/54Fe單斜輝石,基本處于理論計(jì)算的Fe同位素平衡分餾范圍內(nèi)(Zhaoetal., 2012)。

    圖1 華北地幔橄欖巖和輝石巖捕虜體的Fe同位素組成地幔捕虜體數(shù)據(jù)來(lái)源:華北地幔橄欖巖、異剝橄欖巖及輝石巖Fe同位素?cái)?shù)據(jù)來(lái)自Zhao et al.(2010, 2012, 2015, 2017a, b);華南地幔橄欖巖Fe同位素?cái)?shù)據(jù)來(lái)自Huang et al.(2011)和Zhao et al.(2021);其他地區(qū)地幔橄欖巖、異剝橄欖巖和輝石巖數(shù)據(jù)來(lái)自Beard and Johnson(2004),Williams et al.(2005),Schoenberg and von Blanckenburg(2006),Weyer and Ionov(2007),Poitrasson et al.(2013),Williams and Bizimis(2014),Macris et al.(2015)和An et al.(2017). 黑色虛線及其棕色條帶代表球粒隕石平均值(δ57/54Fe=-0.01±0.01‰,2SE;n=146;Johnson et al., 2020). 圖2數(shù)據(jù)來(lái)源同此圖Fig.1 Fe isotopic compositions of mantle peridotite and pyroxene xenoliths from the North China Craton

    圖2 華北地幔橄欖巖和輝石巖中共存單礦物的Fe同位素組成變化Fig.2 Fe isotopic compositions of coexisting minerals in mantle peridotites and pyroxene xenoliths from the North China Craton

    輝石巖是巖石圈地幔僅次于地幔橄欖巖的重要組成部分,它包含著豐富的巖石圈地幔演化信息,是認(rèn)識(shí)地幔不均一性和玄武巖源區(qū)組成的關(guān)鍵對(duì)象。筆者對(duì)華北漢諾壩地區(qū)不同類(lèi)型輝石巖捕虜體進(jìn)行了詳細(xì)的Fe同位素地球化學(xué)研究,發(fā)現(xiàn)綠色型Cr-透輝石輝石巖和黑色型Al-普通輝石二輝巖的Fe同位素組成變化范圍較小(δ57/54Fe=-0.01‰~0.09‰),其平均值(δ57/54Fe=0.03±0.04‰,2SE;n=5;圖1;Zhaoetal., 2017b)與估算的飽滿上地幔的Fe同位素組成在誤差范圍內(nèi)一致(δ57/54Fe=0.03±0.03‰,2SD; Weyer and Ionov, 2007)。與全巖相似,Cr-透輝石輝石巖和黑色型Al-普通輝石二輝巖中單礦物的Fe同位素組成變化范圍也比較小,礦物間的δ57/54Fe橄欖石<δ57/54Fe輝石<δ57/54Fe尖晶石(圖2;Zhaoetal., 2017b),不同礦物之間Fe同位素分餾均達(dá)到平衡。而漢諾壩石榴子石輝石巖和金云母單斜輝石巖具有比綠色型Cr-透輝石輝石巖和黑色型Al-普通輝石二輝巖顯著偏重的Fe同位素組成(δ57/54Fe=0.12‰~0.30‰;圖1;Zhaoetal., 2017b),其平均值(δ57/54Fe=0.20±0.04‰, 2SE;n=7 ;圖1;Zhaoetal., 2017b),與華北新生代玄武巖的Fe同位素平均值接近(δ57/54Fe=0.20±0.02‰,2SE;n=9;Zhaoetal., 2012)。這些石榴子石輝石巖中單礦物的Fe同位素組成具有較大的變化范圍,礦物間的δ57/54Fe橄欖石<δ57/54Fe輝石<δ57/54Fe尖晶石<δ57/54Fe石榴子石(圖2;Zhaoetal., 2017b)。值得注意的是,漢諾壩地區(qū)的金云母單斜輝石巖具有非常重的Fe同位素組成(δ57/54Fe=0.99‰),是目前發(fā)表的幔源巖石中具有最重Fe同位素組成的樣品(圖1;Zhaoetal., 2017b)。與全巖相似,該樣品的組成單礦物單斜輝石和金云母也具有非常重的Fe同位素組成(單斜輝石的δ57/54Fe=1.00‰,金云母的δ57/54Fe=0.95‰;圖2;Zhaoetal., 2017b)。以上研究表明,華北地區(qū)不同類(lèi)型輝石巖捕虜體也具有明顯不均一的Fe同位素組成,礦物間存在明顯的Fe同位素分餾。

    圖3 華北地幔橄欖巖和輝石巖捕虜體Mg同位素組成地幔捕虜體數(shù)據(jù)來(lái)源:華北地幔橄欖巖、異剝橄欖巖及輝石巖Mg同位素?cái)?shù)據(jù)來(lái)自Yang et al.(2009),Teng et al.(2010),Liu et al.(2011),Xiao et al.(2013),Hu et al.(2016),Zhao et al.(2017b)和Hu et al.(2019);華南地幔橄欖巖Mg位素?cái)?shù)據(jù)來(lái)自Huang et al.(2011);其他地區(qū)地幔橄欖巖和輝石巖數(shù)據(jù)來(lái)自Handler et al.(2009),Bourdon et al.(2010),Teng et al.(2010),Pogge von Strandmann et al.(2011),Lai et al.(2015),An et al.(2017)和Stracke et al.(2018). 黑色虛線及其棕色條帶代表硅酸鹽地球平均值(δ26/24Mg=-0.25±0.07‰,2SD;n=139;Teng et al., 2010). 圖4數(shù)據(jù)來(lái)源同此圖Fig.3 Mg isotopic compositions of mantle peridotite and pyroxene xenoliths from the North China Craton

    圖4 華北地幔橄欖巖和輝石巖中共存單礦物的Mg同位素組成變化Fig.4 Mg isotopic composition of coexisting minerals in mantle peridotites and pyroxene xenoliths from the North China Craton

    2.2 華北地幔捕虜體的Mg同位素特征

    與鐵元素相同,Mg也是地幔的主要組成元素,地球中絕大多數(shù)(>99%)Mg存在地幔中。因此,地幔的平均Mg同位素組成基本上代表了地球的Mg同位素組成。前人研究結(jié)果表明,來(lái)自世界不同地區(qū)的地幔橄欖巖捕虜體及其單礦物的Mg同位素組成變化范圍較小(δ26/24Mg=-0.44‰~0.00‰;圖3、圖4),其平均值(δ26/24Mg=-0.24±0.03‰,2SE;n=172)與球粒隕石平均值在誤差范圍內(nèi)相同(δ26/24Mg=-0.28±0.01‰,2SE;n=38;Tengetal., 2010)。然而值得注意的是我國(guó)華北克拉通新生代玄武巖攜帶的不同類(lèi)型地幔橄欖巖和輝石巖捕虜體及其單礦物的Mg同位素組成是明顯不均一的(δ26/24Mg=-1.53‰~-0.05‰;圖3、圖4;Yangetal., 2009; Tengetal., 2010; Liuetal., 2011; Xiaoetal., 2013; Huetal., 2016, 2019; Zhaoetal., 2017b; Huetal., 2019)。例如華北地區(qū)(漢諾壩、三義堂、寬甸以及昌樂(lè)地區(qū)的)鎂橄欖巖具有比較小的Mg同位素組成變化范圍(δ26/24Mg=-0.39‰~-0.19‰;圖3;Yangetal., 2009; Tengetal., 2010; Liuetal., 2011; Xiaoetal., 2013; Huetal., 2016),與華南地區(qū)二輝橄欖巖(δ26/24Mg=-0.38‰~-0.20‰;Huangetal., 2011)以及世界其他地區(qū)的二輝橄欖巖Mg 同位素組成在誤差范圍類(lèi)相似(δ26/24Mg=-0.34‰~-0.13‰;Bourdonetal., 2010; Tengetal., 2010; Pogge von Strandmannetal., 2011; Laietal., 2015; Anetal., 2017),并且這些華北鎂橄欖巖的Mg同位素組成平均值(δ26/24Mg=-0.26±0.01‰,2SE;n=51)與硅酸鹽地球Mg同位素組成相同(δ26/24Mg=-0.25±0.07‰,2SD;n=139;Tengetal., 2010)。同時(shí)這些華北鎂橄欖巖中共生的斜方輝石(δ26/24Mg=-0.28‰~-0.17‰)和橄欖石具有相似的Mg同位素組成(δ26/24Mg=-0.33‰~-0.18‰),而單斜輝石具有比共生的橄欖石系統(tǒng)更重的Mg同位素組成(δ26/24Mg=-0.29‰~-0.08‰),石榴子石具有比其它硅酸巖礦物更輕的Mg同位素組成(δ26/24Mg=-0.53‰~-0.52‰;圖4)。與共生的硅酸鹽礦物相比,尖晶石具有顯著更重的同位素組成,并且與溫度具有很好的相關(guān)性(δ26/24Mg=0.01‰~0.28‰;圖4;Yangetal., 2009; Tengetal., 2010; Liuetal., 2011; Xiaoetal., 2013; Huetal., 2016)。與華北鎂橄欖巖不同,華北昌樂(lè)地區(qū)的富Fe橄欖巖具有比鎂橄欖巖系統(tǒng)偏重的Mg同位素組成(δ26/24Mg=-0.33‰~-0.05‰,平均值δ26/24Mg=-0.18±0.05‰,2SE;n=10;圖3;Yangetal., 2009; Xiaoetal., 2013),這些富Fe橄欖巖礦物之間Mg同位素分餾較大(圖4),遠(yuǎn)遠(yuǎn)超出了這些礦物之間Mg同位素平衡分餾理論計(jì)算值和實(shí)際測(cè)量值,表明這些單礦物間的Mg同位素分餾沒(méi)有達(dá)到平衡(Yangetal., 2009; Xiaoetal., 2013)。Yangetal.(2009)和Xiaoetal.(2013)研究還發(fā)現(xiàn)華北昌樂(lè)地區(qū)的異剝橄欖巖僅有3個(gè)異剝橄欖巖的Mg同位素組成(δ26/24Mg=-0.19‰~0.09‰)與世界其他地區(qū)異剝橄欖巖的Mg同位素組成相似(δ26/24Mg=-0.28‰~0.00‰;圖3),其余12個(gè)異剝橄欖巖的Mg同位素組成(δ26/24Mg=-0.44‰~-0.26‰;圖3)整體比硅酸鹽地球Mg同位素組成偏輕(δ26/24Mg=-0.25±0.07‰,2SD;n=139;Tengetal., 2010),但是其共生礦物之間Mg同位素分餾落在了平衡Mg同位素理論計(jì)算值和實(shí)際測(cè)量值范圍內(nèi),表明這些單礦物間的Mg同位素分餾達(dá)到平衡(Yangetal., 2009; Xiaoetal., 2013)。

    Huetal.(2016)、Zhaoetal.(2017b)以及Huetal.(2019)研究發(fā)現(xiàn)華北漢諾壩地區(qū)的輝石巖及其單礦物的Mg同位素組成非常不均一(δ26/24Mg=-1.54‰~0.06‰;圖3、圖4),其中不含石榴子石的Cr-透輝石輝石巖和黑色型Al-普通輝石二輝巖的Mg同位素組成變化范圍比較小(δ26/24Mg=-0.37‰~-0.17‰;Liuetal., 2011; Huetal., 2016, 2019; Zhaoetal., 2017b; 圖3),其平均值(δ26/24Mg=-0.26‰± 0.03‰,2SE;n=11)與球粒隕石(δ26/24Mg=-0.28±0.01‰,2SE;n=38;Tengetal., 2010)以及大洋玄武巖具有相似的Mg同位素組成(δ26/24Mg=-0.25 ±0.01‰,2SE;n=110;Tengetal., 2010)。與全巖相似,這兩類(lèi)輝石巖礦物間的Mg同位素組成變化范圍也比較小,不同礦物之間Mg同位素分餾均達(dá)到平衡(圖4;Liuetal., 2011; Huetal., 2016; Zhaoetal., 2017b)。與之相比,華北含石榴子石的輝石巖的Mg同位素組成范圍變化較大(δ26/24Mg=-0.48‰~-0.10‰;圖3;Huetal., 2016; Zhaoetal., 2017b),其平均值(δ26/24Mg=-0.30±0.05‰,2SE;n=16)比硅酸鹽地球Mg同位素組成(δ26/24Mg=-0.25±0.07‰,2SD;n=139;Tengetal., 2010)稍微偏輕。相對(duì)于全巖,漢諾壩含石榴子石輝石巖中的單礦物的Mg同位素組成顯示了更高程度的不均一性(圖 4;Huetal., 2016; Zhaoetal., 2017b),其中橄欖石(平均值δ26/24Mg=-0.20±0.06‰,2SE;n=25)、斜方輝石(平均值δ26/24Mg=-0.20±0.06‰,2SE;n=19)和單斜輝石(δ26/24Mg=-0.20±0.05‰,2SE;n=24)之間Mg同位素組成接近,與世界其他地區(qū)的輝石巖中單礦物的Mg同位素組成在誤差范圍內(nèi)一致。另一方面,石榴子石具有較共存的橄欖石和輝石明顯偏輕且高度變化的Mg同位素組成(δ26/24Mg=-0.76‰~-0.37‰;圖4;Huetal., 2016; Zhaoetal., 2017b)。值得注意的是,華北地區(qū)單斜輝石巖及其單礦物具有顯著低于硅酸巖地球的Mg同位素組成(δ26/24Mg=-1.53‰~-1.20‰;圖3、圖4;Huetal., 2016, 2019; Zhaoetal., 2017b),這是目前已報(bào)道幔源捕虜體的最低Mg同位素值。

    圖5 華北地幔橄欖巖和輝石巖捕虜體Ca同位素組成地幔捕虜體數(shù)據(jù)來(lái)源:華北地幔橄欖巖及輝石巖Ca同位素?cái)?shù)據(jù)來(lái)自Kang et al.(2016),Zhao et al.(2017a),Chen et al.(2018)和Dai et al.(2020);華南地幔橄欖巖Ca同位素?cái)?shù)據(jù)來(lái)自Kang et al.(2016,2020);其他地區(qū)地幔橄欖巖和輝石巖數(shù)據(jù)來(lái)自Kang et al.(2017,2019),Ionov et al.(2019),Chen et al.(2019)和Zhu et al.(2021). 黑色虛線及其棕色條帶代表硅酸鹽地球平均值(δ44/40Ca=0.94±0.05‰,2SD;n=14;Kang et al., 2017). 圖6數(shù)據(jù)來(lái)源同此圖Fig.5 Ca isotopic composition of mantle peridotite and pyroxene xenoliths from the North China Craton

    圖6 華北地幔橄欖巖和輝石巖中共存單礦物的Ca同位素組成Fig.6 Ca isotopic compositions of coexisting minerals in mantle peridotites and pyroxene xenoliths from the North China Craton

    2.3 華北地幔捕虜體的Ca同位素特征

    鈣是地幔和地殼的主要組成元素之一。由于沉積碳酸鹽(平均值δ44/40Ca≈0.6‰;Fantle and Tipper, 2014及其參考文獻(xiàn))與硅酸鹽地球的Ca同位素組成(δ44/40Ca=0.94±0.05‰; 2SD; Kangetal., 2017; Chenetal., 2019)存在顯著差異,使得再循環(huán)碳酸鹽有可能顯著地改變地幔的Ca同位素組成,因此Ca同位素被認(rèn)為是示蹤地幔交代和深部碳循環(huán)的重要手段之一。目前,研究者對(duì)華北克拉通地區(qū)的地幔捕虜體進(jìn)行了Ca同位素分析,結(jié)果發(fā)現(xiàn)華北克拉通巖石圈地幔Ca同位素組成不均一(δ44/40Ca=-0.08‰~1.20‰;圖5、圖6;Kangetal., 2016; Zhaoetal., 2017b; Chenetal., 2018; Daietal., 2020)。其中,華北漢諾壩、繁峙和陽(yáng)原地區(qū)的鎂橄欖巖和輝石巖的Ca同位素組成(δ44/40Ca=0.76‰~1.20‰;平均值δ44/40Ca=0.92±0.03‰,2SE;n=39)與華南地區(qū)鎂橄欖巖(δ44/40Ca=0.85‰~1.14‰;平均值δ44/40Ca=0.98±0.04‰,2SE;n=14;Kangetal., 2016, 2020)以及世界其他地區(qū)的鎂橄欖巖和輝石巖的Ca同位素組成接近(δ44/40Ca=0.72‰~1.18‰;平均值δ44/40Ca=0.91±0.02‰,2SE;n=114;圖5;Kangetal., 2017, 2019; Chenetal., 2019; Ionovetal., 2019; Zhuetal., 2021)。但是華北陽(yáng)原地區(qū)的富Fe橄欖巖具有比華北、華南以及世界其他地區(qū)的鎂橄欖巖和輝石巖顯著偏輕且不均一的Ca 同位素組成(δ44/40Ca=-0.08‰~ 0.70‰;圖5;Zhaoetal., 2017b)。斜方輝石和單斜輝石是上地幔的兩個(gè)主要含Ca礦物,它們的Ca同位素組成對(duì)確定上地幔的Ca同位素組成和演化至關(guān)重要。華北地區(qū)鎂橄欖巖和輝石巖中的單斜輝石(δ44/40Ca=0.71‰~1.17‰;圖6)與其全巖具有相似的Ca同位素組成(δ44/40Ca=0.76‰~1.20‰),而斜方輝石具有比共生的單斜輝石系統(tǒng)偏重的Ca同位素組成和更大的變化范圍(δ44/40Ca=0.82‰~1.77‰;圖6;Kangetal., 2016; Zhaoetal., 2017b; Chenetal., 2018; Daietal., 2020),與華南地區(qū)(單斜輝石δ44/40Ca=0.83‰~1.14‰;斜方輝石δ44/40Ca=0.98‰~2.16‰;圖6;Kangetal., 2016, 2020)以及世界其他地區(qū)的鎂橄欖巖中單礦物的Ca同位素組成(單斜輝石δ44/40Ca=0.98‰~1.05‰;斜方輝石δ44/40Ca=1.40‰~ 1.73‰;圖6;Huangetal., 2010; Zhuetal., 2021)在誤差范圍內(nèi)一致。與全巖相同,華北陽(yáng)原地區(qū)的富Fe橄欖巖中單斜輝石和斜方輝石的Ca同位素組成(單斜輝石δ44/40Ca=-0.09‰~0.69‰;斜方輝石δ44/40Ca=-0.24‰~0.70‰;圖6;Zhaoetal., 2017b)明顯比華北、華南以及世界其他地區(qū)的鎂橄欖巖和輝石巖中的單礦物Ca同位素組成偏輕,其中樣品YY09-24中的斜方輝石δ44/40Ca低至-0.24‰,這是目前發(fā)現(xiàn)的地幔中最輕的Ca同位素組分。

    3 華北地幔捕擄體Fe-Mg-Ca同位素組成變化成因

    通過(guò)對(duì)比分析上述華北不同地區(qū)地幔捕虜體及其單礦物的Fe-Mg-Ca同位素組成特征,我們不難發(fā)現(xiàn)華北巖石圈地幔具有高度不均一的Fe-Mg-Ca同位素組成(Fe同位素有~2‰的變化,圖1、圖2;Mg同位素有~1.5‰的變化,圖3、圖4;Ca同位素有~2‰的變化,圖5、圖6)。是什么過(guò)程導(dǎo)致華北巖石圈地幔的Fe-Mg-Ca同位素組成不均一?可能的原因包括:(1)礦物結(jié)構(gòu)控制的礦物間同位素平衡分餾(Huangetal., 2010; Fengetal., 2014; Macrisetal., 2015; Youngetal., 2015; Wangetal., 2017; Antonellietal., 2019);(2)地幔部分熔融過(guò)程(Weyer and Ionov, 2007; Williams and Bizimis, 2014; Kangetal., 2017; Zhuetal., 2018);(3)橄欖巖-熔體反應(yīng)(Weyer and Ionov, 2007; Zhaoetal., 2010, 2012, 2015, 2017a, b; Poitrassonetal., 2013; Xiaoetal., 2013; Kangetal., 2016; Kangetal., 2020);(4)化學(xué)擴(kuò)散引起的動(dòng)力學(xué)分餾(Weyer and Ionov, 2007; Tengetal., 2011; Poitrassonetal., 2013; Kangetal., 2016; Dauphasetal., 2017; Antonelli and Simon, 2020)。下面我們將逐一總結(jié)和討論導(dǎo)致華北巖石圈地幔的Fe-Mg-Ca同位素組成不均一的可能因素。

    理論上,穩(wěn)定同位素平衡分餾程度的大小取決于該元素在不同礦物相中化學(xué)鍵的鍵能強(qiáng)度,而鍵強(qiáng)通常受到鍵長(zhǎng)和配位數(shù)(CN)的影響。通常鍵長(zhǎng)越短、配位數(shù)越小、鍵能越強(qiáng),而重的同位素會(huì)優(yōu)先富集在鍵能最強(qiáng)的位置(Bigeleisen and Mayer, 1947; Urey, 1947; Schauble, 2004; Youngetal., 2015)。由于橄欖石和斜方輝石中Fe2+均為6次配位,它們之間的Fe同位素分餾程度較?。欢鴨涡陛x石礦物結(jié)構(gòu)中既有Fe2+,又有少量Fe3+,因此單斜輝石的Fe-O鍵能相對(duì)較強(qiáng),理論上比橄欖石和斜方輝石略為富集Fe的重同位素;尖晶石中Fe2+以4次配位,F(xiàn)e3+以6次配位充填于八面體空隙中,而Fe2+在石榴子石中的配位數(shù)為8,因此尖晶石和石榴子分別具有顯著重和輕的Fe同位素組成。Macrisetal.(2015)采用離子半徑模型模擬計(jì)算相同溫度下地幔各個(gè)礦物間Fe同位素平衡分餾,計(jì)算結(jié)果進(jìn)一步支持δ57/54Fe尖晶石>δ57/54Fe單斜輝石>δ57/54Fe斜方輝石≥δ57/54Fe橄欖石>δ57/54Fe石榴子石。由于Mg2+和Fe2+的離子半徑相似,它們?cè)诘蒯5V物中通常占據(jù)相同的晶格位置,所以理論上Fe和Mg同位素的分餾方向是一致的。Huangetal.(2013)利用第一性原理計(jì)算了地幔各個(gè)礦物間Mg同位素平衡分餾表現(xiàn)為δ26/24Mg尖晶石>δ26/24Mg單斜輝石>δ26/24Mg斜方輝石≥δ26/24Mg橄欖石>δ26/24Mg石榴子石。單斜輝石和石榴子石中的Ca均是8次配位,但石榴子石中的Ca-O鍵長(zhǎng)(2.27~2.43?)短于單斜輝石中Ca-O鍵長(zhǎng)(2.45~2.50?;Magnaetal., 2015),因此理論上石榴子石相對(duì)單斜輝石富集重的Ca同位素。而斜方輝石和橄欖石中的中的Ca均是6次配位,重的Ca同位素會(huì)優(yōu)先富集在斜方輝石和橄欖石中。Fengetal.(2014)、Wangetal.(2017)和Huangetal.(2019)利用第一性原理計(jì)算了地幔各個(gè)礦物間Ca同位素平衡分餾表現(xiàn)為δ44/40Ca石榴子石>δ44/40Ca橄欖石≥δ44/40Ca斜方輝石>δ44/40Ca單斜輝石,這些理論計(jì)算工作對(duì)地幔深部過(guò)程如部分熔融、結(jié)晶分異和地幔熔體交代過(guò)程中Fe-Mg-Ca同位素分餾行為提供了重要的理論參考。

    巖石圈地幔通常被認(rèn)為是原始地幔經(jīng)歷不同熔融程度抽取玄武質(zhì)巖漿之后的殘留物,而后期可能受到了熔/流體的復(fù)雜改造,使得大陸巖石圈的地球化學(xué)組成高度不均一(Zhang, 2005)。學(xué)者們已經(jīng)對(duì)Fe-Mg-Ca同位素在地幔部分熔融過(guò)程中的分餾進(jìn)行了大量研究,不同的同位素分餾尺度不一。例如已有的研究顯示,洋中脊玄武巖(MORB)、洋島玄武巖(OIB)和島弧玄武巖(IAB)的Fe同位素組成平均值(δ57/54Fe=0.13±0.1‰,2SE;n=369;Johnsonetal., 2020)明顯重于地幔橄欖巖和輝石巖捕虜體的Fe同位素組成平均值(δ57/54Fe=0.03±0.03‰,2SE;n=233;Johnsonetal., 2020及其參考文獻(xiàn)),表明地幔部分熔融過(guò)程中Fe同位素發(fā)生了分餾。這一分餾是由于Fe是典型的變價(jià)元素(Fe2+、Fe3+),在地幔橄欖巖部分熔融過(guò)程中,F(xiàn)e3+比Fe2+更不相容,F(xiàn)e3+優(yōu)先進(jìn)入熔體相,使得熔體相富集Fe3+,難熔的殘留相富集Fe2+,而Fe3+的化學(xué)鍵較強(qiáng),傾向于富集重同位素(Dauphasetal., 2017及其參考文獻(xiàn))。但是由于在地幔部分熔融過(guò)程中,體系中絕大多數(shù)Fe還保留在地幔殘留相中,相對(duì)于玄武質(zhì)熔體,部分熔融作用對(duì)地幔源區(qū)殘留相的Fe同位素影響有限(Fe同位素分餾不到0.15‰;Weyer and Ionov, 2007; Dauphasetal., 2009; Poitrassonetal., 2013)。與Fe同位素不同,前人研究結(jié)果表明,來(lái)自世界不同地區(qū)的地幔橄欖巖捕虜體及其單礦物的Mg同位素組成變化范圍較小(δ26/24Mg=-0.44‰~0.00‰;圖3、圖4),其平均值(δ26/24Mg=-0.24±0.03‰,2SE;n=172)與球粒隕石(δ26/24Mg=-0.28±0.01‰,2SE;n=38;Tengetal., 2010)以及大洋玄武巖具有相似的Mg同位素組成(δ26/24Mg=-0.25±0.01‰,2SE;n=110;Tengetal., 2010),表明地幔部分熔融及巖漿結(jié)晶分異過(guò)程沒(méi)有發(fā)生顯著的Mg同位素分餾。Kangetal.(2017)研究發(fā)現(xiàn),中高程度虧損(>20%)的地幔橄欖巖的Ca同位素組成較硅酸鹽地球的值略微偏重0.1‰,并且大洋玄武巖的Ca同位素較硅酸鹽地球的值略微偏輕 ~0.15‰(Zhuetal., 2018)。部分熔融過(guò)程模擬結(jié)果也表明,在沒(méi)有石榴子石參與的情況下,部分熔融過(guò)程不會(huì)造成明顯可觀測(cè)的Ca同位素組分變化(<0.1‰);而在石榴子石參與的情況下,會(huì)使得Ca同位素分餾值達(dá)到0.15‰(Zhangetal., 2018; Chenetal., 2019; Lietal., 2022)。而華北地幔捕虜體Fe同位素分餾達(dá)到~2‰(圖1、圖2),Mg同位素分餾~1.5‰(圖3、圖4),Ca同位素分餾~2‰(圖5、圖6),因此很難通過(guò)地幔部分熔融過(guò)程來(lái)解釋華北地幔捕擄體Fe-Mg-Ca同位素組成變化。

    圖7 華北陽(yáng)原地區(qū)富鐵橄欖巖中含交代礦物金云母(Phl)和角閃石(Amp)及交代熔蝕結(jié)構(gòu)Ol-橄欖石;Opx-斜方輝石;Cpx-單斜輝石;Sp-尖晶石 ;Spong rim-篩狀邊Fig.7 Fe-rich lherzolites bearing-phlogopite (Phl) and -amphibole (Amp) (plane-polarized light) with resorbed orthopyroxene relics inside newly formed sieve-textured clinopyroxene

    地幔交代作用(橄欖巖-熔體反應(yīng))不僅是造成巖石圈地幔組成不均一性的重要方式,而且是造成巖石圈地幔組成轉(zhuǎn)變的重要原因(Zhang, 2005)。根據(jù)目前累積的數(shù)據(jù)來(lái)看,華北地幔捕虜體Fe-Mg-Ca同位素組成變化與地幔交代作用密不可分。例如華北漢諾壩、三義堂、陽(yáng)原、繁峙、鶴壁以及昌樂(lè)地區(qū)沒(méi)有受到地幔熔體交代的鎂橄欖巖,具有一致的Mg#(~90),單斜輝石具有低的稀土元素總量和虧損或者平坦的稀土元素配分型式,并具有虧損的Sr-Nd同位素組成,可能代表了低程度部分熔融的殘留巖石圈地?;蝾?lèi)似于原始地幔橄欖巖中的新增生的巖石圈地幔,而沒(méi)有受到后期熔/流體的交代作用影響,其Fe-Mg-Ca同位素組成比較均一(趙新苗,2019及其參考文獻(xiàn)),平均值與硅酸鹽地球的Fe-Mg-Ca同位素組成在誤差范圍內(nèi)一致(δ57/54Fe=0.04±0.04‰;δ26/24Mg=-0.25±0.07‰;δ44/40Ca=0.94±0.05‰;Weyer and Ionov, 2007; Tengetal., 2010; Kangetal., 2017)。同時(shí)這類(lèi)橄欖巖的不同共生礦物之間Fe-Mg-Ca同位素分餾達(dá)到了平衡,主要受晶體結(jié)構(gòu)(配位數(shù))控制。而經(jīng)受不同類(lèi)型地幔交代熔體影響的地幔橄欖巖,F(xiàn)e-Mg-Ca同位素組成存在明顯不同(趙新苗,2019及其參考文獻(xiàn))。例如上述地區(qū)的部分鎂橄欖巖的Fe-Mg-Ca同位素組成(δ57/54Fe=-0.28‰~0.19‰;δ26/24Mg=-0.39‰~-0.19‰;δ44/40Ca=0.76‰~1.20‰;圖1、圖3、圖5)相對(duì)沒(méi)有受到地幔熔體交代的鎂橄欖巖具有更大的變化范圍。這類(lèi)橄欖巖通常被認(rèn)為是原始地幔經(jīng)不同熔融程度抽取玄武質(zhì)巖漿之后的殘留物,經(jīng)歷了后期不同程度,不同期次的小規(guī)模熔/流體交代作用的改造,使得橄欖巖的微量元素變富集(如富集輕稀土LREE),甚至伴隨有金云母、角閃石、磷灰石等交代礦物的出現(xiàn)(趙新苗,2019及其參考文獻(xiàn))。然而經(jīng)受這種類(lèi)型地幔熔體交代的地幔橄欖巖,F(xiàn)e-Mg-Ca同位素組成變化范圍有限,它們與其主量(Mg#)、微量(單斜輝石(La/Yb)N)以及相互之間均沒(méi)有明顯的相關(guān)性(趙新苗,2019及其參考文獻(xiàn))。同時(shí)這些鎂橄欖巖不同礦物間Fe-Mg-Ca同位素分餾有限,礦物間的Fe-Mg-Ca同位素分餾均與同位素平衡分餾理論計(jì)算值接近。因此,盡管華北有些地區(qū)的鎂橄欖巖經(jīng)歷了小規(guī)模的硅酸鹽熔體交代(即熔/巖比較小),但是僅僅造成這些鎂橄欖巖微量元素富集,而沒(méi)有導(dǎo)致其Fe-Mg-Ca同位素組成顯著變化,礦物間Fe-Mg-Ca同位素分餾也達(dá)到了平衡(趙新苗,2019及其參考文獻(xiàn))。而受到強(qiáng)烈地幔熔體交代的橄欖巖(例如華北地區(qū)的富Fe橄欖巖和異剝橄欖巖)和輝石巖(例如漢諾壩石榴子石輝石巖和單斜輝石巖)具有與硅酸鹽地球明顯不同的Fe-Mg-Ca同位素組成(δ57/54Fe=-1.04‰~1.00‰;δ26/24Mg=-1.54‰~0.06‰;δ44/40Ca=-0.08‰~0.70‰;圖1、圖3、圖5)并且礦物間Fe-Mg-Ca同位素分餾均未完全達(dá)到平衡。華北這些富Fe橄欖巖和異剝橄欖巖富Fe、具有低的Mg#、富含交代礦物金云母、角閃石和磷灰石,其礦物顆粒內(nèi)部存在單斜輝石取代斜方輝石的現(xiàn)象(圖7),為橄欖巖-熔體反應(yīng)提供了明確的巖石學(xué)證據(jù)。更為重要的是,這些華北地區(qū)的富Fe橄欖巖和異剝橄欖巖的Fe同位素與其Mg#和Ca同位素呈很好的正相關(guān)性,而與其Mg同位素呈反相關(guān)關(guān)系(趙新苗,2019及其參考文獻(xiàn)),表明橄欖巖-熔體反應(yīng)是導(dǎo)致巖石圈地幔Fe-Mg-Ca同位素組成高度不均一及地幔礦物之間存在顯著Fe-Mg-Ca同位素分餾的重要機(jī)制。通過(guò)二端元混合及化學(xué)擴(kuò)散理論模擬計(jì)算(Zhaoetal., 2017a),筆者證實(shí)了華北陽(yáng)原地區(qū)富Fe橄欖巖低的δ44/40Ca和δ57/54Fe值以及礦物間的Fe-Ca同位素不平衡分餾現(xiàn)象是橄欖巖與富Fe的硅酸鹽熔體反應(yīng)過(guò)程中兩者之間的元素交換所引起的動(dòng)力學(xué)分餾的體現(xiàn),而不是碳酸鹽熔體交代的結(jié)果。該研究實(shí)例表明橄欖巖輕的Ca同位素組成不一定是碳酸鹽熔體交代的直接表現(xiàn),這更新了長(zhǎng)期以來(lái)人們對(duì)硅酸鹽Ca同位素組成變化范圍較小的認(rèn)知,為認(rèn)識(shí)地幔Ca同位素組成不均一性以及更好地用Ca同位素示蹤深部碳再循環(huán)提供了重要的約束條件。Chenetal.(2018)報(bào)道了華北繁峙地區(qū)橄欖巖的Ca同位素組成,認(rèn)為碳酸鹽化沉積物的俯沖可能會(huì)導(dǎo)致地幔Ca同位素組成的不均一,而Ionovetal.(2019)認(rèn)為現(xiàn)階段單純使用Ca同位素很難對(duì)硅酸鹽熔體與碳酸鹽熔體進(jìn)行區(qū)分。Huetal.(2016)和Zhaoetal.(2017b)研究發(fā)現(xiàn)石榴子石輝石巖的Mg和Fe同位素組成變化范圍很大,并且存在負(fù)相關(guān)性,很可能是地幔橄欖巖與俯沖大洋板片釋放熔體反應(yīng)的產(chǎn)物。值得注意的是,黑色單斜輝石捕虜體具有非常輕的Mg同位素組成(δ26/24Mg=-1.53‰~-1.04‰;圖3),被認(rèn)為是再循環(huán)碳酸鹽熔體在地幔滲透的結(jié)果(Huetal., 2016, 2019)。上述研究結(jié)果表明,華北地幔捕虜體的Fe-Mg-Ca同位素組成是不均一的,橄欖巖-熔體反應(yīng)可以顯著改變巖石圈地幔的Fe-Mg-Ca同位素提供了新的地球化學(xué)佐證。

    4 總結(jié)與展望

    金屬穩(wěn)定同位素地球化學(xué)是當(dāng)前同位素地球化學(xué)研究的前沿領(lǐng)域,我國(guó)學(xué)者在華北地幔捕虜體Fe,Mg和Ca同位素研究方面做出了長(zhǎng)足的進(jìn)展,證實(shí)Fe、Ca和Mg等金屬穩(wěn)定同位素具有研究地幔不均一性和示蹤大陸巖石圈地幔改造過(guò)程中熔體的性質(zhì)和來(lái)源的潛力(Yangetal., 2009; Tengetal., 2010; Zhaoetal., 2010, 2012, 2015, 2017a, b; Huangetal., 2011; Liuetal., 2011; Xiaoetal., 2013; Huetal., 2016, 2019; Kangetal., 2016; Chenetal., 2018; Daietal., 2020)。然而,現(xiàn)有的研究工作主要是基于華北克拉通個(gè)別地區(qū)的地幔橄欖巖和輝石巖捕虜體的Fe、Ca和Mg同位素?cái)?shù)據(jù),雖然結(jié)合了傳統(tǒng)的巖石學(xué)、元素和同位素地球化學(xué)信息的相互制約和驗(yàn)證,獲得了一些新的認(rèn)識(shí),但是目前獲得的認(rèn)識(shí)尚不夠深入,仍有一些關(guān)鍵科學(xué)問(wèn)題亟待回答。主要表現(xiàn)在:

    (1)對(duì)華北巖石圈地幔的Fe、Ca和Mg同位素組成特征的了解還是局部的、零碎的,尚不能較為全面細(xì)致地揭示華北巖石圈地幔Fe、Ca和Mg同位素組成的時(shí)空演化規(guī)律;

    (2)以往的研究多聚焦于單一的Fe、Ca或Mg同位素研究,對(duì)示蹤地幔改造過(guò)程中熔體性質(zhì)和來(lái)源存在多解性和局限性?,F(xiàn)有的研究已觀察到巖石圈地幔Fe、Ca和Mg同位素組成的變化與橄欖巖-熔體反應(yīng)密切相關(guān),并且這三個(gè)同位素體系之間存在一定的關(guān)聯(lián),它們的聯(lián)合示蹤可以相互檢驗(yàn),這有利于甄別巖石圈地幔演化和改造過(guò)程熔體的性質(zhì)和來(lái)源,為正確認(rèn)識(shí)克拉通破壞的機(jī)制和過(guò)程提供必要條件,但這方面的理論基礎(chǔ)尚不夠完善;

    (3)導(dǎo)致巖石圈地幔組成轉(zhuǎn)變的熔體的來(lái)源和性質(zhì)復(fù)雜多樣,F(xiàn)e、Ca和Mg同位素在不同性質(zhì)熔體遷移過(guò)程中的分餾尺度和具體機(jī)制尚不完全明確。這些目前迫切需要解決的科學(xué)問(wèn)題客觀上要求我們開(kāi)展更大范圍的采樣和更深入的對(duì)比研究,以期對(duì)華北巖石圈地幔的Fe、Ca和Mg同位素組成特征及時(shí)空演化規(guī)律有更全面和清晰的認(rèn)識(shí);進(jìn)一步深化對(duì)Fe、Ca和Mg同位素在地幔熔體交代過(guò)程中的分餾機(jī)制及相互之間的聯(lián)系的認(rèn)識(shí);構(gòu)建Fe-Ca-Mg同位素聯(lián)合示蹤地幔改造過(guò)程中熔體的性質(zhì)和來(lái)源的理論框架。在此基礎(chǔ)上,實(shí)現(xiàn)對(duì)大陸巖石圈地幔改造過(guò)程中熔體的性質(zhì)和來(lái)源的準(zhǔn)確示蹤,為更深層次的制約巖石圈地幔屬性的轉(zhuǎn)變過(guò)程和機(jī)理提供重要依據(jù)。

    謹(jǐn)以此文祝賀周新華老師八十華誕和感恩一直以來(lái)對(duì)第一作者的指導(dǎo)和教誨!

    致謝感謝評(píng)審人對(duì)本文提出的建設(shè)性修改意見(jiàn)。

    猜你喜歡
    橄欖巖巖石圈輝石
    蛇綠巖中識(shí)別出不同類(lèi)型的方輝橄欖巖及其巖相分帶
    ——來(lái)自丁青蛇綠巖專(zhuān)項(xiàng)地質(zhì)調(diào)查的證據(jù)*
    第四章 堅(jiān)硬的巖石圈
    新疆和巴西薔薇輝石玉的對(duì)比研究
    粉色薔薇輝石的寶石學(xué)及礦物學(xué)特征
    河南西峽縣洋淇溝橄欖巖礦床地質(zhì)特征及開(kāi)發(fā)利用
    不同溫度、壓強(qiáng)、氧逸度條件下斜方輝石含水性的實(shí)驗(yàn)研究
    西藏羅布莎地幔橄欖巖礦物學(xué)初探
    巖石圈磁場(chǎng)異常變化與巖石圈結(jié)構(gòu)的關(guān)系
    地震研究(2017年3期)2017-11-06 21:54:14
    2014年魯?shù)?—5級(jí)地震相關(guān)斷裂的巖石圈磁異常分析
    地震研究(2017年3期)2017-11-06 01:58:51
    利用石榴橄欖巖重建大陸俯沖帶的古動(dòng)力學(xué)環(huán)境及其演化過(guò)程*
    国产单亲对白刺激| 操出白浆在线播放| 日韩高清综合在线| 在线看a的网站| 亚洲中文av在线| 欧美老熟妇乱子伦牲交| 日本 av在线| 国产xxxxx性猛交| 亚洲国产精品999在线| 伊人久久大香线蕉亚洲五| 在线永久观看黄色视频| 这个男人来自地球电影免费观看| 国产亚洲欧美98| 欧美最黄视频在线播放免费 | 亚洲精品一卡2卡三卡4卡5卡| 韩国精品一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 国产高清视频在线播放一区| av欧美777| 99国产综合亚洲精品| 两性夫妻黄色片| 免费av中文字幕在线| 女人爽到高潮嗷嗷叫在线视频| 身体一侧抽搐| 岛国视频午夜一区免费看| a级毛片在线看网站| 精品一区二区三区四区五区乱码| 久久 成人 亚洲| 欧美成人午夜精品| 精品乱码久久久久久99久播| 18禁裸乳无遮挡免费网站照片 | 亚洲国产精品一区二区三区在线| 精品久久久久久电影网| 亚洲av片天天在线观看| 人人妻,人人澡人人爽秒播| 亚洲精品美女久久久久99蜜臀| 18美女黄网站色大片免费观看| svipshipincom国产片| xxxhd国产人妻xxx| 午夜视频精品福利| 亚洲精品国产精品久久久不卡| 欧美成狂野欧美在线观看| www.熟女人妻精品国产| 精品国产亚洲在线| 欧美激情 高清一区二区三区| а√天堂www在线а√下载| 久久这里只有精品19| 一级毛片女人18水好多| 亚洲国产欧美日韩在线播放| 99国产精品免费福利视频| 一边摸一边抽搐一进一小说| 伦理电影免费视频| 在线观看www视频免费| 日韩欧美一区二区三区在线观看| 电影成人av| 99久久综合精品五月天人人| 激情在线观看视频在线高清| 精品人妻在线不人妻| 亚洲一码二码三码区别大吗| 欧美日韩亚洲国产一区二区在线观看| 成年版毛片免费区| 久久精品国产亚洲av香蕉五月| 最近最新中文字幕大全免费视频| 老司机午夜福利在线观看视频| 91在线观看av| 日韩 欧美 亚洲 中文字幕| 日本免费a在线| 性少妇av在线| 女性被躁到高潮视频| 久久国产乱子伦精品免费另类| 欧美黑人欧美精品刺激| 青草久久国产| av国产精品久久久久影院| 在线观看www视频免费| 日本精品一区二区三区蜜桃| 国产av又大| 久久人人爽av亚洲精品天堂| 亚洲成人久久性| 中文字幕人妻丝袜一区二区| 99久久国产精品久久久| 午夜福利影视在线免费观看| 日本欧美视频一区| 国产成人精品无人区| 国产精品亚洲av一区麻豆| 丝袜美足系列| 在线观看一区二区三区激情| 90打野战视频偷拍视频| 成年女人毛片免费观看观看9| 午夜福利影视在线免费观看| 日韩精品免费视频一区二区三区| 日韩高清综合在线| 久久 成人 亚洲| av中文乱码字幕在线| 久久精品人人爽人人爽视色| 欧美中文日本在线观看视频| av免费在线观看网站| av在线播放免费不卡| 国产成人精品无人区| svipshipincom国产片| 精品少妇一区二区三区视频日本电影| 欧美日本中文国产一区发布| 精品一区二区三区四区五区乱码| 好看av亚洲va欧美ⅴa在| 国产亚洲精品第一综合不卡| 69av精品久久久久久| 亚洲男人的天堂狠狠| 亚洲精品一卡2卡三卡4卡5卡| 他把我摸到了高潮在线观看| 国产精品美女特级片免费视频播放器 | 欧美人与性动交α欧美软件| 日韩欧美一区二区三区在线观看| 一级作爱视频免费观看| 国产精品久久久久成人av| 男人舔女人下体高潮全视频| av网站免费在线观看视频| √禁漫天堂资源中文www| 又紧又爽又黄一区二区| 国产精品九九99| 嫁个100分男人电影在线观看| 精品久久久久久电影网| 中文字幕人妻熟女乱码| 久久久久国内视频| 1024视频免费在线观看| 国产精品秋霞免费鲁丝片| 久久精品国产99精品国产亚洲性色 | av视频免费观看在线观看| 国产欧美日韩一区二区三| 国产三级黄色录像| 69av精品久久久久久| 精品熟女少妇八av免费久了| www.精华液| 丁香六月欧美| 人人妻人人添人人爽欧美一区卜| 免费一级毛片在线播放高清视频 | 亚洲午夜理论影院| 又黄又粗又硬又大视频| 在线视频色国产色| 色综合站精品国产| 叶爱在线成人免费视频播放| 亚洲aⅴ乱码一区二区在线播放 | 麻豆一二三区av精品| 又紧又爽又黄一区二区| 一级a爱片免费观看的视频| 日韩三级视频一区二区三区| 久久久久久久久久久久大奶| 日本vs欧美在线观看视频| 高潮久久久久久久久久久不卡| 色老头精品视频在线观看| 国产av在哪里看| 99国产精品免费福利视频| 亚洲精华国产精华精| 性欧美人与动物交配| 美女福利国产在线| 国产无遮挡羞羞视频在线观看| 黄色女人牲交| 久久久久久久午夜电影 | 免费观看精品视频网站| 女性被躁到高潮视频| 搡老乐熟女国产| 免费一级毛片在线播放高清视频 | 欧美+亚洲+日韩+国产| 他把我摸到了高潮在线观看| 男女下面插进去视频免费观看| 在线av久久热| 亚洲av片天天在线观看| 丰满饥渴人妻一区二区三| 国产高清激情床上av| 亚洲精品粉嫩美女一区| 他把我摸到了高潮在线观看| 精品国产乱码久久久久久男人| 麻豆av在线久日| 人人妻人人澡人人看| 日本vs欧美在线观看视频| 婷婷丁香在线五月| 丰满迷人的少妇在线观看| 国产精品一区二区免费欧美| 91在线观看av| 激情视频va一区二区三区| 99久久人妻综合| 亚洲精品一二三| 国产精品久久久久久人妻精品电影| 亚洲五月婷婷丁香| 99久久综合精品五月天人人| 国产亚洲精品一区二区www| 午夜视频精品福利| 亚洲久久久国产精品| 亚洲人成电影观看| 亚洲国产中文字幕在线视频| 无限看片的www在线观看| av天堂在线播放| 视频区欧美日本亚洲| 成人黄色视频免费在线看| 成人18禁高潮啪啪吃奶动态图| 国产精品免费视频内射| 国产不卡一卡二| 狠狠狠狠99中文字幕| 一级片'在线观看视频| 一区二区三区国产精品乱码| 欧美精品亚洲一区二区| 亚洲成人免费av在线播放| 午夜福利免费观看在线| 69精品国产乱码久久久| 久久人妻av系列| 久热爱精品视频在线9| 久热这里只有精品99| 国产成人精品久久二区二区91| 两人在一起打扑克的视频| 一区二区三区国产精品乱码| 日韩免费av在线播放| 亚洲熟妇中文字幕五十中出 | 欧美成人午夜精品| 亚洲视频免费观看视频| 18美女黄网站色大片免费观看| svipshipincom国产片| 老司机深夜福利视频在线观看| 国产精品久久久av美女十八| 亚洲av第一区精品v没综合| 久久久国产成人免费| 在线播放国产精品三级| 久久国产精品人妻蜜桃| 国产亚洲精品一区二区www| 国产三级在线视频| 亚洲熟妇熟女久久| 母亲3免费完整高清在线观看| 精品久久久久久电影网| 日韩大尺度精品在线看网址 | 热re99久久国产66热| 一本大道久久a久久精品| 亚洲五月色婷婷综合| 97碰自拍视频| 成年人免费黄色播放视频| 在线观看午夜福利视频| 黑人猛操日本美女一级片| 在线观看免费高清a一片| 久久婷婷成人综合色麻豆| 久久狼人影院| 一夜夜www| 久久草成人影院| 9色porny在线观看| 亚洲精品在线美女| 国产一卡二卡三卡精品| 男人的好看免费观看在线视频 | 中文字幕高清在线视频| 18禁观看日本| 国产精品爽爽va在线观看网站 | 欧美+亚洲+日韩+国产| 日韩视频一区二区在线观看| 成年女人毛片免费观看观看9| 一本大道久久a久久精品| 国产又爽黄色视频| 久久久久久久精品吃奶| 黄色丝袜av网址大全| 国产精品免费视频内射| 99热只有精品国产| 操出白浆在线播放| 美国免费a级毛片| 亚洲国产毛片av蜜桃av| 亚洲av熟女| 夫妻午夜视频| 麻豆成人av在线观看| 免费高清视频大片| 免费不卡黄色视频| 9热在线视频观看99| 久久精品国产清高在天天线| 久久久国产成人免费| 国产1区2区3区精品| 99re在线观看精品视频| 99riav亚洲国产免费| 国产成人av激情在线播放| 自线自在国产av| 亚洲九九香蕉| √禁漫天堂资源中文www| 侵犯人妻中文字幕一二三四区| 午夜免费观看网址| 黄片大片在线免费观看| 午夜免费鲁丝| 黄色 视频免费看| 亚洲中文av在线| 日本精品一区二区三区蜜桃| 啦啦啦 在线观看视频| 亚洲国产精品一区二区三区在线| 精品久久久久久成人av| 无限看片的www在线观看| 午夜福利欧美成人| 18禁美女被吸乳视频| 欧美成人免费av一区二区三区| 成年版毛片免费区| 欧美老熟妇乱子伦牲交| av福利片在线| 韩国av一区二区三区四区| 在线观看一区二区三区激情| 交换朋友夫妻互换小说| 搡老乐熟女国产| 国产成人欧美| 免费高清在线观看日韩| 91麻豆av在线| 欧美成人午夜精品| 黄色视频不卡| 少妇的丰满在线观看| 国产精品国产高清国产av| 亚洲久久久国产精品| 亚洲avbb在线观看| 美女午夜性视频免费| 日本wwww免费看| 亚洲精品在线美女| 久久久久国内视频| 亚洲人成网站在线播放欧美日韩| 国产精品 国内视频| 国产免费男女视频| 国产精品av久久久久免费| 精品国内亚洲2022精品成人| 国产精品久久久人人做人人爽| 久久久国产成人精品二区 | 午夜老司机福利片| 热re99久久精品国产66热6| 亚洲国产欧美一区二区综合| 国产精品国产高清国产av| 国内毛片毛片毛片毛片毛片| 9191精品国产免费久久| 91成年电影在线观看| 人妻丰满熟妇av一区二区三区| 久久精品aⅴ一区二区三区四区| 性欧美人与动物交配| www.自偷自拍.com| 热99re8久久精品国产| 热re99久久精品国产66热6| 久久精品91无色码中文字幕| 亚洲午夜精品一区,二区,三区| 极品教师在线免费播放| 亚洲av片天天在线观看| 国产一区二区三区在线臀色熟女 | 欧美在线一区亚洲| 久久久久国产精品人妻aⅴ院| 久久中文字幕人妻熟女| 精品卡一卡二卡四卡免费| 欧美日本中文国产一区发布| 看免费av毛片| 欧美av亚洲av综合av国产av| 十八禁人妻一区二区| 制服人妻中文乱码| 欧美激情久久久久久爽电影 | 女同久久另类99精品国产91| 亚洲av美国av| 国产麻豆69| 成年版毛片免费区| 天堂√8在线中文| 国产区一区二久久| 一二三四在线观看免费中文在| 亚洲久久久国产精品| 国产97色在线日韩免费| 欧美日韩精品网址| 日本黄色日本黄色录像| 欧美黄色淫秽网站| 午夜两性在线视频| 日韩视频一区二区在线观看| 国产高清videossex| 日韩欧美国产一区二区入口| 18美女黄网站色大片免费观看| 91精品国产国语对白视频| 免费在线观看完整版高清| 午夜成年电影在线免费观看| 国产av一区二区精品久久| 精品午夜福利视频在线观看一区| 啦啦啦免费观看视频1| 亚洲免费av在线视频| 男人舔女人的私密视频| 久久久精品国产亚洲av高清涩受| 国产精品影院久久| 国产亚洲欧美98| av中文乱码字幕在线| 国产麻豆69| 日韩精品青青久久久久久| 中文字幕人妻熟女乱码| 大型黄色视频在线免费观看| 精品无人区乱码1区二区| 老熟妇仑乱视频hdxx| 国产av又大| 国产无遮挡羞羞视频在线观看| 男人的好看免费观看在线视频 | 日韩一卡2卡3卡4卡2021年| 亚洲 国产 在线| 18禁黄网站禁片午夜丰满| 老熟妇仑乱视频hdxx| 久久香蕉激情| www.自偷自拍.com| av片东京热男人的天堂| 精品国产一区二区三区四区第35| 高清av免费在线| 精品国产一区二区三区四区第35| 夫妻午夜视频| 亚洲第一青青草原| xxxhd国产人妻xxx| 久久久国产精品麻豆| 丝袜美足系列| 精品国产乱码久久久久久男人| 香蕉久久夜色| 久久人妻熟女aⅴ| 高清欧美精品videossex| 性色av乱码一区二区三区2| 欧美人与性动交α欧美软件| 国产成人精品无人区| 日韩精品中文字幕看吧| 麻豆av在线久日| 国产亚洲av高清不卡| 两性夫妻黄色片| 色哟哟哟哟哟哟| 欧美日韩视频精品一区| 欧美+亚洲+日韩+国产| 操美女的视频在线观看| 精品一区二区三卡| 国产人伦9x9x在线观看| 日韩三级视频一区二区三区| 在线观看免费高清a一片| 极品教师在线免费播放| 美女国产高潮福利片在线看| 1024视频免费在线观看| 黄色 视频免费看| 国产成人一区二区三区免费视频网站| 久久国产乱子伦精品免费另类| 色综合欧美亚洲国产小说| 两个人免费观看高清视频| 最近最新中文字幕大全电影3 | 大型黄色视频在线免费观看| 久久久国产成人精品二区 | 一级毛片女人18水好多| 免费高清在线观看日韩| 手机成人av网站| 一边摸一边抽搐一进一小说| 一进一出抽搐gif免费好疼 | 国产精品免费一区二区三区在线| 国产成人免费无遮挡视频| av在线天堂中文字幕 | 久久中文字幕人妻熟女| 欧美日韩精品网址| 亚洲色图 男人天堂 中文字幕| 日本免费一区二区三区高清不卡 | 美国免费a级毛片| 久久久久久久精品吃奶| 亚洲精品在线美女| aaaaa片日本免费| 在线观看一区二区三区激情| 成人影院久久| 亚洲中文av在线| 久久精品亚洲精品国产色婷小说| 麻豆av在线久日| 欧美亚洲日本最大视频资源| 免费在线观看日本一区| 久久中文字幕一级| 国产精品免费视频内射| av福利片在线| 天天影视国产精品| 久久热在线av| 亚洲av日韩精品久久久久久密| 无限看片的www在线观看| 电影成人av| 超碰97精品在线观看| 国产精品秋霞免费鲁丝片| 亚洲第一av免费看| 最新美女视频免费是黄的| 国产三级黄色录像| 曰老女人黄片| 午夜免费激情av| 一进一出抽搐gif免费好疼 | 88av欧美| 窝窝影院91人妻| 精品国产乱码久久久久久男人| 成年版毛片免费区| 麻豆国产av国片精品| 国产深夜福利视频在线观看| 一级片'在线观看视频| 亚洲激情在线av| 久久 成人 亚洲| 51午夜福利影视在线观看| 国产视频一区二区在线看| 国产精品美女特级片免费视频播放器 | av在线播放免费不卡| 国产日韩一区二区三区精品不卡| 久久久久久久午夜电影 | 午夜精品在线福利| 日本五十路高清| 夫妻午夜视频| 琪琪午夜伦伦电影理论片6080| 免费日韩欧美在线观看| 香蕉国产在线看| 国产精品影院久久| 成人18禁在线播放| 国产精品九九99| 亚洲va日本ⅴa欧美va伊人久久| 国产成人影院久久av| www.自偷自拍.com| 多毛熟女@视频| 天堂√8在线中文| 韩国av一区二区三区四区| 久久亚洲真实| 999久久久精品免费观看国产| 操美女的视频在线观看| 亚洲视频免费观看视频| 亚洲午夜理论影院| 在线观看午夜福利视频| 纯流量卡能插随身wifi吗| 性欧美人与动物交配| 日韩欧美一区视频在线观看| 日日爽夜夜爽网站| 日韩av在线大香蕉| 超碰97精品在线观看| 视频区图区小说| 日本五十路高清| 国内久久婷婷六月综合欲色啪| av有码第一页| 国产av又大| 亚洲精品一区av在线观看| 国产成人欧美在线观看| 99re在线观看精品视频| 成人18禁高潮啪啪吃奶动态图| 国产一区二区三区在线臀色熟女 | 午夜福利在线观看吧| 国产精品98久久久久久宅男小说| 精品熟女少妇八av免费久了| 啦啦啦免费观看视频1| 亚洲国产中文字幕在线视频| 久久精品影院6| 久久精品亚洲熟妇少妇任你| 欧美乱妇无乱码| 国产欧美日韩一区二区三区在线| 精品久久久久久成人av| 国产成人系列免费观看| 丰满人妻熟妇乱又伦精品不卡| 久久久国产精品麻豆| av免费在线观看网站| 热99re8久久精品国产| 亚洲成人免费av在线播放| 久久国产精品人妻蜜桃| 国产av精品麻豆| 老司机深夜福利视频在线观看| 香蕉国产在线看| 天天影视国产精品| 男女高潮啪啪啪动态图| 老汉色av国产亚洲站长工具| 日韩国内少妇激情av| 男人操女人黄网站| 国产深夜福利视频在线观看| 天堂影院成人在线观看| 香蕉丝袜av| 国产精品国产高清国产av| 免费在线观看黄色视频的| 一区二区日韩欧美中文字幕| 丁香六月欧美| 日本免费a在线| 国产成人欧美| 午夜久久久在线观看| 另类亚洲欧美激情| 香蕉久久夜色| 真人一进一出gif抽搐免费| 日韩成人在线观看一区二区三区| 午夜成年电影在线免费观看| 亚洲精品美女久久av网站| 久久 成人 亚洲| 亚洲精品美女久久av网站| 精品国内亚洲2022精品成人| 国产xxxxx性猛交| 热99re8久久精品国产| 悠悠久久av| 日韩大码丰满熟妇| 亚洲熟女毛片儿| 日韩av在线大香蕉| 两性午夜刺激爽爽歪歪视频在线观看 | 久久人妻av系列| 99热只有精品国产| 黄色片一级片一级黄色片| 亚洲精品在线观看二区| 亚洲欧美激情综合另类| 久久婷婷成人综合色麻豆| 午夜福利在线观看吧| 少妇 在线观看| 亚洲黑人精品在线| 女同久久另类99精品国产91| 超色免费av| 久久精品亚洲熟妇少妇任你| 国产成人精品久久二区二区91| 久久久久久久久中文| 国产一区二区三区在线臀色熟女 | 日韩三级视频一区二区三区| 欧美不卡视频在线免费观看 | 欧美最黄视频在线播放免费 | 久久精品亚洲精品国产色婷小说| 日日干狠狠操夜夜爽| 嫩草影视91久久| xxx96com| 久久久久久久午夜电影 | 一进一出抽搐动态| 亚洲精品国产区一区二| 国产成人啪精品午夜网站| 亚洲av电影在线进入| 久久久久久免费高清国产稀缺| 美女 人体艺术 gogo| 美女高潮喷水抽搐中文字幕| 一边摸一边抽搐一进一出视频| 交换朋友夫妻互换小说| 欧美日韩亚洲高清精品| 正在播放国产对白刺激| 亚洲欧美一区二区三区黑人| 欧美激情久久久久久爽电影 | 国产成人精品久久二区二区免费| 丝袜美腿诱惑在线| 少妇粗大呻吟视频| 国产高清国产精品国产三级| 久久精品人人爽人人爽视色| 精品一区二区三区视频在线观看免费 | 日日摸夜夜添夜夜添小说| 久久精品91无色码中文字幕| 又紧又爽又黄一区二区| 欧美日韩福利视频一区二区| 巨乳人妻的诱惑在线观看| av国产精品久久久久影院| 大型av网站在线播放|