Non-destructive wood density assessment of Scots pine (Pinus sylvestris L.) using Resistograph and Pilodyn

Abstract

We tested two methods for non-destructive assessment of wood density of Scots pine standing trees: one based on penetration depth of a steel pin (Pilodyn) and the other on micro-drilling resistance (Resistograph). As a benchmark we used wood density data from x-ray analysis (SilviScan). We assessed in total 622 trees of 175 full-sib families growing in a single progeny test. Pilodyn was applied with bark (PIL) and without bark (PILB). Raw Resistograph drilling profiles (RES) were adjusted (RESTB) in order to eliminate increasing trend caused by needle friction. Individual narrow-sense heritability of benchmark SilviScan density (DEN; 0.46) was most closely approached by that of adjusted RESTB (0.43). Heritabilities were lower for unadjusted RES (0.35) as well as for PIL and PILB (both 0.32). Additive genetic correlations of the benchmark DEN with RES, RESTB, PIL and PILB were 0.89, 0.96, 0.59 and 0.71, respectively. Our results suggest that Resistograph is a more reliable tool than Pilodyn for wood density assessment of Scots pine; however, we highly recommend adjusting Resistograph drilling profiles prior to further analyses.

About the Authors:

Irena Fundova

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Software, Validation, Visualization, Writing – original draft

Affiliations Umeå Plant Science Centre, Department of Forest Genetics and Plant Physiology, Swedish University of Agricultural Sciences, Umeå, Sweden, Skogforsk (The Forestry Research Institute of Sweden), Sävar, Sweden

ORCID logo http://orcid.org/0000-0002-3817-1599

Tomas Funda

Roles Conceptualization, Investigation, Methodology, Validation, Writing – review & editing

Affiliation: Umeå Plant Science Centre, Department of Forest Genetics and Plant Physiology, Swedish University of Agricultural Sciences, Umeå, Sweden

Harry X. Wu

Roles Funding acquisition, Project administration, Supervision, Writing – review & editing

* E-mail: [email protected]

Affiliation: Umeå Plant Science Centre, Department of Forest Genetics and Plant Physiology, Swedish University of Agricultural Sciences, Umeå, Sweden

Introduction

Scots pine (Pinus sylvestris L.) is a commercially significant forest tree species native to Eurasia. It is mainly used for production of sawn timber, pulp, and furniture [1]. In Sweden, Scots pine breeding started in 1950s with phenotypic selection of superior trees, followed by establishment of seed orchards. Emphasis was mainly put on growth, vitality, branching habit and stem straightness, whereas wood quality traits such as density or stiffness were not considered [2]. Unfortunately, negative relationship between growth and wood quality traits has been reported in a number of conifer tree species including Scots pine [3, 4], radiata pine [5], maritime pine [6], black spruce [7] or Douglas-fir [8] and it is therefore of vital importance to include wood quality traits in breeding programs too. This task brings about a need to find reliable tools for non-destructive wood quality assessment that would be capable of a rapid evaluation of a large number of standing trees.

Wood density is considered to be the best single predictor of wood quality. It is well correlated with other wood quality traits such as strength or stiffness [9] and it significantly affects wood suitability for different end uses. For instance, higher wood density is more appropriate for constructional lumber, is associated with greater pulp yield [6] and increases wood machinability, particularly boring, shaping and mortising [10].

Wood density is defined as the mass of wood per unit volume at a given moisture content. It can be accurately determined by either traditional volumetric method [11] or x-ray densitometry [12]. X-ray densitometry combined with x-ray diffraction and image analysis is incorporated in SilviScan technology, which enables efficient measuring of numerous wood and fiber properties such as wood density, stiffness, microfibril angle or fiber dimensions [13]. However, both the volumetric and x-ray methods are enormously time, cost and labor intensive and, in addition, require extracted increment core samples; thus, they cannot be classified as truly non-destructive.

