譚著明，沈?qū)毭?，譚 云，許彥明，馬豐豐，申愛榮，黃 帆，徐清乾

(1. 湖南省林業(yè)科學(xué)院，湖南 長沙 410004；2.湖南省林下特色生物資源培育與利用研究工程中心，湖南 長沙410004；3. 湖南省菌根性食用菌種質(zhì)資源保護(hù)利用中心，湖南 長沙 410004)

發(fā)展人工林是我國森林覆蓋率提升和資源增長的主要途徑。近四十年來，長江流域及以南大型工程造林項(xiàng)目，包括速生豐產(chǎn)用材林、世界銀行貸款造林、長江中上游防護(hù)林建設(shè)、退耕還林、國家儲備林建設(shè)等工程項(xiàng)目中，馬尾松均是重要造林樹種。在林業(yè)由產(chǎn)業(yè)經(jīng)濟(jì)向生態(tài)保護(hù)轉(zhuǎn)型過程中，人工林的傳統(tǒng)利用方式也在發(fā)生改變。以短周期速生材為目標(biāo)的經(jīng)營模式，逐步向用材與生態(tài)效益并舉的多目標(biāo)經(jīng)營模式轉(zhuǎn)變。馬尾松大徑材培育即是這一轉(zhuǎn)變的重要標(biāo)志。培養(yǎng)大徑材林分，有利于充分發(fā)揮地力，順應(yīng)自然生長規(guī)律，發(fā)揮最大生產(chǎn)力，滿足市場對高端木材產(chǎn)品的需求。而且，由大徑級高大個體構(gòu)成的森林環(huán)境，更有利于昆蟲、鳥獸等動物、林下各類植物、菌物棲息、共生，從而發(fā)揮包容復(fù)雜生物的森林系統(tǒng)的良好生態(tài)作用。

對人工馬尾松林生長規(guī)律的研究是經(jīng)營馬尾松林的重要基礎(chǔ)。相關(guān)研究始于二十世紀(jì)六十年代，主要有關(guān)于馬尾松苗期、林分幼林期、中齡期年生長節(jié)律、不同階段和不同生物氣候帶生長規(guī)律的研究[1-5]。時至今日，基于各類生長模型的關(guān)于馬尾松大徑材生長規(guī)律的研究仍備受關(guān)注[6-11]。關(guān)于林分生長的研究，多以林分平均木或徑階平均木代表林分總體的方式展開，即假設(shè)同齡林分各個體間的生長是同步的，模式是一致的。雖然這是人工林分生長規(guī)律研究通行的方式，但畢竟立地與個體營養(yǎng)空間等的異質(zhì)性是客觀存在的，探討林分個體間生長模式的差異，對于深刻理解林分個體發(fā)育過程，指導(dǎo)目標(biāo)樹經(jīng)營具有一定意義。以往，馬尾松研究對象大多是樹齡30年左右的個體[12-14]，50年以上的馬尾松生長過程資料較為少見，這對于生長期較長的馬尾松大徑材生長規(guī)律與成材規(guī)律的認(rèn)識顯然不足。為此，我們針對湖南省60年左右馬尾松林分個體進(jìn)行了研究，試圖解析馬尾松大徑材生長規(guī)律。

1 研究區(qū)概況

樣地1，位于湖南省靖州苗族侗族自治縣排牙山國有林場(109°31′45″E，26°29′53″N)近同齡馬尾松純林內(nèi)。海拔440m。母巖為石灰?guī)r，土壤為山地黃壤，厚度80cm左右。馬尾松系飛籽成林，經(jīng)3次人工間伐。樣地2，位于慈利縣零陽鎮(zhèn)雙溪村女兒溪(111°12′42.836″E，29°25′27.582″N)馬尾松與杉木、楓香等混交林中，海拔230m。母巖為砂礫巖，黃紅壤，土層厚度40cm。坡向西北，坡度15°。馬尾松亦飛籽而成，經(jīng)2次人工間伐。

2 研究方法

2.1 調(diào)查方法

樣地1內(nèi)隨機(jī)設(shè)立3塊標(biāo)準(zhǔn)地，每樣方大小為20 m×30 m。樣地2設(shè)置標(biāo)準(zhǔn)地1個，樣方大小為30 m×40m。對樣方內(nèi)林木個體進(jìn)行每木檢尺、統(tǒng)計(jì)后，按徑級大小，分別選擇大、中、小徑級個體2～3株，伐倒，按常規(guī)方法進(jìn)行樹干解析。整理數(shù)據(jù)時剔除樣地1中年齡過小的個體2株，剔除樣地2中年齡過大或過小個體各1株。用于數(shù)據(jù)分析的個體中，樣地1取13株，樣地2取8株，使得同一樣地內(nèi)用于分析的個體年齡差異均在6年以下。垂直樹干縱軸截取圓盤，鋸平、拋光后，掃描年輪，或設(shè)定比例尺后對圓盤進(jìn)行微距攝影(Cannon EOS 5D)，圖像數(shù)據(jù)以JEPG格式文件存儲，再用萬深 LA-S 系列植物圖像分析系統(tǒng)(年輪分析獨(dú)立版)(杭州萬深檢測科技有限公司，2015 )將各圓盤年輪圖像數(shù)字化，建立Excel文件。隨后，建立各解析木按年齡時序的胸高斷面積生長過程列表。

