999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

陽澄西湖浮游藻類功能群演替特征及其與環境因子的關系

2019-07-31 06:54:42周彥鋒王晨赫
中國環境科學 2019年7期
關鍵詞:優勢環境功能

葛 優,周彥鋒,王晨赫,尤 洋,*

陽澄西湖浮游藻類功能群演替特征及其與環境因子的關系

葛 優1,周彥鋒1,王晨赫2,尤 洋1,2*

(1.中國水產科學研究院淡水漁業研究中心,農業部淡水漁業和種質資源利用重點實驗室,江蘇 無錫 214081;2.南京農業大學無錫漁業學院,江蘇 無錫 214081)

采用浮游藻類功能群、 RDA(冗余分析)、多元回歸分析等方法,對陽澄西湖浮游藻類功能群的演替特征及其影響因子進行了研究.2015年3月~2016年2月逐月對陽澄西湖浮游藻類和環境因子進行調查分析,結果表明,陽澄西湖浮游藻類可劃分為21個功能群,其中10個功能群(B、D、E、G、J、Lo、S1、X1、X2、Y)為該湖的優勢功能群,反映的生境特征表明該湖為分層敏感,擾動較少且混合程度較高的中-富營養淺水水體.陽澄西湖不同季節優勢功能群演替規律為:由春季的B+X2+Y+X1+J+D,經夏季的B+X2+Y+X1+J+G+S1+Lo+E和秋季的B+X2+Y+Lo+E,轉變為B+X2+Y+X1+D.RDA分析結果顯示,陽澄西湖浮游藻類功能群分布受環境影響較為明顯,但不同季節影響因子不同.整體上,水溫()、pH值、溶解氧(DO)、氨氮(NH4+-N)和濁度(Tur)是影響陽澄西湖浮游藻類功能群分布格局的主要因素.

陽澄西湖;藻類功能群;優勢度;環境因子

浮游藻類作為水生生態系統的重要組成部分,是初級生產力的主要貢獻者[1],其種類組成、群落結構和生物量變化既能反映水環境的變化[2],又能在很大程度上決定著水生生態系統的結構和功能[3],是反映水體富營養化程度和水質狀況的重要指標.因此,在水生態研究中此類生物常被用作評估水體健康狀況的指示性生物類群[4].傳統的浮游植物群落研究基于系統分類方法,將種屬作為研究的基本單元[5],該種分類方法難以凸顯浮游植物在水生態系統中的生態作用和環境特征[6-8].此外,以功能群作為浮游藻類群落劃分方法[9],將具有相似生理、形態和生態特征的浮游植物種類劃分為一個功能群.經修正[3],將浮游植物劃分為39個功能群.功能群內的浮游植物具有相同或相近的生態位,不同的功能群可以反映其特定的生境類型.浮游植物功能類群以個體生態學為依據進行分類,能夠更加準確地描述浮游植物的生境特征,并與浮游植物群落演替過程有效結合,最大限度地揭示浮游植物群落與水體環境間的動態關系,較好地預測群落生態功能的變化[10-11].

國外學者對浮游植物功能群已有了全面的研究,被廣泛應用于湖泊、水庫、河流等多種類型水體的水生態系統研究中[12-15].近年來,國內學者對不同水域的浮游植物功能群也展開了研究,主要從不同類型水體的浮游植物功能群結構特征、季節演替和水質評價等發面進行研究[16-19].國內浮游植物功能群的研究主要集中在河流、水庫、濕地和高原湖泊.而浮游植物功能群的研究在長江下游、太湖流域淺水湖泊的報道相對較少.

陽澄湖為太湖流域第三大淺水湖泊.它不僅是蘇州地區的綜合水源地,還具有工業用水、漁業養殖、農業灌溉、旅游、航運和防汛等多種功能[20-21].由于長期圍網養殖的過度開發,導致陽澄湖水環境遭到了破壞.尤其是西湖區,水草已是零星可見,湖泊由原來的草型湖泊退化成藻型湖泊[22].有關陽澄湖浮游植物的研究主要集中在浮游植物的種類組成、群落結構、季節變化等方面[23-25],對于浮游植物功能群的演替特征及其影響因子的研究未見報道.本研究基于2015~2016年間陽澄西湖浮游藻類和水質理化指標的逐月監測數據,運用浮游藻類功能群分類方法,對其功能群進行劃分,探究浮游藻類功能群的演替特征,并通過冗余分析(RDA)和多元回歸分析明確影響浮游藻類功能群演替的主要環境因子,以期為陽澄西湖的生態研究和治理提供理論依據和技術支撐.

1 材料與方法

1.1 研究區域概況

陽澄湖地處江蘇省蘇州市境內,環湖為昆山市、蘇州市相城區和蘇州工業園區.湖中有2條東北-西南走向的狹長半島,把陽澄湖分為東湖(44.08%)、中湖(29.03%)和西湖(26.89%)整個湖體分為3部分[23].研究區域為蘇州市陽澄湖西湖,湖區水位變幅平緩,年變幅一般在1.2m以下,平均水位3.2m[26].年平均降雨量1100mm,汛期降雨多集中在6~9月.該湖所處區域屬北亞熱帶季風氣候,氣候溫暖濕潤,四季分明(春季 3~5月,夏季6~9月,秋季10~11月,冬季12月~次年2月).

1.2 樣品采集

依據陽澄西湖的地理形態和水文條件,對陽澄西湖4個斷面進行了調查,并結合樣點設置規范[27],選取12個采樣點(31°23¢26.8″~31°28¢55.1″N, 120°42¢13.3″~120°45¢32.8″E),編號采樣點分別設置為斷面Ⅰ(S1、S2、S3)、斷面Ⅱ(S4、S5、S6)、斷面Ⅲ(S7、S8、S9)和斷面Ⅳ(S10、S11、S12)(圖1).2015年 3月~2016年2月,逐月對陽澄西湖的12個采樣點進行樣品采集.

圖1 采樣點示意

1.3 樣品分析方法

水質樣品采集方法按照《湖泊富營養化調查規范》進行[28],水體理化監測項目為溫度()、pH值、溶解氧(DO)、濁度(Tur)、透明度(SD)、水深()、總磷(TP)、溶解性總磷(DTP)、總氮(TN)、氨氮(NH4+-N)、亞硝酸鹽氮(NO2--N)、高錳酸鹽指數(CODMn)、葉綠素a(Chl-a)共13項.其中、pH值、DO、Tur、SD和WD采用多功能水質檢測儀進行現場檢測;TP和DTP采用鉬酸銨分光光度法;TN和DTN采用堿性過硫酸鉀消解紫外分光光度法;NH4+-N采用納氏試劑分光光度法;NO2--N采用紫外分光光度法; CODMn采用酸式滴定法; Chl-a采用丙酮提取分光光度法測定.水質數據按照地表水環境質量標準(GB3838-2002)來進行分析[29].

