Physical quality of an Oxisol after different periods of management systems

Management systems may lead to a loss of soil physical quality as a result of removal of the plant cover and excessive agricultural mechanization. The hypothesis of this study was that the soil aggregate stability, bulk density, macro- and microporosity, and the S index and saturated hydraulic conductivity may be used as indicators of the soil physical quality. The aim was to study the effects of different periods and managements on the physical attributes of a medium-textured Red Oxisol under soybean and corn for two growing seasons, and determine which layers are most susceptible to variations. A completely randomized experimental design was used with split plots (five treatments and four layers), with four replications. The treatments in 2008/09 consisted of: five years of no-tillage (NTS5), seven years of no-tillage (NTS7), nine years of no-tillage (NTS9), conventional tillage (CTS) and an adjacent area of native forest (NF). The treatments were extended for another year, identified in 2009/10 as: NTS6, NTS8, NTS10, CTS and NF. The soil layers 0-0.05, 0.05-0.10, 0.10-0.20 and 0.20-0.30 m were sampled. The highest S index values were observed in the treatment CTS in the 0-0.05 m layer (0.106) and the 0.05-0.10 m layer (0.099) in 2008/09, and in the 0-0.05 m layer (0.066) in 2009/10. This fact may be associated with soil turnover, resulting in high macroporosity in this treatment. In contrast, in the NTS, limiting macroporosity values were observed in some layers (below 0.10 m³ m-3). Highest aggregate stability as well as the highest saturated hydraulic conductivity (Kθ) values were observed in NF in relation to the other treatments. In 2009/10, the Kθ in NF differed only from NTS10. This study showed that the use of the S index alone cannot be recommended as an absolute indicator of the soil physical quality, even at values greater than 0.035.

Estimation of water percolation by different methods using TDR

Detailed knowledge on water percolation into the soil in irrigated areas is fundamental for solving problems of drainage, pollution and the recharge of underground aquifers. The aim of this study was to evaluate the percolation estimated by time-domain-reflectometry (TDR) in a drainage lysimeter. We used Darcy's law with K(θ) functions determined by field and laboratory methods and by the change in water storage in the soil profile at 16 points of moisture measurement at different time intervals. A sandy clay soil was saturated and covered with plastic sheet to prevent evaporation and an internal drainage trial in a drainage lysimeter was installed. The relationship between the observed and estimated percolation values was evaluated by linear regression analysis. The results suggest that percolation in the field or laboratory can be estimated based on continuous monitoring with TDR, and at short time intervals, of the variations in soil water storage. The precision and accuracy of this approach are similar to those of the lysimeter and it has advantages over the other evaluated methods, of which the most relevant are the possibility of estimating percolation in short time intervals and exemption from the predetermination of soil hydraulic properties such as water retention and hydraulic conductivity. The estimates obtained by the Darcy-Buckingham equation for percolation levels using function K(θ) predicted by the method of Hillel et al. (1972) provided compatible water percolation estimates with those obtained in the lysimeter at time intervals greater than 1 h. The methods of Libardi et al. (1980), Sisson et al. (1980) and van Genuchten (1980) underestimated water percolation.

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.

Effects of pig slurry application on soil physical and chemical properties and glyphosate mobility

Pig slurry applied to soil at different rates may affect soil properties and the mobility of chemical compounds within the soil. The purpose of this study was to evaluate the effects of rates of pig slurry application in agricultural areas on soil physical and chemical properties and on the mobility of glyphosate through the soil profile. The study was carried out in the 12th year of an experiment with pig slurry applied at rates of 0 (control), 50, 100 and 200 m³ ha-1 yr-1 on a Latossolo Vermelho distrófico (Hapludox) soil. In the control, the quantities of P and K removed by harvested grains were replaced in the next crop cycle. Soil physical properties (bulk density, porosity, texture, and saturated hydraulic conductivity) and chemical properties (organic matter, pH, extractable P, and exchangeable K) were measured. Soil solution samples were collected at depths of 20, 40 and 80 cm using suction lysimeters, and glyphosate concentrations were measured over a 60-day period after slurry application. Soil physical and chemical properties were little affected by the pig slurry applications, but soil pH was reduced and P levels increased in the surface layers. In turn, K levels were increased in sub-surface layers. Glyphosate concentrations tended to decrease over time but were not affected by pig slurry application. The concentrations of glyphosate found in different depths show that the pratice of this application in agricultural soils has the potential for contamination of groundwater, especially when the water table is the surface and heavy rains occur immediately after application.

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.

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.

