Hydraulic Conductivity Curves Research Articles

The PDI model system for parameterizing soil hydraulic properties

AbstractThe Peters–Durner–Iden (PDI) model system for describing soil hydraulic properties (SHP) has been developed over a decade. Inspired by Rien van Genuchten's seminal work, the PDI system focuses on an efficient and simple parameterization of water retention curves and hydraulic conductivity curves (HCC) across the entire soil moisture spectrum. By combining capillary and non‐capillary components for water retention and conductivity, it aims to reconcile mathematical simplicity and insights on water adsorption and film flow in soils. Recent developments have reduced the number of free parameters of the conductivity model to zero, enhancing the model's applicability in cases of limited data availability. The first reduction was achieved by a prediction of absolute non‐capillary conductivity based on the consideration of film and corner flow on the pore scale, and the second by a prediction of absolute capillary conductivity by a capillary bundle model. This allows a complete characterization of SHP over the entire moisture range with only four retention curve parameters. The inclusion of a maximum pore size in the capillary conductivity model prevents an unrealistic drop of the HCC near saturation. This paper provides a comprehensive overview of the PDI model system, emphasizing its conceptual features and mathematical details. An Excel sheet and a Python code stored in a repository are provided for accessibility.

Implementation and efficient evaluation of backpropagation network training algorithms in parametric simulations of soil hydraulic conductivity curve

The estimation of the soil hydraulic conductivity curve (SHCC) through machine learning tools for determining pedotransfer functions (PTFs) is essential in various agriculture, hydrological, and environmental management-related studies. The accuracy and reliability of machine learning tools such as backpropagation neural network (BNN) for estimating SHCC mainly depends on the optimization of its architecture, specifically the selection of the network training algorithms. By comparing the performance of network training algorithms, various researchers related to agriculture science, hydrology and environmental can choose the appropriate algorithm that provides the best balance of accuracy, speed, computational efficiency, and robustness for hydrological modeling purposes. Therefore, the present study assessed and compared the performance of 4 training algorithms including, gradient descent with momentum (GDM), Quasi-Newton (QN), Bayesian regularization (BR) and Levenberg-Marquardt (LM) for the estimation of the parametric-based PTFs of the van Genuchten-Mualem hydraulic conductivity equation as novelty of this research, which has not been investigated, so far. In this study four GDM, BR, LM and QN-based PTFs were derived to estimate the SHCC from different combinations of readily available properties of 86 soil samples from UNSODA database. The integral root mean square error (IRMSE) and coefficient of determination (R2) and were used to assess the estimations. The performance of BR algorithm-based PTFs on testing data varied from 4.634 to 6.690 with average 5.476 cm day−1 in terms of IRMSE, and 0.958 to 0.966 with average 0.963 in terms of R2, which provided better simulation than other training algorithms. The outcome of the present study on testing data showed that the BR algorithm was stable and robust to heterogeneous and noisy data compared to other training algorithms. Finally, development of accurate and reliable PTFs using efficient algorithms can play a crucial role in soil resource assessments, water resource management, irrigation planning, and designing effective drainage systems, because accurate knowledge of SHCC is vital for optimizing irrigation and fertilization practices and minimizing wastage of water and nutrients.

The effect of amorphous silica on soil–plant–water relations in soils with contrasting textures

