Next Article in Journal
Evapotranspiration and Quantitative Partitioning of Spring Maize with Drip Irrigation under Mulch in an Arid Region of Northwestern China
Next Article in Special Issue
Selected Issues of Adaptive Water Management on the Example of the Białka River Basin
Previous Article in Journal
Impact of Sediment Layer on Longitudinal Dispersion in Sewer Systems
Previous Article in Special Issue
Hydrological Response of the Kunhar River Basin in Pakistan to Climate Change and Anthropogenic Impacts on Runoff Characteristics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a New 8-Parameter Muskingum Flood Routing Model with Modified Inflows

School of Civil Engineering, Chungbuk National University, Cheongju 28644, Korea
Water 2021, 13(22), 3170; https://doi.org/10.3390/w13223170
Submission received: 4 September 2021 / Revised: 4 November 2021 / Accepted: 5 November 2021 / Published: 10 November 2021
(This article belongs to the Special Issue Hydrology in Water Resources Management)

Abstract

:
Flood routing can be subclassified into hydraulic and hydrologic flood routing; the former yields accurate values but requires a large amount of data and complex calculations. The latter, in contrast, requires only inflow and outflow data, and has a simpler calculation process than the hydraulic one. The Muskingum model is a representative hydrologic flood routing model, and various versions of Muskingum flood routing models have been studied. The new Muskingum flood routing model considers inflows at previous and next time during the calculation of the inflow and storage. The self-adaptive vision correction algorithm is used to calculate the parameters of the proposed model. The new model leads to a smaller error compared to the existing Muskingum flood routing models in various flood data. The sum of squares obtained by applying the new model to Wilson’s flood data, Wang’s flood data, the flood data of River Wye from December 1960, Sutculer flood data, and the flood data of River Wyre from October 1982 were 4.11, 759.79, 18,816.99, 217.73, 38.81 (m3/s)2, respectively. The magnitude of error for different types of flood data may be different, but the error may be large if the flow rate of the flood data is large.

1. Introduction

Water resources from rivers are sources of hydroelectric power generation, agricultural water, and industrial water; however, owing to the large volumes of water, such rivers are prone to floods that have adverse impacts on life and property [1]. To reduce or prevent such damage, engineering measures, such as the construction of flood control dams or flood walls (levees), are necessary. Therefore, the evaluation of engineering measures for flood control is critical, and these measures are generally directly related to flood routing. Flood routing can be defined as a procedure for determining the flood hydrograph at a point downstream from the base flood hydrograph at an upstream point. In other words, flood routing is the process of determining the amount by which a flood wave is reduced and how long it takes for a flood wave to pass through an arbitrary section of a river based on the amount of storage in that section.
There are two types of flood routing methodologies: hydraulic and hydrologic [2]. Hydraulic flood routing is a method for solving the partial differential continuity and momentum equations, the governing equations of an unsteady nonuniform flow, which hydraulically represent the flow of the flood wave according to the initial and boundary conditions [3]. In contrast, the hydrologic flood routing method yields an approximate solution using the storage equation based on the continuity equation of the flood wave [2]. The hydrologic flood routing method can be divided into three categories: reservoir routing, channel routing, and watershed routing. Channel routing allows the measurement of the storage effect of natural rivers on flood waves by calculating how the discharge of a flood changes as it progresses downstream and provides a standard hydrologic quantity for river planning. The Muskingum flood routing model is a representative channel-routing model [2].
The first Muskingum flood routing model proposed was the linear Muskingum flood routing model (LMM) with two parameters [4]. However, the LMM did not include lateral inflow, and a new Muskingum flood routing model with three variables (LMM-L) was thus proposed [5]. Additionally, a study using the nonlinear Muskingum flood routing model (NLMM) considered the nonlinear relationship between storage and outflow as means of improving upon the LMM [6]. Two types of NLMMs determined by the location of nonlinear factor have been proposed to calculate the storage [7]. The Broyden–Fletcher–Goldfarb–Shanno technique based on a mathematical gradient was applied to the NLMM [8]. NLMM incorporating lateral flow (NLMM-L) was developed to complement the existing NLMM [9]. In 2018, a Muskingum flood routing model called the advanced NLMM (ANLMM-L) was used for calculating a continuous inflow [2]. ANLMM-L is a type of nonlinear Muskingum flood routing model that considers lateral inflow and continuous flow with time. Generalized storage equations for the NLMM have also been suggested to apply more degrees of freedom in the suggested model [10].
In addition to the aforementioned studies, various other investigations have focused on recalculating the error between the outflow from flood data and the calculated outflow. Various studies on Muskingum flood routing models were conducted before the 2000s. The two parameters of the LMM including nonlinear relation between the storage and weighted flow were determined using the least-squares method [11]. The least-squares method was used to adjust the two parameters of LMM, K and X. The Muskingum parameter estimation/flood routing system was developed for linear LMMs and NLMMs and their results have been compared [12]. The results of the two different Muskingum flood routing models were compared. The genetic algorithm was used to estimate the parameters of the NLMMs [13]. In order to overcome the limitations of traditional methods used for Muskingum flood routing models, the genetic algorithm, a well-known meta-heuristic optimization algorithm, was applied.
Since the 2000s, studies applying various meta-heuristic optimization algorithms to the Muskingum flood routing models have continued. An immune clonal selection algorithm was suggested to improve the convergence speed and it was applied to estimate parameters of the NLMM [14]. The Nelder–Mead simplex algorithm was introduced to improve the usability, and it was used to estimate parameters of the NLMM [15]. Furthermore, the simulated annealing and shuffled frog leaping algorithms were used to estimate the parameters of the Muskingum flood routing model in two benchmark/real case studies, and they were compared with the results of Tung’s method [16]. The honeybee mating optimization algorithm with past convergence speed has also been applied for the parameter estimation of the NLMM [17]. The elitist-mutated particle swarm optimization and improved gravitational search algorithm were applied to estimate the parameters of LMMs and NLMMs [18]. Particle swarm optimization was applied to the parameter estimation of the NLMM with four parameters to fit the multiple-peak hydrographs [19], and various NLMMs with different storage calculations such as parameterized initial storage have been proposed using a weed optimization algorithm [20]. Various meta-heuristic optimization algorithms, such as the genetic algorithm, evolution, particle swarm, and a harmony search have been used for parameter estimations of the nonlinear Muskingum model and the variable parameter McCarthy–Muskingum model [21]. The adaptive genetic algorithm was used to estimate the various exponent parameters of the NLMM and it was applied to Wilson’s flood data [22]. In addition, genetic expression programming with faster convergence speed than existing genetic programming was developed for parameter estimation in the Muskingum flood routing model [23]. The water cycle algorithm was applied to estimate the parameters of the NLMM and compared with the genetic algorithm, particle swarm optimization, harmony search, and imperialist competitive algorithm [24]. Although various meta-heuristic optimization algorithms have been tested, studies focusing on comparing the results of each algorithm have indicated limited improvements for the Muskingum flood routing model.
Studies have also been conducted on the application of hybrid meta-heuristic optimization algorithms, combining a charged system search and particle swarm optimization, for parameter estimation of the Muskingum flood routing models [25]. For example, a hybrid meta-heuristic optimization algorithm combining particle swarm optimization and the Nelder–Mead simplex method was used to estimate the parameters of the Muskingum flood routing model [26]. Parameter estimation was conducted using a hybrid meta-heuristic optimization algorithm combining the shuffled frog leaping algorithm and the Nelder–Mead simplex method [27]. The improved real-coded adaptive genetic algorithm and the Nelder–Mead simplex algorithm were combined for the parameter estimation of two improved NLMMs [28]. The hybrid meta-heuristic optimization algorithm applied in particle swarm optimization and the bat algorithm were used to reduce the computational time of the Muskingum flood routing model [29]. Although good results were obtained in some studies, it is difficult to accurately compare them with the results of other existing Muskingum flood routing models because they were calculated using additional variables.
General improvements of the Muskingum flood routing model have also been considered. For example, a modified Muskingum flood routing approach, in conjunction with the HEC-RAS model, was implemented to determine floodplain flows [30]. A new NLMM with four parameters has been suggested [31]. A parameter estimation method of the Muskingum flood routing model in ungagged channel reaches has also been suggested [32].
Studies have been conducted to apply the hybrid method to Muskingum flood routing models. A hybrid harmony search combined with local search algorithm such as Broyden–Fletcher–Goldfarb–Shanno technique was developed and was applied to estimate parameters in NLMM [33]. The new hybrid optimization technique was suggested by combining the modified honeybee mating optimization and generalized reduced gradient algorithm for the application of the new Muskingum model with six parameters [34]. The particle swarm optimization hybridized with Nelder–Mead simplex method was proposed to improve precision and convergence speed in Muskingum model [26]. The hybrid algorithm combining the shuffled frog leaping algorithm and Nelder–Mead simplex was applied to NLMM with four parameters and NLMM with five parameters, and it was compared with the genetic algorithm-generalized reduced gradient [27]. The improved NLMM was suggested for flood prediction using the hybrid algorithm of particle swarm optimization and bat algorithm and it was compared with particle swarm optimization and bat algorithm [29]. The parameters in the two types of NLMM were estimated to improve precision using the hybrid algorithm combining the improved real-coded adaptive genetic algorithm and the Nelder–Mead simplex [28]. Most of the previously proposed hybrid methods combine optimization algorithms, but the hybrid method of this study is a method combining the inflow at the previous time and the inflow at the next time in the Muskingum flood routing. The honey bee mating optimization algorithm was combined with the generalized reduced gradient algorithm to estimate parameters of improved Muskingum flood routing model and applied to the single and multi-peak flood hydrographs [35].
In this study, a new Muskingum flood routing model was suggested. The new Muskingum flood routing model, which considers continuous inflow at previous and next time in the storage and inflow calculations, can enable accurate flood routing in various flood data. The self-adaptive vision correction algorithm (SAVCA), a recently developed meta-heuristic optimization algorithm, was applied to calibrate various parameters in the new Muskingum flood routing model. SAVCA can overcome the disadvantages of the previously developed vision correction algorithm (VCA). When applied to mathematical benchmark functions and water distribution problems, SAVCA has previously displayed good performance [36]. Various meta-heuristic optimization algorithms as well as SAVCA can be applied to the new Muskingum flood routing model to show good results. In the previous study, the type of meta-heuristic optimization algorithm did not significantly affect the results of Muskingum flood routing models [37].

