Pedotransfer functions estimating soil hydraulic properties using different soil parameters

Estimates of soil hydraulic properties using pedotransfer functions (PTF) are useful in many studies such as hydrochemical modelling and soil mapping. The objective of this study was to calibrate and test parametric PTFs that predict soil water retention and unsaturated hydraulic conductivity parameters. The PTFs are based on neural networks and the Bootstrap method using different sets of predictors and predict the van Genuchten/Mualem parameters. A Danish soil data set (152 horizons) dominated by sandy and sandy loamy soils was used in the development of PTFs to predict the Mualem hydraulic conductivity parameters. A larger data set (1618 horizons) with a broader textural range was used in the development of PTFs to predict the van Genuchten parameters. The PTFs using either three or seven textural classes combined with soil organic mater and bulk density gave the most reliable predictions of the hydraulic properties of the studied soils. We found that introducing measured water content as a predictor generally gave lower errors for water retention predictions and higher errors for conductivity predictions. The best of the developed PTFs for predicting hydraulic conductivity was tested against PTFs from the literature using a subdata set of the data used in the calibration. The test showed that the developed PTFs gave better predictions (lower errors) than the PTFs from the literature. This is not surprising since the developed PTFs are based mainly on hydraulic conductivity data near saturation and sandier soils than the PTFs from the literature. Copyright © 2007 John Wiley & Sons, Ltd.     

Climate,vegetation, and soil controls on hydraulic redistribution in shallow tree roots

Hydraulic redistribution defined as the translocation of soil moisture by plant root systems in response to water potential gradients is a phenomenon widely documented in different climate, vegetation, and soil conditions. Past research has largely focused on hydraulic redistribution in deep tree roots with access to groundwater and/or winter rainfall, while the case of relatively shallow (i.e., ≈1–2 m deep) tree roots has remained poorly investigated. In fact, it is not clear how hydraulic redistribution in shallow root zones is affected by climate, vegetation, and soil properties. In this study, we developed a model to investigate the climate, vegetation, and soil controls on the net direction and magnitude of hydraulic redistribution in shallow tree root systems at the growing season to yearly timescale. We used the model to evaluate the effect of hydraulic redistribution on the water stress of trees and grasses. We found that hydraulic lift increases with decreasing rainfall frequency, depth of the rooting zone, root density in the deep soil and tree leaf area index; at the same time for a given rainfall frequency, hydraulic lift increases with increasing average rainstorm depth and soil hydraulic conductivity. We propose that water drainage into deeper soil layers can lead to the emergence of vertical water potential gradients sufficient to explain the occurrence of hydraulic lift in shallow tree roots without invoking the presence of a shallow water table or winter precipitation. We also found that hydraulic descent reduces the water stress of trees and hydraulic lift reduces the water stress of grass with important implications on tree–grass interactions.     

Determining hydraulic properties of a loam soil by alternative infiltrometer techniques

Testing infiltrometer techniques to determine soil hydraulic properties is necessary for specific soils. For a loam soil, the water retention and hydraulic conductivity predicted by the BEST (Beerkan Estimation of Soil Transfer parameters) procedure of soil hydraulic characterization was compared with data collected by more standard laboratory and field techniques. Six infiltrometer techniques were also compared in terms of saturated soil hydraulic conductivity, K_s. BEST yielded water retention values statistically similar to those obtained in the laboratory and K_s values practically coinciding with those determined in the field with the pressure infiltrometer (PI). The unsaturated soil hydraulic conductivity measured with the tension infiltrometer (TI) was reproduced satisfactorily by BEST only close to saturation. BEST, the PI, one‐potential experiments with both the TI and the mini disk infiltrometer (MDI), the simplified falling head (SFH) technique and the bottomless bucket (BB) method yielded statistically similar estimates of K_s, differing at the most by a factor of three. Smaller values were obtained with longer and more soil‐disturbing infiltration runs. Any of the tested infiltration techniques appears usable to obtain the order of magnitude of K_s at the field site, but the BEST, BB and PI data appear more appropriate to characterize the soil at some stage during a rainfall event. Additional investigations on both similar and different soils would allow development of more general procedures to apply infiltrometer techniques for soil hydraulic characterization. Copyright © 2015 John Wiley & Sons, Ltd.     

