Improving simulation of soil water balance using lysimeter observations in a semiarid climate.

Water balance simulation in cropping systems is a very useful tool to study how water can be used efficiently. However this requires that models simulate an accurate water balance. Comparing model results with field observations will provide information on the performance of the models. The objective of this study was to test the performance of DSSAT model in simulating the water balance by comparing the simulations with observed measurements. The soil water balance in DSSAT uses a one dimensional ?tipping bucket? soil water balance approach where available soil water is determined by the drained upper limit (DUL), lower limit (LL) and saturated water content (SAT). A continuous weighing lysimeter was used to get the observed values of drainage and evapotranspiration (ET). An automated agrometeorological weather station close to the lisymeter was also used to record the climatic data. The model simulated accurately the soil water content after the optimization of the soil parameters. However it was found the inability of the model to capture small changes in daily drainage and ET. For that reason simulated cumulative values had larger errors as the time passed by. These results suggested the need to compare outputs of DSSAT and some hydrological model that simulates soil water movement with a more mechanistic approach. The comparison of the two models will allow us to find which mechanism can be modified or incorporated in DSSAT model to improve the simulations.

An alternative approach for the determination of soil water mobility

A new laboratory method was proposed to establish an easily performed standard for the determination of mobile soil water close to real conditions during the infiltration and redistribution of water in a soil. It consisted of applying a water volume with a tracer ion on top of an undisturbed ring sample on a pressure plate under a known suction or pressure head. Afterwards, soil water mobility was determined by analyzing the tracer-ion concentration in the soil sample. Soil water mobility showed to be a function of the applied water volume. No relation between soil water mobility and applied pressure head could be established with data from the present experiment. A simple one- or two-parameter equation can be fitted to the experimental data to parameterize soil water mobility as a function of applied solute volume. Sandy soils showed higher mobility than loamy soils at low values of applied solute volumes, and both sandy and loamy soils showed an almost complete mobility at high applied solute volumes.

Extraction of soil water by plants: development and validation of a model

A quantitative model of water movement within the immediate vicinity of an individual root is developed and results of an experiment to validate the model are presented. The model is based on the assumption that the amount of water transpired by a plant in a certain period is replaced by an equal volume entering its root system during the same time. The model is based on the Darcy-Buckingham equation to calculate the soil water matric potential at any distance from a plant root as a function of parameters related to crop, soil and atmospheric conditions. The model output is compared against measurements of soil water depletion by rice roots monitored using γ-beam attenuation in a greenhouse of the Escola Superior de Agricultura "Luiz de Queiroz"/Universidade de São Paulo(ESALQ/USP) in Piracicaba, State of São Paulo, Brazil, in 1993. The experimental results are in agreement with the output from the model. Model simulations show that a single plant root is able to withdraw water from more than 0.1 m away within a few days. We therefore can assume that root distribution is a less important factor for soil water extraction efficiency.

Modeling water movement in horizontal columns using fractal theory

Fractal mathematics has been used to characterize water and solute transport in porous media and also to characterize and simulate porous media properties. The objective of this study was to evaluate the correlation between the soil infiltration parameters sorptivity (S) and time exponent (n) and the parameters dimension (D) and the Hurst exponent (H). For this purpose, ten horizontal columns with pure (either clay or loam) and heterogeneous porous media (clay and loam distributed in layers in the column) were simulated following the distribution of a deterministic Cantor Bar with fractal dimension H" 0.63. Horizontal water infiltration experiments were then simulated using Hydrus 2D software. The sorptivity (S) and time exponent (n) parameters of the Philip equation were estimated for each simulation, using the nonlinear regression procedure of the statistical software package SAS®. Sorptivity increased in the columns with the loam content, which was attributed to the relation of S with the capillary radius. The time exponent estimated by nonlinear regression was found to be less than the traditional value of 0.5. The fractal dimension estimated from the Hurst exponent was 17.5 % lower than the fractal dimension of the Cantor Bar used to generate the columns.

