Deficit and Constant Loss Infiltration Method and Parameter Estimation

Deficit and Constant Loss Infiltration Method and Parameter Estimation

The Deficit and Constant Loss Infiltration Method can be implemented in HEC-HMS, HEC-RAS (2D rain-on-grid), and InfoWorks to account for precipitation loss through infiltration and percolation. In order to facilitate the discussion, some key terms are explained below.

  • Porosity or Total Porosity: the fraction of the total soil volume that is taken up by the void space (Unit: dimensionless, V-void/V-bulk soil).
  • Soil Water Content or Moisture Content: the fraction of the total soil volume is taken up by water. The corresponding value at the commencement of a simulation is called initial water content (Unit: dimensionless, V-water/V-bulk soil). The saturated water content can be estimated as total porosity.
  • Field Capacity or Water Holding Capacity: the water remaining in a soil after it has been thoroughly saturated and then allowed to drain freely by gravity (taking one to several days). Usually it is measured as water content retained in soil at −0.33 bar of hydraulic head or suction pressure (Unit: dimensionless, V-water/V-bulk soil).
  • Wilting Point:  the soil water content at which a plant can not absorb any more water and it starts to wilt since the water is held so tightly by the soil particles. Usually it is measured as water content at a soil matric potential of −15 bar (Unit: dimensionless, V-water/V-bulk soil).
  • Available Water Content: water held in soil between its field capacity and wilting point (Unit: dimensionless, V-water/V-bulk soil) – available water content is usually used for irrigation and agriculture.
  • Residual Soil Water Content or Residual Moisture Content: the volumetric water content of a soil where a further increase in negative pore-water pressure does not produce significant changes in water content. It represents the remaining water content after a saturated soil is allowed to drain thoroughly for an extended period of time (Unit: dimensionless, V-water/V-bulk soil).
  • Soil Water Deficit or Soil Moisture Deficit: the amount of water that is required to saturate the soil. The corresponding value at the commencement of a simulation is called initial deficit (Unit: inch).
  • Effective Porosity: usually the soil water content can not be zero: a realistic minimum soil water content can be defined as residual moisture content and thus, effective porosity = porosity – residual moisture content (Unit: dimensionless, V-void/V-bulk soil).
  • Maximum Deficit: the total or maximum amount of water that is required to saturated the soil when the soil is totally “dry” or at its minimum moisture content. The maximum deficit can be estimated as the product of active soil layer depth and effective porosity (Unit: inch).

The concept of the Deficit and Constant Loss Infiltration Method can be explained in Figure 1. At the beginning of a rainfall event, the soil usually is not saturated and the precipitation (after surface/depression storage and canopy/evapotranspiration ) will penetrate the soil to fill the voids (Phase 1 and 2 in Figure 1). When the soil becomes saturated (soil water deficit =0), the infiltrated water will percolate through the bottom of the active soil layer and get lost. The Deficit and Constant Loss Method can be used for continuous simulation if a canopy/evapotranspiration method is used at the same time to extract water from the soil during “dry” periods.

Figure 1

The three required parameters of the deficit and constant loss infiltration method (refer to this post for the steps to create a soils layer and an infiltration layer in RAS Mapper), namely initial deficit (inch), maximum deficit (inch), and percolation rate (in/hr), need to be calibrated whenever it is possible.

The active soil layer depth (D) is required to calculate initial deficit and maximum deficit as explained in Figure 1. Usually the root zone depth can be used to approximate the active soil layer depth. Robert Pitt in the paper titled of “Evapotranspiration and Related Calculations for Bioretention Devices” shows root depth varies from 1.0 ft to 6.0 ft (Table 1). Ying Fan summarized 2200 observations of root depths and concluded that root depth varies from 0.01 meter to 70 meter with a distribution peak at 1.0 meter (Figure 2). Mean root depths of different soil textures and vegetations are also provided by Ying Fan in the same study (Figure 3 and Figure 4).

Table 1
Figure 2
Figure 3
Figure 4

