林业科学  2014, Vol. 50 Issue (11): 96-103   PDF    
DOI: 10.11707/j.1001-7488.20141113
0

文章信息

易敏, 赖猛, 孙晓梅, 张守攻, 赵鲲
Yi Min, Lai Meng, Sun Xiaomei, Zhang Shougong, Zhao Kun
无损检测在日本落叶松材性育种中的应用前景探讨
Application Prospect on Applying Nondestructive Testing Technology in Wood Property Breeding of Larix kaempferi
林业科学, 2014, 50(11): 96-103
Scientia Silvae Sinicae, 2014, 50(11): 96-103.
DOI: 10.11707/j.1001-7488.20141113

文章历史

收稿日期：2013-12-10
修回日期：2014-01-08

作者相关文章

易敏
赖猛
孙晓梅
张守攻
赵鲲

引用本文
易敏, 赖猛, 孙晓梅, 张守攻, 赵鲲. 2014. 无损检测在日本落叶松材性育种中的应用前景探讨[J]. 林业科学, 50(11): 96-103.
Yi Min, Lai Meng, Sun Xiaomei, Zhang Shougong, Zhao Kun. 2014. Application Prospect on Applying Nondestructive Testing Technology in Wood Property Breeding of Larix kaempferi. Scientia Silvae Sinicae, 50(11): 96-103.  DOI: 10.11707/j.1001-7488.20141113
无损检测在日本落叶松材性育种中的应用前景探讨
易敏1, 赖猛1, 孙晓梅1 , 张守攻1, 赵鲲2    
1. 中国林业科学研究院林业研究所 林木遗传育种国家重点实验室 北京 100091；
2. 河南洛阳林业科学研究所 洛阳 471002
收稿日期：2013-12-10；修回日期：2014-01-08
基金项目：国家“十二五” 科技支撑项目“北方针叶树种高世代育种技术研究与示范”(2012BAD01B01)。
通讯作者：孙晓梅
摘要：使用Pilodyn和微秒计对20个日本落叶松无性系进行无损检测，研究无性系间无损检测指标的表型和遗传变异，结果表明： 应力波速(SWV)的表型变异系数和遗传变异系数分别为5.09%和3.57%，低于Pilodyn南向测定值(Ps)和北向测定值(Pn)； 3个无损检测指标无性系重复力的分布范围为0.34~0.48，说明它们受中等强度的遗传控制。研究无损检测指标与各材质性状间的关系，结果表明： Ps和Pn与生长轮密度、早材密度、生长轮壁腔比和早材壁腔比呈极显著负相关； SWV与生长轮微纤丝角、早材微纤丝角和晚材微纤丝角呈极显著负相关，但与生长轮弹性模量、早材弹性模量和晚材弹性模量呈极显著正相关。多元回归分析表明，对Ps，Pn和SWV贡献率排名前2位的材质性状分别为生长轮密度和早材壁腔比、生长轮密度和生长轮壁腔比以及早材微纤丝角和生长轮弹性模量。建立无损检测指标与相对应材质性状的线性回归方程，发现Ps与生长轮密度间方程的决定系数(0.340)略小于Pn与生长轮密度的决定系数(0.360)； SWV与早材微纤丝角、生长轮弹性模量间方程的决定系数分别为0.348和0.316。
关键词： 日本落叶松    无损检测    材质性状    相关关系    线性回归    
Application Prospect on Applying Nondestructive Testing Technology in Wood Property Breeding of Larix kaempferi
Yi Min1, Lai Meng1, Sun Xiaomei1, Zhang Shougong1, Zhao Kun2    
1. State Key Laboratory of Tree Genetics and Breeding Research Institute of Forestry, CAF Beijing 100091;
2. Forestry Research Institute of Luoyang, Henan Province Luoyang 471002
Abstract: The phenotypic and genetic variation of nondestructive testing traits among 20 clones of the living trees of Larix kaempferi by using Pilodyn and microsecond timer were studied. Variant analysis revealed that the phenotypic and genetic variation coefficients for SWV(stress wave velocity) were 5.09% and 3.57%, which higher than those of Ps and Pn. The clonal repeatability for the nondestructive testing traits ranged from 0.34 to 0.48 in this study, which meant that variation in nondestructive testing traits of L. kaempferi were controlled genetically at medium level. The correlations between nondestructive testing traits and wood properties were analysed, and the results showed that there were significant negative correlations between Ps, Pn and ring density, early wood density, ring wall thickness to lumen area and early wood wall thickness to lumen area. The SWV had significant negative correlations with ring microfibrillar angle, early wood microfibrillar angle and late wood microfibrillar angle, but had significant positive correlations with ring modulus of elasticity, early wood modulus of elasticity and late wood modulus of elasticity. The results of the multiple regression analysis showed that the top two contribution rates of wood property to the Ps, Pn and SWV were ring density and early wood wall thickness to lumen area, ring density and ring wall thickness to lumen area, and early wood microfibrillar angle and ring modulus of elasticity. Linear regression equations between nondestructive testing traits and corresponding wood properties were constituted, and the results showed that the coefficient of determination of linear regression equation between Ps and ring density (0.340) was slightly smaller than those of Pn (0.360). The coefficients of determination of linear regression equations between SWV and early wood microfibrillar angle, ring modulus of elasticity were 0.348 and 0.316.
Key words: Larix kaempferi    nondestructive testing    wood property    correlation coefficient    linear regression    

