Next Article in Journal
Stormwater Runoff Treatment Using Rain Garden: Performance Monitoring and Development of Deep Learning-Based Water Quality Prediction Models
Next Article in Special Issue
The Effect of Soil Texture on the Conversion Factor of 1:5 Soil/Water Extract Electrical Conductivity (EC1:5) to Soil Saturated Paste Extract Electrical Conductivity (ECe)
Previous Article in Journal
Using Freshwater Heads to Analyze Flow Directions in Saline Aquifers of the Pingtung Plain, Taiwan
Previous Article in Special Issue
Effects of Drying-Rewetting Cycles on Ferrous Iron-Involved Denitrification in Paddy Soils
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Infiltration under Ponded Conditions

Laboratory of Agricultural Hydraulics, Department of Natural Resources Management and Agricultural Engineering, School of Environment & Agricultural Engineering, Agricultural University of Athens, 75 Iera Odos Street, 11855 Athens, Greece
*
Authors to whom correspondence should be addressed.
Water 2021, 13(24), 3492; https://doi.org/10.3390/w13243492
Submission received: 6 October 2021 / Revised: 20 November 2021 / Accepted: 3 December 2021 / Published: 7 December 2021
(This article belongs to the Special Issue Study of the Soil Water Movement in Irrigated Agriculture Ⅱ)

Abstract

:
Ponded infiltration processes occur in agricultural lands irrigated by flooding of their soil surface or under insufficient drainage conditions. The existing equations describing the phenomenon of vertical infiltration under ponded conditions have not considered the actual contribution of the pressure head gradient to the flow. In this study, simple equations are proposed to describe the horizontal and vertical infiltration under various ponding heads incorporating the actual contribution of the pressure head gradient to the flow. Six soils with known hydraulic properties, covering a wide range of soil textures, were used. Horizontal and vertical infiltration data are obtained by numerical simulation for all soils studied using the Hydrus-1D code. To validate the accuracy of the proposed equations, the solutions of horizontal and vertical infiltrations provided by the proposed equations were compared with numerically simulated ones provided by the Hydrus 1-D. The analysis of the results showed a very good agreement in all soils studied. The proposed vertical infiltration equation was also compared to a simple and accurate equation which does not incorporate the actual contribution of the pressure head gradient to the flow and differences between them were observed in all soils studied.

1. Introduction

