999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

渾太河魚(yú)類(lèi)群落多樣性及生境適宜性量化分析

2019-11-28 10:50:54盧奧然殷旭旺徐宗學(xué)
中國(guó)環(huán)境科學(xué) 2019年11期
關(guān)鍵詞:物種環(huán)境

吳 丹,盧奧然,田 爽,殷旭旺*,張 遠(yuǎn),徐宗學(xué)

渾太河魚(yú)類(lèi)群落多樣性及生境適宜性量化分析

吳 丹1,盧奧然1,田 爽1,殷旭旺1*,張 遠(yuǎn)2,徐宗學(xué)3

(1.大連海洋大學(xué)水產(chǎn)與生命學(xué)院,遼寧省水生生物學(xué)重點(diǎn)實(shí)驗(yàn)室,遼寧 大連 116023;2.中國(guó)環(huán)境科學(xué)研究院,環(huán)境基準(zhǔn)與風(fēng)險(xiǎn)評(píng)估國(guó)家重點(diǎn)實(shí)驗(yàn)室,流域水生態(tài)保護(hù)技術(shù)研究室,北京 100012;3.北京師范大學(xué)水科學(xué)研究院,北京 100875)

以遼寧渾太河為范例,評(píng)價(jià)物種、多樣性水平和功能群組成3種群落結(jié)構(gòu)特征下, 魚(yú)類(lèi)群落與棲息地環(huán)境因子的定量響應(yīng)關(guān)系.結(jié)果表明:流速、BOD5、氨氮、電導(dǎo)率和底質(zhì)指數(shù)是顯著影響渾太河魚(yú)類(lèi)群落空間分布的棲境因子;鯽等耐污種的氨氮和電導(dǎo)率的最適值和正響應(yīng)閾值較高,洛氏鱥和東北七鰓鰻等適應(yīng)較高的流速和底質(zhì)指數(shù),東北七鰓鰻的最適流速為0.59m/s;多樣性水平的棲境因子最適值和閾值隨香農(nóng)多樣性指數(shù)的增加呈現(xiàn)先上升后下降的趨勢(shì),(2-3)區(qū)間的電導(dǎo)率最適值最大,為378.07μS/cm;雜食性功能群魚(yú)類(lèi)對(duì)BOD5、氨氮和電導(dǎo)率有較高的最適值和正響應(yīng)閾值,其最適值分別為1.20mg/L、0.63mg/L和383.37μS/cm.定量分析魚(yú)類(lèi)群落與棲境因子的關(guān)系,為開(kāi)展重要魚(yú)類(lèi)物種保護(hù)、生物多樣性恢復(fù)以及關(guān)鍵棲息地的生境修復(fù)等流域生態(tài)管理決策提供基礎(chǔ)數(shù)據(jù)和科學(xué)依據(jù).

魚(yú)類(lèi)群落;棲息地環(huán)境;最適值;生態(tài)閾值;渾太河

棲息地環(huán)境在整個(gè)河流生態(tài)系統(tǒng)中發(fā)揮著重要作用[1],河流棲息地特征影響當(dāng)?shù)厣锶郝涞慕Y(jié)構(gòu),與生物群落多樣性緊密相關(guān)[2].不同的棲息地環(huán)境因子均與水生生物群落有較強(qiáng)的響應(yīng)關(guān)系[3-5].水體電導(dǎo)率可以綜合底質(zhì)類(lèi)型、溫度、海拔和流速等流域特征,影響河流中硅藻群落的分布[6];水溫、流量、污染物含量、懸浮物含量等可影響?hù)~(yú)類(lèi)物種豐度及組成[7].此外,不同的魚(yú)類(lèi)群落所適宜的棲息地環(huán)境也不盡相同,敏感種洛氏鱥、珍稀瀕危物種東北七鰓鰻和雷氏七鰓鰻、植食性功能群及昆蟲(chóng)食性功能群魚(yú)類(lèi)等主要棲息于底質(zhì)類(lèi)型以石塊為主的清潔水質(zhì)和復(fù)雜生境[5],其河流生境多樣性較高;但鯽等耐污種以及沙棲功能群魚(yú)類(lèi)則多分布于渾濁水體且底質(zhì)為泥沙的單一生境水體中[5],河流生境多樣性也較低.

魚(yú)類(lèi)占據(jù)著河流生態(tài)系統(tǒng)食物鏈的頂端,是水生態(tài)系統(tǒng)監(jiān)測(cè)的重要指示生物[8].國(guó)內(nèi)外關(guān)于魚(yú)類(lèi)群落多樣性與棲息地環(huán)境因子關(guān)系的研究,多是通過(guò)多元統(tǒng)計(jì)分析,定性判斷兩者之間的響應(yīng)程度[9-11];李捷等[12]根據(jù)聚類(lèi)分析和非度量多維標(biāo)度排序的分析方法將連江魚(yú)類(lèi)劃分為3個(gè)類(lèi)群,與歷史數(shù)據(jù)對(duì)比后發(fā)現(xiàn)連江魚(yú)類(lèi)群落多樣性的巨大變化;而有關(guān)魚(yú)類(lèi)群落與棲息地環(huán)境因子的定量分析較少報(bào)道.本研究以遼寧省渾太河為研究范例,選取魚(yú)類(lèi)群落的物種組成、多樣性水平和攝食功能群類(lèi)型為評(píng)價(jià)指標(biāo),借助加權(quán)平均回歸分析和指示物種分析的方法,定量得出它們與棲息地環(huán)境因子最適值和閾值[13-15],為在本流域進(jìn)行重要魚(yú)類(lèi)物種(如東北七鰓鰻和雷氏七鰓鰻)的資源保護(hù)、闡明魚(yú)類(lèi)群落多樣性的維持機(jī)制以及關(guān)鍵棲息地的生態(tài)恢復(fù)等流域生態(tài)管理決策提供基礎(chǔ)數(shù)據(jù)和科學(xué)依據(jù).

1 材料與方法

1.1 研究區(qū)域及采樣點(diǎn)位設(shè)置

渾太河位于遼寧省境內(nèi)(40.45°N~42.30°N, 122.00°E~125.30°E),面積為2.73′104km2[16-18].本文在渾太河開(kāi)展3次獨(dú)立調(diào)查研究(2018年4月、2016年9月和2012年5月),獲得162個(gè)站位的魚(yú)類(lèi)群落、水環(huán)境質(zhì)量、水文特征和物理生境數(shù)據(jù)(圖1).

圖1 渾太河采樣點(diǎn)分布

