Abstract

A method based on a specific power-law relationship between the hydraulic head and the Boltzmann variable, presented using a similarity hypothesis, was recently generalized to a range of powers to satisfy the Bruce and Klute equation exactly. Here, considerations are presented on the proposed similarity assumption, and new analytical support is given to estimate the water density flux into and inside the soil, based on the concept of sorptivity and on Buckingham-Darcy's law. Results show that the new analytical solution satisfies both theories in the calculation of water density fluxes and is in agreement with experimental results of water infiltrating horizontally into sand. However, the utility of this analysis still needs to be verified for a variety of different textured soils having a diverse range of initial soil water contents.

1. Introduction

Based on physical laws of similarity applied to the rate of work required for water to wet and move through a soil, Prevedello et al. [1] presented an analytic solution of a Boltzmann transformed equation of continuity for horizontal infiltration which is derived without invoking the concept or use of the soil water diffusivity function. The derivation assumes that a similarity exists between the shapes of the soil water retention function and the Boltzmann transformation , and the solution successfully described soil water content profiles experimentally measured for different infiltration times into a homogeneous sand. More recently, an extension of this theory generalized the solution to a range of powers to include the saturated zone [2], to satisfy the Bruce and Klute equation exactly.

With a similar assumption, but not exactly as expressed by Prevedello et al. [1], Prevedello et al. [3] obtained a new analytic solution of the Richards equation for the infiltration into the same sand that holds for all infiltration times from zero to infinity, including vertical directions without making use of empirical constants. In this case, the derivation assumed that a similarity exists between the soil water retention function and the soil water content distribution within the soil profile during infiltration.

Although the similarity assumptions used by [1, 3] seem different, we now show that the mathematical development of both leads to the same equation for the description of the soil water content profiles for horizontal infiltration. We also present an analysis to obtain analytical equations to estimate the water density flow into and inside the soil, according to the sorptivity concept and the Buckingham-Darcy law. These analytical equations are compared with experimental results of horizontal infiltration into the same homogeneous marine sand.

2. Theoretical Considerations

2.1. Extended Similarity Hypothesis

Two identical solutions for the horizontal water infiltration into a soil were presented by [1, 3], coincidently (39) in both papers, but using two apparently different initial postulates. The explanation for that is given by a generalization of the similarity theory proposed in [1]. A generalization of the similarity between the Boltzmann transformation and the soil water retention function can be written as follows: in which is a positive integer, the soil water content (m3m-3), and the matric water potential head (). With (1) being the fundamental hypothetical equation, it follows that or and introducing the horizontal position coordinate and the time (s), we have note that .

When , From (2) it follows that

Substituting (7) into the Buckingham-Darcy equation written in the -direction in terms of the soil hydraulic conductivity as a function of , we have where the subscript refers to the model hypothesis of (4). Knowing (8), changing the variables so that the dependent variable becomes the distance , and remembering that because the partial derivative of is taken with respect to , we have which leads to For a fixed time, based on (7), and multiplying by , we have or Holding constant and integrating, the above equation becomes Evaluating at when , Recalling (3) in the form of and equating the two above equations that define , we have Returning to (4), we know that Substituting into the above equation, we obtain Recalling that , the solution is

For the cases in which the soil is initially very dry (i.e., ), the above equation reduces to which is (39) presented in [1, 3].

From (5) we have According to (21), the generalized solution of is But the fundamental hypothesis that gave rise to (22) is The derivative of (22) with respect to is Dividing (22) by ,

Equating (24) and (25), since both follow the postulate, the result is . This shows that any value of different from 2 cannot be accepted in this theory of similarity, as originally postulated and confirmed in [1], since it is with the square root of the time that the advancement of the wetting front remains linear in the process of horizontal infiltration.

2.2. Explanation Notes

Considering and since , [1, equation (19)], for the cases in which at and when the soil is initially very dry , can be written as Therefore, we know that and we can write (26) as follows: which shows that the gradient formulated as an initial postulate in [1] is simplified during the theoretical development to . This statement can also be verified starting with (11). In this equation, holding constant and making , we can write where we can see that the initially formulated gradient becomes . The integration of the equation above leads to (22) in [1, 3].

