Acessibilidade / Reportar erro

Fertilization of Eucalyptus Stands at Advanced Ages in Minas Gerais, Brazil

ABSTRACT

This study evaluated the response to fertilization at advanced ages in Eucalyptus stands. Fertilization at advanced ages (FA) was applied in stands managed according to tall shaft and coppice regimes. All stands received operational fertilization from the company. It was used as the basis to the network of continuous forestry inventory plots (P-IFC) of the company. Connected to the P-IFC, twin-plots received FA throughout four age classes, 30 to 35, 42 to 47, 54 to 59 and 66 to 70 months. The FA presented a biological response and provided volumetric gains of 17 and 51% for stands with tall shaft and coppice, respectively. The most significant responses to the FA was seen at younger ages of the forest in the tall shaft regime. Despite the biological response, the volumetric gains did not provide financial returns in either regime due to the elevated doses applied.

Keywords:
silviculture; forest nutrition; forest management

1. INTRODUCTION

In recent decades in Brazil the average gains in productivity of wood from 10 to 39 m3ha–1year–1 in Eucalyptus stands have been observed. These gains are mainly attributed to an evolution in forestry practices and to genetic improvements (Gonçalves et al., 2008Gonçalves JLM, Stape JL, Laclau JP, Bouillet JP, Ranger J. Assessing the effects of early silvicultural management on long-term site productivity of fastgrowing eucalypt plantations: the Brazilian experience. Southern Forests 2008; 70(2): 105-118. http://dx.doi.org/10.2989/SOUTH.FOR.2008.70.2.6.534.
http://dx.doi.org/10.2989/SOUTH.FOR.2008...
). According to Stape et al. (2010)Stape JL, Binkley D, Ryan MG, Fonseca S, Loos RA, Takahashi EN et al. The Brazil eucalyptus potential productivity project: Influence of water, nutrients and stand uniformity on wood production. Forest Ecology and Management 2010; 259(9): 1684-1694. http://dx.doi.org/10.1016/j.foreco.2010.01.012.
http://dx.doi.org/10.1016/j.foreco.2010....
, the use of mineral fertilizers corresponds, on average, to 30% of the gains in productivity.

Fertilization in Eucalyptus plantations is normally carried out between the pre-planting phase and the canopy closure. After two years of age, fertilization is reduced or stopped completely (Silva et al., 2013Silva PHM, Poggiani F, Libardi PL, Gonçalves AN. Fertilizer management of eucalypt plantations on sandy soil in Brazil: Initial growth and nutrient cycling. Forest Ecology and Management 2013; 301: 67-78. http://dx.doi.org/10.1016/j.foreco.2012.10.033.
http://dx.doi.org/10.1016/j.foreco.2012....
). A direct relation between biomass production and nutrient accumulation was observed by Santana et al. (2008)Santana RC, Barros NF, Novais RF, Leite HG, Comerford NB. Alocação de nutrientes em plantios de eucalipto no Brasil. Revista Brasileira de Ciência do Solo 2008; 32(spe): 2723-2733. http://dx.doi.org/10.1590/S0100-06832008000700016.
http://dx.doi.org/10.1590/S0100-06832008...
. These authors observed that 70% of the total nutrient demand of Eucalyptus plantations in Brazil is accumulated up until 4.5 years of age. Between planting and the maximum leaf area index, the stand shows high levels of nutrient absorption and elevated metabolic activity, factors that explain the greater response to fertilization at a younger age (Gonçalves et al., 2013Gonçalves JLM, Alvares CA, Higa AR, Silva LD, Alfenas AC, Stahl J et al. Integrating genetic and silvicultural strategies to minimize abiotic and biotic constraints in Brazilian eucalypt plantations. Forest Ecology and Management 2013; 301: 6-27. http://dx.doi.org/10.1016/j.foreco.2012.12.030.
http://dx.doi.org/10.1016/j.foreco.2012....
). After the canopy closure phase, the biochemical cycle becomes gradually more established and tends to meet part of the nutritional demand of the stand. The nutrient cycling associated with the fertilization applied during the initial years is considered adequate to satisfy the total demand of the stand because soluble and relatively insoluble fertilizers are used with the intention of providing nutrients throughout the whole cycle. Consequently, it is expected that the plantations are not responsive to the application of nutrients at advanced ages.

