Next Article in Journal
Exceeding Turbidity versus Karst Spring Discharge during Single Rainfall Events: The Case of the Jadro Spring
Previous Article in Journal
Effect of Soil Texture on Water and Salt Transport in Freeze—Thaw Soil in the Shallow Groundwater Area
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Criteria Decision Analyses for the Selection of Hydrological Flood Routing Models

by
Abdurrahman Ufuk Şahin
1,* and
Arzu Özkaya
2
1
Department of Civil Engineering, Hacettepe University, 06800 Ankara, Turkey
2
Department of Civil Engineering, Ankara Yıldırım Beyazıt University, 06010 Ankara, Turkey
*
Author to whom correspondence should be addressed.
Water 2023, 15(14), 2588; https://doi.org/10.3390/w15142588
Submission received: 12 May 2023 / Revised: 7 July 2023 / Accepted: 14 July 2023 / Published: 16 July 2023

Abstract

:
In this study, a framework to circumvent the difficulties in selecting a proper flood routing method was established by employing two different multi-criteria decision analysis (MCDA) tools, namely, TOPSIS and PROMETHEE, with definite decisive criteria such as the error metrics, the number of model parameters, and the model background, under three scenarios. For eight distinct flood datasets, the parameters of 10 different Muskingum models were determined using the water cycle optimization algorithm (WCOA) and the performance of each model was ranked by both MCDA tools considering the hydrograph types of flood datasets, labeled as smooth single peak, non-smooth single peak, multi-peak, and irregular. The results indicate that both tools were compatible by giving similar model results in the rankings of almost all scenarios that include different weights in the criteria. The ranking results from both tools also showed that the routing application in single-peak hydrographs was examined better with empirical models that have a high number of parameters; however, complex hydrographs that have more than one peak with irregular limps can be assessed better using the physical-based routing model that has fewer parameters. The proposed approach serves as an extensive analysis in finding a good agreement between measured and routed hydrographs for flood modelers about the estimation capabilities of commonly used Muskingum models considering the importance of correlation, model complexity, and hydrograph characteristics.

1. Introduction

Since flood routing can help to mitigate the possible casualties and infrastructure damage caused by flood events, it is a crucial concern in the fields of hydrology and water resources engineering. Flood routing can be a complex and challenging task, as it involves many variables such as the amount of precipitation, the rate of runoff, the topography of the landscape, and the capacity of the drainage system [1]. In addition to these variables, the selection of proper methods and flow parameters increases the complexity behind the flood routing process. There are a variety of methods that can be used for flood routing, including hydraulic modeling and hydrological modeling with statistical and data-driven models [2]. In these models, the routing flow problem can be solved by considering either distributed or lumped system approaches. Using the distributed routing approach, the continuity and momentum equations, such as the Saint-Venant equation and kinematic or diffusive wave models for unsteady flow, can be used to portray the flow behavior or the movement of the floodwater [3]. Alternatively, lumped routing approaches can be constructed via the continuity equation together with a storage relation expressed as a function of time. The accuracy and the simplicity in use are trade-offs between the aforementioned approaches. While distributed-based approaches are more accurate, lumped-based solutions are simpler to use. However, due to their simplicity, lumped-based routing methods are viable for preliminary routing designs [4].
Being one of the earliest lumped-based routing models, the Muskingum model, proposed by [5], has been used in practice frequently. The linear storage relation is typically formulated as
S = K w I + 1 w Q
where S is the storage of the control section, I is the inflow and Q is the outflow, K is the storage-time constant, and w is a weighting factor. Over the years, several non-linear Muskingum models have been devised [2,6,7,8,9,10,11]. The extensive literature on the identification of the featured parameters for the models is summarized in Table 1. The classical gradient-based solvers, including non-linear least square (e.g., [12]), Lagrange multiplier (e.g., [13]), the Broyden–Fletcher–Goldfarb–Shanno method (e.g., [14]), and derivative-free optimization methods such as the Nelder–Mead simplex (e.g., [15]) are readily available methods in order to estimate the parameters of non-linear Muskingum models. In addition to these techniques, genetic algorithms (e.g., [16,17]), harmony search (e.g., [18]), particle swarm optimization (e.g., [19]), differential evolution (e.g., [20]), and immune clonal selection (e.g., [21]) are some examples of widely used population-based global optimization algorithms in flood routing. Similarly, the hybridization of population-based algorithms (e.g., [22,23,24,25,26]) is a frequently visited strategy to improve the quality of the parameter estimation procedure.
The aforementioned methods are strong and powerful, there are, however, some issues associated with flood routing procedures as comprehensively discussed in [27]. Some of the storage relations in Table 1 may be grounded on the physical parameters that involve flow characteristics, whereas some include empirical coefficients employed for improving the fit relation between inflow and outflow data. For instance, when m = 3 / 5 , the M4 form could be obtained from kinematic wave theory for a rectangular channel, as shown by [3]. However, the remaining forms are dimensionally inconsistent, which eventually yields different units of K. Dimensionally consistent models may ensure the attaining of global optima [28]. The routing approach is another important point that may directly affect the estimation performance of the implemented algorithms. Some numerical remedies for the time-marching solution could provide better results [27,28,29,30,31].
As thoroughly summarized, the capability of present models to explain a flood event may be improved by adding either physical or empirical parameters. Comparison assessments for the studied algorithms or the proposed non-linear models by researchers were generally undertaken with the very well-known benchmark dataset provided by [32]. It can be understood from the visited literature that selecting a particular storage model is as significant as the algorithm used in the parameter identification process. However, the decision making about which model works better (at least with an acceptable error level) for any kind of flood data seems to be unanswered yet. Moreover, the results given for the implemented flood data with a specific optimization algorithm may not be adequate to conduct a fair comparison. In addition to the error metrics, which are used as the performance indicators, the model complexity (number of model parameters) and the flood data characteristics (hydrographs with single peak, multi-peak, or irregular, etc.) could allow the researchers to perform an unbiased comparison assessment. This study is intended to fill the mentioned research gap via explaining the performance of different storage models with different number of parameters by employing various flood datasets. This goal was achieved by employing two well-known multi-criteria decision analysis (MCDA) tools. To the best of the authors’ knowledge, such an attempt at evaluating hydrological flood routing models with MCDA tools is one of the first studies in the literature.

2. Data and Methods

In this section, we first explain the routing approach followed by the flood routing optimization and flood datasets. The parameters of the Muskingum models implemented in this study are obtained by using the water cycle optimization algorithm (WCOA) whose details are given in the following subsections. Then, the datasets collected from previous studies are introduced. Finally, multi-criteria decision analysis tools are discussed in detail. Figure 1 summarizes the entire work flow behind this study.

2.1. Routing Approach

The mass balance equation is
d S t d t = 1 + γ I t Q t
where S is the storage of the control section, I is the inflow and Q is the outflow, and γ is the lateral flow coefficient. The general non-linear Muskingum model can be written as
S = K w 1 + γ I α + 1 w Q β m
in which K is the storage-time constant, w is the weighting factor, α is the inflow exponent, β is the outflow exponent, and m is the general model exponent. Equation (3) can be rearranged for estimating the outflow as
Q = 1 1 w S K 1 m w 1 w 1 + γ I α 1 β
For any particular time level of i + 1 , Equation (4) can be expressed as
Q i + 1 = 1 1 w S i + 1 K 1 m w 1 w 1 + γ I i + 1 α 1 β
for i = 1 , 2 , , N 1 , where N is the total number of observed data. As discussed in [27], the use of weighted I i and I i + 1 could give better performance than using I i + 1 only. Therefore, the I i + 1 term in Equation (5) is replaced with
I i + 1 = η I i + 1 η I i + 1
in which η is the weight coefficient of the explicit–implicit scheme. Equation (5) is then
Q i + 1 = 1 1 w S i + 1 K 1 m w 1 w 1 + γ η I i + 1 η I i + 1 α 1 β
When the finite difference formulation is applied to Equation (2), it becomes
S i + 1 S i t = 1 + γ I i 1 1 w S i K 1 m w 1 w 1 + γ I i α 1 β
where t denotes the time interval. If the right-hand side terms of Equation (8) are denoted by F i , it will be then
S i + 1 = S i + F i t
The flood routing process starts with the assumed initial conditions ( I 0 = I 1 and Q 0 = Q 1 ). Thus, the initial storage will be
S 1 = K w 1 + γ I 0 α + 1 w Q 0 β m
As a summary, the following computation steps complete the routing process:
Step 1:
Compute the initial storage ( i = 1 ) using Equation (10).
Step 2:
Compute S i + 1 using Equation (9).
Step 3:
Estimate Q i + 1 using Equation (7).
Step 4:
Go to step 2. Repeat the remaining computations for i = 2 , 3 , , N 1 data points.

2.2. Flood Routing Optimization

