HEC-HMS Subbasin Elements: Defining Your Watershed
Subbasins are the fundamental building blocks of HEC-HMS hydrologic models. They represent land areas that receive precipitation and generate runoff. Understanding how to properly configure subbasin elements is essential for accurate watershed simulation.
What Is a Subbasin Element?
A subbasin in HEC-HMS represents a land surface area that:
- Receives precipitation from the meteorologic model
- Computes losses (infiltration, interception, depression storage)
- Transforms excess precipitation into a runoff hydrograph
- May include baseflow from groundwater contributions
- Discharges to a downstream element (junction, reach, reservoir, or sink)
Subbasin Properties Overview
Every subbasin element requires configuration of several key properties:
| Property | Description | Required |
|---|---|---|
| Area | Drainage area of the subbasin | Yes |
| Downstream | Element receiving discharge | Yes |
| Loss Method | How infiltration is calculated | Yes |
| Transform Method | How runoff hydrograph is generated | Yes |
| Baseflow Method | Groundwater contribution | Optional |
| Canopy | Interception by vegetation | Optional |
| Surface | Surface storage and routing | Optional |
Area and Units
The area property defines the size of the drainage area contributing to the subbasin.
Area Units by Unit System
| Unit System | Available Units |
|---|---|
| US Customary | Square miles (MI2), Acres (AC) |
| SI (Metric) | Square kilometers (KM2), Hectares (HA) |
Determining Subbasin Area
Methods for calculating drainage area:
- GIS Analysis: Delineate watershed from DEM using HEC-GeoHMS or ArcGIS
- Topographic Maps: Trace drainage divides on contour maps
- Survey Data: Use site-specific survey information
- Planimeter: Measure area on scaled maps
Loss Methods Overview
Loss methods determine how much precipitation infiltrates into the ground, is intercepted by vegetation, or is stored in surface depressions. The remaining precipitation becomes excess precipitation that generates runoff.
HEC-HMS offers several loss methods:
| Method | Complexity | Data Requirements | Best For |
|---|---|---|---|
| Initial and Constant | Low | 2 parameters | Quick estimates |
| SCS Curve Number | Low | 1-3 parameters | Most common applications |
| Green-Ampt | Medium | 4 parameters | Physical soil modeling |
| Deficit and Constant | Medium | 3+ parameters | Continuous simulation |
| Soil Moisture Accounting | High | 10+ parameters | Detailed continuous modeling |
Initial and Constant Loss
The simplest loss method with two parameters:
Where:
- Ia = Initial abstraction (inches or mm)
- f = Constant loss rate (inches/hour or mm/hour)
- t = Time
| Parameter | Description | Typical Values |
|---|---|---|
| Initial Loss | Depth before runoff begins | 0.1 - 0.5 inches |
| Constant Rate | Ongoing infiltration rate | 0.1 - 1.0 in/hr |
| Impervious | Percent directly connected | 0 - 100% |
Use when:
- Limited soil data available
- Rough preliminary analysis
- Calibrating to observed data
SCS Curve Number Loss
The most widely used loss method, based on NRCS (formerly SCS) methodology.
Where:
- Pe = Excess precipitation (runoff)
- P = Total precipitation
- Ia = Initial abstraction (typically 0.2S)
- S = Potential maximum retention
| Parameter | Description | Range |
|---|---|---|
| Curve Number | Runoff potential index | 30 - 98 |
| Initial Abstraction | Depth before runoff | 0 - 5 inches |
| Impervious | Directly connected impervious area | 0 - 100% |
Curve Number Reference Table
| Land Use | Soil Group A | Soil Group B | Soil Group C | Soil Group D |
|---|---|---|---|---|
| Open space (good condition) | 39 | 61 | 74 | 80 |
| Residential (1/2 acre) | 54 | 70 | 80 | 85 |
| Residential (1/4 acre) | 61 | 75 | 83 | 87 |
| Commercial/Industrial | 89 | 92 | 94 | 95 |
| Parking lots, roads | 98 | 98 | 98 | 98 |
| Wooded (good condition) | 30 | 55 | 70 | 77 |
| Pasture (good condition) | 39 | 61 | 74 | 80 |
| Row crops (contoured) | 65 | 75 | 82 | 86 |
Green-Ampt Loss
A physically-based infiltration method using soil hydraulic properties:
Where:
- f = Infiltration rate
- Ks = Saturated hydraulic conductivity
- psi = Wetting front suction head
- theta = Soil moisture deficit
- F = Cumulative infiltration
| Parameter | Description | Typical Values |
|---|---|---|
| Saturated Content | Porosity | 0.35 - 0.50 |
| Initial Content | Starting moisture | 0.10 - 0.30 |
| Suction | Wetting front head | 2 - 12 inches |
| Conductivity | Saturated K | 0.01 - 5.0 in/hr |
Use when:
- Soil physical properties are known
- Modeling infiltration physics is important
- Site-specific soil data is available
Deficit and Constant Loss
An extension of Initial and Constant for continuous simulation:
| Parameter | Description |
|---|---|
| Initial Deficit | Starting soil moisture deficit |
| Maximum Deficit | Maximum potential storage |
| Constant Rate | Long-term infiltration rate |
| Impervious | Directly connected impervious area |
| Recovery Rate | Rate of soil moisture recovery |
Use when:
- Running continuous simulations over weeks or months
- Soil moisture variation between storms matters
- Recovery of infiltration capacity is important
Transform Methods Overview
Transform methods convert excess precipitation (after losses) into a runoff hydrograph at the subbasin outlet. The relationship between precipitation and discharge involves timing, attenuation, and hydrograph shape.
| Method | Complexity | Key Parameters | Best For |
|---|---|---|---|
| SCS Unit Hydrograph | Low | Lag time only | Most applications |
| Clark Unit Hydrograph | Medium | Tc and R | Larger watersheds |
| Snyder Unit Hydrograph | Medium | Tp and Cp | Traditional studies |
| ModClark | High | Gridded Tc, R | Distributed modeling |
| User-Specified Unit Hydrograph | Variable | UH ordinates | Calibrated models |
SCS Unit Hydrograph
The most commonly used transform method, based on the SCS dimensionless unit hydrograph.
The SCS dimensionless unit hydrograph has a peak at:
Where:
- Tp = Time to peak
- D = Computational time step
- tlag = Subbasin lag time
| Parameter | Description |
|---|---|
| Lag Time | Time from centroid of excess precipitation to peak discharge |
| Graph Type | Standard (484 peak rate factor) or custom |
Estimating Lag Time
The SCS lag time equation:
Where:
- L = Hydraulic length (feet)
- S = (1000/CN) - 10
- Y = Average watershed slope (percent)
Alternative relationship with time of concentration:
Use our Time of Concentration Calculator to estimate Tc.
Clark Unit Hydrograph
The Clark method uses two parameters to describe watershed response:
- Time of Concentration (Tc): Time for water to travel from the most distant point to the outlet
- Storage Coefficient (R): Linear reservoir coefficient representing storage effects
| Parameter | Description | Typical Range |
|---|---|---|
| Time of Concentration | Maximum travel time | 0.5 - 24 hours |
| Storage Coefficient | Reservoir attenuation | 0.3Tc - 0.8Tc |
| Time-Area Curve | Optional spatial distribution | Default or custom |
Use when:
- Time-area relationship is known or can be derived
- More control over hydrograph shape is needed
- Watersheds with significant storage effects
Snyder Unit Hydrograph
A traditional empirical method using two parameters:
| Parameter | Description | Typical Range |
|---|---|---|
| Lag (Tp) | Basin lag time | Hours |
| Peaking Coefficient (Cp) | Peak rate factor | 0.4 - 0.8 |
Use when:
- Maintaining consistency with previous studies
- Regional Snyder coefficients are available
- Traditional USACE analysis required
ModClark Transform
A distributed version of the Clark method that uses gridded time-area relationships:
- Accounts for spatial distribution of excess precipitation
- Useful with radar rainfall data
- Requires GIS-derived grid cell properties
Use when:
- Using gridded (radar) precipitation
- Spatial variability within subbasin is important
- High-resolution analysis is required
Baseflow Methods Overview
Baseflow represents the contribution of groundwater discharge to stream flow. For event-based modeling, baseflow may be omitted. For continuous simulation or calibration to observed data, baseflow is important.
| Method | Complexity | Description |
|---|---|---|
| None | N/A | No baseflow contribution |
| Constant Monthly | Low | Fixed monthly values |
| Linear Reservoir | Medium | Groundwater storage model |
| Recession | Medium | Exponential recession |
| Bounded Recession | Medium | Recession with limits |
| Nonlinear Boussinesq | High | Physical groundwater model |
Constant Monthly
The simplest baseflow method - specify a constant flow for each month:
| Month | Baseflow (cfs) |
|---|---|
| January | 5.0 |
| February | 6.0 |
| … | … |
Use when:
- Long-term average baseflow is known
- Monthly variation is significant
- Quick estimation is needed
Linear Reservoir
Models groundwater as a linear reservoir that releases water proportional to storage:
Where:
- Qbf = Baseflow discharge
- G = Groundwater coefficient
- Sgw = Groundwater storage
| Parameter | Description |
|---|---|
| Initial Discharge | Starting baseflow rate |
| Recession Constant | Rate of baseflow decline |
| Threshold Ratio | Ratio to peak for recharge |
Use when:
- Groundwater-surface water interaction is important
- Calibrating to observed streamflow
- Continuous simulation with baseflow variation
Recession Method
Models baseflow recession using an exponential decay function:
Where:
- Qt = Baseflow at time t
- Q0 = Initial baseflow
- k = Recession constant (0.8 - 0.98)
| Parameter | Description | Typical Values |
|---|---|---|
| Initial Flow | Starting baseflow | Site-specific |
| Recession Constant | Decay rate | 0.85 - 0.95 |
| Threshold Flow | Recharge trigger | Site-specific |
Typical Parameter Values by Land Use
Urban/Developed Areas
| Land Use | CN (B Soil) | Lag Time | Impervious |
|---|---|---|---|
| Commercial | 92 | 15-30 min | 85% |
| Industrial | 88 | 20-40 min | 72% |
| Residential (1/8 ac) | 85 | 20-45 min | 65% |
| Residential (1/4 ac) | 75 | 30-60 min | 38% |
| Residential (1/2 ac) | 70 | 45-90 min | 25% |
| Residential (1 ac) | 68 | 60-120 min | 20% |
Rural/Agricultural Areas
| Land Use | CN (B Soil) | Lag Time | Impervious |
|---|---|---|---|
| Pasture (good) | 61 | 2-4 hours | 0% |
| Row crops | 78 | 1-3 hours | 0% |
| Small grain | 72 | 1-3 hours | 0% |
| Meadow | 58 | 2-5 hours | 0% |
| Forest (good) | 55 | 3-8 hours | 0% |
Entering Parameters Correctly
Step-by-Step Configuration
- Select the subbasin in the Watershed Explorer
- Set the Area property in the Component Editor
- Choose a Loss Method from the dropdown
- Enter loss parameters in the expanded section
- Choose a Transform Method
- Enter transform parameters
- Choose a Baseflow Method (or None)
- Enter baseflow parameters if applicable
- Set the Downstream connection
Parameter Entry Tips
- All parameters must use consistent units within the project
- Some parameters have typical ranges - values outside these may indicate errors
- Use the ellipsis button (…) to access sub-editors for complex parameters
- Hover over parameter names for tooltips with descriptions
Common Mistakes and How to Avoid Them
Incorrect Curve Number
Problem: CN values not appropriate for land use and soil type Solution: Use NRCS tables or our calculator; verify soil group from soil survey
Wrong Lag Time Units
Problem: Lag time entered in wrong units (hours vs. minutes) Solution: Check unit system; typical urban lag is 15-60 minutes, rural is 1-10 hours
Unrealistic Initial Abstraction
Problem: Ia too large, resulting in no runoff for small storms Solution: Use Ia = 0.2S or Ia = 0.05S (more recent research suggests 0.05)
Missing Downstream Connection
Problem: Subbasin not connected to downstream element Solution: Always verify Downstream property is set to valid element
Impervious Area Double-Counting
Problem: Impervious already in CN but also entered separately Solution: If CN includes impervious, set Impervious parameter to 0
Subbasin Delineation Guidelines
When to Create Separate Subbasins
Create separate subbasins when:
- Land use differs significantly: Urban vs. rural areas
- Soil types vary: Different hydrologic soil groups
- Outlets differ: Different discharge points in the system
- Calibration points exist: Gage locations within watershed
- Timing varies: Significantly different lag times
Subbasin Size Recommendations
| Analysis Type | Recommended Size |
|---|---|
| Urban drainage | 10 - 100 acres |
| Suburban development | 50 - 500 acres |
| Rural watersheds | 1 - 10 square miles |
| Regional planning | 10 - 100 square miles |
Lumping vs. Distributing
Lumped approach (fewer, larger subbasins):
- Simpler model
- Less data requirements
- May miss spatial variation
Distributed approach (more, smaller subbasins):
- Better spatial representation
- More data requirements
- More complex model maintenance
Next Steps
Now that you understand subbasin elements:
- Learn calibration: Model Calibration
- Calculate curve numbers: SCS Curve Number Calculator
- Review your model: Your First HEC-HMS Model
References
-
U.S. Army Corps of Engineers, Hydrologic Engineering Center. (2024). HEC-HMS User’s Manual. Davis, CA: USACE.
-
U.S. Army Corps of Engineers, Hydrologic Engineering Center. (2024). HEC-HMS Technical Reference Manual. Davis, CA: USACE.
-
Natural Resources Conservation Service. (1986). Urban Hydrology for Small Watersheds, TR-55. Washington, DC: USDA.
-
Natural Resources Conservation Service. (2004). National Engineering Handbook, Part 630: Hydrology. Washington, DC: USDA.
Summary
Subbasins are the hydrologic engine of HEC-HMS models. Key takeaways:
- Area defines the drainage area contributing to each subbasin
- Loss methods (SCS CN, Green-Ampt) determine how much precipitation infiltrates
- Transform methods (SCS UH, Clark) convert excess precipitation to a hydrograph
- Baseflow methods add groundwater contribution (optional for event modeling)
- Curve Number is the most important parameter for most analyses
Accurate subbasin configuration is essential for reliable HEC-HMS results. Take time to properly characterize land use, soils, and watershed geometry before moving to routing and simulation.