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

    三峽庫區(qū)典型流域水質(zhì)時(shí)空特征及污染防控策略

    2023-01-16 09:36:52李剛浩范先鵬劉宏斌吳茂前張富林翟麗梅周繼文孔祥瓊程子珍
    關(guān)鍵詞:坡耕地氮磷林地

    李剛浩,范先鵬,夏 穎,劉宏斌,吳茂前,張富林,黃 敏,翟麗梅,周繼文,孔祥瓊,程子珍

    三峽庫區(qū)典型流域水質(zhì)時(shí)空特征及污染防控策略

    李剛浩1,2,范先鵬1,夏 穎1※,劉宏斌3,吳茂前1,張富林1,黃 敏2,翟麗梅3,周繼文4,孔祥瓊4,程子珍1

    (1. 湖北省農(nóng)業(yè)科學(xué)院植保土肥研究所/國家農(nóng)業(yè)環(huán)境潛江觀測實(shí)驗(yàn)站/湖北省農(nóng)業(yè)面源污染防治工程技術(shù)研究中心,武漢 430064;2. 武漢理工大學(xué)資源與環(huán)境工程學(xué)院,武漢 430070;3. 中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,北京 100081;4. 興山縣農(nóng)業(yè)農(nóng)村局,興山 443700)

    為了解三峽庫區(qū)典型小流域不同土地利用類型下水質(zhì)的變化特征,于2015—2020年對(duì)三峽庫區(qū)長坪小流域開展水質(zhì)監(jiān)測,對(duì)不同土地利用類型下水體中的總氮(Total Nitrogen,TN)、硝態(tài)氮(NO3--N)、氨氮(NH4+-N)、總磷(Total Phosphorus,TP)、可溶態(tài)總磷(Dissolved Total Phosphorus,DTP)和顆粒態(tài)總磷(Particulate Phosphorus,PP)濃度的時(shí)空變化規(guī)律進(jìn)行研究,識(shí)別流域內(nèi)主要氮磷污染源,從而提出針對(duì)性的污染防控策略。結(jié)果表明:1)不同土地利用類型氮磷輸出濃度從大到小依次為村莊、坡耕地、林地和水庫,其TN平均濃度分別為8.29、2.88、1.57和1.43 mg/L,TP平均濃度分別為0.25、0.13、0.09和0.07 mg/L。2)不同土地利用類型的水質(zhì)在汛期和非汛期存在差異,村莊氮磷輸出濃度為非汛期高于汛期,坡耕地、林地和水庫則表現(xiàn)為汛期高于非汛期。3)內(nèi)梅羅綜合污染指數(shù)評(píng)價(jià)結(jié)果表明林地和水庫的污染程度均為清潔水平,坡耕地污染程度為輕污染水平,村莊受散養(yǎng)生豬數(shù)量的影響,污染程度從污染水平轉(zhuǎn)變?yōu)橹匚廴舅健?)長坪小流域TN和TP年輸出負(fù)荷總量分別為4 278.59和364.93 kg/a。畜禽養(yǎng)殖源是小流域氮磷最主要的污染源,其TN和TP的輸出負(fù)荷分別占流域總負(fù)荷的45.69%和71.77%。坡耕地的TN和TP單位面積污染負(fù)荷分別是村莊的7.58%和1.79%,與村莊相比,坡耕地具有顯著的低污染特征,因此,以坡耕地消納村莊污水和畜禽糞便,可促進(jìn)流域內(nèi)糞污的就地消納,達(dá)到多源共治的效果?;谛×饔虿煌恋乩妙愋偷乃|(zhì)特征和生態(tài)系統(tǒng)服務(wù)功能,將流域劃分為林草水源涵養(yǎng)區(qū)、村莊污染控制區(qū)和坡耕地水土保持區(qū),并進(jìn)行分區(qū)協(xié)同防控,促進(jìn)流域農(nóng)業(yè)面源污染系統(tǒng)高效治理。

    水質(zhì);土地利用;空間分布;三峽庫區(qū);典型流域;農(nóng)業(yè)面源污染;分區(qū)協(xié)同防控策略

    0 引 言

    農(nóng)業(yè)生產(chǎn)活動(dòng)和農(nóng)村生活產(chǎn)生的氮(N)和磷(P)是造成水體富營養(yǎng)化的重要因素[1]。近年來,隨著點(diǎn)源治理力度的增大,由農(nóng)業(yè)生產(chǎn)活動(dòng)和農(nóng)村生活引起的農(nóng)業(yè)面源污染已成為目前水體氮磷污染的最主要來源[2]。因此,研究不同土地利用類型對(duì)流域水質(zhì)變化的作用,對(duì)流域的污染防控和居民居住環(huán)境改善具有重要的意義。

    土地利用類型作為水質(zhì)變化的主要因素之一,隨著時(shí)空尺度的變化對(duì)河流水質(zhì)產(chǎn)生不同的影響[3]。閔惠琳等[4]的研究發(fā)現(xiàn)各形態(tài)氮的濃度均表現(xiàn)出村莊最高,其次為梯田和河流,林地濃度最低的分布特征;王曉燕等[5]對(duì)石匣小流域不同土地利用類型中氮磷流失特征的研究發(fā)現(xiàn),徑流中氮磷的濃度排序?yàn)榇迩f>坡耕地>林果地>荒草地。林地和草地能夠攔蓄污染物,改善水質(zhì),建設(shè)用地和坡耕地導(dǎo)致水質(zhì)惡化[6-7]。彭勃等[8]對(duì)清源河小流域的研究發(fā)現(xiàn),耕地和建設(shè)用地面積越大,臨近水體污染越嚴(yán)重,林地和草地則對(duì)水體污染物輸出具有抑制作用。吳東等[9]在蘭陵溪小流域的研究發(fā)現(xiàn),氮素輸出濃度與居民區(qū)和坡耕地正相關(guān),與林地負(fù)相關(guān),居民區(qū)和坡耕地是主要的氮素輸出源。季節(jié)性的變化對(duì)流域水質(zhì)有重要影響,華玲玲等[10]對(duì)古夫河小流域的研究表明,流域氮磷輸出濃度在豐水期最高,梁爽等[11]對(duì)鄂西長江小流域的研究發(fā)現(xiàn)氮素濃度呈現(xiàn)春夏高,秋冬低的特征。然而,不同流域的土地利用類型不同,受自然和人類活動(dòng)的影響,土地利用和水質(zhì)之間的關(guān)系并不總是一致,流域水質(zhì)的空間和時(shí)間分布特征仍不確定。

    三峽工程是世界上最大的水利樞紐工程之一,在長江流域經(jīng)濟(jì)開發(fā)與水文治理中發(fā)揮了重要作用,取得了顯著的經(jīng)濟(jì)效益和社會(huì)效益,但同時(shí)也對(duì)生態(tài)環(huán)境造成了影響[12]。三峽工程的建成,使庫區(qū)的土地利用類型發(fā)生變化,增加了對(duì)陡坡地的開發(fā)利用,且?guī)靺^(qū)人口稠密,土地利用強(qiáng)度大,再加上不合理的畜禽養(yǎng)殖和農(nóng)事操作,庫區(qū)農(nóng)業(yè)面源污染問題缺少有效處理[13]。三峽庫區(qū)水土流失嚴(yán)重,大量氮磷營養(yǎng)物質(zhì)隨地表徑流和泥沙流失而流失,導(dǎo)致三峽庫區(qū)水質(zhì)下降和水體的富營養(yǎng)化[14-15]。因此有必要深入了解該地區(qū)不同土地利用類型對(duì)水質(zhì)的影響,根據(jù)不同土地利用類型的水質(zhì)特征,提出系統(tǒng)高效的污染防控策略。

    小流域是一個(gè)相對(duì)獨(dú)立的系統(tǒng),是污染產(chǎn)生的源頭,可以作為面源污染的最佳研究單元[16]。不同土地利用類型水質(zhì)特征具有差異性,對(duì)不同土地利用類型的農(nóng)業(yè)面源污染防控,需要制定針對(duì)性的防控策略。因此,本研究選取三峽庫區(qū)上游具有較好獨(dú)立性和封閉性的長坪小流域作為研究對(duì)象,于2015—2020年在長坪小流域進(jìn)行長期連續(xù)的水質(zhì)監(jiān)測,結(jié)合GIS技術(shù),探討流域內(nèi)不同土地利用類型對(duì)水質(zhì)變化的影響,并識(shí)別流域內(nèi)氮磷的主要污染源,基于不同土地利用類型水質(zhì)的保護(hù)與治理需求,提出流域農(nóng)業(yè)面源污染的分區(qū)協(xié)同防控策略,以期為三峽庫區(qū)氮磷污染的系統(tǒng)高效防控、水環(huán)境改善、土地利用方式優(yōu)化提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)域概況

    研究區(qū)域長坪小流域(110°44'21″~110°46'17″E,31°18'30″~31°20'1″N)位于湖北省宜昌市興山縣,總面積330 hm2,是典型的農(nóng)林復(fù)合小流域。流域氣候?yàn)閬啛釒Ъ撅L(fēng)濕潤氣候,年平均氣溫14.9~18.5 ℃,年均降雨量1 000~1 300 mm,降雨季節(jié)分布不均,降雨主要集中在4-9月,具有夏季暴雨多、強(qiáng)度大、歷時(shí)短的特點(diǎn),根據(jù)當(dāng)?shù)貙?shí)際降雨情況將4-9月確定為豐水期,10月至次年3月確定為枯水期。小流域地形為兩邊高、中間低,屬于典型丘陵山區(qū),集水區(qū)徑流匯聚于唯一出口后流入古夫河,再匯入香溪河。主要土壤類型為黃棕壤,土層淺薄,抗蝕性差。主要土地利用類型為林地、坡耕地和村莊,其中,林地面積最大,占流域總面積的64.11%,坡耕地占流域總面積的28.73%,村莊占流域總面積的5.17%,具有村莊內(nèi)人畜混居普遍、地塊高度分散等特征。主要種植作物為玉米、油菜、馬鈴薯等,主要種植模式為玉米-油菜輪作。該小流域在地形地貌、氣候、土地利用、土壤類型和種植模式方面具有典型的三峽庫區(qū)特征。

    1.2 監(jiān)測點(diǎn)布置

    長坪小流域村莊、坡耕地、林地和水庫的土地面積之和占長坪小流域總面積的99%以上,因此在小流域自上而下選取11處水質(zhì)監(jiān)測點(diǎn)分別代表村莊、坡耕地、林地和水庫對(duì)小流域氮磷污染分布的影響(圖1)。其中,村莊、坡耕地、林地選取各自匯水范圍內(nèi)的支流作為采樣點(diǎn),采樣點(diǎn)位于支流與干流交匯的上方,水庫采樣點(diǎn)位于水庫出水口。

    圖1 流域土地利用類型及采樣點(diǎn)分布

    1.3 水樣采集與分析

    水樣采集時(shí)間為2015年1月-2020年12月,采集頻率為每5d采集一次,采樣位置位于河道中央,在所選斷面采集0.5 m以上表層水混合水樣(至少5點(diǎn)混合),以保證水質(zhì)樣品的代表性。使用500 mL聚乙烯瓶采集水樣,水樣采集后放置于冰箱冷凍保存,并及時(shí)送回實(shí)驗(yàn)室分析測試。測定指標(biāo)為總氮(TN)、硝態(tài)氮(NO3--N)、氨氮(NH4+-N)、總磷(TP)和可溶態(tài)總磷(DTP)。TN采用堿性過硫酸鉀氧化—紫外分光光度法測定(HJ636-2012),NO3--N采用紫外分光光度法測定(HJ/T346-2007),NH4+-N采用靛酚藍(lán)比色法測定,TP采用過硫酸鉀氧化—鉬藍(lán)比色法(GB11893-1989)測定,水樣經(jīng)0.45m微孔濾膜過濾后,測定DTP,測定方法與TP相同;顆粒態(tài)總磷(PP)濃度為TP濃度與DTP濃度之差。

    1.4 研究方法

    1.4.1 單因子污染指數(shù)法

    單因子污染指數(shù)法[17]是用來評(píng)價(jià)單項(xiàng)水質(zhì)指標(biāo)是否達(dá)到水體功能標(biāo)準(zhǔn)及其超標(biāo)程度。計(jì)算公式為

    式中P為單因子污染指數(shù);D為該斷面第項(xiàng)污染指標(biāo)的年平均值;D0為第項(xiàng)污染指標(biāo)的評(píng)價(jià)標(biāo)準(zhǔn)值。

    1.4.2 內(nèi)梅羅污染指數(shù)法

    內(nèi)梅羅污染指數(shù)法最早是由Nemerow[18]提出的,將某單項(xiàng)水質(zhì)指標(biāo)的實(shí)測濃度與地表水環(huán)境質(zhì)量中的標(biāo)準(zhǔn)限值進(jìn)行對(duì)比,計(jì)算出綜合指數(shù),并根據(jù)分級(jí)標(biāo)準(zhǔn)評(píng)價(jià)水體污染程度。內(nèi)梅羅污染指數(shù)兼顧了單因子污染指數(shù)的平均值和最大值,且突出最大污染指標(biāo)對(duì)水體污染的影響,更能反映水體的實(shí)際污染程度[19]。具體公式如下:

    式中為內(nèi)梅羅指數(shù),max是P的最大值,PP的平均值。

    1.4.3 輸出系數(shù)模型法

    輸出系數(shù)模型是通過對(duì)研究區(qū)域內(nèi)不同污染源類型的污染負(fù)荷進(jìn)行求和,最后得到該區(qū)域的污染總負(fù)荷,其模型方程[20]如下:

    式中表示污染物的總輸出負(fù)荷,kg/a;E表示第種污染源的輸出系數(shù),kg/(hm2·a)或kg/(人·a)或kg/(頭·a);A表示第種畜禽的數(shù)量(頭)、農(nóng)村人口數(shù)(人)或第種土地利用類型的面積,hm2。

    農(nóng)業(yè)面源污染來源主要包括農(nóng)村生活源、畜禽養(yǎng)殖源和種植業(yè)源,農(nóng)業(yè)面源污染輸出系數(shù)受農(nóng)村生活方式、養(yǎng)殖模式、種植模式和土地利用類型的影響而有所不同[21]。因此,本研究在綜合考慮小流域農(nóng)村生產(chǎn)生活方式的基礎(chǔ)上,借鑒三峽庫區(qū)農(nóng)業(yè)面源污染相關(guān)研究[22-25],確定各污染源的輸出系數(shù)。具體輸出系數(shù)取值見表1。

    表1 各污染源輸出系數(shù)

    1.4.4 污染防控分區(qū)劃分

    參照《全國生態(tài)功能區(qū)劃(修編版)》劃分方法,從流域面源污染治理角度,根據(jù)流域內(nèi)不同土地利用類型水質(zhì)污染程度和生態(tài)系統(tǒng)服務(wù)功能差異,將流域內(nèi)的土地利用類型劃分為不同的防控分區(qū)。

    1.4.5 數(shù)據(jù)分析

    采用Excel 2020進(jìn)行數(shù)據(jù)處理,采用Origin 2018繪圖,采用ArcGIS 10.2繪制流域土地利用類型及采樣點(diǎn)分布圖。

    2 結(jié)果與分析

    2.1 水質(zhì)的空間變化特征

    2015-2020年長坪小流域不同土地利用類型的水質(zhì)特征統(tǒng)計(jì)分析見表2。TN平均濃度從高到低依次為村莊、坡耕地、林地、水庫,村莊TN濃度最高,平均濃度為8.29 mg/L,坡耕地TN平均濃度次之,為2.88 mg/L,林地和水庫TN平均濃度較低,分別為1.57和1.43 mg/L。NO3--N輸出濃度從高到低依次為村莊、坡耕地、林地、水庫,NO3--N是氮素流失的主要形態(tài),占TN的55.24%~81.60%。NH4+-N平均濃度變化規(guī)律與TN和NO3--N不同,從高到低依次為村莊、水庫、林地、坡耕地,不同土地利用類型的NH4+-N濃度均較低,村莊、水庫、林地、坡耕地NH4+-N平均濃度分別為0.51、0.22、0.15、0.11 mg/L。TP平均濃度從高到低依次為村莊、坡耕地、林地、水庫,平均濃度分別為0.25、0.13、0.09和0.07 mg/L。綜上,長坪小流域各土地利用類型中,村莊的氮磷輸出濃度最高,坡耕地次之,林地和水庫的氮磷輸出濃度均較低。同時(shí)在長期監(jiān)測中發(fā)現(xiàn),氮磷濃度變化較大,原因可能是受土地利用類型、降雨、施肥和時(shí)空變化等因素的影響。

    表2 不同土地利用類型的各形態(tài)氮磷濃度統(tǒng)計(jì)分析

    注:不同字母代表差異顯著(<0.05)。

    Note: Different letters represents significant difference at<0.05.

    2.2 水質(zhì)的時(shí)間變化特征

    長坪小流域不同土地利用類型氮磷濃度表現(xiàn)出較強(qiáng)的時(shí)間異質(zhì)性(圖2)。對(duì)于不同形態(tài)的氮濃度而言,村莊TN、NO3--N和NH4+-N濃度峰值多出現(xiàn)于非汛期的11月至次年1月,村莊非汛期的TN、NO3--N和NH4+-N濃度分別為10.14、8.23和0.75 mg/L,高于汛期的6.44、5.23和 0.26 mg/L。與村莊不同,坡耕地TN、NO3--N和NH4+-N濃度峰值則出現(xiàn)在每年汛期的5-8月。林地和水庫不同形態(tài)氮濃度沒有表現(xiàn)出明顯的季節(jié)性差異。對(duì)于不同形態(tài)的磷濃度而言,TP與TN的季節(jié)變化規(guī)律相似,村莊非汛期TP濃度為0.31 mg/L,高于汛期的0.19 mg/L,坡耕地汛期TP濃度為0.15 mg/L,高于非汛期的0.11 mg/L。林地和水庫TP濃度較低且變化范圍較小。DTP和PP的變化規(guī)律與TP基本一致,村莊DTP和PP濃度峰值多出現(xiàn)在11月至次年2月,與之不同,坡耕地、林地和水庫DTP和PP濃度峰值則多出現(xiàn)在汛期的6-7月。從年際尺度看,村莊氮磷濃度呈波動(dòng)上升趨勢,坡耕地、林地和水庫沒有明顯變化。

    注:DTP為可溶態(tài)總磷;PP為顆粒態(tài)總磷。下同。

    對(duì)不同形態(tài)氮磷濃度和降雨量進(jìn)行Pearson相關(guān)性分析(表3)發(fā)現(xiàn),村莊的NH4+-N、TP、PP濃度和月降雨量呈極顯著負(fù)相關(guān)(<0.01),DTP濃度和月降雨量呈顯著負(fù)相關(guān)(<0.05)。與村莊不同,坡耕地的TN、NO3--N、TP、DTP濃度和月降雨量呈極顯著正相關(guān)(<0.01),林地的TP、PP濃度和月降雨量呈極顯著正相關(guān)(<0.01),水庫的TN、NO3--N、TP、顆粒態(tài)總磷(PP)濃度和月降雨量呈極顯著正相關(guān)(<0.01)。從整體上看,隨著降雨量的增加,村莊氮磷濃度逐漸降低,坡耕地、林地和水庫的氮磷濃度逐漸增加,且坡耕地氮磷濃度與降雨量的正相關(guān)性高于林地和水庫。

    表3 不同形態(tài)氮磷濃度與降雨量的相關(guān)性分析

    注:*和**分別表示在<0.05水平相關(guān)和<0.01水平相關(guān)。

    Note:*and**indicate correlation at the<0.05 level and the<0.01 level, respectively.

    2.3 不同土地利用類型水質(zhì)的綜合評(píng)價(jià)

    基于長坪小流域2015-2020年的水質(zhì)監(jiān)測數(shù)據(jù),采用單因子污染指數(shù)法和內(nèi)梅羅污染指數(shù)法相結(jié)合的方法,執(zhí)行GB 3838-2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》Ⅴ類標(biāo)準(zhǔn)來綜合評(píng)價(jià)不同土地利用類型水體水質(zhì),水質(zhì)評(píng)價(jià)結(jié)果見表4。

    表4 不同土地利用類型年際水質(zhì)評(píng)價(jià)結(jié)果

    注:為內(nèi)梅羅指數(shù)。

    Note:is Nemerow index.

    從TN污染指數(shù)看,2015-2020年間,村莊、坡耕地、林地和水庫的TN污染指數(shù)變化范圍分別是2.88~6.41、1.42~1.87、0.67~0.94、0.50~0.82,村莊和坡耕地TN污染指數(shù)均大于1,林地和水庫的TN污染指數(shù)小于1,村莊和坡耕地的TN超標(biāo)。NH4+-N污染指數(shù)均小于1,TP污染指數(shù)總體小于1,長坪小流域水質(zhì)超標(biāo)項(xiàng)目主要為TN。單因子評(píng)價(jià)結(jié)果表明2015-2020年村莊和坡耕地水質(zhì)均為劣Ⅴ類水。內(nèi)梅羅綜合污染指數(shù)評(píng)價(jià)結(jié)果顯示,林地和水庫的污染程度均為清潔水平,坡耕地污染程度為輕污染水平,村莊的污染程度從污染水平轉(zhuǎn)變?yōu)橹匚廴舅?。水質(zhì)的綜合評(píng)價(jià)結(jié)果表明,水質(zhì)污染程度從重到輕為村莊、坡耕地、林地、水庫。

    從表4可知,村莊內(nèi)梅羅污染指數(shù)呈增加趨勢,內(nèi)梅羅污染指數(shù)從2015年的2.20增加到2020年的4.88,且村莊的內(nèi)梅羅污染指數(shù)顯著高于其他土地利用類型,村莊水質(zhì)的惡化使小流域整體水質(zhì)的污染程度不斷增加。坡耕地內(nèi)梅羅污染指數(shù)呈緩慢增加趨勢,這可能與化肥的使用量增加有關(guān)。林地和水庫的內(nèi)梅羅污染指數(shù)變化穩(wěn)定,水質(zhì)未受到明顯污染。

    2.4 不同污染源的氮磷輸出負(fù)荷分析

    從表5可以看出,長坪小流域TN和TP輸出負(fù)荷總量分別為4 278.59和364.93 kg/a。畜禽養(yǎng)殖源是長坪小流域最主要的污染源,其TN輸出負(fù)荷為1 954.84 kg/a,TP輸出負(fù)荷為261.92 kg/a,分別占流域TN和TP輸出負(fù)荷總量的45.69%和71.77%,其中生豬養(yǎng)殖氮磷輸出負(fù)荷最高,占畜禽養(yǎng)殖源氮磷輸出負(fù)荷的95%以上。因此,有效控制畜禽養(yǎng)殖源,是提高小流域水質(zhì)的關(guān)鍵所在。農(nóng)村生活源的TN和TP的輸出負(fù)荷分別占流域TN和TP輸出負(fù)荷總量的15.30%和9.18%。長坪小流域主要土地利用類型為林地和坡耕地,分別占流域面積的64.11%和28.73%,因此林地和坡耕地的氮磷輸出負(fù)荷也較高,坡耕地和林地的TN輸出負(fù)荷總量分別為1 098.52和570.51 kg/a,分別占流域TN輸出負(fù)荷總量的25.67%和13.34%,坡耕地和林地的TP輸出負(fù)荷總量分別為29.35和40.15 kg/a,分別占流域TP輸出負(fù)荷總量的8.04%和11.01%。畜禽養(yǎng)殖源和種植業(yè)源是小流域的主要污染源,是影響小流域水質(zhì)的關(guān)鍵所在。

    表5 不同污染源氮磷污染負(fù)荷及貢獻(xiàn)率

    2.5 分區(qū)協(xié)同防控

    以三峽庫區(qū)典型流域長坪小流域?yàn)槔?,首先按照不同土地利用類型水質(zhì)污染程度分為清潔、輕污染和重污染;再參照《全國生態(tài)功能區(qū)劃(修編版)》中生態(tài)系統(tǒng)服務(wù)功能類型,將流域內(nèi)不同土地利用類型的生態(tài)系統(tǒng)服務(wù)功能分為水源涵養(yǎng)、人居保障和水土保持;綜合土地利用類型、污染程度和生態(tài)系統(tǒng)服務(wù)功能,最終將流域防控分區(qū)分為林草水源涵養(yǎng)區(qū)、村莊污染控制區(qū)和坡耕地水土保持區(qū)(表6)。

    表6 防控分區(qū)劃分

    結(jié)合流域內(nèi)污染特征和不同分區(qū)的防控目標(biāo),制定分區(qū)協(xié)同防控策略,該策略包括3個(gè)方面內(nèi)容:

    1)建立防控分區(qū)。將流域分為林草水源涵養(yǎng)區(qū)、村莊污染控制區(qū)和坡耕地水土保持區(qū)。林草水源涵養(yǎng)區(qū):林地是長坪小流域最主要的土地利用類型,林地面積占流域總面積的64.11%,林地出水水質(zhì)較好(TN 1.57 mg/L,表2)且常年穩(wěn)定,污染程度為清潔水平(表4),該分區(qū)生態(tài)系統(tǒng)服務(wù)功能為水源涵養(yǎng),對(duì)該區(qū)域產(chǎn)生的清水應(yīng)采取保護(hù)策略,避免清水過村后遭受污染。村莊污染控制區(qū):該分區(qū)面積17.04 hm2,占流域總面積的5.17%,TN和TP的污染負(fù)荷卻分別占流域總污染負(fù)荷的60.99%和80.95%,水質(zhì)較差(污染程度為重污染,表4),同時(shí)該分區(qū)存在著人畜混居普遍、畜禽糞便露天堆置遇雨流失問題突出的現(xiàn)象,該分區(qū)生態(tài)系統(tǒng)服務(wù)功能為人居保障,應(yīng)加強(qiáng)污水和糞污的收集、處理和利用,從源頭上控制氮磷污染的輸出。坡耕地水土保持區(qū):該分區(qū)占流域總面積的28.73%,TN和TP的污染負(fù)荷分別占流域總污染負(fù)荷的25.67%和8.04%(表5),坡耕地周邊水體TN和TP濃度分別為2.88和0.13 mg/L(表2),污染程度屬于輕污染,該分區(qū)生態(tài)系統(tǒng)服務(wù)功能為水土保持,應(yīng)以控制水土流失和氮磷流失為主要防控策略,強(qiáng)化徑流的集蓄和梯級(jí)再利用。

    2)分區(qū)協(xié)同防控。通過促進(jìn)各分區(qū)間的水和污染物的流動(dòng),建立流域農(nóng)業(yè)面源污染物分區(qū)協(xié)同防控機(jī)制。林地作為水源涵養(yǎng)區(qū),一方面可為村莊提供生活用水,另一方面林地所產(chǎn)清水流經(jīng)村莊與污水混合后會(huì)帶來污水量大、處理難等問題,應(yīng)在清水入村前進(jìn)行分流調(diào)控,一部分供應(yīng)村莊用水,其余通過清水通道直接入湖/河;坡耕地單位面積TN和TP污染負(fù)荷僅為村莊的7.58%和1.79%,可以通過坡耕地就近就地消納村莊糞便污水,協(xié)同實(shí)現(xiàn)畜禽養(yǎng)殖糞便污染治理與資源化利用,既可以節(jié)約水肥,還可以降低糞污處理成本。

    3)重點(diǎn)防控關(guān)鍵源和關(guān)鍵期。由于不同土地利用類型氮磷濃度在時(shí)間上存在異質(zhì)性(圖2),在汛期,應(yīng)重點(diǎn)管控坡耕地水土與氮磷流失,不僅要防止氮磷肥的過量施用,還要通過集水池或生態(tài)塘等對(duì)汛期的徑流進(jìn)行攔蓄再利用。在非汛期,應(yīng)重點(diǎn)管控村莊中畜禽養(yǎng)殖糞污污染(畜禽養(yǎng)殖源TN和TP污染負(fù)荷比例分別占流域的45.69%和71.77%),糞污貯存和處理設(shè)施的設(shè)計(jì)容積應(yīng)充分考慮糞污產(chǎn)生與利用的錯(cuò)期問題。

    3 討 論

    3.1 時(shí)間變化對(duì)水質(zhì)的影響

    研究發(fā)現(xiàn),2015—2020年研究區(qū)村莊水質(zhì)呈不斷惡化趨勢,TN和TP平均濃度從2015年的6.01和0.22 mg/L增加到2020年的12.69和0.35 mg/L,內(nèi)梅羅污染指數(shù)評(píng)價(jià)結(jié)果表明2020年村莊污染程度為重污染,分析其原因主要是由于受近年來市場行情影響,長坪小流域生豬養(yǎng)殖量增加50%左右。在本研究中,不同土地利用類型氮磷濃度隨季節(jié)變化差異明顯,村莊非汛期氮磷濃度高于汛期,該時(shí)期TN和TP濃度分別比汛期高57.45%和63.16%,坡耕地非汛期氮磷濃度均低于汛期,不同土地利用類型氮磷濃度的季節(jié)性差異主要是受降雨因素的影響,降雨具有污染和稀釋的雙重作用[26],表3的分析結(jié)果證實(shí),降雨量和村莊氮磷濃度負(fù)相關(guān),和坡耕地、林地和水庫的氮磷濃度正相關(guān)。Xu等[27]的研究也表明,村莊排放污水氮磷濃度在汛期因雨水稀釋濃度降低,農(nóng)田汛期氮磷濃度高于非汛期。

    3.2 土地利用對(duì)水質(zhì)的影響

    土地利用類型是影響流域徑流產(chǎn)生、土壤侵蝕和氮磷輸出負(fù)荷的重要因素[28],研究發(fā)現(xiàn),流域內(nèi)氮磷輸出濃度從高到低依次為村莊、坡耕地、林地、水庫,這與已有的研究結(jié)果相似[4]。本研究結(jié)果表明村莊和坡耕地對(duì)流域水質(zhì)起負(fù)向作用,林地和水庫對(duì)流域水質(zhì)起正向作用,這與范雅雙等[29]的研究結(jié)果相同。村莊靠近河道,且村莊的不透水表面的面積較高,道路、屋頂?shù)仁刮廴疚锔菀讌R聚,生活污水和畜禽養(yǎng)殖糞便等直接排放致使大量氮磷污染物進(jìn)入河道,對(duì)河道造成負(fù)擔(dān)[30]。坡耕地受水土流失以及農(nóng)業(yè)活動(dòng)等因素的影響,向河道輸入大量氮磷營養(yǎng)物,同時(shí)也作為植被或濕地系統(tǒng)吸附,消納污染物,對(duì)進(jìn)入河道的污染物有截留作用[31]。林地植被覆蓋度高,根系發(fā)達(dá),可以有效攔截污染物,減少氮磷的輸出。水庫具有較強(qiáng)的自我凈化能力,受人類活動(dòng)的影響較小,氮磷濃度保持低水平,水庫出水能降低流域河道氮磷濃度。

    3.3 流域氮磷輸出負(fù)荷的影響因素

    在本研究區(qū)中,畜禽養(yǎng)殖源是長坪小流域最主要的氮磷污染源,其TN和TP的輸出負(fù)荷分別占流域TN和TP輸出負(fù)荷的45.69%和71.77%。張維理等[32]的研究發(fā)現(xiàn)即使少量的畜禽養(yǎng)殖廢水進(jìn)入水體,也容易使水體富營養(yǎng)化。李榮剛等[33]對(duì)太湖地區(qū)水體污染物排放量的研究表明,畜禽養(yǎng)殖排放的TP占水體總量的86.4%。連慧姝等[34]的研究也表明,畜禽養(yǎng)殖源產(chǎn)生的TP占全流域TP產(chǎn)生量的77.24%。在本研究中,農(nóng)村生活源TN和TP的輸出量雖然僅占全流域的15.30%和9.18%,但由于村莊距河道位置較近,污水在河道匯合處會(huì)引起氮磷濃度的升高,使得河道的局部污染嚴(yán)重。種植業(yè)源也是小流域重要的污染源,TN輸出負(fù)荷占流域TN輸出負(fù)荷總量的25.67%,長坪小流域坡耕地基本采用順坡種植,在6-7月氮磷流失風(fēng)險(xiǎn)期,作物處于苗期,地表裸露大,水土流失風(fēng)險(xiǎn)大。因此,減少村莊和坡耕地的氮磷污染負(fù)荷輸出,是長坪小流域農(nóng)業(yè)面源污染防控的關(guān)鍵。

    3.4 流域農(nóng)業(yè)面源污染分區(qū)協(xié)同防控策略

    《關(guān)于加快推進(jìn)長江經(jīng)濟(jì)帶農(nóng)業(yè)面源污染治理的指導(dǎo)意見》中指出要統(tǒng)籌推進(jìn)種植業(yè)、養(yǎng)殖業(yè)和農(nóng)村生活污染防治,系統(tǒng)設(shè)計(jì)治理方案。因此,需要改變按單個(gè)污染源治理的傳統(tǒng)模式,本文從流域面源污染治理角度,基于流域內(nèi)土地利用類型、污染特征及功能進(jìn)行防控分區(qū),將流域分為林草水源涵養(yǎng)區(qū)、村莊污染控制區(qū)和坡耕地水土保持區(qū),依據(jù)不同分區(qū)的保護(hù)與治理需求,實(shí)施差別化的管理與協(xié)同防控策略。以流域?yàn)閱卧?,統(tǒng)籌考慮種植業(yè)、畜禽養(yǎng)殖業(yè)和農(nóng)村生活污染排放特征,明確污染治理重點(diǎn)和關(guān)鍵期,注重各分區(qū)間水和污染物的緊密銜接,以最低成本實(shí)現(xiàn)最佳治理效果。

    氮磷對(duì)于水體是污染物,但是對(duì)于作物是生長發(fā)育必需的營養(yǎng)物質(zhì)。對(duì)于農(nóng)業(yè)生態(tài)系統(tǒng)來說,最好的農(nóng)業(yè)面源污染治理方式是將各個(gè)污染源中的氮磷通過土地消納,變廢為寶,實(shí)現(xiàn)雙贏。分散養(yǎng)殖的畜禽糞污經(jīng)過無害化處理后就近就地施用于庭院周邊的坡耕地、果園、菜地中比較經(jīng)濟(jì)可行,既促進(jìn)糞肥利用,又可減少污染排放,實(shí)現(xiàn)生產(chǎn)、生活清潔化。在三峽庫區(qū)小流域,由于溝道短,水流急,導(dǎo)致流域自凈能力差,除了考慮土地利用類型外,還應(yīng)根據(jù)流域?qū)嶋H情況,在坡耕地與消落帶或水體交匯處建立臨水生態(tài)凈化區(qū),通過生態(tài)溝、塘、植被過濾帶或消落帶等凈化后再排入水體,對(duì)坡耕地徑流中的泥沙、氮磷和農(nóng)藥等污染物進(jìn)行有效截留[35-36],提升流域自凈能力。

    4 結(jié) 論

    1)不同土地利用類型氮磷輸出濃度具有較強(qiáng)空間異質(zhì)性,氮磷輸出濃度從高到低依次為村莊、坡耕地、林地、水庫;流域氮磷輸出濃度的季節(jié)變化存在顯著差異,村莊氮磷輸出濃度為非汛期高于汛期,坡耕地、林地和水庫則表現(xiàn)為汛期高于非汛期;畜禽養(yǎng)殖源是該流域最主要的氮磷污染源,TN和TP的輸出負(fù)荷占比分別為流域總負(fù)荷的45.69%和71.77%,對(duì)畜禽養(yǎng)殖污染進(jìn)行合理控制,可有效提高小流域水質(zhì);坡耕地與村莊相比具有顯著的低污染特征,小流域坡耕地的TN和TP單位面積污染負(fù)荷分別是村莊的7.58%和1.79%,因此,可以以坡耕地消納村莊污水和畜禽糞便。

    2)根據(jù)流域內(nèi)不同土地利用類型的水質(zhì)時(shí)空變化特征,可將流域分為林草水源涵養(yǎng)區(qū)、村莊污染控制區(qū)和坡耕地水土保持區(qū),建立分區(qū)協(xié)同防控策略,促進(jìn)各分區(qū)間的水和污染物的流動(dòng),對(duì)關(guān)鍵污染源和關(guān)鍵時(shí)期進(jìn)行重點(diǎn)防控,進(jìn)而實(shí)現(xiàn)對(duì)流域農(nóng)業(yè)面源污染的低成本和高效治理。

    [1] Chen C, Gao M, Xie D, et al. Spatial and temporal variations in non-point source losses of nitrogen and phosphorus in a small agricultural catchment in the Three Gorges Region[J]. Environmental Monitoring & Assessment, 2016, 188(4): 257.

    [2] Francisco J D, Anthony T, Dahlgren R A. Agricultural pollutant removal by constructed wetlands: Implications for water management and design[J]. Agricultural Water Management, 2012, 104: 171-183.

    [3] Wu J, Zeng S, Yang L, et al. Spatiotemporal characteristics of the water quality and its multiscale relationship with land use in the Yangtze river basin[J]. Remote Sensing, 2021, 13(16): 3309.

    [4] 閔惠琳,宋維峰. 全福莊小流域不同土地利用類型地表水氮濃度動(dòng)態(tài)特征[J]. 水土保持研究,2019,26(1):92-97.

    Min Huilin, Song Weifeng. Dynamic characterization of surface water nitrogen concentration under various land uses in Quanfuzhuang watershed[J]. Research of Soil and Water Conservation, 2019, 26(1): 92-97. (in Chinese with English abstract)

    [5] 王曉燕,王一峋,王曉峰,等. 密云水庫小流域土地利用方式與氮磷流失規(guī)律[J]. 環(huán)境科學(xué)研究,2003,16(1):30-33.

    Wang Xiaoyan, Wang Yixun, Wang Xiaofeng, et al. The character of nutrient loss and land use in a small watershed of Miyun reservoir[J]. Research of Environmental Sciences, 2003, 16(1): 30-33. (in Chinese with English abstract)

    [6] 紀(jì)仁婧,洪大林,和玉璞,等. 南方低山丘陵區(qū)小流域不同土地利用方式下面源污染分布特征[J]. 水資源與水工程學(xué)報(bào),2020,31(4):181-185.

    Ji Renjing, Hong Dalin, He Yupu, et al. Non-point source pollution distribution characteristics under different landuses in small watershed in the hilly regions of south China[J]. Journal of Water resources & Water Engineering, 2020, 31(4): 181-185. (in Chinese with English abstract)

    [7] 徐啟渝,王鵬,王濤,等. 土地利用結(jié)構(gòu)與景觀格局對(duì)鄱陽湖流域贛江水質(zhì)的影響[J]. 湖泊科學(xué),2020,32(4):1008-1019.

    Xu Qiyu, Wang Peng, Wang Tao, et al. Investigation of the impacts of land use structure and landscape pattern on water quality in the Ganjiang river, Lake Poyang Basin[J]. Journal of Lake Sciences, 2020, 32(4): 1008-1019. (in Chinese with English abstract)

    [8] 彭勃,付永勝,趙君鳳,等. 小流域土地利用景觀格局對(duì)水質(zhì)的影響[J]. 生態(tài)科學(xué),2019,38(3):90-99.

    Peng Bo, Fu Yongsheng, Zhao Junfeng, et al. Effect of land use landscape pattern on water quality in small watershed[J]. Ecological Science, 2019, 38(3): 90-99. (in Chinese with English abstract)

    [9] 吳東,黃志霖,肖文發(fā),等. 三峽庫區(qū)小流域土地利用結(jié)構(gòu)變化及其氮素輸出控制效應(yīng):以蘭陵溪小流域?yàn)槔齕J]. 環(huán)境科學(xué),2016,37(8):2940-2946.

    Wu Dong, Huang Zhilin, Xiao Wenfa, et al. Land use structure change and its control effect of nitrogen output in a small watershed of Three Gorges Reservoir Area: A case study of Lanlingxi watershed[J]. Environmental Science, 2016, 37(8): 2940-2946. (in Chinese with English abstract)

    [10] 華玲玲,李文超,翟麗梅,等. 三峽庫區(qū)古夫河小流域氮磷排放特征[J]. 環(huán)境科學(xué),2017,38(1):138-146.

    Hua Lingling, Li Wenchao, Zhai Limei, et al. Characteristics of nitrogen and phosphorus emissions in the Gufu river small watershed of the Three Georges Reservoir Area[J]. Environmental Science, 2017, 38(1): 138-146. (in Chinese with English abstract)

    [11] 梁爽,陳敏,肖尚斌,等.鄂西長江喀斯特小流域氮磷輸出特征[J]. 長江流域資源與環(huán)境,2021,30(10):2471-2481.

    Liang Shuang, Chen Min, Xiao Shangbin, et al. Nitrogen and phosphorus output in a small karst basin of the Yangtze River in western Hubei[J]. Resources and Environment in the Yangtze Basin, 2021, 30(10): 2471-2481. (in Chinese with English abstract)

    [12] Wang T, Xiao W, Huang Z, et al. Interflow pattern govern nitrogen loss from tea orchard slopes in response to rainfall pattern in Three Gorges Reservoir Area[J]. Agricultural Water Management, 2022, 269: 107684.

    [13] Zhang T, Ni J, Xie D. Severe situation of rural nonpoint source pollution and efficient utilization of agricultural wastes in the Three Gorges Reservoir Area[J]. Environmental Science and Pollution Research, 2015, 22(21): 16453.

    [14] 郭萍,夏振堯,高峰,等. 香根草植物籬對(duì)三峽庫區(qū)坡地紫色土侵蝕的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2021,37(19):105-112.

    Guo Ping, Xia Zhenyao, Gao Feng, et al. Effects ofhedgerow on the erosion of purple soil of slope land in the Three Gorges Reservoir Area of China[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2021, 37(19): 105-112. (in Chinese with English abstract)

    [15] Zhang Y, Xie D, Ni J, et al. Conservation tillage practices reduce nitrogen losses in the sloping upland of the Three Gorges Reservoir Area: No-till is better than mulch-till[J]. Agriculture, Ecosystems & Environment, 2020, 300: 107003.

    [16] 楊桂山. 流域綜合管理導(dǎo)論[M]. 北京:科學(xué)出版社,2004.

    [17] 王文靜,高鵬程,李捷,等. 丹江口水庫典型入庫支流水質(zhì)評(píng)價(jià)與趨勢分析[J]. 水資源保護(hù),2016,32(3):127-130.

    Wang Wenjing, Gao Pengcheng, Li Jie, et al. Water quality assessment and trend analysis of typical tributaries of Danjiangkou reservoir[J]. Water Resources Protection, 2016, 32(3): 127-130. (in Chinese with English abstract)

    [18] Nemerow N. Scientific Stream Pollution Analysis[M]. New York: McGraw-Hil1, 1974.

    [19] 潘犖,黃曉榮,魏曉玥,等. 三種常用水質(zhì)評(píng)價(jià)方法的對(duì)比分析研究[J]. 中國農(nóng)村水利水電,2019(6):51-55.

    Pan Luo, Huang Xiaorong, Wei Xiaoyue, et al. A comparative analysis of three common water quality evaluation methods[J]. China Rural Water and Hydropower, 2019(6): 51-55. (in Chinese with English abstract)

    [20] Johnes P J, Heathwaite A L. Modelling the impact of land use change on water quality in agricultural catchments[J]. Hydrological Processes, 1997, 11: 269-286.

    [21] 龔世飛,丁武漢,居學(xué)海,等. 典型農(nóng)業(yè)小流域面源污染源解析與控制策略:以丹江口水源涵養(yǎng)區(qū)為例[J]. 中國農(nóng)業(yè)科學(xué),2021,54(18):3919-3931.

    Gong Shifei, Ding Wuhan, Ju Xuehai, et al. Source analysis and control strategies of non-point source pollution in typical agricultural small watershed: A case study of Danjiangkou water conservation area[J]. Scientia Agricultura Sinica, 2021, 54(18): 3919-3931. (in Chinese with English abstract)

    [22] 王萌,劉云,李春蕾,等. 三峽庫區(qū)湖北段非點(diǎn)源污染氮磷排放時(shí)空分布特征[J]. 水土保持通報(bào),2018,38(2):46-52.

    Wang Meng, Liu Yun, Li Chunlei, et al. Spatio-temporal distribution of pollution loads from agricultural non-point sources nitrogen and phosphorus in Three Gorges Reservoir Area of Hubei province[J]. Bulletin of Soil and Water Conservation, 2018, 38(2): 46-52. (in Chinese with English abstract)

    [23] 第一次全國污染源普查資料編纂委員會(huì). 污染源普查產(chǎn)排污系數(shù)手冊[M]. 北京:中國環(huán)境科學(xué)出版社,2011.

    [24] 任天志. 全國農(nóng)田面源污染排放系數(shù)手冊[M]. 北京:中國農(nóng)業(yè)出版社,2015.

    [25] 劉佳昆. 基于改進(jìn)輸出系數(shù)模型的澎溪河流域非點(diǎn)源污染研究[D]. 重慶:西南大學(xué),2019.

    Liu Jiakun. An Improved Export Coefficient Model to Estimate Nitrogen and Phosphorus Loads in Pengxi River Basin, Three Gorges Reservoir[D]. Chongqing: Southwest University, 2019. (in Chinese with English abstract)

    [26] 唐廉,胡曉輝,權(quán)冠中,等. 潭江流域水質(zhì)時(shí)空分布特征及其與土地利用的相關(guān)性分析[J]. 地球與環(huán)境,2018,46(4):364-372.

    Tang Lian, Hu Xiaohui, Quan Guanzhong, et al. Spatial-temporal distribution characteristics of water quality in Tanjiang River Basin and its correlation with land use[J]. Earth and Environment, 2018, 46(4): 364-372. (in Chinese with English abstract)

    [27] Xu J, Jin G, Tang H, et al. Response of water quality to land use and sewage outfalls in different seasons[J]. Science of the Total Environment, 2019, 696: 134014.

    [28] Chen X, Xu G, Zhang W, et al. Spatial variation pattern analysis of hydrologic processes and water quality in Three Gorges Reservoir Area[J]. Water, 2019, 11(12): 2608.

    [29] 范雅雙,于婉晴,張婧,等. 太湖上游水源區(qū)河流水質(zhì)對(duì)景觀格局變化的響應(yīng)關(guān)系:以東苕溪上游為例[J]. 湖泊科學(xué),2021,33(5):1478-1489.

    Fan Yashuang, Yu Wanqing, Zhang Jing, et al. Response of water quality to landscape pattern change in the water source area of upperreaches of Lake Taihu: A case study in the upper reaches of Dongtiaoxi River[J]. Journal of Lake Sciences, 2021, 33(5): 1478-1489. (in Chinese with English abstract)

    [30] Asnake K, Worku H, Argaw M. Assessing the impact of watershed land use on Kebena river water quality in Addis Ababa, Ethiopia[J]. Environmental Systems Research, 2021, 10(1): 1-14.

    [31] Kadir A, Ahmed Z, Uddin M M,et al. Integrated approach to quantify the impact of land use and land cover changes on water quality of Surma river, Sylhet, Bangladesh [J]. Water, 2022, 14(1): 17.

    [32] 張維理,武淑霞,冀宏杰,等. 中國農(nóng)業(yè)面源污染形勢估計(jì)及控制對(duì)策 I.21世紀(jì)初期中國農(nóng)業(yè)面源污染的形勢估計(jì)[J]. 中國農(nóng)業(yè)科學(xué),2004,37(7):1008-1017.

    Zhang Weili, Wu Shuxia, Ji Hongjie, et al. Estimation of agricultural non-point source pollution in China and the alleviating strategies I. estimation of agricultural non-point source pollution in China in early 21 century [J]. Scientia Agricultura Sinica, 2004, 37(7): 1008-1017. (in Chinese with English abstract)

    [33] 李榮剛,夏源陵. 江蘇太湖地區(qū)水污染物及其向水體的排放量[J]. 湖泊科學(xué),2000,12(2):147-153.

    Li Ronggang, Xia Yuanling. Water pollutants and their discharges to water bodies in Taihu lake, Jiangsu province [J]. Journal of Lake Sciences, 2000, 12(2): 147-153. (in Chinese with English abstract)

    [34] 連慧姝,劉宏斌,李旭東,等. 太湖蠡河小流域水質(zhì)的空間變化特征及污染物源解析[J]. 環(huán)境科學(xué),2017,38(9):3657-3665.

    Lian Huishu, Liu Hongbin, Li Xudong, et al. Analysis of spatial variability of water quality and pollution sources in Lihe river watershed, Taihu lake basin[J]. Environmental Science, 2017, 38(9): 3657-3665. (in Chinese with English abstract)

    [35] 張維,孫虹蕾,唐翔宇. 三峽水庫消落帶土壤膠體釋放與遷移特[J]. 農(nóng)業(yè)工程學(xué)報(bào),2021,37(1):188-194.

    Zhang Wei, Sun Honglei, Tang Xiangyu. Soil colloid release and transport in the water level fluctuation zone of the Three Gorges Reservoir[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2021, 37(1): 188-194. (in Chinese with English abstract)

    [36] 李?yuàn)櫇桑嚝h,施鳳寧,等. 水庫消落帶研究進(jìn)展[J]. 濕地科學(xué),2019,17(6):689-696.

    Li Shanze, Deng Yue, Shi Fengning, et al. Research progress on reservoir subsidence zone[J]. Wetland Science, 2019, 17(6): 689-696. (in Chinese with English abstract)

    Spatio-temporal characteristics of water quality and pollution prevention strategies in typical watersheds in the Three Gorges reservoir areas

    Li Ganghao1,2, Fan Xianpeng1, Xia Ying1※, Liu Hongbin3, Wu Maoqian1, Zhang Fulin1, Huang Min2, Zhai Limei3, Zhou Jiwen4, Kong Xiangqiong4, Cheng Zizhen1

    (1.,/,430064,;2.,,430070,; 3.,,100081,; 4.,443700,)

    Water quality has been greatly confined to the nitrogen loss and phosphorus pollutants caused by agricultural production in the Three Gorges Reservoir Area. This study aims to clarify the characteristics of water quality under different land use types in typical small watersheds. The water quality monitoring was conducted in the Changping small watershed in the Three Gorges Reservoir Area from 2015 to 2020. A systematic investigation was made on the change patterns of Total Nitrogen (TN), nitrate nitrogen (NO3--N), ammonia nitrogen (NH4+-N), Total Phosphorus (TP), Dissolved Total Phosphorus (DTP), and Total Particulate Phosphorus (PP) concentrations in the water bodies under land use types, together with their influencing factors. The results showed that: 1) The overall nitrogen and phosphorus output concentrations of land use types were ranked in the order of the village > sloping farmland > forest land > reservoir. The average TN concentrations were 8.29, 2.88, 1.57, and 1.43 mg/L in the villages, sloping farmland, forest land, and reservoirs, respectively. The average TP concentrations were 0.25, 0.13, 0.09, and 0.07 mg/L, respectively. 2) There was a different water quality of land use types during the flood and non-flood periods. The concentration of nitrogen and phosphorus output in the villages was higher in the non-flood period than that in the flood period. The sloping farmland, forest land, and reservoirs showed a higher concentration in the flood period than those in the non-flood period. 3) The Nemerow comprehensive pollution index demonstrated that the pollution level of both forest land and reservoir was at the clean level, and the sloping farmland was at the light pollution level. But, the pollution level of the village deteriorated from the pollution level to a heavy pollution level, due to the increase in pig farming. 4) The annual output loads of TN and TP in the Changping sub-basin were 4 278.59 and 364.93 kg/a , respectively. Livestock and poultry breeding was the most important source of nitrogen and phosphorus pollution in the sub-basin, with the output loads of TN and TP accounting for 45.69% and 71.77%, respectively. The TN and TP pollution load per unit area of sloping land in the sub-basin were 7.58% and 1.79% of those of village, respectively. The sloping farmland presented significantly low pollution characteristics, compared with the village. Therefore, the sloping farmland can be used to absorb the village sewage and livestock manure, in order to promote the local absorption of manure in the basin for the better effect of multi-source co-management. The land use structure of the whole watershed can also be adjusted in an effective way to control agricultural surface pollution. Therefore, the small watershed was divided into the forest and grass water connotation zone, village pollution control zone, as well as sloping farmland soil and water conservation zone, according to the water quality and ecosystem service function of land use types. Zoning collaborative prevention and control strategy should be an effective way of promoting the water resource, nitrogen and phosphorus transfer from one pollution control zone to another pollution control zone in order to preventing the pollutants export to the out of the small watershed. The finding can provide the theoretical basis and technical support for the control of nitrogen and phosphorus pollution in the protection of water sources in the Three Gorges Reservoir Area.

    water quality; land use; spatial distribution; Three Gorges Reservoir Area; typical watershed; agricultural non-point source pollution; zoning collaborative prevention and control strategy

    10.11975/j.issn.1002-6819.2022.17.013

    X832

    A

    1002-6819(2022)-17-0118-09

    李剛浩,范先鵬,夏穎,等. 三峽庫區(qū)典型流域水質(zhì)時(shí)空特征及污染防控策略[J]. 農(nóng)業(yè)工程學(xué)報(bào),2022,38(17):118-126.doi:10.11975/j.issn.1002-6819.2022.17.013 http://www.tcsae.org

    Li Ganghao, Fan Xianpeng, Xia Ying, et al. Spatio-temporal characteristics of water quality and pollution prevention strategies in typical watersheds in the Three Gorges reservoir areas[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(17): 118-126. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.17.013 http://www.tcsae.org

    2022-04-27

    2022-08-16

    國家自然科學(xué)基金聯(lián)合基金項(xiàng)目(U21A2025)

    李剛浩,研究方向?yàn)檗r(nóng)業(yè)面源污染防治。Email:lghdc20@163.com

    夏穎,博士,副研究員,研究方向?yàn)檗r(nóng)業(yè)面源污染防治研究。Email:xiayinghappy105@163.com

    猜你喜歡
    坡耕地氮磷林地
    玻利維亞擬建新的氮磷鉀肥料工廠
    丹東市林地分類研究
    建平縣實(shí)施國家坡耕地治理項(xiàng)目成效及經(jīng)驗(yàn)
    淺談林地保護(hù)及恢復(fù)措施
    資陽市雁江區(qū):防治并重 建管結(jié)合 創(chuàng)建坡耕地水土流失綜合治理示范區(qū)
    中國水利(2015年24期)2015-04-25 08:14:51
    植物需要氮磷鉀
    種植苧麻對(duì)南方坡耕地土壤抗蝕性的影響
    林地流轉(zhuǎn)模式的選擇機(jī)理及其政策啟示
    小型無人飛行器用于林地監(jiān)視的嘗試
    不同氮磷水平下AM真菌對(duì)棉花的侵染率及氮磷含量的影響
    亚洲精品在线美女| 18禁观看日本| 一区福利在线观看| 国产一区二区三区在线臀色熟女 | ponron亚洲| 亚洲第一青青草原| av网站免费在线观看视频| 女人被狂操c到高潮| 日日夜夜操网爽| 国产精品免费视频内射| 国产真人三级小视频在线观看| 国产精品九九99| 两性夫妻黄色片| 国产亚洲欧美98| 男男h啪啪无遮挡| 久久久久久免费高清国产稀缺| av视频免费观看在线观看| 男女下面插进去视频免费观看| 亚洲成国产人片在线观看| 久久久久久久久中文| 女人被躁到高潮嗷嗷叫费观| 免费在线观看完整版高清| 久久久久国内视频| 欧美中文综合在线视频| www日本在线高清视频| www.自偷自拍.com| 亚洲精品美女久久av网站| 一区二区三区精品91| 高清在线国产一区| 曰老女人黄片| 在线观看一区二区三区激情| 99精国产麻豆久久婷婷| 亚洲欧美激情综合另类| 麻豆一二三区av精品| 午夜视频精品福利| 天堂动漫精品| 在线观看免费日韩欧美大片| 国产精品一区二区在线不卡| 黄色成人免费大全| 久久亚洲精品不卡| 日韩av在线大香蕉| 91麻豆av在线| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美黑人欧美精品刺激| 国产成人影院久久av| 日韩欧美在线二视频| 又黄又爽又免费观看的视频| 国产一区在线观看成人免费| 性色av乱码一区二区三区2| 操美女的视频在线观看| 中文字幕av电影在线播放| 午夜视频精品福利| 久久香蕉国产精品| 99精品久久久久人妻精品| 啦啦啦在线免费观看视频4| 在线观看www视频免费| 亚洲成人国产一区在线观看| 久久久久国内视频| 搡老乐熟女国产| 欧美成人午夜精品| 一进一出好大好爽视频| 久久天躁狠狠躁夜夜2o2o| 色尼玛亚洲综合影院| 色老头精品视频在线观看| 丁香欧美五月| 亚洲欧美一区二区三区黑人| 国产单亲对白刺激| 99国产精品一区二区三区| 天天影视国产精品| 狠狠狠狠99中文字幕| 成熟少妇高潮喷水视频| 激情在线观看视频在线高清| 午夜福利欧美成人| 男女下面插进去视频免费观看| 免费久久久久久久精品成人欧美视频| 国产精品国产高清国产av| avwww免费| 老鸭窝网址在线观看| 夫妻午夜视频| 777久久人妻少妇嫩草av网站| 日韩欧美一区二区三区在线观看| 又大又爽又粗| 丰满迷人的少妇在线观看| 妹子高潮喷水视频| 中文字幕人妻丝袜一区二区| 波多野结衣一区麻豆| 脱女人内裤的视频| 午夜精品国产一区二区电影| av有码第一页| 动漫黄色视频在线观看| 1024香蕉在线观看| 两性夫妻黄色片| 日日夜夜操网爽| 成人国产一区最新在线观看| 黄色视频,在线免费观看| 别揉我奶头~嗯~啊~动态视频| 国产99久久九九免费精品| 男人舔女人的私密视频| 国产亚洲欧美在线一区二区| x7x7x7水蜜桃| 91九色精品人成在线观看| av中文乱码字幕在线| 日本免费a在线| 午夜老司机福利片| 精品久久久精品久久久| 无遮挡黄片免费观看| 国产成人精品久久二区二区91| 法律面前人人平等表现在哪些方面| 女同久久另类99精品国产91| 国产精品二区激情视频| tocl精华| 中文字幕av电影在线播放| 婷婷精品国产亚洲av在线| 日韩视频一区二区在线观看| 可以免费在线观看a视频的电影网站| 一边摸一边抽搐一进一出视频| 国产主播在线观看一区二区| 国产黄色免费在线视频| 免费在线观看日本一区| www.自偷自拍.com| 俄罗斯特黄特色一大片| 夜夜夜夜夜久久久久| 国产精品久久久久成人av| 天天影视国产精品| 久久亚洲真实| 国产精品一区二区免费欧美| 女人被躁到高潮嗷嗷叫费观| 天天添夜夜摸| e午夜精品久久久久久久| 久久久久久久久中文| 欧美亚洲日本最大视频资源| 色在线成人网| 午夜亚洲福利在线播放| 大码成人一级视频| 少妇 在线观看| 黄片小视频在线播放| 国产欧美日韩一区二区精品| 国产亚洲精品综合一区在线观看 | 欧美久久黑人一区二区| 国产精品综合久久久久久久免费 | 电影成人av| 最近最新中文字幕大全电影3 | 在线观看舔阴道视频| 国产在线观看jvid| 午夜福利在线观看吧| xxxhd国产人妻xxx| 久久欧美精品欧美久久欧美| 亚洲一码二码三码区别大吗| 搡老乐熟女国产| 三上悠亚av全集在线观看| 俄罗斯特黄特色一大片| √禁漫天堂资源中文www| 黑人巨大精品欧美一区二区mp4| 两个人免费观看高清视频| 久热这里只有精品99| 黄色女人牲交| 久久久久久久久免费视频了| 亚洲中文日韩欧美视频| 精品国产超薄肉色丝袜足j| 啦啦啦免费观看视频1| 亚洲欧美一区二区三区黑人| 国产成人欧美在线观看| 国产激情欧美一区二区| 高清欧美精品videossex| 日本免费一区二区三区高清不卡 | 精品电影一区二区在线| 无遮挡黄片免费观看| 国产一区二区三区综合在线观看| 欧美日韩乱码在线| 制服诱惑二区| 伦理电影免费视频| 岛国在线观看网站| 免费看a级黄色片| 国产精品九九99| 久久午夜亚洲精品久久| 成人18禁在线播放| 国产精品乱码一区二三区的特点 | 久久狼人影院| 国产精品自产拍在线观看55亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产色婷婷电影| 欧美日韩精品网址| 丁香欧美五月| 少妇裸体淫交视频免费看高清 | 桃色一区二区三区在线观看| 国产精品免费视频内射| 亚洲av日韩精品久久久久久密| 校园春色视频在线观看| 12—13女人毛片做爰片一| 国产伦人伦偷精品视频| 91大片在线观看| 波多野结衣高清无吗| 精品福利永久在线观看| 日日夜夜操网爽| 日韩 欧美 亚洲 中文字幕| 91国产中文字幕| 好看av亚洲va欧美ⅴa在| 欧美黑人精品巨大| 国产成人一区二区三区免费视频网站| 曰老女人黄片| 五月开心婷婷网| 成熟少妇高潮喷水视频| 亚洲国产欧美日韩在线播放| а√天堂www在线а√下载| 50天的宝宝边吃奶边哭怎么回事| 无遮挡黄片免费观看| 久久九九热精品免费| 欧美人与性动交α欧美软件| 国产成年人精品一区二区 | 午夜免费观看网址| 操美女的视频在线观看| 亚洲三区欧美一区| 黄色视频,在线免费观看| 最新美女视频免费是黄的| 涩涩av久久男人的天堂| 超色免费av| 国产亚洲精品一区二区www| 精品久久久久久,| 母亲3免费完整高清在线观看| 窝窝影院91人妻| 久久狼人影院| bbb黄色大片| 亚洲欧美精品综合一区二区三区| 成人亚洲精品一区在线观看| 色尼玛亚洲综合影院| 99国产极品粉嫩在线观看| 丝袜人妻中文字幕| 搡老熟女国产l中国老女人| 国产91精品成人一区二区三区| x7x7x7水蜜桃| 丰满饥渴人妻一区二区三| 男女之事视频高清在线观看| 搡老乐熟女国产| 天天躁狠狠躁夜夜躁狠狠躁| 99在线视频只有这里精品首页| 一边摸一边抽搐一进一出视频| 女性生殖器流出的白浆| 亚洲欧洲精品一区二区精品久久久| 亚洲少妇的诱惑av| 中文欧美无线码| 欧美激情久久久久久爽电影 | 很黄的视频免费| 久久香蕉精品热| 男女之事视频高清在线观看| 日日爽夜夜爽网站| 久久精品国产亚洲av香蕉五月| 日本wwww免费看| 欧美成人午夜精品| 狠狠狠狠99中文字幕| 国产成人免费无遮挡视频| 在线永久观看黄色视频| 欧美中文日本在线观看视频| 欧美亚洲日本最大视频资源| 一级片免费观看大全| 99久久99久久久精品蜜桃| 精品乱码久久久久久99久播| 老司机福利观看| 纯流量卡能插随身wifi吗| 国产精品久久视频播放| 中文字幕精品免费在线观看视频| 久久久水蜜桃国产精品网| 国产成人啪精品午夜网站| av中文乱码字幕在线| 后天国语完整版免费观看| a级毛片黄视频| 亚洲性夜色夜夜综合| 人人妻人人添人人爽欧美一区卜| 精品国产一区二区三区四区第35| 视频区图区小说| 日本黄色视频三级网站网址| 久久久久久久精品吃奶| 韩国av一区二区三区四区| 久久精品影院6| 国产精品永久免费网站| 亚洲美女黄片视频| 中文字幕人妻熟女乱码| 黄色成人免费大全| а√天堂www在线а√下载| 热re99久久国产66热| 欧美激情 高清一区二区三区| 国产亚洲欧美精品永久| 涩涩av久久男人的天堂| 日韩有码中文字幕| 亚洲中文字幕日韩| av超薄肉色丝袜交足视频| 国产熟女xx| 国产三级黄色录像| 国产aⅴ精品一区二区三区波| 亚洲国产看品久久| 老司机午夜十八禁免费视频| 999精品在线视频| 最近最新免费中文字幕在线| 中文字幕精品免费在线观看视频| 国产免费现黄频在线看| 免费在线观看日本一区| a级片在线免费高清观看视频| 亚洲第一欧美日韩一区二区三区| 波多野结衣高清无吗| 国产熟女午夜一区二区三区| 亚洲人成网站在线播放欧美日韩| 亚洲黑人精品在线| 欧美日韩国产mv在线观看视频| 十八禁网站免费在线| 日本三级黄在线观看| 久久人人精品亚洲av| 欧美日韩福利视频一区二区| 热re99久久精品国产66热6| av网站免费在线观看视频| 欧美乱色亚洲激情| 精品久久久久久久毛片微露脸| 看免费av毛片| 国产精品二区激情视频| 久久影院123| 日韩精品免费视频一区二区三区| 久久精品国产99精品国产亚洲性色 | 欧美日韩精品网址| 亚洲av第一区精品v没综合| 亚洲国产欧美网| 欧美激情极品国产一区二区三区| 精品国产乱子伦一区二区三区| 亚洲成人国产一区在线观看| 黄色丝袜av网址大全| 男人舔女人下体高潮全视频| 久久久国产成人精品二区 | 欧洲精品卡2卡3卡4卡5卡区| 日韩 欧美 亚洲 中文字幕| 国产成人欧美在线观看| 亚洲美女黄片视频| 两人在一起打扑克的视频| 伊人久久大香线蕉亚洲五| 欧美日韩一级在线毛片| 亚洲精品国产区一区二| 成人手机av| 在线观看66精品国产| 69av精品久久久久久| 欧美 亚洲 国产 日韩一| 日韩视频一区二区在线观看| 在线国产一区二区在线| 日韩三级视频一区二区三区| 精品久久久久久久久久免费视频 | 成人亚洲精品av一区二区 | 我的亚洲天堂| 亚洲av美国av| 国产精品一区二区免费欧美| 精品乱码久久久久久99久播| 最新美女视频免费是黄的| 一进一出好大好爽视频| 亚洲欧洲精品一区二区精品久久久| 婷婷六月久久综合丁香| 欧美一级毛片孕妇| 国产一区在线观看成人免费| 老司机福利观看| 国产一区在线观看成人免费| 人妻久久中文字幕网| 在线观看一区二区三区激情| 狂野欧美激情性xxxx| 波多野结衣av一区二区av| e午夜精品久久久久久久| 国产男靠女视频免费网站| 极品人妻少妇av视频| 成人国产一区最新在线观看| 国产成人一区二区三区免费视频网站| 欧美+亚洲+日韩+国产| 人人妻人人澡人人看| 亚洲专区中文字幕在线| 99热国产这里只有精品6| 999久久久国产精品视频| 麻豆国产av国片精品| 悠悠久久av| 久久香蕉激情| 日本免费一区二区三区高清不卡 | 99久久99久久久精品蜜桃| 午夜福利在线免费观看网站| 真人做人爱边吃奶动态| 淫秽高清视频在线观看| 午夜日韩欧美国产| 极品教师在线免费播放| 亚洲专区字幕在线| 免费女性裸体啪啪无遮挡网站| 老司机深夜福利视频在线观看| 香蕉丝袜av| 一本综合久久免费| 免费观看人在逋| 老司机深夜福利视频在线观看| 18禁美女被吸乳视频| 久久久久国产精品人妻aⅴ院| 男女下面插进去视频免费观看| 老熟妇乱子伦视频在线观看| 久久久久久久久久久久大奶| 侵犯人妻中文字幕一二三四区| 久久精品91无色码中文字幕| 久久久久久久久中文| 淫秽高清视频在线观看| 国产亚洲精品一区二区www| 免费观看人在逋| 欧美精品啪啪一区二区三区| 久久久久久久久久久久大奶| 黑人巨大精品欧美一区二区mp4| 嫁个100分男人电影在线观看| 日韩中文字幕欧美一区二区| 波多野结衣一区麻豆| 午夜精品国产一区二区电影| 极品人妻少妇av视频| 国产三级在线视频| 欧美成人免费av一区二区三区| 啦啦啦免费观看视频1| 国产精品久久电影中文字幕| av在线播放免费不卡| 成人18禁在线播放| 亚洲男人天堂网一区| 日韩大尺度精品在线看网址 | 99久久国产精品久久久| 午夜免费成人在线视频| 久久精品国产99精品国产亚洲性色 | 亚洲 国产 在线| 99久久精品国产亚洲精品| 国产精品久久久av美女十八| 熟女少妇亚洲综合色aaa.| 日韩欧美一区二区三区在线观看| 久久影院123| 亚洲精品美女久久av网站| 国产成人精品久久二区二区免费| 99国产综合亚洲精品| 欧美精品啪啪一区二区三区| 久久精品成人免费网站| 亚洲人成网站在线播放欧美日韩| 欧美成狂野欧美在线观看| 一进一出抽搐gif免费好疼 | 欧美午夜高清在线| 精品久久久久久久久久免费视频 | 亚洲在线自拍视频| 久久精品人人爽人人爽视色| 超碰成人久久| 在线观看午夜福利视频| 亚洲成人久久性| 国产麻豆69| 亚洲自拍偷在线| 亚洲男人的天堂狠狠| 宅男免费午夜| 50天的宝宝边吃奶边哭怎么回事| 露出奶头的视频| 亚洲一码二码三码区别大吗| 岛国在线观看网站| 久久精品亚洲av国产电影网| 国产深夜福利视频在线观看| 欧美丝袜亚洲另类 | 电影成人av| 亚洲精品在线美女| 99久久精品国产亚洲精品| 国产不卡一卡二| 伦理电影免费视频| 国产三级黄色录像| 欧美日韩瑟瑟在线播放| a在线观看视频网站| ponron亚洲| 欧美丝袜亚洲另类 | 波多野结衣一区麻豆| 亚洲精品在线观看二区| 一级片'在线观看视频| 欧美老熟妇乱子伦牲交| 中文字幕av电影在线播放| 91av网站免费观看| 欧美成人午夜精品| 国产麻豆69| 在线观看一区二区三区| 我的亚洲天堂| 精品国产乱码久久久久久男人| 曰老女人黄片| 69精品国产乱码久久久| 久久人人精品亚洲av| 久久伊人香网站| 日本精品一区二区三区蜜桃| 咕卡用的链子| 亚洲精华国产精华精| 丝袜在线中文字幕| 成年女人毛片免费观看观看9| 国产精品久久久人人做人人爽| av国产精品久久久久影院| 国产免费现黄频在线看| 丝袜美足系列| 五月开心婷婷网| 免费一级毛片在线播放高清视频 | 一夜夜www| 两性午夜刺激爽爽歪歪视频在线观看 | 俄罗斯特黄特色一大片| 男男h啪啪无遮挡| 午夜福利免费观看在线| 日韩大尺度精品在线看网址 | 琪琪午夜伦伦电影理论片6080| 91老司机精品| 日韩欧美一区二区三区在线观看| 精品日产1卡2卡| av天堂在线播放| 黄片大片在线免费观看| a级毛片黄视频| 中文字幕另类日韩欧美亚洲嫩草| 色综合婷婷激情| 天堂√8在线中文| 级片在线观看| 一a级毛片在线观看| 国产麻豆69| 91字幕亚洲| 69精品国产乱码久久久| 久热爱精品视频在线9| 免费搜索国产男女视频| 国产精品乱码一区二三区的特点 | 亚洲五月天丁香| 精品福利观看| av网站免费在线观看视频| 老司机在亚洲福利影院| 久久久久久久精品吃奶| 午夜激情av网站| 一区在线观看完整版| 天天影视国产精品| 免费观看精品视频网站| 天堂中文最新版在线下载| 亚洲成人免费电影在线观看| 国产97色在线日韩免费| 午夜免费观看网址| 搡老熟女国产l中国老女人| 久久久久国内视频| 亚洲精品一区av在线观看| 亚洲精品国产色婷婷电影| 90打野战视频偷拍视频| 国产精品免费视频内射| 精品久久久久久电影网| 一进一出好大好爽视频| 欧美乱码精品一区二区三区| 大香蕉久久成人网| 波多野结衣av一区二区av| 黄片播放在线免费| 日本免费a在线| 欧美日韩黄片免| 性欧美人与动物交配| 国产成年人精品一区二区 | 在线观看免费日韩欧美大片| 久久久久亚洲av毛片大全| 国产成人欧美在线观看| 啦啦啦免费观看视频1| 欧美在线一区亚洲| 成人18禁高潮啪啪吃奶动态图| 成年人免费黄色播放视频| 女人爽到高潮嗷嗷叫在线视频| 久久99一区二区三区| 亚洲中文字幕日韩| 国产国语露脸激情在线看| 涩涩av久久男人的天堂| 精品午夜福利视频在线观看一区| 国产精品香港三级国产av潘金莲| 欧美中文综合在线视频| 1024香蕉在线观看| 99国产精品99久久久久| 久久人妻熟女aⅴ| 亚洲熟女毛片儿| 女警被强在线播放| 97碰自拍视频| 成年女人毛片免费观看观看9| 搡老岳熟女国产| 成人影院久久| 久久久久久久午夜电影 | 在线免费观看的www视频| 真人做人爱边吃奶动态| 51午夜福利影视在线观看| 手机成人av网站| 大型av网站在线播放| 成人特级黄色片久久久久久久| 国产高清国产精品国产三级| 午夜福利,免费看| 99精品在免费线老司机午夜| 国产精品爽爽va在线观看网站 | 香蕉丝袜av| 老熟妇乱子伦视频在线观看| 亚洲人成电影观看| 亚洲全国av大片| 精品欧美一区二区三区在线| 久久天堂一区二区三区四区| 亚洲 欧美一区二区三区| 精品国产乱码久久久久久男人| 又黄又爽又免费观看的视频| a级毛片在线看网站| 久久久久久免费高清国产稀缺| 一级作爱视频免费观看| 日本五十路高清| 成人亚洲精品av一区二区 | 高清在线国产一区| 天堂√8在线中文| 国产熟女xx| 老鸭窝网址在线观看| 欧美亚洲日本最大视频资源| 精品一区二区三区四区五区乱码| 巨乳人妻的诱惑在线观看| 欧美大码av| 国产亚洲欧美在线一区二区| 不卡av一区二区三区| 久久精品成人免费网站| 国产高清视频在线播放一区| 首页视频小说图片口味搜索| 精品一品国产午夜福利视频| 夜夜躁狠狠躁天天躁| 99国产综合亚洲精品| 国产精品免费视频内射| 欧美人与性动交α欧美精品济南到| 免费看a级黄色片| 多毛熟女@视频| 人人妻人人爽人人添夜夜欢视频| 在线国产一区二区在线| 久久天躁狠狠躁夜夜2o2o| 亚洲国产欧美日韩在线播放| 国产亚洲欧美98| 国产精品一区二区三区四区久久 | 精品第一国产精品| 人成视频在线观看免费观看| 午夜影院日韩av| 一级毛片精品| 国产精品偷伦视频观看了|