竹筍采挖下克隆整合對毛竹林發筍成竹規律的影響

摘 要：【目的】為研發毛竹林高產高效生態經營技術提供一定參考。【方法】調查浙西南地區采挖毛竹林竹筍后的筍- 幼竹生長過程，通過存活曲線分析毛竹林挖筍后的發筍成竹規律，探討挖筍對竹筍中單寧、粗脂肪、粗蛋白等品質指標的影響。【結果】竹筍采挖提高了毛竹林地的發筍能力，挖筍樣地的發筍數比未挖筍樣地增加了12.8%（P ＜ 0.05），而挖筍后選留竹筍的存活率也顯著高于未挖筍樣地；與未挖筍樣地相比，挖筍樣地內毛竹新分株的平均胸徑并未發生顯著改變（P ＞ 0.05），挖筍不會降低新分株胸徑而使林分衰退；竹筍間的非對稱性競爭有利于竹筍地徑的徑級為Ⅰ - Ⅲ級竹筍的生長與存活，隨著地徑的提高竹筍死亡率顯著降低，而徑級Ⅳ級以上竹筍則有著較高的死亡率。盛筍早期以后，徑級相對較小（Ⅳ級以上）的竹筍開始大量萌發，這類竹筍若不及時采挖，大部分會出現退筍，造成資源浪費；挖筍和未挖筍樣地竹筍的含水率、蛋白質和粗脂肪含量在各筍期均未見顯著差異（P ＞ 0.05），盡管挖筍樣地竹筍的單寧含量在末筍期顯著高于未挖筍樣地，但兩種樣地其他筍期單寧含量均未見顯著差異。【結論】毛竹林竹筍的采挖，可緩解竹筍間的非對稱性的資源競爭，有效避免大量退筍導致的毛竹母株營養消耗，實現母株資源的克隆整合，增加竹筍產量。因此，合理的竹筍采挖有利于毛竹筍用林的高效生產及可持續性經營。

關鍵詞：毛竹林，竹筍采挖，非對稱性競爭，克隆整合

中圖分類號：S644.2；727.39 文獻標志碼：A 文章編號：1003—8981（2024）01—0294—10

資源（光照、養分、空間）短缺是植物生存的重要限制因素，自然群落中同一物種的相鄰個體因具有同樣生長習性，往往發生資源競爭[1]。目前，常以“負密度制約假說”來解釋同種個體間相互作用與共存的機制[2]，在生物侵害、資源競爭等因素的影響下，個體在同種密度高的生境中會出現較低的存活率。毛竹Phyllostachys edulis 屬禾本科竹亞科，為高大喬木狀植物，具有生理整合、形態可塑等克隆生態學特性[3-4]，其依靠地下鞭側芽的萌發完成克隆更新，發筍成竹。毛竹在從竹筍到幼竹的生長過程中很大程度上不直接需要光資源，適度的遮光反而有利于幼竹生長[5]。竹筍生長所需的資源主要來源于母株及地下鞭根系統貯藏的光合同化物、礦質養分及水分，母株通過生理整合將這些資源由鞭根系統傳輸至子代新分株，支持子代分株存活與生長，進而拓展生存空間。有研究表明[6]，毛竹林母株數量越多，能夠提供的整合資源就會越充足，竹筍的發筍數和存活率也就越高。因此竹筍存活率是非密度制約的。但母株所能提供資源畢竟有限，隨著毛竹密度和發筍數的增加，竹筍對資源的爭奪也會加劇。同一鞭根系統相鄰的竹筍間往往存在非對稱性競爭，鄒躍國等[7] 認為，較大個體的竹筍在資源爭奪中更占優勢，并會抑制小筍存活。母株將貯藏的資源優先供給那些生長健壯的竹筍，竹筍個體越強壯競爭能力越強。

竹筍采收是毛竹林培育的一項重要經營活動。挖筍后，竹筍間原有的競爭關系被打破，母株貯藏的養分、水分等資源會發生重新分配。此外，竹筍的品質與養分等環境資源關系密切，而母株對資源的再次整合也會對竹筍品質產生影響。近年來，已有較多關于竹類植物的竹筍- 幼竹生長動態規律的研究[8-11]，而尚未見關于竹筍采收能否緩解竹筍間資源的非對稱競爭、如何影響毛竹無性種群的數量和動態變化規律、是否會改變竹筍品質方面的報道。本文調查了浙西南地區毛竹林挖筍后筍- 幼竹生長動態過程，對出筍時間、出筍數、成竹率、退筍率等發筍成竹規律進行分析，并探討了其與單寧、含水率、粗脂肪、粗蛋白等品質指標的關系，可為毛竹林高產高效生態經營技術提供一定參考。

1 材料和方法

1.1 研究區概況

