楸樹人工林生長(zhǎng)規(guī)律及其最優(yōu)生長(zhǎng)模型研究

關(guān)追追　趙江寧　邱權(quán)　麻文俊　王少明,馮小靜,蘇艷　王軍輝　李吉躍　何茜

摘 要：為更系統(tǒng)地研究楸樹的生長(zhǎng)規(guī)律，進(jìn)一步確定其速生期和最優(yōu)生長(zhǎng)模型，以河南省洛寧縣樹齡18～37年生的楸樹為研究對(duì)象，通過(guò)樹干解析獲取樹木的生長(zhǎng)數(shù)據(jù)，應(yīng)用系列模型對(duì)樹高、胸徑和材積生長(zhǎng)過(guò)程進(jìn)行擬合，篩選出楸樹最優(yōu)生長(zhǎng)模型。研究結(jié)果表明：①樹高總生長(zhǎng)量隨著樹齡的增加呈拋物線型增加形式，胸徑總生長(zhǎng)量呈S型曲線增加形式，材積總生長(zhǎng)量呈J型曲線增加形式;前15 a為樹高速生期，平均生長(zhǎng)速率為0.715 m/a，3～15 a為胸徑快速生長(zhǎng)期，此階段平均生長(zhǎng)速率為1.083 cm/a;材積連年生長(zhǎng)量與平均生長(zhǎng)量曲線在37 a時(shí)未相交，此時(shí)楸樹未達(dá)到數(shù)量成熟，不宜進(jìn)行采伐;②樹高的最優(yōu)生長(zhǎng)模型為修正Weibull與Mitscherlich模型，決定系數(shù)R2達(dá)到0.999;胸徑的最優(yōu)生長(zhǎng)模型為修正Weibull模型，R2達(dá)到0.995;材積的最優(yōu)生長(zhǎng)模型為Gompertz與Richards模型，R2達(dá)到0.998。將樹木生長(zhǎng)指標(biāo)實(shí)測(cè)值與模型擬合值進(jìn)行線性擬合，發(fā)現(xiàn)兩者相關(guān)性十分顯著（P<0.01），R2超過(guò)了0.990。從而得到結(jié)論，河南省洛寧縣楸樹的速生期是在造林后前15 a，輪伐期應(yīng)超過(guò)37 a，利用樹木生長(zhǎng)的最優(yōu)模型能夠準(zhǔn)確預(yù)測(cè)楸樹的生長(zhǎng)規(guī)律。

關(guān)鍵詞：楸樹;樹干解析;生長(zhǎng)模型;生長(zhǎng)規(guī)律;數(shù)量成熟

中圖分類號(hào)：S792??? 文獻(xiàn)標(biāo)識(shí)碼：A?? 文章編號(hào)：1006-8023（2021）02-0001-10

Growth Law and Optimal Growth Model of Catalpa bungei Plantation

——A Case Study of Catalpa bungei in Luoning County， Henan Province

GUAN Zhuizhui1， ZHAO Jiangning1， QIU Quan1， MA Wenjun2， WANG Shaoming3，

FENG Xiaojing4， SU Yan1， WANG Junhui2， LI Jiyue1， HE Qian1*

（1.Guangdong Key Laboratory for Innovative Development and Utilization of Forest Plant Germplasm（College of Forestry

and Landscape Architecture， South China Agricultural University）， Guangzhou 510642， China; 2.Key Laboratory of Tree

Breeding and Cultivation， National Forestry and Grassland Administration（Research Institute of Forestry， Chinese Academy of

Forestry），Beijing 100091， China;3.State-owned Luoning County， Quanbaoshan Forest Farm， He nan Luoyang 471000， China;

4.State-owned Luoning County， Lvcun Forest Farm， He nan Luoyang 471000， China）

Abstract： In order to systematically study the growth law of Catalpa bungei and determine its fast-growing period and the optimal growth model， the Catalpa bungei trees， aged from 18a to 37a， were selected as the research object in Luoning County of Henan Province. The growth data of trees was obtained through stem analysis， and a series of models were applied to fit growth process of tree height， diameter breast height （DBH） and volume， and the optimal growth model of Catalpa bungei was screened. The results showed： ①With the increase of tree age， the total growth of tree height increased in the form of parabola， the total growth of DBH increased in the form of S curve， and the total growth of volume increased in the form of J curve. The high-speed growth period of tree height appeared the first 15 a of tree growth， and the average growth rate was 0.715 m/a. The rapid growth stage of DBH appeared

收稿日期：2020-12-17

基金項(xiàng)目：國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（2017YFD060060404）

第一作者簡(jiǎn)介：關(guān)追追，博士研究生。研究方向?yàn)槿斯ち侄ㄏ蚺嘤?。E-mail： guanzhui9402 @163.com

通信作者：何茜，博士，副教授。研究方向?yàn)榱窒陆?jīng)濟(jì)。E-mail： heqian69@126.com

引文格式：關(guān)追追，趙江寧，邱權(quán)，等. 楸樹人工林生長(zhǎng)規(guī)律及其最優(yōu)生長(zhǎng)模型研究-以河南省洛寧縣楸樹為例[J]. 森林工程，2021，37（2）：1-10.

GUAN Z Z， ZHAO J N， QIU Q， et al. Growth law and optimal growth model of Catalpa bungei plantation-a case study of Catalpa bungei in Luoning County， Henan Province[J]. Forest Engineering，2021，37（2）：1-10.

