Finite formulation for the computation of sorptivity

2020 ◽  
Author(s):  
Pierre-Emmanuel Peyneau ◽  
Laurent Lassabatere ◽  
Joseph Pollacco ◽  
Jesús Fernández-Gálvez ◽  
Borja Latorre ◽  
...  
Keyword(s):  
Hydraulic Conductivity ◽  
Water Content ◽  
Water Retention ◽  
Water Pressure ◽  
Pressure Head ◽  
Concentration Function ◽  
Point Of View ◽  
Water Infiltration ◽  
Unsaturated Hydraulic Conductivity ◽  
Flux Concentration

<p>Soil sorptivity is one of the key hydraulic parameters for modelling water infiltration into soil. It quantifies the capacity of a soil to infiltrate water by capillarity. Several formulations, based on various models, have been proposed to compute it from the water retention and the unsaturated hydraulic conductivity functions. All these formulations use the integration of the product of either the hydraulic conductivity or diffusivity function with the flux concentration function. The integration can be performed either over an interval of water pressure head or water content, yielding two equal values. However, the expression of the integral as a function of water pressure head may involve a huge or even infinite interval, which can be numerically difficult to handle. In opposite, the expression of the integral as a function of water content involves the integration of a diverging function (diffusivity) over a large interval, which is also troublesome from a numerical point of view. In this paper, we provide a new expression for sorptivity by cutting the integral in two parts, in order to involve only the integration of a finite function over a finite interval. The dependency of the integral on the flux concentration function is also investigated.</p>

scholarly journals Simultaneous Determination of Water Retention Curve and Unsaturated Hydraulic Conductivity of Substrates Using a Steady-state Laboratory Method

HortScience ◽  
2010 ◽  
Vol 45 (7) ◽  
pp. 1106-1112 ◽  
Author(s):  
Paraskevi A. Londra
Keyword(s):  
Experimental Data ◽  
Steady State ◽  
Hydraulic Conductivity ◽  
Water Retention ◽  
Water Pressure ◽  
Pressure Head ◽  
Laboratory Method ◽  
Water Retention Curve ◽  
Retention Curve ◽  
Unsaturated Hydraulic Conductivity

For effective irrigation and fertilization management, the knowledge of substrate hydraulic properties is essential. In this study, a steady-state laboratory method was used to determine simultaneously the water retention curve, θ(h), and unsaturated hydraulic conductivity as a function of volumetric water content, K(θ), and water pressure head, K(h), of five substrates used widely in horticulture. The substrates examined were pure peat, 75/25 peat/perlite, 50/50 peat/perlite, 50/50 coir/perlite, and pure perlite. The experimental retention curve results showed that in the case of peat and its mixtures with perlite, there is a hysteresis between drying and wetting branches of the retention curve. Whereas in the case of coir/perlite and perlite, the phenomenon of hysteresis was less pronounced. The increase of perlite proportion in the peat/perlite mixtures led to a decrease of total porosity and water-holding capacity and an increase of air space. Study of the K(θ) and K(h) experimental data showed that the hysteresis phenomenon of K(θ) was negligible compared with the K(h) data for all substrates examined. Within a narrow range of water pressure head (0 to –70 cm H2O) that occurs between two successive irrigations, a sharp decrease of the unsaturated hydraulic conductivity was observed. The comparison of the K(θ) experimental data between the peat-based substrate mixtures and the coir-based substrate mixture showed that for water contents lower than 0.40 m3·m−3, the hydraulic conductivity of the 50/50 coir/perlite mixture was greater. The comparison between experimental water retention curves and predictions using Brooks-Corey and van Genuchten models showed a high correlation (0.992 ≤ R2 ≤ 1) for both models for all substrates examined. On the other hand, in the case of unsaturated hydraulic conductivity, the comparison showed a relatively good correlation (0.951 ≤ R2 ≤ 0.981) for the van Genuchten-Mualem model for all substrates used except perlite and a significant deviation (0.436 ≤ R2 ≤ 0.872) for the Brooks-Corey model for all substrates used.

scholarly journals Examining the effect of pore size distribution and shape on flow through unsaturated peat using computed tomography

