Transient three-dimensional modeling of soil water and solute transport with simultaneous root growth,root water and nutrient uptake

A three-dimensional solute transport model was developed and linked to a three-dimensional transient model for soil water flow and root growth. The simulation domain is discretized into a grid of finite elements by which the soil physical properties are spatially distributed. Solute transport modeling includes passive and active nutrient uptake by roots as well as zero- and first-order source/sink terms. Root water uptake modeling accounts for matric and osmotic potential effects on water and passive nutrient uptake. Root age effects on root water and nutrient uptake activity have been included, as well as the influence of nutrient deficiency and ion toxicity on root growth. Examples illustrate simulations with different levels of model complexity, depending on the amount of information available to the user. At the simplest level, root growth is simulated as a function of mechanical soil strength only. Application of the intermediate level with root water and nutrient uptake simulates the influence of timing and amount of NO3 application on leaching. The most comprehensive level includes simulation of root and shoot growth as influenced by soil water and nutrient status, temperature, and dynamic allocation of assimilate to root and shoot.     

Drying of mucilage causes water repellency in the rhizosphere of maize: measurements and modelling

Background and Aims

Although maize roots have been extensively studied, there is limited information on the effect of root exudates on the hydraulic properties of maize rhizosphere. Recent experiments suggested that the mucilaginous fraction of root exudates may cause water repellency of the rhizosphere. Our objectives were: 1) to investigate whether maize rhizosphere turns hydrophobic after drying and subsequent rewetting; 2) to test whether maize mucilage is hydrophobic; and 3) to find a quantitative relation between rhizosphere rewetting, particle size, soil matric potential and mucilage concentration.

Methods

Maize plants were grown in aluminum containers filled with a sandy soil. When the plants were 3-weeks-old, the soil was let dry and then it was irrigated. The soil water content during irrigation was imaged using neutron radiography. In a parallel experiment, ten maize plants were grown in sandy soil for 5 weeks. Mucilage was collected from young brace roots growing above the soil. Mucilage was placed on glass slides and let dry. The contact angle was measured with the sessile drop method for varying mucilage concentration. Additionally, capillary rise experiments were performed in soils of varying particle size mixed with maize mucilage. We then used a pore-network model in which mucilage was randomly distributed in a cubic lattice. The general idea was that rewetting of a pore is impeded when the concentration of mucilage on the pore surface (g cm^?2) is higher than a given threshold value. The threshold value depended on soil matric potential, pore radius and contract angle. Then, we randomly distributed mucilage in the pore network and we calculated the percolation of water across a cubic lattice for varying soil particle size, mucilage concentration and matric potential.

Results

Our results showed that: 1) the rhizosphere of maize stayed temporarily dry after irrigation; 2) mucilage became water repellent after drying. Mucilage contact angle increased with mucilage surface concentration (gram of dry mucilage per surface area); 3) Water could easily cross the rhizosphere when the mucilage concentration was below a given threshold. In contrast, above a critical mucilage concentration water could not flow through the rhizosphere. The critical mucilage concentration decreased with increasing particle size and decreasing matric potential.

Conclusions

These results show the importance of mucilage exudation for the water fluxes across the root-soil interface. Our percolation model predicts at what mucilage concentration the rhizosphere turns hydrophobic depending on soil texture and matric potential. Further studies are needed to extend these results to varying soil conditions and to upscale them to the entire root system.

     

The gradient between pre-dawn rhizoplane and bulk soil matric potentials, and its relation to the pre-dawn root and leaf water potentials of four species

