Soil Hydraulic Conductivity Functions Research Articles

Slope stability analysis of cement and silica fume stabilized expansive soil slope

The objective of this study was to assess the impact of chemical stabilization on slope stability using water retention test data. The analysis focused on hypothetical expansive soil with a 10 m height and a 34° angle. Two critical unsaturated hydraulic properties were considered: (1) soil water retention curve and (2) hydraulic conductivity function. Chemical stabilizers were introduced at varying percentages relative to cement content, with silica fume (SF) partially replacing portions of the cement. The influence of stabilizers on slope stability was evaluated by examining changes in the factor of safety (FOS) with different stabilizer concentrations. Generally, the FOS increased with higher cement content, with further improvement observed when cement was partially replaced by silica fume. Additionally, the addition of silica fume enhanced water retention capacity, as evidenced by altered pore pressure distribution and water table levels post-rainfall infiltration.

Development and Application of an IoT-Based System for Soil Water Status Monitoring in a Soil Profile.

Soil water content (θ), matric potential (h) and hydraulic conductivity (K) are key parameters for hydrological and environmental processes. Several sensors have been developed for measuring soil θ-h-K relationships. The cost of such commercially available sensors may vary over several orders of magnitude. In recent years, some sensors have been designed in the framework of Internet of Things (i.e., IoT) systems to make remote real-time soil data acquisition more straightforward, enabling low-cost field-scale monitoring at high spatio-temporal scales. In this paper, we introduce a new multi-parameter sensor designed for the simultaneous estimation of θ and h at different soil depths and, due to the sensor's specific layout, the soil hydraulic conductivity function via the instantaneous profile method (IPM). Our findings indicate that a second-order polynomial function is the most suitable model (R2 = 0.99) for capturing the behavior of the capacitive-based sensor in estimating θ in the examined soil, which has a silty-loam texture. The effectiveness of low-cost capacitive sensors, coupled with the IPM method, was confirmed as a viable alternative to time domain reflectometry (TDR) probes. Notably, the layout of the sensor makes the IPM method less labor-intensive to implement. The proposed monitoring system consistently demonstrated robust performance throughout extended periods of data acquisition and is highly suitable for ongoing monitoring of soil water status.

Coupled hydro‐mechanical pore‐scale modeling of biopore‐coated clods for upscaling soil shrinkage and hydraulic properties

AbstractEarthworms and plant roots are vital for macropore formation and stabilization. The organo‐mineral coating of biopore surfaces also regulates macropore‐matrix mass exchange during preferential flow. The influence of finer‐textured burrow coatings on macroscopic soil properties during shrinkage could potentially be assessed by upscaling pore‐scale hydraulic and mechanical simulations. The aim was to investigate the influence of micro parameters (particle size, stiffness, and bond strength) on macro parameters (i.e., shrinkage curve and soil hydraulic properties). Drainage experiments and simulations were carried out using biopore‐coated clod‐size samples compared to those without coating. Simulations were performed using a two‐phase pore‐scale finite volume coupled with discrete element model (DEM‐2PFV). The structural dynamics was characterized by analyzing the pore volume and soil shrinkage curve obtained from numerically determined data. The soil hydraulic parameters were described using uni‐ and bimodal van Genuchten (vG) functions. The drainage simulations revealed hydro‐mechanical dynamics characterized by Braudeau‐shrinkage curve subdomains: The matrix‐only samples, with lower particle bond strength, exhibited relatively higher shrinkage. The coated samples, with higher particle stiffness and bond strength, displayed greater hydro‐mechanical stability. The numerically determined initial value of the saturated hydraulic conductivity (Ks) was about 70 times larger for matrix‐only samples than for coated samples. As expected, for the nonrigid soil structures, constant Ks, α, and n values for bimodal vG model resulted in prediction errors. Upscaling DEM‐2PFV pore‐scale model outcomes quantifies micro‐coating effects on macro hydro‐mechanics. This yields void ratio‐based soil water retention and hydraulic conductivity functions, advancing macroscopic soil hydraulic models and enhancing structured soil flow and transport descriptions.

A framework for estimating soil water characteristic curve and hydraulic conductivity function of permeable reactive media