日本落叶松(Larix kaempferi)是我国北方及南方亚高山地区的重要针叶造林树种，具有早期速生、成林快、易于栽培、适应性广等特点，是速生针叶材纸浆树种，受到国内外广泛重视(马常耕，1992)。随着落叶松木材需求量的增长，木材材性改良已成为落叶松遗传改良的重要目标。木材密度是影响木材材质、纸浆产量以及木材强度的重要物理量(Pliura et al.，2006； Schimleck et al.，1999)； 纤维壁腔比是评定木材纤维形态和纸张质量的重要指标(孙晓梅等，2003)； 木材硬度，又称弹性模量(MOE)，是林木作为建筑用材最重要的力学性质，木材的弹性模量越大，在承受荷载时其变形越小； 微纤丝角为细胞次生壁S2层微纤丝与细胞主轴所形成的夹角，其对木材主要物理、化学和力学性能都有着很大的影响，是重要的木材微观材性指标(潘彪等，2005； Cave et al.，1994； Walker et al.，1995； 1999)。

材质性状常规检测方法耗时耗力、经济成本高、精度低、对林木的损害大，如木材密度、纤维壁腔比、微纤丝角等材质性状的测定均需要木芯样品进行测量。活立木木芯取样不仅在一定程度上损坏树木，而且对于大批量样品材质性状的测定程序繁杂、耗时较长。然而，无损检测技术有效克服了传统检测方法的缺点，其操作简单、破坏性小，目前在各领域均有广泛应用。Pilodyn和微秒计均为便携式无损检测仪器。大量研究表明，Pilodyn与木材密度存在极显著的负相关关系，可用于木材密度的早期选择，加速林木育种进程(朱景乐等，2009a；2009b；Greaves et al.，1996)；微秒计可用于检测应力波在树木中传播速度的大小。有研究表明，在辐射松(Pinus radiata)(Chauha et al.，2006)、云杉(Picea asperata)(Wang et al.，2000)、日本落叶松(Ishiguri et al.，2008)和邓恩桉(Eucalyptus dunnii)(Dickson et al.，2003)等树种上，应力波速与弹性模量存在较高的正相关关系，利用应力波速可较好地预测弹性模量，但影响Pilodyn测定值和应力波速的不仅仅是木材密度和弹性模量，纤维壁腔比、微纤丝角等材质性状与无损检测指标亦有一定的关系。本文通过对日本落叶松木材密度、纤维壁腔比、微纤丝角和弹性模量等材质性状的测定研究，旨在探求Pilodyn测定值、应力波速与各材质性状间的关系，评估Pilodyn和微秒计在预测日本落叶松材质性状中的应用价值，从而为无损检测在日本落叶松材性育种中的应用提供理论依据。

1 材料与方法 1.1 试验林概况

试验林位于河南省嵩县五马寺林场，海拔1 400～1 600 m，年均气温8.4 ℃，最高温度24.7 ℃，最低温度-15.5 ℃，无霜期178天，年降水量800～1 200 mm，相对湿度80％。林地土壤为棕壤，pH为6.0。

