National flood discharge mapping in Austria

This article presents the approach and the results of a study in which 30, 100 and 200 year return period flood discharges were estimated for 26,000 km of Austrian streams. Three guiding principles were adopted: combination of automatic methods and manual assessments by hydrologists to allow speedy processing and account for the local hydrological situation; combination of various sources of information including flood peak samples, rainfall data, runoff coefficients and historical flood data; and involvement of the Hydrographic Services to increase the accuracy and enhance the acceptance of results. The flood discharges for ungauged catchments were estimated by the Top-kriging approach with manual adjustment to the local flood characteristics. The adopted combination approach proved to be very efficient both in terms of the project time required and in terms of the accuracy and acceptability of the estimated flood discharges of given return periods.     

Evaluation of rainfall spatial correlation effect on rainfall-runoff modeling uncertainty,considering 2-copula

Lack of accuracy of rainfall-runoff simulation (RRS) remains critical for some applications. Among various sources of uncertainty, precipitation plays a particular role. Rainfall rates as the main input data of RRS are of the first factors controlling the accuracy. In addition to the depth, spatial and temporal distributions of rainfall impact the flood discharge. Most of the previous studies on RRS uncertainty have ignored rainfall spatial distribution, where in large catchments, it is necessary to be modeled explicitly. Karoon III is one most important basin of the Iran because of the Karoon III dam in the outlet. In the present work, effect of spatial correlation of rainfall on HEC-HMS (SMA) continuous RRS uncertainty is evaluated using 2variate copula (2copula). Monte Carlo simulation (MCS) approach was used to consider the rainfall spatial dependence. To reduce the computational expense, sampling efficiency and convergence for MCS, Latin hypercube sampling (LHS) was used. Copula functions consider wide range of marginal probability distribution functions (PDFs), eliminating limits of regular join PDFs. For this aim, two scenarios were investigated. In the first scenario, sub-basin rainfall was considered independent, and in the second scenario, 2copula was adopted to model spatial correlation of rainfall. Dimensionless rainfall depths were calculated for each sub-basin, and the PDFs were determined. The generated random dimensionless rainfalls were reweighted and multiplied by watershed’s mean rainfall value. Stochastic Climate Library was used to generate continuous daily rainfalls. Sampling from dimensionless rainfalls using LHS algorithm, 100 runs of calibrated model-simulated 100 flows for each day following MCS, and 80 % certainty bound was calculated. Results showed that considering dependence decreased 18 % of the maximum uncertainty bound width, so the methodology could be recommended for decreasing predicted runoff error.     

基于统计理论的产流模型

提出了一个基于统计理论的产流模型,该模型考虑了降雨、土壤下渗能力及土壤蓄水容量的空间变异性。假定每个时段的降雨量在空间上可以用概率密度函数或分布函数描述,根据实测降雨资料通过统计拟合优度途径估计各时段降雨的空间概率分布;采用抛物线型函数分别描述土壤下渗能力和土壤蓄水容量的空间分布。按照超渗产流机制计算地表产流量,通过降雨量和土壤下渗能力的联合分布推导得到地表径流量的统计分布,进而得到平均产流量的解析表达式。下渗水量补充土壤含水量,假定满足田间持水量后形成地下径流,其产流量根据下渗量和土壤蓄水容量的空间分配曲线进行计算。以半湿润的黄河支流伊河东湾流域为例,对模型进行了验证和应用,并与新安江模型的结果进行了对比。结果表明,模型对所研究的半湿润区的洪水模拟预报有较好的模拟效果。     

由实测暴雨推求设计洪水方法的探讨

基于南方湿润地区洪水由暴雨形成的特点，可根据实测暴雨资料分析流域代表性产、汇流参数，将历年最大面暴雨系列转换为对应的洪水系列，并由该洪水系列作为样本来计算洪水统计参数。实例分析结果表明，与实测洪水计算结果相近，该方法不但克服了山区河流雨洪同频率的假定，而且还为研究暴雨形成洪水的产、汇流机制，创造了一定的条件。     

Regional frequency analysis of daily rainfall extremes using L-moments approach

