Canadian Journal of Forest Research

Minimum age for clear -cutting native species with energetic potential in the Brazilian semi-arid

Journal: Canadian Journal of Forest Research

Manuscript ID cjfr-2016-0392.R3

Manuscript Type: Note

Date Submitted by the Author: 02-Dec-2016

Complete List of Authors: Santana, Otacilio; Universidade Federal de Pernambuco, Biophysics and Radiobiology

Keyword: growth rings,Draft firewood, Semi-Arid, rotation age, bioenergy

https://mc06.manuscriptcentral.com/cjfr-pubs Page 1 of 24 Canadian Journal of Forest Research

1 Minimum age for clear-cutting native species with energetic potential in the Brazilian

2 semi-arid region

3

4 Otacilio Antunes Santana 1

5 Av. Prof. Moraes Rego, n° 1235, Departamento de Biofísica e Radiobiologia (DBR),

6 Mestrado Profissional em Rede Nacional para o Ensino das Ciências Ambientais

7 (PROFCIAMB), Centro de Biociências (CB), Universidade Federal de Pernambuco (UFPE),

8 Cidade Universitaria, Recife, Pernambuco, Brazil, 50670901

9 [email protected]

10 [email protected]

11 +55 81 997371266

12 1Author for correspondence Draft

13

14

15

16

17

18

19

20

21

22

23

24

25

1 https://mc06.manuscriptcentral.com/cjfr-pubs Canadian Journal of Forest Research Page 2 of 24

26 Abstract

27 What is the minimum age for clearcutting native species with energetic potential in the

28 Brazilian semiarid region? This was the central question of this study, which aimed at

29 estimating the minimum cutting cycle for native woody species from an area covered by a

30 sustainable forest management plan. Individuals (n = 240) of twelve species native to the

31 semiarid region were measured and wood discs were taken from them to collect data relative

32 to: i) volume, by an industrial 3D scanner; ii) density, using a highfrequency densitometer;

33 iii) age, by a growth ring measuring device; and iv) the lower heating value, by an oxygen

34 bomb calorimeter and an elemental analyzer. Data were fitted to mathematical models by

35 regression analysis to determine the relationship and significance between the variables

36 analyzed. The estimated minimum age for the harvesting of native woody species was 47

37 years, determined by the age at whichDraft there was stabilization of volume growth curves over

38 time, and an increase in density and heating power (wood retraction).

39

40 Keywords: growth rings; firewood; rotation age; bioenergy.

41

42

43

44

45

46

47

48

49

50

2 https://mc06.manuscriptcentral.com/cjfr-pubs Page 3 of 24 Canadian Journal of Forest Research

51 1. Introduction

52 The mean annual increment (MAI) and the current annual increment (CAI) are parameters

53 that allow estimation of the optimal age of timber harvest under sustainable forest

54 management, but their calculation is laborious (costly, expensive and requires a large data set)

55 and time consuming in woody species of the Brazilian semiarid region. Although sustainable

56 forest management plans are in place for this region, such difficulty results in ages for clear

57 cut cycles without a scientific basis to justify these possible cycles. Unlike rapidly growing

58 species of Eucalyptus and Pinus that invest their biomass mainly in shoots, native woody

59 species of the Brazilian semiarid region have slow initial growth, investing most of their

60 biomass initially in roots (high root to shoot ratio); this as a support structure, to reduce water

61 loss and against herbivorous attack (from which organs find refuge in the soil) (Lima and

62 Rodal 2010, Costa et al. 2014, SantanaDraft and Encinas 2016). In addition, the individual

63 architecture of woody in the Brazilian semiarid region is based on broad ramifications

64 of the trunk and branches, which makes it difficult to estimate the independent variables,

65 namely height and diameter, used in volumetric calculations (Barros et al. 2010).

66 In addition to shoot volume, other variables are essential for estimating the energy

67 potential of wood. The bulk density, chemical composition (e.g. cellulose and lignin),

68 percentage of vessels and parenchyma, moisture content and size of fibers are some

69 parameters that influence the heating value of wood. However, these parameters are

70 correlated and represented by the density of the wood, which has a significant direct and

71 proportional relationship with the heating value of wood (Gebreegziabher et al. 2013). The

72 age of the individual is another important factor, because, with advancing age and the

73 stabilization of growth in height, woody individuals increase the carbon accumulation on their

74 trunks and branches, thus increasing the value of the mass/volume ratio of the timber. Adult

75 wood is characterized by high density, long tracheids, thick cell walls, high percentage of

3 https://mc06.manuscriptcentral.com/cjfr-pubs Canadian Journal of Forest Research Page 4 of 24

76 latewood, low percentage of spiral grain, low percentage of nodes, high percentage of

77 cellulose, low percentage of compression wood, high transverse shrinkage, smaller microfibril

78 angle and greater mechanical strength (Rowell 2012).

79 The optimal age for harvesting a planted forest stand for wood production is when the

80 current annual increment curve intersects the mean annual increment curve. At that point, the

81 rate of growth in wood volume begins to reduce over time, and the economic profitability in

82 relation to forest investment also begins to fall (Prodan, 1968). In semiarid regions, in

83 addition to a reduced rate of wood volume growth as trees age, wood density often increases

84 over time as an ecophysiological strategy of native tree species to cope with moisture

85 limitations (Liu et al. 2013, Klein et al. 2014, Santana and Encinas 2016). This wood