among 3～15 a， and the average growth rate was 1.083 cm/a. The continuous annual and average growth curve of volume did not intersect when tree age reached 37 a. At this time， trees were not cut because they did not reach maturity. ②The optimal growth model of tree height was correct-Weibull and Mitscherlich model， and the determination coefficient R2 was 0.999. The optimal growth model of DBH was correct-Weibull model， R2 was 0.995. Gompertz and Richards models were the best growth models of volume， R2 was 0.998. The correlation between the measured value of tree growth index and the model fitting value was very significant （P<0.01）， R2 was more than 0.990. The conclusions are as follows， the fast-growing period of Catalpa bungei was in the first 15 years after afforestation， and the rotation period should be more than 37 years in Luoning County， Henan Province. Using the optimal model of tree growth can accurately predict the growth law of Catalpa bungei.

Keywords：Catalpa bungei; stem analysis; growth model; growth law; quantity maturity

0 引言

樹木生長(zhǎng)規(guī)律主要是指樹高、胸徑和材積等生長(zhǎng)因子隨樹齡的變化過(guò)程。了解樹木生長(zhǎng)規(guī)律，不但有助于認(rèn)識(shí)樹木自身生長(zhǎng)發(fā)育特點(diǎn)，而且對(duì)指導(dǎo)森林經(jīng)營(yíng)工作具有重要意義。目前，林麗平等[1]對(duì)44年生樟樹（Cinnamomum camphora）研究發(fā)現(xiàn)，材積平均生長(zhǎng)量和連年生長(zhǎng)量曲線在44 a時(shí)未相交，此時(shí)未達(dá)到數(shù)量成熟齡，舒馬切爾（Schumacher）、單分子式（Mitscherlich）以及坎派茲（Gompertz）模型是擬合樟樹樹高、胸徑和材積的最優(yōu)生長(zhǎng)模型;鐘連香等[2]對(duì)33年生廣東瓊楠（Beilschmiedia fordii）研究發(fā)現(xiàn)，6～12 a為廣東瓊楠的速生期，材積生長(zhǎng)達(dá)到數(shù)量成熟為30～32 a，樹高、胸徑和材積生長(zhǎng)的最優(yōu)模型分別為舒馬切爾、理查德（Richards）與坎派茲模型;劉小軍等[3]對(duì)麻櫟（Quercus acutissima）天然次生林研究發(fā)現(xiàn)，樹高的速生期為4～12 a，胸徑的速生期為10～42 a，材積的速生期為42～59 a，考爾夫（Korf）與坎派茲模型分別為樹高和胸徑的最優(yōu)生長(zhǎng)模型，邏輯斯蒂（Logistic）、坎派茲以及理查德模型為材積的最優(yōu)生長(zhǎng)模型;郝自遠(yuǎn)等[4]對(duì)24年生北美鵝掌楸（Liriodendron tulipifera）研究發(fā)現(xiàn)，樹高速生期為6～17 a;材積速生期為12～20 a，22 a后逐漸降低，林分接近數(shù)量成熟;陶玉華等[5]對(duì)15年生西南樺（Betula alnoides）研究發(fā)現(xiàn)，樹高和胸徑生長(zhǎng)量在第5年時(shí)達(dá)到高峰，材積連年生長(zhǎng)量與平均生長(zhǎng)量分別在第5年和第13年時(shí)達(dá)到高峰。材積連年生長(zhǎng)量與平均生長(zhǎng)量曲線在14 a時(shí)相交，此時(shí)西南樺達(dá)到數(shù)量成熟齡;梁有祥等[6]對(duì)火力楠（Michelia macclurei）研究發(fā)現(xiàn)，樹高速生期集中在前15 a，此期間出現(xiàn)2次明顯的生長(zhǎng)高峰期;胸徑速生期集中在前8 a;樹高與胸徑的最優(yōu)生長(zhǎng)模型為二次曲線模型。這些研究對(duì)揭示樹木速生期以及合理采伐期至關(guān)重要。

楸樹（Catalpa bungei）為紫葳科梓樹屬植物，為珍貴的園林綠化樹種。目前楸樹相關(guān)研究主要集中在良種選育[7]、抗旱節(jié)水評(píng)價(jià)[8-15]、配方施肥[16]以及水肥耦合[17-20]等方面，有關(guān)楸樹生長(zhǎng)規(guī)律的研究涉及較少。楸樹作為珍貴用材樹種，根據(jù)其不同生長(zhǎng)發(fā)育階段規(guī)律，結(jié)合相應(yīng)的栽培措施對(duì)其大徑材培育至關(guān)重要。因此，本研究以河南省洛寧縣樹齡18～37年生的楸樹為研究對(duì)象，通過(guò)樹干解析方法來(lái)探究樹木生長(zhǎng)規(guī)律，并篩選出最優(yōu)的生長(zhǎng)模型，可通過(guò)生長(zhǎng)模型對(duì)樹木發(fā)育過(guò)程進(jìn)行長(zhǎng)期預(yù)測(cè)，以期為楸樹人工林經(jīng)營(yíng)管理措施提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于河南省洛陽(yáng)市洛寧縣的馬店鄉(xiāng)楊村（34°25′09″ N、111°27′24″ E）和長(zhǎng)水鄉(xiāng)連山村（34°24′14″ N、111°22′24″ E）。洛寧縣（34°05′～34°38′ N、111°08′～111°49′ E）屬于伊洛河流域，總面積約為2 300 km2。地形復(fù)雜，有山地、丘陵和平原等多種類型。該地氣候適宜，四季分明，年均氣溫約為14 ℃，年均降水量為600～800 mm。土壤包括黃土、棕壤等多種類型，自然植被主要為暖溫帶落葉闊葉樹種。

