太白山自然保護(hù)區(qū)典型森林群落數(shù)量分類(lèi)、排序及多樣性格局

鄭天義,王 丹,姬柳婷,康 冰

西北農(nóng)林科技大學(xué)生命科學(xué)學(xué)院,楊凌 712100

群落分類(lèi)是揭示植被間斷性的有效方法,而排序是揭示植被連續(xù)性的方法[1],因而將分類(lèi)與排序方法相結(jié)合使用,有利于了解該研究區(qū)域的群落特征,進(jìn)而更為全面和透徹地了解植被的分布及其成因,揭示植物群落間以及植物與環(huán)境因子的生態(tài)關(guān)系,這對(duì)研究區(qū)域的植物保護(hù)提供一定的理論指導(dǎo)[2- 5]。目前,在群落分類(lèi)研究上通常使用的方法包括聚類(lèi)分析(Cluster analysis)、雙向指示種分析(Two-way indicator species analysis,TWINSPAN)和多元回歸樹(shù)(Multivariate regression trees,MRT)等[6]。其中多元回歸樹(shù)是一種將響應(yīng)變量和解釋變量相結(jié)合的分類(lèi)方法,每次劃分的分割點(diǎn)取自解釋變量,并且保留具有最大預(yù)測(cè)能力的回歸樹(shù)[7]。排序是將某一區(qū)域內(nèi)調(diào)查的樣方或植物物種排列在一定的空間,使得排序軸反映一定的生態(tài)梯度,從而解釋植被分布與環(huán)境之間的關(guān)系[8- 9]。目前,常用的排序方法有兩類(lèi)：一類(lèi)是基于線性模型的主成分分析(Principal components analysis,PCA)及其衍生出來(lái)的冗余分析(Redundancy analysis,RDA),另一類(lèi)是基于非線性模型的對(duì)應(yīng)分析(Correspondence analysis,CA)及其直接梯度分析版本“典范對(duì)應(yīng)分析”(Canonical correspondence analysis, CCA)[10]。CCA雖與CA共享一套基礎(chǔ)算法,但在排序迭代過(guò)程中加入環(huán)境變量的加權(quán)回歸分析,因而可以揭示環(huán)境因子對(duì)物種組成和群落分布格局的影響,從而更加直觀和合理得對(duì)群落、物種的分布做出生態(tài)解釋[11]。物種多樣性是群落結(jié)構(gòu)和功能復(fù)雜性的量度[12- 13],對(duì)物種多樣性的研究可以評(píng)價(jià)群落的組成結(jié)構(gòu)及變化,對(duì)于保護(hù)瀕臨物種有著重要意義[2]。

太白山是秦嶺主峰,為中國(guó)青藏高原以東大陸第一高峰。其自然地理位置特殊,居暖溫帶之南緣、亞熱帶之北界,獨(dú)特的地勢(shì)和氣候孕育了豐富的生物資源,是一個(gè)保存完好的自然綜合體和自然生態(tài)系統(tǒng)[14],因此有利于森林群落和物種多樣性研究的開(kāi)展。對(duì)太白山植被群落的數(shù)量分類(lèi)研究雖有許多,但都是運(yùn)用聚類(lèi)分析和雙向指示種的分類(lèi)方法[15- 16],而運(yùn)用多元回歸樹(shù)方法的研究尚為空白。因此本研究以太白山自然保護(hù)區(qū)的典型森林群落為研究對(duì)象,嘗試運(yùn)用多元回歸樹(shù)的方法對(duì)森林植被進(jìn)行分類(lèi),并且使用CCA排序來(lái)驗(yàn)證分類(lèi)結(jié)果,最后通過(guò)群落的多樣性來(lái)近一步探究太白山自然保護(hù)區(qū)森林群落分布特征。本研究旨在全面系統(tǒng)地探討太白山自然保護(hù)區(qū)的森林群落分類(lèi)、影響因子和群落多樣性格局,可為深入了解研究區(qū)域的森林植被格局和管理保護(hù)該區(qū)森林資源提供理論依據(jù)。

