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

    礦物微量元素組成用于火成巖構(gòu)造背景判別

    2016-12-12 08:38:58趙振華中國(guó)科學(xué)院廣州地球化學(xué)研究所廣東廣州510640
    大地構(gòu)造與成礦學(xué) 2016年5期
    關(guān)鍵詞:單斜火成巖輝石

    趙振華(中國(guó)科學(xué)院 廣州地球化學(xué)研究所, 廣東 廣州 510640)

    礦物微量元素組成用于火成巖構(gòu)造背景判別

    趙振華
    (中國(guó)科學(xué)院 廣州地球化學(xué)研究所, 廣東 廣州 510640)

    火成巖中的礦物, 特別是穩(wěn)定副礦物的微量元素組成可用于火成巖構(gòu)造–巖漿背景的判別。產(chǎn)于不同構(gòu)造背景火成巖中單顆粒鋯石原位微量元素含量的準(zhǔn)確測(cè)定, 積累了大量資料。本文介紹了用鋯石微量元素含量的二維核密度分布投影, 構(gòu)建了 U/Yb-Nb/Yb、Sc/Yb-Nb/Yb和 U/Yb-Sc/Yb構(gòu)造背景判別圖解, 可以區(qū)分洋中脊、地幔柱影響的洋島及俯沖有關(guān)的弧(大陸弧)等不同構(gòu)造背景形成的火成巖。單斜輝石、尖晶石、黑云母和角閃石的微量元素或主量元素組成同樣可用于構(gòu)造背景識(shí)別。

    構(gòu)造背景判別; 微量元素; 鋯石; 單斜輝石; 尖晶石; 黑云母; 角閃石

    盡管用巖石微量元素識(shí)別用于構(gòu)造背景的判別圖存在較多爭(zhēng)議, 但用微量元素判別基性巖漿巖構(gòu)造背景仍取得了較明顯的進(jìn)展(Vermeesch, 2006a, 2006b; 趙振華, 2007; Pearce 2014; Li et al., 2015)。主要包括以下三個(gè)方面:

    第一是大量、精確的主、微量元素?cái)?shù)據(jù)的積累。隨著電子探針?lè)治黾夹g(shù), 特別是激光等離子體質(zhì)譜(LA-ICP-MS)技術(shù)的快速發(fā)展, 基性巖全巖(包括單顆粒礦物, 特別是鋯石原位)微量元素含量的準(zhǔn)確測(cè)定積累了大量數(shù)據(jù), 建立了典型巖石類型的數(shù)據(jù)庫(kù),如 PETDB和 GEOROC數(shù)據(jù)庫(kù)(網(wǎng)址: http: //www. earthchem.org/petdb和 http: //www.rankinsider.com/ georoc.mpch-mainz.gwdg.de), 大大拓寬了全巖微量元素組成判別火成巖形成構(gòu)造背景方法的適應(yīng)范圍, 特別是為遭受蝕變巖石的構(gòu)造背景恢復(fù)提供了新的工具(Grimes et al., 2007, 2015; Gonzalez-Jimenez et al., 2012)。

    第二是更有效地運(yùn)用多變量統(tǒng)計(jì)分析方法。在發(fā)表的有關(guān)玄武巖構(gòu)造背景的傳統(tǒng)判別圖解中, 不同構(gòu)造背景分類之間界線的確定往往是用兩個(gè)或三個(gè)變量數(shù)據(jù), 不是數(shù)值化的, 在統(tǒng)計(jì)上不嚴(yán)格, 由此劃分的界線常常是不準(zhǔn)確的、或有較大重疊。例如, 為了將玄武巖構(gòu)造背景的判別分析建立在更嚴(yán)格的基礎(chǔ)上, Vermeesch (2006a, 2006b)收集了45種共756件已知構(gòu)造背景的、有全巖主、微量元素定量分析數(shù)據(jù)的海洋玄武巖樣品(洋島玄武巖259件、洋中脊玄武巖241件、島弧玄武巖256件), 建立了分類樹。分類和回歸樹(classification and regression tree, CART)是用一個(gè)分段常值函數(shù)采用了更嚴(yán)格的統(tǒng)計(jì)方法——線性判別分析和二次判別分析。所選擇的數(shù)據(jù)在統(tǒng)計(jì)分析前均轉(zhuǎn)換為對(duì)數(shù)比值(不是直接用該成分的含量)。判別分析是假定多變量正態(tài)性(multivariate normality), 如果所有的構(gòu)造背景分類具有相同的協(xié)方差結(jié)構(gòu), 則所確定的分類界線是線性的, 稱為線性判別分析(linear discriminant analysis, LDA), 與此相反則為二次判別分析(quadratic discriminant analysis, QDA)。用上述方法產(chǎn)生了 14190個(gè)三元判別圖解, 對(duì)它們的窮舉分析(exhaustive exploration)獲得了最好的線性和二次判別圖解, 分別是不活動(dòng)元素 Ti-V-Sc和 Ti-V-Sm三元圖解。上述方法減少了判別分析的錯(cuò)判, 判別函數(shù)是線性結(jié)合, 它使不同構(gòu)造背景分類之間相對(duì)于每一類型內(nèi)部的變化最大化。Belousova et al. (2002)根據(jù)鋯石的稀土和微量元素含量和組成, 采用遞歸分割法(recursive partitioning, 簡(jiǎn)稱RP法), 該方法關(guān)鍵是選擇節(jié)點(diǎn), 或稱根節(jié)點(diǎn)、終端節(jié)點(diǎn)、葉節(jié)點(diǎn), 進(jìn)而建立決策樹或決定樹。構(gòu)建判別不同類型巖漿巖的分類樹(classification and regression trees)。本文介紹的不同構(gòu)造?巖漿背景鋯石微量元素判別圖則是用其微量元素的二維核密度分布投影(詳見(jiàn)下述)。與上述分類樹類似, 對(duì)于一些特殊類型的火成巖建立了一組圖解, 使用這些判別圖解時(shí)需按一定流程(flow-chart)進(jìn)行, 如對(duì)于鉀質(zhì)火成巖(Müller and Groves, 1997)。

    第三是用礦物, 特別是穩(wěn)定副礦物的、而不是全巖的微量元素組成構(gòu)建構(gòu)造?巖漿背景判別圖解。

    本文重點(diǎn)介紹火成巖中一些造巖礦物和副礦物(特別是鋯石)微量元素組成在構(gòu)造背景判別中的應(yīng)用。

    1 鋯 石

    鋯石是倍受關(guān)注的遍在性礦物, 從地殼巖石到地幔包體、月巖及隕石中都有鋯石存在。鋯石還是化學(xué)上很穩(wěn)定的礦物, 在風(fēng)化、搬運(yùn)過(guò)程, 甚至在變質(zhì)和深熔過(guò)程都可保持化學(xué)上的穩(wěn)定。因此, 火成巖中鋯石U-Pb被廣泛用于定年, 特別是激光等離子體質(zhì)譜的單顆粒礦物定年技術(shù)的發(fā)展, 使鋯石成為巖漿巖U-Pb定年的最重要對(duì)象, 積累了大量年齡數(shù)據(jù)。而在鋯石用于U-Pb定年的同時(shí), 其單顆粒原位微量元素含量的準(zhǔn)確測(cè)定更提供了探討晶出鋯石的巖漿源區(qū)和成巖過(guò)程的重要地球化學(xué)信息。

    1.1鋯石的稀土及微量元素地球化學(xué)特點(diǎn)

    強(qiáng)烈富集重稀土和明顯的 Ce正異常是巖漿鋯石的典型特征, 這種特點(diǎn)主要受控于晶出鋯石的巖漿巖成分以及鋯石的稀土和微量元素分配系數(shù)。

    不同類型巖漿巖中鋯石的稀土和微量元素含量和組成有明顯區(qū)別, 如: 金伯利巖中的鋯石稀土總量低(<50×10-6), 重稀土富集程度也較低, (Yb/Sm)N= 3~30; 但花崗巖中鋯石稀土含量可達(dá)百分之幾, 重稀土富集程度很高, (Yb/Sm)N>100。因此可根據(jù)鋯石稀土元素組成特點(diǎn)區(qū)分源區(qū)的特征(表1)。

    鋯石稀土分配系數(shù)明顯與形成鋯石的巖漿成分有關(guān), 硅不飽和堿性巖中鋯石的稀土分配系數(shù)最高,而基性巖最低, 中性巖、酸性巖居中(Hanchar and Westrenen, 2007)。

    表1 不同類型火成巖鋯石的稀土元素組成Table 1 REE composition of zircon from different type of igneous rocks

    稀土元素在鋯石(ZrSiO4)中的分配系數(shù)還受控于其晶體化學(xué)(離子半徑及電荷)特點(diǎn)。鋯石中稀土元素的進(jìn)入是基于與鋯石等結(jié)構(gòu)的“磷釔礦式”(YPO4)置換:

    稀土元素中重稀土元素Lu3+離子半徑為0.0977 nm, Y3+為 0.1019 nm; 輕稀土元素 La的離子半徑為0.116 nm, 而8配位的Zr4+離子半徑為0.084 nm, 因此, 重稀土元素離子與 Zr4+離子半徑更接近, 更容易進(jìn)入鋯石晶格, 導(dǎo)致鋯石中強(qiáng)烈富集重稀土, 這明顯反映在分配系數(shù)上。Hanchar and Westrenen (2007) 綜合了自然界火山巖和熔體包裹體樣品及由高溫高壓實(shí)驗(yàn)研究所獲得的鋯石的稀土元素分配系數(shù)(圖1), 其特點(diǎn)是: 重稀土(HREE)分配系數(shù)明顯高于輕稀土(LREE), 在一個(gè)樣品中 DLu比 DLa高 1~5個(gè)數(shù)量級(jí); 不同類型巖石之間鋯石稀土分配系數(shù)差別大, 對(duì)于輕稀土變化達(dá) 4個(gè)數(shù)量級(jí), 而重稀土變化相對(duì)較小, 為2個(gè)數(shù)量級(jí)。對(duì)于其他微量元素, 重要的是U、Th、Nb和Sc。U4+、Th4+與Zr4+電價(jià)相同, 八配位Th4+離子半徑為1.0 ?, 八配位的U更大,離子半徑為1.04 ?, 因此鋯石對(duì)它們的相容性比Hf4+低。Hf、U和 Th與鋯石可形成等結(jié)構(gòu)的固溶體: HfSiO4、USiO4和ThSiO4, 但這種固溶體是非理想的、有較大不混溶間斷(Finch and Hanchar, 2003; Geisler et al., 2007)。鋯石對(duì)U和Yb的分配系數(shù)相近, 分別為254和278。Sc和Nb進(jìn)入鋯石的置換方式是:

    上述不同稀土元素分配系數(shù)之間的差異可從其離子半徑及相關(guān)晶格變化得到解釋。在鋯石中稀土元素分配系數(shù)對(duì)離子半徑的 Onuma圖解中, 從 Lu到La構(gòu)成簡(jiǎn)單拋物線(圖 1), 該曲線頂部是最佳被置換的 Zr4+的位置, 從曲線上稀土元素的排列順序可以看出, 重稀土特別是Lu最先置換Zr4+。從圖上還可以看出Ce明顯偏離拋物線, 反映在鋯石Ce的分配系數(shù)DCe明顯高于其相鄰輕稀土元素La和Pr或Nd, 形成顯著Ce正異常。這是由于在氧化條件下 Ce3+氧化成 Ce4+,它與 Zr4+電價(jià)相同, 在置換時(shí)不需電價(jià)補(bǔ)償, 而且, Ce4+的半徑(6配位為 0.087 nm, 8配位為 0.097 nm)比Ce3+更近于Zr4+, 因而比其他LREE更容易進(jìn)入鋯石晶格而形成 Ce正異常。此外, DCe偏離拋物線的程度可作為鋯石生長(zhǎng)過(guò)程中氧化狀態(tài)的指標(biāo)。

    圖1 鋯石–熔體稀土分配系數(shù)與離子半徑和熔體成分關(guān)系的Onuma圖解(引自Hanchar and Westrenen, 2007)Fig.1 Onuma diagrams showing zircon-melt REE partition coefficients as a function of REE ion radius

    1.2不同源區(qū)巖漿巖鋯石的稀土與微量元素組成

    由于鋯石對(duì)重稀土元素的強(qiáng)烈選擇性, 很多巖石中鋯石的稀土含量發(fā)生明顯重疊, 因此難以區(qū)分開(kāi)不同源區(qū)的鋯石。由上述, U與Yb在鋯石中分配系數(shù)相近, 因此, 鋯石的 U/Yb可反映鋯石晶出時(shí)熔體的特點(diǎn)。但是, Nb和Sc在鋯石中的分配系數(shù)變化范圍很大(DNb: 0.2~21; DSc: 2~70), 可達(dá)≥10倍。因此, Grimes et al. (2007, 2015)用鋯石中微量元素U/Nb比值及它們與稀土元素的比值(U/Yb, U/Nb或 Nb/Yb)進(jìn)行不同類型鋯石比較, 這可以使由溫度造成的分配系數(shù)變化導(dǎo)致的濃度改變最小化。對(duì)采自緩慢擴(kuò)張中的大西洋和南印度洋洋脊的 36個(gè)輝長(zhǎng)巖中的300個(gè)鋯石顆粒, 以及1700多個(gè)采自大陸太古宙到顯生宙花崗巖類、金伯利巖和阿拉斯加Talkeetna 島弧石英閃長(zhǎng)巖和英云閃長(zhǎng)巖巖墻的鋯石進(jìn)行了微量元素和 REE元素分析(用 SHRIMPRG離子探針?lè)治?。用 U-Yb、U/Yb-Hf或Y; Th-Yb; Th/Yb-Hf 或Y、P對(duì)上述鋯石作圖(圖 2), 可以將洋殼輝長(zhǎng)巖的鋯石與大陸巖石中的鋯石明顯區(qū)分開(kāi)。雖然對(duì)于同一樣品中的U、Th、Yb的絕對(duì)含量可以有數(shù)量級(jí)的變化, 但洋殼巖石的鋯石 U/Yb比值平均為 0.18, 大陸巖石鋯石平均為 1.07, 金伯利巖為 2.1。這種變化與含這些鋯石的全巖 U/Yb比值一致, 洋中脊玄武巖 MORB平均為 0.01~0.1 (Klein, 2003), 大陸殼平均為0.7(Rudnick and Gao, 2003), 金伯利巖為4~6(Farmer, 2003)。阿拉斯加島弧巖石的鋯石與洋殼巖石的鋯石明顯不同, 但與大陸巖石鋯石重疊, 一些海洋島弧熔巖全巖與MORB明顯重疊, 由此推測(cè)大陸溢流玄武巖和大洋島弧玄武巖中的鋯石與洋殼巖石鋯石相似。

    由上述可見(jiàn), 在構(gòu)筑判別圖解的參數(shù)選擇中, U/Yb是很重要的參數(shù), MORB中鋯石低的U/Yb比值(<0.10)反映了它們的母巖漿及MORB地幔長(zhǎng)期的不相容元素的虧損, 它明顯不同于其他構(gòu)造背景的鋯石。但U/Yb>0.1則不是唯一的, 它包括了多種構(gòu)造背景的鋯石。因此, 需采用包括 U-Nb-Sc-Yb-Gd-Ce等多種陽(yáng)離子的雙變量圖解(Grimes et al., 2015)。例如, 在 U/Yb-Nb/Yb圖解中(實(shí)質(zhì)是鋯石的 U/Nb比值), 包含了多數(shù)MORB、洋島和金伯利巖鋯石的陰影區(qū), 形成了明顯的地幔鋯石排列, 其上部邊界的Nb/Yb和U/Yb值分別為0.004, 0.02和1.0, 10。而巖漿弧和后碰撞大陸鋯石則形成了排列之外的、高U/Yb比值的與之平行的巖漿弧排列(圖3)。

    圖2 區(qū)分大陸和海洋鋯石的U/Yb-Hf (a)和U/Yb-Y (b)圖解(Grimes et al., 2007)Fig.2 Plots of U/Yb vs. Hf (a), and U/Yb vs. Y (b) for discrimination of continental and ocean crust zircons

    圖3 用U/Nb代表火成巖鋯石的構(gòu)造?巖漿源(Grimes et al., 2015)Fig.3 U/Nb proxy for tectono-magmatic source of igneous zircon

    受地幔柱影響的洋島鋯石與MORB的明顯區(qū)別是U/Yb、U/Nb和Nb/Sc比值高; 對(duì)于U/Yb>0.1的鋯石, 高 Sc/Yb可以將弧背景與低 Sc/Yb的受地幔柱影響的洋島區(qū)分開(kāi)。鋯石中Sc的富集可能反映了熔體分異過(guò)程中角閃石的卷入, 進(jìn)而反映了拉斑質(zhì)與鈣堿質(zhì)巖漿系統(tǒng)源區(qū)成分的差異。

    1.3鋯石稀土與微量元素組成對(duì)不同構(gòu)造?巖漿背景的判別

    基于上述火成巖中鋯石微量元素與其源區(qū)及成巖過(guò)程的密切和復(fù)雜關(guān)系, Grimes et al. (2015)對(duì)鋯石微量元素組成開(kāi)展了更系統(tǒng)的研究, 收集了5438個(gè)普通構(gòu)造?巖漿背景產(chǎn)出的火成巖中鋯石SHRIMP-RG分析的微量元素?cái)?shù)據(jù), 包括洋中脊 MOR(大西洋中脊MAR; 西南印度洋脊SWIR; 東太平洋洋隆EPRS)、巖漿弧、大陸碰撞帶和地幔柱影響的洋島; 時(shí)代主要為中生代?新生代。用數(shù)理統(tǒng)計(jì)方法對(duì)這些資料進(jìn)行了計(jì)算, 統(tǒng)計(jì)參數(shù)包括平均值、中值、均方差(1σ)、偏態(tài)、峰態(tài)、第一四分位數(shù)(quartile)和第三四分位數(shù)。表2中給出了不同構(gòu)造背景火成巖中鋯石微量元素含量及U、Th、Ce、Gd、Nb和Sc與Yb的比值及U/Nb比值的平均值、中值、方差和樣品數(shù)。

    根據(jù)上述統(tǒng)計(jì)資料, Grimes et al. (2015)構(gòu)筑了不同構(gòu)造?巖漿背景的鋯石微量元素判別圖解, 這些圖解包括U/Yb-Nb/Yb、Sc/Yb-Nb/Yb和 U/Yb-Sc/Yb (圖4), 圖中坐標(biāo)均為對(duì)數(shù)值(lg)。圖解中陰影區(qū)是所收集的不同構(gòu)造?巖漿背景鋯石微量元素的二維核密度分布投影。等值線分別為 50%、80%、90%和95%, 它們代表了該等值線內(nèi)樣品所占的比例。該密度投影是用二維核密度函數(shù)(kde2d)的 R程序, 即用軸平行(axis-aligned)二元正態(tài)核(bivariated normal kernel)計(jì)算, 用方格網(wǎng)(square grid)評(píng)價(jià)的二維核密度。所產(chǎn)生的構(gòu)造背景圖解中, 洋中脊、地幔柱影響的洋島及俯沖有關(guān)的弧背景明顯分布在不同區(qū)域,每一種背景中的近 80%~90%數(shù)據(jù)都投影于圖中所限定的相應(yīng)區(qū)域中。

    圖4 不同構(gòu)造?巖漿背景鋯石微量元素地球化學(xué)標(biāo)志的密度分布圖(Grimes et al., 2015)Fig.4 Density distribution plots of tectono-magmatic setting based on trace elements in igneous zircon

    在使用上述圖解時(shí)應(yīng)注意的是, Grimes et al. (2015)用鋯石微量元素成分判別構(gòu)造?巖漿背景的圖解是基于現(xiàn)代到中生代的火成巖中的鋯石, 而且對(duì)古生代火成巖的鋯石也做了檢驗(yàn), 其結(jié)果基本一致,但如果擴(kuò)展到前寒武紀(jì)應(yīng)當(dāng)非常小心。此外, 用于構(gòu)建構(gòu)造背景判別圖解的鋯石微量元素資料均來(lái)自SHRIMP-RG分析方法, 不同實(shí)驗(yàn)室之間的數(shù)據(jù)應(yīng)進(jìn)行所分析標(biāo)樣的校正。判別圖中不同區(qū)域之間界線的絕對(duì)值可能有些移動(dòng), 但其空間關(guān)系不會(huì)變化。此外, 應(yīng)綜合采用有 U、Yb、Nb、Sc、Ce和Gd建立的多變量圖解, 不要僅依據(jù)單一圖解。在鋯石微量元素分析中, 應(yīng)特別注意Nb、Sc的分析必須準(zhǔn)確定量測(cè)定, 因?yàn)檫@兩個(gè)元素容易受占鋯石 50%多的Zr的同量異位效應(yīng)的干擾。

    2 單斜輝石

    單斜輝石的化學(xué)成分隨著其寄生火成巖巖石的化學(xué)成分而變化, 對(duì)于單斜輝石的斑晶更為明顯。不同類型玄武巖的單斜輝石斑晶成分之間的差別比相應(yīng)的基質(zhì)成分之間的差別更大。該特征被用作識(shí)別不同構(gòu)造背景的判別要素, 同時(shí)也提供了判別蝕變巖石構(gòu)造背景的可靠手段, 因?yàn)槲g變巖石的單斜輝石斑晶核心的化學(xué)成分可能沒(méi)有變化(Rollinson, 1993)。Leterrier et al. (1982)對(duì) Nisbet and Pearce (1977)提出的單斜輝石MnO-TiO2-Na2O含量判別圖解進(jìn)行了修改完善, 構(gòu)筑了新的 Ti-Cr-Ca-Al-Na投影圖解, 其圖解是根據(jù)一系列更大數(shù)據(jù)庫(kù)為基礎(chǔ)的Ti、Cr、Ca、Al和Na資料而完成的(圖5)。Leterrier et al. (1982)提出的三個(gè)圖解可以區(qū)別堿性玄武巖、擴(kuò)張中心拉斑玄武巖和島弧玄武巖。單斜輝石分析數(shù)據(jù)按六個(gè)氧進(jìn)行晶體化學(xué)計(jì)算, 換算出陽(yáng)離子數(shù),并且只用那些Ca原子數(shù)大于0.5(每個(gè)單斜輝石分子中)的分析數(shù)據(jù)。為了提供計(jì)量比例, 進(jìn)行了 Fe3+和Fe2+的計(jì)算(見(jiàn)Cameron and Papike, 1981)。第一個(gè)圖解(圖 5a)利用 Ti-(Ca+Na)投影, 區(qū)分堿性玄武巖(A: 大洋島和大陸堿性玄武巖)和拉斑玄武巖及鈣堿性玄武巖(T)。第二個(gè)圖解(圖 5b)是(Ti+Cr)-Ca的投影, 區(qū)分非堿性玄武巖及非造山玄武巖(D: MORB, 大洋島拉斑玄武巖和弧后盆地拉斑玄武巖)與火山弧玄武巖(O)。第三個(gè)圖解(圖5c)是Ti-Al總的投影圖解, 區(qū)分火山弧玄武巖和鈣堿性玄武巖(C)與島弧拉斑玄武巖(I)。

    這些圖解不能只用一個(gè)單斜輝石的分析數(shù)據(jù); Leterrier et al. (1982)推薦的最少不得少于10個(gè)分析數(shù)據(jù)。如用20個(gè)分析數(shù)據(jù)投影到圖解上將會(huì)獲得更可靠的結(jié)果。此方法甚至可以應(yīng)用于綠片巖相的變基性巖,其中的單斜輝石成分可能經(jīng)受了變質(zhì)反應(yīng)的調(diào)整。

    Le Bas (1962)曾用單斜輝石的四面體中Al所占比例Al總與Ti含量關(guān)系劃分鎂鐵巖漿系列。Loucks (1990)認(rèn)為, 單斜輝石的 Al總與Ti關(guān)系可用于區(qū)分不同構(gòu)造環(huán)境的鎂鐵–超鎂鐵巖石, 其依據(jù)是在非造山區(qū)(如洋中脊, 弧后盆地、大陸裂谷、熱點(diǎn))的拉斑玄武質(zhì)和硅不飽和巖石中, 單斜輝石晶格結(jié)合的Al大部分是通過(guò)電荷補(bǔ)償VIMgIVSi2=VITiIVAl2, 以CaTiAl2O6分子形式進(jìn)入單斜輝石。這種置換方式由于在熔體中 Ti的活度比Si高, 在硅不飽和鎂鐵-超鎂鐵巖石中更為典型。在匯聚板塊邊緣(拉斑玄武巖、鈣堿性輝長(zhǎng)巖和玄武巖), 鎂鐵質(zhì)火成巖的特點(diǎn)是大部分Al以CaFe3+AlSiO6分子結(jié)合在單斜輝石中,置換方式是VIMgIVSi2?VIFe3+IVAl, 這是由于在島弧軸部水和氧的逸度高。將輝長(zhǎng)巖和超鎂鐵巖堆積巖的輝石中四面體Al對(duì)八面體中的Ti作圖, 可見(jiàn)與島弧有關(guān)的堆積巖中的單斜輝石 Al/Ti比值比與裂谷有關(guān)的拉斑玄武巖中的單斜輝石高, 約為其兩倍。因此,單斜輝石 Al/Ti比值可用于區(qū)分造山帶中蛇綠巖質(zhì)和非蛇綠巖質(zhì)的鎂鐵–超鎂鐵外來(lái)體。我們用這種方法對(duì)新疆克拉瑪依地區(qū)晚石炭世哈圖玄武巖的形成構(gòu)造背景進(jìn)行了討論(Tang et al., 2012, 圖6), 該玄武巖具有與大洋高原玄武巖相似的地球化學(xué)特點(diǎn), 其輝石的Al總對(duì)TiO2含量作圖分布于島弧和裂谷之間區(qū)域, 并從島弧趨勢(shì)偏向裂谷。結(jié)合該區(qū)早石炭世典型島弧環(huán)境及與哈圖玄武巖同時(shí)期的埃達(dá)克巖的存在, 認(rèn)為哈圖玄武巖形成于島弧和伸展的雙重背景, 與洋內(nèi)弧中洋脊俯沖形成的板片窗有關(guān)。

    圖5 玄武巖的單斜輝石斑晶對(duì)構(gòu)造背景的Ti-(Ca+Na), (Ti+Cr)-Ca和Ti-Al總判別圖解(據(jù)Leterrier et al., 1982)Fig.5 Discrimination diagrams for clinopyroxene phenocrysts in basalts

    3 尖 晶 石

    尖晶石的化學(xué)成分對(duì)于其母巖漿的地球化學(xué)非常靈敏, 是其母巖漿形成的構(gòu)造背景的函數(shù), 因此,尖晶石的Cr、Al、Fe、Mg和Ti的比值廣泛用于識(shí)別其母巖漿性質(zhì)和鉻鐵巖形成的構(gòu)造背景。然而,不同構(gòu)造背景形成的鉻鐵巖中尖晶石在 Cr-Al-Fe2+, (Cr/(Cr+Al))/(Mg/(Mg+Fe2+)), TiO2-Al2O3, TiO2-Cr2O3等圖解中出現(xiàn)大量重疊。Gonzalez-Jimenez et al. (2012)用LA-ICP-MS對(duì)已知構(gòu)造背景的鉻鐵巖中的尖晶石進(jìn)行微量元素分析, 并以MORB中的尖晶石為標(biāo)準(zhǔn)進(jìn)行作圖, 發(fā)現(xiàn)可以清楚地將科馬提巖、層狀鎂鐵質(zhì)侵入體及蛇綠巖區(qū)分開(kāi)。蛇綠巖中的尖晶石以富Ni、Mg而與虧損Ni、Mg的科馬提巖、層狀鎂鐵質(zhì)侵入體中的尖晶石區(qū)分開(kāi), 層狀鎂鐵質(zhì)侵入體中的尖晶石則明顯富集Sc、V而與科馬提巖蛇綠巖中的尖晶石相區(qū)分(圖7)。

    圖6 克拉瑪依哈圖玄武巖單斜輝石Al總-TiO2與構(gòu)造背景關(guān)系(據(jù)Tang et al., 2012修改自Loucks, 1990)Fig.6 Plot of Al總vs. TiO2in clinopyroxene for tectonic settings

    4 黑 云 母

    在Debon and Le Fort (1983)提出的侵入體化學(xué)?礦物分類中, 提出了特征礦物(characteristic mineral diagram)圖解, 圖解由參數(shù) A=Al?(K+Na+2Ca); B= (Fe+Mg+Ti)和礦物分類組成。該圖解包含了對(duì)火成巖含Al特點(diǎn)和對(duì)其礦物學(xué)的解釋, 它將火成巖劃分為鈣鐵鎂質(zhì)(cafemic, 主要或全部源于地幔)、Al質(zhì)(主要或全部源于陸殼深熔)及鈣鐵鎂質(zhì)?Al質(zhì)過(guò)渡三種類型。這種特征礦物圖解包含了巖石所固有的構(gòu)造背景信息。據(jù)此, 對(duì)堿性、鈣堿性和過(guò)鋁質(zhì)巖漿中的325個(gè)黑云母典型成分進(jìn)行了分析(Abdel-Rahman, 1994), 其典型特點(diǎn)是堿性非造山系列巖漿巖中黑云母富鐵, 為鐵云母, FeO*/MgO平均值為7.04; 過(guò)鋁質(zhì)系列巖漿巖的黑云母親鐵, FeO*/MgO平均值為 3.48;鈣堿系列黑云母富鎂, FeO*/MgO平均值為1.76。

    根據(jù)上述不同構(gòu)造背景巖漿巖中黑云母成分的顯著差異, 用FeO*-MgO, FeO*-Al2O3, Al2O3-MgO和Al2O3-MgO-FeO*分別作圖(式中 FeO*=[FeO+(Fe2O3×0.89981)]), 在這些圖解中可劃分出伸展(堿性)、俯沖有關(guān)(鈣堿性)和碰撞(過(guò)鋁質(zhì))等不同構(gòu)造背景的分布區(qū)。對(duì)這些圖解中堿性、鈣堿性和過(guò)鋁質(zhì)巖漿巖分布區(qū)進(jìn)行線性統(tǒng)計(jì)分析, 獲得了兩個(gè)線性判別函數(shù), 分別為:

    式中變量為對(duì)數(shù)比轉(zhuǎn)換(log-ratio transformation), 不是重量百分(有關(guān)對(duì)數(shù)比轉(zhuǎn)換方法原理可參考文獻(xiàn):周蒂, 1988, 1998)。用上述判別函數(shù)F1對(duì)F2作圖(圖8), 圖中劃分出A、C和P三部分, 分別代表伸展背景(堿性)、俯沖有關(guān)背景(鈣堿性)和碰撞背景(過(guò)鋁質(zhì)),相應(yīng)的判別準(zhǔn)確率分別為89%、91%和95%。該圖解與火成巖全巖主、微量元素構(gòu)造背景判別圖相結(jié)合,可以更準(zhǔn)確地判斷火成巖形成的構(gòu)造背景。

    圖7 以MORB中尖晶石為標(biāo)準(zhǔn)的不同構(gòu)造背景鉻鐵巖中的尖晶石主、微量元素圖解(據(jù)Gonzalez-Jimenez et al., 2012)Fig.7 Profiles of major and trace elements of spinel in eh chromitites from different tectonic settings normalized to composition of spinel of MORB

    5 角 閃 石

    不同構(gòu)造背景形成的角閃石對(duì)Nb、Ta的富集程度有明顯差異。Coltorfi et al. (2007)研究了來(lái)自巖石圈不同地幔區(qū)地幔包體的角閃石, 發(fā)現(xiàn)它們的微量元素,特別是高場(chǎng)強(qiáng)元素含量明顯不同: 產(chǎn)于俯沖帶上方角閃石(S 型)相對(duì)虧損 Nb(含量<10×10?6), Ti/Nb 和Zr/Nb比值高于球粒隕石, 而產(chǎn)于板內(nèi)的角閃石(I型)富集 Nb(含量在10×10?6~100×10?6), Ti/Nb, Zr/Nb低于球粒隕石。因此, 角閃石是了解俯沖帶上方流體、熔體交代特點(diǎn)的重要工具, 角閃石 Zr/Nb比值以及Zr/Nb-Ti/Nb, Zr/Nb-Ti/Zr作圖可以區(qū)分與不同構(gòu)造背景有關(guān)的、受多期交代的地幔區(qū)很有效(圖9)。

    圖 8 不同構(gòu)造背景火成巖中黑云母成分判別圖解(Abdel-Rahman, 1994)Fig.8 Discrimination diagram of two linear discriminant functions of biotite in different tectonic settings

    圖 9 區(qū)分板內(nèi)和俯沖帶上方角閃石的 Zr/Nb-Ti/Nb(a)和Zr/Nb-Ti/Zr(b)圖解(據(jù)Coltorfi et al., 2007)Fig.9 Zr/Nb vs. Ti/Nb (a) and Zr/Nb vs. Ti/Zr (b) diagrams for intraplate and suprasubduction amphiboles

    6 結(jié) 語(yǔ)

    火成巖中的礦物, 特別是穩(wěn)定副礦物微量元素地球化學(xué)的研究, 為探討火成巖(特別是蝕變的火成巖)構(gòu)造–巖漿背景提供了新的途徑, 可以預(yù)料, 隨著單顆粒礦物原位微量元素準(zhǔn)確測(cè)定技術(shù)的發(fā)展和資料的大量積累, 火成巖形成構(gòu)造背景的地球化學(xué)研究水平將迅速提高。

    趙振華. 2007. 關(guān)于巖石微量元素構(gòu)造環(huán)境判別圖解使用的有關(guān)問(wèn)題. 大地構(gòu)造與成礦學(xué), 31(1): 92–103.

    周蒂. 1988. 成分?jǐn)?shù)據(jù)的對(duì)數(shù)比統(tǒng)計(jì)方法簡(jiǎn)介. 地質(zhì)科技情報(bào), 7(2): 107–114.

    周蒂. 1998. 地質(zhì)成分?jǐn)?shù)據(jù)統(tǒng)計(jì)分析——困難與探索. 地球科學(xué), 23(2): 147–152.

    Abdel-Rahman A M. 1994. Nature of biotite from alkaline, calk-alkaline and peraluminous magmas. Journal of Petrology, 35: 525–541.

    Belousova E A, Griffin W L, O’Relly S Y and Fisher N J. 2002. Igneous zircon: Trace element composition as an indicator of source rock type. Contributions to Mineralogy and Petrology, 143: 602–623.

    Cameron M and Papike J J. 1981. Structual and chemical variations in pyroxene. American Mineralogist, 66: 1–50.

    Coltorti M, Bonadiman C, Faccini B, Grégoire M, O’Reilly S Y and Powell W. 2007. Amphiboles from suprasubduction and intraplate lithospheric mantle. Lithos, 99: 68–84.

    Debon F and Le Fort P. 1983. A chemical-mineralogical classification of common plutonic rocks and association. Transactions of the Royal Society Edinburgh: Earth Science, 73: 135–149.

    Farmer G L. 2003. Continental basaltic rocks // Holland H D and Turekian K K. Treties on Geochemistry, Vol. 3: The Crust. Oxford: Elsevier-Pergamon: 239–263

    Finch R J and Hanchar J M. 2003. Structure and chemistry of zircon and zircon-group minerals. Mineral Society America Reviews in Mineralogy and Geochemistry, 53: 1–25.

    Geisler T, Schaltergger U and Tomascher F. 2007. Re-equilibration of zircon in aqueous fluids and melts. Elements, 3: 43–50.

    Gonzalez-Jimenez J M G, Griffin W L, Locmelis M, O’Reilly S Y and Pearson N J. 2012. Contrasted minor-and trace-element compositions of spinel in chromitites of different tectonic settings. 22nd VM Goldshmidt Conference, Post.

    Grimes C B, John B E and Kelemen P B. 2007. Trace element chemistry of zircon from oceanic crust: A method for distinguishing detrital zircon provenance. Geology, 35(7): 643–646.

    Grimes C B, Wooden J L, Cheadle M J and John B E. 2015.“Fingerprinting” tectono-magmatic provenance using trace elements in igneous zircon. Contributions to Mineralogy and Petrology, 170: 46. DOI: 10.1007/ s00410-015-1199-3

    Hanchar J M and Westrenen W V. 2007. Rare earth element behavior in zircon-melt systems. Elements, 3: 37–42.

    Klein E M. 2003. Geochemistry of the oceanic crust // Holland H D and Turekian K K. Treties on geochemistry, Vol 3: The Crust. Amsterdam: Elsevier-Pogmon: 433–464.

    Le Bas M J. 1962. The role of aluminum in igneous clinopyroxenes with relation to their parentage. American Journal of Science, 260: 267–288.

    Leterrier J, Maury R C, Thonon P, Griard D and Marchal M. 1982. Clinopyroxene composition as a method of identification of the magmatic affinities of paleo-volcanic series. Earth and Planetary Science Letters, 59: 139–154.

    Li C S, Arndt N T, Tang Q Y and Ripley E M. 2015. Trace element indiscrimination diagrams. Lithos, 232: 76–83.

    Loucks R R. 1990. Discrimination of ophiolitic from nonophiolitic ultramafic-mafic allochthons in orogenic belts by the Al/Ti ratio in clinopyroxene. Geology, 18: 346–349.

    Müller D and Groves D I. 1997. Potassic Igneous Rocks and Associated Gold-copper Mineralization. Berlin: Springer: 11–40.

    Nisbet E G and Pearce J A. 1977. Clinopyroxene composition in mafic lavas from different tectonic setting. Contributions to Mineralogy and Petrology, 63: 149–160.

    Pearce J A. 2014. Immoble element fingerpringting of opgiolites. Elements, 10: 101–108.

    Rollinson H R. 1993. Using geochemical data: Evaluation, presentation, interpretation. New York: Longman Scientific Technical: 160–250.

    Rudnick R L and Gao S. 2003. Composition of the Continental Crust. Amsterdam: Elsevier: 3: 1–64.

    Tang G J, Wyman D A, Wang Q, Li J, Li Z X, Zhao Z H and Sun W D. 2012. Arthenosphere-lithosphere interaction triggered by a slab window during ridge subduction: Trace element and Sr-Nd-Hf-Os isotopeic evidence from late Carboniferous tholiites in the western Junggar area(NW China). Earth and Planetary Science Letters, 329-330: 84–96.

    Vermeesch P. 2006a. Tectonic discrimination of basalts with classification trees. Geochimica et Cosmochimica Acta, 70: 1839–1848.

    Vermeesch P. 2006b. Tectonic discrimination diagrams revisited. Geochemistry Geophysics Geosystems, 7: 1–55.

    Discrimination of Tectonic Settings Based on Trace Elements in Igneous Minerals

    ZHAO Zhenhua
    (Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou 510640, Guangdong, China)

    Tectonic settings of igneous rocks can be discriminated based on trace element geochemistry of their constitute minerals, particular the stable accessory minerals. In-situ analyses of trace element contents of minerals in igneous rocks have been compiled and applied to set up a group of discrimination diagrams for tectonic settings. Density distribution plots were created with two-dimensional kernel of trace elements in igneous zircon which present U/Yb-Nb/Yb, Sc/Yb-Nb/Yb and U/Yb-Sc/Yb diagrams. These diagrams are effective for distinguishing the mid-ocean ridge (MOR), magmatic arc (continental arc) and ocean island (and other plume-influenced) settings. Tectonic setting discrimination using clinopyroxene, spinel, biotite and amphibole are also introduced.

    tectonic setting discrimination; trace elements; zircon; clinopyroxene; spinel; biotite; amphibole

    P595

    A

    1001-1552(2016)05-0986-010

    10.16539/j.ddgzyckx.2016.05.008

    2016-03-11; 改回日期: 2016-03-28

    項(xiàng)目資助: NSFC-新疆聯(lián)合基金重點(diǎn)項(xiàng)目(U1203291)和國(guó)家自然科學(xué)基金(41273056)聯(lián)合資助。

    趙振華(1942–), 男, 研究員, 地球化學(xué)專業(yè)。Email: zhzhao@gig.ac.cn

    猜你喜歡
    單斜火成巖輝石
    新疆和巴西薔薇輝石玉的對(duì)比研究
    粉色薔薇輝石的寶石學(xué)及礦物學(xué)特征
    火成巖研磨性試驗(yàn)研究
    不同溫度、壓強(qiáng)、氧逸度條件下斜方輝石含水性的實(shí)驗(yàn)研究
    火成巖巖脈(墻)侵蝕對(duì)工作面的影響
    準(zhǔn)噶爾盆地西緣石炭系火成巖錄井綜合評(píng)價(jià)技術(shù)
    錄井工程(2017年3期)2018-01-22 08:40:24
    低溫單斜相BiPO4∶Tb3+納米晶的合成及其發(fā)光性能
    10位高速分級(jí)比較型單斜模數(shù)轉(zhuǎn)換器
    雙層厚硬火成巖破斷的力學(xué)分析
    南海新生代玄武巖中單斜輝石地球化學(xué)特征及其地質(zhì)意義
    99久久九九国产精品国产免费| 国产午夜精品一二区理论片| 99久久九九国产精品国产免费| 夜夜看夜夜爽夜夜摸| av在线观看视频网站免费| 精品久久久久久久久亚洲| 亚洲电影在线观看av| 3wmmmm亚洲av在线观看| 一级片'在线观看视频| 精品久久国产蜜桃| 在线观看免费高清a一片| 秋霞在线观看毛片| 最近最新中文字幕大全电影3| 日韩电影二区| 亚洲经典国产精华液单| 国产爱豆传媒在线观看| 亚洲怡红院男人天堂| 床上黄色一级片| 国产精品精品国产色婷婷| 欧美激情在线99| 春色校园在线视频观看| 久久久久久久大尺度免费视频| 夜夜爽夜夜爽视频| 黄色配什么色好看| 国产精品综合久久久久久久免费| 亚洲国产精品专区欧美| 91久久精品国产一区二区成人| 婷婷色综合www| 精品99又大又爽又粗少妇毛片| 亚洲三级黄色毛片| 国产免费一级a男人的天堂| 欧美日韩在线观看h| 国内揄拍国产精品人妻在线| 国产高清有码在线观看视频| 国产人妻一区二区三区在| 久久精品久久久久久久性| 成人一区二区视频在线观看| 成年版毛片免费区| 久久精品人妻少妇| 激情五月婷婷亚洲| 五月天丁香电影| 免费看av在线观看网站| 日韩一区二区视频免费看| 色网站视频免费| 在线观看av片永久免费下载| 免费电影在线观看免费观看| 69av精品久久久久久| 精品久久久久久久末码| 97热精品久久久久久| 中文字幕免费在线视频6| 精品一区二区三卡| 最近中文字幕2019免费版| 亚洲真实伦在线观看| 在线免费观看的www视频| 久久综合国产亚洲精品| 亚洲精品日本国产第一区| 亚洲国产成人一精品久久久| 亚洲欧美成人精品一区二区| 身体一侧抽搐| 久久久久久久久大av| 国产精品爽爽va在线观看网站| 嫩草影院入口| 在线免费观看不下载黄p国产| 精品久久久久久久久久久久久| 日韩在线高清观看一区二区三区| 日本黄大片高清| 3wmmmm亚洲av在线观看| 国产精品女同一区二区软件| 久久久久精品性色| 最近最新中文字幕大全电影3| 国产精品不卡视频一区二区| 一区二区三区高清视频在线| 国产 亚洲一区二区三区 | 美女国产视频在线观看| 日本av手机在线免费观看| 好男人在线观看高清免费视频| 91精品国产九色| 久久久久久久久中文| 夜夜看夜夜爽夜夜摸| 3wmmmm亚洲av在线观看| 一区二区三区四区激情视频| 日本黄色片子视频| 91狼人影院| 只有这里有精品99| 国产精品日韩av在线免费观看| 亚洲国产最新在线播放| 大又大粗又爽又黄少妇毛片口| 亚洲最大成人中文| 亚洲国产精品sss在线观看| 亚洲av一区综合| 三级毛片av免费| 天堂影院成人在线观看| 亚洲经典国产精华液单| 一级片'在线观看视频| 能在线免费看毛片的网站| 日本黄色片子视频| 午夜视频国产福利| 国产三级在线视频| 美女xxoo啪啪120秒动态图| 一本久久精品| 国产伦精品一区二区三区四那| 三级男女做爰猛烈吃奶摸视频| 97人妻精品一区二区三区麻豆| 久久久久久国产a免费观看| 亚洲va在线va天堂va国产| 建设人人有责人人尽责人人享有的 | 国产精品一区二区在线观看99 | 天天一区二区日本电影三级| 亚洲在久久综合| 中文字幕av成人在线电影| 精品人妻视频免费看| 亚洲欧洲日产国产| 国产精品熟女久久久久浪| 91精品国产九色| av在线亚洲专区| 国产真实伦视频高清在线观看| 国产成人福利小说| 高清欧美精品videossex| 大又大粗又爽又黄少妇毛片口| 老司机影院毛片| 美女主播在线视频| 久久韩国三级中文字幕| 嘟嘟电影网在线观看| 岛国毛片在线播放| 亚洲国产av新网站| 日韩制服骚丝袜av| 日日干狠狠操夜夜爽| 国产黄色视频一区二区在线观看| 伦精品一区二区三区| 成人亚洲欧美一区二区av| 亚洲怡红院男人天堂| 69人妻影院| 亚洲一区高清亚洲精品| 一级爰片在线观看| 亚洲av中文av极速乱| 精品久久久久久久久av| 特大巨黑吊av在线直播| 内射极品少妇av片p| 日韩不卡一区二区三区视频在线| 综合色丁香网| 国产乱来视频区| 又爽又黄无遮挡网站| 午夜福利视频精品| 蜜桃久久精品国产亚洲av| 只有这里有精品99| 一个人免费在线观看电影| 51国产日韩欧美| 国产激情偷乱视频一区二区| 久久久久久九九精品二区国产| 亚洲av电影在线观看一区二区三区 | 亚洲精品成人av观看孕妇| 成人一区二区视频在线观看| 亚洲激情五月婷婷啪啪| a级毛片免费高清观看在线播放| 国产亚洲一区二区精品| 色5月婷婷丁香| 男女边摸边吃奶| 99久久中文字幕三级久久日本| 青青草视频在线视频观看| 久久久久精品性色| 免费人成在线观看视频色| 欧美日韩综合久久久久久| 亚洲真实伦在线观看| 床上黄色一级片| 好男人在线观看高清免费视频| 国产亚洲av片在线观看秒播厂 | 色综合色国产| 亚洲欧美精品自产自拍| 中文精品一卡2卡3卡4更新| 午夜老司机福利剧场| 国产亚洲精品久久久com| 精品人妻一区二区三区麻豆| 三级经典国产精品| 亚洲国产av新网站| 久久99热这里只频精品6学生| 干丝袜人妻中文字幕| 午夜久久久久精精品| 亚洲在久久综合| 深夜a级毛片| 舔av片在线| 久久久久精品性色| 国产成人一区二区在线| 国产三级在线视频| 淫秽高清视频在线观看| 久久精品国产鲁丝片午夜精品| 国产精品一区二区在线观看99 | 午夜福利成人在线免费观看| 伦精品一区二区三区| 国产精品99久久久久久久久| 亚洲乱码一区二区免费版| a级毛片免费高清观看在线播放| 国产淫片久久久久久久久| 日韩 亚洲 欧美在线| 久久久久久久久久久丰满| 亚洲最大成人中文| 91久久精品电影网| 欧美日韩视频高清一区二区三区二| 啦啦啦啦在线视频资源| 日韩,欧美,国产一区二区三区| 免费观看在线日韩| 有码 亚洲区| 啦啦啦中文免费视频观看日本| 欧美 日韩 精品 国产| 精品酒店卫生间| 一级毛片 在线播放| 免费观看的影片在线观看| 18禁动态无遮挡网站| 伊人久久精品亚洲午夜| 久久99蜜桃精品久久| 国产精品爽爽va在线观看网站| 狂野欧美白嫩少妇大欣赏| 免费观看精品视频网站| 老师上课跳d突然被开到最大视频| 国产av码专区亚洲av| 亚洲av日韩在线播放| 免费电影在线观看免费观看| 麻豆成人午夜福利视频| 联通29元200g的流量卡| 久久久亚洲精品成人影院| 久久久午夜欧美精品| a级毛片免费高清观看在线播放| 精品久久久久久成人av| 2022亚洲国产成人精品| 免费看美女性在线毛片视频| 国产精品伦人一区二区| 免费黄网站久久成人精品| 大香蕉久久网| 我的女老师完整版在线观看| 亚洲精品影视一区二区三区av| 男人狂女人下面高潮的视频| 免费观看精品视频网站| 国产成人精品福利久久| 草草在线视频免费看| 最后的刺客免费高清国语| 国产久久久一区二区三区| 欧美97在线视频| 亚洲精华国产精华液的使用体验| 国产探花极品一区二区| 精品久久久久久久久久久久久| 看黄色毛片网站| 国产av码专区亚洲av| 免费观看在线日韩| www.av在线官网国产| 老司机影院毛片| 成年版毛片免费区| 寂寞人妻少妇视频99o| www.av在线官网国产| 免费观看在线日韩| 亚洲三级黄色毛片| 国产午夜精品论理片| 午夜免费男女啪啪视频观看| 欧美高清成人免费视频www| 免费大片黄手机在线观看| 国产黄频视频在线观看| 男插女下体视频免费在线播放| 在线播放无遮挡| 欧美日韩精品成人综合77777| 在线免费十八禁| 国产探花在线观看一区二区| 国产精品综合久久久久久久免费| 久久综合国产亚洲精品| 久久精品国产鲁丝片午夜精品| 亚洲av日韩在线播放| 久久久精品94久久精品| 色尼玛亚洲综合影院| 岛国毛片在线播放| 色吧在线观看| 一本久久精品| 中文字幕久久专区| 国产 一区精品| 欧美zozozo另类| 国产一区有黄有色的免费视频 | 成年人午夜在线观看视频 | 男人爽女人下面视频在线观看| 男女国产视频网站| 观看美女的网站| 欧美丝袜亚洲另类| 久热久热在线精品观看| 亚洲欧美一区二区三区国产| 尤物成人国产欧美一区二区三区| 亚洲精品一区蜜桃| 好男人在线观看高清免费视频| 国产乱来视频区| 中文在线观看免费www的网站| 99热全是精品| 国产精品国产三级国产专区5o| 久久久久久久国产电影| 亚洲精华国产精华液的使用体验| 人妻少妇偷人精品九色| 日韩 亚洲 欧美在线| 国产精品久久久久久av不卡| 六月丁香七月| 最近中文字幕高清免费大全6| 亚洲av一区综合| 欧美性感艳星| 在线观看免费高清a一片| 久久草成人影院| 色网站视频免费| 亚洲精品,欧美精品| 久久久久久国产a免费观看| 简卡轻食公司| 久久鲁丝午夜福利片| 99热网站在线观看| 99热这里只有是精品在线观看| 久久精品国产亚洲av涩爱| 精品少妇黑人巨大在线播放| 亚洲av日韩在线播放| 亚洲综合色惰| 亚洲国产精品成人综合色| 热99在线观看视频| 搡老妇女老女人老熟妇| 国产麻豆成人av免费视频| av黄色大香蕉| 97超碰精品成人国产| 岛国毛片在线播放| 日韩av在线免费看完整版不卡| 亚洲精品,欧美精品| 在线观看一区二区三区| 国产男人的电影天堂91| 久久久久九九精品影院| 啦啦啦中文免费视频观看日本| 精品久久久精品久久久| 久久久久久久大尺度免费视频| 夜夜看夜夜爽夜夜摸| 99九九线精品视频在线观看视频| 别揉我奶头 嗯啊视频| 边亲边吃奶的免费视频| 一夜夜www| 国产高清三级在线| 免费观看精品视频网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99热6这里只有精品| 色5月婷婷丁香| 最近的中文字幕免费完整| 国产色婷婷99| 欧美日韩在线观看h| 亚洲精品影视一区二区三区av| 天堂影院成人在线观看| 久久99精品国语久久久| 好男人视频免费观看在线| 乱码一卡2卡4卡精品| 亚洲精品aⅴ在线观看| 久久精品久久精品一区二区三区| 嫩草影院精品99| 日韩视频在线欧美| 国产不卡一卡二| 一级毛片我不卡| 可以在线观看毛片的网站| 成人毛片a级毛片在线播放| 啦啦啦中文免费视频观看日本| 亚洲婷婷狠狠爱综合网| 日韩在线高清观看一区二区三区| 国产一区亚洲一区在线观看| 亚洲在线观看片| 亚洲乱码一区二区免费版| 欧美xxxx黑人xx丫x性爽| 99热6这里只有精品| 免费无遮挡裸体视频| 日日啪夜夜撸| 免费观看无遮挡的男女| 亚洲av中文字字幕乱码综合| 美女大奶头视频| 好男人在线观看高清免费视频| 中文资源天堂在线| 日本wwww免费看| 一级黄片播放器| 夫妻性生交免费视频一级片| 日韩欧美精品免费久久| 免费大片黄手机在线观看| 亚洲熟女精品中文字幕| 91久久精品国产一区二区成人| 久久人人爽人人片av| av线在线观看网站| 国内精品一区二区在线观看| 2021天堂中文幕一二区在线观| 六月丁香七月| 亚洲国产日韩欧美精品在线观看| 国产探花极品一区二区| 国产精品综合久久久久久久免费| 高清欧美精品videossex| 国产精品不卡视频一区二区| 黑人高潮一二区| 麻豆精品久久久久久蜜桃| 国产综合懂色| 神马国产精品三级电影在线观看| 性色avwww在线观看| 亚洲图色成人| 国产永久视频网站| 性插视频无遮挡在线免费观看| 2021少妇久久久久久久久久久| 少妇的逼水好多| 国产永久视频网站| 免费大片18禁| 精品久久久久久久久av| 99久国产av精品| ponron亚洲| 久久99热这里只频精品6学生| 国产免费又黄又爽又色| 男人舔奶头视频| 国产精品久久久久久av不卡| 亚洲精品一二三| 亚洲综合色惰| 色播亚洲综合网| av在线亚洲专区| 国产不卡一卡二| 亚洲av免费高清在线观看| 男女那种视频在线观看| 久久99精品国语久久久| 国产 亚洲一区二区三区 | 天堂√8在线中文| 蜜桃久久精品国产亚洲av| 91午夜精品亚洲一区二区三区| 男女视频在线观看网站免费| 亚洲精品一二三| 黑人高潮一二区| 成人鲁丝片一二三区免费| 亚洲成人av在线免费| 我要看日韩黄色一级片| 99热这里只有精品一区| 日韩人妻高清精品专区| 日本免费a在线| 天堂网av新在线| 国产精品福利在线免费观看| 精品欧美国产一区二区三| 一个人看视频在线观看www免费| 日韩大片免费观看网站| 高清欧美精品videossex| 亚洲欧美成人精品一区二区| 国产三级在线视频| 亚洲最大成人手机在线| 亚洲精华国产精华液的使用体验| av免费在线看不卡| 国产精品美女特级片免费视频播放器| 亚洲经典国产精华液单| 一边亲一边摸免费视频| 国产在线一区二区三区精| 日韩精品青青久久久久久| 亚洲精品成人av观看孕妇| av在线蜜桃| 非洲黑人性xxxx精品又粗又长| 91久久精品电影网| 免费黄色在线免费观看| 三级国产精品欧美在线观看| 国产中年淑女户外野战色| 高清午夜精品一区二区三区| 老司机影院成人| 久久久久久久久久人人人人人人| 18禁在线无遮挡免费观看视频| 免费av观看视频| 亚洲精品乱久久久久久| 啦啦啦啦在线视频资源| 人妻一区二区av| 97人妻精品一区二区三区麻豆| 国产在视频线精品| 亚洲av电影在线观看一区二区三区 | 2021天堂中文幕一二区在线观| 熟女人妻精品中文字幕| 成年女人看的毛片在线观看| 国精品久久久久久国模美| 久久久色成人| 国产亚洲一区二区精品| 日韩av不卡免费在线播放| 国产精品熟女久久久久浪| 最近中文字幕2019免费版| 一级a做视频免费观看| 91精品国产九色| 免费看a级黄色片| 男女视频在线观看网站免费| 久久久a久久爽久久v久久| 日本与韩国留学比较| 久久精品国产自在天天线| 在线观看人妻少妇| 国产亚洲午夜精品一区二区久久 | 天天躁夜夜躁狠狠久久av| 最后的刺客免费高清国语| 黄色一级大片看看| 国产黄色小视频在线观看| 成人美女网站在线观看视频| 亚洲美女视频黄频| 欧美3d第一页| 成人亚洲欧美一区二区av| 久久综合国产亚洲精品| 欧美人与善性xxx| 欧美性感艳星| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩在线观看h| 国产精品美女特级片免费视频播放器| 国产精品一区二区三区四区免费观看| 成人漫画全彩无遮挡| 国产综合懂色| 精品一区二区免费观看| 日韩欧美精品免费久久| 国产精品.久久久| 嫩草影院精品99| 人人妻人人看人人澡| 精华霜和精华液先用哪个| 日韩欧美国产在线观看| 国内揄拍国产精品人妻在线| 免费看av在线观看网站| 亚洲国产精品专区欧美| 国产在视频线在精品| 久久久精品94久久精品| 国产伦精品一区二区三区四那| 国产亚洲91精品色在线| av在线老鸭窝| 国产一区亚洲一区在线观看| 少妇的逼水好多| 又粗又硬又长又爽又黄的视频| 免费看不卡的av| 中文精品一卡2卡3卡4更新| 高清在线视频一区二区三区| 国产一区亚洲一区在线观看| 熟妇人妻不卡中文字幕| 亚洲熟妇中文字幕五十中出| 国产亚洲av片在线观看秒播厂 | 日产精品乱码卡一卡2卡三| 国产老妇伦熟女老妇高清| 精品国产露脸久久av麻豆 | 建设人人有责人人尽责人人享有的 | 哪个播放器可以免费观看大片| 日韩一区二区三区影片| 我的老师免费观看完整版| 亚洲最大成人手机在线| 国产91av在线免费观看| 亚洲精品自拍成人| 亚洲人与动物交配视频| 中文字幕av成人在线电影| 午夜激情欧美在线| 久久精品人妻少妇| 久久久久精品性色| 九草在线视频观看| 日韩伦理黄色片| 一二三四中文在线观看免费高清| 日本熟妇午夜| 日本wwww免费看| 亚洲人成网站在线播| 国产成人一区二区在线| 欧美性感艳星| av在线亚洲专区| 一级二级三级毛片免费看| 黄色欧美视频在线观看| 日韩三级伦理在线观看| 最近手机中文字幕大全| 内地一区二区视频在线| 免费观看在线日韩| 性插视频无遮挡在线免费观看| 中国美白少妇内射xxxbb| 嫩草影院新地址| 岛国毛片在线播放| av一本久久久久| 在线播放无遮挡| 少妇丰满av| 久久久久网色| 岛国毛片在线播放| av播播在线观看一区| 老司机影院成人| av网站免费在线观看视频 | 亚洲欧美中文字幕日韩二区| 国产精品.久久久| 六月丁香七月| 久久久精品免费免费高清| 亚洲精品乱码久久久久久按摩| 日韩亚洲欧美综合| 欧美xxⅹ黑人| 乱系列少妇在线播放| 神马国产精品三级电影在线观看| 久久人人爽人人爽人人片va| 久久精品熟女亚洲av麻豆精品 | 精品久久久噜噜| 国产欧美另类精品又又久久亚洲欧美| 欧美精品国产亚洲| 精品一区二区三卡| 亚洲欧洲日产国产| 在线播放无遮挡| 搡老妇女老女人老熟妇| 日本熟妇午夜| 国产一区二区亚洲精品在线观看| 午夜精品一区二区三区免费看| 亚洲欧美日韩东京热| 最近中文字幕2019免费版| 亚洲av成人av| 99热这里只有是精品在线观看| 色网站视频免费| 男女边吃奶边做爰视频| 午夜老司机福利剧场| 久久国内精品自在自线图片| 内地一区二区视频在线| 久久精品熟女亚洲av麻豆精品 | 简卡轻食公司| 亚洲av不卡在线观看| 内射极品少妇av片p| 国产亚洲91精品色在线| 久久久久久久国产电影| 日韩欧美精品v在线| 国产精品嫩草影院av在线观看| 大香蕉久久网| 国产亚洲午夜精品一区二区久久 | 中文字幕久久专区| 婷婷六月久久综合丁香| 中文字幕人妻熟人妻熟丝袜美| 亚洲av.av天堂| 十八禁国产超污无遮挡网站| 五月伊人婷婷丁香| 精品久久久久久久久久久久久| 97超视频在线观看视频| 成人美女网站在线观看视频| 男女边吃奶边做爰视频| 日韩电影二区| 精品一区在线观看国产| 一本久久精品| 国产熟女欧美一区二区| 一区二区三区乱码不卡18| 不卡视频在线观看欧美| 久久6这里有精品| 国产精品一区二区在线观看99 |