Several, to a lesser or greater extent, non-destructive methods have been developed for rapid field wood density assessment. They include quantification of 1- torque of a borer (torsiometers), 2- withdrawal resistance of a nail, 3- penetration depth of a pin (penetrometers) or 4- micro-drilling resistance (resistometers) [14]. Penetrometers and micro-drilling resistometers use a thin steel probe that penetrates into wood, leaving just a slight hole. They were originally developed for testing the quality of wooden structures, but have proven to be suitable for wood density assessment of standing trees as well. The Pilodyn penetrometer measures penetration depth of a spring-loaded blunt pin (Ø 2 mm) that is shot into the wood with an exact force. Its records give strong genetic correlations with wood density in some conifer species [15–18]; however, the penetration depth is limited to a few centimeters of the outer wood (2–3 cm) and, hence, no information is provided about the rest of a tree’s profile. Moreover, it is also worthwhile to consider removing bark prior the measurement [19], as the bark may constitute a substantial portion of the assessed profile and its uneven thickness and lower density might affect the estimates. On the other hand, the Resistograph micro-drill accounts for the whole stem’s profile at a given height, as it measures drilling resistance, i.e., energy needed for a drill needle (Ø 3 mm) to penetrate wood, from bark to bark at a constant speed [20]. A resistogram, a visualized profile of drilling resistance plotted against penetration depth, reveals density variation inside a stem induced by alternating earlywood and latewood as well as by the presence of knots, cavities, or wood in various stages of physical deterioration. Resistograms closely resemble x-ray density profiles; they have similar sensitivity although their local resolution is a little lower. Their drawback is that they often exhibit an increasing trend caused by accumulated needle friction [21]; therefore, they should be detrended to eliminate bias in wood density estimates [12]. The necessity of such post-measurement adjustment is probably the main reason why the Resistograph has not been widely used for wood density assessment in forest genetics field tests [14]. However, strong genetic correlations between wood density and mean Resistograph density have been reported by several studies [8, 12, 22], implying that the Resistograph could be a suitable tool for wood density assessment also in Scots pine.

The aim of this study is to (1) develop an algorithm suitable for Scots pine to determine mean density from a Resistograph’s drilling profile; (2) evaluate the reliability of Resistograph and Pilodyn measurements for wood density assessment in Scots pine using SilviScan data as a benchmark; (3) estimate inheritance of Resistograph and Pilodyn wood density measurements; and (4) calculate their phenotypic and genetic correlations with growth and wood quality traits.

Materials and methods

Test material

The study was conducted in a Scots pine (Pinus sylvestris L.) full-sib progeny test #S23F 711261 “Grundtjärn” (lat. 63.5556° N, long. 17.4139° E, alt. 320 m, area 3.5 ha) located in central Sweden. The test was established for research purposes by The Forestry Research Institute of Sweden (Skogforsk) in 1971 on silty moraine using completely randomized single tree plot design. It comprised of 7,240 trees representing 179 full-sib families that were generated with a partial diallel mating design using 45 parents. Geographic origin of the parents (plus trees) and mating design for this progeny test was described in [3]. The test site was divided into 181 post-blocks to ease orientation and remove environmental effects during data analyses, and each post-block consisted of 40 trees (4 columns by 10 rows with spacing of 2.2 m in both directions). In total 622 trees representing 175 families of 44 parents were included in this study.

Skogforsk, an organization that has the authority to grant permissions to access and maintain this progeny test as well as collect, analyze and publish any data generated therein, approved this study; no other permissions were required. No protected or endangered species were involved.

In situ wood density measurement

The micro-drill Resistograph IML-RESI PD300 (Instrumenta Mechanic Labor, Germany) and penetrometer Pilodyn 6J Forest (PROCEQ, Switzerland) were used for estimating wood density of standing trees. The Resistograph was used for drilling trees bark to bark in two mutually perpendicular directions (from southeast and from southwest) at height of ca 1.2 m and special attention was paid to avoiding drilling through knots or visible stem damages. Each profile was checked immediately after drilling on the tool’s screen and the measurement was repeated when necessary. Pilodyn was applied on each stem with bark and without bark in the southwest direction at ca 1.3 m above ground. Both the Resistograph and Pilodyn (with bark) were applied on fissures between bark scales to ensure that the tools’ tips were stable during measurements and to avoid high proportion of bark in the records. The reciprocal of Pilodyn’s penetration depth and an average value of adjusted Resistograph’s records were used as indirect wood density estimates.

Processing of resistograms

The Resistograph IML-RESI PD300 is capable of measuring trees up to 30.2 cm in diameter with a drilling resolution of 100 points per centimeter. Each radial measurement generates a resistogram that includes lower-density bark on both sides and a pith in the center (Fig 1). The overall shape of a resistogram also reveals differences in juvenile and mature wood densities, resulting in a bowl-like shape in the center passing into a plateau or slightly decreasing at the edges, which is typical for conifers [21]; a smooth curve fitted through a scatter plot of SilviScan mean ring density and cambial age makes the natural shape of wood density variation more clearly visible (Fig 2). Aside from the natural shape, most of the Resistograph profiles exhibited an increasing trend caused by accumulated needle friction (Fig 1A).

[Figure omitted. See PDF.]

Fig 1.

Unadjusted (A), detrended (B), and detrended & debarked (C) Resistograph profiles.

https://doi.org/10.1371/journal.pone.0204518.g001

[Figure omitted. See PDF.]

Fig 2. SilviScan mean ring density plotted against cambial age.

https://doi.org/10.1371/journal.pone.0204518.g002