試驗地設在浙江省遂昌縣，地理坐標為東經118°41′ ～ 119°30′，北緯28°13′ ～ 28°49′，氣候屬中亞熱帶季風類型，山地垂直氣候差異明顯。年平均氣溫16.8 ℃，年降水量1 510 mm，年降水日數172 d，年太陽總輻射量為101 kcal/cm2，年日照時數1 755 h，年無霜期251 d。試驗樣地的土壤類型為山地黃紅壤，土層厚度＞ 100 cm。

2021 年在遂昌縣三仁鄉選取立地條件、立竹密度和竹林結構相近的毛竹筍竹兩用林（采用大小年經營模式，2020 年為冬筍大年并于2021年春留養新竹）。施肥狀況為冬筍大年5 月施尿素700 kg/hm2、過磷酸鈣150 kg/hm2、氯化鉀200 kg/hm2，連續施肥10 a。試驗樣地土壤的基本理化性質見表1。

1.2 樣地設置與測定方法

設置挖筍（對冬筍全部采挖，春筍除選留成竹外，其余全部采挖）和不挖筍（對冬筍和春筍均不采挖）2 種處理。每種處理6 塊樣地，樣地面積為15 m×25 m。

2021 年3 月初，毛竹林地春筍陸續出土，3 月4 日開始第一次調查統計，此后于每天上午8： 00 進行調查。每次調查時在新出土的竹筍（出土5 cm）旁插標簽并加以編號，并對每一竹筍的情況（存活、死亡或被采挖）進行跟蹤記錄，直到出筍結束。開始退筍后進行樣地退筍數量統計，并測定退筍地徑。存活竹筍長成的新分株測定其胸徑。

于各出筍期在每塊樣地內隨機挖取12 個大小基本一致、無機械損傷、出土20 cm 左右的筍樣分析品質：剝去筍籜，取可食部分切成小塊混合保存。竹筍水分含量按GB/5009.3-2010 的方法測定，粗脂肪含量按GB/T5009.6-2003 的方法測定，蛋白質含量按GB/T5009.5-2010 的方法測定，單寧含量按照Folin － Ciocalteu 比色分光光度法測定。

1.3 數據處理

筍期的劃分： 建立Logistic 生長曲線Y=K/（1+ae-bt），擬合竹筍累計出筍數沿著出筍時間序列的變化趨勢。并用Logistic 生長曲線的3 個關鍵點t1=（lna-1.317）/b、t2=lna/b、t3=（lna+1.317）/b， 將毛竹筍的出筍期劃分為4 個同生群，分別為初筍期、盛筍早期、盛筍后期和末筍期。分別對挖筍和不挖筍樣地不同出筍期的出筍數、存活率及存活的新生分株胸徑進行單因素方差分析，探討挖筍對出筍成竹規律的影響。

存活曲線：按竹筍地徑（D）大小將其分為7級，D ≥ 14 cm（Ⅰ級）、12 ≤ D ＜ 14 cm（Ⅱ級）、10 cm ≤ D ＜ 12 cm（Ⅲ級）、8 cm ≤ D ＜ 10 cm（Ⅳ級）、6 cm ≤ D ＜ 8 cm（Ⅴ級）、4 cm ≤ D ＜6 cm（Ⅵ級）、0 ＜ D ＜ 4 cm（Ⅶ級）。Nx 為在X 徑級時，減去該徑級退筍的數量后樣地內竹筍的存活數量，N0 為初始出筍數，1 000 Nx/ N0 為標準化后的存活數量l x。以l x 為縱坐標、竹筍地徑的徑級為橫坐標，繪制存活曲線[12]。

Logistic 生長曲線的擬合及單因素方差分析均在SPSS16.0 軟件中完成，存活曲線繪制用Origin7.5 軟件完成。

2 結果和分析

2.1 挖筍對出筍成竹規律的影響

毛竹林中竹筍采挖包括冬筍和春筍。毛竹的筍芽萌發后，至初冬，筍體逐漸膨大，筍籜呈黃色，被有絨毛，此時的竹筍進入越冬期，在土中處于休眠狀態，是為冬筍。翌年春季地溫上升，冬筍打破休眠，出土生長，是為春筍。2020 年11 月至2021年1 月中旬，將挖筍處理的樣地內未露土的冬筍全部采挖，平均采挖數為3 628.5±213.5 株/hm2。2021 年3 月初竹筍露土，開始采挖春筍并留養新分株。從表2 可以看到，挖筍樣地春筍的筍期為3 月7 日至4 月26 日，而未挖筍樣地筍期為3 月15 日至4 月21 日，說明毛竹林地冬筍采挖可使春筍的出土時間提前，整個筍期也有所延長。在毛竹林地采挖冬筍不僅沒有抑制春筍萌發，而且還提高了春筍的發筍能力，挖筍樣地的春筍發筍數比未挖筍樣地增加了12.8%（P ＜ 0.05）。