2009 ◽  
Vol 13 (10) ◽  
pp. 1993-2002 ◽  
Author(s):  
F. Rezanezhad ◽  
W. L. Quinton ◽  
J. S. Price ◽  
D. Elrick ◽  
T. R. Elliot ◽  
...  
Keyword(s):  
Computed Tomography ◽  
Hydraulic Conductivity ◽  
Pore Structure ◽  
Pore Size ◽  
Soil Water ◽  
Water Pressure ◽  
Peat Soil ◽  
Pressure Head ◽  
Peat Soils ◽  
Unsaturated Hydraulic Conductivity

Abstract. The hydraulic conductivity of unsaturated peat soil is controlled by the air-filled porosity, pore size and geometric distribution as well as other physical properties of peat materials. This study investigates how the size and shape of pores affects the flow of water through peat soils. In this study we used X-ray Computed Tomography (CT), at 45 μm resolution under 5 specific soil-water pressure head levels to provide 3-D, high-resolution images that were used to detect the inner pore structure of peat samples under a changing water regime. Pore structure and configuration were found to be irregular, which affected the rate of water transmission through peat soils. The 3-D analysis suggested that pore distribution is dominated by a single large pore-space. At low pressure head, this single large air-filled pore imparted a more effective flowpath compared to smaller pores. Smaller pores were disconnected and the flowpath was more tortuous than in the single large air-filled pore, and their contribution to flow was negligible when the single large pore was active. We quantify the pore structure of peat soil that affects the hydraulic conductivity in the unsaturated condition, and demonstrate the validity of our estimation of peat unsaturated hydraulic conductivity by making a comparison with a standard permeameter-based method. Estimates of unsaturated hydraulic conductivities were made for the purpose of testing the sensitivity of pore shape and geometry parameters on the hydraulic properties of peats and how to evaluate the structure of the peat and its affects on parameterization. We also studied the ability to quantify these factors for different soil moisture contents in order to define how the factors controlling the shape coefficient vary with changes in soil water pressure head. The relation between measured and estimated unsaturated hydraulic conductivity at various heads shows that rapid initial drainage, that changes the air-filled pore properties, creates a sharp decline in hydraulic conductivity. This is because the large pores readily lose water, the peat rapidly becomes less conductive and the flow path among pores, more tortuous.

scholarly journals Modified tension infiltrometer and its use to determine the unsaturated hydraulic conductivity of upper vadose zone

1998 ◽  
Vol 18 ◽  
Author(s):  
Golam Shabbir Sattar
Keyword(s):  
Steady State ◽  
Hydraulic Conductivity ◽  
Vadose Zone ◽  
Water Pressure ◽  
Fundamental Property ◽  
Pressure Head ◽  
Unknown Parameters ◽  
Steady State Condition ◽  
Tension Infiltrometer ◽  
Unsaturated Hydraulic Conductivity

A modification of the CSIRO type tension infiltrometer was designed (designated UNCEL type) and tested. This infiltrometer is considerably cheaper, but more versatile, than any other infiltrometer reported so far. The modified infiltrometer was used in a reconstructed vadose zone to measure the unsaturated hydraulic conductivity and effect of compaction associated with the upper vadose zone. Hydraulic conductivity, a function of matric potential, is a fundamental property for water entry into the soil as well as in the upper vadose zone. Tension infiltrometer technique for determining unsaturated hydraulic conductivity, which is also a function of soil water pressure head, was compared with the existing techniques. The analysis of data from tension infiltrometer is based mainly on flow rate at steady-state condition. The steady-state measurement also enables to determine unknown parameters on the basis of at least two negative heads at the same location. This modified infiltrometer imposes pressure potentials at soil surface from 0.01 to 0.15 m of water, as a result macropores with an air entry value of less than applied tension are excluded. This is practically useful when investigating the influence of structure and compaction of upper vadose zone. Measurement of infiltration at various negative pressure head, with a single disc diameter allowed sensitive measurement of hydraulic properties of upper vadose zone with minimal surface disturbance. The newly designed tension infiltrometer and adopted schemes of calculation enable to determine unsaturated hydraulic conductivity of upper vadose zone with relevance to structural effects more accurately.