The unsaturated behavior of permeable reactive barriers (PRB) is a critical component in predicting the removal efficiency through the adsorption of contaminants. This study investigates the framework to estimate the soil water characteristic curve (SWCC) and hydraulic conductivity function (HCF) for iron oxide-coated sand (IOCS) and zeolite, which are common materials used in PRBs. A multistep outflow (MSO) experiment was performed and the results of the MSO experiment were used to optimize associated parameters in Kosugi's SWCC and HCF. In addition, three scenarios of optimization analysis were investigated to evaluate the best-fitting model for estimating SWCC and HCF. The low root mean square error (RMSE) of fitted parameters indicates the Kosugi model well described the observed suction profiles in MSO experiments. In addition, the lowest RMSE and coefficient of variation suggested the inclusion of the additional parameter β provided the best estimation of the three materials (clean sand, IOCS, and zeolite). The physically reasonable estimation of SWCC and HCF of the three materials from the optimized parameters suggests the proposed framework is a reasonable model for the unsaturated behavior of PRBs.

Semi-analytical solution of pore-water pressure in unsaturated ground and infinite slope considering highly nonlinear soil hydraulic properties

A multi-exponential function (ME) is proposed to depict the highly nonlinear soil water retention curve and hydraulic conductivity function. Then, semi-analytical solution of pore-water pressure in an unsaturated infinite slope at steady state is derived using ME to depict soil hydraulic properties. The solution is verified with experimental and numerical results. It is found that ME can better describe the highly nonlinear hydraulic properties of soil than the single exponential function (SE), which is commonly-used previously. At steady state, the adoption of SE results in overestimating negative pore-water pressure (PWP; i.e., matric suction) by about 25 kPa and slope’s factor of safety (FOS) by 100% under rainfall, but underestimating negative PWP by about 50 kPa under evaporation. The difference between PWP and FOS calculated by using SE and ME becomes more significant for soil with bimodal hydraulic conductivity function than that with unimodal one. This difference reduces at a larger slope angle and deeper depth. Moreover, this difference generally increases at a lower SME/SSE, where SME and SSE represent the suction at which the unsaturated hydraulic conductivity expressed by ME and SE equals the applied rainfall intensity or evaporation rate, respectively. When SME/SSE is greater than 0.1, the difference basically diminishes.

Reducing uncertainties in hydromechanical modeling with a recently developed Rosetta 3 podeotransfer function

Stability analysis of unsaturated landslide deposits requires reliable estimates of soil moisture and pore water pressure. However, modeled soil moisture and pore water pressure contain substantial uncertainties due to imperfect information on soil hydraulic properties. Due to the relatively high dimensionality, commonly used parameter optimization strategies can be significantly affected by equifinality problems. This study investigates the effectiveness of reducing parameter estimation dimensionality using soil pedo-transfer functions. Specifically, we first estimated soil hydraulic parameters using the traditional Generalized Likelihood Uncertainty Estimation (GLUE) method, with parameters randomly drawn from the entire space (refer to as GLUE-random). In a second strategy, we use the Rosetta 3 pedotransfer function to constrain soil hydraulic parameters (refer to as GLUE-Rosetta). The two methods were tested in a typical landslide deposit with in-situ measured soil moisture dynamics for inverse modeling. The GLUE-random estimated soil hydraulic parameters contained substantial uncertainties –resulting in poorly constrained soil water retention curves (SWCC) and hydraulic conductivity functions (HCF). As a result, the uncertainty bands of pore water pressure and slope stability can cross values with several orders of magnitudes. In contrast, GLUE-Rosetta provided well-constrained SWCC and HCF, which significantly reduce the uncertainties in pore water pressure and slope stability estimates. These results suggest that the Rosetta 3 pedotransfer function can significantly improve the reliability of soil hydraulic parameters by reducing the dimensionality of the optimization problem and high-quality prior information of soil hydraulic properties. In conclusion, Rosetta 3 can enhance the reliability of soil parameters estimates and the reliability of subsurface hydrology, which may benefit the development of landslide early-warning systems.

Assessment of the interactions between soil–biosphere–atmosphere (ISBA) land surface model soil hydrology, using four closed-form soil water relationships and several lysimeters

