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

    大通湖浮游植物群落結構與環(huán)境因子關系

    2017-01-10 05:12:48高子涵孟憲智周緒申
    水利信息化 2016年5期
    關鍵詞:浮游大通富營養(yǎng)化

    高子涵,孟憲智,周緒申

    (海河流域水環(huán)境檢測中心,天津 300170)

    大通湖浮游植物群落結構與環(huán)境因子關系

    高子涵,孟憲智,周緒申

    (海河流域水環(huán)境檢測中心,天津 300170)

    湖南省大通湖是洞庭湖區(qū)典型的淺水養(yǎng)殖湖泊,為合理地開發(fā)利用其自然資源,保護生態(tài)環(huán)境,于2013 年 4 月至 2014 年 3 月對大通湖浮游生物群落結構進行了調查研究,并通過典范對應分析(CCA),探討浮游植物群落結構與環(huán)境因子之間的關系。調查期間共發(fā)現浮游植物 7 門 80 種,全年優(yōu)勢種為多形裸藻 Euglena polymorphya、卵形隱藻 Cryptomonas ovate、尖尾藍隱藻 Chroomonas acuta、鏈狀小環(huán)藻 Syclotella catenata、小球藻 Chlorella vulgaris、圓形衣藻 Chlamydomonas orbicularis 及四尾柵藻 Scenedesmus quadricauda。浮游植物密度和生物量年平均值分別為 6.58 ± 0.66 × 104ind./L 和 0.16 ± 0.05 mg/L,且空間上均呈現南高北低的分布趨勢。浮游植物多樣性指數及豐富度指數均表明大通湖水體呈富營養(yǎng)狀態(tài)。CCA 表明水溫、總磷及 pH 是影響大通湖浮游植物群落結構的主要環(huán)境因子。

    浮游植物;浮游生物;多樣性;群落結構;環(huán)境因子;典范對應分析

    浮游植物通過光合作用制造有機物,同時為水體提供大量的溶解氧,是水生生態(tài)系統(tǒng)的初級生產者,是食物網的基礎環(huán)節(jié),在生態(tài)系統(tǒng)中有著極為重要的地位,同時還參與了水生生態(tài)系統(tǒng)的能量流動、物質循環(huán)及信息的傳遞[1]。浮游植物的種類組成、群落結構及現存量能對水環(huán)境條件的變化做出復雜而快速的響應,能直接影響水生食物網,并對水質狀況起到指示作用,從而在水體尤其是湖泊、水庫水體監(jiān)測中起極為重要的作用[2-3]。此外,在水產養(yǎng)殖方面,浮游植物是水生經濟動物的直接或間接餌料,浮游植物的群落結構特征能為水產品的合理放養(yǎng)提供依據;在水環(huán)境監(jiān)測方面,利用浮游植物對水環(huán)境的敏感性可以廣泛地用于水質監(jiān)測和評價[4]。浮游植物群落結構的動態(tài)變化是多個環(huán)境因子在時間和空間上共同作用的結果,因此浮游植物群落結構的動態(tài)變化與環(huán)境因子關系是一個相當復雜的問題,而典范相應分析能很好地反映出浮游植物群組成與環(huán)境因子之間對應關系,同時結合多個環(huán)境因子分析,包含的信息量大,結果直觀明顯, 從而更好地反映出主要浮游植物類群與環(huán)境因子之間的對應關系[5-6]。

    大通湖(29°09'~29°15'N,112°26'~112°33'E)位于長江中游南岸,湘中偏北,東鄰漉湖,南與沅江市相連,西北與南縣、華容縣比鄰,20 世紀 50 年代初由于圍墾從洞庭湖隔離開,是洞庭湖區(qū)典型的淺水型養(yǎng)殖湖泊,同時也是湖南省面積最大的養(yǎng)殖湖泊,漁業(yè)資源極其豐富[7]。大通湖水產養(yǎng)殖主要為投放四大家魚,且以鰱、鳙為主,采取人放天養(yǎng)并適當施肥投餌的半精養(yǎng)養(yǎng)殖模式。大通湖水面面積 82.67 km2,東西長 15.75 km,南北寬 13.70 km,湖盆呈三角形,湖底高程為海拔 23.8~26.0 m(以吳淞高程為基準)。該湖泊被人為地用攔網分隔成蜜蜂夾湖、大西湖和尼古湖 3 部分[8]。20 世紀 60 年代,盧奮英等[9]對自然狀態(tài)下大通湖水體理化環(huán)境和水生生物進行了全面的調查分析。之后的 50 年內,大通湖漁業(yè)生產方式經歷了由自然增殖到人放天養(yǎng)再到人工養(yǎng)殖的發(fā)展歷程[10]。近年來,隨著水生生物經濟價值的提高,天然養(yǎng)殖湖泊中的捕撈強度日益加大,大通湖水生生物的捕撈量也在逐年加強。但是與此同時隨著捕撈強度的加大及人工施肥使大通湖水體營養(yǎng)狀態(tài)由中營養(yǎng)向富營養(yǎng)轉變[11-12]。

    本研究于 2013 年 4 月至 2014 年 3 月逐月對大通湖浮游植物群落結構進行周年采樣調查,并通過典范對應分析,探討大通湖浮游植物群落結構特征與環(huán)境因子之間的關系,旨在了解掌握大通湖浮游植物群落結構的時空分布變化,評價大通湖水體的健康狀況,并為今后的資源保護及開發(fā)利用提供理論基礎,以期為大通湖乃至長江中下游同類型水體的漁業(yè)可持續(xù)利用及有效管理提供參考。

    1 材料與方法

    1.1 采樣點設置

    大通湖水體近乎封閉,周邊生態(tài)環(huán)境優(yōu)良,無集中的點源污染,養(yǎng)殖物化產品是潛在的污染源[13]。本研究依據湖泊形態(tài)及主要養(yǎng)殖區(qū)域共設置 9 個固定采樣點,采樣點均設置在圍網、湖內簡易建筑物處,以減少施肥投餌、漁業(yè)捕撈對水質和水生生物資源調查結果造成的影響[14],如圖 1 所示。

    圖1 采樣點分布圖

    1.2 樣品采集與處理

    1.2.1 水質樣品采集與處理

    水深(Depth)和透明度(SD)用透明度盤原位測定。水溫(T)、pH 采用 HANNA 多功能水質分析儀原位監(jiān)測??偟═N)、總磷(TP)、硝酸鹽氮可溶性磷(及葉綠素 a(Chl-a)測定所需水樣用 2.5 L 采水器采集表層 0.5 m 水深處水樣,置于聚乙烯瓶中帶回實驗室進行分析。TN 用堿性過硫酸鉀消解紫外分光光度法測定,用紫外分光光度法測定,用鉬酸銨分光光度法測定。葉綠素 a 含量測定采用丙酮萃取分光光度法測定。水樣采集及測定方法均參照《水和廢水監(jiān)測分析方法》[15]進行。

    1.2.2 浮游植物的采集與處理

    浮游植物定性定量樣品采集參照《湖泊富營養(yǎng)化調查規(guī)范(第二版)》[16]及《內陸水域漁業(yè)自然資源調查手冊》[17]進行。

    浮游植物定性樣品采集使用 25# 浮游生物網,關閉網底收集管閥門,在水體表層至 0.5 m 深處呈“∞”字型在采集水域緩慢循回拖拽采集,速度保持在 30 cm/s,采樣過程中網口與水面始終保持垂直,且網口上端不能露出水面。采集完成收網后用水樣將浮游生物網沖洗干凈,然后將采集到的樣品收集到塑料瓶中,貼好標簽,加入 1.0%~1.5% 體積的魯哥試劑固定樣品。

    浮游生物不僅在水平分布上存在差異,而且在垂直分布上也存在差異,大通湖為淺水湖,因此在本研究中并不涉及垂直分布。浮游植物定量樣品在距水體表面 0.5 m 處采集,用有機玻璃采水器采集水樣,放入 1 000 mL 塑料瓶中,貼好標簽,加入 1.0%~1.5% 體積的魯哥試劑進行現場固定。將定量樣品移入沉淀器中,靜置 24~36 h 后,用虹吸管小心吸取上清液,并重復多次直至濃縮至 20~30 mL。

    1.2.3 浮游植物的鑒定與計數

    定性樣品鑒定直接用膠頭滴管移取樣本滴于載玻片上,蓋上蓋玻片,在 10×40 倍顯微鏡下觀察,做好記錄。浮游植物種類鑒定使用 ZEISS-Imager.M2顯微鏡,參照《中國淡水藻類——系統(tǒng)、分類及生態(tài)》[18]及《中國淡水生物圖譜》[19],定性鏡檢到種。

    定量樣品計數時,首先將樣本充分搖勻,用移液器移取 0.1 mL 樣本于 0.1 mL 計數框中,緩慢蓋上蓋玻片,在 10×40 倍顯微鏡下全片觀察,記下每種浮游植物的個數。同一采樣點水樣觀察計數 2 片,取其平均值,若同一樣品的 2 次計數結果之差大于± 15%,則應再增加計數次數。

    每升水樣中浮游植物數量的計算公式為

    式中:N 為水樣中浮游植物的總個數;Pn為計數所得個數;υ 為計數框容積;V1為濃縮后水樣體積;V2為去掉液體積。

    大多數藻類細胞形態(tài)較為規(guī)則,且細胞的比重接近于 1,因此浮游植物生物量計算依照體積法進行,即在鏡檢過程中對不同大小的藻細胞長度、寬度、直徑或高度等參數進行測量,求得各個參數的平均值后用相關公式算出體積,換算為生物量。對于不規(guī)則的形體可分割為多個部分進行測量,并用相似圖形計算體積,直接換算為生物量[16]。

    1.3 數據處理

    利用 Excel 2013 進行數據的統(tǒng)計,SPSS,Statistica 6.0,Canoco 進行數據的分析處理,Origin Pro 8 進行相關圖形的繪制。

    1.3.1 Shannon 多樣性指數

    Shannon 多樣性指數反映了群落中物種的多樣性,其值越大,表明群落結構越復雜程度, H 值在0~2 為富營養(yǎng),2~3 時為中營養(yǎng),> 3 為貧營養(yǎng)[20]。計算公式為

    式中:S 為樣品中總種類數;N 為樣品中浮游植物總個體數;ni為樣品中第 i 種浮游植物個體數。

    1.3.2 Pielou 均勻度指數

    Pielou 均勻度指數主要描述群落中各物種之間個體數目分配的均勻性,J 值在 0~0.3 為富營養(yǎng)化,0.3~0.5 為中營養(yǎng)化,0.5~0.8 為貧營養(yǎng)化,0.8 ~ 1.0 為清潔水體[22]。計算公式為

    式中: S 為樣品中總種類數。

    1.3.3 Margalef 豐富度指數

    Margalef 豐富度指數主要反映了群落物種的豐富度,當 D > 5 時為貧營養(yǎng)水體,在 4~5 之間時為中營養(yǎng)水體,當 D < 4 之間時為富營養(yǎng)水體[22]。計算公式為

    式中:S 為樣品中總種類數;N 為樣品中浮游植物總個體數。

    1.3.4 優(yōu)勢度

    評價浮游動物優(yōu)勢種的優(yōu)勢度。Y 值 > 0.02 的種類為全年調查的優(yōu)勢種[25]。計算公式為

    式中:N 為樣品中浮游植物總個體數;ni為樣品中第 i 種浮游植物個體數;?i為該種在各站點出現的頻率。

    1.3.5 典范對應分析(CCA)

    典范對應分析(CCA)采用 Canoco for Windows 4.5 軟件對浮游生物物種和環(huán)境因子數據進行分析。水環(huán)境因子以及浮游植物數據矩陣經過 lg(x + 1)轉換[27]。

    2 結果與分析

    2.1 大通湖浮游植物群落結構

    調查顯示,大通湖 9 個監(jiān)測點全年共檢測浮游植物 7 門 80 種,其中綠藻門 (Chlorophyta) 26 種,占浮游植物總種類數的 32.5%,硅藻門 (Bacillariophyta) 21 種,占 26.3%,藍藻門 (Cyanophyta) 19 種,占 23.8%,裸藻門 (Euglenophyta) 6 種,占 7.5%,黃藻門 (Xanthophyta) 4 種,占 5.0%,隱藻門(Cryptophyta) 3 種,占 3.8%,甲藻門 (Pyrrophyta) 1 種,占 1.3%。浮游植物種類組成如圖 2 所示。大通湖浮游植物種類繁多,其中,綠藻門、硅藻門和藍藻門在物種數目上占有絕對優(yōu)勢。綠藻和硅藻的種類組成變化相對較大,是引起大通湖浮游植物群落結構變化的主要原因之一。監(jiān)測期間綠藻屬數和種類數始終占首位,硅藻位居第二,根據浮游植物組成特點劃分,大通湖屬于綠藻-硅藻型湖泊。

    圖2 大通湖浮游植物種類組成

    浮游植物種類組成時間分布如圖 3 所示,2013 年8 月物種最為豐富為 29 種,其次為 2014 年 1 月,共檢出 28 種,而在 2013 年 6 月物種數最低,為 12 種,不同月份浮游植物種類組成如表 1 所示。從全年檢測結果來看,綠藻門和硅藻門物種數目在全年每月均占有絕對優(yōu)勢。

    浮游植物種類組成空間分布較為平均,如圖 4所示,監(jiān)測點 S4 浮游植物種類數最為豐富,共檢出49 種;其次為監(jiān)測點 S6,為 47 種;物種數最少的監(jiān)測點為 S9,檢出浮游植物 31 種。不同監(jiān)測點浮游植物種類組成如表 2 所示。從全年檢測結果來看,綠藻門和硅藻門物種數目在各個監(jiān)測點均占有絕對優(yōu)勢。

    圖3 浮游植物種類組成時間分布

    2.2 大通湖浮游植物現存量

    2.2.1 浮游植物密度

    大通湖浮游植物全年平均密度為 6.58 ± 0.66× 104ind./L。如圖 5 所示,2013 年 12 月浮游植物密度最高,為 18.31 ± 1.82×104ind./L,其次為 2013 年8 月,浮游植物密度為 10.16 ± 5.02×104ind./L,2013 年 7 月浮游植物密度最低,為 2.71 ± 0.60× 104ind./L。如圖 6 所示,監(jiān)測點 S3 年平均浮游植物密度最高,為 9.28 ± 4.06×104ind./L,其次為監(jiān)測點S6,浮游植物密度為 8.32 ± 1.69×104ind./L,監(jiān)測點S9 浮游植物密度最低,為 4.34 ± 0.89×104ind./L。

    表1 不同月份浮游植物種類組成

    圖4 浮游植物空間分布

    2.2.2 浮游植物生物量

    大通湖浮游植物全年平均生物量為 0.16 ± 0.05 mg/L。如圖 7 所示,2013 年 11 月浮游植物生物量最高,為 0.64 ± 0.26 mg/L,其次為 2013 年 4 月,浮游植物生物量為 0.26 ± 0.09 mg/L,2014 年 1 月浮游植物生物量最低,為 0.04 ± 0.01 mg/L。如圖 8 所示,監(jiān)測點 S3 年平均浮游植物生物量最高,為 0.36± 0.21 mg/L,其次為監(jiān)測點 S6,浮游植物生物量為0.22 ± 0.11 mg/L,監(jiān)測點 S2 與 S7 浮游植物生物量最低,均為 0.10 ± 0.02 mg/L。

    2.3 大通湖浮游植物優(yōu)勢種

    大通湖全年調查浮游植物的優(yōu)勢種為裸藻門的多形裸藻 Euglena polymorphya,隱藻門的卵形隱藻Cryptomonas ovate 和尖尾藍隱藻 Chroomonas acuta,硅藻門的鏈狀小環(huán)藻 Syclotella catenata,綠藻門的小球藻 Chlorella vulgaris、圓形衣藻 Chlamydomonas orbicularis 及四尾柵藻 Scenedesmus quadricauda,各種類的優(yōu)勢度(Y)如表 3 所示。此外,尖針桿藻Synsdra acus 雖然優(yōu)勢度值(Y)< 0.02,但其全年出現頻率較高,為61.11%。

    表2 不同采樣點浮游植物種類組成

    圖5 浮游植物密度時間分布

    圖6 浮游植物密度空間分布

    圖7 浮游植物生物量時間分布

    圖8 浮游植物生物量空間分布

    2.4 大通湖浮游植物多樣性指數

    2.4.1 Shannon 多樣性指數

    大通湖全年浮游植物 Shannon 多樣性指數(H)年平均值為 1.43,其中,2013 年 8 月為全年最高值,H 值為 2.08,2013 年 11 月最低,為 0.72,如圖 9 所示。根據 Shannon 多樣性指數評價水體富營養(yǎng)化程度:指數值為 0~2 時,水體呈富營養(yǎng)化;指數值為 2~3 時,為中營養(yǎng)化;指數值 > 3 時,為貧營養(yǎng)化。大通湖浮游植物 Shannon 多樣性指數在0.72~2.08 之間,年平均值為 1.43,呈富營養(yǎng)化。大通湖全年浮游植物 Shannon 多樣性指數(H)呈現 H夏> H冬> H春> H秋的趨勢。

    表3 大通湖浮游植物優(yōu)勢種

    圖9 大通湖浮游植物多樣性指數 H 周年變化

    2.4.2 Pielou 均勻度指數

    大通湖全年浮游植物 Pielou 均勻度指數(J)在0.31~0.84 之間,全年平均值為 0.68。其中 2013 年5 月 J 值最高,為 0.86;2013 年 11 最低,為 0.31,如圖 10 所示。根據 Pielou 均勻度指數評價水體富營養(yǎng)化程度:指數值為 0.8~1.0,水體為清潔型水體;指數值為 0.5~0.8,水體呈貧營養(yǎng)化;指數值為 0.3~0.5,水體呈中營養(yǎng)化;指數值在 0~0.2 之間時,水體呈富營養(yǎng)化。大通湖浮游植物 Pielou 均勻度指數變化范圍為 0.31~0.84,年平均值為 0.68,呈貧營養(yǎng)化。大通湖全年浮游植物 Pielou 均勻度指數(J)呈現 J夏>J春>J冬>J秋的趨勢。

    2.4.3 Margalef 豐富度指數

    大通湖全年浮游植物 Margalef 豐富度指數(D)在 2013 年 8 月達到最高值為 1.15,在 2013 年4 月最低,為 0.38,年平均值為 0.75,如圖 11 所示。根據 Margalef 豐富度指數評價水體富營養(yǎng)化程度:指數值為 0~4,水體呈富營養(yǎng)化;指數值為 4~5,呈中營養(yǎng)化,指數值 > 5 時,呈貧營養(yǎng)化。大通湖浮游植物 Margalef 豐富度指數在 0.38~1.15 之間,年平均值為 0.75,呈富營養(yǎng)化。大通湖全年浮游植物 Margalef 豐富度指數(D)呈現D春的趨勢。

    圖10 大通湖浮游植物均勻度指數 J 周年變化

    圖11 大通湖浮游植物豐富度指數 D

    2.5 大通湖浮游植物與環(huán)境因子關系

    選取大通湖浮游植物全年優(yōu)勢種為典范對應分析中的主要代表物種。各物種在 CCA 分析圖中的代碼如表 4 所示。

    表4 CCA 分析浮游植物種類

    為了探討浮游植物與環(huán)境因子的關系,對浮游植物群落的密度與水溫、水深、pH、總氮、總磷、葉綠素、浮游動物密度及浮游動物生物量這 8 項環(huán)境指標做典范對應相關性分析,分析結果如表 5 所示。

    表5 環(huán)境因子與前兩個排序軸的相關系數

    典范對應分析排序圖 12 中前 2 個排序軸的特征根值分別為 0.186 和 0.051,浮游植物種類與環(huán)境因子前 2 個排序軸的相關系數分別為 0.790 和 0.593。前 2 軸共解釋浮游植物種類與環(huán)境關系的 83.1%,如表 6 所示。在典范對應分析排序中,物種和環(huán)境因子相關系數均為 1,表明該排序結果能很好地反映浮游植物與環(huán)境因子間的關系。物種和環(huán)境因子的 2 個排序軸之間的相關系數均分別為 0,表明該排序結果準確、可靠。

    經典范對應分析可知,水溫、pH 及總磷是第一排序軸的重要影響因子。其中,水溫、總磷與第一排序軸呈負相關,相關系數分別為 -0.528 9 和-0.443 5,pH 與第一排序軸呈正相關,相關系數為0.634 3。而其他因子與第一排序軸相關性較小,對浮游植物群落結構影響較小。因此,對浮游植物群落結構影響較大的主要環(huán)境因子是水溫、pH 及總磷。

    圖12 大通湖浮游植物與環(huán)境因子關系 CCA 分析排序圖

    表6 大通湖浮游植物物種和環(huán)境因子間 CCA 分析的統(tǒng)計信息

    3 分析與討論

    3.1 大通湖浮游植物群落結構

    大通湖浮游植物密度的變化很不穩(wěn)定,變化幅度大, 除了季節(jié)性影響外,可能與近年來大通湖養(yǎng)殖密度提高施肥投餌量加大所致。螺、蚌、蜆、魚對浮游植物消耗量過大也是浮游植物密度不規(guī)律變化的主要原因。此外,根據鐘詩群[28]等報道,大通湖在 5—10 月實行間歇施肥,水體中的氮、磷含量始終維持在較高的水平。因此當光照充分時, 浮游植物迅速繁殖, 數量急劇升高。浮游植物密度的最低值出現鰱、鳙在生長季節(jié),鰱、鳙濾食是導致該現象發(fā)生的主要原因之一。

    浮游植物密度在空間分布上呈現從南到北逐漸下降趨勢,這種分布特征與軟體動物生物量由東北向西南逐漸遞減的空間分布格局相反[29]。研究表明,背角無齒蚌 Anodonta woodiana 能顯著降低浮游植物的現存量,同時改變浮游植物的群落結構,降低藍藻門所占比例,提升綠藻門所占比例[30]。此外,河蜆對浮游植物的群落結構也具有間接影響。朱小龍等[31]研究表明,河蜆濾食能顯著降低水體懸浮物濃度及葉綠素 a 含量。大通湖螺、蚌、蜆資源總量十分巨大,據調查,2009 年全湖共捕撈螺類 33 000 t,河蜆 515 t。大通湖水體 TN,TP 含量較高(TN 年平均含量為 4.93 mg/L,TP 年平均含量為 0.69 mg/L)處于富營養(yǎng)狀況,具備水華發(fā)生的營養(yǎng)鹽條件(TN 為 0.1~0.8 mg/L,TP 為 2.5~3.5 mg/L)[32]。大通湖浮游植物密度較其他富營養(yǎng)水體偏低,并且近年的記錄中沒有出現過典型的水華,巨大的軟體動物資源量從一定程度上抑制了藻類水華的發(fā)生,而且使大通湖浮游植物密度長期維持在低于其他同類富營養(yǎng)水體的水平。

    3.2 大通湖浮游植物多樣性

    Shannon 多樣性指數(H)通常用于反映浮游生物群落結構的復雜程度,Margalef 豐富度指數(D)能充分反映浮游生物種類的分布情況,Pielou 均勻度指數(J)反映各浮游生物物種個體密度分配的均勻程度,在選擇多樣性指數進行評價時應選擇多個指數[33]。但在使用浮游植物 Shannon 多樣性、Margalef豐富度及 Pielou 均勻度等指數評價大通湖水體富營養(yǎng)化程度時,結果出現不一致的現象。在評價水體富營養(yǎng)化程度上,浮游植物的 Shannon 多樣性、Margalef 豐富度及 Pielou 均勻度等指數不存在敏感性和準確性等差異,而是各種指數的評價標準需要根據不同類型的水體進行適當調整,且生物多樣性指數法具有相對穩(wěn)定性、長效性和滯后性,因此不排除 Pielou 均勻度指數并不適用于評價大通湖水體的可能性。此外,大通湖浮游植物 Shannon 多樣性與 Margalef 豐富度指數的變化趨勢基本一致,浮游植物多樣性指數的變動,是物種豐富度的變化,即浮游植物種類數變化的綜合反映。

    3.3 環(huán)境因子對浮游植物群落結構特征的影響

    典范對應分析表明大通湖浮游植物密度的變化與水溫、pH 及總磷含量密切相關。研究表明,溫度是影響浮游植物密度的主要因素[34-35]。溫度是影響浮游植物空間和時間分布的重要因素,不同的浮游植物種類有其生長繁殖的最佳溫度范圍[36]。大通湖浮游植物群落結構主要由耐低溫的硅藻構成,耐低溫的硅藻與溫度呈負相關,因此大通湖浮游植物密度呈現出隨著溫度的上升而下降的趨勢,并在 2013 年12 月達到最高值。而耐高溫綠藻與溫度呈良好的正相關,在夏季迅速生長,且密度隨著溫度的升高而增加。

    研究表明,營養(yǎng)鹽是影響浮游植物現存量最基本的因素[37]。本研究中,CCA 分析中顯示總磷較總氮對浮游植物密度的影響更大,與之前李德亮等人[10]對大通湖浮游植物群落結構與環(huán)境因子關系的研究結果一致。

    水體 pH 對浮游植物現存量的影響可能是由于水體中營養(yǎng)物質含量的增加,導致水生生物群落結構發(fā)生變動,各營養(yǎng)級浮游生物現存量發(fā)生變化,從而反過來對水體理化指標造成影響。

    此外,大通湖為典型的內陸淺水型湖泊,水體面積較小,水深較淺,因此透明度和水深對浮游植物密度分布的影響并不明顯。

    大通湖浮游植物與環(huán)境因子 CCA 排序圖反映出不同種類浮游植物對水環(huán)境的適應性不同。四尾柵藻和多形裸藻在 CCA 排序圖中的位置相對接近,表明其在生態(tài)環(huán)境適應性上較為相似的。此外,根據CCA 排序圖中的位置可以看出,鏈狀小環(huán)藻與浮游動物生物量密切相關,尖尾藍隱藻的分布主要受到水體 pH 影響,原因還有待進一步驗證。

    4 結語

    1)大通湖浮游植物主要以硅藻和綠藻絕對占主導地位,因此大通湖浮游植物群落結構的周年性和季節(jié)性變化趨勢的差異主要由硅藻和綠藻引起的。此外,浮游植物密度的空間分布呈現從南到北逐漸下降趨勢。

    2)大通湖浮游植物 Shannon 多樣性指數與Margalef 豐富度指數均表明,大通湖水體呈和營養(yǎng)化。

    3)典范對應分析表明大通湖浮游植物密度的變化與水溫、pH 及總磷含量密切相關。

    [1] 章宗涉,黃祥飛. 淡水浮游生物研究方法[M]. 北京:科學出版社,1991: 333-371.

    [2] 譚香,夏小玲,程曉莉,等. 丹江口水庫浮游植物群落時空動態(tài)及其多樣性指數[J]. 環(huán)境科學,2011,32 (10): 2875-2882.

    [3] 王曉星. 北京翠湖濕地浮游植物群落特征[J]. 濕地科學與管理,2013,9 (2): 57-60.

    [4] 劉健康. 高級水生生物學[M]. 北京:科學出版社,1999: 52-128.

    [5] 張金屯. 植被與環(huán)境關系的分析Ⅱ:CCA 和 DCCA 限定排序[J]. 山西大學學報(自然科學版),1992,15 (3): 292-298.

    [6] 李秋華, 韓博平. 基于 CCA 的典型調水水庫浮游植物群落動態(tài)特征分析[J]. 生態(tài)學報,2007,27 (6): 2355-2364.

    [7] 張健,高子涵,馬行厚,等. 大通湖河蜆(Corbicula fluminea)干重的換算方法[J]. 湖南農業(yè)大學學報(自然科學版),2014,40 (6): 633-636.

    [8] 李德亮,張婷,余建波,等. 長江中游典型湖泊重金屬分布及其風險評價——以大通湖為例[J]. 長江流域資源與環(huán)境,2010,19(增刊 1): 183-189.

    [9] 盧奮英, 章宗涉, 丘昌強,等. 大通湖的生物相及理化環(huán)境與近代沉積的初步報告[J]. 海洋與湖沼,1964,6 (4): 323-330.

    [10] 李德亮,張婷,肖調義,等. 大通湖浮游植物群落結構及其與環(huán)境因子關系[J]. 應用生態(tài)學報,2012,23 (8): 2017-2113.

    [11] 金宏,向建國,唐文熙. 大通湖漁場浮游生物及其魚產力的調查[J]. 內陸水產, 2005 (7): 22-24.

    [12] 李德亮,張婷,余建波,等. 典型養(yǎng)殖湖泊大通湖軟體動物的時空分布格局[J]. 水生生物學報,2011,35 (6): 946-954.

    [13] 高子涵,張健,皮杰,等. 湖南省大通湖浮游動物群落結構及其與環(huán)境因子關系[J]. 生態(tài)學雜志,2016,35 (3): 733-740.

    [14] 王桂蘋,李雪林,皮杰,等. 大通湖環(huán)棱螺的次級生產力[J]. 水生生物學報, 2014,38 (5): 987-993.

    [15] 國家環(huán)境保護總局. 地表水和廢水監(jiān)測技術規(guī)范:HJ/ T 91-2002 [S]. 北京:中國標準出版社,2002: 1-47.

    [16] 金相燦,屠清瑛. 湖泊富營養(yǎng)化調查規(guī)范[S]. 2 版. 北京:中國環(huán)境科學出版社,1990: 21-25.

    [17] 張覺民,何志輝. 內陸水域漁業(yè)自然資源調查手冊[M].北京:農業(yè)出版社,1991: 12-122.

    [18] 胡鴻鈞, 魏印心. 中國淡水藻類—系統(tǒng)、分類及生態(tài)[M]. 北京:科學出版社,2006: 21-915.

    [19] 韓茂森,束蘊芳. 中國淡水生物圖譜[M]. 北京:海洋出版社,1995: 2-133.

    [20] 林碧琴,謝淑琦. 水生藻類與水體污染監(jiān)測[M]. 沈陽:遼寧大學出版社,1988: 227-234.

    [21] Shannon C.E.. Amathematical theory of communication[J]. Bell System Technical Journal,1948 (27): 379-423.

    [22] 沈韞芬,章宗涉,龔循矩,等. 微型生物監(jiān)測新技術[M].北京:中國建筑工業(yè)出版社,1990: 119-130.

    [23] Pielou E.C.. The measurement of diversity in different types of biological collections[J]. Journal of Theoretical Biology,1966 (13): 131-144.

    [24] Margalef R. Information theory in ecology[J].General Systematics,1958 (3): 36-71.

    [25] 張婷,李林,宋立榮. 熊河水庫浮游植物群落結構的周年變化[J]. 生態(tài)學報,2009,29 (6): 2971-2979.

    [26] Gao Y.,Jiang Z.B.,Liu J.J.,et al.. Seasonal variations of net-phytoplankton community structure in the Southern Yellow Sea[J]. Oceanic and Coastal Sea Research,2013 (12): 557-567.

    [27] 邱小琮,趙紅雪,孫曉雪. 沙湖浮游動物與水環(huán)境因子關系的多元分析[J]. 生態(tài)學雜志,2012,31 (4): 896-901.

    [28] 鐘詩群,陳榮坤,周杰,等. 大通湖浮游植物群落結構與富營養(yǎng)化動態(tài)研究[J]. 生態(tài)科學,2014,33 (3): 586-593.

    [29] 李德亮,張婷,肖調義,等. 湖南省大通湖河蜆(Corbicula fluminea)現存量及其時空分布[J]. 湖泊科學,2013,25 (5): 743-748.

    [30] 楊東妹,陳宇煒,劉正文,等. 背角無齒蚌濾食對營養(yǎng)鹽和浮游藻類結構影響的模擬[J]. 湖泊科學,2008,20 (2): 228-234.

    [31] 朱小龍,谷嬌,靳輝,等. 太湖河蜆(Corbicula fluminea)對富營養(yǎng)水體水質的改善作用[J]. 湖泊科學,2015,27 (3): 486-492.

    [32] 朱廣偉. 太湖富營養(yǎng)化現狀及原因分析[J]. 湖泊科學,2008,20(1): 21-26.

    [33] 邱小琮,趙紅雪. 寧夏沙湖浮游植物群落結構及多樣性研究[J]. 水生態(tài)學雜志,2011,32 (1): 20-26.

    [34] 胡韌,熊江霞,韓博平. 具有短水力滯留的小型富營養(yǎng)化水庫浮游植物群落結構與動態(tài)[J]. 生態(tài)環(huán)境,2008,17 (4): 1319-1326.

    [35] 薄芳芳,楊虹,左倬,等. 上海公園水體夏季浮游植物群落與環(huán)境因子的關系[J]. 生態(tài)學雜志,200,28 (7): 1259-1265.

    [36] 陳立婧,吳竹臣,胡忠軍,等. 上海崇明島明珠湖浮游植物群落結構[J]. 應用生態(tài)學報,2011,22 (6): 1599-1605.

    [37] 吳京洪,楊秀環(huán),唐寶英,等. 大亞灣澳頭增養(yǎng)殖區(qū)赤潮與環(huán)境的關系研究 I. 浮游植物總生物量與環(huán)境因子的關系[J]. 中山大學學報(自然科學版),2001,40 (3): 37-40.

    Phytoplankton Community Structure and Its Relationships with Environmental Factors in Datong Lake

    GAO Zihan, MENG Xianzhi, ZHOU Xushen

    (Haihe River Water Environmental Monitoring Center, Tianjin 300170, China)

    Datong Lake is a typical shallow lake in Dongting Lake of Hunan province. In this study, plankton and environmental factors in Datong Lake of Hunan Province are investigated monthly from April 2013 to March 2014, in order to develop and utilize fishery resources rationally, protect the ecological environment. Plankton community structure and its relationships with environmental factors are analyzed by canonical correspondence analysis (CCA). Total 80 species from 7 genera of phytoplankton are identified. According to the abundance and occurrence ofphytoplankton, Euglena polymorphya, Cryptomonas ovate, Chroomonas acuta, Syclotella catenata, Chlorellavulgaris, Chlamydomonas orbicularis and Scenedesmus quadricauda are considered as the dominant species. The mean annual abundance and biomass of phytoplankton are 6.58±0.66×104ind./L and 0.16 ± 0.05 mg/L respectively. The abundance and biomass of phytoplankton have a spatial pattern of decreasing from South to North. The biodiversity index indicates that the lake is under accelerated eutrophication status. The CCA results indicate that water temperature, total phosphorus and pH are the main environmental factors affecting the phytoplankton community structure in the Datong Lake.

    phytoplankton; plankton; biodiversity; community structure; environmental factor; canonical correspondence analysis

    TV213

    A

    1674-9405(2016)05-0021-10

    10.19364/j.1674-9405.2016.05.004

    2016-08-18

    國家水體污染控制與治理科技重大專項(2012ZX07203-002)

    高子涵(1990-),女,天津人,助理工程師,從事水生態(tài)監(jiān)測工作。

    猜你喜歡
    浮游大通富營養(yǎng)化
    基于臨界點的杭州灣水體富營養(yǎng)化多年變化研究
    海洋通報(2022年4期)2022-10-10 07:41:48
    距離完美只差一點 上汽大通MAXUS G90
    車主之友(2022年4期)2022-08-27 00:56:02
    浞河浮游藻類的調查研究與水質評價
    洪口水庫近年富營養(yǎng)化程度時間分布的研究
    《浮游》
    流行色(2019年10期)2019-12-06 08:13:26
    尋找你心中的完美之車
    ——上汽大通D90……虞展
    洞庭湖典型垸內溝渠水體富營養(yǎng)化評價
    上汽大通:C2B個性化定制未來
    浮游藻類在水體PAHs富集與降解中的研究
    富營養(yǎng)化藻的特性與水熱液化成油的研究
    最好的美女福利视频网| 亚洲国产精品合色在线| 内地一区二区视频在线| 麻豆国产av国片精品| 一个人观看的视频www高清免费观看| 五月玫瑰六月丁香| 18禁国产床啪视频网站| 亚洲国产精品sss在线观看| 熟女电影av网| 激情在线观看视频在线高清| 国产精品亚洲美女久久久| 18+在线观看网站| 精品久久久久久久久久免费视频| 法律面前人人平等表现在哪些方面| 两人在一起打扑克的视频| 成人性生交大片免费视频hd| 国产精品久久久久久人妻精品电影| 国产精品日韩av在线免费观看| 午夜免费观看网址| 精华霜和精华液先用哪个| 国产精品影院久久| 日本黄大片高清| 99视频精品全部免费 在线| 亚洲国产精品合色在线| 内射极品少妇av片p| 成人一区二区视频在线观看| 国产久久久一区二区三区| 成年女人毛片免费观看观看9| 在线观看日韩欧美| 国产精品嫩草影院av在线观看 | 国产一区二区亚洲精品在线观看| 欧美黄色片欧美黄色片| 久久精品国产亚洲av香蕉五月| 韩国av一区二区三区四区| www.熟女人妻精品国产| 一区二区三区免费毛片| 十八禁网站免费在线| 深夜精品福利| 一进一出好大好爽视频| 一a级毛片在线观看| 亚洲国产精品999在线| 成人av一区二区三区在线看| 精品国产超薄肉色丝袜足j| 精品福利观看| 欧美最新免费一区二区三区 | 97超级碰碰碰精品色视频在线观看| 亚洲国产色片| 好男人在线观看高清免费视频| 日本一二三区视频观看| 俺也久久电影网| 18禁裸乳无遮挡免费网站照片| 黄片大片在线免费观看| 蜜桃久久精品国产亚洲av| 成人av一区二区三区在线看| 手机成人av网站| 女同久久另类99精品国产91| 国产黄片美女视频| 国语自产精品视频在线第100页| 色老头精品视频在线观看| 国产伦一二天堂av在线观看| 在线观看舔阴道视频| 内射极品少妇av片p| 一本综合久久免费| 又黄又粗又硬又大视频| 无限看片的www在线观看| 18+在线观看网站| 中文字幕高清在线视频| a级一级毛片免费在线观看| 欧美最黄视频在线播放免费| 精品久久久久久久久久免费视频| 黑人欧美特级aaaaaa片| 国产高清三级在线| 久久久久久久亚洲中文字幕 | 97超级碰碰碰精品色视频在线观看| 国产麻豆成人av免费视频| 可以在线观看的亚洲视频| x7x7x7水蜜桃| 亚洲成a人片在线一区二区| 国产v大片淫在线免费观看| 亚洲精品456在线播放app | 人人妻人人澡欧美一区二区| 日本与韩国留学比较| 午夜免费成人在线视频| 在线视频色国产色| 2021天堂中文幕一二区在线观| av中文乱码字幕在线| 久久国产乱子伦精品免费另类| 少妇人妻一区二区三区视频| 精品久久久久久,| 人人妻,人人澡人人爽秒播| 欧美日韩黄片免| 不卡一级毛片| 久久久久久九九精品二区国产| 亚洲成人久久性| 91字幕亚洲| 午夜影院日韩av| 色噜噜av男人的天堂激情| 国产蜜桃级精品一区二区三区| 他把我摸到了高潮在线观看| 国产高清有码在线观看视频| 偷拍熟女少妇极品色| 精品免费久久久久久久清纯| 免费人成视频x8x8入口观看| 精品久久久久久久末码| 99riav亚洲国产免费| 麻豆成人午夜福利视频| 午夜激情欧美在线| 欧美丝袜亚洲另类 | 午夜福利高清视频| 波多野结衣高清无吗| 亚洲自拍偷在线| 一区福利在线观看| 亚洲av美国av| 国产av不卡久久| 嫩草影院精品99| 狠狠狠狠99中文字幕| 男插女下体视频免费在线播放| 在线观看日韩欧美| 国产69精品久久久久777片| 精品午夜福利视频在线观看一区| 欧美国产日韩亚洲一区| 亚洲欧美日韩高清在线视频| 在线观看舔阴道视频| 一个人看的www免费观看视频| 精品一区二区三区人妻视频| 12—13女人毛片做爰片一| 国产精品久久久久久亚洲av鲁大| 亚洲国产精品999在线| 久久99热这里只有精品18| 波多野结衣高清作品| or卡值多少钱| 日本a在线网址| 一区二区三区免费毛片| tocl精华| 少妇人妻精品综合一区二区 | 在线观看66精品国产| 亚洲中文字幕一区二区三区有码在线看| 欧美一区二区精品小视频在线| 亚洲美女黄片视频| 午夜福利欧美成人| 国产91精品成人一区二区三区| 丰满人妻一区二区三区视频av | 麻豆国产av国片精品| 亚洲一区高清亚洲精品| 国产 一区 欧美 日韩| 99久久精品一区二区三区| 少妇的逼水好多| 欧美激情久久久久久爽电影| 国产高清videossex| 免费在线观看日本一区| 亚洲国产日韩欧美精品在线观看 | 乱人视频在线观看| 亚洲av一区综合| 亚洲精品影视一区二区三区av| 男人和女人高潮做爰伦理| 欧美日韩福利视频一区二区| 色综合欧美亚洲国产小说| 免费观看精品视频网站| 国产免费男女视频| 岛国在线观看网站| 午夜福利在线观看免费完整高清在 | netflix在线观看网站| 天天一区二区日本电影三级| 国产精品1区2区在线观看.| 一区二区三区激情视频| 亚洲精品在线观看二区| 国产蜜桃级精品一区二区三区| 综合色av麻豆| 99国产综合亚洲精品| 性色avwww在线观看| 少妇人妻精品综合一区二区 | 亚洲电影在线观看av| 国产三级中文精品| 久久精品国产亚洲av香蕉五月| 久久久久久久午夜电影| 欧美激情在线99| 欧美色视频一区免费| 亚洲av电影不卡..在线观看| 久久久精品大字幕| 99久久99久久久精品蜜桃| 欧美激情久久久久久爽电影| 成人国产一区最新在线观看| 五月伊人婷婷丁香| 久久香蕉国产精品| 日本免费a在线| 中国美女看黄片| 国产麻豆成人av免费视频| 丰满人妻熟妇乱又伦精品不卡| 国产探花在线观看一区二区| 无限看片的www在线观看| 精品一区二区三区人妻视频| 国产真人三级小视频在线观看| 免费高清视频大片| 国内毛片毛片毛片毛片毛片| 黄片大片在线免费观看| 欧美一级毛片孕妇| 久久久精品大字幕| 欧美日韩精品网址| 婷婷精品国产亚洲av在线| 国产91精品成人一区二区三区| 亚洲无线在线观看| 啪啪无遮挡十八禁网站| 一进一出抽搐动态| 黄色丝袜av网址大全| 免费av观看视频| 动漫黄色视频在线观看| 欧美在线黄色| 亚洲精品国产精品久久久不卡| 嫩草影院入口| 国产日本99.免费观看| 观看免费一级毛片| 久久久久国产精品人妻aⅴ院| 老汉色∧v一级毛片| 亚洲美女视频黄频| 欧美一区二区亚洲| 1024手机看黄色片| 日本在线视频免费播放| 国产亚洲精品av在线| 国产99白浆流出| 国产一区二区在线观看日韩 | 国模一区二区三区四区视频| av片东京热男人的天堂| 啦啦啦观看免费观看视频高清| 欧美不卡视频在线免费观看| 国产精品久久久久久亚洲av鲁大| 宅男免费午夜| 免费人成视频x8x8入口观看| 一级黄片播放器| 久久人妻av系列| 亚洲av五月六月丁香网| 2021天堂中文幕一二区在线观| 97超视频在线观看视频| 嫩草影院入口| 欧美3d第一页| 午夜免费男女啪啪视频观看 | 国产99白浆流出| 色噜噜av男人的天堂激情| 偷拍熟女少妇极品色| 国产精品嫩草影院av在线观看 | 成年版毛片免费区| 亚洲乱码一区二区免费版| 久久人人精品亚洲av| 国产国拍精品亚洲av在线观看 | 婷婷丁香在线五月| 中文资源天堂在线| 悠悠久久av| 午夜免费观看网址| 亚洲五月天丁香| 日本撒尿小便嘘嘘汇集6| 免费在线观看日本一区| 小说图片视频综合网站| 又黄又爽又免费观看的视频| 身体一侧抽搐| 久久中文看片网| 啪啪无遮挡十八禁网站| 成年女人毛片免费观看观看9| 一级毛片高清免费大全| 精品一区二区三区视频在线观看免费| 少妇丰满av| 日日干狠狠操夜夜爽| 国产真实乱freesex| 午夜精品一区二区三区免费看| 欧美不卡视频在线免费观看| 美女免费视频网站| 亚洲狠狠婷婷综合久久图片| 真人一进一出gif抽搐免费| 制服人妻中文乱码| 少妇高潮的动态图| 精华霜和精华液先用哪个| 亚洲av一区综合| 日本撒尿小便嘘嘘汇集6| 真人做人爱边吃奶动态| 久久香蕉精品热| 婷婷六月久久综合丁香| 91在线精品国自产拍蜜月 | 日韩免费av在线播放| 国产高清视频在线观看网站| 国产成人影院久久av| 国产成人啪精品午夜网站| 白带黄色成豆腐渣| 18禁黄网站禁片免费观看直播| 成人特级av手机在线观看| 国产精品嫩草影院av在线观看 | 日本黄色视频三级网站网址| 国产精品一及| 深爱激情五月婷婷| 日本撒尿小便嘘嘘汇集6| 人人妻人人看人人澡| 又粗又爽又猛毛片免费看| 变态另类丝袜制服| 一本久久中文字幕| 日本一本二区三区精品| 中文亚洲av片在线观看爽| 一区福利在线观看| 性色av乱码一区二区三区2| 色视频www国产| 国产伦一二天堂av在线观看| 日本与韩国留学比较| 欧美一区二区亚洲| 国产v大片淫在线免费观看| 日韩国内少妇激情av| 欧美av亚洲av综合av国产av| 我要搜黄色片| 51国产日韩欧美| 精品久久久久久久人妻蜜臀av| 日本免费一区二区三区高清不卡| 欧美黑人欧美精品刺激| 99riav亚洲国产免费| 在线播放国产精品三级| 精品福利观看| 最新中文字幕久久久久| 51午夜福利影视在线观看| 国产97色在线日韩免费| 成人永久免费在线观看视频| 美女大奶头视频| 69人妻影院| 国产精品久久电影中文字幕| 免费av不卡在线播放| 男人舔女人下体高潮全视频| 99国产精品一区二区蜜桃av| 久久草成人影院| 757午夜福利合集在线观看| 国内揄拍国产精品人妻在线| 黄色视频,在线免费观看| 可以在线观看的亚洲视频| 国产欧美日韩精品一区二区| 91久久精品电影网| 久久精品综合一区二区三区| 在线免费观看不下载黄p国产 | 久久99热这里只有精品18| 久久精品综合一区二区三区| 亚洲欧美激情综合另类| 国产精品亚洲美女久久久| 亚洲国产精品999在线| 日韩av在线大香蕉| 国产高清videossex| 亚洲电影在线观看av| 老汉色av国产亚洲站长工具| 亚洲人成电影免费在线| 免费人成在线观看视频色| 欧美+亚洲+日韩+国产| 极品教师在线免费播放| 老司机午夜福利在线观看视频| 欧美精品啪啪一区二区三区| 久久这里只有精品中国| a在线观看视频网站| 人妻丰满熟妇av一区二区三区| 亚洲专区中文字幕在线| 极品教师在线免费播放| 老司机福利观看| 国产成+人综合+亚洲专区| 欧美av亚洲av综合av国产av| 91在线精品国自产拍蜜月 | 国产一区在线观看成人免费| 国产精品久久久久久人妻精品电影| 色综合欧美亚洲国产小说| 国产黄片美女视频| 欧美黄色片欧美黄色片| bbb黄色大片| 18+在线观看网站| 一级a爱片免费观看的视频| 免费搜索国产男女视频| 真实男女啪啪啪动态图| 国产爱豆传媒在线观看| 亚洲性夜色夜夜综合| 国产中年淑女户外野战色| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人系列免费观看| 欧美成人a在线观看| 午夜精品一区二区三区免费看| 一进一出抽搐gif免费好疼| 丰满乱子伦码专区| 91麻豆av在线| 成人特级黄色片久久久久久久| 少妇的逼好多水| 99久久精品国产亚洲精品| 搡老熟女国产l中国老女人| 搡老岳熟女国产| 一区二区三区免费毛片| 国产探花在线观看一区二区| 神马国产精品三级电影在线观看| 日日摸夜夜添夜夜添小说| 男人舔奶头视频| 一个人免费在线观看电影| 亚洲不卡免费看| 在线观看66精品国产| 免费高清视频大片| 免费在线观看成人毛片| 搞女人的毛片| 一进一出好大好爽视频| 亚洲人与动物交配视频| 又黄又粗又硬又大视频| 久久99热这里只有精品18| 久9热在线精品视频| 在线十欧美十亚洲十日本专区| 亚洲av美国av| 一级黄片播放器| 欧美激情久久久久久爽电影| 五月玫瑰六月丁香| 美女高潮的动态| 又紧又爽又黄一区二区| a级毛片a级免费在线| 桃色一区二区三区在线观看| 亚洲色图av天堂| 亚洲av美国av| av片东京热男人的天堂| 日本成人三级电影网站| 波多野结衣高清作品| 免费在线观看亚洲国产| 亚洲最大成人中文| 亚洲av美国av| 又爽又黄无遮挡网站| 国产毛片a区久久久久| 波多野结衣高清作品| 日本黄色片子视频| 国产精品亚洲一级av第二区| 噜噜噜噜噜久久久久久91| 日韩欧美在线二视频| 亚洲成人久久性| 两个人视频免费观看高清| 精品一区二区三区av网在线观看| 国产野战对白在线观看| 一本久久中文字幕| 好看av亚洲va欧美ⅴa在| 久久99热这里只有精品18| bbb黄色大片| 精品福利观看| 久久久久久国产a免费观看| 精品久久久久久成人av| 午夜福利在线观看免费完整高清在 | 女人十人毛片免费观看3o分钟| 成年免费大片在线观看| 岛国在线观看网站| 丁香欧美五月| 国产精品 国内视频| 99久久精品一区二区三区| av欧美777| bbb黄色大片| 亚洲精品一卡2卡三卡4卡5卡| 欧美一区二区国产精品久久精品| 熟女少妇亚洲综合色aaa.| 国产午夜福利久久久久久| 亚洲av第一区精品v没综合| 99热精品在线国产| 美女大奶头视频| 熟女人妻精品中文字幕| 欧美日韩瑟瑟在线播放| 国产成人福利小说| svipshipincom国产片| 精品电影一区二区在线| 啪啪无遮挡十八禁网站| 少妇的逼水好多| 国产精品久久久人人做人人爽| 中文字幕精品亚洲无线码一区| 在线观看免费视频日本深夜| 最近在线观看免费完整版| xxxwww97欧美| 熟女电影av网| 深爱激情五月婷婷| 观看免费一级毛片| 欧美日本亚洲视频在线播放| 最后的刺客免费高清国语| 韩国av一区二区三区四区| 国产精品免费一区二区三区在线| 操出白浆在线播放| 日本成人三级电影网站| 国产一区二区三区在线臀色熟女| 亚洲第一欧美日韩一区二区三区| 不卡一级毛片| 免费人成视频x8x8入口观看| 国产欧美日韩精品一区二区| 丰满的人妻完整版| xxxwww97欧美| 亚洲中文字幕一区二区三区有码在线看| 午夜亚洲福利在线播放| 日本一本二区三区精品| 午夜a级毛片| 美女高潮的动态| 三级国产精品欧美在线观看| 观看美女的网站| a级一级毛片免费在线观看| 精品国产美女av久久久久小说| 亚洲一区高清亚洲精品| 欧美日韩黄片免| 欧美黑人巨大hd| 亚洲午夜理论影院| 宅男免费午夜| 91麻豆av在线| 日韩精品青青久久久久久| 精品久久久久久久末码| 日本熟妇午夜| 国产精品爽爽va在线观看网站| 久久人妻av系列| 一区福利在线观看| 精品日产1卡2卡| 女人高潮潮喷娇喘18禁视频| 亚洲成av人片免费观看| 一个人看的www免费观看视频| 18禁在线播放成人免费| 99久久精品一区二区三区| 国产欧美日韩一区二区三| 亚洲欧美精品综合久久99| 成人国产一区最新在线观看| 欧美乱妇无乱码| 亚洲专区国产一区二区| 观看美女的网站| 在线播放无遮挡| 国产一区二区在线av高清观看| 亚洲精品乱码久久久v下载方式 | 免费看光身美女| 亚洲精品国产精品久久久不卡| 精品久久久久久久久久久久久| 男女午夜视频在线观看| 在线观看日韩欧美| 久久精品91蜜桃| 99久久无色码亚洲精品果冻| 久久精品国产综合久久久| 久久精品综合一区二区三区| 观看免费一级毛片| 亚洲美女视频黄频| 久久草成人影院| 女人高潮潮喷娇喘18禁视频| 欧美日韩福利视频一区二区| 国产午夜福利久久久久久| 亚洲五月天丁香| 国产精品乱码一区二三区的特点| 最近视频中文字幕2019在线8| 中文字幕高清在线视频| 欧美黑人欧美精品刺激| 亚洲不卡免费看| 午夜福利欧美成人| 欧美一级毛片孕妇| 久久人人精品亚洲av| 热99re8久久精品国产| 国内精品美女久久久久久| 国模一区二区三区四区视频| 欧美精品啪啪一区二区三区| 日本免费a在线| 久久精品国产综合久久久| 亚洲内射少妇av| 窝窝影院91人妻| 亚洲18禁久久av| 久久久久久久亚洲中文字幕 | 国产真人三级小视频在线观看| av片东京热男人的天堂| 日本免费一区二区三区高清不卡| 麻豆成人午夜福利视频| 国产免费一级a男人的天堂| 国产精品久久视频播放| 激情在线观看视频在线高清| 成人午夜高清在线视频| 在线十欧美十亚洲十日本专区| 香蕉丝袜av| netflix在线观看网站| 欧美一区二区国产精品久久精品| 免费人成在线观看视频色| 国产精品一区二区三区四区久久| av黄色大香蕉| 国产欧美日韩一区二区三| 尤物成人国产欧美一区二区三区| 欧美bdsm另类| 亚洲第一欧美日韩一区二区三区| 亚洲美女黄片视频| 很黄的视频免费| 亚洲精品亚洲一区二区| 一级作爱视频免费观看| 国内精品美女久久久久久| 国产精华一区二区三区| 成年女人看的毛片在线观看| 美女cb高潮喷水在线观看| 三级男女做爰猛烈吃奶摸视频| 国产精品 国内视频| 两性午夜刺激爽爽歪歪视频在线观看| 欧美一区二区亚洲| 91在线观看av| 性色av乱码一区二区三区2| 亚洲va日本ⅴa欧美va伊人久久| 亚洲七黄色美女视频| 搞女人的毛片| 一夜夜www| 在线十欧美十亚洲十日本专区| 国产伦人伦偷精品视频| 成人高潮视频无遮挡免费网站| 久久6这里有精品| 国产色婷婷99| 国产国拍精品亚洲av在线观看 | 亚洲成人免费电影在线观看| 一边摸一边抽搐一进一小说| 欧美区成人在线视频| 俺也久久电影网| 免费av不卡在线播放| 一边摸一边抽搐一进一小说| 国产亚洲精品久久久久久毛片| 美女cb高潮喷水在线观看| 成人特级黄色片久久久久久久| 毛片女人毛片| 91久久精品电影网| 亚洲欧美日韩高清专用| 亚洲专区中文字幕在线| av片东京热男人的天堂| 亚洲国产精品999在线| 校园春色视频在线观看| 草草在线视频免费看| 国产aⅴ精品一区二区三区波| 五月玫瑰六月丁香| 国产av麻豆久久久久久久| or卡值多少钱| 一本综合久久免费| 最近最新中文字幕大全免费视频| 欧美不卡视频在线免费观看| 老鸭窝网址在线观看| 日本成人三级电影网站| 黄片大片在线免费观看| 色哟哟哟哟哟哟| 国产真实乱freesex|