Non-Destructive Assessment of Wood Sti ﬀ ness in Scots Pine ( Pinus sylvestris L.) and its Use in Forest Tree Improvement

: Wood sti ﬀ ness is an important wood mechanical property that predetermines the suitability of sawn timber for construction purposes. Negative genetic correlations between wood sti ﬀ ness and growth traits have, however, been reported for many conifer species including Scots pine. It is, therefore, important that breeding programs consider wood sti ﬀ ness and growth traits simultaneously. The study aims to (1) evaluate di ﬀ erent approaches of calculating the dynamic modulus of elasticity (MOE, non-destructively assessed sti ﬀ ness) using data from X-ray analysis (SilviScan) as a benchmark, (2) estimate genetic parameters, and (3) apply index selection. In total, we non-destructively measured 622 standing trees from 175 full-sib families for acoustic velocity (VEL) using Hitman and for wood density (DEN) using Resistograph and Pilodyn. We combined VEL with di ﬀ erent wood densities, raw (DEN RES ) and adjusted (DEN RES.TB ) Resistograph density, Pilodyn density measured with (DEN PIL ) and without bark (DEN PIL.B ), constant of 1000 kg · m − 3 (DEN CONST ), and SilviScan density (DEN SILV ), to calculate MOEs and compare them with the benchmark SilviScan MOE (MOE SILV ). We also derived Smith–Hazel indices for simultaneous improvement of stem diameter (DBH) and wood sti ﬀ ness. The highest additive genetic and phenotypic correlations of the benchmark MOE SILV with the alternative MOE measures (tested) were attained by MOE DEN SILV (0.95 and 0.75, respectively) and were closely followed by MOE DEN RES.TB (0.91 and 0.70, respectively) and MOE DEN CONST and VEL (0.91 and 0.65, respectively for both). Correlations with MOE DEN PIL , MOE DEN PIL.B , and MOE DEN RES were lower. Narrow-sense heritabilities were moderate, ranging from 0.39 (MOE SILV ) to 0.46 (MOE DEN SILV ). All indices revealed an opportunity for joint improvement of DBH and MOE. Conclusions: MOE DEN RES.TB appears to be the most e ﬃ cient approach for indirect selection for wood sti ﬀ ness in Scots pine, although VEL alone and MOE DEN CONST have provided very good results too. An index combining DBH and MOE DEN RES.TB seems to o ﬀ er the best compromise for simultaneous improvement of growth, ﬁber, and wood quality traits. acoustic velocity and di ﬀ erent wood density estimates.


