Next Article in Journal
Participatory Modelling of Surface and Groundwater to Support Strategic Planning in the Ganga Basin in India
Next Article in Special Issue
Optimization of Spring Wheat Irrigation Schedule in Shallow Groundwater Area of Jiefangzha Region in Hetao Irrigation District
Previous Article in Journal
Assessment of Retrofitted Ramped Weirs to Improve Passage of Potamodromous Fish
Previous Article in Special Issue
Dynamics and Distribution of Soil Salinity under Long-Term Mulched Drip Irrigation in an Arid Area of Northwestern China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of the Flux–Concentration Relation for Horizontal Flow in Soils

Department of Natural Resources Management and Agricultural Engineering, Agricultural University of Athens, 75 Iera Odos, 11855 Athens, Greece
*
Author to whom correspondence should be addressed.
Water 2019, 11(12), 2442; https://doi.org/10.3390/w11122442
Submission received: 6 October 2019 / Revised: 11 November 2019 / Accepted: 19 November 2019 / Published: 21 November 2019
(This article belongs to the Special Issue Study of the Soil Water Movement in Irrigated Agriculture)

Abstract

:
The objective of the present work is to investigate the flux–concentration (F(Θ)) relation, where Θ is the normalized soil volumetric water content for the case of one-dimensional horizontal flow, subject to constant concentration conditions. More specifically, the possibility of describing F(Θ) by an equation of the form F(Θ) = 1 − (1 − Θ)p+1 is examined. Parameter p is estimated from curve-fitting of the experimentally obtained λ(Θ) data to an analytic expression of the form (1 − Θ)p where λ is the well-known Boltzmann transformation λ = xt−0.5 (x = distance, t = time). The results show that the equation of (1 − Θ)p form can satisfactorily describe the λ(Θ) relation for the four porous media tested. The proposed F(Θ) function was compared with the limiting F(Θ) function for linear and Green–Ampt soils and to the actual F(Θ) function. From the results, it was shown that the proposed F(Θ) function gave reasonably accurate results in all cases. Moreover, the analytical expression of the soil water diffusivity (D(Θ)) function, as it was obtained by using the equation for λ(Θ) of the form (1 − Θ)p, appears to be very close to the experimental D(Θ) data (root mean square error (RMSE) = 0.593 m2min−1).

1. Introduction

