Implications of the variability in soil penetration resistance for statistical analysis

Soil penetration resistance (PR) is a measure of soil compaction closely related to soil structure and plant growth. However, the variability in PR hampers the statistical analyses. This study aimed to evaluate the variability of soil PR on the efficiency of parametric and nonparametric analyses in indentifying significant effects of soil compaction and to classify the coefficient of variation of PR into low, medium, high and very high. On six dates, the PR of a typical dystrophic Red Ultisol under continuous no-tillage for 16 years was measured. Three tillage and/or traffic conditions were established with the application of: (i) no chiseling or additional traffic, (ii) additional compaction, and (iii) chiseling. On each date, the nineteen PR data (measured at every 1.5 cm to a depth of 28.5 cm) were grouped in layers with different thickness. In each layer, the treatment effects were evaluated by variance (ANOVA) and Kruskal-Wallis analyses in a completely randomized design, and the coefficients of variation of all analyses were classified (low, intermediate, high and very high). The ANOVA performed better in discriminating the compaction effects, but the rejection rate of null hypothesis decreased from 100 to 80 % when the coefficient of variation increased from 15 to 26 %. The values of 15 and 26 % were the thresholds separating the low/intermediate and the high/very high coefficient variation classes of PR in this Ultisol.

Soil chemical properties related to acidity under successive pig slurry application

Pig slurry application as soil manure can alter the chemical properties of the soil and affect its acidity, modifying the environment for crop growth and development. The objective of this study was to evaluate the chemical properties related to soil acidity subjected to successive applications of pig slurry. The experiment was conducted in May 2000, in an experimental area of the Federal University of Santa Maria (UFSM) under no-tillage and lasted until January 2008. Nineteen surface applications of 0, 20, 40, and 80 m³ ha-1 of pig slurry were performed, during a period of 100 months and the soil sampled in the end (layers 0-2, 2-4, 4-6, 6-8, 8-10, 10-12, 12-14, 14-16, 16-18, 18-20, 20-25, 25-30, 30-35, 35-40, 40-50 and 50-60 cm). The application of pig slurry increased soil pH values, an effect that could reach the depth of 8 cm without affecting the potential acidity values. The applications also resulted in accumulation of Ca and Mg exchangeable levels in the surface layers, increasing base saturation and reducing Al saturation. Long-term applications induced an increase in organic matter in the deeper layers. However, the effect of this residue on the potential CEC was less significant and restricted to the surface layers.

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.

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.

Soil losses from typic cambisols and red latosol as related to three erosive rainfall patterns

Rainfall erosivity is one of the main factors related to water erosion in the tropics. This work focused on relating soil loss from a typic dystrophic Tb Haplic Cambisol (CXbd) and a typic dystrophic Red Latosol (LVdf) to different patterns of natural erosive rainfall. The experimental plots of approximately 26 m² (3 x 8.67 m) consisted of a CXbd area with a 0.15 m m-1 slope and a LVdf area with 0.12 m m-1 slope, both delimited by galvanized plates. Drainpipes were installed at the lower part of these plots to collect runoff, interconnected with a Geib or multislot divisor. To calculate erosivity (EI30), rainfall data, recorded continuously at a weather station in Lavras, were used. The data of erosive rainfall events were measured (10 mm precipitation intervals, accuracy 0.2 mm, 24 h period, 20 min intervals), characterized as rainfall events with more than 10 mm precipitation, maximum intensity > 24 mm h-1 within 15 min, or kinetic energy > 3.6 MJ, which were used in this study to calculate the rainfall erosivity parameter, were classified according to the moment of peak precipitation intensity in advanced, intermediate and delayed patterns. Among the 139 erosive rainfall events with CXbd soil loss, 60 % were attributed to the advanced pattern, with a loss of 415.9 Mg ha-1, and total losses of 776.0 Mg ha-1. As for the LVdf, of the 93 erosive rainfall events with soil loss, 58 % were listed in the advanced pattern, with 37.8 Mg ha-1 soil loss and 50.9 Mg ha-1 of total soil loss. The greatest soil losses were observed in the advanced rain pattern, especially for the CXbd. From the Cambisol, the soil loss per rainfall event was greatest for the advanced pattern, being influenced by the low soil permeability.

Quality of surface water related to land use: a case study in a catchment with small farms and intensive vegetable crop production in southern Brazil

Water degradation is strongly related to agricultural activity. The aim of this study was to evaluate the influence of land use and some environmental components on surface water quality in the Campestre catchment, located in Colombo, state of Parana, Brazil. Physical and chemical attributes were analyzed (total nitrogen, ammonium, nitrate, total phosphorus, electrical conductivity, pH, temperature, turbidity, total solids, biological oxygen demand, chemical oxygen demand and dissolved oxygen). Monthly samples of the river water were taken over one year at eight monitoring sites, distributed over three sub-basins. Overall, water quality was worse in the sub-basin with a higher percentage of agriculture, and was also affected by a lower percentage of native forest and permanent preservation area, and a larger drainage area. Water quality was also negatively affected by the presence of agriculture in the riparian zone. In the summer season, probably due to higher rainfall and intensive soil use, a higher concentration of total nitrogen and particulate nitrogen was observed, as well as higher electrical conductivity, pH and turbidity. All attributes, except for total phosphorus, were in compliance with Brazilian Conama Resolution Nº 357/2005 for freshwater class 1. However, it should be noted that these results referred to the base flow and did not represent a discharge condition since most of the water samples were not collected at or near the rainfall event.

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.

PEDO-TRANSFER FUNCTIONS FOR ESTIMATING SOIL BULK DENSITY IN CENTRAL AMAZONIA

Under field conditions in the Amazon forest, soil bulk density is difficult to measure. Rigorous methodological criteria must be applied to obtain reliable inventories of C stocks and soil nutrients, making this process expensive and sometimes unfeasible. This study aimed to generate models to estimate soil bulk density based on parameters that can be easily and reliably measured in the field and that are available in many soil-related inventories. Stepwise regression models to predict bulk density were developed using data on soil C content, clay content and pH in water from 140 permanent plots in terra firme (upland) forests near Manaus, Amazonas State, Brazil. The model results were interpreted according to the coefficient of determination (R2) and Akaike information criterion (AIC) and were validated with a dataset consisting of 125 plots different from those used to generate the models. The model with best performance in estimating soil bulk density under the conditions of this study included clay content and pH in water as independent variables and had R2 = 0.73 and AIC = -250.29. The performance of this model for predicting soil density was compared with that of models from the literature. The results showed that the locally calibrated equation was the most accurate for estimating soil bulk density for upland forests in the Manaus region.