The infiltration process is of great importance in hydrology and agricultural sciences since it provides the water available for plants and groundwater recharge and defines water runoff at soil surface. A rainfall or irrigation intensity greater than soil infiltration capacity will lead to water runoff at the soil surface, causing ponded conditions. Additionally, in agricultural lands ponded conditions may be developed under insufficient drainage and in irrigation practices when irrigation water is applied by flooding the soil surface. Further ponded conditions are met in lakes, natural or artificial. Consequently, the study of the vertical infiltration under ponding heads is of great interest.
For this, Philip [1,2] tackled the problem of the vertical infiltration under ponded conditions, and he presented analytical solutions relative to this flow problem. Other researchers, later, investigated the same infiltration case and proposed equations to estimate the cumulative ponded infiltration in a homogeneous soil [3,4,5,6,7,8,9,10,11,12,13,14]. The main difference among these equations is the number of physical or fitting parameters used. In several models, the common physical parameters used are the soil sorptivity, S, and the saturated hydraulic conductivity, Ks, which are often met as the main parameters of the two-parameter models.
Green and Ampt [3] proposed a two-parameter equation assuming a piston-type water content profile with a well-defined wetting front characterized by a constant pressure head and constant pressure head at the surface. Brutsaert [4] proposed a three-parameter equation based on the quasi-analytical time-series infiltration solution of Philip [5] which can be converted into a two-parameter form by fixing the value of the third parameter. Parlange et al. [6] proposed a three-parameter equation based on integration of the water content-based form of Richards’ equation which can be converted into a two-parameter form by fixing the value of the third parameter. Parlange et al. [7] proposed a modified version of Parlange et al. [6] equation with five parameters including ponded conditions. Haverkamp et al. [8] proposed a very accurate six-parameter equation based on that of Parlange et al. [7] which can be converted into a three-parameter one. Swartzendruber [9] proposed a three-parameter equation which is a semi-analytical solution of the Richards’ equation for all times t > 0, which is similar to the Philip’s moderate time series solution. This equation can be converted into the two-parameter equation of Stroosnijder [10] by fixing the third parameter. Later, Valiantzas [11] proposed a simple two-parameter linearized equation which is a specific solution located approximately at the middle of the two limiting behavior soils defined by the Green and Ampt [3] and Talsma and Parlange [12] equations. Valiantzas compared his equation with other common nonlinear ones with two and three parameters [3,4,5,6,7,8,9,10]. The comparison showed that the three-parameter equations may reproduce unrealistic estimates and have convergence or nonuniqueness difficulties of the solution [11]. The equation proposed by Valiantzas, over the other two-parameter equations studied, has the advantage of simplicity and accuracy of sorptivity and saturated hydraulic conductivity estimations, as well as the ability to detect and eliminate abnormalities of the infiltration process [11].
However, all the above-mentioned researchers did not incorporate into their analyses the actual contribution of the pressure head gradient to the flow.
Recently, Poulovassilis and Argyrokastritis [15] determined the exact contribution of the pressure head gradient term in the vertical infiltration process in the case of homogeneous porous media under zero ponding head. Specifically, Poulovassilis and Argyrokastritis [15] showed that the contribution to the flow of the pressure head gradient St1/2, included in Philip’s [16] two term equation, is only applied in the case of horizontal infiltration. They defined analytically that this contribution is in the vertical case is (St1/2αt) and that it is smaller than the horizontal one by a factor αt, being a measurable function of the infiltration time t. In conclusion, they proposed a new two-term analytical equation for vertical infiltration which reproduces exactly a vertical cumulative curve and satisfies the physical rules of the infiltration process.
Up to now, no research has been carried out on the actual contribution of the pressure head gradient to flow during the vertical infiltration process under ponded conditions.
The purpose of the present study is to define the physical characteristics of infiltration under ponded conditions and to develop simple equations for it, conforming to these characteristics, incorporating the actual contribution of the pressure head gradient.

2. Materials and Methods

2.1. Theory

