柳树良种苗期年生长模型的探讨.docx

1、柳树良种苗期年生长模型的探讨柳树良种苗期年生长模型的探讨摘 要：以柳树良种品种为材料，利用Logistic方程对各个品种（系）的1a生扦插苗木的苗高、地径的年生长节律进行拟合。根据因素数列的几何形状发展态势的接近程度，将苗木的生长量与影响苗木生长的空气温度、空气相对湿度、降水量和日照时数组成一个灰色系统，进行关联度分析。结果表明：利用Logistic方程拟合效果显著，相关系数皆在0.90以上。并以此将柳树良种品种（系）扦插苗的生长过程分为成活期、生长前期、速生期、和生长后期4个时期，年生长量表现出明显的“慢快慢”的节律，其中速生期内的苗高生长量占总生长量的比率平均为60.5，地径生长量占总生长

2、量的比率平均为64.0。关联度分析表明空气温度对苗木的高生长影响最大，降水量影响最小；而日照时数对苗木的粗生长影响最大，降水量影响最小。通过对苗木年生长节律和影响苗木生长的主要环境因子进行分析，为加强柳树良种苗木生产管理，提高苗木的质量和土地的生产力，提供了理论依据。关键词：柳树良种；苗期；Logistic方程；生长模型；灰色关联度Study on Annual Growth Model of willow excellent fineseed at Nursery StageAbstract: Logistic equation is applied to model growth incr

3、ement of willow excellent clone fineseed at nursery stage. Influence of climatic factors such as air temperature, relative humidity, precipitation and illuminating Hours on the seedling growth are analyzed. The results show that correlation of the modeling is significant And the correlation coeffici

4、ents are more than 0.90. Based on the equation, annual growth of the seedlings Can be divided into four periods, namely survival period, pre-growth period, fast-growth period and post-growth period. Average height and ground diameter in fast-growth period is 60.5%and 64.0% respectively of the total

5、annual. Air temperature is believed to have the most significant influence on the height growth and precipitation the least, whilst illuminating hours have the most significant influence on the ground diameter growth and precipitation the least. This paper has provided theoretical basis for the prod

6、uction of high quality willow seedlings and enhancement of land productivity.Key words: willow excellent; seedling; Logistic equation; growth model; gray correlation analysis为了解柳树良种各品种（系）苗木在甘肃省中部半干旱地区的生长规律，为扩大繁殖、培育良种壮苗提供理论依据，确定各品种（系）在苗木生产中最佳的经营措施，进行了苗期生长节律研究，利用Logistic曲线方程对一年生苗木的生长进行拟合，对影响苗木生长的空气温度、

7、降水量、空气相对湿度和日照时数4个因素进行了灰色关联度分析。1试验材料与方法1.1 试验材料参试的柳树良种品种（系）有4个，分别是从江苏林科院引进的J172柳（S.jiangsuensis CL.172）、J369柳（S.jiangsuensis CL.369）、J194柳(S.jiangsuensis CL.194)以及从新疆引进的青刚柳（S. viminalis L.）。1.2试验方法1.2.1试验地概况试验地设在甘肃省白银市大坪苗圃，土壤为轻壤土，肥力中等。海拔高度1521m，年平均降雨量260mm，年平均蒸发量1372.5mm，年日照时数2559.8h，年平均气温9，无霜期191d。利

8、用各无性系的一年生扦插苗截成粗细一致、长约15的插穗。采用随机区组设计，3次重复，以30株为小区，扦插株行距为2030，于4月2日扦插。1.2.2插后管理插后灌足生根水。大约13d 左右开始萌芽生长，当插条萌枝达2030cm时抹去多余侧枝,只保留1根健壮枝条。68月在苗木快速生长时期适时追肥3次，每次施尿素200250kg/hm，开沟埋施，8月底停止施肥。在苗木生长期发现食叶害虫危害时，喷施40%氧化乐果11000倍液3次，以控制虫害。1.2.3 观测方法扦插40天后开始进行生长节律调查，每小区固定10株苗木，每隔10d测量苗高和地径，直到9月25日苗木接近停止生长为止，共观测了14次。气象资

