李 元

(遼寧省林業(yè)調(diào)查規(guī)劃監(jiān)測院 遼寧,沈陽 110122)

林分空間結(jié)構(gòu)是林分特征的重要組成部分,是森林經(jīng)營的理論基礎(chǔ)[1,2]。林分空間結(jié)構(gòu)很大程度上決定森林功能,如森林生態(tài)系統(tǒng)穩(wěn)定性、林下植被多樣性維持、野生動物棲息地等,合理的空間結(jié)構(gòu)是充分發(fā)揮森林木材和生態(tài)功能的基礎(chǔ)[3-5]。同時(shí),林分空間結(jié)構(gòu)還決定了種群關(guān)系,如競爭優(yōu)勢及其空間生態(tài)位,反映了森林喬木種群物種的空間關(guān)系[6,7],空間結(jié)構(gòu)研究在現(xiàn)代森林經(jīng)營理論研究和應(yīng)用具有重要的意義[8]。林分空間結(jié)構(gòu)包括多個(gè)方面,但主要基于坐標(biāo)位置的三個(gè)方面進(jìn)行描述(多樣性、競爭、分布格局):樹種(多樣性)在空間隔離程度即混交度,表示樹種多樣性基于位置的分布狀態(tài)。林木相對優(yōu)勢程度即競爭和林木個(gè)體在水平地面上的分布形式,描述了林分的競爭優(yōu)勢或生態(tài)位。林木空間分布格局,量化了林木在林地中的集聚程度[7-9]。

長白山林區(qū)是我國重要的天然林分布區(qū),同時(shí)也時(shí)重要的木材生產(chǎn)基地之一,發(fā)揮著重要的經(jīng)濟(jì)和生態(tài)效應(yīng)。經(jīng)過多年連續(xù)的人為采伐,長白山森林生態(tài)系統(tǒng)受到破壞,天然林資源急劇減少[10,11],大面積的天然林轉(zhuǎn)變?yōu)榇紊?。次生林是原始森林?jīng)歷人為反復(fù)不合理采伐或自然災(zāi)害后,在次生裸地上經(jīng)次生演替或干擾后跡地上殘存林木而形成的次生群落[12],與原始林相比,結(jié)構(gòu)差,功能弱。在林分內(nèi)單株間成熟度、生長差異大,林分的生產(chǎn)力低下,森林生態(tài)系統(tǒng)脆弱。本文通過對現(xiàn)有的長白山區(qū)次生林空間結(jié)構(gòu)狀態(tài)進(jìn)行研究分析,制定合理的森林經(jīng)營建議,優(yōu)化林分空間結(jié)構(gòu),進(jìn)而為長白山天然次生林的可持續(xù)發(fā)展提供科學(xué)的理論基礎(chǔ)。

1 研究區(qū)概況

研究區(qū)位于遼寧省開原市黃旗寨滿族鄉(xiāng),地理坐標(biāo)42°10′-42°15′ N,124°22′-124°26′ E,屬于長白山支脈。該區(qū)植被屬于長白山植物區(qū)系,主要樹種有紅松、樟子松、落葉松、云冷杉等針葉,水曲柳、胡桃楸、黃菠蘿等闊葉。由地帶性頂級群落闊葉紅松林或云冷杉混交林,經(jīng)人為干擾后形成的典型次生林。

2 數(shù)據(jù)來源

本論文所用數(shù)據(jù)來自在黃旗寨滿族鄉(xiāng)設(shè)置的30 m×30 m的樣地。樣地調(diào)查時(shí),采用相鄰格子調(diào)查法,即把樣地均勻的劃分為36個(gè)5 m×5 m的調(diào)查單元。起測徑階5cm,對樣地內(nèi)所有樹木(胸徑≥5 cm)進(jìn)行每木檢尺。測定每5 m×5 m調(diào)查單元中的每木樹種、胸徑(DBH)、樹高(H)、枝下高(BH)和冠幅(CW)等測樹傳統(tǒng)因子,并記錄下每株樹木的平面坐標(biāo)位置,建立數(shù)據(jù)庫進(jìn)行研究分析。調(diào)查統(tǒng)計(jì)結(jié)果如表1所示。