For a rational and productive irrigation water application, one needs to know well how water infiltrates into the soil, and how the soil moisture profiles evolve. This can be achieved by solving the appropriate Richards’ equation, which is feasible only when major soil hydraulic properties are known, together with the initial and boundary conditions imposed in each case.
To solve the unsaturated flow equation, many procedures have been proposed, e.g., analytical, semi analytical, finite differences, and finite elements. Philip [1] introduced the flux–concentration (F(Θ)) relation in order to develop a quasi-analytical technique for solving the unsaturated flow equation. The technique is rather simple and gives a holistic view of the phenomenon [1,2].
F(Θ) expresses the relation of soil water flux–density with the normalized volumetric soil water content, and can be used in an attempt to elucidate further the flow phenomena in unsaturated soils, both in the horizontal (absorption, as well as desorption) and in the vertical cases. The form of the F for the case of horizontal absorption and the normalized volumetric soil water content is given by the following equations:
F = q q 0 ,   Θ = θ θ i n θ 0 θ i n
where θin is the initial volumetric soil water content of the soil column, considered uniform throughout the column; θ0 is the volumetric soil water content at the soil surface (x = 0); q is the soil water flux density at a position x where soil water content is θ; and q0 is the soil water flux density at x = 0, the soil surface, where water is absorbed under constant concentration conditions with θ = θ0 = θs (θs = volumetric soil moisture at saturation) when the soil water pressure head is H = H0 = 0 cm column of water.
The values of F and Θ lie in the range of 0 to 1. A zero value for F and Θ when q = 0 is far enough from the soil infiltration surface, where θ = θin; the value is 1 at the soil surface, where q = q0 and θ = θ0 [1]. In general, the F is a function of θ, θin, θ0, and time (t) and it is, in most cases, unknown a priori [3].
For the prediction of the flow regime in porous media, using the flux–concentration relationship, either an iterative procedure should be applied, or accurate estimates of F(Θ) where feasible are used [2].
A number of researchers have proposed analytical expressions for the F(Θ) relationship to solve the flow equation for the case of the horizontal absorption under constant pressure head conditions, applied at the soil surface (x = 0).
Smiles et al. [4] proposed an empirical expression to describe the measured F(Θ) of a fine sand soil:
F ( Θ ) = 1 ( 1 Θ ) 1.19 .
Vauclin and Haverkampt [5] proposed the following simple expression
F ( Θ ) = 2 Θ Θ + 1 ,
while Evangelides et al. [6] proposed the empirical expression
F ( Θ ) = 2 Θ     [ 1     ( 1     Θ ) m ]
with the value of parameter m lying in the range 0.718–0.867, and the value m = 0.8 considered as the most appropriate, which is very close to the mean value.
Recently, Ma et al. [7] proposed the expression
F ( Θ ) = Θ [ 1   +   a 2 ( 1 Θ 1 a 2 ) ]
with the parameter a2 being a function of the initial soil water saturation and the soil pore structure index that reflects the shape of the soil water retention curve.
From the abovementioned expressions, those where F is a function of Θ with a parameter for soil properties (Equations (4) and (5)), compared to those where F is a function of Θ without any parameters (Equations (2) and (3)), are more flexible and can simulate F(Θ) relationships with higher accuracy for a wide range of texture soils [7].
Philip [1] investigated the behavior of F for two distinct cases of soils, characterized by their D(θ) functions—i.e., soils, where their D(θ) function is a delta function (the so called Green–Ampt soils, where the moisture profiles advance step-wise), as well as soils whose D(θ) function is constant, independent of moisture content θ (the so called linear soils). In real soils, the F(Θ) relationship will lie in between the Green–Ampt and linear soils [1].
Clothier et al. [8] investigated the possibility of estimating the expression λ (Θ), where λ   is the Boltzmann transformation λ ( θ ) = x t 0.5 for each soil from the experimental data of a horizontal absorption experiment, according to the equation below (Philip [9], Table 1, no. 2):
λ ( Θ ) = λ i ( 1     Θ ) p ,
where λi is the maximum value of λ(Θ) (i.e., when Θ = 0). Parameter λi considered to be a characteristic measure of the wetted region in a horizontal absorption experiment, and is influenced only by θin for every soil [10,11,12]; p is a fitting parameter.
Equation (6) has the advantage of an easy estimation of parameters λi and p; by a log-transformation, Equation (6) becomes linear, with the slope equal to p and transect equal to logλi when the λ(Θ) is described by Equation (6). This procedure appears preferable to least squares fitting.
It is to be mentioned here that previous similar studies [6] could obtain the F(Θ) relationship using more complex expressions for it. In the present study, we use a simple, mono-parametric F(Θ) expression using a mono-parametric λ(Θ) expression, which adequately describes the experimental data [8].
In this context, the main objectives of this study are (1) tο test the possibility of the mono-parametric λ(Θ) expression to describe the experimental data, as the equation λ ( Θ ) = λ i ( 1   Θ ) p through the proper selection of a fitting parameter p; and (2) tο investigate whether the experimental F(Θ) relations are described accurately enough from the analytical, mono-parametric function F ( Θ ) = 1 ( 1 Θ ) p + 1 for a wide range of soil types, as well as finding out the physical meaning of the parameter p and the range of its values.

2. Theory

