Soil Water Retention Curve Research Articles

A centrifuge modelling setup for dewatering, rewetting, and characterizing soil deposits

Dewatering in soils serves various purposes in the field of geotechnical engineering as well as other disciplines. In geotechnical engineering, dewatering is effective for improving stability, reducing settlements, and limiting seepage of contaminated pore fluid. This study details the development of an experimental setup that enables dewatering and rewetting soil deposits in-flight in a geotechnical centrifuge while the pore pressures/suctions and moisture contents are measured. Cone penetration tests and shear wave velocity measurements are also used to characterize the material in-flight. Detailed discussion is provided regarding the construction of the tensiometers and calibration of the moisture content probes. Tests are performed on two coal combustion products (CCPs) to exemplify the capabilities of the equipment. Illustrative test results provide time histories and profiles of soil suction and water content at different depths, agreeing well with measurements from laboratory soil-water retention curves. The results also show that dewatering of the CCP deposits leads to significant increases in strength and stiffness, and that after subsequent rewetting, the material maintains a permanent increase in stiffness. The developed system enables an economical and reliable study of the impacts related to dewatering and rewetting cycles to fulfill research purposes in a broad range of engineering disciplines.

FEM Analysis on Isotropic Consolidation of Unsaturated Soils Using Ceramic Disks and Microporous Membrane Filters for Suction Control

Abstract Microporous membrane (MM) filters with a significantly thin thickness have recently been used as an alternative to ceramic disks for matric suction control in the axis-translation technique for element tests on unsaturated soils. Although researchers have highlighted some discrepancies in test results when ceramic disks and MM filters are used, the mechanism leading to these discrepancies has not yet been clarified. In this study, suction-controlled triaxial tests were first conducted on unsaturated completely decomposed granite using two filters. Then, it is natural to believe that the different thicknesses of ceramic disk and MM filters may significantly influence the suction equilibrium process of specimens. To clarify this issue, soil-water-air coupled finite element modeling was performed to simulate the consolidation process in triaxial tests as an initial/boundary value problem. In the numerical calculation, the influence of the void ratio on the water retention curve (WRC) of soil was properly considered by adopting a deformation-dependent WRC model. In particular, the ceramic disk and MM filter, along with the specimens, were modeled, and their differences in thickness were considered. The calculated results show that a nonuniform distribution of suction occurred in the ceramic disk, leading to greater drainage discharge of the specimens. In other words, the influence of the thickness difference between the MM filter and ceramic disk cannot be neglected. However, regardless of whether ceramic discs or MM filters are used, the uneven distribution of the state variables in the specimens may occur in the triaxial tests for unsaturated soils, which emphasizes the importance of assessing the soil tests using the present axis-translation technique as initial/boundary value problems.

Estimation of the Hydraulic Parameters of a Stratified Alluvial Soil in the Region of El Haouareb—Central Tunisia. Experiments, Empirical, Analytical and Inverse Models

Modeling water flow and contaminant transport in the unsaturated zone is a difficult task that relies heavily on good hydrodynamic soil characterization. This article presents a complementarity between experimentation, direct modeling and inverse modeling in order to provide a better estimate of the hydrodynamic parameters of stratified alluvial soil in the El Haouareb region of the Kairouane plain in Tunisia. A field sampling campaign was carried out. The samples collected underwent particle size analysis, bulk density measurements and infiltration tests using a mini-Muntz. In parallel, simple evaporation tests were applied to separate strata. In addition, a 2 m soil column was reconstituted and fitted with sensors to monitor water content, tension, temperature and electrical conductivity. An internal drainage test was performed on this monolith. Three methods were applied using experimental data to estimate soil hydrodynamic parameters. In the first method, pedotransfer functions were used (Rosetta platform) based on granulometric results and bulk density. In the second, water tension and water content monitored during the simple evaporation test were used to plot the soil–water retention curve (SWRC) using SWRC-Fit. In the third method, inverse modeling was applied to the internal drainage test. A comparison of the results showed that the inverse method had the lowest RMSE. Uncertainty analysis has been implemented for both the experimental and numerical set up.

Exploring the Soil Structure Analysis of a Subtropical Watershed by Introducing a Blanket Physico‐Hydraulic Health Index