1.2 魚(yú)類(lèi)樣品采集及棲境因子測(cè)定方法

每個(gè)樣點(diǎn)選用電魚(yú)法和掛網(wǎng)法配合完成樣品采集[5],現(xiàn)場(chǎng)對(duì)魚(yú)類(lèi)樣品進(jìn)行物種鑒定[19-20].其中,東北七鰓鰻和雷氏七鰓鰻屬于圓口綱[21],在本文中歸屬劃分為魚(yú)類(lèi);現(xiàn)場(chǎng)進(jìn)行流域棲息地特征的調(diào)查估計(jì),包括經(jīng)緯度、水深(Depth)、流速(Velocity)、水溫(WT)、電導(dǎo)率(COND)、溶解氧(DO)、底質(zhì)含沙量和底質(zhì)指數(shù)特征記錄及計(jì)算公式見(jiàn)相關(guān)研究[5];化學(xué)需氧量(COD)、葉綠素a、5日生化需氧量(BOD5)、高錳酸鹽指數(shù)(CODMn)、氨氮(NH4+-N)、總氮(TN)和總磷(TP)的測(cè)定方法按照《水和廢水監(jiān)測(cè)方法》第四版執(zhí)行[22].

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

物種層面,排除僅在3個(gè)及以下樣點(diǎn)中出現(xiàn)的物種[14]和環(huán)境因子數(shù)據(jù)出現(xiàn)少于3次的物種;計(jì)算魚(yú)類(lèi)群落的香農(nóng)-維納多樣性指數(shù)并根據(jù)結(jié)果將采樣站位進(jìn)行多樣性水平劃分,不同的分組反應(yīng)不同的多樣性水平;功能群組成根據(jù)魚(yú)類(lèi)的營(yíng)養(yǎng)和食性,分為5個(gè)攝食功能群:雜食性功能群、魚(yú)食性功能群、植食性功能群、昆蟲(chóng)食性功能群和底棲動(dòng)物食性功能群[5].

使用Biodiversity Professional 2.0軟件計(jì)算出香農(nóng)-維納多樣性指數(shù);在CANOCO4.5軟件包上對(duì)魚(yú)類(lèi)調(diào)查數(shù)據(jù)進(jìn)行去趨勢(shì)對(duì)應(yīng)分析(DCA)和典范對(duì)應(yīng)分析(CCA)[23];使用加權(quán)平均回歸分析(WA)計(jì)算不同群落水平棲境因子的最適值[24],指示物種分析法(TITAN)確定棲境因子的閾值[15,25];WA和TITAN均在R軟件中完成,WA計(jì)算見(jiàn)式(1).

式中:x是采樣點(diǎn)中的環(huán)境變量值;y是屬種在采樣點(diǎn)種的百分含量;是環(huán)境資料中的總采樣點(diǎn)數(shù)[26].

2 結(jié)果與分析

2.1 魚(yú)類(lèi)群落多樣性結(jié)構(gòu)及其與棲境因子關(guān)系

渾太河共采集到魚(yú)類(lèi)61種(表1),分8目15科共計(jì)24130尾.優(yōu)勢(shì)種為洛氏鱥()、棒花魚(yú)()、寬鰭鱲(),分別占魚(yú)類(lèi)總個(gè)體數(shù)的32.3%、8.3%、8.0%;洛氏鱥為典型的敏感種;另外,有東北七鰓鰻()、雷氏七鰓鰻()和小杜父魚(yú)()等珍稀瀕危魚(yú)類(lèi),分別占魚(yú)類(lèi)總個(gè)體數(shù)的0.3%、0.18%和0.0097%.

表1 渾太河魚(yú)類(lèi)物種組成

由香農(nóng)-維納多樣性指數(shù)得到魚(yú)類(lèi)群落多樣性水平范圍為0.00~3.63(表2),據(jù)此將魚(yú)類(lèi)群落劃分為(0~1)、(1~2)、(2~3)、和(3~4)4個(gè)多樣性水平區(qū)間,其中(2~3)區(qū)間的點(diǎn)位最多,有64個(gè), (3~4]區(qū)間的點(diǎn)位最少,只有19個(gè);從營(yíng)養(yǎng)結(jié)構(gòu)各類(lèi)功能群所占比例來(lái)看:雜食性功能群和植食性功能群在渾太河所占比例最高,其中雜食性功能群魚(yú)類(lèi)分布最廣泛,分布于渾太河的154個(gè)點(diǎn)位,平均物種數(shù)量為85種;底棲動(dòng)物食性功能群在渾太河點(diǎn)位所占比例最少,僅占34個(gè),平均物種數(shù)量為9種.

表2 渾太河魚(yú)類(lèi)物種棲息地環(huán)境因子的最適值分析

圖2 渾太河采樣站位與棲息地環(huán)境因子的典范對(duì)應(yīng)分析

如圖2所示,流速、電導(dǎo)率、BOD5和氨氮是對(duì)渾太河魚(yú)類(lèi)群落起顯著影響的水文和理化因子;相關(guān)研究中發(fā)現(xiàn)底質(zhì)類(lèi)型同樣顯著影響?hù)~(yú)類(lèi)群落的組成[27],因此在本文中也作為一種關(guān)鍵環(huán)境因子進(jìn)行分析.

2.2 魚(yú)類(lèi)物種棲境因子的最適值與閾值

表3 渾太河魚(yú)類(lèi)物種對(duì)棲息地環(huán)境因子的閾值

續(xù)表3

注:環(huán)境閾值: 各物種的響應(yīng)閾值;頻數(shù):每個(gè)分類(lèi)單元的非零豐度值的數(shù)目即物種出現(xiàn)的頻率;響應(yīng)方向:其中“-”表示負(fù)響應(yīng),“+”表示正響應(yīng).

圖3 渾太河魚(yú)類(lèi)物種負(fù)響應(yīng)種(z-)和正響應(yīng)種(z+)指示總分對(duì)棲息地環(huán)境因子的突變點(diǎn)的響應(yīng)曲線

TITAN結(jié)果見(jiàn)表4,渾太河魚(yú)類(lèi)正(z+)負(fù)(z-)響應(yīng)物種指示總分對(duì)棲息地環(huán)境因子突變點(diǎn)的響應(yīng)曲線見(jiàn)圖3(a~e).洛氏鱥在流速、氨氮、BOD5、電導(dǎo)率和底質(zhì)指數(shù)5種棲息地環(huán)境因子下的出現(xiàn)頻率均最高,是流速和底質(zhì)指數(shù)的正響應(yīng)物種,而黃帶克麗鰕虎魚(yú)對(duì)流速和底質(zhì)指數(shù)的出現(xiàn)頻率最低.