The solution obtained by [3] for horizontal infiltration case, although identical to that obtained by [1], presupposes a relation of the type with , that is, or or , and, therefore, it cannot be taken as formal as that presented in [1]. As explained above, this is due to the fact that independently of the chosen value of in the generalised similarity theory here introduced, the gradient of (7), as initially postulated in [1], is always simplified during the theoretical development from to , so that the B-D equation is always converted to , and that, therefore, leads to a final solution that is identical to that obtained by [1].

In this way it is shown that the same soil water profile , expressed by (22) in both papers [1, 3], can be obtained from both similarity approaches, that is, or as images of the soil water retention curve . However, only as an image of translates formally what was postulated into a solution, such as shown in [1] or in the extended similarity hypothesis shown above, in which only carries out formally this postulate.

2.3. The Similarity Hypothesis and New Analytical Support on the Estimation of Horizontal Infiltration

Since Prevedello et al. [1] found a solution such that is proportional to , so that is proportional to , it follows that for , from the definition of , at both and are equal to zero. At a first glance, one could assume that to have a finite flux density at , should be infinite at . According to this interpretation, Barry et al. [2] introduced a nonzero air-entry value (e.g., Haverkamp et al. [4]), such that for , the soil is unsaturated and saturated for , and generalized the Prevedello et al. [1] solution to a range of powers. We here use the original interpretation of Prevedello et al. [1] to calculate the flux density at and at , without running into the shortcomings mentioned above. From Phillip’s analysis, we know that for the case of horizontal infiltration,

If in (30) we substitute by the Bruce and Klute equation or (1) of Barry et al. [2] are obtained. There are many analytical forms to prove the equality in (30). Here, our option was to obtain an analytical expression for both the members and to compare the results obtained for each one. According to the similarity assumption proposed by [1] we know that From van Genuchten we know that so that (31) becomes so that

The water flux density can also be evaluated through the Buckingham-Darcy law directly, using the Prevedello et al. [1, 3] similarity assumption: But, according to the similarity assumption, each value of of the wetting profile is associated to its position value . Therefore, does not represent directly the total or representative matric gradient at the wetting front. In this way, (35) should be interpreted as the water density flow that enters at in response to a “representative” gradient that occurs at the wetting front. Thus, if is taken as the average matric potential of all the active potentials, then we can find analytically this representative matric potential, as well as the position where it acts at the wetting front, applying the similarity assumption. In this way, the previous equation can be written as where is the representative potential active at the wetting front and the position where this potential acts, defined by or

Notice that this interpretation provides a constant value of independent of time, and it follows that

From the soil surface, where the soil is saturated (), the wetting front advances as a response to an average matric potential , and with a force (average matric gradient), which should yield the same results as (34).

The analytical equation to estimate water density flow inside the soil can be determined simply by substituting for or into (40), because the average gradient remains fixed, but inside the soil we have the soil water content or matric potential for each position :

3. Material and Methods

To validate experimentally the equations proposed above to quantify the water density into and inside the soil during the horizontal infiltration process for homogeneous sand, we used the same data presented in [1, 3].

The van Genuchten [6] model was used to describe the soil water retention curve where , , , (residual soil water content), and are the classical parameters of the soil water retention curve. Using a nonlinear regression program as the one suggested by Boratto [5], (42) was fitted to experimental values to obtain the independent parameters , , , , and .

According to van Genuchten [6], introducing (42) into the Mualem [7] theory, the function is obtained as

Data for the saturated hydraulic conductivity and for the parameters of the water retention curve for the marine sand are presented in Table 1, as previously shown in [1, 3]. With a coefficient of determination greater than 0.99, the van Genuchten [6] model provides an excellent representation of the observed values.

4. Results and Discussion

4.1. The Water Density Flow into the Soil (), according to Philip’s Concept of Sorptivity

From the boundary conditions and the soil parameters (Table 1) and assuming the initial soil water content equal to m3·m-3, (34) reduces to which is shown in Figure 1 as a function of time, according to the analysis of Philip [8].

