Pore-morphology-based simulations of drainage and wetting processes in porous media

2011 ◽  
Vol 42 (2-3) ◽  
pp. 128-149 ◽  
Author(s):  
T. P. Chan ◽  
Rao S. Govindaraju

Soil hydraulic properties relating saturation, water pressure, and hydraulic conductivity are known to exhibit hysteresis. In this paper, we focus on the determination of the water retention curve for a porous medium through a novel pore-scale simulation technique that is based on mathematical morphology. We develop an algorithm that allows for the representation of three-dimensional randomly packed porous media of any geometry (i.e. not restricted to idealized geometries such as spherical or ellipsoidal particles/pore space) so that the connectivity-, tortuosity-, and hysteresis-causing mechanisms are represented in both drainage and wetting processes, and their role in determining macroscopic fluid behavior is made explicit. Using this method, we present simulation results that demonstrate hysteretic behavior of wetting and non-wetting phases during both drainage and wetting cycles. A new method for computing interfacial surface areas is developed. The pore-morphology-based method is critically evaluated for accuracy, sample size effects, and resolution effects. It is found that the method computes interfacial areas more accurately than existing methods and allows for (i) examination of relationships between water pressure, saturation and interfacial area for hysteretic soils, and (ii) comparisons with previously developed theoretical models of soil hydraulic properties. The pore-morphology-based method shows promise for applications in vadose zone hydrology.

2020 ◽  
Author(s):  
Patrizia Hangele ◽  
Katharina Luise Müller ◽  
Hannes Laermanns ◽  
Christina Bogner

<p>The need to study the occurrence and effects of microplastic (MP) in different ecosystems has become apparent by a variety of studies in the past years. Until recently, research regarding MP in the environment has mainly focused on marine systems. Within terrestrial systems, studies suggest soils to be the biggest sink for MP. Some studies now started to explore the presence of MP in soils. However, there is a substantial lack of the basic mechanistic understanding of the behaviour of MP particles within soils.</p><p>This study investigates how the presence of MP in soils affects their hydraulic properties. In order to understand these processes, experiments are set up under controlled laboratory conditions as to set unknown influencing variables to a minimum. Different substrates, from simple sands to undisturbed soils, are investigated in soil cylinders. MP particles of different sizes and forms of the most common plastic types are applied to the surface of the soil cylinders and undergo an irrigation for the MP particles to infiltrate. Soil-water retention curves and soil hydraulic conductivity are measured before and after the application of MP particles. It is hypothesised that the infiltrated MP particles clog a part of the pore space and should thus reduce soil hydraulic conductivity and change the soil-water retention curve of the sample. Knowledge about the influence of MP on soil hydraulic properties are crucial to understand transport and retention of MP in soils.</p>


2020 ◽  
Author(s):  
Joseph Pollacco ◽  
Jesús Fernández-Gálvez ◽  
Sam Carrick

<p>Indirect methods for estimating soil hydraulic properties from particle size distribution have been developed due to the difficulty in accurately determining soil hydraulic properties, and the fact that particle size distribution is one piece of basic soil physical information normally available. The similarity of the functions describing the cumulative distribution of particle size and pore size in the soil has been the basis for relating particle size distribution and the water retention function in the soil. Empirical and semi-physical models have been proposed, but these are based on strong assumptions that are not always valid. For example, soil particles are normally assumed to be spherical, with constant density regardless of their size; and the soil pore space has been described by an assembly of capillary tubes, or the pore space in the soil matrix is assumed to be arranged in a similar way regardless of particle size. However, in a natural soil the geometry of the pores may vary with the size of the particles, leading to a variable relation between particle radius and pore radius.</p><p> </p><p>The current work is based on the hypothesis that the geometry of the pore size and the void ratio depends on the size of the soil particles, and that a physically based model can be generalised to predict the water retention curve from particle size distribution. The rearrangement of the soil particles is considered by introducing a mixing function that modulates the cumulative particle size distribution, while the total porosity is constrained by the saturated water content.</p><p> </p><p>The model performance is evaluated by comparing the soil water retention curve derived from laboratory measurements with a mean Nash–Sutcliffe model efficiency a value of 0.92 and a standard deviation of 0.08. The model is valid for all soil types, not just those with a marginal clay fraction.</p>