试验林营建于1998年春，试验材料包括83个日本落叶松无性系。造林采用带状整地，完全随机区组设计，4次重复，4株小区，双行排列，株行距均为2 m×2 m。

1.2 Pilodyn和微秒计的工作原理

本研究选用瑞士PROCEQ公司产、6-J能量、2.5 mm直径探针型号的Pilodyn。Pilodyn的测定原理是以预先设定好的能量，将一个固定规格的探针射入到木材中，探针射入深度就是检测的结果。这与木材密度有紧密负相关关系，木材密度大，则射入深度浅; 反之，则射入深度深(Wei et al.，1997； Greaves et al.，1996)。

微秒计通过测量应力波在开始探头、结束探头、传感器和树干的时间来测得应力波速(Knowles et al.，2004； Chauhan et al.，2006)。应力波速(m·s-1)为应力波传播的距离除以传播时间，公式如下：

$ SWV = L/t。 $

式中： L =1 m，为开始探头与结束探头的距离； t 为传播时间(μs)。

1.3 无损检测指标测定

在试验林内，从83个无性系中随机选取20个无性系，每个无性系选取3个单株(保证每个单株不在同一区组)，共60株样木。Pilodyn的测定方法是在胸高处分正南方向剥去活立木少量外部树皮，用Pilodyn探测2次。如果2次探测值的误差在2 mm以内，则取2次探测值的平均值作为南向的计算数据； 如果2次探测值的误差大于2 mm，则需要进行第3次南向探测，然后在3个数据中选择2个误差在2 mm以内且数据最接近的2个探测值的平均数作为Ps。用同样的方法测得北向的Pilodyn测定值Pn。

采用六通道Fakopp(FRS-06/00，匈牙利)进行应力波速测定，将2个传感器探针以45°角沿纵向分别钉入样木1.5 m和0.5 m处，深度约为5 mm，随后用专用锤敲击其中一个传感器，记录传播时间，对每株样木测量3次，取其平均值。

1.4 材性取样及测定方法

无损检测完成后，在胸高1.3 m处用内径为12 mm的生长锥由南面从韧皮部直至髓心取一无疪木芯，在室内利用澳大利亚联邦科学与工业研究院(CSIRO)自主研发的木材性质快速测定仪(SilviScan-3)进行材性指标的测定。测定的材性指标包括木材密度、纤维壁腔比、微纤丝角和弹性模量。为得到每项指标的生长轮均值、早材均值和晚材均值，需在测定完成后进行定年和早晚材划分处理(Evans et al.，2001)。

1.5 数据统计分析方法

对无损检测数据进行有关遗传参数的估算，采用的线性固定模型(Hansen et al.，1996)如下：

$ {y_{ij}} = \mu + {\alpha _i} + {\beta _i} + {\varepsilon _{ij}}。 $

式中： yij代表第j区组内第i无性系的观测值； μ代表平均值； αi代表第i无性系的固定效应； βj代表第j区组的固定效应； εij代表随机误差。

无性系重复力计算公式(Hansen et al.，1996)如下：

$ R = \sigma _c^2/\left( {\sigma _c^2 + \sigma _e^2/r} \right)。 $

式中： r 代表区组数； σc2代表无性系方差； σe2代表随机误差。

遗传变异系数公式(Hai et al.，2008)如下：

$ CVG = 100 \times \sqrt {\sigma _c^2} /\bar X。 $

式中： $ {\bar X} $代表表型均值。

性状间表型相关关系如下：

$ {r_p} = {\sigma _{p\left( {xy} \right)}}/\sqrt {\sigma _{p\left( x \right)}^2 \times \sigma _{p\left( y \right)}^2} 。 $

式中： σp(x)2代表x性状表型方差分量； σp(y)2代表y性状表型方差分量； σp(xy)代表性状间表型协方差。

性状间遗传相关关系(Pliura et al.，2007)如下：

$ {r_g} = {\sigma _{c\left( {xy} \right)}}/\sqrt {\sigma _{c\left( x \right)}^2 \times \sigma _{c\left( y \right)}^2} . $

式中： σ2c(x)代表x性状无性系方差分量； σ2c(y)代表y性状无性系方差分量； σc(xy)代表性状间无性系协方差。

方差分量计算用SAS/STAT 9.0软件PROC VARCOMP过程中带限制的最大似然法计算(黄少伟等，2001)。