Uptake of soil water by plants may result in significant gradients between bulk soil and soil in the vicinity of roots. Few experimental studies of water potential gradients in close proximity to roots, and no studies on the relationship of water potential gradients to the root and leaf water potentials, have been conducted. The occurrence and importance of pre-dawn gradients in the soil and their relation to the pre-dawn root and leaf water potentials were investigated with seedlings of four species. Pre-germinated seeds were grown without watering for 7 and lid in a silt loam soil with initial soil matric potentials of -0.02, -0.1 and -0.22 MPa. Significant gradients, independent of the species, were observed only at pre-dawn soil matric potentials lower than -0.25 MPa; the initial soil matric potentials were -0.1 MPa. At an initial bulk soil matric potential of -0.22 MPa, a steep gradient between bulk and rhizoplane soil was observed after 7 d for maize (Zea mays L. cv. Issa) and sunflower (Helianthus annuus L. cv. Nanus), in contrast to barley (Hordeum vulgare L. cv. Athos) and wheat (Triticum aestivum L. cv. Kolibri). Pre-dawn root water potentials were usually about the same as the bulk soil matric potential and were higher than the rhizoplane soil matric potential. Pre-dawn root and leaf water potentials tended to be much higher than rhizoplane soil matric potentials when the latter were lower than -0.5 MPa. It is concluded that plants tend to become equilibrated overnight with the wetter bulk soil or with wetter zones in the bulk soil. Plants can thus circumvent negative effects of localized steep pre-dawn soil matric potential gradients. This may be of considerable importance for water uptake and growth in drying soil.     

Modeling soil water movement with water uptake by roots

Soil water movement with root water uptake is a key process for plant growth and transport of water and chemicals in the soil-plant system. In this study, a root water extraction model was developed to incorporate the effect of soil water deficit and plant root distributions on plant transpiration of annual crops. For several annual crops, normalized root density distribution functions were established to characterize the relative distributions of root density at different growth stages. The ratio of actual to potential cumulative transpiration was used to determine plant leaf area index under water stress from measurements of plant leaf area index at optimal soil water condition. The root water uptake model was implemented in a numerical model. The numerical model was applied to simulate soil water movement with root water uptake and simulation results were compared with field experimental data. The simulated soil matric potential, soil water content and cumulative evapotranspiration had reasonable agreement with the measured data. Potentially the numerical model implemented with the root water extraction model is a useful tool to study various problems related to flow transport with plant water uptake in variably saturated soils. This revised version was published online in June 2006 with corrections to the Cover Date.     

Effect of soil strength on root growth under different water conditions

The objective of this study was to determine the effects of soil water and soil strength on root growth in situations where the individual effects of both of these factors were important. Three grain legumes were grown from pre-germinated seeds for five days on 50-mm compacted columns of two major soils of Sri Lanka. Four or five levels of bulk density (1.1 to 1.8 Mg.m^–3) and five or six levels of matric potential (–0.02 to–2.0 MPa) were used.Soil strength and matric potential effects on root growth were independently significant for most crop and soil combinations. Under high (wet) matric potential (>–0.77 MPa) soil conditions, the effect of soil water on root growth was evident only in its effect on soil strength. Bulk density had a significant effect on root growth independent of soil strength and matric potential in three cases.For all crops and soils, root penetration was 80% of the maximum or greater when the average soil strength (soil water not limiting) was 0.75 MPa or less, and when the average matric potential (soil strength not limiting) was –0.77 MPa or greater (wetter). Root penetration was 20% of the maximum or less when the soil strength was greater than 3.30 MPa (soil water not limiting), and when matric potential (soil strength not limiting) was less than –3.57 MPa. The use of pre-germinated seeds, which contained imbibed water, combined with a lack of water loss from the closed chambers containing the plants is the probable cause for the very low (–3.57 MPa) matric potential that allowed root growth at 20% of the maximum.     

Microtensiometer technique for in situ measurement of soil matric potential and root water extraction from a sandy soil

The suitability of microtensiometers to measure the spatial variation of soil matric potential and its diurnal change was tested in a pot experiment with pearl millet (Pennisetum americanum [L.] Leeke) in a sandy soil as the soil dried out.The temporal and spatial resolution of this technique allowed precise measurement of soil matric potential and thus estimation of soil water extraction from different compartments as well as from the whole rooting zone. The technique also allowed the measurement of rehydration of plants at night and root water uptake rate per unit soil volume or per unit root length. The precision of determination of root water uptake depended greatly on the accuracy of the estimate of hydraulic conductivity, which was derived from a bare soil and might be different for a cropped soil owing to aggregation induced by the root system. A linear relationship between root length and water uptake was found (r²=0.82), irrespective of variation in soil water content between compartments and despite the variation in root age, xylem differentiation and suberin formation expected to exist between different compartments of the rooting zone. As the experiment was carried out in a range of soil matric potentials between –4 and –30 kPa, drought stress did not occur. Further information at lower soil matric potentials are required, to address questions such as the importance of soil resistance for water uptake, or which portion of the root system has to be stressed to induce hormonal signals to the shoot. The microtensiometer technique can be applied to soil matric potentials up to –80 kPa.     