2. Materials and Methodologies

2.1. Overview

The two primary methods used in this study were the SAVCA and new Muskingum flood routing model. The error between the flood outflow data and the calculated outflow in the new Muskingum flood routing model was used as an objective function in the SAVCA, which applied an iterative calculation to minimize the error. The iterative calculation progresses as follows:
  • A group of initial solutions is generated by a random value determined between the lower and upper boundaries for each variable of the new Muskingum flood routing model.
  • One among a group of existing solutions is then selected, or a new solution is generated according to the selected probability.
  • The inflow, storage, and outflow are calculated according to the generated solution, and the error between the flood outflow data and calculated outflow is determined as the objective function.
  • The error is calculated using the sum of squares (SSQ), the Nash–Sutcliffe efficiency (NSE), and the root mean square error (RMSE).
The solution refers to the values of the parameters for Muskingum flood routing models. The calculation in the SAVCA is as follows: If all initial solutions are calculated according to the new Muskingum flood routing model process, the errors of the initial solutions are calculated and sorted in ascending order. SAVCA consists of two types of parameters: self-adaptive and fixed. Division rate 1 (DR1), division rate 2 (DR2), and the compression factor (CF) are self-adaptive parameters. The modulation transfer function rate (MR) and astigmatic rate (AR) are fixed parameters.
DR1 determines whether a new solution should be generated in the range of each variable (global search) or if one solution should be selected from the existing solution group (local search). If the generation of a new solution is determined in DR1, the positive and negative direction searches are determined by DR2. The decision variables of the new solution, generated by global search or selected by local search, are adjusted in detail by MR, CF, and AR. After generating a new solution, the calculation process of the new Muskingum flood routing model is applied. The error (SSQ) is calculated for the inflow, storage, and outflow during each time period. If the error of the new solution is lower than that of the worst solution among the existing solution groups, the new solution is included in the existing solution group. DR1 and DR2 are adjusted according to the calculation process for the new solution. All processes are repeated until a certain number of iterations of SAVCA.
The SSQ was used to calculate the first error value for the Muskingum flood routing models in this study. The SSQ between the observed and calculated outflows was used as the objective function in the optimization process. In the new Muskingum flood routing model, eight parameters were used as decision variables, and the objective function is shown in Equation (1).
Minimize   SSQ = ( O o O s ) 2
where Oo is the observed outflow (m3/s), and Os is the calculated outflow (m3/s). The NSE was used to calculate the second error value for the Muskingum flood routing models in this study. The equation of NSE is shown in Equation (2).
NSE = 1 i = 1 n ( O o O s ) 2 i = 1 n ( O o O ¯ o ) 2
where O ¯ o is the average of observed outflow (m3/s), and n is the number of data. The RMSE was used to calculate the third error value for the Muskingum flood routing models in this study. The equation of RMSE is shown in Equation (3).
RMSE = i = 1 n ( O o O s ) 2 n

2.2. New Muskingum Flood Routing Model

The initial LMM was calculated by assuming the amount of storage in the channel as the sum of prism storage and wedge storage. Prism storage is proportional to the outflow, and wedge storage is proportional to the difference between the inflow and outflow. In the LMM, storage is calculated using Equation (4).
S t = K [ X I t + ( 1 X ) O t ]
where St, It, and Ot are the storage, inflow, and outflow at time t, respectively, and X is the weighted factor. In the NLMM, the nonlinear factor is added in the exponential form of Equation (4). Equation (5) represents the storage in an NLMM.
S t = K [ X I t + ( 1 X ) O t ] m
where m is a nonlinear factor, namely different from 1. The storage calculation in the new Muskingum flood routing model is applied by considering an existing generalized storage. The storage is calculated by considering not only the inflow at the current time (t) but also the inflow at the next time point (t + 1). The reason for considering the inflow at t + 1 instead of t − 1 in the storage calculation is as follows. In the nonlinear Muskingum flood routing models, it is assumed that the storage at time t depends on the upstream storage Sin, and the downstream storage Sout. Inflow (I), outflow (O), Sin, Sout were organized as follows according to water depth [38]. I and Sin are shown in Equations (6) and (7).
I = a 1 y c 1
S i n = b 1 y d 1
where y is the water depth and a1, b1, c1, d1 are coefficients. If c1 and d1 are equal, then Sin is shown in Equation (8).
S i n = b 1 ( I a 1 )
O and Sout are shown in Equations (9) and (10).
O = a 2 y c 2
S o u t = b 2 y d 2
where a2, b2, c2 and d2 are coefficients. If c2 and d2 are equal, then Sout is shown in Equation (11).
S o u t = b 2 ( O a 2 )
In NLMM, inflow, outflow and storage are assumed to be water depth related. The storage can be summarized in Equation (12) [38].
S = X S i n + ( 1 X ) S o u t
The storage with K = b1/a1 = b2/a2 can be expressed as Equation (13).
S = K X I + K ( 1 X ) O
Nonlinear parameter m was applied to Equation (11) and it can be expressed as Equation (14).
S = K [ X I + ( 1 X ) O ] m
Additionally, it is assumed that there is an interdependence between the storage at time t and storage at time t + 1 in generalized storage [10].
S t = X 1 S i n ,   t + X 2 S i n , t + 1 + ( 1 X 1 X 2 ) S o u t , t
Equation (16) can be rearranged as in the process from Equations (12)–(14) and it represents storage in the new Muskingum flood routing model.
S t = K [ X 1 I t + X 2 I t + 1 + ( 1 X 1 X 2 ) O t ] m
where X1 is the weighted factor at time t, and X2 is the weighted factor at time t + 1. In addition, It is the inflow at time t, and It+1 is the inflow at time t + 1. Based on Equation (4), the outflow calculation is summarized in Equation (17).
O t = 1 ( 1 X 1 X 2 ) ( S t K ) 1 m X 1 ( 1 X 1 X 2 ) I t X 2 ( 1 X 1 X 2 ) I t + 1
A weighted inflow, including a continuous inflow has been proposed previously [2]. In this study, the inflow at time t + 1 was included to consider the additional continuous inflow, and the weighted inflow could be calculated as shown in Equation (18).
W t = [ ( 1 θ 1 θ 2 θ 3 ) I t + θ 1 I t 1 + θ 2 I t 2 + θ 3 I t + 1 ]
where Wt is the weighted inflow at time t, and θ1 is the weighted factor of the previous inflow at time t − 1. In addition, θ2 is the weighted factor of the previous inflow at time t − 2, and θ3 is the weighted factor of the next inflow at time t + 1. In previous studies, the inflow at time t − 1 (It−1) and the inflow at time t − 2 (It−2) were considered [2,9]. In this study, all inflows before and after the current time were considered by including the inflow at time t + 1 (It+1). If the weighted inflow of Equation (18) is substituted into Equation (17), the outflow is calculated using Equation (19).
O t = 1 ( 1 X 1 X 2 ) ( S t K ) 1 m X 1 ( 1 X 1 X 2 ) W t X 2 ( 1 X 1 X 2 ) W t + 1
where Wt+1 is the weighted inflow at time t + 1. The storage at time t + 1 can be calculated from the outflow in Equation (19) and the observed inflow. The general storage equation is calculated as Equation (20).
d S d t = I t O t
Equation (21) shows the modified storage equation that considers a change in lateral flow.
d S d t = S t = ( 1 + β ) I t O t
where β is the parameter accounting for the lateral flow. The storage at time t + 1 is shown in Equation (22).
S t + 1 = S t + S
Equation (23) represents the storage at time t + 1, and it can be obtained by substituting Equation (21) into Equation (22).
S t + 1 = S t + [ ( 1 + β ) I t O t ] t
where St+1 is the storage at time t + 1. The outflow in the new Muskingum flood routing model is calculated using eight variables, i.e., K, X1, X2, m, β, θ1, θ2, and θ3. The initial Muskingum flood routing model is based on mass conservation. Therefore, the new Muskingum flood routing model was calculated based on the mass conservation.