表4 魚(yú)類(lèi)群落多樣性水平棲息地環(huán)境因子的最適值

2.3 魚(yú)類(lèi)群落多樣性棲境因子的最適值與閾值

如表2所示,中等多樣性水平(2~3)區(qū)間的棲境因子最適值較高,(3~4)區(qū)間對(duì)流速、氨氮、電導(dǎo)率和底質(zhì)指數(shù)的最適值相比最低;閾值結(jié)果如表5所示,正(z+)負(fù)(z-)響應(yīng)物種指示總分對(duì)棲息地環(huán)境因子的突變點(diǎn)響應(yīng)曲線見(jiàn)圖4(f~j).流速和底質(zhì)指數(shù)最小正響應(yīng)值均出現(xiàn)在(2~3)區(qū)間;BOD5的負(fù)響應(yīng)最大值發(fā)生在(0~1)區(qū)間;(3~4)區(qū)間對(duì)氨氮和電導(dǎo)率的負(fù)響應(yīng)閾值最大.

表5 魚(yú)類(lèi)群落多樣性水平棲息地環(huán)境因子的閾值

注:純度為重抽樣中突變點(diǎn)的響應(yīng)方向與觀察到的響應(yīng)方向相匹配的比例;可靠性為在0.05 的顯著性水平下自舉抽樣的可靠概率.

圖4 渾太河魚(yú)類(lèi)多樣性水平負(fù)響應(yīng)種(z-)和正響應(yīng)種(z+)指示總分對(duì)棲息地環(huán)境因子的突變點(diǎn)響應(yīng)曲線

2.4 魚(yú)類(lèi)群落功能群棲境因子的最適值與閾值

表6 魚(yú)類(lèi)群落功能群棲息地環(huán)境因子的最適值

表7 魚(yú)類(lèi)群落功能群棲息地環(huán)境因子的閾值

續(xù)表7

如表6所示,雜食性魚(yú)類(lèi)對(duì)棲境因子的最適值偏高,昆蟲(chóng)食性功能群對(duì)5種環(huán)境因子的最適值均處于中等水平;閾值結(jié)果如表7所示,正(z+)負(fù)(z-)響應(yīng)種指示總分對(duì)棲息地環(huán)境因子的突變點(diǎn)響應(yīng)曲線見(jiàn)圖5(k-o).底棲動(dòng)物食性、昆蟲(chóng)食性、魚(yú)食性和植食性功能群對(duì)氨氮和電導(dǎo)率均呈現(xiàn)負(fù)響應(yīng),昆蟲(chóng)食性功能群魚(yú)類(lèi)對(duì)流速表現(xiàn)出最小的正響應(yīng)閾值,魚(yú)食性功能群魚(yú)類(lèi)對(duì)底質(zhì)指數(shù)有最小的正響應(yīng)閾值.

圖5 渾太河魚(yú)類(lèi)群落功能群負(fù)響應(yīng)種(z-)和正響應(yīng)種(z+)指示總分對(duì)棲息地環(huán)境因子的突變點(diǎn)響應(yīng)曲線

3 討論

3.1 魚(yú)類(lèi)群落結(jié)構(gòu)特征及棲境因子

渾太河魚(yú)類(lèi)群落的物種組成和多樣性水平等與相關(guān)研究結(jié)果一致[5],表現(xiàn)為上游生物多樣性水平較高,敏感種和魚(yú)食性功能群占優(yōu)勢(shì),中下游耐污種和雜食性功能群分布較多,多樣性水平較低.復(fù)雜的棲息地環(huán)境影響?hù)~(yú)類(lèi)的多樣性空間分布[28];研究表明,渾太河上游森林覆蓋率高,河流底質(zhì)復(fù)雜,自然生境完整,因而魚(yú)類(lèi)多樣性水平較高[29];而中下游河岸兩側(cè)底質(zhì)以泥沙、淤泥為主,多為耕作用地或建筑用地,堤岸受到不同程度的侵蝕,其生物多樣性也較低[5].李艷利[30]對(duì)渾太河不同尺度下人類(lèi)活動(dòng)對(duì)魚(yú)類(lèi)群落特征的影響研究中發(fā)現(xiàn),電導(dǎo)率、TDS、氨氮、總氮等是顯著影響?hù)~(yú)類(lèi)群落分布的環(huán)境因子,李麗娟[5]的相關(guān)研究指出底質(zhì)指數(shù)很大程度上也影響了魚(yú)類(lèi)群落功能群的分布;本研究得到結(jié)果類(lèi)似,結(jié)合實(shí)際調(diào)查情況最終得到流速、氨氮、BOD5、電導(dǎo)率和底質(zhì)指數(shù)是顯著影響渾太河魚(yú)類(lèi)群落的棲息地環(huán)境因子.

研究表明,魚(yú)類(lèi)群落受環(huán)境因子影響的結(jié)果會(huì)因流域、地理位置和人類(lèi)干擾因素等的不同而有所差異[31];灤河流域魚(yú)類(lèi)群落空間異質(zhì)性受到TDS、濁度、pH值、水溫、底質(zhì)與堤岸穩(wěn)定性的顯著影響[32];新疆額敏河干流影響?hù)~(yú)類(lèi)群落變異的主要環(huán)境因素是pH值、DO、溫度與電阻率[23];此外,因研究尺度不同,影響同一流域魚(yú)類(lèi)群落的環(huán)境因子也不盡相同[30];太子河魚(yú)類(lèi)功能群結(jié)構(gòu)與多樣性對(duì)土地利用類(lèi)型的響應(yīng)分析中得到底質(zhì)指數(shù)、DO、電導(dǎo)率等是影響?hù)~(yú)類(lèi)群落結(jié)構(gòu)的環(huán)境因子[5],而丁森等[33]對(duì)太子河流域影響?hù)~(yú)類(lèi)空間分布的研究中得到海拔、河流等級(jí)、速度與深度結(jié)合等級(jí)、底質(zhì)等參數(shù)是重要環(huán)境要素.因此,不同模式下的魚(yú)類(lèi)群落結(jié)構(gòu)與棲息地環(huán)境因子的研究,其結(jié)果也有所不同[31].

3.2 魚(yú)類(lèi)群落對(duì)棲境因子的最適值分析