The infiltration process along the horizontal axis x in a homogeneous semi-infinite porous body is governed by the following initial conditions
θ = θ i h = h i t = 0 x 0 ,
θ = θ 0 h = h p t 0 x = 0 ,
where θ0 is the volumetric water content established at x = 0 at a pressure head hp ≥ 0 and the pair θi < θ0, hi < 0 denote the initial uniform steady state prevailing all along the porous column.
When hp = 0, the horizontal cumulative infiltration is described by the following equation
i x = K 0 0 t λ x = 0   d t = S t 1 2 ,
where K0 is the hydraulic conductivity when θ = θ0 prevailing for all t at x = 0 and the pressure head gradient λ x = 0 at x = 0 [5,17] is
λ x = 0 = ( d h d x ) x = 0 = 1 2 t 1 2 S K 0 ,
The imposition of the ponding head hp > 0 at x = 0, adds the pressure head gradient and results in establishing a zone along which positive pressure heads prevail. The length of this zone is continuously increasing as the infiltration time t increases. Thus, if x0 is the value of x at which the hydrostatic pressure becomes zero at time t, then x 0 0 when t 0 and x 0 as t , while the pressure head gradient ( d h p d x ) λ p 0 as t for all x and therefore h p 0 as x 0 . In this case
i x , p = S t 1 2 + K 0 0 t λ x = 0 , p d t ,
while for all x < x0 the hydraulic conductivity is K0.
The pressure head gradient of the first term at the right-hand side in Equation (5) tends to zero as t and the same is true for the gradient λ x = 0 , p . Therefore, the two terms must tend to become parallel to each other. In such a case, the integral term may be expressed in terms of S t 1 2 so that
K 0 0 t λ x = 0 , p d t = n x S t 1 2 ,
where nx > 0. Thus, we may write,
i x , p = S t 1 2 + n x S t 1 2 = S ( 1 + n x ) t 1 2
It may be mentioned that the imposition of hp affects the shape of the wetting front reducing it to a steeper form. From Equation (7)
i x , p t 1 2 = S ( 1 + n x )
so that the value of the ratio at the left hand of Equation (8) is increasing as nx increases and further that this ratio is a single-valued function of the independent variable nx.
If in conditions described by Equation (3) x is replaced by z then they govern, too, the vertical infiltration process. The gravitation component adds to the pressure head gradient, causes an increase of the infiltration rate and promotes the expansion of the zone with positive hydrostatic pressure. Thus, if z0 is the depth at which the hydrostatic pressure becomes zero at a time t, then z 0 0 as t 0 and z 0 as t , and therefore
( d z 0 d t ) > ( d x 0 d t )   and   z 0 > x 0 ,   ( d K 0 d z ) > ( d K 0 d x )  
for all finite t.
Poulovassilis and Argyrokastritis [15] determined that the contribution of the pressure head gradient to the vertical flow is smaller than that of the horizontal flow by a factor at, which is a variable depending on the infiltration time, but measurable if the soil properties S and K0 are known (Equation no. 20 from Poulovassilis and Argyrokastritis [15]). Thus, we may write,
i z , p = S ( 1 + n z ) t 1 2 a t + K 0 t
from which
i z , p + a t K 0 t t 1 2 = S ( 1 + n z )  
The ratio at the left-hand side of Equation (11) is a single-valued function of the independent variable nz.
Taking into account the analytic form of the variable at (Equation no. 22 from Poulovassilis and Argyrokastritis [15])
a t = S t 1 2 ( 1 e c ( K 0 S ) t ) ,
the Equation (10) can be written as
i z , p = S t 1 2 [ ( 1 + n z ) ( 1 e c ( K 0 S ) t ) ] + K 0 t ,
or
i z , p = S t 1 2 ( n z + e c ( K 0 S ) t ) + K 0 t ,
in which c is a constant fitting parameter characteristic of the porous medium.
Poulovassilis and Argyrokastritis [15] also demonstrated that the variable at may be described by the following equation for small and intermediate infiltration time t
α t = k t ,
where the slope coefficient k is characteristic of the porous medium.
In that case, Equation (2) may be approximated by the Poulovassilis and Argyrokastritis (Equation no. 25 from Poulovassilis and Argyrokastritis [15]) following equation
i z = S t 1 2 + ( K 0 k ) t ,

2.2. Porous Media

Six soils ranged from sand to clay with known hydraulic properties were selected covering a wide range of soil texture. The soils studied were the Yolo light clay (YLC), a clay loam (CL), a silty loam (SiL), a sandy loam (SL), a loamy sand (LS) and a sand (S). For all soils, water retention and hydraulic conductivity curves were described by the Mualem [18]-van Genuchten [19] closed-form model. For the Yolo light clay soil, the soil parameters were obtained from Poulovassilis et al. [20] and for the other studied soils by the database included in the Hydrus-1D software package [21]. The corresponding soil parameters are presented in Table 1.

2.3. Numerical Analysis

Very few experimental cumulative infiltration data are available for hp > 0, which is, perhaps, due to the additional experimental difficulties for obtaining such data. For this reason, infiltration data were obtained by numerical simulation for all porous media studied using the Hydrus-1D software package (Department of Environmental Sciences, University of California Riverside: Riverside, CA, USA) [21].
The Hydrus-1D is one of the most common software packages that has been successfully applied in various studies worldwide for predicting soil water movement under different conditions. The HYDRUS-1D is a finite element model for simulating the one-dimensional movement of water, heat, and multiple solutes in variably saturated media. The program numerically solves the Richards equation for saturated–unsaturated water flow and Fickian-based advection–dispersion equations for both heat and solute transport [21]. The unsaturated soil hydraulic properties can be described using van Mualem-Genuchten [18,19], Brooks and Corey [22], modified van Genuchten [23], lognormal distribution model of Kosugi [24], and a dual-porosity model [25] type analytical functions. In this study, the Mualem-van Genuchten closed-form model was used.
To obtain horizontal and vertical infiltration numerical data, boundary conditions were set. Upper boundary was defined as a constant soil surface ponding depth, hp, of 0, 5, 10, 25, 50, 100, 150 and 200 cm. As lower boundary of the uniform soil column, a zero-pressure head gradient was defined (free drainage). Initial pressure heads were set to −208, −500, −300, −200, −100 and −100 cm for YLC, CL, SiL, SL, LS and S, respectively, according to soil texture.
The value of sorptivity S for each porous medium was obtained using horizontal infiltration numerical data and applying the Philip’s [5] equation ix,n = St1/2 (Table 1).

