Soil sampling intensity and spatial distribution pattern of soils attributes and corn yield in no-tillage system

Taking into account that the sampling intensity of soil attributes is a determining factor for applying of concepts of precision agriculture, this study aims to determine the spatial distribution pattern of soil attributes and corn yield at four soil sampling intensities and verify how sampling intensity affects cause-effect relationship between soil attributes and corn yield. A 100-referenced point sample grid was imposed on the experimental site. Thus, each sampling cell encompassed an area of 45 m² and was composed of five 10-m long crop rows, where referenced points were considered the center of the cell. Samples were taken from at 0 to 0.1 m and 0.1 to 0.2 m depths. Soil chemical attributes and clay content were evaluated. Sampling intensities were established by initial 100-point sampling, resulting data sets of 100; 75; 50 and 25 points. The data were submitted to descriptive statistical and geostatistics analyses. The best sampling intensity to know the spatial distribution pattern was dependent on the soil attribute being studied. The attributes P and K+ content showed higher spatial variability; while the clay content, Ca2+, Mg2+ and base saturation values (V) showed lesser spatial variability. The spatial distribution pattern of clay content and V at the 100-point sampling were the ones which best explained the spatial distribution pattern of corn yield.

Outlining precision boundaries among areas with different variability standards using magnetic susceptibility and geomorphic surfaces

There is an increasing demand for detailed maps that represent in a simplified way the knowledge of the variability of a particular area or region maps. The objective was to outline precision boundaries among areas with different accuracy variability standards using magnetic susceptibility and geomorphic surfaces. The study was conducted in an area of 110 ha, which identified three compartment landscapes based on the geomorphic surfaces model. To determinate pH, organic matter, phosphorus, potassium and magnesium, the total sand and clay, 514 soil samples were collected at depths of 0-0.20 m and 0.60-0.80 m. The sum of base, cationic exchange capacity and base saturation were calculated and the magnetic susceptibility was evaluated in the laboratory using a system based on a balance of analytical precision method. Geomorphic surfaces identification allowed setting specific management areas (locations with maximum homogeneity of soil attributes). The map of spatial variability of magnetic susceptibility can be used to validate the precise boundaries among geomorphic surfaces identified in the field and infer the variability of clay content and soil base saturation.

Spatial influence of physical and chemical parameters on management zone definition in apple orchards

ABSTRACT Precision agriculture adoption in Brazilian apple orchards is still incipient. This study aimed at evaluating the spatial variability of certain soil properties as soil density, soil penetration resistance, electrical conductivity, yield, and fruit quality in an apple orchard through digital mapping, as well as assessing the correlation between these factors by means of geostatistics, establishing management zones. Forty representative points were set within 2.5 hectares of apple orchard, wherein soil samples were collected and analyzed, besides measurements of fruit quality (Brix degree, size or diameter, pulp firmness and color) to generate an overall index quality. We concluded that the fruit quality indexes, when isolated, did not show strong spatial dependence, unlike the index of fruit quality (FQI), derived from a combination of these parameters, allowing orchard planning according to management zones based on quality.

Phytosociological variability of weeds in soybean field

Understanding spatial distribution of weeds in the crop enables to perform localized herbicide applications, increasing the technical and economic efficiency of operations and reducing environmental impacts. This work aimed to characterize the spatial and phytosociological variability of weeds occurring in soybean commercial field. It was conducted in an agricultural area located at the municipality of Boa Vista das Missões - RS, during the 2010/2011 harvest season. The area, that had been managed under no-tillage with soybean monoculture (summer) for five years, was divided in regular squares of 50 x 50 m (0.25 ha), totalizing 356 points. For species identification, 0.5 x 0.5 m sample squares were used. During the survey, 1,739 individuals were identified, distributed in 19 species of 13 families. The weed species Cardiospermum halicacabum, Digitaria horizontalis, Urochloa plantaginea and Raphanus raphanistrum showed the highest population variation in the area; however, only C. halicacabum, U. plantaginea and R. raphanistrum stood out based on the Importance Index Value (IVI). Localized management strategies considering the spatial variability of weed species placed in the Magnoliopsidas and Liliopsidas group show a high potential for use in soybean crop. The results show that the sampling method through regular grid was capable of characterizing the occurrence, population density and spatial variability of weed species in soybean crop.

