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.

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.

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.

USE OF SCALED SEMIVARIOGRAMS IN THE PLANNING SAMPLE OF SOIL CHEMICAL PROPERTIES IN SOUTHERN AMAZONAS, BRAZIL

The lack of information concerning the variability of soil properties has been a major concern of researchers in the Amazon region. Thus, the aim of this study was to evaluate the spatial variability of soil chemical properties and determine minimal sampling density to characterize the variability of these properties in five environments located in the south of the State of Amazonas, Brazil. The five environments were archaeological dark earth (ADE), forest, pasture land, agroforestry operation, and sugarcane crop. Regular 70 × 70 m mesh grids were set up in these areas, with 64 sample points spaced at 10 m distance. Soil samples were collected at the 0.0-0.1 m depth. The chemical properties of pH in water, OM, P, K, Ca, Mg, H+Al, SB, CEC, and V were determined at these points. Data were analyzed by descriptive and geostatistical analyses. A large part of the data analyzed showed spatial dependence. Chemical properties were best fitted to the spherical model in almost all the environments evaluated, except for the sugarcane field with a better fit to the exponential model. ADE and sugarcane areas had greater heterogeneity of soil chemical properties, showing a greater range and higher sampling density; however, forest and agroforestry areas had less variability of chemical properties.

ORGANIC PHOSPHORUS FRACTIONS IN SOIL FERTILIZED WITH CATTLE MANURE

Inorganic phosphorus (Pi) usually controls the P availability in tropical soils, but the contribution of organic P (Po) should not be neglected, mainly in systems with low P input or management systems that promote organic matter accumulation. The aims of this study were to evaluate the changes in the Po fractions over time in soil fertilized and not fertilized with cattle manure and to correlate Po forms with available P extracted by anion exchange resin. The experiment was carried out under field conditions, in a sandy-clay loam Haplustox. The experimental design was a 2 × 9 randomized complete block factorial design, in which the first factor was manure application (20 t ha-1) or absence, and the second the soil sampling times (3, 7, 14, 21, 28, 49, 70, 91, and 112 days) after manure incorporation. Labile, moderately labile and non-labile Po fractions were determined in the soil material of each sampling. Manure fertilization increased the Po levels in the moderately labile and non-labile fractions and the total organic P, but did not affect the Po fraction proportions in relation to total organic P. On average, 5.1 % of total Po was in the labile, 44.4 % in the moderately labile and 50.5 % in the non-labile fractions. Available P (resin P) was more affected by the manure soluble Pi rather than by the labile Po forms. The labile and non-labile Po fractions varied randomly with no defined trend in relation to the samplings; for this reason, the data did not fit any mathematical model.

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.

APPRAISAL OF THE SNAP MODEL FOR PREDICTING NITROGEN MINERALIZATION IN TROPICAL SOILS UNDER EUCALYPTUS

The Soil Nitrogen Availability Predictor (SNAP) model predicts daily and annual rates of net N mineralization (NNM) based on daily weather measurements, daily predictions of soil water and soil temperature, and on temperature and moisture modifiers obtained during aerobic incubation (basal rate). The model was based on in situ measurements of NNM in Australian soils under temperate climate. The purpose of this study was to assess this model for use in tropical soils under eucalyptus plantations in São Paulo State, Brazil. Based on field incubations for one month in three, NNM rates were measured at 11 sites (0-20 cm layer) for 21 months. The basal rate was determined in in situ incubations during moist and warm periods (January to March). Annual rates of 150-350 kg ha-1 yr-1 NNM predicted by the SNAP model were reasonably accurate (R2 = 0.84). In other periods, at lower moisture and temperature, NNM rates were overestimated. Therefore, if used carefully, the model can provide adequate predictions of annual NNM and may be useful in practical applications. For NNM predictions for shorter periods than a year or under suboptimal incubation conditions, the temperature and moisture modifiers need to be recalibrated for tropical conditions.

Contribution of Topography and Incident Solar Radiation to Variation of Soil and Plant Litter at an Area with Heterogeneous Terrain

