Fertilization Response, Light Use, and Growth Efficiency in Eucalyptus Plantations across Soil and Climate Gradients in Brazil

Fertilization increases productivity in Eucalyptus plantations, but losses in productivity associated with soil fertility continue at operational scales. In this study, we evaluated the fertilization response (FR), light use efficiency (LUE) and growth efficiency (GE), i.e., the amount of wood biomass accumulated per unit of light absorbed (LUE) and per unit of leaf area index of Eucalyptus plantations. We used a " twin plot " approach, with 161 blocks representing 52,700 ha of planted forests that spanned a broad range of edaphoclimatic conditions in southeastern Brazil. The normal plots (NP) were part of a permanent inventory network, whereas the twin plots (TP) received extra high levels of fertilization and extra weed control after fertilization. The intensive management (twin plots) led to a large increase of 5.3 Mg¨há1¨year´1 of wood increment. The region without dry periods and with soils with high clay content was most responsive to fertilization, with a 15% increment in the LUE and 10% increase in the GE of the TPs compared with those of the NPs. Our results suggested that water availability was the primary element affecting productivity and potential response to fertilization. With this information, decisions can be made on which regions should receive priority fertilization investments. However, more research is required to determine the most limiting nutrient in each type of environment.


Introduction
In recent decades, the global planted area of Eucalyptus has increased significantly, covering approximately 10% of the planted forest area worldwide (approximately 12 million ha) [1].In Brazil, Eucalyptus forests occupy 5.1 million ha, of which over 1.1 million are in São Paulo [2].The supply chain activities of the wood industry, including production, harvesting, and transportation, contribute to the social, economic, and environmental well-being of society.For example, approximately 0.97 ha of native forest is preserved per 1.0 ha of Eucalyptus plantation in Brazil [2]; thus, millions of hectares of native forests are preserved.Eucalyptus is the most productive planted forest in Brazil, with an average productivity of 40 m 3 ¨ha ´1¨year ´1 [3].This high productivity is linked to the availability of natural resources, genetic characteristics, and silvicultural management [4].Productivity in São Paulo State is Forests 2016, 7, 117 2 of 12 limited by nutrient deficiencies in soils of low fertility, high nutrient export at harvest, and inefficient nutrient use.The lack of nutrients such as nitrogen, phosphorus, boron, potassium, and magnesium or the excess of nutrients such as calcium, magnesium, and copper limits productivity [5].
Among the natural resources required for plant growth (i.e., water, light, and nutrients), nutrient availability is the resource most easily manipulated by foresters, throughout practices of soil preparation and conservation and applications of fertilizers.The efficiency of this management approach is site-specific for each type of nutrient applied [5][6][7][8].
Although knowledge on tree nutrition has advanced, losses in timber productivity associated with nutritional problems continue at operational scales [6].Trees without equilibrated nutrition may lead to a heterogeneous population, compromising the site productivity [9].For industrial forests, understanding the value and spatial distribution of the fertilization response is essential for optimal investments in silviculture.A twin plots (TP) design is one approach to quantify the potential response to fertilization.This method establishes paired plots (control and treated plots) across a number of selected sites that represent the landscape [8].This approach differs from that of the classic experimental design of fertilization experiments because of a greater power of statistical inference across many stands, leading to a broader understanding of fertilization response by covering most of the environmental variability in a short period of time [6].The control plot may be a normal plot of a permanent inventory network used to measure the actual productivity of a forest under the traditional fertilization regime of a particular owner.The treated plot receives an intensive treatment (high levels of fertilization and weed control) for quantitative insights into the factors that limit operational productivity.
The relationship between resource availability and productivity is expressed in the following equation: Production = supply + resource capture efficiency + resource use efficiency.Thus, productivity is strongly influenced by the resource supply and is dependent on the efficiency of the plant in capturing and using these resources [10,11].Resource use efficiency is an indication of how plants use resources (e.g., biomass produced per unit of resource consumed) and is a determining factor in wood growth.Growth efficiency, e.g., the amount of wood growth per unit of leaf area, also express how active are the leaves in converting carbon into biomass; this process is strongly influenced by leaf nutritional status [12].
The use efficiency of different natural resources has been examined in several studies (i.e., water use efficiency-WUE; light use efficiency-LUE; and nutrient use efficiency-NUE; [10,13].With increases in the resource supply, e.g., through fertilization, plants use not only nutrients but also other resources more efficiently [10].For example, an increase in the availability of water can increase LUE, WUE, and nitrogen use efficiency [11].Soil fertility influences LUE by increasing the amount of light absorbed with an increase in leaf area index (LAI) or by increasing plant efficiency in using radiation absorbed with an increase in leaf retention, as occurs in response to potassium fertilization [14].
Even though nutritional management of Eucalyptus leads to significant increases in the productivity and sustainability of forest plantations in Brazil [15][16][17], few studies have quantified light use and growth efficiency in response to fertilization.The objective of our study was to gain insight into fertilization response, light use efficiency, and growth efficiency of Eucalyptus plantations along an edaphoclimatic gradient in southeastern Brazil.Moreover, this study addressed the following questions: What is the difference between actual and potential productivity of Eucalyptus forests in southern Brazil?Is there variability in wood productivity on a temporal and spatial scale that can be adjusted by forestry management?Can the increases in productivity correlated with fertilization be explained by increases in light use and growth efficiencies?

Site Description and Location
A total of 161 plots (400 m 2 , containing 66 trees per plot) were installed in Eucalyptus plantations at Suzano Pulp and Paper Company in São Paulo State.The plots were selected based on a random sample of 1832 inventory plots (ranging between 2 and 4 years of age) that spanned a broad range of soil and climatic conditions separated into three regions (Table 1).The selected plots represented a 52,700-ha area of planted forests in three regions in the state of São Paulo (Figure 1).A total of 161 plots (400 m 2 , containing 66 trees per plot) were installed in Eucalyptus plantations at Suzano Pulp and Paper Company in São Paulo State.The plots were selected based on a random sample of 1832 inventory plots (ranging between 2 and 4 years of age) that spanned a broad range of soil and climatic conditions separated into three regions (Table 1).The selected plots represented a 52,700-ha area of planted forests in three regions in the state of São Paulo (Figure 1).
Each inventory plot (normal plot-NP) was paired with a second plot (twin plot-TP) of identical dimensions separated by approximately 20 m.A paired t-test was used to test for initial differences in wood biomass between normal and twin plots.
Following plot establishment, soils were sampled at 0-15 and 15-30 cm depths at four locations inside each plot.The soils were analyzed using the methodology described by Raij et al [18] for the following chemical variables: pH (CaCl2), organic matter (OM), P, Ca, K, Mg, and micronutrients.Nitrogen was not directly included in the original analysis, but we used organic matter amount as an indicator of nitrogen soil content [19][20][21].Brazilian Eucalyptus plantations typically do not respond to nitrogen fertilization [21].Source: Sentelhas et al [22] All stands received maintenance fertilization with N, P, K, Ca, and Mg according to the conventional prescriptions of the company (Table 2).Twin plots received an additional fertilization that contained 4 t•ha −1 of lime, 2.5 t•ha −1 of NPK 18:08:18, 800 kg•ha −1 of Single Superphosphate, and 300 kg•ha −1 of FTE (Fritted Traced Elements: 1.8% B, 0.8% Cu, 3% Fe, 2% Mn, 9% Zn, and 0.1% Mo); the fertilizer was applied four times from January 2004 to April 2005.To eliminate the effects of other factors on tree growth, pests, diseases, and weeds were effectively controlled in the plot areas.Both treatments were weed-controlled, but the twin plot received an extra weed control when necessary to avoid weed effects after the extra fertilization (three times during the study period).Each inventory plot (normal plot-NP) was paired with a second plot (twin plot-TP) of identical dimensions separated by approximately 20 m.A paired t-test was used to test for initial differences in wood biomass between normal and twin plots.
Following plot establishment, soils were sampled at 0-15 and 15-30 cm depths at four locations inside each plot.The soils were analyzed using the methodology described by Raij et al. [18] for the following chemical variables: pH (CaCl 2 ), organic matter (OM), P, Ca, K, Mg, and micronutrients.Nitrogen was not directly included in the original analysis, but we used organic matter amount as an indicator of nitrogen soil content [19][20][21].Brazilian Eucalyptus plantations typically do not respond to nitrogen fertilization [21].All stands received maintenance fertilization with N, P, K, Ca, and Mg according to the conventional prescriptions of the company (Table 2).Twin plots received an additional fertilization that contained 4 t¨ha ´1 of lime, 2.5 t¨ha ´1 of NPK 18:08:18, 800 kg¨ha ´1 of Single Superphosphate, and 300 kg¨ha ´1 of FTE (Fritted Traced Elements: 1.8% B, 0.8% Cu, 3% Fe, 2% Mn, 9% Zn, and 0.1% Mo); the fertilizer was applied four times from January 2004 to April 2005.To eliminate the effects of other factors on tree growth, pests, diseases, and weeds were effectively controlled in the plot areas.Both treatments were weed-controlled, but the twin plot received an extra weed control when necessary to avoid weed effects after the extra fertilization (three times during the study period).Two primary soil types were included in the study areas.In regions 1 and 2, the typical soil is an Entisol-Psamment, which has a sandy texture, low clay content, low cation exchange capacity, and low base saturation.In region 3, the typical soil is an Ultisol, which has a clay texture (clay content > 40% in the 0-30 cm layer), with twice the organic matter content and higher K content than in regions 1 and 2 (Table 3).The differences in physical and chemical properties between the two soil layers (0-15 and 15-30 cm) were small.The nutrient content of the surface layer (0-15 cm) was higher than that in the deeper layer.

Growth Determination and Fertilization Response
We measured all tree diameters (D) and 20% of tree heights: heights of the other trees were estimated with a site-specific hypsometric regression (H est = lnβ 0 ˆln(´β 1 ˆ(1/D)), where Hest is the estimated height and β 0 and β 1 are the coefficients estimated for each region.
Tree volume (Vol) was estimated from tree diameter (DBH, in cm).Tree height (H, in m) was estimated using a specific taper model developed through the Smalian method and then the Schumacher-Hall [24] model was used to determine the individual volume: Vol = e( ´9.56170571390879 + 1.94813158669387 ˆLn(D) + 0.868952145305574 ˆLn(H est )) (p = 274; R 2 = 0.98).Basic wood density was determined after cubage of trees using the hydrostatic scale method [25], which was also used to calculate stem biomass.
The mean annual increment at 7 years of age for each plot was calculated to compare the paired plots and the fertilization response (FR) of each pair, in each period.The FR was determined using the Forests 2016, 7, 117 5 of 12 equation proposed by Ferreira and Stape [6]: FR = (CAW T /BIOM T ´CAW N /BIOM N ) ˆ((BIOM T + BIOM N )/2), where CAW = current annual woody increment (Mg¨ha ´1¨year ´1); BIOM = stem initial biomass (Mg¨ha ´1); and T = twin plot and N = normal plot.

Leaf Area Index, Light Use Efficiency, and Growth Efficiency
Leaf area index (LAI) was determined (m 2 leaf.m ´2soil ) using an indirect method of hemispherical photography.Photographs of the forest canopy were collected using an Opteka Fisheye hemispheric lens 0.22ˆ(New York, USA) and a Fujifilm S5000 digital camera (Tokyo, Japan).A total of six photos were taken inside each plot (three photos within rows and three between rows).All photographs were taken using a tripod leveled and positioned centrally with the north axis (azimuth = 0 ˝).
The photographs were analyzed with Hemisfer ® software [26] to calculate the vegetation area index (VAI) of the canopy.To estimate LAI, we used the calibration equation (LAI = ln((1.24285ˆVAI) ´0.82229 ) developed by Giunti Neto et al. [27].
Growth efficiency (GE) was calculated [28] using the relationship between current annual increment (CAI) and LAI for a determined period (n): GE = CAI n /LAI n .
LUE was estimated using the equation proposed by Landsberg and Gower [29] (LUE = CAI n ˆ100/APAR n ).Absorption of photosynthetically active radiation was calculated based on the Lambert-Beer law that relates LAI and the coefficient of forest light extinction with the incident radiation (APAR n = PAR ´PAR ´k.LAI ), where k is the extinction coefficient of PAR for Eucalyptus, 0.45.We used this value based on some work performed for clonal Eucalyptus in Brazil [11,30].

Statistical Analyses
Each pair of twin plots was treated as a repetition and the twin and normal plots were the treatments.Tree growth, CAW, biomass, tree density, LAI, GE, and LUE after treatment with fertilization were analyzed using paired t-tests.For comparisons of average growth rate and LAI among different regions, the data were submitted to analysis of variance, hypothesis testing, and Tukey's tests, with a significance level of 0.05.The analysis of variance was based on a completely randomized design.

Productivity Plots
As expected, tree volume was not significantly different between the control and twin plots of each pair before the fertilization treatment (158 m 3 ¨ha ´1 (TP) vs. 157 m 3 ¨ha ´1 (NP), p = 0.241; Figure 2).With this equality, the correct interpretation of comparisons between plots after treatment application was assured.Between 2003 and 2006, the average increase in the annual increment of TP plots was 24% compared with NP plots (Figure 3).In the first year (2003)(2004), after the additional fertilization, a small increase was observed in the twin plots (6%).However, responses to fertilization were more evident in the second year of evaluation (2004)(2005), and the approximate increases were 30% in region 1, 44% in region 2, and 70% in region 3, with an average of 45%.In the final period (2005)(2006), the growth rate in the plantations decreased, and the relative difference between treatments was 30%.At 7 years of age, the average attainable productivity of the twin plots (TP) was 49 m 3 ¨ha ´1¨year ´1 (22 Mg¨ha ´1¨year ´1), which was 11% higher than the productivity of the normal plots (NP), i.e., 44 m 3 ¨ha ´1¨year ´1 (20 Mg¨ha ´1¨year ´1).
The mean annual increment at 7 years of age (MAI7) in NPs (54 m 3 ¨ha ´1¨year ´1) and TPs (62 m 3 ¨ha ´1¨year ´1) was the highest in region 3. Region 1 had the lowest MAI7, of 38 m 3 ¨ha ´1¨year ´1 (NP) and 41 m 3 ¨ha ´1¨year ´1 (TP).The greatest difference between the actual and the attainable productivity was in region 3 (15%).The maximum productivity occurred in region 3, i.e., 82 m 3 ¨ha ´1¨year ´1, followed by regions 2 and 1 with maximum productivity of 78 and 57 m 3 ¨ha ´1¨year ´1, respectively (Table 4). 1 n = Total number of plots.

Fertilization Response (FR)
The average fertilization response (FR) between 2003 and 2006 was 5.3 Mg¨ha ´1¨year ´1.Region 3 was the most responsive region, followed by regions 2 and 1 (Table 5).The FR was normally distributed (Figure 4), with ~85% of the data showing a positive FR and 50% of the plots with an FR equal to or higher than the average (5.3Mg¨ha ´1¨year ´1); similar proportions were found by Ferreira & Stape [6].

Fertilization Response (FR)
The average fertilization response (FR) between 2003 and 2006 was 5.3 Mg•ha −1 •year −1 .Region 3 was the most responsive region, followed by regions 2 and 1 (Table 5).The FR was normally distributed (Figure 4), with ~85% of the data showing a positive FR and 50% of the plots with an FR equal to or higher than the average (5.3Mg•ha −1 •year −1 ); similar proportions were found by Ferreira & Stape [6].Fertilization response was also affected by soil characteristics (Table 6).Considering all regions together, organic matter and clay content were the most important variables.Region 1 was strongly  Fertilization response was also affected by soil characteristics (Table 6).Considering all regions together, organic matter and clay content were the most important variables.Region 1 was strongly affected by physical characteristics, such as clay and sand content.Chemical variables such as potassium and magnesium were more strongly related to FR in region 2. In region 3, the sum of bases and percentage of silt were the most strongly related soil variables influencing the FR.

Leaf Area Index, Light Use Efficiency, and Growth Efficiency
In the absence of nutritional limitation, the LAI of the twin plots was 2.6 m 2 ¨m´2 , compared with the LAI of 2.2 m 2 ¨m´2 in the normal plots, which was an increase of 15% in the LAI of the TPs.However, in region 3, with the highest productivity, the difference in LAI between NPs and TPs was not significant (Figure 5a).affected by physical characteristics, such as clay and sand content.Chemical variables such as potassium and magnesium were more strongly related to FR in region 2. In region 3, the sum of bases and percentage of silt were the most strongly related soil variables influencing the FR.

Leaf Area Index, Light Use Efficiency, and Growth Efficiency
In the absence of nutritional limitation, the LAI of the twin plots was 2.6 m 2 •m −2 , compared with the LAI of 2.2 m 2 •m −2 in the normal plots, which was an increase of 15% in the LAI of the TPs.However, in region 3, with the highest productivity, the difference in LAI between NPs and TPs was not significant (Figure 5a).
In additional, the growth efficiency (GE) of the twin plots was higher (10%) than that of the normal plots, with values of 5.8 and 5.3 Mg•ha −1 •year −1 , respectively, the difference was significant (p > 0.05; Figure 5b).
The light use efficiency (LUE) in the twin plots, 0.65 g•MJ −1 , was 15% higher than that in the normal plots, 0.56 g•MJ −1 (p > 0.05; Figure 6).In additional, the growth efficiency (GE) of the twin plots was higher (10%) than that of the normal plots, with values of 5.8 and 5.3 Mg¨ha ´1¨year ´1, respectively, the difference was significant (p > 0.05; Figure 5b).
The light use efficiency (LUE) in the twin plots, 0.65 g¨MJ ´1, was 15% higher than that in the normal plots, 0.56 g¨MJ ´1 (p > 0.05; Figure 6).

Discussion
The additional fertilization led to large increases in wood increment, with an average response in the annual increment of 5.3 Mg 3 •ha −1 •year −1 .Similar results were found in other studies, with responses of 4.8 Mg 3 •ha −1 •year −1 [8] and 4 Mg 3 •ha −1 •year −1 [6] in Eucalyptus plantations in São Paulo State.
Region 3 was the most responsive to fertilization, which can be explained primarily by the climate, with no dry periods, and by the high clay content in the soils; consequently, the water retention capacity and organic matter content were high.Plant growth is controlled by the most limiting resource, according to Liebig's law of the minimum [32].In tropical regions where water and light are not limiting resources, nutrition starts to be important and normally limits wood growth [29].For region 3, the negative correlation of FR with the sum of bases and the silt content, which is an indicator of the primary source of nutrients [33], revealed nutritional aspects interfere more with FR than water availability.
In region 1, the climate has extended dry periods and the soils are sandy, which resulted in the lowest productivity (Table 4).Soil physical properties, particularly the clay content, are directly correlated with wood quality (lignin and holocellulose content) and the productive capacity of a site [34].In this region, where water is a limiting factor, the extra addition of nutrients did not lead to an extra amount of wood, confirmed by the positive correlation between the FR and sand content.
The largest response to fertilization was in the region that also had the highest productivity (Region 3), which was different from the response observed by Ferreira and Stape [6].Therefore, the area with the highest nutrient limitations on growth was also the most productive.Thus, water, the primary growth factor [17,11], was not limiting at this site, whereas in the other regions, particularly in region 1, the most limiting factor was water.
In the twin plots, LAI increased by 15% and GE increased by 10% compared with the normal plots.The increase in leaf area might explain the increase in productivity.Albaugh et al. [35] found much larger gains, an increase of 101% in the LAI, in response to fertilization in loblolly pine.In Eucalyptus nitens plantations, nitrogen fertilization increased the LAI by up to 3.1 units, i.e., 56% increment [36].
Fertilization increased the LUE and the GE, which was also reported by Binkley et al [37] for Eucalyptus and for other forest trees, including loblolly pine [38] and Liquidambar styraciflua [39].These responses might be a result of higher photosynthetic rates with the addition of nutrients [40] or caused by less carbon partitioning to shoots [10,41].

Discussion
The additional fertilization led to large increases in wood increment, with an average response the annual increment of 5.3 Mg 3 ¨ha ´1¨year ´1.Similar results were found in other studies, with responses of 4.8 Mg 3 ¨ha ´1¨year ´1 [8] and 4 Mg 3 ¨ha ´1¨year ´1 [6] in Eucalyptus plantations in São Paulo State.
Region 3 was the most responsive to fertilization, which can be explained primarily by the climate, with no dry periods, and by the high clay content in the soils; consequently, the water retention capacity and organic matter content were high.Plant growth is controlled by the most limiting resource, according to Liebig's law of the minimum [32].In tropical regions where water and light are not limiting resources, nutrition starts to be important and normally limits wood growth [29].For region 3, the negative correlation of FR with the sum of bases and the silt content, which is an indicator of the primary source of nutrients [33], revealed nutritional aspects interfere more with FR than water availability.
In region 1, the climate has extended dry periods and the soils are sandy, which resulted in the lowest productivity (Table 4).Soil physical properties, particularly the clay content, are directly correlated with wood quality (lignin and holocellulose content) and the productive capacity of a site [34].In this region, where water is a limiting factor, the extra addition of nutrients did not lead to an extra amount of wood, confirmed by the positive correlation between the FR and sand content.
The largest response to fertilization was in the region that also had the highest productivity (Region 3), which was different from the response observed by Ferreira and Stape [6].Therefore, the area with the highest nutrient limitations on growth was also the most productive.Thus, water, the primary growth factor [11,17], was not limiting at this site, whereas in the other regions, particularly in region 1, the most limiting factor was water.
In the twin plots, LAI increased by 15% and GE increased by 10% compared with the normal plots.The increase in leaf area might explain the increase in productivity.Albaugh et al. [35] found much larger gains, an increase of 101% in the LAI, in response to fertilization in loblolly pine.In Eucalyptus nitens plantations, nitrogen fertilization increased the LAI by up to 3.1 units, i.e., 56% increment [36].
Fertilization increased the LUE and the GE, which was also reported by Binkley et al. [37] for Eucalyptus and for other forest trees, including loblolly pine [38] and Liquidambar styraciflua [39].These responses might be a result of higher photosynthetic rates with the addition of nutrients [40] or caused by less carbon partitioning to shoots [10,41].
Additional fertilization reduced nutrient limitations in the Eucalyptus plantations and increased wood productivity, LAI, and LUE.Thus, we can identify and quantify opportunities to improve the current fertility management of these plantation populations (high probability of response in 85% of the experimental blocks), particularly in region 3.Moreover, with nutrient limitations, wood production was reduced by 5.3 Mg¨ha ´1¨year ´1.We identified that for all sites, organic matter and clay content were related to FR.Additionally, for the intermediate and high fertility sites (region 2 and 3, respectively), potassium and magnesium, and the sum of bases and the silt content, respectively, were related to FR (Table 6).All of these responses together indicate these variables can be used as diagnostic tools for managers to decide where to invest in fertilization.
Water availability was the primary element affecting productivity and potential response to fertilization.When water resources are available, nutrient limitations become more apparent as the primary difference among the three regions.For practical applications, forest managers can use these results to determine which regions should receive priority in analyses of fertilization investment.For example, on a regional scale, fertilization investments would have a much higher return in region 3 compared with the other regions because of the greater potential for response and generation of economic returns from the plantation forests.Our findings do not provide exact information of where the responses would occur, but the main drivers are described in Table 6 (quantity of organic matter, clay content, potassium, magnesium, sum of bases and silt content).Fisher and Binkley [42] proposed a decision support method to minimize the risks.According to the authors, the risks are related to the average and variance of the response, being a high negative risk if the response is low and with a high variance.In our case, the high percentage of plots that had a positive FR (85%) could help managers to support a decision to apply fertilizer.
The aim of this study was to quantify the potential wood productivity and the proportional response to fertilization with a complete nutrient supplement when plants were grown across a wide range of soil types and climatic conditions.Future studies should focus on the separate effects of different nutrients following fertilization to determine the most limiting elements in each environment.

Figure 1 .
Figure 1.Geographical location of plots in three regions of São Paulo state, Brazil.The circles represent the regions and the red areas the experimental blocks.

Figure 1 .
Figure 1.Geographical location of plots in three regions of São Paulo state, Brazil.The circles represent the regions and the red areas the experimental blocks.

Figure 2 .
Figure 2. Correlation between initial volume of twin (TP) and normal (NP) plots.

Figure 3 .Figure 2 .
Figure 3. Current annual increment (CAI) for stems in the twin (TP) and normal plots (NP) for the periods between 2003 and 2006 in the three regions.

Figure 2 .
Figure 2. Correlation between initial volume of twin (TP) and normal (NP) plots.

Figure 3 .Figure 3 .
Figure 3. Current annual increment (CAI) for stems in the twin (TP) and normal plots (NP) for the periods between 2003 and 2006 in the three regions.

Figure 5 .
Figure 5. (a) Average LAI for each region and treatment.Bars with identical letters are not significantly different.Lowercase letters represent the comparisons of treatments and uppercase letters represent the differences between regions (p = 0.05); (b) Correlation between growth efficiency (GE) in twin plots (TP) and normal plots (NP).

Figure 5 .
Figure 5. (a) Average LAI for each region and treatment.Bars with identical letters are not significantly different.Lowercase letters represent the comparisons of treatments and uppercase letters represent the differences between regions (p = 0.05); (b) Correlation between growth efficiency (GE) in twin plots (TP) and normal plots (NP).

Table 1 .
Climatic characteristics and classification (Köppen) of the three study regions.

Table 2 .
Total amount of macronutrients (kg•ha −1 ) applied to the conventional fertilization plots and the extra fertilization twin plots.

Table 1 .
Climatic characteristics and classification (Köppen) of the three study regions.

Table 2 .
Total amount of macronutrients (kg¨ha ´1) applied to the conventional fertilization plots and the extra fertilization twin plots.

Table 4 .
The average and maximum annual increment at 7 years of age (MAI7) for twin and normal plots.

Table 4 .
The average and maximum annual increment at 7 years of age (MAI7) for twin and normal plots.

Table 5 .
Response to fertilization in the three regions.
* Values with different letters are significantly different according to Tukey's test at 5% probability.

Table 5 .
Response to fertilization in the three regions.
* Values with different letters are significantly different according to Tukey's test at 5% probability.

Table 6 .
Pearson's coefficient relating fertilization response and soil variables.

Table 6 .
Pearson's coefficient relating fertilization response and soil variables.