Abstract. Soil water drainage is the main source of groundwater recharge and river flow. It is therefore a key process for water resource management. In this study, we evaluate the soil hydrology and the soil water drainage, simulated by the interactions between soil–biosphere–atmosphere (ISBA) land surface model currently used for hydrological applications from the watershed scale to the global scale, where parameters are generally not calibrated. This evaluation is done using seven lysimeters from two long-term model approach sites measuring hourly water dynamics between 2009 and 2019 in northeastern France. These 2 m depth lysimeters are filled with different soil types and are either maintained as bare soil or covered with vegetation. Four closed-form equations describing soil water retention and hydraulic conductivity functions are tested, namely the commonly used equations from Brooks and Corey (1966) and van Genuchten (1980), a combination of the van Genuchten (1980) soil water retention function with the Brooks and Corey (1966) unsaturated hydraulic conductivity function, and, for the very first time in a land surface model (LSM), a modified version of the van Genuchten (1980) equations, with a new hydraulic conductivity curve proposed by Iden et al. (2015). The results indicate good performance by ISBA with the different closure equations in terms of soil volumetric water content and water mass. The drained flow at the bottom of the lysimeter is well simulated, using Brooks and Corey (1966), while some weaknesses appear with van Genuchten (1980) due to the abrupt shape near the saturation of its hydraulic conductivity function. The mixed form or the new van Genuchten (1980) hydraulic conductivity function from Iden et al. (2015) allows the solving of this problem and even improves the simulation of the drainage dynamic, especially for intense drainage events. The study also highlights the importance of the vertical heterogeneity of the soil hydrodynamic parameters to correctly simulate the drainage dynamic, in addition to the primary influence of the parameters characterizing the shape of the soil water retention function.

Reduce uncertainty in soil hydrological modeling: A comparison of soil hydraulic parameters generated by random sampling and pedotransfer function

Numerical simulation of unsaturated soil hydrology relies on calibrated soil hydraulic parameters, which are subject to uncertainty due to imperfect information during the inverse modelling. This study investigates the effectiveness of reducing parameter uncertainty using the recently developed Rosetta 3 pedotransfer function. The GLUE method was employed for numerical modeling using the Darcy-Richards equation under two strategies for sampling Mualem-van Genuchten (MvG) parameters: the first uses conventional random generation of MvG parameters (GLUE-random), while the second adopts Rosetta 3 to transfer soil particle composition to MvG parameter (GLUE-Rosetta). Both approaches were used for inverse modeling of 9 typical soils, each with a recommended parameter set defined as true values and associated soil moisture dynamics as observations. The posterior parameters selected with both GLUE-random and GLUE-Rosetta show an equifinality phenomenon. GLUE-random fails to provide well-constrained posterior parameters to recover the pre-defined true values, and its posterior results of soil water characteristic curve (SWCC) and soil hydraulic conductivity function (HCF) are poorly constrained. In contrast, GLUE-Rosetta significantly improves the accuracy of the inversely-estimated soil hydraulic parameters, and the ensemble of posterior SWCC and HCF also encompasses the predefined true curves. The results demonstrate the effectiveness of using Rosetta 3 to reduce the dimensionality of the optimization problem, which results in reliable estimation of soil hydraulic parameters and soil particle compositions. Moreover, GLUE-Rosetta outperforms GLUE-random in predicting soil moisture dynamics under different rainfall intensities. Overall, it is recommended to integrate Rosetta 3 with existing optimization tools to reduce the uncertainty of soil parameters and support more reliable modeling of unsaturated soil hydrology.

Mixed formulation for an easy and robust numerical computation of sorptivity

Abstract. Sorptivity is one of the most important parameters for the quantification of water infiltration into soils. Parlange (1975) proposed a specific formulation to derive sorptivity as a function of the soil water retention and hydraulic conductivity functions, as well as initial and final soil water contents. However, this formulation requires the integration of a function involving hydraulic diffusivity, which may be undefined or present numerical difficulties that cause numerical misestimations. In this study, we propose a mixed formulation that scales sorptivity and splits the integrals into two parts: the first term involves the scaled degree of saturation, while the second involves the scaled water pressure head. The new mixed formulation is shown to be robust and well-suited to any type of hydraulic function – even with infinite hydraulic diffusivity or positive air-entry water pressure heads – and any boundary condition, including infinite initial water pressure head, h→-∞. Lastly, we show the benefits of using the proposed formulation for modeling water into soil with analytical models that use sorptivity.