9、料由白银市气象局提供。1.2.4数学模型的建立利用小区平均值，采用 SAS 统计模块进行数据处理与分析1。许多研究表明，林木不论是个别器官或是整个植物体，不论是单株或群体，它们的重量、表面积、高度、细胞数量，甚至是蛋白质含量等生长曲线都近似于S型。虽然不同树种的生长规律不同，同一树种在不同经营水平以及不同生态条件下生长规律亦不尽一致，但生长曲线的总趋势是相似的2，杨柳科树木的苗期年生长节律符合S型曲线，表现为“慢一快一慢”的生长规律35。根据这一生长特点，用Logistic 曲线方程模拟柳树苗期苗高、地径的年生长动态。该方程的表达式为： （1）式（1）中y 为苗木累积生长量，K 表示苗木年生长

10、极限值，用苗高或地径的等差三点法求得，a 、b 是参数，用最小二乘法求得，t 为生长天数。 （2）对式（1）求二阶导数得令，则： （3）此时的为苗高、地径日生长速度最大值的日期,称之为速生点，通过表1中的月、旬与ti=i之间的对应关系换算得出。苗期苗高、地径年生长阶段划分，对式（1）求三阶导数得令，则： 解得，即 （4） （5）由式（4）、（5）可求出苗高、地径日生长速度变化最大的二点，即当时，为初期生长阶段,此阶段苗木缓慢生长；当时，为苗木中期生长阶段,此阶段生长迅速，称为速生期；当时，为苗木后期生长阶段，此阶段苗木缓慢生长。苗高、地径年生长规律拟合结果，经显著性检验以确定其相关指数。1.2

11、.5 苗木生长与各气象因素的灰色关联度分析灰色关联度分析的原理是根据因素数列的几何形状发展态势的接近程度来衡量其因素间关联度的大小，关联度大的数列与特征数列最为接近的原则6。各无性系的生长量与各气象因子组成一个灰色系统，其中净生长量为母序列X0，气温、降水量、相对湿度和日照时数为子序列Xi(i=14)。Xij1=Xij/Xij （6） 利用公式(6)对原始数值进行初值化变换，使其变成无量纲数值。 （7）利用公式（7）计算每个时刻点上母序列与各子序列差的绝对值。 （8）其中：为的最大值；为的最小值；为tj时刻值求在各时刻点上母序列X0与各子序列Xi的关联系数；K为分辨率系数，K=0.5。 （9）

12、利用公式（9）求关联度7。2 结果与分析2.1 生长模型的建立根据柳树良种品种（系）生长节律观测结果（见表1）可看出柳树苗期生长节律符合S型曲线，各品种（系）一年生扦插苗的生长进程均表现出明显的“慢快慢”的S型曲线。利用（2）式求出K值，利用最小二乘法求出系数a和b，建立各无性的生长模型。相关系数都达到0.90以上，相关显著。用F值进行显著性检验，都达极显著水平（PrF=0.0010.01）（表2），说明利用Logistic方程拟合一年生柳树扦插苗木的生长节律是可行的。表1：各品种（系）苗高、地径生长节律（cm）Table 1:Each species( department) seedlin

13、g is high, a path grows the stanza( cm)月旬TI=IJ172J194J369青刚柳苗高地径苗高地径苗高地径苗高地径5中1120.48110.47110.47100.46下2210.52190.50190.51170.496上中3464920.650.8056810.640.7753800.630.7750780.600.75下51200.921100.911090.901050.877上中671521811.051.711451731.011.131441711.001.131361620.951.06下82011.302001.251981.231871

14、.178上中9102202311.371.432162291.321.402142281.301.372092111.251.31下112511.462471.442431.432291.379上中12132822911.571.702742821.551.672732811.541.662542591.461.52下142961.742901.732871.732621.59表2：各柳树品种（系）苗高、地径的Logistic模型Table 2:Each willow species( department) seedling is high, the model of Logistic of

15、 a path各品种（系）苗 高地 径abK相关系数F值abK相关系数F值J1724.92470.0562321.30.9869724.363.56130.05372.20.9757752.63J1945.36360.0592311.60.99533249.143.34810.04252.60.9343298.22J3695.81390.0676305.60.99251643.313.21030.04922.30.98471270.15青刚柳6.01420.0728301.20.99392638.653.13520.04862.40.9052472.562.2苗木各生长时期的划分及生长特点利用