建立多元回归模型，由各变量标准化偏回归系数的绝对值占标准化偏回归系数绝对值总和的百分比计算各自贡献率，研究各材质性状对Pilodyn南、北向测定值、应力波速贡献的大小(朱景乐等，2009a)，模型为:

$ Y = {B_0} + {B_1}{X_1} + {B_2}{X_2} + \cdots \cdots + {B_n}{X_n}。 $

采用回归分析方法建立Pilodyn南、北向测定值、应力波速与材质性状的回归方程(王军辉等，2008； 张瑛春等，2010)，模型为：

$ y = ax + b。 $
2 结果与分析 2.1 无损检测指标的变异分析 2.1.1 无损检测指标的表型变异分析

日本落叶松无性系无损检测指标表型变异如表 1所示。由表 1可知，日本落叶松无性系Pilodyn南向测定值(Ps)最低为10.25 mm，最高为17.42 mm，Pilodyn北向测定值(Pn)最低为10.95 mm，最高为17.55 mm； Ps和Pn的无性系均值差异不大，分别为13.93和13.53 mm，Ps变异系数为13.53%，略高于Pn(11.88%)。日本落叶松无性系应力波速(SWV)最低为3 344.48 m·s-1，最高为3 989.36 m·s-1，均值为3 680.56 m·s-1； 应力波速表型变异系数仅为5.09%，显著低于Pilodyn南、北方向探测值的表型变异系数(Ps=13.53%，Pn=11.88%)。

表 1 日本落叶松无性系无损检测指标平均值、变幅及表型变异系数① Tab.1 Mean values, ranges and phenotypic variation coefficients of nondestructive testing traits for the clones of Larix kaempferi
2.1.2 无损检测指标的遗传变异分析

表 2列出了日本落叶松无性系无损检测指标的方差分量、无性系重复力、遗传方差占总方差比率、遗传变异系数等遗传参数估算结果。从表中可看出，3个无损检测指标无性系重复力的分布范围为0.34～0.48，说明这些性状受中等强度的遗传控制，其中SWV重复力要高于Ps和Pn； 遗传方差占总方差比率的分析结果与重复力相同，SWV数值最大为0.235，高于Ps(0.181)和Pn(0.147)。3个无损检测指标中，应力波速遗传变异系数最小仅为3.57%，低于Pilodyn南、北向测定值(Ps=8.54%，Pn=6.92%)。

表 2 日本落叶松无性系无损检测指标遗传参数估计 Tab.2 The estimation of genetic parameters of nondestructive testing traits for the clones of Larix kaempferi
2.2 无损检测指标与材质性状的相关关系

表 3列出了日本落叶松无性系无损检测指标与各材质性状的表型和遗传相关系数。从表中可看出，Ps和Pn与生长轮密度、早材密度、生长轮壁腔比和早材壁腔比之间的表型相关均达极显著水平，表现出极强的负向相关关系(-0.609～-0.406)； Ps和Pn与这4个性状的遗传相关均达-0.800以上。SWV与生长轮微纤丝角、早材微纤丝角、晚材微纤丝角、生长轮弹性模量、早材弹性模量和晚材弹性模量之间的表型相关均达极显著水平，其中与生长轮微纤丝角、早材微纤丝角、晚材微纤丝角为负相关关系(-0.633～-0.573)，与生长轮弹性模量、早材弹性模量和晚材弹性模量为正相关关系(0.475～0.573)； 除晚材弹性模量外，SWV与其余5个材质性状的遗传相关均达0.700以上。

表 3 无损检测指标与各材质性状间的相关关系① Tab.3 The correlations between nondestructive testing traits and wood properties
2.3 材质性状对无损检测指标贡献率剖析

由前面结果可知，生长轮密度、早材密度、生长轮壁腔比和早材壁腔比与Ps和Pn紧密相关； Ps和Pn与这4个性状之间的关系可用多元回归模型表示： Y=B0 + B1X1 + B2X2 +B3X3 + B4X4。SWV与生长轮微纤丝角、早材微纤丝角、晚材微纤丝角、生长轮弹性模量、早材弹性模量和晚材弹性模量相关紧密，其与这6个性状之间的关系可用多元回归模型表示： Y=B0 + B1X1 + B2X2 +B3X3 + B4X4 + B5X5 + B6X6。

