# 2:1.2 Runoff Volume: Green & Ampt Infiltration Method The Green & Ampt equation was developed to predict infiltration assuming excess water at the surface at all times (Green & Ampt, 1911). The equation assumes that the soil profile is homogenous and antecedent moisture is uniformly distributed in the profile. As water infiltrates into the soil, the model assumes the soil above the wetting front is completely saturated and there is a sharp break in moisture content at the wetting front. Figure 2:1-2 graphically illustrates the difference between the moisture distribution with depth modeled by the Green & Ampt equation and what occurs in reality. ![](/files/89YJhx1FOWpC8GOjjWlF) Mein and Larson (1973) developed a methodology for determining ponding time with infiltration using the Green & Ampt equation. The Green-Ampt Mein-Larson excess rainfall method was incorporated into SWAT+ to provide an alternative option for determining surface runoff. This method requires sub-daily precipitation data supplied by the user. The Green-Ampt Mein-Larson infiltration rate is defined as: $$f\_{inf,t}=K\_e\*(1+\frac{\Psi\_{wf}\*\Delta\theta\_v}{F\_{inf,t}})$$ 2:1.2.1 where $$f\_{inf}$$ is the infiltration rate at time $$t$$ (mm/hr), $$K\_e$$ is the effective hydraulic conductivity (mm/hr), $$\Psi\_{wf}$$ is the wetting front matric potential (mm), $$\Delta\theta\_v$$ is the change in volumetric moisture content across the wetting front (mm/mm) and $$F\_{inf}$$ is the cumulative infiltration at time $$t$$ (mm H$$\_2$$O). When the rainfall intensity is less than the infiltration rate, all the rainfall will infiltrate during the time period and the cumulative infiltration for that time period is calculated: $$F\_{inf,t}=F\_{inf,t-1}+R\_{\Delta t}$$ 2:1.2.2 where $$F\_{inf,t}$$ is the cumulative infiltration for a given time step (mm H$$*2$$O), $$F*{inf,t-1}$$ is the cumulative infiltration for the previous time step(mm H$$*2$$O), and $$R*{\Delta t}$$ is the amount of rain falling during the time step (mm H$$\_2$$O). The infiltration rate defined by equation 2:1.2.1 is a function of the infiltrated volume, which in turn is a function of the infiltration rates in previous time steps. To avoid numerical errors over long time steps, $$f\_{inf}$$ is replaced by $$dF\_{inf}/dt$$ in equation 2:1.2.1 and integrated to obtain $$F\_{inf,t}=F\_{inf,t-1}+K\_e\*\Delta t+ \Psi\_{wf}*\Delta\theta\_v*ln\[\frac{F\_{inf,t}+\Psi\_{wf}*\Delta\theta\_v}{F\_{inf,t-1}+\Psi\_{wf}*\Delta\theta\_v}]$$ 2:1.2.3 Equation 2:1.2.3 must be solved iteratively for $$F\_{inf,t}$$, the cumulative infiltration at the end of the time step. A successive substitution technique is used. The Green-Ampt effective hydraulic conductivity parameter, $$K\_e$$, is approximately equivalent to one-half the saturated hydraulic conductivity of the soil, $$K\_{sat}$$ (Bouwer, 1969). Nearing et al. (1996) developed an equation to calculate the effective hydraulic conductivity as a function of saturated hydraulic conductivity and curve number. This equation incorporates land cover impacts into the calculated effective hydraulic conductivity. The equation for effective hydraulic conductivity is: $$K\_e=\frac{56.82*K\_{sat}^{0.286}}{1+0.051*exp(0.062\*CN)}-2$$ 2:1.2.4 where $$K\_e$$ is the effective hydraulic conductivity (mm/hr), $$K\_{sat}$$ is the saturated hydraulic conductivity (mm/hr), and $$CN$$ is the curve number. Wetting front matric potential, $$\Psi\_{wf}$$, is calculated as a function of porosity, percent sand and percent clay (Rawls and Brakensiek, 1985): $$\Psi\_{wf}=10*exp\[6.5309-7.32561*\phi\_{soil}+0.001583*m\_c^2+3.809479*\phi\_{soil}^2+0.000344*m\_s*m\_c-0.049837*m\_s*\phi\_{soil}+0.001608*m\_s^2*\phi\_{soil}^2+0.001602*m\_c^2*\phi\_{soil}^2-0.0000136*m\_s^2*m\_c-0.003479*m\_c^2*\phi\_{soil}-0.000799*m\_s^2*\phi\_{soil}]$$ 2:1.2.5 where $$\phi\_{soil}$$ is the porosity of the soil (mm/mm), $$m\_c$$ is the percent clay content, and $$m\_s$$ is the percent sand content. The change in volumetric moisture content across the wetting front is calculated at the beginning of each day as: $$\Delta\theta\_v=(1-\frac{SW}{FC})*(0.95*\phi\_{soil})$$ 2:1.2.6 where $$\Delta\theta\_v$$ is the change in volumetric moisture content across the wetting front (mm/mm), $$SW$$ is the soil water content of the entire profile excluding the amount of water held in the profile at wilting point (mm H$$\_2$$O), $$FC$$ is the amount of water in the soil profile at field capacity (mm H$$*2$$O), and $$\phi*{soil}$$ is the porosity of the soil (mm/mm). If a rainfall event is in progress at midnight, $$\Delta\theta\_v$$ is then calculated: $$\Delta\theta\_v=0.001\*(0.95\*\phi\_{soil})$$ 2:1.2.7 For each time step, SWAT+ calculates the amount of water entering the soil. The water that does not infiltrate into the soil becomes surface runoff. Table 2:1-2: SWAT+ input variables that pertain to Green & Ampt infiltration calculations. | Definition | Source Name | Input Name | Input File | | ------------------------------------------------------------------------------------- | ----------- | ------------- | ---------- | | Rainfall, runoff, routing option. | | IEVENT | .bsn | | Length of time step (min): $$\Delta t$$=IDT/60 | | IDT | file.cio | | $$R\_{\Delta t}$$: Precipitation during time step (mm H$$\_2$$O) | | PRECIPITATION | .pcp | | $$K\_{sat}$$: Saturated hydraulic conductivity of first layer (mm/hr) | | SOL\_K | .sol | | $$CN$$: Moisture condition II curve number | | CN2 | .mgt | | $$CN$$: Moisture condition II curve number | | CNOP | .mgt | | $$\rho\_b$$: Moist bulk density (Mg/$$m^3$$): $$\Psi\_{soil}$$=1 - $$\rho\_b$$ / 2.65 | | SOL\_BD | .sol | | $$m\_c$$: % clay content | | CLAY | .sol | | $$m\_s$$: % sand content | | SAND | .sol | --- # Agent Instructions: Querying This Documentation If you need additional information that is not directly available in this page, you can query the documentation dynamically by asking a question. Perform an HTTP GET request on the current page URL with the `ask` query parameter: ``` GET https://swatplus.gitbook.io/io-docs/theoretical-documentation/section-2-hydrology/chapter-2-1-surface-runoff/2-1.2-runoff-volume-green-and-ampt-infiltration-method.md?ask= ``` The question should be specific, self-contained, and written in natural language. The response will contain a direct answer to the question and relevant excerpts and sources from the documentation. Use this mechanism when the answer is not explicitly present in the current page, you need clarification or additional context, or you want to retrieve related documentation sections.