In order to eliminate bias in wood density estimates, each profile with this issue was detrended, assuming a linear trend. A slope of increase was calculated for each profile separately as(1)where and are Resistograph density numbers at the tail and head of a raw profile, respectively, and Lmax is the total length of the profile. Then every density number point (100 cm-1) along the whole profile was adjusted as(2)where and are Resistograph density numbers of the raw and detrended profiles at point i, respectively, and Li is the length of the profile from the beginning till position i (Fig 1B). When a Resistograph profile was incomplete, whether due to the stem diameter being too large or due to an abrupt profile ending, the density number was detrended based on the overall-site average slope. Bark, together with the first unfinished growth ring, was removed from each detrended profile in order to obtain a more accurate wood density estimate (Fig 1C). The average density number (Resistograph wood density, RES) was calculated as the arithmetic mean of all density numbers DNi along a profile. Resistograph wood density was calculated for each raw, detrended, and detrended & debarked profile and for each of the two drilling directions separately (the southeastern and southwestern sides being marked as A and B, respectively) and combined.

SilviScan data

Bark-to-bark increment cores with diameter of 10 mm were taken at 1.3 m from southeastern side of each stem in 2011. Prior to the SilviScan (CSIRO, Australia) analysis, pith-to-bark radial strips (2 mm thick and 7 mm wide) were sawn from the cores, soaked in acetone to remove extractives and air-dried under laboratory conditions (23°C and 43% relative humidity). The following traits, obtained from SilviScan analysis, were used in this study: mean wood density (DEN); density of earlywood (EWD), transition wood (TWD) and latewood (LWD); proportion of earlywood (EWP), transition-wood (TWP) and latewood (LWP); microfibril angle (MFA); static modulus of elasticity (MOEs); and fiber wall thickness (FWT), fiber coarseness (FCS) and fiber width in radial (FRW) and tangential (FTW) direction. DEN served as a benchmark for evaluation of the Resistograph and Pilodyn for indirect wood density measurement. Earlywood, transition wood and latewood of each growth ring were defined as sections with densities ranging from 0-20%, 20-80% and 80-100% of the total density range within the annual ring (minimum to maximum), respectively. Since mean values weighted by the ring areas were shown to more accurately represent average wood properties [23], area-weighted values (AWV) for fiber and wood quality traits were calculated as:(3)where di is a value for annual ring i with an area αi [24].

Growth traits

Height (HGT) and diameter at breast height (DBH) were measured in summer 2011 and stem volume (VOL) was calculated according to [25] as a function of height and diameter.

Statistical analysis

The response variables were fitted into the following linear mixed model using statistical package ASReml 4 [26]:(4)where yijkl is the lth observation for an offspring of jth and kth parent growing in ith block, μ is the overall mean of a given variable, Bi is random effect of ith block, Gj and Gk are random general combining ability effects of the jth and kth parent, respectively, Sjk is random specific combining ability effect for the cross between parents j and k, and eijkl is random error term.

Individual-tree narrow-sense heritabilities () for each variable were estimated using variance components from the univariate analysis as(5)where , and are variances for additive genetic, phenotypic, general combining ability, specific combining ability, and residual components, respectively. Standard errors were obtained using Taylor series expansion [26]. Phenotypic and genetic correlations (rxy) were calculated as(6)where and are phenotypic or additive genetic variance components for traits x and y, respectively, and σxy is phenotypic or additive genetic covariance component between traits x and y estimated by fitting a bivariate mixed model (Eq 4) [26]. Significance of the correlation coefficients was examined using t-test(7)where t is a t-value with n-2 degrees of freedom, n is the number of pairs and rxy is the Pearson product-moment correlation. Furthermore, hierarchical cluster analysis based on dissimilarity matrix of additive genetic correlations between traits was performed using hclust function in R program [27] in order to construct an illustrative dendrogram. Genetic gain () for direct selection was calculated as(8)where i is selection intensity (1% = 2.665), is individual-tree narrow-sense heritability for trait x, and CVx is coefficient of variation for trait x calculated as phenotypic standard deviation divided by the mean. Correlated response (CRy) to selection for a target trait y was calculated as(9)where hx and hy are square roots of narrow sense heritabilities for selection trait x and target trait y, respectively, rxy is the genetic correlation between trait x and y, and CVy is the coefficient of variation for target trait y.

Results

Variation in the studied traits

Descriptive statistics for wood density estimated by Resistograph, Pilodyn and SilviScan as well as for other wood, fiber and growth traits are summarized in Table 1. SilviScan area-weighted wood density (DEN), used as a benchmark trait in our evaluation study, varied between 365 and 566 kg.m-3, with mean value of 448 kg.m-3. On average, wood contained a substantially higher proportion of earlywood (55%) than latewood (15%). The Resistograph unadjusted wood density (RES) ranged from 1242 to 2403 density units from side A, and from 1427 to 2492 density units from side B. Measurements from side A showed a slightly higher variation (10.7%) compared to side B (9.7%). Linear detrending of resistograms lowered wood density estimates, while debarking increased them since the lower-density bark was excluded (REST = 1644 < RESTB = 1848 < RES = 1857). Pilodyn penetration depth measured with bark (PIL) ranged from 15 to 31 mm with an average of 21.7 mm, and it decreased after bark removal (PILB) to only between 11 and 22 mm with an average of 16.4 mm. Variation of Pilodyn penetration depth was slightly higher after bark removal (11.4%). VOL exhibited the highest variation among the measured traits (34.5%), followed by MFA (23.8%).

