毛竹林空間結(jié)構(gòu)與更替動(dòng)態(tài)的關(guān)系

唐思嘉， 湯孟平， 沈錢勇， 杜秀芳， 龐春梅

（1.浙江農(nóng)林大學(xué) 省部共建亞熱帶森林培育國(guó)家重點(diǎn)實(shí)驗(yàn)室，浙江 杭州311300；2.浙江農(nóng)林大學(xué) 環(huán)境與資源學(xué)院，浙江 杭州311300；3.浙江天目山國(guó)家級(jí)自然保護(hù)區(qū) 管理局， 浙江 杭州311311）

毛竹Phyllostachys edulis林是中國(guó)南方一種特殊的森林類型。中國(guó)是毛竹的主產(chǎn)地，毛竹林面積達(dá)443萬hm2，占所有竹林面積的74%［1］。毛竹生長(zhǎng)快、周期短、產(chǎn)量高、用途廣［2］，同時(shí)具有較強(qiáng)的固碳能力，在應(yīng)對(duì)氣候變化上有不可低估的作用［3－4］。森林空間結(jié)構(gòu)是森林生長(zhǎng)過程的驅(qū)動(dòng)因子，對(duì)森林未來的發(fā)展具有決定性作用［5－6］。毛竹林的結(jié)構(gòu)在一定程度上決定了其穩(wěn)定性、發(fā)展的可能性與經(jīng)營(yíng)空間的大?。?］。近年來，對(duì)毛竹林經(jīng)營(yíng)管理的認(rèn)識(shí)已從立竹度、年齡結(jié)構(gòu)和樹種組成等非空間結(jié)構(gòu)向毛竹空間分布格局、年齡隔離和競(jìng)爭(zhēng)（或大小空間排列）等空間結(jié)構(gòu)［7－8］逐漸深化。毛竹依賴多年生竹鞭發(fā)筍繁殖，出筍周期2 a；竹筍出土50 d后高度生長(zhǎng)和直徑生長(zhǎng)基本完成［9］，生長(zhǎng)速率大于喬木林。目前，對(duì)喬木林空間結(jié)構(gòu)與其生長(zhǎng)規(guī)律的研究較多［10－12］，而關(guān)于毛竹林空間結(jié)構(gòu)及林分更替動(dòng)態(tài)關(guān)系的研究鮮見報(bào)道。本研究以浙江天目山國(guó)家級(jí)自然保護(hù)區(qū)內(nèi)少受人為干擾的毛竹為研究對(duì)象，基于地理信息系統(tǒng)（geographic information system,GIS）的維諾（Voronoi）圖空間分析功能確定空間結(jié)構(gòu)單元［13］，選用最近鄰竹株數(shù)（nearest bamboo number）， 角尺度（uniform angle index）， 大小比數(shù)（neighborhood comparison）， 競(jìng)爭(zhēng)指數(shù)（competition index）等4個(gè)空間結(jié)構(gòu)指數(shù)，采用主成分分析方法，研究毛竹林空間結(jié)構(gòu)與更替動(dòng)態(tài)的關(guān)系，旨在為毛竹林經(jīng)營(yíng)提供科學(xué)依據(jù)。

1 研究區(qū)概況

浙江天目山國(guó)家級(jí)自然保護(hù)區(qū)位于浙江省西北部杭州市臨安區(qū)境內(nèi)的西天目山（30°18′30″～30°24′55″N,119°23′47″～119°28′27″E）， 總面積為 1 050 hm2， 主峰海拔為 1 506 m。 年平均氣溫為 8.8～14.8 ℃，年積溫為2 500～5 100℃，年降水量為1 390～1 870 mm，相對(duì)濕度為76%～81%。受海洋暖濕氣候影響，該區(qū)森林植被豐茂，呈垂直帶狀分布：海拔870 m以下為常綠闊葉林，870～1 100 m為常綠、落葉闊葉混交林，1 100～1 380 m為落葉闊葉林，1 380～1 500 m為落葉矮林。毛竹林鑲嵌于其他森林類型之間，多分布在海拔350～900 m，其林下植被稀少，主要樹種有豹皮樟Litsea coreana var.sinensis，連蕊茶Camellia fraterna，細(xì)葉青岡Cyclobalanopsis myrsinaefolia，微毛柃Eurya hebeclados和短尾柯Lithocarpus brevicaudatus等［14］。該區(qū)多年來堅(jiān)持 “不動(dòng)一草一木”的方針，區(qū)內(nèi)毛竹林少有人為干擾，毛竹生長(zhǎng)狀況接近自然狀態(tài)［15］。

2 研究方法

2.1 研究數(shù)據(jù)

