Canonical correlation analysis of the characteristics of charcoal from Qualea parviflora Mart.

Análise de correlações canônicas das características do carvão vegetal de Qualea parviflora Mart.

Canonical correlation analysis of the characteristics of charcoal from Qualea parviflora Mart.

Análise de correlações canônicas das características do carvão vegetal de Qualea parviflora Mart.

Thiago de Paula ProtásioI, ^*; Rosalvo Maciel Guimarães Neto^II; João de Deus Pereira de Santana^II; José Benedito Guimarães Júnior^II; Paulo Fernando Trugilho^I

^IUniversidade Federal de Lavras - Lavras, Minas Gerais, Brasil

^IIUniversidade Federal do Piauí - Bom Jesus, Piauí, Brasil

ABSTRACT

This study aimed to examine the relationships between the characteristics of charcoal from Qualea parviflora Mart. using canonical correlation analysis. Five trees were analyzed in such way that 5-cm thick discs were removed from each tree at the base, DBH (1.30 m), middle and top sections. The wood was carbonized in a muffle furnace at a heating rate of 1.67 °C min-1. A canonical correlation analysis was conducted to investigate the relationships between the group formed by fixed carbon, volatile matter, ash, elemental carbon, hydrogen, nitrogen, sulfur and oxygen levels and a second group formed by the gravimetric yield, higher heating value and relative bulk density of the charcoal. A tendency was noted for high levels of fixed carbon and elemental carbon to be associated to low levels of volatile matter, ash and oxygen and to low gravimetric yield. Fixed carbon and elemental carbon levels had a positive relation to higher heating value and to relative bulk density, whereas volatile matter, ash and oxygen levels had a negative relation to such characteristics. The higher the gravimetric yield from carbonization, the higher the volatile matter, ash and oxygen levels will be in the resulting charcoal.

Keywords: charcoal, chemical composition, heating value, multivariate analysis.

RESUMO

Neste trabalho, objetivou-se utilizar a análise de correlação canônica para avaliar as relações existentes entre as características do carvão vegetal de Qualea parviflora Mart. Foram avaliadas cinco árvores das quais foram retirados discos de 5 cm de espessura na base, DAP (1,30 m), meio e topo. A madeira foi carbonizada em uma mufla, considerando a taxa de aquecimento de 1,67 ºC min^-1. Foi realizada a análise de correlação canônica entre o grupo formado pelos teores de carbono fixo, materiais voláteis, cinzas, carbono elementar, hidrogênio, nitrogênio, enxofre e oxigênio com um segundo grupo formado pelo rendimento gravimétrico, poder calorífico superior e a densidade relativa aparente do carvão vegetal. Observou-se uma tendência de elevados teores de carbono fixo e carbono elementar associados ao carvão vegetal de baixos teores de materiais voláteis, cinzas e oxigênio e de baixo rendimento gravimétrico. Os teores de carbono fixo e de carbono elementar possuem relação positiva com o poder calorífico superior do carvão vegetal e com a densidade relativa aparente. Já, os teores de materiais voláteis, cinzas e oxigênio possuem relação negativa com essas características. Quanto maior o rendimento gravimétrico da carbonização, maiores serão os teores de materiais voláteis, cinzas e oxigênio presentes no carvão vegetal.

Palavras-chave: carvão vegetal, composição química, poder calorífico, análise multivariada.

1 INTRODUCTION

Species Qualea parviflora Mart. has high basic density (690 kg m-³) and provides high dry matter yield of wood (524 kg ha^-1) (VALE et al., 2002). Popularly known as Pau-terra-de-folha-miúda or simply Pau-terra, the species belongs to family Vochysiaceae and is most commonly found in the transition zone between the Cerrado and Caatinga biomes. Qualea parviflora Mart. naturally occurs in Bahia, Minas Gerais, Mato Grosso do Sul, Piauí and São Paulo states, in cerrado and campos cerrados savannas (LORENZI, 2002).