This study investigates how amorphous silica (ASi) influences soil–plant–water interactions in distinct soil textures. A sandy loam and silty clay soil were mixed with 0 and 2% ASi, and their impact on soil retention and soil hydraulic conductivity curves were determined. In parallel, tomato plants (Solanum lycopersicum L.) were grown in experimental pots under controlled conditions. When plants were established, the soil was saturated, and a controlled drying cycle ensued until plants reached their wilting points. Soil water content, soil water potential, plant transpiration rate, and leaf water potential were monitored during this process. Results indicate a positive impact of ASi on the sandy loam soil, enhancing soil water content at field capacity (FC, factor of 1.3 times) and at permanent wilting point (PWP, a factor of 3.5 times), while its effect in silty clay loam was negligible (< 1.05 times). In addition, the presence of ASi prevented a significant drop in soil hydraulic conductivity (Kh\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${K}_{h}$$\\end{document}) at dry conditions. The Kh\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${K}_{h}$$\\end{document} of ASi-treated sandy loam and silty clay at PWP were 4.3 times higher than their respective control. Transpiration rates in plants grown in ASi-treated sandy loam soil under soil drying conditions were higher than in the control, attributed to improved soil hydraulic conductivity. At the same time, no significant difference was observed in the transpiration of plants treated with ASi in silty clay soil. This suggests ASi boosts soil–plant–water relationships in coarse-textured soils by maintaining heightened hydraulic conductivity, with no significant effect on fine-textured soils.

Can hydraulic-energy-indices be effectively used to describe the saturated hydraulic conductivity?

The saturated hydraulic conductivity (Ks) and water retention curve (SWRC) parameters are important properties for simulating soil hydrological processes and characterizing soil conservation around the world. Therefore, Ks and SWRC are related with the soil physical quality (SPQ) and several SPQ indices can be derived from SWRC, such as the pore size distribution, relative field capacity, plant available water, drainable porosity, and soil hydraulic-energy indices (SHEI). It is well known that the soil structure can be assessed by using SHEI, but a possible physical relationship between Ks and SHEI was not examined yet. Therefore, the objective of this study was to investigate the behavior of Ks as function of SHEI for several soil textural classes. If this relationship be proved, then SHEI might be applied to improve the Ks prediction by PTF models. In this work, a data set of 395 measured SWRC's were fitted to the vG equation to obtain the SHEI to verify whether they are statistically correlated and physically dependent on Ks. The resulting parametric and non-parametric correlation results were split up according to six textural classes. The significant influence of Ks on at least one of the absolute SHEI (Aa or WRa) was verified on the numerical scale when all textures were grouped and on numerical and pF scales for clayey and silty textures. Ks showed significant impact on Aa and WRa indices in four textural classes. Furthermore, Ks had influence on the sum Aa + WRa denoted in pF scale for five of the six textural classes, with a significant linear correlation in the clayey texture when log (Aa + WRa) was applied. The significant and high correlation of Ks on the ratios WRa/AWC and Aa/ϕD was also observed in four of the six classes, and therefore the use of these indices is recommended for the development of PTFs for Ks prediction.

A scaling-based model to describe temporal variability in soil hydraulic properties

Accurate knowledge of temporal variability in soil hydraulic properties (SHPs) can improve the prediction capability of flow and transport models. Lab and field measurements for determining SHPs at multiple timescales are expensive and time-consuming. Further, the existing Fokker-Plank equation (FPE) based numerical and analytical models for describing temporal variation in SHPs require parameterization. In this study, a scaling model is proposed to describe temporal variation in SHPs of four different soils. As field studies related to temporal variation in SHPs are limited, first synthetic temporally varying SHPs, soil water characteristic curves and hydraulic conductivity curves are generated using the analytical solution of FPE. Functional normalization approach is then used to develop reference curves and to obtain scale factors at different time-steps. The time-varying scale factors and the initial SHPs are related by two-step regression equations. The reference curves and the regression equations together can be used to estimate SHPs at any time given the initial SHPs. The developed model is validated using a 5-fold cross validation method and with experimental data at two study sites. The average percentage error in predicted van Genuchten parameters (θs,n and Ksat) except for α are mostly less than 10 %. The performance of the proposed model is comparable to that of FPE analytical model for experimental data. The proposed method holds promise for describing temporal variation in SHPs using only the initial SHPs.

Hydraulic conductivity over a wide suction range of loess with different dry densities

