Go to ⇒  Presentation format


[Time-area Method]   [Clark Unit Hydrograph]   [Cascade of Linear Reservoirs]   [Kinematic Waves]   [Diffusion Waves]   [Assessment]   [Questions]   [Problems]   [References]     

CHAPTER 10:  
CATCHMENT ROUTING  

"We give a mathematical treatment of the 'competition' between kinematic and dynamic waves in river flow,
in order to show how completely the dynamic waves are subordinated in the case of greatest interest,
that is, when the speed of the river is well below critical."
M. J. Lighthill and G. B. Whitham (1955)


This chapter is divided into six sections. Section 10.1 describes the time-area method of hydrologic catchment routing. Section 10.2 describes the Clark unit hydrograph, a procedure closely related to the time-area method. Section 10.3 deals with the cascade of linear reservoirs, a widely accepted method of hydrologic catchment routing. Sections 10.4 and 10.5 describe two hydraulic methods of catchment routing, based on kinematic and diffusion waves, respectively. Section 10.6 contains a discussion of the capabilities and limitations of catchment routing techniques.


10.1  TIME-AREA METHOD

[Clark Unit Hydrograph]   [Cascade of Linear Reservoirs]   [Kinematic Waves]   [Diffusion Waves]   [Assessment]   [Questions]   [Problems]   [References]      [Top]  

Catchment Routing

Catchment routing refers to the calculation of flows in time and space within a catchment. The objective of catchment routing is to transform effective rainfall into streamflow. This is accomplished either in a lumped mode (e.g., time-area method) or in a distributed mode (e.g., kinematic wave method).

Methods for catchment routing are similar to those of reservoir and stream channel routing. In fact, many techniques used in reservoir and channel routing are also applicable to catchment routing. For instance, the concept of linear reservoir is used in both reservoir and catchment routing. Kinematic wave techniques were originaIly developed for river routing [9], but later were applied to catchment routing [20, 23].

Methods for catchment routing are of two types: (1) hydrologic and (2) hydraulic. Hydrologic methods are based on the storage concept and are spatially lumped to provide a runoff hydrograph at the catchment outlet. Examples of hydrologic catchment routing methods are the time-area method and the cascade of linear reservoirs. Hydraulic methods use kinematic or diffusion waves to simulate surface runoff within a catchment in a distributed context. Unlike hydrologic methods, hydraulic methods can provide runoff hydrographs inside the catchment.

Catchment routing models can use parametric, conceptual, andlor deterministic components. For instance, the hydrograph obtained by the time-area method can be routed through a linear reservoir using a storage constant derived by empirical (i.e., parametric) means. The cascade of linear reservoirs is a typical example of a conceptual model used in catchment routing. Kinematic and diffusion models are examples of deterministic methods used in catchment routing.

The concepts of translation and storage are central to the study of flow routing, whether in catchments, reservoirs, or stream channels. They are particularly important in catchment routing, because they can be studied separately. Translation may be interpreted as the movement of water in a direction parallel to the channel bottom. Storage may be interpreted as the movement of water in a direction perpendicular to the channel bottom. Translation is synonymous with runoff concentration; storage is synonymous with runoff diffusion.

In reservoir routing, storage is the primary mechanism, with translation almost nonexistent. In stream channel routing, the situation is reversed, with translation being the predominant mechanism and storage playing only a minor role. This is the reason why kinematic and diffusion waves are useful models of stream channel routing. In catchment routing, translation and storage are about equally important and, therefore, they are often accounted for separately. The translation effect can be related to runoff concentration, whereas the storage effect can be simulated with linear reservoirs.

Time-Area Method

The time-area method of hydrologic catchment routing transforms an effective storm hyetograph into a runoff hydrograph. The method accounts for translation only and does not include storage. Therefore, hydrographs calculated with the time-area method show a lack of diffusion, resulting in higher peaks than those that would have been obtained if storage had been taken into account. If necessary, the required amount of storage can be incorporated by routing the hydrograph obtained by the time-area method through a linear reservoir. The required amount of storage is determined by calibrating the linear reservoir storage constant K with measured data. Alternatively, suitable values of K can be estimated based on regionally derived formulas.

The time-area method is essentially an extension of the runoff concentration principle used in the rational method (Chapter 4). Unlike the rational method, however, the time-area method can account for the temporal variation of rainfall intensity. Therefore, the applicability of the time-area method is extended to midsize catchments.

The time-area method is based on the concept of time-area histogram, i.e., a histogram of contributing catchment subareas. To develop a time-area histogram, the catchment's time of concentration is divided into a number of equal time intervals. Cumulative time at the end of each time interval is used to divide the catchment into zones delimited by isochrone lines, i.e., the loci of points of equal travel time to the catchment outlet, as shown in Fig. 10-1 (a). For any point inside the catchment, the travel time refers to the time that it would take a parcel of water to travel from that point to the outlet. The catchment subareas delimited by the isochrones are measured and plotted in histogram form as shown in Fig. 10-1 (b).

Time-area method: (a) Isochrone delineation; (b) Time-area histogram

Figure 10-1  Time-area method: (a) Isochrone delineation; (b) Time-area histogram.

The time interval of the effective rainfall hyetograph must be equal to the time interval of the time-area histogram. The rationale of the time-area method is that, according to the runoff concentration principle (Section 2.4), the partial flow q at the end of each time interval is equal to the product of effective rainfall i times contributing subarea a, i.e., q = ia. The lagging and summation of the partial flows results in a runoff hydrograph for the given effective rainfall hyetograph and time-area histogram.

While the time-area method accounts for runoff concentration only, it has the advantage that the catchment shape is reflected in the time-area histogram and, therefore, in the runoff hydrograph. The procedure is illustrated by the following example.

 Example 10-1.

A 100-km2 catchment has a 4-h time of concentration, with isochrones at 1-h intervals resulting in the time-area histogram shown in Fig. 10-1 (b). A 6-h storm has the following effective rainfall hyetograph (Fig. 10-2):

Time (h) 0-1 1-2 2-3 3-4 4-5 5-6
Effective rainfall (cm/h) 0.5 1.0 2.0 1.5 1.0 0.5

Effective rainfall hyetograph

Figure 10-2  Effective rainfall hyetograph:  Example 10-1.

Use the time-area method to calculate the outflow hydrograph.


The routing is shown in Table 10-1.

  • Column 1 shows time in hours.

  • The flows shown in Cols. 2-7 were obtained by multiplying effective rainfall intensity times the confributing partial area. For instance, Col. 2 shows the contribution of the first effective rainfall interval (0.5 cm/h) on each of the subareas (10, 30, 20, and 40 km2). At t = 1 h, the partial flow due to the first effective rainfall interval is: 0.5 cm/h × 10 km2 = 5 km2-cm/h (i.e., the flow contributed by the subarea enclosed within the catchment outlet and the first isochrone takes 1 h to concentrate). Likewise, at t = 2 h, the partial flow due to the first effective rainfall interval is: 0.5 cm/hr × 30 km2 = 15 km2-cm/h (i.e., the flow contributed by the subarea enclosed within the first and second isochrones takes 2 h to concentrate at the catchment outlet).

  • The remaining values in Col. 2 (10 and 20) are calculated in a similar way. Finally, at t = 5 h, the flow is zero because it takes a full time interval (in the absence of runoff diffusion) for the last concentrated partial flow to recede back to zero.

  • Columns 2 to 7 show the partial flows contributed by the six effective rainfall intervals, each appropriately lagged a time interval (because the contribution of the second rainfall interval starts at t = 1, the third rainfall interval at t = 2, and so on).

  • The sum of these partial flows, shown in Col. 8, is the catchment outflow hydrograph.

  • In Col. 9, the hydrograph of Col. 8 is expressed in cubic meters per second (Col. 8 × 2.78).

  • The time base of the outflow hydrograph is 10 h, which is equal to the time of concentration (4 h) plus the effective rainfall duration (6 h).

  • To verify the accuracy of the computations, the sum of Col. 8 is 650 km2-cm/h, which represents 6.5 cm of effective rainfall depth uniformly distributed over the entire catchment area (100 km2). This value (6.5 cm) agrees with the total amount of effective rainfall.

Table 10-1  Time-Area Method of Catchment Routing:  Example 10-1.
(1) (2) (3) (4) (5) (6) (7) (8) (9)
Time
(h)
Partial flows (km2-cm/hr) Flow
(m3/s)
0.5
cm/h
1.0
cm/h
2.0
cm/h
1.5
cm/h
1.0
cm/h
0.5
cm/h
Sum
0 0 0 0.0
1 5 0 5 13.9
2 15 10 0 25 69.4
3 10 30 20 0 60 166.7
4 20 20 60 15 0 115 319.4
5 0 40 40 45 10 0 135 375.0
6 0 80 30 30 5 145 402.8
7 0 60 20 15 95 263.9
8 0 40 10 50 138.9
9 0 20 20 55.6
10 0 0 0.0
Total   650

