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

    河岸緩沖帶植被布局對(duì)氮流失的影響

    2023-01-12 13:24:46陳嘉碩唐思玉
    關(guān)鍵詞:水曲柳緩沖帶銨態(tài)氮

    韓 旭,杜 崇,陳嘉碩,唐思玉

    河岸緩沖帶植被布局對(duì)氮流失的影響

    韓 旭,杜 崇※,陳嘉碩,唐思玉

    (黑龍江大學(xué)水利電力學(xué)院,哈爾濱 150006)

    河岸緩沖帶(Rriver Buffer Strips,RBSs)已被證明能有效攔截水流并去除氮,該研究旨在揭示不同配置緩沖帶對(duì)脫氮效果的影響。實(shí)地布置河岸植被緩沖帶,開(kāi)展地下徑流試驗(yàn)。緩沖帶設(shè)置3種不同植物種類(水曲柳,五葉楓,楊樹(shù)),和不同寬度(0、5、10、20、30 m),河岸坡度為3%,樹(shù)木密度為540株/hm2。在此條件下,研究各類河岸植被緩沖帶對(duì)地下徑流銨態(tài)氮、硝態(tài)氮和總氮的截留效果。結(jié)果表明:20 m寬度河岸植被緩沖帶能很好地截留各形態(tài)氮素。30 m寬度下,河岸植被緩沖帶徑流水中的銨態(tài)氮、硝態(tài)氮和總氮的截留率最高,分別為70.4%、67.7%和69.1%。在不同植物種類緩沖帶比較中,楊樹(shù)可顯著降低徑流銨態(tài)氮和總氮濃度,水曲柳可顯著降低徑流硝態(tài)氮濃度。在寬度與不同植物種類的交互關(guān)系中,20 m寬度楊樹(shù)緩沖帶對(duì)銨態(tài)氮和硝態(tài)氮的截留率最高,30 m寬度楊樹(shù)緩沖帶對(duì)總氮的截留率最高,為62.1%。研究結(jié)果可為東北地區(qū)中小型河流河岸緩沖設(shè)計(jì)最大化截留徑流氮污染物提供參考。

    氮;植被;坡度;河岸;緩沖帶;截留率;寬度;地下徑流深度

    0 引 言

    由化肥和農(nóng)藥的過(guò)量施用引起的農(nóng)業(yè)非點(diǎn)源污染是全世界河流水質(zhì)退化的主要原因之一[1]。河岸緩沖系統(tǒng)是用來(lái)控制農(nóng)業(yè)非點(diǎn)源污染和保護(hù)河流環(huán)境的主要措施之一。河岸生態(tài)森林系統(tǒng)具有重要的生態(tài)意義,因?yàn)榕c非河岸生態(tài)森林系統(tǒng)相比,河岸森林具有更加豐富的植物,可以更加有效地保護(hù)河流水質(zhì)和水生動(dòng)物[2]。在各種非點(diǎn)源污染中,農(nóng)田排放的氮污染物已成為當(dāng)前水環(huán)境治理部門(mén)需要解決的一個(gè)基本難題。河岸植被緩沖帶作為陸地生態(tài)系統(tǒng)與水生生態(tài)系統(tǒng)的連接,可有效地?cái)r截農(nóng)田地表徑流的輸入,并通過(guò)河岸緩沖帶的物理截留作用或微生物轉(zhuǎn)化作用有效地截留徑流中的氮污染物[3-5]。河岸植被緩沖帶寬度是決定河岸緩沖帶氮素截留效率的一個(gè)重要因素。王芳等[6]發(fā)現(xiàn)河岸緩沖帶的寬度決定了緩沖帶是否能夠充分發(fā)揮生態(tài)服務(wù)的功能。一般認(rèn)為河岸緩沖帶攔截氮的能力與寬度成正相關(guān)[7]。一個(gè)5 m寬的河岸緩沖帶從地表徑流中截留總氮和硝態(tài)氮的速率分別為40.2%和37.2%[8]。9.1 m和13 m寬的林草河岸緩沖帶可分別截留78%和85%地表徑流中的懸浮顆粒物[9]。中國(guó)巢湖流域23 m的草本和木本混合的河岸緩沖帶對(duì)總氮的截留率達(dá)到50%[10]。有研究表明[11],12~15 m寬的河岸植被緩沖帶是截留農(nóng)田地表徑流中氮素的最佳寬度。目前,許多研究就河岸緩沖帶最佳寬度問(wèn)題并未達(dá)成一致意見(jiàn)。在農(nóng)田地表徑流輸入河岸緩沖帶的過(guò)程中,河岸植被緩沖帶主要通過(guò)滲透、過(guò)濾、吸收、滯留和沉積等[12]作用來(lái)有效地截留農(nóng)田地表徑流中的氮污染物,具體包括土壤和植物的吸附[13],反硝化和微生物固化等[14]。而最重要的截留機(jī)制之一就是植物攔截截留。河岸緩沖帶的大型植物可以從地表徑流中吸收氮污染物[15]。其中河岸緩沖帶截留污染物的效果很大程度上取決于植物種類,因此如何確定河岸緩沖帶的最佳植被類型尤為重要[16]。研究發(fā)現(xiàn)河岸落葉林緩沖區(qū)可以使地表徑流中的氮含量降低68%[17]。然而,適當(dāng)?shù)闹脖慌渲茖?duì)于有效減少污染物也很重要。

    本文在松花江一級(jí)支流何家溝設(shè)計(jì)5種不同寬度、3種東北常見(jiàn)闊葉林樹(shù)木類型的河岸植被緩沖帶,在恒定的河岸坡度和植被密度下,研究農(nóng)田徑流中氮素截留情況,以獲得使河岸帶對(duì)氮素截留效果最佳的植被配置,以期為東北中小型河流選擇最佳氮素截留的植被河岸緩沖帶模式提供參考。

    1 材料和方法

    1.1 試驗(yàn)地概況

    本試驗(yàn)在黑龍江省哈爾濱市西南部何家溝靠山屯段進(jìn)行(126°15E′~126°34′E,45°42N′~45°44′N),研究地?zé)o霜期為168 d,年平均氣溫4.2 ℃,年均日照數(shù)為2 714.4 h,作物生長(zhǎng)季積溫為2 845.1 ℃。研究區(qū)內(nèi)降水集中在夏秋兩季,夏季為雨季,該地區(qū)平均降水天數(shù)為70.2 d,年平均降水量為524.7 mm。受河道保護(hù)政策影響以及隨著河長(zhǎng)制的全面推廣,2015年春季在試驗(yàn)地先后人工種植了3種不同類型的樹(shù)木,同時(shí)還留有1塊裸地。3個(gè)地塊的樹(shù)木密度相同,為540株/hm2,試驗(yàn)地土壤為黑鈣土,土壤糙度為1.14,在植樹(shù)后立即測(cè)量4個(gè)地塊40 cm深處土壤銨態(tài)氮、硝態(tài)氮和總氮的本底值,測(cè)量結(jié)果分別為1.15 mg/kg、2.84 mg/kg和2.68 g/kg。

    1.2 試驗(yàn)布置

    試驗(yàn)地點(diǎn)共有4塊長(zhǎng)15 m×寬30 m的人工林河岸地塊,如圖1所示。從左往右依次為裸地、水曲柳林、五葉楓林和楊樹(shù)林,其基本特征見(jiàn)表1。4種河岸植被緩沖帶地塊各設(shè)置5個(gè)寬度(0、5、10、20、30 m),各個(gè)寬度與農(nóng)田徑流流入方向平行。淋溶水收集裝置是利用內(nèi)徑為10 cm的聚氯乙烯三通管和60目(0.25 mm)濾網(wǎng)制成的試驗(yàn)淋溶管(圖2),將聚氯乙烯三通管分別截成高35 cm和高55 cm,左側(cè)淋溶水收集管道距離底部15 cm,左側(cè)管出口和底部管出口用蓋子封死,水孔的直徑為4 mm,且水孔用濾網(wǎng)包裹緊密,淋溶管在鋪設(shè)河岸種植樹(shù)木的過(guò)程中進(jìn)行埋設(shè),在每個(gè)地塊不同寬度處分別交叉埋設(shè)2組聚氯乙烯淋溶管,每組2個(gè),每組淋濾管的深度分別為20和40 cm,用于地下徑流水樣的采集。每個(gè)淋濾管之間的距離為3 m,且各組之間用PVC擋水板隔開(kāi),擋水板插入深度為50 cm(圖3)。利用施肥機(jī)在每個(gè)地塊前0~0.5 m處均勻施肥,依據(jù)當(dāng)?shù)剞r(nóng)田施肥量61 kg/hm2,所施氮磷鉀復(fù)合肥比率為24∶6∶35,在降雨前施肥,降雨產(chǎn)流后采樣,施肥與采樣間隔8 d。試驗(yàn)樣品在當(dāng)?shù)赜昙? —9月采集,3次采樣時(shí)間分別為2021年6月23日—6月25日,2021年8月6日—8月7日,2021年9月19日—21日。采樣前8天平均降水量分別為63.1、85.3和76.5 mm。

    圖1 試驗(yàn)地示意圖

    1.3 樣本采集和測(cè)定方法

    1.3.1 徑流收集方法

    使用小型抽水器抽取淋濾管中的徑流水清液,并將其裝入300 mL的塑料瓶中,將樣本帶回實(shí)驗(yàn)室存儲(chǔ)在-4~0 ℃的冰箱內(nèi),快速測(cè)定水樣中不同形態(tài)的氮,在每次取樣后,將淋濾管中的水排干,并在遠(yuǎn)離采樣區(qū)的地方排空。

    表1 林型基本特征

    圖2 PVC淋溶管

    圖3 各地塊試驗(yàn)布置圖

    1.3.2 水樣中氮素分析方法

    水樣先經(jīng)過(guò)0.45m濾膜抽濾預(yù)處理,然后,利用納氏試劑分光光度法[18]測(cè)量水樣銨態(tài)氮(NH4+-N),利用紫外分光光度法[19]測(cè)量水樣硝態(tài)氮(NO3--N),利用堿性過(guò)硫酸鉀消解紫外分光光度法[20]測(cè)量水樣總氮(Total N,TN)。759MC型號(hào)紫外分光光度計(jì)生產(chǎn)廠家為上海高致精密公司,精度±5 nm,EFD-NH型精密銨態(tài)氮測(cè)定儀生產(chǎn)廠家合肥恩帆儀器設(shè)備公司,精度≤±5%。

    氮截留率采用式(1)計(jì)算:

    R=(0-C)/0×100% (1)

    式中是河岸植被緩沖帶寬度,m;R是寬度下的氮累積截留率,%;C是寬度下的氮濃度,mg/L;0是寬度為0處徑流水的氮濃度,mg/L。

    1.4 統(tǒng)計(jì)分析

    取3次數(shù)據(jù)采樣平均值進(jìn)行分析。 Origin2018對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)繪圖,利用SPSS23.0版軟件對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析,采用單因素方差分析(One way ANOVA)并利用最小二乘(Least Significant Difference,LSD)、鄧肯(Duncan’s Multiple Range Test)多重比較方法對(duì)平均值進(jìn)行差異顯著性檢驗(yàn),差異在小于0.05時(shí)具有統(tǒng)計(jì)意義。

    2 結(jié)果與分析

    2.1 植被類型對(duì)土壤理化性質(zhì)的影響

    2015年開(kāi)始試驗(yàn)布置,2021年采樣測(cè)定各小區(qū)土壤理化性質(zhì)如表2。由表可知,種植6 a不同植被后,4個(gè)緩沖帶地塊40cm深土壤的pH差異不顯著(0.05)。水曲柳、五葉楓和楊樹(shù)緩沖帶土壤的容重小于裸地土壤的容重(0.05)。水曲柳、五葉楓和楊樹(shù)緩沖帶土壤的總孔隙度和土壤貯水量大于裸地土壤的總孔隙度和土壤貯水量(0.05)。在裸地緩沖帶中,土壤的銨態(tài)氮和硝態(tài)氮質(zhì)量分別為5.47和7.86 mg/kg,總氮的質(zhì)量為4.73 g/kg,高于銨態(tài)氮、硝態(tài)氮和總氮的本底值,這可能是因?yàn)? a中在施有過(guò)量氮肥的農(nóng)田在多次降雨形成下形成農(nóng)田徑流,含有氮素的農(nóng)田地表徑流在流經(jīng)裸地緩沖帶后,使得裸地緩沖帶土壤中氮素質(zhì)量升高。由于植被類型的不同,植被在吸收氮素方面存在差異。水曲柳、五葉楓和楊樹(shù)河岸緩沖帶土壤中銨態(tài)氮、硝態(tài)氮和總氮質(zhì)量與裸地緩沖帶差異顯著(<0.05),這說(shuō)明在裸地上種植植被可有效降低氮素含量。而楊樹(shù)河岸緩沖帶銨態(tài)氮含量顯著低于水曲柳和五葉楓河岸緩沖帶(<0.05),這可能是由于楊樹(shù)的林齡相較于水曲柳和五葉楓的林齡要小,需要更多銨態(tài)氮含量來(lái)維持自身生長(zhǎng)造成的。

    表2 2021年土壤理化性質(zhì)

    注:同行不同小寫(xiě)字母表示各林型間氮素差異顯著(<0.05)。

    Note: Different lowercase letters in the same row indicate significant differences among forest types (<0.05).

    2.2 植被類型、土壤深度與緩沖帶寬度對(duì)徑流氮的交互影響

    對(duì)植被類型、土壤深度和緩沖帶寬度因素進(jìn)行方差分析,結(jié)果見(jiàn)表3。植被類型、緩沖帶寬度對(duì)銨態(tài)氮、硝態(tài)氮和總氮影響顯著;而各因素交互影響不顯著。又從表中看出,地下徑流深度與河岸緩沖帶寬度和植被類型無(wú)顯著交互影響,因而接下來(lái)的分析中將20和40 cm土壤深度的氮素值進(jìn)行平均處理,重點(diǎn)分析不同寬度、不同植被類型和不同寬度×不同植被類型對(duì)氮素截留影響。

    表3 不同因素對(duì)徑流氮影響的方差分析P值結(jié)果

    注:<0.05,顯著。

    Note: Significance (< 0.05).

    2.3 不同寬度河岸緩沖帶對(duì)氮素截留的差異

    2.3.1 不同寬度河岸帶對(duì)NH4+-N截留的影響

    如圖4所示,銨態(tài)氮濃度隨著河岸植被緩沖帶寬度的增加而降低。寬度0、5、10 、20 m河岸植被緩沖帶徑流銨態(tài)氮濃度分別為0.681、0.375、0.270和0.214 mg/L,差異顯著(<0.05);寬度30 m河岸植被緩沖帶的徑流銨態(tài)氮截留率為70.4%,與20 m河岸植被緩沖帶截留率無(wú)顯著差異(>0.05),二者顯著高于5和10 m寬河岸植被緩沖帶銨態(tài)氮截留效率(<0.05)。5、10和20 m寬度河岸植被緩沖帶銨態(tài)氮截留率分別為48.9%、63.3%和69.6%,差異顯著(<0.05),而在20~30 m寬度河岸緩沖帶中,對(duì)銨態(tài)氮截留率相較于5~20 m來(lái)說(shuō)截留效率放緩。

    a. 銨態(tài)氮濃度a. Ammonium nitrogen concentrationb. 銨態(tài)氮截留率b. Ammonium nitrogen retention rate

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

    Note: Different letters represent significant differences (< 0.05). Same as below.

    圖4 不同緩沖帶寬度對(duì)徑流中銨態(tài)氮的影響

    Fig.4 Effects of different buffer zone widths on ammonium nitrogen in runoff

    2.3.2 不同寬度河岸帶對(duì)NO3--N截留的影響

    由圖5所示,寬度0、5、10和20 m的徑流硝態(tài)氮濃度分別為0.792、0.546、0.371和0.290 mg/L,差異顯著(<0.05),寬度5、10、20和30 m河岸植被緩沖帶的硝態(tài)氮截留率分別為34.7%、53.7%、64.2%和67.7%,30 m與20 m河岸植被緩沖帶硝態(tài)氮截留率無(wú)顯著差異(>0.05),但顯著高于5 和10 m寬河岸植被緩沖帶硝態(tài)氮截留效率(<0.05)。隨著寬度的增加,硝態(tài)氮的濃度降低,截留率上升[21-22],寬度30 m的河岸植被緩沖帶具有最高的地下徑流硝態(tài)氮截留率,這樣的結(jié)果與銨態(tài)氮類似。

    a. 硝態(tài)氮濃度a. Nitrate nitrogen concentrationb. 硝態(tài)氮截留率b. Nitrate nitrogen retention rate

    2.3.3 不同寬度河岸帶對(duì)總氮(TN)截留的影響

    如圖6所示,在徑流水中,總氮濃度隨著河岸植被緩沖帶寬的增加而降低,這說(shuō)明土壤在農(nóng)田徑流水下滲的過(guò)程中截留吸附了徑流中的部分氮污染物,這與黃玲玲[23]研究中的土壤的吸附作用使得徑流中氮濃度降低相一致。寬度0、5、10、20和30 m緩沖帶徑流總氮濃度分別為4.137、2.581、2.343、1.486和1.371 mg/L,其中寬度20和30 m緩沖帶徑流總氮濃度與寬度0、5和10 m緩沖帶徑流總氮濃度差異顯著(<0.05)。寬度5、10、20 和30 m緩沖帶徑流總氮截留率分別為39.8%、42.1%、64.7%和69.1%,并且20和30 m緩沖帶徑流總氮截留率與10和20 m緩沖帶徑流截留率差異顯著(<0.05)。

    a. 總氮濃度a. Total nitrogen concentrationb. 總氮截留率b. Total nitrogen retention rate

    結(jié)合河岸植被緩沖帶3種氮素分析表明:徑流中銨態(tài)氮,硝態(tài)氮和總氮的濃度隨著河岸植被緩沖帶寬度的增加而逐漸減少,這與朱曉成等[24-25]研究的結(jié)果基本一致,這可能是當(dāng)河岸緩沖帶的寬度逐漸增大時(shí),使得徑流氮污染物與植被緩沖帶的作用時(shí)間邊長(zhǎng),促進(jìn)了地表徑流下滲,進(jìn)而促進(jìn)土壤和植物對(duì)氮素的截留和吸收。在氮素的截留率方面,20和30 m寬的河岸植被緩沖帶的銨態(tài)氮,硝態(tài)氮和總氮的截留率差異不顯著(>0.05),可見(jiàn),寬度20 m的河岸植被緩沖帶就能夠?qū)?種氮素進(jìn)行很好截留。

    2.4 不同植被類型對(duì)河岸緩沖帶氮素截留的差異

    如圖7所示,以裸地緩沖帶地塊為對(duì)照,來(lái)分析水曲柳、五葉楓和楊樹(shù)植被緩沖帶地塊在銨態(tài)氮、硝態(tài)氮和總氮截留差異。裸地、水曲柳、五葉楓和楊樹(shù)緩沖帶的銨態(tài)氮濃度分別為0.428、0.372、0.364和0.269 mg/L。水曲柳、五葉楓和楊樹(shù)緩沖帶的銨態(tài)氮濃度與裸地植被緩沖帶差異顯著(<0.05),其中楊樹(shù)緩沖帶的銨態(tài)氮濃度顯著低于水曲柳和五葉楓緩沖帶(<0.05),這可能是因?yàn)橄鄬?duì)于水曲柳與五葉楓來(lái)說(shuō),楊樹(shù)的樹(shù)齡更小,維持自身生長(zhǎng)所需吸收的銨態(tài)氮更多造成的。

    由圖7b可知,裸地、水曲柳、五葉楓和楊樹(shù)緩沖帶的硝態(tài)氮濃度分別為0.476、0.322、0.414和0.392 mg/L。各植被緩沖帶硝態(tài)氮濃度與裸地緩沖帶硝態(tài)氮濃度差異顯著(<0.05),而水曲柳緩沖帶硝態(tài)氮濃度顯著低于五葉楓和楊樹(shù)緩沖帶(<0.05),這說(shuō)明相對(duì)于五葉楓和楊樹(shù)緩沖帶,水曲柳緩沖帶具有吸收更多硝態(tài)氮的特性。

    由圖7c可知,裸地、水曲柳、五葉楓和楊樹(shù)緩沖帶的總氮濃度分別為2.588、1.352、1.896和1.117 mg/L,水曲柳和楊樹(shù)緩沖帶總氮濃度顯著低于五葉楓河岸緩沖帶(<0.05),且水曲柳與楊樹(shù)緩沖帶總氮濃度差異不顯著(>0.05),這說(shuō)明水曲柳與楊樹(shù)緩沖帶相比于五葉楓緩沖帶具有降低徑流總氮濃度的特性。

    結(jié)合上述分析可知,水曲柳、五葉楓和楊樹(shù)河岸植被緩沖帶可有效降低徑流氮素濃度,楊樹(shù)緩沖帶在降低徑流銨態(tài)氮濃度方面效果最好,水曲柳緩沖帶在降低硝態(tài)氮濃度方面效果最好,水曲柳與楊樹(shù)緩沖帶相較于五葉楓緩沖帶在降低徑流總氮濃度方面效果更好,且水曲柳與楊樹(shù)緩沖帶在降低徑流總氮濃度方面效果類似。

    2.5 緩沖帶寬度和植被類型對(duì)河岸緩沖帶氮素截留的差異

    圖8顯示了植被類型(水曲柳,五葉楓,楊樹(shù))的河岸緩沖帶對(duì)地下徑流中銨態(tài)氮的截留效果。當(dāng)徑流經(jīng)過(guò)河岸緩沖區(qū)后,徑流中銨態(tài)氮濃度顯著降低,在不同植被類型的河岸緩沖帶對(duì)銨態(tài)氮截留率存在明顯差異。在寬度5、10和20 m河岸植被緩沖帶下,楊樹(shù)河岸緩沖帶的銨態(tài)氮截留率分別為62.3%、75.7%和78.1%,與水曲柳、五葉楓和裸地河岸緩沖帶的銨態(tài)氮截留率差異顯著(<0.05)。而在30 m寬度的河岸植被緩沖帶下,楊樹(shù)與水曲柳河岸緩沖帶的銨態(tài)氮截留率差異不顯著(>0.05)。在各寬度河岸植被緩沖帶中,楊樹(shù)河岸緩沖帶是銨態(tài)氮截留效率最高的緩沖帶,這主要由于試驗(yàn)樣地的楊樹(shù)處于幼年期,相較于五葉楓林和水曲柳林來(lái)說(shuō),幼年楊樹(shù)林需要吸收更多養(yǎng)分,可能對(duì)徑流水中的銨態(tài)氮有較高的截留率。20 m寬度的楊樹(shù)河岸緩沖帶銨態(tài)氮截留率最大為78.1%,比水曲柳,五葉楓和裸地河岸緩沖帶截留率分別高17.9%,26.8%和29.4%。

    a. 銨態(tài)氮濃度a. Ammonium nitrogen concentrationb. 硝態(tài)氮濃度b. Nitrate nitrogen concentrationc. 總氮濃度c. Total nitrogen concentration

    注:BL為裸地、AT為水曲柳、FLM為五葉楓、PL為楊樹(shù)。

    Note: BL is bare land, AT is ashtree, FLM is five-leaf maple, PL is poplar.

    圖7 不同緩沖帶植被類型對(duì)徑流中氮素濃度的影響

    Fig.7 Effect of different buffer zone vegetation types on nitrogen concentration in runoff

    a. 銨態(tài)氮a. Ammonium nitrogenb. 硝態(tài)氮b. Nitrate nitrogenc. 總氮c. Total nitrogen

    不同植被類型的河岸緩沖帶對(duì)徑流中硝態(tài)氮的截留效果如圖8b所示,在寬度5和30 m的河岸緩沖帶中,水曲柳河岸緩沖帶對(duì)徑流硝態(tài)氮的截留率最高,分別為45%和65.3%,寬度10和20 m的楊樹(shù)河岸緩沖帶對(duì)徑流硝態(tài)氮的截留率分別為58.6%和64.2%,寬度20 m的五葉楓河岸緩沖帶的徑流硝態(tài)氮截留率為52.1%,相較于寬度10 m的五葉楓河岸緩沖帶的硝態(tài)氮截留率53.6%略有下降。楊樹(shù)河岸緩沖帶在20 m寬度時(shí),徑流硝態(tài)氮截留率最大,為64.2%,與五葉楓和裸地河岸緩沖帶截留率差異顯著(<0.05),但與水曲柳差異不顯著(>0.05),其中楊樹(shù)河岸緩沖帶比水曲柳、五葉楓和裸地河岸緩沖帶硝態(tài)氮截留率分別高2.1%、12.1%和23.6%,截留率從高到低依次為楊樹(shù)林、水曲柳林、五葉楓林、裸地。30 m寬度下,水曲柳河岸緩沖帶與五葉楓河岸緩沖帶徑流硝態(tài)氮截留率分別為65.3%和54.1%,差異顯著(<0.05),而在其他寬度下,二者徑流硝態(tài)氮平均截留率并沒(méi)有顯著差異(>0.05)。在5和10 m寬的河岸緩沖帶上,水曲柳緩沖帶與五葉楓緩沖帶的徑流硝態(tài)氮平均截留率差異不顯著(>0.05)。

    對(duì)于總氮來(lái)說(shuō),如圖8c所示,在寬度5 m河岸植被緩沖帶下,水曲柳、五葉楓和楊樹(shù)河岸緩沖帶的總氮截留率分別為39.7%、41.2%和47.6%,與裸地河岸緩沖帶差異顯著(<0.05)。在寬度5、20和30 m河岸植被緩沖帶中,楊樹(shù)河岸緩沖帶有著最高的徑流總氮截留率,分別為47.6%、60.5%和62.1%。在寬度10 m河岸植被緩沖帶下,楊樹(shù)河岸緩沖帶與水曲柳河岸緩沖帶的徑流總氮截留率相近,分別為54.6%和54.3%,并且差異不顯著(>0.05)。在30 m處,水曲柳、五葉楓和楊樹(shù)河岸緩沖帶的總氮截留率與裸地河岸緩沖帶總氮截留率存在顯著差異(<0.05)。在寬度10和20 m的河岸植被緩沖帶中,水曲柳林徑流總氮截留率與五葉楓林差異不顯著(>0.05)。

    王金亮等[26]研究表明,在較緩坡度上,林地坡地相比于旱坡地來(lái)說(shuō),林地坡地阻力較大。所以,在一定坡度下,河岸緩沖帶的植被通過(guò)輸送和吸收溶解氧為微生物提供棲息地,降低徑流流速,疏松土壤等機(jī)制來(lái)預(yù)防和控制非點(diǎn)源污染[27-28]。這種機(jī)制在每個(gè)反應(yīng)中可能會(huì)因河岸植被緩沖帶的條件不同而有所不同,某種機(jī)制可能對(duì)某種的氮污染物非常有效[29]。不同類型的植被造成植物組織組成,根系類型和活動(dòng)強(qiáng)度的差異,導(dǎo)致不同植被類型的河岸緩沖帶對(duì)氮污染物的吸收和轉(zhuǎn)化存在差異[30]。本研究中,3種河岸植被緩沖帶對(duì)銨態(tài)氮,硝態(tài)氮和總氮的截留能力均高于裸地,這表明不同植被類型的緩沖帶對(duì)徑流水中氮污染物有著很好的截留效果。楊樹(shù)林河岸緩沖帶對(duì)徑流水中銨態(tài)氮和總氮的平均截留率較高,除寬度30 m外,在其他寬度下,楊樹(shù)林河岸緩沖帶銨態(tài)氮的平均截留率與水曲柳林和五葉楓林河岸緩沖帶銨態(tài)氮平均截留率差異顯著(<0.05),對(duì)于總氮來(lái)說(shuō),各寬度下楊樹(shù)林總氮平均截留率與水曲柳林總氮平均截留率差異不顯著(>0.05)。各寬度楊樹(shù)林緩沖帶對(duì)徑流水中硝態(tài)氮平均截留率的總和略高于各寬度水曲柳林緩沖帶對(duì)徑流水中硝態(tài)氮平均截留率的總和,這可能是因?yàn)闂顦?shù)的根系比水曲柳的更深,能夠吸收更多淋溶到深層土壤中的硝態(tài)氮,且試驗(yàn)點(diǎn)的樹(shù)種處于幼年期,對(duì)養(yǎng)分的需求量大。

    3 結(jié) 論

    布置3種不同植物種類(水曲柳,五葉楓,楊樹(shù))和不同河岸植被緩沖帶寬度(5、10、20、30 m),開(kāi)展地下徑流試驗(yàn),研究各類河岸植被緩沖帶對(duì)地下徑流銨態(tài)氮、硝態(tài)氮和總氮的截留效果,結(jié)果表明:有植被覆蓋的河岸緩沖帶截留氮素效果顯著高于無(wú)植被覆蓋的河岸緩沖帶(裸地),寬度20 m的河岸緩沖帶對(duì)由模擬農(nóng)田輸入徑流氮的污染物的截留效果顯著,滿足脫氮的預(yù)期目標(biāo)。在3種不同河岸植被緩沖帶中,楊樹(shù)林對(duì)銨態(tài)氮的截留率最高,且與五葉楓林,水曲柳林截留差異顯著(<0.05)。對(duì)于硝態(tài)氮來(lái)說(shuō),寬度10和20 m的楊樹(shù)河岸緩沖帶對(duì)徑流硝態(tài)氮的平均截留效率與五葉楓林于水曲柳林的徑流硝態(tài)氮的平均截留效率差異顯著(<0.05),但在寬度5和30 m下,徑流硝態(tài)氮的平均截留效率差異不顯著(>0.05)。在總氮的情況下,除寬度10 m楊樹(shù)河岸緩沖帶徑流總氮的截留率略低于水曲柳河岸緩沖帶徑流總氮的截留率,其他寬度下,楊樹(shù)河岸緩沖帶徑流總氮的截留率均最高。綜合以上結(jié)果,結(jié)合東北地區(qū)生態(tài)經(jīng)濟(jì)條件來(lái)說(shuō),20 m寬度楊樹(shù)林河岸緩沖帶對(duì)于氮污染物的截留達(dá)到最大化的效果。

    [1] Hénault-Ethier L, Gómes Marcelo Pedrosa, Lucotte M, et al. High yields of riparian buffer strips planted with Salix miyabena 'SX64' along field crops in Québec, Canada[J]. Biomass & Bioenergy, 2017, 10(5): 219-229.

    [2] Luke S H, Slade E M, Gray C L, et al. Riparian buffers in tropical agriculture: Scientific support, effectiveness and directions for policy[J]. Journal of Applied Ecology, 2019, 56(1): 85-92.

    [3] 劉瑞霞,王立陽(yáng),孫菲,等. 以農(nóng)業(yè)面源污染阻控為目標(biāo)的河流生態(tài)緩沖帶研究進(jìn)展[J]. 環(huán)境工程學(xué)報(bào),2022,16(1):25-39.

    Liu Ruixia, Wang Liyang, Sun Fei, et al. Research progress in riverine ecological buffer zone for control of agricultural non-point source pollution[J]. Chinese Journal of Environmental Engineering, 2022, 16(1): 25-39. (in Chinese with English abstract)

    [4] 劉贏男,楊帆,韓輝,等. 阿什河流域河岸植被緩沖帶寬度規(guī)劃[J]. 國(guó)土與自然資源研究,2018(6):46-48.

    Liu Yingnan, Yang Fan, Han Hui, et al. Width planning of riparian vegetation buffer zone in Ashe River Basin[J]. Land and Natural Resources Research, 2018(6): 46-48. (in Chinese with English abstract)

    [5] Dlamini J C, Cárdenas L, Tesfamariam E H, et al. Soil methane (CH4) fluxes in cropland with permanent pasture and riparian buffer strips with different vegetation[J]. Journal of Plant Nutrition and Soil Science, 2021, 185(1): 132-144.

    [6] 王芳,汪耀龍,謝祥財(cái). 生態(tài)學(xué)價(jià)值視角下的城市河流綠道寬度研究進(jìn)展[J]. 中國(guó)城市業(yè),2019,17(1):57-61.

    Wang Fang, Wang Yaolong, Xie Xiangcai. Research progress of urban river greenway width from the perspective of ecological value[J]. China Urban Industry, 2019, 17(1): 57-61. (in Chinese with English abstract)

    [7] Lyu C J, Li X J, Yuan P, et al. Nitrogen retention effect of riparian zones in agricultural areas: A meta-analysis[J]. Journal of CleanerProduction, 2021, 315: 128143.

    [8] 張鴻齡,李天嬌,趙志芳,等. 遼河河岸植被緩沖帶構(gòu)建及其對(duì)固體顆粒物和氮阻控能力[J]. 生態(tài)學(xué)雜志,2020,39(7):2185-2192.

    Zhang Hongling, Li Tianjiao, Zhao Zhifang, et al. Construction of vegetation buffer zone along the Liaohe river bank and its resistance and control ability to solid particulate matter and nitrogen[J]. Chinese Journal of Ecology, 2020, 39(7): 2185-2192. (in Chinese with English abstract)

    [9] 湯家喜,何苗苗,王道涵,等. 河岸緩沖帶對(duì)地表徑流及懸浮顆粒物的阻控效應(yīng)[J]. 環(huán)境工程學(xué)報(bào),2016, 10(5):2747-2755.

    Tang Jiaxi, He Miaomiao, Wang Daohan, et al. Effect of riparian buffer zone on surface runoff and suspended particulate matter[J]. Chinese Journal of Environmental Engineering, 2016, 10(5): 2747-2755. (in Chinese with English abstract)

    [10] Cao X, Song C, Xiao J, et al. The optimal width and mechanism of riparian buffers for storm water nutrient rem oval in the Chinese Eutrophic Lake Chaohu Watershed[J]. Water, 2018, 10(10): 1-11.

    [11] Kiffney P M, Bull R J P. Responses of periphyton and insects to experim ental manipulation of riparian buffer width along forest streams[J]. Journal of Applied Ecology, 2010, 40(6): 1060-1076.

    [12] Li Y, Chen Z, Lou H, et al. Denitrification controls in urban riparian soils: Implications for reducing urban nonpoint source nitrogen pollution[J]. Environmental Science & Pollution Research, 2014, 21(17): 10174-10185.

    [13] 胡海波,鄧文斌,王霞. 長(zhǎng)江流域河岸植被緩沖帶生態(tài)功能及構(gòu)建技術(shù)研究進(jìn)展[J]. 浙江農(nóng)林大學(xué)學(xué)報(bào),2022,39(1):214-222.

    Hu Haibo, Deng Wenbin, Wang Xia. Research progress on ecological function and construction technology of riparian vegetation buffer zone in the Yangtze River Basin[J]. Journal of Zhejiang A & F University, 2022, 39(1): 214-222. (in Chinese with English abstract)

    [14] 孫金偉,許文盛. 河岸植被緩沖帶生態(tài)功能及其過(guò)濾機(jī)理的研究進(jìn)展[J]. 長(zhǎng)江科學(xué)院院報(bào),2017,34(3):40-44.

    Sun Jinwei, Xu Wensheng. Research progress on ecological function and filtration mechanism of riparian vegetation buffer zone[J]. Journal of Yangtze River Scientific Research Institute, 2017, 34(3): 40-44. (in Chinese with English abstract)

    [15] Johnson R K, Burdon F J, et al. Forested riparian buffers change the taxonomic and functional composition of stream invertebrate communities in agricultural catchments[J]. Water, 2021, 13(8): 1-17.

    [16] 紀(jì)欽陽(yáng). 九龍江北溪河岸植被緩沖帶構(gòu)建實(shí)驗(yàn)區(qū)氮、磷削減研究[D]. 福州:福建師范大學(xué),2016.

    Ji Qinyang. Research on Nitrogen and Phosphorus Reduction in the Experimental Area of Constructing the Vegetation Buffer Zone on the North River Bank of Jiulong River[D]. Fuzhou: Fujian Normal University, 2016. (in Chinese with English abstract)

    [17] Li R, Kuo Y M. Effects of shallow water table depth on vegetative filter strips retarding transport of nonpoint source pollution in controlled flume experiments[J]. 2021, 15(1): 163-175.

    [18] 張臻宇. 連續(xù)流動(dòng)分析法和納氏試劑分光光度法測(cè)定地表水氨氮的比較研究[J]. 環(huán)境與發(fā)展,2020,32(11):87-88.

    Zhang Zhenyu. Comparative study on determination of ammonia nitrogen in surface water by continuous flow analysis and Nessler's reagent spectrophotometry[J]. Environment and Development, 2020, 32(11): 87-88. (in Chinese with English abstract)

    [19] 李欣然,何巖,朱瑾,等. 丙烯基硫脲和鉬酸鈉對(duì)紫外分光光度法和酚二磺酸光度法測(cè)定硝態(tài)氮的影響[J]. 化學(xué)世界,2019,60(3):149-155.

    Li Xinran, He Yan, Zhu Jin, et al. Effect of propylene thiourea and sodium molybdate on the determination of nitrate nitrogen by ultraviolet spectrophotometry and phenol disulfonic acid spectrophotometry[J]. Chemical World, 2019, 60(3): 149-155. (in Chinese with English abstract)

    [20] 谷東杰,劉倩. 堿性過(guò)硫酸鉀消解紫外分光光度法測(cè)定水中總氮[J]. 山東化工,2020,49(11):103-105.

    Gu Dongjie, Liu Qian. Determination of total nitrogen in water by ultraviolet spectrophotometry with alkaline potassium persulfate digestion[J]. Shandong Chemical Industry, 2020, 49(11): 103-105. (in Chinese with English abstract)

    [21] 王曉琪,姚媛媛,陳寶成,等. 淹水條件硝態(tài)氮和銨態(tài)氮配施對(duì)水稻生長(zhǎng)與土壤養(yǎng)分的影響[J]. 土壤,2020,52(2):254-261.

    Wang Xiaoqi, Yao Yuanyuan, Chen Baocheng, et al. Effects of combined application of nitrate and ammonium nitrogen on rice growth and soil nutrients under flooding conditions[J]. Soils, 2020, 52(2): 254-261. (in Chinese with English abstract)

    [22] 李懷恩,鄧娜,楊寅群,等. 植被過(guò)濾帶對(duì)地表徑流中污染物的凈化效果[J]. 農(nóng)業(yè)工程學(xué)報(bào),2010,26(7):81-86.

    Li Huaien, Deng Na, Yang Yinqun, et al. Clarification efficiency of vegetative filter strips to several pollutants in surface runoff[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(7): 81-86. (in Chinese with English abstract)

    [23] 黃玲玲. 竹林河岸帶對(duì)氮磷截留轉(zhuǎn)化作用的研究[D]. 北京:中國(guó)林業(yè)科學(xué)研究院,2009.

    Huang Lingling. Study on the Interception An Transformation of Nitrogen and Phosphorus in the Riparian Zone of Bamboo Forest[D]. Beijing: Chinese Academy of Forestry Sciences, 2009. (in Chinese with English abstract)

    [24] 朱曉成,吳永波,余昱瑩,等. 太湖喬木林河岸植被緩沖帶截留氮素效率[J]. 浙江農(nóng)林大學(xué)學(xué)報(bào),2019,36(3):565-572.

    Zhu Xiaocheng, Wu Yongbo, Yu Yuying, et al. Nitrogen interception efficiency of riparian vegetation buffer zone of arbor forest in Taihu Lake[J]. Journal of Zhejiang A & F University, 2019, 36(3): 565-572. (in Chinese with English abstract)

    [25] Lv J, Wu Y. Nitrogen removal by different riparian vegetation buffer strips with different stand densities and widths[J]. Water Science & Technology Water Supply, 2021, 11: 1-16.

    [26] 王金亮,陳成龍,倪九派,等. 小流域農(nóng)業(yè)面源污染阻力評(píng)價(jià)及“源-匯”風(fēng)險(xiǎn)空間格局[J]. 農(nóng)業(yè)工程學(xué)報(bào),2018,34(10):216-224.

    Wang Jinliang, Chen Chenglong, Ni Jiupai, et al. Resistance evaluation and “source-sink” risk spatial pattern of agricultural non-point source pollution in small catchment[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(10): 216-224. (in Chinese with English abstract)

    [27] 楊世琦,邢磊,劉宏元,等. 植物籬埂壟向區(qū)田技術(shù)對(duì)坡耕地水土和氮磷流失控制研究[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(22):209-215.

    Yang Shiqi, Xing Lei, Liu Hongyuan, et al. Effect of reducing runoff, sediment, soil nitrogen and phosphorus losses in sloping farmland based on short ridge of clover hedgerow with ridge tillage[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(22): 209-215. (in Chinese with English abstract)

    [28] 劉燕,夏品華,鄭鈞寧,等. 河岸緩沖帶植物配置模式對(duì)面源污染物的凈化效果[J]. 貴州農(nóng)業(yè)科學(xué),2014(10):248-251.

    Liu Yan, Xia Pinhua, Zheng Junning, et al. Purification effect of surface source pollutants by plant configuration mode in riparian buffer zone[J]. Guizhou Agricultural Science, 2014(10): 248-251. (in Chinese with English abstract)

    [29] 孫東耀. 不同植被河岸緩沖帶系統(tǒng)對(duì)氮磷污染物的攔截過(guò)程研究[D]. 福建:福建師范大學(xué),2018.

    Sun Dongyao. Study on Interception Process of Nitrogen and Phosphorus Pollutants by Different Vegetation Riparian Buffer Zone Systems[D]. Fujian: Fujian Normal University, 2018. (in Chinese with English abstract)

    [30] 孫東耀,仝川,紀(jì)欽陽(yáng),等.不同類型植被河岸緩沖帶對(duì)模擬徑流及總磷的消減研究[J].環(huán)境科學(xué)學(xué)報(bào),2018,38(6):2393-2399.

    Sun Dongyao, Tong Chuan, Ji Qinyang, et al. Reduction of simulated runoff and total phosphorus in different vegetation riparian buffer [J].Acta Scientiae Circumstantiae, 2018,38(6): 2393-2399. (in Chinese with English abstract)

    Effects of riparian buffer zone vegetation layout on nitrogen loss

    Han Xu, Du Chong※, Chen Jiashuo, Tang Siyu

    (,,150006,)

    To study the nitrogen retention efficiency of riparian vegetation buffers of small and medium-sized rivers, the common broad-leaved species in Northeast China were slected on the slope of Hejiagou, Harbin in Northeast China. Four riparian planted woodlots of 15 m long by 30 m wide were prepared for the different species experiment, namely bare ground, water willow or ash, five-leaf maple, and poplar riparian vegetation buffers. The bare ground riparian buffers served as controls. The effects of different riparian buffer widths (0, 5, 10, 20, and 30 m) and different plant species above on nitrogen retention in subsurface runoff at 20 and 40 cm depth were analyzed under the condition of riparian slope (3%) and tree density (540 trees/hm2). At the test sites, PVC pipes with a diameter of 10 cm and a 4-mm partial subsurface runoff intake hole on the left side were used to collect subsurface runoff at soil depths of 20 cm and 40 cm. The PVC pipes were buried when the trees were planted. At the four riparian buffer plots from 0 to 0.5 m, the compound fertilizer was evenly applied using a small fertilizer spreader to simulate nitrogen loss conditions. The fertilizer was applied before rainfall, and the ratio of flow-producing compound fertilizer collected after rainfall was 24:6:35 for nitrogen: phosphorus: potassium, with fertilizer application rate of 61 kg/hm2. When sampling, the upper clear liquid in the tube was pumped out with a small pump, and then all the turbid liquid in the tube was pumped out and drained at the far end of the test site to avoid affecting the experimental results, the clear liquid was sampled into 300 mL plastic bottles and stored in a -4-0 ℃ refrigerator for determination of the nitrogen content of the water samples. The ammonium nitrogen in the runoff was determined spectrophotometrically by the nano reagent, the nitrate nitrogen was determined by a UV spectrophotometry, and the total nitrogen was determined by a UV spectrophotometry using the alkaline potassium persulfate elimination method. The statistical analysis results showed that the depth of subsurface runoff did not significantly interact with riparian buffer width and vegetation type. Therefore, the nitrogen values at 20 and 40 cm soil depth were averaged for the followed analysis on the effects of different widths, different vegetation types and different widths × different vegetation types on nitrogen retention. The results showed that the concentrations of runoff ammonium and nitrate nitrogen gradually decreased in riparian vegetation buffers of 0, 5, 10, and 20 m in width. Thus, the riparian vegetation buffer with a width of 20 m could retain all forms of nitrogen well. The 30 m width riparian vegetation buffer had the highest retention of ammonium nitrogen, nitrate nitrogen, and total nitrogen in runoff with 70.4%, 67.7%, and 69.1%, respectively. in the comparison of different plant species buffer zones, poplar could significantly reduce the runoff ammonium nitrogen and total nitrogen concentrations, and water willow could significantly reduce the runoff, the poplar significantly reduced runoff ammonium and total nitrogen concentrations, and ash significantly reduced runoff nitrate-nitrogen concentrations. In the interaction between width and plant species, poplar buffer strips of 20 m width had the highest retention rate of ammonium nitrogen and nitrate nitrogen, and the poplar buffer strips of 30 m width had the highest retention rate of 62.1% for total nitrogen. The results of the study can provide valuable information for the design of riparian buffers to maximize the retention of runoff nitrogen pollutants in small and medium-sized rivers in Northeast China.

    nitrogen; vegetation; slope; streambank; buffer zone; retention rate; width; depth of subsurface runoff

    10.11975/j.issn.1002-6819.2022.16.019

    S153.6

    A

    1002-6819(2022)-16-0172-08

    韓旭,杜崇,陳嘉碩,等. 河岸緩沖帶植被布局對(duì)氮流失的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2022,38(16):172-179.doi:10.11975/j.issn.1002-6819.2022.16.019 http://www.tcsae.org

    Han Xu, Du Chong, Chen Jiashuo, et al. Effects of riparian buffer zone vegetation layout on nitrogen loss[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(16): 172-179. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.16.019 http://www.tcsae.org

    2022-03-16

    2022-07-26

    黑龍江省高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(2018-KYYWF-1570);國(guó)家自然科學(xué)基金(52109055)

    韓旭,研究方向?yàn)樯鷳B(tài)水文。Email:1329626798@qq.com

    杜崇,博士,副教授,研究方向?yàn)檗r(nóng)業(yè)水土資源與環(huán)境保護(hù)。Email:duchong@hlju.edu.cn

    猜你喜歡
    水曲柳緩沖帶銨態(tài)氮
    寶龍店水曲柳無(wú)性系生長(zhǎng)早期評(píng)價(jià)
    不同質(zhì)地土壤銨態(tài)氮吸附/解吸特征
    不同鹽堿化土壤對(duì)NH+4吸附特性研究
    河岸人工林緩沖帶對(duì)徑流水磷素的截留效果
    淺談水曲柳營(yíng)林的速產(chǎn)豐產(chǎn)技術(shù)要點(diǎn)
    水曲柳和落葉松人工純林與混交林的碳儲(chǔ)量
    森林工程(2018年4期)2018-08-04 03:23:10
    林區(qū)水曲柳營(yíng)林的速生豐產(chǎn)技術(shù)要點(diǎn)
    太湖流域流經(jīng)不同類型緩沖帶入湖河流秋、冬季氮污染特征?
    坡耕地作物緩沖帶對(duì)徑流和氮流失的影響分析
    濱岸緩沖帶去除農(nóng)業(yè)面源污染技術(shù)研究
    国产精品一及| 最近最新免费中文字幕在线| 成年免费大片在线观看| 在线国产一区二区在线| 又粗又爽又猛毛片免费看| 在线免费观看的www视频| 一夜夜www| 1024手机看黄色片| 国产精品电影一区二区三区| 亚洲熟女毛片儿| 精华霜和精华液先用哪个| 欧美性猛交黑人性爽| 久久久色成人| 天堂影院成人在线观看| 老司机在亚洲福利影院| 婷婷六月久久综合丁香| 免费搜索国产男女视频| 欧美中文综合在线视频| 国产精品99久久99久久久不卡| 午夜福利18| 国产亚洲欧美在线一区二区| 一本久久中文字幕| 日本一本二区三区精品| 亚洲一区二区三区色噜噜| 久久婷婷人人爽人人干人人爱| 一夜夜www| 成年免费大片在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品久久久av美女十八| 久久久久久久精品吃奶| 99久久综合精品五月天人人| 国内精品一区二区在线观看| 精华霜和精华液先用哪个| 国产精品av久久久久免费| 可以在线观看的亚洲视频| 国产精品一区二区免费欧美| 亚洲av片天天在线观看| 国产高潮美女av| 他把我摸到了高潮在线观看| 免费观看的影片在线观看| 1024手机看黄色片| 欧美国产日韩亚洲一区| 亚洲av片天天在线观看| 亚洲中文字幕一区二区三区有码在线看 | 免费av毛片视频| 欧美成人性av电影在线观看| 首页视频小说图片口味搜索| 国产精品乱码一区二三区的特点| 欧美日本亚洲视频在线播放| 国产亚洲精品av在线| 两个人的视频大全免费| 亚洲av日韩精品久久久久久密| 午夜精品久久久久久毛片777| 中文字幕最新亚洲高清| 身体一侧抽搐| 长腿黑丝高跟| 99久久国产精品久久久| 国产成人精品久久二区二区免费| 亚洲精品一区av在线观看| 亚洲无线在线观看| 欧美3d第一页| 曰老女人黄片| 午夜免费观看网址| 老鸭窝网址在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲色图av天堂| 国产真人三级小视频在线观看| 91av网站免费观看| 国产精品久久电影中文字幕| 亚洲中文字幕日韩| 国产精品久久久人人做人人爽| 亚洲欧美日韩卡通动漫| 国产av一区在线观看免费| 国内精品一区二区在线观看| 香蕉av资源在线| 精品久久久久久久末码| 亚洲人成网站高清观看| 91字幕亚洲| 琪琪午夜伦伦电影理论片6080| av女优亚洲男人天堂 | 一个人免费在线观看的高清视频| 99久久国产精品久久久| 超碰成人久久| 国产又色又爽无遮挡免费看| 久久久久九九精品影院| 看免费av毛片| 男女午夜视频在线观看| 俄罗斯特黄特色一大片| 亚洲av电影不卡..在线观看| 日韩精品中文字幕看吧| 色综合亚洲欧美另类图片| 一个人观看的视频www高清免费观看 | 超碰成人久久| 国产免费av片在线观看野外av| 久久天躁狠狠躁夜夜2o2o| 国产激情久久老熟女| 网址你懂的国产日韩在线| 久久精品aⅴ一区二区三区四区| 男女视频在线观看网站免费| 国产v大片淫在线免费观看| 老司机深夜福利视频在线观看| 亚洲精品一区av在线观看| 叶爱在线成人免费视频播放| 国产日本99.免费观看| 亚洲av成人一区二区三| 中文字幕人成人乱码亚洲影| 18禁裸乳无遮挡免费网站照片| 欧美大码av| 中文亚洲av片在线观看爽| 日韩欧美三级三区| 午夜激情福利司机影院| 又黄又粗又硬又大视频| 国产精品久久久人人做人人爽| 十八禁人妻一区二区| 中文字幕最新亚洲高清| 成年女人看的毛片在线观看| 色综合站精品国产| 99久国产av精品| 亚洲午夜精品一区,二区,三区| 亚洲片人在线观看| tocl精华| 欧美乱码精品一区二区三区| 亚洲美女视频黄频| 午夜精品久久久久久毛片777| 99国产精品99久久久久| 国产精品久久久av美女十八| 免费观看的影片在线观看| 免费在线观看日本一区| 日韩欧美精品v在线| а√天堂www在线а√下载| 亚洲avbb在线观看| 国产一区二区三区视频了| 床上黄色一级片| 久久这里只有精品19| 亚洲无线观看免费| 国产欧美日韩一区二区三| 免费观看的影片在线观看| 精品久久久久久久久久免费视频| 久久久久久人人人人人| 一夜夜www| 国产精品 欧美亚洲| 亚洲 欧美 日韩 在线 免费| 日韩三级视频一区二区三区| 黄色片一级片一级黄色片| 在线观看舔阴道视频| 日韩成人在线观看一区二区三区| 亚洲五月天丁香| 久久热在线av| 久久久国产欧美日韩av| 亚洲精品久久国产高清桃花| 麻豆一二三区av精品| 午夜影院日韩av| 国产成人精品无人区| 精品久久久久久久久久久久久| 亚洲国产高清在线一区二区三| 亚洲av免费在线观看| 国产淫片久久久久久久久 | 欧美日韩中文字幕国产精品一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 三级男女做爰猛烈吃奶摸视频| 亚洲18禁久久av| 午夜激情福利司机影院| 黄色日韩在线| 成人特级av手机在线观看| 狂野欧美激情性xxxx| 国产黄片美女视频| 国产精品爽爽va在线观看网站| 亚洲精品456在线播放app | 十八禁人妻一区二区| 国产精品野战在线观看| 国产主播在线观看一区二区| 亚洲成人久久爱视频| 国产精品亚洲美女久久久| 国产主播在线观看一区二区| 亚洲成人久久爱视频| 一个人看视频在线观看www免费 | 午夜精品在线福利| www日本在线高清视频| 欧美午夜高清在线| 美女 人体艺术 gogo| 欧美精品啪啪一区二区三区| 日韩三级视频一区二区三区| 国内久久婷婷六月综合欲色啪| 久久久久久人人人人人| 久久精品亚洲精品国产色婷小说| 丰满人妻一区二区三区视频av | 成熟少妇高潮喷水视频| 欧美日韩黄片免| 国内毛片毛片毛片毛片毛片| 人妻丰满熟妇av一区二区三区| 日本五十路高清| 久久亚洲真实| 国产熟女xx| 免费av毛片视频| 久久这里只有精品19| h日本视频在线播放| 国产一区二区激情短视频| 国产v大片淫在线免费观看| 精品国产超薄肉色丝袜足j| 精品久久久久久久毛片微露脸| www.精华液| 999久久久国产精品视频| 久久久久久大精品| 观看免费一级毛片| www日本黄色视频网| 给我免费播放毛片高清在线观看| 白带黄色成豆腐渣| 亚洲狠狠婷婷综合久久图片| 国产成人av激情在线播放| 亚洲乱码一区二区免费版| 中亚洲国语对白在线视频| 中文字幕高清在线视频| 国产精品,欧美在线| bbb黄色大片| 老熟妇乱子伦视频在线观看| a级毛片a级免费在线| 国产久久久一区二区三区| 老熟妇仑乱视频hdxx| 两个人视频免费观看高清| 淫妇啪啪啪对白视频| 亚洲成人中文字幕在线播放| 熟女人妻精品中文字幕| 女同久久另类99精品国产91| 国产一区二区三区在线臀色熟女| 精品国产乱子伦一区二区三区| 精品久久久久久,| 毛片女人毛片| 超碰成人久久| 天天添夜夜摸| 99久久精品一区二区三区| av欧美777| 亚洲成人久久爱视频| 黄色日韩在线| 亚洲欧美一区二区三区黑人| 91在线观看av| 哪里可以看免费的av片| 在线播放国产精品三级| 五月伊人婷婷丁香| 免费在线观看亚洲国产| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品中文字幕一二三四区| 窝窝影院91人妻| 久久精品aⅴ一区二区三区四区| 美女 人体艺术 gogo| 国产精品国产高清国产av| 国产成人福利小说| 国产aⅴ精品一区二区三区波| 两个人的视频大全免费| 性色avwww在线观看| 一级毛片高清免费大全| 亚洲国产看品久久| 免费观看的影片在线观看| 免费看光身美女| 哪里可以看免费的av片| 精品无人区乱码1区二区| 亚洲 欧美 日韩 在线 免费| 在线观看美女被高潮喷水网站 | 亚洲avbb在线观看| 黑人欧美特级aaaaaa片| 欧美不卡视频在线免费观看| 国内精品美女久久久久久| 99热只有精品国产| 麻豆av在线久日| 很黄的视频免费| 亚洲专区字幕在线| 日韩人妻高清精品专区| 欧美性猛交黑人性爽| av女优亚洲男人天堂 | 亚洲欧美日韩无卡精品| 国产一区二区激情短视频| 国产伦精品一区二区三区视频9 | 男插女下体视频免费在线播放| 亚洲专区字幕在线| 国产精品九九99| 国产v大片淫在线免费观看| 精品一区二区三区视频在线观看免费| 国内揄拍国产精品人妻在线| cao死你这个sao货| 1000部很黄的大片| 美女cb高潮喷水在线观看 | 亚洲精华国产精华精| 天堂影院成人在线观看| 色噜噜av男人的天堂激情| 99久久国产精品久久久| 久久久久久九九精品二区国产| 中国美女看黄片| 韩国av一区二区三区四区| 免费高清视频大片| 国产又黄又爽又无遮挡在线| 草草在线视频免费看| 老司机午夜十八禁免费视频| 亚洲成人久久性| 中国美女看黄片| 国内精品一区二区在线观看| 亚洲七黄色美女视频| 9191精品国产免费久久| 精品久久久久久久毛片微露脸| 精品国产乱子伦一区二区三区| 日韩欧美一区二区三区在线观看| 美女扒开内裤让男人捅视频| 亚洲国产看品久久| 成熟少妇高潮喷水视频| 国产一区在线观看成人免费| 国产精品亚洲av一区麻豆| 一区二区三区激情视频| 成人无遮挡网站| 熟妇人妻久久中文字幕3abv| 俄罗斯特黄特色一大片| 国产午夜精品论理片| 丰满人妻一区二区三区视频av | 一级毛片女人18水好多| 99热只有精品国产| 日本一本二区三区精品| 深夜精品福利| 亚洲av成人精品一区久久| 日韩欧美一区二区三区在线观看| 亚洲第一电影网av| 国产精品久久久久久精品电影| 色综合欧美亚洲国产小说| 日韩欧美国产一区二区入口| 久久久色成人| 亚洲成人久久性| 禁无遮挡网站| 亚洲 欧美 日韩 在线 免费| 两性夫妻黄色片| 麻豆成人av在线观看| 搡老妇女老女人老熟妇| 91九色精品人成在线观看| 精品午夜福利视频在线观看一区| a级毛片在线看网站| 国产黄色小视频在线观看| 免费观看人在逋| 99久久综合精品五月天人人| 亚洲国产精品合色在线| 男女床上黄色一级片免费看| 母亲3免费完整高清在线观看| 国产一区二区在线观看日韩 | 国产av在哪里看| 国产一级毛片七仙女欲春2| 1000部很黄的大片| 精品一区二区三区视频在线观看免费| 美女黄网站色视频| 一夜夜www| 国产精品永久免费网站| 成年女人看的毛片在线观看| 麻豆av在线久日| 啦啦啦观看免费观看视频高清| 97人妻精品一区二区三区麻豆| 精品无人区乱码1区二区| 亚洲成av人片在线播放无| 久久精品影院6| 欧美色欧美亚洲另类二区| 悠悠久久av| 亚洲中文字幕日韩| 性色av乱码一区二区三区2| 国产激情欧美一区二区| 婷婷六月久久综合丁香| 色av中文字幕| 国产爱豆传媒在线观看| 18美女黄网站色大片免费观看| 亚洲在线观看片| 一级毛片精品| 人人妻人人看人人澡| 好男人电影高清在线观看| 美女免费视频网站| 久久中文字幕人妻熟女| 国产激情久久老熟女| av黄色大香蕉| 少妇熟女aⅴ在线视频| 99久久久亚洲精品蜜臀av| 黄片小视频在线播放| 亚洲欧美日韩卡通动漫| 国产真实乱freesex| 精品久久久久久久久久久久久| 国产精品电影一区二区三区| 亚洲成av人片在线播放无| 黑人巨大精品欧美一区二区mp4| 国产一区在线观看成人免费| 国产高清视频在线观看网站| 中亚洲国语对白在线视频| 欧美3d第一页| 一区二区三区激情视频| 久久中文字幕人妻熟女| 国产亚洲精品一区二区www| 国产伦精品一区二区三区视频9 | 在线观看美女被高潮喷水网站 | 国语自产精品视频在线第100页| 午夜两性在线视频| 无遮挡黄片免费观看| 19禁男女啪啪无遮挡网站| 欧美色欧美亚洲另类二区| 很黄的视频免费| 又大又爽又粗| 国产精品永久免费网站| 国产主播在线观看一区二区| x7x7x7水蜜桃| 免费在线观看影片大全网站| 欧美成狂野欧美在线观看| 国产精品亚洲美女久久久| 午夜免费成人在线视频| 欧美高清成人免费视频www| 国产高潮美女av| 床上黄色一级片| 1024手机看黄色片| 91老司机精品| 又粗又爽又猛毛片免费看| 99在线人妻在线中文字幕| 色av中文字幕| 曰老女人黄片| 美女午夜性视频免费| 狠狠狠狠99中文字幕| 操出白浆在线播放| 成人高潮视频无遮挡免费网站| 亚洲一区二区三区色噜噜| 国产高清视频在线观看网站| 老司机午夜福利在线观看视频| 91九色精品人成在线观看| 亚洲自拍偷在线| 欧美乱妇无乱码| 999精品在线视频| 亚洲欧美一区二区三区黑人| 精品久久久久久久久久免费视频| 午夜福利18| 久久草成人影院| 欧美一区二区精品小视频在线| 99热6这里只有精品| 久久香蕉精品热| 真人一进一出gif抽搐免费| 日韩欧美 国产精品| 香蕉国产在线看| 成年版毛片免费区| 一二三四社区在线视频社区8| 99精品在免费线老司机午夜| 亚洲国产看品久久| av黄色大香蕉| 亚洲精品粉嫩美女一区| 99热这里只有是精品50| 变态另类丝袜制服| 91老司机精品| 亚洲狠狠婷婷综合久久图片| 欧美丝袜亚洲另类 | 性色av乱码一区二区三区2| 精品久久久久久久人妻蜜臀av| 国产伦精品一区二区三区视频9 | 欧美性猛交黑人性爽| 色综合婷婷激情| 88av欧美| 中文字幕av在线有码专区| 免费高清视频大片| 日本三级黄在线观看| 18禁国产床啪视频网站| 国产激情久久老熟女| 97超视频在线观看视频| 欧美中文综合在线视频| 国产精品av久久久久免费| 亚洲在线自拍视频| 亚洲美女视频黄频| 成年女人毛片免费观看观看9| 最新中文字幕久久久久 | 99久久精品一区二区三区| 1024香蕉在线观看| 高清在线国产一区| 少妇丰满av| 亚洲天堂国产精品一区在线| 久久亚洲精品不卡| 亚洲欧美激情综合另类| 日韩 欧美 亚洲 中文字幕| 美女 人体艺术 gogo| 欧美日韩瑟瑟在线播放| 十八禁人妻一区二区| 亚洲国产精品合色在线| 欧美乱妇无乱码| 熟女电影av网| 一夜夜www| 九九在线视频观看精品| 国产高清有码在线观看视频| 亚洲专区国产一区二区| 69av精品久久久久久| 欧美日韩国产亚洲二区| 免费搜索国产男女视频| 超碰成人久久| 99精品欧美一区二区三区四区| 天堂av国产一区二区熟女人妻| 99精品欧美一区二区三区四区| 成人亚洲精品av一区二区| 2021天堂中文幕一二区在线观| 啪啪无遮挡十八禁网站| 午夜视频精品福利| 一本一本综合久久| 久久精品人妻少妇| 国产亚洲av嫩草精品影院| 国产精品乱码一区二三区的特点| 亚洲美女视频黄频| 日本三级黄在线观看| 精华霜和精华液先用哪个| 日韩欧美国产在线观看| 黄色女人牲交| 国产精品电影一区二区三区| 欧美午夜高清在线| 日日夜夜操网爽| 久久中文看片网| 久久伊人香网站| 女生性感内裤真人,穿戴方法视频| 精品久久久久久久久久免费视频| 极品教师在线免费播放| 少妇的逼水好多| 日韩欧美三级三区| 国产人伦9x9x在线观看| 18禁黄网站禁片免费观看直播| 亚洲天堂国产精品一区在线| 精华霜和精华液先用哪个| 国产精品一区二区免费欧美| 色吧在线观看| 啦啦啦观看免费观看视频高清| АⅤ资源中文在线天堂| 伦理电影免费视频| 欧美一区二区国产精品久久精品| 久久久精品大字幕| av在线天堂中文字幕| 真人做人爱边吃奶动态| 亚洲七黄色美女视频| 国产亚洲精品av在线| 亚洲精品美女久久av网站| 欧美乱色亚洲激情| 久久久久久大精品| 午夜精品一区二区三区免费看| 久久天堂一区二区三区四区| 熟女人妻精品中文字幕| 一级a爱片免费观看的视频| 成年女人毛片免费观看观看9| 88av欧美| 伦理电影免费视频| 男女午夜视频在线观看| 一本精品99久久精品77| 免费电影在线观看免费观看| 亚洲av片天天在线观看| 国产一区二区在线av高清观看| 欧美高清成人免费视频www| 99国产综合亚洲精品| 免费在线观看成人毛片| h日本视频在线播放| 久久久久久久久久黄片| 中文字幕精品亚洲无线码一区| 在线国产一区二区在线| 精品久久久久久,| 国产精品 国内视频| 国内久久婷婷六月综合欲色啪| 久久精品人妻少妇| 美女 人体艺术 gogo| 久久国产精品影院| 日韩欧美三级三区| 日本撒尿小便嘘嘘汇集6| 变态另类成人亚洲欧美熟女| 99热只有精品国产| 波多野结衣巨乳人妻| 视频区欧美日本亚洲| 中文字幕久久专区| 亚洲精品美女久久av网站| 十八禁人妻一区二区| 黑人操中国人逼视频| 舔av片在线| 成人18禁在线播放| 五月伊人婷婷丁香| 色播亚洲综合网| 在线观看66精品国产| 久久久久久久午夜电影| 天堂影院成人在线观看| 成人国产一区最新在线观看| 麻豆国产av国片精品| h日本视频在线播放| 免费看美女性在线毛片视频| 一本综合久久免费| 久久精品影院6| svipshipincom国产片| 村上凉子中文字幕在线| 国产日本99.免费观看| 九九在线视频观看精品| 搡老妇女老女人老熟妇| 一卡2卡三卡四卡精品乱码亚洲| 欧美在线黄色| 亚洲欧美日韩高清专用| 亚洲午夜精品一区,二区,三区| 亚洲av日韩精品久久久久久密| 色综合婷婷激情| 国产精品免费一区二区三区在线| 久久草成人影院| 老熟妇仑乱视频hdxx| 欧美中文综合在线视频| 在线永久观看黄色视频| 亚洲专区字幕在线| 黄频高清免费视频| 人妻夜夜爽99麻豆av| 青草久久国产| 手机成人av网站| 三级毛片av免费| 日韩精品青青久久久久久| 99久久成人亚洲精品观看| 岛国视频午夜一区免费看| 国产毛片a区久久久久| 欧美在线一区亚洲| 国产av在哪里看| 最新在线观看一区二区三区| 久久精品国产清高在天天线| 国产成人av激情在线播放| 美女黄网站色视频| 国产精品国产高清国产av| 色播亚洲综合网| 2021天堂中文幕一二区在线观| 色视频www国产| 亚洲熟妇中文字幕五十中出| 人人妻人人看人人澡| 青草久久国产| www.自偷自拍.com| 最近最新中文字幕大全电影3| 免费高清视频大片|