Understanding the Muskingum Method for Hydrologic Routing

The Muskingum method is a widely used hydrological routing technique for predicting the movement of flood waves through a river channel or reservoir. It relates the storage within a reach to both the inflow and outflow at that reach over a period of time.

Hydrological routing methods like the Muskingum method use the equation of continuity (mass balance) and a relationship between storage and flow within the reach. The continuity equation states that the difference between the inflow and outflow equals the change in storage over a time interval:

$$I - O = \frac{dS}{dt}$$

For application over a discrete time interval $\Delta t$, the continuity equation can be written as:

$$\frac{I_1 + I_2}{2} - \frac{O_1 + O_2}{2} = \frac{S_2 - S_1}{\Delta t}$$

Where $I_1$ and $O_1$ are inflow and outflow at the beginning of the interval (time step 1), and $I_2$ and $O_2$ are inflow and outflow at the end of the interval (time step 2).

The Muskingum Storage-Flow Relationship

The Muskingum method employs a specific linear relationship between the storage ($S$), inflow ($I$), and outflow ($O$) in a channel reach. This relationship is given by:

$$S = K [xI + (1-x)O]$$

Here:

• $K$ is a storage time constant with dimensions of time. It is approximately equal to the travel time of a flood wave through the reach.
• $x$ is a dimensionless weighting factor (ranging from 0 to 0.5) that represents the relative importance of inflow and outflow in determining storage. For reservoir-like storage, $x \approx 0$; for channel-like storage, $x$ is typically between 0.1 and 0.3.

Deriving the Muskingum Routing Equation

By substituting the Muskingum storage relationship into the discrete continuity equation and rearranging, we obtain the Muskingum routing equation, which allows calculating the outflow $O_2$ at the end of the time step based on known values $I_1, I_2, O_1$ and the reach parameters $K$ and $x$ (which determine the coefficients):

$$O_2 = C_0 I_2 + C_1 I_1 + C_2 O_1$$

The coefficients $C_0, C_1,$ and $C_2$ are functions of the parameters $K$, $x$, and the chosen time interval $\Delta t$. They are defined as:

• $$C_0 = \frac{\Delta t - 2Kx}{2K(1-x) + \Delta t}$$
• $$C_1 = \frac{\Delta t + 2Kx}{2K(1-x) + \Delta t}$$
• $$C_2 = \frac{2K(1-x) - \Delta t}{2K(1-x) + \Delta t}$$

Property of Muskingum Coefficients

A fundamental property of these Muskingum coefficients is that their sum is always equal to 1. Let's verify this:

$$C_0 + C_1 + C_2 = \frac{\Delta t - 2Kx}{2K(1-x) + \Delta t} + \frac{\Delta t + 2Kx}{2K(1-x) + \Delta t} + \frac{2K(1-x) - \Delta t}{2K(1-x) + \Delta t}$$

Combine the numerators over the common denominator:

$$C_0 + C_1 + C_2 = \frac{(\Delta t - 2Kx) + (\Delta t + 2Kx) + (2K(1-x) - \Delta t)}{2K(1-x) + \Delta t}$$

Simplify the numerator:

$$C_0 + C_1 + C_2 = \frac{\Delta t - 2Kx + \Delta t + 2Kx + 2K(1-x) - \Delta t}{2K(1-x) + \Delta t}$$

$$C_0 + C_1 + C_2 = \frac{(\Delta t + \Delta t - \Delta t) + (-2Kx + 2Kx) + 2K(1-x)}{2K(1-x) + \Delta t}$$

$$C_0 + C_1 + C_2 = \frac{\Delta t + 2K(1-x)}{2K(1-x) + \Delta t}$$

Since the numerator and denominator are the same, their ratio is 1:

$$C_0 + C_1 + C_2 = 1$$

This equation, $C_0 + C_1 + C_2 = 1$, is always satisfied by the coefficients of the Muskingum routing equation and represents a consistency requirement derived from the continuity principle.

Revision Table: Muskingum Method Key Points

Additional Information: Hydrologic Routing

Hydrologic routing simplifies the complex unsteady flow equations (like the Saint-Venant equations) by focusing on the continuity equation and a simplified storage-discharge relationship. This makes calculations easier, suitable for estimating flow changes over longer river reaches or reservoirs. The Muskingum method is a popular example, known for its simplicity and reasonable accuracy for many river systems. Other hydrologic routing methods exist, but the Muskingum method's linear storage assumption makes it computationally efficient.