[Figure omitted. See PDF.]

Table 1. List of variables and descriptive statistics–minimum, maximum, mean, standard deviation (SD), coefficient of variation (CV) and individual narrow-sense heritabilities () with standard errors in parentheses.

https://doi.org/10.1371/journal.pone.0204518.t001

Heritability

Individual-tree narrow-sense heritability of Resistograph wood density substantially increased after detrending of raw density profiles (from 0.35 to 0.44), reaching a similar value as that of the benchmark variable DEN (0.46) obtained from SilviScan (Table 1). Removing bark however did not bring any further improvement. Pilodyn density heritability PIL (0.32) was close to that attained by the Resistograph using raw data (0.35), but since cutting off bark prior to the measurements did not bring any improvement to the estimate either (PILB remained at 0.32), the Resistograph proved to be a superior tool in this regard. Fiber traits exhibited highest heritabilities, particularly FCS (0.57), FTW (0.55) and FWT (0.53), followed by LWD (0.51). Compared to wood and fiber traits, growth traits had lower heritabilities (0.22–0.35).

Phenotypic and genetic correlations of indirect wood density measurements with SilviScan wood density as a benchmark

Phenotypic correlations of unadjusted density RES with the benchmark DEN were moderate (0.53, 0.56 and 0.59 for side A, B and mean of the two sides, respectively) (Fig 3A, Table 2) but they considerably improved following detrending (REST; 0.63, 0.68 and 0.69, respectively) and even slightly more after subsequent debarking (RESTB; 0.66, 0.71, and 0.72, respectively; Fig 3B). Compared to phenotypic correlations, genetic correlations showed only minor differences between the two sides and their mean values. Genetic correlations of the benchmark DEN with mean RES, REST and RESTB increased from 0.89 through 0.93 up to 0.96, respectively. Phenotypic and genetic correlations between DEN and PILB (0.44 and 0.74, respectively) were higher than between DEN and PIL (only 0.38 and 0.59, respectively) (Fig 4, Table 2).

[Figure omitted. See PDF.]

Fig 3.

Relationship between SilviScan area-weighted mean density (DEN) and (A) unadjusted Resistograph mean density (RESA and RESB on the left, RES on the right) and (B) adjusted Resistograph mean density (RESTB_A and RESTB_B on the left, RESTB on the right).

https://doi.org/10.1371/journal.pone.0204518.g003

[Figure omitted. See PDF.]

Fig 4.

Relationship between SilviScan area-weighted mean density (DEN) and inverse value of Pilodyn penetration depth with bark (left) and without bark (right).

https://doi.org/10.1371/journal.pone.0204518.g004

[Figure omitted. See PDF.]

Table 2. Genetic and phenotypic correlations of SilviScan area-weighted mean density (DEN) with different Resistograph and Pilodyn wood density estimates (standard errors in parentheses).

https://doi.org/10.1371/journal.pone.0204518.t002

Phenotypic and genetic correlations of indirect wood density measurements with other wood, fiber and growth traits

Phenotypic and genetic correlations of the benchmark DEN and the indirect wood density variables (RES, REST, RESTB, PIL and PILB) with other wood, fiber and growth traits are summarized in Table 3. In most cases, the adjusted Resistograph wood density gave stronger correlations associated with lower standard errors compared to unadjusted measurements (). Strong positive genetic correlations between Resistograph wood density and the three wood density components (EWD, TWD and LWD) reflected those with benchmark DEN, although they were a little lower. Genetic correlations of the three components with Pilodyn measurements were moderate, but PILB had a stronger relationship than PIL. Phenotypic correlations of the components were strong with DEN, moderate with Resistograph measurements and weak with Pilodyn measurements. Correlations with FWT followed a similar pattern. The indirect wood density estimates, including DEN, had strong to moderate positive correlations with LWP but their correlations with EWP and TWP were negative, whereas genetic correlations were in many cases non-significant and phenotypic correlations were associated with high standard errors. DEN showed moderate negative genetic and weak negative phenotypic correlations with MFA. Both genetic and phenotypic correlations of MFA with Resistograph and Pilodyn measurements were very weak. Moderate genetic and weak phenotypic correlations between MOEs and Resistograph and Pilodyn measurements were lower than those with the benchmark DEN. Growth traits generally correlated negatively with DEN as well as with the indirect wood density estimates (Fig 5). Genetic correlations with DBH and VOL were moderate, while correlations with HGT were weak. Relationships among growth and wood quality traits based on additive genetic correlations are visualized in Fig 6.