Despite this, studies seeking to analyze the nutritional limitations and potential productivity of Eucalyptus have shown a response to fertilization at ages older than two years when very elevated doses of nutrients are applied (Stape et al., 2006Stape JL, Binkley D, Jacob WS, Takahashi EN. Twin-plot approach to determine nutrient limitation and potential productivity in Eucalyptus plantations at landscape scales in Brazil. Forest Ecology and Management 2006; 223(1-3): 358-362. http://dx.doi.org/10.1016/j.foreco.2005.11.015.
http://dx.doi.org/10.1016/j.foreco.2005....
; Ferreira & Stape, 2010Ferreira JMA, Stape JL. Productivity gains by fertilisation in Eucalyptus urophylla clonal plantations across gradients in site and stand conditions. Southern Forests: a Juornal of Forest Science 2010; 71(4): 253-258. http://dx.doi.org/10.2989/SF.2009.71.4.1.1028.
http://dx.doi.org/10.2989/SF.2009.71.4.1...
; Silva et al., 2016Silva RML, Hakamada RE, Bazani JH, Otto MSG, Stape JL. Fertilization response, light use, and growth efficiency in Eucalyptus plantations across Soil and climate gradients in Brazil. Forests 2016; 7(12): 117. http://dx.doi.org/10.3390/f7060117.
http://dx.doi.org/10.3390/f7060117...
). According to the cited authors, the increase in productivity due to fertilization at advanced ages is influenced, amongst other factors, by the stage of development of the forest, soil and climatic conditions, genetic material and type of forestry management adopted. The knowledge of the probability of response to fertilization and its economic viability (Ferreira & Stape, 2010Ferreira JMA, Stape JL. Productivity gains by fertilisation in Eucalyptus urophylla clonal plantations across gradients in site and stand conditions. Southern Forests: a Juornal of Forest Science 2010; 71(4): 253-258. http://dx.doi.org/10.2989/SF.2009.71.4.1.1028.
http://dx.doi.org/10.2989/SF.2009.71.4.1...
) and the classification of responsive and non-responsive sites (Fisher & Binkley, 2000Fisher RF, Binkley D. Ecology and management of forest soils. Forests 2000; 6: 973-991.) are important factors to be considered in making forestry management decisions.

Different to management by tall shaft regime, there have been very little studies on the responses to fertilization in coppice regimes. Generally, the initial fertilization of the stand in coppice regimes takes place after sprout thinning between eight and eighteen months post-harvest. During the phase prior to fertilization, the growth of the stand depends on the nutrient cycling of forest debris, ground cover and the nutrients available in the soil. Reports regarding a drop in productivity during the second rotation in comparison with the first are common (Miranda et al., 1998aMiranda GA, Barros NF, Leite HG, Couto L, Teixeira JL. Produção de povoamentos de eucalipto em regime de talhadia, em função da adubação e da distribuição da galhada da rotação anterior, no vale do Jequitinhonha – MG. Revista Árvore 1998a; 22(3): 307-314. , bMiranda GA, Barros NF, Leite HG, Couto L, Filho MN. Produtividade de povoamento de eucalipto em regime de talhadia, em função da adubação parcelada no Vale do Jequitinhonha- MG. Revista Árvore 1998b; 22(1): 51-59.). Part of the loss of productivity of the stands in the second rotation is attributed to the management adopted in sprout thinning and later fertilization.

Given this, this study aims to evaluate the response to fertilization at advanced ages of Eucalyptus stands under tall shaft and coppice regimes in the state of Minas Gerais.

2. MATERIALS AND METHODS

The present study was realized in clonal stands of Eucalyptus, planted in six regions (R1; R2; R3; R4; R5; and R6) in the state of Minas Gerais (Table 1).

Table 1
The main characteristics of the study area.

To evaluate the potential for response to fertilization at advanced ages (FA), 142 twin-plots were planted (Stape et al., 2006Stape JL, Binkley D, Jacob WS, Takahashi EN. Twin-plot approach to determine nutrient limitation and potential productivity in Eucalyptus plantations at landscape scales in Brazil. Forest Ecology and Management 2006; 223(1-3): 358-362. http://dx.doi.org/10.1016/j.foreco.2005.11.015.
http://dx.doi.org/10.1016/j.foreco.2005....
) connected to the 142 plots of the network of the continuous forestry inventory (P-IFC). The plots chosen represented stands in four age classes (CLI). The first was between 30 and 35 months, the second between 42 and 47 months, the third between 54 and 59 months, and the fourth between 66 and 70 months. The study evaluated 244 circular parcels, with 60 trees per parcel, planted in spaces varying from 4.9 to 9.5 m2 per plant, with 15 clones under two management regimes, tall shaft and coppice (Table 2).

Table 2
Spacing, clone, age class (CLI) and number of plots per region and management regime.

The twin-plots were planted between March and May of 2011, set in five rows at the end of the continuous forestry inventory, with the same soil type and slope.

The evaluated treatments were organized in the following way: Control (C) – operational fertilization (Table 3) represented by the plots of continuous forestry inventory, which received the best forestry practices of the company; fertilization at advanced ages (FA) – recommended to eliminate possible nutritional limitations, applied on the twin-plots that had already received the operational fertilization, in stands in tall shaft and coppice regimes, and in the CLI 1, 2, 3 and 4, according to the proposal of Stape et al. (2006)Stape JL, Binkley D, Jacob WS, Takahashi EN. Twin-plot approach to determine nutrient limitation and potential productivity in Eucalyptus plantations at landscape scales in Brazil. Forest Ecology and Management 2006; 223(1-3): 358-362. http://dx.doi.org/10.1016/j.foreco.2005.11.015.
http://dx.doi.org/10.1016/j.foreco.2005....
. The FA was divided into four applications, the first (FA1) in the establishment of the twin-plots and the others 6 (FA2), 12 (FA3), and 18 (FA4) months later (Table 4).