A modified multiple tension upward infiltration method to estimate the soil hydraulic properties

Determination of saturated hydraulic conductivity, K_s, and the shape parameters α and n of the water retention curve, θ(h), is of paramount importance to characterize the water flow in the vadose zone. This work presents a modified upward infiltration method to estimate K_s, α and n from numerical inverse analysis of the measured cumulative upward infiltration (CUI) at multiple constant tension lower boundary conditions. Using the HYDRUS‐2D software, a theoretical analysis on a synthetic loam soil under different soil tensions (0, 0–10, 0–50 and 0–100 cm), with and without an overpressure step of 10 cm high from the top boundary condition at the end of the upward infiltration process, was performed to check the uniqueness and the accuracy of the solutions. Using a tension sorptivimeter device, the method was validated in a laboratory experiment on five different soils: a coarse and a fine sand, and a 1‐mm sieved loam, clay loam and silt‐gypseous soils. The estimated α and n parameters were compared to the corresponding values measured with the TDR‐pressure cell method. The theoretical analysis demonstrates that K_s and θ(h) can be simultaneously estimated from measured upward cumulative infiltration when high (>50 cm) soil tensions are initially applied at the lower boundary. Alternatively, satisfactory results can be also obtained when medium tensions (<50 cm) and the K_s calculated from the overpressure step at the end of the experiment are considered. A consistent relationship was found between the α (R² = 0.86, p < 0.02) and n (R² = 0.97, p < 0.001) values measured with the TDR‐pressure cell and the corresponding values estimated with the tension sorptivimeter. The error between the α (in logarithm scale) and n values estimated with the inverse analysis and the corresponding values measured with pressure chamber were 3.1 and 6.1%, respectively. Copyright © 2016 John Wiley & Sons, Ltd.     

Soil water balance scenario studies using predicted soil hydraulic parameters

Pedotransfer functions (PTFs) have become a topic drawing increasing interest within the field of soil and environmental research because they can provide important soil physical data at relatively low cost. Few studies, however, explore which contributions PTFs can make to land‐use planning, in terms of examining the expected outcome of certain changes in soil and water management practices. This paper describes three scenario studies that show some aspects of how PTFs may help improve decision making about land management practices. We use an exploratory research approach using simulation modelling to explore the potential effect of alternative solutions in land management. We: (i) evaluate benefits and risks when irrigating a field, and the impact of soil heterogeneity; (ii) examine which changes can be expected (in terms of soil water balance and supply) if organic matter content is changed as a result of an alternative management system; (iii) evaluate the risk of leaching to deeper horizons in some soils of Hungary. Using this research approach, quantitative answers are provided to ‘what if?’ type questions, allowing the distinction of trends and potential problems, which may contribute to the development of sustainable management systems. Copyright © 2005 John Wiley & Sons, Ltd.     

A comparative study of closed-form hydraulic conductivity prediction models

The Mualem and the Burdine hydraulic conductivity prediction models are considered in combination with the van Genuchten analytical retention curve, as well as the Brooks and Corey prediction model. An equivalence is presented between the retention curves of these models. A comparative study follows between hydraulic conductivities that are based on equivalent retention curves. A unified presentation of prediction models provides a framework for the whole analysis. The treatment of the equivalence problem consists in a minimization procedure characterized by uncoupling of the parameters and analytical evaluation of the objective function. Exact analytical equivalence relations are given for significant parts of the parameter ranges, and, for the remaining parts, analytical approximations are proposed. The comparisons between hydraulic conductivities are carried out via an inequality analysis. It is shown that the hydraulic conductivity of the Burdine model is less than that of the other models for extended ranges of equivalent parameters.     

Finite element computation of shape factors for Kirkham's piezometer tube method for hydraulic conductivity of saturated soil

The Kirkham piezometer tube method of measuring hydraulic conductivity of saturated soil uses a numerical constant (shape factor) to account for the relative geometrical parameters of the cavity and the soil domain. Published shape factors have mostly been determined using an electrical analogue. In this paper the application of the finite element method to shape factor determination is considered. Convergence is studied and methods of extrapolation evaluated. A determination is made of the size of finite region required to represent adequately the radially unbounded domain assumed for the piezometer method. Computed shape factors were found to be in good agreement with analogue data.     

