湿地松活立木弹性模量无损评估技术研究

张帅楠，姜景民，徐永勤，栾启福*

(1.中国林业科学研究院亚热带林业研究所，浙江省林木育种技术研究重点实验室，浙江 杭州 311400；2.浙江省杭州市余杭区长乐林场，浙江 杭州 311123)

湿地松活立木弹性模量无损评估技术研究

张帅楠1，姜景民1，徐永勤2，栾启福1*

(1.中国林业科学研究院亚热带林业研究所，浙江省林木育种技术研究重点实验室，浙江 杭州 311400；2.浙江省杭州市余杭区长乐林场，浙江 杭州 311123)

[目的]开发操作简便、破坏性小且结果可靠的弹性模量无损评估技术，为实现以大群体试验材料为基础的遗传选择和评估，培育速生优质林木良种奠定基础。[方法]以22年生湿地松自由授粉家系测定林为研究对象，采用力学阻力仪(Pilodyn)和应力波速(V)记录仪(AV)等活立木无损评估仪器和方法以及木材物理力学常规检测方法，建立活立木弹性模量评估技术体系。[结果]Pilodyn测定值Pr与实测基本密度(ρ)间呈显著负相关(P<0.01)；弹性模量测定值MOEB与弹性模量(计算值)MOE(MOE=ρV2)之间呈极显著线性相关关系，其线性回归模型为MOEB=2.013 3+0.828 6MOE，决定系数R2=0.854 1；利用ρ和V计算的弹性模量值MOE与利用Pr和V计算的弹性模量相对值MOEP(MOEP=PrV2)间呈极显著正相关关系(P<0.001)，回归模型为MOE= 0.400 7+0.231 1MOEP, 决定系数R2=0.734 0；弹性模量MOE与V值相关极显著(P<0.001)，V与MOE的回归方程为MOE= -5.121 8+3.239 4V, 决定系数R2=0.921 9。[结论]在遗传选择时弹性模量(计算值)完全可以代表弹性模量测定值进行评估与选育；可直接利用MOEP作为各样木的弹性模量值参与遗传分析；同时建立了MOE与V值之间的回归方程，在比较同一测定群体弹性模量值相对大小时也可直接比较V值，大大减少了工作量。建立的模型对其它树种测定具有参考意义，但尚需检验。

湿地松；活立木；弹性模量；无损评估

为培育速生优质林木良种，需要大规模对遗传测定群体活立木材性等重要品质性状开展无损评估。其中，弹性模量(MOE)是林木建筑用材最重要的力学性质，能够反映木材的韧性和抵抗变形的能力，木材的弹性模量越大，在承受荷载时其变形越小[1-2]。目前，已成功应用于活立木弹性模量无损评估的是基于应力波速法的AV技术，它通过测量应力波在开始探头和结束探头之间的传播时间来测得AV值(波速V)，该数值与立木强度具有显著的相关性，可以作为活立木强度的检测工具[3-5]。Ross等[6]采用应力波检测方法得出木材的弹性模量(MOE)、木材密度(ρ)和应力波速度的关系式：MOE=ρV2，这为应力波速技术在木材弹性模量无损检测提供了研究基础[7-9]。活立木木材基本密度(ρ)测定的常规方法是用生长锥在其胸径处取木芯，然后利用饱和含水量法进行测定[10]，虽然该方法获得的数据直观、可靠，但是成本高、耗时长、具破坏性。Pilodyn是应用最为普遍的一种基于阻力的活立木密度无损快速检测仪器[11-12]，有效地克服了传统检测方法的缺点，在活立木木材密度评估中已有一定的进展[13-17]。在实现了应力波速(V)和基本密度(ρ)无损评估的基础上，Chen等[5, 18]综合利用了Pilodyn和AV技术，以Pilodyn探测值估测的木材外部密度作为木材基本密度，对挪威云杉(Piceaabies(L.) H. Karst.)大群体MOE和基本密度(ρ)进行了估测，首次实现了数千株大群体活立木材性性状的遗传变异分析。