2018 ◽  
Vol 66 (2) ◽  
pp. 161-169 ◽  
Author(s):  
Camila R. Bezerra-Coelho ◽  
Luwen Zhuang ◽  
Maria C. Barbosa ◽  
Miguel Alfaro Soto ◽  
Martinus Th. van Genuchten

AbstractMany soil, hydrologic and environmental applications require information about the unsaturated soil hydraulic properties. The evaporation method has long been used for estimating the drying branches of the soil hydraulic functions. An increasingly popular version of the evaporation method is the semi-automated HYPROP©measurement system (HMS) commercialized by Decagon Devices (Pullman, WA) and UMS AG (München, Germany). Several studies were previously carried out to test the HMS methodology by using the Richards equation and the van-Genuchten-Mualem (VG) or Kosugi-Mualem soil hydraulic functions to obtain synthetic data for use in the HMS analysis, and then to compare results against the original hydraulic properties. Using HYDRUS-1D, we carried out independent tests of the HYPROP system as applied to the VG functions for a broad range of soil textures. Our results closely agreed with previous findings. Accurate estimates were especially obtained for the soil water retention curve and its parameters, at least over the range of available retention measurements. We also successfully tested a dual-porosity soil, as well as an extremely coarse medium with a very high van Genuchtennvalue. The latter case gave excellent results for water retention, but failed for the hydraulic conductivity. In many cases, especially for soils with intermediate and highnvalues, an independent estimate of the saturated hydraulic conductivity should be obtained. Overall, the HMS methodology performed extremely well and as such constitutes a much-needed addition to current soil hydraulic measurement techniques.


2020 ◽  
Author(s):  
Andreas Cramer ◽  
Ursula Bundschuh ◽  
Pascal Bernard ◽  
Mohsen Zarebanadkouki ◽  
Andrea Carminati

<p>Soils are the largest sink of microplastic particles (MPP) in terrestrial ecosystems. However, there is little knowledge on the implication of MPP contaminating soils. In particular, we don’t know how MPP move and, on the other hand, how they affect soil hydraulic properties and soil moisture dynamics.</p><p>Among the expected effects of MPP on soil hydraulic properties is the likelihood that MPP enhances soil water repellency. This emerges from (1) the MPP surface chemical properties as well as (2) their surface physical properties like size and shape. Here, we tested mixtures of MPP and a model porous media. The Sessile Drop Method was applied and apparent contact angles were measured. We are able to show enlarged contact angles with rising concentrations of MPP. Already in relatively low concentrations of MPP the contact angels exhibit a steep increase and are rapidly reaching areas of super-hydrophobicity. Furthermore, we provide the physical explanation of the apparent contact angles resulting from the three-phase contact line between solid composite surfaces, water and air. The considered modes of a droplet lying on a surface are Wenzel, Cassie-Baxter and Young. The goal here was to differentiate between the involved surfaces building up the apparent contact angle and to pin down the impact of MPP in these systems.</p><p>Thinking about the implications of these results, an increased water repellency alters soil hydraulic properties towards less water content resulting in a shift in the water retention curve. Less water in soils especially at sites of high MPP concentrations leads to a limitation of degradation of MPP by hydrolysis. Additionally, microorganisms themselves and their enzymes cannot migrate in the liquid phase towards the MPP even elongating the process of natural purification.</p>


2020 ◽  
Author(s):  
Mirko Castellini ◽  
Simone Di Prima ◽  
Anna Maria Stellacci ◽  
Massimo Iovino ◽  
Vincenzo Bagarello