As given in Table 1, each implemented Muskingum model has a different number of model parameters, ranging from 2 to 6, and their values can be found by minimizing the cost function which is the discrepancy between the observed outflow data and estimated outflow values. When minimizing the cost function, negative storage and outflow values may appear during the space search for the optimum model parameters. Since negative values for the storage and outflow are physically impossible, such situations are regarded as the constraints of the optimization problem of interest. Therefore, the cost function with the penalized constraints can be written as
C = i = 1 N Q o b s , i Q i K , w , γ , α , β , m , η 2 + ϕ i + ψ i
where Q is the estimated outflow computed by Equation (7) and the penalty functions ϕ and ψ are
ϕ i = λ 1 S i + 1 i f S i + 1 < 0 0 o t h e r w i s e
and
ψ i = λ 2 Q i + 1 i f Q i + 1 < 0 0 o t h e r w i s e
where λ 1 and λ 2 are the penalty coefficients. In this study, λ 1 and λ 2 are assumed to be 0.5. Note that when the storage and the computed outflow values are negative, their penalty values are assumed to be equal in order to continue the flood routing process as outlined in the previous section.
In order to find the optimum model parameters for the cost function given by Equation (11), the space search could be achieved by the water cycle optimization algorithm (WCOA). This was first proposed by [33] for attaining global optima by mimicking the water cycle processes such as forming rivers, streaming, running to the sea, rainfall, and evaporation. The WCOA is constructed on the analogy that streams reach to rivers, rivers are running to the sea. In other words, the candidate solutions will be eliminated under some certain rules of the algorithm, and the best possible solution set, known as sea, will finally satisfy the global optima during the iterations. Suppose C = C o s t x is the cost function to be minimized with its d-dimensional variable set x = x 1 , x 2 , , x d . For a user-defined population size of N p o p , the WCOA starts by generating random candidate solution sets ( i = 1 , 2 , , N p o p ) for the optimized parameter set ( j = 1 , 2 , , d ) as
x i , j 0 = x j L B + Ω i , j × x j U B x j L B
where x j L b and x j U b are the lower and upper bounds of the optimized parameters, respectively, Ω i , j is a uniform random number generator ( Ω i , j 0 , 1 ), and superscript 0 denotes the initial solution at the beginning of the iterations. Once the initial values of the parameters are generated via Equation (14), the corresponding cost function values are evaluated by Equation (11) and sorted with respect to their minimum to maximum values. Then, the sorted positions are grouped as
X = S e a R i v e r 1 R i v e r 2 · R i v e r N r S t r e a m N s r + 1 S t r e a m N s r + 2 · · S t r e a m N s t r e a m = x 1 1 x 2 1 · · x d 1 x 1 2 x 2 2 · · x d 2 · · · · · · · · · · x 1 N p o p x 2 N p o p . . x d N p o p
with
N s r = N r + 1 N s t r e a m = N p o p N s r
where N r is the number of rivers, N s r is the total number of rivers and the sea, and finally N s t r e a m is the number of streams. Thus,
X s e a = X ( 1 , : ) X r i v e r = X ( 2 : N r , : ) X s t r e a m = X ( N s r + 1 : N s t r e a m , : )
By this notation, X a : b , : shows the entries of matrix X from the ath row to the bth row with all columns. As a next step, some of the streams reach the rivers, rivers are running to the sea. By means of this analogy, some solution sets will be stored through iterations and new solution sets will be created as follows
N S i = r o u n d C i i = 1 N s r C i × N s t r e a m
where N S i is the certain number of streams reaching the rivers or sea and C i = C o s t i C o s t N s r + 1 for i = 1 , 2 , , N s r . Note that Equation (18) creates an array with a length of N s r and i = 1 N s r N S ( i ) = N s t r e a m . The positions of the streams, rivers, and sea for each individual member of a d-dimensional variable set will be updated for the next iteration. The streams may reach to the sea directly. Therefore, the stream positions should be updated according to the position of sea as
X s t r e a m t + 1 i , : = X s t r e a m t i , : + Ω Ψ X s e a t X s t r e a m t i , :
for i = 1 , 2 , , N S 1 , where Ω denotes a uniform random generator vector with a length of d and Ψ is a controlling parameter varying Ψ 1 , 2 [34]. Similarly, the streams may run to the rivers. The update mechanism for this case is
X s t r e a m t + 1 j , : = X s t r e a m t j , : Ω Ψ X r i v e r t k , : X s t r e a m t j , :
for j = N S ( 1 ) + 1 , N S ( 1 ) + 2 , , N s t r e a m and k = 1 , 2 , , N s r 1 . Finally, the rivers may flow to the sea. The river positions will then be updated as
X r i v e r t + 1 i , : = X r i v e r t i , : Ω Ψ X s e a t X r i v e r t i , :
for i = 1 , 2 , , N s r 1 . Then, the evaporation phase, which means that some of the streams or rivers will be discarded from the candidate solution set in order to prevent stacking the local optima, triggers the precipitation phase, which forms new streams or rivers, with certain rules. For precipitation to a river, the new stream will be formed as
i f X s e a t X r i v e r t i , : < d m a x | | r a n d < 0.1 X s t r e a m n e w , t k , j = x j L B + r a n d × x j U B x j L B
for i = 1 , 2 , , N s r 1 , j = 1 , 2 , , d and k = N S ( 1 ) + 1 , N S ( 1 ) + 2 , , N s t r e a m , where r a n d is a uniform random number, . is the Euclidean norm, and d m a x is the user-specified small control number which provides searching intensity around the comparison norm. Similarly, when precipitation comes to the streams, a new stream will be generated as
i f X s e a t X s t r e a m t i , : < d m a x | | r a n d < 0.1 X s t r e a m n e w , t i , j = x j L B + r a n d × x j U B x j L B
for i = 1 , 2 , , N S ( 1 ) . Through the iterations, d m a x can be adjusted as
d m a x t + 1 = d m a x t d m a x t i t m a x
where i t m a x shows the maximum iteration number [35]. In this study, the WCOA is utilized with N p o p = 200 and i t m a x = 1000 .

2.3. Flood Datasets

The outlined WCOA is applied to several flood datasets for identifying the model parameters of each Muskingum model. Each dataset is categorized in terms of its characteristic shape such as smooth single peak (SSP), non-smooth single peak (NSP), multi-peak (MP), and irregular as summarized in Table 2 and shown in Figure 2.
As an SSP-type dataset, Wilson’s 1974 flood [32] is labeled as DS1 in this particular study, which is a classical benchmark dataset previously used in many studies. Brutsaert’s data (DS2) as SSP was tested by [27]. River Wye flood data (DS3 and DS7) have again been examined in several studies (e.g., [29,30,31]) and were collected from the River Wye in England. DS4, provided by [36], is associated with the flood event that occurred in the Karun river in the southwest part of Iran. Viessman and Lewis’s data is a classic example of a two-peak hydrograph, which was used to demonstrate the performance of the proposed algorithms and Muskingum models (e.g., [2,31]). Sütçüler flood (DS6), with two peaks, was gathered from Turkey [11,26]. Finally, DS8, which seems to not have a classical hydrograph shape with an apparent peak, was taken from the south canal of Chenggou and Lingqing river in China. It was also studied by [27,37].
Table 2. Flood datasets, their labels, sources, and hydrograph categories in the study.
Table 2. Flood datasets, their labels, sources, and hydrograph categories in the study.
Flood NameLabelSourceHydrograph Category
Wilson’s 1974 floodDS1 [32]Smooth Single Peak
Brutsaert’s dataDS2 [38]Smooth Single Peak
River Wye 1960 floodDS3in [27]Non-smooth Single Peak
Karun river floodDS4 [36]Non-smooth Single Peak
Viessman and Lewis’s dataDS5 [39]Multi-Peak
Sütçüler floodDS6 [11]Multi-Peak
River Wye 1982 floodDS7in [11]Irregular
Chenggou and Lingqing riverDS8 [37]Irregular

2.4. Comparison of Muskingum Models