[Figure omitted. See PDF.]

Fig 5. Relationships of stem volume (VOL) with SilviScan (DEN), Resistograph (RESTB) and Pilodyn (PIL) wood density.

https://doi.org/10.1371/journal.pone.0204518.g005

[Figure omitted. See PDF.]

Fig 6. Dendrogram of genetic relationships among growth and wood quality traits.

https://doi.org/10.1371/journal.pone.0204518.g006

[Figure omitted. See PDF.]

Table 3. Genetic and phenotypic correlations of indirect wood density estimates (SilviScan, Resistograph and Pilodyn) with growth, fiber and other wood quality traits (standard errors in parentheses).

https://doi.org/10.1371/journal.pone.0204518.t003

Genetic gain and response to selection

Genetic gain and correlated genetic response to indirect selection for target traits using selection intensity of 1% is presented in Table 4. Selection based on VOL and DBH resulted in positive gain for all growth traits VOL, DBH and HGT as well as for MFA (ca 20%, 9%, 4% and 5%, respectively) and in negative gain for DEN and MOEs (ca -3% and -6%, respectively). On the other hand, selection based on wood density traits (DEN, RESTB, PIL and PILB) generated negative gain for growth traits and MFA and positive gain for wood quality traits. Selection for DEN led to its moderate genetic gain (9%) and a relatively big correlated response for MOEs (18%), MFA (-12%) and VOL (-13%), while selection based on RESTB, PIL and PILB had a lower impact on MOEs (12%, 8% and 8%, respectively), MFA (-5%, -4% and -2%, respectively) and DEN (8%, 4% and 6%, respectively).

[Figure omitted. See PDF.]

Table 4. Genetic gain (in parentheses) and correlated genetic response (both in %) for selected traits with selection intensity of 1%.

https://doi.org/10.1371/journal.pone.0204518.t004

Discussion

Predictability of wood density using Resistograph and Pilodyn

The relevance of inclusion of wood quality traits (such as density, stiffness or strength) into breeding programs of different conifer species, especially when genetic correlations between growth and wood quality traits are unfavorable, has been intensively studied in a number of species (e.g. [4, 8, 12, 15, 22, 28]), particularly in radiata pine [29–31]. Breeding strategies aiming at overcoming such adverse correlations have been designed [32, 33] and economic weights for optimal combination of growth and stiffness have also been applied in real breeding selection [34]. However, the time and cost demanding assessment of wood quality on standing trees, along with a low market pressure, often prevents wood quality traits from being taken into account during selection [15]. Nevertheless, in spite of the lack of interest from the current wood processing industry, breeders should anticipate prospective needs already now, as long-term breeding of forest trees is not capable of flexible responding to rapid technological advances and market changes [22]. Wood density is generally considered to be the most important wood quality trait because it greatly influences wood properties and performance and thus significantly affects wood suitability for different end uses. Integration of wood density into a breeding program, however, requires non-destructive screening of a large number of standing trees.

Pilodyn penetrometer was successfully applied for wood density estimation in some studies (e.g. [15, 35]) whereas it failed in other studies (e.g. [6, 36]). Application of Pilodyn is quick and simple but it has a shallow penetration, low sensitivity and gets affected by thick bark [19]. In this study, Pilodyn was applied for Scots pine wood density estimation of standing trees both with and without bark. Phenotypic correlations of inverse Pilodyn penetration (PIL) with SilviScan wood density (DEN) were weak, while the genetic correlation was moderate (rG = 0.59) before and relatively strong (rG = 0.74) after bark removal. Correlations of a similar magnitude were found in black spruce [18] or loblolly pine [37]. Nevertheless, Pilodyn does not seem to be a suitable tool for wood density assessment in Scots pine compared to e.g. Norway spruce or Eucalyptus nitens, for which strong genetic correlations were estimated (rG = 0.96) [15, 35].

In contrast to Pilodyn, Resistograph micro-drill measures a whole stem profile with a high accuracy. Its main drawback is, however, a need for adjustment of drilling profiles that show an increasing trend as an effect of accumulated needle friction. Approaches for adjustment of Resistograph measurements vary among authors. For example, [38] and [39] applied centered moving averages and centered moving minima to smooth and detrend profiles, while [22] and [18] used only a few first centimeters of the profile and thus eliminated part of the profile exhibiting a steep increase.

In the present study, the whole Resistograph drilling profiles were detrended, assuming that their increasing trend was linear, and bark was excluded from both sides because bark can reach considerable thickness in Scots pine and thus can constitute a substantial portion of the drilling profile. A strong genetic correlation (0.89) between Resistograph and SilviScan densities was attained already with raw data, and an even stronger correlation (0.96) was reached after profiles’ adjustment. The effect of detrending & debarking was more pronounced at phenotypic correlations (0.59 → 0.72). The results indicate that an adjustment of raw drilling profiles is an important procedure to obtain more accurate wood density estimates and that the Resistograph is a reliable tool for non-destructive wood density assessment of standing Scots pine trees. Similar results were also reported e.g. for loblolly pine (rP = 0.75, rG = 0.92) [38] or maritime pine (rP = 0.77–0.79, rG = 0.96 – 0.98) [22].