2.3. Self-Adaptive Vision Correction Algorithm

SAVCA has a total of six parameters: DR1, DR2, MR, CF, AR, and AF. Among these, DR1, DR2, and CF are self-adaptive, and MR, AR, and AF are fixed. Table 1 shows the parameter types of SAVCA [36].
In SAVCA, the initial decision variables and decision variables generated by the global search are randomly generated within the range between the upper and lower boundaries. Decision variables in the global search are between the current optimal value and the upper boundary or between the lower boundary and the current optimal value based on the probability of DR2. The decision variable generated between the current optimal decision variable and the upper boundary by a global search is shown in Equation (24).
x n = x b + r a n d o m ( 0 ,   1 ) × ( b u x b )
where xn is the new decision variable, and xb is the current optimal value. In addition, random(0, 1) is a random value generated from 0 to 1, and bu is the upper boundary. The decision variable generated between the lower boundary and the current optimal decision variable through a global search is shown in Equation (25).
x n = b l + r a n d o m ( 0 ,   1 ) × ( x b b l )
where bl is the lower boundary. In SAVCA, each decision variable is adjusted by the parameters used in the local search, such as MR, CF, AR, and AF. Equation (26) shows the calculation of the new decision variable.
x n = x n × { 1 + M T F × r a n d o m ( 1 ,   1 ) × ( 1 c u r r e n t   i t e r a t i o n t o t a l   i t e r a t i o n ) C F }
where MTF is the calculated modulation transfer function value. The term random(−1, 1) is a random value between −1 and 1. In addition, CF is a parameter for lens compression. The calculation of MTF is based on the distance (dx) between the current best decision variable and the selected decision variable. dx can be calculated using Equation (27).
d x = r a n k ( x s ) r a n k ( x b ) r a n k ( x l ) r a n k ( x b )
where dx is the relative distance between each decision variable, rank(xs) is the fitness rank of the selected decision variable (xs), rank(xb) is the fitness rank of the best decision variable (xb), and rank(xl) is the fitness rank of the worst decision variable (xl). The fitness rank is the order in which the values of the objective function are sorted. The worst decision variable has the lowest fitness rank. Figure 1 shows the relative distance of dx.
The calculation of the MTF by applying dx is shown in Equation (28).
M T F s = ( d x s ( i = 1 k d x i 2 ) 0.5 ) 0.5
where MTFs is the MTF value of the selected decision variable, and k is the total number of decision variables. In addition, dxs is the relative distance of the selected decision variable, and dxi is the relative distance of the i-th decision variable. The CF in SAVCA can be calculated as shown in Equation (29).
C F = 10 × { s t a n d a r d   d e v i a t i o n ( x i ) a v e r a g e ( x i ) }
where xi is the i-th decision variable. The probability of applying the astigmatism correction process is determined using the AR. The new decision variable adjusted by the application of the astigmatism correction process on a local search is shown in Equation (30).
x n = x n × { 1 + r a n d o m ( 1 ,   1 ) × sin 2 ( A F ) }
where AF is the angle of the astigmatic axis. The application process of SAVCA is summarized as follows: (1) generation of an initial solution group, (2) calculation of the fitness of the initial solution groups, (3) generation of a new solution, (4) application of MR and AR, (5) decision to replace after comparing the new solution with the worst solution in the existing solution group, and (6) repeating (2)–(5) until the total number of iterations is reached. The worst solution is the solution with the lowest fitness rank among the existing solution group. Figure 2 shows the application process.
An initial solution group was created to apply SAVCA to the new Muskingum flood routing model. The weighted inflow, storage, outflow, and SSQ were then calculated. According to the probability of DR1, a new solution was generated by a global search, or one of the existing solution groups was selected through a local search. When creating a new solution with a global search, a new decision variable was created in the positive and negative directions based on the current best decision variable. Each new decision variable was corrected using the MR and AR. In addition, the weighted inflow, storage, outflow, and SSQ were calculated using the new solution with new decision variables. Whether a new solution should be added to the existing solution group was determined by comparing the SSQ of the new solution with the SSQ of the worst solution among the existing group of solutions.

2.4. Flood Data

Five types of flood data were applied to the Muskingum flood routing models: Wilson’s flood data, Wang’s flood data, flood data for River Wye December in 1960, Sutculer flood data, and flood data for River Wyre October in 1982 [5,39,40,41]. All flood data used in this study have been applied in several Muskingum flood routing models in existing studies. The most important aspect of the Muskingum flood routing model is the range of each parameter. The range of each parameter in the new Muskingum flood routing model applied to the five flood datasets is presented in Table 2.
Various Muskingum flood routing models, i.e., LMM, LMM-L, NLMM, NLMM-L, ANLMM-L, and the new Muskingum flood routing model, were compared herein. The parameters used in each Muskingum flood routing model are listed in Table 3.
Differences were observed in the results of the Muskingum flood routing models proposed in other studies. However, the simulation used in this study was conducted according to the parameters listed in Table 3. The data values from existing studies were considered only up to two decimal points when using them as an input for the models in this study. The parameters of Muskingum flood routing models without results from previous studies were obtained by applying SAVCA. However, the values of each parameter were all calculated differently.

3. Application and Results

