## Print version ISSN 0100-0683On-line version ISSN 1806-9657

### Rev. Bras. Ciênc. Solo vol.42  Viçosa  2018  Epub July 02, 2018

#### http://dx.doi.org/10.1590/18069657rbcs20170343

Division - Soil Processes and Properties

Commission - Soil Physics

An Alert Regarding a Common Misinterpretation of the Van Genuchten α Parameter

(1)Universidade de São Paulo, Centro de Energia Nuclear na Agricultura, Piracicaba, São Paulo, Brasil.

ABSTRACT

Among the equations available to describe the relation between matric potential and soil water content, the soil water retention function, the most commonly used is the equation proposed by Van Genuchten in his 1980 landmark paper. In soil physics literature, especially in Brazil, several authors relate the inverse of the Van Genuchten parameter α to the air-entry pressure. This study aimed to show this common interpretation to be erroneous, as 1/α corresponds to water contents lower than saturation. The deviation depends on the m parameter. In fact, α is merely a scaling parameter relative to the matric potential axis. Recognizing this mathematical fact may improve the interpretation of soil hydraulic properties based on water retention parameters.

Key words: soil water retention; soil physics; empirical equations

INTRODUCTION

Several equations are available to describe the relation between matric potential and soil water content, the soil water retention function. The most frequently applied ones are the Brooks and Corey (1964) equation (BC) and the van Genuchten (1980) equation (VG). When expressing effective saturation Θ, a quantity scaling water content from 0 to 1 between residual and saturated values, the BC equation contains two parameters, one of which (hb) explicitly represents the matric potential corresponding to air-entry. On the other hand, the VG equation, to be applied using the Mualem (1976) or Burdine (1953) parametric restriction, holds also two parameters (α and n), but none of them has a clear physical meaning.

In Brazilian soil physics literature, several authors relate the inverse of parameter α to the air-entry pressure, in analogy to the BC parameter hb. Some use this kind of identification when describing VG parameters (Souza et al., 2008a; Lima et al., 2014; Oliveira Júnior et al., 2014), others when explicitly interpreting results of α values to the air-entry pressure (Souza et al., 2008b; Silva et al., 2009; Mota et al., 2017). In international literature, a similar description can sometimes be found (Pollaco and Mohanty, 2012; Aschonitis and Antonopoulos, 2013; Aschonitis et al., 2015; Dokoohaki et al., 2017).

Here, we demonstrate that this interpretation of the VG α parameter is incorrect and should therefore be avoided.

DEVELOPMENT

The air-entry pressure, or “bubbling pressure” hb (m), of a soil or porous material is defined as the matric potential at which the first (largest) pore starts draining its water (Brooks and Corey, 1964). Considering the Young-Laplace capillary equation (Equation 1), it is determined by the radius of the largest pore rm (m) as:

hb=2σcosφρgrm, Eq. 1

in which σ (J m-2) is the surface tension of water, φ the contact angle between the water surface, the surrounding air, and the pore walls, r (kg m-3) the density of water, and g (m s-2) the gravity. In some water retention models, hb (m) is an explicit fitting parameter, notably in the Brooks and Corey (1964) model (Equation 2):

Θ=hhb-λ for h < hb Θ= 1 for hhb Eq. 2

in which h is the matric potential, Θ = (θ – θr)/(θs – θr) is the effective saturation, θ, θr, and θs are water content, residual water content, and saturated water content, respectively, all on a volume base (m3 m-3). The air-entry pressure corresponds to the onset of water content reduction with further decreasing matric potentials. As such, the water content at the air-entry pressure θb (m3 m-3) equals the saturated water content θs and Θ = 1 at h = hb (Equation 2).

The air-entry pressure is not explicitly present in the frequently used van Genuchten (1980) water retention equation (VG, Equation 3):

Θ=1+αhn-m Eq. 3

in which a, n, and m (function of n) are fitting parameters, a having the inverse dimension of h (e.g. m-1). The VG equation is defined together with the theory presented by Mualem (1976) or Burdine (1953), and when applying the respective parametric restrictions (defining m as a function of n), it can be used to estimate the hydraulic conductivity function from retention parameters. The Mualem restriction is as follows (Equation 4):

m=1-1n and n>1 Eq. 4

while the Burdine restriction is (Equation 5):

m=1-2n and n>2 Eq. 5

Consequently, 0< m <1. As mentioned, many authors assume that α is the inverse of the absolute value of the air-entry pressure hb, i.e. (Equation 6):

hb=1αα=1hb Eq. 6

This assumption has its origin in a comparison between equation 2 and 3. If |h| becomes very large, equation 3 reduces to equation 7:

Θ=αhn-m Eq. 7

Equation 7 is equal to equation 2, with λ = mn (or, l = n - 1) with the Mualem restriction, equation 4, and λ = n - 2 with the Burdine restriction, equation 5 and a given by equation 6. However, this does not justify the interpretation of a as the inverse of the bubbling pressure, as equation 7 is only valid for very large values of |h|, whereas |hb| is, in fact, a relatively small value. Fitting soils with several textures from the Hydrus package (Šimůnek et al., 2016), values of |hb| range between 0.05 and 0.4 m (Figure 1), corresponding to pore diameters of 5.87 ∙ 10-4 and 7.35 ∙ 10-5 m, respectively.

Moreover, if the interpretation of a as the inverse of the bubbling pressure were true, then combining equation 6 to equation 3 would result in the following expression (Equation 8) for the effective saturation, corresponding to the bubbling pressure Θb:

Θb=12m Eq. 8