<p>Testing new experimental procedures to assess the effects of the drops impact on the soil sealing formation is a main topic in soil hydrology.</p><p>In this field investigation, the methodological approach proposed first by Bagarello et al. (2014) was extended to account for a greater soil infiltration surface (i.e., about 3.5 times higher), a higher range and number of heights of water pouring and to evaluate the different impact on soil management. For this purpose, the effects of three water pouring heights (low, L=3 cm; medium, M=100 cm; high, H=200 cm) on both no-tilled (NT) and conventionally tilled (CT) loam soil were investigated by Beerkan infiltration runs and using the BEST-procedure of data analysis to estimate the soil hydraulic properties.</p><p>Final infiltration rate decreased when perturbing runs (i.e., M and H) were carried out as compared with the non-perturbing (L) ones (by a factor of 1.5-3.1 under NT and 3.4-4.4 under CT). Similarly, the water retention scale parameter, h<sub>g</sub>, increased (i.e., higher in absolute terms) by a factor 1.6-1.8 under NT and by a factor 1.7 under CT. Saturated hydraulic conductivity, K<sub>s</sub>, changed significantly as a function of the increase of water pouring height; regardless of the soil management, perturbing runs caused a reduction in soil permeability by a factor 5 or 6. Effects on hydraulic functions (i.e., soil water retention curve and hydraulic conductivity function), obtained with the BEST-Steady algorithm, were also highlighted. For instance, differences in water retention curve at fixed soil pressure head values (i.e., field capacity, FC, and permanent wilting point, PWP) due to perturbing and non-perturbing runs, were estimated as higher under NT (3.8%) than CT (3.4%) for FC, and equal to 2.1% or 1.6% for PWP.</p><p>Main results of this investigation confirm that a recently tilled loamy soil, without vegetation cover, can be less resilient as compared to a no-tilled one, and that tested water pouring heights methodology looks promising to mimic effects of high energy rainfall events and to quantify the soil sealing effects under alternative management of the soil.</p><p><strong>Acknowledgments</strong></p><p>The work was supported by the project “STRATEGA, Sperimentazione e TRAsferimento di TEcniche innovative di aGricoltura conservativA”, funded by Regione Puglia–Dipartimento Agricoltura, Sviluppo Rurale ed Ambientale, CUP: B36J14001230007.</p><p><strong> </strong><strong>References</strong></p><p>Bagarello, V., Castellini, M., Di Prima, S., Iovino, M. 2014. Soil hydraulic properties determined by infiltration experiments and different heights of water pouring. Geoderma, 213, 492–501. https://doi.org/10.1016/j.geoderma.2013.08.032</p>


2020 ◽  
Author(s):  
Mahyar Naseri ◽  
Sascha C. Iden ◽  
Wolfgang Durner

<p>Stony soils are soils that contain a high amount of stones and are widespread all over the world.  The effective soil hydraulic properties (SHP), i.e. the water retention curve (WRC) and the hydraulic conductivity curve (HCC) are influenced by the presence of stones in the soil. This influence is normally neglected in vadose zone modeling due to the considerable measurement challenges in stony soils. The available data on the effect of stones on SHP is scarce and there is not a systematic modeling approach to obtain the effective SHP in stony soils. Most of the past studies are limited to the effect of stones on the WRC and saturated hydraulic conductivity and low and medium stone contents (up to 40 % v/v). We investigated the effect of stone content on the effective SHP of stony soils through a series of evaporation experiments. Two soil materials a) sandy loam and b) silt loam as background soils were packed with different volumetric contents (0, 10, 30 and 60 %) of medium stones were in containers with a volume of 5060 cm<sup>3</sup>. Volumetric stone contents were chosen in a way to present stone-free, moderately stony and highly stony soils. All of the experiments were carried out in two replicate packings with an almost identical bulk density. Packed samples were saturated with water from the bottom and subjected to evaporation in a climate-controlled room. During the evaporation experiments, the pressure head and soil temperature were continuously monitored and the water loss from the soil columns was measured with a balance. The dewpoint method provided additional data on the WRC in the dry soil. The resulting data were evaluated by inverse modeling with the Richards equation to identify effective SHP and to analyze the effect of stone content on the evaporation rate, soil temperature, the effective WRC and the effective HCC. The applied methodology was successful in identifying effective SHP with high precision over the full moisture range. The results reveal a quicker transition from stage I to stage II of evaporation in highly stony soils. Evaporation rate reduces with the increase of the volumetric stone content. The existence of a high amount of stone content shorten stage II of evaporation driven by the vapor diffusion through the restricted soil evaporative surface.</p>


2021 ◽  
Vol 69 (4) ◽  
pp. 378-386
Author(s):  
Hongjie Guan ◽  
Xinyu Liu