各樣地毛竹林春筍累計出筍數量與時間序列均呈S 曲線，符合種群密度控制的生長規律，可較好地用Logistic 生長方程進行擬合（R2=0.99），并根據曲線的拐點將筍期劃分為4 個同生群。由表2 可知，挖筍樣地的發筍高峰在盛筍后期，盡管未挖筍樣地的出筍期較晚，但從初筍期至盛筍早期，發筍數迅速達到峰值，隨后發筍數銳減直至末筍期停止發筍。

春筍出筍后，毛竹會迅速完成高生長，發枝展葉成為毛竹新分株，但在此期間往往也會因養分、水分不足導致退筍。挖筍后選留春筍的存活率顯著高于未挖筍樣地（P ＜ 0.05）。隨著時間的推移，從初筍期直至末筍期，挖筍與未挖筍處理樣地的存活率均表現為逐漸降低的趨勢。

2.2 挖筍對新分株胸徑的影響

圖1 為不同筍期毛竹新分株胸徑大小，由圖1可知，挖筍與未挖筍處理各筍期的新分株胸徑均存在一定差異，且規律一致，表現為初筍期和盛筍早期胸徑較大，隨著筍期推移胸徑逐漸減小。兩種處理下，末筍期胸徑分別比初筍期降低了29.1%（挖筍）、15.8%（未挖筍）（P ＜ 0.05）。與未挖筍樣地相比，挖筍并未使樣地內新分株的平均胸徑發生顯著改變（P ＞ 0.05），說明在毛竹林樣地連續施肥、土壤養分得到保證的情況下，挖筍這一經營行為并不會降低新分株平均胸徑而使林分發生衰退。

2.3 竹筍- 幼竹無性系種群存活曲線與退筍率

存活曲線是種群統計分析方法中一種比較直觀的方法，常用來反映不同齡級植物的存活狀況。本研究以竹筍地徑的徑級為橫坐標繪制存活曲線，分析不同徑級下竹筍的存活情況。從圖2 可知，未挖筍樣地的存活曲線呈凹曲線，該曲線類型的竹筍在徑級較小時（Ⅳ級以上）死亡率高，但隨著竹筍地徑徑級的增大，Ⅰ～Ⅲ級的竹筍死亡率有所降低。挖筍樣地在竹筍的人工選留過程中會盡量采去徑級較小的弱筍，保留壯筍（Ⅰ～Ⅳ級），徑級Ⅳ的竹筍死亡率為100%，而徑級為Ⅰ級的竹筍死亡率降低至21.6%。

圖3 為挖筍和未挖筍樣地不同筍期的竹筍- 幼竹存活曲線與退筍率。由圖可知，挖筍樣地在初筍期，由于養分充足，竹筍死亡數相對較少，不同地徑徑級間的竹筍有強烈的非對稱性競爭，大量徑級為Ⅰ級的竹筍存活，死亡率僅為16.7%，從而導致Ⅲ、Ⅳ級的竹筍死亡率為100%。隨著筍期的推移，養分消耗，死亡數增多，而Ⅰ級竹筍的發筍數減少，因此非對稱性競爭得以緩解。至盛筍后期，Ⅰ級竹筍死亡率升為30.9%，Ⅲ級竹筍死亡率降為33.3%，但Ⅳ級竹筍死亡率仍為100%。

未挖筍樣地中，除末筍期徑級Ⅳ級以上的竹筍死亡率為94.4% ～ 100% 外， 其他筍期均為100%。但存活曲線有所不同，初筍期發筍以Ⅰ - Ⅲ級壯筍為主，盡管Ⅳ級以上的竹筍死亡率為100%，由于數量少，存活曲線相對平緩，Ⅰ - Ⅲ級竹筍死亡率降低，但發筍數量多，存活曲線迅速下降。盛筍早期，Ⅳ - Ⅶ級竹筍增多，至盛筍后期和末筍期，發筍以Ⅳ - Ⅶ級級筍為主，Ⅰ - Ⅲ級壯筍數量較小，存活曲線的類型表現為凹曲線。

2.4 挖筍對竹筍品質的影響

竹筍含水率是竹筍口感與鮮嫩程度的重要指標。與未挖筍樣地相比，挖筍后竹筍含水率未出現顯著改變（P ＞ 0.05），而各筍期間竹筍含水率也未見顯著差異（P ＞ 0.05）。這可能是因為竹筍含水率的高低和降水量及土壤水分狀況等因素的關系更為密切[13]，因此挖筍行為對竹筍含水率的影響較小。