Table 3
Average and standard deviation of the nutrients applied in the operational fertilization by region.
Table 4
Nutrients applied in the fertilization at advanced age (FA).

The trees of the plots (C and FA) were measured in the twin-plots installations every six months subsequently, between 2011 and 2013. In each plot, the diameter at 1.30 m from the ground was measured (DBH) for all 60 trees as well as the total height (HT) of the first seven trees. The height of the remaining 53 trees was estimated using artificial neural networks – ANNs (Özçelik et al., 2013Özçelik R, Diamantopoulou MJ, Crecente-Campo F, Eler U. Estimating crimean juniper tree height using nonlinear regression and artificial neural network models. Forest Ecology and Management 2013; 306: 52-60. http://dx.doi.org/10.1016/j.foreco.2013.06.009.
http://dx.doi.org/10.1016/j.foreco.2013....
).

To estimate the ANNs, the quantitative entry variables were considered to be age (years) and DBH (cm), the categorical entry variables were the rotation and genetic material and the output variable was the HT (m). The trained networks were of the multiple layer perceptron type, also known as MLP (Multilayer Perceptron). An ANN was trained for each treatment (C and FA) utilizing the software Statistica 10 (Statsoft, 2010Statsoft. Statistica – data analysis software system. Version 10.0 [online]. Palo Alto: Statsoft Inc; 2010 [cited 2017 Apr 25]. Available from: http://www.statsoft.com
http://www.statsoft.com...
), which used 50% of the data to train the networks, 25% for the test and 25% for validation.

Based on the DBH and HT values, the individual volumes were estimated (V – Equation 1) for each tree, in each parcel evaluated, utilizing the Shumacher & Hall (1933)Shumacher FX, Hall FS. Logarithmic expression of timber-tree volume. Journal of Agricultural Research 1933; 7(9): 719-734. model,

L n V = β 0 + β 1 * L n D B H + β 2 * ln H T + ε (1)

in which: Ln is the Naperian logarithm; V the volume per tree (m3); β0, β1 and β2 the model parameters; DBH the diameter at 1.3 m from the ground (cm); HT the total height of the trees (m); and ε the random error.

To obtain the initial twin character between the plots C and FA, plots whose volumes were statistically equal according to T-paired test analysis were chosen, with a level of significance of 0.05 (Figure 1) (Stape et al. 2006Stape JL, Binkley D, Jacob WS, Takahashi EN. Twin-plot approach to determine nutrient limitation and potential productivity in Eucalyptus plantations at landscape scales in Brazil. Forest Ecology and Management 2006; 223(1-3): 358-362. http://dx.doi.org/10.1016/j.foreco.2005.11.015.
http://dx.doi.org/10.1016/j.foreco.2005....
).

Figure 1
The twin character between the initial volumes of the control (C) and the fertilization at an advanced age (FA) in the stands in tall shaft (A) and coppice (B) regimes.

In each treatment (C and FA) the volumetric gain (GV – Equation 2), was expressed by the difference between the volume obtained in the last and first measurements:

G V = V f i V 0 i (2)

In which: GV is the volumetric gain (m3 ha–1); V the volume (m3 ha–1); f the last measurement; 0 the first measurement; and i the treatment.

The response to FA, independent of R and CLI, was evaluated by comparison with the GV between the treatments C and FA. For this, the T-paired test was applied, with a level of significance of 0.05 utilizing the Statistic 10 software (Statsoft, 2010Statsoft. Statistica – data analysis software system. Version 10.0 [online]. Palo Alto: Statsoft Inc; 2010 [cited 2017 Apr 25]. Available from: http://www.statsoft.com
http://www.statsoft.com...
). The response to FA (Equation 3) by R and CLI was evaluated using descriptive statistics, obtaining an average of the GV for each treatment (C and FA) and comparing them according to the relative difference between the average of GV of the FA in relation to the average of the volumetric gain of C:

R F A = G V F A G V C (3)

in which: RFA is the response to FA in m3 ha–1; and GV is the volumetric gain (m3 ha–1).

To evaluate the economic viability of the GV stemming from the FA, a cost benefit analysis according to Kimberley et al. (2004)Kimberley MO, Wang H, Wilks PJ, Fisher CR, Magesan GN. Economic analysis of growth response from a pine plantation forest applied with biosolids. Forest Ecology and Management 2004; 186(1-3): 345-351. http://dx.doi.org/10.1016/j.foreco.2003.09.003.
http://dx.doi.org/10.1016/j.foreco.2003....
was carried out for two situations: general analysis (independent of R and CLI) and analysis by CLI. The costs per hectare included expenses for fertilizer and labor (F1 – R$ 3,618.9; F2 – R$ 2,090.0; F3 – R$ 2,090; F4 – R$ 2,090). For the calculation of the income, the value of sales of standing wood was considered to be R$ 55.0 m–3. The current values of the costs were deflated according to the General Price Index – internally available at the Fundação Getúlio Vargas (FGV, 2017Fundação Getúlio Vargas – FGV [online]. 2017 [cited 2017 Apr 25]. Available from: http://www.fgv.br
http://www.fgv.br...
) as suggested by Nascimento et al. (2016)Nascimento V, Arf O, Tarsitano MAA, Garcia NFS, Penteado MS, Sabundjian MT. Economic analysis of rice and common bean production in succession to green manure crops and mecahnical soil decompression in Brazilian Savannah. Revista Ceres 2016; 63(3): 315-322. http://dx.doi.org/10.1590/0034-737X201663030006.
http://dx.doi.org/10.1590/0034-737X20166...
.