16、公式（3）、（4）和（5）求得t、t1和t2，t值是苗木连日生长量最大时的时期，可据此求出各无性系的速生点；依据t1和t2将扦插苗的年生长进程划分为3个时期：生长前期、速生期和生长后期，再加上苗木的成活期（从扦插到生根发芽的时期）共4个时期。表3：由回归方程导出的速生期及速生期内净生长量Table 3:From return to return the square distance to lead to living to expect soon and living the clean growth inside the period to measure soon各品种（系）苗 高地 径

17、速生点t速生期始期t1速生期结束期t2速生期占天数/d总生长期的比率/%速生期内生长量占总生长量比率/%速生点t速生期始期t1速生期结束期t2速生期占天数/d总生长期的比率/%速生期内生长量占总生长量比率/%J17298.372.0124.652.635.361.394.659.5129.870.347.267.6J19493.270.4116.245.830.760.388.358.2118.560.340.563.9J36991.470.0112.942.928.861.787.661.8119.657.838.864.5青刚柳88.369.0107.738.726.058.885.158

18、.0112.254.236.460.1平均92.870.4115.445.030.260.588.959.2121.260.740.764.0结果表明：柳树良种1a生苗高的生长峰值出现于萌芽后的92.8d（7月17日左右），与地径的平均速生点（t=88.9d）差异不大；地径的速生期始期平均在6月13日左右（t1= 59.2d），比苗高6月24日左右（tl = 70.4d）早11.2d；地径的速生期结束期8月14日左右（t2=121.2d），比苗高8月8日左右（t2 =115.4d）晚5.8d；地径的速生期持续时间（平均为60.7d）比苗高的速生期持续时间（平均为45.0d）长15.7d。不同品

19、种（系）的苗高的速生点和速生期始期差异不大；而速生期结束期以青刚柳最早，J172最晚；速生期持续时间以J172最长（52.6d），青刚柳(38.7d)最短，相差13.9d。不同品种（系）的地径的速生点差异较大，最早的品种（系）青刚柳(t=85.1d)与最晚的品种（系）J172（t=94.6d）相差9.5d；速生期始期和结束期差异也较大，皆以青刚柳最早，J172最晚；速生期持续时间以J172最长（70.3d），青刚柳最短(54.2d)，相差16.1d。虽然速生期占整个生长期的比率较小（苗高平均为30.2，地径平均为40.7），但速生期内生长量占总生长量的比率却较大（苗高平均为60.5，地径平均为

20、64.0），说明速生期内的生长量决定着整个生长期的生长量。由于速生期时间较短，生长量却较大，因此在苗木生产中，应重点抓好速生期内的合理追肥、除草、病虫害防治等育苗措施，以促进苗木的生长。在生长后期，应停止施肥，以促进苗木木质化，提高苗木质量，使苗木安全越冬，提高育苗、造林成活率。速生持续期的长短与生物量形成有密切关系，一般说来，速生期越长，生长量越大。在一个生长期内，无性系速生点出现的时间，以及速生期的持续时间，是衡量无性系能否充分适应环境的主要标准8。由表3可知，苗高速生期内的生长量占总生长量比率最高的是J369，达到61.7，其次是J172的为61.3，最低的是青刚柳；地径速生期内的生长量

21、占总生长量比率最高的是J172，达到67.6，其次是J369柳，最低的是青刚柳。由此可见，苏柳系的苗木生长占一定的优势。2.3 年生长量与各气象因子的灰色关联分析灰色关联度分析,其原理是根据因素数列的几何形状发展态势的接近程度来衡量其素间关联度的大小，关联度大的数列与特征数列最为接近的原则。根据各品种（系）苗木生长期气象因子（见表4），将每个品种（系）的苗高或地径与温度、降水量、相对湿度和日照时数看作一个灰色系统，温度、降水量、相对湿度和日照时数看作该系统的4个因素,以研究4个柳树品种（系）苗期年生长与温度、降水量、相对湿度和日照时数的关联程度。表4 各品种（系）苗木生长期气象因子table

22、4 each species(department) seedling wood growth period weather factors月旬ti=i空气温度/降水量/mm相对湿度/日照时数/h5中下1218222.6113036.231.0122.5132.46上中下34525.722.323.13.112.6039.842.543.3109.396.297.67上中下67824.923.725.66.767.327.650.651.057.999.252.496.38上中下9101125.323.124.372.115.1055.146.543.676.866.8105.49上中下121

