Genetic Parameters of Stem and Wood Traits in Full-Sib Silver Birch Families

This study investigated heritability of stem and wood traits to improve Swedish silver birch (Betula pendula Roth.) through breeding. Birch is 12% of Sweden’s forest area but mainly used for low value pulp or firewood. This paper applied non-destructive test (NDT) methods, and estimated traits’ heritability (h2), to help breed birch for high value solid wood products. Two trials of 22 families were assessed at age 19 for stem diameter (DBH), stem straightness, rough brown bark height (BH), grain angle (GA), Pilodyn penetration depth (Pilo) and acoustic velocity (AV). X-ray densitometry was performed on a subsample of radial cores taken at 1.3 m from the ground to get an average benchmark density. The h2 values were moderate for GA (0.20 and 0.21) and Pilo (0.53 and 0.48) at the two sites, but the h2 values for AV were low (0.05 and 0.30). There were moderate genotypic correlations between BH and DBH (0.51–0.54). There were low genotypic and phenotypic correlations between NDT measurements and other traits so including NDT in birch breeding efforts should not inadvertently reduce size, stem or wood quality. The high genetic correlations between sites suggest that GA, Pilo and AV values were determined more by genotype than by environment.


Introduction
Over 12% of Sweden's productive forest volume is composed of birch (Betula spp.) [1]. The majority of harvested birch is used for pulp [2] or firewood production as it is considered to be low quality. The volume of high quality wood is presumably insufficient to sustainably supply a domestic wood products industry. Birch regularly occurs in forests managed for Norway spruce or Scots pine production [3] where there is limited interest in managing birch for high quality wood. These birches are scattered in pure or mixed softwood stands and naturally regenerated from unimproved material, especially during early succession stages, such as after clear-cutting or storm damage [3,4]. The current recommendations are to plant improved material and adopt specialized silvicultural practices to produce high-quality birch wood and increase the value of Sweden's birch stands.
The two main species of birch in Swedish forests are silver birch (B. pendula Roth.) and downy birch (B. pubescens Ehrh.), and are commonly treated as a single species for inventory and harvest [1]. Of these two species, downy birch is able to grow on compact or wet sites, while silver birch requires relatively fertile sites with appropriate levels of water for vigorous growth [5,6]. Silver birch has been found to have higher volume production [7], straighter form [8], and higher wood density [9]. Therefore, silver birch is of interest to the Swedish tree breeding program. The breeding program for silver birch started in the 1980s with the main intention to improve volume production, and then to improve stem heritability values for AV; from 0.16 to 0.74 [28] for E. nitens, 0.20 for E. pilularis [43] and 0.62 for E. dunnii [43].
The aim of this study was to investigate phenotypic and genetic variation in growth and wood properties in two clonal experiments at the mid-harvest age, established with clones from 22 full-sib silver birch families produced in the Swedish birch breeding program. Relationships between birch growth traits, external quality traits, and wood properties measured with NDT, were studied to find if an external trait were highly correlated with wood properties, and could then be used in the breeding program as selection criteria. Wood density was estimated by means of indirect and direct methods to determine if birch wood density can be measured indirectly by means of the Pilodyn. The hypothesis tested were that there is no relationship between average family diameter and wood density, i.e., selection for diameter is not decreasing biomass production, and that there is a strong relationship between Pilodyn measurements and wood densities determined from cores. Determining if NDT values are heritable, and that there is sufficient variation in values, was also a key objective in this study, as this relates to the hypothesis that wood quality traits can be included in the birch breeding programs.

