Numerical modeling of the effect of variation of boundary conditions on vadose zone hydraulic properties

An accurate estimation of hydraulic fluxes in the vadose zone is essential for the prediction of water, nutrient and contaminant transport in natural systems. The objective of this study was to simulate the effect of variation of boundary conditions on the estimation of hydraulic properties (i.e. water content, effective unsaturated hydraulic conductivity and hydraulic flux) in a one-dimensional unsaturated flow model domain. Unsaturated one-dimensional vertical water flow was simulated in a pure phase clay loam profile and in clay loam interlayered with silt loam distributed according to the third iteration of the Cantor Bar fractal object Simulations were performed using the numerical model Hydrus 1D. The upper and lower pressure heads were varied around average values of -55 cm for the near-saturation range. This resulted in combinations for the upper and lower constant head boundary conditions, respectively, of -50 and -60 cm, -40 and -70 cm, -30 and -80 cm, -20 and -90 cm, and -10 and -100 cm. For the drier range the average head between the upper and lower boundary conditions was set to -550 cm, resulting in the combinations -500 and -600 cm, -400 and -700 cm, -300 and -800 cm, -200 and -900 cm, and -100 and -1,000 cm, for upper and lower boundary conditions, respectively. There was an increase in water contents, fluxes and hydraulic conductivities with the increase in head difference between boundary conditions. Variation in boundary conditions in the pure phase and interlayered one-dimensional profiles caused significant deviations in fluxes, water contents and hydraulic conductivities compared to the simplest case (a head difference between the upper and lower constant head boundaries of 10 cm in the wetter range and 100 cm in the drier range).

Method to estimate soil macroporosity and microporosity based on sand content and bulk density

Macroporosity is often used in the determination of soil compaction. Reduced macroporosity can lead to poor drainage, low root aeration and soil degradation. The aim of this study was to develop and test different models to estimate macro and microporosity efficiently, using multiple regression. Ten soils were selected within a large range of textures: sand (Sa) 0.07-0.84; silt 0.03-0.24; clay 0.13-0.78 kg kg-1 and subjected to three compaction levels (three bulk densities, BD). Two models with similar accuracy were selected, with a mean error of about 0.02 m³ m-3 (2 %). The model y = a + b.BD + c.Sa, named model 2, was selected for its simplicity to estimate Macro (Ma), Micro (Mi) or total porosity (TP): Ma = 0.693 - 0.465 BD + 0.212 Sa; Mi = 0.337 + 0.120 BD - 0.294 Sa; TP = 1.030 - 0.345 BD 0.082 Sa; porosity values were expressed in m³ m-3; BD in kg dm-3; and Sa in kg kg-1. The model was tested with 76 datum set of several other authors. An error of about 0.04 m³ m-3 (4 %) was observed. Simulations of variations in BD as a function of Sa are presented for Ma = 0 and Ma = 0.10 (10 %). The macroporosity equation was remodeled to obtain other compaction indexes: a) to simulate maximum bulk density (MBD) as a function of Sa (Equation 11), in agreement with literature data; b) to simulate relative bulk density (RBD) as a function of BD and Sa (Equation 13); c) another model to simulate RBD as a function of Ma and Sa (Equation 16), confirming the independence of this variable in relation to Sa for a fixed value of macroporosity and, also, proving the hypothesis of Hakansson & Lipiec that RBD = 0.87 corresponds approximately to 10 % macroporosity (Ma = 0.10 m³ m-3).

Chronofunctions of Heilu soil developed from Loess in Luochuan, on the chinese Loess plateau

Soil chronofunctions are an alternative for the quantification of soil-forming processes and underlie the modeling of soil genesis. To establish soil chronofunctions of a Heilu soil profile on Loess in Luochuan, selected soil properties and the 14C ages in the Holocene were studied. Linear, logarithmic, and third-order polynomial functions were selected to fit the relationships between soil properties and ages. The results indicated that third-order polynomial function fit best for the relationships between clay (< 0.002 mm), silt (0.002-0.02 mm), sand (0.02-2 mm) and soil ages, and a trend of an Ah horizon ocurrence in the profile. The logarithmic function indicated mainly variations of soil organic carbon and pH with time (soil age). The variation in CaCO3 content, Mn/Zr, Fe/Zr, K/Zr, Mg/Zr, Ca/Zr, P/Zr, and Na/Zr ratios with soil age were best described by three-order polynomial functions, in which the trend line showed migration of CaCO3 and some elements.

