INTRODUCTION
Brazil stands out in terms of energy resources since it has all the primary energy sources, with a favorable scenario for its energy demands in long-term (VENTURA FILHO, 2009). This is due to its large territorial area, almost all located in tropical and rainy regions, which offers excellent conditions for the production and energetic use of biomass in large-scale.
Most of the energy biomass in Brazil is used for the production of charcoal, which makes the country the world’s largest producer and consumer of this product (BARCELLOS, 2007). In this context, Cerrado biome is one of the energy sources of Brazil due to its unique ecological aspect, as well as its peculiar physiognomies (SILVA; VALE; MIGUEL, 2015).
Wood (firewood and wood products) is the oldest source of fuel used to produce energy and it is an essential component in meeting Brazil’s energy demand. This energy source will probably continue to be predominant, with most of the consumption located in the production sectors of charcoal, household, industrial, and agricultural (NASCIMENTO; BIAGGIONI, 2010; BRITO, 2007).
Although the large use of wood for energy purposes, few studies describe the performance characteristics of Brazilian species (CINTRA, 2009). Thus, the knowledge of some physical aspects of wood, such as its basic density, is essential for the assessment of its qualities for producing energy inputs, especially charcoal, since a high wood density implies more mass per unit volume and hence more energy (BARCELLOS, 2007).
In addition, the influence of immediate analysis (ashes, volatile matter, and fixed carbon), higher calorific value, and energy density are very important in order to determine the energy potential of wood (ANDRADE; CARVALHO 1998; QUIRINO et al., 2005). Thus, considering the importance of these variables in the use of wood for energy purposes, the aim of this study was to determine a model for the energy density of Cerrado species.
MATERIAL AND METHODS
Study area
Wood samples were collected from trees in a Cerrado area of 10.15 ha located in the Lajeado State Park, with a total area of 9 thousand hectares of Cerrado (COLEN; SILVA; MARTINS, 2007). The area is located in the central-west region of Tocantins State, Brazil, near the city Palmas.
Methodological procedures
Data from the area inventory were used, which presented 2473 individuals, totaling 80 species. From this total, 34 species were selected for sampling of trunks and branches and submitted to immediate analysis, i.e. ashes (ASH), volatile matter (VM), and fixed carbon (FC), and higher calorific value (HCV) (Table 1). Dry and wet mass data were registered for basic density analysis. The number of individuals of each species cut down was determined from the forest inventory, where approximately 3% of trees with DBH > 5 cm, or at least one individual of each species, were collected with permission of the Environmental Agency of Tocantins State (NATURATINS). Disks of 5 cm thick were collected from each individual cut down in three positions of the trunk (base, middle, and top) and in three diameters of the branches (thick branch, with a diameter higher than 10 cm; medium branch, with a diameter of 3 to 10 cm; and thin branch, with a diameter lower than 3 cm). This material was milled, graded below 60 mesh, stored in properly labeled containers, and oven dried at 103 ± 2 °C.
The analyses were carried out having as reference the following Brazilian standards:
- NBR 8112: Immediate analysis. Determination of the contents of moisture, ashes, volatile matter, and fixed carbon of charcoal;
- NBR 8633: Determination of calorific value. It prescribes the method of determining the higher calorific value of the wood at constant volume in an adiabatic, isothermal, or static calorimetric pump.
Energy density (De) was determined as shown in Equation 1:
Where ρb is the basic density (kg m-3) and HCV is the higher calorific value (GJ ton-1).
Data analysis
Data were analyzed by means of analysis of variance (ANOVA). A simple linear regression was used to verify the relationship between variables. The choice of regression models was based on the coefficient of determination (R2) at 5% probability level, standard error, residual analysis, and significance. The tests were performed by using the package Stats from R (R DEVELOPMENT CORE TEAM, 2015).
The equation best explaining the energy density estimation of trees as a function of basic density followed the traditional methods of linear regression in the sequence of importance according to Draper and Smith (1981). In other words, it was based on the graphical analysis of residuals (%) and behavior of the adjusted model in relation to the real basic density, standard error of estimation (absolute and percentage, i.e. Syx (m3) and Syx (%), respectively), adjusted coefficient of determination (R2 adj.), and Fischer’s F-value.
The proposed model considered was:
Where DE is the average energy density of the tree (Mcal m-3), β0 + β1 are the coefficients to be adjusted, and DBA is the average basic density of the tree (kg m-3).
The data of inventory carried out in the area were used for analyzing the energy density per hectare, totaling 80 species in a sample area of 2.16 ha (MIGUEL et al., 2016). Average values of tree density and the model adjusted for energy density were used.
RESULTS AND DISCUSSION
The Cerrado tree species assessed in this study presented values of VM between 77.4 and 85.5% (mean = 82.88%), ASH between 0 and 1% (mean = 0.27%), and FC between 14 and 23% (mean = 16.84%) (Table 1). The species presenting the highest values for VM were Byrsonima laxiflora Griseb (85.63%), Conarus perrotetti (85.62%), and Qualea parviflora (85.95%) whereas the lowest values were registered for Ferdinandusa elliptica (77.42%) and Tapirira guianensis (77.81%). The highest value of ASH was recorded for Myrcia splendens (0.83%), while the lowest value was found for Bowdichia virgiloides (0.02%). FC ranged from 13.69% (Qualea parviflora) to 22.24% (Ferdinandusa elliptica).
According to Brito and Barrichelo (1978), the immediate analysis of a fuel gives the percentage of material burning in the gaseous (volatile matter) and solid (fixed carbon) states, as well as an indication of the residual material (ashes). Vale, Brasil and Leão (2002), Chaves et al. (2013), and Brito and Barrichelo (1982) recommended for biomass, in general, volatile matter contents between 75 and 85% and fixed carbon between 15 and 25%. Cintra (2009) worked with native forest species in São Paulo State, Brazil and found values of volatile matter between 79.6 and 84.9%. The values found in our study for these variables are in accordance with those found in the literature for wood and their variation is influenced by the species. The average ASH varied according to the species, but with values considered close to those of planted forests in Brazil, indicating acceptable results for energy purposes by the analyzed material. According to Chaves et al. (2013) and Pincelli (2011), woods with a high ash content are not recommended for energy purposes for being a substance that remains in the solid form and it does not combust. Vidaurre et al. (2012) studied the energy properties of Paricá wood (Schizolobium amazonicum) and found an average ash content of 1.3%, which is considered high when compared to the ash content of Eucalypt. Morais, Nascimento and Melo (2005) found an average value of 1.26% for wood of Pinus oocarpa and related this high value to the higher content of inorganic compounds in the wood of this species when compared to others. Barcellos (2007) and Moutinho et al. (2016) found differences in ash content according to the studied species, with values ranging from 0.12 to 1.11%. Vale, Brasil and Leão (2002) reported values between 0.15 and 2.73% for 47 Cerrado species, indicating that the values found in our research were similar to those found in the literature.
According to Chaves et al. (2013), during the burning of plant biomass volatile matter rapidly volatilize, which contributes to a low energy efficiency since fuel residence time decreases due to its rapid volatilization. A strong and inversely proportional relationship is usually found between volatile matter and fixed carbon. On the other hand, fixed carbon is directly related to the calorific value. A higher fixed carbon content in the material implies a longer fuel residence time.
Vale, Brasil and Leão (2002) found an average value of fixed carbon of 20.73% for Cerrado woods, Almeida et al. (2015) reported values between 17 and 18.57%, and Cintra (2009) found values between 14.4 and 19.6%, which is close to those found in our study. Chaves et al. (2013), Silva (2013), and Barcellos (2007) emphasized that fixed carbon content presents a negative relation with volatile matter content.
The higher calorific value was between 18,282 GJ ton-1 for Parkia platycephala and 20,121 GJ ton-1 for Pouteria ramiflora (mean = 19,233 GJ ton-1). Pinheiro, Rendeiro and Pinho (2005) found values for the higher calorific value of vegetal residue ranging from 18,810 to 20,900 GJ ton-1 regardless of particle size and density of biomass or origin wood in the case of sawdust. Brito (1993) points out that the variation in higher calorific value for wood, in general, is between 14,630 and 20,900 GJ ton-1. The average values cited in the literature are in agreement with those found in our study, i.e. between 18,282 and 20,121 GJ ton-1. Couto et al. (2004) observed an HCV value for wood of 18,366 GJ ton-1. Quirino et al. (2005), studying the higher calorific value of forest species, found values between 14,003 and 21,986 GJ ton-1, with an average value of 19,779 GJ ton-1. In addition, Vidaurre et al. (2012) registered an average value of 18,638 GJ ton-1 whereas Vale, Brasil and Leão (2002) observed values ranging from 18,876 to 20,854 GJ ton-1 considering 47 Cerrado species.
Table 1 Immediate analysis of wood of 34 woody species from the Brazilian Cerrado, Tocantins State
Species | VM (%) | ASH (%) | FC (%) |
---|---|---|---|
Aspidosperma subincanum Mart. | 81.52 | 0.30 | 18.18 |
Byrsonima laxiflora Griseb. | 85.63 | 0.09 | 14.28 |
Byrsonima pachyphylla A. Juss | 83.17 | 0.10 | 16.73 |
Byrsonima sericea A.Juss. B | 83.39 | 0.37 | 16.24 |
Bocageopsis multiflora (Mart.) | 81.07 | 0.07 | 18.86 |
Bowdichia virgiloides Kunth. | 85.58 | 0.02 | 14.40 |
Connarus perrottetti (DC.) Planch | 85.62 | 0.06 | 14.32 |
Connarus suberosus Planch | 83.98 | 0.43 | 15.59 |
Dalbergia densiflora Benth. | 85.34 | 0.28 | 14.38 |
Emmotum nitens (Benth.) Miers. | 85.11 | 0.27 | 14.62 |
Ferdinandusa elliptica Pohl. Pl. Bras. | 77.42 | 0.34 | 22.24 |
Erythroxylum daphnites Mart. | 84.71 | 0.20 | 15.09 |
Himatanthus sucuuba (Spruce ex Mull. Arg.) | 82.12 | 0.19 | 17.69 |
Inga alba (Sw.) Willd. | 85.35 | 0.56 | 14.09 |
Licania apetala (E. Meyer) Fritsch. | 82.29 | 0.33 | 17.38 |
Mabea fistulifera Mart. | 83.98 | 0.14 | 15.88 |
Maprounea guianensis Aubl. | 80.86 | 0.18 | 18.97 |
Matayba guianensis Aubl. | 84.30 | 0.32 | 15.13 |
Mezilaurus itauba (Meissn.) Taub. | 83.30 | 0.33 | 16.37 |
Miconia albicans (Swartz) | 82.24 | 0.54 | 17.22 |
Miconia cuspidata Mart. Ex Naudin. | 81.79 | 0.37 | 17.84 |
Myrcia splendens (Sw.) DC. | 81.00 | 0.83 | 18.17 |
Ouratea ovallis (Pohl) Engl. | 82.41 | 0.12 | 17.47 |
Parkia pendula (Willd.) Benth. | 84.44 | 0.04 | 15.52 |
Parkia platycephala Benth. | 82.57 | 0.45 | 16.98 |
Pouteria ramiflora (Mart.) Radlk. | 82.19 | 0.26 | 17.55 |
Protium heptaphyllum Mart. | 81.19 | 0.22 | 18.59 |
Qualea parviflora Mart. | 85.95 | 0.36 | 13.69 |
Sacoglottis guianensisBenth. | 83.73 | 0.17 | 16.10 |
Pouteria ramiflora (Mart.) Radlk. | 82.19 | 0.26 | 17.55 |
Tachigali vulgaris L. G. Silva & H. C. Lima. | 82.68 | 0.31 | 17.00 |
Tapirira guianensis Aubl. | 77.81 | 0.40 | 21.79 |
Virola sebifera Aubl. | 79.70 | 0.10 | 20.20 |
Vochysia gardineri Warm. | 82.60 | 0.11 | 17.29 |
Xylopia aromática (Lam.) Mart. | 82.88 | 0.27 | 16.84 |
VM = volatile matter; ASH = ashes; FC = fixed carbon.
Among the 34 species analyzed, Pouteria ramiflora (15,006 Mcal m-3), Miconia cuspidata (16,177 Mcal m-3), and Emmotum nitens (15,977 Mcal m-3) showed the highest energy density values. On the other hand, the species with the lowest energy density were Vochysia gardineri (6,871 Mcal m-3), Connarus perrottetti (8,380 Mcal m-3), and Byrsonima laxiflora (9,315 Mcal m-3). The average value of energy density for the 34 species was 12,459 Mcal m-3.
According to Vidaurre et al. (2012) and Cunha (1989), the density has no direct relation with the amount of energy contained in the wood. Both studies, as well as the result observed in ours, mentioned no correlation between basic density and higher calorific value. However, these authors emphasize that the higher the density is, the greater the amount of energy contained per unit volume, which stimulates the interest in denser woods for burning. When the relationship between the variables of the immediate analysis (VM, ASH, and FC) and HCV was tested, a correlation was observed between VM and FC (R2 = 0.9927) (Figure 1). The higher calorific value and basic density of the studied species did not present a significant relationship since the high-density species showed a low HCV (Figure 2). Table 2 and Figure 3 show the relationship between energy density and other variables.

