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.

Spatial and temporal variability of crop yield and some Rhodic Hapludox properties under no-tillage

Soil properties are closely related with crop production and spite of the measures implemented, spatial variation has been repeatedly observed and described. Identifying and describing spatial variations of soil properties and their effects on crop yield can be a powerful decision-making tool in specific land management systems. The objective of this research was to characterize the spatial and temporal variations in crop yield and chemical and physical properties of a Rhodic Hapludox soil under no-tillage. The studied area of 3.42 ha had been cultivated since 1985 under no-tillage crop rotation in summer and winter. Yield and soil property were sampled in a regular 10 x 10 m grid, with 302 sample points. Yields of several crops were analyzed (soybean, maize, triticale, hyacinth bean and castor bean) as well as soil chemical (pH, Soil Organic Matter (SOM), P, Ca2+, Mg2+, H + Al, B, Fe, Mn, Zn, CEC, sum of bases (SB), and base saturation (V %)) and soil physical properties (saturated hydraulic conductivity, texture, density, total porosity, and mechanical penetration resistance). Data were analyzed using geostatistical analysis procedures and maps based on interpolation by kriging. Great variation in crop yields was observed in the years evaluated. The yield values in the Northern region of the study area were high in some years. Crop yields and some physical and soil chemical properties were spatially correlated.

Spatial variability of physical attributes of an alfisol under different hillslope curvatures

The influence of relief forms has been studied by several authors and explains the variability in the soil attributes of a landscape. Soil physical attributes depend on relief forms, and their assessment is important in mechanized agricultural systems, such as of sugarcane. This study aimed to characterize the spatial variability in the physical soil attributes and their relationship to the hillslope curvatures in an Alfisol developed from sandstone and growing sugarcane. Grids of 100 x 100 m were delimited in a convex and a concave area. The grids had a regular spacing of 10 x 10 m, and the crossing points of this spacing determined a total of 121 georeferenced sampling points. Samples were collected to determine the physical attributes related to soil aggregates, porosity, bulk density, resistance to penetration and moisture within the 0-0.2 and 0.2-0.4 m depth. Statistical analyses, geostatistics and Student's t-tests were performed with the means of the areas. All attributes, except aggregates > 2 mm in the 0-0.2 m depth and macroporosity at both depths, showed significant differences between the hillslope curvatures. The convex area showed the highest values of the mean weighted diameter, mean geometric diameter, aggregates > 2 mm, 1-2 mm aggregates, total porosity and moisture and lower values of bulk density and resistance to penetration in both depth compared to the concave area. The number of soil attributes with greater spatial variability was higher in the concave area.

Pedotransfer functions related to spatial variability of water retention attributes for lowland soils

The estimation of non available soil variables through the knowledge of other related measured variables can be achieved through pedotransfer functions (PTF) mainly saving time and reducing cost. Great differences among soils, however, can yield non desirable results when applying this method. This study discusses the application of developed PTFs by several authors using a variety of soils of different characteristics, to evaluate soil water contents of two Brazilian lowland soils. Comparisons are made between PTF evaluated data and field measured data, using statistical and geostatistical tools, like mean error, root mean square error, semivariogram, cross-validation, and regression coefficient. The eight tested PTFs to evaluate gravimetric soil water contents (Ug) at the tensions of 33 kPa and 1,500 kPa presented a tendency to overestimate Ug 33 kPa and underestimate Ug1,500 kPa. The PTFs were ranked according to their performance and also with respect to their potential in describing the structure of the spatial variability of the set of measured values. Although none of the PTFs have changed the distribution pattern of the data, all resulted in mean and variance statistically different from those observed for all measured values. The PTFs that presented the best predictive values of Ug33 kPa and Ug1,500 kPa were not the same that had the best performance to reproduce the structure of spatial variability of these variables.

Spatial patterns of preconsolidation pressure and soil moisture along transects in two directions under coffee

