Abstract
The probability distribution function (PDF)-based unit hydrographs (UHs) are gaining momentum in an application for more accurate rainfall-runoff transformation. Employing seven statistical performance indices, R2, NSE, MSE, RMSE, MAE, MAPE, and SE in generalized reduced gradient nonlinear programming (GRG-NLP) optimization, 18 known and 12 adaptable PDF-based UHs were assessed against UHs derived from 18 storms in 7 basins across the United States, Turkey, and India. To this end, 27 Maple codes were proposed for UH-application requiring only peak discharge (qp), time to peak (tp), and time base (tb) for derivation. The introduced PDFs, such as Dagum, Generalized Gamma, Log-Logistic, Gumbel Type-I, and Shifted Gompertz, replicated the observed data-derived UHs more closely than did the known PDFs, like Inverse Gaussian, two-parameter gamma distribution (2-PGD), Log-Normal, Inverse-Gamma, and Nagakami. Among the three-parameter (6 nos.), two-parameter (21 nos.), and single-parameter (3 nos.) PDFs, the Dagum, Log-Logistic, and Poisson consistently outperformed their respective counterparts in replication.
HIGHLIGHTS
The present study evaluates the adequacy of 30 (18 old +12 new) probability density functions (PDFs) for unit hydrograph development.
The novel PDFs are more promising than existing PDFs in the literature.
Proposed PDFs with simple formulaic structures are easily adaptable as executables using a calculator.
The incomplete-Gamma distribution has the highest ability for qp and tp.
Maple codes are provided for 27 PDFs.
LIST OF ABBREVIATIONS
Probability distribution function
- UH
Unit hydrograph
- qp
Peak discharge
- tp
Time to peak
- tb
Time base
- qv
Runoff volume
- W50
Width of UH measured at 50% of peak discharge ordinates
- W75
Width of UH measured at 75% of peak discharge ordinates
- WR50
Width of the rising limb of UH measured at 50% of peak discharge ordinate
- WR75
Width of the rising limb of UH measured at 75% of peak discharge ordinate
- L
Length of the longest water course
- LC
Length along the main stream from the outlet to the centroid of the basin
- β
Non-dimensional parameter
- RSO
Overall shape-based rank
- RS
Event-wise shape-based rank
- RqpO
Overall qp-based rank
- RtpO
Overall tp-based rank
- Rqp
Event-wise qp-based rank
- Rtp
Event-wise tp-based rank
- MSE
Mean square error
- RMSE
Root mean square error
- MAE
Mean absolute error
- MAPE
Mean absolute percent error
- SE
Standard error absolute
- RE
Relative error
- R2
Coefficient of determination
- NSE
Nash–Sutcliffe efficiency
- RNSE
Event-wise NSE-based rank
- RR2
Event-wise R2-based rank
- RMSE
Event-wise MSE-based rank
- RRMSE
Event-wise RMSE-based rank
- RMAE
Event-wise MAE-based rank
- RMAPE
Event-wise MAPE-based rank
- RSE
Event-wise SE-based rank
- RRE
Event-wise RE-based rank
INTRODUCTION
To design and operate effectively the stormwater management facilities, treatment plants, hydraulic structures, flood control projects, and ground water development works, it is essential to accurately measure the complex, dynamic, and nonlinear rainfall-runoff transformation process. This complexity is a result of the considerable variability in rainfall patterns and physical characteristics (Tokar & Markus 2000; Ito et al. 2006; Geetha et al. 2008; Yoshitani et al. 2009; Bhadra et al. 2010; Brooks et al. 2012). However, understanding and quantifying this process remain challenging due to the difficulty in directly measuring various water cycle fluxes, especially over large regions where instrumentation may be limited (McDonnell & Tanaka 2001). Consequently, our understanding of the water cycle relies heavily on mathematical models.
Analyzing the transformation process using the unit hydrograph (UH) and its variants has been a focus in hydrologic studies since Sherman (1932) introduced it. Accurately predicting UH parameters, such as peak discharge (qp), time to peak (tp), time base (tb), and runoff volume (qv), is crucial for hydrologic investigations. The UH, representing a linear system, has been extensively studied by various scholars (Nash 1957; Dooge 1959, 1973; Chow et al. 1988; Singh 1988). Various empirical UH models, ranging from single curves to a specific family of equations exist in the literature (Langbein 1940; Commons 1942; Williams 1945; Mockus 1957; Bender & Roberson 1961; SCS 1972; Singh 2000; Bhunya et al. 2003). Single curves can represent UHs for a region or specific watershed sizes. Fixing the scale (qp or tp) allows a standard non-dimensional curve to define a complete UH under unit volumetric conditions. Alternatively, manually drawing a curve conserving unit volume through parameters like qp, tp, tb, W50, W75, WR50, and WR75 can be error-prone. Flexible 2-parameter models, representing a family of curves, capture UH more accurately than do 1-parameter models. However, empirical equations represent UH shapes more precisely than do the family of curves. The adoption of probability distribution functions (PDFs) as an analytical tool to mimic UHs is inspired by their similarities in shape, unit volume, and positive ordinates. The robust mathematical foundation and enhanced hydrological understanding also advocate for the applicability of PDFs in representing UHs.
Empirical hydrographs have been proposed by many researchers, viz., first empirical UH by Edson (1951); two-parameter gamma distribution (2PGD) or Nash (1957, 1958, 1959, 1960) model, Gray (1961, 1962), Reich (1962), Wu (1963), Lienhard (1964), DeCoursey (1966), Lienhard & Meyer (1967), Lienhard & Davis (1971), Gupta & Moin (1974), Gupta et al. (1974), Sokolov et al. (1976), Haan (1977), Cruise & Contractor (1978), Croley (1980), Aron & White (1982), Singh (1981, 1982), Collins (1983), Rosso (1984), Phien & Jivajirajah (1984), Yevjevich & Obeysekera (1984), Ciepielowski (1987), James et al. (1987), Singh (1987), Koutsoyiannis & Xanthopoulos (1989), Haktanir & Sezen (1990), Meadows & Ramsey (1991a, 1991b), Haan et al. (1994), Gottschalk & Weingartner (1998), Evans et al. (2001), Singh (2000, 2004, 2005, 2006, 2007a, 2007b), Yue et al. (2002), inverse Gaussian instantaneous UH (IUH) by Bardsley (1983), Gamma- and log-normal-UH by Bhattacharjya (2004), Gamma-, Beta- Chi-square- and Weibull-synthetic UHs (SUHs) by Bhunya et al. (2003, 2004, 2005, 2007, 2008, 2009), Sahoo et al. (2006), 11 PDF-based UH along with Maple codes by Nadarajah (2007), Rai et al. (2007, 2008, 2009), Nakagami-m and seven other PDF-based UHs by Rai et al. (2010), entropy-based IUH specialized into Gamma, Lienhard and Nakagami-m PDFs by Singh (2011), Singh et al. (2014), and Gamma, Gumbel, Log-normal, Normal, Weibull, and 3-parameter Pearson PDFs transmuted UH by Ghorbani et al. (2017), etc. which laid the base of hydro-physical perception underlying the integration of PDFs in synthesizing hydrographs.
Nevertheless, the utility of PDFs extends beyond hydrologic studies (Abramowitz et al. 1988; Johnson et al. 1994, 1995; Walck 1996; Forbes et al. 2011; Krishnamoorthy 2016). The presented PDFs find application in various fields, for instance, the Dagum PDF proves valuable in actuarial sciences, survival and reliability analysis, atmospheric modeling, and meteorology (Benjamin et al. 2013; Oluyede & Rajasooriya 2013; Shahzad & Asghar 2016). The Log-Logistic PDF is relevant in survival analysis, finance, economics, networking, hydrology (for fitting long-duration extremes), and diverse human-centric domains such as biology, epidemiology, psychology, technology, and energy (Shoukri et al. 1988; Aryal 2013; Gago-Benítez et al. 2013). The Gumbel PDF finds applications in ocean wave modeling, wind and electrical engineering, earthquakes' thermodynamics, financial risk assessment, probability analysis of floods/rainfall/droughts, survival analysis, meteorology, air pollution, and geology (Kotz & Nadarajah 2000). Gamma family PDFs serve engineering, hydrology, survival and reliability analysis, metallurgy, and Bayesian statistics purposes (Agarwal & Kalla 1996; Aksoy 2000; Barriga et al. 2018). The Shifted Gompertz PDF (Bemmaor 1992) is employed for modeling and forecasting adoption timing and diffusion of innovations in markets, as well as predicting the growth and decline of services or in marketing decision models (Jiménez & Jodrá 2008; Bauckhage & Kersting 2014). Other PDFs have specific applications: binomial PDF is suitable for analyzing repeated independent trials, risk management, and quality assurance (Friedman et al. 1983; Consul 1990). The negative binomial PDF is used in failure analysis, machine reliability, RNA and DNA sequencing, meteorological studies, or as an alternative to the Poisson PDF (Bain & Wright 1982). The Poisson PDF is exercised in ET Engg., internet traffic assessment, mortality rate studies, stock markets, genetics, biotechnology, and physics (Clarke 1946). The Gompertz PDF is used to investigate age-specific mortality (Pollard & Valkovics 1992). The Rice PDF is found worthy for sonar echo- and radar signal-analysis and magnetic resonance imaging (Talukdar & Lawing 1991; Carobbi & Cati 2008). The Erlang PDF has utility in telecommunication traffic load assessment, disease analysis, cell cycle time distribution, and business economics (Gavagnin et al. 2019). Notably, with the exception of a few gamma family PDFs such as Generalized- and Log-Gamma, the introduced PDFs have not been specifically adopted in UH investigations so far. PDFs enhance flood forecasting and agricultural planning by improving risk assessment and resource allocation. However, their effectiveness relies on accurate data, appropriate assumptions, and an understanding of uncertainties. PDFs offer insights into flood probabilities for tailored planning, yet limitations include data dependence, sensitivity to assumptions, and computational complexity. Reliability may be compromised in regions with limited monitoring, but a careful consideration of these factors is crucial for successful real-world application.
The method of moments (MOM) often prioritizes the extremities over the center of a PDF, leading to errors, especially at the peak (Nash 1957). It could be addressed through time-consuming linear programming or least-squares approaches, particularly when dealing with oscillatory UHs where computational effort escalates with the number of storms (Bree 1978; Singh 2007b). Genetic algorithms do not outperform nonlinear optimization methods in enhancing PDF performance (Ghorbani et al. 2017). Therefore, the study aims to (a) evaluate the proposed PDFs for UH fitting; (b) determine optimum PDF parameters using the Excel-based generalized reduced gradient (GRG) nonlinear optimization algorithm (Lasdon et al. 1978); and (3) scrutinize the adequacy of PDFs in defining the UH and its properties. The investigation delves into the suitability of both existing and additional PDFs as UHs for more accurate rainfall-runoff transformation.
STUDY AREA AND DATA USED
Eighteen storm-based UHs of seven basins were used to examine PDFs. The salient UH and watershed characteristics obtained from the literature are given in Table 1.
Salient unit hydrograph and watershed characteristics
Sl. No. . | Watershed . | Catchment area A (km2) . | Length of main stream L (km) . | Length from outlet to centroid LC (km) . | Slope . | Peak discharge qp (m3/s/cm) . | Time to peak tp (h) . | Time base tb (h) . |
---|---|---|---|---|---|---|---|---|
1. | Myntdu-Leska, Meghalaya, India | 350.0 | 51.778 | 27.8 | 4.2 km/m | 11.829 | 5 | 30 |
2. | Bridge Catchment No. 253 Narmada River, Madhya Pradesh, India | 114.22 | 35.42 | 16.9 | 2.5 km/m | 5.46 | 5 | 21 |
3. | Karabalcik Creek, Anatolia, Turkey, Basin 16 | 10.6 | 5 | 2.65 | 2.7% | 12.84 | 1.75 | 14.4 |
4. | Inderesi Creek, Anatolia, Turkey, Basin 18 | 98.0 | 17 | 7.8 | 2.0% | 127.5 | 2.22 | 12.5 |
5. | Gormel, Ermenek Creek, Anatolia, Turkey, Basin 17 | 141.5 | 23.5 | 12 | 1.18% | 54 | 6 | 61 |
6. | Kurtleravsari, Aksu, Anatolia, Turkey, Basin 20 | 3360 | 123 | 60 | 0.43% | 301 | 18 | 113 |
Sl. No. . | Watershed . | Catchment area A (km2) . | Length of main stream L (km) . | Length from outlet to centroid LC (km) . | Slope . | Peak discharge qp (m3/s/cm) . | Time to peak tp (h) . | Time base tb (h) . |
---|---|---|---|---|---|---|---|---|
1. | Myntdu-Leska, Meghalaya, India | 350.0 | 51.778 | 27.8 | 4.2 km/m | 11.829 | 5 | 30 |
2. | Bridge Catchment No. 253 Narmada River, Madhya Pradesh, India | 114.22 | 35.42 | 16.9 | 2.5 km/m | 5.46 | 5 | 21 |
3. | Karabalcik Creek, Anatolia, Turkey, Basin 16 | 10.6 | 5 | 2.65 | 2.7% | 12.84 | 1.75 | 14.4 |
4. | Inderesi Creek, Anatolia, Turkey, Basin 18 | 98.0 | 17 | 7.8 | 2.0% | 127.5 | 2.22 | 12.5 |
5. | Gormel, Ermenek Creek, Anatolia, Turkey, Basin 17 | 141.5 | 23.5 | 12 | 1.18% | 54 | 6 | 61 |
6. | Kurtleravsari, Aksu, Anatolia, Turkey, Basin 20 | 3360 | 123 | 60 | 0.43% | 301 | 18 | 113 |
Indian catchments
Narmada (Gondia-Jabalpur) Railway Bridge Catchment No. 253 (Lohani et al. 2001).
Narmada (Gondia-Jabalpur) Railway Bridge Catchment No. 253 (Lohani et al. 2001).
Turkish watersheds
Basin of Anatolia-Turkey (a) Karabalcik Creek, Basin 16 (b) Inderesi Creek, Basin 18 (c) Gormel, Ermenek Creek, Basin 17 (d) Kurtleravsari, Aksu, Basin 20 (Haktanir & Sezen 1990).
Basin of Anatolia-Turkey (a) Karabalcik Creek, Basin 16 (b) Inderesi Creek, Basin 18 (c) Gormel, Ermenek Creek, Basin 17 (d) Kurtleravsari, Aksu, Basin 20 (Haktanir & Sezen 1990).
U.S.A. catchment
ESTIMATION OF DISTRIBUTION PARAMETERS
The basic equations for estimating 2-parameter PDF, f(x) parameters with known qp and tp are: Mode = tp, f(Mode) = qp, β = qp.tp. For 3-parameter PDF, one of the scale parameters is assumed as tp or tb. The proposed 27 Maple codes (Appendix-Table A.1) estimate unknown PDF parameters by solving the expressions (Table 2) with known qp, tp, and tb.
Probability distribution function and it's time to peak and peak discharge
Sl.No. . | Distribution and reference . | PDF ![]() | tp . | qp . |
---|---|---|---|---|
1 | Shifted Gompertz (Torres 2014) | ![]() | ![]() ![]() | ![]() |
2 | Dagum (Kleiber 2008) | ![]() | ![]() | ![]() |
3 | Rice (Talukdar & Lawing 1991) | ![]() | ![]() | ![]() |
4 | Rayleigh (Cleveland et al. 2006) | ![]() | ![]() | ![]() |
5 | Erlang (Mudasir & Ahmad 2017) | ![]() | ![]() | ![]() |
6 | Incomplete-Gamma (Sade 2001) | ![]() | ![]() | ![]() |
7 | Kumaraswamy (Nadarajah 2007) | ![]() | ![]() | ![]() |
8 | Weibull (Nadarajah 2007) | ![]() | ![]() | ![]() |
9 | Log-Normal (Nadarajah 2007) | ![]() | ![]() | ![]() |
10 | Fréchet (Nadarajah 2007) | ![]() | ![]() | ![]() |
11 | Nagakami (Yacoub et al. 1999) | ![]() | ![]() | ![]() |
12 | Beta (Nadarajah 2007) | ![]() | ![]() | ![]() |
13 | Inverse-Gamma (Nadarajah 2007) | ![]() | ![]() | ![]() |
14 | Hybrid/Nash Modified (Bhunya et al. 2005) | ![]() | ![]() | ![]() |
15 | 3 Parameter 2 Sided Power (Nadarajah 2007) | ![]() | ![]() | ![]() |
16 | Simple 3 Parameter (Singh 2015) | ![]() | ![]() | ![]() |
17 | Binomial (Edwards 1960) | ![]() ![]() | ![]() | ![]() |
18 | Negative- Binominal (Fisher 1941) | ![]() | ![]() | ![]() |
19 | Normal (Ghorbani et al. 2017) | ![]() | ![]() | ![]() |
20 | Gumbel Type-I (Chakraborty & Chakravarty 2014) | ![]() ![]() | ![]() | ![]() |
21 | Poisson (Sadooghi-Alvandi, 1990) | ![]() | ![]() | ![]() |
22 | 2-PGD (Bhunya et al. 2003) | ![]() | ![]() | ![]() |
23 | Logistic (Balakrishnan 1991) | ![]() | ![]() | ![]() |
24 | Log-Logistic/Fisk (Ashkar & Mahdi 2006) | ![]() | ![]() | ![]() |
25 | Log-Gamma (Consul & Jain 1971) | ![]() | ![]() | ![]() |
26 | Gumbel Type II (Gumbel 1960) | ![]() | ![]() | ![]() |
27 | Chi-Square (Nadarajah 2007) | ![]() | ![]() | ![]() |
28 | Generalized- Gamma (Stacy & Mihram, 1965) | ![]() | ![]() | ![]() |
29 | Inverse Gaussian (Nadarajah 2007) | ![]() | ![]() | ![]() |
30 | Gompertz (Dey et al. 2018) | ![]() | ![]() ![]() | ![]() |
Sl.No. . | Distribution and reference . | PDF ![]() | tp . | qp . |
---|---|---|---|---|
1 | Shifted Gompertz (Torres 2014) | ![]() | ![]() ![]() | ![]() |
2 | Dagum (Kleiber 2008) | ![]() | ![]() | ![]() |
3 | Rice (Talukdar & Lawing 1991) | ![]() | ![]() | ![]() |
4 | Rayleigh (Cleveland et al. 2006) | ![]() | ![]() | ![]() |
5 | Erlang (Mudasir & Ahmad 2017) | ![]() | ![]() | ![]() |
6 | Incomplete-Gamma (Sade 2001) | ![]() | ![]() | ![]() |
7 | Kumaraswamy (Nadarajah 2007) | ![]() | ![]() | ![]() |
8 | Weibull (Nadarajah 2007) | ![]() | ![]() | ![]() |
9 | Log-Normal (Nadarajah 2007) | ![]() | ![]() | ![]() |
10 | Fréchet (Nadarajah 2007) | ![]() | ![]() | ![]() |
11 | Nagakami (Yacoub et al. 1999) | ![]() | ![]() | ![]() |
12 | Beta (Nadarajah 2007) | ![]() | ![]() | ![]() |
13 | Inverse-Gamma (Nadarajah 2007) | ![]() | ![]() | ![]() |
14 | Hybrid/Nash Modified (Bhunya et al. 2005) | ![]() | ![]() | ![]() |
15 | 3 Parameter 2 Sided Power (Nadarajah 2007) | ![]() | ![]() | ![]() |
16 | Simple 3 Parameter (Singh 2015) | ![]() | ![]() | ![]() |
17 | Binomial (Edwards 1960) | ![]() ![]() | ![]() | ![]() |
18 | Negative- Binominal (Fisher 1941) | ![]() | ![]() | ![]() |
19 | Normal (Ghorbani et al. 2017) | ![]() | ![]() | ![]() |
20 | Gumbel Type-I (Chakraborty & Chakravarty 2014) | ![]() ![]() | ![]() | ![]() |
21 | Poisson (Sadooghi-Alvandi, 1990) | ![]() | ![]() | ![]() |
22 | 2-PGD (Bhunya et al. 2003) | ![]() | ![]() | ![]() |
23 | Logistic (Balakrishnan 1991) | ![]() | ![]() | ![]() |
24 | Log-Logistic/Fisk (Ashkar & Mahdi 2006) | ![]() | ![]() | ![]() |
25 | Log-Gamma (Consul & Jain 1971) | ![]() | ![]() | ![]() |
26 | Gumbel Type II (Gumbel 1960) | ![]() | ![]() | ![]() |
27 | Chi-Square (Nadarajah 2007) | ![]() | ![]() | ![]() |
28 | Generalized- Gamma (Stacy & Mihram, 1965) | ![]() | ![]() | ![]() |
29 | Inverse Gaussian (Nadarajah 2007) | ![]() | ![]() | ![]() |
30 | Gompertz (Dey et al. 2018) | ![]() | ![]() ![]() | ![]() |
PERFORMANCE EVALUATION








