陳 婷，施擁軍,2，周國模,2，鄭澤睿，李 羽中

毛竹碳匯林營造初期林分非空間結(jié)構(gòu)年際變化特征

陳 婷1，施擁軍1,2，周國模1,2，鄭澤睿1，李 羽中1

（1.浙江農(nóng)林大學(xué) 浙江省森林生態(tài)系統(tǒng)碳循環(huán)與固碳減排重點(diǎn)實(shí)驗(yàn)室，浙江 臨安 311300；2.浙江農(nóng)林大學(xué) 亞熱帶森林培育國家重點(diǎn)實(shí)驗(yàn)室培育基地，浙江 臨安311300）

在新造毛竹Phyllostachys edulis碳匯林內(nèi)，設(shè)置80 m×100 m固定標(biāo)準(zhǔn)地，采用固定樣地長期連續(xù)觀測的調(diào)查方法，研究新造毛竹 碳匯林在成林前（第1～5年）林分非空間結(jié)構(gòu)動(dòng)態(tài)變化規(guī)律及影響因子。結(jié)果表明：①新造毛竹碳匯林成林初期，立竹度、林分平均胸徑及林分平均竹高是林分結(jié)構(gòu)變化的主導(dǎo)因子，且隨著造林后年限的增加，立竹度、林分平均胸徑及林分平均竹高均在不斷增加。②新造毛竹碳匯林在成林過程中，立竹度（N）隨造林后年限（y）的增加而增加，并滿足異速生長方程：N=17.214y2＋309.21y－217.4,R2=0.981 1。③隨著造林后年限的增加，新出毛竹的平均胸徑逐年增大，新竹平均胸徑（DBH）與造林后年限（y）的異速生長方程為：DBH=0.015 7y2＋0.727 9y＋1.598 0,R2=0.756 5。④每年新出毛竹的平均竹高（H）隨著造林后年限（y）的增加而增高，其擬合方程為：H=0.108 5y2＋0.429 9y＋2.155 3,R2=0.723 2。圖4表8參20

森林測計(jì)學(xué)；新造毛竹林；碳匯林；立竹度；胸徑；竹高；異速生長方程

在全球共同應(yīng)對氣候變化，大力發(fā)展低碳經(jīng)濟(jì)的背景下，森林生態(tài)系統(tǒng)具有強(qiáng)大的固碳減排功能，是陸地上最大的碳儲(chǔ)存庫，在低碳經(jīng)濟(jì)中有著不可替代的作用，森林碳匯成為相當(dāng)重要的研究領(lǐng)域。碳匯造林也成為增加碳儲(chǔ)量的一種重要的方式，被中國乃是世界各國廣泛采用。竹林是重要的森林植被類型，也是特殊的生態(tài)系統(tǒng)單元。竹林生態(tài)系統(tǒng)的碳儲(chǔ)量及固碳能力不可忽視，成為除喬木林以外最大的碳匯，且竹林擁有穩(wěn)定的碳匯功能［1－3］。許多學(xué)者研究表明：毛竹Phyllostachys edulis具有很強(qiáng)的固碳能力，是速生階段杉木的1.46倍，是熱帶山地雨林的1.33倍，是蘇南27年生杉木林的2.16倍［4］。同時(shí)毛竹又具有良好的經(jīng)濟(jì)利用價(jià)值，深受山區(qū)農(nóng)民喜愛，因此，毛竹造林成為中國碳匯造林的一種重要方式。碳匯造林是指在確定了基線的土地上，以增加碳匯為主要目的，對造林及林木（分）生長過程實(shí)施碳匯計(jì)量和監(jiān)測而開展的有特殊要求的造林活動(dòng)。與普通的造林相比，碳匯造林突出森林的碳匯功能，強(qiáng)調(diào)森林的多重效益，需要對碳儲(chǔ)量的變化過程進(jìn)行科學(xué)的計(jì)量與監(jiān)測工作，以獲得可測量、可報(bào)告、可核查的碳匯量，以滿足后續(xù)碳匯交易的要求。新造毛竹林生長特殊，分為成林前和成林后2個(gè)階段，成林之前的竹林林分生長較為復(fù)雜，隨著造林后年限的增長，毛竹林林分的立竹度、平均胸徑、平均竹高也會(huì)不斷發(fā)生變化，這與成林后較為穩(wěn)定的林分結(jié)構(gòu)（立竹度、胸徑、竹高）有著本質(zhì)的區(qū)別。雖然，對毛竹林成林后的林分結(jié)構(gòu)已經(jīng)有了較多的研究［5－18］，但是對于毛竹林成林前的研究幾乎為空白。竹林成林前的林分結(jié)構(gòu)和碳儲(chǔ)量變化是當(dāng)前竹林碳匯計(jì)量中面臨的最大困難。為此，本研究擬通過碳匯造林試驗(yàn)和長期的監(jiān)測工作，以探究毛竹成林前的林分結(jié)構(gòu)動(dòng)態(tài)變化特征，建立新造毛竹林成林前的異速生長模型，為毛竹林的碳匯計(jì)量監(jiān)測提供支持。