The first flood dataset used for the application of all Muskingum flood routing models was Wilson’s flood data. The parameters of the LMM for Wilson’s flood data were determined to be 29.164640 for K, and 0.118200 for X1. The parameters of the new Muskingum flood routing model for Wilson’s flood data were determined to be 0.943442 for K, 0.340333 for X1, −0.00102 for X2, 1.744439 for m, −0.02166 for β, 0.758873 for θ1, 0.230779 for θ2, and 0.047773 for θ3. The results, including those obtained using the new Muskingum flood routing model, are compared in Table 4.
Among the results in Table 4, those of the LMM-L, NLMM, NLMM-L, and ANLMM-L were calculated in previous studies [2,5,9,42]. The results of LMM and new Muskingum flood routing model were calculated using Wilson’s flood data by applying SAVCA. Notably, the results of the LMM are better than those of the LMM-L. This is because the results of the LMM-L are the results of a previous study wherein the optimization method was not used, while the results of the LMM were obtained using SAVCA. It should be noted that the results of the new Muskingum flood routing model were better than those of the LMM, LMM-L, NLMM, NLMM-L, and ANLMM-L because the new Muskingum flood routing model showed the smallest error in the initial part from 0 to 24 h and showed the smallest error in the overall results. Because the errors of the NLMM-L and ANLMM-L were small, the new Muskingum flood routing model did not lead to a substantial improvement.
Among the existing Muskingum flood routing models, ANLMM-L showed the best results (smallest error). The difference in SSQ between the ANLMM-L and new Muskingum flood routing model was 0.43 (m3/s)2. The differences between the results of the two models were most notable from 6 to 18 h and from 108 to 126 h. The new Muskingum flood routing model showed the closest outflow to the observed outflow from 6 to 18 h. Among other Muskingum flood routing models, the outflows obtained from LMM-L at 6 h, LMM at 12 h, and NLMM-L at 18 h were close to the observed outflow. The difference between the observed outflow and the outflow obtained from the new Muskingum flood routing model was smaller than that attained using other Muskingum flood routing models.
The second flood dataset was Wang’s flood data. The parameters of the LMM-L for Wang’s flood data were determined to be 1.075331 for K, −0.762101 for X1, and −0.003024 for β. The parameters of the new Muskingum flood routing model for Wang’s flood data were determined to be 0.079266 for K, −1.49742 for X1, 0.011592 for X2, 1.360300 for m, −0.000450 for β, 0.421275 for θ1, 0.044483 for θ2, and 0.261537 for θ3. The results using Wang’s flood data, including those for the new Muskingum flood routing model, are compared in Table 5.
The results of LMM, NLMM, NLMM-L, and ANLMM-L were calculated in previous studies [2,8,9,39]. The results of LMM-L and new Muskingum flood routing model were calculated by applying SAVCA. Notably, the results improved dramatically when the lateral inflow was considered, indicated by the difference between the results of the LMM and LMM-L and between those of the NLMM and NLMM-L. A difference between the results of the ANLMM-L and new Muskingum flood routing model was also observed, although it was not due to a lateral inflow but to differences in the calculation equations of the weighted inflow and storage. The results of the new Muskingum flood routing model were overwhelmingly better than those of other Muskingum flood routing models because the new Muskingum flood routing model showed the smallest error in the latter part from 19 to 29 h and showed the smallest error in the overall results.
Among the existing Muskingum flood routing models, the ANLMM-L showed the best results (smallest error). The difference in SSQ between the ANLMM-L and new Muskingum flood routing model was 149.56 (m3/s)2. The difference between the two results was clear from 228 (19) to 288 (24) h. The differences in the outflow obtained using the Muskingum flood routing models and the observed outflow was small. The outflow calculated by the new Muskingum flood routing model was closest to the observed outflow. Among other Muskingum flood routing models, the outflow obtained using NLMM at 21 h was close to the observed outflow. The difference between the observed outflow and the outflow obtained using the new Muskingum flood routing model was smaller than that when using other Muskingum flood routing models.
The third flood dataset was the flood data of River Wye December in 1960. The parameters for the LMM using these data were determined to be 23.877307 for K, and 0.153174 for X1. The parameters for the new Muskingum flood routing model using these data were determined to be 2.318963 for K, 0.499999 for X1, 0.000390 for X2, 1.359406 for m, 0.057839 for β, 0.805567 for θ1, 0.233550 for θ2, and 6.88×10−11 for θ3. The results of all the models for these data are compared in Table 6.
LMM-L, NLMM, NLMM-L, and ANLMM-L results were calculated in previous studies [2,5,9,42]. The results of LMM and new Muskingum flood routing model were calculated by applying SAVCA. Table 6 displays a notable difference between the results of the LMMs and NLMMs; the difference was large because the error in the flood data of River Wye from December 1960 was large. The new Muskingum flood routing model results were better than those of other Muskingum flood routing models because the new Muskingum flood routing model showed the smallest error from 102 to 198 h including the peak value and showed the smallest error in the overall result.
Among other Muskingum flood routing models, ANLMM-L showed the best results (smallest error). The difference in SSQ between the ANLMM-L and new Muskingum flood routing model was 1677.99 (m3/s)2. The greatest difference between the models occurred at 138–186 h. Among other Muskingum flood routing models, the outflow from ANLMM-L was close to the observed outflow at 180 and 186 h. The difference between the observed outflow and the outflow for the new Muskingum flood routing model was smaller than that for other Muskingum flood routing models.
The fourth flood dataset used was Sutculer flood data which is a flood data with a double-peak. The parameters of the LMM for Sutculer flood data were determined to be 1.0 for K, and −0.006097 for X1. The parameters of the LMM-L for Sutculer flood data were determined to be 1.0 for K, −0.025914 for X1, and −0.041042 for β. The parameters of the NLMM for Sutculer flood data were determined to be 1.0 for K, −0.053787 for X1, and 1.002498 for m. The parameters of the new Muskingum flood routing model for Sutculer flood data were determined to be 0.931599 for K, −0.092988 for X1, 0.009066 for X2, 1.000013 for m, −0.036144 for β, 0.817639 for θ1, 0.214801 for θ2, and 0.745272 for θ3. The results of all models are shown in Table 7.
Among the results in Table 7, those of NLMM-L and ANLMM-L were calculated in previous studies [2,9]. The results of LMM, LMM-L, NLMM, and new Muskingum flood routing model were calculated by applying the SAVCA. Notably, the models that consider the lateral inflow show good results. In addition, the results of LMM-L, NLMM-L, and ANLMM-L were better than those of LMM and NLMM because LMM-L, NLMM-L, and ANLMM-L showed relatively small errors in the overall results. Moreover, the difference between the results of the new Muskingum flood routing model and other Muskingum flood routing models was substantial.
The ANLMM-L showed the best results (smallest error) among the considered models. The difference in SSQ between the ANLMM-L and new Muskingum flood routing model was 63.22 (m3/s)2. The time required to show the difference between the two results ranged from 1 to 3 h. At 1 h, the outflow of most Muskingum flood routing models was close to the observed outflow. At 2 and 3 h, except for the new Muskingum flood routing model, the outflow of most Muskingum flood routing models showed a difference from the observed outflow.
The last flood dataset analyzed was the flood data of River Wyre October in 1982. The parameters of the LMM for the flood data of River Wyre from October 1982 were determined to be 3.950351 for K and 0.295668 for X1. The parameters of the NLMM for the flood data of River Wyre from October 1982 were determined to be 8.248204 for K, 0.284338 for X1, and 0.812821 for m. The parameters of the new Muskingum flood routing model for the flood data of River Wyre from October 1982 were determined to be 0.931599 for K, −0.092988 for X1, 0.009066 for X2, 1.000013 for m, −0.036144 for β, 0.817639 for θ1, 0.214801 for θ2, and 0.745272 for θ3; the results, and their comparison with those of the other models, are shown in Table 8.
Of the results given in Table 8, the LMM-L, NLMM-L, and ANLMM-L results were calculated in previous studies [2,5,9], while those of the LMM and NLMM were calculated by applying SAVCA; the results of the latter two models were the same. The errors were the greatest for the results of both LMM and NLMM, and some of the calculated outflow values obtained were negative. The new Muskingum flood routing model results in Table 8 were also calculated by applying SAVCA.
Notably, the results obtained for the models considering the lateral inflow were good; those of the LMM-L, NLMM-L, and ANLMM-L were significantly better than those of the LMM and NLMM. The difference in the results of all Muskingum flood routing models occurs from 26 to 31 h. The outflows of NLMM-L, ANLMM-L, and new Muskingum flood routing model were the closest to the observed outflow from 26 to 31 h. Although the results of the NLMM-L and ANLMM-L were similar to the observed outflow at 26 h, these differed over time. However, the new Muskingum flood routing model results did not differ significantly from the observed outflow from 26 to 31 h. ANLMM-L showed the best results (smallest error) among the existing Muskingum flood routing models. The difference in SSQ between the ANLMM-L and new Muskingum flood routing model was 1.35 (m3/s)2. Although varying results were obtained when using different flood data for the various models, the overall results of new Muskingum flood routing model were better than those of other Muskingum flood routing models because the error of new Muskingum flood routing model was relatively small in the latter part from 16 to 31 h.
The new flood data in Daechung were applied to calibrate and validate the new Muskingum flood routing model. The flood data in April, 2014 were used for calibration and the flood data in April, 2018 were used for validation. Among various Muskingum flood routing models, LMM-L considering lateral inflow to LMM, NLMM considering nonlinearity to LMM and new Muskingum flood routing model were applied to Daechung flood data and compared.
The parameters of LMM were 3.989981 for K and −0.034950 for X. The parameters of LMM-L were 3.865970 for K, −0.043293 for X, and −0.020080 for β. The parameters of NLMM were 2.225924 for K, −1.5 for X, and 1.0 for m. The parameters of NLMM-L were 2.225076 for K, −1.406061 for X, 1.000000 for m, −0.007072 for β, 1.000000 for θ. The parameters of ANLMM-L were 2.123526 for K, −1.500000 for X, 1.000000 for m, −0.010223 for β, 1.000000 for θ1, and 0.017871 for θ2. The parameters of new Muskingum flood routing model were 2.220910 for K, −1.498329 for X1, 0.094832 for X2, 1.000008 for m, −0.012616 for β, 0.999660 for θ1, 0.000093 for θ2, and 0.000288 for θ3. Each parameter was applied equally for 2014 and 2018 data. A total of 100 simulations were conducted for each Muskingum flood routing model, yielding the best results. Table 9 shows the results of Daechung flood data in 2014.
Based on the variable values determined from the 2014 flood data, it was applied to the 2018 flood data. Table 10 shows the results of Daechung flood data in 2018.
Figure 3 shows the results of calibration and validation in Daechung flood data.
In the 2014 flood data, the SSQs of LMM, LMM-L, NLMM, NLMM-L, ANLMM-L and new Muskingum flood routing model were 88.23, 73.81, 43.79, 42.32, 40.16, and 39.55, respectively. In the 2014 flood data, the error of new Muskingum flood routing model was relatively small and the calibration of new Muskingum flood routing model was relatively accurate. The SSQs of new Muskingum flood routing model in the 2018 flood data was relatively small. In the 2018 flood data, the SSQs of LMM, LMM-L, NLMM, NLMM-L, ANLMM-L and new Muskingum flood routing model were 221.92, 180.41, 171.45, 161.99, 159.84, and 157.64, respectively. In the results of Figure 2, more accurate flood routing was performed by applying the new Muskingum flood routing model compared to LMM, LMM-L, NLMM, NLMM-L and ANLMM-L.

