Clustering stochastic point process model for flood risk analysis

Since the introduction into flood risk analysis, the partial duration series method has gained increasing acceptance as an appealing alternative to the annual maximum series method. However, when the base flow is low, there is clustering in the flood peak or flow volume point process. In this case, the general stochastic point process model is not suitable to risk analysis. Therefore, two types of models for flood risk analysis are derived on the basis of clustering stochastic point process theory in this paper. The most remarkable characteristic of these models is that the flood risk is considered directly within the time domain. The acceptability of different models are also discussed with the combination of the flood peak counted process in twenty years at Yichang station on the Yangtze river. The result shows that the two kinds of models are suitable ones for flood risk analysis, which are more flexible compared with the traditional flood risk models derived on the basis of annual maximum series method or the general stochastic point process theory. Received: September 29, 1997     

Quantifying differences between reservoir inflows and dam site floods using frequency and risk analysis methods

Reservoirs are the most important constructions for water resources management and flood control. Great concern has been paid to the effects of reservoir on downstream area and the differences between inflows and dam site floods due to the changes of upstream flow generation and concentration conditions after reservoir’s impoundment. These differences result in inconsistency between inflow quantiles and the reservoir design criteria derived by dam site flood series, which can be a potential risk and must be quantificationally evaluated. In this study, flood frequency analysis (FFA) and flood control risk analysis (FCRA) methods are used with the long reservoir inflow series derived from a multiple inputs and single output model and a copula-based inflow estimation model. The results of FFA and FCRA are compared and the influences on reservoir flood management are also discussed. The Three Gorges Reservoir (TGR) in China is selected as a case study. Results show that the differences between the TGR inflow and dam site floods are significant which result in changes on its flood control risk rates. The mean values of TGR’s annual maximum inflow peak discharge and 3 days flood volume have increased 5.58 and 3.85% than the dam site ones, while declined by 1.82 and 1.72% for the annual maximum 7 and 15 days flood volumes. The flood control risk rates of middle and small flood events are increased while extreme flood events are declined. It is shown that the TGR can satisfy the flood control task under current hydrologic regime and the results can offer references for better management of the TGR.     

洞庭湖三口洪峰流量和水位变异特性分析

近几十年来,受荆江裁弯、葛洲坝工程运用、三峡水库拦蓄调度、洞庭湖治理以及长江上游水土保持措施等因素的综合影响,洞庭湖三口(松滋口、太平口和藕池口)的水文情势发生了显著变化,给湖区防洪、水资源、水生态、水环境等造成一系列影响.为了深入认识三口洪水发生的复杂变化,本文采用水文变异诊断系统和Zivot-Andrews结构突变单位根检验方法,对三口各水文站点的年最大洪峰流量和年最高洪峰水位序列进行变异诊断,并对其变异特性和变异原因进行分析.结果表明,各站点洪峰流量序列的变异具有较好的一致性;受分流能力变化和上游来流变化的影响,三口各站点的洪峰流量多呈现方向向下的趋势变异;受到洪道冲刷、流量减小、顶托减弱、洪道上下游落差增大的影响,三口各站点洪峰水位在2004年均发生方向向下的跳跃变异.     

Stochastic analysis of flood series

A stochastic model based on the renewal process was developed and used to analyse the characteristics of floods: the volume exceedence, the duration of the flood and the maximum annual flow. The model contains a method for determination of total annual volume exceedence and total annual duration of floods, as well as a method for calculation of maximum annual exceedence, maximum flood duration and maximum flow. The subset of the flood occurrence number in a given time interval is common for all analysed phenomena (volume exceedence, flood duration, maximum flow). The subset of given exceedences is common for total annual volume exceedence, as well as for maximum annual volume exceedence. The same holds for durations of individual floods. The model was then applied to analyse the floods on the Drina River at the Paunci hydrological station and on the Danube River at the Bezdan station.     

Evaluation of dam overtopping probability induced by flood and wind

This study develops a probability-based methodology to evaluate dam overtopping probability that accounts for the uncertainties arising from wind speed and peak flood. A wind speed frequency model and flood frequency analysis, including various distribution types and uncertainties in their parameters, are presented. Furthermore, dam overtopping probabilities based on monthly maximum (MMax) series models are compared with those of the annual maximum (AMax) series models. An efficient sampling scheme, which is a combination of importance sampling (IS) and Latin Hypercube sampling (LHS) methods, is proposed to generate samples of peak flow rate and wind speed especially for rare events. Reservoir routing, which incorporates operation rules, wind setup, and run-up, is used to evaluate dam overtopping probability.     

变化环境下珠江流域洪水频率变化特征、成因及影响(1951-2010年)

