Unsaturated Hydraulic Conductivity Function Research Articles

In-situ Measurement of Unsaturated Hydraulic Conductivity Function by a New Point Source Field Dripper Method (NPSFDM)

In-situ Measurement of Unsaturated Hydraulic Conductivity Function by a New Point Source Field Dripper Method (NPSFDM)

Deriving an unsaturated hydraulic conductivity function based on ergun equation

Computational Fluid Dynamics (CFD) is a widely-used numerical technique for simulating flow phenomena, particularly unsaturated flow, which involves two-phase flow. This type of flow can be effectively modeled using the Volume of Fluid (VOF) method. The key to accurately depicting unsaturated flow is the unsaturated hydraulic function (K-function), which connects flow velocity to the hydraulic gradient. However, employing the VOF method alongside the commonly used Ergun equation for modeling porous media can result in an unrealistic K-function. This, in turn, often leads to unrepresentative models of unsaturated flow. This study focuses on determining the equivalent particle diameter, an essential factor in the Ergun equation, and uses this to derive an improved K-function. The modified K-function exhibits promising results, closely matching experimental data and outperforming other empirical methods. This improved function could enhance the description and simulation of unsaturated flows.

SWRC interpretation of iron ore tailings blends focusing on saturation assessment of dry stacking

ABSTRACT The study addresses challenges in dry stacking of fine-grained materials by blending itabirite iron ore slime tailings with coarser flotation tailings. It focuses on flow dynamics, capillary saturation, and the necessity to prevent large, saturated volumes. Geotechnical tests were conducted on various blends, deriving unsaturated hydraulic conductivity functions using Fredlund and Xing’s methodology. The Soil-Water Retention Curve data, combined with transient 1D analysis, assesses capillary rise. Results show capillary rise up to 50 meters with over 80% saturation, and also influenced by weather conditions penetrating up to 20 meters. Integrating up to 30% slime content in flotation tailings mitigates saturation issues. Empirical and numerical findings enable inferring saturation tendencies in distinct iron ore tailings mixtures, establishing a procedural framework for practical field applications.

Prediction of hydraulic and petrophysical parameters from indirect measurements of electrical resistivity to determine soil-water retention curve – studies in granular soils

The characterization of unsaturated soils using hydromechanical methods is an essential requirement in soil science. However, current laboratory techniques used to obtain soil water retention and unsaturated hydraulic conductivity curves are time-consuming. To address this issue, a method based on indirect measures (electrical resistivity/electrical conductivity) was developed to quantitatively characterize soils. A novel unsaturated semi-empirical hydrogeophysical model of soils was developed by incorporating the hydrodynamic, geophysical, and petrophysical characteristics of soils. The model assumes that the parameters influencing the variation in the volumetric water content with matric suction and electrical resistivity are the same. The electrical resistivity characteristic curve (ERCC) defines a function that correlates environmental variables, electrical resistivity, soil water status, matric suction, hydraulic and petrophysical parameters, and fluid electrical resistivity. Model validation confirmed that the proposed approach can estimate the soil water retention curve (SWRC) via the indirect measures, and the results agreed with the experimental data. This indicates that it is possible to determine the SWRC and unsaturated hydraulic conductivity function of soil using the described approach.

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.

Evaluating mechanism and inconsistencies in hydraulic conductivity of unsaturated soil using newly proposed biochar conductivity factor

In the past few decades, numerous studies have been conducted to promote the use of biochar as a soil amendment and most recently, for compacted geo-engineered soils. In general, the definite trends of biochar effects on water retention and fertility of soils have been confirmed. However, the biochar effects on hydraulic conductivity, particularly unsaturated hydraulic conductivity of soil-biochar mix remain unclear, making it difficult to understand water seepage in both agricultural and geo-engineered infrastructures in semi-arid regions. This study examines the unsaturated hydraulic conductivity function derived based on the measurements of soil water characteristic curves of soil with biochar contents of 0%, 5% and 10%. A new parameter “biochar conductivity factor (BCF)” is proposed to evaluate the inconsistency in reported biochar effects on soil hydraulic conductivity and to interpret it from various mechanisms (inter- and intra- pore space filling, cracking, aggregation, bio-film formation and piping/internal erosion). The impact of biochar content on unsaturated hydraulic conductivity appears to reduce as the soil becomes drier with minimal effect in residual zone. Qualitative comparison of near-saturated hydraulic conductivity with test results in the literature showed that the BCF is generally higher for smaller ratio of sand to fine content (clay and silt). Moreover, the particle size of biochar may have significant influence on soil permeability. Future scope of research has been highlighted with respect to biochar production for its applications in agriculture and geo-environmental engineering. Long term effects such as root decay and growth, aggregation and nutrient supply need to be considered.Graphical