表1 樣地每木調(diào)查結(jié)果統(tǒng)計(jì)表

3 研究方法

3.1 林分空間結(jié)構(gòu)單元

本文通過構(gòu)建Voronoi多邊形和Delaunay三角網(wǎng)可以快速、簡潔、合理的確定對象木的鄰近木,在此基礎(chǔ)上計(jì)算混交度、大小比數(shù)和角尺度。林地中任一單株木和距它最近的相鄰木(基于坐標(biāo)的絕對距離),構(gòu)成了林分空間結(jié)構(gòu)的基本單元,空間結(jié)構(gòu)分析便是基于一個(gè)個(gè)結(jié)構(gòu)單元的特征值組合分析[13]。林分空間結(jié)構(gòu)單元中心或林分內(nèi)任一的那株樹被稱為對象木,一般是經(jīng)營中的目標(biāo)樹或目的樹種,或頂級演替階段樹種,而最近的相鄰木則被稱為相鄰木或鄰近木,通常為目標(biāo)樹周圍的伴生樹[13]。分析時(shí),關(guān)鍵的問題是確定近鄰木的個(gè)數(shù)[14],相鄰木的個(gè)數(shù)決定了空間結(jié)構(gòu)的大小[1]。本文確定n值的原則是:以樣地調(diào)查單株木的相對坐標(biāo)視為平面點(diǎn),平面點(diǎn)即為現(xiàn)實(shí)林地中的樹木,根據(jù)Voronoi圖點(diǎn)構(gòu)造和面空間劃分最近性、鄰接性的特點(diǎn),本文利用ArcGIS中的Voronoi工具,創(chuàng)建Voronoi多邊形和構(gòu)造Delaunay三角網(wǎng)[15]:通?，F(xiàn)實(shí)林地中心木的鄰近木株數(shù)n等于二維平面上Voronoi多邊形的邊數(shù)[16];對象木或中心木與最近鄰木的距離,即現(xiàn)實(shí)林木距離,等于二維平面上Delaunay三角網(wǎng)的邊長,二位平面上Delaunay三角網(wǎng)中兩邊的夾角即為現(xiàn)實(shí)林地中任意2個(gè)鄰近木的夾角;相關(guān)研究Voronoi多邊形和Delaunay三角網(wǎng),與鄰近木取4時(shí)的傳統(tǒng)方法在鄰近木和空間結(jié)構(gòu)單元的比較分析表明:兩者在計(jì)算混交度、大小比數(shù)時(shí)有較高的一致性[17]。

3.2 邊緣效應(yīng)

林分中每株樹競爭區(qū)域的大小并不相同,近似圓形區(qū)域的半徑與對象木及林分中最大優(yōu)勢木的大小有關(guān),在數(shù)值上等于待估對象木影響圈半徑與最大優(yōu)勢木影響圈半徑之和[18,19]。分析時(shí),鄰近木的個(gè)數(shù)存在一定程度的偏差或低估,通常樣地設(shè)置時(shí),某一中心木競爭區(qū)域的半徑大于到樣地邊界的距離,它的競爭區(qū)域便是不完全的[19]。在計(jì)算林木空間結(jié)構(gòu)參數(shù)時(shí),若忽略樣地邊緣樹木對內(nèi)部林木鄰近木選擇時(shí)的影響,即處在樣地邊界附近的林木結(jié)果對林分結(jié)果的誤差,計(jì)算的結(jié)果便會有誤差,因此分析前通常對標(biāo)準(zhǔn)地進(jìn)行邊緣的校正,去除一定范圍的邊界木,使計(jì)算的結(jié)果更準(zhǔn)確[20]。

