賈寶軍，胡靖揚(yáng)，林 寬，馮倩男，劉常富,2，于立忠

（1.沈陽(yáng)農(nóng)業(yè)大學(xué) 林學(xué)院，遼寧 沈陽(yáng) 110866；2.中國(guó)林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護(hù)研究所，北京100091；3.中國(guó)科學(xué)院清原森林生態(tài)系統(tǒng)觀測(cè)研究站，遼寧 沈陽(yáng) 110016）

遼東山區(qū)落葉松林枝條葉面積模型選擇

賈寶軍1，胡靖揚(yáng)1，林 寬1，馮倩男1，劉常富1,2，于立忠3

（1.沈陽(yáng)農(nóng)業(yè)大學(xué) 林學(xué)院，遼寧 沈陽(yáng) 110866；2.中國(guó)林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護(hù)研究所，北京100091；3.中國(guó)科學(xué)院清原森林生態(tài)系統(tǒng)觀測(cè)研究站，遼寧 沈陽(yáng) 110016）

本研究基于遼東山區(qū)實(shí)測(cè)落葉松特征因子，通過(guò)枝條解析法獲得了254組枝條的葉面積、枝條基徑、枝條長(zhǎng)度等屬性數(shù)據(jù)，在枝條水平上采用一元非線性方程和二元及多元非線性方程建立以枝條屬性因子為自變量的枝條葉面積模型。結(jié)果表明：枝條基部斷面積(d2)與枝條葉面積有最高的相關(guān)性，最優(yōu)一元枝條葉面積模型為：y=8.967/(1+50.901e-0.084d)，其R2達(dá)到0.719，測(cè)算精度為86.34%（α=0.05）。引入相對(duì)著枝深度(RDINC)，使得二元及多元模型的決定系數(shù)顯著提高，最優(yōu)二元枝條葉面積模型為：y=0.002(d2)2.260e-1.701RDINC，其R2達(dá)到0.796，測(cè)算精度為88.57%（α=0.05）。

落葉松；葉面積；枝水平；回歸模型

葉是植物的營(yíng)養(yǎng)器官、光合作用的主要場(chǎng)所，是森林生態(tài)系統(tǒng)中能量流動(dòng)、物質(zhì)循環(huán)、特別是水文循環(huán)的重要生物調(diào)節(jié)器，是研究森林生態(tài)系統(tǒng)生產(chǎn)力、蒸騰量、光合效率及林冠截流不可缺少的一個(gè)環(huán)節(jié)，所以葉面積是研究森林生態(tài)系統(tǒng)結(jié)構(gòu)和功能的一個(gè)重要指標(biāo)[1]。由于光合作用的主要發(fā)生器官是葉片，因此，葉面積的大小及其分布直接影響著林分對(duì)光能的截獲及利用，進(jìn)而影響著林分生產(chǎn)力[2]。同時(shí)，葉面積的大小對(duì)樹(shù)木的生長(zhǎng)發(fā)育等眾多生命活動(dòng)影響深遠(yuǎn)，尤其和樹(shù)木生物量的積累密切相關(guān)[3]。所以葉面積估算是研究氣體交換和生態(tài)系統(tǒng)模擬的重要基礎(chǔ)[4]，探究將葉片水平上測(cè)得的數(shù)據(jù)轉(zhuǎn)換成為枝條、單株或林分水平上的結(jié)果，是目前樹(shù)木生理生態(tài)學(xué)研究領(lǐng)域中的一個(gè)重要理論問(wèn)題[5]。

