SWMM 2008 - Stormwater Management Model
SWMM 5, Watersheds, Water Quality, Technology, Hydrology, Hydraulics - Watershed
Parameters in the Horton Equation 3-6 are
f∞ = minimum or ultimate value of fp (infiltration capacity), in./hr [mm/hr],
fo = maximum or initial value of fp, in./hr [mm/hr], and
α = decay coefficient, hr-1.
The equation may be derived theoretically (Eagleson, 1970) but is almost always applied by curve-fitting to infiltrometer data. A regeneration parameter, R, is also required to compute a regeneration rate coefficient, ad = R a (Equation 3-20).
Although the Horton infiltration equation is probably the best-known of the several infiltration equations available, there is little to help the user select values of parameters fo and α for a particular application. (Fortunately, some guidance can be found for the value of fc.) Since the actual values of fo and α (and often fc) depend on the soil, vegetation, and initial moisture content, ideally these parameters should be estimated using results from field infiltrometer tests for a number of sites of the watershed and for a number of antecedent wetness conditions. An example of Horton parameters for Georgia soils is given in Table 3-8 (Rawls et al., 1976). Horton’s (1940) estimates are shown in Table 3-9. Skaggs and Khaleel (1982) provide Horton-type decay curves on the basis of theoretical estimates.
Table 3-8. Horton parameters for selected Georgia soils (Rawls et al., 1976).
|
Soil Type |
f∞ in./hr |
fo in./hr |
α 1/hr |
|
Alpha loamy sand |
1.40 |
19.0 |
38.29 |
|
Carnegie sandy loam |
1.77 |
14.77 |
19.64 |
|
Cowarts loamy sand |
1.95 |
15.28 |
10.65 |
|
Dothan loamy sand |
2.63 |
3.47 |
1.40 |
|
Fuquay pebbly loamy sand |
2.42 |
6.24 |
4.70 |
|
Leefield loamy sand |
1.73 |
11.34 |
7.70 |
|
Robersdale loamy sand |
1.18 |
12.41 |
21.75 |
|
Stilson loamy sand |
1.55 |
8.11 |
6.55 |
|
Tooup sand |
1.80 |
23.01 |
32.71 |
Table 3-9. Horton equation parameters provided by Horton (1940).
|
Soil and Cover |
f∞ in./hr |
fo in./hr |
α 1/hr |
|
Standard agricultural (bare) |
0.24 – 8.9 |
11.4 |
96 |
|
Standard agricultural (turfed) |
8.2 – 11.8 |
36.7 |
48 |
|
Peat |
0.82 – 11.8 |
13.3 |
108 |
|
Fine sandy clay (bare) |
0.82 – 1.0 |
8.6 |
120 |
|
Fine sandy clay (turfed) |
4.1 – 1.2 |
27.4 |
84 |
If it is not possible to use field data to find estimates of fo, f∞, and α for each subcatchment, the following guidelines should be helpful. The NRCS, formerly the Soil Conservation Service or SCS, has classified most soils into Hydrologic Soil Groups, A, B, C, and D, dependening on their limiting infiltration capacities, fc. (Well drained, sandy soils are “A”; poorly drained, clayey soils are “D.”) A listing of the groupings for more than 4000 soil types can be found in the SCS Hydrology Handbook (1972, pp. 7.6-7.26); a similar listing is also given in Chow’s Handbook of Applied Hydrology (Ogrosky and Mockus, 1964, pp. 21.12-21.25), but the former reference also gives alternative groupings for some soil types depending on the degree of drainage of the subsoil. NRCS soil classification data are also available on CD-ROMs from vendors such as Earth Info (Boulder, CO). The soil type itself may be found in the U.S. from county NRCS Soil Survey maps.
The best source of information about a particular soil type is a publication entitled “Soil Survey Interpretations” available from a local NRCS office in the U.S. Information on the soil profile, the soil properties, its suitability for a variety of uses, its erosion and crop yield potential, and other data are included on the sheet for undisturbed soil samples. A copy of the listing for Conestoga silt loam is shown in Figure 3-37. Parameter f∞is essentially equal to the saturated hydraulic conductivity, Ks, which is called “permeability” on the soil survey interpretation sheet. For Conestoga Silt Loam, a range of 0.63-2.0 in./hr (16-51 mm/hr) is shown. Similar data are presented in aggregate for most counties in the United States in the Soil Survey for that county, again, obtainable from local NRCS offices or from a local Soil and Water Conservation District.
Alternatively, values for f∞ according to Musgrave (1955) are given in Table 3-10. To help select a value within the range given for each soil group, the user should consider the texture of the layer of least hydraulic conductivity in the profile. Depending on whether that layer is sand, loam, or clay, the fc value should be chosen near the top, middle, and bottom of the range respectively. For example, the data sheet for Conestoga silt loam identifies it as being in Hydrologic Soil Group B, which puts the estimate of fc into the range of 0.15-0.30 in./hr (3.8-7.6 mm/hr), much lower than the Ks value discussed above. Examination of the texture of the layers in the soil profile indicates that they are silty in nature, suggesting that the estimate of the f∞ value should be in the low end of the range, say 0.15-0.20 in./hr (3.8-5.1 mm/hr). A sensitivity test on the f∞ value will indicate the importance of this parameter to the overall result.
Table 3-10. Values of f∞ for Hydrologic Soil Groups (Musgrave, 1955). (Hydrologic Soil Group is defined in Table 3-15 and discussed in Section 3.8.8.5.)
|
Hydrologic Soil Group |
f∞ (in./hr) |
||
|
A |
0.45 |
- |
0.30 |
|
B |
0.30 |
- |
0.15 |
|
C |
0.15 |
- |
0.05 |
|
D |
0.05 |
- |
0 |
Caution should be used in applying values from Table 3-10 to sandy soils (group A) since reported Ks values are often much higher. For instance, sandy soils in Florida have Ks values from 7 to 18 in./hr (180-450 mm/hr) (Carlisle et al., 1981). Unless the water table rises to the surface, ultimate infiltration capacity will be very high, and rainfall rates will almost always be less than f∞, leading to little or no overland flow from such soils.
When water ponds on the surface, the limiting value of surface infiltration is the saturated hydraulic conductivity, Ks. Thus,
fc∞ ≈ Ks (3-133)
Good, generalized estimates for Ks will be discussed in conjunction with Green-Ampt parameter estimates in the next section and likely are the best source of values for f∞ in the absence of site-specific data.
For any field infiltration test the rate of decrease (or “decay”) of infiltration capacity, α, from the initial value, fo, depends on the initial moisture content. Thus the α-value determined for the same soil will vary from test to test.
Figure 4-19. Soil Conservation Service Soil Survey Interpretation for Conestoga silt loam (found near Lancaster, PA).
It is postulated here that, if fo is always specified in relation to a particular soil moisture condition (e.g., dry), and for moisture contents other than this the time scale is changed accordingly (i.e., time “zero” is adjusted to correspond with the constant fo), then α can be considered a constant for the soil independent of initial moisture content. Put another way, this means that infiltration curves for the same soil, but different antecedent conditions, can be made coincident if they are moved along the time axis. Butler (1957) makes a similar assumption.
Values of α found in the literature (Linsley et al., 1975; Overton and Meadows, 1976; Wanielista, 1978; ASCE, 1996) range from 0.67 to 120 hr-1. Nevertheless most of the values cited appear to be in the range 3-6 hr-1 (0.00083-0.00167 sec-1). [MSOffice1] The evidence is not clear as to whether there is any relationship between soil texture and the α value although several published curves seem to indicate a lower value for sandy soils. If no field data are available, an estimate of 0.00115 sec-1 (4.14 hr-1) could be used. Use of such an estimate implies that, under ponded conditions, the infiltration capacity will fall 98 percent of the way towards its minimum value in the first hour, a not uncommon observation. Rates of decay of infiltration for several values of α are shown in Table 3-11.
Table 3-11. Rate of decay of infiltration capacity for different values of α
|
α-value hr-1, rounded, and (sec-1) |
Percent of decline of infiltration capacity towards limiting value fc after 1 hour |
|
|
2 |
(0.00056) |
76 |
|
3 |
(0.00083) |
95 |
|
4 |
(0.00115) |
98 |
|
5 |
(0.00139) |
99 |
The initial infiltration capacity, fo depends primarily on soil type, initial moisture content, and surface vegetation conditions. For example, Linsley et al. (1982) present data that show, for a sandy loam soil, a 60 to 70 percent reduction in the fo value due to wet initial conditions. They also show that lower fo values apply for a loam soil than for a sandy loam soil. As to the effect of vegetation, Jens and McPherson (1964, pp. 20.20-20.38) list data that show that dense grass vegetation nearly doubles the infiltration capacities measured for bare soil surfaces.
For the assumption to hold that the decay coefficient α is independent of initial moisture content, fo must be specified for the dry soil condition. When SWMM is run in continuous mode, the program automatically calculates the fo value applicable for wetter conditions as part of the moisture accounting routine. However, for any single-event simulation, the user must specify the fo value for the storm in question, which may be less than the value for dry soil conditions. For this reason, it is often useful to let the program “warm up” by simulating 7 – 14 days prior to the event in question, to let the model establish its own initial conditions.
Published values of fo vary depending on the soil, moisture, and vegetation conditions for the particular test measurement. The fo values listed in Table 3-12 can be used as a rough guide. Interpolation between the values may be required.
Table 3-12. Representative values for fo.
|
A. DRY soils (with little or no vegetation): Sandy soils: 5 in./hr Loam soils: 3 in./hr Clay soils: 1 in./hr |
|
B. DRY soils (with dense vegetation): Multiply values given in A by 2 (after Jens and McPherson, 1964) |
|
C. MOIST soils (change from dry fo value required for single event simulation only): Soils which have drained but not dried out (i.e., field capacity): divide values from A and B by 3 Soils close to saturation: Choose value close to fc value. Soils which have partially dried out: divide values from A and B by 1.5-2.5. |
For continuous simulation, infiltration capacity will be regenerated (recovered) during dry weather according to Equation 3-20. The value of R is probably << 1.0, (implying a “longer” drying curve than wetting curve, Figure 3-7).
On well-drained porous soils (e.g., medium to coarse sands), recovery of infiltration capacity is quite rapid and could well be complete in a couple of days. For heavier soils, the recovery rate is likely to be slower, say 7 to 14 days. The choice of the value can also be related to the interval between a heavy storm and wilting of vegetation. The value of αd is then,
αd = 0.02/D (3-134)
where
αd = R α = recovery curve decay coefficient, day-1, and
D = number of days required for the soil to dry out (recover).
The factor of 0.02 in Equation 3-134 assumes 98 percent recovery of infiltration capacity (i.e., e-0.02 = 0.98). The value of R may then be calculated from Equation 3-20. For example, for α = 4.14 day-1 and drying times of 3, 7 and 14 days, values of R are 1.61 x 10-3, 6.9 x 10-4 and 3.45 x 10-4, respectively.
Last updated by Robert E Dickinson Apr 21.
Robert E Dickinson
created this social network on Ning.

© 2008 Created by Robert E Dickinson on Ning. Create your own social network
Badges | Report An Issue | SWMM 5 Blog | Privacy | Terms of Service