Applicability of Wakeby distribution in flood frequency analysis: a case study for eastern Australia

Parametric method of flood frequency analysis (FFA) involves fitting of a probability distribution to the observed flood data at the site of interest. When record length at a given site is relatively longer and flood data exhibits skewness, a distribution having more than three parameters is often used in FFA such as log‐Pearson type 3 distribution. This paper examines the suitability of a five‐parameter Wakeby distribution for the annual maximum flood data in eastern Australia. We adopt a Monte Carlo simulation technique to select an appropriate plotting position formula and to derive a probability plot correlation coefficient (PPCC) test statistic for Wakeby distribution. The Weibull plotting position formula has been found to be the most appropriate for the Wakeby distribution. Regression equations for the PPCC tests statistics associated with the Wakeby distribution for different levels of significance have been derived. Furthermore, a power study to estimate the rejection rate associated with the derived PPCC test statistics has been undertaken. Finally, an application using annual maximum flood series data from 91 catchments in eastern Australia has been presented. Results show that the developed regression equations can be used with a high degree of confidence to test whether the Wakeby distribution fits the annual maximum flood series data at a given station. The methodology developed in this paper can be adapted to other probability distributions and to other study areas. Copyright © 2014 John Wiley & Sons, Ltd.     

Model averaging versus model selection: estimating design floods with uncertain river flow data

ABSTRACT

This study compares model averaging and model selection methods to estimate design floods, while accounting for the observation error that is typically associated with annual maximum flow data. Model selection refers to methods where a single distribution function is chosen based on prior knowledge or by means of selection criteria. Model averaging refers to methods where the results of multiple distribution functions are combined. Numerical experiments were carried out by generating synthetic data using the Wakeby distribution function as the parent distribution. For this study, comparisons were made in terms of relative error and root mean square error (RMSE) referring to the 1-in-100 year flood. The experiments show that model averaging and model selection methods lead to similar results, especially when short samples are drawn from a highly asymmetric parent. Also, taking an arithmetic average of all design flood estimates gives estimated variances similar to those obtained with more complex weighted model averaging.     

Timescale-based flood typing to estimate temporal changes in flood frequencies

ABSTRACT

The flood peak is the dominating characteristic in nearly all flood-statistical analyses. Contrary to the general assumptions of design flood estimation, the peak is not closely related to other flood characteristics. Differentiation of floods into types provides a more realistic view. Often different parts of the probability distribution function of annual flood peaks are dominated by different flood types, which raises the question how shifts in flood regimes would modify the statistics of annual maxima. To answer this, a distinction into five flood types is proposed; then, temporal changes in flood-type frequencies are investigated. We show that the frequency of floods caused by heavy rain has increased significantly in recent years. A statistical model is developed that simulates peaks for each event type by type-specific peak–volume relationships. In a simulation study, we show how changes in frequency of flood event type lead to changes in the quantiles of annual maximum series.     

A new seasonal design flood method based on bivariate joint distribution of flood magnitude and date of occurrence

Abstract

Seasonal design floods which consider information on seasonal variation are very important for reservoir operation and management. The seasonal design flood method currently used in China is based on seasonal maximum (SM) samples and assumes that the seasonal design frequency is equal to the annual design frequency. Since the return period associated with annual maximum floods is taken as the standard in China, the current seasonal design flood cannot satisfy flood prevention standards. A new seasonal design flood method, which considers dates of flood occurrence and magnitudes of the peaks (runoff), was proposed and established based on copula function. The mixed von Mises distribution was selected as marginal distribution of flood occurrence dates. The Pearson Type III and exponential distributions were selected as the marginal distribution of flood magnitude for annual maximum flood series and peak-over-threshold samples, respectively. The proposed method was applied at the Geheyan Reservoir, China, and then compared with the currently used seasonal design flood methods. The case study results show that the proposed method can satisfy the flood prevention standard, and provide more information about the flood occurrence probabilities in each sub-season. The results of economic analysis show that the proposed design flood method can enhance the floodwater utilization rate and give economic benefits without lowering the annual flood protection standard.

