周梅，王新華,李春干,楊承伶

(1.廣西大學(xué) 計(jì)算機(jī)與電子信息學(xué)院，廣西 南寧 530004；2.廣西大學(xué) 林學(xué)院，廣西 南寧 530004； 3.廣西林業(yè)勘測(cè)設(shè)計(jì)院，廣西 南寧 530011)

1 研究區(qū)域與研究方法

1.1 研究區(qū)自然概況

研究區(qū)位于廣西壯族自治區(qū)南寧市北部的高峰林場(chǎng)內(nèi)，為一個(gè)呈東北-西南走向的近矩形區(qū)域，中心地理位置為22°58′33″N，108°23′45″E，長(zhǎng)11.2km，寬4.2km，面積約4 770hm2。低山地貌，海拔90-460m，地形封閉，山不高但坡陡，一般坡度為24°-34°，坡度大于24°的面積占研究區(qū)面積的75%。研究區(qū)位于北回歸線以南，屬濕潤(rùn)亞熱帶季風(fēng)氣候，年均氣溫21.6℃左右，年均降雨量1 300mm，相對(duì)溫度79%，熱量充足、雨量充沛，林木生長(zhǎng)迅速。區(qū)內(nèi)森林約95%為人工林，除桉樹(shù)(Eucalyptus)外，年齡大多在15年以上，主要樹(shù)種為尾葉桉(Eucalyptusurophylla)、巨尾桉(Eucalyptusgrandis×E.urophylla)、馬尾松(Pinusmassoniana)、濕地松(Pinuselliottii)、杉木(Cunninghamialanceolata)、八角(Illiciumverum)、紅錐(Castanopsishystrix)、火力楠(Micheliamacclurei)、米老排(Mytilarialaosensis)、厚莢相思(Acaciacrassicarpa)等，約60%的林分為人工天然混交林。其形成過(guò)程為：由于研究區(qū)的水熱條件好，在杉木、馬尾松、尾葉桉、巨尾桉、八角人工造林后，原生樹(shù)種白楸(Mallotuspaniculatus)、木荷(Schimaargentea)等在林內(nèi)萌芽、生長(zhǎng)，形成人工天然復(fù)層混交林，也有部分為人工混交林，如杉木×馬尾松、厚莢相思×紅錐等。在溝谷地帶有少量天然雜木林。

1.2 研究方法

1.2.1 樣地設(shè)置與調(diào)查方法

將研究區(qū)森林分為4個(gè)優(yōu)勢(shì)樹(shù)種組：杉木、松樹(shù)(馬尾松、濕地松)、桉樹(shù)(尾葉桉、巨尾桉)和一般闊葉樹(shù)(紅錐、火力楠、米老排、八角、厚莢相思等和天然雜木林)，每個(gè)樹(shù)種組設(shè)置25塊以上樣地，共設(shè)置了110個(gè)樣地。樣地設(shè)置與調(diào)查方法：(1)在ArcGIS軟件環(huán)境中在小班空間分布圖上根據(jù)優(yōu)勢(shì)樹(shù)種(組)按相對(duì)均勻分布原則半隨機(jī)地布設(shè)樣地，確定樣地中心點(diǎn)的圖面坐標(biāo)。(2)根據(jù)樣地中心的圖面坐標(biāo)，采用手持雙頻差分GPS導(dǎo)航至樣地中心。觀察周圍林分樹(shù)種組成，若該中心四周30m內(nèi)林分均屬同一類型(杉木林或松林、桉樹(shù)林、一般闊葉林)，則該中心即為樣地中心；否則，移動(dòng)該中心至合適的位置；若無(wú)論如何移動(dòng)，樣地都無(wú)法只包含一個(gè)林分類型，則可涵蓋2個(gè)類型，但不能包含無(wú)林地。(3)在樣地中心西北方向20m處附近，隨機(jī)確定樣地西北角點(diǎn)，用森林羅盤儀和激光測(cè)距儀(LeicaDISTOTMX30)測(cè)設(shè)樣地邊界和分隔樣方(樣地面積為30m×30m，分為9個(gè)10m×10m的樣方)。(4)在每個(gè)樣方內(nèi)，對(duì)于胸徑大于或等于5.0cm的林木進(jìn)行每木檢尺，用超聲波測(cè)高儀(Haglf VERTEX Ⅳ)測(cè)定3株平均木和1株優(yōu)勢(shì)木的樹(shù)高。樣方計(jì)算指標(biāo)包括平均直徑、平均高、斷面積、每公頃株數(shù)、每公頃蓄積量等，其中，每公頃蓄積量根據(jù)斷面積和平均高采用形高表計(jì)算。樣地?cái)?shù)據(jù)基本情況見(jiàn)表1。

