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

    河北小五臺山國家級自然保護(hù)區(qū)森林群落與環(huán)境的關(guān)系

    2017-07-20 13:11:03白曉航張金屯王云泉SehrishSadia
    生態(tài)學(xué)報 2017年11期
    關(guān)鍵詞:草本喬木灌木

    白曉航,張金屯,曹 科,王云泉,Sehrish Sadia,曹 格

    北京師范大學(xué)生命科學(xué)學(xué)院, 北京 100875

    ?

    河北小五臺山國家級自然保護(hù)區(qū)森林群落與環(huán)境的關(guān)系

    白曉航,張金屯*,曹 科,王云泉,Sehrish Sadia,曹 格

    北京師范大學(xué)生命科學(xué)學(xué)院, 北京 100875

    植物與環(huán)境之間的關(guān)系是一個復(fù)雜的演變過程,運(yùn)用數(shù)量生態(tài)學(xué)方法探討森林群落的物種組成、種群的生態(tài)特征、不同植物群落與環(huán)境之間的關(guān)系,有助于保護(hù)該區(qū)森林群落的穩(wěn)定性和生物多樣性。根據(jù)148個森林群落樣方數(shù)據(jù),選用雙向指示種分析(TWINSPAN)和典范對應(yīng)分析(CCA)方法,對河北小五臺山國家級自然保護(hù)區(qū)森林群落進(jìn)行分類和排序研究。結(jié)果表明:(1)TWINSPAN將該區(qū)的森林群落分為20個類型;(2)CCA排序結(jié)果較好地反映出群落分布格局與環(huán)境梯度的關(guān)系,各個森林群落類型在前兩軸分異明顯,在11個環(huán)境因子中,海拔、坡位、凋落層厚度、土壤導(dǎo)電率、土壤溫度、土壤濕度、土壤厚度和干擾程度這8個環(huán)境因子對森林群落的分布起較大的作用,影響森林群落的分布格局,形成不同的植被類型。(3)喬木層優(yōu)勢種的CCA二維排序圖所揭示的環(huán)境梯度與群落類型的分布有很大的相似性;t值雙序圖闡明了海拔、凋落層厚度、土壤溫度、干擾程度等環(huán)境因子對森林群落喬木層優(yōu)勢種有著重要影響。采用TWINSPAN 分類與CCA排序的方法,較好地解釋了森林群落與環(huán)境因子的關(guān)系,為小五臺山地區(qū)森林生態(tài)系統(tǒng)的科學(xué)管理和保護(hù)提供了理論依據(jù),研究結(jié)果也為同類地區(qū)森林生態(tài)系統(tǒng)研究及保護(hù)提供參考和借鑒。

    小五臺山國家級自然保護(hù)區(qū);森林群落;TWINSPAN;CCA;環(huán)境梯度

    森林群落結(jié)構(gòu)反映了植物群落對環(huán)境的適應(yīng)、動態(tài)和機(jī)能,群落類型不同,結(jié)構(gòu)也不同。分類能夠在一定程度上揭示植物群落類型的形成、發(fā)展及其與周圍環(huán)境的關(guān)系,是確定植被間斷性的重要方法。等級分類技術(shù)中的雙向指示種分析法(Two-way indicator species analysis,TWINSPAN)分析時可同時完成樣方和物種的分類,方法精細(xì),適合不同尺度的群落分類,在植被數(shù)量分類方法中占有重要影響[1-2]。排序是研究植被連續(xù)變化的方法,是指用數(shù)學(xué)的方法將樣方或植物種排列在一定的空間,使得排序軸反映一定的生態(tài)梯度,從而解釋植物物種、植物群落的分布與環(huán)境之間的關(guān)系[3-4]。典范對應(yīng)分析(Canonical correspondence analysis,CCA)在研究過程中選用植物物種數(shù)據(jù)和環(huán)境因子數(shù)據(jù),每一步的計算結(jié)果都與環(huán)境因子進(jìn)行回歸,揭示環(huán)境因子對物種組成的影響,從而更好地反映群落與環(huán)境的關(guān)系,對植被類型和環(huán)境因子間的關(guān)系做出合理解釋[5-6]。

    保護(hù)區(qū)有著相對原始的自然環(huán)境和較為完整的生態(tài)系統(tǒng),研究自然保護(hù)區(qū)內(nèi)森林群落類型及其與自然環(huán)境的生態(tài)關(guān)系,對保證自然保護(hù)區(qū)內(nèi)植物多樣性和珍稀野生資源的可持續(xù)發(fā)展有著重要的指導(dǎo)意義[7]。河北小五臺山國家自然保護(hù)區(qū)具有華北地區(qū)保存較為完好的暖溫帶森林生態(tài)系統(tǒng),屬暖溫帶落葉闊葉林植被帶,區(qū)內(nèi)生物多樣性豐富,區(qū)系成分復(fù)雜,具有重要的保護(hù)價值。近年來,有諸多學(xué)者從多方面對該地區(qū)植物進(jìn)行研究,如對小五臺山森林生態(tài)系統(tǒng)主要因子間量化關(guān)系研究、海拔梯度上青楊的種群結(jié)構(gòu)和空間分布的研究、天然白樺林草本多樣性影響因素的研究等[8- 10],但運(yùn)用數(shù)量生態(tài)學(xué)的方法對森林植物群落數(shù)量分類及其與環(huán)境之間關(guān)系的研究鮮有報道。本文以河北小五臺山國家級自然保護(hù)區(qū)森林群落為研究對象,運(yùn)用雙向指示種分析(TWINSPAN)和典范對應(yīng)分析(CCA)方法對其進(jìn)行分類和排序,分析該區(qū)植物群落類型,揭示分布于不同群落類型的植物物種與環(huán)境之間的關(guān)系,從而為保護(hù)和管理該地區(qū)森林植物資源提供科學(xué)依據(jù)。

    1 材料和方法

    1.1 研究地區(qū)概況

    河北小五臺山地處太行山、燕山、恒山交匯處,位于燕山—太行山山系區(qū),為太行山主峰,山體峻峭,擁有134座海拔2000米以上的山峰。自然保護(hù)區(qū)位于河北省西北部,張家口市蔚縣和涿鹿縣境內(nèi),地理坐標(biāo)為114°47′8″—115°28′56″(E),39°50′41″—40°6′30″(N), 東西長約60km,南北寬約28km,總面積21833hm2。該區(qū)屬暖溫帶大陸季風(fēng)型山地氣候,具有四季分明、雨熱同期、冬長夏短、夏季晝夜溫差大、隨海拔變化差異顯著等特點(diǎn)。年均溫6.4℃,1月平均氣溫-12.3℃,7月平均氣溫22.1℃,年降水量400—700mm,無霜期110—140d,氣候垂直分帶性比較明顯,集中表現(xiàn)出北方山區(qū)氣候特點(diǎn)[11]。植被類型豐富,垂直帶普明顯,有次生灌草帶、闊葉林帶、針闊混交林帶、亞高山灌叢帶、亞高山草甸帶等,基點(diǎn)高度1200m左右,1200m以下多為農(nóng)田、裸地等。土壤主要有褐土、棕壤、亞高山灌叢草甸土等類型[8]。

    1.2 樣方設(shè)置與調(diào)查

    2015年7—9月在河北小五臺山國家級自然保護(hù)區(qū)內(nèi)進(jìn)行全面踏查,從海拔1300m到2600m分別從東、西、南三臺沿海拔梯度取樣(小五臺山地區(qū)基點(diǎn)高度為1200m,2600m以上全部為亞高山草甸,調(diào)查林地包括灌木林和喬木林),海拔每升高100m三臺分別設(shè)一條樣帶,在每條樣帶上隨機(jī)設(shè)置4個面積10m×10m的樣方[12]。三臺共計37條樣帶,東臺14條、西臺11條、南臺12條,共計獲得148個森林樣方。在10m×10m的樣方內(nèi)調(diào)查并記錄喬木的株數(shù)、高度、胸徑、蓋度、冠幅等數(shù)據(jù),記錄樣方內(nèi)全部灌木和草本的蓋度和高度等數(shù)據(jù)。用手持GPS和地質(zhì)羅盤測量各個樣地的經(jīng)緯度、海拔、坡度、坡向和坡位,用土釬和鋼卷尺測定土壤深度和枯枝落葉層厚度,用土壤原位pH計(IQ150 pH meter)測定土壤pH,用土壤溫濕度儀(Moisture Meter type HH2)測定土壤溫度、濕度、導(dǎo)電率,同時記錄各個群落的生態(tài)學(xué)特征及樣地的生境條件。

    1.3 數(shù)據(jù)處理

    1.3.1 重要值數(shù)據(jù)

    148個樣方中共記錄到392種植物,數(shù)據(jù)分析采用重要值作為綜合指標(biāo)來反映種的特征,計算公式如下:

    喬木重要值=(相對多度+相對優(yōu)勢度+相對高度)/ 3

    灌木、草本重要值=(相對蓋度+相對高度)/ 2

    1.3.2 環(huán)境數(shù)據(jù)

    環(huán)境數(shù)據(jù)中海拔高度、坡度、凋落層厚度、土壤深度、土壤濕度、土壤溫度、土壤導(dǎo)電率、土壤pH以實(shí)際觀測記錄值表示,坡位、坡向、干擾程度采用數(shù)字等級表示。坡位分為五個等級:谷底為1級,下坡位為2級,中坡位為3級,上坡位為4級,山頂為5級。坡向以數(shù)字等級表示,以北為0°,將坡向按照45°的夾角,順時針方向旋轉(zhuǎn)分為8個坡向等級,以數(shù)字表示各等級:1表示北坡(0°—22.5°,337.5°—360°),2表示東北坡(22.5°—67.5°),3表示西北坡(292.5°—337.5°),4表示東坡(67.5°—112.5°),5表示西坡(247.5°—292.5°),6表示東南坡(112.5°—157.5°),7表示西南坡(202.5°—247.5°),8表示南坡(157.5°—202.5°),數(shù)字越大,表示坡向越向陽[13]。干擾情況分為4個等級:無干擾為1級,輕微干擾為2級,中度干擾為3級,強(qiáng)度干擾為4級。共計11種環(huán)境因子。

    1.4 數(shù)量分類與排序

    采用TWINSPAN分類方法,用WinTWINS for Window Version 2.3軟件分析[14- 15]。植被物種數(shù)據(jù)是397個種(植物種類392種,有5種喬木幼苗在灌木層中起作用,算做假種)在148個樣方中的重要值構(gòu)成,以397×148維重要值矩陣進(jìn)行數(shù)量分類。對群落預(yù)先進(jìn)行DCA排序,結(jié)果表明第一軸梯度長度6.371>3SD(Standard deviation, SD),基于單峰模型的排序方法比較合適,因此,選擇CCA排序方法[16]。用CANOCO Version 4.5軟件對397×148維重要值矩陣和11×148維環(huán)境因子矩陣進(jìn)行CCA排序,并制作樣方-環(huán)境因子二維排序圖、喬木優(yōu)勢種-環(huán)境因子二維排序圖、喬木優(yōu)勢種-環(huán)境t值雙序圖。

    2 結(jié)果與分析

    2.1 TWINSPAN分類

    對森林群落的148個樣方進(jìn)行TWINSPAN數(shù)量分類,TWINSPAN樹狀分類如圖1所示。結(jié)合調(diào)查區(qū)域的實(shí)際生態(tài)意義,采用第七級的劃分結(jié)果,將保護(hù)區(qū)的森林群落劃分為20個群落類型。依據(jù)中國植被的分類原則和系統(tǒng),結(jié)合各層優(yōu)勢種和TWINSPAN劃分的指示種命名群落類型。樣方分類結(jié)果見圖1,TWINSPAN所分類的20個群落類型,其群落結(jié)構(gòu)、物種組成差異明顯,主要環(huán)境特征見表1。

    表1 小五臺山國家級自然保護(hù)區(qū)森林群落的主要特征描述

    圖1 小五臺山國家級自然保護(hù)區(qū)森林群落TWINSPAN分類樹狀圖(圖中數(shù)字表示樣方編號)Fig.1 Dendrogram of the TWINSPAN classification of forest communities in the Xiaowutai Monutain National Nature Reserve (The numbers in the figure indicate quadrats serial numbers)

    20個森林群落類型的群落特征及其主要物種組成成分:

    I 白樺-六道木+榛-糙蘇+瓣蕊唐松草群落(Comm.Betulaplatyphylla-Abeliabiflora+Corylusheterophylla-Phlomisumbrosa+Thalictrumpetaloideum)。該群落主要分布于海拔1900—2000m的南臺和2000—2100m的西臺,主要分布在西臺的西坡、南臺的西南坡,坡度35—50°,處于中、上坡位,棕壤,枯枝落葉層略厚。群落總蓋度為60—85%,喬木層蓋度為48—65%,灌木層蓋度為22%—35%,草本層蓋度為12%—21%。喬木層以白樺(Betulaplatyphylla)為主要優(yōu)勢種,伴生有少量棘皮樺(Betuladahurica)、北京花楸(Sorbusdiscolor)和碩樺(Betulacostata)等;灌木層優(yōu)勢種為六道木(Abeliabiflora)、榛(Corylusheterophylla),伴生有金花忍冬(Lonicerachrysantha)、圓葉鼠李(Rhamnusglobosa)、刺五加(Acanthopanaxsenticosus)等;草本層中糙蘇(Phlomisumbrosa)、瓣蕊唐松草(Thalictrumpetaloideum)為優(yōu)勢種,伴生有歪頭菜(Viciaunijuga)、山尖子(Paraseneciohastata)、短毛獨(dú)活(Heracleummoellendorffii)、藜蘆(Veratrumnigrum)、細(xì)葉沙參(Adenophorapaniculata)、綿團(tuán)鐵線蓮(Clematishexapetala)、黃精(Polygonatumsibiricum)、草問荊(Equisetumpratense)、興安老鸛草(Geraniummaximowiczii)等。

    II紅樺-六道木+東陵繡球-興安升麻+披針薹草群落(Comm.Betulaalbosinensis-Abeliabiflora+Hydrangeabretschneideri-Cimicifugadahurica+Carexlancifolia)。該群落主要分布于海拔1818—2100m的南臺,多處東北坡、東坡,坡度31—42°,中、上坡位,棕壤,枯枝落葉層較厚。群落總蓋度為75—91%,喬木層蓋度為52—75%,灌木層蓋度為13—23%,草本層蓋度為30—49%。喬木層以紅樺(Betulaalbosinensis)為主要優(yōu)勢種,伴生有碩樺苗、水榆花楸(Sorbusalnifolia)。灌木層優(yōu)勢種為東陵繡球(Hydrangeabretschneideri)、六道木,伴生有少量中國黃花柳(Salixsinica)、金花忍冬、河北柳(Salixtaishanensisvar.hebeinica)、柳葉鼠李(Rhamnuserythroxylon)、瘤糖茶藨子(Ribeshimalensevar.verruculosum)等;草本層優(yōu)勢種興安升麻(Cimicifugadahurica)、披針薹草(Carexlancifolia),伴生有少量的東亞唐松草(Thalictrumminusvar.hypoleucum)、歪頭菜、山尖子、串鈴草(Phlomismongolica)、野青茅(Deyeuxiaarundinacea)、短花梗黃耆(Astragalushancockii)、蹄葉橐吾(Ligulariafischeri)、牛扁(Aconitumbarbatumvar.puberulum)等。

    III紅樺-土莊繡線菊-東亞唐松草+青綠薹草群落(Comm.Betulaalbosinensis-Spiraeapubescens-Thalictrumminus.var.hypoleucum+Carexbreviculmis)。該群落主要分布于海拔1907—2218m的南臺,東北坡、西坡,坡度35—40°,上坡位,棕壤,枯枝落葉層略厚。群落總蓋度為80—95%,喬木層蓋度為70—89%,灌木層蓋度為5—12%,草本層蓋度為28—32%。喬木層以紅樺為主要優(yōu)勢種,伴生有華北落葉松(Larixprincipis-rupprechtii);灌木層優(yōu)勢種為土莊繡線菊(Larixprincipis-rupprechtii),伴生有瘤糖茶藨子、東陵繡球、柳葉鼠李等;草本層優(yōu)勢種為東亞唐松草、青綠薹草(Carexbreviculmis),樣方內(nèi)還有興安升麻、山尖子、狹葉蕁麻(Urticaangustifolia)、蛇莓(Duchesneaindica)、寬葉薹草(Carexsiderosticta)、黑柴胡(Bupleurumsmithii)、柳葉蒿(Artemisiaintegrifolia)、陰地堇菜(Violayezoensis)、小紅菊、歪頭菜。

    IV紅樺-六道木+四川忍冬-黃精+披針薹草群落(Comm.Betulaalbosinensis-Abeliabiflora+Loniceraszechuanica-Polygonatumsibiricum+Carexlancifolia)。該群落主要分布于海拔1700—1800m的南臺,東北坡,坡度35—39°,下坡位,棕壤,枯枝落葉層不均勻,部分地段較厚。群落總蓋度為83%—90%,喬木層蓋度為68%—85%,灌木層蓋度為 13%—65%,草本層蓋度為26%—35%。喬木層以紅樺為主要優(yōu)勢種,伴生有白樺、遼東櫟(Quercuswutaishanica)等;灌木層優(yōu)勢種為四川忍冬(Loniceraszechuanica)、六道木,伴生有少量的刺五加(Acanthopanaxsenticosus)、毛榛(Corylusmandshurica)、土莊繡線菊、凍綠(Rhamnusutilis)等;草本層中黃精和披針薹草分布較多,還有少量青綠薹草、興安升麻、白蓮蒿(Artemisiasacrorum)、三脈紫菀(Asterageratoides)、山柳菊(Hieraciumumbellatum)、纈草(Valerianaofficinalis)、蒼術(shù)(AtractylodesLancea)、豬殃殃(Galiumaparinevar.tenerum)、北烏頭(Aconitumkusnezoffii)、蛇床(Cnidiummonnieri)、大齒山芹(Ostericumgrosseserratum)、玉竹(Polygonatumodoratum)、白背鐵線蕨(Adiantumdavidii)等。

    VI白樺+紅樺-紅丁香-青綠薹草群落(Comm.Betulaplatyphylla+Betulaalbosinensis-Syringavillosa-Carexbreviculmis)。該群落主要分布于海拔 2100m左右的西臺,西坡、西南坡、南坡,坡度45—58°,上坡位,棕壤,枯枝落葉層較薄。群落總蓋度為92%,喬木層蓋度為40%—55%,灌木層蓋度為15%—25%,草本層蓋度為45%。喬木層以白樺、紅樺為主要優(yōu)勢種,伴生有少量的小青楊;灌木層優(yōu)勢種為紅丁香(Syringavillosa),伴生有六道木、中國黃花柳、金花忍冬、土莊繡線菊;草本層中青綠薹草分布較多,還有少量地榆(Sanguisorbaofficinalis)、毛建草(Dracocephalumrupestre)、莓葉委陵菜(Potentillafragarioides)、黃芩(Scutellariabaicalensis)、蛇床、興安老鸛草、北烏頭等。

    VII白樺-毛榛+照山白-短毛獨(dú)活+珠芽蓼群落(Comm.Betulaplatyphylla-Corylusmandshurica+Rhododendronmicranthum-Heracleummoellendorffii+Polygonumviviparum)。該群落主要分布于海拔1889—2020m的西臺,西北坡,坡度 40—52°,上坡位,枯枝落葉層較薄。群落總蓋度為85%—93%,喬木層蓋度為 60%—70%,灌木層蓋度為25%—50%,草本層蓋度為 14%—20%。喬木層以白樺為主要優(yōu)勢種,伴生有少量小青楊(Populuspseudo-simonii)、北京花楸;灌木層優(yōu)勢種為毛榛、照山白,伴生有美薔薇(Rosabella)、興安胡枝子(Lespedezadavurica)等;草本層中短毛獨(dú)活(Heracleummoellendorffii)、珠芽蓼(Polygonumviviparum)分布較多,還有少量糙蘇(Phlomisumbrosa)、短尾鐵線蓮(Clematisbrevicaudata)、蚊子草(Filipendulapalmata)、瓣蕊唐松草、小玉竹(Polygonatumhumile)、多歧沙參(Adenophorawawreana)、舞鶴草(Maianthemumbifolium)、落新婦(Astilbechinensis)、紫菀(Astertataricus)、牛扁、藍(lán)花棘豆(Oxytropiscoerulea)、遼藁本(Ligusticumjeholense)、草問荊等。

    VII紅樺+白樺-土莊繡線菊+小葉鼠李-莓葉委陵菜+寬葉薹草群落(Comm.Betulaalbosinensis+ Betulaplatyphylla-Spiraeapubescens+Rhamnusparvifolia-Potentillafragarioides+Carexsiderosticta)。該群落主要分布于海拔1705—2137m的東臺,西北坡、東北坡、東坡,坡度12—55°,中、上坡位,棕壤,枯枝落葉層5—10cm。群落總蓋度為80%—97%,喬木層蓋度為35%—70%,灌木層蓋度為45%—60%,草本層蓋度為50%—63%。喬木層以紅樺、白樺為主要優(yōu)勢種,伴生有華北落葉松,小葉楊(Populussimonii)、青扦(Piceawilsonii)等;灌木層優(yōu)勢種為土莊繡線菊、小葉鼠李(Rhamnusparvifolia),伴生有美薔薇、蒙古莢蒾(Viburnummongolicum)、西北栒子(Cotoneasterzabelii)、細(xì)葉小檗(Berberispoiretii)、東北茶藨子(Ribesmandshuricum)、無梗五加(Acanthopanaxsessiliflorus);草本層中莓葉委陵菜、寬葉薹草分布較多,伴生有水金鳳(Impatiensnoli-tangere)、鵝觀草(Roegneriakamoji)、金蓮花(Trolliuschinensis)、南牡蒿(Artemisiaeriopoda)、地榆、草珠黃耆(Astragaluscapillipes)、老鸛草、野青茅、牛扁、兩色烏頭(Aconitumalboviolaceum)、多歧沙參、寬葉蕁麻、北重樓(Parisverticillata)、黃精、蔭生鼠尾草(Salviaumbratica)、蛇莓、黃花列當(dāng)(Orobanchepycnostachya)、舞鶴草、白花碎米薺(Cardamineleucantha)、華北耬斗菜(Aquilegiayabeana)、粗根老鸛草(Geraniumdahuricum)、大葉野豌豆(Viciapseudorobus)、細(xì)葉鴉蔥(Scorzonerapusilla)等,草本種類非常豐富。

    VIII白樺+油松-鼠李+山刺玫-串鈴草+白蓮蒿群落(Comm.Betulaplatyphylla+Pinustabuliformis-Rhamnusdavurica+Rosadavurica-Phlomismongolica+Artemisiasacrorum)。該群落主要分布于海拔1300—1700m的東臺,東北坡、北坡,坡度15—18°,谷底、下坡位,棕壤,枯枝落葉層極薄,干擾較強(qiáng)。群落總蓋度為83—95%,喬木層蓋度為40—55%,灌木層蓋度為25—65%,草本層蓋度為38—65%。喬木層以白樺和油松(Pinustabuliformis)為主要優(yōu)勢種,伴生有山楊(Populusdavidiana)、榆樹(Ulmuspumila)、大葉白蠟(Fraxinusrhynchophylla)等;灌木層優(yōu)勢種為鼠李(Rhamnusdavurica)和山刺玫(Rosadavurica),伴生有牛疊肚(Rubuscrataegifolius)、胡枝子(Lespedezabicolor)、虎榛子(Ostryopsisdavidiana)、杭子梢(Campylotropismacrocarpa)、河北柳、美薔薇、小葉鼠李、黑彈樹(Celtisbungeana)、藍(lán)靛果忍冬、灰栒子(Cotoneasteracutifolius)等,灌木層植物種類較為豐富;草本層中串鈴草和白蓮蒿(Artemisiasacrorum)分布較多,還有少量歧莖蒿(Artemisiaigniaria)、細(xì)葉豬殃殃、龍芽草(Agrimoniapilosa)、朝天委陵菜(Potentillasupina)、香薷(Elsholtziaciliata)、白芷(Angelicadahurica)、華北復(fù)盆子(Rubusidaeusvar.borealisinensis)、鵝觀草、地榆、野青茅、景天三七(Sedumaizoon)、密花香薷(Elsholtziadensa)、小紅菊、北千里光(Seneciodubitabilis)、狹葉黃精(Polygonatumstenophyllum)、細(xì)葉沙參、飛廉(Carduusnutans)、拳參(Polygonumbistorta)、藜蘆、歪頭菜、五?;?Adoxamoschatellina)、牛扁、短尾鐵線蓮(Clematisbrevicaudata)、茜草(Rubiacordifolia)、穿龍薯蕷(Dioscoreanipponica)、遼藁本等。

    IX油松-美薔薇-老鸛草+寬葉薹草群落(Comm.Pinustabuliformis-Rosabella-Geraniumwilfordii+Carexsiderosticta)。該群落主要分布于海拔1409—1598m的西臺,北坡、西北坡,坡度27—37°,下坡位,棕壤,枯枝落葉層較薄,輕微干擾,有金錢豹出沒。群落總蓋度為92%—95%,喬木層蓋度為20%—50%,灌木層蓋度為40%—75%,草本層蓋度為45%—70%。喬木層以油松為主要優(yōu)勢種,伴生有少量華北落葉松;灌木層優(yōu)勢種為美薔薇,伴生有毛櫻桃(Cerasustomentosa)、毛榛、紅瑞木(Swidaalba)、小葉鼠李、灰栒子、三裂繡線菊等;草本層中老鸛草和寬葉薹草分布較多,還有少量鴉蔥(Scorzoneraaustriaca)、芹葉鐵線蓮(Clematisaethusifolia)、多歧沙參、兩型豆(Amphicarpaeaedgworthii)、歧莖蒿、蛇莓、牛扁、寬葉蕁麻(Urticalaetevirens)、龍芽草等。

    X山楊+白樺-中國黃花柳-五?;?+青綠薹草群落(Comm.Populusdavidiana+Betulaplatyphylla-Salixsinica-Adoxamoschatellina+Carexbreviculmis)。該群落主要分布于海拔1592—1838m的西臺,西北坡、西坡、西南坡,坡度29—50°,中、上坡位,棕壤,枯枝落葉層較薄,輕微干擾。群落總蓋度為65%—72%,喬木層蓋度為25%—40%,灌木層蓋度為10%—30%,草本層蓋度為40—55%。喬木層以山楊和白樺為主要優(yōu)勢種,伴生有少量油松、遼東櫟、甘肅山楂等;灌木層優(yōu)勢種為中國黃花柳,伴生有小葉鼠李、三裂繡線菊等;草本層中五?;?Adoxamoschatellina)、青綠薹草分布較多,還有少量小紅菊、藿香(Agastacherugosa)、玉竹、地榆、陰地堇菜、莓葉委陵菜、牛膝菊(Galinsogaparviflora)、并頭黃芩(Scutellariascordifolia)、飛廉、蒲公英(Taraxacummongolicum)、糖芥(Erysimumbungei)、中華小苦荬(Ixeridiumchinense)、鼠掌老鸛草(Geraniumsibiricum)、瓣蕊唐松草、兩型豆等。

    XI榆樹-小葉鼠李+榆葉梅-歧莖蒿群落(Comm.Ulmuspumila-Rhamnusparvifolia+Amygdalustriloba-Artemisiaigniaria)。該群落主要分布于海拔 1300—1500 m的東臺,西北坡,坡度10—25°,下坡位,山地褐土,枯枝落葉層極薄,輕微干擾。群落總蓋度為50%—90%,喬木層蓋度為30%,灌木層蓋度為36%—90%,草本層蓋度為 29%—84%。喬木層以榆樹為主要優(yōu)勢種,伴生有少量小葉楊、皂柳(Salixwallichiana)、油松等;灌木層優(yōu)勢種為小葉鼠李和榆葉梅,伴生有短梗胡枝子、小花溲疏(Deutziaparviflora)、刺果茶藨子、沙棘(Hippophaerhamnoides)、蒙古繡線菊、虎榛子、美薔薇、太平花(Philadelphuspekinensis)、西北栒子等,灌木層植物種類豐富;草本層中歧莖蒿分布較多,還有少量軸藜(Axyrisamaranthoides)、香青蘭(Dracocephalummoldavica)、狗尾草(Setariaviridis)、景天三七、短尾鐵線蓮、耬斗菜(Aquilegiaviridiflora)、小紅菊、石竹(Dianthuschinensis)、糖芥、龍芽草、野青茅、細(xì)葉婆婆納(Veronicalinariifolia)、大葉野豌豆、巴天酸模(Rumexpatientia)等。

    XII油松+杏-毛櫻桃-白蓮蒿+寬葉薹草群落(Comm.Pinustabuliformis+ Armeniacavulgaris-Cerasustomentosa-Artemisiasacrorum+Carexsiderosticta)。該群落主要分布于海拔1300—1600m的西臺,西坡、西北坡,坡度30—45°,下坡位,棕壤,枯枝落葉層極薄,輕微干擾。群落總蓋度為88—95%,喬木層蓋度為30%—40%,灌木層蓋度為35%—75%,草本層蓋度為18%—39%。喬木層以油松和杏為主要優(yōu)勢種,伴生有少量山楊、榆樹、小葉楊;灌木層優(yōu)勢種為毛櫻桃,伴生有甘草、毛榛、達(dá)烏里胡枝子、小葉鼠李、沙棘、北京丁香;草本層中分白蓮蒿和寬葉薹草布較多,伴生有草木犀(Melilotusofficinalis)、黃耆(Astragalusmembranaceus)、野艾蒿(Artemisialavandulaefolia)、大油芒(Spodiopogonsibiricus)、大葉野豌豆、早熟禾(Poaannua)、鵝觀草、狼針草(Stipabaicalensis)、糙葉敗醬(Patriniarupestrissubsp.scabra)、刺兒菜(Cirsiumsetosum)、紫菀、灰綠藜(Chenopodiumglaucum)、豬殃殃、野韭、小紅菊、芹葉鐵線蓮、蒲公英、朝天委陵菜、藿香等。

    XIII杠柳-豬殃殃群落(Comm.Periplocasepium-Galiumaparinevar.tenerum)。該群落主要分布于海拔1290—1388m左右的南臺,南坡,坡度20—35°,谷底,山地褐土,近乎無枯枝落葉層,人為干擾較強(qiáng)。群落總蓋度為90%—92%,灌木層蓋度為88%,草本層蓋度為10%。無喬木層,灌木層優(yōu)勢種為杠柳(Periplocasepium),伴生有水栒子(Cotoneastermultiflorus)、巧玲花(Syringapubescens)等;草本層中豬殃殃分布較多,還有白蓮蒿、狗尾草、小紅菊、車前(Plantagoasiatica)、地榆、菊葉委陵菜(Potentillatanacetifolia)等。

    XIV色木槭+遼東櫟-虎榛子+巧玲花-龍芽草群落(Comm.Acermomo+Quercuswutaishanica—Ostryopsisdavidiana+Syringapubescens—Agrimoniapilosa)。該群落主要分布于海拔1400—1700m的南臺,南坡、西南坡,坡度35 —43°,下、中坡位,山地褐土,凋落層厚度較薄,人為干擾較大。群落總蓋度為92%—94%,喬木層蓋度為59%—70%,灌木層蓋度為65%—80%,草本層蓋度為12%—40%。喬木層以色木槭和遼東櫟為主要優(yōu)勢種,伴生有少量榆樹、野山楂、遼楊等;灌木層優(yōu)勢種為虎榛子和巧玲花,伴生有胡枝子、草珠黃耆、蒙古繡線菊、水栒子、毛榛等;草本層中龍芽草分布較多,還有少量異穗薹草、穿龍薯蕷、串鈴草、歧莖蒿、翠菊、狗尾草、地榆、茜草、艾蒿、華北烏頭等。

    XV紅樺+華北落葉松-銀露梅-細(xì)稈薹草+青綠薹草群落(Comm.Betulaalbosinensis+Larixprincipis-rupprechtii-Potentillaglabra-Carexcapillaris+Carexbreviculmis)。該群落主要分布于海拔2186—2329m的西臺,北坡、西坡,坡度50—55°,上坡位,棕壤,枯枝落葉層厚度2—5cm。群落總蓋度為90%—92%,喬木層蓋度為40%—70%,灌木層蓋度為12%—15%,草本層蓋度為50%—70%。喬木層以紅樺、華北落葉松為主要優(yōu)勢種,伴生有少量青扦、皂柳等;灌木層以銀露梅為主要優(yōu)勢種,伴生有少量土莊繡線菊;草本層以細(xì)稈薹草和青綠薹草為主要優(yōu)勢種,伴生有野青茅、西伯利亞早熟禾(Poasibirica)、密花巖風(fēng)、火絨草(Leontopodiumleontopodioides)、鈴鈴香青(Anaphalishancockii)、兩色烏頭、叉歧繁縷(Stellariadichotoma)、華北藍(lán)盆花(Scabiosatschiliensis)、中國馬先蒿(Pedicularischinensis)、蒲公英、茖蔥(Alliumvictorialis)、毛蕊老鸛草、紅直獐牙菜(Swertiaerythrosticta)、返顧馬先蒿(Pedicularisresupinata)、花錨(Haleniasibirica)、羽節(jié)蕨(Gymnocarpiumjessoense)、珠芽蓼、黑柴胡、皺葉酸模(Rumexcrispus)等。

    XVI紅樺+華北落葉松-密齒柳+青扦幼苗-密花巖風(fēng)群落(Comm.Betulaalbosinensis+Larixprincipis-rupprechtii-Salixcharacta+Piceawilsonii-Libanotiscondensata)。該群落主要分布于海拔2300—2437m的東臺,東北坡,坡度32—40°,上坡位、山頂,棕壤,凋落層略厚,無干擾。群落總蓋度為 50—70%,喬木層蓋度為 55—65%,灌木層蓋度為25%,草本層蓋度為43—70%。喬木層以紅樺、華北落葉松為主要優(yōu)勢種;灌木層優(yōu)勢種為密齒柳和青扦幼苗,伴生有東北茶藨子、黑果栒子(Cotoneastermelanocarpus)等;草本層中密花巖風(fēng)分布較多,還有少量華北馬先蒿、秦艽、藍(lán)刺頭、紅梗蒲公英(Taraxacumerythropodium)、狹苞橐吾(Ligulariaintermedia)、細(xì)葉鴉蔥、寡毛變種(Solidagovirgaureavar.dahurica)、鵝觀草、臭草(Melicascabrosa)、花錨、狹葉黃精、火絨草、鈴鈴香青、蛇床、野韭(Alliumramosum)、地榆、叉歧繁縷、華北藍(lán)盆花、中國馬先蒿、石竹、三脈紫菀、老鸛草等。

    XVII華北落葉松—密齒柳+銀露梅—野青茅群落(Comm.Larixprincipis-rupprechtii—Salixcharacta+Potentillaglabra—Deyeuxiaarundinacea)。該群落主要分布于海拔2397—2589m的東臺,西北坡,坡度35—50°,山頂,棕壤,凋落層3—5cm。群落總蓋度為70—95%,喬木層蓋度為12%,灌木層蓋度為20—25%,草本層蓋度為 70%—95%。喬木層以華北落葉松為主要優(yōu)勢種,喬木層物種數(shù)量稀少;灌木層優(yōu)勢種為密齒柳和銀露梅,伴生有東北茶藨子、全緣栒子(Cotoneasterintegerrimus)等;草本層中野青茅分布較多,還有少量鵝觀草、臭草、花錨、輪葉黃精、細(xì)葉鴉蔥、密花巖風(fēng)、紫苞風(fēng)毛菊、火絨草、鈴鈴香青、叉歧繁縷、華北藍(lán)盆花、中國馬先蒿、大頭風(fēng)毛菊(Saussureabaicalensis)、小竊衣(Torilisjaponica)、瞿麥(Dianthussuperbus)、雪白委陵菜(Potentillanivea)、阿爾泰狗娃花(Heteropappusaltaicus)、北方還陽參(Crepiscrocea)等。

    XVIII華北落葉松-迎紅杜鵑+密齒柳-紫苞風(fēng)毛菊群落(Comm.Larixprincipis-rupprechtii-Rhododendronmucronulatum+Salixcharacta-Saussureaiodostegia)。該群落主要分布于海拔2194—2398m的東臺,北坡,坡度35°,上坡位,棕壤。群落總蓋度為82—88%,喬木層蓋度為10%,灌木層蓋度為 6%,草本層蓋度為 83%。喬木層以華北落葉松為主要優(yōu)勢種,伴生有少量臭冷杉;灌木層優(yōu)勢種為密齒柳和迎紅杜鵑,伴生有碩樺苗;草本層中紫苞風(fēng)毛菊分布較多,還有巴天酸模(Rumexpatientia)、絹茸火絨草(Leontopodiumsmithianum)、鈴鈴香青、花錨、點(diǎn)地梅(Androsaceumbellata)、石竹、長花馬先蒿、黑柴胡、鵝觀草、早熟禾、胭脂花(Primulamaximowiczii)、虎尾草(Lysimachiabarystachys)、橐吾、野罌粟、委陵菜(Potentillachinensis)、飛蓬、密花巖風(fēng)、問荊(Equisetumarvense)、中國馬先蒿、草珠黃耆等。

    XIX金露梅-雪白委陵菜+紫羊茅+異穗薹草群落(Comm.Potentillafruticosa-Potentillanivea+Festucarubra+Carexheterostachya)。該群落主要分布于海拔 2590—2688m的東臺,西坡,坡度30 —45°,山頂,亞高山灌叢草甸土,凋落層厚度3—5cm。群落總蓋度為 93%—96%,灌木層蓋度為15%,草本層蓋度為 95%。灌木層優(yōu)勢種為金露梅;草本層中雪白委陵菜、紫羊茅、異穗薹草(Carexheterostachya)分布較多,還有少量野青茅、細(xì)稈薹草(Carexcapillaris)、叉歧繁縷、密花巖風(fēng)、火絨草、鈴鈴香青、北烏頭、蒲公英、金蓮花、紫苞風(fēng)毛菊、毛蕊老鸛草、藍(lán)花棘豆、小叢紅景天等。

    XX金露梅+銀露梅-嵩草群落(Comm.Potentillafruticosa+Potentillaglabra-Kobresiamyosuroides)。該群落主要分布于海拔2307—2492m的南臺,西南坡、南坡,坡度33—45°,山頂,亞高山灌叢草甸土,凋落層厚度2.7—8cm。群落總蓋度為92%—95%,灌木層蓋度為10%—55%,草本層蓋度為70%—95%。灌木層優(yōu)勢種為金露梅和銀露梅,伴生有鬼箭錦雞兒(Caraganajubata);草本層中嵩草分布較多,還有細(xì)稈薹草、尖嘴薹草(Carexleiorhyncha)、紫苞風(fēng)毛菊、秦艽、鈴鈴香青、卷耳(Cerastiumarvense)、雪白委陵菜、藍(lán)花棘豆、肋柱花(Lomatogoniumcarinthiacum)、華北藍(lán)盆花、胭脂花、翠雀(Delphiniumgrandiflorum)、卵穗荸薺(Heleocharissoloniensis)、西伯利亞三毛草、紫斑風(fēng)鈴草(Campanulapuncatata)、皺葉酸模、黑柴胡等。

    2.2 CCA排序

    典范對應(yīng)分析(Canonical correspondence analysis,CCA)采用包含物種數(shù)據(jù)和環(huán)境數(shù)據(jù)的兩個數(shù)據(jù)矩陣,采用樣方-物種重要值矩陣和樣方-環(huán)境因子矩陣對小五臺山國家級自然保護(hù)區(qū)森林群落進(jìn)行CCA排序,結(jié)果表明,物種CCA排序軸第一軸和第二軸與環(huán)境因子極顯著相關(guān),特征值為0.657和0.403,物種-環(huán)境相關(guān)性為97.30%和90.10%,蒙特卡羅擬合檢驗表明第一軸和前四軸所代表的環(huán)境變量均與物種變量之間呈顯著相關(guān)關(guān)系(第一軸F=5.084,P=0.001;前四軸:F=2.145,P=0.001),前四軸物種-環(huán)境關(guān)系方差累計貢獻(xiàn)率達(dá)到60.2%,采用前兩軸繪制樣方-環(huán)境因子的二維空間排序圖(圖2)。

    圖2 小五臺山國家級自然保護(hù)區(qū)森林群落148個樣方的CCA二維排序圖Fig.2 Two-dimensional CCA ordination diagram of 148 quadrats of forest community in the Xiaowutai Monutain National Nature Reserve Alti: 海拔, Aspe: 坡向,Slop: 坡度,SlPo: 坡位,LiTh: 凋落層厚度,SoTh: 土壤厚度,SoMo: 土壤濕度,SoTe: 土壤溫度,SoEl: 土壤導(dǎo)電率, SopH: 土壤pH,Dist:干擾情況

    Table 2 Intra-set correlations between environmental variables and CCA ordination axes for forest communities in the Xiaowutai Mountain National Nature Reserve

    環(huán)境因子EnvironmentalfactorsCCA排序軸CCAAxes第一軸Axis1第二軸Axis2第三軸Axis3第四軸Axis4海拔Altitude0.9519***-0.04140.0118-0.0255坡向Aspect-0.1655*-0.2344**-0.0532-0.5545***坡度Slope0.2619**0.2180**0.1827*-0.2123**坡位Slopeposition0.7623***-0.1906*0.0555-0.1894*凋落層厚度Litterlayerthickness0.4412***0.6555***-0.0476-0.2168**土壤厚度Soilthickness0.4783***0.5124***0.0229-0.3441***土壤濕度Soilmoisture0.6486***-0.1856*0.1771*0.0986土壤溫度Soiltemperature-0.4695***-0.4889***0.4311***-0.1071土壤導(dǎo)電率Soilelectricalconductivity-0.6680***-0.0347-0.4529***-0.2096**土壤pHSoilpH-0.3017***-0.0093-0.2004*0.4965***干擾情況Disturbance-0.5448***-0.2641**0.2664***0.0623

    *P<0.05,**P<0.01,***P<0.001

    綜合樣方的CCA二維排序圖(圖2)和環(huán)境因子與排序軸的相關(guān)性大小(表2)可以看出,CCA排序第一軸主要反映了以下多個環(huán)境因子的綜合變化趨勢:與海拔、坡位、土壤濕度、土壤厚度、凋落層厚度極顯著正相關(guān),相關(guān)系數(shù)分別為0.9519、0.7623、0.6486、0.4783、0.4412;與土壤導(dǎo)電率、干擾情況、土壤溫度、土壤pH極顯著負(fù)相關(guān),相關(guān)系數(shù)分別為-0.6680、-0.5448、-0.4695、-0.3017,與坡向顯著負(fù)相關(guān),相關(guān)系數(shù)為-0.1655,即沿CCA排序軸第一軸從左到右海拔逐漸升高,土壤類型由山地褐土向山地棕壤、亞高山草甸土轉(zhuǎn)變,坡位趨向于由谷底向下坡位、中坡位、上坡位、山頂轉(zhuǎn)變,土壤濕度由低向高,山頂雨量豐富,較為濕潤,土壤層厚度和枯枝落葉層厚度呈逐漸增加趨勢,土壤導(dǎo)電率逐漸下降,人為或動物干擾逐漸減弱,土壤溫度逐漸降低,土壤由堿性逐漸變?yōu)樗嵝?這11個因子中,海拔與第一軸的相關(guān)性最大。CCA排序第二軸與凋落層厚度、土壤厚度極顯著正相關(guān),相關(guān)系數(shù)分別為0.6555、0.5124,與坡度顯著正相關(guān),相關(guān)系數(shù)為0.2180;與土壤溫度極顯著負(fù)相關(guān),相關(guān)系數(shù)為-0.4889,與干擾情況、坡向顯著負(fù)相關(guān),相關(guān)系數(shù)分別為-0.2641、-0.2344,與坡位、土壤濕度負(fù)相關(guān),相關(guān)系數(shù)為-0.1906、-0.1856,即沿CCA第二軸,從下到上,土壤厚度、凋落層厚度越來越大,坡度由平緩變?yōu)槎盖?土壤溫度逐漸降低,坡向由陽坡轉(zhuǎn)為陰坡,這8個環(huán)境因子中,凋落層厚度與第二軸相關(guān)性最大。綜合前兩軸,海拔梯度、坡位、凋落層厚度、土壤導(dǎo)電率、土壤溫度、土壤濕度、土壤厚度和干擾情況等對森林群落的分布起較大的作用。

    圖3 小五臺山國家級自然保護(hù)區(qū)喬木層優(yōu)勢種的CCA二維排序圖Fig.3 Two-dimensional CCA ordination diagram of dominant species of tree layer in the Xiaowutai Monutain National Nature Reserve

    TWINPSAN分類所得的20個群落類型在排序圖上呈現(xiàn)有規(guī)律的分布,每種植物都有自己的分布中心和分布區(qū)域,如圖2所示,沿第一軸從左到右,適于低海拔、下坡位、土層厚度、凋落層厚度較薄、較為干燥的的群落XII、群落XIII位于排序圖左最左方;適于高海拔、上坡位、土壤溫度較低、較為濕潤的群落XIX位于排序軸的最右方,土壤呈酸性。沿第二軸從下到上,位于排序圖最下方的群落XX到位于排序圖最上方的群落IV,土壤凋落層厚度、土層厚度趨于逐漸增加,土壤溫度呈下降趨勢,由陽坡轉(zhuǎn)為陰坡,各個群落在排序圖中的分布較好地解釋了物種與生境之間的關(guān)系。

    利用前兩軸繪制喬木層優(yōu)勢種的CCA二維排序圖(圖3)和t值雙序圖(圖4)。群落中喬木層優(yōu)勢種的分布格局所揭示的環(huán)境梯度與群落類型的分布梯度有很大的相似性,各個群落沿海拔梯度依次分布(圖3),華北落葉松、紅樺生長于高海拔、冷涼的生境;白樺、山楊生境條件相似,要求溫度、光照逐漸增高;遼東櫟、色木槭、油松喜溫暖光照充足的氣候;杏生長在溫暖向陽的山坡。綜合圖2和圖3可以看出,以杠柳為優(yōu)勢種的群落XIII、以色木槭、遼東櫟為優(yōu)勢種的群落XIV分布在陽坡、低海拔地區(qū);以油松、杏為優(yōu)勢種的群落XII、以榆樹為優(yōu)勢種的群落XI、以油松為優(yōu)勢種的群落Ⅸ分布在低海拔陰坡和半陰坡;群落VIII、X、VI、I、VII、IV、II、III、V位于中心位置,由白樺、油松、山楊為優(yōu)勢種的群落向白樺純林過渡,后轉(zhuǎn)向紅樺、白樺混交林,轉(zhuǎn)向紅樺、華北落葉松混交林,油松、山楊、白樺主要分布在中海拔地區(qū),紅樺、華北落葉松主要分布在高海拔地區(qū);群落XV和群落XVI為喬木層以紅樺、華北落葉松為優(yōu)勢種,群落XVII和群落XVIII為喬木層以華北落葉松為優(yōu)勢種,隨著海拔的升高,由紅樺林向華北落葉松林過渡。群落XIX和群落XX為亞高山灌叢,樣方內(nèi)無喬木,灌木層多為密齒柳、金露梅、銀露梅等。

    利用植物種與特定環(huán)境因子之間回歸系數(shù)t檢驗結(jié)果來分析單個環(huán)境因子是否顯著影響各個優(yōu)勢種的分布,CCA排序得到的前兩軸繪制t值雙序圖,從而解釋各個環(huán)境因子對群落優(yōu)勢種的影響。由樣方的CCA二維排序圖分析得出海拔梯度、坡位、凋落層厚度、土壤導(dǎo)電率、土壤溫度、土壤濕度、土壤厚度和干擾情況等環(huán)境因子對植物種的影響較大,因此,本文繪制喬木層優(yōu)勢種與相關(guān)環(huán)境因子的t值雙序圖。圖4(A)是海拔與優(yōu)勢種回歸系數(shù)的t值雙序圖,可以看出華北落葉松、紅樺分布與海拔呈顯著正相關(guān),油松、遼東櫟、山杏分布與海拔顯著負(fù)相關(guān),表明隨著海拔的升高華北落葉松、紅樺優(yōu)勢度增大,它們趨向于生長在高海拔地區(qū),油松、遼東櫟、山杏優(yōu)勢度逐漸減小,它們趨向于生長在低海拔地區(qū),色木槭、山楊、白樺與海拔的相關(guān)性不大。從圖4(B)中可以看出,喬木層各個優(yōu)勢種與坡位的相關(guān)性不大,坡位不能直接解釋各個喬木樹種的分布。從凋落層厚度與喬木層優(yōu)勢種t值雙序圖上可以看出(圖4(C)),紅樺與凋落層厚度極顯著正相關(guān),其它優(yōu)勢種的箭頭均不能位于響應(yīng)的圓圈內(nèi),表明其它樹種與凋落層厚度相關(guān)性較弱。從圖4(D)可看出,各個樹種與土壤導(dǎo)電率相關(guān)性較弱。從圖4(E)可以看出,紅樺與土壤溫度呈顯著負(fù)相關(guān),其它樹種與土壤溫度相關(guān)性不大。從圖4(F)可看出,各個樹種與土壤濕度相關(guān)性較弱,并無明確指示。從圖4(G)可看出,各個樹種與土壤厚度相關(guān)性不大。干擾程度影響華北落葉松和白樺的分布,其它樹種與干擾程度無關(guān)(圖4(H))。綜上所述,海拔、凋落層厚度、土壤溫度、干擾程度對喬木層優(yōu)勢種的影響較大。

    3 討論

    本文以小五臺山國家級自然保護(hù)區(qū)內(nèi)森林群落為研究對象,選用雙向指示種分析(TWINSPAN)和典范對應(yīng)分析(CCA)等數(shù)量生態(tài)學(xué)方法對該區(qū)山地植物進(jìn)行分類與排序研究。采用樣帶和樣方相結(jié)合的方法沿海拔梯度取樣,在研究區(qū)內(nèi)調(diào)查的148個樣方中,有維管植物75科,222屬,392種,其中蕨類植物6科,8屬,10種;種子植物69科,214屬,382種。在種子植物中,有裸子植物2科,5屬,6種;被子植物67科,209屬,376種。本文所劃分的20個群落類型,群落特征明顯,涉及有油松林、白樺林、紅樺林、華北落葉松林、榆樹林、遼東櫟林、杠柳灌叢、金露梅灌叢、金露梅和銀露梅灌叢等群落類型。其中,將油松林分為IX、XII兩類,多分布在西臺,群落特征明顯,長勢良好,在低海拔處略有人為干擾;將白樺林分為I、V、VI、VIII、X等類群,I、VI為白樺純林,V為白樺、紅樺混交林,VIII為白樺、油松混交林、X為山楊、白樺混交林;將紅樺林分為II、III、IV、VII、XV、XVI等類群,II、III、IV為紅樺純林,VII為紅樺、白樺混交林,XV、XVI為紅樺華北落葉松混交林;華北落葉松林有群落XVII、XVIII兩類,喬木層以華北落葉松為主,灌木層組成不同;群落XI為榆樹林,生長在東臺,多處西北坡,坡度較為平緩;群落XIV以遼東櫟和色木槭為共建種,多生長在的南臺西坡、西南坡、南坡等地;群落XIII為杠柳灌叢,多生長在南臺低海拔處,土壤堿性;群落XIX、XX為亞高山灌叢,以銀露梅、金露梅為主要優(yōu)勢種。前人研究中調(diào)查到本文的部分群落類型,但并不全面,且沒有明確分類[18- 19]。本文劃分出的這些群落類型基本上代表了小五臺自然保護(hù)區(qū)森林類型的全貌,各個群落類型在物種組成以及群落外貌特征等方面均差異明顯,群落分類結(jié)果客觀。

    圖4 小五臺山國家級自然保護(hù)區(qū)喬木層優(yōu)勢種的t值雙序圖Fig.4 Two-dimensional t-value diagram of dominant species of tree layer in the Xiaowutai Monutain National Nature Reserve圓圈代表重點(diǎn)分析的環(huán)境因子,箭頭代表物種與環(huán)境因子的相關(guān)性,箭頭落在圓圈內(nèi),表示該物種與環(huán)境因子顯著相關(guān),落在實(shí)線圈內(nèi),代表顯著正相關(guān),落在虛線圈內(nèi),代表顯著負(fù)相關(guān)[17]

    各個群落類型在CCA二維排序圖上呈現(xiàn)有規(guī)律的分布,每種植物都有自己的分布中心和分布區(qū)域,二維排序圖表明海拔梯度、坡位、土壤厚度、凋落層厚度、土壤溫度、土壤濕度、土壤厚度、干擾程度等對森林群落的分布起較大的作用,影響森林群落的分布格局,形成不同的植被類型,體現(xiàn)了地帶性植被的垂直分布格局。熱量、水分、養(yǎng)分、光照等環(huán)境因子決定植物種的分布格局,其內(nèi)部關(guān)系極為復(fù)雜多變,CCA第一排序軸主要與海拔高度有關(guān),海拔所反映的是水分和溫度條件的綜合梯度[20],這也與其他地區(qū)相關(guān)的研究結(jié)論比較一致[21]。CCA第二排序軸主要與凋落層厚度有關(guān),凋落層越厚,土壤養(yǎng)分越為豐富[22]。雖然CCA排序能夠較好地展示出群落類型、物種分布與環(huán)境之間的關(guān)系,但它并不能有效解決環(huán)境因子間的共線性等問題[23],在今后的研究中,需進(jìn)一步挖掘更多的環(huán)境信息,從而給出更深入全面的環(huán)境解釋。喬木層優(yōu)勢種的CCA二維排序圖所揭示的環(huán)境梯度與群落類型的分布梯度有很大的相似性,在喬木層優(yōu)勢種CCA二維排序圖中(圖3),華北落葉松、紅樺趨向于生長于高海拔、冷涼的生境,遼東櫟、色木槭、油松喜溫暖光照充足的氣候,山杏多生長在溫暖向陽的山坡,白樺和山楊生境條件相似,與前人研究結(jié)果一致[24-27]。t值雙序圖闡明了CCA排序得出的較為重要的8個環(huán)境因子對森林群落喬木層優(yōu)勢種的影響,隨著海拔的升高華北落葉松、紅樺優(yōu)勢度增大,而油松、遼東櫟、杏優(yōu)勢度隨海拔的升高而減?。坏蚵鋵虞^厚、土壤溫度較低的地區(qū),紅樺的優(yōu)勢度較高;華北落葉松和白樺的分布與干擾程度有一定的相關(guān)性。因此,海拔、凋落層厚度、土壤溫度、干擾程度等環(huán)境因子對森林群落喬木層優(yōu)勢種的影響較大。

    綜上所述,小五臺山森林群落TWINSPAN 分類揭示了該區(qū)群落的間斷性,CCA排序揭示了群落的連續(xù)性和生境異質(zhì)性,研究結(jié)果較好地展現(xiàn)了該區(qū)森林群落與環(huán)境因子的關(guān)系,有助于深入理解保護(hù)區(qū)森林群落的空間分布格局及各物種的生態(tài)適應(yīng)特性,為今后保護(hù)區(qū)管理和保護(hù)提供相應(yīng)的理論依據(jù)。

    [1] Dai X, Page B, Duffy K J. Indicator value analysis as a group prediction technique in community classification. South African Journal of Botany, 2006, 72(4):589- 596.

    [2] Borcard D, Gillet F, Legendre P. Numerical Ecology with R. New York: Springer-Verlag, 2011:58- 77.

    [3] 張金屯. 數(shù)量生態(tài)學(xué)(第二版). 北京: 科學(xué)出版社, 2011:123- 127.

    [4] Virtanen R, Luoto M, R?m? T, Mikkola K, Hjort J, Grytnes J A, Birks H J B. Recent vegetation changes at the high-latitude tree line ecotone are controlled by geomorphological disturbance, productivity and diversity. Global Ecology and Biogeography, 2010,19(6): 810- 821.

    [5] Legendre P, Legendre L. Numerical Ecology. 3rd ed. Amsterdam: Elsevier, 2012: 99- 103.

    [6] Vittoz P, Bayfield N, Brooker R, Elston D A, Duff E I, Theurillat J P, Guisan A. Reproducibility of species lists, visual cover estimates and frequency methods for recording high-mountain vegetation. Journal of Vegetation Science, 2010, 21(6):1035- 1047.

    [7] 劉方正,崔國發(fā). 國內(nèi)外保護(hù)區(qū)管理有效性評價方法比較.世界林業(yè)研究, 2013, 26(6):33- 38.

    [8] 魯少波. 河北小五臺山森林生態(tài)系統(tǒng)主要因子間量化關(guān)系研究[D].北京:北京林業(yè)大學(xué), 2009.

    [9] 王志峰, 羅輔燕, 唐婷, 李霄峰, 楊鵬, 李俊鈺. 河北小五臺山海拔梯度上青楊的種群結(jié)構(gòu)和空間分布.西華師范大學(xué)學(xué)報:自然科學(xué)版, 2011, 32(1): 1- 6.

    [10] 劉斌, 陳波, 楊新兵, 魯紹偉, 李少寧. 河北小五臺山天然白樺林草本多樣性影響因素研究. 生態(tài)科學(xué), 2015, 34(2):87- 93.

    [11] 趙建成, 郭書彬, 李盼威. 小五臺山植物志. 北京:科學(xué)出版社, 2011: 2- 4.

    [12] 方精云, 王襄平, 沈澤昊, 唐志堯, 賀金生, 于丹, 江源, 王志恒, 鄭成洋, 朱江玲, 郭兆迪. 植物群落清查的主要內(nèi)容、方法和技術(shù)規(guī)范. 生物多樣性, 2009,17(6):533- 548.

    [13] Zhang J T, Xu B, Li M. Vegetation patterns and species diversity along elevational and disturbance gradients in the Baihua Mountain reserve, Beijing, China. Mountain Research and Development, 2013, 33(2):170- 178.

    [14] Hill M O,milauer P. TWINSPAN for Windows Version 2.3. South Bohemia: Huntingdon & Ceske Budejovice,Centre for Ecology and Hydrology & University of South Bohemia, 2005:45- 59.

    [15] 李林峰. 蘆芽山自然保護(hù)區(qū)青杄林的數(shù)量生態(tài)學(xué)研究[D]. 北京:北京師范大學(xué), 2014.

    [16] Ter Braak C J F,milauer P. CANOCO Reference Manual and CanoDraw for Windows User′s Guide: Software for Canonical Community Ordination (version 4.5). New York: Microcomputer Power Press, 2002:66- 81.

    [18] 尤海舟,畢君,蔡蕾,王超.小五臺山不同林分類型林地水文效應(yīng)及涵養(yǎng)水源生態(tài)功能價值估算. 中國農(nóng)學(xué)通報, 2013, 29(13):60- 65.

    [19] 李白.小五臺山主要喬木群落類型調(diào)查研究.寧夏農(nóng)林科技, 2012, 53(10):87- 88.

    [20] Stella J C, Rodríguez-González P M, Dufour S, Bendix J. 2013. Riparian vegetation research in Mediterranean-climate regions: common patterns, ecological processes, and considerations for management. Hydrobiologia, 719(1):291- 315.

    [21] 任學(xué)敏.太白山主要植物群落數(shù)量分類及其物種組成和豐富度的環(huán)境解釋[D].楊凌:西北農(nóng)林科技大學(xué), 2012.

    [22] Moriyama A, Yonemura S, Kawashima S, Du M Y, Tang Y H. Environmental indicators for estimating the potential soil respiration rate in alpine zone. Ecological Indicators, 2013, 32: 245- 252.

    [23] Davidson T A, Sayer C D, Perrow M, Bramm M Jeppesen E. The simultaneous inference of zooplanktivorous fish and macrophyte density from sub-fossil cladoceran assemblages: a multivariate regression tree approach. Freshwater Biology, 2010, 55(3):546- 564.

    [24] 李豪. 山西歷山保護(hù)區(qū)青檀天然群落數(shù)量生態(tài)學(xué)研究[D].臨汾:山西師范大學(xué), 2013.

    [25] 蘇日古嘎, 張金屯, 張斌, 程佳佳, 田世廣, 張欽弟, 劉素軍. 松山自然保護(hù)區(qū)森林群落的數(shù)量分類和排序. 生態(tài)學(xué)報, 2010, 30(10):2621- 2629.

    [26] Zhang J T, Shao D. Attributes of forest diversity in the yunmeng mountain national forest park in Beijing, China. Applied Ecology and Environmental Research, 2015, 13(3):769- 782.

    [27] 張先平, 王孟本, 佘波, 肖揚(yáng). 龐泉溝國家自然保護(hù)區(qū)森林群落的數(shù)量分類和排序. 生態(tài)學(xué)報, 2006, 26(3):754- 761.

    Relationship between forest communities and the environment in the Xiaowutai Mountain National Nature Reserve, Hebei

    BAI Xiaohang, ZHANG Jintun*, CAO Ke, WANG Yunquan, Sehrish Sadia, CAO Ge

    CollegeofLifeSciences,BeijingNormalUniversity,Beijing100875,China

    The relationship of plants and the environment is a complex evolution process. We discuss species composition of forest communities, ecological characteristics of the population, and the relationship between different plant communities and the environment with quantitative ecology method to protect the stability of forest communities and biodiversity in the area. The forest communities in the Xiaowutai Mountain National Nature Reserve in Hebei province were investigated with a two-way indicator species analysis (TWINSPAN) and canonical correspondence analysis (CCA). All quadrats (n=148) were classified into 20 groups representing 20 plant communities. The results of the CCA ordination clearly reflected the pattern of community distribution, and it was obvious in the gradients of the first and second axis, which mainly represented altitude, slope position, litter layer thickness, soil electrical conductivity, soil temperature, soil moisture, soil thickness, and disturbance. All environmental factors influenced the forest community distribution pattern and the formation of different vegetation types. The CCA ordination of dominant species in tree layer revealed a high degree of similarity between environmental gradients and the distributions of community types. A two-dimensional t-value diagram clarified that altitude, litter layer thickness, soil temperature, and disturbance played an important role in dominant species in the tree layer of forest communities. In this study, we used TWINSPAN classification and CCA ordination to explain the relationships between forest communities and environmental factors. The purpose was to provide a theoretical basis for scientific management of the Xiaowutai Mountain area and to analyze the results that can be used as a reference for research and protection of forest ecosystems in similar areas.

    Xiaowutai Mountain National Nature Reserve; forest community; TWINSPAN; CCA; environmental gradient

    國家自然科學(xué)基金資助項目(31170494)

    2016- 04- 14; 網(wǎng)絡(luò)出版日期:2017- 02- 22

    10.5846/stxb201604140680

    *通訊作者Corresponding author.E-mail: zhangjt@bnu.edu.cn

    白曉航,張金屯,曹科,王云泉,Sehrish Sadia,曹格.河北小五臺山國家級自然保護(hù)區(qū)森林群落與環(huán)境的關(guān)系.生態(tài)學(xué)報,2017,37(11):3683- 3696.

    Bai X H, Zhang J T, Cao K, Wang Y Q, Sehrish Sadia, Cao G.Relationship between forest communities and the environment in the Xiaowutai Mountain National Nature Reserve, Hebei.Acta Ecologica Sinica,2017,37(11):3683- 3696.

    猜你喜歡
    草本喬木灌木
    毛鋪草本酒
    毛鋪草本酒
    ◆ 園林苗木
    彩墨繪草本
    漯河市常見灌木類苗木的整形與修剪要點(diǎn)
    與世隔絕的人世
    詩潮(2017年2期)2017-03-16 10:42:30
    杜鵑等 5 種灌木對PM 2.5 的凈化作用初探
    周碧華的草本人生
    小說月刊(2014年1期)2014-04-23 08:59:59
    沈陽市政材料價格補(bǔ)充信息
    淺談園林喬木的移植
    河南科技(2014年16期)2014-02-27 14:13:33
    99久久精品热视频| 国产精品爽爽va在线观看网站| 色综合亚洲欧美另类图片| 久久久久久久久大av| 欧美bdsm另类| 国产欧美日韩精品亚洲av| 哪里可以看免费的av片| 免费看美女性在线毛片视频| 亚洲内射少妇av| 人人妻人人看人人澡| 国产激情偷乱视频一区二区| 精品久久久噜噜| 日韩一区二区视频免费看| 免费人成在线观看视频色| 亚洲第一电影网av| 联通29元200g的流量卡| 可以在线观看毛片的网站| 免费在线观看日本一区| 亚洲最大成人中文| 免费观看在线日韩| 搡老岳熟女国产| 日韩欧美三级三区| 天堂动漫精品| a级毛片a级免费在线| 午夜福利在线观看吧| 欧美一区二区亚洲| 99久久无色码亚洲精品果冻| 午夜激情福利司机影院| 亚洲久久久久久中文字幕| 十八禁国产超污无遮挡网站| 久久久久久久久久成人| 动漫黄色视频在线观看| 亚洲精品亚洲一区二区| 精品久久久久久久久亚洲 | 久久精品国产99精品国产亚洲性色| 久久精品国产亚洲av天美| 日日夜夜操网爽| 免费av观看视频| 国产乱人伦免费视频| 日本免费一区二区三区高清不卡| 日韩欧美精品v在线| 欧美日韩黄片免| 男人舔女人下体高潮全视频| 亚洲av日韩精品久久久久久密| 99久久精品国产国产毛片| 久久精品国产99精品国产亚洲性色| 麻豆一二三区av精品| 内射极品少妇av片p| 欧美+亚洲+日韩+国产| 成人国产麻豆网| 久久人人爽人人爽人人片va| 我要看日韩黄色一级片| 亚洲国产色片| 欧美性猛交黑人性爽| 亚洲精品国产成人久久av| 久久久久九九精品影院| 久久久国产成人精品二区| 国产精品野战在线观看| 亚洲最大成人中文| 亚洲av中文av极速乱 | 校园春色视频在线观看| 久久午夜福利片| 嫩草影院入口| 能在线免费观看的黄片| 久久99热6这里只有精品| 色尼玛亚洲综合影院| 国产色婷婷99| 又黄又爽又免费观看的视频| 亚洲精品一卡2卡三卡4卡5卡| 国产探花极品一区二区| 中文亚洲av片在线观看爽| 99久久久亚洲精品蜜臀av| 99视频精品全部免费 在线| 精品人妻1区二区| 香蕉av资源在线| 成人国产麻豆网| 乱人视频在线观看| 久久久色成人| 久久精品国产亚洲av天美| 国产激情偷乱视频一区二区| 色哟哟哟哟哟哟| 成年版毛片免费区| 日韩欧美 国产精品| 一进一出好大好爽视频| 亚洲熟妇中文字幕五十中出| 日韩av在线大香蕉| 国产单亲对白刺激| 久99久视频精品免费| 51国产日韩欧美| 亚洲成人久久爱视频| 国内久久婷婷六月综合欲色啪| 亚洲中文日韩欧美视频| 熟女电影av网| av天堂在线播放| 亚洲内射少妇av| 午夜亚洲福利在线播放| 久久久久九九精品影院| 无人区码免费观看不卡| 亚洲avbb在线观看| 久久久久久久亚洲中文字幕| 欧美日本视频| 日本黄大片高清| 亚洲精品一区av在线观看| 麻豆久久精品国产亚洲av| av在线观看视频网站免费| 国产午夜精品久久久久久一区二区三区 | 国产黄色小视频在线观看| 夜夜爽天天搞| 午夜爱爱视频在线播放| 中出人妻视频一区二区| 国产高清三级在线| 中文字幕av在线有码专区| 亚洲中文字幕日韩| www日本黄色视频网| 国产成人一区二区在线| 最好的美女福利视频网| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久电影中文字幕| 女同久久另类99精品国产91| 五月玫瑰六月丁香| 亚洲精品456在线播放app | 亚洲精品久久国产高清桃花| av黄色大香蕉| 成人国产一区最新在线观看| 夜夜爽天天搞| 日韩一区二区视频免费看| 中文字幕av成人在线电影| 哪里可以看免费的av片| 桃色一区二区三区在线观看| 国产高清不卡午夜福利| 国产欧美日韩精品亚洲av| 美女 人体艺术 gogo| 久久久久免费精品人妻一区二区| 久久精品国产自在天天线| 性插视频无遮挡在线免费观看| 九九爱精品视频在线观看| 全区人妻精品视频| or卡值多少钱| 精品99又大又爽又粗少妇毛片 | 18禁黄网站禁片午夜丰满| 在线国产一区二区在线| 国产一区二区在线观看日韩| 久久久久九九精品影院| 色5月婷婷丁香| 午夜a级毛片| 亚洲最大成人手机在线| 波多野结衣巨乳人妻| 亚洲乱码一区二区免费版| 欧美zozozo另类| 久久久久精品国产欧美久久久| 欧美区成人在线视频| АⅤ资源中文在线天堂| 欧美3d第一页| 亚洲第一区二区三区不卡| 亚洲最大成人av| 亚洲avbb在线观看| 欧美+亚洲+日韩+国产| 最近中文字幕高清免费大全6 | 久久久午夜欧美精品| 亚洲无线在线观看| 精品99又大又爽又粗少妇毛片 | 国内精品久久久久久久电影| 看十八女毛片水多多多| 亚洲av日韩精品久久久久久密| 亚洲欧美精品综合久久99| 国产精品1区2区在线观看.| 国产国拍精品亚洲av在线观看| 美女 人体艺术 gogo| 99久久九九国产精品国产免费| 免费电影在线观看免费观看| 国产精品伦人一区二区| 一进一出抽搐gif免费好疼| 美女xxoo啪啪120秒动态图| 色av中文字幕| 熟女电影av网| 天美传媒精品一区二区| 久久久精品大字幕| 日韩精品有码人妻一区| 久久久国产成人免费| 国产探花极品一区二区| 国产一区二区亚洲精品在线观看| 久久精品国产鲁丝片午夜精品 | 老司机午夜福利在线观看视频| 九九在线视频观看精品| 国产精品女同一区二区软件 | 国国产精品蜜臀av免费| 日本三级黄在线观看| 久久久久久九九精品二区国产| 免费观看人在逋| 日本在线视频免费播放| 欧美日本亚洲视频在线播放| 久久久色成人| 天天一区二区日本电影三级| a级毛片免费高清观看在线播放| 免费黄网站久久成人精品| 欧美日韩黄片免| 我要搜黄色片| 99热精品在线国产| 国产主播在线观看一区二区| 成年免费大片在线观看| 欧美精品啪啪一区二区三区| 99国产极品粉嫩在线观看| 亚洲精品国产成人久久av| 久久九九热精品免费| 国产精品伦人一区二区| 久久国产精品人妻蜜桃| 日本一本二区三区精品| 日日摸夜夜添夜夜添av毛片 | 一区二区三区四区激情视频 | 日日撸夜夜添| 丰满乱子伦码专区| а√天堂www在线а√下载| 日本欧美国产在线视频| 性欧美人与动物交配| 乱码一卡2卡4卡精品| 欧美人与善性xxx| 熟女电影av网| 女生性感内裤真人,穿戴方法视频| 麻豆国产av国片精品| 精品99又大又爽又粗少妇毛片 | 少妇裸体淫交视频免费看高清| 亚洲av免费高清在线观看| 尤物成人国产欧美一区二区三区| 欧美不卡视频在线免费观看| 国产高清激情床上av| 亚洲av一区综合| 日韩欧美在线二视频| 又粗又爽又猛毛片免费看| 精品久久国产蜜桃| 亚洲在线自拍视频| www.色视频.com| 亚洲最大成人手机在线| 免费av不卡在线播放| 亚洲狠狠婷婷综合久久图片| 99精品在免费线老司机午夜| 日韩欧美三级三区| 别揉我奶头 嗯啊视频| 国产高潮美女av| 欧美黑人巨大hd| 亚洲成av人片在线播放无| 日韩欧美一区二区三区在线观看| 久久久国产成人精品二区| 色哟哟哟哟哟哟| 亚洲国产精品成人综合色| 成人美女网站在线观看视频| 久久精品国产清高在天天线| 亚洲一区高清亚洲精品| 亚洲电影在线观看av| 国产高清有码在线观看视频| 尾随美女入室| 波多野结衣高清作品| 成人二区视频| 免费看a级黄色片| 欧美成人一区二区免费高清观看| 欧美日韩瑟瑟在线播放| www.色视频.com| 中亚洲国语对白在线视频| 免费人成在线观看视频色| 男女做爰动态图高潮gif福利片| 精品无人区乱码1区二区| 熟妇人妻久久中文字幕3abv| 亚洲性夜色夜夜综合| 成人鲁丝片一二三区免费| 两个人的视频大全免费| 久久久久九九精品影院| 在线观看66精品国产| 中文字幕熟女人妻在线| 亚洲av成人精品一区久久| 久久热精品热| 久久亚洲真实| 搞女人的毛片| 欧美一区二区国产精品久久精品| 久久人妻av系列| 国产 一区精品| 亚洲五月天丁香| 特级一级黄色大片| 嫩草影院精品99| 免费观看精品视频网站| 91麻豆精品激情在线观看国产| 成人鲁丝片一二三区免费| 欧美一区二区精品小视频在线| 亚洲精品久久国产高清桃花| .国产精品久久| 国产成人福利小说| 成熟少妇高潮喷水视频| 成人av一区二区三区在线看| 给我免费播放毛片高清在线观看| 中文字幕久久专区| 日本精品一区二区三区蜜桃| 国产色爽女视频免费观看| 成人国产麻豆网| 又爽又黄无遮挡网站| 嫩草影视91久久| 一边摸一边抽搐一进一小说| 欧美绝顶高潮抽搐喷水| 我要搜黄色片| 国产精品99久久久久久久久| 国产高清视频在线观看网站| 精品久久久久久久久久久久久| 婷婷精品国产亚洲av在线| 日本一二三区视频观看| 熟妇人妻久久中文字幕3abv| 在线观看美女被高潮喷水网站| 看片在线看免费视频| 国产精品,欧美在线| 内地一区二区视频在线| 国产男靠女视频免费网站| 国产成人aa在线观看| 不卡一级毛片| 最近在线观看免费完整版| 一个人看的www免费观看视频| 亚洲最大成人av| 美女被艹到高潮喷水动态| 国产成年人精品一区二区| 亚洲人成网站在线播| 中文字幕av成人在线电影| 好男人在线观看高清免费视频| 真人做人爱边吃奶动态| АⅤ资源中文在线天堂| 免费人成在线观看视频色| 亚州av有码| 国产精品久久电影中文字幕| 国产高清激情床上av| 老司机福利观看| 91av网一区二区| 国产精品美女特级片免费视频播放器| 一区二区三区高清视频在线| 一本一本综合久久| 99视频精品全部免费 在线| 亚洲精华国产精华精| 欧美高清成人免费视频www| 国产精品三级大全| 精品久久久久久久人妻蜜臀av| 免费av毛片视频| 久久午夜亚洲精品久久| 国产精品98久久久久久宅男小说| www.色视频.com| 亚洲不卡免费看| 精品一区二区三区人妻视频| 国产欧美日韩一区二区精品| 三级男女做爰猛烈吃奶摸视频| 观看美女的网站| 午夜视频国产福利| 床上黄色一级片| av在线老鸭窝| 成人国产麻豆网| 中出人妻视频一区二区| 国产爱豆传媒在线观看| 99九九线精品视频在线观看视频| 麻豆av噜噜一区二区三区| 91在线精品国自产拍蜜月| 国产精品久久久久久av不卡| 亚洲欧美日韩高清在线视频| 此物有八面人人有两片| 舔av片在线| 99久国产av精品| 国产在视频线在精品| 国产亚洲91精品色在线| 久99久视频精品免费| 日韩精品中文字幕看吧| 变态另类丝袜制服| 女人被狂操c到高潮| 12—13女人毛片做爰片一| 欧美成人a在线观看| 国产大屁股一区二区在线视频| 成人毛片a级毛片在线播放| 国产一区二区亚洲精品在线观看| 欧美丝袜亚洲另类 | 亚洲第一区二区三区不卡| 岛国在线免费视频观看| 99热精品在线国产| 少妇丰满av| 精品欧美国产一区二区三| 小蜜桃在线观看免费完整版高清| 免费av毛片视频| 午夜爱爱视频在线播放| 精品福利观看| 女生性感内裤真人,穿戴方法视频| 欧美成人免费av一区二区三区| 欧美一区二区精品小视频在线| 午夜亚洲福利在线播放| 亚洲三级黄色毛片| 日韩中字成人| 国产人妻一区二区三区在| 免费无遮挡裸体视频| 日韩强制内射视频| 中文字幕免费在线视频6| 中文亚洲av片在线观看爽| 女生性感内裤真人,穿戴方法视频| 别揉我奶头 嗯啊视频| 国产 一区精品| 久久99热6这里只有精品| 久久欧美精品欧美久久欧美| 久久久国产成人免费| 永久网站在线| 97热精品久久久久久| 亚洲av二区三区四区| 久久久久久久精品吃奶| 波多野结衣巨乳人妻| 女生性感内裤真人,穿戴方法视频| 国内精品久久久久精免费| 97超视频在线观看视频| 欧美日韩精品成人综合77777| videossex国产| av专区在线播放| 亚洲国产精品sss在线观看| 久久精品国产亚洲网站| 国产免费男女视频| 国产蜜桃级精品一区二区三区| 两人在一起打扑克的视频| 国产精品一区二区性色av| 欧美黑人巨大hd| 国产伦一二天堂av在线观看| 五月玫瑰六月丁香| 亚洲第一电影网av| 91av网一区二区| 性色avwww在线观看| 天堂网av新在线| 国产伦精品一区二区三区视频9| 51国产日韩欧美| 久久精品91蜜桃| 亚洲精华国产精华精| 久久精品国产亚洲av涩爱 | 免费看光身美女| 变态另类丝袜制服| 搡老熟女国产l中国老女人| 亚洲av第一区精品v没综合| 黄色配什么色好看| 欧美成人a在线观看| 日本黄色视频三级网站网址| 狂野欧美白嫩少妇大欣赏| 国内少妇人妻偷人精品xxx网站| 国产欧美日韩一区二区精品| 日本五十路高清| 亚洲成a人片在线一区二区| 搡老岳熟女国产| 欧美国产日韩亚洲一区| 久久精品影院6| 免费av观看视频| 欧美一区二区亚洲| 波野结衣二区三区在线| 国产淫片久久久久久久久| 一区二区三区高清视频在线| 啦啦啦啦在线视频资源| 日韩欧美一区二区三区在线观看| 欧美最黄视频在线播放免费| 嫩草影院精品99| 欧美激情久久久久久爽电影| 一本久久中文字幕| 久久精品人妻少妇| 亚洲四区av| 网址你懂的国产日韩在线| 久久国内精品自在自线图片| 一区二区三区高清视频在线| 欧美中文日本在线观看视频| 草草在线视频免费看| 校园春色视频在线观看| 亚洲18禁久久av| 韩国av一区二区三区四区| 中文字幕熟女人妻在线| 色综合色国产| 久久久久国内视频| 又粗又爽又猛毛片免费看| 男人狂女人下面高潮的视频| 男人舔女人下体高潮全视频| 真实男女啪啪啪动态图| or卡值多少钱| 99久久九九国产精品国产免费| 桃红色精品国产亚洲av| 久久久午夜欧美精品| 欧美极品一区二区三区四区| 亚洲美女黄片视频| 亚洲av不卡在线观看| 精品乱码久久久久久99久播| 日韩在线高清观看一区二区三区 | 中文字幕高清在线视频| 69av精品久久久久久| 免费观看在线日韩| 久久99热6这里只有精品| 久久这里只有精品中国| 国产av麻豆久久久久久久| 制服丝袜大香蕉在线| 国产黄片美女视频| 午夜爱爱视频在线播放| 精品久久久久久久久亚洲 | 美女xxoo啪啪120秒动态图| 婷婷精品国产亚洲av| 国产淫片久久久久久久久| 亚洲自拍偷在线| av黄色大香蕉| 十八禁网站免费在线| 国产精品亚洲美女久久久| 69av精品久久久久久| 国产一区二区三区视频了| 久99久视频精品免费| 97人妻精品一区二区三区麻豆| 国产综合懂色| 九九热线精品视视频播放| 搞女人的毛片| 午夜福利视频1000在线观看| 99国产极品粉嫩在线观看| 久久精品国产亚洲av天美| 国产一区二区在线av高清观看| 免费看光身美女| 毛片女人毛片| 日韩一本色道免费dvd| 国产中年淑女户外野战色| 国产一区二区三区av在线 | 欧美高清性xxxxhd video| АⅤ资源中文在线天堂| 亚洲欧美日韩东京热| 午夜精品在线福利| 91久久精品电影网| xxxwww97欧美| 免费电影在线观看免费观看| av在线亚洲专区| www.www免费av| 神马国产精品三级电影在线观看| 亚洲av不卡在线观看| 长腿黑丝高跟| 国产探花极品一区二区| 亚洲欧美日韩高清专用| 精品久久久久久久久亚洲 | 男女之事视频高清在线观看| 亚洲电影在线观看av| 欧美3d第一页| 亚洲18禁久久av| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人影院久久av| 亚洲精品一区av在线观看| 欧美极品一区二区三区四区| 精品久久国产蜜桃| 一个人看的www免费观看视频| 国产精品美女特级片免费视频播放器| 精品国产三级普通话版| 少妇的逼水好多| 亚洲黑人精品在线| 久久精品人妻少妇| 国产av一区在线观看免费| 国产精品爽爽va在线观看网站| .国产精品久久| 一级黄片播放器| 99久久精品热视频| 亚洲午夜理论影院| 国产精品久久久久久av不卡| 日韩精品青青久久久久久| 国产精品永久免费网站| 最近视频中文字幕2019在线8| 九九在线视频观看精品| 欧美3d第一页| 长腿黑丝高跟| 看十八女毛片水多多多| 99热6这里只有精品| 亚洲av.av天堂| 成人永久免费在线观看视频| 亚洲精华国产精华精| 国产av在哪里看| 九色国产91popny在线| 欧美日韩乱码在线| 国产女主播在线喷水免费视频网站 | 亚洲欧美精品综合久久99| 天堂√8在线中文| 成年版毛片免费区| 国内精品久久久久久久电影| 丰满人妻一区二区三区视频av| 欧美激情国产日韩精品一区| 亚洲最大成人手机在线| 一a级毛片在线观看| 午夜视频国产福利| 最后的刺客免费高清国语| 亚洲美女黄片视频| 国产探花在线观看一区二区| а√天堂www在线а√下载| 国产精品综合久久久久久久免费| 又爽又黄无遮挡网站| 亚洲内射少妇av| 波野结衣二区三区在线| 麻豆av噜噜一区二区三区| 亚洲av成人av| 婷婷精品国产亚洲av| 成年免费大片在线观看| 免费av毛片视频| 国产免费男女视频| 乱系列少妇在线播放| 免费av毛片视频| 国产免费男女视频| 人妻少妇偷人精品九色| 一进一出好大好爽视频| 久久精品国产亚洲av涩爱 | 日韩欧美精品免费久久| 午夜福利18| 国产精品电影一区二区三区| videossex国产| 亚洲欧美日韩无卡精品| av专区在线播放| 午夜老司机福利剧场| 日本黄大片高清| 日本三级黄在线观看| 一级毛片久久久久久久久女| 九九在线视频观看精品| 女的被弄到高潮叫床怎么办 | 成人高潮视频无遮挡免费网站| bbb黄色大片| 精品国产三级普通话版| 日日摸夜夜添夜夜添av毛片 | 成人三级黄色视频| 亚洲av中文字字幕乱码综合| 男人的好看免费观看在线视频| 深夜a级毛片| 黄色欧美视频在线观看| 亚洲精品影视一区二区三区av| 人人妻,人人澡人人爽秒播| 国产色婷婷99| 97热精品久久久久久| 亚洲精品成人久久久久久|