Applying Darcy’s law and the mass conservation principle for the one-dimensional horizontal infiltration case, one easily gets the partial differential equation
θ t = x [ D ( θ ) θ x ] ,
where x is the horizontal axis, t the time, and D(θ) the soil water diffusivity. The solution of Equation (7) under the following initial and boundary conditions:
(a)
t = 0, x > 0, θ = θin
(b)
t > 0, x = 0, θ = θ0
(c)
t > 0, x → ∞, θ = θin
could be obtained in terms of the Boltzmann transformation λ ( θ ) = x t 0.5 , considered a function of θ only as
θ i n θ λ ( θ ) d θ = 2 D ( θ ) d θ d λ ,
following the initial and boundary conditions θ = θin, λ → ∞, and θ = θ0, λ = 0.
Introducing the Boltzmann transformation function λ ( θ ) = x t 0.5 in Darcy’s law, one gets
q = Κ ( θ ) d H d x = D ( θ ) d θ d λ λ x = t 0.5 D ( θ ) d θ d λ ,
and the combination of Equations (1), (8), and (9) gives the following expression for F
F = θ i n θ λ ( θ ) d θ θ i n θ 0 λ ( θ ) d θ ,
From the abovementioned equation, it is obvious that function F, for the case of the horizontal absorption, is independent of time and depends on the values of the volumetric soil water, θ, θin and θ0.
When the D(θ) function is a delta function, the flux concentration relationship F for the horizontal absorption is given by
F = Θ ,
while for the case where the D(θ) function is constant, the flux concentration relationship F is given by
F = exp [ ( inverfc ( Θ ) ) 2 ] ,
White et al. [13] have shown that Equation (12) could be reasonably approximated by the expression
F = sin [ ( π 2 ) Θ π 4 ] .
It is obvious from Equations (8) and (10) that the F(Θ) relationship depends on the form of D(θ).
From Equation (6), the soil sorptivity Sc function could be estimated by analytical expression [8], since
S c = ( θ 0     θ i n ) 0 1 λ i ( 1     Θ ) p d Θ = ( θ 0     θ i n ) λ i p   +   1 ,
From the combination of Equations (6), (9), and (10), it could be shown that the F(Θ) relationship could be estimated analytically, according to the mono-parametrical analytical expression
F ( Θ ) = 1 ( 1 Θ ) p + 1 .
A similar expression for F(Θ) was presented by Smiles et al. [4] in a horizontal absorption experiment in sand, where the parameter p had a constant value p = 0.19.
Applying the Bruce and Klute [14] method, the soil water diffusivity function D(Θ) can be estimated according to the expression
D ( Θ ) = 0 Θ λ ( Θ ) d Θ 2 d Θ d λ .
In this respect, it is rather easy, by introducing the λ(Θ) expression given in Equation (6), to get an analytical expression for D(Θ) ([8]), as
D ( Θ ) = p ( p   +   1 ) S 2 [ ( 1     Θ ) p 1 ( 1     Θ ) 2 p ] 2 ( θ 0 θ i n ) 2 .
In practice, this means that for soils with a λ(Θ) expression, such as the one given by Equation (6), their diffusivity function D(Θ) should be described analytically by Equation (17), and the flux concentration relationship F(Θ) by Equation (15).

3. Materials and Methods

