劉 軍，白天道，陳永力，楊開(kāi)太，蔣維昕

(1.廣西國(guó)有七坡林場(chǎng)，廣西 南寧 210000； 2.廣西大學(xué)林學(xué)院，廣西 南寧 210000； 3.廣西林業(yè)科學(xué)研究院，廣西 南寧 210000)

我國(guó)擁有叢生竹種資源100多種，其中包含品質(zhì)優(yōu)良的鮮食型筍用竹種，如版納甜龍竹(Dendrocalamushamiltonii)、勃氏甜龍竹(D.brandisii)等；加工型筍用竹種，如麻竹(D.latiflorus)；鮮食和加工兼用的優(yōu)良竹種，包括綠竹(Dendrocalamopsisoldhami)、大頭典竹(D.beecheyana)、吊絲竹(Dendrocalamusminor)等[1]。

這些大型叢生竹種主要分布于熱帶、南亞熱帶與中亞熱帶氣候過(guò)渡區(qū)，為我國(guó)重要經(jīng)濟(jì)竹種[2]。竹筍味香質(zhì)脆，而且富含多種人體所必需的氨基酸和微量元素，是一種高蛋白低脂肪且營(yíng)養(yǎng)齊全的綠色保健食品，深受現(xiàn)代人追捧[3]。我國(guó)西南、華南、海南及東南沿海區(qū)域開(kāi)展了一系列大型叢生竹的竹種選擇、高效培育及竹筍品質(zhì)等研究[4-7]，尤以麻竹、綠竹、勃氏甜龍竹為生長(zhǎng)及發(fā)筍表現(xiàn)較佳的優(yōu)良筍材兩用竹種。廣西地處中國(guó)竹林面積最大、竹資源最豐富的長(zhǎng)江—南嶺竹區(qū)，其中，桂林和柳州為廣西最有名的兩大產(chǎn)竹區(qū)，均以散生竹種毛竹(Phyllostachysedulis)冬筍為市場(chǎng)主要食用筍種(53.3%)，其次以桂西北的田林八渡筍(麻竹)最為出名。余慧連等[8]在廣西南寧香蕉跡地研究不同施肥對(duì)綠竹(馬蹄筍)產(chǎn)量的影響，以評(píng)估其在桂南部農(nóng)業(yè)用地的發(fā)展?jié)摿?，結(jié)果表明，施肥對(duì)1年生馬蹄筍竹叢產(chǎn)筍質(zhì)量和數(shù)量均無(wú)顯著影響。而于增金等[9]對(duì)不同地類(lèi)種植麻竹的竹筍品質(zhì)研究表明，林地竹筍的纖維素、木質(zhì)素含量均低于農(nóng)業(yè)用地和河岸沖積地的竹筍。因而桂中、南部大面積低丘林下發(fā)展叢生筍用竹具有較大潛力。加勒比松是廣西南寧七坡林場(chǎng)主栽引種樹(shù)種之一，大面積加勒比松純林的營(yíng)造最終會(huì)影響其人工林生態(tài)系統(tǒng)的長(zhǎng)期穩(wěn)定性和持久生產(chǎn)力。當(dāng)前林下經(jīng)濟(jì)已成為現(xiàn)代林業(yè)產(chǎn)業(yè)發(fā)展的主流之一，對(duì)于鞏固拓展脫貧攻堅(jiān)成果、全面推進(jìn)鄉(xiāng)村振興有著特殊的意義[10-11]。本研究選擇適生于我國(guó)華南竹區(qū)的叢生筍用竹種：綠竹、麻竹、勃氏甜龍竹，研究3個(gè)竹種在廣西南寧低丘松林下的造林效果、發(fā)筍成竹規(guī)律及退筍情況，以期為林下叢生筍用竹套種新模式和經(jīng)營(yíng)管理提供參考。

1 材料與方法

1.1 試驗(yàn)地概況

