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.

Plant-available soil water capacity: estimation methods and implications

The plant-available water capacity of the soil is defined as the water content between field capacity and wilting point, and has wide practical application in planning the land use. In a representative profile of the Cerrado Oxisol, methods for estimating the wilting point were studied and compared, using a WP4-T psychrometer and Richards chamber for undisturbed and disturbed samples. In addition, the field capacity was estimated by the water content at 6, 10, 33 kPa and by the inflection point of the water retention curve, calculated by the van Genuchten and cubic polynomial models. We found that the field capacity moisture determined at the inflection point was higher than by the other methods, and that even at the inflection point the estimates differed, according to the model used. By the WP4-T psychrometer, the water content was significantly lower found the estimate of the permanent wilting point. We concluded that the estimation of the available water holding capacity is markedly influenced by the estimation methods, which has to be taken into consideration because of the practical importance of this parameter.

Quality of surface water related to land use: a case study in a catchment with small farms and intensive vegetable crop production in southern Brazil

Water degradation is strongly related to agricultural activity. The aim of this study was to evaluate the influence of land use and some environmental components on surface water quality in the Campestre catchment, located in Colombo, state of Parana, Brazil. Physical and chemical attributes were analyzed (total nitrogen, ammonium, nitrate, total phosphorus, electrical conductivity, pH, temperature, turbidity, total solids, biological oxygen demand, chemical oxygen demand and dissolved oxygen). Monthly samples of the river water were taken over one year at eight monitoring sites, distributed over three sub-basins. Overall, water quality was worse in the sub-basin with a higher percentage of agriculture, and was also affected by a lower percentage of native forest and permanent preservation area, and a larger drainage area. Water quality was also negatively affected by the presence of agriculture in the riparian zone. In the summer season, probably due to higher rainfall and intensive soil use, a higher concentration of total nitrogen and particulate nitrogen was observed, as well as higher electrical conductivity, pH and turbidity. All attributes, except for total phosphorus, were in compliance with Brazilian Conama Resolution Nº 357/2005 for freshwater class 1. However, it should be noted that these results referred to the base flow and did not represent a discharge condition since most of the water samples were not collected at or near the rainfall event.

Pedotransfer functions to predict water retention for soils of the humid tropics: a review

Over the past three decades, pedotransfer functions (PTFs) have been widely used by soil scientists to estimate soils properties in temperate regions in response to the lack of soil data for these regions. Several authors indicated that little effort has been dedicated to the prediction of soil properties in the humid tropics, where the need for soil property information is of even greater priority. The aim of this paper is to provide an up-to-date repository of past and recently published articles as well as papers from proceedings of events dealing with water-retention PTFs for soils of the humid tropics. Of the 35 publications found in the literature on PTFs for prediction of water retention of soils of the humid tropics, 91 % of the PTFs are based on an empirical approach, and only 9 % are based on a semi-physical approach. Of the empirical PTFs, 97 % are continuous, and 3 % (one) is a class PTF; of the empirical PTFs, 97 % are based on multiple linear and polynomial regression of n th order techniques, and 3 % (one) is based on the k-Nearest Neighbor approach; 84 % of the continuous PTFs are point-based, and 16 % are parameter-based; 97 % of the continuous PTFs are equation-based PTFs, and 3 % (one) is based on pattern recognition. Additionally, it was found that 26 % of the tropical water-retention PTFs were developed for soils in Brazil, 26 % for soils in India, 11 % for soils in other countries in America, and 11 % for soils in other countries in Africa.

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.

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

Impacts of deforestation on water balance components of a watershed on the Brazilian East Coast

