SciELO - Scientific Electronic Library Online

 
vol.36 issue3Onion yield and bulb conservation as affected by rates and splitting of side-dressed nitrogenDistribution of Nitrogen Ammonium Sulfate (15N) soil-plant system in a no-tillage crop succession author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Revista Brasileira de Ciência do Solo

Print version ISSN 0100-0683

Rev. Bras. Ciênc. Solo vol.36 no.3 Viçosa May/June 2012

http://dx.doi.org/10.1590/S0100-06832012000300018 

DIVISÃO 3 - USO E MANEJO DO SOLO
COMISSÃO 3.1 - FERTILIDADE DO SOLO E NUTRIÇÃO DE PLANTAS

 

Recommendation of soil fertility levels for willow in the southern highlands of Santa Catarina1

 

Classes recomendáveis de fertilidade do solo para a produção de vime no Planalto Sul Catarinense

 

 

Tássio Dresch RechI; Flávio ZanetteII; João Claudio ZanattaIII; Névio João NuernbergIV; Dieter BrandesV

IPesquisador da Epagri/Estação Experimental de Lages. R. João José Godinho, s/n. Br. Morro do Posto, Caixa Postal 181, CEP 88502-970 Lages (SC). E-mail: tassior@epagri.sc.gov.br
IIProfessor do Programa de Pós-Graduação em Produção Vegetal, Universidade Federal do Paraná - UFPR. R. dos Funcionários 1540, Br. Juvevê, Caixa Potal 19061, CEP 81531-990 Curitiba (PR). E-mail: flazan@ufpr.br
IIIPesquisador da Epagri/Estação Experimental de Lages. E-mail: zanatta@epagri.sc.gov.br
IVPesquisador aposentado. E-mail: nevio.joao@gmail.com
VPesquisador aposentado. E-mail: dieterb@brturbo.com

 

 


SUMMARY

The species Salix x rubens is being grown on the Southern Plateau of Santa Catarina since the 1940s, but so far the soil fertility requirements of the crop have not been assessed. This study is the first to evaluate the production profile of willow plantations in this region, based on the modified method of Summer & Farina (1986), for the recommendation of fertility levels for willow. By this method, based on the law of Minimum and of Maximum for willow production for the conditions on the Southern Plateau of Santa Catarina, the following ranges could be recommended: pH: 5.0-6.5; P: 12-89 mg dm-3; Mg: 3.2-7.5 mg; Zn: 5.0-8.3 mg dm-3; Cu: 0.8-4.6 mg dm-3; and Mn; 20-164 mg dm-3. The Ca/Mg ratio should be between 1.2 and 2.9. For K and Ca only the lower (sufficiency level), but not the upper threshold (excess) was established, with respectively 114 mg dm-3 and 5.3 cmolc dm-3. It was also possible to determine the upper threshold for Al and the Al/Ca ratio, i.e., 1.7 cmolc dm-3 and 0.28, respectively. For maximum yields, the clay in the soil surface layer should be below 320 g dm-3.

Index terms: Salix x rubens; nutritional diagnosis; cropping system.


RESUMO

A espécie Salix x rubens é cultivada no Planalto Sul Catarinense desde a década de 1940, porém ainda não estão estabelecidas suas exigências quanto da fertilidade de solo. O presente trabalho é o primeiro a avaliar o perfil produtivo dos vimais dessa região, com a aplicação do método adaptado de Summer & Farina (1986), para a recomendação de fertilidade para o cultivo do vimeiro. Esse método, baseado na Lei do Mínimo e na Lei do Máximo, permitiu estabelecer as seguintes faixas recomendadas para a produção de vime nas condições do Planalto Sul Catarinense: pH - entre 5,0 e 6,5; P - 12 a 89 mg dm-3; Mg - entre 3,2 e 7,5; Zn - 5,0 a 8,3 mg dm-3; Cu - 0,8 a 4,6 mg dm-3; e Mn - 20 a 164 mg dm-3. A relação Ca/Mg deve estar entre 1,2 e 2,9. Para o K e o Ca foi possível estabelecer apenas o nível de mínimo, mas não de excesso, respectivamente de 114 mg dm-3 e 5,3 cmolc dm-3. Foi possível ainda definir o nível de excesso para Al e para a relação Al/Ca: 1,7 cmolc dm-3 e 0,28, respectivamente. O solo deve apresentar ainda teores de argila na camada superficial inferiores a 320 g dm-3 para que a expectativa de produção seja máxima.