气候变化和人类活动导致珠江流域水文变化,变化前后洪水频率分布显著不同.运用滑动秩和(Mann-Whitney U test)结合Brown-Forsythe、滑动T、有序聚类和Mann-Kendall检验法,并用累积距平曲线法获取年最大流量序列详细信息,综合确定样本最佳变化节点,并对水文变化成因做了系统分析.在此基础上,对整体序列、变化前后序列用线性矩法推求广义极值分布参数以及不同重现期设计流量.结果表明:(1)西江大部以及北江流域最佳变化节点在1991年左右;东江流域最佳变化节点与该流域内3大控制性水库建成时间基本吻合;(2)变化后,西江、北江年最大流量持续增加,洪峰强度增大,尤其是西江干流年最大流量显著增加;东江流域年最大流量显著减小,洪峰强度降低;(3)变化后,西江与北江洪水风险增加,尤其是下游珠三角地区本身受人类活动显著影响,加之西江与北江持续增加的洪水强度,珠三角地区发生洪水的强度及频次加剧,而东江洪水风险减小.此研究对于珠江流域在变化环境下的洪水风险评估与防洪抗灾具有重要意义.     

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.     

Trend detection in river flow series: 2. Flood and low-flow index series / Détection de tendance dans des séries de débit fluvial: 2. Séries d'indices de crue et d'étiage

Abstract

Major floods in Europe and North America during the past decade have provoked the question of whether or not they are an effect of a changing climate. This study investigates changes in observational data, using up to 100-year-long daily mean river flow records at 21 stations worldwide. Trends in seven flood and low-flow index series are assessed using Mann-Kendall and linear regression methods. Emphasis was on the comparison of trends in these flow index series, particularly in peak-over-threshold (POT) series as opposed to annual maximum (AM) river flow series. There is a larger number of significant trends in the AM than in the POT flood magnitude series, probably relating to the way the series are constructed. Low flood peaks occurring at the beginning or end of a time series with trend may be too low to be selected for the POT analysis. However, one peak per year will always be selected for the AM series, making the slope steeper and/or the series longer, resulting in a more significant trend. There is no general pattern of increasing or decreasing numbers or magnitudes of floods, but there are significant increases in half of the low-flow series.     

Impulse response of a linear diffusion analogy model as a flood frequency probability density function

Abstract

The impulse response of a linear convective-diffusion analogy (LD) model used for flow routing in open channels is proposed as a probability distribution for flood frequency analysis. The flood frequency model has two parameters, which are derived using the methods of moments and maximum likelihood. Also derived are errors in quantiles for these parameter estimation methods. The distribution shows that the two methods are equivalent in terms of producing mean values—the important property in case of unknown true distribution function. The flood frequency model is tested using annual peak discharges for the gauging sections of 39 Polish rivers where the average value of the ratio of the coefficient of skewness to the coefficient of variation equals about 2.52, a value closer to the ratio of the LD model than to the gamma or the lognormal model. The likelihood ratio indicates the preference of the LD over the lognormal for 27 out of 39 cases. It is found that the proposed flood frequency model represents flood frequency characteristics well (measured by the moment ratio) when the LD flood routing model is likely to be the best of all linear flow routing models.     

Real‐time monitoring and estimation of the discharge of flash floods in a steep mountain catchment

Synchronously and accurately estimating the flood discharges and dynamic changes in the fluid density is essential for hydraulic analysis and forecasting of flash floods, as well as for risk assessment. However, such information is rare for steep mountain catchments, especially in regions that are hotspots for earthquakes. Therefore, six hydrological monitoring sites were established in the main stream and tributaries of the 78.3‐km² Longxi River catchment, an affected region of the Wenchuan earthquake region in China. Direct real‐time monitoring equipment was installed to measure the flow depths, velocities, and fluid total pressures of the flood hydrographs. On the basis of field measurements, real‐time mean cross‐sectional velocities during the flood hydrographs could be derived from easily obtainable parameters: cross‐sectional maximum velocities and the calibrated dimensionless parameter K_h . Real‐time discharges were determined on the basis of a noncontact method to establish the effective rating curves of this mountainous stream, ranging from 1.46 to 386.34 m³/s with the root mean square errors of ≤10.22 m³/s. Compared with the traditional point‐velocity method and empirical Manning's formula, the proposed noncontact method was reliable and safe for monitoring whole flood hydrographs. Additionally, the real‐time fluid density during the flood hydrographs was calculated on the basis of the direct monitoring parameters for fluid total pressures and water depths. During the flood hydrograph, transient flow behaviour with higher fluid density generally occurred downstream during the flood peak periods when the flow was in the supercritical flow regime. The observed behaviour greatly increased the threat of damage to infrastructure and human life near the river. Thus, it is important to accurately estimate flood discharge and identify for fluid densities so that people at risk from an impending flash flood are given reliable, advanced warning.     

