代 麗，張 喜，吳良如

（1.貴州大學(xué)，貴州貴陽550025；2.貴州省林業(yè)科學(xué)研究院，貴州貴陽550005；3.國家林業(yè)和草原局叢生竹工程技術(shù)中心，云南昆明650224；4.國家林業(yè)和草原局竹子研究開發(fā)中心，浙江杭州310012）

金佛山方竹（Chimonobambusautilis）屬禾本科竹亞科寒竹屬植物，天然分布于重慶、四川、貴州和云南交界區(qū)域，是我國西南地區(qū)特有的筍用竹種之一［1］。黔北渝南的大婁山區(qū)有地球上至今保存最完整、原生性最強和成片面積最大的方竹天然林。本文系統(tǒng)分析了金佛山方竹近50余a來的研發(fā)文獻，總結(jié)已有成果，把握研究的最新進展和未來動態(tài)，旨在為方竹產(chǎn)業(yè)發(fā)展及現(xiàn)代山地特色農(nóng)業(yè)建設(shè)提供借鑒和啟示。

1 研究方法

（1）文獻收集。通過科睿唯安（WOS）和中國知網(wǎng)（CNKI）系統(tǒng)，輸入“金佛山方竹”關(guān)鍵詞搜索并下載到1964－2020年國內(nèi)外公開發(fā)行刊物的中文文獻130余篇、外文文獻近10篇，以及14篇碩士論文。

（2）文獻分析。應(yīng)用中國知網(wǎng)（CNKI）系統(tǒng)中的文獻導(dǎo)出功能羅列文獻清單，進一步分析文獻內(nèi)容及綜述熱點；可視化功能分析文獻量總體變化趨勢、文獻互引和關(guān)鍵詞共現(xiàn)結(jié)構(gòu)，以及學(xué)科、期刊來源、基金類型、作者及機構(gòu)等信息。

2 結(jié)果分析

2.1 文獻計量分析

2.1.1 文獻量動態(tài)變化 金佛山方竹文獻量的時間變化動態(tài)（圖1）發(fā)現(xiàn)：研究文獻、參考文獻及引證文獻整體呈上升趨勢，進入21世紀(jì)后呈波動性快速增長態(tài)勢。國內(nèi)首篇有關(guān)金佛山方竹的研究文獻是1964年發(fā)表在《林業(yè)科學(xué)》上的“國產(chǎn)33種竹材制漿應(yīng)用上纖維形態(tài)結(jié)構(gòu)的研究”，至1985年的20余a間文獻量占總量的4.72%，至1995年間的30余a間文獻量占總量的10.24%，至2005年的40余a間文獻量占總量的30.47%，至2015年的50余a間文獻量占總量的76.56%，2016－2020年的5 a間文獻量占總量的23.62%、碩士論文量占總量的85.71%，文獻量達(dá)新的高度。參考文獻始于1938年前中央大學(xué)楊銜晉教授在重慶金佛山古佛洞采集到此竹標(biāo)本命名起，引證文獻始于1983年。20世紀(jì)80年代后，有關(guān)金佛山方竹研究的參考文獻、引證文獻和研究文獻量增速加大，呈參考文獻量或引證文獻量大于研究文獻量的變化趨勢。

圖1 文獻量（左）和文獻量、參考文獻量及引證文獻量相互關(guān)系（右）的動態(tài)變化Fig.1 Dynamic changes of amount of literature（left）and the amount relationship between literature，reference literature and cited literature（right）