Natural processes that determine soil and plant litter properties are controlled by multiple factors. However, little attention has been given to distinguishing the effects of environmental factors from the effects of spatial structure of the area on the distribution of soil and litter properties in tropical ecosystems covering heterogeneous topographies. The aim of this study was to assess patterns of soil and litter variation in a tropical area that intercepts different levels of solar radiation throughout the year since its topography has slopes predominantly facing opposing geographic directions. Soil data (pH, C, N, P, H+Al, Ca, Mg, K, Al, Na, sand, and silt) and plant litter data (N, K, Ca, P, and Mg) were gathered together with the geographic coordinates (to model the spatial structure) of 40 sampling units established at two sites composed of slopes predominantly facing northwest and southeast (20 units each). Soil and litter chemical properties varied more among slopes within similar geographic orientations than between the slopes facing opposing directions. Both the incident solar radiation and the spatial structure of the area were relevant in explaining the patterns detected in variation of soil and plant litter. Individual contributions of incident solar radiation to explain the variation in the properties evaluated suggested that this and other environmental factors may play a particularly relevant role in determining soil and plant litter distribution in tropical areas with heterogeneous topography. Furthermore, this study corroborates that the spatial structure of the area also plays an important role in the distribution of soil and litter within this type of landscape, which appears to be consistent with the action of water movement mechanisms in such areas.

MODELING OF SOIL LOAD-BEARING CAPACITY AS A FUNCTION OF SOIL MECHANICAL RESISTANCE TO PENETRATION

Estimation of soil load-bearing capacity from mathematical models that relate preconsolidation pressure (σp) to mechanical resistance to penetration (PR) and gravimetric soil water content (U) is important for defining strategies to prevent compaction of agricultural soils. Our objective was therefore to model the σp and compression index (CI) according to the PR (with an impact penetrometer in the field and a static penetrometer inserted at a constant rate in the laboratory) and U in a Rhodic Eutrudox. The experiment consisted of six treatments: no-tillage system (NT); NT with chiseling; and NT with additional compaction by combine traffic (passing 4, 8, 10, and 20 times). Soil bulk density, total porosity, PR (in field and laboratory measurements), U, σp, and CI values were determined in the 5.5-10.5 cm and 13.5-18.5 cm layers. Preconsolidation pressure (σp) and CI were modeled according to PR in different U. The σp increased and the CI decreased linearly with increases in the PR values. The correlations between σp and PR and PR and CI are influenced by U. From these correlations, the soil load-bearing capacity and compaction susceptibility can be estimated by PR readings evaluated in different U.

LOSS OF EXTRACTION CAPACITY OF MEHLICH-1 AND MONOCALCIUM PHOSPHATE AS A VARIABLE OF REMAINING P AND ITS RELATIONSHIP TO CRITICAL LEVELS OF SOIL PHOSPHORUS AND SULFUR

The Mehlich-1 (M-1) extractant and Monocalcium Phosphate in acetic acid (MCPa) have mechanisms for extraction of available P and S in acidity and in ligand exchange, whether of the sulfate of the extractant by the phosphate of the soil, or of the phosphate of the extractant by the sulfate of the soil. In clayey soils, with greater P adsorption capacity, or lower remaining P (Rem-P) value, which corresponds to soils with greater Phosphate Buffer Capacity (PBC), more buffered for acidity, the initially low pH of the extractants increases over their time of contact with the soil in the direction of the pH of the soil; and the sulfate of the M-1 or the phosphate of the MCPa is adsorbed by adsorption sites occupied by these anions or not. This situation makes the extractant lose its extraction capacity, a phenomenon known as loss of extraction capacity or consumption of the extractant, the object of this study. Twenty soil samples were chosen so as to cover the range of Rem-P (0 to 60 mg L-1). Rem-P was used as a measure of the PBC. The P and S contents available from the soil samples through M-1 and MCPa, and the contents of other nutrients and of organic matter were determined. For determination of loss of extraction capacity, after the rest period, the pH and the P and S contents were measured in both the extracts-soils. Although significant, the loss of extraction capacity of the acidity of the M-1 and MCPa extractants with reduction in the Rem-P value did not have a very expressive effect. A “linear plateau” model was observed for the M-1 for discontinuous loss of extraction capacity of the P content in accordance with reduction in the concentration of the Rem-P or increase in the PBC, suggesting that a discontinuous model should also be adopted for interpretation of available P of soils with different Rem-P values. In contrast, a continuous linear response was observed between the P variables in the extract-soil and Rem-P for the MCPa extractor, which shows increasing loss of extraction capacity of this extractor with an increase in the PBC of the soil, indicating the validity of the linear relationship between the available S of the soil and the PBC, estimated by Rem-P, as currently adopted.

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.