1 材料和方法

1.1 研究區(qū)概況

太白山國(guó)家級(jí)自然保護(hù)區(qū)(33°49′—34°5′ N,107°22′—107°51′ E)位于秦嶺山脈中段,地處太白縣、眉縣和周至縣三縣交界處,保護(hù)區(qū)東西長(zhǎng)45 km,南北寬34.5 km,總面積56325 hm2,海拔3767 m。由于受山地條件及大氣環(huán)流的影響,保護(hù)區(qū)內(nèi)南、北氣候差異明顯,具典型的亞高山氣候特點(diǎn),并形成明晰的垂直變化和氣候帶。區(qū)內(nèi)動(dòng)植物資源豐富,植被垂直分帶明顯。優(yōu)勢(shì)樹(shù)種有巴山冷杉(Abiesfargesii)、太白紅杉(Larixchinensis)、紅樺(Betulaalbosinensis)、櫟類(lèi)等,國(guó)家一級(jí)重點(diǎn)保護(hù)植物有獨(dú)葉草(Kingdoniauniflora)、紅豆杉(Taxuschinensis)等。太白山南北坡的土壤有明顯的差異,其中南坡地帶性土壤為黃棕壤,北坡則為褐土。

1.2 研究方法[17]

2019年7—8月在太白山國(guó)家級(jí)自然保護(hù)區(qū)內(nèi)進(jìn)行樣地調(diào)查,調(diào)查樣點(diǎn)選擇相對(duì)典型的森林群落。從1190 m到3391 m沿海拔設(shè)置樣地,海拔每升高100 m調(diào)查1—2個(gè)25 m×25 m的喬木樣方、每個(gè)喬木樣方設(shè)置5個(gè)2 m×2 m的灌木樣方和5個(gè)1 m×1 m的草本樣方,共設(shè)置30個(gè)喬木樣方。調(diào)查記錄喬木樣方的物種名稱(chēng)、數(shù)量、高度、胸徑、冠幅,灌木、草本樣方的物種名稱(chēng)、高度、數(shù)量、蓋度。利用GPS測(cè)定地理坐標(biāo)和海拔,羅盤(pán)測(cè)定坡向、坡度,同時(shí)調(diào)查地貌類(lèi)型、土壤類(lèi)型、枯落物厚度、干擾情況等因子。

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

1.3.1重要值

分別計(jì)算每個(gè)樣地喬木、灌木和草本的重要值,其公式如下[1]：

喬木層重要值 =(相對(duì)密度+相對(duì)高度+相對(duì)顯著度)/ 3

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

1.3.2多樣性指數(shù)

計(jì)算不同群落各層次的多樣性指數(shù),分別為物種豐富度(S)、Shannon-Wiener多樣性指數(shù)(H′)、Pielou均勻度指數(shù)(J)和群落總體多樣性指數(shù)(D)[18- 19]。

S=物種數(shù);H′=-∑PilnPi;J=H′/lnS;D=W1D1+W2D2+W3D3;

式中,Pi為物種i的相對(duì)重要值;D1、D2、D3分別為各層的多樣性指數(shù);參數(shù)W1、W2、W3分別為0.5、0.25、0.25。

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

坡向以數(shù)字等級(jí)表示,具體方法如下,以正北為0°,將坡向按照45°的夾角,順時(shí)針?lè)较蛐D(zhuǎn)分為8個(gè)坡向等級(jí),以數(shù)字表示各等級(jí)：北坡用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ù)字越大,表示坡向越向陽(yáng),光照條件越好[20]。干擾包括人類(lèi)活動(dòng)和動(dòng)物活動(dòng),干擾情況分為4個(gè)等級(jí)：無(wú)干擾為1級(jí),輕微干擾為2級(jí),中度干擾為3級(jí),強(qiáng)度干擾為4級(jí)。海拔、坡度和枯落物厚度使用實(shí)際記錄值,共計(jì)5種環(huán)境因子[8]。