2.1.1 標(biāo)準(zhǔn)地調(diào)查 2009年7月，在浙江天目山國(guó)家級(jí)自然保護(hù)區(qū)內(nèi)，選擇毛竹林典型地段，建立1塊大小為100 m×100 m的固定標(biāo)準(zhǔn)地。中心海拔為840 m，最高處海拔886 m，最低處海拔810 m；主坡向南偏東30°，平均坡度為37°；下部平緩，上部相對(duì)較急，水平方向相對(duì)一致（圖1）；土壤類型為黃壤，壤質(zhì)，土層厚度中。利用南方全站儀NTS355設(shè)定標(biāo)準(zhǔn)地的坐軸，采用相鄰網(wǎng)格法將整個(gè)標(biāo)準(zhǔn)地劃分為100個(gè)10 m×10 m的樣方；測(cè)定標(biāo)準(zhǔn)地內(nèi)每株毛竹的基部三維坐標(biāo)（x，y，z）。用枝痕法、觀桿法確定毛竹年齡，用號(hào)竹法在胸高位置標(biāo)記毛竹年份和竹號(hào)。每木調(diào)查測(cè)定毛竹胸徑、竹高、年齡、枝下高、冠幅、彎曲狀態(tài)等因子。2010－2015年，每年4月調(diào)查毛竹林固定標(biāo)準(zhǔn)地內(nèi)的新筍地徑，7月復(fù)查。

圖1 標(biāo)準(zhǔn)地三維地形Figure 1 3D map of the sample plot

2.1.2 各大年毛竹數(shù)據(jù)統(tǒng)計(jì) 根據(jù)各年毛竹林標(biāo)準(zhǔn)地復(fù)查數(shù)據(jù)，確定當(dāng)年出筍長(zhǎng)成的毛竹為新竹，當(dāng)年?duì)顟B(tài)為存活而翌年死亡的毛竹為死竹。因毛竹生長(zhǎng)周期存在大小年現(xiàn)象，小年很少出筍或不出筍［16］。本試驗(yàn)以毛竹林各大年（2010，2012和2014年）數(shù)據(jù)為依據(jù)，對(duì)毛竹林更替狀況與空間結(jié)構(gòu)的關(guān)系進(jìn)行研究。各大年活竹、新竹和死竹的統(tǒng)計(jì)數(shù)據(jù)見表1。

表1 各大年毛竹統(tǒng)計(jì)數(shù)據(jù)Table 1 Statistical data of Phyllostachys edulis in 2010,2012 and 2014

2.2 空間結(jié)構(gòu)單元及其邊緣矯正

以標(biāo)準(zhǔn)地內(nèi)任意1株毛竹為對(duì)象竹，基于Voronoi圖確定由對(duì)象竹與相鄰竹組成的空間結(jié)構(gòu)單元［17］。每個(gè)Voronoi多邊形內(nèi)包含的1株毛竹即為對(duì)象竹，與其相鄰Voronoi多邊形內(nèi)的毛竹即為最近相鄰竹。對(duì)標(biāo)準(zhǔn)地采用八鄰域平移法作邊緣校正；計(jì)算空間結(jié)構(gòu)指數(shù)時(shí)，矯正標(biāo)準(zhǔn)地內(nèi)的毛竹為對(duì)象竹。

2.3 空間結(jié)構(gòu)指數(shù)

2.3.1 角尺度 從對(duì)象竹出發(fā)，測(cè)定任意2株相鄰竹之間的2個(gè)夾角；設(shè)小角為α，最近鄰竹均勻分布時(shí)的夾角為標(biāo)準(zhǔn)角α0。根據(jù)GADOW 等［18］和 POMMERENING 等［19－20］研 究 ， 角 尺度標(biāo)準(zhǔn)角隨著最近鄰株數(shù)n的變化而變化，取值為360°/（n+1）。角尺度（Wi）被定義為α角小于標(biāo)準(zhǔn)角的個(gè)數(shù)占所考察的最近鄰株數(shù)n的比例。計(jì)算公式為：

式（1）中：Wi為對(duì)象竹i的Voronoi角尺度，n為最近鄰竹株數(shù)，αij為離散變量，當(dāng)?shù)趈個(gè)α角小于α0時(shí)， αij＝1， 反之， αij＝0。

2.3.2 大小比數(shù) 大小比數(shù)是表達(dá)樹種空間優(yōu)勢(shì)程度的良好指標(biāo)，是對(duì)大小分化度的改進(jìn)。大小比數(shù)（Ui）被定義為大于對(duì)象竹的相鄰竹株數(shù)占最近相鄰竹總株數(shù)的比例［21－22］?？蛇x擇胸徑、樹高和冠幅等作為比較指標(biāo)［23］。本研究中以胸徑作為比較指標(biāo)。大小比數(shù)計(jì)算公式為：