2.2 數(shù)據(jù)分析

考慮到個體和林分?jǐn)嗝娣e具有較高的穩(wěn)定性、預(yù)估性及在林業(yè)調(diào)查和生產(chǎn)實(shí)踐中的易測定性[15]，本文參考有關(guān)研究結(jié)果[5，16-21]，基于樹干解析獲得的生長過程數(shù)據(jù)和Chapman-Richards函數(shù)及Logistic函數(shù)等分別建立胸高斷面積生長過程模型，采用Levenberg-Marquarat+通用全局優(yōu)化算法(LM-UGO)求解模型參數(shù)。求解過程借助1stOpt軟件( 七維高科有限公司，綜合優(yōu)化軟件包1stOpt.)實(shí)現(xiàn)[22-24]。模型的擬合效果依據(jù)均方根誤差( RMSE) 、相關(guān)系數(shù)平方值(R2) 、 決定系數(shù)(DC ) 指標(biāo)評估。

3 結(jié)果與分析

擬合結(jié)果見表1和表2。從相關(guān)系數(shù)平方值(R2)和決定系數(shù)(DC)看，兩模型擬合精度較高，說明模型均能較好地?cái)M合胸高斷面積生長過程數(shù)據(jù)。

雖然兩模型擬合度較高，但比較分析模型參數(shù)特點(diǎn)，可以發(fā)現(xiàn)二者反映的模型特征存在較大差異，影響到馬尾松樹干胸高斷面積生長過程特點(diǎn)的解析。

Chapman-Richards與Logistic模型參數(shù)中，a值均表示漸近值上限。由表1可知，按照Chapman-Richards模型模擬結(jié)果，所有參與模擬的21個體中，15株屬于LIU等[18]描述的Case1類型，這是多數(shù)樹木遵循的生長模式。3株屬于Case2類型，3株屬于Case8類型。Case1類型中，只有C4、 CLM7、CLM9、CLM10、CML16、CLM18共6個體的a值在28573～279667之間，占40%，且其拐點(diǎn)值各有不同，但均在合理區(qū)間。說明其生長過程在生長加速期、平穩(wěn)期和減速期等模式特征上均不一致，表現(xiàn)為生長的異型性。其余各株的a值均達(dá)4.85×109以上，顯然遠(yuǎn)大于自然界馬尾松樹干胸高斷面積的極限值，且其拐點(diǎn)出現(xiàn)年份遠(yuǎn)遠(yuǎn)超出各個體現(xiàn)實(shí)樹齡范疇，失去合理性。此外，屬于Case2類型的3株個體，與上述Case1類型個體的生長模式差別更大，其生長率在初始階段保持最大值，此后漸次減小，沒有拐點(diǎn)存在。歸屬Case8類型的3個體，與前述兩種類型均不同，有正常拐點(diǎn)值，且a值也在正常范圍(28 573～279 667之間)。

考慮到樹干胸高斷面積用Chapman-Richards模型模擬時，多數(shù)個體a值和拐點(diǎn)值存在異?，F(xiàn)象，本文采用Logistic模型重新進(jìn)行模擬，結(jié)果見表2。從表2模型參數(shù)反映的生長過程最大值a值和拐點(diǎn)值看，除B1株a值和拐點(diǎn)值異常外，其余均在正常范圍。說明Logistic模型比Chapman-Richards模型更能不失真地模擬馬尾松胸高斷面積生長過程。

