SciELO - Scientific Electronic Library Online

 
vol.66 número6Disponibilidade de nitrogênio, longevidade foliar e mecanismos de conservação de nitrogênio em folhas de espécies arbóreas tropicaisCultivo de Agaricus blazei ss. Heinemann usando diferentes solos como fonte de material para a camada de cobertura índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Artigo

Indicadores

Links relacionados

Compartilhar


Scientia Agricola

versão impressa ISSN 0103-9016

Sci. agric. (Piracicaba, Braz.) vol.66 no.6 Piracicaba nov./dez. 2009

http://dx.doi.org/10.1590/S0103-90162009000600015 

SOILS AND PLANT NUTRITION

 

Hillslope curvature, clay mineralogy, and phosphorus adsorption in an Alfisol cultivated with sugarcane

 

Curvaturas de relevo, mineralogia da argila e adsorção de fósforo em Argissolo cultivado com cana-de-açúcar

 

 

Diogo Mazza BarbieriI; José Marques JúniorII, *; Luis Reynaldo Ferracciú AlleoniIII; Fernando José GarbuioIV; Livia Arantes CamargoI

IUNESP - Programa de Pós-Graduação em Produção Vegetal, Via de Acesso Prof. Paulo Donato Castellane, S/N - 14884-900 - Jaboticabal, SP - Brasil
IIUNESP - Depto. de Solos e Adubos
IIIUSP/ESALQ - Depto. de Ciência do Solo, Av. Pádua Dias, 11 - 13418-900 - Piracicaba, SP - Brasil
IVUSP/ESALQ - Programa de Pós-Graduação em Solos e Nutrição de Plantas

 

 


ABSTRACT

Hillslope curvatures are associated with specific environments that correlate to chemical and mineralogical attributes of soil, so determining specific management zones. Phosphorus is one of the main limiting factors to the development and longevity of sugarcane. The type and the mineralogical constitution of the clay fraction play an important role in the phosphorus (P) adsorption of soil. High proportion of gibbsite (Gb) in soil may be the major responsible for P adsorption. The relationships among spatial variability as a function of hillslope curvature, the proportion of kaolinite (Kt) and Gb, and phosphorus adsorption were evaluated in an Alfisol cultivated with sugarcane. Two plots of 1 ha of a concave and a convex hillslope area were selected and 121 samples were collected in each area. The maximum P adsorption was determined in six samples taken randomly in each area. Data were submitted to descriptive statistical and geostatistical analysis. The lowest average values of available phosphorus were found in the convex area. In this area, the proportion of gibbsite, expressed by the values of the ratio [Gb/(Gb + Ct)] and the values of maximum adsorption capacity of phosphorus were higher than in the concave area.

Key words: Brazil, gibbsite, kaolinite available phosphorus, geostatistics, soil-landscape relationships


RESUMO

As curvaturas do relevo promovem pedoambientes específicos que condicionam os atributos químicos e mineralógicos do solo e podem auxiliar na definição de zonas específicas de manejo. O fósforo (P) é um dos principais elementos limitantes ao desenvolvimento e longevidade do canavial. O teor e a constituição mineralógica da fração argila assumem papel importante na disponibilidade do P, sendo que a gibbsita (Gb), quando presente em altas proporções no solo, pode ser a principal responsável pela sua adsorção e indisponibilidade. Investigaram-se as relações e a variabilidade espacial da adsorção de P e a ocorrência de caulinita (Ct) e gibbsita na fração argila de um Argissolo Vermelho-Amarelo eutrófico originado de rochas areníticas sob diferentes curvaturas do relevo em área sob cultivo de cana-de-açúcar. Duas malhas de 1 ha foram delimitadas numa área côncava e outra área convexa. Foram coletadas 121 amostras em cada área para realização das análises granulométricas, químicas e mineralógicas. A capacidade máxima de adsorção de P foi obtida em seis amostras escolhidas ao acaso em cada área. Os resultados foram submetidos às análises estatísticas descritiva e geoestatística. Os menores valores médios de P disponível encontraram-se na área convexa. Nesta área, a proporção de gibbsita, expressa pelos valores da razão [Gb/(Gb+Ct)] e os valores de capacidade máxima de adsorção de fósforo foram maiores do que na área côncava.

Palavras chave: Brasil, gibbsita, caulinita, fósforo disponível, geoestatística, relação solo-paisagem


 

 