86 retraction is evident in the aerial part of the plant and can enhance the root:shoot ratio even if

87 the area index continues to developDraft over time (Santana and Encinas 2016). Before of the

88 wood retraction occurs, a reduction of the rate of growth in wood volume and an increasing of

89 the rate of growth in wood density are observed (Paula 1993, Santana and Encinas 2016).

90 Thus, in semiarid regions, the time at which wood volume growth levels off and when the

91 density and LHV curves change their inclination angle in relation to age could represent a

92 parameter defines the best time for harvesting wood destined for bioenergy uses (Liu et al.

93 2013, Klein et al. 2014). There is not an analytical or exact method for calculating this point

94 of levelling off or for detecting a critical change in the inclination angle of wood density and

95 LHV curves, but these trends are supported by observations in the literature (Althoff et al.

96 2016, Santana and Encinas 2016).

97 In 2013, 9.1% of the Brazilian energy matrix was based on wood used as firewood and

98 charcoal, which has increased the demand for this resource by approximately 5% per year

99 (Brasil 2013). In the semiarid region of the state of Pernambuco, the extracted wood is used

100 as an energy source of the Polo Gesseiro do Araripe (Araripina Microregion), which has 39

4 https://mc06.manuscriptcentral.com/cjfr-pubs Page 5 of 24 Canadian Journal of Forest Research

101 gypsum mines and 869 industries to cement manufacture, which fill 73% of their energy

102 requirements with wood, but the small manufacturers use 100% wood to meet their energy

103 needs. The production of a 10,000 kg of plaster requires 0.5 st native wood, which has an

104 annual demand of 5.5·10 4 million kg/year (Silva 2008/2009). The sustainable production of

105 wood fuel would require a total area of 155,000 hectares specifically for the plaster industry,

106 with a cutting cycle ranging from 10 to 15 years. But the Northeastern Plants Association

107 (APNE 2014) has the registration of 94 sustainable forest management plans in the state of

108 Pernambuco, covering an area of 39,748 hectares, insufficient to meet the energy demands of

109 the Polo, which grows 10% (Sindugesso 2011) to 25% per year, in addition to a growing

110 ceramics industry (CPRH 2014).

111 The goal of this study was to estimate the minimum age for a clear cutting cycle of

112 native woody species from an area governedDraft by a sustainable forest management plan in the

113 semiarid region, according to requirements for implementation and maintenance of Annual

114 Production Unit, described by the Normative Instruction of the Brazil Ministry of the

115 Environment (number 1 of 06/25/2009 – MMA 2009).

116

117 2. Materials and Methods

118 The wood analyzed came from a forest management farm (Figure 1A) in the state of

119 Pernambuco (8°35’S and 37°59’W), having a semiarid climate with annual rainfall below

120 1,000 mm (550 mm in 2015, APAC, 2016) and mean annual temperature greater than 27°C,

121 or a BSh climate according to the Köppen classification (Peel et al. 2007). This farm is

122 located on the ‘Sertaneja Meridional e Raso da Catarina’ geomorphological depression. The

123 soils have been classified into four classes: Haplic Cambisol, Salic Gleysol, Yellow Oxisol,

124 Haplic Planosol, and the vegetation has been classified into Closed Arboreal, Open Arboreal

125 and Shrub structural classes (Aguiar et al. 2013).

5 https://mc06.manuscriptcentral.com/cjfr-pubs Canadian Journal of Forest Research Page 6 of 24

126 The vegetation analyzed came from the clearcutting of experimental plots in a

127 sustainably managed forest of the SemiArid Forest Management Network (MMA 2009). The

128 twelve predominant tree species have a cumulative importance value index of 96%, measured

129 in 20 plots of 50 x 100 m (10 ha), with a 1,271 stems/ha density (diameter at breast height ≥ 1

130 cm) (Aguiar et al. 2013), evaluated according to the Permanent Plot Measurements Protocol

131 of the SemiArid Forest Management Network, and Rodal et al. (1992).

132 The species and number of evaluated individuals of each species were: Acacia

133 kallunkiae J.W.Grimes & Barneby, (n = 13); Acacia piauhiensis Benth., Fabaceae

134 (n = 15); Anadenanthera colubrina (Vell.) Brenan. Var. cebil (Griseb.) Reis, Fabaceae (n =

135 27); Aspidosperma pyrifolium Mart., Apocynaceae (n = 35); Poincianella pyramidalis (Tul.)

136 L. P. Queiroz, Fabaceae (n = 12); Croton blanchetianus Baill., Euphorbiaceae (n = 9);

137 Erythrina velutina Willd., PapilionoideaeDraft (n = 17); Jatropha elliptica (Pohl.) Mull.Arg.,

138 Euphorbiaceae (n = 33); Maytenus rigida (Mart.) Benth., Celastraceae (n = 16), Mimosa

139 caesalpiniifolia Benth, Fabaceae (n = 8); Myracrodruon urundeuva All., Anacardiaceae (n =

140 29); and dubium (Spreng.) Taub., Fabaceae (n = 26).

141 Harvested trees were analyzed for volumetry (Figure 1B), densitometry (Figure 1C),

142 dating (Figure 1C) and calorimetry (Figure 1D). The volume of the wood (m 3) was measured

143 by an industrial 3D scanner (3D scanner SKDKFX four Lens, Foshan Shangke

144 Machinery Co., Ltd., Guangdong Province, China), in which all parts (trunk + branches) were