Soil Physical Quality and Soybean Yield as Affected by Chiseling and Subsoiling of a No-Till Soil

ABSTRACT The concept of soil physical quality (SPQ) is currently under discussion, and an agreement about which soil physical properties should be included in the SPQ characterization has not been reached. The objectives of this study were to evaluate the ability of SPQ indicators based on static and dynamic soil properties to assess the effects of two loosening treatments (chisel plowing to 0.20 m [ChT] and subsoiling to 0.35 m [DL]) on a soil under NT and to compare the performance of static- and dynamic-based SPQ indicators to define soil proper soil conditions for soybean yield. Soil sampling and field determinations were carried out after crop harvest. Soil water retention curve was determined using a tension table, and field infiltration was measured using a tension disc infiltrometer. Most dynamic SPQ indicators (field saturated hydraulic conductivity, K0, effective macroporosity, εma, total connectivity and macroporosity indexes [CwTP and Cwmac]) were affected by the studied treatments, and were greater for DL compared to NT and ChT (K0 values were 2.17, 2.55, and 4.37 cm h-1 for NT, ChT, and DL, respectively). However, static SPQ indicators (calculated from the water retention curve) were not capable of distinguishing effects among treatments. Crop yield was significantly lower for the DL treatment (NT: 2,400 kg ha-1; ChT: 2,358 kg ha-1; and DL: 2,105 kg ha1), in agreement with significantly higher values of the dynamic SPQ indicators, K0, εma, CwTP, and Cwmac, in this treatment. The results support the idea that SPQ indicators based on static properties are not capable of distinguishing tillage effects and predicting crop yield, whereas dynamic SPQ indicators are useful for distinguishing tillage effects and can explain differences in crop yield when used together with information on weather conditions. However, future studies, monitoring years with different weather conditions, would be useful for increasing knowledge on this topic.

Integration of local-scale hydrological and regional-scale geophysical based on a nonlinear Bayesian sequential simulation approach

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the regional scale represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed a downscaling procedure based on a non-linear Bayesian sequential simulation approach. The basic objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity, which is available throughout the model space. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariate kernel density function. This method is then applied to the stochastic integration of low-resolution, re- gional-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities. Finally, the overall viability of this downscaling approach is tested and verified by performing and comparing flow and transport simulation through the original and the downscaled hydraulic conductivity fields. Our results indicate that the proposed procedure does indeed allow for obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.

Effects of long-term irrigation with treated wastewater. Part I: Evolution of soil physico-chemical properties

The long-term impact of irrigation on a Mediterranean sandy soil irrigated with Treated wastewater (TWW) since 1980 was evaluated. The main soil properties (CEC, pH, size distribution, exchangeable cations and chloride, hydraulic conductivity) as well as the organic matter and Cu, Cr and Pb speciation in an irrigated soil and a non-irrigated control soil at various soil depths were monitored and compared during a 2 years experiment. In this first part, the evolution of the physico-chemical soil properties was described. The irrigation with TWW was beneficial with regard to water and nutrient supplying. All the exchangeable cations other than K(+) were higher in the irrigated soil than in the reference one. A part of the exchangeable cations was not fixed on the exchange complex but stored as labile salts or in concentrated soil solution. Despite the very sandy soil texture, both saturated and unsaturated hydraulic conductivity exhibited a significant diminution in the irrigated soil, but remained high enough to allow water percolation during rainy periods and subsequent leaching of accumulated salts, preventing soil salinization. In the irrigated soil, exchangeable sodium percentage (ESP) exhibited high values (20% on average) and the soil organic C was lower than in the reference. No significant effect was noticed on soil mineralogical composition due to irrigation. (C) 2010 Published by Elsevier Ltd.