浮游藻類樣品的采集及處理方法參照《淡水浮游生物研究方法》[30],使用5L有機玻璃采水器采集0.5m和1.5m的混合水樣,混勻后,取1L加入福爾馬林和魯哥氏碘液進行固定,并于實驗室靜置48h后,虹吸上清液濃縮至50mL.藻類鑒定主要參考《淡水微型生物與底棲動物圖譜(第二版)》[31]和《中國淡水藻類-系統、分類及生態》[32].

浮游藻類在10×40倍顯微鏡下采用目鏡視野法觀察,定性時優勢種應鑒定到種,其他種類鑒定到屬.定量樣品記錄100~200個視野,記下每種浮游藻類的細胞個數.每升水中浮游藻類的數量應按公式(1)進行計算:

式中:s為計數框面積,mm2;s= 20mm×20mm = 400mm2;s為視野面積,mm2;s= π/42=0.159mm2,=450μm;F為計數視野數;為濃縮體積,mL;為計數框容積,m L;P為計數,個.

1.4 數據分析

1.4.1 浮游藻類功能群優勢度分析 根據功能群分類方法[2-3]對淀山湖浮游藻類進行功能分組.優勢度根據藻類功能群的出現頻率及豐度來確定,用優勢度來表示[33]:

y

=

fi

·

Pi

(2)

式中:為優勢度;f為第功能群的出現頻率;P為第功能群豐度占總藻類豐度的比例,當>0.02時,定為優勢功能群[34].選擇優勢度>0.02的浮游藻類功能群進行冗余分析(RDA)分析.

1.4.2 冗余分析(RDA) 采用Canoco for Windows 4.5軟件對功能群數據和環境數據進行RDA分析.浮游藻類功能群矩陣經過lg(+1)轉換,環境因子數據除pH值外全部進行lg(+1)轉換[35].首先進行去趨勢對應分析(DCA)選擇合適的排序和分析方法,根據分析結果,本文選用了冗余分析(RDA),并采用蒙特卡洛置換檢驗,以逐步迭代的方式篩選出有顯著解釋性的環境因子,進行分析作圖,以揭示顯著環境因子對優勢浮游藻類功能群的影響.

本文其他分析和作圖在Origin 2017、SPSS22.0、CANOCO4.5和Excel中完成.

2 結果與分析

2.1 浮游藻類生物量變化及功能群劃分

從不同月份浮游藻類的分布來看(圖2),陽澄西湖浮游藻類季節分布存在極顯著差異(<0.01).3月藻類生物量變化范圍為0.123~1.158mg/L;4月變化范圍為0.260~15.744mg/L,最大值出現在3號點,最小值出現在12號點;5月藻類生物量變化范圍為0.187~3.253mg/L,最大值出現在2號點,最小值出現在7號點;浮游藻類生物量整體表現為南部>中部>北部.6月藻類生物量變化范圍為1.679~12.655mg/L,最大值出現在9號點,最小值出現在12號點;7月藻類生物量變化范圍為1.450~10.892mg/L,最大值出現在1號點,最小值出現在8號點;8月藻類生物量變化范圍為0.373~4.569mg/L,最大值出現在2號點,最小值出現在9號點;浮游藻類生物量整體表現為南部>北部>中部.9月藻類生物量變化范圍為0.354~3.380mg/L,最大值出現在10號點,最小值出現在5號點;10月藻類生物量變化范圍為0.050~ 1.630mg/L,最大值出現在11號點,最小值出現在6號點;11月藻類生物量變化范圍為0.259~1.194mg/L,最大值出現在2號點,最小值出現在12號點;浮游藻類生物量整體表現為北部>中部>南部.12月藻類生物量變化范圍為0.061-1.753mg/L,最大值出現在3號點,最小值出現在7號點;次年1月藻類生物量變化范圍為0.091~0.595mg/L,最大值出現在1號點,最小值出現在4號點;次年2月藻類生物量變化范圍為0.737~ 6.141mg/L,最大值出現在9號點,最小值出現在7號點;浮游藻類生物量整體表現為南部>中部>北部.

表1 陽澄西湖浮游藻類功能群劃分

續表1

根據功能群分類法[19-20],將具備相同或類似生態屬性的浮游藻類歸在一個功能群內,陽澄西湖浮游藻類可分為21個功能群:A、B、D、E、F、G、H1、J、Lo、M、MP、P、S1、T、TD、W1、W2、X1、X2、X3、Y(表1).

2.2 浮游藻類優勢功能群及分布特征

本次調查中優勢度>0.02的功能群處于優勢地位,為陽澄西湖優勢功能群.從不同時期浮游藻類功能群的優勢度來看(圖3),浮游藻類的優勢功能群3月為 B、Y、X2,4月為B、X2、Y、X1、D,5月為X2、Y、B、J、D,6月為Y、B、X2、J、G,7月為X2、B、G、S1、Y、Lo,8月為 X2、B、E、Lo、Y,9月為X2、Y、Lo、B,10月為X2、B、Y、E,11月為B、Y、X2,12月為B、X2、Y、X1,1月為 B、X2、Y、X1、D,2月為B、X2、Y.其中,B、X2、Y在12次的調查中的優勢度均30.02,為陽澄西湖的絕對優勢功能群.本研究中浮游藻類功能群的時空分布特征及其影響因子分析均在上述優勢功能群中展開.

