Development of Non-Destructive-Testing Based Selection and Grading Strategies for Plantation Eucalyptus nitens Sawn Boards

: Stiffness is considered one of the most important structural properties for sawn timber used in buildings and laminated structures including mass timber elements. There is great potential to use plantation Eucalyptus timber for structural applications, and the successful development of a plantation timber supply chain for structural products will depend on the accurate selection and grading of the resource. In this study we aimed to investigate the suitability of non-destructive testing (NDT) to improve selection and grading of sawn boards sourced from a young E. nitens plantation. We studied 268 sawn boards traced from the tree through to ﬁnal processing stages. We found high and positive correlations between stiffness (measured as dynamic modulus of elasticity) tested at each board processing stage through acoustic wave velocity (AWV) and the static board modulus of elasticity measured through mechanical testing on dressed boards. Position of the board in the stem and sawn board processing treatment signiﬁcantly impacted board modulus of elasticity, indicating that early selection of logs would allow larger yield of stiffer boards. We investigated the grading of boards through the traditional Australian Standards using a visual-grading system and through AWV, ﬁnding a classiﬁcation error of 82.5% and 45.2%, respectively. We developed a linear model which was used to re-classify the boards, obtaining a smaller classiﬁcation error, including fewer boards being over-graded. Our results demonstrate that AWV can be used as an early selection method for structural boards and can also be employed to satisfactorily grade E. nitens plantation boards to be used in building structures and as elements of mass timber.


Introduction
The increased demand for timber products is stimulating the use of alternative sources of raw forest materials, increasing the demand on plantation forestry. In the past two decades there was a rapid expansion of plantation forestry, especially in the southern hemisphere [1]. The dominant hardwoods being planted are from the Eucalyptus genus and, in Australia, eucalypts account for over 884,000 ha, which annually generate over 10 million m 3 of pulp logs [2]. The use of eucalypt timber fast grown in plantations is currently being explored as a complement to the softwood used in structural products and for native forest timber. Although bearing large potential to be used in different timber products, the characteristics of eucalypt plantation timber, its processing, and its use in mass timber elements for buildings is still a novel area of research [3].
Eucalypt timber grown in plantations is routinely harvested on 15-year rotations. Stands are established at high stocking rates to favor competition and fast growth and are managed mostly without silvicultural treatments such as thinning or pruning. Consequently, the timber derived from plantations includes certain characteristics such as knots, derived from the branching habits of the trees, drying features such as checks and 1.
Examine the effect of log position in the stem and board treatment stage on sawn boards' stiffness, measured as dynamic modulus of elasticity (MOE dyn ).

2.
Understand the correlative relationship between MOE dyn measured using NDTacoustic wave velocity at different processing stages and the actual stiffness of sawn boards (static modulus of elasticity, MOE stat ).

3.
Analyze the accuracy of the Visual-Stress-Grading (VSG) method and AWV in estimating the MOE stat of sawn boards. 4.
Improve the accuracy of VSG method in predicting the MOE stat of sawn boards using AWV.

Timber Resource
The material used in this study was sourced from a 21-year-old E. nitens plantation located in southern Tasmania, Australia (latitude 43 • 03 S, longitude 146 • 59 E). Trees were originally planted for the production of fiber and were harvested during a commercial pulpwood operation in the winter of 2018. For the study, a total of 18.5 m 3 of logs were harvested from 15 selected stems and all logs up to a small-end diameter of 185 mm were used as sawlogs for the study. Up to four logs were recovered per tree, recorded as log A, log B, log C, and log D (for the bottom, second, third and top log, respectively). A coloring pattern was applied to the end of the logs to maintain traceability during the processing stages. Logs were sawn into boards of four nominal widths of 75 mm, 100 mm, 125 mm, 150 mm, a nominal 45 mm thickness and average length of 5500 mm. A back-sawing pattern was chosen to maximize timber recovery, while retaining sapwood. A total of 268 boards were cut, block-stacked, tallied, and transported to the drying mill.