145 reduced to 1 m length (Figure 1B) in a mobile studio in locus. The density of the wood and

146 the age (growth rings) were measured from sectioned discs, 2 cm thick, by means of high

147 frequency densitometry (Schinker et al. 2003), in LignoStation TM using LignoScop (Rinntech,

148 Heidelberg, Germany) (Figure 1C) to scan surfaces using cameras attached to the microscope,

149 and LignoScan to scan the wood surface using highfrequency (1600 dpi) with 0.001 mm

150 precision (Shchupakivskyy et al. 2014). Standard procedures (Fritts 1976, Lisi et al. 2008)

6 https://mc06.manuscriptcentral.com/cjfr-pubs Page 7 of 24 Canadian Journal of Forest Research

151 and COFECHA software (The University of Tennessee, Knoxville, GrissinoMayer 2001)

152 were used for crossdating and standardization. The irregularity of interannual rain, or of

153 another climate event, could be revised by the crossdating analysis. This analysis could

154 identify growth zone and false growth rings in an interannual period (Pagotto et al. 2015).

155 The equipment performed 1,200 factorial (1,200!) crossdating analyzes (240 studied

156 individuals with five replicates each).

157 The Higher Heating Value (HHV) of the wood was calculated by collecting data

158 measured with a calorimeter (C6000 global standards 2/10 IKA ®, Staufen, Germany)

159 (Günther et al. 2012). Samples were ground to a maximum size of 25 µm. Lower Heating

160 Value (LHV) was determined by subtracting the heat of vaporization of the water vapor from

161 the higher heating value. The water vapor was determined by variance of the hydrogen

162 content of the sample measured throughDraft the ElementAnalyzer (Vario EL, Elementar

163 Analysesystem GmbH, Langenselbold, Germany) (Figure 1D). These procedures followed the

164 standards DIN EN ISO 17162009, DIN 519001 2000, DIN 519003 2005 (Günther et al.

165 2012) and ABNT/NBR 8633/84.

166 The relationships between variables (Y = volume and X = density; Y= volume and X

167 = age; Y = lower heating value and X = density; Y = density and X = age; Y = lower heating

168 value and X = volume; and Y = lower heating value and X = age) were performed by fitting

169 the data to a wide array of mathematical growth models (Table 1) using regression analysis to

170 calculate a coefficient of determination (R2), rootmeansquare error (RMSE), significance

171 level (p) and the fit curve; and the selection of the best fit model (based on maximizing R2,

172 minimizing RMSE, and minimizing p) for each case (Zar 1999). The regression analysis was