1.2 解析木和圓盤獲取

在2019年楸樹生長(zhǎng)季末期，在試驗(yàn)地隨機(jī)設(shè)置7塊樣地，每塊樣地面積為400 m2。試驗(yàn)地現(xiàn)存密度約為500株/hm2。在每塊樣地中隨機(jī)選取1株平均木作為解析木，共獲取7株解析木，進(jìn)行樹高、胸徑和材積建模。將樣地旁隨機(jī)獲取的5株楸樹生長(zhǎng)數(shù)據(jù)對(duì)模型進(jìn)行檢驗(yàn)。所選樣木應(yīng)具有完整的樹冠，樹干無(wú)分叉、無(wú)傾斜。樣木樹齡為18～37 a，樹高為11.50～19.30 m，胸徑為16.45～33.20 cm。其中，8株樣木來(lái)自純林，其他4株樣木來(lái)自農(nóng)林間作的林分，主要間作農(nóng)作物是小麥和玉米。兩種林型條件下土壤理化性質(zhì)請(qǐng)參見表1。

樣木伐倒前，先確定根頸和胸徑位置，同時(shí)測(cè)定胸徑，在樹干上標(biāo)記北向位置。伐倒后，測(cè)定樹干長(zhǎng)度，根據(jù)中央斷面積區(qū)分求積法，在根頸處（0 m）、胸徑處（1.3 、3.6 、5.6 、7.6 m，…）和樹梢（不足1個(gè)區(qū)分段）底部處標(biāo)記需要截取的圓盤位置。圓盤截取時(shí)油鋸與樹干保持垂直，不可偏斜，圓盤厚度一般為3～5 cm。對(duì)圓盤進(jìn)行編號(hào)，并標(biāo)記北向和斷面高度，具有標(biāo)記的一面為非工作面，另一面為工作面。

在實(shí)驗(yàn)室用砂光機(jī)對(duì)圓盤工作面進(jìn)行打磨，保證年輪清晰，易于觀察。圓盤測(cè)定前，用鉛筆在工作面上沿東西、南北兩個(gè)方向標(biāo)記2條直線（需經(jīng)過(guò)圓盤髓心），用大頭針標(biāo)記出每個(gè)年輪與直線的交點(diǎn)，隨后查數(shù)年輪數(shù)量。用直尺測(cè)定圓盤的東西、南北方向的帶皮直徑以及每個(gè)年輪自東、西、南、北4個(gè)方向至髓心的距離（半徑），計(jì)算各年輪的平均寬度。

1.3 研究方法

根據(jù)孟憲宇[21]的方法，將楸樹生長(zhǎng)以3 a劃分為一個(gè)齡階，擬合出各齡階下樹高和胸徑總生長(zhǎng)量，單株材積采用中央斷面區(qū)分求積公式計(jì)算：

V=∑ni=1gi×li+13g0×l0。

式中：V為單株材積，m3;gi為第i區(qū)分段中央斷面積，m2;li為第i區(qū)分段長(zhǎng)度，m;g0為梢頭底端斷面積，m2;l0為梢頭長(zhǎng)度，m。

楸樹各調(diào)查因子的總平均生長(zhǎng)量計(jì)算公式為：

θ=Vtt。

式中：θ為總平均生長(zhǎng)量;Vt為t年時(shí)樹木調(diào)查因子的總生長(zhǎng)量;t為調(diào)查時(shí)的樹木年齡。

楸樹各調(diào)查因子的連年生長(zhǎng)量采用定期平均生長(zhǎng)量代替，計(jì)算公式為：

θn=Vt-Vt-nn。

式中：θn為定期平均生長(zhǎng)量;Vt為t年時(shí)樹木調(diào)查因子的總生長(zhǎng)量;Vt-n為n年前樹木調(diào)查因子的總生長(zhǎng)量;n為間隔期的年數(shù)。

1.4 樹木生長(zhǎng)模型擬合

參考孟憲宇[21]的方法，選取舒馬切爾（Schumacher）和修正Weibull等2種生長(zhǎng)經(jīng)驗(yàn)方程以及邏輯斯蒂（Logistic）、單分子式（Mitscherlich）、坎派茲（Gompertz）和理查德（Richards）等4種生長(zhǎng)理論方程對(duì)楸樹樹高和胸徑生長(zhǎng)過(guò)程進(jìn)行擬合。選取舒馬切爾、邏輯斯蒂、坎派茲以及理查德等4種方程對(duì)材積生長(zhǎng)過(guò)程進(jìn)行擬合。具體模型如下。

（1）舒馬切爾模型 y=aexp（-b/t）。

（2）修正Weibull模型 y=a[1-exp（-btc）]。

（3）邏輯斯蒂模型 y=A/[1+mexp（-rt）]。

（4）單分子式模型 y=A[1-exp（-rt）]。

（5）坎派茲模型 y=Aexp[-bexp（-rt）]。

（6）理查德模型 y=A[1-exp（-rt）]c。

模型中：y分別為樹高、胸徑和材積;t為樹齡，a、b、c、m、r、A均為模型參數(shù)。

1.5 模型檢驗(yàn)

為了檢驗(yàn)樹木生長(zhǎng)方程的擬合效果，采用殘差平方和（SSR，公式中用SSSR表示）、均方差（MSE，公式中用MMSE表示）、均方根（RMSE，公式中用RRMSE表示）以及決定系數(shù)（R2）對(duì)模型的擬合效果進(jìn)行評(píng)價(jià)，模型評(píng)價(jià)公式為：