試驗(yàn)地位于廣西壯族自治區(qū)南寧市吳圩鎮(zhèn)的廣西國(guó)有七坡林場(chǎng)立新分場(chǎng)，東經(jīng)107°59′ —108°21′、北緯22°28′—22°68′，海拔200～300 m。年均氣溫21～23 ℃，≥10 ℃的年積溫7500 ℃，極端最低溫-2.6 ℃，極端最高溫38.4 ℃，干濕季節(jié)明顯，其中4—9月為雨季，年均降水量1200～1300 mm，年相對(duì)濕度為70%～80%。該區(qū)域?qū)賮啛釒ШＱ笮约撅L(fēng)氣候，溫暖濕潤(rùn)、雨量充沛，極為適合叢生竹的生長(zhǎng)。試驗(yàn)栽培地為低丘陵下間伐后的低效加勒比松人工純林，林齡28 a，林分現(xiàn)存密度為6 m×6 m，郁閉度0.6，林下灌木層零星分布有五指毛桃(Ficussimplicissimavar.hirta)、三叉苦(Evodialepta)、地桃花(Urenalobata)等，草本層優(yōu)勢(shì)種有求米草(Oplismentlsundulatifolius)、五節(jié)芒(Miscanthusfloridulus)、華南毛蕨(Cyclosorusparasiticus)、海金沙(Lygodiumjaponicum)以及人工種植林藥天冬(Asparaguscochinchinensis)，長(zhǎng)勢(shì)良好。林地坡度為20°～30°，土壤以頁(yè)巖、砂巖發(fā)育的赤紅壤為主，土壤呈酸性，表土pH值小于4.5。

1.2 試驗(yàn)材料

綠竹(Dendrocalamopsisoldhami，簡(jiǎn)稱(chēng)L)、麻竹(Dendrocalamuslatiflorus，簡(jiǎn)稱(chēng)M)、勃式甜龍竹(Dendrocalamusbrandisii，簡(jiǎn)稱(chēng)T)3個(gè)竹種均為帶土球的1年生帶蔸主枝扦插苗，地徑約2～2.5 cm，竹苗截稈高度為30～40 cm。

1.3 種植設(shè)計(jì)

采用裂區(qū)設(shè)計(jì)，主處理為竹種(M、T、L)，副處理為坡位(上、中、下坡位)，每個(gè)坡位3個(gè)重復(fù)(區(qū)組)，每個(gè)重復(fù)含3個(gè)裂區(qū)，每個(gè)裂區(qū)為1個(gè)竹種(0.1 hm2)，每個(gè)重復(fù)0.3 hm2，總試驗(yàn)面積約2.7 hm2。種植密度為450 株·hm-2，每種竹子405株，松樹(shù)林地不需整地，以垂直縱向開(kāi)穴種植，種植穴約50 cm×60 cm，深30 cm左右。種植時(shí)間為2019年3月初，5月初進(jìn)行除草撫育，同時(shí)溝施氮磷鉀復(fù)合肥(N∶P∶K=5∶2∶1)200 g·株-1。

1.4 調(diào)查方法

根據(jù)試驗(yàn)地的種植狀況，采用樣方法進(jìn)行調(diào)查，每個(gè)坡位每個(gè)重復(fù)內(nèi)設(shè)置固定的矩形標(biāo)準(zhǔn)地(40 m×30 m)。種植后180 d(2019年9月30日)進(jìn)行標(biāo)準(zhǔn)地內(nèi)母竹成活率調(diào)查。筍期調(diào)查從5月20日開(kāi)始，以30 d為一個(gè)觀(guān)察單元，貫穿整個(gè)筍期。借助游標(biāo)卡尺和測(cè)高器等工具記錄各竹種發(fā)筍情況，以筍尖露出地面≥3 cm為發(fā)筍標(biāo)準(zhǔn)，以竹籜葉失水變干、高生長(zhǎng)停止為退筍標(biāo)準(zhǔn)，主要記錄每個(gè)樣方內(nèi)竹種造林母竹成活數(shù)、發(fā)筍數(shù)、退筍數(shù)、退筍高度等。進(jìn)入發(fā)筍盛期(7月22日)，于下坡位各竹種分別選擇16株出土一致(5～10 cm)且生長(zhǎng)狀況良好的筍掛標(biāo)簽編號(hào)，進(jìn)行連續(xù)90 d(固定同一時(shí)間，每2 d測(cè)量1次)的竹筍—幼竹高生長(zhǎng)調(diào)查。

1.5 數(shù)據(jù)統(tǒng)計(jì)分析

利用Excel進(jìn)行初步計(jì)算，SPSS 23.0進(jìn)行統(tǒng)計(jì)分析。相關(guān)統(tǒng)計(jì)公式為：?jiǎn)沃昴钢癜l(fā)筍率=發(fā)筍數(shù)/母竹總數(shù)；單株母竹退筍率=退筍數(shù)/母竹總數(shù)；發(fā)筍量=樣方內(nèi)累計(jì)發(fā)筍數(shù)量；退筍量=樣方內(nèi)累計(jì)退筍數(shù)量；單株母竹成筍率=單株發(fā)筍率-單株退筍率；造林成活率=(成活株數(shù)/總的母竹株數(shù))×100%；蟲(chóng)害發(fā)生率=(受感染母竹數(shù)量/樣方內(nèi)種植母竹數(shù)量)×100%。