Examination of Measured to Predicted Hydraulic Properties for Low Impact Development Substrates

To counter the impacts of climate change and urbanization, engineers have developed ingenious solutions to reduce flooding and capture stormwater contaminants through the use of Low Impact Developments (LIDs). The soil is generally considered to be completely saturated when designing for the LIDs. However, this may not always be an accurate or realistic approach, as the soil could be variably unsaturated leading to inaccurate designs. To analyse the flow under variably unsaturated conditions, Richards’ equation can be used. To solve the Richards’ equation, two nonlinear hydraulic properties, namely soil water characteristic curve (SWCC) and the unsaturated hydraulic conductivity function are required. Laboratory and field measurements of unsaturated hydraulic properties are cumbersome, expensive and time- consuming. Pedotransfer functions (PTFs) estimate soil hydraulic properties using routinely measured soil properties. This paper presents a comparison between the direct measurement obtained through experimental procedures and the use of PTFs to estimate soil hydraulic properties for two green roof and three bioretention soil medias. Comparison between the measured and estimated soil hydraulic properties was accomplished using two different approaches. Statistical analyses and visual comparisons were used to compare the measured and estimated soil hydraulic properties. Additionally, numerical modelling to predict the water balance at the ground surface was conducted using the measured and estimated soil hydraulic properties. In some instances, the use of predicted hydraulic properties resulted in overestimation of the cumulative net infiltration of as much as 60 % for the green roof substrate, but was considered negligible for the bioretention substrate. Design performance criteria for green roof and bioretention facilities were examined using the measured and estimated soil hydraulic properties under extreme precipitation analysis. Results indicate that there is a high level of uncertainty when using PTFs for LID materials. A percent difference between the measured and predicted properties for the green roof peak time delay under a 2-year storm can be as much as 300%. For the bioretention design criteria of a 25-year storm, the surface runoff was overestimated by 14.7 cm and by 100% for the ponding time percent difference.

Prediction of the absolute hydraulic conductivity function from soil water retention data

Abstract. For modeling flow and transport processes in the soil–plant–atmosphere system, knowledge of the unsaturated hydraulic properties in functional form is mandatory. While much data are available for the water retention function, the hydraulic conductivity function often needs to be predicted. The classical approach is to predict the relative conductivity from the retention function and scale it with the measured saturated conductivity, Ks. In this paper we highlight the shortcomings of this approach, namely, that measured Ks values are often highly uncertain and biased, resulting in poor predictions of the unsaturated conductivity function. We propose to reformulate the unsaturated hydraulic conductivity function by replacing the soil-specific Ks as a scaling factor with a generally applicable effective saturated tortuosity parameter τs and predicting total conductivity using only the water retention curve. Using four different unimodal expressions for the water retention curve, a soil-independent general value for τs was derived by fitting the new formulation to 12 data sets containing the relevant information. τs was found to be approximately 0.1. Testing of the new prediction scheme with independent data showed a mean error between the fully predicted conductivity functions and measured data of less than half an order of magnitude. The new scheme can be used when insufficient or no conductivity data are available. The model also helps to predict the saturated conductivity of the soil matrix alone and thus to distinguish between the macropore conductivity and the soil matrix conductivity.

Pore network modeling approach for simulating soil water retention curve under different stress conditions

The field of unsaturated soil mechanics has recently seen the introduction of pore network models, which attempt to replicate the void structure of porous materials. They are robust physically-based simulation tools and have been used to simulate constitutive relationships like soil water retention curves (SWRC) and unsaturated hydraulic conductivity functions. This work aims to present a pore networkmodeling approach for predicting hysteretic SWRC at various stress states using only grain size distribution and porosity data of the granular soils. The soil sample subjected to given stress conditions is simulated using the Discrete Element Method to obtain a stable packing of spherical particles representing the soil structure. From this packing, an algorithm based on the medial axis is availed to extract a network of pores and throats that describes the geometry and topology of the void structure of the granular soils. Various pore-scale mechanisms like the piston-like advance, corner flow, pore body filling, and snap-off are then used to model fluid displacements at the pore scale to simulate SWRC along the drying and wetting paths. The modeled SWRC under various stress conditions is compared with the measured curves obtained for granular soils from the literature, and the predictions are in good agreement with the experimental results.