The probability of the response to FA, by the GV, was analyzed using five classes according to the frequency of response by plot. The class (CL) was distributed according to the average (X¯) and the standard deviation (S): CL1 – R < X¯2S; CL2 – X¯2S<R<X¯1S ; CL3 – X¯1S<R< X¯+1S ; CL4 – X¯+1S<R< X¯+2S; and CL5 – R > X¯+2S.

3. RESULTS

Stands in a tall shaft regime presented, respectively, volumes (m3 ha–1) initial averages and mean annual increment – MAI (m3 ha–1 year–1) – of 128 and 34 for C and 129 and 36 for FA, and for the coppice regime 73 and 26 for C and 76 and 33 for FA.

The FA provided a positive response in terms of volumetric production in both management regimes (Figure 2). The volumetric gain, independent of age and region, was 17% greater in the tall shaft (FA = 77 m3 ha–1; C = 66 m3 ha–1) and 51% in the coppice (FA = 95 m3 ha–1; C = 63 m3 ha–1).

Figure 2
Relationship between the volumetric gains of the control (C) and fertilization at advanced age (FA) in stands with tall shaft (A) and coppice (B) regimes.

Considering the tall shaft regime, the greatest gain was observed in R6 (39%) and the least in R2 (6%) (Table 5). The response to the FA was different depending on the age at the time of application. Stands in the CLI 1 showed volumetric gains of 21% and those in the CL4 of 11%. In the coppice regime, gains were above those of the tall shaft regime with values reaching 57% in the R5 and 62% in the CLI 1.

Table 5
Average volumetric gains of the control (C) and the fertilization at advanced age (FA) for stands in tall shaft and coppice regimes.

The cost for the FA was the same in all evaluated situations (Table 6). It can be observed that the GV with FA did not present an economic return for stands in tall shaft and coppice regimes, independent of R and CLI, given that income obtained was much less than the cost. In analyses only using CLI, it can be observed that the stands that received the FA at younger ages were those presenting the greatest additional incomes.

Table 6
Cost of fertilization at advance age (FA) and income obtained from volumetric gains derived from FA in stands with tall shaft and coppice regimes.

The probability of the response to FA tended toward normal distribution in the stands in tall shaft and coppice regimes (Figure 3). In the tall shaft regime, the average gain was 11 m3 ha–1 with a standard deviation of 22 m3 ha–1. In percentage terms, 38% of the plots presented responses above the average and 19% presented more than twice the average. In the coppice regime, the average was 32 m3 ha–1 with a deviation of 15 m3 ha–1. In this regime, 37% of the plots showed responses above the average and 7% were more than twice the average.

Figure 3
Distribution of the response to fertilization at advanced age of all the parcels, by response class (R): CL1 – R < X¯2S; CL2 – X¯2S<R<X¯1S ; CL3 – X¯1S<R< X¯+1S ; CL4 – X¯+1S<R< X¯+2S; and CL5 – R > X¯+2S, for the stands in tall shaft (A) and coppice (B) regimes.

4. DISCUSSION

The response to FA in terms of GV showed an average of 11 m3 ha–1 in tall shaft regimes (Figure 2). Studies carried out in the state of São Paulo agree with these observations. The average gains of 10 and 8.1 m3 ha–1 year–1 were observed, respectively, by Stape et al. (2006)Stape JL, Binkley D, Jacob WS, Takahashi EN. Twin-plot approach to determine nutrient limitation and potential productivity in Eucalyptus plantations at landscape scales in Brazil. Forest Ecology and Management 2006; 223(1-3): 358-362. http://dx.doi.org/10.1016/j.foreco.2005.11.015.
http://dx.doi.org/10.1016/j.foreco.2005....
and Ferreira & Stape (2010)Ferreira JMA, Stape JL. Productivity gains by fertilisation in Eucalyptus urophylla clonal plantations across gradients in site and stand conditions. Southern Forests: a Juornal of Forest Science 2010; 71(4): 253-258. http://dx.doi.org/10.2989/SF.2009.71.4.1.1028.
http://dx.doi.org/10.2989/SF.2009.71.4.1...
. In both studies, the same FA utilized in the present study was used. The operational fertilization utilized was based on the nutritional balance method. This method was developed to achieve desired productivity and guarantee the sustainability of the forestry sites (Gonçalves et al., 2008Gonçalves JLM, Stape JL, Laclau JP, Bouillet JP, Ranger J. Assessing the effects of early silvicultural management on long-term site productivity of fastgrowing eucalypt plantations: the Brazilian experience. Southern Forests 2008; 70(2): 105-118. http://dx.doi.org/10.2989/SOUTH.FOR.2008.70.2.6.534.
http://dx.doi.org/10.2989/SOUTH.FOR.2008...
) with fertilization generally focusing on the first two years of life to guarantee growth until nutrient cycling was established.