渾太河魚(yú)類(lèi)物種的棲境因子最適值差異較大,魚(yú)類(lèi)自身特征的差異使其對(duì)棲息地環(huán)境的適應(yīng)狀況大不相同[10,31].如七鰓鰻營(yíng)埋伏和吸著生活,喜歡棲居于有巖石和砂礫的清潔河流環(huán)境,夜晚積極的游動(dòng)于石隙和沙礫間來(lái)覓食[34],洛氏鱥同樣喜棲于底質(zhì)多頁(yè)巖、礫石或山澗多石縫的清潔水體[35];而廣適性耐污種魚(yú)類(lèi)鯽,可適應(yīng)不同深淺的水體,在低氧、酸、堿等污染環(huán)境也可生存[36],另外,鯽產(chǎn)出的粘性卵常附著在草上,喜在靜水中繁殖的特性[37],也積極的解釋了其適應(yīng)低流速水環(huán)境的原因;沙塘鱧對(duì)不同環(huán)境因子的最適值也與其底棲穴居的生活習(xí)性相適應(yīng)[38].

魚(yú)類(lèi)多樣性水平具有明顯的空間異質(zhì)性.河流上游生境污染較小,水體清潔、棲息地多樣性和復(fù)雜性一般最大[39],因此渾太河中上游有較高生物多樣性水平,而下游由于人類(lèi)活動(dòng)干擾較多,引起的入河污染物和顆粒物增加,造成水質(zhì)渾濁,其生物多樣性相對(duì)較低;另有研究表明,在流速較低和泥沙質(zhì)的生態(tài)環(huán)境中,更有利于小型水生生物的附著、生存、棲息和繁衍,魚(yú)類(lèi)群落因此得到更充足的食物[40],具有更高的生物多樣性;同時(shí),中度干擾假說(shuō)[41]為本文中中等水平環(huán)境因子最適值對(duì)應(yīng)較高魚(yú)類(lèi)群落多樣性水平香農(nóng)-維納指數(shù)區(qū)間的結(jié)果闡明了原因.

魚(yú)類(lèi)群落攝食功能群的分布差異與棲息地環(huán)境中食物來(lái)源的空間變化緊密相關(guān)[42-43],河流底質(zhì)類(lèi)型決定了生物群落棲息地環(huán)境的不同[44];上游水質(zhì)清潔,物種豐富,底質(zhì)以碎石塊為主[45],其生物多樣性水平較高,魚(yú)類(lèi)食物來(lái)源廣泛;而中下游水環(huán)境較差,水質(zhì)離子濃度高,多呈富營(yíng)養(yǎng)狀態(tài),處于下游的攝食功能群多如鯽等耐受性強(qiáng)的魚(yú)類(lèi),因此不同攝食功能群的棲境因子最適值也有所差異.

3.3 魚(yú)類(lèi)群落對(duì)棲境因子的閾值分析

渾太河不同魚(yú)類(lèi)群落結(jié)構(gòu)的棲境因子閾值結(jié)果與最適值結(jié)果相一致;研究表明,魚(yú)類(lèi)因自身形態(tài)、生理學(xué)特征、生活史和對(duì)棲息地環(huán)境適宜程度的差異而對(duì)環(huán)境因子指示出不同的響應(yīng)結(jié)果[34-37];沙塘鱧是渾太河中重要的頂級(jí)捕食者,其底棲穴居適應(yīng)清潔水體的生活習(xí)性與本文中較高底質(zhì)指數(shù)、較低流速、氨氮、BOD5和電導(dǎo)率的結(jié)果相符合[38].

物種多樣性的變化除受群落組成特征和群落演替動(dòng)態(tài)的影響外,還與自然或人為干擾因子及生境因子息息相關(guān)[46];渾太河中下游破壞嚴(yán)重,加劇的人類(lèi)活動(dòng)改變了河流中物種組成,降低了物種多樣性[39].研究表明,水質(zhì)較好的河流環(huán)境中,水體微生物、水生植物、浮游動(dòng)植物的物種多樣性和豐度都要好于環(huán)境較差的水體[46-47];因此,嚴(yán)控過(guò)度捕撈、環(huán)境污染等人類(lèi)活動(dòng)干擾情況的發(fā)生,對(duì)維持河流中的高物種多樣性具有重要的作用.

功能群組成的分布差異與以食物為基礎(chǔ)的時(shí)空變化密切相關(guān)[42-43].上游物種豐度及多樣性較高,攝食小型魚(yú)類(lèi)的魚(yú)食性功能群多居于此,良好的生存空間為其提供了豐富的食物來(lái)源;中下游水體富營(yíng)養(yǎng)化和有機(jī)污染較嚴(yán)重,食物的質(zhì)量與數(shù)量都不能得到滿(mǎn)足,雜食性魚(yú)類(lèi)食物組成廣泛,水中浮游生物、底棲動(dòng)物、動(dòng)植物有機(jī)碎屑等都可被攝食.此外,渾太河上游底質(zhì)類(lèi)型以石塊為主,流速較大,適合食物鏈頂級(jí)捕食者生活[48];中下游土地利用方式復(fù)雜,水體富營(yíng)養(yǎng)化嚴(yán)重[5],鯽等雜食性功能群魚(yú)類(lèi)具有更好的耐受性.

3.4 基于魚(yú)類(lèi)群落生境適宜性分析的渾太河生境修復(fù)對(duì)策

棲息地環(huán)境的不同決定了生物多樣性水平的不同,影響了魚(yú)類(lèi)的食物來(lái)源,進(jìn)而導(dǎo)致攝食功能群分布的差異.根據(jù)TITAN和WA進(jìn)行定量分析,可確定關(guān)鍵魚(yú)類(lèi)物種或高物種多樣性水平的棲境因子耐受值和最適值,預(yù)判關(guān)鍵物種的適宜生境,對(duì)珍稀瀕危物種的生境進(jìn)行人工修復(fù)與調(diào)整;如保護(hù)珍稀瀕危物種雷氏七鰓鰻,當(dāng)流速大于0.60m/s、底質(zhì)指數(shù)大于2.74、BOD5大于2.85mg/L、氨氮小于0.05mg/L且電導(dǎo)率小于251μs/cm時(shí),其棲息地環(huán)境良好,同時(shí)當(dāng)流速、氨氮、BOD5、電導(dǎo)率和底質(zhì)指數(shù)分別為0.40m/s、0.11mg/L、0.41mg/L、210.3μs/cm和1.78時(shí),雷氏七鰓鰻的物種豐度最高,為最佳生存環(huán)境;因此可利用減少對(duì)雷氏七鰓鰻棲息的渾太河北支河岸帶開(kāi)發(fā)利用的方式對(duì)棲息地進(jìn)行修復(fù)調(diào)整.