Sugarcane productivity correlated with physical-chemical attributes to create soil management zone

The socioeconomic importance of sugar cane in Brazil is unquestionable because it is the raw material for the production of ethanol and sugar. The accurate spatial intervention in the management of the crop, resulting zones of soil management, increases productivity as well as its agricultural yields. The spatial and Person's correlations between sugarcane attributes and physico-chemical attributes of a Typic Tropustalf were studied in the growing season of 2009, in Suzanápolis, State of São Paulo, Brazil (20°28'10'' S lat.; 50°49'20'' W long.), in order to obtain the one that best correlates with agricultural productivity. Thus, the geostatistical grid with 120 sampling points was installed to soil and data collection in a plot of 14.6 ha with second crop sugarcane. Due to their substantial and excellent linear and spatial correlations with the productivity of the sugarcane, the population of plants and the organic matter content of the soil, by evidencing substantial correlations, linear and spatial, with the productivity of sugarcane, were indicators of management zones strongly attached to such productivity.

Scaling of semivariograms and the kriging estimation of field-measured properties

Two methods were evaluated for scaling a set of semivariograms into a unified function for kriging estimation of field-measured properties. Scaling is performed using sample variances and sills of individual semivariograms as scale factors. Theoretical developments show that kriging weights are independent of the scaling factor which appears simply as a constant multiplying both sides of the kriging equations. The scaling techniques were applied to four sets of semivariograms representing spatial scales of 30 x 30 m to 600 x 900 km. Experimental semivariograms in each set successfully coalesced into a single curve by variances and sills of individual semivariograms. To evaluate the scaling techniques, kriged estimates derived from scaled semivariogram models were compared with those derived from unscaled models. Differences in kriged estimates of the order of 5% were found for the cases in which the scaling technique was not successful in coalescing the individual semivariograms, which also means that the spatial variability of these properties is different. The proposed scaling techniques enhance interpretation of semivariograms when a variety of measurements are made at the same location. They also reduce computational times for kriging estimations because kriging weights only need to be calculated for one variable. Weights remain unchanged for all other variables in the data set whose semivariograms are scaled.

Sampling strategies for tropical forest nutrient cycling studies: a case study in São Paulo, Brazil

The precise sampling of soil, biological or micro climatic attributes in tropical forests, which are characterized by a high diversity of species and complex spatial variability, is a difficult task. We found few basic studies to guide sampling procedures. The objective of this study was to define a sampling strategy and data analysis for some parameters frequently used in nutrient cycling studies, i. e., litter amount, total nutrient amounts in litter and its composition (Ca, Mg, Κ, Ν and P), and soil attributes at three depths (organic matter, Ρ content, cation exchange capacity and base saturation). A natural remnant forest in the West of São Paulo State (Brazil) was selected as study area and samples were collected in July, 1989. The total amount of litter and its total nutrient amounts had a high spatial independent variance. Conversely, the variance of litter composition was lower and the spatial dependency was peculiar to each nutrient. The sampling strategy for the estimation of litter amounts and the amount of nutrient in litter should be different than the sampling strategy for nutrient composition. For the estimation of litter amounts and the amount of nutrients in litter (related to quantity) a large number of randomly distributed determinations are needed. Otherwise, for the estimation of litter nutrient composition (related to quality) a smaller amount of spatially located samples should be analyzed. The determination of sampling for soil attributes differed according to the depth. Overall, surface samples (0-5 cm) showed high short distance spatial dependent variance, whereas, subsurface samples exhibited spatial dependency in longer distances. Short transects with sampling interval of 5-10 m are recommended for surface sampling. Subsurface samples must also be spatially located, but with transects or grids with longer distances between sampling points over the entire area. Composite soil samples would not provide a complete understanding of the relation between soil properties and surface dynamic processes or landscape aspects. Precise distribution of Ρ was difficult to estimate.

