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.

Wheat yield and physical properties of a brown latosol under no-tillage in south-central Paraná

Soil management influences the chemical and physical properties of soil. Chemical conditions have been thoroughly studied, while the role of soil physical conditions regarding crop yield has been neglected. This study aimed to analyze the wheat yield and its relationship with physical properties of an Oxisol under no-tillage (NT). The study was carried out between 2010 and 2011, in Reserva do Iguaçu, State of Paraná, Brazil, on the Campo Bonito farm, after 25 years of NT management. Based on harvest maps of barley (2006), wheat (2007) and maize (2009) of a plot (150 ha), zones with higher and lower yield potential (Z1 and Z2, respectively) were identified. Sampling grids with 16 units (50 x 50 m) and three sampling points per unit were established. The wheat grain yield (GY) and water infiltration capacity (WIC) were evaluated in 2010. Soil samples with disturbed and undisturbed structure were collected from the 0.00-0.10 and 0.10-0.20 m layers. The former were used to determine soil organic carbon (Corg) levels and the latter to determine soil bulk density (BD), total porosity (TP), macroporosity (Mac), and microporosity (Mic). Soil penetration resistance (PR) and water content (SWC) were also evaluated. The wheat GY of the whole plot was close to the regional average and the yield between the zones differed significantly, i.e. 22 % higher in Z1 than in Z2. No significant variation in Mic was observed between zones, but Z1 had higher Corg levels, SWC, TP and Mac and lower BD than Z2 in both soil layers, as well as a lower PR than Z2 in the 0.00-0.10 m layer. Therefore, soil physical conditions were more restrictive in Z2, in agreement with wheat yield and zone yield potential defined a priori, based on the harvest maps. Soil WIC in Z1 was significantly higher (30 %) than in Z2, in agreement with the results of TP and Mac which were also higher in Z1 in both soil layers. The correlation analysis of data of the two layers showed a positive relationship between wheat GY and the soil properties TP, SWC and WIC.

The effect of intrinsic soil properties on soil quality assessments

The assessment of soil quality is based on indicators and indices derived from soil properties. However, intrinsic soil properties may interfere with other soil properties that vary under different land uses and are used to calculate the indices. The aim of this study was to assess the extent to which intrinsic soil properties (clay and iron oxide contents) explain variable soil properties (sum of bases, potential acidity, organic carbon, total porosity, and bulk density) under different land uses (native forest, no-tillage and conventional agriculture) on small family farms in Southern Brazil. The results showed that the five properties evaluated can be included in soil quality assessments and are not influenced by the clay and iron oxide contents. It was concluded that for little weathered 1:1 and 2:1 phyllosilicate rich-soils, if the difference between the maximum and the minimum clay content under the different land uses is less than about 200 g kg-1 and the iron oxide content less than about 15 g kg-1, the physico-chemical soil properties in the surface layer are determined mostly by the land use.

Soil physical properties and grape yield influenced by cover crops and management systems

The use of cover crops in vineyards is a conservation practice with the purpose of reducing soil erosion and improving the soil physical quality. The objective of this study was to evaluate cover crop species and management systems on soil physical properties and grape yield. The experiment was carried out in Bento Gonçalves, RS, Southern Brazil, on a Haplic Cambisol, in a vineyard established in 1989, using White and Rose Niagara grape (Vitis labrusca L.) in a horizontal, overhead trellis system. The treatments were established in 2002, consisting of three cover crops: spontaneous species (SS), black oat (Avena strigosa Schreb) (BO), and a mixture of white clover (Trifolium repens L.), red clover (Trifolium pratense L.) and annual rye-grass (Lolium multiflorum L.) (MC). Two management systems were applied: desiccation with herbicide (D) and mechanical mowing (M). Soil under a native forest (NF) area was collected as a reference. The experimental design consisted of completely randomized blocks, with three replications. The soil physical properties in the vine rows were not influenced by cover crops and were similar to the native forest, with good quality of the soil structure. In the inter-rows, however, there was a reduction in biopores, macroporosity, total porosity and an increase in soil density, related to the compaction of the surface soil layer. The M system increased soil aggregate stability compared to the D system. The treatments affected grapevine yield only in years with excess or irregular rainfall.

Effects of tillage systems on physical properties of a cohesive yellow argisol in the northern state of Espírito Santo, Brazil

