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.

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.

Effects of surface application of calcium-magnesium silicate and gypsum on soil fertility and sugarcane yield

Lime application recommendations for amendment of soil acidity in sugarcane were developed with a burnt cane harvesting system in mind. Sugarcane is now harvested in most areas without burning, and lime application for amendment of soil acidity in this system in which the sugarcane crop residue remains on the ground has been carried out without a scientific basis. The aim of this study was to evaluate the changes in soil acidity and stalk and sugar yield with different rates of surface application of calcium, magnesium silicate, and gypsum in ratoon cane. The experiment was performed after the 3rd harvest of the variety SP 81-3250 in a commercial green sugarcane plantation of the São Luiz Sugar Mill (47º 25' 33" W; 21º 59' 46" S), located in Pirassununga, São Paulo, in southeast Brazil. A factorial arrangement of four Ca-Mg silicate rates (0, 850, 1700, and 3400 kg ha-1) and two gypsum rates (0 and 1700 kg ha-1) was used in the experiment. After 12 months, the experiment was harvested and technological measurements of stalk and sugar yield were made. After harvest, soil samples were taken at the depths of 0.00-0.05, 0.05-0.10, 0.10-0.20, 0.20-0.40, and 0.40-0.60 m in all plots, and the following determinations were made: soil pH in CaCl2, organic matter, P, S, K, Ca, Mg, H+Al, Al, Si, and base saturation. The results show that the application of gypsum reduced the exchangeable Al3+ content and Al saturation below 0.05 m, and increased the Ca2+ concentration in the whole profile, the Mg2+ content below 0.10 m, K+ below 0.4 m, and base saturation below 0.20 m. This contributed to the effect of surface application of silicate on amendment of soil acidity reaching deeper layers. From the results of this study, it may be concluded that the silicate rate recommended may be too low, since the greater rates used in this experiment showed greater reduction in soil acidity, higher levels of nutrients at greater depths and an increase in stalk and sugar yield.

CHARACTERIZATION OF BULK SOIL HUMIN AND ITS ALKALINE-SOLUBLE AND ALKALINE-INSOLUBLE FRACTIONS

Humic substances are the major components of soil organic matter. Among the three humic substance components (humic acid, fulvic acid, and humin), humin is the most insoluble in aqueous solution at any pH value and, in turn, the least understood. Humin has poor solubility mainly because it is tightly bonded to inorganic soil colloids. By breaking the linkage between humin and inorganic soil colloids using inorganic or organic solvents, bulk humin can be partially soluble in alkali, enabling a better understanding of the structure and properties of humin. However, the structural relationship between bulk humin and its alkaline-soluble (AS) and alkaline-insoluble (AIS) fractions is still unknown. In this study, we isolated bulk humin from two soils of Northeast China by exhaustive extraction (25 to 28 times) with 0.1 mol L-1 NaOH + 0.1 mol L-1 Na4P2O7, followed by the traditional treatment with 10 % HF-HCl. The isolated bulk humin was then fractionated into AS-humin and AIS-humin by exhaustive extraction (12 to 15 times) with 0.1 mol L-1 NaOH. Elemental analysis and solid-state 13C cross-polarization magic angle spinning nuclear magnetic resonance (13C CPMAS NMR) spectroscopy were used to characterize and compare the chemical structures of bulk humin and its corresponding fractions. The results showed that, regardless of soil types, bulk humin was the most aliphatic and most hydrophobic, AS-humin was the least aliphatic, and AIS-humin was the least alkylated among the three humic components. The results showed that bulk humin and its corresponding AS-humin and AIS-humin fractions are structurally differed from one another, implying that the functions of these humic components in the soil environment differed.

EFFECT OF ALTERNATIVE MULTINUTRIENT SOURCES ON SOIL CHEMICAL PROPERTIES

The current high price of potassium chloride and the dependence of Brazil on imported materials to supply the domestic demand call for studies evaluating the efficiency of alternative sources of nutrients. The aim of this work was to evaluate the effect of silicate rock powder and a manganese mining by-product, and secondary materials originated from these two materials, on soil chemical properties and on brachiaria production. This greenhouse experiment was conducted in pots with 5 kg of soil (Latossolo Vermelho-Amarelo distrófico - Oxisol). The alternative nutrient sources were: verdete, verdete treated with NH4OH, phonolite, ultramafic rock, mining waste and the proportion of 75 % of these K fertilizers and 25 % lime. Mixtures containing 25 % of lime were heated at 800 ºC for 1 h. These sources were applied at rates of 0, 150, 300, 450 and 600 kg ha-1 K2O, and incubated for 45 days. The mixtures of heated silicate rocks with lime promoted higher increases in soil pH in decreasing order: ultramafic rock>verdete>phonolite>mining waste. Applying the mining waste-lime mixture increased soil exchangeable K, and available P when ultramafic rock was incorporated. When ultramafic rock was applied, the release of Ca2+ increased significantly. Mining subproduct released the highest amount of Zn2+ and Mn2+ to the soil. The application of alternative sources of K, with variable chemical composition, altered the nutrient availability and soil chemical properties, improving mainly plant development and K plant uptake, and are important nutrient sources.

DISTRIBUTION OF ORGANIC CARBON IN DIFFERENT SOIL FRACTIONS IN ECOSYSTEMS OF CENTRAL AMAZONIA

Organic matter plays an important role in many soil properties, and for that reason it is necessary to identify management systems which maintain or increase its concentrations. The aim of the present study was to determine the quality and quantity of organic C in different compartments of the soil fraction in different Amazonian ecosystems. The soil organic matter (FSOM) was fractionated and soil C stocks were estimated in primary forest (PF), pasture (P), secondary succession (SS) and an agroforestry system (AFS). Samples were collected at the depths 0-5, 5-10, 10-20, 20-40, 40-60, 60-80, 80-100, 100-160, and 160-200 cm. Densimetric and particle size analysis methods were used for FSOM, obtaining the following fractions: FLF (free light fraction), IALF (intra-aggregate light fraction), F-sand (sand fraction), F-clay (clay fraction) and F-silt (silt fraction). The 0-5 cm layer contains 60 % of soil C, which is associated with the FLF. The F-clay was responsible for 70 % of C retained in the 0-200 cm depth. There was a 12.7 g kg-1 C gain in the FLF from PF to SS, and a 4.4 g kg-1 C gain from PF to AFS, showing that SS and AFS areas recover soil organic C, constituting feasible C-recovery alternatives for degraded and intensively farmed soils in Amazonia. The greatest total stocks of carbon in soil fractions were, in decreasing order: (101.3 Mg ha-1 of C - AFS) > (98.4 Mg ha-1 of C - FP) > (92.9 Mg ha-1 of C - SS) > (64.0 Mg ha-1 of C - P). The forms of land use in the Amazon influence C distribution in soil fractions, resulting in short- or long-term changes.

EFFECT OF WATER AVAILABILITY ON SOIL MICROBIAL BIOMASS IN SECONDARY FOREST IN EASTERN AMAZONIA

Soil microbial biomass (SMB) plays an important role in nutrient cycling in agroecosystems, and is limited by several factors, such as soil water availability. This study assessed the effects of soil water availability on microbial biomass and its variation over time in the Latossolo Amarelo concrecionário of a secondary forest in eastern Amazonia. The fumigation-extraction method was used to estimate the soil microbial biomass carbon and nitrogen content (SMBC and SMBN). An adaptation of the fumigation-incubation method was used to determine basal respiration (CO2-SMB). The metabolic quotient (qCO2) and ratio of microbial carbon:organic carbon (CMIC:CORG) were calculated based on those results. Soil moisture was generally significantly lower during the dry season and in the control plots. Irrigation raised soil moisture to levels close to those observed during the rainy season, but had no significant effect on SMB. The variables did not vary on a seasonal basis, except for the microbial C/N ratio that suggested the occurrence of seasonal shifts in the structure of the microbial community.

EFFECTS OF COMPACTION DUE TO MACHINERY TRAFFIC ON SOIL PORE CONFIGURATION

Soil compaction has been recognized as a severe problem in mechanized agriculture and has an influence on many soil properties and processes. Yet, there are few studies on the long-term effects of soil compaction, and the development of soil compaction has been shown through a limited number of soil parameters. The objectives of this study were to evaluate the persistence of soil compaction effects (three traffic treatments: T0, without traffic; T3, three tractor passes; and T5, five tractor passes) on pore system configuration, through static and dynamic determinations; and to determine changes in soil pore orientation due to soil compaction through measurement of hydraulic conductivity of saturated soil in samples taken vertically and horizontally. Traffic led to persistent changes in all the dynamic indicators studied (saturated hydraulic conductivity, K0; effective macro- and mesoporosity, εma and εme), with significantly lower values of K0, εma, and εme in the T5 treatment. The static indicators of bulk density (BD), derived total porosity (TP), and total macroporosity (θma) did not vary significantly among the treatments. This means that machine traffic did not produce persistent changes on these variables after two years. However, the orientation of the soil pore system was modified by traffic. Even in T0, there were greater changes in K0 measured in the samples taken vertically than horizontally, which was more related to the presence of vertical biopores, and to isotropy of K0 in the treatments with machine traffic. Overall, the results showed that dynamic indicators are more sensitive to the effects of compaction and that, in the future, static indicators should not be used as compaction indicators without being complemented by dynamic indicators.

Boundary between Soil and Saprolite in Alisols in the South of Brazil

Despite numerous studies conducted on the lower limit of soil and its contact with saprolite layers, a great deal of work is left to standardize identification and annotation of these variables in the field. In shallow soils, the appropriately noting these limits or contacts is essential for determining their behavior and potential use. The aims of this study were to identify and define the field contact and/or transition zone between soil and saprolite in profiles of an Alisol derived from fine sandstone and siltstone/claystone in subtropical southern Brazil and to subsequently validate the field observations through a multivariate analysis of laboratory analytical data. In the six Alisol profiles evaluated, the sequence of horizons found was A, Bt, C, and Cr, where C was considered part of the soil due to its pedogenetic structure, and Cr was considered saprolite due to its rock structure. The morphological properties that were determined in the field and that were different between the B and C horizons and the Cr layer were color, structure, texture, and fragments of saprolite. According to the test of means, the properties that support the inclusion of the C horizon as part of the soil are sand, clay, water-dispersible clay, silt/clay ratio, macroporosity, total porosity, resistance to penetration, cation exchange capacity, Fe extracted by DCB, Al, H+Al, and cation exchange capacity of clay. The properties that support the C horizon as a transition zone are silt, Ca, total organic C, and Fe extracted by ammonium oxalate. Discriminant analysis indicated differences among the three horizons evaluated.

Changes in the Structure of a Nigerian Soil under Different Land Management Practices

Quantification of soil physical quality (SPQ) and pore size distribution (PSD) can assist understanding of how changes in land management practices influence dynamics of soil structure, and this understanding could greatly improve the predictability of soil physical behavior and crop yield. The objectives of this study were to measure the SPQ index under two different land management practices (the continuous arable cropping system and natural bush fallow system), and contrast the effects of these practices on the structure of PSD using soil water retention data. Soil water retention curves obtained from a pressure chamber were fitted to van Genuchten’s equation, setting m (= 1-1/n). Although values for soil bulk density were high, soils under the continuous arable cropping system had good SPQ, and maintained the capacity to support root development. However, soils under the natural bush fallow system had a worse structure than the continuous arable system, with restrictions in available water capacity. These two management systems had different PSDs. Results showed the inferiority of the natural bush fallow system with no traffic restriction (which is the common practice) in relation to the continuous arable cropping system in regard to physical quality and structure.

Organic Matter Fractions and Quality of the Surface Layer of a Constructed and Vegetated Soil After Coal Mining. II - Physical Compartments and Carbon Management Index

Soils constructed after mining often have low carbon (C) stocks and low quality of organic matter (OM). Cover crops are decisive for the recovery process of these stocks, improving the quality of constructed soils. Therefore, the goal of this study was to evaluate the effect of cover crops on total organic C (TOC) stocks, C distribution in physical fractions of OM and the C management index (CMI) of a soil constructed after coal mining. The experiment was initiated in 2003 with six treatments: Hemarthria altissima (T1), Paspalum notatum (T2), Cynodon dactylon (T3), Urochloa brizantha (T4), bare constructed soil (T5), and natural soil (T6). Soil samples were collected in 2009 from the 0.00-0.03 m layer, and the TOC and C stocks in the physical particle size fractions (carbon in the coarse fraction - CCF, and mineral-associated carbon - MAC) and density fractions (free light fraction - FLF; occluded light fraction - OLF, and heavy fraction - HF) of OM were determined. The CMI components: carbon pool index (CPI), lability (L) and lability index (LI) were estimated by both fractionation methods. No differences were observed between TOC, CCF and MAC stocks. The lowest C stocks in FLF and OLF fractions were presented by T2, 0.86 and 0.61 Mg ha-1, respectively. The values of TOC stock, C stock in physical fractions and CMI were intermediate, greater than T5 and lower than T6 in all treatments, indicating the partial recovery of soil quality. As a result of the better adaptation of the species Hemarthria and Brizantha, resulting in greater accumulation of labile organic material, the CPI, L, LI and CMI values were higher in these treatments, suggesting a greater potential of these species for recovery of constructed soils.

SOIL PHOSPHORUS THRESHOLDS IN EVALUATING RISK OF ENVIRONMENTAL TRANSFER TO SURFACE WATERS IN SANTA CATARINA, BRAZIL

The State of Santa Catarina, Brazil, has agricultural and livestock activities, such as pig farming, that are responsible for adding large amounts of phosphorus (P) to soils. However, a method is required to evaluate the environmental risk of these high soil P levels. One possible method for evaluating the environmental risk of P fertilization, whether organic or mineral, is to establish threshold levels of soil available P, measured by Mehlich-1 extractions, below which there is not a high risk of P transfer from the soil to surface waters. However, the Mehlich-1 extractant is sensitive to soil clay content, and that factor should be considered when establishing such P-thresholds. The objective of this study was to determine P-thresholds using the Mehlich-1 extractant for soils with different clay contents in the State of Santa Catarina, Brazil. Soil from the B-horizon of an Oxisol with 800 g kg-1 clay was mixed with different amounts of sand to prepare artificial soils with 200, 400, 600, and 800 g kg-1 clay. The artificial soils were incubated for 30 days with moisture content at 80 % of field capacity to stabilize their physicochemical properties, followed by additional incubation for 30 days after liming to raise the pH(H2O) to 6.0. Soil P sorption curves were produced, and the maximum sorption (Pmax) was determined using the Langmuir model for each soil texture evaluated. Based on the Pmax values, seven rates of P were added to four replicates of each soil, and incubated for 20 days more. Following incubation, available P contents (P-Mehlich-1) and P dissolved in the soil solution (P-water) were determined. A change-point value (the P-Mehlich-1 value above which P-water starts increasing sharply) was calculated through the use of segmented equations. The maximum level of P that a soil might safely adsorb (P-threshold) was defined as 80 % of the change-point value to maintain a margin for environmental safety. The P-threshold value, in mg dm-3, was dependent on the soil clay content according to the model P-threshold = 40 + Clay, where the soil clay content is expressed as a percentage. The model was tested in 82 diverse soil samples from the State of Santa Catarina and was able to distinguish samples with high and low environmental risk.

Filter Paper Method for the Determination of the Soil Water Retention Curve

ABSTRACT High cost and long time required to determine a retention curve by the conventional methods of the Richards Chamber and Haines Funnel limit its use; therefore, alternative methods to facilitate this routine are needed. The filter paper method to determine the soil water retention curve was evaluated and compared to the conventional method. Undisturbed samples were collected from five different soils. Using a Haines Funnel and Richards Chamber, moisture content was obtained for tensions of 2; 4; 6; 8; 10; 33; 100; 300; 700; and 1,500 kPa. In the filter paper test, the soil matric potential was obtained from the filter-paper calibration equation, and the moisture subsequently determined based on the gravimetric difference. The van Genuchten model was fitted to the observed data of soil matric potential versus moisture. Moisture values of the conventional and the filter paper methods, estimated by the van Genuchten model, were compared. The filter paper method, with R2 of 0.99, can be used to determine water retention curves of agricultural soils as an alternative to the conventional method.

Contrasts in Areas of Rubber Tree Clones in Regard to Soil and Biomass Carbon Stocks

ABSTRACT Rubber tree (Hevea brasiliensis) crop may accumulate significant amounts of carbon either in biomass or in the soil. However, a comprehensive understanding of the potential of the C stock among different rubber tree clones is still distant, since clones are typically developed to exhibit other traits, such as better yield and disease tolerance. Thus, the aim of this study was to address differences among different areas planted to rubber clones. We hypothesized that different rubber tree clones, developed to adapt to different environmental and biological constrains, diverge in terms of soil and plant biomass C stocks. Clones were compared in respect to soil C stocks at four soil depths and the total depth (0.00-0.05, 0.05-0.10, 0.10-0.20, 0.20-0.40, and 0.00-0.40 m), and in the different compartments of the tree biomass. Five different plantings of rubber clones (FX3864, FDR 5788, PMB 1, MDX 624, and CDC 312) of seven years of age were compared, which were established in a randomized block design in the experimental field in Rio de Janeiro State. No difference was observed among plantings of rubber tree clones in regard to soil C stocks, even considering the total stock from 0.00-0.40 m depth. However, the rubber tree clones were different from each other in terms of total plant C stocks, and this contrast was predominately due to only one component of the total C stock, tree biomass. For biomass C stock, the MDX 624 rubber tree clone was superior to other clones, and the stem was the biomass component which most accounted for total C biomass. The contrast among rubber clones in terms of C stock is mainly due to the biomass C stock; the aboveground (tree biomass) and the belowground (soil) compartments contributed differently to the total C stock, 36.2 and 63.8 %, respectively. Rubber trees did not differ in relation to C stocks in the soil, but the right choice of a rubber clone is a reliable approach for sequestering C from the air in the biomass of trees.