Citation Chen, L., Guo, S. L., Yan, B. W., Liu, P. & Fang, B. (2010) A new seasonal design flood method based on bivariate joint distribution of flood magnitude and date of occurrence. Hydrol. Sci. J. 55(8), 1264–1280.     

The use of the Akaike Information Criterion in the identification of an optimum flood frequency model

Abstract

For a long time now, the hydrologist has been faced with the problem of finding which of the many possible probability distribution functions can be used most effectively in flood frequency analyses. This problem has been mainly due to the insufficiency of the conventional goodness-of-fit procedures when used with the typically skewed flood probability distributions. In this study, the Akaike Information Criterion (AIC) goodness-of-fit test is used to identify more objectively the optimum model for flood frequency analysis in Kenya from a class of competing models. The class is comprised of (a) seven three-parameter density functions, namely, log-normal, Pearson, log-Pearson, Fisher-Tippet, log-Fisher-Tippet, Walter Boughton and log-Walter Boughton; and (b) two five-parameter density functions, namely, Wakeby and log-Wakeby. The AIC is also used in this study as a method of testing for the existence of outlier peak-flow values in the peak annual data used. A modified version of the chi-square goodness-of-fit test is also used, but only for the sake of comparison with the AIC.     

Hydrological Sciences Journal Contents 2002

Abstract

Flood frequency estimation is crucial in both engineering practice and hydrological research. Regional analysis of flood peak discharges is used for more accurate estimates of flood quantiles in ungauged or poorly gauged catchments. This is based on the identification of homogeneous zones, where the probability distribution of annual maximum peak flows is invariant, except for a scale factor represented by an index flood. The numerous applications of this method have highlighted obtaining accurate estimates of index flood as a critical step, especially in ungauged or poorly gauged sections, where direct estimation by sample mean of annual flood series (AFS) is not possible, or inaccurate. Therein indirect methods have to be used. Most indirect methods are based upon empirical relationships that link index flood to hydrological, climatological and morphological catchment characteristics, developed by means of multi-regression analysis, or simplified lumped representation of rainfall–runoff processes. The limits of these approaches are increasingly evident as the size and spatial variability of the catchment increases. In these cases, the use of a spatially-distributed, physically-based hydrological model, and time continuous simulation of discharge can improve estimation of the index flood. This work presents an application of the FEST-WB model for the reconstruction of 29 years of hourly streamflows for an Alpine snow-fed catchment in northern Italy, to be used for index flood estimation. To extend the length of the simulated discharge time series, meteorological forcings given by daily precipitation and temperature at ground automatic weather stations are disaggregated hourly, and then fed to FEST-WB. The accuracy of the method in estimating index flood depending upon length of the simulated series is discussed, and suggestions for use of the methodology provided.

Editor D. Koutsoyiannis     

Flood frequency analysis of Turkey using L‐moments method

The index flood procedure coupled with the L‐moments method is applied to the annual flood peaks data taken at all stream‐gauging stations in Turkey having at least 15‐year‐long records. First, screening of the data is done based on the discordancy measure (D_i) in terms of the L‐moments. Homogeneity of the total geographical area of Turkey is tested using the L‐moments based heterogeneity measure, H, computed on 500 simulations generated using the four parameter Kappa distribution. The L‐moments analysis of the recorded annual flood peaks data at 543 gauged sites indicates that Turkey as a whole is hydrologically heterogeneous, and 45 of 543 gauged sites are discordant which are discarded from further analyses. The catchment areas of these 543 sites vary from 9·9 to 75121 km² and their mean annual peak floods vary from 1·72 to 3739·5 m³ s^?1. The probability distributions used in the analyses, whose parameters are computed by the L‐moments method are the general extreme values (GEV), generalized logistic (GLO), generalized normal (GNO), Pearson type III (PE3), generalized Pareto (GPA), and five‐parameter Wakeby (WAK). Based on the L‐moment ratio diagrams and the |Z_dist|‐statistic criteria, the GEV distribution is identified as the robust distribution for the study area (498 gauged sites). Hence, for estimation of flood magnitudes of various return periods in Turkey, a regional flood frequency relationship is developed using the GEV distribution. Next, the quantiles computed at all of 543 gauged sites by the GEV and the Wakeby distributions are compared with the observed values of the same probability based on two criteria, mean absolute relative error and determination coefficient. Results of these comparisons indicate that both distributions of GEV and Wakeby, whose parameters are computed by the L‐moments method, are adequate in predicting quantile estimates. Copyright © 2011 John Wiley & Sons, Ltd.     