ABSTRACTMultiple disturbed and undisturbed soil samples were collected to assess the soil structure health in support of water use and management in a subtropical watershed, which is composed of shallow Entisols, near to Canguçu, RS state—Brazil. Various physico‐hydraulic properties related to soil structure, such as bulk density, total, micro and macroporosity, saturated hydraulic conductivity, organic carbon, soil water retention curve (SWRC), and the ranges of plant available water (PAWC) and drainable porosity (DP) were measured to be analyzed. The hydraulic‐energy indices of soil aeration and water retention were determined by integrating the areas under SWRC. Combining these properties with the contents of sand and clay, a unique blanket soil physico‐hydraulic health index (BSPHI) was formulated and applied to investigate the impact of land use and management on the soils of the Ellert Creek Watershed (ECW), thereby supporting soil physical health analysis and, consequently, the analysis of water and human health. The favorable outcomes found in the ECW soils under annual crops indicated good root aeration in the surface layer due to the balance between DP and PAWC. The BSPHI results showed that only about 19% of the areas exhibited unhealthy physical health. The soil structure in pasture areas of ECW was not significantly impacted by the extensive cattle grazing, whereas lands under native forest served as a control for the natural soil quality status. The BSPHI successfully captured small changes in soil structure, regardless of land use and management, highlighting its potential as a promising tool for evaluating soil health by assessing the soil physical quality.

How different soil surface treatments in urban areas affect soil pore structure and associated soil properties and processes

Soil water and temperature regimes in urban environments are greatly affected by different surface treatments, and these may lead even to changes in soil properties. The goal of this study was to find out how soil properties had changed after 8 years of soil surface modification. Five surface treatment scenarios were considered: bare soil (BS), bark chips (BC), concrete paving (CP), mown grass (MG), and unmown grass (UG). X-ray computed tomography and micromorphological analyses showed that character of pores in soils with grass (MG, UG), which were mainly influenced by organisms living in soils and roots, differed from pore character in soil covers BC or CP, which mostly were impacted by organisms living in soils. Both groups differed from BS, which was predominantly affected by the regular treatment consisting in weed removal and soil loosening. Soil under BC was more compact than for other treatments due to decomposition of the bark chips mulch and migration of mulch components. Organic matter content was greatest but its quality lowest in the BC soil, followed by UG, MG, CP, and BS. The highest aggregate stability assessed using the water-stable aggregates (WSA) index was found for UG, followed by MG, BC, CP, and BS. The greatest water retention ability was observed for BC followed by UG, MG, CP, and BS. The unsaturated hydraulic conductivities for pressure head of –2 cm measured for UG, MG, and BS were much higher than were those for CP and BC. Finally, the greatest net CO2 efflux was measured for BC and MG, followed by UG, CP, and BS. CO2 emission correlated negatively with soil physical quality expressed as slope of the soil water retention curve at its inflection point. In general, the best soil conditions were observed for UG. No treatment considerably aggravated soil condition.

Conceptual and Applied Aspects of Water Retention Tests on Tailings Using Columns

The water retention capacity of porous materials is crucial in various geotechnical and environmental engineering applications such as slope stability analysis, landfill management, and mining operations. Filtered tailings stacks are considered an alternative to traditional tailings dams. Nevertheless, the mechanical behaviour and stability of the material under different water content conditions are of concern because these stacks can reach considerable heights. The water behaviour in these structures is poorly understood, particularly the effects of the water content on the stability and potential for liquefaction of the stacks. This study aims to investigate the water retention and flow characteristics of compacted iron ore tailings in high columns to better understand their hydromechanical behaviour. The research used 5 m high columns filled with iron ore tailings from the Quadrilátero Ferrífero region in Minas Gerais, Brazil. The columns were prepared in layers, compacted, and instrumented with moisture content sensors and suction sensors to monitor the water movement during various stages of saturation, drainage, infiltration, and evaporation. The sensors provided consistent data and revealed that the tailings exhibited high drainage capacity. The moisture content and suction profiles were effectively established over time and revealed the dynamic water retention behaviour. The comparison of the data with the theoretical soil water retention curve (SWRC) demonstrated a good correlation which indicates that there was no hysteresis in the material response. The study concludes that the column setup effectively captures the water retention and flow characteristics of compacted tailings and provides valuable insights for the hydromechanical analysis of filtered tailings stacks. These findings can significantly help improve numerical models, calibrate material parameters, and contribute to the safer and more efficient management of tailings storage facilities.