Contribution of macroporosity to water flux of a soil under different tillage systems

In view of the importance of the macroporosity for the water transport properties of soils, its quantitative assessment is a challenging task. Measurements of hydraulic conductivity (K) at different soil water tensions and the quantification of water-conducting macropores (θM) of a soil under different tillage systems could help understand the effects on the soil porous system and related hydraulic properties. The purpose of this study was to assess the effects of Conventional Tillage (CT), Chisel Plow (CP) and No Tillage (NT) on θM and on K; and to quantify the contribution of macroporosity to total water flux in a loam soil. A tension disc infiltrometer was used at two soil water pressure heads (-5 cm, and 0) to infer θM and K, during fallow. Macroporosity was determined based on the flow contribution between 0 and -5 cm water potentials (K0, K5, respectively), according to the Hagen-Poiseuille equation. The K0 values were statistically higher for CT than for NT and CP. The K5 values did not differ statistically among treatments. The mean K values varied between 0.20 and 3.70 cm/h. For CT, θM was significantly greater than for CP and NT, following the same trend as K0. No differences in θM were detected between CP and NT. With CT, the formation of water-conducting macropores with persistence until post-harvest was possible, while under CP preparation, the water-conducting macropores were not persistent. These results support the idea that tillage affects the soil water movement mainly by the resulting water-conducting macropores. Future studies on tillage effects on water movement should focus on macroporosity.

Subsurface drip irrigation emitter spacing effects on soil water redistribution, corn yield and water productivity

Emitter spacings of 0.3 to 0.6 m are commonly used for subsurface drip irrigation (SDI) of corn on the deep, silt loam soils of the U.S. Great Plains. Subsurface drip irrigation emitter spacings of 0.3, 0.6, 0.9 and 1.2 m were examined for the resulting differences in soil water redistribution, corn grain yield, yield components, seasonal water use, and water productivity in a 4‐year field study (2005 through 2008) at the Kansas State University Northwest Research‐Extension Center, Colby, Kansas. The results indicate that there is increased preferential water movement along the dripline (parallel) as compared to perpendicular to the dripline and that this phenomenon partially compensates for wider emitter spacings in terms of soil water redistribution. Corn yield and water productivity (WP) were not significantly affected by the emitter spacing with application of a full irrigation regime

Subsurface drip irrigation emitter spacing effects on soil water redistribution, corn yield and water productivity

Emitter spacings of 0.3 to 0.6 m are commonly used for subsurface drip irrigation (SDI) of corn on the deep, silt loam soils of the U.S. Great Plains. Subsurface drip irrigation emitter spacings of 0.3, 0.6, 0.9 and 1.2 m were examined for the resulting differences in soil water redistribution, corn grain yield, yield components, seasonal water use, and water productivity in a 4‐year field study (2005 through 2008) at the Kansas State University Northwest Research‐Extension Center, Colby, Kansas. The results indicate that there is increased preferential water movement along the dripline (parallel) as compared to perpendicular to the dripline and that this phenomenon partially compensates for wider emitter spacings in terms of soil water redistribution. Corn yield and water productivity (WP) were not significantly affected by the emitter spacing with application of a full irrigation regime

An easily implemented agro-hydrological procedure with dynamic root simulation for water transfer in the crop-soil system: validation and application