The Muskingum models in Table 1 are ranked with certain criteria such as the error metrics with normalized root mean squared error (N-RMSE) and the adjusted coefficient of determination ( A d R 2 ). The model complexity, which is characterized by the number of optimized parameters for each model (NP), and the model background (MB) that shows whether the model parameters are physical or empirical. The N-RMSE and A d R 2 , which are the performance indicators of the WCOA-assisted parameter estimation for each dataset, can be, respectively, calculated by
N R M S E = i = 1 N y i s i 2 y m a x y m i n
and
A d R 2 = 1 1 R 2 N 1 N N P 1
with
R 2 = 1 i = 1 N y i s i 2 i = 1 N y i y ¯ 2
where N is the sample size, N P is the number of optimized parameters used in the model, y i is the observed ith data value, and s i is the calculated ith data value. y m a x , y m i n , and y ¯ are the maximum, minimum, and mean of the observed data values, respectively. As expected, the model with the smallest N-RMSE will be superior. Furthermore, higher A d R 2 values point out the powerful estimation performance. A lower NP is preferred. Finally, the models with physical parameters could be better to explain a flood event by providing the link between theory and practice.
To conduct a fair and unbiased comparison assessment, multi-criteria decision analysis (MCDA) tools including the technique for order of preference by similarity to ideal solution (TOPSIS) [40,41,42] and the preference ranking organization method for enrichment evaluation (PROMETHEE) [43] are employed to rank the alternatives (they are Muskingum models in Table 1) under definite criteria (N-RMSE, A d R 2 , NP, and MB).
Suppose A = a 1 , a 2 , , a p is a p-dimensional alternative set, C = c 1 , c 2 , , c r is an r-dimensional criterion set, and Z is the decision matrix as follows
Z = c 1 a 1 c 2 a 1 · · c r a 1 c 1 a 2 c 2 a 2 · · c r a 2 · · · · · · · · · · c 1 a p c 2 a p · · c r a p = z 11 z 12 · · z 1 r z 21 z 22 · · z 2 r · · · · · · · · · · z p 1 z p 2 · · z p r p × r
Once the decision matrix is constructed as shown in Equation (28), the TOPSIS method consists of the following steps:
Step 1:
The entries of decision matrix Z are normalized as
h i j = z i j i = 1 p z i j 2 , i = 1 , 2 , , p , j = 1 , 2 , , r
Step 2:
Each criterion has a definite weight such that i = 1 r w i = 1 . The weighted normalized matrix is then computed as
v i j = w j × h i j , i = 1 , 2 , , p , j = 1 , 2 , , r
Step 3:
As stated earlier, some criteria can be maximized (a higher value is preferred) or minimized (a lower value is preferred). The ideal best V + and the ideal worst V solutions are, respectively, determined as
V j + = m a x ( v i j ) ,   i f   c j   is   maximized m i n ( v i j ) ,   i f   c j   is   minimized   i = 1 , 2 , , p   j = 1 , 2 , , r
V j = m i n ( v i j ) ,   i f   c j   is   maximized m a x ( v i j ) ,   i f   c j   is   minimized   i = 1 , 2 , , p   j = 1 , 2 , , r
Step 4:
The Euclidean distance for the ideal best and the ideal worst solutions will be
ζ i + = i = 1 r v i j V j + 2   i = 1 , 2 , , p   j = 1 , 2 , , r
and
ζ i = i = 1 r v i j V j 2   i = 1 , 2 , , p   j = 1 , 2 , , r
Step 5:
The closeness to the ideal solution for each alternative is then computed by
χ i = ζ i ζ i + ζ i +   i = 1 , 2 , , p
Step 6:
Finally, χ i values are sorted. A larger χ value implies a better alternative.
For PROMETHEE, the computational process is as follows:
Step 1:
Step 1: The decision matrix Z in Equation (28) is normalized first according to the criteria whether some of them are maximized or minimized as
z i j = c j ( a i ) m i n ( c j ) m a x ( c j ) m i n ( c j ) ,   i f   c j   is   maximized m a x ( c j ) c i ( a i ) m a x ( c j ) m i n ( c j ) ,   i f   c j   is   minimized   i = 1 , 2 , , p   j = 1 , 2 , , r
Step 2:
The pairwise distance is
d k ( a i , a j ) = c k ( a i ) c k ( a j )   i j ,   i , j = 1 , 2 , , p   k = 1 , 2 , , r
where d k is the distance between the normalized ith alternative and jth alternative for the kth criterion. Therefore, a total of r × ( p 1 ) pairwise distances is computed.
Step 3:
The preference degree is then determined by a preference function as
π ( a i , a j ) = k = 1 r P k d k a i , a j w k
where P is the preference function. In this study, the binary preference function is used as
P d k = 0 ,   i f   d k 0 1 ,   i f   d k > 0
Thus, π values can be expressed by a p × p matrix whose diagonal entries are empty as
Γ = π ( a 1 , a 2 ) π ( a 1 , a 3 ) · · π ( a 1 , a p ) π ( a 2 , a 1 ) π ( a 2 , a 3 ) · · π ( a 2 , a p ) · · · · · · π ( a p , a 1 ) π ( a p , a 2 ) · · π ( a p 1 , a p ) p × p
Step 4:
The ideal best (positive flow) and the ideal worst (negative flow) solutions are
ξ i + a i = 1 p 1 i = 1 p π a i , x
and
ξ i a i = 1 p 1 i = 1 p π x , a i
In other words, ξ i + is the mean of the ith row of matrix Γ , whereas ξ i is the mean of the ith column of matrix Γ .
Step 5:
The ideal solution (net flow) is
ξ i ( a i ) = ξ i + ( a i ) ξ i ( a i )
Step 6:
The PROMETHEE provides the complete ranking by ordering alternatives according to the decreasing values of the ideal solution ξ .
As a summary, the comparison assessment of Muskingum models is conducted via the following scheme:
1.
The parameters of each model in Table 1 are estimated by the WCOA for each dataset given in Table 2.
2.
The computed outflow values obtained from (1) are substituted into Equations (25) and (26) in order to obtain the N-RMSE and A d R 2 values of each model for the particular dataset.
3.
For each dataset, the estimation score of each model is noted and sorted. Since there is a total of 10 models, the model with the smallest N-RMSE is given a rank of 1 while the highest one is given a rank of 10. A smaller rank indicates a better performance.
4.
Based on the N-RMSE values, the ranks for each dataset are found and then averaged.
5.
Similarly, the A d R 2 values of each model for each dataset are sorted. The model with the highest A d R 2 value is given a rank of 1 while the smallest one is given a rank of 10. A higher rank score shows a better performance.
6.
With A d R 2 values, the rank scores are averaged.
7.
The decision matrix is formed by NP and MB together with the rank scores of N-RMSE and A d R 2 . A binary code (1 for the models with physical parameters, 0 for the models with empirical parameters) is used to classify the MB.
8.
The obtained decision matrix acquired from the previous steps is then studied by TOPSIS and PROMETHEE to elucidate the most efficient Muskingum model.

3. Results and Discussion

3.1. Estimation Performance of Models

Following the WCOA-based parameter estimation approach detailed in the previous section, the model parameters of the Muskingum models listed in Table 1 were obtained using the datasets of eight distinct flood events. For each flood event, the parameters and prediction performance of the models used are presented in Table A1. During the optimization process, the WCOA was run 10 times for each flood dataset in order to obtain more reliable model parameters with accurate ranges. For instance, when the model parameters of M1 were determined by using dataset 1 (DS1), 10 possible candidate solutions were retrieved by the WCOA. Although the estimated model parameters could be very close to each other for each trial, the estimated parameter set with the smallest cost function value among the candidate solutions was noted. Thus, the error metrics for each model were computed, as given in Table A1. The same process was repeated to obtain the error metrics for the remaining models.
According to the findings presented in Table A1, models M7, M8, M9, and M10 yielded highly accurate predictions considering the error metrics of A d R 2 (>0.999) and N-RMSE (<0.01) for DS1 as the smooth single peak (SSP) dataset. For the dataset DS2, which shares similar attributes to DS1, it is worth noting that, in addition to the models mentioned previously, models M2, M3, and M5 also exhibit considerable predictive performance with respect to the aforementioned metrics of A d R 2 (>0.999) and N-RMSE (<0.01). For the non-smooth single peak (NSP) datasets (DS3 and DS4), when evaluated with both metrics ( A d R 2 > 0.98 and N-RMSE < 0.31), models M9 and M10 provide the best results, but it is observed that the prediction capacity of these models decreases when compared to those in SSP datasets. When examining the DS5 dataset, which is of the multi-peak (MP) type, it was observed that all models, except for the M4 model, provided relatively good prediction results. Another dataset (DS6) of the same type (MP) showed that the models between M5 and M10 performed well. However, in these datasets (DS5 and DS6), the performance of the M5, M9, and M10 models stands out compared to others when considering both error metrics. Finally, for the datasets referred to as irregular (DS7 and DS8) in this study, the M5, M9, and M10 models were consistently ranked as the top three performers among the 10 models evaluated in both A d R 2 ( A d R 2 > 0.99) and N-RMSE (N-RMSE < 0.015 for DS7 and N-RMSE < 0.025 for DS8). This situation was also observed for the estimation performance of the models realized in the MP type.
Figure 3 also depicts the relative percentage error of each model over the given flood duration for each dataset. For the SSP-type datasets (DS1 and DS2), all methods could achieve better parameter estimation performance within a certain error tolerance. The falling and the rising limbs of each hydrograph for M4 show an increase in the error percentages and this error is even more pronounced in the falling limbs of the hydrographs, especially for the SSP type. In addition, a similar pattern can be observed for the NSP-type datasets (DS3 and DS4). Likewise, the relative errors of the models, except for model M4, were very small for DS5, that belongs to the MP-type dataset category. However, all methods, in general, gave reasonable estimation results for DS6, that was categorized as an MP-type dataset. Models M1, M3, and M4 were not able to produce the model parameters satisfying a good performance for DS7 while the remaining ones successfully estimated the model parameters. For DS8 as an irregular dataset type, the model parameters of all methods could be easily identified by the WCOA.
The Wilson river flood dataset (DS1) has been studied by several researchers. The reported SSE values (in (m3/s)2) for this dataset were obtained as 39.8 by [44], 17.55 by [45], 9.82 by [11], 7.67 by [2], 5.124 by [46], 4.11 by [47], 4.04 by [25], 1.92 by [31], 1.092 by [48], 0.799 by [49], and 0.65 by [50]. In this study, the SSE value for DS1 with the WCOA was found to be 7.66. Moreover, the proposed model (M10) achieved an SSE value of 4.09. The range of SSE values reported by different researchers suggests that there may be considerable variability in flood modeling results for this dataset, depending on the specific optimization algorithm and routing approach used. For instance, the SSE value of model M7 in this study (with the four-parameter non-linear model proposed by [2]) was calculated as 7.67, while the same metric was found to be 9.82 for model M9 (with the four-parameter non-linear model by [11]). Although these two models have the same number of model parameters, the reported SSE values are different because the researchers used different routing approaches and optimization algorithms as well. Similar findings can also be drawn for other datasets used in this study. For the Viessman and Lewis flood dataset (DS5), the reported SSE values were 71,708 in [45], 65,324 in [44], 28,855 in [50], and 8449 in [49]. However, the implemented routing model with the WCOA in this study provided notably smaller SSE values using the examined models except for M4, as shown in Table A1. It is quite remarkable that the SSE value of model M1 (with the two-parameter linear model proposed by [5]) was 14.53, while a model proposed by [49] with 12 parameters yielded an SSE value of 8449. However, comparing the reported results in the literature directly (without taking into account the applied routing approaches and optimization methods) can be misleading in interpreting which model exhibited a better estimation performance.
In summary, this exercise reveals that the WCOA is able to show a very competitive estimation performance when compared to the reported literature. Thus, the WCOA can be regarded as a viable algorithm with the outlined routing approach. These findings also suggest that the particular models are effective in capturing the outflow patterns in these datasets.
To interpret the performance of the models in Table A1, the models were ordered according to their A d R 2 and N-RMSE. Considering A d R 2 values, the rank score was then assigned to each model such that ranking starts from 1 for the model having the highest A d R 2 value and successively continues to the rank value of 10. For N-RMSE values, the rank score starts from 1 for the smallest N-RMSE value, as shown in Table 3. Furthermore, the averaged rank (AR) of each model was given for the datasets previously categorized as SSP, NSP, MP, and irregular. Based on the overall AR value of A d R 2 and N-RMSE in Table 3, the best two models were found to be M9 and M10, whereas the worst two models were noted to be M1 and M4. When the results are analyzed according to the hydrograph type, model M4 shows the worst performance with a rate of 75% (six times it was assigned with a rank score of 10 among the eight datasets) for both criteria. Although this model has a physical sense evolving from the Manning equation, forcing the model with a constant exponent (m = 3/5) decreases the performance in the calculations. This remark is also supported by the relative error results, presented in Figure 3. However, the best model is not the same for all hydrograph types. Although half of the datasets, more specifically four datasets for the A d R 2 criterion and five datasets for the N-RMSE criterion, have superior results with model M10; models M8, M5, and M9 are ranked as the best for the hydrograph types of SSP, MP, and irregular, respectively.
Considering Table A1 and Table 3, any comparison between the error metrics could be insufficient to form a solid judgment. Therefore, a comprehensive assessment was needed to compare the capabilities of the models. To apply the methodologies of MCDA outlined previously, the decision matrix for each hydrograph type was set as listed in Table 4. The results presented in Table 5 cover the ranking obtained from the MCDA tools considering four hydrograph types under three scenarios. To this end, the average rank values for A d R 2 and N-RMSE were employed as decisive criteria for each dataset. The performance of a model can be evaluated using two metrics, and in this context, a lower value assigned implies better performance. Specifically, the metrics are inversely proportional to the model’s effectiveness, meaning that a decrease in the assigned value indicates an improvement in the model’s performance. Therefore, the aforementioned two metrics were considered as non-beneficiary (minimized) criteria. The number of parameters (NP) for a model could be again considered as a decisive criterion. As listed in Table 1, the NP could indicate the models’ complexity. A model containing fewer parameters may be more desirable than a model with a larger number of parameters. Thus, NP was a non-beneficiary (minimized) criterion. Finally, the dimensional consistency of a model was labeled as the model background (MB) criterion. A model constructed using physical-based parameters was deemed more realistic compared to a model that employed additional parameters solely for improving the fitting accuracy. In this criterion, a binary entry of 1 was assigned to physically sound models, while models containing non-physical parameters solely used for fitting purposes were assigned a binary value of 0, as shown in Table 4. Therefore, the MB was assigned as a beneficiary (maximized) criterion.