173 preceded by the D’Agostino Normality test (D'Agostino et al. 1990) for each variable to

174 validate statistical premises. The models were chosen as indicated by Kleinbaum et al. (2013).

175 The possible multicollinearity among variables was calculated by the FarrarGlauber test

7 https://mc06.manuscriptcentral.com/cjfr-pubs Canadian Journal of Forest Research Page 8 of 24

176 (Farrar & Glauber 1976). Adjustments, graphics, statistical deviation and coefficient of

177 variation were calculated using Statistica 12 (Statsoft, Dell, Tulsa, USA).

178

179 3. Results

180 The mean age of the study population was 46 years (CV% = 15.2%; Figure 2A), ranging from

181 98 years ( Peltophorum dubium ) to 5 years ( Croton blanchetianus ). The mean volume was

182 0.12 m 3 (CV% = 11.9%; Figure 2B) per individual tree, with individuals ranging in size from

183 a minimum volume of 0.003 m 3 (Croton blanchetianus ) to a maximum volume of 0.40 m 3

184 (Acacia kallunkiae ). The mean wood density was 0.8 g · cm 3 (CV% = 7.1%; Figure 2C),

185 varying between 0.4 g · cm 3 (Jatropha elliptica ) and 1.05 g · cm 3 (Acacia kallunkiae ). The

186 mean lower heating value was 4,863 kcal · kg 1 (CV% = 6.2%; Figure 2D), ranging from

187 3,783 kcal · kg 1 (Acacia piauhiensis )Draft to 5,697 kcal · kg 1 (Erythrina velutina ).

188 Data reflecting the relationship between volume and density, and between volume and

189 age, were well fit to sigmoidal models (5 parameters) significantly (R2 > 0.86; RMSE < 0.05;

190 p < 0.001; Figure 3A and 3B, Table 2). In both relationships, volume began to rise sharply

191 (curve growth > 45° slope relative to the dependent variable) to the point where it became

192 constant (curve growth < 5° slope in relation to the dependent variable) and then slowed

193 relative to growth of the independent variable (density or age). The relationship between the

194 lower heating value (LHV) and density was the most significant (R2 > 0.97; RMSE < 0.03; p

195 < 0.001), directly and significantly increasing with increasing density (Figure 3C, Table 2).

196 Multicollinearity was not found in the models tested (p > 0.800). Considering the relationship

197 between density and LHV with the age of individuals, data fitted to the exponential growth

198 model (double and with five parameters), in which, from a given year, 47 years for density

199 and 49.5 for LHV, there is a levelling off of the curve (Figure 3D and 3F; Table 2), with a

200 more marked increase in the values of these variables (curve growth> 45 °slope relative to the

8 https://mc06.manuscriptcentral.com/cjfr-pubs Page 9 of 24 Canadian Journal of Forest Research

201 dependent variable) than previously (growth curve ≈ 30 ° slope relative to the dependent

202 variable). As for the relationship between LHV and volume, data fitted to the Chapman model

203 (four parameters) (Figure 3E; Table 2), indicating that in the higher volume growth phase, the

204 increase in LHV value is less pronounced (curve growth < 15° slope in relation to the

205 dependent variable) than when the volume growth begins to stabilize (growth curve > 45°

206 slope in relation to the dependent variable).

207 The minimum age for clear cutting cycle in an area of Sustainable Forest Management

208 Plan in the semiarid region, with the destination of wood with energy purposes, suggested by

209 the variables analyzed, was from 46.8 (47) years (Figure 3B). This is the average tree age at

210 which there is a stabilization of the volume growth (curve growth < 5° slope in relation to the

211 dependent variable), and pronounced growth of the LHV (from 49.5 years; Figure 3F) and

212 density (from 47 years; Figure 3D). Draft

213

214 4. Discussion

215 One limitation of this study lies in the impracticality of working with a larger sample size,

216 given the time required for cutting and loading for transportation (about 24 hours), making it

217 difficult to scan other woody individuals in the semiarid region. However, 240 individuals

218 were measured with five replicates of each measurement or analysis, improving the reliability

219 and accuracy of data collected.

220 Our findings corroborate what has already been reported for the area (Aguiar et al.

221 2013), which before the Sustainable Forest Management Plan, had reduced logging (> 20

222 individuals / ha / year). The values of volume, density, and LHV varied within the 95%

223 confidence intervals reported in the literature for the mean values of the species analyzed:

224 volume (from biomass) from 0.001 m 3 (Santana and Souto 2006) to 0.50 m 3 (Silva and

225 Sampaio 2008); density from 0.3 g · cm 3 (Lima and Rodal 2010) to 1.2 g · cm 3 (Paula 1993);

9 https://mc06.manuscriptcentral.com/cjfr-pubs Canadian Journal of Forest Research Page 10 of 24

226 and LHV from 3,800 kcal · kg 1 (Medeiros Neto et al. 2012) to 5,800 kcal · kg 1 (Quirino et al.

227 2004).

228 The rapid growth in the volume of tree species, over the years, and a subsequent

229 stabilization of this growth is commonly observed in the literature, not only for fast growing

230 plantation species (Encinas et al. 2011), but also for native species (Felker 1986, Lima 1986,

231 Vico et al. 2015, Santana and Encinas 2016). This stabilization is considered as the saturation

232 period of the plant individual, both by limited environmental resources, by competition with

233 others in their environment, but also by individual senescence marked by the genetics of the

234 species (Hunt 1982, Felker 1986, Soares et al. 2011). The onset of the saturation, senescence

235 or slow growth process, as shown in Figure 3B from 47 years, is marked in forestry science

236 by the beginning of the cutting cycle for vegetation intended for production of firewood and

237 charcoal (Rode et al. 2014, Althoff Draftet al. 2016). When the volume tends to stabilize, other

238 variables, such as density (Figure 3 D) and heating value (Figure 3F), which have a high

239 correlation (R2 > 0.97), continue to grow in value. While a saturation of those values was not

240 observed in this work, it is, however, emphasized in the literature (Hunt 1982, Felker 1986).

241 This reduction of the rate of growth of the wood volume and an increasing of the rate of

242 growth of the wood density and lower heating value (LHV) infer in the begin of wood

243 retraction. Thus, from the age of 47, the harvesting is strongly indicated.

244 In other dry forests on Earth, the cutting cycle of wood is ranges from 63 to 97 years,

245 much longer than the 47 years derived here. The causes for these longer rotation periods are:

246 i) in India, to avoid urban occupation (Agarwal et al., 2016); ii) in Australia and Africa, to

247 conserve the habitat for large mammals (Bhadouria et al. 2016); iii) in West Africa, to reduce

248 wildfire susceptibility and microclimate changes (Scheitera and Savadogo, 2016); iv) in

249 Cameroon and Panama, to preserve forest stands that are used in popular religious and sacred

250 rituals (Kemeuze et al. 2016, Seijo et al. 2016); v) in Costa Rica, India, Papua New Guinea,

10 https://mc06.manuscriptcentral.com/cjfr-pubs Page 11 of 24 Canadian Journal of Forest Research

251 and Southern Africa, to produce medicinal bioproducts (Schmiedel et al. 2016); vi) in

252 Morocco and Argentina, to perpetuate ethnobotanical and other cultural uses of natural

253 resources (Martínez 2015, Blanco and Carrière 2016); and vii) in China, to maintain the

254 groundwater flow near to surface by hydraulic lift from trees (Xiao and Huang 2016). In all

255 cases, like in this work, biodiversity preservation and firewood sustainability also were

256 objectives of forest management (Althoff et al. 2016, Santana 2016).

257

258 5. Conclusion

259 The estimate of the minimum age for cutting cycle of native woody species in the evaluated

260 area in Brazil’s semiarid region was 47 years, determined by the age at which there was

261 stabilization of volume growth over time, and an increase in density and heating power (wood

262 retraction). Draft

263

264 Acknowledgements

265 The author is grateful to ProReitoria de Pesquisa e PosGraduacao of the Universidade

266 Federal de Pernambuco (PROPESQ/UFPE) for the financial and logistical support, and to

267 Research Group ‘Educometria’ (UFPE/CNPq) by discussion and survey support. The author