State-space approach to evaluate the relation between soil physical and chemical properties

The state-space approach is used to evaluate the relation between soil physical and chemical properties in an area cultivated with sugarcane. The experiment was carried out on a Rhodic Kandiudalf in Piracicaba, State of São Paulo, Brazil. Sugarcane was planted on an area of 0.21 ha i.e., in 15 rows 100 m long, spaced 1.4 m. Soil water content, soil organic matter, clay content and aggregate stability were sampled along a transect of 84 points, meter by meter. The state-space approach is used to evaluate how the soil water content is affected by itself and by soil organic matter, clay content, and aggregate stability of neighboring locations, in different combinations, aiming to contribute to a better understanding of the relation among these variables in the soil. Results show that soil water contents were successfully estimated by this approach. Best performances were found when the estimate of soil water content at locations i was related to soil water content, clay content and aggregate stability at locations i-1. Results also indicate that this state-space model using all series describes the soil water content better than any equivalent multiple regression equation.

Soil fertility, organic carbon and fractions of the organic matter at different distances from eucalyptus stumps

Knowledge on variations in vertical, horizontal and temporal characteristics of the soil chemical properties under eucalyptus stumps left in the soil is of fundamental importance for the management of subsequent crops. The objective of this work was to evaluate the effect of eucalyptus stumps (ES) left after cutting on the spatial variability of chemical characteristics in a dystrophic Yellow Argisol in the eastern coastal plain region of Brazil. For this purpose, ES left for 31 and 54 months were selected in two experimental areas with similar characteristics, to assess the decomposition effects of the stumps on soil chemical attributes. Soil samples were collected directly around these ES, and at distances of 30, 60, 90, 120 and 150 cm away from them, in the layers 0-10, 10-20 and 20-40 cm along the row of ES, which is in-between the rows of eucalyptus trees of a new plantation, grown at a spacing of 3 x 3 m. The soil was sampled in five replications in plots of 900 m² each and the samples analyzed for pH, available P and K (Mehlich-1), exchangeable Al, Ca and Mg, total organic carbon (TOC) and C content in humic substances (HS) and in the free light fraction. The pH values and P, K, Ca2+, Mg2+ and Al3+ contents varied between the soil layers with increasing distance from the 31 and 54-monthold stumps. The highest pH, P, K, Ca2+ and Mg2+ values and the lowest Al3+ content were found in the surface soil layer. The TOC of the various fractions of soil organic matter decreased with increasing distance from the 31 and 54-month-old ES in the 0-10 and 10-20 cm layers, indicating that the root (and stump) cycling and rhizodeposition contribute to maintain soil organic matter. The C contents of the free light fraction, of the HS and TOC fractions were higher in the topsoil layer under the ES left for 31 months due to the higher clay levels of this layer, than in those found under the 54-month-old stumps. However, highest C levels of the different fractions of soil organic matter in the topsoil layer reflect the deposition and maintenance of forest residues on the soil surface, mainly after forest harvest.

Inter-relation between soybean yield and soil compaction under degraded pasture in Brazilian Savannah

The Cerrado (Brazilian Savannah) plays an important economic and financial role in the nation, since the pastures of this biome feed cattle for half of the domestic bovine meat productivity, and its agricultural fields produce a third of the country's grain. The variability and spatial dependence between the soil physical attributes and soybean yield were evaluated in a crop rotation planted on a degraded brachiaria pasture, on a dystroferric Red Latosol of an experimental farm of the State University of São Paulo (UNESP), in the 2005/2006 growing season. The linear and spatial correlations between these attributes were also studied, to determine conditions that would allow increased agricultural productivity. In the above pasture area, a grid was installed with 124 plots, spaced 10.0 x 10.0 m and 5.0 x 5.0 m apart, in a total area of 7,500 m². From the linear and spatial point of view, the high grain yield can be explained by the number of grains per plant and soil macroporosity. The high variability observed for most soil properties indicated that the crop - livestock integration system results in environmental heterogeneity of the soil.