Response of microbial activity and community structure to decreasing soil osmotic and matric potential

Low soil water content (low matric potential) and salinity (low osmotic potential) occur frequently in soils, particularly in arid and semi-arid regions. Although the effect of low matric or low osmotic potential on soil microorganisms have been studied before, this is the first report which compares the effect of the two stresses on microbial activity and community structure. A sand and a sandy loam, differing in pore size distribution, nutrient content and microbial biomass and community structure, were used. For the osmotic stress experiment, salt (NaCl) was added to achieve osmotic potentials from ?0.99 to ?13.13 MPa (sand) and from ?0.21 to 3.41 MPa (sandy loam) after which the soils were pre-incubated at optimal water content for 10d. For the matric stress experiment, soils were also pre-incubated at optimal water content for 10d, after which the water content was adjusted to give matric potentials from ?0.03 and ?1.68 MPa (sand) and from ?0.10 to 1.46 MPa (sandy loam). After amendment with 2% (w/w) pea straw (C/N 26), soil respiration was measured over 14d. Osmotic potential decreased with decreasing soil water content, particularly in the sand. Soil respiration decreased with decreasing water potential (osmotic?+?matric). At a given water potential, respiration decreased to a greater extent in the matric stress experiment than in the osmotic stress experiment. Decreasing osmotic and matric potential reduced microbial biomass (sum of phospholipid fatty acids measured after 14 days) and changed microbial community structure: fungi were less tolerant to decreasing osmotic potential than bacteria, but more tolerant to decreasing water content. It is concluded that low matric potential may be more detrimental than a corresponding low osmotic potential at optimal soil water content. This is likely to be a consequence of the restricted diffusion of substrates and thus a reduced ability of the microbes to synthesise osmolytes to help maintain cell water content. The study also highlighted that it needs to be considered that decreasing soil water content concentrates the salts, hence microorganisms in dry soils are exposed to two stressors.     

Topographic position affects the water regime in a semideciduous tropical forest in Panamá

The effects of topographic position on water regime in a semideciduous tropical forest on Barro Colorado Island in Panamá were assessed by measuring soil matric potential using the filter paper technique and by using measured soil water release characteristics to convert a long-term (20 years) gravimetric water content data-set to matric potential. These were also compared against predictions from a simple water balance model. Soil matric potentials on slope sites were significantly higher than on plateau sites throughout the measurement interval and slopes experienced a shorter duration of drought during the annual dry-season. Measured values of matric potential agreed with those predicted from converting the gravimetric measurements using water release characteristics. Annual duration of drought predicted by the simple water balance model agreed with values determined from the converted long term water content data-set and was able to predict the annual duration of drought on plateau sites. On slope sites, the water balance systematically and significantly overestimated the duration of drought obtained from the water content data-set, suggesting that slope sites were supplied with water from upslope. Predictions of annual drought duration from sites with higher annual rainfall than Barro Colorado Island (BCI), suggest that while plateau sites on BCI experience a water regime consistent with annual rainfall, slopes experience a water regime more similar to that of forests with much higher rainfall. We conclude that such large variations in water regime over small spatial scales may play a role in maintaining high species richness through providing opportunities for niche specialisation and by buffering slopes against possible climate change.     

Seasonal and long-term soil water regime in West African tropical forest

Abstract. Variation of soil matric potential of a Wet Evergreen and a Moist Semi-deciduous West African forest were compared. The two forest types differed strongly in their soil water regime. Wet Evergreen forest experienced matric potentials below ?100 kPa only occasionally, while in Moist Semi-deciduous forest matric potentials were less than ?2.5 MPa for periods of several weeks or more each season. A water balance equation was used to simulate the soil water regime at both sites and predict severity and length of the dry period. The predictions showed good agreement with the field measurements of soil water potential over a 2-yr period. The methodology was used to estimate the occurrence and severity of droughts over longer periods at the two sites. The balance calculations suggest that droughts occur occasionally in the Wet Evergreen forest under study. The potential impact of droughts on species distribution and vegetation disturbance in tropical forests is discussed.     