Ranking of PDFs based on agreement with observed UH shape (S, symmetric, PS, positively skewed, NS, negatively skewed)
No . | Catchment/Event . | RS for Bixler Run Basin (U.S.A) . | RS for Turkish Basin . | RS-Indian . | Average of rank (RS) . | Overall rank (RSO) . | |||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | MM DD YYYY PDF (No. of Parameters) . | 02-07-1965 (PS) . | 02-26-1957 (PS) . | 03-06-1963 (NS) . | 03-14-1956 (S) . | 03-17-1963 (PS) . | 04-16-1961 (PS) . | 05-18-1963 (S) . | 06-06-1964 (PS) . | 06-21-1956 (PS) . | 10-04-1962 (PS) . | 10-18-1967 (PS) . | 12-06-1962 (S) . | Karabalcik . | Kurtleravsari . | Inderesi . | Gormel . | Myntdu . | Br. No.253 . | ||
. | R2 w.r.t. overall rank → . | 0.86 . | 0.75 . | 0.10 . | 0.80 . | 0.83 . | 0.87 . | 0.82 . | 0.82 . | 0.53 . | 0.79 . | 0.80 . | 0.85 . | 0.50 (PS) . | 0.49 (PS) . | 0.67 (PS) . | 0.82 (PS) . | 0.62 (PS) . | 0.76 (PS) . | ||
1 | Shifted Gompertz (2P) | 3 | 2 | 21 | 9 | 7 | 4 | 8 | 3 | 4 | 4 | 5 | 4 | 7 | 6 | 5 | 2 | 7 | 3 | 5.78 | 5 |
2 | 2-PGD (2P) | 10 | 10 | 15 | 2 | 2 | 9 | 5 | 9 | 3 | 10 | 11 | 7 | 11 | 3 | 8 | 4 | 8 | 7 | 7.44 | 7 |
3 | Dagum (3P) | 1 | 4 | 1 | 3 | 3 | 1 | 2 | 1 | 8 | 1 | 1 | 1 | 2 | 4 | 1 | 6 | 1 | 1 | 2.33 | 1 |
4 | Rice (2P) | 14 | 15 | 4 | 10 | 11 | 15 | 14 | 14 | 19 | 17 | 17 | 17 | 18 | 17 | 21 | 17 | 22 | 18 | 15.56 | 13 |
5 | Rayleigh (1P) | 30 | 30 | 30 | 30 | 30 | 30 | 30 | 30 | 18 | 30 | 27 | 30 | 17 | 18 | 23 | 25 | 19 | 23 | 26.11 | 28 |
6 | Erlang (2P) | 15 | 14 | 17 | 19 | 18 | 12 | 15 | 11 | 14 | 18 | 13 | 19 | 25 | 12 | 28 | 29 | 25 | 29 | 18.50 | 19 |
7 | Incomplete Gamma (2P) | 20 | 12 | 27 | 6 | 6 | 17 | 22 | 23 | 25 | 14 | 12 | 13 | 20 | 10 | 13 | 12 | 17 | 24 | 16.28 | 15 |
8 | Kumaraswamy (3P) | 18 | 21 | 7 | 21 | 22 | 20 | 19 | 21 | 20 | 23 | 19 | 20 | 16 | 15 | 19 | 16 | 20 | 19 | 18.67 | 22 |
9 | Log-Normal (2P) | 6 | 6 | 19 | 11 | 10 | 8 | 7 | 10 | 7 | 9 | 8 | 8 | 4 | 7 | 3 | 5 | 3 | 6 | 7.61 | 8 |
10 | Weibull (2P) | 21 | 23 | 6 | 17 | 20 | 23 | 21 | 22 | 16 | 24 | 22 | 23 | 14 | 13 | 17 | 14 | 18 | 15 | 18.28 | 17 |
11 | Fréchet (2P) | 19 | 19 | 26 | 25 | 21 | 14 | 23 | 16 | 22 | 12 | 8 | 21 | 5 | 24 | 18 | 24 | 14 | 21 | 18.44 | 18 |
12 | Nakagami (2P) | 11 | 13 | 10 | 4 | 8 | 11 | 10 | 13 | 10 | 13 | 15 | 10 | 12 | 14 | 16 | 9 | 15 | 9 | 11.28 | 10 |
13 | Beta (3P) | 12 | 20 | 8 | 18 | 15 | 18 | 19 | 20 | 11 | 20 | 16 | 21 | 13 | 8 | 14 | 10 | 13 | 11 | 14.83 | 12 |
14 | Inverse-Gamma (2P) | 8 | 7 | 23 | 16 | 12 | 7 | 12 | 5 | 12 | 2 | 3 | 11 | 1 | 21 | 10 | 11 | 9 | 10 | 10.00 | 9 |
15 | Hybrid/Nash Modified (2P) | 28 | 28 | 28 | 28 | 26 | 26 | 28 | 28 | 1 | 27 | 21 | 28 | 10 | 2 | 9 | 21 | 6 | 13 | 19.89 | 24 |
16 | 3 Para.2 Sided Power (3P) | 26 | 24 | 25 | 23 | 24 | 26 | 25 | 24 | 28 | 25 | 29 | 27 | 30 | 26 | 30 | 30 | 30 | 30 | 26.78 | 29 |
17 | Simple 3 Parameter (3P) | 29 | 29 | 29 | 29 | 29 | 29 | 29 | 29 | 24 | 29 | 30 | 29 | 21 | 11 | 27 | 28 | 23 | 27 | 26.72 | 30 |
18 | Poisson (1P) | 22 | 16 | 12 | 20 | 13 | 21 | 10 | 19 | 15 | 15 | 22 | 14 | 24 | 29 | 21 | 18 | 26 | 17 | 18.56 | 21 |
19 | Binomial (2P) | 17 | 8 | 13 | 22 | 23 | 16 | 18 | 18 | 21 | 18 | 20 | 15 | 23 | 30 | 25 | 23 | 28 | 13 | 19.50 | 23 |
20 | Negative Binominal (2P) | 23 | 22 | 14 | 12 | 16 | 22 | 12 | 14 | 13 | 11 | 24 | 12 | 28 | 26 | 15 | 15 | 12 | 20 | 17.28 | 16 |
21 | Normal (2P) | 16 | 18 | 3 | 13 | 14 | 19 | 17 | 17 | 26 | 21 | 18 | 18 | 19 | 22 | 24 | 22 | 24 | 22 | 18.50 | 19 |
22 | Gumbel Type-I (2P) | 2 | 1 | 20 | 7 | 5 | 3 | 6 | 2 | 5 | 5 | 4 | 3 | 8 | 9 | 7 | 1 | 10 | 3 | 5.61 | 4 |
23 | Logistic (2P) | 13 | 17 | 5 | 15 | 17 | 13 | 16 | 12 | 23 | 16 | 14 | 16 | 15 | 20 | 20 | 19 | 21 | 16 | 16.00 | 14 |
24 | Log-Logistic (2P) | 5 | 9 | 16 | 5 | 4 | 2 | 2 | 4 | 9 | 6 | 2 | 2 | 3 | 5 | 2 | 7 | 2 | 2 | 4.83 | 3 |
25 | Log-Gamma (2P) | 7 | 11 | 22 | 14 | 19 | 10 | 9 | 6 | 17 | 3 | 6 | 9 | 27 | 19 | 12 | 12 | 11 | 12 | 12.56 | 11 |
26 | Generalized Gamma (3P) | 9 | 5 | 2 | 1 | 1 | 6 | 1 | 8 | 2 | 8 | 10 | 5 | 9 | 1 | 4 | 3 | 4 | 5 | 4.67 | 2 |
27 | Inverse Gaussian (2P) | 4 | 3 | 18 | 8 | 9 | 5 | 4 | 7 | 6 | 7 | 7 | 6 | 6 | 16 | 6 | 8 | 5 | 8 | 7.39 | 6 |
28 | Gumbel Type-II (2P) | 24 | 25 | 9 | 27 | 25 | 24 | 24 | 25 | 29 | 26 | 25 | 24 | 22 | 23 | 26 | 26 | 27 | 26 | 24.28 | 25 |
29 | Gompertz (2P) | 25 | 26 | 11 | 26 | 27 | 25 | 25 | 26 | 30 | 28 | 26 | 26 | 26 | 25 | 29 | 27 | 29 | 28 | 25.83 | 27 |
30 | Chi-square (1P) | 27 | 26 | 24 | 24 | 28 | 28 | 25 | 27 | 27 | 22 | 28 | 25 | 29 | 28 | 11 | 20 | 16 | 25 | 24.44 | 26 |
No . | Catchment/Event . | RS for Bixler Run Basin (U.S.A) . | RS for Turkish Basin . | RS-Indian . | Average of rank (RS) . | Overall rank (RSO) . | |||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | MM DD YYYY PDF (No. of Parameters) . | 02-07-1965 (PS) . | 02-26-1957 (PS) . | 03-06-1963 (NS) . | 03-14-1956 (S) . | 03-17-1963 (PS) . | 04-16-1961 (PS) . | 05-18-1963 (S) . | 06-06-1964 (PS) . | 06-21-1956 (PS) . | 10-04-1962 (PS) . | 10-18-1967 (PS) . | 12-06-1962 (S) . | Karabalcik . | Kurtleravsari . | Inderesi . | Gormel . | Myntdu . | Br. No.253 . | ||
. | R2 w.r.t. overall rank → . | 0.86 . | 0.75 . | 0.10 . | 0.80 . | 0.83 . | 0.87 . | 0.82 . | 0.82 . | 0.53 . | 0.79 . | 0.80 . | 0.85 . | 0.50 (PS) . | 0.49 (PS) . | 0.67 (PS) . | 0.82 (PS) . | 0.62 (PS) . | 0.76 (PS) . | ||
1 | Shifted Gompertz (2P) | 3 | 2 | 21 | 9 | 7 | 4 | 8 | 3 | 4 | 4 | 5 | 4 | 7 | 6 | 5 | 2 | 7 | 3 | 5.78 | 5 |
2 | 2-PGD (2P) | 10 | 10 | 15 | 2 | 2 | 9 | 5 | 9 | 3 | 10 | 11 | 7 | 11 | 3 | 8 | 4 | 8 | 7 | 7.44 | 7 |
3 | Dagum (3P) | 1 | 4 | 1 | 3 | 3 | 1 | 2 | 1 | 8 | 1 | 1 | 1 | 2 | 4 | 1 | 6 | 1 | 1 | 2.33 | 1 |
4 | Rice (2P) | 14 | 15 | 4 | 10 | 11 | 15 | 14 | 14 | 19 | 17 | 17 | 17 | 18 | 17 | 21 | 17 | 22 | 18 | 15.56 | 13 |
5 | Rayleigh (1P) | 30 | 30 | 30 | 30 | 30 | 30 | 30 | 30 | 18 | 30 | 27 | 30 | 17 | 18 | 23 | 25 | 19 | 23 | 26.11 | 28 |
6 | Erlang (2P) | 15 | 14 | 17 | 19 | 18 | 12 | 15 | 11 | 14 | 18 | 13 | 19 | 25 | 12 | 28 | 29 | 25 | 29 | 18.50 | 19 |
7 | Incomplete Gamma (2P) | 20 | 12 | 27 | 6 | 6 | 17 | 22 | 23 | 25 | 14 | 12 | 13 | 20 | 10 | 13 | 12 | 17 | 24 | 16.28 | 15 |
8 | Kumaraswamy (3P) | 18 | 21 | 7 | 21 | 22 | 20 | 19 | 21 | 20 | 23 | 19 | 20 | 16 | 15 | 19 | 16 | 20 | 19 | 18.67 | 22 |
9 | Log-Normal (2P) | 6 | 6 | 19 | 11 | 10 | 8 | 7 | 10 | 7 | 9 | 8 | 8 | 4 | 7 | 3 | 5 | 3 | 6 | 7.61 | 8 |
10 | Weibull (2P) | 21 | 23 | 6 | 17 | 20 | 23 | 21 | 22 | 16 | 24 | 22 | 23 | 14 | 13 | 17 | 14 | 18 | 15 | 18.28 | 17 |
11 | Fréchet (2P) | 19 | 19 | 26 | 25 | 21 | 14 | 23 | 16 | 22 | 12 | 8 | 21 | 5 | 24 | 18 | 24 | 14 | 21 | 18.44 | 18 |
12 | Nakagami (2P) | 11 | 13 | 10 | 4 | 8 | 11 | 10 | 13 | 10 | 13 | 15 | 10 | 12 | 14 | 16 | 9 | 15 | 9 | 11.28 | 10 |
13 | Beta (3P) | 12 | 20 | 8 | 18 | 15 | 18 | 19 | 20 | 11 | 20 | 16 | 21 | 13 | 8 | 14 | 10 | 13 | 11 | 14.83 | 12 |
14 | Inverse-Gamma (2P) | 8 | 7 | 23 | 16 | 12 | 7 | 12 | 5 | 12 | 2 | 3 | 11 | 1 | 21 | 10 | 11 | 9 | 10 | 10.00 | 9 |
15 | Hybrid/Nash Modified (2P) | 28 | 28 | 28 | 28 | 26 | 26 | 28 | 28 | 1 | 27 | 21 | 28 | 10 | 2 | 9 | 21 | 6 | 13 | 19.89 | 24 |
16 | 3 Para.2 Sided Power (3P) | 26 | 24 | 25 | 23 | 24 | 26 | 25 | 24 | 28 | 25 | 29 | 27 | 30 | 26 | 30 | 30 | 30 | 30 | 26.78 | 29 |
17 | Simple 3 Parameter (3P) | 29 | 29 | 29 | 29 | 29 | 29 | 29 | 29 | 24 | 29 | 30 | 29 | 21 | 11 | 27 | 28 | 23 | 27 | 26.72 | 30 |
18 | Poisson (1P) | 22 | 16 | 12 | 20 | 13 | 21 | 10 | 19 | 15 | 15 | 22 | 14 | 24 | 29 | 21 | 18 | 26 | 17 | 18.56 | 21 |
19 | Binomial (2P) | 17 | 8 | 13 | 22 | 23 | 16 | 18 | 18 | 21 | 18 | 20 | 15 | 23 | 30 | 25 | 23 | 28 | 13 | 19.50 | 23 |
20 | Negative Binominal (2P) | 23 | 22 | 14 | 12 | 16 | 22 | 12 | 14 | 13 | 11 | 24 | 12 | 28 | 26 | 15 | 15 | 12 | 20 | 17.28 | 16 |
21 | Normal (2P) | 16 | 18 | 3 | 13 | 14 | 19 | 17 | 17 | 26 | 21 | 18 | 18 | 19 | 22 | 24 | 22 | 24 | 22 | 18.50 | 19 |
22 | Gumbel Type-I (2P) | 2 | 1 | 20 | 7 | 5 | 3 | 6 | 2 | 5 | 5 | 4 | 3 | 8 | 9 | 7 | 1 | 10 | 3 | 5.61 | 4 |
23 | Logistic (2P) | 13 | 17 | 5 | 15 | 17 | 13 | 16 | 12 | 23 | 16 | 14 | 16 | 15 | 20 | 20 | 19 | 21 | 16 | 16.00 | 14 |
24 | Log-Logistic (2P) | 5 | 9 | 16 | 5 | 4 | 2 | 2 | 4 | 9 | 6 | 2 | 2 | 3 | 5 | 2 | 7 | 2 | 2 | 4.83 | 3 |
25 | Log-Gamma (2P) | 7 | 11 | 22 | 14 | 19 | 10 | 9 | 6 | 17 | 3 | 6 | 9 | 27 | 19 | 12 | 12 | 11 | 12 | 12.56 | 11 |
26 | Generalized Gamma (3P) | 9 | 5 | 2 | 1 | 1 | 6 | 1 | 8 | 2 | 8 | 10 | 5 | 9 | 1 | 4 | 3 | 4 | 5 | 4.67 | 2 |
27 | Inverse Gaussian (2P) | 4 | 3 | 18 | 8 | 9 | 5 | 4 | 7 | 6 | 7 | 7 | 6 | 6 | 16 | 6 | 8 | 5 | 8 | 7.39 | 6 |
28 | Gumbel Type-II (2P) | 24 | 25 | 9 | 27 | 25 | 24 | 24 | 25 | 29 | 26 | 25 | 24 | 22 | 23 | 26 | 26 | 27 | 26 | 24.28 | 25 |
29 | Gompertz (2P) | 25 | 26 | 11 | 26 | 27 | 25 | 25 | 26 | 30 | 28 | 26 | 26 | 26 | 25 | 29 | 27 | 29 | 28 | 25.83 | 27 |
30 | Chi-square (1P) | 27 | 26 | 24 | 24 | 28 | 28 | 25 | 27 | 27 | 22 | 28 | 25 | 29 | 28 | 11 | 20 | 16 | 25 | 24.44 | 26 |
Ranking of PDFs based on agreement with observed UH peak and time to peak
No . | PDF (No. of parameters) . | Average of event-wise qp-based ranks (Rqp) . | Overall rank (RqpO) . | Average of event-wise tp-based ranks (Rtp) . | Overall rank (RtpO) . |
---|---|---|---|---|---|
1 | Shifted Gompertz (2P) | 9.89 | 8 | 1.39 | 4 |
2 | 2-PGD (2P) | 14.39 | 13 | 1.44 | 5 |
3 | Dagum (3P) | 6.17 | 3 | 1.44 | 5 |
4 | Rice (2P) | 15.33 | 15 | 1.78 | 11 |
5 | Rayleigh (1P) | 27.00 | 29 | 2.00 | 14 |
6 | Erlang (2P) | 19.61 | 21 | 1.67 | 9 |
7 | Incomplete Gamma (2P) | 2.06 | 1 | 1.06 | 1 |
8 | Kumaraswamy (3P) | 19.83 | 22 | 1.78 | 11 |
9 | Log-Normal (2P) | 10.83 | 11 | 1.56 | 7 |
10 | Weibull (2P) | 20.89 | 24 | 1.67 | 9 |
11 | Fréchet (2P) | 7.39 | 6 | 1.72 | 10 |
12 | Nakagami (2P) | 15.50 | 17 | 1.61 | 8 |
13 | Beta (3P) | 16.11 | 18 | 1.61 | 8 |
14 | Inverse-Gamma (2P) | 6.83 | 5 | 1.50 | 6 |
15 | Hybrid/Nash Modified (2P) | 25.39 | 27 | 1.94 | 13 |
16 | 3 Para. 2 Sided Power (3P) | 25.22 | 26 | 1.11 | 2 |
17 | Simple 3 Parameter (3P) | 28.28 | 30 | 1.67 | 9 |
18 | Poisson (1P) | 19.33 | 20 | 1.50 | 6 |
19 | Binomial (2P) | 15.06 | 14 | 1.78 | 11 |
20 | Negative Binominal (2P) | 22.28 | 25 | 1.61 | 8 |
21 | Normal (2P) | 15.44 | 16 | 1.89 | 12 |
22 | Gumbel Type-I (2P) | 9.83 | 7 | 1.33 | 3 |
23 | Logistic (2P) | 10.11 | 9 | 1.67 | 9 |
24 | Log-Logistic (2P) | 5.94 | 2 | 1.44 | 5 |
25 | Log-Gamma (2P) | 6.78 | 4 | 1.50 | 6 |
26 | Generalized Gamma (3P) | 14.28 | 12 | 1.50 | 6 |
27 | Inverse Gaussian (2P) | 10.17 | 10 | 1.50 | 6 |
28 | Gumbel Type-II (2P) | 20.61 | 23 | 2.06 | 15 |
29 | Gompertz (2P) | 18.50 | 19 | 2.17 | 16 |
30 | Chi-square (1P) | 25.94 | 28 | 1.72 | 10 |
No . | PDF (No. of parameters) . | Average of event-wise qp-based ranks (Rqp) . | Overall rank (RqpO) . | Average of event-wise tp-based ranks (Rtp) . | Overall rank (RtpO) . |
---|---|---|---|---|---|
1 | Shifted Gompertz (2P) | 9.89 | 8 | 1.39 | 4 |
2 | 2-PGD (2P) | 14.39 | 13 | 1.44 | 5 |
3 | Dagum (3P) | 6.17 | 3 | 1.44 | 5 |
4 | Rice (2P) | 15.33 | 15 | 1.78 | 11 |
5 | Rayleigh (1P) | 27.00 | 29 | 2.00 | 14 |
6 | Erlang (2P) | 19.61 | 21 | 1.67 | 9 |
7 | Incomplete Gamma (2P) | 2.06 | 1 | 1.06 | 1 |
8 | Kumaraswamy (3P) | 19.83 | 22 | 1.78 | 11 |
9 | Log-Normal (2P) | 10.83 | 11 | 1.56 | 7 |
10 | Weibull (2P) | 20.89 | 24 | 1.67 | 9 |
11 | Fréchet (2P) | 7.39 | 6 | 1.72 | 10 |
12 | Nakagami (2P) | 15.50 | 17 | 1.61 | 8 |
13 | Beta (3P) | 16.11 | 18 | 1.61 | 8 |
14 | Inverse-Gamma (2P) | 6.83 | 5 | 1.50 | 6 |
15 | Hybrid/Nash Modified (2P) | 25.39 | 27 | 1.94 | 13 |
16 | 3 Para. 2 Sided Power (3P) | 25.22 | 26 | 1.11 | 2 |
17 | Simple 3 Parameter (3P) | 28.28 | 30 | 1.67 | 9 |
18 | Poisson (1P) | 19.33 | 20 | 1.50 | 6 |
19 | Binomial (2P) | 15.06 | 14 | 1.78 | 11 |
20 | Negative Binominal (2P) | 22.28 | 25 | 1.61 | 8 |
21 | Normal (2P) | 15.44 | 16 | 1.89 | 12 |
22 | Gumbel Type-I (2P) | 9.83 | 7 | 1.33 | 3 |
23 | Logistic (2P) | 10.11 | 9 | 1.67 | 9 |
24 | Log-Logistic (2P) | 5.94 | 2 | 1.44 | 5 |
25 | Log-Gamma (2P) | 6.78 | 4 | 1.50 | 6 |
26 | Generalized Gamma (3P) | 14.28 | 12 | 1.50 | 6 |
27 | Inverse Gaussian (2P) | 10.17 | 10 | 1.50 | 6 |
28 | Gumbel Type-II (2P) | 20.61 | 23 | 2.06 | 15 |
29 | Gompertz (2P) | 18.50 | 19 | 2.17 | 16 |
30 | Chi-square (1P) | 25.94 | 28 | 1.72 | 10 |
RESULTS AND DISCUSSION
Eighteen UHs were employed, each characterized by a unique time step (Δt) and properties such as qp, tp, qv, and shape, including symmetrical/normal (3 UHs), positively skewed (14 UHs), and negatively skewed (01 UH). Utilizing an objective function ‘Minimization of RMSE’ within an Excel Solver routine, optimal PDF parameters were efficiently determined, facilitating the accurate reproduction of actual UHs.
Demonstration of PDF-based UHs against Bixler Run basin UH (6 December 1962)
A symmetrical, equally fragmented actual UH (tp = 10 h) was distributed over tb (=20 h) (Figure 5(b)). Both event-specific-RS and overall-RSO showed a good agreement with R2 = 0.85 (Table 3). Some PDFs shared the same rank (Rtp) due to identical RE(tp), resulting in RtpO ranging from 1 to 16 (Table 4). This pattern was similarly observed in shape prediction, as reflected by the assignment of identical rank, RSO = 19 (Table 3). The discussion here focuses on the top five ranked (RS) PDFs (Figure 5(b)).
Dagum PDF demonstrates superiority as a UH over other PDFs, supported by its indices [RMSE = 0.0048, NSE = 98.9%, R2 = 0.990, MSE = 0.00, MAE = 0.004, MAPE = 24.09%, SE = 0.005] (Table 5). Additionally, it achieves an RSO of 1 (Table 3). Despite its overall performance, Dagum lags in estimating qp (RqpO = 3) and tp (RtpO = 5, as of 2-PGD and Log-Logistic) (Table 4). In the specific event (Table 5), Dagum UH precisely captures tp [RE(tp) = 0, Rtp = 1] but underestimates qp [RE(qp) = 5.17%, PBIAS = 2.9%, Rqp = 4] due to overestimation at the inflection point and tail end, aiming to conserve unit volume (qv = 0.967) over a given tb (Figure 5(b)).
Parameter estimation, statistical evaluation and ranking of the distributions for Bixler Run basin event (12 Jun 1962)
PDF . | RMSE . | NSE (%) . | R2 . | MSE . | MAE . | PBIAS (%) . | MAPE (%) . | SE . | Vol. (C) . | RANK shape (RS) . | Parameters . | Abs RE (qp) (%) . | Abs RE (tp) (%) . | RANK peak (Rqp) . | RANK time to peak (Rtp) . | |||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Shifted Gompertz | 0.0057 | 98.4 | 0.985 | 0.000 | 0.005 | 2.1 | 29.14 | 0.0059 | 0.997 | 4 | η,b | 0.01 | 30.82 | 6.99 | 0.00 | 10 | 1 | |
2-PGD | 0.0062 | 98.1 | 0.982 | 0.000 | 0.005 | 0.6 | 25.36 | 0.0063 | 0.994 | 7 | α,B | 11.63 | 0.94 | 7.77 | 0.00 | 15 | 1 | |
Dagum | 0.0048 | 98.9 | 0.990 | 0.000 | 0.004 | 2.9 | 24.09 | 0.0050 | 0.967 | 1 | α,b,p | 0.01 | 5.28 | 10.4 | 5.17 | 0.00 | 4 | 1 |
Rice | 0.0091 | 95.9 | 0.960 | 0.000 | 0.007 | 0.1 | 26.62 | 0.0093 | 0.999 | 17 | ν,σ | 9.89 | 3.19 | 9.41 | 0.00 | 18 | 1 | |
Rayleigh | 0.0309 | 52.4 | 0.654 | 0.001 | 0.024 | 6.2 | 135.80 | 0.0309 | 0.938 | 30 | σ | 0.03 | 50.22 | 10.00 | 30 | 2 | ||
Erlang | 0.0093 | 95.7 | 0.960 | 0.000 | 0.007 | 0.4 | 30.64 | 0.0096 | 0.996 | 19 | κ,γ | 0.01 | 11.00 | 0.44 | 0.00 | 2 | 1 | |
Incomplete-Gamma | 0.0076 | 97.1 | 0.984 | 0.000 | 0.006 | 0.1 | 31.60 | 0.0078 | 0.999 | 13 | α | 12.50 | 0.00 | 0.00 | 1 | 1 | ||
Kumaraswamy | 0.0129 | 91.7 | 0.921 | 0.002 | 0.000 | 0.0 | 41.32 | 0.0132 | 1.000 | 20 | p,q,b | 0.01 | 3.51 | 6.8 | 11.38 | 10.00 | 22 | 2 |
Log-Normal | 0.0058 | 98.3 | 0.985 | 0.002 | 0.005 | 1.4 | 28.35 | 0.0059 | 0.980 | 8 | μ,σ | 0.01 | 2.37 | 7.03 | 0.00 | 12 | 1 | |
Weibull | 0.0108 | 94.2 | 0.942 | 0.002 | 0.009 | 0.0 | 34.29 | 0.0111 | 1.000 | 23 | α,b | 1.32 | 0.01 | 12.58 | 10.00 | 24 | 2 | |
Fréchet | 0.0101 | 94.9 | 0.957 | 0.002 | 0.008 | 7.7 | 40.06 | 0.0103 | 0.923 | 21 | α,c | 0.01 | 1.24 | 5.39 | 10.00 | 6 | 2 | |
Nakagami | 0.0075 | 97.2 | 0.973 | 0.000 | 0.006 | 0.2 | 26.09 | 0.0077 | 0.998 | 10 | m,Ω | 3.09 | 123.34 | 8.72 | 0.00 | 16 | 1 | |
Beta | 0.0110 | 94.0 | 0.945 | 0.000 | 0.009 | 0.0 | 36.51 | 0.0113 | 1.000 | 21 | q,b,p | 5.51 | 20.00 | 10.6 | 10.14 | 0.00 | 19 | 1 |
Inverse-Gamma | 0.0065 | 97.9 | 0.981 | 0.000 | 0.005 | 2.7 | 31.84 | 0.0066 | 0.973 | 11 | α,B | 11.24 | 117.56 | 6.57 | 0.00 | 8 | 1 | |
Hybrid/Nash Modified | 0.0268 | 64.3 | 0.758 | 0.001 | 0.020 | 8.2 | 72.59 | 0.0275 | 0.918 | 28 | k1,k2 | 2.93 | 2.93 | 45.65 | 10.00 | 28 | 2 | |
3-Para.2-Sided Power | 0.0196 | 80.9 | 0.875 | 0.000 | 0.014 | 0.0 | 92.01 | 0.0201 | 1.000 | 27 | n,b,m | 2.00 | 20.00 | 10.0 | 28.81 | 0.00 | 26 | 1 |
Simple 3-Parameter | 0.0275 | 62.4 | 0.707 | 0.001 | 0.021 | 0.2 | 127.65 | 0.0282 | 0.998 | 29 | M,N,td | 1.37 | 1.37 | 20.0 | 46.37 | 0.00 | 29 | 1 |
Poisson | 0.0079 | 96.9 | 0.969 | 0.000 | 0.006 | 0.3 | 24.03 | 0.0079 | 0.997 | 14 | μ | 10.59 | 12.42 | 0.00 | 23 | 1 | ||
Binomial | 0.0088 | 96.1 | 0.972 | 0.000 | 0.007 | 0.1 | 29.26 | 0.0090 | 0.999 | 15 | n,p | 50.00 | 0.21 | 2.14 | 0.00 | 3 | 1 | |
Negative Binominal | 0.0079 | 96.9 | 0.969 | 0.000 | 0.006 | 0.4 | 22.56 | 0.0081 | 0.996 | 12 | r,p | 1000.00 | 0.99 | 13.27 | 0.00 | 25 | 1 | |
Normal | 0.0093 | 95.7 | 0.958 | 0.000 | 0.008 | 0.1 | 26.18 | 0.0095 | 0.999 | 18 | μ,σ | 10.39 | 3.11 | 9.39 | 0.00 | 17 | 1 | |
Gumbel Type-I | 0.0057 | 98.4 | 0.985 | 0.000 | 0.005 | 2.1 | 29.10 | 0.0059 | 0.979 | 3 | B,μ | 2.80 | 9.69 | 7.00 | 0.00 | 11 | 1 | |
Logistic | 0.0082 | 96.6 | 0.966 | 0.000 | 0.007 | 0.8 | 29.72 | 0.0084 | 0.992 | 16 | s,μ | 1.89 | 10.31 | 6.49 | 0.00 | 7 | 1 | |
Log-Logistic | 0.0048 | 98.8 | 0.989 | 0.000 | 0.004 | 2.7 | 22.88 | 0.0049 | 0.973 | 2 | α,B | 10.62 | 5.47 | 5.22 | 0.00 | 5 | 1 | |
Log-Gamma | 0.0064 | 98.0 | 0.982 | 0.000 | 0.005 | 2.4 | 31.49 | 0.0065 | 0.976 | 9 | α,B | 62.44 | 0.04 | 6.66 | 0.00 | 9 | 1 | |
Generalized Gamma | 0.0059 | 98.3 | 0.984 | 0.000 | 0.005 | 0.8 | 26.37 | 0.0062 | 0.942 | 5 | d,α,p | 17.18 | 0.08 | 0.7 | 7.43 | 0.00 | 14 | 1 |
Inverse Gaussian | 0.0060 | 98.2 | 0.984 | 0.000 | 0.005 | 1.3 | 28.70 | 0.0061 | 0.980 | 6 | γ,μ | 122.41 | 11.17 | 7.10 | 0.00 | 13 | 1 | |
Gumbel Type-II | 0.0170 | 85.6 | 0.858 | 0.000 | 0.015 | 2.0 | 86.83 | 0.0174 | 0.999 | 24 | B,α | 2.95 | 11.04 | 11.33 | 10.00 | 21 | 2 | |
Gompertz | 0.0171 | 85.5 | 0.856 | 0.000 | 0.015 | −0.4 | 91.29 | 0.0175 | 0.987 | 26 | η,b | 0.02 | 0.03 | 10.88 | 10.00 | 20 | 2 | |
Chi-Square | 0.0196 | 80.9 | 0.888 | 0.000 | 0.014 | 5.8 | 42.99 | 0.0196 | 1.004 | 25 | k | 12.00 | 37.55 | 0.00 | 27 | 1 |
PDF . | RMSE . | NSE (%) . | R2 . | MSE . | MAE . | PBIAS (%) . | MAPE (%) . | SE . | Vol. (C) . | RANK shape (RS) . | Parameters . | Abs RE (qp) (%) . | Abs RE (tp) (%) . | RANK peak (Rqp) . | RANK time to peak (Rtp) . | |||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Shifted Gompertz | 0.0057 | 98.4 | 0.985 | 0.000 | 0.005 | 2.1 | 29.14 | 0.0059 | 0.997 | 4 | η,b | 0.01 | 30.82 | 6.99 | 0.00 | 10 | 1 | |
2-PGD | 0.0062 | 98.1 | 0.982 | 0.000 | 0.005 | 0.6 | 25.36 | 0.0063 | 0.994 | 7 | α,B | 11.63 | 0.94 | 7.77 | 0.00 | 15 | 1 | |
Dagum | 0.0048 | 98.9 | 0.990 | 0.000 | 0.004 | 2.9 | 24.09 | 0.0050 | 0.967 | 1 | α,b,p | 0.01 | 5.28 | 10.4 | 5.17 | 0.00 | 4 | 1 |
Rice | 0.0091 | 95.9 | 0.960 | 0.000 | 0.007 | 0.1 | 26.62 | 0.0093 | 0.999 | 17 | ν,σ | 9.89 | 3.19 | 9.41 | 0.00 | 18 | 1 | |
Rayleigh | 0.0309 | 52.4 | 0.654 | 0.001 | 0.024 | 6.2 | 135.80 | 0.0309 | 0.938 | 30 | σ | 0.03 | 50.22 | 10.00 | 30 | 2 | ||
Erlang | 0.0093 | 95.7 | 0.960 | 0.000 | 0.007 | 0.4 | 30.64 | 0.0096 | 0.996 | 19 | κ,γ | 0.01 | 11.00 | 0.44 | 0.00 | 2 | 1 | |
Incomplete-Gamma | 0.0076 | 97.1 | 0.984 | 0.000 | 0.006 | 0.1 | 31.60 | 0.0078 | 0.999 | 13 | α | 12.50 | 0.00 | 0.00 | 1 | 1 | ||
Kumaraswamy | 0.0129 | 91.7 | 0.921 | 0.002 | 0.000 | 0.0 | 41.32 | 0.0132 | 1.000 | 20 | p,q,b | 0.01 | 3.51 | 6.8 | 11.38 | 10.00 | 22 | 2 |
Log-Normal | 0.0058 | 98.3 | 0.985 | 0.002 | 0.005 | 1.4 | 28.35 | 0.0059 | 0.980 | 8 | μ,σ | 0.01 | 2.37 | 7.03 | 0.00 | 12 | 1 | |
Weibull | 0.0108 | 94.2 | 0.942 | 0.002 | 0.009 | 0.0 | 34.29 | 0.0111 | 1.000 | 23 | α,b | 1.32 | 0.01 | 12.58 | 10.00 | 24 | 2 | |
Fréchet | 0.0101 | 94.9 | 0.957 | 0.002 | 0.008 | 7.7 | 40.06 | 0.0103 | 0.923 | 21 | α,c | 0.01 | 1.24 | 5.39 | 10.00 | 6 | 2 | |
Nakagami | 0.0075 | 97.2 | 0.973 | 0.000 | 0.006 | 0.2 | 26.09 | 0.0077 | 0.998 | 10 | m,Ω | 3.09 | 123.34 | 8.72 | 0.00 | 16 | 1 | |
Beta | 0.0110 | 94.0 | 0.945 | 0.000 | 0.009 | 0.0 | 36.51 | 0.0113 | 1.000 | 21 | q,b,p | 5.51 | 20.00 | 10.6 | 10.14 | 0.00 | 19 | 1 |
Inverse-Gamma | 0.0065 | 97.9 | 0.981 | 0.000 | 0.005 | 2.7 | 31.84 | 0.0066 | 0.973 | 11 | α,B | 11.24 | 117.56 | 6.57 | 0.00 | 8 | 1 | |
Hybrid/Nash Modified | 0.0268 | 64.3 | 0.758 | 0.001 | 0.020 | 8.2 | 72.59 | 0.0275 | 0.918 | 28 | k1,k2 | 2.93 | 2.93 | 45.65 | 10.00 | 28 | 2 | |
3-Para.2-Sided Power | 0.0196 | 80.9 | 0.875 | 0.000 | 0.014 | 0.0 | 92.01 | 0.0201 | 1.000 | 27 | n,b,m | 2.00 | 20.00 | 10.0 | 28.81 | 0.00 | 26 | 1 |
Simple 3-Parameter | 0.0275 | 62.4 | 0.707 | 0.001 | 0.021 | 0.2 | 127.65 | 0.0282 | 0.998 | 29 | M,N,td | 1.37 | 1.37 | 20.0 | 46.37 | 0.00 | 29 | 1 |
Poisson | 0.0079 | 96.9 | 0.969 | 0.000 | 0.006 | 0.3 | 24.03 | 0.0079 | 0.997 | 14 | μ | 10.59 | 12.42 | 0.00 | 23 | 1 | ||
Binomial | 0.0088 | 96.1 | 0.972 | 0.000 | 0.007 | 0.1 | 29.26 | 0.0090 | 0.999 | 15 | n,p | 50.00 | 0.21 | 2.14 | 0.00 | 3 | 1 | |
Negative Binominal | 0.0079 | 96.9 | 0.969 | 0.000 | 0.006 | 0.4 | 22.56 | 0.0081 | 0.996 | 12 | r,p | 1000.00 | 0.99 | 13.27 | 0.00 | 25 | 1 | |
Normal | 0.0093 | 95.7 | 0.958 | 0.000 | 0.008 | 0.1 | 26.18 | 0.0095 | 0.999 | 18 | μ,σ | 10.39 | 3.11 | 9.39 | 0.00 | 17 | 1 | |
Gumbel Type-I | 0.0057 | 98.4 | 0.985 | 0.000 | 0.005 | 2.1 | 29.10 | 0.0059 | 0.979 | 3 | B,μ | 2.80 | 9.69 | 7.00 | 0.00 | 11 | 1 | |
Logistic | 0.0082 | 96.6 | 0.966 | 0.000 | 0.007 | 0.8 | 29.72 | 0.0084 | 0.992 | 16 | s,μ | 1.89 | 10.31 | 6.49 | 0.00 | 7 | 1 | |
Log-Logistic | 0.0048 | 98.8 | 0.989 | 0.000 | 0.004 | 2.7 | 22.88 | 0.0049 | 0.973 | 2 | α,B | 10.62 | 5.47 | 5.22 | 0.00 | 5 | 1 | |
Log-Gamma | 0.0064 | 98.0 | 0.982 | 0.000 | 0.005 | 2.4 | 31.49 | 0.0065 | 0.976 | 9 | α,B | 62.44 | 0.04 | 6.66 | 0.00 | 9 | 1 | |
Generalized Gamma | 0.0059 | 98.3 | 0.984 | 0.000 | 0.005 | 0.8 | 26.37 | 0.0062 | 0.942 | 5 | d,α,p | 17.18 | 0.08 | 0.7 | 7.43 | 0.00 | 14 | 1 |
Inverse Gaussian | 0.0060 | 98.2 | 0.984 | 0.000 | 0.005 | 1.3 | 28.70 | 0.0061 | 0.980 | 6 | γ,μ | 122.41 | 11.17 | 7.10 | 0.00 | 13 | 1 | |
Gumbel Type-II | 0.0170 | 85.6 | 0.858 | 0.000 | 0.015 | 2.0 | 86.83 | 0.0174 | 0.999 | 24 | B,α | 2.95 | 11.04 | 11.33 | 10.00 | 21 | 2 | |
Gompertz | 0.0171 | 85.5 | 0.856 | 0.000 | 0.015 | −0.4 | 91.29 | 0.0175 | 0.987 | 26 | η,b | 0.02 | 0.03 | 10.88 | 10.00 | 20 | 2 | |
Chi-Square | 0.0196 | 80.9 | 0.888 | 0.000 | 0.014 | 5.8 | 42.99 | 0.0196 | 1.004 | 25 | k | 12.00 | 37.55 | 0.00 | 27 | 1 |
Following Dagum, the Log-Logistic PDF maintains a comparable qv = 0.973, achieving RS = 2, Rqp = 5 (RE(qp) = 5.22%) and Rtp = 1 (Table 5). However, Log-Logistic UH exhibits an earlier rise than Dagum UH at the head end, underestimates at the peak (PBIAS = 2.7%), overestimates near the point of inflection, and lies below Dagum and the actual UH at the tail end (Figure 5(b)). Notably, Log-Logistic attains RSO = 3, RqpO = 2 and RtpO = 5 (Tables 3 and 4).
Underestimating Gumbel Type-I UH at the head and peak results in RS = 3, Rqp = 11 (RE(qp) = 7%), Rtp = 1, and PBIAS = 2.1% (Table 5). Despite the peak underestimation, overestimation at the inflection point and tail end aids volume conservation (qv = 0.979) (Figure 5(b)) achieving RSO = 4, RqpO = 7, and RtpO = 3 (Tables 3 and 4).
Having RS = 4, Rqp = 10 (RE(qp) = 6.66%), and Rtp = 1, the Shifted Gompertz UH follows Gumbel Type-I UH (Table 5). The broader derived UH encompasses the actual UH, preserving unit volume despite the underestimated peak (PBIAS = 2.1%, Figure 5(b)). Overall, the ranks for Shifted Gompertz are RSO = 5, RqpO = 8, and RtpO = 4 (Tables 3 and 4).
The broad crest of the Generalized Gamma UH reduces the peak (Rqp = 14, RE(qp) = 7.43%) precisely at tp (Rtp = 1) during a given tb, resulting in minimal volumetric shortfall (qv = 0.941) (Table 5). This leads to overestimation before and after the crest, along with underestimation at the tail end, indicated by RS = 5 and PBIAS = 0.8% (Figure 5(b), Table 5). The Generalized Gamma achieves RSO = 2, RqpO = 12, and RtpO = 6 (Tables 3 and 4).
Furthermore, the Inverse Gaussian, 2-PGD, Log-Normal, Log-Gamma, Nagakami, and Inverse-Gamma, rank 6th to 11th, exhibiting consistent shape-based ranking, RSO range (Tables 3 and 5). Subsequent PDFs are ranked based on their respective indices (Tables 3–5). In the case of the positively skewed Bixler Run basin UH (16 April 1961) with tb = 14 h and tp = 6 h, RS and RSO demonstrate strong agreement with R2 = 0.87 (Table 3). However, for the negatively skewed UH (6 March 1963) with tb = 24 h and tp = 14 h, RS and RSO show poor agreement (R2 = 0.10, Table 3).
General discussion
Among the 12 new PDFs, 5 achieved top 5 rankings, RSO is discussed (Table 3). The Dagum PDF secured the 1st rank, RSO for 10 UHs across three regions while obtaining 2nd to 8th ranks for the remaining 08 UHs. Its lowest (i.e. 8th) rank for positively skewed UH (tb = 14 h and tp = 03 h, 21 June 1956) is justified by reasonable indices [RMSE = 0.0086, NSE = 98.3%, R2 = 0.985, MSE = 0.00, MAE = 0.007, MAPE = 61.56%, SE = 0.009, RE(qp) = 5.71%, RE(tp) = 0.00, PBIAS = 2.8% – excess underestimation] compared to the event-specific top-ranked Hybrid–Nash (H-N) PDF (indices: 0.0072, 98.8%, 0.988, 0.00, 0.005, 26.70%, 0.007, 8.09%, 33.33%, 0.2% – negligible underestimation). Despite H-N PDFs erroneous qp − tp estimation, it exhibits volumetric efficiency (qv ≈ 1). Notably, for this specific event (21 June 1956), RS shows poor agreement (R2 = 0.53) with RSO (Table 3).
RSO = 2 confirms Generalized Gamma's second rank, trailing Dagum, with ranks (RS) ranging from 1st (over 4 UHs) to 5th (RS) across 12 out of 18 UHs (Table 3). Notably, its highest (i.e. 10th) rank (RS) against a positively skewed UH (tb = 14 h and tp = 04 h, 18 October 1967) exhibits indices [RMSE = 0.016, NSE = 95.5%, R2 = 0.963, MSE = 0.00, MAE = 0.013, MAPE = 40.35%, SE = 0.018, RE(qp) = 6.51%, RE(tp) = 0.00]. In comparison to the top-ranked (1st) Dagum with respective indices [0.010, 98.1%, 0.984, 0.00, 0.008, 24.90%, 0.011, 3.49%, 0.00], this deviation is attributed to an attenuated peak expanding UH, resulting in a slight deviation (PBIAS = 0.03%). Excluding this exceptional event, Generalized Gamma demonstrates adaptability in reproducing actual qp and tp, as indicated by RqpO = 12 and RtpO = 6 (Table 4).
The Log-Logistic PDF achieves the second rank (RS) across seven events, displaying exceptional performance in replicating negatively skewed UH (6 March 1963) [RS-16, RMSE = 0.009, NSE = 94.5%, R2 = 0.946, MSE = 0.00, MAE = 0.007, MAPE = 54.14%, SE = 0.009, RE(qp) = 2.58%, RE(tp) = 14.29%]. For other events, its rank ranges up to 9th (Table 3). Despite having the same RE(tp) as the superior Dagum PDF [RS-1, RMSE = 0.003, NSE = 99.4%, R2 = 0.994, MSE = 0.00, MAE = 0.002, MAPE = 10.54%, SE = 0.003, RE(qp) = 1.06%, RE(tp) = 14.29%], the Log-Logistic PDF's smooth peak occurring prior to the actual tapered/negatively skewed peak is attributed to UH deviations (PBIAS = 3.2%). Excluding this event, the Log-logistic PDF justifiably estimates qp (Rqp = 2) and tp (Rtp = 5, equivalent to Dagum and 2-PGD) (Table 4). Consequently, it demonstrates superiority over Dagum and Generalized Gamma in qp − tp estimation.
Gumbel Type-I demonstrates performance comparable to Log-Logistic in handling negatively skewed UH [6 March 1963, RS-20, RMSE = 0.011, NSE = 91.9%, R2 = 0.920, MSE = 0.00, MAE = 0.009, MAPE = 64.55%, SE = 0.011, RE(qp) = 6.59%, RE(tp) = 14.29%]. While it attains the 1st rank (RS) for 2 events, its rank for the remaining 15 events ranges up to 10th (Table 3). The underestimated peak expands both limbs of the UH, resulting in a PBIAS of 2.7%. In total, Gumbel Type-I secures the 7th rank for qp and the 3rd rank for tp (Table 4).
Shifted Gompertz exhibits variability in RS, ranging from 2nd (for 2 events) to 9th (Table 3), closely mirroring the Gumbel Type-I PDF with a slight difference [RS-21, RMSE = 0.011, NSE = 91.9%, R2 = 0.920, MSE = 0.00, MAE = 0.009, MAPE = 64.63%, SE = 0.011, RE(qp) = 6.55%, RE(tp) = 14.29%] over a negatively skewed exceptional UH (6 March 1963). In qp − tp fitting, it aligns with Gumbel Type-I, securing the 8th rank for qp and the 4th rank for tp overall (Table 4).
The Inverse Gaussian performs satisfactorily, securing ranks from 3rd to 8th (RS) over 16 events. However, for two specific events (Kurtleravsari UH and Bixler Run Basin UH-6 March 1963), it attains the 16th (RMSE = 0.002, NSE = 97%) and 18th (RMSE = 0.010, NSE = 92.3%) ranks, respectively (Table 3). Notably, Generalized Gamma and Dagum excel in these particular events, achieving the 1st rank (RS = 1). The 16th rank is attributed to deviations in the positively skewed UH limbs (PBIAS = 58.79%) from the actual UH with an oscillatory crest. The 18th rank results from underestimating a smooth symmetric peak that represents the actual UH with a tapered crest. In overall qp and tp estimation, the Inverse Gaussian secures the 10th and 6th ranks, respectively (Table 4).
With RSO = 7, the 2-PGD (Table 3) appears relatively less efficient [RS-15, RMSE = 0.008, NSE = 95.2%, R2 = 0.953, MSE = 0.00, MAE = 0.007, MAPE = 46.50%, SE = 0.008, RE(qp) = 3.67%, RE(tp) = 14.29%] compared to Dagum (RS-1) in the specified event (6 March 1963). This is attributed to peak attenuation, causing UH expansion and deviations/errors (PBIAS = 0.08%). However, its volumetric accuracy is higher (qv = 0.99) than Dagum (qv = 0.98). Excluding this event, 2-PGD demonstrates its capability across other events, with RS ranging from 2nd to 11th (Table 3). It shares the same RtpO (5th) as Dagum and Log-Logistic while following Generalized Gamma with RqpO = 13th (Table 4).
The Log-Normal PDF, with an RSO of 8th (Table 3), holds ranks, RS from 3rd to 11th, except an exceptional RS = 19th [RMSE = 0.010, NSE = 92.7%, R2 = 0.928, MSE = 0.001, MAE = 0.001, MAPE = 59.20%, SE = 0.010, RE in qp = 3.37%, RE in tp = 21.43%] against the UH (6 March 1963). This deviation is attributed to the early occurrence of an underestimated peak (PBIAS = 1.7%). In the overall assessment, Log-Normal is ranked 11th for qp and 7th for tp (Table 4).
Similarly, the Inverse-Gamma, with an RSO of 9th, exhibits RS ranging from 1st to 23rd (6 March 1963) across 18 events (Table 3). In the overall assessment of qp − tp, it secures the 5th and 6th ranks, respectively (Table 4). The Japanese Nagakami PDF holds the 10th place in RSO, with event-wise ranks (RS) varying from 4th to 16th (Inderesi Creek basin UH) across 18 events (Table 3). However, in qp − tp estimation, it attains the 17th and 8th ranks, respectively (Table 4). The new Log-Gamma PDF achieves an RSO of 11th, with RS varying from 3rd to 27th (Karabalcik Creek basin UH) across 18 events (Table 3). Excluding this, it efficiently estimates qp following Dagum with an RqpO of 4th (Table 4). Its performance aligns with Inverse-Gamma, Poisson, Log-Gamma, Generalized Gamma, and Inverse Gaussian, sharing an RtpO of 6th (Table 4)
The Negative Binomial, Chi-square, Poisson, and Binomial PDFs exhibit unacceptable performance with negative NSE and the highest RE(qp), resulting in RS of 26, 28, 29, and 30, respectively. This is evident in defining a positively skewed UH (tp = 17 h, Kurtleravsari Creek basin) with the longest tb (=110 h) among the 18 observed UHs. The computed peaks, occurring either prior to or later than the actual peak, are two to three times in magnitude, aiming to conserve unit volume over a constrained tb (≈40 h). In contrast, by underestimating the peak [RE(qp) = 44.91%] at the exact tp over an identical tb as observed, the 3-Parameter 2-Sided Power-based triangular UH attains identical but unacceptable rank, RS = 26 as the Negative Binomial, while conserving the approximate volume (qv = 0.980) with a positive NSE (=54.8%).
The novel Negative Binomial, Poisson, and Binomial PDFs hold ranks (RSO) of 16th, 21st, and 23rd, respectively (Table 3). Chi-square (RSO = 26th) and 3-Parameter 2-Sided Power (RSO = 29th, attributed to triangular UH fitted) follow them. The Simple 3-Parameter PDF secures the last (30th) rank based on shape and qp (Tables 3 and 4), owing to its underestimated concave peaks against 10 out of 18 observed UHs. In terms of tp, the Gompertz PDF attains the last (16th) rank among all PDFs (Table 4). When tb is computed, qv is either underestimated or overestimated, justifiable by positive or negative PBIAS, respectively. A PBIAS = 0 indicates unit volume conservation over the given tb. Furthermore, the best PDFs for a specific watershed class (based on drainage area) have been identified (Table 6) by averaging the event-wise PDF ranks for that particular watershed class.
Drainage area based watershed class-wise best shape-defining PDFs
S. N. . | Watershed class and events used . | Top 10 UH shape-defining PDFs . | |
---|---|---|---|
1 | Karabalcik Bixler Run 10.6–38.85 km2 (13) |
|
|
2 | Inderesi Bridge No. 253 Gormel Ermenek 98–114.22–141.5 km2 (3) |
|
|
3 | Myntdu-Leska 350 km2 (1) |
|
|
4 | Kurtleravsari Aksu 3360 km2 (1) |
|
|
S. N. . | Watershed class and events used . | Top 10 UH shape-defining PDFs . | |
---|---|---|---|
1 | Karabalcik Bixler Run 10.6–38.85 km2 (13) |
|
|
2 | Inderesi Bridge No. 253 Gormel Ermenek 98–114.22–141.5 km2 (3) |
|
|
3 | Myntdu-Leska 350 km2 (1) |
|
|
4 | Kurtleravsari Aksu 3360 km2 (1) |
|
|
CONCLUSIONS
We examined both established hydrologic literature PDFs and those not yet tested as UH to identify a reliable PDF. However, the PDFs resembling UH are not restricted to those investigated in this study. There is a significant possibility that additional PDFs may exist in the statistical literature, but they require evaluation as UH before being considered. The following conclusions can be derived from the study:
The diverse selection of available PDFs empowers practitioners to choose the most suitable one for their objectives. While PDFs may resemble UH, a single PDF is insufficient to precisely capture all the characteristics of a UH. Different UHs within a region may be well-fitted by more than one PDF.
The parameters obtained through the Excel-based GRG-NLP algorithm define smooth PDF-based UHs, automatically addressing volume and non-negativity constraints. The proposed Maple codes enhance and expedite hydrologic simulation.
The novel PDFs (Dagum, Generalized Gamma, Log-Logistic, Gumbel Type-I, and Shifted Gompertz) accurately replicate UH compared to established PDFs (2-PGD, Chi-Square, Weibull, Beta, Hybrid Model, Simple 3 Parameter, Rayleigh, etc.) recommended in hydrologic literature.
Among the 18 UHs, 16 have been excellently defined (RSO = 1) by new PDFs (10 by Dagum, 04 by Generalized Gamma, and 02 by Gumbel Type-I) while the remaining 2 were well-fitted by the Hybrid/Modified-Nash model.
Dagum, Generalized Gamma, and Log-Logistic PDFs excel in preserving the rising and recession limbs across all watershed classes. The Incomplete-Gamma PDF is particularly effective in predicting both qp and tp. However, the widely-used 2-PGD falls short in accurately estimating qp while it reasonably estimates tp.
Gompertz, among the novel PDFs, is the sole PDF that fits erratic UHs, such as those modeled by the existing PDFs (Rayleigh – 10 events and Simple 3-Parameter – 7 events), both of which rank last across 17 events.
The novel 3-parameter PDFs, Dagum and Generalized Gamma, can be utilized with known ‘qp and tp’ by fixing the least sensitive shape parameter as unity. In contrast, for the existing 3-parameter Kumaraswamy, Beta, 3-Parameter 2-Sided Power, and Simple 3-Parameter PDFs, the application requires knowledge of tb.
The proposed Dagum, Generalized Gamma, Log-Logistic, Gumbel Type-I, and Shifted Gompertz PDFs are highly reliable for small-size watersheds. For mid-size watersheds, Dagum, Log-Logistic, Log-Normal, Generalized Gamma, and Inverse Gaussian PDFs are most suitable. Sequentially, Generalized Gamma, Hybrid/Nash Modified, 2-PGD, Dagum, and Log-Logistic PDFs are efficient for large-size watersheds.
ACKNOWLEDGEMENTS
The authors express their gratitude to IIT, Roorkee and NIH, Roorkee for furnishing the essential research facilities.
DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.
CONFLICT OF INTEREST
The authors declare there is no conflict.