Soil Water Diffusivity Research Articles

Measurement and estimation of evapotranspiration in a maize field: A new method based on an analytical water flux model

Quantifying evapotranspiration (ET) in rainfed cropping systems can be challenging due to complicated interactions among site-specific soil, plant, and management factors. In Northeast China, ET and soil water status in maize fields often display strong spatial and temporal variations due to the changes in tillage practice, planting pattern, and maize plant density. Previous studies have shown that near-surface soil water content (θ) observations at multiple scales provide the potential to estimate surface soil water fluxes. In this study, we introduced a new method to estimate daily field ET by using a soil water flux model mainly based on the time-series of θ at a depth of 2.5 cm. The new method required a calibration of soil water diffusivity with maximum net water flux in the near-surface soil layer, which was related to precipitation redistribution below the canopy. Finally, the new method was evaluated using observed ET values over a 2-year period in a maize field, where independent measurements of soil water evaporation (E) and transpiration (T) were made with heat-pulse sensors and sap-flow gauges, respectively. Field observations showed that E dominated water loss during the seedling stage (16% of total ET). As the canopy was fully developed, E sharply decreased to a value of 0.4 mm d−1, and T accounted for about 89% of ET since the silking stage. The new method to estimate ET performed well in drying periods, while it tended to underestimate ET in wet periods with substantial infiltration into the surface layer. On rain-free days, the ET values estimated with the new method matched well with the measured E+T values, with R2 and RMSE values of 0.85 and 1.93 mm d−1. Therefore, the new approach provides an effective way to quantify maize ET.

Prediction of the weighted-mean soil water diffusivity in supporting the falling rate stage of evaporation