蛋白質和脂肪是竹筍營養成分中含量較高的組成部分，這些營養成分含量的高低也會影響竹筍品質[14]。食品中的脂肪可降低血液中膽固醇含量，具有預防心血管病，促進腦細胞發育等作用，而蛋白質是多種氨基酸組合而成的高分子化合物，為生物體的主要組成部分，也是人體不可缺少的營養物質。由表3 可知，挖筍和未挖筍樣地的竹筍蛋白質含量在初筍期、盛筍早期、盛筍后期相對穩定，末筍期均顯著降低（P ＜ 0.05）。而挖筍和未挖筍樣地間竹筍的蛋白質含量在各筍期均未見顯著差異（P ＞ 0.05），挖筍及筍期變化對竹筍中粗脂肪的含量均沒有顯著影響。說明挖筍活動對竹筍中蛋白質和脂肪等營養物質含量的影響不大。

單寧是多種多元酚衍生物的混合物，是植物次級代謝產物，受到初級代謝的調節控制。竹筍中所含單寧會引發澀感，影響竹筍的風味和口感。不同筍期單寧含量有所不同，在挖筍和未挖筍兩類樣地中，均表現為初筍期竹筍單寧含量相對較低，但隨著發筍時間的推移，至盛筍后期和末筍期單寧含量顯著提高（P ＜ 0.05）。挖筍樣地的竹筍在末筍期單寧含量顯著高于未挖筍樣地，但這兩類樣地的其他筍期單寧含量未見顯著差異（表3），這可能是因為挖筍活動打破了筍芽的非對稱競爭，增加了發筍的數量，釆筍林地對養分的消耗更多，而林地養分的虧缺會導致植物體內單寧含量的提高[15]。因此在末筍期竹筍體內積累了較多的單寧。

冗余分析（Redundancy analysis RDA）是將線性相關分析和排序分析相結合形成的多元排序分析，在對應迭代中，將每次的排序值與環境因子進行多元線性回歸，從而反映環境因子對排序結果的影響。通過對挖筍和未挖筍樣地不同筍期竹筍的品質指標數據進行排序（直接梯度分析），探討挖筍在不同筍期對竹筍品質指標的影響。從圖4 可知，第一軸和第二軸的貢獻率分別為52.4%、4.1%，可以解釋56.5% 以上原變量的特征。第一軸反映了末筍期與其他筍期（初筍期、盛筍早期、盛筍后期）竹筍品質的差異，蛋白質指標分布在第一軸負方向，且與初筍期、盛筍早期、盛筍后期呈不同程度的正相關關系，而與末筍期呈負相關關系。單寧等品質指標分布在第一軸的正軸方向，與末筍期呈較好的正相關關系（圖中箭頭相對方向表示相關性，當相應指標箭頭指向的方向相近時，就有較大的正相關性，方向相反則意味著負相關）。第二軸反映的是挖筍對竹筍品質的影響，由圖中可知，發筍量、單寧含量與釆筍行為呈正相關，且挖筍行為對盛筍后期發筍量及末筍期單寧含量影響較大。粗脂肪含量與釆筍呈負相關，含水率和蛋白質含量與挖筍的相關性較弱。

3 討 論

3.1 挖筍對出筍成竹規律的影響

竹筍采挖是毛竹林經營過程中的重要活動，挖筍延長了春筍的筍期，并使發筍時間提早（表2），這可能與樣地冬筍的采挖有關。冬筍全面翻挖，一方面刺激了休眠筍芽萌發，另一方面，挖筍過程中對土壤的翻動，使表層土壤疏松，有利于春筍露頭。此外，溫度是毛竹種群竹筍—幼竹生長過程中的重要影響因子，陶芳明等[16] 發現，竹筍生長的初期和中期，溫度對竹筍生長起主要作用，出筍的最適溫度在10 ～ 20 ℃之間，而翻挖后粗糙的土壤表面可以更好地吸收太陽輻射，增加了土壤溫度，進而使出筍期提前。

毛竹林內竹筍到幼竹的生長過程對光的需求較少，光合能力弱。研究表明[5]，毛竹幼竹在林下迅速進入林冠層并不依賴較強的光照強度。竹筍生長所需的養分、水分、碳水化合物等資源主要依賴母株通過克隆整合進行供給。克隆整合是克隆植物特有的生理特性，資源可通過地下鞭根系統在母株與子代分株之間進行傳輸與分享[17-18]。因此，竹類植物竹筍的成活對母株有較強的依賴性，在無母株條件下，根莖發筍率低、成活率低[19]。盡管母竹對幼竹生長的貢獻很大，但其所能提供的資源畢竟有限，同一地下鞭根系統上相鄰的竹筍常因爭奪資源產生競爭，而在爭奪資源的過程中，個體在同種密度高的生境中往往會出現較低的存活率。未挖筍樣地由于冬季未采挖冬筍，隨著春季地溫上升，大量冬筍打破休眠，于初筍期至盛筍早期出土生長，竹筍在空間和資源上均存在較大的競爭，退筍率較高，至發筍中后期，由于竹林要維持前期竹筍的生長，造成營養缺乏，發筍受到抑制。而挖筍樣地合理留養竹筍，在各筍期將易出現退筍現象的竹筍（生長較弱或競爭資源空間時處于不利地位）及時采挖，既能產生一定的經濟效益，也可為中后期發筍提供空間，節約資源，緩解競爭的壓力，因而退筍率顯著低于未挖筍樣地（P ＜ 0.05），而發筍數顯著高于未挖筍樣地（P ＜ 0.05）（表2）。