1 研究區(qū)域概況

浙江農(nóng)林大學(xué)2008年在浙江省臨安市藻溪鎮(zhèn)嚴(yán)家山，通過母竹移栽方式新造了1塊毛竹碳匯林，面積為47.72 hm2。研究試驗(yàn)區(qū)選在新造的毛竹碳匯林內(nèi)，中心地理坐標(biāo)為30°10′20″N,119°29′40″E，平均海拔為550 m。該區(qū)屬中亞熱帶季風(fēng)氣候，年平均氣溫為15.8℃，年降水量為1 500 mm左右，無霜期為236 d，地形地貌為低山丘陵，土類為紅壤，亞類為黃紅壤，土屬為黃泥土，土種為黃泥砂土。成土母質(zhì)主要為砂頁巖類殘坡積物和砂礫類殘坡積物，土壤呈酸性。土壤pH 4.4，土壤有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)為36 g·kg－1，堿解氮254.3 mg·kg－1，速效磷126.3 mg·kg－1，有效磷為17.6 mg·kg－1。

毛竹碳匯林初植密度為450株·hm－2，在毛竹林發(fā)育成林階段，竹林撫育措施為春季和初夏各中耕除草1次·a－1,造林后3 a內(nèi)采用穴施方式施用尿素300～450 kg·hm－2·a－1，除挖取退筍、敗筍外,全部留筍長竹,不進(jìn)行伐竹作業(yè)。

2 野外調(diào)查方法

2.1 樣地設(shè)置

2008年，在浙江省臨安市藻溪鎮(zhèn)嚴(yán)家山新造的毛竹碳匯林地內(nèi)，選擇典型的毛竹林設(shè)置80 m× 100 m固定標(biāo)準(zhǔn)地，固定標(biāo)準(zhǔn)地海拔為450～500 m。用高精度全站儀（NTS-352，測角精度2″，測距精度3 mm＋2 mm/km）把固定標(biāo)準(zhǔn)地劃分為完全等面積的20個(gè)20 m×20 m的樣地（調(diào)查單元），使樣地構(gòu)成相鄰格網(wǎng)。

2.2 坐標(biāo)測定

采用全站儀（NTS-352）精確測定每個(gè)角樁的三維坐標(biāo)（x,y,z）及每株毛竹基部三維坐標(biāo)，其中x和y是平面直角坐標(biāo)，z是海拔高度。

2.3 毛竹林林分因子調(diào)查

在每個(gè)調(diào)查單元內(nèi)，對胸徑DBH≥1 cm的毛竹進(jìn)行每木調(diào)查。從2008年開始，每隔2 a測量新竹基部三維坐標(biāo)，并進(jìn)行每木檢尺，測量并記錄竹編號(hào)、竹齡、竹高、胸徑等因子。對樣地內(nèi)新出毛竹進(jìn)行連續(xù)不間斷監(jiān)測。胸徑的測量采用圍尺，以距地面1.3 m處胸高位置為測定標(biāo)準(zhǔn)，胸高處出現(xiàn)節(jié)疤，凹凸或其他不正常情況時(shí)，可在胸高上下距離相等而形狀較正常處，測2個(gè)直徑取平均數(shù)作為胸徑；竹高的測量采用全站儀懸高測量法（測定空中某點(diǎn)距地面的垂直高度）。到2013年（為期5 a）共測定毛竹1 809株及母竹330株。由于母竹在造林后幾年內(nèi)被清除，所以本研究所述的立竹度為不包括母竹的新竹立竹度。