Termos de indexação: Salix x Rubens, diagnóstico nutricional, sistema de cultivo.


 

 

INTRODUCTION

Since the 1940s or even before, willow has been grown on the Southern Plateau of Santa Catarina, although cultivation for commercial purposes only began in the 1970s (Epagri, 2008). So far, the crop requirements for soil fertility on the Santa Catarina Plateau were not determined. Even in countries where willow cultivation is a tradition, e.g., in the UK, little information on the pH and the nutritional requirements of the commonly used varieties is available (Heaton et al., 2000).

Consequently, cultural practices related to liming and fertilization, planting density, rod length and diameter, soil tillage and others are being used empirically. This may affect the yield, willow quality and longevity of plantations.

The determination of the appropriate nutrient levels in soil and plant tissue is usually based on dose-response curve experiments, testing increasing doses of one or more nutrients (calibration tests) (Raij, 1981). However, these methods are somewhat limited, e.g., experiments have to be carried out with combinations of two or three nutrients together to assess individual responses and interactions between them, with restricted ability to assess other co-called controlled/constant factors or interactions among more than three factors (Sims, 2000; Summer & Farina, 1986).

A graphical method proposed by Cate Jr. & Nelson (1965) separates soils in two classes, i.e., soils with high or low probability response to a particular nutrient. Although the arbitrariness in the definition of these classes can be contained by a statistical procedure proposed by the same authors (Cate Jr. & Nelson, 1971), the results allow no more than a first approximation (Raij, 1981).

The most commonly used procedure to determine critical nutrient levels in the soil is known as calibration (Raij, 1981). The greatest difficulty in comparing results from different calibration experiments is the interaction between the controlled variables, maintained constant, and the non- controlled with the tested nutrients. Two or more variables interact when their individual influence is modified by the presence, concentration, intensity or activity of the other factors. When the joint action of factors results in a growth response that is, for example, greater than the sum of their individual effects, the interaction is positive (synergic), when it is lower than the sum of individual effects, the interaction is negative (antagonistic). When the sum of the individual effects of factors is similar to the combined action, there is no interaction (Summer & Farina, 1986, Wilkinson et al., 2000). This process, however requires extensive experimentation. Thus, when there is little or no information about the fertility requirement of a crop, the search for alternatives to determine the critical levels becomes challenging. In a new context, the Law of the Minimum or the limiting factor principle, or minimum law of Liebig, can be used, not for the isolation and control of factors, but for the identification of interactions and complex responses (Summer & Farina, 1986; Wilkinson et al., 2000).

The Law of the Minimum determines that the growth of a plant is limited by that soil nutrient found in the lowest proportion of the plant´s requirement. The Minimum principle is applied under conditions where the supply of energy, minerals and other factors is balanced, near-optimum or unlimited. Once the deficient nutrient is supplied, it induces additional growth, until another nutrient or factor becomes limiting and so on (Raij, 1981; Sims, 2000, Wilkinson et al., 2000).

The application of the Law of the Minimum is complex because there are often several production-limiting nutrients or factors and interactions. Moreover, as defined by the Law of Mitscherlich, biological responses are normally not linear (Raij, 1981). Thus, the application of a nutrient to the soil in sufficient quantity, as described by Liebig, results in a curvilinear response to additional doses of that nutrient until the production becomes maximum, then a "plateau" is reached, as described by Mitscherlich (Malavolta, 1976).

Voisin (1973) emphasizes the role of interactions and claims that the lack of one nutrient in the soil reduces the effectiveness of other nutrients and, consequently, of crop yields. Therefore, each nutrient is most effective when the others are near-optimum. This author also stated that the excess of one nutrient in the soil reduces the efficiency of others and may therefore decrease yields, which is confirmed by the parabolic shape of the response curves. As foreseen by the Mitscherlich Law, in its first approximation, the yield increases decrease with additional doses of a nutrient until reaching a maximum, a plateau, after which the production decreases. This reduction after the maximum was not predicted by Mitscherlich. For this case, Voisin (1973) proposed the Maximum law, based on the understanding that the excess of a nutrient limits or reduces yields, due to nutrient imbalances or toxic effects.