樹(shù)冠是樹(shù)木進(jìn)行光合作用的場(chǎng)所，冠層內(nèi)葉面積的分布與光能的截獲緊密相關(guān)，影響樹(shù)木生物量的積累和林分的生長(zhǎng)[6]。因此，準(zhǔn)確的估算葉面積有助于更深入地了解樹(shù)冠結(jié)構(gòu)，以及育林生產(chǎn)[7]。直接測(cè)算葉面積雖能客觀準(zhǔn)確的測(cè)定樹(shù)冠葉面積，但是這種方法既費(fèi)時(shí)又費(fèi)力，并且會(huì)對(duì)生態(tài)系統(tǒng)造成嚴(yán)重的破壞。植物分枝作為一個(gè)獨(dú)立的分形結(jié)構(gòu)，與植物地上部分是一個(gè)統(tǒng)計(jì)自相似系統(tǒng)，因此，可以通過(guò)簡(jiǎn)單易測(cè)的樹(shù)木特征值來(lái)模擬較難測(cè)算的更具生態(tài)學(xué)意義的指標(biāo)[8]。近些年，許多學(xué)者對(duì)葉面積模擬進(jìn)行了大量研究。有些研究基于葉長(zhǎng)和葉寬等葉片屬性因子對(duì)單片葉面積進(jìn)行模擬，如劉建偉等利用主脈長(zhǎng)及葉寬模擬楊樹(shù)單葉面積取得良好效果[9]；吳春胤等利用BP神經(jīng)網(wǎng)絡(luò)模型簡(jiǎn)單、快速、經(jīng)濟(jì)地測(cè)定荔枝樹(shù)葉面積[10]。還有一些學(xué)者[11-14]通過(guò)建立基于胸徑和樹(shù)高等因子的葉面積模型對(duì)樹(shù)冠葉面積進(jìn)行估算，然而對(duì)于葉面積分布、葉面積指數(shù)等研究其精度難以達(dá)到要求。因此，近些年建立基于枝條屬性因子的葉面積估算模型已成為一種趨勢(shì)。戚繼忠對(duì)紫椴和紫丁香研究發(fā)現(xiàn)，每個(gè)個(gè)體枝橫斷面參數(shù)中以橫斷面積與葉面積呈現(xiàn)最顯著相關(guān)[15:94]。郭孝玉對(duì)長(zhǎng)白落葉松枝條葉面積進(jìn)行研究，結(jié)果表明：以枝條直徑為變量的葉面積異速生長(zhǎng)模型最優(yōu)[16]。

本文以遼東山區(qū)落葉松林為研究對(duì)象，采用全株收獲法，獲取枝條解析因子及枝條葉面積，選擇多種非線性方程及不同變量組合建立枝條葉面積估算模型，并對(duì)最優(yōu)葉面積模型的預(yù)測(cè)效果進(jìn)行分析，確定最優(yōu)估算模型，以期為遼東山區(qū)落葉松樹(shù)冠葉面積和林分葉面積準(zhǔn)確估算奠定基礎(chǔ)。

1 研究方法

1.1 研究區(qū)概況

研究區(qū)位于遼寧省清原縣中國(guó)科學(xué)院清原森林生態(tài)系統(tǒng)觀測(cè)研究站（41°51.102′N(xiāo)，124°54.543′E），地處清原滿(mǎn)族自治縣南部，海拔252～1 116 m。土壤主要類(lèi)型為暗棕和棕色森林土，多呈酸性和中性。氣候類(lèi)型為溫帶大陸濕潤(rùn)季風(fēng)氣候，年平均相對(duì)濕度65%～75%，年平均氣溫3.9～5.4℃，年平均降水量700～850 mm。區(qū)內(nèi)植被屬于長(zhǎng)白植物區(qū)系，是東北地區(qū)次生林的主要代表林型?，F(xiàn)存主要喬木樹(shù)種有紅松Pinus koraiensis、日本落葉松Larix kaempferi、黃花落葉松L. olgensis等。

1.2 數(shù)據(jù)收集與整理

1.2.1 樣地設(shè)置及標(biāo)準(zhǔn)木選取

于2014年9月份，在清原站試驗(yàn)林場(chǎng)內(nèi)選取3塊不同密度、不同年齡的落葉松林，每塊標(biāo)準(zhǔn)樣地面積為20 m×20 m，所選林分均生長(zhǎng)良好。對(duì)每塊標(biāo)準(zhǔn)地進(jìn)行每木調(diào)查，測(cè)量樹(shù)高（H）、胸徑（D）、冠幅（W）等測(cè)樹(shù)因子，計(jì)算林分平均胸徑及平均樹(shù)高，并在每塊標(biāo)準(zhǔn)地附近選取1株冠形良好、無(wú)病蟲(chóng)害的解析木。樣地概況見(jiàn)表1。

表1 落葉松林標(biāo)準(zhǔn)地概況Table 1 Basic information of different larch forest

1.2.2 枝條解析

