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.     

Statistical analysis of extreme events in a non-stationary context via a Bayesian framework: case study with peak-over-threshold data

Statistical analysis of extremes currently assumes that data arise from a stationary process, although such an hypothesis is not easily assessable and should therefore be considered as an uncertainty. The aim of this paper is to describe a Bayesian framework for this purpose, considering several probabilistic models (stationary, step-change and linear trend models) and four extreme values distributions (exponential, generalized Pareto, Gumbel and GEV). Prior distributions are specified by using regional prior knowledge about quantiles. Posterior distributions are used to estimate parameters, quantify the probability of models and derive a realistic frequency analysis, which takes into account estimation, distribution and stationarity uncertainties. MCMC methods are needed for this purpose, and are described in the article. Finally, an application to a POT discharge series is presented, with an analysis of both occurrence process and peak distribution.     

Appraisal of regional and index flood quantile estimators

Our results illustrate the performance of at-site and regional GEV/PWM flood quantile estimators in regions with different coefficients of variation, degrees of regional heterogeneity, record lengths, and number of sites. Analytic approximations of bias and variance are employed. For realistic GEV distributions and short records, the index-flood quantile estimator performs better than a 2-parameter GEV/PWM quantile estimator with a regional shape parameter, or a 3-parameter at-site GEV/PWM quantile estimator, in both humid and especially in arid regions, as long as the degree of regional heterogeneity is moderate. As regional heterogeneity or record lengths increases, 2-parameter estimators quickly dominate. Flood frequency models that assign probabilities larger than 2% to negative flows are unrealistic; experiments employing such distributions provide questionable results. This appraisal generally demonstrates the value of regionalizing estimators of the shape of a flood distribution, and sometimes the coefficient of variation.     

Alternative PWM-estimators of the Gumbel distribution

Probability weighted moments (PWM) are widely used in hydrology for estimating parameters of statistical distributions, including the Gumbel distribution. The classical PWM-approach considers the moments β_i=E[XFⁱ] with i=0,1 for estimation of the Gumbel scale and location parameters. However, there is no reason why these probability weights (F⁰ and F¹) should provide the most efficient PWM-estimators of Gumbel parameters and quantiles. We explore an extended class of PWMs that does not impose arbitrary restrictions on the values of i. Estimation based on the extended class of PWMs is called the generalized method of probability weighted moments (GPWM) to distinguish it from the classical procedure. In fact, our investigation demonstrates that it may be advantage to use weight functions that are not of the form Fⁱ. We propose an alternative PWM-estimator of the Gumbel distribution that maintains the computational simplicity of the classical PWM method, but provides slightly more accurate quantile estimates in terms of mean square error of estimation. A simple empirical formula for the standard error of the proposed quantile estimator is presented.     

Appraisal of regional and index flood quantile estimators

Our results illustrate the performance of at-site and regional GEV/PWM flood quantile estimators in regions with different coefficients of variation, degrees of regional heterogeneity, record lengths, and number of sites. Analytic approximations of bias and variance are employed. For realistic GEV distributions and short records, the index-flood quantile estimator performs better than a 2-parameter GEV/PWM quantile estimator with a regional shape parameter, or a 3-parameter at-site GEV/PWM quantile estimator, in both humid and especially in arid regions, as long as the degree of regional heterogeneity is moderate. As regional heterogeneity or record lengths increases, 2-parameter estimators quickly dominate. Flood frequency models that assign probabilities larger than 2% to negative flows are unrealistic; experiments employing such distributions provide questionable results. This appraisal generally demonstrates the value of regionalizing estimators of the shape of a flood distribution, and sometimes the coefficient of variation.     

Estimation of the annual maximum distribution from samples of maxima in separate seasons

Quantile estimates of the annual maximum distribution can be obtained by fitting theoretical distributions to the maxima in separate seasons, e.g. to the monthly maxima. In this paper, asymptotic expressions for the bias and the variance of such estimates are derived for the case that the seasonal maxima follow a Gumbel distribution. Results from these expressions are presented for a situation with no seasonal variation and for maximum precipitation depths at Uccle/Ukkel (Belgium). It is shown that the bias is often negligible and that the variance reduction by using seasonal maxima instead of just the annual maxima strongly depends on the seasonal variation in the data. A comparison is made between the asymptotic standard error of quantile estimates from monthlymaxima with those from a partial duration series. Much attention is paid to the effect of model misspecification on the resulting quantile estimates of the annual maximum distribution. The use of seasonal maxima should be viewed with caution when the upper tail of this distribution is of interest.     

Cautionary note on multicomponent flood distributions for annual maxima

