Clay mineralogy and magnetic susceptibility of Oxisols in geomorphic surfaces

Studies analyzing the variability of clay minerals and magnetic susceptibility provide data for the delineation of site-specific management areas since many of their attributes are important to agronomy and the environment. This study aimed to evaluate the spatial variability of clay minerals, magnetic susceptibility, adsorbed phosphorus and physical attributes in Oxisols of sandstones in different geomorphic surfaces. For that purpose, soil samples were collected every 25 m along a transect located within the area where the geomorphic surfaces were identified and mapped. The transect occupied the central portion of 500 ha, where it was also sampled for density purposes with one sample per six hectares. Soil samples were collected at a depth of 0.0-0.2 m. The results of the physical, chemical, mineralogical and magnetic susceptibility analyses were subjected to statistical and geostatistical analyses. The nature of the clay minerals and magnetic susceptibility was dependent on the variation of the soil parent material. High values of magnetic susceptibility were associated with the presence of maghemite and magnetite of coarse size. The spatial variability of crystallinity and the content of Fe oxides, as well as magnetic susceptibility, were dependent on the age of the geomorphic surfaces. The youngest surface had greater spatial variability of these attributes. The iron (goethite and hematite) and aluminum (gibbsite) oxides in the youngest geomorphic surface influenced the low values of soil density and high values of total pore volume, micropores and P adsorption. The characterization of the spatial variability of Fe oxides and susceptibility allowed for the delineation of homogeneous areas.


Introduction
Humid tropical regions have areas of global concern about food production, even though there are extremely competitive areas for grain and fuel such as in the Midwest and Southeast of Brazil.Due to soil diversity, soil characterization is necessary, based on a system that takes into account the soil genesis.In this study, we characterized the clay mineralogy of highly weathered soils and their relationship with the landscape feature, aiming to present the potential of this tool to help delineate mapping units.
Tropical soils in humid areas and old geomorphic surfaces have highly weathered clay minerals with a predominance of 1:1 clay minerals, Fe, Al oxides and hydroxides (denominated in this study as oxides).Furthermore, these soils, also known as laterite, have an unquestionable importance to agronomy and civil engineering, because these soils respond to cement stabilization and, in some cases, lime stabilization (Townsend, 1985).Associations between goethite and hematite are the most common Fe oxides in weathered soils from tropical and subtropical regions (Schwertmann and Taylor, 1989).Several relationships between Fe and Al oxides and physical and chemical properties in the soil are reported in the literature (McNeal et al., 1968;Curi and Franzmeier, 1984;Schwertmann and Kämpf, 1985;Borggaard et al., 2004;Mesquita Filho and Torrent, 1993;Schaefer et al., 2004;Camargo et al., 2013).Camargo et al. (2008a) found that the spatial dependence of crystallinity of Fe and Al oxides and kaolinite is relied to landscape compartments in an Oxisol originated from basalt, enabling the delineation of areas based on soil-landscape models in detailed scales.The influence of Fe oxides and their crystallinity in aggregates (Camargo et al., 2008b) as well as in P spatial availability (Camargo et al., 2013) within the compartments was also observed, showing that the delimitation of Fe oxides allow for specific soil management.
Further studies would help in optimizing and recognizing the degree of homogeneity of mapping units, based on the characterization of Fe oxides and magnetic susceptibility (MS).Therefore, this study aimed to evaluate the spatial variability of clay minerals, MS, P adsorbed and physical attributes in Oxisols of sandstones of the Bauru Group in different geomorphic surfaces.

