Hillslope curvature and P adsorption in an Alfisol 819 Sci. Agric. (Piracicaba, Braz.), v.66, n.6, p.819-826, November/December 2009 HILLSLOPE CURVATURE, CLAY MINERALOGY, AND PHOSPHORUS ADSORPTION IN AN ALFISOL CULTIVATED WITH SUGARCANE Diogo Mazza Barbieri1; José Marques Júnior2*; Luis Reynaldo Ferracciú Alleoni3; Fernando José Garbuio4; Livia Arantes Camargo1 1 UNESP - Programa de Pós-Graduação em Produção Vegetal, Via de Acesso Prof. Paulo Donato Castellane, S/ N - 14884-900 - Jaboticabal, SP - Brasil. 2 UNESP - Depto. de Solos e Adubos. 3 USP/ESALQ - Depto. de Ciência do Solo, Av. Pádua Dias, 11 - 13418-900 - Piracicaba, SP - Brasil. 4 USP/ESALQ - Programa de Pós-Graduação em Solos e Nutrição de Plantas. *Corresponding author 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 CURVATURAS DE RELEVO, MINERALOGIA DA ARGILA E ADSORÇÃO DE FÓSFORO EM ARGISSOLO CULTIVADO COM CANA-DE-AÇÚCAR 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 spe- Barbieri et al.820 Sci. Agric. (Piracicaba, Braz.), v.66, n.6, p.819-826, November/December 2009 cific processes that define the quality, the crystalliza- tion, 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 high- lands 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 num- ber of soil samples necessary to characterize an area (Souza et al., 2006) and the required amount of fertil- izer. 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 re- searchers (Francisco et al., 2007; Nogueira & Cardoso, 2007; Galvani et al., 2008). Most crops in Brazil re- spond positively to the application of phosphorus. How- ever, most of the supplementary phosphorus becomes unavailable due to adsorption reactions onto mineral colloids, precipitation reactions, or conversion to or- ganic forms (Busato et al., 2005). The mineralogical composition and content of the clay fraction assumes an important role in the behav- ior 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). Con- versely, Motta et al. (2002) and Mesquita Filho & Tor- rent (1993) pointed out that its contribution to the ad- sorption 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, demonstrat- ing the importance of Al hydroxides in this phenomenon. The objectives of the present study were to evalu- ate: (i) the relationships between phosphorus adsorp- tion 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 sev- eral 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 Pla- teau (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 sug- arcane 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 con- cave and another convex (Figure 1). The agricultural practices of the two areas along the years were iden- tical; 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 con- tent (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 cal- culated 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 & Jack- son, 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). Hillslope curvature and P adsorption in an Alfisol 821 Sci. Agric. (Piracicaba, Braz.), v.66, n.6, p.819-826, November/December 2009 Six nodes were chosen from each grid for obtain- ing 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 L–1) were added to soils samples, be- yond the control treatment, without the addition of phosphorus. The amount of adsorbed P by the soil was calculated by subtracting the value determined in so- lution 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 submit- ted to descriptive statistical analysis using the statisti- cal 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 differ- ent 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 manage- ment in both areas. The highest clay content was found in the convex area (Table 1), which may be a con- tributing 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 min- eralogical characteristics of the soil when the OM content is low. Comparing the mean areas of the diffractogram peaks of Kt and Gb (Table 2), differ- ences 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 °2q in the concave area, and 0.38 to 0.87 °2q 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 °2q. The WHH values for gibbsite ranged from 0.07 to 0.45 °2q in the concave area, and from 0.09 to 0.33 Figure 1 - Location and digital elevation model showing the two hillslope curvatures divisions (Troeh, 1965). The arrows simulate the surface water flow directions. Barbieri et al.822 Sci. Agric. (Piracicaba, Braz.), v.66, n.6, p.819-826, November/December 2009 °2q in the convex area. Ghidin et al. (2006) found WHH values ranging from 0.28 to 0.39 °2q 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 differ- ent. This demonstrates that kaolinite is better crystal- lized in the concave area, that is, the crystals have a smaller specific surface area and smaller P sorption ca- pacity 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 ka- olinite crystals when compared to the concave area. The convex area favors the formation of gibbsite, dem- onstrating 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 proper- ties 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 con- cave 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 etubittA epolslliH erutavruc naeM naideM niM xaM p1 tneiciffeoC )%(.raV yrtemmysA sisotruK m2.0-0.0 mdgm(P –3) evacnoC a21 11 4 52 70.0 14 46.0 32.0- xevnoC b7 7 1 41 51.0 15 10.0 45.0- mdg(MO –3) evacnoC b31 31 9 02 51.0 81 44.0 60.0- xevnoC a51 41 01 02 51.0 41 61.0 11.0 gkg(yalC –1) evacnoC b981 581 08 523 51.0 72 53.0 72.0- xevnoC a032 922 831 523 51.0 61 63.0 33.0 CEClioS mdclomm( –3) evacnoC b25 35 53 66 51.0 31 61.0- 45.0- xevnoC a16 16 74 76 51.0 01 04.0 30.0 CECyalC mdclomm( –3) evacnoC a572 682 302 734 30.0 93 82.0 80.0 xevnoC a562 662 602 043 51.0 23 03.0 11.0 4.0-2.0 m mdgm(P –3) evacnoC a11 01 3 42 30.0 64 96.0 1.0- xevnoC b6 0.6 1 51 51.0 45 02.0 84.0- mdg(MO –3) evacnoC b21 21 6 91 51.0 22 91.0 80.0- xevnoC a41 41 8 81 51.0 41 80.0 400.0 gkg(yalC 1– ) evacnoC b512 002 101 963 10.0 03 74.0 37.0- xevnoC a552 052 871 093 40.0 51 54.0 56.0- CEClioS mdclomm( –3) evacnoC b25 25 63 07 51.0 31 01.0 71.0 xevnoC 16 a 16 94 77 51.0 01 92.0 65.0- CECyalC mdclomm( –3) evacnoC 242 a 062 981 653 70.0 14 50.0 20.1 xevnoC 932 a 442 791 572 51.0 92 51.0 23.0- 1p = minimum level of significance in the Kolmogorov-Smirnov test; Min = minimum, Max = maximum, Var = variation; OM = organic matter; CEC = cation exchange capacity. Means followed by the same letter are not different (Student t test). Table 1 - Soil attributes of the concave and convex hillslope curvatures at 0.0-0.2-m and 0.2-0.4-m depth. Hillslope curvature and P adsorption in an Alfisol 823 Sci. Agric. (Piracicaba, Braz.), v.66, n.6, p.819-826, November/December 2009 etubirttA epolslliH erutavruc naeM naideM niM xaM p1 tneiciffeoC )%(.raV yrtemmysA sisotruK m2.0-0.0 aerakaeP mc()tK( 2) evacnoC 66.1 a 36.1 11.0 32.3 51.0 23 53.0 96.0 xevnoC a77.1 37.1 08.0 97.2 11.0 52 61.0 07.0- aerakaeP mc()bG( 2) evacnoC b11.0 21.0 10.0 82.0 51.0 35 94.0 405.0- xevnoC a42.0 32.0 20.0 44.0 51.0 63 41.0 23.0- tKHHW º( 2θ) evacnoC b65.0 75.0 33.0 38.0 80.0 81 94.0 60.0- xevnoC a46.0 56.0 83.0 78.0 51.0 61 90.0- 66.0- bGHHW º( 2θ) evacnoC a22.0 12.0 70.0 54.0 10.0 04 25.0 60.0 xevnoC a12.0 02.0 90.0 33.0 21.0 32 42.0 05.0- tKDCM )mn evacnoC a32 00.12 0.31 0.83 10.0 52 52.0 17.0- xevnoC b91 00.81 0.9 0.13 10.0 52 07.0 50.0 bGDCM )mn( evacnoC a171 402 9.91 664 10.0 36 97.0 301.0 xevnoC 091 a 402 00.6 704 10.0 74 43.0 10.0- )tK+bG(/bG evacnoC 60.0 b 60.0 61.0 00.0 01.0 4 56.0- 470.0- xevnoC 21.0 a 21.0 32.0 00.0 51.0 5 30.0- 31.0 1p = minimum significance level in the Kolmogorov-Smirnov test; Min = minimum; Max = maximum; Var = variation; WHH = width at half height; MCD = mean crystallite dimension; Gb/(Gb+Kt) =gibbsite ratio /(gibbsite+ kaolinite). Means followed by the same letter are not different (Student t test). Table 2 - Soil clay mineralogical attributes related to kaolinite (Kt) and gibbsite (Gb) on the concave and the convex hillslope curvatures at 0.2-0.4-m depth. OM content of both areas and the CEC of the con- vex 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 remain- ing properties fitted the spherical model. Several au- thors 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 con- cave 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 depen- dence. The remaining properties presented strong spa- tial 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 variabil- ity of chemical properties were found in the convex area (Table 3), demonstrating the influence of relief on the spatial behavior of these properties (Kravchenko & Bul- lock, 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 phos- phorus (Table 4) and the amount of available phospho- rus (Table 1), one can suggest that gibbsite is prima- rily responsible for higher P adsorption and lower P availability in the convex than in the concave area be- cause it has greater proportion and better crystalliza- tion in this area (Table 2). CONCLUSIONS Convex hillslope curvatures were more favorable environments for the formation of gibbsite than con- Barbieri et al.824 Sci. Agric. (Piracicaba, Braz.), v.66, n.6, p.819-826, November/December 2009 cave 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 spa- tial variability for all of the properties studied as com- pared to the concave area, indicating that the sam- pling strategy used in the two areas must be differ- ent. etubirttA epolslliH erutavruc ledoM )m(egnaR 1 oC 2C 0 C+ 1 3 )1C+oC(/oC % r2 m2.0-0.0 mdgm(P 3– ) evacnoC lacirehpS 101 75.11 94.32 94 9.0 xevnoC lacirehpS 62 10.0 4.11 1.0 8.0 mdg(MO 3– ) evacnoC lacirehpS 22 80.0 17.3 2 8.0 xevnoC lacirehpS 22 61.0 98.3 4 9.0 gkg(yalC 1– ) evacnoC lacirehpS 59 765 3252 22 9.0 xevnoC lacirehpS 32 72 4021 2 8.0 CEClioS mdclomm( 3– ) evacnoC lacirehpS 68 49.81 59.54 14 9.0 xevnoC laitnenopxE 52 8.4 56.63 31 7.0 CECyalC mdclomm( 3– ) evacnoC laitnenopxE 33 801 4741 7 9.0 xevnoC lacirehpS 02 83 2501 4 7.0 m4.0-2.0 mdgm(P 3– ) evacnoC lacirehpS 98 2.9 12.32 93 9.0 xevnoC lacirehpS 32 4.0 10.01 4 8.0 mdg(MO 3– ) evacnoC naissuaG 82 86.1 73.3 94 8.0 xevnoC naissuaG 91 5.0 76.3 31 9.0 gkg(yalC 1– ) evacnoC lacirehpS 58 2791 5493 05 9.0 xevnoC lacirehpS 82 14 3021 3 9.0 CEClioS mdclomm( 3– ) evacnoC laitnenopxE 08 9.51 11.16 62 9.0 xevnoC naissuaG 42 7.02 73.74 34 8.0 CECyalC mdclomm( 3– ) evacnoC laitnenopxE 53 501 2801 01 8.0 xevnoC lacirehpS 03 61 226 3 8.0 aerakaeP mc()tK( 2) evacnoC laitnenopxE 63 120.0 4041.0 51 6.0 xevnoC lacirehpS 81 8500.0 491.0 03 7.0 aerakaeP mc()bG( 2) evacnoC lacirehpS 81 410000.0 74200.0 2 8.0 xevnoC lacirehpS 71 51000.0 27600.0 2.0 7.0 2º(tKHHW θ) evacnoC lacirehpS 22 10000.0 9500.0 3.0 7.0 xevnoC lacirehpS 81 720000.0 44900.0 8 8.0 2º(bGHHW θ) evacnoC laitnenopxE 82 64000.0 075500.0 3 9.0 xevnoC lacirehpS 81 50000.0 29100.0 4 9.0 )mn(tKDCM evacnoC lacirehpS 91 10.0 42.62 30.0 7.0 xevnoC lacirehpS 51 49.0 45.91 5 7.0 )mn(bGDCM evacnoC laitnenopxE 92 0211 0429 21 5.0 xevnoC lacirehpS 61 0041 09946 2 5.0 tK+bG/bG evacnoC laitnenopxE 82 181000.0 3100.0 41 7.0 xevnoC lacirehpS 81 50000.0 97100.0 3 9.0 1C0 = nugget effect; 2C0+C1 = sill; 5 C0/(C0+C1 ) = degree of spatial dependence; OM = organic matter; CEC = cation exchange capacity; WHH = width at half height, MCD = mean crystallite dimension; Kt = kaolinite; Gb = gibbsite. Table 3 - Geostatistical attributes of the soil attributes at 0.0-0.2 m and 0.2-0.4 m depth. Hillslope curvature and P adsorption in an Alfisol 825 Sci. Agric. (Piracicaba, Braz.), v.66, n.6, p.819-826, November/December 2009 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. 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). 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. 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. 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. EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Sistema Brasileiro de Classificação de Solos. 2.ed. Rio de Janeiro: Embrapa, 2006. 412p. 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. 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. 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. 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. 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. 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. 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. Table 4 - Maximum sorption capacity of phosphorus in the two studied areas. Means followed by the same letter are not different (Student t test). etiS aerAepolslliHevacnoC aerAepolslliHxevnoC gkgm------------------ 1– ------------------ 1 534 555 2 625 885 3 005 885 4 534 625 5 454 625 6 614 454 -aeM n a164 b935 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. 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. ISAAKS, E.H.; SRIVASTAVA, R.M. An introduction to applied geoestatistics. New York: Oxford University Press, 1989. 561p. 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. 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. 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. LEAL, J.R.; VELOSO, A.C.X. Phosphate adsorption in latosols under Cerrado vegetation. Pesquisa Agropecuária Brasileira, v.8, p.81-88, 1973. 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. McBRATNEY, A.B.; WEBSTER, R. How many observations are needed for regional estimation of soil properties? Soil Science, v.135, p.177-183, 1983. 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. 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. 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. 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. 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. RAIJ, B. van. Análise química para avaliação da fertilidade de solos tropicais. Campinas, Instituto Agronômico, 2001. 285p. 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. 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. SCHLOTZHAVER, S.D.; LITTELL, R.C. SAS: system for elementary statistical analysis. 2 ed. Cary: SAS, Institute, 1997. 905p. 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. 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). Barbieri et al.826 Sci. Agric. (Piracicaba, Braz.), v.66, n.6, p.819-826, November/December 2009 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. TROEH, F.R. Landform equations fitted to contour maps. American Journal of Science, v.263, p.616-627, 1965. 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. 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. Received November 11, 2008 Accepted Agosto 27, 2009