6.3 NRCS (SCS) Curve Number Method & Runoff Volume Estimation

6.3 NRCS (SCS) Curve Number Method & Runoff Volume Estimation

Quick Reference: The NRCS (SCS) Curve Number (CN) Method is the standard hydrologic modeling approach for computing both peak discharge ($Q_{peak}$) and total storm runoff volume ($V$). Governed by the empirical parameter $CN$ ($0\text{ to }100$), potential soil retention is expressed as $S = (1000 / CN) - 10$, with initial abstraction standardly assumed as $I_a = 0.2S$. Direct runoff depth is calculated as $Q = (P - 0.2S)^2 / (P + 0.8S)$ when rainfall $P > 0.2S$. On construction sites, mass grading converts low-CN soils into bare ground with $CN$ values ranging from 77 (HSG A) to 94 (HSG D), dramatically inflating runoff volume and dictating sediment basin sizing.

---

NRCS TR-55 Methodology & Volumetric Modeling Capabilities

Developed in 1954 by the USDA Soil Conservation Service (now the Natural Resources Conservation Service [NRCS]) and published in the National Engineering Handbook Section 4 (NEH-4) and Technical Release 55 (TR-55), the Curve Number method overcomes the fundamental limitations of the Rational Method:

1. Dual Output (Peak & Volume): While the Rational Method yields only a single instantaneous peak flow rate ($Q_{peak}$), the NRCS method computes both peak discharge ($Q_{peak}$) and total cumulative runoff volume ($V$).
2. Hydrograph Generation: By coupling runoff volume with dimensionless unit hydrograph theory, the NRCS method synthesizes complete discharge hydrographs (flow rate plotted over time), enabling dynamic reservoir flood routing through emergency spillways, orifices, and surface skimmers.
3. Watershed Scale Flexibility: While the Rational Method is constrained to small catchments ($\le 20\text{–}100\text{ acres}$), the NRCS CN method is robust across both small parcels and large, complex regional watersheds spanning thousands of acres.
4. Essential for Sediment Basin Sizing: Sediment settling is governed by Stokes' Law and hydraulic detention residence time. Sizing a sediment basin or sediment trap requires knowing the total volume of water entering the impoundment during a design event (typically the 2-year or 10-year, 24-hour storm) to ensure particles have adequate time to settle before effluent exits the outlet.

---

Curve Number ($CN$) Physics & Hydrologic Soil Groups

The Runoff Curve Number ($CN$) is an empirical, dimensionless index ranging from $0\text{ to }100$ that quantifies the runoff potential of a specific land surface. A theoretical $CN$ of $0$ represents a completely permeable surface with infinite infiltration (zero runoff), whereas a $CN$ of $100$ represents a completely impervious surface or open water body where $100%$ of precipitation is converted into surface runoff.

The Curve Number is dictated by three primary landscape variables:

• Hydrologic Soil Group (HSG: A, B, C, or D)
• Land Use and Surface Cover Type
• Hydrologic Condition and Antecedent Moisture Condition (AMC)

NRCS Hydrologic Soil Groups (HSGs)

The NRCS classifies over 20,000 national soil series into four distinct Hydrologic Soil Groups based on their minimum steady-state infiltration rates after prolonged wetting:

• Group A (Low Runoff Potential): Deep, well-drained to excessively drained sands and gravelly sands. Saturated hydraulic conductivity is high ($K_{sat} > 1.42\text{ in/hr}$ or $> 10.0\ \mu\text{m/s}$), and water transmission through the soil column is rapid.
• Group B (Moderately Low Runoff Potential): Moderately deep to deep, moderately well-drained to well-drained soils with moderately fine to moderately coarse textures (sandy loams, silt loams). Water transmission rates range from $0.57\text{ to }1.42\text{ in/hr}$ ($4.0\text{ to }10.0\ \mu\text{m/s}$).
• Group C (Moderately High Runoff Potential): Soils with restrictive soil horizons that impede downward movement, or soils with moderately fine to fine textures (clay loams, silty clay loams). Water transmission rates are low, ranging from $0.06\text{ to }0.57\text{ in/hr}$ ($0.4\text{ to }4.0\ \mu\text{m/s}$).
• Group D (High Runoff Potential): Heavy clay soils with high swelling potential, soils with a permanent high water table within 2 feet of the surface, soils with a claypan or hardpan near the surface, and shallow soils over nearly impermeable bedrock. Saturated hydraulic conductivity is extremely low ($K_{sat} < 0.06\text{ in/hr}$ or $< 0.4\ \mu\text{m/s}$).