calculator image 

ONLINE CALCULATION. Using ONLINE ROUTING06, the result is esentially the same as that of Col. 9, Table 10-1.


It is readily seen that the time-area method and the rational method (Chapter 4) share a common theoretical basis. However, since the time-area method uses effective rainfall and does not rely on runoff coefficients, it can account only for runoff concentration, with no provision for runoff diffusion. Diffusion can be provided by routing the hydrograph calculated by the time-area method through a linear reservoir with an appropriate storage constant.

Time of concentration. The time-area method leads to an alternate way of calculating time of concentration. Provided there is no runoff diffusion as would be the case of a hydrograph calculated by the time-area method, time of concentration can be calculated as the difference between hydrograph time base and effective rainfall duration. Intuitively, as rainfall ceases, the farthest parcels of water concentrate at the catchment outlet at a time equal to the time of concentration. Therefore:

tc = Tb - tr (10-1)

in which tc = time of concentration, Tb = time base of the translated-only hydrograph, and tr = effective rainfall duration.

Equation 10-1 can also be expressed in a slightly different form. Assuming that the point of inflection (i.e., the point of zero curvature) on the receding limb of a measured (i.e., translated and diffused) hydrograph coincides with the end of the translated-only hydrograph, the time to point of inflection of the measured hydrograph can be used in Eq. 10-1 in lieu of time base. Therefore time of concentration can be defined as the difference between the time to point of inflection I and the effective rainfall duration (see Fig. 10-3):

tc = ti - tr (10-2)

in which ti = time to point of inflection on the receding limb of a measured hydrograph. The advantage of Eq. 10-2 over Eq. 10-1 is that, unlike the time base of the translated-only hydrograph, the point of inflection on the receding limb of a measured hydrograph can be readily ascertained.

Alternate definition of time of concentration

Figure 10-3  Alternate definition of time of concentration.


10.2  CLARK UNIT HYDROGRAPH

[Cascade of Linear Reservoirs]   [Kinematic Waves]   [Diffusion Waves]   [Assessment]   [Questions]   [Problems]   [References]      [Top]   [Time-area Method]  

The procedure to derive a Clark unit hydrograph parallels that of the time-area method [2]. First, it is necessary to determine the catchment isochrones. In the Clark method, however, a unit effective rainfall is used in lieu of the effective storm hyetograph used in the time-area method. This leads to an outflow hydrograph corresponding to a unit runoff depth, that is, a unit hydrograph. Since the unit hydrograph calculated in this way lacks (runoff) diffusion, Clark suggested that it be routed through a linear reservoir.

As with the time-area method, an estimate of the linear reservoir storage constant is required. This can be obtained either from the tail of a measured hydrograph or by using a regionally derived formula. In the latter case, the Clark unit hydrograph can be properly regarded as a synthetic unit hydrograph.

Like the time-area method, the Clark unit hydrograph method has the advantage that the catchment's properties (shape, hydraulic length, surface roughness, and so on) are reflected in the time-area histogram and, therefore, on the shape of the unit hydrograph. This feature has contributed to the popularity of the Clark unit hydrograph in engineering practice [7].

When using the Clark or time-area methods, the storage constant can be estimated from the tail of a measured hydrograph. For this purpose, the differential equation of storage (Eq. 8-4) is evaluated at a time for which inflow equals zero (I = 0), i.e., past the end of the translated-only hydrograph. Alternatively, it can be evaluated at the point of inflection on the receding limb of a measured hydrograph (Fig. 10-3). This leads to:

             dS
- O  =  _____
             dt
(10-3)

and since S = KO, the following expression for K is obtained:

                O
K  =  - ________
               dO
             _____
               dt
(10-4)

in which O and dO/dt are evaluated past the end of the translated-only hydrograph or at (the time to) the point of inflection on the receding limb of a measured hydrograph.

The derivation of the Clark unit hydrograph is illustrated by the following example.

 Example 10-2.

Use the Clark method to derive a 2-h unit hydrograph for the catchment of Example 10-1. To provide storage, route the translated-only hydrograph through a linear reservoir of storage constant K = 2 h. Use Δt = 1 h.


The 2-h unit hydrograph has an effective rainfall intensity of 0.5 cm/h (i.e., 1-cm depth distributed over a 2-h duration). The calculations are shown in Table 10-2.

  • Column 1 shows time in hours.

  • Column 2 shows the contribution of the first hour, with 0.5 cm/h of effective rainfall. The procedure is the same as in Table 10-1, Col. 2.

  • Column 3 shows the contribution of the second hour, with 0.5 cm/h of effective rainfall. Again, the procedure is the same as in Table 10-1, Col. 2; but the partial flows are lagged 1 h.

  • The translated-only unit hydrograph shown in Col. 4 is the sum of the partial flows (Cols. 2 and 3). The translated-only unit hydrograph (Col. 4) is the inflow to the linear reservoir. With Δt/K = 1/2, the routing coefficients (Table 8-1) are C0 = 1/5, C1 = 1/5, and C2 = 3/5.

  • The partial flows of the linear reservoir routing are shown in Cols. 5 to 7, and the translated-and-diffused unit hydrograph (in km2-cm/h) shown in Col. 8 is the sum of Cols. 5 to 7 (See Example 8-1 for details of the linear reservoir routing procedure).

  • Column 9 shows the translated-and-diffused Clark unit hydrograph in cubic meters per second.

  • The sum of the ordinates of the translated-only hydrograph (Col. 4) is 100, which amounts to 1 cm of effective rainfall depth uniformly distributed over 100 km2 of catchment area. Likewise, the sum of the ordinates of the translated-and-diffused hydrograph (Col. 8) is 99.98, which verifies not only that the calculated hydrograph is a unit hydrograph but also that the calculation is mass (i.e., volume) conservative.

Note that the peak of the translated-only unit hydrograph (Col. 4) is 30 km2-cm/h, whereas the peak of the translated-and-diffused unit hydrograph (Col. 8) is 21.05 km2-cm/h. Also, note that the time base of the translated-only unit hydrograph ends sharply at 6 h, whereas the time base of the translated-and-diffused unit hydrograph is much longer, with the receding limb of the unit hydrograph gradually approaching zero. This reveals the substantial amount of runoff diffusion provided by the linear reservoir.

Table 10-2  Derivation of Clark Unit Hydrograph:  Example 10-2.
(1) (2) (3) (4) (5) (6) (7) (8) (9)
Time
(h)
Partial Flows and Sum (km2-cm/hr) Flow
(m3/s)
0.5
cm/h
0.5
cm/h
Sum C0 I2 C1 I1 C2 O1 Sum
0 0 0 - - - 0.00 0.00
1 5 0 5 1 0 0 1.00 2.78
2 15 5 20 4 1 0.60 5.60 15.55
3 10 15 25 5 4 3.36 12.36 34.33
4 20 10 30 6 5 7.42 18.42 51.17
5 0 20 20 4 6 11.05 21.05 58.47
6 0 0 0 0 4 12.63 16.63 46.19
7 0 0 0 0 0 9.98 9.98 27.72
8 0 0 0 0 0 5.99 5.99 16.64
9 3.59 3.59 9.97
10 2.15 2.15 5.97
11 1.29 1.29 3.58
12 0.78 0.78 2.17
13 0.46 0.46 1.28
14 0.28 0.28 0.78
15 0.17 0.17 0.47
16 0.10 0.10 0.28
17 0.06 0.06 0.17
18 0.04 0.04 0.11
19 0.02 0.02 0.06
20 0.01 0.01 0.03
Total 100 99.98

calculator image 

ONLINE CALCULATION. Using ONLINE ROUTING07, the result is esentially the same as that of Col. 9, Table 10-2.


Important Note. In the Clark method, the number of partial flows (two in Example 10-2, shown as Cols. 2 and 3 of Table 10-2) is equal to the duration of the unit hydrograph (2 h) divided by the (chosen) time interval of histogram definition (1 h). Also, in Example 10-2, the unit effective rainfall intensity (0.5 cm/h) is equal to the unit rainfall depth (1 cm) divided by the duration of the unit hydrograph (2 h).