Despite the biological response observed with FA, it did not lead to economic viability, given that the additional income derived from GV presented values well below the cost of the fertilizer and the labor for its application (Table 6). Considering the cost benefit analysis by CLI, it was found that the additional income was greater in stands that received FA at a younger age (CLI 1). This fact agrees with the observations of Santana et al. (2008)Santana RC, Barros NF, Novais RF, Leite HG, Comerford NB. Alocação de nutrientes em plantios de eucalipto no Brasil. Revista Brasileira de Ciência do Solo 2008; 32(spe): 2723-2733. http://dx.doi.org/10.1590/S0100-06832008000700016.
http://dx.doi.org/10.1590/S0100-06832008...
and Gonçalves et al. (2013)Gonçalves JLM, Alvares CA, Higa AR, Silva LD, Alfenas AC, Stahl J et al. Integrating genetic and silvicultural strategies to minimize abiotic and biotic constraints in Brazilian eucalypt plantations. Forest Ecology and Management 2013; 301: 6-27. http://dx.doi.org/10.1016/j.foreco.2012.12.030.
http://dx.doi.org/10.1016/j.foreco.2012....
where a greater demand for nutrients during younger phases of forest development was demonstrated with a consequent greater potential for response to its application at younger ages. It is underlined that the economic returns from FA were not expected since it was primarily recommended to guard against nutritional shortfalls in the stands.

The average of the MAI in the coppice regime was 23% lower than the tall shaft regime for operational fertilization. However, a much greater biological response to FA in terms of the tall shaft regime was observed, despite not observing an economic benefit. Probably the unviable economic outcome occurred due to the elevated doses of fertilizers applied in the FA. The use of high doses aimed to avoid growth restrictions due to nutrient shortfalls. This indicates that studies with fertilizer doses in FA in the coppice regime have a greater potential to achieve economic viability.

The operational fertilization in the coppice regime is only carried out after sprout thinning, which occurs at 12 months post-harvest, applying quantities of nutrients well below those of the tall shaft regime (Table 3). This marked reduction of the doses applied, associated with stress caused by the sprout thinning, eliminating between 50% and 70% of the biomass produced, is one of the factors that could explain the drop in productivity and the greater biological response to FA. Less productivity tends to present more significant responses to fertilizations in Eucalyptus stands (Ferreira & Stape, 2010Ferreira JMA, Stape JL. Productivity gains by fertilisation in Eucalyptus urophylla clonal plantations across gradients in site and stand conditions. Southern Forests: a Juornal of Forest Science 2010; 71(4): 253-258. http://dx.doi.org/10.2989/SF.2009.71.4.1.1028.
http://dx.doi.org/10.2989/SF.2009.71.4.1...
; Silva et al., 2016Silva RML, Hakamada RE, Bazani JH, Otto MSG, Stape JL. Fertilization response, light use, and growth efficiency in Eucalyptus plantations across Soil and climate gradients in Brazil. Forests 2016; 7(12): 117. http://dx.doi.org/10.3390/f7060117.
http://dx.doi.org/10.3390/f7060117...
). Studies in stands managed according to a coppice regime have shown lower productivity in relation to the tall shaft regime (Miranda et al., 1998aMiranda GA, Barros NF, Leite HG, Couto L, Filho MN. Produtividade de povoamento de eucalipto em regime de talhadia, em função da adubação parcelada no Vale do Jequitinhonha- MG. Revista Árvore 1998b; 22(1): 51-59., bMiranda GA, Barros NF, Leite HG, Couto L, Teixeira JL. Produção de povoamentos de eucalipto em regime de talhadia, em função da adubação e da distribuição da galhada da rotação anterior, no vale do Jequitinhonha – MG. Revista Árvore 1998a; 22(3): 307-314.; Faria et al., 2002Faria GE, Barros NF, Novais RF, Lima JC, Teixeira JL. Produção e estado nutricional de povoamentos de Eucalyptus grandis, em segunda rotação, em resposta à adubação potássica. Revista Árvore 2002; 26(5): 577-584. http://dx.doi.org/10.1590/S0100-67622002000500008.
http://dx.doi.org/10.1590/S0100-67622002...
). The authors argue that part of the drop in productivity of these stands results from nutritional deficiencies generated by the exportation of nutrients. These nutritional limitations can be minimized if the harvest is restricted to the wood, leaving all the residue in the area (Santana et al., 2002Santana RC, Barros NF, Neves JCL. Eficiência de utilização de nutrientes e sustentabilidade da produção em procedências de Eucalyptus grandis e Eucalyptus saligna em sítios florestais do estado de São Paulo. Revista Árvore 2002; 26(4): 447-457. http://dx.doi.org/10.1590/S0100-67622002000400007.
http://dx.doi.org/10.1590/S0100-67622002...
).