湿地松(PinuselliottiiEngelman var.elliottii)原产于美国东南部，现已成为我国南方重要的用材树种[19-21]。经过数十年的育种工作，湿地松生长性状获得较大的改良[22-24]，以弹性模量为代表的力学性状由于受测定技术影响改良进程缓慢[25-28]。因此，开发操作简单方便、破坏性小且结果可靠的弹性模量无损评估技术，实现以大群体试验材料为基础的遗传选择和评估，对提高遗传测定的精度与效率、促进遗传改良进程具有重大现实意义。本研究采用力学阻力仪(Pilodyn)和应力波速记录仪(AV)等活立木无损评估仪器和方法以及木材物理力学常规检测方法，初步建立湿地松活立木弹性模量评估技术体系，为湿地松弹性模量等材性性状遗传改良奠定基础，同时对其它树种材性测定具有参考意义。

1 试验林概况

试验地位于浙江省杭州市余杭区长乐林场(30°27′ N，119°48′E)，亚热带湿润季风性气候，年均降水量1 399 mm，年均气温16.1℃，属低山缓坡，土壤以红壤为主。

试验林为湿地松子代测定林，营建于1994年春季，含33个半同胞家系，采取单行6株小区、6个重复的随机完全区组设计。

2 研究方法

2.1 材性性状指标测定方法

2.1.1 Pilodyn测定 2015年12月对上述33个湿地松半同胞家系进行取样测定，随机选取6个重复中的3个重复，其中，每个小区选取1株平均木作为样木，总计99个单株，进行Pilodyn(6J，PROCEQ，Switzerland)单株测定[5, 29]。在样木胸高南、北向分别测定2次(若2次探测值差值大于2 mm，则进行第3次探测，选择3次探测差异最小的2个值作为有效探测值)，取2次测定的平均值作为南向的Pilodyn值(Ps)和北向的Pilodyn值(Pn)，并取南、北向的平均值作为Pilodyn测定值(Pr)。

2.1.2 基本密度生长锥测定 Pilodyn测定后，用内径为5 mm的生长锥在样木胸高处南北向钻取木芯(木芯穿透样木)。基本密度(ρ)采用饱和含水量法测定[10]，其计算公式如下：

(1)

式(1)中:ρ为基本密度(g·cm-3);Gmw为饱和含水率时试样的质量(g);Gh为绝干时试样的质量(g)。

2.1.4 弹性模量评估

2.1.4.1 无损评估 根据生长锥法和AV技术得出的测定值ρ和V，应用公式MOE=ρV2计算出各样木的弹性模量(即弹性模量计算值，本文中以MOE表示)；根据Pilodyn法和AV技术得出的测定值Pr和V，以Pr代替ρ，应用上述公式计算出各样木的弹性模量相对值MOEP(MOEP=PrV2)。

2.1.4.2 常规测定 在99株样木中随机选取30株样木依据国家标准《GB/T 1929-2009 木材物理力学试材锯解及试样截取方法》进行破坏性取样，并委托浙江省林科院质检站按照国家标准《GB/T 1936.

2-2009 木材抗弯弹性模量测定方法》测定其抗弯弹性模量(即弹性模量测定值，本文中以MOEB表示)。

2.2 统计分析方法

