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

    源區(qū)組成與轉(zhuǎn)熔礦物組合的選擇性帶入對(duì)花崗巖成分變化的影響:以南秦嶺寧陜花崗巖體群為例*

    2021-12-29 03:03:22王日香李小偉管琪孫雨沁李睿哲謝沛伶吳斌斌夏文月
    巖石學(xué)報(bào) 2021年12期
    關(guān)鍵詞:獨(dú)居石源區(qū)老城

    王日香 李小偉,** 管琪 孫雨沁 李睿哲 謝沛伶 吳斌斌 夏文月

    1. 地質(zhì)過(guò)程與礦產(chǎn)資源國(guó)家重點(diǎn)實(shí)驗(yàn)室,中國(guó)地質(zhì)大學(xué)地球科學(xué)與資源學(xué)院,北京 1000832. 自然資源部金礦成礦過(guò)程與資源利用重點(diǎn)實(shí)驗(yàn)室, 山東省金屬礦產(chǎn)成礦地質(zhì)過(guò)程與資源利用重點(diǎn)實(shí)驗(yàn)室, 山東省地質(zhì)科學(xué)研究院, 濟(jì)南2500133. 河北地質(zhì)大學(xué)地球科學(xué)學(xué)院,河北省戰(zhàn)略性關(guān)鍵礦產(chǎn)資源重點(diǎn)實(shí)驗(yàn)室,石家莊 050031

    花崗巖作為地球大陸地殼的主要組成部分,記錄了地殼生長(zhǎng)和演化的重要信息(Hawkesworth and Kemp, 2006;Kempetal., 2007;Moyenetal., 2017)。因此,花崗巖源區(qū)性質(zhì)以及對(duì)其巖漿演化過(guò)程的研究,可以為地殼生長(zhǎng)、物質(zhì)再循環(huán)和巖漿分異提供重要的約束(Clemensetal., 2011;Clemens and Stevens,2012;Bailieetal., 2020)。普遍認(rèn)為,花崗巖類(lèi)地球化學(xué)特征的多樣性受控于多種因素,如源區(qū)組成的不同,熔融條件和過(guò)程的差異,同化混染及分離結(jié)晶作用(Patio Douce, 1999;Clemens and Stevens, 2012),不同批次地殼或地幔來(lái)源熔體的混合(Kempetal., 2007;Peruginietal., 2008;Applebyetal., 2010; Wangetal., 2013a)等。其中,不同的源區(qū)組成包含變泥質(zhì)巖、變?cè)颇钙瑤r、變硬砂巖、榴輝巖和變玄武巖(變輝長(zhǎng)巖、角閃巖)、變英云閃長(zhǎng)巖、變安山巖、變英安巖等 (Sissonetal., 2005;Gaoetal., 2016;Zhuetal., 2019a);熔融條件的差異包括不同的溫度、壓力和水含量(Clemens, 2003, 2018;Clemens and Stevens, 2012)等;而花崗巖熔融過(guò)程的差異則包括平衡熔融與非平衡熔融過(guò)程(Tangetal., 2014;Wangetal., 2018;Zhuetal., 2019b)。 源區(qū)轉(zhuǎn)熔礦物的選擇性帶入(peritectic assemblage entrainment,PAE,Stevensetal., 2007)被認(rèn)為是影響花崗巖成分變化的重要控制因素(Stevensetal., 2007;Clemensetal., 2011;Huetal., 2019;Bailieetal., 2020;Zhuetal., 2020, 2021;王洛娟等,2021),如國(guó)外對(duì)南非埃迪卡拉紀(jì)中期至早寒武世開(kāi)普花崗巖套的研究(Cape Granite Suite,CGS;Stevensetal., 2007;Clemensetal., 2011;Clemens and Stevens, 2012)、國(guó)內(nèi)對(duì)揚(yáng)子板塊西部中元古代晚期過(guò)鋁質(zhì)花崗巖(Zhuetal., 2020)和新元古代高鎂鐵質(zhì)I型花崗巖(Zhuetal., 2021)的研究等。這些研究認(rèn)為,PAE模型是造成花崗巖高的鎂鐵指數(shù)(摩爾 Fe+Mg)以及Ti、Ca、Zr、Th、Hf等元素含量變化的重要影響因素。

    南秦嶺構(gòu)造帶中段發(fā)育有大規(guī)模的早中生代花崗質(zhì)巖體群,包括光頭山巖體群、寧陜巖體群和五龍巖體群等(圖1),它們是探究花崗巖源區(qū)組成和PAE模型影響酸性侵入巖成分變化的良好對(duì)象。前人對(duì)寧陜巖體的巖石成因進(jìn)行了很多研究(張宏飛等,1997;駱金誠(chéng)等,2010;Dengetal., 2016;Luetal., 2016, 2017),但仍存在不同的認(rèn)識(shí),歸納起來(lái),主要有:(1)寧陜巖體是碰撞背景下南秦嶺構(gòu)造帶的深部地殼物質(zhì)部分熔融的產(chǎn)物(張宏飛等,1997);(2)寧陜巖體早期巖漿(222~216Ma)由元古代變玄武巖部分熔融的巖漿與虧損地幔巖漿混合形成,而晚期巖漿(~210Ma)形成于沉積巖源區(qū)的部分熔融,該巖體形成于同碰撞至后碰撞的過(guò)渡階段(Yangetal., 2012);(3)寧陜巖體中胭脂壩花崗巖是由具有揚(yáng)子板塊屬性的新元古代地殼物質(zhì)重熔形成的(韋龍猛等,2016;方博文等,2017)??梢钥闯?,目前對(duì)影響寧陜巖體巖漿成分變化的主控因素和成因機(jī)制尚未有定論。如前文所述,除源區(qū)物質(zhì)組成的影響外,PAE模型是否影響巖漿成分的變化尚需評(píng)估。

    圖1 中國(guó)中部區(qū)域地質(zhì)簡(jiǎn)圖(a, 據(jù)Yang et al.,2012)和秦嶺造山帶構(gòu)造單元、沉積層序和早中生代花崗巖類(lèi)的分布(b, 據(jù)Dong et al., 2011修改)年齡數(shù)據(jù)據(jù)胡健民等,2004;張宗清等,2006;Yang et al., 2012;Deng et al., 2016;黃雄飛,2016;韋龍猛等,2016;方博文等,2017. WQ-西秦嶺;SG-松潘-甘孜地體;QD(QDM)-柴達(dá)木地體;QL-祁連山造山帶;NQ-北秦嶺;SQ-南秦嶺;YB-揚(yáng)子板塊;NCB-華北陸塊;QT-羌塘地體;LT-拉薩地體Fig.1 Geological map of main blocks and orogenic belts in central China (a, after Yang et al., 2012) and geological sketch of tectonic units, sedimentary sequences and distribution of Early Mesozoic granitoids in the Qinling orogenic belt (b, modified after Dong et al., 2011)Age data according to Hu et al., 2004; Zhang et al., 2006; Yang et al., 2012; Deng et al., 2016; Huang, 2016; Wei et al., 2016; Fang et al., 2017. WQ-West Qinling; SG-Songpan-Garzê Terrane;QD (QDM)-Qaidam; QL-Qilian Terrane; NQ-North Qinling; SQ-South Qinling; YB-Yangtze Block; NCB-North China Block; QT-Qiangtang Terrane; LT-Lhasa Terrane

    本文選擇寧陜花崗巖中代表性巖體為主要研究對(duì)象,對(duì)其開(kāi)展了詳細(xì)的野外地質(zhì)調(diào)查,通過(guò)鋯石U-Pb-Hf同位素和獨(dú)居石U-Pb同位素分析,以及全巖主、微量元素和Sr-Nd同位素地球化學(xué)分析,對(duì)其侵位時(shí)代、源區(qū)組成、巖漿演化過(guò)程如何影響巖漿成分變化進(jìn)行了較為深入的探討,以期查明上述巖體的起源與演化過(guò)程,進(jìn)一步豐富酸性侵入巖的成因機(jī)制。

    1 地質(zhì)背景

    秦嶺造山帶夾持于太平洋構(gòu)造域、特提斯構(gòu)造域和古亞洲洋構(gòu)造域三者之間,呈近東西向展布。其西側(cè)與祁連-昆侖造山帶相鄰,北側(cè)以靈寶-魯山-舞陽(yáng)斷裂帶為界與華北地塊相鄰,南側(cè)以勉略-巴山-襄樊?dāng)嗔褞檫吔缗c揚(yáng)子板塊相接壤,東接大別-蘇魯超高壓變質(zhì)帶(Dongetal., 2011)。秦嶺造山帶及其鄰區(qū)由北到南可劃分為:華北板塊、北秦嶺構(gòu)造帶、商南-丹鳳縫合帶、南秦嶺構(gòu)造帶、勉縣-略陽(yáng)縫合帶和揚(yáng)子板塊(圖1a;張國(guó)偉等,1995)六個(gè)構(gòu)造單元。秦嶺造山帶主要經(jīng)歷了新元古代至早中生代的多期構(gòu)造-巖漿熱事件, 奠定了現(xiàn)今該造山帶的主體構(gòu)造格架(張成立等,2008)。

    秦嶺造山帶大規(guī)模產(chǎn)出的酸性侵入巖具有多期次、多旋回的特點(diǎn)(王曉霞等,2011;Wangetal., 2013a),它們的巖漿作用時(shí)代主要集中在:新元古代(980~710Ma)、古生代(510~400Ma)、早中生代(250~185Ma)和晚中生代(160~100Ma)(Wangetal., 2013b)。早中生代花崗巖類(lèi)以巖體群的形式發(fā)育于南秦嶺構(gòu)造帶內(nèi)(圖1b),其中包括光頭山巖體群、寧陜巖體群和五龍巖體群等。

    本文研究區(qū)位于南秦嶺構(gòu)造帶南部的勉縣-略陽(yáng)縫合帶北部,該縫合帶是由勉略洋盆于古生代至中生代演化而成(Mattaueretal., 1985;eng?r, 1985;Kr?neretal., 1993; Meng and Zhang, 1999, 2000;張國(guó)偉等,2001)。具體而言,其形成包括2個(gè)階段:古生代中期,勉略洋盆逐漸打開(kāi),并演化出獨(dú)立的秦嶺微地塊;古生代晚期開(kāi)始,勉略古洋盆持續(xù)向北俯沖,至早中生代時(shí)期秦嶺微地塊與南部揚(yáng)子板塊發(fā)生碰撞而形成勉縣-略陽(yáng)縫合帶。區(qū)內(nèi)發(fā)育太古宇佛坪雜巖與陡嶺雜巖等結(jié)晶基底,沉積蓋層主要為震旦紀(jì)-早古生代巨厚的被動(dòng)陸緣沉積體系和晚古生代至中三疊紀(jì)秦嶺微板塊內(nèi)陸表海沉積體系(圖1b;張國(guó)偉等,2001;Dongetal., 2015)。研究區(qū)內(nèi)出露有大面積的酸性侵入巖類(lèi)組合,它們侵位于新元古代到古生代地層中,沿勉略縫合帶北側(cè)展布,共同構(gòu)成了一條東西向的印支期花崗巖帶,該花崗巖帶被認(rèn)為是揚(yáng)子板塊北部與秦嶺造山帶碰撞的產(chǎn)物(Sunetal., 2002)。

    2 巖相學(xué)特征

    寧陜花崗巖位于南秦嶺構(gòu)造帶的蜂腰部位,總體呈東西向展布的不規(guī)則橢圓形,出露面積約為1160km2(圖1b),該巖體侵入至寒武紀(jì)-泥盆紀(jì)沉積蓋層中。寧陜花崗巖主要由胭脂壩、懶板凳和老城巖體組成,其中胭脂壩巖體位于寧陜花崗巖的東南部,北部為懶板凳巖體,西部為老城巖體(圖2)。

    圖2 西秦嶺寧陜巖體群地質(zhì)簡(jiǎn)圖(據(jù)Yang et al., 2012;韋龍猛等,2016)Fig.2 Geological map of Ningshan rock group in West Qinling (after Yang et al., 2012;Wei et al., 2016)

    胭脂壩巖體(~530km2)的主體巖性為中粗?;◢忛W長(zhǎng)巖,主要礦物組成為石英(25%~30%)、斜長(zhǎng)石(40%~45%)和鉀長(zhǎng)石(20%~25%);次要礦物為黑云母(~5%)、白云母 (~1%)和角閃石(1%~5%);副礦物為獨(dú)居石、鋯石、磷灰石及其他不透明礦物等。石英呈他形粒狀,部分顆粒可見(jiàn)波狀消光;斜長(zhǎng)石呈自形-半自形板狀,發(fā)育韻律環(huán)帶;微斜長(zhǎng)石發(fā)育格子雙晶;黑云母呈自形片狀,發(fā)育一組極完全解理,局部發(fā)生蝕變(圖3a-c)。

    圖3 寧陜巖體花崗巖正交偏光顯微鏡下照片(a-c)胭脂壩花崗閃長(zhǎng)巖中具有環(huán)帶結(jié)構(gòu)的斜長(zhǎng)石斑晶、具有格子雙晶的微斜長(zhǎng)石;(d)懶板凳二長(zhǎng)花崗巖中具有格子雙晶的微斜長(zhǎng)石;(e、f)老城二長(zhǎng)花崗巖中具有卡斯巴雙晶的鉀長(zhǎng)石. Pl-斜長(zhǎng)石;Kfs-鉀長(zhǎng)石;Mc-微斜長(zhǎng)石;Amp-角閃石;Bt-黑云母;Ms-白云母;Qtz-石英Fig.3 Microphotographs under CPL showing the mineral assemblages of the granite in the Ningshan granitoids(a-c) plagioclase phenocrysts with zoning texture and microcline with cross-hatch twinning from Yanzhiba granodiorite; (d) microcline with cross-hatch twinning from Langbandeng monzogranite; (e, f) K-feldspar with Carlsbad twinnings from Laocheng monzogranite. Pl-plagioclase; Kfs-K-feldspar; Mc-microcline; Amp-amphibole; Bt-biotite; Ms-muscovite; Qtz-quartz

    懶板凳巖體(~80km2)主體巖性為含二云母二長(zhǎng)花崗巖,主要礦物組成包括石英(25%~30%)、斜長(zhǎng)石(25%~30%)和鉀長(zhǎng)石(30%~40%);次要礦物有黑云母(1%~5%)和白云母(1%~2%);副礦物見(jiàn)獨(dú)居石、鋯石等(圖3d)。

    老城巖體(~550km2)主要由石英閃長(zhǎng)巖、花崗閃長(zhǎng)巖和二長(zhǎng)花崗巖組成,其中二長(zhǎng)花崗巖主要由石英(22%~25%)、斜長(zhǎng)石(30%~35%)和鉀長(zhǎng)石(25%~35%)組成,鉀長(zhǎng)石發(fā)育卡斯巴雙晶;次要礦物為黑云母(~5%)和角閃石(~2%);副礦物與上述兩個(gè)巖體類(lèi)似(圖3e, f)。

    3 樣品及測(cè)試方法

    本次研究共選取23件樣品(胭脂壩巖體10件、懶板凳巖體9件、老城巖體4件)進(jìn)行了全巖主微量元素分析,并對(duì)其中的6件樣品進(jìn)行了全巖Sr-Nd同位素分析。此外,對(duì)3件花崗巖類(lèi)樣品(NS18-8、NS18-16和NS18-29)進(jìn)行了鋯石U-Pb及Lu-Hf同位素分析,并對(duì)7件花崗巖類(lèi)樣品(NS18-2、NS18-7、NS18-8、NS18-12、NS18-16、NS18-27和NS18-29)進(jìn)行了獨(dú)居石U-Pb同位素測(cè)定。

    巖石樣品中鋯石與獨(dú)居石顆粒分選采用了常規(guī)的重選法和磁選法,在雙目顯微鏡下手工挑選無(wú)色透明、沒(méi)有裂隙、不含包體的鋯石、獨(dú)居石顆粒并使用環(huán)氧樹(shù)脂制靶,拋光后拍攝陰極發(fā)光(CL)圖像及透射光與反射光圖像,用以觀察鋯石顆粒的外部和內(nèi)部結(jié)構(gòu)。鋯石顆粒與獨(dú)居石顆粒的分選在河北省廊坊市誠(chéng)信地質(zhì)服務(wù)有限公司完成,在北京中科礦研檢測(cè)技術(shù)有限公司完成鋯石與獨(dú)居石的制靶和陰極發(fā)光(CL)顯微照相,在中國(guó)地質(zhì)大學(xué)(北京)顯微礦物實(shí)驗(yàn)室完成了獨(dú)居石的透、反射光照相。

    3.1 鋯石U-Pb-Hf同位素分析

    樣品的鋯石U-Pb同位素分析由桂林理工大學(xué)廣西隱伏金屬礦產(chǎn)勘查重點(diǎn)實(shí)驗(yàn)室完成,分析所用儀器為New Wave Research 研發(fā)的213nm的激光剝蝕系統(tǒng),聯(lián)用Agilent 7500s型ICP-MS。激光束斑直徑為32μm,激光頻率為5Hz,以He為載氣。Plesovice 鋯石(337Ma)為年齡外標(biāo),鋯石樣品中的微量元素含量利用SRM610作為外標(biāo)、Si做內(nèi)標(biāo),使用單內(nèi)標(biāo)多外標(biāo)法進(jìn)行校正,分析流程和具體參數(shù)設(shè)置可參見(jiàn)Liuetal.(2010)。測(cè)試過(guò)程中每8個(gè)測(cè)試樣品前后,均需測(cè)定2個(gè)鋯石標(biāo)樣Plesovice,并在每30個(gè)樣品中間插入3個(gè)GJ-1標(biāo)準(zhǔn)樣品點(diǎn),以監(jiān)控儀器狀態(tài),以保證測(cè)試的精確度。數(shù)據(jù)處理采用ICPMSDataCal程序(Liuetal., 2010)完成,鋯石U-Pb諧和圖的繪制及加權(quán)平均年齡的計(jì)算采用Isoplot3.0程序完成。

    樣品NS18-8、NS18-16、NS18-29原位微區(qū)鋯石Lu-Hf同位素測(cè)試在山東省地質(zhì)科學(xué)研究院利用激光剝蝕多接收杯電感耦合等離子體質(zhì)譜(LA-MC-ICP-MS)完成。激光剝蝕系統(tǒng)為193nm準(zhǔn)分子激光Coherent GeoLas Pro,MC-ICP-MS為T(mén)hermo Fisher Neptune Plus。前人的研究表明,激光剝蝕過(guò)程中采用氦氣作為載氣,比起氬氣可以將信號(hào)靈敏度提高2倍(Huetal., 2008a);同時(shí),少量氮?dú)獾囊脒€可進(jìn)一步提升大部分元素的靈敏度(Huetal., 2008b)。為了提高Hf同位素的分析靈敏度,本實(shí)驗(yàn)采用氦氣作為載氣,并引入了適量的氮?dú)?8ml/min),同時(shí)使用了Neptune Plus配備的高靈敏度Jet采樣錐和X截取錐組合。另外,激光剝蝕系統(tǒng)配置了一個(gè)信號(hào)平滑裝置,確保即使在激光脈沖頻率低達(dá)1Hz的情況下,也能獲得光滑的分析信號(hào)(Huetal., 2012a)。分析測(cè)試采用單點(diǎn)剝蝕模式,激光的輸出能量密度為5~8J/cm2,束斑直徑為44μm,激光頻率為6~8Hz。詳細(xì)的儀器操作條件和分析方法可參照Huetal. (2012b)。

    采用LA-MC-ICP-MS準(zhǔn)確測(cè)試鋯石Hf同位素的難點(diǎn)在于扣除176Yb和176Lu對(duì)176Hf的同量異位素的干擾。在本次實(shí)驗(yàn)中,筆者實(shí)時(shí)獲取了鋯石樣品自身的βYb用于干擾校正。179Hf/177Hf=0.7325和173Yb/171Yb=1.132685 (Fisheretal., 2014)被用于計(jì)算Hf和Yb的質(zhì)量分餾系數(shù)βHf和βYb。179Hf/177Hf 和173Yb/171Yb的比值被用于計(jì)算Hf (βHf)和Yb (βYb)的質(zhì)量偏差。使用176Yb/173Yb=0.79639(Fisheretal., 2014)來(lái)扣除176Yb對(duì)176Hf的同量異位干擾;176Lu/175Lu=0.02656(Blichert-Toftetal., 1997)被用來(lái)扣除干擾程度相對(duì)較小的176Lu對(duì)176Hf的同量異位干擾。由于Yb和Lu具有相似的物理化學(xué)屬性,因此在本實(shí)驗(yàn)中采用Yb的質(zhì)量分餾系數(shù)βYb來(lái)校正Lu的質(zhì)量分餾行為。鋯石91500作為外標(biāo)用來(lái)校正儀器的時(shí)間漂移。利用軟件ICPMSDataCal (Liuetal., 2010)完成分析數(shù)據(jù)(包括對(duì)樣品和空白信號(hào)的選擇、同位素質(zhì)量分餾校正參數(shù)的選擇)的離線處理。

    3.2 獨(dú)居石U-Pb年代學(xué)分析

    獨(dú)居石微區(qū)原位U-Pb同位素分析在中國(guó)地質(zhì)調(diào)查局天津地質(zhì)礦產(chǎn)研究所同位素實(shí)驗(yàn)室完成,所用儀器為激光剝蝕多接收等離子體質(zhì)譜儀(LA-MC-ICP-MS)。采用的激光剝蝕系統(tǒng)為美國(guó)ESI公司生產(chǎn)的NEW WAVE 193nm FX ArF準(zhǔn)分子激光器。激光剝蝕的束斑直徑為20μm,能量密度為10~11J/cm2,頻率為5Hz,所用載氣為氦氣。測(cè)試過(guò)程中每5個(gè)測(cè)試樣品前后均需測(cè)定2個(gè)獨(dú)居石標(biāo)準(zhǔn)樣品。處理數(shù)據(jù)時(shí),用獨(dú)居石標(biāo)樣44069來(lái)校正U-Pb同位素分餾。由TIMS給出的獨(dú)居石標(biāo)準(zhǔn)樣品44069的206Pb/238U年齡為424.9±0.4Ma(萬(wàn)渝生等,2004)。根據(jù)不同時(shí)間和不同的儀器狀態(tài)下用LA-MC-ICP-MS所測(cè)得的該獨(dú)居石標(biāo)樣的206Pb/238U年齡加權(quán)平均值與用TIMS所測(cè)得的206Pb/238U年齡值的差距來(lái)計(jì)算U-Pb同位素分餾系數(shù),其計(jì)算公式為:K=t1/t2。式中:K為用 LA-MC-ICP-MS測(cè)定時(shí)的U-Pb同位素分餾系數(shù),t1為由LA-MC-ICP-MS所測(cè)得的該獨(dú)居石標(biāo)樣的206Pb/238U 年齡加權(quán)平均值;t2為用TIMS所測(cè)得的該獨(dú)居石標(biāo)樣的206Pb/238U年齡值。研究對(duì)象樣品的206Pb/238U年齡計(jì)算公式為:T=t1/K。式中:T為研究對(duì)象樣品的206Pb/238U年齡值;t1為用LA-MC-ICP-MS測(cè)得的該研究對(duì)象樣品的206Pb/238U年齡加權(quán)平均值;K為用LA-MC-ICP-MS測(cè)定時(shí)的U-Pb同位素分餾系數(shù)。本文報(bào)道的獨(dú)居石樣品采用207Pb校正法和等時(shí)線法對(duì)普通鉛進(jìn)行校正(Andersen,2002)。

    3.3 全巖元素地球化學(xué)分析

    本研究共選取了23件新鮮樣品用于全巖元素地球化學(xué)分析。全巖主量元素分析在中國(guó)地質(zhì)大學(xué)(北京)科學(xué)研究院實(shí)驗(yàn)中心完成,分析儀器為島津X射線熒光光譜儀(XRF-1800),具體分析步驟參考Liuetal. (2021)。在本次分析中,每10件樣品間穿插了一件平衡樣,用以監(jiān)測(cè)數(shù)據(jù)質(zhì)量,分析精度優(yōu)于±5%。

    全巖微量元素含量分析由武漢上譜分析科技有限責(zé)任公司完成,分析儀器為Agilent 7700e ICP-MS,具體分析流程參考Liuetal. (2021)和劉穎等(1996),分析精度優(yōu)于±5%。用于校準(zhǔn)測(cè)量樣品元素濃度的標(biāo)準(zhǔn)樣品為AGV-2,BHVO-2,BCR-2和RGM-2。

    3.4 全巖Sr-Nd同位素分析

    選擇具有代表性的6件樣品進(jìn)行了Sr和Nd同位素組成分析。Rb-Sr和Sm-Nd同位素的化學(xué)分離及同位素測(cè)試在武漢上譜分析科技有限責(zé)任公司完成,采用的儀器為美國(guó)Thermo Fisher Scientific 公司的MC-ICP-MS(Neptune Plus),所有分析樣品的86Sr/88Sr和146Nd/144Nd比值分別采用86Sr/88Sr=0.1194和146Nd/144Nd=0.7219進(jìn)行質(zhì)量分餾校正。分別使用USGS標(biāo)準(zhǔn)樣品NBS987和GSB校正所測(cè)樣品的Sr和Nd同位素比值。本次測(cè)得的NBS98787Sr/86Sr比值為0.710241±0.000012(2σ, N=4),GSB146Nd/144Nd比值為0.512441±0.000007(2σ, N=6),與先前發(fā)表的數(shù)據(jù)(87Sr/86Sr=0.710241±0.000012和143Nd/144Nd=0.512439±0.000007),在誤差范圍內(nèi)一致(Thirlwall, 1991;Lietal., 2017)。

    此外,本次分析的標(biāo)準(zhǔn)物質(zhì)BCR-2(玄武巖)和RGM-2(流紋巖)87Sr/86Sr比值分別為0.705002±0.000008(2σ, N=1)和0.704144±0.000008(2σ, N=1),143Nd/144Nd比值分別為0.512640±0.000007(2σ,N=1)和0.512806±0.000007(2σ,N=1),與之前發(fā)布的數(shù)值(BCR-2:87Sr/86Sr=0.705026±0.000020,143Nd/144Nd=0.512638±0.000015;RGM-2:87Sr/86Sr=0.704184±0.000020,143Nd/144Nd=0.512803±0.000015;Weisetal., 2006;Lietal., 2012)在誤差范圍內(nèi)一致。詳細(xì)的分析步驟參考Yangetal. (2019)。

    4 數(shù)據(jù)分析結(jié)果

    4.1 鋯石U-Pb-Hf同位素分析

    本研究選擇了3件寧陜花崗巖鋯石樣品進(jìn)行U-Pb-Hf同位素分析,樣品分別為NS18-8(胭脂壩巖體)、NS18-16(懶板凳巖體)和NS18-29(老城巖體)。鋯石U-Pb同位素?cái)?shù)據(jù)結(jié)果見(jiàn)表1。3件花崗巖樣品中的鋯石具有相似的結(jié)構(gòu)和形態(tài)特征,均表現(xiàn)為自形到半自形,粒度在100~300μm之間,長(zhǎng)/寬比在1:1.5~1:3之間。在CL圖像中,這些鋯石顆粒均具有明顯的震蕩環(huán)帶以及扇形環(huán)帶,且多數(shù)鋯石顆粒無(wú)繼承核(圖4a),屬于典型的巖漿成因鋯石。

    胭脂壩巖體中鋯石(NS18-8)的Th含量在182×10-6~620×10-6之間,U的含量變化范圍較大(351×10-6~2944×10-6),Th/U比值在0.10~0.58之間(表1)。共獲得30個(gè)測(cè)點(diǎn)數(shù)據(jù),除其中2個(gè)測(cè)點(diǎn)諧和度較低,8顆鋯石具有較老的206Pb/238U年齡(207~234Ma)外,其余20個(gè)測(cè)點(diǎn)的206Pb/238U年齡較為集中,獲得加權(quán)平均年齡為197.2±1.3Ma(N=20,MSWD=0.085)(圖5a)。懶板凳巖體中鋯石(NS18-16)具有較高的Th(1581×10-6~5101×10-6)和U(3211×10-6~7388×10-6)含量,Th/U比值在0.40~1.20之間(表1)。共獲得30個(gè)測(cè)點(diǎn)數(shù)據(jù),除去諧和度較低的6個(gè)測(cè)點(diǎn),及具有較老的206Pb/238U年齡(230~235Ma)的5個(gè)測(cè)點(diǎn)外,其余19個(gè)測(cè)點(diǎn)的206Pb/238U年齡較為集中,加權(quán)平均年齡為221.6±1.5Ma(N=19,MSWD=0.82)(圖5b)。老城巖體中鋯石(NS18-29)的Th(136×10-6~5097×10-6)和U(863×10-6~9360×10-6)含量較高且變化范圍較大,對(duì)應(yīng)的Th/U比值為0.09~0.65(表1)。共獲得30個(gè)測(cè)點(diǎn)數(shù)據(jù),除其中1個(gè)測(cè)點(diǎn)的諧和度較低外,其余29個(gè)測(cè)點(diǎn)的206Pb/238U年齡較為集中,加權(quán)平均年齡為209.3±2.3Ma(N=29,MSWD=1.9)(圖5c)。上述鋯石206Pb/238U加權(quán)平均年齡與對(duì)采自胭脂壩(NS18-8)、懶板凳(NS18-16)和老城巖體(NS18-29)的3件鋯石樣品進(jìn)行了Hf同位素分析,根據(jù)鋯石結(jié)晶年齡分別計(jì)算了初始176Hf/177Hf比值、εHf(t)值和Hf模式年齡,分析結(jié)果見(jiàn)表2。3件樣品中鋯石的176Yb/177Hf比值分別在0.020472~0.044924、0.038339~0.110436和0.029547~0.066995之間,176Lu/177Hf比值分別在0.000832~0.001901、0.001502~0.004371和0.001197~0.002653之間,εHf(t)值分別在-5.0~+0.8、-7.1~-3.5和-5.1~+1.2之間(圖6)。3件鋯石樣品均具有較老的Hf同位素二階段模式年齡(tDM2=1384~1061Ma,N=19;tDM2=1514~1316Ma,N=13;tDM2=1386~1055Ma,N=13)。

    表1 寧陜花崗巖鋯石U-Th-Pb定年數(shù)據(jù)

    續(xù)表1

    表2 寧陜花崗巖鋯石Hf同位素?cái)?shù)據(jù)

    表3 寧陜花崗巖獨(dú)居石U-Th-Pb定年數(shù)據(jù)

    續(xù)表3Continued Table 3測(cè)點(diǎn)號(hào)元素含量(×10-6)PbThUTh/U同位素比值年齡(Ma)207Pb206Pb1σ207Pb235U1σ206Pb238U1σ207Pb235U1σ206Pb238U1σ.172562 183856 14961 12.6 0.04862 0.00057 0.20337 0.00271 0.03034 0.00029 188.0 2.5 192.7 1.9 .181873 124038 15709 8.10 0.04919 0.00058 0.20244 0.00269 0.02985 0.00029 187.2 2.5 189.6 1.8 .191657 118255 11631 10.40.04978 0.00127 0.20591 0.00539 0.03000 0.00031 190.1 5.0 190.6 1.9 .201904 128755 14718 8.97 0.04914 0.00058 0.20318 0.00269 0.02998 0.00029 187.8 2.5 190.5 1.8 .212225 160077 14547 11.3 0.04911 0.00066 0.20047 0.00296 0.02960 0.00029 185.5 2.7 188.1 1.8 .222313 162422 15767 10.60.04911 0.00063 0.20417 0.00291 0.03015 0.00030 188.7 2.7 191.5 1.9 .232735 198402 16976 12.0 0.04941 0.00059 0.20627 0.00276 0.03028 0.00030 190.4 2.5 192.3 1.9 .242569 191600 12576 15.6 0.04916 0.00064 0.20562 0.00300 0.03034 0.00030 189.9 2.8 192.7 1.9 .252705 189636 16522 11.8 0.04947 0.00060 0.20603 0.00281 0.03020 0.00030 190.2 2.6 191.8 1.9 NS.18.8.011197 94877 3551 27.4 0.04925 0.00087 0.20007 0.00378 0.02946 0.00029 185.2 3.5 187.2 1.8 .021167 91502 3345 28.1 0.04730 0.00095 0.19248 0.00410 0.02952 0.00029 178.7 3.8 187.5 1.9 .031238 96941 4716 21.1 0.05028 0.00084 0.20314 0.00370 0.02930 0.00029 187.8 3.4 186.2 1.8 .041131 87087 3602 24.8 0.05012 0.00102 0.20305 0.00443 0.02938 0.00029 187.7 4.1 186.7 1.9 .051038 79830 3483 23.5 0.05012 0.00139 0.19968 0.00568 0.02890 0.00029 184.9 5.3 183.6 1.8 .06966.8 73431 3672 20.5 0.04964 0.00093 0.20332 0.00407 0.02971 0.00030 187.9 3.8 188.7 1.9 .071439 114268 4467 26.2 0.05013 0.00088 0.20435 0.00390 0.02956 0.00029 188.8 3.6 187.8 1.9 .081575 124206 4036 31.6 0.05053 0.00106 0.20597 0.00456 0.02956 0.00029 190.2 4.2 187.8 1.9 .091259 100406 3355 30.7 0.04884 0.00114 0.19996 0.00486 0.02969 0.00029 185.1 4.5 188.6 1.8 .101560 126684 3465 37.5 0.05042 0.00112 0.20126 0.00468 0.02895 0.00028 186.2 4.3 184.0 1.8 .111404 112168 3990 28.8 0.04893 0.00078 0.19651 0.00340 0.02913 0.00028 182.2 3.2 185.1 1.8 .121624 129735 4345 30.6 0.05063 0.00094 0.20305 0.00400 0.02909 0.00028 187.7 3.7 184.8 1.8 .131133 90948 3260 28.6 0.05017 0.00087 0.20254 0.00378 0.02928 0.00029 187.3 3.5 186.0 1.8 .141114 88398 2889 31.4 0.04994 0.00095 0.20701 0.00422 0.03006 0.00030 191.0 3.9 190.9 1.9 .151235 96741 3192 31.1 0.04802 0.00097 0.20157 0.00424 0.03045 0.00030 186.5 3.9 193.3 1.9 .161310 101786 4362 23.9 0.04947 0.00073 0.20332 0.00330 0.02981 0.00029 187.9 3.1 189.4 1.9 .171072 83777 3350 25.7 0.05035 0.00082 0.20559 0.00362 0.02961 0.00029 189.8 3.3 188.1 1.9 .181281 99632 4374 23.4 0.05011 0.00081 0.19998 0.00347 0.02894 0.00028 185.1 3.2 183.9 1.8 .191214 98541 3670 27.5 0.05003 0.00081 0.19840 0.00343 0.02876 0.00028 183.8 3.2 182.8 1.8 .201085 84271 3587 24.1 0.04770 0.00080 0.19820 0.00356 0.03014 0.00030 183.6 3.3 191.4 1.9 .211314 104395 3237 33.1 0.04978 0.00128 0.20190 0.00549 0.02942 0.00030 186.7 5.1 186.9 1.9 .221070 83473 3320 25.8 0.04900 0.00107 0.20153 0.00461 0.02983 0.00030 186.4 4.3 189.5 1.9 .231096 84447 4180 20.7 0.04912 0.00082 0.20299 0.00368 0.02997 0.00029 187.7 3.4 190.4 1.9 .241403 113410 3480 33.4 0.04970 0.00109 0.20492 0.00469 0.02990 0.00030 189.3 4.3 190.0 1.9 .251602 129740 3847 34.6 0.04847 0.00115 0.19915 0.00495 0.02980 0.00030 184.4 4.6 189.3 1.9 NS.18.12.012175 173641 3997 44.6 0.04951 0.00119 0.20754 0.00514 0.03040 0.00030 191.5 4.7 193.1 1.9 .021817 148471 4114 37.0 0.04887 0.00100 0.20429 0.00441 0.03032 0.00030 188.8 4.1 192.5 1.9 .032933 233059 7979 30.0 0.05081 0.00064 0.20867 0.00290 0.02978 0.00029 192.4 2.7 189.2 1.8 .042533 201477 5154 40.1 0.04999 0.00071 0.20487 0.00319 0.02972 0.00029 189.2 2.9 188.8 1.8 .054125 337276 8849 39.1 0.05255 0.00065 0.21710 0.00303 0.02996 0.00029 199.5 2.8 190.3 1.8 .062300 178949 6268 29.3 0.05037 0.00083 0.21296 0.00374 0.03066 0.00030 196.0 3.4 194.7 1.9 .071855 143486 5519 26.7 0.04953 0.00082 0.20927 0.00372 0.03064 0.00030 192.9 3.4 194.6 1.9

    圖4 寧陜巖體鋯石陰極發(fā)光(CL)圖像(a)及獨(dú)居石反射光圖像(b)Fig.4 Zircon cathodoluminescence (CL) images (a) and monazite reflected light images (b) for the Ningshan granitoids

    圖5 寧陜花崗巖鋯石(a-c)與獨(dú)居石(d-j)的U-Pb諧和圖Fig.5 Zircon (a-c) and monazite (d-j) U-Pb concordant diagrams of the Ningshan granitoids

    圖6 南秦嶺寧陜花崗巖中鋯石εHf(t)值與鋯石U-Pb年齡關(guān)系圖(據(jù)Liu et al., 2021)Fig.6 Zircon εHf (t) values versus U-Pb age diagram for the Ningshan granitoids in South Qinling (after Liu et al., 2021)

    4.2 獨(dú)居石U-Pb同位素分析

    本研究對(duì)胭脂壩(NS18-2、NS18-7、NS18-8、NS18-12和NS18-27)、懶板凳(NS18-16)和老城巖體(NS18-29)的7件獨(dú)居石樣品進(jìn)行了LA-MC-ICP-MS U-Pb同位素分析,分析結(jié)果見(jiàn)表3。測(cè)試的獨(dú)居石顆粒自形程度較好,粒徑在50~100μm之間(圖4b)。

    圖7 寧陜花崗巖的巖石分類(lèi)和系列圖解(a)A/NK-A/CNK圖解(據(jù)Maniar and Piccoli, 1989);(b)K2O-SiO2圖解(據(jù)Rickwood, 1989)Fig.7 Rock classification and series diagrams of the Ningshan granitoids(a) A/NK vs. A/CNK diagram (after Maniar and Piccoli,1989);(b) K2O vs. SiO2 diagram (after Peccerillo, 1989)

    采集自胭脂壩巖體的5件樣品NS18-2、NS18-7、NS18-8、NS18-12和NS18-27,加權(quán)平均年齡分別為191.8±1.4Ma(N=25,MSWD=3.0)、192.9±0.9Ma(N=25,MSWD=1.3)、187.0±1.0Ma(N=25,MSWD=1.7)、191.6±1.0Ma(N=24,MSWD=1.5)和182.4±1.4Ma(N=25,MSWD=2.5)(圖5d-g, i);分析的獨(dú)居石顆粒具有較高含量的Th和U,對(duì)應(yīng)的Th/U比值分別為0.29~41.54、6.45~15.24、20.00~36.57、21.69~47.40和14.32~36.42。懶板凳巖體獨(dú)居石樣品(NS18-16)的Th和U含量分別在85276×10-6~300341×10-6和3718×10-6~9382×10-6之間,Th/U比值為16.66~37.39,獲得的加權(quán)平均年齡為185.5±1.5Ma(N=25,MSWD=3.5)(圖5h)。老城巖體獨(dú)居石樣品(NS18-29)的Th和U含量分別為123230×10-6~206899×10-6和4715×10-6~9372×10-6,Th/U比值范圍在18.01~34.79之間,加權(quán)平均年齡為185.0±1.0Ma(N=25,MSWD=1.9)(圖5j)。

    4.3 全巖主量和微量元素地球化學(xué)特征

    本文測(cè)試了10件胭脂壩、9件懶板凳和4件老城巖體樣品的全巖元素地球化學(xué)組成,結(jié)果見(jiàn)表4。

    寧陜花崗巖樣品整體具有富硅(SiO2=68.46%~75.36%)、富鋁(Al2O3=13.96%~16.55%)、高鉀(K2O=3.00%~4.69%)和低Mg#(31~51)的特征,Na2O含量為2.56%~3.76%,CaO含量為1.13%~3.00%。相比之下,前人研究中老城巖體樣品的SiO2含量變化范圍較大,同時(shí)包含中性巖和酸性巖樣品;此外,前人的老城巖體樣品數(shù)據(jù)表現(xiàn)出更高的CaO和MgO含量(Yangetal., 2012;方博文等, 2017)。在A/CNK-A/NK圖解上,寧陜花崗巖樣品均落于過(guò)鋁質(zhì)系列范圍內(nèi)(A/CNK=1.10~1.35)(圖7a);從 SiO2-K2O巖石系列劃分圖(圖7b)上看,大部分樣品落入高鉀鈣堿性系列中。在哈克圖解中(圖8),隨著SiO2含量的增加,TiO2、Al2O3、P2O5、MgO和CaO含量降低,Nb/Ta、Dy/Yb比值降低,而Rb/Sr比值增加。

    在球粒隕石標(biāo)準(zhǔn)化稀土元素配分圖(圖9a,c,e)中,寧陜花崗巖樣品輕、重稀土分異明顯,具有富集輕稀土元素,虧損重稀土元素的特征,(La/Yb)N為2.48~31.63。胭脂壩與懶板凳巖體花崗巖的稀土元素總含量(∑REE)分別為138×10-6~228×10-6和102×10-6~195×10-6, 均具有明顯Eu負(fù)異常,Eu/Eu*分別為0.41~0.60和0.23~0.67;老城巖體花崗巖的稀土元素總含量(∑REE)較低(134×10-6~145×10-6),無(wú)明顯Eu異常。在原始地幔標(biāo)準(zhǔn)化微量元素蜘蛛圖(圖9b,d,f)中,所有寧陜花崗巖樣品均表現(xiàn)出富集Rb、Ba、Th、K、Pb和Sr等元素,而虧損Ta、Nb、P和Ti等元素的特征。

    表4 寧陜花崗巖元素地球化學(xué)組成(主量元素:wt%;微量元素:×10-6)

    續(xù)表4

    續(xù)表4

    續(xù)表4

    圖8 寧陜花崗巖SiO2與主量元素和微量元素比值圖解Fig.8 Diagrams of ratio of SiO2 against major and trace elements of the Ningshan granitoids

    圖9 寧陜花崗巖球粒隕石標(biāo)準(zhǔn)化稀土元素配分圖(a、c、e)和原始地幔標(biāo)準(zhǔn)化微量元素蜘蛛圖(b、d、f)(標(biāo)準(zhǔn)化值據(jù)Sun and McDonough, 1989)Fig.9 Chondrite-normalized REE patterns (a, c, e) and primitive mantle-normalized trace element spider diagrams (b, d, f) for the Ningshan granitoids (normalization values after Sun and McDonough, 1989)

    4.4 全巖Sr-Nd同位素地球化學(xué)特征

    寧陜花崗巖的Sr-Nd同位素分析結(jié)果詳見(jiàn)表5。其中,胭脂壩(~197Ma)、懶板凳(~221Ma)和老城(~209Ma)三個(gè)巖體的鋯石206Pb/238U加權(quán)平均年齡分別用于計(jì)算各自的初始Sr-Nd同位素比值。

    圖10 寧陜花崗巖的εNd(t)與年齡的關(guān)系圖解(據(jù)劉桂萍等,2021)Fig.10 Diagram of εNd(t) vs. age of the Ningshan granitoids (after Liu et al., 2021)

    胭脂壩巖體樣品Sr、Nd含量分別在294.10×10-6~351.85×10-6、28.33×10-6~32.94×10-6之間,具有較為均一的初始Sr同位素比值(87Sr/86Sr(t)=0.705348~0.705483),εNd(t) 值為-4.5~-4.7(圖10),二階段模式年齡范圍為1365~1343Ma。懶板凳巖體的樣品具有變化較大的Sr(111.81×10-6~408.16×10-6)和Nd(16.29×10-6~26.22×10-6)含量,初始87Sr/86Sr同位素比值在0.703337~0.705214之間,εNd(t)值(-1.8~-6.3)變化范圍較大(圖10),二階段模式年齡為1511~1147Ma。老城巖體樣品的Sr(815.15×10-6)含量較高,Nd含量為22.13×10-6,初始87Sr/86Sr同位素比值與εNd(t)值分別為0.704873和-3.9(圖10),二階段模式年齡為1309Ma。

    5 討論

    5.1 巖體形成時(shí)代及其地質(zhì)意義

    造山帶的形成記錄了板塊匯聚和造山過(guò)程的演化歷史(Zhangetal., 2015, 2017;王曉霞等,2015;Dengetal., 2017)。秦嶺造山帶作為我國(guó)中央造山帶的重要組成部分,自中生代以來(lái)經(jīng)歷了一系列重要的地質(zhì)事件,包括大陸裂解-洋盆產(chǎn)生、大洋俯沖-大陸增生、大陸碰撞、板內(nèi)構(gòu)造演化等(張國(guó)偉等,2001;張成立等,2008;Dongetal., 2011;Dengetal., 2016)。這些地質(zhì)事件的發(fā)生往往伴隨著廣泛而強(qiáng)烈的巖漿活動(dòng),使得秦嶺造山帶中發(fā)育大量的早中生代花崗巖(張成立等,2008;Dengetal., 2016)。

    近幾十年來(lái),前人對(duì)秦嶺造山帶的地球動(dòng)力學(xué)背景和演化過(guò)程進(jìn)行了大量研究,但對(duì)南秦嶺構(gòu)造帶三疊紀(jì)巖漿活動(dòng)的構(gòu)造環(huán)境的認(rèn)識(shí)仍然存在較大爭(zhēng)議,具體可以概括如下:(1)張成立等(2008)基于對(duì)早中生代花崗巖地球化學(xué)特征的研究,結(jié)合同時(shí)代煌斑巖和基性巖脈的雙峰式巖漿作用特點(diǎn),認(rèn)為秦嶺早中生代花崗巖(245~200Ma)屬于后碰撞階段的產(chǎn)物;(2)~211Ma古特斯洋盆閉合(Jiangetal., 2010),揚(yáng)子板塊北緣與南秦嶺構(gòu)造帶于晚三疊世發(fā)生碰撞(Liuetal., 2005;劉少峰和張國(guó)偉,2008);(3)252~234Ma期間南秦嶺構(gòu)造帶中的勉略洋仍處于洋內(nèi)俯沖階段(黃雄飛,2016;Dengetal., 2016;Huetal., 2019;Xingetal., 2020;Liuetal., 2021),~230Ma開(kāi)始碰撞造山(黃雄飛,2016),同碰撞階段為228~215Ma,后碰撞階段為215~200Ma(Dongetal., 2011)。然而,秦嶺造山帶的變形變質(zhì)事件表明勉略洋盆閉合時(shí)間應(yīng)早于215Ma(黃雄飛,2016)。王曉霞等(2015)對(duì)秦嶺早中生代花崗巖開(kāi)展了綜合研究,指出南秦嶺到北秦嶺不同地區(qū)花崗質(zhì)巖石的年齡、巖石類(lèi)型以及地球化學(xué)特征均沒(méi)有顯示出俯沖極性;楊朋濤等(2013)對(duì)南秦嶺何家莊花崗閃長(zhǎng)巖(~248Ma)開(kāi)展了詳細(xì)的研究,并提出勉略洋殼的俯沖在~248Ma之前就已啟動(dòng),并且南秦嶺地區(qū)在早三疊世仍處于洋殼俯沖動(dòng)力學(xué)背景。

    胭脂壩、懶板凳和老城巖體獨(dú)居石樣品的U-Pb年齡分別為194~181Ma、187~184Ma和~185Ma,均小于鋯石U-Pb年齡(圖11)。樣品NS18-16和NS18-29的鋯石與獨(dú)居石U-Pb年齡分別相差36Myr和24Myr。獨(dú)居石U-Pb體系的封閉溫度比鋯石低,因此獨(dú)居石可以記錄相對(duì)年輕的結(jié)晶年齡;此外,獨(dú)居石的U-Pb體系易受到后期熱事件的擾動(dòng),因此獨(dú)居石定年結(jié)果可能存在多解性(陳旭等,2009;吳黎光和李獻(xiàn)華,2020)。

    本研究的鋯石及獨(dú)居石的年齡結(jié)果指示寧陜花崗巖為晚三疊世到早侏羅世巖漿作用的產(chǎn)物,結(jié)合前人研究的定年結(jié)果(Yangetal., 2011;韋龍猛等, 2016;方博文等, 2017),本文將寧陜花崗質(zhì)巖漿活動(dòng)大體劃分為2個(gè)期次:早期為懶板凳巖體的二長(zhǎng)花崗巖,形成時(shí)代約為223~220Ma;晚期為胭脂壩和老城巖體的石英閃長(zhǎng)巖-花崗閃長(zhǎng)巖,形成時(shí)代約為211~196Ma。與Dongetal.(2011)、黃雄飛(2016)、Dengetal.(2016)、Huetal.(2019)、Xingetal.(2020)和Liuetal.(2021)等研究者對(duì)南秦嶺早中生代花崗巖的年齡劃分結(jié)果相近,故而本研究認(rèn)為寧陜花崗巖中較早期的花崗巖類(lèi)(懶板凳巖體)是南秦嶺構(gòu)造帶與揚(yáng)子板塊碰撞階段的產(chǎn)物,而寧陜巖體較晚期的花崗巖類(lèi)(胭脂壩巖體和老城巖體)形成于后碰撞階段。

    5.2 巖漿過(guò)程與巖漿源區(qū)

    花崗巖成分變化的控制因素主要包括源區(qū)組成、PAE和巖漿過(guò)程等(Zhengetal., 2008;Gaoetal., 2014)。具體而言,巖漿過(guò)程包括分離結(jié)晶作用、巖漿混合作用和地殼混染作用等。

    胭脂壩與懶板凳巖體的樣品具有明顯的Eu負(fù)異常,且隨著Sr含量的減少Rb/Sr比值增加(圖12a),這表明它們可能經(jīng)歷了斜長(zhǎng)石的分離結(jié)晶作用(Yangetal., 2016)。此外,Eu/Eu*與Sr之間的正相關(guān)性,以及Rb/Ba與Rb/Sr的關(guān)系圖解進(jìn)一步證明斜長(zhǎng)石分離結(jié)晶作用的影響(圖12b, f,Huetal., 2017;Bailieetal., 2020)。對(duì)于角閃石,KD(Nb)/ KD(Ta)和KD(Dy)/ KD(Yb)值均超過(guò)1.0,且胭脂壩與懶板凳巖體的樣品中Nb/Ta與Dy/Yb比值隨著SiO2含量的增加而減少,暗示了角閃石的分離結(jié)晶(圖8f、g,Nash and Crecraft, 1985;Stepanov and Hermann, 2013;Huetal., 2019)。此外,La/Yb與La/Sm的關(guān)系圖解(圖12e),進(jìn)一步證明胭脂壩與懶板凳巖體在巖漿演化過(guò)程中存在角閃石的分離結(jié)晶(Chenetal., 2018)。黑云母的分離結(jié)晶會(huì)導(dǎo)致Nb/Ta比值的減少,以及Ti的虧損(Stepanovetal., 2014),這與本文數(shù)據(jù)結(jié)果是一致的(圖8f、圖9b, d, f);V與Sc在黑云母中具有較高的分配系數(shù),Th的分配系數(shù)較低(Beaetal., 1994),胭脂壩與懶板凳巖體的樣品中Sc/Th和V/Th與SiO2/Al2O3之間呈負(fù)相關(guān)關(guān)系(圖12c, d),表明巖漿組分的變化與黑云母的分離結(jié)晶密切相關(guān)(Yangetal., 2016)。

    圖12 寧陜花崗巖熔體分離結(jié)晶與初始熔體的源區(qū)特征辨別圖(a)Rb/Sr-Sr(×10-6)圖解(據(jù)Yang et al., 2016);(b)Eu/Eu*-Sr(×10-6)圖解;(c)Sc/Th-SiO2/Al2O3圖解(據(jù)Yang et al., 2016);(d)V/Th-SiO2/Al2O3圖解(據(jù)Yang et al., 2016);(e)La/Yb-La/Sm圖解(據(jù)Chen et al., 2018);(f)Rb/Ba-Rb/Sr圖解(據(jù)Bailie et al., 2020). 玄武巖和泥質(zhì)巖衍生熔體之間的混合曲線據(jù)Patio Douce and Harris (1998)和Sylvester (1998)Fig.12 Discrimination diagrams of rock genesis and source characteristics of initial melt in Ningshan rock mass(a) Rb/Sr vs. Sr (×10-6) (after Yang et al., 2016); (b) Eu/Eu*vs. Sr (×10-6); (c) Sc/Th vs. SiO2/Al2O3(after Yang et al., 2016); (d) V/Th vs. SiO2/Al2O3(after Yang et al., 2016); (e) La/Yb vs. La/Sm (after Chen et al., 2018); (f) Rb/Ba vs. Rb/Sr (after Bailie et al., 2020). The mixing curve between basalt- and pelite-derived melts is from Patio Douce and Harris (1998) and Sylvester (1998)

    老城巖體的樣品中無(wú)明顯Eu異常,且Rb/Sr和Eu/Eu*與Sr之間無(wú)明顯相關(guān)性(圖12a, b),表明在巖漿演化過(guò)程中可能長(zhǎng)石未發(fā)生明顯分離結(jié)晶。Nb/Ta與Dy/Yb比值隨著SiO2含量的增加而減少,指示角閃石的分離結(jié)晶(圖8 f, g,Nash and Crecraft, 1985;Stepanov and Hermann, 2013;Huetal., 2019)。La/Yb與La/Sm的關(guān)系圖解(圖12e),也反映出存在角閃石的分離結(jié)晶(Chenetal., 2018)。圖12c和圖12d中Sc/Th和V/Th與SiO2/Al2O3之間呈負(fù)相關(guān)關(guān)系,這表明在老城巖體的巖漿演化過(guò)程中存在黑云母的分離結(jié)晶(Yangetal., 2016)。

    寧陜花崗巖SiO2含量較高且變化范圍較大(SiO2=68.46%~75.36%),同時(shí)它們還具有較低的Cr(0.94×10-6~36.76×10-6)和Ni(0.91×10-6~16.59×10-6)含量,這些特征表明這些高硅花崗巖可能并未受到明顯的幔源組分的影響。

    如F-An-Or圖所示(圖13),Castro (2013)整合了不同中酸性熔體的實(shí)驗(yàn)數(shù)據(jù),不同曲線代表了不同實(shí)驗(yàn)條件下熔體的演化路徑。此外,F(xiàn)-An-Or圖解還可用于指示地殼混染過(guò)程的影響;如果巖漿受到顯著的地殼混染,這將導(dǎo)致熔體成分顯著向遠(yuǎn)離An端員的方向演化(Castro,2013),例如雜砂巖與泥質(zhì)巖類(lèi)圍巖通常具有較低含量的鈣和較高含量的鋁。寧陜花崗巖樣品似乎主要沿著熔體分離結(jié)晶演化趨勢(shì)線分布,指示這些花崗巖可能受到同化混染的影響較小。此外,寧陜花崗巖具有比較均一的Sr-Nd同位素組成,也指示了同化混染作用不顯著。

    圖13 寧陜地區(qū)花崗巖類(lèi)F(FeO+MgO+MnO)-An-Or圖解a-混染;Ru-殘留-未混合(據(jù)Castro, 2013);參考數(shù)據(jù)來(lái)自Xing et al. (2020),Liu and Han (2018),張宏飛等(2005);Luo et al. (2018),Patio Douce and Beard(1995)Fig.13 F (FeO+MgO+MnO) -An-Or diagram of Ningshan granitoidsa-assimilation; Ru-restite unmixing (after Castro, 2013). Reference data from Xing et al. (2020), Liu and Han (2018), Zhang et al. (2005), Luo et al. (2018), and Patio Douce and Beard (1995)

    歸納而言,寧陜花崗巖類(lèi)可能經(jīng)歷了斜長(zhǎng)石、黑云母、角閃石等礦物相的分離結(jié)晶,并且未遭受顯著的地殼混染及幔源組分的影響。

    實(shí)驗(yàn)巖石學(xué)的結(jié)果表明:(1)地殼巖石(如泥質(zhì)變質(zhì)巖、變質(zhì)雜砂巖等)脫水熔融形成的深熔熔體具有高Si,低Fe、Mg、Ca、Ti的地球化學(xué)特征(Stevensetal., 2007);(2)變玄武巖和變安山巖的部分熔融可以產(chǎn)生過(guò)鋁質(zhì)熔體(Sissonetal., 2005);(3)變泥質(zhì)巖與變雜砂巖(硬砂巖)部分熔融亦可產(chǎn)生過(guò)鋁質(zhì)熔體(Vielzeuf and Holloway, 1988;Patio Douce and Johnston, 1991;Stevensetal., 1997;Pickering and Johnston, 1998)。

    本文所研究的高硅巖石樣品(SiO2>68%)具有較高的K2O(>3.00%)和LREE,較低的CaO(1.13%~3.00%)含量(表4),初始Sr同位素比值為0.703337~0.705483,εNd(t)值介于-6.3~-1.8之間。在陸內(nèi)造山帶中,弱過(guò)鋁質(zhì)(l-P區(qū)域)至中等過(guò)鋁質(zhì)花崗巖(m-P區(qū)域)的原巖可能是泥質(zhì)變質(zhì)巖、變質(zhì)火成巖和變質(zhì)雜砂巖(圖14h;Villasecaetal., 1998)。胭脂壩巖體樣品在高硅過(guò)鋁質(zhì)到強(qiáng)過(guò)鋁質(zhì)(f-P至h-P)區(qū)域內(nèi)均有分布,懶板凳巖體樣品主要分布于高硅過(guò)鋁質(zhì)(f-P)區(qū)域內(nèi)(圖14h);而老城巖體樣品在偏鋁質(zhì),中到強(qiáng)過(guò)鋁質(zhì)(m-P至h-P)的區(qū)域內(nèi)均有分布(圖14h)。如圖12f所示,寧陜花崗巖中的高硅樣品(SiO2>68%),源區(qū)組分主要為雜砂巖或砂屑巖,其中胭脂壩巖體的源區(qū)組分中主要含砂屑巖與雜砂巖,懶板巖體的源區(qū)組分中砂屑巖含量較高,老城巖體花崗巖(SiO2>68%,A/CNK>1.1)源區(qū)組分主要為雜砂巖。

    筆者將寧陜花崗巖類(lèi)與各種無(wú)水實(shí)驗(yàn)熔體組分進(jìn)行對(duì)比(圖14,T<950℃,P=0.5~1.5GPa)。在Yangetal.(2012)和方博文等(2017)的研究中,寧陜花崗巖樣品的源區(qū)可能是以變玄武巖或新元古代基性下地殼物質(zhì)為主;而本文研究中低鐵鎂質(zhì)的寧陜高硅花崗巖(SiO2>68%)在地球化學(xué)組分上表現(xiàn)出較高的A/CNK值和K含量,和較低的Ca和Al含量,與基性下地殼物質(zhì)組分差距較大,而與雜砂巖、火山碎屑雜砂巖和泥質(zhì)巖的成分更為匹配(圖14)。因此,本文推測(cè)寧陜高硅花崗巖樣品(SiO2>68%)的源區(qū)物質(zhì)可能主要以雜砂巖或砂屑巖為主(如Patio Douce and Beard (1996)和Montel and Vielzeuf (1997)的實(shí)驗(yàn)熔體;圖12、圖14)。

    此外,胭脂壩、懶板凳與老城巖體花崗巖的εHf(t)值分別在-7.1~-3.5、-5.0~+0.8和-5.1~+1.2之間(圖6),二兩階段的Hf模式年齡分別為tDM2=1514~1316Ma (N=13);tDM2=1384~1061Ma (N=19);和tDM2=1386~1055Ma (N=13),這可能指示了來(lái)自古老地殼源區(qū)的貢獻(xiàn)(Xingetal., 2020;Liuetal., 2021)。結(jié)合同位素分析結(jié)果及源區(qū)物質(zhì)特征,本研究認(rèn)為寧陜花崗巖主要為中元古代地殼來(lái)源。

    5.3 轉(zhuǎn)熔礦物組合的選擇性帶入(PAE)

    花崗巖成分變化的控制因素大體可以分為以下五種類(lèi)型:(1)源區(qū)的不均一性,如Clemens and Stevens(2012)的研究指出源區(qū)物質(zhì)是花崗質(zhì)巖漿化學(xué)組分的最主要控制因素,不同批次花崗巖中高度不均一的同位素組成也可以反映出源區(qū)物質(zhì)的高度不均一性(Wangetal., 2013b);(2)熔融條件 (P-T-H2O-fO2) 的變化, 多種實(shí)驗(yàn)巖石學(xué)結(jié)果表明熔融溫度和H2O含量可以控制花崗巖體成分變化(Weinberg and Hasalová, 2015;Zhaoetal., 2015;Gaoetal., 2016);(3)巖漿混合,多數(shù)學(xué)者認(rèn)為花崗巖是由幔源基性巖漿與殼源酸性巖漿混合形成,其成分受控于二者的混合比例(Collins, 1996;Yangetal., 2007;Shaw and Flood, 2009);(4)分離結(jié)晶和同化混染,花崗巖成分變化受到圍巖的混染比例以及分離結(jié)晶程度的影響(Kempetal., 2007;Castroetal., 2021);(5)轉(zhuǎn)熔礦物組合的選擇性帶入(PAE,Peritectic Assemblage Entrainment)亦可對(duì)花崗質(zhì)巖漿成分造成影響,該過(guò)程主要受控于轉(zhuǎn)熔礦物類(lèi)型和被花崗質(zhì)熔體選擇性帶入的比例(Clemens and Stevens, 2012;Bailieetal., 2020;Zhuetal., 2020,2021)等。

    圖14 南秦嶺寧陜花崗巖源區(qū)組成判別圖解

    PAE模型是指源區(qū)礦物的不一致熔融后形成一些轉(zhuǎn)熔礦物和熔體,二者同時(shí)存在,其中轉(zhuǎn)熔礦物組合是從巖漿源區(qū)被攜帶進(jìn)入花崗巖熔體中的,并且是花崗質(zhì)熔體鎂鐵質(zhì)組分變化的主要控制因素(Clemens and Stevens, 2012;Bailieetal., 2020;Zhuetal., 2020,2021)。PAE模型的地球化學(xué)依據(jù)是主、微量元素(如Ti、Ca)與鎂鐵指數(shù)(摩爾Fe+Mg)的相關(guān)性,轉(zhuǎn)熔礦物相是指由含水礦物相(如地殼源區(qū)中黑云母和角閃石礦物相)不一致熔融產(chǎn)生的鎂鐵礦物和斜長(zhǎng)石(如黑云母+角閃石+石英+斜長(zhǎng)石1=熔體+斜長(zhǎng)石2+單斜輝石+斜方輝石+鈦鐵礦±石榴子石(斜長(zhǎng)石1:反應(yīng)礦物相中的斜長(zhǎng)石;斜長(zhǎng)石2:轉(zhuǎn)熔礦物相的斜長(zhǎng)石),Skjerlie and Johnston, 1992;Stevensetal., 2007;Clemensetal., 2011;Clemens and Stevens, 2012)。

    胭脂壩、懶板凳與老城巖體中的花崗巖具有高的SiO2含量(SiO2>68%),這些花崗巖中Ti和Ca均顯示出與Fe+Mg值的高度正相關(guān)性(圖15a, b),以及A/CNK顯示出與Fe+Mg值的負(fù)相關(guān)性(圖15c),反映出轉(zhuǎn)熔鎂鐵質(zhì)組分增加對(duì)花崗巖成分的顯著影響,故而PAE很可能是寧陜花崗巖全巖地球化學(xué)成分變化的一個(gè)重要控制因素(Clemensetal., 2011;Clemens and Stevens, 2012;Farinaetal., 2012)。前文指出,寧陜花崗巖主要為地殼來(lái)源,而黑云母或角閃石是地殼熔融過(guò)程中的主要含K礦物相,胭脂壩、懶板凳與老城高硅花崗巖(SiO2>68%)中K2O/Na2O比值分別為0.92~1.78、1.00~1.66和0.81~1.31(主要集中在1.00~1.78之間),表明在3個(gè)單元花崗巖的源區(qū)熔融反應(yīng)過(guò)程中黑云母或角閃石是主要的反應(yīng)礦物(Clemensetal., 2011;Farinaetal., 2012)。此外,V與Fe+Mg值之間明顯的正相關(guān)性也表明寧陜花崗巖的地殼源區(qū)中黑云母或/和角閃石的含量相對(duì)較高(圖15d,Bailieetal., 2020)。圖14h表明,寧陜花崗巖熔融反應(yīng)中可能沒(méi)有角閃石的參與(Villasecaetal., 1998;Bailieetal., 2020)。此外,Clemensetal. (2011)的模擬結(jié)果顯示(圖15),地殼源區(qū)中不同黑云母/角閃石的比例可能導(dǎo)致花崗質(zhì)熔體中Ti與Fe+Mg值的斜率范圍存在差異,胭脂壩、懶板凳與老城巖體的高硅花崗巖樣品(SiO2>68%)明顯更偏向于含黑云母的地殼源區(qū)(Clemensetal., 2011;Clemens and Stevens, 2012)。

    圖15 南秦嶺寧陜巖體中生代花崗巖的Ti-(Fe+Mg)(a)、Ca-(Fe+Mg)(b)、A/CNK-(Fe+Mg)(c)、V-Ti(d)、Zr-(Fe+Mg)(e)、Hf-(Fe+Mg)(f)、P-(Fe+Mg)(g)、P-Ca(h)和Ti-Ca(i)的關(guān)系圖(據(jù)Clemens et al.,2011)圖15a(據(jù)Clemens et al., 2011)關(guān)系圖顯示轉(zhuǎn)熔礦物裹帶單斜輝石、斜長(zhǎng)石、鈦鐵礦、鋯石、磷灰石和角閃石,陰影區(qū)域顯示了含有黑云母和角閃石的巖石部分熔融產(chǎn)生的地球化學(xué)變化; 實(shí)線表示變質(zhì)泥質(zhì)巖石中黑云母分解產(chǎn)生轉(zhuǎn)熔石榴石和鈦鐵礦的趨勢(shì); 以長(zhǎng)虛線標(biāo)記的趨勢(shì)表示初始熔體與轉(zhuǎn)熔石榴石、單斜輝石和鈦鐵礦之間的混合,這些礦物的比例由巖石與黑云母和角閃石的比例為3:2的部分熔融產(chǎn)生; 以短虛線標(biāo)記的最終趨勢(shì),表示初始熔體與轉(zhuǎn)熔斜方輝石、單斜輝石、石榴石和鈦鐵礦之間的混合,在含有黑云母和角閃石的巖石部分熔融期間形成,比例為2:3Fig.15 The diagrams of Ti vs. (Fe+Mg) (a), Ca vs. (Fe+Mg) (b), A/CNK vs. (Fe+Mg) (c), V vs. Ti (d), Zr vs. (Fe+Mg) (e), Hf vs. (Fe+Mg) (f), P vs. (Fe+Mg) (g), P vs. Ca (h) and Ti vs. Ca (i) for Mesozoic granites of Ningshan pluton in South Qinling (after Clemens et al., 2011) to display peritectic assemblage entrainment of clinopyroxene, plagioclase, ilmenite, zircon, apatite and spheneIn Fig.5a (modified after Clemens et al., 2011), the shadow area displays the geochemical variation produced by partial melting of a rock with biotite and hornblende; the solid line represents the trend that would result from the breakdown of biotite in a metapelitic rock, producing peritectic garnet and ilmenite; the trend marked with long dashes represents mixing between the initial melt and peritectic garnet, clinopyroxene and ilmenite in the ratio of these minerals produced by partial melting of a rock with biotite and hornblende in the ratio 3:2; the final trend, marked with short dashes, represents mixing between the initial melt and peritectic orthopyroxene, clinopyroxene, garnet and ilmenite formed during partial melting of a rock containing biotite and hornblende in the ratio 2:3

    如圖15b, c所示,單斜輝石作為轉(zhuǎn)熔礦物夾帶將會(huì)導(dǎo)致Ca與Fe+Mg值的呈正相關(guān)關(guān)系,以及向偏鋁質(zhì)和高Fe+Mg值方向演化的趨勢(shì)(Clemensetal., 2011;Farinaetal., 2012)。本文樣品中Ca與Fe+Mg值的高度正相關(guān)性,以及A/CNK與Fe+Mg值之間的明顯負(fù)相關(guān)性,表明寧陜花崗巖(SiO2>68%)可能經(jīng)歷了單斜輝石轉(zhuǎn)熔礦物的裹帶。

    除了主要的轉(zhuǎn)熔礦物相之外,含黑云母源區(qū)的不一致熔融也會(huì)產(chǎn)生轉(zhuǎn)熔副礦物相的裹帶(Clemensetal., 2011;Clemens and Stevens, 2012;Farinaetal., 2012;Bailieetal., 2020)。本文中,寧陜花崗巖(SiO2>68%)的Ti與Fe+Mg值的高度正相關(guān)性可能反映了轉(zhuǎn)熔礦物鈦鐵礦的裹帶(圖15a;Bailieetal., 2020;Zhuetal., 2021)。La Tourretteetal.(1991)和Farinaetal.(2012)的研究表明,V在鈦鐵礦中具有很高的相容性(KD≈80),而寧陜花崗巖中Ti和V的正相關(guān)關(guān)系也反應(yīng)出轉(zhuǎn)熔鈦鐵礦的裹帶(圖15d;Clemens and Stevens, 2012;Farinaetal., 2012;Bailieetal., 2020;Zhuetal., 2021)。胭脂壩與懶板凳巖體中Zr和Hf與Fe+Mg值之間存在正相關(guān)性(圖15e, f),暗示源區(qū)中有轉(zhuǎn)熔鋯石的存在(Stevensetal., 2007;Bailieetal., 2020;Zhuetal., 2020, 2021),而老城巖體(SiO2>68%)中Zr和Hf與Fe+Mg值之間無(wú)明顯正相關(guān)關(guān)系,表明其源區(qū)可能不存在轉(zhuǎn)熔鋯石。此外,圖15g-i中寧陜花崗巖的P與Fe+Mg值、P與Ca以及Ti與Ca的正相關(guān)關(guān)系,可能反映了源區(qū)存在含Ca副礦物相(如磷灰石和榍石)的裹帶(Farinaetal., 2012;Bailieetal., 2020;Zhuetal., 2020, 2021)。

    綜上所述,本研究認(rèn)為,寧陜花崗巖是由含黑云母的殼源巖石部分熔融形成的,母巖漿內(nèi)可能存在單斜輝石和副礦物(如鈦鐵礦、鋯石、磷灰石和榍石)等轉(zhuǎn)熔礦物組合的選擇性帶入,老城巖體的母巖漿中可能沒(méi)有轉(zhuǎn)熔鋯石。

    6 結(jié)論

    根據(jù)鋯石U-Pb年齡結(jié)果,可將寧陜花崗巖的巖漿作用劃分為2個(gè)期次:早期為懶板凳巖體的二長(zhǎng)花崗巖,形成時(shí)代約為223~220Ma;晚期為胭脂壩和老城巖體的石英閃長(zhǎng)巖-花崗閃長(zhǎng)巖,形成時(shí)代約為211~196Ma,為晚三疊世到早侏羅世巖漿作用的產(chǎn)物。推斷較早期花崗巖(懶板凳巖體)形成于南秦嶺構(gòu)造帶與揚(yáng)子板塊同碰撞階段,而較晚期的花崗巖(胭脂壩巖體和老城巖體)形成于后碰撞階段。

    將寧陜花崗巖與多種無(wú)水實(shí)驗(yàn)熔體組分進(jìn)行對(duì)比,并結(jié)合主、微量元素與鎂鐵質(zhì)(摩爾Fe+Mg)的相關(guān)性,本研究認(rèn)為胭脂壩巖體源區(qū)組成以含黑云母砂屑巖和雜砂巖為主,懶板凳巖體的源區(qū)組成以含黑云母砂屑巖為主,二者演化過(guò)程中都伴隨有單斜輝石以及副礦物(如鈦鐵礦+鋯石+磷灰石+榍石)等轉(zhuǎn)熔礦物組合的選擇性帶入,且?guī)r漿演化過(guò)程中都存在斜長(zhǎng)石、角閃石與黑云母的分離結(jié)晶。老城巖體的源巖類(lèi)型主要為含黑云母雜砂巖,演化過(guò)程中伴有單斜輝石以及副礦物(鈦鐵礦+磷灰石+榍石)等轉(zhuǎn)熔礦物組合的選擇性帶入,巖漿演化過(guò)程中存在角閃石和黑云母的分離結(jié)晶。

    致謝感謝周紅英教授級(jí)高級(jí)工程師、陳育曉博士、朱毓博士、謝元惠博士和豆敬兆博士在研究過(guò)程中的幫助以及指導(dǎo);感謝審稿人提出的建議。

    猜你喜歡
    獨(dú)居石源區(qū)老城
    獨(dú)居石資源開(kāi)發(fā)利用現(xiàn)狀及政策建議
    獨(dú)居石堿分解熱力學(xué)研究
    電子探針技術(shù)研究大別造山帶富硫獨(dú)居石地球化學(xué)特征及稀土礦化成因
    復(fù)合捕收劑與組合抑制劑對(duì)微細(xì)粒獨(dú)居石與螢石浮選分離的作用機(jī)理研究
    老城舊影 永平門(mén)
    紅巖春秋(2022年11期)2022-05-30 16:00:02
    冬小麥蒸散源區(qū)代表性分析
    老城舊影·和平路
    紅巖春秋(2019年12期)2020-01-03 06:20:56
    老城活力
    人民交通(2019年16期)2019-12-20 07:03:42
    渭河源區(qū)徑流量變化特征及趨勢(shì)分析
    開(kāi)羅伊斯蘭教老城
    久久久久精品久久久久真实原创| 久久久久国产精品人妻一区二区| av在线蜜桃| 欧美bdsm另类| 国产精品一二三区在线看| 国产av不卡久久| 亚洲av欧美aⅴ国产| 真实男女啪啪啪动态图| 王馨瑶露胸无遮挡在线观看| 一区二区av电影网| 欧美xxⅹ黑人| 亚洲欧美日韩另类电影网站 | 久久久久久久久久久免费av| 美女高潮的动态| 一本色道久久久久久精品综合| 特级一级黄色大片| 久久精品久久久久久久性| 亚洲精品久久午夜乱码| 网址你懂的国产日韩在线| 国产精品久久久久久精品电影小说 | 日韩一本色道免费dvd| 国产欧美日韩一区二区三区在线 | 91久久精品国产一区二区三区| 国产精品久久久久久精品电影小说 | 成年女人在线观看亚洲视频 | 在线a可以看的网站| av在线app专区| 国内精品宾馆在线| 综合色av麻豆| 日韩不卡一区二区三区视频在线| 三级经典国产精品| 国产日韩欧美在线精品| 在现免费观看毛片| 偷拍熟女少妇极品色| 久久久久精品久久久久真实原创| 高清午夜精品一区二区三区| 80岁老熟妇乱子伦牲交| 色视频www国产| 久久6这里有精品| 日韩成人av中文字幕在线观看| 国产一区亚洲一区在线观看| 国产高潮美女av| 99久久人妻综合| 男女啪啪激烈高潮av片| 久久久欧美国产精品| 亚洲精品国产av蜜桃| 春色校园在线视频观看| 青青草视频在线视频观看| 亚洲伊人久久精品综合| 亚洲一区二区三区欧美精品 | 国产老妇伦熟女老妇高清| 高清av免费在线| 91aial.com中文字幕在线观看| 国产精品一区二区在线观看99| 男女无遮挡免费网站观看| 极品教师在线视频| 亚洲精品国产av成人精品| 亚洲欧美日韩卡通动漫| 高清日韩中文字幕在线| 国产精品99久久久久久久久| 熟女人妻精品中文字幕| 亚洲精品影视一区二区三区av| 久久久久久伊人网av| 国产精品无大码| 国产免费一级a男人的天堂| 成人特级av手机在线观看| 男人舔奶头视频| 亚洲经典国产精华液单| 18+在线观看网站| 丝袜喷水一区| 亚洲成人久久爱视频| 欧美一区二区亚洲| 国产精品人妻久久久影院| 91aial.com中文字幕在线观看| 寂寞人妻少妇视频99o| 看十八女毛片水多多多| 国产又色又爽无遮挡免| 亚洲精品一区蜜桃| 久久久久国产精品人妻一区二区| 亚洲最大成人中文| 又爽又黄a免费视频| 51国产日韩欧美| 欧美日韩国产mv在线观看视频 | 一级爰片在线观看| 干丝袜人妻中文字幕| 春色校园在线视频观看| 国产高清国产精品国产三级 | 边亲边吃奶的免费视频| 在线播放无遮挡| 直男gayav资源| 色视频在线一区二区三区| av在线亚洲专区| 亚洲图色成人| 久久精品夜色国产| 日韩免费高清中文字幕av| 国产成人精品婷婷| 亚洲美女视频黄频| 国产精品一区二区性色av| 夫妻性生交免费视频一级片| 亚洲国产精品999| 久久这里有精品视频免费| 成年av动漫网址| 麻豆乱淫一区二区| a级毛色黄片| 成人特级av手机在线观看| 国产精品秋霞免费鲁丝片| 日本黄大片高清| 国产黄片美女视频| 伦精品一区二区三区| 欧美+日韩+精品| 国产在线男女| 别揉我奶头 嗯啊视频| 国产欧美日韩一区二区三区在线 | 日本黄大片高清| 搞女人的毛片| 熟女av电影| 成人二区视频| 国产精品久久久久久av不卡| 国产av不卡久久| 国产成年人精品一区二区| 丰满人妻一区二区三区视频av| 国产 一区精品| 国产亚洲5aaaaa淫片| 22中文网久久字幕| 日本免费在线观看一区| 亚洲成色77777| 免费大片18禁| 亚洲av成人精品一二三区| 亚洲精品aⅴ在线观看| 波野结衣二区三区在线| 欧美一区二区亚洲| 成人毛片60女人毛片免费| 精品人妻视频免费看| 一级毛片电影观看| 中国三级夫妇交换| 亚洲欧洲国产日韩| 亚洲熟女精品中文字幕| 成人毛片a级毛片在线播放| 永久免费av网站大全| 欧美激情久久久久久爽电影| 午夜免费男女啪啪视频观看| 久久久久精品性色| 久久久久久国产a免费观看| 在线观看三级黄色| 国产精品国产三级国产av玫瑰| 久久久久精品久久久久真实原创| 国产一区二区三区av在线| 国产成人freesex在线| 18+在线观看网站| 小蜜桃在线观看免费完整版高清| 老师上课跳d突然被开到最大视频| 亚洲成人久久爱视频| 简卡轻食公司| 好男人在线观看高清免费视频| 欧美xxxx黑人xx丫x性爽| 啦啦啦中文免费视频观看日本| 大片电影免费在线观看免费| 在线播放无遮挡| tube8黄色片| 国产大屁股一区二区在线视频| 久久精品国产亚洲av天美| 高清在线视频一区二区三区| 亚洲三级黄色毛片| 夫妻午夜视频| 爱豆传媒免费全集在线观看| 18禁在线播放成人免费| 亚洲精品国产色婷婷电影| 男女边吃奶边做爰视频| 国产成年人精品一区二区| 国产精品久久久久久精品电影| 日本色播在线视频| 国产 一区精品| 哪个播放器可以免费观看大片| 欧美日韩视频精品一区| 另类亚洲欧美激情| 日韩欧美精品v在线| 22中文网久久字幕| 亚洲av日韩在线播放| 中文天堂在线官网| 我的老师免费观看完整版| 日韩欧美 国产精品| 亚洲天堂国产精品一区在线| 97精品久久久久久久久久精品| 国产精品嫩草影院av在线观看| 亚洲国产精品成人综合色| 国产亚洲最大av| 激情五月婷婷亚洲| 丰满乱子伦码专区| 日本爱情动作片www.在线观看| 日韩精品有码人妻一区| 久久精品久久久久久久性| 亚洲av电影在线观看一区二区三区 | 熟女电影av网| av福利片在线观看| 欧美三级亚洲精品| 成年女人看的毛片在线观看| 高清午夜精品一区二区三区| 亚洲国产高清在线一区二区三| 肉色欧美久久久久久久蜜桃 | 欧美性感艳星| 午夜免费鲁丝| 超碰av人人做人人爽久久| 最近中文字幕高清免费大全6| 一区二区三区四区激情视频| 久久97久久精品| 最近中文字幕高清免费大全6| 久久99热这里只频精品6学生| 91久久精品国产一区二区三区| 噜噜噜噜噜久久久久久91| 偷拍熟女少妇极品色| 男女国产视频网站| 欧美性猛交╳xxx乱大交人| 成人漫画全彩无遮挡| 在线观看美女被高潮喷水网站| 偷拍熟女少妇极品色| 我的老师免费观看完整版| 国产免费一级a男人的天堂| 一级毛片黄色毛片免费观看视频| 欧美日韩亚洲高清精品| 欧美3d第一页| 蜜桃久久精品国产亚洲av| 乱码一卡2卡4卡精品| 国产亚洲av片在线观看秒播厂| 一级爰片在线观看| 看黄色毛片网站| 国产乱来视频区| 建设人人有责人人尽责人人享有的 | 欧美3d第一页| 特大巨黑吊av在线直播| 亚洲精品一区蜜桃| 国产成人一区二区在线| 成人免费观看视频高清| 午夜福利网站1000一区二区三区| 亚洲精品日韩在线中文字幕| 久久这里有精品视频免费| 国产成人a区在线观看| 亚洲精品一二三| 色综合色国产| 交换朋友夫妻互换小说| 日本色播在线视频| 五月玫瑰六月丁香| 黄色怎么调成土黄色| 插阴视频在线观看视频| 亚洲成人一二三区av| 秋霞在线观看毛片| 午夜精品国产一区二区电影 | 国产精品国产av在线观看| 1000部很黄的大片| 国产一区二区在线观看日韩| 国产精品人妻久久久影院| 在线观看国产h片| 六月丁香七月| 国产v大片淫在线免费观看| 一边亲一边摸免费视频| 校园人妻丝袜中文字幕| 一级毛片我不卡| 日韩 亚洲 欧美在线| a级毛色黄片| 人妻 亚洲 视频| 国产 一区 欧美 日韩| 成人综合一区亚洲| 精品久久久久久久人妻蜜臀av| 亚洲最大成人中文| 身体一侧抽搐| 日韩人妻高清精品专区| 深夜a级毛片| 国产精品一及| 一级片'在线观看视频| 国产黄a三级三级三级人| 日本午夜av视频| 亚洲第一区二区三区不卡| 亚洲一级一片aⅴ在线观看| 日本午夜av视频| 免费观看在线日韩| 欧美+日韩+精品| 一级爰片在线观看| 最近最新中文字幕免费大全7| 国产成人freesex在线| 高清在线视频一区二区三区| 美女高潮的动态| 丝袜脚勾引网站| 天天躁日日操中文字幕| 搡老乐熟女国产| 日韩精品有码人妻一区| av在线播放精品| 下体分泌物呈黄色| 免费看av在线观看网站| 亚洲人成网站高清观看| 97超碰精品成人国产| 国产精品一区二区在线观看99| 国产伦精品一区二区三区视频9| 久久精品国产自在天天线| 午夜激情久久久久久久| 一区二区三区四区激情视频| 色婷婷久久久亚洲欧美| 亚洲久久久久久中文字幕| 欧美老熟妇乱子伦牲交| 欧美日韩一区二区视频在线观看视频在线 | 国产欧美日韩一区二区三区在线 | 最近的中文字幕免费完整| 日本爱情动作片www.在线观看| 久久久久九九精品影院| 色视频www国产| 免费人成在线观看视频色| 欧美区成人在线视频| 中文欧美无线码| 国产黄a三级三级三级人| 大香蕉97超碰在线| 狂野欧美激情性bbbbbb| tube8黄色片| 亚洲国产日韩一区二区| 欧美精品一区二区大全| 亚洲精品亚洲一区二区| 久久精品国产亚洲av涩爱| 18禁在线无遮挡免费观看视频| 99久久九九国产精品国产免费| 超碰97精品在线观看| 嫩草影院入口| 国产片特级美女逼逼视频| 一级毛片久久久久久久久女| 有码 亚洲区| 亚洲欧美日韩另类电影网站 | 在线免费观看不下载黄p国产| av在线亚洲专区| 制服丝袜香蕉在线| 精品国产乱码久久久久久小说| 18禁在线播放成人免费| 国产淫片久久久久久久久| 亚洲国产精品国产精品| 国产成人freesex在线| 久久久午夜欧美精品| 99久久中文字幕三级久久日本| 免费大片18禁| 18禁动态无遮挡网站| 色5月婷婷丁香| 身体一侧抽搐| 日日撸夜夜添| 中国国产av一级| 亚洲天堂av无毛| 亚州av有码| 一区二区三区四区激情视频| 狂野欧美激情性bbbbbb| tube8黄色片| 亚洲成色77777| 哪个播放器可以免费观看大片| 91久久精品国产一区二区三区| 美女cb高潮喷水在线观看| 天美传媒精品一区二区| 国产精品女同一区二区软件| 欧美性感艳星| 欧美高清成人免费视频www| 成人鲁丝片一二三区免费| 又爽又黄无遮挡网站| 日本黄色片子视频| 国产精品久久久久久精品古装| 成年人午夜在线观看视频| 精品人妻熟女av久视频| 国产视频内射| 日本一二三区视频观看| 国产色婷婷99| 伦精品一区二区三区| 午夜福利高清视频| 大片电影免费在线观看免费| 亚洲精品日本国产第一区| 天堂网av新在线| 欧美xxⅹ黑人| 少妇的逼好多水| 国产精品久久久久久精品电影小说 | 欧美丝袜亚洲另类| 18禁裸乳无遮挡免费网站照片| 韩国高清视频一区二区三区| 国产在线男女| 亚洲精品aⅴ在线观看| 亚洲自偷自拍三级| 日韩免费高清中文字幕av| 高清毛片免费看| 国产精品.久久久| 少妇人妻久久综合中文| 2021少妇久久久久久久久久久| 中文字幕人妻熟人妻熟丝袜美| 网址你懂的国产日韩在线| 免费播放大片免费观看视频在线观看| 免费观看无遮挡的男女| 国产久久久一区二区三区| 黄色日韩在线| 国产精品一及| 丰满少妇做爰视频| 亚洲国产精品国产精品| 成人亚洲精品一区在线观看 | 国产精品麻豆人妻色哟哟久久| 国产精品一区二区在线观看99| 美女国产视频在线观看| 亚洲四区av| 成人综合一区亚洲| 美女脱内裤让男人舔精品视频| 亚洲自偷自拍三级| 高清欧美精品videossex| 一级av片app| 成年版毛片免费区| 精品久久久久久电影网| 一本久久精品| 国产淫语在线视频| 亚洲精品,欧美精品| 嫩草影院新地址| 国产在线一区二区三区精| 三级男女做爰猛烈吃奶摸视频| 国产一区二区三区av在线| 午夜亚洲福利在线播放| 丝袜美腿在线中文| 特大巨黑吊av在线直播| 色综合色国产| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久久精品性色| 97超碰精品成人国产| 日本色播在线视频| 一本一本综合久久| 久久精品国产鲁丝片午夜精品| 国产免费福利视频在线观看| 亚洲精品国产色婷婷电影| 又爽又黄无遮挡网站| 少妇高潮的动态图| 亚洲美女视频黄频| 日韩亚洲欧美综合| 熟女av电影| 99热国产这里只有精品6| 亚洲综合精品二区| 国产午夜精品一二区理论片| 自拍偷自拍亚洲精品老妇| 毛片一级片免费看久久久久| 国产在线一区二区三区精| 国产女主播在线喷水免费视频网站| 一级a做视频免费观看| 五月天丁香电影| 精品人妻一区二区三区麻豆| 亚洲av日韩在线播放| 一级二级三级毛片免费看| 成年av动漫网址| 我的老师免费观看完整版| 天天躁日日操中文字幕| 成人亚洲精品av一区二区| 三级男女做爰猛烈吃奶摸视频| 在线观看免费高清a一片| 在线亚洲精品国产二区图片欧美 | 熟女av电影| 成人国产麻豆网| av女优亚洲男人天堂| 99re6热这里在线精品视频| 七月丁香在线播放| 少妇人妻精品综合一区二区| 亚洲国产日韩一区二区| 国产精品国产三级国产av玫瑰| 亚洲精品第二区| 国产精品一区二区三区四区免费观看| 国内揄拍国产精品人妻在线| 国产精品三级大全| 干丝袜人妻中文字幕| 欧美区成人在线视频| 成人特级av手机在线观看| 免费黄网站久久成人精品| 能在线免费看毛片的网站| 啦啦啦在线观看免费高清www| 午夜精品国产一区二区电影 | 在线a可以看的网站| 欧美另类一区| 亚洲精品国产av蜜桃| 亚洲精品日韩在线中文字幕| 国产精品久久久久久精品电影小说 | 禁无遮挡网站| 亚洲va在线va天堂va国产| 三级国产精品片| 成人亚洲欧美一区二区av| 国国产精品蜜臀av免费| 欧美日韩亚洲高清精品| 国产有黄有色有爽视频| 舔av片在线| 亚洲美女搞黄在线观看| 男的添女的下面高潮视频| 综合色丁香网| 在线观看人妻少妇| 亚洲精品,欧美精品| 亚洲精品影视一区二区三区av| 国产日韩欧美亚洲二区| 2022亚洲国产成人精品| 黄片无遮挡物在线观看| 国产亚洲最大av| 内地一区二区视频在线| 国内少妇人妻偷人精品xxx网站| 深爱激情五月婷婷| 久久久久久久久久久丰满| 亚洲精品视频女| 成人无遮挡网站| 亚洲无线观看免费| 日韩大片免费观看网站| 国产视频内射| 欧美激情国产日韩精品一区| 国产一区二区在线观看日韩| 晚上一个人看的免费电影| 国产精品不卡视频一区二区| 99久久中文字幕三级久久日本| 韩国高清视频一区二区三区| 大陆偷拍与自拍| 超碰av人人做人人爽久久| 国产色爽女视频免费观看| 深夜a级毛片| 亚洲欧美中文字幕日韩二区| 一个人看的www免费观看视频| 久久久久久国产a免费观看| 亚洲精品乱码久久久v下载方式| 亚洲精品自拍成人| 成人无遮挡网站| 可以在线观看毛片的网站| 国产精品av视频在线免费观看| 少妇熟女欧美另类| 亚洲av二区三区四区| 在线观看国产h片| 乱系列少妇在线播放| 国产欧美日韩精品一区二区| 亚洲美女视频黄频| 韩国av在线不卡| 亚洲欧美一区二区三区黑人 | 久久97久久精品| 高清在线视频一区二区三区| 2021天堂中文幕一二区在线观| 亚洲av免费高清在线观看| 黄色日韩在线| 最近中文字幕高清免费大全6| 美女被艹到高潮喷水动态| 久久久久久伊人网av| 91精品国产九色| 日韩av在线免费看完整版不卡| 天天躁日日操中文字幕| 18+在线观看网站| 日韩伦理黄色片| 麻豆精品久久久久久蜜桃| 91久久精品电影网| 日韩欧美精品免费久久| www.av在线官网国产| 一级毛片久久久久久久久女| 日本一本二区三区精品| 亚洲欧美一区二区三区黑人 | 亚洲,一卡二卡三卡| 亚洲国产精品国产精品| 日本与韩国留学比较| 成年人午夜在线观看视频| 欧美最新免费一区二区三区| 久久97久久精品| 内射极品少妇av片p| 日本猛色少妇xxxxx猛交久久| 男人添女人高潮全过程视频| 国产精品一及| 亚洲成人久久爱视频| 日韩av不卡免费在线播放| 久久国内精品自在自线图片| h日本视频在线播放| 一级二级三级毛片免费看| 赤兔流量卡办理| 日韩av免费高清视频| 精品人妻一区二区三区麻豆| 亚洲精品成人久久久久久| 热re99久久精品国产66热6| 日韩av不卡免费在线播放| 老司机影院成人| 熟女电影av网| 91精品伊人久久大香线蕉| 久久亚洲国产成人精品v| xxx大片免费视频| 老师上课跳d突然被开到最大视频| 美女主播在线视频| 久久人人爽人人片av| 亚洲av欧美aⅴ国产| 婷婷色综合www| 国产老妇女一区| 中文欧美无线码| 免费在线观看成人毛片| 夜夜看夜夜爽夜夜摸| 中国国产av一级| 高清毛片免费看| 人妻制服诱惑在线中文字幕| 2021天堂中文幕一二区在线观| 精品午夜福利在线看| 91午夜精品亚洲一区二区三区| 日韩大片免费观看网站| 国产亚洲av嫩草精品影院| 精品人妻一区二区三区麻豆| 国产老妇女一区| 丰满人妻一区二区三区视频av| 一区二区三区免费毛片| 国产成人a区在线观看| 亚洲精品视频女| 一个人观看的视频www高清免费观看| 免费少妇av软件| 精华霜和精华液先用哪个| 国产高清国产精品国产三级 | 色综合色国产| 一级毛片久久久久久久久女| 精品视频人人做人人爽| xxx大片免费视频| 夜夜爽夜夜爽视频| 亚洲av福利一区| 国产成人午夜福利电影在线观看| 大片电影免费在线观看免费| 免费少妇av软件| 久久人人爽人人爽人人片va| 新久久久久国产一级毛片| 国产精品av视频在线免费观看| 99热这里只有精品一区| 黄片wwwwww| 久久久久久久大尺度免费视频| 久久精品国产自在天天线| 国产v大片淫在线免费观看| 精华霜和精华液先用哪个| 亚洲国产成人一精品久久久| 在线观看一区二区三区| 蜜桃久久精品国产亚洲av| 久久久久精品久久久久真实原创|