A comparison of three methods for regional flood frequency analysis in Saudi Arabia

The objective of the study was to compare the relative accuracy of three methodologies of regional flood frequency analysis in areas of limited flood records. Thirty two drainage basins of different characteristics, located mainly in the southwest region of Saudi Arabia, were selected for the study. In the first methodology, region curves were developed and used together with the mean annual flood, estimated from the characteristics of drainage basin, to estimate flood flows at a location in the basin. The second methodology was to fit probability distribution functions to annual maximum rainfall intensity in a drainage basin. The best fitted probability function was used together with common peak flow models to estimate the annual maximum flood flows in the basin. In the third methodology, duration reduction curves were developed and used together with the average flood flow in a basin to estimate the peak flood flows in the basin. The results obtained from each methodology were compared to the flood records of the selected stations using three statistical measures of goodness-of-fit. The first methodology was found best in a case of having short length of record at a drainage basin. The second methodology produced satisfactory results. Thus, it is recommended in areas where data are not sufficient and/or reliable to utilise the first methodology.     

A theoretically based Wakeby distribution for annual flood series

Abstract

The exact distribution of the ratio of any magnitude to the sum of all magnitudes in an annual flood series satisfying the usual distribution-free assumptions of independence and identical distribution, and the additional parametric assumption of exponential tail behaviour with truncation, is shown to be a beta distribution of the first kind. A two-parameter linear transformation of the beta distribution completes the derivation and yields a Wakeby distribution which has the number of members in a series as a given parameter. The Wakeby distribution is developed to illustrate how, in principle, some perceived deficiencies in current flood frequency analysis may be met: more complex parametric assumptions should lead to distributions of wider application. In particular, the distribution has a secure theoretical basis and is hydrologically more realistic because it bounds the variate and requires the definition of a temporally finite annual series. Analytical expressions are obtained for estimating the two distribution parameters; the quantite standard error and a plotting rule. An example is given of the application of the distribution to the design flood problem and an annual flood series is modelled. It is further suggested that a suitable design value for the largest flood to be withstood by a protection work is a statistic of the largest flood occurring during its lifetime. For the derived Wakeby distribution this criterion specifies risk and probability of non-exceedance of the design flood once a lifetime is selected.     

Seasonal and annual maximum streamflow forecasting using climate information: application to the Three Gorges Dam in the Yangtze River basin,China / Prévision d'écoulements saisonnier et maximum annuel à l'aide d'informations climatiques: application au Barrage des Trois Gorges dans le bassin du Fleuve Yangtze,Chine

Abstract

This paper explores the potential for seasonal prediction of hydrological variables that are potentially useful for reservoir operation of the Three Gorges Dam, China. The seasonal flow of the primary inflow season and the peak annual flow are investigated at Yichang hydrological station, a proxy for inflows to the Three Gorges Dam. Building on literature and diagnostic results, a prediction model is constructed using sea-surface temperatures and upland snow cover available one season ahead of the prediction period. A hierarchical Bayesian approach is used to estimate uncertainty in the parameters of the prediction model and to propagate these uncertainties to the predictand. The results show skill for both the seasonal flow and the peak annual flow. The peak annual flow model is then used to estimate a design flood (50-year flood or 2% exceedence probability) on a year-to-year basis. The results demonstrate the inter-annual variability in flood risk. The predictability of both the seasonal total inflow and the peak annual flow (or a design flood volume) offers potential for adaptive management of the Three Gorges Dam reservoir through modification of the operating policy in accordance with the year-to-year changes in these variables.     

An analysis of change in alpine annual maximum discharges: implications for the selection of design discharges

The paper presents an analysis of 17 long annual maximum series (AMS) of flood flows for Swiss Alpine basins, aimed at checking the presence of changes in the frequency regime of annual maxima. We apply Pettitt's change point test, the nonparametric sign test and Sen's test on trends. We also apply a parametric goodness‐of‐fit test for assessing the suitability of distributions estimated on the basis of annual maxima collected up to a certain year for describing the frequency regime of later observations. For a number of series the tests yield consistent indications for significant changes in the frequency regime of annual maxima and increasing trends in the intensity of annual maximum discharges. In most cases, these changes cannot be explained by anthropogenic causes only (e.g. streamflow regulation, construction of dams). Instead, we observe a statistically significant relationship between the year of change and the elevation of the catchment outlet. This evidence is consistent with the findings of recent studies that explain increasing discharges in alpine catchments with an increase in the temperature controlling the portion of mountain catchments above the freezing point. Finally, we analyse the differences in return periods (RPs) estimated for a given flood flow on the basis of recent and past observations. For a large number of the study AMS, we observe that, on average, the 100‐year flood for past observations corresponds to a RP of approximately 10 to 30 years on the basis of more recent observation. From a complementary perspective, we also notice that estimated RP‐year flood (i.e. flood quantile (FQ) associated with RP) increases on average by approximately 20% for the study area, irrespectively of the RP. Practical implications of the observed changes are illustrated and discussed in the paper. Copyright © 2011 John Wiley & Sons, Ltd.     

