SCS defines three antecedent moisture conditions: I—dry (wilting point), II—average moisture, and III—wet (field capacity). The moisture condition I curve number is the lowest value the daily curve number can assume in dry conditions. The curve numbers for moisture conditions I and III are calculated with the equations:
CN1=CN2−(100−CN2+exp[2.533−0.0636∗(100−CN2)])20∗(100−CN2) 2:1.1.4
CN3=CN2∗exp[0.00673∗(100−CN2)] 2:1.1.5
where CN1 is the moisture condition I curve number, CN2 is the moisture condition II curve number, and CN3 is the moisture condition III curve number.
With SWAT+, users are allowed to select between two methods for calculating the retention parameter. The traditional method is to allow the retention parameter to vary with soil profile water content. An alternative added in SWAT+ allows the retention parameter to vary with accumulated plant evapotranspiration. Calculation of the daily CN value as a function of plant evapotranspiration was added because the soil moisture method was predicting too much runoff in shallow soils. By calculating daily CN as a function of plant evapotranspiration, the value is less dependent on soil storage and more dependent on antecedent climate.
When the retention parameter varies with soil profile water content, the following equation is used:
2:1.1.6
where is the retention parameter for a given day (mm), is the maximum value the retention parameter can achieve on any given day (mm), is the soil water content of the entire profile excluding the amount of water held in the profile at wilting point (mm HO), and
The SCS curve number is a function of the soil’s permeability, land use and antecedent soil water conditions. Typical curve numbers for moisture condition II are listed in tables 2:1-1, 2:1-2 and 2:1-3 for various land covers and soil types (SCS Engineering Division, 1986). These values are appropriate for a 5% slope.
Table 2:1-1: Runoff curve numbers for cultivated agricultural lands (from SCS Engineering Division, 1986)
Table 2:1-2: Runoff curve numbers for other agricultural lands (from SCS Engineering Division, 1986)
Table 2:1-3: Runoff curve numbers for urban areas (from SCS Engineering Division, 1986)
The U.S. Natural Resource Conservation Service (NRCS) classifies soils into four hydrologic groups based on infiltration characteristics of the soils. NRCS Soil Survey Staff (1996) defines a hydrologic group as a group of soils having similar runoff potential under similar storm and cover conditions. Soil properties that influence runoff potential are those that impact the minimum rate of infiltration for a bare soil after prolonged wetting and when not frozen.
These properties are depth to seasonally high water table, saturated hydraulic conductivity, and depth to a very slowly permeable layer.
Soil may be placed in one of four groups, A, B, C, and D, or three dual classes, A/D, B/D, and C/D.
Definitions of the classes are:
A: (Low runoff potential). The soils have a high infiltration rate even when thoroughly wetted. They chiefly consist of deep, well drained to excessively drained sands or gravels. They have a high rate of water transmission.
B: The soils have a moderate infiltration rate when thoroughly wetted. They chiefly are moderately deep to deep, moderately well-drained to well-drained soils that have moderately fine to moderately coarse textures. They have a moderate rate of water transmission.
C: The soils have a slow infiltration rate when thoroughly wetted. They chiefly have a layer that impedes downward movement of water or have moderately fine to fine texture. They have a slow rate of water transmission.
D: (High runoff potential). The soils have a very slow infiltration rate when thoroughly wetted. They chiefly consist of clay soils that have a high swelling potential, soils that have a permanent water table, soils that have a claypan or clay layer at or near the surface, and shallow soils over nearly impervious material. They have a very slow rate of water transmission.
Dual hydrologic groups are given for certain wet soils that can be adequately drained. The first letter applies to the drained condition, the second to the undrained. Only soils that are rated D in their natural condition are assigned to dual classes.
The shape coefficients are determined by solving equation 2:1.1.6 assuming that
the retention parameter for moisture condition I curve number corresponds to wilting point soil profile water content,
the retention parameter for moisture condition III curve number corresponds to field capacity soil profile water content, and
the soil has a curve number of 99 (S = 2.54) when completely saturated.
w1=ln[1−S3∗Smax−1FC−FC]+w2∗FC 2.1.1.7
w2=(SAT−FC)(ln[1−S3∗Smax−1FC−FC]−ln[1−2.54∗Smax−1SAT−SAT]) 2.1.1.8
where w1 is the first shape coefficient, w2 is the second shape coefficient, FC is the amount of water in the soil profile at field capacity (mm H2O), S3 is the retention parameter for the moisture condition III curve number, Smax is the retention parameter for the moisture condition I curve number, SAT is the amount of water in the soil profile when completely saturated (mm H2O), and 2.54 is the retention parameter value for a curve number of 99.
When the retention parameter varies with plant evapotranspiration, the following equation is used to update the retention parameter at the end of every day:
S=Sprev+Eo∗exp(Smax−cncoef∗Sprev)−Rday+Qsurf 2:1.1.9
where S is the retention parameter for a given day (mm), Sprev is the retention parameter for the previous day (mm), Eo is the potential evapotranspiration for the day (mm d−1), cncoef is the weighting coefficient used to calculate the retention coefficient for daily curve number calculations dependent on plant evapotranspiration, Smax is the maximum value the retention parameter can achieve on any given day (mm), Rday is the rainfall depth for the day (mm H2O), and Qsurf is the surface runoff (mm H2O). The initial value of the retention parameter is defined as
S=0.9∗Smax
When the top layer of the soil is frozen, the retention parameter is modified using the following equation:
Sfrz=Smax∗[1−exp(−0.000862∗S)] 2:1.1.10
where Sfrz is the retention parameter adjusted for frozen conditions (mm), Smax is the maximum value the retention parameter can achieve on any given day (mm), and S is the retention parameter for a given moisture content calculated with equation 2:1.1.6 (mm).
The daily curve number value adjusted for moisture content is calculated by rearranging equation 2:1.1.2 and inserting the retention parameter calculated for that moisture content:
CN=(S+254)25400 2:1.1.11
where CN is the curve number on a given day and S is the retention parameter calculated for the moisture content of the soil on that day.
The moisture condition II curve numbers provided in the tables are assumed to be appropriate for 5% slopes. Williams (1995) developed an equation to adjust the curve number to a different slope:
CN2s=3(CN3−CN2)∗[1−2∗exp(−13.86∗slp)]+CN2 2:1.1.12
where CN2s is the moisture condition II curve number adjusted for slope, CN3 is the moisture condition III curve number for the default 5% slope, CN2 is the moisture condition II curve number for the default 5% slope, and slp is the average fraction slope of the subbasin. SWAT+ does not adjust curve numbers for slope. If the user wishes to adjust the curve numbers for slope effects, the adjustment must be done prior to entering the curve numbers in the management input file.
Table 2:1-1: SWAT+ input variables that pertain to surface runoff calculated with the SCS curve number method.
: Moisture condition II curve number
CN2
cn2
: Moisture condition II curve number
CNOP
.mgt
Rainfall, runoff, routing option.
IEVENT
.bsn
Daily curve number calculation method: 0 calculate daily CN value as a function of soil moisture; 1 calculate daily CN value as a function of plant evapotranspiration
ICN
.bsn
cncoef: Weighting coefficient used to calculate the retention coefficient for daily curve number calculations dependent on plant evapotranspiration
latq_co
epco
Rday: Daily precipitation (mm H2O)
pcp
pcp