数据分析在R软件(v 3.1.3，CRAN，http://cran.r-project.org)上进行。

(1)各指标间的相关性分析利用R软件中的corr. test()函数计算Pearson相关系数r并进行显著性检验。

(2)MOEB和MOE、MOE和MOEP以及MOE和V之间的回归分析均基于普通最小二乘法(OLS)，利用R软件中的lm()函数拟合其简单线性回归模型：

MOEB= b1+k1MOE

(2)

MOE = b2+k2MOEP

(3)

(4)

式(2)～(4)中:b1、b2、b3代表回归常数，k1、k2、k3代表回归系数。

3 结果与分析

3.1 各测定指标的基本情况

99个湿地松单株材性性状各测定指标的基本情况见表1。湿地松木材基本密度和弹性模量分别为0.502 0 g·cm-3、4.495 6 GPa，其中，MOEP极差与方差均较大，说明弹性模量值变异丰富，这可能是由于MOEP计算时因子Pr与V2的乘积引起的MOEP值波动较大，可能会有较好的遗传选择成效。基本密度性状测定指标(ρ、Pr)变异系数较小(CV<10%)，弹性模量性状测定指标(MOE、MOEP)变异系数较大(CV>20%)，也说明湿地松弹性模量性状相对于基本密度性状具有更丰富的遗传变异性。

表1 湿地松基本密度和弹性模量性状各测定指标的基本情况Table 1 Basic information of determination of wood density and MOE of slash pine

3.2 各测定指标间相关性分析

对湿地松木材基本密度性状测定指标(ρ、Pr)和弹性模量性状测定指标(V、MOE、MOEP)进行相关分析，结果见表2。从表2中可以看出:ρ与Pr呈显著负相关关系(P<0.01)，而这二者基本密度性状测定指标与3个弹性模量性状测定指标之间相关关系不显著或显著性低(P<0.1)；V与MOE及MOEP达到极显著的高度正相关(r>0.90，P<0.001)。这说明波速是影响湿地松弹性模量估测的关键因子，木材密度对弹性模量影响相对较小,而MOE与MOEP之间也呈极显著高度正相关(r=0.85，P<0.001)，因此，可以认为MOEP也可以对MOE进行预测和评估。

表2 湿地松基本密度和弹性模量性状各测定 指标间相关关系Table 2 The correlations of determination of wood density and MOE indexes of slash pine

注：括号中为该相关系数显著性水平；*代表显著性水平为0.1，**代表显著性水平为0.01，***代表显著性水平为0.001。

Note: Significance probability of the correlations were given in parenthesis; *, **, *** representP< 0.05,P< 0.01 andP< 0.001, respectively.

3.3 弹性模量测定值和弹性模量(计算值)的回归方程

湿地松弹性模量测定值与弹性模量(计算值)的线性回归模型为MOEB=2.013 3+0.828 6MOE，其决定系数R2=0.854 1(图1)。模型中，回归常数和回归系数都达极显著水平(P<0.001)，即MOEB与MOE呈极显著线性相关关系。因此，利用关系式MOE=ρV2得出的湿地松弹性模量计算值MOE虽然不能作为其标准值，却完全可以代表弹性模量测定值进行评估与选育。

图1 湿地松弹性模量测定值与弹性模量(计算值)的
简单线性回归模型
Fig.1 Simple linear regression model between MOEBand MOE

3.4 弹性模量(计算值)和弹性模量相对值的回归方程

湿地松弹性模量与弹性模量相对值的线性回归模型为MOE=0.400 7+0.231 1MOEP，其决定系数R2为0.734 0(图2)。模型中，回归常数(P<0.05)和回归系数(P<0.001)都达显著水平，这说明该线性模型适合于MOE与MOEP之间的回归分析，也就是说利用Pilodyn法和AV技术得出的测定值Pr和V，计算出MOEP并应用以上模型即能快速评估弹性模量。

图2 湿地松弹性模量值与弹性模量相对值的简单
线性回归模型
Fig.2 Simple linear regression model between MOE and MOEP

3.5 弹性模量(计算值)和波速的回归方程

湿地松弹性模量与波速的线性回归模型为MOE=-5.121 8+3.239 4V，其决定系数R2为0.921 9(图3)。模型中，回归常数和回归系数也都达极显著水平(P<0.001)，该模型的方差解释率为91.19%，这说明该线性模型能很好地拟合MOE与V之间的线性回归关系，即仅利用AV技术也能实现弹性模量评估。

图3 湿地松弹性模量值与波速V的简单线性回归模型
Fig.3 Simple linear regression model between MOE and V

4 讨论

弹性模量与树木本身抵抗倒伏能力密切相关[31]。湿地松作为当前我国南方最重要的用材、采脂及荒山绿化优良树种，种植面积超过200万hm2，然而，由于前期研究工作对其抗风雪等力学性质的选择力度不够，部分地区湿地松遇到异常天气情况出现雪压、风倒等灾害，成为制约该树种发展的一个重要方面[32]，因此，从树种健康角度考虑，需要继续加强其弹性模量等物理力学性质的评估与选育研究。随着无损评估技术的发展，全林选择测定湿地松弹性模量成为可能。本研究以力学阻力仪(Pilodyn)和应力波速记录仪(AV)等活立木无损评估仪器和方法，初步建立湿地松活立木弹性模量评估技术体系，为湿地松弹性模量等材性性状遗传改良奠定基础，同时对其它树种测定具有参考意义。

自Ross等[6]研究得出关系式MOE=ρV2以来，基于应力波速法的AV技术在雪松(Cedrusdeodara(Roxb.) G. Don.)[7]、欧洲赤松(PinussylvestrisL.)[8]、火炬松(PinustaedaL.)[33]、杉木(Cunninghamialanceolata(Lamb.) Hook.)[9]等树种活立木MOE评估中的研究应用逐渐成熟，为湿地松大群体活立木弹性模量的快速评估奠定了基础。本研究结果也显示出湿地松弹性模量测定值与弹性模量(计算值)之间存在着极显著的线性相关关系，其线性回归模型为MOEB=2.013 3+0.828 6MOE，决定系数R2=0.854 1。可以看出，在湿地松弹性模量评估中，关系式MOE=ρV2不能得出其标准弹性模量值，但这并不能否定该关系式的应用价值，在湿地松无损、快速评估及其遗传改良研究进程中，其仍具有相当大的优势，完全可以代表弹性模量测定值进行评估与选育。

湿地松木材基本密度性状测定指标和弹性模量性状测定指标相关分析结果显示：基本密度ρ与Pilodyn测定值Pr呈显著负相关关系。直接以Pr值代替ρ计算的弹性模量相对值MOEP与弹性模量MOE之间呈极显著的高度正相关关系，回归模型MOE=0.400 7+0.231 1MOEP，因此，在不要求精确获取各样木弹性模量绝对值时，可直接利用MOEP作为各样木的弹性模量值参与遗传分析。以Pr值代替ρ计算的弹性模量相对值MOEP较之先前弹性模量无损评估研究而言，节省了测定相应木材基本密度的工作量[8-9]，具有一定的便利优势，可在现实工作中应用。同时，相关性研究中还发现,基本密度性状测定指标(ρ、Pr)与弹性模量性状测定指标(MOE、MOEP)之间相关关系不显著或显著性较低，波速V与弹性模量MOE以及弹性模量相对值MOEP则均达到极显著高度正相关，说明波速V是弹性模量间接估测中的关键因子，而木材密度对其影响相对较小。弹性模量MOE与波速V之间的线性回归模型MOE= -5.121 8+3.239 4V，该模型很好地拟合了本研究中MOE与V之间的线性关系。可能受试验材料影响，Wu等[34]以22个4年生桉树(Eucalyptusspp.)无性系为研究对象，结果认为Pilodyn与弹性模量呈显著负相关，而本研究中结果显示，尽管这二者之间存在负相关关系，但相关是不显著的，与Chen等[5]对5 618株挪威云杉材性无损测定的结果相同。因此，木材基本密度与弹性模量之间的相关性受试验对象或环境条件影响较大，不同研究对象评估时需要建立对应的评估模型。

一般认为弹性模量由物体基本密度(ρ)和应力波速(V)共同决定，由经验方程MOE=ρV2可以看出V的影响应该更大，本研究中仅应用V就能很好地预测弹性模量，可以将公式中ρ看作是一个常数，这和Kennedy等[9]的观点一致，为活立木弹性模量的简单快速预测奠定了基础，同时研究中直接用Pilodyn测定值Pr代替ρ计算的弹性模量估计值MOEP也可以作为活立木弹性模量的参考值。

5 结论

(1)湿地松弹性模量测定值MOEB与弹性模量(计算值)MOE之间存在极显著线性相关关系，在遗传选择时MOE完全可以代表弹性模量测定值进行评估与选育。

(2)弹性模量相对值MOEP与弹性模量MOE之间呈极显著高度正相关关系，回归模型MOE=0.400 7+0.231 1MOEP，可直接利用MOEP作为各样木的弹性模量值参与遗传分析。

(3)弹性模量MOE与波速V之间的线性回归模型MOE=-5.121 8+3.239 4V，在比较同一测定群体弹性模量值相对大小时也可直接比较V值，大大减少了工作量。

以上建立的模型对其它树种测定具有参考意义，但尚需检验。

[1] Raymond C A, Kube P D, Pinkard L,etal. Evaluation of non-destructive methods of measuring growth stress inEucalyptusglobulus: relationships between strain, wood properties and stress[J]. Forest Ecology and Management, 2004, 190(2): 187-200.

[2] Lasserre J P, Mason E G, Watt M S. Assessing corewood acoustic velocity and modulus of elasticity with two impact based instruments in 11-year-old trees from a clonal-spacing experiment ofPinusradiataD. Don.[J]. Forest Ecology and Management, 2007, 239(1): 217-221.

[3] Grabianowski M, Manley B, Walker J C F. Acoustic measurements on standing trees, logs and green lumber[J]. Wood Sci Technol, 2006, 40(3): 205-216.

[4] Wessels C B, Malan F S, Rypstra T. A review of measurement methods used on standing trees for the prediction of some mechanical properties of timber[J]. European Journal of Forest Research, 2011, 130(6): 881-893.

[5] Chen Z Q, Karlsson B, Lundqvist S O,etal. Estimating solid wood properties using Pilodyn and acoustic velocity on standing trees of Norway spruce[J]. Annals of Forest Science, 2015, 72(4): 499-508.

[6] Ross R J, Pellerin R F. NDE of wood-based composites with longitudinal stress wave[J]. Forest Products Journal, 1988, 38(5):39-45.

[7] Nanami N, Nakamura N, Arima T,etal. Measuring the properties of standing trees with stress waves, 3: Evaluating the properties of standing trees for some forest stands[J]. Journal of the Japan Wood Research Society, 1993, 39(8):903-909.

[8] Auty D, Achim A. The relationship between standing tree acoustic assessment and timber quality in Scots pine and the practical implications for assessing timber quality from naturally regenerated stands[J]. Forestry, 2008, 81(4): 475-487.

[9] Kennedy S G, Cameron A D, Lee S J. Genetic relationships between wood quality traits and diameter growth of juvenile core wood in Sitka spruce[J]. Can J For Res, 2013, 43(1): 1-6.

[10] 成俊卿. 木材学[M]. 北京: 中国林业出版社，1985: 282-464.

[11] Yanchuk A D, and Kiss G K. Genetic variation in growth and wood specific gravity and its utility in the improvement of interior spruce in British Columbia[J]. Silvae Genet, 1993,42(2-3): 141-148.

[12] Koch L, Fins L. Genetic variation in wood specific gravity from progeny tests of ponderosa pine (PinusponderosaLaws.) in northern Idaho and Western Montana[J]. Silvae Genet, 2000, 49(4-5): 174-181.

[13] Cown D J. Comparison of the Pilodyn and torsiometer methods for the rapid assessment of wood density in living trees[J]. N Z J For Sci, 1978, 8(3): 384-391.

[14] Taylor F W. Rapid determination of southern pine specific gravity with a Pilodyn tester[J]. For Sci,1981, 27(1): 59-61.

[15] 朱景乐，王军辉，张守攻，等. 毛白杨材性指标预测及选择[J]. 林业科学，2008，44(7): 23-28.

[16] 栾启福，卢 萍，井振华，等. Pilodyn评估杂交松活立木的基本密度及其性状相关分析[J]. 江西农业大学学报，2011，33(3): 548-552.

[17] 易 敏，赖 猛，孙晓梅，等. 无损检测在日本落叶松材性育种中的应用前景探讨[J]. 林业科学，2014，50(11): 96-103.

[18] Chen Z Q, Gil M R G, Karlsson B,etal. Inheritance of growth and solid wood quality traits in a large Norway spruce population tested at two locations in southern Sweden[J]. Tree Genetics & Genomes, 2014, 10(5): 1291-1303.

[19] 曲桂林，吴志民. 国外林业技术考察报告选编[M]. 北京: 中国林业出版社，1996: 63-72.

[20] Wen X, Kuang Y, Shi M,etal. Biology of Hylobitelus xiaoi (Coleoptera: Curculionidae), a new pest of slash pine, Pinus elliottii[J]. Journal of Economic Entomology, 2004, 97(6): 1958-1964.

[21] Barnett J P, Sheffield R M. Slash pine: characteristics, history, status and trends[C]// Dickens E D, Barnett J P, Hubbard W G,etal. Slash pine: still growing and growing! Proceedings of the slash pine symposium, Jekyll Island, Georgia, USA, 23-25 April 2002. Gen. Tech. Rep. SRS-76. Asheville, NC: U.S. Department of Agriculture, Forest Service, Southern Research Station, 2005:1-6.

[22] 潘志刚，郑勇奇. 湿地松种源试验研究[J]. 林业科学研究，1990，3(4): 314-322.

[23] White T L, Hodge G R, Powell G L. An advanced-generation tree improvement plan for slash pine in the southeastern United States[J]. Silvae Genetica, 1993(42):359-359.

[24] 姜景民，孙海菁，刘昭息. 湿地松自由授粉家系的早期评定[J]. 林业科学研究，1995，8(5):574-581.

[25] Zobel, B J. Inheritance of wood properties in conifers[J]. Silvae Genet, 1961,10(3): 65-96.

[26] Wei X, Borralho N M G. Genetic control of basic density and bark thickness and their relationships with growth traits ofEucalyptusurophyllain south east China[J]. Silvae Genetica, 1997, 46(4):245-250.

[27] Hansen C P. Application of the Pilodyn in forest tree improvement. DFSC Series of Technical Notes. TN55[R]//Humlebaek, Denmark: DANIDA Forest Seed Centre, 2000.

[28] 姜笑梅，骆秀琴，殷亚方，等. 不同湿地松种源木材材性遗传变异的研究[J]. 林业科学，2002，38(3):130-135.

[29] 朱景乐，王军辉，张守攻，等. Pilodyn在林木遗传改良应用中的研究进展[J]. 浙江林学院学报，2008，25(5): 661-665.

[30] Carter P, Briggs D, Ross R J,etal. Acoustic testing to enhance western forest values and meet customer wood quality needs[R] //Harrington C A, Schoenholtz S H. Productivity of western forests: a forest products focus. General technical report PNW-GTR-642. Portland: U.S. Department of Agriculture, Forest Service, Pacific Northwest Research, 2005: 121-129.

[31] 许秀玉，王明怀，仲崇禄，等. 不同树种木材性质及其抗台风性能[J]. 浙江农林大学学报，2014，31(5): 751-757.

[32] 栾启福，卢 萍，肖复明，等. 雨雪冰冻重灾区湿地松受害情况调查及其原因初步分析[J]. 林业科学，2008，44(11): 50-54.

[33] Mora C R, Schimleck L R, Isik F,etal. Relationships between acoustic variables and different measures of stiffness in standing Pinus taeda trees[J]. Canadian Journal of Forest Research, 2009, 39(8): 1421-1429.

[34] Wu S, Xu J, Li G,etal. Use of the Pilodyn for assessing wood properties in standing trees ofEucalyptusclones[J]. Journal of Forestry Research, 2010, 21(1): 68-72.

(责任编辑：金立新)

Study on the Modulus of Elasticity Non-destructive Evaluation Technique of Slash Pine Standing Tree

ZHANGShuai-nan1,JIANGJing-min1,XUYong-qin2,LUANQi-fu1

(1.Research Institute of Subtropical Forestry, Chinese Academy of Forestry, Key Laboratory of Tree Breeding of Zhejiang Province,Hangzhou 311400, Zhejiang, China; 2.Changle Forestry Farm of Yuhang District, Hangzhou, Hangzhou 311123, Zhejiang, China)

[Objective]To develop the non-destructive evaluation technique of modulus of elasticity (MOE) for realizing the genetic selection and evaluation of forest tree with the traits of fast growing and high quality based on a large population of test materials. [Method]The evaluation systems of MOE of standing tree were established by non-destructive evaluation technique, Pilodyn and acoustic velocity (AV), and conventional testing methods based on the open pollinated family of 22-year-old slash pine (Pinuselliottii).[Result]The Pilodyn value (Pr) and the basic density (ρ) showed a significant negative correlation (P<0.01). The measurement value ofMOE(MOEB), determined by conventional methods, had a very significant linear correlation with theMOEvalue calculated by the formulaMOE=ρV2. The regression model wasMOEB=2.013 3+0.828 6MOE(Coefficient of determinationR2=0.854 1). A significant positive correlation (P<0.001) was observed betweenMOEandMOEP, the relative value ofMOE, (MOEP=PrV2) (Vmeans the value ofAV), and the linear regression model wasMOE= 0.400 7+0.231 1MOEP(coefficient of determinationR2=0.734 0). The positive correlation betweenMOEandVwas extremely significant (P<0.001) and the linear regression model betweenMOEandVwasMOE= -5.121 8+3.239 4V(coefficient of determinationR2=0.921 9).[Conclusion]The modulus of elasticity calculated byMOE=ρV2can represent the determination values ofMOE. TheMOEPcan be used directly as theMOEvalue for each sample tree to participate in genetic analysis. It is feasible to compare only theVwhen evaluating the relative value ofMOEof the sample trees in the same population. Therefore, these linear models are also useful for other tree species.

slash pine; standing tree; modulus of elasticity; non-destructive evaluation

10.13275/j.cnki.lykxyj.2017.01.011

2016-05-12

国家自然科学基金项目(31570668)；“十二·五”国家林业科技支撑计划专题(2012BAD01B0203)；国家林业局“948”引进项目(2012-4-46)

张帅楠(1993—)，男，河南安阳人，硕士研究生，主要从事林木遗传育种研究. * 通讯作者：栾启福，副研究员，博士，从事林木遗传育种研究. E-mail: qifu.luan@caf.ac.cn

S791.246

A

1001-1498(2017)01-0075-06