Assessment of the Effects of the 2021 Caldor Megafire on Soil Physical Properties, Eastern Sierra Nevada, USA

The Caldor fire burned ~222,000 acres of the Eastern Sierra Nevada during summer–fall 2021. We evaluated the effects of this “megafire” on the physical properties of a sandy soil developed from glacial tills to document fire-induced soil modifications in this region. We measured soil water retention and hydraulic conductivity functions as well as the thermal properties of five core samples from control (unburned) areas and eight core samples from burned soil of the same soil unit. Soil water repellency was measured in terms of water drop penetration time (WDPT) in the field and apparent contact angle in the laboratory on control and burned soil as well as ash samples. Soil organic matter (SOM) and particle and aggregate size distributions were determined on control and burned soil samples. Additionally, scanning electron microscopy (SEM) was used to image microaggregates of control and burned soil samples. We found a significant difference in SOM content and sand and silt aggregate size distribution between control and burned samples, which we associated with the disintegration of microaggregates due to the fire. We found no significant difference between soil water retention and hydraulic conductivity functions of control and burned soil but observed greater variation in saturated hydraulic conductivity and systematic shifts in thermal conductivity functions of burned compared to control samples. WDPT and apparent contact angle values were significantly higher for burned soils, indicating the occurrence of fire-induced soil hydrophobicity (FISH). Interestingly, the average apparent contact angle of the control soil was >90°, indicating that even the unburned soil was hydrophobic. However, the ash on top of the burned soil was found to be hydrophilic, having apparent contact angles <10°. Our results indicate that SOM and microaggregates were readily affected by the Caldor fire, even for sandy soil with a weakly developed structure. The fire seemed to have moderated thermal properties, significantly and soil wettability but had only minimal effects on water retention and hydraulic conductivity functions. Our findings demonstrate the complex nature of fire-soil interactions in a natural environment and highlight the need for additional investigation into the causes and processes associated with FISH and structure alterations due to fire to improve our ability to rapidly determine potential problem areas in terms of hazards commonly associated with fire-affected soils.

Temperature effects on the hydraulic properties of unsaturated rooted soils

Soil temperature of vegetated landfill earthen covers varies with time during their service, due to heat released by the degradation of municipal solid waste. In this study, an experimental campaign to capture soil water-retention curve (SWRC) and hydraulic conductivity function (HCF) of unsaturated rooted soils was conducted at different temperatures (25, 45, and 60 °C). Simplified evaporation tests were conducted by subjecting soil specimens to a prolonged drying period at maintained temperature. Based on the responses of matric suction and water content, SWRC and HCF were determined by the instantaneous profile method. It was found that increasing temperature reduced water-retention capacity, due to root shrinkage, thermal expansion of liquid phase, and reduction in forces of capillarity and adsorption. Soil hydraulic conductivity was enhanced by the increase in soil temperature due to reduction of water viscosity and potential preferential flow at root channels induced by root shrinkage. Temperature effects on SWRC and HCF of rooted soils are suction-dependent. A new nonisothermal SWRC model was developed to capture the dual effects of roots and temperature. By using the same set of calibrated parameters, the model demonstrates its ability to predict SWRCs reasonably well within the range of temperature investigated.

Some Approaches to the Prediction of Permeability Parameters in a Finite Element Program for Early Warning