根據(jù)個體C3和B2的Logistic模型參數(shù)a值反推計(jì)算胸高直徑，預(yù)估的最大胸高直徑分別達(dá)5.52 m和4.80 m，其余預(yù)估胸徑最大值均在1.42～4.68 m之間。湖南現(xiàn)存馬尾松直徑最大個體位于炎陵縣鹿原鎮(zhèn)，樹齡680年，胸徑達(dá)1.7 m[25]。模型預(yù)估值雖然是現(xiàn)存最大古馬尾松樹胸徑的1～3倍，吻合度不高，但可作為檢驗(yàn)?zāi)Ｐ秃侠硇缘膮⒖?。解析的馬尾松個體中，共有5株拐點(diǎn)出現(xiàn)于50～60年之間。說明馬尾松胸高斷面積快速生長期持續(xù)時間較長，遠(yuǎn)高于以往認(rèn)為的20年左右[5]。除此之外，拐點(diǎn)在30～39年之間的有6株；拐點(diǎn)在20～29年的有5株；13～19年的有4株。拐點(diǎn)出現(xiàn)于50～60年、20～29年間的株數(shù)各占23.81%，拐點(diǎn)出現(xiàn)于30～39年的株數(shù)占28.57%，13～19年的株數(shù)占19.05%。唯有B1的拐點(diǎn)在123年出現(xiàn)，相比屬于生長模式異常株，占4.76%。拐點(diǎn)大于30年的個體數(shù)占57.14%。即林分一半以上馬尾松個體胸高斷面積最大生長速度出現(xiàn)在30年之后，28.57%以上個體出現(xiàn)于50年之后。

將靖州樣地13棵近同齡馬尾松解析木胸高斷面積數(shù)據(jù)進(jìn)行混合，用Logistic函數(shù)建模，得參數(shù)如下：a=97497.2880，b=17.4541，c=0.064376(決定系數(shù)(DC) 0.625096;F統(tǒng)計(jì)值1236.9732)。據(jù)此，計(jì)算得理論最大胸高斷面積31034.38cm2，最大胸高直徑3.52m，拐點(diǎn)出現(xiàn)于44.4年。依據(jù)彭其龍等[9]馬尾松個體胸高斷面積生長過程Logistic模型參數(shù)，計(jì)算得拐點(diǎn)值出現(xiàn)于55年。根據(jù)劉延男等[6]55年生馬尾松個體材積生長的Logistic模型參數(shù)計(jì)算，生長曲線拐點(diǎn)出現(xiàn)于48.1年。考慮到斷面積與材積生長的相關(guān)性，由上述不同研究結(jié)果獲得的3個拐點(diǎn)值均處于50年前后。因此，將馬尾松大徑材林分經(jīng)營年限設(shè)定在50年左右是合理的。

4 結(jié)論與討論

通過比較Chapman-Richards與Logistic函數(shù)擬合21株馬尾松大徑材胸高斷面積生長過程模型參數(shù)特點(diǎn)，結(jié)合林木生長規(guī)律常識，判斷Chapman-Richards模型擬合馬尾松個體胸高斷面積生長過程數(shù)據(jù)時，部分出現(xiàn)不合常理的異常參數(shù)，如a值和拐點(diǎn)值過大，而Logistic模型模擬參數(shù)基本正常。我們還采用其它模型進(jìn)行過模擬，諸如Gompertz、Weibull、李新運(yùn)自適應(yīng)Logistic模型、崔-Lawson模型、張大勇Logistic模型、宋丁全Logistic模型1、模型2[26-27]等，顯示各模型均有多個異常值，不適合用于對馬尾松個體胸高斷面積生長過程的擬合。本文結(jié)果表明，在擬合馬尾松個體胸高斷面積生長過程數(shù)據(jù)時，Logistic函數(shù)在眾多模型中，是相對較優(yōu)的。這一結(jié)果與陳哲夫等[11]用5種函數(shù)擬合馬尾松個體胸高斷面積時Logistic函數(shù)擬合度最高的結(jié)論相符。

值得注意的是，通過1stOpt軟件，嘗試用上述多種函數(shù)擬合斷面積數(shù)據(jù)，均能獲得良好的擬合效果，即R2值均能達(dá)到0.95以上，有的甚至高達(dá)0.999以上，體現(xiàn)出該軟件很強(qiáng)的參數(shù)擬合能力。但所獲參數(shù)是否具有生物學(xué)意義，則要根據(jù)林學(xué)專業(yè)經(jīng)驗(yàn)來判斷，不可簡單以相關(guān)系數(shù)和殘差平方和等值作為判斷模型優(yōu)劣的依據(jù)。這與FEKEDULEGN等[19]的觀點(diǎn)相似。

各馬尾松林分個體，雖起源相似，年齡相仿，但生長模式表現(xiàn)出較明顯差異。主要是生長過程曲線拐點(diǎn)出現(xiàn)時間不同，以致呈現(xiàn)出明顯的異速生長特點(diǎn)。推測這是馬尾松近同齡林分中一種較普遍的現(xiàn)象，與林分個體遺傳基礎(chǔ)及微域環(huán)境差異有關(guān)[28]。馬尾松胸高斷面積生長曲線的拐點(diǎn)值，反映的是生長模式由較高生長速度向較低生長速度轉(zhuǎn)型的時間點(diǎn)。根據(jù)拐點(diǎn)值推斷，可以將馬尾松大徑材經(jīng)營年限設(shè)定在50 年左右。