Variability of apparently homogeneous soilscapes in São Paulo state, Brazil: II. quality of soil maps

The quality of semi-detailed (scale 1:100.000) soil maps and the utility of a taxonomically based legend were assessed by studying 33 apparently homogeneous fields with strongly weathered soils in two regions in São Paulo State: Araras and Assis. An independent data set of 395 auger sites was used to determine purity of soil mapping units and analysis of variance within and between mapping units and soil classification units. Twenty three soil profiles were studied in detail. The studied soil maps have a high purity for some legend criteria, such as B horizon type (> 90%) and soil texture class (> 80%). The purity for the "trophic character" (eutrophic, dystrophic, allic) was only 55% in Assis. It was 88% in Araras, where many soil units had been mapped as associations. In both regions, the base status of clay-textured soils was generally better than suggested by the maps. Analysis of variance showed that mapping was successful for "durable" soil characteristics such as clay content (> 80% of variance explained) and cation exchange capacity (≥ 50% of variance explained) of 0-20 and 60-80 cm layers. For soil characteristics that are easily modified by management, such as base saturation of the 0-20 cm layer, the maps had explained very little (< 15%) of the total variance in the study areas. Intermediate results were obtained for base saturation of the 60-80 cm layer (56% in Assis; 42% in Araras). Variance explained by taxonomic groupings that formed the basis for the legend of the soil maps was similar to, often even smaller than, variance explained by mapping units. The conclusion is that map boundaries have been very carefully located, but descriptions of mapping units could be improved. In future mappings, this could possibly be done at low cost by (a) bulk sampling to remove short range variation and enhance visualization of spatial patterns at distances > 100 m; (b) taking advantage of correlations between easily measured soil characteristics and chemical soil properties and, (c) unbending the link between legend criteria and a taxonomic system. The maps are well suited to obtain an impression of land suitability for high-input farming. Additional field work and data on former land use/management are necessary for the evaluation of chemical properties of surface horizons.

Variability of apparently homogeneous soilscapes in São Paulo state, Brazil: I. spatial analysis

The spatial variability of strongly weathered soils under sugarcane and soybean/wheat rotation was quantitatively assessed on 33 fields in two regions in São Paulo State, Brazil: Araras (15 fields with sugarcane) and Assis (11 fields with sugarcane and seven fields with soybean/wheat rotation). Statistical methods used were: nested analysis of variance (for 11 fields), semivariance analysis and analysis of variance within and between fields. Spatial levels from 50 m to several km were analyzed. Results are discussed with reference to a previously published study carried out in the surroundings of Passo Fundo (RS). Similar variability patterns were found for clay content, organic C content and cation exchange capacity. The fields studied are quite homogeneous with respect to these relatively stable soil characteristics. Spatial variability of other characteristics (resin extractable P, pH, base- and Al-saturation and also soil colour), varies with region and, or land use management. Soil management for sugarcane seems to have induced modifications to greater depths than for soybean/wheat rotation. Surface layers of soils under soybean/wheat present relatively little variation, apparently as a result of very intensive soil management. The major part of within-field variation occurs at short distances (< 50 m) in all study areas. Hence, little extra information would be gained by increasing sampling density from, say, 1/km² to 1/50 m². For many purposes, the soils in the study regions can be mapped with the same observation density, but residual variance will not be the same in all areas. Bulk sampling may help to reveal spatial patterns between 50 and 1.000 m.

Comparison of estimation methods of soil strength in five soils