However, Wallace (1990) proposed the Law of the Maximum in contrast to the Law of the Minimum, somewhat different from the focus and application of the proposal of Voisin. Wallace's Law determines that when the needs for each factor involved in the process are fully met, the rate of this process can reach its full potential, which is greater than the sum of its parts, due to sequential positive interactions. In other words, for each nutrient content or value of a production factor there is a maximum possible yield, which occurs when all other factors and interactions are optimal for the content of this nutrient. In other words, for each level or value of a production factor there is a maximum possible production. This production will only be achieved when all interactions and levels/values of other production factors are not less restrictive than the factor in analysis. Thus, the maximum yield obtained with a certain content or amount of a nutrient or production factor represents the constraint condition for this factor only. All yields below the maximum represent another factor or interaction contributing to the limitation.

Based on these concepts, Summer & Farina (1986) proposed an interpretation model, by which the increase of the levels of a certain nutrient "X" in the soil or of the ratio "X/Y" in the plant tissue up to a maximum critical value would induce yield increases. This region of plant response is interpreted as the region of X insufficiency or of Y excess. The region in which the increases in the levels of this factor X do not affect the yield, was called region of nutrient balance by Summer and Farina. Increases in the levels of the factor, beyond the region of balance, result in a progressive yield decline. This area is defined as the X excess or Y insufficiency region.

This method allows the identification of regions of insufficiency, balance and excess levels and ratios between nutrients based on yield graphs in response to nutrient contents or the ratio between nutrients studied (Wilkinson et al., 2000).

The purpose of this study was to recommend soil fertility levels, based on field data from willow plantations on the southern plateau of Santa Catarina, using the method proposed by Summer & Farina (1986). This is the first study addressing the recommendation of soil fertility levels for willow.

 

MATERIALS AND METHODS

In 2000, 48 sampling units were installed on willow plantations in the counties Bocaina do Sul, Bom Retiro, Lages, Rio Rufino, and Urubici, in the state of Santa Catarina. Plantation segments with no dead or diseased plants were selected as sampling units. The size of each sample unit was adjusted to the spacing, the possibilities of sampling and the availability of the producer; the goal was to evaluate at least six trees and 1 m2 per unit, so the variation was 1-11 m² and 3-27 trees. The branches of each plot were weighed and sampled to determine the moisture content. The samples were oven-dried at 60 °C for 72 h.

In each sampling unit, the water table depth, the presence of a hydromorphic horizon to a depth of 120 cm, and the presence of soil layers impermeable to root penetration were observed.

Soil samples were collected from the 0-20 cm layer, mixed to a composite sample from nine sampling points along a diagonal drawn between the willow rows, of which the first and last points were at the end of opposite rows and the others distributed evenly in-between. Soil samples were subjected to routine analyses (pHwater, P, K, Al3+, Ca2+, and Mg2+) and micronutrients Zn2+, Cu2+ and Mn3+, in soil laboratories of Epagri (Empresa de Pesquisa Agropecuária e Extensão Rural de Santa Catarina, 2008) in Chapecó and São Joaquim, and a laboratory of animal nutrition of Epagri/EELages, according to routine methods of the official network of laboratories in Southern Brazil (Tedesco et al., 1995).

In 2001 and 2002, the procedure of soil sampling and analysis and willow production was repeated, wherever possible on the same plots as in the previous year. The plots that had been harvested or could not be repeated for any other reason, were replaced. At each soil sampling, the willow branches in the plot were cut at the trunk and weighed.

The data were analyzed for normality using the Kolmogorov-Smirnov test (p < 5 %) (Johnson & Wichern, 1998).

Based on the method proposed by Summer & Farina (1986), the soil nutrient levels and pH values as well as nutrient levels in plant tissue were diagrammed with the respective values of willow yield; each variable was placed on the x-axis, and the willow yield on the y-axis (Figure 1).

 

 