Tillage systems are a key element of the technology of crop production, both with a view to crop yield and from the perspective of soil conservation and sustainability of the production system. The aim of this paper was to evaluate the effects of five tillage systems on the physical properties of a cohesive Yellow Argisol. The experiment was installed in the field on January 21, 2011 and lasted 260 days, in an area previously used as pasture with Brachiaria grass without liming or fertilization, but irrigated by a low pressure spray system. The treatments, in five replications and in a randomized block design, consisted of: 1) disk plow (twice) + disk harrow + ridge-furrow tillage (raising a ridge along the planting row), 135 days after transplanting (DP + RID); 2) disk plow (twice) + disk harrow (DP no RID); 3) subsoiler (SB); 4) disk plow (twice) + disk harrow + scarification with three shanks along the plant row (DP + SPR); and 5) disk plow (twice) + disk harrow + scarification with three shanks in the total area (DP + STA). In all tillage systems, furrows were mechanically opened for the papaya plants. After the treatments, the mechanical resistance to penetration was determined, followed by soil moisture, mean weight diameter (MWD), geometric mean diameter (GMD), bulk density (BD), macroporosity (Ma), microporosity (Mi), and number of fruits per plant. There were differences in penetration resistance (PR) between treatments. The subsoiler was more effective to decrease RP to a distance of 0.35 m from the plants, perpendicular to the plant row. The scarifier resulted in a lower PR than DP or SB, even at the depth of 0.40 m, and it was more effective at greater distances perpendicular to the plant. All tillage systems induced a PR between 2.0 and 3.0 MPa at the depth with the highest concentration of papaya tree roots (0-0.25 m), improving the physical conditions to this depth. There was no statistical difference among the treatments for BD, Ma, Mi, MWD, and GMD at a depth of 0.20 m. The disk plow changed the physical properties of the soil most intensely to a depth of 0.20 m. The use of scarification, reduced tillage with a forest subsoiler, or ridge-furrow tillage did not improve the physical properties in the rhizosphere. Reduced tillage with a forest subsoiler resulted in a lower number of fruits per plant than all other treatments, which did not differ from each other.

Chemical and biological properties of phosphorus-fertilized soil under legume and grass cover (Cerrado region, Brazil)

The use of cover crops has been suggested as an effective method to maintain and/or increase the organic matter content, while maintaining and/or enhancing the soil physical, chemical and biological properties. The fertility of Cerrado soils is low and, consequently, phosphorus levels as well. Phosphorus is required at every metabolic stage of the plant, as it plays a role in the processes of protein and energy synthesis and influences the photosynthetic process. This study evaluated the influence of cover crops and phosphorus rates on soil chemical and biological properties after two consecutive years of common bean. The study analyzed an Oxisol in Selvíria (Mato Grosso do Sul, Brazil), in a randomized block, split plot design, in a total of 24 treatments with three replications. The plot treatments consisted of cover crops (millet, pigeon pea, crotalaria, velvet bean, millet + pigeon pea, millet + crotalaria, and millet + velvet bean) and one plot was left fallow. The subplots were represented by phosphorus rates applied as monoammonium phosphate (0, 60 and 90 kg ha-1 P2O5). In August 2011, the soil chemical properties were evaluated (pH, organic matter, phosphorus, potential acidity, cation exchange capacity, and base saturation) as well as biological variables (carbon of released CO2, microbial carbon, metabolic quotient and microbial quotient). After two years of cover crops in rotation with common bean, the cover crop biomass had not altered the soil chemical properties and barely influenced the microbial activity. The biomass production of millet and crotalaria (monoculture or intercropped) was highest. The biological variables were sensitive and responded to increasing phosphorus rates with increases in microbial carbon and reduction of the metabolic quotient.

Physical properties of a Red-Yellow Latosol and productivity of a signalgrass pasture fertilized with increasing nitrogen doses