3.2. Multi-Criteria Decision Analysis (MCDA) Results

To apply MCDA tools effectively, it is essential to identify the weight of each criterion. To achieve this, three distinct weight scenarios were designed. In scenario 1, the weights assigned to the A d R 2 , N-RMSE, NP, and MB criteria were 0.275, 0.275, 0.225, and 0.225, respectively. This scenario implies that the weights for the error metrics criteria (55%) are slightly higher than those for the remaining criteria (45%). Therefore, the results from scenario 1 were expected to provide a balanced ranking of the employed models taking into account the implemented criteria. In scenario 2, weights of 0.3, 0.3, 0.2, and 0.2 were assigned to the A d R 2 , N-RMSE, NP, and MB criteria, respectively. The objective of scenario 2 was to observe the impact of increasing the weights for the error metrics on the ranking. Finally, scenario 3 was designed to assign more weight to the error metric criteria (0.35 each), while the remaining NP and MB criteria were given weights of 0.15 each. Through this assessment, these scenarios can provide a better understanding of how different the weight distribution to each criterion may affect the ranking of the models implemented.
Both TOPSIS and PROMETHEE identified model M8 (a non-linear model with five parameters proposed by [10]) as the most suitable model for the hydrograph type of SSP in all scenarios. With the NSP hydrograph type, model M5 (a linear model with three physical-based parameters) was found to be the best model according to the aforementioned criteria for scenario 1, as determined by both tools. However, in scenario 3, where the error performance was dominant, model M10 (non-linear model with six parameters) was identified as the best alternative model among the comparison poll. It should be noted that model M10 was identified as the best performing model based on the A d R 2 and N-RMSE criteria as shown in Table 4. However, the ranking order may shift due to the influence of the NP and MB criteria, which can cause model M5 to move from a lower position to the top, as observed in scenario 1. This underscores the importance of using a multi-criteria decision analysis approach to ensure a fair and comprehensive evaluation of the capabilities of the examined models. For the MP and irregular hydrograph types, the application of both MCDA tools consistently identified the M5 model and its non-linear variant M9 (with four parameters) as the best and second-best models, respectively, across all scenarios. However, a highly competitive model, M10, was observed to have superior performance based on its A d R 2 and N-RMSE values as listed in Table 4 for the MP and irregular types. Notably, the ranking of the models was observed to be sensitive to the weight scenarios utilized in this study, as demonstrated by the varying influence of the NP and MB criteria. As observed from Table 5, M1, the classical linear model with two parameters, showed a notable performance when model simplicity and background were preferable, as aimed for in scenario 1. Both MCDA tools determined the fourth ranking place for the SSP, and third place for the MP and irregular types. However, fourth place can be obtained by TOPSIS, whereas PROMETHEE put it in sixth place when the NSP-type datasets are employed. All in all, model M1 is still a viable model if the complexity as a decisive criterion is taken into account by the practitioners. Furthermore, the data characteristic may affect the efficiency of the models. Model M8 was found to be the best model for the SSP, whereas its capability dramatically varied for the MP and irregular hydrograph types under the scenarios implemented. As an additional evaluation, a new decision matrix was created using the overall average AR scores of eight flood events (shown in the last column of Table 3) to eliminate any potential impact of hydrograph type, as given in Table 6. Once again, the M5 model and its non-linear variant M9 were found to be the superior models for scenarios 1 and 2, as indicated in the decision matrix presented in Table 3. In scenario 3, where error metrics were dominant, TOPSIS ranked model M9 in the third place, while PROMETHEE placed it in second place.
As mentioned earlier, Table 5 presents the results of four different hydrograph types, while Table 6 displays the outcomes based on the average AR values of eight distinctive flood events. To facilitate the interpretation of the results for each model, the ranking places of each model can be categorized into five main groups: the first or the second place (denoted as 1||2) in group 1 (G1), the third or the fourth place (3||4) in group 2 (G2), the fifth or the sixth place (5||6) in group 3 (G3), the seventh or the eighth place (7||8) in group 4 (G4), and the ninth or the tenth place (9||10) in group 5 (G5), as illustrated in Figure 4. In other words, there are 15 possible cases (four for hydrograph types and one aggregated evaluation for all datasets by three scenarios) provided by each MCDA tool for each model. For instance, in the case of SSP using TOPSIS, the ranking places of the M1 model were found to be the fourth place in scenario 1, sixth place in scenario 2, and ninth place in scenario 3. This procedure was carried out for each model and the implemented cases. The ranking places were noted for each model among the implemented cases and subsequently categorized as explained. For example, the performance of model M1 was categorized as 0% in G1, 46.6% (7 occurrences out of 15) in G2, 20% (3 out of 15) in G3, 26.6% (4 out of 15) in G4, and 6.6% (1 out of 15) in G5, as shown in Figure 4. As also seen in Figure 4, model M5 was able to produce a superior performance to complete the comparison within the first and the second ranking places (G1) with scores of 73.3% by TOPSIS and 66.6% by PROMETHEE. Similarly, model M9 can be evaluated as an identical competitor model regarding its scores in G1, which were 66.6% by TOPSIS and 80% by PROMETHEE. Apart from these models, model M10 can give the highest success rate for the third and fourth places (G2), which were 66.6% by TOPSIS and 80% by PROMETHEE. For the G2 category, model M1 showed a notable performance, with equal success scores of 46.6% by TOPSIS and PROMETHEE. For an alternative interpretation of Figure 4, the success rate for each ranking place category is cumulatively shown in Figure 5 in order to understand which model stands out. As a result of the analyses conducted using both MCDA tools, models M5, M9, and M10 emerge as the top performers among the flood routing models in the comparison pool, as illustrated in Figure 5.

3.3. Comparison of MCDA Tools

