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

    三峽水庫典型支流消落帶泥沙顆粒態(tài)磷復(fù)合指紋示蹤研究

    2019-12-19 01:14:12陳太麗史忠林王永艷嚴(yán)冬春文安邦陳佳村
    農(nóng)業(yè)工程學(xué)報 2019年20期
    關(guān)鍵詞:懸移質(zhì)溪河物源

    陳太麗,史忠林,王永艷,嚴(yán)冬春,文安邦,陳佳村

    三峽水庫典型支流消落帶泥沙顆粒態(tài)磷復(fù)合指紋示蹤研究

    陳太麗1,2,史忠林1※,王永艷1,2,嚴(yán)冬春1,文安邦1,陳佳村1,2

    (1. 中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所,成都 610041;2. 中國科學(xué)院大學(xué),北京 100049)

    三峽水庫消落帶是庫區(qū)水域與周邊陸地環(huán)境的關(guān)鍵過渡地帶,周期性反季節(jié)干濕交替使其具有強烈的物質(zhì)交換特征。辨析消落帶泥沙及其吸附的顆粒態(tài)磷的來源對消落帶土壤污染防治和環(huán)境效應(yīng)評估以及三峽水環(huán)境保護具有重要意義。以三峽庫區(qū)汝溪河支流不同高程(145~155、>155~165、>165~175 m)消落帶為研究對象,運用復(fù)合指紋技術(shù)查明,消落帶泥沙中顆粒態(tài)磷的主要來源是長江干流懸移質(zhì)和汝溪河上游懸移質(zhì)。淹水期間長江干流江水頂托引起的泥沙沉積是顆粒態(tài)磷的主要來源,在>165~175 m高程帶對顆粒態(tài)磷的貢獻達(dá)到最大(54.5%)。雨季初期支流上游懸移質(zhì)對145~155 m高程消落帶的顆粒態(tài)磷貢獻最大(51.6%),而隨高程的增加貢獻率減少。消落帶上方的土壤侵蝕產(chǎn)沙主要堆積在>155~165和>165~175 m高程范圍內(nèi),導(dǎo)致消落帶上方土壤對泥沙和顆粒態(tài)磷的貢獻率都隨高程的增加而增加。

    磷;泥沙;元素;消落帶;復(fù)合指紋技術(shù)

    0 引 言

    三峽水庫的蓄水運行,在庫區(qū)干支流兩岸形成了落差30 m,總面積約349 km2的消落帶[1]。蓄水期間,由于水位加深,流速減緩,河流泥沙及其攜帶的污染物在消落帶沉積;而出露期正值雨季,頻繁的降雨和暴曬引發(fā)嚴(yán)重的水土流失,沉積的泥沙及污染物可能在徑流沖刷作用下再次進入水體。《長江三峽工程生態(tài)與環(huán)境監(jiān)測公報》[2-3]顯示,三峽蓄水以來,盡管庫區(qū)長江干流水質(zhì)總體為良,但受蓄水后干流回水頂托作用影響,很多支流水體富營養(yǎng)化加劇,部分支流回水區(qū)在春季蓄水期水華現(xiàn)象突出。

    磷是水體富營養(yǎng)化的主要限制因子,查明水體磷的來源與遷移是緩解水質(zhì)惡化的重要前提。水環(huán)境中磷的形態(tài)可按物理態(tài)和化學(xué)態(tài)進行劃分,按物理態(tài)可分為溶解態(tài)和顆粒態(tài);按化學(xué)態(tài)可分為無機態(tài)(正磷酸鹽、縮合磷酸鹽等)和有機態(tài)(三磷酸腺苷、磷酸酯等)。而流域內(nèi)磷的主要輸出形態(tài)為顆粒態(tài)[4-6],可占磷輸出總量的90%以上[7]。泥沙是顆粒態(tài)磷的主要載體,泥沙來源對磷的來源有重要影響[8-9]。消落帶作為泥沙的堆積場所,同時也是顆粒態(tài)磷的重要匯,但目前關(guān)于消落帶磷素來源的研究尚未見報道。本研究以庫區(qū)典型支流汝溪河消落帶為研究對象,運用復(fù)合指紋示蹤技術(shù),判別不同高程消落帶沉積泥沙中顆粒態(tài)磷的來源及其相對貢獻,以期為三峽工程蓄水后消落帶污染特征及環(huán)境效應(yīng)評估以及三峽水環(huán)境保護提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    汝溪河為長江一級支流,地處三峽庫區(qū)腹心地帶,行政上主要隸屬于重慶忠縣(圖1)。流域面積720 km2,主河道長54.5 km,多年平均徑流量14.9億m3。汝溪河流域?qū)賮啛釒駶櫦撅L(fēng)氣候區(qū),年均溫18.5 ℃,年均降水1 140 mm,主要集中于5—9月。土地利用以林地、耕地和果園為主,土壤主要為紫色土。三峽庫區(qū)山高坡陡,而消落帶大多在入江口或地勢低洼(海拔145~175 m)的狹小范圍內(nèi)形成,因此消落帶長度與流域長度相比相差較大。汝溪河流域消落帶總長約6 km(距長江匯合口距離),面積6.4 km2,以緩坡型為主,約占80%,蓄水期泥沙沉積明顯。

    1.2 泥沙物源劃分

    應(yīng)用指紋識別技術(shù)首先需在綜合考慮研究目的、時空尺度、物源間差異性等多種因素的基礎(chǔ)上對潛在泥沙物源進行劃分。流域泥沙物源主要分為類型物源(如侵蝕類型或土地利用類型[10-13])和空間物源(如次級流域或地質(zhì)亞區(qū)[14-17])2類,其中以類型物源研究居多。實地調(diào)查表明,對于三峽支流消落帶,受蓄水后長江干流回水頂托影響,河流流速減緩,干流輸送的泥沙會在消落帶發(fā)生沉積;同時,洪水季節(jié)和消落帶淹水期支流所在流域上游來沙也是潛在泥沙來源之一;此外,消落帶上方(高程175 m以上)的農(nóng)用地在雨季的坡面侵蝕也可能對消落帶泥沙有所貢獻。因此,本研究將支流消落帶沉積泥沙源地劃分為消落帶上方土壤(S1)、長江干流懸移質(zhì)(S2)和汝溪河上游懸移質(zhì)(S3)3類(圖1)。

    圖1 研究區(qū)地理位置與采樣點示意圖

    1.3 采樣與分析

    1.3.1 沉積泥沙樣品采集

    汝溪河消落帶在蓄水前主要為農(nóng)耕地或水田,地勢較緩,泥沙淤積明顯。通過實地考察,沿程選擇6個沉積斷面(圖1),每個斷面劃分145~155、>155~165、>165~175 m 3個高程。于2015年8月消落帶出露期,在不同高程帶選擇泥沙易于沉積的平整地塊鋪設(shè)1 m× 1 m的人工草皮收集泥沙。泥沙收集至2016年7月結(jié)束,涵蓋1個完整的蓄水—落干周期。

    1.3.2 泥沙物源樣采集與歸類

    在沉積物采樣斷面175 m高程以上靠近消落帶且易發(fā)生侵蝕的地塊采集0~2 cm表層土壤,每一斷面多點采集后將樣品混合(樣本數(shù)為6),作為消落帶上方產(chǎn)沙物源樣(S1);在距離汝溪河與長江干流匯合口約1 km處的干流河道上布設(shè)浮動式泥沙收集桶,采集干流懸移質(zhì)(S2),采樣時間為2015年9月—2016年6月,頻率為每月1次;同時,在距匯合口約5 km處的汝溪河河道上放置泥沙收集桶,收集汝溪河上游懸移質(zhì)(S3),采樣時間和頻率與干流懸移質(zhì)一致。

    由于不同高程消落帶淹水時間不同(圖2),沉積泥沙的物源組成具有差異,消落帶出露期采集的干、支流懸移質(zhì)不作為其潛在物源。因此對于145~155 m高程消落帶,物源S2和S3的采樣時間為2015年9月—2016年6月(=9),而對于>155~165和>165~175 m高程消落帶,S2和S3的采樣時間則分別為2015年9月—2016年3月(=6)和2015年9月—2016年1月(=4)。

    1.3.3 樣品測試

    消落帶沉積泥沙包裹完整后帶回實驗室沖洗、沉淀、過濾,于60 ℃下烘干至恒質(zhì)量。物源樣經(jīng)自然風(fēng)干,手工去除植物根系等雜質(zhì),研磨后過2 mm篩。取5 g左右沉積泥沙樣,用10% H2O2和10% HCl分別去除有機質(zhì)和碳酸鹽,采用MasterSizer2000激光粒度儀(Malvern Instruments, Malvern, England)測定顆粒組成(超聲波分散2 min)。結(jié)果顯示,不同高程消落帶沉積泥沙以粉粒和黏粒為主,>90%的泥沙粒徑小于63m。因此將所有物源樣和沉積泥沙樣過63m篩,隨后的分析測試均限于<63m組分。測試指標(biāo)包括TOC、TN、TP、K、Mg、Na、Ca、Fe、Al、Cd、Co、Cr、Cu、Mn、Ni、Pb、Ti、Zn共18個,測試方法詳見文獻[18]。

    圖2 不同高程消落帶淹水時間

    1.4 指紋因子篩選及模型求解

    1.4.1 復(fù)合指紋組合構(gòu)建

    對各指標(biāo)測值進行Shapiro-Wilk正態(tài)檢驗[18],結(jié)果顯示不同高程的物源樣和目標(biāo)樣的指紋特征大部分都通過檢驗,數(shù)據(jù)滿足正態(tài)分布(表1)。由于樣品數(shù)量相對較少,而均值相較中值對極端值更敏感,因此使用物源樣和目標(biāo)樣的中值作進一步分析。

    指紋技術(shù)假設(shè)用于模型計算的指紋因子在泥沙遷移過程中具有穩(wěn)定性,從而使泥沙和源地的指紋特征直接可比[19],但這一穩(wěn)定性目前很難直接界定。泥沙輸移過程可能伴隨著粒徑的分選與有機質(zhì)的富集/消耗,這些過程可能會對指紋因子的濃度產(chǎn)生影響。一些研究嘗試通過采用顆粒校正或有機質(zhì)校正來降低這一影響[20-21],但該方法未考慮不同因子對顆粒/有機質(zhì)的敏感程度,采用單一的校正系數(shù)可能帶來過度校正[22]。因此本研究未進行顆粒/有機質(zhì)校正[23],而采用范圍檢驗,即若目標(biāo)泥沙指紋中值落在物源對應(yīng)指紋中值的范圍之內(nèi),則該指紋因子可視為穩(wěn)定[23]。

    為篩選在泥沙源地間具有顯著差異的指紋因子,對通過范圍檢驗的指紋因子進行Kruskal-Wallis H檢驗(KW-H);然后對通過KW-H檢驗的指紋因子(=0.05)進行多元判別分析(discriminant function analysis,DFA),從中篩選出可以區(qū)分不同物源的最佳指紋組合。

    1.4.2 混合模型求解

    采用多元混合模型,求3種潛在物源對不同高程消落帶沉積泥沙的相對貢獻[24]。

    式中R為殘差平方和;C為沉積泥沙指紋因子的濃度,g/kg或mg/kg;C為泥沙源地中指紋因子的濃度,g/kg或mg/kg;P為源地的泥沙貢獻百分比,%;為指紋因子數(shù)量;為泥沙源地數(shù)量。模型運行的邊界條件是所有參與計算的各源地泥沙貢獻百分比P非負(fù)且總和為100%。在滿足上述條件下,當(dāng)R最小時,即可得到各物源的泥沙貢獻百分比。

    利用Crystal Ball軟件的優(yōu)化求解功能,對物源和泥沙各指紋因子進行分布擬合;采用拉丁超立方抽樣方法,進行1 000次重復(fù)求解,尋求最優(yōu)解。誤差分析采用擬合優(yōu)度(goodness of fit,GOF)檢測方法[25],檢驗混合模型對樣本觀測值的擬合程度(式(2))。當(dāng)GOF>0.8時,認(rèn)為模型計算結(jié)果可以接受。

    表1 各指標(biāo)Shapiro-Wilk正態(tài)檢驗結(jié)果

    注:*,真實顯著水平的下限。S1:消落帶上方土壤;S2:長江干流懸移質(zhì);S3:汝溪河上游懸移質(zhì)。下同。

    Note: *, lower bound of the true significance. S1: topsoil above riparian zone; S2: suspended sediment from Yangtze River; S3: suspended sediment from Ruxi River. Same as below.

    1.5 不同源地顆粒態(tài)磷貢獻計算

    忽略泥沙輸移過程中磷的吸附與解析,不同源地對消落帶泥沙中磷的貢獻可由下式求算[9]:

    式中L為泥沙源地的顆粒態(tài)磷貢獻率,%;P為源地的泥沙貢獻百分比,%;C為源地磷素平均濃度,g/kg。

    2 結(jié)果與分析

    2.1 最佳指紋因子組合確定

    不同高程消落帶沉積泥沙及其對應(yīng)物源的測試指標(biāo)值列于表2,指紋判別因子逐步篩選結(jié)果見表3。范圍檢驗表明,18個測試指標(biāo)中TN、TOC、TP、Na、Ti、Mn、Cd 7個因子在3個高程均通過檢驗,可認(rèn)為這些因子在隨泥沙輸移的過程中保持了較好的穩(wěn)定性;Co和Ni在不同高程沉積泥沙中的濃度均超過其對應(yīng)物源濃度的范圍,可能是在輸移過程中因泥沙顆粒的變化而發(fā)生了一定程度的富集,因此將其視為不穩(wěn)定元素予以剔除。

    進一步采用KW-H檢驗和DFA構(gòu)建的145~155 m高程泥沙來源最佳指紋組合為Cd、Ti、Mn,>155~165 m為Mg、Mn、TP、TOC,>165~175 m為Ca、TN和TP,不同指紋組合對源地樣品的正確判別率最低達(dá)到95.8%,判別效果較好(表3)。

    表2 不同高程消落帶泥沙及其物源測試指標(biāo)中值

    表3 不同高程消落帶泥沙來源復(fù)合指紋因子

    2.2 消落帶沉積泥沙來源解析

    不同高程消落帶沉積泥沙物源相對貢獻率如圖3a所示,擬合優(yōu)度檢驗GOF值介于0.88~0.95。模型求解結(jié)果表明,145~155 m高程消落帶泥沙主要來源于長江干流懸移質(zhì)(S2)和汝溪河上游懸移質(zhì)(S3),其貢獻率分別為45.8%和38.4%,消落帶上部土壤(S1)貢獻率僅為15.8%;對于>155~165和>165~175 m高程消落帶,3種物源對沉積泥沙的貢獻比例分布基本一致。對比不同高程消落帶泥沙來源發(fā)現(xiàn),長江干流懸移質(zhì)的貢獻均為最大,平均貢獻率為44.1%,物源貢獻的差異主要表現(xiàn)在消落帶上方土壤和汝溪河上游懸移質(zhì)。

    3種物源隨高程變化對泥沙的貢獻率有所不同。消落帶上方土壤對沉積泥沙的貢獻隨高程的增加從15.8%增加到30.3%,汝溪河上游懸移質(zhì)的貢獻率隨高程增加從38.4%減少至26.3%,而長江干流懸移質(zhì)對不同高程泥沙的貢獻率變化不大,介于43%~45.8%之間。

    不同高程消落帶沉積物來源及其貢獻差異主要與三峽水庫運行調(diào)度引起的消落帶淹水季節(jié)(圖2)有關(guān)。高程145~155 m受淹時間長達(dá)9個月,其中涵蓋了部分雨季(2016年4-6月)。雨季時支流上游來沙增加,洪水過程中低海拔消落帶(145~155 m)的季節(jié)性淹沒帶來的泥沙沉積致使其對該高程帶的泥沙貢獻較大,而>155~165和>165~175 m消落帶淹水時間均處旱季,因此支流懸移質(zhì)的貢獻有所降低。消落帶上部土壤侵蝕對不同高程的泥沙貢獻變化與支流懸移質(zhì)相反,這可能是由于消落帶上部土壤被侵蝕后在向下坡輸移的過程中發(fā)生了近源沉積,大部分泥沙優(yōu)先堆積在較高高程[26]。

    圖3 不同高程消落帶沉積泥沙及其顆粒態(tài)磷的物源相對貢獻率

    2.3 沉積泥沙磷素來源分析

    3種物源對不同高程消落帶沉積泥沙中顆粒態(tài)磷的相對貢獻比例見圖3b。隨著高程的變化,不同物源對磷的貢獻大小變化與泥沙貢獻趨勢基本一致,但相對比例有所不同,這與周慧平等[9]在九鄉(xiāng)河流域的研究結(jié)果相似。消落帶上方土壤(S1)和長江干流來沙(S2)對消落帶的磷素貢獻率隨高程的增加而增加,而汝溪河上游來沙(S3)的貢獻率則隨高程增加呈減少趨勢。

    不同高程消落帶沉積泥沙中顆粒態(tài)磷的來源貢獻主要受泥沙來源影響。長江干流和汝溪河上游來沙對不同高程帶顆粒態(tài)磷的貢獻率之和介于88.6%~95.4%,對應(yīng)的泥沙貢獻率之和為69.7%~84.2%,均占主導(dǎo)地位。長江干流來沙對顆粒態(tài)磷的貢獻率在>165~175 m高程帶達(dá)到最大(54.5%),汝溪河上游來沙對顆粒態(tài)磷的貢獻率在145~155 m高程帶達(dá)到最大(51.6%)。消落帶上方土壤對沉積物磷的貢獻較小,平均為8.5%,低于其泥沙平均貢獻率(25.1%),這是由于該源地土壤中磷的質(zhì)量分?jǐn)?shù)(0.45 g/kg)顯著低于干支流懸移質(zhì)中磷的質(zhì)量分?jǐn)?shù)(表2)。不同源地樣品磷素含量的差異可能與其顆粒粒徑有關(guān)。物源樣品顆粒分析結(jié)果表明(圖4),消落帶上方土壤的顆粒組成較干、支流泥沙偏粗,河流泥沙經(jīng)過一定距離的搬運分選后顆粒變細(xì),顆粒越細(xì)對磷素的吸附能力越強[27]。

    圖4 不同物源粒徑分布特征

    本文局限性分析:

    1)需要指出的是,不同于利用137Cs、210Pb等放射性核素來區(qū)分土壤表層和次表層泥沙來源時具有較為明確的理論基礎(chǔ)[28],基于地球化學(xué)元素的指紋技術(shù)目前很難建立特定指紋參數(shù)與環(huán)境因子之間的直接聯(lián)系,僅僅根據(jù)統(tǒng)計檢驗方法篩選得到的指紋組合可能具有較大的不確定性,不同的參數(shù)組合可能對泥沙來源的定量具有較大影響[29-30],因此應(yīng)結(jié)合其他方法進一步驗證指紋判別結(jié)果的可靠性。

    2)消落帶的侵蝕沉積過程非常復(fù)雜,受水文節(jié)律、地形地貌等多因素影響。本研究只基于1個蓄水周期的采樣分析,樣本量較少,具有一定的局限性;同時,本研究是在定量泥沙來源的基礎(chǔ)上分析磷的來源,因此重點關(guān)注被泥沙吸附的顆粒態(tài)磷,而忽略了泥沙輸移過程和蓄水期間可能發(fā)生的磷素吸附、解析以及形態(tài)轉(zhuǎn)化等。

    3 結(jié) 論

    不同高程,消落帶泥沙來源最佳指紋組合不同。不同指紋組合對源地樣品的正確判別效果均較好(>95.8%)。受近源侵蝕的影響,不同高程消落帶沉積泥沙物源相對貢獻率的差異主要表現(xiàn)在消落帶上方土壤和汝溪河上游懸移質(zhì),3個高程帶長江干流懸移質(zhì)的貢獻率均為最大(43%~45.8%)。隨著高程的變化,不同物源對磷的貢獻大小變化與泥沙貢獻趨勢基本一致。

    長江干流和汝溪河上游來沙是消落帶沉積泥沙及顆粒態(tài)磷的主要來源。受水庫運行調(diào)度引起的消落帶淹水時間影響,隨著高程的增加消落帶上方的岸坡侵蝕和長江干流來沙對顆粒態(tài)磷的貢獻相應(yīng)增多,而汝溪河支流對顆粒態(tài)磷的貢獻率隨著高程的增加呈減少趨勢。干支流泥沙顆粒組成較細(xì),對磷等污染物具有明顯的富集效應(yīng),將對消落帶及庫區(qū)水生態(tài)環(huán)境造成潛在風(fēng)險,而富含營養(yǎng)物質(zhì)的支流泥沙在消落帶堆積,可能隨暴雨侵蝕再次進入水體,加重污染負(fù)荷。蓄水后泥沙沉積引起的消落帶磷素長期富集可能會對水庫消落帶及水環(huán)境產(chǎn)生重要影響。

    [1] Bao Yuhai, Gao Peng, He Xiubin. The water-level fluctuation zone of Three Gorges Reservoir: A unique geomorphological unit[J]. Earth-Science Reviews, 2015, 150: 14-24.

    [2] 中華人民共和國環(huán)境保護部. 長江三峽工程生態(tài)與環(huán)境監(jiān)測公報[R]. 2014:27-29.

    [3] 中華人民共和國環(huán)境保護部. 長江三峽工程生態(tài)與環(huán)境監(jiān)測公報[R]. 2015:25-27.

    [4] Pionke H B, Kunishi H M. Phosphorus status and content of suspended sediment in a Pennsylvania watershed[J]. Soil Sci, 1992, 153: 452-62.

    [5] Vanni M J, Renwick W, Headworth J L, et al. Dissolved and particulate nutrient flux from three adjacent agricultural watersheds: A five-year study[J]. Biogeochemistry, 2001, 54: 85-114.

    [6] Foster I D L, Chapman A S, Hodgkinson R M, et al. Changing suspended sediment and particulate phosphorus loads and pathways in under drained lowland agricultural catchments: Herefordshire and Worcestershire, UK[J]. Hydrobiologia, 2003, 494: 119-126.

    [7] Simard R R, Beauchemin S, Haygarth P M. Potential for pathways of phosphorus transport[J]. J Environ Qual, 2000, 29(1): 97-105.

    [8] Walling D E, Collins A L, Stroud R W. Tracing suspended sediment and particulate phosphorus sources in catchments[J]. Journal of Hydrology, 2008, 350(3/4): 274-289.

    [9] 周慧平,常維娜,張龍江. 基于泥沙指紋識別的小流域顆粒態(tài)磷來源解析[J]. 農(nóng)業(yè)工程學(xué)報,2015,31(13):251-256.

    Zhou Huiping, Chang Weina, Zhang Longjiang. Particulate phosphorus sources apportionment in small watershed based on sediment fingerprinting[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(13): 251-256. (in Chinese with English abstract)

    [10] 楊明義,徐龍江. 黃土高原小流域泥沙來源的復(fù)合指紋識別法分析[J]. 水土保持學(xué)報,2010,24(2):30-34.

    Yang Mingyi, Xu Longjiang. Fingerprinting suspended sediment sources in a small catchment on the Loess Plateau[J]. Journal of Soil and Water Conservation, 2010, 24(2): 30-34. (in Chinese with English abstract)

    [11] 郭進,文安邦,嚴(yán)冬春,等. 復(fù)合指紋識別技術(shù)定量示蹤流域泥沙來源[J]. 農(nóng)業(yè)工程學(xué)報,2014,30(2):94-104.

    Guo Jin, Wen Anbang, Yan Dongchun, et al. Quantifying catchment scale sediment source using composite fingerprinting technique[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(2): 94-104. (in Chinese with English abstract)

    [12] Pulley S, Foster I, Antunes P. The uncertainties associated with sediment fingerprinting suspended and recently deposited fluvial sediment in the Nene river basin[J]. Geomorphology, 2015, 228: 303-319.

    [13] Chen Fangxin, Zhang Fengbao, Fang Nufang, et al. Sediment source analysis using the fingerprinting method in a small catchment of the Loess Plateau, China[J]. Journal of Soils and Sediments, 2016, 16(5): 1655-1669.

    [14] Walling D E, Woodward J C. Tracing sources of suspended sediment in river basins: A case study of the River Culm, Devon, UK[J]. Marine and Freshwater Research, 1995, 46: 327-336.

    [15] Haddadchi A, Olley J, Pietsch T. Quantifying sources of suspended sediment in three size fractions[J]. Journal of Soils and Sediments, 2015, 15(10): 2086-2100.

    [16] Vale S S, Fuller I C, Procter J N, et al. Characterization and quantification of suspended sediment sources to the Manawatu River, New Zealand[J]. Science of The Total Environment, 2016, 543:171-186.

    [17] Olley J, Caitcheon G. Major element chemistry of sediments from the Darling–Barwon river and its tributaries: implications for sediment and phosphorus sources[J]. Hydrological Processes, 2000, 14(7): 1159-1175.

    [18] Shi Zhonglin, Yan Dongchun, Wen Anbang, et al. Sedimentary Processes: Examples from Asia, Turkey and Nigeria[M]. London: IntechOpen, 2019.

    [19] Walling D E, Foster I D L. Using environmental radionuclides, mineral magnetism and sediment geochemistry for tracing and dating fluvial sediments[M]//Kondolf G M, Piégay H. Tools in Fluvial Geomorphology: Second Edition. Chichester: Wiley, 2016: 183-209.

    [20] Collins A L, Walling D E, Leeks G J L. Sediment sources in the upper severn catchment: A fingerprinting approach[J]. Hydrology & Earth System Sciences, 1997, 1(3): 509-521.

    [21] Koiter A J, Owens P N, Petticrew E L, et al. Assessment of particle size and organic matter correction factors in sediment source fingerprinting investigations: An example of two contrasting watershed on Canada[J]. Geoderma, 2018, 325: 195-207.

    [22] Smith H G, Blake W H. Sediment fingerprinting in agricultural catchments: A critical re-examination of source discrimination and data corrections[J]. Geomorphology, 2014, 204: 177-191.

    [23] Wilkinson S N, Hancock G J, Bartley R, et al. Using sediment tracing to assess processes and spatial patterns of erosion in grazed rangelands, Burdekin River basin, Australia[J]. Agriculture, Ecosystems & Environment, 2013, 180: 90-102.

    [24] Collins A L, Walling D E, Leeks G J L. Source type ascription for fluvial suspended sediment based on a quantitative composite fingerprinting technique[J]. Catena, 1997, 29: 1-27.

    [25] Motha J A, Wallbrink P J, Hairsine P B, et al. Determining the sources of suspended sediment in a forested catchment in southeastern Australia[J]. Water Resources Research, 2003, 39(3): 1056-1069.

    [26] 閻丹丹,鮑玉海,賀秀斌,等. 三峽水庫蓄水后長江干支流及消落帶泥沙顆粒特征分析[J]. 水土保持學(xué)報,2014,28(4):289-292,329.

    Yan Dandan, Bao Yuhai, He Xiubin, et al. Particle size characteristics of sediment in draw down area of upper Yangtze River and its major tributaries of Three Gorges Reservoir[J]. Journal of Soil and Water Conservation, 2014, 28(4): 289-292, 329. (in Chinese with English abstract)

    [27] 秦宇,王紫薇,韓超. 懸移質(zhì)泥沙粒徑對磷吸附的影響[J].中國給水排水,2017,33(7):80-83.

    Qin Yu, Wang Ziwei, Han Chao. Effect of different size of suspended sediment on adsorption of phosphorus[J]. China Water & Wastewater, 2017, 33(7): 80-83. (in Chinese with English abstract)

    [28] 張信寶,賀秀斌,文安邦,等. 川中丘陵區(qū)小流域泥沙來源的137Cs和210Pb雙同位素法研究[J]. 科學(xué)通報,2004,49(15):1537-1541.

    [29] 周慧平,周宏偉,陳玉東. 指紋技術(shù)識別泥沙來源:不確定性研究進展[J]. 土壤學(xué)報,2019,56(6):1279-1289.

    Zhou Huiping, Zhou Hongwei, Chen Yudong. Sediment source fingerprinting: Progress in uncertainty analysis[J]. Acta Pedologica Sinica, 2019, 56(6): 1279-1289. (in Chinese with English abstract)

    [30] 周慧平,陳玉東,常維娜. 指紋技術(shù)識別泥沙來源:進展與展望[J]. 水土保持學(xué)報,2018,32(5):4-10.

    Zhou Huiping, Chen Yudong, Chang Weina. Sediments source fingerprinting: Progresses and prospects[J]. Journal of Soil and Water Conservation, 2018, 32(5):4-10. (in Chinese with English abstract)

    Fingerprinting particulate phosphorus absorbed by sediments for riparian zone deposits in tributary of Three Gorges Reservoir

    Chen Taili1,2, Shi Zhonglin1※, Wang Yongyan1,2, Yan Dongchun1, Wen Anbang1, Chen Jiacun1,2

    (1.,,610041,; 2.,100049,)

    Information on deposited sediment and associated particulate phosphorus (PP) is of great importance for development of effective management strategies in the riparian zone of the Three Gorges Reservoir. In this study, a composite fingerprinting approach was used to identify the sources of deposited sediment and associated PP in the riparian zone with different elevations (145-155, >155-165, >165-175 m) along the Ruxi tributary channel. Three potential sediment sources were defined: the topsoil above the 175 m elevation level of the riparian zone, and the suspended sediments transported by both the mainstream of Yangtze River and the upstream Ruxi River. The target deposited sediments were collected at different elevation levels of the riparian zone along the Ruxi River and a total of 18 geochemical properties factors were analyzed for the size fraction of sediment smaller than 63m. A commonly used 3-step (range test, KW-H test and DFA analysis) statistical methodology to element selection was applied to select the optimum subset of tracer properties that could best discriminate sources. The relative contribution of the 3 potential sources to the deposited sediment in the riparian zone was estimated using a frequentist-based multivariate mixing model. The relative importance of those sources to the overall PP was then determined by combining the information obtained from deposited sediment sources with information on the P content of those sources. Source apportionment results showed that the suspended sediments transported by both the Yangtze mainstream and the upstream Ruxi tributary represent the dominant sources of PP deposits in the riparian zone, which account for 88.6%-95.4% of the PP input. The corresponding sediment contribution from these 2 sources ranged from 69.7% to 84.2%, which varied between elevation levels. The PP contribution from the Yangtze mainstream was 54.5% for the >165-175 m level. The relative PP contribution from the Ruxi River was the most important during the early wet season, which was 51.6% (145-155 m). In contrast, the PP input from the topsoil above the 175 m level of the riparian zone exhibited to be less important (mean of 8.5%) during the study period. The information on sediment contribution of this source reveals that most of the eroded sediment from the topsoil above the 175 m level was preferentially deposited in the elevation levels of <155-165 m and >165-175 m. The findings of this research highlight the important impacts of impoundment on the riparian zone of the tributaries of the reservoir in terms of sediment accumulation and pollution.

    phosphorus; sediments; elements; riparian zone; composite fingerprinting technique

    陳太麗,史忠林,王永艷,嚴(yán)冬春,文安邦,陳佳村. 三峽水庫典型支流消落帶泥沙顆粒態(tài)磷復(fù)合指紋示蹤研究[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(20):118-124.doi:10.11975/j.issn.1002-6819.2019.20.015 http://www.tcsae.org

    Chen Taili, Shi Zhonglin, Wang Yongyan, Yan Dongchun, Wen Anbang, Chen Jiacun. Fingerprinting particulate phosphorus absorbed by sediments for riparian zone deposits in tributary of Three Gorges Reservoir[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(20): 118-124. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.20.015 http://www.tcsae.org

    2019-04-02

    2019-09-10

    國家重點研發(fā)計劃“構(gòu)建源頭削減-生物隔離-濕地消納相結(jié)合的高效生態(tài)攔截技術(shù)體系”(2017YFD0800505);國家自然科學(xué)基金“三峽庫區(qū)支流消落帶土-水界面磷素遷移過程與通量”(41430750)

    陳太麗,博士生,主要從事土壤侵蝕與水土保持研究。Email:chentaili1994@126.com

    史忠林,助理研究員,博士。主要從事侵蝕泥沙示蹤技術(shù)研究。Email:shizl@imde.ac.cn

    10.11975/j.issn.1002-6819.2019.20.015

    O613.62

    A

    1002-6819(2019)-20-0118-07

    猜你喜歡
    懸移質(zhì)溪河物源
    人類活動影響下全球河流懸移質(zhì)泥沙通量快速變化研究
    安徽裕溪河流域發(fā)現(xiàn)的石制品研究
    東方考古(2020年0期)2020-11-06 05:34:02
    強震區(qū)泥石流物源演化指標(biāo)選取及規(guī)律分析
    崖門出海航道回淤分析
    珠江水運(2018年21期)2018-12-20 23:17:38
    頭屯河流域河流懸移質(zhì)泥沙分析
    美麗的泗溪河
    美麗的泗溪河
    南海北部陸架表層沉積物重礦物分布特征及物源意義
    懸移質(zhì)含沙量垂線分布
    水道港口(2014年1期)2014-04-27 14:14:35
    物源對子洲氣田山2段儲層的影響
    久久人妻福利社区极品人妻图片| 咕卡用的链子| 黄网站色视频无遮挡免费观看| 国产不卡一卡二| 极品教师在线免费播放| 午夜老司机福利片| 黄色丝袜av网址大全| 无限看片的www在线观看| 精品国产亚洲在线| 91麻豆精品激情在线观看国产 | 99国产精品一区二区蜜桃av | 热99久久久久精品小说推荐| 老司机影院毛片| 青草久久国产| 后天国语完整版免费观看| 久久99一区二区三区| 一个人免费在线观看的高清视频| 久久精品人人爽人人爽视色| 一个人免费在线观看的高清视频| 久久久久久久午夜电影 | 亚洲第一青青草原| 午夜精品在线福利| 国精品久久久久久国模美| 激情视频va一区二区三区| 欧美黑人精品巨大| 亚洲人成电影观看| 成人影院久久| 欧美日韩黄片免| 国产精品免费一区二区三区在线 | 高清欧美精品videossex| 老熟女久久久| 免费高清在线观看日韩| 国产亚洲欧美98| 一区二区三区国产精品乱码| videos熟女内射| 天天添夜夜摸| 可以免费在线观看a视频的电影网站| 在线观看舔阴道视频| 亚洲精品在线观看二区| 男女之事视频高清在线观看| 国产精品免费大片| 热99国产精品久久久久久7| 久久国产乱子伦精品免费另类| 俄罗斯特黄特色一大片| 自拍欧美九色日韩亚洲蝌蚪91| 国产在视频线精品| 91成人精品电影| 久久国产精品大桥未久av| 国产麻豆69| 国产1区2区3区精品| 黄色 视频免费看| 中文字幕人妻丝袜一区二区| 无人区码免费观看不卡| 欧美 亚洲 国产 日韩一| 午夜福利视频在线观看免费| 欧美精品人与动牲交sv欧美| 精品第一国产精品| 久久久精品区二区三区| 欧美老熟妇乱子伦牲交| 亚洲五月婷婷丁香| 精品高清国产在线一区| 丝瓜视频免费看黄片| 女性生殖器流出的白浆| 99久久精品国产亚洲精品| 91在线观看av| 日本五十路高清| 丝袜美腿诱惑在线| 美女扒开内裤让男人捅视频| 妹子高潮喷水视频| 露出奶头的视频| x7x7x7水蜜桃| 亚洲一区中文字幕在线| 欧美黑人欧美精品刺激| 精品国产一区二区久久| 中文字幕精品免费在线观看视频| 人妻 亚洲 视频| 啪啪无遮挡十八禁网站| 黄色a级毛片大全视频| 每晚都被弄得嗷嗷叫到高潮| 久久中文字幕人妻熟女| 国产精品 欧美亚洲| 激情在线观看视频在线高清 | 一级毛片女人18水好多| 免费看a级黄色片| 少妇裸体淫交视频免费看高清 | xxx96com| av天堂在线播放| 精品国产一区二区三区久久久樱花| aaaaa片日本免费| 天天躁日日躁夜夜躁夜夜| 亚洲美女黄片视频| 中亚洲国语对白在线视频| 亚洲欧洲精品一区二区精品久久久| 一级黄色大片毛片| 视频区欧美日本亚洲| 极品人妻少妇av视频| 咕卡用的链子| 他把我摸到了高潮在线观看| 日韩有码中文字幕| 黑人操中国人逼视频| 欧美精品高潮呻吟av久久| av网站免费在线观看视频| 男女高潮啪啪啪动态图| 亚洲全国av大片| 亚洲人成77777在线视频| 女人精品久久久久毛片| 国产精品免费视频内射| 日本精品一区二区三区蜜桃| 99热网站在线观看| 久久国产精品影院| videosex国产| 欧美中文综合在线视频| 身体一侧抽搐| 一夜夜www| 日韩欧美在线二视频 | 美女 人体艺术 gogo| 亚洲一区二区三区不卡视频| 欧美午夜高清在线| 久久久精品国产亚洲av高清涩受| 女人爽到高潮嗷嗷叫在线视频| 日本五十路高清| а√天堂www在线а√下载 | 欧美日韩亚洲综合一区二区三区_| 狂野欧美激情性xxxx| 97人妻天天添夜夜摸| 久久中文看片网| 一级作爱视频免费观看| av不卡在线播放| 成年动漫av网址| 又大又爽又粗| netflix在线观看网站| 国产成人系列免费观看| 青草久久国产| 亚洲自偷自拍图片 自拍| 国产精品电影一区二区三区 | 国产精品 欧美亚洲| 在线观看www视频免费| 久久香蕉国产精品| 国产激情久久老熟女| 亚洲成人免费电影在线观看| 丝袜美腿诱惑在线| 成人18禁高潮啪啪吃奶动态图| 国产成人欧美| 亚洲国产精品一区二区三区在线| 久久热在线av| 日韩成人在线观看一区二区三区| 桃红色精品国产亚洲av| 亚洲av成人一区二区三| 国产精品国产av在线观看| 欧美成狂野欧美在线观看| 热re99久久国产66热| 日本欧美视频一区| 国产熟女午夜一区二区三区| 亚洲一区中文字幕在线| 国产一区二区激情短视频| 色婷婷av一区二区三区视频| 曰老女人黄片| 亚洲中文av在线| 国产成人欧美在线观看 | 视频区图区小说| 欧美人与性动交α欧美软件| 亚洲久久久国产精品| 精品无人区乱码1区二区| 99国产精品一区二区三区| 久久精品国产清高在天天线| 性少妇av在线| 亚洲中文字幕日韩| 久久人妻熟女aⅴ| 久久久久精品国产欧美久久久| 一进一出抽搐gif免费好疼 | 国产精品永久免费网站| 亚洲精品在线观看二区| 国产精华一区二区三区| 欧美精品啪啪一区二区三区| 久久久水蜜桃国产精品网| 最新在线观看一区二区三区| 999精品在线视频| 久久人妻福利社区极品人妻图片| 国产成人av激情在线播放| 女人高潮潮喷娇喘18禁视频| 亚洲精品自拍成人| 搡老熟女国产l中国老女人| 久热这里只有精品99| 在线观看日韩欧美| 久热爱精品视频在线9| 五月开心婷婷网| 亚洲成国产人片在线观看| 在线观看免费视频日本深夜| 国产免费现黄频在线看| 国产激情久久老熟女| 国产高清视频在线播放一区| av超薄肉色丝袜交足视频| 啦啦啦在线免费观看视频4| 一级a爱片免费观看的视频| 桃红色精品国产亚洲av| 成年人午夜在线观看视频| 国产精品二区激情视频| 亚洲国产中文字幕在线视频| 日韩欧美国产一区二区入口| 久久精品亚洲av国产电影网| 美女国产高潮福利片在线看| 久久国产精品大桥未久av| 男女高潮啪啪啪动态图| 精品人妻1区二区| 久久久久国内视频| 热99国产精品久久久久久7| 91精品三级在线观看| 亚洲成国产人片在线观看| 热re99久久精品国产66热6| 狠狠婷婷综合久久久久久88av| 国产精品一区二区在线不卡| 国产一区在线观看成人免费| 亚洲欧美激情综合另类| 国产亚洲av高清不卡| 久久国产精品大桥未久av| 中文字幕最新亚洲高清| 91成人精品电影| 男人的好看免费观看在线视频 | 午夜福利视频在线观看免费| 女警被强在线播放| 久久久精品区二区三区| 色精品久久人妻99蜜桃| 亚洲自偷自拍图片 自拍| 亚洲精品美女久久久久99蜜臀| 成人三级做爰电影| 99久久人妻综合| 日韩大码丰满熟妇| 欧美日韩精品网址| 免费观看人在逋| 怎么达到女性高潮| 他把我摸到了高潮在线观看| 人妻久久中文字幕网| 国产精品综合久久久久久久免费 | av视频免费观看在线观看| 久久这里只有精品19| 少妇 在线观看| 叶爱在线成人免费视频播放| 亚洲 欧美一区二区三区| 成年动漫av网址| 精品一区二区三区四区五区乱码| 欧美亚洲日本最大视频资源| 日韩大码丰满熟妇| 成年版毛片免费区| 90打野战视频偷拍视频| av国产精品久久久久影院| 精品一品国产午夜福利视频| 老司机午夜福利在线观看视频| 欧美 日韩 精品 国产| 欧美精品人与动牲交sv欧美| 国产免费男女视频| 国产成人精品在线电影| av国产精品久久久久影院| 国产主播在线观看一区二区| 老司机午夜福利在线观看视频| 国产高清视频在线播放一区| 欧美精品一区二区免费开放| 侵犯人妻中文字幕一二三四区| 欧美激情高清一区二区三区| 色播在线永久视频| 久久久水蜜桃国产精品网| 国内久久婷婷六月综合欲色啪| 亚洲美女黄片视频| 99久久精品国产亚洲精品| 999久久久精品免费观看国产| 99国产精品一区二区蜜桃av | 一级片免费观看大全| 女人爽到高潮嗷嗷叫在线视频| 99国产精品免费福利视频| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久av美女十八| 夜夜夜夜夜久久久久| 老司机靠b影院| 亚洲国产欧美一区二区综合| 热99国产精品久久久久久7| 女性生殖器流出的白浆| 国产激情欧美一区二区| 99国产精品免费福利视频| 国产精品欧美亚洲77777| 新久久久久国产一级毛片| 后天国语完整版免费观看| 久热爱精品视频在线9| 男女之事视频高清在线观看| 国产成+人综合+亚洲专区| 美女高潮到喷水免费观看| 久久人妻av系列| 精品视频人人做人人爽| 大片电影免费在线观看免费| 国产又爽黄色视频| 99久久精品国产亚洲精品| 少妇粗大呻吟视频| 香蕉国产在线看| 亚洲欧美精品综合一区二区三区| 夜夜爽天天搞| 欧美不卡视频在线免费观看 | 69av精品久久久久久| 久久精品国产综合久久久| 亚洲久久久国产精品| 大码成人一级视频| 国产一区有黄有色的免费视频| 精品国产一区二区久久| www.自偷自拍.com| 日韩有码中文字幕| 18禁观看日本| 欧美日韩一级在线毛片| 美女视频免费永久观看网站| 在线观看免费日韩欧美大片| 亚洲五月色婷婷综合| 狠狠狠狠99中文字幕| 亚洲中文日韩欧美视频| 女人久久www免费人成看片| 久久久久久久久免费视频了| 美女高潮喷水抽搐中文字幕| 成人18禁高潮啪啪吃奶动态图| 午夜福利乱码中文字幕| 色综合欧美亚洲国产小说| 精品人妻1区二区| 操出白浆在线播放| 久久久久久久国产电影| 99香蕉大伊视频| 欧美成人午夜精品| 亚洲五月婷婷丁香| 国产99白浆流出| a在线观看视频网站| 欧美精品亚洲一区二区| 日韩人妻精品一区2区三区| 女性被躁到高潮视频| 一个人免费在线观看的高清视频| 99久久99久久久精品蜜桃| 国产精品.久久久| 亚洲精品久久午夜乱码| 亚洲精品美女久久av网站| 国内久久婷婷六月综合欲色啪| 日本欧美视频一区| 欧美日韩亚洲国产一区二区在线观看 | 日韩精品免费视频一区二区三区| 成熟少妇高潮喷水视频| 99re在线观看精品视频| 国产激情久久老熟女| 国产精品久久久av美女十八| 亚洲综合色网址| 色在线成人网| 国产97色在线日韩免费| 国产在线观看jvid| 一级片'在线观看视频| 日韩欧美国产一区二区入口| 老司机午夜十八禁免费视频| www.自偷自拍.com| 精品无人区乱码1区二区| 亚洲 国产 在线| 18在线观看网站| 久久国产精品影院| 国产成人影院久久av| 国产精品1区2区在线观看. | 在线国产一区二区在线| 黑丝袜美女国产一区| 黑人操中国人逼视频| 国产精品一区二区在线观看99| 看片在线看免费视频| 一二三四社区在线视频社区8| 成人影院久久| 欧美激情 高清一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 黄片大片在线免费观看| 亚洲精品中文字幕在线视频| 中文字幕人妻熟女乱码| 国产亚洲精品久久久久久毛片 | 亚洲精品在线美女| 久久午夜亚洲精品久久| 久久婷婷成人综合色麻豆| 欧美色视频一区免费| 国产成人欧美| av福利片在线| 亚洲精品在线观看二区| 国产亚洲欧美98| 免费黄频网站在线观看国产| 精品久久蜜臀av无| 精品一品国产午夜福利视频| 国产人伦9x9x在线观看| 亚洲熟妇中文字幕五十中出 | 免费在线观看影片大全网站| 久久久国产成人免费| 国产成人欧美| 国产麻豆69| 亚洲精品在线观看二区| 久久草成人影院| 国产又爽黄色视频| 成人18禁在线播放| 欧美黄色片欧美黄色片| 看免费av毛片| 亚洲av电影在线进入| 91在线观看av| 最近最新中文字幕大全免费视频| 一级a爱视频在线免费观看| 国产黄色免费在线视频| 国产无遮挡羞羞视频在线观看| 一边摸一边抽搐一进一小说 | 成人手机av| 免费久久久久久久精品成人欧美视频| 欧美日韩av久久| 国产精品九九99| 在线看a的网站| 国产成人精品久久二区二区免费| 天天添夜夜摸| 亚洲一区二区三区不卡视频| 19禁男女啪啪无遮挡网站| xxx96com| 亚洲成人国产一区在线观看| 国内久久婷婷六月综合欲色啪| 国产在线观看jvid| 免费久久久久久久精品成人欧美视频| 18禁国产床啪视频网站| 久久热在线av| www.999成人在线观看| 亚洲情色 制服丝袜| 两个人免费观看高清视频| 色尼玛亚洲综合影院| 另类亚洲欧美激情| 亚洲情色 制服丝袜| 成年人免费黄色播放视频| xxxhd国产人妻xxx| 91成年电影在线观看| 王馨瑶露胸无遮挡在线观看| 欧美日韩乱码在线| 色婷婷久久久亚洲欧美| 中亚洲国语对白在线视频| 又紧又爽又黄一区二区| 欧美精品一区二区免费开放| 午夜精品国产一区二区电影| 日本vs欧美在线观看视频| 亚洲美女黄片视频| 夫妻午夜视频| 美女高潮喷水抽搐中文字幕| 免费观看精品视频网站| 人人妻人人澡人人爽人人夜夜| 91在线观看av| 在线看a的网站| 少妇 在线观看| 久久久久久久久久久久大奶| 亚洲中文字幕日韩| 日韩有码中文字幕| 午夜免费成人在线视频| 91麻豆av在线| 国产不卡av网站在线观看| 欧美久久黑人一区二区| 欧美激情极品国产一区二区三区| 亚洲午夜理论影院| 亚洲国产精品一区二区三区在线| 一级黄色大片毛片| 精品熟女少妇八av免费久了| 欧美精品人与动牲交sv欧美| 女同久久另类99精品国产91| 大香蕉久久成人网| 999久久久国产精品视频| 很黄的视频免费| 久久国产乱子伦精品免费另类| 国产欧美日韩一区二区精品| 悠悠久久av| 三级毛片av免费| 午夜影院日韩av| 久久久国产一区二区| 午夜福利一区二区在线看| 欧美久久黑人一区二区| 亚洲人成77777在线视频| 又黄又粗又硬又大视频| 中亚洲国语对白在线视频| 一区二区三区激情视频| 满18在线观看网站| 真人做人爱边吃奶动态| 日韩精品免费视频一区二区三区| 久9热在线精品视频| 好看av亚洲va欧美ⅴa在| 91麻豆精品激情在线观看国产 | 丝袜人妻中文字幕| 99精品在免费线老司机午夜| 国产一卡二卡三卡精品| 日本撒尿小便嘘嘘汇集6| 老司机亚洲免费影院| 99re6热这里在线精品视频| 亚洲精品久久午夜乱码| 极品教师在线免费播放| 国产精品98久久久久久宅男小说| 在线观看日韩欧美| 日韩成人在线观看一区二区三区| svipshipincom国产片| 男女床上黄色一级片免费看| 黄色毛片三级朝国网站| 少妇裸体淫交视频免费看高清 | 亚洲国产欧美日韩在线播放| 色婷婷av一区二区三区视频| 精品一品国产午夜福利视频| 91老司机精品| 天天添夜夜摸| 一级黄色大片毛片| 美女 人体艺术 gogo| 日韩 欧美 亚洲 中文字幕| 久久人人爽av亚洲精品天堂| 午夜福利在线免费观看网站| 手机成人av网站| 久久九九热精品免费| 日韩欧美国产一区二区入口| 免费av中文字幕在线| www.999成人在线观看| 香蕉久久夜色| 在线观看午夜福利视频| 亚洲全国av大片| 老司机午夜福利在线观看视频| 精品国产美女av久久久久小说| 亚洲 欧美一区二区三区| 精品国产乱码久久久久久男人| 水蜜桃什么品种好| 午夜精品久久久久久毛片777| 男女午夜视频在线观看| av天堂久久9| 日韩熟女老妇一区二区性免费视频| av不卡在线播放| a在线观看视频网站| 久久久久久人人人人人| 大码成人一级视频| 免费在线观看完整版高清| 热re99久久国产66热| a级毛片在线看网站| 亚洲全国av大片| 国产成人啪精品午夜网站| 老汉色av国产亚洲站长工具| 欧美乱色亚洲激情| 欧美亚洲 丝袜 人妻 在线| 自线自在国产av| 欧美国产精品一级二级三级| 大香蕉久久网| 国产乱人伦免费视频| 久久久国产成人免费| 亚洲国产精品sss在线观看 | 午夜福利影视在线免费观看| 日韩大码丰满熟妇| 亚洲伊人色综图| 黑人猛操日本美女一级片| 亚洲精品粉嫩美女一区| 亚洲国产精品合色在线| 欧美色视频一区免费| 国产真人三级小视频在线观看| 国产在线精品亚洲第一网站| 欧美日韩亚洲综合一区二区三区_| 侵犯人妻中文字幕一二三四区| 国产精品美女特级片免费视频播放器 | 亚洲一区高清亚洲精品| 国产一卡二卡三卡精品| 一本大道久久a久久精品| 成年人黄色毛片网站| 精品久久久久久,| xxx96com| 男女免费视频国产| 久久婷婷成人综合色麻豆| 可以免费在线观看a视频的电影网站| 久久中文看片网| 天天躁日日躁夜夜躁夜夜| 又大又爽又粗| 亚洲少妇的诱惑av| 国产精品一区二区在线不卡| 免费观看人在逋| 丰满饥渴人妻一区二区三| 十分钟在线观看高清视频www| 久久精品国产清高在天天线| 午夜91福利影院| 国产成人影院久久av| 国产男靠女视频免费网站| 欧美大码av| 亚洲九九香蕉| 在线观看www视频免费| 18禁裸乳无遮挡动漫免费视频| 在线免费观看的www视频| 久久精品亚洲精品国产色婷小说| 久久精品亚洲av国产电影网| 欧美在线黄色| 国产精品电影一区二区三区 | 亚洲精品中文字幕在线视频| av网站在线播放免费| 国产欧美亚洲国产| 亚洲专区国产一区二区| 很黄的视频免费| 欧美黄色片欧美黄色片| 欧美 亚洲 国产 日韩一| 欧美日韩视频精品一区| 欧美黑人精品巨大| 看黄色毛片网站| 大香蕉久久成人网| 精品国产一区二区三区四区第35| 99国产精品99久久久久| 十八禁人妻一区二区| 在线观看一区二区三区激情| 脱女人内裤的视频| av网站免费在线观看视频| 伊人久久大香线蕉亚洲五| 黄色丝袜av网址大全| 狠狠婷婷综合久久久久久88av| 国产无遮挡羞羞视频在线观看| 精品午夜福利视频在线观看一区| 视频区欧美日本亚洲| 国产亚洲av高清不卡| 中文字幕色久视频| 国产欧美亚洲国产| 黑人巨大精品欧美一区二区mp4| 操美女的视频在线观看| 人妻一区二区av| 下体分泌物呈黄色| 亚洲人成伊人成综合网2020| 黄片小视频在线播放| 欧美性长视频在线观看| 国产成人一区二区三区免费视频网站| 国产欧美日韩一区二区三| 国产在视频线精品| 免费观看a级毛片全部| 国产高清激情床上av| 久久精品国产综合久久久| 一区在线观看完整版| av欧美777|