2.1.2 文獻來源及時間動態(tài) 公開發(fā)表文獻的統(tǒng)計分析發(fā)現(xiàn)（圖2）：金佛山方竹近40 a研究成果主要發(fā)表在20類科技期刊中，排名前5位的是《世界竹藤通訊》（13.16%）、《貴州林業(yè)科技》（11.40%）、《竹子學(xué)報》（9.65%）、《林業(yè)科技通訊》（6.14%）、《四川林業(yè)科技》《現(xiàn)代農(nóng)業(yè)科技》《種子》（3.51%）；近10 a研究成果主要發(fā)表在10類科技期刊中，排名前4位的期刊序列未變、文獻量比例增大，第5位增加了《貴州大學(xué)學(xué)報》（自然科學(xué)版）。表明文獻量向主流期刊匯聚，權(quán)威期刊文獻量仍然較少。文獻的學(xué)科分類中，林業(yè)占66.46%、園藝占7.45%、農(nóng)業(yè)經(jīng)濟占7.45%、輕（手）工業(yè)或生物學(xué)占3.73%、農(nóng)業(yè)基礎(chǔ)科學(xué)占3.11%，表明林業(yè)行業(yè)對金佛山方竹的關(guān)注度較高。

圖2 文獻來源及40年（左）和10年（右）動態(tài)變化Fig.2 Literature sources and dynamic changes in 40 years（left）and 10 years（right）

2.1.3 論文作者及機構(gòu)分布 公開發(fā)表文獻的統(tǒng)計分析發(fā)現(xiàn)（圖3）：論文作者文獻量呈4個梯隊分布：①文獻量＞3.8%。有張喜（6.98%）、張佐玉（5.43%）和茍光前（3.88%）等，占總量的7.50%；②文獻量在3.10%左右。有江龍和易思榮等8位，占總量的20.00%；③文獻量在2.33%左右。有董文淵和王福升等18位，占總量的45.00%；④文獻量≤1.55%。有劉國華和杜凌等11位，占作者總量的27.50%。論文機構(gòu)文獻量也呈4個梯隊分布：①文獻量＞10%。有貴州大學(xué)（14.79%）和貴州省林業(yè)科學(xué)研究院（11.27%）2家單位，占總量的6.45%；②文獻量5% ～10%。有貴州省桐梓縣林業(yè)局（9.86%）、重慶市南川區(qū)林業(yè)局（9.15%）、南京林業(yè)大學(xué)（5.63%）等5家單位，占總量的16.13%；③文獻量1% ～5%。有西南林業(yè)大學(xué)（4.93%）、貴州省正安縣林業(yè)局（3.52%）、重慶市藥物種植研究所（2.82%）等10家單位，占總量的32.26%；④文獻量＜1%。有云南省林業(yè)科學(xué)研究院和浙江大學(xué)等14家單位，占總量的45.16%。文獻量按?。ㄊ校┘壟判?，貴州?。?3.38%）、重慶市（14.79%）和云南省（7.75%）的文獻量較高，占總量的85.92%。

圖3 文獻作者（左）及機構(gòu)（右）排序Fig.3 Literature ranking of author（left）and institution（right）

2.1.4 基金類型及主題詞分布 按省級及以上資助課題統(tǒng)計（圖4）：金佛山方竹研究主要受國家級科技支撐計劃及國家自然科學(xué)基金資助，占總量的61.76%；省級科技支撐計劃及自然科學(xué)基金資助度較高的貴州省占總量的26.47%，其次是四川省占總量的5.88%，重慶市占總量的2.94%。主題詞頻率可劃分為4類：①頻率＞3.0%。包括金佛山方竹（59.17%）、方竹筍（7.69%）、竹產(chǎn)業(yè)發(fā)展（5.92%）和低產(chǎn)林改造（3.55%），占主題詞總量的21.05%；②頻率2.0% ～3.0%。包括筍用林、母竹栽培、幼竹生長或營養(yǎng)成分等，占主題詞總量的21.05%；③頻率1.0%～2.0%。包括外生菌根真菌、栽培管理或施肥試驗等，占主題詞總量的47.37%；④頻率＜1.0%。包括選擇試驗和信息動態(tài)等，占主題詞總量的10.53%。

圖4 課題資助類型及文獻主題詞分布Fig.4 Types of project funding and the subject term distributions of literature