4. Discussion

Because the calculation process differs for each Muskingum flood routing model, the time required to find the parameters when applying a meta-heuristic optimization algorithm is different for each method. The time required to apply SAVCA was summarized to determine the parameters of each Muskingum flood routing model for Wilson’s flood data. The parameters of SAVCA were set at a constant, and the simulation was conducted 10 times. In addition, the number of iterations was set to 100,000. Table 11 shows the time taken for SAVCA when using Wilson’s flood data.
Depending on the parameters of the SAVCA and flood data, the results over time using the Muskingum flood routing models differ from the results in Table 9. As the number of parameters of the Muskingum flood routing models increased, the time required also increased. The time required by LMMs, LMM and LMM-L, was the shortest. NLMM and NLMM-L required a greater amount of time compared to the LMMs. The new Muskingum flood routing model required more time than the ANLMM-L, which in turn required more time than the NLMM-L. The time required by the new Muskingum flood routing model was approximately 1.5-times that required by LMM and LMM-L. In conclusion, the new Muskingum flood routing model produced more accurate results but took more time owing to the greater number of parameters and calculations.
An analysis was conducted on how each parameter of new Muskingum flood routing model affects the results. Daechung flood data in April, 2014 was applied to analyze the sensitivity of each parameter in the new Muskingum flood routing model. Figure 4 showed the results of the sensitivity analysis for the parameters in the new Muskingum flood routing model.
SSQ decreases and then increases as parameter K increases, and SSQ increases as parameter X1 increases. SSQ increases as parameter X2 increases and SSQ increases as parameter m increases. SSQ decreases and then increases as parameter β increases, and SSQ decreases as the parameter θ1 increases. SSQ increases as parameter θ2 increases, and SSQ increases as parameter θ3 increases. As each parameter changed, the change of the results was not constant. It is difficult to find a uniform pattern in all the results. However, what can be confirmed from the results of sensitivity analysis is that a meta-heuristic optimization algorithm such as SAVCA is required to produce results with a low SSQ.

5. Conclusions

The Muskingum flood routing model is a representative hydrologic flood routing model that is widely used owing to its easy applicability. The proposed Muskingum flood routing model in this study is a simple model that can be applied by researchers that use the existing Muskingum models for accurate flood routing.
In this study, the new Muskingum flood routing model was applied to various flood data, and the results obtained were compared with those of previously developed Muskingum flood routing models. As an index for comparison, the error was calculated between the observed and simulated outflows using the SSQ, NSE and RMSE. In addition, SAVCA, a meta-heuristic optimization algorithm, was applied to adjust the parameters of the new Muskingum flood routing model.
In the sensitivity analysis, the changes of the eight parameters in the new Muskingum flood routing model are different. There are parameters whose results are improved as the value (θ1) increases, some parameters (m, θ2, θ3, x1, x2) whose results are improved as the value decreases, and some parameters (K, β) whose results are changed (improved and then deteriorated) as the value increases. The eight parameters of the new Muskingum flood routing model are decision variables of SAVCA and are calculated through the optimization process.
Muskingum flood routing models considering the lateral inflow are capable of relatively sophisticated simulations, which corroborates that the influence of lateral inflow on the outflow can be considered. Among the existing models, the ANLMM-L showed the highest accuracy, although the difference between its results and those of the other Muskingum flood routing models was insignificant. In the new Muskingum flood routing model, the improved calculation method of the inflow at previous time and next time reflected the trend of the observed outflow.
Since the Muskingum flood routing model proposed in this study has eight parameters, the calculation process is more complicated than the existing Muskingum flood routing models. Accurate flood prediction is possible due to the complicated calculation process, but the calculation time is long.
Many studies, including this study, have been performed by applying various meta-heuristic optimization algorithms to the Muskingum flood routing models. Since various optimization algorithms cannot show advantages in all problems, the results can be improved by appropriately selecting the operators of each meta-heuristic optimization algorithm. In addition, deep learning techniques have been widely used to apply various flood prediction methods including Muskingum flood routing models. By replacing the optimizer in deep learning techniques with meta-heuristic optimization algorithms, it would be possible to produce improved results in flood prediction.

Funding

This research was funded by the National Research Foundation (NRF) of Korea (NRF-2019R1I1A3A01059929).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

This work was supported by a grant from The National Research Foundation (NRF) of Korea (NRF-2019R1I1A3A01059929).

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

ANLMM-LAdvanced nonlinear Muskingum flood routing model considering continuous inflow
NLMM-LNonlinear Muskingum flood routing model incorporating lateral flow
NLMMNonlinear Muskingum method
LMM-LLinear Muskingum method incorporating lateral flow
LMMLinear Muskingum method
SAVCASelf-adaptive vision correction algorithm
DR1Division rate 1
DR2Division rate 2
MTFModulation transfer function
CFCompression factor
ARAstigmatic rate
AFAstigmatic angle
SSQSum of squares
NSENash–Sutcliffe efficiency
RMSERoot mean square error

