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

    小秦嶺自然保護(hù)區(qū)苔蘚植物群落數(shù)量分類、排序及多樣性垂直格局

    2017-06-26 11:28:26馮佳偉韓軍旺葉永忠袁志良
    生態(tài)學(xué)報(bào) 2017年8期
    關(guān)鍵詞:苔蘚秦嶺樣地

    陳 云,馮佳偉,牛 帥,許 寧,韓軍旺,葉永忠,袁志良,*

    1 河南農(nóng)業(yè)大學(xué)林學(xué)院,鄭州 450002 2 河南農(nóng)業(yè)大生命科學(xué)學(xué)院,鄭州 450002 3 河南小秦嶺國家級(jí)自然保護(hù)區(qū)管理局,靈寶 472500

    ?

    小秦嶺自然保護(hù)區(qū)苔蘚植物群落數(shù)量分類、排序及多樣性垂直格局

    陳 云1,馮佳偉2,牛 帥2,許 寧2,韓軍旺3,葉永忠2,袁志良2,*

    1 河南農(nóng)業(yè)大學(xué)林學(xué)院,鄭州 450002 2 河南農(nóng)業(yè)大生命科學(xué)學(xué)院,鄭州 450002 3 河南小秦嶺國家級(jí)自然保護(hù)區(qū)管理局,靈寶 472500

    苔蘚物種多樣性沿海拔梯度呈什么樣的分布格局?基于53個(gè)20 m×20 m樣地,采集并鑒定了樣地內(nèi)所有生境中的苔蘚物種;采用多元回歸樹(multivariate regression trees, MRT)對(duì)苔蘚植物進(jìn)行分類;采用典范對(duì)應(yīng)分析(canonical correspondence analysis, CCA)進(jìn)行排序分析;用廣義線性模型(Generalized Linear Model, GLM)研究苔蘚物種多樣性沿海拔的分布格局。結(jié)果如下:(1)采集的1378份苔蘚植物標(biāo)本經(jīng)過鑒定,共有33科89屬240種,其中蘚類226種,苔類14種;(2)經(jīng)交叉驗(yàn)證認(rèn)為本區(qū)苔蘚植物可分為4類,Ⅰ長肋青蘚(Brachytheciumpopuleum)+大葉匐燈蘚(Plagiomniumsucculentum)+圓葉匐燈蘚(Plagiomnium vesicatum) 群落,Ⅱ灰白青蘚(Brachytheciumalbicans)+寬葉青蘚(Brachytheciumcurtum)+短肋羽蘚(Thuidiumkanedae)+平肋提燈蘚(Mniumlaevinerve) 群落,Ⅲ短肋羽蘚(Thuidiumkanedae)+光萼葉苔(Jungermannialeiantha)+薄羅蘚(Leskeapolycarpa)+叉肋蘚(Trachyphylluminflexum) 群落,Ⅳ寬葉青蘚(Brachytheciumcurtum)+垂蒴棉蘚(Plagiotheciumnemorale)+全緣匐燈蘚(Plagiomniumintegrum) 群落;(3)海拔、坡度和坡向?qū)μμ\物種分布的解釋量為7.29%;(4)苔蘚物種多樣性隨海拔的變化成明顯的上升格局。通過研究了解了小秦嶺國家級(jí)自然保護(hù)區(qū)苔蘚物種組成、群叢劃分以及每種苔蘚植物的海拔分布范圍;苔蘚物種多樣性隨海拔的上升呈顯著的上升格局,與喬木、灌木和草本的分布格局不同。希望此研究能為苔蘚植物生態(tài)學(xué)研究和物種多樣性保護(hù)提供參考。

    多元回歸樹;典范對(duì)應(yīng)分析;海拔;物種多樣性;苔蘚

    苔蘚是森林生態(tài)系統(tǒng)的重要組分,在森林動(dòng)態(tài)變化過程中扮演著重要的作用[1]。在林業(yè)上,苔蘚是劃分森林類型的輔助指標(biāo),并且被作為監(jiān)測環(huán)境的指示植物[2]。然而苔蘚植物個(gè)體相對(duì)較小,缺乏廣泛的經(jīng)濟(jì)價(jià)值,因此很少成為人們關(guān)注的焦點(diǎn),其相關(guān)生態(tài)學(xué)研究也相當(dāng)薄弱[3-4]。

    植被群落分類和排序是植被生態(tài)學(xué)研究中最基本的問題[5]。目前,常用的植被群落數(shù)量分類方法有聚類分析、雙向指示種分析和多元回歸樹[6- 7]。其中多元回歸樹是一種較新的群落數(shù)量分類方法,較聚類分析和雙向指示種分析更加客觀[8]。目前圍繞維管植物群落分類進(jìn)行了較多的研究[9-11]。但對(duì)于苔蘚植物分類,通過查閱文獻(xiàn)僅發(fā)現(xiàn),郭水良等首次在國內(nèi)采用雙向指示種分析,對(duì)長白山苔蘚植物進(jìn)行了分類[2]。隨后張?jiān)鞯纫膊捎秒p向指示種分析,對(duì)新疆博格達(dá)山地地面生苔蘚進(jìn)行了分類[12]。苔蘚植物分類研究較少,并且沒有發(fā)現(xiàn)采用多元回歸樹對(duì)苔蘚植物進(jìn)行劃分的相關(guān)研究。

    海拔作為地理梯度,包含了溫度、濕度和光照等各種環(huán)境因子,被認(rèn)為是植物物種、植物區(qū)系、植物群落和生態(tài)系統(tǒng)的決定因素[13]。山地植被物種多樣性隨海拔梯度的變化規(guī)律也一直是生態(tài)學(xué)家感興趣的問題[14-15]。維管植物多樣性沿海拔梯度的分布格局已有了大量的研究,分布格局主要有5種形式,分別是隨海拔升高先降低后升高、先升高后降低、單調(diào)升高、單調(diào)下降和沒有明顯格局[14,16]。然而苔蘚物種多樣性沿海拔梯度呈什么樣的格局?在國內(nèi)沒有發(fā)現(xiàn)相關(guān)報(bào)道。

    本研究基于53個(gè)20 m×20 m樣地,采集并鑒定了樣地內(nèi)所有生境中的苔蘚物種,采用MTR、CCA和GLM等分析方法,探討了以下問題:(1)小秦嶺地區(qū)苔蘚植物組成,以及優(yōu)勢科和優(yōu)勢種;(2)綜合分析苔蘚植物之間及苔蘚植物與地形因子之間的關(guān)系;(3)苔蘚物種多樣性沿海拔梯度呈什么樣的分布格局?與喬木、灌木和草本的分布格局是否一致?

    1 材料與方法

    1.1 研究區(qū)域概況

    小秦嶺國家級(jí)自然保護(hù)區(qū)位于河南省靈寶市西部(34°23′—34°31′N,110°23′—110°44′E),保護(hù)區(qū)東西長 31 km,南北寬12 km,全區(qū)面積15160 hm2。森林覆蓋率81.2%,其主峰老鴉岔垴海拔2413.8 m,為河南省最高峰。本區(qū)位于暖溫帶大陸性氣候區(qū),屬季風(fēng)型半干旱氣候,年均氣溫11.2—14.2℃,極端最低溫度-17℃,極端最高溫度42.7℃,無霜期170—215 d,年降水量620 mm左右。

    1.2 樣地設(shè)置

    2012年,在保護(hù)區(qū)主峰老鴉岔垴,從海拔1020 m到2413 m設(shè)置樣地。海拔每上升50 m,選擇代表性的群落分別設(shè)置樣地,樣地大小為20 m×20 m,共計(jì)53個(gè)樣地。在樣地內(nèi),從地面到地上2 m范圍內(nèi),對(duì)樣地內(nèi)所有巖石、枯倒木、地面及木本植物個(gè)體上等所有苔蘚進(jìn)行采樣,樣地中同一棵樹或同一塊石頭的不同高度和方位的苔蘚都需要采樣,而同一點(diǎn)位置上同種苔蘚物種只采標(biāo)本一份,同一點(diǎn)位置不同種苔蘚都要采集,標(biāo)本存放于信封中并記錄樣方號(hào)和生境類型。采集過程中仔細(xì)檢查樣地每個(gè)角落,盡最大努力做到不遺漏,共收集到苔蘚標(biāo)本1378份。所有標(biāo)本帶回實(shí)驗(yàn)室通過顯微鏡鑒定到種,憑證標(biāo)本存放于河南農(nóng)業(yè)大學(xué)標(biāo)本館苔蘚標(biāo)本室。同時(shí)通過GPS測量每個(gè)樣地的海拔,地質(zhì)羅盤測量坡度和坡向。

    1.3 苔蘚植物物種多度的計(jì)算

    物種多度一般指的是群落中某物種個(gè)體數(shù)量。傳統(tǒng)的經(jīng)典取樣方法通常是在10 cm×10 cm、20 cm×20 cm或50 cm×50 cm等小尺度上進(jìn)行采樣[17- 20]。然后通過不同的方法計(jì)算物種個(gè)體在樣方內(nèi)的物種個(gè)體數(shù),從而得到多度。這種方式對(duì)于研究小尺度上苔蘚物種多樣性、苔蘚結(jié)皮及樹木與附生苔蘚之間的關(guān)系等方面效果較好。然而這種方法限制了在大尺度上苔蘚生態(tài)學(xué)的發(fā)展。

    為了有利于從宏觀的角度來研究苔蘚植物,研究參考了賈宏汝對(duì)苔蘚物種多度的處理方式[21],即某一苔蘚物種而言,只關(guān)注該物種在樣地中哪一棵樹、哪一塊石頭或者哪一斑塊地面出現(xiàn)了,不關(guān)注它在該棵樹、該塊石頭或者該斑塊地面上同種個(gè)體的數(shù)量信息,也就說該苔蘚物種在某個(gè)點(diǎn)上有或沒有,從而得到苔蘚的多度信息。

    1.4 物種多樣性指數(shù)計(jì)算

    計(jì)算每個(gè)樣地內(nèi)苔蘚植物的生物多樣性。

    (1)Patrick豐富度指數(shù)

    D=S

    (1)

    (2)β多樣性指數(shù)

    β=1-2c/(a+b)

    (2)

    式中,S為樣地中的物種數(shù)目,a和b是兩個(gè)樣地各自所單獨(dú)擁有的物種數(shù),c是兩樣地間共有的物種數(shù)。

    1.5 多元回歸樹

    本研究采用53個(gè)樣地海拔、坡向和坡度為自變量,240種苔蘚植物在53個(gè)樣地內(nèi)的多度為因變量進(jìn)行群叢分類。MRT采用R語言“mvpart”程序包進(jìn)行運(yùn)算[22]。

    1.6 排序分析

    CCA 分析中的物種矩陣選擇苔蘚物種在樣地中的多度組成(多度大于3的物種)。環(huán)境矩陣由海拔、坡向和坡度組成,分析過程中對(duì)數(shù)據(jù)進(jìn)行了標(biāo)準(zhǔn)化處理。通過R軟件中envfit函數(shù)檢驗(yàn)每個(gè)環(huán)境因子與苔蘚物種分布的顯著性,模擬次數(shù)選擇999次[23]。同時(shí)對(duì)CCA的最終結(jié)果進(jìn)行蒙特卡羅置換檢驗(yàn),分析模型是否達(dá)到顯著水平(P<0.05)。CCA分析及顯著性檢驗(yàn)通過R語言軟件中“Vegan”程序包實(shí)現(xiàn)[24]。

    1.7 廣義線性回歸

    苔蘚物種多樣性作為響應(yīng)變量,地形因子作為解釋變量。通過模型線性擬合來描述苔蘚物種多樣性隨地形因子的變化格局。線性模型和作圖通過R語言軟件中“car”程序包實(shí)現(xiàn)[25]。

    2 結(jié)果

    2.1 苔蘚物種多樣性

    采集的1378份苔蘚植物標(biāo)本經(jīng)過室內(nèi)鑒定,共有33科89屬240種。其中蘚類24科79屬226種,占總物種數(shù)的94.2%;苔類9科10屬14種,僅占總物種數(shù)的5.8%,豐富度較低。53個(gè)樣地苔蘚物種多樣性概況如表1所示。小秦嶺自然保護(hù)區(qū)苔蘚植物的優(yōu)勢科和多度排前10位的物種概況如表2所示。不同生境中的苔蘚物種多樣性差異如圖1所示。附表為小秦嶺自然保護(hù)區(qū)苔蘚植物物種名錄和區(qū)系分布(附表)。

    表2 優(yōu)勢科的屬、種統(tǒng)計(jì)及多度排前10位的物種

    2.2 分類結(jié)果

    MRT能兼顧誤差,又能保證樹的規(guī)模足夠小[8],以二岐式分割法將53個(gè)樣地分類,經(jīng)過交叉驗(yàn)證,結(jié)合宋永昌關(guān)于群叢的命名原則[26],認(rèn)為可以分為4類(圖2)。

    Ⅰ 長肋青蘚(Brachytheciumpopuleum)+大葉匐燈蘚(Plagiomniumsucculentum)+圓葉匐燈蘚(Plagiomniumvesicatum) 群落,包括1—35號(hào)樣地,該群落主要分布在海拔1020—1923 m。

    Ⅱ 灰白青蘚(Brachytheciumalbicans)+寬葉青蘚(Brachytheciumcurtum)+短肋羽蘚(Thuidiumkanedae)+平肋提燈蘚(Mniumlaevinerve) 群落,包括42—53號(hào)樣地,該群落主要分布在海拔1970—2080 m的地帶。

    Ⅲ 短肋羽蘚(Thuidiumkanedae)+光萼葉苔(Jungermannialeiantha)+薄羅蘚(Leskeapolycarpa)+叉肋蘚(Trachyphylluminflexum) 群落,包括36,37,38和40號(hào)樣地,該群落主要分布在海拔1970—2075 m。

    Ⅳ 寬葉青蘚(Brachytheciumcurtum)+垂蒴棉蘚(Plagiotheciumnemorale)+全緣匐燈蘚(Plagiomniumintegrum) 群落,包括39和41號(hào)樣地,該群落分布在海拔2020—2080 m。

    圖1 不同生境中苔蘚物種多樣性差異Fig.1 Bryophytes diversity between different habitats白色代表只在該生境中出現(xiàn)的物種數(shù),灰色代表在多個(gè)生境中出現(xiàn)的物種數(shù)

    圖2 小秦嶺自然保護(hù)區(qū)苔蘚植物MRT分類 Fig.2 Dendrogram of the multivariate regression trees classification of the bryophyte community in Xiaoqinling Nature Reserve

    2.3 排序結(jié)果

    對(duì)排序結(jié)果進(jìn)行蒙特卡羅檢驗(yàn),顯著性小于0.05,說明排序的結(jié)果可以接受地形因子對(duì)物種分布的解釋量,排序結(jié)果如圖3所示。每個(gè)環(huán)境因子的顯著性如表3所示。海拔、坡度和坡向三者的解釋量為7.29%,其中第一軸的特征值為21.85%,第二軸的特征值為15.83%。苔蘚在排序軸上的分布基本反映了植物空間分布特征隨地形因子的變化趨勢。

    表3 環(huán)境因子的顯著性檢驗(yàn)

    圖3 苔蘚物種分布與地形因子的關(guān)系Fig.3 Relationship between topography factors and bryophyte species distribution圖中圓點(diǎn)和三角形分別代表物種和樣地

    2.4 苔蘚物種多樣性沿海拔梯度分布格局

    苔蘚物種豐富度與海拔表現(xiàn)出極顯著相關(guān)(P<0.001),而與坡度和坡向未能達(dá)到顯著水平;苔蘚β多樣性與海拔、坡度和坡向均未達(dá)到顯著水平。苔蘚物種多樣性和β多樣性隨海拔的變化格局分別如圖4所示。苔蘚物種多樣性隨海拔的變化成明顯的上升格局,而苔蘚β多樣性隨海拔的變化呈水平變化。

    每種苔蘚物種(個(gè)體數(shù)>3)的海拔分布范圍如圖5所示。長肋青蘚、灰白青蘚、圓葉匐燈蘚和寬葉青蘚等物種的海拔分布極為廣泛,在海拔1020—2378 m都有分布,而叢本蘚(Anoectangiumaestivum)、狹葉小羽蘚(Haplocladiumangustifolium)和灰白青蘚(Brachytheciumalbicans)等物種海拔分布范圍較窄。

    圖4 苔蘚物種多樣性沿海拔分布格局Fig.4 distribution pattern of bryophyte species diversity along elevation gradients

    圖5 小秦嶺自然保護(hù)區(qū)苔蘚物種的海拔分布范圍Fig.5 The elevation range of bryophyte species distribution in Xiaoqinling nature reserve

    3 討論

    通過對(duì)小秦嶺地區(qū)苔蘚植物的調(diào)查研究表明,該保護(hù)區(qū)的苔蘚植物比較豐富。苔類植物總數(shù)遠(yuǎn)小于蘚類,這與汪岱華等的研究結(jié)果是一致的[27]。青蘚科(Bryaceae)、提燈蘚科(Mniaceae)、灰蘚科(Hypnaceae)和羽蘚科(Thuidiaceae),該4科的物種數(shù)占總物種數(shù)的47%,是小秦嶺國家自然保護(hù)區(qū)的主要類群。長肋青蘚(Brachytheciumpopuleum)、灰白青蘚(Brachytheciumalbicans)、圓葉匐燈蘚(Plagiomniumvesicatum)和寬葉青蘚(Brachytheciumcurtum)為小秦嶺自然保護(hù)區(qū)的常見種。

    本次研究結(jié)果顯示,不同的生境中苔蘚多樣性差異較大,在樣地中約有57.5%物種只傾向于分布某一類生境中,這就表明生境異質(zhì)性有助于豐富苔蘚物種多樣性,這也是在有些森林實(shí)施管理中要保持生境類型多樣性的主要原因。

    基于地形因子和物種組成數(shù)據(jù)的MRT技術(shù)將小秦嶺自然護(hù)區(qū)53塊樣地的240種苔蘚分成了4類,結(jié)果顯示4個(gè)群落都有各自的物種組成和群落特征。通過CCA分析表明海拔、坡度和坡向顯著影響著苔蘚物種分布。不同的海拔梯度上溫度、降水及空氣濕度等存在很大差異[28],不同坡度和坡向的土壤水分和養(yǎng)分含量也存在一定的差異[29-30]。因此,海拔、坡度和坡向3個(gè)地形因子顯著影響著苔蘚植物的分布。但三者的解釋量只有7.29%,解釋量較低,并且影響苔蘚物種分布的環(huán)境因子也較多,今后的研究應(yīng)逐步考慮氣候條件、樹皮的理化性質(zhì)及光照強(qiáng)度等更多的環(huán)境因子。

    本研究結(jié)果表明小秦嶺地區(qū)苔蘚物種多樣性隨海拔的上升呈顯著的上升格局。Stehn等在阿巴拉契亞山區(qū)研究也表明苔蘚物種多樣性隨海拔的上升呈上升格局[31],本研究結(jié)果與其是一致。陳云等在小秦嶺地區(qū)研究了喬木、灌木和草本隨海拔的變化格局,結(jié)果顯示喬木層物種豐富度呈單峰曲線,灌木層的豐富度在中海拔呈下降趨勢,草本層物種豐富度沿海拔梯度呈下降趨勢[32]。苔蘚物種與喬木、灌木和草本隨海拔的變化呈現(xiàn)出不一樣的格局,這就表明苔蘚物種生長所需要的環(huán)境條件和喬木、灌木、草本是有很大差異。海拔對(duì)苔蘚β多樣性的影響未能達(dá)到顯著水平,這可能是因?yàn)樘μ\個(gè)體小,受局部微生境影響較大。

    以上分析已表明海拔顯著影響著苔蘚植物的分布,那么每種苔蘚物種在海拔上是怎樣分布的?結(jié)果表明不同苔蘚的海拔生長范圍存在較大差異,特別是一些海拔分布范圍較窄的物種,說明這些物種對(duì)環(huán)境適應(yīng)能力差,因此苔蘚物種保護(hù)過程中應(yīng)多關(guān)注這些物種,并且要針對(duì)不同物種和不同海拔段實(shí)施不同的保護(hù)措施。

    總的來說,通過本研究了解了小秦嶺國家級(jí)自然保護(hù)區(qū)苔蘚物種組成、分類以及每種苔蘚植物的海拔分布范圍;苔蘚物種多樣性隨海拔的上升呈顯著的上升格局,與喬木、灌木和草本的變化格局不同。希望此研究能為苔蘚植物生態(tài)學(xué)研究和物種多樣性保護(hù)提供參考資料。然而,維管植物和苔蘚共存于一個(gè)生態(tài)系統(tǒng)中,本文未能考慮它們對(duì)苔蘚群落的影響,當(dāng)然,維管植物與苔蘚群落之間的關(guān)系也是非常復(fù)雜的,維管植物究竟怎樣影響著苔蘚群落?例如,喬木、灌木和草本分別是如何影響苔蘚物種分布的?不同生境中的苔蘚(樹生苔蘚、地面生苔蘚、石生苔蘚等)與維管植物之間存在怎樣的關(guān)系?這些問題都有待于深入研究。

    致謝:中國科學(xué)院植物研究所賈渝老師幫助苔蘚鑒定,特此致謝。

    [1] 雷波, 包維楷, 賈渝. 6種人工針葉幼林下地表苔蘚植物層片的物種多樣性與結(jié)構(gòu)特征. 植物生態(tài)學(xué)報(bào), 2004, 28(5): 594- 600.

    [2] 郭水良, 曹同. 長白山地區(qū)森林生態(tài)系統(tǒng)樹附生苔蘚植物群落分布格局研究. 植物生態(tài)學(xué)報(bào), 2000, 24(4): 442- 450.

    [3] 葉吉, 郝占慶, 于德永, 閆海冰, 封德全. 苔蘚植物生態(tài)功能的研究進(jìn)展. 應(yīng)用生態(tài)學(xué)報(bào), 2004, 15(10): 1939- 1942.

    [4] 姜炎彬, 邵小明. 苔蘚植物分布及其物種多樣性的研究評(píng)述. 武漢植物學(xué)研究, 2010, 28(3): 385- 390.

    [5] 吳征鎰. 中國植被. 北京: 科學(xué)出版社, 1980: 1- 1375.

    [6] 韓軍旺, 王進(jìn), 張旭, 賈宏汝, 袁志良. 小秦嶺國家級(jí)自然保護(hù)區(qū)植物群落數(shù)量分類與排序. 河南科學(xué), 2015, 33(4): 547- 552.

    [7] 張文靜, 張欽弟, 王晶, 馮飛, 畢潤成. 多元回歸樹與雙向指示種分析在群落分類中的應(yīng)用比較. 植物生態(tài)學(xué)報(bào), 2015, 39(6): 586- 592.

    [8] 賴江山, 米湘成, 任海保, 馬克平. 基于多元回歸樹的常綠闊葉林群叢數(shù)量分類-以古田山24公頃森林樣地為例. 植物生態(tài)學(xué)報(bào), 2010, 34(7): 761- 769.

    [9] De′ath G. Multivariate regression trees: a new technique for modeling species-environment relationships. Ecology, 2002, 83(4): 1105- 1117.

    [10] 張榮, 劉彤. 古爾班通古特沙漠南部植物多樣性及群落分類. 生態(tài)學(xué)報(bào), 2012, 32(19): 6056- 6066.

    [11] 黃永濤, 姚蘭, 艾訓(xùn)儒, 呂世安, 丁易. 鄂西南兩個(gè)自然保護(hù)區(qū)亞熱帶常綠落葉闊葉混交林類型及其常綠和落葉物種組成結(jié)構(gòu)分析. 植物生態(tài)學(xué)報(bào), 2015, 39(10): 990- 1002.

    [12] 張?jiān)? 曹同, 潘伯榮. 新疆博格達(dá)山地面生苔蘚植被的數(shù)量分類與排序研究. 植物生態(tài)學(xué)報(bào), 2002, 26(1): 10- 16.

    [13] Pescador D S, de Bello F, Valladares F, Escudero A. Plant trait variation along an altitudinal gradient in Mediterranean high mountain grasslands: controlling the species turnover effect. PLoS One, 2015, 10(3): e0118876.

    [14] 朱源, 康慕誼, 江源, 劉全儒. 賀蘭山木本植物群落物種多樣性的海拔格局. 植物生態(tài)學(xué)報(bào), 2008, 32(3): 574- 581.

    [15] Koenig I, Feldmeyer-Christe E, Mitchell E A D. Comparative ecology of vascular plant, bryophyte and testate amoeba communities in fourSphagnumpeatlands along an altitudinal gradient in Switzerland. Ecological Indicators, 2015, 54: 48- 59.

    [16] He J S, Chen W L. A review of gradient changes in species diversity of land plant communities. Acta Ecologica Sinica, 1997, 17(1): 91- 99.

    [17] 藺菲, 郝占慶, 李步杭, 葉吉, 戴冠華, 張健, 倪偉東. 長白山暗針葉林苔蘚植物群落特征與林木更新的關(guān)系. 生態(tài)學(xué)報(bào), 2007, 27(4): 1308- 1314.

    [18] 劉蔚秋, 戴小華, 王永繁, 雷純義. 影響廣東黑石頂樹附生苔蘚分布的環(huán)境因子. 生態(tài)學(xué)報(bào), 2008, 28(3): 1080- 1088.

    [19] 王中生, 安樹青, 方炎明. 苔蘚植物生殖生態(tài)學(xué)研究. 生態(tài)學(xué)報(bào), 2003, 23(11): 2444- 2452.

    [20] 福英, 白學(xué)良, 張樂, 畢庚辰, 馮超, 寇瑾, 薩如拉. 五大連池火山熔巖地貌苔蘚植物對(duì)土壤養(yǎng)分積累的作用. 生態(tài)學(xué)報(bào), 2015, 35(10): 3288- 3297.

    [21] 賈宏汝. 環(huán)境空間異質(zhì)性對(duì)暖溫帶-北亞熱帶過渡區(qū)植物群落構(gòu)建的影響[D]. 鄭州: 鄭州大學(xué), 2015.

    [22] De′Ath G. mvpart: multivariate partitioning. R package version 1.2- 4. 2006.

    [24] MHH W H. Vegan: Community Ecology Package.Rpackageversion, 2-0. 2011.

    [25] Fox J, Weisberg S. Using car functions in other functions. 2013.

    [26] 宋永昌. 植被生態(tài)學(xué). 上海: 華東師范大學(xué)出版社, 2001.

    [27] 汪岱華, 王幼芳, 左勤, 李敏, 魏倩倩, 李曉芹, 趙明水. 浙江西天目山30年來苔蘚物種多樣性變化. 生物多樣性, 2013, 21(2): 170- 176.

    [28] Lieberman D, Lieberman M, Peralta R, Hartshorn G S. Tropical forest structure and composition on a large-scale altitudinal gradient in Costa Rica. Journal of Ecology, 1996, 84(2): 137- 152.

    [29] Auslander M, Nevo E, Inbar M. The effects of slope orientation on plant growth, developmental instability and susceptibility to herbivores. Journal of Arid Environments, 2003, 55(3): 405- 416.

    [30] Ferrer-Castán D, Vetaas O R. Pteridophyte richness, climate and topography in the Iberian Peninsula: comparing spatial and nonspatial models of richness patterns. Global Ecology and Biogeography, 2005, 14(2): 155- 165.

    [31] Stehn S E, Webster C R, Glime J M, Jenkins M A. Elevational gradients of bryophyte diversity, life forms, and community assemblage in the southern Appalachian Mountains. Canadian Journal of Forest Research, 2010, 40(11): 2164- 2174.

    [32] 陳云, 王海亮, 韓軍旺, 韋博良, 賈宏汝, 葉永忠, 袁志良. 小秦嶺森林群落數(shù)量分類、排序及多樣性垂直格局. 生態(tài)學(xué)報(bào), 2014, 34(8): 2068- 2075.

    Numerical classification, ordination, and species diversity of bryophytes along elevation gradients in Xiaoqinling

    CHEN Yun1, FENG Jiawei2, NIU Shuai2, XU Ning2, HAN Junwang3, YE Yongzhong2, YUAN Zhiliang2,*

    1CollegeofForestry,HenanAgriculturalUniversit,Zhengzhou450002,China2CollegeofLifeSciences,HenanAgriculturalUniversit,Zhengzhou450002,China3AdministrationBureauofXiaoqinlingNationalNatureReserveofHenan,Lingbao472500,China

    Bryophytes are widely distributed and constitute an important part in vegetation diversity. More than 20,000 species of bryophytes are found worldwide. Bryophytes play important roles in the ecosystem, particularly in water balance and nutrient accumulation. However, bryophytes are smaller than vascular plants and are often ignored in ecological studies. Spatial patterns of species diversity along the elevation gradient are of great significance in illuminating community ecology. What is the distribution pattern of bryophyte species diversity along an elevation gradient? We do not find the related research in China. The Xiaoqinling National Nature Reserve is located in Lingbao (34°23′—34°31′N, 110°23′—110°44′E, Henan Province, China). This nature reserve covers 31 km from east to west and 12 km from north to south, and spreads over a total area of 15,160 km2. The reserve harbors the highest peak in Henan, with an elevation of 2413.8 m. The area has rich biological resources with a diverse and complex vegetation composition. We constructed 53 plots in a 400m2area (20 m × 20 m) along the altitudinal gradient in Xiaoqinling, collected all bryophyte species encountered in these 53 plots, and recorded habitat types such as ground, tree trunks, roots, and rotten wood. Investigators carefully checked every corner of each plot to collect all bryophyte species growing on the plot. A total of 1378 specimens were collected. All specimens were transported back to the laboratory and identified. The bryophyte community was classified using multiple regression trees. Ordination analysis was conducted by canonical correspondence analysis, and the generalized linear model was used to study the distribution pattern of species diversity along elevation gradients. Results were as follows. (1) The collected 1378 bryophyte specimens belonged to 33 families, 89 genera, and 240 species, including 226 species of mosses and 14 species of liverworts. (2) The bryophyte community was divided into four categories by cross validation: I:Brachytheciumpopuleum+Plagiomniumsucculentum+Plagiomniumvesicatum; II:Brachytheciumalbicans+Brachytheciumcurtum+Thuidiumkanedae+Mniumlaevinerve; III:Thuidiumkanedae+Jungermannialeiantha+Leskeapolycarpa+Trachyphylluminflexum; IV:Brachytheciumcurtum+Plagiotheciumnemorale+Plagiomniumintegrum. (3) The interpretation of the elevation, slope, and aspect to the distribution of the bryophyte species was 7.29%. (4) The bryophyte species diversity increased with the increase in altitude. In this study, we examined the species composition, classification, and distribution range along elevation gradients of the bryophyte species in Xiaoqinling. The diversity of the bryophyte species along the altitude showed an increasing tendency and differed from the distribution patterns of trees, shrubs, and herbs. Thus, this study provides a reference for research on bryophyte species coexistence and protection of species diversity. Vascular plants and bryophytes coexist in an ecosystem and the relationship between vascular plants and the bryophyte community is complex. However, we did not consider the effect of vascular plants on the bryophyte community and these factors need further analysis. Environmental factors, such as climatic conditions, physical and chemical properties of bark, and light intensity, must also be considered in future studies.

    MRT; CCA; elevation; species diversity; bryophyte

    附表 小秦嶺自然保護(hù)區(qū)苔蘚植物物種名錄和區(qū)系分布

    河南省科技攻關(guān)(132102110133);河南省教育廳科學(xué)技術(shù)研究重點(diǎn)項(xiàng)目(14A180013)

    2015- 12- 24; 網(wǎng)絡(luò)出版日期:2016- 10- 19

    10.5846/stxb201512242568

    *通訊作者Corresponding author.E-mail: zlyuan163@163.com

    陳云,馮佳偉,牛帥,許寧,韓軍旺,葉永忠,袁志良.小秦嶺自然保護(hù)區(qū)苔蘚植物群落數(shù)量分類、排序及多樣性垂直格局.生態(tài)學(xué)報(bào),2017,37(8):2653- 2664.

    Chen Y, Feng J W, Niu S, Xu N, Han J W, Ye Y Z, Yuan Z L.Numerical classification, ordination, and species diversity of bryophytes along elevation gradients in Xiaoqinling.Acta Ecologica Sinica,2017,37(8):2653- 2664.

    猜你喜歡
    苔蘚秦嶺樣地
    The mysterious “glacier mice”
    暑期秦嶺游
    洞穿秦嶺
    “小矮人”苔蘚的大用途
    軍事文摘(2021年18期)2021-12-02 01:28:36
    額爾古納市興安落葉松中齡林植被碳儲(chǔ)量研究
    昆明市主要綠化樹種閾值測定與分析
    基于角尺度模型的林業(yè)樣地空間結(jié)構(gòu)分析
    好忙好忙的秦嶺
    苔蘚能“對(duì)付”空氣污染嗎?
    黏液水腫性苔蘚1例
    中国美女看黄片| 成人午夜精彩视频在线观看| 成人av在线播放网站| 成年女人永久免费观看视频| 黄片wwwwww| 国内少妇人妻偷人精品xxx网站| 欧美激情久久久久久爽电影| 蜜桃亚洲精品一区二区三区| 少妇熟女aⅴ在线视频| 亚洲美女搞黄在线观看| 欧美丝袜亚洲另类| 亚洲性久久影院| 热99re8久久精品国产| 男女视频在线观看网站免费| 淫秽高清视频在线观看| 亚洲熟妇中文字幕五十中出| 99久久精品热视频| 亚洲欧美精品专区久久| 12—13女人毛片做爰片一| 岛国毛片在线播放| 国产亚洲精品av在线| 成人综合一区亚洲| 亚洲图色成人| 女人被狂操c到高潮| 国产片特级美女逼逼视频| 日韩av在线大香蕉| 青春草亚洲视频在线观看| 久久中文看片网| 久久婷婷人人爽人人干人人爱| 中文欧美无线码| 日韩av不卡免费在线播放| 三级男女做爰猛烈吃奶摸视频| 伦精品一区二区三区| 国产av麻豆久久久久久久| 久久久a久久爽久久v久久| 校园春色视频在线观看| 久久99热6这里只有精品| 一本一本综合久久| 一本久久精品| 亚洲av熟女| or卡值多少钱| 日韩强制内射视频| 国产亚洲91精品色在线| eeuss影院久久| 日日啪夜夜撸| 成人毛片a级毛片在线播放| 欧美一区二区精品小视频在线| 在线a可以看的网站| 岛国毛片在线播放| 伦理电影大哥的女人| 中文字幕av在线有码专区| 精品久久久久久久久久免费视频| 欧美3d第一页| 精品无人区乱码1区二区| 日本欧美国产在线视频| 男女视频在线观看网站免费| 精品免费久久久久久久清纯| 欧美性猛交╳xxx乱大交人| 国产精品三级大全| 爱豆传媒免费全集在线观看| 少妇被粗大猛烈的视频| 国产一区二区在线观看日韩| 国产真实伦视频高清在线观看| 亚洲成人久久性| 久久精品91蜜桃| 国产亚洲精品久久久久久毛片| 黄色日韩在线| .国产精品久久| 一个人看视频在线观看www免费| 亚洲成av人片在线播放无| 亚洲欧美日韩高清在线视频| 观看免费一级毛片| 真实男女啪啪啪动态图| 精品久久久噜噜| 亚洲第一电影网av| 真实男女啪啪啪动态图| 亚洲av成人精品一区久久| av在线老鸭窝| 国产色爽女视频免费观看| 日本三级黄在线观看| 午夜激情福利司机影院| 中国美女看黄片| 五月玫瑰六月丁香| 久久久成人免费电影| 中文欧美无线码| 国产午夜精品久久久久久一区二区三区| 九草在线视频观看| 成人高潮视频无遮挡免费网站| 午夜免费男女啪啪视频观看| 国产在线男女| 欧美日本视频| 国产精品不卡视频一区二区| 国产精品99久久久久久久久| 亚洲最大成人中文| 熟妇人妻久久中文字幕3abv| 免费无遮挡裸体视频| 一级黄色大片毛片| 一区二区三区高清视频在线| 男女啪啪激烈高潮av片| 久久欧美精品欧美久久欧美| 午夜a级毛片| 黄片wwwwww| 国产精品免费一区二区三区在线| 99在线人妻在线中文字幕| 色哟哟·www| av在线天堂中文字幕| 日日干狠狠操夜夜爽| 欧美zozozo另类| 欧美zozozo另类| 久久久色成人| 国产中年淑女户外野战色| 日本黄色片子视频| 九九爱精品视频在线观看| 国产精品麻豆人妻色哟哟久久 | 国产精品99久久久久久久久| 在线观看66精品国产| 午夜福利在线观看免费完整高清在 | a级一级毛片免费在线观看| 久久久欧美国产精品| 亚洲欧美成人精品一区二区| 久久久久国产网址| 高清午夜精品一区二区三区 | 天天躁日日操中文字幕| 久久久国产成人免费| 久久久国产成人精品二区| 成年av动漫网址| 国产欧美日韩精品一区二区| 国产探花极品一区二区| 亚洲成人av在线免费| 国产成人精品久久久久久| 少妇高潮的动态图| 日韩av不卡免费在线播放| 亚洲乱码一区二区免费版| 亚洲美女视频黄频| 天堂√8在线中文| 综合色av麻豆| 免费人成视频x8x8入口观看| 麻豆成人午夜福利视频| 男人狂女人下面高潮的视频| 久久精品国产鲁丝片午夜精品| 中文字幕av在线有码专区| av在线蜜桃| 男人舔女人下体高潮全视频| 欧美精品国产亚洲| а√天堂www在线а√下载| 1024手机看黄色片| 国产69精品久久久久777片| 久久久a久久爽久久v久久| 欧美日本视频| 夫妻性生交免费视频一级片| av在线蜜桃| 亚洲欧美精品专区久久| 日韩欧美一区二区三区在线观看| 国产成人福利小说| 女人十人毛片免费观看3o分钟| 亚洲国产色片| 亚洲精品456在线播放app| 久久人妻av系列| 日韩视频在线欧美| 国产探花极品一区二区| 成人毛片60女人毛片免费| 亚洲成人中文字幕在线播放| 色吧在线观看| 九九热线精品视视频播放| 日本免费a在线| 一本久久精品| 免费看日本二区| 成人美女网站在线观看视频| 国产探花极品一区二区| 三级国产精品欧美在线观看| 国产精品福利在线免费观看| 丰满的人妻完整版| 丰满乱子伦码专区| 天天躁日日操中文字幕| 99热6这里只有精品| 久久久久久久久久久丰满| 欧美高清成人免费视频www| 久久久久久久久久久免费av| 国内少妇人妻偷人精品xxx网站| 3wmmmm亚洲av在线观看| 国产成人福利小说| 亚洲五月天丁香| 少妇人妻精品综合一区二区 | 超碰av人人做人人爽久久| 成人午夜高清在线视频| 三级毛片av免费| 色尼玛亚洲综合影院| 麻豆成人av视频| 亚洲成人久久性| 免费电影在线观看免费观看| 床上黄色一级片| 一本一本综合久久| 午夜精品国产一区二区电影 | 18禁裸乳无遮挡免费网站照片| a级毛色黄片| 国产综合懂色| 麻豆成人av视频| 两个人的视频大全免费| 亚洲av男天堂| 亚洲av熟女| 啦啦啦啦在线视频资源| 免费观看精品视频网站| 色5月婷婷丁香| 欧美一区二区精品小视频在线| 免费av毛片视频| 青青草视频在线视频观看| 国产淫片久久久久久久久| 男女那种视频在线观看| 成人高潮视频无遮挡免费网站| 精品人妻偷拍中文字幕| 国产私拍福利视频在线观看| 男女视频在线观看网站免费| 国产69精品久久久久777片| 可以在线观看毛片的网站| 国产精品一区二区在线观看99 | 精品免费久久久久久久清纯| 51国产日韩欧美| 精品久久久久久久久久免费视频| 国产高清不卡午夜福利| 亚洲精品日韩av片在线观看| 九九热线精品视视频播放| 久久国内精品自在自线图片| 中文字幕制服av| 午夜免费男女啪啪视频观看| 国产一级毛片七仙女欲春2| 国产精品久久久久久精品电影| 美女xxoo啪啪120秒动态图| 久久久国产成人免费| 午夜久久久久精精品| 99久久九九国产精品国产免费| 国产午夜精品久久久久久一区二区三区| 色综合站精品国产| 国产真实乱freesex| 亚洲在线自拍视频| 免费av不卡在线播放| 97人妻精品一区二区三区麻豆| 噜噜噜噜噜久久久久久91| 久久久久网色| 日本在线视频免费播放| 成人永久免费在线观看视频| 日韩av不卡免费在线播放| 国产成年人精品一区二区| 精品久久久久久久末码| 一级av片app| 中文字幕av在线有码专区| 亚洲中文字幕一区二区三区有码在线看| 亚洲美女视频黄频| 日本免费a在线| 韩国av在线不卡| 看片在线看免费视频| 精品熟女少妇av免费看| 男女下面进入的视频免费午夜| 99九九线精品视频在线观看视频| 精品一区二区三区视频在线| 在线播放国产精品三级| 黄色配什么色好看| 神马国产精品三级电影在线观看| 亚洲欧洲日产国产| 国产精华一区二区三区| 色尼玛亚洲综合影院| 九草在线视频观看| 久久国产乱子免费精品| 亚洲内射少妇av| 最后的刺客免费高清国语| 午夜精品国产一区二区电影 | 国产一级毛片在线| 97超碰精品成人国产| 1000部很黄的大片| 97超视频在线观看视频| 一区二区三区免费毛片| 亚洲美女搞黄在线观看| 日本熟妇午夜| 亚洲成人精品中文字幕电影| 国产一区二区在线av高清观看| 毛片女人毛片| av在线老鸭窝| 亚洲人成网站在线播放欧美日韩| 久久亚洲精品不卡| 亚洲国产精品成人综合色| 国产精品伦人一区二区| 中国美白少妇内射xxxbb| 日韩av不卡免费在线播放| 久久久久久国产a免费观看| av女优亚洲男人天堂| 中文在线观看免费www的网站| 国产精品1区2区在线观看.| 3wmmmm亚洲av在线观看| 观看美女的网站| 性欧美人与动物交配| av在线亚洲专区| 国产真实伦视频高清在线观看| 亚洲无线在线观看| 日本黄色片子视频| 日韩欧美在线乱码| 大型黄色视频在线免费观看| 床上黄色一级片| 简卡轻食公司| 国产在线男女| 国产精品久久久久久亚洲av鲁大| 有码 亚洲区| 国产蜜桃级精品一区二区三区| 搡老妇女老女人老熟妇| 91精品国产九色| www.av在线官网国产| 看非洲黑人一级黄片| 最近中文字幕高清免费大全6| 人妻久久中文字幕网| av天堂中文字幕网| 美女国产视频在线观看| 国产成人a∨麻豆精品| 亚洲图色成人| 成人亚洲欧美一区二区av| 91麻豆精品激情在线观看国产| 亚洲av免费高清在线观看| 成人一区二区视频在线观看| 一进一出抽搐gif免费好疼| 国国产精品蜜臀av免费| 欧美色视频一区免费| 日韩一区二区视频免费看| 精品一区二区三区人妻视频| 国产黄a三级三级三级人| 麻豆国产av国片精品| a级一级毛片免费在线观看| 日本黄大片高清| 国产视频首页在线观看| 桃色一区二区三区在线观看| 内射极品少妇av片p| 日韩欧美精品v在线| 国产黄色小视频在线观看| 欧美极品一区二区三区四区| 日韩,欧美,国产一区二区三区 | 亚洲精品日韩在线中文字幕 | 99久久精品一区二区三区| 欧美又色又爽又黄视频| 18禁在线播放成人免费| a级毛色黄片| 日韩,欧美,国产一区二区三区 | 色噜噜av男人的天堂激情| 69av精品久久久久久| 日本在线视频免费播放| 欧美日韩一区二区视频在线观看视频在线 | 国产精品一二三区在线看| av在线亚洲专区| eeuss影院久久| 嫩草影院入口| 男女那种视频在线观看| 永久网站在线| 久久久久久国产a免费观看| 久久精品国产亚洲网站| 国产成人一区二区在线| 婷婷色综合大香蕉| 一区二区三区免费毛片| 丝袜喷水一区| 久久久久久九九精品二区国产| 一级毛片电影观看 | 久久人人爽人人爽人人片va| 人妻少妇偷人精品九色| 国产精品人妻久久久久久| www日本黄色视频网| 国产成人影院久久av| 国产欧美日韩精品一区二区| 日产精品乱码卡一卡2卡三| 国产视频首页在线观看| 国产一级毛片在线| 国产高清不卡午夜福利| 国产精品国产三级国产av玫瑰| 欧美xxxx性猛交bbbb| 国产探花极品一区二区| 婷婷六月久久综合丁香| 亚洲精品国产av成人精品| 99精品在免费线老司机午夜| 久久99精品国语久久久| 亚洲欧美日韩卡通动漫| 国产大屁股一区二区在线视频| 1000部很黄的大片| 听说在线观看完整版免费高清| 成熟少妇高潮喷水视频| 国产黄片美女视频| 精品午夜福利在线看| 一级黄片播放器| 一个人看的www免费观看视频| 丰满的人妻完整版| 老司机影院成人| 91精品一卡2卡3卡4卡| 搡女人真爽免费视频火全软件| 亚洲四区av| 欧美人与善性xxx| 丰满人妻一区二区三区视频av| 久久久国产成人免费| 亚洲,欧美,日韩| 2021天堂中文幕一二区在线观| 亚洲av成人av| 欧美激情在线99| 毛片女人毛片| 卡戴珊不雅视频在线播放| 国产成人影院久久av| 欧美成人a在线观看| 美女xxoo啪啪120秒动态图| 美女内射精品一级片tv| 久久草成人影院| 亚洲国产高清在线一区二区三| 综合色av麻豆| 国产v大片淫在线免费观看| 日韩国内少妇激情av| 国产精品,欧美在线| 欧美xxxx性猛交bbbb| 国产精品乱码一区二三区的特点| 九九久久精品国产亚洲av麻豆| 亚洲国产精品成人久久小说 | 久久精品影院6| 日韩一区二区三区影片| 亚洲成人av在线免费| 综合色av麻豆| 蜜桃亚洲精品一区二区三区| 日本-黄色视频高清免费观看| 久久久久久大精品| 中文字幕制服av| 在线观看66精品国产| 国国产精品蜜臀av免费| 精品日产1卡2卡| 欧洲精品卡2卡3卡4卡5卡区| 在线a可以看的网站| 久久精品综合一区二区三区| 免费看日本二区| 美女脱内裤让男人舔精品视频 | 成年版毛片免费区| 日韩欧美一区二区三区在线观看| av女优亚洲男人天堂| 国语自产精品视频在线第100页| 麻豆国产97在线/欧美| 欧美成人免费av一区二区三区| 国产午夜福利久久久久久| 日日摸夜夜添夜夜爱| 可以在线观看毛片的网站| 亚洲精品自拍成人| 青春草亚洲视频在线观看| 国产91av在线免费观看| 欧美精品一区二区大全| 综合色丁香网| 十八禁国产超污无遮挡网站| 国产av在哪里看| 69av精品久久久久久| 午夜精品一区二区三区免费看| 国产精品一区二区性色av| 欧美性感艳星| 一本精品99久久精品77| 亚洲精品久久久久久婷婷小说 | 身体一侧抽搐| 日本免费a在线| 欧美最黄视频在线播放免费| 美女大奶头视频| 久99久视频精品免费| 久久久久国产网址| 亚洲18禁久久av| 亚洲av电影不卡..在线观看| 天堂中文最新版在线下载 | 男女下面进入的视频免费午夜| 亚洲成人久久爱视频| 国产免费男女视频| 五月玫瑰六月丁香| 一级av片app| 久久热精品热| 国产成人影院久久av| 亚洲性久久影院| 一边亲一边摸免费视频| 亚洲精品456在线播放app| 国产免费一级a男人的天堂| 嫩草影院新地址| 国产av麻豆久久久久久久| 男女下面进入的视频免费午夜| 美女内射精品一级片tv| 最近的中文字幕免费完整| 国产精品国产高清国产av| 亚洲国产欧美在线一区| 国产久久久一区二区三区| 精品人妻偷拍中文字幕| 国产成人影院久久av| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久久久久黄片| 午夜福利成人在线免费观看| 久久欧美精品欧美久久欧美| 男插女下体视频免费在线播放| 日韩成人伦理影院| 丰满人妻一区二区三区视频av| 青春草国产在线视频 | 搡女人真爽免费视频火全软件| 成年免费大片在线观看| 边亲边吃奶的免费视频| 男女视频在线观看网站免费| kizo精华| 亚州av有码| 午夜福利成人在线免费观看| 日韩一本色道免费dvd| 午夜老司机福利剧场| 亚洲在久久综合| 色哟哟·www| 中文在线观看免费www的网站| 精品久久久久久久久亚洲| 中文字幕精品亚洲无线码一区| 亚洲中文字幕日韩| 亚洲av熟女| 欧美激情国产日韩精品一区| 有码 亚洲区| 中文在线观看免费www的网站| 亚洲自拍偷在线| 久久久久性生活片| 国内少妇人妻偷人精品xxx网站| 国产成人a∨麻豆精品| 国产精品乱码一区二三区的特点| 99九九线精品视频在线观看视频| 99热精品在线国产| av免费在线看不卡| 中文字幕免费在线视频6| 神马国产精品三级电影在线观看| 在线观看午夜福利视频| 国产精品乱码一区二三区的特点| 国产成人福利小说| 日本三级黄在线观看| 五月伊人婷婷丁香| 人妻系列 视频| 中文字幕av成人在线电影| 日本欧美国产在线视频| 国产精品国产高清国产av| 免费看美女性在线毛片视频| 国内久久婷婷六月综合欲色啪| 美女大奶头视频| 久久精品夜色国产| 日本黄色视频三级网站网址| 亚洲五月天丁香| 伦理电影大哥的女人| 一级二级三级毛片免费看| 一级毛片aaaaaa免费看小| 日韩三级伦理在线观看| 国产大屁股一区二区在线视频| 婷婷六月久久综合丁香| 99热这里只有精品一区| 91aial.com中文字幕在线观看| 99热6这里只有精品| 日韩欧美 国产精品| 能在线免费看毛片的网站| 九草在线视频观看| 超碰av人人做人人爽久久| 日韩在线高清观看一区二区三区| 又粗又硬又长又爽又黄的视频 | kizo精华| 九色成人免费人妻av| 日韩大尺度精品在线看网址| 亚洲一级一片aⅴ在线观看| 午夜久久久久精精品| 18禁黄网站禁片免费观看直播| 桃色一区二区三区在线观看| 亚洲精品自拍成人| 亚洲国产色片| 麻豆久久精品国产亚洲av| 国产亚洲精品久久久com| 美女脱内裤让男人舔精品视频 | 日韩 亚洲 欧美在线| 12—13女人毛片做爰片一| 日韩亚洲欧美综合| 熟女电影av网| 六月丁香七月| 九草在线视频观看| 免费搜索国产男女视频| 久久久久久久久久久免费av| av视频在线观看入口| 天天一区二区日本电影三级| 天美传媒精品一区二区| 日本熟妇午夜| 69av精品久久久久久| 大香蕉久久网| 亚洲婷婷狠狠爱综合网| 99久久九九国产精品国产免费| 亚洲aⅴ乱码一区二区在线播放| 伊人久久精品亚洲午夜| 少妇裸体淫交视频免费看高清| 国产成人一区二区在线| 午夜亚洲福利在线播放| 日韩国内少妇激情av| 亚洲av熟女| 看十八女毛片水多多多| 一边摸一边抽搐一进一小说| 中文亚洲av片在线观看爽| 免费观看的影片在线观看| 亚洲色图av天堂| 亚洲五月天丁香| 成人特级黄色片久久久久久久| 中国美女看黄片| 成年女人看的毛片在线观看| 亚洲欧美日韩高清专用| 成人毛片a级毛片在线播放| 久99久视频精品免费| 亚洲精品自拍成人| 日本黄色片子视频| av免费观看日本| 校园春色视频在线观看| 男人舔女人下体高潮全视频| 欧洲精品卡2卡3卡4卡5卡区| 国产成人91sexporn| 尤物成人国产欧美一区二区三区| 99热只有精品国产| 日韩一区二区视频免费看| 舔av片在线| 99久久九九国产精品国产免费| 午夜老司机福利剧场| 一区二区三区高清视频在线| 欧美色欧美亚洲另类二区| 哪个播放器可以免费观看大片| 如何舔出高潮| 丰满的人妻完整版| 啦啦啦观看免费观看视频高清| 久久久久久大精品| 日本熟妇午夜| 亚洲高清免费不卡视频| 熟女电影av网| 乱系列少妇在线播放|