Antecedent Moisture Conditions (AMC)

Curve Numbers in standard engineering tables represent Antecedent Moisture Condition II (AMC II), which reflects average, normal soil moisture conditions prior to the design storm. The NRCS recognizes two boundary conditions:

• AMC I (Dry): Soils are dry but not to the permanent wilting point; 5-day total antecedent rainfall is $< 0.5\text{ inches}$ in the dormant season or $< 1.4\text{ inches}$ in the growing season. Infiltration is higher, and the effective Curve Number is lower.
• AMC III (Wet / Saturated): Soils are virtually saturated from recent heavy precipitation; 5-day antecedent rainfall is $> 1.1\text{ inches}$ in the dormant season or $> 2.1\text{ inches}$ in the growing season. Infiltration capacity is heavily depressed, yielding significantly elevated Curve Numbers and near-immediate runoff generation.

---

Potential Maximum Soil Retention ($S$) & Initial Abstraction ($I_a$)

The fundamental physical premise of the NRCS Curve Number method is that the total precipitation ($P$) falling on a watershed is partitioned into three discrete components:

$P = I_a + F + Q$

Where:

• $P$ = Total storm precipitation depth (inches)
• $I_a$ = Initial abstraction (inches)
• $F$ = Actual retention (cumulative infiltration into the soil after runoff begins, inches)
• $Q$ = Direct surface runoff depth (inches)

┌─────────────────────────────────────────────────────────────┐
│                     Total Precipitation (P)                 │
└─────────────────────────────────────────────────────────────┘
        │                          │                     │
        ▼                          ▼                     ▼
┌─────────────────┐      ┌─────────────────┐   ┌─────────────────┐
│ Initial Abstr.  │      │ Actual Infiltr. │   │ Direct Runoff   │
│     (Ia)        │      │       (F)       │   │       (Q)       │
│ Surface storage,│      │ Water absorbed  │   │ Overland flow   │
│ canopy wetting, │      │ into soil after │   │ reaching        │
│ initial soaking │      │ runoff begins   │   │ streams/basins  │
└─────────────────┘      └─────────────────┘   └─────────────────┘

Mathematical Formulation of Soil Retention ($S$)

The potential maximum soil retention ($S$) represents the total upper-bound storage capacity of the soil profile (combined maximum infiltration and depression storage) after runoff begins. It is related directly to the dimensionless Curve Number by:

$S = \frac{1000}{CN} - 10 \quad [\text{in inches}]$

Metric Formulation: In SI metric units, potential maximum retention is expressed as: $S = \frac{25400}{CN} - 254 \quad [\text{in millimeters}]$

Initial Abstraction ($I_a$)

Initial abstraction ($I_a$) comprises all water losses that occur before surface runoff begins. It includes surface ponding in micro-depressions, interception by vegetative canopy and ground litter, evapotranspiration during the storm, and the initial rapid capillary infiltration into unsaturated topsoil.

Based on historical empirical calibrations across hundreds of experimental agricultural watersheds, the NRCS established the standard default relationship:

$I_a = 0.2 \times S$

Substituting this relationship into the retention formula reveals that as $CN$ approaches $100$, $S$ approaches $0$, and $I_a$ approaches $0$ (zero losses, $100%$ runoff). Conversely, lower $CN$ values produce large $S$ and $I_a$ values, meaning significant rainfall depths must occur before any runoff is initiated.

---

The NRCS Runoff Equation & Mathematical Derivation

The NRCS derived its governing rainfall-runoff equation from a simple continuity hypothesis: the ratio of actual retention ($F$) to potential maximum retention ($S$) equals the ratio of direct runoff ($Q$) to potential maximum runoff ($P - I_a$):

