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

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

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

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

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

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

    2013年3月、6月、9月及12月采集了洞庭湖11個斷面的浮游藻類,并對所采集的藻類進行了功能群劃分,結(jié)果顯示,洞庭湖的浮游藻類可以分為24個功能群: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)為對分層敏感、頻繁擾動且混合程度較高的渾濁型中-富營養(yǎng)型水體.不同時期調(diào)查的藻類優(yōu)勢功能群存在明顯差異,其中,B(中營養(yǎng)水體、對分層敏感)、D(較渾濁的淺水水體)、J(混合型高富營養(yǎng)淺水水體)、MP(擾動頻繁的渾濁型淺水水體)、P(混合程度較高中富營養(yǎng)淺水水體)、Y(廣適性)在四次調(diào)查中的優(yōu)勢度均>0.02,成為洞庭湖的絕對優(yōu)勢功能群,洞庭湖藻類優(yōu)勢功能群不同時期的演替規(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é)果顯示,洞庭湖藻類功能群分布受水環(huán)境因子影響較為明顯.整體上,T、pH值、DO、CODCr、NH3-N及TP是影響洞庭湖藻類功能群分布格局的主要因素.

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

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

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

    1 材料與方法

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

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

    1.2 水樣測試方法

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

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

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

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

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

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

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

    2 結(jié)果分析

    2.1 浮游藻類豐度變化及功能群劃分

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

    對所檢出的浮游藻類進行功能分類,如表 1所示,洞庭湖的浮游藻類可以分為 24個功能群: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年浮游藻類生物量的時空變化Fig.2 Temporal and spatial variation of phytoplankton of Dongting Lake in 2013

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

    續(xù)表1

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

    從不同時期藻類功能群的優(yōu)勢度來看(圖3),3月浮游藻類的優(yōu)勢功能群為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)勢度均>0.02,成為洞庭湖的絕對優(yōu)勢功能群.本研究中浮游藻類功能群的時空分布特征及其影響因子分析均在上述優(yōu)勢功能群中展開.

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

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

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

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

    為了使每個主成分的意義更加明確,采用Varimax最大方差法對因子進行旋轉(zhuǎn),按照特征值λ大于1的原則提取出3個主成分,如表2所示,3月、6月、9月和12月的主成分累計貢獻率分別達到81%、94.5%、78.5%及92.8%.該主成分分析所得指標的載荷超過 0.7的即為主要影響因子,3月的主要環(huán)境因子為T、pH值、DO、CODMn、CODCr、BOD、NH3-N、TP、E.coli,6月的主要環(huán)境因子為T、pH值、DO、CODMn、CODCr、BOD、NH3-N、TN、E.coli,9月的主要環(huán)境因子為T、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 浮游藻類優(yōu)勢功能群與環(huán)境因子的相關(guān) (CCA)分析

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

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

    3 討論

    3.1 浮游藻類豐度變化及功能群劃分

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

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

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

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

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

    4 結(jié)論

    4.1 2013年對洞庭湖 11個采樣斷面所采集的浮游藻類的分類結(jié)果顯示,洞庭湖的浮游藻類可以分為24個功能群: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 洞庭湖不同時期的藻類優(yōu)勢功能群存在明顯差異,其中,B、D、J、MP、P、Y在4次調(diào)查中的優(yōu)勢度均>0.02,成為洞庭湖的絕對優(yōu)勢功能群;洞庭湖藻類優(yōu)勢功能群不同季節(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é)果顯示,洞庭湖藻類功能群分布受水環(huán)境因子影響較為明顯.整體上,T、pH值、DO、CODCr、NH3-N及TP是影響洞庭湖藻類功能群分布格局的主要因素.

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

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

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

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

    [15] 李 興,李建茹,徐效清,等.烏梁素海浮游植物功能群季節(jié)演替規(guī)律及影響因子 [J]. 生態(tài)環(huán)境學(xué)報, 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] 汪 星,鄭丙輝,劉錄三,等.洞庭湖典型斷面藻類組成及其與環(huán)境因子典范對應(yīng)分析 [J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2012,31(5):995-1002.

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

    [19] 國家環(huán)??偩帧端蛷U水監(jiān)測分析方法》編委會.水和廢水監(jiān)測分析方法 [M]. 4版.北京:中國環(huán)境出版社, 2002.

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

    [21] 徐金森,鄭天凌,郭清華,等.兩種海洋細菌對赤潮藻的細胞生物量的影響研究 [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] 劉 瀟,薛 瑩,紀毓鵬,等.基于主成分分析法的黃河口及其鄰近水域水質(zhì)評價 [J]. 中國環(huán)境科學(xué), 2015,35(10):3187-3192.

    [25] 王 華,楊樹平,房晟忠,等.滇池浮游植物群落特征及與環(huán)境因子的典范對應(yīng)分析 [J]. 中國環(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] 李利強,張建波.洞庭湖浮游植物調(diào)查與水質(zhì)評價 [J]. 江蘇環(huán)境科技, 1999,12(4):14-16.

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

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

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

    [31] 周 泓,歐伏平,劉 妍.“十一五”期間洞庭湖水環(huán)境質(zhì)量狀況及變化趨勢分析 [J]. 湖南理工學(xué)院學(xué)報(自然科學(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] 汪 星,劉錄三,李 黎,等.鏡泊湖浮游藻類組成及其與環(huán)境因子的相關(guān)分析 [J]. 中國環(huán)境科學(xué), 2015,35(11):3403-3413.

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

    [36] 劉鴻雁,徐云麟.鏡泊湖藻類生長和湖泊富營養(yǎng)化預(yù)測初探 [J].生態(tài)學(xué)報, 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] 潘繼征,熊 飛,李文朝,等.云南撫仙湖透明度的時空變化及影響因子分析 [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] 楊麗標,韓小勇,孫 璞,等.巢湖藻類組成與環(huán)境因子典范對應(yīng)分析 [J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2011,30(5):952-958.

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

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

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

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

    [45] 陳曉江,楊 劼,杜桂森,等.海子水庫浮游植物功能群季節(jié)演替及其驅(qū)動因子 [J]. 水資源保護, 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-),男,湖北荊州人,副研究員,博士,主要從事水生生物評價的研究.發(fā)表論文30余篇.

    2016-06-02

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

    * 責任作者, 研究員, zhengbinghui@craes.org.cn

    猜你喜歡
    優(yōu)勢環(huán)境功能
    也談詩的“功能”
    中華詩詞(2022年6期)2022-12-31 06:41:24
    長期鍛煉創(chuàng)造體內(nèi)抑癌環(huán)境
    一種用于自主學(xué)習(xí)的虛擬仿真環(huán)境
    矮的優(yōu)勢
    趣味(語文)(2020年3期)2020-07-27 01:42:46
    孕期遠離容易致畸的環(huán)境
    環(huán)境
    關(guān)于非首都功能疏解的幾點思考
    畫與話
    談“五老”的五大特殊優(yōu)勢
    中國火炬(2014年11期)2014-07-25 10:31:58
    中西醫(yī)結(jié)合治療甲狀腺功能亢進癥31例
    日本一二三区视频观看| 国产精品电影一区二区三区| 久久精品91无色码中文字幕| 一二三四社区在线视频社区8| 精品久久久久久,| 亚洲精品一卡2卡三卡4卡5卡| 精品国内亚洲2022精品成人| 精品99又大又爽又粗少妇毛片 | 一本精品99久久精品77| 欧美中文日本在线观看视频| 一级黄片播放器| 美女cb高潮喷水在线观看| av视频在线观看入口| 有码 亚洲区| 91麻豆精品激情在线观看国产| 国产午夜精品久久久久久一区二区三区 | 最新中文字幕久久久久| 大型黄色视频在线免费观看| 波野结衣二区三区在线 | 精品电影一区二区在线| 久久久久国产精品人妻aⅴ院| 国产激情欧美一区二区| 精品日产1卡2卡| 成人av在线播放网站| а√天堂www在线а√下载| 久久久久久久午夜电影| 美女黄网站色视频| 亚洲精品成人久久久久久| 亚洲人成网站在线播放欧美日韩| 欧美高清成人免费视频www| 91九色精品人成在线观看| 国内精品久久久久精免费| 色综合欧美亚洲国产小说| 成人鲁丝片一二三区免费| 最后的刺客免费高清国语| 在线观看美女被高潮喷水网站 | 久久人人精品亚洲av| 99国产极品粉嫩在线观看| 在线观看av片永久免费下载| 亚洲专区国产一区二区| 在线观看舔阴道视频| 国产精品一区二区免费欧美| 成年女人毛片免费观看观看9| 久99久视频精品免费| 一级作爱视频免费观看| 欧美在线一区亚洲| 精品久久久久久久毛片微露脸| 长腿黑丝高跟| 天堂√8在线中文| 五月玫瑰六月丁香| www.999成人在线观看| 国产野战对白在线观看| 伊人久久大香线蕉亚洲五| 日本黄大片高清| av天堂中文字幕网| 精品一区二区三区人妻视频| 欧美三级亚洲精品| 一级黄色大片毛片| 搡老熟女国产l中国老女人| 亚洲国产日韩欧美精品在线观看 | 色噜噜av男人的天堂激情| 欧洲精品卡2卡3卡4卡5卡区| 日韩有码中文字幕| 一本一本综合久久| 日日干狠狠操夜夜爽| 九九热线精品视视频播放| 欧美黄色片欧美黄色片| 国产精品免费一区二区三区在线| 亚洲美女黄片视频| 国产乱人视频| 波野结衣二区三区在线 | 我的老师免费观看完整版| 免费无遮挡裸体视频| 亚洲成av人片在线播放无| 免费观看精品视频网站| 成人18禁在线播放| 色视频www国产| 毛片女人毛片| 嫩草影院入口| 偷拍熟女少妇极品色| 久久精品影院6| 草草在线视频免费看| 久久国产乱子伦精品免费另类| 国产99白浆流出| 老汉色∧v一级毛片| 网址你懂的国产日韩在线| 午夜免费男女啪啪视频观看 | 欧美另类亚洲清纯唯美| 欧美大码av| 麻豆一二三区av精品| 欧美日本视频| 久久久久国产精品人妻aⅴ院| 黄色成人免费大全| 久久精品国产综合久久久| 国产在线精品亚洲第一网站| 男女那种视频在线观看| 精品久久久久久久久久免费视频| 亚洲天堂国产精品一区在线| 国产高清videossex| 蜜桃久久精品国产亚洲av| 精品一区二区三区视频在线观看免费| 欧美日韩亚洲国产一区二区在线观看| 亚洲av日韩精品久久久久久密| 在线免费观看的www视频| 国产不卡一卡二| 国产真人三级小视频在线观看| 18禁国产床啪视频网站| 床上黄色一级片| 天天躁日日操中文字幕| 亚洲精华国产精华精| 免费看十八禁软件| 欧美zozozo另类| 黄片大片在线免费观看| 国产精品美女特级片免费视频播放器| av片东京热男人的天堂| 国产精品久久电影中文字幕| 中文字幕人妻熟人妻熟丝袜美 | 老汉色∧v一级毛片| 又粗又爽又猛毛片免费看| 久久精品国产自在天天线| 91字幕亚洲| 身体一侧抽搐| 午夜福利高清视频| 男女午夜视频在线观看| 亚洲片人在线观看| 欧美黄色片欧美黄色片| 久久精品综合一区二区三区| 欧美日韩乱码在线| 天天躁日日操中文字幕| 狂野欧美白嫩少妇大欣赏| 国产单亲对白刺激| 国内久久婷婷六月综合欲色啪| 免费一级毛片在线播放高清视频| 网址你懂的国产日韩在线| 看片在线看免费视频| 国产色婷婷99| 内射极品少妇av片p| 成年女人看的毛片在线观看| 最好的美女福利视频网| 在线观看一区二区三区| 国产av不卡久久| 叶爱在线成人免费视频播放| 一个人观看的视频www高清免费观看| 国产亚洲av嫩草精品影院| 午夜激情欧美在线| 亚洲第一欧美日韩一区二区三区| 91麻豆av在线| 亚洲不卡免费看| 精品人妻一区二区三区麻豆 | 丝袜美腿在线中文| 日韩av在线大香蕉| 久久亚洲真实| 亚洲精品影视一区二区三区av| 日韩有码中文字幕| 欧美国产日韩亚洲一区| 99久久精品一区二区三区| 国语自产精品视频在线第100页| 国产精品亚洲一级av第二区| 波多野结衣高清无吗| 精品久久久久久,| 热99在线观看视频| 国产国拍精品亚洲av在线观看 | 老司机午夜十八禁免费视频| 亚洲av不卡在线观看| 中文字幕人妻熟人妻熟丝袜美 | 90打野战视频偷拍视频| 久久6这里有精品| 午夜福利在线观看吧| 精品熟女少妇八av免费久了| 美女 人体艺术 gogo| 特大巨黑吊av在线直播| 夜夜躁狠狠躁天天躁| 国产亚洲精品综合一区在线观看| 小蜜桃在线观看免费完整版高清| 国产精品98久久久久久宅男小说| 亚洲精品在线美女| 搡老熟女国产l中国老女人| 手机成人av网站| 一进一出抽搐gif免费好疼| 欧美一区二区精品小视频在线| 成人国产一区最新在线观看| 国产精品久久久久久亚洲av鲁大| www日本黄色视频网| 日本精品一区二区三区蜜桃| АⅤ资源中文在线天堂| 欧美一区二区国产精品久久精品| 亚洲性夜色夜夜综合| 亚洲欧美日韩卡通动漫| 久久久久国内视频| 成年人黄色毛片网站| 无人区码免费观看不卡| 色视频www国产| 国产精品 欧美亚洲| 天堂网av新在线| 小蜜桃在线观看免费完整版高清| 熟妇人妻久久中文字幕3abv| 久久精品国产亚洲av涩爱 | 一边摸一边抽搐一进一小说| 岛国在线免费视频观看| 嫁个100分男人电影在线观看| 亚洲人成伊人成综合网2020| 亚洲人成网站高清观看| 两性午夜刺激爽爽歪歪视频在线观看| 色哟哟哟哟哟哟| 最新中文字幕久久久久| 19禁男女啪啪无遮挡网站| 青草久久国产| 天堂网av新在线| 国内精品美女久久久久久| 亚洲美女视频黄频| 51午夜福利影视在线观看| 日韩有码中文字幕| 日本三级黄在线观看| 免费av观看视频| 老熟妇乱子伦视频在线观看| 日韩人妻高清精品专区| www日本黄色视频网| 国产高清视频在线观看网站| 天美传媒精品一区二区| 久久午夜亚洲精品久久| 特大巨黑吊av在线直播| av专区在线播放| 在线观看66精品国产| 精品人妻偷拍中文字幕| 搡女人真爽免费视频火全软件 | 欧美又色又爽又黄视频| 久久久久久久久久黄片| 天天躁日日操中文字幕| 色综合欧美亚洲国产小说| 午夜日韩欧美国产| 日韩欧美国产一区二区入口| 国产一区二区在线av高清观看| 国产精品久久视频播放| 少妇的丰满在线观看| 一区二区三区国产精品乱码| 免费观看人在逋| 岛国在线免费视频观看| 欧美3d第一页| 在线国产一区二区在线| 此物有八面人人有两片| 十八禁网站免费在线| 十八禁人妻一区二区| 日本在线视频免费播放| 欧美激情久久久久久爽电影| 久久久久亚洲av毛片大全| 亚洲内射少妇av| 亚洲欧美日韩高清专用| 两个人看的免费小视频| 欧美激情久久久久久爽电影| 人人妻人人澡欧美一区二区| 日韩欧美三级三区| 国模一区二区三区四区视频| svipshipincom国产片| 黄色日韩在线| 成人三级黄色视频| 欧美性猛交╳xxx乱大交人| 黄色片一级片一级黄色片| 免费观看精品视频网站| 成人av一区二区三区在线看| 国产国拍精品亚洲av在线观看 | 中文在线观看免费www的网站| 亚洲人成网站在线播放欧美日韩| 热99在线观看视频| 最近最新中文字幕大全免费视频| 亚洲成人久久爱视频| www.熟女人妻精品国产| 国产亚洲精品一区二区www| 全区人妻精品视频| 美女大奶头视频| 国产一区二区在线av高清观看| 18禁裸乳无遮挡免费网站照片| 操出白浆在线播放| 国产综合懂色| 成人精品一区二区免费| 亚洲中文字幕一区二区三区有码在线看| 精品久久久久久,| 亚洲熟妇熟女久久| 男插女下体视频免费在线播放| 国产成人aa在线观看| 久久草成人影院| 天堂av国产一区二区熟女人妻| 国产伦精品一区二区三区四那| 中文资源天堂在线| 国产精品99久久99久久久不卡| 嫩草影视91久久| 美女黄网站色视频| 一个人免费在线观看的高清视频| 18禁裸乳无遮挡免费网站照片| 国产麻豆成人av免费视频| 亚洲人成网站高清观看| 757午夜福利合集在线观看| 99国产极品粉嫩在线观看| 他把我摸到了高潮在线观看| 欧美成人一区二区免费高清观看| 精品午夜福利视频在线观看一区| 天堂av国产一区二区熟女人妻| 黄色日韩在线| 亚洲成人中文字幕在线播放| 麻豆国产97在线/欧美| 一区二区三区免费毛片| 国产精品 欧美亚洲| 91九色精品人成在线观看| 99久久久亚洲精品蜜臀av| 亚洲精品国产精品久久久不卡| 很黄的视频免费| 成人一区二区视频在线观看| 精品国产亚洲在线| 十八禁网站免费在线| 成人av在线播放网站| 国产av麻豆久久久久久久| 日韩欧美 国产精品| 国产亚洲欧美在线一区二区| 黄色片一级片一级黄色片| 99久久精品一区二区三区| 欧美又色又爽又黄视频| 淫秽高清视频在线观看| 免费在线观看亚洲国产| 99久久精品一区二区三区| 亚洲电影在线观看av| 国产欧美日韩一区二区三| 欧美在线一区亚洲| 婷婷丁香在线五月| 国产伦在线观看视频一区| 亚洲美女黄片视频| 国产精品久久久久久人妻精品电影| 亚洲不卡免费看| 色哟哟哟哟哟哟| 国产精品一及| 精品欧美国产一区二区三| 国产探花在线观看一区二区| 成人av一区二区三区在线看| 国产一区二区激情短视频| 亚洲人成电影免费在线| 亚洲国产欧美网| 国产精品香港三级国产av潘金莲| 美女cb高潮喷水在线观看| 欧美一级a爱片免费观看看| 午夜福利18| 欧美成狂野欧美在线观看| 欧美日本视频| 中文字幕人成人乱码亚洲影| 99久久无色码亚洲精品果冻| 国产色爽女视频免费观看| 午夜福利高清视频| 日本a在线网址| 国产成人a区在线观看| 国产乱人视频| 天美传媒精品一区二区| 国产精品永久免费网站| 国产精品久久久久久久电影 | 亚洲第一电影网av| 日本黄大片高清| 精品久久久久久久毛片微露脸| 特大巨黑吊av在线直播| 成年女人看的毛片在线观看| 午夜福利免费观看在线| 久久久国产成人精品二区| 国产精品久久久久久久电影 | 久久久久久久亚洲中文字幕 | 亚洲乱码一区二区免费版| 亚洲欧美日韩卡通动漫| 夜夜躁狠狠躁天天躁| 91九色精品人成在线观看| 国产真人三级小视频在线观看| 18禁裸乳无遮挡免费网站照片| 99热精品在线国产| 亚洲精品亚洲一区二区| 午夜免费观看网址| 久久久色成人| 中文资源天堂在线| 亚洲美女视频黄频| 精品一区二区三区视频在线 | 国产av不卡久久| av片东京热男人的天堂| 欧美精品啪啪一区二区三区| 久久中文看片网| 亚洲人成伊人成综合网2020| 99国产精品一区二区三区| 美女黄网站色视频| 非洲黑人性xxxx精品又粗又长| 欧美日韩精品网址| 九色成人免费人妻av| 亚洲欧美精品综合久久99| 一级黄色大片毛片| 19禁男女啪啪无遮挡网站| 色综合欧美亚洲国产小说| 不卡一级毛片| 丝袜美腿在线中文| 丁香六月欧美| 男女下面进入的视频免费午夜| 可以在线观看的亚洲视频| 久久精品国产亚洲av涩爱 | 欧美最黄视频在线播放免费| 五月伊人婷婷丁香| 三级毛片av免费| 网址你懂的国产日韩在线| 九九热线精品视视频播放| 日本三级黄在线观看| 不卡一级毛片| 听说在线观看完整版免费高清| 国产成人啪精品午夜网站| 精品电影一区二区在线| 两个人视频免费观看高清| 99热这里只有是精品50| 男女午夜视频在线观看| 欧美xxxx黑人xx丫x性爽| 一区二区三区国产精品乱码| 日韩 欧美 亚洲 中文字幕| 亚洲欧美日韩高清专用| 欧美乱色亚洲激情| 97超视频在线观看视频| 亚洲成a人片在线一区二区| 网址你懂的国产日韩在线| 一个人免费在线观看的高清视频| 亚洲第一电影网av| av天堂在线播放| 深爱激情五月婷婷| 观看免费一级毛片| 夜夜看夜夜爽夜夜摸| 嫩草影院精品99| 啦啦啦观看免费观看视频高清| 国产三级在线视频| 十八禁人妻一区二区| 少妇的逼水好多| 国产主播在线观看一区二区| 国产美女午夜福利| 国产精品1区2区在线观看.| 国产精品一区二区三区四区久久| 夜夜夜夜夜久久久久| 国产乱人伦免费视频| 丰满乱子伦码专区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 最近最新中文字幕大全免费视频| 精品久久久久久久久久免费视频| 一个人免费在线观看电影| 在线看三级毛片| 日本成人三级电影网站| 最近视频中文字幕2019在线8| 制服人妻中文乱码| 美女大奶头视频| 毛片女人毛片| 国产精品亚洲美女久久久| 亚洲熟妇熟女久久| 九色国产91popny在线| 精品久久久久久,| 桃红色精品国产亚洲av| 少妇高潮的动态图| 少妇的逼好多水| 黑人欧美特级aaaaaa片| 少妇人妻一区二区三区视频| 国产亚洲精品久久久com| 色老头精品视频在线观看| 亚洲av成人av| 亚洲不卡免费看| 亚洲乱码一区二区免费版| 国产精华一区二区三区| 日韩国内少妇激情av| 亚洲精品美女久久久久99蜜臀| 国产精品久久久久久久久免 | 精品99又大又爽又粗少妇毛片 | 热99re8久久精品国产| 国产精品国产高清国产av| 天堂动漫精品| 老司机福利观看| 国产男靠女视频免费网站| 一本一本综合久久| 波野结衣二区三区在线 | 中文字幕精品亚洲无线码一区| 在线a可以看的网站| 18禁美女被吸乳视频| 久久久精品大字幕| eeuss影院久久| 午夜精品久久久久久毛片777| 伊人久久精品亚洲午夜| 国产精品一区二区三区四区久久| 精品日产1卡2卡| 成人特级黄色片久久久久久久| 天美传媒精品一区二区| 日日摸夜夜添夜夜添小说| 日韩欧美一区二区三区在线观看| 国产午夜精品论理片| 亚洲精品日韩av片在线观看 | 国产精品久久电影中文字幕| 伊人久久大香线蕉亚洲五| 在线观看午夜福利视频| 国产亚洲精品一区二区www| 国产一区二区激情短视频| 亚洲av成人精品一区久久| 精品乱码久久久久久99久播| 一卡2卡三卡四卡精品乱码亚洲| 男女午夜视频在线观看| 久久久色成人| 欧美最新免费一区二区三区 | 男女午夜视频在线观看| 麻豆成人午夜福利视频| 全区人妻精品视频| 夜夜爽天天搞| 欧美乱色亚洲激情| 中文字幕人妻丝袜一区二区| 成人18禁在线播放| 啦啦啦免费观看视频1| 热99re8久久精品国产| 久久香蕉精品热| 女人十人毛片免费观看3o分钟| 国产av麻豆久久久久久久| 美女高潮喷水抽搐中文字幕| 丰满人妻熟妇乱又伦精品不卡| 午夜影院日韩av| 给我免费播放毛片高清在线观看| 一个人看的www免费观看视频| av视频在线观看入口| 好看av亚洲va欧美ⅴa在| 国产一区二区激情短视频| 亚洲人成网站高清观看| 欧美bdsm另类| 性欧美人与动物交配| 又粗又爽又猛毛片免费看| 亚洲国产色片| 精品国产亚洲在线| 夜夜爽天天搞| 少妇高潮的动态图| 日本a在线网址| 亚洲无线观看免费| 极品教师在线免费播放| 人人妻人人看人人澡| 欧美一区二区国产精品久久精品| 99热这里只有是精品50| 九九久久精品国产亚洲av麻豆| 成年女人毛片免费观看观看9| 波多野结衣高清无吗| 最近在线观看免费完整版| 国产亚洲精品久久久久久毛片| 中亚洲国语对白在线视频| 亚洲最大成人中文| 日本免费a在线| 在线观看舔阴道视频| 精品久久久久久久毛片微露脸| 可以在线观看毛片的网站| 99精品欧美一区二区三区四区| 日本三级黄在线观看| 在线观看一区二区三区| 老司机午夜福利在线观看视频| 国产野战对白在线观看| 久久精品国产99精品国产亚洲性色| 欧美在线黄色| 99热这里只有精品一区| 精华霜和精华液先用哪个| aaaaa片日本免费| av天堂在线播放| 亚洲成a人片在线一区二区| 日本五十路高清| 女人被狂操c到高潮| 午夜福利18| 首页视频小说图片口味搜索| 最好的美女福利视频网| 又爽又黄无遮挡网站| 美女免费视频网站| 国产单亲对白刺激| 搡老岳熟女国产| 国产精品久久久久久久久免 | a级一级毛片免费在线观看| 长腿黑丝高跟| 久久人妻av系列| 男女做爰动态图高潮gif福利片| 色综合欧美亚洲国产小说| 日韩成人在线观看一区二区三区| 最新在线观看一区二区三区| 午夜免费成人在线视频| 日韩有码中文字幕| 好男人电影高清在线观看| 操出白浆在线播放| 91字幕亚洲| 国产v大片淫在线免费观看| av天堂在线播放| 最新在线观看一区二区三区| 午夜免费成人在线视频| 首页视频小说图片口味搜索| 欧美av亚洲av综合av国产av| 久久草成人影院| 美女免费视频网站| 91麻豆av在线| 一a级毛片在线观看| 国产精品久久电影中文字幕| 男人舔女人下体高潮全视频| av女优亚洲男人天堂| 91字幕亚洲| 成人av在线播放网站| 亚洲专区中文字幕在线| 国产中年淑女户外野战色| 人人妻,人人澡人人爽秒播| 桃色一区二区三区在线观看| 欧美av亚洲av综合av国产av| 婷婷丁香在线五月| 午夜a级毛片| 成人18禁在线播放| 久久中文看片网| 亚洲在线观看片| 亚洲国产欧美人成| 国产亚洲精品综合一区在线观看| 国产精品女同一区二区软件 | 欧美性感艳星| 9191精品国产免费久久| 免费观看人在逋| 激情在线观看视频在线高清| 一区二区三区免费毛片| 精品一区二区三区视频在线观看免费| 成年人黄色毛片网站| 我的老师免费观看完整版| 亚洲人与动物交配视频| 2021天堂中文幕一二区在线观| 首页视频小说图片口味搜索| 久久亚洲精品不卡| 夜夜躁狠狠躁天天躁| 俄罗斯特黄特色一大片| 欧美+日韩+精品|