Information on the spatial structure of soil physical and structural properties is needed to evaluate the soil quality. The purpose of this study was to investigate the spatial behavior of preconsolidation pressure and soil moisture in six transects, three selected along and three across coffee rows, at three different sites under different tillage management systems. The study was carried out on a farm, in Patrocinio, state of Minas Gerais, in the Southeast of Brazil (18 º 59 ' 15 '' S; 46 º 56 ' 47 '' W; 934 m asl). The soil type is a typic dystrophic Red Latosol (Acrustox) and consists of 780 g kg-1 clay; 110 g kg-1 silt and 110 g kg-1 sand, with an average slope of 3 %. Undisturbed soil cores were sampled at a depth of 0.10-0.13 m, at three different points within the coffee plantation: (a) from under the wheel track, where equipment used in farm operations passes; (b) in - between tracks and (c) under the coffee canopy. Six linear transects were established in the experimental area: three transects along and three across the coffee rows. This way, 161 samples were collected in the transect across the coffee rows, from the three locations, while 117 samples were collected in the direction along the row. The shortest sampling distance in the transect across the row was 4 m, and 0.5 m for the transect along the row. No clear patterns of the preconsolidation pressure values were observed in the 200 m transect. The results of the semivariograms for both variables indicated a high nugget value and short range for the studied parameters of all transects. A cyclic pattern of the parameters was observed for the across-rows transect. An inverse relationship between preconsolidation pressure and soil moisture was clearly observed in the samples from under the track, in both directions.

Local influence for spatial analysis of soil physical properties and soybean yield using student's t-distribution

The modeling and estimation of the parameters that define the spatial dependence structure of a regionalized variable by geostatistical methods are fundamental, since these parameters, underlying the kriging of unsampled points, allow the construction of thematic maps. One or more atypical observations in the sample data can affect the estimation of these parameters. Thus, the assessment of the combined influence of these observations by the analysis of Local Influence is essential. The purpose of this paper was to propose local influence analysis methods for the regionalized variable, given that it has n-variate Student's t-distribution, and compare it with the analysis of local influence when the same regionalized variable has n-variate normal distribution. These local influence analysis methods were applied to soil physical properties and soybean yield data of an experiment carried out in a 56.68 ha commercial field in western Paraná, Brazil. Results showed that influential values are efficiently determined with n-variate Student's t-distribution.

Spatial relationships between soil attributes and corn yield in no-tillage system

Soil properties play an important role in spatial variability of crop yield. However, a low spatial correlation has generally been observed between maps of crop yield and of soil properties. The objectives of the present investigation were to assess the spatial pattern variability of soil properties and of corn yield at the same sampling intensity, and evaluate its cause-and-effect relationships. The experimental site was structured in a grid of 100 referenced points, spaced at 10 m intervals along four parallel 250 m long rows spaced 4.5 m apart. Thus, points formed a rectangle containing four columns and 25 rows. Therefore, each sampling cell encompassed an area of 45 m² and consisted of five 10 m long crop rows, in which the referenced points represented the center. Samples were taken from the layers 0-0.1 m and 0.1-0.2 m. Soil physical and chemical properties were evaluated. Statistical analyses consisted of data description and geostatistics. The spatial dependence of corn yield and soil properties was confirmed. The hypothesis of this study was confirmed, i.e., when sampling the soil to determine the values of soil characteristics at similar to sampling intensity as for crop yield assessments, correlations between the spatial distribution of soil characteristics and crop yield were observed. The spatial distribution pattern of soil properties explained 65 % of the spatial distribution pattern of corn yield. The spatial distribution pattern of clay content and percentage of soil base saturation explained most of the spatial distribution pattern of corn yield.

Spatial and linear correlations between soil and corn

