Stochastic Rational Method for Estimation of Flood Peak Uncertainty in Arid Basins: Comparison between Monte Carlo and First Order Second Moment Methods with a Case Study in Southwest Saudi Arabia
Abstract
:1. Introduction
2. Study Area
3. Methodology
3.1. Rational Method
- A = the basin area;
- C = the runoff coefficient;
- i = the intensity of the rainfall; and
- Q = the peak flood.
3.2. Joint Probability Distribution of the Rational Method Parameters
- X = the p-vector of variables;
- = the vector of the mean of the variables;
- = the covariance matrix of the variables;
- p = the number of variables; and
- T = transpose operation of the matrix.
- = the mean of the variable Xj and j = 1, 2, and 3;
- = the variance of the variable Xj;
- = the covariance between variables Xj and Xk where k = 1, 2, and 3.
3.3. Uncertainty Quantification
- = variance of the logarithms of the peak discharge;
- = variance of the logarithms of the runoff coefficient;
- = variance of the logarithms of the rainfall intensity;
- = variance of the logarithms of the area;
- = covariance of the logarithms of runoff coefficient and the logarithms of the rainfall intensity;
- = covariance of the logarithms of runoff coefficient and the logarithms of the area; and
- = covariance of the logarithms of the area and the logarithms of the rainfall intensity.
- = variance of the logarithms of the runoff volume;
- = variance of the logarithms of the rainfall depth;
- = covariance of the logarithms of runoff coefficient and the logarithms of the rainfall depth; and
- = covariance of the logarithms of the area and the logarithms of the rainfall depth.
3.4. Uncertainty Analysis
4. Results and Discussion
4.1. Comparison between the Data and the Generated Realizations of the Rational Method Input Parameters
4.2. Comparison between the Data and the Generated Realizations of the Rational Method Output Variables
4.3. Comparison between the Correlation of the Parameters in the Rational Method and the Generated Realizations
4.4. Comparison between Probability Distributions of the Peak Flow, Runoff Volumes, and the Realizations
4.5. Comparison between the Quantiles of the Data and the Realizations
4.6. Comparison between Statistics of the Data and the Statistics of FOSM and MC
5. Summary and Conclusions
- The correlation coefficient between the rational method parameters (C, A, i, and R) was relatively weak; it also showed a negative correlation except between C and i. The correlation coefficient between A and i was the strongest (0.46), while that between C and R was the weakest.
- Although the correlation between the parameters was weak, the model was capable of simulating the rational model parameters and estimating the Qp and V reasonably well. The reason is that the relations for the generation process do not depend only on the correlations, but also depend on the mean and variance of the parameters.
- The log-normal distribution fit both the data and the generated peak flow well using the rainfall intensity generated from the time of concentration equation developed by some researchers in the literature. Furthermore, the log-normal distribution fit both the data and the generated runoff volumes well. Therefore, the log-normal distribution can be used as a model for the generation of the peak flow and runoff volume in ungagged basins.
- The use of the FOSM method underestimates the data in comparison with the MC method. Therefore, we recommend using the MC method with an intensity calculated with the tc equation [32], since the CV = 1.6, which was the closest to the data (CV = 1.4). However, in terms of runoff volume, FOSM and MC provided a CV = 2.1 and 1.9, respectively, which were in the order of magnitude of the CV of the data (2.4). Therefore, both methods provide an acceptable estimation from a practical point of view.
- The generated realizations fell within the confidence levels, except for a few marginal cases, which are expected due to the long tail of the log-normal distribution, and consequently, an extreme event may occur.
- The model can be used to generate peak flows and the associated confidence limits in ungagged basins from the statistics of the rainfall, basin area, and runoff coefficient based on the equations developed in this study.
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Acknowledgments
Conflicts of Interest
References
- Almazroui, M. Simulation of present and future climate of Saudi Arabia using a regional climate model (PRECIS). Int. J. Clim. 2013, 33, 2247–2259. [Google Scholar] [CrossRef]
- Farquharson, F.; Meigh, J.; Sutcliffe, J. Regional flood frequency analysis in arid and semi-arid areas. J. Hydrol. 1992, 138, 487–501. [Google Scholar] [CrossRef]
- Cordery, I.; Fraser, J. Theme 4-6-Some Hydrological Characteristics of the Australian Arid Zone; National Institute of Hydrology: Roorkee, India, 2000. [Google Scholar]
- Nemec, J.; Rodier, J.A. Streamflow characteristics in areas of low precipitation (with special reference to low and high flows). In The Hydrology of Areas of Low Precipitation; International Association of Hydrological Sciences: Paris, France, 1979; pp. 125–140. [Google Scholar]
- McMahon, T. Hydrological characteristics of arid zones. In The Hydrology of Areas of Low Precipitation; International Association of Hydrological Sciences: Paris, France, 1979; pp. 105–123. [Google Scholar]
- Pilgrim, D.H.; Chapman, T.G.; Doran, D.G. Problems of rainfall-runoff modelling in arid and semiarid regions. Hydrol. Sci. J. 1988, 33, 379–400. [Google Scholar] [CrossRef]
- Majone, B.; Avesani, D.; Zulian, P.; Fiori, A.; Bellin, A. Analysis of high streamflow extremes in climate change studies: How do we calibrate hydrological models? Hydrol. Earth Syst. Sci. 2022, 26, 3863–3883. [Google Scholar] [CrossRef]
- Bedient, P.B.; Huber, W.; Vieux, B. Hydrology and Floodplain Analysis; Pearson: New York, NY, USA, 2008; Volume 816. [Google Scholar]
- Hua, J.; Liang, Z.; Yu, Z. A modified rational formula for flood design in small basins 1. JAWRA J. Am. Water Resour. Assoc. 2003, 39, 1017–1025. [Google Scholar] [CrossRef]
- Wang, S.; Wang, H. Extending the Rational Method for assessing and developing sustainable urban drainage systems. Water Res. 2018, 144, 112–125. [Google Scholar] [CrossRef]
- Pilgrim, D.; Cordery, I. Flood runoff. In Handbook of Hydrology; Chapter 9; Maidment, D., Ed.; McGraw-Hill, Inc.: New York, NY, USA, 1993. [Google Scholar]
- Titmarsh, G.W.; Cordery, I.; Pilgrim, D.H. Calibration Procedures for Rational and USSCS Design Flood Methods. J. Hydraul. Eng. 1995, 121, 61–70. [Google Scholar] [CrossRef]
- Chin, D.A. Estimating Peak Runoff Rates Using the Rational Method. J. Irrig. Drain. Eng. 2019, 145, 04019006. [Google Scholar] [CrossRef]
- Al-Amri, N.S.; Ewea, H.A.; Elfeki, A.M. Revisit the rational method for flood estimation in the Saudi arid environment. Arab. J. Geosci. 2022, 15, 532. [Google Scholar] [CrossRef]
- Dhakal, N.; Fang, X.; Cleveland, T.G.; Thompson, D.B. Revisiting modified rational method. In World Environmental and Water Resources Congress 2011: Bearing Knowledge for Sustainability; American Society of Civil Engineers: Reston, WV, USA, 2011; pp. 751–762. [Google Scholar]
- Young, C.B.; McEnroe, B.M.; Rome, A.C. Empirical Determination of Rational Method Runoff Coefficients. J. Hydrol. Eng. 2009, 14, 1283–1289. [Google Scholar] [CrossRef]
- Efstratiadis, A.; Koussis, A.D.; Koutsoyiannis, D.; Mamassis, N. Flood design recipes vs. reality: Can predictions for ungauged basins be trusted? Nat. Hazards Earth Syst. Sci. 2014, 14, 1417–1428. [Google Scholar] [CrossRef] [Green Version]
- Kottek, M.; Grieser, J.; Beck, C.; Rudolf, B.; Rubel, F. World map of the Köppen-Geiger climate classification updated. Meteorol. Z. 2006, 15, 259–263. [Google Scholar] [CrossRef] [PubMed]
- Sharon, D. The spottiness of rainfall in a desert area. J. Hydrol. 1972, 17, 161–175. [Google Scholar] [CrossRef]
- Dames, M. Representative Basins Study for Wadi: Yiba, Habwnah, Tabalah, Liyyah and Al-Lith (Main Report) Kingdom of Saudi Arabia, Ministry of Agriculture and Water; Water Resource Development Department: Riyadh, Saudi Arabia, 1988. [Google Scholar]
- Ewea, H.A.; Elfeki, A.; Al-Amri, N. Development of intensity–duration–frequency curves for the Kingdom of Saudi Arabia. Geomat. Nat. Hazards Risk 2017, 8, 570–584. [Google Scholar] [CrossRef] [Green Version]
- Ewea, H.A.; Al-Amri, N.S.; Elfeki, A.M. Analysis of maximum flood records in the arid environment of Saudi Arabia. Geomat. Nat. Hazards Risk 2020, 11, 1743–1759. [Google Scholar] [CrossRef]
- Abdulrazzak, M.J.; Sorman, A.U. Transmission Losses from Ephemeral Stream in Arid Region. J. Irrig. Drain. Eng. 1994, 120, 669–675. [Google Scholar] [CrossRef]
- Elfeki, A.M.M.; Ewea, H.A.R.; Bahrawi, J.A.; Al-Amri, N.S. Incorporating transmission losses in flash flood routing in ephemeral streams by using the three-parameter Muskingum method. Arab. J. Geosci. 2015, 8, 5153–5165. [Google Scholar] [CrossRef]
- Farran, M.M.; Elfeki, A.; Elhag, M.; Chaabani, A. A comparative study of the estimation methods for NRCS curve number of natural arid basins and the impact on flash flood predications. Arab. J. Geosci. 2021, 14, 121. [Google Scholar] [CrossRef]
- Al-Shanti, A.; El-Mahdy, O.; Hassan, M.; Hussein, A. A Comparative Study of Five Volcanic-Hosted Sulfide Mineralizations in the Arabian Shield. Earth Sci. 1993, 6, 1–33. [Google Scholar] [CrossRef]
- Wheater, H.; Butler, A.; Stewart, E.; Hamilton, G. A multivariate spatial-temporal model of rainfall in southwest Saudi Arabia. I. Spatial rainfall characteristics and model formulation. J. Hydrol. 1991, 125, 175–199. [Google Scholar] [CrossRef]
- Eagleson, P.S.; Fennessey, N.M.; Qinliang, W.; Rodriguez-Iturbe, I. Application of spatial Poisson models to air mass thunderstorm rainfall. J. Geophys. Res. Atmos. 1987, 92, 9661. [Google Scholar] [CrossRef]
- Wheater, H.; Onof, C.; Butler, A.; Hamilton, G. A multivariate spatial-temporal model of rainfall in southwest Saudi Arabia. II. Regional analysis and long-term performance. J. Hydrol. 1991, 125, 201–220. [Google Scholar] [CrossRef]
- Chow, V.; Maidment, D.; Mays, L. Applied Hydrology; International Editions; McGraw-Hill: New York, NY, USA, 1988. [Google Scholar]
- Mood, A.M. Introduction to the Theory of Statistics; McGraw-Hill: New York, NY, USA, 1950. [Google Scholar]
- Albishi, M.; Bahrawi, J.; Elfeki, A. Empirical equations for flood analysis in arid zones: The Ari-Zo model. Arab. J. Geosci. 2017, 10, 51. [Google Scholar] [CrossRef]
Parameter | Mean | SD | CV | Mean [ln()] | SD [ln()] | CV [ln()] |
---|---|---|---|---|---|---|
Area (m2) | 1.23 × 109 | 9.58 × 108 | 0.78 | 20.64 | 0.80 | 0.04 |
Runoff Coefficient, C | 0.07 | 0.07 | 0.99 | −3.18 | 1.09 | −0.35 |
Rainfall depth, R (m) | 0.0150 | 0.0150 | 2.55 | −4.54 | 0.85 | 0.37 |
Average rainfall intensity, I (m/s) | 1.05 × 10−6 | 1.2 × 10−6 | 1.22 | −14.15 | 0.83 | −0.06 |
Peak discharge, Qp (m3/s) | 84.08 | 117.08 | 2.65 | 3.79 | 1.16 | 0.32 |
Runoff volume, V (m3) | 8.04 × 105 | 1.92 × 106 | 2.08 | 12.92 | 1.20 | 0.10 |
ρ(lnC, lni)= | 0.12 | |||||
ρ(lni, lnA)= | −0.46 | |||||
ρ(lnC, lnA)= | −0.30 | |||||
ρ(lnC, lnR)= | −0.10 | |||||
ρ(lnR, lnA)= | −0.29 |
Item | Data | FOSM | MC (Intensity Based on tc Calculated by Albishi et al. Equation (22)) | MC (Based on Max Intensity within the Storm) | MC (Based on Average Intensity) |
---|---|---|---|---|---|
Mean Q (m3/s) | 84.1 | 59.3 | 118.5 | 120.8 | 75.1 |
SD (Q) (m3/s) | 117.1 | 133.2 | 186.2 | 203.8 | 137.2 |
CV (Q) | 1.4 | 2.2 | 1.6 | 1.7 | 1.8 |
Mean V (m3) | 804,009.0 | 409,373.5 | NA | NA | 974,432.4 |
SD (V) (m3) | 1,920,805.4 | 843,894.8 | NA | NA | 1,837,235.7 |
CV (V) | 2.4 | 2.1 | NA | NA | 1.9 |
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. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Al-Amri, N.S.; Ewea, H.A.; Elfeki, A.M. Stochastic Rational Method for Estimation of Flood Peak Uncertainty in Arid Basins: Comparison between Monte Carlo and First Order Second Moment Methods with a Case Study in Southwest Saudi Arabia. Sustainability 2023, 15, 4719. https://doi.org/10.3390/su15064719
Al-Amri NS, Ewea HA, Elfeki AM. Stochastic Rational Method for Estimation of Flood Peak Uncertainty in Arid Basins: Comparison between Monte Carlo and First Order Second Moment Methods with a Case Study in Southwest Saudi Arabia. Sustainability. 2023; 15(6):4719. https://doi.org/10.3390/su15064719
Chicago/Turabian StyleAl-Amri, Nassir S., Hatem A. Ewea, and Amro M. Elfeki. 2023. "Stochastic Rational Method for Estimation of Flood Peak Uncertainty in Arid Basins: Comparison between Monte Carlo and First Order Second Moment Methods with a Case Study in Southwest Saudi Arabia" Sustainability 15, no. 6: 4719. https://doi.org/10.3390/su15064719
APA StyleAl-Amri, N. S., Ewea, H. A., & Elfeki, A. M. (2023). Stochastic Rational Method for Estimation of Flood Peak Uncertainty in Arid Basins: Comparison between Monte Carlo and First Order Second Moment Methods with a Case Study in Southwest Saudi Arabia. Sustainability, 15(6), 4719. https://doi.org/10.3390/su15064719