Fuzzy logic applied to the modeling of water dynamics in an Oxisol in northeastern Brazil

Modeling of water movement in non-saturated soil usually requires a large number of parameters and variables, such as initial soil water content, saturated water content and saturated hydraulic conductivity, which can be assessed relatively easily. Dimensional flow of water in the soil is usually modeled by a nonlinear partial differential equation, known as the Richards equation. Since this equation cannot be solved analytically in certain cases, one way to approach its solution is by numerical algorithms. The success of numerical models in describing the dynamics of water in the soil is closely related to the accuracy with which the water-physical parameters are determined. That has been a big challenge in the use of numerical models because these parameters are generally difficult to determine since they present great spatial variability in the soil. Therefore, it is necessary to develop and use methods that properly incorporate the uncertainties inherent to water displacement in soils. In this paper, a model based on fuzzy logic is used as an alternative to describe water flow in the vadose zone. This fuzzy model was developed to simulate the displacement of water in a non-vegetated crop soil during the period called the emergency phase. The principle of this model consists of a Mamdani fuzzy rule-based system in which the rules are based on the moisture content of adjacent soil layers. The performances of the results modeled by the fuzzy system were evaluated by the evolution of moisture profiles over time as compared to those obtained in the field. The results obtained through use of the fuzzy model provided satisfactory reproduction of soil moisture profiles.

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.

Variability of root traits in common bean genotypes at different levels of phosphorus supply and ontogenetic stages

Selection of common bean (Phaseolus vulgaris L.) cultivars with enhanced root growth would be a strategy for increasing P uptake and grain yield in tropical soils, but the strong plasticity of root traits may compromise their inclusion in breeding programs. The aim of this study was to evaluate the magnitude of the genotypic variability of root traits in common bean plants at two ontogenetic stages and two soil P levels. Twenty-four common bean genotypes, comprising the four growth habits that exist in the species and two wild genotypes, were grown in 4 kg pots at two levels of applied P (20 and 80 mg kg-1) and harvested at the stages of pod setting and early pod filling. Root area and root length were measured by digital image analysis. Significant genotype × P level and genotype × harvest interactions in analysis of variance indicate that the genotypic variation of root traits depended on soil nutrient availability and the stage at which evaluation was made. Genotypes differed for taproot mass, basal and lateral root mass, root area and root length at both P levels and growth stages; differences in specific root area and length were small. Genotypes with growth habits II (upright indeterminate) and III (prostrate indeterminate) showed better adaptation to limited P supply than genotypes of groups I (determinate) and IV (indeterminate climbing). Between the two harvests, genotypes of groups II and III increased the mass of basal and lateral roots by 40 and 50 %, respectively, whereas genotypes of groups I and IV by only 7 and 19 %. Values of the genotypic coefficient of determination, which estimates the proportion of phenotypic variance resulting from genetic effects, were higher at early pod filling than at pod setting. Correlations between shoot mass and root mass, which could indicate indirect selection of root systems via aboveground biomass, were higher at early pod filling than at pod setting. The results indicate that selection for root traits in common bean genotypes should preferentially be performed at the early pod-filling stage.

Spatial estimation of foliar phosphorus in different species of the genus Coffea based on soil properties

Information underlying analyses of coffee fertilization systems should consider both the soil and the nutritional status of plants. This study investigated the spatial relationship between phosphorus (P) levels in coffee plant tissues and soil chemical and physical properties. The study was performed using two arabica and one canephora coffee variety. Sampling grids were established in the areas, and the points georeferenced. The assessed properties of the soil were levels of available phosphorus (P-Mehlich), remaining phosphorus (P-rem) and particle size, and of the plant tissue, phosphorus levels (foliar P). The data were subjected to descriptive statistical analysis, correlation analysis, cluster analysis, and probability tests. Geostatistical and trend analyses were only performed for pairs of variables with significant linear correlation. The spatial variability for foliar P content was high for the variety Catuai and medium for the other evaluated plants. Unlike P-Mehlich, the variability in P-rem of the soil indicated the nutritional status of this nutrient in the plant.

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.

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.

SPATIAL DEPENDENCE OF ELECTRICAL CONDUCTIVITY AND CHEMICAL PROPERTIES OF THE SOIL BY ELECTROMAGNETIC INDUCTION