從不同時期浮游藻類的功能群空間組成來看(圖4),3月浮游藻類主要功能群Y+X2在各采樣點均有分布,其中,Y在S1、S12占據絕對優勢;X2在S1占據絕對優勢.4月浮游藻類主要功能群由3月Y+X2轉變為B+X2+Y,其中,B在S3占據絕對優勢;X2在S5占據絕對優勢.5月浮游藻類主要功能群由4月B+X2+Y轉變為X2+D,其中,X2在S10占據絕對優勢;D在S11占據絕對優勢.整體來看,X2在陽澄西湖自南向北呈現逐漸上升再下降的趨勢;Y在陽澄西湖自南向北呈現逐漸下降再上升的趨勢.6月浮游藻類主要功能群由5月X2+D轉變為B+X2+J,其中,J在S5、S9均占據絕對優勢.7月浮游藻類主要功能群由6月B+X2+J轉變為B+X2+G+S1+Lo,其中,G在S1占據絕對優勢,Lo在S12占據絕對優勢.8月浮游藻類主要功能群由7月B+X2+G+S1+Lo轉變為B+X2+E+Y+ Lo,其中,X2在S1占據絕對優勢.整體來看,陽澄西湖浮游藻類功能群組成均衡度較高,J在陽澄西湖自南向北呈現逐漸上升再下降的趨勢.9月浮游藻類主要功能群由8月B+X2+E+Y+Lo轉變為X2+Y+Lo,主要優勢功能群在各采樣點均未占據絕對優勢.10月浮游藻類主要功能群由9月X2+Y+Lo轉變為X2+B+E,主要優勢功能群在各采樣點均未占據絕對優勢.11月浮游藻類主要功能群由10月X2+B+E轉變為X2+B+Y,主要優勢功能群在各采樣點均未占據絕對優勢.整體來看,X2在各采樣點均有分布,但未能占據絕對優勢種,各功能群在采樣點所占據的比例相對均衡.12月浮游藻類主要功能群與11月相同,均為X2+ B+Y,其中,X2在S1、S2占據絕對優勢.1月浮游藻類主要功能群由12月X2+B+Y轉變為X2+B+ Y+D,其中,D在S1、S2占據絕對優勢.2月浮游藻類主要功能群由1月X2+B+Y+D轉變為X2+Y,其中,Y在S11占據絕對優勢.整體來看,X2、Y在各采樣點均有分布,X2、D在陽澄西湖自南向北呈現逐漸下降的趨勢;Y在陽澄西湖自南向北呈現逐漸下降再上升的趨勢.

2.3 浮游藻類優勢功能群與環境因子的RDA分析

選取藻類優勢功能群與主要環境因子進行RDA分析,全年分析結果顯示(圖5),環境因子Chl- a、、pH值、DO、NH4+-N和Tur與浮游藻類顯著相關,環境因子第一軸解釋了29.9%的群落結構變化,第二軸解釋了37.8%的群落結構變化,前兩軸共解釋了群落結構的67.7%.B和Lo與Chl-a呈強正相關,與NH4+-N和Tur呈負相關;X2、D和E與DO和Chl-a呈正相關;X1、Y、G和J與pH、和DO呈正相關,與NH4+-N呈強負相關;S1與pH值和呈強正相關,與NH4+-N呈強負相關.

表2 環境因子與排序軸1、排序軸2的相關系數

圖5 陽澄西湖全年浮游藻類優勢功能群與主要環境因子的RDA分析

春季浮游藻類分析結果顯示(圖6),環境因子Chl-a、TP、NO2--N和Tur與浮游藻類顯著相關.其中,B、X1和Lo與Chl-a和NO2--N呈正相關;E、X2、D、Y和J與Chl-a、TP、NO2--N和Tur呈正相關.夏季環境因子pH值、DTN、DO、和Tur與浮游藻類顯著相關,其中,B和D與pH值、DO和Tur呈正相關,與呈負相關;Y和G與pH值、DTN和DO呈正相關;X1和J與、pH值、DO和DTN呈正相關;S1與pH值、DTN、DO和呈正相關,與Tur呈負相關;Lo與和DTN呈正相關,與Tur呈負相關;E和X2與Tur呈正相關,與DTN和呈負相關.秋季環境因子DTN、NH4+-N、NO2--N、DO和與浮游藻類顯著相關,其中,B與DO和NH4+-N呈正相關,與DTN呈負相關;X2與DTN、DO和NH4+-N呈正相關,與NO2--N和呈負相關;Y、E、J、Lo、S1和D與NO2--N和呈正相關,與DTN呈負相關;X1與呈強正相關,與NO2--N呈正相關,與DTN、NH4+-N和DO呈負相關.冬季環境因子DO和與浮游藻類顯著相關,其中,B、G和X2與DO呈強正相關,與呈正相關;D、E和Y與和DO呈正相關;X1和J與和DO呈負相關.

表3 浮游藻類功能群與環境因子的冗余分析(RDA)統計特征

圖6 陽澄西湖不同季節浮游藻類優勢功能群與主要環境因子的RDA分析

3 討論

3.1 浮游藻類功能群季節演替特征

浮游藻類功能群的演替受環境因子變化的影響,環境因子的變化會打破水體原來的生態平衡,浮游藻類可以通過群落演替的方式來達到新的生態平衡[36].一般條件下,浮游藻類對季節變化敏感性較強[37].陽澄西湖調查結果表明,不同時期浮游藻類功能群存在顯著差異.其中,功能群B、X2和Y在12個月調查中均有發現,是本次調查的絕對優勢功能群.研究表明,不同營養水平的熱帶水庫中,低營養鹽的水庫中主要出現的功能群組有B、X2、Y、E、M、Lo和 MP[17,38].B功能群代表性種類為小環藻,適應生境為中營養、淺水水體,并且耐受低溫和弱光[39],陽澄西湖不同時期的溫躍層以上水體溫度相對穩定,適宜該類藻類的生長.X2代表種衣藻和隱藻,在水體相對穩定的情況下,通過鞭毛運動到真光層捕獲光照[40];此外,以綠藻門為主的功能類群X2屬于投機性藻類,一旦營養鹽濃度和溫度適宜,便會迅速生長[41].Y功能群代表種為卵形隱藻,在低捕食壓力下能適宜生活在幾乎所有的靜水生態系統中[42],此外,Y功能群能快速吸收營養并迅速成長,并且可通過鞭毛運動在光照充足的表層與營養鹽充足的底部進行垂向遷移[43].

