Subscriber access provided by University of Newcastle, Australia
Article
POTENTIAL OF NEAR INFRARED SPECTROSCOPY FOR DISTINGUISHING CHARCOAL PRODUCED FROM PLANTED AND NATIVE WOOD FOR ENERGY PURPOSE Fernanda Maria Guedes Ramalho, Paulo Ricardo Gherardi Hein, Jéssica Moreira Andrade, and Alfredo Napoli Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.6b02446 • Publication Date (Web): 03 Jan 2017 Downloaded from http://pubs.acs.org on January 4, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
1
POTENTIAL OF NEAR INFRARED SPECTROSCOPY FOR DISTINGUISHING
2
CHARCOAL PRODUCED FROM PLANTED AND NATIVE WOOD FOR ENERGY
3
PURPOSE
4 5
*Fernanda Maria Guedes Ramalhoa, Paulo Ricardo Gherardi Heinb, Jéssica Moreira
6
Andradeb, Alfredo Napolic
7
a
8
b
9
c
10
Department of Forest Sciences. Federal University of Lavras, Lavras 37200, Brazil. Department of Forest Sciences. Federal University of Lavras, Lavras 37200, Brazil.
Center for International Cooperation in Agronomic Research for Development, Montpellier,
France.
11 12
*Corresponding Author
13
E-mail:
[email protected] 14
Phone: +55 38 988397064
15 16
ABSTRACT
17
The objective of this study was to evaluate the potential of near infrared (NIR) spectroscopy
18
associated with multivariate statistics to distinguish charcoal produced from wood of planted
19
and native forests in Brazil. Timber forest species from the Cerrado (Cedrela sp.,
20
Aspidosperma sp., Jacaranda sp. and Apuleia sp.) and Eucalyptus clones from forestry
21
companies (Vallourec steel producer and Cenibra pulp producer) were pyrolysed under well
22
controlled laboratory scale conditions at the final temperatures of 300 (573,15), 500 (773,15)
23
and 700°C (973,15°K) respectively. Fifteen charcoals of each species were produced for each
24
temperature leading to heighten controlled pyrolysis treatments and finally 270 charcoal
25
samples (3 treatments x 15 repetitions x 6 materials). Principal Component Analysis (PCA)
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
26
and Partial Least Squares Regression (PLS-R) were carried out in the spectra recorded from
27
charcoal specimens. NIR spectroscopy associated with PCA was not able to differentiate the
28
charcoals produced from native and planted woods if the 270 samples were considered in the
29
same analysis. However, the separation of native and planted charcoal was achieved when the
30
samples were analyzed separately by final pyrolysis temperature. Thus, the prediction of
31
native or planted classes by PLS-R presented better performance for samples pyrolysed at
32
300°C, followed by those at 500°C, 700°C and all together.
33
Keywords: NIR; charcoal identification; Cerrado wood; Multivariate analysis; PLS-R; PCA
34 35 36
1. INTRODUCTION Charcoal, originating from charring or slow pyrolysis of wood, with absence or
37
controlled presence of oxygen1 is an important source of energy. Among the major producers
38
of charcoal, Brazil ranks first with a production of 7.24 million tons per year.2 The produced
39
charcoal consumption is concentrated in the domestic market. The charcoal produced from
40
planted forests is a major source of energy used in reducing iron ore in the steel industry
41
mainly in Brazil.3 Brazil is the only country that produces pig iron using charcoal reduction.
42
Of the 7.74 million hectares of trees planted in Brazil, 34% are allocated to pulp and paper
43
and 15.2% for charcoal production for the steel industry.4
44
Deforestation and illegal logging in natural forests are a global concern because of the
45
threats they pose to the rich biodiversity of some regions and because of their contributions
46
towards climate change. In many tropical countries, such as Brazil, despite improved
47
enforcement against illegal logging, deforestation is still increasing. There is no official
48
information concerning the proportion of use of legal and illegal wood for energy purposes.
49
One of the great challenges is to identify illegal charcoal from legal. The distinction of
50
charcoals from planted and native forests is difficult to conduct by environmental control
ACS Paragon Plus Environment
Page 2 of 23
Page 3 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
51
agencies, making it difficult to identify fraud or false documentation. The identification of the
52
material is possible through anatomical characterization which is time consuming and
53
requires specialized professionals. In order to improve the monitoring of charcoal source, it
54
would be useful to develop fast and efficient techniques for charcoal classification.
55
Near infrared (NIR) spectroscopy is a suitable technique for characterization of various
56
materials, including wood and its products. This technique can be applied to any biological
57
material, the analysis is rapid, non-destructive and may require little or no preparation of the
58
sample to be analyzed.5,6 Few studies have applied NIR spectroscopy for evaluating charcoal.
59
In regard to quantitative approaches, Barcellos7 has used NIR spectroscopy to measure the
60
calorific value, the fixed carbon content and the volatiles content of the charcoal produced
61
from the genus Eucalyptus. Andrade et al.8 have estimated volatile material, fixed carbon and
62
gravimetric yield of charcoal obtained from Eucalyptus wood by NIR spectroscopy.
63
Concerning the qualitative analyses, Monteiro et al.9 have associated NIR spectroscopy with
64
principal component analysis (PCA) for discriminating carbonization processes (experimental
65
and industrial) used to produce Eucalyptus charcoal. Nisgoski et al.10 have applied NIR
66
spectroscopy to differentiate wood and charcoal of Brosimum acutifolium and Ficus citrifolia
67
(family Moraceae) and Hieronyma laxiflora and Sapium glandulosum (family
68
Euphorbiaceae). They successfully distinguish those species and families by PCA from NIR
69
spectra recorded from wood; However, using the NIR spectra taken from charcoal, Nisgoski
70
et al.10 could distinguish families but not species. Davrieux et al.11 have verified the feasibility
71
of using spectroscopic reflection in the mid-infrared (MIR) and NIR region for discriminating
72
charcoal produced from Tabebuia serratifolia and Eucalyptus grandis. Their results allowed a
73
successfully discriminant function for two species from their MIR and NIR spectra; where the
74
differentiation was more pronounced using MIR data. Muniz et al.12 also have applied the
75
NIR associated with PCA for discrimination of wood and charcoals of four forest species:
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
76
three native species (Mimosa scabrella, Tabebuia capitat and Hymenaea aurea) and one
77
species considered exotic (Eucalyptus alba) in Brazil. The authors have reported the need to
78
test the efficiency of this technique in different forest species for developing robust statistical
79
models able to identify tree species and, perhaps, the precursor wood used for producing the
80
charcoal.
81
In all of the studies mentioned above, there was no variation in the final wood
82
carbonization temperatures. As final temperature is important in charcoal quality, this aspect
83
deserves further investigation.13 Moreover, it would be useful to include charcoals produced
84
from different wood species in the analyses in order to verify the efficiency of the proposed
85
approach. Thus, the aim of this study was to evaluate the potential of NIR spectroscopy
86
coupled with multivariate analyses to distinguish charcoal produced from planted and native
87
wood and the influence of final carbonization temperature on classifications estimated by
88
NIR-based models.
89
2. MATERIALS AND METHODS
90
2.1. Plant material. Wood species from the Cerrado biome and reforestation were used
91
in this study. The Cerrado species were Cedrela sp. (Cedar, labelled as “C”), Aspidosperma
92
sp. (Peroba labelled as “P”), Jacaranda sp. (Rosewood, labelled as “J”) and Apuleia sp.
93
(Garapa labelled as “T”). In regard to the commercial plantation, two Eucalyptus clones
94
coming from two forestry companies, Vallourec (clones of Eucalyptus grandis x E. urophylla
95
hybrids, 6.5 years old) with a focus on production of charcoal for energy and Cenibra Nipo-
96
Brasileira (clones of Eucalyptus grandis x E. urophylla hybrids, 6 years old) focused on the
97
production of pulp and paper. The four native species occur very frequently in the Brazilian
98
Cerrado while the two types of Eucalyptus materials were chosen in order to cover a large
99
genetic variation existent among reforestation materials in Brazil.
ACS Paragon Plus Environment
Page 4 of 23
Page 5 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
100
2.2. Specimen preparation. Central planks were removed from trees and 45 specimens
101
(defect free) were obtained for each material. Specimens presenting the dimensions of 3.5 cm
102
x 3.5 cm x 10 cm (R x T x L) and 2.5 cm x 2.5 cm x 10 cm (R x T x L) were produced from
103
the native and reforested wood, respectively. All specimens were properly identified using a
104
special pencil (labeling did not disappear after pyrolysis). Before pyrolysis, specimens were
105
kept in an acclimatized room and until reaching 12% moisture.
106
2.3. Pyrolysis process. Wood samples were charred in a Macro ATG oven (Fig. 1)
107
developed by a researcher from the Center for International Cooperation in Agronomic
108
Research for Development (CIRAD, France) in partnership with Federal University of Lavras
109
(UFLA, Brazil). The Macro ATG is a prototype equipped with an electric oven where the
110
temperature can reach up to 1,000 degrees Celsius. The system is provided with a pyrolysis
111
reactor pressure controller, condenser pyroligneous liquid, load cell, and gas chromatograph
112
flow meter. The experiments can be performed with the introduction of various gases such as
113
N2, O2, CO2, CO and H2, simulating different conditions of partial combustion or complete
114
pyrolysis in the presence of an inert atmosphere. A control panel and software were
115
developed by CIRAD, specifically for use in trials in the Macro ATG oven14. Fifteen (15)
116
specimens of each plant material (4 native and 2 Eucalyptus types) were carbonized at three
117
different final temperatures (300°C/ 573,15°K, 500°C/ 773,15°K and 700°C/ 973,15°K)
118
resulting in 270 charcoal samples (6 materials x 3 temperatures x 15 repetitions). The
119
pyrolysis conditions were:
120
• Initial temperature: 40°C.
121
• Heating rate: 5°C min-1.
122
• Final temperature: 300ºC, 500ºC and 700ºC.
123
• Residence time at the final pyrolysis temperature: 1 hour.
124
• Atmosphere: nitrogen gas - N2.
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
125
• Gas flow (N2): 20 NL.min-1 (Normal Litre/min)
126
• Cooling period by natural convection: 15 hours.
127
Figure 1
128
Three temperatures were used for representing the thermal variation in an industrial
129
oven. For each carbonization, the specimens were placed in the crucible located within the
130
pyrolysis reactor. Four thermocouples were placed around the specimens to verify the desired
131
temperature within the system. The condensable gases produced during the thermal
132
decomposition process were collected through a condenser connected to the oven. After
133
cooling the furnace, the produced charcoals were removed from the crucible and taken to an
134
air-conditioned room to stabilize the humidity.
135
2.4. NIR spectra acquisition. A Bruker FT-NIR (model MPA, Bruker Optik GmbH,
136
Ettlingen, Germany) was used in diffuse reflectance mode. This Fourier transform
137
spectrometer is designed for reflectance analysis of solids with an integrating sphere that
138
measures the diffuse reflected infrared energy from a 150 mm² spot. Spectral analysis was
139
performed within the 12,500-3,600 cm-1 range at 8 cm-1 resolution (each spectrum consisted
140
of 2,307 absorption values). Each NIR spectrum was obtained with 32 scans; means were
141
calculated and compared to the sintered gold standard used as background to obtain the
142
absorption spectrum of the sample. Two NIR spectra were recorded in the center of the each
143
side of the transversal surfaces of wood and charcoal samples. The spectrometer was
144
connected to a computer which stored the spectra data collected through the OPUS program,
145
version 7.5. The NIR spectra were considered without pre-treatment and after mathematical
146
treatment by first derivative for these multivariate analyses.
147
2.5. Multivariate statistical analysis. Principal Component Analysis (PCA) and Partial
148
Least Squares Regression (PLS-R) were calculated by means of the Unscrambler (CAMO
149
AS, Norway, v. 9.7) software. First, PCA was used for compressing the main information in a
ACS Paragon Plus Environment
Page 6 of 23
Page 7 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
150
set of variables into a lower number of new variables and for exploring the dependence of the
151
studied materials by means of clusters. The PCA were calculated using full model size,
152
maximum of 6 latent variables (LV’s) and a cross-validation model was adopted (6 segments
153
of 30 NIR spectra per segment). Outliers were identified from the Residuals X-variances and
154
leverage value plot analyses. The number of LV’s adopted for each model corresponded to the
155
first minimal residual X-variance.
156
In a second approach, a discriminant analysis was carried out using partial least squares
157
regression (PLS-R) to estimate the category (native or planted) of the charcoal samples. For
158
the PLS regressions, two variables were attributed to charcoal samples: 0 for native charcoals
159
and 1 for Eucalyptus charcoals. PLS regressions were then carried out in four different data
160
sets: NIR spectra taken from charcoals produced at 300, 500 and 700°C individually (90
161
samples each, 6 x 15) and an another data set including all temperatures together (6 x 15 x 3;
162
270 samples, 540 NIR spectra). In a first step, four PLS-R models were generated and cross-
163
validated. Thus the category (0 or 1) was estimated. The estimated values (rational number)
164
the PLS regression were rounded to 0 or 1 and used as categorical variables. In a second step,
165
an independent validation of the PLS-R models was performed using 20% of total samples
166
(36 samples selected at random) for validating the models of the charcoal produced at 300°C,
167
500°C, 700°C and all together.
168
As estimates of the classes also include numbers less than zero and greater than one, the
169
estimates were rounded and samples presenting estimated values above 0.5 were considered
170
in Category 1 (charcoal from planted woods) while samples whose values were estimated
171
below 0.5 were considered as 0, belonging to the class of native woods. Each PLS-R model
172
was evaluated by the coefficient of determination (R²), the mean square error (RMSE) and the
173
percentage of correct classification.
174
3. RESULTS
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
175
3.1. NIR spectra of charcoals produced at 300, 500 and 700°C. The normal
176
(untreated) NIR spectra measured on charcoal samples produced at 300, 500 and 700°C (Fig.
177
2A) show three distinct spectral signature patterns. It is possible to see that there is noise in
178
the signal at the third overtone region of the untreated spectrum while the first derivative of
179
the NIR spectra (Fig. 2B) improved the signal-to-noise ratio. Delwiche and Reeves15 have
180
demonstrated that the application of pre-treatments on the NIR spectra helps to optimize the
181
NIR information and derive models from them.
182
Figure 2
183
3.2. Principal component analysis (PCA). Fig. 3 shows the scores of the PCA carried
184
out from carbonized samples in the three heat treatments. The three groups formed in this
185
analysis were organized according to the final carbonizing temperature, indicating that the
186
NIR is capable of detecting parameters of the carbonization process of the charcoal
187
specimens. The first principal component (PC) explained 92% of the spectral variability of
188
charcoal and the second PC explained the remaining 6%, taking into account 98% of the total
189
variation (Fig. 3). Monteiro et al.9 also investigated charcoal of seven Eucalyptus species and
190
twenty native species from the cerrado of Minas Gerais, Brazil. They presented clusters from
191
PCA of the NIR spectra showing that it is possible to distinguish the carbonization processes
192
of the samples, but no cluster was formed according to the species.
193
Figure 3
194
It was possible to clearly distinguish two clusters in the PCA two-dimensional plots
195
carried out from the wood (Fig. 4A) and charcoal (Fig. 4B, C and D) sample spectra by final
196
pyrolyses temperature. In the four charts, the two groups are represented by the PCA scores of
197
the samples coming from planted forests (EV and EC) and from native species (J, P, C and T).
198
In regard to charcoal analyses, the lower the final carbonization temperature, the clearer the
199
separation of the two groups (charcoal from planted and native wood). Few samples of
ACS Paragon Plus Environment
Page 8 of 23
Page 9 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
200
Cedrela charcoal (C) fall within the Eucalyptus group, especially in analyses of charcoal from
201
500 and 700 degrees.
202
Figure 4
203
The PCA performed from the first derivative of the NIR spectra was also able to
204
differentiate the charcoal samples into two major groups (planted or native) when only
205
charcoals produced at 300°C were considered (not shown). However, it was not possible to
206
clearly separate the NIR scores of charcoal from 500 and 700 degrees into groups, because the
207
points were randomly distributed along the graphics.
208
3.3. Classification of charcoal origin by multivariate regression. Partial Least
209
Squares Regression (PLS-R) was held from NIR spectral data to develop models for
210
predicting the classification of charcoal samples in the native or planted group. As dependent
211
variable, the value 0 (zero) was assigned to charcoal of native wood and the value one (1) to
212
charcoal of planted wood. Four regressions were obtained to estimate the class of charcoals
213
from 300, 500, 700 degrees and all temperatures. The statistics associated with models and
214
cross and independent validations are presented in Table 1.
215
Table 1
216
PLS-R models developed from the charcoal sample set of 300°C yielded the lowest root
217
mean square error and highest correlation values between real and estimated classifications
218
for cross-validation, as well as for independent validation (Table 1). It can be seen that the
219
R²cv and R²p decrease and RMSE of the models increases as the carbonization temperature
220
increases.
221
Charcoal classifications were performed from the estimates of charcoal classes
222
calculated by cross- and test set-validation models (Table 1). The estimation of class was
223
made on charcoals produced at 300, 500 and 700 degrees and all charcoal together. In this
224
approach, the samples with estimates lower than 0.5 were classified as belonging to the
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
225
category of native wood charcoal while samples whose estimates were equal to or greater than
226
0.5 were classified as reforestation charcoal. Fig. 5 (upper chart) shows the classifications
227
according to the estimates of the cross-validated model. This classification was performed
228
using 15 charcoal samples for plant material produced at the three thermal treatments and for
229
the 45 samples of each plant material taking into account all treatments. The lower chart of
230
Fig. 5 shows the classification estimates of the 36 samples used in the independent validation
231
of each model and the 108 samples of the global model.
232
Higher percentage of correct classifications occurred for the samples carbonized at
233
300°C both for models obtained by cross-validation (Fig. 5A, upper chart), as well as by
234
independent validation (Fig. 5A, lower chart). Only one incorrect classification (red point)
235
was done from the model generated with samples of 300 degrees. However, the estimates of
236
classification obtained from the models of charcoal of 500°C and 700°C and the global model
237
showed greater number of incorrect classifications.
238
Figure 5
239
The details of the number of correct classifications and the charcoal samples more
240
difficult to correctly classify are shown in Table 2. The cross-validated model developed with
241
charcoals produced at 300°C presented 99.2% of correct classifications for native wood
242
charcoals and 100% for reforestation wood charcoal. The percentage of correct classifications
243
decreases as the final carbonization temperature increases. The classifications made from the
244
models validated by an independent validation also showed this trend, but the estimates
245
present a higher percentage of correct classifications compared to class estimates of the cross-
246
validated models.
247 248
Table 2 4. DISCUSSION
ACS Paragon Plus Environment
Page 10 of 23
Page 11 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
249
Energy & Fuels
This study showed that it was possible to distinguish charcoal originating from native
250
and planted wood by combining NIR spectroscopy and multivariate statistics, since the
251
analysis is performed on samples carbonized under the same final temperature. NIR
252
spectroscopy associated with PCA was not able to differentiate charcoals of native or
253
reforested wood simultaneously using charcoal samples produced at different final
254
temperatures. On the other hand, NIR spectra can successfully be used for detecting the
255
temperature at which the charcoal was produced, indicating some information about its
256
quality as an energy source.
257
4.1. NIR spectra of charcoal varies during pyrolysis process. The NIR spectra
258
obtained from charcoal from 300°C (Fig. 2) presented a spectral signature similar to that of
259
wood. In comparison with the studies based on NIR spectra and wood6 the spectra of
260
vegetable charcoal show little absorption in the NIR region. 7, 11, 12, 10
261
As the final carbonization temperature increases, the NIR spectra of the charcoal
262
samples changes and acquires a signature that does not resemble the spectrum of wood that
263
originate the charcoal (not shown). This reflects the thermochemical conversion of wood
264
constituents during pyrolysis.10 At 300°C, there are polymers, like cellulose and lignin, which
265
have not undergone partial or complete conversion into charcoal.16 Due to this fact,
266
interaction between radiation and the chemical constituents of charcoal from 300°C is more
267
informative and it is easier to use its spectral information to distinguish the precursor raw
268
materials. Therefore, it is difficult to differentiate charcoal samples in NIR spectroscopy when
269
the charcoal is produced at high temperatures, losing or homogenizing the chemical
270
information of the chemical carbon link in charcoals. This point may be explained by the fact
271
that the conversion of the main components of wood (extractives, hemicellulose, cellulose and
272
lignin) are converted to charcoal at different temperatures. Pyrolysis at 300°C is mainly
273
responsible for the hemicellulose conversion and a large part of the extractives products. On
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
274
the other hand, cellulose pyrolysis reactions start at 300°C with a fast kinetic conversion.
275
Lignin conversion is reached at a wide range of temperature. 4.2. Analysing charcoal samples produced under different final temperatures at
276 277
the same time. The groupings resulting from PCA of the original (untreated) NIR spectra
278
recorded on charcoal showed better results than those obtained from first derivative NIR
279
spectra (not shown). However, the PCA made from all samples together (carbonized at the
280
three final temperature levels) was not adequate to distinguish the plant material, since the
281
effect of the spectral change caused by the final carbonization temperature in the NIR spectra
282
is greater than the effect of chemical composition variation in the precursor raw material (Fig.
283
3).
284
Taking into account the sensitivity of NIR spectroscopy in clearly separating charcoal
285
samples by final pyrolysis temperature (Fig. 3), it is possible to use this model to
286
preliminarily identify the final temperature at which the charcoal was produced. In a
287
subsequent step, the inspection agent could compare the spectrum of the unknown sample
288
with the database of the temperature that the sample (supposedly) belongs. Thus, depending
289
on the source group in which the charcoal sample spectrum falls, it will be possible to classify
290
it as originating from native or planted wood. The ability of NIR spectroscopy in predicting
291
the final temperature at which the charcoal is produced can be useful as a tool for ranking
292
charcoal in terms of quality. Many studies have shown that charcoal produced at high
293
temperatures present higher fixed carbon content17 and stiffness18. These properties are
294
important for charcoal to be used as bioreducing agents in blast furnaces for steel
295
production19.
296 297
4.3. Discrimination of wood and charcoal specimens by thermal treatment. The graph in Fig. 4 A indicates that it was possible to clearly distinguish samples of native wood
ACS Paragon Plus Environment
Page 12 of 23
Page 13 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
298
from planted wood. This separation shows the sensitivity of NIR spectroscopy in separating
299
wood types according to the plant material.
300
The separating of charcoal samples produced from these woods into groups was also
301
possible (Fig. 4B, C and D). However, as the final carbonization temperature increases, the
302
variations in the charcoal chemical composition are reduced and therefore the technique's
303
ability to separate the groups is reduced.
304
This may be explained by the degradation of the polymers that make up the samples
305
according to their thermal decomposition during pyrolysis. On the other hand, the polymers
306
are preserved in wood and they interact with the electromagnetic radiation generating NIR
307
spectra with information useful to distinguish between materials (Fig. 4A).
308
The distribution of samples within each temperature (Fig. 4) presents differences in
309
terms of homogeneity. The sample carbonized at 300°C are distributed more heterogeneously
310
compared with samples from 500°C and 700°C, that present homogeneous sample
311
dispersion. This can be explained by the degradation of polymers. The degradation is more
312
advanced at higher temperatures and thus the material becomes more chemically
313
homogeneous (greater concentration of carbon). Therefore, it is not possible to clearly
314
separate the carbonized plant materials by temperature.
315
The PLS regressions generated for classifying the origin of charcoals of 300°C showed
316
better statistical performance (higher R² and lower RMSE, Table 1) and the percentage of
317
correct classifications (Table 2) when compared to other models. However, all predictive
318
models showed potential to be applied in the classification of charcoal samples by regulatory
319
agencies.
320
4.4. Limitations and need for further research. The specimens of Cedrela sp. had a
321
greater number of samples incorrectly predicted from all NIR based models. This behavior
322
shows that the charcoal of some native species can be confused with charcoal planted planted
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
323
wood. The similarity between some native and reforested wood can be an obstacle in
324
implementing this approach in real situations, since the chemical composition of some woods
325
are similar. In the PCA, some spectra of Cedar charcoal were also mistakenly grouped in the
326
planted charcoal group.
327 328 329
The approach of this work is a preliminary study to develop a system to identify the source of charcoal that is produced, transported and sold illegally in many countries. Further studies using this approach may include more native and planted forest species
330
in the models. Moreover, models may be developed with charcoal samples produced at
331
different temperatures and from different carbonization kilns in order to generate a
332
comprehensive database that can be effectively used in identifying the source of charcoal.
333
4.5. Perspectives. Forestry companies are seeking genetic improvement of their clones,
334
so that these trees present greater uniformity in their properties and higher quality for the
335
production of charcoal. As can be seen in Fig. 4, the homogeneity of charcoal was more
336
pronounced depending on the final temperature, not only because of genetic variation of
337
materials. This result suggests that companies should pay more attention to the parameters of
338
the carbonization process, such as the final temperature, than just the quality of its raw
339
material, if they want to get a more homogeneous charcoal quality.
340
From this study, it is suggested that the charcoal source be evaluated using a sequential
341
analysis: initially, the final charcoal carbonization temperature is predicted, then the type of
342
wood used to produce the charcoal can be more easily estimated by means of a model
343
calibrated for such final temperature. This methodology may be useful in the actions of
344
regulatory agencies for identifying illegal commerce, contributing to the sustainability of
345
areas not (yet) deforested and the preservation of biodiversity.
346
5. CONCLUSIONS
ACS Paragon Plus Environment
Page 14 of 23
Page 15 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
The association of NIR spectroscopy with multivariate statistics allowed the
347 348
identification of the temperature at which the charcoal was produced. NIR spectroscopy also
349
could successfully detect the type of wood (native and planted) used for producing the
350
charcoal, since charcoals produced under the same final temperature are considered.
351
The prediction of classes, native or planted, by PLS regression yielded a higher
352
percentage of correct classifications for charcoal produced at 300°C, followed by the
353
estimations calculated from NIR models of charcoals from 500°C and 700°C and from the
354
model done simultaneously with all charcoal samples.
355
The higher the final carbonization temperature of the charcoal, the lower the efficiency
356
of predictive models based on NIR spectroscopy for distinguishing charcoals in regard to the
357
wood which originated them.
358
ACKNOWLEDGEMENT
359
The author expresses your special thanks to the Department of Wood Science and Technology
360
of the Universidade Federal de Lavras (UFLA, Brazil) for supporting the experimental work
361
and to the Centre de Cooperation Internationale en Recherche Agronomique pour le
362
Development (UPR114 of CIRAD, Montpellier, France) for laboratory facilities. This study
363
was funded by CNPq (National Council for Scientific and Technological Development,
364
Brazil), CAPES (Higher Education Personnel Improvement Coordination, Brazil) and
365
FAPEMIG (Foundation for Research Support of the State of Minas Gerais, Brazil).
366 367
REFERENCES
368
(1) Wenzl, H. F. J. The chemical technology of wood, 1th ed.; Academic Press: New York;
369
1970.
370
(2) FOOD AND AGRICULTURE ORGANIZATION OF THE UNITED NATIONS.
371
Forestry: forestry production and trade. Rome, 2014.
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
372
(3) Producer of brazilian association of forests planted. Statistical Yearbook of ABRAF 2013
373
years based 2012. Brasília, 2013.
374
(4) INDÚSTRIA BRASILEIRA DE ÁRVORES. Relatório IBÁ 2015. São Paulo, 2015.
375
(5) Pasquini, C. Near infrared spectroscopy: fundamentals, practical aspects and analytical
376
applications. Journal of the Brazilian Chemical Society 2003, 14, 198-219.
377
(6) Tsuchikawa, S.; Schwanninger, M. A review of recent near-infrared research for wood and
378
paper. Applied Spectrosc Reviews 2013, 48, 560-587.
379
(7) Barcellos, D. C. Characterization of the charcoal through the use of spectroscopy in the
380
near infrared 2007. 129 p. Thesis (PhD in Forest Science) - Federal University of Viçosa,
381
Viçosa, MG, 2007.
382
(8) Andrade, C. R.; Trugilho, P. F.; Hein, P. R. G.; Lima, J. T. ; Napoli, A. Near infrared
383
spectroscopy for estimating Eucalyptus charcoal properties. Journal of Near Infrared
384
Spectroscopy 2012, 20, 657-666.
385
(9) Monteiro, T. C.; Da Silva, R. V.; Lima, J. T.; Hein, P. R. G.; Napoli, A. Use of near
386
infrared spectroscopy to distinguish carbonization processes and charcoal sources. Cerne
387
2010, 16, 16-381.
388
(10) Nisgoski, S. ; De Muniz, G. I. B. ; Marrone, S. R.; Schardosin, F. Z., França, R. F. NIR
389
and anatomy of wood and charcoal from Moraceae and Euphorbiaceae species. Ciência da
390
Madeira 2015, 6, 183-190.
391
(11) Davrieux, F.; Rousset, P. L. A.; Pastore, T. C. M.; De Macedo, L. A.; Quirino, W. F.
392
Discrimination of native wood charcoal by infrared spectroscopy. Química Nova 2010, 33,
393
1093-1097.
394
(12) De Muniz, G. I. B.; Carneiro, M. E.; Nisgoski, S.; Ramirez, M. G. L.; Magalhães, W. L.
395
E. SEM and NIR characterization of four forest species charcoal. Wood Science and
396
Technology 2013, 47, 815-823.
ACS Paragon Plus Environment
Page 16 of 23
Page 17 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
397
(13) Kinney, T. J.; Masiello, C. A.; Dugan, B.; Hockaday, W. C.; Dean, M. R.; Zygourakis,
398
K.; Barnes, R. T. Hydrologic properties of biochars produced at different temperatures.
399
Biomass and Bioenergy 2012, 41, 34-43.
400
(14) Jesus, M. S.; Napoli, A.; Andrade, F. W. C.; Trugilho, P. F.; Rocha, M. F. V; Gallet, P.;
401
Boutahar, N. Macro ATG Kiln: gaseous flow study in the pyrolysis process of Eucalyptus
402
Brazilian. Journal of Wood Science 2015, 6, 269-274.
403
(15) Delwiche, S. R.; Reeves, J. B. The effect of spectral pretreatments on the PLS modeling
404
of agricultural products. Near Infrared Spectroscopy Journal 2004, 12, 177-182.
405
(16) Ponder, G. P.; Richards, G. N. A review of some recent studies on mechanisms of
406
pyrolysis of polysaccharides. Biomass and Bioenergy 1994, 7, 1-24.
407
(17) Couto, A. M.; Trugilho, P. F.; Napoli, A.; Lima, J. T.; Silva, J. R. M.; Protásio, T. P.
408
Quality of charcoal from Corymbia and Eucalyptus produced at different final carbonization
409
temperatures. Scientia Forestalis 2015, 43, 817-831.
410
(18) Assis, M. R.; Brancheriau, L.; Napoli, A.; Trugilho, P. F. Factors affecting the mechanics
411
of carbonized wood: Literature review. Wood Science and Technology 2016, 50, 519-536.
412
(19) Griessacher, T.; Antrekowitsch, J.; Steinlechner, S.; Charcoal from agricultural residues
413
as alternative reducing agent in metal recycling. Biomass and Bioenergy 2012, 39, 139-146.
414 415 416 417 418 419 420 421
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
422 423 424
425 426
LIST OF FIGURES
Figure 1 - Macro ATG Oven.
427 428
429 430
Figure 2 - Original (A) and first derivative NIR spectra (B) recorded on charcoal samples
431
produced at 300°C, 500°C and 700°C from 6 types of wood (4 native and 2 Eucalyptus).
ACS Paragon Plus Environment
Page 18 of 23
Page 19 of 23
432 433 434 EC
EV
T
C
J
P
8 6
300°C 700°C
4 Principal Component 2 (6%)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
2 0 -2 -4 -6 -8 500°C
-10 -25
-20
-15
-10 -5 0 5 Principal Component 1 (92%)
10
15
20
435 436
Figure 3 - Principal Component Analysis of the NIR spectra of charcoal samples produced at
437
300°C, 500°C and 700°C.
ACS Paragon Plus Environment
Energy & Fuels
3
3 Native
2 2 1 1 PC 2 (6%)
PC 2 (17%)
0 -1
0
Planted
-2
-1 -3
Planted
Native
-2
-4
A
Wood
B
300°C
-5
-3 -5
-3
-1
1 PC1 (81%)
3
5
7
-15
-10
-5
0
5 10 PC 1 (94%)
15
20
25
2 6
Native
4
1 Planted
PC 2 (6%)
2 PC 2 (22%)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 20 of 23
0
-2
0 Planted
-1
Native
-4
-2 -6
C
500°C -8
438
D
700°C -3
-11
-6
-1 PC 1 (77%)
4
9
-10
-5
0 5 PC 1 (93%)
10
15
439
Figure 4 – Bi-dimensional plot of the principal component (PC) analysis scores of NIR
440
spectra recorded on wood (A) and charcoal samples produced at 300°C (B), 500°C (C) and
441
700°C (D) from native (T, C, J and P) and planted (EC and EV) wood.
442 443
ACS Paragon Plus Environment
Page 21 of 23
444
Cross validated PLS Models 1,5 700°C
All
Planted
500°C
1,0
0,5
Native
Prediction of class by NIR
300°C
0,0
-0,5
-1,0
Test set validated PLS Models 1,5 500°C
700°C
All
Planted
300°C 1,0
0,5
Native
Prediction of class by NIR
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
0,0
-0,5
A
B
C
D
-1,0 Native
Planted
Native
Planted
Native
Planted
Native
Planted
445 446
Figure 5 – Classification of specimens as charcoal of native or planted wood through cross-
447
validated (above) and test set validated (below) PLS-models. Red circles represent specimens
448
incorrectly classified as planted, while yellow circles are the specimens mistakenly classified
449
as native. Classifications were estimated by models developed for charcoal produced at 300°C
450
(A), 500°C (B), 700°C (C) and from all specimens (D). Reference values: zero (0) for
451
charcoal of native wood and one (1) for charcoal of planted wood.
452 453
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 22 of 23
454
List of Tables
455 456
Table 1 – Statistics associate to PLS-R models for estimating the classification of charcoal
457
samples according to its origin by thermal treatment and analyzing all treatments together. Model 300°C 500°C 700°C All
R²cv 0.817 0.729 0.561 0.479
RMSECV 0.204 0.246 0.314 0.342
R²p 0.796 0.744 0.562 0.475
RMSEP 0.227 0.244 0.306 0.342
LV 6 4 3 4
458
R²cv: Coefficient of determination for cross-validation. RMSECV: root mean square error for
459
cross-validation. R²p: Coefficient of determination for validation using a independent test set.
460
RMSEP: root mean square error for test set validation. LV: Latent variables.
ACS Paragon Plus Environment
Page 23 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
461
Energy & Fuels
Table 2 – Prediction of sample as being native or planted, by means of PLS-R models of each heat treatment and for all samples together.
Validationmethod
Sample set
Type
N
NIR prediction Native
300°C 500°C Cross-validation 700°C All 300°C 500°C Test set 700°C All
462
Correct classification
Incorrect classification by species
%
Planted
EC
Native
120
119
1
119
99.2
Planted
60
0
60
60
100
Native
120
114
6
114
95.0
Planted
60
4
56
56
93.3
Native
120
111
7
111
94.2
Planted
60
13
48
48
78.3
Native
360
339
21
339
94.2
Planted
180
61
119
119
66.1
Native
21
20
1
20
95.2
Planted
15
0
15
15
100
Native
23
21
2
21
91.3
Planted
13
1
12
12
92.3
Native
25
23
2
23
92.0
Planted
11
1
10
10
90.9
Native
72
62
10
62
86.0
Planted
36
8
28
28
77.7
EV
C
J
P
T
1 6 4 5 10
3 12
25
2 1
8
36 1 2
1 1
1
1 4 3
1
5
5
EC: Clones of Eucalyptus sp. from Cenibra. EV: Clones of Eucalyptus sp. from Vallourec. C: Cedrela sp. J: Jacaranda sp. P: Aspidosperma sp
ACS Paragon Plus Environment