不同階段發(fā)筍率及退筍率采用R軟件進(jìn)行方差分析。使用R語(yǔ)言自編程序?qū)Πl(fā)筍—幼竹高生長(zhǎng)進(jìn)行擬合，基于AIC信息準(zhǔn)則(Akaike information criterion，AIC)衡量模型擬合優(yōu)良性(Goodness of fit)。通常情況下，AIC定義為：AIC=2k-2ln(L)，其中k為模型參數(shù)個(gè)數(shù)，L為似然函數(shù)。-2ln(L)反映模型的擬合情況，當(dāng)2個(gè)模型之間存在較大差異時(shí)，差異主要體現(xiàn)在似然函數(shù)項(xiàng)-2ln(L)，當(dāng)似然函數(shù)差異不顯著時(shí)，模型參數(shù)的懲罰項(xiàng)2k則起作用，隨著模型中參數(shù)個(gè)數(shù)增加，2k增大，AIC增大，從而參數(shù)個(gè)數(shù)少的模型是較好的選擇。

2 結(jié)果與分析

2.1 3個(gè)筍用竹種在廣西南寧的造林成活及蟲(chóng)害發(fā)生率

由表1可知，3個(gè)叢生筍用竹種在坡地種植的成活率都達(dá)到了60%以上，綜合不同坡位，勃氏甜龍竹的成活率高達(dá)83.15%，明顯高于綠竹(76.6%)和麻竹(62.13%)。不同坡位對(duì)竹種成活率的影響不明顯。進(jìn)入發(fā)筍期(6—8月)，南方多雨濕潤(rùn)，大量蚜蟲(chóng)、介殼蟲(chóng)及竹螟蟲(chóng)危害竹葉，幾個(gè)叢生筍用竹種均遭不同程度損害，以綠竹受蟲(chóng)害感染最為嚴(yán)重，高達(dá)57.92%；勃氏甜龍竹和麻竹蟲(chóng)害發(fā)生率接近(32.32～33.56%)。

表1 不同筍用叢生竹造林成活率及蟲(chóng)害發(fā)生率 %

2.2 不同階段發(fā)筍及退筍情況

發(fā)筍率及成筍率在不同階段不同筍用竹種間，以及退筍率在不同階段均存在極顯著性差異(P<0.001)(表2)。勃氏甜龍竹發(fā)筍時(shí)間最早(5月中旬)，綠竹的發(fā)筍時(shí)間最晚(6月上旬)，發(fā)筍持續(xù)時(shí)間歷時(shí)最長(zhǎng)的為勃氏甜龍竹(5月中旬—11月中旬)，持續(xù)時(shí)間最短的為麻竹(5月下旬—10月下旬)。3個(gè)竹種均以6—7月為發(fā)筍初期，8—9月為發(fā)筍盛期，10—11月為發(fā)筍晚期。由圖1可看出，勃氏甜龍竹(T)、綠竹(L)、麻竹(M)平均每株母竹發(fā)筍率/退筍率，發(fā)筍初期分別為0.468/0.021、0.256/0.009、0.483/0.045；發(fā)筍盛期分別為1.630/0.089、1.399/0.083、2.186/0.115；發(fā)筍晚期分別為1.667/0.032、1.542/0.091、1.911/0.059。勃氏甜龍竹(T)和綠竹(L)晚期的單株母竹發(fā)筍率及成筍率均較發(fā)筍盛期趨于穩(wěn)定，而麻竹(M)發(fā)筍率進(jìn)入發(fā)筍晚期后相對(duì)于盛期降低，這與麻竹在盛期、晚期仍然保持較高的退筍率有關(guān)。另外，勃氏甜龍竹(T)和麻竹(M)在發(fā)筍盛期，退筍率達(dá)到高峰；進(jìn)入發(fā)筍晚期，退筍率也隨之下降，而綠竹退筍率在發(fā)筍晚期仍然持續(xù)上升。

表2 3個(gè)竹種不同發(fā)筍階段單株母竹發(fā)筍率及退筍率方差分析

2.3 不同坡位發(fā)筍及退筍情況