Brazilian soils have natural high chemical variability; thus, apparent electrical conductivity (ECa) can assist interpretation of crop yield variations. We aimed to select soil chemical properties with the best linear and spatial correlations to explain ECa variation in the soil using a Profiler sensor (EMP-400). The study was carried out in Sidrolândia, MS, Brazil. We analyzed the following variables: electrical conductivity - EC (2, 7, and 15 kHz), organic matter, available K, base saturation, and cation exchange capacity (CEC). Soil ECa was measured with the aid of an all-terrain vehicle, which crossed the entire area in strips spaced at 0.45 m. Soil samples were collected at the 0-20 cm depth with a total of 36 samples within about 70 ha. Classical descriptive analysis was applied to each property via SAS software, and GS+ for spatial dependence analysis. The equipment was able to simultaneously detect ECa at the different frequencies. It was also possible to establish site-specific management zones through analysis of correlation with chemical properties. We observed that CEC was the property that had the best correlation with ECa at 15 kHz.

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.

Growth and yield of common bean cultivars at two soil phosphorus levels under biological nitrogen fixation

The genotypic differences on growth and yield of common bean (Phaseolus vulgaris L.) in response to P supply were evaluated in a field experiment under biological N2 fixation. Eight cultivars were grown at two levels of applied P (12 and 50 kg ha-1 of P -- P1 and P2 respectively), in randomized block design in factorial arrangement. Vegetative biomass was sampled at three ontogenetic stages. The effects of genotype and phosphorus were significant for most traits, but not the genotype ´ phosphorus interaction. The cultivars presented different patterns of biomass production and nutrient accumulation, particularly on root system. At P1, P accumulation persisted after the beginning of pod filling, and P translocation from roots to shoots was lower. The nodule senescence observed after flowering might have reduced N2 fixation during pod filling. The responses of vegetative growth to the higher P supply did not reflect with the same magnitude on yield, which increased only 6% at P2; hence the harvest index was lower at P2. The cultivars with highest yields also presented lower grain P concentrations. A sub-optimal supply of N could have limited the expression of the yield potential of cultivars, reducing the genotypic variability of responses to P levels.

A Study to Correlate Soil Consistency Limits with Soil Moisture Tensions, HR-27, 1965

The liquid and plastic limits of a soil are consistency limits that were arbitrarily chosen by Albert Atterberg in 1911. Their determination is by strictly empirical testing procedures. Except for the development of a liquid limit device and subsequent minor refinements the method has remained basically unchanged for over a half century. The empirical determination of an arbitrary limit would seem to be contrary to the very foundations of scientific procedures. However, the tests are relatively simple and the results are generally acceptable and valuable in almost every conceivable use of soil from an engineering standpoint. Such a great volume of information has been collected and compiled by application of these limits to cohesive soils, that it would be impractical and virtually impossible to replace the tests with a more rational testing method. Nevertheless, many believe that the present method is too time consuming and inconsistent. Research was initiated to investigate the development of a rapid and consistent method by relating the limits to soil moisture tension values determined by porous plate and pressure membrane apparatus. With the moisture tension method, hundreds of samples may be run at one time, operator variability is minimal, results are consistent, and a high degree of correlation to present liquid limit tests is possible.

Inclusion of soil erosion impacts in life cycle assessment on a global scale: application to energy crops in Spain

Purpose: Despite the fundamental role of ecosystem goods and services in sustaining human activities, there is no harmonized and internationally agreed method for including them in life cycle assessment (LCA). The main goal of this study was to develop a globally applicable and spatially resolved method for assessing land-use impacts on the erosion regulation ecosystem service.Methods: Soil erosion depends much on location. Thus, unlike conventional LCA, the endpoint method was regionalized at the grid-cell level (5 arc-minutes, approximately 10×10 km2) to reflect the spatial conditions of the site. Spatially explicit characterization factors were not further aggregated at broader spatial scales. Results and discussion: Life cycle inventory data of topsoil and topsoil organic carbon (SOC) losses were interpreted at the endpoint level in terms of the ultimate damage to soil resources and ecosystem quality. Human health damages were excluded from the assessment. The method was tested on a case study of five three-year agricultural rotations, two of them with energy crops, grown in several locations in Spain. A large variation in soil and SOC losses was recorded in the inventory step, depending on climatic and edaphic conditions. The importance of using a spatially explicit model and characterization factors is shown in the case study.Conclusions and outlook: The regionalized assessment takes into account the differences in soil erosion-related environmental impacts caused by the great variability of soils. Taking this regionalized framework as the starting point, further research should focus on testing the applicability of the method trough the complete life cycle of a product and on determining an appropriate spatial scale at which to aggregate characterization factors, in order to deal with data gaps on location of processes, especially in the background system. Additional research should also focus on improving reliability of the method by quantifying and, insofar as it is possible, reducing uncertainty.