Application of the Groenevelt - Grant soil water retention model to predict the hydraulic conductivity

Soil Research ◽  
2010 ◽  
Vol 48 (5) ◽  
pp. 447 ◽  
Author(s):  
C. D. Grant ◽  
P. H. Groenevelt ◽  
N. I. Robinson
Keyword(s):  
Hydraulic Conductivity ◽  
Water Content ◽  
Water Retention ◽  
Soil Water Retention ◽  
Retention Model ◽  
Unsaturated Hydraulic Conductivity ◽  
Incomplete Gamma Functions ◽  
Relative Water ◽  
Van Genuchten Model ◽  
Conductivity Function

We outline several formulations of the Groenevelt–Grant water retention model of 2004 to show how it can be anchored at different points. The model is highly flexible and easy to perform multiple differentiations and integrations on. Among many possible formulations of the model we choose one anchored solely at the saturated water content, θs, to facilitate comparison with the van Genuchten model of 1980 and to obtain a hydraulic conductivity function through analytical integration: where, k0, k1, and n are fitting parameters. We divided this formulation by θs to obtain the relative water content, θr(h), and inverted the function to produce a form required for integration, namely: in which the parameter β is introduced to accommodate both the ‘Burdine’ and ‘Mualem’ models. The integrals are identified as incomplete gamma functions and are distinctly different from the incomplete beta functions embodied in the van Genuchten–Mualem models. Rijtema’s data from 1969 for 20 Dutch soils are used to demonstrate the procedures involved. The water retention curves produced by our Groenevelt–Grant model are virtually indistinguishable from those produced by the van Genuchten model. Relative hydraulic conductivities produced by our Mualem and Burdine models produced closer estimates of Rijtema’s measured values than those produced by the van Genuchten–Mualem model for 19 of his 20 soils. This work provides an alternative to the widely used van Genuchten–Mualem approach and represents a preamble for the, as yet unsatisfactory, treatment of the tortuosity component of the unsaturated hydraulic conductivity function.

Comparison of the Instantaneous Profile Method and inverse modelling for the prediction of effective soil hydraulic properties

Soil Research ◽  
2005 ◽  
Vol 43 (5) ◽  
pp. 599 ◽  
Author(s):  
Oagile Dikinya
Keyword(s):  
Hydraulic Conductivity ◽  
Water Content ◽  
Sandy Soil ◽  
Hydraulic Properties ◽  
Water Retention ◽  
Pressure Head ◽  
Inverse Modelling ◽  
Soil Hydraulic Properties ◽  
Loamy Soil ◽  
Soil Hydraulic

Soil hydraulic conductivity K(θ) and soil water retention θ(h) have been determined from a drainage experiment. Two lysimeters, one filled with a sandy soil and the other with a loamy soil, were set up for a 1-dimensional transient flow experiment. The data were collected after flooding the lysimeters with water. Soil water contents were measured by time domain reflectrometry (TDR) and pressure heads were measured by tensiometers with mercury manometers. The experimental data determined by the instantaneous profile method (IPM) were compared with the results obtained by inverse modelling. The inverse modelling proved to be superior to the IPM methodology in effective prediction of hydraulic properties. The measurable properties water content and pressure head were optimised for the following datasets: water content (WC), pressure head (P-h), and a combination of WC and P-h. For both soils the optimisation of the dataset with both WC and P-h resulted in parameters that corresponded closely to the soil hydraulic data generated by the IPM method. The correspondence for the water retention data was better than for the hydraulic conductivity data. The datasets with WC only or P-h only did not contain enough information to accurately estimate the soil hydraulic properties. In most cases the results indicated that the sandy soil gave better agreement than the loamy soil. This was attributed to the faster drainage of the sandy than the loamy soil.

scholarly journals Examining the effect of pore size distribution and shape on flow through unsaturated peat using 3-D computed tomography