Models for water transfer in the crop-soil system are key components of agro-hydrological models for irrigation, fertilizer and pesticide practices. Many of the hydrological models for water transfer in the crop-soil system are either too approximate due to oversimplified algorithms or employ complex numerical schemes. In this paper we developed a simple and sufficiently accurate algorithm which can be easily adopted in agro-hydrological models for the simulation of water dynamics. We used a dual crop coefficient approach proposed by the FAO for estimating potential evaporation and transpiration, and a dynamic model for calculating relative root length distribution on a daily basis. In a small time step of 0.001 d, we implemented algorithms separately for actual evaporation, root water uptake and soil water content redistribution by decoupling these processes. The Richards equation describing soil water movement was solved using an integration strategy over the soil layers instead of complex numerical schemes. This drastically simplified the procedures of modeling soil water and led to much shorter computer codes. The validity of the proposed model was tested against data from field experiments on two contrasting soils cropped with wheat. Good agreement was achieved between measurement and simulation of soil water content in various depths collected at intervals during crop growth. This indicates that the model is satisfactory in simulating water transfer in the crop-soil system, and therefore can reliably be adopted in agro-hydrological models. Finally we demonstrated how the developed model could be used to study the effect of changes in the environment such as lowering the groundwater table caused by the construction of a motorway on crop transpiration. (c) 2009 Elsevier B.V. All rights reserved.

FRACTAL CHARACTERISTICS OF THE HORIZONTAL MOVEMENT OF WATER IN SOILS

Observed deviations from traditional concepts of soil-water movement are considered in terms of fractals. A connection is made between this movement and a Brownian motion, a random and self-affine type of fractal, to account for the soil-water diffusivity function having auxiliary time dependence for unsaturated soils. The position of a given water content is directly proportional to t(n), where t is time, and exponent n for distinctly unsaturated soil is less than the traditional 0.50. As water saturation is approached, n approaches 0.50. Macroscopic fractional Brownian motion is associated with n < 0.50, but shifts to regular Brownian motion for n = 0.50.

Invariant Solutions of Richards' Equation for Water Movement in Dissimilar Soils

Scaling methods allow a single solution to Richards' equation (RE) to suffice for numerous specific cases of water flow in unsaturated soils. During the past half-century, many such methods were developed for similar soils. In this paper, a new method is proposed for scaling RE for a wide range of dissimilar soils. Exponential-power (EP) functions are used to reduce the dependence of the scaled RE on the soil hydraulic properties. To evaluate the proposed method, the scaled RE was solved numerically considering two test cases: infiltration into relatively dry soils having initially uniform water content distributions, and gravity-dominant drainage occurring from initially wet soil profiles. Although the results for four texturally different soils ranging from sand to heavy clay (adopted from the UNSODA database) showed that the scaled solution were invariant for a wide range of flow conditions, slight deviations were observed when the soil profile was initially wet in the infiltration case or deeply wet in the drainage case. The invariance of the scaled RE makes it possible to generalize a single solution of RE to many dissimilar soils and conditions. Such a procedure reduces the numerical calculations and provides additional opportunities for solving the highly nonlinear RE for unsaturated water flow in soils.

Eucalyptus sp seedling response to potassium fertilization and soil water

A considerable portion of Brazil's commercial eucalypt plantations is located in areas Subjected to periods of water deficit and grown in soils with low natural fertility, particularly poor In potassium. Potassium is influential in controlling water relations of plants. The objective of this study was to verify the influence of potassium fertilization and soil water potential (psi(w)) oil the dry matter production and oil water relations Of eucalypt seedlings grown under greenhouse conditions. the experimental units were arranged in 4x4x2 randomized blocks factorial design, as follow: four species of Eucalyptus (Eucalyptus grandis, Eucalyptus urophylla, Eucalyptus camaldulensis and hybrid Eucalyptus grandis x Eucalyptus urophylla), four dosages of K (0, 50, 100 and 200 mg dm(-3)) and two soil water potentials (-0.01 M Pa and -0.1 M Pa). Plastic containers with 15 cm diameter and 18 cm height, with Styrofoam base, containing 3.0 dm(3) of soil and two plants per container were used. Soil water potential was kept at -0.01 MPa for 40 days after seeding. Afterward, the experimental units were divided into two groups: in one group the potential was kept at 0.01 MPa, and in the other one, at -0.10 MPa. Sol I water potential was control led gravimetrically twice a day with water replacement until the desired potential was reestablished. A week before harvesting, the leaf water potential (psi), the photosynthetic rate (A), the stomatal conductance (gs) and the transpiration rate were evaluated. The last week before harvesting, the mass of the containers was recorded daily before watering to determine the consumption of water by the plants. After harvesting, total dry matter and leaf area were evaluated. the data were Submitted to analysis of variance, to Tukey's tests and regression analyses. The application of K influenced A, gs and the transpiration rate. Plants deficient in K showed lower A and higher Us and transpiration rates. There were no statistical differences in A, gs and transpiration rates ill plants with and Without water deficit. The addition of K reduced the consumption of water per unit of leaf area and, in general, plants submitted to water deficit presented a lower consumption of water.