Estimate of Groundwater Recharge Based on Water Balance in The Unsaturated Soil Zone

ABSTRACT Groundwater management depends on the knowledge on recharge rates and water fluxes within aquifers. The recharge is one of the water cycle components most difficult to estimate. As a result, despite the chosen method, the estimates are subject to uncertainties that can be identified by means of comparison with other approaches. In this study, groundwater recharge estimates based on the water balance in the unsaturated zone is assessed. Firstly, the approach is evaluated by comparing the results with those of another method. Then, the estimates are used as inputs in a transient groundwater flow model in order to assess how the water table would respond to the obtained recharges rates compared to measured levels. The results suggest a good performance of the adopted approach and, despite some inherent limitations, it has advantages over other methods since the data required are easier to obtain.

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.

NUTRIENT RECOMMENDATION MODEL FOR CARROT CROP – FERTICALC CARROT

ABSTRACT The literature on fertilization for carrot growing usually recommends nutrient application rates for yield expectations lower than the yields currently obtained. Moreover, the recommendation only considers the results of soil chemical analysis and does not include effects such as crop residues or variations in yield levels. The aim of this study was to propose a fertilizer recommendation system for carrot cultivation (FERTICALC Carrot) which includes consideration of the nutrient supply by crop residues, variation in intended yield, soil chemical properties, and the growing season (winter or summer). To obtain the data necessary for modeling nutritional requirements, 210 carrot production stands were sampled in the region of Alto Paranaíba, State of Minas Gerais, Brazil. The dry matter content of the roots, the coefficient of biological utilization of nutrients in the roots, and the nutrient harvest index for summer and winter crops were determined for these samples. To model the nutrient supply by the soil, the literature was surveyed in regard to this theme. A modeling system was developed for recommendation of macronutrients and B. For cationic micronutrients, the system only reports crop nutrient export and extraction. The FERTICALC which was developed proved to be efficient for fertilizer recommendation for carrot cultivation. Advantages in relation to official fertilizer recommendation tables are continuous variation of nutrient application rates in accordance with soil properties and in accordance with data regarding the extraction efficiency of modern, higher yielding cultivars.

Potential of Spectroradiometry to Classify Soil Clay Content

ABSTRACT Diffuse reflectance spectroscopy (DRS) is a fast and cheap alternative for soil clay, but needs further investigation to assess the scope of application. The purpose of the study was to develop a linear regression model to predict clay content from DRS data, to classify the soils into three textural classes, similar to those defined by a regulation of the Brazilian Ministry of Agriculture, Livestock and Food Supply. The DRS data of 412 soil samples, from the 0.0-0.5 m layer, from different locations in the state of Rio Grande do Sul, Brazil, were measured at wavelengths of 350 to 2,500 nm in the laboratory. The fitting of the linear regression model developed to predict soil clay content from the DRS data was based on a R2 value of 0.74 and 0.75, with a RMSE of 7.82 and 8.51 % for the calibration and validation sets, respectively. Soil texture classification had an overall accuracy of 79.0 % (calibration) and 80.9 % (validation). The heterogeneity of soil samples affected the performance of the prediction models. Future studies should consider a previous classification of soil samples in different groups by soil type, parent material and/or sampling region.