In agriculture, the soil strength is used to describe the susceptibility to deformation by pressure caused by agricultural machine. The purpose of this study was to compare different methods for estimating the inherent soil strength and to identify their suitability for the evaluation of load support capacity, compaction susceptibility and root growth. The physical, chemical, mineralogical and intrinsic strength properties of seven soil samples, collected from five sampling pits at different locations in Brazil, were measured. Four clay (CS) and three sandy clay loam (SCL) soils were used. The clay soils were collected on a farm in Santo Ângelo, RS (28 º 16 ' 16 '' S; 54 º 13 ' 11 '' W 290 m); A and B horizons at the Universidade Federal de Lavras, Lavras, MG (21 º 13 ' 47 '' S; 44 º 58 ' 6'' W; 918 m) and on the farm Sygenta, in Uberlandia, MG (18 º 58 ' 37 '' S; 48 º 12 ' 05 '' W 866 m). The sandy clay loam soils were collected in Aracruz, ES (19 º 47 ' 10 '' S; 40 º 16 ' 29 '' W 81 m), and on the farm Xavier, Lavras, MG (21 º 13 ' 24 '' S; 45 º 05 ' 00 '' W; 844 m). Soil strength was estimated based on measurements of: (a) a pneumatic consolidometer, (b) manual pocket (non-rotating) penetrometer; and (c) automatic (rotating) penetrometer. The results of soil strength properties were similar by the three methods. The soil structure had a significant influence on soil strength. Results of measurements with both the manual pocket and the electric penetrometer were similar, emphasizing the influence of soil texture. The data showed that, to enhance the reliability of predictions of preconsolidation pressure by penetrometers, it is better to separate the soils into the different classes, rather than analyze them jointly. It can be concluded that the consolidometer method, although expensive, is the best when evaluations of load support capacity and compaction susceptibility of soil samples are desired.

Strategy of specification of management areas: rice grain yield as related to soil fertility

It is well-known nowadays that soil variability can influence crop yields. Therefore, to determine specific areas of soil management, we studied the Pearson and spatial correlations of rice grain yield with organic matter content and pH of an Oxisol (Typic Acrustox) under no- tillage, in the 2009/10 growing season, in Selvíria, State of Mato Grosso do Sul, in the Brazilian Cerrado (longitude 51º24' 21'' W, latitude 20º20' 56'' S). The upland rice cultivar IAC 202 was used as test plant. A geostatistical grid was installed for soil and plant data collection, with 120 sampling points in an area of 3.0 ha with a homogeneous slope of 0.055 m m-1. The properties rice grain yield and organic matter content, pH and potential acidity and aluminum content were analyzed in the 0-0.10 and 0.10-0.20 m soil layers. Spatially, two specific areas of agricultural land management were discriminated, differing in the value of organic matter and rice grain yield, respectively with fertilization at variable rates in the second zone, a substantial increase in agricultural productivity can be obtained. The organic matter content was confirmed as a good indicator of soil quality, when spatially correlated with rice grain yield.

Management zones in agriculture acording to the soil and landscape variables

The study of spatial variability of soil and plants attributes, or precision agriculture, a technique that aims the rational use of natural resources, is expanding commercially in Brazil. Nevertheless, there is a lack of mathematical analysis that supports the correlation of these independent variables and their interactions with the productivity, identifying scientific standards technologically applicable. The aim of this study was to identify patterns of soil variability according to the eleven physical and seven chemical indicators in an agricultural area. It was used two multivariate techniques: the hierarchical cluster analysis (HCA) and the principal component analysis (PCA). According to the HCA, the area was divided into five management zones: zone 1 with 2.87ha, zone 2 with 0.8ha, zone 3 with 1.84ha, zone 4 with 1.33ha and zone 5 with 2.76ha. By the PCA, it was identified the most important variables within each zone: V% for the zone 1, CTC in the zone 2, levels of H+Al in the zone 4 and sand content and altitude in the zone 5. The zone 3 was classified as an intermediate zone with characteristics of all others. According to the results it is concluded that it is possible to separate into groups (management zones) samples with the same patterns of variability by the multivariate statistical techniques.

Multiple Geotechnological Tools Applied to Digital Mapping of Tropical Soils

ABSTRACT In recent years, geotechnologies as remote and proximal sensing and attributes derived from digital terrain elevation models indicated to be very useful for the description of soil variability. However, these information sources are rarely used together. Therefore, a methodology for assessing and specialize soil classes using the information obtained from remote/proximal sensing, GIS and technical knowledge has been applied and evaluated. Two areas of study, in the State of São Paulo, Brazil, totaling approximately 28.000 ha were used for this work. First, in an area (area 1), conventional pedological mapping was done and from the soil classes found patterns were obtained with the following information: a) spectral information (forms of features and absorption intensity of spectral curves with 350 wavelengths -2,500 nm) of soil samples collected at specific points in the area (according to each soil type); b) obtaining equations for determining chemical and physical properties of the soil from the relationship between the results obtained in the laboratory by the conventional method, the levels of chemical and physical attributes with the spectral data; c) supervised classification of Landsat TM 5 images, in order to detect changes in the size of the soil particles (soil texture); d) relationship between classes relief soils and attributes. Subsequently, the obtained patterns were applied in area 2 obtain pedological classification of soils, but in GIS (ArcGIS). Finally, we developed a conventional pedological mapping in area 2 to which was compared with a digital map, ie the one obtained only with pre certain standards. The proposed methodology had a 79 % accuracy in the first categorical level of Soil Classification System, 60 % accuracy in the second category level and became less useful in the categorical level 3 (37 % accuracy).

