Remoção de nutrientes pela colheita de biomassa de Eucalyptus urograndis em diferentes tipos solos: macronutrients

ambiagua Revista Ambiente & Água Rev. Ambient. Água 1980-993X Instituto de Pesquisas Ambientais em Bacias Hidrográficas Resumo O manejo intensivo de povoamentos florestais pode aumentar a produção de biomassa, bem como aumentar a remoção de nutrientes do local. Assim, o objetivo deste estudo foi simular diferentes intensidades de colheita e calcular a eficiência de utilização de nutrientes para Eucalyptus urograndis em diferentes tipos de solo. O estudo foi realizado em um plantio do híbrido E. urograndis de sete anos de idade, no município de Telêmaco Borba, Paraná, Brasil. O local de estudo incluiu duas subáreas com solo arenoso e solo argiloso (Cambissolo Háplico e Latossolo Vermelho, respectivamente). Com os dados de estoque de biomassa e nutrientes foi simulada a remoção de nutrientes em cinco cenários de colheita diferentes. A eficiência de utilização de nutrientes foi obtida a partir da relação entre a quantidade de biomassa e os nutrientes de cada componente da árvore. A colheita da árvore inteira resultou na retirada de aproximadamente 61% dos nutrientes do local em solo arenoso, enquanto em solo argiloso, 57% dos nutrientes foram removidos. Com a colheita apenas da madeira do fuste comercial, apenas 22% dos nutrientes foram retirados do solo arenoso e 21% do solo argiloso. A madeira do caule foi o componente que apresentou os maiores valores da eficiência de utilização de nutrientes para todos os nutrientes analisados. Para alcançar a sustentabilidade nutricional dos povoamentos de E. urograndis, o melhor sistema de colheita envolve a remoção apenas da madeira do tronco comercial. Para a produção de madeira, os solos arenosos apresentam maior eficiência biológica de cálcio e magnésio quando comparados aos solos argilosos. 1. INTRODUCTION Use of renewable energy was boosted by the oil crisis in the 1970s, increasing the demand for energy biomass, which led to the search for available waste after forest harvesting operations (Egnell, 2017). However, intensively managed forest stands can increase biomass production, but may also increase the removal of nutrients from the site (Viera et al., 2011). In the case of energetic biomass, the adopted management can involve harvesting all parts of the tree, including the roots, and, in extreme cases, the removal of accumulated litter and understory (Viera et al., 2015). Viera et al. (2013) emphasize that intensively managed plantations with the use of short rotations, without predicting a minimum period required for nutrient replacement, have been found to be the main reason for the chemical exhaustion of soils. However, for long-term forest management, knowledge of the relationships between the amount of nutrients removed by harvests and the bioavailability of nutrients at the site is imperative in order to achieve sustainable energy standards under various rotations (Santana et al., 2008). New plantations have been implanted into different types of soils and climatic and management conditions (Kumaraswamy et al., 2015). There are differences between plantations and natural ecosystems, which are self-sustainable. Natural ecosystems are able to maintain their productivity and nutrient dynamics. However, forest plantations are a diversified form of land use, and in order to develop plantations that do not have significant negative impacts on the ecosystem, plantations need to be in harmony with the ecological properties of their environment. In this context, knowledge about the growth and nutrient allocation to the different biomass components of trees is essential for crop optimization and the nutritional management of plantations (Schumacher et al., 2011; Yan et al., 2017), because the removal of nutrients in the harvest depends on the productivity of the stands and the concentrations of nutrients in the components of biomass (Silva et al., 2013). The effects of harvesting for forest management are mainly related to the collection system adopted; but other factors, such as climate, soil type, and topography, may also be important (Olsson et al., 2017). The research hypothesis is that the removal of nutrients will occur more intensely in the site with sandy soil and this will result in a reduced number of rotations when compared to clayey soil. With the development of new hybrids and genotypes, it is essential to assess their performance under different edaphoclimatic conditions. The region of the present study does not have a historical silvicultural tradition. Considering this aspect, few studies evaluating the removal of nutrients in different scenarios of harvest and soil texture have been published. The objective of this study was to estimate the removal of macronutrients simulating different harvest intensities and to calculate the nutrient-use efficiency of nutrients in hybrid eucalyptus stands on different soil types. 2. MATERIAL AND METHODS 2.1. Characterization of the experimental area The study was carried out in stands of seven-year-old hybrid Eucalyptus urograndis (Eucalyptus grandis Hill ex Maiden x Eucalyptus urophylla S.T. Blake) at Farm Monte Alegre in the municipality of Telêmaco Borba, Paraná, Brazil. The climate of the region is characterized as Cfa, according to the climatic classification of Köppen. The average annual temperature is approximately 18.6°C and the average annual precipitation reaches 1,443 mm (Alvares et al., 2013). Within this study site, subareas with two distinct soil types were chosen. The first subarea had soil with a sandy texture (Cambisols Inceptisols), and the second had soil with a clayey texture (Ferralsols Oxisols). The determination of the physical and chemical attributes of the soil were carried out with the collection of 3 samples in each plot, totaling 12 samples per soil type, in the depths of 0 - 20 cm, 20 - 40 cm and 40 - 60 cm (Table 1). In each soil type, four sample plots of 2,550 m2 were plotted, with tree spacing of 3.0 m x 2.5 m, composed of 340 plants. The planting was done manually, with a spacing of 3.0 m x 2.5 m and an initial density of 1,333 plants per hectare. For the planting, a soil subsoiling was carried out in the planting line, with a depth of 45 cm, where a dosage of 200 kg ha-1 of natural rock phosphate was incorporated. After the planting, two other fertilizations were carried out; the first was basic fertilization with 15 kg ha-1 of nitrogen, 35 kg ha-1 of phosphor, 15 kg ha-1 of potassium, and the second was cover fertilization with 40 kg ha-1 of nitrogen, 5 kg ha-1 of phosphor, 65 kg ha-1 of potassium + 1,5 kg ha-1 of boron. 2.2. Biomass and nutrients Based on the diameter at breast-height (DBH) survey of all the trees present in the plots, 12 trees were selected in each plot for the determination of the above-ground biomass (leaves, branches, stembark and stemwood) through the destructive method, thus totaling in each type of soil. The biomass of the roots was analyzed from the selection of four medium DBH trees in each soil type. The root system (stump, fine roots < 2 mm, medium roots of 2.1 to 10 mm, and large roots > 10 mm) was extracted using a backhoe and manual excavation (shovels and spades), in an area of 7.5 m² surrounding the tree (according to spacing) to the depth of 1 m. In order to evaluate the understory (all naturally occurring vegetation, native or exotic) were collected all the biomass contained within the useful area (7 m²) of each sectioned tree for the sampling of the tree biomass, thus totaling 12 samples in each soil type, from which the average diameter tree minus a standard deviation, the average diameter tree and the average diameter tree plus one were selected. standard deviation. The evaluation of litter on the soil was carried out through the collection of 10 samples in an area of 30 cm x 30 cm distributed randomly within each plot, totaling 40 samples in each soil type. The selected trees were sectioned at ground level and divided into the following components: leaves, branches, stembark, stemwood and roots. The amount of green biomass in these components of the tree was determined by means of weighing with a hook scale in the field. Table 1. Chemical and physical attributes of distinct soils planted with E. urograndis in the region of Telêmaco Borba, Paraná, Brazil. Attribute Sandy Clayey 0-20 cm 20-40 cm 40-60 cm 0-20 cm 20-40 cm 40-60 cm Organic matter (%) 1.79(0.14) 1.32(0.22) 1.42(0.43) 3.39(0.12) 2.45(0.15) 1.72(0.05) pH (H2O) 3.97(0.07) 3.97(0.10) 3.95(0.02) 3.98(0.04) 4.19(0.04) 4.41(0.01) P*(mg dm-3) 1.61(0.18) 1.12(0.13) 0.89(0.31) 0.86(0.20) 0.68(0.20) 0.68(0.04) K*(mg dm-3) 30.92(1.48) 20.08(1.73) 35.11(7.71) 45.04(1.16) 32.59(1.28) 27.63(0.26) S (mg dm-3) 9.08(1.58) 10.28(2.04) 13.03(2.02) 26.13(2.33) 23.36(2.41) 11.26(3.90) B (mg dm-3) 0.55(0.03) 0.63(0.06) 0.67(0.05) 0.61(0.07) 0.61(0.08) 0.56(0.08) Cu (mg dm-3) 1.25(0.27) 1.22(0.29) 1.19(0.23) 2.65(0.14) 2.06(0.20) 1.40(0.09) Zn (mg dm-3) 0.54(0.07) 0.48(0.11) 0.43(0.09) 0.58(0.12) 0.25(0.01) 0.19(0.02) Ca (cmolc dm-3) 0.09(0.02) 0.06(0.01) 0.06(0.01) 0.18(0.05) 0.05(0.01) 0.03(0.01) Mg (cmolc dm-3) 0.07(0.01) 0.05(0.01) 0.05(0.01) 0.45(0.06) 0.11(0.02) 0.03(0.01) Effective cation exchange capacity(cmolc dm-3) 2.30(0.33) 2.41(0.16) 2.47(0.15) 3.72(0.17) 2.85(0.18) 2.29(0.07) Saturation per exchangeable base (%) 2.54(0.59) 1.53(0.40) 2.68(0.55) 2.89(0.40) 1.00(0.08) 0.75(0.10) Exchangeable aluminum saturation (%) 89.92(2.79) 96.19(0.87) 93.36(1.52) 80.02(1.99) 91.76(0.37) 94.25(0.77) Coarse sand: 2.0 - 0.2 mm (%) 39.89(6.41) 40.33(4.07) 40.72(3.88) 14.68(2.77) 16.51(0.59) 16.18(0.97) Fine sand: 2.0 - 0.05 mm (%) 40.36(7.24) 40.50(3.83) 37.05(4.96) 6.87(2.24) 5.79(1.80) 6.73(0.66) Silt: 0.05 - 0.002 mm (%) 4.19(0.27) 2.11(0.63) 3.67(0.81) 28.23(1.94) 31.48(2.78) 26.87(1.09) Clay: < 0.002 mm (%) 15.56(1.29) 17.06(1.26) 18.56(1.41) 50.22(1.50) 46.22(1.26) 50.22(0.96) * Determination of nutrients using the Mehlich-1 extractor (HCI + H2SO4). P: Phosphorus determined by atomic absorption spectrophotometry (660 nm); K: Potassium determined by flame photometry; S: Sulfur extracted with Ca (H2PO4)2 and determined by turbidimetry (440 nm); B: Boron extracted with curcumin and determined by spectrophotometry (460 nm); Cu: Copper and Zn: Zinc extracted with HCL (0.1 mol L-1) and determined by spectrophotometry atomic absorption (324.75 and 213.86 nm, respectively); Ca: Calcium and Mg: Magnesium, extracted with KCL (1.0 mol L-1) and determined by spectrophotometry atomic absorption (422.70 and 285.20 nm, respectively); ECEC= Effective Cation Exchange Capacity. Values in parentheses represent the standard errors. Representative samples were collected of each component of the biomass of trees, with the exception of stembark and stemwood, which were obtained by 4-cm-thick wooden discs (bark and wood) in the following positions depending on commercial height: base, 25%, 50%, 75%, 100% and tree tops, being considered the fraction of 100% (minimum diameter of 8 cm). After being weighed in the field, the samples were sent to the laboratory to be dried in an air-circulating oven at 70ºC for 72 hours. Dried samples were subsequently milled in a Wiley-type mill, followed by chemical analysis to determine the macronutrients according to Tedesco et al. (1995) and Miyazawa et al. (2009). Sulfuric digestion (H2SO4 + H2O2) was used for the determination of nitrogen through the Kjeldahl method, nitric-perchloric digestion (HNO3 + HCIO4 - [3:1]) for determination Ca (422.67 nm), Mg (285.21nm), with atomic absorption spectrophotometry, K with flame photometry, P (660.00 nm) with spectrophotometry and S (420.00 nm) by turbidimetry. The total tree biomass per hectare was determined based on the dry weight of each sample multiplied by the total wet biomass of each component of the tree, which was extrapolated based on the number of trees in one hectare. In the case of the understory and litter, the estimate was made by extrapolation of the sample area. For the estimation of the total amount of macronutrient, we computed the product between the nutrient content and the dry biomass for each component. A detailed description of the methodology adopted for biomass determination can also be found in Salvador et al. (2015). 2.3. Simulation of harvest intensities Estimation of nutrient removal was calculated considering the accumulated amount in the biomass and the outputs under five simulations of harvest intensities, with the commercial stemwood being considered with a minimum diameter of 8 cm, as follows: 1st simulation - Whole tree = harvest of whole biomass tree (stemwood with stembark, tree tops, leaves, branches and roots, leaving understory and litter); 2nd simulation - Full tree = harvest aboveground biomass (stemwood with stembark, tree tops, leaves and branches, leaving roots, understory and litter); 3rd simulation - Harvest of the stemwood with stembark and tree tops (leaving leaves, branches, roots, understory and litter); 4th simulation - Harvest of the commercial stemwood with stembark (leaving tree tops, leaves, branches, roots, understory and litter); 5th simulation - Harvest of the commercial stemwood (leaving stembark, leaves, branches, tree tops, roots, understory and litter). 2.4. Potential number of rotations and rate of nutrient removal The calculation of the potential number of rotation (PNR) for seven-year-old E. urograndis was performed using the ratio of the amount of nutrients remaining at the site (kg) to the amount of nutrients removed from the site (kg) according to the different harvest intensities for both types of soil (Equation 1). P N R = A m o u n t o f n u t r i e n t r e m a i n i n g a f t e r h a r v e s t i n k g h a - 1 A m o u n t o f n u t r i e n t h a r v e s t e d i n k g h a - 1 (1) The nutrient removal rate (NRR) was obtained by the ratio between the amount of nutrients (kg) removed with the harvest and the amount of biomass harvested (Mg) for each harvest intensity in the two types of soil (Equation 2). N R R = A m o u n t o f n u t r i e n t h a r v e s t e d i n k g h a - 1 B i o m a s s h a r v e s t e d i n M g h a - 1 (2) 2.5. Nutrient-use efficiency The nutrient-use efficiency was obtained according to a calculation proposed by Barros et al. (1986) using the relation between the amount of biomass and the amount of nutrients of each component, both with the same unit (Equation 3). N R R = A m o u n t o f b i o m a s s i n M g h a - 1 A m o u n t o f n u t r i e n t i n M g h a - 1 (3) 2.6. Statistical analysis The data were plotted in box plots for the identification and exclusion of outliers. After, Shapiro-Wilk (to prove the normality of the data) were applied as assumptions to the use of analysis of variance. Results that did not meet the assumptions were transformed by applying the natural logarithm or Box Cox. The amount of nutrients and nutrient-use efficiency were subjected to analysis of variance (ANOVA) for comparison of means, between the same components, in different soil types (P < 0.05). 3. RESULTS AND DISCUSSION 3.1. Harvest simulation based on different harvesting intensities Quantification of the nutrients in different soil types is essential to be able to establish sustainable harvest systems. Sandy soil had greater amounts of nutrients in the treetop wood when compared to clayey soil, with the exception of calcium and sulfur, which were not significant (Table 2). However, the clayey soil had a higher biomass and amount of nutrients than the sandy soil, except for sulfur. Harvesting tree tops in sandy soils can significantly increase nutrient exports when compared to harvesting this component in clayey soils. Table 2. Amounts of nutrients available in biomass and soil in stands of seven-year-old E. urograndis in sandy and clayey soils in the region of Telêmaco Borba, Paraná, Brazil. Solo Component Biomass N P K Ca Mg S Mg ha-1 kg ha-1 Sandy Leaves 2.61 ns* 52.73 ns 2.90 ns 28.15 ns 11.06 b 6.20 ns 2.40 ns Branches 6.95 ns 19.81 b 1.81 ns 33.65 ns 18.57 b 5.16 ns 2.64 ns Tree tops bark 1.42 ns 5.99 ns 1.04 ns 8.92 ns 12.06 a 4.12 ns 0.33 ns Tree tops wood 12.42 ns 27.56 a 2.44 a 42.28 a 6.07 ns 4.89 a 3.84 ns Stembark 12.07 ns 46.23 ns 6.44 ns 115.19 ns 100.67 ns 28.31 ns 4.84 ns Stemwood 180.26 ns 186.18 ns 12.02 ns 169.22 ns 87.37 ns 22.44 ns 40.45 ns Roots 43.05 ns 166.39 ns 7.78 ns 63.57 b 57.30 b 24.11 b 17.98 a Understore 1.08 ns 8.17 ns 0.72 ns 10.19 ns 2.61 ns 1.55 ns 0.92 ns Litter 17.35 b 142.29 b 5.59 b 23.09 b 96.52 b 28.77 b 11.89 a Soil -- 222.96 ns 7.24 ns 172.22 ns 84.17 ns 41.34 ns 64.78 ns Clayey Leaves 4.01 ns 88.98 ns 4.82 ns 46.44 ns 25.07 a 11.08 ns 4.18 ns Branches 11.44 ns 41.35 a 4.05 ns 42.30 ns 98.25 a 12.04 ns 2.87 ns Tree tops bark 1.28 ns 5.55 ns 1.21 ns 7.35 ns 9.28 b 5.91 ns 0.38 ns Tree tops wood 12.52 ns 16.86 b 1.44 b 19.07 b 8.25 ns 2.99 b 3.06 ns Stembark 12.58 ns 46.96 ns 5.98 ns 90.60 ns 157.94 ns 37.61 ns 3.49 ns Stemwood 211.21 ns 213.66 ns 13.30 ns 218.63 ns 132.86 ns 39.27 ns 52.50 ns Roots 36.95 ns 137.57 ns 8.64 ns 149.60 a 87.59 a 40.32 a 9.63 b Understore 0.65 ns 10.68 ns 0.60 ns 6.37 ns 2.42 ns 1.50 ns 1.48 ns Litter 20.93 a 169.12 a 6.42 a 29.31 a 165.11 a 41.61 a 9.70 b Soil -- 377.89 ns 4.44 ns 210.52 ns 104.21 ns 143.44 ns 121.50 ns * Mean of each component in the different soil types followed by ns, did not differ significantly at the level of 5% error. Source: Salvador et al. (2019). Each element should be analyzed based on the nutrients of the amount in the different components and as amount removed at harvest for reposition. Such importance is evident when observing that the leaves, in spite of presenting low biomass production, have high amounts of nutrients, mainly nitrogen, and the same analysis is valid for calcium in the stembark. The harvesting only of stem (third simulation), a system adopted by large companies to exploit the stembark and tree tops (wood and bark) as an energy source and stemwood for cellulose production, resulted in the removal of more than 43-50% of the total phosphor, calcium and potassium in sandy soil. However, in clayey soil, this harvesting system resulted in the removal of 39-43% of phosphorus, potassium and calcium (Figure 1). Not harvesting the stembark can mitigate the effects of biomass harvest on soil fertility. Corroborating this result, Achat et al. (2015), after studying different biomass harvesting systems, found that harvesting only the stemwood was able to reduce the loss of calcium by 56%. Santos et al. (2020), evaluating the estimation the of nutrient export in different eucalypts genotypes in southern Brazil, observed that the harvesting only the stemwood maintained Ca in the other biomass components, which varied from 71 to 82% (E. dunnii and E. benthamii (P1)). In contrast, considering the harvest of the wood with the stem bark, the permanence of the same nutrient in the area reduces to 26 to 37% (E. benthamii (P2) and E. uro globulus). Therefore, when removing the trunk, debarking of the trees in the stands should be carried out to reduce the export of Ca. In this way, soils with lower concentrations of this nutrient could have their productivity preserved for a longer period (Santos et al., 2019a). Figure 1. Removal of nutrients at different harvest intensities of seven-year-old E. urograndis in sandy and clayey soils in the region of Telêmaco Borba, Paraná, Brazil. Whole tree: harvest of whole biomass tree; Full tree: harvest above ground biomass; Stem: harvest of stemwood with stembark; SW+SB: commercial stemwood with stembark; SW: commercial stemwood. *Mean of each harvest’s intensities in the different soil types, followed by ns, did not differ significantly at the level of 5% error. Considering the fourth harvesting system (SW+SB without tree tops, leaves, branches, roots, understory and litter) in sandy site, nutrient removal reached 43% for potassium, 39% for calcium and 38% for phosphor, resulting in a decrease in nutrient loss of 26, 23 and 34%, respectively, compared to the first harvesting system. For the clayey soil, this harvesting system resulted in a higher phosphor, potassium and calcium removal (38%, 38% and 37%, respectively), and a lower removal of magnesium (23%) (Figure 1). In a study on nutrient removal rates under different harvesting intensities of the biomass of E. urophylla x E. globulus, Viera et al. (2015) observed that harvesting the stembark resulted in a removal of more than 47% of the total aboveground accumulated biomass for all macronutrients, with less than 35% of the stemwood being harvested. As such, it was observed that a biomass harvest of the whole tree resulted in greater removal of nutrients from the site. In the fifth harvesting system (SW), considering only the removal of commercial stemwood, nutrient exports in the sandy soil would not exceed 27% for s sulfur and lower than that for the other macronutrients. However, for clayey soil, the maximum nutrient removal was found to be 27% for potassium and the minimum was found to be 12% for magnesium (Figure 1). Under the first harvesting system, 48-72% and 36-77% of the nutrients were removed for sandy and clayey soils, respectively. Under the second harvesting system, this was 36-60% for the sandy soil and 32-61% for the clayey soil. Under the fifth harvesting system, nutrient removal was limited to 27% in sandy soils (sulfur) and 27% (potassium) in clayey soils, on average 39% and 36% less than under the first harvesting system, respectively. As such, we recommend using the third harvest system, which includes the only removal of the stemwood, as this is the most sustainable for the site, and will result in the best nutritional conditions for subsequent rotations. A study by Witschoreck (2008) reported that some of the lost nutrients can be replenished through the application of fertilizers, but the nutrients that can be replenished are usually limited to the triad of elements nitrogen, phosphorus and potassium. Given the distribution of biomass and nutrients in each component of the trees, it is possible to evaluate these during the development of the stand and crop, as this information can be useful for the future maintenance of the site. In general, for all harvesting systems, the sandy soil presented a higher percentage of removal of nitrogen, potassium and magnesium and sulfur (on average 35%, 49%, 36% and 35%, respectively). The clayey soil has greater phosphorus removal, an average of 49% among all harvesting systems. This result confirms our hypothesis that sandy soils would present the greatest removal of nutrients by the eucalyptus harvest. Sandy soils naturally present a low capacity for cation exchange, thus presenting the largest stock of nutrients in the eucalyptus biomass in this system. With the harvest, the nutritional impact on this soil is greater than on clay soils. 3.2. Potential number of rotations and rate of nutrient removal Considering an average of all the macronutrients for commercial stemwood harvesting (SW) results in an increase in the number of rotations 5.8 times greater than the whole tree harvest for the sandy soil and 4.8 times for the clayey soil (Table 3). Comparing the number of rotations between SW + SB and SE, the calcium increases by 2.9 times the number of rotations in the SW for both types of soil (Table 3). The first harvesting system (whole tree) allows one more rotation only for sulfur, while the clay soil, in this harvesting system, would support another rotation for nitrogen, magnesium and sulfur. Considering only the stemwood harvest (SW), phosphor, potassium and sulfur are the ones that most limit the growth of the next rotations (Table 3). Table 3. Estimation of the potential number of rotations based on the amount of nutrients in the biomass and available in the soil, and the nutrient removal rates according to the harvest intensity of E. urograndis biomass, at seven years old in sandy and clayey soils in the region of Telêmaco Borba, Paraná, Brazil. Harvest intensities Sandy Clayey N P K Ca Mg S N P K Ca Mg S Whole tree 0.74b* 0.39ns 0.45ns 0.63ns 0.75b 1.07b 1.01a 0.29ns 0.43ns 0.52ns 1.25a 1.74a Full tree 1.59ns 0.80ns 0.68b 0.90ns 1.30b 1.75ns 1.68ns 0.65ns 0.93a 0.83ns 2.08a 2.14ns Stem 2.30b 1.19ns 0.99b 1.31ns 1.79b 2.03ns 2.92a 1.32ns 1.44a 1.57ns 2.91a 2.51ns SW+SB 2.78b 1.60ns 1.34b 1.53ns 2.29b 2.31ns 3.25a 1.64ns 1.65a 1.72ns 3.37a 2.73ns SW 3.72b 2.99ns 2.94ns 4.45ns 6.44ns 2.71ns 4.19a 2.83ns 2.75ns 4.95ns 7.55ns 2.98ns Whole tree: harvest of whole biomass tree. Full tree: harvest above-ground biomass; Stem: harvest of stemwood with stembark; SW+SB: commercial stemwood with stembark; SW: commercial stemwood. *Mean of each harvest intensities, in the different soil types, followed by ns did not differ significantly at the level of 5% error. The results of the present study indicate that harvest intensity planning should be adopted locally, and should be based on the nutrient removal of each tree section harvested and the type of soil, as the clayey soil was suitable for a higher number of rotations of seven-year-old E. urograndis, in commercial harvesting systems of the stemwood. Viera et al. (2015) verified that in E. urophylla x E. globulus stands, phosphorus and calcium were the main nutrients to become limiting in subsequent rotations under a harvest system in which the stemwood and stembark were harvested. Considering the harvest of stemwood only, the limiting nutrients for both soil types were phosphorus, potassium and sulfur, with numbers less than three rotations in this harvesting system. If you consider harvesting stemwood + stembark, the limiting nutrients are phosphorus, potassium and calcium, with less than two rotations. Magnesium is the least limiting growth nutrient, with a potential of six rotations for sandy soil and seven rotations for clay soil. In the harvest of stemwood + stembark, a very common system by companies in the cellulose sector, the sandy soil has the lowest number of rotations for all the analyzed nutrients. This shows the fragility of this soil and the attention given to fertilizing the plantations implanted in this type of soil to guarantee the sustainability of the system. Just as the estimated number of rotations varies according to the biomass harvest intensity, the rate of nutrient removal also depends on which components will be harvested (Table 4). The SW harvest system has the lowest nutrient removal rate, as it has the lowest amount of nutrients per biomass. The lowest rates of nutrient removal were observed with the harvesting of commercial stemwood, independent on the soil type, due to the high amount of biomass and the low amount of nutrients of this component when compared to the others. When comparing the two types of soil, we found that the removal rate is analogous to the harvesting of the commercial stemwood without the stembark, since the harvest of the stembark with tree tops resulted in higher rates of nutrient removal in the clayey soil, which was possibly related to higher biomass production of this soil. Studies on the removal of nutrients through biomass harvesting provide a basis for the understanding of nutritional dynamics in stands, and are vital for sustainable production, as such studies enable us to predict situations that are critical to both productivity and chemical characteristics of the soil (Viera et al., 2015). Nitrogen, potassium and calcium are the nutrients that present the greatest reduction in the rate of removal of nutrients according to the harvest intensity for both types of soil. There is a 50% reduction in the rate of calcium removal if the harvesting system leaves the bark in place and removes only the stemwood. Regardless of the harvesting system, the rate of sulfur removal in the clayey soil is practically unchanged, since the largest amount of this element is found in the soil followed by stemwood, a component that is removed in all harvesting systems. As in the sandy soil, the sulfur stock in the soil is lower than in the clayey soil, thus it presents a decrease in the removal rate according to the harvest intensity. Table 4. Nutrient removal rate based on nutrient and biomass removed, according to harvest intensity of E. urograndis biomass at seven years old in sandy and clayey soils in the region of Telêmaco Borba, Paraná, Brazil. Harvest intensities Sandy Clayey N P K Ca Mg S N P K Ca Mg S Whole tree 1.96ns* 0.13ns 1.79b 1.14b 0.37b 0.28ns 1.83ns 0.13ns 1.91a 1.72a 0.50a 0.25ns Full tree 1.57ns 0.12ns 1.82ns 1.10b 0.33b 0.25ns 1.56ns 0.12ns 1.67ns 1.63a 0.41a 0.25ns Stem 1.29ns 0.11ns 1.63ns 1.00b 0.29b 0.24ns 1.14ns 0.09ns 1.35ns 1.24a 0.34a 0.24ns SW+SB 1.21ns 0.10ns 1.48ns 0.98b 0.26b 0.24ns 1.11ns 0.08ns 1.32ns 1.24a 0.33a 0.24ns SW 1.03ns 0.07ns 0.94ns 0.48ns 0.12b 0.22ns 1.01ns 0.06ns 1.04ns 0.63ns 0.19a 0.25ns Whole tree: harvest of whole biomass tree. Full tree: harvest above ground biomass; Stem: harvest of stemwood with stembark; SW+SB: commercial stemwood with stembark; SW: commercial stemwood. *Mean of each harvest intensities in the different soil types followed by ns did not differ significantly at the level of 5% error. 3.3. Nutrient-use efficiency The nutrient-use efficiency (NUE) corresponds to the rate of conversion of nutrients into biomass. This value demonstrates how many units of biomass are formed per unit of nutrient, and the higher the value, the more efficient the conversion of nutrients into biomass. By improving the NUE, stemwood productivity can be increased. The NUE is one of the fundamental parameters for the definition of the best management techniques, and consequently for the maintenance of the productive capacity of forest stands (Santana et al., 2002). Analyzing the NUE between the different components of the biomass, we found that the stemwood was the component that had the highest values of NUE for all nutrients analyzed, both in sandy and clayey soil, and the inverse was observed in the leaves for nitrogen, phosphorus, potassium, and sulfur. The stembark presented the lowest NUE for calcium, independent of the soil type (Table 5). The high NUE of wood indicates that the xylem used the lowest quantity of nutrients for each unit of biomass produced (Medeiros et al., 2020). NUE decreased in the following order: phosphorus > magnesium > sulfur > calcium > nitrogen > potassium for sandy and clayey soil. This result is similar to the order found by Santos et al. (2017), for E. urograndis in Rio de Janeiro, Brazil, where the authors observed NUE magnitudes in the order of phosphorus > magnesium> calcium > nitrogen > potassium. However, Schumacher et al. (2019), found inverse values for Eucalyptus spp. in Rio Grande do Sul, Brazil, and reported NUE based on total biomass in the following order: phosphorus > sulfur > magnesium> potassium > nitrogen > calcium. For the production of all components of the biomass, the elements phosphorus, magnesium and sulfur were the macronutrients most efficiently used in both types of soil analyzed (Table 5). However, considering only the stemwood component of the trees, since this is the component of main commercial interest in the stands, the NUE decreased in the following order: phosphor > magnesium > sulfur > calcium > nitrogen > potassium for sandy and clayey soil. Similar results, although with inversion in the distribution of some nutrients, were reported by Santos et al. (2019b) studying the nutritional efficiency of different eucalyptus genotypes in Southern Brazil (P > S > Mg > Ca > N > K). Table 5. Nutrient-use efficiency for a plantation of seven-year-old E. urograndis on sandy and clayey soils, in the region of Telêmaco Borba, Paraná, Brazil. Soil Component N P K Ca Mg S Sandy Leaves 48ns* 852ns 88ns 233a 407a 1069ns Branches 334a 298ns 195ns 377a 1045a 2575b Tree tops bark 219ns 1575ns 152ns 109b 321a 4136ns Tree tops wood 466b 5180b 297b 2157a 2828b 3116b Stembark 258ns 1831ns 101b 125a 439a 2462b Stemwood 979ns 14765ns 1089ns 221a 8354a 4323ns Roots 255ns 5385ns 671a 806ns 1778a 2374b Clayey Leaves 45ns 834ns 84ns 154b 348b 951ns Branches 227b 2952ns 238ns 12b 575b 3799a Tree tops bark 235ns 1321ns 180ns 137a 268b 3713ns Tree tops wood 731a 8885a 666a 1966b 4232a 4642a Stembark 254ns 2017ns 133a 83b 362b 3464a Stemwood 984ns 16888ns 1015ns 1640b 5851b 4217ns Roots 269ns 4275ns 247b 437ns 930b 3831a * Mean of each component, in the different soil types, followed by ns, did not differ significantly at the level of 5% error. The sandy soil has better NUE for calcium and magnesium for the SW harvest system when compared to the clay soil, possibly explained by a luxury consumption, since the clay soil had higher concentrations of these elements in the soil. Considering only the stemwood, the sandy soil showed a higher NUE for calcium and magnesium when compared to clayey soil. The use of fertilization in forest stands is based on the quantification of the nutrients exported from the site by the harvested biomass, and the relation between the biomass production and the nutrients contained in each component (Lafetá et al., 2018). In addition, Santana et al. (2002) point out that nutrient utilization efficiency may be related to the characteristics of each species and to water availability, and nutritional non-equilibrium in the soil-plant system may cause a limitation or excess of one or more available nutrients. Sandy soil had a higher biological efficiency of calcium and magnesium for the production of stemwood when compared to clayey soil. Magnesium was not a growth-limiting nutrient as observed in the simulation of the number of rotations. However, if the stemwood + stembark harvest is considered, calcium was a limiting nutrient to growth. In sandy soils, E. urograndis showed a higher nutrient-use efficiency than clayey soil for both stemwood and stembark, showing the luxury consumption of the species, since the clay soil had higher concentrations of this element in the soil. This study demonstrated the adaptation of eucalyptus to the different types of soils in Brazil, showing the main limiting nutrients for each type of soil in different harvesting systems. Therefore, the NUE is a parameter of great utility in the selection of species to be used in reforestation, especially in nutrient poor soils (Silva et al., 1983). 4. CONCLUSIONS The best harvesting system considers the removal of only commercial stemwood, as this will reduce the removal of nutrients and result in a better nutritional balance of the sites, independently of the soil type. Clay soils had a higher biological efficiency of calcium and magnesium for the production of stemwood when compared to clayey soil, possibly explained by a luxury consumption, since the clay soil had higher concentrations of these elements in the soil. 5. ACKNOWLEDGEMENTS We thank the Klabin Company, and the Monte Alegre farm, for their support during the development of this project, for allowing us to conduct our research in the study area, and for their logistical and financial support. The fourth and fifth authors thank the Coordination for the Improvement of Higher Education Personnel (CAPES) for the scholarship. 6. REFERENCES ACHAT, D. L.; DELEUZE, C.; LANDMANN, G.; POUSSE, N.; RANGER, J.; AUGUSTO, L. Quantifying consequences of removing harvesting residues on forest soils and tree growth - A meta-analysis. Forest Ecology and Management, v. 348, n. 15, p. 124-141, 2015. https://doi.org/10.1016/j.foreco.2015.03.042 ACHAT D. L. DELEUZE C. LANDMANN G. POUSSE N. RANGER J. AUGUSTO L. Quantifying consequences of removing harvesting residues on forest soils and tree growth - A meta-analysis Forest Ecology and Management, 348 15 124 141 2015 https://doi.org/10.1016/j.foreco.2015.03.042 ALVARES, C. A.; STAPE, J. L.; SENTELHAS, P. C.; GONÇALVES, J. L. M.; SPAROVEK, G. Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift, v. 22, n. 6, p. 1-18, 2013. https://doi.org/10.1127/0941-2948/2013/0507 ALVARES C. A. STAPE J. L. SENTELHAS P. C. GONÇALVES J. L. M. SPAROVEK G. Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift, 22 6 1 18 2013 https://doi.org/10.1127/0941-2948/2013/0507 BARROS, N. F.; NOVAIS, R. F.; CARMO, D. N.; NEVES, J. C. L. Classificação nutricional de sítios florestais - descrição de uma metodologia. Revista Árvore, v. 10, p. 112-120, 1986. BARROS N. F. NOVAIS R. F. CARMO D. N. NEVES J. C. L. Classificação nutricional de sítios florestais - descrição de uma metodologia Revista Árvore, 10 112 120 1986 EGNELL, G. A review of Nordic trials studying effects of biomass harvest intensity on subsequent forest production. Forest Ecology and Management, v. 383, n. 1, p. 27-36, 2017. https://doi.org/10.1016/j.foreco.2016.09.019 EGNELL G. A review of Nordic trials studying effects of biomass harvest intensity on subsequent forest production Forest Ecology and Management, 383 1 27 36 2017 https://doi.org/10.1016/j.foreco.2016.09.019 KUMARASWAMY, S.; MENDHAM, D. S.; GROVE, T. S.; O’CONNELL, A. M.; SANKARAN, K. V.; RANCE, S. J. Harvest residue effects on soil organic matter, nutrients and microbial biomass in eucalypt plantations in Kerala, India. Forest Ecology and Management, v. 328, n. 15, p. 140-149, 2015. https://doi.org/10.1016/j.foreco.2014.05.021 KUMARASWAMY S. MENDHAM D. S. GROVE T. S. O’CONNELL A. M. SANKARAN K. V. RANCE S. J. Harvest residue effects on soil organic matter, nutrients and microbial biomass in eucalypt plantations in Kerala, India Forest Ecology and Management, 328 15 140 149 2015 https://doi.org/10.1016/j.foreco.2014.05.021 LAFETÁ, B. O.; SANTANA, R. C.; NOGUEIRA, G. S.; NEVES, J. C. L.; PENIDO, T. M. A. Eficiência de utilização de macronutrientes em eucalipto por método não destrutivo estimados por redes neurais artificiais. Ciência Florestal, v. 28, n. 2, p. 613-623, 2018. https://doi.org/10.1016/j.foreco.2014.05.021 LAFETÁ B. O. SANTANA R. C. NOGUEIRA G. S. NEVES J. C. L. PENIDO T. M. A. Eficiência de utilização de macronutrientes em eucalipto por método não destrutivo estimados por redes neurais artificiais Ciência Florestal, 28 2 613 623 2018 https://doi.org/10.1016/j.foreco.2014.05.021 MEDEIROS, P. L.; SILVA, G. G. C.; OLIVEIRA, E. M. M.; RIBEIRO, C. O.; SILVA, J. M. S.; PIMENTA, A. S. Efficiency of nutrient use for biomass production of a Eucalyptus clone as a function of planting density in short-rotation cropping. Australian forestry, v. 83, n. 2, p. 66-74, 2020. https://doi.org/10.1080/00049158.2020.1774958 MEDEIROS P. L. SILVA G. G. C. OLIVEIRA E. M. M. RIBEIRO C. O. SILVA J. M. S. PIMENTA A. S. Efficiency of nutrient use for biomass production of a Eucalyptus clone as a function of planting density in short-rotation cropping. Australian forestry, 83 2 66 74 2020 https://doi.org/10.1080/00049158.2020.1774958 MIYAZAWA, M.; PAVAN, M. A.; MURAOKA, T.; CARMO, C. A. F. S.; MELO, W. J. Análises químicas de tecido vegetal. In: SILVA, F. C. da (ed.). Manual de análises químicas de solos, plantas e fertilizantes. Rio de Janeiro: Embrapa, 2009. p. 190-233. MIYAZAWA M. PAVAN M. A. MURAOKA T. CARMO C. A. F. S. MELO W. J. Análises químicas de tecido vegetal SILVA F. C. da Manual de análises químicas de solos, plantas e fertilizantes. Rio de Janeiro Embrapa, 2009 190 233 OLSSON, B. A.; ÂKERBLOM, S.; BISHOP, K.; EKLÖF, K.; RING, E. Does the harvest of logging residues and wood ash application affect the mobilization and bioavailability of trace metals? Forest Ecology and Management, v. 383, n. 1, p. 61-72, 2017. https://doi.org/10.1016/j.foreco.2016.09.017 OLSSON B. A. ÂKERBLOM S. BISHOP K. EKLÖF K. RING E. Does the harvest of logging residues and wood ash application affect the mobilization and bioavailability of trace metals? Forest Ecology and Management, 383 1 61 72 2017 https://doi.org/10.1016/j.foreco.2016.09.017 SALVADOR, S. M.; SCHUMACHER, M. V.; STAHL, J.; SANTOS, J. C. Biomassa arbórea e de subbosque e da serapilheira acumulada em povoamentos de Eucalyptus urophylla x Eucalyptus grandis. Ecologia e Nutrição Florestal, v. 3, n. 3, p. 82-93, 2015. https://doi.org/10.5902/2316980X19414 SALVADOR S. M. SCHUMACHER M. V. STAHL J. SANTOS J. C. Biomassa arbórea e de subbosque e da serapilheira acumulada em povoamentos de Eucalyptus urophylla x Eucalyptus grandis Ecologia e Nutrição Florestal, 3 3 82 93 2015 https://doi.org/10.5902/2316980X19414 SALVADOR, S. M.; SCHUMACHER, M. V.; STHAL, J.; LUDVICHAK, A. A.; MOMOLLI, D. R.; CONSENSA, C. B. Effects of soil type in nutrient amount in eucalyptus urograndis: Macronutrients. Journal of Experimental Agriculture International, v. 35, n. 3, p. 1-9, 2019. https://doi.org/10.9734/jeai/2019/v35i330204 SALVADOR S. M. SCHUMACHER M. V. STHAL J. LUDVICHAK A. A. MOMOLLI D. R. CONSENSA C. B. Effects of soil type in nutrient amount in eucalyptus urograndis: Macronutrients Journal of Experimental Agriculture International, 35 3 1 9 2019 https://doi.org/10.9734/jeai/2019/v35i330204 SANTANA, R. C.; BARROS, N. F.; NEVES, J. C. L. 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, v. 26, n. 4, p. 447-457, 2002. https://doi.org/10.1590/S0100-67622002000400007 SANTANA R. C. BARROS N. F. NEVES J. C. L. 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, 26 4 447 457 2002 https://doi.org/10.1590/S0100-67622002000400007 SANTANA, R. C.; BARROS, N. F.; NOVAIS, R. F.; LEITE, H. G.; COMERFORD, N. B. Alocação de nutrientes em plantios de eucalipto no Brasil. Revista Brasileira de Ciência do Solo, v. 32, n. spe, p. 2723-2733, 2008. https://doi.org/10.1590/S0100-06832008000700016 SANTANA R. C. BARROS N. F. NOVAIS R. F. LEITE H. G. COMERFORD N. B. Alocação de nutrientes em plantios de eucalipto no Brasil Revista Brasileira de Ciência do Solo, 32 spe 2723 2733 2008 https://doi.org/10.1590/S0100-06832008000700016 SANTOS, F. M.; CHAER, M.; DINIZ, A. R.; BALIEIRO, F. C. Nutrient cycling over five years of mixed-species plantations of Eucalyptus and Acacia on a sandy tropical soil. Forest Ecology and Management, v. 384, n. 1, p. 110-121, 2017. https://doi.org/10.1016/j.foreco.2016.10.041 SANTOS F. M. CHAER M. DINIZ A. R. BALIEIRO F. C. Nutrient cycling over five years of mixed-species plantations of Eucalyptus and Acacia on a sandy tropical soil Forest Ecology and Management, 384 1 110 121 2017 https://doi.org/10.1016/j.foreco.2016.10.041 SANTOS, K. F.; LUDVICHAK, A. A.; QUEIROZ, T. B.; SCHUMACHER, M. V.; ARAÚJO, E. F. Biomass production and nutrient content in different Eucalyptus genotypes in Pampa Gaúcho, Brazil. Revista Brasileira de Ciências Agrárias, v. 14, n. 4, e6575, 2019a. https://doi.org/10.5039/agraria.v14i4a6575 SANTOS K. F. LUDVICHAK A. A. QUEIROZ T. B. SCHUMACHER M. V. ARAÚJO E. F. Biomass production and nutrient content in different Eucalyptus genotypes in Pampa Gaúcho, Brazil Revista Brasileira de Ciências Agrárias, 14 4 e6575 2019a https://doi.org/10.5039/agraria.v14i4a6575 SANTOS, K. F.; SCHUMACHER, M. V.; LUDVICHAK, A. A.; ARAÚJO, E. F. Biomass and stock of nutrients in different genotypes of eucalypts in Southern Brazil. Journal of Experimental Agriculture International, v. 34, n. 4, p. 1-12, 2019b. https://doi.org/10.9734/jeai/2019/v34i430178 SANTOS K. F. SCHUMACHER M. V. LUDVICHAK A. A. ARAÚJO E. F. Biomass and stock of nutrients in different genotypes of eucalypts in Southern Brazil Journal of Experimental Agriculture International, 34 4 1 12 2019b https://doi.org/10.9734/jeai/2019/v34i430178 SANTOS, K. F.; LUDVICHAK, A. A.; QUEIROZ, T. B.; SCHUMACHER, M. V.; ARAÚJO, E. F.; MOMOLLI, D. R. Estimation of nutrient export in eucalyptus genotypes under different harvest intensities in southern Brazil. Idesia, v. 38, n. 1, p. 33-38, 2020. http://dx.doi.org/10.4067/S0718-34292020000100033 SANTOS K. F. LUDVICHAK A. A. QUEIROZ T. B. SCHUMACHER M. V. ARAÚJO E. F. MOMOLLI D. R. Estimation of nutrient export in eucalyptus genotypes under different harvest intensities in southern Brazil Idesia, 38 1 33 38 2020 http://dx.doi.org/10.4067/S0718-34292020000100033 SCHUMACHER, M. V.; WITSCHORECK, R.; CALIL, F. N.; LOPES, V. G.; VIERA, M. Produção de biomassa no corte raso em plantio de Araucaria angustifolia (Bertol.) Kuntze de 27 anos de idade em Quedas Do Iguaçu, PR. Ciência Florestal, v. 21, n. 1, p. 53-62, 2011. https://doi.org/10.5902/198050982747 SCHUMACHER M. V. WITSCHORECK R. CALIL F. N. LOPES V. G. VIERA M. Produção de biomassa no corte raso em plantio de Araucaria angustifolia (Bertol.) Kuntze de 27 anos de idade em Quedas Do Iguaçu, PR Ciência Florestal, 21 1 53 62 2011 https://doi.org/10.5902/198050982747 SCHUMACHER, M. V.; WITSCHORECK, R.; CALIL, F. N.; LOPES, V. G. Manejo da biomassa e sustentabilidade nutricional em povoamentos de Eucalyptus spp. em pequenas propriedades rurais. Ciência Florestal, v. 29, n. 1, p. 144-156, 2019. http://dx.doi.org/10.5902/198050985135 SCHUMACHER M. V. WITSCHORECK R. CALIL F. N. LOPES V. G. Manejo da biomassa e sustentabilidade nutricional em povoamentos de Eucalyptus spp. em pequenas propriedades rurais Ciência Florestal, 29 1 144 156 2019 http://dx.doi.org/10.5902/198050985135 SILVA, H. D.; POGGIANI, F.; COELHO, L. C. Eficiência de utilização de nutrientes em cinco espécies de Eucalyptus. Boletim de Pesquisa Florestal, n. 6/7, p.1-8, 1983. SILVA H. D. POGGIANI F. COELHO L. C. Eficiência de utilização de nutrientes em cinco espécies de Eucalyptus Boletim de Pesquisa Florestal, 6/7 1 8 1983 SILVA, P. H. M.; POGGIANI, F.; LIBARDI, P. L.; GONÇALVES, N. A. Fertilizer management of eucalypt plantations on sandy soil in Brazil: Initial growth and nutrient cycling. Forest Ecology and Management, v. 301, n. 1, p. 67-78, 2013. https://doi.org/10.1016/j.foreco.2012.10.033 SILVA P. H. M. POGGIANI F. LIBARDI P. L. GONÇALVES N. A. Fertilizer management of eucalypt plantations on sandy soil in Brazil: Initial growth and nutrient cycling Forest Ecology and Management, 301 1 67 78 2013 https://doi.org/10.1016/j.foreco.2012.10.033 TEDESCO, M. J.; GINELLO, C.; BISSANI, C. A.; BOHNEN, H.; VOLKWEISS, S. J. Análise de solo, plantas e outros materiais. Porto Alegre: UFRGS, 1995. 174 p. TEDESCO M. J. GINELLO C. BISSANI C. A. BOHNEN H. VOLKWEISS S. J. Análise de solo, plantas e outros materiais. Porto Alegre UFRGS, 1995 174 VIERA, M.; SCHUMACHER, M. V.; BONACINA, D. M. Biomassa e nutrientes removidos no primeiro desbaste de um povoamento de Pinus taeda L. em Cambará do Sul, RS. Revista Árvore, v. 35, n. 3, p. 371-379, 2011. https://doi.org/10.1590/S0100-67622011000300001 VIERA M. SCHUMACHER M. V. BONACINA D. M. Biomassa e nutrientes removidos no primeiro desbaste de um povoamento de Pinus taeda L. em Cambará do Sul, RS Revista Árvore, 35 3 371 379 2011 https://doi.org/10.1590/S0100-67622011000300001 VIERA, M.; SCHUMACHER, M. V.; KLEINPAUL, I. S. Estoque de nutrientes em consórcios de Eucalyptus urograndis, Acacia mearnsii e Zea mays. Ciência Florestal, v. 23, n. 3, p. 317-327, 2013. https://doi.org/10.5902/1980509810543 VIERA M. SCHUMACHER M. V. KLEINPAUL I. S. Estoque de nutrientes em consórcios de Eucalyptus urograndis, Acacia mearnsii e Zea mays Ciência Florestal, 23 3 317 327 2013 https://doi.org/10.5902/1980509810543 VIERA, M.; SCHUMACHER, M. V.; TRÜBY, P.; ARAÚJO, E. F. Implicações nutricionais com base em diferentes intensidades de colheita da biomassa de Eucalyptus urophylla x Eucalyptus globulus. Ciência Rural, v. 45, n. 3, p. 432-439, 2015. https://doi.org/10.1590/0103-8478cr20120367 VIERA M. SCHUMACHER M. V. TRÜBY P. ARAÚJO E. F. Implicações nutricionais com base em diferentes intensidades de colheita da biomassa de Eucalyptus urophylla x Eucalyptus globulus Ciência Rural, 45 3 432 439 2015 https://doi.org/10.1590/0103-8478cr20120367 WITSCHORECK, R. Biomassa e nutrientes no corte raso de um povoamento de Pinus taeda L. de 17 anos de idade no município de Cambará do Sul-RS. 2008. 80f. Dissertação (Mestrado em Engenharia Florestal) - Universidade Federal de Santa Maria, Santa Maria, 2008. WITSCHORECK R. Biomassa e nutrientes no corte raso de um povoamento de Pinus taeda L. de 17 anos de idade no município de Cambará do Sul-RS. 2008 80f Mestrado Universidade Federal de Santa Maria, Santa Maria YAN, T.; ZHU, J.; YANG, K.; YU, L.; ZHANG, J. Nutrient removal under different harvesting scenarios for larch plantations in northeast China: Implications for nutrient conservation and management. Forest Ecology and Management, v. 400, n. 15, p. 150-158, 2017. https://doi.org/10.1016/j.foreco.2017.06.004 YAN T. ZHU J. YANG K. YU L. ZHANG J. Nutrient removal under different harvesting scenarios for larch plantations in northeast China: Implications for nutrient conservation and management Forest Ecology and Management, 400 15 150 158 2017 https://doi.org/10.1016/j.foreco.2017.06.004