By using Eq. 10-4, the linear reservoir storage constant can be calculated directly from the tail of a measured hydrograph. To illustrate the procedure, in Table 10-2, Col. 9, the two lines for t = 6 h and t = 7 h show zero outflow in the translated-only unit hydrograph (Col. 4), that is, zero inflow to the linear reservoir. Therefore, Eq. 10-4 can be applied between t = 6 h and t = 7 h.

The average outflow (Col. 9) is: (46.19 + 27.72)/2 = 36.955 m3/s. The rate of change of outflow is: (27.72 - 46.19)/(1 h) = -18.47 (m3/s)/h. Therefore, the storage constant (Eq. 10-4) is:

K = - (36.955) / (-18.47) = 2 h

Likewise, between t = 7 h and t = 8 h:

K = - [(27.72 + 16.64) / 2] / [(16.64 - 27.72) / (1)] = 2 h

In other words, Eq. 10-4 applies at the tail of the outflow hydrograph, after the translated-only (inflow) hydrograph has receded back to zero. When using the Clark (or time-area) method, the time base of the translated-only hydrograph is equal to the sum of concentration time plus the unit hydrograph (or effective storm) duration (See Eq. 10-1).

With the help of regional analysis (Chapter 7), the Clark parameters (time of concentration and linear reservoir storage constant) can be estimated based on catchment characteristics. This effectively qualifies the Clark unit hydrograph as a synthetic unit hydrograph. The Eaton [4], O'Kelly [11], and Cordery [3] models are examples of this approach. Singh [18] has a review of Clark synthetic unit hydrograph parameters.


10.3  CASCADE OF LINEAR RESERVOIRS

[Kinematic Waves]   [Diffusion Waves]   [Assessment]   [Questions]   [Problems]   [References]      [Top]   [Time-area Method]   [Clark Unit Hydrograph]  

As seen in Section 8.2, a linear reservoir has a diffusion effect on the inflow hydrograph. If an inflow hydrograph is routed through a linear reservoir, the outflow hydrograph has a reduced peak and an increased time base. This increase in time base causes a difference in the relative timing of inflow and outflow hydrographs, referred to as the lag. The amount of diffusion (and associated lag) is a function of the ratio Δt/K, a larger diffusion effect corresponding to smaller values of Δt/K.

The cascade of linear reservoirs is a widely used method of hydrologic catchment routing. As its name implies, the method is based on the connection of several linear reservoirs in series. For N such reservoirs, the outflow from the first would be taken as inflow to the second, the outflow from the second as inflow to the third, and so on, until the outflow from the (N - 1)th reservoir, is taken as inflow to the N th reservoir. The outflow from the N th reservoir is taken as the outflow from the cascade of linear reservoirs. Admittedly, the cascade of reservoirs to simulate catchment response is an abstract concept; nevertheless, it has proven to be quite useful in practice.

Each reservoir in the series provides a certain amount of diffusion and associated lag. For a given set of parameters Δt/K and N, the outflow from the last reservoir is a function of the inflow to the first reservoir. In this way, a one-parameter linear reservoir method (Δt/K) is extended to a two-parameter catchment routing method. Moreover, the basic routing formula (Eq. 8-15) and routing coefficients (Eqs. 8-16 to 8-18) remain essentially the same.

The addition of the second parameter (N) provides considerable flexibility in simulating a wide range of diffusion and associated lag effects. However, the conceptual basis of the method restricts its general use, since no direct relation between either of the parameters to the physical problem can be readily envisaged. Notwithstanding this apparent limitation, the method has been widely used in catchment simulation, primarily in applications involving large gaged river basins. Rainfall-runoff data can be used to calibrate the method, i.e., to determine a set of parameters Δt/K and N that produces the best fit to the measured data.

The analytical version of the cascade of linear reservoirs is referred to as the Nash model [10]. The numerical version is featured in several hydrologic simulation models developed in the United States and other countries. Notable among them is the SSARR model (Chapter 13), which uses it in its watershed, stream channel routing, and baseflow modules [19]. To derive the routing equation for the method of cascade of linear reservoirs, Eq. 8-15 is reproduced here in a slightly different form:

Q j+1 n+1 = C0 Q j n+1 + C1 Q j n + C2 Q j+1 n (10-5)

in which Q represents discharge, whether inflow or outflow and j and n are space and time indexes, respectively (Fig. 10-4).

Space-time discretization in the method of cascade of linear reservoirs

Figure 10-4  Space-time discretization in the method of cascade of linear reservoirs.

As with Eq. 8-15, the routing coefficients C0, C1 and C2 are a function of the dimensionless ratio Δt /K. This ratio is properly a Courant number (C = Δt /K). In terms of Courant number, Eqs. 8-16 to 8-18 are expressed as follows:

               C
C0  =  _______
            2 + C
(10-6)

C1  =  C0 (10-7)

            2 - C
C2  =  _______
            2 + C
(10-8)

For application to catchment routing, it is convenient to define the average inflow as follows:

            Q j n  +  Q j n+1
j = __________________
                      2
(10-9)

Substituting Eq. 10-6 to 10-9 into Eq. 10-5 gives the following:

Q j+1n+1 = 2 C1 j  +  C2 Q j+1n (10-10)

or, alternatively, through algebraic manipulation:

                       2 C
Q j+1n+1 =   _________ [ j  -  Q j+1n ]  +  Q j+1n
                     2 + C
(10-11)

Equation 10-11 is the routing equation of the SSARR model [20]. Equations 10-10 and 10-11 are in a form convenient for catchment routing because the inflow is usually a rainfall hyetograph, that is, a constant average value per time interval.

Smaller values of C lead to greater amounts of runoff diffusion. For values of C > 2, the behavior of Eq. 10-10 (or Eq. 10-11) is highly dependent on the type of input. For instance, in the case of a unit impulse (rainfall duration equal to the time interval), Eq. 10-10 (or Eq. 10-11) results in negative outflow values (numerical instability). For this reason, Eq. 10-10 (or Eq. 10-11) are restricted in practice to C ≤ 2.

The method of cascade of linear reservoirs is illustrated by the following example.

 Example 10-3.

Use the method of cascade of linear reservoirs to route the following effective storm hyetograph for a 1000-km2 basin. Use N = 3, Δt = 6 h and K = 12 h.

Time (h) 0-6 6-12 12-18 18-24
Effective rainfall (cm/h) 0.2 1.0 0.8 0.4
Effective rainfall (cm) 1.2 6.0 4.8 2.4


The Courant number is C = Δt/K = 6/12 = 1/2, which results in 2C1 = 2/5 and C2 = 3/5. The computations are shown in Table 10-3.

  • Column 1 shows time in hours.

  • Column 2 shows the inflow to the first reservoir (in km2-cm/h) calculated by multiplying each one of the effective rainfall intensities (0.2, 1.0, 0.8, and 0.4 cm/h) times the basin area (1000 km2).

  • Column 3 is the outflow from the first reservoir.

  • Columns 4 and 5 are the inflow and outflow for the second reservoir.

  • Columns 6 and 7 are the inflow and outflow for the third reservoir.

  • Column 8 is the result of the routing, i.e., the outflow from the third reservoir, in m3/s.

To illustrate the calculations for the first reservoir, following Eq. 10-10, 2/5 of the average inflow for the first time interval [(2/5) × 200 km2-cm/h] plus 3/5 of the outflow at time t = 0 h [(3/5) × 0 km2-cm/h] is equal to the outflow at 6 h: 80 km2cm/h. Likewise, 2/5 of the average inflow for the second time interval [(2/5) × 1000 km2-cm/hr] plus 3/5 of the outflow at time t = 6 h [(3/5) × 80 km2-cm/h] is equal to the outflow at 12 h: 448 km2-cm/ h, and so on.

The (average) inflow to the second reservoir (Col. 4) is the average outflow from the first reservoir (Col. 3). For instance, for the first time interval, 40 km2-cm/h is the average of 0 and 80 km2-cm/h.

The calculations proceed in a recursive fashion until the routing through the three linear reservoirs has been completed. Note that the sum of Cols. 3, 5 and 7 is approximately the same: 2400 km2-cm/hr. Since the time interval is 6 h, this is equivalent to 2400 × 6/1000 = 14.4 cm of effective rainfall depth uniformly distributed over 1000 km2 of basin area.

Note that the peak outflow from the first reservoir is 588.8 km2-cm/h, and it occurs at 18 h; the peak outflow from the second reservoir is 396.03 km2-cm/h, occurring at 30 h; and the peak outflow from the third reservoir is 308.61 km2-cm/h, occurring at 36 h. This shows that the effect of the cascade is to produce a certain amount of runoff diffusion at every step, with a corresponding increase in the lag of catchment response.