Spatial variability of soil chemical properties after coffee tree removal

Assessing the spatial variability of soil chemical properties has become an important aspect of soil management strategies with a view to higher crop yields with minimal environmental degradation. This study was carried out at the Centro Experimental of the Instituto Agronomico, in Campinas, São Paulo, Brazil. The aim was to characterize the spatial variability of chemical properties of a Rhodic Hapludox on a recently bulldozer-cleaned area after over 30 years of coffee cultivation. Soil samples were collected in a 20 x 20 m grid with 36 sampling points across a 1 ha area in the layers 0.0-0.2 and 0.2-0.4 m to measure the following chemical properties: pH, organic matter, K+, P, Ca2+, Mg2+, potential acidity, NH4-N, and NO3-N. Descriptive statistics were applied to assess the central tendency and dispersion moments. Geostatistical methods were applied to evaluate and to model the spatial variability of variables by calculating semivariograms and kriging interpolation. Spatial dependence patterns defined by spherical model adjusted semivariograms were made for all cited soil properties. Moderate to strong degrees of spatial dependence were found between 31 and 60 m. It was still possible to map soil spatial variability properties in the layers 0-20 cm and 20-40 cm after plant removal with bulldozers.

Temporal variability of soil water content under different surface conditions in the semiarid region of the Pernambuco state

Rainfall in the semiarid region of Pernambuco is characterized by irregular distribution in time and space, which significantly hinders the rainfed agriculture in the region. This work aims to evaluate the temporal profile of soil moisture in the semiarid region of the Pernambuco State (Brazil) and the effect of different soil surface conditions on soil water content variation and the yield of rainfed beans. To monitor soil water content, five plots 4.5 m wide by 11 m long were installed in a Yellow Argisol (Ultisol). The following treatments were adopted in the experimental plots: natural vegetation, bean intercropped with cactus, beans planted down the slope, beans planted along contour lines with mulch and rock barriers, and bare soil. In each plot, eight PVC access tubes were installed for monitoring the soil water content profile at depths of 0.20 and 0.40 m using a neutron probe device. The surface condition significantly influenced the soil water content variation, both in the dry and rainy seasons. The use of mulch, associated with rock barriers, provided higher soil water content levels than the other treatments and increased the rainfed beans production.

Implications of the variability in soil penetration resistance for statistical analysis

Soil penetration resistance (PR) is a measure of soil compaction closely related to soil structure and plant growth. However, the variability in PR hampers the statistical analyses. This study aimed to evaluate the variability of soil PR on the efficiency of parametric and nonparametric analyses in indentifying significant effects of soil compaction and to classify the coefficient of variation of PR into low, medium, high and very high. On six dates, the PR of a typical dystrophic Red Ultisol under continuous no-tillage for 16 years was measured. Three tillage and/or traffic conditions were established with the application of: (i) no chiseling or additional traffic, (ii) additional compaction, and (iii) chiseling. On each date, the nineteen PR data (measured at every 1.5 cm to a depth of 28.5 cm) were grouped in layers with different thickness. In each layer, the treatment effects were evaluated by variance (ANOVA) and Kruskal-Wallis analyses in a completely randomized design, and the coefficients of variation of all analyses were classified (low, intermediate, high and very high). The ANOVA performed better in discriminating the compaction effects, but the rejection rate of null hypothesis decreased from 100 to 80 % when the coefficient of variation increased from 15 to 26 %. The values of 15 and 26 % were the thresholds separating the low/intermediate and the high/very high coefficient variation classes of PR in this Ultisol.

Uncertainties in the prediction of spatial variability of soil CO2 emissions and related properties