2.1.5 金佛山方竹和主題詞（左）及次主題詞（右）共現(xiàn)頻率 共現(xiàn)頻率矩陣中，金佛山方竹和其它主題詞及次主題詞關(guān)系發(fā)現(xiàn)（圖5）：主題詞排序為方竹筍（16.33%），產(chǎn)業(yè)發(fā)展或幼竹生長（8.16%），母竹栽培、筍用林或采收加工（6.12%），無性系種群、生物學(xué)特性、栽培管理、筍用混交林、施肥試驗、低產(chǎn)林改造或礦質(zhì)元素（4.08%）；次主題詞排序中，方竹筍（50.00%），喬木層（20.00%），海拔高度、形態(tài)學(xué)鑒定或化合物（10.00%）。

圖5 金佛山方竹和主題詞（左）及次主題詞（右）共現(xiàn)頻率Fig.5 Co-occurrence frequency between Chimonobambusa utilis and subject words（left）and sub-subject words（right）

2.2 主要研究領(lǐng)域

2.2.1 竹林特性研究 金佛山方竹是中亞熱帶山地常綠落葉闊葉混交天然林的重要組成物種。如貴州省桐梓縣柏菁喀斯特臺原包含金佛山方竹的森林類型［2］有厚皮栲櫻－金佛山方竹林、厚皮栲野茉莉－金佛山方竹林、厚皮栲稠李－金佛山方竹林、美脈花楸嶺南石櫟－金佛山方竹林、巴東櫟－金佛山方竹林、珙桐［3］－金佛山方竹林、大葉楊－金佛山方竹林及方竹純林8個類型。用地形地貌、母巖土壤及竹林結(jié)構(gòu)等要素因子可將桐梓縣［4］金佛山方竹林聚類劃分為10個類型，竹林平均胸徑、密度及生物量變幅為1.06～3.58 cm、1.69×104～1.31×105株·hm－2和12.1～193.3 t·hm－2，竹林［5］平均胸徑和海拔高度、喬木層植物種數(shù)、平均胸徑及高度的正相關(guān)性顯著、和竹林密度的負(fù)相關(guān)性顯著。隨著金佛山方竹［6］的自然擴散，喬木層植物種類及平均高度變化微弱、蓋度升降不定，灌木層植物種類及平均蓋度降低、平均高度先降后升，草本層植物種類、平均高度及平均蓋度明顯降低，有竹樣地幼苗種豐富度、數(shù)量、幼苗保存率及年際波動頻率顯著低于無竹樣地，幼苗葉生物量及葉根比值顯著高于無竹樣地。高密度金佛山方竹林［7］降低了亮葉水青岡種群數(shù)量及平均胸徑，提高了平均高度。

生活周期60 a左右，以小面積（＜1 hm2）多批次開花為主［8］。除遺傳因素外，干旱及土壤質(zhì)量低劣、林分及立竹年齡結(jié)構(gòu)失調(diào)也是重要的誘導(dǎo)因素［9］?，F(xiàn)存質(zhì)量較高且原生性較強的貴州省桐梓縣柏菁喀斯特臺原［2，9］金佛山方竹林是20世紀(jì)30－40年代開花后經(jīng)種子自然萌發(fā)或?qū)嵣缛斯ひ浦渤闪值模?0］，林地連續(xù)撫育管理5－10 a即有較高的竹筍產(chǎn)量。

竹筍不發(fā)于春而茂于秋，翌年春季完成新竹生長。黔北［11－12］地區(qū)出筍始于9月上中旬，結(jié)束于10月中旬，歷時30余d；新竹高生長始于9月上中旬，止于11月中下旬，歷時60－70 d；地徑生長和高生長同步，止于10月上旬，歷時25－30 d。出筍初日和末日隨海拔高度降低而延后、出筍日數(shù)降低［13］。竹筍重量及凈筍／毛筍重量比值［14］隨地徑及高度增大而增加，稈重［15］隨胸徑及高度增大而增加。翌年［16］新竹枝及葉自冠層上部至下部開始生長、停止時間相同；不同冠層枝長、節(jié)間長、葉長和葉寬的平均值呈上（中）層＞下層的趨勢。

