- Citado por SciELO
- Citado por Google
- Similares em SciELO
- Similares em Google
versão On-line ISSN 0718-9516
J. Soil Sci. Plant Nutr. vol.12 no.1 Temuco 2012
Journal of Soil Science and Plant Nutrition, 2012, 12 (1), 99-112
Water flow and conductivity into capillary and non-capillary pores of soils
Abdelmonem M. Amer1
1Department of Soil Science, Faculty of Agriculture, Menoufia University, 32511 Shebin EI-Kom, Egypt. E-mail address: firstname.lastname@example.org
Downward water flow in the vadose zone occurs principally through the non-capillary pores, while the redistribution and lateral and upward flow occurs in the capillary pores. The purpose of this study was to propose equations to estimate water flow, Q(θ), and hydraulic conductivity, K(θ), in the capillary and non-capillary soil pores. The equations related K(θ) to soil pore radius (r) were based on soil hydraulic data, including water retention h(q), field basic infiltration rate, water sorptivity (S) and distribution density function f(r) of the soil pore size. Calcareous sandy loam and alluvial clay soils located at the Nile Delta were used to test the validity of the assumed equations. Data showed that the values of K(θ) calculated by the proposed equations were in the common ranges for such soils. The equations are therefore expected to be applicable for both coarse and fine textured soils. Also, an equation was derived to estimate the sorptivity at steady state infiltration. It was observed that S is decreased in going from the un-saturation condition to steady state infiltration by 15.1% and 45.9% in sandy loam and clay soils, respectively.
Key words: Unsaturated hydraulic conductivity, soil pore size, water infiltration and sorptivity.
Mathematical models have become increasingly popular in the research and management of soil water flow and transport processes. The unsaturated hydraulic conductivity, K(θ), and the relationships between water content (θ), pressure head (h) and soil pore size are important to quantify the rate of water flow and transport through the vadose zone. Water moves into soil pores due to gradients in water content and potentials, and this is usually from a wet soil to a dry one. The ability of pores to conduct water is controlled mainly by pore sizes, continuity and distribution of the pores in the soil. The soil pore sizes could be classified into macro (non-capillary) pores, coarse capillary pores and fine capillary pores (FCP) (Baver and Gardner 1972). The non-capillary pores are represented as rapidly draining pores (RDP), while the coarse capillary pores (CCP) are represented by the slowly draining pores (SDP) and water holding pores (WHP) (Amer 2002). The pressure head corresponding with the cutoff between capillary and non-capillary pores varies widely, ranging from 0.1 kPa (Beven and Germann 1982) to 10 kPa (Marshall 1956). In the present study, h = 10 kPa has been selected as corresponding to the limit between capillary and non-capillary pores (Amer et al., 2009).
Theoretically, K(θ) may be predicted using statistical pore size distribution models, which assume water flow through " idealized " cylindrical pores. These models are described with the equations of Poiseuille and Darcy (Mualem 1976; Amer et al., 2009) in which water filled pores can be treated as capillary bundle tubes. The non-capillary pores have been recognized as significant pathways for water and solute movement in the soil. In most cases, macro (non-capillary) pores comprise only a small portion of the total soil voids, but under some conditions vertical flows through the macro pores dominate during infiltration (Chen and Wagenet 1992).
The aim of this work was to propose equations for predicting water flow and unsaturated hydraulic conductivity, K(θ), in capillary and non-capillary pores of soils. The approach is based on data on certain hydraulic parameters such as soil water retention y(q), distribution density function f(r) of pore sizes and sorptivity at steady infiltration rate. These parameters can be measured experimentally and the procedures are not as time-consuming as measuring K(θ) in situ. Moreover, results of measuring K(θ) in situ are often unreliable.
2. Theoretical background
2.1. Soil water retention and pore size classes Water is held in soil pores by cohesive and adhesive capillary forces. The size of the pores in unsaturated state can be determined through the so called hydraulic radius (r) of a section of pore space. The relationship between r and the capillary forces expressed as pressure head potential (h in m) is represented by the capillary rise equation (Hillel 1980):
where, r (in m) is the equivalent cylindrical pore size (hydraulic) radius related to the meniscus curvature radius (R) via the equation; r = R cos α, and cos α =1, (as it was assumed that the contact angle (α) is 0 for a wetted surface), y is the surface tension between water and air (at 20Â°C= 0.0727 kg s-2), g is the acceleration due to gravity (9.8 m s-2), pw is the density of water (at 20Â°C= 998kg m-3). As the soil dries out, increasing suction occurs due to the progressive emptying of capillary pores. Pore size diameters were determined for the ranges of soil matric potentials by applying Eq.(1) with respect to the soil water retention curves.
The equivalent pressure (Ψ) ranges of 0-10, 10-33, 10-1500, 33-1500, and > 1500 kPa, are roughly corresponding to rapidly draining pores (RDP), slowly draining pores (SDP), coarse capillary pores (CCP), water holding pores (WHP), and fine capillary pores (FCP) diameters. Classes can be combined into total draining pores (TDP) (0-33 kPa), and total water-storage pores (WSP) (> 33 kPa), as well as into macro (non-capillary) pores (<10 kPa) and soil matrix (capillary) pores (>10 kPa).
2.2. Waterflow in relation to soil capillary pores
The defined ranges of soil pore sizes indicate that the rapidly drainable pores (RDP) radius (r) may range from 15 - 2500 μm (Chen and Wagenet 1992). Within this range, as the radius increases, the Reynolds number (Re, dimensionless) for a water-filled- cylindrical pore (tube) will increase in value, creating a possible change in flow conditions from laminar to non-laminar flow:
where u :is the average flow velocity in the pores tube (LT-1), r is an equivalent radius, and η is the dynamic viscosity of water (10-3 kg m-1s-1 at 20C) and pw is the density of water (103 kg m-3).
Given Re = 1 (i.e., inertia force = viscosity force), and as gravity acceleration, -then from Eq.2a, the criteria radius (rc) between laminar and non laminar flow can be calculated as (Chen and Wagenet, 1922):
By applying the aforementioned values of the parameters in the right term of the Eq.3, then the critical radius of pore is equal to rc= 95μm. As the water-filled, non-capillary pores (RDP) have diameters ≥28.8μm (r≥ 14.4 μm), so it is expected that rc divides the RDP into laminar (Re1 at r = 14.4 - 95 μm) and non-laminar (Re>1 at r> 95μm) water flow. The Poiseuille equation for average velocity (u) of water can also be applied for a pore radius ≤ rc in the form (Amer et al., 2009):
where πr2 is the cross-sectional area of the soil pores. Soil matrix is assumed to behave as a homogeneous soil with a majority of capillary pores (SDP+WHP+FCP), described by classical water flow equations based on Darcy's law, in which:
If the distribution density function of pore sizes, f (r), is defined for pore sizes from rmin . to rmax within a unit area of a capillary bundle, the total water flow capacity Q (L3T-1) at =1, will be:
For capillary pores, water flow Q(r) can be formulated as follows:
where r is the pore radius as applicable to the water adsorption capacity (Wa) at which water does not move, fine capillary pores (FCP), water holding pores (WHP), and slowly drainable pores (SDP), respectively. The pore radius was taken to be as the largest one for each class.
2.3. Infiltration and hydraulic conductivity into non-capillary pores
Under the conditions of non-laminar flow where Re>1, the average flow velocity (u) and then the hydraulic conductivity, K ,can be estimated from the infiltration rate, i and sorptivity, S. The infiltration rate (i) is the volume of water infiltrating through a horizontal unit area of soil surface at any instant (infinitely small period of time) and has unit of LT-1. The cumulative infiltration (I) is the total volume of water that has infiltrated through a unit of horizontal area of soil surface over a given period of time, measured from the beginning of infiltration.. For many soils, a plot of I as a function of time t is described (in depth unit, cm) by the equation of Kostiakov (1932):
where, c and m are constants for a given soil and a given moisture content, respectively. These constants could be obtained if log I is plotted against log t;
At unit time, I = c, where log c = intercept and m = slope of the linear regression of Eq.11.
The instantaneous infiltration rate (i) at any time can be derived by differentiating the expression for I (Eq.10) with respect to time t;
Philip (1957) showed that the cumulative infiltration, I, cm in soil changes with square root of time (t05), with the t in minutes, and it depending on the sorptiv-ity (S) of the soil. Sorptivity is the capacity of a soil to absorb (suck up) water and is influenced by the antecedent water content of the soil. If I is plotted against t05 then a linear relationship is usually found for the first 1 to 3 minutes of infiltration. The slope of the linear relationship,allows to determine S at the unit time. With respect to the Philip infiltration equation, the cumulative infiltration, I, and the sorp-tivity, S, can be calculated at relatively short periods of time as follows:
where α is a constant ranged as 0≤a≤ 1, and Ks is the saturated hydraulic conductivity.
After long periods of times, t at which the instantaneous infiltration rate (i) reached the basic infiltration rate (ib), the steady state of water flow is predominant, and ib is supposed to be proportional to saturated hydraulic conductivity ( i.e. ib = aKs). Amer (2011) proved that the matching factor a (dimensionless) as a function between i and Ks, is:
where, the parameters in the right term of Eq.15 are the same as mentioned above.
In this study, S, after a long period of time, t hr from the beginning infiltration, is denoted as wet sorptivity or steady state sorptivity, Sw (LT-0.5). To estimate Sw, the Eq.14 should be developed using Eq.12 as follows:
At steady state infiltration, i= ib, hence, substituting Eq.16 into Eq.10 produces:
Incorporating Eq.17 into Eq.14, subsequently, we obtain:
During infiltration, the application rate of water to soil surface is often greater than K(θ)RDP + K(θ)h, where K(θ)RDP is the hydraulic conductivity into the rapidly drainable pores and K(θ)h is the matrix unsaturated hydraulic conductivity. If the addition of water continues until the time of steady state infiltration is reached, it is expected that the macro non-capillary pores (Δθ) RDP will eventually be completely saturated with water before the matrix is saturated (Chen and Wagenet 1992). Then the steady state infiltration rate can then be considered as ib= αK(θ)RDP for rapidly drainable pores. In this case, the matching factor (α) is assumed to represents the inverse of the saturation degree at steady state infiltration, i.e. Hence, combing these relationships with respect to Eq. 20, the K(θ)RDP can be estimated as follows :
2.4. Elapsed time for steady state infiltration
The time, t, that must elapse before the instantaneously infiltration rate, i, becomes approximately constant can be expressed in term of the soil property, m. In this case, ib can be defined by differentiating i with respect to time (Eq.12) and then by assuming that stability is reached when the change in i values with time, t does not exceed 10% of i at the previous hour , as follows (Amer, 2011):
By solving for t in the right-hand two terms of the last equation, we find:
2.5. Q(θ) and K(θ) in relation to soil water filled pores
For rapidly drainable pores (RDP), the existence of a range of pore sizes between rmin and rmax (in Eq.8), al lows definition of Q (r) for capillary and non-capillary pores i.e. from rmin = 14.4 μm (= rSDP) up to the critical radius rc (= 95 μm) and from rc up to rmax (≈1500 μm). The r represents the total pores ofsoil at saturation (pF = 0). Substituting Eq.21 into Eq.8, the water flow into RDP can be obtained as follows:
In practices, since the distribution density function, f (r), of pore sizes is difficult to obtain (Chen and Wagenet 1992), a conceptual method of water filled pores, Δθ, is employed to estimate f (r) function whereas the f (r) is usually taken to be similar as the soil water retention (Ψ) function with r corresponding to:where θ is the volumetric soil water content (L3L-3) corresponds to Ψ and r.
Multiplying the first part of Eq.27 by a matching factor Ks/Kc, i.e. the ratio of the measured saturated K to the calculated Kc for total pores, the K(θ) can be rearranged as:
where Δθ is the change in the ratio of volume fraction (L3.L-3) of moisture, going from soil pores occupied by water from the pores of immobile water (soil moisture adsorption capacity, Wa) up to slowly drainable pores (SDP). The Δθ can be expressed as a fraction of the soil bulk volume (as in Eq.28), or as ratio of total pores (Â¿9/9s). When this fraction is substituted for Δθ
in Eqns 28 and 26, the results are:
The critical pore radius (rc) between laminar and turbulent flow is corresponding to the critical moisture content, and rc which can be estimated indirectly from soil moisture retention curve at rc= 95μm. Regardless of the value of rc , the rapid water flow in the non-capillary pores will result in some transfer of water into the soil matrix through the macro-pore wall due to the difference in water potential between the macro-pore and the matrix. The K cutoff point is for the largest r of each certain class, and this was matched with the corresponding Δθ class, and the Wa was subtracted from Δθ of the >1500 kPa class. Larger classes are accumulated of the K's from smaller classes. The value of hydraulic conductivity K(θ) is recognized as being it depends on the nature of the medium and the physical properties;(pwg/η) of the perfuse water.
3. Material and methods
Two soil profiles were chosen to investigate the applicability of the suggested equations. The first soil profile (I) is a calcareous sandy loam soil located at the Nubaria area (northern west Nile Delta). The second (II) is an alluvial clay soil from the farm of the Faculty of Agriculture, Menoufia University, at Shebin Elkom (middle Nile Delta). Disturbed and undisturbed soil samples were taken from each of three depths for each soil profile. The undisturbed soil samples were used to determine bulk density, and saturated hydraulic conductivity by using the constant head method (Klute and Dirksen, 1986). The de-sorption soil water retention curve, for pressure heads h up to about 100 kPa, was determined using undisturbed samples for clay soil (profile II), and a neutron probe and tensiometers were used in situ for calcareous soil (profile I). The determination of soil moisture content (θ), on volume basis, was extended in the laboratory at different h values up to 1500 kPa for the two soil profiles. The data of pore size were obtained using the soil water retention h(θ) curves, and then classified diametrically upon Eq.1. The pore size distribution was expressed as a percent of soil bulk volume (Δθ%) and of total pores volume, (Δθ/θs)x100, (where θs is saturation degree). Disturbed samples were air-dried, gently crushed, sieved through a 2 mm sieve, and used for analysis of desorption at higher pressure heads, CaCO3 content, texture, saturated hydraulic conductivity (Ks), salinity (EC), and sodium adsorption ratio (SAR) These analyses were done following the procedures as mentioned by Sparks et al. (1996) and Dane and Topp (2002). As given in Amer (1993, 2009), we used the water vapour adsorption isotherm method, with applying the BET theory to estimate the moisture adsorption capacity (Wa), which is considered to be immobile water content in a soil. Wa is equal to three layers of adsorbed water (films); Wa = Wm + 2 Wme where, Wm is the mono-adsorbed layer of water vapor on soil particles, and Wme is the external mono-adsorbed layer of water vapor. Some physical and chemical properties of the studied soils are shown in Table 1.
The unsaturated hydraulic conductivity, K(θ), values were calculated by Eq.29 and compared with estimated data (KG ) by van-Genuchten's equation (1980) in calcareous soil:
whereis the volumetric soil moisture content at the moisture suction h, kPa, θr is the residual moisture content, (in this study θr = PbWa), cþ is total porosity, and n is a constant depending on the fitting of the sigmoid curve.
Field infiltration rate was obtained with the double ring infiltrometer (Reynolds et al., 2002). The sorptivity and infiltration parameters were done for both locations of soil profiles. Measured intake rates for individual infiltration runs were obtained at 2 to 10-minute intervals for the triplicates measurements of the locations. Infiltration rates (i in cm/h), as fitted to the Kostiakov equation (1932) were measured in the experimental field.
4. Results and discussion
4.1. Capillary and non-capillary pores Table 2 shows the pore size distribution; RDP, SDP, WHP, and FCP as derived from the soil moisture retention data of the studied soils applying the Eq.1. The non-capillary rapidly drainable pores (RDP) and ratio of air to water (A/W) indicated that the total drained pores, as compared with the capillary and water storage pores, were higher in the sandy loam soil (profile I) than in the clay soil (Profile II). This is due to the relatively coarse texture and high CaCO3 content in the profile I. The maximum value of A/W ratio was observed in the surface depth of the sandy loam soil, but it was occurred in the subsurface depth of the clay soil. As expected, the numbers of capillary pores were as much as high in the clay soil (fine textured profile) if compared to the sandy loam soil profile.
4.2. Infiltration and sorptivity
Experimentally, the typical trends in cumulative infiltration, I (cm) versus time t (minute) and infiltration rate, i (cm/h) are illustrated by the empirical power functions, I= 1.44 t 0.65, i = 56 t -0.35, and I = 0.6 t 0.Â 53, i = 19.1 t -0.47 in sandy loam (profile I) and clay (profile II) soils, respectively. At the onset of wetting, the moisture gradient was the greatest, hence the more rapid infiltration. The infiltration rate i (LT-1) slowed down gradually with time t and reached the steady state (basic infiltration rate, ib ) after 3.5 and 4.7 hours (Eq.23) from the beginning of infiltration for I and II soil profiles, respectively. At the beginning of infiltration (un-saturation conditions), the sorptivity, S, was calculated using I as a function of time tm where m was adjusted to 0.5. The S values were found as 1.324cm/ min0.5 and 0.618 cm/min0.5 for the profiles I and II, respectively. After a long time of infiltration, the sorptivity may be called steady state sorptivity or wet sorptivity (Sw) (Amer 2011). Using Eq.19 to calculate Sw , it was observed that at steady state infiltration the sorptivity values were 1.124 and 0.334 cm / min0.5, 1.Â e., decreased by 15.1% and by 45.95% for I and II soil profiles, respectively. It means that a dry soil (un-saturated) typically has a much greater sorptivity than a wet soil. Though a soil with larger (non-capillary) pores has a lower sorptivity than a soil with smaller (capillary) pores (Hallett 2008), S was higher in sandy loam (calcareous) soil than in clay soil. This result attributed to the abundance of the CaCO3 fraction which has a great ability to suck up water in such calcareous soil (Ghazy, 1993).
4.3. Waterflow rate
The flow of water under saturated conditions is a function of two opposing factors, namely, the effective pore size and the continuity of water films in the conducting pores. At a given suction, the decrease in pore size facilitates the flow of water by increasing the continuity of the water films in the water conducting pores (Amer 2003). Simultaneously, however, it slows down the flow of water by the increased resistance in pores because from Poiseuille's law and the present derived equations, the water flow rate per unit cross-section area of pores is proportional to the square of the pore radius (r2). The ultimate effect of particle-size/pore size (capillary and non-capillary) distribution on unsaturated water flow is thus the resultant of these two opposing phenomena (Mundra et al., 1989). However, as the values of K(θ) increased with an increase in water content (θ), the water flow rate Q(θ) must obey the same trend of K(θ) for the same radius, since; Q(θ) = πr2 K(θ).
4.4. Unsaturated hydraulic conductivity
Data presented in Table 3 show the values of unsatu-rated hydraulic conductivity, K(θ), as calculated by Eq.29 at different soil pore size classes. It was evident that the values of hydraulic conductivity for fine capillary pores, (K(θ)FCP), were very low in the soil matrix whereas they fluctuated between (3.26-4.17) x10-9 and (0.38-2.30)x10-8 cm/s in the I and II soil profiles, respectively. As expected the values of K(θ) increased with increase in capillary and non-capillary pore sizes from FCP to SDP and up to RDP (Fig.1). The larger the pore size the greater is the hydraulic conductivity (Ahuja et al., 1989). The values of K(Δθ)RDP in the range of non-capillary pores remained smaller than those for the range of soil matrix (capillary) pores (K(Δθ)h), in clay soil (Profile II). On the contrary, the K(Δθ)RDP and K(Δθ)SDP values were higher in sandy loam soil (profile I) than that in soil matrix pores. This is due to the relatively large pores and high RDP fraction percent of all depths of the sandy calcareous soil (profile I) as compared to the clay soil (Table 2). The opposite was observed for WHP and FCP where K(Δθ)h values were higher for all depths of the clay soil than in the sandy loam soil.
In comparison with the van-Genuchten's model, the predicted values of K(O) by the suggested equation (29) were in consistent with those calculated by van-Genuchten's model (Eq.31) (K(θ)g ) for the water holding pores values, WHP. The calculated K(θ)g values were obtained only for sandy calcareous soil (profile I), whereas the values of the constant n (in Eq.31) were 0.1837, 0.1632, and 0.1603 for the calcareous soil depths; 0-20, 20-40, and 40-60 cm respectively.
However, the values of K(θ)g were relatively inconsistent for the other pore size classes. These results may be due to the fact that the calculation of the n constant was done only over the water retention range 0-100 kPa of the soil moisture characteristic curve. This range includes the values of air-entry point and (Δθ)WHP. In general, the values of the predicted K(θ) were less by about 3 times than those of K(θ)g / Ks for calcareous soil (Table 3).
The characterization of water flow properties and un-saturated hydraulic conductivity in the root zone of soils has become an ever-increasing problem in Soil Physics. One of the main problems is the understanding of flow in a constantly changing transforming porous media, where the soil pores and water pathways evolve spatially and temporally, for example due to interaction with the flow, aggregates, compaction, destruction of permeability...etc. To overcome this problem, new equations are presumed to predict the water flow Q(θ) and hydraulic conductivity K(θ) into capillary and non-capillary pore size classes in coarse (sandy loam) and fine (clay) textured soils. Then the total hydraulic conductivity of the studied soils is predicted as a summation of hydraulic conductivity for each pore size class. The pore size classes are explained in this work as rapidly drainable pores (RDP), slowly drainable pores (SDP), water holding pores (WHP), and fine capillary pores (FCP). They were derived upon the soil moisture retention data of the studied sandy loam and clay soils. Also, the progress in this work is how to solve the problems regarding the upper limit of the rapidly drainable pores (RDP) diameter, and the expected predominant of non-laminar (turbulent) flow at saturation condition into the large pore sizes. To solve this problem the conceptual of the critical radius that found to be rc = 95 μm was provided in the present study. Upon the rc value, the RDP range was divided into two parts; (1) ranges from the upper limit of slowly drainable pores (SDP, 14.4 μm ) up to the critical limit or radius, rc (= 95 μm), (2) ranges from rc (= 95 μm) up to the maximum limit of pore size at full saturation. In this respect a new model, depending on water sorptivity (S) at steady state infiltration, was suggested to determine the water flow and hydraulic conductivity into the RDP upon the range of non-laminar (turbulent) water flow, i.e. from rc = 95 Â°m up to the maximum size of RDP ( rmax ). To obtain the effective water content which contribute to Q(θ) and K(θ) values, the immobile adsorbed water (soil moisture adsorption capacity, Wa) was reducted from the water content at fine capillary pores (θ)FCP The values of K(θ) calculated by the assumed equations were in agreement with those calculated by van Genuchten's model (K(θ)G) for WHP only. However, the calculated values of K(θ) were within the range of the commonly measured values for the studied soils, indicating the applicability of the assumed equations for predicting the unsaturated hydraulic conductivity in both coarse and fine textured soils.
I would like to express my gratitude to William De-Tar, retired Agric. Engr., USD of Agric.-ARS, US Cotton Res. Station, 17053 Shafter Avenue, Shafter, CA 93263, for correcting the English language of the manuscript. His assistance is gratefully acknowledged. I am also grateful to the anonymous reviewers for their helpful comments.
Ahuja, L.R. Cassel, D.K., Bruce, R.R., Bames, B.B. 1989. Evaluation of spatial distribution of hydraulic conductivity using effective porosity data. Soil Science 148:404-411. [ Links ]
Amer,A.M. 1993. Surface area measurements as related to water vapor adsorption in arid soils of Egypt. Proc. Int. Conf. Desert Develop.,25-30 July, Mexico City, Mexico, pp.619-627. [ Links ]
Amer, A.M. 2002. Drainable and water-filled pores as related to water storage and conductivity in agricultural soils of the Nile Delta. Verh. Internal Verein. Limnology, Vol.28, part 4, p.1912-1919 Stuttgart-Germany, 28th SIL Congress, 4-10 February, 2001, Manash Univ., Melbourne, Australia. [ Links ]
Amer, A.M. 2003. Soil hydrophysics. Part1. El-Dar Al-Arabia for Publish & Distrib., Cairo, Egypt. 452pp. [ Links ]
Amer, A.M. 2009. Moisture adsorption capacity and surface area as deduced from vapour pressure isotherms in relation to hygroscopic water of soils. Biologia 64, 3: 516-521. [ Links ]
Amer, A.M., Logsdon, S.D., Davis, D. 2009. Prediction of hydraulic conductivity in un- saturated soils. Soil Science 174, 9: 508- 515. [ Links ]
Amer, A. M. 2011. Prediction of hydraulic conductivity and sorptivity in soils at steady state infiltration. Archives of Agronomy and Soil Science. Online http://www.tandfonline.com/
Baver, L.D., Gardner, W.R. 1972. Soil physics. 4th ed. John Wiley & Sons, New York, 498pp. [ Links ]
Beven, K., Germann, P.1982. Macropores and water flow in soils. Water Resour. Res. 18: 1311-1325. [ Links ]
Chen, C., Wagenet, R.J. 1992. Simulation of water and chemicals in macropore soils. Part 1. Representation of the equivalent macropore influence and its effect on soil water flow. J. Hydrlogy 130: 105-126. [ Links ]
Dane, J.H., Topp, G.C. (Eds.). 2002. Methods of soil analysis, Part 4, Physical methods. SSSA, Madison, WI. 1692 pp. [ Links ]
Ghazy, A.E. 1993. Pore size distribution and moisture behavior in highly calcareous soils. Egypt. J. Soil Sci. 22, 57-67. [ Links ]
Hallett, P.D. 2008. A brief overview of the causes, impacts and amelioration of soil water repellency - a review. Soil & Water Res. 3, (Sp.Issue 1): S21 - S29. [ Links ]
Hillel, D. 1980. Fundamentals of soil physics. Academic press, New York. [ Links ]
Klute, A., Dirksen, C. 1986. Hydraulic conductivity and diffusivity: Laboratory methods. In Methods of soil analysis. Part 1, Agron. Monogr. 9, 2nd Ed, A. Klute (ed.). ASA and SSSA, Madison, WI, pp.687-734. [ Links ]
Kostiakov, A.N. 1932. On the dynamics of coefficient of water-percolation in soils and on necessity for studying it from a dynamic point of view for purposes of amelioration. Trans. 6th Comm. Intern. Soil Sci., Russia, Part A, pp 17-21. [ Links ]
Marshall, T. J. 1956. Relation between water and soil. Technical Communication No. 50, Commonwealth Bureau of Soils, Farenham Royal, Bucks, England. [ Links ]
Mualem, Y. 1976. A new model predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 12,513-522. [ Links ]
Mundra, M.C., Raj Pal, Siyag, R.S., Poonia, S.R. 1989. Unsaturated water transmission characteristics of soils in relation to texture, aggregate size and initial moisture content. J. Agric. Sci. Camb., 112, 199-204. [ Links ]
Philip, J.R. 1957. The theory of infiltration:1.The infiltration equation and its solution. Soil Science 83, 345-357. [ Links ]
Philip, J.R., 1969. Theory of infiltration. Adv. Hydro-sci. 5, 215-290. [ Links ]
Reynolds, W.D., Elrick, D.E., Young, E.G. 2002. 3.4.2.a. Single - ring and double - ring or concentric - ring infiltrometers. p. 821-826. In (J.H. Dane and G.C. Topp, Eds. Methods of soil analysis Part 4 Physical methods. SSSA, Madison. WI. [ Links ]
Sparks, D.L. (ed.). 1996. Methods of soil analysis, Part 3, Chemical methods. ASA, SSSA, Madison, WI. [ Links ]
Swartzendruber, D., Young, E.G. 1974. A comparison of physically-based infiltration equation. Trans. ASAE 12:822-828. [ Links ]
Van Genuchten, M. Th. 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44, 892-898. [ Links ]
Zhang, R. 1997. Determination of soil sorptivity and hydraulic conductivity from the disk infiltometer Soil Sci. Soc. Am. J. 61,1024-1030. [ Links ]