The physical properties and fertility of the soil are important factors in the formation and establishment of pasture. Changes in physical properties affect the movement of water, air, nutrients and roots, which, in turn, affect the productivity and longevity of pastures. The objective of this study was to evaluate the physical properties of the soil and the dry matter yield of a pasture with signalgrass cv. Basilisk (Brachiaria decumbens cv. Basilisk), fertilized with increasing nitrogen doses (N), on a dystrophic Red-Yellow Latosol. The experiment was conducted on the Fazenda Rio Manso of the Universidade Federal dos Vales do Jequitinhonha e Mucuri, in Couto de Magalhães de Minas, State of Minas Gerais, Brazil. To evaluate the annual forage yield, a split plot scheme in a randomized block design with four replications was used, with N doses (0, 50, 100, 150, and 200 kg/ha/year) in the plots and growing seasons (first and second) in the subplots. For soil evaluation, a split plot scheme was used with N doses (0, 25, 50, 75 and 100 kg/ha/cut) in the plots and three sampling times (prior to the experiment, at the end of the first growing season and at the end of the second growing season) in the subplots in a randomized block design with four replications. This analysis was performed separately at two soil depths (0-3 and 10-13 cm). Forage samples were analyzed for the annual dry matter yield (DMY), and soil samples were analyzed for pre-consolidation pressure (σp), initial soil bulk density (Bd), total pore volume (TPV) and void index (Vd). Higher nitrogen doses increased the dry matter yield of signalgrass pasture and the pre-consolidation pressure of the soil. The total pore volume and void index decreased, and the initial soil bulk density increased, though without promoting soil compaction.

Sugarcane trash management assessed by the interaction of yield with soil properties

Currently, sugarcane plays an important global role, particularly with a view to alternative energy sources. Thus, in a sugarcane field of the mill Vale do Paraná S/A Álcool e Açúcar, Rubineia, São Paulo State, managed under two green cane harvest systems (cane trash left on and cane trash removed from the soil), Pearson and spatial correlations between the sugarcane yield (variety RB855035 in the third cut) and soil physical and chemical properties were studied to identify the property best correlated with stalk yield and the best harvest method. For this purpose, two geostatistical grids (121 sampling points on 1.30 ha) were installed on a eutrophic Red Argisol (homogeneous slope of 0.065 m m-1), in 2011, to determine the properties: stalk yield and sugarcane plant population, and soil resistance to penetration, gravimetric moisture, bulk density, and carbon stock, in the layers 0-0.20 and 0.20-0.40 m. The data were analyzed by descriptive, linear correlation and geostatistical analysis. In both treatments, the property stand density was best correlated with sugarcane yield (r = 0.725 in the trash mulching treatment - TM and r = 0.769 in the trash removal treatment - TR). However, in relation to the soil properties, bulk density (0-0.20 m) was best correlated (r = 0.305 in TM, r = 0.211 in TR). Similarly, from the spatial point of view, stand density was the property that best explained the sugarcane yield. However, in the TM treatment the density (0.20-0.40 m) was the only soil property spatially correlated with stalk yield. The carbon stock in the soil of the TM was 11.5 % higher than in the TR treatment. Results of the TM treatment were best, also with regard to soil management and conservation.

Influence of field slope and coffee plantation age on the physical properties of a red-yellow Latosol

In modern agriculture, several factors cause changes in the soil physical properties. The time of establishment of a crop (plantation age) and the slope are examples of factors that moderate the impact of mechanized operations on the soil structure. The objective of this study was to analyze the effect of machinery traffic on the physical properties of a Red-Yellow Latosol under coffee plantations with different ages (2, 7, 18, and 33 years) and slope positions (3, 9 and 15 %). Samples were collected from three positions between coffee rows (lower wheel track, inter-row and upper wheel track) and at two depths (surface layer and sub-surface). Changes in the total porosity, macroporosity, microporosity, organic matter, bulk density, and aggregate stability were investigated. Our results showed that the slope influenced the organic matter content, microporosity and aggregate stability. The soil samples under the inter-row were minimally damaged in their structure, compared to those from under the lower and upper wheel track, while the structure was better preserved under the lower than the upper track. The time since the establishment of the crop, i.e., the plantation age, was the main factor determining the extent of structural degradation in the coffee plantation.

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.

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.

Soil uses during the sugarcane fallow period: influence on soil chemical and physical properties and on sugarcane productivity

The planting of diversified crops during the sugarcane fallow period can improve the chemical and physical properties and increase the production potential of the soil for the next sugarcane cycle. The primary purpose of this study was to assess the influence of various soil uses during the sugarcane fallow period on soil chemical and physical properties and productivity after the first sugarcane harvest. The experiment was conducted in two areas located in Jaboticabal, São Paulo State, Brazil (21º 14' 05'' S, 48º 17' 09'' W) with two different soil types, namely: an eutroferric Red Latosol (RLe) with high-clay texture (clay content = 680 g kg-1) and an acric Red Latosol (RLa) with clayey texture (clay content = 440 g kg-1). A randomized block design with five replications and four treatments (crop sequences) was used. The crop sequences during the sugarcane fallow period were soybean/millet/soybean, soybean/sunn hemp/soybean, soybean/fallow/soybean, and soybean. Soil use was found not to affect chemical properties and sugarcane productivity of RLe or RLa. The soybean/millet/soybean sequence improved aggregation in the acric Latosol.