3月浮游藻類優勢功能群為B+Y+X2,4月浮游藻類優勢功能群為B+Y+X2+X1+D,比3月增加了X1和D,X1適應生境為混合程度較高的中-富營養淺水水體,對營養缺乏和濾食作用敏感[39].浮游藻類功能群的變化與光照、水體分層以及牧食壓力有關[44].X1占據優勢位置可能與魚類和浮游動物活動不強,X1牧食壓力小,且光照相對充足有關.D適合在低光照渾濁型水體中生存[9], 4月降水少,風擾動強烈,隨著T的升高,浮游藻類開始大量繁殖,導致水中的SD降低,為D類群提供了良好的生長條件[41].5月浮游藻類優勢功能群為B+Y+X2+J+D,隨著T的上升,魚類與浮游動物的活動增強,其他浮游藻類迅速繁殖,X1的優勢地位在營養缺失和牧食壓力下被J代替.J代表種多為非運動的不含膠質的綠藻,生活在富營養渾濁水體中,喜高溫[45],在夏季易形成優勢功能群.6月浮游藻類優勢功能群為B+Y+X2+J+G,G適宜在溫度14℃以上的穩定水體中生長[46],同時,較高的氮營養鹽也會為G提供適宜的條件,并占據優勢地位[9].7月浮游藻類優勢功能群為X2+B+G+S1+ Y+Lo,J中的綠藻對高光照敏感,隨著光照的逐漸增強,J的優勢地位被S1和Lo代替.S1多為絲狀藍藻,主要生活在中-富營養化水體,在高溫、強光和營養鹽充足的情況下,能夠大量繁殖.大量研究表明,S1在湖北省熊河水庫[47],廣東省的橫崗水濂山水庫以及南屏水庫[40]和江西省的贛江流域[49]夏季均占據優勢地位.Lo能夠在20~27℃的高溫穩定水體中生存,是夏季優勢功能群之一.Lo能夠通過鞭毛主動向光照充足的表層和營養鹽充足的底部運動,大量繁殖,使其處于競爭優勢地位[41].8月浮游藻類優勢功能群為X2+B+G+E+Y+Lo,E生活在貧-中營養型淺水水體,隨著7月陽澄西湖中大量的營養鹽被浮游藻類利用,8月水體中營養鹽的減少,E處于一定的優勢地位.浮游藻類的大量繁殖造成了水體中DO和SD降低,導致S1的生存條件遭到破壞,S1的優勢地位被E代替.9月浮游藻類優勢功能群為X2+B+Y+Lo,陽澄西湖水體中營養鹽處于穩定水平,E不能夠占據優勢地位,Lo通過鞭毛主動運動占據優勢地位.10月浮游藻類優勢功能群為X2+B+Y+E,水體中的溫度和光照均不適宜Lo大量繁殖,Lo不再占據優勢地位.水體中營養鹽的降低使得生活貧-中營養型水體中的E再次占據優勢地位.11月浮游藻類優勢功能群為絕對優勢功能群X2+B+Y.12月浮游藻類優勢功能群為X2+B+Y+X1,隨著T降低,魚類和浮游動物活動的降低,X1牧食壓力降低,X1重新占據優勢地位.1月浮游藻類優勢功能群與4月相同,優勢功能群為B+Y+X2+X1+D.研究表明,D在冬季占優勢,主要是由于冬季T較低,適應較低T的硅藻生長[49].2月浮游藻類優勢功能群與3月和11月相同,只有B+Y+X2占據絕對優勢地位.

3.2 浮游藻類功能群與環境因子的相關性

研究表明,pH值、、DO、SD、TN、TP和NH4+-N等環境因子均能影響浮游藻類群落的分布[50-54].本研究全年分析結果顯示,浮游藻類群落分布的主要影響環境因子有Chl-a、、pH值、DO、NH4+-N和Tur.Chl-a是反映浮游植物的重要指標之一,與浮游植物的種類、密度以及生物量密切相關.一般而言,浮游植物密度越高,Chl-a濃度越高,因此本研究中Chl-a不能作為影響浮游植物分布的主要環境因子[55].研究表明,DO與藻類生長呈顯著正相關[56], DO是B、D、J、MP、P、Y分布的主要限制因素.對寒區安邦河濕地[46]、北方的烏梁素海[57]、深水型金盆水庫[58]、貴州高原的阿哈水庫[59]與三板溪水庫[46]等湖庫的研究表明,是該區域藻類功能群結構變化的最主要影響因素.陽澄西湖浮游藻類絕對優勢功能群為B、X2和Y.多元回歸分析結果顯示,B與NH4+-N呈極顯著負相關(<0.01),相關系數為-0.384;與Chl-a和呈極顯著正相關(<0.01),相關系數分別為0.453和0.245.X2與Chl-a和呈極顯著正相關(<0.01),相關系數分別為0.189和0.207;與SD呈顯著正相關(<0.05),相關系數為0.414.Y與NH4+-N呈顯著負相關(<0.05),相關系數為-0.248;與Chl-a、和CODMn呈顯著正相關(<0.05),相關系數分別為0.114,0.121,0.320.冗余分析(RDA)與多元回歸分析結果基本一致,其中,NH4+-N和是主要影響陽澄西湖全年浮游藻類分布的環境因子.此外,多元回歸分析發現,X2與SD呈顯著正相關(<0.05),這一結果與林巋璇等[60]研究結果一致,SD為X2、S1、J、X1和W2等功能群的主要限制因子.

春季陽澄西湖出現功能群X1、J和D的演替,其中,J和D與Tur呈正相關,5月為春季農業灌溉時期,污染物隨水流進入湖區,造成陽澄西湖水體Tur上升,透明度較低,水下的光照較弱,因此,有利于J和D功能群在春季及夏初形成優勢功能群[59].研究表明,J對沖刷有一定的耐受性,總體來說J和D均適合生長在Tur較高的富營養水體中[55].夏季陽澄西湖出現功能群J、G、S1、Lo和E的演替,pH值、DTN、DO、T和Tur成為主要影響演替的主要環境因子.S1和J均與pH值、DTN、DO和呈正相關.陽澄西湖受人類活動影響較大,夏季升高,豐水期降雨量豐富,地表徑流帶來大量外源物質的輸入,導致水體渾濁,透明度降低,水下光強減弱,有利于S1和J功能群的生長繁殖[47].Romo等[61]發現,S1與J常共生于渾濁的富營養化水體中,這在本研究中得到了驗證.研究表明[62],藻類適宜在偏堿性的水體中生長,pH值成為適宜富營養水體生長的J功能群的限制因子,這一結果與本研究結果一致.秋季DTN、NH4+-N、NO2--N、DO和成為主要影響演替的主要環境因子.其中,J和Lo與呈正相關,是影響浮游植物生長最直接的環境因子[63],降低導致J、G和Lo功能群失去優勢位置.B和X2與NH4+-N和呈正相關;Y和E與和NO2--N呈正相關,結果表明陽澄西湖秋季浮游藻類受到氮的限制.這一結果與汪星等人的研究結果一致[62].冬季陽澄西湖浮游藻類的主要限制因子為DO和,其中,B、Y和X2均與和DO呈正相關,受到和DO的限制.冬季陽澄西湖水體中氮含量相對較高,X1類群作為“先鋒種”可在資源競爭中占據優勢[9,39],較高濃度的氮營養鹽也可以維持纖維藻的藻細胞較長時間的增長[64].

