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

    西藏巴松錯(cuò)夏季浮游植物群落的分布

    2023-04-29 20:08:07劉洋安瑞志王陳潘成梅巴桑
    水生態(tài)學(xué)雜志 2023年6期
    關(guān)鍵詞:浮游植物群落結(jié)構(gòu)

    劉洋 安瑞志 王陳 潘成梅 巴桑

    摘要:揭示西藏巴松錯(cuò)浮游植物群落結(jié)構(gòu)的水平和垂直分布特征及其與環(huán)境間的關(guān)系,為進(jìn)一步研究高寒湖泊生態(tài)系統(tǒng)提供基礎(chǔ)資料。于2020年7月對(duì)巴松錯(cuò)出入河流及其湖體設(shè)置 15 個(gè)水平采樣點(diǎn),并選擇湖中心樣點(diǎn)B15 進(jìn)行垂直分層采樣,同時(shí)測(cè)定水體相關(guān)理化指標(biāo)。調(diào)查共采集水樣275個(gè),分析了浮游植物物種組成、細(xì)胞豐度、生物量和物種多樣性等群落結(jié)構(gòu)特征,并探討浮游植物對(duì)水平和垂直環(huán)境梯度變化的響應(yīng)。結(jié)果表明:(1)巴松錯(cuò)夏季共鑒定浮游植物167種,隸屬于5門(mén)8綱20目29科66屬,物種組成為硅藻-綠藻-藍(lán)藻型;(2)在水平方向上,平均細(xì)胞豐度為1.845×104 個(gè)/L,平均生物量為0.012 mg/L;在垂直方向上,平均細(xì)胞豐度為1.006×104 個(gè)/L,平均生物量為0.006 mg/L;(3)巴松錯(cuò)夏季浮游植物群落結(jié)構(gòu)水平分布較為均勻,群落復(fù)雜程度相對(duì)較高,較垂直分布更穩(wěn)定;(4)RDA分析結(jié)果表明,電導(dǎo)率是影響巴松錯(cuò)夏季浮游植物水平群落變異的主要環(huán)境因子,濁度是影響垂直變異的主要環(huán)境因子。

    關(guān)鍵詞:浮游植物;群落結(jié)構(gòu);水平分布;垂直分布;巴松錯(cuò)

    中圖分類(lèi)號(hào):Q179.1? ? ? ? 文獻(xiàn)標(biāo)志碼:A? ? ? ? 文章編號(hào):1674-3075(2023)06-0104-09

    作為湖泊生態(tài)系統(tǒng)的初級(jí)生產(chǎn)者,浮游植物的動(dòng)態(tài)分布受多因素控制,與水體透明度、濁度、水溫、鹽度、DO和pH等水環(huán)境因子密切相關(guān),同時(shí)受到水平和垂直環(huán)境梯度的影響(Ptacnik et al,2008)。孫萍等(2020)調(diào)查了熱帶東印度洋春季浮游植物群落結(jié)構(gòu)水平和垂向上浮游植物種類(lèi)及豐度的差異性,邢冰偉等(2020)分析了九寨溝國(guó)家級(jí)自然保護(hù)區(qū)長(zhǎng)海夏季浮游植物群落結(jié)構(gòu)的水平和垂直分布,然而有關(guān)青藏高原湖泊浮游植物群落結(jié)構(gòu)水平和垂直分布方面的相關(guān)研究鮮有報(bào)道。

    青藏高原高寒湖泊分布星羅密布,占中國(guó)湖泊總面積一半以上,是高原陸-氣系統(tǒng)的基本組成要素(靳錚等,2020)。相對(duì)于其他湖泊而言,高寒湖泊受人類(lèi)活動(dòng)干擾較小,對(duì)氣候變化的響應(yīng)更快、更敏感,基本保持原始生態(tài)系統(tǒng)。自然狀態(tài)下湖泊水生生物群落組成及其對(duì)環(huán)境的響應(yīng),已引起國(guó)內(nèi)外學(xué)者的廣泛關(guān)注(王捷等,2015;陳立婧等,2013;Li et al,2021)。巴松錯(cuò)位于尼洋河最大支流巴河的寬谷中,是西藏東南部最大的高山冰川堰塞湖之一,在維持區(qū)域生態(tài)平衡、調(diào)節(jié)氣候、保持生物多樣性和提高居民生產(chǎn)、生活等方面起著重要作用。目前,有關(guān)巴松錯(cuò)的水生生物學(xué)研究主要集中在浮游植物(安瑞志等,2020;2021)和原生動(dòng)物(楊欣蘭等,2022)等方面,但針對(duì)巴松錯(cuò)浮游植物的系統(tǒng)和持續(xù)研究仍較少。

    為了解西藏巴松錯(cuò)浮游植物群落結(jié)構(gòu)水平和垂直分布差異及其與環(huán)境的關(guān)系,于2020年7月對(duì)巴松錯(cuò)出入河流及其湖體的浮游植物進(jìn)行了水平和垂直分層定性定量采集,研究分析了浮游植物物種組成、細(xì)胞豐度、生物量和物種多樣性等群落結(jié)構(gòu)特征,掌握巴松錯(cuò)水平和垂直上浮游植物群落分布格局和規(guī)律,并探討浮游植物對(duì)水平和垂直環(huán)境梯度變化的響應(yīng),以期為進(jìn)一步研究高寒湖泊生態(tài)系統(tǒng)提供資料。

    1? ?材料與方法

    1.1? ?研究區(qū)域概況

    巴松錯(cuò)(29°59′~30°04′N(xiāo),93°53~94°01′E)位于念青唐古拉山東端南麓、西藏東南部和林芝市西北部工布江達(dá)境內(nèi),屬尼洋河水系,外觀狹長(zhǎng)呈新月形,分為東、西2個(gè)湖盆,湖面平均海拔3 469 m,最大湖深約120 m,湖泊面積25.9 km2,流域面積 1 290 km2。該湖泊四面環(huán)山,氣候溫和,屬于典型的高原溫帶半濕潤(rùn)季風(fēng)氣候,年均降水量 646 mm,主要集中在4-10月,年平均氣溫6.3 ℃,年無(wú)霜期175 d,年日照時(shí)數(shù)為2 016 h(中國(guó)科學(xué)院青藏高原綜合科學(xué)考察隊(duì),1981)。湖水補(bǔ)給主要為天然降水和仲措弄巴曲、新措弄巴曲、白朗曲、羅結(jié)曲和扎拉河的冰川融雪水,經(jīng)由巴河流出,最終匯入尼洋河(中國(guó)科學(xué)院青藏高原綜合科學(xué)考察隊(duì),1984)。

    1.2? ?樣點(diǎn)設(shè)置

    2020年7月,根據(jù)巴松錯(cuò)的生態(tài)環(huán)境特征、功能分區(qū)與梯度結(jié)構(gòu)等特征,在巴松錯(cuò)出入河流及其湖體設(shè)置15個(gè)水平采樣點(diǎn)(圖1),并選擇湖中心樣點(diǎn)B15進(jìn)行垂直分層采樣,現(xiàn)場(chǎng)使用傳統(tǒng)掛錘式測(cè)得B15樣點(diǎn)水深為109.5 m,使用賽氏盤(pán)測(cè)定水體透明度為0.69 m;現(xiàn)場(chǎng)以1 m為間隔測(cè)定溫度,以≥0.1℃/m的溫度變化梯度為溫躍層深度(黃磊等,2015),以表層水溫改變≥0.5℃時(shí)所對(duì)應(yīng)的深度為混合層深度(黃亞男等,2017),以湖底以上0.5 m水深作為底層深度(中國(guó)科學(xué)院南京地理與湖泊研究所,2015),并繪制相應(yīng)的垂直變化曲線,據(jù)此確定垂直采樣深度。采樣點(diǎn)的基本信息見(jiàn)表1。

    1.3? ?樣品采集及鑒定

    浮游植物定性、定量樣品根據(jù)《淡水浮游生物研究方法》(章宗涉和黃祥飛,1991)進(jìn)行采集,物種鑒定參考相關(guān)文獻(xiàn)(遲若文,1990;中國(guó)科學(xué)院青藏高原綜合科學(xué)考察隊(duì),1992;朱蕙忠和陳嘉佑,2000;胡鴻鈞和魏印心,2006)。水深小于40 m的水樣使用有機(jī)玻璃采水器采集,水深40 m以上的水樣則使用卡蓋式采水器(QCC15)采集,河流水樣直接使用2 L采水器采集。其中1 L水樣用于定性樣品,現(xiàn)場(chǎng)使用25號(hào)浮游生物網(wǎng)過(guò)濾濃縮至50 mL,并立即用4%的甲醛溶液固定保存;1 L水樣用于定量樣品,現(xiàn)場(chǎng)立即用15 mL魯哥氏碘液固定,帶回實(shí)驗(yàn)室沉淀48 h后濃縮至30 mL。計(jì)數(shù)時(shí),移取0.1 mL濃縮勻液至0.1 mL計(jì)數(shù)框,在OLYMPUS CX23顯微鏡下進(jìn)行物種鑒定,每個(gè)樣品計(jì)數(shù)2~3次。浮游植物生物量(濕重) 通過(guò)體積換算間接獲得(Sun & Liu,2003;Hillebrand et al,1999;Kremer et al,2014)。

    pH、電導(dǎo)率(EC)、鹽度(Salt)及總?cè)芙恹}(TDS)使用PCSTestr35防水型多參數(shù)筆試測(cè)試儀(EUTECH,美國(guó)) 測(cè)定;水溫(WT)和電阻率(RES)使用HI98192微電腦電阻率測(cè)定儀測(cè)定;溶解氧(DO)使用HI98193微電腦溶解氧測(cè)定儀(HANNA,意大利)測(cè)定;濁度(TUR)使用HI98703微電腦多量程濁度測(cè)定儀(HANNA,意大利)測(cè)定;氨氮(NH4+-N)使用HI83399水質(zhì)檢測(cè)儀(HANNA,意大利)測(cè)定。

    1.4? ?數(shù)據(jù)處理與分析

    為了避免單一多樣性指數(shù)造成結(jié)果偏差,利用Gleason-Margalef 豐富度指數(shù)D(Margalef,1969)、Simpson優(yōu)勢(shì)度指數(shù)S(Simpson,1997)、Shannon-Wiener多樣性指數(shù)H′(Shannon & Weaver,1949)和Pielou 均勻度指數(shù)E(Pielou,1966)這4個(gè)指數(shù)來(lái)計(jì)算巴松錯(cuò)夏季浮游植物的多樣性。計(jì)算公式如下:

    式中:N為同一樣點(diǎn)中個(gè)體總數(shù);Ns為物種數(shù);Pi為第i種個(gè)體所占比例。

    采用Excel 2016和Origin 2022軟件進(jìn)行數(shù)據(jù)分析、統(tǒng)計(jì)和繪圖,采用ArcGIS 10.6軟件繪制地圖。用R X64 4.0.2 軟件包Vegan進(jìn)行BIOENV分析,環(huán)境因子均進(jìn)行l(wèi)g (x+1)轉(zhuǎn)換。使用CANOCO 5.0軟件對(duì)多樣性水平和垂直分布做去趨勢(shì)對(duì)應(yīng)分析(DCA),排序軸長(zhǎng)度均小于3,故采用冗余分析(RDA)。

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

    2.1? ?水體理化因子

    本研究從 15個(gè)樣點(diǎn)采集獲得25個(gè)湖水樣品,巴松錯(cuò)夏季水體理化因子的水平和垂直變化情況見(jiàn)圖 2和表2。

    2.2? ?巴松錯(cuò)浮游植物群落結(jié)構(gòu)

    2.2.1? ?物種組成? ?巴松錯(cuò)夏季共鑒定浮游植物167種,隸屬于5門(mén)8綱20目29科66屬,水平方向物種數(shù)大于垂直方向,物種組成均以硅藻門(mén)-綠藻門(mén)-藍(lán)藻門(mén)為主,裸藻門(mén)和隱藻門(mén)種類(lèi)較少,僅出現(xiàn)在水平方向(表3)。

    巴松錯(cuò)夏季浮游植物群落多集合韋恩圖如圖3。在水平方向上,各樣點(diǎn)浮游植物的共有種類(lèi)僅有1種,為肘狀針桿藻(Synedra ulna)。相比較于水平方向,在垂直方向上,各樣點(diǎn)浮游植物不存在共有種類(lèi)。

    2.2.2? ?細(xì)胞豐度和生物量? ?巴松措夏季浮游植物細(xì)胞豐度和生物量的水平和垂直變化如圖4。在水平方向上,細(xì)胞豐度的變化范圍是0.195×104~7.095×104 個(gè)/L,平均值為1.845×104 個(gè)/L,不同樣點(diǎn)間細(xì)胞豐度差異性較大,最高值在B13-0.5,為7.095×104 個(gè)/L,最低值在B12-0.5,為0.195×104 個(gè)/L;硅藻門(mén)細(xì)胞豐度占比相對(duì)最高,為22.22%~100%,其次是綠藻門(mén)和藍(lán)藻門(mén),細(xì)胞豐度變化范圍分別為0~0.628×104和0~1.256×104 個(gè)/L,占比為0~55%和0~38.89%;生物量為0.003~0.044 mg/L,平均值為0.012 mg/L,不同樣點(diǎn)間浮游植物生物量差異性較大,最高值在B13-0.5,最低值在B14-0.5;硅藻門(mén)生物量占比相對(duì)最高,變化范圍為10.98%~100%,其次是綠藻門(mén)和藍(lán)藻門(mén),生物量變化范圍分別為0~0.013和0~0.008 mg/L,占比為0~88.94%和0~27.49%;

    在垂直方向上,細(xì)胞豐度的變化范圍是0.21×104~2.79×104 個(gè)/L,平均值為1.006×104 個(gè)/L,不同垂直分層細(xì)胞豐度差異性較大,最高值在B15-100,為2.79×104 個(gè)/L,最低值在B15-109,為0.21×104 個(gè)/L,硅藻門(mén)細(xì)胞豐度占比相對(duì)最高,變化范圍為54.55%~100%,藍(lán)藻門(mén)僅在B15-109樣點(diǎn)出現(xiàn),占該樣點(diǎn)豐度的45.45%,綠藻門(mén)在B15-5和B15-15樣點(diǎn)出現(xiàn),分別占細(xì)胞豐度的6.25%和2.33%;生物量的變化為0.001~0.027 mg/L,平均值為0.006 mg/L,不同垂直分層浮游植物生物量差異性較大,最高值在B15-15,為0.027 mg/L,最低值在B15-25,為0.001 mg/L;硅藻門(mén)生物量占比相對(duì)最高,為31.81%~100%,藍(lán)藻門(mén)僅在B15-109樣點(diǎn)出現(xiàn),占該樣點(diǎn)生物量的68.19%,綠藻門(mén)在B15-5和B15-15樣點(diǎn)出現(xiàn),分別占生物量的11.56%和0.86%;

    2.2.3? ?多樣性指數(shù)? ?巴松措各站點(diǎn)夏季浮游植物物種多樣性的水平和垂直變化如圖5。在水平方向上,Margalef豐富度指數(shù)為0.756~2.807,均值1.734;Simpson優(yōu)勢(shì)度指數(shù)為0.695~0.921,均值為0.832;Shannon-Wiener多樣性指數(shù)為1.487~2.792,均值為2.28;Pielou均勻度指數(shù)為0.715~0.858,均值為0.806;在垂直方向上,Margalef豐富度指數(shù)為0.392~1.507,均值為0.955;Simpson優(yōu)勢(shì)度指數(shù)為0.583~0.881,均值為0.698;Shannon-Wiener多樣性指數(shù)為0.998~2.251,均值為1.643;Pielou均勻度指數(shù)為0.685~0.853,均值為0.739。總體上看,水平方向的Margalef豐富度指數(shù)、Simpson優(yōu)勢(shì)度指數(shù)、Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)均高于垂直方向。

    2.3? ?多樣性與環(huán)境因子的關(guān)系

    BIOENV分析發(fā)現(xiàn),pH、DO、EC、TDS、Salt和RES對(duì)巴松錯(cuò)夏季浮游植物群落多樣性水平分布的影響較大;EC和TUR對(duì)巴松錯(cuò)夏季浮游植物群落多樣性垂直分布的影響較大。

    為進(jìn)一步分析上述最優(yōu)環(huán)境因子對(duì)巴松錯(cuò)夏季浮游植物群落多樣性水平和垂直分布的影響,對(duì)群落多樣性與環(huán)境因子進(jìn)行RDA分析(圖6)。結(jié)果表明,在水平方向上,第1排序軸與第2排序軸對(duì)群落多樣性的方差的解釋率分別為48.45%和0.94%,群落多樣性與環(huán)境因子的相關(guān)性分別為0.7117和0.4762,前2軸累計(jì)方差為49.39%,上述環(huán)境因子對(duì)巴松錯(cuò)夏季浮游植物群落多樣性水平分布的影響大小順序?yàn)镋C、TDS、DO、Salt、pH和RES,其中EC(P=0.038,F(xiàn)=6.1)是顯著解釋性變量,對(duì)群落變異的解釋率為31.9%,是影響巴松錯(cuò)夏季浮游植物水平群落變異的主要環(huán)境因子。在垂直方向上,第1排序軸與第2排序軸對(duì)群落多樣性的方差的解釋率分別為57.67%和0.3%,群落多樣性與環(huán)境因子的相關(guān)性分別為0.7693和0.4424,前2軸累計(jì)方差為57.97%,上述環(huán)境因子對(duì)巴松錯(cuò)夏季浮游植物群落多樣性垂直分布的影響大小順序?yàn)門(mén)UR和EC,其中TUR(P=0.008,F(xiàn)=12.1)是顯著解釋性變量,對(duì)群落變異的解釋率為57.4%,是影響巴松錯(cuò)夏季浮游植物群落垂直變異的主要環(huán)境因子。

    3? ?討論

    3.1? ?浮游植物群落結(jié)構(gòu)特征

    本研究共鑒定出浮游植物167種,隸屬于5門(mén)8綱20目29科66屬,少于本課題組之前報(bào)道的242種(安瑞志等,2021),可能和2次調(diào)查的采樣周期、樣點(diǎn)范圍及分類(lèi)鑒定的詳細(xì)程度等因素有關(guān)。這也反映了巴松錯(cuò)浮游植物群落仍然處于動(dòng)態(tài)演替過(guò)程中。從浮游植物種類(lèi)來(lái)看,物種組成以硅藻占絕對(duì)優(yōu)勢(shì), 其次為綠藻和藍(lán)藻,這與本課題組2017-2018年間對(duì)巴松錯(cuò)浮游植物調(diào)查結(jié)果基本一致(安瑞志等,2021),主要原因可能是大多數(shù)硅藻都是狹冷性物種,適合在較寒冷和貧營(yíng)養(yǎng)的水環(huán)境中生活(潘成梅等,2021),而巴松錯(cuò)是一個(gè)過(guò)水型湖泊,其水位相較于吞吐型或封閉型湖泊更為穩(wěn)定(廖夢(mèng)娜等,2020),浮游植物常年在較低溫水環(huán)境下生存,已形成了動(dòng)態(tài)的適應(yīng)性機(jī)制,這一結(jié)果與色林錯(cuò)(劉洋等,2022)、浪錯(cuò)(扎西拉姆等,2021)的研究結(jié)果報(bào)道基本一致。

    巴松錯(cuò)夏季浮游植物群落的分布具有較明顯的水平和垂直空間差異,主要體現(xiàn)在物種數(shù)、細(xì)胞豐度和生物量等方面。結(jié)果表明,浮游植物垂直方向的物種數(shù)、細(xì)胞豐度和生物量均低于水平方向。在水平方向上,生物量最高值出現(xiàn)在湖心區(qū)域的B13-0.5樣點(diǎn),為0.044 mg/L,這主要是因?yàn)樵摌狱c(diǎn)EC(104.33 μS/cm)最高,高EC與有機(jī)質(zhì)降解和離子的輸入有著極大的聯(lián)系,能夠使藻類(lèi)的光合作用得到較好的補(bǔ)充,為浮游植物的大量繁殖提供了較好的營(yíng)養(yǎng)條件(李鋒,2021)。在垂直方向上,生物量最大值出現(xiàn)在B15-15,為0.027 mg/L,主要原因可能是巴松錯(cuò)為深水湖泊, 在夏季具有明顯的溫躍層,水體混合交換弱,水深15 m是巴松錯(cuò)浮游植物最適的生長(zhǎng)區(qū)域,形成了浮游植物上下層分布的明顯差異,這與潘繼征等(2009)的研究和邢冰偉等(2020)的研究結(jié)果基本一致。

    巴松錯(cuò)夏季浮游植物群落的的水平和垂直分布差異不僅體現(xiàn)在物種數(shù)、細(xì)胞豐度和生物量上,還體現(xiàn)在物種多樣性上。物種多樣性是生物與環(huán)境因子長(zhǎng)期作用和適應(yīng)的結(jié)果,用來(lái)衡量浮游植物群落結(jié)構(gòu)是否穩(wěn)定(邱小琮等,2011)。基于浮游植物對(duì)湖泊生態(tài)系統(tǒng)群落穩(wěn)定性進(jìn)行生物學(xué)評(píng)價(jià)有很多方法,本文選擇了應(yīng)用最為廣泛的Margalef豐富度指數(shù)、Simpson優(yōu)勢(shì)度指數(shù)、Shannon-Wiener多樣性指數(shù)和Pielou 均勻度指數(shù)這4個(gè)指數(shù)進(jìn)行評(píng)價(jià)分析。結(jié)果表明,水平方向的4個(gè)指數(shù)值均高于垂直方向,表明巴松錯(cuò)夏季浮游植物群落結(jié)構(gòu)水平分布較為均勻,群落復(fù)雜程度相對(duì)較高,較垂直分布更穩(wěn)定,這與柯志新等(2011)和籍辰旭等(2021)的研究結(jié)果基本一致。

    3.2? ?浮游植物多樣性與環(huán)境因子的關(guān)系

    環(huán)境因子是浮游植物群落分布的主要影響因素,水體透明度、濁度、溶解氧、EC、pH、水溫和鹽度等環(huán)境因子都對(duì)浮游植物群落分布產(chǎn)生影響(趙秀俠等,2018)。本研究中,BIOENV分析結(jié)果顯示,在水平方向上,pH、DO、EC、TDS、Salt和RES對(duì)巴松錯(cuò)夏季浮游植物物種多樣性的影響較大;在垂直方向上,EC和TUR對(duì)巴松錯(cuò)夏季浮游植物物種多樣性垂直分布的影響較大。但不同環(huán)境因子對(duì)浮游植物群落結(jié)構(gòu)變化的影響也存在差異。王宇飛等(2015)的研究發(fā)現(xiàn)影響漢豐湖夏季浮游植物群落結(jié)構(gòu)特征的主要因素是溫度、DO、pH、TN、TP、 DN和DP;牛遠(yuǎn)等(2016)的研究發(fā)現(xiàn)環(huán)境變量較好地解釋了撫仙湖夏季熱分層時(shí)期浮游植物群落組成的空間變化,WT、EC和pH是最關(guān)鍵影響因子,而營(yíng)養(yǎng)鹽的作用并不顯著。本研究中,RDA分析結(jié)果表明,EC是影響巴松錯(cuò)夏季浮游植物水平群落變異的主要環(huán)境因子,TUR是影響垂直變異的主要環(huán)境因子。EC與水體營(yíng)養(yǎng)鹽含量有關(guān)系,而營(yíng)養(yǎng)鹽為浮游植物提供充足的營(yíng)養(yǎng),EC通過(guò)影響浮游植物的代謝活動(dòng)和繁殖從而影響浮游植物的群落結(jié)構(gòu)(趙秀俠等,2021)。而濁度(TUR)是由于水中含有泥沙、黏土、有機(jī)物、無(wú)機(jī)物、浮游生物等懸浮物質(zhì)所造成的,可使光散射或吸收,影響浮游植物對(duì)光的吸收(聶小芬等,2021)。本研究采樣時(shí)間為夏季,降雨徑流所攜帶的陸源物質(zhì)發(fā)生沉降,也會(huì)引起水體濁度升高,透明度降低,光合作用減弱,在垂向上易于產(chǎn)生生物分層,對(duì)浮游植物的垂直分布具有一定的抑制作用(安瑞志等,2021)。因此,浮游植物多樣性的變化是多種環(huán)境因子在水平與垂直分布共同作用的結(jié)果。

    本研究只涉及到1個(gè)季節(jié)以及水體中的部分理化因子, 不能得出巴松錯(cuò)浮游植物水平與垂直分布的季節(jié)動(dòng)態(tài)變化,今后有待于加強(qiáng)對(duì)其季節(jié)動(dòng)態(tài)變化研究,并進(jìn)行初級(jí)生產(chǎn)力估算的研究,以便為碳中和背景下高寒湖泊固碳潛力的計(jì)量及估算提供數(shù)據(jù)支持。

    參考文獻(xiàn)

    安瑞志,潘成梅,劉洋,等,2020. 巴松措(湖)夏季浮游植物優(yōu)勢(shì)種群分布格局與共存機(jī)制研究[J].高原科學(xué)研究,4(2):27-36.

    安瑞志,潘成梅,塔巴拉珍,等,2021. 西藏巴松錯(cuò)浮游植物功能群垂直分布特征及其與環(huán)境因子的關(guān)系[J].湖泊科學(xué),33(1):86-101.

    陳立婧,楊菲,吳淑賢,等,2013. 西藏那曲地區(qū)鹽湖浮游植物群落結(jié)構(gòu)的特征[J].上海海洋大學(xué)學(xué)報(bào),22(4):577-585.

    遲若文,1990.西藏硅藻圖集[M].拉薩:西藏人民出版社.

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

    黃磊,王君波,朱立平等,2015. 納木錯(cuò)水溫變化及熱力學(xué)分層特征初步研究[J].湖泊科學(xué),27(4):711-718.

    黃亞男,紀(jì)道斌,龍良紅,等,2017. 三峽庫(kù)區(qū)典型支流春季特征及其水華優(yōu)勢(shì)種差異分析[J].長(zhǎng)江流域資源與環(huán)境,26(3):461-470.

    籍辰旭,陳逍遙,陳火榮,等,2021. 三沙灣春季浮游植物群落結(jié)構(gòu)特征[J].漁業(yè)研究,43(1):42-52.

    靳錚,張雪芹,次旦央宗,2020. 藏南羊卓雍錯(cuò)湖面非封凍期溫濕變化及輻射平衡分析[J].干旱區(qū)研究,37(4):947-955.

    柯志新,黃良民,譚燁輝,等,2011. 2007年夏季南海北部浮游植物的物種組成及豐度分布[J].熱帶海洋學(xué)報(bào),30(1):131-143.

    廖夢(mèng)娜,金伊麗,李晨瑜,等,2020. 藏東南巴松錯(cuò)200年沉積過(guò)程及其對(duì)硅藻記錄的影響[J].生態(tài)學(xué)報(bào),40(3):1089-1100.

    李鋒,2021. 二灘水庫(kù)浮游植物群落結(jié)構(gòu)與水環(huán)境因子時(shí)空變化及其相關(guān)性研究[D].重慶:西南大學(xué).

    劉洋,安瑞志,巴桑,2022. 西藏色林錯(cuò)-普若崗日國(guó)家公園潛在建設(shè)區(qū)夏季浮游植物群落生態(tài)特征[J].水生態(tài)學(xué)雜志,43(6):51-58.

    聶小芬,楊正健,王從鋒,等,2021. 瀾滄江梯級(jí)水庫(kù)建設(shè)下水體電導(dǎo)率和濁度時(shí)空分布特征[J].三峽大學(xué)學(xué)報(bào)(自然科學(xué)版),43(1):18-24.

    牛遠(yuǎn),孔祥虹,余輝,等,2016. 撫仙湖夏季熱分層時(shí)期浮游植物空間分布特征[J].生態(tài)學(xué)雜志,35(7):1865-1871.

    潘成梅,劉洋,安瑞志,等,2021.西藏麥地卡濕地的浮游植物——1.優(yōu)勢(shì)種的時(shí)空生態(tài)位[J].湖泊科學(xué),33(6):1805-1819.

    潘繼征,熊飛,李文朝,等,2009. 撫仙湖浮游植物群落結(jié)構(gòu)、分布及其影響因子[J].生態(tài)學(xué)報(bào),29(10):5376-5385.

    邱小琮,趙紅雪,2011.寧夏沙湖浮游植物群落結(jié)構(gòu)及多樣性研究[J].水生態(tài)學(xué)雜志,32(1):20-26.

    孫萍,李艷,潘玉龍,等,2020. 熱帶東印度洋春季浮游植物群落結(jié)構(gòu)空間特征分析[J].海洋學(xué)報(bào),42(8):76-88.

    邢冰偉,徐季雄,曹玥,等,2020. 九寨溝國(guó)家級(jí)自然保護(hù)區(qū)長(zhǎng)海夏季浮游植物群落結(jié)構(gòu)及生態(tài)評(píng)價(jià)[J].湖泊科學(xué),32(4):1088-1099.

    王捷,李博,馮佳,等,2015. 西藏西南部湖泊浮游藻類(lèi)區(qū)系及群落結(jié)構(gòu)特征[J].水生生物學(xué)報(bào),39(4):837-844.

    王宇飛,趙秀蘭,何丙輝,等,2015. 漢豐湖夏季浮游植物群落與環(huán)境因子的典范對(duì)應(yīng)分析[J].環(huán)境科學(xué),36(3):922-927.

    楊欣蘭,潘瑛子,巴桑,2022.巴松措原生動(dòng)物群落結(jié)構(gòu)季節(jié)及垂直分布[J].生態(tài)學(xué)報(bào),42(08):3216-3227.

    趙秀俠,方婷,楊坤,等,2018.安徽沱湖夏季浮游植物群落結(jié)構(gòu)特征與環(huán)境因子關(guān)系[J].植物科學(xué)學(xué)報(bào),36(5):687-695.

    趙秀俠,盧文軒,梁陽(yáng)陽(yáng),等,2021.安徽三座通江湖泊秋季浮游植物群落結(jié)構(gòu)特征及其影響因子[J].生態(tài)學(xué)雜志,40(1):67-75.

    扎西拉姆,陳美群,牟振波,等,2022. 西藏浪錯(cuò)夏季浮游植物組成與環(huán)境因子的關(guān)系[J].水產(chǎn)科學(xué),41(1):137-142.

    中國(guó)科學(xué)院青藏高原綜合科學(xué)考察隊(duì),1981.西藏水利[M].北京:科學(xué)出版社.

    中國(guó)科學(xué)院青藏高原綜合科學(xué)考察隊(duì),1984.西藏河流與湖泊[M].北京:科學(xué)出版社.

    中國(guó)科學(xué)院青藏高原綜合科學(xué)考察隊(duì),1992.西藏藻類(lèi)[M].北京:科學(xué)出版社.

    中國(guó)科學(xué)院南京地理與湖泊研究所,2015. 湖泊調(diào)查技術(shù)規(guī)程[M].北京:科學(xué)出版社.

    朱蕙忠,陳嘉佑,2000.中國(guó)西藏硅藻[M].北京:科學(xué)出版社.

    Hillebrand H,Dürselen C D,Kirschtel D,et al,1999. Biovolume calculation for pelagic and benthic microalgae[J]. Journal of Phycology,35(2): 403-424.

    Kremer C T,Gillette J P,Rudstam L G,et al,2014. A compendium of cell and natural unit biovolumes for >1200 freshwater phytoplankton species[J]. Ecology,95(10): 2984.

    Li Z X,Gao Y,Wang S Y,et al,2021. Phytoplankton community response to nutrients along lake salinity and altitude gradients on the Qinghai-Tibet Plateau[J]. Ecological Indicators,128: 107848.

    Margalef R,1969. Diversity and Stability in Ecological Systems[M]. Upton: Brookhacen National Laboratory:13-26.

    Pielou E C,1966.The measurements of diversity in different types of biological collections[J].Journal of Theoretical Biology,13:131-144.

    Ptacnik R,Lepist L,E Willén,et al,2008. Quantitative responses of lake phytoplankton to eutrophication in Northern Europe[J]. Aquatic Ecology,42(2):227-236.

    Shannon C E,Weaver W,1949.The mathematical theory of communication[M]. Urbana:University of Illinois Press.

    Simpson E H,1997. Measurement of Diversity[J]. Journal of Cardiothoracic and Vascular Anesthesia,11(6):812-812.

    Sun J,Liu D Y,2003. Geometric models for calculating cell biovolume and surface area for phytoplankton. Journal of Plankton Research,25(11): 1331-1346.

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

    Horizontal and Vertical Distributions of the Summer Phytoplankton

    Community in Basomtso Lake of Xizang, China

    LIU Yang1, 2, AN Rui‐zhi1,2, WANG Chen1,2, PAN Cheng‐mei1,2, BA Sang1,2

    (1. Center for Carbon Neutrality in the Earth's Third Pole, Tibet University, Lhasa? ?850000, P.R. China;

    2. Laboratory of Wetland and Catchment Ecosystem in Tibetan Plateau, School of Ecology

    and Environment, Tibet University, Lhasa? ?850000, P.R. China)

    Abstract: In this study, we explored the horizontal and vertical distributions of the phytoplankton community and its relationship with environmental factors in Basomtso Lake of Tibet, aiming to provide basic data for the further study of plateau lake ecosystems. In July 2020, horizontal water samples were collected at 15 sampling sites (B1-15) and vertical samples were collected at 11 depths at B15, the sampling site near the center of Basomtso Lake. Physical and chemical parameters of the lake water were also measured. A total of 275 water samples were collected for phytoplankton analysis, focusing on species composition, cell abundance and biomass, and species diversity. The response of phytoplankton to changes in horizontal and vertical environmental gradients were analyzed by redundancy analysis (RDA). Results show that: (1) A total of 167 phytoplankton species from 66 genera, 29 families, 20 orders, 8 classes and 5 phyla were identified, with dominance by bacillariophyta-chlorophyta-cyanophyta. Species richness of the horizontal distribution was higher than that of the vertical distribution. (2) Horizontally, the average cell abundance and biomass were 1.845×104 cells/L and 0.012 mg/L, and vertically, the average cell abundance and biomass were 1.006×104 cells/L and 0.006 mg/L. (3) The Margalef richness, Simpson dominance, Shannon-Wiener diversity and Pielou evenness indices of the summer phytoplankton community in the horizontal distribution were all higher than those of the phytoplankton community in the vertical distribution. These indices show that the horizontal distribution was more uniform, and that community complexity and stability were higher than in the vertical distribution. (4) RDA showed that electrical conductivity (EC) was the primary environmental factor associated with the horizontal distribution of phytoplankton, and turbidity (TUR) was the primary environmental factor associated with the vertical distribution of phytoplankton in Basomtso Lake. The results of this study provide basic data for further study of alpine lake ecosystems.

    Key words:phytoplankton; community structure;horizontal distribution; vertical distribution;

    Basomtso Lake

    收稿日期:2021-12-17? ? ? 修回日期:2023-09-25

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(32070418,31660620);西藏大學(xué)培育基金(ZDQMJH19-34);2022 年中央財(cái)政支持地方高校改革發(fā)展專(zhuān)項(xiàng)資金 (藏財(cái)預(yù)指[2022]1號(hào))。

    作者簡(jiǎn)介:劉洋,女,1988年生,實(shí)驗(yàn)師,主要從事高寒濕地浮游植物生態(tài)學(xué)研究。E-mail:ly@utibet.edu.cn

    通信作者:巴桑,男,1970年生,教授,主要從事高寒濕地與流域生態(tài)系統(tǒng)研究。E-mail:hbasang2003@aliyun.com

    猜你喜歡
    浮游植物群落結(jié)構(gòu)
    海南省紅樹(shù)林區(qū)浮游動(dòng)物多樣性的初步研究
    海南省紅樹(shù)林區(qū)大型底棲動(dòng)物多樣性的初步研究
    塔里木河阿拉爾段及多浪水庫(kù)浮游生物的調(diào)查研究
    遼東灣網(wǎng)采浮游植物名錄
    特征色素在淡水魚(yú)養(yǎng)殖模式中對(duì)浮游植物群落結(jié)構(gòu)的影響
    渭河陜西段浮游植物群落結(jié)構(gòu)及水質(zhì)評(píng)價(jià)
    濟(jì)南流域春季浮游植物多樣性的研究
    峽江樞紐庫(kù)區(qū)水生生物調(diào)查及水質(zhì)評(píng)價(jià)
    濟(jì)南流域輪蟲(chóng)群落結(jié)構(gòu)及多樣性的研究
    黃河蘭州市區(qū)段浮游動(dòng)物群落結(jié)構(gòu)調(diào)查及水質(zhì)初步評(píng)價(jià)
    亚洲欧美日韩无卡精品| 中文字幕高清在线视频| av超薄肉色丝袜交足视频| 欧美激情 高清一区二区三区| 午夜福利高清视频| 女性被躁到高潮视频| 俄罗斯特黄特色一大片| 最新美女视频免费是黄的| 久久久久久久久久久久大奶| 国产午夜精品久久久久久| 日本欧美视频一区| 亚洲av成人一区二区三| 国产97色在线日韩免费| 国产精品久久久久久人妻精品电影| 在线观看一区二区三区| 国产成人欧美| 一级毛片精品| 亚洲成国产人片在线观看| 大香蕉久久成人网| 欧美日韩中文字幕国产精品一区二区三区 | 久热这里只有精品99| 久久精品国产综合久久久| 午夜久久久在线观看| 久久青草综合色| 国产精品久久久久久人妻精品电影| 亚洲人成网站在线播放欧美日韩| 真人做人爱边吃奶动态| 久久人人爽av亚洲精品天堂| www.精华液| 两性夫妻黄色片| 亚洲av成人av| 久久久国产欧美日韩av| 动漫黄色视频在线观看| 国产免费男女视频| 怎么达到女性高潮| 老司机午夜福利在线观看视频| 日韩视频一区二区在线观看| 亚洲片人在线观看| 亚洲三区欧美一区| 国产99久久九九免费精品| 免费不卡黄色视频| 久久精品91无色码中文字幕| 天堂√8在线中文| 日韩成人在线观看一区二区三区| 国语自产精品视频在线第100页| 亚洲欧美日韩高清在线视频| 在线av久久热| 一区二区三区精品91| 叶爱在线成人免费视频播放| 美女大奶头视频| 美女国产高潮福利片在线看| 亚洲 欧美一区二区三区| 伦理电影免费视频| 一级,二级,三级黄色视频| 久久欧美精品欧美久久欧美| 欧美成人午夜精品| 亚洲人成网站在线播放欧美日韩| 侵犯人妻中文字幕一二三四区| 免费在线观看日本一区| 最新在线观看一区二区三区| 两人在一起打扑克的视频| 欧美日韩福利视频一区二区| 欧美日韩亚洲综合一区二区三区_| 日韩大尺度精品在线看网址 | 丁香六月欧美| 窝窝影院91人妻| 欧美一级a爱片免费观看看 | 在线观看舔阴道视频| 精品午夜福利视频在线观看一区| 久久久久国内视频| 久久久国产精品麻豆| 91精品三级在线观看| 91麻豆精品激情在线观看国产| 亚洲五月色婷婷综合| 好看av亚洲va欧美ⅴa在| 免费在线观看黄色视频的| 叶爱在线成人免费视频播放| 亚洲av熟女| 国产亚洲av嫩草精品影院| 国产成人影院久久av| 午夜福利在线观看吧| 国产熟女午夜一区二区三区| 欧美成人一区二区免费高清观看 | 天天躁夜夜躁狠狠躁躁| 丝袜在线中文字幕| av天堂在线播放| 国产精品电影一区二区三区| 午夜精品久久久久久毛片777| 久久精品91无色码中文字幕| 亚洲欧美一区二区三区黑人| 久久久精品欧美日韩精品| 亚洲av美国av| 99国产极品粉嫩在线观看| 97人妻天天添夜夜摸| 亚洲欧美日韩另类电影网站| 国产精品1区2区在线观看.| 精品少妇一区二区三区视频日本电影| 国产成人av激情在线播放| 女人被狂操c到高潮| 一区在线观看完整版| 欧美在线一区亚洲| 亚洲av日韩精品久久久久久密| 91麻豆精品激情在线观看国产| 久久国产精品男人的天堂亚洲| 亚洲专区中文字幕在线| 在线观看免费视频日本深夜| 天天一区二区日本电影三级 | 真人一进一出gif抽搐免费| 日本免费a在线| 啦啦啦 在线观看视频| 午夜a级毛片| 国产野战对白在线观看| 男人舔女人的私密视频| 麻豆久久精品国产亚洲av| 精品少妇一区二区三区视频日本电影| 色老头精品视频在线观看| 91九色精品人成在线观看| 日日摸夜夜添夜夜添小说| 国产成人av激情在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看免费日韩欧美大片| 国产国语露脸激情在线看| 中文字幕色久视频| 欧美日韩亚洲国产一区二区在线观看| 欧美日韩一级在线毛片| 久久亚洲真实| av网站免费在线观看视频| 制服丝袜大香蕉在线| 一边摸一边抽搐一进一小说| bbb黄色大片| 精品一区二区三区四区五区乱码| 精品熟女少妇八av免费久了| 国产欧美日韩一区二区精品| 老司机午夜十八禁免费视频| 精品久久久久久成人av| 在线播放国产精品三级| 一进一出抽搐动态| 97碰自拍视频| 99久久99久久久精品蜜桃| 久久青草综合色| 国产激情久久老熟女| 在线观看午夜福利视频| 日韩成人在线观看一区二区三区| 午夜两性在线视频| 男男h啪啪无遮挡| 50天的宝宝边吃奶边哭怎么回事| 欧洲精品卡2卡3卡4卡5卡区| 亚洲欧美精品综合久久99| 亚洲五月婷婷丁香| 亚洲aⅴ乱码一区二区在线播放 | 韩国精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 乱人伦中国视频| 一边摸一边做爽爽视频免费| 丰满人妻熟妇乱又伦精品不卡| 日本vs欧美在线观看视频| 无人区码免费观看不卡| 很黄的视频免费| 男人操女人黄网站| 狠狠狠狠99中文字幕| 极品教师在线免费播放| 一a级毛片在线观看| 国产免费av片在线观看野外av| 精品久久久久久,| 亚洲电影在线观看av| 老司机深夜福利视频在线观看| 欧美一级a爱片免费观看看 | 国内久久婷婷六月综合欲色啪| 午夜福利视频1000在线观看 | 日本 欧美在线| 长腿黑丝高跟| 亚洲国产精品sss在线观看| 少妇粗大呻吟视频| 少妇 在线观看| 韩国av一区二区三区四区| 夜夜躁狠狠躁天天躁| 欧美不卡视频在线免费观看 | 国产99白浆流出| 黑人操中国人逼视频| 午夜成年电影在线免费观看| 欧美成狂野欧美在线观看| 99在线视频只有这里精品首页| 制服诱惑二区| 一区在线观看完整版| 一个人免费在线观看的高清视频| 国产三级在线视频| 一区二区三区激情视频| 午夜福利影视在线免费观看| 欧美乱色亚洲激情| 夜夜躁狠狠躁天天躁| 啦啦啦观看免费观看视频高清 | 欧美黑人欧美精品刺激| 一区福利在线观看| 亚洲欧美日韩无卡精品| 9色porny在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产单亲对白刺激| aaaaa片日本免费| 亚洲全国av大片| 国产在线精品亚洲第一网站| 可以在线观看毛片的网站| 天天一区二区日本电影三级 | 亚洲欧美一区二区三区黑人| 久久人妻福利社区极品人妻图片| 亚洲无线在线观看| 久久精品91无色码中文字幕| 国产av又大| 国产成人影院久久av| 欧美人与性动交α欧美精品济南到| www.www免费av| 69av精品久久久久久| 亚洲久久久国产精品| 亚洲五月色婷婷综合| 日韩国内少妇激情av| 久久久久久久久中文| 在线国产一区二区在线| 咕卡用的链子| 亚洲成国产人片在线观看| 久久久久国内视频| 男女下面插进去视频免费观看| 日本a在线网址| 欧美日韩乱码在线| 午夜福利免费观看在线| 亚洲激情在线av| 在线天堂中文资源库| 国产蜜桃级精品一区二区三区| 色在线成人网| 神马国产精品三级电影在线观看 | 在线视频色国产色| 婷婷精品国产亚洲av在线| ponron亚洲| 欧美亚洲日本最大视频资源| 亚洲精品国产区一区二| 亚洲电影在线观看av| 亚洲人成网站在线播放欧美日韩| 好看av亚洲va欧美ⅴa在| 午夜日韩欧美国产| 亚洲中文日韩欧美视频| 黄色毛片三级朝国网站| 国产免费男女视频| 亚洲午夜精品一区,二区,三区| 中文字幕高清在线视频| 国产片内射在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品国产色婷婷电影| 中出人妻视频一区二区| 非洲黑人性xxxx精品又粗又长| 人人妻人人澡人人看| 老司机福利观看| 中文字幕久久专区| 国产精品野战在线观看| 免费在线观看日本一区| 中出人妻视频一区二区| 非洲黑人性xxxx精品又粗又长| 亚洲人成电影观看| 可以在线观看的亚洲视频| 99香蕉大伊视频| 国产极品粉嫩免费观看在线| 亚洲av第一区精品v没综合| 自线自在国产av| 久久久久久人人人人人| 日本黄色视频三级网站网址| 久久精品亚洲精品国产色婷小说| 无限看片的www在线观看| 亚洲精品国产区一区二| 久久人人97超碰香蕉20202| 91字幕亚洲| 久久久久久人人人人人| 精品电影一区二区在线| 日韩精品青青久久久久久| av天堂在线播放| 成人三级黄色视频| 美女午夜性视频免费| 日韩欧美一区视频在线观看| 97人妻精品一区二区三区麻豆 | 亚洲av电影在线进入| 成人亚洲精品一区在线观看| 黄频高清免费视频| 国产熟女xx| 国产亚洲精品一区二区www| 身体一侧抽搐| 日韩精品青青久久久久久| 亚洲精品国产区一区二| 神马国产精品三级电影在线观看 | 国产精品一区二区精品视频观看| 老熟妇乱子伦视频在线观看| 亚洲av成人一区二区三| 午夜免费鲁丝| 欧美国产日韩亚洲一区| 午夜久久久久精精品| 国产成人av激情在线播放| 看免费av毛片| 9191精品国产免费久久| 又黄又粗又硬又大视频| 老汉色∧v一级毛片| 99久久久亚洲精品蜜臀av| 久久久久亚洲av毛片大全| 婷婷精品国产亚洲av在线| 国产成人系列免费观看| 久久精品aⅴ一区二区三区四区| av视频免费观看在线观看| 亚洲成av片中文字幕在线观看| 久久久国产精品麻豆| 久久人妻熟女aⅴ| av超薄肉色丝袜交足视频| x7x7x7水蜜桃| 国产国语露脸激情在线看| 欧美色欧美亚洲另类二区 | 欧美绝顶高潮抽搐喷水| 国产精品一区二区精品视频观看| 精品福利观看| 亚洲熟女毛片儿| 亚洲伊人色综图| av在线播放免费不卡| 亚洲成av片中文字幕在线观看| 亚洲激情在线av| 欧美午夜高清在线| 国产成人精品在线电影| 色精品久久人妻99蜜桃| 国产亚洲欧美98| 丰满人妻熟妇乱又伦精品不卡| 可以在线观看毛片的网站| 欧美成人午夜精品| 国产三级在线视频| 啦啦啦 在线观看视频| 久久精品国产亚洲av高清一级| 久久中文看片网| 精品人妻在线不人妻| 在线观看一区二区三区| 他把我摸到了高潮在线观看| 成在线人永久免费视频| 在线观看免费午夜福利视频| 国产97色在线日韩免费| 久9热在线精品视频| 久久午夜综合久久蜜桃| 在线国产一区二区在线| 国产精品亚洲一级av第二区| 欧美日韩黄片免| 成人三级黄色视频| 亚洲中文av在线| 国产精品久久久久久人妻精品电影| 日韩av在线大香蕉| 视频在线观看一区二区三区| 日本 欧美在线| 可以在线观看的亚洲视频| 国内久久婷婷六月综合欲色啪| 黄片播放在线免费| 欧美日韩精品网址| 波多野结衣高清无吗| 亚洲男人天堂网一区| 天堂动漫精品| 欧美日韩乱码在线| 亚洲精华国产精华精| 国产区一区二久久| 亚洲中文日韩欧美视频| 精品卡一卡二卡四卡免费| 午夜两性在线视频| 好男人电影高清在线观看| 国产一区在线观看成人免费| 欧美黄色片欧美黄色片| 亚洲精品中文字幕一二三四区| 久99久视频精品免费| 亚洲在线自拍视频| 亚洲专区字幕在线| 老汉色∧v一级毛片| 久久狼人影院| or卡值多少钱| 国产熟女午夜一区二区三区| 法律面前人人平等表现在哪些方面| 激情在线观看视频在线高清| 欧美亚洲日本最大视频资源| 人人妻人人澡欧美一区二区 | 免费看美女性在线毛片视频| 十八禁人妻一区二区| 亚洲欧美精品综合久久99| 亚洲 欧美一区二区三区| 999久久久精品免费观看国产| 在线永久观看黄色视频| 母亲3免费完整高清在线观看| 人人澡人人妻人| 亚洲五月色婷婷综合| av欧美777| 中文字幕另类日韩欧美亚洲嫩草| a在线观看视频网站| 久久 成人 亚洲| 男男h啪啪无遮挡| 成人国产一区最新在线观看| 人人妻人人爽人人添夜夜欢视频| 这个男人来自地球电影免费观看| aaaaa片日本免费| 日日夜夜操网爽| 亚洲精品中文字幕在线视频| 欧美激情 高清一区二区三区| 亚洲精华国产精华精| 大香蕉久久成人网| 国产成人啪精品午夜网站| 精品国产乱子伦一区二区三区| 久久精品亚洲熟妇少妇任你| 国产高清视频在线播放一区| 久久中文看片网| 视频在线观看一区二区三区| 精品久久久久久成人av| 国产1区2区3区精品| 手机成人av网站| 老汉色∧v一级毛片| 精品国产美女av久久久久小说| 久久久国产精品麻豆| 大香蕉久久成人网| 欧美久久黑人一区二区| 亚洲精品一卡2卡三卡4卡5卡| 久久国产精品男人的天堂亚洲| 亚洲无线在线观看| 69精品国产乱码久久久| 搡老熟女国产l中国老女人| 女人高潮潮喷娇喘18禁视频| 日韩av在线大香蕉| 两个人看的免费小视频| 精品国产亚洲在线| 久久精品人人爽人人爽视色| 中文字幕色久视频| 国产视频一区二区在线看| 啦啦啦 在线观看视频| 欧美丝袜亚洲另类 | 免费少妇av软件| 淫秽高清视频在线观看| 中国美女看黄片| 欧美国产日韩亚洲一区| 免费无遮挡裸体视频| 大型黄色视频在线免费观看| 国产精品久久久久久精品电影 | avwww免费| 亚洲精品中文字幕一二三四区| 国产色视频综合| 最近最新中文字幕大全电影3 | 欧美黑人精品巨大| 两个人免费观看高清视频| 母亲3免费完整高清在线观看| 午夜精品久久久久久毛片777| 欧洲精品卡2卡3卡4卡5卡区| 欧美乱色亚洲激情| 好男人电影高清在线观看| 亚洲精品国产精品久久久不卡| 嫩草影视91久久| 男女床上黄色一级片免费看| 久久香蕉激情| 欧美日本视频| 国产精品影院久久| 大型av网站在线播放| 国产激情久久老熟女| 欧美 亚洲 国产 日韩一| 亚洲av成人不卡在线观看播放网| 国产一区二区在线av高清观看| 久久久久国内视频| 亚洲午夜理论影院| 人人妻人人澡欧美一区二区 | 老司机午夜十八禁免费视频| 亚洲av片天天在线观看| 午夜视频精品福利| 女生性感内裤真人,穿戴方法视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩精品网址| 在线永久观看黄色视频| www.www免费av| 多毛熟女@视频| 91九色精品人成在线观看| 亚洲全国av大片| 中文字幕人妻丝袜一区二区| 禁无遮挡网站| 制服人妻中文乱码| 久热爱精品视频在线9| av天堂在线播放| 国产午夜精品久久久久久| 国产欧美日韩综合在线一区二区| 人妻丰满熟妇av一区二区三区| 亚洲三区欧美一区| 亚洲片人在线观看| 欧美成人性av电影在线观看| 精品欧美一区二区三区在线| 中文字幕av电影在线播放| 亚洲午夜理论影院| 成人特级黄色片久久久久久久| 成人免费观看视频高清| 精品熟女少妇八av免费久了| 午夜免费鲁丝| 午夜激情av网站| 在线天堂中文资源库| 99在线视频只有这里精品首页| 日本在线视频免费播放| 老熟妇乱子伦视频在线观看| 无遮挡黄片免费观看| 免费无遮挡裸体视频| 久久久久久人人人人人| 99精品欧美一区二区三区四区| 国产亚洲精品久久久久5区| 亚洲成a人片在线一区二区| 99国产精品99久久久久| 黑人巨大精品欧美一区二区蜜桃| 国产精品一区二区三区四区久久 | 免费高清在线观看日韩| 女警被强在线播放| 免费在线观看日本一区| 老汉色∧v一级毛片| 国产精品免费一区二区三区在线| 色精品久久人妻99蜜桃| 国产激情久久老熟女| 亚洲欧美精品综合一区二区三区| 久久狼人影院| 高清黄色对白视频在线免费看| 麻豆国产av国片精品| 久久影院123| 国产亚洲av高清不卡| 99精品欧美一区二区三区四区| 后天国语完整版免费观看| 久久伊人香网站| 国产精品98久久久久久宅男小说| 一进一出好大好爽视频| 又大又爽又粗| 可以免费在线观看a视频的电影网站| 亚洲最大成人中文| 精品人妻在线不人妻| 一边摸一边抽搐一进一小说| 午夜久久久久精精品| ponron亚洲| 免费高清视频大片| 可以在线观看的亚洲视频| 成人永久免费在线观看视频| 国产一区在线观看成人免费| 国产色视频综合| 国产成人精品久久二区二区免费| 日本五十路高清| 欧美成人性av电影在线观看| 欧美另类亚洲清纯唯美| 国产99白浆流出| 精品乱码久久久久久99久播| 欧美色视频一区免费| 欧美日韩瑟瑟在线播放| 99国产精品免费福利视频| 日本一区二区免费在线视频| 欧美性长视频在线观看| 欧美黑人精品巨大| 久久香蕉激情| 精品高清国产在线一区| 精品一品国产午夜福利视频| 在线观看免费日韩欧美大片| 女人爽到高潮嗷嗷叫在线视频| 日韩精品免费视频一区二区三区| 免费看美女性在线毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 午夜两性在线视频| or卡值多少钱| 美女扒开内裤让男人捅视频| 999精品在线视频| 久久中文字幕人妻熟女| 中出人妻视频一区二区| 亚洲自偷自拍图片 自拍| 国产一卡二卡三卡精品| 午夜老司机福利片| 国内精品久久久久精免费| 好男人电影高清在线观看| 99国产综合亚洲精品| 国产成人啪精品午夜网站| 老熟妇仑乱视频hdxx| 满18在线观看网站| 国产aⅴ精品一区二区三区波| 999久久久精品免费观看国产| 国产精品av久久久久免费| 久久人妻av系列| 日韩欧美国产一区二区入口| 久久精品91蜜桃| 99riav亚洲国产免费| 精品久久久久久久人妻蜜臀av | 黑人巨大精品欧美一区二区蜜桃| 亚洲五月天丁香| 91在线观看av| 中文字幕人妻丝袜一区二区| 国产97色在线日韩免费| 国产精品久久久久久亚洲av鲁大| 欧美成人一区二区免费高清观看 | 亚洲精品一区av在线观看| 欧美激情极品国产一区二区三区| 欧美一级毛片孕妇| 国产在线观看jvid| 久久久久国产一级毛片高清牌| 精品久久久久久,| 老司机午夜福利在线观看视频| 国产在线观看jvid| 婷婷丁香在线五月| 精品一品国产午夜福利视频| 国产成人啪精品午夜网站| 两个人免费观看高清视频| 高清在线国产一区| 久久中文看片网| 久久人人爽av亚洲精品天堂| 岛国视频午夜一区免费看| 久久人人爽av亚洲精品天堂| 欧美老熟妇乱子伦牲交| 9色porny在线观看| 亚洲熟妇中文字幕五十中出| 黄片播放在线免费| 深夜精品福利| 精品久久久久久久久久免费视频| 精品久久久久久,| 禁无遮挡网站| 欧美国产精品va在线观看不卡| av在线播放免费不卡| 午夜精品在线福利| 精品欧美一区二区三区在线| 欧美日韩黄片免| 一区二区三区精品91| 欧美激情久久久久久爽电影 | 国产精品一区二区三区四区久久 | 亚洲国产看品久久| 99久久国产精品久久久| 亚洲成人免费电影在线观看| 久久中文字幕一级|