2009 ◽  
Vol 6 (3) ◽  
pp. 3835-3862 ◽  
Author(s):  
F. Rezanezhad ◽  
W. L. Quinton ◽  
J. S. Price ◽  
D. Elrick ◽  
T. R. Elliot ◽  
...  
Keyword(s):  
Computed Tomography ◽  
Hydraulic Conductivity ◽  
Pore Structure ◽  
Pore Size Distribution ◽  
Pore Size ◽  
Size Distribution ◽  
Water Pressure ◽  
Pressure Head ◽  
Peat Soils ◽  
Unsaturated Hydraulic Conductivity

Abstract. The hydraulic conductivity of unsaturated peat soils is controlled by the peat structure which affects the air-filled porosity, pore size distribution and shape. This study investigates how the size and shape of pores affects the flow of water through peat soils. In this study we used X-ray Computed Tomography (CT), at 45 µm resolution under 5 specific soil-water pressure head levels to provide 3-D, high-resolution images that were used to detect the inner pore structure of peat samples under a changing water regime. Pore structure and configuration were found to be irregular, which affected the rate of water transmission through peat soils. The 3-D analysis suggested that pore distribution is dominated by a single large pore-space. At low pressure head, this single large air-filled pore imparted a more effective flowpath compared to smaller pores. Smaller pores were disconnected and the flowpath was more tortuous than in the single large air-filled pore, and their contribution to flow was negligible when the single large pore was active. We quantify the pore structure of peat soil that affects the hydraulic conductivity in the unsaturated condition, and demonstrate the validity of our estimation of peat unsaturated hydraulic conductivity by making a comparison with a standard permeameter-based method. Estimates of unsaturated hydraulic conductivities were made for the purpose of testing the sensitivity of pore shape and geometry parameters on the hydraulic properties of peats and how to evaluate the structure of the peat and its affects on parameterization. We also studied the ability to quantify these factors for different soil moisture contents in order to define how the factors controlling the shape coefficient vary with changes in soil water pressure head. The relation between measured and estimated unsaturated hydraulic conductivity at various heads shows that rapid initial drainage, that changes the air-filled pore properties, creates a sharp decline in hydraulic conductivity. This is because the large pores readily lose water, the peat rapidly becomes less conductive and the flow path among pores, more tortuous.

scholarly journals Scaling procedure for straightforward computation of sorptivity

2021 ◽  
Author(s):  
Laurent Lassabatere ◽  
Pierre-Emmanuel Peyneau ◽  
Deniz Yilmaz ◽  
Joseph Pollacco ◽  
Jesús Fernández-Gálvez ◽  
...  
Keyword(s):  
Hydraulic Conductivity ◽  
Water Retention ◽  
Unsaturated Flow ◽  
Water Pressure ◽  
Water Infiltration ◽  
Integral Parameter ◽  
Residual Water ◽  
Shape Parameters ◽  
Scaling Procedure ◽  
Soil Hydraulic

Abstract. Sorptivity is a parameter of primary importance in the study of unsaturated flow in soils. This integral parameter is often considered for modeling the computation of water infiltration into vertical soil profiles (1D or 3D axisymmetric geometry). Sorptivity can be directly estimated from the knowledge of the soil hydraulic functions (water retention of hydraulic conductivity), using the integral formulation of Parlange (Parlange, 1975). However, it requires the prior determination of the soil hydraulic diffusivity and its numerical integration between the initial and the final saturation degrees, which may be tricky for some instances (e.g., coarse soil with diffusivity functions quasi-infinite close to saturation). In this paper, we present a specific scaling procedure for the computation of sorptivity considering slightly positive water pressure heads at the soil surface and initial dry conditions (corresponding to most water infiltration on the field). The square sorptivity is related to the square dimensionless sorptivity (referred to as cp parameter) corresponding to a unit soil (i.e., unit values of all the scaled parameters and zero residual water content) utterly dry at the initial state and saturated at the final state. The cp parameter was computed numerically and analytically for five current models: delta functions (Green and Ampt model), Brooks and Corey, van Genuchten-Mualem, van Genuchten-Burdine, and Kosugi models as a function of the shape parameters. The values are tabulated and can be easily used to determine any dimensional sorptivity value for any case. We propose brand-new analytical expressions for some of the models and validate previous formulations for the other models. Our numerical results also showed that the relation between the cp  parameters and shape parameters strongly depends on the chosen model, with either increasing or decreasing trends when moving from coarse to fine soils. These results highlight the need for carefully selecting the proper model for the description of the water retention and hydraulic conductivity functions for the rigorous estimation of sorptivity. Present results show the need to understand better the hydraulic model's mathematical properties, including the links between their parameters, and, secondly, to better relate these properties to the physical processes of water infiltration into soils.