研究表明[20-22]，當受到人工采挖等外部干擾時，克隆植物能夠促進受到外部干擾的分株再生長或生成新分株，這是其應對外部干擾的重要補償性機制。竹林中的挖筍行為可能通過毛竹克隆整合或克隆儲存等特性，調整竹筍間的非對稱競爭關系來實現補償性機制，進而增加整個筍期的發筍量。挖筍的樣地對冬筍進行全面采挖，初筍期至盛筍早期竹筍萌發的數量顯著低于未挖筍樣地（P ＜ 0.05）。但冬春筍的采收并未減少林地整個筍期的發筍數量，反而增強了春筍發筍能力，這可能與竹筍間的非對稱性競爭有關。研究發現[23]，當較大植物占有較多競爭資源時，較小植物最終只能因資源枯竭而死亡和被淘汰。而隨著林地部分竹筍的采挖，打破了竹筍間原有的非對稱競爭，促進一些原本受抑制的筍芽開始萌發生長。董文淵等[24] 也發現，連年采筍的筇竹林中，初期和盛期較大的竹筍被采挖后? 無性系會自我調控，地下鞭上萌發更多筍芽以彌補子代分株數量的不足，通過增加竹筍數量，增加種群在生境中駐留和延續的機會。

3.2 挖筍對新分株胸徑的影響

初筍期和盛筍早期毛竹新分株的胸徑較大，而隨著筍期的推移胸徑逐漸減小。這主要是因為，初筍期和盛筍早期的竹筍搶先利用母株貯藏的物質，從而妨礙了中后期竹筍生長。李睿等[25] 發現，較早同生群的竹筍及其長成的幼竹總比較晚同生群的長得高、長得壯。而較為粗壯的竹筍往往競爭力強，易成活，在各筍期的退筍率也小（圖2）。

有研究認為[24]，整個筍期一直挖筍，當挖筍強度超過了竹林無性系恢復能力時，會破壞竹林的生長平衡，竹林營養匱乏及踐踏等原因會使幼竹基徑變小。而本研究發現，未挖筍和挖筍樣地內新分株的平均胸徑間沒有顯著差異（P ＞ 0.05）。挖筍樣地會盡量留養較為粗壯的大筍，將生長較弱或空間位置不利的竹筍挖除，而未挖筍樣地在竹筍的自疏過程中，同樣采取優先存活大筍的生存策略，維持種群的競爭力。兩種處理的林地均長期施肥與養護，為大筍的生長提供養分，保證其在競爭中存活，并成長為新的分株。從圖2 可知，未挖筍和挖筍樣地的存活分株中地徑大于12 cm 的分株占了很大的比例，說明挖筍并沒有對新分株胸徑大小產生不利影響。且研究表明[26]，毛竹1 年新分株的胸徑大小很大程度上受到上一代3 年分株的影響，挖筍、未挖筍樣地竹林結構、毛竹平均胸徑以及立地條件均是一致的，兩種處理樣地中新分株的胸徑也可能受到上一世代分株胸徑的影響，因而沒有顯著差異。

3.3 挖筍對林地退筍的影響

處于競爭狀態的植物之間有時會存在資源的不均衡分配，優勢植物會占據大于相應生物量比例的資源，這種現象稱為非對稱競爭[27-29]。竹類植物的母株在發筍成竹階段向竹筍和幼竹提供生長所需資源時，也存在非對稱性競爭的情況。研究發現[7]，較大個體的竹筍在資源爭奪中更占優勢，并會抑制小筍存活，母株會將資源優先供給那些生長健壯的竹筍，竹筍的個體越強壯，競爭能力越強，獲得的資源越多。這使得未挖筍樣地的竹筍在徑級Ⅳ級以上有較高死亡率，而徑級在Ⅰ～Ⅲ級的竹筍死亡率相對較小（圖2），而挖筍樣地，雖然僅保留了壯筍（Ⅰ～Ⅳ級），但竹筍的死亡率隨著地徑的增大逐漸從46.2% 下降至21.6%。

未挖筍樣地在盛筍早期以后，徑級較小（Ⅳ級以上）的竹筍逐漸增多，這類竹筍缺乏資源競爭優勢，若不進行采挖，死亡率極高，造成資源浪費。且竹筍在退筍過程中仍會從母株獲取并消耗一部分資源，這使得出筍中、后期樣地資源更加匱乏，因此，未挖筍樣地在盛筍后期和末筍期發筍以徑級在Ⅳ級以上的竹筍為主，且大量死亡，而Ⅰ—Ⅲ級的壯筍發筍數量相對較少。