采取全株收獲法，小心伐倒解析木以免損壞枝條，首先測(cè)定樹(shù)高、胸徑、枝下高等因子，在樹(shù)干上標(biāo)明北向，沿樹(shù)干每米做標(biāo)記。然后從樹(shù)冠基部第一個(gè)活枝起，依次向上對(duì)枝條進(jìn)行編號(hào)。用米尺測(cè)量枝條距樹(shù)冠頂部的距離，得到著枝深度（DINC），用DINC/CL得到相對(duì)著枝深度（RDINC），再測(cè)量枝條的長(zhǎng)度，用游標(biāo)卡尺測(cè)量枝條直徑。解析木及枝條概況見(jiàn)表2。

表2 落葉松林解析木和枝條解析因子概況Table 2 Summary of sample tree and branch attributes for larch forest

1.2.3 葉面積測(cè)定

采用分段收獲法將一級(jí)枝上的葉片全部收獲，由枝條基部開(kāi)始按50 cm間隔將一級(jí)枝分段（枝梢不足50 cm標(biāo)記實(shí)際長(zhǎng)度），然后將各段葉片全部收獲，分別用自封袋封裝并編號(hào)。當(dāng)天立即隨機(jī)選取一定數(shù)量的葉片，用掃描儀進(jìn)行掃描，然后將掃描后的圖片用Photoshop CS5軟件自帶的分析功能測(cè)算葉面積。然后將樣品烘干，用樣葉面積除以樣葉干重得到每段樣品的比葉面積。最后將葉片全部置于105℃下烘干至恒重，用每段葉片的干重乘以每段的比葉面積求出每個(gè)枝條的葉面積。

1.3 研究方法

1.3.1 一元枝條葉面積模型

落葉松枝條葉面積與屬性因子的相關(guān)系數(shù)分別為基徑d（0.785），枝長(zhǎng)l（0.713），基部斷面積d2（0.827），變量組合d2l（0.794），著枝深度DINC（0.497），相對(duì)著枝深度RDINC（0.419），由此可見(jiàn)枝條葉面積與著枝深度DINC和相對(duì)著枝深度RDINC相關(guān)性較差，所以選擇變量d、l、d2、d2l作為預(yù)選變量。同時(shí)，根據(jù)散點(diǎn)圖分析可知，4個(gè)變量與落葉松枝條葉面積存在明顯的非線性關(guān)系，采用Schumacher方程、Richards方程、Roeecp方程、修正Weibull方程、Levacovic 方程、Logistic方程、Mitscherlich方程、冪函數(shù)方程、Korf方程、Compertz方程共10個(gè)模型進(jìn)行建模，具體公式見(jiàn)文獻(xiàn)[17]。

1.3.2 二元及多元枝條葉面積模型

初選 5個(gè)自變量d、d2、l、DINC、RDINC，模型選擇多元冪函數(shù)多元線性Roceep方程和以最優(yōu)一元模型為基礎(chǔ)的復(fù)合模型為最優(yōu)一元模型）構(gòu)建二元及多元枝條葉面積模型。

1.3.3 模型評(píng)價(jià)與檢驗(yàn)

本文選取決定系數(shù)（R2）、均方誤差（RMSE）、剩余標(biāo)準(zhǔn)差（SEE）評(píng)價(jià)模型的擬合效果。模型采用總相對(duì)誤差（RS）、平均相對(duì)誤差(EE)、平均相對(duì)誤差絕對(duì)值(RMA)和預(yù)估精度(P)4個(gè)指標(biāo)評(píng)價(jià)模型預(yù)測(cè)能力。同時(shí)，采用參數(shù)變動(dòng)系數(shù)（CV）來(lái)評(píng)價(jià)模型的穩(wěn)定性。各指標(biāo)計(jì)算公式見(jiàn)參考文獻(xiàn)[3,16]。模型檢驗(yàn)參見(jiàn)文獻(xiàn)[18]。隨機(jī)選取80%數(shù)據(jù)建模，20%數(shù)據(jù)對(duì)模型進(jìn)行檢驗(yàn)，模型計(jì)算采用SPSS19.0和Excel統(tǒng)計(jì)分析軟完成。

2 結(jié)果與分析

2.1 一元枝條葉面積模型