Table 10-3  Cascade of Linear Reservoirs:  Example 10-3.
(1) (2) (3) (4) (5) (6) (7) (8)
Time
(h)
Flow (km2-cm/hr)Flow
(m3/s)
N = 1 N = 2 N = 3
Inflow Outflow Inflow Outflow Inflow Outflow
0 0.00 0.00 0.00 0.00
200 40.00 8.00
6 80.00 16.00 3.20 8.89
1000 264.00 65.60
12 448.00 115.20 28.16 78.22
800 518.40 195.84
18 588.80 276.48 95.23 264.53
400 551.04 331.39
24 513.28 386.30 189.69 526.93
419.62 391.16
30 307.97 396.03 270.28 750.77
246.38 366.10
36 184.78 336.17 308.61 857.25
147.83 298.50
42 110.87 260.83 304.57 846.02
88.69 226.40
48 66.52 191.98 273.30 759.57
53.21 164.23
54 39.91 136.47 229.67 637.96
31.93 115.56
60 23.95 94.65 184.03 511.18
19.16 79.55
66 14.37 64.45 142.24 395.10
11.50 53.86
72 8.62 43.27 106.89 296.91
6.89 36.00
78 5.17 28.72 78.53 218.15
4.14 23.80
84 3.10 18.89 56.64 157.33
2.48 15.61
90 1.86 12.33 40.23 111.74
1.49 10.16
96 1.12 7.99 28.20 78.34
0.89 6.57
102 0.67 5.15 19.55 54.30
0.53 4.22
108 0.40 3.30 13.42 37.27
0.32 2.70
114 0.24 2.11 9.13 25.37
0.19 1.73
120 0.14 1.34 6.17 17.14
0.11 1.09
126 0.08 0.85 4.14 11.50
0.08 0.69
132 0.05 0.53 2.76 7.67
0.04 0.43
138 0.03 0.33 1.83 5.08
0.02 0.27
144 0.02 0.21 1.21 3.35
0.01 0.17
150 0.01 0.13 0.79 2.19
0.00 0.10
156 0.00 0.06 0.52 1.44
0.00 0.05
162 0.00 0.04 0.33 0.92
0.00 0.03
168 0.00 0.02 0.21 0.58
Sum 2400 2399.83 2399.83 2399.53

calculator image 

ONLINE CALCULATION. Using ONLINE ROUTING08, the result is esentially the same as that of Col. 8, Table 10-3.


The cascade of linear reservoirs provides a convenient mechanism for simulating a wide range of catchment routing problems. Furthermore, the method can be applied to each runoff component (surface runoff, subsurface runoff, and baseflow) separately, and the catchment response can be taken as the sum of the responses of the individual components.

For instance, assume that a certain basin has 10 cm of runoff, of which 7 cm are surface runoff, 2 cm are subsurface runoff, and 1 cm is baseflow. Since surface runoff is the less diffused process, it can be simulated with a high Courant number, say C = 1, and a small number of reservoirs, say N = 3. Subsurface runoff is much more diffused than surface runoff; therefore, it can be simulated with C = 0.4 and N = 5. Baseflow, being very diffused, can be simulated with C = 0.1 and N = 7. In practice, the parameters C and N are determined by extensive calibration. In this sense, the cascade of linear reservoirs remains essentially a conceptual model [17].


10.4  KINEMATIC WAVES

[Diffusion Waves]   [Assessment]   [Questions]   [Problems]   [References]      [Top]   [Time-area Method]   [Clark Unit Hydrograph]   [Cascade of Linear Reservoirs]  

Hydraulic catchment routing using kinematic waves was introduced by Wooding in 1965 [21, 22, 23]. Since then, the kinematic wave approach has been widely used in deterministic catchment modeling. The approach can be either lumped or distributed, depending on whether the parameters are kept constant or allowed to vary in space. Analytical solutions are suited to lumped modeling, whereas numerical solutions lend themselves to distributed modeling.

Wooding used an open-book geometric configuration (Fig. 4-15) to represent the catchment-stream problem physically. As its name implies, an open-book configuration consists of two rectangular catchments separated by a stream and draining laterally into it; in turn the streamflow drains out of the catchment outlet. Wooding used analytical solutions of kinematic waves and the method of characteristics to formulate his method. Since diffusion is absent from these solutions, the method is strictly applicable only to kinematic waves. Criteria for the applicability of kinematic waves have been developed by Woolhiser and Liggett [24] (Eq. 4-55) for overland flow, and by Ponce et al. [12] for stream channel flow (Eq. 9-44).

Kinematic catchment routing models can be approached in a variety of ways. Methods can be either (1) analytical or numerical, (2) lumped or distributed, (3) linear or nonlinear, or (4) single plane, two-plane, or a cascade of planes [7, 8]. Analytical models take advantage of the nondiffusive properties of kinematic waves, whereas numerical models are usually based on the method of finite differences. Linear models assume a constant wave celerity, but nonlinear models relax this restriction. The feature of variable wave celerity often renders the nonlinear models impractical because of wave steepening and associated kinematic shock development [8, 13]. Single- and two-plane models are used in hydrologic engineering practice [7].

The application of kinematic wave modeling to catchment routing is illustrated here with an example of a two-plane finite difference numerical model. The model could be either lumped or distributed, depending on whether the inputs and parameters are allowed to vary in space or not. For simplicity, this example considers constant input (i.e., constant effective rainfall) and constant parameters (i.e., a linear mode of computation). In practice, a computer-aided solution may relax this restriction.

Two-Plane Linear Kinematic Catchment Routing Model

Assume a catchment configured as two rectangular planes adjacent to each other, draining laterally into a stream channel located between them. Each of the planes is 100 m long by 200 m wide, and the channel is 200 m long (Fig. 10-5). The bottom friction in the planes and channel is such that the average velocity in the planes is 0.0417 m/s and the average velocity in the channel is 0.3 m/s. It is desired to obtain the hydrograph at the catchment outlet resulting from an effective rainfall of 9 cm/h lasting 20 min.

Two-plane linear kinematic catchment routing model

Figure 10-5  Two-plane linear kinematic catchment routing model.

Calculation of flow parameters. Since the model is linear, it is first necessary to calculate the flow parameters on which to base the calculation of the routing parameters and coefficients.

The flow per unit width in the midlength of each plane is equal to the effective rainfall intensity times the contributing area (50 m × 1 m):

             9 cm/h × 50 m × 0.01 m/cm
qp  =  ________________________________  =  0.00125 m2/s
                             3600 s/h
(10-12)

Since the average velocity in the planes is vp = 0.0417 m/s, the average flow depth in the planes is:

dp = qp/vp = 0.00125 / 0.0417 = 0.03 m

Laminar flow in the planes is assumed, with discharge-depth rating exponent βp = 3. Therefore, the wave celerity in the planes is:

cp = βp vp = 3 × 0.0417 = 0.125 m/s

The flow in the midlength of the channel is equal to the effective rainfall intensity times the contributing area (2 planes × 100 m × 100 m):

             9 cm/h × 2 × 100 m × 100 m × 0.01 m/cm
Qc  =  ____________________________________________  =  0.5 m3/s
                                         3600 s/h
(10-13)

Assume a channel top width Tc = 5 m. Therefore, the flow per unit width in the channel is:

qc = Qc /Tc = 0.1 m2/s

Since the average velocity in the channel is vc = 0.3 m/s, the average flow depth in the channel (at midlength) is:

dc = qc /vc = 0.333 m

A wide channel and turbulent Manning friction is assumed, with discharge-area rating exponent βc = 1.67. Therefore, the wave celerity in the channel is:

cc = βc vc = 1.67 × 0.3 = 0.5 m/s

The time of conncentration equal to the travel time in the planes plus the travel time in the channel. The travel time in the planes is: (100 m)/(0.125 m/s) = 800 s. The travel time in the channel is: (200 m)/(0.5 m/s) = 400 s. Therefore, the time of concentration is: 800 + 400 = 1200 s, which is equal to the effective rainfall duration. This assures concentrated flow at the catchment outlet.

The maximum possible (i.e., equilibrium) peak flow is equal to the product of rainfall intensity and catchment area:

             9 cm/h × 2 × 100 m × 200 m × 0.01 m/cm
Qc  =  ____________________________________________  =  1 m3/s
                                       3600 s/h
(10-14)

The total volume of runoff is:

             9 cm/h × 20 min × 2 × 100 m × 200 m
Qc  =  _________________________________________  =  1200 m3
                             100 cm/m × 60 min/h
(10-15)