Daily extreme precipitation values are among environmental events with the most disastrous consequences for human society. Information on the magnitudes and frequencies of extreme precipitations is essential for sustainable water resources management, planning for weather-related emergencies, and design of hydraulic structures. In the present study, regional frequency analysis of maximum daily rainfalls was investigated for Golestan province located in the northeastern Iran. This study aimed to find appropriate regional frequency distributions for maximum daily rainfalls and predict the return values of extreme rainfall events (design rainfall depths) for the future. L-moment regionalization procedures coupled with an index rainfall method were applied to maximum rainfall records of 47 stations across the study area. Due to complex geographic and hydro-climatological characteristics of the region, an important research issue focused on breaking down the large area into homogeneous and coherent sub-regions. The study area was divided into five homogeneous regions, based on the cluster analysis of site characteristics and tests for the regional homogeneity. The goodness-of-fit results indicated that the best fitting distribution is different for individual homogeneous regions. The difference may be a result of the distinctive climatic and geographic conditions. The estimated regional quantiles and their accuracy measures produced by Monte Carlo simulations demonstrate that the estimation uncertainty as measured by the RMSE values and 90% error bounds is relatively low when return periods are less than 100 years. But, for higher return periods, rainfall estimates should be treated with caution. More station years, either from longer records or more stations in the regions, would be required for rainfall estimates above T=100 years. It was found from the analyses that, the index rainfall (at-site average maximum rainfall) can be estimated reasonably well as a function of mean annual precipitation in Golestan province. Index rainfalls combined with the regional growth curves, can be used to estimate design rainfalls at ungauged sites. Overall, it was found that cluster analysis together with the L-moments based regional frequency analysis technique could be applied successfully in deriving design rainfall estimates for northeastern Iran. The approach utilized in this study and the findings are of great scientific and practical merit, particularly for the purpose of planning for weather-related emergencies and design of hydraulic engineering structures.     

Monte Carlo analysis of the effect of spatial distribution of storms on prioritization of flood source areas

Implementation of structural and non-structural flood control measures in flood-prone watersheds is on increasing demand. Different watershed areas are not necessarily hydrologically similar and impose variable effects on the outlet flow hydrograph. Meanwhile, prioritization of watershed areas in terms of flood generation is essential for economic flood control planning. Previous works have focused on the definition of a flood index that quantifies the contribution of each subwatershed unit or grid cell to the outlet flood hydrograph through the application of unit flood response (UFR) approach. In the present research, for the first time, the effect of spatial pattern of storm events on the flood index variation was assessed via a Monte Carlo uncertainty analysis. To do so, the UFR approach was carried out for a large number of randomly generated rainfall spatial pattern. The proposed methodology was adopted to the Tangrah watershed in northern Iran. The watershed is frequently hit by floods that have historically caused loss of life and properties. The results indicated that for the more frequent flood events, the flood index is quite sensitive to the spatial distribution of rainfall such that for the highest ranked subwatershed (SW6), the standardized variation of the flood index values (i.e., the uncertainty range) decreases from 1.0 to 0.5 when the rainfall depth increases from 20 to 150 mm, respectively. The results further revealed that increasing the rainfall depth from 20 to 150 mm would cause the effect of rainfall spatial distribution on subwatersheds’ flood indices to diminish. The implications are that if flood control measures are designed for more frequent floods with lower return periods, an uncertainty analysis is required to identify the range of flood index variations.     

Analysis of a Long Record of Annual Maximum Rainfall in Athens,Greece, and Design Rainfall Inferences

