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

    洞庭湖浮游藻類(lèi)功能群的組成特征及其影響因素研究

    2016-12-23 08:02:30李利強(qiáng)鄭丙輝王麗婧中國(guó)環(huán)境科學(xué)研究院環(huán)境基準(zhǔn)與風(fēng)險(xiǎn)評(píng)估國(guó)家重點(diǎn)實(shí)驗(yàn)室北京100012中國(guó)環(huán)境科學(xué)研究院國(guó)家環(huán)境保護(hù)飲用水水源地保護(hù)重點(diǎn)實(shí)驗(yàn)室北京100012國(guó)家環(huán)境保護(hù)洞庭湖科學(xué)觀測(cè)研究站湖南岳陽(yáng)1000洞庭湖生態(tài)環(huán)境監(jiān)測(cè)中心湖南岳陽(yáng)1000
    中國(guó)環(huán)境科學(xué) 2016年12期
    關(guān)鍵詞:優(yōu)勢(shì)環(huán)境功能

    汪 星,李利強(qiáng),鄭丙輝*,劉 琰,田 琪,王麗婧(1.中國(guó)環(huán)境科學(xué)研究院,環(huán)境基準(zhǔn)與風(fēng)險(xiǎn)評(píng)估國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100012;2.中國(guó)環(huán)境科學(xué)研究院,國(guó)家環(huán)境保護(hù)飲用水水源地保護(hù)重點(diǎn)實(shí)驗(yàn)室,北京100012;3.國(guó)家環(huán)境保護(hù)洞庭湖科學(xué)觀測(cè)研究站,湖南 岳陽(yáng) 1000;.洞庭湖生態(tài)環(huán)境監(jiān)測(cè)中心,湖南 岳陽(yáng)1000)

    洞庭湖浮游藻類(lèi)功能群的組成特征及其影響因素研究

    汪 星1,2,3,李利強(qiáng)4,鄭丙輝1,2,3*,劉 琰1,2,田 琪4,王麗婧1,2,3(1.中國(guó)環(huán)境科學(xué)研究院,環(huán)境基準(zhǔn)與風(fēng)險(xiǎn)評(píng)估國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100012;2.中國(guó)環(huán)境科學(xué)研究院,國(guó)家環(huán)境保護(hù)飲用水水源地保護(hù)重點(diǎn)實(shí)驗(yàn)室,北京100012;3.國(guó)家環(huán)境保護(hù)洞庭湖科學(xué)觀測(cè)研究站,湖南 岳陽(yáng) 414000;4.洞庭湖生態(tài)環(huán)境監(jiān)測(cè)中心,湖南 岳陽(yáng)414000)

    2013年3月、6月、9月及12月采集了洞庭湖11個(gè)斷面的浮游藻類(lèi),并對(duì)所采集的藻類(lèi)進(jìn)行了功能群劃分,結(jié)果顯示,洞庭湖的浮游藻類(lèi)可以分為24個(gè)功能群:A、B、C、D、E、F、G、H1、J、LO、LM、M、MP、N、P、S1、SN、T、W1、W2、X1、X2、X3、Y,所反映的生境特征主要表現(xiàn)為對(duì)分層敏感、頻繁擾動(dòng)且混合程度較高的渾濁型中-富營(yíng)養(yǎng)型水體.不同時(shí)期調(diào)查的藻類(lèi)優(yōu)勢(shì)功能群存在明顯差異,其中,B(中營(yíng)養(yǎng)水體、對(duì)分層敏感)、D(較渾濁的淺水水體)、J(混合型高富營(yíng)養(yǎng)淺水水體)、MP(擾動(dòng)頻繁的渾濁型淺水水體)、P(混合程度較高中富營(yíng)養(yǎng)淺水水體)、Y(廣適性)在四次調(diào)查中的優(yōu)勢(shì)度均>0.02,成為洞庭湖的絕對(duì)優(yōu)勢(shì)功能群,洞庭湖藻類(lèi)優(yōu)勢(shì)功能群不同時(shí)期的演替規(guī)律為:3月MP+P+D+B經(jīng)6月MP+J+P+D+B與9月MP+J+P+D+B+LO+Y轉(zhuǎn)變成12月MP+J+D+Y.CCA分析結(jié)果顯示,洞庭湖藻類(lèi)功能群分布受水環(huán)境因子影響較為明顯.整體上,T、pH值、DO、CODCr、NH3-N及TP是影響洞庭湖藻類(lèi)功能群分布格局的主要因素.

    洞庭湖;藻類(lèi)功能群;優(yōu)勢(shì)度;環(huán)境因子

    同一類(lèi)群中的不同種類(lèi)可能具備相似的生態(tài)屬性,因此,生態(tài)學(xué)家們將具備相似生態(tài)屬性的生物體進(jìn)行歸類(lèi),以此來(lái)簡(jiǎn)化復(fù)雜的生態(tài)系統(tǒng),并將該簡(jiǎn)化后的生態(tài)群落定義為功能群[1].浮游藻類(lèi)是一類(lèi)極其多樣、多源的光合自養(yǎng)型原生動(dòng)物或藍(lán)細(xì)菌,是食物鏈和生物地球化學(xué)循環(huán)過(guò)程的初級(jí)生產(chǎn)者[2],20世紀(jì) 70年代以來(lái),功能群分類(lèi)法在浮游藻類(lèi)中的應(yīng)用越來(lái)越廣泛,且不斷被更新.截至目前,浮游藻類(lèi)功能群的分類(lèi)標(biāo)準(zhǔn)包含尺寸、形狀、生長(zhǎng)策略、生態(tài)功能及生境特征等5個(gè) 方 面 .BSS(Biomass size spectrum)/NBS (Normalised Biomass Sizespectrum)、 TTSS (Traditional Taxonomic Size spectrum)以及PGS (Phytoplankton Geometric Shapes)法均以生物體的尺寸和形狀作為分類(lèi)依據(jù),這三種方法的優(yōu)點(diǎn)在于不需要專(zhuān)業(yè)的分類(lèi)學(xué)基礎(chǔ)且錯(cuò)誤歸類(lèi)的風(fēng)險(xiǎn)極低,缺點(diǎn)則是不能反映生物樣本的生態(tài)屬性以及所處生境的特征等[3-6].CSR(Competitive, Stress-tolerant and Ruderal strategists)法是 r/k (rete of increase/carrying capacity)法的衍生,主要反映藻類(lèi)的生長(zhǎng)屬性及策略,優(yōu)點(diǎn)也是對(duì)分類(lèi)學(xué)基礎(chǔ)要求不高,且錯(cuò)誤歸類(lèi)的風(fēng)險(xiǎn)較低,受到研究者們的廣泛關(guān)注并一直沿用至今[7].隨著研究的不斷深入,FG(Functional Groups)法的提出具有劃時(shí)代的意義,該方法能反映藻類(lèi)的棲息地屬性、對(duì)環(huán)境的耐受能力以及水體的營(yíng)養(yǎng)等級(jí)等各個(gè)方面[8-10],近年來(lái)也成為我國(guó)湖泊研究者們使用最為廣泛的方法[11-16].

    洞庭湖位于湖南省北部(28°30′N(xiāo)~30°20′N(xiāo)、111°40′E~113°40′E),長(zhǎng)江中游南岸,西南有湘、資、沅、澧“四水”入湖,北有長(zhǎng)江的藕池、松滋、太平“三口”流入,湖水經(jīng)北部的城陵磯與長(zhǎng)江相通,是我國(guó)第二大淡水湖泊、國(guó)內(nèi)最大淡水濕地以及國(guó)際重要濕地和珍稀候鳥(niǎo)越冬棲息地,亦是長(zhǎng)江洄游和半洄游魚(yú)類(lèi)的產(chǎn)卵場(chǎng)、索餌場(chǎng),具有調(diào)蓄、滯洪、維持長(zhǎng)江中游水域生態(tài)平衡及下游水資源戰(zhàn)略安全等多項(xiàng)生態(tài)服務(wù)功能.因長(zhǎng)年的泥沙淤積和圍墾,洞庭湖被分割為東洞庭湖、南洞庭湖和西洞庭湖.現(xiàn)有天然湖面2 691km2,洪道1 307km2.在自然和人類(lèi)活動(dòng)的雙重作用下,100多年來(lái),湖泊迅速萎縮,生物多樣性顯著下降,部分水域出現(xiàn)藍(lán)藻水華,水情和生態(tài)環(huán)境日益惡化.在以往的研究中曾提出溶解氧和氮是影響洞庭湖藻類(lèi)分布格局的主要因素[17-18],但浮游藻類(lèi)的生長(zhǎng)繁殖及水華的產(chǎn)生也與生境、水文、水動(dòng)力等條件有著潛在的聯(lián)系[15],因此,對(duì)洞庭湖浮游藻類(lèi)生長(zhǎng)、形態(tài)及生態(tài)功能等的演替規(guī)律及其影響因素的研究顯得尤為重要.本研究采用 FG法對(duì)洞庭湖浮游藻類(lèi)功能群特征進(jìn)行分析,并通過(guò)典范對(duì)應(yīng)分析(CCA)明確影響浮游藻類(lèi)功能群演替的主要環(huán)境因子,以期為洞庭湖的生態(tài)與水質(zhì)監(jiān)管提供理論依據(jù)和技術(shù)支撐.

    1 材料與方法

    1.1 采樣位點(diǎn)設(shè)置與采樣方法

    本研究共設(shè)置11個(gè)采樣斷面,如圖1所示,監(jiān)測(cè)時(shí)間為2013年3、6、9和12月.藻類(lèi)樣品的采集使用國(guó)際標(biāo)準(zhǔn)的25號(hào)篩絹制成的小型浮游生物網(wǎng),于水體自下向上垂直托取定量樣品,用 5%的甲醛溶液固定保存.另取表層水樣 2L,搖勻,一部分原水樣直接分裝于250ml磨口瓶,用于測(cè)定總磷、總氮、氨氮、生化需氧量等常規(guī)水質(zhì)指標(biāo).

    1.2 水樣測(cè)試方法

    水樣的生化需氧量(BOD)、高錳酸鹽指數(shù)(CODMn)、化學(xué)需氧量(CODCr)、總氮(TN)、總磷(TP)、氨氮(NH3-N)以及糞大腸菌群(E.coli)的測(cè)定方法參照《水和廢水監(jiān)測(cè)分析方法》[19].其中,稀釋與接種法測(cè)定 BOD,酸性高錳酸鉀法測(cè)定CODMn;重鉻酸鹽法測(cè)定 CODCr;堿性過(guò)硫酸鉀-紫外分光光度法測(cè)定 TN;鉬酸銨分光光度法測(cè)定TP;納氏試劑分光光度法測(cè)定NH3-N;多管發(fā)酵法測(cè)定E.coli.pH值、水溫(T)、溶解氧(DO)的測(cè)定采用便攜式 pH 計(jì)和水質(zhì)多參數(shù)測(cè)定儀(YSI)現(xiàn)場(chǎng)直接測(cè)定(水面下0.5m);參照《中國(guó)淡水藻志》[20]對(duì)采集的樣品進(jìn)行分類(lèi)、鑒定,藻細(xì)胞生物量計(jì)算等采用顯微鏡法[21].

    圖1 洞庭湖采樣斷面示意Fig.1 Sampling sections of Dongting Lake

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

    1.3.1 浮游藻類(lèi)功能群分類(lèi)及生物量計(jì)算 參照2002年Reynolds等[9]和2009年P(guān)adisák等[10]以浮游植物生理生長(zhǎng)特征及其環(huán)境適應(yīng)性機(jī)制為基礎(chǔ),在浮游植物生長(zhǎng)的C-R-S策略分類(lèi)基礎(chǔ)上,對(duì)具有相同適應(yīng)性特征、易于在相同生境條件下共存的浮游植物種類(lèi)進(jìn)行功能分組.浮游藻類(lèi)功能群生物量指的是功能群的代表性藻種(屬)的細(xì)胞生物量之和,通過(guò)各代表性藻種(屬)的生物量相加得出.比如洞庭湖藻類(lèi)功能群B的代表性藻種(屬)包括小環(huán)藻和冠盤(pán)藻,那么該功能群生物量則為小環(huán)藻和冠盤(pán)藻細(xì)胞生物量之和.

    1.3.2 浮游藻類(lèi)功能群優(yōu)勢(shì)度分析 優(yōu)勢(shì)度根據(jù)藻類(lèi)功能群的出現(xiàn)頻率及豐度來(lái)確定,用優(yōu)勢(shì)度來(lái)表示[22]:y=fi·Pi.式中:y為優(yōu)勢(shì)度,fi為第i功能群的出現(xiàn)頻率,Pi為第 i功能群豐度占總藻類(lèi)豐度的比例,當(dāng) y>0.02時(shí),定為優(yōu)勢(shì)功能群[23].選擇優(yōu)勢(shì)度>0.02的浮游藻類(lèi)功能群進(jìn)行典范對(duì)應(yīng)(CCA)分析.

    1.3.3 環(huán)境因子主成分分析(PCA) 采用主成分分析(PCA)中的因子分析法[24],將各采樣斷面的環(huán)境數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理(其中溶解氧為逆指標(biāo),故將其先進(jìn)行倒數(shù)變換,然后在 SPSS13.0中完成數(shù)據(jù)標(biāo)準(zhǔn)化),為了使每個(gè)主成分的意義更加明確,采用 Varimax最大方差法對(duì)因子進(jìn)行旋轉(zhuǎn),按照特征值λ大于1的原則提取出3個(gè)主成分,選擇主成分的因子載荷量>0.7的環(huán)境因子進(jìn)行CCA分析.所有的數(shù)據(jù)分析均在SPSS 13.0軟件下進(jìn)行.

    1.3.4 典范對(duì)應(yīng)分析(CCA) 采用 Canoco for Windows 4.5 軟件對(duì)功能群數(shù)據(jù)和環(huán)境數(shù)據(jù)進(jìn)行CCA分析.浮游藻類(lèi)功能群矩陣經(jīng)過(guò) lg(x+1)轉(zhuǎn)換,環(huán)境因子數(shù)據(jù)除 pH值外全部進(jìn)行l(wèi)g(x+1)轉(zhuǎn)換[25].CCA分析結(jié)果用功能群—環(huán)境因子關(guān)系的雙序圖表示,圖中環(huán)境因子用帶有箭頭的線段表示,向量長(zhǎng)短表示其在主軸中的作用,線段所處象限代表環(huán)境因子與排序軸間的正負(fù)關(guān)系.分析時(shí),做出某一功能群與環(huán)境因子連線的垂直線,垂直線與環(huán)境因子連線相交點(diǎn)離箭頭越近,表示該種與該類(lèi)生境因子的正相關(guān)性越大,處于另一端的則表示與該類(lèi)環(huán)境因子具有的負(fù)相關(guān)性越大.

    2 結(jié)果分析

    2.1 浮游藻類(lèi)豐度變化及功能群劃分

    從不同時(shí)期藻類(lèi)生物量的分布來(lái)看(圖 2),3月藻類(lèi)的生物量變化范圍為 0.08~2.59mg/L,最大值出現(xiàn)在大小西湖(S11),而萬(wàn)子湖(S5)的生物量最小,整體來(lái)看,浮游藻類(lèi)生物量從大到小依次為東洞庭湖、西洞庭湖、南洞庭湖;6月、9月和12月浮游藻類(lèi)的生物量變化范圍分別為0.05~0.89,0.07~1.76,0.02~2.41mg/L,最大值均出現(xiàn)在大小西湖(S11),而坡頭(S1)的生物量最小,整體來(lái)看,浮游藻類(lèi)生物量從大到小依次為東洞庭湖、南洞庭湖、西洞庭湖.

    對(duì)所檢出的浮游藻類(lèi)進(jìn)行功能分類(lèi),如表 1所示,洞庭湖的浮游藻類(lèi)可以分為 24個(gè)功能群:A、B、C、D、E、F、G、H1、J、LO、LM、M、MP、N、P、S1、SN、T、W1、W2、X1、X2、X3、Y.

    圖2 洞庭湖2013年浮游藻類(lèi)生物量的時(shí)空變化Fig.2 Temporal and spatial variation of phytoplankton of Dongting Lake in 2013

    表1 洞庭湖浮游藻類(lèi)功能群劃分Table 1 Algal function groups division in Dongting Lake

    續(xù)表1

    2.2 浮游藻類(lèi)優(yōu)勢(shì)功能群及其分布特征

    從不同時(shí)期藻類(lèi)功能群的優(yōu)勢(shì)度來(lái)看(圖3),3月浮游藻類(lèi)的優(yōu)勢(shì)功能群為D、B、MP、P、J、Y、X1,6月為J、MP、B、Y、D、P,9月為J、Y、P、MP、LO、D、B、X2,12月為Y、D、J、MP、X2、B 、P、X1,其中,B、D、J、MP、P、Y在4次調(diào)查中的優(yōu)勢(shì)度均>0.02,成為洞庭湖的絕對(duì)優(yōu)勢(shì)功能群.本研究中浮游藻類(lèi)功能群的時(shí)空分布特征及其影響因子分析均在上述優(yōu)勢(shì)功能群中展開(kāi).

    從不同時(shí)期藻的功能群組成來(lái)看(圖4),3月藻類(lèi)功能群 MP、P、B、D在各采樣斷面均有分布,其中,B在南嘴(S3)、鹿角(S8)和洞庭湖出口(S10)均占據(jù)絕對(duì)優(yōu)勢(shì),D在大小西湖(S11)占據(jù)絕對(duì)優(yōu)勢(shì),整體來(lái)看,MP、P自西洞庭湖至東洞庭湖湖呈現(xiàn)逐漸下降的趨勢(shì),而B(niǎo)、D則呈現(xiàn)逐漸升高的趨勢(shì)(除 S3外),值得注意的是,萬(wàn)子湖(S4)以Y占據(jù)優(yōu)勢(shì),而MP、P、B、D等優(yōu)勢(shì)類(lèi)群所占比例較小.6月各采樣斷面的主要功能群由3月的MP+P+D+B轉(zhuǎn)變成MP+J+P+D+B,優(yōu)勢(shì)功能類(lèi)型較3月明顯增多,其中,MP、J分別在東洞庭湖(S9)和大小西湖(S11)占據(jù)絕對(duì)優(yōu)勢(shì),整體來(lái)看,西洞庭湖與南洞庭湖的藻類(lèi)優(yōu)勢(shì)功能群組成均衡度較高,而東洞庭湖則以MP、J成為絕對(duì)優(yōu)勢(shì)群落.9月各采樣斷面的主要功能群由 6月的MP+J+P+D+B轉(zhuǎn)變成MP+J+P+D+B+LO+Y,其中,J、Y分別在大小西湖(S11)和洞庭湖出口(S10)占據(jù)絕對(duì)優(yōu)勢(shì),整體來(lái)看,MP、LO自西洞庭湖至東洞庭湖湖呈現(xiàn)逐漸下降的趨勢(shì),J、Y則呈現(xiàn)逐漸上升的趨勢(shì),P、D、B在各采樣斷面所占比例相對(duì)均衡.12月洞庭湖整體以Y占據(jù)絕對(duì)優(yōu)勢(shì),各采樣斷面的主要功能群由 9月的MP+J+P+D+ B+LO+Y轉(zhuǎn)變成MP+J+D+Y,整體來(lái)看,MP自西洞庭湖至東洞庭湖湖呈現(xiàn)逐漸下降的趨勢(shì),而J則呈現(xiàn)相反的趨勢(shì);D在南洞庭湖所占比例明顯高于西洞庭湖和東洞庭湖,而Y在西洞庭湖和東洞庭湖所占比例則明顯高于南洞庭湖.

    圖3 洞庭湖不同季節(jié)浮游藻類(lèi)功能群優(yōu)勢(shì)度Fig.3 Dominance of algal function groups in different seasons in Dongting Lake

    圖4 洞庭湖不同季節(jié)浮游藻類(lèi)功能群組成Fig.4 Algal function groups composition in different seasons in Dongting Lake

    2.3 環(huán)境因子的主成分分析

    為了使每個(gè)主成分的意義更加明確,采用Varimax最大方差法對(duì)因子進(jìn)行旋轉(zhuǎn),按照特征值λ大于1的原則提取出3個(gè)主成分,如表2所示,3月、6月、9月和12月的主成分累計(jì)貢獻(xiàn)率分別達(dá)到81%、94.5%、78.5%及92.8%.該主成分分析所得指標(biāo)的載荷超過(guò) 0.7的即為主要影響因子,3月的主要環(huán)境因子為T(mén)、pH值、DO、CODMn、CODCr、BOD、NH3-N、TP、E.coli,6月的主要環(huán)境因子為T(mén)、pH值、DO、CODMn、CODCr、BOD、NH3-N、TN、E.coli,9月的主要環(huán)境因子為T(mén)、pH值、DO、SD、CODMn、CODCr、NH3-N、TN、E.coli,12月的主要環(huán)境因子為 T、pH值、DO、CODMn、CODCr、BOD、NH3-N、TN、TP、E.coli.

    表2 Varimax轉(zhuǎn)軸后的因子載荷量Table 2 Rotated component matrix

    2.4 浮游藻類(lèi)優(yōu)勢(shì)功能群與環(huán)境因子的相關(guān) (CCA)分析

    圖5 洞庭湖不同季節(jié)浮游藻類(lèi)優(yōu)勢(shì)功能群與主要環(huán)境因子的CCA分析Fig.5 CCA analysis of algal dominant function groups and principal environment factors in different seasons in Dongting Lake

    選取藻類(lèi)優(yōu)勢(shì)功能群與主要環(huán)境因子進(jìn)行CCA分析,如圖5所示,環(huán)境因子第一、二軸間的相關(guān)系數(shù)為 0,表明分析結(jié)果可信[26].3月的分析結(jié)果顯示,Y與DO呈現(xiàn)強(qiáng)正相關(guān); J與pH值呈現(xiàn)強(qiáng)正相關(guān);整體來(lái)看,DO、pH是影響3月藻類(lèi)功能群分布的主要環(huán)境因子.6月的分析結(jié)果顯示,W2、X1與T正相關(guān)性顯著;J、B、P、MP、D與DO、CODCr均呈現(xiàn)強(qiáng)正相關(guān);整體來(lái)看,T、DO及CODCr是影響6月藻類(lèi)功能群分布的主要環(huán)境因子.9月的分析結(jié)果顯示,J與pH呈現(xiàn)強(qiáng)正相關(guān);P與 CODCr、DO呈現(xiàn)強(qiáng)正相關(guān);整體來(lái)看,pH、CODCr及DO是影響9月藻類(lèi)功能群分布的主要環(huán)境因子.12月的分析結(jié)果顯示,P、X1均與 T呈現(xiàn)強(qiáng)正相關(guān);J與 NH3-N呈現(xiàn)強(qiáng)正相關(guān);D、Y、MP、X2均與TP呈現(xiàn)較強(qiáng)的正相關(guān);整體來(lái)看,T、NH3-N、TP是影響12月藻類(lèi)功能群分布的主要環(huán)境因子.

    3 討論

    3.1 浮游藻類(lèi)豐度變化及功能群劃分

    通過(guò)對(duì)洞庭湖不同時(shí)期浮游藻類(lèi)的調(diào)查顯示,該流域內(nèi)藻類(lèi)功能群存在明顯的時(shí)空差異.從優(yōu)勢(shì)功能群的代表性種類(lèi)來(lái)看,3月浮游藻類(lèi)以硅藻占絕對(duì)優(yōu)勢(shì),其次是綠藻和藍(lán)藻,6月以綠藻與硅藻占絕對(duì)優(yōu)勢(shì),其次是藍(lán)藻和隱藻,9月為以綠藻占絕對(duì)優(yōu)勢(shì),其次是甲藻、隱藻,硅藻所占比例顯著減少,12月則轉(zhuǎn)變?yōu)榧自搴碗[藻占優(yōu),其次是硅藻和綠藻.1999年李利強(qiáng)等[27]、2012年汪星等[17]的調(diào)查結(jié)果均顯示,洞庭湖以綠藻門(mén)與硅藻門(mén)占優(yōu),藍(lán)藻門(mén)次之,陳格君等[28]對(duì)另一長(zhǎng)江的典型通江行湖泊鄱陽(yáng)湖的調(diào)查結(jié)果也顯示,浮游藻類(lèi)的組成以綠藻門(mén)占優(yōu),硅藻和藍(lán)藻門(mén)次之,均與本研究的結(jié)果保持一致.可見(jiàn),近年來(lái),盡管土地利用、水產(chǎn)養(yǎng)殖及種植業(yè)等人類(lèi)活動(dòng)的不斷增多,洞庭湖藻類(lèi)功能群的組成并未發(fā)生明顯改變.

    從優(yōu)勢(shì)功能群的生境特征來(lái)看,3月西洞庭湖以B和MP占絕對(duì)優(yōu)勢(shì),特別地,南嘴(S3)的B所占比例超過(guò)80%,B與MP功能群的生境特征均為頻繁擾動(dòng)渾濁型淺水水體,這與西洞庭湖的幾個(gè)采樣斷面實(shí)際情況相符,S1、S3分別屬于沅水和澧水的入湖監(jiān)測(cè)斷面,上游的開(kāi)江春水夾帶大量泥沙且流速較大,導(dǎo)致水體較為渾濁,同樣地,S2、S4也受益陽(yáng)資水的來(lái)水影響,水體混合度較高;南洞庭湖的萬(wàn)子湖(S5)與橫嶺湖(S6)斷面受入湖河流擾動(dòng)較小,水體流動(dòng)速率相對(duì)較緩,適宜生境為混合程度較高的中-富營(yíng)養(yǎng)型水體的 P型藻類(lèi)生長(zhǎng),這也與洞庭湖水體屬于中-輕度富營(yíng)養(yǎng)狀態(tài)的報(bào)道一致[29-30].S7-S8、S9-S11分別受湘江與長(zhǎng)江回水的作用,水體重回渾濁,適宜B、D及MP等功能群的生長(zhǎng).6月以MP和J在各采樣斷面的優(yōu)勢(shì)較為明顯,受上游來(lái)水的影響,MP在坡頭、南嘴、鹿角、東洞庭湖及洞庭湖出口的比例仍占據(jù)優(yōu)勢(shì); J的適宜生境為混合型高富營(yíng)養(yǎng)水體,其比例則呈現(xiàn)自西洞庭湖到東洞庭湖逐漸增加的趨勢(shì),并在大小西湖達(dá)到最大;鐘振宇等[29]、周泓等[31]、王麗婧等[32]、Wang等[33]的研究結(jié)果均表明,洞庭湖的水體富營(yíng)養(yǎng)化狀態(tài)依次為西洞庭湖<南洞庭湖<東洞庭湖,且大小西湖水體的富營(yíng)養(yǎng)化最為惡劣,這些均與本研究中藻類(lèi)功能群的分布狀況相符.9月MP與J的分布狀況與6月的基本一致,但J的比例有所下降,而Y的比例則隨之增加,Y屬于廣適性中營(yíng)養(yǎng)功能群且在萬(wàn)子湖和洞庭湖出口占據(jù)絕對(duì)優(yōu)勢(shì),究其原因,可能是9月屬于豐水期,湖區(qū)上游來(lái)水量增大,對(duì)水體有一定的稀釋作用,導(dǎo)致整體營(yíng)養(yǎng)狀態(tài)有所下降,廣適性功能群比例增加,值得注意的是,大小西湖仍以 J占據(jù)絕對(duì)優(yōu)勢(shì),其富營(yíng)養(yǎng)化狀態(tài)不容小視.12月MP在東洞庭湖的比例明顯小于西洞庭湖和南洞庭湖,東洞庭湖在該時(shí)期受長(zhǎng)江回水影響明顯減少;Y在各采樣斷面占據(jù)優(yōu)勢(shì)明顯,可見(jiàn)該時(shí)期洞庭湖水體營(yíng)養(yǎng)狀況良好,究其原因,可能是枯水期上游來(lái)水減少,進(jìn)入湖區(qū)的污染物總量降低,加之湖區(qū)的自?xún)糇饔?湖區(qū)適應(yīng)中營(yíng)養(yǎng)和廣適性的藻類(lèi)生長(zhǎng);另外,該時(shí)期J的變化趨勢(shì)仍與其他時(shí)期的保持一致,三個(gè)湖區(qū)的水體營(yíng)養(yǎng)狀態(tài)高低次序保持不變.

    3.2 浮游藻類(lèi)功能群與環(huán)境因子的相關(guān)性

    研究表明,pH、水溫、溶解氧、總氮、透明度、總磷、氨氮、生化需氧量以及化學(xué)需氧量均能成為影響浮游藻類(lèi)群落分布的主要環(huán)境因子

    [34-40],本研究結(jié)果顯示,水溫、pH、DO、CODCr、NH3-N及TP是影響洞庭湖藻類(lèi)功能群分布的主要環(huán)境因子,其中,水溫成為P、W2、X1等多個(gè)功能群分布的主要影響因素,雖然洞庭湖屬于北亞熱帶季風(fēng)濕潤(rùn)氣候,但其在不同季節(jié)以及晝夜的水溫仍存在明顯差異,因此,水溫必然成為該區(qū)域藻類(lèi)生長(zhǎng)的限制因子,對(duì)寒區(qū)安邦河濕地[13]、北方的烏梁素海[15]、深水型金盆水庫(kù)[41]、貴州高原的阿哈水庫(kù)[14]與三板溪水庫(kù)[42]等湖庫(kù)的研究亦表明,水溫是該區(qū)域藻類(lèi)功能群結(jié)構(gòu)變化的最主要影響因素,可見(jiàn)水溫對(duì)湖泊藻類(lèi)生長(zhǎng)及結(jié)構(gòu)演替的重要性.藻類(lèi)適宜在偏堿性的水體中生長(zhǎng),可見(jiàn) pH對(duì)藻類(lèi)生長(zhǎng)的重要性,本研究中,pH則成為適宜富營(yíng)養(yǎng)水體生長(zhǎng)的 J功能群的限制因子.藻類(lèi)的生長(zhǎng)除了光合作用,同時(shí)也需要氧氣來(lái)完成呼吸作用,因此,DO成為B、D、J、MP、P、Y分布的主要限制因素,已有研究表明,DO與藻類(lèi)生長(zhǎng)呈顯著正相關(guān)[43],作為典型的通江湖泊,鏡泊湖的浮游藻類(lèi)也受到 DO的限制[34].CODCr通常作為反映點(diǎn)源污染的常規(guī)檢測(cè)指標(biāo),NH3-N、TP則是反映面源污染的兩項(xiàng)重要監(jiān)測(cè)指標(biāo).南屏水庫(kù)、竹仙洞水庫(kù)[44]、北京海子水庫(kù)[45]的研究表明, 磷是浮游藻類(lèi)功能群演替的主要驅(qū)動(dòng)因子,本研究中,J、B、P、MP、D均受CODCr的限制,可見(jiàn)有機(jī)污染物對(duì)藻類(lèi)生長(zhǎng)的影響之深,同時(shí),NH3-N、TP兩項(xiàng)營(yíng)養(yǎng)鹽指標(biāo)也分別與適應(yīng)富營(yíng)養(yǎng)水體生長(zhǎng)的 J、X2呈現(xiàn)強(qiáng)的正相關(guān);研究表明[18],洞庭湖浮游藻類(lèi)的生長(zhǎng)主要受到氮(尤其是氨態(tài)氮)的限制,添加氨氮能明顯促進(jìn)浮游藻類(lèi)生物量的增加.洞庭湖作為長(zhǎng)江中下游通江型調(diào)蓄湖泊,承接“四水”和“三口”,受上游來(lái)水影響顯著,同時(shí),湖區(qū)周邊的人為活動(dòng)也對(duì)洞庭湖的水質(zhì)狀況帶來(lái)了深遠(yuǎn)的影響,整體來(lái)看,洞庭湖的污染主要來(lái)自于上游的工業(yè)廢水和生活污水、湖區(qū)周邊造紙廠廢水、農(nóng)業(yè)用肥以及湖區(qū)內(nèi)漁民的生活污水與垃圾滲濾液等,其中工業(yè)廢水則是CODCr的主要來(lái)源,生活污水及農(nóng)用肥料則是NH3-N與TP的主要貢獻(xiàn)者,故CODCr、NH3-N、TP能成為洞庭湖藻類(lèi)功能群分布的主要限制因子.當(dāng)然,隨著社會(huì)經(jīng)濟(jì)的繼續(xù)穩(wěn)步發(fā)展,人類(lèi)活動(dòng)會(huì)更加頻繁,洞庭湖的生態(tài)環(huán)境質(zhì)量問(wèn)題應(yīng)引起有關(guān)部門(mén)重視.

    4 結(jié)論

    4.1 2013年對(duì)洞庭湖 11個(gè)采樣斷面所采集的浮游藻類(lèi)的分類(lèi)結(jié)果顯示,洞庭湖的浮游藻類(lèi)可以分為24個(gè)功能群:A、B、C、D、E、F、G、H1、J、LO、LM、M、MP、N、P、S1、SN、T、W1、W2、X1、X2、X3、Y.

    4.2 洞庭湖不同時(shí)期的藻類(lèi)優(yōu)勢(shì)功能群存在明顯差異,其中,B、D、J、MP、P、Y在4次調(diào)查中的優(yōu)勢(shì)度均>0.02,成為洞庭湖的絕對(duì)優(yōu)勢(shì)功能群;洞庭湖藻類(lèi)優(yōu)勢(shì)功能群不同季節(jié)的演替規(guī)律為:3月MP+P+D+B經(jīng)6月MP+J+P+D+B與9月 MP+J+P+D+B+LO+Y 轉(zhuǎn) 變 成 12 月MP+J+D+Y.

    4.3 CCA分析結(jié)果顯示,洞庭湖藻類(lèi)功能群分布受水環(huán)境因子影響較為明顯.整體上,T、pH值、DO、CODCr、NH3-N及TP是影響洞庭湖藻類(lèi)功能群分布格局的主要因素.

    [1] Solbrig O T. Plant traits and adaptive strategies: their role in ecosystem function [M]//: Biodiversity and Ecosystem Function (Eds E.D. Schulze & H.A. Mooney), 1993:97-116. Ecological Studies. Springer-Verlag, Berlin.

    [2] Rousseaux C S, Gregg W W. Interannual variation in phytoplankton primary production at a global scale [J]. Remote Sensing, 2013,6(1):1-19.

    [3] Platt T, Denman K. The structure of pelagic marine ecosystems [J]. Rapports et Proces Verbaux des Reunions, 1978,173:60-65.

    [4] Kamenir Y, Dubinsky Z, Zohary T. Phytoplankton size structure stability in a meso-eutrophic subtropical lake [J]. Hydrobiologia, 2004,520(1-3):89-104.

    [5] Kamenir Y, Dubinsky Z, Zohary T. The long-term patterns of phytoplankton taxonomic size-structure and their sensitivity to perturbation: A Lake Kinneret case study [J]. Aquatic Sciences, 2006,68(4):490-501.

    [6] Stanca E, Cellamare M, Basset A. Geometric shape as a trait to study phytoplankton distributions in aquatic ecosystems [J]. Hydrobiologia, 2013,701:99-116.

    [7] Stanca E, Cellamare M, Basset A. Geometric shape as a trait to study phytoplankton distributions in aquatic ecosystems [J]. Hydrobiologia, 2013,701(1):99-116.

    [8] Reynolds C S. Phytoplankton assemblages and their periodicity in stratifying lake systems [J]. Holarctic Ecology, 1980,3(3):141-159.

    [9] Reynolds C S, Huszar V, Kruk C, et al. Towards a functional classification of the freshwater phytoplankton [J]. Journal of Plankton Research, 2002,24(5):417-428.

    [10] Padisak J, Crossetti L O, Naselli-Flores L. Use and misuse in the application of the phytoplankton functional classification: a critical review with updates [J]. Hydrobiologia, 2009,621(1): 1-19.

    [11] 董 靜,李根保,宋立榮.撫仙湖,洱海,滇池浮游藻類(lèi)功能群1960s以來(lái)演變特征 [J]. 湖泊科學(xué), 2014,26(5):735-742.

    [12] 陳曉江,楊 劼,杜桂森,等.海子水庫(kù)浮游植物功能群季節(jié)演替及其驅(qū)動(dòng)因子 [J]. 水資源保護(hù), 2015,31(6):122-127.

    [13] 武安泉,郭 寧,覃雪波.寒區(qū)典型濕地浮游植物功能群季節(jié)變化及其與環(huán)境因子關(guān)系 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2015,35(5):1341-1349.

    [14] 李 磊,李秋華,焦樹(shù)林,等.阿哈水庫(kù)浮游植物功能群時(shí)空分布特征及其影響因子分析 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2015,35(11):3604-3611.

    [15] 李 興,李建茹,徐效清,等.烏梁素海浮游植物功能群季節(jié)演替規(guī)律及影響因子 [J]. 生態(tài)環(huán)境學(xué)報(bào), 2015,24(10):1668-1675.

    [16] Liu X, Qian K, Chen Y. Effects of water level fluctuations on phytoplankton in a Changjiang River floodplain lake (Poyang Lake): Implications for dam operations [J]. Journal of Great Lakes Research, 2015,41(3):770-779.

    [17] 汪 星,鄭丙輝,劉錄三,等.洞庭湖典型斷面藻類(lèi)組成及其與環(huán)境因子典范對(duì)應(yīng)分析 [J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2012,31(5):995-1002.

    [18] 李利強(qiáng),黃代中,熊 劍,等.洞庭湖浮游植物增長(zhǎng)的限制性營(yíng)養(yǎng)元素研究 [J]. 生態(tài)環(huán)境學(xué)報(bào), 2014,23(2):283-288.

    [19] 國(guó)家環(huán)保總局《水和廢水監(jiān)測(cè)分析方法》編委會(huì).水和廢水監(jiān)測(cè)分析方法 [M]. 4版.北京:中國(guó)環(huán)境出版社, 2002.

    [20] 齊雨藻,李家英.中國(guó)淡水藻志 [M]. 北京:科學(xué)出版社, 2004.

    [21] 徐金森,鄭天凌,郭清華,等.兩種海洋細(xì)菌對(duì)赤潮藻的細(xì)胞生物量的影響研究 [J]. 海洋科學(xué), 2002,26(12):57-60,67.

    [22] Aksnes D L, Wassmann P. Modeling the significance of zooplankton grazing for export production [J]. Limnology and Oceanography, 1993,38(5):978-985.

    [23] Lampitt R S, Wishner K F, Turley C M, et al. Marine snow studies in the Northeast Atlantic Ocean: distribution, composition and role as a food source for migrating plankton [J]. Marine Biology, 1993,116(4):689-702.

    [24] 劉 瀟,薛 瑩,紀(jì)毓鵬,等.基于主成分分析法的黃河口及其鄰近水域水質(zhì)評(píng)價(jià) [J]. 中國(guó)環(huán)境科學(xué), 2015,35(10):3187-3192.

    [25] 王 華,楊樹(shù)平,房晟忠,等.滇池浮游植物群落特征及與環(huán)境因子的典范對(duì)應(yīng)分析 [J]. 中國(guó)環(huán)境科學(xué), 2016,36(2):544-552.

    [26] Kent M. Vegetation description and data analysis: a practical approach [M]. John Wiley & Sons, 2011.

    [27] 李利強(qiáng),張建波.洞庭湖浮游植物調(diào)查與水質(zhì)評(píng)價(jià) [J]. 江蘇環(huán)境科技, 1999,12(4):14-16.

    [28] 陳格君,周文斌,胡春華.鄱陽(yáng)湖五河入湖口浮游藻類(lèi)及營(yíng)養(yǎng)現(xiàn)狀評(píng)價(jià) [J]. 湖北農(nóng)業(yè)科學(xué), 2013,52(9):2048-2052.

    [29] 鐘振宇,陳 燦.洞庭湖水質(zhì)及富營(yíng)養(yǎng)狀態(tài)評(píng)價(jià) [J]. 環(huán)境科學(xué)與管理, 2011,36(7):169-173.

    [30] 黃代中,萬(wàn) 群,李利強(qiáng),等.洞庭湖近20年水質(zhì)與富營(yíng)養(yǎng)化狀態(tài)變化 [J]. 環(huán)境科學(xué)研究, 2013,26(1):27-33.

    [31] 周 泓,歐伏平,劉 妍.“十一五”期間洞庭湖水環(huán)境質(zhì)量狀況及變化趨勢(shì)分析 [J]. 湖南理工學(xué)院學(xué)報(bào)(自然科學(xué)版), 2011, 24(2):88-90.

    [32] 王麗婧,汪 星,劉錄三,等.洞庭湖水質(zhì)因子的多元分析 [J]. 環(huán)境科學(xué)研究, 2013,26(1):1-7.

    [33] Wang X, Zheng B, Liu L, et al. Development and evaluation of the Lake Multi-biotic Integrity Index for Dongting Lake, China [J]. Journal of Limnology, 2015,74(3):594-605.

    [34] 汪 星,劉錄三,李 黎,等.鏡泊湖浮游藻類(lèi)組成及其與環(huán)境因子的相關(guān)分析 [J]. 中國(guó)環(huán)境科學(xué), 2015,35(11):3403-3413.

    [35] 沈會(huì)濤,劉存歧.白洋淀浮游植物群落及其與環(huán)境因子的典范對(duì)應(yīng)分析 [J]. 湖泊科學(xué), 2008,20(1):773-779.

    [36] 劉鴻雁,徐云麟.鏡泊湖藻類(lèi)生長(zhǎng)和湖泊富營(yíng)養(yǎng)化預(yù)測(cè)初探 [J].生態(tài)學(xué)報(bào), 1996,16(2):195-201.

    [37] Jiang Y J, He W, Liu W X, et al. The seasonal and spatial variations of phytoplankton community and their correlation with environmental factors in a large eutrophic chinese lake (Lake Chaohu) [J]. Ecological Indicators, 2014,40:58-67.

    [38] 潘繼征,熊 飛,李文朝,等.云南撫仙湖透明度的時(shí)空變化及影響因子分析 [J]. 湖泊科學(xué), 2008,20(5):681-686.

    [39] Arhonditsis G B, Winder M, Brett M T et al. Patterns and mechanisms of phytoplankton variability in Lake Washington (USA). Water Research, 2004,38:4013-4027.

    [40] 楊麗標(biāo),韓小勇,孫 璞,等.巢湖藻類(lèi)組成與環(huán)境因子典范對(duì)應(yīng)分析 [J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2011,30(5):952-958.

    [41] 盧金鎖,張 穎,胡亞潘.深水型水庫(kù)中藻類(lèi)功能群組演替及其與環(huán)境因子的關(guān)系 [J]. 環(huán)境工程學(xué)報(bào), 2014,(11):4605-4611.

    [42] 黃國(guó)佳,李秋華,陳 椽,等.貴州高原三板溪水庫(kù)浮游植物功能群時(shí)空分布特征 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2015,35(2):418-428.

    [43] 游 亮,崔莉鳳,劉載文,等.藻類(lèi)生長(zhǎng)過(guò)程中 DO、pH與葉綠素相關(guān)性分析 [J]. 環(huán)境科學(xué)與技術(shù), 2007,30(9):42-45.

    [44] 張 怡,胡 韌,肖利娟,等.南亞熱帶兩座不同水文動(dòng)態(tài)的水庫(kù)浮游植物的功能類(lèi)群演替比較 [J]. 生態(tài)環(huán)境學(xué)報(bào), 2012,21(1): 107-117.

    [45] 陳曉江,楊 劼,杜桂森,等.海子水庫(kù)浮游植物功能群季節(jié)演替及其驅(qū)動(dòng)因子 [J]. 水資源保護(hù), 2015,31(6):122-127.

    Composition and influential factors of algal function groups in Dongting Lake.

    WANG Xing1,2,3, LI Li-qiang4, ZHENG Bing-hui1,2,3*, LIU Yan1,2, TIAN Qi4, WANG Li-jing1,2,3(1.State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, Beijing 100012, China;2.State Environmental Protection Key Laboratory of Drink Water Resource Protection, Chinese Research Academy of Environmental Sciences, Beijing 100012, China;3.State Environmental Protection Scientific Observation and Research Station for Lake Dongtinghu, Yueyang 414000, China;4.Dongting Lake Eco-Environment Monitoring Centre, Yueyang 414000, China). China Environmental Science, 2016,36(12):3766~3776

    This research was performed at March, June, September and December in 2013 at 11sections of Dongting Lake, aiming to investigate the algal function groups. A total of 24 kinds of algal function groups were identified: A、B、C、D、E、F、G、H1、J、LO、LM、M、MP、N、P、S1、SN、T、W1、W2、X1、X2、X3、Y, which reflected the habitat characteristics with species sensitive to the onset of stratification, and which are adapt to frequently stirred up, inorganically turbid, and mesotrophic to eutrophic shallow lakes. The algal function groups was dominated by B (mesotrophic waters in small-and mediumsized lakes with species sensitive to the onset of stratification), D (shallow turbid waters), J (well-mixed, enriched shallow lakes), MP (frequently stirred up, inorganically turbid shallow lakes), P (well-mixed, mesotrophic to eutrophic shallow lakes) and Y (refers to a wide range of habitats) based on the dominant degree (>0.02) and significant differences which were found on algal function groups across 11sampling sections in the four investigations. The successional law of algal function groups in different seasons could be summarized as follow: MP+P+D+B in March, MP+J+P+D+B in June, MP+J+P+D+B+LO+Y in September, MP+J+D+Y in December. PCA (principal component analysis) was applied to choose the major potential stressors which influenced the habitat of algae by SPSS 13.0, and CCA (canonical correspondence analysis) was used to explore the relationship between algal function groups and major environmental parameters. It is indicated from the results that water temperature (T), pH, dissolvedoxygen (DO), Chemical Oxygen Demand (CODCr), ammonia nitrogen (NH3-N) and total phosphorus (TP) were the most important factors influencing the distribution of algal function groups across the 11sampling sections.

    Dongting Lake;algal function groups;dominance;environmental factors

    X826

    A

    1000-6923(2016)12-3766-11

    汪 星(1983-),男,湖北荊州人,副研究員,博士,主要從事水生生物評(píng)價(jià)的研究.發(fā)表論文30余篇.

    2016-06-02

    “十二五”國(guó)家科技支撐計(jì)劃課題(2014BAC09B02);“973”項(xiàng)目(2012CB417004);洞庭湖江湖生態(tài)監(jiān)測(cè)重點(diǎn)站項(xiàng)目(JJ2013-023)

    * 責(zé)任作者, 研究員, zhengbinghui@craes.org.cn

    猜你喜歡
    優(yōu)勢(shì)環(huán)境功能
    也談詩(shī)的“功能”
    長(zhǎng)期鍛煉創(chuàng)造體內(nèi)抑癌環(huán)境
    一種用于自主學(xué)習(xí)的虛擬仿真環(huán)境
    矮的優(yōu)勢(shì)
    孕期遠(yuǎn)離容易致畸的環(huán)境
    環(huán)境
    關(guān)于非首都功能疏解的幾點(diǎn)思考
    畫(huà)與話
    談“五老”的五大特殊優(yōu)勢(shì)
    中西醫(yī)結(jié)合治療甲狀腺功能亢進(jìn)癥31例
    精品人妻一区二区三区麻豆 | 国产精品亚洲一级av第二区| 国产熟女xx| 国产亚洲精品综合一区在线观看| 精品福利观看| 国产亚洲精品av在线| 欧美xxxx性猛交bbbb| 亚洲最大成人av| 免费观看人在逋| 淫秽高清视频在线观看| 日韩欧美国产一区二区入口| 国产成人av教育| 老司机午夜十八禁免费视频| 午夜老司机福利剧场| 听说在线观看完整版免费高清| 精品久久久久久久久亚洲 | 精品一区二区三区视频在线观看免费| 欧美区成人在线视频| 人妻制服诱惑在线中文字幕| 九九在线视频观看精品| 精品人妻视频免费看| www日本黄色视频网| 色哟哟·www| 免费观看人在逋| 一本一本综合久久| 亚洲最大成人手机在线| 黄色视频,在线免费观看| 午夜福利免费观看在线| www.www免费av| 日韩欧美一区二区三区在线观看| 日本精品一区二区三区蜜桃| 国产主播在线观看一区二区| 五月玫瑰六月丁香| 久久久久久久久大av| 国产精品三级大全| 黄色视频,在线免费观看| 岛国在线免费视频观看| 99久久精品热视频| 欧美一区二区国产精品久久精品| 免费看美女性在线毛片视频| 99久久无色码亚洲精品果冻| 最后的刺客免费高清国语| 国产精华一区二区三区| 可以在线观看的亚洲视频| 国产亚洲欧美98| 久久精品国产99精品国产亚洲性色| 97超级碰碰碰精品色视频在线观看| 精品欧美国产一区二区三| 亚洲欧美日韩高清专用| 久久久色成人| 午夜激情福利司机影院| 欧美丝袜亚洲另类 | 国产精品嫩草影院av在线观看 | 岛国在线免费视频观看| 亚洲人成电影免费在线| 精品人妻熟女av久视频| 久久久久免费精品人妻一区二区| 婷婷精品国产亚洲av在线| 欧美性猛交黑人性爽| 精品人妻熟女av久视频| 欧美激情国产日韩精品一区| 热99在线观看视频| 午夜福利视频1000在线观看| 国产激情偷乱视频一区二区| 自拍偷自拍亚洲精品老妇| 少妇丰满av| 丁香六月欧美| 欧美激情久久久久久爽电影| 欧美激情久久久久久爽电影| 国产精品久久久久久久电影| 久久人人精品亚洲av| 亚洲 国产 在线| 欧美又色又爽又黄视频| 亚洲精品一区av在线观看| 一级a爱片免费观看的视频| 神马国产精品三级电影在线观看| 色综合亚洲欧美另类图片| 欧美黑人欧美精品刺激| 免费搜索国产男女视频| 国模一区二区三区四区视频| 久久久久久久久大av| 亚洲中文字幕一区二区三区有码在线看| 日本黄大片高清| 国产高潮美女av| 日韩精品青青久久久久久| 宅男免费午夜| 欧美成人性av电影在线观看| 中文字幕精品亚洲无线码一区| 99久久精品国产亚洲精品| 天堂网av新在线| 搡老熟女国产l中国老女人| 一进一出抽搐gif免费好疼| 亚洲成人精品中文字幕电影| 国产精品久久久久久精品电影| 久久婷婷人人爽人人干人人爱| 欧美一区二区国产精品久久精品| 给我免费播放毛片高清在线观看| 久久草成人影院| 亚洲成av人片在线播放无| 国产视频一区二区在线看| avwww免费| 在线观看一区二区三区| 十八禁网站免费在线| 国产精品爽爽va在线观看网站| 69人妻影院| 成人av在线播放网站| 欧美国产日韩亚洲一区| 亚洲va日本ⅴa欧美va伊人久久| 日本成人三级电影网站| avwww免费| 欧美zozozo另类| 18美女黄网站色大片免费观看| 午夜免费激情av| 日韩欧美精品免费久久 | 成人性生交大片免费视频hd| 97超级碰碰碰精品色视频在线观看| 韩国av一区二区三区四区| 成人毛片a级毛片在线播放| 丰满的人妻完整版| 乱人视频在线观看| 欧美黄色淫秽网站| 国产伦在线观看视频一区| 91av网一区二区| 亚洲av熟女| 成年免费大片在线观看| 中文字幕久久专区| 欧美三级亚洲精品| 99热只有精品国产| 精品久久久久久久末码| 欧美另类亚洲清纯唯美| 亚洲欧美日韩东京热| 搡老岳熟女国产| 亚洲av熟女| 国产麻豆成人av免费视频| 国产在视频线在精品| 在线天堂最新版资源| 国产蜜桃级精品一区二区三区| 最近中文字幕高清免费大全6 | 亚洲精品久久国产高清桃花| 中出人妻视频一区二区| 赤兔流量卡办理| 欧美一区二区精品小视频在线| 一二三四社区在线视频社区8| 麻豆成人午夜福利视频| av专区在线播放| 一进一出抽搐gif免费好疼| 国产黄a三级三级三级人| 国产麻豆成人av免费视频| 9191精品国产免费久久| 国产av在哪里看| 男女之事视频高清在线观看| 99热这里只有是精品在线观看 | 51午夜福利影视在线观看| 精品日产1卡2卡| 国产探花极品一区二区| 国产伦精品一区二区三区视频9| 日韩av在线大香蕉| aaaaa片日本免费| 中文字幕av在线有码专区| 在线观看av片永久免费下载| 哪里可以看免费的av片| 69av精品久久久久久| 欧美xxxx性猛交bbbb| 美女被艹到高潮喷水动态| 我的老师免费观看完整版| 国产一区二区亚洲精品在线观看| 精品欧美国产一区二区三| 国产精品一区二区免费欧美| 亚洲av二区三区四区| 久久精品国产自在天天线| 日本免费a在线| 亚洲av不卡在线观看| 国内精品一区二区在线观看| 国产免费av片在线观看野外av| 丝袜美腿在线中文| 国产aⅴ精品一区二区三区波| 国产一区二区三区视频了| 哪里可以看免费的av片| 国内精品久久久久精免费| 99久国产av精品| 怎么达到女性高潮| 久久欧美精品欧美久久欧美| 最近视频中文字幕2019在线8| 热99re8久久精品国产| 久久久久免费精品人妻一区二区| 最好的美女福利视频网| 日韩欧美免费精品| 欧美xxxx性猛交bbbb| 日日摸夜夜添夜夜添av毛片 | 色噜噜av男人的天堂激情| 亚洲性夜色夜夜综合| 国产精品久久久久久人妻精品电影| 日韩欧美国产在线观看| 精品一区二区三区人妻视频| 亚洲自偷自拍三级| 日韩 亚洲 欧美在线| 少妇的逼水好多| 久久久久亚洲av毛片大全| 哪里可以看免费的av片| 热99re8久久精品国产| 国产午夜精品论理片| 少妇人妻一区二区三区视频| 91久久精品国产一区二区成人| 国产精品久久电影中文字幕| 国产探花极品一区二区| 国产爱豆传媒在线观看| 床上黄色一级片| 中文字幕精品亚洲无线码一区| 久久人妻av系列| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩中文字幕国产精品一区二区三区| 很黄的视频免费| 日韩国内少妇激情av| 免费av毛片视频| 国产淫片久久久久久久久 | 色哟哟哟哟哟哟| 精品久久久久久久久av| 热99re8久久精品国产| 一本综合久久免费| 91av网一区二区| 久久午夜福利片| 嫩草影院精品99| 好男人在线观看高清免费视频| 亚洲中文日韩欧美视频| 国内少妇人妻偷人精品xxx网站| 舔av片在线| 99热这里只有是精品在线观看 | 岛国在线免费视频观看| 欧美中文日本在线观看视频| 韩国av一区二区三区四区| 午夜免费成人在线视频| 色哟哟·www| 国产成人欧美在线观看| 老熟妇乱子伦视频在线观看| 日本 av在线| 日韩高清综合在线| 亚洲不卡免费看| 国产精品一区二区性色av| 国产精品自产拍在线观看55亚洲| 色综合欧美亚洲国产小说| 在线观看免费视频日本深夜| www.色视频.com| 黄色视频,在线免费观看| 亚洲美女搞黄在线观看 | 亚洲不卡免费看| 1000部很黄的大片| 一夜夜www| 国产在视频线在精品| 欧美潮喷喷水| 人人妻人人澡欧美一区二区| 嫩草影院入口| 亚洲美女搞黄在线观看 | 伊人久久精品亚洲午夜| 久久精品91蜜桃| 免费av毛片视频| 1024手机看黄色片| 午夜视频国产福利| 小说图片视频综合网站| 精华霜和精华液先用哪个| 天天躁日日操中文字幕| 欧美激情在线99| 中文在线观看免费www的网站| 老司机深夜福利视频在线观看| 无遮挡黄片免费观看| 哪里可以看免费的av片| 可以在线观看的亚洲视频| 丁香六月欧美| 永久网站在线| 午夜福利欧美成人| 啪啪无遮挡十八禁网站| 久久久久性生活片| 如何舔出高潮| 麻豆成人av在线观看| 精品乱码久久久久久99久播| 国产亚洲精品久久久com| 99热这里只有是精品50| 99久久久亚洲精品蜜臀av| 长腿黑丝高跟| 国产成人影院久久av| 直男gayav资源| 一区二区三区激情视频| 久久久成人免费电影| 噜噜噜噜噜久久久久久91| 国产精品久久久久久人妻精品电影| 淫妇啪啪啪对白视频| 精品人妻1区二区| 最近在线观看免费完整版| 亚洲精品一卡2卡三卡4卡5卡| 国内毛片毛片毛片毛片毛片| 在线国产一区二区在线| 麻豆国产97在线/欧美| 国产又黄又爽又无遮挡在线| 国产主播在线观看一区二区| 99热这里只有精品一区| 一a级毛片在线观看| 在线看三级毛片| 91久久精品电影网| 1024手机看黄色片| 麻豆久久精品国产亚洲av| 午夜福利在线观看吧| 欧美精品啪啪一区二区三区| 最近中文字幕高清免费大全6 | 亚洲欧美日韩东京热| 最近视频中文字幕2019在线8| 成人一区二区视频在线观看| 久久精品影院6| 久久人妻av系列| 99久久精品热视频| 国产av一区在线观看免费| 老司机福利观看| 99国产极品粉嫩在线观看| 麻豆国产av国片精品| 午夜两性在线视频| 午夜老司机福利剧场| 国产乱人伦免费视频| 观看美女的网站| 9191精品国产免费久久| 国产精品自产拍在线观看55亚洲| 国产亚洲精品av在线| 久久精品综合一区二区三区| 丁香六月欧美| 亚洲无线在线观看| 看片在线看免费视频| 欧美黑人欧美精品刺激| 成人美女网站在线观看视频| 美女黄网站色视频| 国产成人a区在线观看| 宅男免费午夜| 国产又黄又爽又无遮挡在线| 久久久久久久久久黄片| 99热这里只有是精品50| 亚洲熟妇中文字幕五十中出| 日本一本二区三区精品| 欧美精品国产亚洲| 成人欧美大片| 成人特级av手机在线观看| 又黄又爽又刺激的免费视频.| eeuss影院久久| 中出人妻视频一区二区| 在线观看一区二区三区| 国产久久久一区二区三区| 国产不卡一卡二| 久久精品国产99精品国产亚洲性色| av中文乱码字幕在线| 麻豆成人午夜福利视频| 日韩欧美 国产精品| 久久九九热精品免费| 国产视频一区二区在线看| 欧美成人一区二区免费高清观看| 午夜两性在线视频| 日韩亚洲欧美综合| 国产爱豆传媒在线观看| 精品无人区乱码1区二区| 好男人电影高清在线观看| 精品熟女少妇八av免费久了| 欧美+亚洲+日韩+国产| bbb黄色大片| 欧美一区二区精品小视频在线| 99在线视频只有这里精品首页| 午夜福利在线观看免费完整高清在 | 日本 欧美在线| 少妇高潮的动态图| 亚洲精华国产精华精| 亚洲av五月六月丁香网| 色精品久久人妻99蜜桃| 亚洲国产高清在线一区二区三| 桃红色精品国产亚洲av| 色哟哟哟哟哟哟| 欧美色视频一区免费| 三级国产精品欧美在线观看| 精品一区二区三区视频在线| 亚洲精品亚洲一区二区| 夜夜夜夜夜久久久久| 国产野战对白在线观看| 91麻豆av在线| 久久这里只有精品中国| 一本精品99久久精品77| 国产精品免费一区二区三区在线| 国产在视频线在精品| 色吧在线观看| 精品国内亚洲2022精品成人| 90打野战视频偷拍视频| 一本综合久久免费| 欧美在线一区亚洲| 看十八女毛片水多多多| 男插女下体视频免费在线播放| 亚洲第一区二区三区不卡| 国产精品嫩草影院av在线观看 | 久久国产乱子免费精品| 中文字幕免费在线视频6| 一a级毛片在线观看| 亚洲精品乱码久久久v下载方式| 国产精品av视频在线免费观看| 黄色丝袜av网址大全| 高清日韩中文字幕在线| 脱女人内裤的视频| 欧美性猛交黑人性爽| 欧美日本视频| 少妇裸体淫交视频免费看高清| 搞女人的毛片| 欧美成人性av电影在线观看| 三级男女做爰猛烈吃奶摸视频| 最后的刺客免费高清国语| 亚洲第一区二区三区不卡| 国产欧美日韩一区二区精品| 国产精品影院久久| 欧美成人免费av一区二区三区| 午夜a级毛片| 国产av一区在线观看免费| 精品乱码久久久久久99久播| 国产av麻豆久久久久久久| 一个人免费在线观看的高清视频| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲国产精品成人综合色| 一级黄色大片毛片| 国产淫片久久久久久久久 | 男人和女人高潮做爰伦理| 国产高清视频在线播放一区| 夜夜夜夜夜久久久久| 日韩中字成人| 亚洲天堂国产精品一区在线| 国产国拍精品亚洲av在线观看| 小蜜桃在线观看免费完整版高清| 搡老妇女老女人老熟妇| 变态另类成人亚洲欧美熟女| 精品欧美国产一区二区三| 亚洲专区中文字幕在线| 国产精品伦人一区二区| 黄色女人牲交| 高潮久久久久久久久久久不卡| 欧美一区二区国产精品久久精品| 国产免费男女视频| 欧美成人一区二区免费高清观看| 欧美日韩黄片免| 欧美日韩福利视频一区二区| 欧美日本亚洲视频在线播放| 别揉我奶头 嗯啊视频| 久久久成人免费电影| 久久久久久久亚洲中文字幕 | 在线十欧美十亚洲十日本专区| 亚洲五月天丁香| 国产午夜福利久久久久久| 国产成人aa在线观看| 日本与韩国留学比较| 亚洲avbb在线观看| 日韩欧美三级三区| 在线国产一区二区在线| 91久久精品国产一区二区成人| 久久香蕉精品热| 精品熟女少妇八av免费久了| 亚洲国产精品999在线| 中文字幕免费在线视频6| 国产色爽女视频免费观看| 亚洲熟妇中文字幕五十中出| 9191精品国产免费久久| 亚洲av熟女| 国产精品久久久久久久电影| 18禁黄网站禁片午夜丰满| 亚洲成人精品中文字幕电影| 直男gayav资源| 波多野结衣高清作品| 午夜免费男女啪啪视频观看 | 一个人免费在线观看的高清视频| 亚洲欧美日韩东京热| 亚洲精品久久国产高清桃花| ponron亚洲| 美女高潮的动态| 亚洲美女视频黄频| 人妻丰满熟妇av一区二区三区| 亚洲最大成人av| 又黄又爽又免费观看的视频| 午夜免费激情av| 欧美+亚洲+日韩+国产| 亚洲欧美精品综合久久99| 91麻豆精品激情在线观看国产| av专区在线播放| 国产黄色小视频在线观看| 国产视频一区二区在线看| 国产精品美女特级片免费视频播放器| 丁香欧美五月| av在线老鸭窝| 亚洲最大成人av| 又爽又黄a免费视频| 精品久久久久久久久av| 日韩国内少妇激情av| 99riav亚洲国产免费| 一区二区三区高清视频在线| 成年女人毛片免费观看观看9| 他把我摸到了高潮在线观看| 国产在线男女| 色播亚洲综合网| 免费在线观看影片大全网站| 国产成人欧美在线观看| 国产精品免费一区二区三区在线| 国产精品人妻久久久久久| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产高清在线一区二区三| 精品无人区乱码1区二区| 午夜a级毛片| 欧美最黄视频在线播放免费| 天美传媒精品一区二区| 亚洲第一区二区三区不卡| 亚洲成a人片在线一区二区| 草草在线视频免费看| 可以在线观看毛片的网站| 狠狠狠狠99中文字幕| 999久久久精品免费观看国产| x7x7x7水蜜桃| 99久久精品热视频| 中文字幕高清在线视频| 91久久精品国产一区二区成人| 免费在线观看影片大全网站| 亚洲第一欧美日韩一区二区三区| 亚洲自偷自拍三级| 一个人看视频在线观看www免费| 男插女下体视频免费在线播放| 人人妻,人人澡人人爽秒播| 亚洲久久久久久中文字幕| 亚洲内射少妇av| 少妇熟女aⅴ在线视频| 国产亚洲av嫩草精品影院| 亚洲美女黄片视频| 日本一本二区三区精品| 日本免费a在线| 欧美三级亚洲精品| 一级av片app| 男女床上黄色一级片免费看| 内射极品少妇av片p| 国产亚洲精品久久久com| 无遮挡黄片免费观看| 亚洲人成网站高清观看| 精品免费久久久久久久清纯| 国产精品一区二区性色av| 亚洲欧美精品综合久久99| 免费高清视频大片| 婷婷亚洲欧美| 欧美xxxx性猛交bbbb| 男女做爰动态图高潮gif福利片| 真人做人爱边吃奶动态| 搡老妇女老女人老熟妇| 两个人的视频大全免费| 午夜福利成人在线免费观看| 亚洲国产欧美人成| 亚洲专区国产一区二区| 国内精品久久久久久久电影| 国产私拍福利视频在线观看| 欧美乱色亚洲激情| 亚洲av免费高清在线观看| 国产真实伦视频高清在线观看 | 国产爱豆传媒在线观看| 亚洲精品乱码久久久v下载方式| 丰满乱子伦码专区| 国产精品亚洲一级av第二区| 欧美黄色淫秽网站| 十八禁网站免费在线| 欧美激情久久久久久爽电影| 日韩 亚洲 欧美在线| www.色视频.com| 91狼人影院| 久久久色成人| 男人舔女人下体高潮全视频| 中文字幕高清在线视频| 久久久久九九精品影院| 丁香欧美五月| 看片在线看免费视频| 国模一区二区三区四区视频| 村上凉子中文字幕在线| 亚洲三级黄色毛片| 老鸭窝网址在线观看| 日本五十路高清| 国产色爽女视频免费观看| 91在线精品国自产拍蜜月| 国产高清视频在线播放一区| 国产精品永久免费网站| 偷拍熟女少妇极品色| 国产精品一区二区性色av| 国产精品久久电影中文字幕| 嫩草影院精品99| 永久网站在线| 亚洲男人的天堂狠狠| a级毛片a级免费在线| 亚洲成人免费电影在线观看| 91麻豆av在线| 九九热线精品视视频播放| 全区人妻精品视频| 极品教师在线视频| 哪里可以看免费的av片| 在线观看66精品国产| 性欧美人与动物交配| 久久人人爽人人爽人人片va | 在线观看美女被高潮喷水网站 | 国产精华一区二区三区| 少妇被粗大猛烈的视频| 午夜福利高清视频| 国产成人欧美在线观看| 高清日韩中文字幕在线| 波多野结衣巨乳人妻| 成人午夜高清在线视频| 国产极品精品免费视频能看的| 成年免费大片在线观看| 久久精品国产清高在天天线| av专区在线播放| 91久久精品国产一区二区成人| 免费在线观看影片大全网站| 亚洲精品456在线播放app | 极品教师在线免费播放| 少妇高潮的动态图| 给我免费播放毛片高清在线观看| 日韩欧美国产一区二区入口| 男人狂女人下面高潮的视频| 久久99热这里只有精品18| 99热精品在线国产| 男女之事视频高清在线观看| 日韩中文字幕欧美一区二区| 51国产日韩欧美|