Recently, landslides often occurred in natural soil slopes in the tropical region, which correlate with the rainy season. Rainfall infiltration leads to groundwater level fluctuations. The increased positive pore-water pressures due to rainfall influence have affected the properties and behavior of the unsaturated soil slope. In this research, the Finite Element Method of SEEP/W and SLOPE/W analyzes the factor safety of the slope affected by pore water pressure change due to rainfall. The Soil Water Characteristic Curve (SWCC) and Hydraulic Conductivity function were obtained from sieve analysis and Atterberg's limit. In addition, unsaturated soil properties from the UNSODA code are estimated based on grain-size distribution using the SWRC program. The study area is in Khanom District, southern Thailand. The results show that the soil slope at the site became unstable on November 18, 2021, with F.S. = 1.0, which agrees well with the date of the disaster. In conclusion, the slope stability analysis without the parameters from the unsaturated soil hydraulic database (UNSODA) leads to the F.S. value being higher than the actual value, and the alarm estimation would be inaccurate. Doi: 10.28991/CEJ-2022-08-12-014 Full Text: PDF

Inverse Estimation of Soil Hydraulic Parameters in a Landslide Deposit Based on a DE-MC Approach

Extreme rainfall is a common triggering factor of landslide disasters, for infiltration and pore water pressure propagation can reduce suction stress and shear strength at the slip surface. The subsurface hydrological model is an essential component in the early-warning system of rainfall-triggered landslides, whereas soil moisture and pore water pressure simulated by the Darcy–Richards equation could be significantly affected by uncertainties in soil hydraulic parameters. This study conducted an inverse analysis of in situ measured soil moisture in an earthquake-induced landslide deposit, and the soil hydraulic parameters were optimized with the Differential Evolution Markov chain Monte Carlo method (DE-MC). The DE-MC approach was initially validated with a synthetic numerical experiment to demonstrate its effectiveness in finding the true soil hydraulic parameters. Besides, the soil water characteristic curve (SWCC) and hydraulic conductivity function (HCF) described with optimized soil hydraulic parameter sets had similar shapes despite the fact that soil hydraulic parameters may be different. Such equifinality phenomenon in inversely estimated soil hydraulic parameters, however, did not affect the performance of simulated soil moisture dynamics in the synthetic numerical experiment. The application of DE-MC to a real case study of a landslide deposit also indicated satisfying model performance in terms of accurate match between the in situ measured soil moisture content and ensemble of simulations. In conclusion, based on the satisfying performance of simulated soil moisture and the posterior probability density function (PDF) of parameter sets, the DE-MC approach can significantly reduce uncertainties in specified prior soil hydraulic parameters. This study suggested the integration of the DE-MC approach with the Darcy–Richards equation for an accurate quantification of unsaturated soil hydrology, which can be an essential modeling strategy to support the early-warning of rainfall-triggered landslides.

Determining Soil Available Water Capacity and Reasonable Irrigation Volume in Dryland Pasture

Among the methods used to describe soil-available water in terms of yield, field water capacity and permanent wilting coefficient were considered as the upper and lower boundaries, respectively, for normal plant growth. This study investigated the soil properties and groundwater resources of the Hailiutu pasture in Inner Mongolia, China. Soil water characteristics and hydraulic conductivity function curves of the pasture soil were measured based on filter paper method. Experimental results showed that air-entry and residual value were determined to be 34.6% and 13.1%, corresponding to the field water capacity and permanent wilting coefficient, respectively, and the available water capacity was 218.0 mm. The downward velocity of the soil wetting front in the silt is 3.0 mm/h. Site investigation revealed that the pasture was a natural capillary barrier, and its critical suction could optimize the upper boundary of irrigation. Considering the limited underground water resources, shallow water level, and the transpiration of Caragana Korshinskii in the pasture, an available drip irrigation method was designed. The spacing between drip irrigation pipes is 60 cm, the spacing between drip holes is 60 cm, the drip irrigation volume per hole is 475.2~838.8 mL/d, the drip irrigation time is 110~200 min at each night. Each hectare pasture exhausts 13.5~22.5 m3 water each day. Field tests showed that the proposed method is feasible and can be used for irrigation in dryland pastures.

Verification on the mode of moisture transfer in the vadose zone of thick loess