（1）SSSR=∑ni=1（yi-y）2。

（2）MMSE=1n∑ni=1（yi-y）2。

（3）RRMSE=1n∑ni=1（yi-y）2。

（4）R2=1-∑ni=1（yi-y）2∑ni=1（yi-y-）2。

式中：y是樣本某因子的實(shí)測(cè)值;yi是y的無(wú)偏估計(jì)值;y-是y的平均值;n是樣本量。

通常認(rèn)為，SSR、MSE和RMSE越小，R2的值越接近1，模型擬合效果越好。

1.6 統(tǒng)計(jì)分析

采用Microsoft Office Excel 2010軟件對(duì)樹木生長(zhǎng)數(shù)據(jù)進(jìn)行分類整理，采用SPSS 20.0軟件進(jìn)行統(tǒng)計(jì)分析，采用Origin2018軟件進(jìn)行制圖。

2 結(jié)果與分析

2.1 樹高生長(zhǎng)隨樹齡的變化

楸樹高總生長(zhǎng)量隨著樹齡的增加呈現(xiàn)拋物線型增加形式，即樹高先快速增加，達(dá)到一定樹齡后，樹高增加速度開始減慢。具體而言，前15 a是樹高的速生期，平均生長(zhǎng)速率為0.715 m/a，隨后樹高生長(zhǎng)進(jìn)入緩慢生長(zhǎng)期，當(dāng)樹齡超過(guò)33年生時(shí)，樹高總生長(zhǎng)量基本不再增加，如圖1所示。樹高平均生長(zhǎng)量隨著樹齡的增加先迅速增加后逐漸減少，平均生長(zhǎng)量范圍為0.48～0.89 m。樹高連年生長(zhǎng)量總體上隨著樹齡的增加先快速增加后波動(dòng)式降低，尤其在3～21 a，連年生長(zhǎng)量下降十分明顯;在21～37 a，連年生長(zhǎng)量波動(dòng)幅度較小，為0.31～0.37 m。此外，樹高連年生長(zhǎng)量和平均生長(zhǎng)量曲線的交點(diǎn)出現(xiàn)在3～6 a，6 a后平均生長(zhǎng)量高于連年生長(zhǎng)量，如圖2所示。

2.2 胸徑生長(zhǎng)隨樹齡的變化

楸樹胸徑總生長(zhǎng)量隨著樹齡增加總體上呈現(xiàn)S型曲線增加形式，即“緩慢生長(zhǎng)→快速生長(zhǎng)→緩慢生長(zhǎng)”。具體而言，前3 a為緩慢生長(zhǎng)期，平均生長(zhǎng)速率為0.263 cm/a;3～15 a為快速生長(zhǎng)期，平均生長(zhǎng)速率為1.083 cm/a;15 a后生長(zhǎng)速度逐漸減慢，如圖3所示。胸徑平均生長(zhǎng)量隨著樹齡的增加總體上先迅速增加后緩慢降低，在1～12 a，平均生長(zhǎng)量呈現(xiàn)S型曲線增加形式，在12 a時(shí)達(dá)到最大為0.93 cm，隨后平均生長(zhǎng)量緩慢減少。胸徑連年生長(zhǎng)量隨著樹齡的增加總體上呈現(xiàn)“升→降”的趨勢(shì)，其中在3～6 a，連年生長(zhǎng)量增加最快，增加了0.84 cm;在12～15 a，連年生長(zhǎng)量下降最快，降低了0.26 cm;在21～37 a，連年生長(zhǎng)量出現(xiàn)了波動(dòng)現(xiàn)象，變化幅度為0.43～0.73 cm。胸徑連年生長(zhǎng)量和平均生長(zhǎng)量曲線相交于12～15 a，前15 a連年生長(zhǎng)量高于平均生長(zhǎng)量，15 a后連年生長(zhǎng)量低于平均生長(zhǎng)量，平均生長(zhǎng)量緩慢降低，而連年生長(zhǎng)量相對(duì)下降較快，在37 a胸徑連年生長(zhǎng)量低至0.43 cm，如圖4所示。

2.3 材積生長(zhǎng)隨樹齡的變化

楸樹材積總生長(zhǎng)量隨著樹齡增加總體上呈現(xiàn)J型曲線增加形式，即“緩慢生長(zhǎng)→快速生長(zhǎng)”。具體表現(xiàn)為，前15 a材積生長(zhǎng)量增加較緩慢，15 a后材積生長(zhǎng)量迅速增加，在37 a時(shí)達(dá)到0.521 m3，如圖5所示。材積平均生長(zhǎng)量隨著樹齡增加總體上呈現(xiàn)線性增加趨勢(shì)，在37 a時(shí)達(dá)到0.014 m3。材積連年生長(zhǎng)量隨著樹齡增加呈現(xiàn)波動(dòng)式上升的趨勢(shì)，在36 a時(shí)達(dá)到0.027 m3?？傮w而言，各樹齡材積連年生長(zhǎng)量高于平均生長(zhǎng)量，如圖6所示。

2.4 生長(zhǎng)模型的擬合與篩選

根據(jù)最優(yōu)生長(zhǎng)模型的篩選條件[22-27]，發(fā)現(xiàn)修正Weibull、Mitscherlich以及Richards模型的R2達(dá)到了0.999，但Richards模型的RSS、MSE與RMSE均高于其他2種模型，見表2。雖然6種模型的R2均超過(guò)了0.970，但修正Weibull與Mitscherlich模型的RSS、MSE與RMSE的值最小，R2卻最大。故適合樹高生長(zhǎng)的最優(yōu)模型為修正Weibull與Mitscherlich模型。