Multicomponent probability distributions such as the two‐component Gumbel distribution are sometimes applied to annual flood maxima when individual floods are seen as belonging to different classes, depending on physical processes or time of year. However, hydrological inconsistencies may arise if only nonclassified annual maxima are available to estimate the component distribution parameters. In particular, an unconstrained best fit to annual flood maxima may yield some component distributions with a high probability of simulating floods with negative discharge. In such situations, multicomponent distributions cannot be justified as an improved approximation to a local physical reality of mixed flood types, even though a good data fit is achieved. This effect usefully illustrates that a good match to data is no guarantee against degeneracy of hydrological models. Copyright © 2016 John Wiley & Sons, Ltd.     

Probability Estimates of Extreme Temperature Events: Stochastic Modelling Approach vs. Extreme Value Distributions

The paper deals with the probability estimates of temperature extremes (annual temperature maxima and heat waves) in the Czech Republic. Two statistical methods of probability estimations are compared; one based on the stochastic modelling of time series of the daily maximum temperature (TMAX) using the first-order autoregressive (AR(1)) model, the other consisting in fitting the extreme value distribution to the sample of annual temperature peaks.The AR(1) model is able to reproduce the main characteristics of heat waves, though the estimated probabilities should be treated as upper limits because of deficiencies in simulating the temperature variability inherent to the AR(1) model. Theoretical extreme value distributions do not yield good results when applied to maximum annual lengths of heat waves and periods of tropical days (TMAX 30°C), but it is the best method for estimating the probability and recurrence time of annual one-day temperature extremes. However, there are some difficulties in the application: the use of the two-parameter Gumbel distribution and the three-parameter generalized extreme value (GEV) distribution may lead to different results, particularly for long return periods. The resulting values also depend on the chosen procedure of parameter estimation. Based on our findings, the shape parameter testing for the GEV distribution and the L moments technique for parameter estimation may be recommended.The application of the appropriate statistical tools indicates that the heat wave and particularly the long period of consecutive tropical days in 1994 were probably a more rare event than the record-breaking temperatures in July 1983 exceeding 40°C. An improvement of the probability estimate of the 1994 heat wave may be expected from a more sophisticated model of the temperature series.     

On the informative value of the largest sample element of log-Gumbel distribution

Extremes of stream flow and precipitation are commonly modeled by heavytailed distributions. While scrutinizing annual flow maxima or the peaks over threshold, the largest sample elements are quite often suspected to be low quality data, outliers or values corresponding to much longer return periods than the observation period. Since the interest is primarily in the estimation of the right tail (in the case of floods or heavy rainfalls), sensitivity of upper quantiles to largest elements of a series constitutes a problem of special concern. This study investigated the sensitivity problem using the log-Gumbel distribution by generating samples of different sizes (n) and different values of the coefficient of variation by Monte Carlo experiments. Parameters of the log-Gumbel distribution were estimated by the probability weighted moments (PWMs) method, method of moments (MOMs) and maximum likelihood method (MLM), both for complete samples and the samples deprived of their largest elements. In the latter case, the distribution censored by the non-exceedance probability threshold, F _T , was considered. Using F _T instead of the censored threshold T creates possibility of controlling estimator property. The effect of the F _T value on the performance of the quantile estimates was then examined. It is shown that right censoring of data need not reduce an accuracy of large quantile estimates if the method of PWMs or MOMs is employed. Moreover allowing bias of estimates one can get the gain in variance and in mean square error of large quantiles even if ML method is used.     

The PWM large quantile estimates of heavy tailed distributions from samples deprived of their largest element / Estimation des grands quantiles de distributions à queue à décroissance lente par la méthode des moments pondérés par les probabilités à partir d'échantillons amputés de leur plus grande valeur

Abstract

Extremes of streamflow are usually modelled using heavy tailed distributions. While scrutinising annual flow maxima or the peaks over threshold, the largest elements in a sample are often suspected to be low quality data, outliers or values corresponding to much longer return periods than the observation period. In the case of floods, since the interest is focused mainly on the estimation of the right-hand tail of a distribution function, sensitivity of large quantiles to extreme elements of a series becomes the problem of special concern. This study investigated the sensitivity problem using the log-Gumbel distribution by generating samples of different sizes and different values of the coefficient of L-variation by means of Monte Carlo experiments. Parameters of the log-Gumbel distribution were estimated by the probability weighted moments (PWM) method, both for complete samples and the samples deprived of their largest element. In the latter case Hosking's concept of the “A” type PWM with Type II censoring was employed. The largest value was censored above the random threshold T corresponding to the non-exceedence probability F _T. The effect of the F _T value on the performance of the quantile estimates was then examined. Experimental results show that omission of the largest sample element need not result in a decrease in the accuracy of large quantile estimates obtained from the log-Gumbel model by the PWM method.     

