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

    云南5個典型高原湖泊沉水植物群落組成及其與環(huán)境因子的關(guān)系

    2024-01-01 00:00:00朱利明張道熙史方鄒曦張志永胡蓮
    水生態(tài)學(xué)雜志 2024年5期
    關(guān)鍵詞:環(huán)境因子空間分布沉水植物

    摘要:通過對云南中部5個典型高原湖泊沉水植物和水環(huán)境狀況進(jìn)行野外調(diào)查,研究該區(qū)域主要沉水植物群落組成及其與環(huán)境因子的關(guān)系,為云南高原湖泊沉水植被保護(hù)和管理提供理論參考。綜合營養(yǎng)狀態(tài)指數(shù)計算表明,所選取的5個湖泊中,異龍湖、杞麓湖和星云湖屬于富營養(yǎng)湖泊,撫仙湖和陽宗海屬于貧-中營養(yǎng)湖泊。調(diào)查共發(fā)現(xiàn)沉水植物6科10屬15種,微齒眼子菜(Potamogeton maackianus)和穿葉眼子菜(P. perfoliatus)在貧-中營養(yǎng)湖泊中為主要優(yōu)勢種,穗花狐尾藻(Myriophyllum spicatum)和篦齒眼子菜(P. pectinatus)在各營養(yǎng)類型湖泊中均能成為優(yōu)勢種或亞優(yōu)勢種,輪葉黑藻(Hydrilla verticillata)、菹草(P. crispus)和金魚藻(Ceratophyllum demersum)僅在富營養(yǎng)湖泊中成為優(yōu)勢種。沉水植物多樣性分析表明,撫仙湖的香農(nóng)-威納多樣性指數(shù)(H')最高,平均值達(dá)到1.63;杞麓湖的H'最低,平均值僅0.79,總體來說,貧-中營養(yǎng)湖泊的沉水植物多樣性要高于富營養(yǎng)湖泊。沉水植物群落物種與水環(huán)境因子的冗余分析(redundancy analysis, RDA)結(jié)果顯示,透明度(SD)、pH、硝氮(NO3--N)是影響這些不同營養(yǎng)水平湖泊沉水植物分布特征的關(guān)鍵環(huán)境因子,pH、氨氮(NH4+-N)、水溫(WT)是影響貧-中營養(yǎng)湖泊沉水植物分布特征的關(guān)鍵環(huán)境因子,WT、電導(dǎo)率(SPC)、葉綠素a(Chl-a)是影響富營養(yǎng)湖泊沉水植物分布特征的關(guān)鍵環(huán)境因子。研究對于恢復(fù)、保護(hù)和管理該區(qū)域湖泊沉水植物具有重要指導(dǎo)意義。

    關(guān)鍵詞:高原湖泊;沉水植物;空間分布;環(huán)境因子;冗余分析

    中圖分類號:Q178.1" " " " 文獻(xiàn)標(biāo)志碼:A" " " " 文章編號:1674-3075(2024)05-0039-10

    沉水植物是淡水湖泊生態(tài)系統(tǒng)主要的初級生產(chǎn)者,在湖泊物質(zhì)循環(huán)和能量流動中發(fā)揮著重要作用,是維持湖泊清水穩(wěn)態(tài)的關(guān)鍵(Wang et al,2016)。沉水植物生命活動對湖泊環(huán)境因素的變化非常敏感,水體營養(yǎng)鹽濃度、底泥類型、水文水動力條件、濁度、透明度、溶解氧、水深等諸多環(huán)境因子都與沉水植物的群落組成與空間分布密切相關(guān)(Nakamura et al,2008;Bornette amp; Puijalon,2011;孔祥虹等,2015)。近些年來,隨著社會經(jīng)濟(jì)的快速發(fā)展,湖泊富營養(yǎng)化日益加重,我國多個湖泊沉水植物群落演替劇烈,群落組成和優(yōu)勢物種發(fā)生了明顯的改變,最終導(dǎo)致生物量銳減、多樣性減少,湖泊生態(tài)系統(tǒng)功能受到嚴(yán)峻挑戰(zhàn)(Sachse et al,2014)。

    湖泊富營養(yǎng)化會導(dǎo)致水生植物群落結(jié)構(gòu)發(fā)生變化,這與富營養(yǎng)水體本身對水生植物的脅迫、富營養(yǎng)促進(jìn)浮游植物與水生植物之間的競爭以及富營養(yǎng)促進(jìn)單一物種的大量繁殖擴(kuò)散等因素有關(guān)(Kalff,2002)。富營養(yǎng)化引起水生植物群落退化、生物多樣性降低的內(nèi)在機(jī)理是不同水生植物對水環(huán)境的適應(yīng)范圍存在差異(吳愛平等,2005;游灝,2006),耐受范圍小、生態(tài)位寬度窄的物種對水環(huán)境的變化更為敏感。有沉水植物存在的水域水質(zhì)通常較好。有研究表明,沉水植物具有富集水體污染物的功能,在富營養(yǎng)化水體修復(fù)過程中起到凈化水質(zhì)、維持水生態(tài)系統(tǒng)健康的作用(Carpenter amp; Lodge,1986)。當(dāng)前,國內(nèi)外眾多學(xué)者在沉水植被恢復(fù)重建領(lǐng)域開展了大量研究,光照強(qiáng)度(李文朝和連光華,1996)、營養(yǎng)鹽(Horppila amp; Nurminen,2005;金相燦等,2007)、底質(zhì)(謝貽發(fā)等,2007)、濁度(王文林等,2006)、溫度(Barko amp; Smart,1981;閆志強(qiáng)等,2014)等因素對沉水植物的生長至關(guān)重要,共同影響沉水植物的群落組成與分布(王華等,2008)。因此,研究沉水植物群落組成及分布格局與水環(huán)境因子的關(guān)系具有重要意義。

    云南是我國淡水湖泊資源最為集中的地區(qū)之一,其湖泊主要分布在斷裂帶或各大水系的分水嶺地帶,具有海拔較高、湖岸陡峻、面積小而水深的特點(diǎn)(曾理等,2009),兼具調(diào)節(jié)氣候、水量,提供生產(chǎn)生活用水、水產(chǎn)養(yǎng)殖產(chǎn)品,美化環(huán)境和發(fā)展旅游的功能(藍(lán)紅林,2001)。這些湖泊多為源頭型湖泊,其生態(tài)狀況同時影響著當(dāng)?shù)睾拖掠纹皆貐^(qū)的生態(tài)環(huán)境(楊柳和黃耀志,2002)。然而,隨著該地區(qū)工農(nóng)業(yè)的發(fā)展以及長期不合理的湖泊資源開發(fā),云南高原湖泊生態(tài)環(huán)境遭到嚴(yán)重破壞,富營養(yǎng)化嚴(yán)重,水體營養(yǎng)梯度日益明顯,生態(tài)系統(tǒng)日益退化,這一系列問題已引起眾多學(xué)者和社會的關(guān)注(顏昌宙等,2005)。以往的調(diào)查多由不同研究者單獨(dú)進(jìn)行,對多個湖泊同時期沉水植物分布情況了解不足。為了深入了解云南高原湖泊沉水植物群落組成及分布格局,本文對云南5個典型高原湖泊沉水植物進(jìn)行了系統(tǒng)調(diào)查,比較分析該區(qū)域沉水植物的區(qū)系特征,探明沉水植被現(xiàn)狀,對影響沉水植物分布的水環(huán)境因子進(jìn)行解析,辨識關(guān)鍵影響因子,以期為云南高原湖泊沉水植被保護(hù)和管理提供理論參考。

    1" "材料與方法

    1.1" "研究區(qū)域

    分別于2020年和2021年9-10月對云南中部5個典型高原湖泊進(jìn)行2個年度的沉水植物調(diào)查(圖1)。根據(jù)所選湖泊大小、受干擾程度以及沉水植物實(shí)際分布情況設(shè)置采樣點(diǎn),對沉水植物的群落組成及水環(huán)境因子進(jìn)行調(diào)查。其中,撫仙湖(FXH)和陽宗海(YZH)沉水植物主要局限于環(huán)湖湖岸帶一周,因此調(diào)查時沿湖岸帶設(shè)置采樣點(diǎn)(沈亞強(qiáng)等,2010;高弋明等,2021)。星云湖(XYH)、杞麓湖(QLH)和異龍湖(YLH)由于水生態(tài)環(huán)境較差,沉水植物分布面積和生物量已大為減少,并呈現(xiàn)向岸邊萎縮和分布不均勻的趨勢(沈亞強(qiáng)等,2010),同樣采用環(huán)湖設(shè)置采樣點(diǎn)的方法進(jìn)行調(diào)查。其中,撫仙湖設(shè)置4個采樣點(diǎn),星云湖、杞麓湖、陽宗海和異龍湖各設(shè)置3個采樣點(diǎn)。各湖泊的基本形態(tài)學(xué)特征見表1。

    1.2" "環(huán)境因子的測定

    現(xiàn)場采用透明度盤測定透明度(SD),用YSI多參數(shù)水質(zhì)監(jiān)測儀測定水體水溫(WT)、pH、電導(dǎo)率(SPC)、氧化還原電位(ORP)、溶解氧(DO)、鹽度(Sal)等參數(shù),采集表層水樣帶回實(shí)驗(yàn)室分析總氮(TN)、總磷(TP)、氨氮(NH4+-N)、硝氮(NO3--N)、磷酸鹽(PO43-)、葉綠素a(Chl-a)和高錳酸鹽指數(shù)(CODMn),水體化學(xué)指標(biāo)的測定參考《水和廢水監(jiān)測分析方法》(國家環(huán)境保護(hù)總局,2002)。

    1.3" "沉水植物的采集及處理

    深水區(qū)使用0.5 m×0.5 cm的采草器隨機(jī)采集3次;岸邊淺水區(qū)使用面積為 0.5 m×0.5 m 的鋁制樣方框確定采樣區(qū)域,徒手采集樣方框內(nèi)所有植物,采集3次。采集到的沉水植物及時沖洗干凈、分類、稱重,以濕重代表生物量。

    以往研究多使用植物株數(shù)密度來計算物種多樣性指數(shù),但沉水植物由于其特殊的生物學(xué)特性,絕大部分種類多通過分枝或分蘗的方式進(jìn)行繁殖,而不是簡單地通過增加株數(shù)來擴(kuò)大種群生物量。因此,本文參考Jeppesen等(2000)的方法,使用生物量來計算沉水植物的優(yōu)勢度和多樣性指數(shù),即用相對生物量(RB)與相對頻度(RF)來計算物種的優(yōu)勢度(DV):

    DV=(RB + RF)/2×100%" " " ①

    用香農(nóng)-威納多樣性指數(shù)(H')計算公式如下:

    H' = -[i=1Sbib lnbib]" " " " "②

    式中:bi為群落中i物種生物量;b為群落中所有沉水植物的生物量之和;S為該群落總物種數(shù)。

    1.4" "湖泊營養(yǎng)程度評價方法

    采用綜合營養(yǎng)指數(shù)(TLI)法評價湖泊富營養(yǎng)化(金相燦,1995;王明翠等,2002)。計算公式為:

    TLI(Chl-a)=10(2.5+1.086 ln Chl-a)" "③

    TLI(TP)=10(9.436+1.624 ln TP)" " "④

    TLI(TN)=10(5.453+1.694 ln TN)" " ⑤

    TLI(SD)=10(5.118-1.94 ln SD)" " ⑥

    TLI(∑)=W(Chl-a)·TLI(Chl-a)+W(TP)·TLI(TP)+W(TN)·TLI(TN)+W(SD)·TLI(SD)" " " " "⑦

    式中:W(Chl-a)、W(TP)、W(TN)和W(SD)分別為0.3260、0.2300、0.2192和0.2246。采用0~100的一系列連續(xù)數(shù)字對湖泊營養(yǎng)狀態(tài)進(jìn)行分級:TLI(∑)≤30為貧營養(yǎng);30lt;TLI(∑)≤50為中營養(yǎng);50lt;TLI(∑)≤60為輕度富營養(yǎng);60lt;TLI(∑)≤70為中度富營養(yǎng);70lt;TLI(∑)≤100為高度富營養(yǎng)。

    1.5nbsp; "數(shù)據(jù)分析

    利用One-way ANOVA對不同湖泊環(huán)境因子變量均值進(jìn)行分析(Plt;0.05),在進(jìn)行方差分析前,對數(shù)據(jù)進(jìn)行了方差齊性分析和正態(tài)分布分析,若方差齊性,選擇LSD檢驗(yàn),反之使用Tamhane's檢驗(yàn)(Plt;0.05),用SPSS19.0軟件進(jìn)行分析,origin 9.0進(jìn)行作圖。運(yùn)用R 4.1.3中的vegan包分析沉水植物優(yōu)勢度與環(huán)境因子的關(guān)系。在數(shù)據(jù)進(jìn)行排序分析前,對物種數(shù)據(jù)和環(huán)境因子數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理(ErO\"s et al,2012),再進(jìn)行去趨勢對應(yīng)分析(detrended correspondence analysis, DCA)以確定排序模型(Braak amp; Smilauer,2002)。采用方差膨脹因子(variance inflation factor,VIF)分析,剔除共線性強(qiáng)的環(huán)境因子(VIFgt;10);Pearson相關(guān)分析進(jìn)一步剔除強(qiáng)相關(guān)性環(huán)境因子(rgt;0.7);采用蒙特卡洛檢測方法檢驗(yàn)環(huán)境因子對植物影響的顯著性,保留Plt;0.05的環(huán)境因子(Blanchet et al,2008)。

    2" "結(jié)果與分析

    2.1" "水環(huán)境現(xiàn)狀

    不同營養(yǎng)類型湖泊水環(huán)境因子特征如圖2。從中可以看出,異龍湖、杞麓湖、星云湖和陽宗海4個湖泊的WT差異顯著(Plt;0.05),異龍湖最高,星云湖最低,撫仙湖與異龍湖WT具有顯著差異(Plt;0.05),其他湖泊差異不顯著(圖2A)。Sal、SPC、TDS 3個指標(biāo)在不同湖泊間具有顯著差異(Plt;0.05),表現(xiàn)為杞麓湖gt;星云湖gt;異龍湖gt;陽宗海gt;撫仙湖(圖2F、2G、2H)。陽宗海和撫仙湖的SD顯著大于杞麓湖、星云湖和異龍湖(P<0.05)(圖2B)。杞麓湖、星云湖和異龍湖TN、NH4+-N、TP顯著大于撫仙湖(Plt;0.05)(圖2J、2K、2M),pH、DO、ORP、PO43-、CODMn等水環(huán)境因子在各湖泊之間差異不大(圖2C、2D、2E、2N、2O)。綜合營養(yǎng)狀態(tài)指數(shù)法計算表明,異龍湖、杞麓湖和星云湖屬于富營養(yǎng)湖泊,營養(yǎng)指數(shù)分別為63.55±6.28、67.26±1.89和72.08±2.04;撫仙湖和陽宗海屬于貧-中營養(yǎng)湖泊,營養(yǎng)指數(shù)分別為28.01±4.82和41.47±1.96(圖3)。

    2.2" "沉水植物群落現(xiàn)狀

    本次共調(diào)查到沉水植物6科10屬15種,包括苦草(Vallisneria natans)、大茨藻(Najas marina)、篦齒眼子菜(Potamogeton pectinatus)、穿葉眼子菜(P. perfoliatus)、光葉眼子菜(P. lucens)、角果藻(Pteridophyta)、金魚藻(Ceratophyllum demersum)、麗藻屬(Nitella)、輪葉黑藻(Hydrilla verticillata)、輪藻屬(Chara)、馬來眼子菜(P. wrightii)、穗花狐尾藻(Myriophyllum spicatum)、微齒眼子菜(P. maackianus)、伊樂藻(Elodea nuttallii)和菹草(P. crispus)。

    云南高原湖泊中沉水植物的主要優(yōu)勢種為穗花狐尾藻、微齒眼子菜、輪葉黑藻、穿葉眼子菜、菹草、篦齒眼子菜和苦草。不同營養(yǎng)類型湖泊沉水植物優(yōu)勢種有所不同,微齒眼子菜和穿葉眼子菜為貧-中營養(yǎng)湖泊的優(yōu)勢種,穗花狐尾藻和篦齒眼子菜在各營養(yǎng)類型湖泊中均能成為優(yōu)勢種或亞優(yōu)勢種,輪葉黑藻、菹草和金魚藻僅在富營養(yǎng)湖泊中成為優(yōu)勢種(表2)。沉水植物多樣性分析表明,撫仙湖的香農(nóng)-威納多樣性H'最高,平均值達(dá)到1.63;杞麓湖的H'最低,平均值僅0.79??傮w來說,貧-中營養(yǎng)湖泊的沉水植物多樣性要大于富營養(yǎng)湖泊。在16個樣點(diǎn)中,出現(xiàn)頻率最高的是穗花狐尾藻,為87.5%;最低是大茨藻、角果藻和麗藻屬,為6.25%(圖4)。

    2.3" "沉水植物群落物種與環(huán)境因子的關(guān)系

    DCA結(jié)果表明,各排序軸的最大梯度長度(LOG)為:貧-中營養(yǎng)湖泊1.87、富營養(yǎng)湖泊3.67、整體湖泊2.99,都lt;4,符合線性模型,因此本研究統(tǒng)一選用RDA分析方法。分析顯示,貧-中營養(yǎng)湖泊、富營養(yǎng)湖泊、整體湖泊的第1、2軸共解釋物種與環(huán)境因子之間累計方差為86.34%、93.92%、94.61%,說明排序結(jié)果可以很好地解釋沉水植物群落組成與環(huán)境因子的關(guān)系(圖5)。

    蒙特卡洛置換檢驗(yàn)分析表明,貧-中營養(yǎng)湖泊中15個環(huán)境變量中只有3個變量對解釋沉水植物群落組成與環(huán)境因子的關(guān)系做出了顯著貢獻(xiàn)(即Plt; 0.05)。其中,穿葉眼子菜、輪葉黑藻、金魚藻與pH呈顯著正相關(guān),與NH4+-N呈顯著負(fù)相關(guān);菹草、角果藻、輪藻屬和微齒眼子菜與NH4+-N呈顯著正相關(guān),與pH呈顯著負(fù)相關(guān);馬來眼子菜和穗花狐尾藻與WT呈顯著正相關(guān),篦齒眼子菜、苦草和輪藻屬與WT呈顯著負(fù)相關(guān)(圖5a)。環(huán)境因子與坐標(biāo)軸的相關(guān)系數(shù)越高則表明環(huán)境因子對物種的相關(guān)性越強(qiáng),影響越大,由表3可以看出,pH、NH4+-N、WT是影響貧-中營養(yǎng)湖泊沉水植物分布特征的主要因子。

    富營養(yǎng)湖泊RDA分析中,通過蒙特卡洛置換檢驗(yàn)預(yù)選出符合條件的環(huán)境因子為WT、SPC、Chl-a,其中,輪葉黑藻、菹草和馬來眼子菜與SPC呈顯著負(fù)相關(guān),穗花狐尾藻和篦齒眼子菜與SPC呈顯著正相關(guān);金魚藻與Chl-a呈顯著正相關(guān),與WT呈顯著負(fù)相關(guān);穗花狐尾藻、馬來眼子菜、菹草與WT呈顯著正相關(guān),與Chl-a呈顯著負(fù)相關(guān)(圖5b)。環(huán)境因子與坐標(biāo)軸的相關(guān)系數(shù)表明,WT、SPC、Chl-a是影響富營養(yǎng)湖泊沉水植物分布特征的主要因子(表3)。

    整體湖泊RDA分析表明,SD、NO3--N、pH為蒙特卡洛置換檢驗(yàn)預(yù)選出的環(huán)境因子。其中,穿葉眼子菜、苦草、篦齒眼子菜、微齒眼子菜、輪藻屬、伊樂藻和麗藻屬分別與pH和SD呈顯著正相關(guān);菹草、輪葉黑藻、金魚藻以及馬來眼子菜分別與pH和SD呈顯著負(fù)相關(guān);穗花狐尾藻和馬來眼子菜與NO3--N呈顯著正相關(guān);穿葉眼子菜、苦草、篦齒眼子菜、微齒眼子菜、輪藻屬、伊樂藻、麗藻屬、金魚藻和輪葉黑藻都與NO3--N呈顯著負(fù)相關(guān)(圖5c)。環(huán)境因子與坐標(biāo)軸的相關(guān)系數(shù)表明,SD、NO3--N、pH是影響整體湖泊沉水植物分布特征的主要因子(表3)。

    3" "討論

    3.1" "沉水植物群落組成變化

    本研究通過綜合營養(yǎng)狀態(tài)指數(shù)法計算表明,異龍湖、杞麓湖和星云湖屬于富營養(yǎng)湖泊,撫仙湖和陽宗海屬于貧-中營養(yǎng)湖泊(圖3)。調(diào)查共發(fā)現(xiàn)沉水植物6科10屬15種,其中優(yōu)勢種包括穗花狐尾藻、微齒眼子菜、篦齒眼子菜、穿葉眼子菜、輪葉黑藻、菹草和苦草。微齒眼子菜和穿葉眼子菜在貧-中營養(yǎng)湖泊中可以成為優(yōu)勢種,穗花狐尾藻和篦齒眼子菜在各營養(yǎng)湖泊中均能成為優(yōu)勢種或亞優(yōu)勢種,輪葉黑藻、菹草和金魚藻僅能在富營養(yǎng)湖泊中成為優(yōu)勢種(表2)。導(dǎo)致這種分布現(xiàn)象的原因可能包括:(1)環(huán)境因素。已有研究表明水生植物的生長受水深、水溫、光照強(qiáng)度、營養(yǎng)鹽、水文水動力條件、底質(zhì)類型等眾多環(huán)境因子共同作用,水生生境條件對水生植物的物種組成及分布具有顯著的影響(秦伯強(qiáng),2002)。(2)人為因素。社會、經(jīng)濟(jì)、工農(nóng)業(yè)的發(fā)展以及長期以來對湖泊資源的不合理開發(fā),造成云南地區(qū)湖泊生態(tài)環(huán)境嚴(yán)重破壞,湖泊富營養(yǎng)化嚴(yán)重,湖泊生態(tài)系統(tǒng)日益退化(顏昌宙等,2005)。

    已有研究發(fā)現(xiàn),異龍湖水環(huán)境在21世紀(jì)90年代后逐漸惡化,2009年異龍湖水質(zhì)類別甚至達(dá)到劣Ⅴ類,水生生態(tài)系統(tǒng)發(fā)生劇烈變化,平均透明度為0.22 m,綜合營養(yǎng)狀態(tài)指數(shù)為78.13,同年進(jìn)行的該湖泊水生植物調(diào)查顯示沉水植物已完全消失(魏翔和唐光明,2014)。本研究調(diào)查異龍湖的綜合營養(yǎng)狀態(tài)指數(shù)為63.55±6.28,仍屬于中度富營養(yǎng)狀態(tài),沉水植物雖有所恢復(fù),但多以耐污種(穗花狐尾藻和篦齒眼子菜)為主,這與方馨等(2021)的研究結(jié)果一致。1970年代以來,受周邊圍湖造田以及農(nóng)業(yè)面源污染的長期影響,杞麓湖水生生態(tài)系統(tǒng)遭到嚴(yán)重破壞,對水質(zhì)要求非常嚴(yán)格的海菜花(Ottelia acuminata)、水鱉(Hydrocharis dubia)群落、荇菜(Nymphoides peltatum)群落等已相繼消失,對水質(zhì)要求較高的竹葉眼子菜群落分布面積巨減,穗花狐尾藻、菹草等耐污種群落逐漸增多,成為優(yōu)勢類群(沈亞強(qiáng)等,2010),本研究調(diào)查也發(fā)現(xiàn)穗花狐尾藻、篦齒眼子菜、菹草為杞麓湖的主要優(yōu)勢種。已有資料顯示,隨著星云湖周邊畜禽養(yǎng)殖業(yè)生產(chǎn)規(guī)模的擴(kuò)大,畜禽糞便污染已成為星云湖水環(huán)境的重要的污染源(邱成和葉春麗,2015)。另有研究認(rèn)為,星云湖周邊農(nóng)田施用的化學(xué)用劑污染及大量食草魚類放養(yǎng)等因素影響其水質(zhì)質(zhì)量,導(dǎo)致湖內(nèi)植物群落類型單一,僅有少數(shù)對環(huán)境要求較低的植物得以生存(李恒,1980)。本次調(diào)查發(fā)現(xiàn)星云湖綜合營養(yǎng)狀態(tài)指數(shù)為72.08±2.04,屬于高度富營養(yǎng)狀態(tài),間接導(dǎo)致該湖泊沉水植物生物多樣性較低,輪葉黑藻、金魚藻、菹草等耐污種為主要優(yōu)勢種(表2),這與楊加林等(2012)的研究結(jié)果一致。

    3.2" "沉水植物與環(huán)境因子的關(guān)系

    本研究結(jié)果發(fā)現(xiàn),在水體營養(yǎng)鹽水平差異不大的情況下(中營養(yǎng)湖泊和富營養(yǎng)湖泊),WT成為影響沉水植物分布的主要因子(圖5a和圖5b),且在富營養(yǎng)湖泊中,沉水植物物種多樣性隨著溫度的升高而增加(表2)。已有研究表明,水溫和沉積物溫度都會影響水生植物分布及其群落結(jié)構(gòu)特征(Rooney amp; Kalff,2000;Spencer et al,2000;Heino,2002)。溫度可對水體鹽度、營養(yǎng)鹽含量和氣體含量產(chǎn)生明顯的影響,在水生態(tài)系統(tǒng)的物理、化學(xué)和生物過程中發(fā)揮著重要作用(Welch,1952)。較高的水溫會增加水體的蒸發(fā)量,降低水體溶解氧,從而增加許多水生物種的含鹽量,增加其生存壓力(Al-Thahaibawi et al,2021)。本研究發(fā)現(xiàn)富營養(yǎng)湖泊異龍湖,其較高的水溫導(dǎo)致了較低的水體溶解氧和鹽度(圖2),這也是導(dǎo)致該湖泊沉水植物物種多以耐污種為主的原因。

    SPC是影響富營養(yǎng)湖泊沉水植物分布的另一主要環(huán)境因子(圖5b)。水體SPC的變化與水位、蒸發(fā)量都有顯著的關(guān)系(Al-Kenzawi et al,2011),此外,農(nóng)業(yè)活動面源污染也可能導(dǎo)致夏季水體SPC過高(Al-Shawi,2006)。已有研究表明,自1970年代以來,杞麓湖一直受周邊圍湖造田和農(nóng)業(yè)面源污染的影響,使得竹葉眼子菜群落分布面積銳減,穗花狐尾藻、菹草群落等耐污種逐漸蔓延(沈亞強(qiáng)等,2010;邱成和葉春麗,2015)。本研究調(diào)查發(fā)現(xiàn),杞麓湖水體SPC顯著大于其他4個湖泊(Plt;0.05)(圖2),耐鹽能力較強(qiáng)的穗花狐尾藻(劉萌萌等,2019)、篦齒眼子菜(王衛(wèi)紅和季民,2007)為主要優(yōu)勢物種。

    pH是沉水植物生存與繁衍的重要環(huán)境因子,其通過改變水體溶解無機(jī)碳的平衡關(guān)系,進(jìn)而影響水生植物的生長和分布(任南等,1996;何俊等,2008)。已有研究表明,穗花狐尾藻、篦齒眼子菜在pH值為7.5~8.8時對HCO3-形式的碳源都表現(xiàn)出較強(qiáng)的利用能力,從而能適應(yīng)高pH的水環(huán)境(蘇睿麗和李偉,2005;劉玲玲,2011)。本次調(diào)查的5個云南高原湖泊的pH值平均范圍為7.14~7.77,穗花狐尾藻和篦齒眼子菜為主要優(yōu)勢種和亞優(yōu)勢種,也證實(shí)了這一結(jié)論。以往研究還發(fā)現(xiàn),當(dāng)pH為7~9時,輪葉黑藻更易在富營養(yǎng)狀態(tài)水體中生長(張瀏等,2007),本次調(diào)查結(jié)果與其一致,即在富營養(yǎng)湖泊(異龍湖、星云湖)中輪葉黑藻更易成為優(yōu)勢類群。

    由于沉水植物整個生活史大部分位于水中,底質(zhì)和水體中氮元素可直接被其根系和葉片吸收,其中NH4+-N和NO3--N兩類基礎(chǔ)無機(jī)氮源最易被吸收利用(段娜等,2015)。本研究中湖泊NH4+-N和NO3-的平均濃度為0.28 mg/L和0.37 mg/L,其中富營養(yǎng)湖泊NH4+-N為0.40 mg/L、NO3--N為0.53 mg/L,貧-中營養(yǎng)湖泊NH4+-N為0.12 mg/L、NO3--N為0.17 mg/L。已有研究表明,在一定范圍內(nèi),較高的氮素濃度有利于馬來眼子菜的生長(蘇睿麗和李偉,2005)。本研究發(fā)現(xiàn)馬來眼子菜與NO3--N濃度呈顯著正相關(guān)(圖5c),其在富營養(yǎng)湖泊異龍湖中更容易成為優(yōu)勢物種。

    水下光照條件是沉水植物進(jìn)行正常光合作用的重要前提,是影響沉水植物在天然水體中生長的主要因素(王華等,2008)。已有研究表明,沉水植物可正常生長的水深通常為水體SD的2~2.5倍(邱成和葉春麗,2015)。本研究發(fā)現(xiàn),水體SD是影響云南高原湖泊沉水植物分布最關(guān)鍵的環(huán)境因子(圖5c);貧-中營養(yǎng)湖泊(撫仙湖和陽宗海)的SD要顯著大于富營養(yǎng)湖泊(Plt;0.05),其中撫仙湖的SD為(7.77±0.03)m(圖2)。由于沉水植物對光和空間的競爭能力不同,可粗分為草甸型植物和冠層型植物2類(James et al,2004)。其中,草甸型植物具有莖或分枝短、植冠低矮、生長較慢的生物學(xué)特征,絕大部分生活在沉積物附近,對光和生存空間具有相對較弱的競爭能力,如輪藻類、苦草、篦齒眼子菜和微齒眼子菜等物種;相反,冠層型植物具有較長的莖或分枝,其植冠高大、生長迅速,可生長到水面附近,對光和空間的具有較強(qiáng)的競爭能力,如穗花狐尾藻、黑藻和金魚藻等物種(Gopal amp; Usha,1993)。由于撫仙湖、陽宗海具有較高的SD,故苦草、篦齒眼子菜、穿葉眼子菜和微齒眼子菜等草甸型沉水植物可以在該類型湖泊中成為優(yōu)勢種。

    沉水植物的生長分布與水環(huán)境因子是相互影響的。草甸型植物比冠層型植物更有利于維持湖泊的清水穩(wěn)態(tài)(Carpenter amp; Lodge,1986;Rehman et al,2016)。草甸型沉水植物多具有相對稠密的須根,且絕大部分生活在沉積物附近,能有效阻止沉積物向水體釋放營養(yǎng)鹽,顯著提升沉積物滯留營養(yǎng)鹽的能力(Kufel amp; Kufel,2002;James et al,2004)。冠層型沉水植物生長較快,其衰敗腐爛時會釋放大量腐殖質(zhì)和營養(yǎng)鹽進(jìn)入水體,導(dǎo)致浮游植物快速繁殖,引起透明度下降(Engelhardt amp; Ritchie,2001;Hilt et al,2006)。因此,當(dāng)沉水植物類型的轉(zhuǎn)變,也預(yù)示著湖泊營養(yǎng)水平的轉(zhuǎn)變(da Silva et al,2014)。另外,水環(huán)境因子之間也是相互影響、相互制約的,單個環(huán)境因子可以通過影響其他環(huán)境因子進(jìn)而影響沉水植物的生長和分布。因此,在研究影響水生植物生長分布的環(huán)境因子時,應(yīng)當(dāng)綜合考慮各環(huán)境因子之間的內(nèi)在聯(lián)系,全面系統(tǒng)分析水生植物與環(huán)境因子的關(guān)系。

    綜合本次調(diào)查研究結(jié)果,提出以下建議:(1)針對類似本研究的貧-中營養(yǎng)類型湖泊,要遵循沉水植物的生物學(xué)特征和生長規(guī)律,及時收割清除高大的冠層型沉水植物,創(chuàng)造草甸型沉水植物的生長條件,促進(jìn)沉水植物群落向貧營養(yǎng)化湖泊的群落結(jié)構(gòu)方向發(fā)展;(2)針對類似本研究的富營養(yǎng)類型湖泊,其沉水植物恢復(fù)的關(guān)鍵要控制水體營養(yǎng)鹽濃度,抑制藻類的生長,采用不同的恢復(fù)方法,選取不同的植物種類來制定沉水植物的恢復(fù)方案。

    參考文獻(xiàn)

    段娜, 章堯想, 劉芳, 等, 2015. 植物氮素吸收及其轉(zhuǎn)運(yùn)蛋白研究進(jìn)展[J].分子植物育種, 13(2): 461-468.

    方馨, 趙風(fēng)斌, 柴曉利, 等, 2021. 異龍湖沉水植物分布格局與水環(huán)境因子相關(guān)性研究[J].長江流域資源與環(huán)境," 30(3): 636-643.

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

    高弋明, 殷春雨, 劉霞, 等, 2021. 撫仙湖近60年來沉水植物群落變化趨勢分析[J].湖泊科學(xué), 33(4): 1209-1219.

    何俊, 谷孝鴻, 劉國鋒, 2008. 東太湖水生植物及其與環(huán)境的相互作用[J]. 湖泊科學(xué), 20(6): 47-49.

    金相燦, 1995. 中國湖泊環(huán)境[M]. 北京: 海洋出版社.

    金相燦, 楚建周, 王圣瑞, 2007. 水體氮濃度、形態(tài)對黑藻和狐尾藻光合特征的影響[J].應(yīng)用與環(huán)境生物學(xué)報, 13(2): 200-204.

    孔祥虹, 肖蘭蘭, 蘇豪杰, 等, 2015. 長江下游湖泊水生植物現(xiàn)狀及與水環(huán)境因子的關(guān)系[J].湖泊科學(xué), 27(3): 385-391.

    藍(lán)紅林, 2001. 云南高原湖泊的治理與保護(hù)初探[J].云南環(huán)境科學(xué), 20(4): 26-28.

    李恒, 1980. 云南高原湖泊水生植被的研究[J].云南植物研究, 2(2): 113-141.

    李文朝, 連光華, 1996. 幾種沉水植物營養(yǎng)繁殖體萌發(fā)的光需求研究[J]. 湖泊科學(xué)," 8(增刊1): 25-29.

    劉玲玲, 2011. 三種沉水植物無機(jī)碳利用機(jī)制研究[D].武漢:華中師范大學(xué).

    劉萌萌, 劉巧," 楊娜," 等, 2019. 沉水植物穗花狐尾藻耐鹽性與生長[J].生態(tài)學(xué)雜志, 38(3): 778-784.

    秦伯強(qiáng), 2002. 長江中下游淺水湖泊富營養(yǎng)化發(fā)生機(jī)制與控制途徑初探[J].湖泊科學(xué), 14(3): 193-202.

    邱成, 葉春麗, 2015. 淺析江川縣發(fā)展生態(tài)畜牧業(yè)與星云湖水環(huán)境保護(hù)[C]. 中國環(huán)境科學(xué)學(xué)會環(huán)境經(jīng)濟(jì)學(xué)分會2015年學(xué)術(shù)年會.

    任南, 嚴(yán)國安, 馬劍敏, 等, 1996. 環(huán)境因子對東湖幾種沉水植物生理的影響研究[J].武漢大學(xué)學(xué)報(自然科學(xué)版), 42(2): 213-218.

    沈亞強(qiáng), 王海軍, 劉學(xué)勤, 2010. 滇中五湖水生植物區(qū)系及沉水植物群落特征[J]. 長江流域資源與環(huán)境, 19(增刊1): 111-119.

    蘇睿麗, 李偉, 2005. 沉水植物光合作用的特點(diǎn)與研究進(jìn)展[J].植物學(xué)通報, 22(增刊1): 128-138.

    王華, 逄勇, 劉申寶, 等, 2008. 沉水植物生長影響因子研究進(jìn)展[J].生態(tài)學(xué)報," 28(8): 3958-3968.

    王明翠, 劉雪芹, 張建輝, 2002. 湖泊富營養(yǎng)化評價方法及分級標(biāo)準(zhǔn)[J].中國環(huán)境監(jiān)測, 18(5): 47-49.

    王衛(wèi)紅, 季民, 2007. 9種沉水植物的耐鹽性比較[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 26(4): 1259-1263.

    王文林, 王國祥, 李強(qiáng), 等, 2006. 水體濁度對菹草(Potamogeton cripus)幼苗生長發(fā)育的影響[J].生態(tài)學(xué)報, 26(11): 3586-3593.

    魏翔, 唐光明, 2014. 異龍湖近20年來營養(yǎng)鹽與水生生態(tài)系統(tǒng)變化[J].環(huán)境科學(xué)導(dǎo)刊, 33(2): 9-14.

    吳愛平, 吳世凱, 倪樂意, 2005. 長江中游淺水湖泊水生植物氮磷含量與水柱營養(yǎng)的關(guān)系[J].水生生物學(xué)報," 29(4): 406-412.

    謝貽發(fā), 李傳紅, 劉正文, 等, 2007. 基質(zhì)條件對苦草(Vallisneria natans)生長和形態(tài)特征的影響[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 26(4): 1269-1272.

    楊加林, 李杰, 李經(jīng)緯, 等, 2012. 撫仙湖-星云湖水生生物與水環(huán)境研究綜述[J].云南地理環(huán)境研究, 24(2): 98-102.

    楊柳, 黃耀志, 2002. 云南高原湖區(qū)綜合發(fā)展研究[J].重慶大學(xué)學(xué)報(社會科學(xué)版)," 8(3): 19-21.

    顏昌宙, 金相燦, 趙景柱, 等, 2005. 云南洱海的生態(tài)保護(hù)及可持續(xù)利用對策[J]. 環(huán)境科學(xué), 26(5): 38-42.

    閆志強(qiáng), 劉黽, 吳小業(yè), 等, 2014. 溫度對五種沉水植物生長和營養(yǎng)去除效果的影響[J].生態(tài)科學(xué), 33(5): 839-844.

    游灝, 2006. 五種沉水植物對富營養(yǎng)化水體的生態(tài)適應(yīng)性研究[D].南京:南京農(nóng)業(yè)大學(xué).

    張瀏, 陳燦, 高倩, 等," 2007. 兩種營養(yǎng)狀態(tài)下pH對輪葉黑藻(Hydrilla verticillata)生長和抗氧化酶活性的影響[J].生態(tài)環(huán)境, 16(3): 748-752.

    曾理, 吳豐昌, 萬國江, 等, 2009. 中國地區(qū)湖泊沉積物中137Cs分布特征和環(huán)境意義[J].湖泊科學(xué), 21(1): 1-9.

    Al-Kenzawi M," Al-Haidary M," Talib A H," et al," 2011. Environmental study of some water characteristics at Um-Al-Naaj Marsh," south of Iraq[J]. Baghdad Science Journal," 8(1): 531-538.

    Al-Thahaibawi B," Al-Abbawy D," Al-Mayaly I," et al," 2021. The effect of some environmental variables on the abundance and distribution of aquatic macrophytes in Al-Hawizeh marsh," Iraq[J]. IOP Conference Series: Earth and Environmental Science," 722(1): 012001.

    Al-Shawi I," 2006. Comparative study of some physico-chemical characteristics for Northern Al-Hammar marsh waters before destroyed and after Rehabilitation 2004[J]. Marsh Bulletin," 2:127-133.

    Barko J W," Smart R M," 1981. Comparative influences of light and temperature on the growth and metabolism of selected submersed freshwater macrophytes[J]. Ecological Monographs," 51(2): 219-235.

    Blanchet F G," Legendre P," Borcard D," 2008. Forward selection of explanatory variables[J]. Ecology," 89(9): 2623-2632.

    Bornette G," Puijalon S," 2011. Response of aquatic plants to abiotic factors: a review[J]. Aquatic Sciences," 73(1): 1-14.

    Carpenter S R," Lodge D M," 1986. Effects of submersed macrophytes on ecosystem processes[J]. Aquatic Botany," 26: 341-370.

    da Silva S C A D," Cervi A C," Bona C," et al," 2014. Aquatic macrophyte community varies in urban reservoirs with different degrees of eutrophication[J]. Acta Limnol Bras," 26(2): 129-142.

    Engelhardt K," Ritchie M E," 2001. Effects of macrophyte species richness on wetland ecosystem functioning and services[J]. Nature," 411(6838): 687-689.

    EeO\"s T," S[a][']ly P," Tak[a][']cs P," et al," 2012. Temporal variability in the spatial and environmental determinants of functional metacommunity organization-stream fish in a human-modified landscape[J]. Freshwater Biology," 57(9): 1914-1928.

    Gopal B," Usha G," 1993. Competition and allelopathy in aquatic plant communities[J]. The Botanical Review," 59(3): 155-210.

    Heino J," 2002. Concordance of species richness patterns among multiple freshwater taxa: a regional perspective[J]. Biodiversity amp; Conservation," 11(1):137-147.

    Hilt S," Gross E," Hupfer M," et al," 2006. Restoration of submerged vegetation in shallow eutrophic lakes: a guideline and state of the art in Germany[J]. Limnologica- Ecology and Management of Inland Waters," 36(3): 155-171.

    Horppila J," Nurminen L," 2005. Effects of different macrophyte growth forms on sediment and P resuspension in a shallow lake[J]. Hydrobiologia," 545(1): 167-175.

    James W F," Barko J W," Butler M G," 2004. Shear stress and sediment resuspension in relation to submersed macrophyte biomass[J]. Hydrobiologia," 515(1/2/3): 181-191.

    Jeppesen E," Jensen J P," M S?ndergaard," et al," 2000. Trophic structure," species richness and biodiversity in Danish lakes: changes along a phosphorus gradient[J]. Freshwater Biology," 45(2): 201-218.

    Kalff Jacob," 2002. Limnology: Inland Water Ecosystems[M].Prentice Hall," NJ: Benjamin Cummings.

    Kufel L," Kufel I," 2002. Chara beds acting as nutrient sinks in shallow lakes: a review[J]. Aquatic Botany," 72(3/4): 249-260.

    Wang H," Wang Q," Bowler P A," et al," 2016. Invasive aquatic plants in China[J]. Aquatic Invasions," 11(1): 1-9.

    Nakamura K," Kayaba Y," Nishihiro J," et al," 2008. Effects of submerged plants on water quality and biota in large-scale experimental ponds[J]. Landscape amp; Ecological Engineering," 4(1): 1-9.

    Rehman F," Pervez A," Khattak B N," et al," 2016. Constructed wetlands: perspectives of the oxygen released in the rhizosphere of macrophytes[J]. Clean-Soil," Air," Water," 45(1).DOI:10.1002/clen.201600054.

    Rooney N," Kalff J," 2000. Inter-annual variation in submerged macrophyte community biomass and distribution: the influence of temperature and lake morphometry[J]. Aquatic Botany," 68(4): 321-335.

    Sachse R," Petzoldt T," Blumstock M," et al," 2014. Extending one-dimensional models for deep lakes to simulate the impact of submerged macrophytes on water quality[J]. Environmental Modelling amp; Software," 61: 410-423.

    Spencer D F," Ksander G G," Madsen J D," et al," 2000. Emergence of vegetative propagules of Potamogeton nodosus," Potamogeton pectinatus," Vallisneria americana," and Hydrilla verticillata based on accumulated degree-days[J]. Aquatic Botany," 67(3): 237-249.

    Welch P S," 1952. Limnology[M]. New York: McGraw Hill Book Company.

    (責(zé)任編輯" "鄭金秀)

    Community Composition of Submerged Macrophytes and Relationship

    with Environmental Factors in Five Yunnan Plateau Lakes

    ZHU Li‐ming, ZHANG Dao‐xi, SHI Fang, ZOU Xi, ZHANG Zhi‐yong, HU Lian

    (Key Laboratory of Ecological Impacts of Hydraulic-Projects and Restoration of Aquatic Ecosystem

    of Ministry of Water Resources, Institute of Hydroecology, Ministry of Water Resources amp; Chinese

    Academy of Sciences, Wuhan" "430079, P.R. China)

    Abstract:Many studies have suggested that restoration of the submerged plant community is an important means for controlling lake eutrophication. Therefore, it is important to study the composition and distribution of the submerged plant communities in lakes and the relationship with the water environment factors. In this study, five typical plateau lakes in central Yunnan Province were selected for research, including Yilong Lake (YLH), Qilu Lake (QLH), Xingyun Lake (XYH), Fuxian Lake (FXH) and Yangzonghai lake (YZH). We investigated the composition and distribution of submerged macrophytes communities in the five lakes, and the relationship between the distribution of submerged macrophytes and water environmental factors. A two-year investigation of submerged plants and water environment factors were conducted in the five lakes during September and October of 2020 and 2021. The comprehensive nutrient status indices of Yilong Lake, Qilu Lake and Xingyun Lake were 63.55±6.28, 67.26±1.89 and 72.08±2.04, respectively, indicating an eutrophic state, and the comprehensive nutrient status indices of Fuxian Lake (FXH) and Yangzonghai Lake (YZH) were 28.01±4.82 and 41.47±1.96, respectively, indicating an oligotrophic-mesotrophic state. A total of 15 submerged macrophyte species belonging to 6 families and 10 genera were identified during the investigation. Myriophyllum spicatum, P. pectinatus were dominant species in all lakes, Potamogeton maackianus and P. perfoliatus were dominant in the oligotrophic-mesotrophic lakes, and Hydrilla verticillata, P. crispus and Ceratophyllum demersum were dominant in the eutrophic lakes. The Shannon-Weiner diversity (H') index of the submerged macrophytes in FXH was highest, with an average value of 1.63, while the diversity index of QLH was lowest, with an average value of only 0.79. In general, submerged macrophyte diversity was higher in the oligotrophic-mesotrophic lakes than in the eutrophic lakes. Redundancy analysis of the submerged plants and the environmental factors shows that transparency, pH and NO3--N were the key factors affecting the distribution of submerged macrophytes in lakes with different nutrient levels; pH, NH4+-N and WT were the key factors in oligotrophic-mesotrophic lakes; and WT, conductivity and Chl-a were the key factors in eutrophic lakes. In conclusion, we identified the key water environmental factors affecting the distribution of submerged macrophytes in lakes on the Yunnan Plateau, which will provide guidance for the restoration, protection and management of submerged macrophytes in lakes in Yunnan plateau lakes.

    Key words: plateau lakes; submerged macrophytes; spatial distribution; environmental factors; redundancy analysis

    猜你喜歡
    環(huán)境因子空間分布沉水植物
    沉水植物生長影響因子研究進(jìn)展
    花卉(2020年12期)2020-01-09 07:12:52
    沉水植物種植工程實(shí)施后物種多樣性的變化
    ——以武漢紫陽湖為例
    綠色科技(2019年2期)2019-05-21 07:43:18
    白龜山濕地重金屬元素分布特征及其來源分析
    綠色科技(2016年20期)2016-12-27 17:34:13
    阿克蘇地區(qū)杏李與歐洲李氣體交換特征對比研究
    基于GIS技術(shù)的福建省柳葉白前資源適宜性空間分布研究
    江蘇省臭氧污染變化特征
    科技視界(2016年18期)2016-11-03 23:51:58
    鐵路客流時空分布研究綜述
    桉樹人工林胸徑變化特征與環(huán)境因子的關(guān)系研究
    中國市場(2016年36期)2016-10-19 05:25:57
    沉水植物在修復(fù)富營養(yǎng)湖體的研究進(jìn)展
    大科技(2016年1期)2016-07-17 02:42:07
    秦皇島海域夜光藻種群密度與環(huán)境因子的關(guān)系
    国产乱人偷精品视频| 国产精品免费大片| 亚洲国产欧美日韩在线播放| 国产精品国产三级专区第一集| 久热爱精品视频在线9| 日韩人妻精品一区2区三区| 国产精品99久久99久久久不卡 | 亚洲国产精品一区三区| 午夜免费男女啪啪视频观看| 黑人猛操日本美女一级片| 久久国产精品男人的天堂亚洲| 欧美日韩亚洲国产一区二区在线观看 | 一级片免费观看大全| 国产精品99久久99久久久不卡 | 亚洲国产欧美网| 亚洲成色77777| 日韩欧美精品免费久久| 一级毛片黄色毛片免费观看视频| av电影中文网址| 免费女性裸体啪啪无遮挡网站| 老司机亚洲免费影院| 欧美精品一区二区大全| 不卡av一区二区三区| 亚洲一码二码三码区别大吗| 女人爽到高潮嗷嗷叫在线视频| 天天添夜夜摸| 老司机深夜福利视频在线观看 | 久久久国产一区二区| 两个人免费观看高清视频| 国产亚洲最大av| 亚洲,欧美,日韩| 久久 成人 亚洲| 国产精品av久久久久免费| 我要看黄色一级片免费的| 最近最新中文字幕大全免费视频 | 精品一区二区三区av网在线观看 | 91老司机精品| 看非洲黑人一级黄片| 1024视频免费在线观看| 亚洲久久久国产精品| 涩涩av久久男人的天堂| 久久人人爽人人片av| 好男人视频免费观看在线| 欧美精品一区二区大全| 国产精品久久久久久久久免| 色吧在线观看| 国产午夜精品一二区理论片| 最近中文字幕高清免费大全6| 国产男女超爽视频在线观看| av女优亚洲男人天堂| 黄片无遮挡物在线观看| 日日啪夜夜爽| 国产一级毛片在线| av电影中文网址| 亚洲情色 制服丝袜| 十分钟在线观看高清视频www| 国产极品天堂在线| 大陆偷拍与自拍| 妹子高潮喷水视频| 精品久久蜜臀av无| 亚洲国产日韩一区二区| 在线 av 中文字幕| 国产 一区精品| 国产成人系列免费观看| 中国国产av一级| 黑丝袜美女国产一区| 久久久久国产精品人妻一区二区| 日本午夜av视频| 美女午夜性视频免费| 精品少妇久久久久久888优播| 美女脱内裤让男人舔精品视频| 亚洲五月色婷婷综合| 日韩一卡2卡3卡4卡2021年| 亚洲精品一二三| 精品国产国语对白av| 在线观看免费日韩欧美大片| 91成人精品电影| 九九爱精品视频在线观看| 久久久久精品久久久久真实原创| 一级爰片在线观看| 久久99热这里只频精品6学生| 欧美xxⅹ黑人| 悠悠久久av| 国产精品三级大全| 国产一区二区激情短视频 | 午夜福利乱码中文字幕| 91精品国产国语对白视频| a级毛片在线看网站| 人人妻,人人澡人人爽秒播 | 久久久久精品人妻al黑| 中文乱码字字幕精品一区二区三区| 国产精品一国产av| 搡老岳熟女国产| 黄色毛片三级朝国网站| 麻豆av在线久日| 高清在线视频一区二区三区| 中文字幕最新亚洲高清| 一级毛片黄色毛片免费观看视频| 啦啦啦啦在线视频资源| av又黄又爽大尺度在线免费看| 欧美黑人欧美精品刺激| 老司机影院成人| 少妇人妻 视频| 国产成人一区二区在线| 色婷婷久久久亚洲欧美| 99热国产这里只有精品6| 热99久久久久精品小说推荐| 国语对白做爰xxxⅹ性视频网站| 亚洲美女黄色视频免费看| avwww免费| 亚洲欧美清纯卡通| 午夜老司机福利片| 国产成人一区二区在线| 国产精品一区二区精品视频观看| 性少妇av在线| 不卡视频在线观看欧美| 一级a爱视频在线免费观看| 中国三级夫妇交换| 女性生殖器流出的白浆| 麻豆精品久久久久久蜜桃| 日韩av免费高清视频| 大陆偷拍与自拍| 伊人久久大香线蕉亚洲五| 国产高清国产精品国产三级| 成人毛片60女人毛片免费| 丁香六月天网| 午夜福利视频精品| 十八禁人妻一区二区| 免费黄频网站在线观看国产| 日韩大片免费观看网站| 久久人人97超碰香蕉20202| 久久久久久久久免费视频了| 日韩av不卡免费在线播放| 国产日韩欧美在线精品| 亚洲欧美一区二区三区久久| 69精品国产乱码久久久| 男女床上黄色一级片免费看| 亚洲欧美一区二区三区国产| av网站在线播放免费| 国产熟女欧美一区二区| 久久久久久人妻| 日日爽夜夜爽网站| 哪个播放器可以免费观看大片| 久久av网站| 国产日韩欧美视频二区| √禁漫天堂资源中文www| 久久精品国产亚洲av高清一级| 精品人妻一区二区三区麻豆| 成人漫画全彩无遮挡| 男女免费视频国产| 午夜影院在线不卡| 超碰成人久久| 18禁裸乳无遮挡动漫免费视频| 美女脱内裤让男人舔精品视频| 国产乱人偷精品视频| kizo精华| 国产又色又爽无遮挡免| av视频免费观看在线观看| 成人午夜精彩视频在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲第一av免费看| 欧美日韩av久久| 国产片内射在线| 两性夫妻黄色片| 成年人免费黄色播放视频| 国产欧美日韩一区二区三区在线| 国精品久久久久久国模美| 日韩不卡一区二区三区视频在线| 波野结衣二区三区在线| 国产成人系列免费观看| 中文乱码字字幕精品一区二区三区| 久久ye,这里只有精品| www.熟女人妻精品国产| 国产人伦9x9x在线观看| 国产亚洲最大av| 久久影院123| 久久久久久久久免费视频了| 看免费av毛片| 爱豆传媒免费全集在线观看| 亚洲激情五月婷婷啪啪| 男女床上黄色一级片免费看| 亚洲av成人不卡在线观看播放网 | 97在线人人人人妻| 免费日韩欧美在线观看| 欧美中文综合在线视频| 又大又爽又粗| 亚洲中文av在线| 十八禁高潮呻吟视频| 国产精品嫩草影院av在线观看| 男女高潮啪啪啪动态图| 亚洲美女黄色视频免费看| 成人免费观看视频高清| 日本91视频免费播放| 一级a爱视频在线免费观看| 色视频在线一区二区三区| 久久久精品国产亚洲av高清涩受| 人体艺术视频欧美日本| 午夜福利乱码中文字幕| 免费看av在线观看网站| 男男h啪啪无遮挡| 少妇的丰满在线观看| 一级,二级,三级黄色视频| 国产在线视频一区二区| av在线app专区| 777久久人妻少妇嫩草av网站| 亚洲情色 制服丝袜| 久久久精品国产亚洲av高清涩受| 国产免费福利视频在线观看| 亚洲欧美成人精品一区二区| 亚洲av日韩在线播放| 日本色播在线视频| av在线播放精品| 中文字幕亚洲精品专区| 免费少妇av软件| 亚洲三区欧美一区| av卡一久久| 熟女少妇亚洲综合色aaa.| 香蕉国产在线看| 99国产综合亚洲精品| 另类亚洲欧美激情| 丝瓜视频免费看黄片| 老司机亚洲免费影院| 亚洲情色 制服丝袜| av片东京热男人的天堂| 麻豆av在线久日| 成人亚洲精品一区在线观看| 美女午夜性视频免费| 飞空精品影院首页| 五月天丁香电影| av国产精品久久久久影院| 成人三级做爰电影| 制服诱惑二区| 久久99一区二区三区| 99国产精品免费福利视频| 久久久国产欧美日韩av| 亚洲三区欧美一区| 搡老乐熟女国产| 亚洲伊人久久精品综合| 日韩av不卡免费在线播放| 波野结衣二区三区在线| 免费看不卡的av| 精品亚洲成a人片在线观看| 在线免费观看不下载黄p国产| 好男人视频免费观看在线| 亚洲欧洲精品一区二区精品久久久 | 看十八女毛片水多多多| 久久婷婷青草| 嫩草影视91久久| 亚洲熟女毛片儿| av不卡在线播放| 啦啦啦中文免费视频观看日本| 啦啦啦在线观看免费高清www| 99精品久久久久人妻精品| 国产女主播在线喷水免费视频网站| 看免费成人av毛片| 赤兔流量卡办理| 国产日韩欧美亚洲二区| 另类精品久久| 亚洲国产精品成人久久小说| 嫩草影院入口| 精品一区二区三区av网在线观看 | 日本欧美视频一区| 自线自在国产av| 啦啦啦啦在线视频资源| 一边摸一边抽搐一进一出视频| 亚洲精品av麻豆狂野| 亚洲欧洲日产国产| 中国三级夫妇交换| 日本一区二区免费在线视频| 国产黄色免费在线视频| 亚洲七黄色美女视频| 国产免费现黄频在线看| 亚洲,欧美,日韩| 成人午夜精彩视频在线观看| 波多野结衣av一区二区av| 午夜精品国产一区二区电影| 亚洲五月色婷婷综合| www日本在线高清视频| 亚洲婷婷狠狠爱综合网| 精品少妇久久久久久888优播| 丰满乱子伦码专区| 精品久久蜜臀av无| 性少妇av在线| 人人妻人人澡人人看| 波多野结衣一区麻豆| 国产精品国产三级专区第一集| 大话2 男鬼变身卡| 99国产综合亚洲精品| 性高湖久久久久久久久免费观看| 亚洲,欧美精品.| 国产精品国产av在线观看| 国产精品香港三级国产av潘金莲 | 亚洲精品第二区| 国产老妇伦熟女老妇高清| 日韩欧美精品免费久久| 精品少妇内射三级| 日日摸夜夜添夜夜爱| 日本爱情动作片www.在线观看| 久久人人爽人人片av| 国产日韩欧美在线精品| 国产一区二区激情短视频 | 777久久人妻少妇嫩草av网站| 一个人免费看片子| 国产xxxxx性猛交| 免费少妇av软件| 日韩视频在线欧美| 精品一区二区三卡| 亚洲国产看品久久| 成年动漫av网址| 久热这里只有精品99| 日韩精品免费视频一区二区三区| 男女无遮挡免费网站观看| 欧美97在线视频| 久久精品熟女亚洲av麻豆精品| 男女无遮挡免费网站观看| 日韩电影二区| 日韩熟女老妇一区二区性免费视频| 国产精品亚洲av一区麻豆 | 人人妻人人澡人人爽人人夜夜| 久久久久国产一级毛片高清牌| 国产精品久久久久久精品古装| 看免费成人av毛片| 国产激情久久老熟女| 精品少妇一区二区三区视频日本电影 | 亚洲国产欧美日韩在线播放| 久久久久久久久久久免费av| 色综合欧美亚洲国产小说| 青春草国产在线视频| 精品午夜福利在线看| 亚洲熟女毛片儿| 精品少妇黑人巨大在线播放| 18禁裸乳无遮挡动漫免费视频| 最近2019中文字幕mv第一页| 成年人免费黄色播放视频| 一本色道久久久久久精品综合| 激情视频va一区二区三区| 精品酒店卫生间| 免费久久久久久久精品成人欧美视频| 久久久久久免费高清国产稀缺| 国产高清不卡午夜福利| 高清av免费在线| 老汉色∧v一级毛片| 欧美日韩亚洲综合一区二区三区_| 欧美精品av麻豆av| 亚洲av综合色区一区| 欧美激情高清一区二区三区 | 亚洲少妇的诱惑av| 亚洲av综合色区一区| 久久国产精品大桥未久av| 欧美国产精品一级二级三级| 蜜桃国产av成人99| 午夜福利免费观看在线| 国产一区二区在线观看av| 亚洲国产精品成人久久小说| 成人午夜精彩视频在线观看| 国产极品天堂在线| netflix在线观看网站| 久久国产精品男人的天堂亚洲| 制服丝袜香蕉在线| 久久久久人妻精品一区果冻| 啦啦啦啦在线视频资源| 性高湖久久久久久久久免费观看| 久久久久国产一级毛片高清牌| 国产乱来视频区| 国产成人精品福利久久| 交换朋友夫妻互换小说| 亚洲综合色网址| 99久国产av精品国产电影| 伊人久久国产一区二区| 国产精品久久久久久人妻精品电影 | 少妇的丰满在线观看| 亚洲成人一二三区av| 精品免费久久久久久久清纯 | 日韩大片免费观看网站| 在线观看www视频免费| 成人亚洲精品一区在线观看| 男女边吃奶边做爰视频| 人人妻,人人澡人人爽秒播 | 狠狠婷婷综合久久久久久88av| 我的亚洲天堂| 久久久精品免费免费高清| 国产男女内射视频| 国产又色又爽无遮挡免| 综合色丁香网| videosex国产| 我要看黄色一级片免费的| 肉色欧美久久久久久久蜜桃| 老司机在亚洲福利影院| 久久99精品国语久久久| 欧美日韩亚洲高清精品| 国产成人91sexporn| 男女午夜视频在线观看| av网站在线播放免费| 精品人妻在线不人妻| 多毛熟女@视频| 蜜桃在线观看..| 搡老乐熟女国产| 国产成人一区二区在线| 国产1区2区3区精品| 国产精品.久久久| 国产免费一区二区三区四区乱码| 免费观看人在逋| 久久精品国产综合久久久| 亚洲成人av在线免费| 我要看黄色一级片免费的| 国产成人精品久久二区二区91 | 中文字幕人妻丝袜制服| 久久久久网色| 欧美亚洲 丝袜 人妻 在线| 国产无遮挡羞羞视频在线观看| 69精品国产乱码久久久| 亚洲少妇的诱惑av| 可以免费在线观看a视频的电影网站 | 国产一区二区在线观看av| 91成人精品电影| 自拍欧美九色日韩亚洲蝌蚪91| 男女无遮挡免费网站观看| 国产亚洲最大av| 国产精品二区激情视频| 午夜日韩欧美国产| 51午夜福利影视在线观看| 丝袜脚勾引网站| 日韩大码丰满熟妇| 老汉色∧v一级毛片| 午夜老司机福利片| 日本欧美视频一区| 制服人妻中文乱码| 国产精品.久久久| 国语对白做爰xxxⅹ性视频网站| 成年美女黄网站色视频大全免费| 少妇 在线观看| 天天躁夜夜躁狠狠久久av| 免费在线观看视频国产中文字幕亚洲 | 午夜福利网站1000一区二区三区| 国产精品欧美亚洲77777| 久久精品国产a三级三级三级| 男女下面插进去视频免费观看| netflix在线观看网站| 高清视频免费观看一区二区| 国产精品免费大片| 国产野战对白在线观看| 制服诱惑二区| 黄色视频在线播放观看不卡| 欧美最新免费一区二区三区| 亚洲,一卡二卡三卡| 国产不卡av网站在线观看| 亚洲国产欧美在线一区| 天堂俺去俺来也www色官网| 99精国产麻豆久久婷婷| 91老司机精品| 赤兔流量卡办理| 综合色丁香网| 精品国产露脸久久av麻豆| av女优亚洲男人天堂| 亚洲精品,欧美精品| 男人舔女人的私密视频| 99热网站在线观看| 成人午夜精彩视频在线观看| 亚洲av电影在线观看一区二区三区| 黄色视频不卡| 热99久久久久精品小说推荐| xxxhd国产人妻xxx| 在线天堂最新版资源| 日韩 亚洲 欧美在线| 亚洲色图 男人天堂 中文字幕| 亚洲av在线观看美女高潮| 在线观看三级黄色| 黄片小视频在线播放| 午夜福利视频精品| 亚洲精品aⅴ在线观看| 黄网站色视频无遮挡免费观看| 99热国产这里只有精品6| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久精品国产欧美久久久 | 9热在线视频观看99| 中文欧美无线码| 亚洲国产精品国产精品| 国产成人啪精品午夜网站| 欧美人与性动交α欧美精品济南到| 三上悠亚av全集在线观看| 最近最新中文字幕免费大全7| 在线观看免费视频网站a站| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产欧美一区二区综合| 悠悠久久av| 亚洲成国产人片在线观看| 久久久精品免费免费高清| 人人澡人人妻人| 夫妻性生交免费视频一级片| 国产黄色视频一区二区在线观看| 免费观看a级毛片全部| av在线app专区| 男女国产视频网站| 国产免费现黄频在线看| 777久久人妻少妇嫩草av网站| 国产熟女午夜一区二区三区| 亚洲成人免费av在线播放| 成人影院久久| 黑丝袜美女国产一区| 狂野欧美激情性bbbbbb| 少妇人妻 视频| 国产精品香港三级国产av潘金莲 | 纯流量卡能插随身wifi吗| 可以免费在线观看a视频的电影网站 | 精品午夜福利在线看| 欧美黄色片欧美黄色片| 一区二区三区乱码不卡18| 国产黄色免费在线视频| 伊人久久国产一区二区| 国产片内射在线| 90打野战视频偷拍视频| 一级a爱视频在线免费观看| 最黄视频免费看| 色网站视频免费| 我要看黄色一级片免费的| 久久青草综合色| 国产精品二区激情视频| 黄片播放在线免费| 亚洲中文av在线| 婷婷色综合www| 男的添女的下面高潮视频| 国产精品国产三级国产专区5o| 波多野结衣一区麻豆| 老熟女久久久| 少妇的丰满在线观看| 午夜福利一区二区在线看| 国产视频首页在线观看| 80岁老熟妇乱子伦牲交| 成人漫画全彩无遮挡| 女人久久www免费人成看片| 亚洲婷婷狠狠爱综合网| 欧美日韩综合久久久久久| 女性被躁到高潮视频| 国产精品国产三级专区第一集| 女人被躁到高潮嗷嗷叫费观| 日本wwww免费看| 美女中出高潮动态图| 制服人妻中文乱码| 亚洲中文av在线| 国产精品久久久久成人av| 国产成人精品无人区| 亚洲熟女毛片儿| 久久精品亚洲av国产电影网| 欧美 日韩 精品 国产| 亚洲一级一片aⅴ在线观看| 欧美变态另类bdsm刘玥| 观看美女的网站| 老司机影院毛片| 亚洲一区中文字幕在线| 在线观看免费日韩欧美大片| 亚洲av成人不卡在线观看播放网 | 久久精品国产a三级三级三级| 狂野欧美激情性xxxx| 亚洲精品久久成人aⅴ小说| 黑人欧美特级aaaaaa片| 成年人免费黄色播放视频| 精品久久久精品久久久| 男女之事视频高清在线观看 | 天天躁日日躁夜夜躁夜夜| 国产女主播在线喷水免费视频网站| 天堂俺去俺来也www色官网| 日本一区二区免费在线视频| 精品人妻一区二区三区麻豆| 国产男女内射视频| 蜜桃在线观看..| 国产一区二区 视频在线| 精品少妇黑人巨大在线播放| 国产成人免费无遮挡视频| 欧美精品人与动牲交sv欧美| 亚洲图色成人| 青春草亚洲视频在线观看| 成年美女黄网站色视频大全免费| 欧美 日韩 精品 国产| 国产精品 欧美亚洲| 国产av精品麻豆| 国产精品久久久久久精品古装| 久久久国产一区二区| 久久久久精品人妻al黑| 一区二区三区乱码不卡18| 亚洲激情五月婷婷啪啪| 婷婷成人精品国产| 少妇人妻 视频| 美女午夜性视频免费| 中文字幕最新亚洲高清| 国产成人精品福利久久| 我要看黄色一级片免费的| 日韩人妻精品一区2区三区| 亚洲国产精品国产精品| 色视频在线一区二区三区| av网站免费在线观看视频| 纯流量卡能插随身wifi吗| 亚洲欧洲精品一区二区精品久久久 | 观看av在线不卡| 亚洲精品久久成人aⅴ小说| 久久精品国产综合久久久| 国产一区二区三区av在线| 成年动漫av网址| 亚洲欧美一区二区三区久久| 亚洲熟女精品中文字幕| 日韩欧美一区视频在线观看| 久久精品国产综合久久久| 国产 一区精品| 考比视频在线观看| 国产精品久久久久成人av| 久久久久精品久久久久真实原创| 亚洲欧美成人精品一区二区| 免费在线观看视频国产中文字幕亚洲 | 国产欧美亚洲国产| 侵犯人妻中文字幕一二三四区| 最新在线观看一区二区三区 | 亚洲七黄色美女视频| 国产深夜福利视频在线观看| 精品一区二区三区av网在线观看 | 99久国产av精品国产电影| 亚洲欧美成人综合另类久久久|