2.4. Tests of the Formulated Equations

To validate the proposed Equations (7) and (14), the solutions of the ix,p and iz,p provided by these equations, respectively, were compared with numerically simulated ones, ix,p,n and iz,p,n, provided by the Hydrus-1D program for horizontal and vertical infiltration, respectively.
In order to compare the proposed Equation (14) with those proposed for vertical infiltration by other researchers, the Valiantzas equation (denoted as Val) was selected and its algebraic form is given by Equation (17). It is a simple and accurate two-parameter equation, among others more sophisticated needing three or more parameters, as presented in the introduction section.
i = 0.5 K s t + S t 1 2 [ 1 + ( 0.5 K s S ) 2 t ] 0.5 ,
The adjustment of soil sorptivity for various soil surface ponding depths was estimated by the following Equation [11]
S p = S 0 2 + 2 K s h p ( θ s θ i ) ,
where Sp is the soil sorptivity for ponding depth hp, S0 for hp = 0, θs and θi are the water content at saturation and the initial soil water content, respectively.

3. Results and Discussion

In Figure 1, the horizontal cumulative infiltration curves ix,p,n obtained numerically for all soils examined are compared to the predicted ones ix,p using the Equation (7), for various hp values. Correspondingly, the comparison between the vertical cumulative infiltration curves iz,p,n obtained numerically and the predicted ones iz,p using the Equation (14), for various hp values is depicted in Figure 2. The fitting parameter c in Equations (12)–(14) and the slope coefficient k in Equations (15) and (16) for each soil are presented in Table 1.
As shown in Figure 1 and Figure 2, the ix,p and iz,p predicted by applying Equations (7) and (14), for the horizontal and vertical infiltration, respectively, are in close agreement with those produced numerically by Hydrus-1D. Furthermore, the proposed equations (Equations (7) and (14)) were tested for their accuracy with the numerical data using the indices of Root Mean Square Error (RMSE) and Index of Agreement (d) [26]:
d = 1 j = 1 N ( i n , j i p r e d , j ) 2 j = 1 N ( | i p r e d , j i n ¯ | + | i n , j i n ¯ | ) 2 ,   0 d 1 ,
R M S E = j = 1 N ( i n , j i p r e d , j ) 2 N ,
where in,j and ipred,j are the numerical cumulative infiltration values and the predicted ones by the proposed Equation (7) or (14), respectively, i n ¯ is the mean of numerical values and N is the number of values.
As shown in Table 2, both the small values of RMSE and the values of d approaching 1 indicate the accurate match of the proposed equations (Equations (7) and (14)) to the numerical data for all porous media studied.
It may be noted that the hydraulic properties S and K0 present in Equations (7) and (14) can be determined from an available cumulative infiltration curve for vertical infiltration under zero ponding head [15]. Thus, Equation (7) contains one unknown, nx, which depends on the value of ponding head hp (Figure 1a–f) and Equation (14) contains the unknown variable nz which also depends on the value of ponding head hp as shown in Figure 2a–f. It is seen from these figures that nz > nx, and the rather small difference between them may be explained on the basis of the relationships (7) and (10). We have not reached an analytical expression relating hp and n for both cases of infiltration (horizontal and vertical). However, it is found that the plotting of loghp versus logn produces, approximately, a straight line (Figure 3). Thus, if two pairs (hpn)1 and (hpn)2 are determined, then the relationship (hp, n) may be approached for all such pairs.
The vertical cumulative infiltration data iz,p with time t obtained by the proposed Equation (14) and the Valiantzas equation (Equation (17)) for hp = 0 and various soil ponding heads hp, are presented in Figure 4 and the index of RMSE values between the two equations are presented in Table 3.
From Figure 4 and the RMSE values, it is obvious that there are differences in the estimation of vertical cumulative data iz,p between the two procedures. This difference is not systematical, since for hp = 0 the Val (Equation (17)) compared to Equation (14) overestimated the vertical cumulative infiltration data iz,p while for hp > 0 the Val generally underestimated the iz,p values.
Specifically, the Val systematically overestimated the iz,p data for hp = 0 for all the tested porous media. There is a higher difference between the two equations for the fine-textured soils (YCL, SiL) compared to the coarse-textured ones (LS, S). The overestimation of the iz,p data for hp = 0, as has been shown by Poulovassilis and Argyrokastritis [15], is attributed to the non-incorporation of the actual contribution of the pressure head gradient term to this model.
For all the other ponding heads hp > 0, the Val generally underestimated the iz,p data for all the tested porous media. The underestimation of Val compared to Equation (14) presents higher differences for the coarse-textured soils (S, LS), for which there is also a gradual increase of this difference as the hp values increase in the same porous media. For the finer porous media, the Val generally underestimated the iz,p but for an intermediate value of hp, different for each porous medium, the Val gives similar results to Equation [14]. The higher differences between the two equations (Equation (14) and Val) for the ponding depths hp > 0 are observed for coarser materials compared to the finer ones. Specifically, the underestimation of Val compared to Equation (14) for hp > 0 may be attributed to the estimation of Sp for the different hp and also in the way that the needed parameter of θi is estimated.