In this study, two outranking MCDA tools, TOPSIS and PROMETHEE, were employed to evaluate the performance of different Muskingum models under the aforementioned criteria. TOPSIS and PROMETHEE are widely used in decision-making processes due to their ability to handle complex decision problems [51]. In general, any MCDA method can develop a different series of rankings for the same problem depending on the strategy and mathematical background of MCDA implemented. TOPSIS is sensitive to the weights assigned to the criteria and ignores the interrelationships between the criteria [52]. PROMETHEE, on the other hand, generates a pairwise comparison matrix to assess the relative importance of each criterion [53]. However, it is sensitive to the choice of preference function and decision threshold used [54]. Therefore, selecting the appropriate MCDA method depends on the specific problem being addressed and the characteristics of the data being evaluated [55,56,57].
To analyze the differences in rankings obtained from TOPSIS and PROMETHEE, a consistency rate (CR), which is the percentage of cumulative differences of the model ranked from the same case for both MCDA tools, was defined. As observed from Table 5, both MCDA tools ranked models M8, M1, M10, M2, and M3 in the 1st, 4th, 5th, 9th, and 10th ranking places, respectively, in scenario 1. Thus, the CR of both tools for the SSP was realized as 50% (5 out of 10) in scenario 1, 100% in scenario 2 and 70% (7 out of 10) in scenario 3, which yields an average of 73.3% overall. For the NSP type, the CR was observed to be 40%, 30%, and 50% for scenarios 1, 2, and 3, respectively. It is noteworthy that the CR is 70% in scenario 1, 60% in scenario 2, and 60% in scenario 3 for the MP type, whereas the CR values of the irregular type were 80%, 50%, and 30% in scenarios 1, 2, and 3, respectively. This analysis is denoted CR(0) and depicted in Figure 6. The plus (+) values in Figure 6 indicate that TOPSIS ranks the model higher compared to PROMETHEE, while the minus (−) values indicate the opposite condition. Furthermore, the CR was recomputed in such a way that ± 1 rank difference between the MCDA tools was tolerated and denoted as CR(1), shown in a gray color in Figure 5. Considering the CR(1) values, both MCDA tools are in a high accordance for almost all cases implemented in this study. Thus, this analysis highlights the effectiveness of these MCDA tools in providing reliable model rankings.

4. Summary and Conclusions

In this research, an extensive analysis was conducted on the estimation capabilities of several commonly used Muskingum models in flood routing using two different MCDA methods, namely, TOPSIS and PROMETHEE, with the help of eight distinct flood datasets reported in the literature. The estimation performance of the flood routing models may depend on diverse factors such as the model parameters, which can be either physical-based or empirical ones that minimize the error between the output flow calculated by the flood routing methods and the observed output flow, the effectiveness of the applied optimization method, and the numerical flood routing procedure with implicit or explicit schemes. The results demonstrated that not every proposed model performed well on every dataset. Apart from utilizing error metrics to assess the estimation performance of the models, it is also crucial to thoroughly examine how each model performs with varying degrees of success on different datasets. The key findings can be summarized as follows:
(1)
The WCOA used in this study showed a viable ability in parameter estimation when compared to other algorithms used in similar studies.
(2)
The criteria, A d R 2 and N-RMSE, revealed that the non-linear models M9 and M10 (with four and six parameters) showed better predictive performance than the rest of the models on all datasets. However, the temporal distribution of the relative error showed diversity in defining the best routing models overall. Such findings indicate the importance of the criteria in the selection of the proper model and parameters.
(3)
The model performances were also assessed by clustering the hydrograph types, smooth single peak (SSP), non-smooth single peak (NSP), multi-peak (MP), and irregular. According to the average of the ranked criteria overall, models M9 and M10 showed the best performance according to the error metric criteria, including A d R 2 and N-RMSE. For hydrograph types, the definition of the best model showed some variability, the proposed model M10 showed better performance for the majority of the datasets by considering these criteria.
(4)
To obtain a more rigid conclusion about defining the best model, two different multi-criteria decision analysis (MCDA) tools were applied by using two additional criteria: the number of model parameters (NP), which indicates the model complexity; and the model background (MB), which characterizes the model parameters as either physical-based or empirical, with the purpose of improving the fitting.
(5)
The effects of the relative importance of each criterion were investigated by assigning three different weight scenarios in both tools, TOPSIS and PROMETHEE. The most outstanding remark is the compatibility of both tools in all scenarios, giving similar model results in the rankings. However, it was seen that the governing factor in defining the best model depends on the scenario selection. In scenario 1, where additional criterion weights are balanced when compared to the other scenarios, physical-based models were more successful. However, there is a variety of ways to define the best model for hydrograph types. Model M5, a physically sound model, showed the best results for MP and irregular hydrograph types, whereas models M8 and M10, empirical models with a high number of parameters, showed better results compared to the rest of the models for the SSP and NSP types. Such outcomes may indicate that the routing application in complex hydrographs that have more than one peak with irregular limps can be assessed better using physical-based routing models that have fewer parameters.
(6)
The MCDA in a single pool ignoring the clustering of hydrograph types showed that model M5 had the greatest score in the rankings. Moreover, the cumulative score evaluation of 10 Muskingum models that cover the cumulative percentages of the model ranking place shows that models M5, M9, and M10 have the leading results with scores of 85% in the top four ranking places for both tools, TOPSIS and PROMETHEE. Such a result can be beneficial for researchers in choosing the proper routing models in their area of interest.
(7)
There are a variety of Muskingum models applied in the literature. Within this study, the success of the models was assessed concerning the model parameters (complexity), model constructions (empirical or physical), error metrics, and hydrograph types. The most suitable model, model M5, showed that the use of physical-based models with fewer parameters among the possible model sets studied here, may be beneficial.
(8)
Any MCDA method can inherently generate a different ranking for the problem of interest. Although the tools used in this study have different capabilities in terms of relating the criteria and being sensitive to the weights of the criteria, their results are compatible with each other. This outcome reveals the independence of the criteria and stability of the weights in application. This finding was also testified by the consistency rates (CR) that involve the proximity of model rankings for each scenario and hydrograph types performed for both tools.
In this study, although a comprehensive assessment was completed by involving the factors aforementioned above, there is a limited number of datasets. Future studies may improve on the presented results by using new flow datasets. The findings in this work are the result of eight different flood case studies reported previously in the literature. Since the workflow given here does not depend on a specific region, it can be applied to other regions as well. Considering all the events, it can be generalized that achieving a successful running model depends on the purpose of the work and the hydrograph characteristics. Although the purpose of the work that covers the model complexity and the importance of correlation can be decided before a flood happens, estimating the hydrograph type is a challenging task. Such an issue, at least the prediction of the inflow hydrograph, may be handled by implementing hydrologic models that can work with numerical weather prediction products. In addition to remote sensing products with proper hydrologic models, field and laboratory experiments should be made in the future that may provide a broad generalization in the selection of the best routing model.

Author Contributions

Conceptualization, A.U.Ş.; methodology, A.U.Ş.; formal analysis, A.U.Ş.; investigation, A.U.Ş. and A.Ö.; writing—original draft preparation, A.U.Ş. and A.Ö.; writing—review and editing, A.U.Ş. and A.Ö. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All data used in this study can be provided upon by the request.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ARAveraged rank
MBModel background
MCDAMulti-criteria decision analysis
MPMulti-peak
NPNumber of parameters
N-RMSENormalized root mean squared error
NSPNon-smooth single peak
PROMETHEEThe preference ranking organization method for enrichment evaluation
TOPSISThe technique for order of preference by similarity to ideal solution
SSESum of square errors
SSPSmooth single peak
WCOAWater cycle optimization algorithm

Appendix A. Results of Parameter Estimation for Models