References

  1. Go, I.H. Development of watershed integrated water resource management technique. J. Korea Water Resour. Assoc. 2004, 37, 10–15. [Google Scholar]
  2. Lee, E.H.; Lee, H.M.; Kim, J.H. Development and application of advanced Muskingum flood routing model considering continuous flow. Water 2018, 10, 760. [Google Scholar] [CrossRef] [Green Version]
  3. Absi, R. Reinvestigating the parabolic-shaped eddy viscosity profile for free surface flows. Hydrology 2021, 8, 126. [Google Scholar] [CrossRef]
  4. McCarthy, G.T. The unit hydrograph and flood routing. In Proceedings of Conference of North Atlantic Division; U.S. Army Corps of Engineers: New London, CT, USA, 1938. [Google Scholar]
  5. O’donnell, T. A direct three-parameter Muskingum procedure incorporating lateral inflow. Hydrol. Sci. J. 1985, 30, 479–496. [Google Scholar] [CrossRef] [Green Version]
  6. Tung, Y.K. River flood routing by nonlinear Muskingum method. J. Hydraul. Eng. 1985, 111, 1447–1460. [Google Scholar] [CrossRef] [Green Version]
  7. Das, A. Parameter estimation for Muskingum models. J. Irrig. Drain. Eng. 2004, 130, 140–147. [Google Scholar] [CrossRef]
  8. Geem, Z.W. Parameter estimation for the nonlinear Muskingum model using the BFGS technique. J. Irrig. Drain. Eng. 2006, 132, 474–478. [Google Scholar] [CrossRef]
  9. Karahan, H.; Gurarslan, G.; Geem, Z.W. A new nonlinear Muskingum flood routing model incorporating lateral flow. Eng. Optim. 2015, 47, 737–749. [Google Scholar] [CrossRef]
  10. Bozorg-Haddad, O.; Abdi-Dehkordi, M.; Hamedi, F.; Pazoki, M.; Loáiciga, H.A. Generalized storage equations for flood routing with nonlinear Muskingum models. Water Resour. Manag. 2019, 33, 2677–2691. [Google Scholar] [CrossRef]
  11. Gill, M.A. Flood routing by the Muskingum method. J. Hydrol. 1978, 36, 353–363. [Google Scholar] [CrossRef]
  12. Yoon, J.; Padmanabhan, G. Parameter estimation of linear and nonlinear Muskingum models. J. Water Resour. Plann. Manag. 1993, 119, 600–610. [Google Scholar] [CrossRef]
  13. Mohan, S. Parameter estimation of nonlinear Muskingum models using genetic algorithm. J. Hydraul. Eng. 1997, 123, 137–142. [Google Scholar] [CrossRef]
  14. Luo, J.; Xie, J. Parameter estimation for nonlinear Muskingum model based on immune clonal selection algorithm. J. Hydrol. Eng. 2010, 15, 844–851. [Google Scholar] [CrossRef]
  15. Barati, R. Parameter estimation of nonlinear Muskingum models using Nelder-Mead simplex algorithm. J. Hydrol. Eng. 2011, 16, 946–954. [Google Scholar] [CrossRef]
  16. Orouji, H.; Bozorg Haddad, O.; Fallah-Mehdipour, E.; Marino, M.A. Estimation of Muskingum parameter by meta-heuristic algorithms. In Proceedings of the Institution of Civil Engineers-Water Management; Thomas Telford Ltd.: London, UK, 2013; Volume 166, pp. 315–324. [Google Scholar]
  17. Niazkar, M.; Afzali, S.H. Assessment of modified honey bee mating optimization for parameter estimation of nonlinear Muskingum models. J. Hydrol. Eng. 2015, 20, 04014055. [Google Scholar] [CrossRef]
  18. Kang, L.; Zhang, S. Application of the elitist-mutated PSO and an improved GSA to estimate parameters of linear and nonlinear Muskingum flood routing models. PLoS ONE 2017, 11, e0147338. [Google Scholar] [CrossRef]
  19. Moghaddam, A.; Behmanesh, J.; Farsijani, A. Parameters estimation for the new four-parameter nonlinear Muskingum model using the particle swarm optimization. Water Resour. Manag. 2016, 30, 2143–2160. [Google Scholar] [CrossRef]
  20. Hamedi, F.; Bozorg-Haddad, O.; Pazoki, M.; Asgari, H.R.; Parsa, M.; Loáiciga, H.A. Parameter estimation of extended nonlinear Muskingum models with the weed optimization algorithm. J. Irrig. Drain. Eng. 2016, 142, 04016059. [Google Scholar] [CrossRef]
  21. Perumal, M.; Tayfur, G.; Rao, C.M.; Gurarslan, G. Evaluation of a physically based quasi-linear and a conceptually based nonlinear Muskingum methods. J. Hydrol. 2017, 546, 437–449. [Google Scholar] [CrossRef]
  22. Zhang, S.; Kang, L.; Zhou, B. Parameter estimation of nonlinear Muskingum model with variable exponent using adaptive genetic algorithm. In Environmental Conservation, Clean Water, Air & Soil (CleanWAS); IWA Publishing: London, UK, 2017; pp. 231–237. [Google Scholar]
  23. Bagatur, T.; Onen, F. Development of predictive model for flood routing using genetic expression programming. J. Flood Risk Manag. 2018, 11, S444–S454. [Google Scholar] [CrossRef] [Green Version]
  24. Akbarifard, S.; Qaderi, K.; Alinnejad, M. Parameter estimation of the nonlinear Muskingum flood-routing model using water cycle algorithm. J. Watershed Manag. Res. 2018, 8, 34–43. [Google Scholar] [CrossRef]
  25. Talatahari, S.; Sheikholeslami, R.; Farahmand Azar, B.; Daneshpajouh, H. Optimal parameter estimation for Muskingum model using a CSS-PSO method. Adv. Mech. Eng. 2013, 5, 480954. [Google Scholar] [CrossRef]
  26. Ouyang, A.; Li, K.; Truong, T.K.; Sallam, A.; Sha, E.H.M. Hybrid particle swarm optimization for parameter estimation of Muskingum model. Neural Comput. Appl. 2014, 25, 1785–1799. [Google Scholar] [CrossRef]
  27. Haddad, O.B.; Hamedi, F.; Fallah-Mehdipour, E.; Orouji, H.; Marino, M.A. Application of a hybrid optimization method in Muskingum parameter estimation. J. Irrig. Drain. Eng. 2015, 141, 04015026. [Google Scholar] [CrossRef] [Green Version]
  28. Kang, L.; Zhou, L.; Zhang, S. Parameter estimation of two improved nonlinear Muskingum models considering the lateral flow using a hybrid algorithm. Water Resour. Manag. 2017, 31, 4449–4467. [Google Scholar] [CrossRef]
  29. Ehteram, M.; Binti Othman, F.; Mundher Yaseen, Z.; Abdulmohsin Afan, H.; Falah Allawi, M.; Najah Ahmed, A.; Shaid, S.; Singh, V.; El-Shafie, A. Improving the Muskingum flood routing method using a hybrid of particle swarm optimization and bat algorithm. Water 2018, 10, 807. [Google Scholar] [CrossRef] [Green Version]
  30. O’Sullivan, J.J.; Ahilan, S.; Bruen, M. A modified Muskingum routing approach for floodplain flows: Theory and practice. J. Hydrol. 2012, 470, 239–254. [Google Scholar] [CrossRef] [Green Version]
  31. Easa, S.M. New and improved four-parameter non-linear Muskingum model. In Proceedings of the Institution of Civil Engineers-Water Management; Thomas Telford Ltd.: London, UK, 2014; Volume 167, pp. 288–298. [Google Scholar]
  32. Yoo, C.; Lee, J.; Lee, M. Parameter estimation of the Muskingum channel flood-routing model in ungauged channel reaches. J. Hydrol. Eng. 2017, 22, 05017005. [Google Scholar] [CrossRef]
  33. Karahan, H.; Gurarslan, G.; Geem, Z.W. Parameter estimation of the nonlinear Muskingum flood-routing model using a hybrid harmony search algorithm. J. Hydrol. Eng. 2013, 18, 352–360. [Google Scholar] [CrossRef]
  34. Niazkar, M.; Afzali, S.H. Application of new hybrid optimization technique for parameter estimation of new improved version of Muskingum model. Water Resour. Manag. 2016, 30, 4713–4730. [Google Scholar] [CrossRef]
  35. Niazkar, M.; Afzali, S.H. Parameter estimation of an improved nonlinear Muskingum model using a new hybrid method. Hydrol. Res. 2017, 48, 1253–1267. [Google Scholar] [CrossRef]
  36. Lee, E.H. Application of self-adaptive vision-correction algorithm for water-distribution problem. KSCE J. Civ. Eng. 2021, 25, 1106–1115. [Google Scholar] [CrossRef]
  37. Kim, Y.N. Development of Exponential Bandwidth Harmony Search with Centralized Global Search. Master’s Thesis, Chungbuk National University, Cheongju, Korea, 2021. (In Korean). [Google Scholar]
  38. Chow, V.T. Open-Channel Hydraulics; McGraw-Hill Civil Engineering Series; McGraw-Hill: New York, NY, USA, 1959. [Google Scholar]
  39. Wang, W.; Xu, Z.; Qiu, L.; Xu, D. Hybrid chaotic genetic algorithms for optimal parameter estimation of Muskingum flood routing model. In Proceedings of the International Joint Conference on Computational Sciences and Optimization, Sanya, China, 24–26 April 2009; pp. 215–218. [Google Scholar]
  40. Bajracharya, K.; Barry, D.A. Accuracy criteria for linearised diffusion wave flood routing. J. Hydrol. 1997, 195, 200–217. [Google Scholar] [CrossRef]
  41. Ulke, A. Flood Routing Using Muskingum Method. Master’s Thesis, Suleyman Demirel University, Kaskelen, Kazakhstan, 2003. (In Turkish). [Google Scholar]
  42. Karahan, H.; Gurarslan, G. Modelling flood routing problems using kinematic wave approach: Sutculer example. In Proceedings of the 7th National Hydrology Congress, Isparta, Turkey, 26–27 September 2013. [Google Scholar]