Single station and regional analysis of daily rainfall extremes

A peaks over threshold (POT) method of analysing daily rainfall values is developed using a Poisson process of occurrences and a generalised Pareto distribution (GPD) for the exceedances. The parameters of the GPD are estimated by the method of probability weighted moments (PWM) and a method of combining the individual estimates to define a regional curve is proposed.     

On seasonal approach to flood frequency modelling. Part I: Two‐component distribution revisited

The annual peak flow series of the Polish rivers are mixtures of summer and winter flows. In the Part I of a sequence of two papers, theoretical aspects of applicability of seasonal approach to flood frequency analysis (FFA) in Poland are discussed. A testing procedure is introduced for the seasonal model and the data overall fitness. Conditions for objective comparative assessment of accuracy of annual maxima (AM) and seasonal maxima (SM) approaches to FFA are formulated and finally Gumbel (EV1) distribution is chosen as seasonal distribution for detailed investigation. Sampling properties of AM quantile x(F) estimates are analysed and compared for the SM and AM models for equal seasonal variances. For this purpose, four estimation methods were used, employing both asymptotic approach and sampling experiments. Superiority of the SM over AM approach is stated evident in the upper quantile range, particularly for the case of no seasonal variation in the parameters of Gumbel distribution. In order to learn whether the standard two‐ and three‐parameter flood frequency distributions can be used to model the samples generated from the Two‐Component Extreme Value 1 (TCEV1) distribution, the shape of TCEV1 probability density function (PDF) has been tested in terms of bi‐modality. Then the use of upper quantile estimate obtained from the dominant season of extreme floods (DEFS) as AM upper quantile estimate is studied and respective systematic error is assessed. The second part of the paper deals with advantages and disadvantages of SM and AM approach when applied to real flow data of Polish rivers. Copyright © 2011 John Wiley & Sons, Ltd.     

Beta-�� distribution and its application to hydrologic events

The beta-κ distribution is a distinct case of the generalized beta distribution of the second kind. In previous studies, beta-p and beta-κ distributions have played important roles in representing extreme events, and thus, the present paper uses the beta-κ distribution. Further, this paper uses the method of moments and the method of L-moments to estimate the parameters from the beta-κ distribution, and to demonstrate the performance of the proposed model, the paper presents a simulation study using three estimation methods (including the maximum likelihood estimation method) and beta-κ and non beta-κ samples. In addition, this paper evaluates the performance of the beta-κ distribution by employing two widely used extreme value distributions (i.e., the GEV and Gumbel distributions) and two sets of actual data on extreme events.     

A nonstationary extreme value distribution for analysing the cessation of karst spring discharge

The effects of climate change and population growth in recent decades are leading us to consider their combined and potentially extreme consequences, particularly regarding hydrological processes, which can be modeled using a generalized extreme value (GEV) distribution. Most of the GEV models were based on a stationary assumption for hydrological processes, in contrast to the nonstationary reality due to climate change and human activities. In this paper, we present the nonstationary generalized extreme value (NSGEV) distribution and use it to investigate the risk of Niangziguan Springs discharge decreasing to zero. Rather than assuming the location, scale, and shape parameters to be constant as one might do for a stationary GEV distribution analysis, the NSGEV approach can reflect the dynamic processes by defining the GEV parameters as functions of time. Because most of the GEV model is designed to evaluate maxima (e.g. flooding, represented by positive numbers), and spring discharge cessation is a ?minima’, we deduced an NSGEV model for minima by applying opposite numbers, i.e. negative instead of positive numbers. The results of the model application to Niangziguan Springs showed that the probability of zero discharge at Niangziguan Springs will be 1/80 in 2025, and 1/10 in 2030. After 2025, the rate of decrease in spring discharge will accelerate, and the probability that Niangziguan Springs will cease flowing will dramatically increase. The NSGEV model is a robust method for analysing karst spring discharge. Copyright © 2013 John Wiley & Sons, Ltd.     

Depth distribution of interrill overland flow and the formation of rills