4 結(jié)語(yǔ)

流速、BOD5、氨氮、電導(dǎo)率和底質(zhì)指數(shù)是顯著影響渾太河魚(yú)類(lèi)群落空間分布的棲境因子;洛氏鱥和東北七鰓鰻等適應(yīng)較高的流速和底質(zhì)指數(shù),東北七鰓鰻的最適流速為0.59m/s;多樣性水平棲境因子的最適值和閾值隨香農(nóng)多樣性指數(shù)的增加呈現(xiàn)先上升后下降的趨勢(shì),(2~3]區(qū)間的電導(dǎo)率最適值最大,為378.07mS/cm;雜食性功能群魚(yú)類(lèi)對(duì)BOD5、氨氮和電導(dǎo)率有較高的最適值和正響應(yīng)閾值,其最適值分別為1.20mg/L、0.63mg/L和383.37μS/cm.

研究魚(yú)類(lèi)群落與棲息地環(huán)境因子的定量關(guān)系,可以相對(duì)準(zhǔn)確的評(píng)估重要魚(yú)類(lèi)物種(如雷氏七鰓鰻)、高物種多樣性和關(guān)鍵功能群類(lèi)型所對(duì)應(yīng)的棲息地環(huán)境特征,為在本流域進(jìn)行重要魚(yú)類(lèi)物種及資源保護(hù)、維持高物種多樣性及關(guān)鍵棲息地的生境修復(fù)等流域生態(tài)管理決策和環(huán)境保護(hù)提供基礎(chǔ)數(shù)據(jù)和科學(xué)依據(jù).

[1] 鄭丙輝,張 遠(yuǎn),李英博.遼河流域河流棲息地評(píng)價(jià)指標(biāo)與評(píng)價(jià)方法研究 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2007,27(6):928-936. Zheng Bing-Hui, Zhang Yuan, Li Ying-Bo. Study on river habitat evaluation index and evaluation method in Liaohe River Basin [J]. Journal of Environmental Science, 2007,27(6):928-936.

[2] Meffe G K, Sheldon A L. The influence of habitat structure on fish assemblage composition in Southeastern Blackwater Streams [J]. American Midland Naturalist, 1988,120(2):225-240.

[3] 殷旭旺,徐宗學(xué),高 欣,等.渭河流域大型底棲動(dòng)物群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 2013,24(1):218-226. Yin Xu-Wang, Xu Zong-Xue, Gao Xin, et al. The community structure of macrofauna in Weihe River and its relationship with environmental factors [J]. Chinese Journal of Applied Ecology, 2013, 24(1):218-226.

[4] 李麗娟,金 文,王博涵,等.太子河河岸帶土地利用類(lèi)型與硅藻群落結(jié)構(gòu)的關(guān)系 [J]. 環(huán)境科學(xué)研究, 2015,28(11):1662-1669. Li Li-Juan, Jin-Wen, Wang Bo-Han, et al. Relationship between Land Use Types within Riparian Zones and Community Structure of Diatom in Taizi River. China [J]. Research of Environmental Sciences, 2015, 28(11):1662-1669.

[5] 李麗娟,張 吉,吳 丹,等.太子河流域魚(yú)類(lèi)功能群結(jié)構(gòu)與多樣性對(duì)土地利用類(lèi)型的響應(yīng) [J]. 生態(tài)學(xué)報(bào), 2017,37(20):6863-6874. Li Li-Juan, Zhang Ji, Wu Dan, et al. Relationships between structure and diversity of fish functional groups and land use in the Taizi River [J]. Acta Ecologica Sinica, 2017,37(20):6863-6874.

[6] Belore M L, Winter J G, Duthie H C. Use of Diatoms and Macroinvertebrates as Bioindicators of Water Quality in Southern Ontario Rivers [J]. Canadian Water Resources Journal, 2002,27(4): 457-484.

[7] Sickle J V, Baker J, Herlihy A, et al. Projecting the biological condition of streams under alternative scenarios of human land use [J]. Ecological Applications, 2004,14(2):368-380.

[8] 高 欣,丁 森,張 遠(yuǎn),等.魚(yú)類(lèi)生物群落對(duì)太子河流域土地利用、河岸帶棲息地質(zhì)量的響應(yīng) [J]. 生態(tài)學(xué)報(bào), 35(21):7198-7206. Gao Xin, Ding Sen, Zhang Yuan, et al. Exploring the relationship among land-use, riparian habitat quality, and biological integrity of a fish community.Acta Ecologica Sinica, 35(21):7198-7206.

[9] 王 華,楊樹(shù)平,房晟忠,等.滇池浮游植物群落特征及與環(huán)境因子的典范對(duì)應(yīng)分析[J]. 中國(guó)環(huán)境科學(xué), 2016,36(2):544-552. Wang Hua, Yang Shu-Ping, Fang Sheng-Zhong, et al. Canonical correspondence analysis of relationship between characteristics of phytoplankton community and environmental factors in Dianchi Lake. China Environmental Sciencece, 2016,36(2):544-552.

[10] Phillips S J, Schapire R E. A maximum entropy approach to species distribution modeling [C]//International Conference on Machine Learning. ACM, 2004:83.

[11] Lluís Brotons, Wilfried Thuiller, Miguel B A, et al. Presence-absence versus presence-only modelling methods for predicting bird habitat suitability [J]. Ecography, 2004,27(4):437-448.

[12] 李 捷,李新輝,賈曉平,等.連江魚(yú)類(lèi)群落多樣性及其與環(huán)境因子的關(guān)系[J]. 生態(tài)學(xué)報(bào), 2011,32(18):5795-5805. Li Jie, Li Xin-Hui, Jia Xiao-Ping, et al. Diversity of fish communities in Lianjiang and its relationship with environmental factors [J]. ACTA Ecologica Sinica, 2011,32(18):5795-5805.

[13] 閔文武,韓 潔,王培培,等.太子河流域硅藻群落與驅(qū)動(dòng)因子的定量關(guān)系 [J]. 環(huán)境科學(xué)研究, 2016,29(5):672-679. Min Wen-Wu, Han Jie, Wang Pei-Pei, et al. Quantitative relationships among the diatom community and driving factors in Taizi River Basin, China [J]. Research of Environmental Sciences, 2016,29(5):672-679.