Table A1. Results for the parameter estimation for each model and each dataset ( α is the flow exponent, m is the general model exponent, γ is the lateral flow, β is the flow exponent, and η is the coefficient of the explicit/implicit scheme).
Table A1. Results for the parameter estimation for each model and each dataset ( α is the flow exponent, m is the general model exponent, γ is the lateral flow, β is the flow exponent, and η is the coefficient of the explicit/implicit scheme).
DatasetModelKw γ α β m η SSE Ad R 2 N-RMSEMAE
DS 1M129.81060.2388 1415.550.96240.06593.77
M22.76370.2287 1.5012 1245.580.97660.05062.76
M30.51750.2869 1.8681136.770.99650.01961.07
M4235.83400.1661 0.60001752.000.93200.08864.96
M529.74960.2477−0.0091 1411.690.96070.06553.74
M61.62470.0003 3.04181.5676 1184.320.98140.04392.45
M70.83380.2956 0.4331 4.079017.670.99920.00890.47
M80.86460.0395 0.74380.36794.35000.93074.960.99950.00720.29
M90.53420.3005−0.0216 1.86421.09.820.99900.01010.57
M100.66630.0246−0.00770.98500.49533.38660.86574.090.99950.00650.31
DS 2M12.08390.1102 115,895.600.99870.010916.73
M21.08090.1332 1.0839 0.960211,846.220.99900.009414.73
M30.81890.0388 1.1107012,144.800.99900.009514.95
M451.92610.0000 0.60000.5320155,580.390.98740.034153.38
M52.00770.02510.0129 111,913.600.99900.009415.04
M6586.12000.9992 0.23141.1572 0.95548631.990.99920.008012.78
M71.59620.1616 1.7865 0.58231.00008244.860.99930.007812.31
M80.91370.0001 5.00004.09570.26870.51037711.480.99930.007611.77
M91.12350.05600.0117 1.07461.00009043.470.99920.008212.48
M1013.65080.96750.00410.82431.48510.71800.87238510.250.99920.008012.45
DS 3M126.65730.4274 188,387.010.94310.056536.16
M21.63700.4024 1.4050 144,478.760.97040.040125.01
M30.46240.4091 1.5856135,194.510.97660.035622.65
M4491.33090.3475 0.60001200,258.510.87110.085058.54
M526.44130.39920.0465 182,057.200.94540.054435.47
M60.67730.0022 2.24971.4641 137,352.540.97430.036720.56
M70.45020.4074 1.2022 1.3229129,474.960.97970.032620.19
M80.95400.7331 1.00381.18881.3742129,132.200.97920.032420.41
M90.36910.38300.0547 1.6141125,915.270.98220.030617.09
M100.75000.04520.13230.54290.26765.0000120,048.640.98520.026914.51
DS 4M112.56060.2424 189,417.760.97350.054431.20
M210000.1846 0.4563 165,726.230.98010.046627.22
M310000.1749 0.4570164,593.540.98040.046227.07
M4297.39670.1989 0.6000168,814.880.97960.047728.12
M512.35210.3608−0.0314 163,325.140.98080.045827.55
M610000.1101 0.51590.4462 165,666.310.97960.046627.19
M710000.1146 2.8707 0.1600162,607.580.98060.045527.00
M810000.0000 5.00003.49670.1301160,243.110.98080.044626.12
M9980.26460.2431−0.0198 0.4564154,947.640.98290.042625.64
M1010000.0000−0.02495.00002.86880.1551147,663.450.98450.039722.80
DS 5M13.81580.4489 114.531.00000.00060.62
M23.83760.4490 0.9993 113.981.00000.00060.60
M33.80030.4489 1.0005114.331.00000.00060.62
M483.39250.3243 0.60001118,587.650.97250.051357.10
M53.81740.4497−0.0007 17.571.00000.00040.21
M63.75460.4300 1.00740.9978 110.561.00000.00050.42
M73.79600.4490 0.9952 1.005518.961.00000.00040.36
M83.79180.4476 0.99600.99531.005218.951.00000.00040.36
M93.79770.4496−0.0007 1.000717.241.00000.00040.22
M103.77210.4411−0.00071.00531.00100.998616.981.00000.00040.22
DS 6M11.04780.0000 0561.300.99030.02173.12
M210000.6980 5.0000 0.9836327.350.99410.01661.86
M30.98670.0000 1.01060.9995560.620.98990.02173.14
M411.47630.0000 0.60000.98951979.650.96570.04086.19
M50.51970.9999−1.0001 0.9802283.390.99490.01541.68
M610000.7460 5.00005.0535 0.9877322.010.99400.01651.92
M70.01000.6980 5.0000 4.78680.9836327.350.99390.01661.86
M81.29500.0091 0.69120.17985.00001319.370.99380.01641.94
M9126.2182−1.00000.8970 5.00000.9945292.350.99450.01571.82
M1010000.3908−2.62091.68071.70521.99280.9996280.560.99430.01541.81
DS 7M1131.57740.9917 0.956136,896.67−0.35630.375226.36
M2307.88550.7268 0.4827 0.29515272.040.79930.141810.04
M310000.9622 0.26560.454524,252.220.07660.304216.28
M4325.59260.9808 0.60000.824025,624.600.05800.312717.65
M55.15230.22932.5276 154.860.99790.01450.95
M6341.81200.9086 0.42490.6222 0.73613222.910.87270.11098.69
M731.62920.7387 0.5061 4.09490.49613296.050.86990.11217.43
M8125.80950.8845 0.46400.61212.430113244.810.86700.11137.88
M95.67650.22712.5298 0.9800153.660.99790.01430.90
M106.84490.33472.50031.19991.24970.7737144.520.99810.01300.84
DS 8M1116.34200.9914 14582.580.99030.030210.18
M20.12680.0000 1.3211 14542.590.99000.03019.80
M30.12680.0000 1.321114542.590.99000.03019.80
M419.42750.0000 0.60000.89155898.410.98750.034311.24
M50.03210.9994−1.0003 0.1399992.930.99780.01414.70
M60.12680.0000 0.49081.3211 14542.590.98950.03019.80
M70.12680.0000 0.2642 4.99990.57494542.590.98950.03019.80
M82.54890.8092 0.10000.50003.052614534.570.98910.03019.81
M90.2358−1.0000−0.0027 1.21340.1516918.240.99790.01354.51
M1010000.9450−1.18175.05.19150.16370.93903160.900.99210.02518.48