Introduction
Wood is a versatile, renewable, and environmentally sustainable material with a wide range of utilization. It has been used as a building material for thousands of years as it provides good insulation, has good machinability, and is exceptionally strong in proportion to its weight [1]. Stiffness, a non-permanent deformation of a sample when a load is applied, represents an important wood mechanical parameter that predetermines the suitability of sawn boards for construction purposes [2]. Since utilization of wood for construction is foreseen to remain important in the future [3], inclusion of wood stiffness into breeding programs appears to be inevitable. It is particularly crucial for species that exhibit adverse negative correlations between growth and stiffness, as the traditional emphasis on volume maximization has been found to result in stiffness degradation, e.g., [4][5][6][7].
Therefore, a reliable technique suitable for rapid screening of large numbers of trees is needed. The destructive bending stress testing, i.e., quantifying stiffness as the static modulus of elasticity, is commonly replaced by acoustic sensing technology, which is based on the determination of the propagation velocity of stress waves induced by a mechanical force [8,9]. In standing trees, acoustic velocity is recorded between two probes hammered into a stem ca 1 m apart and 2-3 cm deep. Wood stiffness, expressed as the dynamic modulus of elasticity (MOE d , GPa) calculated from acoustic velocity (VEL, km·s −1 ) and green wood density (DEN, kg·m −3 ) according to Equation (1) as has proven to be a useful proxy for the static modulus of elasticity [10]. Several different wood density estimates have been tested as potential surrogates for green wood density in the MOE d calculation (Equation (1)). When static MOE (MOE S ) determined from destructive testing or SilviScan MOE (MOE SILV ) estimated through X-ray density combined with X-ray diffraction [11] were used as benchmark variables, a number of studies reported moderate to very strong correlations with MOE d that utilized X-ray density estimated through SilviScan or other X-ray apparatus, e.g., [12][13][14][15] and moderate correlations when volumetric green density was used [16][17][18]. X-ray or volumetric wood densities are, however, rather inconvenient to employ as their estimation requires time-consuming and expensive laboratory measurements and are therefore often replaced by a constant wood density (usually 1000 kg·m −3 ). Such simplification did not compromise the results, as moderate to strong correlations were attained for MOE d using the constant wood density as well, e.g., [6,7,19]. This alternative approach assumes very little variation in wood density [20] and/or a strong relationship between acoustic velocity and wood stiffness. The latter statement is supported by the attainment of moderate to strong correlations of acoustic velocity itself with benchmark variables observed in several studies [19,21,22]. Nevertheless, these correlations often reach even higher levels when some wood density information is included [12,13,23]. Focusing on fast and non-destructive alternatives of wood density measurement, wood density estimated based on Pilodyn penetration depth appears to be ideal as its MOE d gave an almost perfect genetic correlation with the benchmark MOE SILV (r G = 0.99) in a Norway spruce (Picea abies (L.) H. Karst.) study conducted by [12]. Another candidate approach could be density derived from micro-drilling resistance measured by Resistograph, as strong genetic correlations between adjusted Resistograph density and a benchmark density were reported, e.g., for Scots pine (Pinus sylvestris L.) [24] or loblolly pine (Pinus taeda L.) [25,26].
The identification of a reliable method for wood stiffness assessment that would be suitable for screening large numbers of trees is just the first step towards the inclusion of a new trait into an ongoing tree breeding program. Since stem volume is a prioritized trait in many programs, an appropriate method for multi-trait selection has to be applied. Independent culling, which lies in setting independent truncation thresholds for each trait of interest [27], tandem selection, which considers only a single trait in each breeding cycle [28], or two-stage selection, which focuses on one trait at a time but covers all traits of interest during one breeding cycle [29], are possible options to choose from. Nevertheless, the most efficient method designed for simultaneous improvement of several traits is index selection, which is performed through assigning a weight for each trait according to its economic importance [30,31]. Index selection is suitable for traits that differ in economic importance, heritability, variability, and/or correlation among their phenotypic and genotypic values [27], particularly if the correlation is negative [32]. Estimation of economic weights for the index selection is, however, very complicated, mainly due to the uncertain future of end-products and prerequisites for their production as well as unwillingness of the wood processing industry to fully disclose their financial flows. Nevertheless, economic weights were successfully estimated and applied for radiata pine (Pinus radiata D. Don) structural timber in Australia [33].
The aim of this study is to (1) evaluate different approaches of calculating the dynamic modulus of elasticity (MOE d ) in Scots pine using SilviScan-, Pilodyn-, Resistograph-based, and constant wood density employing SilviScan modulus of elasticity as a benchmark, (2) estimate heritabilities of the best MOE d approaches and calculate their phenotypic and genetic correlations with growth, fiber, and wood quality traits, and (3) apply index selection using non-monetary weights to explore possibilities of simultaneous improvement of stem diameter and wood stiffness under the presence of unfavorable genetic correlation.