由表 4可知，对Ps贡献率最大的是生长轮密度，达47.79%； 其次为早材壁腔比，为27.85%；生长轮壁腔比排名最后，为9.51%。Pn贡献率的排序与南向相比稍有不同，排名第一的依然为生长轮密度，生长轮壁腔比排名第二，为33.79%； 早材密度(8.41%)和早材壁腔比(7.86%)贡献率相差不大。

表 4 Pilodyn南、北向测定值与材质性状的多元回归分析 Tab.4 The multiple regression analysis for Pilodyn south penetration, Pilodyn north penetration and wood properties

由表 5可知，早材微纤丝角对SWV贡献率最大，达28.09%，其次为生长轮弹性模量，为25.94%；生长轮微纤丝角贡献率最低，仅为1.75%。

表 5 应力波速与材质性状的多元回归分析 Tab.5 The multiple regression analysis for stress wave velocity and wood properties
2.4 无损检测指标和材质性状线性回归模型的建立

对Ps，Pn和SWV贡献率排名前2位的材质性状分别为生长轮密度和早材壁腔比、生长轮密度和生长轮壁腔比以及早材微纤丝角和生长轮弹性模量。建立相应的线性回归方程，各参数值和方程的决定系数如表 6所示。由表 6可知，Ps，Pn和SWV与相对应材质性状间均存在极显著的回归关系(P＜0.000 1)。

表 6 无损检测指标与相关材质性状的线性回归分析 Tab.6 The linear regression analysis for nondestructive testing traits and corresponding wood properties

使用Ps和Pn预测木材生长轮密度的线性方程如图 1和图 3所示。由图 1和图 3可知，Ps与生长轮密度间方程的决定系数为0.340，略小于Pn与生长轮密度的决定系数(0.360)； Ps与早材壁腔比、Pn与生长轮壁腔比间方程的决定系数分别为0.228和0.245(图 2、图 4)。SWV与早材微纤丝角、生长轮弹性模量间关系的线性方程如图 5和图 6所示，由图可知，SWV与早材微纤丝角、生长轮弹性模量间方程的决定系数分别为0.348和0.316。

图 1 Pilodyn南向测定值(Ps)与生长轮密度(RD)的线性回归方程 Fig. 1 The linear regression equation between Ps and RD
图 2 Pilodyn南向测定值(Ps)与早材壁腔比(ET/D) 的线性回归方程 Fig. 2 The linear regression equation between Ps and ET/D
图 3 Pilodyn北向测定值(Pn)与生长轮密度(RD)的线性回归方程 Fig. 3 The linear regression equation between Pn and RD
图 4 Pilodyn北向测定值(Pn)与生长轮壁腔比(RT/D)的线性回归方程 Fig. 4 linear regression equation between Pn and RT/D
图 5 应力波速(SWV)与早材微纤丝角(EMFA) 的线性回归方程 Fig. 5 The linear regression equation between SWV and EMFA
图 6 应力波速(SWV)与生长轮弹性模量(RMOE)的线性回归方程 Fig. 6 The linear regression equation between SWV and RMO
3 结论与讨论 3.1 日本落叶松无损检测指标变异分析

Pilodyn南、北向无性系均值分别为13.93和13.53 mm，此结果要小于已发表的日本落叶松(朱景乐等，2009a；王军辉等，2008)和其他树种(Wu et al.，2010；朱景乐等，2009b)的数值。张瑛春等(2010)对河南嵩县79株日本落叶松Pilodyn测定值进行分析后发现，Pilodyn南向均值为15.3 mm，北向均值为14.6 mm。Kien等(2010)对建立在越南北部、中北部和南部的172个赤桉(Eucalyptus camaldulensis)子代进行分析，结果表明，Pilodyn平均值在14.6～15.51 mm之间。本研究中，应力波速无性系均值为3 680 m·s-1，该结果要高于Ishiguri等(2008)对日本落叶松的研究结果(3 180 m·s-1)。应力波在不同树种间的传播速度差异较大，如Matheson等(2008)对2片辐射松子代测定林进行研究，结果表明，2片试验林应力波速均值分别为1 927和2 160 m·s-1；Yin等(2010)利用3种无损检测方法评估杉木(Cunninghamia lanceolata)的材质性状，发现杉木应力波速均值为3 940 m·s-1；Wang等(2001)对168棵加州铁杉(Tsuga heterophylla)和西加云杉(Picea sitchensis)的研究表明，其应力波速均值分别为3 373和3 528 m·s-1。Pilodyn南、北向测定值表型变异系数差异不大(Ps =13.53%，Pn=11.88%)，该结果与张瑛春等(2010)对日本落叶松的研究结论一致。应力波速表型变异系数仅为5.09%，该结果要低于Chauhan等(2006)对辐射松的研究结果。

