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.

N2O emissions from a cultivated mollisol: optimal time of day for sampling and the role of soil temperature

The correct use of closed field chambers to determine N2O emissions requires defining the time of day that best represents the daily mean N2O flux. A short-term field experiment was carried out on a Mollisol soil, on which annual crops were grown under no-till management in the Pampa Ondulada of Argentina. The N2O emission rates were measured every 3 h for three consecutive days. Fluxes ranged from 62.58 to 145.99 ∝g N-N2O m-2 h-1 (average of five field chambers) and were negatively related (R² = 0.34, p < 0.01) to topsoil temperature (14 - 20 ºC). N2O emission rates measured between 9:00 and 12:00 am presented a high relationship to daily mean N2O flux (R² = 0.87, p < 0.01), showing that, in the study region, sampling in the mornings is preferable for GHG.

Soil thermal diffusivity estimated from data of soil temperature and single soil component properties

Under field conditions, thermal diffusivity can be estimated from soil temperature data but also from the properties of soil components together with their spatial organization. We aimed to determine soil thermal diffusivity from half-hourly temperature measurements in a Rhodic Kanhapludalf, using three calculation procedures (the amplitude ratio, phase lag and Seemann procedures), as well as from soil component properties, for a comparison of procedures and methods. To determine thermal conductivity for short wave periods (one day), the phase lag method was more reliable than the amplitude ratio or the Seemann method, especially in deeper layers, where temperature variations are small. The phase lag method resulted in coherent values of thermal diffusivity. The method using properties of single soil components with the values of thermal conductivity for sandstone and kaolinite resulted in thermal diffusivity values of the same order. In the observed water content range (0.26-0.34 m³ m-3), the average thermal diffusivity was 0.034 m² d-1 in the top layer (0.05-0.15 m) and 0.027 m² d-1 in the subsurface layer (0.15-0.30 m).

Building predictive models of soil particle-size distribution

Is it possible to build predictive models (PMs) of soil particle-size distribution (psd) in a region with complex geology and a young and unstable land-surface? The main objective of this study was to answer this question. A set of 339 soil samples from a small slope catchment in Southern Brazil was used to build PMs of psd in the surface soil layer. Multiple linear regression models were constructed using terrain attributes (elevation, slope, catchment area, convergence index, and topographic wetness index). The PMs explained more than half of the data variance. This performance is similar to (or even better than) that of the conventional soil mapping approach. For some size fractions, the PM performance can reach 70 %. Largest uncertainties were observed in geologically more complex areas. Therefore, significant improvements in the predictions can only be achieved if accurate geological data is made available. Meanwhile, PMs built on terrain attributes are efficient in predicting the particle-size distribution (psd) of soils in regions of complex geology.

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.

Spectroscopic quantification of soil phosphorus forms by 31p-nmr after nine years of organic or mineral fertilization

Long-standing applications of mineral fertilizers or types of organic wastes such as manure can cause phosphorus (P) accumulation and changes in the accumulated P forms in the soil. The objective of this research was to evaluate the forms of P accumulated in soils treated with mineral fertilizer or different types of manure in a long-term experiment. Soil was sampled from the 0-5 cm layer of plots fertilized with five different nutrient sources for nine years: 1) control without fertilizer; 2) mineral fertilizer at recommended rates for local conditions; 3) 5 t ha-1 year-1 of moist poultry litter; 4) 60 m³ ha-1 year-1 of liquid cattle manure and 5) 40 m³ ha-1 year-1 of liquid swine manure. The 31P-NMR spectra of soil extracts detected the following P compounds: orthophosphate, pyrophosphate, inositol phosphate, glycerophosphate, and DNA. The use of organic or mineral fertilizer over nine years did not change the soil P forms but influenced their concentration. Fertilization with mineral or organic fertilizers stimulated P accumulation in inorganic forms. Highest inositol phosphate levels were observed after fertilization with any kind of manure and highest organic P concentration in glycerophosphate form in after mineral or no fertilization.

Residual effect of soil tillage on water erosion from a Typic Paleudalf under long-term no-tillage and cropping systems