A Gumbel distribution for maxima is proposed as a model for the depths of interrill overland flow. The model is tested against three sets of field measurements of interrill overland flow depths obtained on shrubland and grassland hillslopes at Walnut Gulch Experimental Watershed, southern Arizona. The model is found to be a satisfactory fit to 81 of the 90 measured distributions. The shape δ and location λ parameters of all fitted distributions are strongly correlated with discharge. However, whereas a common relationship exists between discharge and δ for all depth distributions, the relationships with λ vary systematically downslope. Using the Gumbel distribution as a model for the distribution of overland flow depths, a probabilistic model for the initiation of rills is developed, drawing upon the previous work of Nearing. As an illustration of this approach, we apply this model to the shrubland and grassland hillslopes at Walnut Gulch. It is concluded that the presence of rills on the shrubland, but not on the grassland, is due to the greater runoff coefficient for the shrubland and/or the greater propensity of the shrubland for soil disturbance compared with the grassland. Finally, a generalized conceptual model for rill initiation is proposed. This model takes account of the depth distribution of overland flow, the probability of flow shear stress in excess of local soil shear strength, the spatial variability in soil shear strength and the diffusive effect of soil detachment by raindrops. Copyright © 2005 John Wiley & Sons, Ltd.     

Bivariate flood frequency analysis using the copula function: a case study of the Litija station on the Sava River

As an alternative to the commonly used univariate flood frequency analysis, copula frequency analysis can be used. In this study, 58 flood events at the Litija gauging station on the Sava River in Slovenia were analysed, selected based on annual maximum discharge values. Corresponding hydrograph volumes and durations were considered. Different bivariate copulas from three families were applied and compared using different statistical, graphical and upper tail dependence tests. The parameters of the copulas were estimated using the method of moments with the inversion of Kendall's tau. The Gumbel–Hougaard copula was selected as the most appropriate for the pair of peak discharge and hydrograph volume (Q‐V). The same copula was also selected for the pair hydrograph volume and duration (V‐D), and the Student‐t copula was selected for the pair of peak discharge and hydrograph duration (Q‐D). The differences among most of the applied copulas were not significant. Different primary, secondary and conditional return periods were calculated and compared, and some relationships among them were obtained. Copyright © 2014 John Wiley & Sons, Ltd.     

Simple scaling properties of Canadian annual average streamflow

The spatial scaling properties of Canadian annual average streamflow (abbreviated as AASF) are assessed using both the product moments (PMs) and the probability weighted moments (PWMs) of AASF across the entire country and in its sub-climatic regions. By the PMs, the log relationship between the kth moments of AASF and the drainage area can be almost represented by a perfect straight line across the entire country and in its sub-climatic regions, whose regression parameters are a linear function of the moment order. By the PWMs, the logarithm of the kth PWM is a linear function of the logarithm of drainage area for the entire country and its sub-climatic regions, where its slope (or scale exponent) in a region is constant and is independent of the order. These results indicate that Canadian AASF exhibits simple scaling and drainage area alone may describe most of the variability in the moments of AASF. The third approach, based on the log linearity between quantiles and drainage area, is applied to Region 2, also demonstrate simple scaling of AASF in that region, as concluded from using PMs and PWMs methods, which indicates that all three methods are consistent. The simple scaling results provide a basis for using the index flood method to conduct regional frequency analysis of AASF in Canada.     

Estimation of the GEV distribution from censored samples by method of partial probability weighted moments

The concept of partial probability weighted moments (PPWM), which can be used to estimate a distribution from censored samples, is introduced. Unbiased estimators of PPWM are derived. An application is made to estimating parameters and quantiles of the generalized extreme value (GEV) distribution from censored samples. Censored samples yield high quantile estimates which are almost as efficient as those obtained from uncensored samples. This could be a very useful technique for dealing with the undesirable effects of low outliers which occur in semiarid and arid zones.     

Frequency analysis of climate extreme events in Zanjan, Iran

In this study, generalized extreme value distribution (GEV) and generalized Pareto distribution (GPD) were fitted to the maximum and minimum temperature, maximum wind speed, and maximum precipitation series of Zanjan. Maximum (minimum) daily and absolute annual observations of Zanjan station from 1961 to 2011 were used. The parameters of the distributions were estimated using the maximum likelihood estimation method. Quantiles corresponding to 2, 5, 10, 25, 50, and 100 years return periods were calculated. It was found that both candidate distributions fitted to extreme events series, were statistically reasonable. Most of the observations from 1961 to 2011 were found to fall within 1–10 years return period. Low extremal index (θ) values were found for excess maximum and minimum temperatures over a high threshold, indicating the occurrence of consecutively high peaks. For the purpose of filtering the dependent observations to obtain a set of approximately independent threshold excesses, a declustering method was performed, which separated the excesses into clusters, then the de-clustered peaks were fitted to the GPD. In both models, values of the shape parameters of extreme precipitation and extreme wind speed were close to zero. The shape parameter was less negative in the GPD than the GEV. This leads to significantly lower return period estimates for high extremes with the GPD model.