应力波速重复力(0.48)要高于Pilodyn南、北向测定值(Ps =0.40，Pn=0.34)，遗传方差占总方差比率的分析结果与重复力相同，说明应力波速受遗传的控制作用要强于Pilodyn南、北向测定值。3个无损检测指标中，应力波速遗传变异系数最小(仅为3.57%)，低于Ps(8.54%)和Pn(6.92%)，说明日本落叶松Pilodyn南、北向测定值无性系间的选择潜力要高于应力波速。

3.2 日本落叶松无损检测指标与各材质性状的相关关系

Pilodyn南、北向测定值与各材质性状的相关关系比较复杂，其与木材微纤丝角、弹性模量这2类性状相关性不显著； 与木材密度、木材壁腔比2类性状的生长轮均值、早材均值存在显著负相关关系，而与其晚材均值无显著相关。造成此现象的原因可能是早、晚材质地不同，早材的质地疏松，晚材的质地紧密，阻力不同而使它们之间紧密关系不同。

Pilodyn南、北向测定值与木材生长轮密度、早材密度存在极显著负相关关系，此结果与其他学者的研究结果一致(Ishiguri et al.，2008；Knowles et al.，2004；Hansen，2000；Desmond et al.，2010；栾启福等，2011)。Pilodyn南、北向测定值与生长轮密度、早材密度的表型相关系数范围为-0.425～-0.609，小于Wei等(1997)对尾叶桉(Eucalyptus urophylla)、Wang等(1999)对美国黑松(Pinus contorta)、Cown等(1978)对辐射松的研究结果，这可能是由于取样总数较少造成的。Pilodyn南、北向测定值与木材生长轮壁腔比、早材壁腔比亦存在极显著的负相关关系，这为Pilodyn测定值快速预测纤维壁腔比提供了依据。

应力波速与木材密度、木材壁腔比2类性状相关性不显著，这与前人的研究结果有所不同。有研究认为，应力波速和木材密度并不是相互独立的2个性状，它们之间存在紧密的相关关系(Wang et al.，2000；Wu et al.，2011a)。Wu等(2011b)对华南地区3个试验点的51月生杂种桉无性系分析后认为，应力波速与木材密度遗传相关系数分别为0.80，0.87和0.92，二者之间紧密相关。应力波速与木材弹性模量呈紧密正相关，该结论与Chauhan等(2006)对8，16和25年生辐射松的研究结论一致，即不同年龄辐射松应力波速与弹性模量均存在极显著正相关。

本研究中，无损检测指标与材质性状之间的表型相关和遗传相关的强度趋势基本一致，除个别性状外，遗传相关的绝对值大于表型相关的绝对值，证明无损检测指标与材质性状的相关关系主要由参试无性系的遗传性所决定。例如，Pilodyn南、北向测定值与生长轮密度、早材密度、生长轮壁腔比和早材壁腔比4个性状的遗传相关系数绝对值均达0.800以上，应力波速与生长轮微纤丝角、早材微纤丝角、晚材微纤丝角、生长轮弹性模量和早材弹性模量遗传相关系数绝对值均达0.700以上，远高于其表型相关系数绝对值，这为无损检测指标估测材质性状的遗传参数提供了依据。

3.3 日本落叶松无损检测指标和材质性状回归模型的建立