從胸徑生長(zhǎng)方程擬合結(jié)果來(lái)看（表3），6種生長(zhǎng)模型的R2為0.975～0.995。其中，修正Weibull模型擬合效果最好，其次是Richards模型，Logistic模型擬合效果相對(duì)較差。因此，修正Weibull模型是預(yù)測(cè)楸樹胸徑生長(zhǎng)的最優(yōu)模型。

由表4可知，Gompertz與Richards模型對(duì)楸樹材積生長(zhǎng)過(guò)程擬合效果最好，R2達(dá)到了0.998，RSS、MSE與RMSE在4種模型中均最低。Schumacher與Logistic模型擬合效果相對(duì)較差。

2.5 模型的檢驗(yàn)

模型的檢驗(yàn)采用楸樹樹高、胸徑和材積最優(yōu)生長(zhǎng)模型的擬合值與實(shí)測(cè)值進(jìn)行線性擬合，最后利用RSS、MSE、RMSE與R2進(jìn)行預(yù)測(cè)精度比較。結(jié)果表明，與楸樹生長(zhǎng)指標(biāo)實(shí)際測(cè)量值相比，模型擬合的樹高、胸徑和材積結(jié)果均十分顯著（P<0.01）（表5和圖7）。樹高、胸徑和材積的R2分別達(dá)到了0.999 1、0.994 7和0.997 8。

3 討論

樹木的生長(zhǎng)會(huì)隨著樹齡的增加而發(fā)生變化。本研究中，楸樹樹高生長(zhǎng)表現(xiàn)為前期生長(zhǎng)快，后期生長(zhǎng)慢，這與木荷（Schima superba）[28]、云南松（Pinus yunnanensis）[29]以及濕地松（Pinus elliottii）[30]的生長(zhǎng)規(guī)律類似。具體而言，楸樹高在前15 a生長(zhǎng)速度較快，15 a后生長(zhǎng)速度逐漸降低。由此說(shuō)明，楸樹造林后前15 a是其樹高的速生期。楸樹胸徑生長(zhǎng)總體上呈現(xiàn)S型曲線增加形式，這與臭松（Abies nephrolepis）[31]和頂果木（Acrocarpus fraxinifolius）[32]的研究結(jié)果類似。具體而言，前3 a為胸徑的緩慢生長(zhǎng)期，3～15 a為快速生長(zhǎng)期，15 a后生長(zhǎng)速度逐漸減慢。由此說(shuō)明，3～15 a是楸樹胸徑的速生期，此階段做好施肥撫育工作，可以提供樹木快速生長(zhǎng)所需的養(yǎng)分。因此，本研究認(rèn)為造林后前15 a楸樹正處于幼齡時(shí)期，此階段胸徑連年生長(zhǎng)量與平均生長(zhǎng)量均隨著樹齡的增加而增加，但連年生長(zhǎng)量增加的速度較快，高于平均生長(zhǎng)量，且連年生長(zhǎng)量達(dá)到最高峰的時(shí)間比平均生長(zhǎng)量早，但在平均生長(zhǎng)量達(dá)到最高峰后，連年生長(zhǎng)量小于平均生長(zhǎng)量。造林15 a后楸樹進(jìn)入中齡林階段，樹高和胸徑平均生長(zhǎng)量呈現(xiàn)出線性下降趨勢(shì)，可能是由于樹種間對(duì)溫度、光照和水分的競(jìng)爭(zhēng)加劇，從而降低了樹高和胸徑平均生長(zhǎng)量。此階段可通過(guò)一些營(yíng)林措施，如撫育間伐、施肥或修枝等措施來(lái)調(diào)節(jié)樹木生長(zhǎng)。

值得注意的是，楸樹樹高和胸徑連年生長(zhǎng)量在降低過(guò)程中出現(xiàn)了波動(dòng)現(xiàn)象（21～37 a），可能是年際氣候變化或人為干擾等因素造成的，具體原因還有待進(jìn)一步研究。楸樹材積總生長(zhǎng)量隨著樹齡增加總體上呈現(xiàn)J型曲線增加形式，與華榛（Corylus chinensis）[33]的材積生長(zhǎng)變化規(guī)律類似，說(shuō)明在樹木發(fā)育前期，材積生長(zhǎng)緩慢，但隨著高度和直徑生長(zhǎng)進(jìn)入速生期，材積生長(zhǎng)也迅速增加。隨著樹木生長(zhǎng)逐漸步入老齡階段，材積生長(zhǎng)會(huì)逐漸降低最后趨于平穩(wěn)狀態(tài)，但本研究中在樹齡達(dá)到37 a前，材積一直處于不斷增加的過(guò)程，說(shuō)明此時(shí)楸樹生長(zhǎng)未步入老齡階段。

此外，研究發(fā)現(xiàn)若材積平均生長(zhǎng)量與連年生長(zhǎng)量曲線相交，則表明樹木材積已達(dá)到數(shù)量成熟。本研究中楸樹材積連年生長(zhǎng)量與平均生長(zhǎng)量曲線在37 a時(shí)仍未相交，說(shuō)明在洛寧縣地區(qū)楸樹生長(zhǎng)緩慢，在37 a時(shí)還未達(dá)到數(shù)量成熟，林木材積仍處于速生期。與其他樹種如禿杉（Taiwania flousiana）（數(shù)量成熟齡為35 a）[34]和巒大杉（Cunninghamia konishii）（數(shù)量成熟齡為32 a）[35]相比，楸樹數(shù)量成熟齡更高，說(shuō)明楸樹輪伐期應(yīng)超過(guò)37 a，該樹種的數(shù)量成熟齡有待進(jìn)一步研究。

