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.

Assessment of pedotransfer functions for estimating soil water retention curves for the amazon region

Knowledge of the soil water retention curve (SWRC) is essential for understanding and modeling hydraulic processes in the soil. However, direct determination of the SWRC is time consuming and costly. In addition, it requires a large number of samples, due to the high spatial and temporal variability of soil hydraulic properties. An alternative is the use of models, called pedotransfer functions (PTFs), which estimate the SWRC from easy-to-measure properties. The aim of this paper was to test the accuracy of 16 point or parametric PTFs reported in the literature on different soils from the south and southeast of the State of Pará, Brazil. The PTFs tested were proposed by Pidgeon (1972), Lal (1979), Aina & Periaswamy (1985), Arruda et al. (1987), Dijkerman (1988), Vereecken et al. (1989), Batjes (1996), van den Berg et al. (1997), Tomasella et al. (2000), Hodnett & Tomasella (2002), Oliveira et al. (2002), and Barros (2010). We used a database that includes soil texture (sand, silt, and clay), bulk density, soil organic carbon, soil pH, cation exchange capacity, and the SWRC. Most of the PTFs tested did not show good performance in estimating the SWRC. The parametric PTFs, however, performed better than the point PTFs in assessing the SWRC in the tested region. Among the parametric PTFs, those proposed by Tomasella et al. (2000) achieved the best accuracy in estimating the empirical parameters of the van Genuchten (1980) model, especially when tested in the top soil layer.

Compared performance of penetrometers and effect of soil water content on penetration resistance measurements

Modern agriculture techniques have a great impact on crops and soil quality, especially by the increased machinery traffic and weight. Several devices have been developed for determining soil properties in the field, aimed at managing compacted areas. Penetrometry is a widely used technique; however, there are several types of penetrometers, which have different action modes that can affect the soil resistance measurement. The objective of this study was to compare the functionality of two penetrometry methods (manual and automated mode) in the field identification of compacted, highly mechanized sugarcane areas, considering the influence of soil water volumetric content (θ) on soil penetration resistance (PR). Three sugarcane fields on a Rhodic Eutrudrox were chosen, under a sequence of harvest systems: one manual harvest (1ManH), one mechanized harvest (1MH) and three mechanized harvests (3MH). The different degrees of mechanization were associated to cumulative compaction processes. An electronic penetrometer was used on PR measurements, so that the rod was introduced into the soil by hand (Manual) and by an electromechanical motor (Auto). The θ was measured in the field with a soil moisture sensor. Results showed an effect of θ on PR measurements and that regression models must be used to correct data before comparing harvesting systems. The rod introduction modes resulted in different mean PR values, where the "Manual" overestimated PR compared to the "Auto" mode at low θ.

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.

Effect of management systems and cover crops on organic matter dynamics of soil under vegetables

Vegetable production in conservation tillage has increased in Brazil, with positive effects on the soil quality. Since management systems alter the quantity and quality of organic matter, this study evaluated the influence of different management systems and cover crops on the organic matter dynamics of a dystrophic Red Latosol under vegetables. The treatments consisted of the combination of three soil tillage systems: no-tillage (NT), reduced tillage (RT) and conventional tillage (CT) and of two cover crops: maize monoculture and maize-mucuna intercrop. Vegetables were grown in the winter and the cover crops in the summer for straw production. The experiment was arranged in a randomized block design with four replications. Soil samples were collected between the crop rows in three layers (0.0-0.05, 0.05-0.10, and 0.10-0.30 m) twice: in October, before planting cover crops for straw, and in July, during vegetable cultivation. The total organic carbon (TOC), microbial biomass carbon (MBC), oxidizable fractions, and the carbon fractions fulvic acid (C FA), humic acid (C HA) and humin (C HUM) were determined. The main changes in these properties occurred in the upper layers (0.0-0.05 and 0.05-0.10 m) where, in general, TOC levels were highest in NT with maize straw. The MBC levels were lowest in CT systems, indicating sensitivity to soil disturbance. Under mucuna, the levels of C HA were lower in RT than NT systems, while the C FA levels were lower in RT than CT. For vegetable production, the C HUM values were lowest in the 0.05-0.10 m layer under CT. With regard to the oxidizable fractions, the tillage systems differed only in the most labile C fractions, with higher levels in NT than CT in the 0.0-0.05 m layer in both summer and winter, with no differences between these systems in the other layers. The cabbage yield was not influenced by the soil management system, but benefited from the mulch production of the preceding maize-mucuna intercrop as cover plant.

