Abstract
This study theoretically evaluated the basin concentration time and storage coefficient with their empirical formulas available worldwide. The evaluation results were also validated in the application to major dam basins in Korea. The findings are summarized as follows. As a result of analytical analysis, the concentration time was found to be proportional to the main channel length under laminar flow conditions and to the square of it under turbulent flow conditions, but inversely proportional to the channel slope. It was also found that the storage coefficient and the concentration time are linearly but loosely related. Most empirical formulas for the concentration time concurred with the basic equation form, but just a few for the storage coefficient. Applications to major dam basins in Korea also showed that the concentration time agrees well with the result of theoretical analysis. However, the behavior of the storage coefficient varied much, basin by basin, indicating that additional factors may be needed to explain it.
INTRODUCTION
The runoff characteristics of a basin can be quantified by interpreting the components of a flood hydrograph. Tc and K are critical components when determining the peak flow and peak time of a flood hydrograph. The Clark instantaneous unit hydrograph (IUH), which is generally used for basin flood routing in Korea, can also be derived by these two parameters, Tc and K (Clark 1945). That is, under the assumption of linear system theory, these two parameters can sufficiently represent the characteristics of the rainfall–runoff process in a basin.
However, even though only two parameters are involved in the construction of the Clark IUH, it is not simple to estimate them. Actually, it may be impossible to determine the unique set of parameters that can be applied to various rainfall–runoff events. This is simply because the nonlinear basin system is assumed to be linear. It is well known that the rainfall–runoff process in a basin is nonlinear (Kundzewicz & Napiórkowski 1986; Sinha et al. 2015). Furthermore, as these two parameters are correlated with each other, their estimation procedure can be very complex. In practice, even for an observed rainfall–runoff event, it is practically impossible to estimate a unique set of Tc and K.
On the other hand, in ungauged basins, empirical formulas are used to determine Tc and K. However, as empirical formulas are strongly dependent on the basin or channel characteristics, and sometimes provide very different estimates, the background, such as basin area, basin slope, land use, land cover, etc., should be carefully checked when applying these formulas (Jeong & Yoon 2007; Yoo 2009).
The basin or channel characteristics involved in the empirical formulas for Tc and K, such as basin area, shape factor, channel length, channel slope, rainfall intensity, etc., are diverse. Nevertheless, it has been argued that the components of those formulas are very similar each other (Yoo 2009). This finding indicates that there exist major contributing factors to determine Tc and K, as well as their basic equation forms. If this assumption is true, the regionalization for the consistent application of an empirical formula to a large basin will be practical. However, major contributing factors of empirical formulas for Tc and K have not been investigated in depth by previous researchers.
The main objective of this study is to evaluate the fundamental factors of Tc and K with their theoretical background and empirical formulas available worldwide. Ultimately this study is going to show if a valid basic equation form exists for Tc and K. The result will then be applied to major dam basins in Korea to evaluate its applicability.
MATERIALS AND METHODS
Theory for concentration time



Furthermore, the definition of concentration time, which is the travel time of a rain drop to travel from the furthest location to the outlet of a basin, supports the mathematical expression as a function of L/S or L2/S. Simply the travel time is closely related to the channel length and slope. For these reasons, many empirical formulas for concentration time are expressed as a function of L/S or L2/S (Kirpich 1940; Izzard & Hicks 1947; Johnstone & Cross 1949; Kerby 1959; Carter 1961; Morgali & Linsley 1965; Espey & Winslow 1974; SCS 1975).
Relationship between concentration time and storage coefficient