4 結論

4.1 2015~2016年陽澄西湖浮游藻類可分為21個功能群:A、B、D、E、F、G、H1、J、Lo、M、MP、P、S1、T、TD、W1、W2、X1、X2、X3和Y.

4.2 陽澄西湖不同時期浮游藻類功能群存在差異,其中,B、X2和Y在12次調查中優勢度均>0.02,為陽澄西湖絕對優勢功能群;陽澄西湖不同季節優勢功能群演替規律為:由春季的B+X2+Y+X1+J+D,經夏季的B+X2+Y+X1+J+G+S1+Lo+E,秋季B+X2+ Y+Lo+E,轉變為冬季的B+X2+Y+X1+D.

4.3 RDA分析結果顯示,陽澄西湖浮游藻類功能群分布受環境影響較為明顯,但不同季節影響因子不同.整體上,、pH值、DO、NH4+-N和Tur是影響陽澄西湖浮游藻類功能群分布格局的主要因素.

[1] Sommer U, Adrian R, De Senerpont Domis L,et al. Beyond the Plankton Ecology Group (PEG) model: mechanisms driving plankton succession [J]. Annual Review of Ecology, Evolution, and Systematics, 2012,43:429-448.

[2] Tippett Roger. The ecology of freshwater phytoplankton (C. S. Reynolds) [J]. Limnology and Oceanography, 1987,32(3):779-780.

[3] Padisak J, Borics G, Grigorszky I, et al. Use of phytoplankton assemblages for monitoring ecological status of lakes within the Water Framework Directive: the assemblage index [J]. Hydrobiologia, 2006,553(1):1-14.

[4] 沈韞芬,章宗涉,龔循矩.微型生物監測新技術[M]. 北京:中國建筑工業出版社, 1990:119. Shen Y F, Zhang Z S, Gong X Z. New technology of microbiological monitoring [M]. Beijing: China Building Industry Press, 1990:119.

[5] Bonilla S, Conde D, Aubriot L, et al. Influence of hydrology on phytoplankton species composition and life strategies in a subtropical coastal lagoon periodically connected with the Atlantic Ocean [J]. Estuaries, 2005,28(6):884-895.

[6] Bonilla S, Conde D, Aubriot L, et al. Influence of hydrology on phytoplankton species composition and life strategies in a subtropical coastal lagoon periodically connected with the Atlantic Ocean [J]. Estuaries, 2005,28(6):884-895.

[7] Yu H X, Wu J H, Ma C X, et al. Seasonal dynamics of phytoplankton functional groups and its relationship with the environment in river: a case study in northeast China [J]. Journal of Freshwater Ecology, 2012,27(3):429-441.

[8] 董 靜,李艷暉,李根保,等.東江水系浮游植物功能群季節動態特征及影響因子[J]. 水生生物學報, 2013,37(5):836-843. Dong J, Li Y H, Li G B, et al. Seasonal dynamics and influencing factors of phytoplankton functional groups in the Dongjiang River system [J]. Journal of Hydrobiology, 2013,37(5):836-843.

[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] Brantes N, Antunes S C, Pereira M J,et al. Seasonal succession of cladocerans and phytoplankton and their interactions in a shallow eutrophic lake (Lake Vela, Portugal) [J]. Acta Oecologica, 2006, 29(1):54-64.

[11] 陳曉江,楊 劼,杜桂森,等.官廳水庫浮游植物功能群季節演替及其驅動因子 [J]. 中國環境監測, 2016,32(3):74-81. Chen X J, Yang J, Du G S, et al. Seasonal succession of phytoplankton functional groups in Guanting Reservoir and its driving factors [J]. China Environmental Monitoring, 2016,32(3):74-81.

[12] Becker V, Caputo L, Ordó?ez J, et al. Driving factors of the phytoplankton functional groups in a deep Mediterranean reservoir [J]. Water Research, 2010,44(11):3345-3354.

[13] Pintotd S, Becker V. Diel dynamic of phytoplankton functional groups in a tropical water supply, Extremoz Lake, northeastern Brazil [J]. Acta Limnol Bras, 2014,26(4):356-366.

[14] Costa L S, Huszar V L M, Ovalle A R. Phytoplankton functional groups in a tropical estuary: hydrological control and nutrient limitation [J]. Estuaries and Coasts, 2009,32(3):508-521.

[15] Stankovi? I, Vlahovi? T, UDOvi? M G, et al. Phytoplankton functional and morpho-functional approach in large floodplain rivers [J]. Hydrobiologia, 2012,698(1):217-231.

[16] 田澤斌,劉德富,姚緒姣,等.分層對香溪河庫灣浮游植物功能群季節演替的影響[J]. 長江流域資源與環境, 2014,23(5):700-707. Tian Z B, Liu D F, Yao X Z, et al. Effects of stratification on seasonal succession of phytoplankton functional groups in Xiangxihe Reservoir Bay [J]. Resources and Environment in the Yangtze Basin, 2014,23(5):700-707.

[17] 高國敬,肖利娟,林秋奇,等.海南省典型水庫浮游植物功能類群的結構特征與水質評價[J]. 生態科學, 2013,32(2):144-150. Gao G J, Xiao L J, Lin Q Q, et al. Structural characteristics and water quality evaluation of phytoplankton functional groups in typical reservoirs of Hainan province [J]. 2013,32(2):144-150.

[18] 尹子龍,程李芳,于洪賢,等.扎龍濕地夏季浮游植物功能群特征及水質評價[J]. 北方農業學報, 2016,44(4):56-62. Yin Z L, Cheng L F, Yu H X, et al. Characteristics of phytoplankton functional groups and water quality evaluation in Zhalong wetland in summer [J]. Northern Agricultural Journal, 2016,44(4):56-62.

[19] 董 靜,李根保,宋立榮.撫仙湖、洱海、滇池浮游藻類功能群1960s以來演變特征[J]. 湖泊科學, 2014,26(5):735-742. Dong J, Li G B, Song L R. Evolution characteristics of phytoplankton functional groups in Fuxian Lake, Bohai Sea and Dianchi Lake since 1960s [J]. Journal of Lake Science, 2014,26(5):735-742.