Rice husk ash as corrective of soil acidity

Rice husk ash (RHA) is a by-product from the burning of rice husk that can have favorable effects on the soil in terms of acidity correction. The objectives of this study were to determine the effective calcium carbonate equivalent (ECC) of RHA under field conditions, and establish technical criteria as a basis for estimating the overall ECC of RHA. The 12 treatments of the experiment consisted of 10 RHA dosages (0, 10, 20, 30, 40, 60, 80, 100, 120, and 140 Mg ha-1) and two references, one of which was an absolute control (AC) and the other a plot limed and fertilized according to official recommendations (recommended fertilization - RF). The soil was sampled twice (15 and 210 days after incorporating RHA), in the layers 0.00-0.10 and 0.10-0.20 m, to determine the pH(H2O) and base saturation (V%). The ECC and neutralizing value (NV) of RHA were also determined. The results showed that RHA neutralizes soil acidity, in a faster reaction than conventional limestone, despite a low ECC (around 3 %).

Three-dimensional interpolation of soil data: fertility and pedomorphological features in southern Brazil

The graphical representation of spatial soil properties in a digital environment is complex because it requires a conversion of data collected in a discrete form onto a continuous surface. The objective of this study was to apply three-dimension techniques of interpolation and visualization on soil texture and fertility properties and establish relationships with pedogenetic factors and processes in a slope area. The GRASS Geographic Information System was used to generate three-dimensional models and ParaView software to visualize soil volumes. Samples of the A, AB, BA, and B horizons were collected in a regular 122-point grid in an area of 13 ha, in Pinhais, PR, in southern Brazil. Geoprocessing and graphic computing techniques were effective in identifying and delimiting soil volumes of distinct ranges of fertility properties confined within the soil matrix. Both three-dimensional interpolation and the visualization tool facilitated interpretation in a continuous space (volumes) of the cause-effect relationships between soil texture and fertility properties and pedological factors and processes, such as higher clay contents following the drainage lines of the area. The flattest part with more weathered soils (Oxisols) had the highest pH values and lower Al3+ concentrations. These techniques of data interpolation and visualization have great potential for use in diverse areas of soil science, such as identification of soil volumes occurring side-by-side but that exhibit different physical, chemical, and mineralogical conditions for plant root growth, and monitoring of plumes of organic and inorganic pollutants in soils and sediments, among other applications. The methodological details for interpolation and a three-dimensional view of soil data are presented here.

Iron oxides and quality of organic matter in sugarcane harvesting systems

Improvements in working conditions, sustainable production, and competitiveness have led to substantial changes in sugarcane harvesting systems. Such changes have altered a number of soil properties, including iron oxides and organic matter, as well as some chemical properties, such as the maximum P adsorption capacity of the soil. The aim of this study was to characterize the relationship between iron oxides and the quality of organic matter in sugarcane harvesting systems. For that purpose, two 1 ha plots in mechanically and manually harvested fields were used to obtain soil samples from the 0.00-0.25 m soil layer at 126 different points. The mineralogical, chemical, and physical results were subjected to descriptive statistical analyses, such as the mean comparison test, as well as to multivariate statistical and principal component analyses. Multivariate tests allowed soil properties to be classified in two different groups according to the harvesting method: manual harvest with the burning of residual cane, and mechanical harvest without burning. The mechanical harvesting system was found to enhance pedoenvironmental conditions, leading to changes in the crystallinity of iron oxides, an increase in the humification of organic matter, and a relative decrease in phosphorus adsorption in this area compared to the manual harvesting system.

Competitive sorption and desorption of phosphate and citrate in clayey and sandy loam soils