Materials and Methods
The material used for this study was collected in two clonal trials with silver birch located in Hyssna (57 • 29 40 N, 12 • 33 55 E) and Brunsberg (61 • 28 45 N, 13 • 92 77 E), south to mid Sweden ( Figure 1). Hyssna is generally considered to be the superior site for growth as average monthly temperatures range from −4 to 21 • C and the average monthly rainfall varies from 39 to 80 mm. Brunsberg has a slightly shorter growing season, average monthly temperatures ranging from −8 to 22 • C, and the average monthly rainfall fluctuates between 19 and 82 mm [64]. Both trials were fenced before planting to eliminate browsing. Both trials were planted in a 2 × 2 m spacing in spring 1999 with one-year old material. The experimental design was a single-tree plot design where each tree is randomly assigned to a planting spot. Each tree's ID code, genotype and parent information was then associated with its location: tree number within row, row number within block, and block number. Not all genotypes were present in each block due to the incomplete blocking design. There was a buffer zone surrounding the trial. A total of 3534 and 2479 trees were planted in Hyssna and Brunsberg, respectively.
The planted material included clones from 22-full-sib families produced by crossing trees in the Swedish elite birch population. The elite population contains the 20 "best birch" clones for a deployment zone below 59.5 • latitude. The elite population was established to increase the genetic gain over four birch breeding populations delimited in southern Sweden. There were seven unique mother trees and ten unique father trees in the crosses, with no reciprocal crosses (see Appendix A Table A1). Any tree could be used as a mother and a father tree. One mother tree was used five times, two mothers were used four times, one mother was used three times, and three mothers were used twice. Three fathers were used three times, six fathers two times, and one father one time. There were between 28 and 35 planted genotypes per family/cross in Hyssna and between 8 and 34 in Brunsberg. There were between 2 and 9 individuals (ramets) per clone with an average of 4.6 individuals per clone. In addition, there were seedlings from two half-sib families and seedlings from three reference locations planted in the experiment.
The field measurements were performed on all stems at ages 6, 12 and 19 years for both sites. Diameter (DBH) was measured in the north-south direction at a height of 1.3 m from the stem base (breast height) for all trees. Height was measured at age 6 for all trees in both experiments. Hyssna was thinned between at 12 years and 19 years, reducing the number of stems by 30% during the first thinning. A summary of all measured traits are in Table 1, alongside units and trait descriptions. Stem properties of all remaining trees were assessed in May 2018 when stems were 19 years old from the date of planting. Straightness was ranked on a scale from one to nine, where a score of nine indicated the best quality. This method has been a standard assessment tool for straightness in the birch breeding program [10,53,65]. The assessment of straightness considered any deviation from a cylinder in the bottom 8 m stem section, with a score of 5 considered by the assessor to be a stem representative of the site. A total of 1123 stems at Hyssna and 1410 stems at Brunsberg were ranked for straightness. The height of brown rough bark was measured from the stump base up to the height at which smooth silver bark dominated. Any defects or damages to the stem were recorded, and included features such as scars, breakages and stem lean. In addition, ramicorns (also referred to as epicornic branches) and multileaders were counted for each stem. A subsample of trees was selected from the 22 control pollinated families, which had more than 20 individuals at both sites. Trees with severe defects, a DBH below 60 mm, or with a rough bark height above 140 cm, were excluded from the subsample. The 22 families then had 20 individuals selected randomly. The identity of the trees was tracked. If the tree has been damaged since the original assessment, it was not sampled, and then the next suitable stem was sampled if available. In this sense the sampling was not random, and some families had less than 20 trees sampled (at least 15 at both sites).
The subsample of stems was measured using NDT methods between August and September 2018. The grain angle, Pilodyn penetration depth, and AV along the stems were measured ( Table 1). The Pilodyn 6J Forest (PROCEQ, Zurich, Switzerland) densitometer was applied twice to the north and south stem faces without bark removal [27] at breast height. The second Pilodyn measurement was taken a few centimeters to the left of the initial measurement, and both were averaged to get north and south values. A wedge grain angle gauge (grain angle meter) was used at breast height to measure grain angle on the north and south sides [52,53]. These values were averaged to get a single grain angle value per stem. Stem AV in the longitudinal direction was measured using the Fakopp TreeSonic Timer (FAKOPP Enterprise, Agfalva, Hungary) and converted from ToF (µs m −1 ) to AV (m s −1 ) [66]. The TreeSonic timer was applied to the southern stem face between 0.8 and 1.8 m, so that the 1 m length was centered on breast height. If there were any issues or uncertainties for any of the NDT tools, measurements were repeated and sometimes shifted a few centimeters to avoid defects and reduce operator error.
Average family diameters were used to select twelve families to take cores from; four families with the highest average, four with lowest average, and the four families with mid-range diameters. For each family, five trees per site were randomly selected. A 5 mm diameter core was taken at breast height from the south through to the north face of these trees. In total, 58 cores were collected from Hyssna and 72 cores from Brunsberg, and prepared using the method outlined by Hallingbäck et al. [67].
The Itrax scanning densitometer (Cox Analytical Systems, Mölndal, Sweden) was used to create x-ray images of all the cores as described in Bergsten et al. [68]. Later, the images were analysed using the computer program LIGNOVISION™ (Rinntech, Heidelberg, Germany). The program provided average light intensity (LI) values for a user defined area, which included a single annual ring. LI was then converted to wood density by using a function that included temperature, relative humidity, LI values for a calibration sample, and sample thickness [67]. Average wood density (DENS Itrax ) was used for the remaining analysis, and is the average density value of all rings identified in LIGNOVISION™ for a given stem.
Variance and covariance components for the genetic analysis were estimated separately for each site using ASReml 4.0 [69] with a univariate linear mixed model in the following equation: where y is a vector of observed values, b is a vector of fixed effects (i.e., blocks, plots), α is a vector of random additive genetic effects, c is a vector of random effects of a clone within a full-sib family, f is a vector of random effects of a full-sib family and ε is a vector of random residuals. X, Z α , Z c and Z f , are design matrices linking phenotypic observations to fixed and random effects. The first-order autocorrelation for parameters was tested in the model and kept if the model converged and the log likelihood improved.
Estimates of the individual tree narrow-sense heritability (h 2 ) were calculated for each trait from the variance components obtained in the univariate analysis. Standard errors were estimated with the Taylor series expansion method.
where for any trait i: h 2 i is the narrow-sense heritability,σ 2 A is the additive genetic variance, σ 2 P is the total phenotypic variance (σ 2 A +σ 2 c +σ 2 f ), andσ 2 ε is the total error variance. This value is often used to indicate the proportion of variation in a trait that cannot be attributed to site or phenotypic factors.
Phenotypic coefficients of variation (PCV) and genotypic coefficients of variation (GCV) were calculated for each trait using the following formulas: where x is the general mean, σ 2 p is the phenotypic variance and σ 2 A is the additive genetic variance.
To estimate genetic correlations, the trait variance components and covariance components between pairs of traits were estimated using bivariate analysis. The bivariate model was an expansion of the univariate model to include a vector of values for two traits and matrices corresponding to fixed effect design and random effect design for each trait, respectively.
Variances and covariances of the random effects were estimated as: where A is the numerator relationship matrix, G is the additive genetic (co)variance matrix, C is the clone (co)variance matrix, F is the family co(variance) matrix, R is the residual (co)variance matrix, and I is an identity matrix. The n × n symmetric numerator relationship matrix A in was structured as: where p ij is the additive pedigree-based relationship between any pair of tree individuals i and j for all n individuals included in the analysis. The 2 × 2 order symmetric matrices G, C, F and R, are in turn structured as: are additive genetic, clone within family, fullsib family and residual variances for trait 1 and 2, respectively. The σ a12 , σ c12 , σ f 12 , and σ ε12 , are additive genetic, clonal within family, full-sib family and residual covariances between trait 1 and 2.
Phenotypic (cor P ) and genotypic (cor G ) correlations between traits and age-age genotypic correlations (R (x, y) ) were calculated as: whereσ 2 (1) andσ 2 (2) are the estimated phenotypic or genotypic variances for traits 1 and 2 or the same trait variances at two different ages, respectively, andĉov (1, 2) is the estimated phenotypic or genotypic covariance (σ a12 ) between traits 1 and 2 or between two different ages. These correlations indicate how similar two traits are, so that as one trait changes the other is likely to also change if the phenotype or genotype are linked. Values were recorded with a negative symbol if the correlation was negative.
The between site genetic correlations were estimated by compound variances for the linear mixed model and heterogeneous genetic variance within each site. The model used covariance parameterization [70] and was of the same general form as Equation (5) with the addition of random genetic site effects Z s s; where s is a vector of site random effects, and Z s is a design matrix linking the phenotypic observation to fixed and random effects.
A between-site genetic correlation of one is expected if the genetic influences on the trait are identical between sites, whereas a value of zero implies that the genetic effects on the trait are independent of the site. These are consistently referred to as genetic correlations to avoid confusion with the above genotypic correlations (cor G ).
Graphs were produced using the statistics software R (version 3.6.3, R Core Team, Vienna, Austria) for Pilodyn north against Pilodyn south measurements, Itrax density against stem diameter, and Itrax density against Pilodyn penetration depth. R was also used to plot linear trend lines between the variables. Correlations between variables were calculated using the cor() function and the Pearson's correlation method, which assumes a linear relationship between variables and indicates the strength and direction of this relationship (from −1 to +1).

Descriptive Statistics, Heritability and Variability of Traits
The average height of trees measured at age six was 71 cm for Hyssna and 47 cm for Brunsberg, with respective average diameters of 56 and 36 mm (Table 2). At 19 years old the diameters in Hyssna and Brunsberg were 110 and 98 mm, respectively. These size traits (height and diameters) had heritability estimates around 0.2 for Hyssna and below 0.1 for Brunsberg. PCVs were higher for Brunsberg (maximum 42%) than for Hyssna (maximum 27%), yet Hyssna had higher GCVs (maximum 13.5%) than Brunsberg (maximum 7.8%). Height at age 6 had the lowest coefficients of variation for any size trait, with PCV values of 15% at Hyssna and 28% at Brunsberg, and a GCV of 6.3% at both sites. Wood quality trait (GA, Pilo, AV and DENS Itrax ) average values were poorer at Hyssna than Brunsberg: GA was −1.8 • and −1.3 • , Pilo was 17.6 and 17.2 mm, AV was 4184 and 4237 m s −1 , and DENS Itrax was 431 and 444 kg m −3 (for Hyssna and Brunsberg, respectively). GA had PCVs above 60% and GCVs above 30% at both sites, and heritabilities were 0.20 at Hyssna and 0.21 at Brunsberg. Pilo had a PCV of 7.7% and GCV of 5.8% at Hyssna, and a PCV of 10.2% and GCV of 7.1% at Brunsberg. Heritability for Pilo was around 0.50 at both sites. Heritability for AV was non-existant at Hyssna (0.05) but moderate at Brunsberg (0.30), which also had PCVs around 5% and GCVs below 2%, the lowest of all wood quality traits. DENS Itrax had too few measurements to calculate genetic parameters.