Watershed scale temporal stability of soil water content

The recognition of temporally stable locations with respect to soil water content is of importance for soil water management decisions, especially in sloping land of watersheds. Neutron probe soil water content (0 to 0.8 m), evaluated at 20 dates during a year in the Loess Plateau of China, in a 20 ha watershed dominated by Ust-Sandiic Entisols and Aeolian sandy soils, were used to define their temporal stability through two indices: the standard deviation of relative difference (SDRD) and the mean absolute bias error (MABE). Specific concerns were (a) the relationship of temporal stability with soil depth, (b) the effects of soil texture and land use on temporal stability, and (c) the spatial pattern of the temporal stability. Results showed that temporal stability of soil water content at 0.2 m was significantly weaker than those at the soil depths of 0.6 and 0.8 m. Soil texture can significantly (P<0.05) affect the stability of soil water content except for the existence of an insignificant difference between sandy loam and silt loam textures, while temporal stability of areas covered by bunge needlegrass land was not significantly different from those covered by korshinsk peashrub. Geostatistical analysis showed that the temporal stability was spatially variable in an organized way as inferred by the degree of spatial dependence index. With increasing soil depth, the range of both temporal stability indices showed an increasing trend, being 65.8-120.5 m for SDRD and 148.8-214.1 m for MABE, respectively. This study provides a valuable support for soil water content measurements for soil water management and hydrological applications on sloping land areas. (C) 2010 Elsevier B.V. All rights reserved.

Using a New Criterion to Identify Sites for Mean Soil Water Storage Evaluation

Establishing a few sites in which measurements of soil water storage (SWS) are time stable significantly reduces the efforts involved in determining average values of SWS. This study aimed to apply a new criterion the mean absolute bias error (MABE)-to identify temporally stable sites for mean SWS evaluation. The performance of MABE was compared with that of the commonly used criterion, the standard deviation of relative difference (SDRD). From October 2004 to October 2008, SWS of four soil layers (0-1.0, 1.0-2.0,2.0-3.0, and 3.0-4.0 m) was measured, using a neutron probe, at 28 sites on a hillslope of the Loess Plateau, China. A total of 37 SWS data sets taken over time were divided into two subsets, the first consisting of 22 dates collected during the calibration period from October 2004 to September 2006, and the second with 15 dates collected during the validation period from October 2006 to October 2008. The results showed that if a critical value of 5% for MABE was defined, more than half the sites were temporally stable for both periods, and the number of temporally stable sires generally increased with soil depth. Compared with SDRD, MABE was more suitable for the identification of time-stable sites for mean SS prediction. Since the absolute prediction error of drier sites is more sensitive to changes in relative difference in terms of mean SWS prediction, the sites of wet sectors should be preferable for mean SWS prediction for the same changes in relative difference.

Impacts of land-use pattern on soil water-content variability on the Loess Plateau of China