Optimizing Pavement Base, Subbase, and Subgrade Layers for Cost and Performance of Local Roads, TR-640

This report is one of two products for this project with the other being a design guide. This report describes test results and comparative analysis from 16 different portland cement concrete (PCC) pavement sites on local city and county roads in Iowa. At each site the surface conditions of the pavement (i.e., crack survey) and foundation layer strength, stiffness, and hydraulic conductivity properties were documented. The field test results were used to calculate in situ parameters used in pavement design per SUDAS and AASHTO (1993) design methodologies. Overall, the results of this study demonstrate how in situ and lab testing can be used to assess the support conditions and design values for pavement foundation layers and how the measurements compare to the assumed design values. The measurements show that in Iowa, a wide range of pavement conditions and foundation layer support values exist. The calculated design input values for the test sites (modulus of subgrade reaction, coefficient of drainage, and loss of support) were found to be different than typically assumed. This finding was true for the full range of materials tested. The findings of this study support the recommendation to incorporate field testing as part of the process to field verify pavement design values and to consider the foundation as a design element in the pavement system. Recommendations are provided in the form of a simple matrix for alternative foundation treatment options if the existing foundation materials do not meet the design intent. The PCI prediction model developed from multi-variate analysis in this study demonstrated a link between pavement foundation conditions and PCI. The model analysis shows that by measuring properties of the pavement foundation, the engineer will be able to predict long term performance with higher reliability than by considering age alone. This prediction can be used as motivation to then control the engineering properties of the pavement foundation for new or re-constructed PCC pavements to achieve some desired level of performance (i.e., PCI) with time.

Guidance for Improving Foundation Layers to Increase Pavement Performance on Local Roads, TR-640, 2014

This document is the second of two deliverables for the project Optimizing Pavement Base, Subbase, and Subgrade Layers for Cost and Performance on Local Roads (TR-640). The first deliverable is the 454-page Final Field Data Report. The field data report describes test results and comparative analysis from 16 different portland cement concrete (PCC) pavement sites on local city and county roads in Iowa. At each site the surface conditions of the pavement (i.e., crack survey) and foundation layer strength, stiffness, and hydraulic conductivity properties were documented. The field test results were used to calculate in situ parameters used in pavement design methodologies for AASHTO (1993) and Iowa’s Statewide Urban Design and Specifications (SUDAS). Overall, the results of the study demonstrate how in situ and lab testing can be used to assess the support conditions and design values for pavement foundation layers and how the measurements compare to the assumed design values. This guide summarizes the study results and outlines general guidelines for applying them to optimize pavement bases, subbases, and subgrade layers of local roads with PCC pavements and thus their performance.

Optimization and Management of Materials in Earthwork Construction, TR-501, 2010

As a result of forensic investigations of problems across Iowa, a research study was developed aimed at providing solutions to identified problems through better management and optimization of the available pavement geotechnical materials and through ground improvement, soil reinforcement, and other soil treatment techniques. The overall goal was worked out through simple laboratory experiments, such as particle size analysis, plasticity tests, compaction tests, permeability tests, and strength tests. A review of the problems suggested three areas of study: pavement cracking due to improper management of pavement geotechnical materials, permeability of mixed-subgrade soils, and settlement of soil above the pipe due to improper compaction of the backfill. This resulted in the following three areas of study: (1) The optimization and management of earthwork materials through general soil mixing of various select and unsuitable soils and a specific example of optimization of materials in earthwork construction by soil mixing; (2) An investigation of the saturated permeability of compacted glacial till in relation to validation and prediction with the Enhanced Integrated Climatic Model (EICM); and (3) A field investigation and numerical modeling of culvert settlement. For each area of study, a literature review was conducted, research data were collected and analyzed, and important findings and conclusions were drawn. It was found that optimum mixtures of select and unsuitable soils can be defined that allow the use of unsuitable materials in embankment and subgrade locations. An improved model of saturated hydraulic conductivity was proposed for use with glacial soils from Iowa. The use of proper trench backfill compaction or the use of flowable mortar will reduce the potential for developing a bump above culverts.