A coupled model of stomatal conductance, photosynthesis and transpiration

A model that couples stomatal conductance, photosynthesis, leaf energy balance and transport of water through the soil–plant–atmosphere continuum is presented. Stomatal conductance in the model depends on light, temperature and intercellular CO₂ concentration via photosynthesis and on leaf water potential, which in turn is a function of soil water potential, the rate of water flow through the soil and plant, and on xylem hydraulic resistance. Water transport from soil to roots is simulated through solution of Richards’ equation. The model captures the observed hysteresis in diurnal variations in stomatal conductance, assimilation rate and transpiration for plant canopies. Hysteresis arises because atmospheric demand for water from the leaves typically peaks in mid‐afternoon and because of uneven distribution of soil matric potentials with distance from the roots. Potentials at the root surfaces are lower than in the bulk soil, and once soil water supply starts to limit transpiration, root potentials are substantially less negative in the morning than in the afternoon. This leads to higher stomatal conductances, CO₂ assimilation and transpiration in the morning compared to later in the day. Stomatal conductance is sensitive to soil and plant hydraulic properties and to root length density only after approximately 10 d of soil drying, when supply of water by the soil to the roots becomes limiting. High atmospheric demand causes transpiration rates, LE, to decline at a slightly higher soil water content, θ_s, than at low atmospheric demand, but all curves of LE versus θ_s fall on the same line when soil water supply limits transpiration. Stomatal conductance cannot be modelled in isolation, but must be fully coupled with models of photosynthesis/respiration and the transport of water from soil, through roots, stems and leaves to the atmosphere.     

Bioengineering ground improvement considering root water uptake model

Bioengineering features of native vegetation are currently being evolved to enhance soil stiffness, slope stabilisation and erosion control. The effects of tree roots on soil moisture content and ground settlement are discussed in this paper. Matric suction induced by tree roots is a key factor, governing the properties of unsaturated soils, directly imparting stability to slopes and resistance for yielding behaviour. A mathematical model for the rate of root water uptake that considers ground conditions, type of vegetation and climatic parameters has been developed. This study highlights the inter-related parameters contributing to the development of a conceptual evapo-transpiration and root moisture uptake equilibrium model that is then incorporated in a comprehensive numerical finite element model. The developed model considers fully coupled-flow-deformation behaviour of soil. Field measurements obtained by the Authors from a site in Victoria, South of Australia, are used to validate the model. In this study, the active tree root distribution has been predicted by measuring soil organic content distribution. The predicted results show acceptable agreement with the field data in spite of the assumptions made for simplifying the effects of soil heterogeneity and anisotropy. The results prove that the proposed root water uptake model can reliably predict the region of the maximum matric suction away from the tree axis.     

Parametric studies on bioengineering effects of tree root-based suction on ground behaviour

Using native vegetation to improve soil stiffness, stabilise slopes and control erosion is a rapidly evolving process. A theoretical model previously developed by the authors for the rate of tree root water uptake together with an associated numerical simulation is used to study the effects of a wide range of soil, tree, and atmospheric parameters on partially saturated ground. The influence of different parameters on the maximum initial rate of root water uptake is investigated through parametric and sensitivity analyses. Field measurements taken from previously published literature are compared with numerical predictions for validation. The rate of selected parameters such as potential transpiration and its distribution, suction at wilting point, the coefficient of permeability and the distribution of root length density are studied in detail. The analysis shows that the rate of potential transpiration increases the soil matric suction and ground settlement, while the potential transpiration rate has an insignificant effect on the distribution of soil suction. Root density distribution factors affect the size of the influence zone. Suction at the wilting point increases the soil matric suction and ground settlement, whereas the saturation permeability decreases the maximum soil matric suction generated. The analysis confirms that the most sensitive parameters, including the coefficients of the tree root system, the transpiration rate, the permeability of the soil and its suction at the wilting point should be measured or estimated accurately for an acceptable prediction of ground conditions in the vicinity of trees.     

A highly conductive drainage extension to control the lower boundary condition of lysimeters