對(duì)不同坡位3個(gè)筍用竹種發(fā)筍盛期(8月下旬)的單株發(fā)筍率及退筍率進(jìn)行分析(圖2)。結(jié)果表明，3個(gè)竹種的單株發(fā)筍率均以下坡位略高于中、上坡位。對(duì)于退筍率，不同竹種均以上坡位的退筍現(xiàn)象最為明顯，其中麻竹、綠竹的單株退筍率>0.1，單株退筍率分別為0.138、0.13。可能與上坡位土層相對(duì)瘠薄，石礫含量較高，水分不足，導(dǎo)致竹筍發(fā)育不良而退敗。

橫坐標(biāo)中T、L和M分別為勃氏甜龍竹、綠竹和麻竹圖1 3個(gè)竹種不同發(fā)筍階段單株母竹發(fā)筍率及退筍率

圖2 不同坡位3種筍用竹種發(fā)筍率及退筍率

2.4 筍—幼竹高生長(zhǎng)規(guī)律

由圖3可看出，3個(gè)竹種的筍—幼竹高生長(zhǎng)擬合曲線(xiàn)整體上都呈“S”型，竹種生長(zhǎng)符合“慢—快—慢”的生長(zhǎng)趨勢(shì)。以始末高生長(zhǎng)量占高生長(zhǎng)總生長(zhǎng)量的10%為界限，整個(gè)筍期可分為初期、盛期、末期3個(gè)時(shí)期。3個(gè)竹筍的高生長(zhǎng)初期均在17 d內(nèi)，該階段生長(zhǎng)緩慢，每日高生長(zhǎng)為2 cm；筍出土后的17～61 d為盛期，這一階段生長(zhǎng)較快，平均每日生長(zhǎng)11 cm；筍出土后的62～90 d為生長(zhǎng)末期，平均每日生長(zhǎng)2 cm。不同竹種發(fā)筍—幼竹的高生長(zhǎng)持續(xù)時(shí)間有所差異，勃氏甜龍竹、麻竹、綠竹竹筍平均高生長(zhǎng)持續(xù)時(shí)間從高至低依次為82、79、74 d。勃氏甜龍竹和麻竹的高生長(zhǎng)速率相接近，生長(zhǎng)盛期凈高生長(zhǎng)量平均分別為520、547 cm，日均高增長(zhǎng)量為12 cm，高于綠竹的生長(zhǎng)盛期凈高生長(zhǎng)量(458 cm，日均高增長(zhǎng)量10 cm)。

圖中虛線(xiàn)為樣條函數(shù)擬合曲線(xiàn)，實(shí)線(xiàn)為最佳生長(zhǎng)方程擬合曲線(xiàn)圖3 3種叢生竹的筍—幼竹高生長(zhǎng)規(guī)律