Selection of discrete intervals. For simplicity, a space interval Δx = 100 m is chosen for the planes. This amounts to one spatial increment in the planes. In an actual application using a computer, a smaller value of Δx would be indicated. The time interval is chosen as Δt = 10 min (the rainfall duration is tr = 20 min). This leads to a Courant number in the planes:

Cp = cptx) = 0.75

In the case of the channel, a space interval Δy = 200 m is chosen, that is, one spatial increment in the channel. This leads to a Courant number in the channel:

Cc = cctx) = 1.5

Selection of routing scheme. There are many possible choices for routing scheme. Either first- or second-order schemes may be used (Section 9.2). In practice, first-order schemes are preferred because they are more stable than second-order schemes (compare the results of a first order scheme, Table 9-5, with those of a second order scheme, Table 9-4).

Two first-order schemes are chosen here: (1) Scheme I, forward-in-time, backward-in-space, stable for Courant numbers C ≤ 1 (similar to the convex method, see Example 9-5), and (2) Scheme II, forward-in-space, backward-in-time, stable for Courant numbers C ≥ 1 (exact opposite of the convex method). The use of these two schemes guarantees that the solution will remain stable because scheme I is used for Courant numbers C ≤ 1, whereas scheme II is used for Courant numbers C > 1 [7]. In the present application, scheme I is used for routing in the planes (Cp = 0.75), and scheme II is used for routing in the channel (Cc = 1.5).

Lateral inflows are an integral part of catchment routing. For routing in the planes, lateral inflow is the effective rainfall; for channel routing, lateral inflow is the lateral contribution from the planes. Therefore, it is necessary to discretize the kinematic wave equation with lateral inflow, Eq. 9-43.

The discretization of Eq. 9-43 in a forward-in-time, backward-in-space linear scheme, as shown in Fig. 10-6 (a), leads to:

Q j+1 n+1 = C1 Q j n + C2 Q j+1n + C3 QL (10-16)

in which

C1 = C (10-16a)

C2 = 1 - C (10-16b)

C3 = C (10-16c)

with the Courant number C = βv Δts, with Δs either Δx (planes) or Δy (channel). The term QL is the lateral inflow in cubic meters per second. For routing in the planes, the lateral inflow is equal to the effective rainfall (centimeters per hour) times the applicable area (square meters). For channel routing, the lateral inflow is the average distributed lateral inflow (cubic meters per second per meter) multiplied by the channel length (meters).

Space-time discretization of first-order schemes of kinematic wave equation

Figure 10-6  Space-time discretization of first-order schemes of kinematic wave equation with lateral inflow:
(a) forward-in-time, backward-in-space; (b) forward-in-space, backward-in-time.

The discretization of Eq. 9-43 in a forward-in-space, backward-in-time linear scheme, as shown in Fig. 10-6 (b), leads to:

Q j+1 n+1 = C0 Q j n+1  +  C1 Q j n  +  C3 QL (10-17)

in which

C0 = (C - 1) /C (10-17a)

C1 = 1/C (10-17a)

C3 = 1 (10-17a)

The catchment routing is shown in Table 10-4.

Table 10-4  Two-plane Catchment Routing With Kinematic Wave Technique.
(1) (2) (3) (4) (5) (6) (7)
Time
(min)
Flow (m3/s)
QLp Ip Op QLc Ic Oc
0 0 0 0 0
0.5 0.3750
10 0 0.3750 0 0.3750
0.5 0.8437
20 0 0.4687 0 0.8437
0 0.5859
30 0 0.1172 0 0.5859
0 0.1465
40 0 0.0293 0 0.1465
0 0.0366
50 0 0.0073 0 0.0366
0 0.0091
60 0 0.0018 0 0.0091
0 0.0023
70 0 0.0005 0.0023
0 0.0006
80 0 0.0001 0.0006
0 0.0001
90 0 0 0.0001
Sum 1.00 0.9999 1.9998

Explanation of Table 10-4.

  • Column 1 shows time in minutes,

  • Cols. 2-4 show the plane routing, and Cols. 5-7 show the channel routing.

  • Column 2 shows the lateral inflow to each plane, calculated as follows:

                 9 cm/h × 100 m × 200 m × 0.01 m/cm
    QLp  =  _________________________________________  =  0.5 m3/s
                                         3600 s/h
    (10-18)

    Note that the lateral inflow is an average value for a time interval, and it lasts 20 min (i.e., the effective rainfall duration).

  • Column 3 shows the upstream inflow to the plane, that is, zero (this example does not consider upstream inflow to the planes).

  • Column 4 is obtained by routing with Eq. 10-16, with Courant number in the planes Cp = 0.75. Column 4 is the outflow hydrograph from each plane. The first value of Col. 4 (for t = 0) is: 0; the second value is: (0.75 × 0) + (0.25 × 0) + (0.75 × 0.5) = 0.3750; the third value is: (0.75 × 0) + (0.25 × 0.3750) + (0.75 × 0.5) = 0.4687; and so on.

  • The sum of Col. 2 is 1.0; likewise, the sum of Col. 4 is 0.9999, which confirms that the volume under the outflow hydrograph from each plane is 1 m3/s × 10 min × 60 s/min = 600 m3.

  • Column 5 shows the average lateral inflow to the channel, obtained by multiplying by 2 the average lateral inflow from each plane (Col. 4) (to account for two planes of the same dimensions). The first value of Col. 5 is: [(0 + 0.3750)/2] × 2 = 0.3750; the second value is: [(0.3750 + 0.4687)/2] × 2 = 0.8437; and so on.

  • Column 6 is the upstream inflow to the channel, that is, zero (this example does not consider upstream inflow to the channel).

  • Column 7 is obtained by routing with Eq. 10-17, with Courant number in the channel Cc = 1.5. The first value of Col. 7 (for t = 0) is: 0; the second value is: (0.33 × 0) + (0.67 × 0) + (1.00 × 0.3750) = 0.3750; the third value is: (0.33 × 0) + (0.67 × 0) + (1.00 × 0.8437) = 0.8437; and so on.

  • The sum of Col. 7 is 1.9998, confirming that the total runoff volume is:

    1.9998 m3/s × 10 min × 60 s/min = 1200 m3.

  • Column 7 is the outflow hydrograph from the catchment.

Assessment of kinematic wave method. The calculated outflow hydrograph peak is 0.8437 m3/s, and it occurs at 20 min. This value is less than the maximum peak flow, Eq. 10-14: 1 m3/s. Since the rainfall duration is equal to the time of concentration this implies that the hydrograph has undergone a certain amount of runoff diffusion. This diffusion is really numerical diffusion, due primarily to the coarse grid size and secondarily to the Courant numbers (of planes and channel) being different than 1.

To prove this assertion, it is necessary to reduce the grid size and test the convergence of the kinematic wave schemes, Eqs. 10-16 and 10-17. Convergence refers to the ability of the numerical scheme to approach the analytical solution as the grid is refined. Due to the large number of calculations involved, the procedure is better accomplished with the aid of a computer program. Table 10-5 shows the results obtained by successive grid refinement, using program EH1000B.

Table 10-5  Kinematic Wave Catchment Routing:
Effect of Grid Resolution.

Number of
Increments
Δx
(m)
Δy
(m)
Δt
(s)
Peak Flow
(m3/s)
1 100 200 600 0.8437
2 50 100 300 0.9063
4 25 50 150 0.9490
8 12.5 25 75 0.9776
16 6.25 12.5 37.5 0.9899

It is seen that the results are a function of grid size and that the peak flow value converges to the maximum possible value (1 m3/s) as the grid is refined.

It is concluded that a kinematic wave numerical solution for catchment routing is grid dependent. If necessary, numerical diffusion can be eliminated by successive grid refinement (while keeping the Courant number as close to 1 as possible). However, in this case the calculated hydrograph would be translated only, with no diffusion. This may be adequate for catchments with negligible runoff diffusion (e.g., small catchments with slopes on the order of 1% or more), but is generally not adequate for catchments showing substantial amounts of runoff diffusion (e.g., midsize catchments of mild relief). For the latter, the diffusion wave technique may be used as a viable alternative to the kinematic wave.


10.5  DIFFUSION WAVES

[Assessment]   [Questions]   [Problems]   [References]      [Top]   [Time-area Method]   [Clark Unit Hydrograph]   [Cascade of Linear Reservoirs]   [Kinematic Waves]  

Catchment routing with diffusion waves is applicable to cases where both translation and diffusion are important, that is, for routing in midsize catchments where catchment slope is such that the kinematic wave criterion is not satisfied. Although the concept of diffusion waves and catchment routing dates back to the work of Dooge [1], actual numerical applications are relatively more recent [14, 16]. Diffusion wave routing can provide grid-independent results, and is therefore regarded as an improvement over grid-dependent techniques.

