安吉銀坑塢毛竹林長期跟蹤調查分析

周文偉，童曉青，周仁愛，華錫奇，梁一品，俞宏慶，任柏青

（1. 浙江省林業科學研究院，浙江 杭州 310023；2. 浙江省竹類研究重點實驗室，浙江 杭州 310023；3. 浙江省龍泉市城區林業工作站，浙江 龍泉 323700；4. 安吉縣天荒坪林業服務總公司，浙江 安吉 313302）

安吉銀坑塢毛竹林長期跟蹤調查分析

周文偉1,2，童曉青1,2，周仁愛3，華錫奇1,2，梁一品4，俞宏慶4，任柏青4

（1. 浙江省林業科學研究院，浙江 杭州 310023；2. 浙江省竹類研究重點實驗室，浙江 杭州 310023；3. 浙江省龍泉市城區林業工作站，浙江 龍泉 323700；4. 安吉縣天荒坪林業服務總公司，浙江 安吉 313302）

采用1979-2011年（其中缺1989、1991、1993）的抽樣調查數據對安吉縣天荒坪鎮銀坑塢毛竹（Phyllastachys heterocycla var. pubescens）林的密度、胸徑及竹材產量進行了統計分析。結果表明：單位面積每度毛竹株數有較大幅度的增加，從1979年的600株/hm2提高到2011年的1 350株/hm2，其中最高的是2009年的1 440株/hm2；毛竹新竹有增粗的趨勢，其平均胸徑從1979年的8.75 cm逐步增加到2011年的9.08 cm，其中最高的是2009年，達9.47 cm；毛竹新竹材產量也逐年遞增，從1979年的12 945 kg/hm2增加到2011年的25 680 kg/hm2，其中最高的為2009年，達29 790 kg/hm2。

毛竹林；跟蹤調查；統計分析

1 調查區概況

調查區位于我國著名竹鄉安吉縣天荒坪鎮，118° 88″ E，30° 68″ N，屬中亞熱帶海洋性季風氣候區。氣候溫暖濕潤，雨量充沛，年平均氣溫15.6℃，年平均降水量1 378.6 mm，年平均日照1 849 h，無霜期180 d。

安吉縣天荒坪鎮（原安吉縣山河鄉）在1980-1990年，曾經實施了“六五”國家科技攻關項目“安吉縣以竹類資源為主的綜合開發”和“七五”國家科技星火計劃“安吉縣竹類加工為主的配套技術開發—山河鄉中、低產毛竹林資源綜合開發”項目。自1993年起，天荒坪鎮政府為適應林業發展需要，成立了天荒坪林業服務總公司，下設五個林場，經營管理著400余hm2毛竹山。

調查毛竹（Phyllastachys heterocycla var. pubescens）林設在天荒坪鎮的銀坑塢，為毛竹純林，海拔180 ~ 250 m。竹林地土壤為黃紅壤，pH值4.3 ~ 6.0，有機質含量為4.4%，土壤深度40 ~ 100 cm。為了試驗數據的可比性和延續性，只選取了銀坑第一林場和第二林場，總面積120 hm2。所采取的技術措施為：

（1）在“六五”和“七五”期間，主要采取了以下措施：①科學疏筍：出筍期間，除加強護養外，及時組織竹農挖去青頭筍、退筍，以增加竹農收入；②適度鉤梢：留枝15 ~ 18檔，以改進竹冠結構，提高葉面積指數；③合理砍伐：在全林100多戶竹農中，嚴格砍伐制度，控制砍伐量，全林實行留養4度竹，使立竹保持在3 750株/hm2以上，禁砍冷刀；④留養小年竹：改變大小年竹林為均年竹林；⑤全面劈山：在7、8月間全林劈山，齊地劈除雜草灌木，并均勻覆蓋林地；⑥適地松土：在緩坡進行松土，挖去柴蔀頭、竹蔸、死鞭、老鞭，以擴大竹鞭生長發育的有效土壤空間；⑦按需施肥：根據毛竹生長的不同時期，采用溝施、穴施、竹蔸施與立竹稈基根際施相結合的施肥方法和不同種類的肥料及施肥量；⑧竹林保護，防治竹螟等病蟲害。

1993年以來，天荒坪林業服務總公司采取的技術措施為：①全毛竹林從1995年開始，實行“留二砍三”（即所有毛竹林地上的竹只保留4年生，在毛竹生長到第5年時，等當年的新竹成竹后，就開始采伐，俗稱“砍楊梅紅”，并根據市場的需求，分次分批地全部伐去三度的竹子）經營方式；②適地挖山松土，1995-2002年在緩坡毛竹林地先進行了3次挖山松土，而后打了2次草甘磷除草，2003年以后，既不挖山松土，也不打草甘磷；③全面劈山，在全林地進行劈山，先前在7-8月劈山，現考慮作業人員安全和勞動用工缺少的問題，改在冬季進行；④施肥：全毛竹林施毛竹專用復合肥750 kg/hm2；⑤在毛竹林密度較大的地方，對林地上的竹蔸進行一些劈裂和挖除處理。