Sampling system
The area studied is located in the northeast of São Paulo State, Brazil, in the municipality of Guariba (21°24'17" S; 48º09'46" W; average altitude 640 m.The climate is mesothermal with dry winters (Cwa) (Köppen classification) with average precipitation of 1,400 mm, concentrated between Nov and Feb.The natural vegetation consists of tropical semideciduous forest and the soils are cropped with sugarcane plantations with a mechanical harvesting system.
The geomorphic surfaces (GS) were identified and mapped using topographic and stratigraphic criteria, based on visual observations of unevenness or breaks in the relief shape, as recommended by Ruhe (1969) and Daniels et al. (1971).Three GS's were identified (I, II and III).GS I is approximately 400-m long with a slope ranging from 0 to 4 %.It is considered the oldest surface of the study area with a depositional characteristic.The boundary between GS II and III is approximately at 1,500 m at the top.GS II and III are considered erosional and GS II has the shape of a gentle slope plan.The GS III has a characteristically steeper slope (7 %).Older GS are more stable, usually found in water dividers and, are, thus, depositional environments, since they have gentle slopes.Younger GS, on the other hand, are more variable and occurs in erosional environments (Graham et al., 1990;Phillips, 2004).
To evaluate the spatial variability of the attributes studied, soil samples were collected along a transect where the GS were identified.The direction of the transect was selected from the high topography following the lower elevations.Soil samples were collected every 25 m (at a depth of 0.0 to 0.2 m).We collected 40 samples near the transect in each GS to support the spatial variability analysis.The strip of land which is the transect as well as the points of sample collection of geomorphic surfaces near the transect measures 145 ha (Figure 1).
We delimited and georeferenced an area of 500 ha where the transect occupied a central position.A soil sample was collected every 6 ha at a depth of 0.0 to 0.2 m, totaling 88 samples.These samples were added to those collected from the transect and its side (Figure 1).The soil was classified (Soil Survey Staff, 1999) as Typic Hapludox in GS I, Typic Eutrudox and Typic Hapludox in GS II and Typic Eutrudox in GS III (Table 1).The parent material in the study area is Cretaceous sandstones of Adamantina Formation.

Physical analyses
The particle size analysis was performed using a 0.1 M NaOH solution as a dispersant in a low-speed stirring apparatus (Day, 1965).The clay was determined using a pipette, the fine and coarse sands were separated by sieving and silt was calculated by difference.Undisturbed samples were collected to determine soil porosity (porosity, macropores and micropores).The samples were saturated for 48 h in a tray with water up to 2/3 of the ring height.After saturation, the samples were drained to the potential equivalent to -0.006 MPa, using a tension table and bulk density was determined by the known volume method (Grossman and Reinsch, 2002)).These analyses were performed on samples collected in the 145 ha area.

Mineralogical analyses
Clay mineralogy was determined by x-ray diffraction (XRD).The analyses were conducted in non oriented slides.The clay was separated by centrifugation according to the Jackson method (1985).The characterization of hematite (Hm) and goethite (Gt) was performed after treating the clay fraction with 5 M NaOH (1g clay 100 mL -1 solution), according to the Norrish and Taylor method (1961) with modification by Kämpf and Schwertmann (1982).The ratio Gt/ (Gt+Hm) was obtained after calculating the reflection areas of Hm (012) and Gt (110).The peak area of Gt (110) was multiplied by 0.35, due to the intensity of 35 % of Hm (012) (Kämpf and Schwertmann, 1998).The mean crystallite dimension (MCD) was calculated from the width at half height (WHH) reflection minerals Hm (110 and 012) and Gt (110 and 111) using the Scherrer equation (Klug and Alexander, 1974).
For kaolinite (Kt) and gibbsite (Gb) characterization, Fe oxides were removed from the clay fraction using the dithionite-citrate-bicarbonate method (DCB) according to Mehra and Jackson (1960), and screened through a 0.10-mm mesh.The ratio Kt/(Kt + Gb) was calculated using the reflection areas of the Kt (001) and Gb (002) of the XRD patterns.The MCD of Gb and Kt was calculated from the WHH mineral reflexes Kt (001) and Gb (002) using the Scherrer equation (Klug and Alexander, 1974).
The diffractometer used was the Mini-Flex Rigaku-II, using a cathode copper with a nickel filter and ka radiation (20mA, 30kV).The scan speed used was 1 º2q/minute amplitude 23-49 °2q for the Hm and Gt characterization, and from 11 to 19 °2q for the Kt and Gb characterization.
To calculate the Gt content, the crystalline Fe content was multiplied by the ratio Gt/(Gt + Hm) and by 1.59.For the Hm content, the crystalline Fe content was   multiplied by 1.43 after the Fe content corresponding to Gt had been subtracted from this value, according to the equations:

Magnetic Susceptibility (MS)
MS was measured in a diffuse reflectance spectrometer coupled to an MS2B laboratory sensor of two frequencies.This sensor allowed the MS to be measured using a low and alternating magnetic field at an intensity of 80 A m -1 (10 -4 T) of varying frequency.