Long-term vegetation restoration carried out on the slopes of the Loess Plateau of China employed different spatial and temporal land-use patterns but very little is known about the effects of these patterns on soil water-content variability. For this study the small Donggou catchment was selected to investigate soil water-content distributions for three spatial scales, including the entire catchment area, sampling transects, and land-use systems. Gravimetric soil water contents were determined incrementally to a soil depth of 1.20 m, on 10 occasions from April to October, 2007, at approximately 20-day intervals. Results indicated that soil water contents were affected by the six land-use types, resulting in four distinct patterns of vertical distribution of soil moisture (uniform, increasing, decreasing, and fluctuating with soil depth). The soil water content and its variation were also influenced in a complex manner by five land-use patterns distributed along transects following the gradients of five similar slopes. These patterns with contrasting hydrological responses in different components, such as forage land (alfalfa)-cropland-shrubland or shrubland-grassland (bunge needlegrass)-cropland-grassland, showed the highest soil water-content variability. Soil water at the catchment scale exhibited a moderate variability for each measurement date, and the variability of soil water content decreased exponentially with increasing soil water content. The minimum sample size for accurate data for use in a hydrological model for the catchment, for example, required many more samples for drier (69) than for wet (10) conditions. To enhance erosion and runoff control, this study suggested two strategies for land management: (i) to create a mosaic pattern by land-use arrangement that located units with higher infiltration capacities downslope from those with lower soil infiltrabilities; and (ii) raising the soil-infiltration capacity of units within the spatial mosaic pattern where possible.

Time stability of soil water storage measured by neutron probe and the effects of calibration procedures in a small watershed

The knowledge of soil water storage (SWS) of soil profiles is crucial for the adoption of vegetation restoration practices. With the aim of identifying representative sites to obtain the mean SWS of a watershed, a time stability analysis of neutron probe evaluations of SWS was performed by the means of relative differences and Spearman rank correlation coefficients. At the same time, the effects of different neutron probe calibration procedures were explored on time stability analysis. mean SWS estimation. and preservation of the spatial variability of SWS. The selected watershed, with deep gullies and undulating slopes which cover an area of 20 ha, is characterized by an Ust-Sandiic Entisol and an Aeolian sandy soil. The dominant vegetation species are bunge needlegrass (Stipa bungeana Trim) and korshinsk peashrub (Carugano Korshinskii kom.). From June 11, 2007 to July 23,2008, SWS of the top1 m soil layer was evaluated for 20 dates, based on neutron probe data of 12 sampling sites. Three calibration procedures were employed: type 1, most complete, with each site having its own linear calibration equation (TrE); type II. with TrE equations extended over the whole field: and type III, with one single linear calibration curve for the whole field (UnE) and also correcting its intercept based on site specific relative difference analysis (RdE) and on linear fitting of data (RcE), both maintaining the same slope. A strong time stability of SWS estimated by TrE equations was identified. Soil particle size and soil organic matter content were recognized as the influencing factors for spatial variability of SWS. Land use influenced neither the spatial variability nor the time stability of SWS. Time stability analysis identified one site to represent the mean SWS of the whole watershed with mean absolute percentage errors of less than 10%, therefore. this site can be used as a predictor for the mean SWS of the watershed. Some equations of type II were found to be unsatisfactory to yield reliable mean SWS values or in preserving the associated soil spatial variability. Hence, it is recommended to be cautious in extending calibration equations to other sites since they might not consider the field variability. For the equations with corrected intercept (type III), which consider the spatial variability of calibration in a different way in relation to TrE, it was found that they can yield satisfactory means and standard deviation of SWS, except for the RdE equations, which largely leveled off the SWS values in the watershed. Correlation analysis showed that the neutron probe calibration was linked to soil bulk density and to organic matter content. Therefore, spatial variability of soil properties should be taken into account during the process of neutron probe calibration. This study provides useful information on the mean SWS observation with a time stable site and on distinct neutron probe calibration procedures, and it should be extended to soil water management studies with neutron probes, e.g., the process of vegetation restoration in wider area and soil types of the Loess Plateau in China. (C) 2009 Elsevier B.V. All rights reserved.