The Brazilian East coast was intensely affected by deforestation, which drastically cut back the original biome. The possible impacts of this process on water resources are still unknown. The purpose of this study was an evaluation of the impacts of deforestation on the main water balance components of the Galo creek watershed, in the State of Espírito Santo, on the East coast of Brazil. Considering the real conditions of the watershed, the SWAT model was calibrated with data from 1997 to 2000 and validated for the period between 2001 and 2003. The calibration and validation processes were evaluated by the Nash-Sutcliffe efficiency coefficient and by the statistical parameters (determination coefficient, slope coefficient and F test) of the regression model adjusted for estimated and measured flow data. After calibration and validation of the model, new simulations were carried out for three different land use scenarios: a scenario in compliance with the law (C1), assuming the preservation of PPAs (permanent preservation areas); an optimistic scenario (C2), which considers the watershed to be almost entirely covered by native vegetation; and a pessimistic scenario (C3), in which the watershed would be almost entirely covered by pasture. The scenarios C1, C2 and C3 represent a soil cover of native forest of 76, 97 and 0 %, respectively. The results were compared with the simulation, considering the real scenario (C0) with 54 % forest cover. The Nash-Sutcliffe coefficients were 0.65 and 0.70 for calibration and validation, respectively, indicating satisfactory results in the flow simulation. A mean reduction of 10 % of the native forest cover would cause a mean annual increase of approximately 11.5 mm in total runoff at the watershed outlet. Reforestation would ensure minimum flows in the dry period and regulate the maximum flow of the main watercourse of the watershed.

Water infiltration in an ultisol after cultivation of common bean

Water infiltration in the soil is an important hydrological process that occurs at the interface of the soil-atmosphere system; thus, the soil management practice used has a strong influence on this process. The aim of this study was to evaluate water infiltration in the soil and compare equations for estimating the water infiltration rate in an Ultisol after harvesting common bean (Phaseolus vulgaris L.) under simulated rainfall. Field tests with a rainfall simulator were carried out in three soil management systems: minimum tillage (MT), conventional tillage (CT), and no tillage (NT). In NT, four levels of plant residue on the soil surface were evaluated: 0, 3, 6, and 9 t ha-1. The models of Kostiakov-Lewis, Horton, and Philip were used to estimate the infiltration rate. In the MT system, the final infiltration rate was 54 mm h-1, whereas in the CT and NT systems with up to 3 t ha-1 of plant residue on the soil surface, the rate was near 17 mm h-1. In addition, the results indicated that in the NT system the infiltration rate increased with plant residue coverage greater than 6 t ha-1, i.e., there was a positive correlation between plant cover and the water infiltration rate. The Horton model was the most suitable in representing the water infiltration process in the soil. Therefore, this model can be recommended for estimation of this variable regardless of the soil tillage system used.

Quantification of the least limiting water range in an oxisol using two methodological strategies

The least limiting water range (LLWR) has been used as an indicator of soil physical quality as it represents, in a single parameter, the soil physical properties directly linked to plant growth, with the exception of temperature. The usual procedure for obtaining the LLWR involves determination of the water retention curve (WRC) and the soil resistance to penetration curve (SRC) in soil samples with undisturbed structure in the laboratory. Determination of the WRC and SRC using field measurements (in situ ) is preferable, but requires appropriate instrumentation. The objective of this study was to determine the LLWR from the data collected for determination of WRC and SRC in situ using portable electronic instruments, and to compare those determinations with the ones made in the laboratory. Samples were taken from the 0.0-0.1 m layer of a Latossolo Vermelho distrófico (Oxisol). Two methods were used for quantification of the LLWR: the traditional, with measurements made in soil samples with undisturbed structure; and in situ , with measurements of water content (θ), soil water potential (Ψ), and soil resistance to penetration (SR) through the use of sensors. The in situ measurements of θ, Ψ and SR were taken over a period of four days of soil drying. At the same time, samples with undisturbed structure were collected for determination of bulk density (BD). Due to the limitations of measurement of Ψ by tensiometer, additional determinations of θ were made with a psychrometer (in the laboratory) at the Ψ of -1500 kPa. The results show that it is possible to determine the LLWR by the θ, Ψ and SR measurements using the suggested approach and instrumentation. The quality of fit of the SRC was similar in both strategies. In contrast, the θ and Ψ in situ measurements, associated with those measured with a psychrometer, produced a better WRC description. The estimates of the LLWR were similar in both methodological strategies. The quantification of LLWR in situ can be achieved in 10 % of the time required for the traditional method.