立竹空間分布格局［17］隨研究尺度增大由聚集加隨機性分布向隨機性分布轉(zhuǎn)化。不同年齡立竹［18］葉長、葉寬、葉面積、葉干質(zhì)量和比葉面積5指標(biāo)隨年齡增大而增加，3年竹的相應(yīng)指標(biāo)顯著高于1年竹；海拔高度和喬木冠層［19］對立竹及葉片形態(tài)有明顯影響。密度等級［20］和年齡株數(shù)結(jié)構(gòu)配比試驗發(fā)現(xiàn)密度等級、密度等級和年齡結(jié)構(gòu)交互作用對立竹生物量影響不顯著，年齡結(jié)構(gòu)對立竹生物量影響顯著，立竹生物量在2年竹和3年竹間差異不顯著、在2－3年竹和1年竹及4年竹間差異顯著。

2.2.2 培育技術(shù)研究 金佛山方竹林生長受氣候、地貌及土壤特性的綜合影響。建立的出筍期8－9月溫度、降水及日照量和經(jīng)度、緯度及海拔高度間的回歸模型［21］推算貴州省桐梓縣金佛山方竹氣候適宜區(qū)占國土面積的42.76%，進一步［22］發(fā)現(xiàn)適宜區(qū)、中等適宜區(qū)和勉強適宜區(qū)面積分別為1.45×104hm2、3.96×104hm2和6.85×104hm2，占國土面積的38.40%。海拔高度［23］1 600～2 000 m為適宜區(qū)，＜1 600 m和＞2 000 m為一般區(qū)。

自然落地萌發(fā)的種子當(dāng)年苗高可達(dá)20 cm以上，采摘的種子有半月后熟期，宜陰涼通風(fēng)環(huán)境中短期鮮藏。相同溫度梯度［24］下不同種源種子萌發(fā)率差異顯著，最高萌發(fā)率溫度不同，18～22℃種子萌發(fā)率高且長勢好，22℃以上易誘發(fā)幼苗白化。種子［25］經(jīng)50 mg·L－1萘乙酸、或100 mg·L－1吲哚丁酸或300 mg·L－1吲哚丁酸12 h處理后，在15%聚乙二醇脅迫下的發(fā)芽率、苗高和根冠比等綜合評價值較高。種子苗［26］高度、地徑、株／叢、筍芽數(shù)量、根幅及總鮮重間相關(guān)顯著，據(jù)此以苗高和地徑為主要指標(biāo)可劃分1年苗或2年留床苗的3級苗木質(zhì)量標(biāo)準(zhǔn)。

竹苗包括實生苗［23，26］和母竹苗［23］2種類型。海拔高度［5］和單株母竹造林成活率的正相關(guān)性顯著，呈截稈＞整株及1年竹＞2年竹＞3年竹的趨勢。母竹造林［27］3年后的叢成竹株數(shù)、地徑、高度及走鞭長度呈3措施（松土埋青＋除草＋復(fù)合肥0.5 kg·穴－1）＞2措施（松土埋青＋復(fù)合肥0.5 kg·穴－1，松土埋青＋除草）＞1措施（松土埋青）＞對照的變化趨勢，3措施的竹林郁閉度值第3年可達(dá)0.5、第5年可達(dá)0.8。

低產(chǎn)筍用林改造關(guān)鍵技術(shù)［28］包括調(diào)整竹木結(jié)構(gòu)、竹林密度及年齡結(jié)構(gòu)，以及留筍蓄竹、林地覆蓋和間伐物歸林等。密度調(diào)整是主要的經(jīng)營措施之一、密度受林分平均胸徑的影響，平均胸徑0.86 cm的天然林適宜經(jīng)營密度為4.5×104株·hm－2，平均地徑［29］1.41 cm的人工林適宜經(jīng)營密度為2.0×104株·hm－2。竹木結(jié)構(gòu)調(diào)整是重要的經(jīng)營措施之一，平均胸徑［30］1.06～1.51 cm的天然闊葉混交林內(nèi)，新竹平均胸徑和高度及筍竹比例值呈喬木郁閉度值0.5＞0.3＞0.7＞0.1，地上生物量及總生物量生產(chǎn)量值呈喬木郁閉度值0.5＞0.7＞0.3＞0.1的變化趨勢；人工種植［31］的金佛山方竹闊葉混交林（楓香或鵝掌楸）、針葉混交林（柳杉或日本落葉松）和純林郁閉度值達(dá)0.5及以上的時間分別為3 a、4 a和5 a，竹林平均胸徑及高度值呈闊葉混交林＞針葉混交林＞方竹純林的變化趨勢。天然林［32］試驗發(fā)現(xiàn)1措施（施肥、林地覆蓋或密度調(diào)整）平均胸徑提高率8.37%，2措施組合的相應(yīng)值為16.09%，3措施組合的相應(yīng)值為24.07%。