Evaluation of statistical and geostatistical models of digital soil properties mapping in tropical mountain regions

Soil properties have an enormous impact on economic and environmental aspects of agricultural production. Quantitative relationships between soil properties and the factors that influence their variability are the basis of digital soil mapping. The predictive models of soil properties evaluated in this work are statistical (multiple linear regression-MLR) and geostatistical (ordinary kriging and co-kriging). The study was conducted in the municipality of Bom Jardim, RJ, using a soil database with 208 sampling points. Predictive models were evaluated for sand, silt and clay fractions, pH in water and organic carbon at six depths according to the specifications of the consortium of digital soil mapping at the global level (GlobalSoilMap). Continuous covariates and categorical predictors were used and their contributions to the model assessed. Only the environmental covariates elevation, aspect, stream power index (SPI), soil wetness index (SWI), normalized difference vegetation index (NDVI), and b3/b2 band ratio were significantly correlated with soil properties. The predictive models had a mean coefficient of determination of 0.21. Best results were obtained with the geostatistical predictive models, where the highest coefficient of determination 0.43 was associated with sand properties between 60 to 100 cm deep. The use of a sparse data set of soil properties for digital mapping can explain only part of the spatial variation of these properties. The results may be related to the sampling density and the quantity and quality of the environmental covariates and predictive models used.

Multivariate analysis applied to the study of the relationship between soil and plant properties in a peach orchard

In the State of Rio Grande do Sul, the municipality of Pelotas is responsible for 90 % of peach production due to its suitable climate and soil conditions. However, there is the need for new studies that aim at improved fruit quality and increased yield. The aim of this study was to evaluate the relationship that exists between soil physical properties and properties in the peach plant in the years 2010 and 2011 by the technique of multivariate canonical correlation. The experiment was conducted in a peach orchard located in the municipality of Morro Redondo, RS, Brazil, where an experimental grid of 101 plants was established. In a trench dug beside each one of the 101 plants, soil samples were collected to determine silt, clay, and sand contents, soil density, total porosity, macroporosity, microporosity, and volumetric water content in the 0.00-0.10 and 0.10-0.20 m layers, as well as the depth of the A horizon. In each plant and in each year, the following properties were assessed: trunk diameter, fruit size and number of fruits per plant, average weight of the fruit per plant, fruit pulp firmness, Brix content, and yield from the orchard. Exploratory analysis of the data was undertaken by descriptive statistics, and the relationships between the physical properties of the soil and of the plant were assessed by canonical correlation analysis. The results showed that the clay and microporosity variables were those that exhibited the highest coefficients of canonical cross-loading with the plant properties in the soil layers assessed, and that the variable of mean weight of the fruit per plant was that which had the highest coefficients of canonical loading within the plant group for the two years assessed.

The effects of land use and soil management on the physical properties of an Oxisol in Southeast Brazil

Soils of the tropics are prone to a decrease in quality after conversion from native forest (FO) to a conventional tillage system (CT). However, the adoption of no-tillage (NT) and complex crop rotations may improve soil structural quality. Thus, the aim of this study was to evaluate the physical properties of an Oxisol under FO, CT, and three summer crop sequences in NT: continuous corn (NTcc), continuous soybean (NTcs), and a soybean/corn rotation (NTscr). Both NT and CT decreased soil organic carbon (SOC) content, SOC stock, water stable aggregates (WSA), geometric mean diameter (GMD), soil total porosity (TP), macroporosity (MA), and the least limiting water range (LLWR). However they increased soil bulk density (BD) and tensile strength (TS) of the aggregates when compared to soil under FO. Soil under NT had higher WSA, GMD, BD, TS and microporosty, but lower TP and MA than soil under CT. Soil under FO did not attain critical values for the LLWR, but the lower limit of the LLWR in soils under CT and NT was resistance to penetration (RP) for all values of BD, while the upper limit of field capacity was air-filled porosity for BD values greater than 1.46 (CT), 1.40 (NTscr), 1.42 (NTcc), and 1.41 (NTcs) kg dm-3. Soil under NTcc and NTcs decreased RP even with the increase in BD because of the formation of biopores. Furthermore, higher critical BD was verified under NTcc (1.62 kg dm-3) and NTcs (1.57 kg dm-3) compared to NTscr and CT (1.53 kg dm-3).