[20] 翁建中,李繼影,徐恒省,等.陽澄湖浮游植物研究及其富營養化評價[J]. 環境科學與管理, 2009,34(10):41-44,49. Wen J Z, Li J Y, Xu H S, et al. Phytoplankton research and eutrophication evaluation in Yangcheng Lake [J]. Environmental Science and Management, 2009,34(10):41-44,49.

[21] 駱東玲.淺水湖泊富營養化的機理與應對策略-以陽澄湖為例[J]. 農業資源與環境學報, 2007,24(3):14-18.Luo D L. Mechanism and coping strategies of eutrophication in shallow lakes: A case study of Yangcheng Lake [J]. Journal of Agricultural Resources and Environment, 2007,24(3):14-18.

[22] 葛 優,周彥鋒,呂大偉,等.陽澄西湖葉綠素a的時空分布及其與環境因子的關系[J]. 長江流域資源與環境, 2017,26(7):1068-1075.Ge Y, Zhou Y F, Lv D W, et al. Temporal and spatial distribution of chlorophyll-a in Yangcheng West Lake and its relationship with environmental factors [J]. Resources and Environment in the Yangtze Basin, 2017,26(7):1068.

[23] 丁 娜,周彥鋒,宋江騰,等.陽澄湖浮游植物群落結構特征及其與水環境因子的關系[J]. 上海海洋大學學報, 2015,24(5):702-711. Ding N, Zhou Y F, Song J T, et al. Characteristics of phytoplankton community structure and its relationship with water environment factors in Yangcheng Lake [J]. Journal of Shanghai Ocean University, 2015,24(5):702-711.

[24] 景 明,李繼影,劉孟宇,等.陽澄湖浮游植物群落結構與環境因子的典范對應研究[J]. 環境科學與管理, 2014,39(11):145-148. Jing M, Li J Y, Liu M Y, et al. Canonical correspondence between phytoplankton community structure and environmental factors in Yangcheng Lake [J]. Environmental Science and Management, 2014,39(11):145-148.

[25] 蔡曉鈺,姜 宇,張占恩,等.陽澄湖浮游植物調查及其與水體營養化關系分析[J]. 治淮, 2015,(6):15-17.Cai X Y, Jiang Y, Zhang Z E, et al. Investigation of phytoplankton in Yangcheng Lake and its relationship with nutrientization of water body [J]. Journal of Management, 2015,(6):15-17.

[26] 呂 文,孫瑞瑞,王 誠,等.陽澄湖西湖水源地富營養化時空變化特征分析[J]. 水資源與水工程學報, 2018,(3):87-90.Lv W, Sun R R, Wang C, et al. Temporal and spatial variation characteristics of eutrophication in the water source of West Lake in Yangcheng Lake [J]. Journal of Water Resources and Water Engineering, 2018,(3):87-90.

[27] 陳偉民,黃祥飛,周萬平.湖泊生態系統觀測方法[M]. 北京:中國環境科學出版社, 2005. Chen W M, Huang X F, Zhou W P. Lake ecosystem observation methods [M]. Beijing: China Environmental Science Press, 2005.

[28] 金相燦,屠清英.湖泊富營養化調查規范[M]. 北京:中國環境科學出版社, 1990:198-210. Jin X C, Tu Q Y. Regulations for Lake eutrophication survey [M]. Beijing: China Environmental Science Press, 1990:198-210.

[29] 國家環保局本書編委會.水和廢水監測分析方法 [M]. 中國環境科學出版社, 2002:243-280. Editorial Board of the National Environmental Protection Agency. Water and wastewater monitoring and analysis methods [M]. China Environmental Science Press, 2002:243-280.

[30] 章宗涉,黃祥飛.淡水浮游生物研究方法[M]. 北京:科學出版社, 1995. Zhang Z S, Huang X F. Research methods of freshwater plankton [M]. Beijing: Science Press, 1995.

[31] 周鳳霞.淡水微型生物與底棲動物圖譜[M]. 北京:化學工業出版社, 2011.Zhou F X. Freshwater micro-bio and benthic maps [M]. Beijing: Chemical Industry Press, 2011.

[32] 胡鴻鈞,魏印心.中國淡水藻類-系統、分類及生態[J]. 2006.Hu H J, Wei Y X. Freshwater algae in China-system, classification and ecology [J]. 2006.

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

[34] 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.

[35] 王 華,楊樹平,房晟忠,等.滇池浮游植物群落特征及與環境因子的典范對應分析[J]. 中國環境科學, 2016,36(2):544-552.Wang H, Yang S P, Fang S Z, et al. Canonical correlation analysis of phytoplankton community characteristics and environmental factors in Dianchi Lake [J]. China Environmental Science, 2016,36(2):544-552.

[36] Reynolds C S. What factors influence the species composition of phytoplankton in lakes of different trophic status? [J]. Hydrobiologia, 1998,369(1):11-26.

[37] 杜 虹,黃長江,陳善文,等.2001~2002年粵東柘林灣浮游植物的生態學研究 [J]. 海洋與湖沼, 2003,34(6):604-617. Du H, Huang C J, Chen S W, et al. Ecological study of phytoplankton in Zhelin Bay, Eastern Guangdong Province, 2001~2002 [J]. Ocean and Lake Journal, 2003,34(6):604-617.

[38] 黃享輝,胡 韌,雷臘梅,等.南亞熱帶典型中小型水庫浮游植物功能類群季節演替特征[J]. 生態環境學報, 2013,22(2):311-318. Huang X H, Hu R, Lei L M, et al. Seasonal succession characteristics of phytoplankton functional groups in typical small and medium-sized reservoirs in South Asia [J]. Journal of Eco-Environment, 2013,22(2): 311-318.

[39] 胡 韌,藍于倩,肖利娟,等.淡水浮游植物功能群的概念、劃分方法和應用[J]. 湖泊科學, 2015,27(1):11-23. Hu R, Lan Y Q, Xiao L J, et al. Concept, division method and application of freshwater phytoplankton functional groups [J]. Journal of Lake Science, 2015,27(1):11-23.

[40] 張 怡,胡 韌,肖利娟,等.南亞熱帶兩座不同水文動態的水庫浮游植物的功能類群演替比較[J]. 生態環境學報, 2012,21(1):107-117. Zhang Y, Hu R, Xiao L J, et al. Comparison of functional group succession of two phytoplankton reservoirs in different hydrological dynamics in South Asia [J]. Chinese Journal of Eco-Environment, 2012,21(1):107-117.