Variation in Resistograph density numbers

Resistograph measurements taken in the southeast direction (denoted as A) showed a higher variation and a lower minimum value compared with measurements taken from the southwest direction (denoted as B). This could be caused by wood distortion and/or wood deterioration in some of the assessed trees, as the ten-millimeter-cores extracted for SilviScan analysis were taken in the southeast direction four years earlier. The Resistograph was applied in a reasonably close proximity of healed boring holes with the intention to obtain profiles comparable to those from SilviScan but, at the same time, some distance (ca 5–10 cm) was kept to avoid drilling through wood affected by the increment core extraction.

Heritability

Narrow-sense individual-tree heritability of Resistograph wood density varies between 0.14 and 0.64 [8, 18, 40]. In our study, heritability for adjusted Resistograph density was moderate (0.43), similar to that reported by [22] for maritime pine. Estimates of heritability for Pilodyn penetration range from very low to high. [41] reported individual narrow-sense heritabilities between 0.11 and 0.90 for Scots pine growing at different sites, while [37] between 0.06 and 0.46 for loblolly pine measured by different types of Pilodyn. A high heritability (> 0.6) was determined by [35] for Eucalyptus nitens. In this study, individual narrow-sense heritability of 0.32 was estimated for both Pilodyn penetration with bark and without bark. The results are similar to those reported for Norway spruce by [42] and [15].

Correlations

We have determined negative genetic correlations between wood density and growth traits (-0.47 and -0.44 for DEN vs. DBH and VOL, respectively), which is well in line with earlier studies in Scots pine [3, 24, 43]. Genetic and phenotypic correlations of DBH and VOL with indirect wood density measurements (RESTB, PIL, PILB) followed the same pattern as those with benchmark SilviScan. The only exceptions were phenotypic correlations of PIL, which were double compared to phenotypic correlations of PILB or DEN. We assume that trees with larger diameter have thicker bark and thus Pilodyn’s penetration is deeper as the pin has to penetrate through lower-density bark first before it reaches wood. That is the most plausible explanation of Pilodyn’s higher scores in larger trees, which, using the inverse of Pilodyn penetration depth, resulted in a stronger negative phenotypic correlation with DBH than for Pilodyn without bark (PILB). Nevertheless, such a big difference between PIL and PILB did not emerge in genetic correlations. The effect of Resistograph profiles’ adjustment was clearly visible in the magnitude of genetic correlations and their standard errors. Genetic correlations of DBH and VOL with RES were weak, whereas they were moderate with RESTB and very close to the benchmark values. A similar pattern between genetic correlations of growth traits with wood density estimated by Resistograph or Pilodyn and genetic correlations of growth traits with benchmark wood density were also published in other studies; either the correlations were moderately negative [12, 15, 18 for a mixed forest only] or near zero [38].

Conclusion

In order to find a reliable tool for rapid and non-destructive wood density assessment of standing Scots pine trees, we tested penetrometer Pilodyn and micro-drill Resistograph using SilviScan wood density as a benchmark. As Pilodyn’s measurement is limited to a few centimeters of outer wood, we applied Pilodyn with and without bark, so that the influence of bark could be evaluated. According to our results, it seems that removing bark leads to more accurate estimates; nevertheless, the narrow-sense heritability remained unchanged. Phenotypic correlations of stem diameter with Pilodyn penetrations suggest that conducting the measurements with bark can result in underestimating of wood density as the lower-density bark might constitute a substantial portion of the measured profile. Unlike Pilodyn, Resistograph considers the whole stem profile and provides a detailed scan of wood density variation within it. Resistograms that exhibit increasing trend should be adjusted in order to improve accuracy of wood density estimates. We proposed linear detrending followed by bark removal as the accuracy of Resistograph wood density estimates substantially increased after the detrending and even slightly more after subsequent debarking. We also compared Resistograph’s measurements taken in two mutually perpendicular directions and concluded that a single measurement would be sufficient if a proper distance from knots or visible stem damages was kept. In summary, having compared Pilodyn and Resistograph wood density estimates with the benchmark wood density of SilviScan, it is obvious that Resistograph is a more suitable tool for wood density assessment in Scots pine standing trees (even before drilling profiles’ adjustment) than Pilodyn.

Acknowledgments

The authors thank Skogforsk (The Forestry Research Institute of Sweden) for providing information about and access to the progeny trial and Zhou Hong, Anders Fries, Zhiqiang Chen, Henrik Hallingbäck and David Hall for help with field work.