Spatial patterns of preconsolidation pressure and soil moisture along transects in two directions under coffee

Information on the spatial structure of soil physical and structural properties is needed to evaluate the soil quality. The purpose of this study was to investigate the spatial behavior of preconsolidation pressure and soil moisture in six transects, three selected along and three across coffee rows, at three different sites under different tillage management systems. The study was carried out on a farm, in Patrocinio, state of Minas Gerais, in the Southeast of Brazil (18 º 59 ' 15 '' S; 46 º 56 ' 47 '' W; 934 m asl). The soil type is a typic dystrophic Red Latosol (Acrustox) and consists of 780 g kg-1 clay; 110 g kg-1 silt and 110 g kg-1 sand, with an average slope of 3 %. Undisturbed soil cores were sampled at a depth of 0.10-0.13 m, at three different points within the coffee plantation: (a) from under the wheel track, where equipment used in farm operations passes; (b) in - between tracks and (c) under the coffee canopy. Six linear transects were established in the experimental area: three transects along and three across the coffee rows. This way, 161 samples were collected in the transect across the coffee rows, from the three locations, while 117 samples were collected in the direction along the row. The shortest sampling distance in the transect across the row was 4 m, and 0.5 m for the transect along the row. No clear patterns of the preconsolidation pressure values were observed in the 200 m transect. The results of the semivariograms for both variables indicated a high nugget value and short range for the studied parameters of all transects. A cyclic pattern of the parameters was observed for the across-rows transect. An inverse relationship between preconsolidation pressure and soil moisture was clearly observed in the samples from under the track, in both directions.

Effects of erosion and deposition on particle size distribution of deposited farmland soils on the chinese loess plateau

Particle size distribution (PSD) in the soil profile is strongly related to erosion, deposition, and physical and chemical processes. Water cycling and plant growth are also affected by PSD. Material sedimented upstream of the dam constructions formed large areas of deposited farmland (DF) soils on the Chinese Loess Plateau (CLP), which has been the site of the most severe soil erosion in the world. Two DFs without tillage on the CLP were chosen to study the combined effect of erosion and check dams on PSD. Eighty-eight layers (each 10 cm thick) of filled deposited farmland (FDF) soils and 22 layers of silting deposited farmland (SDF) soils of each studied soil profile were collected and 932 soil samples were investigated using laser granulometry. The particle sizes were stratified in both DFs based on soil properties and erosion resistance. The obtained results of clay and silt fractions showed similar horizontal distribution, indicating parallel characteristics of erosion and deposition processes. Fine sand represented the largest fraction, suggesting the preferential detachment of this fraction. The most erodible range of particle sizes was 0.25-0.5 mm, followed by 0.2-0.25 mm in the studied soil profiles. The correlation between particle size and soil water contents tended to increase with increasing water contents in FDF. Due to the abundant shallow groundwater, the relationship between particle size and soil water content in SDF was lost. Further studies on PSD in the DF area are needed to enhance the conservation management of soil and water resources in this region.

Modelling of soil penetration resistance for an oxisol under no-tillage