Maps of flood statistics for regional flood frequency analysis in New Zealand

Abstract

The two-parameter EV1 distribution adequately describes New Zealand's flood series. Contour maps of [Qbar]/A^0.8 and Q₁₀₀[Qbar] are presented, where [Qbar] is the mean annual flood, A is the basin area and Q₁₀₀ is the 1% annual exceedance probability flood. The maps are based directly on measured discharge series from a large sample of river recording stations. Thus when basins are ungauged, or have just a short record, an estimate of a design flood Q_T with specified annual exceedance probability (1/T) can be obtained using map estimates of [Qbar]/A^0.8 and Q₁₀₀[Qbar], without having first to estimate rainfall statistics for the basin, a particularly difficult task in sparsely instrumented mountainous areas. These maps succinctly summarize a great deal of hydrological information and permit improved flood frequency estimates.     

A distribution-free ordinal classification of floods based on moments

ABSTRACT

Classification of floods is often based on return periods of their peaks estimated from probability distributions and hence depends on assumptions. The choice of an appropriate distribution function and parameter estimation are often connected with high uncertainties. In addition, limited length of data series and the stochastic characteristic of the occurrence of extreme events add further uncertainty. Here, a distribution-free classification approach is proposed based on statistical moments. By using robust estimators the sampling effects are reduced and time series of different lengths can be analysed together. With a developed optimization procedure, locally and regionally consistent flood categories can be defined. In application, it is shown that the resulting flood categories can be used to assess the spatial extent of extreme floods and their coincidences. Moreover, groups of gauges, where simultaneous events belong to the same classes, are indicators for homogeneous groups of gauges in regionalization.     

LH-moment estimation of Wakeby distribution with hydrological applications

This article discusses the method of higher-order L-moment (LH-moment) estimation for the Wakeby distribution (WAD), and describes and formulates details of parameter estimation using LH-moments for WAD. Monte Carlo simulation is performed, to illustrate the performance of the LH-moment method via heavy-tail quantiles (over all quantiles) using WAD. The LH-moment method proves as useful and effective as the L-moment approach in handling data that follow WAD, and it is then applied to annual maximum flood and wave height data.     

Homogeneous region delineation based on annual flood generation mechanisms

Abstract

Flood distributions can have unimodal or multimodal densities due to different flood generation mechanisms such as snowmelt and rainfall in the annual flood series. When applying nonparametric frequency analysis to annual flood data from the province of New Brunswick in Canada, unimodal, bimodal and heavy-tailed distribution shapes were found. By grouping basins with similarly-shaped densities on a geographical basis, homogeneous regions were delineated. Regional equations derived for a homogeneous region gave lower integral square errors than those of province-wide equations.     

Robust flood statistics: comparison of peak over threshold approaches based on monthly maxima and TL-moments

ABSTRACT

Flood quantile estimation based on partial duration series (peak over threshold, POT) represents a noteworthy alternative to the classical annual maximum approach since it enlarges the available information spectrum. Here the POT approach is discussed with reference to its benefits in increasing the robustness of flood quantile estimations. The classical POT approach is based on a Poisson distribution for the annual number of exceedences, although this can be questionable in some cases. Therefore, the Poisson distribution is compared with two other distributions (binomial and Gumbel-Schelling). The results show that only rarely is there a difference from the Poisson distribution. In the second part we investigate the robustness of flood quantiles derived from different approaches in the sense of their temporal stability against the occurrence of extreme events. Besides the classical approach using annual maxima series (AMS) with the generalized extreme value distribution and different parameter estimation methods, two different applications of POT are tested. Both are based on monthly maxima above a threshold, but one also uses trimmed L-moments (TL-moments). It is shown how quantile estimations based on this “robust” POT approach (rPOT) become more robust than AMS-based methods, even in the case of occasional extraordinary extreme events.

