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

    以鰱、鳙養(yǎng)殖為主的長(zhǎng)江中下游武山湖浮游植物群落結(jié)構(gòu)特征*

    2020-10-29 11:13:20張佳敏張澤慧廖明軍趙以軍王和云楊雨慧
    湖泊科學(xué) 2020年6期
    關(guān)鍵詞:硅藻門(mén)武山綠藻

    張佳敏,高 健,楊 誠(chéng),張澤慧,潘 超,廖明軍,趙以軍,王和云,楊雨慧

    (湖北工業(yè)大學(xué)河湖生態(tài)修復(fù)與藻類(lèi)利用湖北省重點(diǎn)實(shí)驗(yàn)室,土木建筑與環(huán)境學(xué)院,武漢 430068)

    武山湖地處湖北省武穴市北郊,位于長(zhǎng)江北岸,距長(zhǎng)江約4 km,屬長(zhǎng)江中游中型淺水通江湖泊. 1949年,武山湖水位16.0 m時(shí),水面面積25.0 km2[1],現(xiàn)實(shí)測(cè)面積15.2 km2. 武山湖漁場(chǎng)屬?lài)?guó)營(yíng)漁場(chǎng),專(zhuān)業(yè)養(yǎng)殖已歷時(shí)60余年. 20世紀(jì)80年代初,武山湖水質(zhì)清澈、水生高等植物多樣性豐富,但自1987年以來(lái)其周?chē)?家工廠向其排污,以及湖泊本身放養(yǎng)量逐漸增大,再加上施肥和投餌養(yǎng)殖,使湖泊水生生態(tài)系統(tǒng)日趨退化,水生高等植物種類(lèi)相繼消失,全湖變?yōu)楦粻I(yíng)養(yǎng)化藻型生態(tài)系統(tǒng)[2-3],其魚(yú)類(lèi)養(yǎng)殖品種逐漸變成以單一的鰱、鳙為主. 近年來(lái),為貫徹落實(shí)國(guó)家和湖北省對(duì)長(zhǎng)江大保護(hù)的重大決策部署,針對(duì)武山湖水污染實(shí)施了一系列基礎(chǔ)建設(shè)工程,對(duì)流域?qū)嵤┝丝匚邸⒐S搬遷和禁止施肥、投餌養(yǎng)殖等措施. 但由于控污不徹底(表1),武山湖水質(zhì)長(zhǎng)期處于富營(yíng)養(yǎng)化水平,水體主要超標(biāo)因子為總磷、總氮以及高錳酸鹽指數(shù),同時(shí)水體透明度低,藻類(lèi)水華暴發(fā)的風(fēng)險(xiǎn)較大.

    近年武山湖冬捕漁獲量中鰱(Hypophthalmichthysmolitrix)和鳙(Aristichthysnobilis)占90%以上(表2),由于鰱、鳙主要以浮游生物為食,鰱、鳙養(yǎng)殖是否能夠?qū)ξ渖胶∮沃参锏倪^(guò)度增長(zhǎng)或藍(lán)藻水華的發(fā)生起到控制作用備受湖泊管理部門(mén)關(guān)注. 鰱、鳙控藻技術(shù)即“非經(jīng)典生物操縱”已經(jīng)廣為研究者和湖泊管理部門(mén)熟知[4-6],但對(duì)鰱、鳙控藻問(wèn)題在富營(yíng)養(yǎng)化湖泊中的應(yīng)用效果仍存在爭(zhēng)議[7]. 本研究基于對(duì)武山湖水體營(yíng)養(yǎng)鹽濃度持續(xù)維持在較高水平和入湖河流污染仍然較高的認(rèn)知前提下,于2017-2018年對(duì)武山湖浮游植物群落結(jié)構(gòu)及其環(huán)境因子首次開(kāi)展較系統(tǒng)的調(diào)查,分析了武山湖浮游植物物種組成及對(duì)其產(chǎn)生影響的關(guān)鍵環(huán)境因子. 研究假設(shè)以鰱、鳙養(yǎng)殖為主的武山湖浮游植物種群在濾食性魚(yú)類(lèi)牧食壓力下其優(yōu)勢(shì)種主要為小型藻類(lèi),但其密度和豐度仍將維持在較高水平. 本研究結(jié)果將為合理評(píng)估鰱、鳙放養(yǎng)在富營(yíng)養(yǎng)湖泊中的控藻作用增添原位資料.

    圖1 武山湖采樣點(diǎn)位置Fig.1 Location of sampling sites in Lake Wushan

    1 材料與方法

    1.1 研究地點(diǎn)和采樣時(shí)間

    武山湖(29°53′03″~29°56′23″N, 115°33′57″~115°36′40″E)水域面積為15.2 km2,常水位湖容約2690×104m3,最大水深為2.46 m,平均水深為1.27 m. 武山湖位于華陽(yáng)水系源頭,入湖河流主要有三部分:梅川河、大金河和鐵石河在近湖處(即三河入湖口)匯水,城區(qū)東港來(lái)水和西港來(lái)水. 武山湖出口在同一地方分叉,一是經(jīng)豐收大港連通太白湖,經(jīng)梅濟(jì)港流入黃梅縣龍感湖,再經(jīng)安徽省華陽(yáng)閘入長(zhǎng)江;二是經(jīng)官司橋大港通過(guò)閘入長(zhǎng)江(圖1).

    采樣時(shí)間為2017年9月-2018年8月,水樣及浮游植物樣品1個(gè)月采集1次. 武山湖3個(gè)主要入湖河口水樣3個(gè)月采集1次.

    2017-2018年對(duì)武山湖三大入湖河口水質(zhì)監(jiān)測(cè)的結(jié)果表明(表1),武山湖來(lái)水的氮、磷污染嚴(yán)重. 東港入湖口在2017年8月、11月和2018年2月總氮(TN)濃度均大于10.0 mg/L,總磷(TP)濃度大于1.30 mg/L. 西港入湖口在2017年8月和2018年2月TN濃度大于10.0 mg/L,TP濃度大于0.80 mg/L. 2018年2月三河入湖口的氮、磷污染最嚴(yán)重,TN和TP濃度分別為3.0和0.30 mg/L. 根據(jù)查得的3個(gè)入湖河口的年均流量可計(jì)算3個(gè)入湖口匯入武山湖的氮、磷總量,其中東港一年排入武山湖的TN總量約為54.96 t,TP約為5.79 t;西港匯入的TN總量約為36.6 t,TP約為3.09 t. 三河入湖口匯入氮、磷最多,TN約為240.22 t,TP約為32.11 t.

    表1 2017年8月-2018年5月武山湖3個(gè)主要入湖河口不同季節(jié)的水質(zhì)特征

    武山湖養(yǎng)殖投放魚(yú)苗主要是鳙和鰱(表2);近四年漁獲量約在1300~2000 t/a,其中鰱約占20%,鳙約占70%(表2),其他魚(yú)類(lèi)包括鯉(Cyprinuscarpio)、鳊(Parabramispekinensis)、鯽(Carassiusauratus)、青魚(yú)(Mylopharyngodonpiceus)、鲌(Culters)、草魚(yú)(Ctenopharyngodonidella)、黃顙魚(yú)(Pelteobagrusfulvidraco)、鳑鲏(RhodeussinensisGunther)、條(HemicculterLeuciclus)等.

    表2 武山湖2014-2018年魚(yú)類(lèi)捕撈量統(tǒng)計(jì)*

    *投放鳙魚(yú)苗約0.25~0.5 kg/尾,鰱魚(yú)苗約0.1~0.15 kg/尾. 投放量和捕撈量的單位均為萬(wàn)kg;數(shù)據(jù)來(lái)源于武山湖武黃漁場(chǎng). “-”表明無(wú)詳細(xì)統(tǒng)計(jì)數(shù)據(jù).

    1.2 樣本采集及分析

    透明度(SD)采用塞氏盤(pán)(Secchi-disk)現(xiàn)場(chǎng)測(cè)定. 水溫(WT)、pH值和溶解氧(DO)等參數(shù)采用YSI ProPlus多參儀(美國(guó))現(xiàn)場(chǎng)測(cè)定.

    用5 L有機(jī)玻璃采水器采集水表面下0.5 m處水樣,共計(jì)20 L水樣放置在用湖水洗凈的塑料桶內(nèi),混合后取約2.5 L水樣低溫保存,帶回實(shí)驗(yàn)室分析水質(zhì). TN、TP和高錳酸鹽指數(shù)(CODMn)等水質(zhì)指標(biāo)的檢驗(yàn)方法依據(jù)《水和廢水監(jiān)測(cè)分析方法(第四版)》[8].

    取上述混合水樣1 L置于廣口塑料瓶中用于浮游植物定量分析,現(xiàn)場(chǎng)加入約水樣體積1%的魯哥試劑固定,帶回實(shí)驗(yàn)室靜置沉淀24 h以上后濃縮并在尼康倒置顯微鏡下計(jì)數(shù),浮游植物鑒定參照《中國(guó)淡水藻類(lèi)——系統(tǒng),分類(lèi)及生態(tài)》[9]盡量鑒定到種,其中部分硅藻鑒定到屬. 計(jì)數(shù)方法為目鏡視野法,用0.1 mL浮游植物計(jì)數(shù)框在10×40倍光學(xué)顯微鏡下進(jìn)行,一般計(jì)數(shù)50個(gè)視野,使所得到的細(xì)胞數(shù)在300以上,每個(gè)樣品至少計(jì)數(shù)2片,對(duì)量小而個(gè)體大的種類(lèi)在10×10倍下全片計(jì)數(shù). 根據(jù)濃縮倍數(shù)換算為每升水樣中的細(xì)胞數(shù)(cells/L),即浮游植物的豐度. 由于浮游植物的比重接近1,故可以直接由浮游植物的體積換算為生物量(鮮重),即生物量為浮游植物的豐度乘以各自的平均體積[10],單位為mg/L.

    1.3 數(shù)據(jù)處理及分析

    1.3.1 綜合營(yíng)養(yǎng)狀態(tài)指數(shù)評(píng)價(jià) 參考湖泊富營(yíng)養(yǎng)化評(píng)價(jià)方法及分級(jí)標(biāo)準(zhǔn)[11],依據(jù)Chl.a、TP、TN、SD和CODMn監(jiān)測(cè)結(jié)果,運(yùn)用綜合營(yíng)養(yǎng)狀態(tài)指數(shù)法對(duì)武山湖生態(tài)系統(tǒng)健康狀況進(jìn)行評(píng)價(jià). 綜合營(yíng)養(yǎng)狀態(tài)指數(shù)計(jì)算公式如下:

    TLI(∑)=∑Wj·TLI(j)

    (1)

    式中,TLI(∑)表示綜合營(yíng)養(yǎng)狀態(tài)指數(shù);TLI(j)代表第j種參數(shù)的營(yíng)養(yǎng)狀態(tài)參數(shù);Wj為第j種參數(shù)的營(yíng)養(yǎng)狀態(tài)指數(shù)的相關(guān)權(quán)重,詳細(xì)計(jì)算公式見(jiàn)參考文獻(xiàn)[11].

    參考《中國(guó)水資源公報(bào)》[12]中湖泊、水庫(kù)富營(yíng)養(yǎng)化評(píng)分與分類(lèi)標(biāo)準(zhǔn)將湖泊營(yíng)養(yǎng)狀態(tài)分成5級(jí),同一營(yíng)養(yǎng)狀態(tài)下,指數(shù)值越高,其營(yíng)養(yǎng)程度越嚴(yán)重:貧營(yíng)養(yǎng)TLI(∑)<30,中營(yíng)養(yǎng)30≤TLI(∑)≤ 50,輕度富營(yíng)養(yǎng)50 70.

    1.3.2 優(yōu)勢(shì)種 浮游植物優(yōu)勢(shì)種根據(jù)種的Mcnaughton優(yōu)勢(shì)度指數(shù)(Y)來(lái)確定[13]:

    (2)

    式中,ni為第i種的個(gè)體數(shù);N為所有種類(lèi)總個(gè)體數(shù);fi為出現(xiàn)頻率;Y值>0.02的種類(lèi)視為優(yōu)勢(shì)種. 優(yōu)勢(shì)種的變化在一定程度上可直接反映水體污染程度和環(huán)境條件的改變[14-15].

    1.3.3 統(tǒng)計(jì)分析 浮游植物豐度和生物量與TN、TP、WT、DO、pH、SD和TLI指數(shù)的Pearson相關(guān)性分析用SPSS Statistics 22.0軟件進(jìn)行. 所有數(shù)據(jù)分析前首先進(jìn)行正態(tài)分布檢驗(yàn),不符合正態(tài)分布的數(shù)據(jù)用ln(x)處理后分析.

    冗余分析通過(guò)CANOCO 4.5進(jìn)行,在分析前,先對(duì)物種豐度矩陣進(jìn)行DCA去趨勢(shì)分析,用于排序的浮游植物物種滿足以下3個(gè)條件:(1)Y≥0.02的優(yōu)勢(shì)屬;(2)該物種在不同采樣點(diǎn)的不同時(shí)期出現(xiàn)頻率>12.5%;(3)該物種至少在1個(gè)采樣點(diǎn)的相對(duì)密度≥1%. 對(duì)篩選后的武山湖浮游植物物種進(jìn)行DCA分析,本研究中最大軸小于3,因此采用冗余分析(RDA)[16]探討武山湖浮游植物豐度與環(huán)境因子(包括WT、DO、pH、TN、TP和SD,而CODMn值受浮游植物濃度的影響而被排除在外)的關(guān)系. 為了獲得數(shù)據(jù)正態(tài)分布,將浮游植物和環(huán)境因子(除pH)均進(jìn)行l(wèi)n(x+1)轉(zhuǎn)換. 物種及環(huán)境因子選取12個(gè)月7個(gè)采樣點(diǎn)共84個(gè)數(shù)據(jù)進(jìn)行RDA分析. 在線性排序圖,從樣方的點(diǎn)在物種箭頭及其延長(zhǎng)線的投影點(diǎn)可以表示它們之間的關(guān)系. 由樣方到物種箭頭線的垂點(diǎn)可以判斷樣方物種的多度值,垂點(diǎn)在物種箭頭方向越遠(yuǎn)位置代表其物種多度越高,反之越低[16].

    聚類(lèi)分析利用大型多元統(tǒng)計(jì)軟件PRIMER 5.0進(jìn)行,對(duì)武山湖浮游植物優(yōu)勢(shì)種分別進(jìn)行時(shí)間和空間的聚類(lèi)分析.

    2 結(jié)果與分析

    2.1 武山湖水體理化特征和季節(jié)變化規(guī)律

    采樣期間,武山湖WT范圍在7.8~34℃. 水體DO濃度均值在5.57~15.03 mg/L,其中2、3月各采樣點(diǎn)DO波動(dòng)較大. 水體pH值的范圍在7.41~9.12之間,其中2、3和5月各采樣點(diǎn)pH波動(dòng)較大. 武山湖水體SD 12月最高,均值為132.7 cm;而在4月最低,均值為27.7 cm,其中有8個(gè)月的SD在40 cm以下,水體除12月外均渾濁(圖2).

    武山湖TP濃度均值范圍為0.17~1.05 mg/L,其值冬季低、夏季高,最低值在12月,最高值在8月. TN濃度均值變化范圍在1.34~3.14 mg/L之間,其值冬季高、夏季低,最低值在9月,最高值在2月,與TP季節(jié)變化相反;盡管2月TN均值最高,但各采樣點(diǎn)濃度有明顯差異,其中S6最高(4.51 mg/L),S2最低(2.09 mg/L). 武山湖CODMn12月最低,均值為4.44 mg/L,而2月最高,均值為8.88 mg/L. 武山湖不同采樣點(diǎn)Chl.a濃度差別很大,12月值最低,平均6.22 μg/L;7月Chl.a濃度最高,達(dá)125.08 μg/L(圖2). 武山湖富營(yíng)養(yǎng)指數(shù)TLI均值變化范圍在51.09~73.1之間,全湖全年處于輕度富營(yíng)養(yǎng)到重度富營(yíng)養(yǎng)水平,其中12月和1月處于輕度富營(yíng)養(yǎng),7月、8月和11月為重度富營(yíng)養(yǎng)化,其他月份為中度富營(yíng)養(yǎng)化.

    圖2 武山湖水體理化指標(biāo)的時(shí)空變化Fig.2 Spatial-temporal change of physicochemical parameters in Lake Wushan

    2.2 武山湖浮游植物特征

    2.2.1 浮游植物種類(lèi)組成及優(yōu)勢(shì)種 采樣期間武山湖共鑒定出浮游植物100種(屬),隸屬于7門(mén)62屬. 其中綠藻門(mén)種類(lèi)數(shù)最多,有28屬50種,占總種數(shù)的50%. 硅藻門(mén)其次,有11屬19種,占總數(shù)的19%. 藍(lán)藻門(mén)有10屬14種,甲藻門(mén)和隱藻門(mén)分別有4屬5種和3屬6種,裸藻門(mén)和金藻門(mén)最少,分別有4屬4種和2屬2種.

    浮游植物優(yōu)勢(shì)種共有23種,其中藍(lán)藻門(mén)有9種,綠藻門(mén)有8種,硅藻門(mén)和隱藻門(mén)各有3種. 其中綠藻門(mén)的四尾柵藻(Scenedesmusquadricauda)在12次調(diào)查中作為優(yōu)勢(shì)種出現(xiàn)了11次,硅藻門(mén)的小環(huán)藻屬一種(Cyclotellasp.)出現(xiàn)了9次,藍(lán)藻門(mén)的偽魚(yú)腥藻屬一種(Pseudanabaenasp.)與細(xì)小平裂藻(Merismopediatenuissima)均出現(xiàn)7次. 夏季和秋季藍(lán)藻門(mén)優(yōu)勢(shì)種最多且優(yōu)勢(shì)度高,秋季和春季綠藻門(mén)和硅藻門(mén)優(yōu)勢(shì)種多且優(yōu)勢(shì)度高,隱藻門(mén)的種類(lèi)僅在2018年1-3月為優(yōu)勢(shì)種. 從秋季、冬季、春季到夏季,武山湖最優(yōu)勢(shì)種的演替過(guò)程為細(xì)小平裂藻、小環(huán)藻一種(Cyclotellasp. )、小球藻一種(Chlorellasp.)再到微囊藻一種(Microcystissp.)、細(xì)小平裂藻(表3).

    以武山湖浮游植物優(yōu)勢(shì)種群的豐度進(jìn)行時(shí)間聚類(lèi)分析(圖3a),在相似水平47%時(shí),2018年1月和其他月份分為兩類(lèi);在相似水平63%時(shí),2018年1月為一類(lèi),2018年2-4月和2017年12月為一類(lèi),2017年9-11月和2018年5-8月為一類(lèi). 以武山湖浮游植物優(yōu)勢(shì)種群的豐度進(jìn)行空間聚類(lèi)分析(圖3b),結(jié)果表明7個(gè)采樣點(diǎn)相似度很高,在相似水平為93.5%時(shí)僅分為兩大類(lèi),一類(lèi)為S1、S2和S3,一類(lèi)為S4、S5、S6和S7.

    圖3 武山湖浮游植物優(yōu)勢(shì)種豐度的時(shí)間(a)和空間(b)聚類(lèi)分析Fig.3 Temporal(a) and spatial(b) cluster analysis of phytoplankton dominant species in Lake Wushan

    2.2.2 浮游植物豐度和生物量 武山湖采樣期間各采樣點(diǎn)浮游植物豐度在0.41×106~341.46×106cells/L之間,其中2017年12月全湖浮游植物平均豐度最低,為1.10×106cells/L;2018年7月浮游植物平均豐度最高,為230.10×106cells/L. 在5-11月,藍(lán)藻門(mén)的豐度占比最高,占總豐度的65.14%~88.76%. 綠藻門(mén)的豐度在2017年12月、2018年1月和4月較高,其中2018年1月浮游植物豐度突然增加,尤其是綠藻門(mén)豐度達(dá)100.11×106cells/L,占總豐度的93.01%. 2018年2月,浮游植物優(yōu)勢(shì)種逐漸從綠藻門(mén)演替成硅藻門(mén),其豐度(15.72×106cells/L)占浮游植物總豐度的75.17%. 硅藻門(mén)豐度在2017年12月和2018年3-4月也比較高,分別占總豐度的37.68%、35.03%和30.49%. 隱藻門(mén)只在2018年3月豐度比較高,為6.08×106cells/L,占總豐度的20.14%. 其他門(mén)的浮游植物豐度均很低(圖4).

    武山湖采樣期間浮游植物生物量在0.061~34.77 mg/L之間,其中12月浮游植物生物量最低,6月浮游植物生物量最高. 7月藍(lán)藻門(mén)的生物量達(dá)4.83 mg/L,占比最多(32.55%). 11月硅藻門(mén)生物量最高,達(dá)10.71 mg/L,全年有9個(gè)月硅藻門(mén)生物量占總生物量的比例最高(圖4). 綠藻門(mén)1月生物量占比最多(達(dá)35.46%),而在6月份生物量最高(3.49 mg/L). 隱藻門(mén)生物量3月最高(3.26 mg/L),裸藻門(mén)生物量最高出現(xiàn)在6月(2.73 mg/L)(圖4).

    從豐度數(shù)據(jù)來(lái)看,夏、秋季藍(lán)藻門(mén)占優(yōu),冬、春季以硅藻和綠藻為主;從生物量數(shù)據(jù)來(lái)看,春、秋、冬硅藻門(mén)占優(yōu),夏季硅藻、藍(lán)藻和綠藻均較高(圖4).

    2.3 武山湖浮游植物與環(huán)境因子的關(guān)系

    Pearson相關(guān)性分析結(jié)果(表4)表明,TP濃度分別與浮游植物豐度和生物量、藍(lán)藻門(mén)豐度、綠藻門(mén)生物量呈顯著正相關(guān),與藍(lán)藻門(mén)生物量呈極顯著正相關(guān). 水溫與藍(lán)藻門(mén)的生物量呈極顯著正相關(guān),與浮游植物總豐度、總生物量以及綠藻門(mén)生物量呈顯著正相關(guān). pH與浮游植物總豐度、總生物量及綠藻門(mén)和硅藻門(mén)的豐度、生物量均呈極顯著正相關(guān). 透明度與浮游植物總豐度和總生物量、硅藻門(mén)的豐度和生物量呈極顯著負(fù)相關(guān),與綠藻門(mén)生物量顯著負(fù)相關(guān). 除綠藻門(mén)豐度外,浮游植物及其主要類(lèi)群豐度和生物量均與綜合營(yíng)養(yǎng)狀態(tài)指數(shù)TLI呈顯著或極顯著正相關(guān). 藍(lán)藻豐度與DO濃度呈顯著負(fù)相關(guān),浮游植物及其類(lèi)群豐度和生物量與TN濃度均不相關(guān).

    圖4 武山湖浮游植物豐度和生物量變化Fig.4 Variations of phytoplankton abundance and biomass in Lake Wushan

    表4 武山湖浮游植物豐度、生物量與環(huán)境因子相關(guān)性分析(n=84)

    通過(guò)冗余分析(RDA)對(duì)浮游植物各優(yōu)勢(shì)屬的豐度和生物量與環(huán)境因子之間的關(guān)系進(jìn)行分析. 基于豐度數(shù)據(jù),主成分軸1和軸2特征值分別是0.361和0.095,即前兩個(gè)軸共解釋了45.6%的浮游植物各優(yōu)勢(shì)屬的豐度變化(圖5a). 主成分軸1和軸2的物種與環(huán)境因子相關(guān)性分別為86.6%和74.9%,兩軸累計(jì)百分比變化率為88.4%,表明前兩個(gè)軸的物種與環(huán)境因子相關(guān)性高. RDA蒙特卡羅前向選擇檢驗(yàn)表明,TP濃度和WT對(duì)浮游植物豐度變化影響最大,這兩個(gè)變量共解釋了36.0%(圖5a). TP和WT與藍(lán)藻門(mén)各優(yōu)勢(shì)屬(Z1~Z6)豐度、綠藻門(mén)Z9、Z10和硅藻門(mén)Z13、Z14呈正相關(guān);而TN濃度與藍(lán)藻門(mén)所有優(yōu)勢(shì)屬豐度均呈負(fù)相關(guān). 透明度與各優(yōu)勢(shì)屬豐度均呈負(fù)相關(guān). TN、TP、pH、WT、DO和SD與軸1相關(guān)性系數(shù)分別為-0.33、0.81、0.38、0.76、-0.39 和-0.53,而與軸2的相關(guān)性系數(shù)分別為0.08、-0.07、-0.53、-0.04、-0.49和0.39.

    基于生物量數(shù)據(jù),主成分軸1和軸2特征值分別是0.346和0.114,即前兩個(gè)軸共解釋了46.0%的浮游植物各優(yōu)勢(shì)屬的生物量變化(圖5b). 主成分軸1和軸2的物種與環(huán)境因子相關(guān)性分別為91.4%和77.7%,兩軸累計(jì)百分比變化率為88.4%,表明前兩個(gè)軸的物種與環(huán)境因子相關(guān)性高. RDA蒙特卡羅檢驗(yàn)前向選擇檢驗(yàn)表明,WT和TP濃度對(duì)浮游植物生物量變化影響最大,這兩個(gè)變量共解釋了36.0%(圖5b). WT、TP和pH與藍(lán)藻門(mén)各優(yōu)勢(shì)屬(Z1~Z6)生物量呈正相關(guān),與綠藻門(mén)Z9、Z10、Z11和硅藻門(mén)Z13、Z14呈正相關(guān);而TN濃度與藍(lán)藻門(mén)所有優(yōu)勢(shì)屬生物量均呈負(fù)相關(guān). 透明度與各優(yōu)勢(shì)屬均呈負(fù)相關(guān). 浮游植物各優(yōu)勢(shì)屬豐度和生物量與環(huán)境因子的相互關(guān)系分析結(jié)果相似. TN、TP、pH、WT、DO和SD與軸1相關(guān)性系數(shù)分別為0.33、-0.77、-0.62、-0.87、0.12和0.69,而與軸2的相關(guān)性系數(shù)分別為-0.02、0.16、-0.47、0.12、-0.70和0.27.

    從RDA三元圖(圖5)可以看出浮游植物樣方四季分布明顯,夏季和冬季樣方分布軸2兩側(cè),其中夏季樣方更集中,而冬季樣方分布在一、四象限(豐度)或二、三象限. 秋季和春季的樣方多分布在軸2附近. 藍(lán)藻門(mén)夏季樣方的多度值靠近物種箭頭的延長(zhǎng)線上,表明其多度值在夏季高,而綠藻門(mén)、硅藻門(mén)和隱藻門(mén)在10月和11月以及4月的樣方多度值較其他月份高.

    圖5 武山湖浮游植物、環(huán)境變量和采樣點(diǎn)的冗余分析(RDA)三元圖Fig.5 Redundancy analysis (RDA) triplot (phytoplankton, environmental variables and sampling sites) in Lake Wushan

    3 討論

    藻類(lèi)大量增長(zhǎng)特別是藍(lán)藻水華的暴發(fā)以及死亡,會(huì)嚴(yán)重影響湖泊水質(zhì)和生態(tài)系統(tǒng)的健康,因此成為湖泊治理的一個(gè)重要問(wèn)題. 根據(jù)相關(guān)湖泊富營(yíng)養(yǎng)評(píng)價(jià)標(biāo)準(zhǔn),盡管武山湖各采樣點(diǎn)由于距離入湖河口不同導(dǎo)致?tīng)I(yíng)養(yǎng)鹽指標(biāo)存在差異(圖2),但從整體上來(lái)講,武山湖TN(1.34~3.14 mg/L)、TP(0.17~1.05 mg/L)、Chl.a濃度(峰值達(dá)125.08 μg/L)以及藻類(lèi)豐度和生物量均可指示武山湖水體處于富營(yíng)養(yǎng)化水平[17],具有藻類(lèi)大量生長(zhǎng)的營(yíng)養(yǎng)條件,再加上武山湖藻類(lèi)優(yōu)勢(shì)種主要是藍(lán)藻和綠藻,藻類(lèi)水華尤其是藍(lán)藻水華發(fā)生的風(fēng)險(xiǎn)高. 本次調(diào)查中12次采樣共鑒定出23種優(yōu)勢(shì)種,其中藍(lán)藻門(mén)占39.1%、綠藻門(mén)占34.9%,夏季和秋季藍(lán)藻門(mén)優(yōu)勢(shì)種最多且優(yōu)勢(shì)度高,冬季和春季綠藻門(mén)和硅藻門(mén)優(yōu)勢(shì)種多且優(yōu)勢(shì)度高. 在長(zhǎng)江中下游類(lèi)似富營(yíng)養(yǎng)水平湖泊中,浮游植物通常都是以綠藻和藍(lán)藻為主. 如長(zhǎng)江中游的長(zhǎng)湖[18]、沙湖[19]、太平湖[20]、巢湖[21]、青山湖[22]以及下游的太湖梅梁灣[23]、滆湖[24]、淀山湖[25]等湖泊,常見(jiàn)的優(yōu)勢(shì)種有微囊藻、平裂藻、小球藻、柵藻等屬的種類(lèi),其中藍(lán)藻主要以群體性如微囊藻、魚(yú)腥藻等占優(yōu)勢(shì),綠藻以柵藻等為主. 與上述這些湖泊相比,武山湖的氮、磷水平接近甚至略高,其藻優(yōu)勢(shì)種及其季節(jié)變化與這些已研究報(bào)道的湖泊呈現(xiàn)類(lèi)似特征. 然而,從豐度和生物量上來(lái)看,武山湖浮游植物豐度在0.41×106~341.46×106cells/L之間,峰值已經(jīng)遠(yuǎn)高于太平湖(6.41×106cells/L)[20]、淀山湖(102×106cells/L)[25]浮游植物豐度峰值,但小于青山湖(726×106cells/L)[22]和巢湖浮游植物豐度峰值(約800×106cells/L)[21]. 根據(jù)報(bào)道,淀山湖、青山湖和巢湖均曾暴發(fā)成片藻類(lèi)水華,而武山湖目前Chl.a濃度、豐度等均達(dá)到很高水平,但是在監(jiān)測(cè)過(guò)程中,沒(méi)有成片藻類(lèi)聚集現(xiàn)象發(fā)生.

    RDA分析表明武山湖的WT、TP和pH是影響浮游植物增長(zhǎng)及其種群分布的關(guān)鍵理化因子. 溫度是影響浮游植物的最重要理化因子之一,溫度不僅會(huì)影響浮游植物生長(zhǎng)率,而且會(huì)影響其群落結(jié)構(gòu)變化[26-28]. 在武山湖,水溫是影響浮游植物季節(jié)變化的最重要環(huán)境因子,與藍(lán)藻、綠藻和硅藻均呈正相關(guān),并可獨(dú)立解釋28%的浮游植物豐度和32%的生物量變化. 武山湖夏季(6-8月)浮游植物以藍(lán)藻門(mén)微囊藻、細(xì)小平裂藻等藍(lán)藻為主,高溫和高磷濃度可能是導(dǎo)致其大量生長(zhǎng)的主要原因,大量研究也表明夏季藻類(lèi)種群分布的主要影響因子是溫度和磷[21,23]. 聚類(lèi)分析(圖3a)和RDA分析(圖5)可清晰表明,武山湖浮游植物群落結(jié)構(gòu)的變化具有明顯的季節(jié)特征,當(dāng)水溫升高時(shí)綠藻和藍(lán)藻的豐度和生物量均顯著增高,原因?yàn)樗{(lán)藻、綠藻喜高溫,其主要出現(xiàn)在溫度較高的季節(jié)[29-32]. 一般來(lái)說(shuō),藍(lán)藻、綠藻的生長(zhǎng)最適溫度為25~35℃[33],15℃以下其生長(zhǎng)受到限制[34]. 本研究中藍(lán)藻的生長(zhǎng)與上述研究結(jié)果一致,但綠藻部分優(yōu)勢(shì)種的生長(zhǎng)與其不一致,如綠藻門(mén)小球藻一種曾在1月份成為優(yōu)勢(shì)種,可能受到冬季捕撈等因素的影響. 小球藻在溫度較低的季節(jié)易成為優(yōu)勢(shì)種[19,23],武山湖冬季捕撈后浮游植物面臨的牧食壓力大大降低,而小球藻生長(zhǎng)繁殖速度又快,這可能是其在1月份成為最優(yōu)勢(shì)種的原因. 武山湖冬季和春季隱藻、硅藻占優(yōu)勢(shì)與多數(shù)研究結(jié)果一致[34].

    而對(duì)于氮磷來(lái)說(shuō),其往往是浮游植物生長(zhǎng)和分布的限制因子[35]. 本研究中磷濃度的增加有利于藍(lán)藻、綠藻和硅藻的生長(zhǎng),其對(duì)藍(lán)藻增長(zhǎng)的影響尤為明顯. 磷是浮游植物存活的必要條件,浮游植物的物種多樣性、豐度和生物量與營(yíng)養(yǎng)鹽的組成和濃度是密不可分的,浮游植物的增長(zhǎng)一般與磷呈正相關(guān)[35]. 但對(duì)于氮來(lái)說(shuō),RDA分析表明武山湖TN與浮游植物生長(zhǎng)呈負(fù)相關(guān)且相關(guān)性不顯著,出現(xiàn)負(fù)相關(guān)的原因可能是因?yàn)槲渖胶N濃度因外源污染排入所出現(xiàn)的峰值與藻類(lèi)的生長(zhǎng)峰值季節(jié)上不同步,同時(shí)溫度升高在促進(jìn)藻類(lèi)生長(zhǎng)的同時(shí)會(huì)造成水體脫氮能力增加[36]. 在一些研究中也發(fā)現(xiàn)TN濃度與浮游植物生長(zhǎng)呈負(fù)相關(guān)[15,37-38],水體中氮濃度峰值一般出現(xiàn)在初春,由于溫度低限制了藻類(lèi)生長(zhǎng),使得氮濃度峰值和藻類(lèi)生物量峰值產(chǎn)生了季節(jié)錯(cuò)位. 例如,在太湖水體中氮、磷濃度就存在類(lèi)似武山湖的季節(jié)變化,氮在春季和冬季濃度較高,而在夏季和秋季濃度較低;磷濃度的季節(jié)變化與氮相反,春季和冬季濃度較低,而夏季和秋季相對(duì)較高[39]. pH對(duì)浮游植物群落也有著重要影響,當(dāng)pH>8.0時(shí),藍(lán)藻往往成為主要優(yōu)勢(shì)種[40]. 同時(shí),藍(lán)藻可通過(guò)一系列的生理生化反應(yīng)來(lái)調(diào)節(jié)水體pH,使其維持在適宜藻類(lèi)生長(zhǎng)的范圍內(nèi)[41]. 藻類(lèi)種群分布與水體pH值的相關(guān)性也可能是由于不同環(huán)境條件下藻類(lèi)種群的生長(zhǎng)繁殖變化導(dǎo)致水體pH值變化所致[42].

    此外,在武山湖,浮游植物群落結(jié)構(gòu)及現(xiàn)存量對(duì)以鰱、鳙尤其是鳙為主的魚(yú)類(lèi)結(jié)構(gòu)的響應(yīng)值得討論. 采樣期間武山湖放養(yǎng)的魚(yú)類(lèi)主要是鰱(約21%)、鳙(約70%),其中鰱現(xiàn)存量約為12.6 g/m3、鳙約為33.8 g/m3,二者比值約3:7. 一般認(rèn)為鰱主要攝食浮游植物,而鳙主要攝食浮游動(dòng)物[43];但劉恩生等研究發(fā)現(xiàn)在以水華藍(lán)藻為主要初級(jí)生產(chǎn)力的富營(yíng)養(yǎng)化湖泊中,鰱、鳙的食性會(huì)根據(jù)環(huán)境的變化而發(fā)生改變,它們均以微囊藻為主要食物[44];王媛等基于穩(wěn)定同位素技術(shù)也證實(shí)鰱、鳙的主要食物來(lái)源為懸浮顆粒有機(jī)物,而懸浮顆粒有機(jī)物的主要來(lái)源為水華藍(lán)藻[45]. 本研究同期采集的魚(yú)類(lèi)同位素測(cè)定結(jié)果也顯示鰱和鳙的δ15N值均接近13.0‰[46],無(wú)顯著區(qū)別,表明武山湖鰱、鳙魚(yú)食性可能類(lèi)似. 在調(diào)查中發(fā)現(xiàn),武山湖2018年夏季浮游植物最優(yōu)勢(shì)種為微囊藻,豐度也很高,但是從2017-2018年周年監(jiān)測(cè)期間并沒(méi)有出現(xiàn)大規(guī)模的水華及藻類(lèi)聚集在水面的現(xiàn)象. 李林春對(duì)南灣水庫(kù)投放鰱、鳙研究表明,當(dāng)鰱與鳙的比值在4:6的情況下,能顯著降低浮游植物的數(shù)量和生物量,同時(shí)也大大降低了藍(lán)藻在藻類(lèi)群落中所占的比例及其生物量,有效地避免了南灣水庫(kù)藍(lán)藻水華的暴發(fā)[5]. 在1998和1999年,千島湖連續(xù)暴發(fā)了2次大面積的藍(lán)藻水華,據(jù)研究,鰱、鳙生物量的減少可能是其主要原因[47]. 馮超群等于2016年7月,在太湖三國(guó)城水域設(shè)置孔徑為1.5 cm的大型尼龍網(wǎng)圍欄,面積約866.7 hm2,鰱、鳙比例為1:1,初始鰱、鳙約2.9 g/尾,生物量約3.27 g/m3,捕撈時(shí)圍網(wǎng)內(nèi)鰱、鳙生物量約82.85 g/m3,結(jié)果表明單純放養(yǎng)鰱、鳙適合控制藍(lán)藻水華而非藻類(lèi)總量,試驗(yàn)期間圍欄內(nèi)大型藻類(lèi)如束絲藻、假魚(yú)腥藻屬等生物量顯著下降,到試驗(yàn)結(jié)束時(shí)幾乎消失[48]. 劉其根等綜合分析眾多鰱、鳙控藻研究結(jié)果,認(rèn)為無(wú)論是基于穩(wěn)定性同位素或PCR-DGGE食性分析,還是通過(guò)合理的圍隔實(shí)驗(yàn),以及大量的養(yǎng)殖實(shí)踐,均證實(shí)了鰱、鳙是能夠控制藍(lán)藻過(guò)度生長(zhǎng),但其應(yīng)用效果受多種因素影響[7]. 因此,盡管對(duì)鰱、鳙放養(yǎng)是否能把浮游植物的數(shù)量和生物量控制在較低水平目前還存在很多爭(zhēng)議,但是對(duì)當(dāng)鰱、鳙養(yǎng)殖密度達(dá)到一定程度時(shí)可有效控制富營(yíng)養(yǎng)水體中藍(lán)藻的過(guò)度生長(zhǎng),且鰱、鳙的濾食作用會(huì)阻止藻類(lèi)在水面聚集成片這一觀點(diǎn)得到較多研究驗(yàn)證.

    4 結(jié)論

    1)綜合理化因子、營(yíng)養(yǎng)狀態(tài)指數(shù)及浮游植物豐度、生物量等指標(biāo),采樣期間武山湖水質(zhì)營(yíng)養(yǎng)程度處于中度-重度富營(yíng)養(yǎng)水平.

    2)武山湖夏季和秋季藍(lán)藻門(mén)優(yōu)勢(shì)種最多且優(yōu)勢(shì)度高,冬季和春季綠藻門(mén)和硅藻門(mén)優(yōu)勢(shì)種多且優(yōu)勢(shì)度高,每月優(yōu)勢(shì)度最大的種類(lèi)主要有藍(lán)藻門(mén)的微囊藻一種和細(xì)小平裂藻、綠藻門(mén)的小球藻一種與硅藻門(mén)的小環(huán)藻一種,表明武山湖浮游植物優(yōu)勢(shì)種主要是小型藻類(lèi).

    3)武山湖浮游植物豐度和生物量在夏季均較高,夏季發(fā)生大規(guī)模藻類(lèi)水華的風(fēng)險(xiǎn)較高. RDA和線性相關(guān)分析結(jié)果表明影響浮游植物豐度和生物量的主要營(yíng)養(yǎng)元素是磷,磷濃度增加與藍(lán)藻和綠藻的豐度和生物量呈正相關(guān),表明外源磷輸入和內(nèi)源磷釋放的控制應(yīng)引起重視. 在武山湖目前外源污染量還較大的情況下,鰱、鳙的放養(yǎng)對(duì)浮游植物過(guò)度生長(zhǎng)的控制具有積極意義.

    猜你喜歡
    硅藻門(mén)武山綠藻
    又吟武山湖
    石佛寺水庫(kù)浮游植物時(shí)間分布相關(guān)研究
    地下水(2022年2期)2022-05-19 12:49:52
    大連松樹(shù)水庫(kù)浮游植物特征及對(duì)環(huán)境因子的響應(yīng)
    香榧綠藻的生物學(xué)特性及物種鑒定
    無(wú)綠藻的微生物學(xué)特性和菌種鑒定
    我的家鄉(xiāng)在武山
    全球近10年無(wú)綠藻病例報(bào)道文獻(xiàn)的回顧
    鐵甲水庫(kù)浮游植物種類(lèi)組成及其易引發(fā)水華浮游植物種類(lèi)探討
    以成長(zhǎng)為主題解讀《窗燈》
    高位循環(huán)水池塘與普通池塘高溫時(shí)節(jié)浮游植物群落的比較
    日韩视频一区二区在线观看| 久久久久久国产a免费观看| 天天躁夜夜躁狠狠躁躁| 97人妻天天添夜夜摸| 久久伊人香网站| 最近最新中文字幕大全电影3 | 国产欧美日韩一区二区精品| 视频在线观看一区二区三区| 亚洲av日韩精品久久久久久密| 精品国产美女av久久久久小说| 久久久久久免费高清国产稀缺| 热re99久久国产66热| 国产精品永久免费网站| 中出人妻视频一区二区| 久久性视频一级片| 操美女的视频在线观看| 免费高清视频大片| 国产一区二区三区综合在线观看| 日韩精品青青久久久久久| 美女免费视频网站| 亚洲自拍偷在线| 亚洲中文日韩欧美视频| 91大片在线观看| 两个人看的免费小视频| 少妇的丰满在线观看| 亚洲avbb在线观看| 婷婷丁香在线五月| 久久久久久亚洲精品国产蜜桃av| 久久伊人香网站| 精品免费久久久久久久清纯| 一个人观看的视频www高清免费观看 | 欧美老熟妇乱子伦牲交| 国产精品野战在线观看| www日本在线高清视频| 欧美日韩精品网址| 婷婷丁香在线五月| 成人18禁在线播放| 一卡2卡三卡四卡精品乱码亚洲| 欧美成人免费av一区二区三区| 18美女黄网站色大片免费观看| 亚洲五月婷婷丁香| 桃色一区二区三区在线观看| 日韩欧美免费精品| 69av精品久久久久久| 宅男免费午夜| 天天添夜夜摸| 午夜福利免费观看在线| 国产高清激情床上av| 在线免费观看的www视频| 久久精品国产清高在天天线| 亚洲第一av免费看| 在线观看一区二区三区| 乱人伦中国视频| 1024香蕉在线观看| 欧美日本亚洲视频在线播放| 免费看a级黄色片| 国产视频一区二区在线看| www.自偷自拍.com| 精品国产超薄肉色丝袜足j| 女人被躁到高潮嗷嗷叫费观| 欧美中文综合在线视频| 亚洲午夜理论影院| 成人精品一区二区免费| 国产精品久久电影中文字幕| 欧美久久黑人一区二区| 日本a在线网址| 久久香蕉激情| 国产亚洲欧美98| 俄罗斯特黄特色一大片| 中文字幕另类日韩欧美亚洲嫩草| 亚洲成av片中文字幕在线观看| 女生性感内裤真人,穿戴方法视频| 99re在线观看精品视频| 中文字幕人成人乱码亚洲影| 亚洲精品在线美女| 久久国产精品男人的天堂亚洲| 久久精品国产亚洲av香蕉五月| 国产三级在线视频| 国产不卡一卡二| 啦啦啦免费观看视频1| 亚洲中文av在线| av视频在线观看入口| 国产av一区二区精品久久| 亚洲精品一区av在线观看| 国产精品久久久久久精品电影 | 国产免费av片在线观看野外av| 亚洲男人天堂网一区| 日韩欧美一区二区三区在线观看| 嫩草影院精品99| 国产精品av久久久久免费| 亚洲人成电影观看| 在线观看午夜福利视频| 一进一出抽搐动态| 一边摸一边抽搐一进一小说| 国产成人影院久久av| bbb黄色大片| 亚洲欧美激情综合另类| 成年女人毛片免费观看观看9| 九色国产91popny在线| 久久青草综合色| 又黄又粗又硬又大视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲中文字幕一区二区三区有码在线看 | 国产av一区在线观看免费| 男人的好看免费观看在线视频 | 国产精品1区2区在线观看.| 日韩精品中文字幕看吧| 久久久久久人人人人人| 亚洲三区欧美一区| 免费看a级黄色片| 欧美老熟妇乱子伦牲交| 黄色 视频免费看| 男女下面进入的视频免费午夜 | 国产精品亚洲一级av第二区| 看黄色毛片网站| 妹子高潮喷水视频| 老司机深夜福利视频在线观看| 真人做人爱边吃奶动态| 国产视频一区二区在线看| 国产一级毛片七仙女欲春2 | 国产又爽黄色视频| 亚洲熟妇中文字幕五十中出| 黄色视频,在线免费观看| 久久精品国产清高在天天线| 亚洲午夜精品一区,二区,三区| 国产午夜福利久久久久久| 午夜激情av网站| av天堂在线播放| 日韩欧美一区视频在线观看| 国产又色又爽无遮挡免费看| 97碰自拍视频| 天堂√8在线中文| 亚洲精品粉嫩美女一区| 亚洲性夜色夜夜综合| 给我免费播放毛片高清在线观看| 国内精品久久久久久久电影| 亚洲一卡2卡3卡4卡5卡精品中文| 桃色一区二区三区在线观看| 国产亚洲欧美精品永久| 丝袜美足系列| 精品国产乱子伦一区二区三区| 国产视频一区二区在线看| 国产精品电影一区二区三区| 国产免费男女视频| 伦理电影免费视频| 9热在线视频观看99| 在线观看免费视频日本深夜| 日韩精品青青久久久久久| 国产三级黄色录像| 如日韩欧美国产精品一区二区三区| 最近最新中文字幕大全免费视频| 女人精品久久久久毛片| 日韩av在线大香蕉| 制服诱惑二区| 90打野战视频偷拍视频| 在线观看日韩欧美| 老司机午夜福利在线观看视频| 午夜日韩欧美国产| 老司机在亚洲福利影院| 中文字幕av电影在线播放| 在线观看免费视频日本深夜| 一a级毛片在线观看| 欧美黑人精品巨大| 极品教师在线免费播放| 两个人视频免费观看高清| 男人操女人黄网站| 99精品久久久久人妻精品| 精品欧美一区二区三区在线| 精品久久久久久久人妻蜜臀av | 午夜a级毛片| 黑丝袜美女国产一区| x7x7x7水蜜桃| 免费不卡黄色视频| 最近最新中文字幕大全电影3 | 亚洲欧美日韩高清在线视频| 两个人看的免费小视频| 国产午夜福利久久久久久| 国产av一区二区精品久久| 色综合站精品国产| 久99久视频精品免费| 久久草成人影院| 国产精品亚洲美女久久久| 精品久久久久久久久久免费视频| 手机成人av网站| 亚洲成人免费电影在线观看| 国产在线观看jvid| 狠狠狠狠99中文字幕| 亚洲五月天丁香| 欧美色欧美亚洲另类二区 | 动漫黄色视频在线观看| 91成人精品电影| 丝袜美足系列| 女人被躁到高潮嗷嗷叫费观| 久久国产精品影院| 日韩 欧美 亚洲 中文字幕| 欧美日本中文国产一区发布| av网站免费在线观看视频| 无人区码免费观看不卡| 亚洲av电影不卡..在线观看| 欧美日韩瑟瑟在线播放| 岛国在线观看网站| 美女大奶头视频| 色综合站精品国产| 免费观看人在逋| 国产欧美日韩一区二区精品| 精品久久蜜臀av无| 日日爽夜夜爽网站| 老司机在亚洲福利影院| 自线自在国产av| 国产一区二区三区在线臀色熟女| 日韩大码丰满熟妇| 首页视频小说图片口味搜索| 午夜福利影视在线免费观看| 黑人操中国人逼视频| 亚洲中文日韩欧美视频| 国产精品久久久久久亚洲av鲁大| 岛国视频午夜一区免费看| 国产av精品麻豆| 欧美成狂野欧美在线观看| 亚洲人成电影观看| 国产日韩一区二区三区精品不卡| 不卡av一区二区三区| 乱人伦中国视频| 99热只有精品国产| 两性夫妻黄色片| av片东京热男人的天堂| 国产日韩一区二区三区精品不卡| 宅男免费午夜| 成人手机av| 多毛熟女@视频| 欧美国产精品va在线观看不卡| 日韩中文字幕欧美一区二区| 亚洲人成网站在线播放欧美日韩| 国产精品久久久人人做人人爽| 91成年电影在线观看| 午夜免费激情av| 国内精品久久久久精免费| 激情在线观看视频在线高清| 国产一区在线观看成人免费| 久久中文字幕人妻熟女| 中文字幕另类日韩欧美亚洲嫩草| 在线观看日韩欧美| 色综合婷婷激情| 国产精品免费视频内射| 午夜福利,免费看| 亚洲七黄色美女视频| 精品久久蜜臀av无| 国产伦一二天堂av在线观看| 伦理电影免费视频| 一本大道久久a久久精品| 悠悠久久av| 女生性感内裤真人,穿戴方法视频| 久久中文字幕一级| 欧美黄色淫秽网站| 欧美不卡视频在线免费观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 首页视频小说图片口味搜索| 老熟妇乱子伦视频在线观看| 99久久99久久久精品蜜桃| 久久精品91蜜桃| 午夜福利,免费看| 在线观看免费午夜福利视频| 久久久久精品国产欧美久久久| 国产成人免费无遮挡视频| 国产三级在线视频| av视频在线观看入口| 男人的好看免费观看在线视频 | 精品国产超薄肉色丝袜足j| 国产不卡一卡二| 午夜福利欧美成人| 午夜福利在线观看吧| av片东京热男人的天堂| 香蕉久久夜色| 夜夜躁狠狠躁天天躁| 伊人久久大香线蕉亚洲五| 少妇被粗大的猛进出69影院| 国产高清有码在线观看视频 | 变态另类成人亚洲欧美熟女 | 欧美日韩亚洲国产一区二区在线观看| 看黄色毛片网站| 大型黄色视频在线免费观看| 久久人妻熟女aⅴ| 久久中文字幕人妻熟女| 午夜成年电影在线免费观看| 色综合婷婷激情| 中文亚洲av片在线观看爽| 国产亚洲精品第一综合不卡| 日韩三级视频一区二区三区| 精品一区二区三区四区五区乱码| 国产成人av激情在线播放| 岛国在线观看网站| 日本 av在线| or卡值多少钱| 欧美中文日本在线观看视频| 欧美性长视频在线观看| 好男人电影高清在线观看| 免费女性裸体啪啪无遮挡网站| 日本 欧美在线| 日日摸夜夜添夜夜添小说| 国产男靠女视频免费网站| 在线视频色国产色| 这个男人来自地球电影免费观看| 精品国内亚洲2022精品成人| 色播在线永久视频| 日韩 欧美 亚洲 中文字幕| 制服诱惑二区| 三级毛片av免费| 国语自产精品视频在线第100页| 日韩欧美一区二区三区在线观看| 精品国产乱码久久久久久男人| 91大片在线观看| 精品高清国产在线一区| 熟妇人妻久久中文字幕3abv| 日韩欧美一区二区三区在线观看| 国产精品美女特级片免费视频播放器 | 国产精品精品国产色婷婷| 久久午夜亚洲精品久久| 一级毛片精品| 午夜影院日韩av| 99精品欧美一区二区三区四区| 丝袜人妻中文字幕| 欧美乱妇无乱码| 熟妇人妻久久中文字幕3abv| 日本免费一区二区三区高清不卡 | 精品久久久精品久久久| 老汉色∧v一级毛片| 欧美日韩乱码在线| 国产极品粉嫩免费观看在线| 91在线观看av| 国产亚洲av高清不卡| 午夜福利免费观看在线| 午夜精品久久久久久毛片777| 在线十欧美十亚洲十日本专区| www.999成人在线观看| 久久久精品欧美日韩精品| 多毛熟女@视频| 美国免费a级毛片| 久久狼人影院| 老鸭窝网址在线观看| 久久 成人 亚洲| 一级a爱视频在线免费观看| 国产一区在线观看成人免费| 久久影院123| 国产视频一区二区在线看| 一级a爱视频在线免费观看| 天堂√8在线中文| 欧美国产日韩亚洲一区| 每晚都被弄得嗷嗷叫到高潮| 两个人视频免费观看高清| 91老司机精品| 精品一区二区三区视频在线观看免费| 深夜精品福利| 久久中文字幕人妻熟女| 国产日韩一区二区三区精品不卡| 成人特级黄色片久久久久久久| 女生性感内裤真人,穿戴方法视频| 欧美一级毛片孕妇| 少妇熟女aⅴ在线视频| 亚洲国产精品999在线| 免费女性裸体啪啪无遮挡网站| 啦啦啦观看免费观看视频高清 | 精品国产一区二区三区四区第35| 少妇裸体淫交视频免费看高清 | 亚洲一区中文字幕在线| 天堂动漫精品| 婷婷丁香在线五月| 熟妇人妻久久中文字幕3abv| 日韩三级视频一区二区三区| 久久久国产成人精品二区| 亚洲中文字幕一区二区三区有码在线看 | 又黄又爽又免费观看的视频| 午夜福利,免费看| 免费人成视频x8x8入口观看| 真人做人爱边吃奶动态| 国产亚洲精品第一综合不卡| 欧美成狂野欧美在线观看| 很黄的视频免费| 国产一级毛片七仙女欲春2 | 久久香蕉国产精品| svipshipincom国产片| 国产亚洲精品久久久久久毛片| 亚洲无线在线观看| 黄片播放在线免费| 国产精品亚洲美女久久久| 亚洲成av片中文字幕在线观看| 久9热在线精品视频| 午夜久久久久精精品| 男女午夜视频在线观看| 国产欧美日韩一区二区精品| 啦啦啦韩国在线观看视频| 一区二区日韩欧美中文字幕| 成在线人永久免费视频| 高潮久久久久久久久久久不卡| 国产一区二区三区视频了| 夜夜爽天天搞| 精品第一国产精品| 无遮挡黄片免费观看| 亚洲,欧美精品.| 色av中文字幕| 成人国产综合亚洲| 欧美日韩乱码在线| 欧美久久黑人一区二区| 最好的美女福利视频网| av超薄肉色丝袜交足视频| 精品国产一区二区久久| 精品不卡国产一区二区三区| 午夜福利一区二区在线看| 国产麻豆成人av免费视频| 不卡av一区二区三区| 午夜免费成人在线视频| 久久性视频一级片| 两性午夜刺激爽爽歪歪视频在线观看 | 九色亚洲精品在线播放| 久久人妻福利社区极品人妻图片| 高清在线国产一区| 12—13女人毛片做爰片一| 性色av乱码一区二区三区2| 精品国产一区二区久久| 欧美激情久久久久久爽电影 | 国产成人一区二区三区免费视频网站| 中文字幕精品免费在线观看视频| 啦啦啦 在线观看视频| 国产极品粉嫩免费观看在线| 中文字幕精品免费在线观看视频| 不卡av一区二区三区| 亚洲性夜色夜夜综合| 丰满人妻熟妇乱又伦精品不卡| 神马国产精品三级电影在线观看 | 亚洲伊人色综图| 欧美乱妇无乱码| 亚洲成人免费电影在线观看| 男女床上黄色一级片免费看| 国产成人啪精品午夜网站| 国产成人精品久久二区二区免费| 亚洲精品美女久久久久99蜜臀| 国产成人av教育| 欧美日韩亚洲国产一区二区在线观看| 欧美最黄视频在线播放免费| 99久久99久久久精品蜜桃| 国产欧美日韩一区二区三| www.精华液| 777久久人妻少妇嫩草av网站| 日日干狠狠操夜夜爽| 黄色a级毛片大全视频| 欧美午夜高清在线| 一本久久中文字幕| 亚洲色图av天堂| 亚洲av成人不卡在线观看播放网| 精品欧美国产一区二区三| 一级黄色大片毛片| 精品人妻在线不人妻| 欧美日韩福利视频一区二区| 少妇熟女aⅴ在线视频| 久久国产精品男人的天堂亚洲| 亚洲欧美激情综合另类| 大型黄色视频在线免费观看| 日韩av在线大香蕉| 欧美激情久久久久久爽电影 | 久久这里只有精品19| 麻豆成人av在线观看| 麻豆av在线久日| av天堂久久9| 黄色女人牲交| 极品人妻少妇av视频| 此物有八面人人有两片| 色综合亚洲欧美另类图片| 自拍欧美九色日韩亚洲蝌蚪91| 日韩欧美国产一区二区入口| 99国产精品一区二区蜜桃av| 欧美一级a爱片免费观看看 | 久久久久久久久免费视频了| 亚洲成人国产一区在线观看| 亚洲色图 男人天堂 中文字幕| 精品少妇一区二区三区视频日本电影| 色播亚洲综合网| av天堂久久9| 国产精品,欧美在线| 美女扒开内裤让男人捅视频| 人人妻人人澡欧美一区二区 | 国产欧美日韩一区二区三区在线| 男人的好看免费观看在线视频 | av福利片在线| 成人特级黄色片久久久久久久| 亚洲五月天丁香| 一级毛片精品| 亚洲av日韩精品久久久久久密| 国产国语露脸激情在线看| 纯流量卡能插随身wifi吗| 精品免费久久久久久久清纯| 久久人人爽av亚洲精品天堂| 国产高清激情床上av| www.精华液| 国产蜜桃级精品一区二区三区| 精品久久蜜臀av无| 纯流量卡能插随身wifi吗| 99久久国产精品久久久| 久久人人97超碰香蕉20202| 黄色视频,在线免费观看| netflix在线观看网站| 久久久精品欧美日韩精品| 在线av久久热| 日韩成人在线观看一区二区三区| 久久午夜亚洲精品久久| 色哟哟哟哟哟哟| 欧美激情高清一区二区三区| 搞女人的毛片| 最近最新中文字幕大全免费视频| 国产精品 欧美亚洲| 亚洲激情在线av| 亚洲精品国产一区二区精华液| 一区二区日韩欧美中文字幕| 日韩有码中文字幕| 国产亚洲精品久久久久久毛片| 男人的好看免费观看在线视频 | 精品一品国产午夜福利视频| 多毛熟女@视频| 欧美午夜高清在线| 国产在线精品亚洲第一网站| 久久国产精品男人的天堂亚洲| 欧美另类亚洲清纯唯美| 亚洲av成人不卡在线观看播放网| 老司机午夜十八禁免费视频| 亚洲最大成人中文| 亚洲五月天丁香| 午夜激情av网站| 久久精品91无色码中文字幕| 欧美国产精品va在线观看不卡| 免费人成视频x8x8入口观看| 一级黄色大片毛片| 国产精品,欧美在线| 黄色a级毛片大全视频| 宅男免费午夜| 夜夜躁狠狠躁天天躁| 国产一级毛片七仙女欲春2 | 成人三级做爰电影| 黑人巨大精品欧美一区二区蜜桃| 欧美在线一区亚洲| 国产1区2区3区精品| 黑丝袜美女国产一区| 久久精品国产综合久久久| 欧美激情高清一区二区三区| 欧美日韩精品网址| 久久久精品欧美日韩精品| 亚洲成a人片在线一区二区| 久久人人精品亚洲av| 国产国语露脸激情在线看| 丝袜在线中文字幕| 久久久精品国产亚洲av高清涩受| 韩国精品一区二区三区| 国产单亲对白刺激| 久久欧美精品欧美久久欧美| 可以免费在线观看a视频的电影网站| 国产又色又爽无遮挡免费看| 高清在线国产一区| 午夜a级毛片| 给我免费播放毛片高清在线观看| 少妇熟女aⅴ在线视频| 视频在线观看一区二区三区| 免费在线观看日本一区| 午夜福利免费观看在线| 亚洲自偷自拍图片 自拍| 免费人成视频x8x8入口观看| 曰老女人黄片| 午夜福利视频1000在线观看 | 欧美成狂野欧美在线观看| 久久人人97超碰香蕉20202| 精品欧美国产一区二区三| 伊人久久大香线蕉亚洲五| 91av网站免费观看| 激情在线观看视频在线高清| 极品教师在线免费播放| 又黄又爽又免费观看的视频| 国产精品久久久久久人妻精品电影| 黄网站色视频无遮挡免费观看| 18禁黄网站禁片午夜丰满| 国产99白浆流出| 国内精品久久久久久久电影| 好男人电影高清在线观看| 久久久久国内视频| 12—13女人毛片做爰片一| 国产亚洲欧美98| 亚洲精品国产区一区二| 成人手机av| 日本在线视频免费播放| АⅤ资源中文在线天堂| 一本综合久久免费| 亚洲性夜色夜夜综合| 国产亚洲精品av在线| 亚洲av第一区精品v没综合| 日日爽夜夜爽网站| 少妇裸体淫交视频免费看高清 | 亚洲成人免费电影在线观看| 操出白浆在线播放| 久久国产精品男人的天堂亚洲| 日韩欧美国产在线观看| 日韩成人在线观看一区二区三区| www.熟女人妻精品国产| 色综合亚洲欧美另类图片| 色综合欧美亚洲国产小说| 欧美午夜高清在线| 91大片在线观看| 淫秽高清视频在线观看| 看片在线看免费视频| 日韩欧美免费精品| 在线观看免费视频日本深夜| 久久久久九九精品影院| 自线自在国产av| 亚洲免费av在线视频| 少妇被粗大的猛进出69影院| 亚洲欧美日韩高清在线视频| 亚洲av成人av| 夜夜夜夜夜久久久久| av中文乱码字幕在线| 亚洲第一电影网av| 亚洲久久久国产精品|