By the original method of Summer & Farina (1986), the cloud of points in the graph of each variable is circumscribed with a line joining the points of maximum production for each value on the x-axis. However, considering the limited information available to evaluate willow requirements in this study, two straight regression lines were adopted as alternative to minimize the error or characterize the inconsistency of the outline of the points of maximum. The ascending line joins the points of maximal responses for each value on the abscissa in the X insufficiency region, where increasing values on the x-axis result in increases in willow yield, on the y-axis (Figure 1). The descending line is to join the points of maximum responses for each value on the abscissa, - X excess region - where increasing values on the x-axis results in decreases in willow yield, on the y-axis.

For each variable evaluated in this study (pHwater, P, K, Al3+, Ca2+, Mg2+, Zn2+, Cu2+ , and Mn3+), a graph was created with the variable on the x-axis and willow yield on the y-axis. In each graph, a straight line on the left - ascending line - was generated from the leftmost value on the yield graph and the value immediately to the right and above it. The values above and to the right of this segment of the line were gradually added to or removed from the equation in order to obtain the regression with the most significant coefficient of determination R2 (Table of significance of R2 values, Steel & Torrie, 1960), without disregarding values to the left and above the generated straight line in the calculation. The absence of significant R2 in at least 95 % with more than two points in the equation was considered as no adjustment.

To obtain the regression line on the right side of the data set - excess region - the procedure of establishing the descending line was the same. That is, a straight right - descending line - was generated from the rightmost value in the yield graph and the value immediate to the left and above it. The values above and left of that line segment were gradually added to or removed from the equation to obtain the regression with the most significant coefficient of determination R2 (Significance table of R2 values, Steel & Torrie, 1960) , without disregarding values to right and above the generated straight line in the calculation. The absence of significant R2 in at least 95 % with more than two points in the equation was considered as no adjustment.

Once the ascending and descending straight lines of the target variable were established, the crossing point of both lines was calculated, i.e., the maximum point of the variable. This procedure was repeated for each variable. The value of the studied factor required to achieve the maximum point was considered the optimum value of the factor. That is: ya = yd , where: ya = aax  + ba and ya = adx  + bd: where: ya = willow yield; aa = slope of the ascending line, ba = constant of the ascending line, ad = slope of the descending line, bd = constant of the descending line; we have:

aax + ba  = adx + bd; so that: xoptimum = bd  - ba / aa - ad

In the cases where one of the lines was not consistent, i.e., where only a trend to the right or left was observed in the data set, the maximum yield observed in the field was used to calculate the optimum point of the target factor under study.

In addition to the optimal value, a critical range was also defined, limited by the amount of the factor under study needed for yields equivalent to the maximum production observed minus a standard deviation. The value in the ascending line was established as sufficiency threshold, while the resulting value of the descending line was called the excess threshold.

Where: ya = aax + ba and ya = adx + bd; for: ymaximum-Σ ; where: ymaximum = maximum yield, Σ = standard deviation, so we have: xsufficiency = (ymaximum - Σ - ba)/aa; and: xexcess = (ymaximum - Σ - ba)/aa; where: xsufficiency = sufficiency threshold of the factor under study, and: xexcess = excess threshold of the factor under study. Thus, for every variable in the study a balance/sufficiency range was established, where the minimum is the sufficiency value and the maximum the excess value. Values of the target trait below this range in the willow plantations were considered insufficient, and those with values beyond the range in excess for each of the study variables.

 

RESULTS AND DISCUSSION

The willow yield was responsive to changes in soil pH (Figure 2). The optimal pH values of 5.8 and pH ranges from 5.0 to 6.5 are compatible with the recommendations of Abrahamson et al. (2002) for willow cultivation aiming at biomass production and with the recommendations of the SBCS (2004) for forest species such as black wattle. Since the soils of the region are acidic (Ernani et al., 2000), very high lime doses must be used to adjust the conditions of the soils in the willow plantations.

 

 

The highest willow yields were observed in a soil P availability range of 12-89 mg dm-3, with a lower limit similar to the range considered sufficient by the SBCS (2004) for "class three" of clay content, i.e., valid for sandy, as well as for medium or clay soils (Figure 3). The upper limit is much higher than the threshold of SBCS (2004) for most crops, however, economic and environmental aspects most likely restrict fertilization of crops(Keoleian & Volk, 2005; Oldham, 2008, Sharpley et al., 2010). This wide range between lower and upper thresholds when another variable interferes with the response, in this case, soil texture and the determination of available P by Mehlich, as observed by Bortolon & Gianello (2008), may be an undesirable feature of the method.

 

 

