歐建德 歐家琳 康永武

（1. 明溪縣林業(yè)局，福建 三明 365200；2. 福建理工學校，福建 福州 350000；3. 沙縣林業(yè)局，福建 三明 365500）

森林生物量是森林生態(tài)系統(tǒng)中最基礎(chǔ)的數(shù)量特征[1]，是影響并決定著森林經(jīng)營決策的基礎(chǔ)。單木生物量模型是估算森林生物量的基礎(chǔ)[2]，提高單木生物量模型的預(yù)測精度一直是林業(yè)科技的重點與難點。模型選擇與優(yōu)化[1?2]、自變量的選取[3]是提高生物量模型預(yù)測精度的主要途徑。鑒于當前植物生物量模型的建立大多基于生長性狀[4?6]變量進行，未見導入總體樹冠結(jié)構(gòu)特征自變量。考慮到樹冠結(jié)構(gòu)改變植物光合作用、同化能力[7?8]，影響生物量積累，選取導入樹冠特征作為自變量，以提高生物量模型的預(yù)測精度十分重要。

云南紅豆杉（Taxus yunnanensis）（以下簡稱紅豆杉）的紫杉醇半合成前體10?脫乙酰巴卡亭Ⅲ（10?DAB）累積含量高[9]，是紫杉烷原料林的主栽樹種和珍貴用材樹種，建立并優(yōu)化單木生物量模型，精準預(yù)測原料林生物量（產(chǎn)量）十分必要。當前原料林研究多關(guān)注遺傳規(guī)律[10?11]、培育[12?13]、采收季節(jié)[14?15]等，尚未見單木生物量模型報道。冠幅、冠形率和樹冠率是樹冠形態(tài)的二維指標體系[7,16?18]，亦是調(diào)控樹冠形態(tài)結(jié)構(gòu)的指標性狀[7,16?18]。因此，以福建省明溪縣3年生云南紅豆杉幼林為研究對象，在生長性狀（地徑、樹高）自變量基礎(chǔ)上，全面導入冠幅、冠形率、樹冠率（樹冠形態(tài)結(jié)構(gòu)）等自變量，運用逐步回歸和獨立樣本檢驗的方法，旨在建立、評價并優(yōu)化紅豆杉系列單木生物量模型，明確理想的樹冠形態(tài)結(jié)構(gòu)與調(diào)控措施，以期為紫杉烷原料林經(jīng)營提供技術(shù)與理論支持。

1 研究區(qū)概況

研究地在福建省西北部的明溪縣（26 °8′～26 °39 ′ N、117 °4 ′～118 °47 ′ E），主要以低山丘陵地貌為主，地處武夷山脈南面，氣候?qū)僦衼啛釒ШＱ笮约撅L氣候，年均溫18.1 ℃，年均降水量1786 mm，5—6月最多，全年平均蒸發(fā)量1364 mm，年均日照時數(shù)1750.7 h，無霜期283 d，年均相對濕度81%[19]。供試材料取自明溪縣王橋的3年生云南紅豆杉原料林林分，株行距25 cm×25 cm，林分郁閉度0.7，平均地徑1.2 cm，平均樹高97 cm。海拔300 m，平均坡度5°左右，山地紅壤，立地肥沃。

2 研究方法

2.1 性狀測定

2019年12 月在原料林基地內(nèi)，隨機布設(shè)1 m×1 m的樣方8個，逐株標志并測量紅豆杉地徑、樹高、枝下高、東西和南北冠幅等，分別精確至0.1、1、1、1 cm，共計128株植株。采用全挖法，挖取樣方內(nèi)全部植株，帶回實驗室，測定枝葉生物量、莖干生物量、根系生物量。生物量測定參考文獻[20]的方法執(zhí)行。樣木基本情況見表1。

表1 樣本基本情況Table 1 Basic information of samples

樹冠形態(tài)特征及生物量按以下公式計算 ：

2.2 生物量數(shù)學模型的擬合與優(yōu)選

地徑、樹高、冠幅均勻分布，按照典型抽樣的方法抽取，選擇80%（103棵）紅豆杉樣本數(shù)據(jù)作為建模數(shù)據(jù)。生物量模型表達式如下：

