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

    渾河流域2010—2014年的魚類群落和水生態(tài)健康變化分析

    2016-11-12 06:40:47賈曉波馬淑芹劉思思
    水生生物學(xué)報(bào) 2016年5期
    關(guān)鍵詞:渾河完整性魚類

    劉 偉 張 遠(yuǎn) 高 欣 賈曉波 馬淑芹 劉思思

    (1. 中國環(huán)境科學(xué)研究院流域水生態(tài)保護(hù)技術(shù)研究室,北京 100012; 2. 環(huán)境標(biāo)準(zhǔn)與風(fēng)險(xiǎn)評估國家重點(diǎn)實(shí)驗(yàn)室,北京 100012)

    渾河流域2010—2014年的魚類群落和水生態(tài)健康變化分析

    劉偉1,2張遠(yuǎn)1,2高欣1,2賈曉波1,2馬淑芹1,2劉思思1,2

    (1. 中國環(huán)境科學(xué)研究院流域水生態(tài)保護(hù)技術(shù)研究室,北京 100012; 2. 環(huán)境標(biāo)準(zhǔn)與風(fēng)險(xiǎn)評估國家重點(diǎn)實(shí)驗(yàn)室,北京 100012)

    為了解渾河流域魚類群落結(jié)構(gòu)的變化趨勢和水生態(tài)健康變化,分別于2010年和2014年對渾河流域46個(gè)采樣點(diǎn)進(jìn)行了魚類和環(huán)境因子調(diào)查。其中2010年采集到魚類15039尾,隸屬于6目9科32種; 2014年采集到魚類10483尾,隸屬于6目10科41種。Mann-Whitney U檢驗(yàn)表明從2010年到2014年魚類總物種數(shù)、香農(nóng)多樣性指數(shù)、底層物種數(shù)百分比、植食性和肉食性個(gè)體數(shù)百分比等魚類特征參數(shù)顯著上升。典型對應(yīng)分析(CCA)結(jié)果表明,2010年顯著影響魚類群落結(jié)構(gòu)的環(huán)境因子為流速、電導(dǎo)率、河流等級和鈣離子濃度; 而2014年顯著影響魚類群落結(jié)構(gòu)的環(huán)境因子則為電導(dǎo)率、河流等級和氨氮濃度。2010年和2014年的魚類完整性指數(shù)(FIBI)評分分別為48.76±24.82和50.41±17.35,與2010年的評價(jià)結(jié)果相比,2014年F-IBI結(jié)果稍有改善,其中極好和好的點(diǎn)位數(shù)分別少3個(gè)和1個(gè),一般的點(diǎn)位數(shù)多9個(gè),而極差和差的點(diǎn)位分別少4個(gè)和1個(gè)。

    魚類;群落結(jié)構(gòu);生物完整性指數(shù);渾河流域

    魚類作為河流生態(tài)系統(tǒng)食物鏈的頂級群落,對各種形式的人類干擾較為敏感,鑒定分類信息完善,壽命較長能提供時(shí)間連續(xù)性的評價(jià),因此在河流生態(tài)健康評價(jià)中占有重要的地位[1,2]。研究發(fā)現(xiàn),魚類的群落結(jié)構(gòu)受水文條件[3]、地形地貌[4]、水環(huán)境質(zhì)量[5]、生物作用[6]等多種因素的共同作用,因此魚類群落結(jié)構(gòu)與一系列的環(huán)境變量存在可預(yù)測的關(guān)系。但由于不同時(shí)期的干擾強(qiáng)度不同且干擾具有時(shí)間累加效應(yīng),魚類群落結(jié)構(gòu)會隨著時(shí)間而發(fā)生變化,影響魚類群落結(jié)構(gòu)的環(huán)境因子也隨之變化[7]。因此,分析魚類群落結(jié)構(gòu)的時(shí)間變化趨勢及不同時(shí)期影響魚類群落結(jié)構(gòu)的主要因子對河流生態(tài)系統(tǒng)的恢復(fù)和管理具有重要意義。

    隨著國家水專項(xiàng)的實(shí)施,遼河流域有機(jī)污染物得到了顯著削減,水質(zhì)逐漸改善,水生態(tài)健康水平不斷提高[8]。研究者針對遼河水生態(tài)系統(tǒng)展開了細(xì)致而廣泛的調(diào)查工作,取得了重要的成果[9,10],基于魚類和藻類的生物完整性評價(jià)也在這一區(qū)域廣泛開展[11,12]。目前,國內(nèi)針對河流水生態(tài)的調(diào)查工作,一般都是一次評價(jià)[13]或不同季節(jié)的比較[14],即在相對集中的1—2年內(nèi),對某一流域進(jìn)行水生生物調(diào)查,從而構(gòu)建生物完整性評價(jià)體系對該流域進(jìn)行健康評價(jià)。然而由于國家環(huán)境保護(hù)措施的開展,某一流域的魚類群落和健康狀況也會隨之發(fā)生變化,導(dǎo)致先前的健康評價(jià)具有落后性和不適性。通過不同年份的多次調(diào)查可以為管理者提供科學(xué)及時(shí)有效的流域健康狀況信息從而有效地進(jìn)行受損流域生態(tài)系統(tǒng)的恢復(fù)。因此,本文基于2010年和2014年在渾河流域46個(gè)樣點(diǎn)魚類調(diào)查的基礎(chǔ)上,分析4年中由于國家水專項(xiàng)的開展,渾河流域魚類群落結(jié)構(gòu)的變化及影響魚類群落因子的變化,并基于魚類完整性指數(shù)進(jìn)行河流健康評價(jià),分析渾河流域生態(tài)健康變化趨勢,以期為渾河流域水生態(tài)系統(tǒng)的保護(hù)和管理提供數(shù)據(jù)支撐。

    1 材料與方法

    1.1研究區(qū)域

    渾河流域是遼河重要支流之一,位于遼寧省東部地區(qū)(40.71°-42.17°N,122.13°-125.21°E),全長為415 km,流域面積達(dá)1.14×104km2,發(fā)源于撫順市清源縣長白山支脈滾馬嶺。渾河為不對稱水系,東側(cè)支流密集,坡陡谷深,水量豐富; 西側(cè)支流很少,水量不大。渾河流經(jīng)遼寧中部城市群,重工業(yè)發(fā)達(dá),人口稠密,為沿岸城市廢水排放的主要通道。渾河上游建有大伙房水庫,集水面積5437 km2,庫容量為21.87×108m3。大伙房水庫以上區(qū)域?yàn)榈蜕角鹆陞^(qū),森林覆蓋率高; 中下游為平原地區(qū),工業(yè)區(qū)密集,土地開發(fā)程度高。

    1.2數(shù)據(jù)獲取

    采樣點(diǎn)設(shè)置依據(jù)渾河水系的分布并綜合考慮野外的可達(dá)性及可操作性,在渾河干流及10條支流(蘇子河、紅河、英額河、社河、章黨河、古城河、李石河、白塔堡河、細(xì)河和蒲河)設(shè)置46個(gè)采樣點(diǎn)(圖1),并分別于2010年6月和2014年8月對這46個(gè)采樣點(diǎn)進(jìn)行魚類調(diào)查。

    樣品采集對于可涉水河段(<1.5 m),沿河岸在采樣點(diǎn)上下游各設(shè)置200 m長度的調(diào)查范圍,采用電魚法,每個(gè)樣點(diǎn)魚類采集時(shí)間為30min。對于不可涉水的河段(≥1.5 m),選用電魚法和掛網(wǎng)法(網(wǎng)徑6 cm×6 cm,12 cm×12 cm和20 cm×20 cm三種)配合進(jìn)行,電魚法只在靠近河岸的淺水區(qū)進(jìn)行; 掛網(wǎng)在深水區(qū)進(jìn)行,在不同類型的生境下布設(shè)網(wǎng)具,掛網(wǎng)時(shí)間為60min。待樣品采集后,在野外進(jìn)行測/稱量和鑒定,記錄魚的全長、體長、體重和數(shù)量等。對于現(xiàn)場不能鑒定的樣品,用5%的甲醛保存魚類樣品,帶回實(shí)驗(yàn)室鑒定。魚類鑒定主要依據(jù)《東北地區(qū)淡水魚類》[15]。

    環(huán)境因子調(diào)查環(huán)境因子的調(diào)查包括自然地理因子、水化學(xué)因子、物理因子和棲息地質(zhì)量(QHEI)調(diào)查。

    自然地理因子包括:海拔和河流等級。其中,海拔通過GPS(Trimble-Juno SB)測得; 河流等級用Strahler(1957)提出的方法從1∶25萬水系圖獲取。

    水化學(xué)因子包括:pH、溶解氧含量(DO)、電導(dǎo)率(EC)、總?cè)芙庑怨腆w(TDS)、懸浮物固體含量(SS)、總氮(TN)、總磷(TP)、高錳酸鹽指數(shù)(CODMn)、鈣離子(Ca2+)、鎂離子(Mg2+)、堿度(Alk)、氨氮(NH+4-N)、硝酸鹽氮(NO-3-N)、活性磷酸鹽(PO34-)。其中,pH、DO、EC和TDS利用手持式Y(jié)SI 多參數(shù)水質(zhì)分析儀(YSI ProPlus)現(xiàn)場測定,其他指標(biāo)通過采集水樣,低溫運(yùn)送至實(shí)驗(yàn)室,48h內(nèi)按照《水和廢水分析檢測方法》中規(guī)定方法進(jìn)行測定。

    圖1 渾河流域采樣點(diǎn)分布Fig. 1 Location of the sample sites in the Hun River Basin

    物理因子包括水深(Depth)、流速(Flow rate)和河寬(Width)。具體測量方法為:分別在采樣點(diǎn)和采樣點(diǎn)上下游50 m河段內(nèi)取一個(gè)斷面,采用激光測距儀(Leupold RX-IV)測量斷面濕寬,同時(shí)每個(gè)斷面從河岸開始每隔1 m用流速儀(FP 201)測量一次水深和流速直到對面河岸,流速采用五點(diǎn)法計(jì)算,最終結(jié)果取3個(gè)斷面的平均值。

    生境質(zhì)量調(diào)查評價(jià)參照鄭丙輝等[16]的研究,選取10個(gè)評價(jià)指標(biāo)(底質(zhì)、棲境復(fù)雜性、速度和深度結(jié)合、堤岸穩(wěn)定性、河道變化、河水水量狀況、植被多樣性、水質(zhì)狀況、人類活動強(qiáng)度、河岸邊土壤利用類型),每項(xiàng)20分,共計(jì)4個(gè)得分等級(好16—20、較好11—15、一般6—10、差0—5)。

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

    魚類群落結(jié)構(gòu)參照以往的研究[11,13]并結(jié)合魚類的習(xí)性、棲息地等信息,將池沼公魚(Hypomesus olidus McAllister)、馬口魚(Opsariichthys bidens Günther)、洛氏鱥(Phoxinus lagowskii Dybowsky)、北方花鰍(Cobitis granoei Rendahl)、點(diǎn)紋銀鮈(Squalidus wolterstorffi Regan)、沙塘鱧(Odontobutis obscurus Temminck et Schlegel)視為敏感種; 將鯽(Carassius auratus Linnaeus)、麥穗魚(Pseudorasbora parva Temminck et Schlegel)、(Hemiculter leucisculus Basilewsky)、泥鰍(Misgurnus anguillicaudatus Cantor)視為耐受種。

    采用Mann-Whitney U檢驗(yàn)分析不同年份魚類特征參數(shù)和環(huán)境因子的差異顯著性。對2010年和2014年各樣點(diǎn)魚類多度數(shù)據(jù)進(jìn)行去趨勢對應(yīng)分析(Detrended correspondence analysis,DCA)發(fā)現(xiàn),最長軸的梯度長度分別為3.1和3.6,介于3.0—4.0,均可以采用典型對應(yīng)分析(CCA)[17]。在對魚類數(shù)據(jù)和環(huán)境因子數(shù)據(jù)進(jìn)行典型對應(yīng)分析時(shí),先對所有數(shù)據(jù)進(jìn)行l(wèi)og10(x+1)轉(zhuǎn)換,然后將通過999次蒙特卡羅檢驗(yàn)后單個(gè)因子貢獻(xiàn)率的P值小于0.05的環(huán)境因子作為顯著影響魚類群落的環(huán)境因子。此外,在進(jìn)行DCA和CCA分析前,需剔除出現(xiàn)頻率低于5%的魚類,以避免稀有種的干擾。Mann-Whitney U檢驗(yàn)于SPSS 19.0中進(jìn)行,DCA和CCA分析在Canoco 5.0中進(jìn)行。

    魚類生物完整性(F-IBI)評價(jià) (1)核心參數(shù)的篩選:根據(jù)渠曉東等建立的標(biāo)準(zhǔn)化參照點(diǎn)和受損點(diǎn)的選取方法[18],綜合考慮水質(zhì)和棲息地質(zhì)量評分篩選敏感點(diǎn)(水質(zhì)達(dá)到Ⅱ類水體,棲息地質(zhì)量評分>120,人類活動強(qiáng)度和河岸邊土壤利用類型兩項(xiàng)得分之和≥24)和受損點(diǎn)(水質(zhì)等級低于Ⅳ級,棲息地質(zhì)量評分<90,人類活動強(qiáng)度和河岸邊土壤利用類型兩項(xiàng)得分之和≤18)。參考裴雪嬌等[11]、黃亮亮等[19]和張浩等[20]確立的候選指標(biāo),結(jié)合本研究實(shí)際采集魚類的情況和數(shù)據(jù)的可獲得性,確立了5類23個(gè)候選指標(biāo)(表1)。參照王備新等[21]的研究方法,對這23個(gè)指標(biāo)分別進(jìn)行分布范圍、判別能力和相關(guān)性分析(|R|<0.9)的檢驗(yàn),最終篩選出F-IBI的核心參數(shù)。(2)比值法計(jì)算F-IBI值:王備新等[21]研究發(fā)現(xiàn)比值法的準(zhǔn)確率明顯優(yōu)于3分法和4分法,因此本研究采用比值法計(jì)算F-IBI核心參數(shù)的分值,比值法具體計(jì)算公式見裴雪嬌等[11]。(3)IBI指標(biāo)體系評價(jià)標(biāo)準(zhǔn):參照渠曉東等[18]研究中的方法,將所有樣點(diǎn)得分的95%分位數(shù)值作為健康評價(jià)標(biāo)準(zhǔn),若樣點(diǎn)的F-IBI值大于95%分位數(shù)值,則該樣點(diǎn)為極好。將小于95%分位數(shù)的分布范圍進(jìn)行4等分,確定出好、一般、差、極差4個(gè)等級的劃分標(biāo)準(zhǔn)。

    2 結(jié)果

    2.1渾河流域魚類群落結(jié)構(gòu)和環(huán)境因子的變化趨勢

    于2010年在渾河流域46個(gè)采樣點(diǎn)共采集到魚類15039尾,隸屬于6目9科32種,其中種類最多的為鯉科(Cyprinidae,20種),占62.5%; 其次為鰍科(Cobitidae,4種),占12.5%; 塘鱧科(Eleotridae,2種),占6.25%,漁獲量較大的物種分別為洛氏鱥(35.52%)、棒花魚(Abbottina rivularis Basiewsky,20.42%)、寬鰭鱲(Zacco platypus Temminck et Schlegel,13.03%)、鯽(6.80%)和麥穗魚(4.04%),占總漁獲量的79.81%。而在2014年在相同的采樣點(diǎn)采集到魚類10483尾,隸屬于6目10科41種,種類數(shù)最多的為鯉科(21種),占51.22%; 其次為鰍科(6種),占14.63%;鰕虎魚科(Gobiidae,4種),占9.76%,漁獲量較大的物種分別為洛氏鱥(22.35%)、北方條鰍(Nemachilus nudes Bleeker,18.10%)、棒花魚(9.49%)、寬鰭鱲(9.44%)和波氏吻蝦虎(Rhinogobius cliffordpopei Nichols,9.33%),占總漁獲量的68.71%。

    通過對魚類群落結(jié)構(gòu)的Mann-Whitney U檢驗(yàn)分析(表2),可以發(fā)現(xiàn),從2010年到2014年,魚類種類數(shù)和香農(nóng)多樣性指數(shù)顯著上升,而魚類個(gè)體總數(shù)則無明顯差異。鯉科物種數(shù)百分比顯著下降,鰕虎魚科和塘鱧科物種數(shù)百分比則顯著升高。此外,中下層魚類物種數(shù)百分比顯著下降,而底層魚類物種數(shù)百分比顯著上升。植食性和肉食性個(gè)體數(shù)百分比都有顯著上升,其余參數(shù)的變化則不顯著。

    通過Mann-Whitney U檢驗(yàn)分析2010年和2014年渾河流域的環(huán)境因子(表3),發(fā)現(xiàn)pH、電導(dǎo)率、鈣離子濃度和堿度這些表征無機(jī)性污染的指標(biāo)顯著性上升,而水深、流速、溶解氧和表征有機(jī)污染的總氮、高錳酸鹽指數(shù)、氨氮及磷酸根濃度這些指標(biāo)顯著下降,即相對于2010年,2014年的水深變淺,流速變小,有機(jī)性污染減輕而無機(jī)性污染加重。

    表1 F-IBI候選指標(biāo)列表Tab. 1 List of F-IBI candidate metrics for the Hun River Basin

    表2 渾河流域兩個(gè)年份魚類特征參數(shù)比較Tab. 2 The comparison of fish parameters between two years in the Hun River Basin

    2.2影響渾河流域魚類群落結(jié)構(gòu)的環(huán)境因子隨時(shí)間的變化

    分別對2010年和2014年渾河流域的環(huán)境因子進(jìn)行CCA分析(圖2),結(jié)果表明:在2010年顯著影響渾河流域魚類群落結(jié)構(gòu)的環(huán)境因子為河流等級、鈣離子濃度、電導(dǎo)率和流速; 而在2014年顯著影響渾河流域魚類群落結(jié)構(gòu)的因子為河流等級、電導(dǎo)率和氨氮濃度。其中,兩個(gè)CCA圖都表現(xiàn)出電導(dǎo)率和河流等級與第一軸正相關(guān),且位于干流下游的電導(dǎo)率明顯高于上游干流和周邊支流,位于大伙房水庫上游的支流點(diǎn)位都處在第二、三象限,與河流等級呈負(fù)相關(guān)。在2010年的樣點(diǎn)與環(huán)境因子的CCA圖中,鈣離子濃度與河流等級相關(guān)性較高,表明隨著河流等級的升高,鈣離子濃度也不斷升高。此外,位于中上游點(diǎn)位的流速明顯高于下游點(diǎn)位的流速。而在2014年樣點(diǎn)與環(huán)境因子的CCA圖表明,氨氮濃度與第一軸正相關(guān),且下游區(qū)域的濃度明顯高于中上游區(qū)域。

    表3 渾河流域兩個(gè)年份環(huán)境因子比較Tab. 3 The comparison of environmental factors between two years in the Hun River Basin

    圖2 渾河流域2010年(a)和2014年(b)環(huán)境因子與魚類群落的典型對應(yīng)分析圖Fig. 2 CCA biplot of environmental variables and species composition of 2010(a) and 2014(b) in the Hun River Basin

    2.3魚類完整性評價(jià)結(jié)果隨時(shí)間的變化

    構(gòu)建渾河流域魚類完整性的核心參數(shù)分布范圍和評分方法如表4所示,基于F-IBI的渾河健康評價(jià)標(biāo)準(zhǔn)如表5所示。渾河流域2010年和2014年的魚類完整性(F-IBI)評價(jià)結(jié)果如圖3所示。

    渾河流域2010年和2014年兩個(gè)年份的魚類完整性評價(jià)結(jié)果大致相同,即渾河上游的支流蘇子河、英額河、紅河和社河健康狀況較好,渾河中游干流和周邊支流的健康狀況一般,而位于渾河下游的渾河干流及支流健康狀況較差。經(jīng)Pearson相關(guān)性分析發(fā)現(xiàn),2010年和2014年的F-IBI的相關(guān)系數(shù)為0.756(P<0.01)。2010年和2014年的F-IBI評分分別為48.76±24.82和50.41±17.35,2010年評價(jià)結(jié)果為極好的點(diǎn)為4個(gè),占總樣點(diǎn)的8.7%; 好的點(diǎn)為14個(gè),占30.4%; 一般的點(diǎn)為11個(gè),占23.9%; 差的點(diǎn)為11個(gè),占23.9%; 極差的點(diǎn)為6個(gè),占13.1%。而2014年評價(jià)結(jié)果為極好、好、一般、差和極差的點(diǎn)位數(shù)分別為1(2.2%)、13(28.3%)、20(43.5%)、10(21.7%)和2(4.3%)。與2010年的F-IBI評價(jià)結(jié)果相比,2014年魚類完整性評價(jià)結(jié)果中極好和好的點(diǎn)位數(shù)分別少3個(gè)和1個(gè),一般的點(diǎn)位數(shù)多9個(gè),而極差和差的點(diǎn)位分別少4個(gè)和1個(gè)。

    表4 渾河流域魚類生物完整性核心參數(shù)的分布范圍和評分方法Tab. 4 Distribution and scoring criteria of core F-IBI metrics in the Hun River Basin

    表5 基于F-IBI的渾河流域健康評價(jià)標(biāo)準(zhǔn)Tab. 5 Assessment criteria of health for F-IBI of two years in the Hun River Basin

    3 討論

    3.1不同年份魚類群落結(jié)構(gòu)和環(huán)境因子的差異

    渾河作為遼河水系的水資源最豐富的河流,是遼寧省中心城市圈的重要水源地,但是由于近年來受到市區(qū)排放的污廢水及沿岸農(nóng)田面源污染的影響[22],渾河水質(zhì)嚴(yán)重惡化,導(dǎo)致魚類群落結(jié)構(gòu)發(fā)生急劇變化。根據(jù)遼河地區(qū)魚類歷史資料[15],20世紀(jì)80年代渾太河水系共鑒定出魚類62種,其中鯉科魚類37種,鰍科魚類3種(渾河是渾太河水系的主要分支,其魚類種類與渾太河水系相差不大)。而20世紀(jì)90年代后期渾河共鑒定出魚類55種,其中鯉科魚類36種,鰍科4種[23]。對比此次調(diào)查,2010年和2014年分別鑒定出魚類32種和41種??梢姡?0年來由于水質(zhì)惡化,渾河流域的魚類種數(shù)明顯減少。

    魚類群落結(jié)構(gòu)會隨著年份而發(fā)生變化[24],對比2010年和2014年的魚類群落結(jié)構(gòu),發(fā)現(xiàn)2014年鑒定出的魚類比2010年多9種,樣點(diǎn)的種類數(shù)和香農(nóng)多樣性指數(shù)明顯增大,且植食性和肉食性個(gè)體百分比都較2010年顯著上升,底層魚類物種數(shù)百分比由于底棲性鰕虎魚科的出現(xiàn)而顯著上升,表明魚類群落功能由單一向復(fù)雜化發(fā)展,水生態(tài)質(zhì)量在緩慢恢復(fù)中。這與長江靖江段魚類物種多樣性的年際變化類似,由于禁漁期制度的全面實(shí)施,從2002—2011年,魚類物種豐富度和多樣性呈顯著上升的趨勢[25]。對比渾河兩個(gè)年份的環(huán)境因子,2014年有機(jī)污染狀況較2010年顯著改善,而無機(jī)污染狀況則顯著加劇。這主要是由于水專項(xiàng)的開展,遼河流域作為重點(diǎn)流域進(jìn)行了水污染治理,控制重點(diǎn)主要集中在COD和NH+4-N的控制減排上[8],因此有機(jī)污染顯著減輕,而無機(jī)污染情況由于挖沙、開山、采礦等活動的廣泛開展而加重。由于有機(jī)污染物排放的消減,渾河水環(huán)境明顯改善,魚類群落功能復(fù)雜化,種類明顯增多,底層魚類種類數(shù)也因?yàn)橛袡C(jī)物沉積減少而顯著上升。

    圖3 渾河流域2010年(a)和2014年(b)的魚類生物完整性評價(jià)Fig. 3 The assessment of environmental condition of 2010(a) and 2014(b) in the Hun River Basin

    3.2不同年份影響魚類群落結(jié)構(gòu)的環(huán)境因子差異性

    魚類群落結(jié)構(gòu)受到環(huán)境因子影響的機(jī)理是非常復(fù)雜的,不同河流、不同季節(jié)影響魚類群落結(jié)構(gòu)的環(huán)境因子都不相同[10,26,27]。本調(diào)查發(fā)現(xiàn)不同年份影響魚類群落結(jié)構(gòu)的環(huán)境因子也不相同,對兩個(gè)年份魚類群落結(jié)構(gòu)都有顯著影響的環(huán)境因子為河流等級和電導(dǎo)率,此外在2010年影響魚類群落結(jié)構(gòu)的影響因子還有鈣離子濃度和流速,而在2014年影響魚類群落的顯著性因子還包括氨氮濃度。從環(huán)境與樣點(diǎn)的CCA圖可以看出,在2010年和2014年電導(dǎo)率較高的點(diǎn)均集中在渾河下游干支流,這里的人口集中,人為干擾大,污染集中,魚類物種數(shù)相對較少。河流等級對魚類群落結(jié)構(gòu)具有較大影響,隨著河流等級的增大,棲息地的復(fù)雜性和多樣性也隨之增加,魚類物種數(shù)也出現(xiàn)線性或非線性的增加[28]。盡管河流等級在2010年和2014年均為魚類群落結(jié)構(gòu)的顯著環(huán)境因子,但其作用機(jī)制不盡相同。經(jīng)Kruskal-Wallis檢驗(yàn)發(fā)現(xiàn),在2010年,渾河流域四級河流的魚類物種數(shù)和個(gè)體數(shù)最大,其次為二級河流,五級河流最低; 而在2014年,渾河流域四級河流魚類物種數(shù)和香農(nóng)多樣性指數(shù)最大,其次為二級河流,一級河流最低,這可能是由于受到了中度干擾理論的影響[28]。在2014年,一級河流物種數(shù)最少,主要原因是相比較于2010年,渾河流域中下游污染控制措施較嚴(yán),水質(zhì)得到了改善,而上游源頭區(qū)由于旅游開發(fā)和礦產(chǎn)資源的需求而不斷受到人為干擾。流速僅在2010年是顯著影響魚類群落結(jié)構(gòu)的環(huán)境因子,這可能主要是因?yàn)樵?010年水流速度較高,對底質(zhì)沖刷較嚴(yán)重,導(dǎo)致底層魚類物種數(shù)較少,從而魚類群落結(jié)構(gòu)受到影響。此次調(diào)查中也發(fā)現(xiàn)鈣離子僅是2010年影響群落結(jié)構(gòu)的顯著性因子,從2010年到2014年,鈣離子濃度顯著升高,經(jīng)過逐漸的適應(yīng)過程,可能魚體對鈣離子的耐受程度加強(qiáng),從而導(dǎo)致2014年鈣離子已經(jīng)不再是顯著影響魚類群落結(jié)構(gòu)的環(huán)境因子。氨氮是水體有機(jī)物污染的重要指標(biāo),其對魚類有毒害作用。由于渾河流域水體污染具有氨氮濃度較高的特點(diǎn)[22],而Amisah和Cowx[29]的研究表明當(dāng)氨氮嚴(yán)重超標(biāo)后會導(dǎo)致魚類物種數(shù)和個(gè)體數(shù)大幅減少。2010年的渾河流域的氨氮濃度均值大于1.5 mg/L,在此濃度下,魚類物種數(shù)大幅減少,相對耐污的鯉科魚類比例較高。而在2014年氨氮濃度由于減排措施的展開而大幅消減,鯉科魚類比例顯著下降,出現(xiàn)對氨氮敏感的物種。此外,氨氮作為藻類的重要營養(yǎng)物質(zhì),會影響藻類的生長和繁殖[12],從而間接影響植食性魚類的數(shù)量。在2014年,植食性魚類顯著上升,對氨氮耐受性較高的鯉科魚類比例下降,因此盡管2014年氨氮濃度較低,但其是影響魚類群落結(jié)構(gòu)的顯著性因子。因而當(dāng)流域進(jìn)行污染物的減排時(shí),由于水體的自凈作用和污染物來源的減少,引起水體環(huán)境因子的變化,魚類群落結(jié)構(gòu)也隨之變化,同時(shí)由于魚類適應(yīng)性地形成、新物種的出現(xiàn)、耐受性鯉科物種的減少及魚類群落結(jié)構(gòu)的復(fù)雜化,導(dǎo)致影響魚類群落結(jié)構(gòu)的顯著性因子也發(fā)生變化,因此分析不同年份影響魚類群落結(jié)構(gòu)的顯著性因子能及時(shí)有效地為管理部門提供數(shù)據(jù)支撐,形成一個(gè)良性的反饋機(jī)制,對于魚類群落結(jié)構(gòu)對環(huán)境因子的響應(yīng)機(jī)制和對受損的生態(tài)系統(tǒng)進(jìn)行恢復(fù)具有重要借鑒作用。

    3.3不同年份水生態(tài)健康狀況的差異

    年跟蹤評價(jià)是完整性評價(jià)中的重要方法,可用來分析人為干擾對生物完整性的影響[30]。Zhu和Chang[31]比較了長江上游1997—2002年的魚類完整性變化,發(fā)現(xiàn)隨著人為活動干擾的加強(qiáng),魚類完整性出現(xiàn)了明顯的退化。此次在渾河流域發(fā)現(xiàn)2014年的IBI的平均評分為50.42,略高于2010年的IBI平均評分48.76,說明兩年的評價(jià)結(jié)果總體上相差不太大,但2014年的健康狀況稍有提高。從每一個(gè)點(diǎn)的健康狀況變化來看,2010年到2014年,健康評價(jià)上升一個(gè)等級的點(diǎn)位12個(gè),上升兩個(gè)等級的點(diǎn)為1個(gè),下降一個(gè)等級的點(diǎn)位有8個(gè),下降兩個(gè)等級的點(diǎn)為2個(gè),評價(jià)結(jié)果沒變化的點(diǎn)位為23個(gè)。這表明2014年的河流健康狀況較2010年有了小幅的提升,且河流評價(jià)等級上升的點(diǎn)位主要分布在渾河流域的中下游區(qū)域和上游支流的末端區(qū)域,下降的點(diǎn)位主要分布在中上游支流的源頭區(qū)域,主要是水專項(xiàng)開展以來加強(qiáng)了渾河流域污染物的減排工作,因此2014年渾河流域的有機(jī)污染物相對于2010年顯著減少。而渾河流域的中下游區(qū)域企業(yè)集中,人口密集,開展減排后效果最明顯,隨著污染物濃度的減少,河流健康等級也隨之上升。而渾河流域上游建有大伙房水庫,由于水庫保護(hù)措施的開展,上游支流入大伙房水庫點(diǎn)位的健康等級也因此提升。而源頭區(qū)域由于一般路況復(fù)雜,監(jiān)管比較困難,隨著人口的增長,源頭區(qū)域開始受到人為干擾,此外由于礦石開挖和挖沙等活動的開展,這些區(qū)域的生境受到破壞,同時(shí)水體電導(dǎo)率、硬度和堿度顯著升高,導(dǎo)致這些區(qū)域健康等級下降。

    [1]Harris J H. The use of fish in ecological assessments [J]. Australian Journal of Ecology,1995,20(1):65—80

    [2]Trautwein C,Schinegger R,Schmutz S. Cumulative effects of land use on fish metrics in different types of running waters in Austria [J]. Aquatic Sciences,2012,74(2):329—341

    [3]Marchetti M P,Moyle P B. Effects of flow regime on fish assemblages in a regulated California stream [J]. Ecological Applications,2001,11(2):530—539

    [4]Bhatt J P,Manish K,Pandit M K. Elevational gradients in fish diversity in the Himalaya:Water discharge is the key driver of distribution patterns [J]. PloS One,2012,7(9):e46237

    [5]Figuerola B,Maceda-Veiga A,De Sostoa A. Assessing the effects of sewage effluents in a Mediterranean creek:fish population features and biotic indices [J]. Hydrobiologia,2012,694(1):75—86

    [6]Súarez Y R,Petrere Jr M,Catella A C. Factors determining the structure of fish communities in Pantanal lagoons(MS,Brazil) [J]. Fisheries Management and Ecology,2001,8(2):173—186

    [7]Yu H C. Fish community structure in the Changjiang Estuary and adjacent waters [D]. Thesis for Master of Science. Institute of Oceanology,Chinese Academy of Sciences,Qingdao. 2008 [于海成. 長江口及鄰近海域魚類群落結(jié)構(gòu)分析. 青島:中國科學(xué)院海洋研究所. 2008]

    [8]Meng W. Construction of technology system for water-shed water pollution treatment and management- exploration and practice of the major water program in Liaohe River Basin [J]. Engineering Science,2013,15(3):4—10[孟偉. 遼河流域水污染治理和水環(huán)境管理技術(shù)體系構(gòu)建-國家重大水專項(xiàng)在遼河流域的探索與實(shí)踐. 中國工程科學(xué),2013,15(3):4—10]

    [9]Wang W,Wang B,He X Y,et al. Study of zoning and distribution characteristics of fish in Taizi River [J]. Research of Environmental Sciences,2013,26(5):494—501[王偉,王冰,何旭穎,等. 太子河魚類群落結(jié)構(gòu)的空間分布特征. 環(huán)境科學(xué)研究,2013,26(5):494—501]

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

    [11]Pei X J,Niu C J,Gao X,et al. The ecological health assessment of Liao River Basin,China,based on biotic integrity index of fish [J]. Acta Ecologica Sinica,2010,30(21):5736—5746 [裴雪嬌,牛翠娟,高欣,等. 應(yīng)用魚類完整性評價(jià)體系評價(jià)遼河流域健康. 生態(tài)學(xué)報(bào),2010,30(21):5736—5746]

    [12]Yin X W,Zhang Y,Qu X D,et al. Community structure and biological integrity of periphyton in Hunhe River water system of Liaoning Province,Northeast China [J]. Chinese Journal of Applied Ecology,2011,22(10):2732—2740 [殷旭旺,張遠(yuǎn),渠曉東,等. 渾河水系著生藻類的群落結(jié)構(gòu)與生物完整性. 應(yīng)用生態(tài)學(xué)報(bào),2011,22(10):2732—2740]

    [13]Wu W,Xu Z X,Yin X W,et al. Fish community structure and the effect of environmental factors in the Wei River basin [J]. Acta Scientiae Circumstantiae,2014,34(5):1298—1308 [武瑋,徐宗學(xué),殷旭旺,等. 渭河流域魚類群落結(jié)構(gòu)特征及其受環(huán)境因子的影響分析. 環(huán)境科學(xué)學(xué)報(bào),2014,34(5):1298—1308]

    [14]Song M J,Deng H T,Zhu F Y,et al. Seasonal dynamics of macroinvertebrates community structure in Daning river after a 175 m depth of impoundment in the Three Gorges Reservoir [J]. Acta Hydrobiologica Sinica,2015,39(5):1046—1053 [宋明江,鄧華堂,朱峰躍,等. 三峽水庫175 m 蓄水后大寧河底棲動物群落結(jié)構(gòu)季節(jié)動態(tài). 水生生物學(xué)報(bào),2015,39(5):1046—1053]

    [15]Xie Y H. Freshwater Fishes in Northeast Region of China[M]. Shenyang:Liaoning Science and Technology Press. 2007,1—460 [解玉浩. 東北地區(qū)淡水魚類. 沈陽:遼寧科學(xué)技術(shù)出版社. 2007,1—460]

    [16]Zheng B H,Zhang Y,Li Y B. Study of indicators and methods for river habitat assessment of Liao River Basin[J]. Acta Scientiae Circumstantiae,2007,27(6):928—936 [鄭丙輝,張遠(yuǎn),李英博. 遼河流域河流棲息地評價(jià)指標(biāo)與評價(jià)方法研究. 環(huán)境科學(xué)學(xué)報(bào),2007,27(6):928—936]

    [17]Lep? L,?milauer P. Multivariate Analysis of Ecological Data Using CANOCO [M]. New York:Cambridge University Press. 2003,184—198

    [18]Qu X D,Liu Z G,Zhang Y. Discussion on the standardized method of reference sites selection for establishing the Benthic-Index of Biotic Integrity [J]. Acta Ecologica Sinica,2012,32(15):4661—4672 [渠曉東,劉志剛,張遠(yuǎn). 標(biāo)準(zhǔn)化方法篩選參照點(diǎn)構(gòu)建大型底棲動物生物完整性指數(shù). 生態(tài)學(xué)報(bào),2012,32(15):4661—4672]

    [19]Huang L L,Wu Z Q,Jiang K,et al. Development and application of IBI based on fish to assess the river's health in the East Tiaoxi River [J]. China Environmental Science,2013,33(7):1280—1289 [黃亮亮,吳志強(qiáng),蔣科,等. 東苕溪魚類生物完整性評價(jià)河流健康體系的構(gòu)建與應(yīng)用. 中國環(huán)境科學(xué),2013,33(7):1280—1289]

    [20]Zhang H,Ding S,Zhang Y,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 [張浩,丁森,張遠(yuǎn),等. 西遼河流域魚類生物完整性指數(shù)評價(jià)及其與環(huán)境因子的關(guān)系. 湖泊科學(xué),2015,27(5):829—839]

    [21]Wang B X,Yang L F,Hu B J,et al. A preliminary study on the assessment of stream ecosystem health in south of Anhui Province using Benthic-Index of Biotic Integrity[J]. Acta Ecologica Sinica,2005,25(6):1481—1490 [王備新,楊蓮芳,胡本進(jìn),等. 應(yīng)用底棲動物完整性指數(shù)BIBI評價(jià)溪流健康. 生態(tài)學(xué)報(bào),2005,25(6):1481—1490]

    [22]Zhang H L,Sun L N,Luo Q,et al. Seasonal variation and sources of pollutants in surface water of Hunhe River [J]. Chinese Journal of Ecology,2011,30(1):119—125 [張宏齡,孫麗娜,羅慶,等. 渾河流域水體污染的季節(jié)性變化及來源. 生態(tài)學(xué)雜志,2011,30(1):119—125]

    [23]Shi Y Q,Tong J,Han T X. Hydrobios diversity in Mid and Upper Stream of Hun River and its protection [J]. Environmental Protection Science,2002,28(110):7—9 [史玉強(qiáng),童君,韓天雪. 渾河中、上游水生生物多樣性及其保護(hù). 環(huán)境保護(hù)科學(xué),2002,28(110):7—9]

    [24]Yu H C,Xian W W. The environment effect on fish assemblage structure in waters adjacent to the Changjiang(Yangtze) River estuary(1998—2001) [J]. Chinese Journal of Oceanology and Limnology,2009,27(3):443—456

    [25]Sun S S. The research of the temporal pattern of fish assemblage in coastal Jingjiang section of the Yangtze River over last decade [D]. Thesis for Master of Science. Shanghai Ocean University,Shanghai. 2013 [孫莎莎. 最近10年長江靖江沿岸魚類群聚特征的時(shí)間格局研究.碩士學(xué)位論文,上海海洋大學(xué),上海. 2013]

    [26]Liu B,Zhang Y,Qu X D,et al. Community structure and diversity of fishes in the Liao River Nature Reserve [J]. Freshwater Fisheries,2013,43(3):49—53 [劉斌,張遠(yuǎn),渠曉東,等. 遼河干流自然保護(hù)區(qū)魚類群落結(jié)構(gòu)及其多樣性變化. 淡水漁業(yè),2013,43(3):49—55]

    [27]Gao X,Zhang Y,Ding S,et al. Response of fish communities to environmental changes in an agriculturally dominated watershed(Liao River Basin) in northeastern China [J]. Ecological Engineering,2015,76:130—141

    [28]Yan Y Z,Zhan Y J,Chu L,et al. Effects of stream size and spatial position on stream-dwelling fish assemblages [J]. Acta Hydrobiologica Sinica,2010,34(5):1022—1030[嚴(yán)云志,占姚軍,儲玲,等. 溪流大小及其空間位置對魚類群落結(jié)構(gòu)的影響. 水生生物學(xué)報(bào),2010,34(5):1022—1030]

    [29]Amisah S,Cowx I G. Response of the fish populations of the River Don in South Yorkshire to water quality and habitats improvements [J]. Environmental Pollution,2000,108(2):191—199

    [30]Teels B M,Mazanti L E,Rewa C A. Using an IBI to assess effectiveness of mitigation measures to replace loss of a wetland-stream ecosystems. Wetlands,2004,24(2):375—384

    [31]Zhu D,Chang J. Annual variations of biotic integrity in the upper Yangtze River using an adapted index of biotic integrity(IBI) [J]. Ecological Indicators,2008,8(5):564—572

    VARIATION OF FISH COMMUNITY STRUCTURE AND WATER SYSTEM CONDITIONS FROM 2010 TO 2014 IN HUN RIVER OF LIAONING PROVINCE,NORTHEAST PART OF CHINA

    LIU Wei1,2,ZHANG Yuan1,2,GAO Xin1,2,JIA Xiao-Bo1,2,MA Shu-Qin1,2and LIU Si-Si1,2
    (1. Laboratory of Riverine Ecological Conservation and Technology,Chinese Research Academy of Environmental Sciences,Beijing 100012,China; 2. State Key Laboratory of Environmental Criteria and Risk Assessment,Beijing 100012,China)

    To understand the structure variation of fish assemblage of Hun River Basin in Liaoning Province of Northeast China,46 sampling sites were investigated in 2010 and 2014. A total of 32 species and 15039 fish,belonging to 9 families and 6 orders,and 41 species and 10483 fish,belonging to 10 families and 6 orders,were identified in 2010 and 2014,respectively. The Mann-Whitney U test supported a significant increase in fish assemblage parameters like the number of fish species,Shannon-Wiener index,proportion of benthic species,proportion of herbivores individuals and proportion of carnivores individuals. Canonical correspondence analysis showed that in 2010,flow rate,conductivity,river rank and calcium ion were the primary environmental factors affecting fish assemblage,while in 2014,conductivity,river rank and ammonia significantly affect fish assemblage. The score of F-IBI were 48.76±24.82 and 50.41±17.35 in 2010 and 2014 respectively. By contrast in 2014,the health condition slightly increased and the sampling sites of ‘very poor',‘poor',‘good' and ‘excellent' decreased by 4,1,1 and 3 sites than those in 2010,respectively. Meanwhile,the sampling sites of ‘moderate' increased by 9 sites than those in 2010.

    Fish assemblage; Community structure; Biotic integrity index; Hun River Basin

    Q145+.2

    A

    1000-3207(2016)05-0968-10

    10.7541/2016.125

    2015-09-06;

    2016-01-05

    水體污染控制與治理科技重大專項(xiàng)(2012ZX07501-001-02)資助 [Supported by Major Science and Technology Program for Water Pollution Control and Treatment in China(Grant No. 2012ZX07501-001-02)]

    劉偉(1991—),男,湖北黃岡人; 碩士; 主要從事流域生態(tài)學(xué)研究。E-mail:weizaitjdx@126.com

    張遠(yuǎn),E-mail:zhangyuan@craes.org.cn

    猜你喜歡
    渾河完整性魚類
    渾河三好橋
    基于MFCC和ResNet的魚類行為識別
    稠油熱采水泥環(huán)完整性研究
    云南化工(2021年9期)2021-12-21 07:44:00
    流域尺度的渾河水污染防治環(huán)境經(jīng)濟(jì)政策
    基于隨機(jī)森林的渾河沈陽段水質(zhì)評價(jià)
    奇妙的古代動物 泥盆紀(jì)的魚類
    莫斷音動聽 且惜意傳情——論音樂作品“完整性欣賞”的意義
    魚類是怎樣保護(hù)自己的
    精子DNA完整性損傷的發(fā)生機(jī)制及診斷治療
    魚類怎樣呼吸
    美女高潮到喷水免费观看| 1024视频免费在线观看| 亚洲av成人不卡在线观看播放网| 国产精品久久电影中文字幕| 十八禁人妻一区二区| videosex国产| 亚洲成av片中文字幕在线观看| 久久午夜综合久久蜜桃| 午夜精品在线福利| 黑丝袜美女国产一区| 日本 av在线| 欧美一区二区精品小视频在线| 欧美国产精品va在线观看不卡| 黄色a级毛片大全视频| 色av中文字幕| 久久午夜亚洲精品久久| 男女视频在线观看网站免费 | 天天一区二区日本电影三级| 国产精品永久免费网站| 国产视频内射| 午夜精品久久久久久毛片777| 男女视频在线观看网站免费 | 色综合站精品国产| 国产不卡一卡二| 午夜福利18| 男人舔女人的私密视频| 侵犯人妻中文字幕一二三四区| 精品卡一卡二卡四卡免费| 怎么达到女性高潮| 欧美一级a爱片免费观看看 | 视频区欧美日本亚洲| 国产免费av片在线观看野外av| 亚洲国产精品999在线| 久久婷婷成人综合色麻豆| 777久久人妻少妇嫩草av网站| 欧美乱色亚洲激情| 9191精品国产免费久久| 亚洲欧洲精品一区二区精品久久久| 国产精品美女特级片免费视频播放器 | 亚洲中文字幕日韩| 成人永久免费在线观看视频| 国产99久久九九免费精品| 怎么达到女性高潮| 亚洲天堂国产精品一区在线| 国产成人欧美在线观看| 在线观看免费午夜福利视频| tocl精华| 午夜福利在线观看吧| 国产精品久久久久久人妻精品电影| 亚洲中文av在线| 88av欧美| 一区二区日韩欧美中文字幕| 亚洲自拍偷在线| 悠悠久久av| 丝袜在线中文字幕| 亚洲成人免费电影在线观看| 久久国产精品影院| 天天一区二区日本电影三级| 韩国av一区二区三区四区| 欧美久久黑人一区二区| 黄色丝袜av网址大全| 很黄的视频免费| 国产欧美日韩一区二区精品| 国产成人影院久久av| 国产黄色小视频在线观看| 露出奶头的视频| 国产三级黄色录像| АⅤ资源中文在线天堂| 日韩欧美三级三区| 日本 av在线| 天天躁狠狠躁夜夜躁狠狠躁| www.自偷自拍.com| 国产精品影院久久| 精品国产超薄肉色丝袜足j| 午夜成年电影在线免费观看| 美女扒开内裤让男人捅视频| 亚洲成国产人片在线观看| 久99久视频精品免费| 日韩高清综合在线| 嫁个100分男人电影在线观看| 日韩欧美免费精品| 国产精品日韩av在线免费观看| 日韩欧美三级三区| 十八禁网站免费在线| 久久久久久久久久黄片| 岛国视频午夜一区免费看| 88av欧美| 国产av不卡久久| 88av欧美| 日韩欧美三级三区| 欧美黑人精品巨大| 精品国产一区二区三区四区第35| www.自偷自拍.com| 国内少妇人妻偷人精品xxx网站 | 19禁男女啪啪无遮挡网站| 国产蜜桃级精品一区二区三区| 男女午夜视频在线观看| 最好的美女福利视频网| 欧美三级亚洲精品| 在线观看www视频免费| 国产视频一区二区在线看| 亚洲中文字幕一区二区三区有码在线看 | 超碰成人久久| 91九色精品人成在线观看| 91九色精品人成在线观看| 韩国av一区二区三区四区| 精品一区二区三区av网在线观看| 欧美中文综合在线视频| 亚洲欧美激情综合另类| 午夜激情福利司机影院| 中文字幕精品亚洲无线码一区 | 99热这里只有精品一区 | 人妻丰满熟妇av一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 成人国语在线视频| 欧美成人一区二区免费高清观看 | 好男人在线观看高清免费视频 | 亚洲自偷自拍图片 自拍| 国产片内射在线| 老司机福利观看| 国产极品粉嫩免费观看在线| 日韩有码中文字幕| 最近在线观看免费完整版| 18禁裸乳无遮挡免费网站照片 | 给我免费播放毛片高清在线观看| 亚洲黑人精品在线| 满18在线观看网站| 精品免费久久久久久久清纯| 叶爱在线成人免费视频播放| 亚洲熟妇中文字幕五十中出| 欧美激情极品国产一区二区三区| 中文字幕人妻丝袜一区二区| 国产精品,欧美在线| 色综合亚洲欧美另类图片| 一区福利在线观看| 久久精品国产清高在天天线| 亚洲熟妇中文字幕五十中出| 国产成年人精品一区二区| 国产一区二区在线av高清观看| 淫秽高清视频在线观看| 国产亚洲精品久久久久5区| 国产伦人伦偷精品视频| 90打野战视频偷拍视频| 成人亚洲精品一区在线观看| 久久久精品欧美日韩精品| 成人免费观看视频高清| 国产精品爽爽va在线观看网站 | 色精品久久人妻99蜜桃| 搡老妇女老女人老熟妇| 亚洲精品国产精品久久久不卡| 一二三四社区在线视频社区8| 国产一区二区三区在线臀色熟女| 免费观看人在逋| 欧美成人性av电影在线观看| 嫁个100分男人电影在线观看| 丝袜人妻中文字幕| 黄色片一级片一级黄色片| 久久草成人影院| 欧美久久黑人一区二区| 欧美乱码精品一区二区三区| 黄片播放在线免费| 久久久久国产一级毛片高清牌| 男人操女人黄网站| 在线观看免费午夜福利视频| www国产在线视频色| 国产精品久久电影中文字幕| 中文字幕久久专区| 又黄又粗又硬又大视频| 欧美激情久久久久久爽电影| 亚洲五月色婷婷综合| 亚洲人成网站在线播放欧美日韩| 国产精品综合久久久久久久免费| 一区二区三区高清视频在线| 一级a爱片免费观看的视频| 国产一卡二卡三卡精品| 成人免费观看视频高清| cao死你这个sao货| 1024手机看黄色片| 女人被狂操c到高潮| 美女扒开内裤让男人捅视频| xxx96com| 色播在线永久视频| 亚洲精品一卡2卡三卡4卡5卡| 国产午夜福利久久久久久| 蜜臀久久99精品久久宅男| 久久亚洲精品不卡| 搡女人真爽免费视频火全软件 | av在线天堂中文字幕| 免费搜索国产男女视频| 国产精华一区二区三区| 亚洲成a人片在线一区二区| 九九爱精品视频在线观看| 欧美最黄视频在线播放免费| 小蜜桃在线观看免费完整版高清| 成年版毛片免费区| 麻豆久久精品国产亚洲av| 自拍偷自拍亚洲精品老妇| 久久99热这里只有精品18| 国产精品人妻久久久影院| or卡值多少钱| 久久韩国三级中文字幕| 一边摸一边抽搐一进一小说| 午夜激情欧美在线| 在线播放国产精品三级| 国产又黄又爽又无遮挡在线| 免费观看精品视频网站| 亚洲最大成人中文| 人人妻人人澡人人爽人人夜夜 | 久久久a久久爽久久v久久| 成人一区二区视频在线观看| 一个人看的www免费观看视频| 亚洲精品日韩av片在线观看| 香蕉av资源在线| 久久人人爽人人爽人人片va| 我要搜黄色片| 性色avwww在线观看| 少妇人妻精品综合一区二区 | 少妇的逼好多水| 国产高清有码在线观看视频| 成人av在线播放网站| 九九爱精品视频在线观看| a级一级毛片免费在线观看| 变态另类成人亚洲欧美熟女| 亚洲精品日韩在线中文字幕 | 成人毛片a级毛片在线播放| 男人舔女人下体高潮全视频| 99热网站在线观看| 久久久精品欧美日韩精品| 日韩亚洲欧美综合| 亚洲国产精品久久男人天堂| 国产在线男女| 网址你懂的国产日韩在线| 亚洲欧美精品自产自拍| 亚洲美女黄片视频| 99热只有精品国产| а√天堂www在线а√下载| 乱系列少妇在线播放| av黄色大香蕉| 亚洲av中文字字幕乱码综合| 白带黄色成豆腐渣| 欧美+亚洲+日韩+国产| 少妇被粗大猛烈的视频| 一级毛片电影观看 | 午夜激情欧美在线| 欧洲精品卡2卡3卡4卡5卡区| 小蜜桃在线观看免费完整版高清| 小说图片视频综合网站| 男人舔奶头视频| 亚洲精品色激情综合| 色播亚洲综合网| a级毛色黄片| 丝袜喷水一区| 波野结衣二区三区在线| 欧美日本亚洲视频在线播放| 十八禁网站免费在线| 女人十人毛片免费观看3o分钟| 久99久视频精品免费| 尤物成人国产欧美一区二区三区| 亚洲无线观看免费| 91午夜精品亚洲一区二区三区| 国产精品,欧美在线| 亚洲aⅴ乱码一区二区在线播放| 久99久视频精品免费| 亚洲天堂国产精品一区在线| 桃色一区二区三区在线观看| 老司机福利观看| 国产激情偷乱视频一区二区| 亚洲欧美精品自产自拍| 亚洲美女视频黄频| 日韩欧美免费精品| 亚洲国产精品成人久久小说 | 亚洲国产精品久久男人天堂| 国内精品一区二区在线观看| 99国产极品粉嫩在线观看| 亚洲国产欧洲综合997久久,| 老司机影院成人| 老师上课跳d突然被开到最大视频| 国产精品乱码一区二三区的特点| 亚洲精品在线观看二区| 久久久久久久午夜电影| 能在线免费观看的黄片| 日韩大尺度精品在线看网址| 国产精品久久久久久久电影| 久久久久久伊人网av| 人人妻人人澡欧美一区二区| 亚洲精品一卡2卡三卡4卡5卡| 3wmmmm亚洲av在线观看| 美女xxoo啪啪120秒动态图| 久久久国产成人免费| 久久久精品大字幕| 看十八女毛片水多多多| 亚洲欧美日韩无卡精品| 我要看日韩黄色一级片| 51国产日韩欧美| 国产在视频线在精品| 精品福利观看| 91av网一区二区| 偷拍熟女少妇极品色| 国产真实乱freesex| 午夜福利成人在线免费观看| 国产伦一二天堂av在线观看| 国产精品久久久久久亚洲av鲁大| 男女边吃奶边做爰视频| 人人妻人人澡欧美一区二区| 国产高清有码在线观看视频| 99久久精品国产国产毛片| 午夜激情福利司机影院| 中文字幕av在线有码专区| 久久久久久九九精品二区国产| 国产精品久久久久久精品电影| 天天躁日日操中文字幕| 高清毛片免费看| 熟妇人妻久久中文字幕3abv| 国产精品福利在线免费观看| 少妇裸体淫交视频免费看高清| 久久久久性生活片| 精品久久久久久久人妻蜜臀av| 日韩一本色道免费dvd| 国产一区亚洲一区在线观看| 一个人观看的视频www高清免费观看| 性欧美人与动物交配| 国产精品99久久久久久久久| 日本与韩国留学比较| 国产精品精品国产色婷婷| 亚洲成人久久爱视频| 波多野结衣巨乳人妻| 久久欧美精品欧美久久欧美| 国产av不卡久久| 久久久国产成人免费| 波多野结衣巨乳人妻| av中文乱码字幕在线| 深爱激情五月婷婷| 校园春色视频在线观看| 三级经典国产精品| 久久精品国产清高在天天线| 日韩高清综合在线| 中文在线观看免费www的网站| 九九在线视频观看精品| 日韩欧美一区二区三区在线观看| 色播亚洲综合网| 尤物成人国产欧美一区二区三区| 免费观看在线日韩| 天堂动漫精品| 午夜激情福利司机影院| 欧美激情国产日韩精品一区| 村上凉子中文字幕在线| 久久久精品欧美日韩精品| 久久精品国产99精品国产亚洲性色| 亚洲,欧美,日韩| 深夜a级毛片| 女人十人毛片免费观看3o分钟| 亚洲精品456在线播放app| 少妇高潮的动态图| 午夜精品一区二区三区免费看| 成人特级黄色片久久久久久久| 国产一区亚洲一区在线观看| 可以在线观看的亚洲视频| 看片在线看免费视频| 婷婷精品国产亚洲av在线| 天美传媒精品一区二区| 欧美成人免费av一区二区三区| 国产精品不卡视频一区二区| 婷婷亚洲欧美| 国产高清激情床上av| 日日摸夜夜添夜夜爱| 啦啦啦啦在线视频资源| 亚洲一区高清亚洲精品| 午夜视频国产福利| 亚洲国产精品合色在线| 亚洲成人精品中文字幕电影| 午夜亚洲福利在线播放| 成熟少妇高潮喷水视频| 国语自产精品视频在线第100页| 国产精品免费一区二区三区在线| 亚洲最大成人手机在线| 男女那种视频在线观看| 一区二区三区高清视频在线| 久久人人爽人人爽人人片va| 久久精品夜色国产| 国产精品人妻久久久久久| 欧美日韩在线观看h| 99国产精品一区二区蜜桃av| 尤物成人国产欧美一区二区三区| 我的老师免费观看完整版| 国产一区二区激情短视频| 成人鲁丝片一二三区免费| 国语自产精品视频在线第100页| 国产精品久久久久久精品电影| 亚洲最大成人手机在线| 久久久欧美国产精品| 日韩一本色道免费dvd| 亚洲七黄色美女视频| 一卡2卡三卡四卡精品乱码亚洲| av福利片在线观看| 国产精品人妻久久久影院| 日本欧美国产在线视频| 国产精品电影一区二区三区| 99国产极品粉嫩在线观看| 一夜夜www| 国内久久婷婷六月综合欲色啪| 成年女人看的毛片在线观看| 国产一区二区亚洲精品在线观看| 成人二区视频| 看十八女毛片水多多多| 五月玫瑰六月丁香| 又爽又黄无遮挡网站| 亚洲专区国产一区二区| 乱码一卡2卡4卡精品| 亚洲成人av在线免费| 97在线视频观看| 精品99又大又爽又粗少妇毛片| 插逼视频在线观看| 色哟哟·www| 精品午夜福利视频在线观看一区| 在线观看免费视频日本深夜| 亚洲av免费在线观看| 国产淫片久久久久久久久| 国产精品不卡视频一区二区| www日本黄色视频网| 性色avwww在线观看| 老女人水多毛片| 亚洲真实伦在线观看| 国产黄色小视频在线观看| 免费黄网站久久成人精品| 97超视频在线观看视频| 搡老妇女老女人老熟妇| 国产在线精品亚洲第一网站| 91午夜精品亚洲一区二区三区| av视频在线观看入口| 淫妇啪啪啪对白视频| 国产三级中文精品| 91在线精品国自产拍蜜月| 99热这里只有是精品50| 男人舔奶头视频| 亚洲婷婷狠狠爱综合网| 黄色欧美视频在线观看| 97人妻精品一区二区三区麻豆| 高清毛片免费观看视频网站| a级一级毛片免费在线观看| 色噜噜av男人的天堂激情| 麻豆国产av国片精品| 又黄又爽又免费观看的视频| 国产高清视频在线观看网站| 91精品国产九色| 精品一区二区三区视频在线观看免费| 一级毛片久久久久久久久女| 亚洲欧美日韩高清在线视频| 五月伊人婷婷丁香| 99久国产av精品国产电影| 欧美激情国产日韩精品一区| 婷婷六月久久综合丁香| 国产亚洲av嫩草精品影院| 久久鲁丝午夜福利片| 精品久久国产蜜桃| 老司机午夜福利在线观看视频| 中文亚洲av片在线观看爽| 偷拍熟女少妇极品色| 国产av一区在线观看免费| 国产男靠女视频免费网站| 高清毛片免费观看视频网站| 亚洲精华国产精华液的使用体验 | 91久久精品电影网| 女人十人毛片免费观看3o分钟| 99热精品在线国产| 久久亚洲精品不卡| 最近的中文字幕免费完整| 黄色一级大片看看| 精品久久国产蜜桃| 蜜桃久久精品国产亚洲av| 性色avwww在线观看| 亚洲精品一区av在线观看| 欧美激情在线99| 午夜免费男女啪啪视频观看 | 亚洲国产精品成人综合色| 欧美一级a爱片免费观看看| 搡老岳熟女国产| 别揉我奶头~嗯~啊~动态视频| 又黄又爽又免费观看的视频| 久久久久久久午夜电影| 少妇高潮的动态图| av黄色大香蕉| 国产男靠女视频免费网站| 亚洲va在线va天堂va国产| 日韩欧美在线乱码| 亚洲国产欧洲综合997久久,| 中出人妻视频一区二区| 永久网站在线| 丝袜美腿在线中文| 一级毛片aaaaaa免费看小| 亚洲欧美日韩无卡精品| 给我免费播放毛片高清在线观看| 久久鲁丝午夜福利片| 一级毛片电影观看 | 亚洲天堂国产精品一区在线| 日本与韩国留学比较| 日韩亚洲欧美综合| 婷婷精品国产亚洲av| 99九九线精品视频在线观看视频| av在线老鸭窝| 久久亚洲精品不卡| 干丝袜人妻中文字幕| 国内少妇人妻偷人精品xxx网站| 国内精品宾馆在线| 精品福利观看| 噜噜噜噜噜久久久久久91| 国产精品电影一区二区三区| 热99在线观看视频| 三级毛片av免费| 麻豆久久精品国产亚洲av| 六月丁香七月| 免费看光身美女| 18禁在线播放成人免费| 小蜜桃在线观看免费完整版高清| www.色视频.com| 免费黄网站久久成人精品| 国产精品久久久久久亚洲av鲁大| 一本一本综合久久| 小说图片视频综合网站| 免费人成视频x8x8入口观看| 欧美人与善性xxx| 嫩草影院精品99| 日韩一区二区视频免费看| 久久久久久久久久黄片| 欧美成人精品欧美一级黄| 此物有八面人人有两片| 精品一区二区三区视频在线| 亚洲色图av天堂| 嫩草影院新地址| 日韩在线高清观看一区二区三区| 一进一出好大好爽视频| 成年女人毛片免费观看观看9| 神马国产精品三级电影在线观看| 最近在线观看免费完整版| 寂寞人妻少妇视频99o| 春色校园在线视频观看| 国产黄片美女视频| 搡老妇女老女人老熟妇| 欧美日本视频| 亚洲国产欧美人成| 成人一区二区视频在线观看| 国产精品国产三级国产av玫瑰| 直男gayav资源| 亚洲欧美成人综合另类久久久 | 91av网一区二区| 久久久国产成人精品二区| 神马国产精品三级电影在线观看| av女优亚洲男人天堂| 精品一区二区三区视频在线观看免费| 插阴视频在线观看视频| 亚洲专区国产一区二区| 国模一区二区三区四区视频| 波野结衣二区三区在线| 18禁在线无遮挡免费观看视频 | 亚洲av熟女| 两性午夜刺激爽爽歪歪视频在线观看| 一级毛片电影观看 | 国产 一区精品| 国产单亲对白刺激| 美女大奶头视频| 久久草成人影院| 少妇被粗大猛烈的视频| 蜜桃久久精品国产亚洲av| 精华霜和精华液先用哪个| 蜜桃亚洲精品一区二区三区| 欧美又色又爽又黄视频| 99热6这里只有精品| 在线观看av片永久免费下载| 亚洲欧美日韩高清专用| 日本-黄色视频高清免费观看| 91久久精品国产一区二区成人| 天堂av国产一区二区熟女人妻| 日韩制服骚丝袜av| 欧美国产日韩亚洲一区| 日韩欧美精品v在线| 蜜桃亚洲精品一区二区三区| 亚洲精品亚洲一区二区| 丝袜美腿在线中文| 欧美人与善性xxx| 亚洲最大成人av| 国产在线精品亚洲第一网站| www日本黄色视频网| 成熟少妇高潮喷水视频| 一区福利在线观看| 欧美国产日韩亚洲一区| 日韩欧美一区二区三区在线观看| 尤物成人国产欧美一区二区三区| 国产精品一及| 丝袜喷水一区| 少妇的逼好多水| 看黄色毛片网站| 国产精品精品国产色婷婷| 国产午夜精品论理片| 天美传媒精品一区二区| 人妻制服诱惑在线中文字幕| 99国产精品一区二区蜜桃av| 婷婷亚洲欧美| av卡一久久| 中文亚洲av片在线观看爽| 久久久国产成人免费| 日韩精品青青久久久久久| 成人鲁丝片一二三区免费| 亚洲成人中文字幕在线播放| 老司机福利观看| 久久婷婷人人爽人人干人人爱| 亚洲天堂国产精品一区在线| 亚洲18禁久久av| av专区在线播放| 国产精品爽爽va在线观看网站| 国产伦在线观看视频一区| 日本a在线网址| 午夜日韩欧美国产| 国产av不卡久久| 可以在线观看毛片的网站| 日本a在线网址| 看片在线看免费视频|