以落葉松枝條葉面積為因變量，選用10類(lèi)模型、4個(gè)自變量分別進(jìn)行回歸擬合共得到40個(gè)模型，以決定系數(shù)為評(píng)價(jià)指標(biāo)，選出R2最高的10個(gè)枝條葉面積模型擬合結(jié)果見(jiàn)表3。從擬合結(jié)果看出，模型6的R2最高達(dá)到0.719，其次依次是模型10（0.715）、模型3（0.713）、模型4、5、8（0.708）、模型7（0.707）、模型2（0.706），其余均小于0.7。然而在保證高決定系數(shù)的同時(shí)，還要考慮參數(shù)的穩(wěn)定性，模型2、4、5、7、10中的CV%大于50%相對(duì)不穩(wěn)定，因此，初步選定1、3、6、8、9作為備選模型，進(jìn)一步進(jìn)行誤差統(tǒng)計(jì)量分析。

表3 一元枝條葉面積估計(jì)模型的擬合結(jié)果Table 3 Fitting result of branch leaf area models for unitary non-linearity regress equation

表4列出了模型擬合時(shí)5種備選模型的誤差統(tǒng)計(jì)量，模型6的RS絕對(duì)值（0.042%）最小，依次為模型3、模型8、模型9、模型1。模型3的EE（4.436%）最小，依次為模型6、模型9、模型8、模型1。模型6的RMA（45.213%）最小，其次為模型3、模型9、模型8、模型1。模型6的P（90.793%）最高。

從擬合結(jié)果和模型的誤差統(tǒng)計(jì)量綜合來(lái)看，模型6可以作為一元落葉松枝條葉面積模型。然而決定系數(shù)較高尚不足以判定所建立的模型是否有偏，僅對(duì)綜合統(tǒng)計(jì)量進(jìn)行分析，并不能發(fā)現(xiàn)擬合結(jié)果數(shù)據(jù)樣式中的差異，需要進(jìn)行殘差分析。模型6（見(jiàn)圖1）殘差呈均勻分布，沒(méi)有明顯的圖樣，其中大多數(shù)分布在±1.5之間，表明擬合效果是良好的，極個(gè)別樣本點(diǎn)點(diǎn)落在±1.5之外，是因?yàn)檫@些樣本的枝條處在樹(shù)冠底層，由于遮蔭和林木之間競(jìng)爭(zhēng)引起，導(dǎo)致枝條的葉片分布不均，葉量較少。綜上所述，模型6即y=8.967/(1+50.901e-0.084d)為最優(yōu)一元枝條葉面積模型，模型精度檢驗(yàn)為86.34%（a=0.05）（見(jiàn)表5）。

表4 備選一元枝條葉面積估計(jì)模型的誤差統(tǒng)計(jì)量Table 4 Statistics of models for selected unitary regress equation

圖1 最優(yōu)一元枝條葉面積模型殘差Fig.1 The residual distribution of models for the optimal unitary non-linearity regress equation

表5 最優(yōu)模型的模擬精度Table 5 Accuracy test of optimal regress equation

2.2 二元及多元枝條葉面積模型

本研究選取d、l、DINC、RDINC、d2為自變量，分別用多元線性、多元冪函數(shù)、二元Roeecp方程和以Logistic為基礎(chǔ)的復(fù)合模型進(jìn)行回歸，共得到46個(gè)模型（12個(gè)線性方程、12個(gè)冪函數(shù)方程、16個(gè)Roeecp方程及6個(gè)復(fù)合模型）。經(jīng)過(guò)初步擬合，多元線性方程的誤差統(tǒng)計(jì)量均較大，將其舍去。以決定系數(shù)為標(biāo)準(zhǔn)，選取數(shù)值較高的3類(lèi)函數(shù)，共9個(gè)模型見(jiàn)表6。

在3類(lèi)函數(shù)中，以d、RDINC為變量的二元模型的擬合優(yōu)度均優(yōu)于以d、DINC和以d、l為變量的二元模型。模型5的R2最高達(dá)到0.796，依次為模型1和2 （0.787）、模型3（0.770）、模型8（0.763）、模型6（0.760）、模型7（0.750），其余均小于0.750。按照枝條葉面模型的參數(shù)變動(dòng)系數(shù)小于50%選擇模型5、7、8作為備選模型。