The diffusion wave catchment routing approach is illustrated here by using the same example as in the previous section. The Muskingum-Cunge method (Chapter 9) is used as the routing scheme of the diffusion wave method [14].

Two-Plane Linear Diffusion Catchment Routing Model

This example is similar to that of the previous section. Assume a catchment configured as two rectangular planes adjacent to each other, draining laterally into a stream channel located between them. Each of the planes is 100 m long by 200 m wide, and the channel is 200 m long (Fig. 10-5). The slopes of planes and channel (in the direction of the flow) are Sop = 0.01; and Soc = 0.01, respectively. The bottom friction in the planes and channel is such that the average velocity in the planes is vp = 0.0417 m/s, and the average velocity in the channel is vc = 0.3 m/s. It is desired to calculate the runoff hydrograph at the catchment outlet resulting from an effective rainfall of 9 cm/h lasting 20 min.

Calculation of flow parameters. Since the model is linear, it is first necessary to calculate the flow parameters on which to base the calculation of the routing parameters and coefficients. As described in the preceding section, the flow per unit width in the mid length of each plane is equal to qp = 0.00125 m2/s. Since the average velocity in the planes is vp = 0.0417 m/s, the average flow depth in the planes is dp = qp/vp = 0.03 m. Laminar flow in the planes is assumed, with a discharge-depth rating exponent βp = 3. Therefore, the wave celerity in the planes is cp = βpvp = 0.125 m/s.

The flow in the midlength of the channel is equal to 0.5 m3/s. Assume a channel top width Tc = 5 m. Therefore, the flow per unit width in the channel is qc = Qc /Tc = 0.1 m2/s. Since the average velocity in the channel is vc = 0.3 m/s, the average flow depth in the channel (at midlength) is dc = qc / vc = 0.333 m. A wide channel with turbulent Manning friction is assumed, with a discharge-area rating exponent βc = 1.67. Therefore, the wave celerity in the channel is cc = βc vc = 0.5 m/s.

The time of concentration is equal to 1200 s, which is equal to the rainfall duration. The maximum possible peak flow is 1.0 m3/s. The total volume of runoff is 1200 m3.

For simplicity, a space interval of Δx = 100 m is chosen for the planes, and Δy = 200 m for the channel. This amounts to one spatial increment in planes and channel. The time interval is chosen as Δt = 10 minutes. In an actual computer application, a finer grid size would be indicated.

Selection of routing scheme. The chosen routing scheme is the Muskingum-Cunge method, Eq. 9-62, with the routing coefficients calculated by Eqs. 9-74 to 9-76, modified with the addition of lateral inflow. For this purpose, Eq. 9-43 is discretized in the same way as Eq. 9-61, leading to:

Q j+1n+1 = C0 Q j n+1  +  C1 Q j n  +  C2 Q j+1n  +  C3 QL (10-19)

which has the same meaning as Eq. 9-62 with Eqs. 9-74 to 9-76, except for the addition of the lateral inflow term, with routing coefficient [14]:

                  2 C
C3  =  _____________
             1 + C + D
(10-20)

The term QL is the lateral inflow in cubic meters per second. For overland flow routing, the lateral inflow is equal to the effective rainfall (centimeters per hour) times the applicable area (square meters). For channel routing, the lateral inflow is the average distributed lateral inflow (cubic meters per second per meter) multiplied by the channel length (meters).

The grid size and physical parameters allow the calculation of the routing parameters. The Courant numbers in the planes and channel (Eq. 9-69) are: Cp = 0.75 and Cc = 1.5, respectively. The cell Reynolds numbers in the planes and channel (Eq. 9-71) are: Dp = 0.01 and Dc = 0.1, respectively.

With Eqs. 9-74 to 9-76 and 10-20, the routing coefficients in the planes are: C0 = - 0.136, C1 = 0.988, C2 = 0.148, and C3 = 0.852. Likewise, the routing coefficients in the channel are: C0 = 0.231, C1 = 0.923, C2 = - 0.154, and C3= 1.154.

The catchment routing is shown in Table 10-6.

Table 10-6  Two-plane Catchment Routing With Diffusion Wave Technique.
(1) (2) (3) (4) (5) (6) (7)
Time
(min)
Flow (m3/s)
QLp Ip Op QLc Ic Oc
0 0 0 0 0
0.5 0.4260
10 0 0.4260 0 0.4916
0.5 0.9150
20 0 0.4890 0 0.9802
0 0.5614
30 0 0.0724 0 0.4969
0 0.0831
40 0 0.0107 0 0.0194
0 0.0123
50 0 0.0016 0 0.0112
0 0.0018
60 0 0.0002 0 0.0004
0 0.0002
70 0 0 0.0002
Sum 1.00 0.9999 1.9999

Explanation of Table 10-6.

  • Column 1 shows time in minutes.

  • Columns 2-4 show the plane routing and Cols. 5-7 show the channel routing.

  • As in the previous example, the lateral inflow to the plane is 0.5 m3/s (Eq. 10-18).

  • Column 3 shows the upstream inflow to each plane, that is, zero (this example does not consider upstream inflow to the planes).

  • Column 4 is the outflow from each one of the planes, obtained by routing with Eq. 10-19. The first value of Col. 4 (for t = 0) is: 0. The second value of Col. 4 is: (- 0.136 × 0) + (0.988 × 0) + (0.148 × 0) + (0.852 × 0.5) = 0.4260. The third value of Col. 4 is: (- 0.136 × 0) + (0.988 × 0) + (0.148 × 0.4260) + (0.852 × 0.5) = 0.4890. The remaining values in Col. 4 are calculated in a similar way.

  • Column 5 shows the average lateral inflow to the channel, obtained by multipying by 2 the average inflow from each plane (Col. 4) (two planes, each of the same dimensions).

  • Column 6 is the upstream inflow to the channel, that is, zero (this example does not consider upstream inflow to the channel).

  • Column 7 is the outflow from the catchment, obtained by routing with Eq. 10-19. The first value of Col. 6 (for t = 0) is: 0. The second value of Col. 6 is: (0.231 × 0) + (0.923 × 0) + (- 0.154 × 0) + (1.154 × 0.4260) = 0.4916. The third value of Col. 6 is: (0.231 × 0) + (0.923 × 0) + (- 0.154 × 0.4916) + (1.154 × 0.9150) = 0.9802. The remaining values in Col. 6 are calculated in a similar way.

  • As in the previous example, the sums of Cols. 2, 4 and 7 confirm that the runoff volumes are appropriately conserved.

  • Column 7 is the outflow hydrograph from the catchment.

Assessment of Diffusion Wave Method

The calculated outflow hydrograph peak is 0.9802 m3/s, and it occurs at 20 min. This value is very close to the maximum possible peak flow (1 m3/s), revealing that the amount of physical diffusion for this particular example is relatively small.

To study the effect of grid size on the hydrograph calculated by the diffusion wave method, a test similar to that of the previous section is performed with the aid of a computer program. Table 10-7 shows the results obtained by successive grid refinement, using program EH1000C. It is shown that the results are essentially independent of grid size. A coarse grid (one space increment in both planes and channel) results in a peak flow of 0.9802, whereas a fine grid (16 increments in both planes and channel) results in a value of 0.9845. However, it should be noted that the coarse-grid solution (one space increment in planes and channel) exhibits a small but perceptible amount of numerical dispersion, as demonstrated by its peak (0.9802) being somewhat greater than the peak obtained with two increments (0.9716). This is caused by the negative C0 in the planes (C0 = - 0.316).

Table 10-7  Diffusion Wave Catchment Routing:
Effect of Grid Resolution.

Number of
Increments
Δx
(m)
Δy
(m)
Δt
(s)
Peak Flow
(m3/s)
1 100 200 600 0.9802
2 50 100 300 0.9716
4 25 50 150 0.9766
8 12.5 25 75 0.9814
16 6.25 12.5 37.5 0.9845

Unlike in the kinematic wave method, in the diffusion wave technique the numerical diffusion is matched to the physical problem. As in the case of stream channel routing (Section 10.4), this procedure works best when numerical dispersion is minimized, that is, when Courant numbers are kept reasonably close to 1. In practice, substantial deviations from this condition may lead to increased numerical dispersion and associated numerical instability. As with other finite difference schemes, a test of grid independence similar to that shown in Table 10-7 is necessary to verify the convergence of the numerical scheme.


10.6  ASSESSMENT