2 調查方法

在“六五”至“七五”期間，即1980-1990年（其中1979年為非當年新竹調查數據），結合項目實施，根據地形地勢，分別采用完全隨機和隨機區組試驗設計，固定樣地與臨時樣地相結合，除對照樣地面積為20 m ×40 m以外，其余樣地面積均為18.25 m×18.25 m。逐年度調查每一樣地的底竹株數、胸徑、枝下高，并編號登記，建立檔案。每年（度）筍期分別樣地統計出筍支數和筍重。新竹成竹后，每一樣地進行新竹每竹的胸徑、枝下高調查。

1995-2011年，每一度新竹成竹后，采用完全隨機定點方法，抽取一定數量的樣地，樣地面積均為10 m× 40 m，對樣地內竹子分年度進行株數和眉圍的每竹調查。

3 結果與分析

3.1 毛竹新竹株數的變化

從1979年以來毛竹新竹留養株數的動態變化趨勢來看，單位面積新竹株數是遞增的，從1979年的平均600株/hm2，增加到最高的2009年平均1 440株/hm2。其變化大致可以分為4個階段：第一階段為1980年代前，每度新竹 < 900株/hm2；第二階段為1981-1995年，每度新竹915 ~ 1035株/hm2；第三階段為1997-2007年，每度新竹達到1 065 ~ 13 20株/hm2；第四階段為2009-2011年，每度新竹達1 350株/hm2以上，盡管2011年春旱比較嚴重，在進行新竹調查時，還發現山坡上部有新竹枯死的現象，但2011年仍然獲得了新竹1 350株/hm2。

通過對1979-2011年新竹調查數據的統計、方差分析，得出了各年度間毛竹新竹株數有極顯著的差異，如表1和表2。

從表1、表2可知，年度間毛竹新竹株數的變化差異達到極顯著水平。為了進一步分析各年度之間的變化情況，把各年度之間新竹株數的平均數、樣地數及標準差列于表3。并對表3各年度之間的歷年數據分別右上和左下部分進行了歷年兩兩的F檢驗和t檢驗。通過右上半部分的F檢驗，發現大部分兩兩年度間的方差相近，但也有一些相差很大。基于F檢驗的結果，本文在進行t檢驗時，采用兩個總體平均數方差不知道和不相等的公式進行t檢驗。左下數字表示的是歷年每兩年間的t值檢驗。從兩兩比較來看，1999年是一條明顯的分界線，與往年比較都有極顯著的差異。

表1 歷年毛竹新竹株數統計Table 1 New culms

表2 毛竹新竹株數單因子方差分析Table 2 ANOVA on new culms number

表3 毛竹林各年度新竹株數變化情況t值、F值檢驗Table 3 Comparisons of new culms numbers every two years (t test, F test)株/hm2

為了更直觀地表達新竹株數的變動情況，把歷年新竹的平均株數繪制成圖1。從表3和圖1可知，在1979-2011年，毛竹林新竹株數有顯著地增加。

3.2 毛竹新竹胸徑的變化

從表4可以看出：在毛竹林新竹株數增加的同時，毛竹林新竹的平均胸徑不但沒有下降，反而有加粗的趨勢（圖2）。

圖1 歷年新竹株數變動情況Figure 1 Change of new culm number every two years

表4 1979—2011年毛竹林每一度新竹胸徑變化情況（t值、F值檢驗）Table 4 Comparisons of DBH of new culms of every two years (t-test, F-test)

從表5、表6可知：1979-2011年間，毛竹林的新竹平均胸徑變化達到了極顯著水平。從歷年的新竹調查原始資料來看，總的趨勢是隨著時間的推移，新竹胸徑的眾位數都在向增粗方向移動。

在表4的右上和左下部分，分別進行了歷年兩兩的F值和t檢驗。左下數字列出了歷年每兩年間的t檢驗，1987年與所有年份都具有極顯著差異，可能的原因之一是在所調查分析的14度竹里，1986年是唯一全面留養小年竹的年份，由于小年竹平均胸徑較小，從而降低了這一度毛竹的平均胸徑。

圖2 歷年新竹胸徑的變化情況Figure 2 Change of DBH of new culm every two years

表5 歷年毛竹新竹胸徑統計Table 5 DBH of new culm

表6 毛竹新竹胸徑單因素方差分析Table 6 ANOVA on DBH of new culm

對歷年的新竹竹材產量作統計和方差分析，結果如表7、表8。

表7 歷年新竹產量統計Table 7 Yield of new culms

表8 竹材產量單因子方差分析Table 8 ANOVA on yield of new culms