為避免邊緣效應(yīng)對鄰近木選取及林分空間結(jié)構(gòu)參數(shù)計(jì)算結(jié)果的影響,本研究根據(jù)數(shù)據(jù)特性和樣地面積,采用緩沖區(qū)校正法,即在楊迪邊界處設(shè)置一定距離的緩沖區(qū)。校正時(shí),邊緣效應(yīng)帶的寬度一般為平均優(yōu)勢木影響圈半徑的2倍,并且林分的樹冠的半徑一般很少超過5m[19,20]。根據(jù)這一研究結(jié)論,分析時(shí)通常事先在軟件上劃分出5 m寬的界限,作為緩沖區(qū),沖緩沖中的林木不作為中心木或?qū)ο竽?而只作為鄰近木。這一思想與“8”鄰域邊緣校正類似[20,21]。因此本研究中,標(biāo)準(zhǔn)地邊界設(shè)立5 m的邊界,作為緩沖區(qū)[20]。邊緣校正如圖1。

圖1 邊緣效應(yīng)緩沖區(qū)

3.3 林分空間結(jié)構(gòu)指數(shù)的計(jì)算

本研究選擇體現(xiàn)樹種多樣性基于空間位置的混交度(M)[21]、反映林木個(gè)體大小分化即競爭、生態(tài)位的大小比數(shù)(U)[22]以及描述林木個(gè)體在林地中的分布格局的角尺度(W)[23]三個(gè)參數(shù)作為研究對象。

(1)混交度

混交度用來說明樹種多樣性基于空間位置的隔離程度,被定義為與對象木不屬同種的個(gè)體的鄰近木占所有鄰近木(通常為4)的比例[24]。計(jì)算公式為:

(1)

式中:Mi為對象木i的混交度。Vij是一個(gè)離散性的變量(不連續(xù)變量)。Mi=0表示中心木周圍的鄰近木均為同種樹木,目標(biāo)樹或中心的樹種隔離高。相反,Mi=1則表示鄰近木均屬他種[2]。一般認(rèn)為樹種隔離程度越高,林分的結(jié)構(gòu)越穩(wěn)定混交度量化了任意一株樹的最近相鄰木為同種或它種的概率程度大小[25]。當(dāng)對象木i與第j株相鄰木非同種時(shí),Vij=1,反之,Vij=0,顯然,0≤Mi≤1[16]。

(2)大小比數(shù)

大小比數(shù)通?？梢詫α帜久磕緶y定指標(biāo)(胸徑、樹高、冠幅等)的狀態(tài)進(jìn)行描述,表示基于位置的中心木和鄰近木每木測定指標(biāo)的相對大小[26],大小比數(shù)被定義為大于參照樹的相鄰木數(shù)占所有鄰近木的比例[23],公式如下:

(2)

式中:Ui為中心木i的大小比數(shù)。大小比數(shù)量化了參照樹與其相鄰木,基于空間位置的競爭關(guān)系或中心木周圍比中心木每木特征值(樹高、胸徑、官府等)大的個(gè)數(shù)[2]。Ui越低,說明比中心木周圍的大的鄰近木越少,中心木的競爭優(yōu)勢越大[1]。本文中,選定大小比數(shù)的分析參數(shù)值為胸徑。大小比數(shù)的值越大,表示比中心木胸徑大的相鄰木越多,則中心木競爭處于劣勢Kij是一個(gè)離散型的變量,當(dāng)相鄰木j比參照樹i小,Kij=0;否則Kij=1[16]。一般認(rèn)為大小比數(shù)越高,目標(biāo)樹或目的樹種相鄰木個(gè)體越占優(yōu)勢,而相應(yīng)地參照樹越不占優(yōu)勢。

(3)角尺度

角尺度用來描述中心木周圍鄰近木的分布均勻性或聚集性[6]。通常角尺度(Wi)被定義為a角小于標(biāo)準(zhǔn)角a0的個(gè)數(shù)占所有角的比例[2],用下式表示:

