Comparison of Methods for Estimating Mechanical Properties of Wood by NIR Spectroscopy
Near-infrared (NIR) spectroscopy provides a rapid alternative to traditional methods of wood property assessment. For organizations who assess wood properties on a large scale, multisite, multispecies calibrations are of practical interest. We examined NIR spectroscopy for the estimation of density (at 12% moisture content), modulus of elasticity (MOE), and modulus of rupture (MOR) using clear wood samples obtained from several pine species (Pinus caribaea var. bahamensis, var. hondurensis, and var. caribaea, P. chiapensis, P. maximinoi, P. oocarpa, P. taeda, and P. tecunumanii). We compared different methodologies for collecting spectra, that is, benchtop instrument versus benchtop fiber-optic probe and field portable fiber-optic probe, and different wood surfaces (radial and transverse). Calibrations based on the benchtop instrument were superior to those obtained using the fiber-optic probe systems. Difficulty with adequately representing the sample when collecting spectra using a fiber-optic probe and lower quality spectra explain the differences among the data sets. Spectra collected from radial and transverse surfaces provided similar calibration statistics. The calibrations obtained for density (R2 = 0.81, SECV = 38.5 kg/m3) and MOE (R2 = 0.81, SECV = 1124 GPa) using benchtop instrument spectra demonstrate that it is possible to obtain general calibrations for estimating the wood properties of a number of tropical, subtropical, and temperate pine species.
For organizations who assess wood properties on a large scale, the cost of traditional methods of analysis limits the number of samples that can be analyzed. The practical implications of such limitations are far reaching in forestry and wood product-related industries. An understanding of wood variability for a given species, that is, within trees, and within and among stands, and among different species is severely limited as is our understanding of how different environments affect wood properties and quality. Similarly, the estimation of genetic parameters for wood properties, critical to the success of tree improvement programs, is limited by the cost and time involved in determining such properties in large numbers [1, 2].
Efforts to lower the cost and time involved in wood property analysis have seen the development of rapid, nondestructive assessment methods. Many studies have explored near-infrared (NIR) spectroscopy as an alternative approach to traditional methods of wood property analysis [3, 4]. Reports can generally be divided into two groups, that is, those related to wood chemistry, for example, cellulose or lignin content and pulp yield, and those related to physical-mechanical properties such as density, modulus of elasticity (MOE), and modulus of rupture (MOR).
In terms of understanding the performance of lumber when used in construction, density (mass divided by volume at a given moisture content), MOE (a measure of the elastic deformation that wood undergoes when subject to an applied load), and MOR (a measure of the breaking strength of wood when subject to an applied load) are critical. Current measurement standards require destructive sampling and careful sample preparation and testing. The samples used for destructive testing, referred to as short, clear samples or static bending samples, have been utilized in a range of NIR-based research over the past twenty years. In these studies, many different NIR instruments have been utilized; however, approaches to spectral acquisition differ in two key aspects. The first is the spectroscopic method utilized to collect NIR spectra; either a fiber-optic probe (referred to in the following text as “Probe”) was used or the surface of the wood sample was placed directly over the NIR beam (“Direct”). The second is the wood surface used for analysis. Owing to the structure of wood, a straight grain, clear wood sample has three distinct surfaces (the radial, tangential, and transverse or cross section) (Figure 1) available for testing. In addition to these differences, test sample size varied as did the size of the area analyzed by NIR. The area analyzed varied both in terms of the number of locations measured for a given piece of wood and the size of the NIR beam, as determined in the Direct approach by the aperture of the NIR instrument or a window placed over the beam and in the Probe method by the diameter of the fiber-optic probe. Table 1 provides a summary of these reports and includes species examined, wood surface utilized, size of the short, clear sample tested, method for collecting NIR spectra, number of spectra collected, and size of the NIR beam.
Comparison of different methodologies (Probe versus Direct or wood surface used for analysis) on the same set of samples is rare as are multiple species calibrations for physical-mechanical properties such as MOE and MOR. Therefore, the objectives of this study were as follows: (i)To develop density, MOE, and MOR calibrations using a number of tropical, subtropical, and temperate pine species(ii)To compare wood property calibrations based on NIR spectra collected using a benchtop instrument, a benchtop instrument fitted with a fiber-optic probe, and field portable instrument fitted with a fiber-optic probe(iii)To compare wood property calibrations based on NIR spectra collected from two different wood surfaces (radial and transverse) using the three different NIR instrument options
2. Material and Methods
2.1. Sample Origin
The samples utilized in this study were from the plantations managed by Brazilian companies in the Center for Genetic Conservation and Improvement of Tropical Pine and Camcore (http://www.camcore.org). Pinus caribaea var. hondurensis (Sénéclauze) W.H. Barrett & Golfari, Pinus chiapensis (Martínez) Andersen, Pinus maximinoi H.E. Moore, Pinus oocarpa Schiede ex Schltdl, Pinus taeda L., and Pinus tecunumanii F. Schwerdtf. ex Eguiluz & J.P. Perry samples were collected from a site located in Ventania, Paraná, Brazil, while a site located in the Municipality of Itararé, São Paulo, Brazil, provided samples from the remaining species (P. caribaea var. bahamensis (Grisebach) W.H. Barrett & Golfari and P. caribaea var. caribaea Morelet). Figure 2 shows the locations of the two plantations, while Table 2 provides a summary of site characteristics.
2.2. Sample Collection and Preparation of Samples
Five medium-sized trees of each species were selected, and diameter at breast height and total height were measured (Table 3). After felling, the length from the base of the tree to the merchantable top was measured and the tree cut into four logs of equal length (2.5 m). At a sawmill, a 100 mm thick bark-bark slab 2.5 m in length was cut through the pith from each log.
At the Federal University of Paraná (Curitiba, Brazil) Wood Technology Laboratory, the slabs were cut to provide clear wood samples from pith to bark with dimensions 20 mm (T) × 20 mm (R) × 300 mm (L). The samples were then conditioned to an equilibrium moisture content (MC) of 12%, and MOE and MOR were determined using an EMIC Universal Testing Machine (EMIC Equipamentos e Sistemas de Ensaio LTDA, São José dos Pinhais, PR, Brasil) according to the COPANT standard . Density of the samples was determined using the dimensions and weight of each piece at 12% MC. After completion of wood property determination, a block with dimensions 20 mm (T) × 20 mm (R) × 20 mm (L) was cut from the end of each sample for NIR analysis.
2.3. Near-Infrared Spectroscopy
Diffuse reflectance NIR spectra were collected from the radial and transverse faces of each static bending sample using a Foss NIRSystems (Laurel MD, USA) Model 5000 scanning spectrometer and an ASD AgriSpec backpack portable NIR spectrometer fitted with a contact probe (approximately 25 mm in diameter, but with a spot size of 10 mm). Spectra collected using the Foss instrument were obtained using two different methods: (1) Foss Static—samples were held in a custom-made holder that had been fitted with a 12.5 mm × 10 mm Teflon mask (the mask was used to ensure that an area of constant size was analyzed) and (2) Foss Probe—NIR spectra were collected using an Optiprobe analyzer fiber-optic probe approximately 5 mm in diameter. Owing to the size of the ASD probe, a single spectrum was collected from each surface, while for the two Foss methods, two spectra were collected and averaged (Figure 3).
Spectra obtained by the Foss NIRSystems instrument were collected at 2 nm intervals over the wavelength range 1100 to 2498 nm, while for the ASD instrument, the sampling interval was 2 nm for 1000 to 2500 nm. The full wavelength range available was used for the Foss Direct method. Owing to noise for the fiber-optic probe systems and also for comparative purposes, data collected by both spectrometers were limited to a wavelength range of 1100–2300 nm.
All measurements were made in a controlled environment of 40% relative humidity and a temperature of 20°C.
2.4. Wood Property Calibrations
Wood property calibrations for three properties: density at 12% MC, MOE, and MOR, were developed using raw spectra, that is, no spectral pretreatments were employed. The Unscrambler software (version 9.2) and partial least squares (PLS) regression (with four cross-validation segments) were used. The software recommended the number of factors reported for each calibration. A total of 334 samples were analyzed; however, two that had very high densities were identified as outliers and were omitted from the final density models.
The standard error of calibration (SEC) (determined from the residuals from the final calibration), standard error of cross-validation (SECV) (determined from the residuals of each cross-validation phase), and the coefficient of determination (R2) were used to assess calibration performance. The ratio of performance to deviation (RPDc) , calculated as the ratio of the standard deviation of the reference data to the SECV, was also used to assess calibration performance. Determination of the RPDc allowed comparison of calibrations developed for different wood properties that have different ranges in values. As this study was only concerned in comparing the relative performance of calibrations obtained using different spectroscopic approaches, it was decided to maximize the samples available for calibration and not split the datasets for calibration and prediction.
3.1. Measured Wood Properties
A summary of measured wood properties for each species is given in Table 4.
3.2. Wood Property Calibrations
Density, MOE, and MOR calibrations obtained using transverse and radial face spectra collected using the three instrument options (Foss Direct, Foss Probe, and ASD Probe) are summarized in Tables 5, 6, and 7, respectively. All calibrations reported used the number of factors indicated by the Unscrambler software, hence comparisons are not exact; however, the number of factors used are generally similar for each property and having the number of factors identical throughout did not change the ranking of instrument options or comments related to the two wood faces examined. Some calibrations required up to 10 factors, which may be considered excessive, and is possibly the result of using untreated NIR data.
For all properties, calibrations based on Foss Direct spectra provided the strongest statistics. Density (Figure 4(a)) and MOE (Figure 4(b)) calibrations had similar R2 (0.78 to 0.81) and RPDc (2.1 to 2.2) values while those for MOR (0.69 to 0.73 and 1.7 to 1.9, resp.) were lower (Figure 4(c)).
Of the two probe options, NIR spectra obtained using the ASD spectrophotometer consistently provided stronger calibration statistics than those obtained using spectra collected using the Foss Probe. For the ASD probe, R2 and RPDc values ranged from 0.64 to 0.75 and 1.6 to 2.0, respectively, while for the Foss Probe they ranged from 0.50 to 0.68 and 1.3 to 1.7. Again, density and MOE calibrations were similar, while those for MOR had weaker statistics.
For a given property, and method of collecting spectra, the transverse face may be marginally better at capturing sample variability than the radial face as RPDc values are often equivalent even though fewer factors were used for the transverse face calibrations. The only instance where a calibration based on transverse face NIR spectra was clearly inferior to the corresponding radial face spectra was for the density calibration based on data provided by the ASD probe.
Calibrations for density, MOE, and MOR were obtained using multiple pine species grown on two sites in southern Brazil. Calibrations were based on NIR spectra collected using three different approaches (Foss Direct, Foss Probe, and ASD Probe) and from two different wood surfaces (radial and transverse). The strong calibration statistics obtained for models using the Foss Direct and ASD Probe methods demonstrates the potential for general calibrations for estimating wood properties of a number of tropical, subtropical, and temperate pine species. The multiple species calibrations may not be as accurate as those based on a single site or species but are sufficient to rank samples .
Comparison of the different spectroscopic approaches used to collect spectra indicated that the Direct approach provided stronger calibration statistics than either of the Probe methods. The advantage of using the Direct approach is that the sample is exposed directly to the NIR beam (whether it is placed directly on the instrument aperture or against a Teflon window) rather than having to transmit the beam through a fiber-optic cable which causes noise in the signal beyond 2300 nm. However, to utilize a Direct approach, subsamples were often removed from the end of the short clear sample to allow the sample to be tested by the NIR instrument [13, 16, 21]. In comparison, utilization of the Probe approach provides much greater flexibility for testing samples [12, 15], and cutting the static bending sample is not required.
When comparing the ASD and Foss Probe spectra, the ASD instrument provided stronger calibrations for all properties. A critical factor in the success of the ASD Probe method was the greater surface area covered by the ASD probe compared to the Foss probe (Figure 3). Owing to its spot size, the ASD probe was able to capture more sample variability, that is, variation in earlywood and latewood components of individual rings, while it was not possible to achieve the same coverage when collecting two spectra using the Foss probe. Using a probe with a large diameter that approaches the width of the static bending sample also offers greater consistency in NIR measurements. For the smaller Foss probe, identification of sampling points was problematic as it always had the potential to over or under selection of latewood, which is critical in determining density and MOE .
An examination of related NIR studies shows that many authors either attempted to collect spectra across the complete surface of the sample, for example, Thumm and Meder , using a 20 mm wide window for samples having radial and tangential dimensions of 20 mm, or were mindful of the cumulative area they were testing, for example, Gindl et al. , noting that three spectra from different locations on the radial face of their samples covered two-thirds of the radial width. In other studies [10–12, 14, 19], radial width of static bending samples varied as the authors wished to include at least three growth rings in their samples introducing the possibility that NIR spectra collected from wider samples underrepresented the sample compared to the same number of spectra collected from narrower samples. While its significance is unclear, a related issue was encountered by Lestander et al.  who noted that because of variation in the width of the target rings, spectral information from the nearest adjacent rings were occasionally included in the collected spectra.
The transverse face had marginally better results regardless of the method used to collect spectra. The only other study that attempted to compare mechanical/physical property calibrations based on spectra collected from radial and transverse surfaces was Schimleck et al.’s study  who also reported similar results for the two surfaces. Other studies have compared calibrations based on radial and tangential surface NIR spectra. Thumm and Meder  reported MOE calibrations for radial and tangential face spectra (moving face) and found much better statistics for radial face data as did Zhao et al.  in their study of E. pellita. In comparison, Kothiyal and Raturi  reported similar calibrations for both radial and tangential face spectra obtained from E. tereticornis static bending samples. In a more recent study, Kothiyal et al.  reported calibrations based on radial and tangential surface spectra and calibrations based on an average of the two. Radial face spectra provided the strongest calibrations; however, the authors concluded that “the composite model is just as good as the radial and tangential face models and can adequately address the problem where it is difficult to obtain perfect radial or tangential surface, particularly in small girth plantation timber or where surface recognition is difficult.” One solution to the issue noted by Kothiyal et al.  is to utilize the transverse surface. In this study, tangential surface spectra were not collected owing to doubts about how well the surface would represent samples with straight grain (i.e., clear radial and tangential surfaces, as shown in Figure 1). As noted by Dahlen et al. , tangential face NIR spectra could easily under- or overestimate properties based on whether an earlywood or latewood component of a ring was present where the spectrum was collected.
Despite the strength of relationships reported in the various studies over the last 20+ years, it is unlikely that NIR-estimated MOE and MOR will replace the standard method used to determine these properties. However, NIR technology can complement the standard methods, for example, to measure applied loads , or in applications where the accuracy of a single value is not critical. Examples include studies of within-tree variation in MOE  or regional variation in MOE and MOR . In such studies, the use of NIR-estimated values can greatly improve approximations of how properties vary spatially. The success of research based on static bending samples has also led to studies aimed at estimating properties of lumber [18, 28, 31–34]. The issue of collecting representative spectra again exists, but results are promising and owing to its ability to estimate a range of properties simultaneously  NIR could have an important role in sawmills in the future. In related research, the assessment of MOE and load for wood in service has recently been demonstrated [35–37].
NIR-based calibrations for the estimation of density, modulus of elasticity (MOE), and modulus of rupture (MOR) were developed using clear wood samples obtained from several tropical, subtropical, and temperate pine species growing on two plantation sites in Brazil. NIR spectra were collected using a benchtop instrument, a benchtop instrument with fiber-optic probe, and a field portable fiber-optic probe instrument and from two different wood surfaces (radial and transverse).
Comparison of calibrations found that the benchtop instrument provided superior results, while for the fiber-optic probe systems, the field portable NIR outperformed the benchtop probe system. Difficulty with adequately representing the sample when collecting spectra using a fiber-optic probe and lower quality spectra explain the differences among the data sets.
NIR spectra collected from the radial and transverse faces provided similar calibration statistics.
The work described in this paper was presented at the 65th Forest Products Society International Convention, Portland, Oregon, USA, 19–21 June 2011.
Conflicts of Interest
The authors declare that they have no competing interests.
The authors acknowledge the financial support for this research provided by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Brazil.
L. R. Schimleck, “Near infrared spectroscopy: a rapid, non-destructive method for measuring wood properties and its application to tree breeding,” New Zealand Journal of Forestry Science, vol. 38, no. 1, pp. 14–35, 2008.View at: Google Scholar
B. K. Via, T. Shupe, L. Groom, M. Stine, and C.-L. So, “Multivariate modelling of density, strength and stiffness from near infrared spectra for mature, juvenile and pith wood of longleaf pine (Pinus palustris),” Journal of Near Infrared Spectroscopy, vol. 11, no. 1, pp. 365–378, 2003.View at: Publisher Site | Google Scholar
L. R. Schimleck, P. D. Jones, A. Clark, R. F. Daniels, and G. F. Peter, “Near infrared spectroscopy for the nondestructive estimation of clear wood properties of Pinus taeda L. from the southern United States,” Forest Products Journal, vol. 55, no. 12, pp. 21–28, 2005.View at: Google Scholar
B. K. Via, C.-L. So, T. Shupe, L. Eckhardt, M. Stine, and L. Groom, “Prediction of wood mechanical and chemical properties in the presence and absence of blue stain using two near infrared instruments,” Journal of Near Infrared Spectroscopy, vol. 13, no. 1, pp. 201–212, 2005.View at: Publisher Site | Google Scholar
E. Y. Haartveit and P. O. Flæte, “Rapid prediction of basic wood properties by near infrared spectroscopy,” New Zealand Journal of Forestry Science, vol. 36, no. 2/3, p. 393, 2006.View at: Google Scholar
T. Fujimoto, Y. Kurata, K. Matsumoto, and S. Tsuchikawa, “Application of near infrared spectroscopy for estimating wood mechanical properties of small clear and full length lumber specimens,” Journal of Near Infrared Spectroscopy, vol. 16, no. 6, pp. 529–537, 2008.View at: Publisher Site | Google Scholar
T. A. Lestander, J. Lindeberg, D. Eriksson, and U. Bergsten, “Prediction of Pinus sylvestris clear-wood properties using NIR spectroscopy and biorthogonal partial least squares regression,” Canadian Journal of Forest Research, vol. 38, no. 7, pp. 2052–2062, 2008.View at: Publisher Site | Google Scholar
R. J. Zhao, X. M. Huo, and L. Zhang, “Estimation of modulus of elasticity of Eucalyptus pellita wood by near infrared spectroscopy,” Spectroscopy and Spectral Analysis, vol. 29, no. 9, pp. 2392–2395, 2009.View at: Google Scholar
L. R. Schimleck, J. L. M. de Matos, J. T. da Silva Oliveira, and G. I. B. Muniz, “Non-destructive estimation of pernambuco (Caesalpinia echinata) clear wood properties using near infrared spectroscopy,” Journal of Near Infrared Spectroscopy, vol. 19, no. 5, pp. 411–419, 2011.View at: Publisher Site | Google Scholar
V. Kothiyal, A. Raturi, Jaideep, and Y. M. Dubey, “Enhancing the applicability of near infrared spectroscopy for estimating specific gravity of green timber from Eucalyptus tereticornis by developing composite calibration using both radial and tangential face of wood,” European Journal of Wood and Wood Products, vol. 72, no. 1, pp. 11–20, 2014.View at: Publisher Site | Google Scholar
COPANT, 555 Método de ensayo de flexion estatica, Comisión Panamericana de Normas Técnicas (Pan American Standards Commission), 1972.
C. R. Mora and L. R. Schimleck, “Determination of within-tree variation of Pinus taeda wood properties by near infrared spectroscopy. Part 2: whole-tree wood property maps,” Appita Journal, vol. 62, no. 3, pp. 232–238, 2009.View at: Google Scholar
F. Antony, L. Jordan, L. R. Schimleck, A. Clark III, R. A. Souter, and R. F. Daniels, “Regional variation in wood modulus of elasticity (stiffness) and modulus of rupture (strength) of planted loblolly pine in the United States,” Canadian Journal of Forest Research, vol. 41, no. 7, pp. 1522–1533, 2011.View at: Publisher Site | Google Scholar
H. Kobori, T. Inagaki, T. Fujimoto, T. Okura, and S. Tsuchikawa, “Fast online NIR technique to predict MOE and moisture content of sawn lumber,” Holzforschung, vol. 69, no. 3, pp. 329–335, 2015.View at: Google Scholar