Applying a Comprehensive Model for Single-Ring Infiltration: Assessment of Temporal Changes in Saturated Hydraulic Conductivity and Physical Soil Properties

Modeling agricultural systems, from the point of view of saving and optimizing water, is a challenging task, because it may require multiple soil physical and hydraulic measurements to investigate the entire crop cycle. The Beerkan method was proposed as a quick and easy approach to estimate the saturated soil hydraulic conductivity, Ks. In this study, a new complete three-dimensional model for Beerkan experiments recently proposed was used. It consists of thirteen different calculation approaches that differ in estimating the macroscopic capillary length, initial (θi) and saturated (θs) soil water contents, use transient or steady-state infiltration data, and different fitting methods to transient data. A steady-state version of the simplified method based on a Beerkan infiltration run (SSBI) was used as the benchmark. Measurements were carried out on five sampling dates during a single growing season (from November to June) in a long-term experiment in which two soil management systems were compared, i.e., minimum tillage (MT) and no tillage (NT). The objectives of this work were (i) to test the proposed new model and calculation approaches under real field conditions, (ii) investigate the impact of MT and NT on soil properties, and (iii) obtain information on the seasonal variability of Ks and other main soil physical properties (θi, soil bulk density, ρb, and water retention curve) under MT and NT. The results showed that the model always overestimated Ks compared to SSBI. Indeed, the estimated Ks differed by a factor of 11 when the most data demanding (A1) approach was considered by a factor of 4–8, depending on the transient or steady-state phase use, when A3 was considered and by a practically negligible factor of 1.0–1.9 with A4. A relatively higher seasonal variability was detected for θi at the MT than NT system. Under both MT and NT, ρb did not change between November and April but increased significantly until the end of the season. The selected calculation approaches provided substantially coherent information on Ks seasonal evolution. Regardless of the approach, the results showed a temporal stability of Ks at least from early April to June under NT; conversely, the MT system was, overall, more affected by temporal changes with a relative stability at the beginning and middle of the season. These findings suggest that a common sampling time for determining Ks could be set at early spring. Soil management affected the soil properties, because the NT system was significantly wetter and more compact than MT on four out of five dates. However, only NT showed a significantly increasing correlation between Ks and the modal pore diameter, suggesting the presence of a relatively smaller and better interconnected pore network in the no-tilled soil. This study confirms the need to test infiltration models under real field conditions to evaluate their pros and cons. The Beerkan method was effective for intensive soil sampling and accurate field investigations on the temporal variability of Ks.

Estimating soil water retention curves over the entire saturation range: A thermal conductivity-based approach

The relationship between soil water content (θ) and suction (h, referring to the absolute value of pressure head), is described by the soil water retention curve (SWRC). Our earlier research (Fu et al. (2021, 2023a) [J. Hydrol.127171]; [J. Hydrol.129898]) recognized underlying correlations between SWRCs and soil thermal conductivity (λ) versus θ curves, and developed methodologies to ascertain the parameters of the van Genuchten (vG) equation using λ(θ) measurements, described by the Ghanbarian & Daigle (GD) equation, and basic soil characteristics. Limitations intrinsic to the van Genuchten equation restrict the GD-vG approach to generate precise estimates only in the wet and medium suction range, specifically h ranging from 0 to 150 mH2O. The validity of these approaches in the dry region remains uncertain. In this study, we associated the Peters-Durner-Iden (PDI) model parameters to those of the GD model. An initial examination was performed on the linearization processes needed to derive the hydraulic continuity water content (θhc) from the capillary water component as characterized by the PDI model and to choose the suction at oven dryness (h0) based on PDI model performance. Subsequently, two piecewise functions and two pedo-transfer functions were formulated to compute the PDI model parameters utilizing soil porosity, particle size distribution, and GD parameters, based on a calibration dataset comprising 25 different soils. The new GD-PDI approach was subsequently assessed with six independent soils and juxtaposed with the previous GD-vG approach. The GD-PDI approach outperformed the GD-vG approach, particularly within the dry range.

A pore network modeling approach to bridge void ratio‐dependent soil water retention and unsaturated hydraulic conductivity curves