1.4 數(shù)據(jù)分析

MRT以環(huán)境因子梯度作為劃分節(jié)點(diǎn),將樣地劃分為盡可能同質(zhì)的類(lèi)別,同時(shí)采用交叉驗(yàn)證來(lái)確定分類(lèi)結(jié)果,因而能夠很好地反映異質(zhì)環(huán)境下群落的結(jié)構(gòu)特征[7,21]。MRT采用R語(yǔ)言mvpart程序包進(jìn)行運(yùn)算,本研究的30個(gè)樣地共有311種植物,其中有45種喬木或灌木由于生長(zhǎng)周期不同而分別記入喬木層和灌木層,因而假設(shè)響應(yīng)變量為356個(gè)。

本研究將物種重要值矩陣(356×30)與環(huán)境因子矩陣(5×30)相結(jié)合來(lái)進(jìn)行分析,同時(shí)采用蒙特卡羅置換檢驗(yàn)對(duì)CCA結(jié)果進(jìn)行顯著性分析,包括CCA分析的有效性和環(huán)境因子影響的顯著度[22]。CCA排序采用R軟件中的Vegan軟件包進(jìn)行分析。

2 結(jié)果與分析

2.1 MRT數(shù)量分類(lèi)

圖1 太白山自然保護(hù)區(qū)森林群落分類(lèi)相對(duì)誤差和交叉驗(yàn)證相對(duì)誤差變化圖Fig.1 Resubstitution and cross-validation relative error for the multivariate regression tree of forest community in Taibai Mountain Nature Reservea：交叉驗(yàn)證相對(duì)誤差變化趨勢(shì),Change trend of cross-validated relative error;b：相對(duì)誤差變化趨勢(shì),Change trend of relative error;丨：標(biāo)準(zhǔn)誤差,Standard error;A：根據(jù)“1-SE”規(guī)則確定的分類(lèi)樹(shù)規(guī)模點(diǎn),The tree size suggested by “1-SE”;Min + 1SE：交叉驗(yàn)證相對(duì)誤差最小值加上一個(gè)標(biāo)準(zhǔn)誤差線,One standard error above the minimum cross-validated relative error;Inf：無(wú)窮大,Infinity

圖2 太白山自然保護(hù)區(qū)森林的多元回歸樹(shù)分類(lèi)樹(shù)狀圖 Fig.2 Dendrogram of multivariate regression trees classification in Taibai Mountain Nature Reserve forest communityError：相對(duì)誤差,Relative error;CV Error：交叉檢驗(yàn)相對(duì)誤差,Relative error of cross validation;SE：標(biāo)準(zhǔn)差,Standard error;Alti：海拔,Altitude;LiTh：枯落物厚度,Litter thickness;Ⅰ—Ⅵ：群落類(lèi)型,Community types;Ⅰ：栓皮櫟+蒙古莢蒾+野棉花群落,Quercus variabilis + Viburnum mongolicum + Anemone vitifolia community;Ⅱ：油松+馬桑+紫草群落,Pinus tabuliformis + Coriaria nepalensis + Lithospermum erythrorhizon community;Ⅲ：華山松+巴山木竹+川西鱗毛蕨群落,Pinus armandii + Bashania fargesii + Dryopteris rosthornii community;Ⅳ：銳齒槲櫟+毛櫻桃+華蟹甲群落,Quercus aliena var. acutiserrata + Cerasus tomentosa + Sinacalia tangutica community;Ⅴ：巴山冷杉+太白杜鵑+橐吾群落,Abies fargesii + Rhododendron purdomii + Ligularia sibirica community;Ⅵ：太白紅杉+冠果忍冬+野青茅群落,Larix chinensis + Lonicera stephanocarpa + Deyeuxia arundinacea community;n代表每個(gè)群落包含的樣地?cái)?shù)