Experimental research into the hydraulic conductivity curve (HCC) of unsaturated soil is limited due to the inherent challenge associated with labor, cost, and time. Typically, the HCC is estimated using the soil water characteristic curve (SWCC) based models and saturated hydraulic conductivity (SHC). However, the efficiency of the SWCC-based model is rarely assessed, and the influence of soil density and pore structure on HCC remains incomplete due to limited experimental data. To address this gap, this study employs an innovative filter-paper-based column method, which can measure the HCC over a wide suction range (e.g. 0−105 kPa), to capture the HCCs of both intact and compacted specimens with varying dry densities. The efficiency of two typical SWCC-based models is assessed using the measured data. Meanwhile, the mercury intrusion porosity (MIP) technique is employed to obtain the pore characteristic (i.e. pore size distribution (PSD)) and a method of predicting the HCC using the PSD data is proposed, emphasizing the dominant role of the pore structure in shaping the HCC. The results reveal that the dry density's influence on the HCC is primarily observed within the low suction range, corresponding to variations in the dominant and large pores. In the high suction range, the HCCs align along a linear trajectory when plotted in a log-log format. A notable finding is the overestimation of the HCC obtained from the SWCC-based models using the measured SHC. When the SHC is regarded as a fitting parameter, good agreement is achieved. The adjusted SHC value is typically 0-1 order of magnitude lower than the measured value, and this discrepancy diminishes as dry density increases. On the other hand, the proposed PSD-based model performs well with the measured SHC data. Caution is exercised when using the SHC to estimate the HCC for modeling water movement in partially saturated soil.

Averaging or adding domain conductivities to calculate the unsaturated soil hydraulic conductivity

AbstractRecent models of the unsaturated hydraulic conductivity curve (UHCC) are the sum of separate UHCCs for domains of capillary water, film water, and water vapor. This requires parallel, noninteracting domains. A theoretical framework for aggregating domain conductivities to a bulk soil UHCC is presented to identify and possibly relax implicit assumptions about domain configuration. The paper develops arithmetic, harmonic, and geometric averages of the liquid‐water conductivities that can be arithmetically averaged with the vapor conductivity. However, current models for capillary and film conductivities are intrinsic, that is, valid within their respective domain. The vapor conductivity is a bulk conductivity, that is, it gives the conductivity of the gaseous domain as it manifests itself in the soil. Conversion relationships use the domain volume fractions as approximations of the as‐yet unknown weighting factors to convert between intrinsic and bulk conductivities. This facilitates consistent averaging of domain conductivities. Even with consistent averaging, a truly physically accurate model of the UHCC based on domain conductivities is fundamentally elusive. Nevertheless, models based on the three averages and the unweighted sum of the domain conductivities produce good fits to data for two soils but diverge in the dry range. The fitted curves for the capillary and film water depend on the averaging (or adding) method. Hence, they are not strictly characteristic of their respective domains. The true intrinsic domain conductivity functions may be impossible to determine.

Vulnerability of Xylem Embolism in Maize Cultivars with Different Drought Tolerance under Water and Salt Stress

Water deficit and soil salinization are the primary abiotic stress factors hindering maize growth. To assess the effect of water and salt stress on xylem embolism in maize and investigate the relationship between drought resistance and xylem vulnerability, a greenhouse experiment was designed using two maize cultivars, Zhengdan 958 (drought-resistant) and Denghai 605 (drought-sensitive). Four treatments were included: control (CK), water deficit (WD), salt stress (SS), and combined water and salt stress (WS). Various hydraulic characteristic indicators, such as stem xylem water potential, leaf xylem water potential, the specific hydraulic conductivity (Ks) and percentage loss of conductivity (PLC), were analyzed. Specific hydraulic conductivity curves and vulnerability curves were constructed, and the hydraulic safety margin (HSM) of the xylem was determined based on stomatal conductance (Gs). The results indicated that the hydraulic conductivity and embolism resistance of maize xylem were not correlated. Compared to Denghai 605, Zhengdan 958 had lower maximum specific hydraulic conductivity Ksmax and P50 values (xylem water potential at 50% PLC) in all treatments, indicating lower water transport capacity but stronger resistance to embolism. Under single-cultivar conditions, salt stress had a greater inhibitory effect on Ksmax and HSM in maize xylem compared to water deficit; thus, more severe embolism was found under salt stress. Under different treatment conditions, Zhengdan 958 had a larger HSM than Denghai 605, showing a wider water transport safety range and overall superior water transport security. To summarize, water and salt stress inhibited the water transport efficiency of the xylem in maize stems, and stronger drought-resistant cultivars showed greater resistance to embolism and larger hydraulic safety margins.