References

  1. Afzali, S.H. Variable-parameter Muskingum model. Iran. J. Sci. Technol. Trans. Civ. Eng. 2016, 40, 59–68. [Google Scholar] [CrossRef]
  2. Easa, S.M. New and Improved Four-Parameter Non-Linear Muskingum Model. Proc. Inst. Civ. Eng. Water Manag. 2014, 167, 288–298. [Google Scholar] [CrossRef]
  3. Gąsiorowski, D.; Romuald, S. Dimensionally Consistent Nonlinear Muskingum Equation. J. Hydrol. Eng. 2018, 23, 04018039. [Google Scholar] [CrossRef]
  4. Bai, T.; Wei, J.; Yang, W.; Huang, Q. Multi-objective parameter estimation of improved Muskingum model by wolf pack algorithm and its application in Upper Hanjiang River, China. Water 2018, 10, 1415. [Google Scholar] [CrossRef] [Green Version]
  5. McCarthy, G.T. The unit hydrograph and flood routing. In Conference of North Atlantic Division; US Army Corps of Engineers: Washington, DC, USA, 1938. [Google Scholar]
  6. Chow, V.T. Open-Channel Hydraulics; McGraw-Hill: New York, NY, USA, 1959. [Google Scholar]
  7. Gill, M.A. Flood routing by the Muskingum method. J. Hydrol. 1978, 36, 353–363. [Google Scholar] [CrossRef]
  8. O’donnell, T. A Direct Three-Parameter Muskingum Procedure Incorporating Lateral Inflow. Hydrol. Sci. J. 1985, 30, 479–496. [Google Scholar] [CrossRef] [Green Version]
  9. Gavilan, G.; Houck, M.H. Optimal Muskingum river routing. In Proceedings of the ASCE Computer Applications in Water Resources, Buffalo, NY, USA, 10–12 June 1985; American Society of Civil Engineers: New York, NY, USA, 1985; pp. 1294–1302. [Google Scholar]
  10. Vatankhah, A.R. Discussion of parameter estimation of the nonlinear Muskingum flood routing model using a hybrid harmony search algorithm by Halil Karahan, Gurhan Gurarslan, and Zong Woo Geem. J. Hydrol. Eng. 2014, 19, 839–842. [Google Scholar] [CrossRef]
  11. 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]
  12. Yoon, J.W.; Padmanabhan, G. Parameter estimation of linear and nonlinear Muskingum models. ASCE J. Water Resour. Plan. Manag. 1993, 119, 600–610. [Google Scholar] [CrossRef]
  13. Das, A. Parameter estimation for Muskingum models. ASCE J. Irrig. Drain. Eng. 2004, 130, 140–147. [Google Scholar] [CrossRef]
  14. Geem, Z.W. Parameter estimation for the nonlinear Muskingum model using the BFGS technique. ASCE J. Irrig. Drain. Eng. 2006, 132, 474–478. [Google Scholar] [CrossRef]
  15. Barati, R. Parameter estimation of nonlinear Muskingum models using Nelder–Mead simplex algorithm. ASCE J. Hydrol. Eng. 2011, 16, 946–954. [Google Scholar] [CrossRef]
  16. Mohan, S. Parameter estimation of nonlinear Muskingum models using genetic algorithm. ASCE J. Hydraul. Eng. 1997, 123, 137–142. [Google Scholar] [CrossRef]
  17. Zhang, S.; Ling, K.; Liwei, Z.; Xiaoming, G. A New Modified Nonlinear Muskingum Model and Its Parameter Estimation Using the Adaptive Genetic Algorithm. Hydrol. Res. 2017, 48, 17–27. [Google Scholar] [CrossRef]
  18. Kim, J.H.; Geem, Z.W.; Kim, E.S. Parameter estimation of the nonlinear Muskingum model using harmony search. J. Am. Water Resour. Assoc. 2001, 37, 1131–1138. [Google Scholar] [CrossRef]
  19. Chu, H.J.; Chang, L.C. Applying particle swarm optimization to parameter estimation of the nonlinear Muskingum model. ASCE J. Hydrol. Eng. 2009, 14, 1024–1027. [Google Scholar] [CrossRef]
  20. Xu, D.; Qiu, L.; Chen, S. Estimation of nonlinear Muskingum model parameters using differential evolution. ASCE J. Hydrol. Eng. 2015, 16, 348–353. [Google Scholar] [CrossRef]
  21. Luo, J.; Xie, J. Parameter estimation for nonlinear Muskingum model based on immune clonal selection algorithm. ASCE J. Hydrol. Eng. 2010, 15, 844–851. [Google Scholar] [CrossRef]
  22. Ehteram, M.; Binti Othman, F.; Mundher Yaseen, Z.; Abdulmohsin Afan, H.; Falah Allawi, M.; Najah Ahmed, A.; Shahid, S.; Singh, V.P.; 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]
  23. 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]
  24. Bozorg-Haddad, O.; Hamedi, F.; Orouji, H.; Pazoki, M.; Loáiciga, H.A. A re-parameterized and improved nonlinear Muskingum model for flood routing. Water Resour. Manag. 2015, 29, 3419–3440. [Google Scholar] [CrossRef]
  25. 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]
  26. Farzin, S.; Singh, V.P.; Karami, H.; Farahani, N.; Ehteram, M.; Kisi, O.; Allawi, M.F.; Mohd, N.S.; El-Shafie, A. Flood routing in river reaches using a three-parameter Muskingum model coupled with an improved bat algorithm. Water 2018, 10, 1130. [Google Scholar] [CrossRef] [Green Version]
  27. Geem, Z.W. Issues in Optimal Parameter Estimation for the Nonlinear Muskingum Flood Routing Model. Eng. Optim. 2014, 46, 328–339. [Google Scholar] [CrossRef]
  28. Gąsiorowski, D.; Romuald, S. Identification of Parameters Influencing the Accuracy of the Solution of the Nonlinear Muskingum Equation. Water Resour. Manag. 2020, 34, 3147–3164. [Google Scholar] [CrossRef]
  29. Easa, S.M. Improved nonlinear Muskingum model with variable exponent parameter. ASCE J. Hydrol. Eng. 2013, 18, 1790–1794. [Google Scholar] [CrossRef]
  30. Karahan, H. Discussion of ‘‘Improved nonlinear Muskingum model with variable exponent parameter” by Said M Easa. J. Hydrol. Eng. 2014, 19, 1–9. [Google Scholar] [CrossRef]
  31. Ayvaz, M.T.; Gurarslan, G. A new partitioning approach for nonlinear Muskingum flood routing models with lateral flow contribution. J. Hydrol. 2017, 553, 142–159. [Google Scholar] [CrossRef]
  32. Wilson, E.M. Engineering Hydrology; Macmillan Education: Basingstoke, UK, 1974. [Google Scholar]
  33. Eskandar, H.; Sadollah, A.; Bahreininejad, A.; Hamdi, M. Water cycle algorithm-A novel metaheuristic optimization method for solving constrained engineering optimization problems. Comput. Struct. 2012, 110–111, 151–166. [Google Scholar] [CrossRef]
  34. Sahin, A.U. Automatic Shifting Method for the Identification of Generalized Radial Flow Parameters by Water Cycle Optimization. Water Resour. Manag. 2021, 35, 5205–5223. [Google Scholar] [CrossRef]
  35. Nasir, M.; Sadollah, A.; Choi, Y.H.; Kim, J.H. A comprehensive review on water cycle algorithm and its applications. Neural Comput. Appl. 2020, 32, 17433–17488. [Google Scholar] [CrossRef]
  36. Orouji, H.; Bozorg, H.O.; Fallah-Mehdipour, E.; Marino, M.A. Estimation of muskingum parameter by meta-heuristic algorithms. Proc. Inst. Civ. Eng. Water Manag. 2013, 166, 315–324. [Google Scholar] [CrossRef]
  37. 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]
  38. Brutsaert, W. Hydrology: An Introduction; Cambridge University Press: New York, NY, USA, 2005. [Google Scholar]
  39. Viessman, W.; Lewis, G.L. Introduction to Hydrology; Pearson Education, Inc.: Upper Saddle River, NJ, USA, 2003. [Google Scholar]
  40. Hwang, C.L.; Yoon, K. Multiple Attribute Decision Making: Methods and Applications; Springer: New York, NY, USA, 1981. [Google Scholar]
  41. Yoon, K. A reconciliation among discrete compromise situations. J. Oper. Res. Soc. 1987, 38, 277–286. [Google Scholar] [CrossRef]
  42. Hwang, C.L.; Lai, Y.J.; Liu, T.Y. A new approach for multiple objective decision making. Comput. Oper. Res. 1993, 20, 889–899. [Google Scholar] [CrossRef]
  43. Brans, J.P. L’ingénierie de la décision: Elaboration d’instruments d’aide à la décision. La méthode PROMETHEE. In L’aide à la Décision: Nature, Instruments et Perspectives d’Avenir; Presses de l’Université Laval: Quebec, QC, Canada, 1982. [Google Scholar]
  44. Vatankhah, A.R. The lumped muskingum flood routing model revisited: The storage relationship. Hydrol. Sci. J. 2021, 66, 1625–1637. [Google Scholar] [CrossRef]
  45. Lu, C.; Ji, K.; Wang, W.; Zhang, Y.; Ealotswe, T.K.; Qin, W.; Lu, J.; Liu, B.; Shu, L. Estimation of the interaction between groundwater and surface water based on flow routing using an improved nonlinear muskingum-cunge method. Water Resour. Manag. 2021, 35, 2649–2666. [Google Scholar] [CrossRef]
  46. Farahani, N.; Karami, H.; Farzin, S.; Ehteram, M.; Kisi, O.; El Shafie, A. A new method for flood routing utilizing four-parameter nonlinear muskingum and shark algorithm. Water Resour. Manag. 2019, 33, 4879–4893. [Google Scholar] [CrossRef]
  47. Lee, E.H. Development of a new 8-parameter muskingum flood routing model with modified inflows. Water 2021, 13, 3170. [Google Scholar] [CrossRef]
  48. Akbari, R.; Hessami-Kermani, M.-R.; Shojaee, S. Flood routing: Improving outflow using a new non-linear muskingum model with four variable parameters coupled with PSO-GA algorithm. Water Resour. Manag. 2020, 34, 3291–3316. [Google Scholar] [CrossRef]
  49. Moradi, E.; Yaghoubi, B.; Shabanlou, S. A new technique for flood routing by nonlinear Muskingum model and artificial gorilla troops algorithm. Appl. Water Sci. 2023, 13, 49. [Google Scholar] [CrossRef]
  50. Bozorg-Haddad, O.; Mohammad-Azari, S.; Hamedi, F.; Pazoki, M.; Loáiciga, H.A. Application of a new hybrid non-linear Muskingum modelto flood routing. Proc. Inst. Civ. Eng. Water Manag. 2020, 173, 109–120. [Google Scholar] [CrossRef]
  51. Zlaugotne, B.; Zihare, L.; Balode, L.; Kalnbalkite, A.; Khabdullin, A.; Blumberga, D. Multi-Criteria Decision Analysis Methods Comparison. Environ. Clim. Technol. 2020, 24, 454–471. [Google Scholar] [CrossRef]
  52. Kaya, G.K.; Ozturk, F. A Comparison of the Multi-criteria Decision-Making Methods for the Selection of Researchers. In Industrial Engineering in the Internet-of-Things World: Selected Papers from the Virtual Global Joint Conference on Industrial Engineering and Its Application Areas, GJCIE 2020, Virtually, 14–15 August 2020; Calisir, F., Ed.; Lecture Notes in Management and Industrial Engineering Series; Springer: Cham, Switzerland, 2020. [Google Scholar] [CrossRef]
  53. Vassoney, E.; Mammoliti Mochet, A.; Desiderio, E.; Negro, G.; Pilloni, M.G.; Comoglio, C. Comparing Multi-Criteria Decision-Making Methods for the Assessment of FlowRelease Scenarios From Small Hydropower Plants in the Alpine Area. Front. Environ. Sci. 2021, 9, 635100. [Google Scholar] [CrossRef]
  54. Khan, I.; Pintelon, L.; Martin, H. The Application of Multicriteria Decision Analysis Methods in Health Care: A Literature Review. Med. Decis. Mak. 2022, 42, 262–274. [Google Scholar] [CrossRef] [PubMed]
  55. Cinelli, M.; Burgherr, P.; Kadziński, M.; Słowiński, R. Proper and improper uses of MCDA methods in energy systems analysis. Decis. Support Syst. 2022, 163, 113848. [Google Scholar] [CrossRef]
  56. Zavadskas, E.K.; Kaklauskas, A. The new method of multicriteria evaluation of projects. In Deutsch-Litauisch-Polnisches Kolloquim zum Baubetriebswesen; 3 Jahrgang, Sonderheft; Hochschule fur Technik, Wirtschaft und Kultur: Leipzig, Germany, 1996; pp. 3–8. [Google Scholar]
  57. Zavadskas, E.K.; Kaklauskas, A.; Vilutiene, T. Multicriteria evaluation of apartment blocks maintenance contractors: Lithuanian case study. Int. J. Strateg. Prop. Manag. 2009, 13, 319–338. [Google Scholar] [CrossRef]