3.4 挖筍對竹筍品質的影響

竹筍的營養價值與其蛋白質的含量密切相關，挖筍和未挖筍樣地內竹筍蛋白質含量在初筍期、盛筍早期、盛筍后期差異不大，但末筍期均顯著降低。這是因為竹筍進入春季生長階段，竹筍快速生長需要大量消耗毛竹體內蛋白質，從而使得末筍期竹筍中蛋白質貯存量減少[13]。與未挖筍林地相比，挖筍后林地竹筍單寧含量在初筍期至盛筍后期沒有顯著差異，兩種林地均表現為初筍期竹筍單寧含量相對較低，隨著發筍時間推移，盛筍后期和末筍期單寧含量顯著提高，但末筍期采筍樣地的竹筍單寧含量顯著高于未采筍樣地。末筍期單寧含量的增加可能與毛竹林地大量發筍導致末筍期林地土壤速效養分的虧缺有關，土壤養分的不足往往會導致植物體內單寧的含量增加。研究表明[15，30]，溫度、光照強度、土壤含量水、土壤肥力等外界因子都會影響植物體內單寧的積累。邱永華等[31] 也認為，當土壤相對貧瘠時，林地雖然通過施肥補充養分，但迅速增加的竹筍產量會使土壤養分虧缺狀態加劇，從而使竹筍中的單寧含量反而提高。總的來說，挖筍不僅有效地提高發筍數量，竹筍釆挖量的增大也并未對竹筍品質產生明顯不良影響。

4 結 論

采挖毛竹林內的竹筍可緩解竹筍間非對稱性的資源競爭，增加竹筍產量，提高竹筍存活率，有效避免大量退筍導致的毛竹母株營養消耗，使母株資源得以重新整合，減少資源浪費。而竹筍釆挖量的增加也并未對新分株胸徑及竹筍品質產生不良影響。因此，合理的竹筍采挖有利于毛竹筍用林的高效生產及可持續性經營。本研究主要探討的是竹筍采挖對毛竹林竹筍產量以及營養品質的影響。但隨著生活水平的提高，人們在毛竹林經營過程中，不僅會關注竹林生產力的提升[32-33]，對竹筍的食味品質（甜味、酸味、鮮味、芳香味等）也會有所要求。因此，食味品質是體現竹筍經濟價值的重要因素，而竹筍食味是多種物質綜合作用的結果，采挖過程中是否會對竹筍內與食味相關的芳香類氨基酸及醛、醇類等物質的含量產生影響，尚待后續的研究，以期實現毛竹林高產高效生態經營模式的進一步優化。

參考文獻：

[1] COMITA L S， HUBBELL S P. Local neighborhood and species’shade tolerance influence survival in a diverse seedling bank[J].Ecology，2009，90（2）：328-334.

[2] 賈美玉， 張金屯， 梁鈺， 等. 五臺山林線過渡帶華北落葉松格局過程和尺度分析[J]. 生態學報，2018，38（20）：7499-7507.

JIA M Y， ZHANG J T， LIANG Y， et al. Spatial pattern andprocess of Larix principis-rupprechti at different scales in thetimberline ecotone of Mount Wutai[J]. Acta Ecologica Sinica，2018，38（20）：7499-7507.

[3] 黃彪， 劉廣路， 范少輝， 等. 毛竹向闊葉林擴展過程細根可塑性變化[J]. 中南林業科技大學學報，2021，41（10）：11-19.

HUANG B， LIU G L， FAN S H， et al. The plasiticity of fine rootPhyllostachys edulis expanding into adjacent broadleaved forest[J].Journal of Central South University of Forestry amp; Technology，2021，41（10）：11-19.

[4] 蔡宗明， 鄧智文， 李東寶， 等. 帶狀采伐對毛竹生物量和根系非結構性碳水化合物的影響[J]. 中南林業科技大學學報，2023，43（4）：33-42.

CAI Z M， DENG Z W， LI D B， et al. Effects of strip logging on thebiomass and root non-structural carbohydrates of Phyllostachysedulis[J]. Journal of Central South University of Forestry amp;Technology，2023，43（4）：33-42.

[5] 方飛燕， 白尚斌， 周國模， 等. 光照強度對毛竹竹筍到幼竹生長過程的影響[J]. 東北林業大學學報，2012（3）：11-13，27.

FANG F Y， BAI S B， ZHOU G M， et al. Effects of shading ongrowth of Phyllostaehys pubescens[J]. Journal of Northeast ForestryUniversity，2012（3）：11-13，27.