由表8方差分析結果表明：毛竹新竹竹材產量在各年度之間存在著極顯著差異。

從以上的新竹株數和新竹的胸徑變化，可以推斷出新竹竹材產量是增產了，事實上增產的幅度也是相當大的。

表9 歷年毛竹新竹產量t檢驗、F檢驗Table 9 Comparisons of new culm yield every two years (t-test, F-test)

從歷年新竹竹材產量變化圖（圖3）可以看到：自1999年至今，每一度新竹竹材產量雖有起伏，但都維持在21 000 kg/hm2以上，明顯地高于往年，達到了高產穩產。

圖3 歷年毛竹新竹竹材產量變化Figure 3 Change of new culm yield every two years

4 小結與討論

經過30余a的培育，安吉縣天荒坪鎮銀坑塢毛竹林的毛竹株數、胸徑、產量都有較大幅度的增長。這說明各個時期的毛竹林培育技術措施與當時的毛竹林狀態是相適應的、是有效的。1980年代初，鑒于當時留養三度竹，毛竹立竹較稀（2 400株/hm2），鉤梢過度等情況，相應地提出了“留養四度竹，嚴格控制鉤梢強度”的技術措施，要求做到“四度留養，五度環空，禁砍‘冷刀’，鉤梢留枝15 ~18盤”。通過這些措施的實施，增加了毛竹林葉面積指數，調整了毛竹林分結構，同時采用劈山、松土、施肥等技術，為毛竹林的生長營造了良好的立地條件。到1980年代末，立竹達3 000株/hm2以上，隨著立竹的逐年遞增，毛竹林密度進一步加大，1995年天荒坪林業服務總公司適時提出了“留二砍三”的毛竹林采伐新技術，使竹林的密度保持在一個合適的范圍內，同時使竹林的結構更趨年輕化，保持了竹林的高產穩產。

[1] 潘春霞，李雪濤，呂玉龍. 安吉縣毛竹資源及其生物量研究[J]. 浙江林業科技，2010，30（1）：82-84.

[2] 楊校生，吳良如，李正才，等. 毛竹林經濟和生態公益價值綜合評價[J]. 竹子研究匯刊，2007，26（1）：1-5.

[3] 周國模，姜培坤．毛竹林的碳密度和碳貯量及其空間分布[J]. 林業科學，2004，40（2）：20-24.

[4] 周芳純. 竹林培育和利用[M]. 北京：中國林業出版社，1998.

[5] 黃伯惠，胡正堅. 毛竹周年采伐及經營管理配套技術研究[J]. 竹子研究匯刊，1996，15（4）：52-67.

[6] 朱劍秋，曾新宇. 毛竹林不同留養度數對生長和產量的影響[J]. 竹子研究匯刊，1990，9（1）：55-62.

[7] 浙江省林業志編纂委員會. 浙江省林業志[M]. 北京：中華書局，2001. 648-719

[8] 周文偉，龔立忠，陳永華，等. 毛竹林松土效益分析[J]. 浙江林業科技，1992，12（4）：42-46.

[9] Gopal K. Kanji. 100 Statistical Tests[M]. London: SAGE Publications Inc, 1999.

Long-term Follow-up Investigation on Phyllastachys heterocycla var. pubescens Stands in Anji

ZHOU Wen-wei1,2，TONG Xiao-qing1,2，ZHOU Ren-ai3，HUA Xi-qi1,2，LIANG Yi-pin4，YU Hong-qing4，REN Bo-qing4
（1. Zhejiang Forestry Academy, Hangzhou 310023, China; 2. Zhejiang Provincial Key Laboratory of Bamboo Research, Hangzhou 310023, China; 3. Longquan Urban F orestry Station of Zhe jiang, Longquan 323700, China; 4. Anji Tianhuan ping Forestry Company of Zhejiang, Anji 313302, China）

Analysis was made on data of sampling investigation on Phyllastachys heter ocycla var. pubes cens stands from 1979-2011(missing 1989,1991 and 1993) at Yinkengwu of Tianhuangping town, Anji county of Zhejiang province. The result demonstrated that unit area density was increase, from 600 culms/ha in 1979 to 1350 culms/ha in 2011, topped to 1440 culms/ha in 2009. Mean DBH of bamboo increased from 8.75cm in 1979 to 9.08cm in 2011, and topped to 9.47cm in 2009. Culm yield increased steadily, from 12945 kg/ ha in 1979 to 25 680kg/ ha in 2011, with the highest of 29 790kg/ ha in 2009.

Phyllastachys heterocycla var. pubescens stand; follow-up investigation; statistical analysis

SS757.2+9

A

1001-3776（2012）02-0006-06

2011-11-21；

2012-02-18

周文偉（1959-），男，浙江麗水人，副研究員，從事竹類資源、培育研究。