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

    西藏巴松錯夏季浮游植物優(yōu)勢種生態(tài)位及其種間聯(lián)結性

    2024-05-23 19:14:39劉洋安瑞志王陳潘成梅巴桑
    水生態(tài)學雜志 2024年2期
    關鍵詞:浮游植物生態(tài)位

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

    摘要:為了解青藏高原高寒湖泊不同空間維度浮游植物群落優(yōu)勢種的生態(tài)位及種間聯(lián)結性,選取青藏高原東南部最大的高寒冰山堰塞湖之一的巴松錯為研究地點,2020年7月進行了浮游植物定性和定量采集,共采集水平和垂直維度浮游植物水樣275個,結合生態(tài)位寬度、生態(tài)位重疊值、生態(tài)響應速率、總體聯(lián)結性、[χ]2 檢驗、聯(lián)結系數(shù)( AC) 、共同出現(xiàn)百分比( PC) 以及點相關系數(shù)(φ) 等方法,分析優(yōu)勢種的生態(tài)位及種間聯(lián)結性。結果顯示:(1)共鑒定出浮游植物優(yōu)勢種5屬8種,均屬于硅藻門;(2)生態(tài)位寬度值差異較大,可分為3大生態(tài)位類群,水平維度對二維生態(tài)位寬度影響較明顯;(3)在水平維度上,生態(tài)位重疊度達到有意義的種對比例較高,具有明顯的同質性;在垂直和水平-垂直維度上,生態(tài)位重疊度達到有意義的種對比例較低,具有明顯的異質性;(4)夏季浮游植物群落優(yōu)勢種生態(tài)響應速率在水平維度上形成正反饋,在垂直和水平-垂直維度上形成負反饋,表明巴松錯夏季浮游植物群落存在多穩(wěn)態(tài),垂直和水平-垂直維度群落較水平維度更穩(wěn)定;(5)夏季浮游植物群落總體處于演替初期,有正向演替趨勢,物種間聯(lián)結性不強,種間關系較為松散。研究結果為進一步研究高寒湖泊優(yōu)勢種生態(tài)位和群落多穩(wěn)態(tài)提供了基礎資料。

    關鍵詞:浮游植物;優(yōu)勢種;生態(tài)位;種間聯(lián)結;巴松錯

    中圖分類號:Q178.1? ? ? ? 文獻標志碼:A? ? ? ? 文章編號:1674-3075(2024)02-0111-10

    浮游植物是指在水體中營浮游生活的微型藻類,在水生態(tài)系統(tǒng)的食物網(wǎng)和質能流動中起著十分重要的作用(安瑞志等,2020)。因其體積小、結構簡單、生命周期短,對水環(huán)境變化反應敏感,能夠反映水生態(tài)環(huán)境的變化規(guī)律,是水生態(tài)健康的重要指示物種(李港等,2021)。相對于浮游植物其他種群,優(yōu)勢種是研究生態(tài)過程的關鍵(Avolio et al,2019;劉洋等,2022)。

    生態(tài)位反映了一個物種在群落和生態(tài)系統(tǒng)中的地位和狀態(tài),是分析和評價物種間的生存機制、種間關系和群落動態(tài)演替等方面的重要理論(潘成梅等,2021)。種間關系是指不同物種在空間分布上相互作用所形成的關系,是水生態(tài)系統(tǒng)的本質特征之一。種間關聯(lián)性研究有助于分析群落發(fā)育過程中種間相互關系及表征優(yōu)勢種分布,從而研究群落結構、功能及其演替趨勢(趙萍萍等,2012;張苗苗等,2018)。將生態(tài)位理論與種間聯(lián)結性結合,能夠客觀有效地反映自然種群,了解種間資源利用情況和競爭機制,對湖泊生態(tài)系統(tǒng)水生生物多樣性保護、管理和恢復具有重要的指導意義(陶敏等,2021)。迄今將生態(tài)位與種間聯(lián)結性結合的研究主要集中于植物群落(陳龍濤等,2016;江煥等,2019),以浮游植物(吳佳夢等,2019)、原生動物(安瑞志等,2021a)和魚類(董靜瑞等,2017)等水生生物為對象的相關研究近年來才出現(xiàn)。然而,從水平、垂直和水平-垂直等多維度分析浮游植物群落優(yōu)勢種的生態(tài)位及種間聯(lián)結性的研究卻鮮有報道。

    青藏高原高寒湖泊分布星羅棋布,占中國湖泊總面積一半以上,是高原陸-氣系統(tǒng)的基本組成部分(靳錚等,2020)。與其他湖泊相比,高寒湖泊受人類活動干擾較小,對氣候變化的響應更快、更敏感,基本保持原始生態(tài)系統(tǒng),有利于研究自然狀態(tài)下湖泊水生生物群落組成及其對環(huán)境的響應,因而備受國內外學者關注(王捷等,2015;Li et al,2021)。巴松錯位于尼洋河北岸最大支流巴河的高峽深谷中,是青藏高原東南部最大的高寒冰山堰塞湖之一,對維持區(qū)域生態(tài)平衡等起著重要作用。目前,有關巴松錯的水生生物學研究主要集中在浮游植物(安瑞志等,2020;2021b)和原生動物(楊欣蘭等,2022)等方面,但針對其浮游植物的系統(tǒng)研究較少(劉洋等,2023)。

    為了解西藏巴松錯夏季不同維度浮游植物群落優(yōu)勢種的生態(tài)位及種間聯(lián)結性,本研究于2020年7月對巴松錯出入河流及其湖體的浮游植物進行了水平和垂直分層定性和定量采集,共采集水樣275個,從水平、垂直和水平-垂直等多維度對巴松錯夏季浮游植物群落優(yōu)勢種的生態(tài)位和種間聯(lián)結性開展研究,揭示浮游植物群落優(yōu)勢種在不同維度的生態(tài)適應性、功能地位及資源利用狀況,以期為進一步研究其群落的穩(wěn)定性和空間分布格局提供基礎數(shù)據(jù)。

    1? ?材料與方法

    1.1? ?研究設計

    選取巴松錯出入河流及其湖體進行水平和垂直采樣(圖1),其中水平采樣點15個,垂直采樣點11個,共采集水樣275個,采樣時間為2020年7月。具體研究區(qū)概況、樣點設置、樣品采集及鑒定和環(huán)境因子的采樣和分析過程等內容見劉洋等(2023)研究。

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

    1.2.1? ?優(yōu)勢種? ?優(yōu)勢種采用Mcnaughton優(yōu)勢度指數(shù)(Y)確定(Mcnaughton,1967; Habib et al,1997),計算公式如下:

    式中:fi為第i種的出現(xiàn)頻率,ni為第i種的個體數(shù)量,n為同一樣點的個體總數(shù);Y>0.02則為優(yōu)勢種,Y>0.10為絕對優(yōu)勢種。

    1.2.2? ?生態(tài)位? ?采用Colwell等(1971)加權修正的Levins指數(shù)計算生態(tài)位寬度(Bi),計算公式如下:

    式中:Bi為物種i的生態(tài)位寬度,Pij為物種i在第j個樣點的個體數(shù)量占該種個體總數(shù)的比例,Pi=Nij/Ni,i表示物種,j為資源,Nij為種群i利用資源狀態(tài)j的個體數(shù)量,Ni為種群i的個體總數(shù),r為采樣點數(shù)。

    采用Pinaka指數(shù)(Pianka,1973)計算生態(tài)位重疊值(Oik),計算公式如下:

    式中:Oik為生態(tài)位重疊值,表示種i與種k的資源利用曲線的重疊指數(shù),取值區(qū)間為[0,1];Pij和Pkj分別表示種i和種k(k為不同于i的另一物種)在第j個站點中的個體數(shù)量占該種個體總數(shù)的比例,N為樣點總數(shù);R表示生態(tài)響應速率,∑k=1Oik指種群i占用其他種群的總資源量,∑i=1Oik指種群k被其他種群侵占的總資源量,當i=k時,?Oik>0,說明該種處于發(fā)展期;?Oik<0,該種處于衰退期;當?Oik=0,該種是中性的。

    根據(jù)Cody(1974)和May(1975)提出的用于來估計多維資源軸上的生態(tài)位特征值的2個原則,巴松錯夏季浮游植物優(yōu)勢種的水平-垂直二維生態(tài)位寬度和重疊由浮游植物水平和垂直生態(tài)位特征值的乘積來表示:

    水平-垂直生態(tài)位寬度值=水平生態(tài)位寬度×垂直生態(tài)位寬度? ? ? ? ⑥

    水平-垂直生態(tài)位重疊值=水平生態(tài)位重疊×垂直生態(tài)位重疊? ? ? ? ? ⑦

    1.2.3? ?種間聯(lián)結性? ?方差比率法(VR)用于確定總體關聯(lián)性,并使用統(tǒng)計量W值進一步檢驗顯著程度(Schluter,1984),計算公式如下:

    式中:σ2T表示樣點的總體方差,Pi為物種i出現(xiàn)的頻率,ni表示物種i出現(xiàn)的樣點數(shù)量,N表示總樣點數(shù);S2T表示物種的總體方差,S表示優(yōu)勢種的種群總數(shù),Tj表示樣點j內優(yōu)勢種出現(xiàn)的物種數(shù),t表示樣點中優(yōu)勢種的平均物種數(shù)(t=(T1+T2+…+Tn)/N);當VR=1時,物種間沒有關聯(lián),說明物種間是獨立的;當VR>1時,物種間呈現(xiàn)正相關;當VR<1時,表明物種間存在負關聯(lián);統(tǒng)計量W值用于測試VR值是否顯著偏離1;如果物種間總體無顯著關聯(lián),則[χ]20.95(N)

    [χ]2檢驗是基于2×2列聯(lián)表進行種間關聯(lián)的定性研究;由于這項研究是一個不連續(xù)的樣本,因此通過Yates的連續(xù)校正公式對[χ]2進行了校正;當某一物種的出現(xiàn)頻率為100%時,b和d值的權重為1,以避免分母為0時出現(xiàn)不可計算的情況,目的是獲得更客觀的結果 (張金屯, 2004),計算公式如下:

    式中:N是總樣點數(shù),a是包含物種i和k的樣點數(shù),b是包含物種i但不包含k的樣點數(shù),c是包含物種k但不包含i的樣點數(shù),d是物種i和k均不包含的樣點數(shù);當[χ]2<3.841,表示種間無顯著聯(lián)結性(P>0.05);當3.841≤[χ]2≤6.635,表示種間聯(lián)結顯著(0.01≤P≤0.05),若[χ]2>6.635,則表示種對間聯(lián)結性極顯著(P<0.01);因[χ]2計算結果無法反映負值,種對間的正、負聯(lián)結性利用ad-bc的符號進行判斷,若其值為正數(shù),可判定為正聯(lián)結,反之則判定為負聯(lián)結,相等則判定為無聯(lián)結。

    聯(lián)結系數(shù)(AC)用于測定種間聯(lián)結的程度(張金屯, 2004),計算公式如下:

    式中:聯(lián)結系數(shù)AC在[-1,1]范圍內有效,如果AC=0,說明種間完全獨立,AC值越趨向1,表示種間正聯(lián)結程度越高,AC值越趨向-1,則種間的負聯(lián)結性越強。

    共同出現(xiàn)百分率(PC)用于分析種間的正聯(lián)結程度(張金屯, 2004),計算公式如下:

    式中:PC在[0,1]范圍內有效,當PC取值越偏向1,表示種間聯(lián)結程度越強,當種間聯(lián)結性不存在時,PC值等于0。

    點相關系數(shù)(φ) 用于解釋種間聯(lián)結的程度(張金屯, 2004),計算公式如下:

    式中:φ取值在[-1,1]范圍內有效;φ值越偏向1,種間的正聯(lián)結性越強;φ值越偏向-1,種間的負聯(lián)結性越強。

    利用ArcGis對研究區(qū)域經(jīng)緯度數(shù)據(jù)進行可視化標注處理;通過Excel2019軟件計算優(yōu)勢度指數(shù)、生態(tài)位寬度指數(shù)、生態(tài)位重疊指數(shù)和種間聯(lián)結性;使用Origin Pro 2021對優(yōu)勢種水平、垂直和水平-垂直維度的生態(tài)位寬度進行聚類分析及繪圖。

    2? ?結果與分析

    2.1? ?優(yōu)勢種

    巴松錯夏季浮游植物群落共鑒定出優(yōu)勢種5屬8種,均屬于硅藻門。不同維度上優(yōu)勢種的種類和優(yōu)勢度均有差異(表1)。在水平維度上優(yōu)勢種較少,僅有4種,小環(huán)藻屬一種(Cyclotella sp.)和肘狀針桿藻(Synedra ulna)為絕對優(yōu)勢種;在垂直維度上優(yōu)勢種較多,為6種,尖針桿藻(Synedra acus)為絕對優(yōu)勢種。

    2.2? ?生態(tài)位

    2.2.1? ?生態(tài)位寬度? ?巴松錯夏季浮游植物優(yōu)勢種生態(tài)位寬度如表2,在水平維度上,優(yōu)勢種的生態(tài)位寬度為0.067~0.648,平均0.248。在垂直維度上,優(yōu)勢種的生態(tài)位寬度為0.091~0.431,平均0.238;在水平-垂直維度上,生態(tài)位寬度為0.013~0.154,平均0.058。

    根據(jù)浮游植物優(yōu)勢種在不同維度的生態(tài)位寬度聚類分析(圖2),巴松錯夏季浮游植物優(yōu)勢種可分為3大類,第1類為廣生態(tài)位種類,生態(tài)位寬度值較大,包含小環(huán)藻屬一種,這類優(yōu)勢種對水環(huán)境適應能力強,生態(tài)適應幅度廣;第2類為中生態(tài)位種類,生態(tài)位寬度值較小,包含尖針桿藻,這類優(yōu)勢種對生境有一定要求,對水環(huán)境資源有較強選擇性,在適宜條件下能充分利用資源;第3類為窄生態(tài)位種類,生態(tài)位寬度值小,包含纖細等片藻、雙月短縫藻、鈍脆桿藻、脆桿藻屬一種、針桿藻屬一種和肘狀針桿藻,對環(huán)境因子有一定依賴性,對水環(huán)境條件變化較敏感,適應性較弱。

    2.2.2? ?生態(tài)位重疊值? ?巴松錯夏季浮游植物優(yōu)勢種不同維度生態(tài)位重疊值如表3、表4和表5。不同維度物種間對資源利用或生態(tài)適應能力的相似程度存在較大差異,說明不同優(yōu)勢種之間對環(huán)境資源的利用既相互獨立,又相互重疊。

    水平維度上,生態(tài)位重疊值為[0,0.814](表3)。28個種對中,重疊度較高(Oik>0.6)的有9對,占比32.14%;重疊度較低(0.30.3,為有意義的重疊,約占總種對數(shù)的73.91%,表明該湖泊浮游植物優(yōu)勢種在水平維度上具有明顯的同質性。

    垂直維度上,生態(tài)位重疊值為[0,913](表4)。28 個種對中,重疊度較高(Oik>0.6)的有4 對,占比14.29%;重疊度較低(0.30.3,約占總種對數(shù)的17.86%,表明該湖泊浮游植物優(yōu)勢種在垂直維度上達到有意義的種對比例較低,具有明顯的異質性。

    水平-垂直維度上,生態(tài)位重疊值為[0,0.744](表5)。28個種對中,重疊度較高(Oik>0.6)的有1對,占總種對數(shù)的3.57%;重疊度較低(0.30.3,約占22.4%,表明該湖泊浮游植物優(yōu)勢種在水平-垂直維具有明顯的異質性。

    2.2.3? ?生態(tài)響應速率? ?巴松錯夏季浮游植物群落優(yōu)勢種不同維度的相對資源占有量(△Oik)和生態(tài)響應速率(R)如表6。在水平維度上,雙月短縫藻的相對資源占有量(△Oik)變化幅度最大;垂直和水平-垂直維度上,均以尖針桿藻的相對資源占有量(△Oik)變化幅度最大。在水平維度上,小環(huán)藻屬一種、雙月短縫藻和針桿藻屬一種屬于衰退型,其中雙月短縫藻有最大的負△Oik,表明該種衰退空間最大,與其他優(yōu)勢種在資源競爭上處劣勢;纖細等片藻、鈍脆桿藻、脆桿藻屬一種、尖針桿藻和肘狀針桿藻屬于發(fā)展型,其中肘狀針桿藻有最大的正△Oik,表明該種發(fā)展空間最大,與其他優(yōu)勢種在資源競爭中占據(jù)優(yōu)勢;在垂直維度上,纖細等片藻、鈍脆桿藻和尖針桿藻屬于發(fā)展型,其他各優(yōu)勢種處于衰退型,小環(huán)藻屬一種有最大的衰退空間,雙月短縫藻衰退空間最??;在水平-垂直維度上,小環(huán)藻屬一種、雙月短縫藻、脆桿藻屬一種、針桿藻屬一種和肘狀針桿藻屬于衰退型,其他各優(yōu)勢種處于發(fā)展型,尖針桿藻有最大的發(fā)展空間,其次是鈍脆桿藻,纖細等片藻的發(fā)展空間最小。

    基于生態(tài)響應速率(R)可進一步確定優(yōu)勢種在不同維度上對生境條件的生態(tài)響應狀況。在水平維度上,脆桿藻屬一種的發(fā)展趨勢最為強勁(R=0.161),肘狀針桿藻的發(fā)展趨勢最弱(R=0.027);小環(huán)藻屬一種的衰退趨勢最明顯(R=-0.111),針桿藻屬一種的衰退趨勢最弱(R=-0.005)。在垂直維度上,尖針桿藻的發(fā)展趨勢最為強勁(R=0.059),鈍脆桿藻的發(fā)展趨勢最弱(R=0.041);雙月短縫藻的衰退趨勢最為強勁(R=-0.430),小環(huán)藻屬一種的衰退趨勢最弱(R=-0.035);在水平-垂直維度上,尖針桿藻的發(fā)展趨勢最為強勁(R=0.014),纖細等片藻的發(fā)展趨勢最弱(R=0.005);肘狀針桿藻的衰退趨勢最為強勁(R=-0.127),針桿藻屬一種的衰退趨勢最弱(R=-0.128)。

    從總體上看,巴松錯夏季浮游植物優(yōu)勢種在水平維度上生態(tài)響應速率之和為正數(shù)(0.178),在垂直維度上生態(tài)響應速率之和為負數(shù)(-0.73),在水平-垂直維度上生態(tài)響應速率之和為負數(shù)(-0.154),表明巴松錯夏季浮游植物群落優(yōu)勢種在水平維度上形成正反饋,在垂直和水平-垂直維度上形成負反饋調節(jié)。

    2.3? ?種間聯(lián)結性

    2.3.1? ?總體關聯(lián)性? ?巴松錯夏季不同維度浮游植物優(yōu)勢種的種間總體關聯(lián)性如表7,[χ]20.95 (25)=14.611,[χ]20.05 (25)=37.652,VR值均大于1,統(tǒng)計量W為[[χ]20.95(N),[χ]20.05 (N)],表明巴松錯夏季浮游植物優(yōu)勢種總體上呈正關聯(lián)不顯著。

    2.3.2? ?種間聯(lián)結測度? ?巴松錯夏季浮游植物群落優(yōu)勢種的[χ]2檢驗、聯(lián)結系數(shù)(AC)、共同出現(xiàn)百分率(PC)和點相關系數(shù)(φ)結果如圖3。[χ]2檢驗結果顯示,在28個種對內,聯(lián)結性不顯著[χ]2<3.841(P>0.05)的種對有26對,占92.86%;聯(lián)結性顯著(3.841≤[χ]2≤6.635,0.01≤P≤0.05)的種對有1對,占3.57%;聯(lián)結性極顯著([χ]2>6. 635,P<0. 01) 的種對有1對,占 3.57%;聯(lián)結系數(shù)(AC)結果顯示,正聯(lián)結種對有13對,負聯(lián)結種對有14對,無聯(lián)結1對,群落正負聯(lián)結比為0.93,負聯(lián)結占優(yōu)勢。正聯(lián)結對中,有2對正聯(lián)結系數(shù)程度較高(AC≥0.6),約占總對數(shù)的7.14%;5對正聯(lián)結系數(shù)程度一般(0.2≤AC<0.6),約占17.86%;6對正聯(lián)結系數(shù)程度趨向獨立(0

    3? ?討論

    3.1? ?優(yōu)勢種的生態(tài)位寬度和重疊值

    生態(tài)位寬度是衡量物種的分布狀況、利用資源能力均勻程度的綜合指標,反映了物種在特定生境中對環(huán)境資源的利用情況(汪志聰?shù)龋?010)??臻g生態(tài)位包括垂直生態(tài)位和水平生態(tài)位;水平生態(tài)位寬度反映物種在水平空間維度上的分布范圍和均勻程度,垂直生態(tài)位寬度反映物種與資源在垂直空間上的分化和配置方式(肖紅等,2003)。然而,水平或垂直生態(tài)位研究僅反映了物種對一維資源的利用和競爭可能性,而物種都是生存棲息在水平和垂直的空間中,因此水平-垂直二維生態(tài)位更能準確反映物種對水平空間和垂直空間資源的利用程度(王剛等,1984)。本研究中,巴松錯夏季浮游植物群落優(yōu)勢種在垂直維度上生態(tài)位寬度值差別不大(極差:0.255),故水平-垂直二維生態(tài)位寬度的變化幅度(極差:0.241)受水平維度生態(tài)位寬度(極差:0.409)影響更明顯,即水平維度對巴松錯夏季浮游植物優(yōu)勢種的分布范圍和均勻程度影響更大。不同維度生態(tài)位寬度聚類分析結果表明,小環(huán)藻屬一種為廣生態(tài)位種,主要是因為小環(huán)藻屬一種數(shù)量較多、出現(xiàn)頻率較高,具有較強的資源利用和生態(tài)適應能力,對維持群落的結構和群落內部環(huán)境的形成以及生態(tài)系統(tǒng)的穩(wěn)定起著重要作用。

    生態(tài)位重疊值反映了物種在資源利用或環(huán)境適應能力上的相似性特征和競爭關系(安瑞志等,2020)。在水平維度上,優(yōu)勢種的生態(tài)位重疊值達到有意義的種對比例較高,具有明顯的同質性,表明浮游植物優(yōu)勢種在水平空間分布上相似性極高,種內和種間都存在強烈的競爭關系,這主要是因為巴松錯夏季浮游植物群落大部分優(yōu)勢種對水平維度資源的利用接近程度較高,因此呈現(xiàn)出水平生態(tài)位重疊較高的現(xiàn)象。在垂直維度上,優(yōu)勢種的生態(tài)位重疊度達到有意義的種對比例較低,具有明顯的異質性,這可能是因為溫躍層的存在致使上下水層的物理化學性質差異明顯,形成了不同的生境(安瑞志等,2021a),種間競爭不激烈,物種組成和優(yōu)勢種等群落結構在短期內可能不會有明顯變動,處于相對穩(wěn)定狀態(tài),這與王祥福等(2008)的研究和姚和金等(2008)的研究結果較一致。在水平-垂直維度上,優(yōu)勢種的生態(tài)位重疊達到有意義的種對比例較低,具有明顯的異質性,表明種群間的生態(tài)相似性較大,利用資源的相似性程度較高,這主要是因為水平-垂直二維生態(tài)位重疊度主要受垂直維度分布影響,因此呈現(xiàn)出浮游植物優(yōu)勢種的水平-垂直二維生態(tài)位重疊較低的現(xiàn)象,這與本研究中水平-垂直生態(tài)位寬度值的變化主要受水平生態(tài)位寬度值影響的結果相反。

    3.2? ?優(yōu)勢種的生態(tài)響應速率和群落多穩(wěn)態(tài)

    種群呈發(fā)展或衰退趨勢是種間競爭的直接結果,是生物與非生物因子之間的復雜和多維的互作關系的生態(tài)響應(龐春花等,2015)。在水平維度上,尖針桿藻有較大的發(fā)展空間,發(fā)展趨勢較弱;而在垂直和水平-垂直維度上,尖針桿藻的發(fā)展空間最大,發(fā)展趨勢最強;這表明水平維度上,尖針桿藻與其他優(yōu)勢種在資源競爭上占據(jù)相對優(yōu)勢地位,但發(fā)展趨勢變緩,為了獲得更多的生態(tài)資源可能會通過在垂直維度上擴大生態(tài)位寬度來增強種間競爭力,從而使其在水平-垂直維度獲取最大的發(fā)展空間。同樣,衰退型種群,如針桿藻屬一種在水平、垂直以及水平-垂直維度上均為衰退優(yōu)勢種,但是其生態(tài)響應速率較小,表明其衰退趨勢也會變緩。在水平維度上,脆桿藻屬一種有最小的發(fā)展空間,但發(fā)展趨勢最大;小環(huán)藻屬一種有最小的衰退空間,但衰退趨勢最強;垂直維度上,雙月短縫藻有最小的衰退空間,但衰退趨勢最強;在水平-垂直維度上,肘狀針桿藻有最大的衰退空間,衰退趨勢最強;在垂直和水平-垂直維度上,尖針桿藻發(fā)展空間最大,發(fā)展趨勢最大。以上結果表明,物種的生態(tài)位寬度和生態(tài)響應速率并無直接聯(lián)系,生態(tài)位較寬的物種呈現(xiàn)明顯的衰退趨勢,如小環(huán)藻屬一種在水平、垂直和水平維度生態(tài)位寬度均較大,物種豐度接近飽和,由于自疏或他疏作用使其種群數(shù)量增長緩慢甚至有所減少,早期發(fā)展起來的優(yōu)勢種將逐漸衰退;而生態(tài)位較窄的物種,如尖針桿藻為了獲得更多的資源,可能會通過擴大生態(tài)位寬度來提高自身的生態(tài)適應,其他衰退型的物種將逐漸被替代,這與楊文煥等(2020)和馬一明等(2021)的研究結果較一致。

    生態(tài)響應主要是通過反饋調節(jié)完成的,分為負反饋和正反饋2種,負反饋作用能夠使群落達到和保持平衡或穩(wěn)態(tài),而正反饋作用往往使群落遠離平衡狀態(tài)或穩(wěn)態(tài)(尚玉昌,2010)。本研究中,通過計算生態(tài)響應速率之和,可對群落的反饋機制進行判斷,以此來監(jiān)測群落的穩(wěn)定狀況。巴松錯夏季浮游植物群落優(yōu)勢種生態(tài)響應速率在水平維度上形成正反饋,表明群落遠離平衡狀態(tài)或穩(wěn)態(tài),可能的原因是巴松錯周圍旅游活動頻繁和人為干擾等,加上雨水淋溶和下游筑壩形成水庫的影響(安瑞志等,2021a),致使優(yōu)勢種獲得更多的優(yōu)勢,其他物種受到排擠與去除,多樣性減少,導致群落穩(wěn)定性降低。優(yōu)勢種生態(tài)響應速率在垂直和水平-垂直維度上形成負反饋,表明群落達到和保持平衡或穩(wěn)態(tài),優(yōu)勢種通過自身的響應減緩群落內的壓力,以維持群落的穩(wěn)態(tài)。因此,巴松錯夏季浮游植物群落存在群落多穩(wěn)態(tài),群落在垂直和水平-垂直維度上較水平維度更穩(wěn)定。

    3.3? ?優(yōu)勢種的種間聯(lián)結性

    種間聯(lián)結性大小反映物種間的關系和群落演替的相互作用,總體關聯(lián)性反映了群落演替階段的穩(wěn)定性(劉潤紅等,2020)。本研究采用方差比率法(VR)對巴松錯夏季浮游植物群落優(yōu)勢種的總體關聯(lián)性進行測定,結合[χ]2檢驗定性判定種間聯(lián)結性,并以聯(lián)結系數(shù)(AC)和共同出現(xiàn)百分率(PC)為輔助參數(shù),并經(jīng)過點相關系數(shù)(φ)的檢驗測定,獲得更加精確和可信的結果(程瑞梅等,2013;董靜瑞等,2017;戴冬旭等,2019)。巴松錯夏季浮游植物群落優(yōu)勢種總體關聯(lián)性呈不顯著正關聯(lián),表明群落總體處于演替的初期階段且存在正向演替的趨勢。[χ]2檢驗判斷了巴松錯夏季浮游植物群落8個優(yōu)勢種的種間聯(lián)結顯著程度。從整體上看,巴松錯夏季浮游植物群落的正負聯(lián)結比為0.93,負聯(lián)結占優(yōu)勢,說明巴松錯夏季浮游植物群落比較脆弱,結構簡單,仍處于演替的初期階段。但大部分種間聯(lián)結性均為未達到顯著水平,表明巴松錯夏季浮游植物優(yōu)勢種種間聯(lián)結性不強,種間關系較為松散,可能是因為選取的優(yōu)勢種物種數(shù)相對較少,且各自占據(jù)有利的資源空間,同域率低,相互依賴程度較低引起的;另一方面可能與物種的生物學特性、空間分布的異質性及多樣性有關,浮游植物個體較小,生命周期短,營浮游生活,繁殖快,對環(huán)境變化敏感(胡鴻鈞和魏印心, 2016),加之由于人為活動、極端氣候等因素的干擾,導致浮游植物群落呈現(xiàn)不穩(wěn)定的波動狀態(tài)并伴有演替現(xiàn)象發(fā)生,從而使優(yōu)勢種的種間聯(lián)結性也并不顯著,這與楊文煥等(2020)和陶敏等(2021)的研究結果較一致。

    在[χ]2檢驗的基礎上,結合聯(lián)結系數(shù)(AC)和共同出現(xiàn)百分率(PC)分析,不僅能反映種間聯(lián)結的真實情況,而且是對[χ]2檢驗中不顯著的種間聯(lián)結性的進一步驗證,從而得到更精確的結果(安瑞志等,2021a)。聯(lián)結系數(shù)(AC)結果顯示,2對正聯(lián)結系數(shù)程度較高(AC≥0.6),5對負聯(lián)結系數(shù)程度較高(AC<-0.6)。這個結果進一步驗證了[χ]2檢驗中不顯著的種間聯(lián)結性,表明大部分優(yōu)勢種具有相似的生境需求或生物學特性,并且生態(tài)位重疊值也較高,這與申涵(2020)和安瑞志等(2021a)的研究結果基本一致。根據(jù)共同出現(xiàn)百分率(PC)可知,聯(lián)結程度較強的種對僅有1對,為小環(huán)藻屬一種和肘狀針桿藻形成的種對,這可能是因為盡管它們的空間重合度高,對生境有共同資源需求,但資源利用方式存在不同,利用方式表現(xiàn)為互補而非種間競爭,趨向于物種間的穩(wěn)定共存(李丘霖等,2017)。在種間聯(lián)結性中,聯(lián)結系數(shù)(AC)和共同出現(xiàn)百分率(PC)能體現(xiàn)出那些由[χ]2檢驗判斷為不顯著種對間的聯(lián)結性及大小,但共同出現(xiàn)百分率(PC)夸大了2個物種均出現(xiàn)時的聯(lián)結顯著程度,聯(lián)結系數(shù)(AC)夸大了2個物種均不出現(xiàn)時的聯(lián)結顯著程度, 而點相關系數(shù)(φ)可以降低這2種情況的影響(潘琪等,2015)。經(jīng)過點相關系數(shù)(φ)的測定,小環(huán)藻屬一種與尖針桿藻形成的種對聯(lián)結系數(shù)(AC)與共同出現(xiàn)百分率(PC)均有明顯降低。因此,點相關系數(shù)(φ)可以對[χ]2檢驗、聯(lián)結系數(shù)(AC)和共同出現(xiàn)百分率(PC)測定結果對種間聯(lián)結性進行校正,使研究結果更精確。

    綜合比較x2檢驗、聯(lián)結系數(shù)(AC)、共同出現(xiàn)百分率(PC)及點相關系數(shù)(φ)的研究結果,由于不同的測定方法,所得結果之間存在一定差異,但整體上是趨于一致的(申涵,2020)。比如小環(huán)藻屬一種和肘狀針桿藻形成的種對,[χ]2=0.931 (ad>bc),表明種對間呈不顯著正聯(lián)結性;聯(lián)結系數(shù)(AC)值為0.219,表明種對間負聯(lián)結系數(shù)程度一般;共同出現(xiàn)百分率(PC)值呈聯(lián)結程度較強,表明二者之間同時出現(xiàn)的幾率較大;點相關系數(shù)(φ)值為0.3,則表明種對間正聯(lián)結程度一般。因此,結合不同的指標綜合分析,可以更科學地判定種間聯(lián)結性。

    參考文獻

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

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

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

    陳龍濤,石曉東,高潤梅,2016.山西陵川南方紅豆杉群落種間聯(lián)結與生態(tài)位特征研究[J].植物科學學報,34(4):521-529.

    程瑞梅,王瑞麗,劉澤彬,等,2013.三峽庫區(qū)栲屬群落主要喬木種群的種間聯(lián)結性[J]. 林業(yè)科學,49(5):36-42.

    戴冬旭,俞存根,劉惠,等,2019.南鹿列島東側海域主要游泳動物生態(tài)位及種間聯(lián)結性[J]. 應用生態(tài)學報,30(12):4249-4258.

    董靜瑞,水柏年,胡成業(yè),等,2017.溫州南部沿岸海域主要魚類的生態(tài)位及種間聯(lián)結性[J].應用生態(tài)學報,28(5):1699-1706.

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

    李港,林妙麗,陳誠,等,2021.京杭運河蘇北段秋冬季浮游植物群落結構特征及其與環(huán)境因子的關系[J].水生態(tài)學雜志,42(5):119-126.

    江煥,張輝,龍文興,等,2019.金鐘藤入侵群落的種間聯(lián)結及生態(tài)位特征[J].生物多樣性,27(4):388-399.

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

    李丘霖,宗秀虹,鄧洪平,等,2017.赤水桫欏群落喬木層優(yōu)勢物種生態(tài)位與種間聯(lián)結性研究[J].西北植物學報,37(7):1422-1428.

    劉潤紅,陳樂,涂洪潤,等,2020.桂林巖溶石山青岡群落灌木層主要物種生態(tài)位與種間聯(lián)結[J].生態(tài)學報,40(6) : 2057-2071.

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

    劉洋,安瑞志,王陳,等,2023.西藏巴松錯夏季浮游植物群落的分布[J].水生態(tài)學雜志,44(6):104-112.

    馬一明,李秋華,潘少樸,等,2021.貴州高原花溪水庫浮游植物優(yōu)勢種生態(tài)位及種間聯(lián)結性動態(tài)分析[J].湖泊科學,33(3):785-796.

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

    潘琪,梁軍,倪楊,等,2015.昆崳山腮扁葉蜂與喬木層樹種的聯(lián)結性分析[J]. 林業(yè)科學,51(10) : 85-92.

    龐春花,范曉,張峰,等,2015.不同資源維度上汾河流域下游優(yōu)勢種的生態(tài)位[J].生態(tài)學雜志,34(2):380-386.

    尚玉昌,2010.普通生態(tài)學[M].北京:北京大學出版社.

    申涵,2020. 寒旱區(qū)湖泊浮游植物生態(tài)位及種間聯(lián)結特征分析[D].包頭:內蒙古科技大學.

    陶敏,岳興建,岳珊,等,2021.四川丘陵區(qū)水庫浮游植物群落結構與藍藻水華風險——基于優(yōu)勢種生態(tài)位與種間聯(lián)結研究[J].生態(tài)學報,(23):1-13.

    王剛,趙松嶺,張鵬云,等,1984.關于生態(tài)位定義的探討及生態(tài)位重疊計測公式改進的研究[J].生態(tài)學報,(2):119-127.

    王捷,李博,馮佳,等,2015.西藏西南部湖泊浮游藻類區(qū)系及群落結構特征[J].水生生物學報,39(4):837-844.

    王祥福,郭泉水,巴哈爾古麗,等,2008.崖柏群落優(yōu)勢喬木種群生態(tài)位[J].林業(yè)科學(4):6-13.

    汪志聰,吳衛(wèi)菊,左明,等, 2010.巢湖浮游植物群落生態(tài)位的研究[J]. 長江流域資源與環(huán)境,19(6):685-691.

    吳佳夢,徐娜娜,張文珺,等,2019.浙江舟山定海護城河浮游植物優(yōu)勢種生態(tài)位與種間聯(lián)結性季節(jié)性分析[J].湖泊科學,31(2):429-439.

    肖紅,周啟星,曹瑩,等,2003.沈陽地區(qū)水田主要雜草種群的消長動態(tài)及生態(tài)位分析[J].農(nóng)村生態(tài)環(huán)境, (3):9-13.

    楊文煥,申涵,周明利,等,2020.包頭南海湖浮游植物優(yōu)勢種生態(tài)位及種間聯(lián)結性季節(jié)分析[J].中國環(huán)境科學,40(1):383-391.

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

    姚和金,金宗來,楊偉斌, 等, 2008.浙西南紅黃壤果園雜草種群消長動態(tài)及生態(tài)位的研究[J].江蘇農(nóng)業(yè)學報,(5):649-655.

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

    張苗苗,王詠雪,田闊,等,2018.臺州玉環(huán)北部沿岸海域主要游泳動物生態(tài)位和種間聯(lián)結性[J].應用生態(tài)學報,29(11):3867-3875.

    趙萍萍,衛(wèi)志宏,吳秋婷,等,2012.洱海浮游橈足類群落結構對季節(jié)性休漁的響應[J].應用與環(huán)境生物學報,18(3) : 421-425.

    Avolio M L, Forrestel E J, Chang C C, et al, 2019.Demystifying dominant species[J].New Phytologist, 223(3):1106-1126.

    Cody M L, 1974. Competition and Structure of Bird Communities[M]. Princeton: Princeton University Press.

    Colwell R K,F(xiàn)utuyma D J, 1971. On the measurement of Niche Breadth and Overlap[J]. Ecology,52(4):567-576.

    Habib O A , Tippett R , Murphy K J, 1997. Seasonal changes in phytoplankton community structure in relation to physicochemical factors in Loch Lomond, Scotland[J]. Hydrobiologia,350(1/2/3):63-79.

    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.

    May R M, 1975. Some notes on estimating the competition matrix[J]. Ecology, 56: 737-74.

    Mcnaughton S J. 1967. Relationships among functional properties of Californian grassland[J]. Nature, 216(5111) : 168-169.

    Pianka E R, 1973. The structure of lizard communities[J]. Annual Review of Ecology and Systematics,4:53-74.

    Schluter D,1984. A variance test for detecting species associations, with some example applications[J]. Ecology,65(3):998-1005.

    (責任編輯? ?鄭金秀)

    收稿日期:2022-03-14? ? ? 修回日期:2023-12-19

    基金項目:國家自然科學基金(31660620,32070418);2022 年中央財政支持地方高校改革發(fā)展專項資金 (藏財預指[2022]1號) 。

    作者簡介:劉洋,女,1988年生,實驗師,主要從事高寒水體浮游植物生態(tài)學研究。E-mail:ly@utibet.edu.cn

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

    猜你喜歡
    浮游植物生態(tài)位
    生態(tài)位視域下大學生就業(yè)支持研究
    生態(tài)位視角下物流產(chǎn)業(yè)競爭力評價
    塔里木河阿拉爾段及多浪水庫浮游生物的調查研究
    基于生態(tài)位的企業(yè)成長理論研究
    時代金融(2016年27期)2016-11-25 17:07:08
    株洲市房地產(chǎn)業(yè)住宅產(chǎn)品結構的生態(tài)位診斷及應對
    商(2016年31期)2016-11-22 22:17:08
    遼東灣網(wǎng)采浮游植物名錄
    生態(tài)位視野下純網(wǎng)綜藝優(yōu)勢初探
    今傳媒(2016年5期)2016-06-01 00:10:32
    特征色素在淡水魚養(yǎng)殖模式中對浮游植物群落結構的影響
    渭河陜西段浮游植物群落結構及水質評價
    濟南流域春季浮游植物多樣性的研究
    av电影中文网址| 亚洲中文字幕日韩| 欧美在线黄色| 精品久久久久久久毛片微露脸| 一区二区三区精品91| av网站在线播放免费| 黄色毛片三级朝国网站| 黄色成人免费大全| 国产欧美日韩一区二区三| 精品免费久久久久久久清纯 | 久久久久国产一级毛片高清牌| 久久国产亚洲av麻豆专区| 麻豆乱淫一区二区| 18禁国产床啪视频网站| 嫁个100分男人电影在线观看| 欧美另类亚洲清纯唯美| 美国免费a级毛片| 久久久久久久国产电影| 19禁男女啪啪无遮挡网站| 免费人妻精品一区二区三区视频| 精品第一国产精品| 一级a爱视频在线免费观看| 搡老岳熟女国产| 国产一区二区三区综合在线观看| 精品熟女少妇八av免费久了| 淫妇啪啪啪对白视频| 亚洲五月色婷婷综合| 久久精品熟女亚洲av麻豆精品| 日韩欧美免费精品| 精品久久蜜臀av无| a级毛片黄视频| 午夜福利视频在线观看免费| 精品免费久久久久久久清纯 | 国产成+人综合+亚洲专区| 狠狠精品人妻久久久久久综合| 18禁国产床啪视频网站| 亚洲人成电影观看| 国产成人av激情在线播放| 久热爱精品视频在线9| 男女无遮挡免费网站观看| 午夜福利免费观看在线| 18禁裸乳无遮挡动漫免费视频| 少妇猛男粗大的猛烈进出视频| 久久婷婷成人综合色麻豆| 国产伦理片在线播放av一区| 后天国语完整版免费观看| 最新美女视频免费是黄的| 在线观看66精品国产| 一区福利在线观看| 亚洲一码二码三码区别大吗| 夫妻午夜视频| 黄网站色视频无遮挡免费观看| 国产av一区二区精品久久| 欧美精品一区二区大全| 中文字幕最新亚洲高清| 日本黄色视频三级网站网址 | 18在线观看网站| 免费av中文字幕在线| 亚洲午夜精品一区,二区,三区| 在线观看一区二区三区激情| 母亲3免费完整高清在线观看| 高清av免费在线| 岛国在线观看网站| 男女无遮挡免费网站观看| 国产精品久久久久久人妻精品电影 | www.熟女人妻精品国产| 精品国产超薄肉色丝袜足j| 国产一卡二卡三卡精品| 波多野结衣一区麻豆| 亚洲欧美日韩高清在线视频 | 免费不卡黄色视频| 久久精品亚洲熟妇少妇任你| 午夜福利一区二区在线看| av超薄肉色丝袜交足视频| 午夜成年电影在线免费观看| 激情视频va一区二区三区| 国产不卡av网站在线观看| 1024视频免费在线观看| 成人18禁在线播放| 变态另类成人亚洲欧美熟女 | 一级毛片精品| 丝袜喷水一区| 久久人妻福利社区极品人妻图片| 18禁美女被吸乳视频| 久久精品国产99精品国产亚洲性色 | 成人国产一区最新在线观看| 国产亚洲精品久久久久5区| 一本色道久久久久久精品综合| 一区二区三区乱码不卡18| 精品人妻1区二区| 免费观看人在逋| 一本色道久久久久久精品综合| 精品少妇内射三级| 99re在线观看精品视频| 一区二区三区精品91| 高潮久久久久久久久久久不卡| 妹子高潮喷水视频| 亚洲av电影在线进入| 性少妇av在线| 九色亚洲精品在线播放| 91精品三级在线观看| 久9热在线精品视频| 日日夜夜操网爽| 在线亚洲精品国产二区图片欧美| 精品久久久久久久毛片微露脸| 日韩有码中文字幕| 欧美人与性动交α欧美软件| 在线观看一区二区三区激情| 成人手机av| 亚洲色图综合在线观看| 伊人久久大香线蕉亚洲五| 亚洲av电影在线进入| 18禁裸乳无遮挡动漫免费视频| 在线十欧美十亚洲十日本专区| 国产主播在线观看一区二区| 波多野结衣一区麻豆| 亚洲中文日韩欧美视频| 国产淫语在线视频| 午夜精品久久久久久毛片777| 91精品国产国语对白视频| 精品人妻1区二区| 日本a在线网址| 亚洲avbb在线观看| 免费在线观看影片大全网站| 男女之事视频高清在线观看| 最新美女视频免费是黄的| 999精品在线视频| 岛国毛片在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 水蜜桃什么品种好| 国产精品.久久久| 日本vs欧美在线观看视频| 免费在线观看视频国产中文字幕亚洲| 精品乱码久久久久久99久播| 黄色片一级片一级黄色片| 欧美乱码精品一区二区三区| 免费少妇av软件| 久久亚洲真实| 黄色毛片三级朝国网站| 精品一品国产午夜福利视频| 男人操女人黄网站| 亚洲中文日韩欧美视频| 日韩欧美一区二区三区在线观看 | 日本撒尿小便嘘嘘汇集6| 一本大道久久a久久精品| 男女无遮挡免费网站观看| 一区二区三区精品91| 丁香六月欧美| 久久久国产欧美日韩av| 国产色视频综合| 国产免费视频播放在线视频| 欧美变态另类bdsm刘玥| 国产一区二区在线观看av| 亚洲自偷自拍图片 自拍| 午夜福利乱码中文字幕| 国产欧美日韩综合在线一区二区| av电影中文网址| 国产av一区二区精品久久| 一区在线观看完整版| 丝袜美足系列| 一进一出好大好爽视频| 日韩一卡2卡3卡4卡2021年| 国产老妇伦熟女老妇高清| 久久毛片免费看一区二区三区| 国产熟女午夜一区二区三区| 欧美乱码精品一区二区三区| 日本精品一区二区三区蜜桃| 日韩大片免费观看网站| 天堂俺去俺来也www色官网| 黄色视频不卡| 国产成人啪精品午夜网站| 另类精品久久| 日韩欧美免费精品| 国产精品电影一区二区三区 | 叶爱在线成人免费视频播放| 精品一区二区三区av网在线观看 | 国产熟女午夜一区二区三区| 又黄又粗又硬又大视频| 99国产精品免费福利视频| 久久久精品国产亚洲av高清涩受| 国产日韩欧美视频二区| 免费观看av网站的网址| 两个人看的免费小视频| 啦啦啦免费观看视频1| 亚洲七黄色美女视频| 国产1区2区3区精品| 国产野战对白在线观看| 制服诱惑二区| 欧美日韩av久久| 十八禁网站免费在线| 一本—道久久a久久精品蜜桃钙片| 国产亚洲av高清不卡| 亚洲人成77777在线视频| 丰满人妻熟妇乱又伦精品不卡| 每晚都被弄得嗷嗷叫到高潮| a级毛片在线看网站| 亚洲色图综合在线观看| 久久天堂一区二区三区四区| 国产麻豆69| 九色亚洲精品在线播放| 久久中文字幕一级| 免费在线观看日本一区| 两个人免费观看高清视频| 99久久国产精品久久久| 久久久久久免费高清国产稀缺| 久久人妻福利社区极品人妻图片| 一级毛片精品| 国产高清激情床上av| 99精品久久久久人妻精品| 超碰成人久久| 午夜福利免费观看在线| 日本av手机在线免费观看| 男女午夜视频在线观看| 亚洲人成电影免费在线| 亚洲五月色婷婷综合| 国产精品自产拍在线观看55亚洲 | 亚洲精品成人av观看孕妇| 女人久久www免费人成看片| 一级片'在线观看视频| 免费在线观看影片大全网站| 99riav亚洲国产免费| 天天躁日日躁夜夜躁夜夜| 欧美老熟妇乱子伦牲交| 人妻久久中文字幕网| 国产成人欧美| 超色免费av| 亚洲午夜理论影院| 天堂动漫精品| 国产高清国产精品国产三级| 天天影视国产精品| 日本撒尿小便嘘嘘汇集6| 一区二区三区激情视频| 交换朋友夫妻互换小说| 一本久久精品| 十八禁网站网址无遮挡| 美女扒开内裤让男人捅视频| 新久久久久国产一级毛片| 国产黄频视频在线观看| 高清欧美精品videossex| 国产欧美日韩精品亚洲av| 精品午夜福利视频在线观看一区 | 汤姆久久久久久久影院中文字幕| 欧美精品啪啪一区二区三区| 国产欧美亚洲国产| 国产高清videossex| 亚洲伊人色综图| 国产成人影院久久av| 精品久久久精品久久久| 欧美黄色片欧美黄色片| 精品久久蜜臀av无| 18禁观看日本| 制服诱惑二区| 黄网站色视频无遮挡免费观看| 精品免费久久久久久久清纯 | 妹子高潮喷水视频| 国产成人精品无人区| 国产99久久九九免费精品| 国产高清激情床上av| 亚洲成a人片在线一区二区| 一区二区三区激情视频| 1024视频免费在线观看| 老熟妇乱子伦视频在线观看| 国产免费福利视频在线观看| 9191精品国产免费久久| 视频区欧美日本亚洲| 亚洲精品国产色婷婷电影| 午夜免费成人在线视频| 在线观看免费午夜福利视频| 动漫黄色视频在线观看| 麻豆av在线久日| 中文字幕精品免费在线观看视频| 91成年电影在线观看| 成人18禁高潮啪啪吃奶动态图| 免费少妇av软件| 国产深夜福利视频在线观看| 午夜免费成人在线视频| 99久久国产精品久久久| 狠狠精品人妻久久久久久综合| 欧美亚洲 丝袜 人妻 在线| 欧美 亚洲 国产 日韩一| 女人久久www免费人成看片| 久久影院123| av在线播放免费不卡| 色婷婷av一区二区三区视频| 日韩制服丝袜自拍偷拍| 别揉我奶头~嗯~啊~动态视频| 亚洲精品一卡2卡三卡4卡5卡| 国产av精品麻豆| 免费在线观看日本一区| 精品少妇黑人巨大在线播放| 国产在线精品亚洲第一网站| 亚洲成人国产一区在线观看| 成人18禁在线播放| 久热这里只有精品99| 一本色道久久久久久精品综合| 久久久久久久国产电影| 欧美日韩福利视频一区二区| 婷婷成人精品国产| 青草久久国产| 一区二区三区激情视频| 亚洲人成电影观看| 国产精品一区二区精品视频观看| 99国产精品一区二区蜜桃av | 亚洲欧美精品综合一区二区三区| 亚洲av第一区精品v没综合| 午夜福利视频精品| 19禁男女啪啪无遮挡网站| 国产精品一区二区免费欧美| 日韩欧美三级三区| 国产精品一区二区免费欧美| 国产精品二区激情视频| 人人妻人人添人人爽欧美一区卜| 日韩欧美三级三区| 国产精品一区二区免费欧美| 男女床上黄色一级片免费看| 国产一区二区 视频在线| 两个人免费观看高清视频| 777米奇影视久久| 精品亚洲乱码少妇综合久久| 免费不卡黄色视频| 精品亚洲成a人片在线观看| 欧美日韩一级在线毛片| 丰满饥渴人妻一区二区三| 老鸭窝网址在线观看| 国产精品熟女久久久久浪| 亚洲九九香蕉| 久久精品亚洲精品国产色婷小说| 成在线人永久免费视频| 一本综合久久免费| 久久人妻福利社区极品人妻图片| 一本大道久久a久久精品| 国产精品自产拍在线观看55亚洲 | 国产在线一区二区三区精| 免费观看a级毛片全部| 亚洲精品在线美女| 亚洲天堂av无毛| 亚洲熟女毛片儿| 中国美女看黄片| 国产精品98久久久久久宅男小说| 久久久久久久久久久久大奶| 国产精品98久久久久久宅男小说| 久久精品成人免费网站| 精品少妇一区二区三区视频日本电影| 超色免费av| 日韩一区二区三区影片| 久久精品亚洲av国产电影网| 亚洲精品在线美女| 69精品国产乱码久久久| 老熟妇仑乱视频hdxx| 欧美乱妇无乱码| 日韩一区二区三区影片| 日本av免费视频播放| 亚洲中文字幕日韩| 后天国语完整版免费观看| 久久久久久免费高清国产稀缺| 大香蕉久久成人网| 亚洲国产欧美日韩在线播放| av福利片在线| 成人国产一区最新在线观看| 女人精品久久久久毛片| 国产日韩欧美在线精品| 午夜久久久在线观看| 国产真人三级小视频在线观看| 国产精品亚洲av一区麻豆| 涩涩av久久男人的天堂| 一本色道久久久久久精品综合| 色在线成人网| 高清av免费在线| 一区福利在线观看| 国产高清国产精品国产三级| 精品少妇黑人巨大在线播放| 亚洲精品久久午夜乱码| 亚洲精品国产精品久久久不卡| 国精品久久久久久国模美| 国产一区二区激情短视频| 日日夜夜操网爽| 亚洲中文av在线| 啦啦啦在线免费观看视频4| 十八禁网站网址无遮挡| 在线观看免费高清a一片| 午夜福利免费观看在线| 无限看片的www在线观看| 18禁美女被吸乳视频| 色婷婷久久久亚洲欧美| 国产精品 国内视频| 变态另类成人亚洲欧美熟女 | 亚洲少妇的诱惑av| 在线av久久热| 精品一品国产午夜福利视频| 国产aⅴ精品一区二区三区波| 久久av网站| 人人澡人人妻人| 91精品三级在线观看| 日日夜夜操网爽| 日韩熟女老妇一区二区性免费视频| 中文字幕高清在线视频| 久久人妻av系列| 人人妻人人爽人人添夜夜欢视频| 亚洲精品av麻豆狂野| 亚洲精品美女久久久久99蜜臀| 国精品久久久久久国模美| 久久精品人人爽人人爽视色| 亚洲自偷自拍图片 自拍| 高清黄色对白视频在线免费看| 色老头精品视频在线观看| 老鸭窝网址在线观看| 99国产精品一区二区三区| 亚洲精品美女久久av网站| 中文字幕av电影在线播放| 国产一区二区激情短视频| 91精品国产国语对白视频| 国产一区二区三区在线臀色熟女 | 一区二区三区国产精品乱码| 久久中文字幕人妻熟女| 色婷婷久久久亚洲欧美| 最新的欧美精品一区二区| 超碰97精品在线观看| av福利片在线| 一个人免费看片子| svipshipincom国产片| 久久毛片免费看一区二区三区| 欧美久久黑人一区二区| 亚洲中文av在线| 女人爽到高潮嗷嗷叫在线视频| 免费观看人在逋| 下体分泌物呈黄色| 高清毛片免费观看视频网站 | 精品少妇内射三级| 精品免费久久久久久久清纯 | 国产欧美日韩一区二区三| 91老司机精品| 亚洲av日韩精品久久久久久密| 99久久99久久久精品蜜桃| 精品国产国语对白av| 午夜福利在线免费观看网站| 老司机深夜福利视频在线观看| 天天操日日干夜夜撸| 麻豆成人av在线观看| 成人黄色视频免费在线看| 午夜精品国产一区二区电影| 少妇的丰满在线观看| 国产精品香港三级国产av潘金莲| 亚洲情色 制服丝袜| 亚洲精品粉嫩美女一区| 欧美乱码精品一区二区三区| 久久久久久久久久久久大奶| 婷婷成人精品国产| 亚洲成a人片在线一区二区| 日韩 欧美 亚洲 中文字幕| 中文字幕最新亚洲高清| 性高湖久久久久久久久免费观看| 王馨瑶露胸无遮挡在线观看| 天天操日日干夜夜撸| 天堂俺去俺来也www色官网| 欧美日韩福利视频一区二区| 在线观看一区二区三区激情| 久久影院123| 美女高潮到喷水免费观看| 99九九在线精品视频| 18禁美女被吸乳视频| 变态另类成人亚洲欧美熟女 | 亚洲美女黄片视频| 久久久精品94久久精品| 国产一区二区激情短视频| 一本大道久久a久久精品| 亚洲精华国产精华精| 精品一区二区三区视频在线观看免费 | 国产精品影院久久| 国产精品免费大片| 精品人妻熟女毛片av久久网站| 精品卡一卡二卡四卡免费| 日本a在线网址| 亚洲成人免费电影在线观看| 女警被强在线播放| 亚洲第一av免费看| 最新美女视频免费是黄的| 成人亚洲精品一区在线观看| 成年版毛片免费区| 欧美人与性动交α欧美精品济南到| 人人妻,人人澡人人爽秒播| 2018国产大陆天天弄谢| 精品少妇黑人巨大在线播放| 久久精品亚洲精品国产色婷小说| a在线观看视频网站| 久久人人97超碰香蕉20202| 精品午夜福利视频在线观看一区 | 午夜成年电影在线免费观看| 亚洲精品粉嫩美女一区| 成年人午夜在线观看视频| 香蕉国产在线看| videos熟女内射| 国产成人啪精品午夜网站| 国产一区二区三区综合在线观看| 超碰成人久久| 一区二区三区精品91| 久久精品国产a三级三级三级| 亚洲综合色网址| 男女床上黄色一级片免费看| 久久精品亚洲精品国产色婷小说| a在线观看视频网站| 亚洲av成人不卡在线观看播放网| 2018国产大陆天天弄谢| 成年人免费黄色播放视频| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩欧美三级三区| 精品久久久精品久久久| 18在线观看网站| av线在线观看网站| 欧美中文综合在线视频| av有码第一页| 老司机在亚洲福利影院| 国产极品粉嫩免费观看在线| 在线观看免费视频网站a站| 亚洲色图综合在线观看| 大型av网站在线播放| 老司机靠b影院| 亚洲国产欧美在线一区| 老汉色av国产亚洲站长工具| 免费在线观看视频国产中文字幕亚洲| 成人手机av| 丝袜喷水一区| 淫妇啪啪啪对白视频| 91av网站免费观看| 成人特级黄色片久久久久久久 | 欧美日韩亚洲高清精品| 欧美成人免费av一区二区三区 | 国产免费福利视频在线观看| 久久婷婷成人综合色麻豆| 老司机午夜福利在线观看视频 | 伦理电影免费视频| 天天添夜夜摸| 亚洲熟女毛片儿| 久久久久久亚洲精品国产蜜桃av| 一进一出好大好爽视频| 丰满少妇做爰视频| 中国美女看黄片| 精品亚洲乱码少妇综合久久| 91大片在线观看| 叶爱在线成人免费视频播放| 在线观看www视频免费| e午夜精品久久久久久久| 久久天堂一区二区三区四区| 国产亚洲精品一区二区www | 一级毛片女人18水好多| 男女之事视频高清在线观看| 悠悠久久av| 欧美成人免费av一区二区三区 | 在线看a的网站| 女人被躁到高潮嗷嗷叫费观| 国产成人啪精品午夜网站| videosex国产| 久热爱精品视频在线9| 老司机午夜十八禁免费视频| 1024香蕉在线观看| 无限看片的www在线观看| 亚洲色图 男人天堂 中文字幕| 国产xxxxx性猛交| 免费在线观看日本一区| 欧美精品人与动牲交sv欧美| 99久久人妻综合| 正在播放国产对白刺激| 99精国产麻豆久久婷婷| 超碰97精品在线观看| 成人三级做爰电影| 五月开心婷婷网| 欧美精品亚洲一区二区| 久久久精品国产亚洲av高清涩受| 91麻豆精品激情在线观看国产 | 真人做人爱边吃奶动态| 1024视频免费在线观看| 亚洲美女黄片视频| 看免费av毛片| 亚洲精华国产精华精| 涩涩av久久男人的天堂| 精品国产乱子伦一区二区三区| 亚洲欧美一区二区三区黑人| 亚洲专区中文字幕在线| 在线播放国产精品三级| 免费观看人在逋| 女人被躁到高潮嗷嗷叫费观| 欧美国产精品一级二级三级| 精品久久久久久久毛片微露脸| 建设人人有责人人尽责人人享有的| 成人国产av品久久久| 精品第一国产精品| 狠狠狠狠99中文字幕| 亚洲 国产 在线| 国产片内射在线| 久久人妻熟女aⅴ| 又大又爽又粗| 欧美av亚洲av综合av国产av| 啦啦啦中文免费视频观看日本| 亚洲精华国产精华精| 亚洲美女黄片视频| 窝窝影院91人妻| 又大又爽又粗| 国产日韩欧美在线精品| 午夜福利视频精品| 91国产中文字幕| 国产成人系列免费观看| 午夜福利视频精品| 91国产中文字幕| 在线观看免费视频网站a站| 亚洲 欧美一区二区三区| 少妇裸体淫交视频免费看高清 | 两性夫妻黄色片| 国产主播在线观看一区二区| av电影中文网址| 日韩有码中文字幕| 久久中文看片网| 国产精品熟女久久久久浪| 精品国产一区二区久久| 男女床上黄色一级片免费看| 国产精品98久久久久久宅男小说| 成人国语在线视频|