[6] 李睿， 鐘章成，M J A 維爾格. 中國亞熱帶高大竹類植物毛竹竹筍克隆生長的密度調節[J]. 植物生態學報，1997，21（1）：9-18.

LI R， ZHONG Z C， WERGER M J A. Density regulation of theclonal growth of new shoots in the giant bamboo Phyllostachyspubescens in subtropical China[J]. Acta Phytoecologica Sinica，1997，21（1）：9-18.

[7] 鄒躍國. 毛竹林種群結構特征變化的海拔效應研究[J]. 竹子研究匯刊，2012，31（4）：17-21.

ZHOU Y G. Influence of altitude on population structurecharacteristic of Moso Bamboo forest[J]. Journal of BambooResearch，2012，31（4）：17-21.

[8] 鄭祥乾， 崔永忠， 陳凌娜， 等. 屏邊空竹四季出筍及幼竹生長發育規律研究[J]. 林業科學研究，2018，31（5）：131-136.

ZHENG X Q， CUI Y Z， CHEN L N， et al. Study on bambooshooting and shoot growth of Cephalostachyum pingbianense[J].Forest Research，2018，31（5）：131-136.

[9] 張瑋， 林振清， 楊前宇， 等. 椽竹出筍與幼竹生長規律[J]. 浙江農林大學學報，2015，32（3）：478-482.

ZHANG W， LIN Z Q， YANG Q Y， et al. Bambusa textilis var.fasca shooting and young bamboo growth rhythm[J]. Journal ofZhejiang A amp; F University，2015，32（3）：478-482.

[10] 張國學， 輝朝茂， 石明， 等. 甜龍竹發筍和退筍規律稈高生長特性研究[J]. 世界竹藤通訊，2020，18（1）：30-33.

ZHANG G X， HUI C M， SHI M， et al. Shoot germinationand degradation law and culm growth characteristics ofDendrocalamus brandisii[J]. World Bamboo and Rattan，2020，18（1）：30-33.

[11] 李露雙， 陳凌娜， 李斌， 等. 屏邊空竹四季出筍動態及其與氣象因子相關性分析[J]. 西北林學院學報，2021，36（4）：66-72.

LI L S， CHEN L N， LI B， et al. Seasonal variations of bambooshooting and its correlation with meteorological factors onCephalostachyum pingbianense[J]. Journal of Northwest ForestryUniversity，2021，36（4）：66-72.

[12] 劉慶， 鐘章成. 斑苦竹無性系種群的數量統計[J]. 西南師范大學學報（ 自然科學版），1995，20（2）：176-182.

LIU Q， ZHONG Z C. The demography of pleioblastus maculateclone population[J]. Journal of Southwest China NormalUniversity （Natural Science），1995，20（2）：176-182.

[13] 劉耀榮. 毛竹筍期的營養動態[J]. 林業科學研究，1990，3（4）：363-367.

LIU Y R. Nutrient dynamics in the shoot stage of MaoBamboo[J]. Forest Research，1990，3（4）：363-367.

[14] 徐森， 谷瑞， 陳雙林， 等. 毛竹春筍營養、食味和基于重金屬的安全品質的海拔效應[J]. 江西農業大學學報，2021，43（1）：144-152.

XU S， GU R， CHEN S L， et al. Altitude effect on nutrition，tasteand safety quality based on heavy metals of Phyllostachys edulisspring shoots[J]. Acta Agriculturae Universitatis Jiangxiensis，2021，43（1）：144-152.

[15] HANSEN A H， JONASSON S， Michelsen A， et al. Long-termexperimental warming， shading and nutrient addition affect theconcentration of phenolic compounds in arctic-alpine deciduousand evergreen dwarf shrubs[J]. Oecologia，2006，147（1）：1-11.

[16] 陶芳明， 洪傳禮. 毛竹林火燒跡地竹筍生長與氣候因子的關系[J]. 竹子研究匯刊，1990，9（2）：65-71.

TAO F M， HONG C L. Relationship between bamboo shootgrowth and climate factor in fire-slash of bamboo forest[J].Journal of Bamboo Research，1990，9（2）：65-71.

[17] YU F H， WANG N， ALPERT P， et al. Physiological integrationin an introduced， invasive plant increases its spread intoexperimental communities and modifies their structure[J].American Journal of Botany，2009，96（11）：1983-1989.

[18] YU F H， WANG N， HE W M， et al. Effects of clonal integrationon species composition and biomass of sand dune communities[J].Journal of Arid Environment，2010，74（6）：632-637.

[19] 王瓊， 蘇智先， 雷濘菲， 等. 慈竹母株大小對克隆生長的影響[J]. 植物生態學報，2005，29（1）：116-121.

WANG Q， SU Z X， LEI N F， et al. Effects of parent ramet sizeon clonal growth in neosinocalam affinis[J]. Acta PhytoecologicaSinica，2005，29（1）：116-121.

