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.

Estimation of soil moisture in the root-zone from remote sensing data

Field-based soil moisture measurements are cumbersome. Thus, remote sensing techniques are needed because allows field and landscape-scale mapping of soil moisture depth-averaged through the root zone of existing vegetation. The objective of the study was to evaluate the accuracy of an empirical relationship to calculate soil moisture from remote sensing data of irrigated soils of the Apodi Plateau, in the Brazilian semiarid region. The empirical relationship had previously been tested for irrigated soils in Mexico, Egypt, and Pakistan, with promising results. In this study, the relationship was evaluated from experimental data collected from a cotton field. The experiment was carried out in an area of 5 ha with irrigated cotton. The energy balance and evaporative fraction (Λ) were measured by the Bowen ratio method. Soil moisture (θ) data were collected using a PR2 - Profile Probe (Delta-T Devices Ltd). The empirical relationship was tested using experimentally collected Λ and θ values and was applied using the Λ values obtained from the Surface Energy Balance Algorithm for Land (SEBAL) and three TM - Landsat 5 images. There was a close correlation between measured and estimated θ values (p<0.05, R² = 0.84) and there were no significant differences according to the Student t-test (p<0.01). The statistical analyses showed that the empirical relationship can be applied to estimate the root-zone soil moisture of irrigated soils, i.e. when the evaporative fraction is greater than 0.45.

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.

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.

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.

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.

SPATIAL VARIABILITY OF SOIL PROPERTIES IN AN AGRARIAN REFORM SETTLEMENT

ABSTRACT The study of soil chemical and physical properties variability is important for suitable management practices. The aim of this study was to evaluate the spatial variability of soil properties in the Malhada do Meio settlement to subsidize soil use planning. The settlement is located in Chapadinha, MA, Brazil, and has an area of 630.86 ha. The vegetation is seasonal submontane deciduous forest and steppe savanna. The geology is formed of sandstones and siltstones of theItapecuru Formation and by colluvial and alluvial deposits. The relief consists of hills with rounded and flat tops with an average altitude of 67 m, and frequently covered over by ferruginous duricrusts. A total of 183 georeferenced soil samples were collected at the depth of 0.00-0.20 m inPlintossolos, Neossolo andGleissolo. The following chemical variables were analyzed: pH(CaCl2), H+Al, Al, SB, V, CEC, P, K, OM, Ca, Mg, SiO2, Al2O3, and Fe2O3; along with particle size variables: clay, silt, and sand. Descriptive statistical and geostatistical analyses were carried out. The coefficient of variation (CV) was high for most of the variables, with the exception of pH with a low CV, and of sand with a medium CV. The models fitted to the experimental semivariograms of these variables were the exponential and the spherical. The range values were from 999 m to 3,690 m. For the variables pH(CaCl2), SB, and clay, there are three specific areas for land use planning. The central part of the area (zone III), where thePlintossolos Pétricos and Neossolos Flúvicos occur, is the most suitable for crops due to higher macronutrient content, organic matter and pH. Zones I and II are indicated for environmental preservation.

Effect of soil-water tension on herbaceous cotton yield

A field experiment was conducted during two years, 1990/91, in an alluvial soil, in the State of Paraíba, Brazil, to study the effect of the levels of soil-water tension, 50, 100, 200, 300, 400 and 600 kPa, at 20 cm depth, on upland cotton (Gossypium hirsutum L.r. latifolium Hutch, cv. CNPA-6H) yield. The experimental design was a complete randomized block with six treatments and four repetitions. There was an effect of the treatments on plant height, leaf area index and cotton yield, but the precocity index was not modified. Water should be applied when the soil-water tension, measured at 20 cm depth, reaches values around 200 kPa. There was a quadratic (R² = 0.893**) response of cotton yields to soil water tension, with the maximum when water was applied at 52% of soil water depletion.

Bacterial diversity of soil under eucalyptus assessed by 16S rDNA sequencing analysis

Studies on the impact of Eucalyptus spp. on Brazilian soils have focused on soil chemical properties and isolating interesting microbial organisms. Few studies have focused on microbial diversity and ecology in Brazil due to limited coverage of traditional cultivation and isolation methods. Molecular microbial ecology methods based on PCR amplified 16S rDNA have enriched the knowledge of soils microbial biodiversity. The objective of this work was to compare and estimate the bacterial diversity of sympatric communities within soils from two areas, a native forest (NFA) and an eucalyptus arboretum (EAA). PCR primers, whose target soil metagenomic 16S rDNA were used to amplify soil DNA, were cloned using pGEM-T and sequenced to determine bacterial diversity. From the NFA soil 134 clones were analyzed, while 116 clones were analyzed from the EAA soil samples. The sequences were compared with those online at the GenBank. Phylogenetic analyses revealed differences between the soil types and high diversity in both communities. Soil from the Eucalyptus spp. arboretum was found to have a greater bacterial diversity than the soil investigated from the native forest area.

Bioindicators of soil quality in coffee organic cultivation systems