式中：a、b、c、d、e、f是待解的參數(shù)；Y表示枝葉生物量（Ly）、莖干生物量（Sy），地上生物量（Ay），總生物量（Ty）；D表示地徑，H表示苗高，CW表示冠幅，CSR表示冠形率，CR表示樹冠率。

采用逐步回歸方法擬合數(shù)學模型，根據(jù)擬合模型判定系數(shù)（R2）的大小，篩選出最優(yōu)數(shù)學模型[1,23]。

2.3 優(yōu)選模型檢驗

采用獨立樣本檢驗評價模型預(yù)測水平和精度，即以剩余的20%（25棵）紅豆杉樣本數(shù)據(jù)作為檢驗數(shù)據(jù)，利用系列優(yōu)選生物量模型，分別計算出檢驗數(shù)據(jù)的生物量預(yù)測值后[1,23]，對預(yù)測值與實測值進行配對樣本t檢驗[1,23]、計算平均偏差（ME）、平均絕對偏差（MAE）、平均相對偏差（MPE）、平均相對偏差絕對值（MAPE）等系列偏差指標，估算預(yù)測模型精度（P）評價模型預(yù)測能力[1,24?26]；偏差檢驗與模型預(yù)測精度估算按參考文獻[24?26]的方法執(zhí)行。

2.4 數(shù)據(jù)處理

采用Excel 2003軟件進行數(shù)據(jù)預(yù)處理與統(tǒng)計；采用SPSS 21.0軟件中逐步回歸方法擬合數(shù)學模型，配對樣本t檢驗方法檢驗預(yù)測值、實測值間差異。

3 結(jié)果與分析

3.1 單木枝葉生物量模型擬合與比較

紅豆杉單木枝葉生物量模型的逐步回歸擬合結(jié)果見表2。在模型的擬合過程中，變量D、H、CW、CSR依次進入，始終未出現(xiàn)CR變量，說明CR單獨作用不顯著。在擬合模型1～3過程，模型參數(shù)均通過參數(shù)t檢驗（P<0.01），說明D、H、CW對枝葉生物量均有顯著直接作用。模型調(diào)整判定系數(shù)（R2）隨變量D、H、CW的依次進入而持續(xù)增大，導入CW變量可提高擬合效果。在模型3基礎(chǔ)上導入CSR變量，導致模型4中H變量未通過參數(shù)t檢驗（P=0.158>0.05），為此繼續(xù)逐步回歸擬合得模型5，模型5的系列參數(shù)（常量、D、CW、CSR）均通過參數(shù)t檢驗（P<0.01）。由表2可知，模型5中所有參數(shù)均通過t檢驗且R2最大（0.812），為最優(yōu)模型；模型中CW、CSR變量系數(shù)分別為1.181、18.372，均大于0，表明冠幅和冠形率有顯著促進枝葉生物量效應(yīng)。結(jié)果（表2）還顯示，模型2為導入樹冠形態(tài)特征因子前（以下簡稱導入前）最優(yōu)模型（R2最大為0.762）。

表2 單木枝葉生物量數(shù)學模型擬合結(jié)果Table 2 Fitting results of the single tree branches and leaves biomass mathematical model of T. yunnanensis

為分析導入樹冠形態(tài)特征因子前后（以下簡稱導入前后）系列生物量模型的預(yù)測效果，采取獨立樣本檢驗方法，對系列導入樹冠形態(tài)特征前后的2種最優(yōu)模型（R2最大且模型參數(shù)均通過t檢驗）進行配對t檢驗、偏差檢驗和預(yù)測精度估算比較，結(jié)果見表3。

表3 優(yōu)化的單木紅豆杉生物量模型檢驗結(jié)果Table 3 Test results of optimized single tree biomass model of T. yunnanensis

由表2可知，導入前后最優(yōu)模型（模型2、5）的R2分別0.762、0.812，說明導入后提高了模型的擬合效果。由表3可知，導入前后模型配對t檢驗的P值分別為0.040、0.620，意味著導入前模型的實測值、擬合值間差異顯著，而導入后模型實測值、擬合值間差異不顯著，表明導入后提高了模型的預(yù)測效果。導入前后模型的預(yù)測精度分別達89.67%、92.50%，導入后提高了模型預(yù)測精度。導入后模型的ME、MAE、MPE、MAPE分別為0.81 g、5.94 g、?5.22%、18.64%，總體較導入前模型理想。