Phosphorus adsorption
P adsorption was assessed in accordance with the methodology of Fassbender and Igue (1967) using samples consisting of 2.0 g of air dried soil transferred to 50 ml centrifuge tube.20 mL of a 0.005 M KH 2 PO 4 solution containing CaCl 2 were added at a P concentration of 100 mg L -1 .The soil:solution ratio was 1:10.Samples containing phosphate solutions were stirred for 16 h at constant temperature.After stirring, centrifugation was performed at 2,500 rpm for 15 min.The supernatant was collected to determine the P concentration in equilibrium by means of a molecular absorption spectrophotometer.To determine P in the spectrophotometer, we transferred the supernatant to 1 ml aliquot in the test tube, adding 50 mL of a sulfomolybdic acid solution and three similar solutions of 200 g L -1 of ascorbic acid.This solution was stirred individually for 30 s and left to rest for 15 min.The spectrophotometer wavelength was selected for 650 nm, which corresponds to the maximum absorption peak of the complex formed.The amount of P adsorbed (P ads ) by the soil was calculated by subtracting the value determined in solution from the total amount added.This analysis was performed in the samples collected from the 145 ha area.

Treatment of selective dissolution of maghemite
The treatment of maghemite selective dissolution was applied in accordance with the methodology proposed by Schwertmann and Fechter (1984), which consists of treating 1 g air-dried soil with 40 ml of 1.8 M H 2 SO 4 75-80 °C for 2 h and eliminating part of the Hm and/or Gt existing in the samples.Therefore, we carried out a test with two samples of each GS changing the exposure time (from 45 min to 2 h) and the H 2 SO 4 concentration (0.8 to 1.8 M) to determine the most suitable conditions.The results of the test showed that the best time was 1.2 h and best concentration of H 2 SO 4 was 1.3 M (unpublished data).The Fe content extracted by H 2 SO 4 was determined colorimetrically using a spectrophotometer following the method described by Olson and Ellis (1982).The waste dissolution was subjected to MS analysis after washing with deionized water and dried.

Statistical and geostatistical analyses
The variability of soil attributes was evaluated by exploratory data analysis, calculating the mean, median, coefficient of variation and the coefficient of skewness and kurtosis.The mean difference between the GS was tested by Tukey test.The analysis of spatial dependence of data was performed by geostatistics (Matheron, 1963;Isaaks and Srivastava, 1989;Vieira et al., 1983) using semivariograms.To determine spatial dependence, we used GS + software, and to model and set the semivariogram, the models considered were spherical, exponential and gaussian (McBratney and Webster, 1983).

