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

    鄱陽湖大型底棲動(dòng)物時(shí)空演變特征及驅(qū)動(dòng)因素

    2021-07-22 12:11:48鄒亮華張慶吉龔志軍張艷杰魯順保蔡永久
    中國環(huán)境科學(xué) 2021年6期
    關(guān)鍵詞:搖蚊底質(zhì)湖區(qū)

    鄒亮華,鄒 偉,張慶吉,3,李 穎,4,龔志軍,張艷杰,魯順保*,蔡永久**

    鄱陽湖大型底棲動(dòng)物時(shí)空演變特征及驅(qū)動(dòng)因素

    鄒亮華1,2,鄒 偉2,張慶吉2,3,李 穎2,4,龔志軍2,張艷杰1,魯順保1*,蔡永久2**

    (1.江西師范大學(xué)生命科學(xué)學(xué)院,江西 南昌 330027;2.中國科學(xué)院南京地理與湖泊研究所,中國科學(xué)院流域地理學(xué)重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210008;3.南京大學(xué)地理與海洋科學(xué)學(xué)院,江蘇 南京 210023;4.西南大學(xué)地理科學(xué)學(xué)院,重慶 400715)

    于2019年7月對鄱陽湖133個(gè)樣點(diǎn)大型底棲動(dòng)物進(jìn)行調(diào)查,分析其群落結(jié)構(gòu)與環(huán)境因子之間的關(guān)系.結(jié)果表明,共記錄底棲動(dòng)物48種,隸屬7綱16目22科38屬,全湖底棲動(dòng)物的平均密度和生物量分別為158.38ind/m2和173.76g/m2,現(xiàn)階段優(yōu)勢種主要為河蜆、銅銹環(huán)棱螺和大沼螺.將全湖分為5個(gè)湖區(qū):北部湖區(qū)、西北部湖區(qū)、中部湖區(qū)、東部湖區(qū)和南部湖區(qū),單因素方差分析表明不同湖區(qū)水體理化指標(biāo)差異顯著(ANOVA,<0.05),沉積物環(huán)境指標(biāo)差異不顯著.相似性分析結(jié)果表明不同湖區(qū)底棲動(dòng)物群落結(jié)構(gòu)差異性顯著,北部湖區(qū)生物多樣性顯著高于其它湖區(qū).典范對應(yīng)分析表明水深、溶解氧、濁度、總磷、葉綠素a、燒失量和底質(zhì)類型與底棲動(dòng)物群落結(jié)構(gòu)顯著相關(guān).與歷史研究相比,鄱陽湖底棲動(dòng)物多樣性下降明顯,優(yōu)勢種由大型軟體動(dòng)物逐漸演變成小型軟體動(dòng)物和昆蟲類.鄱陽湖采砂、水文情勢變化和水生植被衰退是影響鄱陽湖大型底棲動(dòng)物群落結(jié)構(gòu)演替的主要環(huán)境因素.

    鄱陽湖;大型底棲動(dòng)物;群落結(jié)構(gòu);生物多樣性;環(huán)境因子

    底棲動(dòng)物是指生活史的全部或大部分時(shí)間生活于水體底部的無脊椎動(dòng)物,是維系水生生物多樣性、水生態(tài)系統(tǒng)結(jié)構(gòu)和功能的重要組成部分[1].底棲動(dòng)物可加速水底碎屑的分解、調(diào)節(jié)泥水界面的物質(zhì)交換及促進(jìn)水體的自凈等作用,同時(shí)還是水生態(tài)系統(tǒng)食物鏈的重要環(huán)節(jié)[2].底棲動(dòng)物具有分布范圍廣、生命周期長和遷移能力弱、不同物種對環(huán)境變化的敏感程度不同等特點(diǎn)[3],因此被廣泛運(yùn)用到水生態(tài)評價(jià)中.隨著我國對流域水生態(tài)管理的日益重視,水生態(tài)系統(tǒng)健康的研究也更加深入,底棲動(dòng)物作為生態(tài)監(jiān)測的關(guān)鍵類群,需要開展進(jìn)一步的研究[4].

    鄱陽湖位于長江中游,江西省北部,是中國第一大淡水湖,正常水位(14~15m)情況下面積達(dá)3914km2.鄱陽湖是典型的吞吐性、季節(jié)性的湖泊,與贛江、撫河、信江、饒河、修水相接,調(diào)蓄后經(jīng)湖口注入長江,水位漲落變化較大[5].近年來受上游大型水利工程和流域來水影響,改變了長江和鄱陽湖之間的江湖關(guān)系,從而導(dǎo)致鄱陽湖生態(tài)系統(tǒng)格局出現(xiàn)巨大變化[6].此外,近年來鄱陽湖氮、磷污染越來越嚴(yán)重,水體富營養(yǎng)化也進(jìn)一步加劇,水環(huán)境質(zhì)量不斷下降[7],而底棲動(dòng)物對環(huán)境的變化十分敏感,其群落結(jié)構(gòu)往往會(huì)隨著水環(huán)境的變化發(fā)生改變[8].

    對鄱陽湖底棲動(dòng)物的調(diào)查與研究已有一些報(bào)道,20世紀(jì)80年代初,謝欽銘等[9]對鄱陽湖底棲動(dòng)物有過定性定量調(diào)查,調(diào)查樣點(diǎn)共20個(gè),記錄底棲動(dòng)物95種;20世紀(jì)90年代,Wang等[10]對鄱陽湖底棲動(dòng)物進(jìn)行定量調(diào)查,調(diào)查樣點(diǎn)共18個(gè),記錄底棲動(dòng)物58種;歐陽珊等[11]于2007、2008年在鄱陽湖湖區(qū)分別調(diào)查了43個(gè)和32個(gè)樣點(diǎn),采集到大型底棲動(dòng)物35種;Cai等[12]在2012年對鄱陽湖底棲動(dòng)物進(jìn)行4個(gè)季度調(diào)查,調(diào)查樣點(diǎn)共15個(gè),采集到大型底棲動(dòng)物41種,并結(jié)合歷史研究分析表明底棲動(dòng)物優(yōu)勢種發(fā)生明顯變化.然而以往研究均表現(xiàn)出空間分辨率低的特點(diǎn),難以全面反映鄱陽湖底棲動(dòng)物群落結(jié)構(gòu)特征,為此本研究根據(jù)鄱陽湖豐水期水域面積布設(shè)了133個(gè)樣點(diǎn),系統(tǒng)調(diào)查鄱陽湖大型底棲動(dòng)物群落及環(huán)境因子,探討鄱陽湖底棲動(dòng)物群落結(jié)構(gòu)的空間異質(zhì)性及其驅(qū)動(dòng)因子.并與歷史研究進(jìn)行比較,揭示其底棲動(dòng)物群落結(jié)構(gòu)的長期變化,以期為評價(jià)鄱陽湖水生態(tài)系統(tǒng)健康狀況及影響因素提供科學(xué)支持.

    1 材料與方法

    1.1 研究區(qū)域概括

    鄱陽湖(28°22′~29°45′N和115°47′~116°45′E)是我國最大的淡水湖,也是目前長江中游僅存的2個(gè)大型自然通江湖泊之一,位于江西省北部長江南岸.鄱陽湖地處我國亞熱帶濕潤氣候區(qū),氣候溫和,降雨充沛,日照充足,年均氣溫16.5~17.8℃,7月平均氣溫28.4~29.8℃,多年平均年降水量1570mm.鄱陽湖承載贛江、撫河、信江、饒河、修水及漳田河等來水,經(jīng)調(diào)蓄后由湖口注入長江,是一個(gè)過水性、吞吐型、季節(jié)性的湖泊,其具有“高水湖相,低水河相”的獨(dú)特自然地理景觀,其豐、枯水位面積相差10倍[13].鄱陽湖的水位受到支流和湖盆、北部長江的影響,鄱陽湖有明顯的季節(jié)性水位階段,12月到次年3月鄱陽湖星子站水位低于10m(吳淞),湖區(qū)面積通常減少到不到1000km2.4~6月,水位增加至14m,豐水期7~9月,水位高于14m,湖泊面積超過4000km2,10~11月,水位急劇下降至10m左右,并在12月至次年3月保持在低水位階段.作為少數(shù)幾個(gè)自由連接長江的湖泊之一,鄱陽湖在維持長江洪泛生態(tài)系統(tǒng)獨(dú)特生物多樣性方面發(fā)揮著重要作用[14].

    圖1 鄱陽湖采樣點(diǎn)分布

    根據(jù)鄱陽湖的湖泊形態(tài)、水動(dòng)力和生境條件等,把鄱陽湖分為5個(gè)湖區(qū),共設(shè)置了133個(gè)樣點(diǎn)(圖1),分別為北部湖區(qū)(=17)、西北部湖區(qū)(=19)、中部湖區(qū)(=43)、東部湖區(qū)(=30)和南部湖區(qū)(=24).北部湖區(qū)為通江水道區(qū)域,湖面狹窄,流速較快,水較深,受采砂活動(dòng)影響較多;西北部湖區(qū)為半封閉寬闊湖面,流速不均勻,河流入湖口,部分區(qū)域有采砂活動(dòng);中部湖區(qū)湖面寬闊,流速較快;東部湖區(qū)為湖灣區(qū)域,流速緩慢,水面相對靜止,撮萁湖區(qū)域有大量的養(yǎng)殖活動(dòng);南部湖區(qū)流速緩慢,河流入湖口較多,河流帶來的污染物較多.

    1.2 環(huán)境因子測定和樣本采集

    2019年鄱陽湖豐水期時(shí),進(jìn)行野外環(huán)境因子測定和樣品采集.現(xiàn)場使用測深儀測定水深(WD),使用多參數(shù)水質(zhì)監(jiān)測儀(YSI 6600V2)測定pH值、電導(dǎo)率(Cond)、溶解氧(DO)、濁度(Tur),同時(shí)采集水樣和表層沉積物樣品并冷藏保存帶回實(shí)驗(yàn)室.現(xiàn)場使用0.0625m2的改良彼得森采泥器,在每個(gè)采樣點(diǎn)附近3個(gè)不同位置分別抓取1次合成1個(gè)樣品,所采的泥樣經(jīng)60目尼龍篩網(wǎng)篩洗干凈后,剩余物冷藏帶回實(shí)驗(yàn)室進(jìn)行分揀.在實(shí)驗(yàn)室中將剩余物放在白色解剖盤中,將底棲動(dòng)物逐一挑出,樣本用10%福爾馬林保存.利用解剖鏡和顯微鏡將樣本都鑒定至盡可能低的分類單元[15-17],然后用濾紙吸去表面固定液,置于萬分之一電子天平上稱重,將結(jié)果換算成單位面積的密度(ind/m2)和生物量(g/m2).

    在實(shí)驗(yàn)室內(nèi)測定水體銨氮(NH4+-N)、硝酸鹽氮(NO3--N)、正磷酸鹽(PO43-P)、溶解性有機(jī)碳(DOC)、總氮(TN)、總磷(TP)、高錳酸鹽指數(shù)(CODMn)和葉綠素a(Chl-a)濃度;測定沉積物總氮(TNs)、總磷(TPs)、燒失量(LOI)和底質(zhì)類型組成,其中底質(zhì)類型組成按黏土(<4μm)、細(xì)粉砂(4~16μm)、粗粉砂(16~63μm)、砂(>63μm)的所占質(zhì)量百分比劃分.營養(yǎng)鹽參數(shù)的測定方法參考《水和廢水監(jiān)測分析方法》[18].

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

    1.3.1 優(yōu)勢種 優(yōu)勢種是指群落中占優(yōu)勢的且對整個(gè)群落具有控制性影響的種類[19],本文采用相對重要性指數(shù)(IRI)進(jìn)行優(yōu)勢種的確定,其計(jì)算公式為:

    式中:為某一種類的生物量占各點(diǎn)大型底棲動(dòng)物總生物量的比例,%;為該種類的密度占各點(diǎn)總密度的比例,%;為該種類在各湖區(qū)出現(xiàn)的相對頻率.

    1.3.2 多樣性指數(shù) 本文運(yùn)用Simpson優(yōu)勢度指數(shù)()、Margalef豐富度指數(shù)(ma)、Shannon-Wiener多樣性指數(shù)()及Pielou均勻度指數(shù)()這4個(gè)多樣性指數(shù)計(jì)算各樣點(diǎn)的底棲動(dòng)物多樣性.

    式中:為大型底棲動(dòng)物群落總物種數(shù);為大型底棲動(dòng)物群落總個(gè)體數(shù);n為第個(gè)物種的個(gè)體數(shù)

    1.3.3 統(tǒng)計(jì)分析 為提高數(shù)據(jù)的正態(tài)性和方差齊性,將密度、生物量與所有水環(huán)境因子(除pH值外)數(shù)據(jù)進(jìn)行l(wèi)og(+1)轉(zhuǎn)換.使用單因素方差分析(one-way ANOVA)比較各湖區(qū)環(huán)境指標(biāo)間是否存在顯著差異.采用主成分分析(PCA)分析理化數(shù)據(jù),將各環(huán)境因子重新進(jìn)行降維,提取出可以盡可能多的反映原來變量信息的綜合變量,尋找環(huán)境因子的潛在梯度.

    利用相似性分析(ANOSIM)檢驗(yàn)不同湖區(qū)的動(dòng)物群落組成的差異顯著性,利用相似百分比(SIMPER)找出對底棲動(dòng)物群落差異貢獻(xiàn)最大的物種,找出各湖區(qū)的特征種及其相應(yīng)的貢獻(xiàn)率.采用典范對應(yīng)分析(CCA)分析底棲動(dòng)物群落與環(huán)境因子之間的關(guān)系,首先使用向前引入法對鄱陽湖環(huán)境因子進(jìn)行篩選,保留能通過蒙特卡羅置換檢驗(yàn)的顯著因子(<0.05).

    單因素方差分析在SPSS 22中完成,ANOSIM分析和SIMPER在PRIMER 5.0中完成,PCA和CCA在R語言Vegan包中完成,底棲動(dòng)物的空間分布圖使用ArcGIS 10.7繪制.

    2 結(jié)果與分析

    2.1 環(huán)境因子特征

    由圖2、表1可見,One-way ANOVA表明水體理化指標(biāo)在各湖區(qū)均表現(xiàn)出顯著差異(<0.05).北部湖區(qū)的平均水深為11.36m,濁度為20.04NTU, NO3--N濃度為0.66mg/L,均顯著高于其它湖區(qū)(<0.05),但電導(dǎo)率和pH值顯著低于其它湖區(qū)(<0.001).南部湖區(qū)DO平均值為9.83mg/L,顯著高于其它湖區(qū)(<0.01).西北部湖區(qū)的PO43-P、NH4+-N、CODMn、DOC濃度顯著高于其它湖區(qū)(<0.05).

    圖2 鄱陽湖水體和沉積物理化參數(shù)空間格局

    表1 鄱陽湖各湖區(qū)水體和沉積物理化參數(shù)

    中部湖區(qū)的TN濃度為1.40mg/L,顯著高于其它湖區(qū)(<0.05).南部湖區(qū)的TP濃度為0.047mg/L, Chl-a濃度為34.62μg/L,顯著高于其它湖區(qū)(< 0.001).沉積物理化指標(biāo)方面,不同湖區(qū)間無顯著差異(>0.05),南部湖區(qū)的TNs、TPs平均值最高,中部湖區(qū)的LOI平均值最高.各湖區(qū)底質(zhì)粒徑組成中粗粉砂占比較為均勻,西北部湖區(qū)和東部湖區(qū)黏土、細(xì)粉砂占比較高,而北部湖區(qū)、中部湖區(qū)和南部湖區(qū)砂占比較高,接近于其它2個(gè)湖區(qū)的2倍(圖3).

    主成分分析顯示第一主成分與電導(dǎo)率、DOC、CODMn、DO、pH值、TNs、黏土、細(xì)粉砂、粗粉砂呈負(fù)相關(guān),與NO3--N、濁度、水深、砂呈正相關(guān).第二主成分與Chl-a、TP、TN、PO43--P呈正相關(guān),與LOI、TPs呈負(fù)相關(guān)(圖4).不同湖區(qū)樣點(diǎn)在主成分排序圖上能較好地區(qū)分,東部湖區(qū)主要分布于第三象限,沉積物氮磷和有機(jī)質(zhì)含量高,黏土和淤泥含量高.北部湖區(qū)的采樣點(diǎn)主要分布于右側(cè),具有較高的水深和濁度,較低的電導(dǎo)率、有機(jī)污染和營養(yǎng)水平,沉積物砂含量高.西北湖區(qū)的洲灘區(qū)域樣點(diǎn)集中分布在第二象限,水體中DOC、CODMn、NH4+-N濃度較高,表明有機(jī)質(zhì)濃度較高.南部湖區(qū)樣點(diǎn)主要分布于排序圖的中間區(qū)域,表明各項(xiàng)指標(biāo)值從全湖來看處于中等水平.中部湖區(qū)樣點(diǎn)在排序圖中分布較為分散,表明湖區(qū)內(nèi)部理化因子變異較大.

    2.2 種類組成與優(yōu)勢種

    圖3 鄱陽湖底質(zhì)粒徑組成空間格局

    共鑒定到大型底棲動(dòng)物48種,隸屬7綱16目22科38屬,其中軟體動(dòng)物腹足綱4科6屬8種(16.7%),雙殼綱5科10屬12種(25%),環(huán)節(jié)動(dòng)物寡毛綱1科2屬2種(4.2%),多毛綱4科4屬4種(8.3%),蛭綱2科2屬2種(4.2%),節(jié)肢動(dòng)物甲殼綱2科2屬2種(4.2%),昆蟲綱中搖蚊科幼蟲9屬15種(31.3%),其他類3種(6.3%).北部湖區(qū)物種豐富度相對較高,共采集到底棲動(dòng)物38種,其次為中部湖區(qū)和南部湖區(qū),為29種,東部湖區(qū)物種數(shù)為24種,而西北部湖區(qū)物種數(shù)最少,為21種.在所有底棲動(dòng)物中,出現(xiàn)率最高的是銅銹環(huán)棱螺(42.65%),其次分別是河蜆(31.62%)、大沼螺(23.53%)、紋沼螺(16.91%)、長角涵螺(13.97%)、寡鰓齒吻沙蠶(13.24%).根據(jù)鄱陽湖大型底棲動(dòng)物的IRI值表明,河蜆、銅銹環(huán)棱螺和大沼螺在各湖區(qū)均占據(jù)優(yōu)勢地位,為現(xiàn)階段的優(yōu)勢種(表2).

    圖4 鄱陽湖環(huán)境因子主成分分析

    表2 鄱陽湖不同湖區(qū)大型底棲動(dòng)物的優(yōu)勢種及其相對重要性指數(shù)得分

    注:表中數(shù)值為相對重要性指數(shù)(IRI),且僅列出IRI>1的物種.

    2.2 密度和生物量

    全湖大型底棲動(dòng)物的平均密度和平均生物量分別為158.38ind/m2和173.76g/m2.大型底棲動(dòng)物密度空間差異較大(圖5),最高值出現(xiàn)在西北湖區(qū)的樣點(diǎn),達(dá)到1897.14ind/m2,西北湖區(qū)的平均密度要高于其它湖區(qū),為333.34ind/m2,其次為北部湖區(qū),平均密度為288.21ind/m2,而中部湖區(qū)平均密度最低,僅為85.03ind/m2,東部湖區(qū)與南部湖區(qū)平均密度較為接近,分別為110.03和119.74ind/m2.大型底棲生物量空間變化也較大,最高值出現(xiàn)在南部湖區(qū)的一個(gè)樣點(diǎn),達(dá)到4678.11g/m2,與密度分布相同,西北湖區(qū)的平均生物量最高,為354.41g/m2,中部湖區(qū)最低,僅為65.67g/m2,而南部湖區(qū)、北部湖區(qū)、東部湖區(qū)分別為256.96,220.94,120.99g/m2.從不同類群所占比重可以看出,密度方面,腹足綱和雙殼綱是優(yōu)勢類群,平均密度分別為78.29,38.53ind/m2,昆蟲綱也占據(jù)一定比重,平均密度為23.07ind/m2.生物量方面,雙殼綱和腹足綱主導(dǎo)了全湖生物量的空間變化,平均生物量分別為119.06g/m2和54.00g/m2,主要由于軟體動(dòng)物個(gè)體相對較大.

    軟體動(dòng)物在各湖區(qū)均有出現(xiàn),主要分布在湖區(qū)內(nèi)部區(qū)域,而沿岸區(qū)域分布較少.銅銹環(huán)棱螺、大沼螺、河蜆和淡水殼菜是軟體動(dòng)物的優(yōu)勢種,占軟體動(dòng)物總密度和總生物量的80.14%和64.33%.其中銅銹環(huán)棱螺、大沼螺全湖區(qū)廣泛分布,在西北湖區(qū)平均密度和平均生物量最高,中部湖區(qū)較低.河蜆與前二者分布相似,在中部湖區(qū)分布最低,平均密度和平均生物量僅分別為3.45ind/m2、7.06g/m2,在其它湖區(qū)則分布較高.淡水殼菜主要集中分布在北部湖區(qū),其它湖區(qū)相對較少.

    昆蟲綱主要以搖蚊幼蟲為主,全湖搖蚊幼蟲平均密度為19.94ind/m2,其優(yōu)勢種為中國長足搖蚊、德永雕翅搖蚊和黃色羽搖蚊,平均密度為8.95ind/m2,占搖蚊幼蟲總密度的44.87%.搖蚊幼蟲空間分布上差異較大,其中中部湖區(qū)、南部湖區(qū)、北部湖區(qū)搖蚊幼蟲種類較多,平均密度較高,分別為29.53,24.54, 22.24ind/m2,東部湖區(qū)發(fā)現(xiàn)種類較少,搖蚊幼蟲平均密度為12.50ind/m2,西北部湖區(qū)僅發(fā)現(xiàn)中國長足搖蚊,平均密度為3.86ind/m2.除搖蚊幼蟲外,昆蟲綱其它種類較少,主要以楊氏桑嘎蜉為主,北部湖區(qū)和西北部湖區(qū)分布較高,平均密度分別為10.60和9.85ind/m2,而中部湖區(qū)和南部湖區(qū)未發(fā)現(xiàn).

    圖5 鄱陽湖大型底棲動(dòng)物各類群密度和生物量的空間格局

    寡毛綱、多毛綱、蛭綱和甲殼綱在鄱陽湖大型底棲動(dòng)物總密度中所占比例很低(9.95%).寡毛綱僅采集到蘇氏尾鰓蚓和霍甫水絲蚓2種,全湖平均密度分別為2.39和2.68ind/m2,主要分布在西北部湖區(qū)、中部湖區(qū)和東部湖區(qū),其他湖區(qū)密度較低.多毛綱主要分布在各湖區(qū)沿岸帶部分樣點(diǎn)中,其優(yōu)勢種為寡鰓齒吻沙蠶,主要分布在北部湖區(qū)和西北部湖區(qū),其他湖區(qū)密度較低.蛭綱主要分布在中部湖區(qū)及北部湖區(qū)的沿岸帶,東部湖區(qū)并未發(fā)現(xiàn)該類群.甲殼綱在全湖區(qū)出現(xiàn)率較低,優(yōu)勢種為秀麗白蝦,北部湖區(qū)相對分布較多,其它湖區(qū)密度很低.

    2.3 群落結(jié)構(gòu)差異和多樣性

    表3 不同湖區(qū)大型底棲動(dòng)物特征種平均密度(ind/m2)及其對組內(nèi)相似性貢獻(xiàn)百分比(%)

    ANOSIM檢驗(yàn)結(jié)果表明各湖區(qū)之間的底棲動(dòng)物群落存在顯著差異(=0.002),各湖區(qū)的相似性百分比介于80.80%~89.62%之間(表3).北部湖區(qū)的特征種為河蜆、銅銹環(huán)棱螺、大沼螺、寡鰓齒吻沙蠶、楊氏桑嘎蜉、淡水殼菜,SIMPER分析表明這6個(gè)種類對組內(nèi)相似性貢獻(xiàn)率最大(累積達(dá)到90.30%).西北部湖區(qū)的特征種為銅銹環(huán)棱螺、河蜆、大沼螺、中國長足搖蚊、楊氏桑嘎蜉,其對組內(nèi)相似性的累積貢獻(xiàn)率為90.45%.中部湖區(qū)的特征種為銅銹環(huán)棱螺、紋沼螺、大沼螺、背蚓蟲、扁蛭屬一種、河蜆、中國長足搖蚊,其對組內(nèi)相似性的累積貢獻(xiàn)率為90.12%.東部湖區(qū)的特征種為河蜆、銅銹環(huán)棱螺、長角涵螺、大沼螺、紅裸須搖蚊、寡鰓齒吻沙蠶,其對組內(nèi)相似性的累積貢獻(xiàn)率為92.19%.南部湖區(qū)的特征種為銅銹環(huán)棱螺、中國長足搖蚊、紋沼螺、河蜆、蘇氏尾鰓蚓、大沼螺、黃色羽搖蚊,其對組內(nèi)相似性的累積貢獻(xiàn)率為93.26%.

    4個(gè)多樣性指數(shù)在不同湖區(qū)間均存在顯著差異(<0.05,圖6).Simpson指數(shù)全湖區(qū)平均值為0.43,北部湖區(qū)平均值最高(0.68),東部湖區(qū)平均值最低(0.32).Margalef指數(shù)全湖區(qū)平均值為1.03,北部湖區(qū)平均值最高(1.82),東部湖區(qū)平均值最低(0.68). Shannon-Wiener指數(shù)全湖區(qū)平均值為0.80,北部湖區(qū)平均值最高(1.39),東部湖區(qū)平均值最低(0.56). Pielou指數(shù)全湖區(qū)平均值為0.62,北部湖區(qū)平均值最高(0.88),西北部湖區(qū)平均值最低(0.50).總體而言,北部湖區(qū)相較于其它湖區(qū)的底棲動(dòng)物多樣性更高.

    不同字母表示湖區(qū)間具有顯著差異(<0.05)

    2.4 關(guān)鍵環(huán)境影響因子識(shí)別

    DCA分析結(jié)果表明第一軸的軸長5.24,因此選擇單峰模型CCA排序分析.與沉積物參數(shù)作CCA分析篩選出2個(gè)環(huán)境因子,第一排序軸與LOI相關(guān)性較高,第二排序軸與砂(>63μm)相關(guān)性較高(圖7),CCA分析的第1軸和第2軸的特征值分別為0.241和0.161,分別解釋10.4%的物種數(shù)據(jù)方差變異和44.07%的物種-環(huán)境關(guān)系變異(表4).底棲動(dòng)物密度與水體理化因子CCA分析最終篩選出5個(gè)環(huán)境因子.第1排序軸與水深、Chl-a、TP、DO相關(guān)性較高,第2排序軸與濁度相關(guān)性較高,CCA分析的第1軸和第2軸的特征值分別為0.318和0.293,分別解釋11.2%的物種數(shù)據(jù)方差變異和33.95%的物種-環(huán)境關(guān)系變異.

    表4 底棲動(dòng)物群落與環(huán)境因子的典范對應(yīng)分析結(jié)果

    3 討論

    3.1 底棲動(dòng)物群落結(jié)構(gòu)空間異質(zhì)性及影響因素

    不同的區(qū)域環(huán)境因素對底棲動(dòng)物群落結(jié)構(gòu)和現(xiàn)存量會(huì)產(chǎn)生顯著影響.本文研究表明,底質(zhì)類型、水深、溶解氧、濁度、總磷、葉綠素a是影響鄱陽湖大型底棲動(dòng)物分布的主要影響因子.富營養(yǎng)化對湖泊底棲動(dòng)物的影響前人已有充分的研究[20-21],此外,對于過水性湖泊而言,鄱陽湖平均換水周期約為30d[22],湖泊中的底棲藻類將會(huì)更容易被水流帶走,因此對底棲動(dòng)物而言,將可能產(chǎn)生食物限制.而CCA分析結(jié)果表明Chla、TP為驅(qū)動(dòng)底棲動(dòng)物群落結(jié)構(gòu)的主要影響因子,這2個(gè)指標(biāo)除了可以衡量水體的富營養(yǎng)水平外,也可以從一定程度上表征底棲動(dòng)物食物量多少,這一現(xiàn)象在同為過水性湖泊洞庭湖的研究中也有報(bào)道[23].此外,LOI是表征沉積物中有機(jī)質(zhì)含量的指標(biāo),也可以表征底棲動(dòng)物食物數(shù)量和質(zhì)量.由于LOI和Chla、TP均是驅(qū)動(dòng)底棲動(dòng)物群落結(jié)構(gòu)的影響顯著因子,可以推測食物資源可能主要影響著鄱陽湖底棲動(dòng)物群落結(jié)構(gòu).鄱陽湖東部湖區(qū)底棲動(dòng)物均勻度與多樣性顯著低于其它湖區(qū),該區(qū)域底質(zhì)粒徑偏細(xì),粒徑范圍較窄,底質(zhì)異質(zhì)性較低難以為多類底棲動(dòng)物提供棲息場所.北部湖區(qū)與南部湖區(qū)底質(zhì)粒徑范圍較廣,生境多以淤泥+沙混合底質(zhì)為主,為底棲動(dòng)物提供了生存空間和庇護(hù)場所,適合軟體動(dòng)物棲息,同時(shí)也有利于寡毛類和搖蚊幼蟲生存,從而表現(xiàn)出更高的多樣性.底質(zhì)粒徑大小是底質(zhì)類型的重要決定因素.一般而言,底棲動(dòng)物多樣性隨底質(zhì)粒徑增大而發(fā)生顯著變化,且一般不同類型底質(zhì)中的優(yōu)勢類群不同[24].底質(zhì)粒徑大小與有機(jī)質(zhì)和有機(jī)污染物的富集有關(guān),粒徑較小的底質(zhì)與有機(jī)質(zhì)更易結(jié)合,低有機(jī)質(zhì)含量的底質(zhì)無法支撐復(fù)雜的底棲動(dòng)物食物網(wǎng)結(jié)構(gòu)和較高生物量,造成湖泊底棲動(dòng)物群落結(jié)構(gòu)的改變和群落多樣性的顯著降低[25].不同類群底棲動(dòng)物對底質(zhì)類型的喜好差異較大,一般而言,粒徑偏小的底質(zhì)不利于濾食者的攝食,影響軟體動(dòng)物的棲息,如河蜆喜歡砂質(zhì)淤泥底質(zhì),而顫蚓類和搖蚊幼蟲喜好棲居于淤泥底質(zhì)中,因此底質(zhì)類型直接影響底棲動(dòng)物的群落組成和多樣性[26].一般認(rèn)為,從沿岸帶生境多樣性較高,能夠容納更多的動(dòng)物生存,因此底棲動(dòng)物多樣性會(huì)隨著沿岸帶到湖心逐漸降低.本研究中雖然水深也是驅(qū)動(dòng)鄱陽湖底棲動(dòng)物群落結(jié)構(gòu)的主要影響因子,但結(jié)果表明,鄱陽湖深水區(qū)域相較于沿岸洲灘與湖灣淺水區(qū)域底棲動(dòng)物豐富度較高,如軟體動(dòng)物主要分布在湖區(qū)內(nèi)部深水區(qū)域,而沿岸區(qū)域分布較少.出現(xiàn)這種截然相反的情況是因?yàn)檑蛾柡煌竟?jié)水位波動(dòng)較大,在干旱季節(jié),沿岸帶長期出露水面,從而導(dǎo)致其生境不穩(wěn)定,抑制底棲動(dòng)物生存.濁度和溶解氧能夠直接影響底棲動(dòng)物生存,尤其是雙殼類,對濁度和溶解氧要求較高.

    3.2 底棲動(dòng)物群落結(jié)構(gòu)歷史演變

    近年來,受人類活動(dòng)、江湖關(guān)系逐漸改變等影響,鄱陽湖湖泊富營養(yǎng)化程度嚴(yán)重[27],底棲動(dòng)物群落組成與結(jié)構(gòu)也發(fā)生了顯著變化.本次調(diào)查中,共采集底棲動(dòng)物48種,平均密度和生物量為158.38ind/m2和173.76g/m2,優(yōu)勢種為河蜆、銅銹環(huán)棱螺和大沼螺.歷史上對鄱陽湖底棲動(dòng)物進(jìn)行過多次調(diào)查(表5),與20世紀(jì)90年代相比,物種數(shù)、平均密度與生物量下降明顯.另外大型底棲動(dòng)物群落結(jié)構(gòu)變化顯著,優(yōu)勢種呈現(xiàn)小型化的趨勢,從豐富的軟體動(dòng)物演變?yōu)樾⌒蛙涹w動(dòng)物河蜆和耐污能力較強(qiáng)的銅銹環(huán)棱螺,出現(xiàn)了蚌類資源衰退的現(xiàn)象.張屹等[28]對洞庭湖大型底棲動(dòng)物群落演變的研究表明,近30a大型底棲動(dòng)物的密度呈波動(dòng)式下降趨勢,相對較為耐污的種類成為優(yōu)勢種,這一演變特征與鄱陽湖較為相似.此外,由于本次調(diào)查在鄱陽湖布設(shè)了133個(gè)監(jiān)測點(diǎn),因此與歷史研究相比,本次調(diào)查中沿岸帶監(jiān)測點(diǎn)位更多,而螺類在沿岸帶分布較多,這也可能是本次調(diào)查中優(yōu)勢種多數(shù)為螺類的主要原因.

    表5 鄱陽湖大型底棲動(dòng)物的物種數(shù)、優(yōu)勢種、平均密度和平均生物量的年代變化

    底棲動(dòng)物群落結(jié)構(gòu)的變化,與棲息環(huán)境變化和人為干擾程度的加劇密不可分.引起鄱陽湖底棲動(dòng)物群落結(jié)構(gòu)變化可能主要有以下方面原因:一方面是鄱陽湖存在大量采砂作業(yè)的現(xiàn)象,2001~2010年鄱陽湖采砂面積約為260.40km2,挖砂深度平均4.95m, 2005~2006年采砂量為2.36億m3,占長江中下游地區(qū)消費(fèi)量的40%[29-30].采砂造成的底質(zhì)破壞以及水體濁度升高使底棲動(dòng)物棲息地喪失和食源嚴(yán)重不足[31].劉振元等[32]對南洞庭湖軟體動(dòng)物多樣性評估表明,采砂區(qū)物種組成最為單一,采砂是影響軟體動(dòng)物生存的重要制約因素.Meng等[33]對洞庭湖的研究表明采砂活動(dòng)直接導(dǎo)致水深增加,透明度降低,沉積物粒徑改變,軟體動(dòng)物類群喪失賴以生存的棲息地,進(jìn)而導(dǎo)致湖泊底棲動(dòng)物群落結(jié)構(gòu)變得簡單化.劇烈的采砂活動(dòng)導(dǎo)致湖泊底質(zhì)生境穩(wěn)定性較低,部分底棲動(dòng)物直接被清除或死亡.另一方面是鄱陽湖水位變化明顯,研究表明,1990~2016年期間,長江干流水量減少導(dǎo)致鄱陽湖豐水期平均水位下降明顯,而地形改變導(dǎo)致枯水期水位下降,水域面積低至500km2[34].2005年后鄱陽湖水位與之前相比降幅可達(dá)9.4%,枯水期提前且持續(xù)時(shí)間長,10月份水位異常偏低[35].高變幅的水位對底棲動(dòng)物群落將產(chǎn)生變化,特別對湖泊沿岸帶底棲動(dòng)物產(chǎn)生顯著影響[36]. Aroviita等[37]對芬蘭23個(gè)湖泊的對比研究發(fā)現(xiàn),沿岸帶底棲動(dòng)物在不同水位變幅情況下具有顯著差異,物種豐富度隨變幅的增大而降低.一般認(rèn)為,水草是底棲動(dòng)物分布的重要環(huán)境因子.大型沉水植物為底棲動(dòng)物提供了更加多樣的生境,利于腹足類、昆蟲類等底棲動(dòng)物生活、攝食和繁殖.茂盛水草的存在亦為底棲動(dòng)物躲避魚類及其它大型動(dòng)物的掠食以及水流、風(fēng)力的不利影響提供了避難場所.近30a來鄱陽湖長期持續(xù)低枯湖水位,水生植被生存空間被逐漸壓縮,導(dǎo)致水生植被大幅縮減[38],可能也是鄱陽湖底棲動(dòng)物群落結(jié)構(gòu)變化的重要原因.

    4 結(jié)論

    4.1 本次鄱陽湖調(diào)查共采集大型底棲動(dòng)物48種,平均密度為158.38ind/m2,平均生物量為173.76g/m2,底棲動(dòng)物的主要優(yōu)勢種為河蜆、銅銹環(huán)棱螺、大沼螺,不同湖區(qū)的種類數(shù)、優(yōu)勢種均存在顯著差異.

    4.2 空間尺度上,水深、溶解氧、濁度、總磷、葉綠素a、燒失量和底質(zhì)類型是鄱陽湖底棲動(dòng)物群落結(jié)構(gòu)的顯著影響因子.

    4.3 時(shí)間尺度上,與20世紀(jì)90年代相比鄱陽湖底棲動(dòng)物的物種數(shù)呈下降的趨勢,優(yōu)勢種從大型軟體動(dòng)物逐漸演變成小型軟體動(dòng)物,引起鄱陽湖底棲動(dòng)物群落結(jié)構(gòu)變化可能與采砂、水文情勢變化、水生植被衰退等因素有關(guān).

    [1] 王備新,楊蓮芳.大型底棲無脊椎動(dòng)物水質(zhì)快速生物評價(jià)的研究進(jìn)展[J]. 南京農(nóng)業(yè)大學(xué)學(xué)報(bào), 2001,24(4):107-111.

    Wang B X, Yang L F. Advances in rapid bio-assessment of water quality using benthic macroinvertebrates [J]. Journal of Nanjing Agricultural University, 2001,24(4):107-111.

    [2] Covich A P, Palmer M A, Crowl T A, et al. The role of benthic invertebrate species in freshwater ecosystems: zoobenthic species influence energy flows and nutrient cycling [J].BioScience, 2004,54(8):767-775.

    [3] 張 又,劉 凌,蔡永久,等.太湖流域河流及溪流大型底棲動(dòng)物群落結(jié)構(gòu)及影響因素[J]. 中國環(huán)境科學(xué), 2015,35(5):1535-1546.

    Zhang Y, Liu L, Cai Y J, et al. Benthic macroinvertebrate community structure in rivers and streams of Lake Taihu Basin and environmental determinants [J]. China Environmental Science, 2015,35(5):1535-1546.

    [4] 金小偉,王業(yè)耀,王備新,等.我國流域水生態(tài)完整性評價(jià)方法構(gòu)建[J]. 中國環(huán)境監(jiān)測, 2017,33(1):75-81.

    Jin X W, Wang Y Y, Wang B X, et al. Methods development for monitoring and assessment of ecological integrity of surface waters in China [J]. Environmental Monitoring in China, 2017,33(1):75-81.

    [5] Li B, Yang G S, Wan R R, et al. Comparison of random forests and other statistical methods for the prediction of lake water level: a case study of the Poyang Lake in China [J]. Hydrology Research, 2016, 47(S1):69-83.

    [6] Zhang Z X, Chen X, Xu C Y, et al. Examining the influence of river–lake interaction on the drought and water resources in the Poyang Lake basin [J]. Journal of Hydrology, 2015,522:510-521.

    [7] 張 琍,陳曉玲,張 媛,等.水文地貌分區(qū)下鄱陽湖豐水期水質(zhì)空間差異及影響機(jī)制[J]. 中國環(huán)境科學(xué), 2014,34(10):2637-2645.

    Zhang L, Chen X L, Zhang Y, et al. Spatial distribution of water quality and its impacting factor in the wet season of Poyang Lake using the hydro-geomorphological partitions [J]. China Environmental Science, 2014,34(10):2637-2645.

    [8] 王偉莉,閆振廣,何 麗,等.五種底棲動(dòng)物對優(yōu)控污染物的敏感性評價(jià) [J]. 中國環(huán)境科學(xué), 2013,33(10):1856-1862.

    Wang W L, Yan Z G, He L, et al. Sensitivity evaluation of five zoobenthos to priority pollutants [J]. China Environmental Science, 2013,33(10):1856-1862.

    [9] 謝欽銘,李 云,熊國根.鄱陽湖底棲動(dòng)物生態(tài)研究及其底層魚產(chǎn)力的估算[J]. 江西科學(xué), 1995,13(3):161-170.

    Xie Q M, Li Y, Xiong G G. The ecological study of macrozoobenthic and the estimation of demersal fish productivity in Poyang Lake [J]. Jiangxi Science, 1995,13(3):161-170.

    [10] Wang H Z, Xu Q Q, Cui Y D, et al. Macrozoobenthic community of Poyang Lake, the largest freshwater lake of China, in the Yangtze floodplain. [J]. Limnology, 2007,8(1):65-71.

    [11] 歐陽珊,詹 誠,陳堂華,等.鄱陽湖大型底棲動(dòng)物物種多樣性及資源現(xiàn)狀評價(jià)[J]. 南昌大學(xué)學(xué)報(bào)(工科版), 2009,31(1):9-13.

    Ou Y S, Zhan C, Chen T H, et al. Species diversity and resource assessment of macrozoobenthos in Poyang Lake [J]. Journal of Nanchang University(Engineering & Technology), 2009,31(1):9-13.

    [12] Cai Y J, Lu Y J, Wu Z S, et al. Community structure and decadal changes in macrozoobenthic assemblages in Lake Poyang, the largest freshwater lake in China [J]. Knowledge and Management of Aquatic Ecosystems, 2014,414:1-18.

    [13] 楊桂山,馬榮華,張 路,等.中國湖泊現(xiàn)狀及面臨的重大問題與保護(hù)策略[J]. 湖泊科學(xué), 2010,22(6):799-810.

    Yang G S, Ma R H, Zhang L, et al. Lake status, major problems and protection strategy in China [J]. Journal of Lake Sciences, 2010,22(6): 799-810.

    [14] 錢奎梅,劉寶貴,陳宇煒.鄱陽湖浮游植物功能群的長期變化特征(2009~2016年) [J]. 湖泊科學(xué), 2019,31(4):1035-1044.

    Qian K M, Liu B G, Chen Y W. Long term dynamics of phytoplankton functional groups in Lake Poyang during 2009~2016 [J]. Journal of Lake Sciences, 2019,31(4):1035-1044.

    [15] 劉月英,張文珍,王躍先,等.中國經(jīng)濟(jì)動(dòng)物志-淡水軟體動(dòng)物[M]. 北京:科學(xué)出版社, 1979.

    Liu Y Y, Zhang W Z, Wang Y X, et al. Economic fauna of China (freshwater mollusk) [M]. Beijing:Science Press, 1979.

    [16] Morse J C, Yang L F, Tian L X. Aquatic insects of China useful for monitoring water quality [M]. Nanjing: Hohai University Press, 1994.

    [17] 王俊才,王新華.中國北方搖蚊幼蟲[M]. 北京:中國言實(shí)出版社, 2011.

    Wang J C, Wang X H. Northern China Chironomus larva [M]. Beijing:Yanshi Press in China, 2011.

    [18] 國家環(huán)境保護(hù)總局《水和廢水監(jiān)測分析方法》編委會(huì).水和廢水監(jiān)測分析方法[M]. 4版.北京:中國環(huán)境科學(xué)出版社, 2002.

    “Water and wastewater monitoring and analysis method”editorial board of State Environmental Protection Administration of China. Monitoring and analysis methods of water and wastewater [M]. fourth edition.Beijing:China Environmental Science Press, 2002.

    [19] 韓 潔,張志南,于子山.渤海中、南部大型底棲動(dòng)物的群落結(jié)構(gòu)[J]. 生態(tài)學(xué)報(bào), 2004,24(3):531-537.

    Han J, Zhang Z N, Yu Z S. Macrobethic community structure in the southern and central Bohai Sea, China [J]. Acta Ecologica Sinica, 2004,24(3):531-537.

    [20] 龔志軍,謝 平,唐匯涓,等.水體富營養(yǎng)化對大型底棲動(dòng)物群落結(jié)構(gòu)及多樣性的影響[J]. 水生生物學(xué)報(bào), 2001,25(3):210-216.

    Gong Z J, Xie P, Tang H J, et al. The influence of eutrophycation upon community structure and biodiversity of macrozoobenthos [J]. Acta Hydrobiologica Sinica, 2001,25(3):210-216.

    [21] 蔡永久,龔志軍,秦伯強(qiáng).太湖大型底棲動(dòng)物群落結(jié)構(gòu)及多樣性[J]. 生物多樣性, 2010,18(1):50-59.

    Cai Y J, Gong Z J, Qin B Q. Community structure and diversity of macrozoobenthos in Lake Taihu, a large shallow eutrophic lake in China [J]. Biodiversity Science, 2010,18(1):50-59.

    [22] 李云良,姚 靜,李夢凡,等.鄱陽湖換水周期與示蹤劑傳輸時(shí)間特征的數(shù)值模擬[J]. 湖泊科學(xué), 2017,29(1):32-42.

    Li Y L, Yao J, Li M F, et al. Numerical simulations of residence time and tracer travel time in Lake Poyang [J]. Journal of Lake Sciences, 2017,29(1):32-42.

    [23] 王小毛,歐伏平,王丑明,等.洞庭湖底棲動(dòng)物長期演變特征及影響因素分析[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2016,35(2):336-345.

    Wang X M, Ou F P, Wang C M, et al. Long-term evolution and influencing factors of macrozoobenthos in Dongting Lake [J]. Journal of Agro-Environment Science, 2016,35(2):336-345.

    [24] 段學(xué)花,王兆印,田世民.河床底質(zhì)對大型底棲動(dòng)物多樣性影響的野外試驗(yàn)[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版), 2007,47(9):1553-1556.

    Duan X H, Wang Z Y, Tian S M. Field experiment on the effect of streambed substrate on macroinvertebrate diversity [J]. Journal of Tsinghua University(Science and Technology), 2007,47(9):1553-1556.

    [25] Thompson B, Lowe S. Assessment of macrobenthos response to sediment contamination in the San Francisco Estuary,California,USA [J]. Environmental Toxicology and Chemistry, 2004,23(9):2178-2187.

    [26] 鄒 偉,李太民,劉 利,等.蘇北駱馬湖大型底棲動(dòng)物群落結(jié)構(gòu)及水質(zhì)評價(jià)[J]. 湖泊科學(xué), 2017,29(5):1177-1187.

    Zou W, Li T M, Liu L, et al. Macrozoobenthic community structure and water quality assessment of Lake Luoma, Jiangsu Province,China [J]. Journal of Lake Sciences, 2017,29(5):1177-1187.

    [27] 胡春華,周文斌,王毛蘭,等.鄱陽湖氮磷營養(yǎng)鹽變化特征及潛在性富營養(yǎng)化評價(jià)[J]. 湖泊科學(xué), 2010,22(5):723-728.

    Hu C H, Zhou W B, Wang M L, et al. norganic nitrogen and phosphate and potential eutrophication assessment in Lake Poyang [J]. Journal of Lake Sciences, 2010,22(5):723-728.

    [28] 張 屹,王丑明,石慧華,等.洞庭湖近30年大型底棲動(dòng)物的群落演變[J]. 生態(tài)環(huán)境學(xué)報(bào), 2015,24(8):1348-1353.

    Zhang Y, Wang C M, Shi H H, et al. Macrobenthic community succession during last thirty years in Dongting Lake [J]. Ecology and Environmental Sciences, 2015,24(8):1348-1353.

    [29] 江 豐,齊述華,廖富強(qiáng),等.2001~2010年鄱陽湖采砂規(guī)模及其水文泥沙效應(yīng)[J]. 地理學(xué)報(bào), 2015,70(5):837-845.

    Jiang F, Qi S H, Liao F Q, et al. Hydrological and sediment effects from sand mining in Poyang Lake during 2001~2010 [J]. Acta Geographica Sinica, 2015,70(5):837-845.

    [30] Lai X J, David S, Claire H, et al. Sand mining and increasing Poyang Lake’s discharge ability: A reassessment of causes for lake decline in China [J]. Journal of Hydrology, 2014,519:1698-1706.

    [31] Zou W, Tolonen K T, Zhu G W, et al. Catastrophic effects of sand mining on macroinvertebrates in a large shallow lake with implications for management [J]. Science of the Total Environment, 2019,695: 133706.DOI:10.1016/j.scitotenv.2019.133706.

    [32] 劉振元,孟星亮,李正飛,等.南洞庭湖區(qū)軟體動(dòng)物物種多樣性評估及保護(hù)對策[J]. 生物多樣性, 2020,28(2):155-165.

    Liu Z Y, Meng X L, Li F F, et al. Diversity assessment and protection strategies for the mollusk community in the southern Dongting Lake [J]. Biodiversity Science, 2020,28(2):155-165.

    [33] Meng X L, Jiang X M, Li Z F,et al.Responses of macroinvertebrates and local environment to short-term commercial sand dredging practices in a flood-plain lake [J]. Science of the Total Environment, 2018,631-632:1350-1359.DOI:10.1016/j.scitotenv.2018.03.086.

    [34] 孫芳蒂,馬榮華.鄱陽湖水文特征動(dòng)態(tài)變化遙感監(jiān)測[J]. 地理學(xué)報(bào), 2020,75(3):544-557.

    Sun F D, Ma R H. Hydrologic changes of Poyang Lake based on radar altimeter and optical sensor [J]. Acta Geographica Sinica, 2020,75(3): 544-557.

    [35] 鄧志民,張 翔,肖 洋,等.鄱陽湖水位演變及其影響因子分析[J]. 武漢大學(xué)學(xué)報(bào)(工學(xué)版), 2015,48(5):615-621.

    Deng Z M, Zhang X, Xiao Y, et al. Study of evolution of water level in Poyang Lake and impact factors [J]. Engineering Journal of Wuhan University, 2015,48(5):615-621.

    [36] 彭建華,鄭志偉,鄒 曦,等.漢豐湖底棲動(dòng)物與環(huán)境因子間相關(guān)性初步研究[J]. 環(huán)境影響評價(jià), 2015,37(5):63-68.

    Peng J H, Zheng Z W, Zou X, et al. Study of correlation between benthos animals and environmental factors in Hanfeng Lake [J]. Environmental Impact Assessment, 2015,37(5):63-68.

    [37] Aroviita J, H?m?l?inen H. The impact of water-level regulation on littoral macroinvertebrate assemblages in boreal lakes [J]. Hydrobiologia, 2008,613(1):45-56.

    [38] 胡振鵬,林玉茹.鄱陽湖水生植被30年演變及其驅(qū)動(dòng)因素分析[J]. 長江流域資源與環(huán)境, 2019,28(8):1947-1955.

    Hu Z P, Lin Y R. Analysis of evolution process and driving factors for aquatic vegetations of Poyang Lake in 30years [J]. Resources and Environment in the Yangtze Basin, 2019,28(8):1947-1955.

    Characteristics and driving factors of spatiotemporal succession of macrozoobenthos in Poyang Lake.

    ZOU Liang-hua1,2, ZOU Wei2, ZHANG Qing-ji2,3, LI Ying2,4, GONG Zhi-jun2, ZHANG Yan-jie1, LU Shun-bao1*, CAI Yong-jiu2**

    (1.College of Life Sciences, Jiangxi Normal University, Nanchang 330027, China;2.Key Laboratory of Watershed Geographic Sciences, Nanjing Institute of Geography & Limnology, Chinese Academy of Sciences, Nanjing 210008, China;3. School of Geography & Ocean Science, Nanjing University , Nanjing 210023, China;4.School of Geography science, Southwest University, Chongqing 400715, China)., 2021,41(6):2881~2892

    Macrozoobenthos community and environmental variables were investigated at 133 sites in Poyang Lake in July 2019. The relationships between macrozoobenthic community structure and environmental factors were examined. A total of 48 species belonging to 7 classes, 16 orders, 22 families and 38 genera were recorded. The average density and biomass of macrozoobenthos were 158.38 ind/m2 and 173.76g/m2, respectively, mainly dominated by,and. In terms of spatial patterns, the whole lake was divided into five areas, including northern, northwest, central, eastern and southern lake areas. One-way ANOVA analysis detected significant differences of physicochemical variables among the five lake areas, while no significant differences in sediment parameters. The community structure of benthic animals among the five lake areas also varied significantly as indicated by one-way analysis of similarity (ANOSIM). Regarding to biodiversity indexes, significantly higher values in northern lake were found than those in other lake areas. Canonical correspondence analyses (CCA) showed that water depth, dissolved oxygen, turbidity, total phosphorus, chlorophyll a, loss on ignition, and sediment type were significantly related to the community structure of macrozoobenthos. In the long term, the biodiversity of macrozoobenthos in Poyang Lake has declined significantly compared with those in their historical condition, with the dominated species changed from large molluscs to small molluscs and insects. Sand mining, changes in hydrological conditions and decline of aquatic vegetation were the main environmental factors regulating the spatiotemporal succession of the macrozoobenthic community structure in Poyang Lake.

    Poyang Lake;macrozoobenthos;community structure;biodiversity;environmental factors

    X171

    A

    1000-6923(2021)06-2881-12

    2020-09-20

    中國科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)(A類)(XDA23040203);國家自然科學(xué)基金資助項(xiàng)目(32071572,31660072,31760136);江西省自然科學(xué)基金資助項(xiàng)目(20201BAB204237);江西省教育廳基金(GJJ160312, GJJ170198)

    * 責(zé)任作者, 教授, luxunbao8012@126.com; ** 副研究員, caiyj@niglas.ac.cn

    鄒亮華(1994-),男,江西贛州人,碩士研究生,主要從事濕地生態(tài)學(xué)研究.

    猜你喜歡
    搖蚊底質(zhì)湖區(qū)
    不同發(fā)育階段中華絨螯蟹對底質(zhì)的喜好性研究
    一種ReliefF和隨機(jī)森林模型組合的多波束海底底質(zhì)分類方法
    大通湖區(qū)河蟹產(chǎn)業(yè)發(fā)展綜述
    用于海底目標(biāo)識(shí)別與底質(zhì)分類的多波束水體波形預(yù)處理
    云南省雙江縣4種多足搖蚊翅膀形態(tài)特征及COI基因序列分析
    沉睡搖蚊的生存智慧
    生活在湖區(qū)
    海峽旅游(2018年4期)2018-06-01 11:20:00
    湖區(qū)航道風(fēng)速預(yù)警監(jiān)測點(diǎn)布設(shè)研究
    江西建材(2018年4期)2018-04-10 12:37:24
    文蛤的底質(zhì)選擇性及潛沙能力研究
    常用基因序列在搖蚊科昆蟲系統(tǒng)發(fā)育研究中的應(yīng)用進(jìn)展
    国产高清videossex| 最新在线观看一区二区三区| 男人舔女人下体高潮全视频| 国产亚洲精品综合一区在线观看 | 亚洲自拍偷在线| 欧美成人免费av一区二区三区| 国产主播在线观看一区二区| 免费一级毛片在线播放高清视频 | 亚洲熟女毛片儿| 欧美乱色亚洲激情| 国产亚洲精品av在线| 黄色a级毛片大全视频| av天堂久久9| 亚洲国产精品久久男人天堂| 校园春色视频在线观看| 欧美色欧美亚洲另类二区 | 久热爱精品视频在线9| 老司机深夜福利视频在线观看| 欧美一级a爱片免费观看看 | 精品国产美女av久久久久小说| 欧美在线黄色| 大型黄色视频在线免费观看| 激情视频va一区二区三区| 九色亚洲精品在线播放| 搞女人的毛片| 99国产精品一区二区三区| 真人做人爱边吃奶动态| 日日夜夜操网爽| 日韩大尺度精品在线看网址 | 伊人久久大香线蕉亚洲五| 亚洲 国产 在线| 亚洲伊人色综图| 久久婷婷成人综合色麻豆| 精品国产超薄肉色丝袜足j| 成人特级黄色片久久久久久久| 亚洲 欧美 日韩 在线 免费| 国产亚洲精品久久久久久毛片| 精品少妇一区二区三区视频日本电影| 亚洲精品一卡2卡三卡4卡5卡| 不卡av一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 亚洲最大成人中文| 日韩国内少妇激情av| 欧美日本视频| 精品国产国语对白av| 色播在线永久视频| 香蕉久久夜色| 99国产精品99久久久久| 黄色毛片三级朝国网站| 韩国av一区二区三区四区| 久久精品国产综合久久久| 黑人操中国人逼视频| 波多野结衣巨乳人妻| 精品卡一卡二卡四卡免费| 香蕉久久夜色| 国产熟女xx| 韩国av一区二区三区四区| 国产免费男女视频| 老司机深夜福利视频在线观看| 久久精品亚洲熟妇少妇任你| 搡老熟女国产l中国老女人| 亚洲av片天天在线观看| 一区二区三区国产精品乱码| 久9热在线精品视频| 亚洲自拍偷在线| 亚洲国产日韩欧美精品在线观看 | 免费看美女性在线毛片视频| 日韩精品青青久久久久久| 亚洲av成人一区二区三| 男人操女人黄网站| 国产精品一区二区三区四区久久 | 国产成人影院久久av| 1024香蕉在线观看| 久久久久久人人人人人| 伊人久久大香线蕉亚洲五| 日本五十路高清| 日本欧美视频一区| 国产蜜桃级精品一区二区三区| 宅男免费午夜| 亚洲精品国产精品久久久不卡| 亚洲av电影在线进入| 精品国产乱码久久久久久男人| 久久久久久久精品吃奶| 嫁个100分男人电影在线观看| 91国产中文字幕| 成年人黄色毛片网站| 99香蕉大伊视频| 无人区码免费观看不卡| www国产在线视频色| 成人18禁在线播放| 九色国产91popny在线| 看黄色毛片网站| 国产精品亚洲美女久久久| 手机成人av网站| 欧美性长视频在线观看| 国语自产精品视频在线第100页| 国产精品98久久久久久宅男小说| 9191精品国产免费久久| 看黄色毛片网站| 色婷婷久久久亚洲欧美| 亚洲午夜理论影院| 好看av亚洲va欧美ⅴa在| 欧美 亚洲 国产 日韩一| 欧美精品啪啪一区二区三区| 国语自产精品视频在线第100页| 激情视频va一区二区三区| 国产欧美日韩一区二区三区在线| 午夜福利18| 在线观看免费视频网站a站| 精品国产美女av久久久久小说| АⅤ资源中文在线天堂| 99国产综合亚洲精品| a级毛片在线看网站| 中文字幕人妻熟女乱码| 欧美丝袜亚洲另类 | 婷婷六月久久综合丁香| 色综合亚洲欧美另类图片| 国产免费av片在线观看野外av| 亚洲激情在线av| 变态另类丝袜制服| 一夜夜www| 日韩大码丰满熟妇| 久久九九热精品免费| 真人一进一出gif抽搐免费| 后天国语完整版免费观看| 国产精品美女特级片免费视频播放器 | 美女 人体艺术 gogo| 国产精品美女特级片免费视频播放器 | 亚洲成人精品中文字幕电影| 男人舔女人下体高潮全视频| 18禁国产床啪视频网站| 国产黄a三级三级三级人| 黄色视频,在线免费观看| 国产成人精品无人区| 国产精品久久视频播放| 在线观看午夜福利视频| 麻豆国产av国片精品| 国产欧美日韩综合在线一区二区| 97超级碰碰碰精品色视频在线观看| 老司机福利观看| 中文字幕久久专区| 亚洲中文日韩欧美视频| 欧美日韩瑟瑟在线播放| 国产午夜福利久久久久久| 大型黄色视频在线免费观看| 国产又色又爽无遮挡免费看| 久久天堂一区二区三区四区| 亚洲精品粉嫩美女一区| 国产精品爽爽va在线观看网站 | 激情视频va一区二区三区| 亚洲欧美激情在线| 神马国产精品三级电影在线观看 | www.精华液| 一级毛片女人18水好多| 好看av亚洲va欧美ⅴa在| 欧美精品啪啪一区二区三区| 国产午夜福利久久久久久| 国产精品久久久av美女十八| 美女免费视频网站| 少妇粗大呻吟视频| 日韩欧美国产一区二区入口| 亚洲国产精品sss在线观看| 亚洲熟女毛片儿| 国产精品免费视频内射| 欧美 亚洲 国产 日韩一| 亚洲五月色婷婷综合| 久久欧美精品欧美久久欧美| 日本 av在线| 欧美日韩中文字幕国产精品一区二区三区 | 一区福利在线观看| 美女高潮到喷水免费观看| 两个人视频免费观看高清| 久9热在线精品视频| 丰满人妻熟妇乱又伦精品不卡| 精品高清国产在线一区| av免费在线观看网站| 涩涩av久久男人的天堂| 亚洲成人国产一区在线观看| 亚洲国产毛片av蜜桃av| 一区福利在线观看| 淫秽高清视频在线观看| 亚洲精华国产精华精| 国产成人av教育| svipshipincom国产片| 国产成人精品久久二区二区免费| 一边摸一边抽搐一进一出视频| 午夜视频精品福利| 亚洲少妇的诱惑av| 亚洲全国av大片| 一级a爱片免费观看的视频| 日本精品一区二区三区蜜桃| 久久人妻福利社区极品人妻图片| 日韩高清综合在线| 欧美丝袜亚洲另类 | 欧美av亚洲av综合av国产av| 国产xxxxx性猛交| 99国产精品免费福利视频| 18禁美女被吸乳视频| 99精品久久久久人妻精品| 美女扒开内裤让男人捅视频| 黑人操中国人逼视频| 亚洲 欧美 日韩 在线 免费| 黄色女人牲交| or卡值多少钱| 一级作爱视频免费观看| 美女扒开内裤让男人捅视频| 99国产精品99久久久久| 变态另类成人亚洲欧美熟女 | 真人一进一出gif抽搐免费| 亚洲av电影不卡..在线观看| 麻豆一二三区av精品| 亚洲熟女毛片儿| 精品人妻在线不人妻| 成人免费观看视频高清| 禁无遮挡网站| 亚洲av电影不卡..在线观看| 男男h啪啪无遮挡| 黑人巨大精品欧美一区二区mp4| 国产精品影院久久| 99国产精品免费福利视频| 久久久久久久精品吃奶| 两人在一起打扑克的视频| 免费观看精品视频网站| 国产精品98久久久久久宅男小说| 日韩成人在线观看一区二区三区| 久久久水蜜桃国产精品网| 在线观看午夜福利视频| 中出人妻视频一区二区| 亚洲中文日韩欧美视频| 身体一侧抽搐| 国产精品久久久人人做人人爽| 亚洲无线在线观看| 给我免费播放毛片高清在线观看| 免费久久久久久久精品成人欧美视频| 91在线观看av| 午夜精品国产一区二区电影| 午夜免费观看网址| 亚洲av熟女| 亚洲一区高清亚洲精品| 国产成人免费无遮挡视频| 老汉色∧v一级毛片| 男女做爰动态图高潮gif福利片 | 如日韩欧美国产精品一区二区三区| 老司机午夜福利在线观看视频| 91在线观看av| 国产精品久久久久久亚洲av鲁大| 亚洲熟妇中文字幕五十中出| 国产成人影院久久av| 成人三级做爰电影| 欧美乱码精品一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 中文亚洲av片在线观看爽| 巨乳人妻的诱惑在线观看| 少妇裸体淫交视频免费看高清 | 一级黄色大片毛片| 久久狼人影院| 国产成人精品无人区| 国产亚洲欧美98| 看黄色毛片网站| 波多野结衣高清无吗| 叶爱在线成人免费视频播放| 老司机午夜十八禁免费视频| 亚洲av成人一区二区三| 嫩草影视91久久| 国产精品九九99| bbb黄色大片| 国产亚洲精品一区二区www| 免费少妇av软件| 久久天躁狠狠躁夜夜2o2o| 欧美国产日韩亚洲一区| 午夜福利18| 又黄又爽又免费观看的视频| 国产真人三级小视频在线观看| 免费在线观看黄色视频的| 美女国产高潮福利片在线看| 国产av精品麻豆| 久久人人爽av亚洲精品天堂| avwww免费| 性少妇av在线| 自线自在国产av| 嫁个100分男人电影在线观看| 欧美绝顶高潮抽搐喷水| 久久久国产精品麻豆| 九色国产91popny在线| 巨乳人妻的诱惑在线观看| bbb黄色大片| 亚洲欧美激情在线| 中国美女看黄片| 老司机深夜福利视频在线观看| 19禁男女啪啪无遮挡网站| 国产精品 欧美亚洲| av天堂在线播放| 成人手机av| 国产精品乱码一区二三区的特点 | 欧美日韩亚洲国产一区二区在线观看| 97人妻天天添夜夜摸| 日韩精品免费视频一区二区三区| 十八禁人妻一区二区| 不卡av一区二区三区| 亚洲av日韩精品久久久久久密| 亚洲天堂国产精品一区在线| 午夜福利一区二区在线看| 最新在线观看一区二区三区| 日韩精品青青久久久久久| 欧美乱妇无乱码| 国产精品1区2区在线观看.| 欧美日韩亚洲综合一区二区三区_| 啦啦啦韩国在线观看视频| 久久精品国产清高在天天线| 亚洲国产中文字幕在线视频| 丁香欧美五月| 首页视频小说图片口味搜索| 久久欧美精品欧美久久欧美| 一进一出抽搐动态| 日本免费一区二区三区高清不卡 | 国产精品影院久久| 人人妻人人澡人人看| 久久精品91蜜桃| 亚洲精华国产精华精| 亚洲最大成人中文| 日韩欧美免费精品| 精品久久久精品久久久| 露出奶头的视频| 亚洲色图综合在线观看| 免费高清视频大片| 国内精品久久久久精免费| 精品人妻在线不人妻| 老司机午夜福利在线观看视频| 波多野结衣高清无吗| 一区在线观看完整版| 午夜激情av网站| 亚洲欧美一区二区三区黑人| 亚洲va日本ⅴa欧美va伊人久久| 嫩草影院精品99| 日日摸夜夜添夜夜添小说| 岛国视频午夜一区免费看| 国内久久婷婷六月综合欲色啪| 欧美av亚洲av综合av国产av| 十分钟在线观看高清视频www| 成人国产一区最新在线观看| 国产成人av教育| 12—13女人毛片做爰片一| 午夜福利视频1000在线观看 | 久99久视频精品免费| 久久久久久免费高清国产稀缺| 亚洲国产欧美一区二区综合| 亚洲国产精品合色在线| 亚洲aⅴ乱码一区二区在线播放 | 午夜日韩欧美国产| 亚洲人成77777在线视频| 成人欧美大片| 久久久久久久午夜电影| 日日摸夜夜添夜夜添小说| 99国产精品免费福利视频| 久久影院123| 国产欧美日韩一区二区三| 亚洲男人的天堂狠狠| 国产xxxxx性猛交| www国产在线视频色| 女人高潮潮喷娇喘18禁视频| 欧美丝袜亚洲另类 | 亚洲成国产人片在线观看| 精品人妻1区二区| 少妇 在线观看| 亚洲欧美一区二区三区黑人| 日韩精品中文字幕看吧| 亚洲黑人精品在线| 午夜激情av网站| 啪啪无遮挡十八禁网站| 久久久久国产精品人妻aⅴ院| 久久人妻熟女aⅴ| 亚洲一区二区三区不卡视频| 久久人妻熟女aⅴ| 久久精品亚洲精品国产色婷小说| 日本a在线网址| 午夜福利影视在线免费观看| 精品免费久久久久久久清纯| 亚洲一区中文字幕在线| 亚洲精华国产精华精| 九色亚洲精品在线播放| www.www免费av| 国语自产精品视频在线第100页| 99国产精品一区二区蜜桃av| 亚洲色图av天堂| 最近最新中文字幕大全免费视频| 免费无遮挡裸体视频| 日韩av在线大香蕉| 久热爱精品视频在线9| 久久久久久久久久久久大奶| 国产亚洲av高清不卡| 国产成人影院久久av| 亚洲欧美激情在线| 高清毛片免费观看视频网站| 亚洲精品国产区一区二| 可以在线观看的亚洲视频| 999久久久国产精品视频| 国产精品美女特级片免费视频播放器 | 欧美国产精品va在线观看不卡| 搡老妇女老女人老熟妇| 一边摸一边抽搐一进一出视频| 亚洲精品中文字幕一二三四区| 中文字幕精品免费在线观看视频| 大型av网站在线播放| 老熟妇乱子伦视频在线观看| a级毛片在线看网站| 亚洲专区国产一区二区| 国产色视频综合| 麻豆久久精品国产亚洲av| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲专区中文字幕在线| 亚洲人成伊人成综合网2020| 欧美日本亚洲视频在线播放| 少妇的丰满在线观看| 啦啦啦免费观看视频1| 91麻豆av在线| 亚洲第一av免费看| 99久久精品国产亚洲精品| 亚洲国产欧美日韩在线播放| 国产成人一区二区三区免费视频网站| 国产野战对白在线观看| 校园春色视频在线观看| 少妇被粗大的猛进出69影院| 欧美一级a爱片免费观看看 | 国产精品免费一区二区三区在线| 免费在线观看日本一区| 成人亚洲精品一区在线观看| 欧美绝顶高潮抽搐喷水| 最近最新免费中文字幕在线| 不卡一级毛片| 国内精品久久久久久久电影| 日本欧美视频一区| 国产人伦9x9x在线观看| 欧美日本中文国产一区发布| 欧美黄色淫秽网站| 激情视频va一区二区三区| 制服诱惑二区| 亚洲精品在线观看二区| 欧美成人免费av一区二区三区| 亚洲精品久久成人aⅴ小说| 天堂影院成人在线观看| 啦啦啦 在线观看视频| 亚洲国产欧美网| 男女之事视频高清在线观看| 侵犯人妻中文字幕一二三四区| 午夜免费激情av| 少妇粗大呻吟视频| 欧美黄色淫秽网站| 日韩欧美一区二区三区在线观看| 亚洲在线自拍视频| 淫秽高清视频在线观看| 欧美日韩瑟瑟在线播放| 这个男人来自地球电影免费观看| 身体一侧抽搐| 97超级碰碰碰精品色视频在线观看| 免费少妇av软件| 久久久久精品国产欧美久久久| 亚洲精品美女久久久久99蜜臀| 看黄色毛片网站| 久久人妻熟女aⅴ| 九色亚洲精品在线播放| 久9热在线精品视频| 人人妻,人人澡人人爽秒播| 久久久久国产一级毛片高清牌| av视频免费观看在线观看| 我的亚洲天堂| 国产精品精品国产色婷婷| 在线十欧美十亚洲十日本专区| 一级,二级,三级黄色视频| 欧美日韩黄片免| 精品无人区乱码1区二区| 欧美av亚洲av综合av国产av| 亚洲欧洲精品一区二区精品久久久| www.999成人在线观看| 亚洲国产欧美网| 亚洲一码二码三码区别大吗| 国产午夜福利久久久久久| 三级毛片av免费| 亚洲国产精品999在线| 午夜成年电影在线免费观看| 国产精品二区激情视频| 欧美日本视频| av网站免费在线观看视频| 在线观看www视频免费| 日本一区二区免费在线视频| 国产又爽黄色视频| 制服人妻中文乱码| 色综合站精品国产| 国产精品乱码一区二三区的特点 | 嫩草影视91久久| 99久久综合精品五月天人人| 亚洲第一青青草原| 中文字幕高清在线视频| 国产成人精品久久二区二区91| 一区二区三区激情视频| x7x7x7水蜜桃| 日本 av在线| 50天的宝宝边吃奶边哭怎么回事| 亚洲九九香蕉| 一级黄色大片毛片| 美女 人体艺术 gogo| 国产一区在线观看成人免费| 国产成人av激情在线播放| 国产av精品麻豆| 亚洲欧美精品综合一区二区三区| 一级作爱视频免费观看| 亚洲欧美日韩无卡精品| xxx96com| av在线天堂中文字幕| 黄色成人免费大全| 99久久国产精品久久久| avwww免费| 免费在线观看日本一区| 免费观看人在逋| 久久香蕉精品热| 18禁黄网站禁片午夜丰满| 午夜福利欧美成人| 久久青草综合色| 91av网站免费观看| 久久精品91蜜桃| 桃色一区二区三区在线观看| 国产一卡二卡三卡精品| 真人做人爱边吃奶动态| 久久亚洲真实| 天天躁夜夜躁狠狠躁躁| 午夜成年电影在线免费观看| 美女高潮到喷水免费观看| 不卡av一区二区三区| 成人手机av| 亚洲成人精品中文字幕电影| 亚洲情色 制服丝袜| 国内毛片毛片毛片毛片毛片| 在线观看免费视频日本深夜| aaaaa片日本免费| 色综合亚洲欧美另类图片| 热re99久久国产66热| 最近最新免费中文字幕在线| 黑人巨大精品欧美一区二区蜜桃| 法律面前人人平等表现在哪些方面| 国产黄a三级三级三级人| 亚洲国产精品合色在线| 午夜精品久久久久久毛片777| 国产又爽黄色视频| 精品国内亚洲2022精品成人| 搡老妇女老女人老熟妇| 国产成人欧美| 成人三级做爰电影| 久久香蕉国产精品| 97人妻天天添夜夜摸| 18禁黄网站禁片午夜丰满| 亚洲精品久久成人aⅴ小说| 久久国产亚洲av麻豆专区| av网站免费在线观看视频| 男人舔女人的私密视频| √禁漫天堂资源中文www| 欧美日韩中文字幕国产精品一区二区三区 | 搡老岳熟女国产| 啦啦啦观看免费观看视频高清 | 国产精品亚洲一级av第二区| 欧美色欧美亚洲另类二区 | 国产成人精品久久二区二区免费| av超薄肉色丝袜交足视频| 亚洲av成人不卡在线观看播放网| 久久久久精品国产欧美久久久| 久久久精品欧美日韩精品| 麻豆成人av在线观看| 久久这里只有精品19| 欧美午夜高清在线| 在线观看舔阴道视频| 在线观看免费日韩欧美大片| 美女大奶头视频| 看免费av毛片| 老鸭窝网址在线观看| 久久人人97超碰香蕉20202| www.熟女人妻精品国产| 黄色视频,在线免费观看| 在线观看www视频免费| 99久久综合精品五月天人人| 久久香蕉国产精品| 久久国产精品影院| 国产区一区二久久| 免费在线观看完整版高清| 欧美 亚洲 国产 日韩一| 69精品国产乱码久久久| 一区二区日韩欧美中文字幕| 亚洲av成人一区二区三| 久久性视频一级片| 日本欧美视频一区| 看片在线看免费视频| 在线观看午夜福利视频| 波多野结衣av一区二区av| 怎么达到女性高潮| 久久久久久久久久久久大奶| 亚洲精品在线观看二区| 9色porny在线观看| 99国产精品一区二区蜜桃av| 黑人巨大精品欧美一区二区蜜桃| 国产日韩一区二区三区精品不卡| 日本a在线网址| 久久人人97超碰香蕉20202| 制服人妻中文乱码| 亚洲精品中文字幕一二三四区| 999久久久国产精品视频| 日本在线视频免费播放| 国产精品一区二区三区四区久久 | 久久影院123| 亚洲欧美一区二区三区黑人| 男人舔女人的私密视频| 91大片在线观看| 搡老妇女老女人老熟妇| 午夜福利高清视频| 制服人妻中文乱码| 亚洲一码二码三码区别大吗| 国产一区二区三区视频了| 婷婷丁香在线五月| 91成年电影在线观看|