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

    華北克拉通中東部新太古代晚期變質(zhì)火山巖及動力學體制

    2021-02-06 05:35:18劉樹文包涵高磊孫國正王偉郭榮榮郭博然付敬浩胡雅璐白翔胡方泱
    巖石學報 2021年1期
    關鍵詞:鎂質(zhì)安山巖火山巖

    劉樹文 包涵 高磊 孫國正 王偉 郭榮榮 郭博然 付敬浩 胡雅璐 白翔 胡方泱

    LIU ShuWen1, BAO Han1, GAO Lei1, SUN GuoZheng1, WANG Wei2, GUO RongRong3, GUO BoRan4, FU JingHao5, HU YaLu1, BAI Xiang6, 7 and HU FangYang8

    1. 北京大學造山帶與地殼演化教育部重點實驗室,地球與空間科學學院,北京 100871

    2. 中國地質(zhì)大學(北京)地質(zhì)過程與礦產(chǎn)資源國家重點實驗室,地球科學與資源學院,北京 100083

    3. 東北大學深部金屬礦山安全開采教育部重點實驗室,資源與土木工程學院地質(zhì)系,沈陽 110819

    4. 北京礦產(chǎn)地質(zhì)研究院,北京 100012

    5. 西南石油大學地球科學與技術學院,成都 610500

    6. 中國地震局地質(zhì)研究所,吉林長白山火山國家野外科學觀測研究站,北京 100029

    7. 中國地震局地震與火山災害重點實驗室,北京 100029

    8. 中國科學院地質(zhì)與地球物理研究所巖石圈演化國家重點實驗室,北京 100029

    1. MOE Key Laboratory of Orogenic Belts and Crustal Evolution, School of Earth and Space Sciences, Peking University, Beijing 100871, China

    2. State Key Laboratory of Geological Processes and Mineral Resources, School of Earth Sciences and Resources, China University of Geosciences, Beijing 100083, China

    3. MOE Key Laboratory on Safe Mining of Deep Metal Mines, Department of Geology, School of Resources and Civil Engineering, Northeastern University, Shenyang 110819, China

    4. Beijing Institute of Geology for Mineral Resources, Beijing 100012, China

    5. School of Geoscience and Technology, Southwest Petroleum University, Chengdu 610500, China

    6. National Observation and Research Station of Jilin Changbaishan Volcano, Institute of Geology, China Earthquake Administration, Beijing 100029, China

    7. Key Laboratory of Seismic and Volcanic Hazards, China Earthquake Administration, Beijing 100029, China

    8. State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China

    冥古宙和太古宙占據(jù)了地球一半以上的形成與演化歷史,其殼幔動力學體制是地球科學的關鍵問題。長期以來研究者們基于地球早期比太古宙以后更高的地幔潛能溫度和大陸長英質(zhì)地殼起源于早期玄武質(zhì)地殼的認識,提出各種各樣的太古宙殼幔動力學體制模型,其代表性模型包括:(1)滯留-蓋層模型(Stagnant-lid),該模型認為冥古宙到太古宙時期,隨著地球的冷卻,早期地殼表面玄武質(zhì)巖石漸漸冷卻形成滯留蓋層,蓋層之下不同規(guī)模的地幔對流導致對流室匯聚、沉降邊緣對接增厚向下拖曳,發(fā)生高壓變質(zhì)甚至榴輝巖化、部分熔融形成最早的TTG巖漿。與玄武質(zhì)巖漿結(jié)晶分異形成最早的長英質(zhì)地殼的認識一樣,該模型需要地殼中絕大部分的地質(zhì)記錄是拉斑玄武巖,而TTG片麻巖只沿著匯聚邊界分布,然后這一匯聚邊界后退形成更多的長英質(zhì)片麻巖(Moyen and Laurent, 2018);(2)地幔柱模型認為鐵鎂質(zhì)和超鐵鎂質(zhì)巖漿起源于上下地幔邊界或者核幔邊界未虧損地幔的高溫部分熔融,可以形成多級地幔柱體系,產(chǎn)生地球早期具有代表性的科馬提巖、科馬提質(zhì)玄武巖、苦橄巖和富集的拉斑玄武巖等地質(zhì)記錄,其中地幔柱頭部和邊部的拉斑玄武巖成分略有差別。該模型以巴布頓綠巖帶為典型代表,大量的鐵鎂質(zhì)-超鐵鎂質(zhì)巖出現(xiàn)是該模型的關鍵地質(zhì)記錄(Bédard, 2018; Zhao and Zhai, 2013)。這種模型的關鍵是超高溫和干的伸展構(gòu)造體系,基本不虧損的高度熔融的幔源巖漿,出現(xiàn)鐵鎂質(zhì)巖墻群、高溫A型鉀質(zhì)花崗質(zhì)巖石和少量殼?;旌铣梢虻闹行詭r漿巖石記錄。這種體制下長英質(zhì)巖石通常形成于干體系高溫和低壓熱動力學過程;(3)重力沉降模型(Sagduction tectonics)主要是基于太古宙的地幔溫度高于現(xiàn)代地幔溫度250℃以上,地殼巖石處于高熱和高度塑性狀態(tài),地殼表層高密度鐵鎂質(zhì)巖石因重力不穩(wěn)定而下沉,到地殼底部熔融形成TTG巖漿(Lin and Beakhouse, 2013)。這種構(gòu)造體制主要受控于地殼不同層次的巖石密度,主要下沉的是高密度的鐵鎂質(zhì)巖石和BIF,形成重力卵型底劈構(gòu)造和高角度構(gòu)造面理和線理,當沉降深度的圍巖與沉降巖石密度接近時運動停止。由于其它變質(zhì)沉積巖與圍巖的密度差很小,難于發(fā)生大規(guī)模沉降,尤其是變質(zhì)泥質(zhì)巖和長英質(zhì)巖石幾乎不可能發(fā)生重力沉降,且被重力沉降的巖石很難折返回到地表。沉降后形成一些類拗拉槽型上部沉積充填帶;(4)滴落構(gòu)造模型(Dripping tectonics)是結(jié)合滯留蓋層思路,通過熱-力學數(shù)值模擬,獲得巖石圈下部在地幔對流粘滯力帶動下物質(zhì)粘性流動,巖石圈上部保持基本穩(wěn)定的滯留層,形成巖石圈底部不規(guī)則向下突出、沉降到拆沉-滴落,導致下部巖石圈物質(zhì)再循環(huán)進入地幔,熔體通過沉降管道上升形成TTG巖漿(Capitanioetal., 2019a, b; Cawoodetal., 2018; Nebeletal., 2018),但是地表物質(zhì)尤其是長英質(zhì)沉積物很難通過這種模型進入到地幔,且這種構(gòu)造機制形成再循環(huán)型巖漿的過程需要非常長的時間周期;(5)包括我們在內(nèi)的很大一部分研究者將太古宙構(gòu)造體制,尤其是新太古代構(gòu)造體制與現(xiàn)代主體構(gòu)造體制類比,用板塊構(gòu)造體制來解釋新太古代殼幔作用的記錄。這種觀點通常認為從古太古代時就存在殼幔再循環(huán)(Wangetal., 2020)。這種殼幔再循環(huán)導致的地幔含水量和硅-堿-鋁質(zhì)組分有所增加,流體-殼源熔體交代地幔,導致其降溫并部分熔融,形成新太古代巖漿記錄的多樣性和明顯的中新太古代地幔的不均一性,且氧逸度明顯增加。這種作用的機制與現(xiàn)代板塊構(gòu)造體制有明顯差別,從中太古代晚期到新太古代晚期,以存在較為剛性的塊體之間的擠壓碰撞、熱的板片俯沖和流體-熔體交代地幔產(chǎn)生的相關殼幔巖漿作用記錄為主要標志。這種俯沖機制包括高地熱梯度條件下塊體對接(Buckling)、中等地熱梯度下小規(guī)模高角度淺俯沖然后俯沖板片翻轉(zhuǎn)和斷離(Rolling back and break-off)(Baietal., 2014, 2016; Fuetal., 2017, 2018; Guoetal., 2013, 2015a, b, 2017a, b, 2018; Kusky, 2020; Liuetal., 2002, 2004, 2019; Wangetal., 2013, 2015; Zhai and Santosh, 2011; 劉樹文等, 2018a)。近年來的研究進展表明,至少前三種模式可能從冥古宙、經(jīng)始太古代、古太古代到中太古代都存在,在不同演化階段可能不同的動力學體制都起到過主導作用。不同階段都存在多種方式的動力學體制的聯(lián)合作用,從地球早期地幔柱與滯留蓋層、與沉降構(gòu)造聯(lián)合動力學體制到中太古代滯留-蓋層或者地幔柱和板塊俯沖聯(lián)合動力學體制逐漸演化過渡為太古代晚期以熱俯沖的動力學體制為主(Capitanioetal., 2019a, b; Cawoodetal., 2018; Hawkesworthetal., 2019; Nebeletal., 2018)。

    本文整理了我們研究組未發(fā)表的和眾多同行們已經(jīng)發(fā)表的華北克拉通中東部新太古代晚期變質(zhì)火山巖資料,我們認為這些資料較好地記錄了新太古代晚期的殼幔作用過程。結(jié)合我們研究組最近關于華北克拉通中東部新太古代熱狀態(tài)和地殼厚度的研究,這些火山巖特征、成因及其反映的殼幔動力學過程對于理解新太古代晚期的殼幔動力學體制可能具有重要意義,供同行們批評和討論。

    1 新太古代晚期變質(zhì)火山巖大數(shù)據(jù)收集和分類原則

    本研究使用樣品的選取原則是:(1)樣品需具有配套的常量元素、微量元素(包括稀土元素)數(shù)據(jù);(2)必須是常量元素主要為XRF方法分析(少量濕法分析),微量元素(包括稀土元素)使用電感耦合等離子體質(zhì)譜(ICP-MS)進行分析,且分析數(shù)據(jù)符合精度要求;(3)每一地區(qū)同一批次樣品有可靠的定年數(shù)據(jù)和年齡;(4)同一批次樣品有與年齡配套的鋯石Lu-Hf同位素或者全巖Sm-Nd同位素數(shù)據(jù)進行成因限定。因此,我們共收集了華北克拉通中東部ca. 2.6~2.49Ga的有主量、微量和稀土元素分析結(jié)果的變質(zhì)火山巖樣品726件,其中232件為我們項目組尚未發(fā)表的資料,其它樣品數(shù)據(jù)引自近年來在國內(nèi)外學術刊物上已經(jīng)發(fā)表的有效分析資料。所有的726件分析樣品中的地區(qū)分布如下:膠東地區(qū)20件(Shanetal., 2015; Tangetal., 2007),吉林南部39件(Guoetal., 2016; 李承東等, 2014),遼北地區(qū)212件(Li and Wei, 2017; Pengetal., 2015, 2019; Wangetal., 2017a; 我們未發(fā)表的資料),遼南地區(qū)66件(Gaoetal., 2020b; Guoetal., 2017b; Hanetal., 2014; Zhuetal., 2015),魯西地區(qū)48件(Gaoetal., 2019, 2020a; Lietal., 2016; Pengetal., 2013; Shietal., 2019; Yuetal., 2019),阜新地區(qū)42件(Wangetal., 2011, 2015),冀東青龍南-遷安-灤縣地區(qū)50件(Lvetal., 2012; 郭榮榮等, 2014; 劉樹文等, 2018a;我們未發(fā)表的資料),青龍-上營-灑河橋-遵化地區(qū)91件(Guoetal., 2013, 2015b, 2017c; Wangetal., 2019a),冀北平泉-承德-赤城-宣化-懷安地區(qū)68件(Geetal., 2015; Liouetal., 2017; Wangetal., 2017b),五臺山-云中山地區(qū)54件(Gao and Santosh, 2019; Liuetal., 2016; Lüetal., 2006; Polatetal., 2005; Wangetal., 2004, 2014, 2019b),贊皇地區(qū)16件(Dengetal., 2013; Zhangetal., 2019),登封地區(qū)20件(Diwuetal., 2011; Zhangetal., 2018)。

    對這些變質(zhì)火山巖的分類和命名是一個很大的問題。長期以來研究者們主要是依據(jù)一些統(tǒng)計現(xiàn)代火山巖獲得的微量元素分類圖對所研究的變質(zhì)火山巖進行分類,忽略了這些變質(zhì)火山巖的主量元素特征,尤其是能夠直接反映礦物組成特征的SiO2、Al2O3、MgO、FeO(Fe2O3)、堿質(zhì)成分和TiO2等關鍵組分,導致一些命名的變質(zhì)火山巖與其實際巖石種類不完全一致。由于一些富鎂貧硅的巖石常常表現(xiàn)明顯蝕變,導致較高的燒失量,這也會導致巖石化學分類時巖石名稱的一些不確定性。鑒于目前這種研究現(xiàn)狀,本文的巖石化學分類主要參考IUGS推薦的Pearce and Robinson (2010)分類和Pearce and Reagan (2019)研究全球玻安巖時的火山巖分類,即將全巖主量巖石化學分析成分去掉燒失量,按照全球范圍內(nèi)太古宙鐵鎂質(zhì)巖石和遼東地區(qū)新太古代巖石中Fe3+:Fe2+=2:8(離子比, Pearce and Reagan, 2019; Pengetal., 2015),將全鐵成分換算成Fe2O3和FeO,然后歸一化到100%,重新核算每一個主量元素的組分含量(用下角標C加以注釋)。在此基礎上,為了大數(shù)據(jù)統(tǒng)計處理方便起見,將SiO2C≤45%的巖石命名為變質(zhì)超鐵鎂質(zhì)巖石、將45%≤SiO2C<52%且MgOC≥18%時的科馬提巖和12%≤MgOC<18%的變質(zhì)苦橄巖也歸于超鐵鎂質(zhì)巖石。45%1.3,因此選用(La/Yb)N=1.3將變質(zhì)玄武巖類劃分為稀土未分異型(相當于TH1型)和稀土分異型(相當于TH2型,包括拉斑玄武巖、鈣堿性玄武巖和二者過渡型)變質(zhì)玄武巖。變質(zhì)安山巖類則定義為在52%52%、MgOC≥8%和TiO2C<0.5%(Pearce and Reagan, 2019)。對既符合高鎂安山巖又符合玻安巖分類標準的巖石一律稱為變質(zhì)玻安巖,對除變質(zhì)玻安巖外、剩余滿足SiO2C>52%和MgOC≥5%的安山巖在統(tǒng)計分析時放進高鎂安山巖類一起討論。當SiO2C>63%時稱為變質(zhì)英安巖、流紋巖或者變質(zhì)長英質(zhì)火山巖,未細分。

    2 新太古代晚期變質(zhì)火山巖類型及巖石組合與分布

    從現(xiàn)有分析資料不難看出華北克拉通中東部新太古代晚期變質(zhì)火山巖巖石組合非常復雜,保存有很少量的變質(zhì)超鐵鎂質(zhì)巖石和變質(zhì)玻安巖記錄,最主要的巖石分布是拉斑玄武巖和玄武安山巖,其次是安山巖、英安巖和流紋巖(圖1a)。一些變質(zhì)火山巖也表現(xiàn)了非常高的總堿質(zhì)(Na2O+K2O)變化量,從小于1%到超過9%,尤其是遼北和冀東北部地區(qū)的少量玄武質(zhì)-玄武安山質(zhì)巖石表現(xiàn)了很高的總堿質(zhì)含量,落入了偏堿性巖石系列范圍(圖1a, b)。我們假設收集到的樣品分析結(jié)果沒有問題,那么這種高度變化的總堿質(zhì)含量可能一方面與原巖初始化學組分變化有關,這一點也被微量元素分類圖解上部分樣品落入堿性玄武巖、粗面玄武巖、粗面安山巖樣品所證實(圖1b);另一個重要原因就是變質(zhì)變形過程中部分樣品受到了較為強烈的堿質(zhì)流體交代作用。在MgOC-SiO2C分類圖上(圖1c),這些樣品表現(xiàn)為兩個重要的巖石組合,一是分布量最大的高鐵低鎂巖石組合,包括很少量的苦橄質(zhì)玄武巖、最大量的玄武巖和玄武安山巖-安山巖和少量的長英質(zhì)火山巖,另一個就是高鎂巖石組合,包括科馬提巖、苦橄巖、硅質(zhì)高鎂玄武巖、玻安巖和高鎂安山巖(圖1c黃色范圍)。與顯生宙玄武巖相比,所有的這些變質(zhì)鐵鎂質(zhì)和超鐵鎂質(zhì)巖石總體上表現(xiàn)為低鈦特征,只有很少量的樣品TiO2含量大于2%,尤其是高鎂巖石系列表現(xiàn)了從拉斑玄武巖到鈣堿性巖石系列的漸變過渡關系(圖1d)。

    按照上述分類命名原則,華北克拉通中東部新太古代晚期各地區(qū)的巖石組合特征如下(圖2b):膠東地區(qū)主要為稀土未分異型變質(zhì)玄武巖(30%)和稀土分異型變質(zhì)玄武巖(65%),含有少量的高鎂安山巖;吉林南部地區(qū)主要由稀土分異型變質(zhì)玄武巖(31%)、高鎂安山巖和安山巖(46%)、變質(zhì)長英質(zhì)火山巖(21%)和少量超鐵鎂質(zhì)巖石組成;遼北地區(qū)出現(xiàn)超鐵鎂質(zhì)巖石(12%),包括變質(zhì)蝕變橄欖巖、輝橄巖、蛇紋巖和異剝鈣榴巖(圖2a, c),這些巖石是堆晶巖還是洋殼殘片尚不清楚,最主要的巖石類型是變質(zhì)玄武巖和安山巖,包含稀土未分異玄武巖(14%)和稀土分異型變質(zhì)玄武巖(25%)、變質(zhì)高鎂安山巖和安山巖(44%),含有少量的變質(zhì)玻安巖和長英質(zhì)火山巖;遼南地區(qū)超鐵鎂質(zhì)巖石15%,稀土未分異和稀土分異型變質(zhì)玄武巖為32%,變質(zhì)玻安巖占9%,變質(zhì)高鎂安山巖和安山巖占27%,但變質(zhì)長英質(zhì)火山巖相對較多(17%)。魯西地區(qū)主要由稀土未分異和稀土分異型變質(zhì)玄武巖(44%)、變質(zhì)高鎂安山巖、安山巖(42%)為主,含有較多的變質(zhì)玻安巖(10%),并含有少量的變質(zhì)超鐵鎂質(zhì)巖石(4%);阜新地區(qū)缺少變質(zhì)超鐵鎂質(zhì)巖石和玻安巖,以稀土未分異和稀土分異型變質(zhì)玄武巖(50%)、變質(zhì)高鎂安山巖和安山巖(42%)為主,含有少量的變質(zhì)長英質(zhì)火山巖(7%);冀東青龍南-遷安-灤縣地區(qū)以大量的變質(zhì)長英質(zhì)火山巖(48%),稀土未分異和稀土分異型變質(zhì)玄武巖(22%)及變質(zhì)高鎂安山巖和安山巖(22%)為主,含有少量玻安巖(6%)和少量的變質(zhì)超鐵鎂質(zhì)巖石;冀東北部青龍-上營-灑河橋-遵化地區(qū)以稀土未分異和稀土分異型變質(zhì)玄武巖(58%)、變質(zhì)高鎂安山巖、安山巖(20%)和變質(zhì)長英質(zhì)火山巖(13%)為主,出現(xiàn)少量變質(zhì)玻安巖和變質(zhì)超鐵鎂質(zhì)巖石,包括變質(zhì)角閃石巖和易熔巖等;與冀東北部地區(qū)相似,冀北平泉-承德-赤城-宣化-懷安地區(qū)以稀土未分異和稀土分異型變質(zhì)玄武巖(35%)、變質(zhì)高鎂安山巖和安山巖(47%)、變質(zhì)長英質(zhì)火山巖(9%)為主,并少量變質(zhì)玻安巖和變質(zhì)超鐵鎂質(zhì)巖石;五臺山-云中山地區(qū)主要火山巖巖性為大量稀土分異型變質(zhì)玄武巖(19%)、變質(zhì)高鎂安山巖和安山巖(71%),含有少量稀土未分異變質(zhì)玄武巖、變質(zhì)玻安巖和長英質(zhì)火山巖;與五臺地區(qū)類似,贊皇地區(qū)主要火山巖巖性為稀土分異型變質(zhì)玄武巖(63%)、變質(zhì)高鎂安山巖和安山巖(37%);登封地區(qū)缺少變質(zhì)超鐵鎂質(zhì)巖石和玻安巖,主要為稀土未分異變質(zhì)玄武巖(55%)、稀土分異型變質(zhì)玄武巖(15%)和變質(zhì)高鎂安山巖(20%),含有少量變質(zhì)安山巖和長英質(zhì)火山巖。

    圖1 變質(zhì)火山巖的地球化學分類(a) IUGS火山巖TAS分類(Le Bas et al., 1986);(b)微量元素分類圖(Winchester and Floyd, 1976);(c) MgOC-SiO2C分類圖(Pearce and Reagan, 2019);(d) La-Yb巖石系列分布圖(Ross and Bédard, 2009). LSB-低硅玻安巖;HSB-高硅玻安巖;SHMB-硅質(zhì)高鎂玄武巖;HMA-高鎂安山巖;A-安山巖;D-英安巖;Ol-橄欖石;Opx-斜方輝石Fig.1 Geochemical classification diagrams for metavolcanic rocks(a) International Union of Geological Sciences (IUGS) total alkalis-silica (TAS) diagram of volcanic rocks (Le Bas et al., 1986); (b) Zr/TiO2×0.0001 vs. Nb/Y diagram (Winchester and Floyd, 1976); (c) MgOC-SiO2C diagram (Pearce and Reagan, 2019); (d) La vs. Yb diagram (Ross and Bédard, 2009). LSB, low-Si boninite; HSB, high-Si boninite; SHMB, siliceous high-Mg basalt; HMA, high-Mg andesite; A, andesite; D, dacite; Ol, olivine; Opx, orthopyroxene

    從上述各地區(qū)不同巖石類型組合和各類巖石分布量比(圖2)不難看出,新太古代變質(zhì)超鐵鎂質(zhì)巖石在遼北、遼南、冀東、冀北和吉南地區(qū)均有少量分布,所占火山巖樣品總量比僅有~7%,稀土未分異型玄武巖占~14%,稀土分異型玄武巖占~27%,變質(zhì)玻安巖只少量分布(~4%)。而高鎂安山巖和安山巖比例為~38%,變質(zhì)長英質(zhì)火山巖僅占~10%。

    魯西、遼南、遼北、冀東北部、冀東南部(遷安-灤縣)、冀北和五臺-云中地區(qū)可大致劃歸為同一種巖石組合類型,以稀土未分異型玄武巖(REE-undifferentiated basalts)、稀土分異型玄武巖(REE-differentiated basalts)、高鎂安山巖(high-Mg andesites)、安山巖(andesites)組合為主體,出現(xiàn)少量超鐵鎂質(zhì)巖(ultramafic rocks)、玻安巖(boninites)和長英質(zhì)火山巖(felsic volcanic rocks),我們以Ultramafic rock+Basalt+Boninite+Andesite作為其標志性巖石組合,稱為UBBA型(圖2),最突出的特點是出現(xiàn)變質(zhì)玻安巖,且出現(xiàn)大量的稀土未分異和稀土分異型玄武巖,高鎂和低鎂安山巖,并組合有少量英安巖,存在超鐵鎂質(zhì)巖石。膠東、登封和阜新地區(qū)主要保存了大量的稀土未分異、稀土分異型玄武巖和不等量安山巖,可以稱為Tholeiite-Basalt+Andesite(TBA)型組合(圖2b)。吉林南部和贊皇地區(qū)巖石組合為稀土分異型拉斑玄武巖、高鎂和低鎂安山巖,只有少量的長英質(zhì)火山巖,我們可以稱其為Basalt-Andesite(BA)型組合(圖2b)。

    圖2 華北克拉通中東部陸塊變質(zhì)火山巖區(qū)域分布及巖石組合特征(a)華北克拉通地質(zhì)簡圖(據(jù)Zhao et al., 2005修改);(b)不同區(qū)域變質(zhì)火山巖巖石組合. DF-登封;FX-阜新;JD-膠東;N-EH-冀東北部;NH-冀北;NL-遼北;S-EH-冀東南部(青龍南部-遷安-灤縣地區(qū));SJ-吉南;SL-遼南;TH-太華;WS-魯西;WT-YZ-五臺-云中;ZH-贊皇;ZT-中條Fig.2 Regional distributions and lithological assemblages of metamorphic volcanic rocks in the basement terranes of Central and Eastern Block of NCC(a) geological sketch map of the North China Craton (modified after Zhao et al., 2005); (b) the lithological assemblages of metamorphic volcanic rocks in different areas. DF, Dengfeng; FX, Fuxin; JD, eastern Shandong; N-EH, northern part of eastern Hebei; NH, northern Hebei; NL, northern Liaoning; S-EH, southern part of eastern Hebei (including Qianan, luanxian and southern Qinglong areas); SJ, southern Jilin; SL, southern Liaoning; TH, Taihua; WS, western Shandong; WT-YZ, Wutai and Yunzhong areas; ZH, Zanhuang; ZT, Zhongtiao

    3 變質(zhì)火山巖成因

    為了進一步探討這些新太古代晚期變質(zhì)火山巖所代表的動力學意義,我們必須了解各類巖石組合的巖漿作用特征和大致巖石成因。

    3.1 UBBA變質(zhì)火山巖

    這種類型變質(zhì)火山巖在華北克拉通不僅發(fā)育量大而且分布范圍廣,巖石類型復雜。該類火山巖中含48個超鐵鎂質(zhì)巖石樣品,包括蛇紋巖、蛇紋石化輝橄巖、富鐵的易熔巖、富鎂的透閃石角閃石巖、普通角閃石巖、透輝石巖,水鋁榴石異剝鈣榴巖等,因為各種類型的蝕變改造和交代作用導致該類巖石巖石學特征異常復雜。這些超鐵鎂質(zhì)巖石表現(xiàn)了類似于SSZ型蛇綠巖的巖石組合特征。Pengetal. (2015)通過大量的巖石學、巖石地球化學和Sr-Nd同位素研究認為其中超鐵鎂質(zhì)巖石直接來源于地幔,其它相關的火山-侵入巖形成于俯沖相關的弧巖漿作用。按照稀土元素特征,這些超鐵鎂質(zhì)巖石樣品中10個具有稀土未分異型拉斑玄武巖的稀土配分樣式,大部分樣品具有70~90的Mg#值,一些富鐵巖石Mg#變化范圍較寬。個別樣品表現(xiàn)出Ce負異常(圖3a, Pengetal., 2015),一些大離子親石元素和P濃度變化較大并呈不明顯的Nb和Ti負異常(圖3b),巖石成因研究表明其巖漿來源于原始地幔到虧損地幔較高程度的部分熔融,且地幔源區(qū)受到了來自俯沖相關的沉積物的熔體交代作用(圖3e, f),巖石受到了海底流體或者后期流體的強烈交代。另一類超鐵鎂質(zhì)巖(38個樣品)則表現(xiàn)為明顯右斜式稀土圖譜和Th、Nb、Ta的負異常,Sr、P和Ti從正異常到明顯負異常(圖3c, d),大部分樣品TiO2含量小于2%。巖石成因研究表明該類巖石的巖漿起源于原始地幔較高程度部分熔融,源區(qū)明顯受到了俯沖板片熔體的交代,并存在OIB型地幔端元,只是貢獻量很小(圖3e, f)。

    圖3 超鐵鎂質(zhì)巖石的微量元素特征和成因稀土未分異型超鐵鎂質(zhì)巖石球粒隕石標準化稀土圖譜(a)和原始地幔標準化多元素圖譜(b);稀土分異型超鐵鎂質(zhì)巖石球粒隕石標準化的稀土圖譜(c)和原始地幔標準化多元素圖譜(d);球粒隕石和原始地幔值引自Sun and McDonough (1989);(e) La/Yb-Nb/Yb巖石成因鑒別圖,反映原始巖漿的地幔端元和熔體交代作用(據(jù)Pearce, 2008修改);(f) Th/Yb-Nb/Yb值圖(Pearce, 2008). DM虧損地幔;N-MORB正常洋中脊玄武巖;E-MORB富集洋中脊玄武巖;CC大陸地殼;OIB洋島玄武巖;AC平均地殼Fig.3 Chondrite-normalized REE patterns and primitive mantle-normalized spider diagrams, and petrogenetic discrimination diagrams for the ultramafic rocksChondrite-normalized REE patterns (a) and primitive mantle-normalized spider diagrams (b) for REE-undifferentiated ultramafic rocks; chondrite-normalized REE patterns (c) and primitive mantle-normalized spider diagrams (d) for REE-differentiated ultramafic rocks; the chondrite and primitive mantle values are after Sun and McDonough (1989); (e) La/Yb vs. Nb/Yb diagram (modified after Pearce, 2008); (f) Th/Yb vs. Nb/Yb diagram (Pearce, 2008). DM, depleted mantle; N-MORB, normal-type mid-ocean ridge basalt; E-MORB, enriched-type mid-ocean ridge basalt; CC, continental crust; OIB, ocean island basalt; AC, average crust

    該巖石組合中的稀土未分異型高鎂拉斑玄武巖表現(xiàn)輕稀土虧損到平坦型稀土配分模式,絕大部分巖石沒有明顯的Nb、Ta、Ti負異常(圖4a, b),巖石成因研究表明該類巖石巖漿起源于尖晶石橄欖巖熔融程度在5%~30%范圍內(nèi)的部分熔融,與稀土未分異的超鐵鎂質(zhì)巖石有明顯一致的巖石成因特征,且其巖漿源區(qū)同樣受到了俯沖沉積物熔體的交代作用(圖4e, f)。與稀土未分異型高鎂拉斑玄武巖相比,稀土分異型高鎂玄武巖呈現(xiàn)明顯高的稀土總量、(La/Yb)N比值和右斜式稀土圖譜,從不明顯到明顯的Nb、Ta和Ti負異常(圖4c, d),TiO2在0.2%到2.1%范圍內(nèi)變化,清楚地表現(xiàn)了島弧型拉斑玄武巖的地球化學特征。巖石成因研究表明該類巖石巖漿起源于兩個初始巖漿端元,一個和稀土未分異型巖漿源完全一致,但是主要受到了俯沖板片流體、熔體交代的地幔尖晶石橄欖巖的部分熔融。另一個端元巖漿為OIB型巖漿源,源區(qū)物質(zhì)中石榴石和尖晶石的比例超過50:50,反映了其來源深度相對較大,該類巖漿源同樣受到了來自板片熔體的交代作用或者形成于兩者的熔體混合(圖4e, f)。

    圖4 變質(zhì)玄武巖類的微量元素特征和巖石成因稀土未分異型玄武巖球粒隕石標準化稀土圖譜(a)和原始地幔標準化多元素圖譜(b) (其中粉紅色圖標為該類樣品中高鎂玄武巖樣品);稀土分異型玄武巖球粒隕石標準化稀土圖譜(c)和原始地幔標準化多元素圖譜(d) (其中深藍色圖標為該類樣品中高鎂玄武巖樣品);球粒隕石和原始地幔值引自Sun and McDonough (1989);(e) Sm/Yb-Sm玄武質(zhì)巖漿巖石成因鑒別圖(Aldanmaz et al., 2000);(f)兩類玄武巖源區(qū)成分和俯沖熔體交代特征(據(jù)Pearce, 2008修改). PM-原始地幔;其余英文縮寫含義同圖3Fig.4 Chondrite-normalized REE patterns and primitive mantle-normalized spider diagrams, and petrogenetic discrimination diagrams for the meta-basaltsChondrite-normalized REE patterns (a) and primitive mantle-normalized spider diagrams (b) for REE-undifferentiated basalts (blue) and REE-undifferentiated high-Mg basalts (pink); chondrite-normalized REE patterns (c) and primitive mantle-normalized spider diagrams (d) for REE-differentiated basalts (red) and REE-differentiated high-Mg basalts (dark blue); the chondrite and primitive mantle values are after Sun and McDonough (1989); (e) Sm/Yb vs. Sm diagram (Aldanmaz et al., 2000); (f) La/Yb vs. Nb/Yb diagram (modified after Pearce, 2008). PM, primitive mantle; the other abbreviations are the same as Fig.3

    與高鎂拉斑玄武巖一樣,低鎂拉斑玄武巖同樣也分為稀土未分異和稀土分異型拉斑玄武巖,它們的微量元素特征與對應的高鎂拉斑玄武巖特征類似,主要形成于對應的稀土未分異和稀土分異型高鎂拉斑玄武巖的初始巖漿結(jié)晶分異,少量樣品反映了大陸地殼混染的地球化學特征,反映了它們形成演化過程中存在同化混染和分離結(jié)晶(AFC)作用。與高鎂拉斑玄武巖相比,低鎂拉斑玄武巖類缺少明確的OIB端元,但是存在E-MORB到OIB之間的樣品分布,不能完全排除OIB端元巖漿作用的存在。

    少量玻安巖(28個樣品)和大量高鎂安山巖的出現(xiàn)是該類巖石組合的標志性特征。這些巖石表現(xiàn)為高鎂含量(圖5a, b)和稀土分異型拉斑玄武巖的稀土圖譜特征,呈Nb、Ta和Ti負異常和部分樣品的Th負異常(圖5c-f),屬于鈣堿性巖石系列,只有少數(shù)樣品表現(xiàn)為稀土未分異型的稀土圖譜特征,表現(xiàn)為拉斑玄武巖系列(圖5c)。其中高鎂安山巖屬于低溫高鎂安山巖(圖5a)。低鎂安山巖表現(xiàn)為與高鎂安山巖一致的微量元素地球化學特征(圖5g, h)。在La/Sm對La巖石成因鑒別圖(圖5i)上,變質(zhì)玻安巖樣品表現(xiàn)了高的正斜率線性分布,變質(zhì)高鎂安山巖樣品表現(xiàn)了與變質(zhì)玻安巖樣品類似的正斜率線性分布但有少數(shù)樣品呈近水平分布,而低鎂安山巖較多的樣品呈近水平分布,因此變質(zhì)玻安巖和高鎂安山巖主要為部分熔融成因,變質(zhì)低鎂安山巖主要為結(jié)晶分異成因(圖5i)。在La/Yb對Nb/Yb的對數(shù)坐標圖上,所有變質(zhì)玻安巖、低鎂安山巖樣品和絕大部分高鎂安山巖樣品都表現(xiàn)了俯沖板片熔體交代和島弧型火山巖的分布趨勢,部分高鎂安山巖樣品表現(xiàn)了俯沖沉積物熔體交代的分布趨勢(圖5j)。變質(zhì)玻安巖樣品表現(xiàn)了與IBM型島弧區(qū)近似的地球化學分布,在MgOC對SiO2C分類、成因和構(gòu)造環(huán)境鑒別圖上,所有的玻安巖樣品都落在了方輝橄欖巖低壓下部分熔融的熔體范圍內(nèi),并且落在了Setouchi proto-arc、arc-basin體系和典型太古代Whudu型島弧區(qū)的巖石成分分布區(qū)內(nèi),反映了它們的低壓下方輝橄欖巖部分熔融的成因和匯聚型板塊邊緣的構(gòu)造背景(圖5b)。其中拉斑系列巖石高鎂安山巖落在了接近DM和N-MORB范圍,體現(xiàn)了它們與高鎂拉斑玄武巖的近緣關系,并且其源區(qū)受到了俯沖板片流體-熔體的交代作用,落在島弧-盆地構(gòu)造環(huán)境(圖5b),與澳大利亞太古代Whundu型玻安巖一致。低鎂安山巖類巖石主要形成于各類變質(zhì)玄武巖和高鎂安山巖巖漿的分離結(jié)晶,即可能存在AFC巖漿演化過程。

    圖5 變質(zhì)玻安巖和各類安山巖的地球化學特征(a)變質(zhì)安山巖分類圖(據(jù)鄧晉福等, 2018修改);(b)玻安巖成因分類-成因圖(Pearce and Reagan, 2019);變質(zhì)玻安巖球粒隕石標準化稀土圖譜(c)和原始地幔標準化多元素圖譜(d);變質(zhì)高鎂安山巖球粒隕石標準化稀土圖譜(e)和原始地幔標準化多元素圖譜(f);變質(zhì)低鎂安山巖球粒隕石標準化稀土圖譜(g)和原始地幔標準化多元素圖譜(h);球粒隕石和原始地幔值引自Sun and McDonough (1989);(i) La/Sm-La巖石成因鑒別圖(Treuil and Joron, 1975);(j)巖漿源區(qū)成分和俯沖熔體交代特征(據(jù)Pearce, 2008修改).縮寫同圖1和圖3Fig.5 Geochemical characteristics of the meta-boninites and andesites(a) MgOC vs. SiO2C diagram for the andesites (modified after Deng et al., 2018); (b) MgOC vs. SiO2C diagram for the boninites (Pearce and Reagan, 2019); Chondrite-normalized REE patterns and primitive mantle-normalized spider diagrams for meta-boninites (c, d), for meta high-Mg andesites (e, f) and for meta low-Mg andesites (g and h); the chondrite and primitive mantle values are after Sun and McDonough (1989); (i) La/Sm vs. La petrogenetic discrimination diagram (Treuil and Joron, 1975); (j) La/Yb vs. Nb/Yb diagram (modified after Pearce, 2008). The abbreviations are the same as Fig.1 and Fig.3

    該巖石組合中還有少量的變質(zhì)英安巖-流紋巖樣品,其中絕大部分樣品表現(xiàn)了較高的(La/Yb)N比值(>25)。這些樣品分布在冀北、冀東北部和遼北地區(qū),Mg#值在20~70之間變化,其中符合埃達克巖地球化學特征的樣品很少,主要分布在遼北地區(qū),形成于俯沖板片部分熔融,且熔體上升過程中受到了地幔楔物質(zhì)的混染(Pengetal., 2015)。絕大部分樣品盡管Mg#較高,變化范圍大,但是Mg#值與較低濃度的Cr和V等過渡族元素沒有明顯相關性,說明它們起源于下部地殼鐵鎂質(zhì)巖石的部分熔融,經(jīng)歷了分離結(jié)晶作用(Gaoetal., 2019, 2020a, b; Guoetal., 2013, 2015b, 2017c; Li and Wei, 2017; Liuetal., 2012; Lüetal., 2006; Pengetal., 2013; Wangetal., 2004)。

    3.2 TBA型變質(zhì)火山巖

    該類變質(zhì)火山巖主要分布在膠東地區(qū)、登封地區(qū)和阜新地區(qū),該類巖石組合最突出的特征是稀土未分異型和稀土分異型拉斑玄武巖占有絕對優(yōu)勢,組合有少量安山質(zhì)巖石和長英質(zhì)火山巖(圖2)。其中稀土未分異型拉斑玄武巖表現(xiàn)為近于平坦的稀土配分模式,沒有明顯的Eu異常。在原始地幔標準化微量元素圖譜中,大部分元素表現(xiàn)出較為平坦的圖譜特征,沒有Nb、Ta和Ti負異常(Diwuetal., 2011; Zhangetal., 2018; 王偉等, 2015),其La/Yb、Th/Yb和Nb/Yb比值均接近于N-MORB,并表現(xiàn)為N-MORB-E-MORB-OIB排列趨勢,指示其地幔源區(qū)沒有受到明顯的外來物質(zhì)的擾動,形成于尖晶石二輝橄欖巖10%~15%的部分熔融,即軟流圈地幔的低壓部分熔融(圖4a, b, e, f)。

    稀土分異型拉斑玄武巖,具有輕微右斜式稀土配分模式,沒有明顯的Eu異常,表現(xiàn)為微弱的Nb、Ta和Ti的負異常,具有相對較高的La/Yb和Nb/Yb比值,表現(xiàn)了拉斑玄武巖向鈣堿性玄武巖過渡的特征(圖1d、圖4c, d)。它們形成于受到俯沖流體-熔體交代富集,并有不同比例石榴石參與的虧損地幔石榴石尖晶石二輝橄欖巖的部分熔融,形成深度應該高于稀土未分異型拉斑玄武巖(圖4e, f)。它們的原巖可能形成于洋內(nèi)初始俯沖作用階段,主要起源于受到微弱俯沖流體作用改造的虧損-弱富集大洋巖石圈地幔的部分熔融(Diwuetal., 2011; Zhangetal., 2018; 王偉等, 2015)。

    本巖石組合中的高鎂和低鎂安山巖,表現(xiàn)為較寬范圍變化的總堿質(zhì)含量(圖1a),輕重稀土強烈分餾的稀土配分模式,無明顯Eu異常,通常(La/Yb)N比值較高,在2.5~25.2之間變化。在原始地幔標準化微量元素圖譜中,它們具有明顯Nb、Ta和Ti負異常,具有明顯高的La/Yb、Nb/Yb、Th/Yb和Nb/Y比值(圖5e-j),以至于部分樣品落入堿性玄武巖和粗面安山巖的范圍(圖1a, b)。這些特征表現(xiàn)了與贊岐巖類巖石類似的地球化學特征。一些高鎂安山巖類主要形成于受俯沖板片流體、俯沖板片和沉積物熔體交代的地幔楔的部分熔融,而低鎂安山巖的一部分與高鎂安山巖結(jié)晶分異有關,另一部分表現(xiàn)了與稀土分異型拉斑玄武巖的親緣性,經(jīng)歷了AFC演化過程(圖5i, j;Diwuetal., 2011; Zhangetal., 2018; 王偉等, 2015)。

    本巖石組合中存在很少量的變質(zhì)英安巖,它們表現(xiàn)出明顯高的Sr/Y(>300)和(La/Yb)N(>20)比值,及明顯高的MgO含量和Mg#值,輕稀土強烈分餾,重稀土相對分餾較弱,并呈明顯的正Eu異常。這些特征表明這類巖石形成過程中可能有幔源物質(zhì)的卷入,或者形成于俯沖板片熔體受到了地幔物質(zhì)污染,或者形成于高鎂安山質(zhì)巖漿的結(jié)晶分異。

    3.3 BA型變質(zhì)火山巖

    就目前已有的資料看,該變質(zhì)火山巖組合最重要的特征是缺少稀土未分異型變質(zhì)火山巖組合,保存了大量稀土分異型拉斑玄武巖、大量高鎂和低鎂安山巖及少量的英安巖(圖1a、圖2)。其中吉林南部夾皮溝地區(qū)變質(zhì)高鎂安山巖形成于2.49Ga,具有與日本Stouchi地區(qū)經(jīng)典的高鎂安山巖一致的高鎂鈣堿性地球化學特征和類似于贊岐巖的微量元素地球化學特征,具有高的Mg#值(68~71)和Cr、Ni、Co等過渡族元素含量,虧損Nb、Ta和Ti等高場強元素,形成于板片流體-熔體和俯沖沉積物熔體交代的地幔楔物質(zhì)的部分熔融(李承東等, 2014)。變質(zhì)稀土分異型拉斑玄武巖形成于2.59~2.54Ga,Mg#的范圍為38~64,具有虧損Nb、Ta和Ti的地球化學特征和較寬的鋯石Hf同位素范圍,并保存了古老的~2.7Ga的捕獲鋯石,形成于俯沖流體-熔體交代地幔楔物質(zhì)部分熔融過程,且經(jīng)歷了AFC演化過程(Guoetal., 2016)。與吉林南部新太古代晚期變質(zhì)火山巖類似,贊皇地區(qū)新太古代晚期變質(zhì)火山巖以稀土分異型變質(zhì)拉斑玄武巖和安山質(zhì)巖石為主體,表現(xiàn)為Nb、Ti等高場強元素虧損,形成于大洋島弧和弧前火山作用(Dengetal., 2013)。在前人分析數(shù)據(jù)基礎上,本文對這些巖石的進一步巖石成因分析表明,其初始未演化玄武質(zhì)巖漿起源于俯沖板片流體-熔體交代的地幔楔尖晶石石榴石二輝橄欖巖5%~10%的部分熔融,并經(jīng)歷了分離結(jié)晶和大陸地殼物質(zhì)的混染過程,形成了玄武質(zhì)、玄武安山質(zhì)和安山質(zhì)巖石。該殼幔巖漿作用發(fā)生于島弧和活動大陸邊緣構(gòu)造背景,與前人的認識基本一致。

    4 新太古代動力學體制

    在太古宙滯留蓋層、地幔柱還是板塊構(gòu)造體制研究中,太古宙不同階段與現(xiàn)代地幔溫度差(ΔT)是關鍵。當ΔT<100℃以內(nèi)產(chǎn)生規(guī)模較大的類現(xiàn)代板塊俯沖構(gòu)造體制;當100℃≤ΔT<175℃時產(chǎn)生規(guī)模較小、高角度的淺俯沖,然后俯沖板片發(fā)生翻轉(zhuǎn)和斷離,導致頻繁發(fā)生板片俯沖和后退,形成以側(cè)向增生為主導的殼幔動力學體制;當175℃≤ΔT<250℃時形成前俯沖構(gòu)造體制(Pre-subduction),巖石圈以塑性變形增厚和對接為主,形成擠壓增厚的構(gòu)造體制;當ΔT≥250℃不可能形成俯沖構(gòu)造,以地幔柱、滯留層和重力沉降(Sagduction or dripping tectonics)為主(Capitanioetal., 2019a, b; Cawoodetal., 2018; Gerya, 2014; Nebeletal., 2018; Sizovaetal., 2010; Van Hunen and Moyen, 2012)。

    4.1 變質(zhì)火山巖指示的動力學體制

    華北克拉通中東部新太古代變質(zhì)火山巖的少量全巖Sm-Nd的同位素資料(Wuetal., 2005)和大量的鋯石Lu-Hf同位素資料表明,新太古代經(jīng)歷了~2.5Ga和2.7Ga兩期主要的地殼生長事件,鐵鎂質(zhì)巖漿起源于虧損地幔,大部分樣品都不同程度地受到了老地殼物質(zhì)混染的影響(Gaoetal., 2019, 2020a, b; Guoetal., 2013, 2015a, b, 2016, 2017b, 2018; Liuetal., 2011; Wangetal., 2013, 2015; 萬渝生等, 2015; 王偉等, 2015)。本文和前人關于華北克拉通中東部新太古代變質(zhì)火山巖的巖石成因研究表明,這些來源于地幔的巖漿作用受到了不同程度的俯沖流體、板片熔體和沉積物熔體交代。本文和Pengetal. (2015)研究表明華北克拉通中東部新太古代晚期變質(zhì)火山巖中包括很少量超鐵鎂質(zhì)巖石樣品,主要是蛇紋石化的輝橄巖、異剝鈣榴巖、易熔巖和各種角閃石巖。超鐵鎂質(zhì)和鐵鎂質(zhì)巖石,按照Pearce and Reagan (2019)分類,其中6個樣品屬于科馬提巖,15個樣品屬于變質(zhì)苦橄巖,9個為苦橄質(zhì)玄武巖,此外還有3個MgO含量超過30%樣品在硅質(zhì)高鎂玄武巖范圍(蛇紋石化輝橄巖)。這些鐵鎂質(zhì)和超鐵鎂質(zhì)巖石(占火山巖樣品總量的4.6%)主要分布在遼北、遼南、冀東北部青龍-上營-遵化等地區(qū)。

    為了確定新太古代晚期地幔巖漿作用的特點,我們主要用MgOC≥8%的樣品討論其巖漿起源,以盡量減少結(jié)晶分異和地殼混染作用的影響。結(jié)果顯示初始巖漿有四個端元,最主要的是稀土未分異型端元和稀土分異型拉斑玄武巖端元(圖1d),以及少量的玻安巖端元和OIB型玄武巖端元。其中OIB型玄武巖樣品表現(xiàn)出較高的TiO2含量和Nb/Yb比值?!?%的典型玻安巖和大量高鎂玄武巖類樣品的存在,表現(xiàn)為明顯的Nb、Ta和Ti虧損的大量稀土未分異和稀土分異型玄武巖,巖石成因研究揭示這些巖石主要形成于俯沖板片流體-熔體和沉積物熔體交代地幔的部分熔融,而不是殼幔巖漿混合作用(Tangetal., 2007; Wangetal., 2017b)。這些火山巖成分和成因表明在新太古代晚期板片俯沖相關的殼幔巖漿作用已經(jīng)在華北克拉通中東部地區(qū)成為主導的動力學體制,但是明顯存在多元化幔源巖漿作用,尤其是OIB巖漿作用、科馬提質(zhì)玄武巖和苦橄質(zhì)玄武巖類的存在表明高溫-高壓地幔巖漿作用在局部是存在的,并且與俯沖板片體制下的巖漿作用相伴產(chǎn)出,最有可能是板塊構(gòu)造體制和地幔柱構(gòu)造體制聯(lián)合作用的產(chǎn)物(Gaoetal., 2019; Polat, 2009)。

    圖6 華北克拉通中東部新太古代地殼厚度和地熱梯度P-T空間圖(底圖據(jù)Brown, 2014修改)Lg-低程度變質(zhì), 包括沸石相;Gr-綠片巖相;Am-角閃巖相;Gn-麻粒巖相;UHTM-超高溫變質(zhì)作用;Bl-藍片巖相;E-HPG-榴輝巖-高壓麻粒巖相;UHPM-超高壓變質(zhì)作用;DF-登封;EH-冀東;JD-膠東;NL-遼北;SJ-吉南;TH-太華;WL-遼西;WS-魯西;ZT-中條Fig.6 P-T space for Neoarchean crustal thickness and geothermal gradient of Central and Eastern Block of North China Craton (base map modified after Brown, 2014)Lg, low-grade metamorphism, includes the zeolite facies; Gr, greenschist facies; Am, amphibolite facies; Gn, granulite facies; UHT, ultrahigh temperature metamorphism; Bl, blueschist facies; E-HPG, eclogite to high-pressure granulite facies; UHPM, ultrahigh pressure metamorphism; DF, Dengfeng; EH, eastern Hebei; JD, eastern Shandong; NL, northern Liaoning; SJ, southern Jilin; TH, Taihua; WS, western Shandong; WL, western Liaoning; ZT, Zhongtiao

    4.2 新太古代晚期地殼厚度和熱狀態(tài)

    近年來數(shù)值模擬實驗(Number experiments of thermal mechanics)技術的發(fā)展為太古宙殼幔動力學體制研究注入了活力。這些實驗基于假設的太古宙熱狀態(tài)參數(shù)獲得的模擬結(jié)果,提供了動力學體制、構(gòu)造樣式和動態(tài)演進模式。到目前為止,早前寒武紀地球科學研究的基本假設是太古宙地幔溫度和莫霍面地熱梯度比現(xiàn)代高得多(Gerya, 2014; Sizovaetal., 2010)。但是該假設缺少有效的約束,因為太古宙殼幔作用的地球物理記錄基本上被后期長期的改造所抹掉,所以要限定早期地殼的熱狀態(tài)只有地質(zhì)學、巖石學和地球化學方法是目前基本可用的手段,因此我們希望通過這些可用的方法,初步判斷地球早期的熱狀態(tài)。目前的一個關鍵科學問題是太古宙不同時期地殼究竟有多厚,地熱梯度究竟有多高,地幔溫度究竟比現(xiàn)代高多少?到目前為止這些還是未知數(shù),不同的作者基于不同的模型有各種各樣的假設或者推斷?;谀壳暗倪@種情況,我們研究組發(fā)展了利用純殼源鐵鎂質(zhì)巖石部分熔融產(chǎn)生的TTG巖漿實驗、微量元素理論模擬和一維熱傳導,結(jié)合長英質(zhì)地殼放射性熱源,初步獲得了華北克拉通中東部陸塊中-新太古代時期的地殼厚度和莫霍面地熱梯度(具體方法、步驟和成果將另文發(fā)表)。研究結(jié)果表明~2.5Ga華北克拉通中東部各區(qū)地殼厚度為35~56km,其中中條山44~47km、魯西地區(qū)40~43km、膠東地區(qū)35~38km、遼北-吉南40~43km、登封-太華地區(qū)50~53km和冀東-遼西地區(qū)53~56km(Sunetal., 2019b, 2020)?!?.5Ga各地區(qū)的莫霍面地熱梯度分別為中條山13.2℃/km、魯西地區(qū)13.2℃/km、膠東地區(qū)18℃/km、遼北-吉南15.1℃/km、登封-太華10.1℃/km和冀東-遼西8.7℃/km(Sunetal., 2019b, 2020)。上述估算結(jié)果表明華北克拉通中東部新太古代晚期地殼厚度和地熱梯度具有明顯的空間分帶性,南緣的登封-太華地區(qū)和北部的冀東-遼西地區(qū)地殼厚度大于50km,而地熱梯度低(8.7~10.1℃/km),與現(xiàn)代地熱梯度相當或略有不同(按照均變估計約11.7℃/km,比現(xiàn)代地幔溫度高不足100K),完全滿足現(xiàn)代板塊構(gòu)造起主導作用的物理條件。而中條、魯西地區(qū)略低于13.5℃/km,相當于比現(xiàn)代地幔溫度高約100~150K(低于175K),板塊構(gòu)造體制仍可以起主導作用,但是應表現(xiàn)為熱俯沖,即規(guī)模小、頻次高、俯沖角度大,導致俯沖板片回轉(zhuǎn)和斷離深度明顯低于現(xiàn)代板片俯沖。而遼北-吉南地區(qū)和膠東地區(qū)明顯較高的地熱梯度(>14.2℃/km),可能反映了弧后盆地的熱狀態(tài)特征(圖6; Gerya, 2014; Sizovaetal., 2010)。圖6反映了這樣的熱狀態(tài)均發(fā)育在高壓麻粒巖和高壓-中壓麻粒巖相分界線上。空間上的這種熱狀態(tài)與現(xiàn)代板塊構(gòu)造從NW向SE向的俯沖和弧后盆地熱狀態(tài)特征類似(Baietal., 2016; Fuetal., 2018, 2019; Guoetal., 2013, 2015b, 2017b; Wangetal., 2013, 2015; 劉樹文等, 2018b)。這種~2.5Ga時期的熱狀態(tài)說明,新太古代晚期華北克拉通中東部滿足類現(xiàn)代板塊構(gòu)造和熱俯沖的熱狀態(tài)條件,不能排除東部地區(qū)吉南-遼東-膠東等地區(qū)局部存在地幔柱(Plume)或者沉降(Sagduction or Dripping)構(gòu)造體制的可能性。

    結(jié)合上述的華北克拉通中東部新太古代晚期的熱狀態(tài)和地殼厚度研究結(jié)果,區(qū)域上遼北、遼南、冀東和魯西地區(qū)與這些變質(zhì)火山巖伴生有大量的TTG片麻巖和贊岐巖質(zhì)侵入巖,高鎂安山巖廣泛存在,典型的玻安巖在部分地區(qū)存在,明顯與俯沖相關的流體、熔體的地幔交代作用有關(Gaoetal., 2019, 2020a, b; Guoetal., 2017c; Sunetal., 2019a, 2020; Wangetal., 2015),說明新太古代晚期地球的殼幔動力學已經(jīng)由早期地幔柱體制(mantle plume regime)為主轉(zhuǎn)化為以早期板塊俯沖動力學體制為主,地幔柱高熱未虧損深部軟流圈地幔來源的巖漿產(chǎn)物已經(jīng)很少,說明這種動力學作用在局部還存在,但是已經(jīng)不是主導的殼幔動力學體制。因此,華北克拉通在新太古代晚期總體表現(xiàn)為以早期板塊體制為主導,局部為地幔柱和板塊構(gòu)造聯(lián)合作用的殼幔動力學體制(Gaoetal., 2019)。

    5 結(jié)語

    結(jié)合上述地殼厚度和熱狀態(tài)研究,新太古代晚期華北克拉通中東部的地熱梯度發(fā)生了明顯分化,冀東-遼西等地區(qū)表現(xiàn)為厚的地殼厚度和低的地熱梯度,遼北-吉南地區(qū)和膠東地區(qū)表現(xiàn)為相對較薄的地殼厚度和較高的地熱梯度,局部超出了前板片俯沖的地熱梯度范圍,而其它地區(qū)在現(xiàn)代地熱梯度和熱俯沖的地幔溫度范圍之間。結(jié)合新太古代晚期變質(zhì)火山巖組合和成因研究,我們認為在新太古代晚期,板塊構(gòu)造體制已經(jīng)是最主要的殼幔動力學體制,地幔柱構(gòu)造體制、板塊構(gòu)造與地幔柱聯(lián)合作用體制仍然在局部地區(qū)存在,但是其作用的范圍和強度已經(jīng)明顯減小。因此中太古代到新太古代晚期,隨著地幔潛能溫度的下降,地幔柱和板片俯沖的動力學體制可能是相互關聯(lián)的此消彼長的動力學過程,而不是一個災變事件導致的突變過程。

    謹以此文祝賀沈其韓院士百年華誕,祝愿沈先生健康長壽!

    猜你喜歡
    鎂質(zhì)安山巖火山巖
    制備工藝對鎂質(zhì)膨脹劑水化活性和微觀特性影響研究
    碳化作用對鎂質(zhì)膠凝材料微觀結(jié)構(gòu)演變過程的影響
    接財接福
    寶藏(2021年4期)2021-05-27 08:10:50
    含廢鎂碳磚中間包鎂質(zhì)干式料的研制
    昆鋼科技(2021年6期)2021-03-09 06:10:22
    渤中34-9油田古近系火山巖巖相特征與分布預測
    西藏拿若銅礦床安山巖元素地球化學特征研究?
    王官屯油田安山巖油藏滾動評價探索
    激發(fā)極化法尋找固體礦產(chǎn)的應用
    二連油田安山巖裂縫性儲層壓裂技術研究
    牛東火山巖油藏壓裂技術研究與應用
    一级作爱视频免费观看| 精品国产乱子伦一区二区三区| 免费观看精品视频网站| 亚洲色图av天堂| 欧美+亚洲+日韩+国产| 国产精品一及| 成人手机av| 国产精品久久电影中文字幕| 色老头精品视频在线观看| 亚洲乱码一区二区免费版| 亚洲精品粉嫩美女一区| 日本一区二区免费在线视频| 九色国产91popny在线| 成年版毛片免费区| 婷婷六月久久综合丁香| 这个男人来自地球电影免费观看| 日韩欧美国产在线观看| 久久久国产欧美日韩av| 丝袜美腿诱惑在线| 国产亚洲欧美98| 精品第一国产精品| 日本成人三级电影网站| 黄片小视频在线播放| a级毛片在线看网站| 99久久综合精品五月天人人| 国产亚洲精品一区二区www| 亚洲全国av大片| 十八禁网站免费在线| 中国美女看黄片| 国产激情久久老熟女| 99国产精品一区二区三区| 亚洲精品一区av在线观看| 亚洲av日韩精品久久久久久密| 99在线视频只有这里精品首页| 国产在线精品亚洲第一网站| 成人国产综合亚洲| 美女午夜性视频免费| 久久久国产精品麻豆| 久99久视频精品免费| 啦啦啦观看免费观看视频高清| 日韩中文字幕欧美一区二区| 一二三四社区在线视频社区8| 日韩大码丰满熟妇| 怎么达到女性高潮| 人人妻,人人澡人人爽秒播| 国产精品日韩av在线免费观看| 欧美成狂野欧美在线观看| 精品第一国产精品| 久久久久久人人人人人| 一个人免费在线观看的高清视频| 成年女人毛片免费观看观看9| 国产熟女xx| 国产av在哪里看| 90打野战视频偷拍视频| 99热6这里只有精品| 天堂av国产一区二区熟女人妻 | 美女高潮喷水抽搐中文字幕| 欧美色视频一区免费| 99热这里只有是精品50| 午夜视频精品福利| 久久久久久久午夜电影| 99久久精品国产亚洲精品| 99久久99久久久精品蜜桃| 久久精品人妻少妇| 制服丝袜大香蕉在线| 性欧美人与动物交配| 熟妇人妻久久中文字幕3abv| 久久婷婷人人爽人人干人人爱| 一级片免费观看大全| 国产成人啪精品午夜网站| 欧美日韩精品网址| 久久中文看片网| 高清毛片免费观看视频网站| 99精品在免费线老司机午夜| 国产爱豆传媒在线观看 | 欧美日韩国产亚洲二区| 亚洲自拍偷在线| 在线观看美女被高潮喷水网站 | 男人舔女人的私密视频| 国产精品爽爽va在线观看网站| 美女 人体艺术 gogo| 中文字幕高清在线视频| 一级a爱片免费观看的视频| 国产成人精品无人区| 亚洲精品国产一区二区精华液| 老熟妇乱子伦视频在线观看| 1024香蕉在线观看| 久久九九热精品免费| 91老司机精品| 人人妻,人人澡人人爽秒播| 欧美日韩亚洲综合一区二区三区_| 最新在线观看一区二区三区| 亚洲色图av天堂| 国产精品久久久人人做人人爽| 18美女黄网站色大片免费观看| 亚洲欧洲精品一区二区精品久久久| 熟女电影av网| 亚洲中文av在线| bbb黄色大片| 黑人操中国人逼视频| 国产男靠女视频免费网站| 在线观看美女被高潮喷水网站 | 人人妻人人澡欧美一区二区| 91av网站免费观看| 国产高清videossex| 欧美成人午夜精品| 波多野结衣高清作品| 在线观看日韩欧美| 一本大道久久a久久精品| www日本黄色视频网| 99久久无色码亚洲精品果冻| 精品第一国产精品| 午夜免费成人在线视频| 一级片免费观看大全| 观看免费一级毛片| 亚洲真实伦在线观看| 男女之事视频高清在线观看| 最新美女视频免费是黄的| 身体一侧抽搐| 免费观看人在逋| 成熟少妇高潮喷水视频| 国产伦在线观看视频一区| 亚洲一区二区三区色噜噜| 夜夜看夜夜爽夜夜摸| 日韩免费av在线播放| 日韩成人在线观看一区二区三区| 日韩欧美精品v在线| 草草在线视频免费看| 亚洲熟妇熟女久久| 男女视频在线观看网站免费 | 久久婷婷成人综合色麻豆| 午夜福利欧美成人| or卡值多少钱| 制服诱惑二区| 精品久久蜜臀av无| 国产成人精品久久二区二区免费| 国产精品1区2区在线观看.| 精品人妻1区二区| 国产精品一区二区三区四区久久| 男女视频在线观看网站免费 | 人人妻人人澡欧美一区二区| 一本综合久久免费| 母亲3免费完整高清在线观看| 在线观看午夜福利视频| 国产亚洲av嫩草精品影院| x7x7x7水蜜桃| 亚洲欧洲精品一区二区精品久久久| 韩国av一区二区三区四区| 国产99白浆流出| 久久精品夜夜夜夜夜久久蜜豆 | 午夜免费激情av| 亚洲av电影不卡..在线观看| netflix在线观看网站| 日本黄色视频三级网站网址| 国产不卡一卡二| 99国产综合亚洲精品| 别揉我奶头~嗯~啊~动态视频| 亚洲av电影不卡..在线观看| 国产伦一二天堂av在线观看| 成人国产一区最新在线观看| 久久香蕉国产精品| 91麻豆av在线| 少妇粗大呻吟视频| 免费搜索国产男女视频| 看黄色毛片网站| 男女下面进入的视频免费午夜| 欧美乱妇无乱码| 国产亚洲精品综合一区在线观看 | 脱女人内裤的视频| 欧美在线黄色| 久久香蕉激情| 日韩精品中文字幕看吧| 99在线人妻在线中文字幕| 国产区一区二久久| 国产精品一区二区精品视频观看| 国产精品爽爽va在线观看网站| 成人三级做爰电影| 精品国产超薄肉色丝袜足j| videosex国产| 久久热在线av| 在线观看一区二区三区| 欧美成狂野欧美在线观看| 国产亚洲精品一区二区www| 人人妻人人澡欧美一区二区| av在线天堂中文字幕| 国产激情久久老熟女| 悠悠久久av| 男女之事视频高清在线观看| 国产激情偷乱视频一区二区| 国产一区二区在线观看日韩 | АⅤ资源中文在线天堂| 午夜福利成人在线免费观看| 亚洲成人久久性| 欧美日韩乱码在线| xxx96com| 日本 av在线| 午夜精品久久久久久毛片777| 国产成+人综合+亚洲专区| 久久精品国产清高在天天线| 久久久久国产一级毛片高清牌| 精品国产超薄肉色丝袜足j| 免费看美女性在线毛片视频| 91老司机精品| 国产高清激情床上av| 国产真实乱freesex| 日韩国内少妇激情av| 91麻豆精品激情在线观看国产| x7x7x7水蜜桃| 亚洲精品中文字幕在线视频| 国产亚洲精品久久久久5区| 观看免费一级毛片| 2021天堂中文幕一二区在线观| 在线观看午夜福利视频| 少妇熟女aⅴ在线视频| 欧美黑人欧美精品刺激| www.精华液| 日韩免费av在线播放| 欧美黑人巨大hd| 黄色毛片三级朝国网站| 免费电影在线观看免费观看| 中出人妻视频一区二区| 不卡一级毛片| www.www免费av| 中文字幕最新亚洲高清| 国产亚洲精品一区二区www| 日韩欧美国产在线观看| 日本成人三级电影网站| 亚洲精品美女久久久久99蜜臀| 桃红色精品国产亚洲av| 精品午夜福利视频在线观看一区| 亚洲自偷自拍图片 自拍| 亚洲国产精品999在线| 久久中文字幕一级| 免费在线观看亚洲国产| 91在线观看av| 欧美日本亚洲视频在线播放| 99国产综合亚洲精品| 欧美3d第一页| 日韩国内少妇激情av| 国产真人三级小视频在线观看| 免费人成视频x8x8入口观看| 精品午夜福利视频在线观看一区| 欧美不卡视频在线免费观看 | 亚洲成人久久爱视频| 9191精品国产免费久久| 亚洲自拍偷在线| 男女视频在线观看网站免费 | 色精品久久人妻99蜜桃| 国产精品乱码一区二三区的特点| 男女床上黄色一级片免费看| 国产高清videossex| 亚洲乱码一区二区免费版| 香蕉国产在线看| 岛国在线观看网站| 午夜亚洲福利在线播放| 色综合婷婷激情| 精品乱码久久久久久99久播| 亚洲第一电影网av| 一区福利在线观看| 色综合欧美亚洲国产小说| 亚洲专区中文字幕在线| 最近在线观看免费完整版| 成年版毛片免费区| 欧美成人性av电影在线观看| 国产亚洲精品av在线| 亚洲av美国av| 亚洲第一电影网av| 亚洲精品久久成人aⅴ小说| 欧美色欧美亚洲另类二区| 久久精品91蜜桃| 两个人视频免费观看高清| 国产激情偷乱视频一区二区| 亚洲精品美女久久av网站| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕人成人乱码亚洲影| 久久草成人影院| 欧美+亚洲+日韩+国产| 亚洲成av人片免费观看| 成人特级黄色片久久久久久久| 成人午夜高清在线视频| 欧美性长视频在线观看| 我要搜黄色片| 香蕉久久夜色| 国产精品自产拍在线观看55亚洲| 成人欧美大片| 床上黄色一级片| 正在播放国产对白刺激| 两个人免费观看高清视频| 国产亚洲欧美98| 美女高潮喷水抽搐中文字幕| 国产精品永久免费网站| 久久亚洲精品不卡| 欧美另类亚洲清纯唯美| 午夜福利欧美成人| 成人三级做爰电影| 国产一区二区在线观看日韩 | av有码第一页| 日韩大尺度精品在线看网址| 中文字幕人妻丝袜一区二区| 国产在线精品亚洲第一网站| 亚洲av美国av| 精品久久久久久久久久免费视频| 欧美午夜高清在线| 日韩欧美在线二视频| 村上凉子中文字幕在线| 亚洲五月婷婷丁香| 中文字幕最新亚洲高清| 国产v大片淫在线免费观看| 国产伦在线观看视频一区| 亚洲自拍偷在线| 国产精品久久久久久人妻精品电影| 日本在线视频免费播放| 亚洲精品一卡2卡三卡4卡5卡| 国产高清视频在线观看网站| 三级毛片av免费| 岛国在线观看网站| 久久精品国产清高在天天线| 国产成人aa在线观看| 成人高潮视频无遮挡免费网站| 国产欧美日韩一区二区三| 香蕉国产在线看| 中文资源天堂在线| 国内揄拍国产精品人妻在线| 欧美色视频一区免费| 亚洲专区字幕在线| 免费电影在线观看免费观看| 国内精品久久久久久久电影| 一二三四在线观看免费中文在| 99在线人妻在线中文字幕| 老司机午夜福利在线观看视频| 999久久久精品免费观看国产| 亚洲性夜色夜夜综合| 欧美国产日韩亚洲一区| 亚洲精品中文字幕在线视频| 国产精品自产拍在线观看55亚洲| 伦理电影免费视频| 人人妻,人人澡人人爽秒播| 岛国在线免费视频观看| 亚洲第一电影网av| 日本黄色视频三级网站网址| 麻豆国产97在线/欧美 | 精品一区二区三区av网在线观看| 婷婷精品国产亚洲av在线| 老鸭窝网址在线观看| 日日爽夜夜爽网站| 一a级毛片在线观看| 久9热在线精品视频| 国内精品一区二区在线观看| 88av欧美| 欧美中文日本在线观看视频| 亚洲电影在线观看av| 国内毛片毛片毛片毛片毛片| 国产aⅴ精品一区二区三区波| av天堂在线播放| 中国美女看黄片| 日韩欧美在线乱码| 午夜免费观看网址| 国产亚洲欧美98| 国产99白浆流出| 精品国产乱码久久久久久男人| 一卡2卡三卡四卡精品乱码亚洲| 久久国产精品影院| 天堂√8在线中文| 国产97色在线日韩免费| 亚洲色图av天堂| 伦理电影免费视频| 欧美极品一区二区三区四区| 99精品在免费线老司机午夜| 少妇人妻一区二区三区视频| 一本久久中文字幕| xxx96com| 美女免费视频网站| 最近在线观看免费完整版| 色老头精品视频在线观看| 国产午夜精品久久久久久| 国产久久久一区二区三区| 成人国产一区最新在线观看| 免费看十八禁软件| 日韩欧美在线乱码| 50天的宝宝边吃奶边哭怎么回事| 色老头精品视频在线观看| 日本撒尿小便嘘嘘汇集6| 热99re8久久精品国产| 韩国av一区二区三区四区| 久久久久久亚洲精品国产蜜桃av| 99在线视频只有这里精品首页| 又大又爽又粗| av福利片在线观看| 午夜亚洲福利在线播放| 欧美最黄视频在线播放免费| 男女做爰动态图高潮gif福利片| 97超级碰碰碰精品色视频在线观看| 国产精品久久久久久久电影 | 变态另类丝袜制服| 成人三级做爰电影| 国产不卡一卡二| e午夜精品久久久久久久| 国产av在哪里看| 成人av在线播放网站| 精品久久久久久成人av| 少妇的丰满在线观看| 国产精品一及| 久久精品影院6| 精品免费久久久久久久清纯| 很黄的视频免费| 精品高清国产在线一区| 神马国产精品三级电影在线观看 | 两性夫妻黄色片| 可以在线观看的亚洲视频| 午夜精品一区二区三区免费看| 成人高潮视频无遮挡免费网站| 777久久人妻少妇嫩草av网站| 午夜日韩欧美国产| 亚洲va日本ⅴa欧美va伊人久久| 欧美精品亚洲一区二区| 18禁观看日本| 久久精品国产清高在天天线| 日本免费a在线| 精品熟女少妇八av免费久了| 丰满人妻熟妇乱又伦精品不卡| 国产精品野战在线观看| 久热爱精品视频在线9| 日韩精品青青久久久久久| 亚洲激情在线av| 琪琪午夜伦伦电影理论片6080| 国产成人一区二区三区免费视频网站| 国产爱豆传媒在线观看 | 男人舔女人的私密视频| 中亚洲国语对白在线视频| 在线国产一区二区在线| 一个人免费在线观看的高清视频| 99国产精品99久久久久| 久久久久久大精品| 中文字幕高清在线视频| 久久精品国产清高在天天线| 91麻豆av在线| 亚洲片人在线观看| 老司机福利观看| 国产精品综合久久久久久久免费| 午夜福利视频1000在线观看| 999精品在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 成人av一区二区三区在线看| 国产成人av教育| 18禁观看日本| 国产熟女午夜一区二区三区| 免费看日本二区| 国产高清videossex| 中文字幕最新亚洲高清| 三级男女做爰猛烈吃奶摸视频| 亚洲激情在线av| 欧美色视频一区免费| 国产高清videossex| 成人午夜高清在线视频| 亚洲成a人片在线一区二区| 国产精品亚洲美女久久久| 无遮挡黄片免费观看| 在线观看美女被高潮喷水网站 | 久久久国产精品麻豆| 成年女人毛片免费观看观看9| 国产亚洲av高清不卡| 国产蜜桃级精品一区二区三区| cao死你这个sao货| 国产精品亚洲av一区麻豆| 欧美久久黑人一区二区| 欧美日韩国产亚洲二区| 国产免费男女视频| 人人妻人人澡欧美一区二区| 天堂动漫精品| 男人的好看免费观看在线视频 | 99国产极品粉嫩在线观看| 久久精品国产清高在天天线| 一个人观看的视频www高清免费观看 | 亚洲成人久久爱视频| 50天的宝宝边吃奶边哭怎么回事| 日韩精品青青久久久久久| 欧美色欧美亚洲另类二区| 色在线成人网| 99热这里只有精品一区 | 91成年电影在线观看| 欧美乱码精品一区二区三区| 免费搜索国产男女视频| 亚洲一区二区三区色噜噜| 日本一区二区免费在线视频| 美女免费视频网站| 亚洲欧洲精品一区二区精品久久久| 久久久久久久久久黄片| 真人做人爱边吃奶动态| 91av网站免费观看| 精品福利观看| 日本一二三区视频观看| 18禁观看日本| 久久久久久亚洲精品国产蜜桃av| 久久精品aⅴ一区二区三区四区| 国产aⅴ精品一区二区三区波| 天堂动漫精品| 国产成+人综合+亚洲专区| 久久久久久久精品吃奶| 亚洲精品中文字幕在线视频| 黄色视频,在线免费观看| 久久精品国产99精品国产亚洲性色| 欧美精品啪啪一区二区三区| 一本大道久久a久久精品| avwww免费| 成年人黄色毛片网站| 免费一级毛片在线播放高清视频| 免费看a级黄色片| 老鸭窝网址在线观看| 怎么达到女性高潮| av中文乱码字幕在线| 99国产精品一区二区蜜桃av| 久久午夜综合久久蜜桃| 大型黄色视频在线免费观看| 欧美性猛交╳xxx乱大交人| 欧美丝袜亚洲另类 | 精品久久久久久成人av| 欧美日本视频| 亚洲精品在线观看二区| 久久中文字幕人妻熟女| 最新在线观看一区二区三区| 99国产精品一区二区三区| 亚洲国产看品久久| 18禁裸乳无遮挡免费网站照片| 精品久久久久久久久久免费视频| 亚洲 欧美 日韩 在线 免费| 99精品欧美一区二区三区四区| 亚洲精品国产一区二区精华液| 精品乱码久久久久久99久播| 久久久久久久久中文| 日日干狠狠操夜夜爽| 久久久久久免费高清国产稀缺| 美女 人体艺术 gogo| 美女午夜性视频免费| 高潮久久久久久久久久久不卡| 久久久精品欧美日韩精品| 国产精品国产高清国产av| 好男人在线观看高清免费视频| 亚洲真实伦在线观看| 正在播放国产对白刺激| 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧美日韩东京热| 男人的好看免费观看在线视频 | 丰满人妻熟妇乱又伦精品不卡| 女生性感内裤真人,穿戴方法视频| 久久精品亚洲精品国产色婷小说| 国产成人影院久久av| 成人国产一区最新在线观看| 亚洲av电影不卡..在线观看| 床上黄色一级片| 国产精品国产高清国产av| 久久香蕉激情| 亚洲成人中文字幕在线播放| 人妻丰满熟妇av一区二区三区| 久久中文字幕一级| 国产视频一区二区在线看| 免费在线观看亚洲国产| 午夜免费观看网址| 中文资源天堂在线| 久久精品成人免费网站| 欧美+亚洲+日韩+国产| svipshipincom国产片| 麻豆国产97在线/欧美 | 免费看美女性在线毛片视频| 女同久久另类99精品国产91| 午夜激情福利司机影院| 亚洲美女视频黄频| 麻豆一二三区av精品| 欧美在线一区亚洲| 精品电影一区二区在线| 熟女少妇亚洲综合色aaa.| 十八禁网站免费在线| 搞女人的毛片| 午夜精品在线福利| 国产一区二区在线观看日韩 | 无人区码免费观看不卡| 国产精品香港三级国产av潘金莲| 国内精品久久久久精免费| 黄频高清免费视频| 亚洲中文字幕一区二区三区有码在线看 | 一夜夜www| 欧美日韩黄片免| 又爽又黄无遮挡网站| 99精品欧美一区二区三区四区| 国产精品乱码一区二三区的特点| 真人一进一出gif抽搐免费| 国产在线观看jvid| 国产av麻豆久久久久久久| 免费人成视频x8x8入口观看| 欧美一区二区精品小视频在线| 色av中文字幕| 国产午夜精品久久久久久| 国产精品一区二区三区四区免费观看 | 国产三级黄色录像| 国产精品98久久久久久宅男小说| 美女 人体艺术 gogo| 免费在线观看完整版高清| 日韩国内少妇激情av| 一进一出抽搐动态| 美女黄网站色视频| 久久久久国产精品人妻aⅴ院| 操出白浆在线播放| 国产精品久久久久久精品电影| 欧美日本视频| 久9热在线精品视频| 国产精品久久久人人做人人爽| 91九色精品人成在线观看| 欧美乱妇无乱码| 亚洲国产看品久久| 91在线观看av| 国产熟女午夜一区二区三区| 他把我摸到了高潮在线观看| 欧美性长视频在线观看| 午夜激情福利司机影院| 男女床上黄色一级片免费看| 欧美日韩瑟瑟在线播放|