Board Treatment Stages and Testing
Prior to drying, each board was measured on 5 points along its length for width and thickness to assess individual board volume, and the board stiffness was measured using AWV with the acoustic resonance device Director HM200 (Fibre-gen, New Zealand). Each single board was tested alone on insulated trestles, to avoid transmission of the waves to the supports or to other timber. The test consisted of tapping with a hammer at one end of the board and reading the AWV value with the hand-held tool. Each board was weighed on a calibrated scale and the mass was used with the board volume to calculate the green density (GD, kg/m 3 ) of the boards. The dynamic MOE (MOE dyn, GPa) of each sawn board was then calculated as This set of measurements is referred to as 'green-stage', with moisture content (MC) of the boards being more than 25%, due to the large amount of water contained in the freshly cut timber. The boards were then air-dried for a period of fourteen months and reconditioned prior to final kiln drying. The same measurements performed on the green boards were repeated on the air-dried boards to record volume, mass, and AWV, referred to as 'air-dry-stage', in which the MC was tested at the mill to be at 12%. The boards were reconditioned and kiln dried, employing the operational drying schedule used for eucalypt material in the mill, to a nominal MC of 12%. After kiln-drying, the same measurements were performed on the boards to record volume, mass and AWV at the 'kiln-dry-stage'. The dried boards were square dressed (planed on the width) to final widths of 70, 90, 120, 140, 165 mm (five boards were over-sized as green), planed to a final thickness of 35 mm, and to an average board length of 5.5 m, maintaining the boards' identity. The boards were transported to the engineering laboratory at the University of Tasmania, where volume, mass, and AWV were measured, and this stage is referred to as 'dressed stage', where the boards' MC was at 12%. Each board was visually assessed and characteristics were recorded according to the Visual-Stress-Grading (VSG) method adopted for structural boards in Australia following AS 2082 [16]. Important features to estimate the board grade through the VSG were also noted as part of the investigation, and the total number of knots, knot type (sound/unsound), number of major knots, number of knot clusters, presence of board checking (checks deeper than 3 mm or longer than 1 4 the length of the board), presence of pith in the boards, presence of rot, wane, gum vein, and insect traces were all recorded.
Using the VSG method, we placed the sawn boards into structural grades and the relevant F-grade (stress grade) could be allocated, which corresponded to the expected MOE stat listed in AS 1720.1 [10]. The relevant grades are outlined in Table 1, reporting the expected minimum MOE values per each F-grade. After VSG, structural tests were performed. The static modulus of elasticity (MOE stat , GPa) of the boards was tested in an edge-wise four-point static bending-test, in accordance with the test procedures outlined in AS 4063.1 [17], and calculated according to Equation (2): with b and d being the thickness and the width of the board (mm) and l the span length, corresponding to 18 times the width (mm). a corresponds to 6 times the board width (mm). F 2 and F 1 correspond to 40% and 10% of the maximum load at failure point (Fmax). ϕ 2 and ϕ 1 are respectively the maximum displacement (mm) at F 2 and F 1 loads. Two samples of timber were recovered from each tested board, consistently cut from the top and bottom ends, to measure Basic Density (BD) and MC according to the procedure described in AS 1080.3 and AS 1080.1 [18,19] and the following Equations (3) and (4). The test was performed to record the density of the timber at the time of the mechanical testing as required by the Australian Standards. The average density of the timber was 578.9 kg/m 3 .
where m 1 is the mass of the sample at the time of the testing (kg), V is the volume of the sample before oven-drying (m 3 ), and m 0 is the mass of the sample after oven drying (kg). The MC and BD of each board were calculated as an average of the samples at the bottom and top ends of each board. The MOE stat values obtained were adjusted based on the MC of each board, according to AS 2878 [20].

Statistical Analysis
Statistical analyses were performed using R studio statistical software [21]. We used the Kolmogorov-Smirnov test and Levene's tests to verify the normality and homogeneity of variance in the data. We inspected the dataset for outliers using Tukey's fences method [22] and visual inspection of plots. We investigated differences in MOE dyn of the boards sourced from different positions in the stem (A, B, C, D) at different treatment stages (green, air-dry, kiln-dry, and dressed) through two-way repeated measures analysis of variance (ANOVA), comparing the means of MOE dyn of the boards at different treatment stages and between logs. We used tree as a random effect to account for non-independence of logs coming from the same tree. The repeated-measures model considered log position (four fixed levels), treatment stage (four fixed levels), and their interaction (nine fixed levels). Significance levels were kept at p = 0.05. We used multiple pairwise comparisons as post-hoc tests to identify differences between boards MOE dyn in the treatment stages and TukeyHSD test to investigate differences in board MOE dyn among log positions. We investigated differences in modulus of elasticity of the boards as measured through AWV at the green and dressed stage and modulus of elasticity measured mechanically at the dressed stage, thus comparing the MOE dyn of green or dressed boards and MOE stat of dressed boards. We used two-way ANOVA to account for log position, obtaining a repeated measure model consisting of log position (four fixed levels), measurement stage (three fixed levels), and their interaction (six fixed levels). We used multiple pairwise comparisons as post-hoc tests to identify differences between board modulus of elasticity among the two measurements. We then modelled the relationship between MOE dyn of the boards measured at different stages (green, air-dry, kiln-dry, and dressed) and MOE stat using linear regression.
To analyze the results of the visual grading (VSG) and non-destructive grading (NDT) we compared the grade classification of the dressed boards from the two system with the actual board MOE stat tested through mechanical testing. We used one-way ANOVA to compare the means of dressed board modulus of elasticity (MOE dyn and MOE stat ) using the F-grade determined through VSG as the factor variable.
We used Pearsons correlation coefficient at a probability level of 0.05 to test the correlation among board features and board MOE stat . If significantly correlated, those features could be included as variables into multiple regression modelling with MOE dyn to predict MOE stat . We used linear regression to model the relationship between MOE dyn and MOE stat and used the predicted values to classify the boards and then compared that classification with the classifications presented for VSG and AWV.