[20] QIN R R，WANG X L. Effects of crown height on thecompensatory growth of Italian ryegrass based on combined effectsof stored organic matter and cytokinin[J]. Grassland Science，2020，66（1）：29-39.

[21] ZHANG D M， ZHAO W Z， LUO W C. Effect of the populationdensity on belowground bud bank of a rhizomatous clonal plantLeymus secalinus in Mu Us sandy land[J]. Journal of Plant Research，2019，132（1） ：69-80.

[22] 葉學華， 薛建國， 謝秀芳， 等. 外部干擾對根莖型克隆植物甘草自然種群植株生長及主要藥用成分含量的影響[J]. 植物生態學報，2020，44（9）：951-961.

YE X H， XUE J G， XIE X F， et al. Effects of differentdisturbances on plant growth and content of main medicinalingredients of rhizomatous clonal plant Glycyrrhiza uralensisin a natural population[J]. Chinese Journal of Plant Ecology，2020，44（9）：951-961.

[23] SCHWINNING S， WWINER J. Mechanisms determiningthe degree of size asymmetry in competition among plants[J].Oecologia，1998，113（4）：447-455.

[24] 董文淵， 黃寶龍， 謝澤軒， 等. 密度調節與輪閑制采筍對筇竹林竹筍- 幼竹生長的影響[J]. 林業科學，2002，38（5）：78-82.

DONG W Y， HUANG B L， XIE Z X， et al. The Effects of densityregulation and rotational harvesting shoot on the shoot and youngbamboo growth of qiongzhuea tumidinoda forest[J]. ScientiaSilvae Sinicae，2002，38（5）：78-82.

[25] 李睿， 鐘章成，M J A 維爾格. 毛竹（Phyllostachys pubescens）竹筍群動態的研究[J]. 植物生態學報，1997，21（1）：53-59.

LI R， ZHONG Z C， WERGER M J A. Studies on the dynamicsof the bamboo shoots in phyllostachys pubescens[J]. ActaPhytoecologica Sinica，1997，21（1）：53-59.

[26] 王意錕， 金愛武， 朱強根， 等. 施肥對毛竹種群不同年齡分株間胸徑大小關系的影響[J]. 植物生態學報，2014，38（3）：289-297.

WANG Y K， JIN A W， ZHU Q G， et al. Effects of fertilization onthe relations of diameter at breast height between different-agedramets of Phyllostachys edulis population[J]. Chinese Journal ofPlant Ecology，2014，38（3）：289-297.

[27] VELAZQUEZ E， WIEGAND T. Competition for light andpersistence of rare light-demanding species within tree-fall gapsin a moist tropical forest[J]. Ecology，2020，101（7）：e03034.

[28] PICARD N. Asymmetric competition can shape the sizedistribution of trees in a natural tropical forest[J]. Forest Science，2019，65（5）：562-569.

[29] PRETZSCH H， BRAVO-OVIEDO A， HILMERS T， et al. Withincreasing site quality asymmetric competition and mortalityreduces Scots pine （Pinus sylvestris L.） stand structuring acrossEurope[J]. Forest Ecology and Management，2022，520：120365.

[30] ANURAGA M， DUARSA P， LOVETT J V， et al. Soil moisture andtemperature affect condensed tannin concentrations and growth inLotus corniculatus and Lotus pedunculatus[J]. Australian Journalof Agricultural Research，1995，44（7）：1667-1681.

[31] 邱永華， 金愛武， 張四海， 等. 不同施肥方式對竹筍品質的影響[J]. 竹子學報，2017，36（1）：41-48.

QIU Y H， JIN A W， ZHANG S H， et al. Effect of differentfertilization methods on quality of bamboo shoots[J]. Journal ofBamboo Research，2017，36（1）：41-48.

[32] 羅天磊， 范少輝， 劉廣路， 等. 閩中毛竹立地質量評價及分級研究[J]. 中南林業科技大學學報，2020，40（7）：48-57.

LUO T L， FAN S H， LIU G L， et al. Study on site quality evaluationand classification of Phyllostachys edulis in central Fujian[J].Journal of Central South University of Forestry amp; Technology，2020，40（7）：48-57.

[33] 孫慧， 曹志華， 吳中能， 等. 覆蓋對毛竹林土壤養分、酶活性及微生物群落的影響[J]. 經濟林研究，2023，41（2）：223-233.

SUN H， CAO Z H， WU Z N， et al. Effects of mulching onsoil nutrients， enzyme activities and microbial community inPhyllostachys edulis forest[J]. Non-wood Forest Research，2023，41（2）：223-233.

[ 本文編校：趙 坤]

基金項目：浙江省科技廳領雁計劃科技攻關項目（2022C03119）；麗水市公益性技術應用研究項目（2022GYX19）；中央財政林業科技推廣示范資金項目（〔2019〕TS16）。