268 is very grateful to the reviewers for their careful and meticulous reading of the paper.

269

270 References

271 Agarwal, S., Nagendra, H., and Ghate, R. 2016. The Influence of Forest Management

272 Regimes on Deforestation in a Central Indian Dry Deciduous Forest Landscape. Land. 5: 27

273 43. doi: 10.3390/land5030027

274 Aguiar, M. M. B., Santana, O. A., Inácio, E. dos S. B., Amorim, L. B. de, and Almeida

275 Cortez, J. S. 2013. Tree resilience after clearcutting in sustainable forest management of

11 https://mc06.manuscriptcentral.com/cjfr-pubs Canadian Journal of Forest Research Page 12 of 24

276 semiarids areas. In : 2nd United Nations Convention to Combat Desertification Scientific

277 Conference, Edited by Walter J. Ammann, UN, Bonn, p. 157158. Available from

278 [accessed 12 January 2016].

279 Althoff, T. D., Menezes, R. S. C., Carvalho, A. L., Pinto, A. S., Santiago, G. A. C. F., Ometto,

280 J. P. H. B., Randow, C., and Sampaio, E. V. D. S. B. 2016. Climate change impacts on the

281 sustainability of the firewood harvest and vegetation and soil carbon stocks in a tropical dry

282 forest in Santa Teresinha Municipality, Northeast Brazil. Forest Ecology and Management.

283 360 : 367375. doi: 10.1016/j.foreco.2015.10.001

284 APAC – Agência Pernambucana de Águas e Clima. 2016. Monitoramento Pluviometrico

285 [online]. Available from < http://www.apac.pe.gov.br/> [accessed 11 November 2016].

286 APNE – Associação de Plantas do Nordeste. 2014. Centro Nordestino de Informações sobre

287 Plantas. Planos de Manejo SustentáveisDraft da Caatinga [online]. Available from

288 [accessed 12 January 2014].

289 Associação Brasileira de Normas Técnicas – ABNT, 1984. NBR 8633: Carvão vegetal –

290 determinação do poder calorífico. ABNT, Rio de Janeiro. 12p.

291 Barros, B. C., Silva, J. A. A., Ferreira, R. L. C., and Rebouças, A. C. M. N. 2010. Volumetria

292 e sobrevivência de espécies nativas e exóticas no Pólo Gesseiro do AraripePE. Ciência

293 Florestal. 20 : 641647. doi: 10.5902/198050982422

294 Bhadouria, R., Singh, R., Srivastava, P., and Raghubanshi, A. S. 2016. Understanding the

295 ecology of treeseedling growth in dry tropical environment: a management perspective.

296 Energy, Ecology and Environment. 1: 296–309. doi: 10.1007/s4097401600383

297 Blanco, J., and Carrière, S. M. 2016. Sharing local ecological knowledge as a human

298 adaptation strategy to arid environments: Evidence from an ethnobotany survey in Morocco.

299 Journal of Arid Environments. 127 : 3043. doi: 10.1016/j.jaridenv.2015.10.021

12 https://mc06.manuscriptcentral.com/cjfr-pubs Page 13 of 24 Canadian Journal of Forest Research

300 Brasil. Empresa de Pesquisa Energética, 2013. Balanço Energético Nacional 2013 – Ano base

301 2012: Relatório Síntese. EPE, Rio de Janeiro. 55p.

302 CPRH – Agência Ambiental do Meio Ambiente de Pernambuco, 2014. Madeira Ilegal –

303 Buscas [online]. Available from

304 MADEIRA%20ILEGAL>. [accessed 01 June 2014].

305 Costa, T. L., Sampaio, E. V. S. B., Sales, M. F., Accioly, L. J. O., Althoff, T. D., Pareyn, F.

306 G. C., Albuquerque, E. R. G. M., and Menezes R. S. C. 2014. Root and shoot biomasses in

307 the tropical dry forest of semiarid Northeast Brazil. Plant and Soil. 378 : 113123. doi:

308 10.1007/s1110401320091

309 D'Agostino, R. B., Belanger, A., and D'agostino Jr., R. B. 1990. A suggestion for using

310 powerful and informative tests of normality. American Statistician. 44 : 316321. doi:

311 10.2307/2684359 Draft

312 Encinas, J. I., Santana, O. A., and Imaña, C. R. 2011. Volumetric and Economic optimal

313 rotations for firewood production of Eucalyptus urophylla in Ipamery, State of Goias.

314 Floresta. 41 : 905912. doi: 10.5380/rf.v41i4.25353

315 Farrar, D., and Glauber, R. 1976. Multicollinearity in Regression Analysis: the Problem

316 Revisited. Review of Economics and Statistics. 49 : 92107. doi: 10.2307/1937887

317 Felker, P. 1986. Establishment and productivity of tree plantings in semiarid regions. Elsevier,

318 Amsterdam. 444 p.

319 Fritts, H. C. 1976. Tree Rings and Climate. Academic Press, London. 566p.

320 GrissinoMayer, H. D. 2001. Evaluating crossdating accuracy: a manual and tutorial for the

321 computer program COFECHA. TreeRing Research. 57 : 205–221.

322 Gebreegziabher, T., Oyedun, A. O., and Hui, C. W. 2013. Optimum biomass drying for

323 combustion – A modeling approach. Energy. 53 : 6773. doi: 10.1016/j.energy.2013.03.004