3 結(jié)果與分析

3.1 新造毛竹碳匯林林分結(jié)構(gòu)年際變化特征

竹齡、立竹度、胸徑和竹高等是林分結(jié)構(gòu)變化的重要指標(biāo)。新造毛竹林主要生長因子動(dòng)態(tài)變化見表1。從表1可以看出：在研究樣地內(nèi)，新造毛竹林造林后第1年成竹數(shù)為48株，到造林后第5年立竹度達(dá)到1 809株，呈急速上升趨勢；新造毛竹林內(nèi)毛竹林林分平均胸徑從第1年的2.6 cm增加到第5年的4.4 cm；新造毛竹林林分平均竹高也發(fā)生著明顯的變化，造林后第1年林分平均竹高3.1 m，到造林后第5年林分平均竹高達(dá)到5.3 m。

綜上可知：新造毛竹林在整個(gè)成林過程中，隨著造林后年限的增加，立竹度在不斷增加，林分平均胸徑也在增大，林分平均竹高在不斷增高?？梢钥闯觯盒略烀窳质且环N隨著造林年限不斷動(dòng)態(tài)變化著的，立竹度、林分平均胸徑以及林分平均竹高與造林后年限之間存在異速生長模型。

表1 新造毛竹碳匯林林分特征Table 1 The characteristics of moso bamboo samples

3.2 新造毛竹林各年齡非空間結(jié)構(gòu)指標(biāo)的方差分析

為探究新造毛竹林非空間結(jié)構(gòu)的各指標(biāo)年際變化的影響及其顯著性水平，采用SPSS 19.0軟件，對造林后不同年限的毛竹的立竹度、新竹胸徑、新竹竹高等做單因素方差分析，結(jié)果見表2?？梢钥闯觯盒略烀窳种胁煌炝趾竽晗蓿ㄖ颀g）新竹的立竹度、胸徑、竹高等因子均呈極顯著差異（P＜0.01），說明造林后年限的變化是新造毛竹林成林前的林分結(jié)構(gòu)變化的重要因子。

為探究大樣地內(nèi)部是否存在差異性，對20個(gè)子樣地做方差分析，結(jié)果見表3。表3表明：樣地之間立竹度具有顯著差異（P=0.000），但不管是對胸徑還是對竹高而言，樣地之間的P都遠(yuǎn)大于0.05，說明樣地之間胸徑和竹高沒有顯著差異；

表2 新造毛竹碳匯林不同年齡非空間結(jié)構(gòu)各指標(biāo)變化的顯著性檢驗(yàn)Table 2 Each index’s significant test of non spatial structurewith different ages

表3 新造毛竹碳匯林各子樣地間非空間結(jié)構(gòu)各指標(biāo)變化的顯著性檢驗(yàn)Table3 Each index’s significant test of non spatial structurewith different sub-plots

3.3 新造毛竹碳匯林不同造林后年限毛竹立竹度變化

立竹度是竹林出筍前單位面積內(nèi)正常生長的竹子株數(shù)，是新造毛竹林林分結(jié)構(gòu)變化的重要指標(biāo)，在毛竹林成林前期，會(huì)隨著造林年限呈急速增加趨勢。表4為新造毛竹林新竹株數(shù)的年際變化?？梢姡谘芯繕拥貎?nèi)，新造毛竹林造林后第1年新竹的成竹數(shù)48株，造林后第2年新竹的成竹數(shù)為555株，為出筍成竹盛期，第3年和第4年新竹分別成竹229株和374株，第5年新竹出筍成竹數(shù)為603株。第1年新竹出筍成竹數(shù)較少，第2年和第5年新竹出筍成竹數(shù)較高，到達(dá)出筍成竹盛期。但林分內(nèi)，隨著造林年限的增加，立竹度隨造林后年限的增加呈線性上升趨勢（表1）。

