Civil, Structural & Geotechnical Engineering
Bagnold Stream Power Formula: Cross-Sectional Omega (W/m) & Unit Stream Power omega (W/m^2)
Quantify river geomorphic energy dissipation and channel incision/avulsion potential using Ralph A. Bagnold’s (1966, USGS PP 422-I) Cross-Sectional Stream Power Omega = rho * g * Q * S (W/m) and Specific Unit Stream Power omega = Omega / W (W/m^2).
Treat an alluvial river as a gravitational machine expending potential energy per meter and per square meter!
In his landmark USGS Professional Paper 422-I, desert and fluvial physicist Brigadier Ralph A. Bagnold (1966) reframed sediment transport and channel instability around thermodynamic power expenditure: as discharge Q (m^3/s) of water density rho (~= 1,000 kg/m^3) descends energy slope S (m/m) under gravity g = 9.80665 m/s^2, Cross-Sectional Stream Power Omega = rho * g * Q * S (Watts per meter of channel length, W/m) measures total gravitational power released per unit reach length!
Dividing Omega by active water-surface width W (m) yields Specific (Unit) Stream Power omega = Omega / W = (rho * g * Q * S) / W = tau_0 * V (Watts per square meter of river bed, W/m^2)! Across gravel- and cobble-bed rivers worldwide, Brookes (1988) and Magilligan (1992) established classic geomorphic stability benchmarks: channels with omega < 35 W/m^2 tend to remain stable or aggrade after straightening, whereas extreme floods exceeding omega >= 300 W/m^2 trigger catastrophic floodplain stripping, bank avulsion, and bedrock/alluvial incision!
| Specific Unit Stream Power omega (W/m^2) | Geomorphic Channel Response Regime | Engineering Design Implication |
|---|---|---|
| omega < 35 W/m^2 | Low-Energy Depositional / Stable | Prone to silt/sand aggradation if over-widened |
| 35 <= omega < 300 W/m^2 | Active Bedload & Meander Migration | Requires bank toe protection & equilibrium sizing |
| omega >= 300 W/m^2 | Catastrophic Flood Geomorphic Threshold | Major scour, avulsion, and boulder mobilization risk |
Bagnold Cross-Sectional Stream Power Omega (W/m) = 1000 * 9.80665 * Q * S
Evaluates 9806.65 * Q * S across river discharge Q (B2 m^3/s) and longitudinal energy slope S (B3 m/m).
=IF(OR(B2 <= 0, B3 <= 0), "Requires Discharge Q > 0 m^3/s and Slope S > 0", ROUND(1000 * 9.80665 * B2 * B3, 1))Specific Unit Stream Power omega (W/m^2) = Omega / Width & Stability Tier
Divides Omega by active channel top width W (B4 m) to output bed power flux in W/m^2.
=IF(OR(B2 <= 0, B3 <= 0, B4 <= 0), "Requires Positive Discharge Q, Slope S, and Channel Width W", ROUND((1000 * 9.80665 * B2 * B3) / B4, 2))Try the example
Use the original values or change them to test your understanding.