Soil penetration resistance is an important property that affects root growth and elongation and water movement in the soil. Since no-till systems tend to increase organic matter in the soil, the purpose of this study was to evaluate the efficiency with which soil penetration resistance is estimated using a proposed model based on moisture content, density and organic matter content in an Oxisol containing 665, 221 and 114 g kg-1 of clay, silt and sand respectively under annual no-till cropping, located in Londrina, Paraná State, Brazil. Penetration resistance was evaluated at random locations continually from May 2008 to February 2011, using an impact penetrometer to obtain a total of 960 replications. For the measurements, soil was sampled at depths of 0 to 20 cm to determine gravimetric moisture (G), bulk density (D) and organic matter content (M). The penetration resistance curve (PR) was adjusted using two non-linear models (PR = a Db Gc and PR' = a Db Gc Md), where a, b, c and d are coefficients of the adjusted model. It was found that the model that included M was the most efficient for estimating PR, explaining 91 % of PR variability, compared to 82 % of the other model.

Runoff and sediment transport in a degraded area

Gully erosion occurs by the combined action of splash, sheetwash and rill-wash (interrill and rill erosion). These erosion processes have a great capacity for both sediment production and sediment transport. The objectives of this experiment were to evaluate hydrological and sediment transport in a degraded area, severely dissected by gullies; to assess the hydraulic flow characteristics and their aggregate transport capacity; and to measure the initial splash erosion rate. In the study area in Guarapuava, State of Paraná, Brazil (lat 25º 24' S; long 51º24' W; 1034 m asl), the soil was classified as Cambissolo Húmico alumínico, with the following particle-size composition: sand 0.116 kg kg-1; silt 0.180 kg kg-1; and clay 0.704 kg kg-1. The approach of this research was based on microcatchments formed in the ground, to study the hydrological response and sediment transport. A total of eight rill systems were simulated with dry and wet soil. An average rainfall of 33.7 ± 4.0 mm was produced for 35 to 54 min by a rainfall simulator. The equipment was installed, and a trough was placed at the end of the rill to collect sediments and water. During the simulation, the following variables were measured: time to runoff, time to ponding, time of recession, flow velocity, depth, ratio of the initial splash and grain size. The rainsplash of dry topsoil was more than twice as high as under moist conditions (5 g m-2 min-1 and 2 g m-2 min-1, respectively). The characteristics of the flow hydraulics indicate transition from laminar to turbulent flow [Re (Reynolds number) 1000-2000]. In addition, it was observed that a flow velocity of 0.12 m s-1 was the threshold for turbulent flow (Re > 2000), especially at the end of the rainfall simulation. The rill flow tended to be subcritical [Fr (Froude Number) < 1.0]. The variation in hydrological attributes (infiltration and runoff) was lower, while the sediment yield was variable. The erosion in the rill systems was characterized as limited transport, although the degraded area generated an average of 394 g m-2 of sediment in each simulation.

Comparison of the methane oxidation rate in four media

Landfill gas emissions are one of the main sources of anthropogenic methane (CH4), a major greenhouse gas. In this paper, an economically attractive alternative to minimize greenhouse gas emissions from municipal solid waste landfills was sought. This alternative consists in special biofilters as landfill covers with oxidative capacity in the presence of CH4. To improve the quality/cost ratio of the project, compost was chosen as one of the cover substrates and soil (Typic red yellow-silt-clay Podzolic) as the other. The performance of four substrates was studied in laboratory experiments: municipal solid waste (MSW) compost, soil, and two soil-compost at different proportions. This study aimed to evaluate the suitability and environmental compatibility as a means of CH4 oxidation in biofilters. Four biofilters were constructed in 60 cm PVC tubes with an internal diameter of 10 cm. Each filter contained 2.3 L of oxidizing substrate at the beginning of the experiment. The gas used was a mixture of CH4 and air introduced at the bottom of each biofilter, at a flow of 150 mL min-1, by a flow meter. One hundred days after the beginning of the experiment, the best biofilter was the MSW compost with an oxidation rate of 990 g m-3 day-1 , corresponding to an efficiency of 44 %. It can be concluded that the four substrates studied have satisfactory oxidative capacity, and the substrates can be used advantageously as cover substrate of MSW landfills.

Performance of a reciprocal shaker in mechanical dispersion of soil samples for particle-size analysis

The dispersion of the samples in soil particle-size analysis is a fundamental step, which is commonly achieved with a combination of chemical agents and mechanical agitation. The purpose of this study was to evaluate the efficiency of a low-speed reciprocal shaker for the mechanical dispersion of soil samples of different textural classes. The particle size of 61 soil samples was analyzed in four replications, using the pipette method to determine the clay fraction and sieving to determine coarse, fine and total sand fractions. The silt content was obtained by difference. To evaluate the performance, the results of the reciprocal shaker (RSh) were compared with data of the same soil samples available in reports of the Proficiency testing for Soil Analysis Laboratories of the Agronomic Institute of Campinas (Prolab/IAC). The accuracy was analyzed based on the maximum and minimum values defining the confidence intervals for the particle-size fractions of each soil sample. Graphical indicators were also used for data comparison, based on dispersion and linear adjustment. The descriptive statistics indicated predominantly low variability in more than 90 % of the results for sand, medium-textured and clay samples, and for 68 % of the results for heavy clay samples, indicating satisfactory repeatability of measurements with the RSh. Medium variability was frequently associated with silt, followed by the fine sand fraction. The sensitivity analyses indicated an accuracy of 100 % for the three main separates (total sand, silt and clay), in all 52 samples of the textural classes heavy clay, clay and medium. For the nine sand soil samples, the average accuracy was 85.2 %; highest deviations were observed for the silt fraction. In relation to the linear adjustments, the correlation coefficients of 0.93 (silt) or > 0.93 (total sand and clay), as well as the differences between the angular coefficients and the unit < 0.16, indicated a high correlation between the reference data (Prolab/IAC) and results obtained with the RSh. In conclusion, the mechanical dispersion by the reciprocal shaker of soil samples of different textural classes was satisfactory. The results allowed recommending the use of the equipment at low agitation for particle size- analysis. The advantages of this Brazilian apparatus are its low cost, the possibility to simultaneously analyze a great number of samples using ordinary, easily replaceable glass or plastic bottles.

Pedotransfer functions to estimate water retention parameters of soils in northeastern Brazil

Pedotransfer functions (PTF) were developed to estimate the parameters (α, n, θr and θs) of the van Genuchten model (1980) to describe soil water retention curves. The data came from various sources, mainly from studies conducted by universities in Northeast Brazil, by the Brazilian Agricultural Research Corporation (Embrapa) and by a corporation for the development of the São Francisco and Parnaíba river basins (Codevasf), totaling 786 retention curves, which were divided into two data sets: 85 % for the development of PTFs, and 15 % for testing and validation, considered independent data. Aside from the development of general PTFs for all soils together, specific PTFs were developed for the soil classes Ultisols, Oxisols, Entisols, and Alfisols by multiple regression techniques, using a stepwise procedure (forward and backward) to select the best predictors. Two types of PTFs were developed: the first included all predictors (soil density, proportions of sand, silt, clay, and organic matter), and the second only the proportions of sand, silt and clay. The evaluation of adequacy of the PTFs was based on the correlation coefficient (R) and Willmott index (d). To evaluate the PTF for the moisture content at specific pressure heads, we used the root mean square error (RMSE). The PTF-predicted retention curve is relatively poor, except for the residual water content. The inclusion of organic matter as a PTF predictor improved the prediction of parameter a of van Genuchten. The performance of soil-class-specific PTFs was not better than of the general PTF. Except for the water content of saturated soil estimated by particle size distribution, the tested models for water content prediction at specific pressure heads proved satisfactory. Predictions of water content at pressure heads more negative than -0.6 m, using a PTF considering particle size distribution, are only slightly lower than those obtained by PTFs including bulk density and organic matter content.

Incubation methods for assessing mineralizable nitrogen in soils under sugarcane

Considering nitrogen mineralization (N) of soil organic matter is a key aspect for the efficient management of N fertilizers in agricultural systems. Long-term aerobic incubation is the standard technique for calibrating the chemical extraction methods used to estimate the potentially mineralizable N in soil. However, the technique is laborious, expensive and time-consuming. In this context, the aims of this study were to determine the amount of soil mineralizable N in the 0-60 cm layer and to evaluate the use of short-term anaerobic incubation instead of long-term aerobic incubation for the estimation of net N mineralization rates in soils under sugarcane. Five soils from areas without previous N fertilization were used in the layers 0-20, 20-40 and 40-60 cm. Soil samples were aerobically incubated at 35 ºC for 32 weeks or anaerobically incubated (waterlogged) at 40 ºC for seven days to determine the net soil N mineralization. The sand, silt and clay contents were highly correlated with the indexes used for predicting mineralizable N. The 0-40 cm layer was the best sampling depth for the estimation of soil mineralizable N, while in the 40-60 cm layer net N mineralization was low in both incubation procedures. Anaerobic incubation provided reliable estimates of mineralizable N in the soil that correlated well with the indexes obtained using aerobic incubation. The inclusion of the pre-existing NH4+-N content improved the reliability of the estimate of mineralizable N obtained using anaerobic incubation.

Impacts of land leveling on lowland soil physical properties

The practice of land leveling alters the soil surface to create a uniform slope to improve land conditions for the application of all agricultural practices. The aims of this study were to evaluate the impacts of land leveling through the magnitudes, variances and spatial distributions of selected soil physical properties of a lowland area in the State of Rio Grande do Sul, Brazil; the relationships between the magnitude of cuts and/or fills and soil physical properties after the leveling process; and evaluation of the effect of leveling on the spatial distribution of the top of the B horizon in relation to the soil surface. In the 0-0.20 m layer, a 100-point geo-referenced grid covering two taxonomic soil classes was used in assessment of the following soil properties: soil particle density (Pd) and bulk density (Bd); total porosity (Tp), macroporosity (Macro) and microporosity (Micro); available water capacity (AWC); sand, silt, clay, and dispersed clay in water (Disp clay) contents; electrical conductivity (EC); and weighted average diameter of aggregates (WAD). Soil depth to the top of the B horizon was also measured before leveling. The overall effect of leveling on selected soil physical properties was evaluated by paired "t" tests. The effect on the variability of each property was evaluated through the homogeneity of variance test. The thematic maps constructed by kriging or by the inverse of the square of the distances were visually analyzed to evaluate the effect of leveling on the spatial distribution of the properties and of the top of the B horizon in relation to the soil surface. Linear regression models were fitted with the aim of evaluating the relationship between soil properties and the magnitude of cuts and fills. Leveling altered the mean value of several soil properties and the agronomic effect was negative. The mean values of Bd and Disp clay increased and Tp, Macro and Micro, WAD, AWC and EC decreased. Spatial distributions of all soil physical properties changed as a result of leveling and its effect on all soil physical properties occurred in the whole area and not specifically in the cutting or filling areas. In future designs of leveling, we recommend overlaying a cut/fill map on the map of soil depth to the top of the B horizon in order to minimize areas with shallow surface soil after leveling.

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.

Multivariate analysis applied to the study of the relationship between soil and plant properties in a peach orchard

In the State of Rio Grande do Sul, the municipality of Pelotas is responsible for 90 % of peach production due to its suitable climate and soil conditions. However, there is the need for new studies that aim at improved fruit quality and increased yield. The aim of this study was to evaluate the relationship that exists between soil physical properties and properties in the peach plant in the years 2010 and 2011 by the technique of multivariate canonical correlation. The experiment was conducted in a peach orchard located in the municipality of Morro Redondo, RS, Brazil, where an experimental grid of 101 plants was established. In a trench dug beside each one of the 101 plants, soil samples were collected to determine silt, clay, and sand contents, soil density, total porosity, macroporosity, microporosity, and volumetric water content in the 0.00-0.10 and 0.10-0.20 m layers, as well as the depth of the A horizon. In each plant and in each year, the following properties were assessed: trunk diameter, fruit size and number of fruits per plant, average weight of the fruit per plant, fruit pulp firmness, Brix content, and yield from the orchard. Exploratory analysis of the data was undertaken by descriptive statistics, and the relationships between the physical properties of the soil and of the plant were assessed by canonical correlation analysis. The results showed that the clay and microporosity variables were those that exhibited the highest coefficients of canonical cross-loading with the plant properties in the soil layers assessed, and that the variable of mean weight of the fruit per plant was that which had the highest coefficients of canonical loading within the plant group for the two years assessed.