An annual series of maximum dailyrainfall extending through 1860–1995, i.e., 136 years,was extracted from the archives of a meteorologicalstation in Athens. This is the longest rainfall recordavailable in Greece and its analysis is required forthe prediction of intense rainfall in Athens, wherecurrently major flood protection works are under way.Moreover, the statistical analysis of this long recordcan be useful for investigating more generalisedissues regarding the adequacy of extreme valuedistributions for extreme rainfall analysis and theeffect of sample size on design rainfall inferences.Statistical exploration and tests based on this longrecord indicate no statistically significant climaticchanges in extreme rainfall during the last 136 years.Furthermore, statistical analysis shows that theconventionally employed Extreme Value Type I (EV1 orGumbel) distribution is inappropriate for the examinedrecord (especially in its upper tail), whereas thisdistribution would seem as an appropriate model iffewer years of measurements were available (i.e., partof this sample were used). On the contrary, theGeneral Extreme Value (GEV) distribution appears to besuitable for the examined series and its predictionsfor large return periods agree with the probablemaximum precipitation estimated by the statistical(Hershfield's) method, when the latter is consideredfrom a probabilistic point of view. Thus, the resultsof the analysis of this record agree with a recently(and internationally) expressed scepticism about theEV1 distribution which tends to underestimate thelargest extreme rainfall amounts. It is demonstratedthat the underestimation is quite substantial (e.g.,1 : 2) for large return periods and this fact must beconsidered as a warning against the widespread use ofthe EV1 distribution for rainfall extremes.     

Variation of uncertainty of drainage density in flood hazard mapping assessment with coupled 1D–2D hydrodynamics model

Coupled 1D–2D hydrodynamic models are widely utilized in flood hazard mapping. Previous studies adopted conceptual hydrological models or 1D hydrodynamic models to evaluate the impact of drainage density on river flow. However, the drainage density affects not only river flow, but also the flooded area and location. Therefore, this work adopts the 1D–2D model SOBEK to investigate the impact of drainage density on river flow. The uncertainty of drainage density in flood hazard mapping is assessed by a designed case and a real case, Yanshuixi Drainage in Tainan, Taiwan. Analytical results indicate that under the same return period rainfall, reduction in tributary drainages in a model (indicating a lower drainage density) results in an underestimate of the flooded area in tributary drainages. This underestimate causes higher peak discharges and total volume of discharges in the drainages, leading to flooding in certain downstream reaches, thereby overestimating the flooded area. The uncertainty of drainage density decreases with increased rainfall. We suggest that modeling flood hazard mapping with low return period rainfalls requires tributary drainages. For extreme rainfall events, a lower drainage density could be selected, but the drainage density of local key areas should be raised.

     

七大流域水文特性分析

我国东西和南北方向地理跨度大,水文特性差异明显。从气候特征、降雨和径流、暴雨和洪水等方面,分析七大流域的水文特性,总结降雨和径流在上下游、干支流等空间上以及年内和年际等时间上的分布规律,并分析暴雨和洪水在成因、发生时间、空间分布、过程等方面的规律,为径流和洪水成果计算提供了基础,也为水文成果的合理性检验提供了参考。     

Hydrologic behavior and flood probability for selected arid basins in Makkah area, western Saudi Arabia

In arid regions, flash floods often occur as a consequence of excessive rainfall. Occasionally causing major loss of property and life, floods are large events of relatively short duration. Makkah area in western Saudi Arabia is characterized by high rainfall intensity that leads to flash floods. This study quantifies the hydrological characteristics and flood probability of some major wadis in western Saudi Arabia, including Na’man, Fatimah, and Usfan. Flood responses in these wadis vary due to the nature and rainfall distribution within these wadis. Rainfall frequency analysis was performed using selected annual maximums of 24-h rainfall from eight stations located in the area. Two of the most applied methods of statistical distribution, Gumbel’s extreme value distribution and log Pearson type III distribution, were applied to maximum daily rainfall data over 26 to 40 years. The Gumbel’s model was found to be the best fitting model for identifying and predicting future rainfall occurrence. Rainfall estimations from different return periods were identified. Probable maximum floods of the major wadis studied were also estimated for different return periods, which were extrapolated from the probable maximum precipitation.     

基于SoilGrids的栅格新安江模型参数空间分布估算

为实现中小流域降雨径流过程精细化模拟, 合理估算水文模型参数的空间分布具有重要意义。基于新版全球数字土壤制图系统(SoilGrids)构建栅格新安江模型(GXM)参数化方案, 对陕西省陈河流域2003—2012年16场洪水进行模拟, 与新安江模型计算结果进行对比, 开展基于洪水过程划分的自由水蓄水容量敏感性及空间分布特征量化分析。结果表明: GXM模拟的峰现时间误差水平降低约0.31 h, 洪峰和洪量模拟精度较高, 模型能够对土壤水饱和度等水文要素的动态空间分布进行较合理的模拟; 自由水蓄水容量参数对洪峰和涨洪过程的确定性系数以及涨洪段的洪量相对误差影响较大, 对退水过程影响小; 自由水蓄水容量在陈河流域河谷和山脊附近较大, 坡段中部较小。     

降雨重现期及其用于滑坡概率分析的探讨

极端气候条件往往会诱发各种地质灾害,而降雨型滑坡的发生则与极值降雨关系密切。为了有效预防和控制汛期滑坡灾害的发生,定量评估滑坡灾害造成的人员与经济风险,讨论了汛期极值降雨条件下滑坡概率的分析方法。利用Gumbel极值分布理论,以三峡库区巴东县1990~2006年日降雨量为基础资料,采用统计分析方法,求取研究区在汛期(6月中旬至9月)最大一次连续降雨量、多日累积最大降雨量的极值及其分布曲线; 在此基础上,以研究区一个滑坡实例为对象,综合采用渗流模拟、稳定性分析和基于蒙特卡罗的滑坡概率分析方法,讨论了降雨极值及其重现期成果分别在降雨新生型滑坡和降雨复活型滑坡概率分析中的应用思路与方法。结果表明,随着重现期的增加,一次降雨过程的降雨量也增加; 随着降雨日数的增加,具有不同重现期降雨事件的累积降雨量均会增加,且重现期越长,累积降雨量值会越大; 降雨极值曲线分别应用于降雨新生型和降雨复活型滑坡概率分析的思路是可行并有效的。     

江苏省无资料山丘区洛阳河流域径流模拟方法探讨

近年来,无资料或资料缺乏地区汇流计算成为水文研究热点。基于地形参数的地貌瞬时单位线(GIUH)方法凭借其对历史降雨径流实测资料的不依赖性,已然成为径流模拟研究中被普遍采用且效果较好的方法之一。为探讨江苏省无资料山丘区的径流模拟方法,选取具有典型地形地貌特征的镇江通胜洛阳河小流域展开研究。根据地貌特征值等数据提取出该流域的地貌单位线,采用地貌单位线计算径流过程,利用2016、2017年汛期实测的4场洪水进行方法验证,将计算径流与实测数据进行对比分析,结果表明:基于地貌特征的地貌单位线法在该流域洪水模拟中具有较高精度,又因该流域于江苏省乃代表性区域,初步判定可将该方法推广至其他无资料山丘区流域。     

地区线性矩法在太湖流域暴雨频率分析中的应用

介绍了基于水文气象途径的地区线性矩法的概念,通过基于次序统计量的线性矩进行参数估计与基于水文气象一致区的地区分析法相结合,以太湖流域1d时段的年极值降雨资料为例,进行暴雨频率分析。应用水文气象一致区的判别准则,将太湖流域划分为8个水文气象一致区;综合考虑三种拟合优度检测方法,选择1~8区的最优分布线型分别为:GEV、GLO、GEV、GEV、GNO、GNO、GEV、GNO;根据地区分析法原理,估算各雨量站的暴雨频率设计值。分析表明:太湖流域各重现期下的年极值降雨空间分布形态基本一致,西南山区是太湖流域的暴雨高风险区,应该在地区防洪规划中引起重视。结果表明:地区线性矩法具有很高的学术和实用价值,建议在全国范围内推广,作为防洪规划的顶层设计和基础工作,以满足工程防洪设计、地区防洪规划、山洪预警和城市防涝防洪规划等方面的需求。     

Combination of GIS and mathematical modeling to predict floods in semiarid areas: case of Rheraya watershed (Western High Atlas,Morocco)

The Tinitine sub-watershed of Rheraya is renowned for its semiarid climate, irregular supply of water flow, and its historical floods; for these reasons, it is classified as vulnerable areas during a rainfall event. We conducted this study to propose flood prediction models adapted to this risky zone. For this, a hydrological method of flood forecasting rainfall-runoff used to determine a model conforms to the semiarid basin. This model is based on the articulation of the series production and transfer function. The goal of the production function is to determine the portion of the rainfall runoff, which is performed by two approaches that differ in their structures: (1) the first approach is based on the mathematical model of Horton; (2) the second approach is based on the determination of the part of the rain contributes to runoff and obtained by a spatial map of runoff coefficient (GIS). The transfer function is based on the two models of Clark and Nash, rainy efficient routes to the catchment outlet from a unit hydrograph. The characteristic parameters of these models are obtained through the application of mathematical optimization methods such as genetic algorithms. Thereby, the coupling function producing Horton (identified by the parameters: initial infiltration?=?15.03, final infiltration?=?0.3, and phi?=?0.45) and Clark transfer function (identified by the parameters C_A?=?0.21 and C_B?=?0.79) was given very satisfying results, mainly the low difference between observed and simulated hydrograph and Nash coefficient which is about 85%. This shows the interest of this coupling model in flood forecasting.     

Rainfall flood hazard at nuclear power plants in India

Flood hazard evaluation is an important input for Nuclear Power Plants external events safety studies. In the present study, flood hazard at various nuclear sites in India due to rainfall has been evaluated. Hazard estimation is a statistical procedure by which rainfall intensity versus occurrence frequency is estimated from historical records of rainfall data and extrapolated with asymptotic extreme value distribution. Rainfall data needed for flood hazard assessment are daily annual maximum rainfall (24?h data). The observed data points have been fitted using Gumbel, power law and exponential distribution, and return period has been estimated. To study the stationarity of rainfall data, a moving window estimate of the parameters has been performed. The rainfall pattern is stationary in both coastal and inland regions over the period of observation. The coastal regions show intense rainfall and higher variability than inland regions. Based on the plant layout, catchment area and drainage capacity, the prototype fast breeder reactor (PFBR) site is unlikely to be flooded.     

20世纪暴雨和洪水极值的年代际变化

利用20世纪100年(特别是后50年)2万多个雨量站和175个大中河流水文站的实测和调查雨量和流量资料,分析了中国大陆(外流区)暴雨和洪水极值的年代际变化,其中包括长期站暴雨的年代均值变化、20世纪50年点雨量极值的前后期平均值对比、最大点雨量发生年代分布和特大洪峰流量的年代分布。分别就5种历时、不同地区的暴雨洪水极值随年代的变化作了分析,揭示了变化的事实,计算了多项年代比率,探讨了影响因子。     

基于非对称Archimedean Copula的三变量洪水风险评估

分析洪峰、洪量和历时三变量联合分布与风险概率及其设计分位数,为水利工程规划设计和风险评估提供参考依据。以珠江流域西江高要站52年洪水数据为例,采用非对称阿基米德M6 Copula函数与Kendall分布函数计算三变量洪水联合分布的“或”重现期、“且”重现期和二次重现期及其最可能的设计分位数。结果表明:“或”重现期的风险率偏高,“且”重现期的风险率偏低,二次重现期更准确地反映了特定设计频率情况下三变量洪水要素遭遇的风险率;按三变量“或”重现期或三变量同频率设计值推算的洪水设计值偏高,以最大可能概率推算的三变量洪水要素的二次重现期设计值可为防洪工程安全与风险管理提供新的选择。     

用相关函数优化法计算离散线性水文模型参数

借助于日降雨和出流资料系列,采用相关函数优化方法确定离散线性水文模型参数是一种新的方法.文中选用4个流域8年日降雨和径流资料系列来检验相关函数优化法(6年作为率定期,2年作为验证期),所得结果与常规最小二乘法计算的结果进行比较表明:新方法计算精度比常规最小二乘法计算的精度有所提高,常规最小二乘法仅是相关函数优化法的一种特例.     

Regional modelling of rainfall-induced runoff using hydrological model by incorporating plant cover effects: case study in Kelantan,Malaysia

The objective of this paper is to develop a spatial temporal runoff modelling of local rainfall patterns effect on the plant cover hilly lands in Kelantan River Basin. Rainfall interception loss based on leaf area index, loss/infiltration on the ground surface, and runoff calculation were considered as the main plant cover effects on the runoff volume. In this regard, a hydrological and geotechnical grid-based regional model (integrated model) was performed using Microsoft Excel^® and GIS framework system for deterministic modelling of rainfall-induced runoff by incorporating plant cover effects. The infiltration process of the current model was integrated with the precipitation distribution method and rainfall interception approach while the runoff analysis of integrated model was employed based on loss/infiltration water on the ground surface with consideration of water interception loss by canopy and the remaining surface water. In the following, the spatial temporal analysis of rainfall-induced runoff was performed using 10 days of hourly rainfall events at the end of December 2014 in Kelantan River Basin. The corresponding changes in pressure head and consequent rate of infiltration were calculated during rainfall events. Subsequently, flood volume is computed using local rainfall patterns, along with water interception loss and the remaining surface water in the study area. The results showed the land cover changes caused significant differences in hydrological response to surface water. The increase in runoff volume of the Kelantan River Basin is as a function of deforestation and urbanization, especially converting the forest area to agricultural land (i.e. rubber and mixed agriculture).