如圖1所示,根據(jù)1-SE(1-標(biāo)準(zhǔn)誤差)規(guī)則：在保證通過(guò)交叉驗(yàn)證誤差(CVRE,通過(guò)交叉驗(yàn)證獲得)盡量小的范圍內(nèi),選取規(guī)模最小的樹(shù),最終把“A”點(diǎn)作為回歸樹(shù)規(guī)模[23- 24]。依據(jù)植物群落分類(lèi)和命名原則,將30個(gè)樣地劃分為6類(lèi)群落。由圖2可知,6類(lèi)群落基本以海拔梯度為節(jié)點(diǎn)劃分,因而群落順序以海拔梯度從低到高排列,群落Ⅱ和群落Ⅲ處于相同梯度,因而以枯落物厚度從低到高依次排列。

群落Ⅰ 栓皮櫟(Quercusvariabilis)+蒙古莢蒾(Viburnummongolicum)+野棉花(Anemonevitifolia)群落。包括樣地：1、2、5共3個(gè)樣地。喬木層優(yōu)勢(shì)種為栓皮櫟,重要值64.83%,伴生種為華山松(Pinusarmandii)、君遷子(Diospyroslotus)等;灌木層優(yōu)勢(shì)種為蒙古莢蒾,重要值15.39%,伴生種為樺葉莢蒾(Viburnumbetulifolium)、黃櫨(Cotinuscoggygria)、木姜子(Litseapungens)等;草本層優(yōu)勢(shì)種為野棉花,重要值13.82%,伴生種為唐松草(Thalictrumaquilegiifoliumvar.sibiricum)、披針薹草(Carexlancifolia)、野青茅(Deyeuxiaarundinacea)等。

群落Ⅱ 油松(Pinustabuliformis)+馬桑(Coriarianepalensis)+紫草(Lithospermumerythrorhizon)群落。包括樣地：10、11、12、13共4個(gè)樣地。喬木層優(yōu)勢(shì)種為油松,重要值50.59%,伴生種為鹽膚木(Rhuschinensis)、漆(Toxicodendronvernicifluum)、刺楸(Kalopanaxseptemlobus)等;灌木層優(yōu)勢(shì)種為馬桑,重要值13.46%,伴生種為栓翅衛(wèi)矛(Euonymusphellomanus)、蒙古莢蒾、木姜子等;草本層優(yōu)勢(shì)種為紫草,重要值7.98%,伴生種為唐松草、荻(Triarrhenasacchariflora)、野青茅等?？萋湮锖穸燃s為2 cm。

群落Ⅲ 華山松+巴山木竹(Bashaniafargesii)+川西鱗毛蕨(Dryopterisrosthornii)群落。包括樣地：16、18、19、20共4個(gè)樣地。喬木層優(yōu)勢(shì)種為華山松,重要值53.9%,伴生種為油松、銳齒槲櫟(Quercusalienavar.acutiserrata)等;灌木層優(yōu)勢(shì)種為巴山木竹,重要值8.95%,伴生種為二翅六道木(Abeliamacrotera)、唐棣(Amelanchiersinica)、粘毛忍冬(Lonicerafargesii)等;草本層優(yōu)勢(shì)種為川西鱗毛蕨,重要值9.71%,伴生種為羽節(jié)蕨(Gymnocarpiumjessoense)、糙蘇(Phlomisumbrosa)、野棉花等?？萋湮锖穸燃s為4 cm。

群落Ⅳ 銳齒槲櫟+毛櫻桃(Cerasustomentosa)+華蟹甲(Sinacaliatangutica)群落。包括3、4、6、14、17、7、8、9共8個(gè)樣地。喬木層優(yōu)勢(shì)種為銳齒槲櫟,重要值17.38%,伴生種為遼東櫟(Quercuswutaishanica)、青楊(Populuscathayana)、華山松等;灌木層優(yōu)勢(shì)種為毛櫻桃,重要值6.71%,伴生種為米面蓊(Buckleyalanceolate)、巴山木竹、野山楂(Crataeguscuneata)等。草本層優(yōu)勢(shì)種為華蟹甲,重要值9.65%,伴生種為翼柄翅果菊(Pterocypselatriangulata)、蟹甲草(Parasenecioforrestii)、升麻(Cimicifugafoetida)等。

群落Ⅴ 巴山冷杉+太白杜鵑(Rhododendronpurdomii)+橐吾(Ligulariasibirica)群落。包括樣地：21、22、29、15、23、24共6個(gè)樣地。喬木層優(yōu)勢(shì)種為巴山冷杉,重要值45.21%,伴生種為紅樺、川柳(Salixhylonoma)等;灌木層優(yōu)勢(shì)種為太白杜鵑,重要值10.26%,伴生種為豪豬刺(Berberisjulianae)、叉子圓柏(Sabinavulgaris)、漸尖茶藨子(Ribestakare)等;草本層優(yōu)勢(shì)種為橐吾,重要值7.39%,伴生種為艾(Artemisiaargyi)、蕨(Pteridiumaquilinumvar.latiusculum)、陜西蛾眉蕨(Lunathyriumgiraldii)等。

群落Ⅵ 太白紅杉+冠果忍冬(Lonicerastephanocarpa)+野青茅群落。包括樣地：25、26、27、28、30共5個(gè)樣地。喬木層優(yōu)勢(shì)種為太白紅杉,重要值p51.35%,伴生種為巴山冷杉、紅樺等;灌木層優(yōu)勢(shì)種為冠果忍冬,重要值45.84%,伴生種為香柏(Juniperuspingiivar.wilsonii)、太白杜鵑、紅脈忍冬(Loniceranervosa)等;草本層優(yōu)勢(shì)種為野青茅,重要值14.91%,伴生種為披針薹草、千里光(Senecioscandens)、獨(dú)葉草等。

2.2 群落分布與環(huán)境因子關(guān)系

通過(guò)對(duì)研究區(qū)域的30個(gè)樣地進(jìn)行典范對(duì)應(yīng)分析,5個(gè)環(huán)境因子的解釋量為24.11%,前4個(gè)排序軸的特征值分別為0.8234、0.5240、0.3899、0.3474,CCA排序結(jié)果有效性的蒙特卡羅置換檢驗(yàn)結(jié)果表明在0.001水平上顯著,因此本次CCA排序在一定程度上能很好地解釋植物群落的分布格局與環(huán)境因子的關(guān)系。

由表2可知,與第一排序軸正相關(guān)的環(huán)境因子依次為：海拔、坡向、干擾情況和枯落物厚度,與第一排序軸負(fù)相關(guān)的環(huán)境因子為坡度;與第二排序軸正相關(guān)的環(huán)境因子依次為：枯落物厚度、干擾情況、坡向和海拔,與第二排序軸負(fù)相關(guān)的環(huán)境因子依次為坡度。從顯著性檢驗(yàn)來(lái)看,第一排序軸基本是海拔主導(dǎo)的,第二排序軸是坡度、枯落物厚度和干擾情況主導(dǎo)的。

表1 太白山自然保護(hù)區(qū)森林群落多元回歸樹(shù)(MRT,multivariate regression trees)分類(lèi)的基本特征