The Technologies setting at Agricultural production system have the main characteristics the vertical productivity, reduced costs, soil physical, chemical and biological improvement to promote production sustainable growth. Thus, the study aimed to determine the variability and the linear and special correlations between the plant and soil attributes in order to select and indicate good representation of soil physical quality for forage productivity. In the growing season of 2006, on the Fazenda Bonança in Pereira Barreto (SP), the productivity of autumn corn forage (FDM) in an irrigated no-tillage system and the soil physical properties were analyzed. The purpose was to study the variability and the linear and spatial correlations between the plant and soil properties, to select an indicator of soil physical quality related to corn forage yield. A geostatistical grid was installed to collect soil and plant data, with 125 sampling points in an area of 2,500 m². The results show that the studied properties did not vary randomly and that data variability was low to very high, with well-defined spatial patterns, ranging from 7.8 to 38.0 m. On the other hand, the linear correlation between the plant and the soil properties was low and highly significant. The pairs forage dry matter versus microporosity and stem diameter versus bulk density were best correlated in the 0-0.10 m layer, while the other pairs - forage dry matter versus macro - and total porosity - were inversely correlated in the same layer. However, from the spatial point of view, there was a high inverse correlation between forage dry matter with microporosity, so that microporosity in the 0-0.10 m layer can be considered a good indicator of soil physical quality, with a view to corn forage yield.

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.

Mineralogy of the clay fraction of Alfisols in two slope curvatures: III - spatial variability

A good knowledge of the spatial distribution of clay minerals in the landscape facilitates the understanding of the influence of relief on the content and crystallographic attributes of soil minerals such as goethite, hematite, kaolinite and gibbsite. This study aimed at describing the relationships between the mineral properties of the clay fraction and landscape shapes by determining the mineral properties of goethite, hematite, kaolinite and gibbsite, and assessing their dependence and spatial variability, in two slope curvatures. To this end, two 100 × 100 m grids were used to establish a total of 121 regularly spaced georeferenced sampling nodes 10 m apart. Samples were collected from the layer 0.0-0.2 m and analysed for iron oxides, and kaolinite and gibbsite in the clay fraction. Minerals in the clay fraction were characterized from their X-ray diffraction (XRD) patterns, which were interpreted and used to calculate the width at half height (WHH) and mean crystallite dimension (MCD) of iron oxides, kaolinite, and gibbsite, as well as aluminium substitution and specific surface area (SSA) in hematite and goethite. Additional calculations included the goethite and hematite contents, and the goethite/(goethite+hematite) [Gt/(Gt+Hm)] and kaolinite/(kaolinite+gibbsite) [Kt/(Kt+Gb)] ratios. Mineral properties were established by statistical analysis of the XRD data, and spatial dependence was assessed geostatistically. Mineralogical properties differed significantly between the convex area and concave area. The geostatistical analysis showed a greater number of mineralogical properties with spatial dependence and a higher range in the convex than in the concave area.

Mineralogy of the clay fraction of alfisols in two slope curvatures: IV - spatial correlation with physical properties

Although the influence of clay mineralogy on soil physical properties has been widely studied, spatial relationships between these features in Alfisols have rarely been examined. The purpose of this work was to relate the clay minerals and physical properties of an Alfisol of sandstone origin in two slope curvatures. The crystallographic properties such as mean crystallite size (MCS) and width at half height (WHH) of hematite, goethite, kaolinite and gibbsite; contents of hematite and goethite; aluminium substitution (AS) and specific surface area (SSA) of hematite and goethite; the goethite/(goethite+hematite) and kaolinite/(kaolinite+gibbsite) ratios; and the citrate/bicarbonate/dithionite extractable Fe (Fe d) were correlated with the soil physical properties through Pearson correlation coefficients and cross-semivariograms. The correlations found between aluminium substitution in goethite and the soil physical properties suggest that the degree of crystallinity of this mineral influences soil properties used as soil quality indicators. Thus, goethite with a high aluminium substitution resulted in large aggregate sizes and a high porosity, and also in a low bulk density and soil penetration resistance. The presence of highly crystalline gibbsite resulted in a high density and micropore content, as well as in smaller aggregates. Interpretation of the cross-semivariogram and classification of landscape compartments in terms of the spatial dependence pattern for the relief-dependent physical and mineralogical properties of the soil proved an effective supplementary method for assessing Pearson correlations between the soil physical and mineralogical properties.

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.

Harvest managements and cultural practices in sugarcane