When comparing the two extremes of response to the same FA by region, in the tall shaft regime we can infer that greater water availability favors production (Tables 1 and 5). To demonstrate this effect, regions R2 and R3 were selected, since they possess all the age classes evaluated. A greater response to fertilization in regions with greater water availability is widely reported in the literature (Smethurst et al., 2003Smethurst P, Baillie C, Cherry M, Holz G. Fertilizer effects on LAI and growth of four Eucalyptus nitens plantations. Forest Ecology and Management 2003; 176(1-3): 531-542. http://dx.doi.org/10.1016/S0378-1127(02)00226-8.
http://dx.doi.org/10.1016/S0378-1127(02)...
; Stape et al., 2010Stape JL, Binkley D, Ryan MG, Fonseca S, Loos RA, Takahashi EN et al. The Brazil eucalyptus potential productivity project: Influence of water, nutrients and stand uniformity on wood production. Forest Ecology and Management 2010; 259(9): 1684-1694. http://dx.doi.org/10.1016/j.foreco.2010.01.012.
http://dx.doi.org/10.1016/j.foreco.2010....
; Silva et al., 2016Silva RML, Hakamada RE, Bazani JH, Otto MSG, Stape JL. Fertilization response, light use, and growth efficiency in Eucalyptus plantations across Soil and climate gradients in Brazil. Forests 2016; 7(12): 117. http://dx.doi.org/10.3390/f7060117.
http://dx.doi.org/10.3390/f7060117...
).

Greater responses to the FA was observed during development in forests with tall shaft regimes, in other words, the younger the age, the greater the response observed. Stands with 2.5 years (CLI 1) are nearer to the canopy closure stage. This stage tends to show a greater leaf area index (Smethurst et al., 2003Smethurst P, Baillie C, Cherry M, Holz G. Fertilizer effects on LAI and growth of four Eucalyptus nitens plantations. Forest Ecology and Management 2003; 176(1-3): 531-542. http://dx.doi.org/10.1016/S0378-1127(02)00226-8.
http://dx.doi.org/10.1016/S0378-1127(02)...
; Laclau et al., 2010Laclau J-P, Ranger J, Gonçalves JLM, Maquère V, Krusche AV, M’Bou AT et al. Biogeochemical cycles of nutrients in tropical Eucalyptus plantations main features shown by intensive monitoring in Congo and Brazil. Forest Ecology and Management 2010; 259(9): 1771-1785. http://dx.doi.org/10.1016/j.foreco.2009.06.010.
http://dx.doi.org/10.1016/j.foreco.2009....
), consequently, there is a greater demand for nutrients. The lower response to FA at more advanced ages can be attributed to the establishment of nutrient cycling (Santana et al., 2008Santana RC, Barros NF, Novais RF, Leite HG, Comerford NB. Alocação de nutrientes em plantios de eucalipto no Brasil. Revista Brasileira de Ciência do Solo 2008; 32(spe): 2723-2733. http://dx.doi.org/10.1590/S0100-06832008000700016.
http://dx.doi.org/10.1590/S0100-06832008...
; Gonçalves et al., 2013Gonçalves JLM, Alvares CA, Higa AR, Silva LD, Alfenas AC, Stahl J et al. Integrating genetic and silvicultural strategies to minimize abiotic and biotic constraints in Brazilian eucalypt plantations. Forest Ecology and Management 2013; 301: 6-27. http://dx.doi.org/10.1016/j.foreco.2012.12.030.
http://dx.doi.org/10.1016/j.foreco.2012....
) and of operational fertilization that adopted the nutritional balance model where one does not expect nutritional restrictions. Some authors report greater responses to FA at advanced ages in studies with twin-plots (Stape et al., 2006Stape JL, Binkley D, Jacob WS, Takahashi EN. Twin-plot approach to determine nutrient limitation and potential productivity in Eucalyptus plantations at landscape scales in Brazil. Forest Ecology and Management 2006; 223(1-3): 358-362. http://dx.doi.org/10.1016/j.foreco.2005.11.015.
http://dx.doi.org/10.1016/j.foreco.2005....
; Ferreira & Stape, 2010Ferreira JMA, Stape JL. Productivity gains by fertilisation in Eucalyptus urophylla clonal plantations across gradients in site and stand conditions. Southern Forests: a Juornal of Forest Science 2010; 71(4): 253-258. http://dx.doi.org/10.2989/SF.2009.71.4.1.1028.
http://dx.doi.org/10.2989/SF.2009.71.4.1...
; Silva et al., 2016Silva RML, Hakamada RE, Bazani JH, Otto MSG, Stape JL. Fertilization response, light use, and growth efficiency in Eucalyptus plantations across Soil and climate gradients in Brazil. Forests 2016; 7(12): 117. http://dx.doi.org/10.3390/f7060117.
http://dx.doi.org/10.3390/f7060117...
) where nutritional restrictions were observed using limited quantities of nutrients applied during operational fertilization.

The probability of the response to the FA can be employed to classify the stands with greater or lesser potential for response (Fisher & Binkley, 2000Fisher RF, Binkley D. Ecology and management of forest soils. Forests 2000; 6: 973-991.). In the present study, approximately 37% of the plots presented a greater potential for biological response to the FA in both regimes evaluated (Figure 3). This highlights that applying the same fertilization indiscriminately without considering regional quality may not generate the best outcomes. It is common in Brazilian Eucalyptus cultivation to apply similar quantities of nutrients in locations with distinct potentials in the response, a fact that can compromise economic returns.