式（2）中：Ui為對(duì)象竹i的大小比數(shù)，n為最近鄰竹株數(shù)，Kij為離散變量，當(dāng)相鄰竹j胸徑大于對(duì)象竹i時(shí)， Kij＝1， 反之， Kij＝0。

2.3.3 競(jìng)爭(zhēng)指數(shù) 競(jìng)爭(zhēng)指數(shù)被廣泛應(yīng)用并被證明能夠很好地解釋植物競(jìng)爭(zhēng)的強(qiáng)度、作用和競(jìng)爭(zhēng)結(jié)果［23］。Hegyi競(jìng)爭(zhēng)指數(shù)在形式上反映了林木個(gè)體競(jìng)爭(zhēng)與生存空間的關(guān)系，方法簡(jiǎn)便易行［24］。因此，本研究采用Hegyi競(jìng)爭(zhēng)指數(shù)計(jì)算毛竹林標(biāo)準(zhǔn)地內(nèi)各對(duì)象竹的競(jìng)爭(zhēng)指數(shù)。計(jì)算公式為：

式（3）中：Ci為對(duì)象竹i的競(jìng)爭(zhēng)指數(shù)，n為最近鄰竹株數(shù)，di，dj分別為對(duì)象竹i和相鄰竹j胸徑，Lij為對(duì)象竹i和相鄰竹j之間的距離。

2.4 方差分析

方差分析（analysis of variance,ANOVA）是用于2個(gè)及2個(gè)以上樣本均數(shù)差別的顯著性檢驗(yàn)。對(duì)毛竹各空間結(jié)構(gòu)指數(shù)中活竹、新竹和死竹的均值比較采用單因素方差分析（one-way ANOVA），且兩兩比較時(shí)選擇最小顯著性差異（least significant difference,LSD）法。

2.5 主成分分析

主成分分析（principal component analysis,PCA）是一種簡(jiǎn)化數(shù)據(jù)集的方法。通過對(duì)毛竹林空間結(jié)構(gòu)指數(shù)作主成分分析可以確定影響毛竹林更替的主導(dǎo)因子。對(duì)毛竹林空間結(jié)構(gòu)指數(shù)進(jìn)行標(biāo)準(zhǔn)化；然后計(jì)算空間結(jié)構(gòu)指數(shù)的相關(guān)矩陣特征值、貢獻(xiàn)率、累積貢獻(xiàn)率以及其各主成分；最后對(duì)各主成分進(jìn)行解釋，并進(jìn)行重要性排序［25］。本研究采用SPSS for Windows軟件計(jì)算各空間結(jié)構(gòu)指數(shù)的相關(guān)矩陣特征值、貢獻(xiàn)率、累積貢獻(xiàn)率和載荷矩陣。

3 結(jié)果與分析

3.1 最近鄰竹株數(shù)與毛竹更替的關(guān)系

對(duì)象竹的最近鄰竹對(duì)其更替有最直接的影響［25］?；贕IS的Voronoi圖可以確定各對(duì)象竹的最近鄰竹株數(shù)。結(jié)果表明（表2）：在各大年，活竹、新竹和死竹的最近鄰竹株數(shù)頻率分布中5～7株的頻率最多；各年新竹或死竹的平均最近鄰竹株數(shù)與活竹平均最近鄰竹株數(shù)接近，均為6株左右。單因素方差分析發(fā)現(xiàn)，各年活竹、新竹和死竹在最近鄰竹株數(shù)上不存在顯著差異（P＞0.05），表明三者具有一致的分布特征，即最近鄰竹株數(shù)分布不影響竹林更替，毛竹林分具有相對(duì)穩(wěn)定性。

表2 各大年毛竹活竹、新竹和死竹最近鄰竹株數(shù)統(tǒng)計(jì)Table 2 Statistics of nearest bamboo number of live bamboo,new bamboo and dead bamboo in 2010,2012 and 2014

對(duì)比新竹與活竹、死竹與活竹的最近鄰竹株數(shù)頻率分布（圖2）可以看出：新竹在 “3株” “4株”或 “7株以上”的最近鄰竹株數(shù)頻率要高于活竹；而死竹在 “7株以上”的最近鄰竹株數(shù)頻率要高于活竹。說明新竹更容易在較低密度或較高密度區(qū)域內(nèi)發(fā)筍生長(zhǎng)，而死竹主要出現(xiàn)于較高密度區(qū)域。因此，為減少毛竹生長(zhǎng)過程中的自然損耗，在毛竹林經(jīng)營(yíng)過程中應(yīng)對(duì)毛竹林高密度（最近鄰竹株數(shù)＞7株）區(qū)域進(jìn)行擇伐。

圖2 各大年新竹與活竹、死竹與活竹的最近鄰竹株數(shù)分布差異Figure 2 Distribution of the nearest bamboo number of live bamboo,new bamboo and dead bamboo