4. Conclusions

The study of infiltration under ponded conditions is of great interest in agricultural lands. However, there are no studies that have considered the actual contribution of the pressure head gradient to the flow in their analyses.
To fill that scientific gap, two simple equations are proposed, to describe the horizontal and vertical infiltration under various ponding heads, incorporating the actual contribution of the pressure head gradient to the flow.
The accuracy of the proposed equations was validated using horizontal and vertical numerical infiltration data for six soils, covering a wide range of soil textures.
The soil sorptivity and saturated hydraulic conductivity, which are included in the proposed equations, can be determined from an available vertical cumulative infiltration under zero ponding head. Thus, both the horizontal and vertical proposed equations contain one unknown parameter each, nx and nz respectively, which depend on the value of ponding head hp. The plotting of loghp versus logn produces, approximately, a straight line.
Finally, the proposed vertical infiltration equation was compared to that of Valiantzas which does not incorporate the actual contribution of the pressure head gradient to the flow and differences were observed in all porous media studied. The vertical cumulative infiltration data obtained by Valiantzas equation are overestimated for hp = 0 and are underestimated for hp > 0.

Author Contributions

Conceptualization, I.A.; methodology, I.A.; software, I.A., M.P. and P.A.L.; validation, I.A., M.P. and P.A.L.; formal analysis, I.A., M.P. and P.A.L.; investigation, I.A., M.P. and P.A.L.; data curation, I.A., M.P. and P.A.L.; writing—review and editing, I.A., M.P. and P.A.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Philip, J.R. The theory of infiltration: 6 Effect of water depth over soil. Soil Sci. 1958, 85, 278–286. [Google Scholar] [CrossRef]
  2. Philip, J.R. The theory of infiltration: 7. The infiltration equation and its solution. Soil Sci. 1958, 85, 333–337. [Google Scholar] [CrossRef]
  3. Green, W.A.; Ampt, G.A. Studies on soils physics: 1. The flow of air and water through soils. J. Agric. Sci. 1911, 4, 1–24. [Google Scholar]
  4. Brutsaert, W. Vertical infiltration in dry soil. Water Resour. Res. 1977, 13, 363–368. [Google Scholar] [CrossRef]
  5. Philip, J.R. The theory of infiltration: 1. The infiltration equation and its solution. Soil Sci. 1957, 83, 345–357. [Google Scholar] [CrossRef]
  6. Parlange, J.-Y.; Lisle, I.; Braddock, R.D.; Smith, R.E. The three-parameter infiltration equation. Soil Sci. 1982, 133, 337–341. [Google Scholar] [CrossRef]
  7. Parlange, J.-Y.; Haverkamp, R.; Touma, J. Infiltration under ponded conditions: 1. Optimal analytical solution and comparison with experimental observations. Soil Sci. 1985, 139, 305–311. [Google Scholar] [CrossRef]
  8. Haverkamp, R.; Parlange, J.-Y.; Starr, J.L.; Schmitz, G.; Fuentes, C. Infiltration under ponded conditions: 3. A predictive equation based on physical parameters. Soil Sci. 1990, 149, 292–300. [Google Scholar] [CrossRef]
  9. Swartzendruber, D. A quasi-solution of Richards’ equation for the downward infiltration of water into soil. Water Resour. Res. 1987, 23, 809–817. [Google Scholar] [CrossRef]
  10. Stroosnijder, L. Infiltratie en Herverdeling van Water in Gronde (Infiltration and Redistribution of Water in Soils). Versl. Landbouwkd. Onderz. 847. 1976. Available online: https://core.ac.uk/download/pdf/29386787.pdf (accessed on 6 September 2021).
  11. Valiantzas, J.D. New linearized two-parameter infiltration equation for direct determination of conductivity and sorptivity. J. Hydrol. 2010, 384, 1–13. [Google Scholar] [CrossRef]
  12. Talsma, T.; Parlange, J.-Y. One-dimensional vertical infiltration. Aust. J. Soil Res. 1972, 10, 143–150. [Google Scholar] [CrossRef]
  13. Haverkamp, R.; Kutilek, M.; Parlange, J.-Y.; Rendon, L.; Krejca, M. Infiltration under ponded conditions: 2. Infiltration equations tested for parameter time-dependence and predictive use. Soil Sci. 1988, 145, 317–329. [Google Scholar] [CrossRef]
  14. Barry, D.A.; Parlange, J.-Y.; Haverkamp, R.; Ross, J.P. Infiltration under ponded conditions: 4. An explicit predictive infiltration formula. Soil Sci. 1995, 160, 8–17. [Google Scholar] [CrossRef]
  15. Poulovassilis, A.; Argyrokastritis, I. A new approach for studying vertical infiltration. Soil Res. 2020, 58, 509–518. [Google Scholar] [CrossRef]
  16. Philip, J.R. The theory of infiltration: 4. Sorptivity and algebraic infiltration equations. Soil Sci. 1957, 84, 257–264. [Google Scholar] [CrossRef]
  17. Philip, J.R. Numerical solution of equations of the diffusion type with diffusivity concentration dependent. Trans. Faraday Soc. 1955, 51, 885–892. [Google Scholar] [CrossRef]
  18. Mualem, Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 1976, 12, 513–522. [Google Scholar] [CrossRef] [Green Version]
  19. Van Genuchten, M.T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef] [Green Version]
  20. Poulovassilis, A.; Elmaloglou, S.; Kerkides, P.; Argyrokastritis, I. A variable sorptivity infiltration equation. Water Resour. Manag. 1989, 3, 287–298. [Google Scholar] [CrossRef]
  21. Šimůnek, J.; Šejna, M.; Saito, H.; Sakai, M.; van Genuchten, M.T. The HYDRUS-1D Software Package for Simulating the One-Dimensional Movement of Water, Heat, and Multiple Solutes in Variably-Saturated Media, Version 4.0, Hydrus Series 3; Department of Environmental Sciences, University of California Riverside: Riverside, CA, USA, 2008. [Google Scholar]
  22. Brooks, R.H.; Corey, A.T. Hydraulic Properties of Porous Media; Hydrology Paper No. 3; Colorado State University: Fort Collins, CO, USA, 1964. [Google Scholar]
  23. Vogel, T.; Císlerová, M. On the reliability of unsaturated hydraulic conductivity calculated from the moisture retention curve. Transp. Porous Media 1988, 3, 1–15. [Google Scholar] [CrossRef]
  24. Kosugi, K. Lognormal distribution model for unsaturated soil hydraulic properties. Water Resour. Res. 1996, 32, 2697–2703. [Google Scholar] [CrossRef]
  25. Durner, W. Hydraulic conductivity estimation for soils with heterogeneous pore structure. Water Resour. Res. 1994, 32, 211–223. [Google Scholar] [CrossRef]
  26. Willmott, C.J. Some comments on the evaluation of model performance. Bull. Am. Meteorol. Soc. 1982, 62, 1309–1313. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Comparative presentation of the horizontal cumulative infiltration numerical data ix,p,n and predicted ones ix,p using Equation (7) with time t for various ponding heads hp, for Yolo light clay (a), clay loam (b), silty loam (c), sandy loam (d), loamy sand (e) and sand (f).
