1. Introduction
The future prosperity of the Canadian forest sector is increasingly dependent on its ability to embrace value-based management given that the economic viability of the traditional volumetric yield maximization proposition is becoming more challenging. Principally, this is due to a combination of factors which include increasing global competiveness, accessibility of economically-viable fibre sources, and evolving market demands for increased end-product diversity and value (e.g., [
1,
2,
3]). Operationally, this transitional shift has generated a renewed focus on deploying intensive silvicultural-based crop plans that can result in improvements in wood quality and the production of a more diverse stream of end-products throughout the rotation [
4]. Such a representative crop plan would consist of species, genotype and initial spacing control via plantation establishment of genetically-improved stock on well-prepared (scarified) sites followed by early vegetation management treatments and subsequent maintenance of optimal site occupancy levels through density management treatments (e.g., precommercial and commercial thinning; [
5]). Transitioning to a value-based management paradigm also requires improved operational intelligence for decision-making, particularly in relation to segregation and merchandising efficiency within the upstream portion of the forest products supply chain.
The quality and associated economic value of manufactured wood-based end-products derived from the merchantable stem portion of a harvested tree, such as pulp and paper products, dimensional lumber, engineered wood composites and utility poles, are largely dependent on the characteristics of internal fibre attributes (
Table 1). The ability to estimate these internal attributes and by extension classify standing trees according to their end-product potential before harvest, could provide the prerequisite information for optimal segregation decision-making (e.g., directing the trees to the most appropriate conversion facility upon harvest (pole yards, pulp and paper mills or saw mills)).
Relatively recently, a suite of innovative non-destructive operational survey tools have been developed for estimating the end-product potential of standing trees. These tools include: (1) time-of-flight acoustic velocity instruments for indirectly providing an estimate of wood stiffness such as the Director ST300 developed by Fibre-gen Inc. of Christchurch, New Zealand (e.g., [
8]) and the TreeSonic microsecond timer developed by Fakopp Enterprise, Ágfalva, Hungary (e.g., [
9]); and (2) micro-drill resistance and impact tools for indirectly providing a non-destructive estimate of wood density such as the Resistograph developed by Instrumenta Mechanik Labor GmbH of Wiesloch, Germany, and the Pilodyn developed by PROCEQ of Zurich, Switzerland, respectively (see [
10] for a comparative review).
Traditionally, one of the more important attributes associated with lumber quality and associated value as evident from its use in machine stress lumber grading [
11], is the degree of bending stiffness as quantified by the static modulus of elasticity. This attribute has been historically determined through destructive, laborious and expensive bending stress tests of the extracted end-product (e.g., dimensional lumber). However, its dynamic analogue (dynamic modulus of elasticity) has been shown to be a useful surrogate measure of static wood stiffness that can be non-destructively determined through its relationship to acoustic velocity and wood density (e.g., [
8,
12,
13]). The dynamic and static modulus of elasticity estimates have been shown to be highly correlated with each other (correlation coefficients >0.95) when compared within the same solid wood product (lumber), with the dynamic estimate being approximately 10% greater than that of the static value for a given wood sample (sensu [
14]). Furthermore, based on comprehensive reviews, the acoustic-based approach has been shown to be of consequential utility in the non-destructive estimation of wood stiffness of both standing trees and harvested logs (e.g., [
15]).
Conceptually, the velocity of a dilatational stress wave arising from a mechanically-induced impact that propagates through a standing tree, is related to dynamic modulus of elasticity according to Equation (1) (sensu [
16,
17]):
where
me is the dynamic modulus of elasticity or MOE
dyn (GPa),
P is a species- or sample-specific transverse/axial stain ratio (Poisson ratio) estimate that is commonly treated as an unknown constant when parameterizing the relationship,
wd(g) is a species- or sample-specific green wood density (kg/m
3) estimate, and
vd is the speed of the mechanically-induced dilatational stress wave (km/s) that propagates between a lower and upper probe positioned approximately at stem heights of 0.3 and 1.5 m, respectively. Deployments of simplified representations of this relationship in which acoustic velocity is singularly assumed to be an indirect measure of wood stiffness have been used as a surrogate response metric for evaluating thinning effects on wood quality (e.g., in Loblolly pine (
Pinus taeda L.; [
18]) and Douglas-fir (
Pseudotsuga menziesii (Mirb.) Franco; [
19]). In other simplified applications, a universal constant has been used in place of the species- or sample-specific
wd(g) estimate (e.g., [
8,
20]).
Operationally, in-forest segregation decisions based on acoustic-based stiffness estimates may result in optimal wood allocations and accompanying increases in economic profitability [
21]. However, empirical evaluations of the functional expression given by Equation (1) in which the dynamic modulus of elasticity and wood density are actually measured on clear xylem samples extracted from standing softwood trees, have shown considerable variation among investigations. Particularly, in terms of statistical significance, proportion of variability explained and precision of the resultant estimates (e.g., as reviewed by Hong ([
8]) and Mora ([
9])). Differences among studies in relation to species, locale, silvicultural treatment histories, instrumentation, environmental conditions at the time of sampling, and analytical approaches, have negated a definitive and conclusive determination of the overall merits of the acoustic approach to wood quality estimation. The degree of between-study variability suggests that acoustic relationships may be unique to the experimental approach and sample population utilized, and hence may be best evaluated on a species-, locale- and analytical- specific basis.
Advances in attribute determination via the introduction of the Silviscan system [
22,
23] have provided a means to evaluate the acoustic approach without resorting to costly and time-consuming methods. Briefly, the Silviscan system is an automated fibre analysis system originally developed by CSIRO’s (Commonwealth Scientific and Industrial Research Organisation) Forestry and Forest Products Division, in Australia. Based on a semi-empirical analytical approach, the system was designed to provide rapid and cost-effective estimates of wood quality attributes related to end-use performance [
24]. Utilizing the Silviscan-based modulus of elasticity estimate and the oven-dry wood density estimate (
wd) as a surrogate for
wd(g), Hong ([
8]) established a statistically significant (
p ≤ 0.05) acoustic relationship for 778 standing Scot pine (
Pinus sylvestris L.) trees growing in northern Sweden. In that study, 67% of the variation in the dynamic modulus of elasticity was explained by a simple linear regression model specification. Similarly, using Silviscan-based estimates of
me and
wd in conjunction with a simple linear regression model, Newton ([
25]) established a statistically significant (
p ≤ 0.05) acoustic relationship for 54 standing jack pine (
Pinus banksiana Lamb.) trees growing in the central portion of the Canadian Boreal Forest Region [
26]. Seventy-one percent of the variation in
me was explained by the parameterized model in that study.
In addition to the wood stiffness parameter,
me, wood density, microfibril angle, tracheid dimensions (radial and tangential diameters and wall thickness), fibre coarseness and specific surface area, are important attributes underlying the type, quality and value of end-products produced (
Table 1). Statistically, Silviscan-based
me estimates have been shown to be correlated with these secondary attributes as empirically exemplified in
Table 2 for three representative boreal softwood species: black spruce (
Picea mariana (Mill) B.S.P.), jack pine and red pine (
Pinus resinosa Ait.). As shown, Silviscan-based estimates of wood density, microfibril angle, radial and tangential tracheid diameters, tracheid wall thickness, fibre coarseness and specific surface area derived from a relatively large number of cross-sectional disk (xylem) samples, were mostly significantly (
p ≤ 0.05) correlated with the dynamic modulus of elasticity, across all 3 species. The only exceptions to these generalized correlative inferences were the relationships for tangential diameter for the pines and fibre coarseness for the spruce. Among the significant relationships, correlation coefficients ranged from a high of −0.898 for microfibril angle (black spruce) to a mean low of −0.304 for tangential tracheid diameter (black spruce).
Deploying these correlative relationships has enabled the formulation of a more encompassing acoustic-based inferential framework than that solely described by Equation (1) (sensu [
28]). Assuming
P is an unknown constant and utilizing the oven-dry Silviscan-based wood density estimate (
wd) as a surrogate for
wd(g) and the empirical-based correlative associations between
and microfibril angle (MFA denoted
in this study; °), tracheid wall thickness (
; µm), radial and tangential tracheid diameters (
(µm) and
(µm), respectively), fibre coarseness (
; µg/m) and specific surface area (
; m
2/kg), results in an expanded set of attribute-acoustic velocity associations. Specifically, these include the following:
- (1)
given that is inversely proportional to yields ;
- (2)
given that tracheid wall thickness (; µm) is directly proportional to , and radial (; µm) and tangential (; µm) tracheid diameters are inversely proportional to , yields the relationships and , respectively; and
- (3)
given that fibre coarseness (; µg/m) is directly proportional to , and specific surface area (; m2/kg) is inversely proportional to , yields the relationships and , respectively.
This expanded analytical structure inclusive of an assessment of the predictive ability of the individual relationships has yet to be completed for commercially-important softwood species in Canada. Consequently, the objective of this study was to investigate the nature, strength and predictability of the relationships between acoustic velocity and the Silviscan-based estimates of the dynamic modulus of elasticity, wood density, microfibril angle, tracheid wall thickness, radial and tangential tracheid diameters, fibre coarseness and specific surface area, for standing red pine trees. Furthermore, given that an estimate of wood density is required for the operational deployment of acoustic velocity analytical framework, the relationship between the micro-drill resistance (% amplitude as determined via the Resistograph [
29]) and wood density, was also investigated.
4. Discussion
End-product type and associated quality of an individual softwood tree is a function of its external morphological characteristics (e.g., stem diameter, height, sweep and taper, and the number and size of biotic and abiotic branches), and internal anatomical characteristics of the xylem tissue (e.g., modulus of elasticity, density, microfibril angle, tracheid wall thickness, radial and tangential tracheid diameters, fibre coarseness and specific surface area). Red pine produces a wide array of economically-important end-products which includes appearance-based boards used for interior flooring, exterior decking and wall panelling, dimensional lumber for residential home construction, utility poles used in building electrical transmission grids, veneer logs used in furniture manufacturing, and raw fibre for pulp for paper production and engineered wood composites [
37]. Consequently, estimating end-product potential before harvest could provide the prerequisite knowledge for increasing segregation and merchandizing efficiency within the upstream portion of the forest products supply chain.
More specifically, based on an examination of the relationship between acoustic velocity and a suite of commercially-relevant red pine fibre attributes, the results of this study indicated that 5 of the 8 attributes studied, could be non-destructively estimated from time-of-flight acoustic velocity measurements. Contrary to expectation, the results for the relationship where microfibril angle is expressed as a function of density-weighted acoustic velocity revealed no graphical or statistical support for such a relationship. Microfibril angle and the dynamic modulus of elasticity are inversely proportional as empirically evident by the significance and multitude of the correlation between these attributes for the 54 red pine trees analyzed in this study (i.e.,
r = −0.81 (
p ≤ 0.05);
Table 2). Furthermore, other investigations have reported significant relationships between microfibril angle and acoustic velocity (e.g., [
17,
38]). Thus the results reported here for red pine should be considered tentative and suggest that additional research should be initiated in order to arrive at a more conclusive determination of the microfibril angle - acoustic velocity relationship for this species. (e.g., determining and accounting for potential covariates that could be influencing the relationship).
Collectively, based on a set of goodness-of-fit, lack-of-fit and predictive ability criteria, the results of this study indicated that viable relationships could be obtained for
me,
wd,
wt,
co and
sa. Specifically, based on their statistical significance (
p ≤ 0.05;
Table 4), proportion of variability explained (40%, 14%, 45%, 27% and 43% of the variation in
me,
wd,
wt,
co and
sa explained, respectively;
Table 4), and predictive precision (e.g., 95% of all future errors would be within 12%, 8%, 7%, 8% and 6% of the true value of
me,
wd,
wt,
co and
sa, respectively;
Table 5)). Furthermore, given that a wood density estimate is required for deploying the
me,
wt,
co and
sa relationships, two non-destructive approaches for estimating
wd were also evaluated (acoustic and micro-drill resistance measures). Results from these analyses indicated that both approaches could provide unbiased wood density estimates at moderate levels of precision (e.g., ±8%;
Table 5 and
Table 7). Combining the acoustic-based density estimates with the parameterized functions revealed that
me,
wt,
co and
sa could be unbiasedly predicted at moderate levels of precision at the individual tree level and at relatively high levels of precision for stand-level mean values (
Table 8): (1) the future error arising from
me,
wt,
co or
sa prediction for a newly sampled tree would be expected to fall within 13%, 11%, 10% and 8%, respectively, of their true values; and (2) mean error arising from
me,
wt,
co or
sa predictions for a newly sampled stand of trees would be expected to fall within 3%, 2%, 2% and 2%, respectively, of their true values.
4.1. The Acoustic Velocity-Stiffness Relationship and Associated Inferences
The relationship between the modulus of elasticity and the density-weighted acoustic velocity of a mechanically-induced dilatational stress wave within standing trees was originally derived from engineering principles and subsequently empirically validated through field and laboratory experimentation [
39,
40]. Practically, however, the relationship is difficult to apply in the field given the logistical challenges of estimating wood density non-destructively. Thus, apart from a few studies that have incorporated surrogate measures of density obtained through the use of non-destructive tools such as the Resistograph (this study) or the Pilodyn (PROCEQ, Zurich, Switzerland; [
41]), most previous studies have either omitted the density term or assumed that it is an invariant species or sample specific constant. Analytically comparable studies to this study, such as that completed by Chen [
41], reported a significant (
p ≤ 0.05) but weak relationship (
r2 = 0.28) in this regard for Norway spruce (
Picea abies (L.) Karst.). Contrasting this result with that obtained for red pine, likewise exhibited a significant (
p ≤ 0.05) but slightly higher descriptive relationship (
r2 = 0.42). Other studies employing the simpler density-unweighted acoustic velocity variable, have reported a considerable range in the proportion of variation explained by the stiffness - acoustic velocity relationship: coefficients of determination have ranged from 0.11 to 0.41 [
39]. More complex regression models in which presumed covariates affecting acoustic velocity have been incorporated, have also been proposed. Although one of the most frequent covariates considered is tree size (diameter at breast-height), results have been mixed in terms of its influence in increasing the proportion of variability explained when included in specifications derived from Equation (1) [
15]. Analyzing the data for the 54 red pine trees assessed in this study, indicated no significant correlation between acoustic velocity and breast-height diameter (i.e.,
r = −0.018 (
p > 0.05)), thus providing confirmatory support for the deployment of the more classical representation of the modulus of elasticity − acoustic velocity relationship (sensu Equation (1)).
Previous studies have assessed the direct relationship between the dynamic
me of standing trees and the corresponding static
me of the resultant end-products as determined via the bending stress tests of dimensional lumber or laboratory assessment of clear wood specimens extracted from sawn boards. Overall, the strength of these relationships as measured by the degree of correlation or the proportion of variation explained, has also varied among studies. For example, Amateis and Burkhart [
42] found a non-significant (
p > 0.05) relationship between time-of-flight acoustic velocity impulses within standing loblolly pine (
Pinus taeda L.) trees and the static modulus of elasticity of the resultant sawn boards. Chen [
41] and Fischer [
43] reported significant (
p ≤ 0.05) but relatively weak relationships in this regard for Norway spruce (
r2 values of 0.25 and 0.13, respectively). Conversely, Wang [
12] reported significant and relatively moderately strong relationships for western hemlock (
Tsuga heterophylla (Raf.) Sarg.;
r2 = 0.73) and Sitka spruce (
Picea sitchensis (Bong.) Carr.);
r2 = 0.77).
This wide range of results among studies in terms of the significance and strength of the relationships is partially due to differences among the investigations in terms of the species examined, sampling protocols used, environmental conditions at the time of measurements (e.g., seasonal differences in temperature and moisture), locale, instrumentation (e.g., Director ST300 or TreeSonic acoustic velocity tools; Resistograph or the Pilodyn (PROCEQ, Zurich, Switzerland) wood density estimation tools), model specifications (e.g., density-unweighted or density-weighted acoustic velocity variable), correlation versus regression analysis, simple or multiple regression analyses), and availability of modulus of elasticity measures within clear wood samples from standing trees (e.g., this study) or from derived end-products (dimensional lumber). These differences are problematic in terms of drawing explicit comparisons among and between studies. Nevertheless, on a collective basis, the results presented in this study for standing red pine trees are in general agreement with the range of previous results in terms of statistical significance, explanatory performance and predictive ability. Thus providing further incremental empirical support for the generality of the relationship between density-weighted acoustic velocity and the dynamic modulus of elasticity for standing trees.
4.2. Standing Tree Versus Log Acoustic Relationships and a Poisson Ratio Estimate for Red Pine
Conceptually, the relationship between the dynamic modulus of elasticity and acoustic velocity varies between standing trees and sawn logs because the wave types being generated are different. For standing trees, it is the time-of-flight (velocity) of a mechanically-induced dilatational or quasi-dilatational stress wave that enters from the circumference of the stem just above stump height (0.3 m), progressing vertically through the xylem tissue which transects breast-height (1.3 m), and then exits the stem at a height of approximately 1.5 m. For logs, it is the velocity of a mechanically-induced resonance-based longitudinal stress wave that enters the log at one of its open cross-sectional faces (log ends), progressing horizontally through the xylem tissue until arriving at the opposite log face. Although the wave types differ along with their functional relationship with the modulus of elasticity, the velocity measurements are correlated. Specifically, Wang [
17] reported that the mean ratio between the time-of-flight acoustic velocity estimate for standing trees (
vd) and the resonance acoustic velocity estimate for butt logs (
vl) for the same sample trees across 5 coniferous species, Sitka spruce (
Picea sitchensis (Bong.) Carr.), western hemlock (
Tsuga heterophylla (Raf.) Sarg), jack pine, ponderosa pine (
Pinus ponderosa Dougl. ex Laws.), and radiata pine (
Pinus radiata D. Don), ranged from 1.07 to 1.36 with a mean value of 1.20. Based on tree and log acoustic measures for the red pine trees considered in this study, results revealed an overall mean ratio of 1.39 with individual ratios ranging from a minimum of 1.28 to a maximum of 1.56.
Although tree and log velocities are not equivalent given that they are reflecting 2 different types of stress waves (dilatational for tree and longitudinal for logs), their inter-relationship can be used to empirically estimate the Poisson ratio, which is a principal covariate in the primary acoustic relationship (Equation (1)). Mechanically, the Poisson ratio represents the transverse to axial strain relationship when a wood sample is axially loaded (i.e., ratio of the deformation perpendicular to the direction of the load (transverse strain) is proportional to the deformation parallel to the direction of the load (axial strain)). Given that the Poisson ratio varies within and between species and is affected by the moisture content and specific gravity [
44], it is commonly treated as an unknown constant when acoustically estimating wood stiffness. However, Wang [
17] demonstrated an approach for empirically deriving the Poisson ratio (
P) based on the ratio between
vd and
vl: specifically, by inputting a mean ratio value and subsequently solving for
P in the
relationship. Solving this relationship employing the acoustic ratio obtained for the red pine trees sampled in this study, yields an estimated
P value of 0.39. This value is slightly greater than the largest value reported by Wang [
17] (i.e., 0.38 for Ponderosa pine). More generally, however, a generic mean value of 0.37 is commonly assumed for both hardwoods and softwoods [
45] and hence the acoustic-based empirical
P estimate for red pine is not dissimilar. Irrespective of its numeric value relative to other species, provision of the
P estimate could be of future utility when quantifying acoustic-based relationships for red pine.
4.3. A Suite of Acoustic Velocity-Attribute Relationships and their Potential Operational Utility
The conceptual expansion of the primary
relationship to include secondary relationships, consequentially expands the acoustic-based analytical framework. These additional attributes provide for a much more comprehensive assessment of end-product potential than that based on stiffness alone. The empirical results obtained in this study for standing red pine trees, indicated that the acoustic velocity approach could be used to estimate wood density, tracheid wall thickness, fibre coarseness and specific surface area, in addition to the dynamic modulus of elasticity. This expanded set of attributes are associated with a wider range of end-products which can be used as surrogate indicators of potential end-product quality (
Table 1). Results from a parallel analysis deploying a similar inferential framework for investigating the relationship between the same attributes used in this study and velocity of a mechanically-induced longitudinal stress wave but for red pine logs derived from the same sample trees as used in this study, were in accord with those found in this study [
28]. Specifically, assessment of the relationship between acoustic velocity as measured via the Director HM200 acoustic velocity resonance tool (Fibre-gen Inc., Christchurch, New Zealand;
www.fibre-gen.com) and the same 8 Silviscan-determined attributes as used in this study, revealed viable regression relationships for the same 5 variables based on statistical significance, unbiasedness and predictive ability (i.e., modulus of elasticity, wood density, cell wall thickness, fibre coarseness and specific surface area).
Operationally, the parameterized equations can be used to generate estimates for the dynamic modulus of elasticity, wood density, tracheid wall thickness, fibre coarseness and specific surface area within standing red pine trees (
Table 4). In-forest implementation of the regression relationships will require an end-user to attain a prerequisite wood density estimate using either an acoustic velocity measurement obtained using the Director 300 time-of-flight tool (i.e.,
wd =
;
Table 4;
Figure 1), or from an amplitude measurement obtained using the Resistograph micro-drill tool (i.e.,
wd =
;
Table 5;
Figure 2). In regards to the latter approach, the assessment and quantification of the relationship between drill resistance amplitude and wood density represents an alternative field-based approach for obtaining a density estimate. The Resistograph measures the relative resistance (expressed as the percent amplitude) of a micro-diameter drill bit rotating at a constant rate and being inserted at a constant rate when drilled radially into a standing tree. Although designed for assessing the structural integrity of load-bearing wood-based structures including bridge supports, timber beams and utility poles, the Resistograph as shown in this study and by others (e.g., [
46]), can provide an indirect estimate of wood density. The correlation between the mean amplitude value derived from drill resistance profiles and oven-dry wood density has varied across studies, ranging from a low of 0.29 to a high of 0.89 ([
29,
46], respectively). The result from this study for red pine falls midway within this range (
r = 0.51; derived from
Table 5) and hence is in general agreement with these previous findings. Thus, in addition to providing a species-specific parameterized equation for predicting density within standing red pine trees, the results of this analysis provide additional incremental support for the micro-drilling approach in wood density estimation.
From a logistical perspective, however, the acoustic approach is less complex to implement. Furthermore, deploying acoustic-based wood density estimates would yield unbiased individual and stand-level attribute estimates at generally tolerable levels of precision. For example, at the individual tree-level, there was a 95% probability that a future error arising from the prediction of
me,
wt,
co and
sa using a newly acquired acoustic velocity measurement for an individual red pine tree along with the corresponding acoustic-based density estimates, would be within 13%, 11%, 10% and 8% of their true values, respectively (
Table 8). At the stand-level, there was a 95% probability that the mean error arising from the prediction of
me,
wt,
co and
sa using 30 newly acquired acoustic velocity measurements from a stand of trees along with the corresponding acoustic-based density estimates, would be within 3%, 2%, 2% and 2% of their true values, respectively (
Table 8). More generally, the tolerance intervals indicated that 95% of all future errors would be within 15%, 12%, 11% and 10% of their true
me,
wt,
co and
sa values, respectively (
Table 8).
The utility of the acoustic approach for in-forest segregation of individual trees into end-product categories and associated grade classes based on the
me,
wd,
wt,
co or
sa estimates, is ultimately dependent on the accuracy requirements of the end-user. The precision of the acoustic-based estimates for red pine as quantified by the prediction and tolerance error intervals can provide operational guidance for such a determination. For example, the mean difference for the static
me between the 14 consecutive machine stress-rated lumber grades for conifers is approximately 0.7 GPa, according to the National Lumber Grades Authority [
11]. However, the prediction interval for an acoustic-based dynamic
me estimate for an individual red pine tree is ±1.5 GPa (
Table 8). Thus, even if one assumed a 1-to-1 relationship between dynamic elasticity estimates in standing trees and static elasticity estimates within derived dimensional lumber products, the acoustic-based stiffness estimates would not be precise enough to sort standing trees into grade classes requiring such a ±0.7 GPa precision level. Alternatively, grouping the grade categories into a smaller set of discrete classes may be a viable approach. Based on the NLGA [
47] machine stress-rated lumber specifications for spruce, pine and fir boards, Paradis [
48] used 3
me-based grade classes to represent lumber end-product potentials for black spruce. These 3 classes, denoted low, medium and high grade, were differentiated by approximately 1.5
me units (GPa), and hence if similarly applied to red pine, the acoustic-based dynamic
me estimate would be precise enough to segregate individual trees into one of these 3 classes.
The lack of specific design specifications for the other predictable attributes (wood density, cell wall thickness, fibre coarseness and specific surface area) negates a similar assessment of the potential utility of the estimates to explicitly differentiate standing red pine trees into discrete end-product quality classes. However, assessing individual trees or stands using all the of the estimable attributes collectively, affords the end-user the ability to segregate standing red pine trees into coarse-level end-product-based categories (sensu
Table 1). For example, the solid wood end-product potential of standing trees could be determined from the modulus of elasticity and wood density estimates given that these attributes are directly proportional to lumber stiffness and strength, respectively. Conversely, the pulp and paper end-product potential of standing red pine trees could be inferred from the wall thickness, fibre coarseness and specific surface area estimates, given that these estimates are inversely proportional to the tensile strength, directly proportional to tear strength and inversely proportional to yield of derived paper products, respectively.
The expansion of the acoustic approach to estimate additional internal wood attributes along with provision of the prerequisite prediction equations, provides the foundation for potentially deploying an expanded inferential framework in red pine segregation operations. In terms of further research, it would be advantageous to direct efforts towards identifying and minimizing potential sources of variation influencing the attribute-specific acoustic relationships, and establishing quantitative linkages between within-tree and within-product attribute estimates. For example, establishing species-specific relationships between acoustic-based dynamic modulus of elasticity estimates within standing trees and mechanical-based static modulus of elasticity estimates within derived solid wood products, would improve in-forest end-product forecasts and provide more precise wood property characterization of the raw material before merchandizing (sensu [
49]). Ultimately, such efforts could potentially yield a more comprehensive and precise end-product segregation protocol resulting in improved in-forest decision-making within the upstream portion of the red pine forest products supply chain.