13 https://mc06.manuscriptcentral.com/cjfr-pubs Canadian Journal of Forest Research Page 14 of 24

324 Günther, B., Gebauer, K., Barkowski, R., Rosenthal, M., and Bues C. T. 2012. Calorific value

325 of selected wood species and wood products. European Journal of Wood and Wood Products.

326 70 : 755–757. doi: 10.1007/s001070120613z

327 Hunt, R. 1982. Plant growth curves. The functional approach to plant growth analysis.

328 Edward Arnold Ltd., London, 248 p.

329 Kemeuze, V. A., Sonwa, D. J., Nkongmeneck, B. A., and Mapongmetsem, P. M. 2016.

330 Sacred groves and biodiversity conservation in semiarid area of Cameroon: Case study of

331 Diamare plain. In: Quels botanistes pour le 21e siècle? Métiers, enjeux et opportunités. Edited

332 by N. R. Rakotoarisoa, S. Blackmore and B. Riera. UNESCO, Paris. pp. 171183. Available

333 from [accessed 12 March 2016].

334 Klein, T., Hoch, G., Yakir, D., and Körner, C. 2014. Drought stress, growth and nonstructural

335 carbohydrate dynamics of pine treesDraft in a semiarid forest. Tree Physiology. 71 : 112. doi:

336 10.1093/treephys/tpu071

337 Kleinbaum, D., Kupper, L., Nizam, A., and Rosenberg, E. 2013. Applied regression analysis

338 and other multivariable methods. Cengage Learning, Michigan, 928p.

339 Lima, P. C. F. 1986. Tree productivity in the semiarid zone of Brazil. Forest ecology and

340 Management. 16 : 513. doi: 10.1016/03781127(86)900034

341 Lima, A. L. A., and Rodal, M. J. N. 2010. Phenology and wood density of plants growing in

342 the semiarid region of northeastern Brazil. Journal of Arid Environments. 74 : 13631373.

343 doi: 10.1016/j.jaridenv.2010.05.009

344 Lisi, C. S., Tomazello Filho, M., Botosso, P. C., Roig, F. A., Maria, V. B. R., FerreiraFedele,

345 L., and Voigt, A. R. A. 2008. Treering formation, radial increment periodicity, and

346 phenology of tree species from a seazonal semideciduous forest in southeast Brazil. IAWA

347 Journal. 29 : 189207. doi: 10.1163/2294193290000179

14 https://mc06.manuscriptcentral.com/cjfr-pubs Page 15 of 24 Canadian Journal of Forest Research

348 Liu, H., Park, W. A., Allen, C. D., Guo, D., Wu, X., Anenkhonov, O. A., Liang, E.,

349 Sandanov, D.V., Yin, Y., Qi, Z., and Badmaeva, N. K. 2013. Rapid warming accelerates tree

350 growth decline in semiarid forests of Inner Asia. Global Change Biology. 19 : 25002510.

351 doi: 10.1111/gcb.12217.

352 Martínez, G. J. 2015. Cultural patterns of firewood use as a tool for conservation: A study of

353 multiple perceptions in a semiarid region of Cordoba, Central Argentina. Journal of Arid

354 Environments. 121 : 8499. doi: 10.1016/j.jaridenv.2015.05.004

355 Medeiros Neto, P. N., Oliveira, E., Calegari, L., Almeida, A. M. C., Pimenta, A. S., and

356 Carneiro, A. C. O. 2012. Características FísioQuímicas e energéticas de duas espécies de

357 ocorrência no Semiárido Brasileiro. Ciência Florestal. 22 : 579588. doi:

358 10.5902/198050986624

359 Ministério do Meio Ambiente – MMA,Draft 2009. Dispõe sobre procedimentos técnicos para

360 elaboração, apresentação, execução e avaliação técnica de Planos de Manejo Florestal

361 SustentávelPMFS da Caatinga e suas formações sucessoras, e dá outras providências.

362 Instrução Normativa N° 1, de 25 de junho de 2009. Diário Oficial da União – Sec. 1, 120: 93.

363 Pagotto, M. A., Roig, F. A., Ribeiro, A. de S., and Lisi, C. S. 2015. Influence of regional

364 rainfall and Atlantic sea surface temperature on treering growth of Poincianella pyramidalis,

365 semiarid forest from Brazil. Dendrochronologia. 35 : 1423. doi:

366 10.1016/j.dendro.2015.05.007

367 Paula, J. E., 1993. Madeiras da Caatinga úteis para produção de energia. Pesquisa

368 Agropecuária Brasileira. 28 : 153165.

369 Peel, M. C., Finlayson, B. L., and McMahon, T. A. 2007. Updated world map of the Köppen

370 Geiger climate classification. Hydrology and Earth System Sciences. 11 : 1633–1644. doi:

371 10.5194/hess1116332007

15 https://mc06.manuscriptcentral.com/cjfr-pubs Canadian Journal of Forest Research Page 16 of 24

372 Quirino, W. F., Vale, A. T., Andrade, A. P. A., Abreu, V. L. S., and Azevedo, A. C. S. 2004.

373 Poder Calorífico da Madeira e de Resíduos Lignolcelulósicos. Biomassa & Energia. 1: 173

374 182.

375 Prodan, M. 1968. Forest Biometrics. Pergamon, Oxford. doi: 10.1016/B978008012441

376 4.500018.

377 Rede de Manejo Florestal da Caatinga. 2005. Protocolo de medições de parcelas permanentes

