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

    拉薩河中下游原生動物群落物種分布格局及其驅(qū)動因素

    2023-04-29 00:44:03達珍張鵬安瑞志潘成梅喬楠茜巴桑
    水生態(tài)學雜志 2023年4期
    關(guān)鍵詞:原生動物拉薩河物種多樣性

    達珍 張鵬 安瑞志 潘成梅 喬楠茜 巴桑

    摘要:了解拉薩河中下游原生動物群落物種分布格局及其驅(qū)動因素,為拉薩河中下游生物多樣性保護及可持續(xù)利用提供基礎(chǔ)資料。選取了拉薩河中下游17個樣點,于2015年5月、8月和2016年11月進行原生動物定性和定量樣品采集,共采集水樣153份,分析物種組成、生物多樣性及其與環(huán)境因子之間的關(guān)系,確定該流域原生動物類群組成中稀有種、偶見種和常見種對物種多樣性的貢獻、生物多樣性格局及群落與理化因子相關(guān)性。結(jié)果表明:(1)3個季節(jié)共調(diào)查到稀有種80種,占總物種數(shù)的40%;偶見種95種,占48%;常見種21種,占11%;恒見種1種,占0.5%。(2)拉薩河中下游原生動物偶見種是原生動物群落物種多樣性格局的重要貢獻者。(3)多樣性指數(shù)表明,Richness指數(shù)和Margalef指數(shù)春季和夏季高于秋季,Shannon-Wiener指數(shù)和Simpson指數(shù)春季高于秋季高于夏季。(4)基于Bray-Curtis距離的主坐標分析和相似性分析結(jié)果表明,拉薩河中下游原生動物群落結(jié)構(gòu)在不同季節(jié)和不同流域均具有顯著差異,且流域差異相比季節(jié)差異較大。(5)根據(jù)Mantel檢驗分析原生動物與環(huán)境因子的相關(guān)性,在所測定環(huán)境因子中,水溫(WT)、總?cè)芙庑怨腆w(TDS)、鹽度(Salt)、電導(dǎo)率(EC)、總磷(TP)及溶解氧(DO)是影響原生動物群落結(jié)構(gòu)的主要環(huán)境因子。

    關(guān)鍵詞:原生動物;物種多樣性;環(huán)境因子;拉薩河

    中圖分類號:Q178.1 ? ? ? ?文獻標志碼:A ? ? ? ?文章編號:1674-3075(2023)04-0010-08

    原生動物是動物界中一大類單細胞真核生物的統(tǒng)稱,其個體小,生長發(fā)育快,繁殖能力強,在水生態(tài)系統(tǒng)物質(zhì)循環(huán)、能量傳遞與生物資源利用過程中發(fā)揮著十分重要的作用,其物種組成、優(yōu)勢種以及多樣性能直接反映水環(huán)境的狀況(Pawlowski,2014;安瑞志等,2021)。物種分布格局及其驅(qū)動因素一直是生態(tài)學家關(guān)注的焦點(Lennon et al,2004;王思齊和朱章明,2020)。Moura (2021)根據(jù)物種的出現(xiàn)頻率,將物種劃分為恒有種、常見種、偶見種和稀有種,其中有關(guān)稀有種、偶見種和常見種對物種多樣性格局相對貢獻的研究較多,主要集中于植物群落(劉旻霞等,2020)和浮游植物群落(Moura et al,2021;索雨澤等,2021),然而很少有人將此分析方法應(yīng)用于原生動物群落的相關(guān)研究。

    雅魯藏布江位于青藏高原西南部,是我國水能資源最豐富的河流(彭定志等,2008)。雅魯藏布江流域由于高寒環(huán)境、地形復(fù)雜、人口稀疏、交通不便等因素,保持著相對穩(wěn)定的原始生態(tài)系統(tǒng),因而備受國內(nèi)外研究生物起源、進化、區(qū)系和遷徙的學者的廣泛關(guān)注(Zhang et al,2018;安瑞志等,2020)。拉薩河是雅魯藏布江五大支流之一,在調(diào)節(jié)區(qū)域氣候、水資源供給和維持生物多樣性等方面起著重要的生態(tài)安全屏障作用。目前,有關(guān)拉薩河的水生生物學研究主要集中在浮游植物(君珊等,2019)和原生動物(巴桑等,2017;楊欣蘭等,2019)等方面,但針對拉薩河中下游原生動物的系統(tǒng)和持續(xù)研究仍較少。

    為了解拉薩河中下游原生動物群落物種分布格局及其驅(qū)動因素,于2015年5月、8月和2016年11月對拉薩河中下游原生動物進行了定性和定量樣品采集,分析其物種組成、生物多樣性及其與環(huán)境因子之間的關(guān)系,確定該流域原生動物類群組成中稀有種、偶見種和常見種對物種多樣性的貢獻、生物多樣性格局及群落與理化因子相關(guān)性,以期為拉薩河中下游生物多樣性保護及可持續(xù)利用提供基礎(chǔ)資料。

    1 ? 研究方法

    1.1 ? 區(qū)域概況

    拉薩河位于青藏高原的中南部,是雅魯藏布江長度最長、面積最大的支流,發(fā)源于念青唐古拉山南麓(安寶晟等,2021)。根據(jù)河谷形態(tài)和地質(zhì)情況,分為上游、中游、下游3段(何柳,2019)。自源頭至桑曲縣匯入口為上游段,位于自然區(qū)域,全長256 km,氣候常年寒冷干燥,空氣稀薄,人口稀少,多為牧民,生態(tài)環(huán)境較為原始(劉久潭等,2020)。從桑曲到直孔為中游段,全長138 km,為高原溫帶半干旱季風氣候,日照充足,晝夜溫差大,人類活動較少(陳歆等,2019)。直孔以下為下游段,全長157 km,地勢平緩,河谷寬廣,相比中上游,人類活動干擾較強,沿拉薩市流域城鎮(zhèn)數(shù)量變多,人口密集,且該流域為拉薩市區(qū)主要糧食產(chǎn)區(qū)(陳麗等,2019)。

    1.2 ? 樣點設(shè)置

    依據(jù)《水域生態(tài)系統(tǒng)觀測規(guī)范》,在拉薩河中下游河段共計設(shè)置17個樣點(圖1),樣點基本信息見表1。分別于2015年5月、8月和2016年11月進行采樣,共采集樣本153份。

    1.3 ? 水體理化因子測定

    水體pH、水溫(WT)、總?cè)芙庑怨腆w(TDS)、鹽度(Salt)及電導(dǎo)率(EC)現(xiàn)場用PSCTestr35防水型多參數(shù)筆測試儀(EUTECH,美國)測定;溶解氧(DO)現(xiàn)場用HI9146便攜式溶解氧測定儀(HANNA,德國)測定;水流速度(WS)用FP111直讀式流速儀(HANNA,美國)測定;濁度(Tur)用HI98703高精度數(shù)據(jù)型濁度測定儀(HANNA,意大利)測定;總氮(TN)、總磷(TP)、重鉻酸鹽指數(shù)(CODCr)、氨氮(NH3-N)和硝態(tài)氮(NO3-N)送檢測公司測定,測定方法參照國家環(huán)境保護總局《水和廢水監(jiān)測分析方法》和《地表水環(huán)境質(zhì)量標準》(GB3838-2002)送檢測公司進行測定。

    1.4 ? 樣品采集與處理

    樣品采集標準及實驗室處理方法均參照《淡水浮游生物研究方法》(章宗涉和黃祥飛,1991)進行。定性樣品用25號浮游生物網(wǎng)(孔徑64 μm)在水面下0.8~1.2 m處作“∞”字型撈取5~10 min,將濃縮樣品裝入標本瓶,在室內(nèi)利用Olympus CX23研究型光學顯微鏡,采用活體鏡檢和固定染色法鑒定物種(中國科學院青藏高原綜合科學考察隊,1981;沈韞芬,1999)。定量樣品在水面下0.8~1.2 m采集原水1 L,加10 mL(10%)魯哥氏液固定,實驗室靜置24 h沉淀濃縮至50 mL。移取0.1 mL濃縮液至0.1 mL計數(shù)框,在10×40倍顯微鏡下觀察,全片平行計數(shù)3次(計數(shù)差異不超過15%)。

    1.5 ? 數(shù)據(jù)統(tǒng)計與分析

    參照Moura(2005)計算原生動物在群落中的出現(xiàn)頻率(FO),分為恒見種(FO≥70%)、常見種(30%≤FO<70%)、偶見種(10%≤FO<30%)和稀有種(FO<10%)。

    利用Excel 2016軟件進行原始數(shù)據(jù)的基礎(chǔ)統(tǒng)計處理;利用R語言軟件Vegan包計算Shannon-Wiener指數(shù)(H')、Margalef指數(shù)(dM)、Richness指數(shù)(R)、Simpson指數(shù)(D)、單因素方差分析(one-way ANOVA);Tidyverse、Vegan包進行Mantel檢驗分析;ggplot2、ggcor包進行繪圖;利用Origin pro 2019分析原生動物稀有種、偶見種和常見種對群落物種多樣性的相對貢獻;利用ArcGIS10.7軟件繪制研究區(qū)域地圖。

    2 ? 結(jié)果與分析

    2.1 ? 水環(huán)境因子

    拉薩河中下游春季、夏季、秋季的水環(huán)境因子如表2所示。

    由表2可得,pH、總氮、電導(dǎo)率、總?cè)芙庑怨腆w和氨氮的值秋季>春季>夏季;水流速度、總磷的值夏季>春季>秋季;濁度、溶解氧的值夏季>秋季>春季;鹽度值春季>秋季>夏季;硝態(tài)氮值春季>夏季>秋季。

    2.2 ? 原生動物群落物種組成

    拉薩河中下游3個季節(jié)共鑒定出原生動物197種,其中,春、夏、秋分別鑒定出97、106、107種。各季節(jié)稀有種、偶見種、常見種數(shù)量及占比見表3。

    2.3 ? 原生動物群落多樣性

    春、夏、秋3個季節(jié)多樣性指數(shù)結(jié)果見表4。其中,Richness指數(shù)(R)和Margalef指數(shù)(dM)春季[≥]夏季>秋季,Shannon-Wiener指數(shù)(H')和Simpson指數(shù)(D)春秋季節(jié)>夏季。

    拉薩河中下游各樣點不同季節(jié)多樣性指數(shù)如圖2所示,3個季節(jié)各樣點多樣性指數(shù)比較發(fā)現(xiàn),整體上,中游段各樣點(L1、L5、L6)指數(shù)高于下游段各樣點(L14、L9、L12、L13)。

    2.4 ? 出現(xiàn)頻率對物種多樣性的相對貢獻

    原生動物稀有種、偶見種和常見種對物種多樣性的相對貢獻如圖3所示。從各樣點看,春季去除偶見種引起的Richness指數(shù)、Shannon多樣性指數(shù)和Simpson指數(shù)減小效應(yīng)最強,其次為去除稀有種和常見種。夏季去除偶見種引起的Richness指數(shù)、Shannon多樣性指數(shù)和Simpson指數(shù)減小效應(yīng)最強,其次為去除稀有種,去除常見種減小效應(yīng)最弱。秋季去除偶見種引起的Richness指數(shù)、Shannon多樣性指數(shù)和Simpson指數(shù)減小效應(yīng)最強,其次為去除稀有種,而去除常見種引起減小效應(yīng)并不明顯。

    比較拉薩河中下游不同季節(jié)原生動物去除稀有種、偶見種和常見種對物種多樣性的減小效應(yīng),當去除偶見種時,各樣點物種多樣性減小效應(yīng)表現(xiàn)最強,表明在原生動物物種類群組成中偶見種是物種多樣性格局的重要貢獻者。

    2.5 ? 原生動物群落相似性分析

    基于Bray-Curtis距離的主坐標分析(Principal Coordinate Analysis,PCoA)和相似性分析(Analysis Of Similarities,ANOSIM)結(jié)果如圖4a和圖4b,拉薩河中下游原生動物群落在不同季節(jié)(R=0.18,P=0.001)和不同流域(R=0.42,P=0.001)均具有顯著性差異(P<0.05),且流域間差異(R=0.42)比季節(jié)間差異(R=0.18)更大。為探究拉薩河中下游原生動物群落時空格局的形成機制,對原生動物群落相似性與環(huán)境因子及地理距離的相關(guān)性進行分析,結(jié)果如圖4c和圖4d,拉薩河中下游原生動物群落相似性隨環(huán)境因子和地理距離的差異增大而具有顯著衰減趨勢(P<0.01),這表明拉薩河中下游原生動物群落差異性是由環(huán)境因子差異和地理距離差異共同作用的結(jié)果。

    2.6 ? 群落結(jié)構(gòu)與環(huán)境因子的關(guān)系

    利用Mantel檢驗分析拉薩河中下游原生動物群落結(jié)構(gòu)差異與環(huán)境因子差異之間的相關(guān)性,如圖5。在所測定的13個環(huán)境因子中,DO、EC、Salt、TDS、WT和TP是驅(qū)動拉薩河中下游原生動物群落組成變化的最顯著環(huán)境因子(P<0.01)。

    3 ? 討論

    3.1 ? 原生動物群落結(jié)構(gòu)特征

    生物多樣性反映生物群落組成的復(fù)雜性和穩(wěn)定性,是衡量水生態(tài)系統(tǒng)環(huán)境狀態(tài)的重要指標(Pool et al,2014;陳澤豪等,2021)。從空間尺度上,拉薩河中下游原生動物物種多樣性分析結(jié)果顯示,中游段多樣性指數(shù)高于下游段,表明中游原生動物群落組成較下游復(fù)雜,群落穩(wěn)定性較下游好,這與方艷紅等(2014)對大渡河河口春秋季原生動物的群落結(jié)構(gòu)研究結(jié)果基本一致。這可能是由于中游地區(qū)主要為農(nóng)牧區(qū),生態(tài)環(huán)境較為原始,水環(huán)境受人為干擾影響較少,水質(zhì)較好;而下游地區(qū)主要為拉薩市城鎮(zhèn)區(qū),人口密度大,受人類活動干擾較為嚴重,且河流易受工農(nóng)業(yè)廢水與生活污水等污染,對水生態(tài)系統(tǒng)造成一定影響,導(dǎo)致水質(zhì)較差,多樣性指數(shù)低。研究結(jié)果表明,原生動物物種多樣性在人類活動干擾下,表現(xiàn)出較強的空間異質(zhì)性。如魏洪祥等(2018)對臥龍湖浮游生物群落影響因子的研究發(fā)現(xiàn),臥龍湖濕地生態(tài)系統(tǒng)的穩(wěn)定性受不同程度人類開發(fā)活動和自然環(huán)境因素的影響,采樣面積逐年下降,生物種群急劇減少,致使不同區(qū)段采集鑒定的物種種類組成上呈現(xiàn)明顯的空間異質(zhì)性,物種多樣性指數(shù)上體現(xiàn)空間分布差異。

    從時間尺度上,拉薩河中下游原生動物群落Richness指數(shù)、Margalef指數(shù)春季和夏季高于秋季,可能是由于春季(12.9℃)和夏季(12.2℃)的水溫較秋季(5.4℃)更適合浮游植物生長與繁殖,浮游植物在適宜的溫度范圍內(nèi),生長速度加快,可為原生動物生命活動提供更充足的營養(yǎng)物質(zhì)、食物來源及棲息場所,從而使原生動物物種數(shù)目變多、種類變豐富、多樣性指數(shù)升高(金瓊貝等,1991;方艷紅等,2012)。這與代培等(2019)在探究太湖五里湖湖濱帶浮游動物群落結(jié)構(gòu)特征中原生動物生物多樣性隨水溫的季節(jié)變化規(guī)律基本一致。原生動物作為水生態(tài)系統(tǒng)食物鏈中的初級消費者,同時受初級生產(chǎn)者和次級消費者種內(nèi)種間競爭和捕食的影響導(dǎo)致其多樣性指數(shù)降低(Zhu et al,2012;龔森森等,2021)。原生動物群落Shannon-Wiener指數(shù)(H')季節(jié)變化趨勢均為春季和秋季高于夏季,從水生態(tài)系統(tǒng)食物鏈角度來看,可能是由于夏季水溫逐漸升高,促進以原生動物為食的其他大型浮游生物和魚類等生物的生長代謝,從而加快水體中營養(yǎng)物質(zhì)的消耗速度,加劇物種內(nèi)部競爭與捕食關(guān)系(陳旭淼和徐奎棟,2014)。這與安睿(2016)對三江平原典型濕地浮游生物群落特征研究結(jié)果類似。本研究結(jié)果表明,原生動物在應(yīng)對人類不同程度干擾與水環(huán)境變化等情況時,反應(yīng)靈敏迅速,其群落組成與物種多樣性能夠在一定程度上反映水環(huán)境適應(yīng)性變化與水環(huán)境質(zhì)量狀況。

    3.2 ? 物種組成對物種多樣性的相對貢獻

    物種多樣性格局的形成不僅依賴于群落的構(gòu)建過程,還依賴于群落中物種類群組成(Marquet et al,2004;王世雄等,2016)。群落中不同物種類群組成是決定物種多樣性格局的重要因素,研究不同物種類群組成對物種多樣性相對貢獻是探究物種多樣性格局形成過程的一個重要途徑(Lennon et al,2004)。本研究首次通過去除原生動物稀有種(或偶見種、常見種)形成一個“新”的群落,然后比較“新”群落物種多樣性與原群落物種多樣性間減小效應(yīng)來間接度量稀有種、偶見種和常見種對群落物種多樣性格局的相對貢獻。本研究中,當去除偶見種時,原生動物物種多樣性減小效應(yīng)最強,表明偶見種是原生動物物種多樣性格局的重要貢獻者。這與王世雄等(2018)運用此分析方法對黃土高原遼東櫟群落物種多樣性貢獻的多尺度分析結(jié)果存在一定的差異,這可能是由于2個研究研究定義稀有種、偶見種和常見種的標準不同,導(dǎo)致群落物種類群組成的比例差異,因此造成不同物種類群組成對物種多樣性格局的貢獻不一致。也有可能是因為不同季節(jié)物種類群組成在應(yīng)對環(huán)境差異變化時,具有不同的功能性狀,如擴散能力或者環(huán)境耐受力,它們對物種多樣性的主導(dǎo)作用受其自身比例及其他物種比例等因素的影響(彭少麟等,2003)。如稀有種本身對各種因素的干擾表現(xiàn)更為敏感,生存環(huán)境較為局限,因此在種間種內(nèi)競爭過程中比例逐漸減少,去除稀有種引起的整個群落物種多樣性的變化較小。

    3.3 ? 環(huán)境因子對群落結(jié)構(gòu)的影響

    原生動物群落結(jié)構(gòu)特征不僅與食物供給、競爭及捕食—被捕食的關(guān)系等生物因子影響有關(guān),還與pH、水溫、溶解氧等非生物因子密切相關(guān)(龔森森等,2021)。使用Mantel檢驗方法來確定原生動物群落結(jié)構(gòu)差異與環(huán)境變量差異之間的潛在相關(guān)性,結(jié)果表明,拉薩河中下游流域原生動物群落結(jié)構(gòu)均與水溫、總?cè)芙庑怨腆w、鹽度、電導(dǎo)率、總磷、溶解氧呈顯著相關(guān)性。水溫是影響拉薩河中下游流域原生動物群落生長和數(shù)量變化的一個重要因素,水溫會影響水體的理化性質(zhì)(如DO、Salt、pH),直接或間接地影響原生動物的食物來源、生長代謝和其他生命活動(Wang et al,2019)。在溫度適宜的條件下,豐富的浮游植物為原生動物提供大量餌料與棲息場所,使原生動物多樣性大大提高。此外,溶解氧和總磷也是影響該流域原生動物群落結(jié)構(gòu)的重要因素,這與溶解氧和總磷均是生命體進行正常生命活動的必要物質(zhì)密不可分(魏洪祥等,2018)。浮游生物對溶解氧依賴程度隨溫度升高而增加,當水溫升高時,水中的溶解氧會直接限制生物的豐度,水中含氧量隨溫度而有晝夜和季節(jié)的變化,所以,不同原生動物對氧的需求也有一定的適應(yīng)范圍(沈韞芬,1990)。水體鹽度可改變原生動物的細胞滲透壓和水體各種離子,進而影響原生動物的生命活動。Davidson等(2010)對南海原生生物群落的組成和分布展開調(diào)查,發(fā)現(xiàn)原生動物群落與水溫、鹽度及電導(dǎo)率具有顯著相關(guān)性,與本研究調(diào)查結(jié)果一致??傮w而言,非生物因子是決定某種原生動物能否生活在某一特定生境中的決定性條件,而生物因子則是確定某一種群相對豐度的首要條件。本研究與前人對原生動物的群落結(jié)構(gòu)特征研究結(jié)果類似,表明造成其群落物種豐度、組成、結(jié)構(gòu)等方面的差異是原生動物不同類群生活史特征與生境條件協(xié)同作用的結(jié)果(趙瑞等,2019)。

    參考文獻

    安瑞志,張鵬,達珍,等,2021. 西藏麥地卡濕地不同水文期原生動物優(yōu)勢種生態(tài)位及其種間聯(lián)結(jié)性[J].林業(yè)科學,57(2):126-138.

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

    安寶晟,姚檀棟,郭燕紅,等,2021. 拉薩河流域典型區(qū)域保護、修復(fù)、治理技術(shù)示范體系[J].科學通報,66(22):2775-2784.

    安睿,2016. 三江平原典型濕地浮游生物群落特征及其與環(huán)境因子的關(guān)系[D].哈爾濱:東北林業(yè)大學.

    巴桑,楊欣蘭,黃香,等,2017. 拉薩河中上游春季原生動物群落特征[J].西北師范大學學報(自然科學版),53(5):70-76,83.

    陳旭淼,徐奎棟,2014. 深海原生動物多樣性研究進展[J].海洋科學,38(10):119-126.

    陳歆,靳甜甜,蘇輝東,等,2019. 拉薩河河流健康評價指標體系構(gòu)建及應(yīng)用[J].生態(tài)學報,39(3):799-809.

    陳麗,王東波,君珊,等,2019. 拉薩河流域大型底棲動物群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系[J].生態(tài)學報,39(3):757-769.

    陳澤豪,楊文,王穎,等,2021. 白洋淀大型底棲動物群落結(jié)構(gòu)及其與環(huán)境因子關(guān)系[J].生態(tài)學雜志,40(7):2175-2185.

    代培,劉凱,周彥鋒,等,2019. 太湖五里湖湖濱帶浮游動物群落結(jié)構(gòu)特征[J].水生態(tài)學雜志,40(1):55-63.

    方艷紅,黃道明,王文君,等,2012. 大渡河河口秋季浮游動物的群落結(jié)構(gòu)特征[J].水生態(tài)學雜志,33(3):34-40.

    方艷紅,王崇,王文君,等,2014. 大渡河河口春秋季原生動物的群落結(jié)構(gòu)[J].水生態(tài)學雜志,35(4):29-34.

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

    龔森森,吳嘉偉,柴毅,等,2021. 長湖浮游動物群落結(jié)構(gòu)特征及其季節(jié)變化[J].水產(chǎn)科學,40(3):329-338.

    何柳,2019. 拉薩河流域水文地球化學特征及其風化指示[D].南昌:東華理工大學.

    君珊,王東波,周健華,等,2019. 拉薩河流域浮游植物群落結(jié)構(gòu)特征及與環(huán)境因子的關(guān)系[J].生態(tài)學報,39(3):787-798.

    劉久潭,李穎智,高宗軍,等,2020. 拉薩河流域中下游地區(qū)水化學及地表水-地下水轉(zhuǎn)化關(guān)系研究[J].山東科技大學學報(自然科學版),39(5):10-20.

    劉旻霞,李全弟,蔣曉軒,等,2020. 甘南亞高寒草甸稀有種對物種多樣性和物種多度分布格局的貢獻[J].生物多樣性,28(2):107-116.

    彭定志,徐宗學,鞏同梁,2008. 雅魯藏布江拉薩河流域水文模型應(yīng)用研究[J].北京師范大學學報(自然科學版),(1):92-95.

    彭少麟,殷祚云,任海,等,2003.多物種集合的種-多度關(guān)系模型研究進展[J].生態(tài)學報,(8):1590-1605.

    索雨澤,劉洋,劉宇康,等,2021. 雅魯藏布江中上游浮游植物污染指示種及水質(zhì)評價[J].環(huán)境生態(tài)學,3(3):26-32.

    沈韞芬,1990. 微型生物監(jiān)測新技術(shù)[M].北京:中國建筑工業(yè)出版社.

    沈韞芬,1999.原生動物學[M].北京:科學出版社.

    王思齊,朱章明,2022. 中國薔薇屬植物物種豐富度分布格局及其與環(huán)境因子的關(guān)系[J].生態(tài)學報,42(1):209-219.

    王世雄,趙亮,李娜,等,2016. 稀有種和常見種對植物群落物種豐富度格局的相對貢獻[J].生物多樣性,24(6):658-664.

    王世雄,王孝安,郭華,等,2018. 稀有種和常見種對黃土高原遼東櫟群落物種多樣性貢獻的多尺度分析[J].生態(tài)學報,38(22):8060-8069.

    魏洪祥,石俊艷,王曉光,等,2018. 臥龍湖浮游生物群落的演替、影響因子及其水質(zhì)評價[J].生態(tài)科學,37(1):55-61.

    楊欣蘭,巴桑,黃香,等,2019. 拉薩河中上游夏秋季纖毛蟲群落時空變動及其與環(huán)境的關(guān)系[J].生態(tài)學報,39(9):3121-3132.

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

    趙瑞,趙學敏,茍婷,等,2019.廣西貝江大型底棲動物群落結(jié)構(gòu)時空分布特征[J].中國環(huán)境監(jiān)測,35(4):40-49.

    中國科學院青藏高原綜合科學考察隊,1981.西藏地熱[M].北京:科學出版社.

    Lennon J J, Koleff P, Greenwood J J D, et al,2004. Contribution of rarity and commonness to patterns of species richness[J]. Ecology Letters, 7(2) :81-87.

    Marquet P A, Fernández M, Navarrete S A, et al,2004. Diversity emerging toward a deconstruction of biodiversity patterns[M]//Lomolino M,Heaneyl.Frontiers of Biogeography :New dinections in the Geography of Nature,Sunderland MA:Sinauer Associates:191-209.

    Moura L C S,dos Santos S M, de Souza C A, et al,2021. Phytoplankton richness and abundance in response to seasonality and spatiality in a tropical reservoir[J].Acta Limnologica Breasiliensia,33(1):DOI:10.1590/s2179-975x11419.

    Pool T K, Grenouille G, Villéger S, 2014.Species contribute differently to the taxonomic,functional,and phylogenetic alpha and beta diversity of freshwater fish communities[J]. Diversity and Distributions, 20(11) : 1235-1244.

    Zhang Y J, Duo L, Pang Y Z, et al, 2018.Modern pollen assemblages and their relationships to vegetation and climate in the Lhasa Valley, Tibetan Plateau, China[J]. Quaternary International, 467:210-221.

    Zhu W T, Zhang A K, Qin C X, et al, 2021.Seasonal and spatial variation of protist communities from reef water and open ocean water in patchy coral reef areas of a semi-enclosed bay[J]. Marine Environmental Research Volume,169:105407.

    (責任編輯 ? 鄭金秀)

    Species Distribution and Factors Driving the Protozoan Community

    in the Middle and Lower Lhasa River

    DA Zhen, ZHANG Peng, AN Rui‐zhi, PAN Cheng‐mei, QIAO Nan‐qian, BA Sang

    ( Laboratory of Wetland and Catchments Ecology in Tibetan Plateau, Faculty of Natural Sciences,

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

    Abstract:Lhasa River is one of the five largest tributaries of the Yarlung Zangbo River and plays important roles as an ecological corridor, in regional climate regulation, as a water source and in maintaining biodiversity in Tibet. The current study explored species distribution patterns and factors influencing the protozoan community in the middle and lower reaches of Lhasa River, aiming to provide basic data for conserving biodiversity and water resources. In May (spring), August (summer) of 2015 and November (autumn) of 2016, qualitative and quantitative protozoan samples were collected at 17 sampling sites in the middle and lower reaches of the Lhasa River. Species composition and distribution patterns of the protozoan community and their correlation with physicochemical factors of the river water were then analyzed. Results show that: (1) A total of 197 protozoan species were identified during the three seasons, including 80 rare species, 95 occasional species, 21 common species, and 1 constant species, which accounted for 40%, 48%, 11% and 0.5% of the total species, respectively. (2) Occasional species were important contributors to the species diversity of the protozoan community in Lhasa River. (3) The Richness and Margalef indices were higher in spring and summer than in autumn, while the Shannon-Wiener and Simpson indices were higher in spring and autumn than in summer. (4) The results of principal coordinate analysis and similarity analysis based on Bray-Curtis distance showed that the protozoan community structure in the middle and lower reaches of Lhasa River varied significantly with season and river segment, with larger variations between river segments than between seasons. (5) The Mantel test was used to analyze the correlations between protozoan community parameters and environmental factors. Water temperature (WT), total dissolved solids (TDS), salinity (Salt), electrical conductivity (EC), total phosphorus (TP) and dissolved oxygen (DO) were the primary environmental factors affecting protozoan community structure.

    Key words: protozoan; species diversity; environmental factors; Lhasa River

    猜你喜歡
    原生動物拉薩河物種多樣性
    青海湖流域原生動物群落結(jié)構(gòu)及影響因素研究
    拉薩河流域浮游植物生態(tài)學研究現(xiàn)狀
    DNA條形碼在原生動物系統(tǒng)分類學中的應(yīng)用
    動物學實驗教學中“原生動物的結(jié)構(gòu)觀察”伸縮泡觀察方法的優(yōu)化
    拉薩河,從古遠流來
    布達拉(2018年4期)2018-05-14 08:53:30
    拉薩市拉薩河城區(qū)段綜合整治工程4~#閘設(shè)計
    江西大崗山木荷天然林群落物種多樣性分析
    綠色科技(2016年21期)2016-12-27 10:53:26
    納板河保護區(qū)不同林分群落喬木物種多樣性研究
    用固定樣地法監(jiān)測銀瓶山闊葉林物種多樣性
    基于概率統(tǒng)計和模糊綜合評價法的水質(zhì)評價模型及其應(yīng)用——以拉薩河水質(zhì)評價為例
    757午夜福利合集在线观看| 夜夜夜夜夜久久久久| 高潮久久久久久久久久久不卡| 看十八女毛片水多多多| 黄色女人牲交| 日日干狠狠操夜夜爽| 十八禁人妻一区二区| 久久草成人影院| 亚洲人成网站高清观看| 亚洲国产高清在线一区二区三| 久久亚洲精品不卡| 国产精品久久视频播放| 国产久久久一区二区三区| 窝窝影院91人妻| 不卡一级毛片| 国产一区二区亚洲精品在线观看| 国内精品久久久久久久电影| 丁香六月欧美| 在线观看av片永久免费下载| 成人特级av手机在线观看| 亚洲国产色片| 亚洲中文日韩欧美视频| 欧美日本亚洲视频在线播放| 亚洲国产欧美人成| 成年人黄色毛片网站| 国内久久婷婷六月综合欲色啪| 亚洲国产日韩欧美精品在线观看| 亚洲美女黄片视频| av视频在线观看入口| 国产成人av教育| 99久久精品一区二区三区| 12—13女人毛片做爰片一| 亚洲av成人av| 国产精品嫩草影院av在线观看 | 丁香六月欧美| av中文乱码字幕在线| 亚洲av成人av| 日韩有码中文字幕| 免费av观看视频| 日本免费一区二区三区高清不卡| 此物有八面人人有两片| 99国产极品粉嫩在线观看| 天堂√8在线中文| 99riav亚洲国产免费| 一卡2卡三卡四卡精品乱码亚洲| 18+在线观看网站| 久久热精品热| 最近视频中文字幕2019在线8| 亚洲国产色片| 91狼人影院| 日本与韩国留学比较| 露出奶头的视频| 国内精品久久久久久久电影| 日日夜夜操网爽| 国产高清视频在线播放一区| 国产精品三级大全| 在线看三级毛片| 小蜜桃在线观看免费完整版高清| 嫩草影视91久久| 90打野战视频偷拍视频| 国产v大片淫在线免费观看| 亚洲午夜理论影院| 少妇丰满av| 日日干狠狠操夜夜爽| 国产精品亚洲av一区麻豆| 在线观看舔阴道视频| 日本与韩国留学比较| 亚洲av一区综合| 99热这里只有是精品在线观看 | 久久6这里有精品| 亚洲av不卡在线观看| 免费看光身美女| 少妇被粗大猛烈的视频| 国产精华一区二区三区| 欧美最黄视频在线播放免费| 亚洲精品亚洲一区二区| 欧美一区二区亚洲| 欧洲精品卡2卡3卡4卡5卡区| 国产精品98久久久久久宅男小说| 亚洲熟妇中文字幕五十中出| 热99re8久久精品国产| 赤兔流量卡办理| 精品人妻1区二区| 一级毛片久久久久久久久女| 欧美高清性xxxxhd video| 成年版毛片免费区| 日韩欧美一区二区三区在线观看| 亚洲精品久久国产高清桃花| 国产精品一区二区免费欧美| 好男人电影高清在线观看| 噜噜噜噜噜久久久久久91| 国产69精品久久久久777片| 一本综合久久免费| 午夜福利欧美成人| 国产精品嫩草影院av在线观看 | 国产成+人综合+亚洲专区| 日韩欧美精品v在线| 噜噜噜噜噜久久久久久91| 成年免费大片在线观看| 人妻制服诱惑在线中文字幕| 热99在线观看视频| 在线观看午夜福利视频| 高潮久久久久久久久久久不卡| 成人毛片a级毛片在线播放| 亚洲无线在线观看| 国产成人欧美在线观看| 国产探花极品一区二区| 91麻豆精品激情在线观看国产| 免费人成视频x8x8入口观看| 又爽又黄a免费视频| 中文亚洲av片在线观看爽| 国产精品久久久久久精品电影| 午夜亚洲福利在线播放| 日本一二三区视频观看| 99riav亚洲国产免费| 制服丝袜大香蕉在线| 国内毛片毛片毛片毛片毛片| 99国产精品一区二区蜜桃av| 久久天躁狠狠躁夜夜2o2o| 美女 人体艺术 gogo| 免费黄网站久久成人精品 | 五月玫瑰六月丁香| 制服丝袜大香蕉在线| 能在线免费观看的黄片| 欧美3d第一页| 国产成人啪精品午夜网站| 好看av亚洲va欧美ⅴa在| 日韩国内少妇激情av| 亚洲专区国产一区二区| 亚洲美女视频黄频| 嫩草影院精品99| 一区二区三区免费毛片| 日本在线视频免费播放| 美女被艹到高潮喷水动态| 97碰自拍视频| 床上黄色一级片| 亚洲欧美清纯卡通| 尤物成人国产欧美一区二区三区| 少妇熟女aⅴ在线视频| 久久精品91蜜桃| 男人舔奶头视频| 天堂动漫精品| 欧美乱妇无乱码| 精品久久久久久久久久免费视频| 亚洲无线观看免费| 成人毛片a级毛片在线播放| 欧美精品国产亚洲| av女优亚洲男人天堂| 黄色丝袜av网址大全| 免费看光身美女| 久久欧美精品欧美久久欧美| 精品国产三级普通话版| 夜夜爽天天搞| 欧美高清性xxxxhd video| 成人高潮视频无遮挡免费网站| 日韩欧美国产一区二区入口| 女生性感内裤真人,穿戴方法视频| 桃红色精品国产亚洲av| 精品人妻一区二区三区麻豆 | 一a级毛片在线观看| 免费av毛片视频| 极品教师在线免费播放| 国产成人aa在线观看| 国产精品野战在线观看| 久久精品91蜜桃| 亚洲国产欧洲综合997久久,| 别揉我奶头~嗯~啊~动态视频| 可以在线观看毛片的网站| 亚洲精品在线美女| 欧美一区二区国产精品久久精品| 国产黄色小视频在线观看| 久久久久国内视频| 99热这里只有是精品在线观看 | 国产精品一及| 亚洲av第一区精品v没综合| 欧美xxxx性猛交bbbb| av在线老鸭窝| 夜夜看夜夜爽夜夜摸| www日本黄色视频网| 99在线人妻在线中文字幕| 超碰av人人做人人爽久久| 一级作爱视频免费观看| 中文字幕人妻熟人妻熟丝袜美| 伦理电影大哥的女人| 国产精品av视频在线免费观看| 欧美一区二区亚洲| 99热这里只有是精品50| eeuss影院久久| 1024手机看黄色片| 国产探花极品一区二区| 女生性感内裤真人,穿戴方法视频| 成年免费大片在线观看| 成人特级黄色片久久久久久久| 色哟哟·www| 欧美成人免费av一区二区三区| av天堂中文字幕网| 免费观看人在逋| 最近最新中文字幕大全电影3| 中文字幕av成人在线电影| 两个人视频免费观看高清| 乱码一卡2卡4卡精品| 国产精品久久久久久亚洲av鲁大| 成人无遮挡网站| 国产激情偷乱视频一区二区| 午夜a级毛片| 一本一本综合久久| 国产探花极品一区二区| 欧美性猛交黑人性爽| 国产成人欧美在线观看| 久久精品国产亚洲av涩爱 | 中文字幕久久专区| 国产极品精品免费视频能看的| 日韩欧美一区二区三区在线观看| x7x7x7水蜜桃| 老司机福利观看| 搞女人的毛片| 久久精品夜夜夜夜夜久久蜜豆| 男人的好看免费观看在线视频| 99久久精品一区二区三区| 午夜a级毛片| 国产精品嫩草影院av在线观看 | 亚洲成人免费电影在线观看| 99国产精品一区二区蜜桃av| 中文字幕高清在线视频| 国产精品av视频在线免费观看| 欧美xxxx性猛交bbbb| 国产精品98久久久久久宅男小说| a级毛片a级免费在线| 亚洲久久久久久中文字幕| 亚洲成人久久性| 嫩草影院新地址| 免费观看精品视频网站| 久久精品夜夜夜夜夜久久蜜豆| 国产精品美女特级片免费视频播放器| 在线观看免费视频日本深夜| 成人午夜高清在线视频| 中文字幕高清在线视频| 欧美精品国产亚洲| 亚洲天堂国产精品一区在线| 人人妻,人人澡人人爽秒播| 俺也久久电影网| 国产精品不卡视频一区二区 | 波多野结衣巨乳人妻| 一a级毛片在线观看| 亚洲最大成人手机在线| 欧美激情国产日韩精品一区| 露出奶头的视频| 久久国产精品人妻蜜桃| 亚洲国产精品sss在线观看| 国产亚洲av嫩草精品影院| 日韩欧美在线乱码| 国产成年人精品一区二区| 99久久精品国产亚洲精品| 欧美区成人在线视频| 亚洲av电影在线进入| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一区二区三区免费毛片| 少妇熟女aⅴ在线视频| 亚洲av成人av| 亚洲男人的天堂狠狠| 人妻夜夜爽99麻豆av| 99久久无色码亚洲精品果冻| 在线观看舔阴道视频| 日韩欧美三级三区| 久久伊人香网站| 大型黄色视频在线免费观看| 真人一进一出gif抽搐免费| 国产伦一二天堂av在线观看| 日本精品一区二区三区蜜桃| 高清在线国产一区| 搡女人真爽免费视频火全软件 | 欧美性猛交╳xxx乱大交人| 中文在线观看免费www的网站| 午夜影院日韩av| 欧美乱妇无乱码| 亚洲色图av天堂| 国产一区二区三区在线臀色熟女| 中文字幕久久专区| 国产精品一及| 尤物成人国产欧美一区二区三区| 热99在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| 国内少妇人妻偷人精品xxx网站| 亚洲五月天丁香| 国产亚洲精品综合一区在线观看| 琪琪午夜伦伦电影理论片6080| 日本在线视频免费播放| 美女免费视频网站| xxxwww97欧美| 成人av一区二区三区在线看| 日韩亚洲欧美综合| 黄色日韩在线| 亚洲精华国产精华精| 内射极品少妇av片p| 欧美性猛交黑人性爽| 中文字幕人妻熟人妻熟丝袜美| 一级a爱片免费观看的视频| 精品熟女少妇八av免费久了| 免费在线观看影片大全网站| 欧美色视频一区免费| 少妇熟女aⅴ在线视频| 国产亚洲欧美在线一区二区| av专区在线播放| 美女高潮喷水抽搐中文字幕| 国产精品一区二区免费欧美| 99久久精品热视频| 直男gayav资源| 亚洲成人精品中文字幕电影| 国产成人av教育| 三级国产精品欧美在线观看| 免费观看的影片在线观看| 国产精品98久久久久久宅男小说| 男女视频在线观看网站免费| 精品福利观看| 69av精品久久久久久| 免费av观看视频| 精品久久国产蜜桃| a级一级毛片免费在线观看| 日本精品一区二区三区蜜桃| 国产aⅴ精品一区二区三区波| 欧美bdsm另类| 亚洲久久久久久中文字幕| 91av网一区二区| 色视频www国产| 日韩欧美国产一区二区入口| x7x7x7水蜜桃| 国产精品一区二区性色av| 俺也久久电影网| 欧美+日韩+精品| 麻豆国产97在线/欧美| 国产成人a区在线观看| 免费av毛片视频| 成人鲁丝片一二三区免费| 波多野结衣高清无吗| 欧美成狂野欧美在线观看| 国产老妇女一区| 欧美成人a在线观看| 欧美激情久久久久久爽电影| 18禁在线播放成人免费| 日韩欧美在线乱码| 舔av片在线| 小蜜桃在线观看免费完整版高清| 成人特级av手机在线观看| 欧美最黄视频在线播放免费| 草草在线视频免费看| 国产日本99.免费观看| 三级国产精品欧美在线观看| 国产淫片久久久久久久久 | 中亚洲国语对白在线视频| 国产精品国产高清国产av| 亚洲av电影在线进入| 日日摸夜夜添夜夜添小说| АⅤ资源中文在线天堂| 一本精品99久久精品77| 18禁黄网站禁片免费观看直播| 美女cb高潮喷水在线观看| 欧美日本视频| 精品一区二区三区人妻视频| 精品一区二区三区av网在线观看| 天堂av国产一区二区熟女人妻| 精品一区二区三区人妻视频| 亚洲精品色激情综合| 国产亚洲精品久久久com| 精品久久久久久久人妻蜜臀av| 亚洲真实伦在线观看| 一级毛片久久久久久久久女| 精品久久久久久久末码| 精品久久久久久久久久免费视频| 日韩欧美国产一区二区入口| 亚洲av第一区精品v没综合| 亚洲欧美激情综合另类| 亚洲av中文字字幕乱码综合| 亚洲精品亚洲一区二区| 国产精品爽爽va在线观看网站| 国产aⅴ精品一区二区三区波| 村上凉子中文字幕在线| 亚洲精品久久国产高清桃花| 国产亚洲精品综合一区在线观看| 久久婷婷人人爽人人干人人爱| 在现免费观看毛片| 久久久久久久久中文| 午夜福利欧美成人| 色综合欧美亚洲国产小说| 少妇被粗大猛烈的视频| 亚洲第一电影网av| 午夜福利欧美成人| 久久久久精品国产欧美久久久| 色综合亚洲欧美另类图片| 精品久久久久久成人av| 久久精品夜夜夜夜夜久久蜜豆| 欧美高清性xxxxhd video| 久久这里只有精品中国| 国产精品久久电影中文字幕| 国产一区二区在线观看日韩| 黄色丝袜av网址大全| 国产精品一及| 高潮久久久久久久久久久不卡| 国产av不卡久久| 亚洲人成网站高清观看| 精品午夜福利在线看| 亚洲三级黄色毛片| 亚洲国产欧美人成| 国产伦一二天堂av在线观看| 亚洲av日韩精品久久久久久密| 午夜日韩欧美国产| 国产精品精品国产色婷婷| 色综合婷婷激情| 一个人免费在线观看的高清视频| 欧美日韩国产亚洲二区| 女人十人毛片免费观看3o分钟| 欧美激情久久久久久爽电影| 老鸭窝网址在线观看| 在线观看免费视频日本深夜| 日韩亚洲欧美综合| 男人舔奶头视频| 热99re8久久精品国产| 亚洲精品在线观看二区| 免费看日本二区| 亚洲国产精品999在线| 色视频www国产| 五月伊人婷婷丁香| 欧美国产日韩亚洲一区| 国产精品免费一区二区三区在线| 国产综合懂色| 成年女人永久免费观看视频| 欧美日本视频| 51午夜福利影视在线观看| 白带黄色成豆腐渣| 午夜老司机福利剧场| 夜夜看夜夜爽夜夜摸| 青草久久国产| 亚洲av成人不卡在线观看播放网| 极品教师在线视频| 免费电影在线观看免费观看| 又黄又爽又刺激的免费视频.| 五月玫瑰六月丁香| 免费av观看视频| 51午夜福利影视在线观看| 亚洲经典国产精华液单 | 成人高潮视频无遮挡免费网站| 国产精品久久久久久人妻精品电影| 99久久九九国产精品国产免费| 成人无遮挡网站| 亚洲精品在线美女| 国内精品久久久久久久电影| 午夜免费成人在线视频| 蜜桃久久精品国产亚洲av| 丰满人妻一区二区三区视频av| 可以在线观看毛片的网站| 欧美性猛交黑人性爽| 欧美日本视频| 日日干狠狠操夜夜爽| 亚洲色图av天堂| 美女cb高潮喷水在线观看| 日韩免费av在线播放| 欧美3d第一页| 国产精品久久久久久精品电影| 日韩精品中文字幕看吧| or卡值多少钱| 国产精品乱码一区二三区的特点| 欧美黑人巨大hd| 国产精品亚洲一级av第二区| 一进一出抽搐动态| 3wmmmm亚洲av在线观看| 日韩欧美国产在线观看| 欧美精品国产亚洲| 好看av亚洲va欧美ⅴa在| 搡老妇女老女人老熟妇| 欧美最新免费一区二区三区 | 在线观看66精品国产| 99在线视频只有这里精品首页| 嫩草影视91久久| 精品午夜福利视频在线观看一区| 国产不卡一卡二| 久久久久久久亚洲中文字幕 | 国产精品国产高清国产av| 啪啪无遮挡十八禁网站| 国产精品一区二区性色av| 99热这里只有是精品在线观看 | 国产成人a区在线观看| 国产精品一区二区三区四区久久| 国产精品人妻久久久久久| 91在线精品国自产拍蜜月| 在线a可以看的网站| 欧美日韩福利视频一区二区| 97人妻精品一区二区三区麻豆| 日本三级黄在线观看| 香蕉av资源在线| 免费在线观看日本一区| 搡老妇女老女人老熟妇| 在线观看66精品国产| 可以在线观看的亚洲视频| 午夜福利在线在线| 极品教师在线视频| 国产高清视频在线观看网站| 国产免费男女视频| 久久天躁狠狠躁夜夜2o2o| 久久99热这里只有精品18| 久久草成人影院| 日韩欧美一区二区三区在线观看| 亚洲av成人精品一区久久| 久久久久久久午夜电影| 亚洲,欧美,日韩| 久久午夜福利片| 大型黄色视频在线免费观看| 国产v大片淫在线免费观看| 精品一区二区免费观看| 久久久久久久久中文| 日日干狠狠操夜夜爽| 国产视频一区二区在线看| 国产中年淑女户外野战色| 美女免费视频网站| 国产大屁股一区二区在线视频| 国产三级黄色录像| 白带黄色成豆腐渣| 真人一进一出gif抽搐免费| av中文乱码字幕在线| 亚洲男人的天堂狠狠| 亚洲18禁久久av| 国产伦在线观看视频一区| 91久久精品国产一区二区成人| 精品一区二区三区人妻视频| 国产成人欧美在线观看| 欧美午夜高清在线| 91九色精品人成在线观看| 欧美成人a在线观看| 日本免费一区二区三区高清不卡| xxxwww97欧美| 热99在线观看视频| 久久国产精品人妻蜜桃| 18禁在线播放成人免费| 国产老妇女一区| 亚洲人成伊人成综合网2020| 亚洲激情在线av| 淫妇啪啪啪对白视频| 一夜夜www| 久久精品91蜜桃| 国产成人aa在线观看| www日本黄色视频网| 国产大屁股一区二区在线视频| 琪琪午夜伦伦电影理论片6080| 色5月婷婷丁香| 亚洲欧美清纯卡通| 深夜a级毛片| 午夜精品一区二区三区免费看| 午夜亚洲福利在线播放| 三级男女做爰猛烈吃奶摸视频| 亚洲成人久久爱视频| 此物有八面人人有两片| 国产亚洲欧美98| 亚洲综合色惰| 一本综合久久免费| 一本一本综合久久| 国产真实伦视频高清在线观看 | av在线蜜桃| 如何舔出高潮| 女生性感内裤真人,穿戴方法视频| 日本熟妇午夜| 久久九九热精品免费| 色在线成人网| 亚洲三级黄色毛片| 一个人观看的视频www高清免费观看| 日本成人三级电影网站| 97热精品久久久久久| 精品国产亚洲在线| 国产精品三级大全| 欧美日本视频| 极品教师在线免费播放| 久久久久免费精品人妻一区二区| 亚洲七黄色美女视频| 国产精品乱码一区二三区的特点| 久久久精品欧美日韩精品| 在线观看舔阴道视频| 国产精品女同一区二区软件 | 变态另类丝袜制服| 99热这里只有是精品50| 观看美女的网站| 麻豆久久精品国产亚洲av| 亚洲成a人片在线一区二区| 精品乱码久久久久久99久播| 欧美+日韩+精品| 久久久国产成人免费| 久久这里只有精品中国| 日日摸夜夜添夜夜添av毛片 | 国产色婷婷99| av视频在线观看入口| 12—13女人毛片做爰片一| 亚洲国产精品久久男人天堂| 成年女人永久免费观看视频| 国产老妇女一区| 国产精品伦人一区二区| 亚洲国产欧美人成| 99视频精品全部免费 在线| 成年女人毛片免费观看观看9| 黄片小视频在线播放| 国产探花极品一区二区| 日韩 亚洲 欧美在线| 国产精品一及| 欧美精品啪啪一区二区三区| 无人区码免费观看不卡| 国产精品久久久久久精品电影| 欧美乱色亚洲激情| 亚洲av免费高清在线观看| 精品国产三级普通话版| 亚洲无线观看免费| 别揉我奶头~嗯~啊~动态视频| 色综合亚洲欧美另类图片| 9191精品国产免费久久| 午夜精品在线福利| av国产免费在线观看| 免费一级毛片在线播放高清视频| 一个人观看的视频www高清免费观看| av在线观看视频网站免费| 99热这里只有是精品50| 毛片一级片免费看久久久久 |