Soil water erosion under different cultivation systems and different fertilization rates and forms over 10 years

The action of rain and surface runoff together are the active agents of water erosion, and further influences are the soil type, terrain, soil cover, soil management, and conservation practices. Soil water erosion is low in the no-tillage management system, being influenced by the amount and form of lime and fertilizer application to the soil, among other factors. The aim was to evaluate the effect of the form of liming, the quantity and management of fertilizer application on the soil and water losses by erosion under natural rainfall. The study was carried out between 2003 and 2013 on a Humic Dystrupept soil, with the following treatments: T1 - cultivation with liming and corrective fertilizer incorporated into the soil in the first year, and with 100 % annual maintenance fertilization of P and K; T2 - surface liming and corrective fertilization distributed over five years, and with 75 % annual maintenance fertilization of P and K; T3 - surface liming and corrective fertilization distributed over three years, and with 50 % annual maintenance fertilization of P and K; T4 - surface liming and corrective fertilization distributed over two years, and with 25 % annual maintenance fertilization of P and K; T5 - fallow soil, without liming or fertilization. In the rotation the crops black oat (Avena strigosa ), soybean (Glycine max ), common vetch (Vicia sativa ), maize (Zea mays ), fodder radish (Raphanus sativus ), and black beans (Phaseolus vulgaris ). The split application of lime and mineral fertilizer to the soil surface in a no-tillage system over three and five years, results in better control of soil losses than when split in two years. The increase in the amount of fertilizer applied to the soil surface under no-tillage cultivation increases phytomass production and reduces soil loss by water erosion. Water losses in treatments under no-tillage cultivation were low in all crop cycles, with a similar behavior as soil losses.

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.

WATER INFILTRATION IN TWO CULTIVATED SOILS IN SOUTHERN BRAZIL

Infiltration is the passage of water through the soil surface, influenced by the soil type and cultivation and by the soil roughness, surface cover and water content. Infiltration absorbs most of the rainwater and is therefore crucial for planning mechanical conservation practices to manage runoff. This study determined water infiltration in two soil types under different types of management and cultivation, with simulated rainfall of varying intensity and duration applied at different times, and to adjust the empirical model of Horton to the infiltration data. The study was conducted in southern Brazil, on Dystric Nitisol (Nitossolo Bruno aluminoférrico húmico) and Humic Cambisol (Cambissolo Húmico alumínico léptico) soils to assess the following situations: simulated rains on the Nitisol from 2001 to 2012 in 31 treatments, differing in crop type, sowing direction, type of soil opener on the seeder, amount and type of crop residue and amount of liquid swine manure applied; on the Cambisol, rains were simlated from 2006 to 2012 and 18 treatments were evaluated, differing in crop, seeding direction and crop residue type. The constant of the water infiltration rate into the soil varies significantly with the soil type (30.2 mm h-1 in the Nitisol and 6.6 mm h-1 in the Cambisol), regardless of the management system, application time and rain intensity and duration. At the end of rainfalls, soil-water infiltration varies significantly with the management system, with the timing of application and rain intensity and duration, with values ranging from 13 to 59 mm h-1, in the two studied soils. The characteristics of the sowing operation in terms of relief, crop type and amount and type of crop residue influenced soil water infiltration: in the Nitisol, the values of contour and downhill seeding vary between 27 and 43 mm h-1, respectively, with crop residues of corn, wheat and soybean while in the Cambisol, the variation is between 2 and 36 mm h-1, respectively, in soybean and corn crops. The Horton model fits the values of water infiltration rate into the soil, resulting in the equation i = 30.2 + (68.2 - 30.2) e-0.0371t (R2 = 0.94**) for the Nitisol and i = 6.6 + (64.5 - 6.6) e-0.0537t (R2 = 0.99**) for the Cambisol.