Multivariate analysis and geostatistics of the fertility of a humic rhodic hapludox under coffee cultivation

The spatial variability of soil and plant properties exerts great influence on the yeld of agricultural crops. This study analyzed the spatial variability of the fertility of a Humic Rhodic Hapludox with Arabic coffee, using principal component analysis, cluster analysis and geostatistics in combination. The experiment was carried out in an area under Coffea arabica L., variety Catucai 20/15 - 479. The soil was sampled at a depth 0.20 m, at 50 points of a sampling grid. The following chemical properties were determined: P, K+, Ca2+, Mg2+, Na+, S, Al3+, pH, H + Al, SB, t, T, V, m, OM, Na saturation index (SSI), remaining phosphorus (P-rem), and micronutrients (Zn, Fe, Mn, Cu and B). The data were analyzed with descriptive statistics, followed by principal component and cluster analyses. Geostatistics were used to check and quantify the degree of spatial dependence of properties, represented by principal components. The principal component analysis allowed a dimensional reduction of the problem, providing interpretable components, with little information loss. Despite the characteristic information loss of principal component analysis, the combination of this technique with geostatistical analysis was efficient for the quantification and determination of the structure of spatial dependence of soil fertility. In general, the availability of soil mineral nutrients was low and the levels of acidity and exchangeable Al were high.

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.

Optimization of the sampling scheme for maps of physical and chemical properties estimated by kriging

The sampling scheme is essential in the investigation of the spatial variability of soil properties in Soil Science studies. The high costs of sampling schemes optimized with additional sampling points for each physical and chemical soil property, prevent their use in precision agriculture. The purpose of this study was to obtain an optimal sampling scheme for physical and chemical property sets and investigate its effect on the quality of soil sampling. Soil was sampled on a 42-ha area, with 206 geo-referenced points arranged in a regular grid spaced 50 m from each other, in a depth range of 0.00-0.20 m. In order to obtain an optimal sampling scheme for every physical and chemical property, a sample grid, a medium-scale variogram and the extended Spatial Simulated Annealing (SSA) method were used to minimize kriging variance. The optimization procedure was validated by constructing maps of relative improvement comparing the sample configuration before and after the process. A greater concentration of recommended points in specific areas (NW-SE direction) was observed, which also reflects a greater estimate variance at these locations. The addition of optimal samples, for specific regions, increased the accuracy up to 2 % for chemical and 1 % for physical properties. The use of a sample grid and medium-scale variogram, as previous information for the conception of additional sampling schemes, was very promising to determine the locations of these additional points for all physical and chemical soil properties, enhancing the accuracy of kriging estimates of the physical-chemical 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.

USE OF SCALED SEMIVARIOGRAMS IN THE PLANNING SAMPLE OF SOIL PHYSICAL PROPERTIES IN SOUTHERN AMAZONAS, BRAZIL

There is a great lack of information from soil surveys in the southern part of the State of Amazonas, Brazil. The use of tools such as geostatistics may improve environmental planning, use and management. In this study, we aimed to use scaled semivariograms in sample design of soil physical properties of some environments in Amazonas. We selected five areas located in the south of the state of Amazonas, Brazil, with varied soil uses, such as forest, archaeological dark earth (ADE), pasture, sugarcane cropping, and agroforestry. Regular mesh grids were set up in these areas with 64 sample points spaced at 10 m from each other. At these points, we determined the particle size composition, soil resistance to penetration, moisture, soil bulk density and particle density, macroporosity, microporosity, total porosity, and aggregate stability in water at a depth of 0.00-0.20 m. Descriptive and geostatistical analyses were performed. The sample density requirements were lower in the pasture area but higher in the forest. We concluded that managed-environments had differences in their soil physical properties compared to the natural forest; notably, the soil in the ADE environment is physically improved in relation to the others. The physical properties evaluated showed a structure of spatial dependence with a slight variability of the forest compared to the others. The use of the range parameter of the semivariogram analysis proved to be effective in determining an ideal sample density.