Horizontal infiltration experiments were conducted (Scheme 1) in three porous media with different soil textures: a sandy loam (SL) (13.2% clay, 8% silt, 78.8% sand), a loam (L) (20% clay, 38% silt, 42% sand), and a silty clay loam (SiCL) (36.5% clay, 52% silt, 11.5% sand). The respective soils were air-dried, ground, and passed through a 2 mm sieve. The soil samples were uniformly packed in a rectangular column 5 cm in width, 2 cm in height, and 50 cm long. The upper soil surface was open to the air, thus securing that the soil air pressure at the soil pores was at atmospheric pressure (Scheme 1). The experimental data from horizontal absorption experiments conducted in sand (S) by Poulovassilis [15] were also used.
The volumetric soil water θj(xj, t*) profile, after a time period t* was elapsed from the beginning of the experiment, was determined by cutting, as quickly as possible, the soil column into small rectangular pieces. Specifically, the soil column was sectioned in 0.01 m increments, and the volumetric water content of each rectangular piece θj was determined by using the gravimetric water content and the dry soil bulk density (ρφ). The θj value corresponds to the center of the soil samples, at distance xj from the soil infiltration surface.
The initial values of θ and θin were those associated to the air-dry soil, and the condition of the soil infiltration surface was that of a constant pressure head H value (H = 0 at x = 0, with x = 0 denoting the soil infiltration surface). The zero value of H was maintained by a Mariotte device (Scheme 1). A thin wire mesh at x = 0 was installed to keep the soil at rest, and also provide the least possible resistance to the soil water entry.
The time duration of these experiments varied according to the soil type—i.e., smaller values for coarse-textured soils and larger values for fine-textured soils. For the sand (S), it was 10.4 min [15], for the SL it was 194 min, for the L it was 251 min, and for the SiCL it was 1630 min. During the experimental process, a continuous monitoring of the water volume entering the column was obtained, thus allowing the cumulative infiltration I(t) experimental curve to be determined. The soil sorptivity S is immediately available from the well-known relationship [16]
S = I t 0.5 ,

4. Results and Discussion

