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

    地球深部水循環(huán)的動力學機制與相關(guān)問題研究進展

    2018-06-08 03:27:50,
    西北地質(zhì) 2018年2期
    關(guān)鍵詞:島弧板片水循環(huán)

    ,

    (蘭州大學地質(zhì)科學與礦產(chǎn)資源學院 ,甘肅省西部礦產(chǎn)資源重點實驗室(蘭州大學)甘肅 蘭州 730000)

    地幔轉(zhuǎn)換帶(mantle transition zone,MTZ)是介于上下地幔之間的區(qū)域,是指410km和660km不連續(xù)面之間的地幔部分。地幔轉(zhuǎn)換帶內(nèi)部發(fā)生的地質(zhì)過程深刻地影響著與之相關(guān)的地質(zhì)事件,而地球深部的水循環(huán)在其中起到了重要作用,是解決許多重要地球科學問題的鑰匙(BIRCH,1952)。但是,由于缺乏直接來自于地幔轉(zhuǎn)換帶的樣品,對地幔轉(zhuǎn)換帶的研究主要使用的是一些交叉學科的研究手段,尤其以地球物理為主,由此產(chǎn)生的相關(guān)地球物理的數(shù)據(jù)又由于自生多解性的緣故,在詳細系統(tǒng)的說明有關(guān)地球深部的物質(zhì)結(jié)構(gòu)和地質(zhì)過程的時候仍然存在較大的困難,由此必須藉由高溫高壓的實驗來證實和解釋(金振民,1997)。在各個類型的構(gòu)造背景之中包含有種類豐富的地殼與地幔之間的耦合作用,而其中的俯沖帶地區(qū)是二者相互作用最為強烈的區(qū)域之一(趙子福等,2015)。在俯沖帶,巖漿的來源主要有地幔楔物質(zhì)的熔融以及俯沖洋殼和深海沉積物(THORKELSON,1996)。此外,巖漿產(chǎn)生的一個重要來源包括向地幔進行俯沖的板塊在深度達到一定程度的時候而釋放的水的參與,這是一個不可忽視的機制,這會影響巖漿的形成。例如,地幔物質(zhì)的絕熱減壓的過程,地幔楔物質(zhì)熔融的程度大小以及地幔楔的水通量都在其形成過程中起到了重要的作用(HELFFRICH et al.,2001)。因此,可以通過觀測一些數(shù)據(jù)來解釋目前難以觀察到的地方,而廣泛分布在地表的大陸溢流玄武巖(continental flood basalt,CFB)便是其中很好的一個例子。大陸溢流玄武巖對于殼幔循環(huán)的動力學機制研究意義重大。通過對全球范圍內(nèi)的大陸溢流玄武巖的地球化學特征的分析,結(jié)合其與構(gòu)造環(huán)境的相關(guān)關(guān)系,便可以推知其源區(qū)的一些性質(zhì),借此便能夠初探地球深部水循環(huán)的部分性質(zhì)。目前提出的地幔轉(zhuǎn)換帶水過濾器模型(WANG et al.,2016)對于地幔不均一性的問題給出了一個很好的解釋,地幔轉(zhuǎn)換帶的富水程度遠高于上、下地幔,由俯沖過程帶入的殘余板片下沉進入地幔轉(zhuǎn)換帶,殘余板片在地幔轉(zhuǎn)換帶發(fā)生相變并伴隨脫水作用,釋放出來的水導(dǎo)致部分熔融的發(fā)生,之后引發(fā)相應(yīng)的含水流體的上涌,交代上覆巖石圈地幔,使得巖漿作用產(chǎn)生的大陸溢流玄武巖具有類似島弧玄武巖的地球化學特征。而傳統(tǒng)的玄武巖構(gòu)造判別圖在區(qū)分這2類玄武巖時判別精度不夠高(王金榮等,2017;楊婧等,2016)。可以通過進一步的敏感的元素進行劃分(XIA et al.,2014),所以需謹慎地使用構(gòu)造判別圖。筆者在查閱大量文獻的基礎(chǔ)上,結(jié)合最新的研究進展,對地球深部水循環(huán)的一些理論進行對比總結(jié),旨在更好地理解巖漿成因機理與構(gòu)造動力學過程之關(guān)系。

    1 深部地球水循環(huán)研究歷史及趨向

    1.1 研究歷史回顧

    深部地球水循環(huán)的研究始于20世紀50至60年代(MARTIN et al.,1972)。伴隨著高溫高壓實驗技術(shù)的興起以及測試分析水平的提高,地幔轉(zhuǎn)換帶的深度環(huán)境條件已經(jīng)能夠成功的在實驗環(huán)境中模擬實現(xiàn),這使得我們對于在地幔轉(zhuǎn)換帶之中的礦物巖石學研究、礦物物理學研究、熱動力學狀態(tài)以及流變學諸多方面的研究有了長足的發(fā)展。然而要想掌握地球深部成分和結(jié)構(gòu)的情況,目前的手段是采取對比地球深部物理特征,結(jié)合特殊的礦物組合的性質(zhì),由此對地球深部的物理化學性質(zhì)進行評價和估計,我們需要認可并且采用一種符合實際情況的地幔成分模型,而對于地幔轉(zhuǎn)換帶物質(zhì)組成的分析必須建立在現(xiàn)有的地幔巖石學模型的基礎(chǔ)上進行操作。目前有2個上地幔物質(zhì)組成的相關(guān)模型,主要是pyrolite(地幔巖)模型(RINGWOOD,1962;RINGWOOD,1966;RINGWOOD,1975;RINGWOOD,1979;SUN,1982) 和piclogite(苦橄質(zhì)榴輝巖)模型(ANDERSON et al.,1984; BASS et al.,1984; DUFFY,et al.,1989)。其中,前者更為廣大學者所接受。

    由RINGWOOD首先明確指出的“pyrolite”模型,指的是地幔含有60%的橄欖石成分、30%的輝石成分以及數(shù)量較少的含鋁相礦物組成,其中含鋁相礦物根據(jù)壓力的大小情況可劃分為石榴子石、尖晶石或者斜長石(RINGWOOD,1991)。這種模型對于許多地球物理的相關(guān)現(xiàn)象都能夠很好的進行解釋,然而該模型中存在的最大問題就是其中的橄欖石成分在深度為410km時引起的密度和波速的突然改變都與地震波的實際測量值不同,并且要高于實際測量值(SHEARER et al.,1999;CAMMARANO et al.,2005),這是唯一存在于該模型的一個無法解釋的問題。

    地幔的另一種礦物組成模型是由ANDERSON et al(1984)提出的“piclogite”模型。在這個模型中,橄欖石在地幔之中的含量下降至40%~50%,而輝石和石榴子石含量卻一一上升,然而piclogite模型一直缺乏直接的高溫高壓實驗數(shù)據(jù)的支持。在后來的關(guān)于地震的波速及密度研究中的結(jié)果表明,實際測得的地球物理結(jié)果更加傾向于pyrolite模型。經(jīng)理論計算,2種模型的礦物組成隨著深度的變化而變化,pyrolite模型在實驗中受到廣泛應(yīng)用,加深了我們對于地幔的認識水平。

    想要清楚地了解和掌握地球深部的物質(zhì)狀況就必須借助科學的實驗,這就要求我們結(jié)合各類相關(guān)學科的知識,尤其著重于材料領(lǐng)域的科學、傳統(tǒng)的物理學化學手段以及新興的計算機科學來實現(xiàn)綜合性水平的提升?,F(xiàn)在高溫高壓實驗的技術(shù)已經(jīng)取得了長足的突破,尤其是激光加熱金剛石壓砧(LHDAC)和大體積壓機(LVP)這2種技術(shù)手段對于高溫高壓實驗有著至關(guān)重要的意義,在此基礎(chǔ)上再結(jié)合同步輻射技術(shù)(synchrotron radiation),就可以使得研究人員在高溫高壓的條件下進行原位(in situ)定量地對地球深部物質(zhì)的相變情況和相關(guān)礦物的物理性質(zhì)進行探究(OHTANI et al.,2008)。

    需要指出的是,對于地球深部的研究,關(guān)于地球深部賦存礦物的相變情況和結(jié)構(gòu)的研究是不可忽視的一個方面,因為地幔轉(zhuǎn)換帶中的地幔物質(zhì)的成分和結(jié)構(gòu)受控于相關(guān)礦物的相變狀態(tài),而這也對探究相關(guān)地幔流體的交代作用等一系列地球動力學過程具有指導(dǎo)作用。位于地幔410~660km深度范圍的地幔轉(zhuǎn)換帶(MTZ),處于該區(qū)域內(nèi)的橄欖石會經(jīng)歷連續(xù)的礦物相轉(zhuǎn)變過程位于410km鄰近區(qū)域的橄欖石就會轉(zhuǎn)變?yōu)橥咂澙?wadsleyite),這是對橄欖石而言的高壓礦物相,也是410~520km深度主要的礦物相,常稱為β相橄欖石;接下來在520km的區(qū)域內(nèi)發(fā)生相的轉(zhuǎn)變,轉(zhuǎn)變?yōu)榱治榈率?ringwoodite),學術(shù)界認為這是520~660km區(qū)域內(nèi)的主要礦物相,稱為γ相橄欖石;最后在660km附近分解為鈣鈦礦(perovskite)以及鎂方鐵石(magnesiowustite),為下地幔的主要礦物相。

    地球深部動力學研究離不開對其中水的研究,因為水可以降低地幔熔融溫度,通過流體傳輸元素,加強擴散作用,并且可能改變相轉(zhuǎn)變邊界的位置(WOOD,1995)。

    水對地球深部過程具有重要影響,如礦物的相轉(zhuǎn)變(SMYTH et al.,2002),地震波的性質(zhì)(INOUE et al.,1998) 等,同時也是地球深部發(fā)生流體交代作用的主體。在地球淺部,水主要賦存于俯沖板塊的各種含水礦物中,以及俯沖洋殼帶來的沉積物中,并且通過俯沖從而進入地球深部位置,在俯沖深度不斷增加的過程中,這之中裹攜的大量含水礦物狀態(tài)逐漸不穩(wěn)定從而引發(fā)脫水作用,這個過程中賦存的水的主要形式是名義上無水礦物(Nominal Anhydrous Minerals,NAMs),包括一部分致密的水鎂硅酸鹽之內(nèi)(BELL et al.,1992)。名義上無水礦物通常表述為在化學式中未見到H元素,但是實際中卻可以在晶格中觀察到一定數(shù)目的H元素的一類礦物。KUBO et al.(1998)觀察到水可以加強橄欖石α-β相的轉(zhuǎn)變,同時水可以改變瓦茲利石的流變學性質(zhì),使其攜帶0.02%的水含量。

    地幔轉(zhuǎn)換帶中的主要礦物,瓦茲利石和林伍德石都是名義上無水礦物,地幔轉(zhuǎn)換帶中的的含水數(shù)量都取決于這些名義上無水礦物。研究表明,地幔轉(zhuǎn)換帶中潛在的含水量較之地表含量很高(IONUE et al.,1995)。

    借助于地幔中主要儲水礦物的含量,可以大致估算出包括地幔轉(zhuǎn)換帶在內(nèi)的地幔各個圈層的儲水能力。表1表示地幔中各圈層的儲水能力,通過其中礦物的含水能力以及礦物的比例計算得出,礦物的比例依照RINGWOOD提出的地幔巖模型中的各個礦物的比例。

    表1 地幔各圈層的儲水能力(ANDERSON,2007)Tab.1 Water storage capacity of each mantle layer (ANDERSON,2007)

    1.2 存在的爭議

    目前關(guān)于地球深部水循環(huán)存在許多爭議(GREEN et al.,2008) 。其中之一就是地幔轉(zhuǎn)換帶有無水的參與,也就是轉(zhuǎn)換帶的“干”“濕”之爭,盡管實驗結(jié)果顯示名義上無水礦物都含有一定能力的儲水本領(lǐng),但并不代表轉(zhuǎn)換帶中就一定儲存著大量的水,而且根據(jù)不同的模型的計算,具體含水多少等問題目前還存在很大爭議(OHTANI, 2005) 。

    現(xiàn)已發(fā)現(xiàn)地幔中名義上無水礦物可以含有高達3%的水(KOHLSTEDT et al.,1996)。而且如果礦物含水的話,其物理性質(zhì)與不含水時將會有很大的差異(HIRSCHMANN,2006)。實驗也證明,地幔轉(zhuǎn)換帶中存在的林伍德石可以含有相當多的水。OHTANI et al.(2001) 發(fā)現(xiàn)林伍德石的含水量與溫度的變化呈現(xiàn)負相關(guān)性,伴著溫度的升高而呈現(xiàn)降低的趨勢,原因在于相邊界與林伍德石的含水量有密切的關(guān)系,水參與后的影響與溫度的影響正好相反,尤其是溫度較低時水影響的程度更甚,這是因為水的傳輸能力與板片的熱狀態(tài)有關(guān)。例如,熱的年輕的板片只能將很少水帶入地幔深部,然而冷的板片在臨界條件(600℃,6 GPa)則能夠攜帶大量水進入深部地幔(KUBO et al.,1998)。660 km地震波速的大幅下降也被認為是轉(zhuǎn)換帶中水存在的證據(jù)之一(LITASOV et al.,2005)。

    盡管實驗的結(jié)果說明了地幔轉(zhuǎn)換帶含水的事實是存在的,不過不同實驗得出的具體含水數(shù)目則存在較大差異。HUANG 等(2005)認為,地幔轉(zhuǎn)換帶中含水,但是含水量在不同的地區(qū)存在很大的差異,而且通過全球地幔轉(zhuǎn)換帶電導(dǎo)率的研究表明,其在全球的分布存在著不均一性,而且通過和大地電磁測量實驗的結(jié)果相互比較得出的結(jié)論是,在無水條件下橄欖石的電導(dǎo)率與地球物理測量結(jié)果存在很好的相關(guān)性,由此可以推導(dǎo)得出,水分布的不均一性是全球范圍內(nèi)地幔轉(zhuǎn)換帶的一個普遍現(xiàn)象,尤其在一些構(gòu)造背景下,如筆者討論的俯沖帶環(huán)境,地幔轉(zhuǎn)換帶是富含一定量的水的,但是在有些地區(qū)則呈現(xiàn)出無水的情形,這就使得地幔轉(zhuǎn)換帶是否含水成為一個學術(shù)界研究討論的熱點問題。

    盡管賦存于地幔轉(zhuǎn)換帶中的某些礦物存在含有一定數(shù)量水的可能性,但是直接證據(jù)的缺乏,即目前尚未直接觀察到來自地幔轉(zhuǎn)換帶的樣品來證明其中是否含水。地幔轉(zhuǎn)換帶能否將如此多的水穩(wěn)定地保存其中仍然值得商榷,換句話說,高的儲水能力并不代表著實際的高含水量。

    1.3 地球深部水循環(huán)研究新進展

    隨著板塊構(gòu)造理論逐漸被學術(shù)界廣泛地接受,在20世紀80年代,伴隨著玄武巖構(gòu)造環(huán)境判別方法的發(fā)展,使得玄武巖在地球動力學、大地構(gòu)造背景研究中所體現(xiàn)出的影響和作用愈來愈大(王金榮等,2017)。20世紀90年代所使用判別圖由于時間久遠、研究區(qū)域普遍較小、研究思路以及研究手段的匱乏以及分析技術(shù)的瓶頸所限,從發(fā)展的角度來看,不足之處難以避免,許多不同類型的玄武巖在判別圖上因為重疊區(qū)域太大而難以區(qū)分(LI et al.,2015)。WANG 等(2015,2016)在基于全球6個大火成巖省大陸溢流玄武巖以及與裂谷作用有關(guān)的美國盆地-山脈省的玄武巖數(shù)據(jù)進一步驗證了玄武巖構(gòu)造判別圖的置信度,認為利用主要主量元素、微量元素及其原始地幔標準化蛛網(wǎng)圖等綜合的判別體系是能夠區(qū)別類島弧地球化學特征的大陸溢流玄武巖與島弧玄武巖的,并且提出了一個新的地幔轉(zhuǎn)換帶水過濾器模型(WANG et al.,2015),進一步指出地球深部流體(主要是水)再循環(huán)是大陸溢流玄武巖富集流體活動性元素的重要機制,由此可見地球深部水循環(huán)所起到的重要作用。

    該模型提出,來自地幔轉(zhuǎn)換帶內(nèi)殘余板片熔融產(chǎn)生的含水流體的上涌組分是構(gòu)成大陸溢流玄武巖的主要來源。俯沖板塊產(chǎn)生的流體在形成弧巖漿和新地殼的過程中占有十分重要的地位,殘余板片物質(zhì)繼續(xù)向深部俯沖,因此對深部地幔的地球化學組成產(chǎn)生了一系列的影響,在板內(nèi)玄武巖巖漿中可能伴隨有部分古老板塊的參與(國坤等,2016)。而在地表占絕大多數(shù)的大陸溢流玄武巖往往顯示出類似島弧玄武巖的地球化學特征,其原因也可以通過地球深部水循環(huán)得到很好地解釋,并且從根本上將其與真正的島弧玄武巖區(qū)分開來(WANG et al.,2015;WANG et al.,2016),而傳統(tǒng)意義上的構(gòu)造判別圖在區(qū)分兩者時得到的結(jié)果往往不很理想,所以重新理解大陸溢流玄武巖的起源就顯得尤為重要。因為判別玄武巖的構(gòu)造環(huán)境對于接下來的工作開展尤為關(guān)鍵,如果能夠基于一個正確且合理的結(jié)論,與之相關(guān)的工作也會進行的更加科學和系統(tǒng)。

    玄武巖的地球化學特征取決于其源區(qū)的組成,熔融的條件(溫度、壓力以及有無含水相)以及之后發(fā)生的部分熔融和AFC作用等(DUNCAN,1987)。往往只有這些地質(zhì)過程對于某個特殊的構(gòu)造環(huán)境顯示出獨特的鑒定特征時,這些玄武巖的組成才能作為構(gòu)造環(huán)境的指示器(XIA, 2014)。因此,這些弧形特征是如何產(chǎn)生的就成為了區(qū)分弧型和在俯沖環(huán)境下產(chǎn)生的真正的島弧玄武巖的關(guān)鍵。

    1.4 與地球深部水循環(huán)有關(guān)的CFB的起源

    目前一致的看法是CFB的類似島弧玄武巖地球化學特征都可以用直接或間接地來自從俯沖板片中釋放出來的水(包括其他流體相)來解釋(IVANOV et al., 2013)。地球化學和巖石學都說明了類似島弧玄武巖的大陸玄武巖的地球化學特征和富集的放射性同位素特征都來自于交代過的、長期遠離地幔對流的地幔源區(qū)(MERLE et al.,2014)。但是,俯沖板塊產(chǎn)生的流體是如何影響和控制這些類似島弧玄武巖特征仍是一個未解決的問題。

    有的學者認為部分熔融過程發(fā)生于無水的大陸巖石圈地幔(SCLM)橄欖巖固相線條件下,顯然干的SCLM不可能是CFB的源區(qū),因此地幔轉(zhuǎn)換帶中的含水能力能夠解釋這個問題,這也與GALLAGHER(GALLAGHER et al.,1992)等人的觀點相一致。如果巖石圈地幔含有足夠多的水(0.3%H2O),橄欖巖固相線可以降低500℃,在該條件下,巖石圈地幔的部分熔融就可以作為CFB的主要源區(qū)。

    FUKAO et al.(2009)結(jié)合數(shù)學建模與地球物理的研究說明俯沖過程產(chǎn)生的大量水可以到達地幔轉(zhuǎn)換帶。因此,類似島弧玄武巖地球化學特征或許與地幔轉(zhuǎn)換帶的深部地球流體循環(huán)有關(guān)(IVANOV et al.,2008)。

    對于大陸溢流玄武巖來說,這種大規(guī)模噴發(fā)的玄武巖需要一個足夠大的且連續(xù)的含水富集源區(qū),并且需要一個相當長的穩(wěn)定時間來建立同位素異常儲庫。因此,首先可排除SCLM作為CFB源區(qū)的可能性。如果需要一個更深的地幔來源作為CFB的物質(zhì)來源和噴發(fā)基礎(chǔ),那么這種來源毫無疑問應(yīng)該是地幔轉(zhuǎn)換帶。最近觀察到的10 μm林伍德石的顆粒提供了地幔轉(zhuǎn)換帶含水本質(zhì)的第一手證據(jù)(PEARSON et al., 2014) 。

    在深部水循環(huán)模型中,當俯沖進入地幔轉(zhuǎn)換帶內(nèi)的殘余板塊由于重力原因繼續(xù)下沉時,發(fā)生了礦物相轉(zhuǎn)變,因此使得密度發(fā)生了變化,由于浮力的原因,便會產(chǎn)生含水流體的上涌,上涌的流體交代上覆地幔橄欖巖,由于不相容元素主要存在于角閃石等礦物中,當含水流體上升交代地幔橄欖巖時,礦物在流體中的溶解率受陽離子從晶體中擴散速率的影響,如金紅石這樣的富Ti氧化物,具有高電荷的陽離子,溶解速率低于硅酸鹽礦物,高度交代的地幔巖石中含有大量的金紅石,它們將導(dǎo)致LREE,LILE優(yōu)先于HFSE而進入流體中,這一過程有實驗結(jié)果所支持(MUNKER, 1998)。除了礦物相本身的影響,LREE、LILE、HFSE在流體中的溶解度還存在差異,流體活動性元素容易在流體中遷移,而HFSE不易溶解,因此交代地幔部分熔融產(chǎn)物具有高的LILE/HFSE值。LREE富集的原因是因為部分熔融程度低,也可能源于源區(qū)的富集,此外,當洋殼快速俯沖至地幔發(fā)生板片斷裂,下部的富集地幔通過裂隙上涌發(fā)生減壓熔融,也能夠形成不相容元素和LREE的富集。

    長石、輝石、橄欖石等主要造巖礦物中含有的H2O和Ce都是不相容元素,并且其不相容性非常相似,而且H2O/Ce值在部分熔融和結(jié)晶分異兩種過程中都不具有明顯的分異現(xiàn)象,所以巖漿的H2O/Ce值在排除了脫氣作用的影響后就代表了源區(qū)的H2O/Ce值(夏群科等,2015)。因此,巖漿的水含量和H2O/Ce值可以用來指示源區(qū)的組成(THY et al.,2006)。因為水的行為表現(xiàn)的像高度不相容元素,如Ce和K2O。在地-幔熔融過程中,H2O/Ce值可以作為水流體在熔融中的貢獻的指標(HAURI et al.,2006)。DIXON 等(2002)基于全球島弧玄武巖數(shù)據(jù),結(jié)合RUSCITTO 等(2012)的研究,認為流體活動性元素與非流體活動性元素(La,Nb)的比值與H2O/Ce高度相關(guān)。這就意味著流體活動性元素的富集可以用來限定初始的H2O/Ce值。地殼中的脫氣作用導(dǎo)致H2O的損失,但對REE沒有影響,因此會降低H2O/Ce值(PLANK et al.,2013)。具有較低的H2O/Ce的玄武巖源區(qū)可能有發(fā)生過強烈脫水作用的洋殼玄武巖的參與,而顯示出高的H2O/Ce值的玄武巖源區(qū)由2部分組成,既有俯沖洋殼的玄武巖的參與,同時還包括沉積物組分(夏群科等,2015)。在類似島弧玄武巖的陸內(nèi)玄武巖和島弧玄武巖中流體活動性元素的富集和HFSE的虧損的相似性都說明了深部地球水循環(huán)在形成地球大火成巖省中的重要作用。

    1.5 地球深部水循環(huán)模式

    一個新的模型建立的意義是因為它具有很好的解釋預(yù)測和指導(dǎo)作用,可以用來解釋許多地球動力學過程。傳統(tǒng)的地幔柱理論及其模式得到廣泛承認的原因正是因為它們很好地解釋了大洋板內(nèi)玄武巖的相關(guān)問題。在21世紀初期,在地幔柱模式的指導(dǎo)下,全世界范圍內(nèi)取得了一系列卓越的研究成果(CIVIERO et al.,2015)。然而目前這些傳統(tǒng)理論及模式并未像人們所希望的那樣在陸內(nèi)造山帶地區(qū)產(chǎn)生重大影響,這嚴重地撼動了地幔柱模式在地球動力學研究中的地位(FROMM et al.,2015)。于是以WANG 等(2015)為代表的一批學者,通過對世界范圍內(nèi)多個地區(qū)的CFB進行的觀察和研究,提出了以地球深部水循環(huán)為指導(dǎo)的地球深部水過濾器模型。該模型將地球深部水循環(huán)、殘余板片、含水流體地幔上涌、大規(guī)模板內(nèi)巖漿作用和超大陸的拼合與裂解結(jié)合在了一起(圖1)。由圖1可知,當俯沖板片進入地球深部的地幔轉(zhuǎn)換帶,許多殘余板片滯留在地幔轉(zhuǎn)換帶之中,在此階段地幔對流呈現(xiàn)分層性,滯留在地幔轉(zhuǎn)換帶的殘余板片最終俯沖進入下地幔,有一部分發(fā)生礦物相轉(zhuǎn)變回到上地幔,在下地幔頂部發(fā)生部分熔融,導(dǎo)致殘余板片釋放出流體和熔體成分,這會導(dǎo)致大規(guī)模的含水流體上涌,并且在地幔轉(zhuǎn)換帶引發(fā)交代作用。從俯沖板塊回到上覆地殼和軟流圈中的元素是地球上各個位置的元素物質(zhì)再循環(huán)的重要過程之一(國坤等,2016)。聚集在下地幔頂部的地幔柱物質(zhì)也會加強地球表面以及地幔轉(zhuǎn)換帶之間水的循環(huán),并且會促進含水流體的上涌。

    圖1 地幔轉(zhuǎn)換帶水過濾器模型圖(據(jù)WANG,2015)Fig.1 Mantle transition zone water-filting model(After WANG, 2015)

    與地幔柱模式不同,地幔轉(zhuǎn)換帶水過濾器模型將地球深部水流體循環(huán)、大規(guī)模陸內(nèi)巖漿作用和超大陸拼合與裂解組成了一個體系。這個模型說明了超大陸的聚合潛在地導(dǎo)致了含水大洋板片在地幔轉(zhuǎn)換帶內(nèi)的積累,并且可能改變了地幔對流的模式(FUKAO et al.,2015)。

    Pangea古陸和Rodinia超大陸的重建模型支持了這一模型,它們的周圍圍繞著環(huán)超級大陸俯沖帶(LI et al., 2009)。這些特征也適用于其他超級大陸,并且有重要的地球動力學啟示(LI et al.,2008)。地幔轉(zhuǎn)換帶內(nèi)的板片的黏度和礦物相轉(zhuǎn)變是決定滯留板片命運的關(guān)鍵因素(NAKAKUKI et al.,2010)。緩慢的輝石-石榴子石的相轉(zhuǎn)變說明滯留板片在地幔轉(zhuǎn)換帶內(nèi)存留超過了100 Ma(LI et al.,2008),大量的基于板片黏度的熱動力學模型同樣說明滯留板片的存在周期為100~130 Ma(LI et al, 2009)。這也與Pangea古陸和Rodinia超大陸的最終拼合與初始裂解之間的140~150 Ma的時間相一致(LI et al.,2008)。這種相似性可能對板片崩塌和超大陸的裂解的聯(lián)系有所啟示。地幔轉(zhuǎn)換帶內(nèi)相對穩(wěn)定的板片最終會俯沖進入下地幔。這會產(chǎn)生大規(guī)模的含水流體上涌并交代上覆大陸巖石圈,由于重力不穩(wěn)定以及礦物相轉(zhuǎn)變的原因,導(dǎo)致大陸巖石圈最下部的不穩(wěn)定(WINDLEY et al.,2010)。更重要的是,含水的玄武質(zhì)流體會匯聚在巖石圈和軟流圈的邊界(LAB),為加速板塊移動提供一個潤滑的作用(SAKAKMAKI et al.,2013)。殘余板片崩塌導(dǎo)致的上涌流體以及含水玄武質(zhì)熔體的活動性將會促進超級大陸的裂解。超級大陸裂解相應(yīng)的也可能會促進殘余板片的崩塌。由于含水玄武質(zhì)熔體的物理性質(zhì)偏向于集中在LAB,所以對于理解與克拉通破壞有關(guān)的玄武質(zhì)巖石成因十分關(guān)鍵。例如,中國華北克拉通地區(qū)。因此,地幔轉(zhuǎn)換帶內(nèi)部的殘余板片崩塌導(dǎo)致的上涌含水流體在超大陸的裂解與在地球表面和地幔轉(zhuǎn)換帶之間的水循環(huán)過程中起到了關(guān)鍵作用。如此大規(guī)模的含水流體上涌引發(fā)并控制了大規(guī)模的陸內(nèi)巖漿作用。

    這里再次強調(diào)起源于核幔邊界的地幔柱并不是產(chǎn)生類似島弧玄武巖CFB的重要條件,但是地幔柱的存在必然將會加強深部流體循環(huán)。WANG 等(2015)認為地幔柱和超級地幔柱產(chǎn)生的物質(zhì)加速了礦物相轉(zhuǎn)變,并且隨后加強了地球深部的流體循環(huán)。新的模型認為板片的崩塌導(dǎo)致了上涌的流體可能會引發(fā)最初的大陸裂解。這與Pangea古陸最初裂解產(chǎn)生的顯著的類似島弧玄武巖的大陸玄武巖相一致。

    2 相關(guān)的大陸溢流玄武巖實例

    2.1 CFB的鑒定特征

    隨著地球科學的發(fā)展,地質(zhì)學家對于大陸溢流玄武巖和島弧玄武巖的研究取得了一系列進展,但與此同時也產(chǎn)生了許多問題。由于來自地球的深部流體對于大陸地殼或巖石圈的混染作用可以給出似消減帶的信號,從而會使得人們將大陸玄武巖誤判為島弧玄武巖 (夏林圻等,2007)。因為大陸溢流玄武巖往往顯示出類似島弧玄武巖的地球化學特征(Arc-like signature),所以如何將大陸玄武巖與島弧玄武巖區(qū)分開來,就成為了大陸火山巖研究中一個十分重要的課題(ERNST et al.,2005;XIA,2014)。

    大陸溢流玄武巖(CFB)也稱為板內(nèi)玄武巖,大多集中于裂谷帶上,保存于被動大陸邊緣,部分與板內(nèi)深大斷裂有關(guān),廣泛的分布于大陸之上,是一類研究詳細的玄武巖類,巨大規(guī)模的大陸溢流玄武巖被稱為大火成巖省(王金榮等,2016)。前人的研究表明,大陸溢流玄武巖通常具有平坦的REE分配形式,并以缺乏Nb,Ta和Ti的負異常為特征,沒有受到混染的大陸溢流玄武巖通常是以(Th/Nb)n<1,Nb/La>>1,低87Sr/86Sr值(0.703 4~0.706 4),高εNd(t)值為特征(HAWKESWORTH et al.,1995)。而受到混染的大陸溢流玄武巖則具有較低的Nb/La(<1.0),較高的87Sr/86Sr值(0.704 3~0.713 7),較低的εNd(t)值,并且Nb,Ta和Ti負異常程度顯著高于俯沖帶的玄武巖。板片來源的流體/熔體會導(dǎo)致低的Nb/La值,巖石圈的混染也會使Nb/La值降低,但是會產(chǎn)生負的εNd(t)值,從而造成類似島弧玄武巖的地球化學特征(XIA,2014)。

    CFB通常分為高Ti和低Ti兩類,這反映了地幔源區(qū)的不均一性(SAUNDERS et al.,1992),這也與之前提到的地幔轉(zhuǎn)換帶內(nèi)部水的分布相對應(yīng)。在不相容元素中,Pb富集特征差異最大,Pb鑒于是一種俯沖遷移元素,俯沖組分中的Pb可能會掩蓋地幔中Pb的組分。Pb如果有明顯的正異常,那么CFB可能反映出板塊俯沖作用有陸源物質(zhì)的加入,而且富集地幔的不同還表現(xiàn)在Pb的異常方面,這也是地幔不均一性的表現(xiàn)(WEAVER, 1991)。Hf是一種非俯沖遷移元素,對揭示俯沖帶巖漿初始物源具有重要作用(JORDAN et al.,2007)。

    CFB擁有較寬的ε(Nd)和ε(Sr),表明具有不同的Sm/Nd、Rb/Sr 值長期存在的獨立源區(qū),這也從側(cè)面說明了對于CFB這種大規(guī)模噴發(fā)的玄武巖需要一個相當長的并且獨立穩(wěn)定的時間來建立同位素異常儲庫(PUFFER,2001)。

    2.2 研究實例

    根據(jù)PANG 等(2016)在內(nèi)蒙古的本巴圖組玄武巖的研究結(jié)果顯示,該地區(qū)玄武巖具有類似弧型島弧玄武巖的地球化學特征,具有較低的Sr87/Sr86初始值(0.704 2~0.704 8),εNd(t)(+8.99~+9.24)以及εHf(t)(+15.38~+15.65),在蛛網(wǎng)圖上富集Rb,U,Pb,Zr和Hf;虧損Nb,Ta,Ti和Sr,顯示出富集LILE,LREE,虧損HFSE的趨勢,并且具有平緩的REE配分形式圖,但是其中的Zr與Hf元素的富集以及Sr的虧損與典型的島弧玄武巖區(qū)分了開來。筆者認為造成該地區(qū)玄武巖具有類似島弧玄武巖地球化學特征的原因可能是地幔源區(qū)發(fā)生的交代作用導(dǎo)致的,而這一作用則發(fā)生在地幔轉(zhuǎn)換帶之中,通過俯沖作用進入的板片釋放出的含水流體引起的。

    在實際研究工作中,如果不結(jié)合實際的地質(zhì)情況和地球化學特征加以認真詳細地分析判斷,只是通過某些使用頻率較高卻又最容易與島弧玄武巖混淆的元素數(shù)值特征和地球化學圖解,或者利用構(gòu)造環(huán)境判別圖解,就會將表現(xiàn)出島弧型信號的大陸溢流玄武巖誤判為島弧玄武巖(夏林圻等,2007)。所以當利用玄武巖構(gòu)造環(huán)境判別圖解來鑒別有弧型特征的CFB時,必須謹慎地使用相關(guān)圖解,尤其是在面對被卷入造山帶中形成的環(huán)境未知的古老火山系時。在這種情況下,夏林圻等(2007)認為應(yīng)使用Zr-Zr/Y圖解,因為大陸玄武巖具有高的Zr(>70×10-6),Zr/Y(>3.4),這與是否被混染無關(guān),因此這些具有弧型特征的CFB仍然會落在板內(nèi)玄武巖(WPB)之內(nèi),這是因為Zr和Y在玄武巖巖漿上升時未受到影響的緣故,它們?nèi)詫⒈3制湓紳舛?TURNER et al.,2003)。但是,如果使用與Nb、Ta、Ti有關(guān)的判別圖時,由于地殼巖石中含有較低的TiO2含量和Nb、Ta含量,將會落入島弧玄武巖(IAB),因為類似島弧玄武巖的CFB與IAB的多元素原始地幔標準化分配型式近似程度高,在判斷時比較困難(KEPPLER,1996)。

    典型的島弧玄武巖(IAB)與大陸溢流玄武巖相比,總體上是以具有較低的Zr含量(<130×10-6)和Zr/Y值(<3.3)為特征;而大陸溢流玄武巖不論是否有弧型特征,它們都具有較高的Zr含量(>70×10-6)和Zr/Y值(>3.4)(PEARCE et al.,1979)。

    3 結(jié)論

    (1)地球深部研究是地學研究中的前沿學科,因此對于地球深部水循環(huán)的研究對地球內(nèi)部的許多動力學過程都具有重要意義。近幾十年來,地球深部水循環(huán)過程的研究成為地球科學研究的熱點。但是依然存在許多未解決的問題(地幔轉(zhuǎn)換帶中水的實際含量等)。

    (2)對地幔轉(zhuǎn)換帶結(jié)構(gòu)的研究是清楚了解地球深部水循環(huán)的重要前提,因此有必要在中國發(fā)展高溫高壓實驗技術(shù),深入探索地球深部結(jié)構(gòu)和性質(zhì),探討地球深部水循環(huán)與地球動力學之內(nèi)在的聯(lián)系。

    (3)大陸溢流玄武巖(CFB)的成因比較復(fù)雜,有可能與地球深部的水循環(huán)有關(guān),而CFB常因具有類似島弧玄武巖特征常被誤判為島弧玄武巖,所以正確理解CFB的起源對于正確區(qū)分二者具有重要的意義。

    (4)以地幔對流為基礎(chǔ)的地幔柱模式常用于解釋與大洋玄武巖有關(guān)的地質(zhì)過程,而以地球深部水循環(huán)為基礎(chǔ)的地幔轉(zhuǎn)換帶水過濾器模型則為復(fù)雜的超大陸旋回、深部流體循環(huán)、陸內(nèi)巖漿作用共存模式。地球深部水循環(huán)模型尚在探索階段,這2種模式不會區(qū)分絕對的對錯,只是在不同的情況下適用性不同。

    (5)傳統(tǒng)的玄武巖構(gòu)造判別圖在區(qū)分2類玄武巖(CFB、IAB)時判別精度不夠高,可以通過進一步的敏感的元素進行劃分,所以需謹慎地使用構(gòu)造判別圖。

    參考文獻(References):

    金振民. 我國高溫高壓實驗研究進展和展望[J]. 地球物理學報, 1997, 40(增刊Ⅰ):71-81.

    JIN Zhenmin. The progress and perspectives of high-T and high-p experimental study in China[J]. Chinese Journal of Geophysics, 1997, 40(SupplⅠ): 71-81.

    趙子福,戴立群,鄭永飛. 大陸俯沖帶兩類殼幔相互作用[J]. 中國科學(地球科學), 2015, 45: 900-915.

    ZHAO Zifu, DAI Liqun, ZHENG Yongfei. Two types of the crust-mantle interaction incontinental subduction zones[J]. Science China: Earth Sciences, 2015, 45: 900-915.

    王金榮,陳萬峰,張旗,等. N-MORB和E-MORB數(shù)據(jù)挖掘——玄武巖判別圖及洋中脊源區(qū)地幔性質(zhì)的討論[J]. 巖石學報, 2017, 33(3):993-1005.

    WANG Jinrong, CHEN Wanfeng, ZHANG Qi, et al. Preliminary research on data mining of N-MORB and E-MORB: Disdussion on method of basalt discrimination diagrams and the character of MORB’s mantle sourse[J]. Acta Petrologica Sinica,2017, 33(3):993-1005.

    夏林圻,夏祖春,徐學義,等.利用地球化學方法判別大陸玄武巖和島弧玄武巖[J]. 巖石礦物學雜志,2007,26(1):77-89.

    XIA Linqi,XIA Zuchun,XU Xueyi, et al. The discrimination between continental basalt and island arc basalt based on geochemical method[J]. Acta Petrologica et Mineralogica,2007,26(1):77-89.

    楊婧,王金榮,張旗,等. 弧后盆地玄武巖( BABB) 數(shù)據(jù)挖掘:與MORB 及IAB的對比[J]. 地球科學進展,2016,31(1): 66-77.

    YANG Jing,WANG Jinrong,ZHANG Qi,et al. Back-arc basin basalt( BABB) data mining: Comparison with MORB and IAB[J]. Advances in Earth Science,2016,31(1):66-77.

    周春銀,金振民,章軍鋒.地幔轉(zhuǎn)換帶:地球深部研究的重要方向[J]. 地學前緣,2010,17(3):90-113.

    ZHOU Chunyin,JIN Zhenmin,ZHANG Junfeng. Mantle transition zone: An important field in the studies of Earth’s deep interior[J]. Earth Science Frontiers, 2010,17(3):90-113.

    楊翠平,金振民,吳耀.地幔轉(zhuǎn)換帶中的水及其地球動力學意義[J]. 地學前緣,2010,17(3): 114-126.

    YANG Cuiping,JIN Zhenmin,WU Yao. Water in the mantle transition zone and its geodynamic implications[J]. Earth Science Frontiers,2010,17(3): 114-126.

    張鴻翔,徐志方,馬英軍,等.大陸溢流玄武巖的地球化學特征及起源[J]. 地球科學,2001,26(3): 261-268.

    ZHANG Hongxiang,XU Zhifang,MA Yingjun,et al. Geochemical features and origin of continental flood basalts[J]. Earth Science,2001,26(3): 261-268.

    國坤,翟世奎,于增慧,等.板塊俯沖對巖漿作用影響的同位素地球化學示蹤研究[J].海洋科學,2016,40(6):126-132.

    GUO Kun,ZHAI Shihui,YU Zenghui,et al.Advances in isotopic geochemistry tracing for the influence of subduction over magmatism[J].Marine Sciences,2016,40(6):126-132.

    夏群科,劉佳,陳歡,等.大陸玄武巖原始水含量的測定及其對源區(qū)組分的制約[J].巖石礦物學雜志,2015,34(3):371-381.

    XIA Qunke,LIU Jia,CHEN Huan,et al.Estimation of water content of primary magma for continental basalts and itsconstraint on source components[J].Acta Petrologica et Mineralogica,2015,34(3):371-381.

    王金榮,陳萬峰,張旗,等.MORB數(shù)據(jù)挖掘:玄武巖判別圖反思[J].大地構(gòu)造與成礦學,2017,41(2):420-431.

    WANG Jinrong,CHEN Wanfeng,ZHAN Qi,et al.MORB Data Mining:Reflection of Basalt Discrimination Diagram[J].Geotectonica et Metallogenia,2017,41(2):420-431.

    王金榮, 潘振杰, 張旗,等. 大陸板內(nèi)玄武巖數(shù)據(jù)挖掘:成分多樣性及在判別圖中的表現(xiàn)[J]. 巖石學報, 2016, 32(7):1919-1933.

    WANG Jinrong,PAN Zhenjie,ZHANG Qi,et al.Intra-continental basalt data mining: The diversity of their constituents and the performance in basalt discrimination diagrams[J].Acta Petrologica Sinica,2016,32(7):1919-1933.

    BIRCH F. Elasticity and constitution of the Earth’s interior[J]. Journal of Geophysical Research, 1952, 57(2): 277-286.

    THORKELSON D J. Subduction of diverging plates and the principles of slab window formation[J]. Tectonophysics, 1996, 255:47-63.

    HELFFRICH G R, WOOD B J. The Earth’s mantle[J]. Nature, 2001, 412:501-507.

    WANG, SIMON A W, XU B, et al. Origin of arc-like continental basalts: Implications for deep-Earth fluid cycling and tectonic discrimination[J]. Lithos,2015,http://dx.doi.org/10.1016/j.lithos.2015.12.014.

    MARTIN R F, DONNAY G. Hydroxyl in the mantle[J]. American Mineralogist, 1972, 57:554-570.

    ANDERSON D L. Theory of Earth[M]. Boston: Blackwell Scientific, 1989: 34-44.

    RINGWOOD A E. A model for upper mantle[J]. Journal of Geophysical Research, 1962, 67(2):857-867.

    RINGWOOD A E. Mineralogy of the mantle[M]. Hurley P M. Advances in Earth Science. Cambridge: Massachusetts Institute of Technology Press, 1966: 357-398.

    RINGWOOD A E. Compsition and Petrology of the Earth’s Mantle[M]. New York: McGraw-Hill,1975.

    RINGWOOD A E. Origin of the Earth and Moon[M]. New York: Springer-Verlag, 1979.

    SUN S S. Chemical composition and origin of the Earth’s primitive mantle[J]. Geochemical et Cosmochimica Acta, 1982, 46(2): 179-192.

    ANDERSON D L, BASS J D. Mineralogy and composition of the upper mantle[J]. Geophysical Research Letters, 1984, 11(7): 637-640.

    BASSJ D, ANDERSON D L. Composition of the upper mantle:Geophysical tests of two petrological models[J]. Geophysical Research Letters, 1984, 11(3): 229-232.

    DUFFY T S, ANDERSON D L. Seismic velocities in mantle minerals and the mineralogy of the upper mantle[J]. Journal of Geophysical Research, 1989, 94(B2): 1895-1912.

    RINGWOOD A E. Phase transformations and their bearing in the constitution and dynamics of the mantle[J]. Geochimica et Cosmochimca Acta, 1991, 55(8): 31-51.

    SHEARER P M, FLANAGAN M P. Seismic velocity and density jumps across the 410-and 660-kilometer discontinuities[J]. Sciences,1999, 285: 1545-1548.

    CAMMARANO F, GOES S, et al. Is a pyrolitic adiabatic mantle compatible with seismic data[J]. Earth and Planetary Science Letters, 2005, 232:227-243.

    OHTANI E, SAKAI T. Recent advances in the study of mantle phase transitions[J]. Physics of the Earth and Planetary Interiors,2008,170(3/4): 240-247.

    SMYTH J R, FROST D J. The effect of water on the 410 km discontinuity: An experimental study[J]. Geophysical Research Letters, 2002,29. Doi: 101029/2001G L014418.

    INOUE T, WEIDNER D J, et al. Elastic properties of hydrous ringwoodite in Mg2SiO4[J]. Earth and Planetary Sciences Letters,1998, 160(1/2): 107-113.

    BELL D R, ROSSMAN G R. Water in Earth’s mantle: The role of nominally anhydrous minerals[J]. Sciences,1992,255:1391-1397.

    INOUE T, YURIMOTO H, et al. Hydrous modified spinel, Mg1.75SiH0.5O4: A new water reservoir in the mantle transition region[J]. Geophysical Research Letters, 1995,22(2): 117-120.

    ANDERSON D L. New Theory of the Earth[M]. Cambridge: Cambridge University Press, 2007.

    GREEN H W, CHEN W, et al. Water is not Recycled into the Deep Mantle in Subducting Lithosphere[A]. American Geophysical Union[C], Fall Meeting, 2008: T13C1966G.

    OHTANI E. Water in the mantle[J]. Elements,2005,1(1):25-30.

    KOHLSTEDT D L, KEPPLER H, et al. Solubility of water in the α,β and γ phases of (Mg,Fe)2SiO4[J]. Contributions to Mineralogy and Petrology, 1996, 123(4): 345-357.

    HIRSCHMANN M M. Water,melting and the deep Earth H2O cycle[J]. Annual Review of Earth and Planetary Sciences, 2006, 34(1): 629-653.

    OHTANI E, TOMA M, et al. Stability of dense hydrous magnesium silicate phases and water storage capacity in the transition zone and lower mantle[J]. Physics of the Earth and Planetary Interiors, 2001, 124(1/2): 105-117.

    LITASOV K D, OHTANI E, et al. Wet subduction versus cold subduction[J]. Geophysical Research Letters, 2005, 32:L13312,doi: 10.1029/2005GL022921.

    HUANG X, XU Y, et al. Water content in the transition zone from electrical conductivity of wadsleyite and ringwoodite[J]. Natures,2005, 434:746-749.

    LI CS, TANG QY, et al. Trace element in discrimination diagrams. Lithos, 2015,232:76-83.

    WANG XC, SIMON A W, et al. Continental flood basalts derived from the hydrous mantle transition zone[J]. Nat.Commun.6:7700 doi:10.1038/ncomms8700(2015).

    ARNDT N T, CZAMANSKE G K, et al. Mantle and crustal contributions to continental flood vocanism[J]. Tectonphysics,1993, 223:39-52.

    DUNCAN A R. The Karoo igneous province - a problem area for inferring tectonic setting from basalt geochemistry[J]. Journal of Volcanology and Geothermal Research, 1987, 32: 13-34.

    XIA LQ. The geochemical criteria to distinguish continental basalts from arc related ones[J]. Earth-Science Reviews, 2014, 139: 195-212.

    IVANOV A V, LITASOV K D. The deep water cycle and blood basalt volcanism[J]. International Geology Review,2013,56: 1-14.

    MEREL R, MARZOLI A, et al. Sr, Nd, Pb and Os isotope systematic of CAMP tholeiites from Eastern North America(ENA): evidence of a subduction-enriched,mantle source[J]. Journal of Petrology,2014,55: 133-180.

    GALLAGHER K, HAWKESWORTH C. Dehydration melting and the generation of continental flood basalts[J]. Natures, 1992, 358: 57-59.

    FUKAO Y, OBAYASHI M, et al.Stagnant slabs: a review[J]. Annual Review of Earth and Planetary Sciences, 2009, 37: 19-46.

    IVANOV A, DEMONTEROVA E, et al. Low-Ti melts from the southeastern Siberian Traps Large Igneous Province: evidence for a water-rich mantle source[J]. Journal of Earth System Science,2008, 117:1-21.

    PEARSON D G, BRENKER F E, et al. Hydrous mantle transition zone indicated by ringwoodite included within diamond[J]. Nature, 2014,507: 221-224.

    MUNKER C. Nb/Ta fractionation in a Cambrian arc/back arc system, New Zealand: source constrains and application of refined ICPMS techniques[J]. Chem Geol, 1998, 144(1-2):23-45.

    THY P, LESHER C E, et al. Experimental constrains on the skaergaard liquid line of descent[J]. Lithos, 2006, 92: 154-180.

    HAURI E H, GAETANI G A, et al. Partitioning of water during melting of the Earth’s upper mantle at H2O-undersaturated conditions[J]. Earth and Planetary Science Letters.,2006, 248: 715-734.

    DIXON J E, LEIST L, et al. Recycled dehydrated lithosphere observed in plume-influenced mid-ocean-ridge basalt[J]. Nature,2002, 420:385-389.

    RUSCITTO D M, WALLACE P J, et al. Global variations in H2O/Ce: 2 relationships to arc magma geochemistry and volatile fluxes[J]. Geochemistry, Geophysics, Geosystems,2012, 13, Q03025.

    PLANK T, KELLEY K A, et al. Why do mafic arc magmas contain ~ 4wt% water on average[J]. Earth and Planetary Science Letters, 2013,364: 168-179.

    CIVIERO C, HAMMOND J O S, et al. Multiple mantle upwellings in the transition zone beneath the northern EastAfrican Rift system from relative P-wave travel-time tomography[J].Geochemistry,Geophysics,Geosystems,2015,16:2949-2968.

    FROMMT, PLANERT L, et al. South Atlantic opening: a plume-induced breakup[J]. Geology,http//dx.doi.org/10.1130/g36936.1.

    LI Z X, ZHONG S. Supercontinent-superplume coupling, true polar wander and plume mobility: Plate dominance in whole-mantle tectonics[J]. Phys, Earth Planet, Inter, 2009, 176: 143-156.

    LI Z X et al. Assembly, configuration, and break-up history of Rodinia: a synthesis[J]. Precambrian Research, 2008, 160: 179-210.

    NAKAKUKI T, TAGAWA M, et al. Dynamical mechanisms controlling formation and avalanche of a stagnant slab[J]. Earth Planet, 2010, 183: 309-320.

    WINDLEY B F, MARUYAMA S, et al. Delamination/thinning of sub-continental lithospheric mantle under Eastern China: the role of water and multiple subduction[J]. American Journal of Science,2010, 310: 1250-1293.

    SAKAKMAKI T et al. Ponded melt at the boundary between the lithosphere and asthenosphere[J]. Nature Geoscience, 2013, 6: 1041-1044.

    ERNST R E, BUCHAN K L, et al. Frontiers in large igneous province research[J]. Lithos, 2005,79: 271-297.

    HAWKESWORTH C, TURNER S, et al. Calc-alkaline magmatism, lithospheric thinning and extension in the Basin and Range[J]. Geophysical Research, 1995, 100(10):271-286.

    SAUNDERS A D, STOREY M, et al. Consequences of plume-lithossphere interaction[A]. London, 1992,68: 41-60.

    WEAVER B L. The origin of ocean island basalt end-member compositions: trace element and isotopic constraints[J]. Earth Planet Letters, 1991, 104: 381-397.

    JORDAN F, BERTRAND H, et al. Major and Trace element and Sr,Nd,Hf, and Pb isotope compositions of the Karoo Large Igneous Province, Botswana-Zimbabwe: lithosphere vs mantle plume contuibution[J]. Journal of Petrology, 2007, 48: 1043-1077.

    PUFFER J H. Contrasting high field strength elements contents of continental flood basalts from plume versus reactivated-arc sources[J]. Geology, 2001, 29: 675-678.

    PANG C J, WANG X C, et al. Late Carboniferous N-MORB-type basalts in central Inner Mongolia, China: Products of hydrous melting in an intraplate setting[J]. Lithos, 2016,http://dx.doi.org/10.1016/j.lithos.2016.05.005.

    TURNER S, FODEN J, et al. Rates and peocesses of potassic magma evolution beneath Sangeang Apivolcano, east Sunda arc, Indonesia[J]. Petrol, 2003, 44(3): 491-515.

    KEPPLER H. Constraints from partitioning experiments on the composition of subduction-zone fluids[J]. Nature, 1996, 380: 237-240.

    PEARCE J, NORRY M. Petrogenetic implications of Ti,Zr,Y,and Nb variations in volcanic rocks[J]. Mineralogy and Petrology, 1979, 69: 33-47.

    WOOD B J. The effect of H2O on the 410-kilometer seismic discontinuity[J]. Science,1995,268:74-76.

    KUBO T, OHTANI E, et al. Effects of water on the α-β transformation kinetics in San Carlos olivine[J]. Science,1998,281:85-87.

    猜你喜歡
    島弧板片水循環(huán)
    用袋裝水模擬水循環(huán)
    核島板式換熱器缺陷原因分析及改進建議
    硫酸法鈦白粉水洗及水循環(huán)膜處理工藝
    化工管理(2021年7期)2021-05-13 00:46:24
    新型平行板反應(yīng)器板片動力響應(yīng)研究
    化工機械(2021年1期)2021-04-19 02:56:34
    伊豆-博寧-馬里亞納島弧地殼厚度分布及其對巖漿活動的指示
    東天山帕爾塔格西銅礦床地質(zhì)特征及找礦方向
    青海多彩銅多金屬礦區(qū)火山巖地球化學特征及其對島弧環(huán)境的限定
    板片斷裂失效分析
    科技視界(2020年17期)2020-07-30 14:03:27
    水循環(huán)高效礦井乏風熱泵系統(tǒng)分析與應(yīng)用
    同煤科技(2015年4期)2015-08-21 12:51:02
    40MN板片液壓機液壓系統(tǒng)設(shè)計
    国产精品.久久久| 亚洲性久久影院| 毛片女人毛片| 在线观看一区二区三区| 日本色播在线视频| 亚洲av中文av极速乱| 欧美日本亚洲视频在线播放| 一区二区三区高清视频在线| 干丝袜人妻中文字幕| 最近的中文字幕免费完整| 亚洲精品国产成人久久av| av国产免费在线观看| 中文欧美无线码| 日本免费一区二区三区高清不卡| 亚洲精品亚洲一区二区| 亚州av有码| 卡戴珊不雅视频在线播放| 欧美色视频一区免费| 亚洲精品456在线播放app| 成人综合一区亚洲| 97热精品久久久久久| 综合色丁香网| 69人妻影院| 一区二区三区四区激情视频| 国产精品久久久久久av不卡| 亚洲婷婷狠狠爱综合网| 亚洲av成人av| 亚洲av男天堂| 成人三级黄色视频| 国产国拍精品亚洲av在线观看| 天美传媒精品一区二区| 毛片一级片免费看久久久久| 欧美成人一区二区免费高清观看| 亚洲国产精品成人久久小说| 久久精品久久久久久久性| 色尼玛亚洲综合影院| 一级毛片aaaaaa免费看小| 欧美+日韩+精品| 人妻制服诱惑在线中文字幕| 日本黄色视频三级网站网址| 亚洲精品乱码久久久久久按摩| 国产精品福利在线免费观看| 亚洲伊人久久精品综合 | 亚洲国产精品专区欧美| 久久久久久久午夜电影| 精品免费久久久久久久清纯| 久99久视频精品免费| 成人综合一区亚洲| 久久久久久九九精品二区国产| 精品午夜福利在线看| 毛片女人毛片| 国产色婷婷99| 午夜福利高清视频| 狂野欧美激情性xxxx在线观看| 免费黄色在线免费观看| 国产激情偷乱视频一区二区| 日韩av在线免费看完整版不卡| 少妇被粗大猛烈的视频| 亚洲自偷自拍三级| 国产黄a三级三级三级人| 亚洲成人中文字幕在线播放| 一卡2卡三卡四卡精品乱码亚洲| 老司机福利观看| 欧美高清性xxxxhd video| 亚洲av成人精品一区久久| 啦啦啦韩国在线观看视频| 一级爰片在线观看| 久久久精品大字幕| 少妇高潮的动态图| 欧美三级亚洲精品| 亚洲精品乱码久久久久久按摩| 黄片无遮挡物在线观看| 五月伊人婷婷丁香| 亚洲伊人久久精品综合 | 天美传媒精品一区二区| 亚洲最大成人手机在线| 欧美丝袜亚洲另类| 久久久久久久国产电影| av专区在线播放| 三级经典国产精品| 男人舔奶头视频| 菩萨蛮人人尽说江南好唐韦庄 | 三级毛片av免费| 老师上课跳d突然被开到最大视频| 特级一级黄色大片| 日韩av在线免费看完整版不卡| 精品午夜福利在线看| 欧美性感艳星| 国产69精品久久久久777片| 麻豆成人av视频| 亚洲国产精品sss在线观看| 我要搜黄色片| 别揉我奶头 嗯啊视频| 夫妻性生交免费视频一级片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | eeuss影院久久| 一二三四中文在线观看免费高清| 又爽又黄a免费视频| 亚洲中文字幕日韩| 中文精品一卡2卡3卡4更新| 午夜爱爱视频在线播放| 亚洲性久久影院| 3wmmmm亚洲av在线观看| 一区二区三区乱码不卡18| 真实男女啪啪啪动态图| 国产高潮美女av| 亚洲国产欧美人成| 国产单亲对白刺激| 国产白丝娇喘喷水9色精品| www.色视频.com| 国语对白做爰xxxⅹ性视频网站| 国产又黄又爽又无遮挡在线| 欧美精品一区二区大全| 久久久精品94久久精品| 国产视频内射| 黄色日韩在线| 免费av毛片视频| 久久久久久国产a免费观看| 日日摸夜夜添夜夜爱| 国模一区二区三区四区视频| 亚洲精品456在线播放app| 99久久人妻综合| 国产精品人妻久久久久久| 在线天堂最新版资源| 日韩精品有码人妻一区| 免费观看性生交大片5| 久久精品国产亚洲网站| 视频中文字幕在线观看| 淫秽高清视频在线观看| 国产黄片视频在线免费观看| 精品人妻熟女av久视频| 亚洲国产精品久久男人天堂| 久久久久久久久中文| 午夜福利在线观看免费完整高清在| 亚洲av男天堂| 能在线免费看毛片的网站| 国产精品久久久久久久久免| 麻豆成人av视频| 日日啪夜夜撸| 亚洲国产欧美在线一区| 91在线精品国自产拍蜜月| 少妇熟女欧美另类| 联通29元200g的流量卡| www.av在线官网国产| 日本免费一区二区三区高清不卡| 久久婷婷人人爽人人干人人爱| 深爱激情五月婷婷| 久久人人爽人人片av| 国产淫语在线视频| 亚洲中文字幕日韩| 最新中文字幕久久久久| 舔av片在线| 国产精品一区www在线观看| 欧美一级a爱片免费观看看| 久久鲁丝午夜福利片| 日产精品乱码卡一卡2卡三| 丰满人妻一区二区三区视频av| 国产亚洲午夜精品一区二区久久 | 国产成人aa在线观看| 女的被弄到高潮叫床怎么办| 天天一区二区日本电影三级| 3wmmmm亚洲av在线观看| 91精品伊人久久大香线蕉| 蜜桃久久精品国产亚洲av| 少妇高潮的动态图| 国产一区二区在线观看日韩| 最近中文字幕高清免费大全6| 蜜桃亚洲精品一区二区三区| 国产伦一二天堂av在线观看| 村上凉子中文字幕在线| 国产精品爽爽va在线观看网站| 国语对白做爰xxxⅹ性视频网站| 午夜福利成人在线免费观看| 亚洲丝袜综合中文字幕| 在线观看66精品国产| av国产免费在线观看| 听说在线观看完整版免费高清| 亚洲欧美中文字幕日韩二区| 国产色婷婷99| 天堂中文最新版在线下载 | 韩国av在线不卡| 麻豆一二三区av精品| 久久久久久九九精品二区国产| 亚洲av福利一区| 美女cb高潮喷水在线观看| 天堂影院成人在线观看| 插逼视频在线观看| www.色视频.com| 一本久久精品| 国产精品99久久久久久久久| 亚洲精品自拍成人| 久久久久久国产a免费观看| 大香蕉97超碰在线| 亚洲美女搞黄在线观看| 婷婷色av中文字幕| 尾随美女入室| av在线观看视频网站免费| 欧美性猛交黑人性爽| 亚洲成人久久爱视频| 欧美97在线视频| 亚洲电影在线观看av| 国产成人免费观看mmmm| 色播亚洲综合网| 特大巨黑吊av在线直播| 欧美人与善性xxx| 少妇人妻一区二区三区视频| 少妇裸体淫交视频免费看高清| 超碰97精品在线观看| 久久久久久久久久成人| 美女大奶头视频| 一区二区三区四区激情视频| 99久久人妻综合| 日韩av在线免费看完整版不卡| av卡一久久| 麻豆一二三区av精品| 久久这里只有精品中国| 麻豆精品久久久久久蜜桃| 老司机影院毛片| 91久久精品国产一区二区三区| 色尼玛亚洲综合影院| 你懂的网址亚洲精品在线观看 | 一级黄片播放器| 欧美97在线视频| 亚洲精品乱码久久久久久按摩| 亚洲av一区综合| 国产麻豆成人av免费视频| 美女被艹到高潮喷水动态| 永久免费av网站大全| 国产成年人精品一区二区| 乱系列少妇在线播放| 美女被艹到高潮喷水动态| 纵有疾风起免费观看全集完整版 | 纵有疾风起免费观看全集完整版 | 免费观看a级毛片全部| 在线观看美女被高潮喷水网站| 99热这里只有是精品在线观看| 国产精品野战在线观看| 国产三级中文精品| 久久久精品欧美日韩精品| 成年免费大片在线观看| 日日摸夜夜添夜夜添av毛片| 国产亚洲5aaaaa淫片| 九九热线精品视视频播放| 99热6这里只有精品| 级片在线观看| 一区二区三区乱码不卡18| 成年免费大片在线观看| 中文资源天堂在线| 亚洲av中文av极速乱| 免费大片18禁| 一级av片app| 不卡视频在线观看欧美| av播播在线观看一区| 97在线视频观看| 国产淫片久久久久久久久| 亚洲不卡免费看| 91精品伊人久久大香线蕉| 日本黄色片子视频| 日韩国内少妇激情av| 亚洲成人中文字幕在线播放| 久久韩国三级中文字幕| 国产精品永久免费网站| 看免费成人av毛片| 欧美精品一区二区大全| 国产视频首页在线观看| a级毛色黄片| 不卡视频在线观看欧美| 免费看a级黄色片| 成人漫画全彩无遮挡| 久久精品国产99精品国产亚洲性色| 国产精品乱码一区二三区的特点| 久久久久精品久久久久真实原创| 亚洲欧美清纯卡通| 国产精品美女特级片免费视频播放器| 亚洲欧美成人精品一区二区| 日韩欧美精品v在线| 久久韩国三级中文字幕| 成人美女网站在线观看视频| 亚洲第一区二区三区不卡| 亚洲精品成人久久久久久| 99久久无色码亚洲精品果冻| 一卡2卡三卡四卡精品乱码亚洲| 精品国产一区二区三区久久久樱花 | 久久久久久国产a免费观看| 久久久国产成人精品二区| 亚洲av二区三区四区| 国产欧美另类精品又又久久亚洲欧美| 久99久视频精品免费| 日本熟妇午夜| 高清日韩中文字幕在线| 亚洲欧美日韩高清专用| 午夜视频国产福利| 国产真实伦视频高清在线观看| 岛国在线免费视频观看| 老师上课跳d突然被开到最大视频| 人妻夜夜爽99麻豆av| 亚洲欧洲日产国产| 麻豆av噜噜一区二区三区| 亚洲最大成人中文| 国产成人免费观看mmmm| 久久久色成人| 波多野结衣巨乳人妻| 欧美又色又爽又黄视频| 久久久久久大精品| 超碰97精品在线观看| 国产午夜精品一二区理论片| 亚洲精品aⅴ在线观看| 我要搜黄色片| 级片在线观看| 少妇熟女欧美另类| 联通29元200g的流量卡| 午夜日本视频在线| 最近视频中文字幕2019在线8| 亚洲欧美日韩卡通动漫| 级片在线观看| 国产乱人视频| 搡老妇女老女人老熟妇| 亚洲天堂国产精品一区在线| 又粗又爽又猛毛片免费看| 国产黄色小视频在线观看| 国产乱人偷精品视频| 日韩人妻高清精品专区| 岛国在线免费视频观看| 天堂网av新在线| 村上凉子中文字幕在线| 一级毛片我不卡| 国产精品人妻久久久久久| 日韩制服骚丝袜av| 日韩av在线大香蕉| av专区在线播放| 我的老师免费观看完整版| 亚洲性久久影院| 国产三级在线视频| 久久久久九九精品影院| 啦啦啦观看免费观看视频高清| 2021天堂中文幕一二区在线观| 91在线精品国自产拍蜜月| 久久久久久久久久久丰满| 亚洲真实伦在线观看| 毛片女人毛片| 国产精品乱码一区二三区的特点| 日韩大片免费观看网站 | 色视频www国产| 亚洲中文字幕日韩| 日韩强制内射视频| 国产 一区 欧美 日韩| 亚洲精品日韩在线中文字幕| 毛片一级片免费看久久久久| 欧美一级a爱片免费观看看| 少妇人妻一区二区三区视频| 国产欧美另类精品又又久久亚洲欧美| 久久久国产成人精品二区| 亚洲欧美精品自产自拍| 国产男人的电影天堂91| 久久久久久久久久久丰满| 自拍偷自拍亚洲精品老妇| 18禁裸乳无遮挡免费网站照片| 国产伦在线观看视频一区| 国产一区二区亚洲精品在线观看| 精品人妻一区二区三区麻豆| 五月伊人婷婷丁香| 亚洲欧美中文字幕日韩二区| 欧美最新免费一区二区三区| 91久久精品国产一区二区三区| 九九爱精品视频在线观看| 永久网站在线| 免费观看性生交大片5| av专区在线播放| 99久久中文字幕三级久久日本| 国产v大片淫在线免费观看| 国产精品嫩草影院av在线观看| 人妻制服诱惑在线中文字幕| 床上黄色一级片| 在线播放国产精品三级| 亚洲精品日韩av片在线观看| 乱系列少妇在线播放| 高清午夜精品一区二区三区| 免费大片18禁| 成人欧美大片| 亚洲欧美成人综合另类久久久 | 波多野结衣巨乳人妻| 99热网站在线观看| 三级国产精品片| 三级男女做爰猛烈吃奶摸视频| 亚洲国产欧美人成| 欧美极品一区二区三区四区| av.在线天堂| 日本av手机在线免费观看| 又黄又爽又刺激的免费视频.| 久久99蜜桃精品久久| 亚洲最大成人手机在线| 日韩国内少妇激情av| 亚洲内射少妇av| 免费黄网站久久成人精品| 久久热精品热| 欧美成人一区二区免费高清观看| 国产亚洲精品av在线| 亚洲国产欧美在线一区| 国产精品久久电影中文字幕| 日韩一本色道免费dvd| 欧美成人一区二区免费高清观看| 特级一级黄色大片| 九色成人免费人妻av| 午夜福利成人在线免费观看| 国内精品一区二区在线观看| 国产成人午夜福利电影在线观看| 夫妻性生交免费视频一级片| 综合色av麻豆| 国产乱人偷精品视频| 日韩精品青青久久久久久| 国产一区二区在线av高清观看| 中文字幕亚洲精品专区| 中文字幕av成人在线电影| 午夜福利在线观看吧| 精品少妇黑人巨大在线播放 | 国产精品国产三级国产专区5o | 国产亚洲av嫩草精品影院| 少妇高潮的动态图| 青青草视频在线视频观看| 成人欧美大片| 日韩欧美精品v在线| 一区二区三区四区激情视频| 精品久久久噜噜| 插阴视频在线观看视频| 又爽又黄a免费视频| av免费观看日本| 中文乱码字字幕精品一区二区三区 | 亚洲自偷自拍三级| 午夜激情福利司机影院| 日本熟妇午夜| 亚洲最大成人av| 亚洲精品乱码久久久v下载方式| 中文字幕av成人在线电影| 一区二区三区免费毛片| 亚洲在久久综合| 男女国产视频网站| 精品久久久噜噜| 99热网站在线观看| 亚洲国产成人一精品久久久| 国产精品国产三级专区第一集| 亚洲欧美中文字幕日韩二区| 在线天堂最新版资源| 国产免费视频播放在线视频 | 欧美不卡视频在线免费观看| 嫩草影院精品99| 波多野结衣巨乳人妻| 日本一二三区视频观看| 久久久久九九精品影院| 一个人看视频在线观看www免费| 久久久精品94久久精品| 日韩视频在线欧美| 国产精品一及| 久久久a久久爽久久v久久| 国产精品,欧美在线| 国产精品久久电影中文字幕| 欧美日韩精品成人综合77777| 99久国产av精品| 中文天堂在线官网| av女优亚洲男人天堂| 国内精品美女久久久久久| 日韩亚洲欧美综合| 久久精品91蜜桃| 能在线免费看毛片的网站| 在线a可以看的网站| 天天躁日日操中文字幕| 搞女人的毛片| 欧美性猛交黑人性爽| 成人美女网站在线观看视频| 男女下面进入的视频免费午夜| 欧美3d第一页| 国产精品综合久久久久久久免费| 简卡轻食公司| 91午夜精品亚洲一区二区三区| 熟女人妻精品中文字幕| 国产爱豆传媒在线观看| 在线免费十八禁| 欧美激情久久久久久爽电影| 少妇裸体淫交视频免费看高清| 看片在线看免费视频| 亚洲国产欧美人成| 麻豆一二三区av精品| 久久精品综合一区二区三区| 中文字幕av在线有码专区| 国国产精品蜜臀av免费| 天堂av国产一区二区熟女人妻| 国产成人精品婷婷| 亚洲三级黄色毛片| 97在线视频观看| kizo精华| 久久婷婷人人爽人人干人人爱| av线在线观看网站| 男女下面进入的视频免费午夜| 午夜福利成人在线免费观看| 黄片wwwwww| 99久国产av精品| 26uuu在线亚洲综合色| 99久久人妻综合| 在线免费观看的www视频| 大又大粗又爽又黄少妇毛片口| 变态另类丝袜制服| 老女人水多毛片| 国产欧美另类精品又又久久亚洲欧美| 亚洲av成人精品一区久久| АⅤ资源中文在线天堂| 国产精品野战在线观看| 乱人视频在线观看| av在线亚洲专区| 一级二级三级毛片免费看| 国产午夜精品论理片| 嫩草影院新地址| 少妇熟女aⅴ在线视频| 国产精品精品国产色婷婷| 亚洲性久久影院| 岛国在线免费视频观看| 九草在线视频观看| 男女视频在线观看网站免费| 午夜免费男女啪啪视频观看| 精品久久国产蜜桃| 午夜免费男女啪啪视频观看| 日韩强制内射视频| 一级爰片在线观看| 99在线人妻在线中文字幕| 精品国内亚洲2022精品成人| 日韩欧美三级三区| 欧美一级a爱片免费观看看| 亚洲国产欧洲综合997久久,| 亚洲av中文字字幕乱码综合| 狂野欧美激情性xxxx在线观看| 国产伦理片在线播放av一区| 视频中文字幕在线观看| 黄色配什么色好看| 一区二区三区免费毛片| 日本黄色视频三级网站网址| 国产老妇伦熟女老妇高清| 国产精品综合久久久久久久免费| 国产中年淑女户外野战色| 国产乱人视频| 青春草国产在线视频| 国产亚洲5aaaaa淫片| 女人被狂操c到高潮| 日韩欧美精品v在线| 啦啦啦啦在线视频资源| 国产精品无大码| 极品教师在线视频| 97人妻精品一区二区三区麻豆| 久热久热在线精品观看| 日本午夜av视频| 成年av动漫网址| 国产成人免费观看mmmm| 成人美女网站在线观看视频| 亚洲内射少妇av| 久久久欧美国产精品| 免费大片18禁| 白带黄色成豆腐渣| 国产精华一区二区三区| 纵有疾风起免费观看全集完整版 | 在线观看av片永久免费下载| 久久久久久久久大av| 在线播放国产精品三级| 亚洲精品乱久久久久久| 久久精品熟女亚洲av麻豆精品 | 黄色配什么色好看| 91精品伊人久久大香线蕉| 亚洲成色77777| 成人特级av手机在线观看| 日本欧美国产在线视频| 午夜免费男女啪啪视频观看| 综合色av麻豆| 精品少妇黑人巨大在线播放 | 天堂影院成人在线观看| 国产精品,欧美在线| 午夜精品一区二区三区免费看| 国产亚洲午夜精品一区二区久久 | 亚洲va在线va天堂va国产| 日本黄大片高清| 亚洲成人久久爱视频| 日产精品乱码卡一卡2卡三| 嘟嘟电影网在线观看| 国产精品无大码| 亚洲国产成人一精品久久久| 黄色日韩在线| 日韩三级伦理在线观看| 国产 一区 欧美 日韩| 国产欧美另类精品又又久久亚洲欧美| 亚洲在线自拍视频| 九色成人免费人妻av| 久久久久久久久中文| av在线天堂中文字幕| 久久久久国产网址| 最后的刺客免费高清国语| 永久免费av网站大全| 桃色一区二区三区在线观看| 日本午夜av视频| 乱码一卡2卡4卡精品| 色综合站精品国产| 欧美日韩国产亚洲二区| 少妇猛男粗大的猛烈进出视频 | 国产亚洲av嫩草精品影院| 22中文网久久字幕| 亚洲一区高清亚洲精品| 久久婷婷人人爽人人干人人爱| 成人综合一区亚洲| 一级二级三级毛片免费看| 成年版毛片免费区| 久久国内精品自在自线图片| 国产三级中文精品| 亚洲av电影不卡..在线观看| av免费在线看不卡| 国产淫语在线视频| 国产高清国产精品国产三级 | a级毛色黄片| 亚洲av成人av| 三级经典国产精品| 女人被狂操c到高潮| 久久久久久久午夜电影|