2.2.3 竹筍產(chǎn)量及質(zhì)量研究 竹筍產(chǎn)量受立地、林分和經(jīng)營措施的綜合影響。海拔［33］1 400～2 000 m變幅內(nèi)，竹筍產(chǎn)量隨海拔高度升高而增加，變幅0.50～1.97 t·hm－2·a－1。金佛山方竹天然林中，平均胸徑［28］0.86 cm及密度4.5×104株·hm－2的林分連續(xù)4 a改造，平均竹筍產(chǎn)量達(dá)到3.74 t·hm－2·a－1；平均胸徑［20］2.01 cm及密度2.25×104株·hm－2的林分（1年竹∶2年竹∶3年竹∶4年竹＝1∶1∶1∶0）竹筍產(chǎn)量達(dá)4.06 t·hm－·a－1；竹筍產(chǎn)量［30］及日均出筍量和喬木郁閉度值關(guān)系呈0.5＞0.3＞0.7＞0.1的變化趨勢；平均胸徑［34］0.92 cm及密度8 200株·hm－2的林分林地覆蓋糠殼20 cm和10 cm后，當(dāng)年竹筍產(chǎn)量為對照的357.06%、177.79%，達(dá)到4.01 t·hm－2·a－1和2.00 t·hm－2·a－1。金佛山方竹人工林中，純林、針葉混交林和闊葉混交林［31］第5年竹筍產(chǎn)量分別為1.14、1.17和2.00 t·hm－2·a－1，第9年的相應(yīng)值為1.87、1.92和2.83 t·hm－2·a－1；平均地徑1.41 cm及2.0×104株·hm－2的10年生林分，連續(xù)4年密度調(diào)整并輔施0.38 t·hm－2復(fù)合肥，竹筍產(chǎn)量持續(xù)增加、變幅0.98～3.54 t·hm－2·a－1。方竹天然林中，平均胸徑2.20 cm及2.7×104株·hm－2的林分［32］連續(xù)3 a定位觀測發(fā)現(xiàn)1措施、2措施組合和3措施組合竹筍產(chǎn)量為對照的135.08%～163.73%、172.81%～223.17%、250.93%，分別達(dá)到3.1、4.2和4.7 t·hm－2·a－1；平均胸徑2.88 cm及2.4×104株·hm－2的林分中，不同肥種及肥量組合的竹筍增產(chǎn)量和海拔高度變化無關(guān)，相同海拔高度的竹筍增產(chǎn)量隨施肥量增加而增大，為對照的112.75%～127.96%，達(dá)到4.76～6.72 t·hm－2·a－1；平均胸徑1.90 cm林分［35］］中，林分密度（2.3×104～2.5×104株·hm－2）和復(fù)合肥量（1.0×103～1.5×103t·hm－2）組合的竹筍產(chǎn)量較高，達(dá)到4.36 t·hm－2·a－1。