Potassium availability was less limiting than would be expected, as farmers are not used to applying K or organic fertilizers (Figure 2). The extraction at harvest of this nutrient was estimated at 29 kg ha-1 year for average dry matter yields of 10 t ha-1 for the regional conditions of this study (Rech et al., 2005). Adegbidi et al. (2001) found a removal of 7 to 44 kg ha-1 K in studies conducted in New York State, USA.

The soils had an average K content of 129 mg dm-3. As similarly observed for P, the minimum level of K sufficiency (114 mg dm-3) was similar to the high and very high levels of the SBCS recommendations (SBCS, 2004) for soil with medium CEC (5.1-15.0 cmolc dm-3 ). However, it was not possible to establish the upper threshold for the surface layer, probably due to the low density of samples with high enough K levels to limit willow yields, be it by nutritional imbalance, or by toxicity or salinity.

The mean Ca content in the soil was 8.4 and the sufficiency threshold 5.3 cmol dm-3 for willow production. This value is similar to the limit of the sufficiency range of the SBCS classification (SBCS, 2004). It is also compatible with the levels reported by Adegbidi et al. (2003) and Labrecque et al. (1998) in experimental areas with willow. Similarly as for K, no upper threshold could be established for Ca either.

For Mg it was possible to establish a sufficiency range of 3.2-7.5 cmol dm-3, while the mean value was 4.4 cmol dm-3 (Figure 2). For the Ca/Mg ratio, the average was 2.3, while the sufficiency range was determined from 1.2 to 2.9. Willow yields were higher in soils with narrower Ca/Mg ratios (Figure 2) than in soils growing Salix for biomass production (Labrecque et al., 1998). Still, the value was within an adequate range for many agricultural crops (Camberato & Pan, 2000). This narrow Ca/Mg ratio does not seem to be a problem, because the limestone used in the region is usually dolomitic lime, and therefore, appropriate for the plant requirements. Perhaps for this reason, only few of the sampled plantations were classified in the insufficiency range.

Although willow yields were higher in the absence of Al, the level tolerated for yields of over 90 % of the maximum were obtained with Al contents of up to 1.7 cmol dm-3 (Figure 2). However, Cronan & Grigal (1995) claimed that the relationship between Ca and Al can be a more appropriate indicator for the evaluation of stress in the forest development. Based on the data and procedure adopted in this study, the maximum value of Al/Ca was 0.28. Above this threshold, the Al/Ca ratio may be the main yield-limiting factor.

For the micronutrients (Figure 3), the optimum Zn content was determined at 6.8 mg dm-3 and the limits of the sufficiency range 5.0 mg and 8.3 mg dm-3. The average Zn content found in willow plantations was 5.8 mg dm-3. These values are high, compared to the classification of the (SBCS, 2004). However, they are similar to levels reported by Labrecque et al. (1998) in soils of willow plantations for biomass production which were assessed for organic fertilization.

The optimal Cu amount determined based on willow yield was 2.8 mg dm-3, which is considered high (SBCS, 2004). Even the lower threshold of the sufficiency range (0.8-4.6 mg dm-3), was considered high according to this classification. On the other hand, Labrecque et al. (1998) reported higher Cu than observed in this study.

The optimal level of Mn, 34.0 mg dm-3, is far higher than the ranges defined by the classification of SBCS (2004), so that the lower sufficiency limit (20.0-163.4 mg dm-3) is almost 10 times greater than the lower limit of the upper class of that classification (Figure 3) and higher than the levels presented by Labrecque et al. (1998).

The average clay content in the samples was 270 g kg-1, in both soil layers. It was not possible to determine a threshold of minimum clay content. This may be the result of the presence of organic soils and areas where there is constant replenishment of nutrients through water flow, situations in which the contribution of clay minerals to nutrient replacement may be negligible. However, soils with higher clay content seem limiting, perhaps due to the difficulties of drainage and soil aeration. According to the model, soils with a clay content > 320 g dm-3 could limit willow yields.

 

CONCLUSIONS