Estimation of unsaturated hydraulic conductivity function using extrapolated saturated front widths of line source field dripper

Estimation of unsaturated hydraulic conductivity function using extrapolated saturated front widths of line source field dripper

Field-scale assessment of the unsaturated hydraulic properties of residual soils in southeastern Brazil

Abstract Field tests were carried out to estimate effective unsaturated soil hydraulic properties of layered residual soils in Rio de Janeiro, southeastern Brazil. Data of this type are important for understanding the initiation of rainstorm-induced soil landslides, which often occur in the state of Rio de Janeiro as well as other areas having similar geologic settings and climate conditions. Tests were carried out using a simplified field approach, referred to as the Monitored Infiltration Test, which requires only a tensiometer to measure pressure heads below the wetting front, triggered by flow from a Mariotte bottle which maintains a constant pressure at the top edge of the soil profile. The data can then be analyzed by numerical inversion using the HYDRUS-2D software package. The test is relatively fast since no steady-state flow conditions are needed, and versatile since the test can be carried out quickly on steep slopes with the help of a manual auger. Soil water retention and the unsaturated hydraulic conductivity functions were obtained for a range of young, mature and saprolitic residual soils. The effective hydraulic properties of the distinct residual soil layers can be quite large, reflecting a need to provide a careful analysis of field-scale hydraulic heterogeneity in geotechnical analyses.

Measurement of the hydraulic properties of chalk using centrifuge permeameter; the study of chalk hydraulic properties under accelerated gravitational force

In this study, a geotechnical centrifuge has been used to study the unsaturated hydraulic properties of rock samples from the chalk aquifer, SE England. This method allows rapid measurement of hydraulic properties in a controlled environment, in contrast to previous studies on the chalk unsaturated zone, which required either an extended period (years) of data monitoring in the field or extended experimental periods (weeks) in the laboratory. Three types of specially built sensors were used to monitor water flow through chalk samples: a water pressure transducer to measure matric potential, frequency domain reflectometry probes to measure the volumetric water content and pressure transducers to measure the volume of water passing through the sample. Chalk samples were tested during wetting and draining processes to understand any hysteresis occurring during periodic recharge of the aquifer. Before undertaking physical tests of chalk samples in the centrifuge, a theoretical model of chalk hydraulic behaviour under centrifugal force was developed. This model was used to define and justify the instrumentation plan of the physical model and predict the shape of the soil moisture characteristic (SMC) curve and unsaturated hydraulic conductivity (Ku) function. The results were then evaluated and compared with the experimental results. These results show that chalk samples can be successfully tested under centrifuge conditions and hydraulic properties can be measured, including soil moisture characteristic curves and hydraulic conductivity.Supplementary material:Calculation spread sheets are available athttps://doi.org/10.17033/DATA.00000283

Merenje hidrauličkih osobina podloga za uzgoj sa HYPROP sistemom

Knowledge of hydro-physical properties is an essential prerequisite for assessing the suitability and quality of growing media. The method used for sample preparation is important for the measurement results. Three different sample preparation methods were compared. The methods differed in terms of the way the 250°cm3 steel cylinder was filled and the height of preloading. Measurements on loosely filled cylinders were included. The comparison was carried out on 15 growing media using the HYPROP device. HYPROP enables a complex analysis of the hydro-physical properties with high accuracy and reproducibility. The water retention curve, the unsaturated hydraulic conductivity function, the dry bulk density, the shrinkage and the rewetting properties can be measured simultaneously. The air capacity and the amount of plant-available water in pots depend on the height of the pot. In the field, it is related to the field capacity. The quality assessment was carried out both for flowerpots of different height and for field conditions with free drainage. Loosely filled samples consolidated hydraulically shortly after the start of the measurement. These geometric changes can be taken into account with the HYPROP. The sample preparation method - preloading or loose filling - yielded significantly different results for the pore volume, dry bulk density, plant available water and air capacity. The total pore volume of the loosely filled cylinders varied between 86.8 and 95.2°% by vol. (preloaded 81.3 and 87.7°% by vol.). The most critical factor was the air capacity. Loosely filled substrate samples achieved the highest air capacities, but also did not reach the critical value of 10°% by volume in shallow flowerpots, e.g. in 10 cm pots with 5.8°% by volume. The sample preparation method, measurement and quality assessment of the hydro-physical properties of growing media should be adapted to the conditions of use - whether they are used in a field with free drainage or in pots or containers in greenhouses.

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.