竹筍質(zhì)量受種源及地形的影響，大量元素（僅指氮、磷、鉀、鈣和鎂總量）和微量元素（僅指鐵、錳、銅和鋅總量）、必需氨基酸（僅指賴氨酸、纈氨酸、異亮氨酸、蘇氨酸和苯丙氨酸總量）和呈味氨基酸（僅指谷氨酸、丙氨酸、組氨酸、甘氨酸、天門冬氨酸、精氨酸和脯氨酸總量）含量變化的趨勢性各異。重慶南川區(qū)［36］竹筍含水率、粗蛋白、粗纖維和灰分量較高，貴州桐梓縣［37－38］竹筍粗脂肪含量較高。必需氨基酸量呈貴州金沙縣［39］＞重慶南川區(qū)［36］＞貴州正安縣［13］，呈味氨基酸量呈貴州金沙縣［39］＞重慶南川區(qū)［36］＞貴州正安縣［13］的變化趨勢，大量元素量呈貴州桐梓縣［37－38］＞重慶南川區(qū)［36］的變化趨勢，部分微量元素量呈貴州桐梓縣［37－38］＞重慶南川區(qū)［36］的變化趨勢。山脊、山坡和山腳方竹筍［37］變動系數(shù)（標(biāo)準(zhǔn)差與平均值的%）＞10%的有粗脂肪、粗纖維、可溶性糖和維生素C含量，＜2%的有粗蛋白、氮、鉀和含水量；竹筍［13］必需氨基酸和呈味氨基酸量有低海拔高度（1 300～1 450 m）＜中海拔高度（1 450～1 750 m）＜高海拔高度（1 750～1 900 m）的變化趨勢。

3 討論與建議

3.1 筍用林培育

金佛山方竹已形成較為完整的培育技術(shù)體系，如下方面值得進一步研究：⑴適生區(qū)規(guī)劃及異地引種試驗。金佛山方竹［2，5］在渝南黔北海拔1 300 m以上可見天然林，成片分布于1 600 m以上，2 000 m左右和箭竹呈交錯分布。近年來人工林栽培已至1 000 m左右，也有引種到海拔500 m以上并初步成功的報道［40］。依據(jù)金佛山方竹生物學(xué)主要特性和海拔、氣候、地質(zhì)地貌及土壤等多要素［21－22］規(guī)劃適宜區(qū)并評價生產(chǎn)力，是擴大種植區(qū)域?qū)崿F(xiàn)產(chǎn)業(yè)規(guī)?；闹匾侄?。⑵竹筍產(chǎn)量評估。竹筍產(chǎn)量有年際［28］變化，受林分平均胸徑、密度及經(jīng)營措施的綜合影響。應(yīng)用筍用林定位樣地中長期觀測數(shù)據(jù)，建立竹筍產(chǎn)量和氣候、立地及林分結(jié)構(gòu)等要素間的關(guān)系模型，是科學(xué)測算竹筍產(chǎn)量的重要途徑。⑶竹筍質(zhì)量評估。竹筍質(zhì)量形成受產(chǎn)地區(qū)位、氣候、地質(zhì)地貌、土壤和林分質(zhì)量的綜合影響，土壤重金屬和微量元素有利害之分、有量變及質(zhì)變的變化，建議制定筍用林基地建設(shè)主要環(huán)境元素背景值閥值和竹筍產(chǎn)地溯源機制，適應(yīng)金佛山方竹竹筍產(chǎn)業(yè)高質(zhì)化發(fā)展。⑷新技術(shù)引進。高新技術(shù)引進可加速金佛山方竹林培育提質(zhì)增效，如改進林地覆蓋措施［34］提高竹筍產(chǎn)量，實生苗［41］接種AM真菌后提高苗木生物量、主要酶活性及葉綠素含量，等。

3.2 以竹筍產(chǎn)業(yè)鏈為主導(dǎo)引領(lǐng)全產(chǎn)業(yè)鏈發(fā)展

對照竹產(chǎn)業(yè)生態(tài)系統(tǒng)［42］演化的一般規(guī)律，勾繪出金佛山方竹全產(chǎn)業(yè)鏈圖譜（圖6）。竹筍產(chǎn)業(yè)是主要目標(biāo)，竹筍加工制品系列化及高端化有待加速推進，竹材［43］、竹炭［44］系列產(chǎn)品開發(fā)和竹林［45］綜合利用不可或缺。由此形成全林利用的復(fù)合產(chǎn)業(yè)鏈，進一步鞏固脫貧攻堅成果及助力鄉(xiāng)村振興。

圖6 金佛山方竹全產(chǎn)業(yè)鏈圖譜Fig.6 Whole industry chain map of Chimonobambusa utilis