$\frac{F}{S} = \frac{Q}{P - I_a}$

From water balance continuity, actual retention after runoff initiates is the precipitation remaining after subtracting initial abstraction and direct runoff: $F = (P - I_a) - Q$. Substituting this into the continuity ratio yields:

$\frac{(P - I_a) - Q}{S} = \frac{Q}{P - I_a}$

Multiplying both sides by $S \times (P - I_a)$:

$(P - I_a)^2 - Q(P - I_a) = Q \times S$

$(P - I_a)^2 = Q \times S + Q(P - I_a) = Q \times [S + P - I_a]$

Solving explicitly for direct runoff depth ($Q$):

$Q = \frac{(P - I_a)^2}{(P - I_a) + S}$

Substituting the standard empirical assumption $I_a = 0.2S$:

$Q = \frac{(P - 0.2S)^2}{(P - 0.2S) + S} = \frac{(P - 0.2S)^2}{P + 0.8S}$

Governing Runoff Threshold Condition

This equation is valid only when precipitation exceeds initial abstraction:

$Q = \begin{cases} \frac{(P - 0.2S)^2}{P + 0.8S} & \text{for } P > 0.2S \\ 0 & \text{for } P \le 0.2S \end{cases}$

If total rainfall $P$ is less than or equal to $0.2S$, then $Q = 0$. Every drop of rain is completely absorbed by depression storage, interception, and initial infiltration, producing zero surface runoff.

---

Calculating Total Storm Runoff Volume ($V$)

Once direct runoff depth ($Q$, in inches) is determined from the NRCS equation, calculating the total cumulative storm runoff volume ($V$) over a drainage catchment of area $A$ (in acres) is straightforward.

Volume in Acre-Feet

$\text{Runoff Volume } (V_{\text{ac-ft}}) = \frac{Q \text{ (inches)} \times A \text{ (acres)}}{12\text{ inches/foot}}$

Volume in Cubic Feet

Because $1\text{ acre} = 43,560\text{ sq ft}$:

$\text{Runoff Volume } (V_{\text{cu ft}}) = V_{\text{ac-ft}} \times 43,560\text{ ft}^2/\text{acre} = \frac{Q \times A \times 43,560}{12} = Q \times A \times 3,630$

Regulatory Significance for EPA Construction General Permit (CGP)

Where a state permit requires a temporary sediment basin (commonly at 5 or 10 or more disturbed acres draining to a common point), both the state permit and EPA CGP Part 2.2.12 accept a basin sized to the larger of two storage volumes:

1. $3,600\text{ cubic feet}$ of storage per acre drained; OR
2. The calculated volume of runoff from a 2-year, 24-hour storm event from all contributing upland drainage area.

The NRCS Curve Number method is the universal tool used to calculate this 2-year, 24-hour runoff volume ($V$) to verify whether a proposed sediment basin meets this critical regulatory volume standard.

---

Curve Numbers for Construction Sites vs. Pre-Development

Table 2-2 of NRCS TR-55 provides standard Curve Numbers across cover types and hydrologic soil groups for average moisture conditions (AMC II). On construction sites, the removal of vegetation and disturbance of topsoil elevates Curve Numbers significantly:

The Dramatic Hydrologic Shift from Woods to Bare Soil

Consider an undisturbed woodland site in Hydrologic Soil Group C ($CN = 70$) undergoing clearing and mass grading into bare soil ($CN = 91$):

• Pre-Development ($CN = 70$): $S = \frac{1000}{70} - 10 = 4.29\text{ inches}; \quad I_a = 0.2 \times 4.29 = 0.86\text{ inches}$
• Post-Grading ($CN = 91$): $S = \frac{1000}{91} - 10 = 0.99\text{ inches}; \quad I_a = 0.2 \times 0.99 = 0.20\text{ inches}$

Potential soil water retention ($S$) plummets by $77%$, and the initial abstraction threshold required to initiate surface runoff drops from $0.86\text{ inches}$ to a mere $0.20\text{ inches}$. Light rain events that previously generated zero runoff now produce immediate, sediment-laden discharges.

