1. INTRODUCTION
Assessing wood quality is a key stage in genetic improvement programs and to determine the best use of this material. The chemical and anatomical components of wood present great variability, both among species and within the same individual. Determining and quantifying these components implies a laborious and time consuming process from the preparation of the samples to the results ( Hein et al., 2016 ).
Near-infrared spectroscopy (NIRS) is an analytical technique that employs reading the molecular vibrations resulting from the incidence of an infrared ray on some compound or sample performed at a wavelength of 750 to 2500 nm ( Pasquini, 2003 ). Variations in spectrum absorbance due to changes in the chemical composition of the elements can be quantified and used in several studies ( Campos et al., 2009 ).
NIRS has proven to be an efficient and reliable method in the most varied forms of use, such as for characterization studies of wood quality ( Duca et al., 2016 ). The technique can be used to identify molecular species in complex chemical mixtures, such as wood ( Viana et al., 2009 ; Terdwongworakul et al., 2005 ). At the beginning of its application in the forestry sector, this application was more common in the pulp and paper industry segment due to the advancement and promotion of studies in this area ( Nisgoski et al., 2016 ; Samistraro et al., 2009 ). However, many studies have been produced in several areas of wood technology such as charcoal and wood chemistry, among others ( Ramalho et al., 2017 ; Duca et al., 2016 ; Meder et al., 2017 ).
Readings of the electromagnetic radiation in the NIRS range can be performed in samples with almost no previous preparation, allowing for qualitative and quantitative analyses of the wood without the need of felling the tree, thus constituting a non-destructive method ( Schimleck et al., 2004 ). However, the application of this technique requires a database with reference values obtained by traditional methods for calibrating and validating the results found by this technique.
Due to the variability of the reading spectra, they need to undergo corrections through mathematical models in order to minimize reading errors ( Williams & Norris, 2001 ). In this context, the use of near-infrared spectroscopy (NIRS) represents a fast, reliable and inexpensive analytical method. When combined with multivariate analysis, it can correlate a variable of interest with the obtained spectra. This method is fast, taking about 90 seconds to carry out an analysis, in addition to having the advantage of not using any chemical substances and an almost total independence of operator errors.
The application of this technique to Brazilian cerrado species lies is important due to the large amount of commercialized charcoal coming from this biome. It is also important for the identification of cerrado species with potential to produce charcoal. This will contribute to the implantation of pure plantations to reduce pressure on species that do not have characteristics for this purpose.
Wood from the Sclerolobium paniculatum Vogel and Qualea grandiflora Mart. species are widely used for charcoal production and present rapid growth in regions of low soil fertility and large water deficit. Due to these characteristics, it is justifiable to study them as energy potential.
Thus, this study aimed to evaluate the feasibility of using NIRS for determining gravimetric yield, ash content, volatile materials and fixed carbon from the coal obtained from Sclerolobium paniculatum Vogel and Qualea grandiflora Mart., which are typical pioneering species of the sensu stricto cerrado.
2. MATERIAL AND METHODS
Sclerolobium paniculatum Vogel and Qualea grandiflora Mart. trees were collected from a suppression area for pasture maintenance in the municipalities of Gurupi-TO and Dueré-TO, with an average age of 7 years. Age was estimated as a function of the last suppression performed for pasture maintenance. Four individuals of each species were randomly selected taking into account the presence of heartwood (detected macroscopically) and DBH between 15 and 19 cm.
We used Eucalyptus wood as comparison, since there are many studies on this species. Therefore, four (4) samples of Eucalyptus sp. obtained from a commercial plantation located in the municipality of Gurupi-TO, with 5 years of age and mean DBH of 23 cm were used as controls.
2.1. Sample acquisition
The trees were felled at ground level and a 5.0 cm thick disc was removed from the base of each tree trunk.
The discs were divided into three positions - external, internal and central; next, 2.0 cm x 2.0 cm x 5.0 cm specimens were removed from each position of the disc ( Figure 1 ).