Around 35% to 50% of the charcoal consumed by Brazilian steel companies is sourced from residues from management and exploration of natural forests, particularly of the Cerrado biome, since homogeneous forest plantations have proved unable to meet the entire market demand (ASSOCIAÇÃO BRASILEIRA DE PRODUTORES DE FLORESTAS PLANTADAS - ABRAF, 2012; BRASIL, 2011; CALAIS, 2009). Therefore, research studies become increasingly critical not only to evaluate the quality of charcoal from native species with high potential for use in energy and steel companies, as is the case of Qualea parviflora Mart. (PROTÁSIO et al., 2011a; VALE et al., 2002), but also to investigate the relationships or correlations between the characteristics of interest.

In order to evaluate charcoal intended for use in energy and/or steel production, it is necessary to quantify characteristics such as carbonization yield, heating value, bulk density and immediate chemical composition, and learning the relations between such characteristics can be very useful for controlling the pyrolysis process, optimizing charcoal quality and ensuring genetic improvement of the species of interest.

With that in mind, multivariate techniques such as canonical correlation are a viable alternative not only for analysis and interpretation of data in studies exploring wood and charcoal quality, but also for identifying quality coefficients and/or indices (TRUGILHO et al., 1997a).

Canonical correlation analysis enables verifying the existing linear correlations between two groups or sets of variables (X and Y) and is intended to maximize the correlations between such variables. The analysis consists of obtaining a pair of latent variables, known as canonical statistical variables, that are linear combinations of the variables of two vectors (X and Y), and the information contained in the parameters should be focused on the correlation between these new variables. Each pair of canonical statistical variables is termed canonical function (ABREU; VETTER, 1978; CRUZ; REGAZZI, 1994; FERREIRA, 2008; HAIR JÚNIOR et al., 2005, 2009; MINGOTI, 2005).

The above technique is used in exploratory studies, as in some situations researchers may have a large set of variables and yet be interested only in studying a few linear combinations of the variables within such set (TRUGILHO et al., 1997b, 2003).

Ferreira (2008) argued that knowledge of the interrelation between distinct variable vectors can allow researchers to identify possible dependence structures between the different sets of variables, which would otherwise not be easily identifiable by simple linear correlation analysis. And yet canonical correlation analysis is the multivariate statistical technique least explored by users (MINGOTI, 2005).

Given the above, the aim of this study was to use canonical correlation analysis to evaluate the relationships between the characteristics of charcoal from Qualea parviflora Mart.

2 MATERIAL AND METHODS

2.1 Carbonization procedure and analysis of the resulting charcoal

The wood material was obtained from a Cerrado remnant sitting between the municipalities of Bom Jesus and Currais, in Piauí state. Five trees were randomly sampled, and from each tree 5-cm thick discs were removed from the base, DBH (1.30 m), middle and top sections, with a total of five replicates being run per lengthwise position, noting that 5.38 m was defined as the average total height. The average DBH outside bark was 51.0 cm, against 40.4 cm inside bark. For purposes of sampling, height was defined as including only tree trunks.

The wood material was carbonized in an electric muffle furnace at a heating rate of 1.67 °C min-1. The initial temperature was 100 °C and the final temperature was 450 °C, remaining stable for 30 minutes (ASSIS et al., 2012; NEVES et al., 2011; TRUGILHO et al., 1997a, 1997b, 2005). The total carbonization time was four hours. In each trial around 800 gr. of wood was used, which had been previously oven dried at 103 ± 2 ºC. After each carbonization procedure, the gravimetric charcoal yield was computed.

The relative bulk density of the resulting charcoal was determined by the hydrostatic method, considering water as the medium and as described in NBR 11941 (ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS - ABNT, 2003).

The charcoal was subjected to immediate chemical analysis for quantification of fixed carbon, volatile matter and ash levels, according to the guidelines set out in standard NBR 8112 (ABNT, 1983).

For the quantification of higher heating value, an IKA C-200 digital calorimeter was used following standard NBR 8633 (ABNT, 1984).