Results and Discussion
Higher clay and silt contents and smaller sand contents were found in the GS III area (Table 2).This can be explained by a variation of the sandstone strata, sandy on top, with increasing clay content at the lower elevations.The sandstones of Adamantina Formation (Bauru Group) have fine granulation, usually with cross bedding, alternating with clay-rich materials such as mudstones.In the site studied, the higher clay content and the higher Fe 2 O 3 levels (Table 1), indicate that the GS III is a region under the influence of mudstones, which are characterized by a higher concentration of clay.Soils affected by mudstones have peculiarities that affect their physical and chemical attributes and, therefore, should be studied in order to implement the appropriate management procedures.
The crystallinity of Gt (111) was different for GS III from the other surfaces (p < 0.05).In GS III, Gt (111) had the highest degree of crystallinity, corroborated by high MCD values.The opposite occurred with the Hm mineral, which was different from WHH (with the exception of Hm ( 012)) and from MCD between GS I and the other surfaces (p < 0.05).The Hm presented a lower degree of crystallinity in GS I and II (Table 3).
The AS in Gt and Hm had the highest values where the crystals showed a lower degree of crystallinity (Table 4).This is attributed to the isomorphic substitution of Fe by Al in the crystal structure, causing a contraction of the unit cell, thus decreasing the size of the crystals in these oxides, lowering the degree of crystallinity Norrish and Taylor, 1961).
The Gt and Hm contents were higher (p < 0.05) in GS III in relation to the other surfaces.In all GSs, Hm contents are greater than Gt contents; however, there was a greater ratio Gt/(Gt+Hm) in GS I and II (Table 4) though the difference in GS III was less pronounced.This result was consistent with the more yellowish color of the soil in GS I and II (Table 1).Higher levels of Fed and Feo were found in GS III, which is consistent with the higher Fe 2 O 3 values in this surface, thus confirming the variability in the parent material.Soils at the landscape summit tend to have less water content and higher temperatures (Curi and Franzmeier, 1984), favoring the dehydration step of the ferrihydrite to Hm (Torrent et al., 1982;Kampf and Schwertmann, 1983).In GS I, the pH was around 7.4 at the depth sampled for the mineral study (Table 1).A pH near neutrality favors ferrihydrite aggregation and, thereby, the Hm formation (Schwertmann and Murad, 1983).This is consistent with a higher degree of crystallinity of Hm in soils that occur in GS I, which is older and more stable.
In GS III, the source material may provide higher Fe content for Hm formation but with a lower degree of crystallinity.However, in this GS, conditions were favorable for increasing the degree of crystallinity of Gt.Organic matter (OM) and Gt formation over Hm (Schwertmann and Taylor, 1989).The role of OM consists of complexing Fe, inhibiting the formation of ferrihydrite, which is an obligatory Hm precursor (Curi and Franzmeier, 1984).However, in this study, the influence of organic matter may be affecting the degree of crystallinity of the mineral, not its content.
In addition to geochemical processes, the position of soils in the landscape affects the crystallinity of Fe oxides (Curi and Franzmeier, 1984).These authors found Gt with a lower degree of crystallinity in the Standard deviation; (2) Coefficient of variation (%); (3) g kg -1 .GS= geomorphic surface; Var.= Variance; Max.= Maximum; Min.= Minimum; TS = Total Sand.Values in a column followed by the same letter were not different (Tukey test, p < 0.05).Standard deviation; (2) Coefficient of variation (%).GS = geomorphic surface; Min.= Minimum; Max.= Maximum; Var.= Variance; Hm = Hematite, Gt = Goethite.Values in a column followed by the same letter were not different (Tukey test, p < 0.05).
lower landscape, which disagrees with the current study.However, the source material in the study site used by Curi and Franzmeier (1984) was the product of changes in the basaltic rocks.These results indicate that the variation of the source material indicates that generalizations about the formation of Fe oxides in the soil-landscape relationships should be approached with caution.
The MS was also higher in GS III, which is consistent with the higher clay and Fe oxide contents in the soil surface.The mineralogical origin of MS in soils with low levels of ferrimagnetic minerals cannot be identified by using only the XRD technique.The magnetic behavior is more evident in soils with greater clay fraction (Fontes et al., 1996).In the sand fraction, magnetite is oxidized directly in Hm, while in the clay fraction magnetite is oxidized in maghemite (Singer et al., 1995).The soils used in this study were sandstone and the MS values were, on average, relatively lower than in maghemite soils (2.86 to 4.40 × 10 -8 m 3 kg -1 ) (Thompson and Oldfield, 1986).However, the highest MS value (2.91 × 10 -8 m 3 kg -1 ) in GS III shows the need for further studies to investigate the presence of ferrimagnetic minerals in these soils.
At 1.2 h of exposure to H 2 SO 4 (1.3 M), the MS in the residues decreased drastically and remained constant at low values, except for GS III (Figure 2A), but the Fe content extracted (Figure 2B) continued to increase in time.These results confirm that MS was affected by the presence of maghemite that had dissolved in the soil after 1.2 h, because the iron content extracted increased thereafter due to the dissolution of other Fe oxides such as Hm and Gt.Costa et al. (1999) found that for soils with a high content of Fe oxides (Fe 2 O 3 content higher than 180 g kg -1 ) the period of exposure to H 2 SO 4 (1.8 M) for maghemite dissolution was 2 h.The difference in relation to this study may be attributed to intermediate content of Fe oxides in the soil compared to the works cited and consequently with intermediate maghemite contents.These results indicate that the period of exposure to H 2 SO 4 should take into account the Fe oxide content in each soil.
The MS measures in GS III remained constant at high values after successive Fe extractions with H 2 SO 4 over several exposure periods (Figure 2A).This result can be attributed to the presence of magnetite from the rock with very small size in addition to goethite, hematite and maghemite, since soils originating from this parent material may contain its traces.Moreover, magnetite of coarse size is not dissolved by H 2 SO 4 , which may have affected the measurements.Standard deviation; (2) Coefficient of variation (%); (3) g kg -1 ; (4) 10 -8 m 3 kg -1 .GS= geomorphic surface; Min.= Minimum; Max.= Maximum; Var.= Variance; Hm = Hematite; Gt = Goethite; Fed = iron oxides extracted by dithionite-citrate-bicarbonate; Feo = iron extracted by ammonium oxalate; Hm = Hematite, Gt = Goethite.Values in a column followed by the same letter were not different (Tukey test, p < 0.05).Kt and Gb showed opposite behaviors.Kt had a lower degree of crystallinity in surface III, while Gb had a greater degree of crystallinity in this area (Table 5).Regarding the Kt/(Kt + Gb) ratio, we observed a difference between all GS and lower values, that is , the lowest Kt rate over Gb was found in GS III (p < 0.05).However, in this study, a higher degree of crystallinity and a greater Kt rate with respect to Gb in older surfaces are related to higher sand contents in these surfaces with a consequent greater release of silica.The quartz of fine sand fraction may act as a continuous source of silica (Lepsch et al., 1977;Furian et al., 2002) keeping the high silicon activity inside favoring the formation of Kt (Lindsay, 1979) in relation to Gb.Nevertheless, in GS III, in addition to the smaller sand contents, the highest Fe 2 O 3 levels favor Gb formation, since Fe contributes to prevent Kt crystallization favoring Gb.