(3)

式中:Wi為對象木i的角尺度。對中心木i的n個(gè)鄰近木而言,中心木處于均勻分布狀態(tài)時(shí),其分布角應(yīng)各為360°/n[15]。Zij是一個(gè)離散型的變量,當(dāng)?shù)趈個(gè)a角小于標(biāo)準(zhǔn)角時(shí),Zij=1,相反Zij=0[16]。Wi越接近0.5,林木水平分布格局越理想。

4 結(jié)果與分析

4.1 Voronoi圖和Delaunay三角網(wǎng)構(gòu)建

在研究區(qū)域內(nèi),樣地中共有195株林木,其中在核心區(qū)域內(nèi),共測試113株,將調(diào)查的樣地每株樹木點(diǎn)位置數(shù)據(jù)加載進(jìn)ArcGIS軟件,每株樹木的平面坐標(biāo)可以在ArcGIS中確定一個(gè)點(diǎn),設(shè)置緩沖區(qū),創(chuàng)建基于空間位置的樣地中林木點(diǎn)的Voronoi多邊形和構(gòu)造Delaunay三角網(wǎng)(圖2);根據(jù)Voronoi多邊形特征,確定對象木的最近鄰木株數(shù)N,最近鄰木數(shù)取值4、5和6較多,平均取值為5(圖3)。

圖2 Voronoi圖與Delaunay三角網(wǎng)

圖3 最近鄰木株數(shù)分布

4.2 空間結(jié)構(gòu)指數(shù)的計(jì)算

利用ArcGIS軟件構(gòu)建樹木點(diǎn)的Voronoi圖和Delaunay三角網(wǎng),確定鄰近木株數(shù)。利用式(1)、(2)、(3)計(jì)算林分個(gè)體的空間結(jié)構(gòu)指數(shù),計(jì)算結(jié)果如表2所示。林分的平均混交度只有0.38,屬于弱度混交,這是由于樣地中樹種種類比較單一,并且花曲柳株數(shù)所占的比例較大;林分的平均大小比數(shù)為0.47,根據(jù)大小比數(shù)的定義,處于中庸?fàn)顟B(tài),這與樣地為天然次生的實(shí)際情況相符;平均角尺度在[0.475,0.517]為隨機(jī)分布,依照空間分布格局的判定標(biāo)準(zhǔn)[27],小于0.475為均勻分布,大于0.517為聚集分布,樣地的平均角尺度為0.55,所以樣地的林木水平分布格局屬于聚集分布。

表2 基于Voronoi圖與Delaunay三角網(wǎng)空間結(jié)構(gòu)指數(shù)結(jié)果

5 結(jié)論

利用ArcGIS的空間分析功能,對標(biāo)準(zhǔn)地的每木調(diào)查數(shù)據(jù),構(gòu)建Voronoi圖和Delaunay三角網(wǎng),可以快速、直觀、靈活、合理的確定對象木的最近鄰接木,并確定出最近鄰接木的株數(shù)N的取值,在此基礎(chǔ)上計(jì)算空間結(jié)構(gòu)指數(shù)。

由計(jì)算結(jié)果得出,標(biāo)準(zhǔn)地角尺度為0.55,林分的空間分布格局為聚集分布;林分的混交度僅為0.38,林分處于弱度混交;林分的大小比數(shù)為0.47,林分內(nèi)的林木個(gè)體大小分化程度屬于中庸?fàn)顟B(tài)。我們的研究結(jié)果表明,本文所研究的林分為天然次生林,其呈現(xiàn)的是一個(gè)林木個(gè)體團(tuán)狀分布、弱度混交狀態(tài)的森林群落。建議在森林資源經(jīng)營中,根據(jù)林分空間結(jié)構(gòu)分析結(jié)果制定相應(yīng)的營林措施,以期逐步優(yōu)化林分的空間結(jié)構(gòu),達(dá)到穩(wěn)定健康狀態(tài)。