Editor M.C. Acreman Associate editor A. Viglione     

Partial L-moments for the analysis of censored flood samples / Utilisation des L-moments partiels pour l’analyse d’échantillons tronqués de crues

Abstract

Abstract A parameter estimation method is proposed for fitting the generalized extreme value (GEV) distribution to censored flood samples. Partial L-moments (PL-moments), which are variants of L-moments and analogous to ?partial probability weighted moments?, are defined for the analysis of such flood samples. Expressions are derived to calculate PL-moments directly from uncensored annual floods, and to fit the parameters of the GEV distribution using PL-moments. Results of Monte Carlo simulation study show that sampling properties of PL-moments, with censoring flood samples of up to 30% are similar to those of simple L-moments, and also that both PL-moment and LH-moments (higher-order L-moments) have similar sampling properties. Finally, simple L-moments, LH-moments, and PL-moments are used to fit the GEV distribution to 75 annual maximum flow series of Nepalese and Irish catchments, and it is found that, in some situations, both LH- and PL-moments can produce a better fit to the larger flow values than simple L-moments.     

Comparison between the peaks-over-threshold method and the annual maximum method for flood frequency analysis

Abstract

Flood frequency analysis can be made by using two types of flood peak series, i.e. the annual maximum (AM) and peaks-over-threshold (POT) series. This study presents a comparison of the results of both methods for data from the Litija 1 gauging station on the Sava River in Slovenia. Six commonly used distribution functions and three different parameter estimation techniques were considered in the AM analyses. The results showed a better performance for the method of L-moments (ML) when compared with the conventional moments and maximum likelihood estimation. The combination of the ML and the log-Pearson type 3 distribution gave the best results of all the considered AM cases. The POT method gave better results than the AM method. The binomial distribution did not offer any noticeable improvement over the Poisson distribution for modelling the annual number of exceedences above the threshold.

Editor D. Koutsoyiannis

Citation Bezak, N., Brilly, M., and ?raj, M., 2014. Comparison between the peaks-over-threshold method and the annual maximum method for flood frequency analysis. Hydrological Sciences Journal, 59 (5), 959–977.     

The impact of seasonal flood peak dependence on annual maxima design quantiles

If the maximum annual peak flow series are a mixture of summer and winter flows, a seasonal approach to flood frequency analysis is necessary. While considering seasonal maxima as mutually independent events, the annual maxima distribution is defined as the product of seasonal distributions. However, if the independency assumption does not hold, a bivariate approach with dependent margins should be applied, i.e. the copula approach. The impact of dependency on design quantiles is investigated here in the context of the Fréchet-Hoeffding inequality defining copula bounds and the definition of dependency. The results of the two approaches are compared using six catchments in the San River basin, where in four cases the dependency of seasonal maxima has been identified as positive significant and no strong dominance of any one season is observed. The product model leads to higher estimates of design quantiles than do models where the dependency is taken into account and, therefore, is safe.

EDITOR R. Woods ASSOCIATE EDITOR A. Fiori     

Flood season segmentation based on the probability change-point analysis technique

Abstract

The segmentation of flood seasons has both theoretical and practical importance in hydrological sciences and water resources management. The probability change-point analysis technique is applied to segmenting a defined flood season into a number of sub-seasons. Two alternative sampling methods, annual maximum and peaks-over-threshold, are used to construct the new flow series. The series is assumed to follow the binomial distribution and is analysed with the probability change-point analysis technique. A Monte Carlo experiment is designed to evaluate the performance of proposed flood season segmentation models. It is shown that the change-point based models for flood season segmentation can rationally partition a flood season into appropriate sub-seasons. China's new Three Gorges Reservoir, located on the upper Yangtze River, was selected as a case study since a hydrological station with observed flow data from 1882 to 2003 is located 40 km downstream of the dam. The flood season of the reservoir can be reasonably divided into three sub-seasons: the pre-flood season (1 June–2 July); the main flood season (3 July–10 September); and the post-flood season (11–30 September). The results of flood season segmentation and the characteristics of flood events are reasonable for this region.