Saturated hydraulic conductivity and water retention curve of variably weathered tuff breccia bedrock in a headwater catchment

Although recent studies have emphasized the importance of water flow through bedrock to hydrological processes in headwater catchments, information on the hydraulic properties of bedrock, required for physically characterizing water flow processes, is lacking. In this study, the saturated hydraulic conductivity, Ks, and water retention curve of tuff breccia bedrock core were determined in samples collected from boreholes in a headwater catchment in northern Japan. These core-scale properties, together with core shape and in situ Ks values, were used to define two groups of bedrock: CM- and CL-class (weakly to moderately weathered) and D-class (highly weathered). Little change in the volumetric water content θ of CM- and CL-class bedrock was observed in the dry range at a pressure head ψ < –200 cm, but slight changes occurred in the wet range (ψ > –200 cm), resulting in low core-scale Ks values of 10–8–10–7 cm s−1. Core shape and in situ Ks values much larger than the core-scale values suggested that in situ water flow is restricted to fractures, with little contribution of the matrix to total flow. In D-class bedrock, on the other hand, the change in θ with ψ was somewhat larger than in CM- and CL-class bedrock, especially in the wet range, resulting in moderate core-scale Ks values of 10–5 cm s−1. These results indicated that flow in D-class bedrock is better represented by matrix flow than by preferential flow through fractures. The hydraulic properties presented herein provide a basis for analyzing subsurface water flow, including through bedrock layers. Those analyses will provide a better understanding of water flow processes in headwater catchments.

A stress-controlled soil column test equipment for investigating water seepage-consolidation characteristics of the compacted soil

For deep-filled building site, it is not acceptable to neglect water seepage and land settlement caused by the coupling effect of wetting and loading. However, previous soil column tests employed in investigating the water seepage always failed to consider the vertical stress. In this article, a stress-controlled soil column test was conducted to investigate the effect of vertical stress on water seepage and deformation of the compacted soil. The soil column was equipped with a vertical loading device and a water-recharge device to control the boundary conditions and water sensors, tensiometers, dial indicators and a data collector to monitor the displacement and water movement. The soil column test was utilised to perform seepage tests on compacted loess with vertical stresses of 0 and 400 kPa. The time-history data of wetting front depth, volumetric water content (VWC), suction and vertical displacement have been monitored directly. The VWC and suction profiles, soil–water characteristic curves and hydraulic conductivity curves were then obtained with the monitored data. The test results demonstrated that the vertical stress affected the water seepage and consolidation of the compacted soil, but this effect gradually decreased as the depth of the soil column increased. Moreover, the reliability of the proposed test method was verified by the comparison of some available test results. The findings of this study contribute to a better understanding of water seepage and consolidation characteristics of the compacted soil subjected to vertical stress.

A Simple Model for Estimating the Hydraulic Conductivity of Unsaturated Soil