BEYOND THE “LEAST LIMITING WATER RANGE”: RETHINKING SOIL PHYSICS RESEARCH IN BRAZIL

As opposed to objective definitions in soil physics, the subjective term “soil physical quality” is increasingly found in publications in the soil physics area. A supposed indicator of soil physical quality that has been the focus of attention, especially in the Brazilian literature, is the Least Limiting Water Range (RLL), translated in Portuguese as "Intervalo Hídrico Ótimo" or IHO. In this paper the four limiting water contents that define RLLare discussed in the light of objectively determinable soil physical properties, pointing to inconsistencies in the RLLdefinition and calculation. It also discusses the interpretation of RLL as an indicator of crop productivity or soil physical quality, showing its inability to consider common phenological and pedological boundary conditions. It is shown that so-called “critical densities” found by the RLL through a commonly applied calculation method are questionable. Considering the availability of robust models for agronomy, ecology, hydrology, meteorology and other related areas, the attractiveness of RLL as an indicator to Brazilian soil physicists is not related to its (never proven) effectiveness, but rather to the simplicity with which it is dealt. Determining the respective limiting contents in a simplified manner, relegating the study or concern on the actual functioning of the system to a lower priority, goes against scientific construction and systemic understanding. This study suggests a realignment of the research in soil physics in Brazil with scientific precepts, towards mechanistic soil physics, to replace the currently predominant search for empirical correlations below the state of the art of soil physics.

Water Balance Components in Covered and Uncovered Soil Growing Irrigated Muskmelon

ABSTRACT Knowledge of the terms (or processes) of the soil water balance equation or simply the components of the soil water balance over the cycle of an agricultural crop is essential for soil and water management. Thus, the aim of this study was to analyze these components in a Cambissolo Háplico (Haplocambids) growing muskmelon (Cucumis melo L.) under drip irrigation, with covered and uncovered soil, in the municipality of Baraúna, State of Rio Grande do Norte, Brazil (05º 04’ 48” S, 37º 37’ 00” W). Muskmelon, variety AF-646, was cultivated in a flat experimental area (20 × 50 m). The crop was spaced at 2.00 m between rows and 0.35 m between plants, in a total of ten 50-m-long plant rows. At points corresponding to ⅓ and ⅔ of each plant row, four tensiometers (at a distance of 0.1 m from each other) were set up at the depths of 0.1, 0.2, 0.3, and 0.4 m, adjacent to the irrigation line (0.1 m from the plant row), between two selected plants. Five random plant rows were mulched using dry leaves of banana (Musa sp.) along the drip line, forming a 0.5-m-wide strip, which covered an area of 25 m2 per of plant row with covered soil. In the other five rows, there was no covering. Thus, the experiment consisted of two treatments, with 10 replicates, in four phenological stages: initial (7-22 DAS - days after sowing), growing (22-40 DAS), fruiting (40-58 DAS) and maturation (58-70 DAS). Rainfall was measured with a rain gauge and water storage was estimated by the trapezoidal method, based on tensiometer readings and soil water retention curves. For soil water flux densities at 0.3 m, the tensiometers at the depths of 0.2, 0.3, and 0.4 m were considered; the tensiometer at 0.3 m was used to estimate soil water content from the soil water retention curve at this depth, and the other two to calculate the total potential gradient. Flux densities were calculated through use of the Darcy-Buckingham equation, with hydraulic conductivity determined by the instantaneous profile method. Crop actual evapotranspiration was calculated as the unknown of the soil water balance equation. The soil water balance method is effective in estimating the actual evapotranspiration of irrigated muskmelon; there was no significant effect of soil coverage on capillary rise, internal drainage, crop actual evapotranspiration, and muskmelon yield compared with the uncovered soil; the transport of water caused by evaporation in the uncovered soil was controlled by the break in capillarity at the soil-atmosphere interface, which caused similar water dynamics for both management practices applied.

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.