3.2 角尺度與毛竹更替的關(guān)系

角尺度（Wi）可用于描述林木空間分布格局。根據(jù)Voronoi角尺度判別標(biāo)準(zhǔn)［26］，林木空間分布格局有3種：Wi∈［0.327，0.357］，隨機(jī)分布；Wi＜0.327，均勻分布；Wi＞0.357，聚集分布。對(duì)各大年活竹、新竹、死竹角尺度進(jìn)行統(tǒng)計(jì)（表3）發(fā)現(xiàn)，新竹呈隨機(jī)分布的頻率高于活竹，原因是毛竹受所處生境條件影響，新竹呈現(xiàn)隨機(jī)分布的可能性變大，與王微等［27］對(duì)重慶典型毛竹分株種群結(jié)構(gòu)特征研究結(jié)果一致。死竹呈聚集分布的頻率高于活竹，原因是聚集程度較高的區(qū)域競(jìng)爭(zhēng)也更激烈，毛竹更易死亡。整體來看，各大年毛竹活竹、新竹和死竹平均角尺度均大于0.357，說明毛竹林整體上呈聚集分布，與黃麗霞等［28］研究結(jié)果一致。平均角尺度的單方差分析結(jié)果則表明，三者不存在顯著性差異（P＞0.05）；說明毛竹林整體的空間分布格局少受毛竹更替的影響，具有相對(duì)穩(wěn)定性。

表3 各大年毛竹活竹、新竹和死竹角尺度統(tǒng)計(jì)Table 3 Statistics of uniform angle of live bamboo,new bamboo and dead bamboo in 2010,2012 and 2014

3.3 大小比數(shù)與毛竹更替的關(guān)系

大小比數(shù)描述相鄰毛竹之間的差異，反映立竹在毛竹林中的生長(zhǎng)優(yōu)勢(shì)程度。對(duì)某一對(duì)象竹而言，大小比數(shù)取值為 0.00，（0.00， 0.25］，（0.25， 0.50］，（0.50， 0.75］，（0.75， 1.00］時(shí)， 分別對(duì)應(yīng)優(yōu)勢(shì)、 亞優(yōu)勢(shì)、中庸、劣勢(shì)和絕對(duì)劣勢(shì)。對(duì)各大年活竹、新竹、死竹的大小比數(shù)進(jìn)行統(tǒng)計(jì)（表4）可以看出：新竹處于優(yōu)勢(shì)和亞優(yōu)勢(shì)的頻率都大于活竹，說明新竹在林分生長(zhǎng)中處于有利地位；而死竹處于劣勢(shì)和絕對(duì)劣勢(shì)的頻率高于活竹，說明在受壓環(huán)境下的毛竹更易死亡?？偟膩砜?，毛竹林整體處于中庸?fàn)顟B(tài)，大小比數(shù)的頻率分布基本呈現(xiàn)均衡分布的特征；新竹主要處于亞優(yōu)勢(shì)狀態(tài)，而死竹主要處于劣勢(shì)狀態(tài)。單因素方差分析結(jié)果表明：各大年毛竹活竹、新竹和死竹在大小比數(shù)上存在極顯著差異（P＜0.01），大小比數(shù)排序?yàn)樗乐瘢净钪瘢拘轮瘢从沉肆⒅竦膬?yōu)勢(shì)程度對(duì)毛竹更替有重要影響。因此，在毛竹林經(jīng)營(yíng)過程中，應(yīng)提早采伐處于受壓狀態(tài)的毛竹，釋放生長(zhǎng)空間，為潛在新竹創(chuàng)造有利生長(zhǎng)條件。

表4 各大年毛竹活竹、新竹和死竹大小比數(shù)統(tǒng)計(jì)Table 4 Statistics of neighborhood comparison of live bamboo,new bamboo and dead bamboo in 2010,2012 and 2014

3.4 競(jìng)爭(zhēng)強(qiáng)度與毛竹更替的關(guān)系