Figure 1 Scheme for obtaining the samples: (E) disk external position; (I) disc intermediate position and (C) disk center position.
After preparation, the specimens were submitted to the NIR spectrometer for spectra acquisition. The specimens were subsequently grinded in a Wiley mill, and sieved in a 1 mm mesh sieve. The material that passed through this sieve was also submitted to the spectra acquisition. The spectral acquisition procedures were performed at CIRAD.
2.2. Obtaining the spectra
BRUKER® equipment with an integration sphere was used to obtain the spectra. We collected one spectrum per sample. Each spectrum was obtained from 32 and 64 scans in the range of 1100-2500 nm with a resolution of 8 cm-1, obtained from the radial side of the specimens in their original state (OS) and after being ground (GS) and sieved in a 1 mm mesh.
The specimens were conditioned in air-conditioned room with a relative humidity of 65% and temperature of 21 °C, thus giving the test specimens a moisture content of 12%.
2.3. Energy characterization
For the energy characterization process, 6g of totally dry sawdust (103 ± 2 ºC) from each specimen were placed into a covered crucible in a muffle oven.
With the purpose of simulating the conventional carbonization process as close as possible, the carbonization rate was performed at 4.5 °C/minute, taking 1 hour and 40 minutes to reach the temperature of 450 ºC, and remaining for another 30 minutes at this temperature. After the muffle oven cooled, the weight of the obtained coal was determined.
To determine the gravimetric yield (GY%), the following equation was used:
In which, GY%: gravimetric yield; DWW: dry wood weight; DCW: dry charcoal weight (103 ± 2 °C).
The AFNOR 1860-2 standard ( AFNOR, 2000 ) was used to determine the content of volatile materials (VMC), ash (AC) and fixed carbon (FCC) of the produced charcoal.
2.4. Statistical analysis
The experiment was conducted using a Random Block Design (RBD) under a 3 x 3 factorial arrangement (species x positions) with four replications. Data were submitted to transformation in arcsen√ (x÷100), since they did not comply with the principle of normal distribution. The results were submitted to analysis of variance by the F test at the 5% probability level to verify the effect of variability between species on gravimetric yield, content of volatile materials, ash content and fixed carbon content. When significant, differences were identified (P<0.05) and the comparative test between means at the 5% probability level was applied by the Tukey test. Statistical analyzes were performed using the Sisvar and Microsoft Office Excel 2010 programs.
2.5. Pre-treatment of spectral data
We acquired spectral data (reflectance intensity versus wavelength) using the Unscrambler 9.7 software by CAMO A.S., according to standard ASTME1655-05 ( ASTM, 2000 ). In order to minimize variations of the spectral readings, which can affect the distribution and uniformity, the following normalization pre-treatments were tested. These tests are used to obtain the data with similar scales to obtain a more uniform distribution of the variances, and smoothing applied to the reduction of spectral noise: standard normal variation transformation (SNV) with 20-point window with quadratic adjustment; and Partial least squares regression (PLS) applied to the data matrix of the cross-validated samples (with the exception of the outliers) with the objective of correlating the spectra with the following variables: Gravimetric yield (GY), Volatile Materials Content (VMC), Fixed Carbon Content (FCC) and Ash Content (AC).
The choice for the data pre-treatment was based on the coefficient of determination (R 2) that represents the adjustment measure of a generalized linear statistical model in relation to the observed values, indicating how much the model can explain the observed values. In addition to the determination coefficient, a correlation between the studied variables and the spectra was also adopted for the choice of pre-treatments, which portrays the accuracy of the model and allows for predicted error standardization, making it possible to compare different properties ( Hein et al., 2010 ).
For the construction of the models, the range of all samples was considered to predict the variables for immediate charcoal chemistry and gravimetric yield due to the amount of samples obtained in this study.
3. RESULTS AND DISCUSSION
The following results were obtained from the energy evaluation, as presented in Table 1 .
Table 1 Averages obtained for the variables: GY, FCC, VMC and AC after the energetic evaluation of the studied wood.
Species | %GY | %FCC | %VMC | %AC |
---|---|---|---|---|
Eucalyptus sp. | 24.90Cb | 79.14 Aa | 19.10 Ca | 1.75 Aa |
Qualea grandiflora Mart | 25.96 Ba | 77.79 Aa | 20.65 Ba | 1.54 Ba* |
Sclerolobium paniculatum Vogel | 29.56 Aa | 72.05 Ba | 26.38 Aa | 1.56 Ba |
* Means followed by the same uppercase letter in the column and the same lowercase letter in the row do not differ at the 5% probability level by the Tukey test;
GY: gravimetric yield; FCC: fixed carbon content; VMC: volatile materials content; AC: ashes content.
According to the data presented in Table 1 , we observed a high gravimetric yield percentage of the species studied in relation to Eucalyptus sp.. We conclude that this high gravimetric yield, in addition to the low ash content found, shows superiority of the cerrado species in relation to the traditional species for charcoal production. This relationship was also observed by Vale et al. (2010) , who studied the gravimetric yield of charcoal produced from cerrado species and found an average of 35%. This difference may be associated with the 50% slower heating ramp in this study, which may increase the gravimetric yield. Table 2 shows the results of the calibration models for gravimetric yield.
Table 2 Multivariate calibration model (PLS) for gravimetric yield.
Type of scans* | Reference values** | Without Pretreatment | Pre-treatments | ||
---|---|---|---|---|---|
Smoothing | Normalization | SNV | |||
Ground 64 scans |
R2 | 0.3023 | 0.3327 | 0.3378 | 0.7711 |
Correlation | 0.5415 | 0.5601 | 0.5710 | 0.8738 | |
Original 32 scans |
R2 | 0.7184 | 0.7254 | 0.7352 | 0.7437 |
Correlation | 0.8467 | 0.8509 | 0.8494 | 0.8589 |
*Best results for each treatment of the specimens x number of scans;
** The reference values were determined by the coefficient of determination (R2 ) and by the correlation between the spectra and the variable in question;
PLS: Partial Least Squares; SNV: Standard Normal Variate.
Table 2 shows the gravimetric yield variable versus the calculated values, based on the NIRS spectra obtained from the wood, with the best coefficient of determination presented by standard normal variate transformation (SNV). Ground wood presented the best results.
The best coefficient of determination (R2) found for calibration was higher than 0.7 for the ground wood (0.7711), with 87.38% correlation for powder samples with the GY variable, indicating a strong correlation between the variable and the spectrum. Andrade et al. (2012) found a coefficient of determination (R2) of 0.74 with a correlation of 86% for the gravimetric yield variable in the prediction of wood and coal properties from clonal plantations of Eucalyptus sp. using near-infrared spectroscopy.
Table 3 shows that the volatile materials variable presented the highest coefficient of determination (R2) among all analyzed variables, and the spectral readings of the ground samples showed the best determination coefficients (0.88) by SNV data processing method (standard normal variation transformation). Campos et al. (2008) found a coefficient of determination (R2) varying from 0.83 to 0.96 in the prediction of the volatile material content of Eucalyptus sp. by near-infrared spectroscopy. This indicates a potential of the technique in the determination of volatile material in charcoal samples. Table 4 presents the calibration model of fixed carbon content for the studied species.
Table 3 Multivariate calibration model (PLS) of volatile materials.
Type of scans* | Reference values** | Without Pretreatment | Pre-treatments | ||
---|---|---|---|---|---|
Smoothing | Normalization | SNV | |||
Ground 64 scans |
R2 | 0.2842 | 0.2885 | 0.2940 | 0.8834 |
Correlation | 0.5244 | 0.5307 | 0.5186 | 0.9377 | |
Original 32 scans |
R2 | 0.8547 | 0.8525 | 0.8511 | 0.8661 |
Correlation | 0.9204 | 0.9188 | 0.9189 | 0.9289 |
*Best results for each treatment of the specimens x number of scans;
** The reference values were determined by the coefficient of determination (R2 ) and by the correlation between the spectra and the variable in question;
PLS: Partial Least Squares; SNV: Standard Normal Variate.
Table 4 Multivariate calibration model (PLS) of fixed carbon.
Type of scans* | Reference values** | Without Pretreatment | Pre-treatments | |||
---|---|---|---|---|---|---|
Smoothing | Normalization | SNV | ||||
Ground 64 scans |
R2 | 0.3011 | 0.2636 | 0.2617 | 0.8694 | |
Correlation | 0.5342 | 0.5000 | 0.5012 | 0.9326 | ||
Original 16 scans |
R2 | 0.8511 | 0.8521 | 0.8611 | 0.8638 | |
Correlation | 0.9209 | 0.9227 | 0.9257 | 0.9282 |
*Best results for each treatment of the specimens x number of scans;
** The reference values were determined by the coefficient of determination (R2 ) and by the correlation between the spectra and the variable in question;
PLS: Partial Least Squares; SNV: Standard Normal Variate.
In Table 4 , the best results for predicting fixed carbon contents for the charcoal production were for ground wood samples. The standard normal variation transformation (SNV) was the best pre-treatment. The correlation obtained for this variable was 93%, and the coefficient of determination (R 2) was 0.86. These results corroborate those by Campos et al. (2008) , who obtained a similar coefficient of determination in their study which determined the immediate chemistry of the charcoal varying (R2) between 0.83 and 0.96; while Andrade et al. (2012) found a coefficient of determination (R2) of 0.88, both for the prediction of the fixed carbon contents of Eucalyptus sp.
Despite representing an extremely important variable for charcoal, the ash content prediction was not successful in any of the pre-treatments for the evaluated data, which can be explained due to ash not having organic compounds in its constitution. Studies carried out by Campos et al. (2009) found that ash contents cannot be predicted through near-infrared spectroscopy in conjunction with multivariate tools.
4. Conclusions
Near-infrared spectroscopy presented satisfactory results in the prediction of gravimetric yield, content of volatile materials and fixed carbon for the evaluated species.
The calibrations obtained using the spectra obtained from ground samples presented the best correlations for the gravimetric yield, content of volatile materials and fixed carbon variables when compared to the spectra acquired from solid wood.