[Questions]   [Problems]   [References]      [Top]   [Time-area Method]   [Clark Unit Hydrograph]   [Cascade of Linear Reservoirs]   [Kinematic Waves]   [Diffusion Waves]  

Catchment routing techniques have evolved from the simple time-area methods to the more elaborate physically based kinematic and diffusion wave techniques. The variety of existing methods and techniques reflects the fact that no one method is applicable to all cases. Surface runoff in catchments is a complex phenomenon, and research continues to unveil improved ways of solving the problem.

In nature, catchments can be either small, midsize, or large. Runoff processes are nonlinear and distributed; nonlinear in the sense that the parameters do not increase in the same proportion as the flow, and distributed in the sense that the parameters vary within the catchment. Hydrographs originating in catchment runoff are generally concentrated, diffused, and dispersed. From the mathematical standpoint, concentration is a first-order process, diffusion a second-order process, and dispersion a third-order process [5]. Being a first-order process, runoff concentration is the primary mechanism; it is also referred to as translation or wave travel.

For certain applications, diffusion also plays an important role-for instance, in the modeling of runoff response in catchments with mild slopes. Early description of the diffusion mechanism referred to it as storage. The usage is so widespread that it has been preserved in this book. In a routing context, diffusion and (longitudinal channel) storage have essentially the same meaning. Dispersion is a third order process; therefore, it is usually much smaller than runoff concentration and diffusion. Physical dispersion translates into wave steepening, a concept similar to that of skewness, or the third moment of a statistical distribution (Chapter 6).

The time-area method can calculate concentration but it cannot account for diffusion or dispersion. Therefore, the time-area method should be limited to small and midsize catchments where translation is by far the predominant mechanism. When used indiscriminately, the time-area method always overestimates the peak of the outflow hydrograph. If necessary, diffusion can be added by routing the time-area hydrograph through a linear reservoir. However, the storage constant would have to be determined either from measured data or by synthetic means. While the time-area method is a lumped method, it can be used as a component of larger network models which have a distributed structure [19].

In principle, the method of cascade of linear reservoirs accounts for runoff diffusion only. However, the connection of several linear reservoirs in series provides enough diffusion so that translation is actually being simulated by means of diffusion. The method is linear insofar as the routing parameters are determined by calibration. Nevertheless, by successive calibrations at different flow levels, from low to high, the method can be made to reflect the nonlinearity actually existing in nature. The method is lumped, but it can be used as a component of larger network models that have a distributed structure. The method has been successfully applied to very large basins (in excess of 10,000 km2), which exhibit substantial amounts of diffusion. For such large basins, the distributed models of the kinematic and diffusion type may be impractical due to the prohibitive amount of data required to properly specify the spatial diversity.

The kinematic wave model provides translation and diffusion, the latter, however, due only to the finite grid size. The method can be linear or nonlinear, and lumped or distributed, depending on the numerical scheme and input data. The method is applicable to small catchments with steep slopes where diffusion is small and can be controlled by grid refinement. Theoretically, the method could also be applicable to midsize catchments, as long as physical diffusion remains small. In practice, the larger the catchment, the more unlikely it is that physical diffusion is negligible.

The distributed nature of kinematic wave models results in substantial data needs; the use of average parameters would render the model lumped, with the consequent loss of detail. Another important consideration in kinematic wave models is the validity of the geometric configuration. For instance, two-plane descriptions are adequate as long as the catchment geometry fits the two-plane model configuration. Otherwise, a certain amount of lumping would be introduced by the model's inability to properly account for the physical detail. In practice, the larger the catchment, the more difficult it is to fit catchments within two-plane descriptions. Multiple-plane descriptions are possible but invariably lead to additional complexity [8].

The diffusion wave technique (Muskingum-Cunge scheme) provides translation and diffusion and, unlike the kinematic wave model, its solution is generally independent of grid size. Therefore, it is applicable to catchments with substantial amounts of physical diffusion, either small catchments of mild slope or midsize catchments with average slopes. The method is linear or nonlinear, and lumped or distributed, depending on the numerical scheme and input data. As with kinematic models, the validity of the geometric configuration is important in diffusion wave catchment models. Unless the model's geometric abstraction is a reasonable representation of the catchment's actual geometry, the degree of lumping introduced may tend to mask the distributed nature of the model.

The preceding comments have referred specifically to catchment routing, i.e., the conversion of effective rainfall into runoff. In practice, comparative evaluations of the performance of catchment routing methods/models are hampered by the fact that it is seldom possible to determine effective rainfall with any degree of certainty. Hydrologic abstractions are time-variant, distributed, and nonlinear. Therefore, a proper estimation of hydrologic abstractions is crucial to the performance evaluation of catchment routing models.

From this discussion, it can be concluded that no one method or model is suitable for all applications. All have strengths and weaknesses; they are either simple or complex, and suffer from lack of detail or require a substantial amount of data for their successful operation. In practice, the choice of catchment routing method remains one of individual preference and experience.


QUESTIONS

[Problems]   [References]      [Top]   [Time-area Method]   [Clark Unit Hydrograph]   [Cascade of Linear Reservoirs]   [Kinematic Waves]   [Diffusion Waves]   [Assessment]  

  1. What is the difference between hydrologic and hydraulic methods of catchment routing?

  2. What are catchment isochrones? How are they determined?

  3. How is time of concentration defined when using hydrographs generated by the time-area method?

  4. When can the Clark unit hydrograph be considered synthetic? Explain.

  5. What is the difference between translation and diffusion? How are the time-area and rational methods related?

  6. How can the linear reservoir storage coefficient be determined for runoff data?

  7. What is the principle behind the method of cascade of linear reservoirs used in catchment routing?

  8. Why is it necessary to use two schemes in catchment routing using first-order kinematic wave techniques?

  9. Why is a kinematic wave solution using numerical techniques usually grid-dependent? Why is the diffusion wave solution grid-independent?


PROBLEMS

[References]      [Top]   [Time-area Method]   [Clark Unit Hydrograph]   [Cascade of Linear Reservoirs]   [Kinematic Waves]   [Diffusion Waves]   [Assessment]   [Questions]  

  1. A 45-km2 catchment has a 6-h time fo concentration with isochrones at 2-h intervals, resulting in the following time-area histogram:

    Time (h) 0-2 2-4 4-6
    Area (km2) 9 21 15

    Use the time-area method to calculate the outflow hydrograph from the following effective storm pattern:

    Time (h) 0-2 2-4 4-6 6-8 8-10 10-12
    Effective rainfall (cm/h) 0.5 1.0 2.0 3.0 1.0 0.5

    Use a spreadsheet. Verify your results with ONLINE ROUTING06.

  2. A 68-km2 catchment has a 4-h time of concentration with isochrones at 1-h intervals resulting in the following time-area histogram:

    Time (h) 0-1 1-2 2-3 3-4
    Area (km2) 12 19 26 11

    Use the time-area method to calculate the outflow hydrograph from the following storm:

    Time (h) 0-1 1-2 2-3 3-4 4-5 5-6
    Total rainfall (cm/h) 1 2 4 3 2 1

    The runoff curve number is CN = 75. Use a spreadsheet. Verify your results with ONLINE ROUTING06.

  3. An 82-km2 catchment has the following characteristics:

    Time (h) 0-1 1-2 2-3 3-4 4-5 5-6
    Area (km2) 14 22 29 17 0 0
    Effective rainfall (cm/h) 0.8 1.2 1.3 2.1 0.7 0.5

    Using a spreadsheet, calculate the outflow hydrograph by the time-area method. Verify your results with ONLINE ROUTING06.

  4. Use the Clark method with a spreadsheet to derive a 3-h unit hydrograph for a catchment with the following time-area diagram:

    Time (h) 0-3 3-6 6-9 9-12
    Area (km2) 57 72 39 15

    Use Δt = 3 h and K = 3 h. Verify your results with ONLINE ROUTING07.

  5. Use the Clark method with a spreadsheet to derive a 1-h unit hydrograph for a catchment with the following time-area diagram:

    Time (h) 0-1 1-2 2-3 3-4 4-5 5-6
    Area (km2) 12 20 42 66 30 16

    Use Δt = 1 h and K = 2 h. Verify your results with ONLINE ROUTING07.

  6. Use the Clark method with a spreadsheet to derive a 2-h unit hydrograph for a catchment with the following time-area diagram:

    Time (h) 0-1 1-2 2-3 3-4 4-5 5-6
    Area (km2) 10 20 30 20 12 8

    Use Δt = 1 h and K = 4 h. Verify your results with ONLINE ROUTING07.

  7. The 2-h unit hydrograph for a 92-km2 catchment is the following:

    Time (h) 0 2 4 6 8 10 12 14
    Flow (m3/s) 0 2.778 8.611 14.333 19.433 21.938 20.942 15.897
    Time (h) 16 18 20 22 24 26 28 30
    Flow (m3/s) 9.538 5.722 3.433 2.061 1.236 0.742 0.444 0.267

    Given the following time-area diagram, what is the linear reservoir storage constant in the Clark method?

    Time (h) 0-2 2-4 4-6 6-8 8-10 10-12
    Area (km2) 10 15 18 21 16 12

  8. A 1-h unit hydrograph derived from measured data has the following ordinates:

    Time (h) 0 1 2 3 4 5 6 7
    Flow (m3/s) 0 7 22 48 60 90 74 47
    Time (h) 8 9 10 11 12 13 14
    Flow (m3/s) 28 17 10 6 4 3 2

    Assuming a time of concentration tc = 6 h, calculate the linear reservoir storage constant in the Clark method.

  9. Using ONLINE ROUTING08, route the effective rainfall hyetograph of Example 10-3 (with time interval Δt = 6 h) using: (a) K = 12 h and N = 4; and (b) K = 18 h and N = 3.

  10. Use ONLINE ROUTING08 to route the following storm hyetograph through a 535-km2 catchment.

    Time (h) 0-6 6-12 12-18 18-24
    Effective rainfall depth (cm) 1.0 1.5 2.5 1.2
    Effective rainfall depth (cm/hr) 0.1667 0.2500 0.4167 0.2000

    Set K = 12 h and N = 3. Report peak outflow and time-to-peak.

  11. Using the method of cascade of linear reservoirs, derive a 3-h unit hydrograph for a 432-km2 basin. Assume Δt = 3 h, K = 9 h and two reservoirs (N = 2). Verify the results with ONLINE ROUTING08.

  12. Using the method of cascade of linear reservoirs, derive a 6-h unit hydrograph for a 1235-km2 basin. Assume Δt = 6 h, K = 6 h and N = 4. Verify the results with ONLINE ROUTING08.

  13. Derive Eq. 10-16.

  14. Derive Eq. 10-17.

  15. Derive Eq. 10-20.