Figure 1. Comparative presentation of the horizontal cumulative infiltration numerical data ix,p,n and predicted ones ix,p using Equation (7) with time t for various ponding heads hp, for Yolo light clay (a), clay loam (b), silty loam (c), sandy loam (d), loamy sand (e) and sand (f).
Water 13 03492 g001aWater 13 03492 g001b
Figure 2. Comparative presentation of the vertical cumulative infiltration numerical data iz,p,n and predicted ones iz,p using Equation (14) with time t for various ponding heads hp, for Yolo light clay (a), clay loam (b), silty loam (c), sandy loam (d), loamy sand (e) and sand (f).
Figure 2. Comparative presentation of the vertical cumulative infiltration numerical data iz,p,n and predicted ones iz,p using Equation (14) with time t for various ponding heads hp, for Yolo light clay (a), clay loam (b), silty loam (c), sandy loam (d), loamy sand (e) and sand (f).
Water 13 03492 g002aWater 13 03492 g002b
Figure 3. Relationships between (a) the parameter nx of the Equation (7) and ponding head hp for horizontal infiltration and (b) the parameter nz of the Equation (14) and ponding head hp for vertical infiltration (b), for all porous media studied.
Figure 3. Relationships between (a) the parameter nx of the Equation (7) and ponding head hp for horizontal infiltration and (b) the parameter nz of the Equation (14) and ponding head hp for vertical infiltration (b), for all porous media studied.
Water 13 03492 g003
Figure 4. Comparative presentation of the vertical cumulative infiltration data iz,p predicted by the proposed Equation (14) and the Valiantzas equation [11] with time t for various ponding heads hp, for Yolo light clay (a), clay loam (b), silty loam (c), sandy loam (d), loamy sand (e) and sand (f).
Figure 4. Comparative presentation of the vertical cumulative infiltration data iz,p predicted by the proposed Equation (14) and the Valiantzas equation [11] with time t for various ponding heads hp, for Yolo light clay (a), clay loam (b), silty loam (c), sandy loam (d), loamy sand (e) and sand (f).
Water 13 03492 g004
Table 1. The Mualem-van Genuchten closed-form model parameters θs, θr, α, n, the tortuosity factor l, the saturated hydraulic conductivity Ks, the sorptivity S and the parameters c (Equation (12)) and k (Equation (15)) of the six porous media studied.
Table 1. The Mualem-van Genuchten closed-form model parameters θs, θr, α, n, the tortuosity factor l, the saturated hydraulic conductivity Ks, the sorptivity S and the parameters c (Equation (12)) and k (Equation (15)) of the six porous media studied.
SoilSoil ParametersInfiltration Parameters
θsθrαnlKsSck
(cm3 cm−3)(cm3 cm−3)(cm−1)(-)(-)(cm h−1)(cm h−1/2)(-)(cm h−1)
YLC0.4950.1240.0152.00.50.044280.75120.64680.024
CL0.410.0950.0191.310.50.260.74530.98980.1969
SiL0.4500.0670.021.410.50.451.23951.0660.2986
SL0.410.0650.0751.890.54.420833.67750.731.9899
LS0.410.0570.1242.280.514.59176.00060.58835.9645
S0.430.0450.1452.680.529.708.90920.5069.9463
Table 2. Root Mean Square Error (RMSE) and Index of Agreement (d) for testing the accuracy of the proposed Equations (7) and (14) with the numerical data for various ponding heads hp for the six porous media studied.
Table 2. Root Mean Square Error (RMSE) and Index of Agreement (d) for testing the accuracy of the proposed Equations (7) and (14) with the numerical data for various ponding heads hp for the six porous media studied.
hp (cm)YLCCLSiLSLLSS
RMSEdRMSEdRMSEdRMSEdRMSEdRMSEd
Horizontal infiltration using Equation (7)
00.0360.9990.0250.9990.0110.9990.0150.9990.0250.9990.0830.999
50.0460.9990.1690.9880.0031.0000.0300.9990.0690.9990.1510.999
100.0520.9990.2670.9800.0021.0000.0470.9990.0920.9990.2250.999
250.0640.9990.4640.9700.0011.0000.0750.9990.1270.9990.2420.999
500.0790.9990.6690.9660.0120.9990.1030.9990.1570.9990.2900.999
1000.1010.9990.9160.9660.0720.9990.1320.9990.1810.9990.3320.999
1500.1170.9990.5230.9930.1370.9990.1490.9990.1850.9990.3530.999
2000.1290.9990.6670.9910.0260.9990.1590.9990.1970.9990.3620.999
Vertical infiltration using Equation (14)
00.0450.9990.0420.9980.0590.9990.0170.9990.0720.9990.1970.999
50.0640.9990.0780.9960.0451.0000.1170.9990.1040.9990.1760.999
100.0910.9990.1260.9930.0580.9990.1400.9990.1120.9990.2150.999
250.1270.9990.2330.9870.0720.9990.1710.9990.1310.9990.1920.999
500.1610.9990.3430.9840.0790.9990.1990.9990.1480.9990.2120.999
1000.1850.9990.4520.9840.0890.9990.2290.9990.1550.9990.2380.999
1500.2110.9990.3180.9950.0990.9990.2460.9990.1480.9990.2520.999
2000.1980.9990.2170.9980.1140.9990.2590.9990.1490.9990.2620.999
Table 3. Root Mean Square Error (RMSE) values from comparing the vertical cumulative infiltration data iz,p predicted by the proposed Equation (14) and the Valiantzas Equation [11] for various ponding heads hp for the six porous media studied.
Table 3. Root Mean Square Error (RMSE) values from comparing the vertical cumulative infiltration data iz,p predicted by the proposed Equation (14) and the Valiantzas Equation [11] for various ponding heads hp for the six porous media studied.
hp (cm)YLCCLSiLSLLSS
00.7570.1940.7280.5790.4410.472
50.2880.0430.1590.0250.3790.696
100.1190.1290.0510.1820.6011.319
250.1730.2780.1910.3410.8251.536
500.4070.4230.2820.3980.9261.783
1000.5550.6170.3560.4070.9781.935
1500.6410.6030.4490.3960.9981.992
2000.6110.6060.4870.3790.9942.021
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Argyrokastritis, I.; Psychogiou, M.; Londra, P.A. Infiltration under Ponded Conditions. Water 2021, 13, 3492. https://doi.org/10.3390/w13243492

AMA Style

Argyrokastritis I, Psychogiou M, Londra PA. Infiltration under Ponded Conditions. Water. 2021; 13(24):3492. https://doi.org/10.3390/w13243492

Chicago/Turabian Style

Argyrokastritis, Ioannis, Maria Psychogiou, and Paraskevi A. Londra. 2021. "Infiltration under Ponded Conditions" Water 13, no. 24: 3492. https://doi.org/10.3390/w13243492

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop