Nonstationary Flood Frequency Analysis for Annual Flood Peak and Volume Series in Both Univariate and Bivariate Domain

Flood frequency analysis for practical application is traditionally based on the assumption of stationarity, but this assumption has been open to doubt in recent years. A number of studies have focused on the nonstationary flood frequency analysis, and the associated causes of nonstationarity. In this study, the annual maximum flood peak and flood volume of Wangkuai reservoir watershed were used, and several univariate and bivariate models were established to investigate the nonstationary flood frequency, with the distribution parameters changing over the climate indices (NPO, Niño3) and the check dam indices (CDIp, CDIv). In the univariate models, the Weibull distribution performed best and exhibited an undulate behavior for both flood peak and volume, which tended to describe the nonstationarity reasonably well. The bivariate models were constructed using copulas, of which the optimal Weibull distribution in the univariate flood frequency analysis was considered as marginal distributions within the joint distribution. The results showed that the Gumbel-Hougaard copula offered the best joint distribution, and most of the probability isolines crossed each other, which demonstrated the possibility that the occurrence of combinations of the flood peak and volume may be the same under multiple effects of phase changes in the climate patterns and certain human activities (i.e. soil and water conservation). The most likely events were elaborated in diagrams, and the associated combinations of the flood peak and volume were smaller than that estimated by the fixed parameters (i.e. stationary condition) during most of the study period, while it was the opposite in 1956, 1959 and 1963. The results highlight the necessity of nonstationary flood frequency analysis under various conditions in both univariate and multivariate domains.