A concise parameterization of the hydraulic conductivity of unsaturated soils

The parameter n in the well-known expression for hydraulic conductivity K=K₀S_eⁿ (where K₀ is its value at satiation and S_e the effective saturation) is determined as a function of the exponent in the power form of the soil–water retention relationship. The result is validated with an extensive experimental database comprising some 43 soils, collected by Mualem.     

Multi‐decadal water table manipulation alters peatland hydraulic structure and moisture retention

A peatland complex disturbed by berm construction in the 1950s was used to examine the long‐term impact of water table (WT) manipulation on peatland hydraulic properties and moisture retention at three adjacent sites with increasing depth to WT (WET, INTermediate reference and DRY). Saturated hydraulic conductivity (K_s) was found to decrease with depth by several orders of magnitude over a depth of 1–1.5 m at all sites. The depth dependence of WT response to rainfall was similar across sites: WT response increased from 1 : 1 at the surface, to 5 : 1 at 50 cm depth. While surface specific yield (S_y) values were similar across all sites, it decreased with depth at a rate of 0.014 cm^?1 in hollows and 0.007 cm^?1 in hummocks. Bulk density (ρ_b) exhibited similar depth‐dependent trends as S_y and explains a high amount of variance (r² > 0.69) in moisture retention across a range of pore water pressures (?15 to ?500 cm H₂O). Because of higher ρ_b, hollow peat had greater moisture retention, where site effects were minimal. However, the estimated residual water content for surface Sphagnum samples, while on average lower in hummocks (0.082 m³ m^?3) versus hollows (0.087 m³ m^?3), increased from WET (0.058 m³ m^?3) to INT (0.088 m³ m^?3) to DRY (0.108 m³ m^?3) which has important implications for moisture stress under conditions of persistent WT drawdown. Given the potential importance of microtopographic succession for altering peatland hydraulic structure, our findings point to the need for a better understanding of what controls the relative height and proportional coverage of hummocks in relation to long‐term disturbance‐response dynamics. Copyright © 2014 John Wiley & Sons, Ltd.     

Identification of optimal soil hydraulic functions and parameters for predicting soil moisture

Abstract

The accuracy of six combined methods formed by three commonly-used soil hydraulic functions and two methods to determine soil hydraulic parameters based on a soil hydraulic parameter look-up table and soil pedotransfer functions was examined for simulating soil moisture. A novel data analysis and modelling approach was used that eliminated the effects of evapotranspiration so that specific sources of error among the six combined methods could be identified and quantified. By comparing simulated and observed soil moisture at six sites of the USDA Soil Climate Analysis Network, we identified the optimal soil hydraulic functions and parameters for predicting soil moisture. Through sensitivity tests, we also showed that adjusting only the soil saturated hydraulic conductivity, K_s , is insufficient for representing important effects of macropores on soil hydraulic conductivity. Our analysis illustrates that, in general, soil hydraulic conductivity is less sensitive to K_s than to the soil pore-size distribution parameter.

Editor D. Koutsoyiannis; Associate editor D. Hughes

Citation Pan, F., McKane, R.B. and Stieglitz, M., 2012. Identification of optimal soil hydraulic functions and parameters for predicting soil moisture. Hydrological Sciences Journal, 57 (4), 723–737.     

On estimating the hydraulic properties of soil,Part 2. A new empirical equation for estimating hydraulic conductivity for sands

A new empirical equation to estimate hydraulic conductivity is proposed, based on a large set of measured data for hydraulic properties of soil. The equation is simpler and more accurate than the series-parallel model. Under conditions of insufficient data, the new equation provides a good estimation of hydraulic conductivity for sands. For the same class of soils, another empirical equation is proposed to estimate the power N in the Averjanov-Irmay function.     

Establishing the role of pedogenesis in changing soil hydraulic properties

Chronosequences provide suitable sources for the study of changes in soil hydraulic behaviour as a result of long-term pedogenesis. For a podsol chronosequence in the Scottish Highlands, data are presented to indicate the changes that have occurred over 13000 years in the saturated hydraulic conductivity (K_sat) in each horizon. As the soil profile has evolved into a differentiated sequence of three horizons, the resulting hydrological changes can be both measured and quantified by relating K_sat to textural properties and bulk density. The results are significant for interpretation of changing runoff processes and slope stability.     