Phenotypic and Genotypic Correlations between Traits
Phenotypic correlations (cor P ) were above 0.83 between diameters (DBH 6 , DBH 12 , and DBH 19 ), and above 0.7 between H 6 and diameters (Table 3). Genotypic correlations (cor G ) were above 0.9 between diameters (DBH 6 , DBH 12 , and DBH 19 ), and above 0.55 between H 6 and diameters (Table 3). Diameters had high cor P with BH (above 0.5), and moderate cor G with BH (0.23-0.5). Other cor P values were low between size traits, bark height and wood quality traits, with most cor P values close to ±0.1, and the largest cor P (−0.27) between AV and H 6 at Brunsberg.
The cor P between BH and wood quality traits were low, ranging from 0.04 for Pilo at Brunsberg to 0.25 for AV at Hyssna. BH had slightly higher cor G values with Pilo (−0.17 and −0.36) and AV (−0.11 and 0.29), for Hyssna and Brunsberg, respectively.
Wood quality traits (GA, Pilo, and AV) generally had low cor P (−0.03 to 0.16) and cor G (around 0.1) with other wood quality traits. There was a moderate (−0.46) cor G between GA and AV at Hyssna, but a low cor G value for Brunsberg (−0.07). GA had low to high (0.24 to 0.99) cor G with size traits.

