Specification of the properties of the energy wood pieces (EN 14961-1:2011) .
Pruning of fruit trees produces a great quantity of biomass each year that can be used for energy production. For this purpose, it is necessary to carry out an energy characterization of these pruned wastes, where the determination of heating value is significant. This value is usually measured by an adiabatic or isoperibolic calorimeter, which causes high economic costs and wastes time. The present study is focused on the development of indirect models for heating value prediction of biomass from orange trees Citrus × sinensis Osbeck, almond trees Prunus dulcis (Mill) D.A. Webb, and olive trees Olea europaea L. from an elemental analysis in order to reduce the time of determination as well as the economic costs. Residual biomass was classified and characterized according to CEN regulations such as received, without drying. Also, moisture content wet basis, bark ratio, density, heating value, and elemental composition (carbon, hydrogen, nitrogen, and sulfur) were measured. The influence of these variables on the heating value was analyzed. Finally, mathematical models were developed to predict this value for this studied species. These models showed coefficients of determination between 0.83 and 0.97, being suitable for industrial use.
- economical studies
- wood residues
- higher heating value
A number of researchers have published mathematical models to predict the higher heating value of different biomass materials from the concentration of the main elements present, such as percentage of carbon, percentage of hydrogen, and percentage of nitrogen together with others [1–6]. On the other hand, other models have used proximate analysis [7–10] or structural analysis [9–11]. Indirect calculation of the higher heating value by means of these types of models is justified by the expensive cost of the use of calorimeters [3, 12]. The aim of this chapter is to compare the resources used in direct heat value determination with that used in indirect calculation from elemental analysis by means of prediction models in three common lignocellulosic materials coming from pruning Mediterranean fruit trees. In this work, predictive models of the heating value of biomass from pruning of
The heat value is an invariable parameter for a material with constant composition, defined by the empirical form CH
The capacity of retaining water in the biomass, caused by its porosity, must be considered in the measurement of the calorific value. The moisture content in the material changes its molecular formula, and therefore the gravimetric percentages of C, H, O, and N. For this reason, standards to determine the calorific value for a particular material, such as UNE 164001:2005 EX , require the material to be obtained in the anhydrous state, or with a known moisture content.
The problem in industrial applications is that the biomass materials received for the combustion in boilers show variation in their moisture content, they are sometimes mixed with leaves and even with other materials. Under these conditions, the heat properties are not constant. In order to measure the calorific value instantly, regression models are analyzed in this paper from the percentage contents of C, H, and N of the materials as received. Currently, there are devices capable of measuring these elements in 5 min [5, 6].
When there is variability in the composition, and uncertainty of the type of materials or proportion of the mixture thereof, the mathematical models for the indirect determination of the calorific value are only applicable in the scope where they were developed. This uncertainty is quite common in industry. Thus researchers, such as Francis and Lloyd , Ebeling and Jenkins , and Kathiravale et al. , have provided models for different types of mixtures. However, none of them showed an economical study to justify the use of the model. In this chapter, the development of models specific for the studied material is shown; in addition, an analysis was carried out to compare the cost of direct determination with an indirect mathematical model from elemental analysis.
In industrial facilities, it is very difficult to find biomass received for combustion without some moisture. Because of processes increase production costs, they are rarely used in the production of energy wood. On the other hand, open air drying rarely decreases the moisture content below 20% [17, 18]. Moreover, it is normal that industrial boilers do not work with a well-defined type of material but with mixtures of different types of biomass. These reasons make the composition of the biomass used in industry variable that directly affects the expected calorific value. So, calorific determination before the introduction of the materials in the boiler is useful to understand the energy performance of the combustion. If this direct determination by the bomb calorimeter is more expensive than the indirect determination from their composition, developing predictive models is fully justified. This is studied in this work.
2. Materials and methods
2.1. Vegetal material
The species studied in this work were
The studied orange variety was “Valencia Late”. This is one of the latest varieties most cultivated in Spain. The harvesting begins from March, and it is usually pruned in summer after harvesting. Its main use is fresh consumption due to its high juice content. It also has high chances of industrialization .
The almond variety studied was “Blaquerna”. It is a self-fertile cultivar that is currently replacing other varieties of lower production. Almond pruning can be annual or biennial according to the development of the tree . During the first 4 years, winter pruning is carried out focused on defining the architecture of the vegetation. From the fourth production year, pruning is performed to remove unproductive branches and to improve the fruit quality .
The studied olive tree variety was “Villalonga”, which is the variety most widely used in Valencia, with a total cultivation area of 23,550 ha. Its main use is in the manufacture of oil and it is also used for fresh consumption .
To define these raw materials as biomass for biofuel, the standard UNE EN 14961-1:2011 was used. According to this standard, the analyzed samples in this work were classified by their origin and sources of solid biofuels: 1. Wood biomass → 1.1. Wood biomass from forest or plantation → 1.1.4. Residues from cuttings → 126.96.36.199. Fresh/green, hardwood (including leaves). Therefore, it would be “cutting hardwood waste” (188.8.131.52). Following the mentioned standard, the specifications of the properties of the studied biofuel are defined in Table 1.
Cutting hardwood waste (184.108.40.206)
|Wood biomass (1.1)
|Marketed form||Wood logs, firewood|
|Moisture content, M (% as it is received) CEN/TS 15149-1, CEN/TS 15149-2|
M 45 (≤45%)
M 35 (≤35%)
M 40 (≤40%)
|Volume or weight, m3, kg, loose or piled as received||15–20 kg of each specied.|
|Volume ratio of split logs||Trunk without cuttings, whole branches|
|Cut surfacea||Surface smooth and regular cut|
|Rust and rot||None of the samples has mold|
2.2. Sample preparation
The branches of each species, obtained from pruning, were divided into six size classes (0–1, 1–2, 2–3, 3–4, 4–5, and >5 cm). Five samples of each class were taken for analysis; therefore, it resulted in 30 samples per species. For sample preparation, the methods defined by the UNE-EN 14780:2012  were followed. The main purpose of sample preparation was to reduce the size of the branches in test portions with the same initial composition, being representative of the original sample. Initial wet basis moisture content average in wood was about 42.24%. In Figure 1, the process of the preparation of the sample is shown. First step was to separate the leaves and wood of the 30 branches of each species arrived at the laboratory (Figure 2a–c). These leaves were crushed with hammer mill and stored in airtight jars with identification labels. On the other hand, the wood was milled until their sizes were lower than 3 mm. Special care was taken to avoid loss of fine particles and moisture during milling and other operations. Average wet basis moisture content in tested particles was 29.85%. The devices used for the sample preparation were as follows:
Manual saw. In order to prevent the moisture loss at the border, cuttings of the central part was used for the analytical determinations. The pieces obtained using this device are shown in Figure 2(d).
Hammer mill crusher of stainless steel, equipped with a 3 mm screen.
Once the samples were crushed, leaves and wood were mixed in defined proportions in each group of size class. The different proportions analyzed were from 10% weight leaves/weight wood in the first sample, 20% weight leaves/weight wood in the first sample in the second, etc. up to 50% of mixture leaves/wood to the fifth sample. So, the effect of leaves on the calorific value was measured. This is very important in this study because pruned material used in boilers has always got a high percentage of leaves, and this influences the moisture content and composition. In many publications, the calorific value calculated referred to that of completely dry and bare material, but this condition is far from the actual applications. In this work, the analysis was focused on samples with variable elemental composition which is usually obtained due to variation in the percentage of leaf and bark content (caused by different size class) and the moisture content (obtained after outdoor drying). For this reason, neither the wood nor the leaf fractions were dried rigorously, but they only experience the natural loss of moisture content during the transport and storage prior to analysis.
2.3. Measurement process
Higher heating value, wet basis moisture content, and elemental composition (C, H, and N) of each sample were measured. For this, standards shown in Table 2 were used.
|Reference of standard||Standard|
|CEN/TS 14779||Solid biofuels—sampling—methods for preparing sampling plans and sampling certificates|
|CEN/TS 14780||Solid biofuels—methods for sample preparation|
|EN 14774-3||Solid biofuels—determination of moisture content—Stufe drying method part 3 moisture content analysis for overall sample analysis|
|EN 14918||Solid biofuels—determination of calorific value|
|CEN/TS 15104||Solid biofuels—determination of total carbon, hydrogen, and nitrogen—instrumental methods|
|CEN/TS 15290||Solid biofuels—determination of major elements|
|EN 14961-1||Solid biofuels—fuel specifications and classes—part 1: general requirements|
|Variable||Species||Average||Standard deviation||Coef. skewness||Coef. kurtosis||Max.||Min.|
|Moisture content %||Almond||19.009||4.870||0.967||−0.392||30.314||12.291|
|C||H||N||S||% leaves||% bark||% moisture content||HHV (kJ/kg)|
||C||H||N||S||% leaves||% bark||% moisture content||HHV (kJ/kg)|
||C||H||N||S||% leaves||% bark||% moisture content||HHV (kJ/kg)|
|Mixing three species||HHV = 717.79 + 380.85·C + 7.61% leaves−14.00% bark||97.352||278.77||195.74|
|HHV = 863.61 + 383.93·C − 15.257% bark||96.976||297.90||218.00|
|HHV =769.58 + 371.80·C + 8.02% leaves||97.147||289.35||201.07|
|HHV = 928.66 + 374.30·C||96.733||309.62||223.64|
||HHV = 7624.11 + 237.76·C + 228.75·N – 83.37%w||97.131||166.21||125.38|
|HHV = −1024.04 + 412.99·C + 272.54·N||96.743||177.08||133.28|
|HHV = −1106.84 + 422.77·C||94.630||227.38||173.87|
||HHV = 4323.54 + 301.46·C – 72.74% bark||88.127||331.30||250.40|
|HHV = 1254.39 + 365.79·C||83.17||394.46||303.26|
|HHV = −790.92 + 421.37·C||94.858||231.17||118.44|
The calorific value was measured by means of a LECO AC-500 isoperibolic calorimeter. Before analyzing the samples, the calorimeter was calibrated by the combustion of a reactive standard of a known calorific power (benzoic acid, 1 g). Subsequently, each sample was prepared with a mass between 0.1 and 1 g. This was introduced into a combustion vessel where a fuse wire caused ignition. Note that 10 ml of distilled water was added. Then, combustion vessel was sealed, and oxygen with a pressure of 3000 kPa was introduced inside the calorimeter . This container was placed in a bucket of water which was surrounded by an insulating layer to maintain a constant temperature. During analysis, the water temperature was measured by an electronic thermometer with an accuracy of 1/10,000 degree. In order to control energy exchange, the temperatures of the cuvette and the insulating layer were continuously monitored. With this, the device applies a correction to the result.The water temperature was monitored by a microprocessor in every 6 seconds. The difference between the water temperature before ignition and after ignition was processed through the calorimeter software, for obtaining the calorific value.
The weight percentage of carbon, hydrogen, and nitrogen was measured by means of a LECO TruSpec CHN analyzer. According to EN 14918, samples between 0.1 and 1 g were weighed. Then, they were wrapped in titanium sheets that is completely inorganic. These were placed in a feeding carrousel. The analysis cycle consists of three phases: purging, combustion, and analysis. In the purge phase, the sample is casted into the load compartment, which is sealed and atmospheric gases are removed. In the second phase, the sample is casted into a compartment at 950°C, and oxygen is injected for rapid and complete combustion. The gases pass through an afterburner at 850°C to oxidize and remove particles. The combustion gases are collected in the ballast (gas collection vessel). During the analysis phase, combustion gases are homogenized in the ballast. Subsequently, an aliquot of 3 cm3 is purged with helium through infrared detectors of CO2 and H2O. Another gas loop aliquot is transported through hot copper to remove O2 and transform the NO
In order to calculate the cost of analysis all inputs were counted. Market prices of nine enterprises were checked. Consumables, labor maintenance, and amortization were considered. Cost of technical labor was estimated in 20 €/h. Residual value of the device was considered to be 10% of investment. Time of analysis was measured.
2.4. Percentage of bark
The percentage of bark in the branches was calculated after their classification according to six diameter classes: 0–1, 1–2, 2–3, 3–4, 5–6, and >6 cm:. The diameter influences the bark ratio . For each diameter, class five samples were taken, so 30 samples in total were analyzed in each species. In each branch, diameter with bark and diameter without bark were measured using of a digital caliper with 0.01 mm accuracy, as shown in Figure 3. From these diameters, the percentage of bark was calculated by Eq. (1) , where
2.5. Determination of wood density
Wood density is expressed as the mass of dry wood per unit volume. To calculate the density, first the samples were immersed in a beaker with water, calibrated to 250 mL. The volume was measured by the difference between the water level before and after immersing. Then the dry weight of the samples was determined; for this purpose, samples were placed in a drying oven at a constant temperature of 105 ± 2°C for 24 h. Mean and standard deviation for the densities are obtained by Eq. (2) :
2.6. Prediction model of heat value
To obtain predictive models by regression up to three variables (C, H, and N) have been used. For evaluating the models, the coefficient of determination (
For all equations, 30 data were used to develop the models, and another 15 independent data were used for validation. The statistical program used was Statgraphics Centurion XV©14, even for the calculation of the significance of the variables of the mathematical prediction models by the beta coefficients and Student’s
3. Results and discussion
Several tests were initially applied to check the normal distribution of the data, such as Shapiro-Wilk test [27, 28], Anderson-Darling test, and the Lilliefors test [29, 30]. In Table 3, statistical description of the studied variables is shown for each species. It is observed that the coefficient of skewness and kurtosis are between −2 and +2. This indicates that the observed samples are normally distributed.
The calorific values are clearly influenced by the moisture content and the leaf content in the sample. The obtained calorific value for the three species was between 12 and 16 MJ/kg, which were relatively lower than values cited in the literature. This is because, in this work, we have studied samples of diverse diameters with different percentage of bark, without any drying process, and mixtures of wood and leaves in different proportions, so smaller high heat values (HHVs) were obtained. For example, González et al.  gave values of HHV for biomass from orange tree pruning as 16–18 MJ/kg. Yin  also analyzed mixtures of biomass and got HHV values 18 MJ/kg, but all these were measured on a dry condition basis.
To compare the calorific value of the three species studied, analysis of variance was performed. In Figure 4, LSD intervals are shown at 95% confidence level. It can be seen that the HHV of almond and olive trees were similar, but that of orange was significantly lower, which may be due to the characteristics of wood, leaves, and bark with moisture.
Table 4 shows the correlation between analyzed variables. Significant negative influence of the percentage of moisture in the calorific value (HHV) is observed, with −0.97 Pearson coefficient. On the other hand, it is noted that a higher percentage of C increases the HHV (+0.97). It is observed that the percentage of H is associated with the highest moisture content (+0.63) and obviously decreases the percentage of carbon in the sample and the calorific value (−0.54). It is observed that the percentages of sulfur, leaves, and bark do not have a clear influence on the calorific value of almonds but have a negative influence on the olive and orange trees.
Calorific models from elemental analyses are proposed in Table 5. All models show a high
Vargas-Moreno et al.  conducted a review of models to predict the calorific value for different biomass materials based on elemental analyses. Most of these models reviewed gave determination coefficients between 0.8 and 0.99. For this reason, the correlation coefficients obtained in this work, higher than 0.8, were considered very acceptable (Table 4). However, on the other hand, Acda  proposed models with higher
Velázquez-Martí et al.  already applied this method to obtain prediction models to predict a high heat value on lignocellulosic waste materials from urban tree pruning. In these works, the determination coefficients obtained are similar to those of the present work, between 0.7 and 0.9. This means that the presence of moisture causes decrease in accuracy of the calorific value calculated by prediction models.
It is observed that a higher carbon content (C) provides bigger calorific values. This coincides with many other studies [2, 34]. It is also shown in Table 4 with a Pearson correlation coefficient of 0.97 between HHV and C. Because of this, all proposed equations for predicting the gross calorific value (Table 5), whether univariate or multivariate, present the variable C.
The sulfur (S) is not high in all species studied (Table 3), so it would not be a problem in biomass combustion boilers caused by this element. As it can be seen, in olive and citrus wood the S and N contents are lower than the limits established by the standard EN 14691-part 4 , which fix the conditions for chips used as biofuels in %N < 1% and %S < 0.1%. Almond has not got values excessive high. This allows preparing mixtures with materials with low N and S contents.
This discrepancy may be due to the extractive substances of plant biomass (sugars, tannins, sterols, fatty acids, resin acids, oligomeric terpenics, hydrocarbons, etc.) that influence the gross calorific value. Their contents depend on the species, the part of the plant, season, and the growing stage, among other factors . Another explanation could be the botanical family of the studied species .
Once calculated using the prediction models, we proceeded to the validation of these data by applying them to new samples, whose observations were compared with the predicted values.
3.1. Comparison cost between direct and indirect measurement of HHV
Measurements with the isoperibolic LECO AC-500 calorimeter showed that the time per analyzed sample was 15 min, which includes sample preparation. As consumables, 6 L of pure oxygen measured at standard conditions (25°C and 1 atm of pressure) and a fuse for ignition of the sample are used. Table 6 shows the cost of reagents. Device depreciation and calibration cost are respectively shown in Tables 7 and 8. The total costs for different number of analyzed samples are shown in Table 9.
|Consumable||Price per consumable pack||Consumption per sample||Analysis cost (€/analysis)|
|Standard for calibration ref.
502-208-T Ac. Benzoic
|160.0||€/100 tablets||5 tablets/calibration||1.60|
|Spike for combustion ref.
502-815 Mineral oil
|Fuse for ignition ref. 502-462||38.2||€/375 fuse wire||1||Ud||0.10|
|Rent gas tank O2||94.4||€/year||0.019|
|Residual value 10%||5000||€|
|Time (min/analysis)||Total time (min)||Cost (€)|
|Number of analysis||5||20||100||33.33|
|Standard for calibration ref. 502-208-T Ac. Benzoic||5||8.00|
|Fuse wire, ref. 502–462||5||0.51|
|Rent gas tanks + maintenance||1.294|
|Sample number||Number of calibration||Technical labor time (h)||Technical labor cost (€)||Gasses (€)||Wire fuse (€)||Cost/sample (€/sample)|
The time spent for elemental analysis using a TruSpec CHN analyzer is 5 min. The costs of the reagents, calibration, and amortization are respectively shown in Tables 10–12.
|Consumable||Price per consumable pack||Consumption per sample||Analysis cost (€/analysis)|
|Standard for calibration EDTA ref. 502-092||44.8||€/50 g||0.2 × 8||g/calibration||0.18|
|Tin Foil for solid samples
(Large Tin Foil ref. 502-397-400)
|65.5||€/400 tin foil||1||Unit/sample||0.16|
|Rent tank O2||94.4||€/year||0.019|
|Rent tank N2||47.2||€/year||0.009|
|Rent tank He||47.2||€/year||0.009|
|Rent tank compressed air||47.2||€/year||0.009|
|Gas O2||450.0||€/tank||1||Tank/200 analysis||2.25|
|Gas N2||342.9||€/tank||1||Tank/200 analysis||1.71|
|Gas He||342.9||€/tank||1||Tank/200 analysis||1.71|
|Time (min/analysis)||Total time (min)||Cost|
|Time of calibration calculation||Blanc test||15||5||75|
|Standard for calibration||EDTA ref. 502-092||8||1.43||€|
|Tin foil for solid samples||Large Tin Foil ref. 502-397-400||8||1.31||€|
|Rent gas tanks + maintenance||8||2.30||€|
|Residual value 10%||12,500||€|
In Table 13, the total costs of the elemental analysis with the analyzer LECO CHNS TruSpec are shown. When Tables 9 and 13 are compared, it can be seen that the cost for determining the calorific value indirectly from elemental analysis is 23% cheaper than the direct measurement with AC500 LECO isoperibolic calorimeter for 25 samples. Moreover, time of determination is lower. The possibility to calculate the calorific value of a substance from its elemental composition reduces investment to a single computer, instead of two. This is very important in laboratories with limited resources of small and medium enterprises.
|No. analysis||Calibration||Calibration cost||Analysis time (h)||Technician labor||Gases||Tin foil||Cost/sample (€/sample)|
The advantage of using indirect methods for determining the heat value based on regression models from the analysis of C, H, and N elements has been proven this paper. Time to determine the high heat value indirectly is 40% lower than the time taken using the calorimeter directly. The cost of indirect method is 23% cheaper. Along with the cost savings, reduced analysis time is associated with a lower environmental impact linked to the reagents used.
It is proved that a higher carbon content (C) provides bigger calorific values. However, hydrogen (H) has a negative influence, with a negative Pearson’s coefficient. This is due to the fact that hydrogen is associated with the water content. As it is known, the moisture content decreases the high heat value of the biomass, therefore, hydrogen presents inverse proportionality with the heat obtained from combustion. The colinearity between the moisture content and the hydrogen ratio justifies that both were rarely considered in the same model to predict the heat value.
In this paper, models for determining the calorific value in samples of
According to the thermochemical characterization of plum wood, the residual biomass from pruning can be used as chips for bioenergy. These species did not have significant differences in C and H composition, between 30 and 40% C, between 5 and 7% H. However, small differences exist with respect to N and S. Olive and citrus wood have S and N contents lower than the limits established by the standard EN 14691-part 4 , which fix the conditions for chips used as biofuels in %N < 1% and %S < 0.1%. Almond has not got values excessive high but still it exceeds. This leads to the preparation of mixtures with materials with low N and S contents to decrease their content.
Friedl A, Padouvas E, Rotter H, Varmuza K. Prediction of heating values of biomass fuel from elemental composition. Anal Chim Acta. 2005; 544:191–198. DOI: 10.1016/j.aca.2005.01.0
Sheng C, Azevedo JLT. Estimating the higher heating value of biomass fuels from basic analysis data. Biomass Bioenergy. 2005; 28(5):499–507. DOI: 10.1016/j.biombiec.2004.11.008
Callejón-Ferre AJ, Velázquez-Martí B, López-Martinez JA, Manzano-Agugliaro F. Greenhouse crop residues: energy potential and models for prediction of their higher heating value. Renew Sust Energ Rev. 2011; 15(2):948–955. DOI: 10.1016/j.rser.2010.11.01
Vargas-Moreno JM, Callejón-Ferre AJ, Pérez-Alonso J, Velázquez-Martí B. A review of the mathematical models for predicting the heating value of biomass materials. Renew Sust Energ Rev. 2012; 16:3065–3083. DOI: 10.1016/j.rser.2012.02.054
Velázquez-Martí B, Sajdak M, López-Cortés I, Callejón-Ferre AJ. Wood characterization for energy application proceeding from pruning Morus albaL., Platanus hispanicaMünchh. and Sophora japonicaL. in urban áreas. Renew Energy. 2014; 62:478–483. DOI: 10.1016/j.renene.2013.08.010
Pérez-Arévalo JJ, Callejón-Ferre AJ, Velázquez-Martí B, Suárez-Medina MD. Prediction models based on higher heating value from the elemental analysis of neem, mango, avocado, banana, and carob trees in Guayas (Ecuador). J Renew Sustain Ener. 2015; 7:053122.
Saidur R, Abdelaziz EA, Demirbaş A, Hossain MS, Mekhilef S. A review on biomass as a fuel for boilers. Renew Sust Energ Rev. 2011; 15(5):2262–2289. DOI: 10.1016/j.rser.2011.02.015
Khan AA, Jonga WD, Jansens PJ, Spliethoff H. Biomass combustion in fluidized bed boilers: potential problems and remedies. Fuel Process Technol. 2009; 90:21–50. DOI: 10.1016/j.fuproc.2008.07.012
Telmo C, Lousada J, Moreira N. Proximate analysis, backwards stepwise regression between gross calorific value, ultimate and chemical analysis of wood. Bioresource Technol. 2010; 101(11):3808–3815. DOI: 10.1016/j.biortech.2010.01.021
Demirbas A, Demirbas AH. Estimating the calorific values of lignocellulosic fuels. Energ Explor Exploit. 2004; 22:135–143. DOI: 10.1260/0144598041475198.
Demirbas A. Relationship between heating value and lignin, moisture, ash and extractive contents of biomass fuels. Energy Explor Exploit. 2002; 20:105–111. DOI: 10.1260/014459802760170420
Callejón-Ferre AJ, Carreño-Sánchez J, Suárez-Medina FJ, Pérez-Alonso J, Velázquez-Martí B . Prediction models for higher heating value based on the structural analysis of the biomass of plant remains from the greenhouses of Almería. Fuel. 2014; 116:377–387. DOI: 10.1016/j.fuel.2013.08.023
UNE-EN 164001:2005 EX. Solids Biofuels. Method for the detemination of calorific value. AENOR: Madrid. Spain; 2005.
Francis HE, Lloyd WG. Predicting heating value from elemental composition. J Coal Qual. 1983; 2(2):21–25.
Ebeling JM, Jenkins BM. Physical and chemical properties of biomass. Trans ASAE. 1998; 28:898–902.
Kathiravale S, Yunus MNM, Sopian K, Samsuddin AH, Rahman RA. Modeling the heating value of municipal solid waste. Fuel. 2003; 82(9):1119–1125. DOI: 10.1016/S0016-2361(03)00009-7
Velázquez-Martí B, Fernández-González E, López-Cortes I, Salazar-Hernández DM. Quantification of the residual biomass obtained from pruning of vineyards in Mediterranean area. Biomass Bioenergy. 2011; 35(3):3453–3464. DOI: 10.1016/j.biombioe.2011.04.009
Velázquez-Martí B, Fernández-González E, López-Cortes I, Salazar-Hernández DM. Quantification of the residual biomass obtained from pruning of trees in Mediterranean olive groves. Biomass Bioenergy. 2011; 35(2):3208–3217. DOI: 10.1016/j.biombioe.2011.04.042
INE. Instituto Nacional de Estadística. Censo Agrario 2009 [Internet]. 2009 [Updated: 2014]. Available from: http://www.ine.es/censoagrario/censoag.htm
Avanza MM, Bramardi SJ, Mazza SM. Statistical models to describe the fruit growth pattern in sweet orange “Valencia late”. Spanish J Agric Res. 2008; 6(4):577–585.
Velázquez-Martí B, Fernández-González E, López-Cortes I, Salazar-Hernández DM. Quantification of the residual biomass obtained from pruning of trees in Mediterranean almond groves. Renew Energy. 2011; 36:621–626. DOI: 10.1016/j.renene.2010.08.008
López-Cortés I, Salazar-García DC, Salazar DM, Velázquez-Martí B. Fatty Acid, Sterol and Polyphenol content in most Monovarietal Oils of the Spanish East Areas. In: De Leonardis A. (eds.). Virgin olive oil. Production, composition, uses and benefit for the man. New York: Nova Publishers; 2014.
Barranco D, Rallo L. Olive cultivars in Spain. Horttechnology. 2000; 10(1):107–110.
UNE-EN 14961-1:2011. Biocombustibles sólidos. Especificaciones y clases de combustibles. Parte 1: Requisitos generales. AENOR: Madrid, Spain; 2011.
UNE-EN 14780:2012. Biocombustibles sólidos. Preparación de muestras. AENOR: Madrid, Spain; 2012.
Husch B, Beers TW, Kershaw JA Jr, editors. Forest Mensuration. John Wiley & Sons. Inc.: Hoboken, NJ; 2003.
Shapiro SS, Francia RS. An approximate analysis of variance test for normality. J Amer Statist Assoc. 1972; 67: 215–216.
Shapiro SS, Wilk MB. An analysis of variance test for normality (complete samples). Biometrika. 1965;52:591–611.
Sheskin DJ. Handbook of Parametric and Nonparametric Statistical Procedures: Fifth Edition. London: Chapman & Hall/CRC Taylor & Francis Group 2011.
Razali NM, Wah YB. Power Comparisons of Shapiro-Wilk, Kolmogorov-Smirnov, Lilliefors and Anderson-Darling Tests. J Stat Model Anal. 2011;2(1):21–33.
González Z, Rosal A, Requejo A, Rodríguez A. Production of pulp and energy using orange tree prunings. Bioresource Technol. 2011; 102(19):9330–9334. DOI: 10.1016/j.biortech.2011.07.088
Yin CY. Prediction of higher heating values of biomass from proximate and ultimate analyses. Fuel. 2011; 90(3):1128–1132. DOI: 10.1016/j.fuel.2010.11.031
Acda MN. Physico-chemical properties of wood pellets from coppice of short rotation tropical hardwoods. Fuel. 2015; 160:531–533. doi:10.1016/j.fuel.2015.08.018
Erol M, Haykiri-Acma H, Kücükbayrak S. Calorific value estimation of biomass from their proximate analysis data. Renew Energy. 2010; 35:170–173. DOI:10.1016/j.renene.2009.05.008
EN 14691-part 4. Wood chips—product standard for non-industrial use. CEN-UNE, Madrid, Spain, 2011.
Nhuchhen DR. Prediction of carbon, hydrogen, and oxygen compositions of raw and torrefied biomass using proximate analysis. Fuel. 2016; 180:348–356. doi: 10.1016/j.fuel.2016.04.058
Veiga JPS, Valle TL, Feltran JC, Bizzo WA. Characterization and productivity of cassava waste and its use as an energy source. Renew Energy. 2016; 93:691–699. doi:10.1016/j.renene.2016.02.078