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

    河流過程對小型集水區(qū)沉積物地球化學(xué)與重礦物組成的影響

    2025-03-04 00:00:00趙慧瑩謝遠(yuǎn)云遲云平康春國吳鵬孫磊魏振宇
    沉積學(xué)報(bào) 2025年1期
    關(guān)鍵詞:特征

    摘 要 【目的】河流沉積物記錄了豐富的源區(qū)信息,對于探究流域源區(qū)風(fēng)化、搬運(yùn)—沉積循環(huán)及母巖成分具有重要意義。與整合了不同風(fēng)化機(jī)制并包含復(fù)雜巖石類型的大型流域不同,小型河流具有相似的氣候和巖性背景,對研究沉積物物質(zhì)組成有重要研究價(jià)值?!痉椒ā繛榱颂骄亢恿鬟^程對沉積物物質(zhì)組成的影響,將從巴蘭河上游到下游的邊灘上獲取的8個(gè)樣品進(jìn)行分粒級(jí)處理(lt;63 μm、63~125 μm、125~250 μm),將小于63 μm的樣品進(jìn)行地球化學(xué)測試(常量、微量、稀土和Sr-Nd同位素),63~125 μm和125~250 μm的樣品進(jìn)行重礦物測試?!窘Y(jié)果】巴蘭河沉積物經(jīng)歷了較低的化學(xué)風(fēng)化作用、成熟度較低,處于初次循環(huán)沉積,物源主要來源于長英質(zhì)母巖。物源指標(biāo)(La/Sc、Co/Th、Cr/Th、Th/Sc)和Sr-Nd同位素組成總體變化較小,能夠代表巴蘭河總體流域特征?!窘Y(jié)論】巴蘭河樣品呈現(xiàn)出一致的化學(xué)風(fēng)化程度、成熟度、再循環(huán)特征、母巖特征,這說明河流過程并沒有對巴蘭河沉積物的地球化學(xué)組成造成明顯影響。重礦物經(jīng)過河流的搬運(yùn)后稀釋混合,不同河段及不同粒級(jí)組分的重礦物明顯呈現(xiàn)不同的含量組成,表明短距離的搬運(yùn)很難使重礦物均勻分配到不同粒級(jí)組分。因此,在進(jìn)行河流源—匯系統(tǒng)研究時(shí),少數(shù)的樣品很難代表整個(gè)流域的重礦物組成,即便是在小型河流中。另外,強(qiáng)調(diào)了Sr-Nd同位素與不活動(dòng)元素比值結(jié)合可以有效進(jìn)行河流物源示蹤。

    關(guān)鍵詞 巴蘭河;河流過程;地球化學(xué)組成;重礦物;Sr-Nd同位素

    第一作者簡介 趙慧瑩,女,1998年出生,碩士研究生,第四紀(jì)地質(zhì)與環(huán)境變化,E-mail: 18746917868@163.com

    通信作者 謝遠(yuǎn)云,男,教授,E-mail: xyy0451@hrbnu.edu.cn

    中圖分類號(hào) P595 文獻(xiàn)標(biāo)志碼 A

    0 引言

    河流作為陸源碎屑物向沉積區(qū)搬運(yùn)的主要通道,是陸源物質(zhì)遷移的重要途徑。河流沉積物記錄了豐富的地形、古氣候、古環(huán)境以及構(gòu)造活動(dòng)和地殼演化等多種信息,能很好地反映流域物質(zhì)組成,是沉積物源—匯系統(tǒng)研究的重要載體[1?4]。然而沉積物在搬運(yùn)和沉積過程中常發(fā)生分選作用,具有相同沉降速度的顆粒沉積在一起。沉積物中的礦物因其粒度、形狀和密度的差異而被分選,這在很大程度上導(dǎo)致沉積物礦物組成和化學(xué)成分發(fā)生了改變[5?7]。因此,對河流沉積物的地球化學(xué)元素和重礦物組成進(jìn)行研究,對了解河流地質(zhì)過程具有重要意義。

    大部分對河流物源分析的工作都是針對流域面積廣闊、支流水系發(fā)達(dá)、環(huán)境較為復(fù)雜的大河流域,這些大河甚至貫穿多個(gè)地貌單元和構(gòu)造帶。例如,長江作為鏈接亞洲季風(fēng)系統(tǒng)以及青藏高原構(gòu)造運(yùn)動(dòng)的紐帶,流經(jīng)了廣闊的大陸克拉通區(qū)域,其沉積物不僅記錄了沉積區(qū)域的環(huán)境信息,還記錄了流域多個(gè)地貌體系和構(gòu)造單元的演化信息[8?11];黃河的形成演化則對理解東亞地貌格局發(fā)育歷史和構(gòu)造環(huán)境變化至關(guān)重要[12?14]。然而,與整合了不同風(fēng)化機(jī)制并包含復(fù)雜巖石類型的大型河流流域不同,小型集水區(qū)具有相似的氣候和巖性背景,這使研究的復(fù)雜性大大降低[6]。因此,這些小河可被視為研究河流過程對沉積物物質(zhì)組成影響的天然的沉積學(xué)實(shí)驗(yàn)室。

    松花江作為中國七大水系之一,是東北地區(qū)最重要的河流。20世紀(jì)70—90年代一些學(xué)者對其開展了深入研究[15],近年來學(xué)者對松花江進(jìn)行了重礦物分析來探討松花江的水系演化模式[7,16?17]。巴蘭河作為松花江水系的支流,流域面積較小,長度較短,具有相對單一的氣候和巖性背景,為探究小型流域的河流過程提供了一個(gè)很好的機(jī)會(huì)。目前對于巴蘭河的研究僅局限于與松花江支流的對比分析[17?19],針對巴蘭河這類小型集水區(qū)的系統(tǒng)研究相對較少。本文以巴蘭河為研究對象,闡述其地球化學(xué)特征(常量、微量、稀土和Sr-Nd同位素)以及重礦物特征,揭示河流過程對小型集水區(qū)沉積物物質(zhì)組成產(chǎn)生的影響,完善了松花江水系的研究,也為今后探索其他小河流域奠定基礎(chǔ)。

    1 研究區(qū)概況

    松花江水系位于松嫩平原東北部,是中國七大河流之一,主要流經(jīng)黑龍江、吉林兩省,流域面積為55.68×104 km2。松花江有南北兩源,北源為發(fā)源于大興安嶺支脈伊勒呼里山的嫩江,南源為發(fā)源于吉林省長白山天池的第二松花江,兩江于三岔河匯合后稱松花江[20]。根據(jù)松花江的河道特性以及流域地形特征,可將其分為上、中、下三段。上游為松原市至哈爾濱市,全長240 km;中游為哈爾濱市至佳木斯市,全長432 km;下游為佳木斯市至同江市,全長267 km[17,20]。根據(jù)物源對松花江水系的貢獻(xiàn)不同,可以劃分為依蘭方向河流和松原方向河流。依蘭方向河流物源主要來自小興安嶺和張廣才嶺,松原方向河流物源主要來自大興安嶺和長白山西側(cè)[17]。

    巴蘭河是松花江左岸支流(圖1),位于哈爾濱市東北部依蘭縣(45°50′40″~46°39′20″ N,129°11′50″~130°11′40″ E)內(nèi)。全長108 km,河寬20~60 m,水深0.3~1.0 m,流域總面積為2 075 km2,每年11月中旬至次年4月上旬為結(jié)冰期。巴蘭河發(fā)源于小興安嶺余脈五華頂子北側(cè),流經(jīng)依蘭縣并注入松花江。依蘭縣河網(wǎng)密布,河流較多,年平均氣溫3.1 ℃~3.6 ℃,松花江以北地區(qū)較以南地區(qū)低0.5 ℃。

    2 材料與方法

    在巴蘭河的邊灘,從上游到下游采集8個(gè)樣品(圖1),依次命名為BLH-8~BLH-1。采樣位置盡量避開城市和可能污染源,剝?nèi)ケ砻?0 cm進(jìn)行取樣。樣品在室內(nèi)風(fēng)干后,使用200目(63 μm)標(biāo)準(zhǔn)分樣篩,獲得小于63 μm的粒度組分用于地球化學(xué)分析。使用荷蘭帕納(PANalytica)XRF(X-ray fluorescence)光譜儀進(jìn)行常量元素測試,測量誤差控制在3%以內(nèi)。使用電感耦合等離子體質(zhì)譜儀(ICP-MS)進(jìn)行微量、稀土元素的測試。為了保證實(shí)驗(yàn)結(jié)果的精準(zhǔn)度,測量過程中添加國際標(biāo)樣、空白樣品及重復(fù)樣進(jìn)行分析[21],測試準(zhǔn)確度、精度和不確定度(相對偏差)均小于2%[22?23]。采用標(biāo)準(zhǔn)離子交換方法進(jìn)行Sr和Nd的分離,并分別使用86Sr/88Sr=0.119 4 和146Nd/144Nd=0.721 9,將Sr和Nd同位素比值標(biāo)準(zhǔn)化以校正質(zhì)量分?jǐn)?shù)。通過定期使用國際標(biāo)準(zhǔn)樣品NBS 987和JMC進(jìn)行儀器的準(zhǔn)確度測試,Sr、Nd同位素測試平均值分別為87Sr/86Sr=0.710 250±7(2σ)和143Nd/144Nd=0.512 109±3 (2σ)。εNd(0)值(Nd同位素組成)采用球粒隕石均一源(CHUR)現(xiàn)代值143Nd/144Nd=0.512 638 計(jì)算。整個(gè)流程的化學(xué)分析空白為:Srlt;1 ng,Ndlt;50 pg。元素地球化學(xué)測試和Sr?Nd 同位素的測定在蘭州大學(xué)進(jìn)行。

    此外,這8 個(gè)樣品被干篩分成63~125 μm 和125~250 μm兩個(gè)粒度組分用于重礦物分析。首先將重礦物樣品烘干并稱重后,利用三溴甲烷(密度為2.89 g/cm3)分離重礦物樣品,用酒精反復(fù)沖洗分離得到的樣品,烘干后(60 ℃恒溫)再次稱重[21,24?25],最后用薄紙包住永久磁鐵反復(fù)磁吸,得到重礦物的無磁、電磁、強(qiáng)磁部分,將分離出來的重礦物部分放在雙目鏡下進(jìn)行鑒定,每個(gè)重礦物樣品鑒定顆粒數(shù)在900粒以上[23,26]。重礦物鑒定在廊坊誠信地質(zhì)公司進(jìn)行。

    3 結(jié)果

    3.1 地球化學(xué)元素特征

    巴蘭河沉積物的化學(xué)成分以SiO2(69.40%)、Al2O3(14.64%)和Fe2O3(3.97%)為主,相較于UCC(上地殼),MnO和TiO2輕微富集,MgO、CaO和P2O5呈現(xiàn)虧損狀態(tài),F(xiàn)e2O3 和Na2O 呈現(xiàn)出較大可變性,SiO2、Al2O3、K2O 與UCC 含量相似(圖2)。CaO 的含量在BLH-8樣品以及BLH-6樣品出現(xiàn)高值;Na2O在BLH-7樣品以及BLH-3樣品出現(xiàn)高值;Fe2O3在BLH-2樣品出現(xiàn)高值;其余元素樣品間差異較小。

    巴蘭河沉積物微量元素與UCC相比(圖2),微量元素中Zr、Hf呈現(xiàn)富集狀態(tài),Y、Th、U部分富集,其余元素全部虧損。BLH-2樣品的Sc、V、Co、Ni、Zn、Nb值顯著大于其他樣品。Zr在BLH-1樣品和BLH-2樣品的含量顯著高于其他樣品。Rb在BLH-7樣品含量高于其他樣品。Y、Th、U、Hf在BLH-2樣品、BLH-7樣品中含量高于其他樣品,其余微量元素差異不明顯。

    巴蘭河沉積物稀土總量(ΣREE)為101.15×10-6,低于UCC(平均值146.37×10-6)和后太古宙頁巖標(biāo)準(zhǔn)值(PAAS:184.77×10-6)[27];輕稀土(LREE)富集,平均值為88.45×10-6;重稀土(HREE)豐度較低,平均值為12.70×10-6(圖2)。輕稀土與重稀土比值為5.15~8.78,平均為6.97;Eu/Eu* =0.46×10-6,Ce/Ce* =0.93×10-6,表明巴蘭河存在Eu、Ce負(fù)異常。(La/Sm)N 值介于3.26~3.78,平均為3.52,表明輕稀土間分餾程度較強(qiáng);(Gd/Yb)N 值介于0.87~1.39,平均為1.19,表明重稀土間分餾程度較弱;(La/Yb)N 值介于3.74~8.19,平均為6.16,稀土配分曲線為右傾型,表明輕稀土富集,重稀土虧損。就稀土元素而言,BLH-7 樣品、BLH-6樣品與BLH-2樣品的值明顯高于其他樣品。

    巴蘭河的87Sr/86Sr比值比較穩(wěn)定,介于0.705 130~0.714 083,平均為0.711 931,與松花江依蘭方向支流樣品的87Sr/86Sr 比值相吻合[19]。143Nd/144Nd 值介于0.512 335~0.512 402,平均為0.512 367。εNd(0)值介于-5.914~-4.604,平均為-5.277(表1)。

    3.2 重礦物組合特征

    利用沉積物的重礦物類型和組成來判斷源區(qū)是探討河流過程的重要手段[25,28]。在河流沉積物中,由于重礦物具有化學(xué)性質(zhì)穩(wěn)定、耐磨蝕等特性,因此在風(fēng)化、搬運(yùn)以及成巖過程中受到的影響相對較弱,能夠較多地保留母巖信息[29]。巴蘭河重礦物分為兩個(gè)粒度組分(圖3):63~125 μm、125~250 μm,檢查出的重礦物有16種,其中菱鐵礦僅在個(gè)別樣品中出現(xiàn),為偶見礦物。

    巴蘭河(63~125 μm)的重礦物以角閃石(22.12%~55.06%)、鐵磁性礦物(6.17%~22.16%)和榍石(10.09%~18.37%)為主,占72.96%;鈦鐵礦(3.60%~18.32%)、綠簾石(1.72%~11.30%)、鋯石(0.45%~5.01%)、石榴子石(0.33%~4.27%)、赤褐鐵礦(1.33%~3.96%)次之,占20.57%;其他礦物占6.47%。在125~250 μm組分中(圖3),巴蘭河的重礦物以角閃石(36.21%~53.39%)、榍石(12.04%~35.97%)和綠簾石(2.68%~26.21%)為主,占75.96%;鈦鐵礦(0.21%~14.25%)、鐵磁性礦物(1.08%~17.50%)、赤褐鐵礦(0.27%~3.52%)、石榴子石(0.48%~2.62%)、白鈦石(0.10%~4.83%)次之,占16.79%;其他重礦物均不足1%。

    4 討論

    沉積物的地球化學(xué)組成主要受物源區(qū)母巖性質(zhì)和構(gòu)造背景的影響,但是河流沉積物在搬運(yùn)過程中還會(huì)受到化學(xué)風(fēng)化、搬運(yùn)方式和介質(zhì)、遷移距離、分選作用、再循環(huán)等影響,給利用沉積物的地球化學(xué)特征進(jìn)行物源分析及成巖構(gòu)造背景判別帶來多解性與復(fù)雜性[30]。為了評估小型河流過程是否會(huì)對來自同一物源區(qū)的沉積物地球化學(xué)與重礦物組成造成影響,為此我們從沉積物的化學(xué)風(fēng)化特征、成熟度、再循環(huán)特征以及母巖特征來探究河流過程對地球化學(xué)組成的影響,結(jié)合河流的分選作用探究河流過程對沉積重礦物組成的影響。

    4.1 沉積物化學(xué)風(fēng)化特征與成熟度

    河流沉積物的化學(xué)組成可以敏感地反映流域化學(xué)風(fēng)化程度,而沉積物風(fēng)化指數(shù)不僅會(huì)受到風(fēng)化強(qiáng)度的影響,還會(huì)受到沉積物侵蝕、搬運(yùn)和沉積過程中水力過程的干擾[5,31?32],這些風(fēng)化指數(shù)不僅記錄了沉積物搬運(yùn)和沉積過程中的水力分選效應(yīng),還反映了流域完整的風(fēng)化歷史[6,33]。

    Nesbitt et al.[34]為了評價(jià)沉積物的風(fēng)化程度,提出了化學(xué)蝕變指數(shù)(CIA):

    CIA = 100 × [Al2O3/(Al2O3 + CaO* + Na2O +K2O) ] (1)

    式中:各氧化物均為氧化物分子摩爾數(shù),其中CaO*為不包含碳酸鹽和磷酸鹽內(nèi)CaO的硅酸鹽礦物中的摩爾含量。

    CIA的值越大,沉積物的化學(xué)風(fēng)化程度越強(qiáng),CIA介于50~60代表初等程度的化學(xué)風(fēng)化,CIA介于60~80代表中等程度的化學(xué)風(fēng)化,CIA介于80~100代表強(qiáng)烈的化學(xué)風(fēng)化程度。巴蘭河沉積物的CIA值介于42.77~58.19(圖4),平均為51.16,說明巴蘭河經(jīng)歷了程度較弱的化學(xué)風(fēng)化。

    WIP是Parker[35]最先提出的評價(jià)硅酸鹽巖風(fēng)化強(qiáng)度的指標(biāo),通過計(jì)算堿金屬和堿土金屬元素在風(fēng)化產(chǎn)物中所占的比例來判斷化學(xué)風(fēng)化程度。

    WIP = 100 × (CaO*/0.7 + 2Na2O/0.35 +2K2O/0.25 + MgO/0.9) (2)

    WIP 值越小,風(fēng)化作用越強(qiáng)[33,36?37]。巴蘭河的WIP值介于54.00~99.94(圖5),平均為71.84,反映了巴蘭河風(fēng)化程度較弱。

    在A-CN-K圖中(圖4),根據(jù)元素的遷移順序,風(fēng)化作用可分為以下三個(gè)階段:早期斜長石風(fēng)化,表現(xiàn)為風(fēng)化趨勢線平行于A-CN線指向A;中期鉀長石和伊利石風(fēng)化,斜長石消失,表現(xiàn)為風(fēng)化趨勢線平行于A-K線;晚期殘余風(fēng)化產(chǎn)物(高嶺石、三水鋁石以及綠泥石),表現(xiàn)為產(chǎn)物集中分布在A點(diǎn)[33,38]。巴蘭河樣品分布在斜長石與鉀長石附近,風(fēng)化趨勢線平行于A-CN線指向A點(diǎn),樣品點(diǎn)全部低于PAAS,說明樣品經(jīng)歷了較低的化學(xué)風(fēng)化作用。

    氧化物比值可以用來表征樣品之間的差異,黏土礦物和非黏土硅酸鹽礦物的特點(diǎn)是氧化鋁的比例差別很大。Cox et al. [37]定義成分可變性指數(shù)ICV(Index of Compositional Variability),用于測量巖石或礦物中氧化鋁相對于其他主要陽離子的豐度。ICV也可用于泥巖作為成分成熟度的測量,成分成熟的泥巖缺乏非鈣硅酸鹽或以礦物(如kandite族的礦物)為主,將具有低值。成分不成熟的泥巖往往出現(xiàn)在構(gòu)造活躍的環(huán)境中,將具有高值,這里的Fe2O3為總鐵含量[33,37,39?40]。

    ICV = (CaO + K2O + Na2O + Fe2O3 + MgO +TiO2 + MnO)/Al2O3 (3)

    成分不成熟的沉積物具有高ICV值(gt;1),而成分成熟的沉積物具有低ICV值(lt;1)。巴蘭河的ICV均值為1.08,說明樣品具有高的ICV值,沉積物成熟度低。

    SiO2/Al2O3、Al2O3/Na2O和K2O/Na2O比值亦是判斷沉積物成熟度的指標(biāo),高比值(分別大于6、大于5和大于1)表明沉積物成熟度高[37]。較低的K2O/Na2O比值表明,源區(qū)主要是演化程度較低的富含斜長石的巖石[41?42]。此外,碎屑沉積物中的K2O/Al2O3比值可以指示堿性長石相對于鋁硅酸鹽的豐度。K2O/Al2O3大于0.5表明有大量的堿性長石,而小于0.5表明堿性長石相對于鋁硅酸鹽有所減少[37]。巴蘭河沉積物SiO2/Al2O3 值介于4.51~5.01,平均為4.75;Al2O3/Na2O值介于1.83~5.10,平均為3.56;K2O/Na2O 值介于0.42~0.95,平均為0.73;三者的低比值都表明巴蘭河沉積物的成分成熟度較低。

    綜合巴蘭河沉積物的CIA、WIP、A-CN-K圖、ICV以及SiO2/Al2O3、Al2O3/Na2O和K2O/Na2O比值,表明巴蘭河沉積物經(jīng)歷了較低的化學(xué)風(fēng)化作用,成熟度較低。從巴蘭河沉積物的化學(xué)風(fēng)化特征以及成熟度的指標(biāo)上看,這8個(gè)樣品并沒有明顯變化,表明河流過程對巴蘭河沉積物本身的化學(xué)風(fēng)化特征以及成熟度影響較小。

    4.2 沉積物再循環(huán)特征

    在化學(xué)風(fēng)化水平高的地區(qū),初次循環(huán)下的沉積物CIA/WIP值一般小于10,而多次循環(huán)下的沉積物CIA/WIP比值一般大于10,甚至可以達(dá)到100[5]。巴蘭河的CIA/WIP 值介于0.43~1.08(圖5),平均為0.76,所有值均小于10,遠(yuǎn)低于初次與多次循環(huán)的臨界點(diǎn)(CIA/WIP=10),說明巴蘭河樣品處于初次循環(huán)。

    Ohta[43]提出,SiO2/Al2O3 和Na2O/K2O 比值在分選過程中對石英、鈉長石和伊利石的分餾非常敏感。SiO2/Al2O3與Na2O/K2O二元圖解(圖5)可通過沉積分選作用對沉積物化學(xué)組成的影響來判斷沉積物的來源及沉積再循環(huán)特征。該二元圖解可被用來度量區(qū)分首次循環(huán)和再循環(huán)沉積:直接來源于巖漿巖源區(qū)的沉積物(初次循環(huán))的化學(xué)組成形成一條垂直的線性趨勢線;而直接來源于沉積巖源區(qū)的沉積物(再循環(huán)沉積)的化學(xué)組成會(huì)形成一條水平的線性趨勢線。從巴蘭河樣品的二元圖解來看,樣品分布呈現(xiàn)垂直狀態(tài),說明巴蘭河樣品處于初次循環(huán)。

    綜合巴蘭河沉積物的CIA/WIP比值、SiO2/Al2O3與Na2O/K2O二元圖解,巴蘭河沉積物的再循環(huán)特征指標(biāo)值變化較小,均指示為初次循環(huán)沉積,河流過程對巴蘭河沉積物的再循環(huán)特征影響較小。

    4.3 母巖性質(zhì)

    巴蘭河作為松花江的支流,對于它的母巖性質(zhì)研究僅局限于從重礦物方向闡述,針對其地球化學(xué)方面的研究尚未涉及。由巴蘭河的地理位置可知其為依蘭方向河流,物源主要來自小興安嶺東南側(cè)[7,17,44]。而小興安嶺的巖性以花崗巖、花崗閃長巖和變質(zhì)巖為特征[44]。巴蘭河作為一條氣候與巖性背景相對單一的小型河流,8個(gè)樣品能否準(zhǔn)確判別其母巖類型,能否對應(yīng)實(shí)際的地質(zhì)情況,河流過程是否會(huì)改變沉積物攜帶的母巖性質(zhì),其結(jié)果能為未來探究小型河流源區(qū)時(shí)提供一定的依據(jù)。

    不活動(dòng)元素及其比值在長英質(zhì)和鐵鎂質(zhì)源巖中存在顯著差異,常被用來示蹤物源[23,45?47]。相容元素和不相容元素在母巖中的相對比值可以確定母巖性質(zhì):過渡元素(Ni、Sc、V、Co、Cr)等富集在基性母巖(鐵鎂質(zhì))中[46];而高場強(qiáng)元素(Y、Zr、Hf、Nb、Ta、Th、U)為不相容元素,主要富集在酸性母巖(長英質(zhì))中[48]。因此,元素比率和基于這些元素的圖解對于區(qū)分沉積物來源是非常有效的。在Co/Th-La/Sc圖解中(圖6a)[49],樣品投點(diǎn)在長英質(zhì)火山巖附近,且樣品的Co/Th小于1;在Cr/Th-Th/Sc二元圖解中(圖6b),樣品的Cr/Th值和Th/Sc值,分別介于0.73~2.74(平均為2.06)、0.64~2.22(平均為1.20),從圖上可看出巴蘭河沉積物主要在長英質(zhì)火山巖附近。

    在TiO2-Zr雙變量圖中(圖6c)[50],樣品點(diǎn)全部分布在長英質(zhì)火成巖區(qū)域,表明所研究的巴蘭河沉積物物源巖具有長英質(zhì)火成巖性質(zhì)。Roser et al.[51]提出了一個(gè)用7 種常量元素組成兩組判別函數(shù)(DF1、DF2)用來確定4種不同來源成分的物源鑒別圖(圖6d),在此圖中巴蘭河沉積物全部分布在長英質(zhì)火成巖物源區(qū)。

    在Co/Th-La/Sc、Cr/Th-Th/Sc、TiO2-Zr二元圖圖解以及判別函數(shù)圖上,樣品幾乎全部落入長英質(zhì)火成巖源區(qū)范圍內(nèi),指示巴蘭河沉積物物源主要來源于長英質(zhì)火成巖。將巴蘭河沉積物樣品的母巖特征與小興安嶺的巖性[7]進(jìn)行對比,其母巖特征與現(xiàn)今物源區(qū)出露的巖石特征基本一致。8個(gè)樣品足以判別其母巖類型,河流過程并沒有改變巴蘭河沉積物攜帶的母巖性質(zhì),能夠很好地反映物源。

    4.4 河流地質(zhì)過程對地球化學(xué)組成的影響

    Sr-Nd同位素比值在示蹤沉積物的物源和搬運(yùn)方面具有巨大的潛力,Sr同位素分餾對風(fēng)化條件的變化和遷移過程的差異很敏感,一般認(rèn)為會(huì)受到化學(xué)風(fēng)化和粒度效應(yīng)的影響。而Nd同位素比值僅與源巖的性質(zhì)有關(guān),基本不受粒度分異和化學(xué)風(fēng)化的影響。因此,Sr-Nd同位素比值變化幅度的大小能夠作為河流過程是否影響沉積物地球化學(xué)組成的判斷依據(jù)之一[52]。

    87Sr/86Sr 和143Nd/144Nd 比值被廣泛用于地球化學(xué)示蹤[24,53],87Sr單方向增長并以86Sr作為比較基礎(chǔ),根據(jù)等時(shí)線外推或含Sr礦物測定得到地質(zhì)體形成時(shí)的初始87Sr/86Sr來示蹤其物質(zhì)來源[54]。Nd同位素比較穩(wěn)定,不易受到風(fēng)化、搬運(yùn)、堆積以及后期成巖影響,能真實(shí)反映母巖特征。Nd同位素組成的差異是由于不同環(huán)境中其放射性母體含量的差異長期積累或演化造成的。εNd值指Nd同位素的初始值相對于球粒隕石質(zhì)地幔演化值的偏差[54]。

    從不同樣品之間的87Sr/86Sr比值差異來看(表1),只有樣品之間的差異大于0.000 050時(shí),才被認(rèn)為不同樣品之間的絕對Sr同位素差異具有地質(zhì)意義。也有不同的觀點(diǎn)認(rèn)為,如果樣品之間的87Sr/86Sr比值的差異在0.001 500 以內(nèi),則Sr 同位素組成是相同的[52,55?56]。巴蘭河的87Sr/86Sr 值從0.705 160 614 到0.714 082 914,最大差值為0.008 952 299,最小差值為0.000 400 076。BLH-2樣品與其他樣品的Sr具有明顯差異。相比之下,BLH-8樣品和BLH-7樣品只能與其余樣品發(fā)生部分差異,其余樣品間差異不大,說明其余樣品之間不具有地質(zhì)意義上的差異。總體來說,Sr同位素僅在個(gè)別樣品間產(chǎn)生差異,大部分樣品之間不具有地質(zhì)意義上的差異。

    從Nd同位素的變化范圍來看(表1),143Nd/144Nd值的變化不大,最低值為4.6,最高值為5.9。對于143Nd/144Nd之間εNd(0)值差異來看,目前認(rèn)為只有當(dāng)樣品間εNd(0)值的差異大于1ε 單位時(shí),Nd同位素組成才會(huì)發(fā)生本質(zhì)變化或出現(xiàn)同位素異常[57],這種Nd同位素變化才是由地質(zhì)作用而不是實(shí)驗(yàn)誤差觸發(fā)的[53]。從上游到下游的樣品中,樣品之間εNd(0)最大差值為1.31ε,最小差值為0.03ε,僅有兩個(gè)樣品超過1ε。樣品間εNd(0)差值較小,變化不明顯。所以可以認(rèn)為整條巴蘭河的Nd同位素很穩(wěn)定,Nd同位素不易受到風(fēng)化搬運(yùn)堆積以及后期成巖影響,能真實(shí)反映流域母巖特征,8個(gè)樣品可以代表巴蘭河的Nd同位素特征。

    河流過程對地球化學(xué)組成的影響體現(xiàn)在多個(gè)方面,包括化學(xué)風(fēng)化特征、再循環(huán)特征以及母巖性質(zhì)特征。在巴蘭河樣品中,常量元素、微量元素變化不明顯,不活動(dòng)元素比值La/Sc、Co/Th、Cr/Th、Th/Sc(圖6)和Sr-Nd同位素(表1)總體上變化較小,能夠很好地示蹤物源。樣品呈現(xiàn)出一致的化學(xué)風(fēng)化程度、成熟度、再循環(huán)特征、母巖特征,這說明河流過程對地球化學(xué)組成的影響較小。雖然物源示蹤指標(biāo)ΣLREE、ΣHREE、ΣREE、LREE/HREE 以及化學(xué)風(fēng)化指標(biāo)CIA、CIA/WIP在BLH-2樣品出現(xiàn)最高值,但附近的BLH-1、BLH-3號(hào)樣品卻沒有出現(xiàn)BLH-2樣品的高值現(xiàn)象,推測可能BLH-2樣品的稀土元素與風(fēng)化指標(biāo)的高值受到了人為因素的影響。

    巴蘭河的8個(gè)樣品地球化學(xué)組成變化不大,雖然BLH-2樣品的稀土元素與風(fēng)化指標(biāo)出現(xiàn)高值,但并不影響最后的物源判別,單一樣品能夠代表巴蘭河總體的地球化學(xué)特征,但在選擇樣品數(shù)量時(shí)也要考慮人為因素的影響。對于流域面積較小、氣候相似、巖性背景單一的巴蘭河而言,河流過程沒有對巴蘭河沉積物的地球化學(xué)組成造成明顯影響,可以準(zhǔn)確恢復(fù)源區(qū)母巖性質(zhì)。

    4.5 河流地質(zhì)過程對重礦物組成的影響

    從巴蘭河上游到下游,河床寬度越來越寬,河床沙出現(xiàn)的規(guī)模也越來越大,河床沙顆粒的磨圓度也越來越好。在63~125 μm組分中,角閃石、鐵磁性礦物、榍石、鈦鐵礦富集,重礦物組合為角閃石+鐵磁性礦物+榍石+鈦鐵礦+綠簾石+鋯石+石榴子石,金紅石、電氣石、菱鐵礦僅在個(gè)別樣品出現(xiàn)。在125~250μm組分中,角閃石、榍石、綠簾石、鈦鐵礦富集,重礦物組合為角閃石+榍石+綠簾石+鈦鐵礦+鐵磁性礦物+赤褐鐵礦+石榴子石,鋯石、銳鈦礦、金紅石、電氣石、菱鐵礦、磁鐵礦僅在個(gè)別樣品中出現(xiàn)。不同粒徑組分重礦物組成不同,對比兩個(gè)粒度組分的重礦物含量可以看出,榍石、石榴子石、綠簾石、鈦鐵礦、磁鐵礦存在明顯差別,但兩粒級(jí)組分的重礦物組合類似。

    同一粒度組分不同河段的重礦物含量則波動(dòng)變化(圖3),呈現(xiàn)出明顯差異,特別是金紅石、電氣石只在下游穩(wěn)定出現(xiàn),菱鐵礦僅在個(gè)別樣品(BLH-2)出現(xiàn)。在物源相同的條件下,巴蘭河樣品表現(xiàn)出不同的礦物組成。同一粒徑不同河段出現(xiàn)明顯的重礦物缺失(金紅石、電氣石),同一河段不同粒徑也出現(xiàn)重礦物缺失的現(xiàn)象(鋯石、金紅石、銳鈦礦、磁鐵礦)。當(dāng)某些重礦物只局限分布在某一樣品、某一粒度組分中,單一窗口或者過少樣品都很難準(zhǔn)確得到源區(qū)母巖重礦物組成。因此,為了準(zhǔn)確獲得源區(qū)母巖信息,樣品的數(shù)量非常重要。對于巴蘭河的8個(gè)重礦物樣品而言,單一樣品或者過少樣品都很難代表整條河流的重礦物特征,河流在搬運(yùn)過程中由于水動(dòng)力分選作用對重礦物組成造成較大影響。那么,穩(wěn)定重礦物比值能否消除重礦物粒度效應(yīng)的影響,這對于進(jìn)一步探討小型流域河流過程對重礦物的影響提供幫助。

    重礦物組合在河流搬運(yùn)以及沉積過程中會(huì)受到諸如機(jī)械磨損、化學(xué)溶蝕和水動(dòng)力分選等因素的影響,導(dǎo)致物源判斷準(zhǔn)確性不足[17,58]。在相似的水動(dòng)力條件下,穩(wěn)定重礦物比值能夠在物源示蹤過程中最大限度地消除礦物的粒度效應(yīng)的影響[17],較好地反映沉積物的成熟度并反映物源特征[58?59],這些比值被稱為重礦物特征指數(shù)。如:ZTR指數(shù)、ATi指數(shù)、GZi指數(shù)。公式如下:

    ZTR指數(shù)= 鋯石% + 電氣石% + 金紅石% (4)

    ATi指數(shù)= 100 × 磷灰石%/磷灰石% + 電氣石% (5)

    GZi指數(shù)= 100 × 石榴子石%/石榴子石% + 鋯石% (6)

    隨著粒徑的增大,巴蘭河的ZTR指數(shù)值減小(除BLH-3 樣品)、ATi 值減小、Gzi 值增大(圖7)。63~125 μm 粒度組分的樣品ZTR 指數(shù)明顯高于125~250 μm粒度組分(除BLH-3樣品),這說明細(xì)顆粒組分的成熟度更高。河流分選作用同時(shí)作用于粗細(xì)顆粒,但細(xì)顆粒組分搬運(yùn)距離更遠(yuǎn),更易被風(fēng)化磨蝕,使細(xì)顆粒組分的成熟度更高[24]。而125~250 μm粒度組分的ATi指數(shù)相比于63~125 μm粒度組分較低的原因可能是由于水動(dòng)力分選作用,粗顆粒組分沉積速率較快,風(fēng)化程度弱,ATi指數(shù)較低[24]。石榴子石富集于125~250 μm組分,分選作用下的石榴子石優(yōu)先富集于粗顆粒組分,使GZi指數(shù)高于細(xì)顆粒組分。從巴蘭河沉積物的重礦物特征指數(shù)上看,其值明顯受到粒度效應(yīng)的影響,重礦物特征指數(shù)無法最小化或者消除由粒度效應(yīng)帶來的偏差,這說明河流沉積物從上游到下游的搬運(yùn)和沉積過程中,在水動(dòng)力的作用下進(jìn)行分選,粗、細(xì)顆粒發(fā)生分離,從而導(dǎo)致混合礦物分離進(jìn)入不同粒級(jí)組分,不同河段的重礦物組成存在顯著差異,很難根據(jù)較少樣品代表整條巴蘭河的重礦物特征[17]。兩個(gè)組分的重礦物呈現(xiàn)不同的含量組成,明顯體現(xiàn)河流作用對重礦物組成的影響。

    5 結(jié)論

    (1) 巴蘭河沉積物在河流搬運(yùn)過程中因礦物分選而發(fā)生了顯著變化。河流沉積物在搬運(yùn)和沉積過程中,由于水動(dòng)力分選作用,粗、細(xì)顆粒發(fā)生分離,從而導(dǎo)致混合礦物分離進(jìn)入不同粒級(jí)組分,不同河段不同組分的重礦物組成存在顯著差異,這種差異明顯體現(xiàn)了小型集水區(qū)河流搬運(yùn)作用對重礦物組成的影響。

    (2) 巴蘭河沉積物的87Sr/86Sr和143Nd/144Nd及不活動(dòng)元素比值變化較小,所有指標(biāo)均指示巴蘭河沉積物經(jīng)歷了較弱程度的化學(xué)風(fēng)化作用、成熟度較低,處于初次循環(huán),其物源主要來源于長英質(zhì)火成巖。河流過程沒有對巴蘭河沉積物的地球化學(xué)組成造成明顯影響,可以準(zhǔn)確恢復(fù)源區(qū)母巖性質(zhì)。

    (3) 從巴蘭河整體的重礦物特征以及地球化學(xué)組成來看,重礦物在不同粒度組分的變化較為明顯。對于大型河流而言,長距離搬運(yùn)不能在不同粒徑中均勻混合重礦物,而對于小型河流巴蘭河而言,短距離的搬運(yùn)更難均勻混合重礦物。在地球化學(xué)組成上,Sr-Nd同位素與不活動(dòng)元素比值結(jié)合能夠較好地反映巴蘭河的物源。因此,在研究小型河流時(shí),過少的樣品很難代表整個(gè)流域的重礦物組成。相比于重礦物而言,地球化學(xué)組成能夠較好地進(jìn)行物源示蹤并反映流域情況。

    致謝 對提供寶貴修改意見的評審專家,一并謹(jǐn)致謝忱。感謝編輯部老師的辛勤付出,使論文得以完善。地球化學(xué)組成以及Sr?Nd同位素組成得到蘭州大學(xué)甘肅省西部礦產(chǎn)資源重點(diǎn)實(shí)驗(yàn)室熊聰慧、閆曉麗和陳萬峰老師的支持,重礦物組成得到了河北廊坊誠信地質(zhì)公司李鵬、張?jiān)频ひ约皬埮遢胬蠋煹闹С?,研究生李思琪、孫楊、孫建華、李秋杭、侯心茹、符錦霞參與了論文寫作討論,在此一并表示感謝。

    參考文獻(xiàn)(References)

    [1] McLennan S M, Taylor S R, Kr?ner A. Geochemical evolution ofArchean shales from South Africa. I. The Swaziland and PongolaSupergroups[J]. Precambrian Research, 1983, 22(1/2): 93-124.

    [2] 李明月,孫學(xué)軍,Karki K,等. 喜馬拉雅山中段柯西河跨境流域河流沉積物的礦物和元素特征[J]. 礦物巖石地球化學(xué)通報(bào),2019,38(5):989-997.[Li Mingyue, Sun Xuejun, Karki K, et al.Mineral and elemental characteristics of bed sediments in thetrans-boundary Koshi River basin in the central Himalayas[J].Bulletin of Mineralogy, Petrology and Geochemistry, 2019, 38(5):989-997.]

    [3] 陳丹婷,彭渤,方小紅,等. 洞庭湖“四水”入湖河床沉積物主量元素地球化學(xué)特征及意義[J]. 第四紀(jì)研究,2021,41(5):1267-1280.[Chen Danting, Peng Bo, Fang Xiaohong, et al. Geochemistryof major elements in bed sediments from inlets of the fourrivers to Dongting Lake, China[J]. Quaternary Sciences, 2021, 41(5): 1267-1280.]

    [4] 楊允柳,何云龍,解習(xí)農(nóng),等. 湄公河三角洲第四系沉積物地球化學(xué)特征及其地質(zhì)意義[J]. 地球科學(xué),2022,47(3):1107-1121.[Yang Yunliu, He Yunlong, Xie Xinong, et al. Geochemicalcharacteristics and geological significance of Quaternary sedimentsin the Mekong delta[J]. Earth Science, 2022, 47(3): 1107-1121.]

    [5] Garzanti E, Padoan M, Setti M, et al. Weathering geochemistryand Sr-Nd fingerprints of equatorial Upper Nile and Congo muds[J]. Geochemistry, Geophysics, Geosystems, 2013, 14(2):292-316.

    [6] Guo Y L, Yang S Y, Su N, et al. Revisiting the effects of hydrodynamicsorting and sedimentary recycling on chemical weatheringindices[J]. Geochimica et Cosmochimica Acta, 2018, 227: 48-63.

    [7] Wu P, Xie Y Y, Kang C G, et al. Effects of provenance, transportprocesses and chemical weathering on heavy mineral composition:A case study from the Songhua River drainage, NE China[J]. Frontiers in Earth Science, 2022, 10: 839745.

    [8] Barbour G B. Physiographic history of the Yangtze[J]. The GeographicalJournal, 1936, 87(1): 17-32.

    [9] Zheng H B. Birth of the Yangtze River: Age and tectonicgeomorphicimplications[J]. National Science Review, 2015,2(4): 438-453.

    [10] Wei H H, Wang E, Wu G L, et al. No sedimentary records indicatingsoutherly flow of the paleo-upper Yangtze River from thefirst bend in southeastern Tibet[J]. Gondwana Research, 2016,32: 93-104.

    [11] 魏傳義,尹功明,劉春茹,等. 長江現(xiàn)代河流沉積物石英ESR信號(hào)強(qiáng)度空間分布特征及其物源示蹤意義[J]. 第四紀(jì)研究,2022,42(4):1168-1180. [Wei Chuanyi, Yin Gongming, LiuChunru, et al. Spatial diversity of quartz ESR signal intensity ofthe modern Yangtze River fluvial sediments and its implicationsfor sediment provenance[J]. Quaternary Sciences, 2022, 42(4):1168-1180.]

    [12] 李吉均,方小敏,馬海洲,等. 晚新生代黃河上游地貌演化與青藏高原隆起[J]. 中國科學(xué):地球科學(xué),1996,26(4):316-322.[Li Jijun, Fang Xiaomin, Ma Haizhou, et al. Geomorphologicaland environmental evolution in the upper reaches of the YellowRiver during the Late Cenozoic[J]. Science China: Earth Sciences,1996, 26(4): 316-322.]

    [13] 傅建利,張珂,馬占武,等. 中更新世晚期以來高階地發(fā)育與中游黃河貫通[J]. 地學(xué)前緣,2013,20(4):166-181.[Fu Jianli,Zhang Ke, Ma Zhanwu, et al. The terrace (T5 and T4) formationsince the late Middle Pleistocene and its implication in thethrough cutting of the middle reach of Yellow River[J]. EarthScience Frontiers, 2013, 20(4): 166-181.]

    [14] 趙希濤,楊艷,賈麗云,等. 論晚期共和古湖時(shí)代、演化過程及其與地殼運(yùn)動(dòng)和黃河發(fā)育的關(guān)系[J]. 地球?qū)W報(bào),2021,42(4):451-471.[Zhao Xitao, Yang Yan, Jia Liyun, et al. A discussionon the age and evolution process of the late Gonghe paleolakeand its relations with the crustal movement and the developmentof the Yellow River[J]. Acta Geoscientica Sinica, 2021, 42(4):451-471.]

    [15] 宋國利,劉釗,于桂云. 黑龍江省松花江流域主要巖石類型中若干元素的背景含量及其環(huán)境意義[J]. 哈爾濱師范大學(xué)自然科學(xué)學(xué)報(bào),1986(4):113-121.[Song Guoli, Liu Zhao, Yu Guiyun.Background contents and environmental significance ofsome elements in major rock types in Songhua River basin[J].Natural Science Journal of Harbin Normal University, 1986(4):113-121.]

    [16] Xie Y Y, Kang C G, Chi Y P, et al. Reversal of the middle-upperSonghua River in the late Early Pleistocene, northeast China[J].Geomorphology, 2020, 369: 107373.

    [17] 李思琪,謝遠(yuǎn)云,康春國,等. 物源—河流過程—化學(xué)風(fēng)化對松花江水系沉積物重礦物組成的影響[J]. 地質(zhì)科學(xué),2022,57(1):207-229.[Li Siqi, Xie Yuanyun, Kang Chunguo, et al. Influenceof provenance-river process-chemical weathering onheavy mineral composition of the Songhua River sediment[J].Chinese Journal of Geology, 2022, 57(1): 207-229.]

    [18] 吳鵬,謝遠(yuǎn)云,康春國,等. 早更新世晚期松花江水系襲奪:地球化學(xué)和沉積學(xué)記錄[J]. 地質(zhì)學(xué)報(bào),2020,94(10):3144-3160.[Wu Peng, Xie Yuanyun, Kang Chunguo, et al. The capture ofthe Songhua River system in the late Early Pleistocene: Geochemicaland sedimentological records[J]. Acta Geologica Sinica,2020, 94(10): 3144-3160.]

    [19] 魏振宇,謝遠(yuǎn)云,康春國,等. 早更新世松花江水系反轉(zhuǎn):來自荒山巖芯Sr-Nd 同位素特征指示[J]. 沉積學(xué)報(bào),2020,38(6):1192-1203.[Wei Zhenyu, Xie Yuanyun, Kang Chunguo, et al.The inversion of the Songhua River system in the Early Pleistocene:Implications from Sr-Nd isotopic composition in the HarbinHuangshan cores[J]. Acta Sedimentologica Sinica, 2020, 38(6): 1192-1203.]

    [20] 孫忠,賈長青. 松花江流域蓄滯洪區(qū)建設(shè)有關(guān)問題探討[J]. 東北水利水電,2007,25(10):31-33.[Sun Zhong, Jia Changqing.Discussion on the construction of flood storage and detentionarea in Songhua River basin[J]. Water Resources amp; Hydropowerof Northeast China, 2007, 25(10): 31-33.]

    [21] 魏春艷,謝遠(yuǎn)云,康春國,等. 哈爾濱地區(qū)羅家窩棚組地層的沉積學(xué)、礦物學(xué)及地球化學(xué)特征:對沉積環(huán)境的指示[J]. 地質(zhì)科學(xué),2022,57(1):172-189. [Wei Chunyan, Xie Yuanyun,Kang Chunguo, et al. Sedimentological, mineralogical, and geochemicalcharacteristics of the Luojiawopeng Fm. in Harbin: Implicationsfor the sedimentary environment[J]. Chinese Journalof Geology, 2022, 57(1): 172-189.]

    [22] 袁方,謝遠(yuǎn)云,遲云平. 哈爾濱塵暴天氣沉降物的物質(zhì)組成及其對物源的限制[J]. 中國地質(zhì),2018,45(6):1177-1187.[YuanFang, Xie Yuanyun, Chi Yunping. Material characteristics ofdust fallouts during the dust-storm weather in Harbin: Constrainton the provenance[J]. Geology in China, 2018, 45(6): 1177-1187.]

    [23] 孫建華,謝遠(yuǎn)云,康春國,等. 大興安嶺東麓平安鎮(zhèn)白土山組的地層屬性:對物源和沉積環(huán)境的指示[J]. 地層學(xué)雜志,2022,46(2):196-208.[Sun Jianhua, Xie Yuanyun, Kang Chunguo, etal. Stratigraphic properties of the Baitushan Formation inPing’an town, the eastern foot of the Great Hinggan Mountains:An indication of provenance and sedimentary environment[J].Journal of Stratigraphy, 2022, 46(2): 196-208.]

    [24] 孫磊,謝遠(yuǎn)云,康春國,等. 呼倫貝爾沙地重礦物、Sr-Nd 同位素組成及其對亞洲風(fēng)塵系統(tǒng)的指示[J]. 中國地質(zhì),2021,48(6):1965-1974.[Sun Lei, Xie Yuanyun, Kang Chunguo, et al.Heavy minerals, Sr-Nd isotopic composition of sandy land in HulunBuir, Inner Mongolia and their implications for Asian aeoliandust system[J]. Geology in China, 2021, 48(6): 1965-1974.]

    [25] 康春國,李長安,王節(jié)濤,等. 江漢平原沉積物重礦物特征及其對三峽貫通的指示[J]. 地球科學(xué):中國地質(zhì)大學(xué)學(xué)報(bào),2009,34(3):419-427.[Kang Chunguo, Li Chang’an, Wang Jietao, etal. Heavy minerals characteristics of sediments in Jianghanplain and its indication to the forming of the Three Gorges[J].Earth Science: Journal of China University of Geosciences,2009, 34(3): 419-427.]

    [26] 王嘉新,謝遠(yuǎn)云,康春國,等. 哈爾濱荒山巖心重礦物特征對松花江第四紀(jì)水系演化的指示[J]. 第四紀(jì)研究,2020,40(1):79-94.[Wang Jiaxin, Xie Yuanyun, Kang Chunguo, et al. Theindication of the heavy mineral characteristics of the core in HarbinHuangshan to the Quaternary drainage evolution of SonghuaRiver[J]. Quaternary Sciences, 2020, 40(1): 79-94.]

    [27] Taylor S R, McLennan S M. The continental crust: Its compositionand evolution[M]. Oxford: Blackwell Scientific Publications,1985: 57-72.

    [28] 由文智,向芳,黃恒旭,等. 青藏高原東緣宜賓地區(qū)第四紀(jì)河流沉積物中鐵質(zhì)重礦物特征及物源意義[J]. 地學(xué)前緣,2022,29(4):278-292.[You Wenzhi, Xiang Fang, Huang Hengxu, etal. Characteristics and provenance significance of iron-richheavy minerals in Quaternary fluvial sediments in Yibin area,eastern margin of Tibetan Plateau[J]. Earth Science Frontiers,2022, 29(4): 278-292.]

    [29] 田豹. 重礦物物源分析研究進(jìn)展[J]. 中國錳業(yè),2017,35(1):107-109,115. [Tian Bao. A research progress in provenanceanalysis of heavy minerals[J]. China’s Manganese Industry,2017, 35(1): 107-109, 115.]

    [30] 毛光周,劉池洋. 地球化學(xué)在物源及沉積背景分析中的應(yīng)用[J]. 地球科學(xué)與環(huán)境學(xué)報(bào),2011,33(4):337-348.[Mao Guangzhou,Liu Chiyang. Application of geochemistry in provenanceand depositional setting analysis[J]. Journal of Earch Sciencesand Environment, 2011, 33(4): 337-348.]

    [31] Nesbitt H W, Young G M, McLennan S M, et al. Effects ofchemical weathering and sorting on the petrogenesis of siliciclasticsediments, with implications for provenance studies[J]. TheJournal of Geology, 1996, 104(5): 525-542.

    [32] Laceby J P, Evrard O, Smith H G, et al. The challenges and opportunitiesof addressing particle size effects in sediment sourcefingerprinting: A review[J]. Earth-Science Reviews, 2017, 169:85-103.

    [33] 李緒龍,張霞,林春明,等. 常用化學(xué)風(fēng)化指標(biāo)綜述:應(yīng)用與展望[J]. 高校地質(zhì)學(xué)報(bào),2022,28(1):51-63.[Li Xulong, ZhangXia, Lin Chunming, et al. Overview of the application and prospectof common chemical weathering indices[J]. GeologicalJournal of China Universities, 2022, 28(1): 51-63.]

    [34] Nesbitt H W, Young G M. Early Proterozoic climates and platemotions inferred from major element chemistry of lutites[J]. Nature,1982, 299(5885): 715-717.

    [35] Parker A. An index of weathering for silicate rocks[J]. GeologicalMagazine, 1970, 107(6): 501-504.

    [36] Nicholls G D. Environmental studies in sedimentary geochemistry[J]. Science Progress, 1963, 51(201): 12-31.

    [37] Cox R, Lowe D R, Cullers R L. The influence of sediment recyclingand basement composition on evolution of mudrock chemistryin the southwestern United States[J]. Geochimica et CosmochimicaActa, 1995, 59(14): 2919-2940.

    [38] 陳駿,安芷生,劉連文,等. 最近2. 5Ma 以來黃土高原風(fēng)塵化學(xué)組成的變化與亞洲內(nèi)陸的化學(xué)風(fēng)化[J]. 中國科學(xué):地球科學(xué),2001,31(2):136-145. [Chen Jun, An Zhisheng,Liu Lianwen, et al. Variations in chemical compositions of theeolian dust in Chinese Loess Plateau over the past 2. 5 Ma andchemical weathering in the Asian inland[J]. Science China:Earth Sciences, 2001, 31(2): 136-145.]

    [39] Cullers R L, Podkovyrov V N. Geochemistry of the MesoproterozoicLakhanda shales in southeastern Yakutia, Russia: Implicationsfor mineralogical and provenance control, and recycling[J]. Precambrian Research, 2000, 104(1/2): 77-93.

    [40] Ahmad I, Chandra R. Geochemistry of loess-paleosol sedimentsof Kashmir Valley, India: Provenance and weathering[J]. Journalof Asian Earth Sciences, 2013, 66: 73-89.

    [41] Dey S, Rai A K, Chaki A. Palaeoweathering, composition andtectonics of provenance of the Proterozoic intracratonic Kaladgi-Badami Basin, Karnataka, southern India: Evidence from sandstonepetrography and geochemistry[J]. Journal of Asian EarthSciences, 2009, 34(6): 703-715.

    [42] Armstrong-Altrin J S, Machain-Castillo M L. Mineralogy, geochemistry,and radiocarbon ages of deep sea sediments from thegulf of Mexico, Mexico[J]. Journal of South American Earth Sciences,2016, 71: 182-200.

    [43] Ohta T. Geochemistry of Jurassic to earliest Cretaceous depositsin the Nagato Basin, SW Japan: Implication of factor analysis tosorting effects and provenance signatures[J]. Sedimentary Geology,2004, 171(1/2/3/4): 159-180.

    [44] 李秋杭,謝遠(yuǎn)云,康春國,等. 基于人工和TIMA自動(dòng)化方法的松花江水系重礦物組成:對源—匯物源示蹤的指示[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì),2022,42(3):170-183.[Li Qiuhang, XieYuanyun, Kang Chunguo, et al. Heavy mineral composition ofthe Songhua River system identified by manual and TIMAautomatic methods and implications for provenance tracing[J].Marine Geology amp; Quaternary Geology, 2022, 42(3): 170-183.]

    [45] Cullers R L, Basu A, Suttner L J. Geochemical signature ofprovenance in sand-size material in soils and stream sedimentsnear the tobacco root batholith, Montana, U. S. A. [J]. ChemicalGeology, 1988, 70(4): 335-348.

    [46] Cullers R L. The controls on the major- and trace-element evolutionof shales, siltstones and sandstones of Ordovician to Tertiaryage in the Wet Mountains region, Colorado, U. S. A. [J]. ChemicalGeology, 1995, 123(1/2/3/4): 107-131.

    [47] Condie K C, Wronkiewicz D J. The Cr/Th ratio in Precambrianpelites from the Kaapvaal Craton as an index of craton evolution[J]. Earth and Planetary Science Letters, 1990, 97(3/4):256-267.

    [48] Feng R, Kerrich R. Geochemistry of fine-grained clastic sedimentsin the Archean Abitibi greenstone belt, Canada: Implicationsfor provenance and tectonic setting[J]. Geochimica et CosmochimicaActa, 1990, 54(4): 1061-1081.

    [49] Spalletti L A, Queralt I, Matheos S D, et al. Sedimentary petrologyand geochemistry of siliciclastic rocks from the Upper JurassicTordillo Formation (Neuquén Basin, western Argentina): Implicationsfor provenance and tectonic setting[J]. Journal ofSouth American Earth Sciences, 2008, 25(4): 440-463.

    [50] Hayashi K I, Fujisawa H, Holland H D, et al. Geochemistryof~1. 9 Ga sedimentary rocks from northeastern Labrador, Canada[J]. Geochimica et Cosmochimica Acta, 1997, 61(19): 4115-4137.

    [51] Roser B P, Korsch R J. Determination of tectonic setting ofsandstone-mudstone suites using SiO2 content and K2O/Na2Oratio[J]. The Journal of Geology, 1986, 94(5): 635-650.

    [52] Ling Z Y, Yang S L, Xia D S, et al. Source of the aeolian sedimentsin the Yarlung Tsangpo valley and its potential dust contributionto adjacent oceans[J]. Earth Surface Processes and Landforms,2022, 47(7): 1860-1871.

    [53] Xie Y Y, Liu L, Kang C G, et al. Sr-Nd isotopic characteristicsof the northeast sandy land, China and their implications for tracingsources of regional dust[J]. CATENA, 2020, 184: 104303.

    [54] 鄭海飛,郝瑞霞. 普通地球化學(xué)[M]. 北京:北京大學(xué)出版社,2007:76-79.[Zheng Haifei, Hao Ruixia. General geochemistry[M]. Beijing: Peking University Press, 2007: 76-79.]

    [55] Yang J D, Li G J, Rao W B, et al. Isotopic evidences for provenanceof East Asian dust[J]. Atmospheric Environment, 2009, 43(29): 4481-4490.

    [56] Rao W B, Mao C P, Wang Y G, et al. Using Nd-Sr isotopes andrare earth elements to study sediment provenance of the modernradial sand ridges in the southwestern Yellow Sea[J]. AppliedGeochemistry, 2017, 81: 23-35.

    [57] Bayon G, Toucanne S, Skonieczny C, et al. Rare earth elementsand neodymium isotopes in world river sediments revisited[J].Geochimica et Cosmochimica Acta, 2015, 170: 17-38.

    [58] Morton A C, Hallsworth C R. Processes controlling the compositionof heavy mineral assemblages in sandstones[J]. SedimentaryGeology, 1999, 124(1/2/3/4): 3-29.

    [59] 何夢穎. 長江河流沉積物礦物學(xué)、地球化學(xué)和碎屑鋯石年代學(xué)物源示蹤研究[D]. 南京:南京大學(xué),2014.[He Mengying.The provenance study on the Yangtze River sediments, based onmineralogy, geochemistry and detrital zircon dating[D]. Nanjing:Nanjing University, 2014.]

    基金項(xiàng)目:國家自然科學(xué)基金項(xiàng)目(42171006,41871013)

    猜你喜歡
    特征
    抓住特征巧觀察
    離散型隨機(jī)變量的分布列與數(shù)字特征
    具有兩個(gè)P’維非線性不可約特征標(biāo)的非可解群
    月震特征及與地震的對比
    如何表達(dá)“特征”
    被k(2≤k≤16)整除的正整數(shù)的特征
    不忠誠的四個(gè)特征
    詈語的文化蘊(yùn)含與現(xiàn)代特征
    新聞傳播(2018年11期)2018-08-29 08:15:24
    抓住特征巧觀察
    基于特征篩選的模型選擇
    久久国产乱子伦精品免费另类| 桃色一区二区三区在线观看| 啦啦啦观看免费观看视频高清| 亚洲一区中文字幕在线| 亚洲国产看品久久| 国产成人精品久久二区二区91| 99热6这里只有精品| 99热这里只有是精品50| 日本黄色视频三级网站网址| 免费在线观看完整版高清| 久久亚洲真实| 成人18禁在线播放| 亚洲av五月六月丁香网| 一a级毛片在线观看| 日韩欧美免费精品| 国产单亲对白刺激| 国产亚洲精品一区二区www| 大型黄色视频在线免费观看| 亚洲av电影在线进入| 国产熟女午夜一区二区三区| 日韩大尺度精品在线看网址| 少妇粗大呻吟视频| 亚洲熟妇中文字幕五十中出| 亚洲人成电影免费在线| 亚洲精品av麻豆狂野| 人妻夜夜爽99麻豆av| 1024香蕉在线观看| 国产精品久久视频播放| 色尼玛亚洲综合影院| 757午夜福利合集在线观看| 国产亚洲精品第一综合不卡| 国产黄色小视频在线观看| 久久这里只有精品19| 夜夜躁狠狠躁天天躁| 亚洲中文字幕日韩| 99精品欧美一区二区三区四区| 欧美成人一区二区免费高清观看 | 国产又黄又爽又无遮挡在线| 国产午夜福利久久久久久| 精品国产美女av久久久久小说| 免费人成视频x8x8入口观看| 成人18禁在线播放| 国产99白浆流出| 18禁国产床啪视频网站| 麻豆国产av国片精品| 黑人欧美特级aaaaaa片| 免费在线观看影片大全网站| www.999成人在线观看| 成年人黄色毛片网站| 色尼玛亚洲综合影院| 亚洲欧美精品综合一区二区三区| 亚洲美女视频黄频| 久久天躁狠狠躁夜夜2o2o| 美女扒开内裤让男人捅视频| 在线a可以看的网站| 99re在线观看精品视频| 成人午夜高清在线视频| 正在播放国产对白刺激| 国产午夜福利久久久久久| 欧美黑人欧美精品刺激| 免费在线观看黄色视频的| 啦啦啦韩国在线观看视频| 我的老师免费观看完整版| 五月玫瑰六月丁香| 一本一本综合久久| 欧美最黄视频在线播放免费| 久久欧美精品欧美久久欧美| 午夜福利免费观看在线| 99久久综合精品五月天人人| 成人国产综合亚洲| 国产av一区二区精品久久| 母亲3免费完整高清在线观看| 一卡2卡三卡四卡精品乱码亚洲| 久久久久国产精品人妻aⅴ院| 久久人人精品亚洲av| 俄罗斯特黄特色一大片| 亚洲av电影不卡..在线观看| 中文亚洲av片在线观看爽| 国产免费av片在线观看野外av| 国产亚洲精品久久久久5区| x7x7x7水蜜桃| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品一区二区三区四区久久| 亚洲人成电影免费在线| 超碰成人久久| 国产欧美日韩一区二区三| 桃红色精品国产亚洲av| 欧美一级a爱片免费观看看 | 熟女少妇亚洲综合色aaa.| 亚洲一区中文字幕在线| 深夜精品福利| 最近在线观看免费完整版| 午夜久久久久精精品| 黑人巨大精品欧美一区二区mp4| 香蕉av资源在线| 免费人成视频x8x8入口观看| 久久香蕉激情| 午夜两性在线视频| 麻豆久久精品国产亚洲av| 三级男女做爰猛烈吃奶摸视频| 一边摸一边做爽爽视频免费| 国产99白浆流出| 手机成人av网站| 我的老师免费观看完整版| 无遮挡黄片免费观看| 丰满人妻一区二区三区视频av | 免费在线观看成人毛片| 亚洲第一电影网av| 国产午夜精品论理片| 在线观看www视频免费| 国产精品久久久av美女十八| 国产三级黄色录像| 国产成人影院久久av| 精品一区二区三区av网在线观看| 婷婷精品国产亚洲av在线| 一夜夜www| 久久伊人香网站| 日本成人三级电影网站| 国产av不卡久久| x7x7x7水蜜桃| 一边摸一边做爽爽视频免费| 极品教师在线免费播放| 色综合欧美亚洲国产小说| 国产精品免费视频内射| 国产一区二区激情短视频| 免费在线观看亚洲国产| 久久国产精品影院| 久久久久久大精品| 婷婷丁香在线五月| 色哟哟哟哟哟哟| 久久久久久大精品| 欧美国产日韩亚洲一区| 日本免费一区二区三区高清不卡| 人成视频在线观看免费观看| 一级毛片女人18水好多| 人人妻,人人澡人人爽秒播| 成在线人永久免费视频| 亚洲男人的天堂狠狠| 国产视频一区二区在线看| 十八禁网站免费在线| 国产精品美女特级片免费视频播放器 | 久久精品亚洲精品国产色婷小说| 黄片小视频在线播放| 一级作爱视频免费观看| 亚洲午夜精品一区,二区,三区| 国产精品久久久久久精品电影| 国产欧美日韩一区二区精品| 成年版毛片免费区| 国产精品日韩av在线免费观看| 日本撒尿小便嘘嘘汇集6| www日本在线高清视频| 在线永久观看黄色视频| 久9热在线精品视频| 国产又黄又爽又无遮挡在线| 国产1区2区3区精品| av有码第一页| 日韩有码中文字幕| 人妻久久中文字幕网| 久久久精品大字幕| 精品不卡国产一区二区三区| 黄色女人牲交| 色精品久久人妻99蜜桃| 午夜老司机福利片| 久久久久免费精品人妻一区二区| 国产三级中文精品| 亚洲熟妇熟女久久| 18禁裸乳无遮挡免费网站照片| 国内久久婷婷六月综合欲色啪| 日本熟妇午夜| 人人妻人人澡欧美一区二区| 国产av不卡久久| 最近最新免费中文字幕在线| 非洲黑人性xxxx精品又粗又长| 国产aⅴ精品一区二区三区波| 不卡av一区二区三区| 天天躁夜夜躁狠狠躁躁| 久久天堂一区二区三区四区| 一个人免费在线观看电影 | www.www免费av| 国产v大片淫在线免费观看| 真人做人爱边吃奶动态| 久久久久久免费高清国产稀缺| 日本 av在线| 狠狠狠狠99中文字幕| 精品久久蜜臀av无| 色哟哟哟哟哟哟| 午夜亚洲福利在线播放| 亚洲精品av麻豆狂野| 美女午夜性视频免费| 国产亚洲av嫩草精品影院| 久久久久免费精品人妻一区二区| www.999成人在线观看| 欧美3d第一页| 成人特级黄色片久久久久久久| 亚洲无线在线观看| 两个人看的免费小视频| av在线播放免费不卡| 国产欧美日韩一区二区精品| 老熟妇乱子伦视频在线观看| 国产欧美日韩精品亚洲av| 亚洲第一欧美日韩一区二区三区| 天天躁夜夜躁狠狠躁躁| 在线观看免费日韩欧美大片| 亚洲一区中文字幕在线| av在线天堂中文字幕| 精品久久久久久久久久久久久| 亚洲全国av大片| 大型av网站在线播放| 国产探花在线观看一区二区| 黄色女人牲交| 亚洲第一欧美日韩一区二区三区| 亚洲av成人不卡在线观看播放网| 亚洲五月天丁香| 日日摸夜夜添夜夜添小说| 国产不卡一卡二| 欧美人与性动交α欧美精品济南到| 欧美乱色亚洲激情| 女生性感内裤真人,穿戴方法视频| 欧美日韩亚洲国产一区二区在线观看| 身体一侧抽搐| 久久久久久久久久黄片| 啪啪无遮挡十八禁网站| 性欧美人与动物交配| 91九色精品人成在线观看| 又粗又爽又猛毛片免费看| 黄片小视频在线播放| 一边摸一边抽搐一进一小说| 成人亚洲精品av一区二区| www日本在线高清视频| 久久人人精品亚洲av| 搡老妇女老女人老熟妇| 美女 人体艺术 gogo| 9191精品国产免费久久| 婷婷精品国产亚洲av在线| 日本免费a在线| 日韩欧美国产在线观看| 精品第一国产精品| 别揉我奶头~嗯~啊~动态视频| 久久久久精品国产欧美久久久| 国产视频一区二区在线看| 波多野结衣高清无吗| 一区福利在线观看| 国产黄色小视频在线观看| 欧美色欧美亚洲另类二区| 一级毛片高清免费大全| 欧美+亚洲+日韩+国产| av有码第一页| 国产日本99.免费观看| 无限看片的www在线观看| 五月伊人婷婷丁香| 日韩av在线大香蕉| 日韩欧美在线二视频| 国产视频一区二区在线看| 99精品欧美一区二区三区四区| 麻豆国产av国片精品| 精品久久久久久久人妻蜜臀av| av天堂在线播放| 免费观看人在逋| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩黄片免| 淫秽高清视频在线观看| 国产精品电影一区二区三区| 国产成年人精品一区二区| 色综合欧美亚洲国产小说| 伊人久久大香线蕉亚洲五| 男人的好看免费观看在线视频 | 亚洲人成77777在线视频| 制服人妻中文乱码| 少妇被粗大的猛进出69影院| 欧美极品一区二区三区四区| 日韩欧美精品v在线| 好男人在线观看高清免费视频| 手机成人av网站| 听说在线观看完整版免费高清| 久久热在线av| 国产私拍福利视频在线观看| 欧美国产日韩亚洲一区| 国产主播在线观看一区二区| 91麻豆av在线| 美女免费视频网站| 免费看日本二区| 午夜福利高清视频| 国产欧美日韩一区二区精品| 999久久久精品免费观看国产| 日韩av在线大香蕉| 又黄又爽又免费观看的视频| 久久精品国产亚洲av高清一级| 在线国产一区二区在线| 亚洲精品粉嫩美女一区| 欧美最黄视频在线播放免费| 天堂动漫精品| 黄色a级毛片大全视频| 在线观看一区二区三区| 色综合欧美亚洲国产小说| av天堂在线播放| 一级作爱视频免费观看| 麻豆久久精品国产亚洲av| 在线观看午夜福利视频| 一级黄色大片毛片| 国产精品久久视频播放| 一二三四在线观看免费中文在| 免费观看人在逋| 国产成人精品久久二区二区免费| 19禁男女啪啪无遮挡网站| 两个人看的免费小视频| 免费在线观看完整版高清| 欧洲精品卡2卡3卡4卡5卡区| 国产真人三级小视频在线观看| aaaaa片日本免费| 久久久久久人人人人人| 久久香蕉激情| 一进一出抽搐gif免费好疼| a在线观看视频网站| 国产乱人伦免费视频| 欧美日韩乱码在线| 日韩国内少妇激情av| 国产一区二区三区在线臀色熟女| 国产成人精品久久二区二区免费| 91成年电影在线观看| 国产探花在线观看一区二区| 变态另类丝袜制服| 欧美最黄视频在线播放免费| 久久香蕉国产精品| 成人一区二区视频在线观看| 19禁男女啪啪无遮挡网站| 日本一本二区三区精品| 亚洲男人天堂网一区| 精品久久久久久久久久久久久| 国产精品香港三级国产av潘金莲| 亚洲男人的天堂狠狠| 精华霜和精华液先用哪个| 久久久久久亚洲精品国产蜜桃av| 精品久久久久久,| 日本撒尿小便嘘嘘汇集6| 亚洲 国产 在线| 欧美色视频一区免费| 亚洲电影在线观看av| 久久久久国内视频| 欧美成人午夜精品| 给我免费播放毛片高清在线观看| 亚洲电影在线观看av| 18禁裸乳无遮挡免费网站照片| 黄色成人免费大全| 国产成人精品无人区| 欧美乱色亚洲激情| 香蕉丝袜av| 欧美高清成人免费视频www| 九色成人免费人妻av| 国产精品一区二区免费欧美| 中文字幕人妻丝袜一区二区| 中文字幕久久专区| 亚洲精品色激情综合| 国产又黄又爽又无遮挡在线| 五月玫瑰六月丁香| 国产av麻豆久久久久久久| 亚洲国产欧洲综合997久久,| 一级毛片女人18水好多| 国产亚洲欧美在线一区二区| 亚洲成av人片在线播放无| 90打野战视频偷拍视频| 国产亚洲精品一区二区www| 99国产综合亚洲精品| 中文资源天堂在线| 18禁美女被吸乳视频| 日日爽夜夜爽网站| 国产高清激情床上av| 欧美黄色片欧美黄色片| 久久欧美精品欧美久久欧美| 国内揄拍国产精品人妻在线| 欧美zozozo另类| 成人18禁高潮啪啪吃奶动态图| 成人av一区二区三区在线看| 亚洲七黄色美女视频| 国产av一区在线观看免费| 男人舔女人的私密视频| 久久精品91无色码中文字幕| 99国产精品99久久久久| 日韩欧美一区二区三区在线观看| 亚洲成a人片在线一区二区| 久久精品91蜜桃| 制服丝袜大香蕉在线| 国内精品久久久久精免费| 啦啦啦免费观看视频1| 亚洲一卡2卡3卡4卡5卡精品中文| 手机成人av网站| 90打野战视频偷拍视频| 精品久久久久久久久久久久久| 女人高潮潮喷娇喘18禁视频| 青草久久国产| 色老头精品视频在线观看| 国产精品自产拍在线观看55亚洲| 99热这里只有是精品50| 男女下面进入的视频免费午夜| 国产av又大| 一区二区三区激情视频| 欧美日韩一级在线毛片| 男人舔奶头视频| 成年免费大片在线观看| 在线视频色国产色| 日本三级黄在线观看| 国产免费av片在线观看野外av| 中文字幕最新亚洲高清| www.999成人在线观看| 成人欧美大片| 国产一区二区三区在线臀色熟女| 两性午夜刺激爽爽歪歪视频在线观看 | 成年免费大片在线观看| av在线播放免费不卡| 亚洲avbb在线观看| 久久天堂一区二区三区四区| 怎么达到女性高潮| 亚洲人成77777在线视频| 国产精品亚洲美女久久久| 欧美成人性av电影在线观看| 看黄色毛片网站| 男女床上黄色一级片免费看| 99re在线观看精品视频| 一边摸一边抽搐一进一小说| 天堂av国产一区二区熟女人妻 | 亚洲熟妇中文字幕五十中出| 成人av在线播放网站| 亚洲欧美日韩东京热| 变态另类丝袜制服| 欧美日韩亚洲国产一区二区在线观看| 又粗又爽又猛毛片免费看| 99久久无色码亚洲精品果冻| 国产亚洲精品一区二区www| 天堂√8在线中文| 亚洲欧洲精品一区二区精品久久久| 国产野战对白在线观看| 在线观看免费日韩欧美大片| 亚洲精品国产精品久久久不卡| 国产成人啪精品午夜网站| 欧美日本亚洲视频在线播放| 亚洲人成电影免费在线| 一级黄色大片毛片| 午夜免费激情av| 一进一出抽搐gif免费好疼| 在线看三级毛片| 五月伊人婷婷丁香| 久久国产精品人妻蜜桃| videosex国产| 在线十欧美十亚洲十日本专区| 亚洲五月天丁香| 精品国产亚洲在线| 久久久精品国产亚洲av高清涩受| 巨乳人妻的诱惑在线观看| www.熟女人妻精品国产| АⅤ资源中文在线天堂| 亚洲专区字幕在线| 精品久久久久久久末码| 日本黄大片高清| x7x7x7水蜜桃| 日本 av在线| 亚洲av成人精品一区久久| 欧美一区二区国产精品久久精品 | 99在线人妻在线中文字幕| 亚洲成av人片在线播放无| 两个人视频免费观看高清| 亚洲专区国产一区二区| 在线观看免费午夜福利视频| 999久久久国产精品视频| 午夜精品久久久久久毛片777| 两个人看的免费小视频| 国产乱人伦免费视频| 丝袜美腿诱惑在线| 国产精品野战在线观看| 欧美黑人精品巨大| 欧美中文日本在线观看视频| 99国产极品粉嫩在线观看| 日本一二三区视频观看| 国内精品一区二区在线观看| 亚洲人成77777在线视频| 国产精品电影一区二区三区| 欧美黄色片欧美黄色片| 久久伊人香网站| 一区福利在线观看| 婷婷丁香在线五月| 亚洲国产精品999在线| 国产97色在线日韩免费| av天堂在线播放| 亚洲人成伊人成综合网2020| 亚洲专区中文字幕在线| 国产真实乱freesex| 老司机靠b影院| 午夜免费激情av| 在线观看免费视频日本深夜| 久久天躁狠狠躁夜夜2o2o| 免费无遮挡裸体视频| 欧美zozozo另类| 欧美激情久久久久久爽电影| 搡老妇女老女人老熟妇| 免费在线观看视频国产中文字幕亚洲| av在线天堂中文字幕| 亚洲中文日韩欧美视频| 免费人成视频x8x8入口观看| 成人av一区二区三区在线看| 精品午夜福利视频在线观看一区| 天堂动漫精品| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产中文字幕在线视频| 国产精品一区二区三区四区久久| 亚洲专区字幕在线| 禁无遮挡网站| 久久国产乱子伦精品免费另类| 天堂影院成人在线观看| 精品久久久久久久毛片微露脸| 亚洲色图av天堂| 后天国语完整版免费观看| 无遮挡黄片免费观看| 久久午夜综合久久蜜桃| 欧美性猛交╳xxx乱大交人| 国语自产精品视频在线第100页| 亚洲一卡2卡3卡4卡5卡精品中文| 在线视频色国产色| 精品一区二区三区av网在线观看| 欧美日韩乱码在线| 999久久久精品免费观看国产| 欧美3d第一页| 欧美黑人欧美精品刺激| 欧美绝顶高潮抽搐喷水| 成人手机av| 无遮挡黄片免费观看| 中出人妻视频一区二区| 精华霜和精华液先用哪个| 美女午夜性视频免费| 国产精品亚洲av一区麻豆| 天天躁狠狠躁夜夜躁狠狠躁| 18禁黄网站禁片午夜丰满| 国产午夜福利久久久久久| 亚洲五月婷婷丁香| 老汉色∧v一级毛片| tocl精华| 一区二区三区激情视频| 一边摸一边抽搐一进一小说| 成年版毛片免费区| 亚洲国产欧美一区二区综合| 亚洲欧美精品综合一区二区三区| 中亚洲国语对白在线视频| 在线a可以看的网站| a在线观看视频网站| 日本黄色视频三级网站网址| 香蕉国产在线看| 国产精品综合久久久久久久免费| 国产精品美女特级片免费视频播放器 | 亚洲男人的天堂狠狠| 他把我摸到了高潮在线观看| 久久久久久久久免费视频了| 88av欧美| 国产精品野战在线观看| 精品久久久久久久毛片微露脸| 欧美乱色亚洲激情| 99国产精品一区二区蜜桃av| 国产黄片美女视频| 一本久久中文字幕| 亚洲成人久久性| 99久久国产精品久久久| 人人妻人人澡欧美一区二区| 精品午夜福利视频在线观看一区| 国产私拍福利视频在线观看| 国产伦一二天堂av在线观看| 国产av又大| 午夜福利欧美成人| 国产1区2区3区精品| 最近最新免费中文字幕在线| 精品人妻1区二区| 九色国产91popny在线| 性色av乱码一区二区三区2| 免费在线观看影片大全网站| 法律面前人人平等表现在哪些方面| 天堂√8在线中文| 欧美不卡视频在线免费观看 | 在线观看一区二区三区| 一级毛片高清免费大全| 99riav亚洲国产免费| 亚洲欧美日韩无卡精品| 日日夜夜操网爽| 久久亚洲真实| av中文乱码字幕在线| 精品一区二区三区av网在线观看| 大型av网站在线播放| 嫩草影视91久久| 最新美女视频免费是黄的| 日韩欧美免费精品| 少妇粗大呻吟视频| 成人18禁在线播放| 香蕉丝袜av| 国产高清视频在线播放一区| 国产精品 欧美亚洲| 俄罗斯特黄特色一大片| 我的老师免费观看完整版| 欧洲精品卡2卡3卡4卡5卡区| 听说在线观看完整版免费高清| 国产一区二区三区视频了| 欧美中文综合在线视频| av在线天堂中文字幕| 我的老师免费观看完整版| 久久伊人香网站| 日韩欧美免费精品| 老熟妇仑乱视频hdxx| 精品日产1卡2卡| 丰满人妻一区二区三区视频av | 免费在线观看黄色视频的| 不卡一级毛片| 日本一区二区免费在线视频| 免费看a级黄色片| 欧美丝袜亚洲另类 | 人妻丰满熟妇av一区二区三区| 黄色视频不卡| 免费在线观看视频国产中文字幕亚洲| 欧美zozozo另类| 麻豆成人av在线观看| 给我免费播放毛片高清在线观看| 久久久久久九九精品二区国产 | 亚洲成人久久性|