The soil CO2 emission has high spatial variability because it depends strongly on soil properties. The purpose of this study was to (i) characterize the spatial variability of soil respiration and related properties, (ii) evaluate the accuracy of results of the ordinary kriging method and sequential Gaussian simulation, and (iii) evaluate the uncertainty in predicting the spatial variability of soil CO2 emission and other properties using sequential Gaussian simulations. The study was conducted in a sugarcane area, using a regular sampling grid with 141 points, where soil CO2 emission, soil temperature, air-filled pore space, soil organic matter and soil bulk density were evaluated. All variables showed spatial dependence structure. The soil CO2 emission was positively correlated with organic matter (r = 0.25, p < 0.05) and air-filled pore space (r = 0.27, p < 0.01) and negatively with soil bulk density (r = -0.41, p < 0.01). However, when the estimated spatial values were considered, the air-filled pore space was the variable mainly responsible for the spatial characteristics of soil respiration, with a correlation of 0.26 (p < 0.01). For all variables, individual simulations represented the cumulative distribution functions and variograms better than ordinary kriging and E-type estimates. The greatest uncertainties in predicting soil CO2 emission were associated with areas with the highest estimated values, which produced estimates from 0.18 to 1.85 t CO2 ha-1, according to the different scenarios considered. The knowledge of the uncertainties generated by the different scenarios can be used in inventories of greenhouse gases, to provide conservative estimates of the potential emission of these gases.

Genesis of textural contrasts in subsurface soil horizons in the Northern Pantanal-Brazil

The Pantanal region can be characterized as a quaternary floodplain with predominant sedimentation in the form of alluvial fans. In the geomorphologic and sedimentary evolution, the avulsion process is inherent to this depositional system and its dynamics, together with surface water floods, influence soil sedimentation on this plain. The knowledge and differentiation of these two events can contribute to a better understanding of the variability of soil properties and distribution under the influence of these sedimentation processes. Therefore, this study investigated the genesis of soils in the Northern Pantanal with textural contrasts in deeper horizons and their relationship with the depositional system dynamics. We analyzed four soil profiles in the region of Barão de Melgaço, Mato Grosso State, Brazil (RPPN SESC Pantanal). Two profiles were sampled near the Rio Cuiabá (AP1 and AP4) and two near the Rio São Lourenço (AP10 and AP11). In AP11, the horizons contrast in particle size between the profile basis and the surface. In AP1, AP4 and AP10, the horizons overlaying the sand layer have similar particle size properties, mainly in terms of sand distribution. In the first case, floods (surface water) seem to have originated the horizons and layers with contrasting texture. In the second case, avulsion is the most pronounced process. Therefore, the two modes can form soils with contrasting texture that are discriminable by soil morphology, based on the distinct features associated to the specific sedimentation processes.

Water storage variability in a vineyard soil in the southern highlands of Santa Catarina state

In the subtropical regions of southern Brazil, rainfall distribution is uneven, which results in temporal variability of soil water storage. For grapes, water is generally available in excess and water deficiency occurs only occasionally. Furthermore, on the Southern Plateau of Santa Catarina, there are differences in soil properties, which results in high spatial variability. These two factors affect the composition of wine grapes. Spatio-temporal analyses are therefore useful in the selection of cultural practices as well as of adequate soils for vineyards. In this way, well-suited areas can produce grapes with a more appropriate composition for the production of quality wines. The aim of this study was to evaluate the spatio-temporal variability of water storage in a Cambisol during the growth cycle of a Cabernet Sauvignon vineyard and its relation to selected soil properties. The experimental area consisted of a commercial 8-year-old vineyard in São Joaquim, Santa Catarina, Brazil. A sampling grid with five rows and seven points per row, spaced 12 m apart, was outlined on an area of 3,456 m². Soil samples were collected with an auger at these points, 0.30 m away from the grapevines, in the 0.00-0.30 m layer, to determine gravimetric soil moisture. Measurements were taken once a week from December 2008 to April 2009, and every two weeks from December 2009 to March 2010. In December 2008, undisturbed soil samples were collected to determine bulk density, macro- and microporosity, and disturbed samples were used to quantify particle size distribution and organic carbon content. Results were subjected to descriptive analysis and semivariogram analysis, calculating the mean relative difference and the Pearson correlation. The average water storage in a Cambisol under grapevine on ridges had variable spatial dependence, i.e., the lower the average water storage, the higher the range of spatial dependence. Water storage had a stable spatial pattern during the trial period, indicating that the points with lower water storage or points with higher water storage during a certain period maintain these conditions throughout the experimental period. The relative difference is a simple method to identify positions that represent the average soil water storage more adequately at any time for a given area.