采用Pearson相關(guān)性分析方法對新造毛竹林立竹度與造林后年限之間做相關(guān)性分析，呈顯著相關(guān)性（R=0.991）。對不同造林后年限和立竹度進(jìn)行方差分析（表2），F(xiàn)分布的顯著性概率為0.000，達(dá)到了極顯著水平，有統(tǒng)計(jì)學(xué)意義。這說明因變量和自變量的關(guān)系極顯著。在成林過程中，立竹度（N）隨造林后年限（y）變化的異速生長方程為：N=17.214y2+309.210y－217.400（R2=0.981 1）。據(jù)此模型，可以通過了解新造毛竹林的造林后年限計(jì)算出新造毛竹林各林分的立竹度。

表4 新造毛竹碳匯林新竹株數(shù)年際變化Table 4 Changes of bamboo stand culms with year

3.4 新造毛竹碳匯林各年齡毛竹胸徑變化

新造毛竹林在成林過程中，隨著造林后年限的增加，除了立竹度增加以外，新竹平均胸徑也會(huì)逐步增大，與成林后的穩(wěn)定竹林表現(xiàn)出完全不一致的特征［19］。本研究表明（表5）：新造毛竹林在造林后第1年新竹平均胸徑為2.6 cm，到成林后的第5年新竹平均胸徑增加到5.8 cm，說明隨著年齡的增大新造毛竹林新竹平均胸徑在快速的增加，各年份最小胸徑差異性不大，為1.0～1.4 cm；胸徑的極大值變化較大，變化幅度為5.4～10.2 cm。隨著造林年限的增加，毛竹林個(gè)體之間的差異也在不斷變大，標(biāo)準(zhǔn)差從第1年的0.88 cm增加到第5年的1.56 cm。

表5 新造毛竹碳匯林新竹胸徑年際變化Table 5 Change of new bamboo stand diameter in breast height with year

采用Pearson相關(guān)性分析方法對新造毛竹林新竹胸徑與造林后年限之間做相關(guān)性分析，呈顯著相關(guān)性（R=0.992）。對新造毛竹林各子樣地內(nèi)不同年齡新竹胸徑求平均，得到散點(diǎn)圖（圖1），并進(jìn)行曲線擬合，得到不同新造毛竹林新竹胸徑年際變化異速生長模型（表6），結(jié)果表明：新造毛竹林在成林過程中，新竹的平均胸徑逐年增大。通過模型比較，多項(xiàng)式模型DBH=0.015 7y2＋0.727 9y＋1.598 0,R2=0.756 5，能夠最好地反映新造毛竹林新竹平均胸徑隨年齡的遞增狀況。其次是線性模型D=0.824 8y＋1.477 9,R2= 0.756 1。

從圖2可以看出：新造毛竹林立竹度按徑階的分布符合威布爾分布，株數(shù)最多的集中在1.5～6.5 cm，隨著造林后年限的增加，新竹的徑階分布將向大徑階移動(dòng)，未來將逐漸趨向于正態(tài)分布。

同時(shí)，對研究樣地內(nèi)對不同造林年限的毛竹林立竹度和林分平均胸徑進(jìn)行相關(guān)性分析，發(fā)現(xiàn)新造毛竹林林分平均胸徑和立竹度（林分密度）之間呈正相關(guān)（R2=0.974 0，P=0.005），對新造毛竹林內(nèi)各年份林分平均胸徑和立竹度進(jìn)行曲線擬合，得到線性方程： DBH=0.000 9N＋2.366 9，R2=0.950 0，函數(shù)曲線見圖3。此模型說明:新造毛竹林在成林初期，隨著造林后年限的增加，新竹平均胸徑在年際增大的同時(shí)，立竹度也在不斷增加。竹林立竹度的大小，與竹林新竹產(chǎn)量有密切關(guān)系。由此看來，應(yīng)該是立竹度影響翌年的新竹大小。相對于毛竹胸徑的測量，立竹度數(shù)據(jù)更便于獲得。

3.5 新造毛竹碳匯林各年齡毛竹竹高變化