The presence of trash from the mechanical harvest of green cane on sugarcane plantations promotes changes in the agricultural management, for example, in the mechanical cultural practices of ratoon cane in-between the rows and nitrogen (N) fertilization. The goal of this study was to evaluate the performance of sugarcane in different harvest systems, associated to the mechanical cultural practices in interrows and N rates. The study was carried out on a sugarcane plantation in Sales Oliveira, São Paulo, Brazil, with the sugarcane variety SP81-3250, on soil classified as Acrudox, in a randomized block design with split-split plots and four replications. The main treatments consisted of harvest systems (harvesting green cane or burnt cane), the secondary treatment consisted of the mechanical cultural practices in the interrows and the tertiary treatments were N rates (0, 30, 60, 90, 120 and 160 kg ha-1), using ammonium nitrate (33 % N) as N source. The harvest systems did not differ in sugarcane yield (tons of cane per hectare - TCH), but in burnt cane, the pol percent and total sugar recovery (TSR) were higher. This could be explained by the higher quantity of plant impurities in the harvested raw material in the system without burning, which reduces the processing quality. Mechanical cultural practices in the interrows after harvest had no effect on cane yield and sugar quality, indicating that this operation can be omitted in areas with mechanical harvesting. The application of N fertilizer at rates of 88 and 144 kg ha-1 N, respectively, increased stalk height and TCH quadratically to the highest values for these variables. For the sugar yield per hectare (in pol %), N fertilization induced a linear increase.

Impacts of land leveling on lowland soil physical properties

The practice of land leveling alters the soil surface to create a uniform slope to improve land conditions for the application of all agricultural practices. The aims of this study were to evaluate the impacts of land leveling through the magnitudes, variances and spatial distributions of selected soil physical properties of a lowland area in the State of Rio Grande do Sul, Brazil; the relationships between the magnitude of cuts and/or fills and soil physical properties after the leveling process; and evaluation of the effect of leveling on the spatial distribution of the top of the B horizon in relation to the soil surface. In the 0-0.20 m layer, a 100-point geo-referenced grid covering two taxonomic soil classes was used in assessment of the following soil properties: soil particle density (Pd) and bulk density (Bd); total porosity (Tp), macroporosity (Macro) and microporosity (Micro); available water capacity (AWC); sand, silt, clay, and dispersed clay in water (Disp clay) contents; electrical conductivity (EC); and weighted average diameter of aggregates (WAD). Soil depth to the top of the B horizon was also measured before leveling. The overall effect of leveling on selected soil physical properties was evaluated by paired "t" tests. The effect on the variability of each property was evaluated through the homogeneity of variance test. The thematic maps constructed by kriging or by the inverse of the square of the distances were visually analyzed to evaluate the effect of leveling on the spatial distribution of the properties and of the top of the B horizon in relation to the soil surface. Linear regression models were fitted with the aim of evaluating the relationship between soil properties and the magnitude of cuts and fills. Leveling altered the mean value of several soil properties and the agronomic effect was negative. The mean values of Bd and Disp clay increased and Tp, Macro and Micro, WAD, AWC and EC decreased. Spatial distributions of all soil physical properties changed as a result of leveling and its effect on all soil physical properties occurred in the whole area and not specifically in the cutting or filling areas. In future designs of leveling, we recommend overlaying a cut/fill map on the map of soil depth to the top of the B horizon in order to minimize areas with shallow surface soil after leveling.

Spatial estimation of foliar phosphorus in different species of the genus Coffea based on soil properties

Information underlying analyses of coffee fertilization systems should consider both the soil and the nutritional status of plants. This study investigated the spatial relationship between phosphorus (P) levels in coffee plant tissues and soil chemical and physical properties. The study was performed using two arabica and one canephora coffee variety. Sampling grids were established in the areas, and the points georeferenced. The assessed properties of the soil were levels of available phosphorus (P-Mehlich), remaining phosphorus (P-rem) and particle size, and of the plant tissue, phosphorus levels (foliar P). The data were subjected to descriptive statistical analysis, correlation analysis, cluster analysis, and probability tests. Geostatistical and trend analyses were only performed for pairs of variables with significant linear correlation. The spatial variability for foliar P content was high for the variety Catuai and medium for the other evaluated plants. Unlike P-Mehlich, the variability in P-rem of the soil indicated the nutritional status of this nutrient in the plant.