[14] 湯 婷,任 澤,唐 濤,等.基于附石硅藻的三峽水庫(kù)入庫(kù)支流氮、磷閾值 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 2016,27(8):2670-2678. Tang Ting, Ren Ze, Tang Tao, et al. Total nitrogen and total phosphorus thresholds for epilithic diatom assemblages in inflow tributaries of the Three Gorges Reservoir, China [J]. Chinese Journal of Applied Ecology, 2016,27(8):2670-2678.

[15] Baker M E, King R S. A new method for detecting and interpreting biodiversity and ecological community thresholds [J]. Methods in Ecology & Evolution, 2010,1(1):25-37.

[16] 殷旭旺,渠曉東,李慶南,等.基于著生藻類(lèi)的太子河流域水生態(tài)系統(tǒng)健康評(píng)價(jià) [J]. 生態(tài)學(xué)報(bào), 2012,32(6):1677-1691. Yin Xu-Wang, Qu Xiao-Dong, Li Qing-Nan, et al. Using periphyton assemblages to assess stream conditions of Taizi River Basin [J]. ACTA Ecologica Sinica, 2012,32(6):1677-1691.

[17] 李艷利,李艷粉,徐宗學(xué).影響渾太河流域大型底棲動(dòng)物群落結(jié)構(gòu)的環(huán)境因子分析 [J]. 環(huán)境科學(xué), 2015,(1):94-106. Li Yan-Li, Li Yan-Fen, Xu Zong-Xue. Analysis of environmental factors affecting the community structure of macrobenthic animals in the HunTai River Basin [J]. Environmental Science, 2015,(1):94-106.

[18] 沈洪艷,曹志會(huì),劉軍偉,等.太子河流域大型底棲動(dòng)物功能攝食類(lèi)群與環(huán)境要素的關(guān)系[J]. 中國(guó)環(huán)境科學(xué), 2015,35(2):579-590. Shen Hong-Yan, Cao Zhi-Hui, Liu Jun-Wei, et al. Relationship between functional feeding groups of macroinvertebrates and environmental factors in Taizi River basin. China Environmental Sciencece, 2015,35(2):579-590.

[19] 劉禪馨,秦克靜.遼寧動(dòng)物志:魚(yú)類(lèi) [M]. 沈陽(yáng):遼寧科學(xué)技術(shù)出版社, 1987. Liu Chan-Xin, Qin Ke-Jing. Liaoning Zoology: Fish [M]. Shenyang: Liaoning Science and Technology Press, 1987.

[20] 李思忠,張春光.中國(guó)動(dòng)物志,硬骨魚(yú)綱 [M]. 北京:科學(xué)出版社, 2011. Li Si-Zhong, Zhang Chun-Guang [M]. Beijing: Chinese animal records, hard bone fish class [M]. Science Press, 2011.

[21] 李 軍,董彥嬌,韓英倫,等.東北七鰓鰻和日本七鰓鰻成體形態(tài)性狀對(duì)體質(zhì)量的影響分析 [J]. 海洋學(xué)報(bào), 2017,39(4):61-71. Li Jun, Dong Yan-Jiao, Han Ying-Lun, et al. Effect of morphogenetic characters of the Northeast quail eel and Japanese quail eel on body quality [J]. Marine Journal, 2017,39(4):61-71.

[22] 國(guó)家環(huán)境保護(hù)總局《水和廢水監(jiān)測(cè)分析方法》編委會(huì).水和廢水監(jiān)測(cè)分析方法:第4版[M]. 北京:中國(guó)環(huán)境科學(xué)出版社, 2002. “Monitoring and analysis method of water and wastewater ” editorial board of State Environmental Protection Administration of China ed. Monitoring and analysis method of water and wastewater: Fourth Edition. [M]. Beijing: China Environmental Science Press, 2002.

[23] 劉 鴻,牛建功,劉春池,等.新疆額敏河干流魚(yú)類(lèi)群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系 [J]. 生態(tài)學(xué)雜志, 2017,36(12):3558-3563. Liu Hong, Niu Jian-Gong, Liu Chun-Chi, et al. Fish community structure and its relationship with environmental factors of main stream of Emin River, Xinjiang [J]. Chinese Journal of Ecology, 2017, 36(12):3558-3563.

[24] Birks H J B, Line J M, Juggens S, et al. Diatoms and pH Reconstruction [J]. Royal Society of London Philosophical Transactions, 1990,327(1240):263-278.

[25] King R S, Baker M E. Considerations for analyzing ecological community thresholds in response to anthropogenic environmental gradients [J]. Journal of the North American Benthological Society, 2010,29(3):11.

[26] 鄧培雁.水質(zhì)對(duì)東江流域附生硅藻群落的影響 [J]. 生態(tài)學(xué)報(bào), 2015,35(6):1852-1861. Deng Pei-Yan.The effects of water quality on epilithic diatoms communities of Dongjiang river basin [J]. Acta Ecologica Sinica, 2015,35(6):1852-1861.

[27] Kadye W T, Magadza C H D, Moyo N A G, et al. Stream fish assemblages in relation to environmental factors on a montane plateau (Nyika Plateau, Malawi) [J]. Environmental Biology of Fishes, 2008,83(4):417-428.

[28] 鄧朝陽(yáng),朱 仁,嚴(yán)云志.長(zhǎng)江蕪湖江段魚(yú)類(lèi)多樣性及其群落結(jié)構(gòu)的時(shí)空格局[J]. 淡水漁業(yè), 2016,43(1). DENG Chao-Yang, ZHU Ren, YAN Yun Zhi. Spatial and temporal patterns of fish diversity and assemblage structures in the Wuhu section of the Yangtze River [J]. Freshwater Fisheries, 2016,43(1): 28-31.

[29] 項(xiàng)珍龍.在渾太河流域水生態(tài)系統(tǒng)檢驗(yàn)中度干擾假說(shuō)[D]. 大連:大連海洋大學(xué), 2017. Xiang Zhen-Long. Test of the intermediate disturbance hypothesis in the Huntai River Basin water ecosystem [D]. Dalian: Dalian Ocean University, 2017.

[30] 李艷利,李艷粉.李 科.不同尺度下人類(lèi)活動(dòng)對(duì)渾太河流域魚(yú)類(lèi)和大型底棲動(dòng)物群落特征的影響[J]. 環(huán)境科學(xué)研究, 2016,29(8): 1145-1153. Li Yan-Li, Li Yan-Fen, Li Ke. Relative influence of anthropogenic stressors on fish and macroinvertebrate communities at different scales in Huntai Rivers [J].Research of Environmental Sciences, 2016,29(8): 1145-1153.

