This study examines all of the equations necessary to derive the parameters for seven probability distributions of three parameters typically used in flood frequency research, namely the Pearson III (PE3), the generalized extreme value (GEV), the Weibull (W3), the log-normal (LN3), the generalized Pareto Type II (PG), the Rayleigh (RY) and the log-logistic (LL3) distributions, using the higher-order linear moments method (LH-moments). The analysis represents the expansion of previous research whose results were presented in previous materials, and is part of hydrological research aimed at developing a standard for calculating maximum flows based on L-moments and LH-moments. The given methods for calculating the parameters of the examined distributions are used to calculate the maximum flows on Romania’s Prigor River. For both methods, the criterion for selecting the most suitable distribution is represented by the diagram of the L-skewness–L-kurtosis and LH-skewness–LH-kurtosis. The results for Prigor River show that the PG distribution is the best model for the L-moments method, the theoretical values of the statistical indicators being 0.399 and 0.221. The RY distribution is the best model for the LH-moments technique, with values of 0.398 and 0.192 for the two statistical indicators.

Considering the wide applicability of two-parameter distributions in the frequency analysis of extreme events, this article presents new elements regarding the use of thirteen two-parameter probability distributions, using three parameter estimation methods. All the necessary elements for the application of these distributions are presented using the method of ordinary moments (MOM), the method of linear moments (L-moments) and the method of high order linear moments (LH-moments). Only these three methods are analyzed, because they are usually applied in the analysis regarding the regionalization of extreme events. As a case study, the frequency analysis of the maximum annual flows on the Siret River, Lungoci station, Romania, was made. For the recommended methods (L- and LH-moments), from the thirteen analyzed distributions, the log-normal distribution had the best results, with the theoretical values L-coefficient of variation and L-kurtosis (0.297, 0.192, 0.323, 0.185, 0.336, and 0.185) best approximating the corresponding values of the recorded data (0.339, 0.185, 0.233, 0.199, 0.198, and 0.205).

A direct way to estimate the likelihood and magnitude of extreme events is frequency analysis. This analysis is based on historical data and assumptions of stationarity, and is carried out with the help of probability distributions and different methods of estimating their parameters. Thus, this article presents all the relations necessary to estimate the parameters with the LH-moments method for the family of distributions defined only by the quantile function, namely, the Wakeby distribution of 4 and 5 parameters, the Lambda distribution of 4 and 5 parameters, and the Davis distribution. The LH-moments method is a method commonly used in flood frequency analysis, and it uses the annual series of maximum flows. The frequency characteristics of the two analyzed methods, which are both involved in expressing the distributions used in the first two linear moments, as well as in determining the confidence interval, are presented. The performances of the analyzed distributions and the two presented methods are verified in the following maximum flows, with the Bahna river used as a case study. The results are presented in comparison with the L-moments method. Following the results obtained, the Wakeby and Lambda distributions have the best performances, and the LH-skewness and LH-kurtosis statistical indicators best model the indicators’ values of the sample (0.5769, 0.3781, 0.548 and 0.3451). Similar to the L-moments method, this represents the main selection criterion of the best fit distribution.