Flood routing
A flood wave changes as it travels. Passing through a reservoir or down a long river reach, its peak is lowered and delayed, because water goes into storage on the rising limb and is released on the falling limb. Flood routing calculates the outflow hydrograph from a known inflow hydrograph, so that engineers can size spillways, predict downstream flood levels and plan warnings.
All routing rests on the continuity equation: inflow minus outflow equals the rate of change of storage. What differs is how storage is related to flow. In a reservoir with a free spillway, storage depends only on outflow, through the water level, so level-pool routing, such as the modified Puls method, steps through time using the storage-elevation-discharge relation.
In a river reach, storage depends on both inflow and outflow, because the water surface is not level during a flood. The Muskingum method expresses storage as K times a weighted mix of inflow and outflow, with K close to the travel time through the reach and a weighting factor X between 0 and 0.5. Its routing equation calculates each new outflow from the current and previous inflow and the previous outflow, using three coefficients that add to one.
Hydraulic routing solves the full unsteady flow equations of Saint-Venant and is used for complex river systems, dam break studies and floodplain mapping, where backwater and overbank flow matter.
Check that C₀ + C₁ + C₂ = 1 and that C₀ is not negative; if it is, reduce K X or change Δt.
Continuity: I − O = dS / dtMuskingum: S = K [X I + (1 − X) O]O₂ = C₀ I₂ + C₁ I₁ + C₂ O₁, with D = K − K X + 0.5 ΔtC₀ = (−K X + 0.5 Δt) / D; C₁ = (K X + 0.5 Δt) / D; C₂ = (K − K X − 0.5 Δt) / D
Choosing a time step much longer than the travel time, which gives unstable or negative coefficients.
Problem. A reach has K = 12 h and X = 0.2. With Δt = 6 h, find the Muskingum coefficients. Then route one step: I₁ = 50, I₂ = 120 and O₁ = 50 m³/s.
- D = 12 − 2.4 + 3 = 12.6.
- C₀ = (−2.4 + 3) / 12.6 = 0.0476; C₁ = (2.4 + 3) / 12.6 = 0.4286; C₂ = (12 − 2.4 − 3) / 12.6 = 0.5238. Sum = 1.000.
- O₂ = 0.0476 × 120 + 0.4286 × 50 + 0.5238 × 50 = 5.71 + 21.43 + 26.19.
- O₂ = 53.3 m³/s: the outflow lags well behind the rising inflow.
Answer. C₀ = 0.048, C₁ = 0.429, C₂ = 0.524; O₂ ≈ 53.3 m³/s
Quick quiz
1What equation underlies all flood routing?
Continuity: inflow minus outflow equals rate of change of storage.
2What does K represent in Muskingum?
Approximately the travel time of the flood through the reach.
3Which method suits a reservoir with a free spillway?
Level-pool (modified Puls) routing.
- IS 11223 (India): Fixing spillway capacity: guidelines.
Cited by number and clause only. We do not reproduce clause text, tables or figures (BIS Act 2016, s.10(5) and s.11). Check the current edition and amendments before use.
References and further reading
Standard texts, listed in no order of merit. We do not copy text, figures, tables or solved problems from any book; our explanations and examples are original.
- K. Subramanya, Engineering Hydrology, McGraw Hill Education (India), latest edition.
- V. T. Chow, D. R. Maidment, L. W. Mays, Applied Hydrology, McGraw Hill, 1988.
Free and official sources
Found an error? Write to indiacongroupofengineers [at] gmail [dot] com. For corrections and feedback. We never ask for payment or personal details.