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

    低頭壩養(yǎng)魚對上下游局域棲息地和魚類群落時空格局的影響*

    2016-02-23 10:21:04國家水體污染控制與治理科技重大專項項目2012ZX07501002003安徽省自然科學(xué)基金項目1508085QD75安徽省教學(xué)研究項目2013jyxm150安慶師范學(xué)院博士科研啟動基金項目044K05000130032和安徽省高校自然科學(xué)基金重點項目KJ2015A222聯(lián)合資助20141010收稿20150623
    湖泊科學(xué) 2016年1期
    關(guān)鍵詞:養(yǎng)魚

    國家水體污染控制與治理科技重大專項項目(2012ZX07501002-003)、安徽省自然科學(xué)基金項目(1508085QD75)、安徽省教學(xué)研究項目(2013jyxm150)、安慶師范學(xué)院博士科研啟動基金項目(044-K05000130032)和安徽省高校自然科學(xué)基金重點項目(KJ2015A222)聯(lián)合資助. 2014-10-10收稿;2015-06-23

    收修改稿. 萬安(1980~), 男, 博士研究生, 副教授;E-mail:wanan@aqtc.edu.cn。

    萬 安1,2,3,張曉可3,謝 楓2,3,韓 旭2,3,鄧 潔2,3,鐘 明1,李 寧1*通信作者;E-mail:Lining196@126.com.,安樹青1

    ?

    低頭壩養(yǎng)魚對上下游局域棲息地和魚類群落時空格局的影響*

    *國家水體污染控制與治理科技重大專項項目(2012ZX07501002-003)、安徽省自然科學(xué)基金項目(1508085QD75)、安徽省教學(xué)研究項目(2013jyxm150)、安慶師范學(xué)院博士科研啟動基金項目(044-K05000130032)和安徽省高校自然科學(xué)基金重點項目(KJ2015A222)聯(lián)合資助. 2014-10-10收稿;2015-06-23

    收修改稿. 萬安(1980~), 男, 博士研究生, 副教授;E-mail:wanan@aqtc.edu.cn。

    萬安1,2,3,張曉可3,謝楓2,3,韓旭2,3,鄧潔2,3,鐘明1,李寧1*通信作者;E-mail:Lining196@126.com.,安樹青1

    (1:南京大學(xué)生命科學(xué)學(xué)院,南京 210093)

    (2:安慶師范學(xué)院資源環(huán)境學(xué)院,安慶246000)

    (3:安慶師范學(xué)院水生生物保護(hù)與水生態(tài)修復(fù)安徽省高等學(xué)校工程技術(shù)研究中心,安慶246000)

    摘要:河流魚類的群落特征會隨棲息地環(huán)境的變化而發(fā)生變化.為了解魚類群落對低頭壩養(yǎng)魚的響應(yīng)程度,2013-2014年共6次調(diào)查了烏龍河(長江下游支流)兩個可涉水河段樣點(上游壩和下游壩)的環(huán)境因子、魚類多樣性及其群落結(jié)構(gòu).共采集到5目10科21屬24種魚類,以鯉形目鯉科最多,占總數(shù)的58.3%.通過解析時空動態(tài)變化與局域棲息地條件對魚類群落的影響,發(fā)現(xiàn)季節(jié)因素對魚類群落結(jié)構(gòu)無顯著影響,而魚類群落所在空間位置(上游壩,下游壩)對群落結(jié)構(gòu)差異具有較顯著影響,其中上游中華鳑鲏(Rhodeus sinensis)、 (Hemiculter leucisculus)、棒花魚(Abbottina rivularis)、大鰭鱊(Acheilognathus macropterus),下游食蚊魚(Gambusia affinis)的多度決定了這種差異;而3個棲息地變量(溶解氧濃度、電導(dǎo)率和流速)均對魚類群落特征具有顯著性影響.棲息地環(huán)境差異對河流上、下游魚類群落的影響較大,具體體現(xiàn)在下游樣點魚類的捕獲重量、捕獲數(shù)量、密度、多樣性明顯下降,推測利用低頭壩養(yǎng)魚是造成這些現(xiàn)象的重要原因。

    關(guān)鍵詞:低頭壩; 養(yǎng)魚;局域棲息地;魚類群落特征;時空格局;烏龍河

    魚類作為河流生態(tài)系統(tǒng)中的高等動物類群,它們的出生、生長、繁殖和死亡等生活史事件都是在河流中完成的,河流的生境特征對魚類多樣性及其群落結(jié)構(gòu)有著重要影響.同時由于魚類對河流生態(tài)系統(tǒng)的食物網(wǎng)動態(tài)、營養(yǎng)和物質(zhì)循環(huán)、生態(tài)系統(tǒng)的恢復(fù)能力等具有重要的調(diào)節(jié)作用,因此魚類多樣性對于維持穩(wěn)定的河流生態(tài)系統(tǒng)結(jié)構(gòu)和功能有著十分重要的作用.自然河流的環(huán)境呈現(xiàn)的是一種連續(xù)的線性梯度變化過程,河流魚類的物種分布、豐富度和群落結(jié)構(gòu)會隨環(huán)境梯度的變化而變化[1],但當(dāng)這種環(huán)境梯度被中斷(如水利工程[2])時,河流自然流態(tài)機(jī)制將隨之而改變,使魚類難以正常完成生活史循環(huán);此外受農(nóng)業(yè)灌溉、水體污染[3-4]、土地利用[5]等的影響,河流魚類多樣性和物種豐富度也正在逐漸下降和減少,魚類群落的空間格局也隨之改變[6-8].通過觀測魚類群落空間格局的變化,可以了解河流生態(tài)系統(tǒng)的變化情況[9]。

    圖1 研究樣點位置Fig.1 Location of study area

    在諸多影響魚類群落變化的因素中,水壩無疑是最普遍、最典型的一種人為干擾形式[10].水壩又可分為高頭壩(high-head dam)和低頭壩(low-head dam),其本質(zhì)區(qū)別在于前者壩頭較高并能完全截流,后者壩頭較低因而不能完全截流[11].目前,國內(nèi)外對水壩生態(tài)學(xué)效應(yīng)的研究大多數(shù)聚焦于大型水壩,僅有少量研究有關(guān)低頭壩對河流魚類群落結(jié)構(gòu)及其多樣性的影響,如在魚類多樣性上,Tiemann等發(fā)現(xiàn)緊鄰水壩局部河段中的魚類多度、密度和均勻度受到低頭壩的顯著影響,而魚類物種豐富度無顯著變化[12];Dodd等則發(fā)現(xiàn),壩下沖刷區(qū)的物種豐富度顯著高于其它河段[13].從魚類群落結(jié)構(gòu)水平看,雖然低頭壩顯著影響了少數(shù)物種的優(yōu)勢度,而魚類群落結(jié)構(gòu)卻沒有顯著變化[13-14];但Gillette等[15]、Poulet[16]和Yan等[17]則觀察到與遠(yuǎn)離水壩的河段相比,蓄水區(qū)的急流性魚類顯著減少且緩流性或靜水性魚類顯著增多;而Yan等更進(jìn)一步發(fā)現(xiàn)低頭壩減弱了魚類群落在河流網(wǎng)絡(luò)中的空間自相關(guān)性[18],并且與本地入侵種(native invaders)魚類的多樣性增加有較大相關(guān)[19].上述有關(guān)低頭壩對魚類群落影響的研究結(jié)果之間尚存很大爭議和分歧[15-17, 20-21],推測低頭壩產(chǎn)生的生態(tài)學(xué)效應(yīng)很可能與水壩自身的屬性、人為干擾程度以及數(shù)量累積效應(yīng)密切相關(guān)[11,18,22]。

    為此,本研究以某低頭壩攔河養(yǎng)殖的河段為對象,連續(xù)2年調(diào)查了其上、下游的壩下沖刷區(qū)(plunge area)的魚類群落特征,旨在確定該河段魚類群落的組成與時空格局動態(tài)變化,解析用于養(yǎng)殖目的的低頭壩對上下游魚類群落的影響,從而為區(qū)域內(nèi)魚類資源的保護(hù)和河流的科學(xué)管理提供參考依據(jù),也為生態(tài)水利學(xué)的相關(guān)研究積累資料。

    1 材料與方法

    1.1 研究區(qū)域概況

    研究樣點位于烏龍河,該河匯入石門湖后經(jīng)皖河入長江,屬亞熱帶長江中下游地區(qū)的小型河流,該河段已修筑2座低頭壩攔河形成人工湖泊(圖1所示雙龍湖),用于養(yǎng)殖,常年飼喂鰱、鳙、鯽,湖中也有少量野生、麥穗魚、鳑鲏和鲇等.其上游壩高約0.9 m,下游壩高約1.0 m,兩壩之間無其他阻隔,相距1.2 km.研究于2013年4、7、11月和2014年4、7、10月分6次對上游壩和下游壩的沖刷區(qū)共進(jìn)行12樣次的采集工作(圖1).研究河段周邊地區(qū)實行封閉管理,絕無酷漁濫捕等現(xiàn)象。

    1.2 采樣方法

    采集工具為背負(fù)式超聲波捕魚器(CWB-5000 P),采取“之”字形捕撈路線,一人電魚,兩人以撈網(wǎng)(網(wǎng)目為0.2 cm)尾隨捕撈,采樣河道長30 m,采樣時間20 min,每樣點均采用相同捕撈努力.為避免小生境差異對魚類群落的影響,樣點盡可能包括該河段的各類生境(岸邊、河中央、靜水、流水).所獲標(biāo)本在新鮮狀態(tài)下進(jìn)行常規(guī)鑒別、稱重,疑難種以10%甲醛溶液固定后,帶回實驗室鑒別,其余活體全部放歸[23]。

    1.3 現(xiàn)場環(huán)境因子觀察與測量

    現(xiàn)場觀察水色、氣味和沿岸植被情況。

    在每個樣點的3個等距截面上測量水寬(WW)、透明度(SD),于各截面的3個等距點分別測量水深(WD)、水溫(Tem)、溶解氧(DO,雷磁JPB-607A型)、電導(dǎo)率(Cond,雷磁DDBJ-350型)和 pH值(雷磁PHB-4型),水深60%處測量流速(Speed,榮圣MT-A型).所得數(shù)據(jù)全部取平均值,作為該樣次的環(huán)境因子數(shù)據(jù)。

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

    樣點的物種多樣性指數(shù)采用Shannon-Wiener指數(shù)(H′):

    H′=-∑Pi·lnPi

    (1)

    相關(guān)的魚類多樣性指標(biāo)密度(D)、頻率(F)、相對多度(P)和重要值指數(shù)(IVI)[24-25]的計算公式分別為:

    Di=Ni/A

    (2)

    Fi=100%(Si/S)

    (3)

    Pi=100%(Ni/N)

    (4)

    IVIi=104Fi·Pi

    (5)

    式中,Si和Ni分別代表物種i的采到次數(shù)和累計個體數(shù),S、N和A分別代表全部采樣次數(shù)、全部魚類個體數(shù)和采樣的水面積。

    利用t-test檢驗各樣點的多樣性、密度和捕獲量數(shù)據(jù)的差異顯著性。

    為提高數(shù)據(jù)的正態(tài)性和方差齊性,降低極端數(shù)值的負(fù)面影響,將相關(guān)魚類數(shù)據(jù)和棲息地變量數(shù)據(jù)進(jìn)行對數(shù)lg(x+1)轉(zhuǎn)化,在PRIMER 5.0 軟件中,運(yùn)用單因素相似性分析(one-way crossed ANOSIM),分別檢驗空間因素和時間因素對魚類群落結(jié)構(gòu)的影響;運(yùn)用非參數(shù)多變量排序(nonmetric multidimensional scaling,NMS)方法,構(gòu)建魚類群落結(jié)構(gòu)的二維雙標(biāo)圖;最后,運(yùn)用相似性百分比(SIMPER)分析,解析形成群落差異的關(guān)鍵貢獻(xiàn)物種[19]。

    利用CANOCO 4.5軟件的典范對應(yīng)分析(CCA)解析棲息地變量對魚類群落結(jié)構(gòu)的影響,此外對稀有種(出現(xiàn)樣點數(shù)不超過2次的物種)不予以典范對應(yīng)分析,以降低稀有物種的權(quán)重[23,26]。

    3 結(jié)果與分析

    3.1 棲息地環(huán)境觀察

    調(diào)查期間內(nèi),上游壩水色為澄清,無味,岸邊有菰(Zizanialatifolia)、菖蒲(Acoruscalamus)、蘆葦(Phragmitesaustralis),河岸及河中央有苦草(Vallisnerianatans)、金魚藻(Ceratophyllumdemersum)、狐尾藻(Myriophyllum)等水下植物;下游壩水色為深綠色,有臭味,岸邊有狗牙根(Cynodondactylon)、水花生(Alternantheraphiloxeroides)、水蓼(Polygonumhydropiper),河岸水下有黑藻(Hydrilla),河中央無沉水植物。

    3.2 物種組成與分布

    12樣次共捕獲標(biāo)本626尾,隸屬5目10科21屬24種.以鯉形目最多(15種),占總數(shù)的62.5%,鱸形目有5種,占總數(shù)的20.8%;全部10科中,鯉科魚類最多(14種),分布于6個亞科,占總數(shù)的58.3%,其次是塘鱧科魚類2種,占8.3%,其他各科魚類總計占33.4%;所捕獲的24種魚類都屬于定居型物種,無洄游性或半洄游性物種。

    下游壩常見種為鯽、泥鰍、食蚊魚3種(F>40%),麥穗魚為偶見種(10%100) (表1),均為靜水性魚類[23]。

    表1 烏龍河魚類分布及其頻率、相對多度和重要值指數(shù)

    3.3 上下游魚類群落差異

    統(tǒng)計2年各樣點數(shù)據(jù),發(fā)現(xiàn)下游壩魚類的捕獲重量、捕獲數(shù)量、密度和多樣性(Shannon-Wiener指數(shù))水平低于上游壩的水平;同時上游壩每年秋季(10、11月)的捕獲量和密度均為最大.而下游壩在7月捕獲量和密度均為最大(圖2)。

    圖2 各樣次捕獲重量、捕獲數(shù)量、密度和多樣性的比較(*表示P<0.01)Fig.2 Comparisons of the catch weigh, catch number, density and diversity among the samples(*stands for P<0.01)

    3.4 魚類群落結(jié)構(gòu)時空動態(tài)變化與排序

    圖3 魚類群落結(jié)構(gòu)的空間變化非度量多維標(biāo)度排序Fig.3 Spatial variations of fish assemblages based on non-metric multi-dimensional scaling

    運(yùn)用ANOSIM檢驗空間因素(樣點)和時間因素(季節(jié))對魚類群落結(jié)構(gòu)的影響,結(jié)果表明,在整體上,魚類群落結(jié)構(gòu)無顯著的季節(jié)變化(Global R=-0.45,P<0.05)和顯著的空間變化(Global R=1,P<0.05).而非參數(shù)多變量排序(NMS)也表明上游壩和下游壩的魚類群落結(jié)構(gòu)明顯不同,聚為兩大區(qū)域(圖3)。

    由于上游壩和下游壩間的魚類群落結(jié)構(gòu)存在顯著差異,運(yùn)用相似性百分比(SIMPER)分別解析維持上游壩和下游壩魚類群落內(nèi)不相似性的關(guān)鍵貢獻(xiàn)物種.結(jié)果顯示造成上游壩和下游壩群落結(jié)構(gòu)差異的主要物種為中華鰟鲏、、棒花魚、大鰭鱊、馬口魚、寬鰭鱲、興凱鱊、飄魚、子陵吻鰕虎魚、細(xì)鱗鲴、食蚊魚、銀鮈、沙塘鱧、刺鰍、麥穗魚(累積貢獻(xiàn)率>90%);其中上游壩樣點關(guān)鍵物種是中華鰟鲏、、棒花魚,大鰭鱊 (貢獻(xiàn)值>7%),下游壩樣點為食蚊魚(表2)。

    3.5 魚類物種與棲息地環(huán)境的關(guān)系

    CCA分析發(fā)現(xiàn)3個棲息地變量(溶解氧、電導(dǎo)率和流速)對魚類群落結(jié)構(gòu)具有顯著影響(P<0.05).溶解氧、電導(dǎo)率和流速均主要與軸1相關(guān),其中溶解氧、流速與軸1、軸2均呈負(fù)相關(guān),而電導(dǎo)率與軸1、軸2均呈正相關(guān)(圖4)。

    表2 關(guān)鍵物種對上、下游樣點魚類群落結(jié)構(gòu)差異的貢獻(xiàn)

    圖4 CCA所得魚類群落結(jié)構(gòu)與環(huán)境因素的關(guān)系圖Fig.4 Correlations between fish assemblages and habitat variables after CCA

    4 討論

    水壩是人為干擾河流常見形式之一,由于水壩的修建使得河流流速減緩、大面積的靜水水體出現(xiàn)在水壩上、下游,并導(dǎo)致下游水量減少;河段魚類產(chǎn)卵場縮小,產(chǎn)漂流性卵魚類發(fā)育受限,繁殖季節(jié)滯后.此外,水壩的建立對河流生態(tài)系統(tǒng)中的理化環(huán)境與生物成分也產(chǎn)生較大的影響,由于壩上蓄水區(qū)(impoundment area)流速變緩,水體對于污染物的稀釋、擴(kuò)散、遷移和凈化能力將下降.壩下流水的流速趨于平穩(wěn)或者減小,水位大幅度下降甚至斷流,河底泥沙淤積,減弱了水體的擴(kuò)散稀釋污染物的能力,影響了河流的自凈功能[27],并因此導(dǎo)致河流水質(zhì)的惡化.特別是在人為干擾較大地區(qū),水壩除了導(dǎo)致河流自然流態(tài)的變化外,由于其所承載的社會服務(wù)功能較山區(qū)水壩多,其魚類群落特征變化也更為明顯.近年來各地濕地資源日益枯竭,湖泊面積逐漸減少[28],自然漁業(yè)資源受到嚴(yán)重威脅,淡水養(yǎng)殖業(yè)開始日益發(fā)展.而《中華人民共和國漁業(yè)法》第十條也規(guī)定“國家鼓勵全民所有制單位、集體所有制單位和個人充分利用適于養(yǎng)殖的水域、灘涂發(fā)展養(yǎng)殖業(yè)”.很多地方開始利用河流進(jìn)行漁業(yè)養(yǎng)殖,但是養(yǎng)殖過程產(chǎn)生的殘餌、化肥以及未經(jīng)處理就直接排放的農(nóng)藥易引起水質(zhì)中的氮和磷超標(biāo),滋生病原;特別是河流養(yǎng)殖中,大量增加的人工設(shè)施會使下游流速降低,影響了河流營養(yǎng)物質(zhì)的輸入和輸出,使陸源污染物得不到及時的稀釋擴(kuò)散,從而導(dǎo)致河流水質(zhì)受到不同程度的污染,并有繼續(xù)發(fā)展的趨勢[23].人們也意識到隨著養(yǎng)殖年限增加,養(yǎng)殖水體的富營養(yǎng)化程度也逐漸加重[29],導(dǎo)致養(yǎng)殖水體中污染物蓄積,從而嚴(yán)重危害到下游魚類的生存[30].因此在諸多河流養(yǎng)殖形式中,利用低頭壩攔截河流進(jìn)行養(yǎng)魚,會對河段上下游自然魚類群落造成較大影響。

    魚類群落的動態(tài)變化受到內(nèi)源性和外源性因素聯(lián)合影響[31].內(nèi)源性因素一般指魚類自身繁殖或洄游導(dǎo)致的種群數(shù)量和群落結(jié)構(gòu)的變化[32];外源性因素主要是引起棲息地穩(wěn)定性及其有效性的顯著性季節(jié)動態(tài),從而導(dǎo)致溪流魚類群落的物種組成及數(shù)量發(fā)生相應(yīng)的變化[33-34].Yan等通過對青弋江3座低頭壩的調(diào)查發(fā)現(xiàn)低頭壩對魚類的影響主要體現(xiàn)在:1) 與上下游1 km的參考點相比,壩上蓄水區(qū)的喜靜水魚類增加導(dǎo)致物種多樣性和密度增加;2) 壩下沖刷區(qū)魚類組成情況則與上下游1 km處參考點的物種組成和群落結(jié)構(gòu)相差不大;3) 導(dǎo)致各樣點魚類群落差異的主要原因不是季節(jié)因素,而是低頭壩造成水文條件改變,其中水寬、水深、流速以及底質(zhì)類型對魚類群落有較大影響[35]。

    本研究中所得結(jié)果表明,通過ANOSIM檢驗認(rèn)為魚類群落結(jié)構(gòu)無顯著性的季節(jié)變化(Global R=-0.45,P<0.05),但具有顯著的空間變化(Global R=1,P<0.05),所以養(yǎng)殖河段上下游的魚類群落差異也主要由空間位置(上游壩、下游壩)引起,因此低頭壩的作用明顯(圖3)。

    此外,上游壩的魚類數(shù)量相對多度排列順序為:中華鳑鲏(27.81%)>(10.71%)>鯽(9.67%)>棒花(9.5%),下游壩為鯽(57.14%)>食蚊魚(23.81%)>泥鰍(16.67%)>麥穗魚(2.38%).與表1中所示順序不同,相似性百分比(SIMPER)顯示上游壩的中華鰟鲏、、棒花魚、大鰭鱊,下游壩的食蚊魚對上、下游魚類群落的空間變化貢獻(xiàn)值較大。

    由于魚類物種組成的空間變化與物種的棲息地選擇密切聯(lián)系[36],上游壩中的大鰭鱊(產(chǎn)卵需較高的水體溶解氧濃度)以及下游壩中的食蚊魚(喜食耐污、好低氧的孑孓)在上、下游壩樣點的平均多度較高(表2),提示下游壩較上游壩的水體溶解氧濃度低,污染程度大[23,36],是導(dǎo)致比食蚊魚和大鰭鱊相對多度大的鯽(表1)卻對上、下游的魚類群落空間差異的貢獻(xiàn)較小的原因。

    同樣,CCA分析結(jié)果表明,該河段中魚類(除食蚊魚和鯽外)基本集中在軸1和軸2交叉點的左部(圖4),表明這些物種生境需求基本類似,種間競爭激烈,且生態(tài)幅較小,無論環(huán)境中流速、溶解氧濃度增加和減少,均不利于這些魚類生存,說明該河段中的魚類群落對外界環(huán)境的變化也應(yīng)較為敏感,生境較脆弱[35];而魚類分布差異主要受溶解氧濃度、電導(dǎo)率和流速3個環(huán)境因素影響,特別是食蚊魚和鯽,這兩種較耐受不良環(huán)境的魚,遠(yuǎn)離其他魚類聚集的區(qū)域,喜溫度較高、電導(dǎo)率較高的環(huán)境(圖4);在魚類組成種類上,上游壩共有23種魚,其中寬鰭鱲、馬口魚、亮銀鮈、銀鮈、沙塘鱧5種為喜流性魚類,其余18種為靜水性魚類;而下游壩的魚類群落僅有麥穗魚、鯽、泥鰍和食蚊魚,并且3種優(yōu)勢種均為耐污生命力強(qiáng)的靜水性魚類[23].這與Yan等發(fā)現(xiàn)的壩下沖刷區(qū)與參考區(qū)域的魚類(喜靜水、喜流水均有)群落組成相差不大的結(jié)論有明顯的不同[35].以上結(jié)果表明上下壩沖刷區(qū)環(huán)境條件的改變造成了各樣點間不同的魚類群落組成特點[19,35]。

    同時發(fā)現(xiàn)在上游樣點各月份的捕獲重量、捕獲數(shù)量、密度和多樣性均高于下游樣點,且上游樣點自身秋季(10、11月)的捕獲重量、捕獲數(shù)量、密度和多樣性高于4月(圖2).一般認(rèn)為亞熱帶地區(qū)的平水期(4月)和枯水期(10、11月)恰為魚類的繁殖前期和非繁殖期,而冬季普遍溫度低且降雨量極小,非常不利于魚類的生存,因此經(jīng)過越冬后,4月份魚類種群數(shù)量、生物量均出現(xiàn)下降,而秋季(10、11月)因為有當(dāng)年生的幼體補(bǔ)充,魚類種群數(shù)量和密度反而增加(圖2).這與在同一地區(qū)的青弋江、秋浦河觀察到的結(jié)果相似[23,36].同時與在秋浦河觀察的結(jié)果一致[36],相對于4月,上游壩采樣地7月份的漁獲物數(shù)量和密度降低(圖2),而多樣性指數(shù)則有降(2013年)有升(2014年)(圖2),這很可能與豐水期間雨量充沛,魚類棲息地面積增加,電捕器的捕獲效率下降,導(dǎo)致偶然誤差較多有關(guān).同時觀察到在下游采樣點7月份的捕獲重量、捕獲數(shù)量、密度和多樣性指數(shù)均較其他月份出現(xiàn)增高的趨勢(圖2),這很可能與下游壩在夏季水量較大,水體環(huán)境得到改善,下游魚類上溯有關(guān)。

    總之,與上游壩沖刷區(qū)相比,用于養(yǎng)殖目的的低頭壩對下游壩沖刷區(qū)的局域棲息地及魚類群落結(jié)構(gòu)具有較顯著的影響,其中3個棲息地變量(溶解氧、電導(dǎo)率和流速)均對魚類群落特征具有顯著影響,并且出現(xiàn)魚類物種多樣性、捕獲數(shù)量和重量、密度顯著減少等特點,而這些變化與季節(jié)無顯著相關(guān)性.推測下游低溶解氧濃度、高污染是造成上游中華鰟鲏、、棒花魚、大鰭鱊和下游食蚊魚對魚類群落空間差異貢獻(xiàn)較大的主要原因. 因此,今后在探討低頭壩對魚類分類群影響時,不能僅就其蓄水區(qū)和沖刷區(qū)魚類分類群進(jìn)行分析,還需對低頭壩自身用途以及魚類功能群(棲息地、營養(yǎng)、繁殖)組成等因素加以考慮。

    當(dāng)前,隨著我國經(jīng)濟(jì)建設(shè)的快速發(fā)展,河流所承擔(dān)的服務(wù)功能將越來越多,由此導(dǎo)致的一系列的棲息地退化,水體污染和入侵種威脅的危險也隨之增加[37-38].因此,本文所得結(jié)果將有助于人們了解水利工程對河流生態(tài)系統(tǒng)的影響,為科學(xué)利用河流與河流生態(tài)系統(tǒng)健康、水生生物的保護(hù)等相關(guān)研究提供科學(xué)參考。

    致謝:衷心感謝原中國科學(xué)院南京地理與湖泊研究所朱松泉研究員親自指點魚類鑒別工作,安慶師范學(xué)院生命科學(xué)學(xué)院陶峰、李家磊同學(xué)的辛苦采樣工作。

    參考文獻(xiàn)5

    [1]Vannote RL, Minshall GW, Cumins KWetal. The river continuum concept.CanadianJournalofFisheriesandAquaticSciences, 1980, 37(1): 130-137。

    [2]March JG, Benstead JP, Pringle CMetal. Damming tropical island streams: problems, solutions, and alternatives.BioSicence, 2003, 53(11): 1069-1078。

    [3]Allan JD. Landscapes and riverscapes: the influence of land use on stream ecosystems.AnnualReviewofEcology,Evolution&Systematics, 2004, 35: 257-284. 。

    [4]Lima-Junior SE, Cardone IB, Goitein R. Fish assemblage structure and aquatic pollution in a Brazilian stream: some limitations of diversity indices and models for environmental impact studies.EcologyofFreshwaterFish, 2006, 15(3): 284-290。

    [5]Harding JS, Benfield EF, Bolstad PVetal. Stream biodiversity: the ghost of land use past.ProceedingsoftheNationalAcademyofSciencesoftheUnitedStatesofAmerica, 1998, 95(25): 14843-14847。

    [6]Bonner TH, Wilde GR. Changes in the Canadian River fish assemblage associated with reservoir construction.JournalofFreshwaterEcology, 2000, 15(2): 189-198。

    [7]Meldgaard T, Nielsen EE, Loeschcke V. Fragmentation by weirs in a riverine system: a study of genetic variation in time and space among populations of European grayling (Thymallusthymallus) in a Danish river system.ConservationGenetics, 2003, 4(6): 735-747。

    [8]Galat DL, Berry CR, Garaner WMetal. Spatiotemporal patterns and changes in Missouri River fishes.AmericanFisheriesSocietySymposium, 2005, 45: 249-291。

    [9]Karr JR, Fausch KD, Angermeier PLetal. Assessing biological integrity in running waters: A method and its rationale. Champaign: Illinois Natural History Survey(Special Publication), 1986: 1-31。

    [10]Rosenberg DM, Berks F, Boadaly RAetal. Large-scale impacts of hydroelectric development.EnvironmentalReviews, 1997, 5(1): 27-54。

    [11]Poff NL, Hart DD. How dams vary and why it matters for the emerging science of dam removal.BioScience, 2002, 52(8): 659-668。

    [12]Tiemann JS, Gillette DP, Wildhaber MLetal. Effects of lowhead dams on riffle-dwelling fishes and macro invertebrates in a Midwestern river.TransactionsofAmericanFisheriesSociety, 2004, 133(3): 705-717。

    [13]Dodd HR, Hayes DB, Baylis JRetal. Low-head sea lamprey barrier effects on stream habitat and fish communities in Great Lakes basin.JournalofGreatLakesResearch, 2003, 29(Suppl. 1): 386-402。

    [14]Raborn SW, Schramm SW. Fish assemblage response to recent mitigation of a channelized warmwater stream.RiverResearchandApplications, 2003, 19(4): 289-301。

    [15]Gillette DP, Tiemann JS, Edds DRetal. Spatiotemporal patterns of fish assemblage structure in a river impoundment by low-head dams.Copeia, 2005, 3: 539-549。

    [16]Poulet N. Impact of weirs on fish communities in a piedmont stream.RiverResearchandApplications, 2007, 23(9): 1038-1047。

    [17]Yan YZ, Wang H, Zhu Retal. Influences of low-head dams on the fish assemblages in the headwater streams of the Qingyi watershed, China.EnvironmentalBiologyofFishes, 2013, 96(4): 495-506。

    [18]Yan YZ. Influences of low-head dams on fish assemblages in the headwater streams of the Qingyi River, China. 144th Annual Meeting of the American Fisheries Society. Afs, 2014。

    [19]Chu L, Wang W, Zhu Retal. Variation in fish assemblages across impoundments of low-head dams in headwater streams of the Qingyi River, China: effects of abiotic factors and native invaders.EnvironmentalBiologyofFishes, 2014, 98(1): 101-112。

    [20]Cummings GS. The impact of low-head dams on fish species richness in Wisconsin, USA.EcologicalApplications, 2004, 14(5): 1495-1506。

    [21]Helms BS, Werneke DC, Gangloff MMetal. The influence of low-head dams on fish assemblages in streams across Alabama.JournaloftheNorthAmericanBenthologicalSociety, 2011, 30(4): 1095-1106。

    [22]Ward JV, Stanford JA. The serial discontinuity concept of lotic ecosystems. In: Fontaine TD, Bartell SM eds. Dynamics of lotic ecosystems, Ann Arbor Science, Ann Arbor, MI, 1983: 29-42。

    [23]王紅. 低頭壩對青弋江河源溪流局域棲息地, 魚類多樣性及其群落結(jié)構(gòu)的影響[學(xué)位論文]. 合肥: 安徽師范大學(xué), 2012。

    [24]祁繼英, 阮曉紅. 大壩對河流生態(tài)系統(tǒng)的環(huán)境影響分析. 河海大學(xué)學(xué)報: 自然科學(xué)版, 2005, 33(1): 37-40。

    [25]鄭奕, 劉新春, 楊金龍等. 塔里木河上游地區(qū)阿拉爾段天然退化生態(tài)系統(tǒng)植物群落物種多樣性特征分析. 水土保持研究, 2005, 12(2): 131-134。

    [26]周勁風(fēng), 溫琰茂, 梁志謙. 珠江三角洲密養(yǎng)池塘的水質(zhì)特點分析. 水產(chǎn)科學(xué), 2005, 23(12): 13-17。

    [27]陳生熬, 童華文, 霍佳佳等. 塔里木河流域淡水養(yǎng)殖業(yè)可持續(xù)發(fā)展思考. 塔里木大學(xué)學(xué)報, 2008, 19(4): 103-106。

    [28]安樹青, 李哈濱, 關(guān)保華等. 中國的天然濕地: 過去的問題, 現(xiàn)狀和未來的挑戰(zhàn). AMBIO-人類環(huán)境雜志, 2007, 36(4): 317-324。

    [29]鄧朝陽, 朱仁, 嚴(yán)云志. 長江蕪湖江段魚類多樣性及其群落結(jié)構(gòu)的時空格局. 淡水漁業(yè), 2013, 43(1): 28-36。

    [30]Gauch HG. Multivariate analysis in community ecology. New York: Cambridge University Press, 1982: 23。

    [31]Belyea LR, Lancaster J. Assembly rules within a contingent ecology.Oikos, 1999, 86(3): 402-416。

    [32]Welcomme RL. The biology and ecology of the fishes of a small tropical stream.JournalofZoology, 1969, (4): 485-529。

    [33]Grossman GD, Moyle PB, Whitaker Jr JR. Stochasticity in structural and functional characteristics of an Indiana stream fish assemblage: a test of community theory.TheAmericanNaturalist, 1982, 120(4): 423-454

    [34]Matthews WJ. Fish faunal “breaks” and stream order in the eastern and central United States.EnvironmentalBiologyofFishes, 1986, 17(2): 81-92。

    [35]Yan YZ, Xiang XY, Chu Letal. Influences of local habitat and stream spatial position on fish assemblages in a dammed watershed, the Qingyi Stream, China.EcologyofFreshwaterFish, 2011, 20(2): 199-208。

    [36]王文劍, 儲玲, 司春等. 秋浦河源國家濕地公園溪流魚類群落的時空格局. 動物學(xué)研究, 2013, 34(4): 417-428。

    [37]Fu CZ, Wu JH, Chen JKetal. Freshwater fish biodiversity in the Yangtze River basin of China: patterns, threats and conservation.BiodiversityandConservation, 2003, 12(8): 1649-1685。

    [38]Kang B, He DM, Perrett Letal. Fish and fisheries in the Upper Mekong: current assessment of the fish community, threat and conservation.ReviewsinFishBiologyandFisheries, 2009, 19(4): 465-480。

    J.LakeSci.(湖泊科學(xué)), 2016, 28(1): 178-186

    ?2016 byJournalofLakeSciences

    Low-head-dam fish culture effects on spatial-temporal patterns of local habitat and fish assemblages in the upstream and downstream of rivers

    WAN An1,2,3, ZHANG Xiaoke3, XIE Feng2,3, HAN Xu2,3, DENG Jie2,3, ZHONG Ming1, LI Ning1**& AN Shuqing1***

    (1:SchoolofLifeScience,NanjingUniversity,Nanjing210093,P.R.China)

    (2:SchoolofResourcesandEnvironment,AnqingNormalUniversity,Anqing246000,P.R.China)

    (3:ResearchCenterofAquaticOrganismConservationandWaterEcosystemRestorationofAnhuiEducationDepartment,AnqingNormalUniversity,Anqing246000,P.R.China)

    Abstract:Fish assemblage will change along with the river habitat environment. In order to find the response degree of fish assemblage to fish culture within low-head-dam water, two low-head dams (including the upstream and downstream dams) of the Wulong River (a tributary of the lower reach of the Yangtze River, China) were selected as sampling sites, where the habitat variables, the fish diversity and fish assemblage structure were surveyed for 6 times from 2013 to 2014 in these wadeable reaches of the sampling sites. A total of 24 fish species which belong to 5 orders, 21 genera and 24 species were caught, with dominant Cyprinidae accounting for 58.3% of the total. After an analysis on the influences of temporal changes and local habitat conditions upon the fish assemblage, seasonal factors were found to have no apparent effect on the fish assemblage structure. The spatial location (the upstream dam or the downstream dam) affected the assemblage’s structural differences significantly, which were determined by the abundance of Hemiculter leucisculus, Abbottina rivularis, Rhodeus sinensis, Acheilognathus macropterus in the upstream site, and the abundance of Gambusia affinis in downstream site. Besides, all the three habitat variables including dissolved oxygen, conductivity and current velocity exerted a great impact on the fish assemblage. The differences in habitat environment bring about some great influences on the fish assemblage, showing in the dramatic decrease of capture weight, quantity, density and diversity of fish in the sampling sites of the downstream dam, which are probably caused by the fish culture within the low-head dam。

    Keywords:Low-head dam; fish culture; local habitat; fish assemblage characteristics; spatial-temporal pattern; Wulong River

    通信作者*共同;E-mail:anshq@nju.edu.cn。

    DOI10.18307/2016.0121

    猜你喜歡
    養(yǎng)魚
    養(yǎng)魚記
    為什么自來水不能養(yǎng)魚?
    養(yǎng)魚記
    魚經(jīng)
    銅梁區(qū)建成首個移動式微循環(huán)流水養(yǎng)魚槽項目
    電腦迷(2019年3期)2019-01-05 11:09:33
    二氫吡啶在養(yǎng)魚生產(chǎn)上的應(yīng)用
    廣東飼料(2016年6期)2016-12-01 03:43:29
    利用顆粒飼料養(yǎng)魚好處多
    蓮田養(yǎng)魚效益好
    養(yǎng)魚
    怎樣避免網(wǎng)箱養(yǎng)魚患“營養(yǎng)不良”癥
    乱系列少妇在线播放| 99热网站在线观看| 精品国产三级普通话版| 波多野结衣高清无吗| 久久久久久九九精品二区国产| 看十八女毛片水多多多| 久久久久国产精品人妻aⅴ院| 国产探花极品一区二区| 成人欧美大片| 亚洲三级黄色毛片| 久久久精品大字幕| 久久久精品大字幕| 国产高清视频在线播放一区| 国产精品日韩av在线免费观看| 午夜福利视频1000在线观看| 国产日本99.免费观看| 国产高清三级在线| 欧美最新免费一区二区三区| 日日撸夜夜添| 久久九九热精品免费| 国产精品永久免费网站| 女同久久另类99精品国产91| 日韩欧美三级三区| 18禁裸乳无遮挡免费网站照片| 中文字幕免费在线视频6| 搡老熟女国产l中国老女人| 国产爱豆传媒在线观看| 十八禁网站免费在线| 性色avwww在线观看| 亚洲男人的天堂狠狠| 看免费成人av毛片| 两个人的视频大全免费| 亚洲,欧美,日韩| 亚洲av日韩精品久久久久久密| 狂野欧美激情性xxxx在线观看| 久久6这里有精品| 男人舔女人下体高潮全视频| 真人做人爱边吃奶动态| 欧美成人免费av一区二区三区| 婷婷丁香在线五月| 99热这里只有精品一区| 老司机福利观看| 啦啦啦啦在线视频资源| 欧美3d第一页| 国产久久久一区二区三区| 国产大屁股一区二区在线视频| 亚洲三级黄色毛片| 搡老熟女国产l中国老女人| 老熟妇乱子伦视频在线观看| 精品一区二区免费观看| 免费观看的影片在线观看| 久久久久精品国产欧美久久久| 成人国产综合亚洲| 亚洲综合色惰| 成年人黄色毛片网站| 最近最新中文字幕大全电影3| 变态另类成人亚洲欧美熟女| 97超视频在线观看视频| 在线免费十八禁| 啦啦啦观看免费观看视频高清| 999久久久精品免费观看国产| 美女黄网站色视频| 国产亚洲av嫩草精品影院| 亚洲天堂国产精品一区在线| av在线蜜桃| 国产亚洲91精品色在线| 51国产日韩欧美| 99久久精品一区二区三区| 一区二区三区高清视频在线| 亚洲狠狠婷婷综合久久图片| 99久国产av精品| 国产探花极品一区二区| 久久久久久伊人网av| videossex国产| 女的被弄到高潮叫床怎么办 | 丰满乱子伦码专区| 亚洲成a人片在线一区二区| 窝窝影院91人妻| 窝窝影院91人妻| 国产私拍福利视频在线观看| 男女啪啪激烈高潮av片| 不卡一级毛片| 久久国产乱子免费精品| 欧美一区二区国产精品久久精品| 美女cb高潮喷水在线观看| 一卡2卡三卡四卡精品乱码亚洲| 麻豆一二三区av精品| 欧美潮喷喷水| 午夜视频国产福利| 亚洲欧美日韩东京热| 全区人妻精品视频| 国产精品久久久久久亚洲av鲁大| av国产免费在线观看| 91麻豆精品激情在线观看国产| 国产精品三级大全| 99热只有精品国产| 少妇被粗大猛烈的视频| 国产 一区精品| 一a级毛片在线观看| 欧美在线一区亚洲| 亚洲在线自拍视频| 亚洲精品在线观看二区| 色播亚洲综合网| 欧美极品一区二区三区四区| 伦理电影大哥的女人| 久久这里只有精品中国| 欧美色视频一区免费| 窝窝影院91人妻| 午夜福利在线在线| 亚洲av熟女| 亚洲国产精品久久男人天堂| 国产精品电影一区二区三区| 免费在线观看日本一区| 成人午夜高清在线视频| 久久午夜福利片| 国产欧美日韩一区二区精品| 国产一区二区激情短视频| 久久久久国内视频| 国产人妻一区二区三区在| 欧美日韩亚洲国产一区二区在线观看| 亚洲三级黄色毛片| 99久久九九国产精品国产免费| 国内少妇人妻偷人精品xxx网站| 男女做爰动态图高潮gif福利片| 淫秽高清视频在线观看| 国产亚洲精品久久久com| 蜜桃亚洲精品一区二区三区| 全区人妻精品视频| 日韩欧美免费精品| 国产 一区 欧美 日韩| 国产一级毛片七仙女欲春2| 精品久久久久久久久av| 女的被弄到高潮叫床怎么办 | av视频在线观看入口| 欧美国产日韩亚洲一区| 久久久午夜欧美精品| 制服丝袜大香蕉在线| 久久精品国产鲁丝片午夜精品 | 免费在线观看成人毛片| 成人高潮视频无遮挡免费网站| 久久亚洲精品不卡| 久99久视频精品免费| 午夜日韩欧美国产| 国产91精品成人一区二区三区| 欧美色视频一区免费| 91在线观看av| 国产私拍福利视频在线观看| 啦啦啦啦在线视频资源| 人妻夜夜爽99麻豆av| 深爱激情五月婷婷| 别揉我奶头 嗯啊视频| 免费av不卡在线播放| 他把我摸到了高潮在线观看| 简卡轻食公司| 成人永久免费在线观看视频| 免费av毛片视频| 男人的好看免费观看在线视频| 日本一本二区三区精品| 亚洲无线观看免费| 一个人免费在线观看电影| 18禁黄网站禁片免费观看直播| 久久人人精品亚洲av| 丰满的人妻完整版| 超碰av人人做人人爽久久| 日韩亚洲欧美综合| 免费黄网站久久成人精品| 校园春色视频在线观看| 91麻豆av在线| 久久精品国产亚洲av天美| 韩国av一区二区三区四区| 国产精品一区二区性色av| 国产一区二区在线观看日韩| aaaaa片日本免费| 亚洲av日韩精品久久久久久密| 欧美成人性av电影在线观看| 欧美成人a在线观看| 亚洲av不卡在线观看| 窝窝影院91人妻| 久久精品91蜜桃| 搡女人真爽免费视频火全软件 | 极品教师在线免费播放| 国产麻豆成人av免费视频| 欧美最黄视频在线播放免费| 男女啪啪激烈高潮av片| 久久久久久大精品| 欧美成人免费av一区二区三区| 亚洲午夜理论影院| 国产91精品成人一区二区三区| 在现免费观看毛片| 搡老熟女国产l中国老女人| 国产一区二区在线观看日韩| 最近在线观看免费完整版| 中文字幕久久专区| 男人舔女人下体高潮全视频| 午夜激情欧美在线| 国产高清不卡午夜福利| 欧美性猛交黑人性爽| 91精品国产九色| 变态另类成人亚洲欧美熟女| 亚洲av成人av| 人妻久久中文字幕网| 国产亚洲精品久久久久久毛片| 99热只有精品国产| 免费看日本二区| 亚洲专区国产一区二区| 国产精品精品国产色婷婷| 欧美日韩亚洲国产一区二区在线观看| 三级国产精品欧美在线观看| 一区二区三区激情视频| 69av精品久久久久久| 波多野结衣高清作品| 精品国内亚洲2022精品成人| 精品人妻一区二区三区麻豆 | 男女下面进入的视频免费午夜| 春色校园在线视频观看| 日本熟妇午夜| 亚洲av成人av| 日本在线视频免费播放| 日本与韩国留学比较| 精品久久久久久久末码| 在线观看舔阴道视频| 九九久久精品国产亚洲av麻豆| 好男人在线观看高清免费视频| 国产激情偷乱视频一区二区| 免费av毛片视频| 亚洲成人久久爱视频| 欧美日韩瑟瑟在线播放| 久久久久久久午夜电影| 两人在一起打扑克的视频| 欧美在线一区亚洲| 欧美又色又爽又黄视频| 日日撸夜夜添| 亚洲最大成人av| 天堂动漫精品| 18禁黄网站禁片免费观看直播| 草草在线视频免费看| 日韩欧美三级三区| 久久人人爽人人爽人人片va| 国产伦人伦偷精品视频| 成人欧美大片| 久久久久久伊人网av| 国产成人av教育| 变态另类成人亚洲欧美熟女| 精品免费久久久久久久清纯| a级一级毛片免费在线观看| 日韩中字成人| 中文字幕精品亚洲无线码一区| 亚洲色图av天堂| 精品久久久久久久久av| 亚洲图色成人| 婷婷精品国产亚洲av在线| 国产精品久久视频播放| 一个人看视频在线观看www免费| 国产精品国产三级国产av玫瑰| 亚洲欧美激情综合另类| 麻豆久久精品国产亚洲av| av视频在线观看入口| 午夜精品一区二区三区免费看| 亚洲精品一区av在线观看| 国产成人影院久久av| 国产伦一二天堂av在线观看| 国产一区二区在线av高清观看| 欧美中文日本在线观看视频| 神马国产精品三级电影在线观看| 禁无遮挡网站| 国产精品1区2区在线观看.| 久久久久久大精品| 九色国产91popny在线| 少妇人妻一区二区三区视频| 国产麻豆成人av免费视频| 精品不卡国产一区二区三区| 国产一区二区三区视频了| 春色校园在线视频观看| 亚洲国产精品成人综合色| 日韩强制内射视频| 日本免费a在线| 噜噜噜噜噜久久久久久91| 最新中文字幕久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 欧美绝顶高潮抽搐喷水| 内射极品少妇av片p| 日韩一本色道免费dvd| 亚洲精品一区av在线观看| www日本黄色视频网| 听说在线观看完整版免费高清| 啦啦啦观看免费观看视频高清| 麻豆一二三区av精品| 熟妇人妻久久中文字幕3abv| 人妻夜夜爽99麻豆av| 黄片wwwwww| 在线免费观看的www视频| 亚洲精品日韩av片在线观看| 成人高潮视频无遮挡免费网站| 99视频精品全部免费 在线| 国内少妇人妻偷人精品xxx网站| 免费观看在线日韩| 国产女主播在线喷水免费视频网站 | 88av欧美| 亚洲男人的天堂狠狠| 亚洲av中文av极速乱 | 又粗又爽又猛毛片免费看| 桃红色精品国产亚洲av| 麻豆成人av在线观看| 在线观看美女被高潮喷水网站| 欧美一级a爱片免费观看看| 久久精品久久久久久噜噜老黄 | 草草在线视频免费看| 啦啦啦啦在线视频资源| 波多野结衣高清无吗| 99久久精品国产国产毛片| 亚洲色图av天堂| 免费观看人在逋| 日韩一区二区视频免费看| 精品一区二区免费观看| 成人亚洲精品av一区二区| 日韩欧美国产一区二区入口| 乱系列少妇在线播放| 亚洲电影在线观看av| 啪啪无遮挡十八禁网站| 老司机深夜福利视频在线观看| 床上黄色一级片| 尤物成人国产欧美一区二区三区| 搡老岳熟女国产| 国内精品久久久久久久电影| 悠悠久久av| 欧美xxxx黑人xx丫x性爽| 熟妇人妻久久中文字幕3abv| 国内毛片毛片毛片毛片毛片| 18+在线观看网站| 草草在线视频免费看| 女生性感内裤真人,穿戴方法视频| 国产精品一及| 亚洲人成网站在线播放欧美日韩| 日本色播在线视频| 久久精品综合一区二区三区| a在线观看视频网站| 欧洲精品卡2卡3卡4卡5卡区| 精品福利观看| 久久久久久伊人网av| a级毛片a级免费在线| 亚洲精品成人久久久久久| 美女高潮喷水抽搐中文字幕| 国产男靠女视频免费网站| 69人妻影院| 尤物成人国产欧美一区二区三区| 变态另类成人亚洲欧美熟女| 国产精品无大码| 最近中文字幕高清免费大全6 | 国产成人a区在线观看| 我的老师免费观看完整版| 我要看日韩黄色一级片| 欧美国产日韩亚洲一区| 久久婷婷人人爽人人干人人爱| 国产精品久久电影中文字幕| 色噜噜av男人的天堂激情| 91精品国产九色| 亚洲国产精品sss在线观看| 99热6这里只有精品| ponron亚洲| 久久久久九九精品影院| 亚洲中文字幕一区二区三区有码在线看| 最后的刺客免费高清国语| h日本视频在线播放| 欧美激情在线99| 两性午夜刺激爽爽歪歪视频在线观看| 狠狠狠狠99中文字幕| 搞女人的毛片| 国产伦人伦偷精品视频| 日韩,欧美,国产一区二区三区 | 日本-黄色视频高清免费观看| 在线观看午夜福利视频| 日本黄色视频三级网站网址| 亚洲最大成人手机在线| 精品一区二区三区视频在线观看免费| 成年人黄色毛片网站| 在线天堂最新版资源| videossex国产| 国语自产精品视频在线第100页| 91久久精品国产一区二区成人| 色综合色国产| 国产精品久久视频播放| 成人特级黄色片久久久久久久| 嫁个100分男人电影在线观看| 免费av不卡在线播放| 性插视频无遮挡在线免费观看| 少妇熟女aⅴ在线视频| 免费观看人在逋| 亚洲va日本ⅴa欧美va伊人久久| 少妇裸体淫交视频免费看高清| 国产午夜精品论理片| 我要搜黄色片| 成人国产综合亚洲| 亚洲精品影视一区二区三区av| 成人国产综合亚洲| videossex国产| 在线天堂最新版资源| 99riav亚洲国产免费| 成人特级黄色片久久久久久久| a在线观看视频网站| 一区二区三区激情视频| 简卡轻食公司| av天堂在线播放| 久久99热6这里只有精品| 亚洲真实伦在线观看| 国产69精品久久久久777片| 婷婷精品国产亚洲av在线| 夜夜夜夜夜久久久久| 久久人人精品亚洲av| 亚洲成a人片在线一区二区| 色在线成人网| 日本三级黄在线观看| 日韩精品有码人妻一区| 中文字幕人妻熟人妻熟丝袜美| 香蕉av资源在线| 久久国产精品人妻蜜桃| 琪琪午夜伦伦电影理论片6080| 国产激情偷乱视频一区二区| 超碰av人人做人人爽久久| 天天一区二区日本电影三级| 久久久久久久精品吃奶| 亚洲国产高清在线一区二区三| 嫁个100分男人电影在线观看| 国产一区二区三区视频了| av女优亚洲男人天堂| ponron亚洲| xxxwww97欧美| av国产免费在线观看| 婷婷丁香在线五月| 我的女老师完整版在线观看| 国内精品久久久久精免费| 精品久久久久久久久av| av在线亚洲专区| 91在线精品国自产拍蜜月| .国产精品久久| 国产亚洲精品久久久com| 三级国产精品欧美在线观看| xxxwww97欧美| 欧美bdsm另类| 一本精品99久久精品77| 亚洲内射少妇av| 小蜜桃在线观看免费完整版高清| 国产精品精品国产色婷婷| 成人一区二区视频在线观看| 五月玫瑰六月丁香| 欧美高清成人免费视频www| 国产黄色小视频在线观看| 欧美不卡视频在线免费观看| 午夜福利视频1000在线观看| 国产精品综合久久久久久久免费| 在线看三级毛片| 无遮挡黄片免费观看| 日韩欧美在线二视频| 一进一出好大好爽视频| 国产伦精品一区二区三区四那| 国产一区二区三区在线臀色熟女| 欧美不卡视频在线免费观看| 一本久久中文字幕| 国产三级中文精品| 老司机午夜福利在线观看视频| 男女做爰动态图高潮gif福利片| 老女人水多毛片| 亚洲国产欧洲综合997久久,| 色5月婷婷丁香| 亚洲av第一区精品v没综合| 日韩精品中文字幕看吧| 午夜精品在线福利| 男女啪啪激烈高潮av片| 亚洲av电影不卡..在线观看| 丰满人妻一区二区三区视频av| 成人一区二区视频在线观看| 嫩草影院新地址| 变态另类成人亚洲欧美熟女| 免费观看人在逋| 老司机深夜福利视频在线观看| videossex国产| 三级国产精品欧美在线观看| xxxwww97欧美| 久久99热这里只有精品18| 麻豆成人午夜福利视频| 日本五十路高清| 极品教师在线视频| 国产精品国产三级国产av玫瑰| 日韩一区二区视频免费看| 日本免费a在线| 亚洲精品一卡2卡三卡4卡5卡| 日韩亚洲欧美综合| 欧美日本亚洲视频在线播放| 直男gayav资源| 国产一区二区在线av高清观看| 国产伦在线观看视频一区| 嫩草影院入口| 小说图片视频综合网站| 国产精品久久视频播放| 精品午夜福利视频在线观看一区| 老司机深夜福利视频在线观看| 亚洲最大成人av| 国产老妇女一区| 国产精品国产三级国产av玫瑰| 亚洲avbb在线观看| 深夜a级毛片| 毛片女人毛片| 亚洲人成网站在线播放欧美日韩| 久久午夜福利片| x7x7x7水蜜桃| 国产aⅴ精品一区二区三区波| 国产精品福利在线免费观看| 国模一区二区三区四区视频| 国产精品98久久久久久宅男小说| 成人午夜高清在线视频| 国产伦一二天堂av在线观看| 午夜亚洲福利在线播放| 久久久久久久久久成人| 午夜亚洲福利在线播放| 天天一区二区日本电影三级| 久久久久久久精品吃奶| 国产老妇女一区| 午夜福利视频1000在线观看| 在线观看免费视频日本深夜| 很黄的视频免费| 国内揄拍国产精品人妻在线| 看片在线看免费视频| 91在线精品国自产拍蜜月| 动漫黄色视频在线观看| 免费看av在线观看网站| 美女黄网站色视频| 在线观看舔阴道视频| 亚洲精品久久国产高清桃花| 此物有八面人人有两片| 直男gayav资源| 淫秽高清视频在线观看| 一级黄色大片毛片| 18禁黄网站禁片午夜丰满| 久久亚洲精品不卡| 国产精品电影一区二区三区| 岛国在线免费视频观看| 国内精品美女久久久久久| 啦啦啦啦在线视频资源| 亚洲av美国av| 欧美+日韩+精品| 欧美成人一区二区免费高清观看| 日韩一本色道免费dvd| 国产久久久一区二区三区| 国国产精品蜜臀av免费| 亚洲色图av天堂| 亚洲五月天丁香| 成人特级黄色片久久久久久久| 日韩欧美精品v在线| 99在线视频只有这里精品首页| 国产乱人视频| 一本一本综合久久| 亚洲经典国产精华液单| 色精品久久人妻99蜜桃| 在线a可以看的网站| 国产精品久久视频播放| 九九在线视频观看精品| 日韩高清综合在线| 99视频精品全部免费 在线| 日本一二三区视频观看| 午夜影院日韩av| 淫妇啪啪啪对白视频| 天美传媒精品一区二区| 身体一侧抽搐| 亚洲人成网站在线播| 久久精品综合一区二区三区| 婷婷精品国产亚洲av| 免费在线观看成人毛片| 日日撸夜夜添| 老司机深夜福利视频在线观看| 成年人黄色毛片网站| 五月伊人婷婷丁香| av专区在线播放| 久久天躁狠狠躁夜夜2o2o| 日韩精品青青久久久久久| 色尼玛亚洲综合影院| a级毛片免费高清观看在线播放| 亚洲经典国产精华液单| 久久99热6这里只有精品| 97超视频在线观看视频| 特级一级黄色大片| 男人和女人高潮做爰伦理| 97人妻精品一区二区三区麻豆| 精品一区二区三区av网在线观看| 岛国在线免费视频观看| 久久欧美精品欧美久久欧美| 天堂影院成人在线观看| 欧美性感艳星| 夜夜看夜夜爽夜夜摸| 女的被弄到高潮叫床怎么办 | 老女人水多毛片| 少妇裸体淫交视频免费看高清| 亚洲中文字幕一区二区三区有码在线看| 精品午夜福利视频在线观看一区| 精品久久久久久久久av| 小说图片视频综合网站| 最新在线观看一区二区三区| 国产精品一区二区三区四区久久| 日本熟妇午夜| 能在线免费观看的黄片| 久久久久久伊人网av| 内射极品少妇av片p| 国产人妻一区二区三区在| 午夜a级毛片| 色综合婷婷激情| 久久国内精品自在自线图片| 少妇被粗大猛烈的视频| 久久久久久久久中文| 午夜免费激情av| 在线观看一区二区三区| 小蜜桃在线观看免费完整版高清| 禁无遮挡网站| eeuss影院久久| 亚洲国产日韩欧美精品在线观看| 禁无遮挡网站| 搡女人真爽免费视频火全软件 | 婷婷精品国产亚洲av|