The increase of organic acids in soils can reduce phosphorus sorption. The objective of the study was to evaluate the competitive sorption of P and citrate in clayey and sandy loam soils, using a stirred-flow system. Three experiments were performed with soil samples (0-20 cm layer) of clayey (RYL-cl) and sandy loam (RYL-sl) Red Yellow Latosols (Oxisols). In the first study, the treatments were arranged in a 2 × 5 factorial design, with two soil types and five combinations of phosphorus and citrate application (only P; P + citrate; and citrate applied 7, 22, 52 min before P); in the second, the treatments were arranged in a 2 × 2 factorial design, corresponding to two soils and two forms of P and citrate application (only citrate and citrate + P); and in the third study, the treatments in a 2 × 2 × 6 factorial design consisted of two soils, two extractors (citrate and water) and six incubation times. In the RYL-cl and RYL-sl, P sorption was highest (44 and 25 % of P application, respectively), in the absence of citrate application. Under citrate application, P sorption was reduced in all treatments. The combined application of citrate and P reduced P sorption to 25.8 % of the initially applied P in RYL-cl and to 16.7 % in RYL-sl, in comparison to P without citrate. Citrate sorption in RYL-cl and RYL-sl was highest in the absence of P application, corresponding to 32.0 and 30.2 % of the citrate applied, respectively. With P application, citrate sorption was reduced to 26.4 and 19.7 % of the initially applied citrate in RYL-cl and RYL-sl, respectively. Phosphorus desorption was greater when citrate was used. Phosphorus desorption with citrate and water was higher in the beginning (until 24 h of incubation of P) in RYL-cl and RYL-sl, indicating a rapid initial phase, followed by a slow release phase. This suggests that according to the contact time of P with the soil colloids, the previously adsorbed P can be released to the soil solution in the presence of competing ligands such as citrate. In conclusion, a soil management with continuous input of organic acids is desirable, in view of their potential to compete for P sorption sites, especially in rather weathered soils.

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).

Stabilization of soil hydraulic properties under a long term no-till system

The area under the no-tillage system (NT) has been increasing over the last few years. Some authors indicate that stabilization of soil physical properties is reached after some years under NT while other authors debate this. The objective of this study was to determine the effect of the last crop in the rotation sequence (1st year: maize, 2nd year: soybean, 3rd year: wheat/soybean) on soil pore configuration and hydraulic properties in two different soils (site 1: loam, site 2: sandy loam) from the Argentinean Pampas region under long-term NT treatments in order to determine if stabilization of soil physical properties is reached apart from a specific time in the crop sequence. In addition, we compared two procedures for evaluating water-conducting macroporosities, and evaluated the efficiency of the pedotransfer function ROSETTA in estimating the parameters of the van Genuchten-Mualem (VGM) model in these soils. Soil pore configuration and hydraulic properties were not stable and changed according to the crop sequence and the last crop grown in both sites. For both sites, saturated hydraulic conductivity, K0, water-conducting macroporosity, εma, and flow-weighted mean pore radius, R0ma, increased from the 1st to the 2nd year of the crop sequence, and this was attributed to the creation of water-conducting macropores by the maize roots. The VGM model adequately described the water retention curve (WRC) for these soils, but not the hydraulic conductivity (K) vs tension (h) curve. The ROSETTA function failed in the estimation of these parameters. In summary, mean values of K0 ranged from 0.74 to 3.88 cm h-1. In studies on NT effects on soil physical properties, the crop effect must be considered.

Effects of pig slurry application on soil physical and chemical properties and glyphosate mobility

Pig slurry applied to soil at different rates may affect soil properties and the mobility of chemical compounds within the soil. The purpose of this study was to evaluate the effects of rates of pig slurry application in agricultural areas on soil physical and chemical properties and on the mobility of glyphosate through the soil profile. The study was carried out in the 12th year of an experiment with pig slurry applied at rates of 0 (control), 50, 100 and 200 m³ ha-1 yr-1 on a Latossolo Vermelho distrófico (Hapludox) soil. In the control, the quantities of P and K removed by harvested grains were replaced in the next crop cycle. Soil physical properties (bulk density, porosity, texture, and saturated hydraulic conductivity) and chemical properties (organic matter, pH, extractable P, and exchangeable K) were measured. Soil solution samples were collected at depths of 20, 40 and 80 cm using suction lysimeters, and glyphosate concentrations were measured over a 60-day period after slurry application. Soil physical and chemical properties were little affected by the pig slurry applications, but soil pH was reduced and P levels increased in the surface layers. In turn, K levels were increased in sub-surface layers. Glyphosate concentrations tended to decrease over time but were not affected by pig slurry application. The concentrations of glyphosate found in different depths show that the pratice of this application in agricultural soils has the potential for contamination of groundwater, especially when the water table is the surface and heavy rains occur immediately after application.