AbstractIn geotechnical engineering, understanding the relationship between soil permeability and deformation is essential, particularly for applications like earth dams, where compaction‐induced permeability reduction is crucial for performance optimization. In unsaturated soils, soil moisture content significantly impacts hydraulic conductivity. Traditionally, changes in unsaturated hydraulic conductivity have been linked to soil void ratios. However, a pore network modeling perspective reveals the significance of structural parameters, such as pore and throat size distribution and pore coordination number. This study introduces a pore network model to estimate unsaturated hydraulic conductivity based on void ratio‐dependent soil‐water retention curves. It examines how soil deformation at varying stress levels affects structural parameters and the water phase continuity. The model shows strong potential in predicting unsaturated hydraulic conductivity across different stress levels, aligning well with experimental data and established equations. Notably, the aspect ratio and coordination number parameters are most affected by stress levels. The study also presents relationships to describe changes in pore network structural parameters with soil void ratio, which can be used to predict soil‐water retention curves and unsaturated hydraulic conductivity at various stress levels.

Measuring shrinkage of expansive soils using a novel automated high‐frequency setup

AbstractSoil shrinkage characteristic curves are used to describe the shrinkage behavior and hydraulic properties of unsaturated soils. To construct soil shrinkage characteristic curves, a high‐data‐density measurement method is needed that relates water content to soil volume changes. We present a fully automated soil shrinkage measurement setup, based on the simplified evaporation method, to characterize the shrinkage behavior of undisturbed natural expansive clay soils. The high data density creates the opportunity to produce soil shrinkage characteristic curves without the need for a mathematical model. The technique allows for resaturation of the samples after drying, enabling differentiation between reversible and irreversible shrinkage. The setup consists of the commercialized HYPROP2 apparatus combined with optical distance sensors to measure the horizontal and vertical dimensions of the samples, yielding data on the sample volume, weight, and soil water suction. The measurement frequency is once per 10 min, and the measurement period is up to 4 weeks, providing a detailed time series of the drying and shrinkage characteristics. The setup can capture the different shrinkage phases and offers the opportunity to relate the soil shrinkage characteristic curve to soil water suction. The measurement data acquisition rate and accuracy enable detailed interpretation of soil water retention curves for nonrigid soils and are shown to be essential for understanding and quantifying the shrinkage potential of several types of deposits.

Evaluation of Pedotransfer Functions to Estimate Soil Water Retention Curve: A Conceptual Review

The soil water retention curve (SWRC) is a vital soil property used to evaluate the soil’s water holding capacity, a critical factor in various applications such as determining soil water availability for plants, soil conservation and management, climate change adaptation, and mitigation of flood risks. Estimating SWRC directly in the field and laboratory is a time-consuming and laborious process and requires numerous instruments and measurements at a specific location. In this context, various estimation approaches have been developed, including pedotransfer functions (PTFs), over the past three decades to estimate soil water retention and its associated properties. Despite the efficiencies, PTFs and semi-physical approach-based models often have several limitations, particularly in the dry range of the SWRC. PTFs-based modeling has become a key research topic due to readily available soil data and cost-effective methods for deriving essential soil parameters, which enable more efficient decision-making in sustainable land-use management. Therefore, advancement and adjustment are necessary for reliable estimations of the SWRC from readily available data. This article reviews the evaluation of the current and past PTFs for estimating the SWRC. This study aims to evaluate PTF techniques and semi-physical approaches based on soil texture, bulk density, porosity, and other related factors. Additionally, it also assesses the performance and limitations of various common semi-physical models proposed and developed by Arya and Paris, Haverkamp and Parlange, the Modified Kovács model by Aubertin et al., Chang and Cheng, Meskini-Vishkaee et al., Vidler et al., and Zhai et al. This assessment will be effective for researchers in this field and provide valuable insight into the importance of new PTFs for modeling SWRC.

Effect of salt concentration on osmotic potential in drying soils—Measurement and models