---

Unit Hydrograph Theory & Synthetic 24-Hour Rainfall Distributions

To transform total runoff volume ($Q$, in inches) into a continuous time-discharge hydrograph and determine peak flow ($Q_{peak}$), the NRCS utilizes Dimensionless Unit Hydrograph theory (developed by Victor Mockus in 1957):

$q_p = \frac{484 \times A \times Q}{t_p}$

Where:

• $q_p$ = Peak discharge rate (cfs)
• $484$ = Standard NRCS peak rate factor for curvilinear unit hydrographs (reflecting $37.5%$ of hydrograph volume under the rising limb)
• $A$ = Drainage area (square miles, where $1\text{ sq mi} = 640\text{ acres}$)
• $Q$ = Direct runoff depth (inches)
• $t_p$ = Time to peak (hours), where $t_p = \frac{\Delta t}{2} + 0.6 T_c$

SCS 24-Hour Synthetic Rainfall Distributions

Because rainfall intensity varies dramatically during a 24-hour storm, the NRCS developed four regional synthetic 24-hour rainfall distributions based on National Weather Service records:

• Type I: Maritime climate with wet winters and dry summers (Southern and Central California, Western Oregon and Washington, Hawaii).
• Type IA: Pacific Northwest coastal regions with low-intensity, long-duration persistent winter rains.
• Type II: The most intense convective storm distribution; applicable to the vast interior continental United States, the Midwest, Great Plains, and Northeast. Approximately $50%$ of the entire 24-hour rainfall depth falls within a single 2-hour window centered at hour 12 ($t = 11.0\text{ to }13.0\text{ hours}$), creating extreme peak discharge spikes.
• Type III: Gulf Coast and Southeastern coastal areas characterized by tropical moisture surges, high-volume hurricanes, and intense subtropical squalls.

---

Full Worked NRCS Runoff Volume Calculation

Problem Statement

A 20.0-acre commercial site in Columbus, Ohio is undergoing complete clearing and mass grading (bare soil). Soil survey mapping confirms the underlying soil is Miami silt loam, classified as Hydrologic Soil Group C.

The CPESC practitioner must determine the total runoff volume generated by the 10-year, 24-hour design storm to size an on-site temporary sediment basin. From NOAA Atlas 14, the 10-year, 24-hour design precipitation depth is $P = 4.50\text{ inches}$.

Step 1: Determine the Design Curve Number ($CN$)

From TR-55 Table 2-2, for newly graded areas (bare soil) on Hydrologic Soil Group C:

$CN = 91$

Step 2: Calculate Potential Maximum Soil Retention ($S$)

$S = \frac{1000}{CN} - 10$

$S = \frac{1000}{91} - 10 = 10.989 - 10 = 0.989\text{ inches}$

Step 3: Calculate Initial Abstraction ($I_a$)

$I_a = 0.2 \times S = 0.2 \times 0.989\text{ in} = 0.198\text{ inches}$

Check runoff initiation criterion: $P = 4.50\text{ in} > I_a = 0.198\text{ in} \implies \text{Runoff occurs}$

Step 4: Calculate Direct Runoff Depth ($Q$)

$P - 0.2S = 4.50 - 0.198 = 4.302\text{ inches}$

$(P - 0.2S)^2 = (4.302)^2 = 18.507\text{ in}^2$

$P + 0.8S = 4.50 + (0.8 \times 0.989) = 4.50 + 0.791 = 5.291\text{ inches}$

$Q = \frac{(P - 0.2S)^2}{P + 0.8S} = \frac{18.507}{5.291} = 3.498\text{ inches} \approx 3.50\text{ inches}$

Out of 4.50 inches of rainfall, 3.50 inches (77.8%) becomes direct surface runoff!

Step 5: Calculate Total Storm Runoff Volume ($V$)

Volume in Acre-Feet:

$V_{\text{ac-ft}} = \frac{Q \times A}{12} = \frac{3.498\text{ in} \times 20.0\text{ acres}}{12\text{ in/ft}} = \frac{69.96}{12} = 5.830\text{ acre-feet}$