23、31421.520.618.242.005.150.340.134.3113.0123.6130.7由式（8）、（9）计算各品种（系）的关联度表明（见表4）：各气象因子对苗高和地径生长节律的影响顺序在各品种（系）间有差别。但是总体来讲，气温对苗高影响最大，降水量影响最小；而日照时数对地径影响最大，降水量影响最小。表5 各无性系苗高、地径关联度及排序结果table 5 each asexualities of the seedling is high a path connection and the row prefaces are as a result项目各品种（系）气象因子关联度排序苗高

24、J172空气温度/降水量/mm相对湿度/%日照时数/h0.91210.74250.83320.78111423J194空气温度/降水量/mm相对湿度/%日照时数/h0.80900.75410.79030.77561423J369空气温度/降水量/mm相对湿度/%日照时数/h0.68900.67210.66750.63351234青刚柳空气温度/降水量/mm相对湿度/%日照时数/h0.67810.65520.66270.62561324地径J172空气温度/降水量/mm相对湿度/%日照时数/h0.89080.80810.94450.95063421J194空气温度/降水量/mm相对湿度/%日照时

25、数/h0.85060.80290.87210.90313421J369空气温度/降水量/mm相对湿度/%日照时数/h0.90220.84070.94270.97213421空气温度/降水量/mm相对湿度/%日照时数/h0.92760.84070.94230.972134213 结论与讨论3.1 利用 Logistic 方程拟合各柳树品种（系）苗高、地径的年生长节律，相关指数均在0.90 以上，拟合效果显著。3.2 对Logistic 曲线进行一次、二次求导所确定的速生点及速生期，在无性系选择及管理上具有重要意义10。各品种（系）的苗高、地径生长在生长前期、速生期和生长后期呈现出“慢快慢”的节律

26、，具有明显的阶段性。虽然速生期持续时间占整个生长期的比率较小，但速生期内生长量占总生长量的比率却较大，因此根据苗木的速生期、速生点来制定合理的水肥管理措施、促进苗木生长、提高苗木质量，从而降低生产成本，提高经济效益。3.3 苗高和地径的速生点差异不大，但地径的速生期始期较苗高早且持续时间较苗高长。不同品种（系）间苗高的速生点和速生期始期差异不大；而速生期结束期和速生期持续时间差异较大；地径的速生点、速生期始期、速生期结束期和速生期持续时间也差异较大，苏柳系的苗木生长占一定的优势。3.4在45月要重点做好抗旱工作，确保苗木成活，78月苗木高生长高峰期和89月地径生长高峰期，是管理的关键，应及时降

27、渍排水，追施肥料，创造苗木快速生长的适宜条件，确保齐苗壮苗，提高优质苗出圃率。3.5 各品种（系）的生长量受气象因子影响。气温对苗高影响最大，降水量影响最小；而日照时数对地径影响最大，降水量影响最小。3.6参试的柳树良种品种（系）苏柳新无性系均能适应本地土壤、气候等因子，生长迅速，适应性强，适宜本地大面积推广。参考文献1黄少伟，谢维辉.实用SAS编程与林业试验数据分析M.广州；华南理工大学出版社，20012罗世家，艾训儒，彭 诚.黄山松人工林材积生长规律的研究J.林业科技，2002，27（6）：463张俊生，刘建民，杨久延等.欧美杨类新无性系苗期生长模型及灰色关联度分析J.辽宁林业科技，199

28、7（5）：18-204崔长占，黄传林，于国民.等.小黑杨扦插苗生长规律的研究J.高师理科学刊，1999，19（2）：37395向志民，何 敏.几种杨树生长进程动态分析J.西北林学院学报，1994。9（2）：82866袁嘉祖.灰色系统理论及其应用M.北京：科学出版社，1991.7王学萌，张继忠，王荣.灰色系统分析及实用计算程序M.武昌：华中科技大学出版社，20018邬荣德，王明庥，黄敏仁.黑杨派新无性系研究J.南京林业大学学报，1998（4）：1-139Farance J（金之庆，高亮之译）.农业中的数学模型M.北京：农业出版社，199110徐宏远，陈章水.不同密度I-69杨树生长规律的研究J.林业科学研究，1994，7（1）：616611 秦光华，乔玉玲，孟昭和.美洲黑杨新无性T26和T66林期年高生长节律的研究J.江苏林业科技，2002，29（4）：68