樹木生長(zhǎng)方程作為描述林木大小隨年齡變化的模型，可以反映樹木生長(zhǎng)的規(guī)律性。本研究發(fā)現(xiàn)修正Weibull和Mitscherlich模型對(duì)楸樹樹高的模擬效果最好，R2達(dá)到了0.999;最適合胸徑生長(zhǎng)的模型為修正Weibull模型，R2達(dá)到了0.995;Gompertz與Richards模型是預(yù)測(cè)材積生長(zhǎng)的最優(yōu)模型，R2達(dá)到了0.998。模型的檢驗(yàn)結(jié)果表明，與實(shí)測(cè)值相比，模型擬合的楸樹樹高、胸徑和材積結(jié)果均十分顯著，R2均超過(guò)了0.990。

以上篩選的樹木最優(yōu)生長(zhǎng)模型可以重塑河南省洛寧縣地區(qū)楸樹動(dòng)態(tài)發(fā)育過(guò)程，由于本研究所選樣本量和立地類型的限制，建立的楸樹最優(yōu)生長(zhǎng)模型僅適用于洛寧縣楸樹人工林。今后應(yīng)增加試驗(yàn)樣本量，綜合多種立地類型，從而更加系統(tǒng)地建立與完善楸樹人工林最優(yōu)生長(zhǎng)模型，以期為楸樹撫育措施的制定提供參考依據(jù)。

4 結(jié)論

綜合河南省洛寧縣楸樹測(cè)樹因子生長(zhǎng)量的研究，發(fā)現(xiàn)該樹種在前15 a具有明顯的速生特性。通過(guò)系列方程擬合，認(rèn)為楸樹樹高的最優(yōu)生長(zhǎng)模型為修正Weibull和Mitscherlich模型;胸徑的最優(yōu)生長(zhǎng)模型為修正Weibull模型;材積的最優(yōu)生長(zhǎng)模型為Gompertz與Richards模型。楸樹材積的平均生長(zhǎng)量與連年生長(zhǎng)量曲線在37a時(shí)未相交，此時(shí)未達(dá)到數(shù)量成熟齡，因此不宜進(jìn)行采伐。通過(guò)對(duì)河南省洛寧縣地區(qū)楸樹生長(zhǎng)規(guī)律的研究，了解該樹種在不同年齡下的生長(zhǎng)規(guī)律從而進(jìn)行科學(xué)的撫育管理，對(duì)這一樹種未來(lái)的發(fā)展具有重要意義。

【參考文獻(xiàn)】