Test Material
In total, 622 trees from 175 families of 44 parents were sampled in a Scots pine full-sib progeny test "Grundtjärn" (#S23F 711261, lat. 63.5556 • N, long. 17.4139 • E, alt. 320 m, area 3.5 ha) located in central Sweden. The test was established by the Forestry Research Institute of Sweden (Skogforsk) in 1971 on a slightly west-sloping sandy moraine with the site productivity around 3 m 3 /ha/year [34]. It was comprised of 7240 seedlings representing 179 full-sib families of 45 parents that were crossed following a partial diallel mating design, where each parent participated 8-9 times as either a mother or a father. Phenotypically selected parental trees originated from forest populations between latitudes 63 • N and 64 • N. One-year-old seedlings were planted with 2.2 × 2.2 m spacing following a completely randomized single tree plot design. The experiment was divided into 181 post-blocks, each consisting of 40 trees (4 columns by 10 rows) and was systematically thinned at age 35.

Wood Density Assessment
Wood density of standing trees was assessed non-destructively using micro-drill Resistograph IML-RESI PD 300 (Instrumenta Mechanic Labor, Germany) and penetrometer Pilodyn 6J Forest (PROCEQ, Switzerland) in early summer at tree ages of 45 and 47, respectively. Both tools utilize a slender needle that penetrates into wood and leave just negligible openings. Each tree was drilled with Resistograph bark-to-bark and shot with Pilodyn through bark (DEN PIL ) and after removing a patch of bark (DEN PIL.B ) ca 1.3 m above ground. All wood density measurements were performed with an effort to avoid knots or any visible damages that could cause a bias in wood density estimates. Resistograph's drilling profiles (resistograms) were adjusted (detrended and debarked) according to [24]. Mean values of unadjusted and adjusted Resistograph's records (DEN RES and DEN RES.TB , respectively) and reciprocals of Pilodyn's readings were used in this study.

SilviScan Data
Pith-to-bark increment cores of 10 mm in diameter, extracted ca 1.3 m above ground, were cut into thin strips and analyzed with SilviScan technology (CSIRO, Australia). SilviScan wood stiffness (MOE SILV ) derived from X-ray absorption and diffraction was used as a benchmark for evaluation of six different dynamic moduli of elasticity calculated according to Equation (1). Other traits estimated by SilviScan that are used in this study include: Wood density (DEN SILV ); density of earlywood (EWD), transition wood (TWD), and latewood (LWD); proportion of earlywood (EWP), transition-wood (TWP), and latewood (LWP); microfibril angle (MFA); fiber wall thickness (FWT); fiber coarseness (FCS); and fiber width in radial (FRW) and tangential (FTW) direction. Earlywood, transition wood, and latewood were defined as annual ring segments with densities ranging from 0-20%, 20-80%, and 80-100% of the span from minimum to maximum density within the ring, respectively. All traits measured by SilviScan were obtained for each annual ring; mean values (excluding EWP, TWP, and LWP) were calculated as area-weighted values (AWV), where d i is a value for ith annual ring with cross-sectional area α i [35].

Wood Stiffness Assessment through Acoustic Velocity
Standing-tree acoustic velocity (VEL) was measured in autumn at a tree age of 42 by Hitman ST300 (Fibre-gen, New Zealand) between two probes hammered 0.7-1.3 m apart into the southeastern side of a stem. Time-of-flight of a stress wave induced by tapping the transmitter probe with a steel hammer was recorded and, knowing the exact distance between the two probes, acoustic velocity was calculated. Mean value of two series of eight taps was used for calculating the dynamic modulus of elasticity (MOE d , Equation (1) reflecting the wood stiffness. Different wood densities were used for MOE d calculation: The VEL and MOE d estimates were compared with SilviScan modulus of elasticity (MOE SILV ), which served as a benchmark. Pairwise comparisons between MOE means were performed using Duncan's multiple range test (duncan.test in R program [36]) at the significance level of 0.01. Variation in wood quality traits due to measurements being taken in different seasons and/or at different tree ages was considered negligible.

Growth Data
Height (HGT) and diameter at breast height (DBH) were obtained from Skogforsk. Stem volume (VOL) was calculated as a function of height and diameter [37]:

Statistical Analysis
Using statistical package ASReml 4 [38], the response variables were fitted into the following linear model in order to estimate variance and covariance components: where y ijkl is the variable measurement for lth offspring from full-sib family of parents j and k growing in ith block; µ is the variable overall mean; B, P, F, and e are random effects of block, parent, family, and residuals, respectively. Individual-tree narrow sense heritability (h 2 i ) for each variable was estimated using variance components from the univariate analysis as where σ 2 A , σ 2 P , σ 2 p , σ 2 f , and σ 2 e are variances for additive genetic, phenotypic, parental, family, and residual components, respectively. Taylor series expansion was employed to estimate standard errors [38]. Phenotypic and genetic correlations (r xy ) were calculated as where σ 2 x and σ 2 y are phenotypic or additive genetic variance components for trait x and y, respectively, and σ xy is phenotypic or additive genetic covariance component between traits x and y estimated by fitting a bivariate model (Equation (4)) [38]. Dendrogram of additive genetic correlations depicting relationships among traits was produced by hierarchical cluster analysis based on dissimilarity matrix using hclust function in R program [36].
Index selection was applied to explore the possibility of simultaneous improvement of stem diameter and wood stiffness, in spite of the negative genetic correlation between the two traits. Index selection was also constructed for DBH and wood density as a close proxy for wood stiffness. Different sets of weights for stiffness (GPa) and density (kg·m −3 ), ranging from 0 to 3.6 and from 0 to 0.12, respectively, relative to 1 for DBH (cm) were used. The Smith-Hazel index, which maximizes the efficiency of selection through accounting for heritabilities and genetic and phenotypic correlations among traits, was employed to find optimal weights with desired genetic gains for both traits [30,39].
The index (I) is generally defined as where P s represent phenotypic performance for each trait and b s are their corresponding index coefficients. Index coefficients were calculated as where P is phenotypic variance-covariance matrix, A is additive genetic variance-covariance matrix, and w is vector of weights for each trait. The expected genetic gain (∆A x ) for trait x included in the Smith-Hazel index was estimated following [32] as: where i is the selection intensity (1% refers to i = 2.67), σ 2 A x is additive genetic variance of trait x, σ A xy is additive genetic covariance for traits x and y, and σ I is a square root of phenotypic variance of the index calculated as Genetic gain (G A x ) for direct selection was estimated as where i is selection intensity, h 2 x is heritability for trait x, and CV x is coefficient of variation for trait x calculated as a phenotypic standard deviation divided by mean. Correlated response (CR y ) to selection for a target trait y was determined as CR y = ih x h y r xy CV y (12)  where h x and h y are squared roots of narrow sense heritabilities for selection trait x and target trait y, respectively, r xy is additive genetic correlation between trait x and y, and CV y is coefficient of variation for target trait y. Traits involved in the index selection were included as a new aggregated variable x calculated according to Equation (7). Weights corresponding to intersections of plotted weights versus expected genetic gains were used for calculation of index coefficients.

Variation in Different Wood Stiffness (MOE) Estimates
Different estimates of wood stiffness (MOE) together with their descriptive statistics and narrow-sense heritabilities are summarized in Table 1

Heritability
All individual-tree narrow-sense heritabilities obtained for different estimates of MOE and VEL were moderate ( Table 1). Heritability of the benchmark MOE SILV (0.39) was higher than that of VEL and MOE DEN CONST (both 0.37) but lower than those for all other estimates of MOE d . The highest heritability was attained by MOE DEN SILV (0.46), followed by MOE DEN RES.TB (0.45) and MOE DEN RES (0.44).

Comparison of Moduli of Elasticity
Using MOE SILV as the benchmark, MOE DEN SILV gave the best estimate of modulus of elasticity (r G = 0.95 and r P = 0.75; Table 2, Figure 1). The second-best estimate of MOE was, at the genetic level, attained by VEL, MOE DEN RES.TB , and MOE DEN CONST (r G = 0.91 for all); however, at the phenotypic level, MOE DEN RES.TB gave a better estimate (r P = 0.70 and r P = 0.65, respectively). It means that additional information in the form of adjusted Resistograph density improved the MOE estimate, although the improvement was just minor. On the other hand, neither the two Pilodyn densities nor the unadjusted Resistograph density improved the MOE estimation. Their correlations with the benchmark MOE SILV were the same or lower compared with those of MOE SILV with VEL where no additional wood density information was supplied.

Phenotypic and Additive Genetic Correlations of Selected Moduli of Elasticity with Other Wood, Fiber, and Growth Traits
Phenotypic and additive genetic correlations of the benchmark MOESILV, VEL, and the best performing moduli of elasticity estimated based on VEL (MOEDEN CONST , MOEDEN SILV , MOEDEN RES.TB ) are presented in Table 3. Correlations of VEL and MOEDEN CONST were almost identical. Genetic correlations among moduli of elasticity and different wood density estimates were moderate to strong, while phenotypic correlations were moderate to weak. The strongest relationships with wood density traits were obtained for MOEDEN RES.TB , followed by MOEDEN SILV . On the other hand, relationships with VEL and MOEDEN CONST were the weakest. Genetic correlations of VEL and moduli of elasticity with LWP, TWP, and EWP were moderately positive, weakly negative, and close to zero, respectively. MOESILV was strongly negatively correlated with MFA. Genetic correlations of VEL and MOEDEN CONST with MFA were also strong and negative; other correlations with MFA were moderate.
Strong genetic and moderate phenotypic positive correlations were obtained for moduli of elasticity with FTW, except for MOEDEN CONST that had weaker correlations. Genetic and phenotypic correlation with FCS were moderate, except for weak phenotypic correlation of MOEDEN CONST . Correlations with FRW and FTW were weak and negative for all moduli of elasticity. Genetic and phenotypic correlations of VEL and moduli of elasticity with growth traits were moderately negative, except for week phenotypic correlations with HGT. MOESILV produced the strongest negative correlations with growth traits.
Relationships among variables based on additive genetic correlations are shown in Figure 2. The dendrogram revealed two distinct groups: One consisting of two sub-groups (growth traits; and earlywood and transition wood proportion together with fiber width) and the other consisting of four sub-groups (Pilodyn density along with latewood proportion; moduli of elasticity, acoustic velocity,   Table 3. Correlations of VEL and MOE DEN CONST were almost identical. Genetic correlations among moduli of elasticity and different wood density estimates were moderate to strong, while phenotypic correlations were moderate to weak. The strongest relationships with wood density traits were obtained for MOE DEN RES.TB , followed by MOE DEN SILV . On the other hand, relationships with VEL and MOE DEN CONST were the weakest. Genetic correlations of VEL and moduli of elasticity with LWP, TWP, and EWP were moderately positive, weakly negative, and close to zero, respectively. MOE SILV was strongly negatively correlated with MFA. Genetic correlations of VEL and MOE DEN CONST with MFA were also strong and negative; other correlations with MFA were moderate. Strong genetic and moderate phenotypic positive correlations were obtained for moduli of elasticity with FTW, except for MOE DEN CONST that had weaker correlations. Genetic and phenotypic correlation with FCS were moderate, except for weak phenotypic correlation of MOE DEN CONST . Correlations with FRW and FTW were weak and negative for all moduli of elasticity. Genetic and phenotypic correlations of VEL and moduli of elasticity with growth traits were moderately negative, except for week phenotypic correlations with HGT. MOE SILV produced the strongest negative correlations with growth traits. Relationships among variables based on additive genetic correlations are shown in Figure 2. The dendrogram revealed two distinct groups: One consisting of two sub-groups (growth traits; and earlywood and transition wood proportion together with fiber width) and the other consisting of four sub-groups (Pilodyn density along with latewood proportion; moduli of elasticity, acoustic velocity, and earlywood density; SilviScan, Resistograph, transition wood and latewood density, plus fiber wall thickness; and microfibril angle with fiber coarseness). and earlywood density; SilviScan, Resistograph, transition wood and latewood density, plus fiber wall thickness; and microfibril angle with fiber coarseness).

Index Selection
Index selection, with different sets of weights for MOE and DEN relative to 1 for DBH, was employed in order to evaluate the possibility of simultaneous improvement. For all studied combinations of DBH and a wood quality trait, weights plotted against expected genetic gains (Figure 3 density) compared to the pair with DEN SILV (0.8 cm for DBH and 16 kg·m −3 for density). Weights corresponding to the graphs' intersections were used for calculating index coefficients and a subsequent analysis of the correlated response to selection. between weights 0.2 and 3.2, with the highest expected gain being attained at the intersection (1.3 cm for DHB and 1.3 GPa for MOE). Pairs with MOEDEN SILV and MOEDEN RES.TB showed an almost identical pattern, with a positive region between weights 0.4 and 2.8 and the intersectional gain of 1. The graph of simultaneous selection of DBH and MOESILV appeared to have the shortest region of positive values (0.5-1.8) but the intersectional expected gain was 1 as well. The pair of DBH with DENSILV had a slightly longer positive region (0.035-0.075) than the one with DENRES.TB (0.030-0.065). Nevertheless, the pair with DENRES.TB reached a higher genetic gain for both traits (1 cm for DBH and 20 kg·m -3 for density) compared to the pair with DENSILV (0.8 cm for DBH and 16 kg·m -3 for density). Weights corresponding to the graphs' intersections were used for calculating index coefficients and a subsequent analysis of the correlated response to selection.

Response to Selection
Correlated response to selection for important growth (DBH, VOL), fibre (MFA), and wood quality (MOESILV, DENSILV) traits is summarized in Table 4. Selection for DBH separately resulted in positive genetic gains in growth traits but in negative gains in wood quality traits. Note that a positive gain for MFA means an unfavorable increase in microfibril angle. Contrariwise, selection based on MOESILV led to improvement of wood quality and fiber traits but to a decline in growth traits. Therefore, indices combining DBH with different MOE and DEN estimates have been employed. All indices resulted in desirable simultaneous improvement of both growth and wood quality traits as well as MFA. The highest genetic response of growth traits was attained by selection based on DBH and MOEDEN CONST , followed by DBH and MOEDEN SILV and DBH and DENSILV, while the highest response of stiffness was achieved by DBH and MOEDEN RES.TB and DBH and MOESILV. Selection based on DBH and DENRES.TB and DBH and MOEDEN RES.TB resulted in the highest gain for wood density. The

Response to Selection
Correlated response to selection for important growth (DBH, VOL), fibre (MFA), and wood quality (MOE SILV , DEN SILV ) traits is summarized in Table 4. Selection for DBH separately resulted in positive genetic gains in growth traits but in negative gains in wood quality traits. Note that a positive gain for MFA means an unfavorable increase in microfibril angle. Contrariwise, selection based on MOE SILV led to improvement of wood quality and fiber traits but to a decline in growth traits. Therefore, indices combining DBH with different MOE and DEN estimates have been employed. All indices resulted in desirable simultaneous improvement of both growth and wood quality traits as well as MFA. The highest genetic response of growth traits was attained by selection based on DBH and MOE DEN CONST , followed by DBH and MOE DEN SILV and DBH and DEN SILV , while the highest response of stiffness was achieved by DBH and MOE DEN RES.TB and DBH and MOE SILV . Selection based on DBH and DEN RES.TB and DBH and MOE DEN RES.TB resulted in the highest gain for wood density. The highest decrease in microfibril angle was achieved by selecting for DBH and MOE SILV , followed by MOE DEN SILV and MOE DEN RES.TB , whereas selection combining DBH and DEN had nearly no effect on MFA. Weights (in brackets) correspond to intersections of plotted weights versus expected genetic gains ( Figure 3) and were used for calculating index coefficients.

Estimation of Wood Stiffness
Wood stiffness is an important mechanical parameter determining the suitability of wood for construction purposes. It is, however, negatively correlated with growth traits, which means that a systematic selection for bigger trees leads to deterioration in wood stiffness. In order for wood stiffness to be included in forest tree breeding programs, a reliable method for its assessment, suitable for fast screening of thousands of trees, is needed. Hence, we have evaluated several different approaches of calculating modulus of elasticity from acoustic velocity (Equation (1), Table 2) measured non-destructively on standing trees. Having MOE SILV as a benchmark, MOE DEN SILV , which combines acoustic velocity and wood density derived from X-ray absorption of SilviScan, was found to be the best proxy, providing the highest correlations (r G = 0.95 and r P = 0.75). It is in good accordance with other studies-for instance, strong correlations between MOE SILV and MOE DEN SILV were reported for radiata pine (r P = 0.97) [15] or Norway spruce (r G = 1.01 and r P = 0.66) [12]. This relationship can, however, be slightly overestimated as SilviScan wood density was used for calculating both MOE SILV and MOE DEN SILV . Nevertheless, employing SilviScan or another X-ray-based technology is feasible only for a small-scale assessment, as the wood sample preparation and scanning are laborious and costly. Of the proxies suitable for inclusion in operational breeding programs, the best MOE estimate incorporating acoustic velocity and Resistograph adjusted density (MOE DEN RES.TB ) was closely followed by the MOE estimate based solely on acoustic velocity (MOE DEN CONST ) and acoustic velocity itself (VEL). The MOE DEN RES.TB had a slightly higher phenotypic correlation (r P = 0.75 versus 0.65) with the benchmark MOE SILV but genetic correlations of MOE DEN RES.TB , MOE DEN CONST and VEL were identical (r G = 0.91). Similar correlations between MOE SILV and MOE DEN CONST were obtained, e.g., for lodgepole pine (Pinus contorta Douglas ex Loudon) (r G = 0.91 and 0.90, r P = 0.67 and 0.61) [6]; a strong phenotypic correlation of VEL with a benchmark MOE was published for Scots pine (r P = 0.72) [19] or Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) (r P = 0.79) [13]. Genetic correlations of all other MOE estimates with MOE SILV in this study were lower than the one of VEL, which indicates that the additional information on wood density measured by Pilodyn, either with or without bark, or by Resistograph, without the subsequent adjustment of drilling profiles, did not bring any improvement to the MOE estimation. On the other hand, a combination of Pilodyn wood density and acoustic velocity (MOE DEN PIL ) was identified as an efficient method for indirect selection for MOE in Norway spruce (r G = 0.99, r P = 0.53) [12].

Correlations
Additive genetic correlations of MOE and VEL with DEN SILV were strong and positive in this study, which corresponds with correlations of MOE SILV and MOE DEN PIL reported for Norway spruce [12]. It suggests that selection based on one of the traits only, either MOE or DEN, whichever is more convenient to measure, would lead to improvement of both traits at the same time. On the other hand, weak to moderate positive correlations of VEL or MOE DEN CONST with a benchmark density were observed for other conifers, e.g., [7,13,42]. Strong negative genetic correlations of MOE and VEL with MFA were in good accordance with a number of earlies studies, e.g., [6,12,[42][43][44], implying that selection for MOE or just VEL would result in a desirably lower MFA. Negative correlations of MOE and VEL with growth traits were weak to moderate, of which the weakest were those with HGT. Similar results were obtained, e.g., for Douglas-fir [13] or Norway spruce [12]. Moderate negative genetic correlations between MOE DEN CONST and DBH were also reported for Sitka spruce (Picea sitchensis (Bong.) Carr.) and lodgepole pine, whilst negative genetic correlations between MOE SILV and DBH and between MOE S and DBH were strong in the two species [7] and [6], respectively. In contrast, weakly positive genetic correlations of VEL with HGT and DBH were estimated for black spruce [40] and white spruce (Picea glauca (Moench) Voss) [42].

Factors Affecting Acoustic Velocity Measurements
A number of factors have been reported to affect VEL measurements on standing trees such as the presence of knots, compression wood [45], spiral grain [46], age [47], temperature [48], moisture content [49][50][51], or direction of the measurement on sloping terrain [52]. It has also been observed that VEL can be influenced by insertion depth of probes [53], their imbedding angle [54], and variation in hitting intensity [45]. Moreover, VEL measured on standing trees considers only the outermost wood, which is usually stiffer than the whole log [55]. Nevertheless, in this study, VEL has proven to be an efficient method for wood stiffness assessment in the whole stem profile, being capable of revealing high MFA too.

Response to Selection
Taking into account the negative correlations between growth and wood quality traits, it is obvious that selection for growth results in a lower wood quality and vice versa. Moreover, prioritizing growth traits increases MFA, which negatively influences the quality of end-products, e.g., stiffness of sawn boards [56] or pulp yield [57]. A similar pattern of the expected genetic response after single-trait selection was also observed in lodgepole pine [6] or Norway spruce [12,58].
In order to explore opportunities for joint improvement, we have constructed indices consisting of DBH as an easy-to-measure growth trait and several estimates of MOE and DEN. Different sets of weights for MOE and DEN relative to 1 for DBH were plotted against expected genetic gains ( Figure 3) whereby regions with positive gains for both traits were revealed. Intersections of the two respective curves on each graph were considered as points of a good compromise between the potential gain increase in one trait and decline in the other when each of the traits is selected alone. Weights corresponding to these intersections were used for calculating index coefficients. Weights for the index selection of DBH and MOE SILV determined based on the intersection of their plotted expected genetic gains (1 cm for DHB and 1 GPa for MOE) correspond with weights suggested by [58] for Norway spruce and applied by [59], who studied lodgepole pine growing in northern Sweden (note that DBH in the two studies was expressed in mm, whereas we use cm).
The performance of indices calculated using SilviScan-based estimates of wood stiffness and density (MOE SILV and DEN SILV , respectively) was compared with that of indices calculated using non-destructive techniques for wood quality assessment (acoustic velocity and drilling resistance, Table 4). An index combining DBH and MOE DEN CONST , representing the simplest approach to MOE estimation, resulted in the highest gain in growth but lowest in wood density, which was anticipated as no density information was supplied. This index also gave a rather low gain in wood stiffness compared to other indices that involved MOE. The index incorporating DBH and mean wood density calculated from adjusted drilling profiles (DEN RES.TB ) performed best for wood density improvement only but genetic gains for all other traits were low. The combination of DBH and MOE DEN RES.TB calculated from acoustic velocity and adjusted drilling resistance seems to be the best option, as selection based on this index generated the highest and second highest gain for stiffness and wood density, respectively. Gains for growth traits were higher compared to the benchmark DBH and MOE SILV but lower than those using indices with MOE DEN SILV and MOE DEN CONST . Correlated response of MFA was comparable with DBH and MOE DEN SILV but a little lower than the benchmark DBH and MOE SILV .
The joint selection for DBH and MOE DEN RES.TB appears to offer a reasonable compromise in terms of the attained genetic gain for growth and wood quality traits as well as for MFA, resulting in trees with wood of versatile properties suitable for a great variety of purposes.

Practical Implications
Apart from the accuracy of assessment and purchasing expenses, the time required for conducting the measurements should also be considered when a new technique is evaluated for potential utilization. The Hitman ST300 operated by a team of two people is more than twice as efficient as if it is operated by a single person. However, the IML Resistograph can be efficiently operated by a single person while the measurements are acquired ca 2-3 times faster than those by Hitmen when operated by two people. Furthermore, while the measurement speed of the Resistograph depends on stem diameter only, the speed of Hitman, although constant under ideal conditions, is strongly affected by the presence of branches reaching below the upper probe, which must be removed prior to the measurements. The necessity to repeat measurements is comparable for both tools. Post-measurement adjustment of Resistograph's drilling profiles is recommended in order to get more accurate estimates of wood density [25]. A linear detrending, followed by bark removal proposed by [24], represents a simple and quick fix. When one has to decide for just one of the wood quality assessment methods (Hitman or Resistograph) due to shortage of funding, time, or people, Hitman appears to a better choice than Resistograph, although the measurements will require a considerably higher amounts of time. Nevertheless, measuring with the Resistograph is fast and the inclusion of wood density estimated by the Resistograph in MOE calculations brings a substantial improvement in genetic gains for wood quality traits.

Conclusions
In this study, we focused on wood stiffness as an important mechanical parameter. We evaluated several non-destructive approaches of calculating the dynamic modulus of elasticity (MOE d , non-destructively assessed wood stiffness) based on acoustic velocity and wood density measurements. The following wood density measures were considered: Constant density, SilviScan density, unadjusted and adjusted Resistograph density, and Pilodyn density measured with and without bark. MOE d calculated from acoustic velocity and wood density derived from Resistograph's adjusted drilling profiles (MOE DEN RES.TB ) appears to be the most efficient combination for indirect selection for wood stiffness. Nevertheless, acoustic velocity alone and MOE d derived from acoustic velocity and constant wood density reached nearly the same results as MOE DEN RES.TB . In order to investigate the possibility of a joint improvement of unfavorably correlated growth and wood quality traits, we applied index selection. An index that combines stem diameter and MOE DEN RES.TB seems to offer the best compromise for simultaneous improvement of growth, fiber, and wood quality traits.