新造毛竹林在成林過程中，隨著造林后年限的增加，除了立竹度和新竹平均胸徑在不斷增加，新竹平均竹高也會(huì)逐漸變高，與成林后的穩(wěn)定竹林表現(xiàn)出完全不一致的特征［19］。本研究表明（表7），新造毛竹林在造林后第1年平均新竹竹高為3.1 m，到造林后的第5年新竹平均竹高迅速增加到7.3 m，說明隨著造林后年限的增加，新造毛竹林新竹平均竹高也在快速的高，各年份新竹竹高最小值的差異性不大，為1.0～1.9 m；新竹竹高的變化主要體現(xiàn)在極大值上，變化幅度為6.3～18.0 m。隨著造林年限的增加，毛竹林個(gè)體之間的差異也在不斷變大，標(biāo)準(zhǔn)差從第1年的1.01 m增加到第5年的2.06 m。

表6 新造毛竹碳匯林新竹胸徑年際變化異速生長模型比較Table 6 Different types of allometric models for predicting new moso bamboo by DBHand age

圖1 新造毛竹碳匯林新竹平均胸徑年際變化Figure 1 Change of bamboo stand diameterin bright with year

圖2 新造毛竹碳匯林徑階分布Figure 2 Diameter distribution of mosobamboo

圖3 不同年限立竹度和林分平均胸徑之間的關(guān)系Figure 3 Relationships between mean DBHand culms of new moso bamboo

表7 新造毛竹碳匯林新竹竹高年際變化Table 7 Changes of the height of new moso bamboo stand with year

采用Pearson相關(guān)性分析方法對新造毛竹林新竹竹高和造林后年限之間做相關(guān)性分析，結(jié)果表明:新造毛竹林新竹竹高與造林后年限呈顯著相關(guān)（R2=0.884 0）。對20個(gè)子樣地中的不同造林后年限的新竹竹高分別取平均，生成散點(diǎn)圖（圖4），并進(jìn)行曲線擬合，得到不同新造毛竹林新竹平均竹高年際變化異速生長模型（表8）。結(jié)果表明：新造毛竹林在成林過程中，新竹平均竹高和造林后年限呈正相關(guān)，隨著新造毛竹林造林后年限的增加，新竹的平均竹高也逐年增大。通過模型比較發(fā)現(xiàn)，多項(xiàng)式模型：H= 0.108 5y2＋0.429 9y＋2.155 3,R2=0.723 2能夠最好地反映新造毛竹林新竹平均竹高隨不同造林后年限的遞增狀況；其次是線性模型H=1.099 8y＋1.324 4,R2=0.713 5。

4 結(jié)論與討論

圖4 新造毛竹碳匯林子樣地新竹平均竹高年際變化Figure 4 Changes of bamboo stand culm height with year in sub-plots

本研究在新造毛竹碳匯林地內(nèi)設(shè)置典型固定大樣地，并連續(xù)監(jiān)測，共獲得5 a新造毛竹碳匯林林分?jǐn)?shù)據(jù)，對毛竹林非空間結(jié)構(gòu)特征及動(dòng)態(tài)變化進(jìn)行研究，并對數(shù)據(jù)建模。得到如下結(jié)論。①新造毛竹碳匯林立竹度、造林后年限和林分平均胸徑是新造毛竹碳匯林林分非空間結(jié)構(gòu)變化的主導(dǎo)因子，林分平均竹高是新造毛竹碳匯林林分非空間結(jié)構(gòu)變化的從屬因子；②新造毛竹碳匯林在成林過程中，立竹度（N）隨造林后年限（y）增加呈線性上升，并滿足異速生長方程：N=17.214 0y2＋309.210 0y－217.400 0,R2=0.981 1。③新造毛竹碳匯林在成林過程中，隨著造林年限的增加，新竹的平均胸徑逐年增大，新竹平均胸徑（DBH）與造林后年限（y）滿足異速生長方程為：DBH=0.015 7y2＋0.727 9y＋1.598 0,R2=0.756 5。④新造毛竹碳匯林在成林過程中，新竹平均竹高（H）隨著造林后年限（y）增加逐漸增高，其擬合方程為：H=0.108 5y2＋0.429 9y ＋2.155 3,R2=0.723 2。