Soil fertility, nutrition and yield of maize and barley with gypsum application on soil surface in no-till

Annual crop yield and nutrition have shown differentiated responses to modifications in soil chemical properties brought about by gypsum application. The aim of this study was to evaluate the effect of gypsum application rates on the chemical properties of a Latossolo Bruno (Clayey Oxisol), as well as on the nutrition and yield of a maize-barley succession under no-till. The experiment was set up in November 2009 in Guarapuava, Parana, Brazil, applying gypsum rates of 0.0, 1.5, 3.0, 4.5, and 6.0 Mg ha-1 to the soil surface upon sowing maize, with crop succession of barley. Gypsum application decreased the levels of Al3+ and Mg2+ in the 0.0-0.1 m layer and increased soil pH in the layers from 0.2-0.6 m depth. Gypsum application has increased the levels of Ca2+ in all soil layers up to 0.6 m, and the levels of S-SO4(2-) up to 0.8 m. In both crops, the leaf concentrations of Ca and S were increased while Mg concentrations have decreased as a function of gypsum rates. There was also an effect of gypsum rates on grain yield, with a quadratic response of maize and a linear increase for barley. Yield increases were up to 11 and 12 % in relation to control for the maximum technical efficiency (MTE) rates of 3.8 and 6.0 Mg ha-1 of gypsum, respectively. Gypsum application improved soil fertility in the profile, especially in the subsurface, as well as plant nutrition, increasing the yields of maize and barley.

Artificial neural networks applied for soil class prediction in mountainous landscape of the Serra do Mar¹

Soil information is needed for managing the agricultural environment. The aim of this study was to apply artificial neural networks (ANNs) for the prediction of soil classes using orbital remote sensing products, terrain attributes derived from a digital elevation model and local geology information as data sources. This approach to digital soil mapping was evaluated in an area with a high degree of lithologic diversity in the Serra do Mar. The neural network simulator used in this study was JavaNNS and the backpropagation learning algorithm. For soil class prediction, different combinations of the selected discriminant variables were tested: elevation, declivity, aspect, curvature, curvature plan, curvature profile, topographic index, solar radiation, LS topographic factor, local geology information, and clay mineral indices, iron oxides and the normalized difference vegetation index (NDVI) derived from an image of a Landsat-7 Enhanced Thematic Mapper Plus (ETM+) sensor. With the tested sets, best results were obtained when all discriminant variables were associated with geological information (overall accuracy 93.2 - 95.6 %, Kappa index 0.924 - 0.951, for set 13). Excluding the variable profile curvature (set 12), overall accuracy ranged from 93.9 to 95.4 % and the Kappa index from 0.932 to 0.948. The maps based on the neural network classifier were consistent and similar to conventional soil maps drawn for the study area, although with more spatial details. The results show the potential of ANNs for soil class prediction in mountainous areas with lithological diversity.

Index of soil physical quality of hardsetting soils on the brazilian coast

Many soils have a hard-setting behavior, also known as cohesive or "coesos". In such soils, the penetration resistance increases markedly when dry and decreases considerably when moist, creating serious limitations for plant emergence and growth. To evaluate the level of structure degradation in hard-setting soils with different texture classes and to create an index for assessing soil hardness levels in hard-setting soils, six soil representative profiles were selected in the field in various regions of Brazil. The following indices were tested: S, which measures soil physical quality, and H , which analyzes the degree of hardness and the effective stress in the soil during drying. Both indices were calculated using previously described functions based on data from the water-retention curves for the soils. The hard-setting values identified in different soils of the Brazilian Coastal Tablelands have distinct compaction (hardness) levels and can be satisfactorily measured by the H index. The S index was adequate for evaluating the structural characteristics of the hard-setting soils, classifying them as suitable or poor for cultivation, but only when the moisture level of the soil was near the inflection point. The H index showed that increases in density in hard-setting soils result from increases in effective stress and not from the soil texture. Values for Bd > 1.48 kg dm-3 classify the soil as hard-setting, and the structural organization is considered "poor".