Figure 1. Relative distance of dx.
Figure 1. Relative distance of dx.
Water 13 03170 g001
Figure 2. Application process.
Figure 2. Application process.
Water 13 03170 g002
Figure 3. Results of calibration and validation in Daechung flood data (a) 2014 flood data; (b) 2018 flood data..
Figure 3. Results of calibration and validation in Daechung flood data (a) 2014 flood data; (b) 2018 flood data..
Water 13 03170 g003
Figure 4. Results of sensitivity analysis for the parameters in the new Muskingum flood routing model.
Figure 4. Results of sensitivity analysis for the parameters in the new Muskingum flood routing model.
Water 13 03170 g004
Table 1. Parameter types of SAVCA.
Table 1. Parameter types of SAVCA.
ParametersDR1DR2MRCFARAF
TypesSelf-adaptiveSelf-adaptiveFixedSelf-adaptiveFixedFixed
Table 2. Range of parameters in the new Muskingum flood routing model.
Table 2. Range of parameters in the new Muskingum flood routing model.
ParametersWilson’s Flood DataWang’s Flood DataFlood Data for River Wye December in 1960Sutculer Flood DataFlood Data for River Wyre October in 1982
K0.01–50.000.01–50.000.01–50.000.01–50.000.01–50.00
X1−0.50–0.50−1.50–1.50−0.50–0.50−0.50–0.50−0.50–0.50
X2−0.50–0.50−1.50–1.50−0.50–0.50−0.50–0.50−0.50–0.50
m1.00–3.001.00–3.001.00–3.001.00–3.000.00–1.00
β−0.10–0.10−3.00–3.00−0.10–0.10−0.10–0.10−3.00–3.00
θ10.00–1.000.00–1.000.00–1.000.00–1.000.00–1.00
θ20.00–1.000.00–1.000.00–1.000.00–1.000.00–1.00
θ30.00–1.000.00–1.000.00–1.000.00–1.000.00–1.00
Table 3. Parameters used in each Muskingum flood routing model (○: applied, X: not applied).
Table 3. Parameters used in each Muskingum flood routing model (○: applied, X: not applied).
ParametersLMMLMM-LNLMMNLMM-LANLMM-LThis Study
K
X1
X2
m
β
θ1
θ2
θ3
Table 4. Results when using Wilson’s flood data.
Table 4. Results when using Wilson’s flood data.
Time
(h)
Inflow
(m3/s)
Outflow
(m3/s)
LMM
(m3/s)
LMM-L
(m3/s)
[5]
NLMM
(m3/s)
[42]
NLMM-L
(m3/s)
[9]
ANLMM-L
(m3/s)
[2]
This Study
(m3/s)
0222222.0022.0022.0022.0022.0022.00
6232121.8721.1022.0021.7121.5721.33
12352120.5221.7022.4022.0221.6721.13
18712619.0722.6026.6026.0825.4625.53
241033426.9030.7034.5033.5134.5934.75
301114443.5844.7044.2042.8343.7343.52
361095559.5858.1056.9055.4454.5954.62
421006672.3268.9068.1066.6766.0166.08
48867580.6576.1077.1075.7775.5275.53
54718283.9179.2083.3082.1282.1682.11
60598582.5178.5085.9084.7885.0485.08
66478478.6375.6084.5083.4284.0083.89
72398072.3270.7080.6079.4479.6279.61
78327365.4965.1073.7072.4872.6372.53
84286458.2159.1065.4064.0863.8063.81
90245451.7053.4056.0054.5854.3154.27
96224445.5047.9046.7045.2244.8044.84
102213640.1543.1037.7036.3436.2536.32
108203035.8238.9030.5029.2129.4529.52
114192532.2635.4025.2024.2124.6324.66
120192229.1732.3021.7020.9621.3921.46
126181926.9329.9020.0019.4119.8119.77
SSQ
(m3/s)2
--605.63815.6836.779.824.544.11
Squared root of SSQ (m3/s)--24.6128.566.063.132.132.03
NSE--0.9743220.9743260.9924120.9995830.9998080.999826
RMSE
(m3/s)
--5.3102595.3698852.9194110.6839930.4641240.442254
Table 5. Results when using Wang’s flood data.
Table 5. Results when using Wang’s flood data.
Time
(12 h)
Inflow
(m3/s)
Outflow
(m3/s)
LMM
(m3/s)
[39]
LMM-L
(m3/s)
NLMM
(m3/s)
[8]
NLMM-L
(m3/s)
[9]
ANLMM-L
(m3/s)
[2]
This Study
(m3/s)
1261228228.00228.00228.00228.00228.00228.00
2389300305.19300.19303.80299.74300.92301.75
3462382382.00377.92382.30382.57381.51382.38
4505444442.70440.10442.40442.76443.15442.81
5525490483.60482.17482.40482.16482.69483.63
6543513513.00511.70511.2509.89510.09510.15
7556528534.29532.96532.30530.72530.66530.75
8567543550.44548.97548.50546.77546.62546.79
9577553563.53561.89561.70559.96559.77559.53
10583564573.16571.53571.60569.94569.80569.75
11587573580.02578.38578.70577.07576.95577.89
12595581587.32585.44586.20584.39584.22584.03
13597588592.14590.40591.20589.68589.60589.77
14597594594.59592.93593.90592.34592.30591.61
15589592592.02590.68591.80590.33590.34586.67
16556584574.89574.62575.70574.68574.86576.15
17538566556.85556.15558.50556.41556.23556.07
18516550536.93536.22539.00537.43537.13536.33
19486520512.18511.79514.80513.47513.35521.23
20505504507.96505.60509.60507.07506.51502.72
21477483493.22492.40484.90494.86494.95492.05
22429461462.34462.82464.80464.39464.94463.80
23379420421.87422.73425.10423.97424.15422.09
24320368372.34373.60376.10375.05375.07374.32
25263318318.97320.23322.40321.35321.35322.59
26220271270.39271.06272.50271.42271.40271.68
27182234226.99227.38227.50226.94227.09229.70
28167193197.20196.67195.70194.92195.13194.64
29152178174.87174.28172.60172.46172.76174.61
SSQ
(m3/s)2
--1086.84999.83979.96917.06909.35759.79
Squared root of SSQ (m3/s)--32.9731.6231.3030.2830.1627.56
NSE--0.9982470.9983260.9983590.9984640.9984780.998728
RMSE
(m3/s)
--6.0081115.87116935.8130545.6234235.5987625.118558
Table 6. Results when using flood data of River Wye December in 1960.
Table 6. Results when using flood data of River Wye December in 1960.
Time
(h)
Inflow
(m3/s)
Outflow
(m3/s)
LMM
(m3/s)
LMM-L
(m3/s)
[5]
NLMM
(m3/s)
[42]
NLMM-L
(m3/s)
[9]
ANLMM-L
(m3/s)
[2]
This Study
(m3/s)
0154102102.00102.00102.00102.00102.00102.00
6150140118.15116.00154.00149.50146.52141.89
12219169115.12120.00152.00156.59155.74155.50
18182190152.64147.00181.00191.40194.41185.46
24182209161.35158.00191.00200.79194.19190.53
30192218165.67165.00185.00195.14196.05195.99
36165210178.37176.00187.00197.46198.35196.69
42150194177.11178.00179.00188.48186.83188.20
48128172173.05176.00162.00170.80172.12175.53
54168149152.45164.00141.00148.10150.37157.72
60260136140.42160.00154.00162.59167.56169.06
66471228137.74167.00198.00210.36216.61213.24
72717303192.13218.00264.00281.58294.27287.51
781092366280.05303.00344.00367.75378.29378.89
841145456511.40484.00416.00447.65461.17465.87
90600615797.99690.00599.00629.57612.03609.41
96365830781.75700.00871.00892.78862.51863.65
102277969674.00642.00834.00859.01884.60887.00
108227665565.24572.00689.00719.30737.54730.86
114187519472.11505.00535.00567.50565.33555.56
120161444392.21442.00397.00427.85414.97410.06
126143321326.86386.00283.00308.86297.45300.33
132126208275.37338.00202.00220.90216.14224.40
138115176233.04296.00152.00163.64164.43174.61
144102148200.36260.00124.00131.90134.94143.56
15093125172.80228.00106.00111.93114.46121.64
15688114150.03201.0094.0099.28101.24106.75
16282106132.71179.0088.0092.9094.0097.42
1687697118.75160.0082.0086.1486.9489.67
1747389106.60144.0075.0079.3480.1382.79
180708197.18130.0073.0076.4676.8778.56
186677689.65118.0069.0073.1373.5474.88
192637183.66109.0066.0069.8570.2371.46
198596678.25100.0062.0065.0965.6067.24
SSQ
(m3/s)2
--196,077.12251,802.0037,944.1525,915.2720,494.9818,816.99
Squared root of SSQ (m3/s)--442.81501.80194.79160.98143.16137.18
NSE--0.9166660.9216000.9592080.9889860.9912900.992003
RMSE
(m3/s)
--77.08262574.76575053.93017828.02361224.92107723.879109
Table 7. Results when using Sutculer flood data.
Table 7. Results when using Sutculer flood data.
Time
(h)
Inflow
(m3/s)
Outflow
(m3/s)
LMM
(m3/s)
LMM-L
(m3/s)
NLMM
(m3/s)
NLMM-L
(m3/s)
[9]
ANLMM-L
(m3/s)
[2]
This Study
(m3/s)
07.537.007.007.007.007.007.007.00
19.068.007.597.257.587.247.268.14
228.0023.0010.069.119.949.009.0111.97
379.8025.0029.9527.6629.5627.3527.3525.63
464.3075.0076.0474.9275.8674.8474.8173.93
538.2060.0063.4761.3663.7061.5761.5962.61
641.4040.0039.8437.3339.9637.4037.4137.54
741.3041.0041.3039.6441.3139.6339.6239.37
833.8041.0040.8739.4240.9239.4739.4739.72
932.0032.0034.1032.5534.1532.5732.5832.56
1029.0030.0031.9530.6631.9830.6830.6831.10
1135.0034.0029.5128.0329.4928.0028.0029.48
1263.1035.0036.3034.1036.1033.9333.9336.09
13110.0060.0064.2660.9863.8160.6260.6263.47
14170.00105.00110.82105.81110.12105.25105.25108.08
15216.00160.00169.24162.69168.46162.06162.07157.14
16131.00206.00208.43203.95208.54204.11204.11205.42
17101.00128.00133.73126.88134.55127.58127.61126.32
1865.0097.00100.8196.74101.3397.1497.1098.08
1962.4061.0066.9163.1467.1963.3363.3263.18
2053.8060.0062.1659.7162.2659.7859.7659.18
2136.3050.0053.2751.3753.4451.5151.5051.71
2229.6033.0036.8935.0737.0335.1635.1635.19
2325.0027.0029.7528.4429.8228.4928.4828.49
2421.3023.0025.0624.0025.1124.0324.0324.13
2519.6019.0021.4220.4721.4420.4920.4920.53
2618.0018.0019.6118.8019.6318.8118.8118.90
2717.3017.0018.0517.2818.0617.2917.2917.38
2817.0017.0017.3316.6017.3316.6016.6016.63
2916.0017.0016.9616.2916.9716.2916.2916.53
SSQ
(m3/s)2
--512.87282.89510.18281.11280.95217.73
Squared root of SSQ (m3/s)--22.6516.8222.5916.7716.7614.76
NSE--0.9925570.9958950.9925960.9959210.9959220.996840
RMSE
(m3/s)
--4.1346943.0708024.1238233.0610803.0605932.694028
Table 8. Results when using flood data of River Wyre October in 1982.
Table 8. Results when using flood data of River Wyre October in 1982.
Time
(h)
Inflow
(m3/s)
Outflow
(m3/s)
LMM
(m3/s)
LMM-L
(m3/s)
[5]
NLMM
(m3/s)
NLMM-L
(m3/s)
[9]
ANLMM-L
(m3/s)
[2]
This Study
(m3/s)
02.608.308.308.308.308.308.308.30
14.209.005.588.206.008.518.528.73
212.309.901.688.102.278.799.9410.11
325.4010.200.0012.700.0010.9412.7412.75
424.1018.909.6727.908.6620.2819.7119.51
520.3035.9016.4539.9015.5037.5435.7336.22
623.3051.8016.5845.7016.0249.0748.8749.25
727.7059.4017.1552.2016.9155.1155.9555.83
827.7063.3020.9461.4020.9062.5062.7462.54
926.9069.6023.7168.9023.7871.4471.3571.33
1024.8076.7025.7374.7025.8078.0377.9577.87
1126.9082.0024.5277.2024.5982.0782.6782.68
1233.7085.3022.5279.8022.7783.7285.2785.10
1333.9089.0026.4587.8026.9287.4388.1187.71
1427.8094.6031.6995.5032.0995.4994.7494.61
1520.8098.8033.2397.7033.18100.8899.9099.91
1615.6098.0030.9594.4030.4399.2998.8798.75
1711.9091.8026.9887.9026.2992.0692.0591.82
189.5082.3022.5779.8021.9882.2282.3682.12
197.8072.0018.5971.5018.2471.7571.8871.67
206.5061.9015.2663.6015.1961.9461.9361.80
215.8053.0012.4056.1012.6053.1253.1053.03
225.0045.6010.3749.6010.7445.4745.3745.34
234.8039.208.5243.709.0439.1439.0439.07
244.5033.807.3138.807.8733.7633.6533.68
254.1029.306.4734.607.0329.5529.3929.44
263.7026.205.7830.906.3426.1225.9626.02
273.4023.505.1627.705.7123.2023.0823.14
283.2021.204.6124.805.1420.6720.5920.64
292.9019.204.2322.304.7318.5218.4418.48
302.8017.703.7920.104.2716.7116.6816.72
312.6016.403.5218.203.9615.1215.0915.23
SSQ
(m3/s)2
--53,544.67468.8453,544.9953.6640.1638.81
Squared root of SSQ (m3/s)--231.4021.65231.407.336.346.23
NSE--−0.2139580.989570−0.2139650.9988420.9990900.999120
RMSE
(m3/s)
--40.9056333.79078040.9057551.2635631.1203201.101288
Table 9. Results of Daechung flood data in 2014.
Table 9. Results of Daechung flood data in 2014.
Time
(h)
Inflow
(m3/s)
Outflow
(m3/s)
LMM
(m3/s)
LMM-L
(m3/s)
NLMM
(m3/s)
NLMM-L
(m3/s)
ANLMM-L
(m3/s)
This Study
(m3/s)
00.790.470.470.470.470.470.470.47
32.121.040.750.750.640.650.650.62
63.542.111.791.801.811.801.781.79
950.753.194.664.963.143.122.622.34
12103.9038.3239.9440.4241.8641.6540.8841.01
15112.3390.6886.6986.4389.9989.6889.4389.79
1882.41106.18104.31102.80104.36104.08104.41104.44
2145.0682.3287.1484.8884.2583.9884.3984.04
2423.1251.5955.8353.7052.8052.4652.5852.12
2715.7631.2131.8230.3130.7430.3330.1729.87
3015.0320.8220.1319.2220.6320.2319.9819.86
3317.2217.2416.5015.9717.4117.0716.8616.80
3617.2717.6517.0216.6417.9117.6417.5117.50
3915.0416.0917.1316.7617.5817.3817.3217.30
429.9713.0215.4415.0515.5915.4415.4415.41
456.3810.3911.3410.9811.1611.0411.0410.98
485.678.177.717.437.597.497.457.39
514.457.076.195.996.366.276.236.21
544.235.994.924.774.994.924.894.87
574.184.914.424.304.524.464.434.42
602.254.204.184.074.324.274.274.27
632.334.022.782.692.672.642.632.60
662.243.012.452.382.512.482.462.45
692.113.132.292.242.342.312.302.30
722.833.492.182.142.182.162.142.13
754.254.552.702.672.732.712.682.67
782.834.203.783.723.943.923.923.94
812.152.083.072.992.952.932.942.92
842.131.042.402.332.332.312.302.37
SSQ
(m3/s)2
--88.2373.8143.7942.3240.1639.55
Squared root of SSQ (m3/s)--1.781.621.251.231.201.19
NSE--0.9966120.9971660.9983190.9983750.9984580.998482
RMSE
(m3/s)
--1.7751351.6235771.2505011.2293651.1975841.188440
Table 10. Results of Daechung flood data in 2018.
Table 10. Results of Daechung flood data in 2018.
Time
(h)
Inflow
(m3/s)
Outflow
(m3/s)
LMM
(m3/s)
LMM-L
(m3/s)
NLMM
(m3/s)
NLMM-L
(m3/s)
ANLMM-L
(m3/s)
This Study
(m3/s)
00.530.320.320.320.320.320.320.32
31.861.020.520.520.430.440.440.41
63.302.091.541.551.571.561.541.55
959.343.794.715.082.912.902.301.95
1285.3742.4645.2745.6148.8248.5848.0548.45
15124.7779.8875.7375.5175.7275.5075.0574.93
1888.73115.68110.14108.82113.04112.65112.85113.28
2148.9887.4293.2490.9389.9689.6990.2289.81
2425.2357.0560.2958.0356.8856.5556.7056.21
2717.0925.1934.5432.9333.2632.8532.6832.36
3012.4319.6421.7120.7022.2921.8821.6621.54
3319.0315.3915.1914.6315.6915.3515.0914.96
3616.3920.7417.8917.5119.3119.0218.8818.94
3916.4415.6916.8016.4317.0116.8116.7416.68
428.9712.4516.2815.8916.7116.5416.5516.56
456.8910.1710.9010.5210.4810.3610.3610.26
486.769.507.987.717.977.867.807.76
514.905.937.036.837.287.197.157.15
544.554.985.475.315.495.415.395.36
573.514.504.764.634.884.824.804.79
602.454.053.823.703.853.803.803.88
SSQ
(m3/s)2
--221.92180.41171.45161.99159.84157.64
Squared root of SSQ (m3/s)--2.822.542.472.412.392.37
NSE--0.9907920.9925140.9928860.9932790.9933680.993459
RMSE
(m3/s)
--2.8152922.5383422.4745222.4052712.3892822.372746
Table 11. Time required by Muskingum flood routing models when using Wilson’s flood data.
Table 11. Time required by Muskingum flood routing models when using Wilson’s flood data.
Comparative IndicatorsLMMLMM-LNLMMNLMM-LANLMM-LThis Study
Time (s)634633741753810936
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, E.H. Development of a New 8-Parameter Muskingum Flood Routing Model with Modified Inflows. Water 2021, 13, 3170. https://doi.org/10.3390/w13223170

AMA Style

Lee EH. Development of a New 8-Parameter Muskingum Flood Routing Model with Modified Inflows. Water. 2021; 13(22):3170. https://doi.org/10.3390/w13223170

Chicago/Turabian Style

Lee, Eui Hoon. 2021. "Development of a New 8-Parameter Muskingum Flood Routing Model with Modified Inflows" Water 13, no. 22: 3170. https://doi.org/10.3390/w13223170

APA Style

Lee, E. H. (2021). Development of a New 8-Parameter Muskingum Flood Routing Model with Modified Inflows. Water, 13(22), 3170. https://doi.org/10.3390/w13223170

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop