葉 鵬,湯孟平

(1.浙江農(nóng)林大學(xué) 環(huán)境與資源學(xué)院,杭州 311300；2.浙江省龍泉市林業(yè)局,浙江 麗水 323700;3.浙江農(nóng)林大學(xué) 省部共建亞熱帶森林培育國家重點(diǎn)實(shí)驗(yàn)室,杭州 311300)

常綠闊葉林是我國亞熱帶地區(qū)最復(fù)雜、生產(chǎn)力最高、生物多樣性最豐富的地帶性植被類型之一,對(duì)保護(hù)環(huán)境、維持全球性碳循環(huán)的平衡和人類的持續(xù)發(fā)展等都具有極重要的作用[1]。郁閉度和樹冠重疊度是描述森林系統(tǒng)狀態(tài)與環(huán)境指標(biāo)的重要特征因子,對(duì)森林內(nèi)部環(huán)境、群落結(jié)構(gòu)、林木生長等都有重要的影響；郁閉度還常被作為控制擇伐和撫育采伐強(qiáng)度的指標(biāo),也是區(qū)分未成林造林地、有林地和疏林地的指標(biāo)之一,在森林經(jīng)營中具有極其重要的作用[2]。

目前,對(duì)常綠闊葉林的研究主要關(guān)注物種組成及樹種多樣性[3-6]、空間結(jié)構(gòu)[7-10]、生物量[11-13]等方面；對(duì)郁閉度和樹冠重疊度的研究主要集中于某一樹種,如木荷(SchimasuperbaGardn.et Champ)、馬尾松(PinusmassonianaLamb.)等的郁閉度與樹冠重疊度特征及相互關(guān)系和影響因子[14-16]、冠層覆蓋度和冠層郁閉度的測(cè)量方法和技術(shù)研究[17-19]等方面。但是,對(duì)于常綠闊葉林內(nèi)多樹種的樹冠重疊度和郁閉度分析及其與各生長因子之間的關(guān)系研究還鮮有報(bào)道。

采用GIS與樹冠投影法相結(jié)合的方法,既能保證郁閉度的測(cè)量精度,同時(shí)還能提高工作效率[2]。本研究應(yīng)用GIS的空間分析功能,對(duì)常綠闊葉林樹冠重疊度和郁閉度進(jìn)行計(jì)算與分析,探究常綠闊葉林不同樹種的郁閉度和樹冠重疊度與各生長因子之間的關(guān)系,并確定優(yōu)勢(shì)樹種,分析優(yōu)勢(shì)樹種對(duì)林分郁閉度的貢獻(xiàn)度,旨在為常綠闊葉林的科學(xué)經(jīng)營提供理論依據(jù)。

1 研究區(qū)概況

浙江省天目山國家級(jí)自然保護(hù)區(qū)位于浙江省西北部臨安市境內(nèi)的西天目山,距省會(huì)杭州94km,地理坐標(biāo)為30°18′30″～30°24′55″N,119°23′47″～119°28′27″E。區(qū)內(nèi)年平均氣溫為8.8～14.8℃,≥10℃積溫2 500～5 100℃,年降水量1 390～1 870 mm,相對(duì)濕度76%～81%。自然保護(hù)區(qū)受海洋暖濕氣候影響較深,具有中亞熱帶向北亞熱帶過渡特征,森林植被十分茂盛。研究區(qū)內(nèi)植物區(qū)系古老,成分復(fù)雜,有種子植物1 718種(國家重點(diǎn)保護(hù)植物35種,占浙江省重點(diǎn)保護(hù)植物的64%),有蕨類植物171種,是華東地區(qū)植被保存較完好的地區(qū)之一[20]。

2 研究方法

2.1 標(biāo)準(zhǔn)地設(shè)置

標(biāo)準(zhǔn)地設(shè)置采用相鄰網(wǎng)格調(diào)查法,把100m×100m的固定標(biāo)準(zhǔn)地分割成10m×10m的網(wǎng)格(圖1)。以每個(gè)網(wǎng)格為調(diào)查單元,用全站儀圍測(cè)標(biāo)準(zhǔn)地邊界,并測(cè)定各調(diào)查單元角樁基部坐標(biāo)(x,y,z)。

2.2 每木調(diào)查

在每個(gè)調(diào)查單元內(nèi),對(duì)胸徑(DBH)≥5cm的木本植物進(jìn)行每木調(diào)查,記錄樹木種類,測(cè)定每株樹木的胸徑、樹高、活枝下高、樹冠等因子,并采用激光對(duì)中全站儀(徠卡TCR702Xrange)測(cè)定每株樹木坐標(biāo)(x,y,z)。

2.3 郁閉度和樹冠重疊度計(jì)算方法

分別對(duì)標(biāo)準(zhǔn)地內(nèi)各樹種的樹冠投影面積和各樹種的林冠投影面積進(jìn)行計(jì)算。為了便于區(qū)分,把標(biāo)準(zhǔn)地的林冠投影面積與標(biāo)準(zhǔn)地面積之比稱為林分郁閉度,把某一樹種的林冠投影面積與標(biāo)準(zhǔn)地面積之比稱為樹種郁閉度。

首先,利用GIS的緩沖區(qū)分析功能,以樹冠半徑作為緩沖區(qū)半徑,建立樹冠投影圖層；然后,利用Clip(裁剪)工具生成標(biāo)準(zhǔn)地內(nèi)林冠投影圖層,并計(jì)算標(biāo)準(zhǔn)地內(nèi)林冠投影面積,按(1)式計(jì)算林分郁閉度,按(2)式計(jì)算樹種郁閉度；最后,用擦除法(Erase工具)或用裁剪法(Clip工具)求得標(biāo)準(zhǔn)地內(nèi)樹冠投影總面積,按(3)式計(jì)算樹冠重疊度。

(1)

(2)

(3)

2.4 優(yōu)勢(shì)樹種劃分方法

優(yōu)勢(shì)樹種劃分采用優(yōu)勢(shì)度分析法,具體步驟為:首先,計(jì)算群落喬木層中每個(gè)樹種的相對(duì)胸高斷面積,并將其作為優(yōu)勢(shì)度,各樹種按優(yōu)勢(shì)度從大到小排序；其次,通過(4)式確定優(yōu)勢(shì)樹種數(shù)量。

如果群落只有一個(gè)優(yōu)勢(shì)樹種,則優(yōu)勢(shì)樹種的理想百分比為100%；如果有兩個(gè)優(yōu)勢(shì)樹種,則它們的理想百分比為50%；如果有3個(gè)優(yōu)勢(shì)樹種,則理想百分比為33.3%；依次類推,分別計(jì)算d值,當(dāng)d值最小時(shí)的上位種數(shù)為群落喬木層的優(yōu)勢(shì)樹種數(shù)[8]。

(4)

式中:xi為排序在前的上位種中(T)的相對(duì)斷面積,x為優(yōu)勢(shì)樹種所占的理想百分比,xj為除上位種外的剩余種(U)的相對(duì)胸高斷面積,N為總種數(shù)。

2.5 不同樹種對(duì)林分郁閉度的貢獻(xiàn)度計(jì)算方法

將每一樹種的樹冠投影面積按與其它樹種有無重疊分為兩部分。對(duì)于與其它樹種無重疊的樹冠投影面積,直接納入該樹種對(duì)林分郁閉度的貢獻(xiàn)度計(jì)算之中；對(duì)于與其它樹種相重疊的樹冠投影面積,按重疊面積中不同樹種的平均樹高進(jìn)行調(diào)整,再納入該樹種對(duì)林分郁閉度的貢獻(xiàn)計(jì)算之中。

重疊部分面積中,各樹種納入貢獻(xiàn)度計(jì)算之中的樹冠投影面積根據(jù)(5)式進(jìn)行調(diào)整,調(diào)整后，根據(jù)(6)式計(jì)算不同樹種對(duì)林分郁閉度的貢獻(xiàn)度。

(5)

(6)

2.6 聚集指數(shù)計(jì)算

聚集指數(shù)是檢驗(yàn)種群空間分布格局的常用指數(shù)[21]。計(jì)算公式為:

(7)

式中:ri表示第i株樹到其最近鄰樹的距離；F表示標(biāo)準(zhǔn)地面積；N表示標(biāo)準(zhǔn)地內(nèi)此樹種總株數(shù)。

若R<1,那么這塊標(biāo)準(zhǔn)地的該樹種屬于聚集分布；R=1,屬于隨機(jī)分布；R>1,屬于均勻分布。

3 結(jié)果與分析

