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

    亞熱帶丘陵區(qū)池塘底棲動(dòng)物群落特征與水體面源 污染關(guān)系研究

    2022-06-29 09:31:58彭健李希李裕元王浩孟岑曾睿王棟
    關(guān)鍵詞:生物量池塘物種

    彭健 ,李希 ,李裕元 ,王浩 ,孟岑,曾睿,王棟

    (1. 中國科學(xué)院亞熱帶農(nóng)業(yè)生態(tài)研究所,亞熱帶農(nóng)業(yè)生態(tài)過程重點(diǎn)實(shí)驗(yàn)室,湖南 長沙 410125;2. 中國科學(xué)院大學(xué),北京 100049;3. 湖南艾布魯環(huán)保科技股份有限公司, 湖南 長沙 410000)

    隨著點(diǎn)源污染不斷得到有效控制,農(nóng)業(yè)面源污染問題已成為我國農(nóng)村環(huán)境的最突出問題[1]。面源污染量大面廣,污染源分散度高,對(duì)水體污染持續(xù)性較強(qiáng)[2],容易導(dǎo)致底棲動(dòng)物功能、豐富度和群落結(jié)構(gòu)的改變,因此,水體健康和底棲動(dòng)物的相關(guān)研究備受國內(nèi)外學(xué)者關(guān)注[3],探討池塘底棲動(dòng)物群落特征與水體面源污染的關(guān)系對(duì)于池塘水體面源污染的預(yù)防與修復(fù)治理具有一定的理論價(jià)值和現(xiàn)實(shí)意義。

    目前我國農(nóng)村地區(qū)以污染物遷移方向?yàn)橹骶€、以小流域或集水區(qū)為單元開展綜合生態(tài)治理是治理面源污染的重要方向[4-5]。池塘是我國南方地區(qū)十分常見的小型水體[6-7],作為小流域水文循環(huán)的重要節(jié)點(diǎn),對(duì)于調(diào)蓄地表經(jīng)流、抗旱和農(nóng)業(yè)面源污染物消納轉(zhuǎn)化均具有重要作用[8]。同時(shí),池塘水體水質(zhì)的變化也會(huì)引起包含底棲動(dòng)物在內(nèi)的水生生態(tài)系統(tǒng)特征發(fā)生不同程度的變化,從而進(jìn)一步影響到池塘的生態(tài)功能。

    底棲動(dòng)物作為池塘等濕地生態(tài)系統(tǒng)中物質(zhì)循環(huán)和能量流動(dòng)的重要參與者,是池塘濕地生態(tài)功能健康的重要標(biāo)志之一[9],在水生生態(tài)系統(tǒng)中扮演消費(fèi)者角色,一方面以細(xì)菌、有機(jī)碎屑顆粒和底棲藻類為食,促進(jìn)水體有機(jī)物降解,另一方面可為魚類提供天然餌料,加速池塘生態(tài)系統(tǒng)營養(yǎng)物質(zhì)的遷移轉(zhuǎn)化[10]。由于底棲動(dòng)物對(duì)環(huán)境變化反應(yīng)敏感,當(dāng)水體受到污染時(shí),其群落結(jié)構(gòu)及多樣性會(huì)發(fā)生相應(yīng)的變化,可對(duì)濕地系統(tǒng)水質(zhì)變化起到指示作用[11],因此也受到水環(huán)境研究者的廣泛關(guān)注。

    目前針對(duì)河流水系[12]、湖泊[13]、水庫[14]等大型水體底棲動(dòng)物的研究常見報(bào)道,與農(nóng)業(yè)面源污染治理密切相關(guān)的小型水體中底棲動(dòng)物的研究則相對(duì)集中于人工濕地[15],人工濕地與天然池塘之間的底棲動(dòng)物群落特征既有區(qū)別又有聯(lián)系,進(jìn)一步研究天然池塘底棲動(dòng)物群落特征可以反映一個(gè)地區(qū)池塘水體面源污染狀況,建立具有一定指示作用的生物指標(biāo)體系。為此,本研究于2020 年7 月(夏季)和12月(冬季)對(duì)亞熱帶丘陵區(qū)金井小流域15 口典型池塘底棲動(dòng)物物種、豐度、生物量和水體污染情況進(jìn)行了現(xiàn)場采樣調(diào)查,通過分析兩季池塘水體污染情況和相關(guān)底棲動(dòng)物的群落特征,探討了底棲動(dòng)物群落特征隨季節(jié)和水體污染程度變化的關(guān)系,研究結(jié)果可為池塘水體面源污染的治理修復(fù)提供參考依據(jù)和數(shù)據(jù)支撐。

    1 材料與方法

    1.1 研究區(qū)概況與采樣位點(diǎn)設(shè)置

    研究區(qū)位于湖南省長沙縣金井鎮(zhèn),地理坐標(biāo)范圍28°30′~28°39′ N、113°18′~113°29′ E(圖1),區(qū)內(nèi)平均海拔98.3 m,為典型丘陵地貌,屬于典型亞熱帶濕潤季風(fēng)氣候,多年平均氣溫17.1 ℃,全年無霜期274 d,日照時(shí)數(shù)1 663 h,熱量充足,多年平均降水量1 394.6 mm,主要集中于4—7 月份。當(dāng)?shù)赝恋乩梅绞街饕修r(nóng)田(36.47%)、茶園(3.40%)、林地(49.01%)等[16]。利用研究區(qū)高分辨率影像(2016年12 月,0.5 m 分辨率),采用ArcGIS 軟件完成研究區(qū)內(nèi)森林、農(nóng)田、居民區(qū)、池塘、水庫和道路等土地利用要素的目視解譯與矢量圖層輸出,按影像解譯結(jié)果,該小流域池塘數(shù)量有2 012 個(gè),密度高達(dá)19 個(gè)/km2。課題組定期監(jiān)測其中142 個(gè)池塘的水質(zhì)變化,通過實(shí)地調(diào)查,獲取池塘周邊景觀類型、池塘經(jīng)營管理(是否養(yǎng)魚)和池塘周邊污染物輸入來源、綜合水質(zhì)情況等基本信息,從142 個(gè)池塘中挑選15 個(gè)水質(zhì)狀況、污染物來源、分布區(qū)域和經(jīng)營模式等有顯著差異的典型池塘(平均水深約3 m)進(jìn)行底棲動(dòng)物群落調(diào)查與采樣分析,采樣點(diǎn)分布見圖1。

    1.2 池塘底棲動(dòng)物調(diào)查

    分別于2020 年7 月(夏季)和12 月(冬季)開展兩次底棲動(dòng)物的調(diào)查采樣分析,其中底棲動(dòng)物調(diào)查采用1/16 m2彼得森采泥器,每口池塘采用五點(diǎn)取樣法采集底泥樣,泥樣當(dāng)即用60 目尼龍網(wǎng)篩清洗后倒入采樣瓶,瓶中加入10%甲醛密封保存后帶回實(shí)驗(yàn)室,再用60 目尼龍網(wǎng)篩在自來水下沖洗,將甲醛洗凈后將網(wǎng)篩內(nèi)的剩余物倒入解剖盤中,用鑷子將底棲動(dòng)物樣本逐一撿出,放入盛有80%乙醇溶液的離心管中,隨后在光學(xué)顯微鏡下逐個(gè)進(jìn)行物種鑒定和分類計(jì)數(shù),并用濾紙將底棲動(dòng)物表面水分吸干后在萬分之一天平上稱重,最終匯總記錄結(jié)果并換算成單位面積的個(gè)體數(shù)量即豐度(ind/m2)和鮮重質(zhì)量即生物量(g/m2)[17]。

    1.3 水樣采集與測定

    在進(jìn)行底棲動(dòng)物調(diào)查的同時(shí)采集水樣并帶回室內(nèi)進(jìn)行相關(guān)水質(zhì)指標(biāo)的分析,部分水體理化指標(biāo)現(xiàn)場實(shí)時(shí)測定。用哈希便攜式水質(zhì)分析儀(型號(hào):HQ40d)測定水溫(T)、酸堿度(pH)和溶解氧(DO),用雷磁DDB-303A 便攜式電導(dǎo)率儀測定電導(dǎo)率(EC)。

    水質(zhì)指標(biāo)測定方法按照《水和廢水監(jiān)測分析方法》(第4 版)進(jìn)行。氨氮(NH4+-N)、硝態(tài)氮 (NO3--N)和溶解態(tài)有機(jī)碳(DOC)的測定方法為將水樣抽濾過0.45 μm 膜后,取濾液直接上流動(dòng)注射儀(AA3,德國SEAL 公司)和TOC 分析儀(島津Vwp,日本)測定;總氮(TN)采用堿性過硫酸鉀消解后用流動(dòng)注射儀(AA3,德國SEAL 公司)測定;總磷(TP)采用過硫酸鉀消解后,用鉬酸銨分光光度法測定[18]。

    1.4 數(shù)據(jù)處理與計(jì)算

    底棲動(dòng)物群落特征指數(shù)采用優(yōu)勢(shì)度指數(shù)(Y)確定池塘底棲動(dòng)物的優(yōu)勢(shì)物種;采用單位面積底棲動(dòng)物的數(shù)量(豐度,ind/m2)和生物量(g/m2)來評(píng)價(jià)底棲動(dòng)物的數(shù)量特征;采用Shannon-Weiner 多樣性指數(shù)(H ′)、Margalef 豐富度指數(shù)(D)、Pielou均勻度指數(shù)(J)分析底棲動(dòng)物的多樣性特征。各主要指數(shù)的計(jì)算公式為:

    式中:Y 為優(yōu)勢(shì)度指數(shù),當(dāng)Y >0.02 時(shí),確定該物種為優(yōu)勢(shì)種;Ai為相對(duì)豐度,即種i 的個(gè)體數(shù)量占總物種個(gè)體數(shù)量的比例;Fi為物種i 出現(xiàn)的頻率;S為總的物種數(shù)目。

    水質(zhì)健康評(píng)價(jià)參照有關(guān)監(jiān)測技術(shù)規(guī)范對(duì)相關(guān)水質(zhì)參數(shù)進(jìn)行限定,本研究中pH 的限定范圍取6.50~8.50,處于此區(qū)間內(nèi)為合格;DO 限定值取5 mg/L,高于5 mg/L 為合格;NH4+-N、NO3--N、TN、TP 和DOC 的限定值分別取1.00 mg/L、1.00 mg/L、1.00 mg/L、0.20 mg/L 和6.00 mg/L。評(píng)價(jià)方法依據(jù)農(nóng)用水源環(huán)境質(zhì)量監(jiān)測技術(shù)規(guī)范[19],采用單項(xiàng)污染指數(shù)和負(fù)荷比對(duì)監(jiān)測參數(shù)進(jìn)行單項(xiàng)評(píng)價(jià),再通過綜合污染指數(shù)對(duì)水體環(huán)境質(zhì)量進(jìn)行整體評(píng)價(jià)。單項(xiàng)污染指數(shù)計(jì)算公式為:

    式中:Pi為水環(huán)境中污染物i 的污染指數(shù);Ci為水環(huán)境中污染物i 的實(shí)測值;C0為水環(huán)境中污染物的限量標(biāo)準(zhǔn)值。根據(jù)該計(jì)算方法,為了便于累積計(jì)算比較,需對(duì)Pi數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化計(jì)算:當(dāng)Pi≤1 時(shí)表示水環(huán)境未受污染,指標(biāo)合格,標(biāo)準(zhǔn)化值Ni=Pi;當(dāng)Pi>1 時(shí),表示水環(huán)境受到污染,指標(biāo)不合格,標(biāo)準(zhǔn)化值Ni= 1.0+5×lgPi,由此,在單項(xiàng)污染指數(shù)評(píng)價(jià)的基礎(chǔ)上,采用兼顧單項(xiàng)污染指數(shù)最大值和平均值的綜合污染指數(shù)Pj進(jìn)行評(píng)價(jià),其計(jì)算公式為:

    式中:Nmaxi為標(biāo)準(zhǔn)化的最大單項(xiàng)污染指數(shù),Navei為標(biāo)準(zhǔn)化的平均單項(xiàng)污染指數(shù)。依據(jù)水體環(huán)境綜合污染指數(shù),可將池塘水質(zhì)狀況分為5 個(gè)等級(jí)(表1),對(duì)所研究池塘水質(zhì)的污染程度進(jìn)行評(píng)價(jià)。

    表1 水體綜合污染指數(shù)分級(jí)Table 1 Classification of comprehensive pollution index of waters

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

    運(yùn)用單因素方差分析(One-way ANOVA)來檢驗(yàn)不同類型池塘環(huán)境因子、底棲動(dòng)物群落特征差異的顯著性,采用Duncan’s 多重比較檢驗(yàn)組間差異。采用皮爾遜相關(guān)分析來檢驗(yàn)底棲動(dòng)物種類數(shù)、不同類群密度和生物量與主要環(huán)境因子的相關(guān)性,如果P <0.05,則認(rèn)為有顯著相關(guān)關(guān)系,并對(duì)底棲動(dòng)物與環(huán)境因子間進(jìn)行主成分分析(PCA)。以上所有分析在SPSS 21.0 軟件中完成。使用CANOCO 5 進(jìn)行冗余分析(RDA)。

    2 結(jié)果與分析

    2.1 池塘水體污染程度與分類

    對(duì)15 個(gè)調(diào)查池塘兩個(gè)觀測季節(jié)綜合污染指數(shù)的計(jì)算結(jié)果表明,水體污染均達(dá)到輕污染程度以上(Pj>1,圖2),研究區(qū)池塘水體主要污染物平均值為3.03,超出警戒水平??傮w上,池塘冬季水體污染程度明顯高于夏季,4、6 和13~15 號(hào)池塘冬季污染程度有所降低,區(qū)域內(nèi)夏季高溫多雨,冬季低溫少雨,降雨可能是導(dǎo)致兩季污染程度不同的主要原因,13~15 號(hào)池塘可能由于周邊地理環(huán)境屬于林地而影響了地表徑流,導(dǎo)致夏季污染程度反而高于冬季,而4 號(hào)和6 號(hào)池塘周邊為居民地,受人為因素影響較大。根據(jù)實(shí)際調(diào)查結(jié)果,結(jié)合水體污染指數(shù)分級(jí),可將15 個(gè)池塘劃分為3 類,即輕度污染(綜合污染指數(shù)1~2 級(jí)):夏季7~12 號(hào)塘,冬季7 號(hào)塘和13 號(hào)塘;中度污染(綜合污染指數(shù)3 級(jí)):夏季1~3 號(hào)池塘,冬季11、14 和15 號(hào)塘;重度污染(綜合污染指數(shù)>3 級(jí)):夏季4~6 號(hào)塘和13~15 號(hào)塘,冬季1~6、8、9、10 和12 號(hào)塘。

    圖2 水環(huán)境綜合污染指數(shù)Fig. 2 Comprehensive pollution index of water environment

    根據(jù)池塘分類結(jié)果,參照國家地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB 3838—2002),金井河小流域池塘水體總磷測定結(jié)果一般可達(dá)到地表水Ⅱ類~Ⅲ類之間,兩季平均總磷含量僅為0.09 mg/L,而總氮一般為Ⅴ類或劣Ⅴ類,夏季總氮平均含量為3.98 mg/L,地表水Ⅳ類超標(biāo)率達(dá)到73.3%,冬季總氮和氨氮平均含量甚至高達(dá)6.21 mg/L 和4.59 mg/L,總氮和氨氮地表水Ⅴ類超標(biāo)率分別為100%和86%以上(表2),表明 總氮和氨氮是金井河流域池塘水體污染的主要指標(biāo)。

    表2 不同污染程度池塘氮磷指標(biāo)Table 2 Nitrogen and phosphorus indexes of ponds under different pollution degrees

    2.2 底棲動(dòng)物群落特征

    對(duì)底棲動(dòng)物種類的調(diào)查結(jié)果顯示,2 個(gè)季節(jié)共采集底棲動(dòng)物22 種,包括節(jié)肢動(dòng)物、環(huán)節(jié)動(dòng)物和軟體動(dòng)物三大類群,節(jié)肢動(dòng)物種類14 種,占群落總種數(shù)64%,軟體動(dòng)物和環(huán)節(jié)動(dòng)物種數(shù)均為4 種,群落結(jié)構(gòu)相對(duì)較為簡單(表3)。兩季物種數(shù)略有差異,其中夏季略多為16 種,冬季略少為14 種??傮w而言,研究區(qū)池塘中主要以節(jié)肢動(dòng)物和環(huán)節(jié)動(dòng)物為主,軟體動(dòng)物數(shù)目和種類均較少,耐污性較強(qiáng)的霍甫水絲蚓、黃色羽搖蚊和中國長足搖蚊在研究區(qū)域池塘中占據(jù)絕對(duì)優(yōu)勢(shì)。

    底棲動(dòng)物豐度調(diào)查結(jié)果顯示,夏季和冬季污染程度較低的輕污染池塘均具有較高的底棲動(dòng)物豐度,平均分別為930 ind/m2和1 808 ind/m2(圖3), 隨污染程度的升高底棲動(dòng)物豐度呈顯著下降的趨勢(shì)(P<0.05),重污染池塘豐度平均僅為251 ind/m2(夏季)和352 ind/m2(冬季)??傮w而言,冬季平均豐度(913 ind/m2)顯著高于夏季(544 ind/m2)。底棲動(dòng)物生物量的變化(圖3)與豐度表現(xiàn)為相同的趨勢(shì),冬季平均生物量(6.67 g/m2)總體上顯著高于夏季(1.48 g/m2,P <0.05),且隨污染程度增加底棲動(dòng)物生物量呈顯著增加趨勢(shì)(P <0.05),但冬季中污染和重污染之間生物量差異不顯著(P >0.05)。

    表3 不同池塘冬夏兩季底棲動(dòng)物物種及優(yōu)勢(shì)度(Y)Table 3 Species and dominance (Y) of zoobenthos in different ponds in winter and summer

    圖3 池塘不同污染程度下底棲動(dòng)物豐度和生物量Fig. 3 Abundance and biomass of zoobenthos under different pollution levels in ponds

    圖4 池塘不同污染程度下底棲動(dòng)物的相對(duì)豐度Fig. 4 Relative abundance of zoobenthos under different pollution levels in ponds

    從池塘不同物種底棲動(dòng)物相對(duì)豐度的變化來看(圖4),其受季節(jié)影響較大,冬季黃色羽搖蚊和幽蚊相對(duì)豐度明顯要低于夏季,而霍甫水絲蚓冬季的相對(duì)豐度則顯著增加(P<0.01),庫蠓、蘇氏尾鰓蚓和克拉伯水絲蚓 的相對(duì)豐度在冬季也有一定程度的增加。底棲動(dòng)物相對(duì)豐度受污染程度的影響也較大,總體來看,黃色羽搖蚊在中度污染池塘中相對(duì)豐度最大,水體污染程度升高會(huì)導(dǎo)致其相對(duì)豐度的下降。同樣,霍甫水絲蚓在總體污染程度較低的夏季,其相對(duì)豐度與污染程度成正比,而在總體污染程度較高的冬季,其相對(duì)豐度則與污染程度成反比。蘇氏尾鰓蚓則在冬季重污染程度池塘中相對(duì)豐度顯著增加(P<0.01)。

    池塘底棲動(dòng)物多樣性指數(shù)的分析結(jié)果表明,總體上冬季各項(xiàng)指數(shù)均顯著高于夏季(P<0.05,表4),而多樣性指數(shù)隨水質(zhì)變化的趨勢(shì)基本一致:Shannon- Weiner 多樣性指數(shù)(H′)和Margalef 豐富度指數(shù)(D)均隨著池塘水體污染程度的升高呈降低趨勢(shì),而Pielou 均勻度指數(shù)(J)變化趨勢(shì)不明顯,且三種多樣性指數(shù)中除了冬季的物種豐富度指數(shù)隨水質(zhì)變化的差異達(dá)到統(tǒng)計(jì)顯著水平以外(P<0.05),其余指標(biāo)的差異均不顯著(P>0.05)。

    表4 不同污染程度池塘底棲動(dòng)物多樣性指數(shù)Table 4 Diversity index of zoobenthos under different pollution levels in ponds

    2.3 底棲動(dòng)物物種分布與環(huán)境因子的相關(guān)性分析

    對(duì)池塘水環(huán)境因子的主成分分析(PCA)結(jié)果表明,9 個(gè)環(huán)境因子的累計(jì)解釋方差占總方差的75.94%,其中第一(F1)、第二(F2)和第三主成分 (F3)的貢獻(xiàn)率分別為36.19%、21.95%和17.81% (表5),表明9 種環(huán)境因子均對(duì)底棲動(dòng)物的組成和分布有較為重要的影響。F1 主要為TN、NH4+-N、TP 的組合,且與三者均呈正荷載;F2 主要為T、pH、DO、NO3--N 組合,與T、DO 呈正荷載,與pH、NO3

    --N 呈負(fù)荷載;F3 主要為EC 和DOC,均呈正荷載。

    表5 環(huán)境因子總方差解釋表Table 5 Total variance explained of environmental factors

    通過前選法和蒙特卡羅檢驗(yàn)排除貢獻(xiàn)小的因子,進(jìn)一步對(duì)6 種優(yōu)勢(shì)底棲物種(Y>0.02,且出現(xiàn)頻度超過一次的)霍甫水絲蚓、蘇氏尾鰓蚓、黃色羽搖蚊、中國長足搖蚊、幽蚊、庫蠓的生物量與環(huán)境因子間進(jìn)行RDA 分析,結(jié)果發(fā)現(xiàn),總解釋變量達(dá)到41.39%,其中軸1 和軸2 的解釋率較高,特征值分別為0.150 和0.132,分別解釋了14.99%和13.15%的物種數(shù)據(jù)方差變異以及34.68%和30.43%的物種與環(huán)境關(guān)系變異(表6),表明水體環(huán)境對(duì)底棲動(dòng)物的豐度、生物量以及優(yōu)勢(shì)度有重要影響。從RDA 排序圖中可知:第一排序軸與N 相關(guān)環(huán)境因子呈顯著正相關(guān),與DO、T 呈顯著負(fù)相關(guān)。黃色羽搖蚊與T 呈顯著正相關(guān)關(guān)系,而蘇氏尾鰓蚓和霍甫水絲蚓這兩種耐污能力強(qiáng)的物種則與TN、NH4+-N 呈顯著正相關(guān)。第一排序軸上顯示最大正值的主要有霍甫水絲蚓、蘇氏尾鰓蚓、中國長足搖蚊和幽蚊等,表明這些物種能夠耐受較高程度污染(圖5)。

    3 討論

    3.1 典型池塘底棲動(dòng)物群落特征及其季節(jié)變化

    本研究發(fā)現(xiàn)金井小流域池塘兩季底棲動(dòng)物物種組成相對(duì)簡單,物種多樣性總體上較低,冬季豐度與生物量顯著高于夏季。區(qū)內(nèi)調(diào)查已知的底棲動(dòng)物僅有22 種,物種數(shù)目和多樣性均遠(yuǎn)低于我國相關(guān)區(qū)域河流、湖泊等大型水體的調(diào)查結(jié)果,如韋建福等[20]在廣西南溪河研究發(fā)現(xiàn)的大型底棲動(dòng)物有208種,王丑明等[21]在洞庭湖的研究共鑒定出大型底棲動(dòng)物58 種,導(dǎo)致這一現(xiàn)象的主要原因可能在于池塘屬于相對(duì)破碎化的微小型生境,相互間具有一定的獨(dú)立性,而且池塘與周邊大型水體之間的物種遷移途徑和信息交流機(jī)會(huì)相對(duì)較少[22]。本研究中發(fā)現(xiàn)的軟體動(dòng)物也相對(duì)較少,僅有中華圓田螺、銅銹環(huán)棱螺、膀胱螺、凸旋螺和小土蝸5 種,且出現(xiàn)頻次較低,均不是優(yōu)勢(shì)種,這與大型湖泊內(nèi)底棲動(dòng)物的物種構(gòu)成有顯著區(qū)別,如鄒亮華等[23]對(duì)鄱陽湖的研究結(jié)果表明,河蜆、銅銹環(huán)棱螺和大沼螺在各湖區(qū)均占據(jù)顯著優(yōu)勢(shì)地位。研究區(qū)軟體動(dòng)物偏少的原因是多方面的,除了上述原因以外,可能還與池塘水體NH4+-N 濃度總體偏高(0.3~9.0 mg/L,表2)有密切關(guān)系,據(jù)Ilarri 等[24]的研究,NH4+-N 對(duì)底棲動(dòng)物尤其是軟體動(dòng)物會(huì)產(chǎn)生較高的生物毒性,會(huì)導(dǎo)致軟體動(dòng)物的豐度和生物量下降。

    表6 底棲動(dòng)物群落與環(huán)境因子冗余分析(RDA)結(jié)果Table 6 Results of redundancy analysis (RDA) of zoobenthos community and environmental factors

    圖5 底棲動(dòng)物群落與環(huán)境因子冗余分析Fig. 5 Redundancy analysis (RDA) of zoobenthos community and environmental factors

    金井河流域池塘水體冬季的污染程度總體上明顯高于夏季。研究區(qū)屬于典型亞熱帶濕潤季風(fēng)氣候,夏季高溫多雨,冬季低溫少雨,降雨對(duì)稀釋區(qū)域內(nèi)池塘水體的污染物濃度有重要作用。申雅莉等[25]對(duì)金井河小流域調(diào)查發(fā)現(xiàn),該區(qū)域池塘TN、TP 濃度呈明顯的季節(jié)性變化趨勢(shì),其中雨季較低而旱季較高,這與本研究的觀測結(jié)果基本一致,13~15 號(hào)池塘冬季綜合污染指數(shù)低于夏季,可能是因?yàn)檫@三個(gè)池塘周邊環(huán)境為林地,與農(nóng)田相比,林地能夠更好地控制徑流養(yǎng)分流失,緩解地下水污染,有利于農(nóng)業(yè)面源污染的控制[26]。相應(yīng)地,研究區(qū)冬季底棲動(dòng)物豐度和密度均顯著高于夏季(圖3),其中分布最為廣泛的霍甫水絲蚓和蘇氏尾鰓蚓等耐污種的豐度及生物量均與水體富營養(yǎng)化水平呈顯著正相關(guān)關(guān)系(圖5),對(duì)反映池塘水體的富營養(yǎng)化程度具有明顯的指示作用,其豐度及生物量越大,則水體的富營養(yǎng)化程度就越高,這與國內(nèi)相關(guān)區(qū)域的研究結(jié)果較一致,如超富營養(yǎng)化水體中霍甫水絲蚓豐度最高可達(dá)到10 524 ind/m2[27],而在中、低富營養(yǎng)化水體中僅為27 ind/m2[28],這些均表明面源污染導(dǎo)致的池塘水體富營養(yǎng)化變化對(duì)底棲動(dòng)物物種及生物量的影響程度甚至要高于溫度(季節(jié))變化的影響。水溫主要直接影響搖蚊類、寡毛類水生昆蟲的生長、繁殖和羽化進(jìn)程[29],導(dǎo)致其在夏季更為活躍,冬季則以幼蟲形態(tài)棲息于底泥之中,因此該類底棲動(dòng)物在冬季的豐度相對(duì)較高,這可能主要決定于該類物種生活史與生活習(xí)性的本性。但是郭寧寧等[30]對(duì)亞熱帶丘陵區(qū)淺水人工濕地(水深20 cm 左右)的研究表明,夏、秋季濕地底棲動(dòng)物的多樣性指數(shù)、豐富度指數(shù)和均勻度指數(shù)??傮w上均高于春、冬季,與本文結(jié)果不盡一致,表明水深可能也是影響底棲動(dòng)物多樣性的重要因素,具體影響機(jī)理尚值得進(jìn)一步深入研究。

    3.2 農(nóng)業(yè)面源污染對(duì)池塘底棲動(dòng)物群落特征的影響

    隨污染程度升高,底棲動(dòng)物豐度呈顯著下降趨勢(shì)(P <0.05),多樣性指數(shù)和豐富度指數(shù)均隨著池塘水體污染程度的升高呈降低趨勢(shì)但并未達(dá)到統(tǒng)計(jì)顯著水平(P >0.05),主要面源污染指標(biāo)TN、NH4

    +-N 可能是影響研究區(qū)池塘底棲動(dòng)物群落特征的最主要因素。李麗娟等[31]對(duì)太子河大型底棲動(dòng)物的研究發(fā)現(xiàn),氮磷等無機(jī)污染物進(jìn)入河流后會(huì)導(dǎo)致底棲動(dòng)物攝食功能群以收集者功能群和濾食者功能群為主,表明底棲動(dòng)物多樣性受水體污染影響且與污染程度呈負(fù)相關(guān)。

    對(duì)池塘底棲動(dòng)物豐度、生物量、多樣性指數(shù)和底棲動(dòng)物物種分布與池塘水質(zhì)的相關(guān)性分析結(jié)果均顯示,非季節(jié)性分布的霍甫水絲蚓和蘇氏尾鰓蚓這兩種耐污能力較強(qiáng)的物種的豐度和生物量與TN、NH4

    +-N 呈顯著正相關(guān)關(guān)系,表明水體總氮與氨氮濃度是影響本區(qū)域?qū)Φ讞珓?dòng)物物種多樣性和群落構(gòu)成的主導(dǎo)因子,霍甫水絲蚓和蘇氏尾鰓蚓對(duì)反映池塘水體的富營養(yǎng)化程度具有一定的指示作用。對(duì)于不同的底棲動(dòng)物,對(duì)環(huán)境因素變化的響應(yīng)特征也不盡相同,如郭寧寧等[30]研究發(fā)現(xiàn),水體溶氧(DO)是影響腹足綱底棲動(dòng)物(如螺類)分布的主要因素,而寡毛綱(如霍甫水絲蚓)和昆蟲綱(如搖蚊幼蟲)與水體TN、COD 的關(guān)系更為密切。

    本研究也發(fā)現(xiàn),底棲動(dòng)物生物量冬季較夏季高,但在水體污染程度由中到高變化時(shí)池塘底棲動(dòng)物生物量的增加則不顯著,這可能是因?yàn)楫?dāng)水質(zhì)指標(biāo)超過耐污種的正響應(yīng)閾值時(shí),部分耐污種達(dá)到了耐受極限而開始減少,如程佩瑄等[32]研究發(fā)現(xiàn),大型底棲動(dòng)物群落結(jié)構(gòu)對(duì)水質(zhì)指示的指標(biāo)閾值也有相應(yīng)的變化,不同底棲動(dòng)物均有其適宜的水質(zhì)范圍,超出該范圍,其存活能力就會(huì)顯著降低。富營養(yǎng)化水體雖然能為底棲動(dòng)物提供更為豐富的食物來源,但過高的富營養(yǎng)化程度也會(huì)導(dǎo)致底泥中底棲動(dòng)物群落結(jié)構(gòu)的變化,本研究中主要體現(xiàn)為物種多樣性的顯著降低,同時(shí)PCA 分析中發(fā)現(xiàn)底棲動(dòng)物豐度和生物量與pH 呈負(fù)荷載,這與蘆康樂[33]對(duì)黃河三角洲蘆葦濕地的調(diào)查結(jié)果較為一致,水體的pH 是影響底棲動(dòng)物分布的一個(gè)重要環(huán)境因子,不同的底棲動(dòng)物,其適宜的pH 不同,因而分布也不同。

    3.3 底棲動(dòng)物對(duì)農(nóng)業(yè)面源污染治理的重要意義

    本次調(diào)查的池塘主要分布于土地利用方式為農(nóng)田和茶園的區(qū)域,化肥的施用和降雨可造成農(nóng)業(yè)面源污染,導(dǎo)致池塘水體營養(yǎng)鹽濃度總體偏高。區(qū)域內(nèi)池塘底棲動(dòng)物物種多樣性較低且優(yōu)勢(shì)種均為耐污種,表明了池塘水環(huán)境較差的污染現(xiàn)狀,為防止池塘生態(tài)環(huán)境進(jìn)一步惡化,應(yīng)加強(qiáng)該區(qū)域水體環(huán)境的保護(hù)和污染治理工作。

    底棲動(dòng)物群落特征與水體污染密切關(guān)聯(lián),其作為水生生態(tài)系統(tǒng)食物鏈的重要組成部分,對(duì)營養(yǎng)鹽的流動(dòng)循環(huán)和污染物的去除有重要作用[34]。實(shí)際調(diào)查過程中也發(fā)現(xiàn),污染程度相對(duì)較高的池塘水生植物相對(duì)較少,因此底棲動(dòng)物的棲息環(huán)境也相對(duì)較差。靳聰聰?shù)萚35]發(fā)現(xiàn)沉水植物的構(gòu)建可以有效控制農(nóng)業(yè)地區(qū)水體中農(nóng)藥等污染物的濃度,降低其生態(tài)風(fēng)險(xiǎn),有利于恢復(fù)底棲動(dòng)物的多樣性。底棲動(dòng)物通過移動(dòng)、攝食和筑穴等生物擾動(dòng)作用改變沉積物的物理和化學(xué)性質(zhì),從而影響水生態(tài)系統(tǒng)中的物質(zhì)循環(huán)和能量流動(dòng),田勝艷等[36]研究發(fā)現(xiàn),生物擾動(dòng)可以促進(jìn)沉積物顆粒的遷移混合、釋放沉積物中有機(jī)物質(zhì)、改變沉積環(huán)境并提高沉積物中有機(jī)污染物的生物可利用性,從而提高水體有機(jī)污染物的去除效果,因此通過底棲動(dòng)物添加來強(qiáng)化污染池塘濕地等水體底泥污染的方法也不失為一條有效途徑,一些試驗(yàn)已經(jīng)證明了該技術(shù)的可行性,如陳桐等[37]在野外進(jìn)行的圍隔試驗(yàn)發(fā)現(xiàn),投加螺和蚌的密度為150 g/m2時(shí)對(duì)水體TN 的去除率可達(dá)72.6%,朱明璇等[38]利用水生動(dòng)植物聯(lián)合技術(shù)凈化水庫水試驗(yàn)的研究也表明,水生動(dòng)植物聯(lián)合能取得最佳的水質(zhì)凈化效果。因此,構(gòu)建包含底棲動(dòng)物的復(fù)合人工濕地生態(tài)系統(tǒng)可能是提升人工濕地污染物消納容量的有效途徑,其作用機(jī)理尚值得進(jìn)一步深入研究。

    4 結(jié)論

    對(duì)亞熱帶丘陵區(qū)典型池塘底棲動(dòng)物群落結(jié)構(gòu)的研究表明,該區(qū)域池塘中底棲動(dòng)物多樣性相對(duì)較低,共采集到底棲動(dòng)物22 種,其中以節(jié)肢動(dòng)物居多,耐污性較強(qiáng)的黃色羽搖蚊、霍甫水絲蚓和中華長足搖蚊為主要優(yōu)勢(shì)種。

    池塘底棲動(dòng)物豐度、生物量與水體富營養(yǎng)化程度呈顯著正相關(guān)關(guān)系,同時(shí)也受季節(jié)變化的顯著影響,其中冬季底棲動(dòng)物豐度和生物量顯著高于夏季。

    不同底棲動(dòng)物對(duì)水環(huán)境因子表現(xiàn)出不同的響應(yīng)特征,黃色羽搖蚊生物量與溫度呈顯著正相關(guān),而霍甫水絲蚓和蘇氏尾鰓蚓生物量則與水體總氮、氨氮呈顯著正相關(guān),霍甫水絲蚓和蘇氏尾鰓蚓的豐度和生物量對(duì)于反映該區(qū)域池塘水體面源污染狀況具有一定的指示作用。

    猜你喜歡
    生物量池塘物種
    吃光入侵物種真的是解決之道嗎?
    英語世界(2023年10期)2023-11-17 09:18:18
    輪牧能有效促進(jìn)高寒草地生物量和穩(wěn)定性
    熱鬧的池塘
    回首2018,這些新物種值得關(guān)注
    電咖再造新物種
    汽車觀察(2018年10期)2018-11-06 07:05:26
    池塘
    生物量高的富鋅酵母的開發(fā)應(yīng)用
    瘋狂的外來入侵物種
    夏天的池塘
    基于SPOT-5遙感影像估算玉米成熟期地上生物量及其碳氮累積量
    久久精品久久精品一区二区三区| 亚洲精品一二三| 国产视频一区二区在线看| 老司机亚洲免费影院| 国产成人免费无遮挡视频| 国产精品一区二区在线不卡| svipshipincom国产片| 操美女的视频在线观看| 国产精品一国产av| 亚洲一码二码三码区别大吗| 免费看不卡的av| 国产午夜精品一二区理论片| 天天躁狠狠躁夜夜躁狠狠躁| 午夜免费成人在线视频| 高清欧美精品videossex| 91国产中文字幕| 精品久久久久久久毛片微露脸 | 午夜福利在线免费观看网站| 中文字幕人妻丝袜一区二区| 国产91精品成人一区二区三区 | 赤兔流量卡办理| 亚洲人成网站在线观看播放| 人妻 亚洲 视频| 国产精品.久久久| 国产精品 国内视频| 最黄视频免费看| 久久久国产一区二区| 久久精品久久精品一区二区三区| 老司机亚洲免费影院| 国产欧美日韩一区二区三区在线| 亚洲国产精品999| 国产亚洲精品久久久久5区| 18在线观看网站| 欧美黄色淫秽网站| 人人澡人人妻人| av国产精品久久久久影院| 国产精品一区二区在线不卡| 成人18禁高潮啪啪吃奶动态图| 成人国语在线视频| 色婷婷av一区二区三区视频| 国产成人影院久久av| 男女之事视频高清在线观看 | 国产成人精品久久久久久| 国产熟女欧美一区二区| 啦啦啦 在线观看视频| 男女床上黄色一级片免费看| 亚洲久久久国产精品| 尾随美女入室| 久久亚洲国产成人精品v| 91精品三级在线观看| 大香蕉久久网| 夫妻性生交免费视频一级片| 久久国产精品男人的天堂亚洲| av线在线观看网站| 国产日韩欧美亚洲二区| 成年av动漫网址| 韩国高清视频一区二区三区| 国产欧美日韩一区二区三区在线| a 毛片基地| 精品少妇内射三级| 91精品伊人久久大香线蕉| 久久国产亚洲av麻豆专区| 日日摸夜夜添夜夜爱| 一边亲一边摸免费视频| 成人三级做爰电影| 国产成人欧美在线观看 | 午夜福利一区二区在线看| 国产精品麻豆人妻色哟哟久久| 女人被躁到高潮嗷嗷叫费观| 蜜桃国产av成人99| 久久人人爽人人片av| 老司机深夜福利视频在线观看 | 岛国毛片在线播放| 欧美日本中文国产一区发布| 18禁国产床啪视频网站| 日韩电影二区| 国产精品久久久人人做人人爽| 欧美av亚洲av综合av国产av| 国产成人系列免费观看| av在线app专区| 欧美av亚洲av综合av国产av| 两个人免费观看高清视频| 欧美黄色片欧美黄色片| 黄色一级大片看看| 国产一卡二卡三卡精品| 尾随美女入室| 一本一本久久a久久精品综合妖精| 深夜精品福利| 赤兔流量卡办理| 另类亚洲欧美激情| 国产精品久久久人人做人人爽| 考比视频在线观看| 国产精品久久久人人做人人爽| 国产不卡av网站在线观看| 欧美黑人欧美精品刺激| www.999成人在线观看| 国产在线视频一区二区| 视频区欧美日本亚洲| 国产精品av久久久久免费| 久久av网站| 久久久国产精品麻豆| 夫妻性生交免费视频一级片| 国产av精品麻豆| a级毛片在线看网站| 超色免费av| 69精品国产乱码久久久| 久久这里只有精品19| 国产成人免费观看mmmm| 人人澡人人妻人| 亚洲欧美激情在线| 日韩视频在线欧美| 中文欧美无线码| 久久ye,这里只有精品| 亚洲三区欧美一区| 亚洲视频免费观看视频| 麻豆av在线久日| 亚洲情色 制服丝袜| 亚洲欧美日韩另类电影网站| 欧美人与善性xxx| 亚洲av国产av综合av卡| 黄色怎么调成土黄色| 国产亚洲av高清不卡| 99国产精品免费福利视频| 又大又爽又粗| 欧美精品av麻豆av| 欧美日韩国产mv在线观看视频| 久久久久久人人人人人| 国产免费福利视频在线观看| av福利片在线| 欧美日韩亚洲国产一区二区在线观看 | 97精品久久久久久久久久精品| av国产精品久久久久影院| 高清不卡的av网站| 久久国产亚洲av麻豆专区| 又粗又硬又长又爽又黄的视频| 欧美日韩亚洲国产一区二区在线观看 | xxxhd国产人妻xxx| 精品高清国产在线一区| 国产精品免费视频内射| 丰满饥渴人妻一区二区三| 亚洲精品在线美女| 蜜桃国产av成人99| 亚洲五月婷婷丁香| 亚洲一区二区三区欧美精品| 在线观看免费日韩欧美大片| 免费不卡黄色视频| 成年女人毛片免费观看观看9 | 中文字幕制服av| 亚洲精品中文字幕在线视频| 日韩制服丝袜自拍偷拍| 久久亚洲精品不卡| 国产又色又爽无遮挡免| 黄色一级大片看看| 国产成人精品在线电影| 国产精品一国产av| 免费高清在线观看视频在线观看| 欧美在线黄色| 18禁国产床啪视频网站| 午夜福利免费观看在线| 高清黄色对白视频在线免费看| 成人黄色视频免费在线看| 桃花免费在线播放| 人人澡人人妻人| 久久久精品区二区三区| 日本欧美视频一区| 国产亚洲欧美在线一区二区| 这个男人来自地球电影免费观看| 操出白浆在线播放| 国产成人啪精品午夜网站| 亚洲精品国产av成人精品| 咕卡用的链子| 少妇猛男粗大的猛烈进出视频| 国产黄色视频一区二区在线观看| 晚上一个人看的免费电影| 丰满饥渴人妻一区二区三| 丁香六月天网| 精品人妻在线不人妻| 五月开心婷婷网| 亚洲av电影在线进入| 在线观看免费日韩欧美大片| 两人在一起打扑克的视频| www.精华液| 高清不卡的av网站| 不卡av一区二区三区| 亚洲欧美激情在线| 91成人精品电影| 人妻一区二区av| 老司机影院毛片| 国产精品偷伦视频观看了| 日韩一卡2卡3卡4卡2021年| 亚洲精品国产色婷婷电影| 国产免费一区二区三区四区乱码| 下体分泌物呈黄色| 国产成人欧美| 新久久久久国产一级毛片| 午夜福利,免费看| 天堂中文最新版在线下载| 日韩av在线免费看完整版不卡| 只有这里有精品99| 性少妇av在线| av在线播放精品| 欧美国产精品va在线观看不卡| 欧美性长视频在线观看| 丰满饥渴人妻一区二区三| 欧美亚洲日本最大视频资源| www.自偷自拍.com| av在线老鸭窝| av视频免费观看在线观看| 国产精品一区二区在线观看99| 菩萨蛮人人尽说江南好唐韦庄| 丰满迷人的少妇在线观看| 一区二区三区乱码不卡18| 一边摸一边做爽爽视频免费| av天堂在线播放| 黄片小视频在线播放| cao死你这个sao货| 精品国产一区二区久久| 亚洲欧洲日产国产| 欧美日韩一级在线毛片| 亚洲欧美一区二区三区黑人| 国产成人精品无人区| 欧美日韩成人在线一区二区| 叶爱在线成人免费视频播放| 国产色视频综合| 中国国产av一级| a 毛片基地| 午夜免费成人在线视频| 国产欧美亚洲国产| 久久性视频一级片| 国产亚洲午夜精品一区二区久久| 美国免费a级毛片| 三上悠亚av全集在线观看| 亚洲av电影在线进入| 一级,二级,三级黄色视频| 国产成人91sexporn| 9色porny在线观看| 久久久久国产精品人妻一区二区| 日本欧美视频一区| 精品国产一区二区三区久久久樱花| 国产精品成人在线| 中文字幕色久视频| 亚洲国产最新在线播放| 丝袜脚勾引网站| 精品一区二区三卡| 免费观看人在逋| 国产精品偷伦视频观看了| 爱豆传媒免费全集在线观看| 国产成人免费无遮挡视频| 一级黄片播放器| 欧美精品高潮呻吟av久久| 久久久久精品人妻al黑| 伊人亚洲综合成人网| 观看av在线不卡| 搡老乐熟女国产| 国产精品久久久久久精品电影小说| 十八禁高潮呻吟视频| 老汉色av国产亚洲站长工具| 国产黄色免费在线视频| 精品久久蜜臀av无| 美女福利国产在线| 97人妻天天添夜夜摸| 午夜影院在线不卡| 国产av一区二区精品久久| 丝袜美腿诱惑在线| 亚洲人成电影免费在线| 欧美大码av| 成年人午夜在线观看视频| 在线观看人妻少妇| 美女大奶头黄色视频| 午夜免费鲁丝| 国产在线一区二区三区精| 色网站视频免费| www.熟女人妻精品国产| 蜜桃国产av成人99| 午夜福利免费观看在线| 久久99热这里只频精品6学生| 一级黄片播放器| 欧美黑人精品巨大| 亚洲国产欧美在线一区| 欧美成人精品欧美一级黄| 久久久亚洲精品成人影院| 日本av手机在线免费观看| 91成人精品电影| 精品福利永久在线观看| 大片免费播放器 马上看| 国产淫语在线视频| 国产成人91sexporn| 国产精品一区二区精品视频观看| 久久av网站| 国产男女内射视频| 久久精品国产a三级三级三级| 中文精品一卡2卡3卡4更新| 亚洲三区欧美一区| 男女边摸边吃奶| 国产1区2区3区精品| 蜜桃在线观看..| 午夜老司机福利片| 亚洲成人免费电影在线观看 | 老司机午夜十八禁免费视频| 亚洲三区欧美一区| 国产亚洲欧美在线一区二区| 两性夫妻黄色片| 大片电影免费在线观看免费| 99久久99久久久精品蜜桃| 午夜精品国产一区二区电影| 在线精品无人区一区二区三| 免费看十八禁软件| 免费av中文字幕在线| 少妇 在线观看| 女性生殖器流出的白浆| 一级毛片电影观看| 无遮挡黄片免费观看| 久久性视频一级片| 免费日韩欧美在线观看| 视频区欧美日本亚洲| 久久这里只有精品19| 天天躁日日躁夜夜躁夜夜| 国产高清视频在线播放一区 | 国产片特级美女逼逼视频| 精品人妻1区二区| 天天添夜夜摸| 久久这里只有精品19| 国产成人精品久久二区二区免费| 女人被躁到高潮嗷嗷叫费观| 无遮挡黄片免费观看| 老司机影院毛片| 亚洲人成电影免费在线| 欧美少妇被猛烈插入视频| 免费日韩欧美在线观看| 国产主播在线观看一区二区 | 在线观看国产h片| 女性被躁到高潮视频| 久久精品国产a三级三级三级| 国产精品九九99| 国产极品粉嫩免费观看在线| 天堂8中文在线网| 9热在线视频观看99| 老汉色av国产亚洲站长工具| 波多野结衣一区麻豆| 中国国产av一级| 亚洲精品国产区一区二| av网站免费在线观看视频| 真人做人爱边吃奶动态| 国产又色又爽无遮挡免| 日韩中文字幕视频在线看片| 50天的宝宝边吃奶边哭怎么回事| 免费看av在线观看网站| 亚洲视频免费观看视频| 免费观看av网站的网址| 99久久精品国产亚洲精品| 精品久久久精品久久久| 欧美av亚洲av综合av国产av| 性色av乱码一区二区三区2| 两人在一起打扑克的视频| 伦理电影免费视频| 曰老女人黄片| 成人国产一区最新在线观看 | xxx大片免费视频| 亚洲国产欧美日韩在线播放| 国产成人精品无人区| 欧美国产精品va在线观看不卡| 亚洲国产日韩一区二区| 女人精品久久久久毛片| 在线 av 中文字幕| 一边摸一边做爽爽视频免费| 久热这里只有精品99| 只有这里有精品99| 热99国产精品久久久久久7| 女性生殖器流出的白浆| 久久午夜综合久久蜜桃| 91麻豆精品激情在线观看国产 | 国产免费一区二区三区四区乱码| 国产一区二区激情短视频 | 国产片内射在线| 免费日韩欧美在线观看| 三上悠亚av全集在线观看| 成人午夜精彩视频在线观看| 免费女性裸体啪啪无遮挡网站| 日韩av免费高清视频| 亚洲色图综合在线观看| 国产高清视频在线播放一区 | 国产精品香港三级国产av潘金莲 | 一级a爱视频在线免费观看| 欧美激情极品国产一区二区三区| 久久久久精品国产欧美久久久 | 亚洲精品在线美女| 啦啦啦 在线观看视频| bbb黄色大片| 亚洲精品国产av成人精品| 欧美日韩精品网址| 久久亚洲国产成人精品v| 两性夫妻黄色片| 亚洲精品久久午夜乱码| 嫩草影视91久久| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩视频高清一区二区三区二| av欧美777| 最新在线观看一区二区三区 | 国产亚洲欧美在线一区二区| 亚洲国产精品一区三区| 18禁国产床啪视频网站| 老司机午夜十八禁免费视频| 欧美+亚洲+日韩+国产| 欧美日韩视频精品一区| 国产成人精品久久二区二区免费| 久久精品人人爽人人爽视色| www.熟女人妻精品国产| 飞空精品影院首页| 久久精品人人爽人人爽视色| 国产成人av教育| 99久久99久久久精品蜜桃| 欧美精品高潮呻吟av久久| 国产伦人伦偷精品视频| 婷婷丁香在线五月| 亚洲av综合色区一区| 在线看a的网站| 一本久久精品| 国产成人精品在线电影| 91精品国产国语对白视频| 精品一品国产午夜福利视频| 欧美黄色淫秽网站| 97在线人人人人妻| 欧美日韩亚洲国产一区二区在线观看 | 国产成人免费观看mmmm| 亚洲精品日本国产第一区| 少妇被粗大的猛进出69影院| 国产精品国产三级专区第一集| 91精品国产国语对白视频| 男人添女人高潮全过程视频| 大香蕉久久成人网| 久久青草综合色| 亚洲精品国产av蜜桃| 深夜精品福利| 国产一区二区在线观看av| 一本大道久久a久久精品| 国产主播在线观看一区二区 | 亚洲欧美清纯卡通| 午夜影院在线不卡| 十八禁人妻一区二区| 在线观看免费高清a一片| 啦啦啦啦在线视频资源| 国产午夜精品一二区理论片| 中国美女看黄片| 美女福利国产在线| 91精品三级在线观看| 国产精品欧美亚洲77777| 晚上一个人看的免费电影| 99精国产麻豆久久婷婷| 国产激情久久老熟女| 在线亚洲精品国产二区图片欧美| 国产精品成人在线| 亚洲三区欧美一区| 十八禁人妻一区二区| 亚洲精品国产av成人精品| 国产精品av久久久久免费| 好男人视频免费观看在线| 九草在线视频观看| 国产一区二区 视频在线| 亚洲精品一二三| 最新在线观看一区二区三区 | 亚洲欧美一区二区三区国产| 蜜桃国产av成人99| avwww免费| 只有这里有精品99| xxxhd国产人妻xxx| av欧美777| 中文精品一卡2卡3卡4更新| 久久毛片免费看一区二区三区| 午夜久久久在线观看| 精品一区二区三卡| 一级片'在线观看视频| 国产精品麻豆人妻色哟哟久久| 免费av中文字幕在线| 人妻 亚洲 视频| 夫妻午夜视频| 赤兔流量卡办理| svipshipincom国产片| 97在线人人人人妻| 亚洲欧洲日产国产| 免费看不卡的av| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩免费高清中文字幕av| 另类亚洲欧美激情| 午夜老司机福利片| 国产精品一区二区精品视频观看| 国产亚洲欧美在线一区二区| 亚洲国产中文字幕在线视频| 亚洲精品中文字幕在线视频| 永久免费av网站大全| 男女高潮啪啪啪动态图| 每晚都被弄得嗷嗷叫到高潮| 男女免费视频国产| 人人妻人人添人人爽欧美一区卜| 国产一区有黄有色的免费视频| 久久性视频一级片| 久久精品国产综合久久久| 精品一区二区三卡| 免费观看a级毛片全部| 久久精品亚洲熟妇少妇任你| 亚洲,欧美,日韩| 美女主播在线视频| 精品少妇内射三级| 天天躁夜夜躁狠狠躁躁| 90打野战视频偷拍视频| 中文精品一卡2卡3卡4更新| 老汉色∧v一级毛片| 男女免费视频国产| 日本猛色少妇xxxxx猛交久久| 久久久久久久国产电影| 免费看十八禁软件| 国产老妇伦熟女老妇高清| 久久中文字幕一级| 在线观看免费高清a一片| 国产欧美亚洲国产| 色综合欧美亚洲国产小说| 久久久国产欧美日韩av| 午夜激情久久久久久久| 少妇精品久久久久久久| 青春草亚洲视频在线观看| 人人妻人人爽人人添夜夜欢视频| 脱女人内裤的视频| 亚洲人成77777在线视频| 久久99热这里只频精品6学生| 精品人妻在线不人妻| 桃花免费在线播放| 性色av一级| 男女下面插进去视频免费观看| 啦啦啦在线免费观看视频4| 欧美精品高潮呻吟av久久| 69精品国产乱码久久久| 欧美日韩精品网址| 波野结衣二区三区在线| 黄片播放在线免费| 久久毛片免费看一区二区三区| 欧美大码av| 90打野战视频偷拍视频| 亚洲欧美清纯卡通| 亚洲一码二码三码区别大吗| 国产黄色视频一区二区在线观看| 免费观看人在逋| 国产精品99久久99久久久不卡| 亚洲精品自拍成人| 亚洲精品日本国产第一区| 日本午夜av视频| 亚洲av成人不卡在线观看播放网 | 精品久久久久久久毛片微露脸 | 欧美精品一区二区大全| 大片免费播放器 马上看| 亚洲,欧美精品.| 国产精品免费大片| 99久久人妻综合| 欧美日韩视频精品一区| 老司机在亚洲福利影院| 性色av乱码一区二区三区2| 一级毛片女人18水好多 | 校园人妻丝袜中文字幕| 女人久久www免费人成看片| 99九九在线精品视频| 国产三级黄色录像| 丝袜美足系列| 最近中文字幕2019免费版| 下体分泌物呈黄色| 51午夜福利影视在线观看| 亚洲精品国产色婷婷电影| 午夜影院在线不卡| 国产极品粉嫩免费观看在线| 天天躁日日躁夜夜躁夜夜| 在线天堂中文资源库| 欧美大码av| 久久亚洲精品不卡| 可以免费在线观看a视频的电影网站| 亚洲精品国产区一区二| 日韩大码丰满熟妇| 在线 av 中文字幕| 天天添夜夜摸| 99久久综合免费| 一区二区日韩欧美中文字幕| 欧美日韩亚洲高清精品| 国产精品国产三级专区第一集| 水蜜桃什么品种好| 国产伦人伦偷精品视频| 亚洲av片天天在线观看| 黄片小视频在线播放| 国产激情久久老熟女| 搡老岳熟女国产| 日本vs欧美在线观看视频| 欧美 日韩 精品 国产| 午夜老司机福利片| 免费少妇av软件| 国产高清videossex| 久久亚洲国产成人精品v| 热99国产精品久久久久久7| 日韩一本色道免费dvd| 国产一级毛片在线| 色婷婷av一区二区三区视频| 亚洲国产看品久久| 国产免费视频播放在线视频| 欧美xxⅹ黑人| 一边摸一边做爽爽视频免费| 最新在线观看一区二区三区 | 老鸭窝网址在线观看| 午夜两性在线视频| 国产欧美日韩综合在线一区二区| 蜜桃在线观看..| 亚洲欧洲日产国产| 精品少妇黑人巨大在线播放| 久久精品国产综合久久久| 男女下面插进去视频免费观看| 侵犯人妻中文字幕一二三四区| 久9热在线精品视频| 国产日韩欧美亚洲二区| 啦啦啦在线观看免费高清www| 成人午夜精彩视频在线观看| 校园人妻丝袜中文字幕| 丝袜美腿诱惑在线|