378 / Comitê Técnico Científico. Associação Plantas do Nordeste, Recife, 21 p.

379 Rodal, M. J. N., Sampaio, E. V. S. B., and Figueiredo, M. A. 1992. Manual sobre métodos de

380 estudo florístico e fitossociológico ecossistema Caatinga. Sociedade Botânica do Brasil,

381 Brasília, 24 p.

382 Rode, R., Leite, G. L., Silva, M. L., Ribeiro, C. A. A. S., and Binoti, D. H. B. 2014. The

383 economics and optimal managementDraft regimes of eucalyptus plantations: A case study of

384 forestry outgrower schemes in Brazil. Forest Policy and Economics. 44 : 2633. doi:

385 10.1016/j.forpol.2014.05.001

386 Rowell, R. M. 2012. Handbook of wood chemistry and wood composites. CRC press, Boca

387 Raton, 703 p.

388 Santana, J. A. S., and Souto, J. S. 2006. Diversidade e estrutura fitossociológica da Caatinga

389 na Estação Ecológica do SeridóRN. Revista de Biologia e Ciências da Terra. 6: 232242.

390 Santana, O. A. 2016. Resistência social na Caatinga árida: a narrativa de quem ficou no

391 colapso ambiental. Desenvolvimento e Meio Ambiente. 38 : 419438. doi:

392 http://dx.doi.org/10.5380/dma.v38i0.43574

393 Santana, O. A., and Encinas, J. I. 2016. Dendrophysiological plant strategies of Poincianella

394 pyramidalis (Tul.) L.P. Queiroz after wood herbivory in semiarid region of Paraíba Brazil.

395 Acta Scientiarum. Biological Sciences. 38 : 179186. doi: 10.4025/actascibiolsci.v38i2.29089

16 https://mc06.manuscriptcentral.com/cjfr-pubs Page 17 of 24 Canadian Journal of Forest Research

396 Schneider, P. R. 1998. Análise de regressão aplicada à engenharia florestal. UFSM/CEPEF,

397 Santa Maria, 236 p.

398 Schinker, M. G., Hansen, N., Spiecker H. 2003. Highfrequency densitometry — a new

399 method for the rapid evaluation of wood density variations. IAWA Journal. 24 : 231–239. doi:

400 10.1163/2294193290001592

401 Shchupakivskyy, R., Clauder, L., Linke, N., and Pfriem, A. 2014. Application of high

402 frequency densitometry to detect changes in early and latewood density of oak ( Quercus

403 robur L.) due to thermal modification. European Journal of Wood and Wood Products. 72 : 5

404 10. doi: 10.1007/s001070130744x

405 Scheitera, S., and Savadogo, P. 2016. Ecosystem management can mitigate vegetation shifts

406 induced by climate change in West Africa. Ecological Modelling. 332 : 19–27. doi:

407 10.1016/j.ecolmodel.2016.03.022 Draft

408 Schmiedel, U., Araya, Y., Bortolotto, M. I., Boeckenhoff, L., Hallwachs, W., Janzen, D.,

409 Kolipaka, S. S., Novotny, V., Palm, M., Parfondry, M., Smanis, A., and Toko, P. 2016.

410 Contributions of paraecologists and parataxonomists to research, conservation, and social

411 development. Conservation Biology. 30 : 506–519. doi: 10.1111/cobi.12661

412 Seijo, M. M., Huerta, R. P., Torné, J. M., Torné, C. M., and Vidal, E. A. 2016. Madera

413 Carbonizada en contextos funerários de la jefatura de Río Grande Panamá: Antracología en el

414 sitio de el Caño. Chungará. 48 : 277294. doi: 10.4067/S071773562016005000013

415 Silva, G. C., and Sampaio, E. V. S. B. 2008. Biomassas de partes aéreas em plantas de

416 caatinga. Revista Árvore. 32 : 567575. doi: 10.1590/S010067622008000300017

417 Sindusgesso. 2001. Newsletter [online]. Available from .

418 [accessed 13 June 2014].

17 https://mc06.manuscriptcentral.com/cjfr-pubs Canadian Journal of Forest Research Page 18 of 24

419 Silva, J. A. A. 2008/2009. Potencialidades de florestas energéticas de Eucalyptus no Pólo

420 Gesseiro do AraripePernambuco. Anais da Academia Pernambucana de Ciência

421 Agronômica. 5-6: 301319.

422 Soares, C. P. B., Paula Neto, F., and Souza, A. L. 2011. Dendrometria e inventário florestal.

423 Editora UFV, Viçosa, 272p.

424 Vico, G., Thompson, S. E., Manzoni, S., Molini, A., Albertson, J. D., Almeida‐Cortez, J. S.,

425 Fay, P. A., Feng, X., Guswa, A. J., Liu, H., Wilson, T. G., and Porporato, A. 2015. Climatic,

426 ecophysiological, and phenological controls on plant ecohydrological strategies in seasonally

427 dry ecosystems. Ecohydrology. 8, 660681. doi: 10.1002/eco.1533

428 Xiao, Q., and Huang, M. 2016. Fine root distributions of shelterbelt trees and their water

429 sources in an oasis of arid northwestern China. Journal of Arid Environments. 130 : 3039.

430 doi: 10.1016/j.jaridenv.2016.03.004 Draft

431 Zar, J. 1999. Biostatistical analysis. Prentice Hall, New Jersey.

