• <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)不良”癥
    亚洲av五月六月丁香网| 国产精品久久久久久亚洲av鲁大| 一级毛片高清免费大全| 久久精品国产综合久久久| 亚洲国产精品久久男人天堂| 国产精品永久免费网站| 成熟少妇高潮喷水视频| 婷婷丁香在线五月| 99热这里只有精品一区| 男人和女人高潮做爰伦理| 床上黄色一级片| 18禁国产床啪视频网站| 亚洲av美国av| 人人妻人人澡欧美一区二区| 欧美三级亚洲精品| 在线观看av片永久免费下载| 国产一区二区三区视频了| ponron亚洲| 欧美日本亚洲视频在线播放| 12—13女人毛片做爰片一| 亚洲av五月六月丁香网| 韩国av一区二区三区四区| 一本一本综合久久| 1024手机看黄色片| 精品熟女少妇八av免费久了| 欧美中文日本在线观看视频| 欧美另类亚洲清纯唯美| 国产高清三级在线| 国产成+人综合+亚洲专区| 亚洲中文日韩欧美视频| 中文字幕高清在线视频| 九九在线视频观看精品| 精品人妻偷拍中文字幕| 色噜噜av男人的天堂激情| 一二三四社区在线视频社区8| 亚洲激情在线av| 成年女人永久免费观看视频| 日本一二三区视频观看| 色av中文字幕| 免费人成在线观看视频色| 国产探花极品一区二区| 麻豆久久精品国产亚洲av| 国产精品久久久久久亚洲av鲁大| 欧美bdsm另类| 久久久久久久午夜电影| 国产精品女同一区二区软件 | 最后的刺客免费高清国语| 日韩高清综合在线| 99在线视频只有这里精品首页| 一级作爱视频免费观看| 夜夜夜夜夜久久久久| 欧美另类亚洲清纯唯美| 丁香六月欧美| 欧美激情久久久久久爽电影| 精品福利观看| 免费在线观看影片大全网站| 男女下面进入的视频免费午夜| 有码 亚洲区| 青草久久国产| 婷婷丁香在线五月| 国产av麻豆久久久久久久| 老汉色∧v一级毛片| 欧美黑人欧美精品刺激| 狂野欧美白嫩少妇大欣赏| 中文字幕久久专区| 日本黄色视频三级网站网址| 成人三级黄色视频| 国产精品免费一区二区三区在线| 国产三级黄色录像| 一夜夜www| 免费看美女性在线毛片视频| 免费看十八禁软件| x7x7x7水蜜桃| 国产精品永久免费网站| 亚洲av中文字字幕乱码综合| 一个人看的www免费观看视频| 亚洲国产高清在线一区二区三| 国产一区二区三区视频了| 国产av在哪里看| 午夜免费成人在线视频| 欧美3d第一页| 啦啦啦观看免费观看视频高清| 又黄又爽又免费观看的视频| 又粗又爽又猛毛片免费看| 草草在线视频免费看| 日韩欧美一区二区三区在线观看| 久久久久亚洲av毛片大全| 超碰av人人做人人爽久久 | 日本一二三区视频观看| 在线看三级毛片| 女人被狂操c到高潮| 国产久久久一区二区三区| 看免费av毛片| 黑人欧美特级aaaaaa片| 韩国av一区二区三区四区| 久久精品亚洲精品国产色婷小说| 熟女人妻精品中文字幕| 最新在线观看一区二区三区| 在线看三级毛片| 色综合欧美亚洲国产小说| 亚洲久久久久久中文字幕| 久久久久久久亚洲中文字幕 | 国产高清三级在线| 亚洲av五月六月丁香网| 久久久久亚洲av毛片大全| 欧美午夜高清在线| 淫秽高清视频在线观看| 观看美女的网站| 真人一进一出gif抽搐免费| 精品日产1卡2卡| www.熟女人妻精品国产| 午夜日韩欧美国产| 在线a可以看的网站| 看片在线看免费视频| 午夜福利成人在线免费观看| 国产一区二区三区在线臀色熟女| 午夜免费成人在线视频| 成年版毛片免费区| 成人性生交大片免费视频hd| 亚洲熟妇熟女久久| 亚洲国产精品999在线| 国产老妇女一区| 小蜜桃在线观看免费完整版高清| 天堂影院成人在线观看| 欧美中文综合在线视频| 成人亚洲精品av一区二区| 蜜桃亚洲精品一区二区三区| 国产亚洲欧美98| 中文字幕精品亚洲无线码一区| 亚洲人与动物交配视频| 韩国av一区二区三区四区| 欧美日韩国产亚洲二区| 亚洲国产高清在线一区二区三| 久久精品夜夜夜夜夜久久蜜豆| 最近视频中文字幕2019在线8| 日韩欧美精品免费久久 | 在线观看日韩欧美| 国产在视频线在精品| 国产精品乱码一区二三区的特点| 男人舔奶头视频| 国产av不卡久久| 午夜福利视频1000在线观看| 国产精品 欧美亚洲| 亚洲无线观看免费| 国产高清视频在线播放一区| 国产精品一及| 夜夜爽天天搞| 亚洲欧美精品综合久久99| 欧美日本视频| 91久久精品电影网| 久久精品国产亚洲av香蕉五月| 久久精品影院6| 99久久无色码亚洲精品果冻| 免费在线观看成人毛片| 欧美色欧美亚洲另类二区| 国产色爽女视频免费观看| 手机成人av网站| 又黄又爽又免费观看的视频| 精品人妻一区二区三区麻豆 | 欧美黄色片欧美黄色片| 99国产综合亚洲精品| 欧美日韩福利视频一区二区| 欧美国产日韩亚洲一区| 国产伦一二天堂av在线观看| 99久久久亚洲精品蜜臀av| 国产精品久久久人人做人人爽| 免费大片18禁| 99在线视频只有这里精品首页| 在线看三级毛片| 亚洲中文日韩欧美视频| 久久久久国产精品人妻aⅴ院| 亚洲午夜理论影院| 国产探花在线观看一区二区| 观看免费一级毛片| xxx96com| 一级毛片女人18水好多| 真实男女啪啪啪动态图| 成人亚洲精品av一区二区| 看免费av毛片| 在线观看66精品国产| 最好的美女福利视频网| 久久久久久大精品| 99久久精品国产亚洲精品| 国产精品久久久久久久久免 | 国产亚洲精品一区二区www| 亚洲男人的天堂狠狠| 久久久久性生活片| 成年人黄色毛片网站| 母亲3免费完整高清在线观看| 国产一区二区在线av高清观看| av国产免费在线观看| 12—13女人毛片做爰片一| 男人舔女人下体高潮全视频| 我要搜黄色片| 国产日本99.免费观看| 精品国内亚洲2022精品成人| 一区二区三区激情视频| 一二三四社区在线视频社区8| 国产中年淑女户外野战色| 欧美绝顶高潮抽搐喷水| 亚洲18禁久久av| av片东京热男人的天堂| 看片在线看免费视频| 一区二区三区高清视频在线| 国产主播在线观看一区二区| 国产国拍精品亚洲av在线观看 | 欧美在线一区亚洲| 91字幕亚洲| 人人妻人人澡欧美一区二区| 最近视频中文字幕2019在线8| 1024手机看黄色片| 欧美最黄视频在线播放免费| 欧美日韩精品网址| 成人av一区二区三区在线看| 九色成人免费人妻av| 在线观看午夜福利视频| 久久久久性生活片| 麻豆一二三区av精品| 两个人的视频大全免费| 深夜精品福利| 毛片女人毛片| 丁香六月欧美| 精品乱码久久久久久99久播| 99热这里只有精品一区| 97碰自拍视频| 国产精华一区二区三区| 国产av麻豆久久久久久久| 在线十欧美十亚洲十日本专区| 老司机午夜十八禁免费视频| 制服人妻中文乱码| 一级毛片女人18水好多| 成人av一区二区三区在线看| 国产91精品成人一区二区三区| 欧美成人a在线观看| 99久久精品国产亚洲精品| 看片在线看免费视频| 中文字幕av成人在线电影| 久久精品综合一区二区三区| 麻豆国产av国片精品| 亚洲成av人片在线播放无| 婷婷六月久久综合丁香| 亚洲人成电影免费在线| 国产视频一区二区在线看| 噜噜噜噜噜久久久久久91| www.熟女人妻精品国产| 中文字幕人成人乱码亚洲影| 亚洲人与动物交配视频| 亚洲中文字幕一区二区三区有码在线看| 午夜亚洲福利在线播放| 亚洲av中文字字幕乱码综合| 五月玫瑰六月丁香| 日韩欧美国产在线观看| 亚洲自拍偷在线| 久久午夜亚洲精品久久| 国语自产精品视频在线第100页| 国产精品日韩av在线免费观看| 成年女人毛片免费观看观看9| 不卡一级毛片| 天天躁日日操中文字幕| 无限看片的www在线观看| 久久精品影院6| 欧洲精品卡2卡3卡4卡5卡区| 亚洲自拍偷在线| 麻豆国产av国片精品| 亚洲av五月六月丁香网| 久久久久久久久久黄片| 在线观看日韩欧美| 国产乱人伦免费视频| 露出奶头的视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | aaaaa片日本免费| 少妇的逼水好多| 此物有八面人人有两片| 亚洲av熟女| 法律面前人人平等表现在哪些方面| 亚洲精品一卡2卡三卡4卡5卡| 久久99热这里只有精品18| 99riav亚洲国产免费| 国产又黄又爽又无遮挡在线| 成人国产一区最新在线观看| 亚洲国产精品成人综合色| 尤物成人国产欧美一区二区三区| 日韩人妻高清精品专区| 最近视频中文字幕2019在线8| 欧美乱妇无乱码| 日本撒尿小便嘘嘘汇集6| 日本黄色片子视频| 国内久久婷婷六月综合欲色啪| 可以在线观看毛片的网站| 岛国在线观看网站| 欧美bdsm另类| 国产精品久久久久久人妻精品电影| 亚洲av不卡在线观看| 国产探花极品一区二区| 久久中文看片网| 成人国产综合亚洲| 久久精品国产自在天天线| 一个人观看的视频www高清免费观看| 熟女少妇亚洲综合色aaa.| 亚洲无线观看免费| 国产男靠女视频免费网站| 69人妻影院| 国产aⅴ精品一区二区三区波| 麻豆久久精品国产亚洲av| 欧美性感艳星| 国产精品永久免费网站| 国产爱豆传媒在线观看| 欧美黑人巨大hd| 欧美+日韩+精品| 琪琪午夜伦伦电影理论片6080| 亚洲在线自拍视频| www.999成人在线观看| 在线播放国产精品三级| 国产精品久久久人人做人人爽| 成熟少妇高潮喷水视频| 少妇丰满av| 又爽又黄无遮挡网站| av天堂在线播放| 精品欧美国产一区二区三| 内射极品少妇av片p| 18禁国产床啪视频网站| 亚洲欧美精品综合久久99| 亚洲av成人不卡在线观看播放网| 在线观看日韩欧美| ponron亚洲| 精品欧美国产一区二区三| 亚洲成人免费电影在线观看| 久久天躁狠狠躁夜夜2o2o| 国产精品乱码一区二三区的特点| 两个人视频免费观看高清| 人人妻人人看人人澡| 国产成人福利小说| 成人精品一区二区免费| 中文字幕熟女人妻在线| 在线播放无遮挡| 亚洲国产欧美人成| 一级毛片高清免费大全| 天天添夜夜摸| 18+在线观看网站| 国产三级黄色录像| 色综合婷婷激情| 久久久久久久午夜电影| 欧美高清成人免费视频www| 亚洲熟妇熟女久久| 波多野结衣巨乳人妻| 国产精品三级大全| 人人妻人人看人人澡| 久久久久久久午夜电影| netflix在线观看网站| 久久久久久国产a免费观看| 成人一区二区视频在线观看| 久久久久久国产a免费观看| 可以在线观看毛片的网站| 国产私拍福利视频在线观看| 在线免费观看不下载黄p国产 | 制服丝袜大香蕉在线| 中文字幕av在线有码专区| 久久婷婷人人爽人人干人人爱| 国产欧美日韩精品亚洲av| 国产精品免费一区二区三区在线| 色精品久久人妻99蜜桃| 天堂av国产一区二区熟女人妻| 亚洲精品成人久久久久久| 久久九九热精品免费| 亚洲av成人精品一区久久| 午夜影院日韩av| 啪啪无遮挡十八禁网站| 欧美国产日韩亚洲一区| av天堂在线播放| 久久久久久久久大av| 最近视频中文字幕2019在线8| 欧美日韩福利视频一区二区| 香蕉久久夜色| 午夜免费观看网址| 老司机午夜十八禁免费视频| 日本黄色视频三级网站网址| 天天添夜夜摸| 亚洲在线自拍视频| 亚洲欧美日韩无卡精品| 中国美女看黄片| 国产熟女xx| 99精品欧美一区二区三区四区| 少妇的丰满在线观看| 人妻夜夜爽99麻豆av| 久久亚洲真实| 久久久久免费精品人妻一区二区| 国产一区二区在线观看日韩 | 黄色片一级片一级黄色片| 免费av不卡在线播放| 俄罗斯特黄特色一大片| 日日夜夜操网爽| 一区二区三区免费毛片| 国产高清三级在线| 亚洲美女黄片视频| 亚洲欧美激情综合另类| ponron亚洲| 国产成人福利小说| 国产高清videossex| 老司机在亚洲福利影院| 乱人视频在线观看| 欧美最黄视频在线播放免费| 午夜免费观看网址| 内地一区二区视频在线| 波多野结衣高清无吗| 亚洲成人久久爱视频| 欧美三级亚洲精品| 亚洲精品影视一区二区三区av| 一个人看视频在线观看www免费 | 特大巨黑吊av在线直播| 精品电影一区二区在线| 给我免费播放毛片高清在线观看| 九九在线视频观看精品| 极品教师在线免费播放| 国产伦精品一区二区三区四那| 精品国产美女av久久久久小说| 亚洲国产中文字幕在线视频| 亚洲成人精品中文字幕电影| 亚洲成a人片在线一区二区| 成年女人看的毛片在线观看| 91在线精品国自产拍蜜月 | or卡值多少钱| 男女那种视频在线观看| 91字幕亚洲| 天天添夜夜摸| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲精品乱码久久久v下载方式 | 亚洲七黄色美女视频| 欧美中文综合在线视频| 亚洲人成网站在线播放欧美日韩| 国产精品精品国产色婷婷| 看免费av毛片| 欧美黑人欧美精品刺激| 亚洲av中文字字幕乱码综合| 免费观看人在逋| 少妇熟女aⅴ在线视频| 中亚洲国语对白在线视频| 成人特级黄色片久久久久久久| 男人和女人高潮做爰伦理| 国产精品综合久久久久久久免费| e午夜精品久久久久久久| eeuss影院久久| 日本在线视频免费播放| 久久精品影院6| 免费av观看视频| 69人妻影院| 在线a可以看的网站| 亚洲欧美日韩高清专用| 免费看a级黄色片| 欧美极品一区二区三区四区| 久久久久久久久久黄片| 18+在线观看网站| 日本a在线网址| 757午夜福利合集在线观看| 香蕉丝袜av| 国产午夜精品论理片| 2021天堂中文幕一二区在线观| a级毛片a级免费在线| 亚洲在线自拍视频| 国产成年人精品一区二区| xxx96com| 他把我摸到了高潮在线观看| 成熟少妇高潮喷水视频| 人人妻,人人澡人人爽秒播| 国产精品女同一区二区软件 | 99热6这里只有精品| 91av网一区二区| 精品人妻一区二区三区麻豆 | 国产综合懂色| 九九热线精品视视频播放| 亚洲欧美日韩东京热| 日韩欧美精品免费久久 | 欧美一区二区国产精品久久精品| 国产伦在线观看视频一区| 啪啪无遮挡十八禁网站| 悠悠久久av| 97超级碰碰碰精品色视频在线观看| 午夜福利欧美成人| 亚洲欧美日韩高清在线视频| 最近视频中文字幕2019在线8| 少妇高潮的动态图| 给我免费播放毛片高清在线观看| 女人高潮潮喷娇喘18禁视频| 欧美国产日韩亚洲一区| 有码 亚洲区| 欧美绝顶高潮抽搐喷水| 黄片大片在线免费观看| av中文乱码字幕在线| 精品99又大又爽又粗少妇毛片 | 国产伦精品一区二区三区视频9 | 精品乱码久久久久久99久播| 亚洲午夜理论影院| 久久精品91无色码中文字幕| 亚洲 国产 在线| 成熟少妇高潮喷水视频| 一二三四社区在线视频社区8| 在线观看舔阴道视频| 日本 欧美在线| 国产高清视频在线播放一区| 婷婷丁香在线五月| 一二三四社区在线视频社区8| 国产伦一二天堂av在线观看| 俺也久久电影网| 成人国产一区最新在线观看| 蜜桃久久精品国产亚洲av| 日本精品一区二区三区蜜桃| 国产亚洲精品一区二区www| 免费看十八禁软件| 波多野结衣高清无吗| 亚洲不卡免费看| 国产一区二区在线av高清观看| 天天一区二区日本电影三级| 青草久久国产| 在线视频色国产色| 天堂网av新在线| 欧美日韩乱码在线| 999久久久精品免费观看国产| 日韩有码中文字幕| 久久精品国产自在天天线| 高清在线国产一区| 国产精品久久久久久精品电影| 尤物成人国产欧美一区二区三区| 天天一区二区日本电影三级| 五月玫瑰六月丁香| xxx96com| 88av欧美| 99精品在免费线老司机午夜| 欧美一区二区国产精品久久精品| 国产欧美日韩精品一区二区| 香蕉av资源在线| 国产又黄又爽又无遮挡在线| 久久精品国产清高在天天线| 精品人妻一区二区三区麻豆 | 色在线成人网| 又爽又黄无遮挡网站| 91av网一区二区| 精品日产1卡2卡| 丁香欧美五月| 校园春色视频在线观看| 亚洲欧美日韩无卡精品| 亚洲精品456在线播放app | 国产乱人视频| 日本三级黄在线观看| 色综合站精品国产| h日本视频在线播放| 18禁黄网站禁片午夜丰满| 国内揄拍国产精品人妻在线| 波多野结衣巨乳人妻| 中亚洲国语对白在线视频| 国产三级中文精品| 天堂影院成人在线观看| 香蕉丝袜av| 亚洲一区二区三区色噜噜| 亚洲精品456在线播放app | 成年免费大片在线观看| 91麻豆精品激情在线观看国产| 亚洲国产欧美网| 国产成人影院久久av| 国产午夜精品久久久久久一区二区三区 | 亚洲国产精品成人综合色| 久久久久久九九精品二区国产| 亚洲国产中文字幕在线视频| 90打野战视频偷拍视频| 在线观看66精品国产| 两人在一起打扑克的视频| 国产精品,欧美在线| 又黄又粗又硬又大视频| 国产精品嫩草影院av在线观看 | 亚洲国产精品合色在线| 国产精品美女特级片免费视频播放器| 国产高清激情床上av| 色综合婷婷激情| 嫁个100分男人电影在线观看| 国产午夜精品久久久久久一区二区三区 | 精品欧美国产一区二区三| 丁香欧美五月| 丁香六月欧美| 免费看十八禁软件| 亚洲成人精品中文字幕电影| 看免费av毛片| 国产精品 国内视频| 亚洲av成人不卡在线观看播放网| 久久久国产成人精品二区| 热99在线观看视频| 亚洲成av人片在线播放无| 青草久久国产| 一a级毛片在线观看| 久久久久久久久中文| 黄片大片在线免费观看| www日本黄色视频网| 国产av麻豆久久久久久久| 18美女黄网站色大片免费观看| 久久人人精品亚洲av| 脱女人内裤的视频| 国产精品 国内视频| 2021天堂中文幕一二区在线观| 国产色婷婷99| 亚洲黑人精品在线| 免费看a级黄色片| 一级毛片女人18水好多| 亚洲精品影视一区二区三区av| 最好的美女福利视频网| 国产探花在线观看一区二区| 色精品久久人妻99蜜桃| 久久国产乱子伦精品免费另类| 可以在线观看毛片的网站| 欧美成人一区二区免费高清观看| 日韩欧美国产一区二区入口| 国产综合懂色| 欧美激情在线99| 女人高潮潮喷娇喘18禁视频| 啦啦啦观看免费观看视频高清| 亚洲精华国产精华精| 在线视频色国产色| 真人做人爱边吃奶动态| 禁无遮挡网站|