表7列出了模型擬合時(shí)3種備選模型的誤差統(tǒng)計(jì)量，模型5的RS、EE和RMA均為三種模型中最小，分別為0.171%、4.071%、42.658% 。預(yù)估精度最高，達(dá)到92.166%。模型5殘差（見(jiàn)圖2）呈均勻分布，沒(méi)有明顯的圖樣，多分布于±1.0之間。因此，模型5即y=0.002(d2)2.260e-1.701RDINC為二元枝條葉面積最優(yōu)模型，模型精度檢驗(yàn)為88.57%（a=0.05）（見(jiàn)表7）。

表6 多元枝條葉面積估計(jì)模型的擬合結(jié)果Table 6 Fitting result of branch leaf area models for multivariate non-linearity regress equation

圖2 最優(yōu)多元枝條葉面積模型殘差Fig.2 The residual distribution of models for the optimal multivariate non-linearity regress equation

表7 備選多元枝條葉面積估計(jì)模型的誤差統(tǒng)計(jì)量Table 7 Statistics of models for selected multivariate nonlinearity regress equation

3 結(jié)論與討論

研究表明枝條基部斷面積與枝條葉面積具有最高的相關(guān)性，最優(yōu)一元枝條葉面積模型為：y=8.967/(1+50.901e-0.084d)，其R2達(dá)到 0.719，測(cè)算精度為86.34%（α=0.05）。引入相對(duì)著枝深度，使得二元及多元模型的決定系數(shù)顯著提高，最優(yōu)二元枝條葉面積模型為：y=0.002(d2)2.260e-1.701RDINC，其R2達(dá)到0.796，測(cè)算精度為88.57%（α=0.05）。

以往學(xué)者建立的枝條葉面積模型決定系數(shù)基本在0.7以上[15,19]，本研究中最優(yōu)二元枝條葉面積的決定系數(shù)為0.8，同時(shí)均方差、殘差分布等均表現(xiàn)較好。本研究在以基徑為自變量的基礎(chǔ)上，把相對(duì)著枝深度引入到方程中，模型的決定系數(shù)從0.7提高到0.8，明顯提高了模型的擬合效果。這與Bartelink HH[20]研究結(jié)果一致，其原因是枝條葉面積主要受基徑和著生位置的影響。本研究發(fā)現(xiàn)枝條基部斷面積與枝條葉面積有最高的相關(guān)性，說(shuō)明枝條基部斷面積對(duì)枝條葉面積有顯著的影響，這與劉志國(guó)[21]等研究亞熱帶闊葉林植物時(shí)發(fā)現(xiàn)，小枝的莖截面積與總?cè)~面積之間呈異速生長(zhǎng)關(guān)系以及劉常富[8]和王東良[22]發(fā)現(xiàn)影響樹(shù)木三維綠量的最主要因素是胸高斷面積等研究相吻合。傳統(tǒng)林業(yè)研究焦點(diǎn)主要集中在樹(shù)木基徑（或胸徑）。而眾多研究表明，樹(shù)木基部斷面積（或胸高斷面積）在研究枝條生長(zhǎng)和林木生長(zhǎng)中的重要作用益發(fā)凸顯。

在模型構(gòu)建中，模型中自變量的個(gè)數(shù)并非越多越好，相反地，在估測(cè)精度允許的前提下，應(yīng)使模型所含的自變量個(gè)數(shù)盡可能的少。一方面可以減少因自變量帶來(lái)的工作量，另一方面，自變量多會(huì)使模型參數(shù)變動(dòng)較大，最終導(dǎo)致模型預(yù)測(cè)效果不理想[23]。由于葉片性狀對(duì)海拔高度變化比較敏感[24]，同時(shí)，可能針葉受光照、水分、營(yíng)養(yǎng)及土壤條件影響較大有關(guān)[25]，造成枝條葉面積大小變動(dòng)范圍很大。本研究發(fā)現(xiàn)把枝長(zhǎng)引入到葉面積模型中，模型的穩(wěn)定性明顯變差，其原因可能是基徑和枝長(zhǎng)之間存在明顯的共線性，導(dǎo)致模型不穩(wěn)定。然而考慮著生位置對(duì)枝條葉片的影響，選取枝條基徑和相對(duì)著枝深度構(gòu)建模型，既提高了模型的決定系數(shù)同時(shí)也保證了模型的穩(wěn)定性。快速準(zhǔn)確地估算葉面積，有助于林分、景觀以及地區(qū)尺度上對(duì)碳、能量、水分通量進(jìn)行研究[26]。李懷等通過(guò)用光合速率、蒸騰速率乘以葉面積、年日照時(shí)數(shù)，計(jì)算出綠地年吸二氧化碳量、放氧量，進(jìn)而計(jì)算出樹(shù)木生態(tài)效益價(jià)值[27]；石鐵矛等利用單株葉面積乘以單位葉面積釋氧能力值算出單株植物釋氧量，進(jìn)而研究植物的生態(tài)效益[28]；于佳等通過(guò)光合速率和葉面積估算沈陽(yáng)園林植物的碳匯值，進(jìn)而研究植物的碳匯能力[29]。因此，本研究建立的最優(yōu)枝條葉面積模型，可以為本地區(qū)的林木固碳、釋氧、滯塵等生態(tài)效益研究奠定基礎(chǔ)。