INTRODUCTION

The hillslope curvatures and steepness make up a complex pattern of water and solute transportation in a soil landscape (Marques Júnior & Lepsch, 2000; Franzen et. al., 2006). They are also correlated to specific processes that define the quality, the crystallization, the spatial distribution and the type of the soil clay minerals (Kravchenko & Bullock, 2000; Ghidin et al., 2006). Camargo et al. (2008) verified the influence of small differences of the relief in the crystallinity and distribution of soil clay minerals. Reatto et al. (2008) studied an Oxisols' toposequence in the central highlands of Brazil both, under local and regional scale, and observed higher values for the gibbsite / (kaolinite + gibbsite) in hilly areas in comparison to flatter areas at the local scale of the study.

Similarities between the spatial distribution of the soil clay minerals and their chemical properties may explain the cause/effect relationship between these properties and the relief curvature (Franzen et al., 2006). This relationship is important because it can guide farm management decisions, such as the number of soil samples necessary to characterize an area (Souza et al., 2006) and the required amount of fertilizer.

Reduced P availability in soil is the most common nutritional limitation in agricultural production (Motta et al., 2002; Johnson & Loeppert, 2006), especially in highly weathered soils. Improving the efficiency of phosphate fertilization has been objective of many researchers (Francisco et al., 2007; Nogueira & Cardoso, 2007; Galvani et al., 2008). Most crops in Brazil respond positively to the application of phosphorus. However, most of the supplementary phosphorus becomes unavailable due to adsorption reactions onto mineral colloids, precipitation reactions, or conversion to organic forms (Busato et al., 2005).

The mineralogical composition and content of the clay fraction assumes an important role in the behavior of phosphorus. Al and Fe hydroxides and oxides in the clay fraction of soils are considered primarily responsible for the adsorption of phosphorus (Johnson & Loeppert, 2006; Ulén & Snäll, 2007), and goethite is the most effective form (Motta et al., 2002). Conversely, Motta et al. (2002) and Mesquita Filho & Torrent (1993) pointed out that its contribution to the adsorption of phosphorus can be the equal to or higher than that of goethite when in soils with a high amount of gibbsite. Schaefer et al. (2004) observed that the maximum capacity of P adsorption is correlated better with the gibbsite content than with goethite, demonstrating the importance of Al hydroxides in this phenomenon.

The objectives of the present study were to evaluate: (i) the relationships between phosphorus adsorption and the occurrence of kaolinite and gibbsite in the soil clay fraction from samples of an Alfisol cultivated with sugarcane; and (ii) the spatial variability of several soil attributes as a function of differences in hillslope curvatures.

 

MATERIAL AND METHODS