In Table 1, some soil characteristics, such as θ0, θin, ρφ, λi and S (Equation (18)), together with fitting parameter p and the value of Sc (Equation (14)), are shown for each porous media used in this study. It is easily noticeable that the values of S and λi depend strongly on the soil type, and they tend to decrease as soils become finer in texture. The same trend for the values S and λi for six different soils were presented from McBride and Horton [10].
In Figure 1, the λ(θ) profiles are presented, as these were obtained from the experimental θ(x,t) data and the well-known Boltzmann transformation ( λ ( θ ) = x t 0.5 ), together with the fitted λ(θ) relationship (Equation (6)), after the fitting parameter p was properly selected. From the comparison, it is shown that Equation (6) is suitable for the description of the λ(θ) relationships for the horizontal absorption experiments. In each soil, the p-value was determined as the one where the root mean square error (RMSE) from a series of neighboring p-values was the least.
One could argue, considering that the λ(Θ) relationship is unique for each soil, independent of the duration of the experiment, that the two parameters of Equation (6) (λi and p) do not represent simple fitting parameters, but are related to the soil’s hydraulic properties. Some researchers insist that parameter λi may be considered as an index of the soil’s hydraulic properties [10,11]. Moreover, Shao and Horton [11] correlated parameter p with the parameter n (p = 1/n) of the equation of van Genuchten [17], which describes the soil moisture retention curve.
In what follows, an investigation of the characteristics of parameter p will be carried out. According to Equation (14), the ratio S/λi(θ0θin) is equal to (p + 1)−1. For the Green–Ampt soils (D is a Dirac delta-function of θ), the ratio S/λi(θ0θin) should be unit, thus p would be zero. Similarly, for the linear soils, the expression (p + 1)−1 would be 0.31, and the value of p will be 2.23 [8]. Consequently, the values of p are related to the form of the diffusivity D(θ) function. In order to examine this relationship, Equation (6) is rewritten as in Equation (19):
λ λ i = ( 1 Θ ) p Θ = 1 ( λ λ i ) 1 p .
In Figure 2, the Θ(λ/λi) relationship is shown for various values of the parameter p (0 < p < 2.23). From Figure 2, it is shown that all the Θ(λ/λi) relations lie in the area with its borders defined by the values of the parameter p (i.e., p 0 ; D(θ) Dirac delta function) and p = 2.23 (D(θ) constant). From the investigation of the λ(θ) relations in this study, it is found that p-values fall in the range 0.1 < p < 0.4. Also, Evangelides et al. [6] showed that p-values obtained using data from horizontal infiltration experiments in seven soils fall in the range 0.149 < p < 0.389. In other words, the range of p is narrower than the range 0 p 2.23 .
After the parameter p selection, using also the λi, θ0, and θin values for each soil, a re-estimation of sorptivity Sc was performed using Equation (14). From the values of S and Sc presented in Table 1, Sc values are reasonably close to the experimental values of S (S = I/t0.5) for three out of the four soils examined (absolute values of RE: 1.55% < RE < 3.08%). The relative error value for the SiCL soil appears to be rather high (13.27%). This could be attributed to the long time duration of the experiment (1630 min) and unavoidable soil water evaporation from the upper soil surface. In any case, the overall differences are small, and therefore one may consider that Equation (14) can lead to a quick and reliable way of estimating S from a set of horizontal absorption experimental data.
In Figure 3, the D(θ) data for the loamy soil (L) are shown. Closed circles denote data obtained by the Bruce and Klute [14] method (Equation (16)), while open circles are data obtained from Equation (17). The comparison of the above indicates that Equation (17) could reasonably describe the experimental D(θ) data. One should note that by using the analytical expression (Equation (17)), the problem of differentiating the experimental data d Θ d λ , in which there is scatter, is overcome. Moreover, the problem of estimating the slope d Θ d λ near saturation, where the λ(Θ) relationship is almost parallel to the λ-axis, is overcome by the application of the analytical expression (Equation (17)). Similar results were also obtained for the other soils under present investigation.
In Figure 4, the F(Θ) relationship for the linear Equation (13) and the Green–Ampt soils (Equation (11)) is shown, together with that obtained according to Equation (15) for all soils examined. The experimental F(Θ) points for all soils were obtained by the application of Equation (10) and the measured value of the sorptivity S (Equation (18)). From the results, it is shown that the one-parameter Equation (15) gave practically the same values for the F(Θ) relationships as the ones obtained experimentally, and lies between the two limits (linear and Green–Ampt soils) in all cases that were examined.
It is worth investigating the ability of the equation F(Θ) = 1 − (1 – Θ)p+1 to describe the upper and lower theoretical boundaries of F(Θ). As has already mentioned, when p = 0 the value of F = Θ—i.e., it converges to the lower limit (D(θ) function is a delta function). However, the values of F for p = 2.23 (resulted from (p + 1)−1 = 0.31 [8]) differs from the values of F resulting from Equation (13), which is the upper theoretical limit of F(Θ). Specifically, it gives higher values of F(Θ) than the theoretical curve. The fitting of the F(Θ) = 1 − (1 – Θ)p+1 equation to the theoretical curve (Equation (13)) showed that the equation F(Θ) = 1 − (1 – Θ)2.36 gives very good results of F(Θ) estimation over the entire range of Θ [7]. Therefore, the variation range of the parameter p is between 0 and 1.36, and the resulting shape parameter value for the linear soils ((p + 1)−1 = 0.423) is greater than that presented by Clothier et al. [8].
It can be concluded that F(Θ) = 1 − (1 – Θ)p+1 seems to be appropriate functional form for describing actual flux–saturation curves of general soils, and its parameter p has a physical meaning, i.e., it represents the shape of the soil moisture profile.
Furthermore, one could observe that from the study of Evangelides et al. (Table 2 [6]), where seven different soils are studied, the fitting parameters m (Equation (4) [6]) and n = p + 1 (Equation (15)) are strongly linearly related, as n = p + 1 = 1.483 m + 2.4375 . If one uses m = 4/5, as Evangelides et al. [6] proposed, then this leads to n ¯ = 1.251 which corresponds to our p value 0.25. In this respect, Equation (4) [6], with m = 0.8, becomes equivalent to Equation (15) in the present study, with p = 0.25.

5. Conclusions