綜上，導入樹冠形態(tài)特征因子后明顯提高模型的擬合與預(yù)測效果，與前文中冠幅和冠形率顯著促進枝葉生物量的結(jié)論相吻合。

3.2 單木莖干生物量模型擬合與比較

紅豆杉單木莖干生物量模型的逐步回歸擬合結(jié)果見表4。在模型的擬合過程中，變量D、H、CW、CSR依次進入，始終未出現(xiàn)CR變量，說明CR的單獨作用效應(yīng)不顯著；在模型6～8的擬合過程，模型的參數(shù)均通過參數(shù)t檢驗（P<0.01），說明D、H、CW對莖干生物量均有顯著直接作用；系列模型的R2隨變量D、H、CW的依次進入而持續(xù)增大，說明導入CW變量提高了模型擬合效果。在模型8基礎(chǔ)上導入CSR得模型9，但模型9中H變量未通過參數(shù)t檢驗（P=0.198>0.05），為此繼續(xù)逐步回歸擬合得模型10。模型10參數(shù)（常量、D、CW、CSR）均通過參數(shù)t檢驗（P<0.01）。綜合6～10系列模型的參數(shù)t檢驗和R2值表現(xiàn)（表4），模型10所有參數(shù)均通過參數(shù)t檢驗（P<0.01）且R2值最大（0.818），為最優(yōu)模型。模型中CW、CSR系數(shù)分別為0.540、8.105，均大于0，表明冠幅和冠形率顯著促進莖干生物量效應(yīng)。結(jié)果（表4）還顯示，模型7為導入樹冠形態(tài)特征前最優(yōu)模型（R2最大為0.773）。

結(jié)果（表4）顯示，導入樹冠形態(tài)特征前后最優(yōu)模型（模型7、10）的R2分別為0.773、0.818，說明導入后提高了模型的R2、擬合效果。結(jié)果（表3）顯示，導入前后最優(yōu)模型配對t檢驗的P值分別為0.033、0.527，意味著導入前的擬合值、實測值間差異顯著，而導入后的實測值、擬合值差異不顯著、預(yù)測效果好，表明導入后提高模型的預(yù)測效果。導入前后最優(yōu)模型的預(yù)測精度分別為88.61%、92.36%，即導入后提高了模型的預(yù)測精度。導入后最優(yōu)模型（模型10）的ME、MAE、MPE、MAPE分 別 為0.46 g、2.29 g、?2.73%、18.95%，較導入前（模型7）的相應(yīng)偏差性狀均有更好表現(xiàn)。

表4 單木莖干生物量數(shù)學模型擬合結(jié)果Table 4 Fitting results of the single tree stem biomass mathematical model of T. yunnanensis

綜上，導入樹冠形態(tài)特征因子明顯提高模型的擬合與預(yù)測效果，與前文中冠幅和冠形率顯著促進莖干生物量的結(jié)論相吻合。

3.3 單木地上生物量模型擬合與比較

紅豆杉單木地上生物量模型的逐步回歸擬合結(jié)果見表5。在模型的擬合過程中，變量D、H、CW、CSR依次進入，始終未出現(xiàn)CR變量，說明CR單獨作用不顯著；在模型11～13的擬合過程，模型參數(shù)均通過參數(shù)t檢驗（P<0.01），說明D、H、CW變量對地上生物量均有著顯著的直接作用；模型R2隨著變量D、H、CW的依次進入而持續(xù)增大、擬合效果持續(xù)變好，說明導入CW變量可提高模型擬合效果。在模型13基礎(chǔ)上導入CSR變量得模型14，但模型14中H變量未通過參數(shù)t檢驗（P=0.170>0.05），繼續(xù)回歸擬合得模型15。綜合系列模型11～15的參數(shù)t檢驗和R2值表現(xiàn)（表5），模型15的參數(shù)均通過的參數(shù)t檢驗（P<0.01）、且R2值最大（0.815），為最優(yōu)模型。模型中CW、CSR的系數(shù)分別為1.721、26.427，均大于0，表明冠幅和冠形率對地上生物量有顯著促進效應(yīng)。結(jié)果（表5）還顯示，模型12為導入樹冠形態(tài)特征前最優(yōu)模型（R2最大、0.766）。