Soil samples were collected in Catanduva, state of São Paulo, Brazil. It lies in the Paulista Western Plateau (21º05'S; 49º01'W) geomorphic province. The climate of the region is classified as Aw (Köppen's classification); dry in the winter, with average rainfall of 1,350 mm, and an average annual temperature of 23ºC. The primary vegetation of the Catanduva region is savannah seasonal pluvial forest, cultivated with sugarcane for over 20 years. The soil parent material was identified as a sandstone from the Bauru Group, Adamantina formation (IPT, 1981), and the soil is a Typic Hapludalf (Soil Survey Staff, 1999).

Hillslope curvature were classified according to Troeh (1965). Two areas were delimited: one concave and another convex (Figure 1). The agricultural practices of the two areas along the years were identical; both have been cultivated with the sugarcane, cultivar SP80-1842, and the last replanting was in 2003. Fluid fertilization was performed with 700 L ha-1 of 04-12-10 at replanting and 600 L ha-1 10-00-10 in the ratoon. All applications were performed at 20-cm depth. Filter cake and vinasse were not added in these areas.

 

 

Two 1-ha plots were delimited: one in the concave hillslope area and the other in the adjacent convex hillslope area. A spacing grid of 10 × 10 m was used in each of the plots, and the grid nodes were georeferenced using a geodesic GPS. Soil samples were collected at the grid nodes at 0.0-0.2 m and 0.2-0.4 m, totaling 121 sites per plot (Figure 1).

The potential acidity (H+Al) and organic matter content (OM) were determined according to Raij et al. (2001). Exchangeable Ca, Mg, K, and available P were extracted with an ionic exchange resin (Raij et al., 2001). Clay cation exchange capacity (CEC) was calculated by dividing the soil CEC by the amount of clay (Embrapa, 2006). Particle size analysis was performed by the pipette method (Embrapa, 1997).

The kaolinite-gibbsite system was evaluated in 121 samples of the iron-free clay fraction (Mehra & Jackson, 1960) at 0.0-0.2 depth by X-ray diffraction (XRD) - powder method (Camargo et al., 1986) using a diffractometer HZG - 4/B, with a Cu cathode with Ni filter, ka radiation (20 mA, 30 kV), and scanning speed of 1º2q/min. The gibbsite/(gibbsite+kaolinite) [Gb/(Gb+Kt)] ratio was calculated by using the reflection areas of kaolinite [Kt (001)] and gibbsite [Gb (002)]. The mean crystallite dimension (MCD) of Kt and Gb was calculated from the widths at half height (WHH) and the position of the reflections of the minerals Kt (001) and Gb (002) by using Scherrer's equation (Schulze, 1984).

Six nodes were chosen from each grid for obtaining the P maximum adsorption data according to the method described by Casagrande & Camargo (1997). Concentrations of phosphorus (100, 200, 400, 800, and 1,600 mg L1) were added to soils samples, beyond the control treatment, without the addition of phosphorus. The amount of adsorbed P by the soil was calculated by subtracting the value determined in solution from the total amount added. Linearization was performed using Langmuir's equation to obtain the P maximum sorption capacity.

Chemical and mineralogical soil data were submitted to descriptive statistical analysis using the statistical software, SAS (Schlotzhaver & Littell, 1997). Means of the variables were compared by the Student t test. Next, semivariograms were plotted using geostatistical analysis.

 

RESULTS AND DISCUSSION

The available contents of phosphorus were different in the two areas (Table 1). The mean values of the concave and convex areas were 12 and 7 mg dm-3 at 0-0.2 m, and 11 and 6 mg dm-3 at 0.2-0.4-m depth, respectively. These differences are relevant and may influence the decisions regarding soil use and management in both areas. The highest clay content was found in the convex area (Table 1), which may be a contributing factor to reduce P availability. Rolim Netto et al. (2004) found that the most significant correlations for the maximum P sorption capacity were associated with the presence of gibbsite and the clay content.

 

 

The largest soil CEC values were observed in the convex area, which had the least amount of available P (Table 1). On the other hand, the CEC of the clay fractions of both areas were similar, which suggests that the two areas had the same type of clay, although more clay was found in the convex area. So, the larger CEC values for the convex area are a result of its higher contents of clay and organic matter.

The presence of OM favors P desorption to the soil solution, competing with P for sorption sites (Mesquita Filho & Torrent, 1993; Heredia & Cirelli, 2007). Mesquita Filho & Torrent (1993) pointed out that P adsorption can be predicted based on the mineralogical characteristics of the soil when the OM content is low. Comparing the mean areas of the diffractogram peaks of Kt and Gb (Table 2), differences were only found for gibbsite. These results demonstrate that the convex area soil has more gibbsite and that both areas have very similar amounts of kaolinite (Table 2). The WHH of kaolinite varied from 0.33 to 0.83 º2θ in the concave area, and 0.38 to 0.87 º2θ in the convex area. Ghidin et al. (2006) studied a toposequence of Oxisols and found WHH values for kaolinite ranging from 0.55 to 0.85 º2θ. The WHH values for gibbsite ranged from 0.07 to 0.45 º2θ in the concave area, and from 0.09 to 0.33 º2θ in the convex area. Ghidin et al. (2006) found WHH values ranging from 0.28 to 0.39 º2θ for gibbsite.

 

 

The MCD values of kaolinite for both areas were different, in the same way to the WHH values. On the other hand, those values for gibbsite were not different. This demonstrates that kaolinite is better crystallized in the concave area, that is, the crystals have a smaller specific surface area and smaller P sorption capacity than the crystals in the convex area soil. (Leal & Velloso, 1973; Juo & Fox, 1977; Fontes & Weed, 1996)

The mean Gb/(Gb+Kt) ratio of the concave area was 0.06, while that of the convex area was 0.12, and this difference was significant (Student t test). The crystallographic attributes (Table 2) are different for gibbsite, while they are qualitatively different in the two studied areas for kaolinite. Therefore, the convex area has a highest proportion of gibbsite and smaller kaolinite crystals when compared to the concave area. The convex area favors the formation of gibbsite, demonstrating that the relief may indicate the boundaries between the environments and specific areas for the management and transfer of information (Marques Júnior & Lepsch, 2000; Franzen et al., 2006).

The mean and median values of all of the properties were similar, with symmetrical distributions, as confirmed by the asymmetry values close to zero (Tables 1 and 2). The normality of the data predicts stationarity and shows that the so-called proportional effect was not observed; therefore, the data were suited for geostatistics (Isaaks & Srivastava, 1989). According to the classification proposed by Warrick & Nielsen (1980), low variation coefficients were found for the Gb/(Gb+Kt) ratio in the convex and concave areas. The same was observed for the CEC of the convex area at both depth ranges (Tables 1 and 2). The remaining properties presented mean and high coefficients of variation.

All variables presented spatial dependence at the two depth ranges (Table 3). At 0.0-0.2 m, only the convex soil CEC and the clay CEC of the concave area fitted the exponential model. At 0.2-0.4 m, the OM content of both areas and the CEC of the convex area fitted the Gaussian model. The CEC values of the soil and clay, the peak area of kaolinite, the MCD of gibbsite, and the Gb/(Gb+Kt) ratio of the concave area fitted the exponential model better. The remaining properties fitted the spherical model. Several authors obtained spherical model fitting for most of the soil properties (Grego & Vieira, 2005).

 

 

The degree of spatial dependence of the variables that were studied (Table 3) was evaluated by using the classification of Cambardella et al. (1994). At 0.0-0.2-m depth in the concave area, the P properties and the soil CEC presented moderate dependence. In the concave area, at 0.2-0.4 m, the P, OM, clay properties and the CEC of the soil and the kaolinite peak area of the convex area presented moderate spatial dependence. The remaining properties presented strong spatial dependence.

Ranges were found to vary with depth from 22 to 101 m at 0.0-0.2 m, and 18 to 89 at 0.2-0.4 m in the concave area. In the convex area, ranges varied from 20 to 26 m at 0.0-0.2 m and from 15 to 30 m at 0.2-0.4 m (Table 3). The properties presented different spatial dependence ranges, and ranges for the concave area were larger for all properties.

The lowest range values and highest spatial variability of chemical properties were found in the convex area (Table 3), demonstrating the influence of relief on the spatial behavior of these properties (Kravchenko & Bullock, 2000; Souza et al., 2006; Camargo et al., 2008). Mcbratney & Webster (1983) reported that the efficacy of the soil sampling strategy can be improved by using the spatial variability model, which can define specific management zones when associated to the relief forms (Marques Júnior & Lepsch, 2000).

The maximum sorption capacity of phosphorus in the soil from the concave area varied from 416 to 526 mg kg-1, while the soil from the convex area ranged from 454 to 588 mg kg-1, with mean values of 461 and 539 mg kg-1, respectively (Table 4). This confirms that gibbsite extensively affect P adsorption and is more effective than OM (Mesquita Filho & Torrent, 1993). Motta et al. (2002) found that the maximum capacity of sorption of phosphorus varied from 676 to 1,233 mg kg-1 when the Gb/(Gb+Kt) ratio varied from 0.47 to 0.54, respectively.

 

 

Based in the maximum sorption capacity of phosphorus (Table 4) and the amount of available phosphorus (Table 1), one can suggest that gibbsite is primarily responsible for higher P adsorption and lower P availability in the convex than in the concave area because it has greater proportion and better crystallization in this area (Table 2).

 

CONCLUSIONS

Convex hillslope curvatures were more favorable environments for the formation of gibbsite than concave areas. Thus, hillslope curvatures can be used to define the boundaries of the more homogenous areas.

Gibbsite is a determining factor in the maximum sorption capacity of phosphorus and in its smaller availability in the soil of the convex area.

The soil in the hillside area presented greater spatial variability for all of the properties studied as compared to the concave area, indicating that the sampling strategy used in the two areas must be different.

 

REFERENCES

BUSATO, J.G.; CANELLAS, L.P.; VELLOSO, A.C.X. Fósforo num Cambissolo cultivado com cana-de-açúcar por longo tempo. I. Fracionamento sequencial. Revista Brasileira de Ciência do Solo, v.29, p.935-945, 2005.         [ Links ]

CAMARGO, O.A.; MONIZ, A.C.; JORGE, J.A.; VALADARES, J.M.A.S. Métodos de análise química, mineralógica e física de solos do IAC. Campinas: Instituto Agronômico de Campinas, 1986. 94p. (Boletim Técnico, 106).         [ Links ]

CAMARGO, L.A.; MARQUES JUNIOR, J.; PEREIRA, G.T.; HORVAT, R.A. Variabilidade espacial de atributos mineralógicos de um Latossolo sob diferentes formas do relevo: I - mineralogia da fração argila. Revista Brasileira de Ciência do Solo, v.32, p.2269-2277, 2008.         [ Links ]

CAMBARDELLA, C.A.; MOOMAN, T.B.; NOVAK, J.M.; PARKIN, T.B.; KARLEM, D.L.; TURVO, R.F.; KONOPA, A.E. Field scale variability of soil properties in Central Iowa soil. Soil Science Society of America Journal, v.58, p.1501-1511, 1994.         [ Links ]

CASAGRANDE, J.C.; CAMARGO, O.A. Adsorção de fosfato em solos com caráter ácrico validada por um modelo de complexação de superfície. Revista Brasileira de Ciência do Solo, v.21, p.353-360, 1997.         [ Links ]

EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Sistema Brasileiro de Classificação de Solos. 2.ed. Rio de Janeiro: Embrapa, 2006. 412p.         [ Links ]

EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Manual de métodos de análise de solo. 2 ed. Rio de Janeiro: Ministério da Agricultura, 1997. 212p.         [ Links ]

FONTES, M.P.F.; WEED, S.B. Phosphate adsorption by clays from Brazilian Oxisols: relationships with specific surface area and mineralogy. Geoderma, v.72, p.37-51, 1996.         [ Links ]

FRANCISCO, E.A.B.; PROCHNOW, L.I.; TOLEDO, M.C.M.; FERRARI, V.C.; JESUS, S.L. Thermal treatment of aluminous phosphates of the crandallite group and its effect on phosphorus solubility. Scientia Agricola, v.64, p.269-274, 2007.         [ Links ]

FRANZEN, D.W.; NANNA, T.; NORVELL, W.A. A survey of soil attributes in North Dakota by landscape position. Agronomy Journal, v.98, p.1015-1022, 2006.         [ Links ]

GALVANI, R.; HOTTA, L.F.K.; ROSOLEM, C.A. Phosphorus sources and fractions in an Oxisol under no-tilled soybean. Scientia Agricola, v.65, p.415-421, 2008.         [ Links ]

GHIDIN, A.A.; MELO, V.F.; LIMA, V.C.; LIMA, J.M.J.C. Topossequência de latossolos originados de rochas basálticas no Paraná. I. Mineralogia da fração argila. Revista Brasileira de Ciência do Solo, v.30, p.293-306, 2006.         [ Links ]

GREGO, C.R.; VIEIRA, S.R. Variabilidade espacial de propriedades físicas do solo em uma parcela experimental. Revista Brasileira de Ciência do Solo, v.29, p.169-177, 2005.         [ Links ]

HEREDIA, O.S.; CIRELLI, A.F. Environmental risks of increasing phosphorus addition in relation to soil sorption capacity. Geoderma, v.137, p.426-431, 2007.         [ Links ]

INSTITUTO DE PESQUISAS TECNOLOGICAS DO ESTADO DE SÃO PAULO - IPT. Mapa geomorfológico do Estado de São Paulo. Escala - 1:1.000.000. São Paulo, IPT, 1981. v.1, 94p.         [ Links ]

ISAAKS, E.H.; SRIVASTAVA, R.M. An introduction to applied geoestatistics. New York: Oxford University Press, 1989. 561p.         [ Links ]

JOHNSON, S.E.; LOEPPERT R.H. Role of organic acids in phosphate mobilization from iron oxide. Soil Science Society of America Journal, v.7, p.222-234, 2006.         [ Links ]

JUO, A.S.R.; FOX, R.L. Phosphate sorption characteristics of some bench-mark soils of West Africa. Soil Science, v.124, p.370 376, 1977.         [ Links ]

KRAVCHENKO, A.N.; BULLOCK, D.G. Correlation of corn and soybean grain yield with topography and soil properties. Agronomy Journal, v.92, p.75-83, 2000.         [ Links ]

LEAL, J.R.; VELOSO, A.C.X. Phosphate adsorption in latosols under Cerrado vegetation. Pesquisa Agropecuária Brasileira, v.8, p.81-88, 1973.         [ Links ]

MARQUES JÚNIOR, J.; LEPSCH, I.F. Depósitos superficiais neocenozóicos, superfícies geomórficas e solos em Monte Alto, SP. Geociência, v.19, p.226-231, 2000.         [ Links ]

McBRATNEY, A.B.; WEBSTER, R. How many observations are needed for regional estimation of soil properties? Soil Science, v.135, p.177-183, 1983.         [ Links ]

MEHRA, O.P.; JACKSON, M.L. Iron oxide removal from soils and clay by a dithionite-citrate system bulfered with sodium bicarbonate. Clays and Clay Minerals, v.7, p.317-327, 1960.         [ Links ]

MESQUITA FILHO, M.V.; TORRENT, J. Phosphate sorption as related to mineralogy of a hydrosequence of soils from Cerrado Region (Brazil). Geoderma, v.58, p.107-123, 1993.         [ Links ]

MOTTA, P.E.F. Adsorção e formas de fósforo em Latossolos: influencia da mineralogia e histórico de uso. Revista Brasileira de Ciência do Solo, v.26, p.349-359, 2002.         [ Links ]

NOGUEIRA, M.A.; CARDOSO, E.J.B.N. Phosphorus availability changes the internal and external endomycorrhizal colonization and affects symbiotic effectiveness. Scientia Agricola, v.64, p.295-300, 2007.         [ Links ]

REATTO, A.; BRUAND, A.; MARTINS, E.S.; MULLER, F.; SILVA, E.M.; CARVALHO JR, O.A.; BROSSARD, M. Variation of the kaolinite and gibbsite content at regional and local scale in Latosols of the Brazilian Central Plateau. Comptes Rendus Geoscience, v.340, p.741-748, 2008.         [ Links ]

RAIJ, B. van. Análise química para avaliação da fertilidade de solos tropicais. Campinas, Instituto Agronômico, 2001. 285p.         [ Links ]

ROLIM NETO, F.C.; SCHAEFER, C.E.G.R.; COSTA, L.M.; CORRÊA, M.M.; FERNANDES FILHO, E.I.; IBRAIMO, M.M Adsorção de fósforo, superfície específica atributos mineralógicos em solos desenvolvidos de rochas vulcânicas do Alto Paranaíba (MG). Revista Brasileira de Ciência do Solo, v.28, p.953-964, 2004.         [ Links ]

SCHAEFER, C.E.G.R.; GILKES, R.J.; FERNANDES, R.B.A. EDS/SEM study on microaggregates of Brazilian Latosols in relation to P adsorption and clay fraction attributes. Geoderma, v.123, p.1-33, 2004.         [ Links ]

SCHLOTZHAVER, S.D.; LITTELL, R.C. SAS: system for elementary statistical analysis. 2 ed. Cary: SAS, Institute, 1997. 905p.         [ Links ]

SCHULZE, D.G. The influence of aluminium on iron oxides VIII. Unit-cel dimension of Al-substituted of goethites and estimation of Al from them. Clays and Clay Minerals, v.32, p.36-44, 1984.         [ Links ]

SOIL SURVEY STAFF. Soil taxonomy: a basic system of soil classification for making and interpreting soil surveys. 2 ed. Washington, D.C.: USDA, 1999. 869p. (Agriculture Handbook, v.436).         [ Links ]

SOUZA, Z.M.; MARQUES JR, J.; PEREIRA, G.T.; BARBIERI, D.M. Effects of relief shapes on soil chemical properties of an Oxisol cultivated with sugar cane. Scientia Agricola, v.63, p.161-168, 2006.         [ Links ]

TROEH, F.R. Landform equations fitted to contour maps. American Journal of Science, v.263, p.616-627, 1965.         [ Links ]

ULÉN, B.; SNÄLL, S. Forms and retention of phosphorus in an illite-clay soil profile with a history of fertilisation with pig manure and mineral fertilizers. Geoderma, v.137, p.455-465, 2007.         [ Links ]

WARRICK, A.W.; NIELSEN, D.R. Spatial variability of soil physical properties in the field. In: HILLEL, D. (Ed.) Applications of soil physics. New York: Academic Press, 1980. chap.2, p.319-344.         [ Links ]

 

 

Received November 11, 2008
Accepted Agosto 27, 2009

 

 

* Corresponding author <marques@fcav.unesp.br>