Effect of Unimodal and Bimodal Soil Hydraulic Properties on Slope Stability Analysis

Rainfall infiltration is the primary triggering factor of slope instability. The process of rainfall infiltration leads to changes in the water content and internal stress of the slope soil, thereby affecting slope stability. The soil water retention curve (SWRC) was used to describe the relationship between soil water content, matric suction, and the water retention characteristics of the soil. This characteristic is essential for estimating the properties of unsaturated soils, such as unsaturated hydraulic conductivity function and shear strength. Thus, SWRC is regarded as important information for depicting the properties of unsaturated soil. The SWRC is primarily affected by the soil pore size distribution (PSD) and has unimodal and bimodal features. The bimodal SWRC is suitable for soils with structural or dual-porous media. This model can describe the structure of micropores and macropores in the soil and allow the hydraulic behavior at different pore scales to be understood. Therefore, this model is more consistent with the properties of onsite soil. Few studies have explored the differences in the impact of unimodal and bimodal models on unsaturated slopes. This study aims to consider unimodal and bimodal SWRC to evaluate the impact of unsaturated slope stability under actual rainfall conditions. A conceptual model of the slope was built based on field data to simulate changes in the hydraulic behavior of the slope. The results of seepage analysis show that the bimodal model has a better water retention capacity than the unimodal model, and therefore, its water storage performance is better. Under the same saturated hydraulic conductivity function, the wetting front of the bimodal model moves down faster. This results in changes in the pressure head, water content, and internal stress of the soil. The results show that the water content and suction stress changes of the bimodal model are higher than those of the unimodal model due to the difference in water retention capacity. Based on the stability of the slope, calculated using the seepage analysis, the results indicate that the potential failure depth of the bimodal model is deeper than that of the unimodal model.

Importance of volumetric shrinkage curve (VSC) for determination of soil–water retention curve (SWRC) for low plastic natural soils

The study of unsaturated behavior of soils pivots around the relationship between soil suction and the amount of water content present in the soil mass, which is termed as soil–water retention curve (SWRC). The drying SWRC is generally determined by neglecting the volume change behavior of the soil for low plastic soils. Such a procedure is bound by an inherent assumption that the soil exhibit zero or negligible volume change and can lead to approximations in important SWRC parameters, including saturated volumetric water content, residual water content, and fitting parameters, which are essential for the estimation of unsaturated hydraulic characteristics of the soil. This study deals with the determination of SWRC parameters for a range of soils starting from non-plastic to high plastic soils, with and without considering volume change characteristics. For this purpose, the SWRC of six different natural soils were experimentally measured using a miniature tensiometer (T5) and a dew point potentiometer (WP4C), while the VSC of the used soils were measured using the balloon method. Efforts were made to critically understand the influence of no volume change assumption on the SWRC parameters, unsaturated hydraulic conductivity function, and flow of water through 1D soil column. The experimental results indicated that the no volume change assumption underpredicted the water retention curves of low plastic soils, and this error in SWRC increases with the soil plasticity. The measured VSC of the used soils further revealed that there could be a range of saturated volumetric water content (VWC) depending on the initial state of the soil sample, which also influences the sensitivity of the SWRC parameters. Therefore, the present study proposed a new methodology to determine unambiguous SWRC (in terms of VWC-suction) for a particular soil type combining both the measured VSC data and SWRC data. Based on the obtained results, it was recommended that the details of VSC must be incorporated along with the measured SWRC data for soils having plasticity and clay content higher than 15 and 11%, respectively.

Estimation of unsaturated hydraulic conductivity function: implication of low to high suction measurements

Establishment of the relationship between soil suction and water content, commonly termed as soil–water characteristic curve (SWCC), is of prime importance in the field of unsaturated soil mechanics. There are several instruments available that can be used to measure the SWCC of soil, but every suction measuring device has its own limitations in terms of its suction measurement range. Therefore, the preciseness of the estimated unsaturated soil properties largely depends on the range of suction measurements and the type of instruments used. The primary objective of this study is to quantify the error that can occur during the estimation of unsaturated hydraulic conductivity function (UHCF) from SWCC for low plastic soils. Experiments were performed to investigate the influence of different suction measurement devices on the estimated UHCF for four different soils with varying clay content. A dew point potentiometer (WP4) and a miniature tensiometer (T5) have been used in this study for the suction measurement. The SWCC of the selected soils were predicted mathematically using a commonly used pedo-transfer function (PTF). The experimental results clearly indicated that the sole use of WP4 overestimated the SWCC parameters, as well as UHCF (overestimation in the conductivity value is in order of 104 times). Rather, a combination of T5 and WP4 data, within their accurate range, provides a more precise estimation of UHCF. Further, the accuracy of the PTF was found very effective for low plastic soils with a relatively low percentage of clay (% clay < 10), in the absence of any experimental data.