AbstractPrecipitation is the only source of groundwater recharge in the Chinese loess region. However, rainfall has been the main driver of geological disasters and engineering challenges in China, such as landslides, slope collapse and sinking of foundations. The roles of preferential flow and piston flow in moisture transfer in the thick loess vadose zone remain controversial. This study used a 100‐m thick loess shaft to monitor soil moisture content and soil temperature over 5 years, whereas climate indices were monitored on the Earth surface. The soil water characteristic curves and hydraulic conductivity function of the loess and the palaeosol collected in the profile were measured in a laboratory. The soil profile moisture transfer process was simulated based on the climatic boundary conditions and measured hydraulic parameters. The monitored and simulated results suggest that loess moisture transfer is dominated by piston flow. The thick loess vadose zone above the groundwater level could be divided into two sub‐zones, namely the active and steady zones. The active zone has a thickness of ~7 m and shows a direct response to climate change, with soil moisture content regulated by rainfall infiltration, evaporation and transpiration and soil moisture moving as transient flow. Although preferential flow can move in the active zone along joints, fissures and holes when water accumulates in the lower ground during heavy rainfall, piston flow remains the main mode of water movement. The steady zone has a thickness of <7 m and is above the groundwater level. Soil moisture content and soil temperature in the steady zone show temporal stability and soil moisture moves as a steady flow. Piston flow is the only mode of water movement in the steady zone. Since percolation velocity in the steady zone is estimated to be 73.3 mm/year, 1330 years of rainfall infiltration through the 97.5 m thick vadose zone would be required to recharge groundwater.

Estimating Soil Hydraulic Parameters of Different Textured Soils in Semiarid Duhok Conditions - Iraqi-Kurdistan Region by using RETC Program

RETC is a computer program that can be used to estimate soil hydraulic functions by easily determining soil properties. Using van Genuchten model in the RETC program gave high accuracy for estimating the hydraulic functions e.g., θ (h), k (θ) and D (h), also the mentioned program is considered a good tool for potential optimization of water management especially though calculating irrigation water after consuming 60% of the available water in root zone restoring technique.
 In the current study surface, soil samples (0 – 30 cm) from eight locations with a wide range in texture were collected RETC code is used to refine the water retention curve through minimizing the residual error according to van Genuchten – Mualem approach (1980). the obtained results indicated that the sum square SSQ (residual error) values becomes smaller after refining for eight studied soils which were (Avkendal, Berderash, Chemrash, Kanisark, Keperto, Khanik, Qerwola and Terbsbia), which their SSQ values before RETC refining were ( 0.0470264, 0.0680298, 0.0240080, 0.0190478, 0.0396200, 0.0166960, 0.094400 and 0.0199100) whereas these values decrease as fallowing (0.000333400, 0.000942618, 0.0000997008, 0.000167315, 0.002448800, 0.000243300, 0.00080000, and 0.000099917 ) with ha high coefficient of determination (R2) greater than (0.964) and the maximum refining was (0.000099917) at Chemrash silty clay loam soil with 0.9986 (R2), whereas the minimum value was(0.002448800) at Keperto loamy soil with 0.9720(R2) and this means higher precision in fitting using the retention model.
 RETC code was used to analyze soil water retention curve and hydraulic conductivity functions which are the key parameters in any quantitative description of water flow into the unsaturated zone of soils, van Genuchten et. al. (1991)

An improved method for measuring wetting front advancing velocity

AbstractThe unsaturated soil hydraulic conductivity function is one of the keys to solving unsaturated seepage issues. The wetting front advancing method (WFAM) is a proposed method for quickly measuring the coefficient of permeability of the unsaturated soils, but it has some demerits and limitations. In this study, a series of capillary rise tests were conducted and the coefficient of permeability of six kinds of unsaturated soils were obtained by using the WFAM. Based on the analyzing of characteristic changes in volumetric water content via soil moisture sensors in the capillary rise test, improvements were made to the method of determining the wetting front advancing velocity. The RMSE evaluation index was used to analyze the deviation of the data obtained by the improved method from the original visual method. The RMSE of the wetting front advancing velocity ranged from 0.0007 to 0.0074 cm s–1, and the RMSE of the coefficient of permeability ranged from 0 to 5.2 × 10–6 m s–1. The results indicated the applicability of the improved method.

Evaluation of pedotransfer functions for predicting soil hydraulic properties: A voyage from regional to field scales across Europe