Equation 8 yields values for Θb between 1 (at m = 0) and 0.5 (at m = 1), as shown in figure 2, and in obvious disagreement with the notion that θb = θs; consequently, Θb = 1. The value of m = 0 implies in n = 1 (Mualem restriction, equation 4) or n = 2 (Burdine restriction, equation 5). Such values are physically unrealistic, as m = 0 results in Θ = 1 for any value of h (Equation 7). Equation 8 implies in the fact that the greater the value of m (and, from equations 4 and 5, the greater the value of n), the larger the deviation between the inverse of a and the air-entry pressure. This can also be seen in figure 3, showing the retention curves (Θ as a function of the matric potential) for two values of parameter a and n = 2 (top) and n = 5 (bottom). In this figure, 1/a indicates the supposed values of the air-entry pressure according to equation 6, with corresponding Θb given by equation 8. Figure 2 also clearly demonstrates the effect of a on the shape of the retention curve, with a being a mere scaling parameter relative to the matric potential axis.

CONCLUSION

We showed, mathematically and graphically, that the Van Genuchten retention equation parameter a is not equal to, nor simply correlated to the (inverse of) air-entry matric potential, as frequently alleged. Instead, a is a scaling parameter relative to the matric potential axis. Recognizing this mathematical fact may improve the interpretation of soil hydraulic properties based on water retention parameters and prevent the error of using the relationship shown in equation 6 to correlate parameters from equations 2 and 3.

REFERENCES

Aschonitis VG, Antonopoulos VZ. New equations for the determination of soil saturated hydraulic conductivity using the van Genuchten model parameters and effective porosity. Irrig Drain. 2013;62:537-42. https://doi.org/10.1002/ird.1751Links ]

Aschonitis VG, Salemi E, Colombani N, Mastrocicco M. Comparison of different ‘‘S-index’’ expressions to evaluate the state of physical soil properties. Geotech Geol Eng. 2015;33:1055-66. https://doi.org/10.1007/s10706-015-9887-3Links ]

Brooks RH, Corey AT. Hydraulic properties of porous media: Hydrology Papers. Fort Collins: Colorado State University; 1964. Burdine NT. Relative permeability calculation from pore-size distribution data. J Petrol Technol. 1953;198:71-7. http://dx.doi.org/10.2118/225-gLinks ]

Burdine NT. Relative permeability calculation from pore-size distribution data. J Petrol Technol. 1953;5:71-8. https://doi.org/10.2118/225-GLinks ]

Dokoohaki H, Miguez FE, Laird D, Horton R, Basso AS. Assessing the biochar effects on selected physical properties of a sandy soil: an analytical approach. Commun Soil Sci Plan. 2017;48:1387-98. https://doi.org/10.1080/00103624.2017.1358742Links ]

Lima JRS, Souza ES, Antonino ACD, Silva IF, Corrêa MM, Lira CAB. Atributos físico-hídricos de um Latossolo Amarelo cultivado e sob mata nativa no Brejo Paraibano. Rev Bras Cienc Agrar. 2014;9:599-605. https://doi.org/10.5039/agraria.v9i4a3532Links ]

Mota JCA, Libardi PL, Brito AS, Moraes SO, Nascimento IV, Alencar TL. Variabilidade espacial dos parâmetros da equação de van Genuchten em um Latossolo Vermelho-Amarelo. Revista Agro@mbiente On-line. 2017;11:92-100. https://doi.org/10.18227/1982-8470ragro.v11i2.4023Links ]

Mualem Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour Res. 1976;12:513-22. https://doi.org/10.1029/WR012i003p00513Links ]

Oliveira Júnior JAS, Souza ES, Correa MM, Lima JRS, Souza RMS, Silva Filho LA. Variabilidade espacial de propriedades hidrodinâmicas de um Neossolo Regolítico sob pastagem e caatinga. Rev Bras Eng Agric Ambient. 2014;18:631-9. https://doi.org/10.1590/S1415-43662014000600010Links ]

Pollaco JAP, Mohanty BP. Uncertainties of water fluxes in soil-vegetation-atmosphere transfer models: inverting surface soil moisture and evapotranspiration retrieved from remote sensing. Vadose Zone J. 2012;11:vzj2011.0167. https://doi.org/10.2136/vzj2011.0167Links ]

Silva AP, Leão TP, Tormena CA, Gonçalves ACA. Determinação da permeabilidade ao ar em amostras indeformadas de solo pelo método da pressão decrescente. Rev Bras Cienc Solo. 2009;33:1535-45. https://doi.org/10.1590/S0100-06832009000600003Links ]

Šimůnek J, van Genuchten MT, Šejna M. Recent developments and applications of the HYDRUS computer software packages. Vadose Zone J. 2016;15:vzj2016.04.0033. https://doi.org/10.2136/vzj2016.04.0033Links ]

Souza ES, Antonino ACD, Ângulo-Jaramillo R, Maciel Netto A, Montenegro SMGL, Silva EB. Variabilidade espacial dos parâmetros hidrodinâmicos de duas parcelas agrícolas no estado da Paraíba. Rev Bras Cienc Solo. 2008b;32:1795-804. https://doi.org/10.1590/S0100-06832008000500001Links ]

Souza ES, Antonino ACD, Angulo-Jaramillo R, Maciel Netto A. Caracterização hidrodinâmica de solos: aplicação do método Beerkan. Rev Bras Eng Agric Ambient. 2008a;12:128-35. https://doi.org/10.1590/S1415-43662008000200004Links ]

van Genuchten MT. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci Soc Am J. 1980;44:892-8. https://doi.org/10.2136/sssaj1980.03615995004400050002xLinks ]

Received: October 30, 2017; Accepted: March 23, 2018

* Corresponding author: E-mail: qdjvlier@usp.br

This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.