Lysimeters are used to study and monitor water, fertilizers, salts and other contaminants and are particularly valuable in transpiration and evapotranspiration research. Saturation at the soil bottom boundary in a lysimeter is inherent to its design. A drainage extension made of porous media with high hydraulic conductivity and substantial water holding capacity was devised to extend the lysimeter in order to produce soil moisture conditions mimicking those in the field. Design criteria that assure equal discharge in the soil and in the highly conductive drain (HCD) were established and formulated. Desired matric head at the lysimeter base is determined by HCD extension length. Its value can be manipulated and can range between saturation and the soil's field capacity. Conditions where the HCD is not limiting to flow are obtained through selection of the appropriate cross sectional area ratio between the soil in the lysimeter and the HCD. The validity of these criteria was confirmed with 200 l working lysimeters in the field, with and without plants, and with detailed flow tests utilizing smaller (15 l) lysimeters. Comparison of computed and measured matric head and leachate volume indicates that the proposed method can serve to maintain conditions similar to those in the field.     

Refining and unifying the upper limits of the least limiting water range using soil and plant properties

The least limiting water range (LLWR) was introduced as an integrated soil water content indicator, measuring the impact of mechanical impedance, oxygen and water availability on water uptake and crop growth. However, a rigorous definition of the upper limit of the LLWR using plant physiological and soil physical concepts was not given. We introduce in this study an upper limit of the LLWR, based on soil physical and plant physiological properties. We further evaluate the sensitivity of this boundary to different soil and crop variables, and compare the sensitivity of the upper limit of the LLWR to previous definitions of soil water content at field capacity. The current study confirms that the upper limit of the LLWR can be predicted from knowledge of the soil moisture characteristic curve, plant root depth and oxygen consumption rate. The sensitivity analysis shows further that the upper limit of the LLWR approaches the volumetric soil water content at saturation when the oxygen consumption rate by plants becomes less than 2 µmol m^?3 s^?1. When plants are susceptible to aeration (e.g. potato and avocado), there is a big difference between the upper limit of the LLWR and the soil water content at field capacity, in particular for sandy soils. Results also show that the soil water content at aeration porosity corresponding to 10% cannot be considered as an appropriate upper limit of LLWR because it does not appropriately reflect the crop water requirements. Similar poor results are obtained when considering the soil water content at matric potential ?0.033 MPa or when defining the soil water content at field capacity based on drainage flux rate. It is observed that the upper limit of the LLWR is higher than either soil water content at ?0.033 MPa matric potential or soil water content at field capacity as based on drainage flux rate, especially in sandy soils.     

Calibration of Watermark soil moisture sensors for soil matric potential and temperature

Rapid, accurate, and automated measurement of soil matric potential is desirable. Evidence suggested that the Watermark resistance block might be an appropriate and inexpensive tool, so we conducted an evaluation of its relevant characteristics. A number of these blocks were calibrated under laboratory conditions to determine their individual and aggregate responses to soil matric potential, soil type, and temperature. We found that the temperature response could be expressed as a single equation, valid for all tested blocks, but comparison against matric potential revealed that each block had a characteristic response. Furthermore, block responses were different in two soils and, for a given soil, not necessarily reproducible. Given these limitations, these sensors are probably useful only as relative indicators of soil water status.     

Leaf surface wetness in sorghum and resistance to shoot fly, Atherigona soccata: role of soil and plant water potentials

In experiments with potted plants, the relationships between soil matric potential, plant water potential and production of water droplets (leaf surface wetness) on the folded central whorl leaf of seedlings of sorghum genotypes that are either resistant or susceptible to shoot fly (Atherigona soccata) damage were investigated. Differences in soil matric potentials in the pots affected the plant water status, which in turn had profound effects on the production of water droplets on the central whorl leaf of the sorghum genotype susceptible to shoot fly. There was no consistent variation in the relationship between plant water potential and soil matric potential of resistant and susceptible sorghum genotypes. However, there was very little or practically no water droplets on the central whorl leaf of the resistant genotypes, indicating that the production of water droplets is not solely the result of internal water status of the plant. It is suggested that leaf surface wetness is genetically controlled and that an understanding of the mechanism by which water is transferred to the leaf surface will enhance breeding for resistance to shoot fly.     

The effect of water stress on the mortality ofBetula pendula Roth. andBuddleia davidii Franch. seedlings