REFERENCES

   [Top]   [Time-area Method]   [Clark Unit Hydrograph]   [Cascade of Linear Reservoirs]   [Kinematic Waves]   [Diffusion Waves]   [Assessment]   [Questions]   [Problems]  

  1. Agricultural Research Service, U.S. Department of Agriculture. (1973). "Linear Theory of Hydrologic Systems," Technical Bulletin No. 1468, J. C. I. Dooge, author), Washington, D.C.

  2. Clark, C. O. (1945). "Storage and the Unit Hydrograph," Transactions. ASCE, Vol. 110. pp. 1416-1446.

  3. Cordery, 1. (1968). "Synthetic Unit Graphs for Small Catchments in Eastern New South Wales," Civil Engineering Transactions of the Institution of Engineers, Australia, Vol. 10, pp. 47-58.

  4. Eaton, T. D. (1954) . "The Derivation and Synthesis of Unit Hydrographs When Rainfall Records Are Inadequate," Institution of Engineers, Australia, Vol. 26, pp. 239-246.

  5. Ferrick, M. G., J. Bilmes, and S. E. Long. (1983). "Analysis of a Diffusion Wave Flow Routing Model with Application to Flow in Tailwaters," U.S. Army Corps of Engineers, Cold Regions Research and Engineering Laboratory, CRREL Report 83-7, March.

  6. Hromadka, T. V. , R. H. McCuen, and C. C. Yen. (1987). "Comparison of Overland Flow Hydrograph Models," Journal of Hydraulic Engineering, ASCE. Vol. 113, No. 11, November. pp. 1422-1440.

  7. Hydrologic Engineering Center, U.S. Army Corps of Engineers. (1990). "HEC-l, Flood Hydrograph Package, User's Manual," Version 4.0, September.

  8. Kibler, D. F., and D. A. Woolhiser. (1970). "The Kinematic Cascade as a Hydrologic Model," Colorado State University Hydrology Paper No. 39, Fort Collins, Colorado.

  9. Lighthill, M. J., and G. B. Whitham. (1955). "On Kinematic Waves. 1. Flood Movement in Long Rivers," Proceedings. Royal Society, London, Vol. A229. May, pp. 281-316.

  10. Nash, J. E. (1958). "The Form of the Instantaneous Unit Hydrograph," International Association for Scientific Hydrology, Publication No. 42, No.3, pp. 114-118.

  11. O'Kelly, 1. J. (1955). "The Employment of Unit Hydrograph to Determine the Flows of Irish Arterial Drainage Channels," Proceedings of the Institution of Civil Engineers, Dublin, Ireland, Vol. 4, No.3, pp. 365-412.

  12. Ponce, V. M., R. M. Li, and D. B. Simons. (1978). "Applicability of Kinematic and Diffusion Models," Journal of the Hydraulics Division, ASCE, Vol. 104, No. HY3, March, pp. 353-360.

  13. Ponce, V. M., and D. Windingland. (1985). "Kinematic Shock: Sensitivity Analysis," Journal of Hydraulic Engineering, ASCE. Vol. 111, No. 4, April, pp. 600-611.

  14. Ponce, V. M. (1986). "Diffusion Wave Modeling of Catchment Dynamics," Journal of Hydraulic Engineering, ASCE, Vol. 112, No. 8, August, pp. 716-727.

  15. Ponce, V. M., and P. T. Huston. 1994. "New perspective on the convection-diffusion-dispersion equation," Water Resources Research, Vol. 30, No. 5, May, pp. 1619-1620.

  16. Orlandini, S., and R. Rosso. (1996). "Diffusion Wave Modeling of Distributed Catchment Dynamics," Journal of Hydrologic Engineering, ASCE, Vol. 1, No. 3, July, pp. 103-113.

  17. Rao, R. A., J. W. Delleur, and B. S. Sarma. (1972). "Conceptual Hydrologic Models for Urbanizing Basins," Journal of the Hydraulics Division, ASCE, Vol. 98, No. HY7, July, pp. 1205-1220.

  18. Singh, V. P. (1988). Hydrologic Systems. Rainfall-Runoff Modeling, Vol. 1. Englewood Cliffs, N.J.: Prentice-Hall.

  19. Stall, J. B., and Terstriep, M. L. (1972). "Storm-Sewer Design: An Evaluation of the RRL Method," EPA Technology Series EPA-R2-72-068, October.

  20. U.S. Army Engineer Division, North Pacific Division. (1986). "Program Description and User Manual for SSARR Model, Streamflow Synthesis and Reservoir Regulation," Portland, Oregon, April.

  21. Wooding, R. A. (1965). "A Hydraulic Model for the Catchment-stream Problem, 1. Kinematic Wave Theory," Journal of Hydrology, Vol. 3, pp. 254-267.

  22. Wooding, R. A. (1965). "A Hydraulic Model for the Catchment-stream Problem, II. Numerical Solutions," Journal of Hydrology, Vol. 3, 1965, pp. 268-282.

  23. Wooding, R. A. (1966). "A Hydraulic Model for the Catchment-stream Problem, III. Comparison with Runoff Observations," Journal of Hydrology, Vol. 4, pp. 21-37.

  24. Woolhiser, D. A., and Liggett, J. A. (1967). "Unsteady One-dimensional Flow Over a Plane: The Rising Hydrograph," Water Resources Research, Vol. 3, No.3, pp. 753-771.

SUGGESTED READINGS

  • Agricultural Research Service, U.S. Department of Agriculture. (1973). "Linear Theory of Hydrologic Systems," Technical Bulletin No. 1468, (J. C. I. Dooge, author), Washington, D.C.

  • Clark, C. O. (1945). "Storage and the Unit Hydrograph," Transactions, ASCE, Vol. 110, 1945, pp. 1416-1446.

  • Hydrologic Engineering Center, U.S. Army Corps of Engineers. (1990). "HEC-1, Flood Hydrograph Package, User's Manual," Version 4.0, September.

  • U.S. Army Engineer Division, North Pacific Division. (1986). "Program Description and User Manual for SSARR Model, Streamflow Synthesis and Reservoir Regulation," Portland, Oregon, April.

  • Wooding, R. A. (1965). "A Hydraulic Model for the Catchment-stream Problem, I. Kinematic Wave Theory," Journal of Hydrology, Vol. 3, 1965, pp. 254-267.


http://engineeringhydrology.sdsu.edu
150714