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.

Characterization of soil chemical properties of strawberry fields using principal component analysis

One of the largest strawberry-producing municipalities of Rio Grande do Sul (RS) is Turuçu, in the South of the State. The strawberry production system adopted by farmers is similar to that used in other regions in Brazil and in the world. The main difference is related to the soil management, which can change the soil chemical properties during the strawberry cycle. This study had the objective of assessing the spatial and temporal distribution of soil fertility parameters using principal component analysis (PCA). Soil sampling was based on topography, dividing the field in three thirds: upper, middle and lower. From each of these thirds, five soil samples were randomly collected in the 0-0.20 m layer, to form a composite sample for each third. Four samples were taken during the strawberry cycle and the following properties were determined: soil organic matter (OM), soil total nitrogen (N), available phosphorus (P) and potassium (K), exchangeable calcium (Ca) and magnesium (Mg), soil pH (pH), cation exchange capacity (CEC) at pH 7.0, soil base (V%) and soil aluminum saturation(m%). No spatial variation was observed for any of the studied soil fertility parameters in the strawberry fields and temporal variation was only detected for available K. Phosphorus and K contents were always high or very high from the beginning of the strawberry cycle, while pH values ranged from very low to very high. Principal component analysis allowed the clustering of all strawberry fields based on variables related to soil acidity and organic matter content.

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.

Fuzzy logic applied to the modeling of water dynamics in an Oxisol in northeastern Brazil

Modeling of water movement in non-saturated soil usually requires a large number of parameters and variables, such as initial soil water content, saturated water content and saturated hydraulic conductivity, which can be assessed relatively easily. Dimensional flow of water in the soil is usually modeled by a nonlinear partial differential equation, known as the Richards equation. Since this equation cannot be solved analytically in certain cases, one way to approach its solution is by numerical algorithms. The success of numerical models in describing the dynamics of water in the soil is closely related to the accuracy with which the water-physical parameters are determined. That has been a big challenge in the use of numerical models because these parameters are generally difficult to determine since they present great spatial variability in the soil. Therefore, it is necessary to develop and use methods that properly incorporate the uncertainties inherent to water displacement in soils. In this paper, a model based on fuzzy logic is used as an alternative to describe water flow in the vadose zone. This fuzzy model was developed to simulate the displacement of water in a non-vegetated crop soil during the period called the emergency phase. The principle of this model consists of a Mamdani fuzzy rule-based system in which the rules are based on the moisture content of adjacent soil layers. The performances of the results modeled by the fuzzy system were evaluated by the evolution of moisture profiles over time as compared to those obtained in the field. The results obtained through use of the fuzzy model provided satisfactory reproduction of soil moisture profiles.

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.

SPATIAL VARIABILITY AND VITALITY OF EPIGEOUS TERMITE MOUNDS IN PASTURES OF MATO GROSSO DO SUL, BRAZIL

Epigeous termite mounds are frequently observed in pasture areas, but the processes regulating their population dynamics are poorly known. This study evaluated epigeous termite mounds in cultivated grasslands used as pastures, assessing their spatial distribution by means of geostatistics and evaluating their vitality. The study was conducted in the Cerrado biome in the municipality of Rio Brilhante, Mato Grosso do Sul, Brazil. In two pasture areas (Pasture 1 and Pasture 2), epigeous mounds (nests) were georeferenced and analyzed for height, circumference and vitality (inhabited or not). The area occupied by the mounds was calculated and termite specimens were collected for taxonomic identification. The spatial distribution pattern of the mounds was analyzed with geostatistical procedures. In both pasture areas, all epigeous mounds were built by the same species, Cornitermes cumulans. The mean number of mounds per hectare was 68 in Pasture 1 and 127 in Pasture 2, representing 0.4 and 1 % of the entire area, respectively. A large majority of the mounds were active (vitality), 91 % in Pasture 1 and 84 % in Pasture 2. A “pure nugget effect” was observed in the semivariograms of height and nest circumference in both pastures reflecting randomized spatial distribution and confirming that the distribution of termite mounds in pastures had a non-standard distribution.

SPATIAL CORRELATION BETWEEN PHYSICAL PROPERTIES OF SOIL AND WEEDS IN TWO MANAGEMENT SYSTEMS

The spatial correlation between soil properties and weeds is relevant in agronomic and environmental terms. The analysis of this correlation is crucial for the interpretation of its meaning, for influencing factors such as dispersal mechanisms, seed production and survival, and the range of influence of soil management techniques. This study aimed to evaluate the spatial correlation between the physical properties of soil and weeds in no-tillage (NT) and conventional tillage (CT) systems. The following physical properties of soil and weeds were analyzed: soil bulk density, macroporosity, microporosity, total porosity, aeration capacity of soil matrix, soil water content at field capacity, weed shoot biomass, weed density, Commelina benghalensis density, and Bidens pilosa density. Generally, the ranges of the spatial correlations were higher in NT than in CT. The cross-variograms showed that many variables have a structure of combined spatial variation and can therefore be mapped from one another by co-kriging. This combined variation also allows inferences about the physical and biological meanings of the study variables. Results also showed that soil management systems influence the spatial dependence structure significantly.