Citation: Fundova I, Funda T, Wu HX (2018) Non-destructive wood density assessment of Scots pine (Pinus sylvestris L.) using Resistograph and Pilodyn. PLoS ONE 13(9): e0204518. https://doi.org/10.1371/journal.pone.0204518

References

1. Krakau U-K, Liesebach M, Aronen T, Lelu-Walter M-A, Schneck V. Scots Pine (Pinus sylvestris L.). In: Pâques LE, editor. Forest tree breeding in Europe—current state-of-the-art and perspectives. Vol 25. Managing forest ecosystems. Netherlands: Springer 2013. p. 267–323.

2. Wilhelmsson L, Andersson B. Breeding of Scots pine (Pinus sylvestris) and lodgepole pine (Pinus contorta ssp. latifolia). In: Lee SJ, editor. Corrected repring from: Progeny testing and breeding strategies, proceedings of the Nordic group of tree breeding. Edinburgh: Forestry Commission; 1993. p. 184.

3. Fries A. Genetic parameters, genetic gain and correlated responses in growth, fibre dimensions and wood density in a Scots pine breeding population. Annals of Forest Science. 2012; 69(7):783–94.

4. Hong Z, Fries A, Wu HX. High negative genetic correlations between growth traits and wood properties suggest incorporating multiple traits selection including economic weights for the future Scots pine breeding programs. Annals of Forest Science. 2014; 71(4):463–72.

5. Li L, Wu HX. Efficiency of early selection for rotation-aged growth and wood density traits in Pinus radiata. Canadian Journal of Forest Research. 2005; 35(8):2019–29.

6. Pot D, Chantre G, Rozenberg P, Rodrigues JC, Jones GL, Pereira H, et al. Genetic control of pulp and timber properties in maritime pine (Pinus pinaster Ait.). Annals of Forest Science. 2002; 59(5–6):563–75.

7. Zhang SY, Morgenstern EK. Genetic variation and inheritance of wood density in black spruce (Picea mariana) and its relationship with growth: Implications for tree breeding. Wood Science and Technology. 1995; 30(1):63–75.

8. El-Kassaby YA, Mansfield S, Isik F, Stoehr M. In situ wood quality assessment in Douglas-fir. Tree Genetics & Genomes. 2011; 7(3):553–61.

9. Armstrong JP, Skaar C, Dezeeuw C. The effect of specific gravity on several mechanical properties of some world woods. Wood Science and Technology. 1984; 18(2):137–46.

10. Hernández RE, Bustos C, Fortin Y, Beaulieu J. Wood machining properties of white spruce from plantation forests. Forest Products Journal. 2001; 51(6):82–8.

11. Glass SV, Zelinka SL. Moisture relations and physical properties of wood. Wood handbook: Wood as an engineering material: chapter 4. General Technical Report FPL-GTR-190. Madison, Wisconsin: U.S. Department of Agriculture, Forest Service, Forest Products Laboratory; 2010. p. 4.1–4.19.

12. Isik F, Li BL. Rapid assessment of wood density of live trees using the Resistograph for selection in tree improvement programs. Canadian Journal of Forest Research. 2003; 33(12):2426–35.

13. Evans R. Rapid measurement of the transverse dimensions of tracheids in radial wood sections from Pinus radiata. Holzforschung. 1994; 48(2):168–72.

14. Gao S, Wang X, Wiemann MC, Brashaw BK, Ross RJ, Wang L. A critical analysis of methods for rapid and nondestructive determination of wood density in standing trees. Annals of Forest Science. 2017; 74(2):1–13.

15. Chen Z-Q, Karlsson B, Lundqvist S-O, Gil MRG, Olsson L, Wu HX. Estimating solid wood properties using Pilodyn and acoustic velocity on standing trees of Norway spruce. Annals of Forest Science. 2015; 72(4):499–508.

16. Yanchuk AD, Kiss GK. Genetic variation in growth and wood specific gravity and its utility in the improvement of interior spruce in British Columbia. Silvae Genetica. 1993; 42(2–3):141–8.

17. Fukatsu E, Tamura A, Takahashi M, Fukuda Y, Nakada R, Kubota M, et al. Efficiency of the indirect selection and the evaluation of the genotype by environment interaction using Pilodyn for the genetic improvement of wood density in Cryptomeria japonica. Journal of Forest Research. 2011; 16(2):128–35.

18. Desponts M, Perron M, DeBlois J. Rapid assessment of wood traits for large-scale breeding selection in Picea mariana [Mill.] B.S.P. Annals of Forest Science. 2017; 74(3):1–12.

19. Cown D. Comparison of the pilodyn and torsiometer methods for the rapid assessment of wood density in living trees. New Zealand Journal of Forestry Science. 1978; 8(3):384–91.

20. Rinn F, Schweingruber FH, Schär E. Resistograph and X-ray density charts of wood. Comparative evaluation of drill resistance profiles and X-ray density charts of different wood species. Holzforschung. 1996; 50(4):303–11.