Sci
At larger coefficients of variation (CV) of the Gt/ (Gt+Hm) ratio contents of Fe oxides and MS were found in GS III (Table 4).We also found the presence of two soil classes in the same GS (Figure 1 and Table 1), indicating greater variability of these attributes in this GS.This is because the GS III is in a less stable landscape, corroborating studies carried out by Daniels et al. (1970), Parsons et al. (1970), Gamble et al. (1970) and Lepsch et al. (1977), who stated that soils in younger GS have greater variation in their attributes.For the attributes related to Fe oxides, skewness and kurtosis parameters were above 1.0; however, the mean and median of these attributes were similar without compromising the geostatistical analysis.The MS had different mean and median values and this attribute showed impairment to carry out the geostatistical analysis.
The temporal sequence of GSs was not related to the nature of the mineral fraction, but it influenced the degree of crystallinity and variability, mainly in Fe oxides.The contents of Fe oxides and MS in the soil were related to lithological variations in the area.These results emphasize the importance of knowing the stratigraphy of the ground in order to understand the relationship between cause and effect when applying the concepts and to understand GS soils as natural bodies.Standard deviation; (2) Coefficient of variation (%);GS = geomorphic surface; Min.= Minimum; Max.= Maximum; Var.= Variance; WHH = Width at half height (º2q); MCD = Mean Crystallite Dimension (nm).Values in a column followed by the same letter were not different (Tukey test, p < 0.05).The lowest Ds values and the highest values of total porosity and microporosity were found in GS III because of the clay mineralogy of this GS (Table 6).Fe and Al oxides contribute to soil aggregation, acting as cementing agents (Lutz, 1936;Kroth and Page, 1947;Schwertmann and Kampf, 1985;Camargo et al., 2008b), increasing soil porosity and decreasing Ds.In GS III, we also observed the highest values of P ads when compared with the other GS.The highest values of P ads in GS III were related to higher Gt and Hm contents in the soil surface.Fe oxides adsorb phosphate ions in adsorption sites on the surface of iron oxides (Parfitt, 1978;Stumm, 1992;Persson et al., 1996).These results indicate the potential of clay mineralogy characterization for the understanding of physical and chemical properties in the soil.Thus, GSs that delimit the occurrence and crystallinity of clay minerals can be used as specific management areas, once they also delimit the behavior of soil physical and chemical properties.