計(jì)算各大年活竹、新竹、死竹競(jìng)爭(zhēng)指數(shù)，并根據(jù)毛竹競(jìng)爭(zhēng)指數(shù)（Ci）的大小，將競(jìng)爭(zhēng)強(qiáng)度劃分成3個(gè)等級(jí)：Ci≤4，弱度競(jìng)爭(zhēng)；4＜Ci≤8，中度競(jìng)爭(zhēng)；8＜Ci，強(qiáng)度競(jìng)爭(zhēng)。分析結(jié)果見表5。從表5可知：新竹處于弱度競(jìng)爭(zhēng)的頻率高于活竹，處于中度競(jìng)爭(zhēng)和強(qiáng)度競(jìng)爭(zhēng)的頻率低于活竹，平均競(jìng)爭(zhēng)指數(shù)小于活竹，說明新竹中有較多個(gè)體在弱競(jìng)爭(zhēng)環(huán)境下生長(zhǎng)。而死竹處于強(qiáng)度競(jìng)爭(zhēng)的頻率大于活竹，平均競(jìng)爭(zhēng)指數(shù)也大于活竹，說明在強(qiáng)競(jìng)爭(zhēng)環(huán)境下的毛竹更易死亡。對(duì)各大年毛竹活竹、新竹和死竹的競(jìng)爭(zhēng)指數(shù)進(jìn)行單因素方差分析，結(jié)果表明：活竹和新竹無顯著差異（P＞0.05），而活竹和死竹、新竹和死竹間存在顯著差異（P＜0.05）。推測(cè)原因是新竹在筍期時(shí)已經(jīng)歷競(jìng)爭(zhēng)，保留了較多優(yōu)勢(shì)新竹，因此，在整體呈聚集分布的毛竹林中處于有利競(jìng)爭(zhēng)地位；同時(shí)受環(huán)境條件影響，部分新竹在弱競(jìng)爭(zhēng)環(huán)境下生長(zhǎng)，逐漸加劇這個(gè)區(qū)域的競(jìng)爭(zhēng)強(qiáng)度，表現(xiàn)在競(jìng)爭(zhēng)強(qiáng)度上為無顯著差異；而死竹在強(qiáng)競(jìng)爭(zhēng)環(huán)境下，多處于劣勢(shì)和絕對(duì)劣勢(shì)易被淘汰，表現(xiàn)在競(jìng)爭(zhēng)指數(shù)上為存在顯著差異，說明競(jìng)爭(zhēng)是驅(qū)動(dòng)毛竹更替的重要因素。

表5 各大年毛竹活竹、新竹和死竹競(jìng)爭(zhēng)強(qiáng)度統(tǒng)計(jì)Table 5 Statistics of competition intensity of livel bamboo,new bamboo and dead bamboo in 2010,2012 and 2014

3.5 影響毛竹林更替的重要空間結(jié)構(gòu)因子分析

為確定影響毛竹林更替的主要空間結(jié)構(gòu)因子，對(duì)各大年毛竹林新竹和死竹的空間結(jié)構(gòu)指數(shù)進(jìn)行主成分分析。由表6可知：新竹空間結(jié)構(gòu)指數(shù)的第1主成分中最近鄰竹株數(shù)的因子載荷最大，表明毛竹新竹受最近鄰竹的影響最大；最近鄰竹株數(shù)和角尺度為正載荷，大小比數(shù)和競(jìng)爭(zhēng)指數(shù)為負(fù)載荷，表明前兩者和后兩者對(duì)毛竹新竹的影響不同，前兩者反映毛竹的分布特征，后兩者反映毛竹之間的競(jìng)爭(zhēng)關(guān)系。第2主成分中，競(jìng)爭(zhēng)指數(shù)的載荷最大，反映立竹在林分中的競(jìng)爭(zhēng)強(qiáng)度對(duì)毛竹新竹存在重要影響。第3主成分中主要反映了立竹在林分中的生長(zhǎng)優(yōu)勢(shì)程度對(duì)新竹的影響。主成分分析表明：影響新竹生長(zhǎng)的空間結(jié)構(gòu)指數(shù)重要性排序依次為最近鄰竹株數(shù)＞大小比數(shù)＞競(jìng)爭(zhēng)指數(shù)＞角尺度。

分析表7發(fā)現(xiàn)：第1主成分中最近鄰竹株數(shù)的因子載荷最大，表明毛竹死亡受最近鄰竹的影響最大；最近鄰竹株數(shù)和角尺度為正載荷，大小比數(shù)和競(jìng)爭(zhēng)指數(shù)為負(fù)載荷，表明前兩者和后兩者對(duì)毛竹死亡的影響不同。第2主成分中，大小比數(shù)的載荷最大，反映立竹在林分中的生長(zhǎng)優(yōu)勢(shì)程度對(duì)毛竹死亡存在重要影響。第3主成分中主要反映了競(jìng)爭(zhēng)強(qiáng)度對(duì)毛竹死亡的影響。主成分分析結(jié)果表明：影響毛竹死亡的空間結(jié)構(gòu)指數(shù)重要性排序依次為最近鄰竹株數(shù)＞大小比數(shù)＞競(jìng)爭(zhēng)指數(shù)＞角尺度。

4 結(jié)論與討論