Correlations between Sites and Measurements
The between site correlation for H6 was 0.76 the highest for any size trait ( Table 4). The same correlation for DBH6 was 0.44, DBH12 was 0.47, and DBH19 0.35. The between-site correlations were higher for stem form traits (above 0.6) than for size traits (Table 4). Wood quality traits had high between-site correlations; 0.96 for GA, 0.82 for Pilo, and 1.00 for AV. The average Pilodyn readings for the north and south stem faces had a correlation of 0.809 for Hyssna and 0.823 for Brunsberg ( Figure 2). Given the strong correlation, an average of the north and south Pilodyn readings for each stem (Pilo) is used elsewhere in this paper. DENS Itrax was moderately correlated with DBH19 ( Figure 3a

Discussion
Since the 1980s, the Swedish breeding program has managed to increase volume production per area (m 3 ha −1 ) by 15% and improved overall stem quality, compared to unimproved silver birch regeneration [13]. Unfortunately, the annual birch seedling production is less than 1% of the Norway spruce or Scots pine seedling production [71], reflecting the still limited interest in planting improved birch. This gap can be reduced if opportunities and markets are developed, and product education occurs, to encourage the use of birch wood [72]. Birch is believed to have lower production than the planted spruce and pine [73,74] and be too low quality for solid products, so any activity to improve growth rates, form or perception, could increase the amount of birch planted. As the use of broadleaved species is increasingly promoted in Sweden to meet goals for biodiversity [75], plantation management will need to focus on maximizing production value from broadleaf trees, both in mixtures and pure stands. For birch, this will include silvicultural choices such as improved genetic material and establishment methods.
It is important to test planting material across sites (GxE correlation) located in different conditions [76]. In our study, there were between site differences indicating that variables not related to genetics might influence the performance of genotypes. Phenotypic coefficients of variation were higher at Brunsberg than Hyssna. Brunsberg stem form traits and wood quality traits had higher genotypic coefficients of variation, while the reverse was true for size traits; diameters and height had lower genotypic variation at Hyssna. Higher genotypic variation is desirable in breeding programs as it indicates these traits have potential for improvement, and higher variability can make selection progress rapidly [77]. There were also high genetic correlations between sites for height, stem form, and wood quality traits, while diameters had moderate (below 0.5) between site genetic correlations. High between-site genetic correlations suggest that the best families will perform well at both sites as family performance rankings were stable.
The heritabilities for size traits reported here were slightly lower than in other genetics studies of silver birch in Sweden. In the study by Stener and Jansson [78], average heritibilities for diameter (0.29) and height (0.32) of open-pollinated silver birch trees at age 7 were similar to those in this study (0.04-0.28) and varied between sites. Gailis et al. [79] studied open-pollinated silver birch in Latvia at ages 10 and 14 years [79], and had heritibilities above 0.12 for diameter and 0.30 for height. This paper's heritabilities were much lower than those found by Skrøppa and Solvin [80] at age 16 for height (0.38) and diameter (0.81) for full and half-sib families in Norway. The genotypes used by Skroppa and Solvin were produced from a 9 × 9 diallel mating design [80], which may explain the relatively large differences in heritabilities compared to our study. Here, the moderately high heritability, and reasonable variability, indicate that selection for size traits is possible in the studied population.
Of the stem form traits, straightness has been more frequently assessed [78,79] than any bark metric [81]. For these other studies of silver birch [78,79], heritability for straightness ranged from 0.02 to 0.45 and was generally lower than for size traits. The heritibilities for straightness at Hyssna and Brunsberg, were above 0.2 in our study and would fit within this range of values. Viherä-Aarnio and Velling [81] measured bark thickness of 22-year-old silver birch in a provenance trial. They found that seed origin latitude was related to bark thickness, yet, unfortunately, did not provide heritabilities. In our study, heritability for bark height was very high at both Hyssna and Brunsberg, but genetic variation was very low at Hyssna. The results of this study suggest stem form traits can be improved through breeding. Stem straightness has previously been considered in the breeding program for birch [10,78]. The high genotypic variation and heritability of rough bark height makes this trait interesting and relatively easy to include in breeding, if there is enough commercial interest.
Grain angle heritability using the wedge grain angle meter has been commonly assessed for conifer species [31,32,52,82], yet it is difficult to find any hardwood species genetics trials which considered grain angle [83]. For 9-year-old open-pollinated family trials of Eucalyptus in Australia, E. dunnii had a heritability for grain angle of 1.0 [83] compared to 0.0 for E. pilularis [43], which is a large range of heritability estimates between species. The heritability for silver birch wood density in an 11-year-old Swedish clonal trial [10] was 0.73, which was higher than for Pilo in this study. Liang et al., studied 60 halfsib B. platyphylla families in China [84] and found heritability estimates of 0.30 and 0.37 for air-dry and oven-dry wood density, respectively. The same paper reported heritability values for flexural strength (0.37) and elastic modulus (0.15), which could be comparable to this study's AV measurements but included the influence of wood density on stiffness [37]. It is common for wood quality traits to have higher heritability estimates than size traits [85] as size traits tend to be under strong environmental influence.
The NDT tools used in this study to indirectly assess wood quality traits had fairly high heritabilities. Pilo had high heritability estimates at both sites, but Pilo had less variation than GA, which had low-moderate (0.2) heritability values. AV heritability was moderate at Brunsberg but extremely low at Hyssna, and had the lowest genetic coefficients of variation for any trait. Therefore, the results of this study should be complemented with more observations of B. pendula solid wood quality traits to see if these results are species or population specific, and to confirm that they are not unusual. One benefit of this work, is that the sampled trees are still growing, and could be assessed again later to see how well early performance relates to log mechanical properties, or how relationships vary with age [31,86]. The assessment of wood mechanical properties with NDT tools was also rapid, so the methodology should be easy for other researchers to apply.
The phenotypic correlations between bark height and wood quality traits were moderate, and genotypic correlations were slightly stronger. The rough bark height could affect the NDT measurements despite stems with rough bark at 1.3 m being excluded from the sample. It would be worth exploring these relationships further as rough bark height can be rapidly measured. Viherä-Aarnio and Velling [81] related bark thickness to seed origin latitude in a 22-year-old silver birch progeny trial in Finland. Phenotypic and genotypic correlations between wood quality traits were low to moderate, and varied by site. For wood quality traits this means that a single value is not representative, and that selection for one trait will not impact another.
Pilodyn's north and south values were highly correlated, and a single value should be sufficient for later assessment in silver birch field trails. Practical issues meant that taking two measurements per stem face was helpful, especially as the tool sometimes slipped or misfired. A second value can help to identify errors and provide more accurate estimates. Hyssna had higher growth, and lower Pilo and AV values, in agreement with studies of birch wood properties where faster growth rates reduced density [19,20] and strength [18]. However, the plots of diameter against density showed a poor relationship and a low correlation (phenotypic correlations were also low) meaning this difference is likely due to site conditions. In agreement with other studies, growth rate seemingly had no effect on wood density [87,88], but did vary by site. Itrax density and Pilodyn penetration were moderately correlated, although this relationship could be improved with more observations. The level of accuracy provided by Pilodyn has been sufficient for genetics trials comparing families, but not for comparing individuals within families [40]. Selection for the best trees to clone may be improved by using detailed density information, as provided by core scanning or destructive sampling.
The average DENS Itrax values were 431 kg m −3 at Hyssna and 444 kg m −3 at Brunsberg. Some comparable studies of juvenile silver birch [10,19,81] had average density values of: 478 kg m −3 for 22-year-old Finnish provenance trials [81], 455 kg m −3 for Latvian plantation-grown birch for ages 11 to 15 years [19], and of 429 kg m −3 for 11-year-old Swedish birch clones [10]. The density values obtained in this report are similar to other studies, yet there is variation within families, and often even within a stem [20,89], which is ignored when average values are used.

Conclusions
Pilodyn penetration was under strong genetic control, and grain angle had moderate heritabilities with high genotypic variation. This indicates that there is the potential to improve grain angle through breeding. However, low genotypic variation for Pilodyn measurements indicate that it will be harder to improve Pilodyn penetration through breeding, and consequently harder to influence wood density. Genotypic and phenotypic correlations between growth traits measured at different ages were high, indicating that early selection is possible for size traits. Low genotypic and phenotypic correlations between NDT measurements and growth traits indicate that wood quality traits cannot be improved by indirect selection. The high genetic correlations between sites suggest that grain angle, Pilodyn penetration and acoustic velocity values are determined by genotype rather than by environment. The results of the study indicate that NDT tools can be applied to silver birch genetics programs, yet more studies are needed to confirm their usefulness. The measurements for grain angle seem to have the greatest potential for improvement, while Pilodyn penetration also shows promise for inclusion. Pilodyn penetration was moderately correlated with wood density. The measurement was relatively quick and can be applied to a large number of individuals. No unfavorable correlations were observed between wood density and diameter, as are often observed in coniferous species. A broader study is needed to identify and quantify the effects of variables which could influence NDT measurements in birch genetics trials. Further work should also explore the possible relationships between bark height and size/wood quality traits. High between site genetic correlations support the use of NDT tools for comparing traits between families, but more accurate or destructive methods may be required to select individual trees. Funding: This work was financed by the research program FRAS-The Future Silviculture in Southern Sweden. The project is a collaboration between the Swedish University of Agriculture Science (SLU), Linnaeus University and Skogforsk (the Forestry Research Institute of Sweden). This research was also funded by the Kungl. Skogs-och Lantbruksakademien, grant number GFS2018-0073.

Data Availability Statement:
The measurements used in this report are also being analyzed separately for a study of the effect of stem form on NDT properties and measurements. Data can be made available to interested parties after this report is produced.

Conflicts of Interest:
The authors declare no conflict of interest.