Development of pedotransfer functions for soil hydraulic properties in the critical zone on the Loess Plateau,China

Soil hydraulic properties (SHPs) including the soil water retention curve and saturated soil hydraulic conductivity (Ks) are crucial input data for simulations of soil water and solute transport in the Earth's critical zone. However, obtaining direct measurements of SHPs at a wide range of scales is time consuming and expensive. Pedotransfer functions (PTFs) are employed as an alternative method for indirectly estimating these parameters based on readily measured soil properties. However, PTFs for SHPs for the deep soil layer in the Earth's critical zone are lacking. In this study, we developed new PTFs in the deep soil profile for Ks and soil water retention curve on the Loess Plateau, China, which were fitted with the van Genuchten equation. In total, 206 data sets comprising the hydraulic and basic soil properties were obtained from three typical sites. Samples were collected from the top of the soil profile to the bedrock by soil core drilling. PTFs were developed between the SHPs and basic soil properties using stepwise multiple linear regression. The PTFs obtained the best predictions for Ks (R_adj² = 0.561) and the worst for van Genuchten α (R_adj² = 0.474). The bulk density and sand content were important input variables for predicting Ks, α, and θs, and bulk density, clay content, and soil organic carbon were important for n. The PTFs developed in this study performed better than existing PTFs. This study contains the first set of PTFs of SHPs to be developed for the deep profile on the Loess Plateau, and they may be applicable to other regions.     

Validity limits for the van Genuchten–Mualem model and implications for parameter estimation and numerical simulation

The calculation of the relative hydraulic conductivity function based on water retention data is an attractive and widely used approach, since direct measurements of unsaturated conductivities are difficult. We show theoretically under which conditions an air-entry value for water retention data is definitely required when using the statistical approach of Mualem. Moreover we rigorously specify the conditions for which the classical van Genuchten–Mualem model leads to wrong predictions of relative hydraulic conductivity and, hence, an alternative formulation including an air-entry value should be used. Significant consequences are demonstrated for the inverse parameter estimation based on multistep outflow experiments. Furthermore it is shown that the use of a physically correct formulation of the water retention curve including an air-entry value and the derived hydraulic conductivity function influences not only the stability of numerical simulations but also their final results. This is especially grave as simulations with van Genuchten–Mualem parameters are frequently used to compare experiments and simulations and to draw conclusions on the correctness of Richards’ equation.     

Upscaling hydraulic conductivity based on the topology of the sub-scale structure

The hydraulic conductivity of heterogeneous porous media depends on the distribution function and the geometry of local conductivities at the smaller scale. There are various approaches to estimate the effective conductivity K_eff at the larger scale based on information about the small scale heterogeneity. A critical geometric property in this ‘upscaling’ procedure is the spatial connectivity of the small-scale conductivities. We present an approach based on the Euler-number to quantify the topological properties of heterogeneous conductivity fields, and we derive two key parameters which are used to estimate K_eff. The required coefficients for the upscaling formula are obtained by regression based on numerical simulations of various heterogeneous fields. They are found to be generally valid for various different isotropic structures. The effective unsaturated conductivity function K_eff (ψ_m) could be predicted satisfactorily. We compare our approach with an alternative based on percolation theory and critical path analysis which yield the same type of topological parameters. An advantage of using the Euler-number in comparison to percolation theory is the fact that it can be obtained from local measurements without the need to analyze the entire structure. We found that for the heterogeneous field used in this study both methods are equivalent.     

Use of practical aspects of soil behaviour to evaluate different methods to generate soil hydraulic functions