The objective of this work was to assess the effect of different coffee organic cultivation systems on chemical and biological soil characteristics, in different seasons of the year. The following systems were evaluated: coffee intercropped with one (CJ1), two (CJ2) or three (CJ3) pigeon pea (Cajanus cajan) alleys; coffee planted under full sun (CS); area planted with sweet pepper and snap bean in a conventional tillage system (AC); and secondary forest area (FFR). Row spacing in CJ1, CJ2, CJ3 and CS was 2.0x1.0, 2.8x1.0, 3.6x1.0, and 2.8x1.0 m, respectively. Soil samples were collected at 10-cm depth, during the four seasons of the year. The results were subjected to analysis of variance, principal component analysis, and redundancy analysis. There was an increase in edaphic macrofauna, soil basal respiration, and microbial quotient in the summer. Total macrofauna density was greater in CJ2 followed by CJ3, CS, CJ1, AC and FFR; Coleoptera, Formicidae, and Isoptera were the most abundant groups. There are no significant differences among the areas for soil basal respiration, and the metabolic quotient is higher in CJ1, CJ3, and FFR. Microbial biomass carbon and the contents of K, pH, Ca+Mg, and P show greater values in AC.

Correlation of physical and chemical attributes of soil with sugarcane yield

The objective of this work was to evaluate the correlation between sugarcane yield and some physical and chemical attributes of soil. For this, a 42‑ha test area in Araras, SP, Brazil, was used. Soil properties were determined from samples collected at the beginning of the 2003/2004 harvest season, using a regular 100x100 m grid. Yield assessment was done with a yield monitor (Simprocana). Correlation analyses were performed between sugarcane yield and the following soil properties: pH, pH CaCl2, N, C, cone index, clay content, soil organic matter, P, K, Ca, Mg, H+AL, cation exchange capacity, and base saturation. Correlation coefficients were respectively ‑0.05, ‑0.29, 0.33, 0.41, ‑0.27, 0.22, 0.44, ‑0.24, trace, ‑0.06, 0.01, 0.32, 0.14, and 0.04. Correlations of chemical and physical attributes of soil with sugarcane yield are weak, and, per se, they are not able to explain sugarcane yield variation, which suggests that other variables, besides soil attributes, should be analysed.

Prediction of soil orders with high spatial resolution: response of different classifiers to sampling density

The objective of this work was to evaluate sampling density on the prediction accuracy of soil orders, with high spatial resolution, in a viticultural zone of Serra Gaúcha, Southern Brazil. A digital elevation model (DEM), a cartographic base, a conventional soil map, and the Idrisi software were used. Seven predictor variables were calculated and read along with soil classes in randomly distributed points, with sampling densities of 0.5, 1, 1.5, 2, and 4 points per hectare. Data were used to train a decision tree (Gini) and three artificial neural networks: adaptive resonance theory, fuzzy ARTMap; self‑organizing map, SOM; and multi‑layer perceptron, MLP. Estimated maps were compared with the conventional soil map to calculate omission and commission errors, overall accuracy, and quantity and allocation disagreement. The decision tree was less sensitive to sampling density and had the highest accuracy and consistence. The SOM was the less sensitive and most consistent network. The MLP had a critical minimum and showed high inconsistency, whereas fuzzy ARTMap was more sensitive and less accurate. Results indicate that sampling densities used in conventional soil surveys can serve as a reference to predict soil orders in Serra Gaúcha.

Prediction of soil shear strength in agricultural and natural environments of the Brazilian Cerrado

The objective of this work was to develop uni- and multivariate models to predict maximum soil shear strength (τmax) under different normal stresses (σn), water contents (U), and soil managements. The study was carried out in a Rhodic Haplustox under Cerrado (control area) and under no-tillage and conventional tillage systems. Undisturbed soil samples were taken in the 0.00-0.05 m layer and subjected to increasing U and σn, in shear strength tests. The uni- and multivariate models - respectively τmax=10(a+bU) and τmax=10(a+bU+cσn) - were significant in all three soil management systems evaluated and they satisfactorily explain the relationship between U, σn, and τmax. The soil under Cerrado has the highest shear strength (τ) estimated with the univariate model, regardless of the soil water content, whereas the soil under conventional tillage shows the highest values with the multivariate model, which were associated to the lowest water contents at the soil consistency limits in this management system.

Distribution of chemical compartments of soil organic matter and c stocks of a cambisol from south Brazil as affected by Pinus afforestation

Distribution and stocks of soil organic matter (SOM) compartments after Pinus monoculture introduction in a native pasture area of a Cambisol, Santa Catarina, Brazil, were investigated. Pinus introduction increased soil acidity, content of exchangeable Al+3 and diminished soil nutrients. Nevertheless, soil C stock increased in all humic fractions of the 0-5 cm layer after Pinus afforestation. In the subsurface, the vegetation change only promoted SOM redistribution from the NaOH-extractable humic substances to a less hydrophobic humin fraction. Under Pinus, soil organo-mineral interactions were relevant up to a 15 cm depth, while in pasture environment, this mechanism occurred mainly in the surface layer.

Spatial variability of soil physical properties in two management systems in sugarcane crop

The aim of this study was to characterize the spatial variability of soil bulk density (Bd), soil moisture content (θ) and total porosity (Tp) in two management systems of sugarcane harvesting, with or without burning, in a Haplustox soil, in the 0-0.20 m layer. The study area is located in Rio Brilhante, state of Mato Grosso do Sul, Brazil, in Eldorado Sugar Mill. The plots have presented 180 m length, and 145.6 m width, totaling 90 points distributed in the form of a grid of nine rows by ten columns, with points spaced 20 m from its neighbor. Soil samples were collected at 0-0.20 m layer in 2007/2008 and 2008/2009 crops. The harvest with burning system had a higher density compared to mechanized harvest, in the two study periods. The moisture content as well as the porosity increased proportionally with the decrease of the density of the harvest burning system compared to the mechanized.