Describing the hydraulic conductivity of unsaturated soil is very important in predicting water transport. Most current models have complex forms and generally need to be calibrated by the measured unsaturated hydraulic conductivity curve. A simple model, by which it is possible to conveniently predict the unsaturated hydraulic conductivity, is proposed in this study. The soil–water characteristic curve and hydraulic conductivity curve are separated into three parts. The soil–water characteristic curve is represented by Fredlund and Xing’s equation. A simple model composed of three lines is proposed for estimating the hydraulic conductivity of unsaturated soil. The model parameters can be conveniently calibrated from the measured soil–water characteristic curve and saturated hydraulic conductivity. Finally, the proposed model is validated by the experimental data from different soils. The proposed model provides a simple approach to estimating the hydraulic conductivity of unsaturated soil, which is more convenient for practical application.

Prediction of absolute unsaturated hydraulic conductivity – comparison of four different capillary bundle models

Abstract. To model water, solute, and energy transport in porous media, it is essential to have accurate information about the soil hydraulic properties (SHPs), i.e., the water retention curve (WRC) and the soil hydraulic conductivity curve (HCC). It is important to have reliable data to parameterize these models, but equally critical is the selection of appropriate SHP models. While various expressions for the WRC are frequently compared, the capillary conductivity model proposed by Mualem (1976a) is widely used but rarely compared to alternatives. The objective of this study was to compare four different capillary bundle models in terms of their ability to accurately predict the HCC without scaling the conductivity function by a measured conductivity value. The four capillary bundle models include two simple models proposed by Burdine (1953) and Alexander and Skaggs (1986), which assume a bundle of parallel capillaries with tortuous flow paths, and two more sophisticated models based on statistical cut-and-random-rejoin approaches, namely those proposed by Childs and Collis-George (1950) and the aforementioned model of Mualem (1976a). To examine how the choice of the WRC parameterization affects the adequacy of different capillary bundle models, we utilized four different capillary saturation models in combination with each of the conductivity prediction models, resulting in 16 SHP model schemes. All schemes were calibrated using 12 carefully selected data sets that provided water retention and hydraulic conductivity data over a wide saturation range. Subsequently, the calibrated models were tested and rated by their ability to predict the hydraulic conductivity of 23 independent data sets of soils with varying textures. The statistical cut-and-random-rejoin models, particularly the Mualem (1976a) model, outperformed the simpler capillary bundle models in terms of predictive accuracy. This was independent of the specific WRC model used. Our findings suggest that the widespread use of the Mualem model is justified.

Effect of the Mid-Layer on the Diversion Length and Drainage Performance of a Three-Layer Cover with Capillary Barrier

The capillary barrier is a type of soil cover system commonly used in various geotechnical applications, such as limiting infiltration for slopes or landfills or providing cover for solid waste. It serves to prevent the movement of water through the soil layers by utilizing contrasting particle sizes. This paper focuses on investigating the effect of the granular layer on the performance of a three-layer cover with a capillary barrier, integrating the granular layer within clayey sand. The investigation involved one-dimensional infiltration tests utilizing four uniform granular soils with varying grain sizes. These tests were instrumental in calibrating soil water characteristic curves and hydraulic conductivity curves via back analysis. Subsequently, numerical analyses were conducted using a 15 m long model for each of the four distinct cover types. The results indicated that the fine gravel significantly improved the barrier performance beyond one-dimensional tests, owing to its high permeability and the influence of the slope. After the capillary barrier failure, the intermediate layers transitioned into efficient drainage layers, particularly in the gravel layer with the highest lateral drainage capacity. Clayey sand at the bottom delayed percolation, thereby supporting the conversion of the intermediate layer into an effective drainage component. Overall, the multi-layer system showed superior percolation performance compared to the clayey sand cover lacking a granular layer.

Divergent effects of food waste derived hydrochar on hydraulic properties and infiltration in a sandy soil