Measured soil water content profiles (Figure  5 of Prevedello et al. [1]) were integrated with respect to distance for times , , and equal to 15, 60, and 300 min, resulting in 165, 330, and 732 mm, respectively. We also know that the cumulative infiltration into the soil surface

is or where is the sorptivity.

In light of (46) and substituting the value  mm·(min)-1/2, we have  mm,  mm, and mm, respectively, which are very close to those estimated by integration of experimental data and showed in Figure 2.

4.2. The Water Density Flow into the Soil (), According to Buckingham-Darcy Law

To verify the analytical validity of (40) to quantify the water density flow at , (38) and (39) were evaluated, resulting, respectively, in  m and . These results show that the similarity theory proposed by [1], here represented by (36) or (40), does not imply in a water flux density zero at . Substituting these values and  m/min into (40), we obtain which is exactly the same value obtained in the calculation through Philip’s concept for sorptivity (44) and what therefore also adjusts very well with the experimental results showed in the Figure 2.

Another apparently analogous way to obtain the results of (47) above is to consider the derivative from (20), that is, . But according to the similarity theory proposed in [1, 3], the value of is a constant for each position at which it is only dependent on the time of advancement, and the derivative of a constant is null. The same restriction is applicable to the Green and Ampt theory [9] when we want to obtain the matric gradient through derivative of (the pressure head at the wetting front) with respect to because is always a constant independent of . So, the derivative of with respect to cannot be used here, except if associated to the concept of the similarity, and as exposed previously in item b above, we know that the matric gradient is reduced to , and therefore we can also write (40) as

Comparing (30) with (48), we know that for , for , and for . It is now readily apparent from the above analysis that the distribution of was substituted by the constant in (48) to determine the water flow density at in response to the force that occurs at distance where . Inside of the soil the values of or are consistent with the values of found from the values of in (39) of Prevedello et al. [1, 3] because (39) provides the shape of the soil water content profile at the wetting profile while the average gradient / remains active on the time.

Equation (48) with values of from Table 1 and also yields values very close to the experimental results (Figure 2) and exactly the same results obtained in the calculation through Philip’s concept for sorptivity (44).

4.3. The Density of Flow inside of the Soil, for Any Position and Time

Considering the boundary conditions of soil water content varying between 0.0187 and 0.387 m3·m-3, to use (41) we first estimate the corresponding matric potentials for any stipulated between 0.0187 and 0.387 m3·m-3, to obtain their respective positions from (20) and their respective ) values from (43).

Through (41) we generate the water density flow for each position inside the soil, that is,

Analogously as before in (48), the water density flow inside the soil can also be determined through the following: in which , and the values of and are showed in the Table 1.

Both (49a) and (49b) yield exactly the same values, and Figure 3 shows the water density flow estimated from those equations for each reached at 15, 60 and 300 minutes of the horizontal infiltration into the homogeneous marine sand.

It can be noted that for the water density flow values coincide with those obtained from (44), (47), or (48), as it should be.

From the above, it follows that

or

indicating that the gradient decreases with the time, shown in Figure 4 for times of 15, 60, and 300 minutes.

To verify how accurate these estimated results are, the same boundary value problem was solved numerically from Richard’s equation without any similarity assumption, using the procedure of Philip [10] with appropriate initial and boundary conditions and using the measured soil water properties shown in Section 3 (calculations based on  m3·m-3). Results showed that both water density flows at and inside the soil were very similar in the transmission zone and wetting front for each of the three chosen times of horizontal infiltration. However, the utility of this analysis still needs to be verified for a variety of different textured soils having a diverse range of initial soil water contents.

5. Conclusion

The new analytical supports to estimate the water density into and inside of the soil during horizontal infiltration agree with the experimental results, with the concept of sorptivity proposed by Philip [8] as well as the Buckingham-Darcy law, when a “representative” matric gradient is considered acting at the wetting front.

Acknowledgment

The authors acknowledge Donald R. Nielsen and Klaus Reichardt for their support and guidance in preparing this paper.