Infiltration is a component of the general mass balance hydrologic budget. There are several ways to estimate the volume and water infiltration rate into the soil. The rigorous standard that fully couples groundwater to surface water through a non-homogeneous soil is the numerical solution of
Richards' equation. A newer method that allows 1-D groundwater and surface water coupling in homogeneous soil layers and that is related to the Richards equation is the
Finite water-content vadose zone flow method solution of the
Soil Moisture Velocity Equation. In the case of uniform initial soil water content and deep, well-drained soil, some excellent approximate methods exist to solve the infiltration flux for a single rainfall event. Among these are the Green and Ampt (1911) method, Parlange et al. (1982). Beyond these methods, there are a host of empirical methods such as SCS method, Horton's method, etc., that are little more than curve fitting exercises.
General hydrologic budget The general hydrologic budget, with all the components, with respect to infiltration
F. Given all the other variables and infiltration is the only unknown, simple algebra solves the infiltration question. :F=B_I+P-E-T-ET-S-I_A-R-B_O where :
F is infiltration, which can be measured as a volume or length; :B_I is the boundary input, which is essentially the output watershed from adjacent, directly connected impervious areas; :B_O is the boundary output, which is also related to surface runoff,
R, depending on where one chooses to define the exit point or points for the boundary output; :
P is
precipitation; :
E is
evaporation; :
T is
transpiration; :
ET is
evapotranspiration; :
S is the storage through either
retention or
detention areas; :I_A is the initial abstraction, which is the short-term surface storage such as puddles or even possibly
detention ponds depending on size; :
R is
surface runoff. The only note on this method is one must be wise about which variables to use and which to omit, for doubles can easily be encountered. An easy example of double counting variables is when the evaporation,
E, and the transpiration,
T, are placed in the equation as well as the evapotranspiration,
ET.
ET has included in it
T as well as a portion of
E. Interception also needs to be accounted for, not just raw precipitation.
Richards' equation (1931) The standard rigorous approach for calculating infiltration into soils is
Richards' equation, which is a
partial differential equation with very nonlinear coefficients. The Richards equation is computationally expensive, not guaranteed to converge, and sometimes has difficulty with mass conservation.
Finite water-content vadose zone flow method This method approximates Richards' (1931) partial differential equation that de-emphasizes soil water diffusion. This was established by comparing the solution of the advection-like term of the
Soil Moisture Velocity Equation and comparing against exact analytical solutions of infiltration using special forms of the soil constitutive relations. Results showed that this approximation does not affect the calculated infiltration flux because the diffusive flux is small and that the
finite water-content vadose zone flow method is a valid solution of the equation is a set of three
ordinary differential equations, is guaranteed to converge and to conserve mass. It requires the assumption that the flow occurs in the vertical direction only (1-dimensional) and that the soil is uniform within layers.
Green and Ampt The name was derived from two men: Green and Ampt. The Green-Ampt method of infiltration estimation accounts for many variables that other methods, such as Darcy's law, do not. It is a function of the soil suction head, porosity, hydraulic conductivity, and time. :\int_0^{F(t)} {F\over F+\psi\,\Delta\theta}\, dF = \int_0^t K\,dt where :{\psi} is wetting front soil suction head (L); :\theta is
water content (-); :K is
hydraulic conductivity (L/T); :F(t) is the cumulative depth of infiltration (L). Once integrated, one can easily choose to solve for either volume of infiltration or instantaneous infiltration rate: :F(t)=Kt+\psi \, \Delta\theta \ln \left[1+{F(t)\over \psi \, \Delta\theta}\right]. Using this model one can find the volume easily by solving for F(t). However, the variable being solved for is in the equation itself so when solving for this one must set the variable in question to converge on zero, or another appropriate constant. A good first guess for F is the larger value of either Kt or \sqrt {2\psi \, \Delta\theta Kt}. These values can be obtained by solving the model with a log replaced with its Taylor-Expansion around one, of the zeroth and second order respectively. The only note on using this formula is that one must assume that h_0, the water head or the depth of ponded water above the surface, is negligible. Using the infiltration volume from this equation one may then substitute F into the corresponding infiltration rate equation below to find the instantaneous infiltration rate at the time, t, F was measured. :f(t)=K\left[{\psi \, \Delta\theta\over F(t)}+1\right].
Horton's equation Named after the same
Robert E. Horton mentioned above, Horton's equation :f(t) = akt^{a-1}+f_0\! in integrated form, the cumulative volume is expressed as: :F(t) = kt^{a}+f_0t\! Where :f_0 approximates but does not necessarily equate to the final infiltration rate of the soil.
Darcy's law This method used for infiltration is using a simplified version of
Darcy's law. :{f} Infiltration rate f (mm hour−1)) :K is the
hydraulic conductivity (mm hour−1)); :L is the vague total depth of subsurface ground in question (mm). This vague definition explains why this method should be avoided. :{S_f} is wetting front soil suction head ({-\psi}) = ({-\psi_f}) (mm) :h_0 is the depth of ponded water above the ground surface (mm); ==See also==