本研究以浙江天目山國(guó)家級(jí)自然保護(hù)區(qū)毛竹林為研究對(duì)象，利用最近鄰竹株數(shù)、角尺度、大小比數(shù)和競(jìng)爭(zhēng)指數(shù)等4個(gè)空間結(jié)構(gòu)指數(shù)，研究了空間結(jié)構(gòu)與毛竹林更替的關(guān)系，為提高毛竹林經(jīng)營(yíng)管理水平提供了參考依據(jù)。研究表明：各大年活竹、新竹和死竹的最近鄰竹株數(shù)均為6株，三者在角尺度上無顯著性差異，均具有相對(duì)穩(wěn)定的聚集分布特征。說明雖然毛竹屬散生竹，但其繁殖方式是以母竹為中心向周邊發(fā)筍長(zhǎng)出新竹，新竹分布受竹鞭延伸限制，所以呈現(xiàn)聚集分布；但聚集分布會(huì)造成毛竹種內(nèi)競(jìng)爭(zhēng)和自然死亡，且這種死亡并不改變其聚集分布特征。顯然，毛竹林分布格局與其更替存在密切相關(guān)。

在毛竹林經(jīng)營(yíng)管理中，應(yīng)當(dāng)充分考慮毛竹林更替與空間結(jié)構(gòu)關(guān)系，通過砍竹和留筍等措施優(yōu)化調(diào)控毛竹林的空間結(jié)構(gòu)［25,29］，為毛竹林豐產(chǎn)營(yíng)造最佳密度、胸徑和年齡結(jié)構(gòu)［30］。本研究結(jié)果提示：為減少毛竹生長(zhǎng)過程中的自然損耗，為毛竹生長(zhǎng)創(chuàng)造有利條件，應(yīng)對(duì)毛竹林高密度區(qū)域進(jìn)行疏伐（最近鄰竹株數(shù)＞7株）；對(duì)處于受壓狀態(tài)的毛竹，還應(yīng)提早采伐；維持毛竹林適當(dāng)?shù)母?jìng)爭(zhēng)強(qiáng)度，可避免因種內(nèi)競(jìng)爭(zhēng)激烈而導(dǎo)致毛竹減產(chǎn)。毛竹林在生長(zhǎng)過程中受多因素的影響，實(shí)際上，環(huán)境因子如地形因子、土壤因子、植被因子等都影響毛竹林的更替，因此，結(jié)合立地條件分析毛竹林空間與更替的關(guān)系有待進(jìn)一步研究。

表6 各大年毛竹新竹空間結(jié)構(gòu)指數(shù)主成分分析Table 6 Principal component analysis on spatial structure indices of new bamboo in 2010,2012 and 2014

表7 各大年毛竹死竹空間結(jié)構(gòu)指數(shù)主成分分析Table 7 Principal component analysis on spatial structure indices of dead bamboo in 2010,2012 and 2014

［1］劉珉.多角度解讀第8次全國(guó)森林資源清查結(jié)果［J］.林業(yè)經(jīng)濟(jì)，2014（5）：3－15.LIU Min.Study on the 8th national forest inventory from multiple perspective ［J］.For Econ,2014（5）:3 － 15.

［2］周國(guó)模，姜培坤.毛竹林的碳密度和碳貯量及其空間分布［J］.林業(yè)科學(xué)，2004，40（6）：20－24.ZHOU Guomo,JIANG Peikun.Density,storage and spatial distribution of carbon in Phyllostachys pubescens forest［J］.Sci Silv Sin,2004,40（6）:20 － 24.

［3］肖復(fù)明，范少輝，汪思龍，等.毛竹、杉木人工林生態(tài)系統(tǒng)碳平衡估算［J］.林業(yè)科學(xué)，2010，46（11）：59－65.XIAO Fuming,FAN Shaohui,WANG Silong,et al.Estimation of the carbon balance in Moso bamboo and Chinese fir plantation ecosystem ［J］.Sci Silv Sin,2010,46（11）:59 － 65.

［4］漆良華，劉廣路，范少輝，等.不同撫育措施對(duì)閩西毛竹林碳密度、碳貯量與碳格局的影響［J］.生態(tài)學(xué)雜志，2009， 28（8）： 1482 － 1488.QI Lianghua,LIU Guanglu,FAN Shaohui,et al.Effects of different tending measures on carbon density,storage,and allocation pattern of Phyllostachy edulis forests in western Fujian Province ［J］.Chin J Ecol,2009,28（8）:1482 － 1488.

［5］PRETZSCH H.Analysis and modeling of spatial stand structures:methodological considerations based on mixed beech-larch stands in Lower Saxony ［J］.For Ecol Manage,1997,97（3）:237 － 253.

［6］湯孟平，唐守正，雷相東，等.林分擇伐空間結(jié)構(gòu)優(yōu)化模型研究［J］.林業(yè)科學(xué)，2004，40（5）：25－31.TANG Mengping,TANG Shouzheng,LEI Xiangdong,et al.Study on spatial structure optimizing model of stand selection cutting ［J］.Sci Silv Sin,2004,40（5）:25 － 31.