[41] 鄭 誠,陸開宏,徐 鎮,等.四明湖水庫浮游植物功能類群的季節演替及其影響因子[J]. 環境科學, 2018,39(6):2688-2697. Zheng C, Lu K H, Xu Z, et al. Seasonal succession of phytoplankton functional groups in the Siming Lake Reservoir and its influencing factors [J]. Environmental Science, 2018,39(6):2688-2697.

[42] Padisák 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.

[43] Bovo-Scomparin V M, Train S. Long-term variability of the phytoplankton community in an isolated floodplain lake of the Ivinhema River State Park, Brazil [J]. Hydrobiologia, 2008,610(1): 331-344.

[44] Becker V, Vera Lúcia M. Huszar, Crossetti L O . Responses of phytoplankton functional groups to the mixing regime in a deep subtropical reservoir [J]. Hydrobiologia, 2009,628(1):137-151.

[45] 武安泉,郭 寧,覃雪波.寒區典型濕地浮游植物功能群季節變化及其與環境因子關系 [J]. 環境科學學報, 2015,35(5):1341-1349. Wu A Q, Guo N, Qi X B. Seasonal changes of phytoplankton functional groups in typical wetlands in cold regions and their relationship with environmental factors [J]. Journal of Environmental Science, 2015,35(5):1341-1349.

[46] 黃國佳,李秋華,陳 椽,等.貴州高原紅楓湖水庫浮游植物功能分組及其時空分布特征 [J]. 生態學報, 2015,35(17):5573-5584. Huang G J, Li Q H, Chen W, et al. Functional grouping of phytoplankton in Hongfeng Lake Reservoir of Guizhou Plateau and its spatial and temporal distribution characteristics [J]. Chinese Journal of Ecology, 2015,35(17):5573-5584.

[47] 張 婷,李 林,宋立榮.熊河水庫浮游植物群落結構的周年變化 [J]. 生態學報, 2009,29(6):2971-2979. Zhang T, Li L, Song L R. Annual variation of phytoplankton community structure in Xionghe Reservoir [J]. Journal of Ecology, 2009,29(6):2971-2979.

[48] 劉足根,張 柱,張 萌,等.贛江流域浮游植物群落結構與功能類群劃分[J]. 長江流域資源與環境, 2012,21(3):375. Liu Z G, Zhang Z, Zhang M, et al. Division of phytoplankton community structure and functional groups in Ganjiang River Basin [J]. Resources and environment in the Yangtze River Basin, 2012,21 (3):375.

[49] 代龔圓,李 杰,李 林,等.滇池北部湖區浮游植物時空格局及相關環境因子 [J]. 水生生物學報, 2012,36(5):946-956. Dai G Y, Li J, Li L, et al. Spatial-temporal pattern of phytoplankton and related environmental factors in the Northern Lake area of Dianchi Lake [J]. Journal of Aquatic Biology, 2012,36(5):946-956.

[50] 汪 星,劉錄三,李 黎,等.鏡泊湖浮游藻類組成及其與環境因子的相關分析 [J]. 中國環境科學, 2015,35(11):3403-3413. Wang X, Liu L S, Li L, et al. Composition of phytoplankton in Jingbo Lake and its correlation with environmental factors [J]. China Environmental Science, 2015,35(11):3403-3413.

[51] 沈會濤,劉存歧.白洋淀浮游植物群落及其與環境因子的典范對應分析 [J]. 湖泊科學, 2008,20(1):773-779. Shen H T, Liu C Q. Phytoplankton community in Baiyangdian Lake and its canonical correspondence with environmental factors [J]. Lake Science, 2008,20(1):773-779.

[52] 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.

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

[54] 潘繼征,熊 飛,李文朝,等.云南撫仙湖透明度的時空變化及影響因子分析[J]. 湖泊科學, 2008,20(5):681-686. Pan J Z, Xiong F, Li W C, et al. Temporal and spatial changes and influence factors of transparency in Fuxian Lake, Yunnan [J]. Journal of Lake Science, 2008,20(5):681-686.

[55] 李 磊,李秋華,焦樹林,等.小關水庫夏季浮游植物功能群對富營養化特征的響應[J]. 環境科學, 2015,36(12):4436-4443. Li L, Li Q H, Jiao S L, et al. Response of phytoplankton functional groups to eutrophication characteristics in Xiaoguan reservoir in Summer [J]. Environmental Science, 2015,36(12):4436-4443.

[56] 游 亮,崔莉鳳,劉載文,等.藻類生長過程中DO、pH與葉綠素相關性分析[J]. 環境科學與技術, 2007,30(9):42-45. You L, Cui L F, Liu Z W, et al. Correlation analysis of DO, pH and chlorophyll in algae growth process [J]. Environmental Science and Technology, 2007,30(9):42-45.

[57] 李 興,李建茹,徐效清,等.烏梁素海浮游植物功能群季節演替規律及影響因子[J]. 生態環境學報, 2015,24(10):1668-1675. Li X, Li J R, Xu X Q, et al. Seasonal succession and influencing factors of phytoplankton functional groups in Wuliangsuhai lake [J]. Chinese Journal of Eco-Environment, 2015,24(10):1668-1675.

[58] 盧金鎖,張 穎,胡亞潘.深水型水庫中藻類功能群組演替及其與環境因子的關系[J]. 環境工程學報, 2014,(11):4605-4611. Lu J S, Zhang Y, Hu Y P. The succession of algal functional groups in deep water reservoirs and its relationship with environmental factors [J]. Journal of Environmental Engineering, 2014,(11):4605-4611.

[59] 李 磊,李秋華,焦樹林,等.阿哈水庫浮游植物功能群時空分布特征及其影響因子分析 [J]. 環境科學學報, 2015,35(11):3604-3611. Li L, Li Q H, Jiao S L, et al. Temporal and spatial distribution characteristics of phytoplankton functional groups in Aha Reservoir and their influencing factors [J]. Journal of Environmental Science, 2015,35(11):3604-3611.

[60] 林巋璇,汪 星,王 瑜,等.鏡泊湖浮游藻類功能群的演替特征及其影響因素[J]. 環境科學研究, 2017,30(4):520-528. Lin W X, Wang X, Wang Y et al. Succession characteristics and influencing factors of phytoplankton functional groups in Jingbo Lake [J]. Environmental Science Research, 2017,30(4):520-528.

[61] Romo S. Phytoplankton strategies and diversity under different nutrient levels and planktivorous fish densities in a shallow Mediterranean lake [J]. Journal of Plankton Research, 2005,27(12): 1273-1286.