枝條基礎(chǔ)數(shù)據(jù)的獲取相對(duì)較困難，然而激光雷達(dá)技術(shù)的出現(xiàn)為無(wú)損精準(zhǔn)地獲取枝條基礎(chǔ)數(shù)據(jù)提供了可能。隨著計(jì)算機(jī)可視化技術(shù)與測(cè)繪技術(shù)的不斷發(fā)展，我國(guó)“數(shù)字化林業(yè)”建設(shè)進(jìn)程不斷向前推進(jìn)，地面三維激光掃描技術(shù)作為一種可以快速獲取目標(biāo)三維點(diǎn)云模型的新測(cè)繪技術(shù)，已在林業(yè)資源調(diào)查、林分結(jié)構(gòu)研究以及單木三維重建等方面得到了初步應(yīng)用[30]。王寧寧等利用三維掃描技術(shù)獲取了楊樹(shù)樹(shù)冠結(jié)構(gòu)特征參數(shù)，其精度符合要求[31]；王向玉等利用地面激光雷達(dá)點(diǎn)云數(shù)據(jù)重建了單木三維結(jié)構(gòu)[32]；謝鴻宇等利用地面激光雷達(dá)實(shí)現(xiàn)無(wú)損精確地獲取樹(shù)冠體積和表面積[33]；劉偉樂(lè)等利用地面三維掃描技術(shù)，運(yùn)用快速凸包算法提取樹(shù)木胸徑精度與其他方法相比精度最高[34]。三維激光掃描成像系統(tǒng)在獲取測(cè)樹(shù)因子的同時(shí)，不會(huì)對(duì)樹(shù)木造成任何損傷，便于連續(xù)觀測(cè)，進(jìn)行森林的動(dòng)態(tài)研究[35]。為此，在今后的研究中可以利用激光雷達(dá)獲取樹(shù)冠結(jié)構(gòu)參數(shù)，同時(shí)結(jié)合枝條葉面積模型實(shí)現(xiàn)無(wú)損精準(zhǔn)地估算樹(shù)冠葉面積以及林分的葉面積。

[1]呂 勇.馬尾松林葉面積的預(yù)測(cè)[J].湖北林業(yè)科技,1996,(3):38-41.

[2]尹 婧,邱國(guó)玉,何 凡,等.半干旱黃土丘陵區(qū)人工林葉面積特征[J].植物生態(tài)學(xué)報(bào),2008,32(2):440-447.

[3]刁 軍,國(guó) 紅,盧 軍,等.油松針葉面積估計(jì)模型及比葉面積的研究[J].林業(yè)科學(xué)研究,2013,26(2):174-180.

[4]Pietsch S A,Hasenauer H,Thornton P E. BGC-model parameters for tree species growing in central European forests[J]. Forest Ecology and Management, 2005,211(3):264-295.

[5]朱春全,雷靜品,劉曉東,等.不同經(jīng)營(yíng)方式下楊樹(shù)人工林葉面積分布與動(dòng)態(tài)研究[J].林業(yè)科學(xué),2001,37(1):47-51.

[6]趙 東,楊喜田,樊 巍,等.楊樹(shù)農(nóng)田防護(hù)林帶單木葉面積的變化[J].林業(yè)科學(xué),2011,47(4):107-113.