Hydrothermal carbonization of food waste (FW) and application of the solid product (hydrochar, HC) in agriculture, is expected to be an effective way to promote the circulation in food-energy-water nexus. However, little is known on the effects of applying FW-derived HC as a soil amendment on soil properties and functions. The current study investigated the effects of FW-derived HC on soil hydraulic properties and infiltration process. Two types of FW-derived HC were mixed into the top 0–20 cm soil layer at a rate of 1 % and 2 % (w/w) in a soil column experiment with grass growth. In-situ infiltration experiment was conducted within the soil column using a disc infiltrometer at different negative pressure (−20 cm, −6 cm, and −2 cm water head). Unsaturated hydraulic conductivity (K0) and sorptivity (S0) at each negative pressure head was estimated by fitting experimental cumulative infiltration data to a 4-terms infiltration model. Saturated hydraulic conductivity (Ks) and soil water retention curve (SWRC) were measured using undisturbed soil cores taken from the soil column. Soil pore size distribution and soil water constants were estimated from the SWRCs. The results indicate that application of FW-derived HC significantly (p < 0.05) affected soil structure by increasing meso (0.03–0.0002 mm) and macro (> 0.03 mm) pores, which can be attributed to the promoted soil aggregation. Saturated hydraulic conductivity was significantly increased and water retention characteristics were changed with significantly (p < 0.05) increased soil water holding capacity and plant available water content. Besides, FW-derived HCs were found hydrophobic and can render wettable sandy soil water repellent. The K0 and S0 were significantly (p < 0.05) decreased and cumulative infiltration was reduced when soil was wetted from an initial dry condition. Negative effect was found on grass growth with significantly (p < 0.05) decreased aboveground biomass after HCs application. Our study indicates that agriculture application of FW-derived HC as a soil amendment could effectively improve soil fertility by recycling nutrient elements and improving soil structure. However, attention should be paid to the hydrophobicity and phototoxicity of FW-derived HC to avoid the related detrimental influences on plant growth and environment.

Predicting soil hydraulic properties for binary mixtures – concept and application for constructed Technosols

Abstract. Constructed Technosols are an important means of substituting natural soil material, such as peat and geogenic material, for use in urban green infrastructure. One characteristic of Technosols important to their role in urban green infrastructure, specifically with respect to urban water management, is their soil hydraulic properties (SHPs). The SHPs depend on the composition of the constructed Technosols (e.g. their components and their mixing ratio). The diversity of possible components and the infinite number of mixing ratios practically prohibit the experimental identification of the composition needed to achieve suitable soil hydrological functions. In this study, we propose a compositional model for predicting the water retention curves (WRCs) of any binary mixture based on the measured WRCs of its two pure components only (basic scheme) or with one additional mixture (extended scheme). The unsaturated hydraulic conductivity curves (HCCs) are predicted based on the modelled WRCs. The compositional model is developed from existing methods for estimating the porosity of binary mixtures. The model was tested on four data sets of measured WRCs of different binary mixtures. The distribution of water and air in 50 cm high soil columns filled with these mixtures was predicted under hydrostatic conditions in order to assess their suitability for typical urban applications. The difference between the maxima of the pore size distributions ΔPSDmax (m) of the components indicates the applicability of the compositional approach. For binary mixtures with small ΔPSDmax, the water content deviations between the predicted and the measured WRCs range from 0.004 to 0.039 cm3 cm−3. For mixtures with a large ΔPSDmax, the compositional model is not applicable. The prediction of the soil hydraulic properties of any mixing ratio facilitates the simulation of flow and transport processes in constructed Technosols before they are produced (e.g. for specific urban water management purposes).

Extension of the Gardner exponential equation to represent the hydraulic conductivity curve: Inclusion of macropore flow effects