For the elemental analysis, the samples were crushed and sieved, and the fraction retained by a 270-mesh sieve was used, as described by Paula et al. (2011) and Protásio et al. (2011b). The quantification of carbon, hydrogen, nitrogen and sulfur levels, relative to dry charcoal weight, was performed using a universal Elementar analyzer (model Vario Micro Cube). This analyzer uses helium as carrier gas and oxygen as ignition gas. The 2-mg samples were placed in tin capsules and thoroughly incinerated at 1,200 °C. Oxygen content was obtained by the difference (Eq. 1), following the guidelines provided by Bech et al. (2009).

where: O is oxygen content (%), C is carbon content (%), H is hydrogen content (%), N is nitrogen content (%), S is sulfur content (%) and Cz is ash content (%).

2.2 Canonical correlation analysis

A canonical correlation analysis was performed to establish correlations between the group formed by fixed carbon, volatile matter, ash, elemental carbon, hydrogen, nitrogen, sulfur and oxygen levels and a second group formed by the gravimetric yield, higher heating value and relative bulk density of the charcoal from Qualea parviflora Mart. The first group represents the dependent variables (Y) whereas the second group represents the independent variables (X).

It was possible to determine three canonical functions, or three pairs of canonical statistical variables. Estimations were made of the canonical loadings, that is, the correlations between the original variables and their respective canonical statistical variables, the canonical cross loadings which represent the correlation between an original variable of a given group and the canonical statistical variable of the opposite group, and also the canonical weights (similar to beta coefficients of a multiple regression).

To correct any possible effects of the lengthwise sampling positions on the characteristics of the charcoal being evaluated, the canonical correlations were computed based on the residuals of variance analyzes performed for each variable.

Preliminarily to the canonical correlation analysis, simple linear correlations were also estimated between the variables so as to test for multicollinearity between variables from different groups.

The Hotelling's multivariate test (approximation of the F-distribution) was used to test the significance of the canonical roots jointly. Also, the generalized Shapiro-Wilk test was used to test for multivariate normality in the data used, since the significance test of canonical correlation is only valid when vectors X and Y are multivariate normal (FERREIRA, 2008; JOHNSON; WICHERN, 1992; MINGOTI, 2005).

The amount of explained variance, that is, the percentage of variance in the dependent canonical statistical variable that may be explained by the independent canonical statistical variable, and vice versa, was determined by squaring the canonical correlation (canonical R²).

The amount of explained variance between the observed dependent and independent variables and their respective canonical statistical variables was determined by squaring the canonical loadings. To calculate the amount of variance explained by the canonical statistical variable, an arithmetic average was taken of the canonical loadings squared. The same procedure was used for the canonical cross loadings in order to estimate the shared variance between the observed dependent or independent variable and the opposite canonical statistical variable.

The redundancy index was also computed, it being the average squared canonical loading times the canonical R². The redundancy index expresses the amount of variance in a canonical statistical variable (dependent or independent) explained by the opposite canonical statistical variable in a canonical function. A high redundancy index entails a high canonical correlation and a high degree of shared variance explained by the canonical statistical variable. The redundancy index should be considered in trying to overcome the bias and uncertainty inherent in the use of canonical roots (squared canonical correlations) as a measure of shared variance (HAIR JÚNIOR et al., 2009).

The statistical analysis procedures used were grounded in studies performed by Cruz and Regazzi (1994), Ferreira (2008), Hair Júnior et al. (2005, 2009), Mingoti (2005), Protásio et al. (2012) and Trugilho et al. (1997b, 2003).

All analysis procedures were performed using software R version 2.11.0 (R DEVELOPMENT CORE TEAM, 2008), packages CCA (GONZÁLEZ; DÉJEAN, 2009), CCP (MENZEL, 2009) and yacca (BUTTS, 2009).

3 RESULTS AND DISCUSSION

Table 1 provides values of canonical correlations, canonical R² and Hotelling's multivariate significance test with approximation of the F-distribution.

Canonical functions 1 and 2 were found to be significant and have a high canonical R². Results reveal that 97.32% and 79.42% of the variance between the groups of charcoal characteristics is explained by canonical functions 1 and 2 respectively. This is indicative of interdependence between the groups of variables and can be of great help in modeling, in understanding the pyrolysis process and in ensuring the genetic improvement of the species for bioenergy applications, particularly in the steel industry.