21. Rinn F. Basics of micro-resistance drilling for timber inspection. Holztechnologie. 2012; 53(3):24–9.

22. Bouffier L, Charlot C, Raffin A, Rozenberg P, Kremer A. Can wood density be efficiently selected at early stage in maritime pine (Pinus pinaster Ait.)? Annals of Forest Science. 2008; 65(1):1–8.

23. Gräns D, Hannrup B, Isik F, Lundqvist S-O, McKeand S. Genetic variation and relationships to growth traits for microfibril angle, wood density and modulus of elasticity in a Picea abies clonal trial in southern Sweden. Scandinavian Journal of Forest Research. 2009; 24(6):494–503.

24. Hannrup B, Ekberg I, Persson A. Genetic correlations among wood, growth capacity and stem traits in Pinus sylvestris. Scandinavian Journal of Forest Research. 2000; 15(2):161–70.

25. Brandel G. Volume functions for individual trees: Scots pine (Pinus sylvestris), Norway spruce (Picea abies) and birch (Betula pendula & Betula pubescens). Garpenberg: Swedish University of Agricultural Sciences. Department of Forest Yield Research. Report 26: 183 p. ISSN 0348-7636.; 1990.

26. Gilmour A, Gogel B, Cullis B, Welham S, Thompson R. ASReml user guide. Release 4.1. Structural specification. Hemel Hempstead: VSN International Ltd. 2015.

27. RCoreTeam. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2013.

28. Hayatgheibi H, Fries A, Kroon J, Wu HX. Genetic analysis of lodgepole pine (Pinus contorta) solid-wood quality traits. Canadian Journal of Forest Research. 2017; 47(10):1303–13.

29. Baltunis BS, Wu HX, Powell MB. Inheritance of density, microfibril angle, and modulus of elasticity in juvenile wood of Pinus radiata at two locations in Australia. Canadian Journal of Forest Research. 2007; 37(11):2164–74.

30. Gapare WJ, Baltunis BS, Ivković M, Wu HX. Genetic correlations among juvenile wood quality and growth traits and implications for selection strategy in Pinus radiata D. Don. Annals of Forest Science. 2009; 66(6):1–9.

31. Wu HX, Ivković M, Gapare WJ, Matheson AC, Baltunis BS. Breeding for wood quality and profit in Pinus radiata: a review of genetic parameter estimates and implications for breeding and deployment. New Zealand Journal of Forestry Science. 2008; 38(1):56–87.

32. Wu HX, Sanchez L. Effect of selection method on genetic correlation and gain in a two-trait selection scheme. Australian Forestry. 2011; 74(1):36–42.

33. Hallingbäck HR, Sánchez L, Wu HX. Single versus subdivided population strategies in breeding against an adverse genetic correlation. Tree Genetics & Genomes. 2014; 10(3):605–17.

34. Ivković M, Wu HX, McRae TA, Powell MB. Developing breeding objectives for radiata pine structural wood production. I. Bioeconomic model and economic weights. Canadian Journal of Forest Research. 2006; 36(11):2920–31.

35. Greaves BL, Borralho NMG, Raymond CA, Farrington A. Use of a Pilodyn for the indirect selection of basic density in Eucalyptus nitens. Canadian Journal of Forest Research. 1996; 26(9):1643–50.

36. Livingston AK, Cameron AD, Petty JA, Lee SL. Effect of growth rate on wood properties of genetically improved Sitka spruce. Forestry. 2004; 77(4):325–34.

37. Sprague JR, Talbert JT, Jett JB, Bryant RL. Utility of the Pilodyn in selection for mature wood specific gravity in loblolly pine. Forest Science. 1983; 29(4):696–701.

38. Eckard JT, Isik F, Bullock B, Li BL, Gumpertz M. Selection efficiency for solid wood traits in Pinus taeda using time-of-flight acoustic and micro-drill resistance methods. Forest Science. 2010; 56(3):233–41.

39. Isik F, Mora CR, Schimleck LR. Genetic variation in Pinus taeda wood properties predicted using non-destructive techniques. Annals of Forest Science. 2011; 68(2):283–93.

40. Gwaze D, Stevenson A. Genetic variation of wood density and its relationship with drill resistance in shortleaf pine. Southern Journal of Applied Forestry. 2008; 32(3):130–3.

41. Haapanen M, Velling P, Annala ML. Progeny trial estimates of genetic parameters for growth and quality traits in Scots pine. Silva Fennica. 1997; 31(1):3–12.

42. Hallingbäck HR, Jansson G, Hannrup B. Genetic correlations between spiral grain and growth and quality traits in Picea abies. Canadian Journal of Forest Research. 2010; 40(2):173–83.

43. Fries A, Ericsson T. Estimating genetic parameters for wood density of Scots pine (Pinus sylvestris L.). Silvae Genetica. 2006; 55(2):84–92.