Soil hydraulic functions can be obtained with methods that range from complex and costly to simple and cheap. Decisions as to which is the most appropriate method for a specific application have to be based on a comparison of generated hydraulic functions. This comparison should preferably be based on a statistical comparison of practical applications calculated with the different hydraulic functions rather than on a statistical comparison of the functions themselves. In this study four different methods were used to generate hydraulic functions: (A) direct on-site measurement, (B) measurement in soil horizons in the area, (C) use of a national data set, and (D) use of Van Genuchten parameters correlated with soil texture and organic matter content. The four methods were compared by their effect on two practical aspects of soil behaviour: (1) evapotranspiration deficit and (2) flux through a plane at 30 cm below soil surface. These two aspects are highly relevant for agricultural and environmental use. However, direct measurement is not feasible. A validated simulation model was used for the calculations and results obtained with method A were taken as a reference. Calculations were performed for three soil profiles for a period of seven years. Deficits and fluxes, calculated with the four methods to generate hydraulic functions, were not significantly different using the data of the seven-year period. However, methods were significantly different when rainfall deficits were used as a covariable. This is true with the exception of downward fluxes in the period October until March which are most important for leaching of pollutants. The user has to decide whether differences between methods are sufficiently large to justify repeated, expensive on-site measurements (method A) or whether an investment will be made to make standard series of curves to be used everywhere (methods C and D).     

Delineation of hydrologically similar units in a watershed based on fuzzy classification of soil hydraulic properties

Evaluation of flow and transport processes in a watershed‐scale requires that the watershed be divided into homogenous spatial units referred to as hydrologically similar units (HSUs). Although a few discretization schemes are already in use, a universally acceptable method of obtaining HSUs is yet to emerge. In this study, we developed a fuzzy inference system (FIS) to classify the saturated hydraulic conductivity (K_s) and two water‐retention parameters α and n into fuzzy logic‐based soil hydrologic classes (FSHCs). Analysis of these classes showed that soil properties within an FSHC have less variability and those between two FSHCs have large variability. This result suggested that soils belonging to a specific FSHC may be more similar than those across different FSHCs and may be grouped together to represent an HSU. Soils within a specific hydrologic class were aggregated to delineate HSUs within the watershed. For the Dengei Pahad micro‐watershed (DPW), this approach showed five distinct regions representing a discretized zone having similar soil hydraulic properties. Application of this approach on a larger international database of soil hydraulic properties revealed that the developed hydrologic classes are quite comparable across different databases. The delineated HSUs based on these FSHCs were also better than the soil series map of the watershed in maintaining the soil heterogeneity of the watershed. Moreover, this new discretization scheme using the SWAT modelling environment showed better performance than the soil series‐based discretization approach. Copyright © 2010 John Wiley & Sons, Ltd.     

A dripper-TDR method for in situ determination of hydraulic conductivity and chemical transport properties of surface soils

Field determined hydraulic and chemical transport properties can be useful for the protection of groundwater resources from land-applied chemicals. Most field methods to determine flow and transport parameters are either time or energy consuming and/or they provide a single measurement for a given time period. In this study, we present a dripper-TDR field method that allows measurement of hydraulic conductivity and chemical transport parameters at multiple field locations within a short time period. Specifically, the dripper-TDR determines saturated hydraulic conductivity (K_s), macroscopic capillary length (λ_c), immobile water fraction (θ_im/θ), mass exchange coefficient (α) and dispersion coefficient (D_m). Multiple dripper lines were positioned over five crop rows in a field. Background and step solutions were applied through drippers to determine surface hydraulic conductivity parameters at 44 locations and surface transport properties at 38 locations. The hydraulic conductivity parameters (K_s, λ_c) were determined by application of three discharge rates from the drippers and measurements of the resultant steady-state flux densities at the soil surface beneath each dripper. Time domain reflectometry (TDR) was used to measure the bulk electrical conductivity of the soil during steady infiltration of a salt solution. Breakthrough curves (BTCs) for all sites were determined from the TDR measurements. The K_s and λ_c values were found to be lognormally distributed with average values of 31.4 cm h⁻¹ and 6.0 cm, respectively. BTC analysis produced chemical properties, θ_im/θ, α, and D_m with average values of 0.23, 0.0036 h⁻¹, and 1220 cm² h⁻¹, respectively. The estimated values of the flow and transport parameters were found to be within the ranges of values reported by previous studies conducted at nearby field locations. The dripper TDR method is a rapid and useful technique for in situ measurements of hydraulic conductivity and solute transport properties. The measurements reported in this study give clear evidence to the occurrence of non-equilibrium water and chemical movement in surface soil. The method allows for quantification of non-equilibrium model parameters and preferential flow. Quantifying the parameters is a necessary step toward determining the influences of surface properties on infiltration, runoff, and vadose zone transport.