表1 900m2樣地的主要林分參數(shù)

注：杉木(Cunninghamialanceolata)；松樹(shù)(Pinus)為馬尾松(P.massoniana)、濕地松(P.elliottii)；桉樹(shù)(Eucalyptus)為尾葉桉(E.urophylla)、巨尾桉(E.grandis×E.urophylla)。

1.2.2 樣方組合及其林分參數(shù)計(jì)算方法

在900m2的樣地上，以西北角樣方(1號(hào)樣方)為固定樣方，組合成面積分別為100m2、200m2、300m2、400m2、600m2和900m2的6個(gè)不同面積大小的樣地(圖1)。

圖1 樣地設(shè)置和組合圖

根據(jù)以上組合，每個(gè)900m2樣地可以計(jì)算6個(gè)不同面積大小的樣地?cái)?shù)據(jù)，其方法是對(duì)于每個(gè)給定面積(如400m2等)的樣地，其斷面積、每公頃株數(shù)、每公頃蓄積量分別為其包含樣方相應(yīng)值之和，平均胸徑、平均高分別取其包含樣方相應(yīng)值的斷面積加權(quán)平均值，優(yōu)勢(shì)高為其包含樣方中最高優(yōu)勢(shì)木的高。

1.2.3 林木直徑變動(dòng)分析方法

由于研究區(qū)內(nèi)絕大部分林分為人工林，多經(jīng)撫育間伐等人為干擾，且大部分為復(fù)層混交林，采用SPSS 21軟件作Shapiro-Wilk檢驗(yàn)、Kolmogorov-Smirnov檢驗(yàn)和矩法檢驗(yàn)(U檢驗(yàn))——偏度和峰度檢驗(yàn)[10]，只有不到20%的樣地林木直徑分布服從正態(tài)分布，故采用變動(dòng)系數(shù)對(duì)不同面積樣地的林木直徑變動(dòng)情況進(jìn)行分析。

2 結(jié)果與分析

2.1 不同面積大小樣地的林分參數(shù)平均值比較

不同優(yōu)勢(shì)樹(shù)種(組)、不同面積樣地的林分參數(shù)的算術(shù)平均值及其變動(dòng)系數(shù)見(jiàn)表2。由表2可以看出,(1)無(wú)論是全部樣地還是各個(gè)優(yōu)勢(shì)樹(shù)種，不同面積大小樣地的平均直徑、平均高、每公頃斷面積、每公頃株數(shù)和每公頃蓄積量總平均值十分接近；(2)隨著樣地面積的增大，樣地間各個(gè)林分參數(shù)的變動(dòng)系數(shù)均呈逐漸減小的趨勢(shì)，說(shuō)明樣地面積越大林分參數(shù)的變化越穩(wěn)定。進(jìn)一步分析可以發(fā)現(xiàn),①各個(gè)優(yōu)勢(shì)樹(shù)種中，不同面積大小樣地的平均直徑、平均高、每公頃斷面積、每公頃株數(shù)和每公頃蓄積量的最大值和最小值的相對(duì)相差(最大相差)均在8%以內(nèi)，松樹(shù)、桉樹(shù)和一般闊葉樹(shù)各個(gè)參數(shù)的最大相差絕大部分小于5%(松樹(shù)、桉樹(shù)每公頃蓄積量最大相差分別為6.8%和6.1%)，杉木各個(gè)參數(shù)的最大相差略大，但均小于8.0%；②不同面積大小樣地的林分優(yōu)勢(shì)高最大相差為10.0%-19.8%，遠(yuǎn)大于其他林分參數(shù)；③在110個(gè)900m2樣地中，優(yōu)勢(shì)木出現(xiàn)在1-9樣方的概率分別為10.0%、8.2%、14.5%、8.2%、12.7%、7.3%、10.9%、12.7%、15.5%，作樣方組合后，優(yōu)勢(shì)木出現(xiàn)在100m2、200m2、300m2、400m2、600m2、900m2樣地的概率分別為10.0%、18.2%、38.2%、60.9%、100.0%，說(shuō)明隨著樣地范圍的擴(kuò)大，更高的樣木出現(xiàn)的可能性增大。

表2 不同面積樣地林分參數(shù)的平均值和變動(dòng)系數(shù)

以上結(jié)果與黃清麟等采用125塊杉木人工林樣地的分析結(jié)果相似，在該研究中，面積分別為667m2、534m2、400m2、267m2的樣地每公頃蓄積量總體平均值最大相差為2.2%，其中：幼齡林為3.7%，中齡林為1.8%，近熟林為4.7%，成過(guò)熟林為7.9%[3]。

2.2 不同面積樣地的林分參數(shù)與900m2樣地的比較