AbstractThe water potential in drying soils, comprising both matric potential and osmotic potential components, can be measured using the dew point method (DPM). By combining DPM data with retention curve data acquired from techniques such as the suction plate method or the simplified evaporation method (SEM), it becomes possible to determine the soil water retention curve across the entire moisture spectrum. However, as the latter methods only determine the matric potential, the osmotic potential component in DPM data must either be negligible or known so that osmotic and matric potential components can be separated. This study aims to critically analyse common approaches for calculating the osmotic potential. To achieve this, we measured the water retention properties of a silt loam, a sandy loam and a sand across the entire moisture range by combining SEM and DPM. By using almost salt‐free soil material, we characterized reference water retention curves with negligible osmotic potential components. The impact of salt on water potential was analysed by conditioning soils with MgCl2 solutions of different concentrations, drying them, and measuring the water potential at different water contents using the DPM. The resulting water potentials were compared to the reference potentials and differences were interpreted as the osmotic potential component. The DPM‐measured water potentials in drying soils can be significantly affected by osmotic potential, especially at higher matric potentials (low suctions). Two models accounting for ideal and one model accounting for non‐ideal electrolyte behaviour were used to compare osmotic potential predictions with measurements. At low to medium salt concentrations, all models performed fairly well. At high concentrations, only the model accounting for non‐ideal behaviour predicted the osmotic potential satisfactorily, whereas at very high concentrations, all models underestimated the impact of osmotic potential on water potential. This suggests that the surface properties of the soil matrix, such as the specific surface area and surface charges, may lead to a decrease in osmotic potential beyond what is expected in pure solutions.

A state-dependent bimodal soil water retention model considering evolutions of pore-scale soil–water interaction

Existing soil water retention models can capture the effects of void ratio and/or net mean stress on the changes in soil water retention curve (SWRC). However, most of them, especially those for bimodal soil that is composed of macro- and micro-pore systems, do not consider the simultaneous evolutions of pore structure and particle contacts, and neglected the drainage processes of bulk water and meniscus water when subjected to any external hydro-mechanical loading condition. This paper proposes a state-dependent bimodal SWRC model considering pore structure evolution and pore-scale water drainage. In this model, an elastoplastic relationship is introduced to describe the variations of void ratio with net mean stress and suction, whilst the corresponding changes in the pore structure are modelled by the shift and scaling of the pore-size distribution and the volumetric proportion between the macro- and micro-pore systems. By relating stress-induced changes in void ratio and coordination number to the mean pore radius and total meniscus water volume, respectively, the soil water content at different stages of water drainage can be determined. The model is validated against existing and new data, and can capture the effects of void ratio, isotropic stress, and K0 stress for unimodal and bimodal soils.

The added value of CoreVESS score and penetration resistance in predicting the soil–water retention curve

Soil hydraulic properties, such as the soil–water retention curve (SWRC), play a crucial role in simulating water transport within the vadose zone. They can be directly measured in the laboratory or field, or indirectly predicted by using Pedotransfer Functions (PTFs). It has been advocated that accounting for soil structure could marginally improve the prediction accuracy of PTFs in its wet range. The aim of this study was to test whether the use of an easy-to-determine soil structure-related variable could effectively increase the prediction accuracy of water retention curve PTFs. Additionally, we explored whether including soil strength, another soil property that is easy and quickly to determine in the field in many replicates, could improve PTF performance. Our investigation involved extensive sampling of 252 soil horizons across 42 cropped fields, each exhibiting within-field variations in soil structural degradation. Soil structure was represented by a soil structural quality score (Sq) taken from a semi-quantitative visual soil assessment (CoreVESS), which was included in the PTFs as predictor variable or as discriminator (Sq) for data grouping in addition to or instead of ‘classical’ predictors such as soil organic carbon content, clay and sand content, bulk density and soil layer. Penetration resistance (PR) was taken as a variable for soil strength. Various regression methods from classical regression to machine learning were evaluated on Train and Test subdatasets. To evaluate the effect of data grouping, the dataset was split in two based on structure and on texture. Both point and parametric PTFs were developed, with the latter predicting the parameters of the Van Genuchten (VG) water retention equation. Overall, the k-nearest neighbors PTFs (KNN-PTFs) achieved the best performance showing the highest R2 and lowest RMSE values. Moreover, the point PTFs showed much better and more stable predictive ability than parametric PTFs. The inclusion of PR as an additional predictor variable slightly improved the prediction performance in the wet range of the SWRC. Interestingly, Sq proved beneficial when substituted for bulk density but not when used in addition to it. Grouping data based on texture and structure improved the prediction accuracy of the PTFs, particularly evident in textural grouping and MLR, but not in KNN.

Modeling soil water dynamics of an intensively cultivated histosol