Under the field condition, it is common that the lower portions of the root zone are not active during the hydrology process and therefore the active soil layer depth can range from half of the root zone depth to the full root zone depth (less than 1 foot to several feet). For lack of a better estimation or guidance, in reference to the study results mentioned above, the active soil layer depth may be estimated as a value between 2.0 ft (considering the fact that lower portion of root zone is not active) and 3.0 ft (close to 1.0 meter – the root depth distribution peak) .

The initial deficit is the amount of water that is required to fill the voids in the soil at the commencement of a simulation and it is usually a value between zero (saturated) and maximum deficit (with residual water content). For a model to simulate a historical storm event, a practical way to estimate the initial deficit is to begin the simulation three days after the end of a previous precipitation event which saturated the soil layer. Three days is the approximate time required for saturated soil to drain freely by gravity to a water content known as field capacity – in this case, the initial deficit can be estimated as the product of active soil layer depth and the difference between porosity (not effective porosity!!! – see Figure 5 for illustration) and field capacity, D * (porosity – field capacity). When simulating a historical storm event it is always a good practice to star the model run in advance enough from the time period of interest so that the modeling results are insensitive to the initial condition values.

Under the average antecedent moisture condition (average AMC), Dr. David Maidment in his book of Handbook of Hydrology suggested using wilting point (-1500 kPa water content) as initial water content in the western USA and using field capacity (-33kPa water content) as initial water content in the eastern USA. Correspondingly, the initial deficit will be calculated as D * (porosity – wilting point) for Western USA or D * (porosity – field capacity) for Eastern USA for average AMC. The typical values of (porosity – wilting point) or (porosity – field capacity) can be found in Table 2 of this post about Green-Ampt infiltration method.

Figure 5

The maximum deficit can be estimated as the product of active soil layer depth and effective porosity, or D * effective porosity. Some literatures suggest using D * (porosity – wilting point), however, since the realistic minimum soil moisture content is the residual moisture content, not wilting point – a water content limit meant for plants to absorb water, it is determined that “D * effective porosity” is more reasonable to calculate the maximum deficit.

The percolation rate in the Deficit and Constant Loss method can be initially estimated as hydraulic conductivity. In field, the percolation rate tends to be less than the saturated hydraulic conductivity due to incomplete saturation at the bottom of the soil layer or root zone. For this reason, the percolation rate may be reduced from its initial estimation during calibration. At each time step during a model run, the entered percolation rate in HEC-HMS or HEC-RAS will be compared with the instant precipitation rate and the lesser of the two will be used to calculate the loss.

Some useful soil properties discussed above are summarized in Table 2 for the typical soil texture. The values in the table are taken from USACE EM 1110-2-1417, 1982 paper authored by Rawls et al, and 1983 paper authored by Rawls et al. It is not clear why the two papers authored by Rawls in 1982 and 1983 presented two sets of hydraulic conductivity values. For deficit and constant infiltration method, one way of dealing with the two sets of values is to apply the 1982 paper hydraulic conductivity values for initial model setup and use the 1983 paper values as calibration limits. Note this is different from what is suggested in this post mainly because the 1983 paper is specifically titled for Green-Ampt Infiltration Method.

Table 2

InfoWorks ICM has two options for rainfall-runoff process – the traditional subcatchment hydrology method and the Rain-on-Grid method. Currently as of 2021, the deficit and constant infiltration method can only be applied to the traditional subcatchment hydrology method, and it is not available for Rain-on-Grid.

In XPSWMM, the deficit and constant loss method is named as uniform loss (Figure 6). The maximum deficit is not a required parameter by XPSWMM, since the uniform loss method does not support continuous simulation and thus there is no soil water content regeneration in XPSWMM when using uniform loss as the infiltration method.

Figure 6

For illustration purpose, a HEC-HMS model result is shown in Figure 7 where the loss method is Deficit and Constant.

Figure 7

Leave a Reply

Your email address will not be published. Required fields are marked *