在多元回归模型中，对Pilodyn南向测定值、北向测定值和应力波速贡献率排名前2位的材质性状分别为生长轮密度和早材壁腔比、生长轮密度和生长轮壁腔比以及早材微纤丝角和生长轮弹性模量。为进一步了解无损检测指标对这些性状的预测精度及彼此之间的关系，需建立相应的线性回归模型。Pilodyn南向测定值与生长轮密度间方程的决定系数(0.340)小于Pilodyn北向测定值与生长轮密度的决定系数(0.360)，说明在北向进行测量取得的数据比在南向取得的数据更能反映木材生长轮密度的变化。Pilodyn南向测定值与早材壁腔比、Pilodyn北向测定值与生长轮壁腔比间方程的决定系数分别为0.228和0.245，低于Pilodyn南、北向测定值与木材生长轮密度间方程的决定系数，说明Pilodyn对木材生长轮密度的预测效果要优于纤维壁腔比。应力波速与早材微纤丝角间方程的决定系数为0.348，大于其与生长轮弹性模量间方程的决定系数0.316。

Pilodyn测定值、应力波速与各材质性状间的回归关系均达极显著水平(P＜0.000 1)，说明利用Pilodyn测定值估测生长轮密度、生长轮壁腔比和早材壁腔比，利用应力波速估测早材微纤丝角和生长轮弹性模量是可行的。Pilodyn测定值和应力波速可以有效区分不同群体间木材性状的差异，适合对大量遗传测定林进行比较测定，从而为遗传选择提供依据。然而，Pilodyn测定值和应力波速只是木材性状的间接估测值，不能精确估算每株树木木材性状的具体数值。因此，如何有效结合这2种无损检测方法来预测材性性状还需进一步分析。