Figure 1 provides scores of significant canonical functions by Hotelling's test. A significant interdependence was found between the variables of each group, that is, the characteristics of charcoal encompassed by groups X and Y not only are not independent but have a high correlation.

Table 2 provides canonical weights, canonical loadings and canonical cross loadings for canonical function 1.

An analysis and interpretation of canonical weights entails examining the signal and its magnitude. Variables with relatively higher canonical weights contribute more toward the statistical variables and vice versa. As regards the dependent group, it was noted that fixed carbon, volatile matter and ash levels have the highest canonical weights. And as for the independent group, gravimetric charcoal yield can be considered the most important variable. It should be noted, however, that canonical weights are subject to high variability (instability) and are sensitive to multicollinearity (HAIR JÚNIOR et al., 2009).

In the dependent group, the highest canonical loadings, in module, relate to fixed carbon and volatile matter levels. These variables can thus be considered the most important in the dependent canonical variable set. Hydrogen and sulfur levels had statistically negligible canonical loadings. This probably occurred because the low levels of these chemical components in charcoal from Qualea parviflora are indicative of their low correlations with other parameters being evaluated.

In the dependent group, the highest canonical loadings, in module, relate to fixed carbon and volatile matter levels. These variables can thus be considered the most important in the dependent canonical variable set. Hydrogen and sulfur levels had statistically negligible canonical loadings. This probably occurred because the low levels of these chemical components in charcoal from Qualea parviflora are indicative of their low correlations with other parameters being evaluated.

In the independent canonical variables group, only relative bulk density had low canonical loadings. The same trend was observed for canonical cross loadings.

An analysis of the canonical loadings revealed a tendency for high fixed carbon levels to be associated to low volatile matter, ash and oxygen levels and to high elemental carbon level. The existing relationships between fixed carbon and nitrogen levels were not sufficiently explained by the canonical correlation analysis, and this could be attributed to the low canonical loading observed for nitrogen level.

Again in the independent group, a positive correlation was found between relative bulk density and higher heating value. Gravimetric charcoal yield, in its turn, had a negative correlation with relative bulk density and higher heating value.

Vale et al. (2001) found no correlation between gravimetric charcoal yield and bulk density using simple linear correlation analysis, which reinforces the importance of using canonical correlation analysis for understanding existing correlations between charcoal characteristics.

Trugilho and Silva (2001) found no positive correlation between relative bulk density and higher heating value in charcoal from Jatobá, unlike the findings in this study. These authors found a positive linear correlation between charcoal yield and bulk density. However, similarly to the findings from the canonical correlation analysis, Trugilho and Silva (2001) found a negative correlation between gravimetric charcoal yield and higher heating value.

An analysis of the canonical cross loadings revealed that higher fixed carbon and elemental carbon levels correlate positively with higher heating value and relative bulk density but negatively with gravimetric yield.

It was further noted that the higher the gravimetric charcoal yield the higher the volatile matter, ash and oxygen levels in charcoal will be. Likewise, the higher the volatile matter, ash and oxygen levels, the lower the values of relative bulk density and higher heating value will be and vice versa.

Therefore, when evaluating charcoal for use in the steel industry, a high fixed carbon content is desirable on account of the higher percentage of reducing carbon that can be used in the chemical reactions. Additionally, the higher the fixed carbon content and the lower the volatile matter content, the greater the higher heating value will be for the resulting charcoal, as has been found both in this study and in studies performed by Reis et al. (2012) and by Trugilho and Silva (2001), and the greater the thermal stability of the fuel and its residence time in the furnace will be for thorough burning (BRAND, 2010; PROTÁSIO et al., 2013; VALE et al., 2001).

Also, it is important to consider the influence maximum carbonization temperature has on charcoal characteristics, such as immediate chemical composition and yields evaluated (TRUGILHO; SILVA, 2001; VALENTE et al., 1993). In this study, however, the Qualea parviflora charcoal samples were obtained under the same condition of pyrolysis, that is, the observed correlations allow understanding the typical interrelations between the characteristics of this bioreductant without the effect of temperature.