3.1 林分及各樹種郁閉度和樹冠重疊度分析

3.1.1林分因子及郁閉度和樹冠重疊度計(jì)算結(jié)果

將標(biāo)準(zhǔn)地內(nèi)所有樹木按樹種進(jìn)行區(qū)分,并對(duì)相關(guān)調(diào)查因子進(jìn)行分類匯總。各樹種調(diào)查因子統(tǒng)計(jì)特征如表1所示。經(jīng)調(diào)查,標(biāo)準(zhǔn)地內(nèi)樹木共有83種36科,以殼斗科(Fagaceae)、山茶科(Theaceae)、樟科(Lauraceae)等常綠植物為主。標(biāo)準(zhǔn)地內(nèi)樹木共計(jì)2 625株,其中,細(xì)葉青岡(Cyclobalanopsisgracilis)的株數(shù)最多,達(dá)665株；連蕊茶(Camelliafraterna)的株數(shù)次之,為339株；青岡(Cyclobalanopsisgracilis)、豹皮樟(Litseacoreanavar.sinensis)、短尾柯(Lithocarpusbrevicaudatus)、榧樹(Torreyagrandis)均超過100株。標(biāo)準(zhǔn)地內(nèi)樹木胸高斷面積之和為297 866.32cm2,其中,細(xì)葉青岡胸高斷面積最大,為39 834.12cm2；杉木(Cunninghamialanceolata)次之,為29 435.29 cm2；長柱紫莖(StewartiarostrataSpongb)最小,為3.14cm2。標(biāo)準(zhǔn)地內(nèi)樹冠面積之和為24 689.784 2m2,其中,細(xì)葉青岡樹冠面積最大,為5 035.330 7m2,披針葉回香(IlliciumlanceolatumA.C.Smith)樹冠面積最小,為0.196 4m2。標(biāo)準(zhǔn)地林分郁閉度為0.81。各樹種中:細(xì)葉青岡的樹種郁閉度最高,為0.17,榧樹、短尾柯、小葉青岡(Cyclobalanopsismyrsinaefolia)、青岡、楓香(Liquidambarformosana)、白櫟(Quercusfabri)、杉木(Cunninghamialanceolata)次之,在0.05～0.08之間；其余樹種均低于0.05。標(biāo)準(zhǔn)地樹冠重疊度為203.55%,其中,南京椴(TiliamiquelianaMaxim)的樹冠重疊度最高,為241.43%；紅果榆(UlmusszechuanicaFang)的樹冠重疊度次之,為207.12%；天目樸(Celtischekiangensis)和杜鵑(RhododendronsimsiiPlanch)的樹冠重疊度最低,為0.41%。

(續(xù)表)

(續(xù)表)

3.1.2常綠和落葉樹種郁閉度和樹冠重疊度分析

將不同樹種根據(jù)葉片的生理生態(tài)特征分為常綠樹種和落葉樹種,并對(duì)兩類樹種的調(diào)查因子進(jìn)行差異性分析。對(duì)于兩類樹種的樹種郁閉度、樹冠重疊度進(jìn)行正態(tài)性檢驗(yàn),結(jié)果表明均不符合正態(tài)分布(P<0.05)(表2),所以采用非參數(shù)檢驗(yàn)進(jìn)行差異性分析,結(jié)果如表3所示。常綠樹種和落葉樹種在樹種郁閉度、樹冠重疊度等兩個(gè)方面都沒有顯著性差異(P>0.05)。

表2 常綠和落葉樹種郁閉度和樹冠重疊度正態(tài)性檢驗(yàn)Tab.2 Normality test of canopy density and canopy overlapof evergreen and deciduous tree species

表3 常綠和落葉樹種郁閉度和樹冠重疊度對(duì)比分析Tab.3 Comparative analysis of canopy density and canopy overlap of evergreen and deciduous tree species

3.1.3優(yōu)勢(shì)樹種對(duì)郁閉度的貢獻(xiàn)度分析