AbstractUnderstanding soil–water dynamics in cultivated organic soils (histosols) is crucial for sustainable agriculture, ecosystem preservation, and climate change mitigation. Data on the soil water retention curves (SWRCs) of these histosols in Canada are not readily available in the literature. The Hydrus‐1D model was used to predict SWRCs for a cultivated organic soil in Quebec. The model was validated with matric potential measured in a soil column at 10, 26, and 48 cm. The optimized hydraulic parameters in the Hydrus model resulted in field capacity moisture contents of 0.6, 0.696, and 0.49 cm3 cm−3 at the three depths, respectively. Wilting point moisture contents of 0.165, 0.107, and 0.14 cm3 cm−3 were obtained at the same depths. Model accuracy was confirmed with a Nash–Sutcliffe efficiency of 0.76, 0.91, and 0.78 for the three depths, respectively. These findings can inform decisions regarding water management practices such as irrigation and drainage requirements for intensively cultivated organic soils.

Effective approach to predict soil-water retention curve of bentonites considering adsorption and capillarity

Understanding soil–water retention behavior is a longstanding topic. Water retained in soils can be decomposed into adsorptive and capillary components, controlled by different physicochemical mechanisms. The capillary water retention was frequently discussed in literatures, but the adsorption role was rarely considered, especially for high active clays (e.g., bentonite). In this study, two novel equations for quantifying adsorptive and capillary water retentions are proposed, generating a twofold model to continuously simulate soil–water retention behavior of bentonites. Only 5 parameters, i.e., the maximum adsorption capacity, characteristic adsorptive and capillary suctions, and uniformities of adsorptive and capillary pores, are defined with clear physical meanings. A graphical method was suggested to firstly determine the maximum adsorption capacity, and the remaining parameters are efficiently estimated by non-linear curve fitting. The water retention data for various bentonites, representing a variety of hydration conditions, initial compactness, montmorillonite content and suction range, are used to assess the model performance. The predictions agree well with the measured total, adsorptive and capillary water contents of a Wyoming bentonite, and the fitting curves also match well with test data for other bentonites over the full suction range. Adsorption parameters are distributed within a narrow range, while the characteristic capillary suction is also distributed within a limited range under constant volume condition. The proposed model allows reasonable predictions about the capillary onset and the transition from adsorption to capillarity, revealing obvious superiority by comparison with other hybrid models. This work offers a new pathway to quantitively assess the soil–water retention curve of high active clays.

Implementation of a slope stability method in the CRITERIA-1D agro-hydrological modeling scheme

AbstractThis paper presents the implementation of a slope stability method for rainfall-induced shallow landslides in CRITERIA-1D, which is an agro-hydrological model based on Richards’ equation for transient infiltration and redistribution processes. CRITERIA-1D can simulate the presence and development of roots and canopies over space and time, the regulation of transpiration activity based on real meteorological data, and the evaporation reduction caused by canopies. The slope can be considered composed of a multi-layered soil, leading to the possibility of simulating the bedrock and of setting an initial water table level. CRITERIA-1D can consider different soil horizons characterized by different hydraulic conductivities and soil water retention curves, thus allowing the simulation of capillarity barriers. The validation of the proposed physically based slope stability model was conducted through the simulation of the collected water content and water potential data of an experimental slope. The monitored slope is located close to Montuè, in the north-eastern sector of Oltrepò Pavese (northern Apennines—Italy). Just close to the monitoring station, a shallow landslide occurred in 2014 at a depth of around 100 cm. The results show the utility of agro-hydrological modeling schemes in modeling the antecedent soil moisture condition and in reducing the overestimation of landslides events detection, which is an issue for early warning systems and slope management related to rainfall-induced shallow landslides. The presented model can be used also to test different bioengineering solutions for slope stabilization, especially when data about rooting systems and plant physiology are known.

Root Influences Rhizosphere Hydraulic Properties through Soil Organic Carbon and Microbial Activity.