The FA applied at advanced ages where operational fertilization was established employing the nutritional balance method shows that starting from 2.5 years of age there was no economic benefit, despite the biological response. To better exploit the potential of this type of fertilization, its application, taking into account different levels and dosages in a sequence of ages between six and 30 months, can demonstrate the viability or not of the response to fertilization at advanced ages.

5. CONCLUSION

  • Fertilization at advanced age presented a biological response and provided volumetric gains for stands managed according to tall shaft and coppice regimes.

  • Stands with a coppice regime possess greater potential for biological response to fertilization at an advanced age.

  • Despite the biological response, the volumetric gains did not provide a financial return in either regime due to the elevated doses applied.

ACKNOWLEDGEMENTS

We would like to thank Universidade Federal dos Vales do Jequitinhonha e Mucuri – UFVJM; Conselho Nacional de Pesquisa e Desenvolvimento – CNPq; Fundação de Amparo à Pesquisa de Minas Gerais – FAPEMIG; This study was financed in party by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brasil (CAPES). Gerdau Florestal S/A for the data base for this work.

  • FINANCIAL SUPORT Universidade Federal dos Vales do Jequitinhonha e Mucuri – UFVJM; Conselho Nacional de Pesquisa e Desenvolvimento – CNPq; Fundação de Amparo à Pesquisa de Minas Gerais – FAPEMIG; This study was financed in party by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brasil (CAPES) Finance Code 001. Gerdau Florestal AS for the data base for this work.