結(jié)合群落分類(lèi)和排序結(jié)果可以分析6個(gè)群落的分布規(guī)律。其中,栓皮櫟+蒙古莢蒾+野棉花群落(群落Ⅰ)分布在海拔低,坡度大的區(qū)域,位于排序圖的左下方;油松+馬桑+紫草群落 (群落Ⅱ)和華山松+巴山木竹+川西鱗毛蕨群落(群落Ⅲ)均分布在海拔低且坡度平緩的區(qū)域,位于排序圖的左方,群落Ⅲ的枯落物厚度大于群落Ⅱ;銳齒槲櫟+毛櫻桃+華蟹甲群落(群落Ⅳ)分布在海拔中等,坡度平緩且干擾強(qiáng)度適中的區(qū)域,位于排序圖的上中部和左下部;巴山冷杉+太白杜鵑+橐吾群落(群落Ⅴ)分布在海拔較高、坡度平緩、干擾較大的區(qū)域,位于排序圖的右方中部;太白紅杉+冠果忍冬+野青茅群落(群落Ⅵ)分布在海拔高,坡度較陡且干擾較大的區(qū)域,位于排序圖的右方中部。

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

圖3 太白山自然保護(hù)區(qū)森林群落30個(gè)樣方的典范對(duì)應(yīng)分析(CCA,canonical correspondence analysis)二維排序圖Fig.3 Two-dimensional CCA ordination diagram of 30 plots of Taibai Mountain Nature Reserve forest communityAlti：海拔Altitude;Aspe：坡向Aspect;LiTh：枯落物厚度Litter thickness;Dist：干擾情況Disturbance;Slop：坡度Slope;1—30代表樣方號(hào);Ⅰ—Ⅵ：群落類(lèi)型,Community types

2.3 群落總體多樣性

結(jié)果由圖4、5、6表明,群落Ⅳ的三種多樣性指數(shù)均為最高,群落Ⅰ的物種豐富度和Pielou均勻度最低,群落Ⅵ的Shannon-Wiener多樣性指數(shù)最低。群落Ⅱ和群落Ⅲ處于同一海拔,但群落Ⅱ的三種多樣性指數(shù)均較高于群落Ⅲ。除了群落Ⅱ和群落Ⅲ處于同一梯度海拔以外,其余森林群落基本以海拔梯度分布,因而可以看出位于中海拔區(qū)域森林群落的三種總體多樣性指數(shù)均較高,而位于高海拔和低海拔區(qū)域森林群落的三種總體多樣性指數(shù)均較低。

圖4 太白山自然保護(hù)區(qū)森林群落的物種豐富度差異 Fig.4 Differences in species richness of forest community Taibai Mountain Nature Reserve

圖5 太白山自然保護(hù)區(qū)森林群落的多樣性指數(shù)差異 Fig.5 Differences in diversity index of forest community Taibai Mountain Nature Reserve

圖6 太白山自然保護(hù)區(qū)森林群落的均勻度指數(shù)差異 Fig.6 Differences in evenness index of forest community Taibai Mountain Nature Reserve

3 討論

本次研究使用多元回歸樹(shù)(MRT)將研究區(qū)域的30個(gè)樣地分為6個(gè)森林群落。從圖2可以看出群落Ⅱ和群落Ⅲ由枯落物厚度差異來(lái)劃分,而其余群落均由海拔差異來(lái)劃分,因而在一定程度上反應(yīng)了生境條件的梯度變化。其中群落Ⅱ和群落Ⅲ、群落Ⅴ和群落Ⅵ的物種組成相似度分別較高。群落Ⅳ中包含著較多的樣地和物種,但由于樣地的生境和物種組成結(jié)構(gòu)相似度高,因而將其歸為一個(gè)群落,而有些樣地雖地理位置距離較遠(yuǎn),但仍被劃為一個(gè)群落,這可能是由于太白山地形復(fù)雜進(jìn)而產(chǎn)生的小氣候?qū)е碌?。群落Ⅰ、Ⅱ位于低海拔區(qū)域,熱量充足,而栓皮櫟喜暖,因此在兩個(gè)群落中均有出現(xiàn)。群落Ⅵ位于高海拔區(qū)域,溫度低,因而存在許多耐寒的植物,如太白紅杉、銀露梅(Potentillaglabra)、秦嶺沙參(Adenophorapetiolata)、金蓮花(Trolliuschinensis)、獨(dú)葉草等。在其他關(guān)于太白山森林群落的分類(lèi)研究中,唐志堯等[15]使用TWINSPAN將太白山植被分為銳齒槲櫟林、遼東櫟林、紅樺林、糙皮樺(Betulautilis)林、巴山冷杉林、太白落葉松(太白紅杉)林和高山灌叢7種類(lèi)型,并且沿海拔從下至上依次分布,而張善紅等[25]則采用了地理信息系統(tǒng)將太白山植被分為灌叢草甸、太白紅杉、巴山冷杉、樺林、櫟林和非林地六種植被類(lèi)型?？梢钥闯?位于較高海拔區(qū)域的植被群落分類(lèi)結(jié)果與本研究一致,而位于中低海拔區(qū)域植被復(fù)雜,不同的樣地設(shè)置和研究方法可能導(dǎo)致分類(lèi)結(jié)果與本研究略有不同。