通過優(yōu)勢(shì)度分析法確定常綠闊葉林的優(yōu)勢(shì)樹種。結(jié)果表明,常綠闊葉林的優(yōu)勢(shì)樹種共13種8科(表4),按胸高斷面積從高到低依次為細(xì)葉青岡、杉木、楓香(Liquidambarformosana)、青岡、短尾柯、榧樹、黃連木(Pistaciachinensis)、小葉青岡、藍(lán)果樹(Nyssasinensis)、櫸樹(Zelkovaschneideriana)、豹皮樟、天目木姜子(Litseaauriculata),分別屬于殼斗科、杉科(Taxodiaceae)、金縷梅科(Hamamelidaceae R.Br.)、紅豆杉科(Taxaceae)、漆樹科(Anacardiaceae)、藍(lán)果樹科(Nyssaceae)、榆科(Ulmaceae)、樟科。

分別計(jì)算每個(gè)優(yōu)勢(shì)樹種對(duì)林分郁閉度的貢獻(xiàn)度,結(jié)果如表4所示。13個(gè)優(yōu)勢(shì)樹種對(duì)林分郁閉度的貢獻(xiàn)率達(dá)90.11%,其余70個(gè)樹種對(duì)林分郁閉度的貢獻(xiàn)率僅為9.89%,故只對(duì)13個(gè)優(yōu)勢(shì)樹種的貢獻(xiàn)度進(jìn)行計(jì)算和分析。在優(yōu)勢(shì)樹種中,細(xì)葉青岡對(duì)林分郁閉度的貢獻(xiàn)度最大,為21.19%；榧樹的貢獻(xiàn)度次之,為9.36%；黃連木的貢獻(xiàn)度最小,為2.49%。不同科優(yōu)勢(shì)樹種中,殼斗科樹種對(duì)林分郁閉度的貢獻(xiàn)度最大,達(dá)51.32%。

表4 標(biāo)準(zhǔn)地優(yōu)勢(shì)樹種及其對(duì)林分郁閉度的貢獻(xiàn)度Tab.4 Dominant species of sample plot and their contribution to canopy density

3.2 林分因子與樹種郁閉度和樹冠重疊度的相關(guān)性分析

3.2.1株數(shù)與樹種郁閉度的相關(guān)性分析

對(duì)各樹種的株數(shù)與樹種郁閉度進(jìn)行相關(guān)性分析,結(jié)果如表5所示。由表5可知,株數(shù)與樹種郁閉度的Pearson相關(guān)系數(shù)為0.821,呈顯著的相關(guān)關(guān)系。由圖2可知,株數(shù)和樹種郁閉度兩者的函數(shù)關(guān)系式為y=0.0002x+0.0023,R2=0.6732。

圖2 株數(shù)和樹種郁閉度之間的關(guān)系Fig.2 Relationship between number of trees and canopydensity of tree species

表5 林分因子與樹種郁閉度和樹冠重疊度的相關(guān)性分析Tab.5 Correlation analysis of stand factors and treespecies canopy density and canopy overlap

3.2.2樹冠面積與樹種郁閉度的相關(guān)性分析

對(duì)各樹種的樹冠面積與樹種郁閉度進(jìn)行相關(guān)性分析,結(jié)果如表5所示。由表5可知,樹冠面積與樹種郁閉度的Pearson相關(guān)系數(shù)為0.962,呈顯著的相關(guān)關(guān)系。由圖3可知,樹冠面積和樹種郁閉度兩者的函數(shù)關(guān)系式為y=0.00004x-0.001,R2=0.9245。

圖3 樹冠面積與樹冠郁閉度之間的關(guān)系Fig.3 Relationship between the canopy area and canopydensity of tree species

3.2.3胸高斷面積與樹種郁閉度的相關(guān)性分析

對(duì)各樹種的胸高斷面積與樹種郁閉度進(jìn)行相關(guān)性分析,結(jié)果如表5所示。由表5可知,胸高斷面積與樹種郁閉度的Pearson相關(guān)系數(shù)為0.926,呈顯著的相關(guān)關(guān)系。由圖4可知,胸高斷面積和樹種郁閉度兩者的函數(shù)關(guān)系式為y=0.000003x-0.0018,R2=0.8576。

圖4 胸高斷面積與樹種郁閉度之間的關(guān)系Fig.4 Relationship between the chest height sectional areaand canopy density of tree species

3.2.4株數(shù) 樹冠面積 胸高斷面積 樹種郁閉度與樹冠重疊度的相關(guān)性分析