Volume in Cubic Feet:

$V_{\text{cu ft}} = V_{\text{ac-ft}} \times 43,560\text{ ft}^2/\text{acre} = 5.830\text{ ac-ft} \times 43,560\text{ ft}^2/\text{ac} = 253,955\text{ ft}^3 \approx 254,000\text{ ft}^3$

(Alternatively: $V_{\text{cu ft}} = Q \times A \times 3,630 = 3.498 \times 20.0 \times 3,630 = 253,955\text{ cu ft}$)

Step 6: Engineering Evaluation & Sediment Basin Sizing Comparison

• EPA CGP Standard Minimum Storage Rule: $V_{\text{min}} = 20.0\text{ acres} \times 3,600\text{ ft}^3/\text{acre} = 72,000\text{ cubic feet}$
• Actual 10-Year, 24-Hour Storm Runoff Volume: $V_{10\text{-yr}} = 253,955\text{ cubic feet}$

Comparing these values demonstrates that the actual 10-year storm generates over 3.5 times the default 3,600 cu ft/acre volume standard ($254,000\text{ ft}^3$ vs. $72,000\text{ ft}^3$). If a CPESC professional designs the sediment basin strictly using the static $3,600\text{ ft}^3/\text{acre}$ rule without modeling the 10-year storm volume, the basin will be severely overwhelmed during a 10-year event, overtopping the embankment and causing catastrophic downstream sediment pollution.

---

Isohyetal Maps and the Composite Curve Number

Two named blueprint items sit either side of the curve number calculation: where the rainfall depth comes from, and what to do when the watershed is not uniform.

Isohyetal Maps — Do Not Confuse Them with Isoerodent Maps

An isohyetal map contours lines of equal rainfall depth. There are two distinct uses:

• Design rainfall. NOAA Atlas 14 (and Atlas 2 in the states it still covers) publishes contoured depth-duration-frequency maps: the 2-year 24-hour depth, the 10-year 24-hour depth, and so on. Reading a design depth $P$ off an isohyetal map or the NOAA Precipitation Frequency Data Server (PFDS) is the first step of every curve-number calculation, and the depth may vary meaningfully across a single large or linear project.
• Storm analysis. Contouring gauge totals after an actual event yields the mean areal precipitation over a watershed by the isohyetal method, which weights each inter-contour band by its area — more accurate than a simple arithmetic gauge average when rainfall is spatially uneven.

Mixing these up is a classic exam distractor: a value read off an isohyetal map is inches, a value read off an isoerodent map is a dimensionless erosivity index.

The Composite (Weighted) Curve Number

Real watersheds are mosaics of soil groups and cover types. The composite curve number — also called the weighted curve number — is the area-weighted average:

$CN_{\text{comp}} = \frac{\sum (CN_i \times A_i)}{\sum A_i}$

Worked example. A 20.0-acre catchment during construction: 12.0 ac of bare graded HSG C soil ($CN = 91$), 5.0 ac of undisturbed woods in good condition on HSG C ($CN = 70$), and 3.0 ac of new asphalt ($CN = 98$). $CN_{\text{comp}} = \frac{(91)(12.0) + (70)(5.0) + (98)(3.0)}{20.0} = \frac{1{,}092 + 350 + 294}{20.0} = \frac{1{,}736}{20.0} = 86.8 \approx 87$ Then $S = \frac{1000}{87} - 10 = 1.49$ in, and for $P = 3.5$ in: $Q = \frac{(3.5 - 0.298)^2}{3.5 + 1.195} = \frac{10.25}{4.695} = 2.18$ in of direct runoff.

Three cautions the CPESC must apply:

1. Never average curve numbers across sub-basins that drain to different points. Weight only within a single hydrologic response area.
2. Do not composite impervious areas with pervious areas when total imperviousness exceeds about 30%, or when the impervious area is unconnected. TR-55 provides separate adjustment figures for those cases; a naive weighted average understates the peak.
3. Recompute for each construction phase. The CN of a site changes from pre-development woods, to stripped subgrade, to final landscaped and paved condition, and the basin must be sized on the worst of them — normally the fully stripped condition.