In soil hydraulics, it is crucial to establish an accurate representation of the relative hydraulic conductive curve (rHCC), K_r(h). This paper proposes a simple way to determine K_r(h), called the Modified Gardner Dual model (MGD), using a logarithmic extension of the classical Gardner exponential representation and including macropore flow effects. MGD has five parameters which are hydraulic constants clearly identified in the bilogarithmic representation of K_r(h). Two of them are related to the main inflection point coordinates of rHCC; from them, it is possible to determine the macroscopic capillary length of the infiltration theory. The model was tested in the suction interval 0 < h < 15,000 cm with a total of 249 soil samples from two databases, and employing a flexible representation of the Mualem-van Genuchten (MVG) equation as a reference. Using the RMSE statistics (with log base) to measure the fitting errors, we obtained a 31% reduction in errors (RMSE_MGD = 0.27, RMSE_MVG = 0.39). In 74% of the soils, including samples from the two databases, the reduction was 53% (RMSE_MGD = 0.19, RMSE_MVG = 0.40); the rHCC data fitting of this group was accurate over all the suction h intervals, with RMSE_MGD < 0.32 in each soil sample. In the remaining 26% of the samples, the quality of the MGD fitting degraded due mainly to the presence of multiple rHCC data inflection points. Therefore, in soils without this structural peculiarity, the proposed model revealed to be quite accurate in addition to being analytically simple. Another advantage of MGD is that its parameters depend mainly on the data with h around and lower than the main inflection suction value, which, in turn, never exceeded the 300-cm limit in this study. Hence, in soils that do not have multiple inflections, the extrapolations of the model in drier intervals (1000 cm < h < 15,000 cm) are reliable. The MGD parameter optimization software has been called KUNSAT. It is available in the Supplementary Material or from the corresponding author on request.

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.

The assessment of the effectiveness of the unimodal and bimodal models to evaluate the water flow through nature-based solutions substrates.

Nature-based solutions are popular techniques for managing stormwater. Most of them allow porous media as their main layer. The description of the Soil Water Retention Curve (SWRC) as the Unsaturated Hydraulic Conductivity Curve (UHCC) is often required to run the hydrological simulations with the physically based models. Using the unimodal and bimodal models to assess the SWRC and UHCC of soils is a widespread technique but their evaluation is often present in literature only in terms of curve fitting. Based on these assumptions, this work presents the performance assessment of the van Genuchten unimodal and bimodal models by functional evaluation of them based on the runoff from several substrates. Four substrates were investigated to define the structure, the SWRC, and the UHCC. Results showed that all substrates had a bimodal behaviour with lowest values of RMSE (RMSE_Θ = 0.0023 to 0.0037, RMSE_K = 0.0636 to 0.1284). Finally, a numerical simulation using the HYDRUS-1D model was performed for a three-month data set to check the effectiveness of the unimodal model instead of the bimodal one. The findings have shown that the unimodal model must be preferred instead of the bimodal because it has fewer parameters and assured low discrepancies in runoff volume (ε=0.00% to 6.25%).

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.

Linking Soil Structure, Hydraulic Properties, and Organic Carbon Dynamics: A Holistic Framework to Study the Impact of Climate Change and Land Management

AbstractClimate change and unsustainable land management practices have resulted in extensive soil degradation, including alteration of soil structure (i.e., aggregate and pore size distributions), loss of soil organic carbon, and reduction of water and nutrient holding capacities. Although soil structure, hydrologic processes, and biogeochemical fluxes are tightly linked, their interaction is often unaccounted for in current ecohydrological, hydrological and terrestrial biosphere models. For more holistic predictions of soil hydrological and biogeochemical cycles, models need to incorporate soil structure and macroporosity dynamics, whether in a natural or agricultural ecosystem. Here, we present a theoretical framework that couples soil hydrologic processes and soil microbial activity to soil organic carbon dynamics through the dynamics of soil structure. In particular, we link the Millennial model for soil carbon dynamics, which explicitly models the formation and breakdown of soil aggregates, to a recent parameterization of the soil water retention and hydraulic conductivity curves and to solute and O2 diffusivities to soil microsites based on soil macroporosity. To illustrate the significance of incorporating the dynamics of soil structure, we apply the framework to a case study in which soil and vegetation recover over time from agricultural practices. The new framework enables more holistic predictions of the effects of climate change and land management practices on coupled soil hydrological and biogeochemical cycles.