本研究發(fā)現(xiàn)新造毛竹碳匯林成林初期林分平均胸徑和立竹度呈正比，其線性關(guān)系為： DBH=0.000 9N＋2.366 9，R2=0.950 0，而毛竹碳匯林成林之后立竹度與平均胸徑呈反比［20］，說明毛竹碳匯林成林前的林分結(jié)構(gòu)變動(dòng)特征與成林后的林分結(jié)構(gòu)特征不一致，其可能的原因是：毛竹林成林前競爭較弱，新竹平均胸徑每年不斷增加，并且由于竹鞭的蔓延，大量發(fā)筍，導(dǎo)致立竹度也不斷增加；成林后由于空間和土壤養(yǎng)分的限制，立竹度的大量增加會(huì)導(dǎo)致林分內(nèi)的競爭漸趨激烈，從而影響發(fā)筍的數(shù)量和質(zhì)量，因此隨著立竹度的增加平均胸徑會(huì)減小。

陳雙林等［19］年對新造毛竹林林分結(jié)構(gòu)演替規(guī)律的研究中提出過新造毛竹林的胸徑異速生長方程，但采用的調(diào)查方法是在同一時(shí)間內(nèi)對不同竹齡的進(jìn)行調(diào)查研究，可能存在以時(shí)間換空間的概念，對竹子的年齡的判斷具有不確定性，而且樣地的初植密度不一，在建模時(shí)很容易產(chǎn)生誤差。而本研究是基于相同初植密度的20塊標(biāo)準(zhǔn)子樣地每年連續(xù)調(diào)查的基礎(chǔ)上進(jìn)行建模，數(shù)據(jù)更加可靠。

本研究是基于新造毛竹碳匯林成林初期（第1～5年）所進(jìn)行的研究，毛竹種植后至成林前，采用一般經(jīng)營模式，通過不斷發(fā)筍、留筍和養(yǎng)竹第7年后可以成林，一般毛竹到達(dá)成林穩(wěn)定階段往往需要6～7 a［19］或更長時(shí)間。本研究今后還將持續(xù)監(jiān)測成林?jǐn)U展特征，補(bǔ)充完善相關(guān)的模型。

表8 新造毛竹碳匯林新竹竹高年際變化異速生長模型比較Table 8 Different types of allometric models for predicting moso bamboo by culm height and age

［1］ 郭起榮，楊光耀，杜天真，等．中國竹林的碳素特征［J］．世界竹藤通訊，2005，3（3）:25－28. GUO Qirong,YANG Gangyao,DU Tianzhen,et al.Carbon character of Chinese bamboo forest［J］.World Bamboo Rattan,2005,3（3）:25－28.

［2］ 肖復(fù)明,范少輝,汪思龍，等.毛竹Phyllostachy pubescens杉木Cunninghamia lanceolata人工林生態(tài)系統(tǒng)碳儲(chǔ)量及其分配特征［J］.生態(tài)學(xué)報(bào),2007,27（7）：2794－2801. XIAO Fuming,FAN Shaohui,WANG Silong,et al.Carbon storage and spatial distribution in Phyllostachy pubescens and Cunninghamia lanceolata plantation ecosystem［J］.Acta Ecol Sin,2007，27（7）：2794－2801.

［3］ 陳先剛,張一平,張小全，等．過去50年中國竹林碳儲(chǔ)量變化［J］.生態(tài)學(xué)報(bào),2008,28（11）：5218－5226．CHEN Xiangang,ZHANG Yiping,ZHANG Xiaoquan,et al.Carbon stock changes in bamboo stands in China over the last 50 years［J］.Acta Ecol Sin,2008,28（11）：5218－5226．

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

［5］ 聶道平.江西省大崗山毛竹林的結(jié)構(gòu)特征［J］.林業(yè)科學(xué)研究，1992，5（6）:693－699. NIE Daoping.The structure characteristics of baoboo groves in Dagangshan Region,Jiangxi Province［J］.For Res, 1992，5（6）:693－699.