Processes of water retention and movement and the hydraulic conductivity are altered in the rhizosphere. The aim of this study was to investigate the physical-hydric properties of soil aggregates in the rhizosphere of annual ryegrass (Lolium multiflorum) cropped in a Kandiudalfic Eutrudox, taking into account aspects related to soil aggregate stability. Soil aggregates from rhizosphere soil (RZS) and soil between plant rows (SBP) were used to determine soil water retention curves (SWRCs) and saturated hydraulic conductivity (Ksat). In addition, properties related to soil aggregate stability, such as water-dispersible clay, soil organic carbon (SOC), and microbial activity, were also assessed. The higher microbial activity observed in the RZS was facilitated by increased SOC and microbial activity, resulting in improved soil aggregation (less water-dispersible clay). For nearly all measured matric potentials, RZS had a higher water content than SBP. This was attributed to the stability of aggregates, increase in SOC content, and the root exudates, which improved soil water retention. The increase in total porosity in RZS was associated with improved soil aggregation, which prevents deterioration of the soil pore space and results in higher Ksat and hydraulic conductivity as a function of the effective relative saturation in RZS compared to SBP.

Estimating the soil water retention curve by the HYPROP-WP4C system, HYPROP-based PCNN-PTF and inverse modeling using HYDRUS-1D

Due to the difficulty of direct field measurement, soil hydraulic properties are often obtained in laboratory settings using small undisturbed soil samples or estimated indirectly through pedotransfer functions (PTFs). The pseudo-continuous pedotransfer function (PCNN-PTF) is a neural network-based approach for estimating soil hydraulic properties. The main objective of this study was to use field soil moisture and tension data to assess soil water retention curves obtained from the HYPROP-WP4C system and the HYPROP-based PCNN-PTF. The in-situ soil water retention data were simultaneously acquired using Acclima TDT and MeterGroup MPS-6 sensors (every 30 min, May to September 2020) from 11 hybrid bermudagrass plots under different irrigation treatments in Riverside, California. We utilized extended evaporation and dewpoint methods using HYPROP-WP4C (Meter Group Inc., USA) devices to obtain lab-measured SWRCs. In addition, SWRCs were estimated from Rosetta (Schaap et al., 2001) and by inverse modeling in HYDRUS-1D utilizing in-situ moisture retention data. Although the hysteresis impacted field data, overall, there was a good agreement between the in-situ and lab water retention data for most samples, especially within the pF range of 2–3.5. The PCNN-PTF outperformed Rosetta in estimating laboratory (RMSE=0.034 cm3 cm−3 vs 0.063 cm3 cm−3) and in-situ soil moisture data (RMSE=0.048 cm3 cm−3 vs 0.082 cm3 cm−3). Inverse modeling of in-situ data also performed well in estimating the SWRC (RMSE=0.043 cm3 cm−3); however, further attention is required in dry and saturated soil conditions. We developed a simple, free, and easy-to-access tool called PC-PTF for estimating the SWRC using the PCNN-PTF model evaluated in this study. The PC-PTF can be accessed from the Verdi Water Management Group website: http://www.ucrwater.com/software-and-tools.html.

Insight into the variation of soil hydraulic properties under beech and spruce forest—A case study in the forest of Tharandt, NE Germany

AbstractBackgroundThe increasing vulnerability of forests in the temperate zone due to climate change has led to modification in the forest structure to secure woody raw materials and ecosystem benefits. Such changes will influence hydrological processes both at the stand and catchment scale. Soil hydraulic properties (SHP) play an important role in assessing the water cycle in these ecosystems. Yet, knowledge regarding the effect of forest‐ and site‐specific conditions on SHP in temperate climates is scarce.AimsThis work addresses this research gap by assessing the variation of SHP under two common European forest stands, Fagus sylvatica and Picea abies (1) with comparable site conditions, and (2) across differing site conditions.MethodsWe determined the soil water retention curve (WRC) and the hydraulic conductivity curve (HCC) in several plots with the bimodal Kosugi–Mualem's hydraulic model. These functions were determined using combined field and laboratory measurements, including hydraulic conductivity and water content from soil samples.Results(1) We observed distinct variations in SHP between beech and spruce forest stands with comparable site conditions; however, no clear pattern in the variation was discernible. (2) A noticeable effect of the site‐specific characteristics on the SHP was detected. Moreover, SHP in each analysed forest type presented individual variations.ConclusionsThis study demonstrates that SHP present a wide range of variations in terms of both forest‐ and site‐specific conditions. Hence, due to its heterogeneity, we emphasise the need for more research to better characterise SHP in temperate zone forests. Moreover, this study underlines the urgent use of a minimum set of parameters in studies when addressing SHP (e.g., tree age, soil texture).