The mono-parametric equation of the form (1 – Θ)p can reliably describe the λ(Θ) relationship after the proper selection of the parameter p for a relatively large range of soils. It is also shown that the analytical expressions of the soil hydraulic diffusivity D(θ) and soil sorptivity S approach the experimental ones well. Moreover, for the case of the D(θ) relationship, there is the advantage of obtaining values near saturation, where the classical methodology of Bruce and Klute might be inadequate.
During the laboratory experiment, a new mono-parametric analytical function was used for F(Θ) evaluation, where the parameter p is estimated from curve-fitting of the experimentally obtained λ(Θ) data to an analytical expression of the form (1 – Θ)p. Parameter p seems to be strongly related to soil hydraulic properties, and further investigation is needed to find this exact relationship. The analytical F(Θ) relationship, for all soil types investigated, approaches the experimental ones very well, and lies within the limiting F(Θ) function for linear and delta function soils. In addition, the upper and lower limit curves of F(Θ) calculated by the proposed expression were consistent with theoretical curves.

Author Contributions

All authors contributed equally to this work.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Philip, J.R. On solving the unsaturated flow equation: 1. The flux-concentration relation. Soil Sci. 1973, 116, 328–335. [Google Scholar] [CrossRef]
  2. Philip, J.R.; Knight, J.H. On solving the unsaturated flow equation: 3. New quasi-analytical technique. Soil Sci. 1974, 117, 1–13. [Google Scholar] [CrossRef]
  3. White, I. Measured and approximate flux-concetration lrelations for absorption of water by soils. Soil Sci. Soc. Am. J. 1979, 4, 1074–1080. [Google Scholar] [CrossRef]
  4. Smiles, D.E.; Knight, J.H.; Perroux, K.M. Absorption of water by soil: The effect of a surface crust. Soil Sci. Soc. Am. J. 1982, 46, 476–481. [Google Scholar] [CrossRef]
  5. Vauclin, M.; Haverkampt, R. Solutions quasi analytiques de l’equation d’absorption de l’eau par les sols non satires. I. Analyse critique. Agronomie 1985, 5, 597–606. [Google Scholar] [CrossRef]
  6. Evangelides, C.; Tzimopoulos, C.; Arampatzis, G. Flux-Saturation relationship for unsaturated horizontal infiltration. Soil Sci. 2005, 170, 671–679. [Google Scholar] [CrossRef]
  7. Ma, D.H.; Zhang, J.B.; Lu, Y.X. Derivation and Validation of a New Soil Pore-Structure-Dependent Flux–Saturation Relationship. Vadose Zone J. 2017, 16, 1–15. [Google Scholar] [CrossRef]
  8. Clothier, B.E.; Scotter, D.R.; Green, A.E. Diffusivity and One-Dimensional Absorption Experiments. Soil Sci. Soc. Am. J. 1983, 47, 641–644. [Google Scholar] [CrossRef]
  9. Philip, J.R. General method of exact solutions in concentration-dependent diffusion equation. Aust. J. Phys. 1960, 13, 1. [Google Scholar] [CrossRef]
  10. McBride, J.F.; Horton, R. An empirical function to describe measured water distributions from horizontal infiltration experiments. Water Resour. Res. 1985, 21, 1539–1544. [Google Scholar] [CrossRef]
  11. Shao, M.; Horton, R. Integral method for estimating soil hydraulic properties. Soil Sci. Soc. Am. J. 1998, 62, 585–592. [Google Scholar] [CrossRef]
  12. Ma, D.H.; Wang, Q.; Shao, M.A. Analytical method for estimating soil hydraulic parameters from horizontal absorption. Soil Sci. Soc. Am. J. 2009, 73, 727–736. [Google Scholar] [CrossRef]
  13. White, I.; Smiles, D.E.; Peroux, K.M. Absorption of water by soil: The constant flux boundary condition. Soil Sci. Soc. Am. J. 1979, 43, 659–664. [Google Scholar] [CrossRef]
  14. Bruce, R.R.; Klute, A. The measurement of soil water diffusivity. Soil Sci. Soc. Am. Proc. 1956, 20, 458–462. [Google Scholar] [CrossRef]
  15. Poulovassilis, A. Flow characteristics during infiltration into a horizontal sand column. Water Resour. Res. 1977, 13, 369–374. [Google Scholar] [CrossRef]
  16. Philip, J.R. Theory of infiltration. Adv. Hydrosc. 1969, 5, 215–296. [Google Scholar]
  17. 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]