Hyphal colonization of Rhizophagus irregularis increases unsaturated hydraulic conductivity of a loamy sand distant from roots

2021 ◽  
Author(s):  
Michael Bitterlich ◽  
Richard Pauwels
Keyword(s):  
Hydraulic Conductivity ◽  
Water Content ◽  
Hydraulic Properties ◽  
Volumetric Water Content ◽  
Rhizophagus Irregularis ◽  
Loamy Sand ◽  
Soil Hydraulic Properties ◽  
Unsaturated Hydraulic Conductivity ◽  
Soil Hydraulic ◽  
Root Exclusion

<p>Hydraulic properties of mycorrhizal soils have rarely been reported and difficulties in directly assigning potential effects to hyphae of arbuscular mycorrhizal fungi (AMF) arise from other consequences of AMF being present, i.e. their influence on growth and water consumption rates of their host plants that both also influence soil hydraulic properties.</p><p>We assumed that the typical nylon meshes used for root-exclusion experiments in mycorrhizal research can provide a dynamic hydraulic barrier. It is expected that the uniform pore size of the rigid meshes causes a sudden hydraulic decoupling of the enmeshed inner volume from the surrounding soil as soon as the mesh pores become air-filled. Growing plants below the soil moisture threshold for hydraulic decoupling would minimize plant-size effects on root-exclusion compartments and allow for a more direct assignment of hyphal presence to modulations in soil hydraulic properties.</p><p>We carried out water retention and hydraulic conductivity measurements with two tensiometers introduced in two different heights in a cylindrical compartment (250 cm³) containing a loamy sand, either with or without the introduction of a 20 µm nylon mesh equidistantly between the tensiometers. Introduction of a mesh reduced hydraulic conductivity across the soil volumes by two orders of magnitude from 471 to 6 µm d<sup>-1</sup> at 20% volumetric water content.</p><p>We grew maize plants inoculated or not with Rhizophagus irregularis in the same soil in pots that contained root-exclusion compartments while maintaining 20% volumetric water content. When hyphae were present in the compartments, water potential and unsaturated hydraulic conductivity increased for a given water content compared to compartments free of hyphae. These differences increased with progressive soil drying.</p><p>We conclude that water extractability from soils distant to roots can be facilitated under dry conditions when AMF hyphae are present.</p><p> </p>

Experimental study on dual capillary barrier using recycled asphalt pavement materials

2014 ◽  
Vol 51 (10) ◽  
pp. 1165-1177 ◽  
Author(s):  
F.R. Harnas ◽  
H. Rahardjo ◽  
E.C. Leong ◽  
J.Y. Wang
Keyword(s):  
Water Content ◽  
Pore Water Pressure ◽  
Water Pressure ◽  
Water Storage ◽  
Pressure Head ◽  
Volumetric Water Content ◽  
Capillary Barrier ◽  
Recycled Asphalt ◽  
Fine Grained ◽  
Drainage Rate

The performance of a capillary barrier cover as a cover system is affected by the ability of the capillary barrier to store water. To increase the water storage of a capillary barrier cover, the dual capillary barrier (DCB) concept is proposed. The objective of this paper is to investigate the water storage of the proposed DCB as compared to the storage of a traditional single capillary barrier (SCB). The investigation is conducted using two one-dimensional infiltration column tests under different rainfall conditions. The results show that a DCB stores more water as compared to SCB. The results show that the fine-grained layers of a DCB have higher volumetric water contents during drainage as compared to that of the fine-grained layer of an SCB. The higher volumetric water content is caused by the fact that the thickness of the layers in a DCB corresponds to a pore-water pressure head range where the material has the highest volumetric water content. In addition, a slower drainage rate is resulted from additional layering in a DCB.

Sign in / Sign up

Export Citation Format

Share Document