参考文献(References)
[1] 黄少伟,谢维辉.2001.实用SAS编程与林业试验数据分析.广州:华南农业大学出版社.(1)
[2] 栾启福,卢萍,井振华,等.2011.Pilodyn评估杂交松活立木的基本密度及其性状相关分析.江西农业大学学报,33(3): 548-552.(1)
[3] 马常耕.1992.落叶松种和种源选择.北京:北京农业大学出版社.(1)
[4] 潘彪,徐朝阳,王章荣.2005.杂交鹅掌楸木材解剖性质及其径向变异规律.南京林业大学学报,29(1): 79-82.(1)
[5] 孙晓梅,张守攻,齐力旺,等.2003.日本落叶松自由授粉家系纸浆材材性遗传变异的研究.林业科学研究,16(5): 515-522.(1)
[6] 王军辉,张守攻,张建国,等.2008.Pilodyn在日本落叶松材性育种中应用的初步研究.林业科学研究,21(6): 808-812.(2)
[7] 张瑛春,王军辉,张守攻,等.2010.Pilodyn和日本落叶松材性指标的关系.林业科学,46(7): 114-119.(3)
[8] 朱景乐,王军辉,张守攻,等.2009a.Pilodyn在日本落叶松活立木材性指标预测中的应用.林业科学研究,22(1): 75-79.(3)
[9] 朱景乐,赵宁,刘增喜,等.2009b.Pilodyn在毛白杨材性指标的预测及选择中的应用.河南农业大学学报,43(4): 376-381.(2)
[10] Cave I D,Walker J C W.1994.Stiffness of wood in fast-grown plantation softwoods: the influence of microfibril angle.For Prod,44: 43-48.(1)
[11] Chauhan S S,Walker J C F.2006.Variation in acoustic and density with age,and their interrelationships in radiation pine.Forest Ecology and Management,229: 388-394.(4)
[12] Cown D J.1978.Comparison of the Pilodyn and torsiometer methods for the rapid assessment of wood density in living trees.Newzealand J For Sci,8: 384-391.(1)
[13] Desmond J,Stackpole R,Vaillancourt M,et al. 2010.Age trends in genetic parameters for growth and wood density in Eucalyptus globulus. Tree Genetics & Genomes,6: 179-193.(1)
[14] Dickson R L,Raymond C A,Joe W,et al. 2003.Segregation of Eucalyptus dunnii logs using acoustics.Forest Ecology and Management,179: 243-251(1)
[15] Evans R,Ilic J.2001.Rapid prediction of wood stiffness from microfibril angle and density.For Prod,51: 53-57.(1)
[16] Greaves B L,Borralho N M G,Raymond C A,et al. 1996.Use of a Pilodyn for the indirect selection of basic density in Eucalyptus nitens. Canadian Journal of Forest Research,26(9): 1643-1650.(4)
[17] Hai P H,Jansson G,Harwood C,et al. 2008.Genetic variation in growth,stem straightness and branch thickness in clonal trials of Acacia auriculiformis at three contrasting sites in Vietnam.Forest Ecology and Management,255: 156-167.(1)
[18] Hansen C P.2000.Application of the Pilodyn in forest tree improvement.DFSC Series of Technical Notes.TN55.Danida Forest Seed Centre,Humlebaek,Denmark.(1)
[19] Hansen J K,Roulund H.1996.Genetic parameters for spiral grain,stem form,Pilodyn and growth in 13 year old clones of Sitka Spruce(Picea sitchensis)(Bong.)Carr.Silvae Genetica,46(2/3): 107-113.(2)
[20] Ishiguri F,Matsui R,Lizuka K,et al. 2008.Prediction of the mechanical properties of lumber by stress-wave velocity and Pilodyn penetration of 36-year-old Japanese larch trees.Original Arbelten.Originals,66: 275-280.(3)
[21] Kien N D,Jansson G,Harwood C,et al. 2010.Clonal variation and genotype by environment interactions in growth and wood density in Eucalyptus camaldulensis at three contrasting sites in Vietnam.Silva Genetica,59(1): 17-28.(1)
[22] Knowles L R,Hansen L W,Wendding A,et al. 2004.Evaluation of non-destructive methods for assessing stiffness of Douglas fir trees.New Zealand Journal of Forestry Science,34(1): 87-101.(2)
[23] Matheson A C,Gapare W J,Ilic J,et al. 2008.Inheritance and genetic gain in wood stiffness in Radiata pine assessed acoustically in young standing trees.Silva Genetica,57(2): 56-64.(1)
[24] Pliura A,Zhang S Y,Bqusquet J,et al. 2006.Age trends in genotypic variation of wood density and its intra-ring components in young poplar hybrid crosses.Ann For Sci,63: 673-685.(1)
[25] Pliura A,Zhang S Y,Mackay J,et al. 2007.Genotypic variation in wood density and growth traits of poplar hybrids at four clonal trials.Forest Ecology and Management,238: 92-106.(1)
[26] Schimleck L R,Michell A J,Raymond C A,et al. 1999.Estimation of basic density of Eucalyotus globulus using near-infrared spectroscopy.Canadian Journal of Forest Research,29: 194-202.(1)
[27] Walker J C F,Butterfield B G.1995.The importance of microfibril angle for the processing industries.N.Z.Forestry,40: 34-40.(1)
[28] Walker J C F,Nakada R.1999.Understanding corewood in some soft-woods: a selective review on stiffness and acoustics.Int For Rev,1: 251-259.(1)
[29] Wang T L,Aitken S N,Rozenberg P,et al. 1999.Selection for height growth and Pilodyn pin penetration in lodgepole pine: effects on growth traits,wood properties,and their relationships.Can J For Res,29: 434-445.(1)
[30] Wang X P,Ross R J,Mcclellan M,et al. 2001.Nondestructive evaluation of standing trees with a stress wave method.Wood and Fiber Science,33(4): 522-533.(1)
[31] Wang X,Ross R J,McClellan M, et al.2000.Strength and stiffness assessment of standing trees using a nondestructive stress wave technique.Res.Pap.FPL-RP-585.U.S.Department of Agriculture,Forest Service,Forest Products Laboratory,Madison,WI,197-206.(2)
[32] Wei X,Borralho N M G.1997.Genetic control of basic density and bark thickness and their relationships with growth traits of Eucalyptus urophylla in south east China.Silvae Genetica,46(4): 245-250.(2)
[33] Wu S J,Xu J M,Li G Y,et al. 2010.Use of the Pilodyn for assessing wood properties in standing trees of Eucalyptus clones.J For Res,21(1): 68-72.(1)
[34] Wu S J,Xu J M,Li G Y,et al. 2011a.Estimation of basic density and modulus of elasticity of eucalypts clones using nondestructive methods in southern China.Journal of Tropical Forestry Science,23(1): 51-56.(1)
[35] Wu S J,Xu J M,Li G Y,et al. 2011b.Genotypic variation in wood properties and growth traits of Eucalyptus hybrid clones in southern China.New Forests,42: 35-50.(1)
[36] Yin Y F,Nagao H,Liu X L,et al. 2010.Mechanical properties assessment of Cunninghamia lanceolata plantation wood with three acoustic-based nondestructive methods.J Wood Sci,56: 33-40.(1)