對(duì)各樹種的株數(shù)、樹冠面積、胸高斷面積、樹種郁閉度與樹冠重疊度進(jìn)行相關(guān)性分析,結(jié)果如表5所示。由表5可知,株數(shù)與樹冠重疊度的Pearson相關(guān)系數(shù)為0.047,樹冠面積與樹冠重疊度的Pearson相關(guān)系數(shù)為0.091,胸高斷面積與樹冠重疊度的Pearson相關(guān)系數(shù)為0.034,樹種郁閉度與樹冠重疊度的Pearson相關(guān)性為0.027,相關(guān)關(guān)系均未達(dá)顯著水平。

3.3 種群空間分布格局與樹冠重疊度的相關(guān)性分析

通過相關(guān)性分析,發(fā)現(xiàn)樹冠重疊度與株數(shù)、樹冠面積、胸高斷面積、樹種郁閉度之間無顯著的相關(guān)關(guān)系。為進(jìn)一步探究樹冠重疊度的影響因子,對(duì)各樹種的樹冠重疊度與聚集指數(shù)進(jìn)行相關(guān)性分析,結(jié)果如表6所示。各樹種樹冠重疊度與聚集指數(shù)的Pearson相關(guān)系數(shù)為-0.520,即兩者存在顯著的負(fù)相關(guān)關(guān)系。由圖5可知,隨著聚集指數(shù)的逐漸增大,樹冠重疊度有減小的趨勢(shì)。

圖5 樹冠重疊度與聚集指數(shù)的關(guān)系Fig.5 The relationship between crown overlapand aggregation index

表6 各樹種樹冠重疊度與聚集指數(shù)相關(guān)性分析Tab.6 Correlation analysis of crown overlap and aggregationindex of each tree species

為了驗(yàn)證這一趨勢(shì),對(duì)優(yōu)勢(shì)樹種的聚集指數(shù)進(jìn)行計(jì)算,結(jié)果如表7所示。13個(gè)優(yōu)勢(shì)樹種的聚集指數(shù)均小于1,即優(yōu)勢(shì)樹種皆為聚集分布,這與湯孟平等[8]對(duì)天目山常綠闊葉林優(yōu)勢(shì)種群的空間分布格局研究結(jié)論相一致。優(yōu)勢(shì)樹種的樹冠重疊度隨聚集指數(shù)的增大有減小的趨勢(shì),這與樹冠重疊度與聚集指數(shù)的相關(guān)性分析結(jié)論相一致。

表7 優(yōu)勢(shì)樹種聚集指數(shù)與樹冠重疊度Tab.7 Aggregation index of dominant treespecies and crown overlap

(續(xù)表)

4 結(jié)論與討論

天目山常綠闊葉林中,細(xì)葉青岡的樹種郁閉度最高,為0.17；榧樹、短尾柯、小葉青岡、青岡、楓香、白櫟、杉木次之,在0.05～0.08之間；其余樹種均低于0.05。林分的樹冠重疊度為203.55%,其中南京椴的樹冠重疊度最高,為241.43%；天目樸和杜鵑的樹冠重疊度最低,為0.41%。

天目山常綠闊葉林的優(yōu)勢(shì)樹種共有13種8科1 755株,13個(gè)優(yōu)勢(shì)樹種對(duì)郁閉度的貢獻(xiàn)率達(dá)90.11%,其中細(xì)葉青岡的貢獻(xiàn)度最大,為21.19%。以細(xì)葉青岡為主的殼斗科樹種對(duì)常綠闊葉林林分冠層的影響最大,其對(duì)林分郁閉度的貢獻(xiàn)度達(dá)51.32%。

常綠闊葉林中,株數(shù)與樹種郁閉度、樹冠面積與樹種郁閉度、胸高斷面積與樹種郁閉度之間均具有顯著的正相關(guān)關(guān)系；株數(shù)與樹冠重疊度、樹冠面積與樹冠重疊度、胸高斷面積與樹冠重疊度、樹種郁閉度與樹冠重疊度之間不存在顯著的相關(guān)關(guān)系。

林木空間分布格局是影響常綠闊葉林樹冠重疊度的重要因子,隨著聚集指數(shù)的逐漸增大,樹冠重疊度有減小的趨勢(shì)。但是,個(gè)別樹種如楓香等出現(xiàn)聚集指數(shù)較大且樹冠重疊度較大的現(xiàn)象(表7),說明樹冠重疊度不僅受林木空間分布格局的影響,也可能是多個(gè)影響因子共同作用的結(jié)果,這個(gè)問題有待進(jìn)一步研究。