理論上，隨著樣地面積的增大，樣木株數(shù)的增多，樣地各個(gè)林分參數(shù)與所在林分相應(yīng)的參數(shù)越接近，當(dāng)樣地面積增大至與林分面積相等時(shí)，樣地的各個(gè)林分參數(shù)就與林分的相應(yīng)參數(shù)相等。在本調(diào)查中，900m2是樣地的最大面積，各個(gè)林分參數(shù)的樣地間的變動(dòng)系數(shù)亦最小，因此，可以假定其林分參數(shù)與林分真值最接近，故可用900m2樣地的各個(gè)林分參數(shù)值考察其余不同面積大小樣地相應(yīng)參數(shù)值的偏差。

對(duì)于每個(gè)優(yōu)勢(shì)樹(shù)種(組)的25-31組樣地?cái)?shù)據(jù)，逐組計(jì)算100m2、200m2、300m2、400m2、600m2樣地各項(xiàng)林分參數(shù)與900m2樣地相應(yīng)參數(shù)的相對(duì)偏差，可以發(fā)現(xiàn)：總體上，隨著樣地面積的增大，各個(gè)林分參數(shù)的偏差均呈逐漸減小的趨勢(shì)，600m2樣地的偏差最小。圖2表示4個(gè)優(yōu)勢(shì)樹(shù)種(組)不同面積樣地每公頃蓄積量與900m2樣地每公頃蓄積量的相對(duì)偏差曲線。由該圖可以看出，各個(gè)優(yōu)勢(shì)樹(shù)種(組)中，600m2樣地的偏差曲線與偏差為0的水平直線最接近。

圖2 不同面積樣地與900m2樣地每公頃蓄積量的相對(duì)偏差曲線

將上述相對(duì)偏差取絕對(duì)值后按優(yōu)勢(shì)樹(shù)種、樣地面積分別取平均值，得到各個(gè)優(yōu)勢(shì)樹(shù)種不同面積大小樣地林分參數(shù)與900m2樣地相應(yīng)參數(shù)的平均相對(duì)偏差(表3)。

由表3可以看出，當(dāng)樣地面積達(dá)到600m2時(shí)，平均直徑、平均高、優(yōu)勢(shì)高和每公頃斷面積及桉樹(shù)每公頃株數(shù)、每公頃蓄積量的平均相對(duì)偏差均小于5%，杉木、松樹(shù)和一般闊葉樹(shù)每公頃株數(shù)、每公頃蓄積量的偏差略大，但均小于8%；隨著樣地面積的增大，各項(xiàng)林分參數(shù)與900m2樣地的最大相對(duì)偏差亦呈逐漸減小的趨勢(shì)，尤其是每公頃蓄積量和每公頃株數(shù)，減少的幅度最大，600m2樣地的最大偏差比100m2樣地減小了50%以上。

進(jìn)一步分析發(fā)現(xiàn)，隨著樣地面積的增大，各個(gè)林分參數(shù)與900m2樣地相應(yīng)參數(shù)偏差大的樣地?cái)?shù)量逐漸減少。以每公頃蓄積量為例，當(dāng)樣地面積為100m2時(shí)，杉木、松樹(shù)、桉樹(shù)和一般闊葉樹(shù)偏差大于10%的樣地?cái)?shù)量比重分別為72.0%、65.5%、64.0%和74.2%，當(dāng)樣地面積增大至600m2時(shí)，這一比重分別減小至12.0%、20.7%、8.0%和25.8%。

2.3 不同面積樣地對(duì)其林木直徑變動(dòng)的影響

理論上，當(dāng)樣地面積的增大、林木株數(shù)增多時(shí)，林木直徑的變動(dòng)將趨于穩(wěn)定，變動(dòng)系數(shù)變小。但分析結(jié)果表明，無(wú)論是全部樣地，還是各個(gè)優(yōu)勢(shì)樹(shù)種(組)，林木直徑變動(dòng)系數(shù)并未隨著樣地面積的增大而變小，而是呈現(xiàn)出無(wú)規(guī)律的變化，但總平均值十分接近，見(jiàn)表4。

表4 不同面積樣地的林木直徑平均變動(dòng)系數(shù)

進(jìn)一步分析發(fā)現(xiàn)，盡管不同面積大小樣地的林木直徑變動(dòng)系數(shù)變化的規(guī)律性不強(qiáng)，但總體存在100m2樣地的直徑變化系數(shù)大于900m2樣地的情況，但隨著樣地面積的增大，樣地間的變動(dòng)系數(shù)呈減小的趨勢(shì)。以杉木為例，其林木直徑變動(dòng)系數(shù)的變化情況見(jiàn)圖3。

圖3 不同面積大小杉木樣地林木直徑變動(dòng)系數(shù)的變化