Scheme 1. Experimental device.
Scheme 1. Experimental device.
Water 11 02442 sch001
Figure 1. The λ(θ) relationships, as obtained directly from experimental θ(x,t) data for each soil, and the estimated values according to Equation (6).
Figure 1. The λ(θ) relationships, as obtained directly from experimental θ(x,t) data for each soil, and the estimated values according to Equation (6).
Water 11 02442 g001
Figure 2. The relationship Θ(λ/λi) for different values of the parameter p.
Figure 2. The relationship Θ(λ/λi) for different values of the parameter p.
Water 11 02442 g002
Figure 3. A comparative presentation of the relationship D(θ) as obtained according to the Bruce and Klute [14] method (Equation (16)) and Equation (17) for the loamy (L) soil.
Figure 3. A comparative presentation of the relationship D(θ) as obtained according to the Bruce and Klute [14] method (Equation (16)) and Equation (17) for the loamy (L) soil.
Water 11 02442 g003
Figure 4. Comparative presentation of the relationship F(Θ) for the linear soils and Green–Ampt soils, as well as the F(Θ) calculated according to Equation (15) and the experimental F(Θ) values described in the text for the four soils examined.
Figure 4. Comparative presentation of the relationship F(Θ) for the linear soils and Green–Ampt soils, as well as the F(Θ) calculated according to Equation (15) and the experimental F(Θ) values described in the text for the four soils examined.
Water 11 02442 g004aWater 11 02442 g004b
Table 1. The values of the soil characteristics θ0, θin, ρφ, and S, together with the values of the parameters λi, p, and Sc for each soil examined. The value of the parameter p was estimated from Equation (6). The value of sorptivity S was obtained directly from the experimental I(t) data (Equation (18)), while Sc comes from Equation (14). RE denotes the absolute relative error between actual S and estimated Sc values for each soil examined.
Table 1. The values of the soil characteristics θ0, θin, ρφ, and S, together with the values of the parameters λi, p, and Sc for each soil examined. The value of the parameter p was estimated from Equation (6). The value of sorptivity S was obtained directly from the experimental I(t) data (Equation (18)), while Sc comes from Equation (14). RE denotes the absolute relative error between actual S and estimated Sc values for each soil examined.
Porous Mediaθ0 (m3 m−3)θin (m3 m−3)ρφ (g cm−3) 1λi (cm min−0.5)pS (cm min−0.5)Sc (cm min−0.5)RE (%)
Sand (S)0.2840.06-12.2100.351.9301.9802.59
Sandy loam (SL)0.4180.0151.413.5100.121.3001.263.08
Loam (L)0.4650.0221.121.5450.10.6440.6341.55
Silty clay loam (SiCL)0.5110.0281.220.6560.180.3090.26813.27
1 Dry soil bulk density value for sand soil is not referred to by Poulovassilis [15].

Share and Cite

MDPI and ACS Style

Kargas, G.; Londra, P.; Kerkides, P. Investigation of the Flux–Concentration Relation for Horizontal Flow in Soils. Water 2019, 11, 2442. https://doi.org/10.3390/w11122442

AMA Style

Kargas G, Londra P, Kerkides P. Investigation of the Flux–Concentration Relation for Horizontal Flow in Soils. Water. 2019; 11(12):2442. https://doi.org/10.3390/w11122442

Chicago/Turabian Style

Kargas, George, Paraskevi Londra, and Petros Kerkides. 2019. "Investigation of the Flux–Concentration Relation for Horizontal Flow in Soils" Water 11, no. 12: 2442. https://doi.org/10.3390/w11122442

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