Determining the ash content is critical when evaluating the bioenergetic potential and use potential of charcoal in the steel industry, inasmuch as high ash contents contribute to reducing the heating value, given that minerals do not participate in the combustion process. This statement corroborates the results found both in this work and in studies performed by Brand (2010), Paula et al. (2011), Protásio et al. (2011b) and Reis et al. (2012). Additionally, high ash contents in charcoal may cause impurities to accumulate in the center of solidified metal pieces, leading to variations in the physical, chemical and mechanical properties of the resulting pig iron or iron alloys.

With respect to the influence of oxygen content plus minerals on the heating value of biomass fuels, Protásio et al. (2011b) reported a negative correlation between such variables, which supports the findings for Qualea parviflora charcoal. The relationships found between levels of fixed carbon, volatile matter, higher heating value and gravimetric yield are consistent with results found in literature (TRUGILHO; SILVA, 2001; VALE et al., 2001).

Table 3 provides canonical weights, canonical loadings and canonical cross loadings for canonical function 2. In the dependent group, elemental chemical components had the highest canonical weights, whereas in the independent group, gravimetric yield had the highest canonical weights.

In the independent group, relative bulk density had the highest canonical loading in module, whereas in the dependent group, nitrogen content can be considered the most important variable. Analysis results for canonical function 2 revealed a tendency for nitrogen content to increase with decreasing relative bulk density and vice versa.

Redundancy measures of shared variance for canonical functions 1 and 2 are provided in Table 4.

In canonical function 2, the variables relative bulk density and nitrogen content had the highest percentage of explained variance within and between groups.

Results demonstrate that fixed carbon, volatile matter, higher heating value, gravimetric yield and relative bulk density are the most significant variables in the canonical correlation analysis.

Figure 2 provides redundancy indices for significant canonical functions according to Hotelling's test.

In canonical function 1, 50% of variance in the dependent variables was explained by the independent canonical statistical variable. The dependent canonical statistical variable (Y), on the other hand, explained around 25% of variance in the independent variables (X). In canonical function 2, the independent group explained around 26% of variance in the dependent variables.

4 CONCLUSIONS

The canonical correlation analysis proved effective in determining the existing relationships between the characteristics of charcoal from Qualea parviflora Mart., revealing an interdependence between the groups of interest.

A tendency was noted for high fixed carbon and elemental carbon contents to be associated to charcoal with low volatile matter, ash and oxygen contents and with low gravimetric yield.

Fixed carbon and elemental carbon levels correlated positively with higher heating value of charcoal and with relative bulk density. Volatile matter, ash and oxygen levels, on the other hand, correlated negatively with such characteristics.

The higher the gravimetric yield from carbonization, the higher the volatile matter, ash and oxygen contents present in the resulting charcoal.

The authors hope these results will be valuable in modeling, in understanding the pyrolysis process and in improving the bioenergetic qualities of charcoal from Qualea parviflora Mart.

5 REFERENCES

Recebido: 23 de setembro de 2011

Aceito: 22 de julho de 2013