Variability of soil fertility properties in areas planted to sugarcane in the State of Goias, Brazil

Soil sampling should provide an accurate representation of a given area so that recommendations for amendments of soil acidity, fertilization and soil conservation may be drafted to increase yield and improve the use of inputs. The aim of this study was to evaluate the variability of soil fertility properties of Oxisols in areas planted to sugarcane in the State of Goias, Brazil. Two areas of approximately 8,100 m² each were selected, representing two fields of the Goiasa sugarcane mill in Goiatuba. The sugarcane crop had a row spacing of 1.5 m and subsamples were taken from 49 points in the row and 49 between the row with a Dutch auger at depths of 0.0-0.2 and 0.2-0.4 m, for a total of 196 subsamples for each area. The samples were individually subjected to chemical analyses of soil fertility (pH in CaCl2, potential acidity, organic matter, P, K, Ca and Mg) and particle size analysis. The number of subsamples required to compose a sample within the acceptable ranges of error of 5, 10, 20 and 40 % of each property were computed from the coefficients of variation and the Student t-value for 95 % confidence. The soil properties under analysis exhibited different variabilities: high (P and K), medium (potential acidity, Ca and Mg) and low (pH, organic matter and clay content). Most of the properties analyzed showed an error of less than 20 % for a group of 20 subsamples, except for P and K, which were capable of showing an error greater than 40 % around the mean. The extreme variability in phosphorus, particularly at the depth of 0.2-0.4 m, attributed to banded application of high rates of P fertilizers at planting, places limitations on assessment of its availability due to the high number of subsamples required for a composite sample.

SPATIAL VARIABILITY OF SOIL PROPERTIES IN AN AGRARIAN REFORM SETTLEMENT

ABSTRACT The study of soil chemical and physical properties variability is important for suitable management practices. The aim of this study was to evaluate the spatial variability of soil properties in the Malhada do Meio settlement to subsidize soil use planning. The settlement is located in Chapadinha, MA, Brazil, and has an area of 630.86 ha. The vegetation is seasonal submontane deciduous forest and steppe savanna. The geology is formed of sandstones and siltstones of theItapecuru Formation and by colluvial and alluvial deposits. The relief consists of hills with rounded and flat tops with an average altitude of 67 m, and frequently covered over by ferruginous duricrusts. A total of 183 georeferenced soil samples were collected at the depth of 0.00-0.20 m inPlintossolos, Neossolo andGleissolo. The following chemical variables were analyzed: pH(CaCl2), H+Al, Al, SB, V, CEC, P, K, OM, Ca, Mg, SiO2, Al2O3, and Fe2O3; along with particle size variables: clay, silt, and sand. Descriptive statistical and geostatistical analyses were carried out. The coefficient of variation (CV) was high for most of the variables, with the exception of pH with a low CV, and of sand with a medium CV. The models fitted to the experimental semivariograms of these variables were the exponential and the spherical. The range values were from 999 m to 3,690 m. For the variables pH(CaCl2), SB, and clay, there are three specific areas for land use planning. The central part of the area (zone III), where thePlintossolos Pétricos and Neossolos Flúvicos occur, is the most suitable for crops due to higher macronutrient content, organic matter and pH. Zones I and II are indicated for environmental preservation.

Temporal variability of soil CO2 emission after conventional and reduced tillage described by an exponential decay in time model

To study Assessing the impact of tillage practices on soil carbon losses dependents it is necessary to describe the temporal variability of soil CO2 emission after tillage. It has been argued that large amounts of CO2 emitted after tillage may serve as an indicator for longer-term changes in soil carbon stocks. Here we present a two-step function model based on soil temperature and soil moisture including an exponential decay in time component that is efficient in fitting intermediate-term emission after disk plow followed by a leveling harrow (conventional), and chisel plow coupled with a roller for clod breaking (reduced) tillage. Emission after reduced tillage was described using a non-linear estimator with determination coefficient (R²) as high as 0.98. Results indicate that when emission after tillage is addressed it is important to consider an exponential decay in time in order to predict the impact of tillage in short-term emissions.