Spectral induced polarization: Laboratory measurements on sandy quartz samples of varied pore characteristics

Understanding the influence of pore space characteristics on the hydraulic conductivity and spectral induced polarization (SIP) response is critical for establishing relationships between the electrical and hydrological properties of surficial sedimentary deposits. Here, we present the results of laboratory SIP measurements on saturated quartz samples with granulometric characteristics ranging from fine sand to fine gravel. We alter the pore characteristics using three principal methods: (i) variation of the grain sizes, (ii) changing the degree of compaction, and (iii) changing the level of sorting. We then examine how these changes affect both the SIP response and the hydraulic conductivity. In general, the results indicate a clear connection between the applied changes in pore characteristics and the SIP response. In particular, we observe a systematic correlation between the hydraulic conductivity and the relaxation time of the Cole-Cole model describing the observed SIP effect for the whole range of considered grain sizes.

Development of a Subgrade Drainage Model for Unpaved Roads, TR-654, 2015

With over 68 thousand miles of gravel roads in Iowa and the importance of these roads within the farm-to-market transportation system, proper water management becomes critical for maintaining the integrity of the roadway materials. However, the build-up of water within the aggregate subbase can lead to frost boils and ultimately potholes forming at the road surface. The aggregate subbase and subgrade soils under these gravel roads are produced with material opportunistically chosen from local sources near the site and, many times, the compositions of these sublayers are far from ideal in terms of proper water drainage with the full effects of this shortcut not being fully understood. The primary objective of this project was to provide a physically-based model for evaluating the drainability of potential subbase and subgrade materials for gravel roads in Iowa. The Richards equation provided the appropriate framework to study the transient unsaturated flow that usually occurs through the subbase and subgrade of a gravel road. From which, we identified that the saturated hydraulic conductivity, Ks, was a key parameter driving the time to drain of subgrade soils found in Iowa, thus being a good proxy variable for accessing roadway drainability. Using Ks, derived from soil texture, we were able to identify potential problem areas in terms of roadway drainage . It was found that there is a threshold for Ks of 15 cm/day that determines if the roadway will drain efficiently, based on the requirement that the time to drain, Td, the surface roadway layer does not exceed a 2-hr limit. Two of the three highest abundant textures (loam and silty clay loam), which cover nearly 60% of the state of Iowa, were found to have average Td values greater than the 2-hr limit. With such a large percentage of the state at risk for the formation of boils due to the soil with relatively low saturated hydraulic conductivity values, it seems pertinent that we propose alternative design and/or maintenance practices to limit the expensive repair work in Iowa. The addition of drain tiles or French mattresses my help address drainage problems. However, before pursuing this recommendation, a comprehensive cost-benefit analysis is needed.

Bayesian Markov-chain-Monte-Carlo inversion of time-lapse cross hole ground-penetrating radar data to characterize the vadose zone at the Arrenaes field site, Denmark

The ground-penetrating radar (GPR) geophysical method has the potential to provide valuable information on the hydraulic properties of the vadose zone because of its strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR traveltime data can allow for a significant reduction in uncertainty regarding subsurface van Genuchten-Mualem (VGM) parameters. Much of the previous work on the stochastic estimation of VGM parameters from crosshole GPR data has considered the case of steady-state infiltration conditions, which represent only a small fraction of practically relevant scenarios. We explored in detail the dynamic infiltration case, specifically examining to what extent time-lapse crosshole GPR traveltimes, measured during a forced infiltration experiment at the Arreneas field site in Denmark, could help to quantify VGM parameters and their uncertainties in a layered medium, as well as the corresponding soil hydraulic properties. We used a Bayesian Markov-chain-Monte-Carlo inversion approach. We first explored the advantages and limitations of this approach with regard to a realistic synthetic example before applying it to field measurements. In our analysis, we also considered different degrees of prior information. Our findings indicate that the stochastic inversion of the time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions compared with the corresponding priors, which in turn significantly improves knowledge of soil hydraulic properties. Overall, the results obtained clearly demonstrate the value of the information contained in time-lapse GPR data for characterizing vadose zone dynamics.