由表5可知，導入前后最優(yōu)模型（模型12、15）的R2分別為0.766、0.815，說明導入后提高了模型的R2、擬合效果。由表3可知，導入前后配對t檢驗的P值分別為0.041、0.587，意味著導入前擬合值、實測值間差異顯著，預(yù)測效果差，而導入后擬合值、實測值間差異不顯著，預(yù)測效果好，表明導入樹冠形態(tài)特征可提高模型的預(yù)測效果。導入前后最優(yōu)模型的預(yù)測精度分別89.39%、92.51%，即導入后提高了模型的預(yù)測精度。導入后最優(yōu)模型（模型15）的ME、MAE、MPE、MAPE分別為1.28 g、8.22 g、?4.47%、19.82%，總體較導入前（模型12）的表現(xiàn)理想。

表5 單木地上生物量數(shù)學模型擬合結(jié)果Table 5 Fitting results of the single tree aboveground biomass mathematical model of T. yunnanensis

綜上，導入樹冠形態(tài)特征因子明顯提高模型的擬合與預(yù)測效果，與前文中冠幅和冠形率顯著促進地上生物量的結(jié)論相吻合。

3.4 單木紅豆杉總生物量模型擬合與比較

紅豆杉單木總生物量模型的逐步回歸擬合結(jié)果見表6。在模型的擬合過程中，變量D、H、CR、CW依次進入，始終未出現(xiàn)CSR變量，說明冠形率CSR單獨作用不顯著。模型R2隨變量D、H、CR、CW的依次進入而持續(xù)增大、擬合效果持續(xù)變好，說明導入樹冠形態(tài)特征因子變量提高了模型擬合效果。綜合系列模型16～19的參數(shù)t檢驗、R2值表現(xiàn)（表6），模型19中參數(shù)均通過參數(shù)t檢驗（P<0.01）、且R2值最大（0.831），為最優(yōu)模型。模型中CR、CW的系數(shù)分別為65.347、0.808，均大于0，表明樹冠率和冠幅對紅豆杉總生物量有著顯著性促進效應(yīng)。結(jié)果（表6）還顯示，模型17為導入前最優(yōu)模型（R2最大為0.780）。

由表6可知，導入樹冠形態(tài)特征前后最優(yōu)模型（模型17、19）的R2分別為0.780、0.831，說明導入后提高了模型的R2、擬合效果。由表3可知，導入前后配對t檢驗的P值分別為0.060、0.587，雖然導入前后擬合值、實測值均無顯著性差異，但導入后的配對t檢驗的P=0.587大于導入前的P=0.060，表明導入樹冠形態(tài)特征可提高模型的預(yù)測效果。導入前后最優(yōu)模型的預(yù)測精度分別90.57%、92.58%，即導入后提高了模型的預(yù)測精度。導入后最優(yōu)模型（模型19）的ME、MAE、MPE、MAPE分別為0.31 g、9.64 g、?4.26%、13.96%，總體較導入前（模型17）的表現(xiàn)理想。

表6 單木紅豆杉總生物量數(shù)學模型擬合結(jié)果Table 6 Fitting results of the single tree total biomass mathematical model of T. yunnanensis

綜上，導入樹冠形態(tài)特征因子后明顯提高模型的擬合與預(yù)測效果，與前文中樹冠率和冠幅顯著促進紅豆杉總生物量的結(jié)論相吻合。

4 結(jié)論與討論

本研究全面導入系列樹冠形態(tài)特征因子（冠幅、冠形率、樹冠率）自變量后，建立了紅豆杉系列單木生物量模型，提高了模型的預(yù)測精度與水平，可為同類研究提供參考與借鑒；系列單木生物量優(yōu)化模型預(yù)測精度均在92.36%以上，且偏差表現(xiàn)良好，預(yù)測精度高，可用于紅豆杉原料林生物量與產(chǎn)量預(yù)測預(yù)估。建議在云南紅豆杉原料林培育過程中，首先采用良種壯苗、優(yōu)化密度和幼林撫育以促進地徑生長，其次宜結(jié)合培育模式（收獲物的差異）調(diào)適樹冠結(jié)構(gòu)調(diào)控措施與策略。