Daily soil evaporation from drying soil can be described as isothermal linear diffusion process, resulting in a simple equation for estimating soil evaporation rate. However, determining soil water diffusivity and subsequently the crucial parameter of weighted-mean diffusivity, essential for estimating soil evaporation rate, is challenging and time-consuming. In this paper, the soil water diffusivity was re-evaluated by considering the effect of film flow, a phenomenon often overlooked but recently argued by Wang et al. (2019, https://doi.org/10.1029/2019WR025003) as the dominant process in the falling rate stage of evaporation. Theoretical derivation reveals a universal relationship between soil water diffusivity and matric potential, independent of soil types. Additionally, soil water diffusivity can be expressed in the form of an exponential function, a form commonly assumed empirically in the literature. Testing with various datasets suggests that the theoretically derived soil water diffusivity generally agrees closely with observations for matric potentials less than about −1 m, wherein the diffusivity accounting for capillary flow shows significant underestimation. Furthermore, a simple function is derived for estimating the weighted-mean diffusivity, enabling a direct estimation of soil evaporation rates in the falling rate stage. Testing with various data sets of different soil types and ambient conditions generally demonstrates close agreements between the estimated soil evaporation rates and observations, indicating that the developed method provides a robust estimation of the weighted-mean soil water diffusivity. When the initial matric potential of the linear diffusion process of drying is higher than about −1 m, which usually occurs under very high atmospheric demand conditions or for very coarse-textured soils, however, the effects of capillary flow were found to be important and need consideration. The effects of vapor diffusion were also discussed and found to have a negligible impact on the estimation of weighted-mean soil water diffusivity.

Effect of unsaturated soil water transport function linearization for mulched drip irrigation

Water flow in the unsaturated zone is an important component of the water cycle. Linearized methods have been widely adopted to solve unsaturated flow problems described with Richards’ equation, which is nonlinear. To better understand the effect of unsaturated soil water transport function linearization, this article presents an analytical study of a two-dimensional soil water flow system. Analytical solutions, experimental observations and numerical simulations were compared, and the results indicated that the soil water content variation trends were the same, but the values were a little different. The analytical solutions have a certain reference value for researching unsaturated water flow. However, the soil volumetric water content (SWC) results were sensitive to linearized water flow parameter variations in soil water diffusivity (D). The D changes with the SWC. The results showed that unsaturated water flow function linearization affected the analytical solution results essentially. As the nonlinear math problem was not solved, we recommend the numerical solution when more precise values are needed and complex conditions are considered.

Effects of urea solution concentration on soil hydraulic properties and water infiltration capacity

Elucidating the effect of fertigation on soil hydraulic parameters and water-solute transportation is fundamental to the design of farmland irrigation systems and their sustainable utilization. Few studies have focused on soil hydraulic parameters or water infiltration characteristics or how they are influenced by urea solution concentration. In this study, the clay loam and sandy loam in Yangling District of Shaanxi Province, China, were used as test soil, and experiments involving seven urea solution concentrations (0.2, 0.4, 0.6, 0.8, 1, 3, and 5 g/L) and a control treatment (0 g/L) were conducted to explore the influence of the various urea solution concentrations on soil hydraulic parameters and water infiltration characteristics. The results indicated that the cumulative infiltration and wetting front migration depth increased with urea solution concentration, as accurately estimated using the Kostiakov model and a power function, respectively. In addition, the coefficients of the Kostiakov model and the power function increased with urea solution concentration. Treatment with multiple concentrations of urea solution resulted in an increase in the volume of macro pores in the soil but a reduction in the volume of mesopores and micro pores in the soil, leading to increases in the saturated water content, saturated hydraulic conductivity, soil water diffusivity, and infiltration capacity and a reduction in the water-holding capacity of the soil. The effect of urea solute potential on the inhibition of soil water movement is small, and this inhibitory effect is far weaker than the improvement effect of the urea solution on soil structure, and hence enhance the soil water infiltration capacity. Our results increase the understanding of soil hydrological mechanisms and may be usefully applied for improving the management of fertigation.

ANALYTICAL SOLUTION OF UNSATURATED SOIL WATER FLOW FROM A POINT SOURCE

Water flow into unsaturated porous media is governed by the Richards’ partial differential equation expressing the mass conservation and Darcy’s laws. The Richards’ equation may be written in three forms,where the dependent variable is pressure head or moisture content, and the constitutive relationships between water content and pressure head allow for conversion of one form into the other. In the present paper, the “moisture-based" form of Richards’ equation is linearized by applying Kirchhoff’s transformation, whichcombines the soil water diffusivity and soil water content. Then the similarity method is used to obtain the analytical solution of wetting front position. This exact solution is obtained by means of Lie’s method of infinitesimal transformation groups. The predicted results of the analytical solution agreed well with available results of experiments and numerical solutions.

Soil water diffusivity as function of the pore size distribution and precompression stress

The idea of the study was to combine the saturation-dependent hydraulic conductivity, K(Se), as intensity parameter with the volume fraction related to the fractional capillary potential for each of the characteristic pore size classes that based on the matric flow potential, φ, and the soil water diffusivity, D(θ). The analysis is focused on the soil hydraulic properties of differently textured soils with differences in the pore size distribution, PSD. The soil-moisture retention, θ(h), and K(Se) as function of the effective saturation were determined with undisturbed samples for each soil profile and horizon and fitted by Mualem–van Genuchten type analytical function. The results indicate a strong linear relationship between φ values and the soil water diffusivity, D(θ), of soils with its volume fractions of macropores (r2: 0.65–0.91), wide coarse pores (r2: 0.55–0.77), and medium pores (r2: 0.83–0.85). The precompression stress, Pc, as mechanical soil strength parameter affects the soil water diffusivity of the chosen datasets especially in the wide coarse pore range (r2: 0.85). The results suggest that the PSD’s-related φ and D(θ) values are affected by both textural and structural effects. The soil water diffusivity itself strongly depends on the PSD and is an effective link between soil’s intensity and capacity parameter.

Towards the influences of three types of biocrusts on soil water in drylands: Insights from horizontal infiltration and soil water retention

As a “living skin” of soil in drylands, biocrusts possibly change the retention and movement of soil water, and thus critically influence hydrological and erosion processes as well as soil moisture regime. Contrary to vertical infiltration, little attention was paid to horizontal infiltration, which may provide important information on soil water diffusivity and sorptivity. In this study, the one-dimensional horizontal infiltration experiments for bare soil and three types of biocrusts (cyanobacterial, cyanobacterial-moss mixed, and moss crusts) were conducted with disturbed soil samples in a semiarid dryland. Undisturbed soil samples were also taken for each treatment, and their soil water retention capacity were determined in the laboratory. The substantial differences that were found were analyzed, and the differences were further validated through monitoring and comparing their in situ soil moisture regimes at depth of 0–20 cm. Our results showed that the three types of biocrusts had significantly lower constant infiltration rate (0.11 vs. 0.18 cm min−1) and infiltration water amount (6.75 vs. 10.01 cm) than bare soil, and their wetting front moved much slower (0.60 vs. 1.14 cm min−1) in comparison to bare soil, with bare soil > cyanobacterial crusts > cyanobacterial-moss mixed crusts > moss crusts. Total sorptivity of cyanobacterial, mixed, and moss crusts were 21%, 32%, and 53% lower than that of bare soil, respectively. More importantly, at the same level of water content, all crust types had lower water diffusivity, sorptivity, and conductivity than the bare soil, and exhibited a higher water content than bare soil at a same level of matric potential (across the range of −15,000–0 hPa), indicating an increase in water-holding capacity and water availability of the biocrusts in contrast to bare soil. Furthermore, the in-situ soil water content of all biocrusts increased on average by 30% at depth of 0–10 cm but decreased by 23% at 10–20 cm. Confirming the lower infiltrability and higher water-holding capacity of the biocrusts in comparison to bare soil, biocrusts increase surface soil water retention capacity and water availability at the surface, thus playing a vital role in reshaping the surface soil water balance, which subsequently may affect surface hydropedological and biochemical processes in drylands.

Anisotropy of soil water diffusivity of hillslope soil under spruce forest derived by X-ray CT and lab experiments

The idea of the study is to indicate direction-dependent differences in hydraulic conductivity, K(Se), and soil water diffusivity, D(θ), as function of the volume fraction related to the fractional capillary potential for each of the characteristic pore size classes by extended anisotropy factors. The study is exemplary focused on a BwC horizon of a Dystric Cambisol under spruce forest formed on the weathered and fractured granite bedrock in the mountainous hillslopes Uhlirska catchment (Czech Republic). Thus, undisturbed soil samples were taken in vertical (0°, y = x-axis) and horizontal (90°, z-axis) direction. The D(θ) values and especially the D(θ)-weighted anisotropy ratios showed that anisotropy increases with the volume fraction of macropores, MaP (d > 0.03 mm), with r2 between 0.89 and 0.92. The X-ray computer tomography (CT) based anisotropy ratio (ACT) is larger for the horizontal sampled soil core with 0.31 than for the vertical with 0.09. This underlines the existence of a predominantly horizontally oriented pore network and the fact that weathered bedrock strata can initiate lateral preferential flow. The study results suggest that combining the hydraulic conductivity as intensity and the capacity parameter by means of diffusivity results in an extended anisotropy ratio which unveils the role of the soil hydraulic characteristics in generation of small-scale lateral preferential flow. In future, the small-scale direction-dependent differences in the soil hydraulic capacity and intensity parameter will be used for model-based upscaling for better understanding of preferential flow at the catchment scale.

A new method of estimating soil water diffusivity based on variation in the wetting front position and on the principle of stationary action

This paper presents a novel method for the estimation of soil water diffusivity that solves the nongravitational Richards equation using the variation in the boundary of the advancing wetting front. The solution is based on the principle of stationary action. The position of the wetting front is determined from the continuity of water content at the wetting front interface, as the underlying physical principles of this area are well known and thus it can be described using relatively few assumptions. This also eliminates the need for the Boltzmann-transformation, which is commonly required in traditional methods. The Brooks-Corey empirical water diffusion function is used to analytically determine the significant parameters of diffusivity. The soil water diffusivities of 15 soil samples matched well with their measured soil–water diffusivities, thereby confirming the reliability of the method. All of the measurements required for using this method are rapidly and easily obtained, and the experimental apparatus is simple and inexpensive. Compared with other methods, this novel method is more accurate and convenient for use in laboratories for measuring the soil hydraulic properties.

SOIL HYDRAULIC PROPERTIES OF A CHROMIC LUVISOL IN KATUMANI, KENYA

<p><strong>Background:</strong> Soil hydraulic parameters in non-saturated conditions are crucial for explaining soil water dynamics in the field. It is therefore necessary to understand the link between soil water potential and hydraulic conductivity in the soil in order to estimate plant available water and hence simulate its movement within the soils. However, measurement of such hydraulic properties in the field and laboratory is somehow difficult, laborious and costly. <strong>Objective</strong>: To determine soil hydraulic properties of Ferro-chromic Luvisols in Katumani using the RETC code based on pedo-transfer functions; % sand, silt, clay and soil bulk density. <strong>Methodology:</strong> Undisturbed soil samples were collected from a profile pit at 0-15, 16-30, 31-45 and 46-60 cm depths using core rings for bulk densities and texture determination. Soil water retention curves and saturated hydraulic conductivities (Ksat) were estimated for all the samples using standard suction apparatus and the constant head method, respectively. The air entry suction (α) and pore size distribution (n) were generated using the RETC code. <strong>Results:</strong> The permanent wilting point and field capacity were at 0.081, 0.102, 0.107 and 0.121 and 0.188, 0.225, 0.241, 0.262 m<sup>3</sup>m<sup>-3 </sup>H<sub>2</sub>O, whilst its soil water diffusivity ranged from a low of 6.39, 6.94, 9.03 to a high of 12.5 cm<sup>2</sup>min<sup>-1</sup> in the 0-15, 16-30, 31-45 and 46-60 cm depth, respectively. Ksat values from RETC code ranged from 29 - 48 cm day<sup>-1</sup>, while the total and readily available water within the soil profile were 330.4 and 214.7 mm H<sub>2</sub>O, respectively. <strong>Implication:</strong> The air entry value (α) and pore size distribution (n) implied an almost even distribution from the top and subsequent horizons. The soils ‘field capacity was achieved at pF 2.0 whilst PWP was arbitrary indicated at pF 4.2 reducing the time it takes to calculate irrigation cycles based on the amount of water available to the crops. <strong>Conclusion:</strong> The data indicates that pedo-transfer functions; especially high bulk densities negatively impact on soil hydraulics conductivity.</p>

Physico-chemical and thermal characteristics of sandy loam soils contaminated by single and mixed pollutants (mineral and vegetable oils).

Samples of contaminated top soil (0-30cm) and uncontaminated soil (control) from two locations in Precambrian basement complex area were analyzed to assess the effects of single and mixed oil contaminants on the physico-chemical and thermal properties of soils. Pearson's correlation and analysis of variance (ANOVA) were used to study the interrelationships of the studied parameters as well as variation of studied soil characteristics under the different oil contaminants, respectively. Results showed insignificant impact of pollutant(s) on the textural class of contaminated soils. The highest and lowest mean soil resistivity (SR) values were found in petrol-contaminated and mixed surfactants (shampoo + conditioner) -contaminated soils, respectively. The least values of mean specific heat capacity (SHC), heat capacity (HC), and soil water diffusivity (SWD) were found in soils contaminated by mixed surfactants-contaminants. However, mixed mineral oils (petrol + diesel + engine oil) and mixed vegetable oils (palm oil + groundnut oil) had mean SHC values lower than those of control soils at the two sampling locations. The mixed surfactants-polluted soil is characterized by lowest mean bulk density (BD) and highest mean gravimetric water content (GWC) while lowest mean GWC and highest mean SWD characterized engine oil-contaminated soils. The analysis of variance (ANOVA) result revealed significant variation in % sand at 5% level (p < 0.05) for petroleum-derived contaminants but no significant differences in mean values of all analyzed soil properties under vegetable oil contaminants.

Determination of soil water hydraulic parameters from infiltration data

PurposeTraditional laboratory measurements of soil water diffusivity (D) and soil water retention curve (SWRC) are always time-consuming and labor-intensive. Therefore, this paper aims to present a simple and robust test method for determining D and SWRC without reducing accuracy.Design/methodology/approachIn this study, a D model of unsaturated soil was established based on Gardner–Russo model and then a combination of Gardner–Russo model with one-dimensional horizontal absorption method to obtain n and a parameters of Gardner–Russo model. One-dimensional horizontal absorption experiments on loam, silt loam and sandy clay loam were conducted to obtain the relationships between measured infiltration rate and cumulative infiltration with wetting front distance. Based on the obtained relationships, the measured infiltration data from the one-dimensional horizontal absorption tests were used to calculate n and a parameters and further constructing D and SWRC.FindingsBoth the calculated D and SWRC inversed from the infiltration data were in good agreement with the measured ones that obtained from the traditional horizontal absorption method and the centrifuge method, respectively. Error analysis indicated that only the infiltration data are enough to reliably synchronously determine D and SWRC.Originality/valueA simple and robust method is proposed for synchronous determination of soil water diffusivity and water retention curve.

Transition Between Constitutive Equations and the Mechanics of Water Flow in Unsaturated Soil: Numerical Simulations

Richards’ partial differential equation normally governs water flows in soil. It can be utilized to investigate water infiltration in the soil. A case study of the Haverkamp’s water infiltration simulation into Yolo Light Clay was carried out. A common practice to change between the hydraulic functions (from Haverkamp to van Genuchten) could cause a discrepancy in the simulation results. Hence, the first objective is to modify van Genuchten’s equations to reproduce the water infiltration result that was obtained using the Haverkamp constitutive functions. The method used was able to recreate water infiltration by increasing the fitting parameters of the van Genuchten constitutive functions. The second objective is to identify the mechanism governing the flow of water. The soil water diffusivity and the hydraulic conductivity were responsible for the flux of water in the soil. The Richards’ equation was discretized using a finite difference method, and the algebraic solution was coded into Simply Fortran 2008.

Anisotropy of unsaturated hydraulic properties of compacted mineral capping systems seven years after construction

The mechanical compaction of soil material of mineral landfill systems affects the continuity and connectivity of the complex soil pore network. A horizontally-oriented layering is intended to generate a slope-induced lateral water flow out of mineral capping systems that is sufficient to minimise the statutory required vertical percolation through the underlying waste body and the potential leachate.In this case, soil compaction affects both the porosity and water retention as capacity values and the hydraulic conductivity as intensity parameter. The idea of this study was to combine information on both soil properties in an extended anisotropy factor based on the soil water diffusivity. The analysis is focused on the direction-dependent soil hydraulic properties of a mechanically compacted landfill capping system.In particular, the volume fractions were related to the fractional capillary potential for each of the characteristic pore size classes. Three different soil profiles of top, middle, and bottom slopes of the mineral capping system of the Rastorf landfill in Northern Germany were sampled seven years after construction. Undisturbed soil cores of 100 cm3 and 438 cm3 were extracted in vertical (ver) and horizontal directions (hor) in depths of 20, 50, and 80 cm representing the main layers. The soil water retention and unsaturated hydraulic conductivity, K, functions were determined by suctions plates, permeameter, and the evaporation method. In the coarse pores range (pressure head values of h ≥ −300 hPa), the standard anisotropy ratio, AR, (K(Se)hor/K(Se)ver) as a function of effective saturation, Se, in the sealing liner in 80 cm depth was larger than 1, indicating higher horizontal than vertical K(Se) values. Thus, AR-values above 1 in the range close to water saturation especially in 80 cm depth suggest the tendency of lateral water flow out of mineral capping system due to a sufficient hydraulic potential and thus its reasonable functionality, even seven years after construction.The anisotropy factor was extended in two steps; for AR* and AR**, the pore size class-related matric flux potential, ϕ, and the soil water diffusivity, D(θ), were proposed to combine intensity parameters with capacity-based volume fractions of pore size classes and the fractional capillary potential. The ϕ- and D(θ)-weighted anisotropy ratios, AR* and AR**, indicate that anisotropy increases with the volume fraction of macropores (r2AR* of 0.69−0.77; r2AR** of 0.71−0.80) and wide coarse pores (r2AR* of 0.57−0.78; r2AR** of 0.79−0.89) in both directions. The results suggest that by combining both the intensity and the capacity parameters of the soil hydraulic properties in an extended anisotropy ratio improves the information on compacted mineral capping systems.

Methodology to simulate unsaturated zone hydrology in Storm Water Management Model (SWMM) for green infrastructure design and evaluation.

Hydrologic models such as the USEPA Stormwater Management Model (SWMM) are commonly used to assess the design and performance of green infrastructure (GI). To accurately represent GI performance models used in design need to be able to address both the hydrology/hydraulics of the catchment and the GI unsaturated (vadose) zone hydrology. While hydrologic models, such as SWMM, address the need for catchment hydrology/hydraulics, they often simplify the unsaturated zone hydrology. This paper presents a methodology utilizing existing components of SWMM to represent unsaturated zone hydrology in an accessible format that does not require adjustments to the SWMM source code. The methodology simulated the unsaturated soil water movement by considering flow caused by differences of soil matric head and flow caused by gravity between soil layers with finite depth/length. The flow flux related to the soil matric head is a function of soil water diffusivity (D) and the soil moisture gradient, where D can be represented by a pump curve in SWMM. The flow flux related to gravity was controlled by unsaturated hydraulic conductivity (K) only and was also simulated by a pump. The methodology was compared to another variably saturated model, HYDRUS, with theoretical soils (with single layers of sand, loam, silt, and clay, as well as dual-layer scenarios). Field data was used to compare the methodology to HYDRUS and the SWMM LID (Low Impact Development) module. In all comparisons the presented methodology and HYDRUS delivered similar results for the vadose zone response to a storm event, while the LID module of SWMM exhibited slower water movement. The results showed that under natural conditions, the approximation of the presented methodology yielded satisfactory results to simulate flow through the unsaturated vadose zone.

Quantifying the hydraulic properties of some Egyptian soils using RETC code

Temporal and spatial soil variabilities reduce the accuracy of quantifying the hydraulic properties of the soil, leading to poor irrigation management. Modeling estimation and computer codes (e.g., RETC code) have been used to quantify and predict soil hydraulic properties. One hundred soil samples were collected to represent soil textural classes according to USDA textural triangle. Following the physical, chemical and hydraulic characterization of the soil samples collected, one soil sample was chosen to represent each texture class. The selected soil samples covered six USDA textural classes. Then, Brooks-Corey (B-C) (1964) and van Genuchten (vG) (1980) parametric models were used to describe the functional relationship between soil water tension and water content, i.e. the soil water retention curve, and the Mualem (1976) theoretical models of pore size distribution were used to predict the functions of unsaturated conductivity K (θ) and soil water diffusivity D (θ). The RETC (retention curve) code was used to obtain the parameters of vG and B-C models by simultaneously matching to the observed data. Output file of refining process showed the iteration levels to achieve minimum value of residual summed square (RSSQ).The values of determination coefficient (R2) of the fitted θ (h) as a power function increased after fitting the experimental data either in B-C or in vG models. The values of R2 of the fitted equation of θ (h) with both B-C and vG data increased after refining process through RETC code. The increments of R2 values after refining are more pronounced in coarse textured samples (Sand, Loamy sand, and Sandy loam).The values of determination coefficient (R2) of the fitted functions K (θ) and D (θ) as a power function are highly significant after fitting soil parameters either in B-C or in vG models. While using RETC code refined the obtained soil parameters of both K (θ) and D(θ) increasing R2 values of fitted power function and reducing residual summation square. RETC reduced the value of residual summed square (RSSQ) of the objective function O(b) under using both B-C and vG models. This reduction of RSSQ by using RETC fluctuated between 60.4% and 98.2% with B-C for six soil textural classes. While, the reduction of RSSQ by using RETC fluctuated between 90.4% and 98.6% with vG for the same six soil textural classes. In general, the RETC code is a good tool for obtaining accurate values of hydraulic properties in a variety of soil textures.

Investigation of the Flux–Concentration Relation for Horizontal Flow in Soils

The objective of the present work is to investigate the flux–concentration (F(Θ)) relation, where Θ is the normalized soil volumetric water content for the case of one-dimensional horizontal flow, subject to constant concentration conditions. More specifically, the possibility of describing F(Θ) by an equation of the form F(Θ) = 1 − (1 − Θ)p+1 is examined. Parameter p is estimated from curve-fitting of the experimentally obtained λ(Θ) data to an analytic expression of the form (1 − Θ)p where λ is the well-known Boltzmann transformation λ = xt−0.5 (x = distance, t = time). The results show that the equation of (1 − Θ)p form can satisfactorily describe the λ(Θ) relation for the four porous media tested. The proposed F(Θ) function was compared with the limiting F(Θ) function for linear and Green–Ampt soils and to the actual F(Θ) function. From the results, it was shown that the proposed F(Θ) function gave reasonably accurate results in all cases. Moreover, the analytical expression of the soil water diffusivity (D(Θ)) function, as it was obtained by using the equation for λ(Θ) of the form (1 − Θ)p, appears to be very close to the experimental D(Θ) data (root mean square error (RMSE) = 0.593 m2min−1).

Effect of Organic Matter and Initial Moisture Content on Water Transmission Characteristics of Alfisols of Assam

A study was conducted to study the effect of organic matter and initial moisture content on water transmission behaviour of three texturally different Alfisols. Saturated hydraulic conductivity showed a variable increase with increase in organic matter. Soil water diffusivity D (θ) showed an irregular trend with increasing levels of organic matter while sorptivity (s) increased with increasing levels of organic matter, penetrability (P) of wetting front decreased. Weighted mean diffusivity did not show any particular trend amongst the treatments, however, it was found to be higher in treated soils as compared to untreated control. Specific water capacity increased with increasing levels of organic matter, while capillary conductivity (k) did not show any definite trend. All the water transmission parameters except sorptivity were found to be increased with an increase in initial moisture content. Diffusivity was in general higher at lower initial moisture content indicating that soil water diffusivity was greatly influenced by initial moisture content. Prediction of D (θ) was, by and large, satisfactory in coarse textured soils and inconclusive in case of K (θ) relation in most of the soils.

Diffusivity and sorptivity determination at different soil water contents from horizontal infiltration

Determination of soil hydraulic properties under non-saturated conditions is crucial in order to describe soil water dynamics in the field. The objectives of this work were: i- to determine soil water diffusivity (D) as a function of soil moisture with a simple lab set-up, based on Whisler et al. (1968) methodology, from soil moisture profiles obtained using capacitive sensors in three contrasting soils from Argentinean Pampas Region; and ii- to obtain a function that allows to determine sorptivity (S) values in the whole moisture range from diffusivity data. Additionally, horizontal infiltration data was used to obtain hydraulic conductivity function (K (h)) through inverse parameterization using HYDRUS 1-D. Horizontal infiltration was performed in order to determine D and S at different soil water content in 15 repacked columns from three different textured soils (loam, silty loam and sandy loam). The results indicate that the new approach is a suitable methodology for D (θ) and S (θ) determination. D values measured using the proposed laboratory setup ranged between 1 × 10−5 and 1.67 cm2 s−1, showing an exponential dependence with soil water content (r2 > 0.85). S values for different soil water contents ranged between 0.078 and 0.14 cm s−1/2. The highest value of S corresponded to the lowest soil water content, remaining practically constant in a large range of soil water content, and abruptly decreasing near saturation toward zero. The sandy loam soil showed higher values of D at saturation (D0) and total sorptivity (S0), as compared with the loam and silty loam soils. In conclusion, it was possible to determine D (θ) using a simple and robust setup, based on Whisler et al. (1968) method. Additionally, the proposed setup represents a simple way to measure soil S (θ) from horizontal infiltration, and to obtain K(h) functions through inverse parameterization.

The effect of microbial activity on soil water diffusivity.

In this study, we explored the effects of microbial activity on the evaporation of water from cores of a sandy soil under laboratory conditions. We applied treatments to stimulate microbial activity by adding different amounts of synthetic analogue root exudates. For comparison, we used soil samples without synthetic root exudates as control and samples treated with mercuric chloride to suppress microbial activity. Our results suggest that increasing microbial activity reduces the rate of evaporation from soil. Estimated diffusivities in soil with the largest amounts of added root exudates were one third of those estimated in samples where microbial activity was suppressed by adding mercuric chloride. We discuss the effect of our results with respect to water uptake by roots.Highlights We explored effects of microbial activity on the evaporation of water from cores of a sandy soil.We found the effect of microbial activity on water release characteristic was small.Increasing microbial activity reduced evaporation from soil, while microbial suppression increased it.Effect of microbial activity on root water uptake was estimated to be equivalent to a change in soil structure.