Reach Routing Methods

Overview

Reach routing simulates the movement of a flood wave through a channel or conduit. As a hydrograph travels through a reach, two effects occur: translation (the hydrograph is delayed) and attenuation (the peak is reduced and the hydrograph is spread out). Different routing methods model these effects with varying levels of detail.

Method Comparison

Method	Translation	Attenuation	Parameters	Physical Basis	Computational Cost
Lag	Yes	No	Lag time	None	Very Low
Muskingum-Cunge	Yes	Yes	Channel geometry	Physics-based	Low
Modified Puls	Yes	Yes	Storage-discharge table	Empirical	Medium
Kinematic Wave	Yes	None (numerical only)	Channel geometry	Physics-based	Medium

Lag Method

The simplest routing approach. It shifts the inflow hydrograph forward in time by a fixed lag time, and no attenuation occurs: the outflow hydrograph has the same shape and peak as the inflow, only translated (shifted) in time.

Muskingum-Cunge Method

Derives Muskingum K and X parameters from physical channel properties, eliminating the need for calibration. This method is widely used for open channel routing when channel geometry is known. HydraLink uses the variable-parameter form: K and X are recomputed at every time step from the flow actually being routed, not fixed once for the reach.

Process

1. From the inflow peak, estimate the wave celerity and divide the reach into subreaches so each is traversed in about one time step
2. At each time step and subreach, solve Manning's equation for normal depth at the current flow using the reach cross-section (trapezoidal, rectangular, triangular, circular pipe, or an eight-point station–elevation section)
3. Compute the kinematic wave celerity c = dQ/dA from that section and K = Δx / c
4. Compute X = 0.5 × (1 − Q / (T × c × S × Δx)) where T = top width at normal depth, limited so the Muskingum coefficients stay non-negative
5. Route the subreach with the standard Muskingum coefficients

Modified Puls (Storage-Indication) Method

Uses a known storage-discharge relationship to route the hydrograph.

Storage Indication Curve

Routing

At each time step, interpolate O from the SI curve. This method is suitable for reaches where the storage-discharge relationship is known from field data or detailed hydraulic modeling. The table is entered as storage in acre-feet against discharge in cfs, and the reach can be split into subreaches.

Kinematic Wave Method

Assumes the momentum equation reduces to a balance between friction and gravity (kinematic approximation). Uses the same channel geometry as Muskingum-Cunge.

Best for steep, friction-dominated channels. HydraLink warns when the slope is below 0.002 ft/ft — any apparent attenuation there is numerical, not physical — and recommends Muskingum-Cunge instead. It also reports when the time step is poorly matched to the reach length (Courant number far from 1). Not suitable for backwater effects.

Selecting a Routing Method

The appropriate routing method depends on the site conditions, available data, and local design criteria. The following considerations may help guide the selection:

• Lag is appropriate for relatively short reaches where attenuation is negligible. It provides translation only (no peak reduction). The engineer should consider whether neglecting attenuation is appropriate for the given reach length and conditions.
• Muskingum-Cunge is commonly used when channel geometry is known but calibration data is unavailable. It derives routing parameters from measurable channel properties.
• Modified Puls works well when the storage-discharge relationship is known from field data or detailed hydraulic modeling (e.g., wide floodplain sections).
• Kinematic Wave suits steep channels where inertial and pressure gradient terms are small.

References

• Cunge, J.A. (1969). “On the Subject of a Flood Propagation Computation Method.” Journal of Hydraulic Research.
• USACE (2000). HEC-HMS Technical Reference Manual.
• Chow, V.T. (1959). Open-Channel Hydraulics. McGraw-Hill.