Abstract The presence of biocrusts changes water infiltration in the Mu Us Desert. Knowledge of the hydraulic properties of biocrusts and parameterization of soil hydraulic properties are important to improve simulation of infiltration and soil water dynamics in vegetation-soil-water models. In this study, four treatments, including bare land with sporadic cyanobacterial biocrusts (BL), lichen-dominated biocrusts (LB), early-successional moss biocrusts (EMB), and late-successional moss biocrusts (LMB), were established to evaluate the effects of biocrust development on soil water infiltration in the Mu Us Desert, northwest of China. Moreover, a combined Wooding inverse approach was used for the estimation of soil hydraulic parameters. The results showed that infiltration rate followed the pattern BL > LB > EMB > LMB. Moreover, the LB, EMB, and LMB treatments had significantly lower infiltration rates than the BL treatment. The saturated soil moisture (θs ) and shape parameter (α VG) for the EMB and LMB treatments were higher than that for the BL and LB treatments, although the difference among four treatments was insignificant. Water retention increased with biocrust development at high-pressure heads, whereas the opposite was observed at low-pressure heads. The development of biocrusts influences van Genuchten parameters, subsequently affects the water retention curve, and thereby alters available water in the biocrust layer. The findings regarding the parameterization of soil hydraulic properties have important implications for the simulation of eco-hydrological processes in dryland ecosystems.


1999 ◽  
Author(s):  
Arthur Warrick ◽  
Uri Shani ◽  
Dani Or ◽  
Muluneh Yitayew

The primary information for accurately predicting water and solute movement and their impact on water quality is the characterization of soil hydraulic properties. This project was designed to develop methods for rapid and reliable estimates of unsaturated hydraulic properties of the soil. Particularly, in situ methodology is put forth, based on subsurface point sources. Devices were designed to allow introduction of water in subsurface settings at constant negative heads. The ability to operate at a negative head allows a direct method of finding unsaturated soil properties and a mechanism for eliminating extremely rapid preferential flow from the slow matrix flow. The project included field, laboratory and modeling components. By coupling the measurements and the modeling together, a wider range of designs can be examined, while at the same time realistic performance is assured. The developed methodology greatly expands the possibilities for evaluating hydraulic properties in place, especially for measurements in undisturbed soil within plant rooting zones. The objectives of the project were (i) To develop methods for obtaining rapid and reliable estimates of unsaturated hydraulic properties in situ, based on water distribution from subsurface point sources. These can be operated with a constant flow or at a constant head; (ii) To develop methods for distinguishing between matrix and preferential flow using cavities/permeameters under tension; (iii) To evaluate auxiliary measurements such as soil water content or tensions near the operating cavities to improve reliability of results; and (iv: To develop numerical and analytical models for obtaining soil hydraulic properties based on measurements from buried-cavity sources and the auxiliary measurements. The project began in July 1995 and was terminated in November 1998. All of the objectives were pursued. Three new subsurface point sources were designed and tested and two old types were also used. Two of the three new designs used a nylon cloth membrane (30 mm) arranged in a cylindrical geometry and operating at a negative water pressure (tension). A separate bladder arrangement allowed inflation under a positive pressure to maintain contact between the membrane and the soil cavity. The third new design used porous stainless steel (0.5 and 5 mm) arranged in six segments, each with its own water inlet, assembled to form a cylindrical supply surface when inflated in a borehole. The "old" types included an "off-the-shelf" porous cup as well as measurements from a subsurface drip emitter in a small subsurface cavity. Reasonable measurements were made with all systems. Sustained use of the cloth membrane devices were difficult because of leaks and plugging problems. All of the devices require careful consideration to assure contact with the soil system. Steady flow was established which simplified the analysis (except for the drip emitter which used a transient analysis).


2004 ◽  
Vol 8 (1) ◽  
pp. 2-7 ◽  
Author(s):  
A. Bayer ◽  
H.-J. Vogel ◽  
K. Roth

Abstract. X-ray absorption measurements have been explored as a fast experimental approach to determine soil hydraulic properties and to study rapid dynamic processes. As examples, the pressure-saturation relation θ(Ψ) for a uniform sand column has been considered as has capillary rise in an initially dry sintered glass column. The θ(Ψ)-relation is in reasonable agreement with that obtained by inverting a traditional multi-step outflow experiment. Monitoring the initial phase of capillary rise reveals behaviour that deviates qualitatively from the single-phase, local-equilibrium regime described by Richards’ equation. Keywords: X-ray absorption, soil hydraulic properties, soil water dynamics, Richards’ equation


Sign in / Sign up

Export Citation Format

Share Document