通過(guò)對(duì)5個(gè)解釋因子進(jìn)行蒙特卡羅置換檢驗(yàn),得出海拔、枯落物厚度、坡度和干擾情況4種因子對(duì)太白山自然保護(hù)區(qū)的植被分布格局有著顯著影響。MRT分類(lèi)和CCA排序的結(jié)果均反映了海拔是影響太白山森林群落分布的最主要的原因,6個(gè)群落基本在CCA排序圖上沿著第一排序軸規(guī)律分布。第一排序軸主要與海拔高度有關(guān),這與其他相關(guān)研究結(jié)論一致[26- 27],海拔高度屬于綜合性的因子,其變化會(huì)影響水分、溫度等環(huán)境因子[28]。第二排序軸主要與坡度、枯落物厚度、干擾情況有關(guān)。坡度會(huì)影響土壤的持水能力,從而限制了植物生長(zhǎng)發(fā)育所需的水分[29],枯落物厚度越大土壤養(yǎng)分越豐富[30]。不同的干擾強(qiáng)度也會(huì)影響植物的分布[8,31]。群落Ⅳ雖位于中部海拔區(qū)域,地勢(shì)平緩且物種多樣性高,但有著中等干擾,這可能與附近居民采草藥和旅游等活動(dòng)有關(guān);群落Ⅴ、Ⅵ雖處于高海拔區(qū)域,但依舊有部分干擾,這可能與當(dāng)?shù)芈糜螛I(yè)開(kāi)發(fā)有關(guān),因此在進(jìn)行旅游業(yè)開(kāi)發(fā)的同時(shí),應(yīng)注重對(duì)植被的保護(hù)和管理,防止過(guò)度開(kāi)發(fā)對(duì)太白山植被的破壞。在一些研究中,坡向?qū)χ脖蝗郝涞姆植加兄@著影響[2,32],而本次研究中坡向?qū)ι秩郝涞姆植嫉挠绊戨m表現(xiàn)為與第一排序軸的相關(guān)系數(shù)較高,但是并不顯著,這可能是研究范圍不一致所導(dǎo)致的。在解釋量方面,本次研究所使用的環(huán)境因子未對(duì)剩余的75.89%做出解釋,盡管在群落研究中環(huán)境解釋量允許較低[33],但在以后的研究中應(yīng)考慮更多的環(huán)境因子,如土壤的理化性質(zhì)、土壤厚度和巖石暴露率等。