1. The application of the modified method of Summer & Farina (1986) allows the establishment of recommendations of: a) levels of sufficiency and excess for pH, P, Mg, Zn, Cu, and Mn and the Ca/Mg ratio, b) sufficiency but not excess levels of K and Ca, c) the upper threshold for Al, Al/Ca ratio, and for clay contents.

2. For willow production, soils should have a pH between 5.0 and 6.5 and the following nutrient contents: P 12-89 mg dm-3, K > 114 mg dm-3; Ca > 5.3 cmolc dm-3; Zn 5.0-8.3 mg dm-3; Cu 0.8-4.6 mg dm-3; and Mn 20-164 mg dm-3. The level of Al should not be higher than 1.73 cmolc dm-3, the Ca/Mg ratio should be between 1.2 and 2.9 and the Al/Ca ratio < 0.28. The soil clay content should be < 320 g dm-3.

 

ACKNOWLEDGEMENTS

We thank the agricultural extension workers of EPAGRI (Empresa de Pesquisa Agropecuária e Extensão Rural de Santa Catarina): Antônio Edu Arruda, Saulo Poffo, Pedro Donizete de Souza, Evaldo Roberto Schlemper, and Donizete Cruz de Souza for their help in finding partner producers.

 

LITERATURE CITED

ABRAHAMSON, L.P.; VOLK, R.F.; WHITE, E.H. & BALLARD, J.L. Willow Biomass Producer's Handbook. Short-Rotation Woody Crops Program, Sunny College of Environmental Science and Forestry. Syracuse, 2002. Available at: <http://www.esf.edu/willow/PDF/2001%20finalhandbook.pdf> Accessed: 09 May 2011.         [ Links ]

ADEGBIDI, H.G.; VOLK, T.A.; WHITE, E.H.; ABRAHAMSON, L.P.; BRIGGS, R.D. & BICKELHAUPT, D.H. Biomass and nutrient removal by willow clones in experimental bioenergy plantations in New York State. Biomass Bioenergy, 20:399-411, 2001.         [ Links ]

ADEGBIDI, H.G.; BRIGGS, R.D.; VOLK, T.A.; WHITE, E.H. & ABRAHAMSON, L.P. Effect of organic amendments and slow-release nitrogen fertilizer on willow biomass production and soil chemical characteristics. Biomass Bioenergy, 25:389-398, 2003.         [ Links ]

BORTOLON, L. & GIANELLO, C. Interpretation of the avalilabre phosphorus by Mehlich-1 and Mehlich-3 extractant solution in soils of the Rio Grande do Sul State, Brazil. R. Bras. Ci. Solo, 32:2751-2756, 2008.         [ Links ]

CAMBERATO, J.J. & PAN, W.L. Biovailability of calcium, magnesium, and sulfur. In: SUMMER, M., ed. Handbook of soil science. Boca Raton, CRC Press, 2000. p. D53-D70.         [ Links ]

CATE Jr, R. B. & NELSON, L.A. A rapid method for correlation of soil test analysis with plant response data. Raleigh, International Soil Testing, 1965. 24p.         [ Links ]

CATE Jr, R.B. & NELSON, L.A. A simple statistical procedure for partitioning soil test correlation data into two classes. Soil Sci. Am. Proc., 35:658-659, 1971.         [ Links ]

CRONAN, C.S. & GRIGAL, D.F. Use of calcium/aluminum ratios as indicators of stress in forest ecosystems. J. Environ. Qual., 24:209-226, 1995.         [ Links ]

EMPRESA DE PESQUISA AGROPECUÁRIA E EXTENSÃO RURAL DE SANTA CATARINA - EPAGRI. O cultivo do vime. 2008. Available at <http://www.epagri.sc.gov.br/index.php?option=com_content&view=article&id=207:o-cultivo-do-vime&catid=28:projetos-destaqe&Itemid =43>. Accessed: 08 May 2011.         [ Links ]

ERNANI, P.R.; NASCIMENTO, J.A.L.; CAMPOS, M.L. & CAMILLO, R.J. Influência da combinação de fósforo e calcário no rendimento de milho. R. Bras. Ci. Solo, 24:537-544, 2000.         [ Links ]