Board Treatment Stages
Data followed normality and homogeneity of variance. The dynamic board modulus of elasticity measured at different treatment stages is represented in Figure 1. There was a general increase in MOE dyn from the green to the kiln-dry stage, after which the MOE dyn recorded for the dressed stage decreased. Differences in MOE dyn were due both to the position of the log in the tree and to the treatment; however, the interaction between log position and treatment was not significant ( Table 2). Air-dry and kiln-dry MOEdyn values were significantly higher than green MOEdyn for all log positions ( Figure 1). Dressed MOEdyn was not significantly different between green MOEdyn for boards coming from the first, second, and top logs (log A t(106) = 2.53, p = 0.07,  Air-dry and kiln-dry MOE dyn values were significantly higher than green MOE dyn for all log positions (Figure 1). Dressed MOE dyn was not significantly different between green MOE dyn for boards coming from the first, second, and top logs (log A t(106) = 2.53, p = 0.07, log B t(83) = 2.31, p = 0.14, log D t(14) = 0.09, p = 1), but significantly different for boards coming from third logs (log C t(65) = 3.42, p < 0.05). Dressed MOE dyn was also significantly lower than kiln-dry MOE dyn (log A t(106) = 17.7, p < 0.001, log B t(83) = 13.9, p < 0.001, log C t(65) = 20.5, p < 0.001, log D t(14) = 6.72, p < 0.001).
There were no significant differences between modulus of elasticity measured at the green stage or dressed stage through AWV and actual modulus of elasticity measured at the dressed stage via mechanical testing (F = 1.94, p = 0.15); however, significant differences among log positions were present (F = 17.3, p < 0.001). There was no interaction between log position and measurement type (F = 0.17, p = 0.99). Table 3 reports the values of modulus of elasticity at the green and dressed stages measured via AWV (MOE dyn ) and actual modulus of elasticity measured mechanically (MOE stat ). Log position in the stem influenced the modulus of elasticity of boards ( Figure 1). Boards coming from bottom logs (log A) had consistently lower modulus of elasticity than boards from other positions in all treatment stages, but this relationship was not always significant. Although lower, the modulus of elasticity of boards from bottom logs (A) was not significantly different than that of boards from top logs (D), but significantly lower than that of boards from middle positions in the stems (log B and C).
The correlations between AWV MOE dyn measured at different stages and actual MOE stat were strong and significant at the 0.001 level (Figure 2). Considering specifically the first and last panel, measurement of MOE dyn on green boards through AWV can explain almost 60% of the variability in actual MOE stat , and measurement of MOE dyn on dressed boards through AWV can explain almost 70% of the variability in actual MOE stat . not significantly different than that of boards from top logs (D), but significantly lower than that of boards from middle positions in the stems (log B and C).
The correlations between AWV MOEdyn measured at different stages and actual MOEstat were strong and significant at the 0.001 level (Figure 2). Considering specifically the first and last panel, measurement of MOEdyn on green boards through AWV can explain almost 60% of the variability in actual MOEstat, and measurement of MOEdyn on dressed boards through AWV can explain almost 70% of the variability in actual MOEstat.

Grading Systems
The results of the allocation of boards to F-grades according to both visual grading (VSG) and AWV are presented in Tables 4 and 5, respectively. These are reported alongside the F-grade of the boards determined using the actual modulus of elasticity of the boards measured through mechanical testing (MOEstat). The VSG method had a large error, misclassifying the boards with an 82.5% rate of error. 70.2% of the board F-grades were underestimated, while 12.3% were overestimated, in respect to the actual MOEstat of the boards.

Grading Systems
The results of the allocation of boards to F-grades according to both visual grading (VSG) and AWV are presented in Tables 4 and 5, respectively. These are reported alongside the F-grade of the boards determined using the actual modulus of elasticity of the boards measured through mechanical testing (MOE stat ). The VSG method had a large error, misclassifying the boards with an 82.5% rate of error. 70.2% of the board F-grades were underestimated, while 12.3% were overestimated, in respect to the actual MOE stat of the boards. The classification error was considerably lower when utilizing the non-destructive technique of AWV to classify the boards into F-grades. A total error of 45.2% was recorded, with 22% of the boards being underestimated in their grade and 23.5% being overestimated. The indirect measure of MOE dyn via AWV is correlated with the actual MOE stat of the boards (last panel of Figure 2), which is likely the reason for a better classification of boards compared to the VSG method. Table 6 presents the comparison of AWV MOE dyn and actual MOE stat of the F-grades determined through VSG method. No significant difference was detected between the MOE stat of the grades, except for MOE stat of F17 and F11 which were found to be significantly different (p < 0.05). This further strengthens the lack of reliability of the VSG to be applied to plantation E. nitens timber in detecting the actual stiffness of the boards, either measured as MOE dyn or MOE stat . Table 5. Classification of boards via AWV and relative error. Board grades are as described in Table 4.  We found that the features of boards that are mostly correlated with actual MOE stat were board density (r = 0.66, p < 0.001), number of sound knots (r = −0.40, p < 0.001), number of knots (r = −0.37, p < 0.001), number of major knots (r = −0.27, p < 0.001), presence of pith in the boards (r = −0.38, p < 0.001), and presence of checks deeper than 3 mm (r = −0.18, p < 0.01). Although significantly correlated with the modulus of elasticity of the boards, those features did not contribute to improving the correlation between MOE dyn measured via AWV and actual MOE stat . Therefore, we modelled MOE stat only using MOE dyn as tested with AWV, without accounting further for board features. The predictions of this model were validated against the observed MOE stat values and it was found that the model could explain 69% of the variability in actual MOE stat of the boards with a RMSE of 1.26 GPa (Figure 3a). The residuals of this model were normally distributed and did not show apparent bias with the fitted values (Figure 3b).

AWV
The regression equation showed that for an increase of 0.67 GPa in MOE dyn there would be a corresponding increase in MOE stat of 1 GPa. Using this model, the actual MOE stat was predicted on the overall dataset, obtaining the classification of boards as presented in Table 7. The overall error was less than both the VSG and the classification made directly with the MOE dyn values tested through AWV, for a total of 43.3%. Using the regression equation to predict the MOE stat of the boards, 18.3% of the boards would be over-estimated, while 25% would be underestimated. In respect to the actual values recorded through AWV, the overestimation is lower, and the underestimation is higher. via AWV and actual MOEstat. Therefore, we modelled MOEstat only using MOEdyn as tested with AWV, without accounting further for board features. The predictions of this model were validated against the observed MOEstat values and it was found that the model could explain 69% of the variability in actual MOEstat of the boards with a RMSE of 1.26 GPa (Figure 3a). The residuals of this model were normally distributed and did not show apparent bias with the fitted values (Figure 3b). The regression equation showed that for an increase of 0.67 GPa in MOEdyn there would be a corresponding increase in MOEstat of 1 GPa. Using this model, the actual MOEstat was predicted on the overall dataset, obtaining the classification of boards as presented in Table 7. The overall error was less than both the VSG and the classification made directly with the MOEdyn values tested through AWV, for a total of 43.3%. Using the regression equation to predict the MOEstat of the boards, 18.3% of the boards would be overestimated, while 25% would be underestimated. In respect to the actual values recorded through AWV, the overestimation is lower, and the underestimation is higher. * Percentage of boards in F-grades determined through classification via AWV. ** Total error determined as sum of the percentage of boards misclassified (minus the percentage of boards placed into the correct F-grade). See Table 4 for example.  * Percentage of boards in F-grades determined through classification via AWV. ** Total error determined as sum of the percentage of boards misclassified (minus the percentage of boards placed into the correct F-grade). See Table 4 for example.

Board Treatment Stages
This work investigated sources of variation in stiffness, measured as modulus of elasticity of E. nitens sawn boards, and examined the potential for non-destructive testing (NDT) techniques to depict actual board stiffness and grading of boards. We found that an indirect measure of modulus of elasticity through acoustic wave velocity (AWV) at the earliest stage of log sawing corresponds well to the values that can be recorded from the dressed timber, allowing for the use of this technology as a viable substitute for cumbersome mechanical testing. Furthermore, AWV can be used as a segregation tool to select the stiffer boards early in the production chain, allowing only the best boards to be directed through the lengthy drying and reconditioning process, and re-purposing the lower-grade boards to other processing streams, thus saving resources by processing only high-quality boards.
We found that modulus of elasticity varies due to log position in the stem and board treatment stage, but there was no interaction between log position and treatment stage. The influence of log position on board modulus of elasticity and its independence from the treatments applied during timber processing is a finding of considerable impact for the timber processing industry. It suggests that log selection prior to sawmilling could increase the possibility of obtaining timber of higher stiffness, thus reducing the large costs associated with board drying, storing, treatments, and dressing. In this study, logs sourced from the bottom of the stems delivered boards of lower stiffness and values did not change considerably after board treatment and processing. During the treatment stages there were significant changes in the modulus of elasticity of the boards. There was an increase from the green to the dried stages (air-dry and kiln-dry), but after the dressing stage modulus of elasticity was not significantly different from the values measured at the green stage, when boards were freshly sawn. The apparent increase in modulus of elasticity during the drying stages is due to the impact on moisture content of boards, which affects both AWV and timber density [23].
We tested differences between modulus of elasticity measured through AWV and actual modulus of elasticity measured through mechanical testing, which we used as a benchmark test of the actual structural properties of the boards. We found that modulus of elasticity measured through AWV at the green or dressed stage did not significantly differ from the actual board modulus of elasticity. When modulus of elasticity is measured through AWV it is recorded as dynamic modulus of elasticity (MOE dyn ), which differs from actual modulus of elasticity measured through mechanical testing, recorded as static modulus of elasticity (MOE stat ). Usually, lower values of modulus of elasticity are found when measured as MOE stat in respect to MOE dyn , due to the nature of the measurement, where the static modulus of elasticity is recorded through mechanical bending tests, and the dynamic measurement is through the use of stress waves exerted into the timber [9]. MOE dyn is also influenced by the density of the sample, which in the earliest stages of processing is higher due to a larger moisture content in the timber. At the dressed stage, with the samples being tested at the same moisture content and density, we found very consistent values between the two measurements (MOE dyn and MOE stat ). The accordance between the AWV method and the mechanical testing again demonstrates the reliability of the former in detecting the stiffness of the boards. We modelled the relationship between MOE dyn measured at each stage and MOE stat , finding strong and significant correlations between the two measurements. Our results concur with previous research on differences among measurements for modulus of elasticity on boards [13,24], where dynamic modulus of elasticity measured through acoustic waves provides great reliability in measuring stiffness which would otherwise be tested through more cumbersome mechanical testing.

Grading Systems
Currently, dressed boards require appropriate grading to be classified into structural grades, and our study has compared traditional visual-stress grading (VSG) to the grading of boards via AWV. Visual-Stress Grading was designed to grade native timbers and does not account for species grown in plantation settings. The latter, such as the widely grown E. nitens, present features which are classified as impermissible in the VSG system (as large number of knots or checks) and this renders the boards to be down-graded, while their actual stiffness might be suitable for high-grades. Visual-Stress Grading of the boards in this study led to a large classification error (82.5%) with most of the boards being under-graded.
This result is aligned with previous research on the same species [11] although with a more pronounced classification error. The boards utilized in this study developed rot signs during the processing, which, although not impacting the structural properties, have affected the classification of the majority of the boards as under-grade, according to the Australian Standard, where there is little allowance for signs of rot in structural boards [9]. This led to a large number of boards being allocated to the under-grade class (UG) (39.2%), which according to the actual modulus of elasticity of the boards should have been placed in other categories. The downgrading of the boards to the UG class alone led to a classification error of almost 36%. Other characteristics, including presence of major knots and checking, severely impacted the classification of the boards, placing the majority of stiffer boards into low-grade categories, due to the restrictions in knot size and checking presence in the standard. These features might not actually impact the modulus of elasticity of the timber, and in fact, we found that the average modulus of elasticity of boards as classified through VSG was not different among grades. This result is concerning, as it highlights that applying VSG to plantation timber not only misclassifies boards either by down-grading or over-grading, but does not actually produce any meaningful grade classification.
Boards of higher grades such as F22 and F17 are sought after for construction purposes and need to adhere to a minimum stiffness threshold to be used in buildings. Hence, they are priced differently than lower grades, which might be utilized for other, non-structural purposes. Grading of timber needs to be applied with a standard that takes into account the silvicultural history and the characteristics of the timber, and VSG was originally designed for native timbers grown in remarkably different settings than plantations. Thus, the problem of how to classify plantation timber needs careful consideration to encourage use of these timbers for construction purposes. We tested the use of AWV as an alternative grading method, obtaining a considerably lower classification error (45.2%), with 22% underestimation of grade and 23.5% overestimation. This result is mostly due to the agreement between the modulus of elasticity tested through AWV and the actual stiffness of the boards, of which the relationship explained almost 70% of the variability in boards' actual modulus of elasticity. While 30% of the variation remains unexplained, this may be a suitable way to classify plantation E. nitens sawn boards. We accounted for features which are originally important in VSG, such as the presence of major knots and checks, and although they were correlated with board modulus of elasticity, they did not contribute to an increase in the predictive power of the model of modulus of elasticity measured through AWV and actual modulus of elasticity, a finding noted also on other species [25].
Using the model developed only with modulus of elasticity tested through AWV we re-classified the boards, obtaining a final classification error of 43.3%, with more boards being underestimated in their grade than overestimated (18.3% over-estimation and 25% underestimation). From an operational perspective there might be an advantage to undergrade the boards rather than over-grade, thus avoiding false attribution of boards to a higher stiffness category in respect to the actual stiffness value [11]. This result showed that acoustic measurements not only can be used as a classification tool for structural boards, but will achieve better results than the traditional visual stress grading method. Future research might focus on studying how features important for structural board use can be accounted for in a combined VSG and AWV grading method, including more features or other characteristics, to achieve the best results with minimal classification errors and satisfying the requirements of engineers and designers who would have to use the timber in structures.

Conclusions
This study presents novel findings on grading and structural properties of plantation E. nitens sawn boards. We found that board stiffness is significantly impacted by the position of the log in the stem, as well as by the treatment applied to the boards. Boards originating from logs at the bottom of the stems presented significantly lower stiffness, with boards of higher stiffness found in the middle part of the stems. Stiffness tested through AWV is a reliable measure, being highly and significantly correlated with the actual stiffness measured through mechanical testing, and can be used as a grading tool as early as the stage of log breakdown at the mill. A priori log selection would allow for only the best logs to be allocated for sawing, and after sawing only the best boards could be directed for further processing. Once dried and dressed, the grading of boards needs to take into account the peculiarities of the species grown in plantation settings, which render the traditional visual-stress grading system unsuitable for grading plantation hardwood timber. We found that AWV -based grading can be used as an alternative with remarkably low classification error. This study showed how early timber selection can improve the yield of structurally stiffer timber through fast and reliable NDT methods, which can also be used as grading tools, providing insights on possible methods which can be implemented in future grading standards and upscaled to industrial timber processing on the production lines of saw millers.