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

    向家壩水庫浮游植物群落時空特征及其環(huán)境驅動因子

    2025-03-22 00:00:00田楚銘張地繼米武娟齊青松邢龍孫志峰曹光榮畢永紅
    水生態(tài)學雜志 2025年2期

    摘要:分析向家壩水庫浮游植物群落組成及其與環(huán)境因子的耦合關系,為水庫生態(tài)系統(tǒng)的保護與管理提供基礎數(shù)據(jù)和參考依據(jù)。2019年1月—2020年1月在向家壩水庫干支流20個采樣點逐月采集浮游植物樣品,同時測定水溫、溶解氧、電導率等環(huán)境因子,采用α多樣性指數(shù)分析浮游植物群落特征,并通過RDA分析其與環(huán)境因子的相關性。結果共采集并鑒定出浮游植物119種,隸屬于8門79屬;群落優(yōu)勢種6種,主要是小環(huán)藻(Cyclotella sp.)、變異直鏈藻(Melosira varians)等,群落為綠藻-硅藻-藍藻型,隱藻、甲藻、裸藻、金藻等種類較少,占比僅15.12%。浮游植物群落的時空差異較大,不同季節(jié)的細胞密度大小為夏季gt;春季gt;秋季gt;冬季,且夏季遠高于其他季節(jié);空間上細胞密度表現(xiàn)為支流大于干流。水環(huán)境因子的時空差異明顯,水溫、濁度和流速平均值表現(xiàn)為夏季gt;秋季gt;春季gt;冬季,總氮平均濃度為夏季gt;冬季gt;春季gt;秋季,總磷和硝態(tài)氮表現(xiàn)為春夏季gt;秋冬季,氨氮、溶解氧和pH表現(xiàn)為春冬季gt;夏秋季,高錳酸鹽指數(shù)平均濃度為春季gt;秋季gt;冬季gt;夏季;空間上高錳酸鹽指數(shù)、電導率、濁度、透明度和流速為干流gt;支流,其中流速具有極顯著差異(Plt;0.01),濁度具有顯著差異(Plt;0.05),總氮、硝態(tài)氮、氨氮、水溫、溶解氧和氧化還原電位顯示為支流gt;干流,其中溶解氧具有顯著差異(Plt;0.05)。單因素方差分析表明,向家壩水庫環(huán)境因子間具有顯著的時空異質性(Plt;0.05)。RDA分析表明,影響浮游植物群落時空變化的主要因子為水溫、營養(yǎng)鹽、透明度和電導率,其次為流速和濁度。

    關鍵詞:浮游植物;環(huán)境因子;冗余分析;向家壩水庫

    中圖分類號:Q145+.2" " " " 文獻標志碼:A" " " " 文章編號:1674-3075(2025)02-0194-10

    浮游植物是水域生態(tài)系統(tǒng)重要的初級生產(chǎn)者,食物鏈的基礎環(huán)節(jié)(Herrera-Silveira amp; Morales-Ojeda,2009),參與完成生態(tài)系統(tǒng)內(nèi)物質循環(huán)和能量流動過程,對構建群落環(huán)境、維持生態(tài)平衡、供應能源等諸多方面都具有重要意義(Reynolds,2006)。實質上,水域生態(tài)系統(tǒng)的結構功能與浮游植物群落的種類組成和數(shù)量特征息息相關,浮游植物群落演替能顯著影響水域生態(tài)系統(tǒng)的平衡(Ha et al,1998)。相關研究表明,浮游植物種類、數(shù)量特征及其時空格局的變化與營養(yǎng)鹽、溫度、光照、水動力條件等環(huán)境因子具有耦合關系(van Gremberghe,2008;Becker et al,2010;Vogt et al,2015)。不同水體浮游植物群落結構特征不同,因此,浮游植物因其對不同水環(huán)境的敏感響應和特異適應性而成為水環(huán)境監(jiān)測的重要生物指示物種(張春梅等,2021)。

    水庫中的浮游植物對水庫生態(tài)系統(tǒng)的結構功能具有重要意義,同時,浮游植物群落與環(huán)境因子間存在密切的耦合關系,開展水庫浮游植物群落調查研究對于水庫的運行管理具有重要的支撐作用。向家壩水庫位于四川省宜賓市(左岸)和云南省水富市(右岸)交界的金沙江下游河段,是金沙江下游4座特大型水庫(烏東德、白鶴灘、溪洛渡和向家壩)之一(張忠波,2014);總庫容為51.63×108 m3,調節(jié)庫容為9×108 m3,多年平均徑流量為3 810 m3/s;屬熱帶季風氣候,年平均氣溫18.3 ℃,年平均降水量為908.1 mm,豐水期多年平均降水量占全年降水的91.8%(王耀耀等,2019);屬于典型的山區(qū)河道型水庫(梁秀蓮,2021;米武娟等,2022)。向家壩水庫的運行具有重要的社會經(jīng)濟效益,但同時也使原有的河流生態(tài)系統(tǒng)和環(huán)境發(fā)生了深刻的改變(龍良紅等,2017;王耀耀等,2019;米武娟等,2022)。河流建壩成庫改變了原有的水文情勢進而改變了水域生態(tài)系統(tǒng)的結構功能狀況,向家壩水庫作為新建水庫,其水域生態(tài)系統(tǒng)結構功能備受關注;鑒于浮游植物是重要的環(huán)境指示生物,從浮游植物的角度深入認識向家壩水庫的生態(tài)環(huán)境狀況具有必要性和迫切性。目前,還沒有向家壩水庫浮游植物群落方面的研究報道。為此,本文以向家壩水庫為對象,于2019至2020年對浮游植物進行了逐月調查,分析了向家壩水庫浮游植物群落組成及環(huán)境因子的動態(tài)特征,探討了浮游植物與環(huán)境因子的耦合關系,擬為水庫生態(tài)系統(tǒng)的保護與管理提供基礎數(shù)據(jù)和參考依據(jù)。

    1" "材料與方法

    1.1" "研究區(qū)域與采樣點

    2019年1月至2020年1月在向家壩水庫的干流及庫灣支流等水域(32.68°~39.99° N,111.72°~116.27° E)進行逐月調查,每月采集樣品1次,共計采樣13次。干流采樣點8個(XJB01~08),支流采樣點12個(XB01~12),共設置20個采樣點(圖1)。

    1.2" "樣品的采集與處理

    浮游植物定量樣品:用5 L玻璃采樣器在每個采樣點的水面下0.5 m處取1 L的水樣,立刻加入10 mL 10%魯哥氏液固定,樣品帶回實驗室靜置48 h后,采用虹吸法濃縮至30 mL,再搖勻后吸取0.1 mL至浮游生物計數(shù)框在Olympus CX23光學顯微鏡下進行觀察、鑒定和計數(shù)(章宗涉和黃祥飛,1991)。鑒定工作參考的主要資料為《中國淡水藻類——系統(tǒng)、分類及生態(tài)》(胡鴻鈞和魏印心,2006)及《中國內(nèi)陸水域常見藻類圖譜》(鄧堅,2012)等。

    浮游植物定性樣品:使用25號浮游生物網(wǎng)在水面下0.5 m處呈“∞”形拖拽采集定性樣品,用適量4%甲醛溶液固定樣品(章宗涉和黃祥飛,1991)。

    現(xiàn)場使用YSI多參數(shù)水質分析儀現(xiàn)場測定水溫(WT)、溶解氧(DO)、電導率(SPC)、氧化還原電位(ORP)和pH,采用塞氏盤測定透明度(SD)。按照《水和廢水監(jiān)測分析方法》(國家環(huán)境保護總局, 2002)對各項理化指標進行測定,包括總磷(TP)、總氮(TN)、硝態(tài)氮(NO3--N)、氨氮(NH4+-N)、正磷酸鹽(PO43--P)、高錳酸鹽指數(shù)(CODMn)、葉綠素a(Chl-a)(國家環(huán)境保護總局,2002)。使用流速儀和濁度計測定流速(Vel)和濁度(Tur)。

    1.3" "數(shù)據(jù)處理與分析

    優(yōu)勢度計算公式為:

    Y=(Ni/N)×fi" " " " " " " " " " " " " " " " " ①

    式中:N為同一樣品中浮游植物總個體數(shù),Ni為第i種浮游植物個體數(shù),Y為優(yōu)勢度,fi為第i種浮游植物在樣品中出現(xiàn)的頻度,Y≥0.02的物種定義為優(yōu)勢種(章宗涉和黃祥飛,1991)。

    采用α多樣性指數(shù)Shannon-Wiener多樣性指數(shù)(Shannon)、生態(tài)優(yōu)勢度指數(shù)(Simpson)、豐富度指數(shù)(Margalef)對浮游植物群落特征進行分析(宋高飛等,2021)。

    利用Origin 2021和SPSS statistics 26.0對浮游植物各類群密度數(shù)據(jù)和水體參數(shù)進行了處理和制圖。使用one-way ANOVA進行統(tǒng)計分析。統(tǒng)計數(shù)據(jù)除pH外其他環(huán)境因子及物種矩陣均經(jīng)過lg(x+1)轉換及其Z標準化。

    采用Canoco 5.0軟件對浮游植物物種數(shù)據(jù)與環(huán)境因子數(shù)據(jù)進行排序,并通過前項選擇(forward selection)挑選具有顯著性的因子,通過物種數(shù)據(jù)進行去趨勢對應分析(DCA),確定排序模型;通過蒙特卡洛檢驗分析環(huán)境變量解釋顯著性;通過物種與環(huán)境因子的相關性確定兩者關系(賴江山,2013)。

    2" "結果與分析

    2.1" "浮游植物群落特征

    2.1.1" "物種組成與優(yōu)勢種" 調查期間共鑒定到浮游植物8門79屬119種,以綠藻為主,共有59種,占總種類數(shù)的比例為49.58%;其次為硅藻,有32種,所占比例為26.89%;藍藻10種,所占比例為8.41%;隱藻、甲藻、裸藻、金藻、定鞭藻種類較少,所占比例分別為5.04%、4.20%、2.52%、2.52%和0.84%(圖2)。浮游植物優(yōu)勢種6種,隸屬于3門5屬,為蹄形藻、塔胞藻、小環(huán)藻、變異直鏈藻、隱藻和卵形隱藻(表1)。

    從季節(jié)上看,浮游植物優(yōu)勢種類數(shù)排序為秋季gt;夏季=春季gt;冬季。夏秋兩季綠藻種類最多,硅藻、藍藻次之,春季種類數(shù)排序為硅藻gt;綠藻gt;甲藻,冬季種類數(shù)排序為硅藻gt;綠藻gt;甲藻=隱藻。其中,春季優(yōu)勢種為蹄形藻、小環(huán)藻和隱藻;夏季為小環(huán)藻、隱藻和變異直鏈藻;秋季為小環(huán)藻、變異直鏈藻、塔胞藻和卵形隱藻;冬季為小環(huán)藻和變異直鏈藻,其中,小環(huán)藻在春季、秋季和冬季優(yōu)勢度最大,夏季以變異直鏈藻主(表1);向家壩水庫浮游植物在種類數(shù)上主要為綠藻-硅藻-藍藻型(圖2),而在豐度組成上則為硅藻-綠藻-隱藻型(圖3)。

    2.1.2" "α多樣性指數(shù)" 如圖4所示,向家壩水庫浮游植物Margalef指數(shù)、Shannon-Wiener多樣性指數(shù)和均勻性指數(shù)季節(jié)均值分別為:Margalef指數(shù)(春季:0.67、夏季:0.54、秋季:0.40、冬季:0.24);Shannon-Wiener多樣性指數(shù)(春季:1.50、夏季:0.90、秋季:1.15、冬季:0.66);均勻性指數(shù)(春季:0.63、夏季:0.38、秋季:0.56、冬季:0.33)。Margalef指數(shù)均值在時間尺度上表現(xiàn)為春季gt;夏季gt;秋季gt;冬季的變化趨勢,其余指數(shù)均值季節(jié)變化為春季gt;秋季gt;夏季gt;冬季;且各指數(shù)均值在季節(jié)變化上差異顯著(Plt;0.05)。

    2.2" "細胞密度與葉綠素a濃度的時空動態(tài)

    浮游植物細胞密度為0.30×105~127.26×105個/L,平均密度為(10.23±7.39)×105個/L,季節(jié)密度大小表現(xiàn)為夏季gt;春季gt;秋季gt;冬季(圖5a),單因素方差分析顯示夏季顯著高于其他季節(jié)(Plt;0.05),特別值得注意的是,在冬季壩前水域(XJB01)浮游植物細胞密度明顯高于其他季節(jié)同點位的細胞密度,可能原因為冬季枯水期水位低,水體流速緩慢,營養(yǎng)鹽沉積,導致該區(qū)域浮游植物易于在此區(qū)域富集。在空間尺度上的變化為沿程干流(XJB01~08)→支流(XB01~12)各樣點浮游植物細胞密度呈現(xiàn)先增加后降低再增加再降低的“M”型變化趨勢,冬季變化不明顯,趨近于直線,在支流(XB02~11)處均值較高,在春季XB07、夏季XB11和秋季XB09處達到峰值(圖5b)。

    Chl-a濃度全年變化范圍為0.12~22.46 μg/L,平均值為(3.18±1.59)μg/L,其中春季gt;夏季gt;秋季gt;冬季,單因素方差分析顯示4個季節(jié)間具有顯著差異(Plt;0.05)(圖5c)。在空間上表現(xiàn)為,除冬季變化不明顯外,春季在XJB02、XB01、XB03和XB08出現(xiàn)峰值,夏季在XB03和XB12處出現(xiàn)峰值,秋季在XJB02、XB02和XB09處出現(xiàn)峰值,可以看出除冬季外支流(XB01~XB12)的細胞密度和葉綠素a濃度明顯高于干流(XJB01~XJB08),具有明顯的時空異質性(圖5d)。

    2.3" "環(huán)境因子

    在調查期間向家壩水庫水溫為12.96~27.40 ℃,流速為0.10~5.60 m/s,NO3--N濃度為0.40~2.22 mg/L,NH4+-N濃度為0~0.25 mg/L,TP濃度為0~0.67 mg/L,TN濃度為0.52~2.54 mg/L,PO43--P濃度為0~0.08 mg/L,CODMn濃度為0.40~2.34 mg/L,Chl-a濃度為0.12~22.46 μg/L,DO濃度為7.02~12.82 mg/L,SPC為237.6~543.00 μS/cm,pH為7.32~9.69,均值為8.58,水體整體偏堿性,ORP為42.00~306.70 mV,Tur為0.30~77.60,SD為0.35~4.90 m(圖6)。TN平均濃度為夏季gt;冬季gt;春季gt;秋季,TP平均濃度為春季gt;夏季gt;冬季gt;秋季,NO3--N平均濃度為夏季gt;春季gt;冬季gt;秋季,NH4+-N和pH的平均值表現(xiàn)為春季gt;冬季gt;夏季gt;秋季,CODMn平均濃度為春季gt;秋季gt;冬季gt;夏季,Chl-a平均濃度為春季gt;夏季gt;秋季gt;冬季,WT、Tur和Vel的平均值表現(xiàn)為夏季gt;秋季gt;春季gt;冬季,DO平均濃度為春季gt;冬季gt;秋季gt;夏季,SPC平均值表現(xiàn)為夏季gt;冬季gt;秋季gt;春季,ORP平均值表現(xiàn)為秋季gt;冬季gt;夏季gt;春季,SD平均值表現(xiàn)為冬季gt;春季gt;秋季gt;夏季(表2)。單因素方差分析表明,TN、TP和SPC除夏季和冬季外,其余季節(jié)間均有顯著差異(Plt;0.05),NO3--N和WT 在四季之間均有顯著差異(Plt;0.05),NH4+-N和SD除春季和冬季外,其余季節(jié)間均有顯著差異(Plt;0.05),PO43--P和Chl-a除春季和夏季外,其余季節(jié)間均有顯著差異(Plt;0.05),CODMn春季和夏季、夏季和秋季、冬季間具有顯著差異(Plt;0.05),DO春季和夏季、秋季、冬季,秋季和冬季之間有顯著差異(Plt;0.05),ORP春季和秋季、冬季,夏季和秋季、冬季有顯著差異(Plt;0.05),秋季和冬季之間差異不顯著,Tur和pH除夏季和秋季外,其余季節(jié)間均有顯著差異(Plt;0.05),Vel除春季與冬季、夏季與秋季外,其余季節(jié)間均有顯著差異。干支流TP和PO43--P的平均濃度相同,CODMn、SPC、Tur、SD和Vel顯示為干流gt;支流,其中Vel具有極顯著差異(Plt;0.01),Tur具有顯著差異(Plt;0.05),TN、NO3--N、NH4+-N、WT、DO和ORP顯示為支流gt;干流,其中DO具有顯著差異(Plt;0.05),各環(huán)境因子隨季節(jié)和樣點變化趨勢如圖6所示。

    2.4" "浮游植物群落與環(huán)境因子的排序分析

    分別以4月、7月、10月和1月數(shù)據(jù)代表春、夏、秋、冬季,對向家壩水庫4個季節(jié)的14個環(huán)境變量與對應月份的浮游植物總細胞密度及優(yōu)勢種細胞密度進行冗余分析(RDA),研究向家壩水庫浮游植物群落的分布規(guī)律,并探究驅動浮游植物群落時空格局發(fā)生動態(tài)變化的環(huán)境決定因子。

    DCA分析結果顯示排序軸最長軸長度梯度lt;3,表明采用線性模型RDA能更好地解釋本研究排序分析。RDA排序分析(圖7)結果顯示,第一排序軸和第二排序軸共解釋了66.9%浮游植物群落與環(huán)境因子的關系變量,說明本研究中的排序分析較好地解釋了浮游植物群落與環(huán)境因子的關系。使用Monte carlo置換檢驗前向選擇方法對影響浮游植物群落演替的環(huán)境變量進行識別和篩選,發(fā)現(xiàn)營養(yǎng)鹽、透明度、電導率和水溫是驅動向家壩水庫浮游植物群落時空分布的環(huán)境決定因子,其次是流速、濁度。浮游植物總細胞密度與水溫和氮磷營養(yǎng)鹽呈顯著正相關,與電導率、透明度呈顯著負相關,與流速呈負相關;小環(huán)藻作為四季優(yōu)勢種,與NO3--N、TN、PO43--P、SPC和WT呈顯著負相關,與SD呈顯著正相關,與Vel和Tur呈負相關;變異直鏈藻是夏、秋、冬3季的優(yōu)勢種,與SD和SPC呈顯著正相關,與WT、NO3--N、TN和PO43--P呈顯著負相關,與Vel呈正相關;蹄形藻和隱藻同為春季優(yōu)勢種,與NO3--N、TN和PO43--P呈顯著正相關,與SD和SPC呈顯著負相關,與Vel和Tur呈負相關;塔胞藻和卵形隱藻同為秋季優(yōu)勢種,塔胞藻表現(xiàn)為與SD和營養(yǎng)鹽呈正相關,與WT、Vel、Tur和SPC呈負相關,而卵形隱藻則與SD和Vel呈正相關,與SPC、WT、Tur、NO3--N、TN、PO43--P和TP呈負相關關系。樣方中秋季采樣點與冬季采樣點重合度較高,且優(yōu)勢種中小環(huán)藻為絕對優(yōu)勢種,表明秋冬兩季小環(huán)藻生態(tài)位相似度高,而夏季絕對優(yōu)勢種為變異直鏈藻,其生態(tài)位與秋冬季小環(huán)藻生態(tài)位具有較大差異。

    3" "討論

    3.1" "向家壩水庫水體環(huán)境因子特征

    河流生態(tài)系統(tǒng)作為地球物理化學循環(huán)的重要載體,具有輸水輸沙、供水自凈、航運發(fā)電和生態(tài)水文等多種功能(倪晉仁和劉元元,2006)。在調節(jié)氣候、改善生態(tài)環(huán)境的同時,為維持系統(tǒng)內(nèi)的生物群落和棲息環(huán)境提供了重要保障。而在金沙江下游修建了大壩后,就從天然的河流生態(tài)系統(tǒng)轉變?yōu)榫哂泻恿?湖泊特征的復合生態(tài)系統(tǒng),從而使向家壩水庫具有了河流和湖泊的復雜生境,使得水域生態(tài)系統(tǒng)在物質循環(huán)和能量流動的生態(tài)系統(tǒng)功能發(fā)生了改變。

    水體理化因子(水溫、pH、溶解氧、電導率、營養(yǎng)鹽等)和水文條件是表征水環(huán)境特征、水生態(tài)結構功能的基礎性指標(張斌,2019),為研究向家壩水庫水體環(huán)境提供了直接的證據(jù)。有關研究(王耀耀,2019;米武娟,2022)表明,金沙江TN形態(tài)以NO3--N為主,占TN濃度的67.3%~91.8%。向家壩水庫為中營養(yǎng)型水庫(陳興,2020)。周圍多礦區(qū),氮、磷主要為外源輸入,具有滯留效應。向家壩水庫在蓄水運行過程中有大量泥沙沉積,成為氮、磷營養(yǎng)鹽的主要儲蓄庫。在蓄水前后,當水體的水環(huán)境發(fā)生改變,庫區(qū)底泥中的部分氮、磷會重新釋放到水體中(劉尚武等,2019;米武娟等,2022)。庫區(qū)中水的pH呈中性偏弱堿性,除個別點外,全年變化不大,與本文研究一致。ORP表現(xiàn)為還原性,對磷營養(yǎng)鹽的釋放具有促進作用(米武娟等,2022),表現(xiàn)出顯著的季節(jié)差異性。SD、Tur季節(jié)變化明顯,Vel除個別點外全年變化不明顯,三者具有相關性,表現(xiàn)為流速變大導致水體發(fā)生攪動使得濁度升高,透明度降低。TN、TP和氮磷營養(yǎng)鹽都具有季節(jié)變化特征。通過對相關文獻(張斌,2019;陳興,2020;梁秀蓮,2021)的對比,發(fā)現(xiàn)近年來水體理化指標的變化差異較小。因此,鑒于金沙江流域作為我國一個跨多省的大型水文系統(tǒng),而其參考研究數(shù)據(jù)較少,建議今后對其進行長期水質監(jiān)測,以保證實時掌控水庫的水質情況,防止水質惡化、富營養(yǎng)化的發(fā)生以及探究水庫大壩的修建對流域水體水質變化的影響狀況,為全球氣候環(huán)境變化提供科學依據(jù)。

    3.2" "向家壩水庫浮游植物群落結構時空變化

    浮游植物又稱微藻,是一類在水中以浮游方式生活的低等植物(劉建康,1999)。浮游植物是水生態(tài)系統(tǒng)食物網(wǎng)的基礎,在生態(tài)系統(tǒng)平衡中起到重要作用,作為食物鏈的底層,為其他生物提供食物。由于個體較小、生活史較短,浮游植物在一個周期內(nèi)可以繁殖多代,因而可以迅速地響應水體環(huán)境變化,不同生境中浮游植物的群落組成具有一定的差異(Reynolds,1984)。本研究中,向家壩水庫共鑒定浮游植物119個分類單位,其中綠藻門最多,共59種,占浮游植物種數(shù)49.58%,其次為硅藻門,共32種,占比為26.89%,再次為藍藻門共10種,所占比例為8.41%,其他藻類如隱藻、甲藻、裸藻、金藻、定鞭藻也有,但占比不大。研究期間向家壩水庫浮游植物群落生態(tài)類型呈綠藻-硅藻-藍藻型。

    向家壩地區(qū)屬于熱帶季風氣候(李清雪等,2022),年均水溫為19.48 ℃,夏季高溫多雨,且降雨集中。一般認為,硅藻門種類適合生活在早春深秋低溫的環(huán)境,該環(huán)境下某些藻類如綠藻門和藍藻門的部分種類生長會受到抑制;夏季隨著水溫的升高,綠藻門和藍藻門的種類和數(shù)量有所上升。除此之外,其他門類的藻類在群落構建過程中也起到一定的貢獻作用。從時間看,向家壩水庫浮游植物種類、細胞密度在季節(jié)上也有不同差異。其中春季藻類種類最多,冬季種類最少,細胞密度最高值出現(xiàn)在夏季,其中硅藻在四季中都為優(yōu)勢種且優(yōu)勢度最大。春季為枯水期,當外源輸入減少,水位降低時,硅藻的生物量得到明顯增加(Sarmento et al,2006)。研究表明,自然水體中小環(huán)藻數(shù)量與流速成反比(李波等,2022),春季屬于枯水期,枯水期泥沙含量降低、透明度升高,有利于浮游植物的光合作用。水體流動緩慢,易于藻類的集聚,提高現(xiàn)存量(王敏等,2016),這與本文研究相一致。研究表明向家壩水庫為中營養(yǎng)型水庫(陳興,2020),小環(huán)藻最適宜生境為β-中污染水體,充足的營養(yǎng)鹽使之成為小環(huán)藻的適宜生境。此外,小環(huán)藻與水溫呈顯著的負相關關系,低水溫的持續(xù)時間越長,越利于小環(huán)藻的發(fā)生(鄧春光,2007)。在本研究中,小環(huán)藻在除夏季外的其他季節(jié)中均顯示出最大的優(yōu)勢度,特別是在冬季,浮游植物的細胞密度降至最低,原因為冬季水溫的降低,導致浮游植物細胞密度減少。然而,即便在這樣的條件下,小環(huán)藻依然保持了其最大的優(yōu)勢度,這表明小環(huán)藻對低溫環(huán)境具有較強的適應能力。

    3.3" "浮游植物群落結構與環(huán)境因子的相關性

    浮游植物群落結構分布特征與環(huán)境因子及其本身的生存策略密切相關,隨不同季節(jié)和流域內(nèi)水文條件的變化而變化,浮游植物群落能夠通過調控演替機制而維持群落的生態(tài)平衡,如浮游植物的形態(tài)結構、沉降情況、光合能力、固碳作用和對營養(yǎng)鹽的利用等生態(tài)響應機制(Reynolds et al,1993;馬煜等,2021)。影響浮游植物群落結構的環(huán)境因素主要包括水溫、光照、pH、溶解氧、營養(yǎng)鹽與有機物含量等。不同物種對于不同環(huán)境因子的適應性差異會造成浮游植物群落發(fā)生演替。浮游植物的生長都有對應的耐受溫度范圍,超過這個范圍初級生產(chǎn)力將會下降,甚至無法生存,一般情況下適宜的水溫區(qū)間為10~40 ℃(方麗娟等,2014)。水溫還能通過改變浮游動物的捕食壓力和競爭策略來改變食物的質量和閾值,從而間接影響藻類的群落結構(楊位迪,2007)。光在時間和空間上具有異質性,并且能夠被介質、可溶性懸浮顆粒、浮游植物等吸收,因此可以調控水體中的浮游植物垂直遷移(黃志敏等,2014)。氮、磷等營養(yǎng)鹽與微量元素是浮游植物的基礎資源,如硅藻還需要硅元素。有研究表明,磷是限制水生態(tài)系統(tǒng)初級生產(chǎn)力的主要影響因素,多數(shù)情況下水生態(tài)系統(tǒng)并不缺少氮,過量磷的輸入會造成富營養(yǎng)化(Schindler et al,2008)。

    將物種數(shù)據(jù)的變化分解為與環(huán)境變量有關的方差,用來探索群落物種組成受環(huán)境變量約束的關系(張金屯,2004;Lep? amp; ?milauer,2006)。RDA分析可同時結合多個環(huán)境因子,結果直觀且信息量大,分析效果明顯。本文RDA分析結果表明,影響向家壩水庫浮游植物群落時空分布的環(huán)境決定因子為WT、SPC和NO3--N,其次為流速和濁度。在向家壩水庫研究期間,WT、SPC和NO3--N四季間具有顯著性差異(Plt;0.05),WT具有明顯的時空異質性,NO3--N除秋季外其余季節(jié)間在空間上也存在差異性,SPC除夏季和冬季外,其余季節(jié)間均有顯著差異。說明水溫、電導率及氮磷營養(yǎng)鹽含量為驅動向家壩水庫浮游植物群落動態(tài)分布格局的環(huán)境決定變量。因此,在今后對向家壩水庫的管理中,為防止水質惡化、水體富營養(yǎng)化以及水華的發(fā)生,應管控好氮磷營養(yǎng)鹽的外源性輸入,并充分利用水動力條件(流速、流量和水位等)對浮游植物進行科學調控。

    參考文獻

    陳興, 2020. 水庫“藻-菌”群落變化及其對水體CO2代謝的影響研究[D]. 重慶:重慶交通大學.

    鄧春光, 2007. 三峽庫區(qū)富營養(yǎng)化研究[M]. 北京:中國環(huán)境科學出版社.

    鄧堅, 2012. 中國內(nèi)陸水域常見藻類圖譜[M]. 武漢:長江出版社.

    方麗娟, 劉德富, 楊正健, 等, 2014. 水溫對浮游植物群落結構的影響實驗研究[J]. 環(huán)境科學與技術, 37(增刊2):45-50.

    FANG L J, LIU D F, YANG Z J, et al, 2014. Effects of water temperature on the phytoplankton community structure[J]. Environmental Science amp; Technology, 37(Suppl.2):45-50.

    國家環(huán)境保護總局, 2002. 水和廢水監(jiān)測分析方法[M].4版. 北京:中國環(huán)境科學出版社:243-285.

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

    黃志敏, 陳椽, 劉之威, 等, 2014. 貴州百花湖夏季浮游植物晝夜垂直分布特征[J]. 生態(tài)學報, 34(19):5389-5397.

    HUANG Z M, CHEN C, LIU Z W, et al,2014. Diel vertical distribution patterns of phytoplanktion in Baihua Reservoir, Guizhou Province[J]. Acta Ecologica Sinica, 34(19):5389-5397.

    賴江山, 2013. 生態(tài)學多元數(shù)據(jù)排序分析軟件Canoco5介紹[J]. 生物多樣性," 21(6):765-768.

    LAI J S, 2013. Canoco 5: a new version of an ecological multivariate data ordination program[J]. Biodiversity Science, 21(6): 765-768.

    李波, 張羽珩, 鄧睿, 等, 2022. 流速對梅尼小環(huán)藻生長與光合活性的影響[J]. 生態(tài)學雜志, 41(7):1385-1390.

    LI B, ZHANG Y H, DENG R, et al, 2022. Effects of water velocity on the growth and photosynthetic activity of Cyclotella meneghiniana (diatom)[J]. Chinese Journal of Ecology, 41(7): 1385-1390.

    李清雪, 靳慧慧, 趙海萍, 等, 2022. 向家壩庫區(qū)沉積物磷形態(tài)分布及釋放風險[J]. 環(huán)境科學學報, 42(9):182-190.

    LI Q X, JIN H H, ZHAO H P, et al, 2022. Distribution characteristics of sediments phosphorus species of Xiangjiaba Reservoir and release risk[J]. Acta Scientiae Circumstantiae, 42(9):182-190.

    梁秀蓮, 2021. 向家壩庫區(qū)底泥氮形態(tài)分布特征及微生物群落結構研究[D]. 邯鄲:河北工程大學.

    劉建康, 1999. 高級水生生物學[M]. 北京:科學出版社:176-198.

    劉尚武, 張小峰, 呂平毓, 等, 2019. 金沙江下游梯級水庫對氮、磷營養(yǎng)鹽的滯留效應[J]. 湖泊科學, 31(3):656-666.

    LIU S W, ZHANG X F, LYU P Y, et al, 2019. Effects of cascade reservoirs in the lower reaches of Jinsha River on nitrogen and phosphorus retention[J]. J. Lake Sci., 31(3):656-666.

    龍良紅, 徐慧, 紀道斌, 等, 2017. 向家壩水庫水溫時空特征及其成因分析[J]. 長江流域資源與環(huán)境, 26(5):738-746.

    LONG L H, XU H, JI D B, et al, 2017. Temporal and spatial characteristics of water Temperature in Xiangjiaba reservoir[J]. Resources and Environment in the Yangtze Basin, 26(5):738-746.

    米武娟, 張地繼, 許元釗, 等, 2022. 向家壩水庫沉積物理化特性及其營養(yǎng)釋放風險[J]. 河南師范大學學報(自然科學版), 50(5):117-125.

    MI W J, ZHANG D J, XU Y Z, et al, 2022. Physical and chemical characteristics and nutrient release risk of sediments in Xiangjiaba Reservoir[J]. Journal of Henan Normal University(Natural Science Edition), 50(5):117-125.

    馬煜, 陸欣鑫, 范亞文, 2021. 松花江哈爾濱段浮游植物群落格局及其與環(huán)境因子的相關性[J]. 生態(tài)學報, 41(1):224-234.

    MA Y, LU X X, FAN Y W, 2021. Correlation between phytoplankton community patterns and environmental factors in Harbin section of the Songhua River[J]. Acta Ecologica Sinica, 41(1):224-234.

    倪晉仁, 劉元元, 2006. 論河流生態(tài)修復[J]. 水利學報, 37(9):1029-1037.

    NI J R, LIU Y Y, 2006. Ecological rehabilitation of damaged river system[J]. Journal of Hydraulic Engineering, 37(9):1029-1037.

    宋高飛, 朱宇軒, Anila P Ajayan, 等, 2021. 重金屬污染區(qū)秋季浮游植物群落特征及其與環(huán)境因子的關系[J]. 水生生物學報, 45(5):1082-1092.

    SONG G F, ZHU Y X, ANILA P A, et al, 2021. The community structure of phytoplankton and its relationships with environmental factors of heavy metal pollution districts in autumn[J]. Acta Hydrobiologica Sinica, 45(5):1082-1092.

    王耀耀, 呂林鵬, 紀道斌, 等, 2019. 向家壩水庫營養(yǎng)鹽時空分布特征及滯留效應[J]. 環(huán)境科學, 40(8):3530-3538.

    WANG Y Y, LYU L P, JI D B, et al, 2019. Spatial and Temporal Distribution Characteristics and the Retention Effects of Nutrients in Xiangjiaba Reservoir[J]. Environmental Science, 40(8):3530-3538.

    王敏, 袁紹春, 徐煒, 等, 2016. 嘉陵江重慶出口段小環(huán)藻水華影響因素分析[J]. 重慶交通大學學報(自然科學版), 35(6):86-90.

    WANG M, YUAN S C, XU W, et al, 2016. Influence factor analysis on cyclotella bloom of Jialing river estuary in Chongqing[J]. Journal of Chongqing Jiaotong University(Natural Science), 35(6):86-90.

    楊位迪, 2007. 廈門港浮游動物對浮游植物的攝食壓力研究[D]. 廈門:廈門大學.

    張斌, 2019. 金沙江下游梯級水庫水體CH4、CO2分布與源匯特征研究[D]. 重慶:中國科學院大學(中國科學院重慶綠色智能技術研究院).

    張春梅, 朱宇軒, 宋高飛, 等, 2021. 南水北調中線干渠浮游植物群落時空格局及其決定因子[J]. 湖泊科學, 33(3):675-686.

    ZHANG C M, ZHU Y X, SONG G F, et al, 2021. Spatiotemporal pattern of phytoplankton community structure and its determining factors in the channel of the middle route of South-to-North Water Diversion Project[J]. J. Lake Sci., 33(3):675-686.

    張金屯, 2004. 數(shù)量生態(tài)學[M]. 北京:科學出版社.

    張忠波, 2014. 三峽與金沙江下游梯級水庫群發(fā)電優(yōu)化調度研究[D]. 天津:天津大學.

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

    HA K, KIM H W, JOO G J, 1998. The phytoplankton succession in the lower part of hypertrophic Nakdong River (Mulgum), South Korea[C]// Phytoplankton and Trophic Gradients. 1998. Dordrecht: Springer Netherlands: 217-227.

    HERRERA-SILVEIRA J A, MORALES-OJEDA S M, 2009. Evaluation of the health status of a coastal ecosystem in southeast Mexico: assessment of water quality, phytoplankton and submerged aquatic vegetation[J]. Marine Pollution Bulletin, 59(1/2/3):72-86.

    LEP? J, ?MILAUER P, 2006. Multivariate analysis of ecological data[J]. The Bulletin of the Ecological Society of America, 87(3):193.

    REYNOLDS C S, 2006. Ecology of phytoplankton[M]. Cambridge: Cambridge University Press.

    REYNOLDS C S, 1984. Phytoplankton periodicity: the interactions of form, function and environmental variability. Freshwater Biology,14(2):111-142.

    REYNOLDS C S, PADISáK J, SOMMER U, 1993. Intermediate disturbance in the ecology of phytoplankton and the maintenance of species diversity: a synthesis[J]. Hydrobiologia, 249:183-188.

    SARMENTO H, ISUMBISHO M, DESCY J P, 2006. Phytoplankton ecology of Lake Kivu (eastern Africa)[J]. Journal of Plankton Research, 28(9):815-829.

    SCHINDLER D W, HECKY R E, FINDLAY D L, et al, 2008. Eutrophication of lakes cannot be controlled by reducing nitrogen input: results of a 37-year whole-ecosystem experiment[J]. Proceedings of the National Academy of Sciences of the United States of America, 105(32):11254-11258.

    VAN GREMBERGHE I, VAN WICHELEN J, VAN DER GUCHT K, et al, 2008. Covariation between zooplankton community composition and cyanobacterial community dynamics in Lake Blaarmeersen (Belgium)[J]. FEMS Microbiology Ecology,63(2):222-237.

    VOGT R J, SHARMA S, LEAVITT P R, 2015. Decadal regulation of phytoplankton abundance and water clarity in a large continental reservoir by climatic, hydrologic and trophic processes[J]. Journal of Great Lakes Research, 41:81-90.

    (責任編輯" "熊美華)

    Spatio-temporal Characteristics and Environmental Driving Factors of the Phytoplankton Community in Xiangjiaba Reservoir

    TIAN Chuming1,2, ZHANG Diji3, MI Wujuan2, QI Qingsong2, XING Long3,SUN Zhifeng3, CAO Guangrong3, BI Yonghong2

    (1.The College of Life Sciences and Technology, Harbin Normal University, Key Laboratory of Aquatic Biodiversity Research in Heilongjiang Province, Harbin" "150025, P.R. China;

    2. State Key Laboratory of Fresh Water Ecology and Biotechnology, Institute of Hydrobiology,Chinese Academy of Sciences, Wuhan" "430072, P.R. China;

    3. River Basin Complex Administration Center, China Three Gorges, Corporation,Yichang" "443000, P.R. China)

    Abstract:Xiangjiaba reservoir is one of the four mega-reservoirs in the lower reaches of Jinsha River. In this study, we investigated phytoplankton community structure and water quality in Xiangjiaba reservoir, and analyzed the coupling of the phytoplankton community with environmental factors. From January 2019 to January 2020, phytoplankton and water environment surveys were carried out monthly at 20 sampling sites in the mainstream (XJB01-XJB08) and tributaries (XB01-XB12 ) of Xiangjiaba reservoir. The α-diversity index and redundancy analysis were used to analyze phytoplankton community characteristics and their relationship with environmental factors. A total of 119 phytoplankton species from 79 genera and 8 phyla were identified, with the dominance by Chlorophyta (59 species, 49.58%), followed by Bacillariophyta (32 species, 26.89%), Cyanobacteria (10 species, 8.41%), Cryptophytes (6 species, 5.04%), Dinoflagellates (5 species, 4.2%) Euglenophytes (3 species, 2.52%), Chrysophytes (3 species, 2.52%), and Haptophytes (1 species, 0.84%). There were 6 dominant species, belonging to 5 genera and 3 phyla, namely, Kirchneriella sp., Cyclotella sp., Cryptomonas sp., Melosira varians, Pyramimonas sp., and Cryptomonas ovata. The spatial and temporal differences of phytoplankton community density were significant, and the phytoplankton density in different seasons was in the order summer gt; spring gt; autumn gt; winter, and much higher in summer. Spatially, phytoplankton density in the tributaries was higher than that in the mainstream. Temporal differences in water quality parameters were obvious, with mean values in the following seasonal orders: water temperature, turbidity and flow velocity (summergt;autumngt;springgt;winter); total nitrogen (summergt;wintergt;springgt;autumn); total phosphorus and nitrate nitrogen (spring/summergt;autumn/winter); ammonia nitrogen, dissolved oxygen and pH (spring/wintergt;summer/autumn) and CODMn (springgt;autumngt;wintergt;summer). Spatially, the concentrations of CODMn, conductivity, turbidity, transparency, and flow velocity were higher in the main stream than in tributaries, with flow velocity showing extremely significant differences (Plt;0.01) and turbidity showing significant differences (Plt;0.05). On the other hand, total nitrogen, nitrate nitrogen, ammonia nitrogen, water temperature, dissolved oxygen, and oxidation-reduction potential (ORP)were higher in tributaries than in the main stream, with dissolved oxygen showing significant differences (Plt;0.05). One-way ANOVA showed significant spatial and temporal heterogeneity among environmental factors in Xiangjiaba reservoir (Plt;0.05). Redundancy analysis (RDA) shows that the primary factors affecting spatial and temporal variations within the phytoplankton community were water temperature, nutrients, transparency and conductivity, followed by flow velocity and turbidity. The results of this study provide a basis for scientific decision-making about reservoir ecological management and protection.

    Key words: phytoplankton; environmental factors; redundancy analysis; Xiangjiaba reservoir

    基金項目:中國長江三峽集團有限公司2019年向家壩水庫水華監(jiān)測項目(SXSN/4495)。

    作者簡介:田楚銘,1997年生,男,碩士研究生,專業(yè)方向為藻類生態(tài)學。E-mail:421454328@qq.com

    通信作者:畢永紅,1974年生,男,博士,研究員,研究方向為水域生態(tài)學。E-mail:biyh@ihb.ac.cn

    ponron亚洲| 人人妻人人看人人澡| 亚洲人成网站高清观看| 国产在视频线在精品| 深夜a级毛片| 精品国内亚洲2022精品成人| 国产精品亚洲av一区麻豆| 一个人观看的视频www高清免费观看| 黄片小视频在线播放| 男人舔女人下体高潮全视频| 激情在线观看视频在线高清| 少妇的逼水好多| 99久久精品国产亚洲精品| 可以在线观看的亚洲视频| 国产视频内射| 在线播放国产精品三级| 亚洲av五月六月丁香网| 亚洲av成人不卡在线观看播放网| 亚洲熟妇中文字幕五十中出| 亚洲精品久久国产高清桃花| 国产精品久久电影中文字幕| 国产探花极品一区二区| 日韩中字成人| 色播亚洲综合网| 性欧美人与动物交配| 麻豆一二三区av精品| 国产极品精品免费视频能看的| avwww免费| 性色av乱码一区二区三区2| 中文字幕精品亚洲无线码一区| 中文字幕久久专区| 成人特级黄色片久久久久久久| 日韩欧美精品v在线| 亚洲av日韩精品久久久久久密| 男人和女人高潮做爰伦理| 精品人妻一区二区三区麻豆 | 国产黄a三级三级三级人| 在线观看午夜福利视频| 在线播放国产精品三级| 国语自产精品视频在线第100页| 一级a爱片免费观看的视频| 舔av片在线| 男人和女人高潮做爰伦理| 老女人水多毛片| 日本成人三级电影网站| 97超视频在线观看视频| 亚洲欧美日韩高清专用| 天堂av国产一区二区熟女人妻| 久久香蕉精品热| 男女下面进入的视频免费午夜| 精品免费久久久久久久清纯| 久久久国产成人精品二区| 在线播放无遮挡| 精品一区二区三区视频在线观看免费| 亚洲av.av天堂| 久久天躁狠狠躁夜夜2o2o| 精品久久久久久成人av| 精品福利观看| 国产精品免费一区二区三区在线| 亚洲欧美日韩卡通动漫| 内地一区二区视频在线| 成年女人看的毛片在线观看| 亚洲av美国av| 男人和女人高潮做爰伦理| 国产在线精品亚洲第一网站| 精品国产亚洲在线| 国产成人福利小说| 亚洲欧美激情综合另类| 亚洲国产精品合色在线| 国产精品,欧美在线| 亚洲精华国产精华精| 国产视频内射| 99久久九九国产精品国产免费| 成人三级黄色视频| 毛片一级片免费看久久久久 | 一区二区三区免费毛片| 在线免费观看的www视频| 久久久久免费精品人妻一区二区| 啦啦啦观看免费观看视频高清| 亚洲一区二区三区不卡视频| 欧美区成人在线视频| 亚洲经典国产精华液单 | 欧美成人性av电影在线观看| 最近在线观看免费完整版| 婷婷丁香在线五月| 日本一本二区三区精品| 内射极品少妇av片p| 国产午夜福利久久久久久| 国产三级在线视频| 日韩欧美国产一区二区入口| 色噜噜av男人的天堂激情| 天堂√8在线中文| 全区人妻精品视频| 国产aⅴ精品一区二区三区波| 看黄色毛片网站| 亚洲av五月六月丁香网| 一个人看的www免费观看视频| 日韩欧美在线乱码| 欧美日韩国产亚洲二区| 欧美激情国产日韩精品一区| 亚洲 欧美 日韩 在线 免费| 自拍偷自拍亚洲精品老妇| 亚洲激情在线av| 极品教师在线免费播放| 国产成人影院久久av| 99热这里只有是精品在线观看 | 亚洲经典国产精华液单 | 欧美性猛交╳xxx乱大交人| 性色avwww在线观看| 很黄的视频免费| 国产探花在线观看一区二区| 国产爱豆传媒在线观看| 深夜精品福利| 国产私拍福利视频在线观看| 成人精品一区二区免费| 在线播放国产精品三级| 窝窝影院91人妻| 国产精品乱码一区二三区的特点| 身体一侧抽搐| 一区福利在线观看| or卡值多少钱| 俺也久久电影网| 成熟少妇高潮喷水视频| 天堂av国产一区二区熟女人妻| 白带黄色成豆腐渣| 国产精品久久久久久亚洲av鲁大| 亚洲av.av天堂| 天美传媒精品一区二区| 国产淫片久久久久久久久 | 丁香六月欧美| 日韩中字成人| av女优亚洲男人天堂| 看黄色毛片网站| 亚洲欧美日韩无卡精品| 日本 欧美在线| 国产高清三级在线| 国产午夜精品久久久久久一区二区三区 | 久久精品国产99精品国产亚洲性色| 少妇熟女aⅴ在线视频| 婷婷六月久久综合丁香| 欧美成狂野欧美在线观看| 丁香六月欧美| 无遮挡黄片免费观看| 99国产精品一区二区蜜桃av| 日本一本二区三区精品| 在线观看av片永久免费下载| 男女之事视频高清在线观看| 国产aⅴ精品一区二区三区波| 老熟妇仑乱视频hdxx| 久久久精品欧美日韩精品| 男人狂女人下面高潮的视频| 麻豆成人av在线观看| 一二三四社区在线视频社区8| 亚洲人成伊人成综合网2020| 欧美日韩瑟瑟在线播放| 国产成年人精品一区二区| 午夜激情福利司机影院| 变态另类成人亚洲欧美熟女| АⅤ资源中文在线天堂| 最新在线观看一区二区三区| 国产精品一区二区三区四区免费观看 | 欧美乱色亚洲激情| 99热只有精品国产| 日本黄色片子视频| 国产精品久久久久久久电影| 亚洲无线在线观看| 国内久久婷婷六月综合欲色啪| 成人特级黄色片久久久久久久| 色综合站精品国产| 直男gayav资源| 国产高清视频在线播放一区| 51午夜福利影视在线观看| 亚洲avbb在线观看| 波多野结衣高清无吗| 成人午夜高清在线视频| 国产在线男女| 99久久精品热视频| aaaaa片日本免费| 99久久九九国产精品国产免费| 男女床上黄色一级片免费看| 色视频www国产| 国语自产精品视频在线第100页| 长腿黑丝高跟| 亚洲国产精品久久男人天堂| 一区二区三区高清视频在线| 欧美成人性av电影在线观看| 精品国内亚洲2022精品成人| 国产不卡一卡二| 禁无遮挡网站| 亚洲欧美日韩东京热| 黄片小视频在线播放| 免费在线观看成人毛片| 全区人妻精品视频| av视频在线观看入口| 国产精品一区二区免费欧美| 精品一区二区三区人妻视频| 国产精品久久久久久亚洲av鲁大| 国产综合懂色| 韩国av一区二区三区四区| 欧美日韩黄片免| 久久久久久久久大av| 国产精品久久视频播放| 欧美成人免费av一区二区三区| 简卡轻食公司| 午夜精品一区二区三区免费看| 亚洲中文字幕日韩| 精品一区二区免费观看| 久久天躁狠狠躁夜夜2o2o| 欧洲精品卡2卡3卡4卡5卡区| 黄色一级大片看看| 午夜福利18| 精品久久国产蜜桃| 久久性视频一级片| 淫妇啪啪啪对白视频| 小说图片视频综合网站| 日本免费a在线| 国产大屁股一区二区在线视频| 日韩精品青青久久久久久| 国产高清有码在线观看视频| 国产欧美日韩精品一区二区| 精品久久久久久久久亚洲 | 最近在线观看免费完整版| 欧美激情国产日韩精品一区| 国产精品一区二区三区四区久久| 国产精品伦人一区二区| 亚洲 国产 在线| 国产精品久久久久久人妻精品电影| 午夜精品一区二区三区免费看| 三级国产精品欧美在线观看| 亚洲一区二区三区色噜噜| 高潮久久久久久久久久久不卡| 91麻豆av在线| 欧美高清性xxxxhd video| 美女被艹到高潮喷水动态| 99视频精品全部免费 在线| 国内揄拍国产精品人妻在线| 欧美zozozo另类| 精品日产1卡2卡| 欧美日韩瑟瑟在线播放| 中文字幕人成人乱码亚洲影| 久久6这里有精品| 亚洲国产精品久久男人天堂| 成年版毛片免费区| 久久精品国产清高在天天线| 国产精品久久视频播放| 久久人妻av系列| 999久久久精品免费观看国产| 亚洲av免费在线观看| 亚洲欧美日韩高清专用| 中文字幕高清在线视频| 神马国产精品三级电影在线观看| 露出奶头的视频| 亚洲一区高清亚洲精品| 免费无遮挡裸体视频| 一二三四社区在线视频社区8| 亚洲av电影不卡..在线观看| av天堂在线播放| 桃红色精品国产亚洲av| 看十八女毛片水多多多| 中文字幕熟女人妻在线| 国产aⅴ精品一区二区三区波| 精品人妻视频免费看| 麻豆国产av国片精品| 激情在线观看视频在线高清| av在线天堂中文字幕| 久久久久久久亚洲中文字幕 | 亚洲国产色片| 亚洲成人免费电影在线观看| 国产高清有码在线观看视频| 国产精品不卡视频一区二区 | 搡老熟女国产l中国老女人| 日本精品一区二区三区蜜桃| 99热这里只有精品一区| 日本与韩国留学比较| 日韩亚洲欧美综合| 日韩中文字幕欧美一区二区| 一个人观看的视频www高清免费观看| netflix在线观看网站| 国产精品免费一区二区三区在线| 色综合欧美亚洲国产小说| 久久精品人妻少妇| 真人一进一出gif抽搐免费| 日韩国内少妇激情av| 国产亚洲欧美在线一区二区| 九色国产91popny在线| 色噜噜av男人的天堂激情| 老司机午夜十八禁免费视频| 91久久精品电影网| 久9热在线精品视频| 一级黄色大片毛片| 夜夜爽天天搞| 欧美不卡视频在线免费观看| 亚洲性夜色夜夜综合| 啦啦啦韩国在线观看视频| 久久99热这里只有精品18| 免费看日本二区| 在线观看av片永久免费下载| 国产亚洲精品久久久com| 欧美不卡视频在线免费观看| aaaaa片日本免费| 三级男女做爰猛烈吃奶摸视频| 免费看美女性在线毛片视频| 国产av麻豆久久久久久久| 日韩欧美三级三区| 最近最新中文字幕大全电影3| 国产不卡一卡二| 色视频www国产| 亚洲专区中文字幕在线| 给我免费播放毛片高清在线观看| 九九在线视频观看精品| 久久天躁狠狠躁夜夜2o2o| 日本熟妇午夜| 亚洲熟妇中文字幕五十中出| 一边摸一边抽搐一进一小说| 国产三级中文精品| 我要看日韩黄色一级片| 男人和女人高潮做爰伦理| 国产人妻一区二区三区在| 午夜福利成人在线免费观看| 国产探花极品一区二区| 欧美在线黄色| 性色avwww在线观看| 亚洲男人的天堂狠狠| 757午夜福利合集在线观看| .国产精品久久| 亚洲av.av天堂| 性欧美人与动物交配| 精品一区二区三区视频在线观看免费| 99国产综合亚洲精品| 国产高清三级在线| 我的女老师完整版在线观看| 成年女人毛片免费观看观看9| 白带黄色成豆腐渣| 欧美午夜高清在线| 国产欧美日韩一区二区精品| x7x7x7水蜜桃| 国产成人欧美在线观看| 啦啦啦韩国在线观看视频| 在线观看66精品国产| 欧美bdsm另类| 亚洲欧美日韩卡通动漫| 亚洲不卡免费看| 亚洲电影在线观看av| 久久久精品欧美日韩精品| 尤物成人国产欧美一区二区三区| 色哟哟哟哟哟哟| 90打野战视频偷拍视频| 99久久成人亚洲精品观看| av在线观看视频网站免费| 日韩中字成人| 午夜亚洲福利在线播放| 久久久色成人| 亚洲午夜理论影院| 国产一区二区在线观看日韩| 久久性视频一级片| 国产精品人妻久久久久久| 亚洲精品成人久久久久久| 一本久久中文字幕| 美女被艹到高潮喷水动态| 亚洲成人精品中文字幕电影| 女人十人毛片免费观看3o分钟| 怎么达到女性高潮| 香蕉av资源在线| 蜜桃久久精品国产亚洲av| 国产亚洲精品久久久com| 亚洲国产精品合色在线| 久久亚洲真实| 国语自产精品视频在线第100页| а√天堂www在线а√下载| 亚洲国产日韩欧美精品在线观看| 一个人免费在线观看的高清视频| 午夜日韩欧美国产| 丰满乱子伦码专区| 中文字幕av在线有码专区| 国产高清视频在线观看网站| 精品久久久久久久久久久久久| avwww免费| 色综合欧美亚洲国产小说| 深夜a级毛片| 桃色一区二区三区在线观看| 免费看a级黄色片| 非洲黑人性xxxx精品又粗又长| 身体一侧抽搐| 色精品久久人妻99蜜桃| 夜夜看夜夜爽夜夜摸| 天天躁日日操中文字幕| 国产亚洲精品久久久久久毛片| 最近最新免费中文字幕在线| 免费av不卡在线播放| 91久久精品国产一区二区成人| 久久精品国产亚洲av涩爱 | 国产成人欧美在线观看| 五月伊人婷婷丁香| 91麻豆av在线| 美女免费视频网站| 欧美日韩亚洲国产一区二区在线观看| 午夜a级毛片| 一级av片app| 亚洲精品成人久久久久久| 老司机福利观看| 一级a爱片免费观看的视频| 欧美黄色淫秽网站| 在线观看一区二区三区| 国产精品永久免费网站| 他把我摸到了高潮在线观看| 亚洲第一电影网av| 国产精品,欧美在线| 国产大屁股一区二区在线视频| 久久精品国产清高在天天线| 国产在线精品亚洲第一网站| 天堂√8在线中文| 脱女人内裤的视频| 丰满乱子伦码专区| 特级一级黄色大片| 久久6这里有精品| 成人国产综合亚洲| 精品久久久久久久末码| 亚洲av二区三区四区| 国产精品爽爽va在线观看网站| 91麻豆av在线| 1000部很黄的大片| 日韩av在线大香蕉| 国产探花极品一区二区| 亚洲av熟女| 色5月婷婷丁香| 亚洲av电影在线进入| 国产亚洲精品av在线| 99国产综合亚洲精品| 天天一区二区日本电影三级| 69av精品久久久久久| 在线观看66精品国产| 国产亚洲av嫩草精品影院| 一区二区三区激情视频| 好男人在线观看高清免费视频| 悠悠久久av| 亚洲人成电影免费在线| 人人妻,人人澡人人爽秒播| 国产日本99.免费观看| 久久精品91蜜桃| 亚洲精华国产精华精| 又黄又爽又刺激的免费视频.| 熟女人妻精品中文字幕| 欧美绝顶高潮抽搐喷水| 日韩欧美精品v在线| 精品久久久久久久末码| 亚洲成av人片在线播放无| 午夜福利免费观看在线| 在线十欧美十亚洲十日本专区| 中亚洲国语对白在线视频| 女同久久另类99精品国产91| 欧美色视频一区免费| 国产精品影院久久| 老司机午夜福利在线观看视频| 欧美在线黄色| av在线老鸭窝| 国产三级黄色录像| 亚洲va日本ⅴa欧美va伊人久久| 丝袜美腿在线中文| 999久久久精品免费观看国产| 成人性生交大片免费视频hd| 免费看a级黄色片| 亚洲av电影不卡..在线观看| 国产91精品成人一区二区三区| 亚洲在线观看片| 国产三级黄色录像| 亚洲va日本ⅴa欧美va伊人久久| 一区二区三区免费毛片| 三级国产精品欧美在线观看| 国产中年淑女户外野战色| 亚洲自拍偷在线| 一卡2卡三卡四卡精品乱码亚洲| 丰满的人妻完整版| 看免费av毛片| 欧美性猛交╳xxx乱大交人| 日本熟妇午夜| 99riav亚洲国产免费| 国产免费一级a男人的天堂| 热99re8久久精品国产| 97人妻精品一区二区三区麻豆| 蜜桃久久精品国产亚洲av| 亚洲自偷自拍三级| 午夜日韩欧美国产| 色播亚洲综合网| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美清纯卡通| 国产乱人视频| 国产精品久久久久久久久免 | 中亚洲国语对白在线视频| 国产精品一区二区性色av| 精品无人区乱码1区二区| 成人三级黄色视频| 国产在视频线在精品| 国产一区二区亚洲精品在线观看| 熟妇人妻久久中文字幕3abv| 桃红色精品国产亚洲av| 97碰自拍视频| 亚洲美女搞黄在线观看 | 亚洲美女黄片视频| av黄色大香蕉| 99热只有精品国产| 一本精品99久久精品77| 国产高清激情床上av| 两人在一起打扑克的视频| 真实男女啪啪啪动态图| 午夜日韩欧美国产| 国产精品女同一区二区软件 | 国产真实乱freesex| 两个人的视频大全免费| 亚洲最大成人手机在线| 日本黄色视频三级网站网址| 亚洲精华国产精华精| 小说图片视频综合网站| 久久天躁狠狠躁夜夜2o2o| 黄色日韩在线| 69人妻影院| 老司机深夜福利视频在线观看| 淫秽高清视频在线观看| 免费观看精品视频网站| 久久久久久久久久黄片| 亚洲天堂国产精品一区在线| 男女做爰动态图高潮gif福利片| 校园春色视频在线观看| 国产主播在线观看一区二区| 亚洲人成网站在线播| 中文字幕人成人乱码亚洲影| 国产精品伦人一区二区| 色在线成人网| 久久久国产成人精品二区| 亚洲国产精品sss在线观看| 性欧美人与动物交配| 搡老岳熟女国产| 亚洲av不卡在线观看| 国产不卡一卡二| 国内精品久久久久精免费| www.色视频.com| 国产成人a区在线观看| 天堂√8在线中文| 天天一区二区日本电影三级| 麻豆国产av国片精品| 欧美在线黄色| 国产av麻豆久久久久久久| 狠狠狠狠99中文字幕| 国产极品精品免费视频能看的| 97超级碰碰碰精品色视频在线观看| 香蕉av资源在线| 91九色精品人成在线观看| 欧美xxxx黑人xx丫x性爽| 欧美又色又爽又黄视频| 日本三级黄在线观看| 精品久久久久久久久av| 看十八女毛片水多多多| 中文字幕精品亚洲无线码一区| 欧美bdsm另类| 波多野结衣巨乳人妻| 免费搜索国产男女视频| 日韩免费av在线播放| 欧美另类亚洲清纯唯美| 有码 亚洲区| 久久久久国产精品人妻aⅴ院| 日韩亚洲欧美综合| 五月伊人婷婷丁香| 国产伦精品一区二区三区四那| 欧美黑人巨大hd| 久久欧美精品欧美久久欧美| 如何舔出高潮| 亚洲国产欧美人成| 亚洲国产精品合色在线| 久久久久国内视频| 亚洲中文字幕日韩| 久久久色成人| 久久欧美精品欧美久久欧美| 国产精品伦人一区二区| 搡女人真爽免费视频火全软件 | 欧美日本亚洲视频在线播放| 精品日产1卡2卡| 久久99热这里只有精品18| 国产aⅴ精品一区二区三区波| 久久久久久久久久黄片| 成人一区二区视频在线观看| 淫秽高清视频在线观看| 黄色配什么色好看| 禁无遮挡网站| 精品人妻偷拍中文字幕| 亚洲一区高清亚洲精品| 搡女人真爽免费视频火全软件 | 深夜精品福利| 亚洲av中文字字幕乱码综合| 观看美女的网站| 伊人久久精品亚洲午夜| 免费人成视频x8x8入口观看| 精品一区二区三区视频在线观看免费| 免费在线观看亚洲国产| av在线观看视频网站免费| 校园春色视频在线观看| 精品午夜福利在线看| 成人欧美大片| 国产精品精品国产色婷婷| 熟妇人妻久久中文字幕3abv| 午夜老司机福利剧场| 国产免费一级a男人的天堂| 欧美黄色淫秽网站| 91麻豆精品激情在线观看国产| 男女那种视频在线观看| 成年女人永久免费观看视频| 国产欧美日韩精品一区二区| 日本免费a在线| 在线a可以看的网站| 真人一进一出gif抽搐免费| 99国产极品粉嫩在线观看| 成年女人毛片免费观看观看9| 色综合亚洲欧美另类图片| 怎么达到女性高潮| 乱人视频在线观看| 国产精品1区2区在线观看.| 久久九九热精品免费|