REFERENCES

  • Alvares CA, Stape JL, Sentelhas PC, Gonçalves JLM, Sparovek G. Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift (Berlin) 2013; 22(6): 711-728. http://dx.doi.org/10.1127/0941-2948/2013/0507
    » http://dx.doi.org/10.1127/0941-2948/2013/0507
  • Faria GE, Barros NF, Novais RF, Lima JC, Teixeira JL. Produção e estado nutricional de povoamentos de Eucalyptus grandis, em segunda rotação, em resposta à adubação potássica. Revista Árvore 2002; 26(5): 577-584. http://dx.doi.org/10.1590/S0100-67622002000500008
    » http://dx.doi.org/10.1590/S0100-67622002000500008
  • Ferreira JMA, Stape JL. Productivity gains by fertilisation in Eucalyptus urophylla clonal plantations across gradients in site and stand conditions. Southern Forests: a Juornal of Forest Science 2010; 71(4): 253-258. http://dx.doi.org/10.2989/SF.2009.71.4.1.1028
    » http://dx.doi.org/10.2989/SF.2009.71.4.1.1028
  • Fisher RF, Binkley D. Ecology and management of forest soils. Forests 2000; 6: 973-991.
  • Fundação Getúlio Vargas – FGV [online]. 2017 [cited 2017 Apr 25]. Available from: http://www.fgv.br
    » http://www.fgv.br
  • Gonçalves JLM, Alvares CA, Higa AR, Silva LD, Alfenas AC, Stahl J et al. Integrating genetic and silvicultural strategies to minimize abiotic and biotic constraints in Brazilian eucalypt plantations. Forest Ecology and Management 2013; 301: 6-27. http://dx.doi.org/10.1016/j.foreco.2012.12.030
    » http://dx.doi.org/10.1016/j.foreco.2012.12.030
  • Gonçalves JLM, Stape JL, Laclau JP, Bouillet JP, Ranger J. Assessing the effects of early silvicultural management on long-term site productivity of fastgrowing eucalypt plantations: the Brazilian experience. Southern Forests 2008; 70(2): 105-118. http://dx.doi.org/10.2989/SOUTH.FOR.2008.70.2.6.534
    » http://dx.doi.org/10.2989/SOUTH.FOR.2008.70.2.6.534
  • Instituto de Geociências Aplicadas – IGA [online]. 2015 [cited 2017 May 25]. Available from: http://www.iga.br/siteiga/mapas/cgi/Imagens/Localizacao_MG.pdf
    » http://www.iga.br/siteiga/mapas/cgi/Imagens/Localizacao_MG.pdf
  • Kimberley MO, Wang H, Wilks PJ, Fisher CR, Magesan GN. Economic analysis of growth response from a pine plantation forest applied with biosolids. Forest Ecology and Management 2004; 186(1-3): 345-351. http://dx.doi.org/10.1016/j.foreco.2003.09.003
    » http://dx.doi.org/10.1016/j.foreco.2003.09.003
  • Laclau J-P, Ranger J, Gonçalves JLM, Maquère V, Krusche AV, M’Bou AT et al. Biogeochemical cycles of nutrients in tropical Eucalyptus plantations main features shown by intensive monitoring in Congo and Brazil. Forest Ecology and Management 2010; 259(9): 1771-1785. http://dx.doi.org/10.1016/j.foreco.2009.06.010
    » http://dx.doi.org/10.1016/j.foreco.2009.06.010
  • Meneses VA, Trugilho PF, Calegario N, Leite HG. Efeito da idade e do sítio na densidade básica e produção de massa seca de madeira em um clone do Eucalyptus urophylla. Scientia Forestalis 2015; 43(105): 101-116.
  • Miranda GA, Barros NF, Leite HG, Couto L, Filho MN. Produtividade de povoamento de eucalipto em regime de talhadia, em função da adubação parcelada no Vale do Jequitinhonha- MG. Revista Árvore 1998b; 22(1): 51-59.
  • Miranda GA, Barros NF, Leite HG, Couto L, Teixeira JL. Produção de povoamentos de eucalipto em regime de talhadia, em função da adubação e da distribuição da galhada da rotação anterior, no vale do Jequitinhonha – MG. Revista Árvore 1998a; 22(3): 307-314.
  • Nascimento V, Arf O, Tarsitano MAA, Garcia NFS, Penteado MS, Sabundjian MT. Economic analysis of rice and common bean production in succession to green manure crops and mecahnical soil decompression in Brazilian Savannah. Revista Ceres 2016; 63(3): 315-322. http://dx.doi.org/10.1590/0034-737X201663030006
    » http://dx.doi.org/10.1590/0034-737X201663030006
  • Özçelik R, Diamantopoulou MJ, Crecente-Campo F, Eler U. Estimating crimean juniper tree height using nonlinear regression and artificial neural network models. Forest Ecology and Management 2013; 306: 52-60. http://dx.doi.org/10.1016/j.foreco.2013.06.009
    » http://dx.doi.org/10.1016/j.foreco.2013.06.009
  • Santana RC, Barros NF, Neves JCL. Eficiência de utilização de nutrientes e sustentabilidade da produção em procedências de Eucalyptus grandis e Eucalyptus saligna em sítios florestais do estado de São Paulo. Revista Árvore 2002; 26(4): 447-457. http://dx.doi.org/10.1590/S0100-67622002000400007
    » http://dx.doi.org/10.1590/S0100-67622002000400007
  • Santana RC, Barros NF, Novais RF, Leite HG, Comerford NB. Alocação de nutrientes em plantios de eucalipto no Brasil. Revista Brasileira de Ciência do Solo 2008; 32(spe): 2723-2733. http://dx.doi.org/10.1590/S0100-06832008000700016
    » http://dx.doi.org/10.1590/S0100-06832008000700016
  • Shumacher FX, Hall FS. Logarithmic expression of timber-tree volume. Journal of Agricultural Research 1933; 7(9): 719-734.
  • Silva PHM, Poggiani F, Libardi PL, Gonçalves AN. Fertilizer management of eucalypt plantations on sandy soil in Brazil: Initial growth and nutrient cycling. Forest Ecology and Management 2013; 301: 67-78. http://dx.doi.org/10.1016/j.foreco.2012.10.033
    » http://dx.doi.org/10.1016/j.foreco.2012.10.033
  • Silva RML, Hakamada RE, Bazani JH, Otto MSG, Stape JL. Fertilization response, light use, and growth efficiency in Eucalyptus plantations across Soil and climate gradients in Brazil. Forests 2016; 7(12): 117. http://dx.doi.org/10.3390/f7060117
    » http://dx.doi.org/10.3390/f7060117
  • Smethurst P, Baillie C, Cherry M, Holz G. Fertilizer effects on LAI and growth of four Eucalyptus nitens plantations. Forest Ecology and Management 2003; 176(1-3): 531-542. http://dx.doi.org/10.1016/S0378-1127(02)00226-8
    » http://dx.doi.org/10.1016/S0378-1127(02)00226-8
  • Stape JL, Binkley D, Jacob WS, Takahashi EN. Twin-plot approach to determine nutrient limitation and potential productivity in Eucalyptus plantations at landscape scales in Brazil. Forest Ecology and Management 2006; 223(1-3): 358-362. http://dx.doi.org/10.1016/j.foreco.2005.11.015
    » http://dx.doi.org/10.1016/j.foreco.2005.11.015
  • Stape JL, Binkley D, Ryan MG, Fonseca S, Loos RA, Takahashi EN et al. The Brazil eucalyptus potential productivity project: Influence of water, nutrients and stand uniformity on wood production. Forest Ecology and Management 2010; 259(9): 1684-1694. http://dx.doi.org/10.1016/j.foreco.2010.01.012
    » http://dx.doi.org/10.1016/j.foreco.2010.01.012
  • Statsoft. Statistica – data analysis software system Version 10.0 [online]. Palo Alto: Statsoft Inc; 2010 [cited 2017 Apr 25]. Available from: http://www.statsoft.com
    » http://www.statsoft.com

Publication Dates

  • Publication in this collection
    01 July 2019
  • Date of issue
    2019

History

  • Received
    27 Oct 2017
  • Accepted
    12 July 2018
Instituto de Florestas da Universidade Federal Rural do Rio de Janeiro Rodovia BR 465 Km 7, CEP 23897-000, Tel.: (21) 2682 0558 | (21) 3787-4033 - Seropédica - RJ - Brazil
E-mail: floram@ufrrj.br