• ABREU, M. A.; VETTER, D. A análise de relação entre conjuntos de variáveis na matriz geográfica: correlação canônica. In: FAISSOL, S. (Ed.). Tendências atuais na geografia urbana/regional: teorização e quantificação Rio de Janeiro: IBGE, 1978. p. 133-144.
• ASSIS, M. R.; PROTÁSIO, T. P.; ASSIS, C. O.; TRUGILHO, P. F.; SANTANA, W. M. S. Qualidade e rendimentos do carvão vegetal de um clone híbrido de Eucalyptus grandis x Eucalyptus urophylla Pesquisa Florestal Brasileira, Colombo, v. 32, n. 71, p. 291-302, 2012.
• ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 11941: madeira: determinação da densidade básica Rio de Janeiro, 2003. 6 p.
• ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 8112: análise imediata: material volátil, cinzas, carbono fixo Rio de Janeiro, 1983. 6 p.
• ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 8633: carvão vegetal: determinação do poder calorífico Rio de Janeiro, 1984. 13 p.
• ASSOCIAÇÃO BRASILEIRA DE PRODUTORES DE FLORESTAS PLANTADAS. Anuário estatístico da ABRAF 2012: ano base 2011 Brasília, 2012. 150 p.
• BECH, N.; JENSEN, P. A.; DAM-JOHANSEN, K. Determining the elemental composition of fuels by bomb calorimetry and the inverse correlation of HHV with elemental composition. Biomass and Bioenergy, Oxford, v. 33, n. 3, p. 534-537, 2009.
• BRAND, M. A. Energia de biomassa florestal Rio de Janeiro: Interciência, 2010. 131 p.
• BRASIL. Ministério do Meio Ambiente. Plano de ação para prevenção e controle do desmatamento e das queimadas: cerrado Brasília, 2011. 200 p.
• BUTTS, C. T. Yet another canonical correlation analysis package R package version 2.11.0. Disponível em: <http://r-forge.r-project.org>. Acesso em: 24 nov. 2009.
• CALAIS, D. Florestas energéticas no Brasil: demanda e disponibilidade Belo Horizonte: Associação Mineira de Silvicultura, 2009. 23 p. Disponível em: <http://www.silviminas.com.br/arquivo/publicacoes.aspx?ano=23>. Acesso em: 12 jul. 2011.
• CRUZ, C. D.; REGAZZI, A. J. Modelos biométricos aplicados ao melhoramento genético Viçosa, MG: UFV, 1994. 390 p.
• FERREIRA, D. F. Estatística multivariada Lavras: UFLA, 2008. 662 p.
• GONZÁLEZ, I.; DÉJEAN, S. T. Canonical correlation analysis R package version 2.11.1. Disponível em: <http://r-forge.r-project.org>. Acesso em: 24 nov. 2009.
• HAIR JÚNIOR, J. F.; ANDERSON, R. E.; TATHAM, R. L.; BLACK, W. C. Análise multivariada de dados Porto Alegre: Bookman, 2005. 593 p.
• HAIR JÚNIOR, J. F.; BLACK, W. C.; BABIN, B. J.; ANDERSON, R. E.; TATHAM, R. L. Análise multivariada de dados Porto Alegre: Bookman, 2009. 688 p.
• JOHNSON, R. A.; WICHERN, D. W. Applied multivariate statistical analysis 3rd ed. New Jersey: Prentice-Hall, 1992. 642 p.
• LORENZI, H. Árvores brasileiras: manual de identificação e cultivo de plantas arbóreas do Brasil 2. ed. Nova Odessa: Plantarum, 2002. v. 2, 368 p.
• MENZEL, U. Significance tests for Canonical Correlation Analysis (CCA) R package version 2.11.0. Disponível em: <http://r-forge.r-project.org>. Acesso em: 24 nov. 2009.
• MINGOTI, S. A. Análise de dados através de métodos de estatística multivariada: uma abordagem aplicada Belo Horizonte: UFMG, 2005. 297 p.
• NEVES, T. A.; PROTÁSIO, T. P.; COUTO, A. M.; TRUGILHO, P. F.; SILVA, V. O.; VIEIRA, C. M. M. Avaliação de clones de Eucalyptus em diferentes locais visando à produção de carvão vegetal. Pesquisa Florestal Brasileira, Colombo, v. 31, n. 68, p. 319-330, 2011.
• PAULA, L. E. R.; TRUGILHO, P. F.; NAPOLI, A.; BIANCHI, M. L. Characterization of residues from plant biomass for use in energy generation. Cerne, Lavras, v. 17, p. 237-246, 2011.
• PROTÁSIO, T. P.; BUFALINO, L.; TONOLI, G. H. D.; COUTO, A. M.; TRUGILHO, P. F.; GUIMARÃES JÚNIOR, M. Relação entre o poder calorífico superior e os componentes elementares e minerais da biomassa vegetal. Pesquisa Florestal Brasileira, Colombo, v. 31, n. 66, p. 122-133, 2011a.
• PROTÁSIO, T. P.; MELO, I. C. N. A.; GUIMARÃES JUNIOR, M.; MENDES, R. F.; TRUGILHO, P. F. Thermal decomposition of torrefied and carbonized briquettes of residues from coffee grain processing. Ciência e Agrotecnologia, Lavras, v. 37, n. 3, p. 221-228, maio/jun. 2013.
• PROTÁSIO, T. P.; SANTANA, J. D. P.; GUIMARÃES NETO, R. M.; GUIMARÃES JÚNIOR, J. B.; TRUGILHO, P. F.; RIBEIRO, I. B. Avaliação da qualidade do carvão vegetal de Qualea parviflora Pesquisa Florestal Brasileira, Colombo, v. 31, n. 68, p. 295-307, 2011b.
• PROTÁSIO, T. P.; TONOLI, G. H. D.; GUIMARÃES JÚNIOR, M.; BUFALINO, L.; COUTO, A. M.; TRUGILHO, P. F. Correlações canônicas entre as características químicas e energéticas de resíduos lignocelulósicos. Cerne, Lavras, v. 18, n. 3, p. 433-439, 2012.
• R DEVELOPMENT CORE TEAM. R: a language and environment for statistical computing Vienna: R Foundation for Statistical Computing, 2008. Disponível em: <http://www.R-project.org>. Acesso em: 10 dez. 2010.
• REIS, A. A.; PROTÁSIO, T. P.; MELO, I. C. N. A.; TRUGILHO, P. F.; CARNEIRO, A. C. Composição da madeira e do carvão vegetal de Eucalyptus urophylla em diferentes locais de plantio. Pesquisa Florestal Brasileira, Colombo, v. 32, n. 71, p. 277-290, 2012.
• TRUGILHO, P. F.; LIMA, J. T.; MORI, F. A. Correlação canônica das características químicas e físicas da madeira de clones de Eucalyptus grandis e Eucalyptus urophylla Cerne, Lavras, v. 9, n. 1, p. 66-80, 2003.
• TRUGILHO, P. F.; REGAZZI, A. J.; VITAL, B. R.; GOMIDE, J. L. Aplicação de algumas técnicas multivariadas na avaliação da qualidade da madeira de Eucalyptus e seleção de genótipos superiores para produção de carvão vegetal. Revista Árvore, Viçosa, v. 21, n. 1, p. 113-130, 1997a.
• TRUGILHO, P. F.; SILVA, D. A. Influência da temperatura final de carbonização nas características físicas e químicas do carvão vegetal de jatobá (Himenea courbaril L.). Scientia Agraria, Curitiba, v. 2, n. 1/2, p. 45-53, 2001.
• TRUGILHO, P. F.; SILVA, J. R. M.; MORI, F. A.; LIMA, J. T.; MENDES, L. M.; MENDES, L. F. B. Rendimentos e características do carvão vegetal em função da posição radial de amostragem em clones de Eucalyptus Cerne, Lavras, v. 11, n. 2, p. 178-186, 2005.
• TRUGILHO, P. F.; VITAL, B. R.; REGAZZI, A. J.; GOMIDE, J. L. Aplicação da análise de correlação canônica na identificação de índices de qualidade da madeira de eucalipto para a produção de carvão vegetal. Revista Árvore, Viçosa, v. 21, n. 2, p. 259-267, 1997b.
• VALE, A. T.; COSTA, A. F.; GONÇALEZ, J. C.; NOGUEIRA, M. Relações entre a densidade básica da madeira, o rendimento e a qualidade do carvão vegetal de espécies do cerrado. Revista Árvore, Viçosa, v. 25, n. 8, p. 89-95, 2001.
• VALE, A. T.; FIEDLER, N. C.; SILVA, G. F. Avaliação energética da biomassa do cerrado em função do diâmetro das árvores. Ciência Florestal, Santa Maria, v. 12, n. 2, p. 115-126, 2002.
• VALENTE, O. F.; BARBOSA, J. A.; VITAL, B. R.; LELLES, J. G. Rendimento e qualidade do carvão produzido da madeira de algaroba (Prosopis juliflora DC.). Revista Árvore, Viçosa, v. 17, n. 1, p. 38-44, 1993.

Publication Dates

• Publication in this collection
  11 June 2014
• Date of issue
  Mar 2014

History

• Received
  23 Sept 2011
• Accepted
  22 July 2013