[7]Zellers C E, Saunders M R, Morrissey R C,et al.Development of allometric leaf area models for intensively managed black walnut(Juglans nigraL.)[J]. Annals of Forest Science, 2012,69(8):907-913.

[8]劉常富,何興元,陳 瑋,等.沈陽(yáng)城市森林三維綠量模擬及其影響因子[J].應(yīng)用生態(tài)學(xué)報(bào),2008,19(6):1173-1178.

[9]劉建偉,劉雅榮,朱春全,等.楊樹(shù)不同無(wú)性系葉面積模型的選擇及其建立[J].林業(yè)科學(xué),1994,30(6):481-485.

[10]吳春胤,張文昭,歐陽(yáng)慶,等.基于BP神經(jīng)網(wǎng)絡(luò)模型的荔枝樹(shù)葉面積測(cè)定方法[J].農(nóng)業(yè)工程學(xué)報(bào),2007,23(7):166-169.

[11]李軒然,劉琪璟,蔡 哲,等.千煙洲針葉林的比葉面積及葉面積指數(shù)[J].植物生態(tài)學(xué)報(bào),2007,31(1):93-101.

[12]董學(xué)軍.幾種沙生植物葉面積的估算與經(jīng)驗(yàn)公式[J].干旱區(qū)研究,1993,10(4):33-36.

[13]樊 巍,趙 東,高喜榮,等.‘中林46楊’農(nóng)田防護(hù)林單木葉面積特征[J].中國(guó)農(nóng)學(xué)通報(bào),2012, 28(4): 36-40.

[14]Bartelink H H. Allometric relationships for biomass and leaf area of beech[J].Ann Science Foresty,1997,(54):39-50.

[15]戚繼忠.枝徑與葉面積協(xié)調(diào)生長(zhǎng)關(guān)系的比較研究[J].南京林業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版,2004,28(1):92-94.

[16]郭孝玉.長(zhǎng)白落葉松人工林樹(shù)冠結(jié)構(gòu)及生長(zhǎng)模型研究[D].北京:北京林業(yè)大學(xué),2013.35-37.

[17]孟憲宇. 測(cè)樹(shù)學(xué)[M].北京:中國(guó)林業(yè)出版社,2006.187-188.

[18]劉常富,何興元,陳 瑋,等.沈陽(yáng)城市森林三維綠量測(cè)算[J].北京林業(yè)大學(xué)學(xué)報(bào),2006,28(3):32-37.

[19]Weiskittel A R, Kershaw J A, Hofmeyer P V,et al.Species differences in total and vertical distribution of branch-and tree-level leaf area for the five primary conifer species in Maine, USA[J].Forest Ecology and Management, 2009, 258(7): 1695-1703.

[20]Bartelink H H. Allometric relationships on biomass and needle area of Douglas- fir[J]. Forest Ecology and Management, 1996,86(1):193-203.

[21]劉志國(guó),蔡永立,李 愷.亞熱帶常綠闊葉林植物葉-小枝的異速生長(zhǎng)[J].植物生態(tài)學(xué)報(bào),2008,32(2):363-369.

[22]王東良,金荷仙,范麗琨,等.療養(yǎng)院人工綠地三維綠量分布特征及影響因子[J].浙江農(nóng)林大學(xué)學(xué)報(bào),2013,30(4):529-535.

[23]徐小軍,杜華強(qiáng),周?chē)?guó)模,等.基于遙感植被生物量估算模型自變量相關(guān)性分析綜述[J].遙感技術(shù)與應(yīng)用,2004,23(2):240-247.

[24]孫忠林,蔡雪穎,孫令娟,等. 黃花落葉松葉片生長(zhǎng)性狀[J].東北林業(yè)大學(xué)學(xué)報(bào),2011,39(12):24-26.

[25]Bond-Lamberty B,Wang C,Gower S T. Aboveground and belowground biomass and sapwood area allometric equations for six boreal tree species of northern Manitoba[J].Canadian Journal of Forest Research, 2002, 32 (8):1441-1450.

[26]陳 廈,桑衛(wèi)國(guó).暖溫帶地區(qū)3種森林群落葉面積指數(shù)和林冠開(kāi)闊度的季節(jié)動(dòng)態(tài)[J].植物生態(tài)學(xué)報(bào),2007,31(3):431-436.