[62] 汪 星,李利強,鄭丙輝,等.洞庭湖浮游藻類功能群的組成特征及其影響因素研究[J]. 中國環境科學, 2016,36(12):3766-3776. Wang X, Li L Q, Zheng B H, et al. Composition characteristics and influencing factors of phytoplankton functional groups in Dongting Lake [J]. China Environmental Science, 2016,36(12):3766-3776.

[63] 田永強.淡水浮游植物功能類群劃分方法及其生態學應用研究進展[J]. 亞熱帶植物科學, 2015,44(4):349-354. Tian Y Q. Progress in the classification of functional groups of freshwater phytoplankton and its application in ecology [J]. Subtropical Plant Science, 2015,44(4):349-354.

[64] 王英英.不同氮磷質量濃度對太浦河四種優勢藻類生長影響的研究[D]. 上海:上海師范大學, 2016. Wang Y Y. Study on the effects of different nitrogen and phosphorus concentrations on the growth of four dominant algae in Taipu River [D]. Shanghai: Shanghai Normal University, 2016.

Succession patterns of phytoplankton functional groups in western area of Yangcheng Lake and their relationship with environmental factors.

GE You1, ZHOU Yan-feng1, WANG Cheng-he2, YOU Yang1,2*

(1.Key Laboratory of Freshwater Fisheries and Utilization of Germplasm Resources, Ministry of Agriculture, Freshwater Fisheries Research Center, Chinese Academy of Fishery Sciences,Wuxi 214081, China;2.Wuxi Fishery College, Nanjing Agricultural University, Wuxi 214081, China)., 2019,39(7):3027~3039

The functional classification of phytoplankton assemblages in western area of Yangcheng Lake and the influencing factors were explored based upon their community composition structure analysis and by using the redundancy analysis and multiple regression methods. Samples were collected monthly from March 2015 to February 2016 from 12 sites distributed in western area of Yangcheng Lake, and the phytoplankton and environmental variables were investigated. A total of 21 phytoplankton functional groups were identified as A、B、D、E、F、G、H1、J、Lo、M、MP、P、S1、T、TD、W1、W2、X1、X2、X3 and Y. And they were dominated by the 10groups of B、D、E、G、J、Lo、S1、X1、X2 and Y. The western water of Yangcheng Lake was a shallow medium-eutrophic water body with less disturbance and high mixing degree. Meanwhile, the dominant functional groups of phytoplankton in western area of Yangcheng Lake followed the seasonal succession pattern, that was B+X2+Y+X1+J+D in spring, B+X2+Y+X1+J+G+S1+Lo+E in summer, B+X2+Y+Lo+E in autumn and B+X2+Y+X1+D in winter. As well, the distribution of these functional groups in west area of Yangcheng Lake was closely related to the environmental variables. Water temperature, pH, DO, ammonia, and turbidity were the main influencing factors.

western area of Yangcheng Lake;phytoplankton functional groups;dominance;environmental factors

X826

A

1000-6923(2019)07-3027-13

葛 優(1991-),男,江蘇高郵人,碩士,主要從事生態漁業方面的研究.發表論文1篇.

2018-12-13

國家重點研發計劃資助(2018YFD0900905);國家現代農業產業技術體系(CARS-46)

* 責任作者, 研究員, youy@ffrc.cn

猜你喜歡
優勢環境功能
也談詩的“功能”
中華詩詞(2022年6期)2022-12-31 06:41:24
長期鍛煉創造體內抑癌環境
一種用于自主學習的虛擬仿真環境
矮的優勢
趣味(語文)(2020年3期)2020-07-27 01:42:46
孕期遠離容易致畸的環境
環境
關于非首都功能疏解的幾點思考
畫與話
談“五老”的五大特殊優勢
中國火炬(2014年11期)2014-07-25 10:31:58
中西醫結合治療甲狀腺功能亢進癥31例
主站蜘蛛池模板: 97亚洲色综久久精品| 99久久精品视香蕉蕉| 久久久久国产一区二区| 99热国产在线精品99| 国产精品女人呻吟在线观看| 久草视频精品| 啪啪啪亚洲无码| 自拍中文字幕| 亚洲天堂高清| 国产二级毛片| 亚洲无线视频| 91无码视频在线观看| 亚洲浓毛av| 香蕉视频国产精品人| 欧美α片免费观看| 91视频免费观看网站| 国产在线一二三区| 亚洲AV一二三区无码AV蜜桃| 欧美日本在线播放| 国内精自视频品线一二区| 国产无码网站在线观看| 中文字幕波多野不卡一区| 亚洲第一中文字幕| 激情视频综合网| 丰满人妻被猛烈进入无码| 99久久精品免费看国产电影| 91小视频版在线观看www| 伊人天堂网| 激情综合网址| 久久国产黑丝袜视频| 国产精欧美一区二区三区| 91免费国产高清观看| 久久成人国产精品免费软件| 婷婷色婷婷| 在线va视频| 精品久久人人爽人人玩人人妻| 蝌蚪国产精品视频第一页| 2021国产v亚洲v天堂无码| 亚洲午夜天堂| 中日韩欧亚无码视频| 日韩高清中文字幕| 日韩色图在线观看| 伊人国产无码高清视频| 国产伦片中文免费观看| 亚洲日韩精品欧美中文字幕| 91香蕉视频下载网站| 亚洲an第二区国产精品| 亚洲天堂精品在线观看| 亚洲精选无码久久久| 波多野结衣国产精品| 思思热在线视频精品| 亚洲视频黄| 久久毛片网| 免费视频在线2021入口| 2021国产精品自拍| 天天摸夜夜操| 欧美影院久久| 天堂av综合网| 久久国产精品电影| 婷婷综合亚洲| 免费在线观看av| 毛片网站在线看| 国产一区二区三区精品久久呦| 在线观看国产网址你懂的| 亚洲香蕉伊综合在人在线| 毛片在线播放a| 免费女人18毛片a级毛片视频| 中文字幕无码中文字幕有码在线| 久久美女精品国产精品亚洲| 人妻丰满熟妇αv无码| 亚洲有无码中文网| 波多野结衣在线se| 三级欧美在线| 国产精品女同一区三区五区| 日韩在线成年视频人网站观看| 九九九九热精品视频| 国产成人无码Av在线播放无广告| 精品小视频在线观看| 欧美特黄一级大黄录像| 综合色区亚洲熟妇在线| 久久久精品无码一二三区| 国产欧美日韩视频怡春院|