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

    城市小型湖泊浮游植物群落結(jié)構(gòu)特征及演替規(guī)律
    ——以武漢沙湖為例

    2017-11-03 06:39:49吳亦瀟張維昊
    生態(tài)學(xué)報(bào) 2017年18期
    關(guān)鍵詞:物種優(yōu)勢

    魯 蕾,吳亦瀟,張維昊,2,*

    1 武漢大學(xué)資源與環(huán)境科學(xué)學(xué)院,武漢 430072 2 水資源安全保障湖北省協(xié)同創(chuàng)新中心,武漢 430072

    城市小型湖泊浮游植物群落結(jié)構(gòu)特征及演替規(guī)律
    ——以武漢沙湖為例

    魯 蕾1,吳亦瀟1,張維昊1,2,*

    1 武漢大學(xué)資源與環(huán)境科學(xué)學(xué)院,武漢 430072 2 水資源安全保障湖北省協(xié)同創(chuàng)新中心,武漢 430072

    城市小型湖泊具有對環(huán)境變化的高敏感性,其浮游植物群落特征和演替規(guī)律對湖泊生態(tài)系統(tǒng)的監(jiān)測與管理有重要意義。本研究于2015年3—12月,對武漢市沙湖進(jìn)行了生態(tài)和水質(zhì)監(jiān)測,共鑒定出藻類10門312種,綠藻、硅藻和藍(lán)藻種類數(shù)占優(yōu)。為探究富營養(yǎng)化城市小型湖泊浮游植物群落時空異質(zhì)性及其演替驅(qū)動力,選取了4個多樣性指數(shù)(Shannon-Wiener指數(shù)、Simpson指數(shù)、Pielou均勻度指數(shù)和Whittaker指數(shù)),并進(jìn)行了多樣性指數(shù)與環(huán)境因子的相關(guān)性分析、主成分分析和聚類分析,計(jì)算了優(yōu)勢種在監(jiān)測時段的3個種間聯(lián)結(jié)性指數(shù)和校正χ2檢驗(yàn)系數(shù)。結(jié)果表明,夏秋季節(jié)浮游植物群落處于演替的穩(wěn)定狀態(tài),環(huán)境因子驅(qū)動作用明顯,其中離子濃度的影響大于溫度的影響,鈉、鎂離子濃度為主要影響因素;藍(lán)藻門和綠藻門優(yōu)勢種聯(lián)結(jié)性較高,硅藻門優(yōu)勢種內(nèi)部聯(lián)結(jié)性較高而與其他門類相關(guān)性較低,從季節(jié)上看,夏秋季節(jié)的優(yōu)勢種種間聯(lián)結(jié)性較高,春冬季優(yōu)勢種與其他時段優(yōu)勢種相關(guān)性低,可能為群落在冬春季不穩(wěn)定的原因。

    城市湖泊; 浮游植物; 優(yōu)勢種; α多樣性; 常量離子; 種間聯(lián)結(jié)性

    近年來,城市湖泊成為各界關(guān)注的熱點(diǎn)之一。其中,城市小型湖泊由于其具有易受城市生態(tài)系統(tǒng)干擾的特性,受到了廣泛關(guān)注。研究者發(fā)現(xiàn)人類通過農(nóng)業(yè)和氣候?qū)丛斐沙掷m(xù)干擾,會使浮游植物群落受到明顯影響[1];同時,研究表明浮游植物群落的初級生產(chǎn)力變化與水體生態(tài)系統(tǒng)的穩(wěn)定性有相互影響的關(guān)系[2, 3]。因此,城市小型湖泊浮游植物群落結(jié)構(gòu)特征和演替規(guī)律的研究對于其生態(tài)保護(hù)有重要的意義。

    在小型湖泊中易發(fā)生多種污染,小型湖泊浮游動植物群落對有機(jī)污染物的抗性與群落結(jié)構(gòu)有明顯的相關(guān)性,均呈現(xiàn)出季節(jié)性變化的特點(diǎn)[4],群落中的物種變化可能比表觀水體性質(zhì)參數(shù)如葉綠素濃度、總磷要更敏感[5]。目前對于氮磷營養(yǎng)鹽與浮游植物群落相關(guān)性的研究較多[6],但是湖泊水體成分較為復(fù)雜,除氮磷外,浮游植物群落結(jié)構(gòu)變化與水體中的其他組分存在一定關(guān)系,Tilman[7]等認(rèn)為常量和微量金屬元素在一些情況下會成為浮游植物生長的限制因素;Holopainen[8]等研究發(fā)現(xiàn)礦地廢水流入造成陽離子濃度顯著升高的淺水湖泊中,浮游植物群落生物量達(dá)到了對照無污染湖泊的兩倍。因此對于水體中的常量離子對浮游植物群落的影響應(yīng)該給予關(guān)注。

    目前,富營養(yǎng)化是城市小型湖泊面臨的主要環(huán)境問題。在富營養(yǎng)化湖泊中,營養(yǎng)鹽不再是限制因素,浮游動植物群落結(jié)構(gòu)和群落內(nèi)種間關(guān)系成為群落變化的原因[9]。優(yōu)勢種和優(yōu)勢類群很大程度上可以成為湖泊水質(zhì)的評判標(biāo)準(zhǔn)[10],并具有一定的時序變化性[11]。目前研究者普遍進(jìn)行的是群落物種表觀數(shù)量關(guān)系研究[12-13],優(yōu)勢種種間關(guān)聯(lián)研究缺乏。而優(yōu)勢種作為群落結(jié)構(gòu)的指示,其種間關(guān)聯(lián)對群落演替具有潛在的驅(qū)動作用,并有較高的研究價值。

    圖1 研究區(qū)采樣點(diǎn)設(shè)置Fig.1 Sample location of area under study

    綜上所述,本研究以沙湖為例,研究了典型富營養(yǎng)化城市小型湖泊浮游植物群落的時空異質(zhì)性,并從環(huán)境因素和優(yōu)勢種更替兩個方面探究浮游植物群落的演替規(guī)律。其中,環(huán)境因子通過參與復(fù)雜的生化反應(yīng)影響群落的演替過程(即群落演替外驅(qū)動力),而優(yōu)勢種更替則直接影響群落的演替方向(即群落演替內(nèi)驅(qū)動力)。

    1 實(shí)驗(yàn)方法

    1.1 采樣點(diǎn)設(shè)置

    采樣區(qū)域?yàn)樯鷳B(tài)功能處于退化狀態(tài)的城市小型湖泊——沙湖,位于武漢市武昌東北部,水域面積3.197km2,于2013年完成一期公園建設(shè),通過人工挖掘的楚河與東湖連通。設(shè)置的采樣點(diǎn)見圖1。

    采樣點(diǎn)坐標(biāo)見表1。采樣點(diǎn)分布的湖區(qū)存在明顯差異,5號采樣點(diǎn)為碼頭;3、4號采樣點(diǎn)位于挺水植物生長區(qū)邊緣;1、2、6號采樣點(diǎn)鄰近公園,其中1號點(diǎn)鄰近公園道路,2號點(diǎn)鄰近公園樹林,6號點(diǎn)鄰近公園草坡。

    1.2 群落樣品采集和水質(zhì)參數(shù)測定

    由于武漢1、2月多雨雪天氣,水溫偏低,表層浮游植物群落采集困難,因此設(shè)計(jì)的浮游植物群落采集時間為2015年3—12月,樣品采集使用1L有機(jī)玻璃采水器,采集表層0.5m處水樣,使用魯哥氏液進(jìn)行固定濃縮后,于高倍顯微鏡下使用0.1mL計(jì)數(shù)框分類計(jì)數(shù)。

    常量陽離子監(jiān)測時間段為2015年5—12月。陽離子濃度測定使用EP- 600便攜式離子色譜儀,陽離子保護(hù)柱:Thermo scientific Dionex IonPacTMCG12A RFICTMGuard 4×50 mm,陽離子交換柱:Thermo scientific Dionex IonPacTMCG12A RFICTMAnalytical 4×250 mm,流動相:18mmol/L甲烷磺酸。

    表1 采樣點(diǎn)位

    1.3 多樣性指數(shù)計(jì)算

    4個多樣性指數(shù),分別為Shannon-Wiener指數(shù)、Simpson指數(shù)、Pielou均勻度指數(shù)[14]和Whittaker指數(shù)[15]。計(jì)算公式如下:

    Shannon-Wiener指數(shù):

    (1)

    Pi=ni/N

    (2)

    Simpson指數(shù):

    (3)

    Pielou均勻度指數(shù):

    J=H/lnS

    (4)

    Whittaker指數(shù):

    β=S/A-1

    (5)

    式中,S為每月采樣出現(xiàn)的物種總數(shù),A為每月各采樣點(diǎn)的平均物種數(shù),Pi為第i種的相對重要值,ni為采樣點(diǎn)第i種物種的個體數(shù),N為采樣點(diǎn)各物種的總個體數(shù)。H值范圍大于等于3為清潔;2—3為輕度污染;1—2為中度污染;小于1為重污染[16]。J值范圍0—1時,J值大體現(xiàn)種間個體分布均勻,污染程度增加J值減小[17]。

    1.4 優(yōu)勢種與種間聯(lián)結(jié)性計(jì)算

    優(yōu)勢種使用優(yōu)勢度進(jìn)行判斷[17],種間聯(lián)結(jié)性選取了3個聯(lián)結(jié)性指數(shù),Ochiai指數(shù)、Dice指數(shù)和Jaccard指數(shù),以及校正χ2檢驗(yàn)系數(shù)[19]進(jìn)行描述。

    優(yōu)勢度:

    Y=fi(ni/N)

    (6)

    Ochiai指數(shù):

    (7)

    Dice指數(shù):

    (8)

    Jaccard指數(shù):

    (9)

    校正χ2檢驗(yàn)系數(shù):

    (10)

    式中,fi為物種i的出現(xiàn)頻率;a為兩種物種均出現(xiàn)的采樣月數(shù),b、c為只有其中一種物種出現(xiàn)的采樣月數(shù),d為兩種物種均未出現(xiàn)的采樣月數(shù),n為采樣總月數(shù)。取Y≥0.02為優(yōu)勢種[17]。當(dāng)ad>bc時為正聯(lián)結(jié),當(dāng)ad3.841時為顯著聯(lián)結(jié),當(dāng)χ2>6.635時為極顯著聯(lián)結(jié)[18]。

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

    數(shù)據(jù)分析使用軟件Matlab2014a,聚類分析使用K-均值法[19]進(jìn)行計(jì)算。數(shù)據(jù)繪圖使用Originlab 8.5。

    2 結(jié)果

    2.1 浮游植物群落時空異質(zhì)性

    2.1.1 群落結(jié)構(gòu)變化規(guī)律

    沙湖浮游植物群落的物種組成和相對豐度均表現(xiàn)出明顯的季節(jié)性。

    從群落物種組成來看(圖2),2015年3—12月,綠藻門(Chlorophyta)和硅藻門(Bacillariopgyta)物種數(shù)超過沙湖的浮游植物群落總體的70%;其中,春季硅藻門物種數(shù)超過總體的30%,夏秋季綠藻門物種數(shù)超過總體的50%;藍(lán)藻門(Cyanophyta)物種數(shù)占比范圍為7%—14%,5—7月以及11—12月藍(lán)藻門物種數(shù)占比超過10%;裸藻門(Euglenophyta)物種數(shù)占比(3%—9%)呈現(xiàn)波形變化,5月達(dá)最小值,8月、12月達(dá)到最大值;黃藻門(Xanthophyta)物種數(shù)占比(3%—6%)4月達(dá)最小值,8月達(dá)最大值,其他各月占比穩(wěn)定。

    由群落相對豐度變化可知(圖3),2015年3—12月綠藻門占比均超過總豐度的50%;藍(lán)藻門(1%—20%)和硅藻門(9%—30%)生物量比例呈波形變化,最大值分別出現(xiàn)在8月和4月,最小值分別出現(xiàn)在3月和6月;裸藻門占比大幅增加出現(xiàn)在4月和9—10月,黃藻門占比大幅增加出現(xiàn)在5—8月;甲藻門(Dinophyta)在3月、金藻門(Crysophyceae)在12月比重明顯增大。隱藻門(Cryptophyta)、定鞭藻門(Haptophyta)和紅藻門(Rhodophyta)物種數(shù)和相對豐度占群落比例均非常小,為稀有種或偶見種。

    圖2 不同采樣月沙湖浮游植物群落物種組成 Fig.2 Taxonomic composition of phytoplankton in different sampling months in Sand Lake

    圖3 不同采樣月沙湖浮游植物群落物種相對豐度 Fig.3 Relative abundance of phytoplankton in different sampling months in Sand Lake

    2.1.2 群落演替規(guī)律

    (1)α和β多樣性時間異質(zhì)性

    β多樣性指數(shù)為Whittaker指數(shù);其值降低說明生境異質(zhì)性增加[17]。由圖4(a)可知,β多樣性分布有明顯季節(jié)性。沙湖各采樣點(diǎn)浮游植物群落的異質(zhì)性在7—9月最高,此時間段,采樣點(diǎn)內(nèi)的物種數(shù)較多而采樣點(diǎn)間的物種差異最大,群落處于物種更新緩慢的穩(wěn)定狀態(tài);3—7月群落異質(zhì)性逐漸增強(qiáng),表明此時段沙湖整體浮游植物群落在向穩(wěn)定狀態(tài)演替;9—12月群落異質(zhì)性逐漸下降,表明此時段浮游植物群落開始向不穩(wěn)定狀態(tài)演替;12月浮游植物群落分布為監(jiān)測時段最均勻的月份(Whittaker指數(shù)達(dá)最大值),說明此時沙湖各采樣點(diǎn)的共有物種明顯減少,整體物種更新速率加快。

    圖4 監(jiān)測時段沙湖浮游植物群落α、β多樣性指數(shù)Fig.4 α and β diversity indexes of phytoplankton in Sand Lake from March to December

    圖5 監(jiān)測時段沙湖浮游植物群落三種α多樣性聚類圖 Fig.5 Cluster result of three kinds of α diversity indexes of phytoplankton in Sand Lake in monitoring period

    α多樣性指數(shù)分別為Shannon-Wiener指數(shù)H、Simpson指數(shù)D和Pielou均勻度指數(shù)J;其中,H值和J值反映水體受污染的狀態(tài),H值與豐富度的關(guān)系最密切,D值和J值對均勻度的敏感性較高[14]。由圖4(b)可知,α多樣性呈現(xiàn)明顯波形變化規(guī)律??傮w來說,監(jiān)測時段內(nèi)H值均集中在1—2之間,J值大都集中在0.3—0.4之間,說明水體保持在中度污染狀態(tài);3、6、12月J值低于0.3說明此時水體處于重污染狀態(tài)。α多樣性分布規(guī)律β多樣性基本一致,7—9月沙湖的平均α多樣性指數(shù)最高,說明此時各采樣點(diǎn)群落物種豐富度和均勻度均達(dá)到穩(wěn)定狀態(tài);3—7月為α多樣性增加段,9—12月為α多樣性下降段,總的來說,3—12月沙湖浮游植物群落的演替規(guī)律為由不穩(wěn)定向穩(wěn)定再向不穩(wěn)定群落發(fā)展的過程。

    (2) 浮游植物群落空間異質(zhì)性

    分別以Shannon-Wiener指數(shù)、Simpson指數(shù)和Pielou均勻度指數(shù)為每月每采樣點(diǎn)的x,y,z坐標(biāo)進(jìn)行聚類分析(圖5),最優(yōu)得到4類,各類中心信息見表2。

    表2 α多樣性指數(shù)聚類中心信息

    由表2可知,4類的綜合α多樣性排序?yàn)椋?>1>2>4。對監(jiān)測時段各采樣點(diǎn)α多樣性根據(jù)聚類進(jìn)行分類可知(表3),沙湖群落存在明顯的時空異質(zhì)性。從空間分布來看,2號點(diǎn)α多樣性最高,群落狀態(tài)最穩(wěn)定;1、3號點(diǎn)以及4、6號點(diǎn)差異不大,其群落演替狀態(tài)較接近;5號點(diǎn)α多樣性最低,群落受到擾動最大;總體來看,沙湖的浮游植物群落空間異質(zhì)性表現(xiàn)為湖泊中段群落較兩側(cè)穩(wěn)定,同一側(cè)湖岸的兩端群落演替狀態(tài)相似,湖泊中人類頻繁活動區(qū)域群落多樣性最低。

    表3 沙湖浮游植物群落α多樣性聚類結(jié)果

    2.2 浮游植物群落演替外驅(qū)動力

    2.2.1 環(huán)境因子監(jiān)測結(jié)果

    本研究考慮的環(huán)境因子包括溫度、鈉、鉀、鎂、鈣的陽離子濃度,除溫度外環(huán)境因子測定時間為5—12月。由圖6可知,水質(zhì)參數(shù)結(jié)果分布呈明顯季節(jié)性,總體上為3—7月逐漸增大,9—12月逐漸減小,8月出現(xiàn)夏季最低值。鈣離子濃度明顯高于其他3種離子,9月各離子濃度達(dá)到最大值,7月溫度達(dá)到最大值,4—10月表層水溫度高于20℃,12月水溫低于10℃。

    圖6 監(jiān)測時段沙湖水質(zhì)因子變化圖Fig.6 The results of water quality factors of Sand Lake from March to December

    2.2.2 外驅(qū)動力分析

    (1)多樣性與環(huán)境因子相關(guān)性分析

    由表4可知,Whittaker指數(shù)與鉀、鈣離子濃度和溫度有極強(qiáng)負(fù)相關(guān)關(guān)系,說明當(dāng)其濃度增大或溫度增加,各采樣點(diǎn)間的異質(zhì)性增強(qiáng)。Shannon-Wiener指數(shù)與鎂離子濃度和溫度有強(qiáng)相關(guān)關(guān)系,與鈉離子濃度有極強(qiáng)相關(guān)關(guān)系;Simpson指數(shù)與鈉離子濃度有強(qiáng)相關(guān)關(guān)系,與鎂離子濃度有極強(qiáng)相關(guān)關(guān)系;Pielou均勻度指數(shù)與鈉、鎂離子濃度有強(qiáng)相關(guān)關(guān)系。相關(guān)性分析結(jié)果說明,環(huán)境因子對浮游植物群落狀態(tài)有明顯的影響,其中常量離子濃度對群落演替有比溫度更強(qiáng)的影響力,鈉、鎂離子濃度是α多樣性的主要影響因素,溫度主要影響群落物種豐富度,鉀、鈣離子濃度主要影響沙湖浮游植物群落β多樣性。

    表4 多樣性指數(shù)與水質(zhì)因子相關(guān)性分析結(jié)果

    圖7 監(jiān)測時段沙湖生物多樣性與水質(zhì)因子主成分分析聚類圖 Fig.7 Cluster results of principal component analysis of biodiversity indexes and water quality factors from March to December主成分3對應(yīng)的數(shù)值為原數(shù)值減去最小值,以使點(diǎn)均落在z軸正向; 藍(lán)色點(diǎn)為類別1,綠色點(diǎn)為類別2,紅色點(diǎn)為類別3

    (2) 主成分分析

    主成分分析考慮的因素包括α多樣性指數(shù)與環(huán)境因子。得到的主成分中,主成分1、主成分2和主成分3分別能夠解釋60.01%、 24.44%和11.28%(累積95.73%)的數(shù)據(jù)量(表5)。根據(jù)監(jiān)測數(shù)據(jù)對應(yīng)3個主成分空間的變換值進(jìn)行聚類,最優(yōu)可得3類結(jié)果,對3類結(jié)果進(jìn)行整理見表6,主成分分析聚類結(jié)果見圖7。

    主成分分析是從原本存在相關(guān)性的各成分中提取出同樣個數(shù)的完全線性無關(guān)的主成分進(jìn)行分析的過程。由表5可知,主成分1綜合性較強(qiáng);主成分2多樣性指數(shù)為主要因子;主成分3鎂離子濃度、溫度為主要因子。

    由表6類別中心可知,類別3主成分1為主要因素,該類點(diǎn)α多樣性和環(huán)境因子數(shù)值均較高;類別2處于3個成分的負(fù)象限,說明該類點(diǎn)α多樣性和環(huán)境因子數(shù)值均較低;類別1主成分1和2為主要因素,主成分1為負(fù)象限,說明該類點(diǎn)的α多樣性較高而離子濃度等環(huán)境因子數(shù)值較低。由表6分可知,對各采樣點(diǎn)位來說,6月均為2類,7—10月均為3類,此時段內(nèi)群落狀態(tài)明顯受到環(huán)境因子的驅(qū)動,表現(xiàn)為兩者分布趨勢一致;5、11、12月大多數(shù)采樣點(diǎn)為1類,說明此時段內(nèi),除了環(huán)境因子的影響外,可能存在其他作用較強(qiáng)的影響因素。結(jié)合表6和表3的聚類結(jié)果可知,夏秋季節(jié)群落狀態(tài)較穩(wěn)定,環(huán)境因子有明顯的驅(qū)動性,冬春季節(jié)群落處于演替狀態(tài),環(huán)境因子的驅(qū)動力減弱,內(nèi)驅(qū)動力可能為主要作用因素。

    表5 主成分系數(shù)

    表6 主成分分析聚類結(jié)果

    每一類中心為:類別1(-1.642,1.154,0.227),類別2(-2.384,-1.847,-0.678),類別3(1.889,-0.154,0.0745)

    2.3 浮游植物群落演替內(nèi)驅(qū)動力

    2.3.1 優(yōu)勢種變化規(guī)律

    每月的優(yōu)勢種由優(yōu)勢度確定,3—12月出現(xiàn)的優(yōu)勢種共40種(表7)。其中綠藻門21種,硅藻門10種,藍(lán)藻門4種,黃藻門和裸藻門各2種,金藻門和甲藻門各1種。由表7可知,每月的優(yōu)勢種數(shù)量呈現(xiàn)先增多再減少的規(guī)律,最大值為9月,達(dá)到17種,最小為3月,有5種;其他各月基本上呈現(xiàn)出季節(jié)上的穩(wěn)定性;監(jiān)測時段主要優(yōu)勢種有5個,為S9、S10、S17、S23、S33,主要優(yōu)勢種集中在春季和夏季,且均為綠藻。

    表7 監(jiān)測時段沙湖各月優(yōu)勢種匯總信息

    *:主要優(yōu)勢種

    2.3.2 內(nèi)驅(qū)動力分析

    優(yōu)勢種每月分布的聯(lián)結(jié)性計(jì)算結(jié)果見表8。Ochiai、Dice、Jaccard指數(shù)確定的完全相關(guān)的物種分別為S4和S3、S20和S19、S21和S19、S21和S20、S25和S24、S30和S29;這6對物種同時作為優(yōu)勢種出現(xiàn)的時間為9月,且其優(yōu)勢度均在0.02-0.04中間,不是主要優(yōu)勢種。Jaccard指數(shù)強(qiáng)相關(guān)的物種分別為S9和S10,S32和S31;Dice指數(shù)強(qiáng)相關(guān)的物種有5對,分別為S9和S10,S32和S31,以及S27和S28、S33和S31、S39和S40;Ochiai指數(shù)強(qiáng)相關(guān)的物種有8對,分別為S9和S10,S32和S31、S27和S28、S33和S31、S39和S40,以及S33和S32、S11和S9、S29和S27。

    表8 監(jiān)測時段優(yōu)勢種各聯(lián)結(jié)性指數(shù)半矩陣分析結(jié)果

    校正χ2檢驗(yàn)結(jié)果見表9。χ2檢驗(yàn)結(jié)果中顯著正聯(lián)結(jié)的物種有6對,分別為S22和S19、S28和S27、S29和S27、S30和S27、S33和S31、S40和S39;極顯著正聯(lián)結(jié)的物種有7對,分別為S4和S3、S20和S19、S21和S19、S21和S20、S25和S24、S30和S29以及S32和S31;不相關(guān)的結(jié)果中,S9的與編號在其后的優(yōu)勢種不相關(guān),2對不相關(guān)的物種分別為S26和S22、S39和S35。

    表9 監(jiān)測時段優(yōu)勢種校正χ2檢驗(yàn)半矩陣分析結(jié)果

    種間聯(lián)結(jié)性測定值本身有一定的生態(tài)學(xué)意義,在一定程度上衡量了物種的相互關(guān)系和物種對環(huán)境綜合因子反應(yīng)的差異[18];χ2檢驗(yàn)?zāi)鼙容^準(zhǔn)確度量種對的關(guān)聯(lián)顯著程度[20]。3個種間聯(lián)結(jié)性指數(shù)結(jié)果略有不同,Jaccard、Dice指數(shù)得到的信息量小于Ochiai指數(shù)。優(yōu)勢種種間聯(lián)結(jié)性表現(xiàn)為同屬相關(guān)性較強(qiáng),包括三個指數(shù)均為強(qiáng)相關(guān)的2種小球藻和2種平裂藻,以及Ochiai和Dice指數(shù)強(qiáng)相關(guān)的直鏈藻,不同屬也存在強(qiáng)相關(guān)的現(xiàn)象,包括三角四角藻和普通黃絲藻、綠球藻和橢圓小球藻;綠藻門和藍(lán)藻門優(yōu)勢種顯著相關(guān),包括集球藻和2種平裂藻;硅藻門和其他藻類相關(guān)性均較低,表現(xiàn)為同屬的高相關(guān)性。χ2檢驗(yàn)結(jié)果和聯(lián)結(jié)性指數(shù)計(jì)算結(jié)果基本一致,但χ2檢驗(yàn)結(jié)果出現(xiàn)不相關(guān)的物種,包括模糊直鏈藻和綠脆桿藻、三角四角藻和月牙藻,其中,χ2檢驗(yàn)中橢圓小球藻與其他物種相關(guān)性均較低。

    總體來說,9月群落處于監(jiān)測時段最穩(wěn)定狀態(tài),4月和9月優(yōu)勢種相關(guān)性較高,夏季和秋季優(yōu)勢種相關(guān)性較高,該結(jié)果與表6和表3結(jié)果一致,說明群落物種結(jié)構(gòu)的強(qiáng)關(guān)聯(lián)性對夏秋季節(jié)沙湖浮游植物群落演替的穩(wěn)定狀態(tài)有一定貢獻(xiàn);5和11月優(yōu)勢種僅2種(S31,S32)與7、8月優(yōu)勢種S33有強(qiáng)相關(guān)性,其他優(yōu)勢種表現(xiàn)為月內(nèi)優(yōu)勢種極強(qiáng)相關(guān)性(S9和S10)和不與其他時段優(yōu)勢種相關(guān)的情況,12月優(yōu)勢種和其他優(yōu)勢種相關(guān)性也普遍較低,出現(xiàn)的是耐低溫的物種;該結(jié)果與表6結(jié)果一致,說明5、11、12月浮游植物群落演替主要由群落內(nèi)物種構(gòu)成變化驅(qū)動。

    3 討論和結(jié)論

    3.1 城市小型湖泊浮游植物群落時空分布異質(zhì)性

    沙湖的浮游植物群落結(jié)構(gòu)與分布存在明顯的時空變化規(guī)律。沙湖是典型的富營養(yǎng)化湖泊,主要優(yōu)勢類群為綠藻-硅藻-藍(lán)藻(累積物種數(shù)占比超過80%),研究者普遍認(rèn)為優(yōu)勢類群的占比可以用于指示湖泊的污染狀態(tài)[10],硅藻的物種數(shù)占比從春季(約30%)到夏季(低于20%)有明顯下降,說明沙湖受污染程度從夏季開始加深;藍(lán)藻、裸藻物種數(shù)占比在夏秋季達(dá)到最大值,表明沙湖在夏秋季處于富營養(yǎng)化加劇的狀態(tài)。整個監(jiān)測時段中,7—9月浮游植物群落整體α多樣性最高,β多樣性最低(圖4),群落處于演替的穩(wěn)定階段,該時段的群落表現(xiàn)為耐污種占主要(表7),各采樣點(diǎn)在該時段群落多樣性差異不大。α多樣性較低的3—7月和9—12月,多樣性分布存在明顯空間差異(表3),5號點(diǎn)碼頭的多樣性總體來說最小,表明人為擾動對小型湖泊生物多樣性有明顯的影響;劉春光[21]等認(rèn)為水生植物能夠通過多種機(jī)制影響湖泊生態(tài)系統(tǒng),3、4號點(diǎn)為人工種植挺水、沉水植物的水域邊緣帶,由采樣點(diǎn)間對比發(fā)現(xiàn),挺水、沉水植物對小型湖泊浮游植物群落的影響較??;6月各點(diǎn)生物多樣性均明顯下降,而1號點(diǎn)多樣性下降最明顯,由1、2、6號采樣點(diǎn)對比結(jié)果可知,小型湖泊的浮游植物群落的穩(wěn)定性與湖濱帶的環(huán)境有相關(guān)關(guān)系。

    3.2 浮游植物群落演替外驅(qū)動力

    外驅(qū)動力主要作用因子為鈉、鎂離子濃度,溫度為輔助作用因子,外驅(qū)動力主要作用時段為夏秋季節(jié)。由主成分分析結(jié)果可知,6—10月α多樣性與環(huán)境因子變化規(guī)律一致,該時段內(nèi)浮游植物群落主要驅(qū)動力為外驅(qū)動力,可能原因?yàn)槎鄻有暂^高的7—9月浮游植物群落處于較穩(wěn)定狀態(tài),物種替換減緩,以物種更新為主的內(nèi)驅(qū)動力減小,環(huán)境因子的驅(qū)動作用相對增大;多樣性較低的6月浮游植物群落處于較敏感狀態(tài),加之水體整體營養(yǎng)水平在該時段較低[22],環(huán)境因子成為群落更新生長的限制因素。通過浮游植物群落多樣性和環(huán)境因子的相關(guān)性分析可知,鈉、鎂離子濃度與群落α多樣性有強(qiáng)相關(guān)性,主要原因可能為鎂離子是藻類的葉綠素的組成元素之一且與葉綠體基粒堆垛直接相關(guān),能夠?qū)θ~綠體熒光產(chǎn)率和光合作用原初過程造成影響[23],該結(jié)果與其他研究結(jié)果類似[24]。本研究發(fā)現(xiàn)溫度與浮游藻類群落的Shannon-Wiener指數(shù)有強(qiáng)相關(guān)性,但與Simpson指數(shù)相關(guān)性較低,與Pielou均勻度指數(shù)基本無相關(guān)性,與大中型湖泊水庫研究結(jié)果[25,26]有明顯差異,表明大中型湖泊與小型湖泊浮游植物群落演替外驅(qū)動力因素存在差異,溫度影響降低可能由于小型湖泊層間溫度和湖區(qū)溫度差異化不明顯。Whittaker指數(shù)與Shannon-Wiener指數(shù)和其他環(huán)境因子均有強(qiáng)負(fù)相關(guān)關(guān)系,說明較高的生境內(nèi)多樣性與較高的常量離子濃度有利于浮游植物群落維持演替的穩(wěn)定狀態(tài)。

    3.3 浮游植物群落演替內(nèi)驅(qū)動力

    浮游植物優(yōu)勢種種類及其數(shù)量對浮游植物群落結(jié)構(gòu)的穩(wěn)定性有重要影響,優(yōu)勢種種類數(shù)越多且優(yōu)勢度越小,則群落結(jié)構(gòu)越復(fù)雜、穩(wěn)定[26],沙湖每月優(yōu)勢種數(shù)量變化表明其群落演替規(guī)律為由不穩(wěn)定到穩(wěn)定再逐步向不穩(wěn)定變化的過程。從優(yōu)勢種分類來看,硅藻和藍(lán)綠藻對環(huán)境因素的響應(yīng)存在明顯差異[27]。藍(lán)藻大量生長的時間段為夏季和冬季,既有相對生物量的增加,也有種類數(shù)的明顯增加;夏季爆發(fā)可能由于該季節(jié)水體狀態(tài)更適宜藍(lán)藻生長[28],冬季爆發(fā)可能由于2015年11月降雨量增大,出現(xiàn)適宜藍(lán)藻生長的面源輸入[29]。硅藻大量生長的時間為春季和秋季,Sarmento[30]等通過對于淺水湖泊的研究發(fā)現(xiàn),外源輸入減少水位降低時,硅藻生物量得到明顯增加;硅藻夏季種類數(shù)比重變化穩(wěn)定,但是生物量比重明顯下降,硅藻與藍(lán)藻占總密度百分比為極顯著負(fù)相關(guān)[31]。硅藻和藍(lán)藻優(yōu)勢種的負(fù)關(guān)聯(lián)性可能是5、11月浮游植物群落結(jié)構(gòu)演替與夏秋季不同的主要驅(qū)動因素。甲藻的生物量最大值出現(xiàn)在3月,達(dá)到6%,胡韌[32]等研究發(fā)現(xiàn),在溫度不是主要的限制因素時,春季甲藻的生物量較高,與硅藻門和綠藻門構(gòu)成優(yōu)勢類群;12月湖水溫度低于10℃,優(yōu)勢種以金藻門為主,其特點(diǎn)為耐低溫性,其與其他優(yōu)勢種幾乎無種間關(guān)聯(lián),可能為冬季群落演替的主要驅(qū)動力。

    [1] Lemma B. Ecological changes in two ethiopian lakes caused by contrasting human intervention. Limnologica-Ecology and Management of Inland Waters, 2003, 33(1): 44- 53.

    [2] Mukankomeje R, Plisnier P D, Descy J P, Massaut L. Lake muzahi, rwanda: limnological features and phytoplankton production. Hydrobiologia, 1993, 257(2): 107- 120.

    [3] Sommer U, Sommer F. Cladocerans versus copepods: The cause of contrasting top-down controls on freshwater and marine phytoplankton. Oecologia, 2006, 147(2): 183- 194.

    [4] Willis K J, Van den Brink P J, Green J D. Seasonal variation in plankton community responses of mesocosms dosed with pentachlorophenol. Ecotoxicology, 2004, 13(7): 707- 720.

    [5] Huszar V L M, Silva L H S, Domingos P, Marinho M, Melo S. Phytoplankton species composition is more sensitive than OECD criteria to the trophic status of three Brazilian tropical lakes. Hydrobiologia, 1998, 369: 59- 71.

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

    [7] Tilman D, Kilham S S, Kilham P. Phytoplankton community ecology:the role of limiting nutrients. Annual Review of Ecology and Systematics, 1982, 13(1): 349- 372.

    [8] Holopainen I J, Holopainen A L, H?m?l?inen H, Rahkola-Sorsa M, Tkatcheva V, Viljanen M. Effects of mining industry waste waters on a shallow lake ecosystem in karelia, north-west russia. Hydrobiologia, 2003, 506(1/3): 111- 119.

    [9] Mataloni G, Tesolín G, Sacullo F, Tell G. Factors regulating summer phytoplankton in a highly eutrophic antarctic lake. Hydrobiologia, 2000, 432(1/3): 65- 72.

    [10] 龐科, 姚錦仙, 王昊, 劉松濤, 李翀, 呂植. 額爾古納河流域秋季浮游植物群落結(jié)構(gòu)特征. 生態(tài)學(xué)報(bào), 2011, 31(12): 3391- 3398.

    [11] 田永強(qiáng), 俞超超, 王磊, 黃邦欽. 福建九龍江北溪浮游植物群落分布特征及其影響因子. 應(yīng)用生態(tài)學(xué)報(bào), 2012, 23(9): 2559- 2565.

    [12] 王麗卿, 施榮, 季高華, 范志鋒, 程婧蕾. 淀山湖浮游植物群落特征及其演替規(guī)律. 生物多樣性, 2011, 19(1): 48- 56.

    [13] 王靜雅, 汪志聰, 李翀, 葉少文, 連玉喜, 劉家壽, 張?zhí)昧? 李鐘杰. 三峽水庫壩前水域浮游植物群落時空動態(tài)研究. 水生態(tài)學(xué)雜志, 2015, 39(5): 877- 884.

    [14] 許晴, 張放, 許中旗, 賈彥龍, 尤建民. Simpson指數(shù)和Shannon-Wiener指數(shù)若干特征的分析及“稀釋效應(yīng)”. 草業(yè)科學(xué), 2011, 28(4): 527- 531.

    [15] 周洪華, 陳亞寧, 李衛(wèi)紅. 塔里木河下游綠洲——荒漠過渡帶植物多樣性特征及優(yōu)勢種群分布格局. 中國沙漠, 2009, 29(4): 688- 696.

    [16] 劉文盈, 張秋良, 高潤宏, 鐵牛, 崔琥. 鹽沼濕地浮游植物多樣性與環(huán)境關(guān)系研究. 內(nèi)蒙古林業(yè)科技, 2012, 38(1): 5- 7, 22- 22.

    [17] 郝媛媛, 孫國鈞, 張立勛, 龔雪平, 許莎莎, 劉慧明, 張芬. 黑河流域浮游植物群落特征與環(huán)境因子的關(guān)系. 湖泊科學(xué), 2014, 26(1): 121- 130.

    [18] 徐敏, 高翔, 呂光輝, 傅德平, 楊曉東. 克拉瑪依典型植物群落種間關(guān)聯(lián)分析. 新疆農(nóng)業(yè)科學(xué), 2010, 47(5): 1041- 1046.

    [19] 周世兵, 徐振源, 唐旭清. 新的k-均值算法最佳聚類數(shù)確定方法. 計(jì)算機(jī)工程與應(yīng)用, 2010, 46(16): 27- 31.

    [20] 陳家長, 孟順龍, 尤洋, 胡庚東, 瞿建宏, 吳偉, 范立民, 馬曉燕. 太湖五里湖浮游植物群落結(jié)構(gòu)特征分析. 生態(tài)環(huán)境學(xué)報(bào), 2009, 18(4): 1358- 1367.

    [21] 劉春光, 邱金泉, 王雯, 莊源益. 富營養(yǎng)化湖泊治理中的生物操縱理論. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2004, 23(1): 198- 201.

    [22] Li H M, Tang H J, Shi X Y, Zhang C S, Wang X L. Increased nutrient loads from the Changjiang (Yangtze) River have led to increased Harmful Algal Blooms. Harmful Algae, 2014, 39: 92- 101.

    [23] 孫谷疇. 鎂離子對柱孢魚腥藻(Anabaena cylindrica)類囊體膜吸收光譜和熒光光譜的影響. 植物生理學(xué)報(bào), 1983, 9(4): 339- 345.

    [24] 劉文盈, 鐵牛, 張秋良, 高潤宏. 藻類在鹽沼濕地分布的關(guān)鍵限制因子研究. 內(nèi)蒙古林業(yè)科技, 2012, 38(4): 16- 18, 26- 26.

    [25] 鄧文麗, 劉均平, 王曉星, 杜桂森, 洪劍明. 北京野鴨湖浮游植物群落結(jié)構(gòu)與水質(zhì)關(guān)系研究. 濕地科學(xué), 2013, 11(1): 27- 34.

    [26] 裴國霞, 郭琦, 魏樂, 李汗青. 金海水庫浮游植物調(diào)查及群落結(jié)構(gòu)研究. 水生態(tài)學(xué)雜志, 2014, 35(2): 68- 73.

    [27] 王宇飛, 趙秀蘭, 何丙輝, 黃琪. 漢豐湖夏季浮游植物群落與環(huán)境因子的典范對應(yīng)分析. 環(huán)境科學(xué), 2015, 36(3): 922- 927.

    [28] 鄭吳柯, 劉憲斌, 趙興貴, 陳曦. 于橋水庫浮游植物群落特征. 中國環(huán)境監(jiān)測, 2015, 31(1): 35- 40.

    [29] Rakko A, Laugaste R, Ott I. Algal blooms in estonian small lakes // Evangelista V, Barsanti L, Frassanito A M, Passarelli V, Gualtieri P, eds. Algal Toxins: Nature, Occurrence, Effect and Detection. Netherlands: Springer, 2008: 211- 220.

    [30] Sarmento H, Isumbisho M, Descy J P. Phytoplankton ecology of Lake Kivu (Eastern Africa). Journal of Plankton Research, 2006, 28(9): 815- 829.

    [31] 李娣, 李旭文, 牛志春, 王霞, 師偉, 于紅霞. 太湖浮游植物群落結(jié)構(gòu)及其與水質(zhì)指標(biāo)間的關(guān)系. 生態(tài)環(huán)境學(xué)報(bào), 2014, 23(11): 1814- 1820.

    [32] 胡韌, 雷臘梅, 韓博平. 南亞熱帶大型貧營養(yǎng)水庫浮游植物群落結(jié)構(gòu)與季節(jié)變化——以新豐江水庫為例. 生態(tài)學(xué)報(bào), 2008, 28(10): 4652- 4664.

    Distributionandsuccessionpatternofphytoplanktoncommunitiesinasmallurbanlake,SandLakeinWuhanCity,China

    LU Lei1, WU Yixiao1, ZHANG Weihao1,2,*

    1SchoolofResourceandEnvironmentalSciences,WuhanUniversity,Wuhan430072,China2HubeiCollaborativeInnovationCenterforWaterResourceSecurity,Wuhan430072,China

    Small urban lakes are highly sensitive to changes in their environment and are closely associated with the surrounding urban systems. The distribution and succession pattern of phytoplankton can significantly influence the entire ecosystem in a small urban lake. From March to December 2015, phytoplankton and water quality were monitored in Sand Lake, Wuhan city, China. In total, 312 species were detected representing 10 phyla, with Chlorophyta, Bacillariophyta, and Cyanophyta being the major phyla represented. Biodiversity indices were calculated using the methods of Shannon-Wiener, Simpson, Pielou, and Whittaker. The biodiversity indices and environmental factors were subjected to correlation analysis, principal component analysis, and cluster analysis. Three interspecific association indices and the calibrated coefficients ofχ2test were calculated. The results showed that the phytoplankton reached a stable state in summer and autumn, during which the environmental factors were the main influence; the effect of constant ion concentration was greater than that of temperature, particularly those of sodium and magnesium. Analysis of the data for the dominant species showed that Chlorophyta and Cyanophyta were significantly correlated, whereas Bacillariophyta showed no significant correlation with other phyla. A significant correlation occurred between summer and autumn dominant species, while dominant species in winter and spring showed no significant correlation with those in other seasons, which could be attributed to the low stability of the phytoplankton during this period.

    urban lake; phytoplankton; dominant species; α diversity; constant ion; inter-specific association

    國家自然科學(xué)基金(41331174);江西省水利廳重大科技項(xiàng)目(KT201412);水利部科技推廣項(xiàng)目(TG1520)

    2016- 06- 21; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 04- 25

    *通訊作者Corresponding author.E-mail: zhangwh@whu.edu.cn

    10.5846/stxb201606211215

    魯蕾,吳亦瀟,張維昊.城市小型湖泊浮游植物群落結(jié)構(gòu)特征及演替規(guī)律——以武漢沙湖為例.生態(tài)學(xué)報(bào),2017,37(18):5993- 6004.

    Lu L, Wu Y X, Zhang W H.Distribution and succession pattern of phytoplankton communities in a small urban lake, Sand Lake in Wuhan City, China.Acta Ecologica Sinica,2017,37(18):5993- 6004.

    猜你喜歡
    物種優(yōu)勢
    物種大偵探
    物種大偵探
    優(yōu)勢 等
    吃光入侵物種真的是解決之道嗎?
    英語世界(2023年10期)2023-11-17 09:18:18
    矮的優(yōu)勢
    趣味(語文)(2020年3期)2020-07-27 01:42:46
    回首2018,這些新物種值得關(guān)注
    電咖再造新物種
    汽車觀察(2018年10期)2018-11-06 07:05:26
    畫與話
    發(fā)揚(yáng)優(yōu)勢 有所作為
    瘋狂的外來入侵物種
    十八禁网站免费在线| 天天操日日干夜夜撸| 国产精品成人在线| 高潮久久久久久久久久久不卡| 国产精品国产高清国产av | 免费人成视频x8x8入口观看| 久久久国产一区二区| 午夜日韩欧美国产| 国产精品影院久久| 久久精品人人爽人人爽视色| 免费在线观看亚洲国产| 9191精品国产免费久久| 新久久久久国产一级毛片| 亚洲av成人av| 精品熟女少妇八av免费久了| 精品一区二区三区av网在线观看| 婷婷成人精品国产| 宅男免费午夜| 两性夫妻黄色片| 伦理电影免费视频| 午夜两性在线视频| 久热这里只有精品99| 亚洲av第一区精品v没综合| 免费在线观看亚洲国产| 免费在线观看视频国产中文字幕亚洲| av在线播放免费不卡| 欧美日韩乱码在线| 成熟少妇高潮喷水视频| 久热爱精品视频在线9| 午夜影院日韩av| 亚洲人成77777在线视频| 后天国语完整版免费观看| 搡老乐熟女国产| 9热在线视频观看99| 亚洲国产毛片av蜜桃av| 高清黄色对白视频在线免费看| 亚洲一区高清亚洲精品| 热99国产精品久久久久久7| 午夜福利视频在线观看免费| 天天添夜夜摸| 可以免费在线观看a视频的电影网站| 搡老熟女国产l中国老女人| 韩国精品一区二区三区| 色94色欧美一区二区| 在线观看日韩欧美| 国产1区2区3区精品| 国产高清视频在线播放一区| 亚洲熟妇熟女久久| 日韩有码中文字幕| 韩国av一区二区三区四区| 久久人人97超碰香蕉20202| 啦啦啦 在线观看视频| 一区在线观看完整版| 精品久久久久久久久久免费视频 | 18禁观看日本| 99在线人妻在线中文字幕 | 色在线成人网| 高清黄色对白视频在线免费看| 老司机福利观看| 久久精品亚洲熟妇少妇任你| 另类亚洲欧美激情| 亚洲一区二区三区欧美精品| svipshipincom国产片| 老熟妇乱子伦视频在线观看| 日韩 欧美 亚洲 中文字幕| 69av精品久久久久久| 亚洲综合色网址| x7x7x7水蜜桃| 香蕉久久夜色| 国产一区在线观看成人免费| 久久精品成人免费网站| 亚洲欧美日韩高清在线视频| 久久久久久久国产电影| 极品教师在线免费播放| 老司机午夜福利在线观看视频| 一本大道久久a久久精品| 男女午夜视频在线观看| 欧美日韩亚洲综合一区二区三区_| e午夜精品久久久久久久| 9191精品国产免费久久| 亚洲精品久久成人aⅴ小说| 精品福利永久在线观看| 丰满迷人的少妇在线观看| 国产无遮挡羞羞视频在线观看| 高清视频免费观看一区二区| 一区在线观看完整版| 啦啦啦视频在线资源免费观看| 嫩草影视91久久| 欧美日本中文国产一区发布| 亚洲五月天丁香| 女人精品久久久久毛片| 久久久国产一区二区| 亚洲中文日韩欧美视频| 亚洲欧美精品综合一区二区三区| 国产高清视频在线播放一区| 国产精品免费一区二区三区在线 | 亚洲免费av在线视频| 午夜成年电影在线免费观看| 国产亚洲一区二区精品| 99国产精品免费福利视频| 一区在线观看完整版| 男男h啪啪无遮挡| 黑人巨大精品欧美一区二区蜜桃| 满18在线观看网站| 国产成人av激情在线播放| 老汉色∧v一级毛片| 久久精品人人爽人人爽视色| 免费久久久久久久精品成人欧美视频| 精品国产亚洲在线| 亚洲精品美女久久av网站| 亚洲人成电影免费在线| 欧美人与性动交α欧美软件| 久久中文看片网| 如日韩欧美国产精品一区二区三区| 人妻丰满熟妇av一区二区三区 | 国产成人影院久久av| 国产精品一区二区免费欧美| 精品人妻1区二区| 免费观看人在逋| 自拍欧美九色日韩亚洲蝌蚪91| 9热在线视频观看99| 久久中文看片网| 欧美日韩福利视频一区二区| 亚洲va日本ⅴa欧美va伊人久久| 国产不卡一卡二| 亚洲一区二区三区欧美精品| 欧美最黄视频在线播放免费 | 久久精品aⅴ一区二区三区四区| 免费看十八禁软件| 看免费av毛片| 在线观看一区二区三区激情| 成年人黄色毛片网站| 成人精品一区二区免费| 亚洲欧美色中文字幕在线| 精品福利永久在线观看| 国精品久久久久久国模美| 国产亚洲精品第一综合不卡| 国产黄色免费在线视频| 在线观看舔阴道视频| 一级a爱片免费观看的视频| av网站免费在线观看视频| 久久热在线av| 老熟女久久久| 麻豆av在线久日| 99精品在免费线老司机午夜| 亚洲美女黄片视频| 一区二区三区精品91| 真人做人爱边吃奶动态| 国产又爽黄色视频| 亚洲色图 男人天堂 中文字幕| 真人做人爱边吃奶动态| 久久精品熟女亚洲av麻豆精品| 亚洲aⅴ乱码一区二区在线播放 | 国产成人啪精品午夜网站| 亚洲av成人一区二区三| 久久精品aⅴ一区二区三区四区| 久久精品成人免费网站| 久久久国产一区二区| 精品国内亚洲2022精品成人 | 在线视频色国产色| 99re6热这里在线精品视频| 日韩视频一区二区在线观看| 久久精品国产亚洲av高清一级| 一级,二级,三级黄色视频| 久久九九热精品免费| 很黄的视频免费| 国产高清国产精品国产三级| 成人精品一区二区免费| 亚洲精品中文字幕在线视频| 亚洲第一av免费看| 免费在线观看完整版高清| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人精品无人区| 中亚洲国语对白在线视频| 人人妻人人澡人人看| 日日摸夜夜添夜夜添小说| 午夜两性在线视频| 一二三四在线观看免费中文在| 九色亚洲精品在线播放| 亚洲国产毛片av蜜桃av| 老司机影院毛片| 午夜福利在线观看吧| 99热网站在线观看| 亚洲精品久久午夜乱码| 国产亚洲一区二区精品| 99热只有精品国产| 一区二区三区精品91| 欧美大码av| 又紧又爽又黄一区二区| 一a级毛片在线观看| 国产伦人伦偷精品视频| 精品国产一区二区久久| 麻豆乱淫一区二区| 久久国产精品男人的天堂亚洲| 欧美激情极品国产一区二区三区| 亚洲国产欧美一区二区综合| www.熟女人妻精品国产| 欧美日韩黄片免| 欧美午夜高清在线| 国产精品久久视频播放| 久久久久久亚洲精品国产蜜桃av| 少妇猛男粗大的猛烈进出视频| 丝袜在线中文字幕| 黄片小视频在线播放| 电影成人av| 亚洲精品国产一区二区精华液| 在线十欧美十亚洲十日本专区| 日韩成人在线观看一区二区三区| 亚洲性夜色夜夜综合| 黄色片一级片一级黄色片| 脱女人内裤的视频| 成人特级黄色片久久久久久久| 欧美午夜高清在线| 国产成人影院久久av| 色尼玛亚洲综合影院| 国产乱人伦免费视频| 免费少妇av软件| 日韩欧美免费精品| 十分钟在线观看高清视频www| 少妇猛男粗大的猛烈进出视频| 国产精品综合久久久久久久免费 | 国产熟女午夜一区二区三区| 久久人妻av系列| 亚洲欧美激情在线| 欧美黑人欧美精品刺激| 在线观看日韩欧美| 亚洲欧美一区二区三区久久| 无遮挡黄片免费观看| 日韩视频一区二区在线观看| 免费在线观看亚洲国产| 色综合婷婷激情| 欧美日韩黄片免| 国产一区二区三区综合在线观看| 亚洲熟女精品中文字幕| svipshipincom国产片| 欧美日韩精品网址| 制服人妻中文乱码| 免费在线观看视频国产中文字幕亚洲| 精品亚洲成国产av| 日本精品一区二区三区蜜桃| 欧美在线一区亚洲| 久久久久久亚洲精品国产蜜桃av| 天天躁日日躁夜夜躁夜夜| 国产野战对白在线观看| 黑丝袜美女国产一区| 免费在线观看亚洲国产| 黄色丝袜av网址大全| 欧美精品啪啪一区二区三区| 啦啦啦免费观看视频1| 国产成人精品无人区| 午夜久久久在线观看| 麻豆成人av在线观看| 一进一出好大好爽视频| 国产精品久久视频播放| 久久影院123| 亚洲三区欧美一区| 99精品欧美一区二区三区四区| 制服诱惑二区| 一级毛片精品| 999久久久国产精品视频| 日韩视频一区二区在线观看| 99热网站在线观看| 中文字幕高清在线视频| 日韩成人在线观看一区二区三区| 国产成人精品久久二区二区91| 国产成人啪精品午夜网站| 国产视频一区二区在线看| 久久中文看片网| 久久香蕉激情| 午夜91福利影院| 美女午夜性视频免费| 90打野战视频偷拍视频| 亚洲欧美激情综合另类| 欧美 亚洲 国产 日韩一| 久久狼人影院| 麻豆av在线久日| 黄片播放在线免费| 每晚都被弄得嗷嗷叫到高潮| 久久久久久久午夜电影 | 最新在线观看一区二区三区| 国产精品久久久久成人av| 美女午夜性视频免费| 亚洲五月天丁香| 我的亚洲天堂| 久久婷婷成人综合色麻豆| 国产一区二区三区综合在线观看| 黄色片一级片一级黄色片| 欧美 亚洲 国产 日韩一| 丁香欧美五月| 色精品久久人妻99蜜桃| 国产在线精品亚洲第一网站| 九色亚洲精品在线播放| 亚洲第一青青草原| 男人操女人黄网站| 一级毛片精品| 日韩免费av在线播放| 欧美激情高清一区二区三区| 国产亚洲欧美精品永久| 757午夜福利合集在线观看| 日韩欧美一区二区三区在线观看 | 50天的宝宝边吃奶边哭怎么回事| 身体一侧抽搐| 99精国产麻豆久久婷婷| 女人被狂操c到高潮| 国产精品98久久久久久宅男小说| 精品乱码久久久久久99久播| 99精品久久久久人妻精品| 成人永久免费在线观看视频| 久久精品aⅴ一区二区三区四区| 成人手机av| 王馨瑶露胸无遮挡在线观看| 久久狼人影院| 久久午夜亚洲精品久久| 91在线观看av| 一区二区日韩欧美中文字幕| 男人的好看免费观看在线视频 | av天堂在线播放| 久久精品人人爽人人爽视色| 人人妻人人澡人人看| 无限看片的www在线观看| 国产精品国产av在线观看| 欧美人与性动交α欧美精品济南到| 久久精品国产综合久久久| 亚洲精品久久成人aⅴ小说| 国产精品偷伦视频观看了| 亚洲成人手机| 国产又色又爽无遮挡免费看| 亚洲一区二区三区欧美精品| 中文字幕av电影在线播放| 色在线成人网| 又大又爽又粗| 国产99久久九九免费精品| 午夜影院日韩av| 在线观看午夜福利视频| 在线观看免费高清a一片| 国产又色又爽无遮挡免费看| 欧美乱妇无乱码| 女人爽到高潮嗷嗷叫在线视频| 国产主播在线观看一区二区| 久久午夜亚洲精品久久| 看黄色毛片网站| 亚洲国产看品久久| 男人操女人黄网站| 亚洲av片天天在线观看| 国产又爽黄色视频| 国产av又大| 少妇猛男粗大的猛烈进出视频| 女人高潮潮喷娇喘18禁视频| 中国美女看黄片| 中文字幕高清在线视频| 91成人精品电影| 桃红色精品国产亚洲av| 一级毛片高清免费大全| 人人妻,人人澡人人爽秒播| а√天堂www在线а√下载 | 日韩精品免费视频一区二区三区| 久久精品91无色码中文字幕| 亚洲 欧美一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品91无色码中文字幕| 免费一级毛片在线播放高清视频 | www.熟女人妻精品国产| 搡老岳熟女国产| 精品无人区乱码1区二区| 国产一区二区三区视频了| 一级,二级,三级黄色视频| 久久精品亚洲精品国产色婷小说| 国产精品九九99| 久久ye,这里只有精品| 精品久久久久久,| 如日韩欧美国产精品一区二区三区| 777久久人妻少妇嫩草av网站| av视频免费观看在线观看| 精品欧美一区二区三区在线| 国产在线一区二区三区精| 亚洲人成电影观看| 亚洲av片天天在线观看| 宅男免费午夜| 免费不卡黄色视频| 国产亚洲一区二区精品| 美女国产高潮福利片在线看| 女人被躁到高潮嗷嗷叫费观| 免费观看人在逋| 亚洲综合色网址| 丰满饥渴人妻一区二区三| 午夜福利在线免费观看网站| 国产精品一区二区在线不卡| 午夜影院日韩av| 久久精品亚洲av国产电影网| 成年人午夜在线观看视频| 精品国产美女av久久久久小说| 国产成人一区二区三区免费视频网站| 亚洲专区国产一区二区| 国产免费现黄频在线看| 午夜免费鲁丝| 91麻豆精品激情在线观看国产 | 久热爱精品视频在线9| 这个男人来自地球电影免费观看| 亚洲精品中文字幕在线视频| 国产成人一区二区三区免费视频网站| 日韩成人在线观看一区二区三区| 欧美中文综合在线视频| 一进一出好大好爽视频| 精品国内亚洲2022精品成人 | 国产在线精品亚洲第一网站| 国产精品九九99| 亚洲av成人一区二区三| 中文字幕另类日韩欧美亚洲嫩草| 9热在线视频观看99| 免费看十八禁软件| 老汉色∧v一级毛片| 欧美日韩福利视频一区二区| 人人妻,人人澡人人爽秒播| 在线av久久热| 国产精品香港三级国产av潘金莲| 国产高清国产精品国产三级| 亚洲精品美女久久av网站| 男男h啪啪无遮挡| 久久中文字幕一级| 久久久久国产精品人妻aⅴ院 | 午夜精品国产一区二区电影| 性少妇av在线| 成年动漫av网址| 欧美精品av麻豆av| 精品人妻熟女毛片av久久网站| 在线观看免费高清a一片| 国产av精品麻豆| 色婷婷久久久亚洲欧美| 一级,二级,三级黄色视频| 天堂俺去俺来也www色官网| 午夜精品久久久久久毛片777| 国产成人欧美| 操美女的视频在线观看| 午夜福利乱码中文字幕| 50天的宝宝边吃奶边哭怎么回事| 午夜福利在线免费观看网站| 亚洲午夜精品一区,二区,三区| 亚洲精华国产精华精| 精品国产亚洲在线| 国产麻豆69| 日韩成人在线观看一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 国产一卡二卡三卡精品| 国产欧美日韩一区二区精品| 中文亚洲av片在线观看爽 | 涩涩av久久男人的天堂| 国产单亲对白刺激| 黄色女人牲交| 搡老乐熟女国产| 精品人妻在线不人妻| 久久久国产成人精品二区 | 欧美人与性动交α欧美精品济南到| 麻豆乱淫一区二区| 国产精品久久视频播放| 亚洲五月色婷婷综合| 久久99一区二区三区| 久久精品国产综合久久久| 99精国产麻豆久久婷婷| 9191精品国产免费久久| 老司机影院毛片| 女人爽到高潮嗷嗷叫在线视频| 国产精品av久久久久免费| 亚洲自偷自拍图片 自拍| 亚洲九九香蕉| 韩国av一区二区三区四区| 午夜亚洲福利在线播放| 黑人欧美特级aaaaaa片| 国产高清激情床上av| 法律面前人人平等表现在哪些方面| 国产不卡一卡二| 欧美日韩亚洲国产一区二区在线观看 | 国产精品香港三级国产av潘金莲| 欧美另类亚洲清纯唯美| avwww免费| √禁漫天堂资源中文www| 村上凉子中文字幕在线| 中文字幕色久视频| 丁香欧美五月| 高清av免费在线| 久热这里只有精品99| 成人18禁在线播放| 一区二区日韩欧美中文字幕| 成熟少妇高潮喷水视频| 巨乳人妻的诱惑在线观看| 女人被狂操c到高潮| 午夜影院日韩av| 操美女的视频在线观看| 无遮挡黄片免费观看| 久久久久国产一级毛片高清牌| 一级a爱片免费观看的视频| 久久久国产一区二区| 亚洲精品久久成人aⅴ小说| 黑人猛操日本美女一级片| 他把我摸到了高潮在线观看| 国产精品免费大片| 午夜福利一区二区在线看| 黄色丝袜av网址大全| 99热只有精品国产| 精品久久久久久久久久免费视频 | 飞空精品影院首页| 99热国产这里只有精品6| 欧美日韩黄片免| 久久久国产成人免费| 午夜视频精品福利| 国产精品美女特级片免费视频播放器 | 一级,二级,三级黄色视频| 精品久久久精品久久久| 大陆偷拍与自拍| 午夜91福利影院| 大陆偷拍与自拍| 91九色精品人成在线观看| 日韩熟女老妇一区二区性免费视频| 久久久久久久久免费视频了| 亚洲,欧美精品.| 最新美女视频免费是黄的| 人人妻人人爽人人添夜夜欢视频| 久久热在线av| 亚洲成av片中文字幕在线观看| 久久青草综合色| 久久香蕉国产精品| 91精品国产国语对白视频| 久久久久久久午夜电影 | 久久精品国产综合久久久| 99国产综合亚洲精品| 无限看片的www在线观看| 黄频高清免费视频| 在线十欧美十亚洲十日本专区| 久久久久久久精品吃奶| 黄色成人免费大全| 亚洲精华国产精华精| 亚洲三区欧美一区| 如日韩欧美国产精品一区二区三区| 国产一区二区三区视频了| 久9热在线精品视频| 天天躁夜夜躁狠狠躁躁| 免费在线观看黄色视频的| 国产又爽黄色视频| 窝窝影院91人妻| 亚洲国产精品一区二区三区在线| 12—13女人毛片做爰片一| 色综合欧美亚洲国产小说| 日韩三级视频一区二区三区| 9热在线视频观看99| 黄色片一级片一级黄色片| 黄色女人牲交| 韩国av一区二区三区四区| 男人的好看免费观看在线视频 | 人成视频在线观看免费观看| 水蜜桃什么品种好| 欧美日韩成人在线一区二区| 国产精品综合久久久久久久免费 | 韩国精品一区二区三区| 国产97色在线日韩免费| 成年动漫av网址| 国产aⅴ精品一区二区三区波| 国产欧美日韩精品亚洲av| 亚洲午夜理论影院| 51午夜福利影视在线观看| 久热爱精品视频在线9| tocl精华| 亚洲 国产 在线| 可以免费在线观看a视频的电影网站| 午夜影院日韩av| 一级片'在线观看视频| 国产高清激情床上av| 最新美女视频免费是黄的| 在线天堂中文资源库| 真人做人爱边吃奶动态| 精品久久久久久久毛片微露脸| 搡老熟女国产l中国老女人| 天堂中文最新版在线下载| 亚洲一区中文字幕在线| 一边摸一边抽搐一进一小说 | 国产色视频综合| ponron亚洲| 国产精品亚洲av一区麻豆| 久久狼人影院| 久久国产精品人妻蜜桃| 大香蕉久久网| av电影中文网址| 一级片免费观看大全| 久久久精品区二区三区| 国产精品美女特级片免费视频播放器 | 久久精品国产亚洲av香蕉五月 | 久久久久精品人妻al黑| 亚洲成人手机| 在线看a的网站| 999久久久精品免费观看国产| 亚洲第一欧美日韩一区二区三区| 亚洲五月色婷婷综合| 精品国产一区二区久久| 亚洲av美国av| 岛国在线观看网站| 久久这里只有精品19| 人人妻人人爽人人添夜夜欢视频| 美女国产高潮福利片在线看| 免费在线观看亚洲国产| 亚洲久久久国产精品| 下体分泌物呈黄色| 99久久精品国产亚洲精品| 国产欧美日韩精品亚洲av| 亚洲欧美精品综合一区二区三区| 在线av久久热| 美女高潮到喷水免费观看| 亚洲熟女毛片儿| 高清av免费在线| 极品教师在线免费播放| 悠悠久久av| 国产精品欧美亚洲77777| 精品国产乱码久久久久久男人| 老熟妇仑乱视频hdxx| 久久ye,这里只有精品| 男女免费视频国产| ponron亚洲| 国产精品一区二区免费欧美|