Summary A simple container is described whereby small seedlings may be grown at controlled levels of water stress. The water stress was induced in the soil by an osmoticum which is separated from the soil by a semi-permeable membrane. The mortality ofBetula pendula seedlings was markedly increased at a matric potential of –1.6 bars whereas the mortality ofBuddleia davidii was only affected below –2.8 bars. This difference in tolerance to water stress at the seedling stage might not be reflected in the distribution of the species in the colonisation of chalk and sand pits in England unless there is a dry spring.     

土壤-植物-大气连续体水流阻力分布规律的研究

本文依据田间实测资料,分析了土壤-植物-大气连续体水流阻力的相对重要性,结果表明在连续体中的水流阻力主要分布于从叶气孔腔到大气的扩散过程和根系的吸水过程。叶-气之间的水流阻力比土-根之间要大50倍。最后,讨论了控制连续体水流运动的气孔阻力的变化规律及其与环境因素之间的关系。     

Colonization of Wheat Roots by an Exopolysaccharide-Producing Pantoea agglomerans Strain and Its Effect on Rhizosphere Soil Aggregation

The effect of bacterial secretion of an exopolysaccharide (EPS) on rhizosphere soil physical properties was investigated by inoculating strain NAS206, which was isolated from the rhizosphere of wheat (Triticum durum L.) growing in a Moroccan vertisol and was identified as Pantoea aglomerans. Phenotypic identification of this strain with the Biotype-100 system was confirmed by amplified ribosomal DNA restriction analysis. After inoculation of wheat seedlings with strain NAS206, colonization increased at the rhizoplane and in root-adhering soil (RAS) but not in bulk soil. Colonization further increased under relatively dry conditions (20% soil water content; matric potential, −0.55 MPa). By means of genetic fingerprinting using enterobacterial repetitive intergenic consensus PCR, we were able to verify that colonies counted as strain NAS206 on agar plates descended from inoculated strain NAS206. The intense colonization of the wheat rhizosphere by these EPS-producing bacteria was associated with significant soil aggregation, as shown by increased ratios of RAS dry mass to root tissue (RT) dry mass (RAS/RT) and the improved water stability of adhering soil aggregates. The maximum effect of strain NAS206 on both the RAS/RT ratio and aggregate stability was measured at 24% average soil water content (matric potential, −0.20 MPa). Inoculated strain NAS206 improved RAS macroporosity (pore diameter, 10 to 30 μm) compared to the noninoculated control, particularly when the soil was nearly water saturated (matric potential, −0.05 MPa). Our results suggest that P. agglomerans NAS206 can play an important role in the regulation of the water content (excess or deficit) of the rhizosphere of wheat by improving soil aggregation.     

Hydrological and biogeochemical controls on the timing and magnitude of nitrous oxide flux across an agricultural landscape

Anticipated increases in precipitation intensity due to climate change may affect hydrological controls on soil N₂O fluxes, resulting in a feedback between climate change and soil greenhouse gas emissions. We evaluated soil hydrologic controls on N₂O emissions during experimental water table fluctuations in large, intact soil columns amended with 100 kg ha^?1 KNO₃‐N. Soil columns were collected from three landscape positions that vary in hydrological and biogeochemical properties (N= 12 columns). We flooded columns from bottom to surface to simulate water table fluctuations that are typical for this site, and expected to increase given future climate change scenarios. After the soil was saturated to the surface, we allowed the columns to drain freely while monitoring volumetric soil water content, matric potential and N₂O emissions over 96 h. Across all landscape positions and replicate soil columns, there was a positive linear relationship between total soil N and the log of cumulative N₂O emissions (r²= 0.47; P= 0.013). Within individual soil columns, N₂O flux was a Gaussian function of water‐filled pore space (WFPS) during drainage (mean r²= 0.90). However, instantaneous maximum N₂O flux rates did not occur at a consistent WFPS, ranging from 63% to 98% WFPS across landscape positions and replicate soil columns. In contrast, instantaneous maximum N₂O flux rates occurred within a narrow range (?1.88 to ?4.48 kPa) of soil matric potential that approximated field capacity. The relatively consistent relationship between maximum N₂O flux rates and matric potential indicates that water filled pore size is an important factor affecting soil N₂O fluxes. These data demonstrate that matric potential is the strongest predictor of the timing of N₂O fluxes across soils that differ in texture, structure and bulk density.