HEATON, R.J.; RANDERSON, P.F. & SLATER, F.M. The silviculture, nutrition and economics of short rotation willow coppice in the uplands of mid-Wales. 2000. Cardiff University for ETSU New and Renewable Energy Progamme. available at: <http://www.ecodyfi.org.uk/pdf/coppice.pdf>. Accessed 08 May 2011.         [ Links ]

JONHSON, R.A. & WICHERN, D.W. Applied multivariate statistical analysis. Madison, Prentice Hall International, 1998. 607p.         [ Links ]

KEOLEIAN, G.A. & VOLK, T.A. Renewable energy from willow biomass crops: Life cycle energy, environmental and economic performance. Crit. Rev. Plant Sci., 24:385-406, 2005.         [ Links ]

LABRECQUE, M.; TEODORESCU, T.I. & DAIGLE, S. Early performance and nutrition of two willow species in short-rotation intensive culture fertilized with wastewater sludge and impact on the soil characteristics. Can. J. For. Res., 28:1621-1638, 1998.         [ Links ]

MALAVOLTA, E. Manual de química agrícola. São Paulo, Agronômica Ceres, 1976. 528p.         [ Links ]

OLDHAM, L. Phosphorus in Mississippi soils. 2008. Available at: <http://msucares.com/pubs/infosheets/is0871.pdf>. Accessed: 08 May 2011.         [ Links ]

RAIJ, B.van. Avaliação da fertilidade do solo. Piracicaba, Instituto da Potassa e Fosfato/Instituto Internacional da Potassa, 1981. 142p.         [ Links ]

RECH, T.D.; NUERNBERG, N.; BRANDES, D.; ARRUDA, A.E. & ZANETTE, F. Remoção de nutrientes pelo cultivo do vime no Vale do Rio Canoas, SC. In: CONGRESSO INTERNACIONAL DE DESENVOLVIMENTO RURAL E AGROINDÚSTRIA FAMILIAR, 1., São Luiz Gonzaga-RS, 2005. Anais... São Luiz Gonzaga, 2005. p.263-267.         [ Links ]

SOCIEDADE BRASILEIRA DE CIÊNCIA DO SOLO - SBCS. Manual de adubação e calagem para os Estados do Rio Grande do Sul e de Santa Catarina. 10.ed. Porto Alegre, SBCS/ Núcleo Regional Sul; Comissão de Química e Fertilidade do Solo - RS/SC, 2004. 400p.         [ Links ]

SHARPLEY, A.; DANIELS, M.; VANDEVENDER, K. & SLATON, N. Soil phosphorus: Management and recommendations. 2010. Available at <http://www.uaex.edu/Other_Areas/publications/pdf/FSA-1029.pdf>. Accessed: 08 May 2011.         [ Links ]

SIMS, J.T. Soil fertility evaluation. In: SUMMER, M. ed. Handbook of soil science. Boca Raton, CRC Press, 2000. p.D113-D153.         [ Links ]

STEEL, R.G.D. & TORRIE, J.H. Principles and procedures of statistics. New York, McGraw, 1960.         [ Links ]

SUMMER, M.E. & FARINA, M.P.W. Phosphorus interactions with other nutrients and lime in field cropping systems. Adv. Soil Sci., 5:201-236, 1986.         [ Links ]

TEDESCO, M.J.; VOLKWEISS, S.J. & BOMMEN, H. Análises de solo, plantas e outros materiais. Porto Alegre, Universidade Federal do Rio Grande do Sul, 1995, 174p.         [ Links ]

VOISIN, A. Adubos - Novas leis científicas de sua aplicação. São Paulo, MestreJou, 1973. 130p.         [ Links ]

WALLACE, Z. Crop improvement through multi-disciplinary approaches to different types of stresses: Law of maximum. J. Plant Nutr., 13:313-325, 1990.         [ Links ]

WILKINSON, S.R.; GRUNES, D.L. & SUMMER, M.E. Nutrient interactions in soil and plant nutrition. In: SUMMER, M., ed. Handbook of soil science. Boca Raton, CRC Press, 2000. p.D89-D112.         [ Links ]

 

 

1 Trabalho financiado pela Epagri, Sebrae e Fapesc, extraído da Tese de Doutorado do primeiro autor junto ao Programa de Pós-graduação em Produção Vegetal, Universidade Federal do Paraná - UFPR. Received for publication in May 23, 2011 and approved in February 28, 2012.