研究結(jié)果顯示，單木枝葉、莖干和地上生物量模型的實測值與擬合值差異顯著（P<0.05），且預(yù)測精度均小于90%，表明擬合紅豆杉單木枝葉、莖干和地上生物量模型僅導入地徑、苗高等生長自變量是不夠的，該系列單木生物量模型中導入樹冠形態(tài)特征變量是必要的，驗證了生物量模型主要以與冠幅相關(guān)的因子為自變量[27]以及含有樹冠形態(tài)特征變量的單木生物量預(yù)測精度表現(xiàn)好[28]的研究結(jié)論一致。

系列單木生物量模型的逐步回歸擬合過程顯示，率先進入的變量均為地徑，其次是樹高，且均表現(xiàn)出顯著促進作用，表明地徑在直接促進系列生物量的作用超過樹高，是主導影響因子。提高原料林生物量（產(chǎn)量），首先是促進地徑生長。進入系列最優(yōu)生物量模型的樹冠形態(tài)特征與先后順序明顯變化，表明樹冠形態(tài)特征作用系列生物量規(guī)律，以及利于系列單木生物量積累的理想樹冠形態(tài)結(jié)構(gòu)與調(diào)控措施明顯不同，建議生產(chǎn)中宜根據(jù)目標收獲物的差異調(diào)適樹冠結(jié)構(gòu)調(diào)控措施與策略。這與因調(diào)控目標改變調(diào)整樹冠調(diào)控措施[7,17?18]的研究結(jié)論相互驗證。先后依次進入單木枝葉、莖干和地上生物量最優(yōu)模型的樹冠形態(tài)特征是冠幅、冠形率，且均有顯著促進生物量積累作用。這意味著促進該系列生物量積累的理想樹冠形態(tài)結(jié)構(gòu)首先是冠幅寬大、其次為狹長的冠形。促進生物量積累的樹冠形態(tài)結(jié)構(gòu)調(diào)控重點與方向是促進冠幅寬大，建議在以采收枝葉、地上生物量為目標的原料林矮林作業(yè)模式中，宜選擇保留或采用冠幅寬大的個體、品種為主，或減少密度以促進冠幅生長。先后依次進入單木總生物量最優(yōu)模型的樹冠形態(tài)特征為樹冠率、冠形率，且均有顯著促進總生物量積累作用，意味著促進總生物量理想樹冠形態(tài)結(jié)構(gòu)首先是大的樹冠率，其次是冠幅寬大。促進總生物量積累的樹冠形態(tài)調(diào)控重點與方向是提高樹冠率，建議在全挖式采收的短周期原料林模式中，宜保留或采用樹冠率大的個體、品種為主，以及經(jīng)營過程減少降低樹冠率的行為措施。

導入樹冠形態(tài)結(jié)構(gòu)特征，采用逐步回歸的方法，明顯提升了云南紅豆杉單木生物量模型的預(yù)測精度與水平，與前人采用逐步回歸的方法優(yōu)化預(yù)測模型[4?6,20,29]以及在自變量不斷增加情況，模型的擬合精度有不斷提升的趨勢[28]結(jié)論相互驗證。結(jié)果表明樹冠結(jié)構(gòu)顯著影響原料林單木生物量積累，這是原料林樹冠既是枝葉生產(chǎn)器官，又是光合作用的同化器官，其樹冠結(jié)構(gòu)改變了利用地上空間能力、影響同化能力與生產(chǎn)力所致。這與前人關(guān)于樹冠結(jié)構(gòu)影響同化能力、地上空間利用能力[7?8]和林木生產(chǎn)力[7?8,30]的研究結(jié)論相互驗證。王軼夫等[31]認為冠長、冠幅等樹冠自變量與馬尾松（Pinus massoniana）樹冠、樹枝、葉花果等生物量相關(guān)關(guān)系不明顯；造成結(jié)論差異的原因，與樹種、培育模式以及年齡不同有關(guān)。同時亦驗證了生物量模型的自變量選取因樹種明顯不同[27]的結(jié)論。

本研究僅從導入樹冠形態(tài)特征自變量方面，擬合并優(yōu)化單木生物量模型，有關(guān)模型的自變量類型、模型選擇等優(yōu)化研究尚待后續(xù)。