Unsaturated Soil Hydraulic Properties according to Double-Ring Infiltration of Saline Water

Inverse modeling is a relatively complex procedure that allows quick estimation of soil hydraulic properties, yielding parameters for both soil water retention and unsaturated hydraulic conductivity function from a single experiment. The purpose of this research was to evaluate saline water effect on unsaturated hydraulic properties of clay loamy nonsaline-sodic soil layers and estimate these properties inversely from infiltration data using a double-ring infiltrometer.Three levels of water salinity (EC) 1.1, 2.2 and 5.8 dS/m were used for each treatment (T1, T2, and T3, respectively). Soil samples at three depths (0–20, 20–40 and 40–60 cm) were taken to measure Mualem–van Genuchten (MVG) hydraulic parameters (θr, residual water content, θs, saturated water content, α and n, shape parameters, and Ks, lab saturated hydraulic conductivity). Soil hydraulic parameters were estimated with an inverse solution using HYDRUS-1D model and infiltration data. Besides, the most sensitive MVG parameters to infiltration data (n, θs, and Ks) were inversely estimated due to difficulty in estimating several parameters simultaneously. Results showed that salinity had no significant effect on soil water infiltration and soil hydraulic conductivity changes in the specified salinity range according to the variance analyses (ANOVA).The Nash–Sutcliffe coefficient of efficiency (Ce) between simulated (HYDRUS) and measured values of cumulative infiltration data was about 0.998, indicating an excellent match. Soil water retention, as well as unsaturated soil hydraulic conductivity, were optimized with reasonable accuracy when the soil profile was assumed to be homogenous. In addition, simulated soil-water contents were quite similar to the measured values in which the values of Nash–Sutcliffe efficiency (Ce) for T1, T2, and T3 were 0.83, 0.85, and 0.85, respectively. In conclusion, the soil hydraulic parameters can be inversely estimated using the measured infiltration data through the double ring method if the soil is homogeneous.

Simulation of vadose zone flow processes via inverse modeling of modified multistep outflow for fine‐grained soils

AbstractCharacterization and measurement of the hydraulic properties of unsaturated porous media is still a challenge in natural environments, although exact knowledge of the soil's hydraulic properties [unsaturated soil hydraulic conductivity (K), volumetric water content (θ), soil matric head (h)] is crucial for solving many soil, hydrological, and environmental issues. The main purpose of this study was to establish a modified multistep outflow technique that facilitates laboratory operations and reduces time and costs. We also focused on inverse parameter optimization of the Mualem–van Genuchten (MVG) model with laboratory‐provided data for three fine‐grained soils. This modified version was assessed by inverse modeling of the MVG model and the generalized reduced gradient algorithm. The results represent the best estimate of the adjusted parameters and the very low uncertainty associated with the fitted values. The soil hydraulic functions were simulated well up to pressures of 1,500 kPa and their uncertainties were extremely low. The R2 values were 88 to 97 and 81 to 96% for the soil water characteristics and the unsaturated hydraulic conductivity functions, respectively. The simultaneous optimization of the nonlinear hydraulic functions from a single transient flow experiment agreed well with their observed data for each soil examined. Besides, by eliminating the use of tensiometers and having no independent measurements of saturated soil hydraulic conductivity and saturated volumetric water content, the modified version was faster and easier than the conventional method. The proposed protocol appears to be suited for the vadose zone flow simulation of fine‐grained soils.

Applicability of wetting front advancing method in the sand to silty clay soils

The unsaturated hydraulic conductivity function is one of the keys to solve unsaturated seepage issues, and it is difficult to be obtained via the laboratory test considering the time involved. In this study, a series of laboratory tests, include capillary rise test and downward infiltration test, were conducted to obtain the unsaturated hydraulic conductivity function of six kinds of soils varies from sand to silty clay by using the wetting front advancing method. The similarities and differences between instantaneous profile method and wetting front advancing method were discussed. The interpretation of data, demand details and the aspect that notice in test process, and application of the capillary rise test and downward infiltration test were summarized and analyzed. The demerits and its countermeasures for the method were given.