Soil erosion is one of the chief causes of agricultural land degradation. Practices of conservation agriculture, such as no-tillage and cover crops, are the key strategies of soil erosion control. In a long-term experiment on a Typic Paleudalf, we evaluated the temporal changes of soil loss and water runoff rates promoted by the transition from conventional to no-tillage systems in the treatments: bare soil (BS); grassland (GL); winter fallow (WF); intercrop maize and velvet bean (M+VB); intercrop maize and jack bean (M+JB); forage radish as winter cover crop (FR); and winter cover crop consortium ryegrass - common vetch (RG+CV). Intensive soil tillage induced higher soil losses and water runoff rates; these effects persisted for up to three years after the adoption of no-tillage. The planting of cover crops resulted in a faster decrease of soil and water loss rates in the first years after conversion from conventional to no-tillage than to winter fallow. The association of no-tillage with cover crops promoted progressive soil stabilization; after three years, soil losses were similar and water runoff was lower than from grassland soil. In the treatments of cropping systems with cover crops, soil losses were reduced by 99.7 and 66.7 %, compared to bare soil and winter fallow, while the water losses were reduced by 96.8 and 71.8 % in relation to the same treatments, respectively.

Critical limits of soil penetration resistance in a rhodic Eutrudox

Soil penetration resistance is an important indicator of soil physical quality and the critical limit of 2 MPa has been widely used to characterize the soil physical quality, in both no-tillage and conventional systems. The aim of this study was to quantify the influence of different tillage and cropping systems on the soil penetration resistance in a Rhodic Eutrudox. The experiment was carried out in a 5 × 2 factorial, completely randomized block design (tillage systems vs cropping systems), with four replications. The tillage systems consisted of: conventional tillage disk harrow; minimum tillage with annual chiseling; minimum tillage with chiseling every three years; no-tillage for 11 consecutive years; and no-tillage for 24 consecutive years. The factor cropping systems was represented by: crop rotation and crop succession. The soil penetration resistance (SPR) was determined in 20 soil samples per treatment and layer (0.0-0.10; 0.10-0.20 and 0.20-0.30 m) for each soil matric potential: -6, -10, -33, -100, -500 kPa. The SPR was determined at a volumetric soil water content equivalent to the fraction of plant-available water of 0.7. There were no differences of soil penetration resistance between the two cropping systems. Differences in soil penetration resistance among tillage systems were related to the matric potential at which the samples were equilibrated. The critical SPR limit of 2 MPa normally used for conventional tillage should be maintained. However, this value of 2 MPa is inappropriate for the physical quality characterization of Rhodic Eutrudox under no-tillage and/or minimum tillage with chiseling. Regardless of the cropping systems, the critical SPR limit should be raised to 3 MPa for minimum tillage with chiseling and to 3.5 MPa for no-tillage.

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

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

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.

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

Comparing the artificial neural network with parcial least squares for prediction of soil organic carbon and pH at different moisture content levels using visible and near-infrared spectroscopy

Visible and near infrared (vis-NIR) spectroscopy is widely used to detect soil properties. The objective of this study is to evaluate the combined effect of moisture content (MC) and the modeling algorithm on prediction of soil organic carbon (SOC) and pH. Partial least squares (PLS) and the Artificial neural network (ANN) for modeling of SOC and pH at different MC levels were compared in terms of efficiency in prediction of regression. A total of 270 soil samples were used. Before spectral measurement, dry soil samples were weighed to determine the amount of water to be added by weight to achieve the specified gravimetric MC levels of 5, 10, 15, 20, and 25 %. A fiber-optic vis-NIR spectrophotometer (350-2500 nm) was used to measure spectra of soil samples in the diffuse reflectance mode. Spectra preprocessing and PLS regression were carried using Unscrambler® software. Statistica® software was used for ANN modeling. The best prediction result for SOC was obtained using the ANN (RMSEP = 0.82 % and RPD = 4.23) for soil samples with 25 % MC. The best prediction results for pH were obtained with PLS for dry soil samples (RMSEP = 0.65 % and RPD = 1.68) and soil samples with 10 % MC (RMSEP = 0.61 % and RPD = 1.71). Whereas the ANN showed better performance for SOC prediction at all MC levels, PLS showed better predictive accuracy of pH at all MC levels except for 25 % MC. Therefore, based on the data set used in the current study, the ANN is recommended for the analyses of SOC at all MC levels, whereas PLS is recommended for the analysis of pH at MC levels below 20 %.

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.