以上樣地林木直徑的變動(dòng)系數(shù)無(wú)規(guī)律的變化，可能與研究區(qū)大多為混交林，且除桉樹(shù)外都經(jīng)歷過(guò)撫育間伐有關(guān)。

3 結(jié)論與討論

在本研究中，除優(yōu)勢(shì)高外，不同面積大小樣地的林分參數(shù)總平均值十分接近，全部樣本各個(gè)參數(shù)的最大相差均小于4%，每公頃蓄積量的最大相差只有2.9%，各個(gè)優(yōu)勢(shì)樹(shù)種全部林分參數(shù)的最大相差均小于8%，大多數(shù)參數(shù)的最大相差小于5%，說(shuō)明對(duì)于總體平均數(shù)的估計(jì)而言，100m2的樣地與900m2的樣地具有相同或相近的結(jié)果。因此，在采用抽樣方法進(jìn)行森林資源調(diào)查時(shí)，可考慮設(shè)置較小面積的樣地，如300-600m2，在確保抽樣精度的前提下，盡可能減少調(diào)查工作量，提高工作效率。上述結(jié)果與黃清麟等的研究中建議采用400m2樣地的結(jié)論[3]相近。

雖然不同面積大小樣地的林分參數(shù)總平均值十分接近，但在同一組樣地中，不同面積大小樣地林分參數(shù)仍然存在較大差異，并且隨著樣地面積的增大，樣地林分參數(shù)與900m2樣地的偏差逐漸減小，600m2樣地與900m2樣地的林分參數(shù)的平均絕對(duì)相對(duì)相差均不超過(guò)7%，大部分小于5%。此外，不同面積樣地的林木直徑的變動(dòng)系數(shù)表現(xiàn)為無(wú)規(guī)律的變化。因此，在作林分參數(shù)建模樣地設(shè)置時(shí)，600m2的樣地也可能取得較好結(jié)果。這一結(jié)果與司洪生[2]編制杉木立地指數(shù)表、龍紅等[4]編制濕地松數(shù)表建議采用的800m2樣地有較大差異。

本研究的樣地?cái)?shù)量偏少，樣地最大面積亦不足夠大，抽樣方法亦非完全隨機(jī)抽樣，因此，其結(jié)論的可靠性仍有待進(jìn)一步的驗(yàn)證。

致謝：廣西壯族自治區(qū)林業(yè)勘測(cè)設(shè)計(jì)院梁耀、黃孝發(fā)、何黃之園等15人，廣西大學(xué)林學(xué)院陸道調(diào)等4人參與樣地調(diào)查，在此謹(jǐn)表謝意。

參考文獻(xiàn)：

[1]雷相東,洪玲霞,陸元昌,等.國(guó)家級(jí)森林資源清查地面樣地設(shè)計(jì)[J].世界林業(yè)研究,2008,21(4):35-40.

[2]司洪生.對(duì)編制杉木立地指數(shù)表時(shí)樣地面積和優(yōu)勢(shì)木株數(shù)的探討[J].林業(yè)科技通訊,1981(2):22-24.

[3]黃清麟,陳平留,林金國(guó).杉木樣地調(diào)查面積初探[J].中南林業(yè)調(diào)查規(guī)劃,1992(2):18-21.

[4]龍紅,周錫祥,沈家智,等.濕地松數(shù)表研究中樣地面積大小的確定[J].江西農(nóng)業(yè)大學(xué)學(xué)報(bào),1993(S3):119-122.

[5]周紅敏,惠剛盈,趙中華,等.森林結(jié)構(gòu)調(diào)查中最適樣方面積和數(shù)量的研究[J].林業(yè)科學(xué)研究,2009,22(4):492-495.

[6]魏士忠,竇宏海,劉學(xué)東,等.立木調(diào)查中樣地的面積、數(shù)量與抽樣精度的關(guān)系[J].河北林果研究,2013,28(1):41-43.

[7]湯孟平,周國(guó)模,施擁軍,等.天目山常綠闊葉林群落最小取樣面積與物種多樣性[J].浙江林學(xué)院學(xué)報(bào),2006,23(4):357-361.

[8]陳宏偉,劉永剛,馮弦.西南樺人工林群落取樣面積探討[J].云南林業(yè)科技,1999(3):24-27.

[9]孫麗娟,張國(guó)蓉,劉濤,等.林分平均直徑100倍圓法與角規(guī)、標(biāo)準(zhǔn)地法測(cè)樹(shù)特征的對(duì)比分析[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2014,34(6):1-6.

[10]劉慶武,胡志艷.如何用SPSS、SAS統(tǒng)計(jì)軟件進(jìn)行正態(tài)性檢驗(yàn)[J].湘南學(xué)院學(xué)報(bào)(自然科學(xué)版),2005,7(3):56-58.