［6］ 聶道平.毛竹林結(jié)構(gòu)的動(dòng)態(tài)特性［J］．林業(yè)科學(xué)，1994，30（3）：201－208．NIE Daoping.Structural dynamics of bamboo forest stands［J］.Sci Silv Sin，1994，30（3）:201－207.

［7］ 何東進(jìn)，洪偉，吳承禎.毛竹林林分平均胸徑模擬預(yù)測模型的研究［J］.林業(yè)科學(xué)，2000，36（1）:148－153. HE Dongjin,HONG Wei,WU Chengzhen.A study on simulating predictive model of mean DBH for bamboo stands［J］.Sci Silv Sin,2000,36（1）:148－153.

［8］ 朱錦懋，黃茂提，陳由強(qiáng)，等.筍材兩用毛竹林林分結(jié)構(gòu)數(shù)量關(guān)系研究［J］.植物生態(tài)學(xué)報(bào),2000，24（4）：483 －488. ZHU Jinmao,HUANG Maoti,CHEN Youqiang,at al.The structure of a culm and shoot producing stand of Phyllostachys pubescens［J］.Acta Phytoecol Sin,2000,24（4）：483－488.

［9］ 湯孟平，徐文兵，陳永剛，等.地形條件對毛竹林分結(jié)構(gòu)和植被碳儲(chǔ)量的影響［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．

［10］ 鄧英英，湯孟平，徐文兵，等.天目山近自然毛竹純林的竹稈空間結(jié)構(gòu)特征［J］.浙江農(nóng)林大學(xué)學(xué)報(bào)，2011，28（2）:173－179. DENG Yingying，TANG Mengping，XU Wenbing,et al.Spatial structure of bamboo culm of an almost natural，pure Phyllostachys pubescens forest in Mount Tianmu［J］.J Zhejiang A＆F Univ，2011，28（2）：173－179.

［11］ 白尚斌，周國模，王懿祥,等.天目山國家級(jí)自然保護(hù)區(qū)毛竹擴(kuò)散過程的林分結(jié)構(gòu)變化研究［J］.西部林業(yè)科學(xué)，2012，41（1）:77－82. BAI Shangbin,ZHOU Guomo,WANG Yixiang,et al.Stand structure change of Phyllostachys pubescens forest expansion in Tianmushan National Nature Reserve［J］.J West China For Sci，2012，41（1）:77－82.

［12］ 劉恩斌,施擁軍,李永夫,等.浙江毛竹林分非空間結(jié)構(gòu)特征及其動(dòng)態(tài)變化［J］.林業(yè)科學(xué),2013,49（9）:1－7. LIU Enbin,SHI Yongjun,LI Yongfu,et al.Non spatial structural characteristic of moso bamboo forest and its dynamics in Zhejiang Province［J］.Sci Silv Sin，2013，49（9）:1－7.

［13］ 徐振國,郭起榮,馮云,等.淡竹林分結(jié)構(gòu)與林分因子相關(guān)研究［J］.中南林業(yè)科技大學(xué)學(xué)報(bào),2014,34（2）:72－76. XU Zhenguo,GUO Qirong,FENG Yun,et al.Related research on structure characteristics and stand factors of Phyllostachys glauca［J］.J Cent South Univ For&Technol,2014,34（2）:72－76.

［14］ 范葉青,周國模,施擁軍,等.地形條件對毛竹林分結(jié)構(gòu)和植被碳儲(chǔ)量的影響［J］.林業(yè)科學(xué),2013,49（11）:178 －182. FAN Yeqing,ZHOU Guomo,SHI Yongjun,et al.Effects of terrain on stand structure and vegetation carbon storage of Phyllostachys edulis forest［J］.Sci Silv Sin，2013，49（11）:178－182.

［15］ ARUN J N，GITASREE D，ASHESH K D.Above ground standing biomass and carbon storage in village bamboos in North East India［J］.Biom Bioenergy，2009，33（9）：1188－1196.

［16］ DU Huaqiang，CUI Ruirui，ZHOU Guomo，et al.The responses of moso bamboo （Phyllostachys heterocycla var. pubescens）forest aboveground biomass to Landsat TM spectral reflectance and NDVI［J］.Acta Ecol Sin，2010，30 （5）：257－263.