采用數(shù)學(xué)模型模擬竹子的高生長(zhǎng)已被廣泛采用。本文以3種叢生筍用竹竹筍—幼竹株高為因變量，發(fā)筍后的天數(shù)為自變量，采用Richards、Logistic和Gomperz.T方程分別對(duì)不同竹種萌發(fā)筍的高生長(zhǎng)過(guò)程進(jìn)行模擬，基于最小AIC值獲得最優(yōu)模型(表3)。結(jié)果表明，Richards較好地模擬勃氏甜龍竹和綠竹的高生長(zhǎng)過(guò)程，Gomperz.T則較好地模擬麻竹的高生長(zhǎng)過(guò)程。最佳模擬方程分別為：勃氏甜龍竹：yT=722.2519(1-0.22076exp(1.418142-0.050794t))4.529806；綠竹：yL=590.5062(1+0.45790exp(2.586007-0.071649t))-2.18388；麻竹：yM=683.68380exp(-exp(1.932864-0.06847t)。

表3 3個(gè)叢生竹種高生長(zhǎng)曲線(xiàn)擬合模型比較

3 討論與結(jié)論

本試驗(yàn)研究表明，不同竹種發(fā)筍周期有所差異。3種叢生竹分別從5月中旬—6月上旬開(kāi)始發(fā)筍，以勃氏甜龍竹的發(fā)筍持續(xù)時(shí)間歷時(shí)最長(zhǎng)(180 d)，麻竹的發(fā)筍持續(xù)時(shí)間最短(150 d)。根據(jù)報(bào)道，勃氏甜龍竹在滇西—南5個(gè)縣的發(fā)筍歷期平均在120 d左右[12]；在云南雙江、勐海、彌勒的發(fā)筍歷期分別為153、182、124 d[7]；而引種至溫州地區(qū)發(fā)筍期持續(xù)150 d左右[13]；綠竹在其最適宜區(qū)福建福安的發(fā)筍期在5月下旬—10月初[14]。本試驗(yàn)中，綠竹在廣西南寧的發(fā)筍時(shí)間稍晚(6月初)，但發(fā)筍結(jié)束期在11月中旬，歷時(shí)相對(duì)較長(zhǎng)(165 d)。可見(jiàn)，相同的竹種，其發(fā)筍規(guī)律受地理分布、氣候因素、不同年份和立地條件的影響。從整體來(lái)看，筍期可以分為初期(6—7月)、盛期(8—9月)、末期(10—11月)3個(gè)階段。這與大部分叢生竹的發(fā)筍節(jié)律一致[7，12-13]。退筍量最多的時(shí)期是在8—9月，因其發(fā)筍量也多，故退筍也多。不同竹種間單株發(fā)筍率/退筍率存在顯著差異(P<0.05)，以麻竹的單株發(fā)筍率/退筍率最高，發(fā)筍初期、盛期、末期分別為0.483/0.045、2.186/0.115、1.911/0.059；其次為勃氏甜龍竹，分別為0.468/0.021、1.630/0.089、1.667/0.032；綠竹分別為0.256/0.009、1.399/0.083、1.542/0.091。

不同坡位中，3個(gè)竹種的單株發(fā)筍率均以下坡位略高于中、上坡位；而退筍率，不同竹種均以上坡位的退筍現(xiàn)象最為明顯。從整體來(lái)看，上坡位的成筍率相對(duì)較低，下坡位由于土層肥沃，水分充足，成筍率相對(duì)較高。3個(gè)竹種間，以麻竹的發(fā)筍對(duì)坡位響應(yīng)最為敏感，表明在山頂(上坡位)等土層瘠薄的地段較不適宜培育麻竹筍用林。已有研究指出，無(wú)論叢生竹還是散生竹，土層肥沃，水分、養(yǎng)分充足是培育優(yōu)質(zhì)筍用林的重要保證[2]。在竹林管理過(guò)程中，可針對(duì)不同坡位采取不同經(jīng)營(yíng)模式，如坡下部及夾溝處，采取集約經(jīng)營(yíng)，實(shí)現(xiàn)高效筍竹培育；對(duì)于土層瘠薄、水肥較差的中上坡位，可采取近自然經(jīng)營(yíng)，打造山區(qū)生態(tài)竹筍精品，亦可提升品牌效益[15]。

勃氏甜龍竹、麻竹、綠竹幼竹在南寧的生長(zhǎng)歷期74～82 d，與譚宏超等[7]在云南(80 d以后高生長(zhǎng)趨緩至停)觀(guān)察到的幾種叢生竹種的生長(zhǎng)周期接近。而夏海濤等[13]研究幾種叢生竹種在浙江溫州的發(fā)筍及高生長(zhǎng)規(guī)律，觀(guān)測(cè)到勃氏甜龍竹的高生長(zhǎng)在40 d后基本停止，表明幼竹的高生長(zhǎng)速率與氣候條件密切相關(guān)[7，13]。勃氏甜龍竹、綠竹、麻竹幼竹在南寧的高生長(zhǎng)均呈現(xiàn)慢—快—慢“S”型規(guī)律，這與大多數(shù)竹種[6-7，16-18]類(lèi)似。其中，勃氏甜龍竹與麻竹的高生長(zhǎng)速率相接近，生長(zhǎng)盛期凈高生長(zhǎng)量約為綠竹的1.2 倍。有研究表明，Logistic[16-17]及Doseresp[18]曲線(xiàn)模型可以較好地模擬部分散生竹種的高生長(zhǎng)過(guò)程。本試驗(yàn)中，采用多項(xiàng)模型(Richards、Logistic、Gomperz.T)進(jìn)行評(píng)估，依據(jù)AIC權(quán)衡估計(jì)模型復(fù)雜度和擬合數(shù)據(jù)優(yōu)良性。根據(jù)AIC值顯示，Richards較好地模擬廣西南寧地區(qū)種植的勃氏甜龍竹和綠竹高生長(zhǎng)過(guò)程，Gomperz.T則較好地模擬該區(qū)域麻竹的高生長(zhǎng)過(guò)程。