［7］AGUIRRE O,HUI Gangying,GADOW K,et al.An analysis of spatial forest structure using neighborhood-based variables ［J］.For Ecol Manage,2003,183（1/3）:137 － 145.

［8］湯孟平.森林空間結(jié)構(gòu)研究現(xiàn)狀與發(fā)展趨勢(shì)［J］.林業(yè)科學(xué)，2010，46（1）：117－122.TANG Mengping.Advances in study of forest spatial structure ［J］.Sci Silv Sin,2010,46（1）:117 － 122.

［9］鄭進(jìn)烜，董文淵，陳沖，等.海子坪天然毛竹種群生長(zhǎng)規(guī)律研究［J］.竹子研究匯刊，2008，27（2）：32－37.ZHENG Jinxuan,DONG Wenyuan,CHEN Chong,et al.Studies on the growth and rhythm of natural moso bamboo population in Haiziping ［J］.J Bamboo Res,2008,27（2）:32 － 37.

［10］賈秀紅，鄭小賢.長(zhǎng)白山過伐林區(qū)云冷杉針闊混交林空間結(jié)構(gòu)分析［J］.華中農(nóng)業(yè)大學(xué)學(xué)報(bào)，2006，25（4）：436－440.JIA Xiuhong,ZHENG Xiaoxian.Analysis on the spatial structure of the spruce-fir mixed stands of coniferous and broad-leaved trees in the over-logged forest［J］.J Huazhong Agric Univ,2006,25（4）:436 － 440.

［11］毛磊，王冬梅，楊曉暉，等.樟子松幼樹在不同林分結(jié)構(gòu)中的空間分布及其更新分析［J］.北京林業(yè)大學(xué)學(xué)報(bào)， 2008， 30（6）： 71 － 77.MAO Lei,WANG Dongmei,YANG Xiaohui,et al.Spatial patterns of young Pinus sylvestris var.mongolica saplings and their regeneration analysis in different stands of Inner Mongolia,northern China ［J］.J Beijing For Univ,2008,30（6）:71 － 77.

［12］席青虎，鐵牛，淑梅，等.寒溫帶興安落葉松林天然更新研究［J］.林業(yè)資源管理，2009（1）：44－48.XI Qinghu,TIE Niu,SHU Mei,et al.Study on natural regeneration of Larix gmelinii in cold-temperate zone ［J］.For Resour Manage,2009（1）:44-48.

［13］湯孟平，周國(guó)模，陳永剛，等.基于Voronoi圖的天目山常綠闊葉林混交度［J］.林業(yè)科學(xué)，2009，45（6）：1－5.TANG Mengping,ZHOU Guomo,CHEN Yonggang,et al.Mingling of evergreen broad-leaved forests in Tianmu Mountain based on Voronoi diagram ［J］.Sci Silv Sin,2009,45（6）:1-5.

［14］仇建習(xí)，湯孟平，沈利芬，等.天目山近自然毛竹林空間結(jié)構(gòu)與胸徑的關(guān)系［J］.生態(tài)學(xué)報(bào)，2015，35（12）：4081－4088.QIU Jianxi,TANG Mengping,SHEN Lifen,et al.The relationship between spatial structure and DBH of close-to-nature Phyllostachys edulis stands in Tianmu Mountain ［J］.Acta Ecol Sin,2015,35（12）:4081 － 4088.

［15］丁麗霞，王祖良，周國(guó)模，等.天目山國(guó)家級(jí)自然保護(hù)區(qū)毛竹林?jǐn)U張遙感監(jiān)測(cè)［J］.浙江林學(xué)院學(xué)報(bào)，2006， 23（3）： 297 － 300.DING Lixia,WANG Zuliang,ZHOU Guomo,et al.Monitoring Phyllostachys pubescens stands expansion in National Nature Reserve of Mount Tianmu by remote sensing ［J］.J Zhejiang For Coll,2006,23（3）:297 － 300.

［16］陳新安.毛竹林大小年生長(zhǎng)規(guī)律探討［J］.中南林業(yè)調(diào)查規(guī)劃，2010，29（1）：21－23.CHEN Xin’an.Discussion on the growth regulation of on-year and off-year moso bamboo ［J］.Cent South For Inv Plan,2010,29（1）:21 － 23.

［17］湯孟平，陳永剛，施擁軍，等.基于Voronoi圖的群落優(yōu)勢(shì)樹種種內(nèi)種間競(jìng)爭(zhēng)［J］.生態(tài)學(xué)報(bào)，2007，27（11）：4707－4716.TANG Mengping,CHEN Yonggang,SHI Yongjun,et al.Intraspecific and interspecific competition analysis of community dominant plant populations based on Voronoi diagram ［J］.Acta Ecol Sin,2007,27（11）:4707 － 4716.