[31] 劉 偉,張 遠(yuǎn),高 欣,等.渾河流域魚(yú)類(lèi)群落特征及其與環(huán)境因子的關(guān)系 [J]. 環(huán)境工程技術(shù)學(xué)報(bào), 2016,6(3):266-274. Liu Wei, Zhang Yuan, Gao Xin, et al. Fish community structure and its relationship with environmental variables in Hun River Basin of Liaoning Province [J]. Journal of Environmental Engineering Technology, 2016,6(3):266-274.

[32] 王曉寧,彭世賢,張 亞,等.灤河流域魚(yú)類(lèi)群落結(jié)構(gòu)空間異質(zhì)性與影響因子分析 [J]. 環(huán)境科學(xué)研究, 2018,31(2):273-282. Wang Xiao-Ning, Peng Shi-Xian, Zhang Ya, et al. Analysis of fish community spatial heterogeneity and its driving factors in the Luanhe River Basin [J]. Research of Environmental Sciences, 2018, 31(2):273-282.

[33] 丁 森,張 遠(yuǎn),渠曉東,等.影響太子河流域魚(yú)類(lèi)空間分布的不同尺度環(huán)境因子分析 [J]. 環(huán)境科學(xué), 2012,33(7):2272-2280. Ding Sen, Zhang Yuan, Qu Xiao-Dong, et al. Influence on the spatial distribution of fish in Taizi River Basin by environmental factors at multiple scales [J]. Chinese Journal of Environmental Science, 2012, 33(7):2272-2280.

[34] 馬常夫,于春林.渾江七鰓鰻的初步觀察 [J]. 動(dòng)物學(xué)雜志, 1959, (3):29-31. Ma Chang-Fu, Yu Chun-Lin. Preliminary observation on the seven pelican eels in Pujiang River [J]. Zoology Magazine, 1959,(3):29-31.

[35] 康 鑫,張 遠(yuǎn),張 楠,等.太子河洛氏鱥幼魚(yú)棲息地適宜度評(píng)估 [J]. 生態(tài)毒理學(xué)報(bào), 2011,6(3):310-320. Kang Xin, Zhang Yuan, Zhang Nan, et al. Assessment of habitat suitability ofin Taizi River [J]. Asian Journal of Ecotoxicolog, 2011,6(3):310-320.

[36] Ornellas A B, Coutinho R. Spatial and temporal patterns of distribution and abundance of a tropical fish assemblage in a seasonal Sargassum bed, Cabo Frio Island, Brazil [J]. Journal of Fish Biology, 2010,53(sA): 198-208.

[37] 饒發(fā)祥.鯽魚(yú)種類(lèi)生物學(xué)特性及其特殊的繁殖方式[J]. 北京水產(chǎn), 1996,(Z1):21-23. RAO Fa-Xiang. Biological characteristics of carp species and its special breeding methods [J]. Journal of Beijing fisheries, 1996,(Z1): 21-23.

[38] 胡廷尖,章文敏,周 穎.沙塘鱧的生物學(xué)習(xí)性與繁養(yǎng)技術(shù)的研究 [J]. 北京水產(chǎn), 2004,(5):50-50. Hu Ting-Jian, Zhang Wen-Min, Zhou Ying. Study on biological learning and breeding techniques of Snakehead [J]. Journal of Beijing Fisheries, 2004,(5):50-50.

[39] 何美峰.汀江中上游魚(yú)類(lèi)多樣性及其影響因子 [J]. 福建農(nóng)業(yè)學(xué)報(bào), 2016,31(6):566-574. He Mei-Feng. Factors affecting diversity of fish species in mid-and upper-streams of Tingjiang River [J]. Fujian Agricultural Journal, 2016,31(6):566-574.

[40] 白曉慧,張曉紅,許彭鵬.城市景觀河道不同駁岸界面水生生物多樣性[J]. 南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2008,32(1):111-114. Bai Xiao-Hui, Zhang Xiao-Hong, Xu Peng-Peng. Diversity of aquatic organisms at the interface of different river embankments [J]. Journal of Nanjing Forestry University(Natural Science Edition, 2008,32(1): 111-114.

[41] Sommer U. An experimental test of the intermediate disturbance hypothesis using cultures of marine phytoplankton [J]. Limnology and Oceanography, 1995,40(7):1271-1277.

[42] 張 波,唐啟升,金顯仕.東海高營(yíng)養(yǎng)層次魚(yú)類(lèi)功能群及其主要種類(lèi) [J]. 中國(guó)水產(chǎn)科學(xué), 2007,14(6):939-949. Zhang Bo, Tang Qi-Sheng, Jin Xian-Shi.Functional groups of fish assemblages and their major species at high trophic level in the East China Sea [J].Journal of Fishery Sciences of China, 2007,14(6): 939-949.

[43] 韓 潔,張 杰,宋蒙蒙,等.濟(jì)南市水域魚(yú)類(lèi)功能群結(jié)構(gòu)及其季節(jié)性變化 [J]. 環(huán)境科學(xué)研究, 31(9). Han Jie, Zhang Jie, Song Meng-Meng, et al. Structure and seasonal changes of fish functional groups in Jinan Rgion [J]. Research of Environmental Sciences,2018,(9):1537-1544.

[44] 張 浩,丁 森,張 遠(yuǎn),等.遼河流域類(lèi)生物完整性指數(shù)評(píng)價(jià)及與環(huán)境因子的關(guān)系 [J]. 湖泊科學(xué), 2015,27(5):829-839. Zhagn Hao, Ding Sen, Zhang Yuan. et al. Assessment of the fish index of biotic integrity and its relationship with environmental factors in the Xiliao River Basin [J]. Journal of Lake Sciences, 2015,27(5):829-839.

[45] Beugly J, Pyron M. Temporal and spatial variation in the long-term functional organization of fish assemblages in a large river [J]. Hydrobiologia, 2010,654(1):215-226.

[46] 夏智勇,王 力,齊代華.水生植物群落物種多樣性與水環(huán)境相關(guān)性研究 [J]. 環(huán)境保護(hù), 2011,(5):39-41.Xia Zhi-Yong, Wang Li, Qi Dai-Hua. Relationship between species diversity and water environment of aquatic plant communities [J]. Environmental Protection, 2011,(5):39-41.

[47] 孫寓姣,陳 程,丁愛(ài)中,等.官?gòu)d水庫(kù)水質(zhì)特征及水體微生物多樣性的響應(yīng) [J]. 中國(guó)環(huán)境科學(xué), 2015,35(5):1547-1553. Sun Yu-Jiao, Chen Cheng, Ding Ai-Zhong. et al. The corresponding of microbial diversity on water quality and environmental variables of Guanting Reservoir [J]. China Environmental Sciencece, 2015,35(5): 1547-1553.