Citation Liu, P., Guo, S., Xiong, L. & Chen, L. (2010) Flood season segmentation based on the probability change-point analysis technique. Hydrol. Sci. J. 55(4), 540–554.     

基于SWAT模型的淮河上游流域设计洪水修订

变化环境下洪水序列的一致性遭到破坏,引发基于统计原理计算的设计洪水可靠性下降,亟需开展非一致性条件下的设计洪水修订研究.以淮河上游流域为研究区域,运用Pettitt检验法和滑动t检验法综合检测年最大洪峰流量序列突变点,在此基础上,采用SWAT分布式水文模型对变异前的洪峰与洪量序列进行还现,利用径流深的模拟结果修订设计洪...     

Reassessment of flood frequency using historical information for the River Ouse at York,UK (1200–2000)

Abstract

The reassessment of flood risk at York, UK, is pertinent in light of major flooding in November 2000, and heightened concerns of a perceived increase in flooding nationally. Systematic flood level readings from 1877 and a wealth of documentary records dating back as far as 1263 AD give the City of York a long and rich history of flood records. This extended flood record provides an opportunity to reassess estimates of flood frequency over a time scale not normally possible within flood frequency analysis. This paper re-evaluates flood frequency at York, considering the strengths and weaknesses in estimates resulting from four contrasting methods of analysis and their corresponding data: (a) single-site analysis of gauged annual maxima; (b) pooled analysis of multi-site gauged annual maxima; (c) combined analysis of systematic annual maxima augmented with historical peaks, and (d) analysis of only the very largest peaks using a Generalized Pareto Distribution. Use of the historical information was found to yield risk estimates which were lower and considered to be more credible than those achieved using gauged records alone.

Citation Macdonald, N. & Black, A. R. (2010) Reassessment of flood frequency using historical information for the River Ouse at York, UK (1200–2000). Hydrol. Sci. J. 55(7), 1152–1162.     

Flood frequency analysis of southern Africa: I. Delineation of homogeneous regions

Abstract

A graphical test is presented to check if recorded annual maximum flood data for a group of gauging stations in a region belong to a common parent distribution (P). The test compares the observed at site L-coefficient of variation (Lcv) with its sampling distribution. The latter is obtained by generating synthetic sequences from an assumed parent distribution (P). A group of sites is deemed to be homogeneous if the observed Lcv, treated as an order statistic, lies within its sampling distribution. The proposed test has been applied to annual maximum flood data from Tanzania to delineate the country into 12 homogeneous regions.     

Non-identical models for seasonal flood frequency analysis

Abstract

The aim is to build a seasonal flood frequency analysis model and estimate seasonal design floods. The importance of seasonal flood frequency analysis and the advantages of considering seasonal design floods in the derivation of reservoir planning and operating rules are discussed, recognising that seasonal flood frequency models have been in use for over 30 years. A set of non-identical models with non-constant parameters is proposed and developed to describe flows that reflect seasonal flood variation. The peak-over-threshold (POT) sampling method was used, as it is considered to provide significantly more information on flood seasonality than annual maximum (AM) sampling and has better performance in flood seasonality estimation. The number of exceedences is assumed to follow the Poisson distribution (Po), while the peak exceedences are described by the exponential (Ex) and generalized Pareto (GP) distributions and a combination of both, resulting in three models, viz. Po-Ex, Po-GP and Po-Ex/GP. Their performances are analysed and compared. The Geheyan and the Baiyunshan reservoirs were chosen for the case study. The application and statistical experiment results show that each model has its merits and that the Po-Ex/GP model performs best. Use of the Po-Ex/GP model is recommended in seasonal flood frequency analysis for the purpose of deriving reservoir operation rules.