Figure 1 Relationship between fixed carbon, ashes, and volatile matter and higher calorific value of Cerrado tree species from Tocantins State, Brazil.
Table 2 Relationship between basic density (ρb), higher calorific value (HCV), ashes (ASH), and volatile matter (VM) of Cerrado tree species from Tocantins State, Brazil.
Relation | R2 | F | p | Residual |
---|---|---|---|---|
Db × HCV | 6.03−05 | 0.004705 | 0.94549 | 0.013738 |
DE × VM | 0.001653 | 0.129144 | 0.720292 | 2.9211 |
DE × FC | 0.006186 | 0.485523 | 0.488003 | 2.9E11 |
DE × ASH | 0.036694642 | 0.08872 | 0.08872 | 2.8111 |
DE × Db | 0.985211 | 5129.459 | 0.735273 | 3.3511 |

Figure 2 Relationship between higher calorific value and basic density of 34 Cerrado tree species from Tocantins State, Brazil.

Figure 3 Relationship between energy density and calorific value, volatile material content (VM), fixed carbon content (FC), ash content (ASH) of 34 Cerrado tree species from Tocantins State, Brazil.
In a study on the energy potential of six-year-old eucalyptus trees, Lima et al. (2007) found 9,168 Mcal m-3 of energy density. Pinheiro, Rendeiro and Pinho (2005) studied the energy density of plant residues (fruit peels) and found values of energy density from 400 to 2,600 Mcal m-3. When analyzing six-year-old eucalyptus trees, Lima (2011) observed a value of energy density of 11,103 Mcal m-3 whereas Brito, Barrichelo and Seixas (1983), found values ranging from 7,843 to 12,038 Mcal m-3 for eucalyptus species. According to Lima (2011), energy density considers the energy contained in a given wood volume. The variation of values in our study was close to the values found in the literature, ranging from 6,871 to 16,177 Mcal m-3, with an average of 12,459 Mcal m-3. Although the average energy density in the literature does not present so much variation, few studies with native species can be found.
Brito and Barrichelo (1977) concluded that in the choice of wood to obtain charcoal with better chemical properties (higher contents of fixed carbon and lower of volatile substances and ashes), those with high lignin contents and high basic density should be sought. This provides an increase in the amount of dry matter in the oven. In a study with 108 forest species, Quirino et al. (2005) observed that the highest calorific value does not match the species with the highest basic density. In addition, Quirino et al. (2012) reported that low-density woods could be used for energy if converted into briquettes.
Fixed carbon is directly related to the calorific value, especially the charcoal, being the variation higher, ranging from 6% in a dry base of wood to up to 90% in charcoal. This may be the reason for the low relationship between the higher calorific value and fixed carbon content found in our study. Another possible explanation may be the presence of extractives. Vale, Sarmento and Almeida (2005) reported that the best use of wood for energy production as heat is linked to a higher specific mass, higher fixed carbon content, and lower volatile matter content. Souza (2010) observed that the parameter density might influence the calorific value of a material as the heat generated by combustion is related to the volume. Thus, wood and denser wood products presented higher calorific value per unit volume.
The coefficient of determination (R2) presented in Table 3 showed that 98% of the energy density variation registered can be explained by the average basic density of the tree, with mean percentage errors of less than 2.1%, as the standard error of the estimate (Syx). The model adjusted for energy density estimation as a function of the basic density is shown in Table 3. The residual distribution did not exceed a variation of 6% either to underestimate or overestimate the energy density of the tree by the basic density of the 34 species (Figure 4).
Table 3 Model adjustment to estimate the energy density of the tree (DE) as a function of its basic density.
Model | R2 adj. | Syx | Syx (%) | F |
---|---|---|---|---|
|
||||
DE = −1.25645 + 19.21955 × DBA | 0.98 | 277.06 | 2.15 | 4034.15 |
DE = average energy density of the tree; β0 + β1 = estimated coefficients; DBA = average basic density of the tree; R2 adj. = coefficient of determination; Syx = standard error of estimation; F = F value from the analysis of variance.

Figure 4 Residual distribution according to the adjusted model to determine the energy density as a function of the basic density of the tree.
The estimated energy density of the area was 1,378,541 Mcal ha-1. The proposed model for calculating the energy density as a function of the basic density facilitates the energy surveys in areas to be explored since only the values of basic density need to be obtained, without the need of calorific value analysis. The knowledge of the energy potential of the analyzed native species allows advancing in the studies for possible improvements in the quality of wood species directed to a sustainable forest energy production.