Acessibilidade / Reportar erro

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

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.

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)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.03615995004400050002x
https://doi.org/10.2136/sssaj1980.036159...
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)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/WR012i003p00513
https://doi.org/10.1029/WR012i003p00513...
or Burdine (1953)Burdine NT. Relative permeability calculation from pore-size distribution data. J Petrol Technol. 1953;5:71-8. https://doi.org/10.2118/225-G
https://doi.org/10.2118/225-G...
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., 2014Lima 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.v9i4a3532
https://doi.org/10.5039/agraria.v9i4a353...
; 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., 2009Silva 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-06832009000600003
https://doi.org/10.1590/S0100-0683200900...
; Mota et al., 2017Mota 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.4023
https://doi.org/10.18227/1982-8470ragro....
). In international literature, a similar description can sometimes be found (Pollaco and Mohanty, 2012Pollaco 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.0167
https://doi.org/10.2136/vzj2011.0167...
; Aschonitis and Antonopoulos, 2013Aschonitis 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.1751
https://doi.org/10.1002/ird.1751...
; Aschonitis et al., 2015Aschonitis 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-3
https://doi.org/10.1007/s10706-015-9887-...
; Dokoohaki et al., 2017Dokoohaki 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.1358742
https://doi.org/10.1080/00103624.2017.13...
).

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:

h b = 2 σcosφ ρgr m , 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):

Θ = h h b - λ for h < h b Θ = 1 for h h b 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)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.03615995004400050002x
https://doi.org/10.2136/sssaj1980.036159...
water retention equation (VG, Equation 3):

Θ = 1 + α h n - 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)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/WR012i003p00513
https://doi.org/10.1029/WR012i003p00513...
or Burdine (1953)Burdine NT. Relative permeability calculation from pore-size distribution data. J Petrol Technol. 1953;5:71-8. https://doi.org/10.2118/225-G
https://doi.org/10.2118/225-G...
, 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 - 1 n and n > 1 Eq. 4

while the Burdine restriction is (Equation 5):

m = 1 - 2 n 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):

h b = 1 α α = 1 h b 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:

Θ = α h n - 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Š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.0033
https://doi.org/10.2136/vzj2016.04.0033...
), 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.

Figure 1
Representative values of |hb| for soils from the Hydrus package (Šimůnek et al., 2016Š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.0033
https://doi.org/10.2136/vzj2016.04.0033...
).

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 = 1 2 m 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.

Figure 2
Effective saturation Θb as a function of the Van Genuchten parameter m, assuming parameter a to be the inverse of the air-entry pressure.

Figure3
Effective saturation Θ as a function of the matric potential for two values of parameter a and n = 2 (top) and n = 5 (bottom). The lines 1/ a indicate the supposed values of the air-entry pressure according to equation 6, with corresponding Θb given by equation 8.

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.1751
    » https://doi.org/10.1002/ird.1751
  • 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-3
    » https://doi.org/10.1007/s10706-015-9887-3
  • 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-g
    » http://dx.doi.org/10.2118/225-g
  • Burdine NT. Relative permeability calculation from pore-size distribution data. J Petrol Technol. 1953;5:71-8. https://doi.org/10.2118/225-G
    » https://doi.org/10.2118/225-G
  • 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.1358742
    » https://doi.org/10.1080/00103624.2017.1358742
  • 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.v9i4a3532
    » https://doi.org/10.5039/agraria.v9i4a3532
  • 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.4023
    » https://doi.org/10.18227/1982-8470ragro.v11i2.4023
  • 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/WR012i003p00513
    » https://doi.org/10.1029/WR012i003p00513
  • 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-43662014000600010
    » https://doi.org/10.1590/S1415-43662014000600010
  • 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.0167
    » https://doi.org/10.2136/vzj2011.0167
  • 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-06832009000600003
    » https://doi.org/10.1590/S0100-06832009000600003
  • Š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.0033
    » https://doi.org/10.2136/vzj2016.04.0033
  • 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-06832008000500001
    » https://doi.org/10.1590/S0100-06832008000500001
  • 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-43662008000200004
    » https://doi.org/10.1590/S1415-43662008000200004
  • 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.03615995004400050002x
    » https://doi.org/10.2136/sssaj1980.03615995004400050002x

Publication Dates

  • Publication in this collection
    02 July 2018
  • Date of issue
    2018

History

  • Received
    30 Oct 2017
  • Accepted
    23 Mar 2018
Sociedade Brasileira de Ciência do Solo Secretaria Executiva , Caixa Postal 231, 36570-000 Viçosa MG Brasil, Tel.: (55 31) 3899 2471 - Viçosa - MG - Brazil
E-mail: sbcs@ufv.br