Macroscopic root water uptake distribution using a matric flux potential approach

Hydrological models featuring root water uptake usually do not include compensation mechanisms such that reductions in uptake from dry layers are compensated by an increase in uptake from wetter layers. We developed a physically based root water uptake model with an implicit compensation mechanism. Based on an expression for the matric flux potential (M) as a function of the distance to the root, and assuming a depth-independent value of M at the root surface, uptake per layer is shown to be a function of layer bulk M, root surface M, and a weighting factor that depends on root length density and root radius. Actual transpiration can be calculated from the sum of layer uptake rates. The proposed reduction function (PRF) was built into the SWAP model, and predictions were compared to those made with the Feddes reduction function (FRF). Simulation results were tested against data from Canada (continuous spring wheat [(Triticum aestivum L.]) and Germany (spring wheat, winter barley [Hordeum vulgare L.], sugarbeet [Beta vulgaris L.], winter wheat rotation). For the Canadian data, the root mean square error of prediction (RMSEP) for water content in the upper soil layers was very similar for FRF and PRF; for the deeper layers, RMSEP was smaller for PRF. For the German data, RMSEP was lower for PRF in the upper layers and was similar for both models in the deeper layers. In conclusion, but dependent on the properties of the data sets available for testing,the incorporation of the new reduction function into SWAP was successful, providing new capabilities for simulating compensated root water uptake without increasing the number of input parameters or degrading model performance.

Measurement systems of soil water matric potential and evaluation of soil moisture under different irrigation depths

The development of new methodologies and tools that enable to determine the water content in soil is of fundamental importance to the practice of irrigation. The objective of this study was to evaluate soil matric potential using mercury tensiometer and puncture digital tensiometer, and to compare the gravimetric soil moisture values obtained by tensiometric system with gravimetric soil moisture obtained by neutron attenuation technique. Four experimental plots were maintained with different soil moisture by irrigation. Three repetitions of each type of tensiometer were installed at 0.20 m depth. Based on the soil matric potential and the soil water retention curve, the corresponding gravimetric soil moisture was determined. The data was then compared to those obtained by neutron attenuation technique. The results showed that both tensiometric methods showed no difference under soil matric potential higher than -40 kPa. However, under drier soil, when the water was replaced by irrigation, the soil matric potential of the puncture digital tensiometer was less than those of the mercury tensiometer.

Contributions of matric and osmotic potentials to the unfrozen water content of frozen soils

Recent reports show that biogeochemical processes continue when the soil is frozen, but are limited by water availability. However, there is little knowledge about the interactive effects of soil and environmental variables on amounts of unfrozen water in frozen soils. The aims of this study were to determine the contributions of matric and osmotic potentials to the unfrozen water content of frozen soil. We determined the effects of matric and osmotic potential on unfrozen water contents of frozen mineral soil fractions (ranging from coarse sand to fine silt) at -7 degrees C, and estimated the contributions of these potentials to liquid water contents in samples from organic surface layers of boreal soils frozen at -4 degrees C. In the mineral soil fractions the unfrozen water contents appeared to be governed solely by the osmotic potential, but in the humus layers of the sampled boreal soils both the osmotic and matric potentials control unfrozen water content, with osmotic potential contributing 20 to 69% of the total water potential. We also determined pore size equivalents, where unfrozen water resides at -4 degrees C, and found a strong correlation between these equivalents and microbial CO2 production. The larger the pores in which the unfrozen water is found the larger the microbial activity that can be sustained. The osmotic potential may therefore be a key determinant of unfrozen water and carbon dynamics in frozen soil. (C) 2008 Elsevier B.V. All rights reserved.

Measurement systems of soil water matric potential and evaluation of soil moisture under different irrigation depths

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Effect of solute and matric potential on growth rate of fungal species isolated from cheese

The effect of water potential ( J w ) on the growth of 15 fungal species isolated from cheeses was analysed. The species, identified mainly by analysis of DNA sequences, belonged to genera Penicillium , Geotrichum , Mucor , Aspergillus , Microascus and Talaromyces . Particularly, the effect of matric potential ( J m ), and ionic (NaCl) and non-ionic (glycerol) solute potentials ( J s ) on growth rate was studied. The response of strains was highly dependent on the type of J w . For J s , clear profiles for optimal, permissive and marginal conditions for growth were obtained, and differences in growth rate were achieved comparing NaCl and glycerol for most of the species. Conversely, a sustained growth was obtained for J m in all the strains, with the exception of Aspergillus pseudoglaucus , whose growth increased proportionally to the level of water stress. Our results might help to understand the impact of environmental factors on the ecophysiology and dynamics of fungal populations associated to cheeses.

Effect of solute and matric potential on growth rate of fungal species isolated from cheese

The effect of water potential ( J w ) on the growth of 15 fungal species isolated from cheeses was analysed. The species, identi fi ed mainly by analysis of DNA sequences, belonged to genera Penicillium, Geotrichum, Mucor , Aspergillus , Microascus and Talaromyces . Particularly, the effect of matric potential ( J m ), and ionic (NaCl) and non-ionic (glycerol) solute potentials ( J s ) on growth rate was studied. The response of strains was highly dependent on the type of J w . For J s, clear profiles for optimal, permissive and marginal conditions for growth were obtained, and differences in growth rate were achieved comparing NaCl and glycerol for most of the species. Conversely, a sustained growth was obtained for J m in all the strains, with the exception of Aspergillus pseudoglaucus, whose growth increased proportionally to the level of water stress. Our results might help to understand the impact of environmental factors on the ecophysiology and dynamics of fungal populations associated to cheeses.

Unsaturated Infinite Slope Stability Considering Surface Flux Conditions

A slope stability model is derived for an infinite slope subjected to unsaturated infiltration flow above a phreatic surface. Closed form steady state solutions are derived for the matric suction and degree of saturation profiles. Soil unit weight, consistent with the degree of saturation profile, is also directly calculated and introduced into the analyzes, resulting in closed-form solutions for typical soil parameters and an infinite series solution for arbitrary soil parameters. The solutions are coupled with the infinite slope stability equations to establish a fully realized safety factor function. In general, consideration of soil suction results in higher factor of safety. The increase in shear strength due to the inclusion of soil suction is analogous to making an addition to the cohesion, which, of course, increases the factor of safety against sliding. However, for cohesive soils, the results show lower safety factors for slip surfaces approaching the phreatic surface compared to those produced by common safety factor calculations. The lower factor of safety is due to the increased soil unit weight considered in the matric suction model but not usually accounted for in practice wherein the soil is treated as dry above the phreatic surface. The developed model is verified with a published case study, correctly predicting stability under dry conditions and correctly predicting failure for a particular storm.

New Analytic Solution of Boltzmann Transform for Horizontal Water Infiltration into Sand

The use of the Boltzmann transform function, lambda(theta), to solve the Richards equation when the diffusivity, D, is a function of only soil water content,., is now commonplace in the literature. Nevertheless, a new analytic solution of the Boltzmann transform lambda(h) as a function of matric potential for horizontal water infiltration into a sand was derived without invoking the concept or use of D(theta). The derivation assumes that a similarity exists between the soil water retention function and the Boltzmann transform lambda(theta). The solution successfully described soil water content profiles experimentally measured for different infiltration times into a homogeneous sand and agrees with those presented by Philip in 1955 and 1957. The applicability of this solution for all soils remains open, but it is anticipated to hold for soils whose air-filled pore-size distribution before wetting is sufficiently narrow to yield a sharp increase of water content at the wetting front during infiltration. It also improves and provides a versatile alternative to the well-known analysis pioneered by Green and Ampt in 1911.

Pupal Development of Ceratitis capitata (Diptera: Tephritidae) and Diachasmimorpha longicaudata (Hymenoptera: Braconidae) at Different Moisture Values in Four Soil Types

This study aimed to evaluate adult emergence and duration of the pupal stage of the Mediterranean fruit fly, Ceratitis capitata (Wiedemann), and emergence of the fruit fly parasitoid, Diachasmimorpha longicaudata (Ashmead), under different moisture conditions in four soil types, using soil water matric potential Pupal stage duration in C capitata was influenced differently for males and females In females, only soil type affected pupal stage duration, which was longer in a clay soil In males, pupal stage duration was individually influenced by moisture and soil type, with a reduction in pupal stage duration in a heavy clay soil and in a sandy clay, with longer duration in the clay soil As allude potential decreased, duration of the pupal stage of C capitata males increased, regardless of soil type C capitata emergence was affected by moisture, regardless of soil type, and was higher in drier soils The emergence of D longicaudata adults was individually influenced by soil type and moisture factors, and the number of emerged D longicaudata adults was three times higher in sandy loam and lower in a heavy clay soil Always, the number of emerged adults was higher at higher moisture conditions C capitata and D longicaudata pupal development was affected by moisture and soil type, which may facilitate pest sampling and allow release areas for the parasitoid to be defined under field conditions.

Root Water Extraction under Combined Water and Osmotic Stress

Using a numerical implicit model for root water extraction by a single root in a symmetric radial flow problem, based on the Richards equation and the combined convection-dispersion equation, we investigated some aspects of the response of root water uptake to combined water and osmotic stress. The model implicitly incorporates the effect of simultaneous pressure head and osmotic head on root water uptake, and does not require additional assumptions (additive or multiplicative) to derive the combined effect of water and salt stress. Simulation results showed that relative transpiration equals relative matric flux potential, which is defined as the matric flux potential calculated with an osmotic pressure head-dependent lower bound of integration, divided by the matric flux potential at the onset of limiting hydraulic conditions. In the falling rate phase, the osmotic head near the root surface was shown to increase in time due to decreasing root water extraction rates, causing a more gradual decline of relative transpiration than with water stress alone. Results furthermore show that osmotic stress effects on uptake depend on pressure head or water content, allowing a refinement of the approach in which fixed reduction factors based on the electrical conductivity of the saturated soil solution extract are used. One of the consequences is that osmotic stress is predicted to occur in situations not predicted by the saturation extract analysis approach. It is also shown that this way of combining salinity and water as stressors yields results that are different from a purely multiplicative approach. An analytical steady state solution is presented to calculate the solute content at the root surface, and compared with the outputs of the numerical model. Using the analytical solution, a method has been developed to estimate relative transpiration as a function of system parameters, which are often already used in vadose zone models: potential transpiration rate, root length density, minimum root surface pressure head, and soil theta-h and K-h functions.

Validation of a root water uptake model to estimate transpiration constraints

Experimental results obtained from a greenhouse trial with common bean (Phaseolus vulgaris L) plants performed to test model hypotheses regarding the onset of limiting hydraulic conditions and the shape of the transpiration reduction curve in the falling rate phase are presented. According to these hypotheses based on simulations with an upscaled single-root model, the matric flux potential at the onset of limiting hydraulic conditions is as a function of root length density and potential transpiration rate, while the relative transpiration in the falling rate phase equals the relative matric flux potential. Transpiration of bean plants in water stressed pots with four different soils was determined daily by weighing and compared to values obtained from non-stressed pots. This procedure allowed determining the onset of the falling rate phase and corresponding soil hydraulic conditions. At the onset of the falling rate phase, the value of matric flux potential M(I) showed to differ in order of magnitude from the model predicted value for three out of four soils. This difference between model and experiment can be explained by the heterogeneity of the root distribution which is not considered by the model. An empirical factor to deal with this heterogeneity should be included in the model to improve predictions. Comparing the predictions of relative transpiration in the falling rate phase using a linear shape with water content, pressure head or matric flux potential, the matric flux potential based reduction function, in agreement with the hypothesis, showed the best performance, while the pressure head based equation resulted in the highest deviations between observed and predicted values of relative transpiration rates. (C) 2010 Elsevier B.V. All rights reserved.

Modeling of transpiration reduction in van Genuchten-Mualem type soils

We derive an analytic expression for the matric flux potential (M) for van Genuchten-Mualem (VGM) type soils which can also be written in terms of a converging infinite series. Considering the first four terms of this series, the accuracy of the approximation was verified by comparing it to values of M estimated by numerical finite difference integration. Using values of the parameters for three soils from different texture classes, the proposed four-term approximation showed an almost perfect match with the numerical solution, except for effective saturations higher than 0.9. Including more terms reduced the discrepancy but also increased the complexity of the equation. The four-term equation can be used for most applications. Cases with special interest in nearly saturated soils should include more terms from the infinite series. A transpiration reduction function for use with the VGM equations is derived by combining the derived expression for M with a root water extraction model. The shape of the resulting reduction function and its dependency on the derivative of the soil hydraulic diffusivity D with respect to the soil water content theta is discussed. Positive and negative values of dD/d theta yield concave and convex or S-shaped reduction functions, respectively. On the basis of three data sets, the hydraulic properties of virtually all soils yield concave reduction curves. Such curves based solely on soil hydraulic properties do not account for the complex interactions between shoot growth, root growth, and water availability.

Evaluating responses of maize (Zea mays L.) to soil physical conditions using a boundary line approach

The functional relation between the decline in the rate of a physiological process and the magnitude of a stress related to soil physical conditions is an important tool for uses as diverse as assessment of the stress-related sensitivity of different plant cultivars and characterization of soil structure. Two of the most pervasive sources of stress are soil resistance to root penetration (SR) and matric potential (psi). However, the assessment of these sources of stress on physiological processes in different soils can be complicated by other sources of stress and by the strong relation between SR and psi in a soil. A multivariate boundary line approach was assessed as a means of reducing these cornplications. The effects of SR and psi stress conditions on plant responses were examined under growth chamber conditions. Maize plants (Zea mays L.) were grown in soils at different water contents and having different structures arising from variation in texture, organic carbon content and soil compaction. Measurements of carbon exchange (CE), leaf transpiration (ILT), plant transpiration (PT), leaf area (LA), leaf + shoot dry weight (LSDW), root total length (RTL), root surface area (RSA) and root dry weight (RDW) were determined after plants reached the 12-leaf stage. The LT, PT and LA were described as a function of SR and psi with a double S-shaped function using the multivariate boundary line approach. The CE and LSDW were described by the combination of an S-shaped function for SR and a linear function for psi. The root parameters were described by a single S-shaped function for SR. The sensitivity to SR and psi depended on the plant parameter. Values of PT, LA and LSDW were most sensitive to SR. Among those parameters exhibiting a significant response to psi, PT was most sensitive. The boundary line approach was found to be a useful tool to describe the functional relation between the decline in the rate of a physiological process and the magnitude of a stress related to soil physical conditions. (C) 2009 Elsevier B.V. All rights reserved.