通過(guò)對(duì)不同群落的總體多樣性指數(shù)進(jìn)行分析,得出位于中海拔地區(qū)的群落Ⅳ的物種豐富度、Shannon-Wiener多樣性和Pielou均勻度指數(shù)均為最高,這可能與中海拔地區(qū)的溫度、光照、水分都相對(duì)適中,并且坡度較為平緩有關(guān);而位于低海拔的群落Ⅰ、Ⅱ、Ⅲ的總體多樣性偏低,這可能是因?yàn)樵搮^(qū)域雖熱量充足,但缺乏充足的光照和水分導(dǎo)致的。而群落Ⅴ、Ⅵ的總體多樣性隨著海拔的繼續(xù)升高而逐漸減少,這與高海拔區(qū)域溫度較低有關(guān),許多植物因耐寒能力不足而無(wú)法生存。從物種豐富度、Shannon-Wiener多樣性上看,群落總體多樣性在海拔高度層面呈現(xiàn)出單峰模型,這與之前學(xué)者的研究一致[17],而Pielou均勻度指數(shù)在海拔高度層面所呈現(xiàn)的規(guī)律不明顯。本次研究群落總體多樣性是通過(guò)經(jīng)驗(yàn)加權(quán)的方式進(jìn)行計(jì)算,并將加權(quán)和未加權(quán)的多樣性指數(shù)進(jìn)行相關(guān)性分析檢驗(yàn),結(jié)果表明兩種方法所得到的群體多樣性指數(shù)在0.01水平上顯著,因而本次研究所使用的加權(quán)系數(shù)是合理的。張峰等[18]認(rèn)為經(jīng)驗(yàn)權(quán)重系數(shù)是比較合理的,它沒(méi)有影響客觀生態(tài)意義反映,同時(shí)參數(shù)值也突出體現(xiàn)了群落中不同結(jié)構(gòu)、不同層次對(duì)于群落總體多樣性具有不同的貢獻(xiàn)的特點(diǎn)。

太白山植物種類(lèi)豐富,有著明顯的植被垂直分帶,而中低部海拔區(qū)域生境優(yōu)越,孕育著豐富的物種,復(fù)雜的物種分布使得群落的劃分較為困難,因而本次研究嘗試將MRT分類(lèi)和CCA排序相結(jié)合的方法對(duì)其進(jìn)行研究,將分類(lèi)的間斷性和排序的連續(xù)性相結(jié)合,能夠更合理地探究群落組成、揭示群落分布和環(huán)境的關(guān)系。通過(guò)對(duì)比不同群落的總體多樣性差異,進(jìn)一步研究森林群落多樣性在海拔水平的變化格局,這對(duì)揭示太白山生物多樣性環(huán)境梯度變化規(guī)律具有重要意義。本次研究有助于深入了該保護(hù)區(qū)森林群落的空間分布規(guī)律,掌握不同物種的生活習(xí)性,保護(hù)珍稀物種,為今后保護(hù)區(qū)森林的管理、利用和保護(hù)提供一定的理論依據(jù)。

4 結(jié)論

本研究通過(guò)多元回歸樹(shù)(MRT)對(duì)太白山自然保護(hù)區(qū)典型森林群落進(jìn)行數(shù)量分類(lèi),采用典范對(duì)應(yīng)分析(CCA)進(jìn)行排序,并對(duì)分類(lèi)得到的森林群落的總體多樣性進(jìn)行分析。主要結(jié)論如下：

(1)將太白山自然保護(hù)區(qū)典型森林群落分為6類(lèi)：栓皮櫟+蒙古莢蒾+野棉花群落(群落Ⅰ)、油松+馬桑+紫草群落(群落Ⅱ)、華山松+巴山木竹+川西鱗毛蕨群落(群落Ⅲ)、銳齒槲櫟+毛櫻桃+華蟹甲群落(群落Ⅳ)、巴山冷杉+太白杜鵑+橐吾群落(群落Ⅴ)和太白紅杉+冠果忍冬+野青茅群落(群落Ⅵ)。

(2)CCA排序反映了太白山森林群落與環(huán)境因子的關(guān)系,結(jié)果表明：海拔、坡度、枯落物厚度和干擾情況4個(gè)變量對(duì)本研究區(qū)群落分布有顯著影響。

(3)位于中海拔區(qū)域森林群落的3種總體多樣性指數(shù)(物種豐富度、Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù))均較高,而位于高海拔和低海拔區(qū)域森林群落的3種總體多樣性指數(shù)相對(duì)較低。