Detection of trends in magnitude and frequency of flood peaks across Europe

This study analyses the differences in significant trends in magnitude and frequency of floods detected in annual maximum flood (AMF) and peak over threshold (POT) flood peak series, for the period 1965–2005. Flood peaks are identified from European daily discharge data using a baseflow-based algorithm and significant trends in the AMF series are compared with those in the POT series, derived for six different exceedence thresholds. The results show that more trends in flood magnitude are detected in the AMF than in the POT series and for the POT series more significant trends are detected in flood frequency than in flood magnitude. Spatially coherent patterns of significant trends are detected, which are further investigated by stratifying the results into five regions based on catchment and hydro-climatic characteristics. All data and tools used in this study are open-access and the results are fully reproducible.     

On seasonal and semi-annual approach for flood frequency analysis

As a supplementary method to the conventional flood frequency analysis based on annual peak flows, we propose an approach in this paper to infer the flood frequency distribution on quarterly and semi-annual time scale, which are then converted to annual time scale to obtain the floods corresponding to return periods in unit of year. Two criteria for test of data independence, namely, minimum 7 and 15-day interval between two consecutive peak flows, are tested. The proposed approach was applied to Des Moines River at Fort Dodge, Iowa, USA using its 62 years of observation daily flows. The results show that the estimated floods for given return periods from quarter-annual data series are in general higher than the corresponding estimated floods from semi-annual data series, which is further larger than estimated floods from annual peak flows. The floods estimated from semi-annual data series agree well with the results of previous US Geological Survey study.     

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

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

Hydrological risk analysis of dam overtopping using bivariate statistical approach: a case study from Geheyan Reservoir,China

Hydrological risk analysis is essential and provides useful information for dam safety management and decision-making. This study presents the application of bivariate flood frequency analysis to risk analysis of dam overtopping for Geheyan Reservoir in China. The dependence between the flood peak and volume is modelled with the copula function. A Monte Carlo procedure is conducted to generate 100,000 random flood peak-volume pairs, which are subsequently transformed to corresponding design flood hydrographs (DFHs) by amplifying the selected annual maximum flood hydrographs (AMFHs). These synthetic DFHs are routed through the reservoir to obtain the frequency curve of maximum water level and assess the risk of dam overtopping. Sensitive analysis is performed to investigate the influence of different AMFH shapes and correlation coefficients of flood peak and volume on estimated overtopping risks. The results show that synthetic DFH with AMFH shape characterized by a delayed time to peak results in higher risk, and therefore highlight the importance of including a range of possible AMFH shapes in the dam risk analysis. It is also demonstrated that the overtopping risk is increased as the correlation coefficient of flood peak and volume increases and underestimated in the independence case (i.e. traditional univariate approach), while overestimated in the full dependence case. The bivariate statistical approach based on copulas can effectively capture the actual dependence between flood peak and volume, which should be preferred in the dam risk analysis practice.     

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.     

Quantifying the sources of simulation uncertainty in natural catastrophe models

The risk from natural catastrophes is typically estimated using complex simulation models involving multiple stochastic components in a nested structure. This risk is principally assessed via the mean annual loss, and selected quantiles of the annual loss. Determining an appropriate simulation strategy is important in order to achieve satisfactory convergence of these statistics, without excessive computation time and data storage requirements. This necessitates an understanding of the relative contribution of each of the stochastic components to the total variance of the statistics. A simple framework using random effects models and analysis of variance is used to partition the variance of the annual loss, which permits calculation of the variance of the mean annual loss with varying numbers of samples of each of the components. An extension to quantiles is developed using the empirical distribution function in combination with bootstrapping. The methods are applied to a European flood model, where the primary stochastic component relates to the frequency and severity of flood events, and three secondary components relate to defence levels, exposure locations and building vulnerability. As expected, it is found that the uncertainty due to the secondary components increases as the size of the portfolio of exposures decreases, and is higher for industrial and commercial business, compared with residential for all statistics of interest. In addition, interesting insights are gained as to the impact of flood defences on convergence.