［18］GADOOW K V,HUI Gangying,ALBERT M.Das Winkelmaβ-ein strukturparameter zur beschreibung der individualverteilung in Waldbest?nden ［J］.Centr Ges For,1998,115（1）:1 － 9.

［19］POMMERENING A.Evaluating structural indices by reversing forest structural analysis ［J］.For Ecol Manage,2006,224（3）:266 － 277.

［20］POMMERENING A,STOYAN D.Edge-correction needs in estimating indices of spatial forest structure ［J］.Can J For Res,2006,36（7）:1723 － 1739.

［21］惠剛盈，von GADOW K，ALBERT M.一個(gè)新的林分空間結(jié)構(gòu)參數(shù)：大小比數(shù)［J］.林業(yè)科學(xué)研究，1999，12（1）： 1 － 6.HUI Gangying,von GADOW K,ALBERT M.A new parameter for stand spatial structure:neighborhood comparison［J］.For Res,1999,12（1）:1 － 6.

［22］趙中華，惠剛盈，胡艷波，等.基于大小比數(shù)的林分空間優(yōu)勢(shì)度表達(dá)方法及其應(yīng)用［J］.北京林業(yè)大學(xué)學(xué)報(bào)，2014， 36（1）： 78 － 82.ZHAO Zhonghua,HUI Gangying,HU Yanbo,et al.Method and application of stand spatial advantage degree based on the neighborhood comparison ［J］.J Beijing For Univ,2014,36（1）:78 － 82.

［23］辛營(yíng)營(yíng)，韋新良.青山湖針闊混交林優(yōu)勢(shì)樹種競(jìng)爭(zhēng)的數(shù)量研究［J］.浙江農(nóng)林大學(xué)學(xué)報(bào)，2011，28（4）：601－606.XIN Yingying,WEI Xinliang.Dominance in a mixed conifer and broadleaved forest of Qingshan Lake,Zhejiang ［J］.J Zhejiang A&F Univ,2011,28（4）:601 － 606.

［24］HEGYI F.A simulation model for managing Jack-pine stands ［G］//FRIES J.Growth Models for Tree and Stand Simulation.Stockholm:Royal College of Forestry,1974:74－90.

［25］湯孟平，徐文兵，陳永剛，等.天目山近自然毛竹林空間結(jié)構(gòu)與生物量的關(guān)系［J］.林業(yè)科學(xué)，2011，47（8）：1－6.TANG Mengping,XU Wenbing,CHEN Yonggang,et al.Relationship between spatial structure and biomass of a close-to-nature Phyllostachys edulis stand in Tianmu Mountain ［J］.Sci Silv Sin,2011,47（8）:1 － 6.

［26］李際平，封堯，趙春燕，等.基于Voronoi圖的角尺度分析方法探討：以湖南省平江縣福壽林場(chǎng)為例［J］.林業(yè)資源管理， 2015（4）： 33－38.LI Jiping,FENG Yao,ZHAO Chunyan,et al.Discussion on the analysis method of uniform angle index based on Voronoi diagram:take Fushou Forest Farm as an example ［J］.For Resour Manage,2015（4）:33 － 38.

［27］王微，胡凱，吳冬琴，等.重慶地區(qū)毛竹分株種群結(jié)構(gòu)特征研究［J］.竹子研究匯刊，2008，27（2）：16－20.WANG Wei,HU Kai,WU Dongqin,et al.Studies on the ramet population structure of the clonal moso bamboo in Chongqing ［J］.J Bamboo Res,2008,27（2）:16 － 20.

［28］黃麗霞，袁位高，黃建花，等.不同經(jīng)營(yíng)方式下毛竹林的林分空間結(jié)構(gòu)比較研究［J］.浙江林業(yè)科技，2008， 28（3）： 48 － 51.HUANG Lixia,YUAN Weigao,HUANG Jianhua,et al.Comparative study on spatial structure of Phyllostachys heterocycla var.pubescens stand with different management［J］.J Zhejiang For Sci Technol,2008,28（3）:48 － 51.

［29］曹流清，李曉鳳.毛竹大徑材培育技術(shù)研究［J］.竹子研究匯刊，2003，22（4）：34－41.CAO Liuqing,LI Xiaofeng.Research on cultivation technology of big-diameter bamboo ［J］.J Bamboo Res,2003,22（4）:34 － 41.

［30］鄭郁善，洪偉.毛竹林豐產(chǎn)年齡結(jié)構(gòu)模型與應(yīng)用研究［J］.林業(yè)科學(xué)，1998，34（3）：32－39.ZHENG Yushan,HONG Wei.A study on age structure model of bamboo stand and its application ［J］.Sci Silv Sin,1998,34（3）:32 － 39.