[48] 王博涵,吳 丹,張 吉,等.濟(jì)南河流大型底棲動(dòng)物攝食功能群多樣性及時(shí)空動(dòng)態(tài) [J]. 生態(tài)學(xué)報(bào), 2017,(21):115-126. Wang Bo-Han, Wu Dan, Zhang Ji, et al. Diversity and temporal- spatial dynamics of macroinvertebrate functional feeding groups in the rivers of the Jinan Region [J]. Acta Ecologica Sinica, 2017,(21): 115-126.

Quantitative analysis of the diversity and habitat suitability of the Huntai River fish community.

WU dan1, LU Ao-ran1, TIAN Shuang1, YIN Xu-wang1*, ZHANG Yuan2, XU Zong-xue3

(1.Liaoning Provincial Key Laboratory for Hydrobiology, College of Fisheries and Life Science, Dalian Ocean University, Dalian 116023, China;2.State Key Laboratory of Environmental Benchmarking and Risk Assessment, Chinese Academy of Environmental Sciences, Laboratory of Watershed Ecological Protection Technology, Beijing 100012, China;3.Institute of Water Science, Beijing Normal University, Beijing 100875, China)., 2019,39(11):4875~4885

In the Huntai River, Liaoning Province, the optimum and threshold values were calculated for several environmental factors for habitats of different fish community structure characteristics. Three kinds of community organization characteristics were considered: the level of community species diversity, the Shannon diversity index, and functional groups. The results showed that velocity, BOD5, NH4+-N, conductivity, and IOS were the environmental factors that significantly influenced the spatial distribution of fish community structure in the Huntai River basin. The optimum value and positive response threshold of NH4+-N and conductivity were higher in pollution resistant species, such asand.andwere adapted to a higher velocity and IOS, with 0.59 m/s being the optimal velocity forFor the level of fish community diversity, the optimum value and threshold for the environmental factors of the habitat first increased and then decreased, in accordance with an increase in the Shannon diversity index. For example, the maximum conductivity range of (2~3) was 378.07μS/cm. In the middle and lower reaches of the Huntai River, fish in the omnivorous functional group had higher optimum values and positive response thresholds for BOD5, NH4+-N, and conductivity. The optimum values were 1.20 mg/L, 0.63 mg/L, and 383.37 μS/cm, respectively. This quantitative analysis of the community and environmental factors will assist in protecting resources for important fish species and conserving fish biodiversity. This research is necessary to inform management decisions regarding the restoration of key habitats, provide basic data, and provide a scientific basis for environmental protection.

fish communities;habitat environment;optimum value;ecological threshold;Huntai River

X176

A

1000-6923(2019)11-4875-11

吳 丹(1993-)遼寧營(yíng)口人,大連海洋大學(xué)碩士研究生,主要從事河流生態(tài)學(xué)方面的研究.

2019-04-26

遼寧省優(yōu)秀人才支持計(jì)劃(LR2015009);國(guó)家自然科學(xué)基金項(xiàng)目(41977193);遼寧省“興遼英才計(jì)劃”項(xiàng)目(XLYC1807228)

* 責(zé)任作者, 教授, yinxuwang@dlou.edu.cn

猜你喜歡
物種環(huán)境
物種大偵探
物種大偵探
吃光入侵物種真的是解決之道嗎?
長(zhǎng)期鍛煉創(chuàng)造體內(nèi)抑癌環(huán)境
一種用于自主學(xué)習(xí)的虛擬仿真環(huán)境
孕期遠(yuǎn)離容易致畸的環(huán)境
不能改變環(huán)境,那就改變心境
回首2018,這些新物種值得關(guān)注
環(huán)境
孕期遠(yuǎn)離容易致畸的環(huán)境
主站蜘蛛池模板: 女人av社区男人的天堂| 国产区福利小视频在线观看尤物| 国产草草影院18成年视频| 99久久性生片| 成人毛片免费在线观看| 国产亚洲精品91| 欧洲高清无码在线| 亚洲区第一页| 欧美一区二区三区国产精品| 国产清纯在线一区二区WWW| 青青青国产在线播放| a亚洲视频| 女同国产精品一区二区| 亚洲精品男人天堂| 国产亚洲成AⅤ人片在线观看| 亚洲日韩AV无码一区二区三区人 | 视频一本大道香蕉久在线播放| 日韩精品一区二区三区视频免费看| 极品尤物av美乳在线观看| 91免费片| 看国产一级毛片| 欧美日韩免费在线视频| 97se亚洲综合在线韩国专区福利| 成年人福利视频| 亚洲第一成网站| 亚洲区欧美区| 亚洲狼网站狼狼鲁亚洲下载| 天天色天天综合网| 91网址在线播放| 国产成人1024精品| 国产精品免费电影| 久久久久夜色精品波多野结衣| 久久精品无码一区二区国产区| 久久久久免费看成人影片| 欧美激情成人网| 亚洲欧美人成人让影院| 国产精品欧美在线观看| 国产毛片高清一级国语| 五月天综合婷婷| 亚洲AV人人澡人人双人| 99国产在线视频| 亚洲精品无码抽插日韩| 99热最新网址| 国产性生交xxxxx免费| 在线视频97| AV片亚洲国产男人的天堂| 日韩成人高清无码| 99re66精品视频在线观看| 免费可以看的无遮挡av无码| 亚洲AV无码久久精品色欲| 亚洲精品中文字幕午夜| 在线网站18禁| 免费国产一级 片内射老| 麻豆精品在线| 第一区免费在线观看| 成年片色大黄全免费网站久久| 黄色污网站在线观看| 久久国产精品国产自线拍| 91免费观看视频| 71pao成人国产永久免费视频| 亚洲a免费| 中文字幕佐山爱一区二区免费| 亚洲五月激情网| 97视频精品全国在线观看| 亚洲无线一二三四区男男| 亚洲一级无毛片无码在线免费视频 | 国产黄色爱视频| 亚洲精品国产自在现线最新| 国产97公开成人免费视频| 久热re国产手机在线观看| 日韩 欧美 国产 精品 综合| 免费大黄网站在线观看| 国产精品一区在线观看你懂的| 激情综合网址| 亚洲天堂网2014| 99免费视频观看| 天堂在线视频精品| 亚洲天堂网2014| 四虎永久免费地址| 国产女人爽到高潮的免费视频 | 亚洲区第一页| 国产喷水视频|