Fly Ash Soil Stabilization for Non-Uniform Subgrade Soils, Volume II: Influence of Subgrade Non-Uniformity on PCC Pavement Performance, TR-461, 2005

To provide insight into subgrade non-uniformity and its effects on pavement performance, this study investigated the influence of non-uniform subgrade support on pavement responses (stress and deflection) that affect pavement performance. Several reconstructed PCC pavement projects in Iowa were studied to document and evaluate the influence of subgrade/subbase non-uniformity on pavement performance. In situ field tests were performed at 12 sites to determine the subgrade/subbase engineering properties and develop a database of engineering parameter values for statistical and numerical analysis. Results of stiffness, moisture and density, strength, and soil classification were used to determine the spatial variability of a given property. Natural subgrade soils, fly ash-stabilized subgrade, reclaimed hydrated fly ash subbase, and granular subbase were studied. The influence of the spatial variability of subgrade/subbase on pavement performance was then evaluated by modeling the elastic properties of the pavement and subgrade using the ISLAB2000 finite element analysis program. A major conclusion from this study is that non-uniform subgrade/subbase stiffness increases localized deflections and causes principal stress concentrations in the pavement, which can lead to fatigue cracking and other types of pavement distresses. Field data show that hydrated fly ash, self-cementing fly ash-stabilized subgrade, and granular subbases exhibit lower variability than natural subgrade soils. Pavement life should be increased through the use of more uniform subgrade support. Subgrade/subbase construction in the future should consider uniformity as a key to long-term pavement performance.

Spatio-temporal variability of sugarcane yield efficiency in the state of São Paulo, Brazil

The objective of this work was to assess the spatial and temporal variability of sugarcane yield efficiency and yield gap in the state of São Paulo, Brazil, throughout 16 growing seasons, considering climate and soil as main effects, and socioeconomic factors as complementary. An empirical model was used to assess potential and attainable yields, using climate data series from 37 weather stations. Soil effects were analyzed using the concept of production environments associated with a soil aptitude map for sugarcane. Crop yield efficiency increased from 0.42 to 0.58 in the analyzed period (1990/1991 to 2005/2006 crop seasons), and yield gap consequently decreased from 58 to 42%. Climatic factors explained 43% of the variability of sugarcane yield efficiency, in the following order of importance: solar radiation, water deficit, maximum air temperature, precipitation, and minimum air temperature. Soil explained 15% of the variability, considering the average of all seasons. There was a change in the correlation pattern of climate and soil with yield efficiency after the 2001/2002 season, probably due to the crop expansion to the west of the state during the subsequent period. Socioeconomic, biotic and crop management factors together explain 42% of sugarcane yield efficiency in the state of São Paulo.

Effects of nitrogen supply on wheat and on soil nitrate

Field experiments were conducted at two locations during two growing seasons in the Ebro Valley (Spain), to evaluate the effects of N fertilization on yield and quality of Mediterranean-type wheat in irrigated conditions. Seven N treatments and a control were investigated. The average grain yields ranged from 2117 to 5551 kg ha-1 depending on the year and location. Grain protein ranged from 14.25 to 16.9%, and other quality parameters such as the dough strength (W) also varied with year and location, confirming the suitability of Mediterranean-type wheat and the climate for the production of good bread-making quality wheat. However, grain yields are normally low and both yields and quality can be greatly affected by the variability of this type of climate, even under irrigation. Under these conditions, grain yield increases were mainly due to an increase in the number of grains per m2 without a reduction in the N content per spike, suggesting that N in the grain was not source-limited, possibly due to the lower grain yields and relatively high soil nitrate concentrations. In soils with lower initial soil NO-3N contents, better grain yields could be achieved by applying a N fertilizer rate of about 100 kg N ha-1, whereas in soils with high initial NO-3N contents, no N or a maximum rate of 50 kg N ha-1 is needed to obtain a good grain quality, showing the possibility of producing high-quality wheat with a low amount of N fertilizer and thus increasing the sustainability of the cropping system.