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

    木論喀斯特森林木本植物多樣性垂直格局

    2016-02-07 00:44:20蘭斯安曾馥平宋同清彭晚霞覃文更
    生態(tài)學(xué)報(bào) 2016年22期
    關(guān)鍵詞:喀斯特樣方樣地

    蘭斯安, 宋 敏,3, 曾馥平, 宋同清,*, 彭晚霞, 杜 虎, 韓 暢,4, 陳 莉,4, 覃文更

    1 中國(guó)科學(xué)院亞熱帶農(nóng)業(yè)生態(tài)研究所亞熱帶農(nóng)業(yè)生態(tài)過(guò)程重點(diǎn)實(shí)驗(yàn)室, 長(zhǎng)沙 410125 2 中國(guó)科學(xué)院環(huán)江喀斯特生態(tài)系統(tǒng)觀測(cè)研究站, 環(huán)江 547100 3 湖南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院, 長(zhǎng)沙 410128 4 湖南農(nóng)業(yè)大學(xué)生物科學(xué)技術(shù)學(xué)院, 長(zhǎng)沙 410128 5 廣西壯族自治區(qū)木論國(guó)家級(jí)自然保護(hù)區(qū)管理局,環(huán)江 547200

    木論喀斯特森林木本植物多樣性垂直格局

    蘭斯安1,2, 宋 敏1,2,3, 曾馥平1,2, 宋同清1,2,*, 彭晚霞1,2, 杜 虎1,2, 韓 暢1,2,4, 陳 莉1,2,4, 覃文更5

    1 中國(guó)科學(xué)院亞熱帶農(nóng)業(yè)生態(tài)研究所亞熱帶農(nóng)業(yè)生態(tài)過(guò)程重點(diǎn)實(shí)驗(yàn)室, 長(zhǎng)沙 410125 2 中國(guó)科學(xué)院環(huán)江喀斯特生態(tài)系統(tǒng)觀測(cè)研究站, 環(huán)江 547100 3 湖南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院, 長(zhǎng)沙 410128 4 湖南農(nóng)業(yè)大學(xué)生物科學(xué)技術(shù)學(xué)院, 長(zhǎng)沙 410128 5 廣西壯族自治區(qū)木論國(guó)家級(jí)自然保護(hù)區(qū)管理局,環(huán)江 547200

    基于木論國(guó)家級(jí)自然保護(hù)區(qū)內(nèi)西南坡一條垂直樣帶 16個(gè)樣方的調(diào)查資料,運(yùn)用豐富度指數(shù)、物種多樣性指數(shù)和均勻度指數(shù)對(duì)木論喀斯特森林木本植物多樣性的垂直格局進(jìn)行研究。結(jié)果表明:(1) 在16個(gè)樣方中共記錄植株5089個(gè),分屬48科90屬120種。(2) 喀斯特森林木本植物群落具有明顯的垂直分布格局,采用等級(jí)聚類(lèi)和非度量多維標(biāo)度排序方法可將植物群落沿海拔梯度劃分為 3個(gè)類(lèi)型,分別是:下坡位以小果厚殼桂(Cryptocaryaaustrokweichouensis)和灰?guī)r棒柄花(Cleidionbracteosum) 為主要優(yōu)勢(shì)種的群落;中坡位以長(zhǎng)管越南茜(Rubovietnamiaaristata)和羅傘(Brassaiopsis glomerulata)為主要優(yōu)勢(shì)種的群落;上坡位以化香樹(shù)(Platycaryastrobilacea)和密花樹(shù)(Rapaneaneriifolia)為主要優(yōu)勢(shì)種的群落。沿海拔梯度升高,水熱組合發(fā)生變化,地帶性植被以及群落內(nèi)的建群種和優(yōu)勢(shì)種都發(fā)生了變化。(3) 物種生活型方面,常綠樹(shù)種所占的比例高于落葉闊葉樹(shù)種,針葉樹(shù)種僅出現(xiàn)在山頂。(4) 沿海拔梯度物種豐富度及α多樣性指數(shù)中的 Shannon-Wiener 指數(shù)曲線(xiàn)均符合單峰格局,峰值出現(xiàn)在坡腰位置。(5)β多樣性指數(shù)起伏較大,總體上β多樣性隨著海拔的升高呈遞減的趨勢(shì)。

    垂直格局;α多樣性;β多樣性;木本植物;喀斯特

    植物多樣性的空間分布格局是植被生態(tài)學(xué)研究的核心,可以進(jìn)一步揭示生物多樣性的維持機(jī)制。例如,山地植物群落的物種多樣性如何隨海拔梯度的變化一直是生態(tài)學(xué)熱點(diǎn)研究領(lǐng)域之一[1- 3]。而且由于環(huán)境因子沿海拔梯度的變化要比沿緯度梯度快1000倍[4]。沿環(huán)境梯度不同生活型的物種多樣性差異較大,木本植物具有比草本植物更為明確的分布格局[5]。因此,研究木本植物多樣性的海拔梯度格局對(duì)于揭示生物多樣性隨環(huán)境梯度的變化規(guī)律具有重要意義[6- 7]。

    桂西北喀斯特山地既是廣西生物多樣性研究熱點(diǎn)區(qū)域之一,也是我國(guó)生物多樣性3個(gè)特有分布中心之一[8- 9]。木論國(guó)家自然保護(hù)區(qū)位于廣西壯族自治區(qū)環(huán)江毛南族自治縣西北部,與貴州省茂蘭國(guó)家自然保護(hù)區(qū)連片構(gòu)成當(dāng)今世界上分布面積最廣、保護(hù)最好、原生性最強(qiáng)的喀斯特非地帶性植被區(qū),受地球內(nèi)動(dòng)力、強(qiáng)烈的地質(zhì)運(yùn)動(dòng)、高溫多雨且分布不均、碳酸鹽巖溶蝕性強(qiáng)、水文二維結(jié)構(gòu)明顯以及其適生植物具有嗜鈣性、耐旱性和石生性等限制特點(diǎn)的綜合影響,其生境和植被具有高度異質(zhì)性。同時(shí),該保護(hù)區(qū)位于我國(guó)植物區(qū)系的相匯交錯(cuò)區(qū)和交接過(guò)渡的中心地帶[10]。從地貌類(lèi)型來(lái)看,該區(qū)域?qū)儆诳λ固胤鍏餐莸氐孛病G叭搜芯勘砻?峰叢洼地集中分布在西南喀斯特地區(qū)南部斜坡地帶,面積約12.5萬(wàn)km2[11],其特點(diǎn)為絕對(duì)高差小、水熱條件差異大、植被垂直變化明顯、景觀異質(zhì)性強(qiáng)、巖石裸露率高、土層淺薄、土壤稀少而干燥且以石生植物為主,頂極群落具備生境復(fù)雜、植被生產(chǎn)力低、生物種類(lèi)多樣等總體特征[12]。喀斯特原生性森林是喀斯特植被演替的頂級(jí)群落,對(duì)深入研究喀斯特的植被演替或恢復(fù)具有重要參考價(jià)值。

    盡管諸多學(xué)者已對(duì)云南哈巴雪山、秦嶺太白山、廣東鼎湖山、海南尖峰嶺以及神龍架地區(qū)的植物多樣性海拔梯度分布格局進(jìn)行了研究[13- 18],但其研究結(jié)果不一定適用于亞熱帶喀斯特森林。雖然已有一些研究報(bào)導(dǎo)了喀斯特森林植被的空間格局[12,19- 22],但由于受到樣地大小或樣地形狀的限制,尚不能全面反映其多樣性與海拔梯度的聯(lián)系。為此,本研究通過(guò)建立一條從洼地到山頂?shù)莫M長(zhǎng)樣帶來(lái)分析喀斯特森林的群落類(lèi)型、結(jié)構(gòu)及其物種多樣性沿海拔梯度的變化趨勢(shì),以期為深入了解喀斯特原生性森林植物群落分布格局、闡明喀斯特退化植物群落的穩(wěn)定性機(jī)制提供科學(xué)依據(jù)。

    1 研究區(qū)概況與研究方法

    1.1 研究區(qū)概況

    木論國(guó)家自然保護(hù)區(qū)(107°54′01″—108°05′51″E,25°07′01″—25°12′22″N)位于廣西環(huán)江毛南族自治縣的西北部(圖1),東西寬19.80 km,南北長(zhǎng)10.75 km,林區(qū)總面積89.69 km2,屬于中亞熱帶石灰?guī)r區(qū)常綠落葉闊葉混交林生態(tài)系統(tǒng),是我國(guó)生物區(qū)系相匯交錯(cuò)區(qū)和交接過(guò)渡的中心,生境異質(zhì)性極高,區(qū)系成分復(fù)雜,生物種類(lèi)豐富,生態(tài)環(huán)境脆弱,是目前世界上喀斯特地區(qū)保存最完好、面積最大的原生林,森林覆蓋率達(dá)94.8%,其典型景觀單元為峰叢洼地。研究區(qū)屬中亞熱帶季風(fēng)氣候,年均氣溫19.3°C,極端最高溫36°C,極端最低溫- 5°C,1月平均3.4—8.7℃,7月平均23.0—26.7°C,≥10°C年活動(dòng)積溫6260°C,無(wú)霜期為310 d。年均日照1451.1 h,占1年可照時(shí)數(shù)(4422 h)的32.8%,月平均日照以7月份最多(190.7 h),2月份最少(61.5 h),年均降水量1529.2 mm (海拔420 m),降水集中在4—8月,占全年降水量的73.7%,林內(nèi)陰濕,年均相對(duì)濕度79%,干燥度為0.69,濕潤(rùn)度為1.46(K>1.0)?;鶐r裸露達(dá)80%—90%,土被面積不足20%,土壤發(fā)育于白云巖、石灰?guī)r的石灰土[23]。

    圖1 研究區(qū)及樣地的地理位置Fig.1 The location of study area and the investigated plots

    1.2 研究方法

    1.2.1 樣地設(shè)置

    本研究旨在反映喀斯特森林植物多樣性的海拔梯度格局,因此,在保護(hù)區(qū)內(nèi),從洼地底部至山頂建立一條較為整齊的森林群落樣帶,樣帶寬20 m,長(zhǎng)320 m,樣帶最低點(diǎn)海拔448.6 m,最高點(diǎn)海拔633.8 m。坡度范圍11.8—57.5 °,主坡向?yàn)槲髂?。整個(gè)樣帶包括16個(gè)投影面積20 m×20 m的連續(xù)樣方,樣方的4個(gè)角用水泥樁做永久標(biāo)記,在每個(gè)20 m×20 m的樣方用插值法細(xì)分為4個(gè)10 m×10 m和16個(gè)5 m×5 m的小樣方。以5 m ×5 m 小樣方為基本測(cè)量單元,按順時(shí)針掛牌標(biāo)記樣方內(nèi)每個(gè)胸徑(DBH)≥1 cm的木本植物個(gè)體,記錄物種名稱(chēng)、測(cè)定樹(shù)木的胸徑、樹(shù)高和冠幅、坐標(biāo)和植物生長(zhǎng)狀況(斜生、斷頂、倒伏等);利用GPS(Trimble XT6000)、手持羅盤(pán)測(cè)定每一調(diào)查樣地的經(jīng)緯度、海拔高度、坡度和坡向。調(diào)查樣地基本情況見(jiàn)表1。

    表1 調(diào)查樣地基本情況

    1.2.2 物種豐富度及其結(jié)構(gòu)成分

    (1) 物種豐富度 統(tǒng)計(jì)樣帶中維管束植物的科、屬、種數(shù),計(jì)算各樣方內(nèi)屬數(shù)/科數(shù)、種數(shù)/屬數(shù)、種數(shù)/科數(shù)的比值。

    (2) 生態(tài)成分 根據(jù)物種的生活型、落葉性劃分不同的生態(tài)成分,并計(jì)算其比值。

    1.2.3 群落數(shù)量分類(lèi)和排序

    采用等級(jí)聚類(lèi)分析(Hierarchical Cluster Analysis)方法,依據(jù)喬木層優(yōu)勢(shì)種重要值對(duì)樣方進(jìn)行聚類(lèi)分析;同時(shí),采用非度量多維標(biāo)度排序(Non-metric Multidimensional Scaling,NMS)法進(jìn)行矩陣排序。其中重要值計(jì)算方法為:

    IV=(A+D+F)/3

    式中,A為相對(duì)多度,指每個(gè)種的多度占所有種多度之和的百分比;D為相對(duì)優(yōu)勢(shì)度,指每個(gè)種的胸高斷面積之和占所有種個(gè)體胸高斷面積之和的百分比;F為相對(duì)頻度,指每個(gè)種的頻度占所有種頻度之和的百分比。

    1.2.4 多樣性測(cè)度

    豐富度指數(shù)S

    S=樣方內(nèi)物種數(shù)

    多樣性指數(shù)Pielou均勻度指數(shù)

    Jsw=-(∑PilnPi)/lnS

    Shannon-Wiener指數(shù)

    Simpson指數(shù)

    D=1-∑(Ni/N)2

    式中,Pi=Ni/N為第i個(gè)物種的個(gè)體在樣地中所占的比例;N為所有種的個(gè)體總數(shù);s為種數(shù);Ni為第i種的個(gè)體數(shù)。

    Bray-Curtis 指數(shù)

    CN= 2jN/(Na+Nb)

    式中,Na為樣地a的所有物種個(gè)體數(shù)目和,Nb為樣地b的所有物種個(gè)體數(shù)目和,jN為樣地a和b共有種中個(gè)體數(shù)目較小者之和,即jN=∑min(jNa,jNb)。

    所有統(tǒng)計(jì)分析均在R軟件(R Development CoreTeam,2009)中進(jìn)行。

    2 結(jié)果與分析

    2.1 物種組成與分類(lèi)

    2.1.1 群落種類(lèi)組成

    本次調(diào)查的16個(gè)森林群落樣方中喬木層共記錄到5089個(gè)胸徑≥1.0 cm的植株,喬木樹(shù)種120 種,隸屬于48科90屬,其中有落葉樹(shù)種21科37屬49種,常綠樹(shù)種34科58屬70種,針葉樹(shù)種僅有1種,落葉與常綠樹(shù)種分別占總數(shù)的40.83%和58.33%。樣地中重要值>10.00的物種共有25種,占總種數(shù)的20.8%,占總重要值的28.4%;樣帶中屬于單屬、單種的科有17種,占總科數(shù)的35.4%,而含有2個(gè)屬以下的科共有34個(gè),占總科數(shù)的70.8%;種數(shù)大于5的科有樟科(Lauraceae)、蕓香科(Rutaceae)、薔薇科(Rosaceae)、茜草科(Rubiaceae)和大戟科(Euphorbiaceae)5個(gè)科。種數(shù)最多的科為大戟科(14種)。

    由表2可知,洼地中多以小果厚殼桂(Cryptocaryaaustrokweichouensis)、梔子皮 (Itoaorientalis)為主,巖石裸露率高的樣方以灰?guī)r棒柄花(Cleidionbracteosum)、長(zhǎng)管越南茜(Rubovietnamiaaristata)、羅傘(Brassaiopsisglomerulata)為主,其次還有密榴木 (Miliusachunii);隨海拔上升,陸續(xù)出現(xiàn)廣西密花樹(shù)(Rapaneakwangsiensis)、化香樹(shù)(Platycaryastrobilacea)、密花樹(shù) (Rapaneaneriifolia)、羅城鵝耳櫪(Carpinusluochengensis),山頂則以清香木(Pistaciaweinmannifolia)、化香樹(shù)、密花樹(shù)和米念芭(Tirpitziaovoidea)為主。

    表2 各樣地重要值大于5的物種

    不同分類(lèi)等級(jí)之間的數(shù)量比(圖2) 基本上反映了該樣帶不同植物群落類(lèi)型在物種組成方面的差異。種數(shù)/科數(shù)、種數(shù)/屬數(shù)、屬數(shù)/科數(shù)的比值是對(duì)植物區(qū)系分化強(qiáng)度的反映,其值低則表示寡種屬多,區(qū)系的古老和殘遺性強(qiáng);其值高,則表明多種屬比例大,反映區(qū)系的物種分化強(qiáng)烈。由圖2可知各比例的最大值出現(xiàn)在y11樣方(海拔502.1 m)。

    2.1.2 群落生活型

    從圖 3可以看出,針葉樹(shù)種僅出現(xiàn)在13和14樣方,這可能與較高海拔的低溫有關(guān)??傮w來(lái)看,在整個(gè)海拔梯度上,常綠樹(shù)種的重要值比重始終遠(yuǎn)高于落葉樹(shù)種,表明該區(qū)域植物群落構(gòu)成以常綠樹(shù)種為主。

    圖2 不同系統(tǒng)等級(jí)分類(lèi)單元的多樣性梯度Fig.2 Altitudinal gradient of diversity of units at different taxonomic levels

    圖3 不同生活型植物重要值百分比沿海拔梯度的變化Fig.3 The altitudinal patterns of IV percentage in different life forms

    2.1.3 群落分類(lèi)與排序

    綜合等級(jí)聚類(lèi)分析和NMS排序結(jié)果,16個(gè)樣地被明顯地分為了3個(gè)類(lèi)型:第Ⅰ類(lèi),是以小果厚殼桂和灰?guī)r棒柄花為主要優(yōu)勢(shì)種的群落,包括1、2、3、4、5、6號(hào)樣地。第II類(lèi),是以長(zhǎng)管越南茜和羅傘為主要優(yōu)勢(shì)種的群落、包括7、8、9、10、11號(hào)樣地;第III類(lèi),是以化香樹(shù)和密花樹(shù)為主要優(yōu)勢(shì)種的群落、包括12、13、14、15、16號(hào)樣地。等級(jí)聚類(lèi)分析和NMS排序方法得出結(jié)果的一致,說(shuō)明分類(lèi)結(jié)果的合理性。另外,數(shù)量分類(lèi)的結(jié)果與現(xiàn)地植被類(lèi)型的分布格局相符,同樣印證了分類(lèi)結(jié)果的合理性。

    圖4 群落樣方的聚類(lèi)分析Fig.4 Cluster analysis classification of the community

    圖5 非度量多維標(biāo)度(NMS)排序Fig.5 Non-metric multidimensional scaling (NMS) of the species × attributes matrix

    2.2 物種豐富度

    圖6表示了物種豐富度沿海拔梯度的變化趨勢(shì)??梢钥闯?在低海拔地區(qū),隨海拔升高,木本植物的豐富度先緩慢增加,從7號(hào)樣方開(kāi)始迅速增大,在11號(hào)樣方達(dá)到峰值后開(kāi)始下降。也就是說(shuō),11號(hào)樣方的分類(lèi)群豐富度最高。表明該區(qū)域物種總數(shù)在中等海拔最大,符合“中間高度膨脹(mid-altitude bulge)”理論。

    2.3 群落多樣性

    2.3.1 α多樣性

    Shannon-Wiener多樣性指數(shù)與Simpson優(yōu)勢(shì)度指數(shù)沿海拔梯度表現(xiàn)出一致的趨勢(shì)(圖 7),隨著海拔升高,基本呈現(xiàn)出“低-高-低”的單峰分布格局,并均在12號(hào)樣方達(dá)到峰值。

    由圖8可知,Pielou均勻度指數(shù)在研究區(qū)域內(nèi)總體趨勢(shì)也是隨著海拔梯度先升高后逐漸降低,峰值同樣出現(xiàn)在12樣方。低海拔區(qū)域內(nèi),隨著海拔梯度的升高,外界干擾活動(dòng)逐漸減少,均勻度開(kāi)始逐漸升高,在中高海拔區(qū)域內(nèi),由于生境的異質(zhì)性和特殊性,只有特定的物種才能生存,物種類(lèi)型較為單一,均勻度也隨之下降。

    3個(gè)指數(shù)都反映出基本一致的變化趨勢(shì):即隨著海拔的升高先增加后減少,而3個(gè)指數(shù)均在 6號(hào)樣方突然減少,可能是由于該樣方巖石裸露率高,藤本覆蓋面積較大導(dǎo)致木本植物數(shù)量較少,且分支較多所致。

    圖6 豐富度隨海拔梯度變化趨勢(shì)圖Fig.6 Changes in species richness along an altitudinal gradient

    圖7 多樣性指數(shù)隨海拔梯度變化趨勢(shì)圖Fig.7 Species diversity along an altitudinal gradient

    2.3.2 β多樣性

    Bray-Curtis指數(shù)屬于相似性測(cè)度,不同地段共有種越多,其β多樣性越低,從圖9可以看出,1號(hào)樣地與其它不同海拔樣地喬木層的Bray-Curtis指數(shù)在海拔梯度上大致呈現(xiàn)先升高后下降再升高的趨勢(shì),表明隨著海拔落差的增加,其他樣方與1號(hào)樣方的共有物種先增加后減少再增加,β多樣性隨著海拔的升高減小。出現(xiàn)這種趨勢(shì)的原因可能是隨著海拔升高,植物群落進(jìn)一步發(fā)育,群落物種組成差異性不斷增加,群落之間的生態(tài)距離逐漸變遠(yuǎn),到達(dá)一定高度后干擾減少,群落所處的局部小環(huán)境更加穩(wěn)定。

    圖8 均勻度指數(shù)隨海拔梯度變化趨勢(shì)圖Fig.8 Pielou′s index along an altitudinal gradient

    圖9 1 號(hào)樣地與其它樣地不同尺度喬木層β樣多樣性Fig.9 The first quadrat′s beta diversity of tree layer vs that of each of the other quadrats′

    圖10 相鄰樣地喬木層β多樣性 Fig.10 The beta diversity in tree layer′ of each two adjacent quadrats

    圖10顯示了沿海拔梯度相鄰樣地間的β多樣性變化趨勢(shì),結(jié)果表明,相鄰樣地間種類(lèi)組成隨海拔梯度變化差異比較明顯,1—2、4—5及6—7樣地間種類(lèi)組成差異較小,相似性較高,3—4和8—9樣地間差異較大,相似性較低。

    3 討論

    喀斯特峰叢洼地景觀異質(zhì)性強(qiáng),巖石裸露率達(dá)80%以上,土層淺薄,土壤稀少而干燥,植物多生長(zhǎng)在巖石上或石縫中,既要有石生性、喜鈣性和耐旱性的特點(diǎn),又必須具有發(fā)達(dá)而強(qiáng)壯的根系才能扎根和生長(zhǎng)[12]。木論喀斯特森林群落物種組成復(fù)雜多樣,根據(jù)本文研究結(jié)果,在16個(gè)樣地0.64 hm2內(nèi)統(tǒng)計(jì)到48科90屬120種植物,與宋同清等[12]在該保護(hù)區(qū)0.8 hm2內(nèi)調(diào)查到的43科91屬123種木本植物的結(jié)果基本一致。相比而言:遠(yuǎn)高于長(zhǎng)白山25 hm2闊葉紅松林樣地的52種[24]和東靈山20 hm2樣地的58種[25],低于古田山24 hm2樣地[26]、鼎湖山[27]和西雙版納[28]20 hm2樣地分別有的159、210和468種。

    植物多樣性與海拔梯度的關(guān)系錯(cuò)綜復(fù)雜,與群落類(lèi)型、水熱條件、所處的土壤基質(zhì)和取樣尺度等不同有關(guān)[1]。先前的研究一般認(rèn)為:物種多樣性具有多種海拔梯度分布格局,其中最普遍的觀點(diǎn)是物種多樣性的海拔梯度格局與緯度格局相似,即隨著海拔上升物種多樣性降低[15,29- 30],但也有研究表明隨著海拔升高,物種多樣性呈現(xiàn)單峰分布格局,即物種多樣性隨著海拔上升先增加、后減少[31- 33]或呈現(xiàn)出更為復(fù)雜的變化規(guī)律[34- 35]。

    相對(duì)海拔是影響喀斯特山地水熱條件分布的主要因素,它通過(guò)影響光照、水分及土壤資源的再分配進(jìn)而影響植物群叢的分布及結(jié)構(gòu)[36]。本研究所選樣帶的α多樣性與常見(jiàn)的隨海拔升高而下降的模式[30]不同,表現(xiàn)出與佛坪自然保護(hù)區(qū)[37]和五鹿山國(guó)家級(jí)自然保護(hù)區(qū)物種多樣性海拔格局[38]一致的“單峰格局”。主要原因可能是:由于樣地選址位于喀斯特峰叢洼地,處于亞熱帶區(qū)域,海拔最高不超過(guò)700 m,隨著海拔升高,溫度下降不明顯,仍然屬于適合植物生長(zhǎng)的范圍內(nèi),喀斯特山地地形復(fù)雜、坡位變化及其它環(huán)境資源因子分配差異大,低海拔洼地,土層較深厚、濕度條件較好,但太陽(yáng)直射光照時(shí)間不長(zhǎng),一般僅限于某些喜蔭和耐蔭性強(qiáng)的種類(lèi),因而較相對(duì)于光照充足的相鄰的中等海拔山坡多樣性和物種種類(lèi)組成都低。而在高海拔地段山體陡峭、少土、日照時(shí)間長(zhǎng)導(dǎo)致土壤蒸發(fā)強(qiáng)烈,造成嚴(yán)重干旱缺水等不利植物定居的惡劣環(huán)境,因而多樣性隨海拔變化而呈現(xiàn)下降趨勢(shì)。另外海拔梯度實(shí)際上在一定程度反映了水熱狀況的梯度變化,因此海拔的影響實(shí)際反映了氣候因子的綜合影響[39]。很多情況下,在低海拔地段,植物生長(zhǎng)季節(jié)時(shí)所需熱量充足而水分常常不足;在高海拔地段,水分充足而熱量可能成為植物分布的限制因子。相對(duì)于低海拔和高海拔,在中海拔地段可能是熱量和水分組合最好的生境,從而使得資源的可利用率較高[32]。6號(hào)樣方出現(xiàn)的多樣性突然降低是由于該樣地的特殊地形,因此在以后的研究中應(yīng)該更多的考慮局部小環(huán)境的影響。

    β多樣性表示沿某一環(huán)境梯度物種替代速率,能夠較直觀地反映不同群落間物種的組成差異和不同生境間的梯度變化[40]。在對(duì)β多樣性隨海拔梯度變化格局的研究中,數(shù)量數(shù)據(jù)測(cè)度和二元屬性數(shù)據(jù)法各有利弊,不同學(xué)者也得出了不同結(jié)論:鼎湖山群落β多樣性隨海拔呈無(wú)規(guī)律變化[16];南嶺國(guó)家級(jí)自然保護(hù)區(qū)森林群落β多樣性隨海拔升高而下降[41];長(zhǎng)白山北坡植物群落β多樣性隨群落間海拔差的增加而增大[42]。但生態(tài)位理論認(rèn)為群落相似性的降低是由于環(huán)境異質(zhì)性的增加[33],中性理論則認(rèn)為擴(kuò)散限制產(chǎn)生了β多樣性[30]。所以本研究區(qū)域的森林群落β多樣性符合隨海拔升高而減小的常見(jiàn)格局[16]。

    本研究選取的木論保護(hù)區(qū)屬于典型的喀斯特峰叢洼地地貌,從洼地到山頂?shù)慕^對(duì)高度不大,水熱條件差異大,植被垂直變化明顯。從低海拔往高海拔16個(gè)樣地被明顯地分為了3個(gè)類(lèi)型:第一類(lèi)是洼地及其邊緣,由于光照時(shí)間不長(zhǎng),比較蔭蔽,濕度條件較好,故分布著以小果厚殼桂和灰?guī)r棒柄花為主要優(yōu)勢(shì)種的喜濕和耐蔭性強(qiáng)的植物群落。第二類(lèi),是山坡山腰地段,這里因太陽(yáng)直射時(shí)間比洼地長(zhǎng),土壤較少,土層較淺薄,保水性能較差,濕度條件不如洼地,群落主要優(yōu)勢(shì)種為長(zhǎng)管越南茜和羅傘;第三類(lèi),是山頂及其周?chē)?這里環(huán)境條件惡劣,全天受到太陽(yáng)直射,溫度高,土壤極少,巖石裸露,加之山頂上風(fēng)大、保水性能和濕度條件最差,群落主要優(yōu)勢(shì)種為化香樹(shù)和密花樹(shù)等耐旱植物。因此,在喀斯特山地植被恢復(fù)時(shí),要根據(jù)植物群落物種豐富度的海拔梯度規(guī)律,科學(xué)、合理地進(jìn)行立地類(lèi)型劃分,選擇適宜栽種物種,制訂相應(yīng)技術(shù)和管理措施,開(kāi)展高效的植物群落營(yíng)建和封育管護(hù),才能達(dá)到費(fèi)省效宏的效果。

    致謝:樣地調(diào)查得到木論國(guó)家級(jí)自然保護(hù)區(qū)管理局工作人員和中南林業(yè)科技大學(xué)2012級(jí)生態(tài)學(xué)專(zhuān)業(yè)實(shí)習(xí)同學(xué)的大力支持,特此致謝。

    [1] 方精云. 探索中國(guó)山地植物多樣性的分布規(guī)律. 生物多樣性, 2004, 12(1): 1- 4.

    [2] Brown J H. Mammals on mountainsides: elevational patterns of diversity. Global Ecology and Biogeography, 2001, 10(1): 101- 109.

    [3] Hegazy A K, El-Demerdash M A, Hosni H A. Vegetation, species diversity and floristic relations along an altitudinal gradient in south-west Saudi Arabia. Journal of Arid Environments, 1998, 38(1): 3- 13.

    [4] Küchler A W. Review: vegetation of the Earth. Ecology, 1980, 61(1): 206- 206.

    [5] O′Brien E M. Climatic gradients in woody plant species richness: towards an explanation based on an analysis of southern Africa′s woody flora. Journal of Biogeography, 1993, 20(2): 181- 198.

    [6] 唐志堯, 方精云. 植物物種多樣性的垂直分布格局. 生物多樣性, 2004, 12(1): 20- 28.

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

    [8] 蘇宗明, 李先琨. 廣西巖溶植被類(lèi)型及其分類(lèi)系統(tǒng). 廣西植物, 2003, 23(4): 289- 293.

    [9] Hou M F, López-Pujol J, Qin H N, Wang L S, Liu Y. Distribution pattern and conservation priorities for vascular plants in Southern China: Guangxi Province as a case study. Botanical Studies, 2010, 51(3): 377- 386.

    [10] 鄭穎吾. 木論喀斯特林區(qū)概論. 北京: 科學(xué)出版社, 1999: 36-38.

    [11] 蔣忠誠(chéng), 李先琨, 曾馥平. 巖溶峰叢洼地生態(tài)重建. 北京: 地質(zhì)出版社, 2007: 16-17.

    [12] 宋同清, 彭晚霞, 曾馥平, 王克林, 曹洪麟, 李先琨, 覃文更, 譚衛(wèi)寧, 劉璐. 喀斯特峰叢洼地不同類(lèi)型森林群落的組成與生物多樣性特征. 生物多樣性, 2010, 18(4): 355- 364.

    [13] 沈澤昊, 胡會(huì)峰, 周宇, 方精云. 神農(nóng)架南坡植物群落多樣性的海拔梯度格局. 生物多樣性, 2004, 12(1): 99- 107.

    [14] 陶晶, 臧潤(rùn)國(guó), 余昌元. 云南哈巴雪山植物群落和植物多樣性海拔梯度分布格局. 林業(yè)科學(xué), 2011, 47(7): 1- 6.

    [15] 唐志堯, 方精云, 張玲. 秦嶺太白山木本植物物種多樣性的梯度格局及環(huán)境解釋. 生物多樣性, 2004, 12(1): 115- 122.

    [16] 林國(guó)俊, 黃忠良, 竺琳, 歐陽(yáng)學(xué)軍. 鼎湖山森林群落β多樣性. 生態(tài)學(xué)報(bào), 2010, 30(18): 4875- 4880.

    [17] 李林, 周小勇, 黃忠良, 魏識(shí)廣, 史軍輝. 鼎湖山植物群落α多樣性與環(huán)境的關(guān)系. 生態(tài)學(xué)報(bào), 2006, 26(7): 2301- 2307.

    [18] 吳裕鵬, 許涵, 李意德, 駱土壽, 陳德祥, 林明獻(xiàn). 海南尖峰嶺熱帶林喬灌木層物種多樣性沿海拔梯度分布格局. 林業(yè)科學(xué), 2013, 49(4): 16- 23.

    [19] 蔣國(guó)芳. 木論喀斯特森林區(qū)生物多樣性初報(bào). 生物多樣性, 1995, 3(2): 91- 94.

    [20] 宋同清, 彭晚霞, 曾馥平, 王克林, 覃文更, 譚衛(wèi)寧, 劉璐, 杜虎, 鹿士楊. 木論喀斯特峰叢洼地森林群落空間格局及環(huán)境解釋. 植物生態(tài)學(xué)報(bào), 2010, 34(3): 298- 308.

    [21] 韓文衡, 向悟生, 葉鐸, 呂仕洪, 丁濤, 李先琨. 廣西木論保護(hù)區(qū)喀斯特常綠落葉闊葉混交林優(yōu)勢(shì)種空間格局及其相關(guān)性. 應(yīng)用生態(tài)學(xué)報(bào), 2010, 21(11): 2769- 2776.

    [22] 彭晚霞, 宋同清, 曾馥平, 王克林, 劉璐. 喀斯特常綠落葉闊葉混交林植被的空間異質(zhì)性. 西北植物學(xué)報(bào), 2011, 31(4): 815- 822.

    [23] 宋同清. 西南喀斯特植物與環(huán)境. 北京: 科學(xué)出版社, 2015:115-123.

    [24] 郝占慶, 李步杭, 張健, 王緒高, 葉吉, 姚曉琳. 長(zhǎng)白山闊葉紅松林樣地(CBS): 群落組成與結(jié)構(gòu). 植物生態(tài)學(xué)報(bào), 2008, 32(2): 238- 250.

    [25] 劉海豐, 李亮, 桑衛(wèi)國(guó). 東靈山暖溫帶落葉闊葉次生林動(dòng)態(tài)監(jiān)測(cè)樣地: 物種組成與群落結(jié)構(gòu). 生物多樣性, 2011, 19(2): 232- 242.

    [26] 祝燕, 趙谷風(fēng), 張儷文, 沈國(guó)春, 米湘成, 任海保, 于明堅(jiān), 陳建華, 陳聲文, 方騰, 馬克平. 古田山中亞熱帶常綠闊葉林動(dòng)態(tài)監(jiān)測(cè)樣地——群落組成與結(jié)構(gòu). 植物生態(tài)學(xué)報(bào), 2008, 32(2): 262- 273.

    [27] 葉萬(wàn)輝, 曹洪麟, 黃忠良, 練琚愉, 王志高, 李林, 魏識(shí)廣, 王章明. 鼎湖山南亞熱帶常綠闊葉林20公頃樣地群落特征研究. 植物生態(tài)學(xué)報(bào), 2008, 32(2): 274- 286.

    [28] 蘭國(guó)玉, 胡躍華, 曹敏, 朱華, 王洪, 周仕順, 鄧曉保, 崔景云, 黃建國(guó), 劉林云, 許海龍, 宋軍平, 何有才. 西雙版納熱帶森林動(dòng)態(tài)監(jiān)測(cè)樣地—樹(shù)種組成與空間分布格局. 植物生態(tài)學(xué)報(bào), 2008, 32(2): 287- 298.

    [29] 郝占慶, 于德永, 楊曉明, 丁之慧. 長(zhǎng)白山北坡植物群落α多樣性及其隨海拔梯度的變化. 應(yīng)用生態(tài)學(xué)報(bào), 2002, 13(7): 786- 789.

    [30] Vázquez G J A, Givnish T J. Altitudinal gradients in tropical forest composition, structure, and diversity in the Sierra de Manantlán. Journal of Ecology, 1998, 86(6): 999- 1020.

    [31] Rahbek C. The role of spatial scale and the perception of large-scale species-richness patterns. Ecology Letters, 2005, 8(2): 224- 239.

    [32] 王國(guó)宏. 祁連山北坡中段植物群落多樣性的垂直分布格局. 生物多樣性, 2002, 10(1): 7- 14.

    [33] Wang G H, Zhou G S, Yang L M, Li Z Q. Distribution, species diversity and life-form spectra of plant communities along an altitudinal gradient in the northern slopes of Qilianshan Mountains, Gansu, China. Plant Ecology, 2003, 165(2): 169- 181.

    [34] 張華, 鄭楠, 武晶, 祝業(yè)平, 李劍源, 張寶財(cái), 孫衛(wèi)東. 遼東山地老禿頂子北坡植物物種多樣性的垂直分布格局. 遼寧師范大學(xué)學(xué)報(bào): 自然科學(xué)版, 2009, 32(4): 484- 489.

    [35] Dole?al J,tek M. Altitudinal changes in composition and structure of mountain-temperate vegetation: a case study from the Western Carpathians. Plant Ecology, 2002, 158(2): 201- 221.

    [36] 袁道先, 章程. 巖溶動(dòng)力學(xué)的理論探索與實(shí)踐. 地球?qū)W報(bào), 2008, 29(3): 355- 365.

    [37] 王宇超, 王得祥. 佛坪自然保護(hù)區(qū)植物群落物種多樣性和復(fù)雜性的海拔格局研究. 西北植物學(xué)報(bào), 2013, 33(1): 169- 176.

    [38] 何艷華, 閆明, 張欽弟, 張直峰, 苗艷明, 畢潤(rùn)成. 五鹿山國(guó)家級(jí)自然保護(hù)區(qū)物種多樣性海拔格局. 生態(tài)學(xué)報(bào), 2013, 33(8): 2452- 2462.

    [39] Whittaker R J, Willis K J, Field R. Scale and species richness: towards a general, hierarchical theory of species diversity. Journal of Biogeography, 2001, 28(4): 453- 470.

    [40] 馬克平, 劉燦然, 劉玉明. 生物群落多樣性的測(cè)度方法——Ⅱβ多樣性的測(cè)度方法. 生物多樣性, 1995, 3(1): 38- 43.

    [41] 張璐, 蘇志堯, 李鎮(zhèn)魁. 南嶺國(guó)家級(jí)自然保護(hù)區(qū)森林群落β多樣性隨海拔梯度的變化. 熱帶亞熱帶植物學(xué)報(bào), 2007, 15(6): 506- 512.

    [42] 郝占慶, 于德永, 吳鋼, 鄧紅兵, 姜萍, 王慶禮. 長(zhǎng)白山北坡植物群落β多樣性分析. 生態(tài)學(xué)報(bào), 2001, 21(12): 2018- 20

    Altitudinal pattern of woody plant species diversity in the karst forest in Mulun, China

    LAN Sian1,2,SONG Min1,2,3,ZENG Fuping1,2,SONG Tongqing1,2,*,PENG Wanxia1,2,DU Hu1,2,HAN Chang1,2,4, CHEN Li1,2,4,QIN Wengeng5

    1KeyLaboratoryofAgro-ecologicalProcessesinSubtropicalRegion,InstituteofSubtropicalAgriculture,ChineseAcademyofSciences,Changsha410125,China2HuanjiangObservationandResearchStationofKarstEcosystem,ChineseAcademyofSciences,Huanjiang547100,China3CollegeofBioscienceandBiotechnology,HunanAgriculturalUniversity,Changsha410128,China4CollegeofAgriculture,HunanAgriculturalUniversity,Changsha410128,China5AdministrativeBureauofMulunNationalNatureReserve,GuangxiZhuangAutonomousRegion,Huanjiang547200,China

    Variation in species diversity along ecological gradients is the main focus in recent biodiversity studies. Among multiple interacting resource gradients, altitude gradient may be a decisive factor determining spatial patterns of species diversity. Forest community classification, ordination, and species diversity patterns along altitudinal gradients can reveal relationships between species distribution and environmental factors, and provide an important basis for understanding the function of mountain ecosystems. For studying the species diversity along ecological gradients in karst ecosystems, we conducted a study in Mulun National Nature Reserve (MNNR), northwest of Huanjiang Maonan Autonomous County in Guangxi, China. This reserve and Maolan National Nature Reserve together possesses the largest and most typical karst forests in the world, and play an important role in plant biodiversity research due to their unique geography, climate, and vegetation conditions. The climax community in this region is an evergreen and deciduous broadleaf mixed forest on limestone in the central subtropical zone in MNNR. The area of this forest is 89.69 km2with 94.8% of vegetation coverage. The soil type is calcareous lithosols (limestone soil). Sixteen sample plots (20 m × 20 m) were established in a typical depression landscape between karst hills on a southwestern slope in MNNR. All woody plants [DBH (diameter at breast height) ≥ 1 cm] from the 16 plots were investigated. We analyzed the pattern of woody plant species diversity along an altitudinal gradient in MNNR by measuring the indices of species richness, diversity, and evenness. The results showed that: (1) there were 5089 individuals in the studied plots, which belonged to 120 species, 48 families and 90 genera; (2) there was a distinct altitudinal pattern of woody plant species in the karst forest communities. Hierarchical cluster analysis and Non-metric Multidimensional Scaling analysis showed that the forest communities could be classified into three types, including those in lower slope, middle slope and upper slope respectively. In the lower slope, the dominant species wereCryptocaryaaustrokweichouensisandCleidionbracteosumare; in the middle slope, the dominant species wereRubovietnamiaaristataandBrassaiopsisglomerulata;PlatycaryastrobilaceaandRapaneaneriifoliaarewere the dominant species in the upper slope. As the altitude increases, the temperature and moisture also changed. Aridity and low temperature respectively limited plant growth at low and high altitude; the precipitation and sunshine were plentiful in middle elevations, which supports a high diversity of plant species; (3) among the woody species, there was a greater proportion of evergreen species than deciduous broadleaf species. In addition, coniferous species only grew on the top of the mountain; (4) Species richness and the Shannon-Wiener index, along the altitudinal gradient, exhibited a pattern of unimodal distribution, with the peak appearing in the mid slope; (5) theβdiversity index fluctuated considerably, and decreased with rising elevation in the ecosystem. Our results will be useful for developing strategies of vegetation restoration in the karst region.

    altitudinal gradient;αdiversity;βdiversity; woody species; karst

    中國(guó)科學(xué)院科技服務(wù)網(wǎng)絡(luò)計(jì)劃STS資助項(xiàng)目(KFJ-EW-STS-092); 國(guó)家自然科學(xué)基金資助項(xiàng)目(31370485, 31370623, 31400412, 31460135); 廣西科技惠民資助項(xiàng)目(桂科轉(zhuǎn)1599001- 6);廣西特聘專(zhuān)家資助項(xiàng)目; 湖南省研究生科研創(chuàng)新資助項(xiàng)目(CX2015B231)

    2016- 01- 08;

    2016- 05- 06

    10.5846/stxb201601080055

    *通訊作者Corresponding author.E-mail: songtongq@isa.ac.cn

    蘭斯安, 宋敏, 曾馥平, 宋同清, 彭晚霞, 杜虎, 韓暢, 陳莉, 覃文更.木論喀斯特森林木本植物多樣性垂直格局.生態(tài)學(xué)報(bào),2016,36(22):7374- 7383.

    Lan S A,Song M,Zeng F P,Song T Q,Peng W X,Du H,Han C, Chen L,Qin W G.Altitudinal pattern of woody plant species diversity in the karst forest in Mulun, China.Acta Ecologica Sinica,2016,36(22):7374- 7383.

    猜你喜歡
    喀斯特樣方樣地
    藥王谷多種因素對(duì)漏蘆、野韭生長(zhǎng)的影響
    昆明金殿國(guó)家森林公園林業(yè)有害生物調(diào)查研究
    額爾古納市興安落葉松中齡林植被碳儲(chǔ)量研究
    典型野生刺梨(Rosa roxburghii Tratt.)灌草叢植物多樣性研究
    昆明市主要綠化樹(shù)種閾值測(cè)定與分析
    基于角尺度模型的林業(yè)樣地空間結(jié)構(gòu)分析
    大自然的鬼斧神工:“中國(guó)南方喀斯特”
    云南九龍山自然保護(hù)區(qū)水青岡群落物種多樣性及群落特征
    “別有洞天”的喀斯特王國(guó)
    喀斯特
    ————水溶蝕巖石的奇觀
    家教世界(2019年4期)2019-02-26 13:44:20
    亚洲中文字幕一区二区三区有码在线看| 夜夜夜夜夜久久久久| 国产高清三级在线| www.av在线官网国产| 国产精品久久久久久亚洲av鲁大| 国产成年人精品一区二区| 一个人观看的视频www高清免费观看| 国产视频内射| 午夜视频国产福利| 日日啪夜夜撸| 中国国产av一级| 国产人妻一区二区三区在| 国产 一区 欧美 日韩| 一个人看的www免费观看视频| 99国产极品粉嫩在线观看| 国产av不卡久久| 午夜免费激情av| 久久精品国产亚洲av天美| 一个人看视频在线观看www免费| av天堂在线播放| av免费在线看不卡| 亚洲内射少妇av| 亚洲无线在线观看| 观看美女的网站| 99久国产av精品国产电影| 欧美丝袜亚洲另类| 女的被弄到高潮叫床怎么办| 午夜激情欧美在线| 麻豆av噜噜一区二区三区| 欧美色欧美亚洲另类二区| 可以在线观看毛片的网站| 日本黄色视频三级网站网址| 久久久久久久久大av| 卡戴珊不雅视频在线播放| 又爽又黄无遮挡网站| 日韩欧美在线乱码| 国产午夜精品久久久久久一区二区三区| 亚洲无线在线观看| 午夜激情福利司机影院| 亚洲人成网站在线观看播放| 亚洲国产精品sss在线观看| 观看美女的网站| 亚洲成a人片在线一区二区| 精品久久久久久久人妻蜜臀av| 国产免费一级a男人的天堂| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产亚洲av片在线观看秒播厂 | 99久久成人亚洲精品观看| 看免费成人av毛片| 成人三级黄色视频| 一区二区三区免费毛片| 99久久成人亚洲精品观看| 夫妻性生交免费视频一级片| 亚洲18禁久久av| 午夜免费激情av| 丝袜喷水一区| 不卡一级毛片| 舔av片在线| 99热6这里只有精品| 亚洲av中文字字幕乱码综合| av在线观看视频网站免费| 噜噜噜噜噜久久久久久91| 最近2019中文字幕mv第一页| 亚洲无线观看免费| 久久国产乱子免费精品| 毛片一级片免费看久久久久| 亚洲av成人精品一区久久| 亚洲人成网站在线播放欧美日韩| 又粗又硬又长又爽又黄的视频 | www.av在线官网国产| 国产麻豆成人av免费视频| 亚洲久久久久久中文字幕| 欧美不卡视频在线免费观看| 免费av毛片视频| 高清毛片免费观看视频网站| 精品久久久久久成人av| 不卡视频在线观看欧美| 美女内射精品一级片tv| 超碰av人人做人人爽久久| 麻豆av噜噜一区二区三区| 精品少妇黑人巨大在线播放 | 哪里可以看免费的av片| 好男人在线观看高清免费视频| 成人综合一区亚洲| 1000部很黄的大片| 最近视频中文字幕2019在线8| 秋霞在线观看毛片| 小蜜桃在线观看免费完整版高清| 美女xxoo啪啪120秒动态图| 亚洲四区av| 成人美女网站在线观看视频| 国产午夜福利久久久久久| 国产一区亚洲一区在线观看| 99热这里只有是精品在线观看| 久久精品国产亚洲av涩爱 | 久久人妻av系列| 精品熟女少妇av免费看| 国产精品一区二区三区四区久久| 99热全是精品| 国产精品.久久久| 久久精品夜夜夜夜夜久久蜜豆| 欧美在线一区亚洲| 精品久久久久久久久久免费视频| 日韩精品青青久久久久久| 欧美性猛交黑人性爽| 99热全是精品| 免费无遮挡裸体视频| 99久久精品热视频| 最近2019中文字幕mv第一页| 国产精品麻豆人妻色哟哟久久 | 观看美女的网站| 一本久久精品| 午夜爱爱视频在线播放| 欧美极品一区二区三区四区| 免费观看精品视频网站| 亚洲在线观看片| 可以在线观看毛片的网站| 久久鲁丝午夜福利片| 欧美三级亚洲精品| 啦啦啦啦在线视频资源| 日本爱情动作片www.在线观看| 亚洲婷婷狠狠爱综合网| 亚洲av免费高清在线观看| 国产69精品久久久久777片| 亚洲不卡免费看| 99久久九九国产精品国产免费| 一区福利在线观看| 性色avwww在线观看| 联通29元200g的流量卡| 亚洲精品自拍成人| 69人妻影院| 校园春色视频在线观看| 国产女主播在线喷水免费视频网站 | 免费av不卡在线播放| 国产精品av视频在线免费观看| 男人的好看免费观看在线视频| 国产日韩欧美在线精品| 亚洲激情五月婷婷啪啪| 午夜亚洲福利在线播放| 久久久a久久爽久久v久久| 悠悠久久av| 亚洲最大成人av| 国产一级毛片在线| 黄色日韩在线| 97在线视频观看| 亚洲精品456在线播放app| 国产精品.久久久| 国产精品人妻久久久影院| 国产精品久久久久久精品电影小说 | 国产毛片a区久久久久| 免费电影在线观看免费观看| 热99re8久久精品国产| 91狼人影院| 免费观看a级毛片全部| 亚洲五月天丁香| 亚洲天堂国产精品一区在线| 男的添女的下面高潮视频| 青春草亚洲视频在线观看| 特大巨黑吊av在线直播| 亚洲人成网站在线播放欧美日韩| 精品免费久久久久久久清纯| 亚洲av不卡在线观看| 亚洲成人av在线免费| 长腿黑丝高跟| 成人亚洲精品av一区二区| 国产亚洲精品久久久久久毛片| 麻豆国产97在线/欧美| 国产精品一二三区在线看| 老熟妇乱子伦视频在线观看| 成年女人永久免费观看视频| 又爽又黄a免费视频| 天天躁夜夜躁狠狠久久av| 我要搜黄色片| 老师上课跳d突然被开到最大视频| 12—13女人毛片做爰片一| 久久综合国产亚洲精品| 精品久久久久久久久久久久久| av在线蜜桃| 中国国产av一级| 国产av一区在线观看免费| 在线观看免费视频日本深夜| 国产爱豆传媒在线观看| 欧美日韩在线观看h| 国产私拍福利视频在线观看| 淫秽高清视频在线观看| 99久久精品国产国产毛片| 丝袜喷水一区| 色噜噜av男人的天堂激情| 久久久色成人| 看十八女毛片水多多多| 综合色av麻豆| 国产又黄又爽又无遮挡在线| 少妇人妻一区二区三区视频| 高清毛片免费观看视频网站| 国产伦一二天堂av在线观看| 美女国产视频在线观看| 日本黄色片子视频| 成人毛片a级毛片在线播放| 欧美日韩乱码在线| 亚洲av中文字字幕乱码综合| 日日撸夜夜添| 国产精品蜜桃在线观看 | 女人十人毛片免费观看3o分钟| 黄色欧美视频在线观看| 在线观看一区二区三区| 精品欧美国产一区二区三| 特级一级黄色大片| 国产精品免费一区二区三区在线| 久久人人爽人人片av| 欧美最新免费一区二区三区| 午夜亚洲福利在线播放| 精品国内亚洲2022精品成人| 国产精品av视频在线免费观看| 国产成人aa在线观看| 黄色配什么色好看| 综合色av麻豆| 青春草视频在线免费观看| 深夜a级毛片| 国内精品宾馆在线| 伦精品一区二区三区| av免费在线看不卡| 精品99又大又爽又粗少妇毛片| av在线蜜桃| 一区二区三区四区激情视频 | 我要搜黄色片| 国产成人精品一,二区 | 一级毛片aaaaaa免费看小| av在线观看视频网站免费| 女人十人毛片免费观看3o分钟| 国产伦精品一区二区三区四那| 国产不卡一卡二| 久久久国产成人免费| 亚洲最大成人中文| 一级毛片电影观看 | 亚洲欧美中文字幕日韩二区| 久久久久久久久久成人| av免费观看日本| 好男人在线观看高清免费视频| 亚洲国产精品合色在线| 成人特级av手机在线观看| 最近的中文字幕免费完整| 欧美另类亚洲清纯唯美| 日日撸夜夜添| 亚洲成人久久爱视频| 大香蕉久久网| 中文精品一卡2卡3卡4更新| 国产亚洲欧美98| 看免费成人av毛片| 亚洲国产欧美在线一区| 欧美xxxx黑人xx丫x性爽| 国产精品久久电影中文字幕| av.在线天堂| 听说在线观看完整版免费高清| 亚洲精品日韩av片在线观看| 色播亚洲综合网| 精品久久久噜噜| av在线观看视频网站免费| 国产精品久久视频播放| 天堂影院成人在线观看| 在线观看免费视频日本深夜| 日日撸夜夜添| 午夜精品在线福利| 国产大屁股一区二区在线视频| 99久久九九国产精品国产免费| 一级毛片电影观看 | 色尼玛亚洲综合影院| 亚洲精品日韩在线中文字幕 | 日韩精品青青久久久久久| 男插女下体视频免费在线播放| 熟妇人妻久久中文字幕3abv| 淫秽高清视频在线观看| a级一级毛片免费在线观看| 我的女老师完整版在线观看| 精品日产1卡2卡| 男女啪啪激烈高潮av片| 最近手机中文字幕大全| 亚洲成人av在线免费| 国产黄片视频在线免费观看| 日韩成人伦理影院| 中文精品一卡2卡3卡4更新| 亚洲图色成人| 国产不卡一卡二| 国产精品一二三区在线看| 边亲边吃奶的免费视频| 99久久精品国产国产毛片| 网址你懂的国产日韩在线| 日韩强制内射视频| 一级毛片aaaaaa免费看小| 高清日韩中文字幕在线| 久久精品国产亚洲av天美| 国产一区二区三区av在线 | eeuss影院久久| 热99re8久久精品国产| 亚洲一区二区三区色噜噜| av卡一久久| 精品少妇黑人巨大在线播放 | 成年免费大片在线观看| 男人舔奶头视频| 午夜福利视频1000在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产一级毛片在线| 亚洲国产精品sss在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲电影在线观看av| 亚洲高清免费不卡视频| 欧美不卡视频在线免费观看| 日韩av在线大香蕉| 看十八女毛片水多多多| 国产成人a∨麻豆精品| 天天一区二区日本电影三级| 日日干狠狠操夜夜爽| АⅤ资源中文在线天堂| 国国产精品蜜臀av免费| 国产成人a区在线观看| 男女做爰动态图高潮gif福利片| 啦啦啦啦在线视频资源| 亚洲自拍偷在线| 久久久久性生活片| 人人妻人人澡欧美一区二区| 成年免费大片在线观看| 亚洲欧美精品专区久久| 性欧美人与动物交配| 亚洲精品自拍成人| 麻豆乱淫一区二区| 欧美又色又爽又黄视频| 少妇人妻精品综合一区二区 | 夜夜爽天天搞| 波野结衣二区三区在线| 夜夜爽天天搞| 成人亚洲欧美一区二区av| 免费人成视频x8x8入口观看| 国产毛片a区久久久久| 日本av手机在线免费观看| 国产精品人妻久久久久久| 成年女人永久免费观看视频| 国产在线男女| 嫩草影院精品99| 国产在线男女| 中国国产av一级| 18+在线观看网站| 成人特级av手机在线观看| 国产欧美日韩精品一区二区| 成人特级av手机在线观看| 日韩一区二区视频免费看| 亚洲18禁久久av| 久久久久久久亚洲中文字幕| 亚洲欧美日韩无卡精品| 亚洲欧美清纯卡通| 久久综合国产亚洲精品| 国产视频内射| 男的添女的下面高潮视频| 欧美xxxx黑人xx丫x性爽| av专区在线播放| 亚洲精品乱码久久久久久按摩| 美女内射精品一级片tv| 亚洲国产日韩欧美精品在线观看| 亚洲精品久久久久久婷婷小说 | 亚洲成人久久性| 国产精华一区二区三区| 三级国产精品欧美在线观看| 一本一本综合久久| 久久热精品热| 国产精品一二三区在线看| 久久99精品国语久久久| 嫩草影院新地址| 日韩大尺度精品在线看网址| 人妻久久中文字幕网| 欧美变态另类bdsm刘玥| 日韩欧美精品v在线| 国产精华一区二区三区| 国产伦一二天堂av在线观看| 色哟哟·www| 国产精品永久免费网站| ponron亚洲| 亚洲精品国产成人久久av| 国产一区二区在线观看日韩| 中国美女看黄片| 欧美日本视频| 美女国产视频在线观看| 国产美女午夜福利| 久久久久国产网址| 麻豆国产97在线/欧美| 久久精品国产亚洲av天美| 亚洲激情五月婷婷啪啪| 国产黄片视频在线免费观看| 国产 一区 欧美 日韩| 亚洲乱码一区二区免费版| 观看美女的网站| 人人妻人人看人人澡| 国产一区二区三区在线臀色熟女| 久久精品国产亚洲av涩爱 | 亚洲,欧美,日韩| 深爱激情五月婷婷| 日韩,欧美,国产一区二区三区 | 精品久久久久久成人av| 免费av毛片视频| av在线老鸭窝| АⅤ资源中文在线天堂| 亚洲精华国产精华液的使用体验 | 亚洲经典国产精华液单| 最后的刺客免费高清国语| 国产老妇女一区| 18禁裸乳无遮挡免费网站照片| 99视频精品全部免费 在线| 高清毛片免费观看视频网站| 99久久久亚洲精品蜜臀av| 亚洲国产色片| 亚洲真实伦在线观看| 国产单亲对白刺激| 精品国产三级普通话版| 日本一本二区三区精品| 久久精品综合一区二区三区| 亚洲自拍偷在线| 久久精品国产亚洲网站| 欧美日韩在线观看h| 亚洲av不卡在线观看| 久久婷婷人人爽人人干人人爱| 高清在线视频一区二区三区 | 国产精品爽爽va在线观看网站| av在线亚洲专区| 日本免费a在线| 国产精品.久久久| 99riav亚洲国产免费| 国产精品一区二区性色av| 国产麻豆成人av免费视频| 自拍偷自拍亚洲精品老妇| 晚上一个人看的免费电影| 久久精品国产亚洲av香蕉五月| 久久人妻av系列| 床上黄色一级片| 国产精品不卡视频一区二区| 最近中文字幕高清免费大全6| 亚洲经典国产精华液单| 国产精品嫩草影院av在线观看| 国内少妇人妻偷人精品xxx网站| 人妻系列 视频| www日本黄色视频网| 内地一区二区视频在线| 在线免费观看的www视频| 亚洲欧美日韩高清在线视频| a级毛片免费高清观看在线播放| 亚洲性久久影院| 亚洲七黄色美女视频| 亚洲成人av在线免费| 日韩亚洲欧美综合| 亚洲精品日韩在线中文字幕 | 亚洲18禁久久av| 赤兔流量卡办理| 免费看美女性在线毛片视频| 久久人妻av系列| 国产男人的电影天堂91| 国产午夜精品一二区理论片| 两个人视频免费观看高清| 亚洲av二区三区四区| 高清毛片免费观看视频网站| 又粗又爽又猛毛片免费看| 国产精品伦人一区二区| 99久久精品热视频| 久久这里有精品视频免费| 成人漫画全彩无遮挡| 亚洲欧美日韩东京热| 一级毛片我不卡| 男的添女的下面高潮视频| 亚洲成a人片在线一区二区| 久久精品夜夜夜夜夜久久蜜豆| 日本一本二区三区精品| 最近视频中文字幕2019在线8| 九草在线视频观看| 久久午夜亚洲精品久久| 午夜福利在线在线| 校园春色视频在线观看| 国产男人的电影天堂91| 国产精品99久久久久久久久| 国产成人影院久久av| 中文字幕av成人在线电影| 蜜桃亚洲精品一区二区三区| 久久精品国产清高在天天线| 国产在线男女| 卡戴珊不雅视频在线播放| 国内精品美女久久久久久| 秋霞在线观看毛片| 久久精品国产亚洲av香蕉五月| 国产探花在线观看一区二区| 搡老妇女老女人老熟妇| 精品熟女少妇av免费看| 久久99热6这里只有精品| 男女做爰动态图高潮gif福利片| 国产色婷婷99| 嫩草影院新地址| 成人特级黄色片久久久久久久| 美女 人体艺术 gogo| 最新中文字幕久久久久| 99热这里只有是精品50| 乱人视频在线观看| 国产老妇女一区| 热99re8久久精品国产| 少妇丰满av| 99热这里只有是精品50| 亚洲av一区综合| 在线观看66精品国产| 亚洲国产色片| eeuss影院久久| 亚洲在线观看片| 青青草视频在线视频观看| 九九在线视频观看精品| 欧美色欧美亚洲另类二区| 久久99热这里只有精品18| 精品国内亚洲2022精品成人| 69人妻影院| 国产一级毛片七仙女欲春2| 日本黄色片子视频| 日韩欧美三级三区| 麻豆成人av视频| 午夜免费激情av| 久久99精品国语久久久| 亚洲天堂国产精品一区在线| 中国国产av一级| 小说图片视频综合网站| 国产日本99.免费观看| 国产精品久久久久久精品电影| 赤兔流量卡办理| 六月丁香七月| 波多野结衣高清无吗| 欧美高清性xxxxhd video| 又粗又硬又长又爽又黄的视频 | 亚洲最大成人av| 在线观看66精品国产| 一级黄片播放器| 日韩欧美一区二区三区在线观看| 18+在线观看网站| videossex国产| 国产极品天堂在线| 亚洲第一电影网av| 亚洲精品国产成人久久av| 波野结衣二区三区在线| 免费在线观看成人毛片| 久久精品国产清高在天天线| 91狼人影院| 给我免费播放毛片高清在线观看| 国产午夜精品久久久久久一区二区三区| 国产色爽女视频免费观看| 免费观看人在逋| 国产成人精品久久久久久| 国产精品伦人一区二区| 天堂√8在线中文| 黄色一级大片看看| 欧美xxxx性猛交bbbb| 精品一区二区免费观看| 国产伦理片在线播放av一区 | 亚洲国产精品国产精品| 亚洲美女搞黄在线观看| av免费观看日本| 国产乱人偷精品视频| 波多野结衣高清作品| 亚洲av成人精品一区久久| 精品一区二区免费观看| 久久久久久久亚洲中文字幕| 国产精品日韩av在线免费观看| 伦精品一区二区三区| 我要搜黄色片| 国产白丝娇喘喷水9色精品| 亚洲精品久久国产高清桃花| 亚洲第一区二区三区不卡| 亚洲四区av| 久久久久久九九精品二区国产| 在线免费观看不下载黄p国产| 亚洲中文字幕日韩| 99在线视频只有这里精品首页| 久久精品人妻少妇| 桃色一区二区三区在线观看| 村上凉子中文字幕在线| 国产中年淑女户外野战色| 欧美又色又爽又黄视频| 精品人妻熟女av久视频| 2022亚洲国产成人精品| 成人亚洲精品av一区二区| 99久久精品热视频| 国产高清三级在线| 日韩强制内射视频| 免费av不卡在线播放| 成人av在线播放网站| www日本黄色视频网| 成人美女网站在线观看视频| 中国国产av一级| 一个人看的www免费观看视频| 1024手机看黄色片| 亚洲欧美精品专区久久| 老司机福利观看| 亚洲欧洲国产日韩| 精品久久久久久久久久久久久| 亚洲自偷自拍三级| 成人欧美大片| 亚洲精品国产成人久久av| 搞女人的毛片| 久久6这里有精品| 国产精品一区二区在线观看99 | 国产伦精品一区二区三区四那| 欧美性感艳星| 看免费成人av毛片| 99riav亚洲国产免费| 22中文网久久字幕| 性欧美人与动物交配| 亚洲综合色惰| 国产精品电影一区二区三区| 欧美一区二区亚洲| av福利片在线观看| 午夜爱爱视频在线播放| 亚洲国产精品sss在线观看| 级片在线观看| 丰满乱子伦码专区| 亚洲欧美日韩东京热| 亚洲欧美清纯卡通| 亚洲av中文av极速乱| 午夜精品国产一区二区电影 | 如何舔出高潮|