Figure 1. Schematic summary of the entire computational process (blue dotted lines show the implementation and output of WCOA, red dotted lines show the details of the WCOA).
Figure 1. Schematic summary of the entire computational process (blue dotted lines show the implementation and output of WCOA, red dotted lines show the details of the WCOA).
Water 15 02588 g001
Figure 2. Flowrate datasets. (a) Wilson’s 1974 flood (SSP), (b) Brutsaert’s data (SSP), (c) River Wye 1960 flood (NSP), (d) Karun river flood (NSP), (e) Viessman and Lewis’s data (MP), (f) Sütçüler flood (MP), (g) River Wye 1982 flood (irregular), and (h) Chenggou and Lingqing river (irregular).
Figure 2. Flowrate datasets. (a) Wilson’s 1974 flood (SSP), (b) Brutsaert’s data (SSP), (c) River Wye 1960 flood (NSP), (d) Karun river flood (NSP), (e) Viessman and Lewis’s data (MP), (f) Sütçüler flood (MP), (g) River Wye 1982 flood (irregular), and (h) Chenggou and Lingqing river (irregular).
Water 15 02588 g002
Figure 3. Estimation errors of the models for each dataset.
Figure 3. Estimation errors of the models for each dataset.
Water 15 02588 g003
Figure 4. Percentage of ranking for both methods ( | | represents OR function; in the overall evaluation of PROMETHEE, model M9 takes 1st or 2nd place with 80%).
Figure 4. Percentage of ranking for both methods ( | | represents OR function; in the overall evaluation of PROMETHEE, model M9 takes 1st or 2nd place with 80%).
Water 15 02588 g004
Figure 5. Cumulative score of each model by TOPSIS and PROMETHEE.
Figure 5. Cumulative score of each model by TOPSIS and PROMETHEE.
Water 15 02588 g005
Figure 6. Consistency in ranking for each hydrograph type and scenario between TOPSIS and PROMETHEE (CR is consistency rate, CR(0) shows the exact coherency, CR(1) shows the coherency in the gray zone. For scenario 1 in SSP, 50% for CR(0) shows that 5 out of 10 models obtain the same rankings from both methods, whereas 90% for CR(1) shows that 9 out of 10 models obtain close rankings, located in the gray zone).
Figure 6. Consistency in ranking for each hydrograph type and scenario between TOPSIS and PROMETHEE (CR is consistency rate, CR(0) shows the exact coherency, CR(1) shows the coherency in the gray zone. For scenario 1 in SSP, 50% for CR(0) shows that 5 out of 10 models obtain the same rankings from both methods, whereas 90% for CR(1) shows that 9 out of 10 models obtain close rankings, located in the gray zone).
Water 15 02588 g006
Table 1. Muskingum models with number of parameters and references.
Table 1. Muskingum models with number of parameters and references.
Model No.Storage Equation *# of ParametersStudied by
M1 S = K [ w I + ( 1 w ) Q ] 2[5]
M2 S = K [ w I α + ( 1 w ) Q α ] 3 [6]
M3 S = K [ w I + ( 1 w ) Q ] m 3 [7]
M4 S = K [ w I + ( 1 w ) Q ] 3 / 5 2 [3]
M5 S = K [ w ( 1 + γ ) I + ( 1 w ) Q ] 3 [8]
M6 S = K [ w I α + ( 1 w ) Q β ] 4[9]
M7 S = K [ w I α + ( 1 w ) Q α ] m 4 [2]
M8 S = K [ w I α + ( 1 w ) Q β ] m 5 [10]
M9 S = K [ w ( 1 + γ ) I + ( 1 w ) Q ] m 4 [11]
M10 S = K [ w ( 1 + γ ) I α + ( 1 w ) Q β ] m 6Proposed
* where α is the flow exponent, m is the general model exponent, γ is the lateral flow, and β is the flow exponent.
Table 3. Ranking results of models for criteria A d R 2 and N-RMSE.
Table 3. Ranking results of models for criteria A d R 2 and N-RMSE.
A d R 2
SSPNSPMPIrregularOverall
ModelDS1DS2ARDS3DS4ARDS5DS6ARDS7DS8ARAR
M1898.59109.5787.510478.125
M2766.57778467566.375
M3586.5565.59998677
M41010101099.51010109109.59.75
M59788461112224.25
M6634.5687655.5475.55.625
M7322.5354465586.54.5
M8211.5433.5576697.54.625
M94442222223122.5
M101531113331322.25
N-RMSE
M19999109.59991099.59.25
M2766.5787.5777766.56.875
M3586.5565.5888877.56.875
M41010101099.51010109109.59.75
M5877.5856.5322.5322.54.75
M6645676.5655.54865.75
M7322.5444565.5655.54.375
M8211.5333444544.53.25
M9454.5222232.5211.52.625
M101321111111321.5
Table 4. Decision matrix for four hydrograph types: SSP, NSP, MP, and irregular.
Table 4. Decision matrix for four hydrograph types: SSP, NSP, MP, and irregular.
SSPNSP
Model A d R 2 N-RMSENPMBModel A d R 2 N-RMSENPMB
M18.5921M19.5921
M26.56.530M27730
M36.56.530M35.5530
M4101021M49.5921
M587.531M56631
M64.5540M67640
M72.52.540M74440
M81.51.550M83.5350
M944.540M92240
M103260M101160
MPIrregular
Model A d R 2 N-RMSENPMBModel A d R 2 N-RMSENPMB
M17.5921M179.521
M26730M266.530
M39830M377.530
M4101021M49.59.521
M512.531M522.531
M65.55.540M65.5640
M755.540M76.55.540
M86450M87.54.550
M922.540M921.540
M103160M102260
Table 5. Ranking results (starting from the best model) and scores of TOPSIS and PROMETHEE for three scenarios and four hydrograph types.
Table 5. Ranking results (starting from the best model) and scores of TOPSIS and PROMETHEE for three scenarios and four hydrograph types.
Dataset TypeScenario 1Scenario 2Scenario 3
TOPSISPROMETHEETOPSISPROMETHEETOPSISPROMETHEE
Model χ Model ξ Model χ Model ξ Model χ Model ξ
SSP80.547580.125180.597180.169380.697080.2577
50.539870.115770.577970.146470.673070.2078
70.531150.0507100.5477100.0353100.6481100.1245
10.511610.041350.495450.012490.565690.0477
100.498110−0.009390.491890.009260.526660.0020
40.46849−0.010110.46221−0.010550.40875−0.0642
90.45366−0.046160.46136−0.030120.38501−0.1141
60.42684−0.048540.41794−0.108530.38502−0.1165
20.35252−0.109420.36372−0.111810.36263−0.1165
30.35253−0.109430.36373−0.111840.31584−0.2284
NSP50.593950.1586100.581290.1856100.6834100.2618
100.530790.151690.5771100.152990.672390.2536
90.5303100.098550.555550.130180.590480.0975
10.469370.007880.503780.032070.545950.0730
40.46938−0.000770.476170.028850.483670.0706
80.45931−0.012610.41883−0.033330.42933−0.0250
70.43964−0.012640.41881−0.069310.31666−0.1810
30.38043−0.037530.39734−0.069340.31661−0.1827
20.28492−0.163360.28352−0.170660.30234−0.1827
60.27306−0.189920.28146−0.186920.27472−0.1851
MP50.885350.467150.885650.466750.886350.4657
90.536190.120790.582190.151990.674390.2142
10.516110.0883100.5612100.0593100.6577100.1525
100.5125100.012710.468710.040780.47927−0.0451
40.45564−0.030680.41997−0.070470.46701−0.0543
70.38987−0.083070.41744−0.088960.44548−0.0651
80.38656−0.100040.40546−0.088910.37576−0.0667
60.37462−0.105460.40008−0.107420.37082−0.1114
20.34508−0.128520.35412−0.107440.30474−0.2056
30.24443−0.241230.23193−0.255630.20723−0.2843
Irregular50.896650.469750.899550.469450.905150.4690
90.543390.195490.591290.233390.690090.3093
10.534910.0611100.5528100.1014100.6574100.2016
100.5020100.051310.487210.011160.42131−0.0889
40.48964−0.040740.43872−0.097270.39866−0.0910
20.35152−0.096160.37434−0.100080.39792−0.0995
60.34946−0.119120.36346−0.109710.39197−0.1185
70.33427−0.140770.35667−0.133320.38728−0.1634
80.31933−0.175080.34753−0.183340.33443−0.2000
30.29658−0.205830.29568−0.191730.29394−0.2185
Table 6. MCDA results of all flood events. (a) Decision matrix obtained from 8 flood events and (b) corresponding scores of TOPSIS and PROMETHEE for three scenarios.
Table 6. MCDA results of all flood events. (a) Decision matrix obtained from 8 flood events and (b) corresponding scores of TOPSIS and PROMETHEE for three scenarios.
(a)
Model A d R 2 N-RMSENPMB
M18.1259.2521
M26.3756.87530
M376.87530
M49.759.7521
M54.254.7531
M65.6255.7540
M74.54.37540
M84.6253.2550
M92.52.62540
M102.251.560
(b)
Scenario 1Scenario 2Scenario 3
TOPSISPROMETHEETOPSISPROMETHEETOPSISPROMETHEE
Model χ Model ξ Model χ Model ξ Model χ Model ξ
50.753950.291250.734650.274750.700850.2418
90.527190.128790.574290.1606100.672290.2244
10.5201100.0556100.5687100.106190.6711100.2071
100.517910.029280.494670.001080.583380.0431
40.48217−0.017670.48048−0.014670.555870.0382
80.44988−0.043510.47071−0.023760.44736−0.0849
70.44204−0.055640.43136−0.104520.37171−0.1295
60.36866−0.114460.39594−0.116210.37012−0.1352
20.34292−0.124120.35272−0.127830.34203−0.1676
30.32483−0.149530.33063−0.155640.32784−0.2374
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Şahin, A.U.; Özkaya, A. Multi-Criteria Decision Analyses for the Selection of Hydrological Flood Routing Models. Water 2023, 15, 2588. https://doi.org/10.3390/w15142588

AMA Style

Şahin AU, Özkaya A. Multi-Criteria Decision Analyses for the Selection of Hydrological Flood Routing Models. Water. 2023; 15(14):2588. https://doi.org/10.3390/w15142588

Chicago/Turabian Style

Şahin, Abdurrahman Ufuk, and Arzu Özkaya. 2023. "Multi-Criteria Decision Analyses for the Selection of Hydrological Flood Routing Models" Water 15, no. 14: 2588. https://doi.org/10.3390/w15142588

APA Style

Şahin, A. U., & Özkaya, A. (2023). Multi-Criteria Decision Analyses for the Selection of Hydrological Flood Routing Models. Water, 15(14), 2588. https://doi.org/10.3390/w15142588

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