18 https://mc06.manuscriptcentral.com/cjfr-pubs Page 19 of 24 Canadian Journal of Forest Research

Draft

Figure 1. Map showing the origin of the tree individuals studied (A); and methodological steps: volumetry (B); densitometry and dating (C) and calorimetry (D).

714x851mm (72 x 72 DPI)

https://mc06.manuscriptcentral.com/cjfr-pubs Canadian Journal of Forest Research Page 20 of 24

Draft

Figure 2. Age (A), volume (B), density (C) and lower heating value (LHV) (D) of twelve native species analyzed in the semi-arid region.

1103x1446mm (72 x 72 DPI)

https://mc06.manuscriptcentral.com/cjfr-pubs Page 21 of 24 Canadian Journal of Forest Research

Draft

Figure 3. Relationship between the analyzed variables: volume, lower heating value (LHV), density and age of the 240 woody individuals measured.

1431x1577mm (72 x 72 DPI)

https://mc06.manuscriptcentral.com/cjfr-pubs Canadian Journal of Forest Research Page 22 of 24

Table 1. Mathematical models used for data fit.

Model Equation*

Polynomial (linear) Yi=β 0 + β 1i ⋅ X + ε i

2 Polynomial (Quadratic) Yi=ββ 01i +⋅+⋅ X β 2i X + ε i

2 3 Polynomial (Cubic) Yi=+⋅+⋅ββ 01i2i XX β +⋅ β 3i X + ε i

β1 Yi= + ε i Sigmoidal (3 parameters) X1−β 2  −  1+ e β3 

β1 Yi=β 0 + + ε i Sigmoidal (4 parameters) Xi−β 2  −  1+ e β3 

β1 Yi=β 0 +β + ε i X −β  4 Sigmoidal (5 parameters) −i 2   1+ e β3    

β1 Yi=β + ε i Logistic (3 parameters) X  3 Draft1+ i  β2 

β1 Yi=β 0 +β + ε i Logistic (4 parameters) X  3 1+ i  β2 

β 1  4  X−β + β ln2 β4  −i 2 3   β  5     Weibull (4 parameters) Yi=β 1  1 − e  + ε i      

β 1  4  X−β + β ln2 β4  −i 2 3   β  5     Weibull (5 parameters) Yi=+−β 0 β 1  1 e  + ε i      

Xi−β 2    Gompertz (3 parameters) −eβ3  Yi=β 1 e + ε i

Xi−β 2    Gompertz (4 parameters) −eβ3  Yi01=β + β e + ε i

β X β2 Hill (3 parameters) Y =1 i + ε iβ2 β 2 i β3+ X i

1 https://mc06.manuscriptcentral.com/cjfr-pubs Page 23 of 24 Canadian Journal of Forest Research

β X β2 Hill (4 parameters) Y =β +1 i + ε i 0β2 β 2 i β3+ X i

β3 Chapman (3 parameters) −β2X i Yi=β 1( 1 − e ) + ε i

β3 Chapman (4 parameters) −β2X i Yi01=+β β( 1 − e ) + ε i

β1⋅X i Exponential Growth (Simple 1 parameter) Yi= e + ε i

β2⋅X i Exponential Growth (Simple 2 parameters) Yi=β 1 ⋅ e + ε i

β2⋅X i Exponential Growth (Simple 3 parameters) Yi01=β +⋅ β e + ε i

β2i⋅X β 4i ⋅ X Exponential Growth (Double 1 parameter) Yi1=⋅β e +⋅ β 3 e + ε i

β2i⋅X β 4i ⋅ X Exponential Growth (Double 2 parameter) Yi01=+⋅ββ e +⋅ β 3 e + ε i

*Y i = observed value of the dependent variable; X i = observed value of the independent variable; βi = regression equation

parameters; εi = fit error. Draft

2 https://mc06.manuscriptcentral.com/cjfr-pubs Canadian Journal of Forest Research Page 24 of 24

Table 2. Significant fit of data to models for six analyzed relationships: volume vs. density (A); volume vs. age (B); lower heating value vs. density (C); density vs. age (D); lower heating value vs. volume (E); and lower heating value vs. age (F).

Model Function R2 RMSE p

0.509 Sigmoidal (5 Yi = − 0.12 + 0.005 Xi − 0.423  A −   0.86 0.044 < 0.001 0.001  parameters) 1+ e   

0.031 Sigmoidal (5 Yi = 0.007 + 3.799 Xi − 23.861  B −   0.89 0.037 < 0.001 4.287  parameters) 1+ e   

Polynomial

C Yi= 3.736 + 1.738 ⋅ X i 0.97 0.026 < 0.001 (linear)

Exponential

0.001⋅ Xi 0.019 ⋅ X i D Growth (Double 5 Yi =− 0.193 + 0.248 ⋅ e + 0.210 ⋅ e 0.91 0.034 < 0.001

parameters)

Chapman (4 6.438 E Y= 0.004 + 0.001Draft 1 − e −6.297 ⋅ X i 0.90 0.048 < 0.001 i ( ) parameters)

Exponential

0.001⋅ Xi 0.044 ⋅ X i F Growth (Double 5 Yi =− 1415.94 + 5142.98 ⋅ e + 76.118 ⋅ e 0.94 0.050 < 0.012

parameters)

Source: Author.

3 https://mc06.manuscriptcentral.com/cjfr-pubs