[27]李 懷,付曉云,郭廷旭.朝陽(yáng)地區(qū)公路綠化樹(shù)種生態(tài)經(jīng)濟(jì)效益研究[J].西北林學(xué)院學(xué)報(bào),2007,22(4): 199-202.

[28]石鐵矛,潘續(xù)文,高 暢,等.城市綠地釋氧能力研究[J].沈陽(yáng)建筑大學(xué)學(xué)報(bào):自然科學(xué)版,2013,29(2):349-354.

[29]于 佳,陳宏偉,閆紅偉,等.沈陽(yáng)市常用園林植物碳匯功能研究[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2015,35(8):94-97.

[30]趙 陽(yáng),余新曉,信忠保,等.地面三維激光掃描技術(shù)在林業(yè)中的應(yīng)用與展望[J].世界林業(yè)研究,2010,23(4):41-45.

[31]王寧寧,尹文廣,黃秦軍,等.三維掃描技術(shù)在獲取楊樹(shù)樹(shù)冠結(jié)構(gòu)特征參數(shù)上的應(yīng)用[J].林業(yè)科學(xué), 2015,51(5):108-116.

[32]王向玉,謝東輝,汪 艷,等.基于地面激光雷達(dá)點(diǎn)云數(shù)據(jù)的單木三維重建[J].遙感技術(shù)與應(yīng)用,2015,30(3):455-460.

[33]謝鴻宇,趙耀龍,楊木壯,等.基于地面Lidar的樹(shù)冠體積和表面積測(cè)量方法研究[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2015,35(4): 1-6.

[34]劉偉樂(lè),林 輝,孫 華,等.基于地面三維激光掃描技術(shù)的樹(shù)木胸徑提取算法分析[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào),2014,34(11): 111-115.

[35]羅 旭,馮仲科,鄧向瑞,等.三維激光掃描成像系統(tǒng)在森林計(jì)測(cè)中的應(yīng)用[J].北京林業(yè)大學(xué)學(xué)報(bào),2007,29(增刊):82-87.

Leaf area simulation of branches for larch forests in the montane regions of eastern Liaoning province

JIA Bao-jun1, HU Jing-yang1, LIN Kuan1, FENG Qian-nan1, LIU Chang-fu1,2, YU Li-zhong3
(1.Forestry College, Shenyang Agricultural University, Shenyang 110866, Liaoning, China; 2. Research Institute of Forest Ecology,Environment and Protection, Chinese Academy of Forestry, Beijing 100091, China; 3.Qingyuan Forest CERN,Chinese Academy of Sciences, Shenyang 110016, Liaoning, China)

Based on the characteristics measured in montane region of eastern Liaoning province, 254 branches’ leaf area, basal branch diameter, branch length and any more branch attributes were obtained by branch analytical method. At branch level the models for leaf area and branch attributes were developed by unitary nonlinear regression and binary nonliner regression and multivariate nonlinear regression respectively. The results showed that the basal area of branch (d2) was highly correlated to branch leaf area. The modle ofy=8.967/(1+50.901e-0.084d) was the optimal unitary non-linearity regress equation. It’s coef ficient of determination (R2) was 0.719 and the precision was 86.34% (α=0.05). The coef ficient of determination(R2)was improved signi ficantly after we brought the variable of relative depth into crown (RDINC) into the equation. The modle ofy=0.002(d2)2.260e-1.701RDINCwas the optimal binary non-linearity regress equation. It’s coef ficient of determination (R2) was up to 0.796 and the precision was 88.57%(α=0.05).

Larixspp.; leaf area; branch level; regress model

S791.22

A

1673-923X(2016)10-0054-06

10.14067/j.cnki.1673-923x.2016.10.010

2015-12-20

引進(jìn)國(guó)際先進(jìn)林業(yè)科學(xué)技術(shù)（2013-4-59）；國(guó)家自然科學(xué)基金項(xiàng)目（31270490）

賈寶軍，碩士研究生

劉常富，教授；E-mail：liucf898@163.com

賈寶軍，胡靖揚(yáng)，林 寬，等. 遼東山區(qū)落葉松林枝條葉面積模型選擇[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2016, 36(10): 54-59.

[本文編校：吳 彬]