Sci
The degree of crystallinity of clay minerals, the content of Fe oxides, AS in Gt and MS were spatially dependent (Table 7).The spherical model was fitted to the semivariogram of all attributes that showed spatial dependence, except for the Gt content , the Gt /(Hm + Gt) ratio, Fed, Feo and MS that showed the gaussian model fitted to the semivariograms.For attributes that showed spatial dependence at the site of 145 ha, we constructed spatial distribution maps (Figures 3,4 and 5).The maps of the DMC Gt (111 and 110) (Figures 3A  and 3D) had similar patterns, where a sharp transition occured between surfaces II and III.
The Hm showed a higher degree of crystallinity (higher MCD) in GS I and the limit between GS I and GS II corresponded to the limit found in the field (Figure 2C).The maps for the contents of Gt, Hm, Fed, Feo and MS showed similar spatial distributions, lower values for GSs I and II and higher values in GS III (Figures 4A, 4B,  4D, 4E, 4F).Kt had the highest degree of crystallinity in GSs I and II, and the Gb showed the highest degree of crystallinity in GS III.These results corroborate the mean difference for both minerals (Figures 5A and 5B).
The maps confirmed the age of GS III (youngest) and its unstable position in the landscape, mainly with regard to the contents of Fe oxides (Gt, Hm, Fed and Feo) and SM.In this GS, we found a greater spatial variability of attributes, agreeing with the variability found in simple statistics (Table 4).In GS III, more than one class of values was found for these attributes whereas in GSs I and II, only one class was found.
To evaluate the effect of scale variability of clay mineralogy, statistical and geostatistical analyses were performed over an area of 500 ha, where the transect occupied the central position as we analyzed the contents of Gt and Hm, ratio Gt/(Gt + Hm), Fed, Feo and MS (Table 8 and Figure 5).The same attributes that presented spatial dependence in the 145 ha area also presented spatial dependence in the 500 ha area.
A correspondence of the spatial distribution pattern between contents of Gt Hm, Fed, Feo and SM and boundaries of GS (Figure 6) was also observed in the 500-ha area.GSs conditioned the variability of Fe oxides on a large scale, indicating that these mappings can be used in homogeneous areas, since these oxides have many attribute covariates of agronomic interest.Standard deviation; (2) Coefficient of variation (%); (3) g kg -1 .Var. = Variance; Min.= Minimum; Max.= Maximum; Hm = Hematite; Gt = Goethite; Fed = iron oxides extracted by dithionite-citrate-bicarbonate; Feo = iron extracted by ammonium oxalate; MS = magnetic susceptibility (10 -8 m 3 kg -1 ); Hm = Hematite, Gt = Goethite.
More than one class of values for the contents of Gt, Hm, Fed, Feo and SM (Figure 6) was found within a soil class (Figure 1, Table 1), mainly in GS III.This shows variability within mapping units that may be relevant in delimiting site-specific management areas for agricultural purposes expressed by the influence of Fe oxides in soil attributes.
The techniques used in this study demonstrated the bilateralism of the GS and the influence of these Fe oxides and MS, using univariate statistical and geostatistical analyses.They also allowed for the identification of variability within the minimum mapping unit that may help in the establishment of site-specific management areas on a detailed scale even in large areas.The level of characterization proposed here can be used for areas with Fe contents and the presence of geomorphic surfaces similar to those studied, that recognized features in large tropical areas that control

Conclusions
1.The nature of the clay minerals and magnetic susceptibility were dependent on the variation of the soil parent material.
2. The high values of magnetic susceptibility were associated with the presence of maghemite and magnetite of coarse size.3. The spatial variability of crystallinity of the content of both Fe oxides, as well as the magnetic susceptibility, were dependent on the age of the geomorphic surfaces.Younger surfaces (III) showed greater spatial variability of these attributes.4. The iron (goethite and hematite) and aluminum

Figure 2 -
Figure 2 -A) Magnetic susceptibility of the sample after Fe extraction by H 2 SO 4 (1.3 M) regarding exposure period and B) Concentration of Fe extracted with H 2 SO 4 (1.3 M) in terms of exposure time.
The use of MS, as noted, enables this characterization to replace costly and time-consuming chemical and mineralogical soil analysis.

Table 1 -
Chemical and physical characteristics of the profiles by geomorphic surface.

Table 2 -
Descriptive statistics of the texture attributes in the area of 145 ha.

Table 3 -
Descriptive statistics of Width at Half Height [WHH (º2q)] and of Mean Crystallite Dimension [MCD (nm)] of Fe oxides, area of 145 ha at 0.0 to 0.2 m.

Table 4 -
Descriptive statistics of Al substitution (AS), contents of Fe oxides and magnetic susceptibility (MS) of the area of 145 ha, at 0.0 to 0.2 m.

Table 5 -
Descriptive statistics of kaolinite (Kt) and gibbsite (Gb) of area of 145 ha, at 0.0 to 0.2 m.

Table 6 -
Mean values of density (Ds), soil porosity and P adsorbed (Pads) at 0.00 to 0.20 m.

Table 7 -
Parameters of the models fitted to the semivariograms of attributes studied in the area of 145 ha.

Table 8 -
Descriptive statistics of the attributes the study area (area 500 ha) at 0.0 to 0.2 m.