• <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
    av在线老鸭窝| 亚洲内射少妇av| 性色av乱码一区二区三区2| 日韩有码中文字幕| 国产av麻豆久久久久久久| 九九在线视频观看精品| 国产私拍福利视频在线观看| 男插女下体视频免费在线播放| 中文字幕免费在线视频6| 啦啦啦韩国在线观看视频| 女同久久另类99精品国产91| 免费电影在线观看免费观看| 成年人黄色毛片网站| 国产在视频线在精品| 欧美不卡视频在线免费观看| 久久久成人免费电影| 日本a在线网址| 日本黄色片子视频| 黄色配什么色好看| www日本黄色视频网| 一级毛片久久久久久久久女| 国产日本99.免费观看| 美女被艹到高潮喷水动态| 国产极品精品免费视频能看的| 久久精品夜夜夜夜夜久久蜜豆| 最近中文字幕高清免费大全6 | 女人被狂操c到高潮| 亚洲午夜理论影院| 亚洲av成人不卡在线观看播放网| 欧美一区二区亚洲| 一级黄片播放器| 午夜视频国产福利| 国产午夜精品久久久久久一区二区三区 | av专区在线播放| 麻豆久久精品国产亚洲av| 黄片小视频在线播放| 久久午夜福利片| 内地一区二区视频在线| 变态另类成人亚洲欧美熟女| 亚洲精品色激情综合| 十八禁人妻一区二区| www.熟女人妻精品国产| 亚洲专区国产一区二区| 美女黄网站色视频| 午夜精品久久久久久毛片777| 国产色婷婷99| 亚洲精品在线观看二区| 最近中文字幕高清免费大全6 | 国产精品一区二区免费欧美| 亚洲专区中文字幕在线| 深爱激情五月婷婷| 99热这里只有是精品在线观看 | 色哟哟哟哟哟哟| 欧美xxxx黑人xx丫x性爽| 毛片一级片免费看久久久久 | 成人性生交大片免费视频hd| 国产精品人妻久久久久久| 大型黄色视频在线免费观看| 亚洲精品粉嫩美女一区| 国产成人欧美在线观看| 欧美日韩瑟瑟在线播放| 草草在线视频免费看| 老鸭窝网址在线观看| 免费黄网站久久成人精品 | 午夜精品在线福利| av专区在线播放| 亚洲精品亚洲一区二区| 亚洲无线在线观看| 99热6这里只有精品| 在线国产一区二区在线| 一级a爱片免费观看的视频| 嫩草影院新地址| 如何舔出高潮| 一区二区三区激情视频| 色噜噜av男人的天堂激情| 国产精品日韩av在线免费观看| 最后的刺客免费高清国语| 成人av在线播放网站| 天堂影院成人在线观看| 不卡一级毛片| 淫秽高清视频在线观看| 人妻久久中文字幕网| 成人永久免费在线观看视频| 免费一级毛片在线播放高清视频| av天堂在线播放| 久久精品人妻少妇| 91麻豆精品激情在线观看国产| 国产伦精品一区二区三区视频9| 久久99热这里只有精品18| 久久99热这里只有精品18| 亚洲美女黄片视频| 亚洲真实伦在线观看| 国产精品av视频在线免费观看| 床上黄色一级片| 久久6这里有精品| 给我免费播放毛片高清在线观看| 99久久精品国产亚洲精品| 在线播放国产精品三级| 久久久久免费精品人妻一区二区| 变态另类成人亚洲欧美熟女| 美女大奶头视频| 十八禁网站免费在线| 欧美zozozo另类| xxxwww97欧美| 国产精品电影一区二区三区| 国产免费男女视频| 最新中文字幕久久久久| 真人一进一出gif抽搐免费| 日日夜夜操网爽| 又粗又爽又猛毛片免费看| 听说在线观看完整版免费高清| 欧美黄色淫秽网站| 99视频精品全部免费 在线| 99热只有精品国产| 久久99热6这里只有精品| 欧美日韩瑟瑟在线播放| 精品国内亚洲2022精品成人| 国产野战对白在线观看| 亚洲精品久久国产高清桃花| 天堂av国产一区二区熟女人妻| 国产精品久久久久久亚洲av鲁大| 日韩av在线大香蕉| 窝窝影院91人妻| 日本五十路高清| 精品免费久久久久久久清纯| 国产 一区 欧美 日韩| 日本黄大片高清| 欧美日韩中文字幕国产精品一区二区三区| 精品日产1卡2卡| 国产视频一区二区在线看| 亚洲av成人av| 久久天躁狠狠躁夜夜2o2o| av女优亚洲男人天堂| 国产精品一区二区三区四区免费观看 | 一级a爱片免费观看的视频| 国产免费男女视频| 三级男女做爰猛烈吃奶摸视频| 成人国产一区最新在线观看| 久久中文看片网| 亚洲七黄色美女视频| 午夜激情福利司机影院| 看黄色毛片网站| 成人特级黄色片久久久久久久| 热99在线观看视频| 国产伦人伦偷精品视频| 亚洲人成伊人成综合网2020| 中文字幕久久专区| 欧美区成人在线视频| 色噜噜av男人的天堂激情| 亚洲av第一区精品v没综合| 在线免费观看不下载黄p国产 | 91九色精品人成在线观看| 亚洲精品一卡2卡三卡4卡5卡| 村上凉子中文字幕在线| 最新中文字幕久久久久| 日韩有码中文字幕| 日本三级黄在线观看| 午夜亚洲福利在线播放| 性色avwww在线观看| 国产淫片久久久久久久久 | 男人和女人高潮做爰伦理| 狠狠狠狠99中文字幕| 国产亚洲av嫩草精品影院| 毛片一级片免费看久久久久 | 亚洲黑人精品在线| 成人av在线播放网站| www.熟女人妻精品国产| 两性午夜刺激爽爽歪歪视频在线观看| 国产一区二区在线观看日韩| 久久久久久久精品吃奶| 毛片一级片免费看久久久久 | 亚洲最大成人av| 最后的刺客免费高清国语| 夜夜躁狠狠躁天天躁| 国产视频一区二区在线看| 五月玫瑰六月丁香| 久9热在线精品视频| 亚洲avbb在线观看| 91九色精品人成在线观看| 国产不卡一卡二| 一本久久中文字幕| 免费观看精品视频网站| 久久久久亚洲av毛片大全| 十八禁网站免费在线| 麻豆成人午夜福利视频| 中文字幕久久专区| 少妇熟女aⅴ在线视频| 色哟哟哟哟哟哟| 18禁在线播放成人免费| 99国产综合亚洲精品| 简卡轻食公司| 国产激情偷乱视频一区二区| 欧美日本亚洲视频在线播放| 久久久久久久精品吃奶| 国产又黄又爽又无遮挡在线| 精品福利观看| 久久久久久久久大av| 亚洲最大成人手机在线| 亚洲一区高清亚洲精品| 啦啦啦观看免费观看视频高清| 人人妻,人人澡人人爽秒播| 久久久久久久久大av| 最后的刺客免费高清国语| 中文在线观看免费www的网站| 激情在线观看视频在线高清| 国内揄拍国产精品人妻在线| 国产亚洲欧美在线一区二区| 亚洲成av人片免费观看| 99久久成人亚洲精品观看| 亚洲久久久久久中文字幕| 天堂av国产一区二区熟女人妻| av福利片在线观看| 国产主播在线观看一区二区| 国产黄色小视频在线观看| 久久九九热精品免费| 国产亚洲av嫩草精品影院| 成人精品一区二区免费| 2021天堂中文幕一二区在线观| 亚洲七黄色美女视频| 18禁黄网站禁片午夜丰满| 最新在线观看一区二区三区| 天堂av国产一区二区熟女人妻| 欧美+亚洲+日韩+国产| 亚洲精品粉嫩美女一区| 欧美xxxx性猛交bbbb| 日韩欧美在线乱码| 国产视频一区二区在线看| 九九久久精品国产亚洲av麻豆| 日韩欧美精品免费久久 | 国产主播在线观看一区二区| 欧美潮喷喷水| 亚洲中文字幕日韩| 欧美在线黄色| 国产精品,欧美在线| 精品熟女少妇八av免费久了| 亚洲专区国产一区二区| 亚洲av一区综合| 国产在线精品亚洲第一网站| 无遮挡黄片免费观看| 成年人黄色毛片网站| 黄色一级大片看看| 深爱激情五月婷婷| 看黄色毛片网站| 99久久精品热视频| 亚洲第一欧美日韩一区二区三区| av中文乱码字幕在线| 97热精品久久久久久| 亚洲无线观看免费| 淫秽高清视频在线观看| 成人一区二区视频在线观看| 色综合亚洲欧美另类图片| 97超级碰碰碰精品色视频在线观看| 国产亚洲精品久久久久久毛片| 亚洲国产欧美人成| 国产大屁股一区二区在线视频| 国产午夜精品久久久久久一区二区三区 | 国产精品综合久久久久久久免费| 99riav亚洲国产免费| 亚洲va日本ⅴa欧美va伊人久久| 变态另类丝袜制服| 国产又黄又爽又无遮挡在线| 中文字幕人成人乱码亚洲影| 国产单亲对白刺激| 午夜免费成人在线视频| 国产成人av教育| 久久精品91蜜桃| 久久久久久久精品吃奶| 日本熟妇午夜| 成人特级av手机在线观看| 午夜久久久久精精品| av国产免费在线观看| 中文字幕精品亚洲无线码一区| 一区二区三区激情视频| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品综合一区在线观看| 国产午夜精品久久久久久一区二区三区 | 最近中文字幕高清免费大全6 | 一级黄片播放器| 亚洲经典国产精华液单 | 国产色爽女视频免费观看| 国产探花极品一区二区| 一a级毛片在线观看| 97超视频在线观看视频| 简卡轻食公司| 高潮久久久久久久久久久不卡| ponron亚洲| 久久这里只有精品中国| 中亚洲国语对白在线视频| 国产大屁股一区二区在线视频| 成人特级黄色片久久久久久久| 色哟哟·www| 国产欧美日韩一区二区精品| 久久久久久久久中文| 一进一出抽搐gif免费好疼| 国产蜜桃级精品一区二区三区| 两个人视频免费观看高清| 亚洲精品乱码久久久v下载方式| 国产白丝娇喘喷水9色精品| 亚洲美女黄片视频| 一区二区三区四区激情视频 | 在线国产一区二区在线| 不卡一级毛片| 精华霜和精华液先用哪个| 岛国在线免费视频观看| 国产成人a区在线观看| 亚洲 国产 在线| 黄片小视频在线播放| 久久午夜亚洲精品久久| 亚洲无线观看免费| av欧美777| 久久午夜亚洲精品久久| 国产成+人综合+亚洲专区| 欧美+亚洲+日韩+国产| 好男人电影高清在线观看| 欧美午夜高清在线| 日本 欧美在线| 国产亚洲av嫩草精品影院| 欧美日韩亚洲国产一区二区在线观看| 国产精品1区2区在线观看.| 国产av在哪里看| 欧美一区二区亚洲| 男女那种视频在线观看| 好男人电影高清在线观看| 18禁裸乳无遮挡免费网站照片| 国产午夜精品论理片| 亚洲精品色激情综合| 毛片一级片免费看久久久久 | 久久久久九九精品影院| 天堂av国产一区二区熟女人妻| 欧美日韩综合久久久久久 | 神马国产精品三级电影在线观看| av在线天堂中文字幕| 日韩大尺度精品在线看网址| 好男人电影高清在线观看| 日本a在线网址| www.熟女人妻精品国产| 亚洲精品在线观看二区| 一夜夜www| 免费在线观看日本一区| 国产精品av视频在线免费观看| 最近在线观看免费完整版| 欧美激情久久久久久爽电影| 亚洲精品粉嫩美女一区| 中文在线观看免费www的网站| av天堂在线播放| 午夜精品在线福利| 少妇高潮的动态图| 黄色视频,在线免费观看| 国产在线精品亚洲第一网站| 18禁黄网站禁片免费观看直播| 亚洲国产色片| 免费看光身美女| 我的女老师完整版在线观看| 免费在线观看成人毛片| 日韩免费av在线播放| 国产黄色小视频在线观看| 国产精品野战在线观看| 日韩欧美一区二区三区在线观看| 午夜激情福利司机影院| 国产伦一二天堂av在线观看| 精品国产三级普通话版| 能在线免费观看的黄片| 亚洲av二区三区四区| 亚洲,欧美精品.| 亚洲av一区综合| 久久国产精品人妻蜜桃| 91av网一区二区| 午夜精品在线福利| 久久久久国产精品人妻aⅴ院| 九色成人免费人妻av| 国产精品伦人一区二区| 亚洲成人久久性| 18美女黄网站色大片免费观看| 麻豆一二三区av精品| 91av网一区二区| 成人无遮挡网站| 国产欧美日韩精品一区二区| 高清日韩中文字幕在线| 亚洲欧美激情综合另类| 亚洲成人免费电影在线观看| 久久精品综合一区二区三区| 国产精品永久免费网站| 日韩精品青青久久久久久| 一夜夜www| 日本一本二区三区精品| 热99在线观看视频| xxxwww97欧美| 欧美日韩黄片免| 九九热线精品视视频播放| 麻豆av噜噜一区二区三区| 国产精品亚洲一级av第二区| 偷拍熟女少妇极品色| 亚洲第一电影网av| 中文亚洲av片在线观看爽| 琪琪午夜伦伦电影理论片6080| 特大巨黑吊av在线直播| 动漫黄色视频在线观看| 俄罗斯特黄特色一大片| 免费一级毛片在线播放高清视频| 十八禁网站免费在线| 久久伊人香网站| 亚洲美女黄片视频| 简卡轻食公司| 女同久久另类99精品国产91| 精品免费久久久久久久清纯| 怎么达到女性高潮| 国产蜜桃级精品一区二区三区| 亚洲人成伊人成综合网2020| 夜夜看夜夜爽夜夜摸| a级毛片免费高清观看在线播放| 国产久久久一区二区三区| 亚洲自拍偷在线| www.www免费av| 少妇高潮的动态图| 久久久久久久午夜电影| 欧美国产日韩亚洲一区| 神马国产精品三级电影在线观看| 十八禁人妻一区二区| 国产免费男女视频| 久久国产乱子免费精品| 在线看三级毛片| 日韩欧美 国产精品| 一区二区三区免费毛片| 欧美激情久久久久久爽电影| 欧美一区二区亚洲| 窝窝影院91人妻| 日韩欧美三级三区| 首页视频小说图片口味搜索| 国产视频一区二区在线看| 一区二区三区四区激情视频 | 国产精品亚洲美女久久久| 国产成人福利小说| 丰满人妻熟妇乱又伦精品不卡| 人人妻,人人澡人人爽秒播| 欧美色视频一区免费| 欧美区成人在线视频| 午夜老司机福利剧场| 少妇熟女aⅴ在线视频| 岛国在线免费视频观看| 欧美乱色亚洲激情| 亚洲自偷自拍三级| 日韩欧美国产在线观看| 男人舔女人下体高潮全视频| 一区二区三区免费毛片| 免费一级毛片在线播放高清视频| 国产高清激情床上av| 免费在线观看成人毛片| 日本 欧美在线| 国产伦在线观看视频一区| 国产一区二区亚洲精品在线观看| 十八禁国产超污无遮挡网站| 国产蜜桃级精品一区二区三区| 色av中文字幕| 亚洲国产精品sss在线观看| aaaaa片日本免费| 亚洲 国产 在线| 国产真实乱freesex| 国产午夜福利久久久久久| 一本久久中文字幕| 一个人看的www免费观看视频| 欧美在线一区亚洲| 国产一级毛片七仙女欲春2| 在线免费观看不下载黄p国产 | 欧美一区二区亚洲| 免费av不卡在线播放| 一区二区三区免费毛片| 久久久久久九九精品二区国产| 亚洲一区高清亚洲精品| 欧美午夜高清在线| 中文字幕高清在线视频| 91久久精品国产一区二区成人| 精品国产亚洲在线| 99热这里只有是精品在线观看 | 99热精品在线国产| 亚洲av一区综合| 免费av不卡在线播放| av女优亚洲男人天堂| 亚洲成人精品中文字幕电影| 国产av一区在线观看免费| 亚洲精品粉嫩美女一区| 757午夜福利合集在线观看| 男人的好看免费观看在线视频| 嫩草影视91久久| 国产精品久久久久久久久免 | 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美三级三区| 91狼人影院| 99久久无色码亚洲精品果冻| 一区二区三区高清视频在线| 如何舔出高潮| 一二三四社区在线视频社区8| 91狼人影院| 一边摸一边抽搐一进一小说| 精品人妻偷拍中文字幕| 校园春色视频在线观看| 在线免费观看的www视频| 成人精品一区二区免费| 夜夜爽天天搞| 成人av一区二区三区在线看| 久久99热6这里只有精品| 色综合亚洲欧美另类图片| 午夜精品一区二区三区免费看| 久久精品影院6| 久久精品国产亚洲av涩爱 | 亚洲乱码一区二区免费版| 久久久久精品国产欧美久久久| 久久久久久久午夜电影| 国产乱人伦免费视频| 成人性生交大片免费视频hd| 欧美bdsm另类| 国产av麻豆久久久久久久| 欧美激情在线99| 每晚都被弄得嗷嗷叫到高潮| 精品日产1卡2卡| 长腿黑丝高跟| 国产精品乱码一区二三区的特点| 久久久色成人| 婷婷色综合大香蕉| 免费人成在线观看视频色| 狠狠狠狠99中文字幕| 精品一区二区三区视频在线观看免费| 首页视频小说图片口味搜索| 国产成年人精品一区二区| 亚洲,欧美,日韩| 国产亚洲精品久久久com| 亚洲乱码一区二区免费版| 最近最新中文字幕大全电影3| 亚洲性夜色夜夜综合| 国产视频一区二区在线看| 亚洲欧美激情综合另类| 宅男免费午夜| 一本综合久久免费| 亚洲欧美清纯卡通| 久久久久九九精品影院| 人人妻人人澡欧美一区二区| 真实男女啪啪啪动态图| 欧美一区二区亚洲| 999久久久精品免费观看国产| 国产69精品久久久久777片| 色吧在线观看| 搞女人的毛片| 日韩亚洲欧美综合| av国产免费在线观看| 免费人成在线观看视频色| 99国产精品一区二区蜜桃av| 久久久久久九九精品二区国产| 亚洲男人的天堂狠狠| 国产亚洲欧美在线一区二区| 欧美最新免费一区二区三区 | 成人av一区二区三区在线看| 成人亚洲精品av一区二区| 国产精品三级大全| 亚洲av五月六月丁香网| 99热6这里只有精品| 91麻豆精品激情在线观看国产| 久久久久免费精品人妻一区二区| 青草久久国产| 3wmmmm亚洲av在线观看| 午夜视频国产福利| www.www免费av| 精品久久久久久,| 人妻丰满熟妇av一区二区三区| 精品免费久久久久久久清纯| 成人av一区二区三区在线看| 内地一区二区视频在线| 久久久精品欧美日韩精品| 一区二区三区激情视频| 99久国产av精品| 精品久久国产蜜桃| 麻豆成人av在线观看| 亚洲一区高清亚洲精品| АⅤ资源中文在线天堂| 婷婷精品国产亚洲av| 国产视频内射| 波多野结衣高清作品| 少妇高潮的动态图| 亚洲人成伊人成综合网2020| 亚州av有码| 婷婷亚洲欧美| 老司机午夜十八禁免费视频| 欧美黑人欧美精品刺激| 又黄又爽又免费观看的视频| 欧美xxxx性猛交bbbb| 久久中文看片网| 757午夜福利合集在线观看| 无人区码免费观看不卡| 精品乱码久久久久久99久播| 免费看光身美女| bbb黄色大片| 久久久久久久久大av| 伊人久久精品亚洲午夜| 一进一出抽搐动态| 亚洲综合色惰| 日本熟妇午夜| 在线免费观看不下载黄p国产 | 欧美色视频一区免费| 亚洲美女视频黄频| 国产熟女xx| 别揉我奶头~嗯~啊~动态视频| 亚洲美女视频黄频| 国产精品亚洲美女久久久| 简卡轻食公司| 国产私拍福利视频在线观看| 亚洲欧美日韩东京热| 51国产日韩欧美| 亚洲av美国av| 国产乱人伦免费视频| 欧美黄色淫秽网站| 女生性感内裤真人,穿戴方法视频| 免费av不卡在线播放| 国产单亲对白刺激| 国产一区二区在线av高清观看| 99视频精品全部免费 在线| 色综合欧美亚洲国产小说| 欧美一区二区国产精品久久精品|