Study regionEurope. A total of 660, 522, and 4940 soil samples belonging to GRIZZLY, HYPRES, and EU-HYDI databases, respectively, were used for parametric evaluation. Study focusThe soil water retention and hydraulic conductivity functions are crucial input information for land surface models. Determining these functions by using direct methods is hampered by excessive time and unaffordable costs required for field activities and laboratory analyses. Pedotransfer functions (PTFs) are widely-used indirect techniques enabling soil hydraulic properties to be predicted by using easily-retrievable soil information. In a parametric evaluation, the predictive capability of PTFs is examined by comparing measured and estimated soil water retention parameters and saturated hydraulic conductivity. Yet information about the performance of PTFs for specific modeling applications is mandatory to evaluate PTF effectiveness in greater depth. This approach is commonly defined as functional evaluation. New hydrological insights for the regionThe best performing four PTFs selected in the parametric evaluations are tested under two functional evaluations. The first encompasses a spatial interpolation with a geostatistical technique, whereas the second employs Hydrus-1D to simulate the water balance components along an experimental transect. Our results reinforce and integrate the insights of previous studies about the use of a PTF, and highlight the ability, or inability, of this technique to adequately reproduce the observed spatial variability of soil hydraulic properties and simulated water fluxes.

Hydraulic properties of unsaturated calcarenites by means of a new integrated approach

Most infiltration and water retention models for unsaturated porous media proposed in the literature often are based on simplified representations of the pore system and use the matric head as the main state variable.An integrated study on the hydraulic properties of homogeneous calcarenites with a complex pore network is presented, aimed at assessing water infiltration and retention mechanisms, starting from dry to wet conditions under high water inflow.Medium-grained grainstones and fine-grained packstones were used as material for laboratory testing and model set up. Petrophysical characterization and bimodal distribution of pore volume was obtained by standard and unconventional procedures. Permeability test under saturated conditions was performed with constant and falling head tests. Falling head infiltration tests at bench scale were carried out to estimate the hydraulic properties of the materials under unsaturated conditions.Unsaturated hydraulic relationships such as the soil water retention curve (SWRC) and hydraulic conductivity function were examined. The advection–diffusion form of Richard’s equation, together with the pore bundle model and dual porosity concept under equilibrium conditions represent the theoretical framework used to investigate the unsaturated hydraulic behavior of the calcarenites in the falling head infiltration test. The SWRC, relative hydraulic conductivity and celerity were determined through the capillary theory and pore bundle concept considering the bi-modal nature of the pore-size distribution. Theoretical solutions and model predictions were compared with experimental data, which demonstrated the dominant control of macro-porosity on water flow and infiltration rate. In particular, our experimental results indicate that the wetting front propagates mainly through the macropores with a gravity driven flow velocity higher than the diffusive flow through the micropores. The effect of texture, degree of saturation and boundary conditions on infiltration dynamics is discussed.

Rapid and accurate measurement methods for determining soil hydraulic properties: A review

Abstract The determination of soil hydraulic properties is important in several environmental sciences but may be expensive and time consuming. Therefore, during the last decades, a great effort has been made in soil sciences to develop relatively easy, robust, and inexpensive methods for soil hydraulic characterization. In this manuscript, we reviewed and discussed different infiltrometer techniques in light of the available experimental applications. More specifically, we considered the simplified falling head (SFH) infiltrometer technique and the single-ring infiltration experiment of the Beerkan type. Concerning this latter method, we considered different algorithms for data analysis: two simplified methods based on the analysis of transient (TSBI) and steady (SSBI) Beerkan infiltration data, and the Beerkan Estimation of Soil pedoTransfer parameters algorithm (BEST), that allows to estimate the soil characteristics curves, i.e., the soil water retention curve and hydraulic conductivity functions. For a given method, after dealing briefly theory and practice, available literature references were reported to account for specific applications in order to provide findings on method validation and application. With the aim to provide practical information on available tools for a simpler application of the reviewed methods, several video tutorials were reported to show i) how to conduct correctly field experiments and ii) how to calculate saturated hydraulic conductivity or soil hydraulic functions using user-friendly tools for data analysis. Finally, details on a new automated single-ring infiltrometer for Beerkan infiltration experiments (i.e., construction, assembly and field use) were presented.