Study basins
A total of five dam basins in Chungju, Namgang, Andong, Imha and Hapcheon were considered. Among these five dam basins, only Chungju dam basin belongs to the Han river basin covering the central part of the Korean Peninsula. The remaining dam basins in Namgang, Andong, Imha and Hapcheon are located within the Nakdong river basin covering the south-eastern part of the Korean Peninsula. Figure 1 shows the location of the selected basins within the Korean Peninsula.
Each basin contains many stream gauge stations, among which this study selected only those less affected by dam or other hydraulic structures to secure the accurate rainfall–runoff data for reasonable parameter estimation. As a result, 18 stream gauge stations within Chungju, nine within Namgang, six within Andong, four within Imha and four within Hapcheon were selected. Table 1 summarizes the topographic characteristics of the sub-basins considered in this study. In Table 1, A is the basin area (km2), L is the channel length (km), S is the basin slope.
Topographic characteristics of sub-basins considered in this study
Basin . | Stream gauge station . | A (km²) . | L (km) . | S . | Basin . | Stream gauge station . | A (km²) . | L (km) . | S . |
---|---|---|---|---|---|---|---|---|---|
Chungju dam | Imokjeong | 55.8 | 16.6 | 0.016 | Andong dam | Jangseong | 123.4 | 21.6 | 0.013 |
Jangpyeong-gyo | 105.1 | 26.0 | 0.015 | Dosan | 227.6 | 36.0 | 0.010 | ||
Baekokpo | 143.9 | 23.0 | 0.014 | Buncheon | 502.1 | 59.0 | 0.0078 | ||
Anheung | 187.0 | 31.8 | 0.0049 | Socheon | 642.8 | 68.4 | 0.0069 | ||
Songcheon | 349.5 | 62.8 | 0.0089 | Ungok | 1,144.2 | 107.8 | 0.0052 | ||
Sanganmi | 393.3 | 44.2 | 0.0085 | Andong dam | 1,584.0 | 139.1 | 0.0039 | ||
Najeon | 452.6 | 59.1 | 0.0081 | Namgang dam | Hamyang | 123.6 | 22.2 | 0.023 | |
Bangrim | 527.2 | 51.9 | 0.0077 | Taesu | 243.2 | 28.3 | 0.037 | ||
Jucheon | 533.2 | 71.0 | 0.0042 | Samga | 101.0 | 13.3 | 0.0068 | ||
Shincheon | 598.3 | 84.1 | 0.0040 | Aneui | 215.5 | 31.7 | 0.015 | ||
Pyeongchang | 695.7 | 74.8 | 0.0061 | Changchon | 328.3 | 40.0 | 0.018 | ||
Panun | 879.1 | 90.3 | 0.0049 | Macheon | 315.4 | 31.3 | 0.0093 | ||
Yeongwol 1 | 1,524.1 | 125.7 | 0.0038 | Imcheon | 459.0 | 47.2 | 0.0093 | ||
Jeongseon 2 | 1,688.1 | 103.5 | 0.0050 | Sancheong | 1,122.3 | 67.2 | 0.0067 | ||
Yeongwol 2 | 2,283.4 | 181.5 | 0.0035 | Namgang dam | 2,285.0 | 111.2 | 0.0041 | ||
Yeongwol | 2,440.4 | 186.5 | 0.0033 | Imha dam | Yoengyang | 334.0 | 36.6 | 0.0090 | |
Youngchun | 4,775.0 | 189.1 | 0.0031 | Cheongsong | 305.0 | 41.0 | 0.0050 | ||
Chungju dam | 6,648.0 | 282.2 | 0.0017 | Giran | 411.7 | 54.7 | 0.0049 | ||
Imha dam | 1,361.0 | 90.7 | 0.0035 | ||||||
Hapcheon dam | Jisan | 159.9 | 26.0 | 0.019 | |||||
Geochang 2 | 180.7 | 29.5 | 0.015 | ||||||
Geochang 1 | 227.0 | 30.7 | 0.014 | ||||||
Hapcheon dam | 925.0 | 64.2 | 0.0054 |
Basin . | Stream gauge station . | A (km²) . | L (km) . | S . | Basin . | Stream gauge station . | A (km²) . | L (km) . | S . |
---|---|---|---|---|---|---|---|---|---|
Chungju dam | Imokjeong | 55.8 | 16.6 | 0.016 | Andong dam | Jangseong | 123.4 | 21.6 | 0.013 |
Jangpyeong-gyo | 105.1 | 26.0 | 0.015 | Dosan | 227.6 | 36.0 | 0.010 | ||
Baekokpo | 143.9 | 23.0 | 0.014 | Buncheon | 502.1 | 59.0 | 0.0078 | ||
Anheung | 187.0 | 31.8 | 0.0049 | Socheon | 642.8 | 68.4 | 0.0069 | ||
Songcheon | 349.5 | 62.8 | 0.0089 | Ungok | 1,144.2 | 107.8 | 0.0052 | ||
Sanganmi | 393.3 | 44.2 | 0.0085 | Andong dam | 1,584.0 | 139.1 | 0.0039 | ||
Najeon | 452.6 | 59.1 | 0.0081 | Namgang dam | Hamyang | 123.6 | 22.2 | 0.023 | |
Bangrim | 527.2 | 51.9 | 0.0077 | Taesu | 243.2 | 28.3 | 0.037 | ||
Jucheon | 533.2 | 71.0 | 0.0042 | Samga | 101.0 | 13.3 | 0.0068 | ||
Shincheon | 598.3 | 84.1 | 0.0040 | Aneui | 215.5 | 31.7 | 0.015 | ||
Pyeongchang | 695.7 | 74.8 | 0.0061 | Changchon | 328.3 | 40.0 | 0.018 | ||
Panun | 879.1 | 90.3 | 0.0049 | Macheon | 315.4 | 31.3 | 0.0093 | ||
Yeongwol 1 | 1,524.1 | 125.7 | 0.0038 | Imcheon | 459.0 | 47.2 | 0.0093 | ||
Jeongseon 2 | 1,688.1 | 103.5 | 0.0050 | Sancheong | 1,122.3 | 67.2 | 0.0067 | ||
Yeongwol 2 | 2,283.4 | 181.5 | 0.0035 | Namgang dam | 2,285.0 | 111.2 | 0.0041 | ||
Yeongwol | 2,440.4 | 186.5 | 0.0033 | Imha dam | Yoengyang | 334.0 | 36.6 | 0.0090 | |
Youngchun | 4,775.0 | 189.1 | 0.0031 | Cheongsong | 305.0 | 41.0 | 0.0050 | ||
Chungju dam | 6,648.0 | 282.2 | 0.0017 | Giran | 411.7 | 54.7 | 0.0049 | ||
Imha dam | 1,361.0 | 90.7 | 0.0035 | ||||||
Hapcheon dam | Jisan | 159.9 | 26.0 | 0.019 | |||||
Geochang 2 | 180.7 | 29.5 | 0.015 | ||||||
Geochang 1 | 227.0 | 30.7 | 0.014 | ||||||
Hapcheon dam | 925.0 | 64.2 | 0.0054 |
Preparation of storm event data
Hourly rainfall data were used in this study. The observation period of the data varies from 1 year to 36 years depending on the rain gauge station. In cases where several rain gauges were available, the basin average rainfall data were prepared by the Thiessen polygon method. Major storm events were then separated by applying some conditions like the mean rainfall intensity 10 mm/hr or maximum rainfall intensity 50 mm/hr. Additionally, only those storm events satisfying the AMC (Antecedent Moisture Condition) III condition were selected to make sure of enough runoff volume. As a result, each stream gauge station could secure storm events from a minimum 2 to a maximum 47, and they were used for the estimation of Tc and K.
EVALUATION OF EMPIRICAL FORMULAS
Concentration time
Various empirical formulas for Tc have been proposed globally and are summarized in Table 2. In the table, Tc is the concentration time (hours), A is the basin area (km2), L is the channel length (km), S is the basin slope or channel slope, N is the retardance coefficient, V is the mean velocity (m/sec), i is the rainfall intensity (in/hr), and n is the Manning's roughness coefficient.
Empirical formulas for Tc
Name . | Formula . | Name . | Formula . |
---|---|---|---|
Kraven (I) (JSCE 1999) | ![]() | Espey & Winslow (1974) | ![]() |
Kraven (II) (JSCE 1999) | ![]() | SCS (1975) | ![]() |
Rziha (1876) | ![]() | SCS lag (1985) | ![]() |
Kirpich (1940) | ![]() | Ahn & Lee (1986) | ![]() |
Johnstone & Cross (1949) | ![]() | FAA (1970) | ![]() |
Kerby (1959) | ![]() | USGS (2000) | ![]() |
Kerby-Hathaway (Hathaway 1945) | ![]() | Yoon & Park (2002) | ![]() |
Carter (1961) | ![]() | Jung (2005) | ![]() |
Morgali & Linsley (1965) | ![]() | California DoT (1955) | ![]() |
Izzard & Hicks (1947) | ![]() | Kim (2015) | ![]() |
MOCT (1974) | ![]() |
Name . | Formula . | Name . | Formula . |
---|---|---|---|
Kraven (I) (JSCE 1999) | ![]() | Espey & Winslow (1974) | ![]() |
Kraven (II) (JSCE 1999) | ![]() | SCS (1975) | ![]() |
Rziha (1876) | ![]() | SCS lag (1985) | ![]() |
Kirpich (1940) | ![]() | Ahn & Lee (1986) | ![]() |
Johnstone & Cross (1949) | ![]() | FAA (1970) | ![]() |
Kerby (1959) | ![]() | USGS (2000) | ![]() |
Kerby-Hathaway (Hathaway 1945) | ![]() | Yoon & Park (2002) | ![]() |
Carter (1961) | ![]() | Jung (2005) | ![]() |
Morgali & Linsley (1965) | ![]() | California DoT (1955) | ![]() |
Izzard & Hicks (1947) | ![]() | Kim (2015) | ![]() |
MOCT (1974) | ![]() |
The empirical formulas for Tc summarized in Table 2 are compared with those in the theoretical analysis. Most of the empirical formulas are expressed as a function of both L and S. As, theoretically, Tc should be represented by a function of L2/S for turbulent flow and by a function of L/S for laminar flow, those empirical formulas in Table 2 can be compared effectively on the plane of exponents of L and S, as shown in Figure 2(a). Each empirical formula for Tc and K is expressed as a single value in Figure 2. Since the exponents of L and S in the empirical formulas considered in this study do not exceed the value of 2, the ranges of the x-axis and y-axis are set from 0 to 2.
Comparison of empirical formulas on the plane of exponents of S and L.
The dotted lines in Figure 2 are the reference lines for the functions with slopes 2 and 1, respectively, representing the ratios between exponents of L and S for L2/S and L/S. These reference lines are based on the theoretical review that Tc and K should follow a function of L2/S for the turbulent flow and a function of L/S for the laminar flow. That is, the nearer to the reference line the point (indicating an empirical formula) is located, the better it is matched to its own theoretical background.
As can be seen in Figure 2(a), formulas like Kirpich (1940), Kraven (I) (JSCE 1999), Kraven (II) (JSCE 1999), California DoT (1955), Kerby (1959), SCS (1975), Morgali & Linsley (1965), Carter (1961), MOCT (1974), Espey & Winslow (1974), and Kerby-Hathaway (Hathaway 1945) are found to be matched well to the theoretical result for turbulent flow. That is, the ratio of exponents of L and S is 2. Clearly, these formulas were developed for the condition of heavy rainfall in a steep upstream basin. On the other hand, the empirical formulas like Johnstone & Cross (1949) and Izzard & Hicks (1947) are found to follow the case of laminar flow; that is, the ratio of exponents of L and S is 1. Some other formulas like Rziha (1876), FAA (1970) and SCS lag (1985) are found to be in between the turbulent and laminar flow conditions.
Storage coefficient
The empirical formulas for K identified in this study are summarized in Table 3, where K represents the storage coefficient (hour), C is the runoff coefficient, L is the channel length (km), A is the basin area (km2), S is the channel slope, and b and a are the correction factors.
Empirical formulas for K
Name . | Formula . | Name . | Formula . |
---|---|---|---|
Clark (1945) | ![]() | USGS (2000) | ![]() |
Linsley (1945) | ![]() | Yoon & Park (2002) | ![]() |
Laurenson (1962) | ![]() | Jung (2005) | ![]() |
Russell (Russell et al. 1979) | ![]() | Go (2014) | ![]() |
Sabol (1988) | ![]() | Kim (2015) | ![]() |
Yoon et al. (1994) | ![]() |
Name . | Formula . | Name . | Formula . |
---|---|---|---|
Clark (1945) | ![]() | USGS (2000) | ![]() |
Linsley (1945) | ![]() | Yoon & Park (2002) | ![]() |
Laurenson (1962) | ![]() | Jung (2005) | ![]() |
Russell (Russell et al. 1979) | ![]() | Go (2014) | ![]() |
Sabol (1988) | ![]() | Kim (2015) | ![]() |
Yoon et al. (1994) | ![]() |
Based on the theoretical analysis in the previous section, it was found that K can be proportional to Tc. Thus, the forms of empirical formulas can be analyzed similarly to those for Tc. Like Figure 2(a), Figure 2(b) compares empirical formulas of K in Table 3 on the plane of exponents of L and S. Also, each point in Figure 2(b) indicates an empirical formula located by the exponents of S and L.
In Figure 2(b), formulas like Clark (1945), Russell (Russell et al. 1979), Sabol (1988) and Jung (2005) closely match the theoretical result. The empirical formulas like Clark (1945) and Jung (2005) show a ratio of exponents of L and S of 2. In the Russell (Russell et al. 1979) and Sabol (1988) formulas, the relationship between Tc and K is only given by a proportionality constant. The formula by Yoon & Park (2002) shows a slightly higher ratio of exponents (approximately 2.5), and the Kim (2015) formula is between turbulent and laminar flow (approximately 1.5). However, the ratio of exponents is 4 for the Linsley (1945) formula and is even higher for the Yoon et al. (1994) and Go (2014) formulas.
Relationship between concentration time and storage coefficient
When comparing the empirical formulas for Tc and K, it is obvious that the empirical formulas for Tc more closely concur with the theoretical result. More variation was found in the empirical formulas for K. For example, A is considered important in many empirical equations including Linsley (1945), Laurenson (1962), Kim (2015), etc. The exponents used for A in the Linsley (1945), Laurenson (1962), and Kim (2015) formulas are rather high at 0.5, 0.27, and 0.238, respectively. Additional consideration of L to A also indicates that the basin shape factor plays an important role, which can be found in the formulas by Linsley (1945), Yoon et al. (1994), Go (2014), and Kim (2015).
Also, based on Figure 2(b), it can be easily conjectured that K in a basin is not that simple like Tc. A linear relationship between the two may not be satisfied in most river basins. Also, K in a basin may not be simplified as a function of channel characteristics. Characteristics of basin shape and channel network, as adopted in many empirical formulas, may improve the accuracy of their application. At this moment, it may not be concluded that K in a basin is simply proportional to Tc.
EVALUATION FOR MAJOR DAM BASINS IN KOREA
Estimation of concentration time and storage coefficient
This study estimated Tc and K of an observed storm event using the method proposed by Yoo et al. (2013). This method is a recursive approach considering the structure of the Nash IUH (Nash 1957). Yoo et al. (2013) showed that the Nash IUH has a distinct relationship between Tc and K, as in Equation (7). Also, as shown in Equations (5) and (6), Tc and K are correlated nonlinearly. In fact, Tc shows more sensitivity to x than K does. This difference causes the strong nonlinear behavior of these two parameters when they are estimated from rainfall–runoff measurements. K can show rather stable behavior without significant variation, but Tc cannot. However, the method by Yoo et al. (2013) was found to overcome this problem when estimating Tc and K using the observed rainfall–runoff data.
Results and discussion
Empirical equations for Tc and K were derived for each basin using those estimated using the observed data. Figure 3 compares the resulting empirical equations on log–log paper.
Comparison of empirical equations of Tc and K derived for the dam basins in this study.
Comparison of empirical equations of Tc and K derived for the dam basins in this study.
As can be seen in Figure 3, both derived empirical equations for Tc and K are well matched to the observed values. Especially, Tc seems very well explained by L2/S in all basins. All the empirical equations are found to have similar slopes on log–log paper. On the other hand, the empirical equations for K show different behavior; especially the slopes and intercepts are found to be different basin by basin. That is, K seems to vary much basin by basin. In fact, this phenomenon was also found in the previous section where various forms of empirical formulas were analyzed. Differently from Tc, more various forms of empirical formulas exist, and they also consider other basin characteristics like A and shape factor. Simply L and S may not be enough to explain K.
Finally, this study evaluated the relationship between Tc and K, which was found to be linear. The coefficients of determination of those regression lines are from 0.932 to 0.995. However, it also true that each line is different from the others. The difference may be quantified by the Russell parameter (Russell et al. 1979), which is the relation between K and Tc or simply the slope of the line in this figure. The highest Russell parameter was found in Hapcheon dam basin, which is 1.262, and in Chungju dam basin it was the smallest at 0.986. Even though the Russell parameter in Hapcheon dam basin is a bit higher than 1.2, all basins are included in the range 0.8 ∼ 1.2 of natural basins in Korea (Jeong & Yoon 2007).
CONCLUSIONS
This study tried the theoretical evaluation of Tc and K with more than 20 empirical formulas available worldwide. The evaluation result was also confirmed in the application to major dam basins in Korea. The findings of this study are as follows:
- (1)
Tc could be expressed as a function of channel characteristics: L and S. Especially, Tc was found to be expressed as a function of L/S under laminar flow conditions and L2/S under turbulent flow conditions.
- (2)
Most empirical formulas for Tc are found to follow the basic equation form considered in this study, however, those for K are found to vary considerably. It seems that K may not be explained fully by L2/S. Additional factors may be needed to effectively explain the storage effect of a basin.
- (3)
In the application to several dam basins in Korea, it was found that both Tc and K could be well modeled by the factor L2/S. Especially, the relationship between Tc and L2/S was very consistent in all basins, but the behavior of K varied much basin by basin.
ACKNOWLEDGEMENTS
This work is supported by the Korea Agency for Infrastructure Technology Advancement (KAIA) grant funded by the Ministry of Land, Infrastructure and Transport (Grant 18AWMP-B083066-05).