[1]林麗平，徐期瑚，羅勇，等.廣東省樟樹立木生長(zhǎng)規(guī)律和生長(zhǎng)模型研究[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2018，38（6）：23-29.

LIN L P， XU Q H， LUO Y， et al. Study on the growth laws and models of Cinnamomum camphora in Guangdong province[J]. Journal of Central South University of Forestry & Technology， 2018， 38（6）：23-29.

[2]鐘連香，米偉，施福軍.廣東瓊楠人工林生長(zhǎng)規(guī)律研究[J].西南農(nóng)業(yè)學(xué)報(bào)，2018，31（12）：2666-2671.

ZHONG L X， MI W， SHI F J. Research on growth regularity of Beilschmiedia fordii plantation[J]. Southwest China Journal of Agricultural Sciences， 2018， 31（12）：2666-2671.

[3]劉小軍，劉平宇，趙苗苗，等.天然麻櫟次生林生長(zhǎng)規(guī)律研究[J].西北林學(xué)院學(xué)報(bào)，2019，34（1）：175-179.

LIU X J， LIU P Y， ZHAO M M， et al. Growth law of natural Quercus acutissima secondary forest[J]. Journal of Northwest Forestry University， 2019， 34（1）：175-179.

[4]郝自遠(yuǎn)， 李火根， 康昊， 等. 北美鵝掌楸人工林生長(zhǎng)規(guī)律及早期選擇可行性探究[J]. 林業(yè)科學(xué)研究， 2017， 30（5）： 878-885.

HAO Z Y， LI H G， KANG H， et al. Growth pattern and early selection of Liriodendron tulipifera[J]. Forest Research， 2017， 30（5）： 878-885.

[5]陶玉華， 趙峰， 龐正轟， 等. 西南樺人工林生長(zhǎng)規(guī)律研究[J]. 西北林學(xué)院學(xué)報(bào)， 2019， 34（1）： 180-183.

TAO Y H， ZHAO F， PANG Z H， et al. Growth law of Betula alnoides plantation[J]. Journal of Northwest Forestry University， 2019， 34（1）： 180-183.

[6]梁有祥，秦武明，玉桂成，等.桂東南地區(qū)火力楠人工林生長(zhǎng)規(guī)律研究[J].西北林學(xué)院學(xué)報(bào)，2011，26（2）：150-154.

LIANG Y X， QIN W M， YU G C， et al. Growth regularity of Michelia macclurei plantation in southeastern of Guangxi[J]. Journal of Northwest Forestry University， 2011， 26（2）：150-154.

[7]王秋霞.‘中林5號(hào)楸樹優(yōu)良品種選育研究[J].河南林業(yè)科技，2019，39（2）：11-14.

WANG Q X. Study on breeding of the excellent seedling tree of Zhonglin Sorbus No. 5[J]. Journal of Henan Forestry Science and Technology， 2019， 39（2）：11-14.

[8]董蕾， 陳博， 李吉躍， 等. 5個(gè)楸樹無(wú)性系抗旱性的綜合評(píng)價(jià)[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào)， 2014， 34（4）： 30-35.

DONG L， CHEN B， LI J Y， et al. Drought resistance evaluation of 5 Catalpa bungei clones[J]. Journal of Central South University of Forestry & Technology， 2014， 34（4）： 30-35.

[9]董蕾， 李吉躍， 王軍輝， 等. 干旱脅迫下楸樹無(wú)性系水分生理表征關(guān)系初探[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào)， 2015， 35（10）： 26-32.

DONG L， LI J Y， WANG J H， et al. Comparation of water consumption and transpiration of 3 Catalpa bungei clones under drought[J]. Journal of Central South University of Forestry & Technology， 2015， 35（10）： 26-32.

[10]董蕾， 李吉躍， 謝堃， 等. 5個(gè)楸樹（Catalpa bungei）無(wú)性系耗水特性比較[J]. 山東農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2016， 47（1）： 31-36.

DONG L， LI J Y， XIE K， et al. Comparison among characteristics of water consumption in 5 clones of Catalpa bungei[J]. Journal of Shandong Agricultural University （Natural Science Edition）， 2016， 47（1）： 31-36.

[11]QIU Q， WANG J H， SU Y， et al. Organ-level evaluation of the carbon starvation hypothesis in deciduous broad-leaved Catalpa bungei plants undergoing drought-induced mortality[J]. Dendrobiology， 2018， 80： 48-60.

[12]SHI H L， MA W J， SONG J Y， et al. Physiological and transcriptional responses of Catalpa bungei to drought stress under sufficient and deficient-nitrogen conditions[J]. Tree Physiology， 2017， 37（11）： 1457-1468.

[13]LUO Z B， LUO J. Uncovering the physiological mechanisms that allow nitrogen availability to affect drought acclimation in Catalpa bungei[J]. Tree Physiology， 2017， 37（11）： 1453-1456.

[14]ZHENG H F， ZHANG X， MA W J， et al. Morphological and physiological responses to cyclic drought in two contrasting genotypes of Catalpa bungei[J]. Environmental and Experimental Botany， 2017， 138： 77-87.

[15]WU J W， LI J Y， SU Y， et al. A morphophysiological analysis of the effects of drought and shade on Catalpa bungei plantlets[J]. Acta Physiologiae Plantarum， 2017， 39（3）： 80.

[16]白晶晶， 吳俊文， 何茜， 等. 不同配方施肥對(duì)楸樹幼苗生物量分配及養(yǎng)分利用的影響[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報(bào)， 2015， 36（6）： 91-97.

BAI J J， WU J W， HE Q， et al. Effects of different fertilization formulas on Catalpa bungei seedling biomass allocation and nutrient use efficiency[J]. Journal of South China Agricultural University， 2015， 36（6）： 91-97.

[17]邱權(quán)， 李吉躍， 王軍輝， 等. 水肥耦合效應(yīng)對(duì)楸樹苗木葉片δ13C的影響[J]. 西北植物學(xué)報(bào)， 2015， 35（11）： 2290-2298.

QIU Q， LI J Y， WANG J H， et al. Coupling effects of soil water and fertilizer application on leaf δ13C of Catalpa bungei seedlings[J]. Acta Botanica Boreali-Occidentalia Sinica， 2015， 35（11）： 2290-2298.

[18]邱權(quán)， 李吉躍， 王軍輝， 等. 水肥耦合效應(yīng)對(duì)楸樹苗期葉片凈光合速率和SPAD值的影響[J]. 生態(tài)學(xué)報(bào)， 2016， 36（11）： 3459-3468.

QIU Q， LI J Y， WANG J H， et al. Interactive effects of soil water and fertilizer application on leaf net photosynthetic rate and SPAD readings of Catalpa bungei seedlings[J]. Acta Ecologica Sinica， 2016， 36（11）： 3459-3468.

[19]QIU Q， HE Q， WANG J H， et al. Interactions between soil water and fertilizer on leaf photosynthesis and δ13C of Catalpa bungei seedlings[J]. Dendrobiology， 2016， 76： 25-34.

[20]QIU Q， LI J Y， WANG J H， et al. Coupling effects of water and fertilizer on the growth characteristics of Catalpa bungei seedlings[J]. Pakistan Journal of Botany， 2015， 47（3）： 889-896.

[21]孟憲宇. 測(cè)樹學(xué)[M]. 北京： 中國(guó)林業(yè)出版社， 2006.

MENG X Y. Dendrometry[ [M]. Beijing： China Forestry Publishing House， 2006.

[22]李貴， 童方平， 劉振華， 等. 黧蒴栲生長(zhǎng)過(guò)程及數(shù)量成熟、工藝成熟的初步研究[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào)， 2013， 33（12）： 53-56.

LI G， TONG F P， LIU Z H， et al. Preliminary study on growth process， quantitative maturity and technical maturity of Castanopsis fissa[J]. Journal of Central South University of Forestry & Technology， 2013， 33（12）： 53-56.

[23]白文斌， 廖超英， 康樂(lè)， 等. 西藏昌都地區(qū)川西云杉林木生長(zhǎng)規(guī)律研究[J]. 西北林學(xué)院學(xué)報(bào)， 2012， 27（5）： 158-162.

BAI W B， LIAO C Y， KANG L， et al. Growth regularity of Picea likiangensis var. balfouriana in Changdu of Tibet[J]. Journal of Northwest Forestry University， 2012， 27（5）： 158-162.

[24]韋善華， 覃靜， 朱賢良， 等. 南寧地區(qū)灰木蓮人工林生長(zhǎng)規(guī)律研究[J]. 西北林學(xué)院學(xué)報(bào)， 2011， 26（5）： 174-178.

WEI S H， QIN J， ZHU X L， et al. Growth regularity of Manglietia glauce plantation growing in Nanning[J]. Journal of Northwest Forestry University， 2011， 26（5）： 174-178.

[25]陳英， 楊華， 李偉， 等. 北京地區(qū)側(cè)柏人工林標(biāo)準(zhǔn)木生長(zhǎng)過(guò)程研究[J]. 西北林學(xué)院學(xué)報(bào)， 2012， 27（5）： 153-157.

CHEN Y， YANG H， LI W， et al. Standard tree growth of Platycladus orientalis plantation in Beijing[J]. Journal of Northwest Forestry University， 2012， 27（5）： 153-157.

[26]何遠(yuǎn)洋， 段劼， 魏永忠， 等. 遼東山區(qū)長(zhǎng)白落葉松人工林生長(zhǎng)規(guī)律的研究[J]. 西北林學(xué)院學(xué)報(bào)， 2014， 29（2）： 184-190.

HE Y Y， DUAN J， WEI Y Z， et al. Growth regularity of Larix olgensis plantation in mountainous regions of eastern Liaoning Province， China[J]. Journal of Northwest Forestry University， 2014， 29（2）： 184-190.

[27]李鵬飛，高曉陽(yáng)，李紅嶺，等.基于3D掃描溫室葡萄葉片生長(zhǎng)建模及可視化[J].林業(yè)機(jī)械與木工設(shè)備，2020，48（2）：20-25.

LI P F， GAO X Y， LI H L， et al. Modeling and visualization of grape leaf growth in greenhouse based on 3D scanning[J]. Forestry Machinery & Woodworking Equipment， 2020， 48（2）：20-25.

[28]韋昌幸. 48年生木荷人工林生產(chǎn)力及經(jīng)濟(jì)效益評(píng)價(jià)研究[D]. 長(zhǎng)沙： 中南林業(yè)科技大學(xué)， 2014.

WEI C X. Productivity analysis and economic benefit evaluation of 48a Schima superba plantation[D]. Changsha： Central South University of Forestry & Technology， 2014.

[29]束傳林， 廖聲熙， 李昆， 等. 一平浪林場(chǎng)云南松單株木生長(zhǎng)規(guī)律初步研究[J]. 陜西林業(yè)科技， 2008，36（3）： 1-3.

SHU C L， LIAO S X， LI K， et al. Preliminary study on the growth of single plant of Pinus yunnanensis in Yipinglang forest farm[J]. Shaanxi Forest Science and Technology， 2008， 36（3）： 1-3.

[30]周義罡， 王瑞輝， 謝裕宏， 等. 架空線下速生超高樹木樹高生長(zhǎng)規(guī)律研究[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào)， 2016， 36（2）： 75-78.

ZHOU Y G， WANG R H， XIE Y H， et al. Tree height growth rhythm of fast-growing ultra-high trees under overhead line[J]. Journal of Central South University of Forestry & Technology， 2016， 36（2）： 75-78.

[31]劉玉. 臭松次生林結(jié)構(gòu)規(guī)律和生長(zhǎng)模擬研究[D]. 長(zhǎng)沙： 中南林業(yè)科技大學(xué)， 2015.

LIU Y. Research on stand structure and growth of Adiessibirica ledeb[D]. Changsha： Central South University of Forestry & Technology， 2015.

[32]呂曼芳， 梁乃鵬， 秦武明， 等. 頂果木人工林生長(zhǎng)規(guī)律的研究[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào)， 2013， 33（8）： 43-49.

LV M F， LIANG N P， QIN W M， et al. Study on growth rhythm of Acrocarpus fraxinifolius plantation[J]. Journal of Central South University of Forestry & Technology， 2013， 33（8）： 43-49.

[33]楊亞莉. 華榛人工林生長(zhǎng)規(guī)律及開花結(jié)實(shí)特性研究[D]. 長(zhǎng)沙： 中南林業(yè)科技大學(xué)， 2017.

YANG Y L. Growth pattern， blossoms and solid characteristic of Corylus chinensis plantation[D]. Changsha： Central South University of Forestry & Technology， 2017.

[34]孫昀. 禿杉人工林生長(zhǎng)規(guī)律及生產(chǎn)力研究[D]. 長(zhǎng)沙： 中南林業(yè)科技大學(xué)， 2016.

SUN Y. The study on growth rule and productivity of Taiwania flousiana[D]. Changsha： Central South University of Forestry & Technology， 2016.

[35]歐建德， 吳志莊， 康永武. 福建省巒大杉人工林生長(zhǎng)模擬及其數(shù)量成熟研究[J]. 東北林業(yè)大學(xué)學(xué)報(bào)， 2018， 46（4）： 18-22.

OU J D， WU Z Z， KANG Y W. Growth simulation and quantity maturity of Cun ninghamia konishii plantation in Fujian Province[J]. Journal of Northeast Forestry University， 2018， 46（4）： 18-22.