［17］ YEN T M，JI Y J，LEE J S.Estimating biomass production and carbon storage for a fast-growing makino bamboo （Phyllostachys makinoi）plant based on the diameter distribution model［J］.For Ecol Manage，2010，260（3）:339 －344.

［18］ FU Weijun,JIANG Peikun，ZHAO Keli,et al.The carbon storage in moso bamboo plantation and its spatial variation in Anji County of southeastern China［J］.J Soils Sediment，2014，14（2）:320－329.

［19］ 陳雙林,吳柏林,吳明,等.新造毛竹林林分結(jié)構(gòu)年際演替規(guī)律及影響因子［J］.浙江林學(xué)院學(xué)報(bào),2004,21（4）: 393－397. CHEN Shuanglin,WU Bolin,WU Ming,et al.A study of the interannual succession rule and influential factors of young stands structures of Phyllostachys pubescens［J］.J Zhejiang For Coll,2004,21（4）:393－397.

［20］ YEN T M,LEE J S.Comparing aboveground carbon sequestration between moso bamboo（Phyllostachys heterocycla）and China fir（Cunninghamia lanceolata）forests based on the allometric model［J］.For Ecol Manage,2011,261（6）: 995－1002.

Interannual variation characteristics for stand structures in the early stages of new moso bamboo carbon sink stands

CHEN Ting1,SHI Yongjun1,2,ZHOU Guomo1,2,ZHENG Zerui1,LI Chong1
（1.Zhejiang Provincial Key Laboratory of Carbon Cycling in Forest Ecosystems and Carbon Sequestration,Zhejiang A &F University,Lin’an 311300,Zhejiang,China;2.The Nurturing Station for the State Key Laboratory of Subtropical Silviculture,Zhejiang A&F University,Lin’an 311300,Zhejiang,China）

To determine a rule for interannual variation characteristics in non-spatial,structural stands of new moso bamboo （Phyllostachys edulis）based on carbon sink （1－5 a）.An 80 m×100 m fixed plot was established at a long-term,afforestation,fixed observation research site for carbon sequestration.Linear regression and correlation were used for the analysis.Results showed that 1）stand density,DBH,and culm height were the dominating factors for interannual component changes,and these increased with bamboo age.2）In the young stand,stand density（N）increased with stand year（y）with a fitted regression equation of:N=17.214y2＋309.21y－217.4,R2=0.981 1.3）The allometric equation for mean DBH and stand year（y）fitted with linear regression was DBH=0.015 7y2＋0.727 9y＋1.598 0,R2=0.756 5.4）The linear regression for the relationship between culm height（H）and stand year（y）was projected to be H=0.108 5y2＋0.429 9y＋2.155 3, R2=0.723 2.Therefore,in the young moso bamboo stand,stand density,DBH,and culm height all had a positive relationship with age.［Ch,4 fig.8 tab.20 ref.］

forest mensuration;new moso bamboo forest;carbon sequestration forest;stand density;DBH;culmheight;allometric model

S758

A

2095-0756（2015）02-0181-07

浙 江 農(nóng) 林 大 學(xué) 學(xué) 報(bào)，2015，32（2）：181-187

Journal of Zhejiang A＆F University

10.11833/j.issn.2095-0756.2015.02.003

2014-07-26；

2014-08-28

“十二五”國家科技支撐計(jì)劃項(xiàng)目（2012BAD22B0503）；國家林業(yè)局引進(jìn)國際先進(jìn)農(nóng)業(yè)科學(xué)技術(shù)計(jì)劃（“948”計(jì)劃）項(xiàng)目（2013-4-71）；國家自然科學(xué)基金資助項(xiàng)目（31370637）；浙江省重點(diǎn)科技創(chuàng)新團(tuán)隊(duì)項(xiàng)目（2010R50030）；浙江省自然科學(xué)基金資助項(xiàng)目（Y5110145）

陳婷，從事森林碳匯計(jì)量與監(jiān)測研究。E-mail：ctzun_1990@126.com。通信作者：施擁軍，副教授，從事森林碳匯計(jì)量與監(jiān)測研究。E-mail：syjwwh@163.com